Nucleotide sequences of HS-alpha satellite DNA from kangaroo rat Dipodomys ordii and characterization of similar sequences in other rodents

Cytogenomics of Frieseomelitta varia (Hymenoptera: Apidae) and the Sharing of a Satellite DNA Family in Several Neotropical Meliponini Genera

BACKGROUND/OBJECTIVES

A striking feature of the karyotypes of stingless bees is the large amount of heterochromatin present in most species. Cytogenomic studies performed in some Meliponini species have suggested that evolutionary events related to the diversification and amplification of satellite DNA families in the heterochromatin may reflect the structuring of phylogenetic clades in this tribe. In this study, we performed a genomic analysis in Frieseomelitta varia to characterize different satDNA families in its genome. We also investigated the presence of the most abundant satDNA family of F. varia in its own chromosomes, in two other Frieseomelitta species, and in other Meliponini genera encompassing the three main clades of Neotropical Meliponini, according to the available molecular phylogeny.

METHODS

Genomic analyses were performed using RepeatExplorer2 on the Galaxy platform, and chromosomal investigations were conducted using fluorescent in situ hybridization.

RESULTS

Seven satDNA families were recovered, which together totaled an abundance of 11.223% of the analyzed F. varia genomic fraction. The most abundant satDNA family, FvarSat01-306, predominates in the analyzed repetitive fraction (representing around 89%) and was recently amplified and homogenized in almost all the heterochromatin of F. varia. In addition, the data revealed an unprecedented sharing of this satDNA family in the centromeric/pericentromeric heterochromatin among different Meliponini genera, with independent amplifications and loss of this sequence in some taxa.

CONCLUSIONS

One family of satellite DNA makes up most of the heterochromatin in this species and is shared with other Meliponini.

Satellite DNA shapes dictate pericentromere packaging in female meiosis

The abundance and sequence of satellite DNA at and around centromeres is evolving rapidly despite the highly conserved and essential process through which the centromere directs chromosome inheritance1-3. The impact of such rapid evolution is unclear. Here we find that sequence-dependent DNA shape dictates packaging of pericentromeric satellites in female meiosis through a conserved DNA-shape-recognizing chromatin architectural protein, high mobility group AT-hook 1 (HMGA1)4,5. Pericentromeric heterochromatin in two closely related mouse species, M. musculus and M. spretus, forms on divergent satellites that differ by both density of narrow DNA minor grooves and HMGA1 recruitment. HMGA1 binds preferentially to M. musculus satellites, and depletion in M. musculus oocytes causes massive stretching of pericentromeric satellites, disruption of kinetochore organization and delays in bipolar spindle assembly. In M. musculus × spretus hybrid oocytes, HMGA1 depletion disproportionately impairs M. musculus pericentromeres and microtubule attachment to their kinetochores. Thus, DNA shape affects both pericentromere packaging and the segregation machinery. We propose that rapid evolution of centromere and pericentromere DNA does not disrupt these essential processes when the satellites adopt DNA shapes recognized by conserved architectural proteins (such as HMGA1). By packaging these satellites, architectural proteins become part of the centromeric and pericentromeric chromatin, suggesting an evolutionary strategy that lowers the cost of megabase-scale satellite expansion.

The highly dynamic satellitomes of cultivated wheat species

BACKGROUND AND AIMS

Durum wheat, Triticum turgidum, and bread wheat, Triticum aestivum, are two allopolyploid species of very recent origin that have been subjected to intense selection programmes during the thousands of years they have been cultivated. In this paper, we study the durum wheat satellitome and establish a comparative analysis with the previously published bread wheat satellitome.

METHODS

We revealed the durum wheat satellitome using the satMiner protocol, which is based on consecutive rounds of clustering of Illumina reads by RepeatExplorer2, and we estimated abundance and variation for each identified satellite DNA (satDNA) with RepeatMasker v4.0.5. We also performed a deep satDNA family characterization including chromosomal location by fluorescence in situ hybridization (FISH) in durum wheat and its comparison with FISH patterns in bread wheat. BLAST was used for trailing each satDNA in the assembly of durum wheat genome through NCBI's Genome Data Viewer and the genome assemblies of both species were compared. Sequence divergence and consensus turnover rate between homologous satDNA families of durum and bread wheat were estimated using MEGA11.

KEY RESULTS

This study reveals that in an exceedingly short period, significant qualitative and quantitative changes have occurred in the set of satDNAs of both species, with expansions/contractions of the number of repeats and the loci per satellite, different in each species, and a high rate of sequence change for most of these satellites, in addition to the emergence/loss of satDNAs not shared between the two species analysed. These evolutionary changes in satDNA are common between species but what is truly remarkable and novel is that these processes have taken place in less than the last ~8000 years separating the two species, indicating an accelerated evolution of their satDNAs.

CONCLUSIONS

These results, together with the relationship of many of these satellites with transposable elements and the polymorphisms they generate at the level of centromeres and subtelomeric regions of their chromosomes, are analysed and discussed in the context of the evolutionary origin of these species and the selection pressure exerted by humans throughout the history of their cultivation.

Evidence for the Transcription of a Satellite DNA Widely Found in Frogs

BACKGROUND

The satellite DNA (satDNA) PcP190 has been identified in multiple frog species from seven phylogenetically distant families within Hyloidea, indicating its broad distribution. This satDNA consists of repeats of approximately 190 bp and exhibits a highly conserved region (CR) of 120 bp, which is similar to the transcribed region of 5S ribosomal DNA (rDNA), and a hypervariable region (HR) that varies in size and nucleotide composition among and within species. Here, to improve our understanding of PcP190 satDNA, we searched for evidence of its transcription in the available transcriptomes of Rhinella marina (Bufonidae) and Engystomops pustulosus (Leptodactylidae), two phylogenetically distantly related species.

METHODS

We first characterized the 5S rDNA and PcP190 sequences in these species by searching for them in available genome assemblies. Next, we used the PcP190 (CR and HR) and 5S rDNA sequences of each species as queries to search for these sequences in RNA-seq libraries.

RESULTS

We identified two types of 5S rDNA in each analyzed species, with a new type found in E. pustulosus. Our results also revealed a novel type of PcP190 sequence in R. marina and a new subtype of PcP-1 in E. pustulosus. Transcriptome analyses confirmed the expected transcription of the 5S rRNA gene and showed transcription of both the CR and HR of the PcP190 satDNA in both species and in different tissues.

CONCLUSIONS

As the entire repeat of this satDNA is susceptible to transcription, the high variability observed in the HR cannot be attributed to transcriptional activity confined to the CR.

The burst of satellite DNA in Leptidea wood white butterflies and their putative role in karyotype evolution

