Differential rates of local and global homogenization in centromere satellites from Arabidopsis relatives

Cycles of satellite and transposon evolution in Arabidopsis centromeres

Centromeres are critical for cell division, loading CENH3 or CENPA histone variant nucleosomes, directing kinetochore formation and allowing chromosome segregation^1,2. Despite their conserved function, centromere size and structure are diverse across species. To understand this centromere paradox^3,4, it is necessary to know how centromeric diversity is generated and whether it reflects ancient trans-species variation or, instead, rapid post-speciation divergence. To address these questions, we assembled 346 centromeres from 66 Arabidopsis thaliana and 2 Arabidopsis lyrata accessions, which exhibited a remarkable degree of intra- and inter-species diversity. A. thaliana centromere repeat arrays are embedded in linkage blocks, despite ongoing internal satellite turnover, consistent with roles for unidirectional gene conversion or unequal crossover between sister chromatids in sequence diversification. Additionally, centrophilic ATHILA transposons have recently invaded the satellite arrays. To counter ATHILA invasion, chromosome-specific bursts of satellite homogenization generate higher-order repeats and purge transposons, in line with cycles of repeat evolution. Centromeric sequence changes are even more extreme in comparison between A. thaliana and A. lyrata. Together, our findings identify rapid cycles of transposon invasion and purging through satellite homogenization, which drive centromere evolution and ultimately contribute to speciation.

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.

Hybrid seed incompatibility in Capsella is connected to chromatin condensation defects in the endosperm

Hybridization of closely related plant species is frequently connected to endosperm arrest and seed failure, for reasons that remain to be identified. In this study, we investigated the molecular events accompanying seed failure in hybrids of the closely related species pair Capsella rubella and C. grandiflora. Mapping of QTL for the underlying cause of hybrid incompatibility in Capsella identified three QTL that were close to pericentromeric regions. We investigated whether there are specific changes in heterochromatin associated with interspecific hybridizations and found a strong reduction of chromatin condensation in the endosperm, connected with a strong loss of CHG and CHH methylation and random loss of a single chromosome. Consistent with reduced DNA methylation in the hybrid endosperm, we found a disproportionate deregulation of genes located close to pericentromeric regions, suggesting that reduced DNA methylation allows access of transcription factors to targets located in heterochromatic regions. Since the identified QTL were also associated with pericentromeric regions, we propose that relaxation of heterochromatin in response to interspecies hybridization exposes and activates loci leading to hybrid seed failure.

Elucidating the major hidden genomic components of the A, C, and AC genomes and their influence on Brassica evolution

Decoding complete genome sequences is prerequisite for comprehensive genomics studies. However, the currently available reference genome sequences of Brassica rapa (A genome), B. oleracea (C) and B. napus (AC) cover 391, 540, and 850 Mbp and represent 80.6, 85.7, and 75.2% of the estimated genome size, respectively, while remained are hidden or unassembled due to highly repetitive nature of these genome components. Here, we performed the first comprehensive genome-wide analysis using low-coverage whole-genome sequences to explore the hidden genome components based on characterization of major repeat families in the B. rapa and B. oleracea genomes. Our analysis revealed 10 major repeats (MRs) including a new family comprising about 18.8, 10.8, and 11.5% of the A, C and AC genomes, respectively. Nevertheless, these 10 MRs represented less than 0.7% of each assembled reference genome. Genomic survey and molecular cytogenetic analyses validates our insilico analysis and also pointed to diversity, differential distribution, and evolutionary dynamics in the three Brassica species. Overall, our work elucidates hidden portions of three Brassica genomes, thus providing a resource for understanding the complete genome structures. Furthermore, we observed that asymmetrical accumulation of the major repeats might be a cause of diversification between the A and C genomes.

Centromere location in Arabidopsis is unaltered by extreme divergence in CENH3 protein sequence

During cell division, spindle fibers attach to chromosomes at centromeres. The DNA sequence at regional centromeres is fast evolving with no conserved genetic signature for centromere identity. Instead CENH3, a centromere-specific histone H3 variant, is the epigenetic signature that specifies centromere location across both plant and animal kingdoms. Paradoxically, CENH3 is also adaptively evolving. An ongoing question is whether CENH3 evolution is driven by a functional relationship with the underlying DNA sequence. Here, we demonstrate that despite extensive protein sequence divergence, CENH3 histones from distant species assemble centromeres on the same underlying DNA sequence. We first characterized the organization and diversity of centromere repeats in wild-type Arabidopsis thaliana. We show that A. thaliana CENH3-containing nucleosomes exhibit a strong preference for a unique subset of centromeric repeats. These sequences are largely missing from the genome assemblies and represent the youngest and most homogeneous class of repeats. Next, we tested the evolutionary specificity of this interaction in a background in which the native A. thaliana CENH3 is replaced with CENH3s from distant species. Strikingly, we find that CENH3 from Lepidium oleraceum and Zea mays, although specifying epigenetically weaker centromeres that result in genome elimination upon outcrossing, show a binding pattern on A. thaliana centromere repeats that is indistinguishable from the native CENH3. Our results demonstrate positional stability of a highly diverged CENH3 on independently evolved repeats, suggesting that the sequence specificity of centromeres is determined by a mechanism independent of CENH3.

