TAREAN: a computational tool for identification and characterization of satellite DNA from unassembled short reads

Variants of a major DNA satellite discriminate parental subgenomes in a hybrid parthenogenetic lizard Darevskia unisexualis (Darevsky, 1966)

Hybrid parthenogenetic animals are an exceptionally interesting model for studying the mechanisms and evolution of sexual and asexual reproduction. A diploid parthenogenetic lizard Darevskia unisexualis is a result of an ancestral cross between a maternal species Darevskia raddei nairensis and a paternal species Darevskia valentini and presents a unique opportunity for a cytogenetic and computational analysis of a hybrid karyotype. Our previous results demonstrated a significant divergence between the pericentromeric DNA sequences of the parental Darevskia species; however, an in-depth comparative study of their pericentromeres is still lacking. Here, using target sequencing of microdissected pericentromeric regions, we reveal and compare the repertoires of the pericentromeric tandem repeats of the parental Darevskia lizards. We found species-specific sequences of the major pericentromeric tandem repeat CLsat, which allowed computational prediction and experimental validation of fluorescent DNA probes discriminating parental chromosomes within the hybrid karyotype of D. unisexualis. Moreover, we have implemented a generalizable computational method, based on the optimization of the Levenshtein distance between tandem repeat monomers, for finding species-specific fluorescent probes for pericentromere staining. In total, we anticipate that our comparative analysis of Darevskia pericentromeric repeats, the species-specific fluorescent probes that we found and the pipeline that we developed will form a basis for the future detailed cytogenomic studies of a wide range of natural and laboratory hybrids.

Genome-wide sequence divergence of satellite DNA could underlie meiotic failure in male hybrids of bighead catfish and North African catfish (Clarias, Clariidae)

Hybrid sterility, a hallmark of postzygotic isolation, arises from parental genome divergence disrupting meiosis. While chromosomal incompatibility is often implicated, the underlying mechanisms remain unclear. This study investigated meiotic behavior and genome-wide divergence in bighead catfish (C. macrocephalus), North African catfish (C. gariepinus), and their sterile male hybrids (important in aquaculture). Repetitive DNA analysis using bioinformatics and cytogenetics revealed significant divergence in satellite DNA (satDNA) families between parental species. Notably, one hybrid exhibited successful meiosis and spermatozoa production, suggesting potential variation in sterility expression. Our findings suggest that genome-wide satDNA divergence, rather than chromosome number differences, likely contributes to meiotic failure and male sterility in these catfish hybrids.

The holocentricity in the dioecious nutmeg (Myristica fragrans) is not based on major satellite repeats

Holocentric species are characterized by the presence of centromeres throughout the length of the chromosomes. We confirmed the holocentricity of the dioecious, small chromosome-size species Myristica fragrans based on the chromosome-wide distribution of the centromere-specific protein KNL1, α-tubulin fibers, and the cell cycle-dependent histone H3 serine 28 phosphorylation (H3S28ph) mark. Each holocentromere is likely composed of, on average, ten centromere units, but none of the identified and in situ hybridized high-copy satellite repeats is centromere-specific. No sex-specific major repeats are present in the high-copy repeat composition of male or female plants, or a significant difference in genome size was detected. Therefore, it is unlikely that M. fragrans possesses heteromorphic sex chromosomes.

Reinforcement of repressive marks in the chicken primordial germ cell epigenetic signature: divergence from basal state resetting in mammals

BACKGROUND

In mammals, primordial germ cells (PGCs), the embryonic precursors of the germline, arise from embryonic or extra-embryonic cells upon induction by the surrounding tissues during gastrulation, according to mechanisms which are elucidated in mice but remain controversial in primates. They undergo genome-wide epigenetic reprogramming, consisting of extensive DNA demethylation and histone post-translational modification (PTM) changes, toward a basal, euchromatinized state. In contrast, chicken PGCs are specified by preformation before gastrulation based on maternally-inherited factors. They can be isolated from the bloodstream during their migration to the genital ridges. Our prior research highlighted differences in the global epigenetic profile of cultured chicken PGCs compared with chicken somatic cells and mammalian PGCs. This study investigates the acquisition and evolution of this profile during development.

RESULTS

Quantitative analysis of global DNA methylation and histone PTMs, including their distribution, during key stages of chicken early development revealed divergent PGC epigenetic changes compared with mammals. Unlike mammalian PGCs, chicken PGCs do not undergo genome-wide DNA demethylation or exhibit a decrease in histone H3 lysine 9 dimethylation. However, chicken PGCs show 5‑hydroxymethylcytosine loss, macroH2A redistribution, and chromatin decompaction, mirroring mammalian processes. Chicken PGCs initiate their epigenetic signature during migration, progressively accumulating high global levels of H3K9me3, with preferential enrichment in inactive genome regions. Despite apparent global chromatin decompaction, abundant heterochromatin marks, including repressive histone PTMs, HP1 variants, and DNA methylation, persists in chicken PGCs, contrasting with mammalian PGCs.

CONCLUSIONS

Chicken PGCs' epigenetic signature does not align with the basal chromatin state observed in mammals, suggesting a departure from extensive epigenetic reprogramming. Despite disparities in early PGC development, the persistence of several epigenetic features shared with mammals implies their involvement in chromatin-regulated germ cell properties, with the distinctive elevation of chicken-specific H3K9me3 potentially participating in these processes.

Ancient hybridization and repetitive element proliferation in the evolutionary history of the monocot genus Amomum (Zingiberaceae)

Genome size variation is a crucial aspect of plant evolution, influenced by a complex interplay of factors. Repetitive elements, which are fundamental components of genomic architecture, often play a role in genome expansion by selectively amplifying specific repeat motifs. This study focuses on Amomum, a genus in the ginger family (Zingiberaceae), known for its 4.4-fold variation in genome size. Using a robust methodology involving PhyloNet reconstruction, RepeatExplorer clustering, and repeat similarity-based phylogenetic network construction, we investigated the repeatome composition, analyzed repeat dynamics, and identified potential hybridization events within the genus. Our analysis confirmed the presence of four major infrageneric clades (A-D) within Amomum, with clades A-C exclusively comprising diploid species (2n = 48) and clade D encompassing both diploid and tetraploid species (2n = 48 and 96). We observed an increase in the repeat content within the genus, ranging from 84% to 89%, compared to outgroup species with 75% of the repeatome. The SIRE lineage of the Ty1-Copia repeat superfamily was prevalent in most analyzed ingroup genomes. We identified significant difference in repeatome structure between the basal Amomum clades (A, B, C) and the most diverged clade D. Our investigation revealed evidence of ancient hybridization events within Amomum, coinciding with a substantial proliferation of multiple repeat groups. This finding supports the hypothesis that ancient hybridization is a driving force in the genomic evolution of Amomum. Furthermore, we contextualize our findings within the broader context of genome size variations and repeatome dynamics observed across major monocot lineages. This study enhances our understanding of evolutionary processes within monocots by highlighting the crucial roles of repetitive elements in shaping genome size and suggesting the mechanisms that drive these changes.