Satellite DNAs (satDNAs) are abundant components of eukaryotic genomes, playing pivotal roles in chromosomal organization, genome stability, and evolution. Here, we combined cytogenetic and genomic methods to characterize the satDNAs in the genomes of Leptidea butterflies. Leptidea is characterized by the presence of a high heterochromatin content, large genomes, and extensive chromosomal reshuffling as well as the occurrence of cryptic species. We show that, in contrast to other Lepidoptera, satDNAs constitute a considerable proportion of Leptidea genomes, ranging between 4.11% and 11.05%. This amplification of satDNAs, together with the hyperactivity of transposable elements, contributes to the substantial genome expansion in Leptidea. Using chromosomal mapping, we show that, particularly LepSat01-100 and LepSat03-167 satDNAs, are preferentially localized in heterochromatin exhibiting variable distribution that may have contributed to the highly diverse karyotypes within the genus. The satDNAs also exhibit W-chromosome accumulation, suggesting their involvement in sex chromosome evolution. Our results provide insights into the dynamics of satDNAs in Lepidoptera genomes and highlight their role in genome expansion and chromosomal organization, which could influence the speciation process. The high proportion of repetitive DNAs in the genomes of Leptidea underscores the complex evolutionary dynamics revealing the interplay between repetitive DNAs and genomic architecture in the genus.

Differential amplification and contraction of satellite DNAs in the distinct lineages of the beetle Euchroma gigantea

Satellite DNA (satDNA) consists of tandem repeat sequences that typically evolve rapidly through evolutionary mechanisms, including unequal crossover, transposition events, and others. The evolutionary history of Euchroma gigantea is marked by complex chromosomal evolution between lineages, making this species an interesting model for understanding satDNA evolution at intraspecies level. Therefore, our aim was to comprehend the potential contribution of satDNAs to the greater chromosomal differentiation of evolutionary lineages in E. gigantea by investigating the differential patterns of amplification and contraction of the repeats. To achieve this, we employed de novo identification of satDNA using RepeatExplorer and TAREAN, allowing the satellitome characterization between lineages. A total of 26 satDNA families were identified, ranging from 18 to 1101 nucleotides in length, with most families being shared between individuals/lineages, as predicted by the library hypothesis, except for the satDNA EgiSat21-168 that was absent for Northeast Lineage. The total satDNA content of the individuals was less than 11.2%, and it appeared to increase in two directions following the chromosomal evolution model. Thirteen satDNAs exhibited different patterns of amplification, and nine ones were contracted among individuals. Additionally, most repeats showed a divergence of about 10% for these satDNAs, indicating satellitome differentiation for each lineage/individual. This scenario suggests that the expansion of the satellitome occurred differentially among individuals/lineages of E. gigantea, with the contribution of various DNA turnover mechanisms after geographical isolation, and that they could be involved with karyotype evolution.

Highly divergent satellitomes of two barley species of agronomic importance, Hordeum chilense and H. vulgare

In this paper, we have performed an in-depth study of the complete set of the satellite DNA (satDNA) families (i.e. the satellitomes) in the genome of two barley species of agronomic value in a breeding framework, H. chilense (H1 and H7 accessions) and H. vulgare (H106 accession), which can be useful tools for studying chromosome associations during meiosis. The study has led to the analysis of a total of 18 satDNA families in H. vulgare, 25 satDNA families in H. chilense (accession H1) and 27 satDNA families in H. chilense (accession H7) that constitute 46 different satDNA families forming 36 homology groups. Our study highlights different important contributions of evolutionary and applied interests. Thus, both barley species show very divergent satDNA profiles, which could be partly explained by the differential effects of domestication versus wildlife. Divergence derives from the differential amplification of different common ancestral satellites and the emergence of new satellites in H. chilense, usually from pre-existing ones but also random sequences. There are also differences between the two H. chilense accessions, which support genetically distinct groups. The fluorescence in situ hybridization (FISH) patterns of some satDNAs yield distinctive genetic markers for the identification of specific H. chilense or H. vulgare chromosomes. Some of the satellites have peculiar structures or are related to transposable elements which provide information about their origin and expansion. Among these, we discuss the existence of different (peri)centromeric satellites that supply this region with some plasticity important for centromere evolution. These peri(centromeric) satDNAs and the set of subtelomeric satDNAs (a total of 38 different families) are analyzed in the framework of breeding as the high diversity found in the subtelomeric regions might support their putative implication in chromosome recognition and pairing during meiosis, a key point in the production of addition/substitution lines and hybrids.

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.

In silico Characterization of Satellitomes and Cross-Amplification of Putative satDNAs in Two Species of the Hypostomus ancistroides Complex (Siluriformes, Loricariidae)

INTRODUCTION

The mapping of the satellite DNA on chromosomes is vital to understanding the distribution and evolution of repetitions in the genome since these chromosomal studies have shown the origin, evolutionary mode, and function of repetitive sequences. This study aimed to prospect the satellitome and determine its location in the genome of two cryptic species of Hypostomus, H. aff. ancistroides and H. ancistroides, with and without XX/XY sexual chromosome system.

METHODS

Mitotic chromosomes and DNA extraction were obtained according to protocols. After the whole genome sequencing, the satDNAs were retrieved, amplified, and hybridized in chromosome preparations for male and female individuals.

RESULTS

We found 30 satellite families (47 variants, two superfamilies) in H. ancistroides and 38 satellite families (45 variants, four superfamilies) in H. aff. ancistroides. The sequences varied from 14 bp to 2,662 bp in H. ancistroides and from 14 bp to 2,918 bp in H. aff. ancistroides. We did not observe any tandem repeats that were exclusive to each of the libraries; however, many sequences showed very different abundances and copy numbers between the libraries. Four satDNAs did not hybridize on the chromosomes of either species. Conversely, one satDNA hybridized in both species, HxySat1-80. However, the phenotypes found varied among species, populations, and in the same individual. There was no sign of HanSat3-464 and HanSat11-335 in any individuals of H. aff. ancistroides, but markings were in the chromosomes of H. ancistroides. HxySat12-1127 and HxySat8-52, on the other hand, were only hybridized in H. aff. ancistroides, while H. ancistroides had a negative sign. No hybridization of satDNAs was found in the X and Y sex chromosomes as they were mostly composed of euchromatin.

CONCLUSION

We distinguish H. aff. ancistroides as genetically different from H. ancistroides, recognizing that such characteristics go far beyond morphological, karyotypic, and molecular data. Our data support the differential abundance and location of satellite DNAs and confirm that many organisms, including fish, have repetitive sequences that validate the library hypothesis. All found and validated satDNAs and the characterization of the satellitomes of the two species represent important contributions to cytogenomic studies of the genus Hypostomus.

Interspecific cytogenomic comparison reveals a potential chromosomal centromeric marker in Proceratophrys frog species