Satellite DNA in Plants: More than Just Rubbish

For decades, satellite DNAs have been the hidden part of genomes. Initially considered as junk DNA, there is currently an increasing appreciation of the functional significance of satellite DNA repeats and of their sequences. Satellite DNA families accumulate in the heterochromatin in different parts of the eukaryotic chromosomes, mainly in pericentromeric and subtelomeric regions, but they also span the functional centromere. Tandem repeat sequences may spread from subtelomeric to interstitial loci, leading to the formation of chromosome-specific loci or to the accumulation in equilocal sites in different chromosomes. They also appear as the main components of the heterochromatin in the sex-specific region of sex chromosomes. Satellite DNA, required for chromosome organization, also plays a role in pairing and segregation. Some satellite repeats are transcribed and can participate in the formation and maintenance of heterochromatin structure and in the modulation of gene expression. In addition to the identification of the different satellite DNA families, their characteristics and location, we are interested in determining their impact on the genomes, by identifying the mechanisms leading to their appearance and amplification as well as in understanding how they change over time, the factors affecting these changes, and the influence exerted by the evolutionary history of the organisms. On the other hand, satellite DNA sequences are rapidly evolving sequences that may cause reproductive barriers between organisms and promote speciation. The accumulation of experimental data collected in recent years and the emergence of new approaches based on next-generation sequencing and high-throughput genome analysis are opening new perspectives that are changing our understanding of satellite DNA. This review examines recent data to provide a timely update on the overall information gathered about this part of the genome, focusing on the advances in the knowledge of its origin, its evolution, and its potential functional roles.

Diversity and evolution of centromere repeats in the maize genome

Centromere repeats are found in most eukaryotes and play a critical role in kinetochore formation. Though centromere repeats exhibit considerable diversity both within and among species, little is understood about the mechanisms that drive centromere repeat evolution. Here, we use maize as a model to investigate how a complex history involving polyploidy, fractionation, and recent domestication has impacted the diversity of the maize centromeric repeat CentC. We first validate the existence of long tandem arrays of repeats in maize and other taxa in the genus Zea. Although we find considerable sequence diversity among CentC copies genome-wide, genetic similarity among repeats is highest within these arrays, suggesting that tandem duplications are the primary mechanism for the generation of new copies. Nonetheless, clustering analyses identify similar sequences among distant repeats, and simulations suggest that this pattern may be due to homoplasious mutation. Although the two ancestral subgenomes of maize have contributed nearly equal numbers of centromeres, our analysis shows that the majority of all CentC repeats derive from one of the parental genomes, with an even stronger bias when examining the largest assembled contiguous clusters. Finally, by comparing maize with its wild progenitor teosinte, we find that the abundance of CentC likely decreased after domestication, while the pericentromeric repeat Cent4 has drastically increased.

Identification and characterization of functional centromeres of the common bean

In higher eukaryotes, centromeres are typically composed of megabase-sized arrays of satellite repeats that evolve rapidly and homogenize within a species' genome. Despite the importance of centromeres, our knowledge is limited to a few model species. We conducted a comprehensive analysis of common bean (Phaseolus vulgaris) centromeric satellite DNA using genomic data, fluorescence in situ hybridization (FISH), immunofluorescence and chromatin immunoprecipitation (ChIP). Two unrelated centromere-specific satellite repeats, CentPv1 and CentPv2, and the common bean centromere-specific histone H3 (PvCENH3) were identified. FISH showed that CentPv1 and CentPv2 are predominantly located at subsets of eight and three centromeres, respectively. Immunofluorescence- and ChIP-based assays demonstrated the functional significance of CentPv1 and CentPv2 at centromeres. Genomic analysis revealed several interesting features of CentPv1 and CentPv2: (i) CentPv1 is organized into an higher-order repeat structure, named Nazca, of 528 bp, whereas CentPv2 is composed of tandemly organized monomers; (ii) CentPv1 and CentPv2 have undergone chromosome-specific homogenization; and (iii) CentPv1 and CentPv2 are not likely to be commingled in the genome. These findings suggest that two distinct sets of centromere sequences have evolved independently within the common bean genome, and provide insight into centromere satellite evolution.

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.

Conserved DNA Motifs, Including the CENP-B Box-like, Are Possible Promoters of Satellite DNA Array Rearrangements in Nematodes

Tandemly arrayed non-coding sequences or satellite DNAs (satDNAs) are rapidly evolving segments of eukaryotic genomes, including the centromere, and may raise a genetic barrier that leads to speciation. However, determinants and mechanisms of satDNA sequence dynamics are only partially understood. Sequence analyses of a library of five satDNAs common to the root-knot nematodes Meloidogyne chitwoodi and M. fallax together with a satDNA, which is specific for M. chitwoodi only revealed low sequence identity (32-64%) among them. However, despite sequence differences, two conserved motifs were recovered. One of them turned out to be highly similar to the CENP-B box of human alpha satDNA, identical in 10-12 out of 17 nucleotides. In addition, organization of nematode satDNAs was comparable to that found in alpha satDNA of human and primates, characterized by monomers concurrently arranged in simple and higher-order repeat (HOR) arrays. In contrast to alpha satDNA, phylogenetic clustering of nematode satDNA monomers extracted either from simple or from HOR array indicated frequent shuffling between these two organizational forms. Comparison of homogeneous simple arrays and complex HORs composed of different satDNAs, enabled, for the first time, the identification of conserved motifs as obligatory components of monomer junctions. This observation highlights the role of short motifs in rearrangements, even among highly divergent sequences. Two mechanisms are proposed to be involved in this process, i.e., putative transposition-related cut-and-paste insertions and/or illegitimate recombination. Possibility for involvement of the nematode CENP-B box-like sequence in the transposition-related mechanism and together with previously established similarity of the human CENP-B protein and pogo-like transposases implicate a novel role of the CENP-B box and related sequence motifs in addition to the known function in centromere protein binding.