Revealing the Satellite DNA Content in Ancistrus sp. (Siluriformes: Loricariidae) by Genomic and Bioinformatic Analysis

INTRODUCTION

Eukaryotic genomes are composed of simple, repetitive sequences, including satellite DNAs (satDNA), which are noncoding sequences arranged in tandem arrays. These sequences play a crucial role in genomic functions and innovations, influencing processes such as the maintenance of nuclear material, the formation of heterochromatin and the differentiation of sex chromosomes. In this genomic era, advances in next-generation sequencing and bioinformatics tools have facilitated the exhaustive cataloging of repetitive elements in genomes, particularly in non-model species. This study focuses on the satDNA content of Ancistrus sp., a diverse species of fish from the Loricariidae family. The genus Ancistrus shows significant karyotypic evolution, with extensive variability from the ancestral diploid number.

METHODS

By means of bioinformatic approaches, 40 satDNA families in Ancistrus sp., constituting 5.19% of the genome were identified. Analysis of the abundance and divergence landscape revealed diverse profiles, indicating recent amplification and homogenization of these satDNA sequences.

RESULTS

The most abundant satellite, AnSat1-142, constitutes 2.1% of the genome, while the least abundant, AnSat40-52, represents 0.0034%. The length of the monomer repeat varies from 16 to 142 base pairs, with an average length of 61 bp. These results contribute to understanding the genomic dynamics and evolution of satDNAs in Ancistrus sp.

CONCLUSION

The study underscores the variability of satDNAs between fish species and provides valuable information on chromosome organization and the evolution of repetitive elements in non-model organisms.

Phased Assembly of Neo-Sex Chromosomes Reveals Extensive Y Degeneration and Rapid Genome Evolution in Rumex hastatulus

Y chromosomes are thought to undergo progressive degeneration due to stepwise loss of recombination and subsequent reduction in selection efficiency. However, the timescales and evolutionary forces driving degeneration remain unclear. To investigate the evolution of sex chromosomes on multiple timescales, we generated a high-quality phased genome assembly of the massive older (<10 MYA) and neo (<200,000 yr) sex chromosomes in the XYY cytotype of the dioecious plant Rumex hastatulus and a hermaphroditic outgroup Rumex salicifolius. Our assemblies, supported by fluorescence in situ hybridization, confirmed that the neo-sex chromosomes were formed by two key events: an X-autosome fusion and a reciprocal translocation between the homologous autosome and the Y chromosome. The enormous sex-linked regions of the X (296 Mb) and two Y chromosomes (503 Mb) both evolved from large repeat-rich genomic regions with low recombination; however, the complete loss of recombination on the Y still led to over 30% gene loss and major rearrangements. In the older sex-linked region, there has been a significant increase in transposable element abundance, even into and near genes. In the neo-sex-linked regions, we observed evidence of extensive rearrangements without gene degeneration and loss. Overall, we inferred significant degeneration during the first 10 million years of Y chromosome evolution but not on very short timescales. Our results indicate that even when sex chromosomes emerge from repetitive regions of already-low recombination, the complete loss of recombination on the Y chromosome still leads to a substantial increase in repetitive element content and gene degeneration.

Satellite DNAs, heterochromatin, and sex chromosomes of the wattled jacana (Charadriiformes; Jacanidae): a species with highly rearranged karyotype

Charadriiformes, which comprises shorebirds and their relatives, is one of the most diverse avian orders, with over 390 species showing a wide range of karyotypes. Here, we isolated and characterized the whole collection of satellite DNAs (satDNAs) at both molecular and cytogenetic levels of one of its representative species, named the wattled jacana (Jacana jacana), a species that contains a typical ZZ/ZW sex chromosome system and a highly rearranged karyotype. In addition, we also investigate the in situ location of telomeric and microsatellite repeats. A small catalog of 11 satDNAs was identified that typically accumulated on microchromosomes and on the W chromosome. The latter also showed a significant accumulation of telomeric signals, being (GA)10 the only microsatellite with positive hybridization signals among all the 16 tested ones. These current findings contribute to our understanding of the genomic organization of repetitive DNAs in a bird species with high degree of chromosomal reorganization contrary to the majority of bird species that have stable karyotypes.

Cytogenomic Characterization of Transposable Elements and Satellite DNA in Passiflora L. Species

The species Passiflora alata, P. cincinnata, and P. edulis have great economic value due to the use of their fruits for human consumption. In this study, we compared the repetitive genome fractions of these three species. The compositions of the repetitive DNA of these three species' genomes were analyzed using clustering and identification of the repetitive sequences with RepeatExplorer. It was found that repetitive DNA content represents 74.70%, 66.86%, and 62.24% of the genome of P. alata, P. edulis, and P. cincinnata, respectively. LTR Ty3/Gypsy retrotransposons represent the highest genome proportions in P. alata and P. edulis, while Ty1/Copia comprises the largest proportion of P. cincinnata genome. Chromosomal mapping by Fluorescent In Situ Hybridization (FISH) showed that LTR retrotransposons have a dispersed distribution along chromosomes. The subtelomeric region of chromosomes is where 145 bp satellite DNA is located, suggesting that these elements may play important roles in genome structure and organization in these species. In this work, we obtained the first global characterization of the composition of repetitive DNA in Passiflora, showing that an increase in genome size is related to an increase in repetitive DNA, which represents an important evolutionary route for these species.

Heterochromatin Is Not the Only Place for satDNAs: The High Diversity of satDNAs in the Euchromatin of the Beetle Chrysolina americana (Coleoptera, Chrysomelidae)