Among the repetitive elements, satellite DNA (SatDNA) emerges as extensive arrays of highly similar tandemly repeated units, spanning megabases in length. Given that the satDNA PboSat01-176, previously characterized in P. boiei, prompted our interest for having a high abundance in P. boiei and potential for centromeric satellite, here, we employed various approaches, including low coverage genome sequencing, followed by computational analysis and chromosomal localization techniques in four Proceratophrys species and, investigating the genomic presence and sharing, as well as its potential for chromosomal centromere marker in Proceratophrys frog species. Our findings demonstrate that PboSat01-176 exhibits high abundance across all four Proceratophrys species, displaying distinct characteristics that establish it as the predominant repetitive DNA element in these species. The satellite DNA is prominently clustered in the peri/centromeric region of the chromosomes, particularly in the heterochromatic regions. The widespread presence of PboSat01-176 in closely related Proceratophrys species reinforces the validity of the library hypothesis for repetitive sequences. Thus, this study highlighted the utility of the satDNA family PboSat01-176 as a reliable centromeric marker in Proceratophrys species, with potential to be applied in other species of anuran amphibians. The observed sharing and maintenance of this sequence within the genus suggest possibilities for future research, particularly through expanded sampling to elucidate parameters that underlie the library hypothesis and the evolutionary dynamics of satDNA 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.

Evolution of bird sex chromosomes: a cytogenomic approach in Palaeognathae species

BACKGROUND

Different patterns of sex chromosome differentiation are seen in Palaeognathae birds, a lineage that includes the ratites (Struthioniformes, Rheiformes, Apterygiformes, Casuariiformes, and the sister group Tinamiformes). While some Tinamiform species have well-differentiated W chromosomes, both Z and W of all the flightless ratites are still morphologically undifferentiated. Here, we conducted a comprehensive analysis of the ZW differentiation in birds using a combination of cytogenetic, genomic, and bioinformatic approaches. The whole set of satDNAs from the emu (Dromaius novaehollandiae) was described and characterized. Furthermore, we examined the in situ locations of these satDNAs alongside several microsatellite repeats and carried out Comparative Genomic Hybridizations in two related species: the greater rhea (Rhea americana) and the tataupa tinamou (Crypturellus tataupa).

RESULTS

From the 24 satDNA families identified (which represent the greatest diversity of satDNAs ever uncovered in any bird species), only three of them were found to accumulate on the emu's sex chromosomes, with no discernible accumulation observed on the W chromosome. The W chromosomes of both the greater rhea and the emu did not exhibit a significant buildup of either C-positive heterochromatin or repetitive DNAs, indicating their large undifferentiation both at morphological and molecular levels. In contrast, the tataupa tinamou has a highly differentiated W chromosome that accumulates several DNA repeats.

CONCLUSION

The findings provide new information on the architecture of the avian genome and an inside look at the starting points of sex chromosome differentiation in birds.

Evolution of ancient satellite DNAs in extant alligators and caimans (Crocodylia, Reptilia)

BACKGROUND

Crocodilians are one of the oldest extant vertebrate lineages, exhibiting a combination of evolutionary success and morphological resilience that has persisted throughout the history of life on Earth. This ability to endure over such a long geological time span is of great evolutionary importance. Here, we have utilized the combination of genomic and chromosomal data to identify and compare the full catalogs of satellite DNA families (satDNAs, i.e., the satellitomes) of 5 out of the 8 extant Alligatoridae species. As crocodilian genomes reveal ancestral patterns of evolution, by employing this multispecies data collection, we can investigate and assess how satDNA families evolve over time.

RESULTS

Alligators and caimans displayed a small number of satDNA families, ranging from 3 to 13 satDNAs in A. sinensis and C. latirostris, respectively. Together with little variation both within and between species it highlighted long-term conservation of satDNA elements throughout evolution. Furthermore, we traced the origin of the ancestral forms of all satDNAs belonging to the common ancestor of Caimaninae and Alligatorinae. Fluorescence in situ experiments showed distinct hybridization patterns for identical orthologous satDNAs, indicating their dynamic genomic placement.

CONCLUSIONS

Alligators and caimans possess one of the smallest satDNA libraries ever reported, comprising only four sets of satDNAs that are shared by all species. Besides, our findings indicated limited intraspecific variation in satellite DNA, suggesting that the majority of new satellite sequences likely evolved from pre-existing ones.

Bread wheat satellitome: a complex scenario in a huge genome

In bread wheat (Triticum aestivum L.), chromosome associations during meiosis are extremely regulated and initiate at the telomeres and subtelomeres, which are enriched in satellite DNA (satDNA). We present the study and characterization of the bread wheat satellitome to shed light on the molecular organization of wheat subtelomeres. Our results revealed that the 2.53% of bread wheat genome is composed by satDNA and subtelomeres are particularly enriched in such DNA sequences. Thirty-four satellite DNA (21 for the first time in this work) have been identified, analyzed and cytogenetically validated. Many of the satDNAs were specifically found at particular subtelomeric chromosome regions revealing the asymmetry in subtelomere organisation among the wheat subgenomes, which might play a role in proper homologous recognition and pairing during meiosis. An integrated physical map of the wheat satellitome was also constructed. To the best of our knowledge, our results show that the combination of both cytogenetics and genome research allowed the first comprehensive analysis of the wheat satellitome, shedding light on the complex wheat genome organization, especially on the polymorphic nature of subtelomeres and their putative implication in chromosome recognition and pairing during meiosis.

Conserved satellite DNA motif and lack of interstitial telomeric sites in highly rearranged African Nothobranchius killifish karyotypes

Using African annual killifishes of the genus Nothobranchius from temporary savannah pools with rapid karyotype and sex chromosome evolution, we analysed the chromosomal distribution of telomeric (TTAGGG)n repeat and Nfu-SatC satellite DNA (satDNA; isolated from Nothobranchius furzeri) in 15 species across the Nothobranchius killifish phylogeny, and with Fundulosoma thierryi as an out-group. Our fluorescence in situ hybridization experiments revealed that all analysed taxa share the presence of Nfu-SatC repeat but with diverse organization and distribution on chromosomes. Nfu-SatC landscape was similar in conspecific populations of Nothobranchius guentheri and Nothobranchius melanospilus but slightly-to-moderately differed between populations of Nothobranchius pienaari, and between closely related Nothobranchius kuhntae and Nothobranchius orthonotus. Inter-individual variability in Nfu-SatC patterns was found in N. orthonotus and Nothobranchius krysanovi. We revealed mostly no sex-linked patterns of studied repetitive DNA distribution. Only in Nothobranchius brieni, possessing multiple sex chromosomes, Nfu-SatC repeat occupied a substantial portion of the neo-Y chromosome, similarly as formerly found in the XY sex chromosome system of turquoise killifish N. furzeri and its sister species Nothobranchius kadleci-representatives not closely related to N. brieni. All studied species further shared patterns of expected telomeric repeats at the ends of all chromosomes and no additional interstitial telomeric sites. In summary, we revealed (i) the presence of conserved satDNA class in Nothobranchius clades (a rare pattern among ray-finned fishes); (ii) independent trajectories of Nothobranchius sex chromosome differentiation, with recurrent and convergent accumulation of Nfu-SatC on the Y chromosome in some species; and (iii) genus-wide shared tendency to loss of telomeric repeats during interchromosomal rearrangements. Collectively, our findings advance our understanding of genome structure, mechanisms of karyotype reshuffling, and sex chromosome differentiation in Nothobranchius killifishes from the genus-wide perspective.