Tandem repeats derived from centromeric retrotransposons

Background

Tandem repeats are ubiquitous and abundant in higher eukaryotic genomes and constitute, along with transposable elements, much of DNA underlying centromeres and other heterochromatic domains. In maize, centromeric satellite repeat (CentC) and centromeric retrotransposons (CR), a class of Ty3/gypsy retrotransposons, are enriched at centromeres. Some satellite repeats have homology to retrotransposons and several mechanisms have been proposed to explain the expansion, contraction as well as homogenization of tandem repeats. However, the origin and evolution of tandem repeat loci remain largely unknown.

Results

CRM1TR and CRM4TR are novel tandem repeats that we show to be entirely derived from CR elements belonging to two different subfamilies, CRM1 and CRM4. Although these tandem repeats clearly originated in at least two separate events, they are derived from similar regions of their respective parent element, namely the long terminal repeat (LTR) and untranslated region (UTR). The 5^′ ends of the monomer repeat units of CRM1TR and CRM4TR map to different locations within their respective LTRs, while their 3^′ ends map to the same relative position within a conserved region of their UTRs. Based on the insertion times of heterologous retrotransposons that have inserted into these tandem repeats, amplification of the repeats is estimated to have begun at least ~4 (CRM1TR) and ~1 (CRM4TR) million years ago. Distinct CRM1TR sequence variants occupy the two CRM1TR loci, indicating that there is little or no movement of repeats between loci, even though they are separated by only ~1.4 Mb.

Conclusions

The discovery of two novel retrotransposon derived tandem repeats supports the conclusions from earlier studies that retrotransposons can give rise to tandem repeats in eukaryotic genomes. Analysis of monomers from two different CRM1TR loci shows that gene conversion is the major cause of sequence variation. We propose that successive intrastrand deletions generated the initial repeat structure, and gene conversions increased the size of each tandem repeat locus.

Concerted evolution of satellite DNA in Sarcocapnos: a matter of time

SarkOne is a genus-specific satellite-DNA family, isolated from the genomes of the species of the genus Sarcocapnos. This satellite DNA is composed of repeats with a consensus length of 855 bp and a mean G+C content of 52.5%. We have sequenced a total of 189 SarkOne monomeric repeats belonging to a total of seven species of the genus Sarcocapnos. The comparative analysis of these sequences both at the intraspecific and the interspecific levels have revealed divergence patterns between species are proportional to between-species divergence according to the phylogeny of the genus. Our study demonstrates that the molecular drive leading to the concerted-evolution pattern of this satellite DNA is a time-dependent process by which new mutations are spreading through genomes and populations at a gradual pace. However, time is a limiting factor in the observation of concerted evolution in some pairwise comparisons. Thus, pairwise comparisons of species sharing a recent common ancestor did not reveal nucleotide sites in transitional stages higher than stage III according to the Strachan's model. By contrast, there was a gradation in the percentage of upper transition stages (IV, V, VI) the more phylogenetically distant the species were. In addition, closely related species shared a high number of polymorphic sites, but these types of sites were not common when comparing more distant species. All these data are discussed in the light of current life-cycle models of satellite-DNA evolution.

Organization and molecular evolution of CENP-A--associated satellite DNA families in a basal primate genome

Centromeric regions in many complex eukaryotic species contain highly repetitive satellite DNAs. Despite the diversity of centromeric DNA sequences among species, the functional centromeres in all species studied to date are marked by CENP-A, a centromere-specific histone H3 variant. Although it is well established that families of multimeric higher-order alpha satellite are conserved at the centromeres of human and great ape chromosomes and that diverged monomeric alpha satellite is found in old and new world monkey genomes, little is known about the organization, function, and evolution of centromeric sequences in more distant primates, including lemurs. Aye-Aye (Daubentonia madagascariensis) is a basal primate and is located at a key position in the evolutionary tree to study centromeric satellite transitions in primate genomes. Using the approach of chromatin immunoprecipitation with antibodies directed to CENP-A, we have identified two satellite families, Daubentonia madagascariensis Aye-Aye 1 (DMA1) and Daubentonia madagascariensis Aye-Aye 2 (DMA2), related to each other but unrelated in sequence to alpha satellite or any other previously described primate or mammalian satellite DNA families. Here, we describe the initial genomic and phylogenetic organization of DMA1 and DMA2 and present evidence of higher-order repeats in Aye-Aye centromeric domains, providing an opportunity to study the emergence of chromosome-specific modes of satellite DNA evolution in primate genomes.