The satellitome of the beetle Chrysolina americana Linneo, 1758 has been characterized through chromosomal analysis, genomic sequencing, and bioinformatics tools. C-banding reveals the presence of constitutive heterochromatin blocks enriched in A+T content, primarily located in pericentromeric regions. Furthermore, a comprehensive satellitome analysis unveils the extensive diversity of satellite DNA families within the genome of C. americana. Using fluorescence in situ hybridization techniques and the innovative CHRISMAPP approach, we precisely map the localization of satDNA families on assembled chromosomes, providing insights into their organization and distribution patterns. Among the 165 identified satDNA families, only three of them exhibit a remarkable amplification and accumulation, forming large blocks predominantly in pericentromeric regions. In contrast, the remaining, less abundant satDNA families are dispersed throughout euchromatic regions, challenging the traditional association of satDNA with heterochromatin. Overall, our findings underscore the complexity of repetitive DNA elements in the genome of C. americana and emphasize the need for further exploration to elucidate their functional significance and evolutionary implications.

Divergence of 10 satellite repeats in Artemisia (Asteraceae: Anthemideae) based on sequential fluorescence in situ hybridization analysis: evidence for species identification and evolution

Artemisia is a large genus encompassing about 400 diverse species, many of which have considerable medicinal and ecological value. However, complex morphological information and variation in ploidy level and nuclear DNA content have presented challenges for evolution studies of this genus. Consequently, taxonomic inconsistencies within the genus persist, hindering the utilization of such large plant resources. Researchers have utilized satellite DNAs to aid in chromosome identification, species classification, and evolutionary studies due to their significant sequence and copy number variation between species and close relatives. In the present study, the RepeatExplorer2 pipeline was utilized to identify 10 satellite DNAs from three species (Artemisia annua, Artemisia vulgaris, Artemisia viridisquama), and fluorescence in situ hybridization confirmed their distribution on chromosomes in 24 species, including 19 Artemisia species with 5 outgroup species from Ajania and Chrysanthemum. Signals of satellite DNAs exhibited substantial differences between species. We obtained one genus-specific satellite from the sequences. Additionally, molecular cytogenetic maps were constructed for Artemisia vulgaris, Artemisia leucophylla, and Artemisia viridisquama. One species (Artemisia verbenacea) showed a FISH distribution pattern suggestive of an allotriploid origin. Heteromorphic FISH signals between homologous chromosomes in Artemisia plants were observed at a high level. Additionally, the relative relationships between species were discussed by comparing ideograms. The results of the present study provide new insights into the accurate identification and taxonomy of the Artemisia genus using molecular cytological methods.

How diverse a monocentric chromosome can be? Repeatome and centromeric organization of Juncus effusus (Juncaceae)

Juncus is the largest genus of Juncaceae and was considered holocentric for a long time. Recent findings, however, indicated that 11 species from different clades of the genus have monocentric chromosomes. Thus, the Juncus centromere organization and evolution need to be reassessed. We aimed to investigate the major repetitive DNA sequences of two accessions of Juncus effusus and its centromeric structure by employing whole-genome analyses, fluorescent in situ hybridization, CENH3 immunodetection, and chromatin immunoprecipitation sequencing. We showed that the repetitive fraction of the small J. effusus genome (~270 Mbp/1C) is mainly composed of Class I and Class II transposable elements (TEs) and satellite DNAs. Three identified satellite DNA families were mainly (peri)centromeric, with two being associated with the centromeric protein CENH3, but not strictly centromeric. Two types of centromere organization were discerned in J. effusus: type 1 was characterized by a single CENH3 domain enriched with JefSAT1-155 or JefSAT2-180, whereas type 2 showed multiple CENH3 domains interrupted by other satellites, TEs or genes. Furthermore, while type 1 centromeres showed a higher degree of satellite identity along the array, type 2 centromeres had less homogenized arrays along the multiple CENH3 domains per chromosome. Although the analyses confirmed the monocentric organization of J. effusus chromosomes, our data indicate a more dynamic arrangement of J. effusus centromeres than observed for other plant species, suggesting it may constitute a transient state between mono- and holocentricity.

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.

Insights into the genome of the 'Loco' Concholepas concholepas (Gastropoda: Muricidae) from low-coverage short-read sequencing: genome size, ploidy, transposable elements, nuclear RNA gene operon, mitochondrial genome, and phylogenetic placement in the family Muricidae

BACKGROUND

The Peruvian 'chanque' or Chilean 'loco' Concholepas concholepas is an economically, ecologically, and culturally important muricid gastropod heavily exploited by artisanal fisheries in the temperate southeastern Pacific Ocean. In this study, we have profited from a set of bioinformatics tools to recover important biological information of C. concholepas from low-coverage short-read NGS datasets. Specifically, we calculated the size of the nuclear genome, ploidy, and estimated transposable elements content using an in silico k-mer approach, we discovered, annotated, and quantified those transposable elements, we assembled and annotated the 45S rDNA RNA operon and mitochondrial genome, and we confirmed the phylogenetic position of C. concholepas within the muricid subfamily Rapaninae based on translated protein coding genes.

RESULTS

Using a k-mer approach, the haploid genome size estimated for the predicted diploid genome of C. concholepas varied between 1.83 Gbp (with kmer = 24) and 2.32 Gbp (with kmer = 36). Between half and two thirds of the nuclear genome of C. concholepas was composed of transposable elements. The most common transposable elements were classified as Long Interspersed Nuclear Elements and Short Interspersed Nuclear Elements, which were more abundant than DNA transposons, simple repeats, and Long Terminal Repeats. Less abundant repeat elements included Helitron mobile elements, 45S rRNA DNA, and Satellite DNA, among a few others.The 45S rRNA DNA operon of C. concholepas that encodes for the ssrRNA, 5.8S rRNA, and lsrRNA genes was assembled into a single contig 8,090 bp long. The assembled mitochondrial genome of C. concholepas is 15,449 bp long and encodes 13 protein coding genes, two ribosomal genes, and 22 transfer RNAs.

CONCLUSION

The information gained by this study will inform the assembly of a high quality nuclear genome for C. concholepas and will support bioprospecting and biomonitoring using environmental DNA to advance development of conservation and management plans in this overexploited marine snail.

Sex-biased gene content is associated with sex chromosome turnover in Danaini butterflies

