High conservation of the differentially amplified MPA2 satellite DNA family in parthenogenetic root-knot nematodes

Unzipped genome assemblies of polyploid root-knot nematodes reveal unusual and clade-specific telomeric repeats

Using long-read sequencing, we assembled and unzipped the polyploid genomes of Meloidogyne incognita, M. javanica and M. arenaria, three of the most devastating plant-parasitic nematodes. We found the canonical nematode telomeric repeat to be missing in these and other Meloidogyne genomes. In addition, we find no evidence for the enzyme telomerase or for orthologs of C. elegans telomere-associated proteins, suggesting alternative lengthening of telomeres. Instead, analyzing our assembled genomes, we identify species-specific composite repeats enriched mostly at one extremity of contigs. These repeats are G-rich, oriented, and transcribed, similarly to canonical telomeric repeats. We confirm them as telomeric using fluorescent in situ hybridization. These repeats are mostly found at one single end of chromosomes in these species. The discovery of unusual and specific complex telomeric repeats opens a plethora of perspectives and highlights the evolutionary diversity of telomeres despite their central roles in senescence, aging, and chromosome integrity.

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.

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.

The Centromere Histone Is Conserved and Associated with Tandem Repeats Sharing a Conserved 19-bp Box in the Holocentromere of Meloidogyne Nematodes

Although centromeres have conserved function, centromere-specific histone H3 (CenH3) and centromeric DNA evolve rapidly. The centromere drive model explains this phenomenon as a consequence of the conflict between fast-evolving DNA and CenH3, suggesting asymmetry in female meiosis as a crucial factor. We characterized evolution of the CenH3 protein in three closely related, polyploid mitotic parthenogenetic species of the Meloidogyne incognita group, and in the distantly related meiotic parthenogen Meloidogyne hapla. We identified duplication of the CenH3 gene in a putative sexual ancestral Meloidogyne. We found that one CenH3 (αCenH3) remained conserved in all extant species, including in distant Meloidogyne hapla, whereas the other evolved rapidly and under positive selection into four different CenH3 variants. This pattern of CenH3 evolution in Meloidogyne species suggests the subspecialization of CenH3s in ancestral sexual species. Immunofluorescence performed on mitotic Meloidogyne incognita revealed a dominant role of αCenH3 on its centromere, whereas the other CenH3s have lost their function in mitosis. The observed αCenH3 chromosome distribution disclosed cluster-like centromeric organization. The ChIP-Seq analysis revealed that in M. incognita αCenH3-associated DNA dominantly comprises tandem repeats, composed of divergent monomers which share a completely conserved 19-bp long box. Conserved αCenH3-associated DNA is also confirmed in the related mitotic Meloidogyne incognita group species suggesting preservation of both centromere protein and DNA constituents. We hypothesize that the absence of centromere drive in mitosis might allow for CenH3 and its associated DNA to achieve an equilibrium in which they can persist for long periods of time.

Long-time evolution and highly dynamic satellite DNA in leptodactylid and hylodid frogs

BACKGROUND

Satellite DNA sequences are the most abundant components of heterochromatin and are repeated in tandem hundreds to thousands of times in the genome. However, the number of repeats of a specific satellite family can vary even between the genomes of related species or populations. The PcP190 satellite DNA family was identified in the genome of the leptodactylid frog Physalaemus cuvieri, which showed to be derived most likely from the 5S rDNA in an ancestral species. In this study, we investigate the presence of the PcP190 satellite DNA in several P. cuvieri populations and in four closely related species at the chromosomal and molecular level. Furthermore, we investigate the occurrence of this satellite DNA in the genomes of P. marmoratus as well as in representative species of the leptodactylid genus Leptodactylus (L. latrans) and the hylodid family (Crossodactylus gaudichaudii), all with the aim of investigating if the PcP190 satellite DNA presents or not a restricted distribution.

RESULTS