Characterization of satellite CentC repeats from heterochromatic regions on the long arm of maize B-chromosome

The B-chromosome of maize contains an A-chromosome centromere-specific satellite CentC repeat in its centromere region (CENB) and at multiple locations in its distal heterochromatic regions (BDHs). Because CentC is highly repetitive, it is a challenge to study CentC sequences within individual centromeres or chromosome regions. The combined structure of CentC and a BDH-specific CL-repeat has allowed us to isolate CentC sequences from BDHs. In the study described herein, we have used a PCR method to amplify 13 CL-CentC variant products that were specifically mapped to A-centromeres (CENAs), the CENB, and BDHs via the tertiary trisomes and hypoploids of five B-10L translocations. Cloning and sequence analyses of these CL-CentC products have revealed a local CentC homogenization within the three CentC-containing regions. Phylogenetic analysis has indicated that the CentC sequences of BDHs are more closely related to those of CENAs in comparison to that of the CENB. Furthermore, the CentC monomers that are within the CENB are more diverse than those within BDHs and CENAs. These results shed light on the evolution of CentC repeats on the B-chromosome and provide a better understanding of B-chromosome evolution.

Rapid divergence of repetitive DNAs in Brassica relatives

Centromeric, subtelomeric, and telomeric repetitive DNAs were characterized in Brassica species and the related Raphanus sativus and Arabidopsis thaliana. In general, rapid divergence of the repeats was found. The centromeric tandem satellite repeats were differentially distributed in the species studied, suggesting that centromeric repeats have diverged during the evolution of the A/C and B genome lineages. Sequence analysis of centromeric repeats suggested rapid evolution. Pericentromere-associated retrotransposons were identified and showed divergence during the evolution of the lineages as centromeric repeats. A novel subtelomeric tandem repeat from B. nigra was found to be conserved across the diploid Brassica genomes; however, this sequence was not identified in the related species. In contrast to previous studies, interstitial telomere-like repeats were identified in the pericentromeres of Brassica chromosomes, and these repeats may be associated with genomic stability. These results provide insight into genome evolution during polyploidization in Brassica and divergence within the Brassicaceae.

Fall and rise of satellite repeats in allopolyploids of Nicotiana over c. 5 million years

Allopolyploids represent natural experiments in which DNA sequences from different species are combined into a single nucleus and then coevolve, enabling us to follow the parental genomes, their interactions and evolution over time. Here, we examine the fate of satellite DNA over 5 million yr of divergence in plant genus Nicotiana (family Solanaceae). We isolated subtelomeric, tandemly repeated satellite DNA from Nicotiana diploid and allopolyploid species and analysed patterns of inheritance and divergence by sequence analysis, Southern blot hybridization and fluorescent in situ hybridization (FISH). We observed that parental satellite sequences redistribute around the genome in allopolyploids of Nicotiana section Polydicliae, formed c. 1 million yr ago (Mya), and that new satellite repeats evolved and amplified in section Repandae, which was formed c. 5 Mya. In some cases that process involved the complete replacement of parental satellite sequences. The rate of satellite repeat replacement is faster than theoretical predictions assuming the mechanism involved is unequal recombination and crossing-over. Instead we propose that this mechanism occurs with the deletion of large chromatin blocks and reamplification, perhaps via rolling circle replication.

Molecular and chromosomal evidence for allopolyploidy in soybean

Recent studies have documented that the soybean (Glycine max) genome has undergone two rounds of large-scale genome and/or segmental duplication. To shed light on the timing and nature of these duplication events, we characterized and analyzed two subfamilies of high-copy centromeric satellite repeats, CentGm-1 and CentGm-2, using a combination of computational and molecular cytogenetic approaches. These two subfamilies of satellite repeats mark distinct subsets of soybean centromeres and, in at least one case, a pair of homologs, suggesting their origins from an allopolyploid event. The satellite monomers of each subfamily are arranged in large tandem arrays, and intermingled monomers of the two subfamilies were not detected by fluorescence in situ hybridization on extended DNA fibers nor at the sequence level. This indicates that there has been little recombination and homogenization of satellite DNA between these two sets of centromeres. These satellite repeats are also present in Glycine soja, the proposed wild progenitor of soybean, but could not be detected in any other relatives of soybean examined in this study, suggesting the rapid divergence of the centromeric satellite DNA within the Glycine genus. Together, these observations provide direct evidence, at molecular and chromosomal levels, in support of the hypothesis that the soybean genome has experienced a recent allopolyploidization event.

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.

Evolutionary dynamics and sites of illegitimate recombination revealed in the interspersion and sequence junctions of two nonhomologous satellite DNAs in cactophilic Drosophila species