Sex chromosomes play an outsized role in adaptation and speciation, and thus deserve particular attention in evolutionary genomics. In particular, fusions between sex chromosomes and autosomes can produce neo-sex chromosomes, which offer important insights into the evolutionary dynamics of sex chromosomes. Here, we investigate the evolutionary origin of the previously reported Danaus neo-sex chromosome within the tribe Danaini. We assembled and annotated genomes of Tirumala septentrionis (subtribe Danaina), Ideopsis similis (Amaurina), Idea leuconoe (Euploeina) and Lycorea halia (Itunina) and identified their Z-linked scaffolds. We found that the Danaus neo-sex chromosome resulting from the fusion between a Z chromosome and an autosome corresponding to the Melitaea cinxia chromosome (McChr) 21 arose in a common ancestor of Danaina, Amaurina and Euploina. We also identified two additional fusions as the W chromosome further fused with the synteny block McChr31 in I. similis and independent fusion occurred between ancestral Z chromosome and McChr12 in L. halia. We further tested a possible role of sexually antagonistic selection in sex chromosome turnover by analysing the genomic distribution of sex-biased genes in I. leuconoe and L. halia. The autosomes corresponding to McChr21 and McChr31 involved in the fusions are significantly enriched in female- and male-biased genes, respectively, which could have hypothetically facilitated fixation of the neo-sex chromosomes. This suggests a role of sexual antagonism in sex chromosome turnover in Lepidoptera. The neo-Z chromosomes of both I. leuconoe and L. halia appear fully compensated in somatic tissues, but the extent of dosage compensation for the ancestral Z varies across tissues and species.

The Genome of the Yellow Mealworm, Tenebrio molitor: It's Bigger Than You Think

BACKGROUND

Insects are a sustainable source of protein for human food and animal feed. We present a genome assembly, CRISPR gene editing, and life stage-specific transcriptomes for the yellow mealworm, Tenebrio molitor, one of the most intensively farmed insects worldwide.

METHODS

Long and short reads and long-range data were obtained from a T. molitor male pupa. Sequencing transcripts from 12 T. molitor life stages resulted in 279 million reads for gene prediction and genetic engineering. A unique plasmid delivery system containing guide RNAs targeting the eye color gene vermilion flanking the muscle actin gene promoter and EGFP marker was used in CRISPR/Cas9 transformation.

RESULTS

The assembly is approximately 53% of the genome size of 756.8 ± 9.6 Mb, measured using flow cytometry. Assembly was complicated by a satellitome of at least 11 highly conserved satDNAs occupying 28% of the genome. The injection of the plasmid into embryos resulted in knock-out of Tm vermilion and knock-in of EGFP.

CONCLUSIONS

The genome of T. molitor is longer than current assemblies (including ours) due to a substantial amount (26.5%) of only one highly abundant satellite DNA sequence. Genetic sequences and transformation tools for an insect important to the food and feed industries will promote the sustainable utilization of mealworms and other farmed insects.

Chromosomal dynamics in Senna: comparative PLOP-FISH analysis of tandem repeats and flow cytometric nuclear genome size estimations

Introduction

Tandem repeats (TRs) occur abundantly in plant genomes. They play essential roles that affect genome organization and evolution by inducing or generating chromosomal rearrangements such as duplications, deletions, inversions, and translocations. These impact gene expression and chromosome structure and even contribute to the emergence of new species.

Method

We investigated the effects of TRs on speciation in Senna genus by performing a comparative analysis using fluorescence in situ hybridization (FISH) with S. tora-specific TR probes. We examined the chromosomal distribution of these TRs and compared the genome sizes of seven Senna species (estimated using flow cytometry) to better understand their evolutionary relationships.

Results

Two (StoTR03_159 and StoTR04_55) of the nine studied TRs were not detected in any of the seven Senna species, whereas the remaining seven were found in all or some species with patterns that were similar to or contrasted with those of S. tora. Of these studies species, only S. angulata showed significant genome rearrangements and dysploid karyotypes resembling those of S. tora. The genome sizes varied among these species and did not positively correlate with chromosome number. Notably, S. angulata had the fewest chromosomes (2n = 22) but a relatively large genome size.

Discussion

These findings reveal the dynamics of TRs and provide a cytogenetic depiction of chromosomal rearrangements during speciation in Senna. To further elucidate the dynamics of repeat sequences in Senna, future studies must include related species and extensive repeatomic studies, including those on transposable elements.

Repeat-based phylogenomics shed light on unclear relationships in the monocentric genus Juncus L. (Juncaceae)

The repetitive fraction (repeatome) of eukaryotic genomes is diverse and usually fast evolving, being an important tool for clarify plant systematics. The genus Juncus L. comprises 332 species, karyotypically recognized by having holocentric chromosomes. However, four species were recently described as monocentric, yet our understanding of their genome evolution is largely masked by unclear phylogenetic relationships. Here, we reassess the current Juncus systematics using low-coverage genome skimming data of 33 taxa to construct repeats, nuclear rDNA and plastome-based phylogenetic hypothesis. Furthermore, we characterize the repeatome and chromosomal distribution of Juncus-specific centromeric repeats/CENH3 protein to test the monocentricity reach in the genus. Repeat-base phylogenies revealed topologies congruent with the rDNA tree, but not with the plastome tree. The incongruence between nuclear and plastome chloroplast dataset suggest an ancient hybridization in the divergence of Juncotypus and Tenageia sections 40 Myr ago. The phylogenetic resolution at section level was better fitted with the rDNA/repeat-based approaches, with the recognition of two monophyletic sections (Stygiopsis and Tenageia). We found specific repeatome trends for the main lineages, such as the higher abundances of TEs in the Caespitosi and Iridifolii + Ozophyllum clades. CENH3 immunostaining confirmed the monocentricity of Juncus, which can be a generic synapomorphy for the genus. The heterogeneity of the repeatomes, with high phylogenetic informativeness, identified here may be correlated with their ancient origin (56 Mya) and reveals the potential of comparative genomic analyses for understanding plant systematics and evolution.

Repositioning of centromere-associated repeats during karyotype evolution in Oryzias fishes

The karyotype, which is the number and shape of chromosomes, is a fundamental characteristic of all eukaryotes. Karyotypic changes play an important role in many aspects of evolutionary processes, including speciation. In organisms with monocentric chromosomes, it was previously thought that chromosome number changes were mainly caused by centric fusions and fissions, whereas chromosome shape changes, that is, changes in arm numbers, were mainly due to pericentric inversions. However, recent genomic and cytogenetic studies have revealed examples of alternative cases, such as tandem fusions and centromere repositioning, found in the karyotypic changes within and between species. Here, we employed comparative genomic approaches to investigate whether centromere repositioning occurred during karyotype evolution in medaka fishes. In the medaka family (Adrianichthyidae), the three phylogenetic groups differed substantially in their karyotypes. The Oryzias latipes species group has larger numbers of chromosome arms than the other groups, with most chromosomes being metacentric. The O. javanicus species group has similar numbers of chromosomes to the O. latipes species group, but smaller arm numbers, with most chromosomes being acrocentric. The O. celebensis species group has fewer chromosomes than the other two groups and several large metacentric chromosomes that were likely formed by chromosomal fusions. By comparing the genome assemblies of O. latipes, O. javanicus, and O. celebensis, we found that repositioning of centromere-associated repeats might be more common than simple pericentric inversion. Our results demonstrated that centromere repositioning may play a more important role in karyotype evolution than previously appreciated.