The PcP190 satellite DNA was detected in all the analyzed species. Some of them exhibited particular sequence differences, allowing the identification of species-specific groups of sequences, but in other species, the sequences were more conserved. However, in a general analysis, conserved and variable domains have been recognized within the PcP190 monomer. The chromosomal analysis performed on P. cuvieri populations and closely related species revealed high variability of the satellite DNA amount and its chromosomal location, which has always been coincident with regions of centromeric/pericentromeric heterochromatin.

CONCLUSION

The PcP190 satellite DNA was found in representatives of two families, Leptodactylidae and Hylodidae, indicating that these sequences are widely distributed and conserved in these frogs. There is a pattern of non-random variation within the repeating units, indicating interplay between stochastic events and selective pressure along the PcP190 sequences. Karyotypic differences involving the PcP190 satellite DNA prove to be highly dynamic on the chromosomes of the Physalaemus and its differential accumulation has contributed to the differentiation process of the Z and W sex chromosomes in P. ephippifer.

An ancient satellite DNA has maintained repetitive units of the original structure in most species of the living fossil plant genus Zamia

ZpS1 satellite DNA is specific to the genus Zamia and presents repetitive units organized as long arrays and also as very short arrays dispersed in the genome. We have characterized the structure of the ZpS1 repeats in 12 species representative of the whole geographic distribution of the genus. In most species, the clone most common sequences (cMCS) were so similar that a general most common sequence (GMCS) of the ZpS1 repetitive unit in the genus could be obtained. The few partial variations from the GMCS found in cMCS of some species correspond to variable positions present in most other species, as indicated by the clone consensus sequences (cCS). Two species have an additional species-specific variety of ZpS1 satellite. The dispersed repeats were found to contain more mutations than repeats from long arrays. Our results indicate that all or most species of Zamia inherited the ZpS1 satellite from a common ancestor in Miocene and have maintained repetitive units of the original structure till present. The features of ZpS1 satellite in the genus Zamia are poorly compatible with the model of concerted evolution, but they are perfectly consistent with a new model of satellite evolution based on experimental evidences indicating that a specific amplification-substitution repair mechanism maintains the homogeneity and stability of the repeats structure in each satellite DNA originally present in a species as long as the species exists.

Sequence variability of the MspI satellite DNA family of the pinewood nematode Bursaphelenchus xylophilus at different geographic scales

Tandemly repeated sequences known as satellite DNA (satDNA) generally exhibit complex evolutionary patterns of concerted evolution in which mutations are homogenized and fixed in a stochastic process of molecular drive. Here, the nucleotidic variability of the MspI satDNA family of the pinewood nematode Bursaphelenchus xylophilus is analyzed in order to understand the evolutionary dynamics of satDNA at the intraspecific level. A total of 425 MspI monomer units, either PCR-amplified from isolates of local (Peninsula of Setúbal, Portugal) or worldwide origin, or retrieved from the B. xylophilus genome sequence, were characterized and compared. Whatever their origin, sliding window analysis of sequence variability patterns among monomers revealed low, moderate and highly variant domains, indicating that variable levels of evolutionary constraint may act upon the entire monomers. The phylogenetic inference based on the different sets of MspI satDNA family for this species shows a broad polymorphism of the individual monomers, which were distributed into four main clusters. However, such clustering appeared independent from the geographic origin of the nematodes, and could not discriminate isolates or groups of geographically close isolates. Rather, the formation of different phylogenetic groups within this satDNA family suggests an a priori embodying of a set of diverging repeats from a common ancestor satDNA library, which have been differently amplified along the evolutionary pathway of this species. The present work improves knowledge on the evolutionary dynamics of satDNA at the intraspecific level, and provides new information on satDNA sequence variability among natural populations sampled at a local geographic scale.

Long-term conservation vs high sequence divergence: the case of an extraordinarily old satellite DNA in bivalve mollusks