Karyotypes of water frogs from the Pelophylax esculentus complex: results of cross-species chromosomal painting

Amphibian species have the largest genome size enriched with repetitive sequences and relatively similar karyotypes. Moreover, many amphibian species frequently hybridize causing nuclear and mitochondrial genome introgressions. In addition, hybridization in some amphibian species may lead to clonality and polyploidization. All such events were found in water frogs from the genus Pelophylax. Among the species within the genus Pelophylax, P. esculentus complex is the most widely distributed and well-studied. This complex includes two parental species, P. ridibundus and P. lessonae, and their hybrids, P. esculentus, reproducing hemiclonally. Parental species and their hybrids have similar but slightly polymorphic karyotypes, so their precise identification is still required. Here, we have developed a complete set of 13 chromosome painting probes for two parental species allowing the precise identification of all chromosomes. Applying chromosomal painting, we identified homologous chromosomes in both parental species and orthologous chromosomes in their diploid hemiclonal hybrids. Comparative painting did not reveal interchromosomal exchanges between the studied water frog species and their hybrids. Using cross-specific chromosome painting, we detected unequal distribution of the signals along chromosomes suggesting the presence of species-specific tandem repeats. Application of chromosomal paints to the karyotypes of hybrids revealed differences in the intensity of staining for P. ridibundus and P. lessonae chromosomes. Thus, both parental genomes have a divergence in unique sequences. Obtained chromosome probes may serve as a powerful tool to unravel chromosomal evolution in phylogenetically related species, identify individual chromosomes in different cell types, and investigate the elimination of chromosomes in hybrid water frogs.

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.

Investigating the diversification of holocentromeric satellite DNA Tyba in Rhynchospora (Cyperaceae)

BACKGROUND AND AIMS

Satellite DNAs (satDNAs) are repetitive sequences composed by tandemly arranged, often highly homogenized units called monomers. Although satDNAs are usually fast evolving, some satDNA families can be conserved across species separated by several millions of years, probably because of their functional roles in the genomes. Tyba was the first centromere-specific satDNA described for a holocentric organism, until now being characterized for only eight species of the genus Rhynchospora Vahl. (Cyperaceae). Here, we characterized Tyba across a broad sampling of the genus, analysing and comparing its evolutionary patterns with other satDNAs.

METHODS

We characterized the structure and sequence evolution of satDNAs across a robust dadated phylogeny based on Hybrid Target-Capture Sequencing (hyb-seq) of 70 species. We mined the repetitive fraction for Tyba-like satellites to compare its features with other satDNAs and to construct a Tyba-based phylogeny for the genus.

KEY RESULTS

Our results show that Tyba is present in the majority of examined species of the genus, spanning four of the five major clades and maintaining intrafamily pairwise identity of 70.9% over 31 Myr. In comparison, other satellite families presented higher intrafamily pairwise identity but are phylogenetically restricted. Furthermore, Tyba sequences could be divided into 12 variants grouped into three different clade-specific subfamilies, showing evidence of traditional models of satDNA evolution, such as the concerted evolution and library models. Besides, a Tyba-based phylogeny showed high congruence with the hyb-seq topology. Our results show structural indications of a possible relationship of Tyba with nucleosomes, given its high curvature peaks over conserved regions and overall high bendability values compared with other non-centromeric satellites.

CONCLUSIONS

Overall, Tyba shows a remarkable sequence conservation and phylogenetic significance across the genus Rhynchospora, which suggests that functional roles might lead to long-term stability and conservation for satDNAs in the genome.

The Low-Copy-Number Satellite DNAs of the Model Beetle Tribolium castaneum

The red flour beetle Tribolium castaneum is an important pest of stored agricultural products and the first beetle whose genome was sequenced. So far, one high-copy-number and ten moderate-copy-number satellite DNAs (satDNAs) have been described in the assembled part of its genome. In this work, we aimed to catalog the entire collection of T. castaneum satDNAs. We resequenced the genome using Illumina technology and predicted potential satDNAs via graph-based sequence clustering. In this way, we discovered 46 novel satDNAs that occupied a total of 2.1% of the genome and were, therefore, considered low-copy-number satellites. Their repeat units, preferentially 140-180 bp and 300-340 bp long, showed a high A + T composition ranging from 59.2 to 80.1%. In the current assembly, we annotated the majority of the low-copy-number satDNAs on one or a few chromosomes, discovering mainly transposable elements in their vicinity. The current assembly also revealed that many of the in silico predicted satDNAs were organized into short arrays not much longer than five consecutive repeats, and some of them also had numerous repeat units scattered throughout the genome. Although 20% of the unassembled genome sequence masked the genuine state, the predominance of scattered repeats for some low-copy satDNAs raises the question of whether these are essentially interspersed repeats that occur in tandem only sporadically, with the potential to be satDNA "seeds".

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.

Satellite DNA evolution in Corvoidea inferred from short and long reads

Satellite DNA (satDNA) is a fast-evolving portion of eukaryotic genomes. The homogeneous and repetitive nature of such satDNA causes problems during the assembly of genomes, and therefore it is still difficult to study it in detail in nonmodel organisms as well as across broad evolutionary timescales. Here, we combined the use of short- and long-read data to explore the diversity and evolution of satDNA between individuals of the same species and between genera of birds spanning ~40 millions of years of bird evolution using birds-of-paradise (Paradisaeidae) and crow (Corvus) species. These avian species highlighted the presence of a GC-rich Corvoidea satellitome composed of 61 satellite families and provided a set of candidate satDNA monomers for being centromeric on the basis of length, abundance, homogeneity and transcription. Surprisingly, we found that the satDNA of crow species rapidly diverged between closely related species while the satDNA appeared more similar between birds-of-paradise species belonging to different genera.

A Satellite-Free Centromere in Equus przewalskii Chromosome 10

In mammals, centromeres are epigenetically specified by the histone H3 variant CENP-A and are typically associated with satellite DNA. We previously described the first example of a natural satellite-free centromere on Equus caballus chromosome 11 (ECA11) and, subsequently, on several chromosomes in other species of the genus Equus. We discovered that these satellite-free neocentromeres arose recently during evolution through centromere repositioning and/or chromosomal fusion, after inactivation of the ancestral centromere, where, in many cases, blocks of satellite sequences were maintained. Here, we investigated by FISH the chromosomal distribution of satellite DNA families in Equus przewalskii (EPR), demonstrating a good degree of conservation of the localization of the major horse satellite families 37cen and 2PI with the domestic horse. Moreover, we demonstrated, by ChIP-seq, that 37cen is the satellite bound by CENP-A and that the centromere of EPR10, the ortholog of ECA11, is devoid of satellite sequences. Our results confirm that these two species are closely related and that the event of centromere repositioning which gave rise to EPR10/ECA11 centromeres occurred in the common ancestor, before the separation of the two horse lineages.

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.