Centromere innovations within a mouse species

Mammalian centromeres direct faithful genetic inheritance and are typically characterized by regions of highly repetitive and rapidly evolving DNA. We focused on a mouse species, Mus pahari, that we found has evolved to house centromere-specifying centromere protein-A (CENP-A) nucleosomes at the nexus of a satellite repeat that we identified and termed π-satellite (π-sat), a small number of recruitment sites for CENP-B, and short stretches of perfect telomere repeats. One M. pahari chromosome, however, houses a radically divergent centromere harboring ~6 mega-base pairs of a homogenized π-sat-related repeat, π-satB, that contains >20,000 functional CENP-B boxes. There, CENP-B abundance promotes accumulation of microtubule-binding components of the kinetochore and a microtubule-destabilizing kinesin of the inner centromere. We propose that the balance of pro- and anti-microtubule binding by the new centromere is what permits it to segregate during cell division with high fidelity alongside the older ones whose sequence creates a markedly different molecular composition.

De novo reconstruction of satellite repeat units from sequence data

Satellite DNA are long tandemly repeating sequences in a genome and may be organized as high-order repeats (HORs). They are enriched in centromeres and are challenging to assemble. Existing algorithms for identifying satellite repeats either require the complete assembly of satellites or only work for simple repeat structures without HORs. Here we describe Satellite Repeat Finder (SRF), a new algorithm for reconstructing satellite repeat units and HORs from accurate reads or assemblies without prior knowledge on repeat structures. Applying SRF to real sequence data, we show that SRF could reconstruct known satellites in human and well-studied model organisms. We also find satellite repeats are pervasive in various other species, accounting for up to 12% of their genome contents but are often underrepresented in assemblies. With the rapid progress in genome sequencing, SRF will help the annotation of new genomes and the study of satellite DNA evolution even if such repeats are not fully assembled.

Genome-wide analysis of transposable elements and satellite DNA in Humulus scandens, a dioecious plant with XX/XY1Y2 chromosomes

Transposable elements (TEs) and satellite DNAs, two major categories of repetitive sequences, are expected to accumulate in non-recombining genome regions, including sex-linked regions, and contribute to sex chromosome evolution. The dioecious plant, Humulus scandens, can be used for studying the evolution of the XX/XY1Y2 sex chromosomes. In this study, we thoroughly examined the repetitive components of male and female H. scandens using next-generation sequencing data followed by bioinformatics analysis and florescence in situ hybridization (FISH). The H. scandens genome has a high overall repetitive sequence composition, 68.30% in the female and 66.78% in the male genome, with abundant long terminal repeat (LTR) retrotransposons (RTs), including more Ty3/Gypsy than Ty1/Copia elements, particularly two Ty3/Gypsy lineages, Tekay and Retand. Most LTR-RT lineages were found dispersed across the chromosomes, though CRM and Athila elements were predominately found within the centromeres and the pericentromeric regions. The Athila elements also showed clearly higher FISH signal intensities in the Y1 and Y2 chromosomes than in the X or autosomes. Three novel satellite DNAs were specifically distributed in the centromeric and/or telomeric regions, with markedly different distributions on the X, Y1, and Y2 chromosomes. Combined with FISH using satellite DNAs to stain chromosomes during meiotic diakinesis, we determined the synapsis pattern and distinguish pseudoautosomal regions (PARs). The results indicate that the XY1Y2 sex chromosomes of H. scandens might have originated from a centric fission event. This study improves our understanding of the repetitive sequence organization of H. scandens genome and provides a basis for further analysis of their chromosome evolution process.

Comparative analysis of repetitive DNA in dysploid and non-dysploid Phaseolus beans

Structural karyotype changes result from ectopic recombination events frequently associated with repetitive DNA. Although most Phaseolus species present relatively stable karyotypes with 2n = 22 chromosomes, the karyotypes of species of the Leptostachyus group show high rates of structural rearrangements, including a nested chromosome fusion that led to the dysploid chromosome number of the group (2n = 20). We examined the roles of repetitive landscapes in the rearrangements of species of the Leptostachyus group using genome-skimming data to characterize the repeatome in a range of Phaseolus species and compared them to species of that group (P. leptostachyus and P. macvaughii). LTR retrotransposons, especially the Ty3/gypsy lineage Chromovirus, were the most abundant elements in the genomes. Differences in the abundance of Tekay, Retand, and SIRE elements between P. macvaughii and P. leptostachyus were reflected in their total amounts of Ty3/gypsy and Ty1/copia. The satellite DNA fraction was the most divergent among the species, varying both in abundance and distribution, even between P. leptostachyus and P. macvaughii. The rapid turnover of repeats in the Leptostachyus group may be associated with the several rearrangements observed.

The nature and organization of satellite DNAs in Petunia hybrida, related, and ancestral genomes

Introduction

The garden petunia, Petunia hybrida (Solanaceae) is a fertile, diploid, annual hybrid species (2n=14) originating from P. axillaris and P. inflata 200 years ago. To understand the recent evolution of the P. hybrida genome, we examined tandemly repeated or satellite sequences using bioinformatic and molecular cytogenetic analysis.

Methods

Raw reads from available genomic assemblies and survey sequences of P. axillaris N (PaxiN), P. inflata S6, (PinfS6), P. hybrida (PhybR27) and the here sequenced P. parodii S7 (PparS7) were used for graph and k-mer based cluster analysis of TAREAN and RepeatExplorer. Analysis of repeat specific monomer lengths and sequence heterogeneity of the major tandem repeat families with more than 0.01% genome proportion were complemented by fluorescent in situ hybridization (FISH) using consensus sequences as probes to chromosomes of all four species.

Results