Satellite DNA (satDNA) is a major component of genomes but relatively little is known about the fine-scale organization of unrelated satDNAs residing at the same chromosome location, and the sequence structure and dynamics of satDNA junctions. We studied the organization and sequence junctions of two nonhomologous satDNAs, pBuM and DBC-150, in three species from the neotropical Drosophila buzzatii cluster (repleta group). In situ hybridization to microchromosomes, interphase nuclei and extended DNA fibers showed frequent interspersion of the two satellites in D. gouveai, D. antonietae and, to a lesser extent, D. seriema. We isolated by PCR six pBuM x DBC-150 junctions: four are exclusive to D. gouveai and two are exclusive to D. antonietae. The six junction breakpoints occur at different positions within monomers, suggesting independent origin. Four junctions showed abrupt transitions between the two satellites, whereas two junctions showed a distinct 10 bp tandem duplication before the junction. Unlike pBuM, DBC-150 junction repeats are more variable than randomly cloned monomers and showed diagnostic features in common to a 3-monomer higher-order repeat seen in the sister species D. serido. The high levels of interspersion between pBuM and DBC-150 repeats suggest extensive rearrangements between the two satellites, maybe favored by specific features of the microchromosomes. Our interpretation is that the junctions evolved by multiples events of illegitimate recombination between nonhomologous satDNA repeats, with subsequent rounds of unequal crossing-over expanding the copy number of some of the junctions.

Satellite DNA in insects: a review

The study of insect satellite DNAs (satDNAs) indicates the evolutionary conservation of certain features despite their sequence heterogeneity. Such features can include total length, monomer length, motifs, particular regions and/or secondary and tertiary structures. satDNAs may act as protein-binding sites, structural domains or sites for epigenetic modifications. The selective constraints in the evolution of satDNAs may be due to the satDNA sequence interaction with specific proteins important in heterochromatin formation and possible a role in controlling gene expression. The transcription of satDNA has been described in vertebrates, invertebrates and plants. In insects, differential satDNA expression has been observed in different cells, developmental stages, sex and caste of the individuals. These transcription differences may suggest their involvement in gene-regulation processes. In addition, the satDNA or its transcripts appear to be involved in heterochromatin formation and in chromatin-elimination processes. The importance of transposable elements to insect satDNA is shown by their presence as a constituent of satDNA in several species of insects (including possible active elements). In addition, they may be involved in the formation of centromeres and telomeres and in the homogenization and expansion of satDNA.

Sequence analysis, chromosomal distribution and long-range organization show that rapid turnover of new and old pBuM satellite DNA repeats leads to different patterns of variation in seven species of the Drosophila buzzatii cluster

We aimed to study patterns of variation and factors influencing the evolutionary dynamics of a satellite DNA, pBuM, in all seven Drosophila species from the buzzatii cluster (repleta group). We analyzed 117 alpha pBuM-1 (monomer length 190 bp) and 119 composite alpha/beta (370 bp) pBuM-2 repeats and determined the chromosome location and long-range organization on DNA fibers of major sequence variants. Such combined methodologies in the study of satDNAs have been used in very few organisms. In most species, concerted evolution is linked to high copy number of pBuM repeats. Species presenting low-abundance and scattered distributed pBuM repeats did not undergo concerted evolution and maintained part of the ancestral inter-repeat variability. The alpha and alpha/beta repeats colocalized in heterochromatic regions and were distributed on multiple chromosomes, with notable differences between species. High-resolution FISH revealed array sizes of a few kilobases to over 0.7 Mb and mutual arrangements of alpha and alpha/beta repeats along the same DNA fibers, but with considerable changes in the amount of each variant across species. From sequence, chromosomal and phylogenetic data, we could infer that homogenization and amplification events involved both new and ancestral pBuM variants. Altogether, the data on the structure and organization of the pBuM satDNA give insights into genome evolution including mechanisms that contribute to concerted evolution and diversification.

Meiotic transmission of an in vitro-assembled autonomous maize minichromosome

Autonomous chromosomes are generated in yeast (yeast artificial chromosomes) and human fibrosarcoma cells (human artificial chromosomes) by introducing purified DNA fragments that nucleate a kinetochore, replicate, and segregate to daughter cells. These autonomous minichromosomes are convenient for manipulating and delivering DNA segments containing multiple genes. In contrast, commercial production of transgenic crops relies on methods that integrate one or a few genes into host chromosomes; extensive screening to identify insertions with the desired expression level, copy number, structure, and genomic location; and long breeding programs to produce varieties that carry multiple transgenes. As a step toward improving transgenic crop production, we report the development of autonomous maize minichromosomes (MMCs). We constructed circular MMCs by combining DsRed and nptII marker genes with 7-190 kb of genomic maize DNA fragments containing satellites, retroelements, and/or other repeats commonly found in centromeres and using particle bombardment to deliver these constructs into embryogenic maize tissue. We selected transformed cells, regenerated plants, and propagated their progeny for multiple generations in the absence of selection. Fluorescent in situ hybridization and segregation analysis demonstrated that autonomous MMCs can be mitotically and meiotically maintained. The MMC described here showed meiotic segregation ratios approaching Mendelian inheritance: 93% transmission as a disome (100% expected), 39% transmission as a monosome crossed to wild type (50% expected), and 59% transmission in self crosses (75% expected). The fluorescent DsRed reporter gene on the MMC was expressed through four generations, and Southern blot analysis indicated the encoded genes were intact. This novel approach for plant transformation can facilitate crop biotechnology by (i) combining several trait genes on a single DNA fragment, (ii) arranging genes in a defined sequence context for more consistent gene expression, and (iii) providing an independent linkage group that can be rapidly introgressed into various germplasms.

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.