Tandem Repeat DNA Provides Many Cytological Markers for Hybrid Zone Analysis in Two Subspecies of the Grasshopper Chorthippus parallelus

Recent advances in next generation sequencing (NGS) have greatly increased our understanding of non-coding tandem repeat (TR) DNA. Here we show how TR DNA can be useful for the study of hybrid zones (HZ), as it serves as a marker to identify introgression in areas where two biological entities come in contact. We used Illumina libraries to analyse two subspecies of the grasshopper Chorthippus parallelus, which currently form a HZ in the Pyrenees. We retrieved a total of 152 TR sequences, and used fluorescent in situ hybridization (FISH) to map 77 families in purebred individuals from both subspecies. Our analysis revealed 50 TR families that could serve as markers for analysis of this HZ, using FISH. Differential TR bands were unevenly distributed between chromosomes and subspecies. Some of these TR families yielded FISH bands in only one of the subspecies, suggesting the amplification of these TR families after the geographic separation of the subspecies in the Pleistocene. Our cytological analysis of two TR markers along a transect of the Pyrenean hybrid zone showed asymmetrical introgression of one subspecies into the other, consistent with previous findings using other markers. These results demonstrate the reliability of TR-band markers for hybrid zone studies.

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.

Satellitome Analysis on Talpa aquitania Genome and Inferences about the satDNAs Evolution on Some Talpidae

In the genus Talpa a new species, named Talpa aquitania, has been recently described. Only cytogenetic data are available for the nuclear genome of this species. In this work, we characterize the satellitome of the T. aquitania genome that presents 16 different families, including telomeric sequences, and they represent 1.24% of the genome. The first satellite DNA family (TaquSat1-183) represents 0.558%, and six more abundant families, including TaquSat1-183, comprise 1.13%, while the remaining 11 sat-DNAs represent only 0.11%. The average A + T content of the SatDNA families was 50.43% and the median monomer length was 289.24 bp. The analysis of these SatDNAs indicated that they have different grades of clusterization, homogenization, and degeneration. Most of the satDNA families are present in the genomes of the other Talpa species analyzed, while in the genomes of other more distant species of Talpidae, only some of them are present, in accordance with the library hypothesis. Moreover, chromosomal localization by FISH revealed that some satDNAs are localized preferentially on centromeric and non-centromeric heterochromatin in T. aquitania and also in the sister species T. occidentalis karyotype. The differences observed between T. aquitania and the close relative T. occidentalis and T. europaea suggested that the satellitome is a very dynamic component of the genomes and that the satDNAs could be responsible for chromosomal differences between the species. Finally, in a broad context, these data contribute to the understanding of the evolution of satellitomes on mammals.

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

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

Satellitome analyses in nematodes illuminate complex species history and show conserved features in satellite DNAs

Background

Satellite DNAs (satDNAs) are tandemly repeated non-coding DNA sequences that belong to the most abundant and the fastest evolving parts of the eukaryotic genome. A satellitome represents the collection of different satDNAs in a genome. Due to extreme diversity and methodological difficulties to characterize and compare satDNA collection in complex genomes, knowledge on their putative functional constraints and capacity to participate in genome evolution remains rather elusive. SatDNA transcripts have been detected in many species, however comparative studies of satDNA transcriptome between species are extremely rare.

Results

We conducted a genome-wide survey and comparative analyses of satellitomes among different closely related Meloidogyne spp. nematodes. The evolutionary trends of satDNAs suggest that each round of proposed polyploidization in the evolutionary history is concomitant with the addition of a new set of satDNAs in the satellitome of any particular Meloidogyne species. Successive incorporation of new sets of satDNAs in the genome along the process of polyploidization supports multiple hybridization events as the main factor responsible for the formation of these species. Through comparative analyses of 83 distinct satDNAs, we found a CENP-B box-like sequence motif conserved among 11 divergent satDNAs (similarity ranges from 36 to 74%). We also found satDNAs that harbor a splice leader (SL) sequence which, in spite of overall divergence, shows conservation across species in two putative functional regions, the 25-nt SL exon and the Sm binding site. Intra- and interspecific comparative expression analyses of the complete satDNA set in the analyzed Meloidogyne species revealed transcription profiles including a subset of 14 actively transcribed satDNAs. Among those, 9 show active transcription in every species where they are found in the genome and throughout developmental stages.

Conclusions

Our results demonstrate the feasibility and power of comparative analysis of the non-coding repetitive genome for elucidation of the origin of species with a complex history. Although satDNAs generally evolve extremely quickly, the comparative analyses of 83 satDNAs detected in the analyzed Meloidogyne species revealed conserved sequence features in some satDNAs suggesting sequence evolution under selective pressure. SatDNAs that are actively transcribed in related genomes and throughout nematode development support the view that their expression is not stochastic.

Identification and characterization of a new family of long satellite DNA, specific of true toads (Anura, Amphibia, Bufonidae)

Amphibians have some of the most variable genome sizes among vertebrates. Genome size variation has been attributed to repetitive and noncoding DNA, including satellite repeats, transposable elements, introns, and nuclear insertions of viral and organelle DNA. In vertebrates, satellite DNAs have been widely described in mammals, but few molecular studies have been carried out in amphibians. Here, we provide a detailed characterization of a new family of satellite DNA, present in all 15 examined species of the family Bufonidae. Southern-blot analysis and PCR reveal that this satellite is formed by monomers of 807 bp, is organized in tandem arrays, and has an AT-content of 57.4%. Phylogenetic analyses show that most clades exhibit species-specific variances, indicating that this satellite DNA has evolved by concerted evolution. The homogenization/fixation process is heterogeneous in Bufonidae, where the genera Bufo and Bufotes do not show species-specific differences, while populations from Rhinella marina exhibit population-specific changes. Additionally, variants of this satellite DNA have been identified in Duttaphrynus melanostictus and R. marina, supporting the 'library hypothesis' (a set, 'library', of satellite DNAs is shared by a species group). Physical mapping in Bufo bufo, Bufo spinosus, Epidalea calamita and Bufotes viridis provides evidence that this repetitive DNA is not dispersed in the karyotype, but accumulated in pericentromeric regions of some chromosomal pairs. This location, together with its presence in the transcriptomes of bufonids, could indicate a role in centromere function or heterochromatin formation and maintenance.

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.

Study and Physical Mapping of the Species-Specific Tandem Repeat CS-237 Linked with 45S Ribosomal DNA Intergenic Spacer in Cannabis sativa L