Seven repeat families, PSAT1, PSAT3, PSAT4, PSAT5 PSAT6, PSAT7 and PSAT8, shared high consensus sequence similarity and organisation between the four genomes. Additionally, many degenerate copies were present. FISH in P. hybrida and in the three wild petunias confirmed the bioinformatics data and gave corresponding signals on all or some chromosomes. PSAT1 is located at the ends of all chromosomes except the 45S rDNA bearing short arms of chromosomes II and III, and we classify it as a telomere associated sequence (TAS). It is the most abundant satellite repeat with over 300,000 copies, 0.2% of the genomes. PSAT3 and the variant PSAT7 are located adjacent to the centromere or mid-arm of one to three chromosome pairs. PSAT5 has a strong signal at the end of the short arm of chromosome III in P. axillaris and P.inflata, while in P. hybrida additional interstitial sites were present. PSAT6 is located at the centromeres of chromosomes II and III. PSAT4 and PSAT8 were found with only short arrays.

Discussion

These results demonstrate that (i) repeat families occupy distinct niches within chromosomes, (ii) they differ in the copy number, cluster organization and homogenization events, and that (iii) the recent genome hybridization in breeding P. hybrida preserved the chromosomal position of repeats but affected the copy number of repetitive DNA.

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

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

Repetitive DNA sequence detection and its role in the human genome

Repetitive DNA sequences playing critical roles in driving evolution, inducing variation, and regulating gene expression. In this review, we summarized the definition, arrangement, and structural characteristics of repeats. Besides, we introduced diverse biological functions of repeats and reviewed existing methods for automatic repeat detection, classification, and masking. Finally, we analyzed the type, structure, and regulation of repeats in the human genome and their role in the induction of complex diseases. We believe that this review will facilitate a comprehensive understanding of repeats and provide guidance for repeat annotation and in-depth exploration of its association with human diseases.

Abundance and Diversification of Repetitive Elements in Decapoda Genomes

Repetitive elements are a major component of DNA sequences due to their ability to propagate through the genome. Characterization of Metazoan repetitive profiles is improving; however, current pipelines fail to identify a significant proportion of divergent repeats in non-model organisms. The Decapoda order, for which repeat content analyses are largely lacking, is characterized by extremely variable genome sizes that suggest an important presence of repetitive elements. Here, we developed a new standardized pipeline to annotate repetitive elements in non-model organisms, which we applied to twenty Decapoda and six other Crustacea genomes. Using this new tool, we identified 10% more repetitive elements than standard pipelines. Repetitive elements were more abundant in Decapoda species than in other Crustacea, with a very large number of highly repeated satellite DNA families. Moreover, we demonstrated a high correlation between assembly size and transposable elements and different repeat dynamics between Dendrobranchiata and Reptantia. The patterns of repetitive elements largely reflect the phylogenetic relationships of Decapoda and the distinct evolutionary trajectories within Crustacea. In summary, our results highlight the impact of repetitive elements on genome evolution in Decapoda and the value of our novel annotation pipeline, which will provide a baseline for future comparative analyses.

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.

miR-430 microRNA Family in Fishes: Molecular Characterization and Evolution

The miR-430 microRNA family has been described in multiple fish species as one of the first microRNAs expressed by the zygote. It has been suggested that this family is implicated in maternal mRNA elimination, but may also play a role in steroidogenesis, sexual differentiation, and flatfish metamorphosis. The miR-430 sequences have been found in multiple-copy tandem clusters but evidence of their conservation outside of teleost fishes is scarce. In the present study, we have characterized the tandem repeats organization of these microRNAs in different fish species, both model and of interest in aquaculture. A phylogenetic analysis of this family has allowed us to identify that the miR-430 duplication, which took place before the Chondrostei and Neopterygii groups' divergence, has resulted in three variants ("a", "b", and "c"). According to our data, variant "b" is the most closely related to the ancestral sequence. Furthermore, we have detected isolated instances of the miR-430 repeat subunit in some species, which suggests that this microRNA family may be affected by DNA rearrangements. This study provides new data about the abundance, variability, and organization of the miR-430 family in fishes.

Holocentromeres can consist of merely a few megabase-sized satellite arrays

The centromere is the chromosome region where microtubules attach during cell division. In contrast to monocentric chromosomes with one centromere, holocentric species usually distribute hundreds of centromere units along the entire chromatid. We assembled the chromosome-scale reference genome and analyzed the holocentromere and (epi)genome organization of the lilioid Chionographis japonica. Remarkably, each of its holocentric chromatids consists of only 7 to 11 evenly spaced megabase-sized centromere-specific histone H3-positive units. These units contain satellite arrays of 23 and 28 bp-long monomers capable of forming palindromic structures. Like monocentric species, C. japonica forms clustered centromeres in chromocenters at interphase. In addition, the large-scale eu- and heterochromatin arrangement differs between C. japonica and other known holocentric species. Finally, using polymer simulations, we model the formation of prometaphase line-like holocentromeres from interphase centromere clusters. Our findings broaden the knowledge about centromere diversity, showing that holocentricity is not restricted to species with numerous and small centromere units.

Cytogenetics Meets Genomics: Cytotaxonomy and Genomic Relationships among Color Variants of the Asian Arowana Scleropages formosus

Scleropages formosus (Osteoglossiformes, Teleostei) represents one of the most valued ornamental fishes, yet it is critically endangered due to overexploitation and habitat destruction. This species encompasses three major color groups that naturally occur in allopatric populations, but the evolutionary and taxonomic relationships of S. formosus color varieties remain uncertain. Here, we utilized a range of molecular cytogenetic techniques to characterize the karyotypes of five S. formosus color phenotypes, which correspond to naturally occurring variants: the red ones (Super Red); the golden ones (Golden Crossback and Highback Golden); the green ones (Asian Green and Yellow Tail Silver). Additionally, we describe the satellitome of S. formosus (Highback Golden) by applying a high-throughput sequencing technology. All color phenotypes possessed the same karyotype structure 2n = 50 (8m/sm + 42st/a) and distribution of SatDNAs, but different chromosomal locations of rDNAs, which were involved in a chromosome size polymorphism. Our results show indications of population genetic structure and microstructure differences in karyotypes of the color phenotypes. However, the findings do not clearly back up the hypothesis that there are discrete lineages or evolutionary units among the color phenotypes of S. formosus, but another case of interspecific chromosome stasis cannot be excluded.

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.

Centromere Innovations Within a Mouse Species