Meiotic transmission of an in vitro-assembled autonomous maize minichromosome

Autonomous chromosomes are generated in yeast (yeast artificial chromosomes) and human fibrosarcoma cells (human artificial chromosomes) by introducing purified DNA fragments that nucleate a kinetochore, replicate, and segregate to daughter cells. These autonomous minichromosomes are convenient for manipulating and delivering DNA segments containing multiple genes. In contrast, commercial production of transgenic crops relies on methods that integrate one or a few genes into host chromosomes; extensive screening to identify insertions with the desired expression level, copy number, structure, and genomic location; and long breeding programs to produce varieties that carry multiple transgenes. As a step toward improving transgenic crop production, we report the development of autonomous maize minichromosomes (MMCs). We constructed circular MMCs by combining DsRed and nptII marker genes with 7-190 kb of genomic maize DNA fragments containing satellites, retroelements, and/or other repeats commonly found in centromeres and using particle bombardment to deliver these constructs into embryogenic maize tissue. We selected transformed cells, regenerated plants, and propagated their progeny for multiple generations in the absence of selection. Fluorescent in situ hybridization and segregation analysis demonstrated that autonomous MMCs can be mitotically and meiotically maintained. The MMC described here showed meiotic segregation ratios approaching Mendelian inheritance: 93% transmission as a disome (100% expected), 39% transmission as a monosome crossed to wild type (50% expected), and 59% transmission in self crosses (75% expected). The fluorescent DsRed reporter gene on the MMC was expressed through four generations, and Southern blot analysis indicated the encoded genes were intact. This novel approach for plant transformation can facilitate crop biotechnology by (i) combining several trait genes on a single DNA fragment, (ii) arranging genes in a defined sequence context for more consistent gene expression, and (iii) providing an independent linkage group that can be rapidly introgressed into various germplasms.

Global repeat discovery and estimation of genomic copy number in a large, complex genome using a high-throughput 454 sequence survey

Background

Extensive computational and database tools are available to mine genomic and genetic databases for model organisms, but little genomic data is available for many species of ecological or agricultural significance, especially those with large genomes. Genome surveys using conventional sequencing techniques are powerful, particularly for detecting sequences present in many copies per genome. However these methods are time-consuming and have potential drawbacks. High throughput 454 sequencing provides an alternative method by which much information can be gained quickly and cheaply from high-coverage surveys of genomic DNA.

Results

We sequenced 78 million base-pairs of randomly sheared soybean DNA which passed our quality criteria. Computational analysis of the survey sequences provided global information on the abundant repetitive sequences in soybean. The sequence was used to determine the copy number across regions of large genomic clones or contigs and discover higher-order structures within satellite repeats. We have created an annotated, online database of sequences present in multiple copies in the soybean genome. The low bias of pyrosequencing against repeat sequences is demonstrated by the overall composition of the survey data, which matches well with past estimates of repetitive DNA content obtained by DNA re-association kinetics (Cot analysis).

Conclusion

This approach provides a potential aid to conventional or shotgun genome assembly, by allowing rapid assessment of copy number in any clone or clone-end sequence. In addition, we show that partial sequencing can provide access to partial protein-coding sequences.

Global repeat discovery and estimation of genomic copy number in a large, complex genome using a high-throughput 454 sequence survey

BACKGROUND

Extensive computational and database tools are available to mine genomic and genetic databases for model organisms, but little genomic data is available for many species of ecological or agricultural significance, especially those with large genomes. Genome surveys using conventional sequencing techniques are powerful, particularly for detecting sequences present in many copies per genome. However these methods are time-consuming and have potential drawbacks. High throughput 454 sequencing provides an alternative method by which much information can be gained quickly and cheaply from high-coverage surveys of genomic DNA.

RESULTS

We sequenced 78 million base-pairs of randomly sheared soybean DNA which passed our quality criteria. Computational analysis of the survey sequences provided global information on the abundant repetitive sequences in soybean. The sequence was used to determine the copy number across regions of large genomic clones or contigs and discover higher-order structures within satellite repeats. We have created an annotated, online database of sequences present in multiple copies in the soybean genome. The low bias of pyrosequencing against repeat sequences is demonstrated by the overall composition of the survey data, which matches well with past estimates of repetitive DNA content obtained by DNA re-association kinetics (Cot analysis).

CONCLUSION

This approach provides a potential aid to conventional or shotgun genome assembly, by allowing rapid assessment of copy number in any clone or clone-end sequence. In addition, we show that partial sequencing can provide access to partial protein-coding sequences.

Chromosome-specific satellite sequences in Turritis glabra

Two novel repetitive sequence families were isolated from Turritis glabra (2n = 2x = 12). These two repeat families are similar to those of centromeric repeats in Arabidopsis thaliana, are co-localized on one chromosome pair, and differ by about 20% from each other. Phylogenetic analysis revealed that the two repeat families of T. glabra are more similar to each other than to the centromeric repeat families of other Arabidopsis and related species. The relationships of satellite sequences reflected the species phylogeny, indicating that the replacement of satellite sequences has occurred in each species lineage independently, and shared variants could not have existed for a long time between species.