Hemp (Cannabis sativa L.) is a valuable crop and model plant for studying sex chromosomes. The scientific interest in the plant has led to its whole genome sequencing and the determination of its cytogenetic characteristics. A range of cytogenetic markers (subtelomeric repeat CS-1, 5S rDNA, and 45S rDNA) has been mapped onto hemp's chromosomes by fluorescent in situ hybridization (FISH). In this study, another cytogenetic marker (the tandem repeat CS-237, with a 237 bp monomer) was found, studied, and localized on chromosomes by FISH. The signal distribution and karyotyping revealed that the CS-237 probe was localized in chromosome 6 with one hybridization site and in chromosome 8 with two hybridization sites, one of which colocalizes with the 45S rDNA probe (with which a nucleolus organizer region, NOR, was detected). A BLAST analysis of the genomic data and PCR experiments showed that the modified CS-237 monomers (delCS-237, 208 bp in size) were present in the intergenic spacers (IGSs) of hemp 45S rDNA monomers. Such a feature was firstly observed in Cannabaceae species. However, IGS-linked DNA repeats were found in several plant species of other families (Fabaceae, Solanaceae, and Asteraceae). This phenomenon is discussed in this article. The example of CS-237 may be useful for further studying the phenomenon as well as for the physical mapping of hemp chromosomes.

Satellitome of the Red Palm Weevil, Rhynchophorus ferrugineus (Coleoptera: Curculionidae), the Most Diverse Among Insects

The red palm weevil, Rhynchophorus ferrugineus, is the most harmful species among those pests affecting palm trees. Its impact causes important economic losses around the World. Nevertheless, the genetic information of Rh. ferrugineus is very scarce. Last year, the first genome assembly was published including a rough description of its repeatome. However, no information has been added about one of the main components of repeated DNA, the satellite DNA. Herein, we presented the characterization of the satellitome of this important species that includes 112 satellite DNA families, the largest number in an insect genome. These satellite DNA families made up around 25% of the genome while the most abundant family, RferSat01-169, alone represented 20.4%. Chromosomal location of most abundant satellite DNA families performed by fluorescence in situ hybridization showed that all of them are dispersed in the euchromatin on all chromosomes but some of them are also specifically accumulated either on the pericentromeric heterochromatic regions of all chromosomes or on specific chromosomes. Finally, the transcription of satellitome families was analyzed through Rh. ferrugineus development. It was found that 55 out of 112 satellite DNA families showed transcription, some families seemed to be transcribed across all stages while a few appeared to be stage-specific, indicating a possible role of those satellite DNA sequences in the development of this species.

The Singular Evolution of Olea Genome Structure

The current view of plant genome evolution proposes that genome size has mainly been determined by polyploidisation and amplification/loss of transposons, with a minor role played by other repeated sequences, such as tandem repeats. In cultivated olive (Olea europaea subsp. europaea var. europaea), available data suggest a singular model of genome evolution, in which a massive expansion of tandem-repeated sequences accompanied changes in nuclear architecture. This peculiar scenario highlights the importance of focusing on Olea genus evolution, to shed light on mechanisms that led to its present genomic structure. Next-generation sequencing technologies, bioinformatics and in situ hybridisation were applied to study the genomic structure of five related Olea taxa, which originated at different times from their last common ancestor. On average, repetitive DNA in the Olea taxa ranged from ~59% to ~73% of the total genome, showing remarkable differences in terms of composition. Among repeats, we identified 11 major families of tandem repeats, with different abundances in the analysed taxa, five of which were novel discoveries. Interestingly, overall tandem repeat abundance was inversely correlated to that of retrotransposons. This trend might imply a competition in the proliferation of these repeat classes. Indeed, O. paniculata, the species closest to the Olea common ancestor, showed very few tandem-repeated sequences, while it was rich in long terminal repeat retrotransposons, suggesting that the amplification of tandem repeats occurred after its divergence from the Olea ancestor. Furthermore, some tandem repeats were physically localised in closely related O. europaea subspecies (i.e., cultivated olive and O. europaea subsp. cuspidata), which showed a significant difference in tandem repeats abundance. For 4 tandem repeats families, a similar number of hybridisation signals were observed in both subspecies, apparently indicating that, after their dissemination throughout the olive genome, these tandem repeats families differentially amplified maintaining the same positions in each genome. Overall, our research identified the temporal dynamics shaping genome structure during Olea speciation, which represented a singular model of genome evolution in higher plants.

Satellitome comparison of two oedipodine grasshoppers highlights the contingent nature of satellite DNA evolution

BACKGROUND

The full catalog of satellite DNA (satDNA) within a same genome constitutes the satellitome. The Library Hypothesis predicts that satDNA in relative species reflects that in their common ancestor, but the evolutionary mechanisms and pathways of satDNA evolution have never been analyzed for full satellitomes. We compare here the satellitomes of two Oedipodine grasshoppers (Locusta migratoria and Oedaleus decorus) which shared their most recent common ancestor about 22.8 Ma ago.

RESULTS

We found that about one third of their satDNA families (near 60 in every species) showed sequence homology and were grouped into 12 orthologous superfamilies. The turnover rate of consensus sequences was extremely variable among the 20 orthologous family pairs analyzed in both species. The satDNAs shared by both species showed poor association with sequence signatures and motives frequently argued as functional, except for short inverted repeats allowing short dyad symmetries and non-B DNA conformations. Orthologous satDNAs frequently showed different FISH patterns at both intra- and interspecific levels. We defined indices of homogenization and degeneration and quantified the level of incomplete library sorting between species.

CONCLUSIONS

Our analyses revealed that satDNA degenerates through point mutation and homogenizes through partial turnovers caused by massive tandem duplications (the so-called satDNA amplification). Remarkably, satDNA amplification increases homogenization, at intragenomic level, and diversification between species, thus constituting the basis for concerted evolution. We suggest a model of satDNA evolution by means of recursive cycles of amplification and degeneration, leading to mostly contingent evolutionary pathways where concerted evolution emerges promptly after lineages split.

The extensive amplification of heterochromatin in Melipona bees revealed by high throughput genomic and chromosomal analysis

Satellite DNAs (satDNAs) and transposable elements (TEs) are among the main components of constitutive heterochromatin (c-heterochromatin) and are related to their functionality, dynamics, and evolution. A peculiar case regarding the quantity and distribution of c-heterochromatin is observed in the genus of bees, Melipona, with species having a low amount of heterochromatin and species with high amount occupying almost all chromosomes. By combining low-pass genome sequencing and chromosomal analysis, we characterized the satDNAs and TEs of Melipona quadrifasciata (low c-heterochromatin) and Melipona scutellaris (high low c-heterochromatin) to understand c-heterochromatin composition and evolution. We identified 15 satDNA families and 20 TEs for both species. Significant variations in the repeat landscapes were observed between the species. In M. quadrifasciata, the repetitive fraction corresponded to only 3.78% of the genome library studied, whereas in M. scutellaris, it represented 54.95%. Massive quantitative and qualitative changes contributed to the differential amplification of c-heterochromatin, mainly due to the amplification of exclusive repetitions in M. scutellaris, as the satDNA MscuSat01-195 and the TE LTR/Gypsy_1 that represent 38.20 and 14.4% of its genome, respectively. The amplification of these two repeats is evident at the chromosomal level, with observation of their occurrence on most c-heterochromatin. Moreover, we detected repeats shared between species, revealing that they experienced mainly quantitative variations and varied in the organization on chromosomes and evolutionary patterns. Together, our data allow the discussion of patterns of evolution of repetitive DNAs and c-heterochromatin that occurred in a short period of time, after separation of the Michmelia and Melipona subgenera.