Mammalian centromeres direct faithful genetic inheritance and are typically characterized by regions of highly repetitive and rapidly evolving DNA. We focused on a mouse species, Mus pahari, that we found has evolved to house centromere-specifying CENP-A nucleosomes at the nexus of a satellite repeat that we identified and term π-satellite (π-sat), a small number of recruitment sites for CENP-B, and short stretches of perfect telomere repeats. One M. pahari chromosome, however, houses a radically divergent centromere harboring ~6 Mbp of a homogenized π-sat-related repeat, π-satB, that contains >20,000 functional CENP-B boxes. There, CENP-B abundance drives accumulation of microtubule-binding components of the kinetochore, as well as a microtubule-destabilizing kinesin of the inner centromere. The balance of pro- and anti-microtubule-binding by the new centromere permits it to segregate during cell division with high fidelity alongside the older ones whose sequence creates a markedly different molecular composition.

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".

Analysis in Proceratophrys boiei genome illuminates the satellite DNA content in a frog from the Brazilian Atlantic forest

Satellite DNAs (satDNAs) are one of the most abundant elements in genomes. Characterized as tandemly organized sequences that can be amplified into multiple copies, mainly in heterochromatic regions. The frog P. boiei (2n = 22, ZZ♂/ZW♀) is found in the Brazilian Atlantic forest and has an atypical pattern of heterochromatin distribution when compared to other anuran amphibians, with large pericentromeric blocks on all chromosomes. In addition, females of Proceratophrys boiei have a metacentric sex chromosome W showing heterochromatin in all chromosomal extension. In this work, we performed high-throughput genomic, bioinformatic, and cytogenetic analyses to characterize the satellite DNA content (satellitome) in P. boiei, mainly due to high amount of C-positive heterochromatin and the highly heterochromatic W sex chromosome. After all the analyses, it is remarkable that the satellitome of P. boiei is composed of a high number of satDNA families (226), making P. boiei the frog species with the highest number of satellites described so far. Consistent with the observation of large centromeric C-positive heterochromatin blocks, the genome of P. boiei is enriched with high copy number of repetitive DNAs, with total satDNA abundance comprising 16.87% of the genome. We successfully mapped via Fluorescence in situ hybridization the two most abundant repeats in the genome, PboSat01-176 and PboSat02-192, highlighting the presence of certain satDNAs sequences in strategic chromosomal regions (e.g., centromere and pericentromeric region), which leads to their participation in crucial processes for genomic organization and maintenance. Our study reveals a great diversity of satellite repeats that are driving genomic organization in this frog species. The characterization and approaches regarding satDNAs in this species of frog allowed the confirmation of some insights from satellite biology and a possible relationship with the evolution of sex chromosomes, especially in anuran amphibians, including P. boiei, for which data were not available.

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.

Maps of Constitutive-Heterochromatin Distribution for Four Martes Species (Mustelidae, Carnivora, Mammalia) Show the Formative Role of Macrosatellite Repeats in Interspecific Variation of Chromosome Structure

Constitutive-heterochromatin placement in the genome affects chromosome structure by occupying centromeric areas and forming large blocks. To investigate the basis for heterochromatin variation in the genome, we chose a group of species with a conserved euchromatin part: the genus Martes [stone marten (M. foina, 2n = 38), sable (M. zibellina, 2n = 38), pine marten (M. martes, 2n = 38), and yellow-throated marten (M. flavigula, 2n = 40)]. We mined the stone marten genome for the most abundant tandem repeats and selected the top 11 macrosatellite repetitive sequences. Fluorescent in situ hybridization revealed distributions of the tandemly repeated sequences (macrosatellites, telomeric repeats, and ribosomal DNA). We next characterized the AT/GC content of constitutive heterochromatin by CDAG (Chromomycin A3-DAPI-after G-banding). The euchromatin conservatism was shown by comparative chromosome painting with stone marten probes in newly built maps of the sable and pine marten. Thus, for the four Martes species, we mapped three different types of tandemly repeated sequences critical for chromosome structure. Most macrosatellites are shared by the four species with individual patterns of amplification. Some macrosatellites are specific to a species, autosomes, or the X chromosome. The variation of core macrosatellites and their prevalence in a genome are responsible for the species-specific variation of the heterochromatic blocks.

Complex polyploid and hybrid species in an apomictic and sexual tropical forage grass group: genomic composition and evolution in Urochloa (Brachiaria) species

BACKGROUND AND AIMS

Diploid and polyploid Urochloa (including Brachiaria, Panicum and Megathyrsus species) C4 tropical forage grasses originating from Africa are important for food security and the environment, often being planted in marginal lands worldwide. We aimed to characterize the nature of their genomes, the repetitive DNA and the genome composition of polyploids, leading to a model of the evolutionary pathways within the group including many apomictic species.

METHODS

Some 362 forage grass accessions from international germplasm collections were studied, and ploidy was determined using an optimized flow cytometry method. Whole-genome survey sequencing and molecular cytogenetic analysis were used to identify chromosomes and genomes in Urochloa accessions belonging to the 'brizantha' and 'humidicola' agamic complexes and U. maxima.

KEY RESULTS

Genome structures are complex and variable, with multiple ploidies and genome compositions within the species, and no clear geographical patterns. Sequence analysis of nine diploid and polyploid accessions enabled identification of abundant genome-specific repetitive DNA motifs. In situ hybridization with a combination of repetitive DNA and genomic DNA probes identified evolutionary divergence and allowed us to discriminate the different genomes present in polyploids.

CONCLUSIONS

We suggest a new coherent nomenclature for the genomes present. We develop a model of evolution at the whole-genome level in diploid and polyploid accessions showing processes of grass evolution. We support the retention of narrow species concepts for Urochloa brizantha, U. decumbens and U. ruziziensis, and do not consider diploids and polyploids of single species as cytotypes. The results and model will be valuable in making rational choices of parents for new hybrids, assist in use of the germplasm for breeding and selection of Urochloa with improved sustainability and agronomic potential, and assist in measuring and conserving biodiversity in grasslands.

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.

Autopolyploid inheritance and a heterozygous reciprocal translocation shape chromosome genetic behavior in tetraploid blueberry (Vaccinium corymbosum)