The ubiquity of satellite DNA (satDNA) sequences has raised much controversy over the abundance of divergent monomer variants and the long-time nucleotide sequence stability observed for many satDNA families. In this work, we describe the satDNA BIV160, characterized in nine species of the three main bivalve clades (Protobranchia, Pteriomorphia and Heteroconchia). BIV160 monomers are similar in repeat size and nucleotide sequence to satDNAs described earlier in oysters and in the clam Donax trunculus. The broad distribution of BIV160 satDNA indicates that similar variants existed in the ancestral bivalve species that lived about 540 million years ago; this makes BIV160 the most ancient satDNA described so far. In the species examined, monomer variants are distributed in quite a complex pattern. This pattern includes (i) species characterized by a specific group of variants, (ii) species that share distinct group(s) of variants and (iii) species with both specific and shared types. The evolutionary scenario suggested by these data reconciles sequence uniformity in homogenization-maintained satDNA arrays with the genomic richness of divergent monomer variants formed by diversification of the same ancestral satDNA sequence. Diversified repeats can continue to evolve in a non-concerted manner and behave as independent amplification-contraction units in the framework of a 'library of satDNA variants' representing a permanent source of monomers that can be amplified into novel homogeneous satDNA arrays. On the whole, diversification of satDNA monomers and copy number fluctuations provide a highly dynamic genomic environment able to form and displace satDNA sequence variants rapidly in evolution.

Satellite DNAs between selfishness and functionality: structure, genomics and evolution of tandem repeats in centromeric (hetero)chromatin

Satellite DNAs (tandemly repeated, non-coding DNA sequences) stretch over almost all native centromeres and surrounding pericentromeric heterochromatin. Once considered as inert by-products of genome dynamics in heterochromatic regions, recent studies showed that satellite DNA evolution is interplay of stochastic events and selective pressure. This points to a functional significance of satellite sequences, which in (peri)centromeres may play some fundamental functional roles. First, specific interactions with DNA-binding proteins are proposed to complement sequence-independent epigenetic processes. The second role is achieved through RNAi mechanism, in which transcripts of satellite sequences initialize heterochromatin formation. In addition, satellite DNAs in (peri)centromeric regions affect chromosomal dynamics and genome plasticity. Paradoxically, while centromeric function is conserved through eukaryotes, the profile of satellite DNAs in this region is almost always species-specific. We argue that tandem repeats may be advantageous forms of DNA sequences in (peri)centromeres due to concerted evolution, which maintains high intra-array and intrapopulation sequence homogeneity of satellite arrays, while allowing rapid changes in nucleotide sequence and/or composition of satellite repeats. This feature may be crucial for long-term stability of DNA-protein interactions in centromeric regions.

Interplay of selective pressure and stochastic events directs evolution of the MEL172 satellite DNA library in root-knot nematodes

According to the library model, related species can have in common satellite DNA (satDNA) families amplified in differing abundances, but reasons for persistence of particular sequences in the library during long periods of time are poorly understood. In this paper, we characterize 3 related satDNAs coexisting in the form of a library in mitotic parthenogenetic root-knot nematodes of the genus Meloidogyne. Due to sequence similarity and conserved monomer length of 172 bp, this group of satDNAs is named MEL172. Analysis of sequence variability patterns among monomers of the 3 MEL172 satellites revealed 2 low-variable (LV) domains highly reluctant to sequence changes, 2 moderately variable (MV) domains characterized by limited number of mutations, and 1 highly variable (HV) domain. The latter domain is prone to rapid spread and homogenization of changes. Comparison of the 3 MEL172 consensus sequences shows that the LV domains have 6% changed nucleotide positions, the MV domains have 48%, whereas 78% divergence is concentrated in the HV domain. Conserved distribution of intersatellite variability might indicate a complex pattern of interactions in heterochromatin, which limits the range and phasing of allowed changes, implying a possible selection imposed on monomer sequences. The lack of fixed species-diagnostic mutations in each of the examined MEL172 satellites suggests that they existed in unaltered form in a common ancestor of extant species. Consequently, the evolution of these satellites seems to be driven by interplay of selective constraints and stochastic events. We propose that new satellites were derived from an ancestral progenitor sequence by nonrandom accumulation of mutations due to selective pressure on particular sequence segments. In the library of particular taxa, established satellites might be subject to differential amplification at chance due to stochastic mechanisms of concerted evolution.