Chromosome and Genome Diversity in the Genus Trifolium (Fabaceae)

Trifolium L. is an economically important genus that is characterized by variable karyotypes relating to its ploidy level and basic chromosome numbers. The advent of genomic resources combined with molecular cytogenetics provides an opportunity to develop our understanding of plant genomes in general. Here, we summarize the current state of knowledge on Trifolium genomes and chromosomes and review methodologies using molecular markers that have contributed to Trifolium research. We discuss possible future applications of cytogenetic methods in research on the Trifolium genome and chromosomes.

High-Throughput Genomic Data Reveal Complex Phylogenetic Relationships in Stylosanthes Sw (Leguminosae)

Allopolyploidy is widely present across plant lineages. Though estimating the correct phylogenetic relationships and origin of allopolyploids may sometimes become a hard task. In the genus Stylosanthes Sw. (Leguminosae), an important legume crop, allopolyploidy is a key speciation force. This makes difficult adequate species recognition and breeding efforts on the genus. Based on comparative analysis of nine high-throughput sequencing (HTS) samples, including three allopolyploids (S. capitata Vogel cv. “Campo Grande,” S. capitata “RS024” and S. scabra Vogel) and six diploids (S. hamata Taub, S. viscosa (L.) Sw., S. macrocephala M. B. Ferreira and Sousa Costa, S. guianensis (Aubl.) Sw., S. pilosa M. B. Ferreira and Sousa Costa and S. seabrana B. L. Maass & 't Mannetje) we provide a working pipeline to identify organelle and nuclear genome signatures that allowed us to trace the origin and parental genome recognition of allopolyploids. First, organelle genomes were de novo assembled and used to identify maternal genome donors by alignment-based phylogenies and synteny analysis. Second, nuclear-derived reads were subjected to repetitive DNA identification with RepeatExplorer2. Identified repeats were compared based on abundance and presence on diploids in relation to allopolyploids by comparative repeat analysis. Third, reads were extracted and grouped based on the following groups: chloroplast, mitochondrial, satellite DNA, ribosomal DNA, repeat clustered- and total genomic reads. These sets of reads were then subjected to alignment and assembly free phylogenetic analyses and were compared to classical alignment-based phylogenetic methods. Comparative analysis of shared and unique satellite repeats also allowed the tracing of allopolyploid origin in Stylosanthes, especially those with high abundance such as the StyloSat1 in the Scabra complex. This satellite was in situ mapped in the proximal region of the chromosomes and made it possible to identify its previously proposed parents. Hence, with simple genome skimming data we were able to provide evidence for the recognition of parental genomes and understand genome evolution of two Stylosanthes allopolyploids.

Chromosomal Differentiation of Deschampsia (Poaceae) Based on Four Satellite DNA Families

Diverse families of satellite DNA (satDNA) were detected in heterochromatin regions of Deschampsia. This kind of repetitive DNA consists of tandem repeat sequences forming big arrays in genomes, and can contribute to lineages differentiation. The differentiation between types of satDNA is related to their sequence identity, the size and number of monomers forming the array, and their chromosomal location. In this work, four families of satDNA (D2, D3, D12, D13), previously isolated by genomic analysis, were studied on chromosomal preparations of 12 species of Deschampsia (D. airiformis, D. antarctica, D. cespitosa, D. cordillerarum, D. elongata, D. kingii, D. laxa, D. mendocina, D. parvula, D. patula, D. venustula, and Deschampsia sp) and one of Deyeuxia (D. eminens). Despite the number of satDNA loci showing interspecific variation, the general distribution pattern of each satDNA family is maintained. The four satDNA families are AT-rich and associated with DAPI + heterochromatin regions. D2, D3, and D12 have mainly subterminal distribution, while D13 is distributed in intercalary regions. Such conservation of satDNA patterns suggests a not random distribution in genomes, where the variation between species is mainly associated with the array size and the loci number. The presence of satDNA in all species studied suggests a low genetic differentiation of sequences. On the other hand, the variation of the distribution pattern of satDNA has no clear association with phylogeny. This may be related to high differential amplification and contraction of sequences between lineages, as explained by the library model.

The Genomics of Plant Satellite DNA

The twenty-first century began with a certain indifference to the research of satellite DNA (satDNA). Neither genome sequencing projects were able to accurately encompass the study of satDNA nor classic methodologies were able to go further in undertaking a better comprehensive study of the whole set of satDNA sequences of a genome. Nonetheless, knowledge of satDNA has progressively advanced during this century with the advent of new analytical techniques. The enormous advantages that genome-wide approaches have brought to its analysis have now stimulated a renewed interest in the study of satDNA. At this point, we can look back and try to assess more accurately many of the key questions that were left unsolved in the past about this enigmatic and important component of the genome. I review here the understanding gathered on plant satDNAs over the last few decades with an eye on the near future.

First Description of a Satellite DNA in Manatees' Centromeric Regions

Trichechus manatus and Trichechus inunguis are the two Sirenia species that occur in the Americas. Despite their increasing extinction risk, many aspects of their biology remain understudied, including the repetitive DNA fraction of their genomes. Here we used the sequenced genome of T. manatus and TAREAN to identify satellite DNAs (satDNAs) in this species. We report the first description of TMAsat, a satDNA comprising ~0.87% of the genome, with ~684bp monomers and centromeric localization. In T. inunguis, TMAsat showed similar monomer length, chromosome localization and conserved CENP-B box-like motifs as in T. manatus. We also detected this satDNA in the Dugong dugon and in the now extinct Hydrodamalis gigas genomes. The neighbor-joining tree shows that TMAsat sequences from T. manatus, T. inunguis, D. dugon, and H. gigas lack species-specific clusters, which disagrees with the predictions of concerted evolution. We detected a divergent TMAsat-like homologous sequence in elephants and hyraxes, but not in other mammals, suggesting this sequence was already present in the common ancestor of Paenungulata, and later became a satDNA in the Sirenians. This is the first description of a centromeric satDNA in manatees and will facilitate the inclusion of Sirenia in future studies of centromeres and satDNA biology.

Comparative Repeat Profiling of Two Closely Related Conifers (Larix decidua and Larix kaempferi) Reveals High Genome Similarity With Only Few Fast-Evolving Satellite DNAs

In eukaryotic genomes, cycles of repeat expansion and removal lead to large-scale genomic changes and propel organisms forward in evolution. However, in conifers, active repeat removal is thought to be limited, leading to expansions of their genomes, mostly exceeding 10 giga base pairs. As a result, conifer genomes are largely littered with fragmented and decayed repeats. Here, we aim to investigate how the repeat landscapes of two related conifers have diverged, given the conifers' accumulative genome evolution mode. For this, we applied low-coverage sequencing and read clustering to the genomes of European and Japanese larch, Larix decidua (Lamb.) Carrière and Larix kaempferi (Mill.), that arose from a common ancestor, but are now geographically isolated. We found that both Larix species harbored largely similar repeat landscapes, especially regarding the transposable element content. To pin down possible genomic changes, we focused on the repeat class with the fastest sequence turnover: satellite DNAs (satDNAs). Using comparative bioinformatics, Southern, and fluorescent in situ hybridization, we reveal the satDNAs' organizational patterns, their abundances, and chromosomal locations. Four out of the five identified satDNAs are widespread in the Larix genus, with two even present in the more distantly related Pseudotsuga and Abies genera. Unexpectedly, the EulaSat3 family was restricted to L. decidua and absent from L. kaempferi, indicating its evolutionarily young age. Taken together, our results exemplify how the accumulative genome evolution of conifers may limit the overall divergence of repeats after speciation, producing only few repeat-induced genomic novelties.

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.

Sequence, Chromatin and Evolution of Satellite DNA

Satellite DNA consists of abundant tandem repeats that play important roles in cellular processes, including chromosome segregation, genome organization and chromosome end protection. Most satellite DNA repeat units are either of nucleosomal length or 5–10 bp long and occupy centromeric, pericentromeric or telomeric regions. Due to high repetitiveness, satellite DNA sequences have largely been absent from genome assemblies. Although few conserved satellite-specific sequence motifs have been identified, DNA curvature, dyad symmetries and inverted repeats are features of various satellite DNAs in several organisms. Satellite DNA sequences are either embedded in highly compact gene-poor heterochromatin or specialized chromatin that is distinct from euchromatin. Nevertheless, some satellite DNAs are transcribed into non-coding RNAs that may play important roles in satellite DNA function. Intriguingly, satellite DNAs are among the most rapidly evolving genomic elements, such that a large fraction is species-specific in most organisms. Here we describe the different classes of satellite DNA sequences, their satellite-specific chromatin features, and how these features may contribute to satellite DNA biology and evolution. We also discuss how the evolution of functional satellite DNA classes may contribute to speciation in plants and animals.

Origin and evolution of highly polymorphic rDNA sites in Alstroemeria longistaminea (Alstroemeriaceae) and related species

Alstroemeria (Alstroemeriaceae) displays a conserved and highly asymmetric karyotype, where most rDNA sites can be properly recognized by the size and morphology of the chromosomes. We analyzed the intraspecific variation of rDNA sites in A. longistaminea and compared with their distribution in other species (A. caryophyllaea and A. piauhyensis) and a representative of a sister genus, Bomarea edulis. All three species of Alstroemeria presented 2n = 16, and one to six B chromosomes were found in some individuals of A. longistaminea. There was a set of 12 conserved rDNA sites (four 5S and eight 35S) and up to 11 variable sites. B chromosomes were almost entirely covered by 35S signals, coupled with tiny 5S sites. Noteworthy, most rDNA sites found in A. caryophyllaea and A. piauhyensis were localized in chromosome positions similar to those in A. longistaminea, suggesting the existence of conserved hotspots for rDNA accumulation. Some of these hotspots were absent in Chilean Alstromeria as well in B. edulis. We propose that insertions of rDNA sequences on chromosomes do not occur randomly but rather on preferential sites or hotspots for insertions. The maintenance of these arrays, however, may be favored/constrained by different factors, resulting in stable or polymorphic sites.

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.

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.

Satellite DNA in Neotropical Deer Species

The taxonomy and phylogenetics of Neotropical deer have been mostly based on morphological criteria and needs a critical revision on the basis of new molecular and cytogenetic markers. In this study, we used the variation in the sequence, copy number, and chromosome localization of satellite I-IV DNA to evaluate evolutionary relationships among eight Neotropical deer species. Using FISH with satI-IV probes derived from Mazama gouazoubira, we proved the presence of satellite DNA blocks in peri/centromeric regions of all analyzed deer. Satellite DNA was also detected in the interstitial chromosome regions of species of the genus Mazama with highly reduced chromosome numbers. In contrast to Blastocerus dichotomus, Ozotoceros bezoarticus, and Odocoileus virginianus, Mazama species showed high abundance of satIV DNA by FISH. The phylogenetic analysis of the satellite DNA showed close relationships between O. bezoarticus and B. dichotomus. Furthermore, the Neotropical and Nearctic populations of O. virginianus formed a single clade. However, the satellite DNA phylogeny did not allow resolving the relationships within the genus Mazama. The high abundance of the satellite DNA in centromeres probably contributes to the formation of chromosomal rearrangements, thus leading to a fast and ongoing speciation in this genus, which has not yet been reflected in the satellite DNA sequence diversification.

Comparative analysis of morabine grasshopper genomes reveals highly abundant transposable elements and rapidly proliferating satellite DNA repeats

BACKGROUND

Repetitive DNA sequences, including transposable elements (TEs) and tandemly repeated satellite DNA (satDNAs), collectively called the "repeatome", are found in high proportion in organisms across the Tree of Life. Grasshoppers have large genomes, averaging 9 Gb, that contain a high proportion of repetitive DNA, which has hampered progress in assembling reference genomes. Here we combined linked-read genomics with transcriptomics to assemble, characterize, and compare the structure of repetitive DNA sequences in four chromosomal races of the morabine grasshopper Vandiemenella viatica species complex and determine their contribution to genome evolution.

RESULTS

We obtained linked-read genome assemblies of 2.73-3.27 Gb from estimated genome sizes of 4.26-5.07 Gb DNA per haploid genome of the four chromosomal races of V. viatica. These constitute the third largest insect genomes assembled so far. Combining complementary annotation tools and manual curation, we found a large diversity of TEs and satDNAs, constituting 66 to 75% per genome assembly. A comparison of sequence divergence within the TE classes revealed massive accumulation of recent TEs in all four races (314-463 Mb per assembly), indicating that their large genome sizes are likely due to similar rates of TE accumulation. Transcriptome sequencing showed more biased TE expression in reproductive tissues than somatic tissues, implying permissive transcription in gametogenesis. Out of 129 satDNA families, 102 satDNA families were shared among the four chromosomal races, which likely represent a diversity of satDNA families in the ancestor of the V. viatica chromosomal races. Notably, 50 of these shared satDNA families underwent differential proliferation since the recent diversification of the V. viatica species complex.

CONCLUSION

This in-depth annotation of the repeatome in morabine grasshoppers provided new insights into the genome evolution of Orthoptera. Our TEs analysis revealed a massive recent accumulation of TEs equivalent to the size of entire Drosophila genomes, which likely explains the large genome sizes in grasshoppers. Despite an overall high similarity of the TE and satDNA diversity between races, the patterns of TE expression and satDNA proliferation suggest rapid evolution of grasshopper genomes on recent timescales.