Understanding chromosome recombination behavior in polyploidy species is key to advancing genetic discoveries. In blueberry, a tetraploid species, the line of evidences about its genetic behavior still remain poorly understood, owing to the inter-specific, and inter-ploidy admixture of its genome and lack of in depth genome-wide inheritance and comparative structural studies. Here we describe a new high-quality, phased, chromosome-scale genome of a diploid blueberry, clone W85. The genome was integrated with cytogenetics and high-density, genetic maps representing six tetraploid blueberry cultivars, harboring different levels of wild genome admixture, to uncover recombination behavior and structural genome divergence across tetraploid and wild diploid species. Analysis of chromosome inheritance and pairing demonstrated that tetraploid blueberry behaves as an autotetraploid with tetrasomic inheritance. Comparative analysis demonstrated the presence of a reciprocal, heterozygous, translocation spanning one homolog of chr-6 and one of chr-10 in the cultivar Draper. The translocation affects pairing and recombination of chromosomes 6 and 10. Besides the translocation detected in Draper, no other structural genomic divergences were detected across tetraploid cultivars and highly inter-crossable wild diploid species. These findings and resources will facilitate new genetic and comparative genomic studies in Vaccinium and the development of genomic assisted selection strategy for this crop.

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

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

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

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

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

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

Conservation of Major Satellite DNAs in Snake Heterochromatin

Repetitive DNA sequences constitute a sizeable portion of animal genomes, and tandemly organized satellite DNAs are a major part of them. They are usually located in constitutive heterochromatin clusters in or near the centromeres or telomeres, and less frequently in the interstitial parts of chromosome arms. They are also frequently accumulated in sex chromosomes. The function of these clusters is to sustain the architecture of the chromosomes and the nucleus, and to regulate chromosome behavior during mitosis and meiosis. The study of satellite DNA diversity is important for understanding sex chromosome evolution, interspecific hybridization, and speciation. In this work, we identified four satellite DNA families in the genomes of two snakes from different families: Daboia russelii (Viperidae) and Pantherophis guttatus (Colubridae) and determine their chromosomal localization. We found that one family is localized in the centromeres of both species, whereas the others form clusters in certain chromosomes or subsets of chromosomes. BLAST with snake genome assemblies showed the conservation of such clusters, as well as a subtle presence of the satellites in the interspersed manner outside the clusters. Overall, our results show high conservation of satellite DNA in snakes and confirm the "library" model of satellite DNA evolution.

Analysis of 5S rDNA Genomic Organization Through the RepeatExplorer2 Pipeline: A Simplified Protocol

The ribosomal RNA genes (rDNA) are universal genome components with a housekeeping function, given the crucial role of ribosomal RNA in the synthesis of ribosomes and thus for life-on-Earth. Therefore, their genomic organization is of considerable interest for biologists, in general. Ribosomal RNA genes have also been largely used to establish phylogenetic relationships, and to identify allopolyploid or homoploid hybridization.Here, we demonstrate how high-throughput sequencing data, through graph clustering implemented in RepeatExplorer2 pipeline ( https://repeatexplorer-elixir.cerit-sc.cz/galaxy/ ), can be helpful to decipher the genomic organization of 5S rRNA genes. We show that the linear shapes of cluster graphs are reminiscent to the linked organization of 5S and 35S rDNA (L-type arrangement) while the circular graphs correspond to their separate arrangement (S-type). We further present a simplified protocol based on the paper by (Garcia et al., Front Plant Sci 11:41, 2020) about the use of graph clustering of 5S rDNA homoeologs (S-type) to identify hybridization events in the species history. We found that the graph complexity (i.e., graph circularity in this case) is related to ploidy and genome complexity, with diploids typically showing circular-shaped graphs while allopolyploids and other interspecific hybrids display more complex graphs, with usually two or more interconnected loops representing intergenic spacers. When a three-genomic comparative clustering analysis from a given hybrid (homoploid/allopolyploid) and its putative progenitor species (diploids) is performed, it is possible to identify the corresponding homoeologous 5S rRNA gene families, and to elucidate the contribution of each putative parental genome to the 5S rDNA pool of the hybrid. Thus, the analysis of 5S rDNA cluster graphs by RepeatExplorer, together with information coming from other sources (e.g., morphology, cytogenetics) is a complementary approach for the determination of allopolyploid or homoploid hybridization and even ancient introgression events.

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.

qPCR as a Selective Tool for Cytogenetics

qPCR is widely used in quantitative studies of plant genomes and transcriptomes. In this article, this method is considered as an auxiliary step in the preparation and selection of markers for FISH analysis. Several cases from the authors' research on populations of the same species were reviewed, and a comparison of the closely related species, as well as the adaptation of the markers, based on satellite tandem repeats (TRs) using quantitative qPCR data was conducted. In the selected cases, TRs with contrast abundance were identified in the cases of the Dasypyrum, Thinopyrum and Aegilops species, and the transfer of TRs between the wheat and related species was demonstrated. TRs with intraspecific copy number variation were revealed in Thinopyrum ponticum and wheat-wheatgrass partial amphidiploids, and the TR showing predominant hybridization to the sea buckthorn Y chromosome was identified. Additionally, problems such as the absence of a reference gene for qPCR, and low-efficiency and self-complementary primers, were illustrated. In the cases considered here, the qPCR results clearly show high correlation with the subsequent results of the FISH analysis, which confirms the value of this method for cytogenetic studies.

Identification of Iguania Ancestral Syntenic Blocks and Putative Sex Chromosomes in the Veiled Chameleon (Chamaeleo calyptratus, Chamaeleonidae, Iguania)

The veiled chameleon (Chamaeleo calyptratus) is a typical member of the family Chamaeleonidae and a promising object for comparative cytogenetics and genomics. The karyotype of C. calyptratus differs from the putative ancestral chameleon karyotype (2n = 36) due to a smaller chromosome number (2n = 24) resulting from multiple chromosome fusions. The homomorphic sex chromosomes of an XX/XY system were described recently using male-specific RADseq markers. However, the chromosomal pair carrying these markers was not identified. Here we obtained chromosome-specific DNA libraries of C. calyptratus by chromosome flow sorting that were assigned by FISH and sequenced. Sequence comparison with three squamate reptiles reference genomes revealed the ancestral syntenic regions in the C. calyptratus chromosomes. We demonstrated that reducing the chromosome number in the C. calyptratus karyotype occurred through two fusions between microchromosomes and four fusions between micro-and macrochromosomes. PCR-assisted mapping of a previously described Y-specific marker indicates that chromosome 5 may be the sex chromosome pair. One of the chromosome 5 conserved synteny blocks shares homology with the ancestral pleurodont X chromosome, assuming parallelism in the evolution of sex chromosomes from two basal Iguania clades (pleurodonts and acrodonts). The comparative chromosome map produced here can serve as the foundation for future genome assembly of chameleons and vertebrate-wide comparative genomic studies.