Satellite DNA junctions identify the potential origin of new repetitive elements in the beetle Tribolium madens

Two related satellite DNA families (satellite I and satellite II) with complex higher-order repeat (HOR) monomers represent major DNA components equilocated in the pericentromeric heterochromatin of all Tribolium madens chromosomes. Fragments obtained upon genomic DNA restriction revealed two subfamilies of satellite II monomers, and also identified regions of transition between satellite I and satellite II sequences. The two subfamilies differ not only in diagnostic nucleotides, but also in flipped orientation of constituent subunits. Hybrid genomic fragments comprise directly linked satellite I and satellite II monomers that cannot be distinguished from randomly cloned monomers of corresponding families. An exception is the most proximal satellite I monomer in the hybrid fragment named TMADhinf, which shows sequence divergence typical for repeats evolving at array ends, in zones of low homogenization efficiency. This pattern points to the extensive rearrangement processes generating abrupt transitions between satellite arrays combined with array maintenance by unequal crossover. Switching points between adjacent satellites as well as the edges of flipped subunits are localized within a short sequence segment, indicating a preferential site of recombination within satellite subunits. Multiple copies of TMADhinf junction fragment support the hypothesis that sites of evolutionary origin of novel satellite repeat (sub)families can be localized at array ends, in regions of enhanced sequence divergence.

Sequence of events leading to near-complete genome turnover in allopolyploid Nicotiana within five million years

Analyses of selected bacterial artificial chromosomes (BACs) clones suggest that the retrotransposon component of angiosperm genomes can be amplified or deleted, leading to genome turnover. Here, Nicotiana allopolyploids were used to characterize the nature of sequence turnover across the whole genome in allopolyploids known to be of different ages. Using molecular-clock analyses, the likely age of Nicotiana allopolyploids was estimated. Genomic in situ hybridization (GISH) and tandem repeat characterization were used to determine how the parental genomic compartments of these allopolyploids have diverged over time. Paternal genome sequence losses, retroelement activity and intergenomic translocation have been reported in early Nicotiana tabacum evolution (up to 200,000 yr divergence). Here it is shown that within 1 million years of allopolyploid divergence there is considerable exchange of repeats between parental chromosome sets. After c. 5 million years of divergence GISH fails. This GISH failure may represent near-complete genome turnover, probably involving the replacement of nongenic sequences with new, or previously rare sequence types, all occurring within a conserved karyotype structure. This mode of evolution may influence or be influenced by long-term diploidization processes that characterize angiosperm polyploidy-diploid evolutionary cycles.

Patterns of DNA variation among three centromere satellite families in Arabidopsis halleri and A. lyrata

We describe patterns of DNA variation among the three centromeric satellite families in Arabidopsis halleri and lyrata. The newly studied subspecies (A. halleri ssp. halleri and A. lyrata ssp. lyrata and petraea), like the previously studied A. halleri ssp. gemmifera and A. lyrata ssp. kawasakiana, have three different centromeric satellite families, the older pAa family (also present in A. arenosa) and two families, pAge1 and pAge2, that probably evolved more recently. Sequence variability is high in all three satellite families, and the pAa sequences do not cluster by their species of origin. Diversity in the pAge2 family is complex, and different from variation among copies of the other two families, showing clear evidence for exchange events among family members, especially in A. halleri ssp. halleri. In A. lyrata ssp. lyrata there is some evidence for recent rapid spread of pAge2 variants, suggesting selection favoring these sequences.

Centromeres put epigenetics in the driver's seat

A defining feature of chromosomes is the centromere, the site for spindle attachment at mitosis and meiosis. Intriguingly, centromeres of plants and animals are maintained by both sequence-specific and sequence-independent (epigenetic) processes. Epigenetic inheritance might enable kinetochores (the structures that attach centromeres to spindles) to maintain an optimal size. However, centromeres are susceptible to the evolution of "selfish" DNA repeats that bind to kinetochore proteins. We argue that such sequence-specific interactions are evolutionarily unstable because they enable repeat arrays to influence kinetochore size. Changes in kinetochore size could affect the interaction of kinetochores with the spindle and, in principle, skew Mendelian segregation. We propose that key kinetochore proteins have adapted to disrupt such sequence-specific interactions and restore epigenetic inheritance.

Chromosome arrangement and nuclear architecture but not centromeric sequences are conserved between Arabidopsis thaliana and Arabidopsis lyrata

In contrast to the situation described for mammals and Drosophila, chromosome territory (CT) arrangement and somatic homologous pairing in interphase nuclei of Arabidopsis thaliana (n = 5) are predominantly random except for a more frequent association of the chromosomes bearing a homologous nucleolus organizer region. To find out whether this chromosome arrangement is also characteristic for other species of the genus Arabidopsis, we investigated Arabidopsis lyrata ssp. lyrata (n = 8), one of the closest relatives of A. thaliana. First, we determined the size of each chromosome and chromosome arm, the sequence type of centromeric repeats and their distribution between individual centromeres and the position of the 5S/45S rDNA arrays in A. lyrata. Then we demonstrated that CT arrangement, homologous pairing and sister chromatid alignment of distinct euchromatic and/or heterochromatic regions within A. lyrata interphase nuclei are similar to that in A. thaliana nuclei. Thus, the arrangement of interphase chromosomes appears to be conserved between both taxa that diverged about 5 million years ago. Since the chromosomes of A. lyrata resemble those of the presumed ancestral karyotype, a similar arrangement of interphase chromosomes is also to be expected for other closely related diploid species of the Brassicaceae family.

Transcription and Evolutionary Dynamics of the Centromeric Satellite Repeat CentO in Rice

Satellite DNA is a major component of centromeric heterochromatin in most multicellular eukaryotes, where it is typically organized into megabase-sized tandem arrays. It has recently been demonstrated that small interfering RNAs (siRNAs) processed from centromeric satellite repeats can be involved in epigenetic chromatin modifications which appear to underpin centromere function. However, the structural organization and evolution of the centromeric satellite DNA is still poorly understood. We analyzed the centromeric satellite repeat arrays from rice chromosomes 1 and 8 and identified higher order structures and local homogenization of the CentO repeats in these 2 centromeres. We also cloned the CentO repeats from the CENH3-associated nucleosomes by a chromatin immunoprecipitation (ChIP)-based method. Sequence variability analysis of the ChIPed CentO repeats revealed a single variable domain within the repeat. We detected transcripts derived from both strands of the CentO repeats. The CentO transcripts are processed into siRNA, suggesting a potential role of this satellite repeat family in epigenetic chromatin modification.

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.

Sequence homogenization and chromosomal localization of VicTR-B satellites differ between closely related Vicia species

Satellite sequences of the VicTR-B family are specific for the genus Vicia (Leguminosae), but their abundance varies among the species, being the highest in Vicia sativa and Vicia grandiflora. In this study, we have sequenced multiple randomly cloned VicTR-B fragments from these two species and analyzed their sequence variability, periodicity, and chromosomal localization. We have found that V. sativa VicTR-B sequences are homogeneous with respect to their nucleotide sequences and periodicity (monomers of 38 bp), whereas V. grandiflora repeats are considerably more variable, occurring in at least four distinct sequence subfamilies. Although the periodicity of 38 bp was conserved in most of the V. grandiflora sequences, one of the subfamilies was composed of higher-order repeats of 186 bp, which originated from a pentamer of the basic repeated unit. Individual VicTR-B subfamilies were preferentially located in either intercalary or subtelomeric regions of chromosomes. Interestingly, two V. grandiflora subfamilies with the highest similarity to V. sativa VicTR-B sequences were located in intercalary heterochromatic bands, showing similar chromosomal distribution as the majority of VicTR-B repeats in V. sativa. The other two V. grandiflora subfamilies showing a considerable divergence from V. sativa sequences were found to be accumulated at subtelomeric regions of V. grandiflora chromosomes.

The origin, meiotic behavior, and transmission of a novel minichromosome in Arabidopsis thaliana

A plant carrying a small extra chromosome was found in Landsberg erecta ecotype of Arabidopsis thaliana. Fluorescence in situ hybridization revealed that this minichromosome was derived from the short arm of chromosome 4. The size of this "mini4S" chromosome was estimated to be approximately 7.5 Mb on the basis of previously reported data and the amount of the centromeric major satellite (180-bp family) present, which was determined to be about 1 Mb, or about one third of that in the normal chromosome 4. No pairing between mini4S and its original chromosome 4 was observed at pachytene and metaphase I stages. The transmission of mini4S through pollen was limited, but about 30% of selfed progeny carried the mini4S chromosomes. The transmission rates considerably increased when the mini4S chromosomes were transferred to plants with a Columbia background by successive backcrosses. This suggests that the stability of the minichromosomes is controlled genetically by factors that can vary between ecotypes.

Dynamic evolution at pericentromeres

Pericentromeres are exceptional genomic regions: in animals they contain extensive segmental duplications implicated in gene creation, and in plants they sustain rearrangements and insertions uncommon in euchromatin. To examine the mechanisms and patterns of plant pericentromere evolution, we compared pericentromere sequence from four Brassicaceae species separated by <15 million years (Myr). This flowering plant family is ideal for studying relationships between genome reorganization and pericentromere evolution-its members have undergone recent polyploidization and hybridization, with close relatives changing in genome size and chromosome number. Through sequence and hybridization analyses, we examined regions from Arabidopsis arenosa, Capsella rubella, and Olimarabidopsis pumila that are homologous to Arabidopsis thaliana pericentromeres (peri-CENs) III and V, and used FISH to demonstrate they have been maintained near centromere satellite arrays in each species. Sequence analysis revealed a set of highly conserved genes, yet we discovered substantial differences in intergenic length and species-specific changes in sequence content and gene density. We discovered that A. thaliana has undergone recent, significant expansions within its pericentromeres, in some cases measuring hundreds of kilobases; these findings are in marked contrast to euchromatic segments in these species that exhibit only minor length changes. While plant pericentromeres do contain some duplications, we did not find evidence of extensive segmental duplications, as has been documented in primates. Our data support a model in which plant pericentromeres may experience selective pressures distinct from euchromatin, tolerating rapid, dynamic changes in structure and sequence content, including large insertions of mobile elements, 5S rDNA arrays and pseudogenes.