Possible role of natural selection in the formation of tandem-repetitive noncoding DNA

The Structure of Simple Satellite Variation in the Human Genome and Its Correlation With Centromere Ancestry

Although repetitive DNA forms much of the human genome, its study is challenging due to limitations in assembly and alignment of repetitive short-reads. We have deployed k-Seek, software that detects tandem repeats embedded in single reads, on 2,504 human genomes from the 1,000 Genomes Project to quantify the variation and abundance of simple satellites (repeat units <20 bp). We find that the ancestral monomer of Human Satellite 3 makes up the largest portion of simple satellite content in humans (mean of ∼8 Mb). We discovered ∼50,000 rare tandem repeats that are not detected in the T2T-CHM13v2.0 assembly, including undescribed variants of telomericand pericentromeric repeats. We find broad homogeneity of the most abundant repeats across populations, except for AG-rich repeats which are more abundant in African individuals. We also find cliques of highly similar AG- and AT-rich satellites that are interspersed and form higher-order structures that covary in copy number across individuals, likely through concerted amplification via unequal exchange. Finally, we use pericentromeric polymorphisms to estimate centromeric genetic relatedness between individuals and find a strong predictive relationship between centromeric lineages and pericentromeric simple satellite abundances. In particular, ancestral monomers of Human Satellite 2 and Human Satellite 3 abundances correlate with clusters of centromeric ancestry on chromosome 16 and chromosome 9, with some clusters structured by population. These results provide new descriptions of the population dynamics that underlie the evolution of simple satellites in humans.

The implications of satellite DNA instability on cellular function and evolution

Abundant tandemly repeated satellite DNA is present in most eukaryotic genomes. Previous limitations including a pervasive view that it was uninteresting junk DNA, combined with challenges in studying it, are starting to dissolve - and recent studies have found important functions for satellite DNAs. The observed rapid evolution and implied instability of satellite DNA now has important significance for their functions and maintenance within the genome. In this review, we discuss the processes that lead to satellite DNA copy number instability, and the importance of mechanisms to manage the potential negative effects of instability. Satellite DNA is vulnerable to challenges during replication and repair, since it forms difficult-to-process secondary structures and its homology within tandem arrays can result in various types of recombination. Satellite DNA instability may be managed by DNA or chromatin-binding proteins ensuring proper nuclear localization and repair, or by proteins that process aberrant structures that satellite DNAs tend to form. We also discuss the pattern of satellite DNA mutations from recent mutation accumulation (MA) studies that have tracked changes in satellite DNA for up to 1000 generations with minimal selection. Finally, we highlight examples of satellite evolution from studies that have characterized satellites across millions of years of Drosophila fruit fly evolution, and discuss possible ways that selection might act on the satellite DNA composition.

Conservation Genetics of the Asian Giant Soft-Shelled Turtle (Pelochelys cantorii) with Novel Microsatellite Multiplexes

To understand the genetic structure of the protected turtle species Pelochelys cantorii we used transcriptome data to design more than 30,000 tri- and tetranucleotide repeat microsatellite primer pairs, of which 230 pairs were used for laboratory experiments. After two screenings, only 10 microsatellite markers with good specificity, high amplification efficiency, and polymorphisms were obtained. Using the selected primers, two multiplex PCR systems were established to compare and analyze the genetic diversity of artificially assisted breeding generations from four parents (two females and two males) continuously bred over two years. A total of 25 alleles were detected among the 10 microsatellite loci of the offspring. The polymorphism information content (PIC) was 0.313-0.674, with an average of 0.401, among which two loci were highly polymorphic (PIC ≥ 0.5). The number of alleles was 2-5 and the number of effective alleles was 1.635-3.614. The observed heterozygosity was 0.341-0.813, with an average of 0.582, whereas the average expected heterozygosity was 0.389-0.725, with an average of 0.493. Moreover, on the basis of Nei's genetic distance and the Bayesian clustering algorithm, the 182 offspring were divided into two subgroups, and the results corresponded to the two maternal lines. This is the first study to investigate the molecular markers of P. cantorii, and the results obtained can be used to protect genetic resources and provide a genetic basis for the design of population recovery plans.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Deciphering the Mechanism of Glyphosate Resistance in Amaranthus palmeri by Cytogenomics

In agriculture, various chemicals are used to control the weeds. Out of which, glyphosate is an important herbicide invariably used in the cultivation of glyphosate-resistant crops to control weeds. Overuse of glyphosate results in the evolution of glyphosate-resistant weeds. Evolution of glyphosate resistance (GR) in Amaranthus palmeri (AP) is a serious concern in the USA. Investigation of the mechanism of GR in AP identified different resistance mechanisms of which 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS) gene amplification is predominant. Molecular analysis of GR AP identified the presence of a 5- to >160-fold increase in copies of the EPSPS gene than in a glyphosate-susceptible (GS) population. This increased copy number of the EPSPS gene increased the genome size ranging from 3.5 to 11.8%, depending on the copy number compared to the genome size of GS AP. FISH analysis using a 399-kb EPSPS cassette derived from bacterial artificial chromosomes (BACs) as probes identified that amplified EPSPS copies in GR AP exist in extrachromosomal circular DNA (eccDNA) in addition to the native copy in the chromosome. The EPSPS gene-containing eccDNA having a size of ∼400 kb is termed EPSPS-eccDNA and showed somatic mosacism in size and copy number. EPSPS-eccDNA has a genetic mechanism to tether randomly to mitotic or meiotic chromosomes during cell division or gamete formation and is inherited to daughter cells or progeny generating copy number variation. These eccDNAs are stable genetic elements that can replicate and exist independently. The genomic characterization of the EPSPS locus, along with the flanking regions, identified the presence of a complex array of repeats and mobile genetic elements. The cytogenomics approach in understanding the biology of EPSPS-eccDNA sheds light on various characteristics of EPSPS-eccDNA that favor GR in AP.

Variation and Evolution of Human Centromeres: A Field Guide and Perspective

We are entering a new era in genomics where entire centromeric regions are accurately represented in human reference assemblies. Access to these high-resolution maps will enable new surveys of sequence and epigenetic variation in the population and offer new insight into satellite array genomics and centromere function. Here, we focus on the sequence organization and evolution of alpha satellites, which are credited as the genetic and genomic definition of human centromeres due to their interaction with inner kinetochore proteins and their importance in the development of human artificial chromosome assays. We provide an overview of alpha satellite repeat structure and array organization in the context of these high-quality reference data sets; discuss the emergence of variation-based surveys; and provide perspective on the role of this new source of genetic and epigenetic variation in the context of chromosome biology, genome instability, and human disease.

Copy number evolution in simple and complex tandem repeats across the C57BL/6 and C57BL/10 inbred mouse lines

Simple sequence tandem repeats are among the most rapidly evolving compartments of the genome. Some repeat expansions are associated with mammalian disease or meiotic segregation distortion, yet the rates of copy number change across generations are not well known. Here, we use 14 distinct sublineages of the C57BL/6 and C57BL/10 inbred mouse strains, which have been evolving independently over about 300 generations, to estimate the rates of copy number changes in genome-wide tandem repeats. Rates of change varied across repeats and across lines. Notably, CAG, whose expansions in coding regions are associated with many neurological and genetic disorders, was highly stable in copy number, likely indicating stabilizing selection. Rates of change were positively correlated with copy number, but the direction and magnitude of changes varied across lines. Some mouse lines experienced consistent losses or gains across most simple repeats, but this did not correlate with copy number changes in complex repeats. Rates of copy number change were similar between simple repeats and the more abundant complex repeats after normalization by copy number. Finally, the Y-specific centromeric repeat had a fourfold higher rate of change than the homologous centromeric repeat on other chromosomes. Structural differences in satellite complexity, or restriction to the Y chromosome and elevated mutation rates of the male germline, may explain the higher rate of change. Overall, our work underscores the mutational fluidity of long tandem arrays of repeats, and the correlations and constraints between genome-wide tandem repeats, which suggest that turnover is not a completely neutral process.

Sequence, Chromatin and Evolution of Satellite DNA

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

Functional Significance of Satellite DNAs: Insights From Drosophila

Since their discovery more than 60 years ago, satellite repeats are still one of the most enigmatic parts of eukaryotic genomes. Being non-coding DNA, satellites were earlier considered to be non-functional “junk,” but recently this concept has been extensively revised. Satellite DNA contributes to the essential processes of formation of crucial chromosome structures, heterochromatin establishment, dosage compensation, reproductive isolation, genome stability and development. Genomic abundance of satellites is under stabilizing selection owing of their role in the maintenance of vital regions of the genome – centromeres, pericentromeric regions, and telomeres. Many satellites are transcribed with the generation of long or small non-coding RNAs. Misregulation of their expression is found to lead to various defects in the maintenance of genomic architecture, chromosome segregation and gametogenesis. This review summarizes our current knowledge concerning satellite functions, the mechanisms of regulation and evolution of satellites, focusing on recent findings in Drosophila. We discuss here experimental and bioinformatics data obtained in Drosophila in recent years, suggesting relevance of our analysis to a wide range of eukaryotic organisms.

Characterization of repeat arrays in ultra-long nanopore reads reveals frequent origin of satellite DNA from retrotransposon-derived tandem repeats

Amplification of monomer sequences into long contiguous arrays is the main feature distinguishing satellite DNA from other tandem repeats, yet it is also the main obstacle in its investigation because these arrays are in principle difficult to assemble. Here we explore an alternative, assembly-free approach that utilizes ultra-long Oxford Nanopore reads to infer the length distribution of satellite repeat arrays, their association with other repeats and the prevailing sequence periodicities. Using the satellite DNA-rich legume plant Lathyrus sativus as a model, we demonstrated this approach by analyzing 11 major satellite repeats using a set of nanopore reads ranging from 30 to over 200 kb in length and representing 0.73× genome coverage. We found surprising differences between the analyzed repeats because only two of them were predominantly organized in long arrays typical for satellite DNA. The remaining nine satellites were found to be derived from short tandem arrays located within LTR-retrotransposons that occasionally expanded in length. While the corresponding LTR-retrotransposons were dispersed across the genome, this array expansion occurred mainly in the primary constrictions of the L. sativus chromosomes, which suggests that these genome regions are favourable for satellite DNA accumulation.

Dynamic turnover of centromeres drives karyotype evolution in Drosophila

Centromeres are the basic unit for chromosome inheritance, but their evolutionary dynamics is poorly understood. We generate high-quality reference genomes for multiple Drosophila obscura group species to reconstruct karyotype evolution. All chromosomes in this lineage were ancestrally telocentric and the creation of metacentric chromosomes in some species was driven by de novo seeding of new centromeres at ancestrally gene-rich regions, independently of chromosomal rearrangements. The emergence of centromeres resulted in a drastic size increase due to repeat accumulation, and dozens of genes previously located in euchromatin are now embedded in pericentromeric heterochromatin. Metacentric chromosomes secondarily became telocentric in the pseudoobscura subgroup through centromere repositioning and a pericentric inversion. The former (peri)centric sequences left behind shrunk dramatically in size after their inactivation, yet contain remnants of their evolutionary past, including increased repeat-content and heterochromatic environment. Centromere movements are accompanied by rapid turnover of the major satellite DNA detected in (peri)centromeric regions.

Function of Junk: Pericentromeric Satellite DNA in Chromosome Maintenance

Satellite DNAs are simple tandem repeats that exist at centromeric and pericentromeric regions on eukaryotic chromosomes. Unlike the centromeric satellite DNA that comprises the vast majority of natural centromeres, function(s) for the much more abundant pericentromeric satellite repeats are poorly understood. In fact, the lack of coding potential allied with rapid divergence of repeat sequences across eukaryotes has led to their dismissal as "junk DNA" or "selfish parasites." Although implicated in various biological processes, a conserved function for pericentromeric satellite DNA remains unidentified. We have addressed the role of satellite DNA through studying chromocenters, a cytological aggregation of pericentromeric satellite DNA from multiple chromosomes into DNA-dense nuclear foci. We have shown that multivalent satellite DNA-binding proteins cross-link pericentromeric satellite DNA on chromosomes into chromocenters. Disruption of chromocenters results in the formation of micronuclei, which arise by budding off the nucleus during interphase. We propose a model that satellite DNAs are critical chromosome elements that are recognized by satellite DNA-binding proteins and incorporated into chromocenters. We suggest that chromocenters function to preserve the entire chromosomal complement in a single nucleus, a fundamental and unquestioned feature of eukaryotic genomes. We speculate that the rapid divergence of satellite DNA sequences between closely related species results in discordant chromocenter function and may underlie speciation and hybrid incompatibility.

Double insertion of transposable elements provides a substrate for the evolution of satellite DNA

Eukaryotic genomes are replete with repeated sequences in the form of transposable elements (TEs) dispersed across the genome or as satellite arrays, large stretches of tandemly repeated sequences. Many satellites clearly originated as TEs, but it is unclear how mobile genetic parasites can transform into megabase-sized tandem arrays. Comprehensive population genomic sampling is needed to determine the frequency and generative mechanisms of tandem TEs, at all stages from their initial formation to their subsequent expansion and maintenance as satellites. The best available population resources, short-read DNA sequences, are often considered to be of limited utility for analyzing repetitive DNA due to the challenge of mapping individual repeats to unique genomic locations. Here we develop a new pipeline called ConTExt that demonstrates that paired-end Illumina data can be successfully leveraged to identify a wide range of structural variation within repetitive sequence, including tandem elements. By analyzing 85 genomes from five populations of Drosophila melanogaster, we discover that TEs commonly form tandem dimers. Our results further suggest that insertion site preference is the major mechanism by which dimers arise and that, consequently, dimers form rapidly during periods of active transposition. This abundance of TE dimers has the potential to provide source material for future expansion into satellite arrays, and we discover one such copy number expansion of the DNA transposon hobo to approximately 16 tandem copies in a single line. The very process that defines TEs—transposition—thus regularly generates sequences from which new satellites can arise.

High-throughput analysis of satellite DNA in the grasshopper Pyrgomorpha conica reveals abundance of homologous and heterologous higher-order repeats

Satellite DNA (satDNA) constitutes an important fraction of repetitive DNA in eukaryotic genomes, but it is barely known in most species. The high-throughput analysis of satDNA in the grasshopper Pyrgomorpha conica revealed 87 satDNA variants grouped into 76 different families, representing 9.4% of the genome. Fluorescent in situ hybridization (FISH) analysis of the 38 most abundant satDNA families revealed four different patterns of chromosome distribution. Homology search between the 76 satDNA families showed the existence of 15 superfamilies, each including two or more families, with the most abundant superfamily representing more than 80% of all satDNA found in this species. This also revealed the presence of two types of higher-order repeats (HORs), one showing internal homologous subrepeats, as conventional HORs, and an additional type showing non-homologous internal subrepeats, the latter arising by the combination of a given satDNA family with a non-annotated sequence, or with telomeric DNA. Interestingly, the heterologous subrepeats included in these HORs showed higher divergence within the HOR than outside it, suggesting that heterologous HORs show poor homogenization, in high contrast with conventional (homologous) HORs. Finally, heterologous HORs can show high differences in divergence between their constituent subrepeats, suggesting the possibility of regional homogenization.

The effects of transcription and recombination on mutational dynamics of short tandem repeats

Short tandem repeats (STR) are ubiquitous components of the genomic architecture of most living organisms. Recent work has highlighted the widespread functional significance of such repeats, particularly around gene regulation, but the mutational processes underlying the evolution of these highly abundant and highly variable sequences are not fully understood. Traditional models assume that strand misalignment during replication is the predominant mechanism, but empirical data suggest the involvement of other processes including recombination and transcription. Despite this evidence, the relative influences of these processes have not previously been tested experimentally on a genome-wide scale. Using deep sequencing, we identify mutations at >200 microsatellites, across 700 generations in replicated populations of two otherwise identical sexual and asexual Saccharomyces cerevisiae strains. Using generalized linear models, we investigate correlates of STR mutability including the nature of the mutation, STR composition and contextual factors including recombination, transcription and replication origins. Sexual capability was not a significant predictor of microsatellite mutability, but, intriguingly, we identify transcription as a significant positive predictor. We also find that STR density is substantially increased in regions neighboring, but not within, recombination hotspots.

Lack of satellite DNA species-specific homogenization and relationship to chromosomal rearrangements in monitor lizards (Varanidae, Squamata)

BACKGROUND

Satellite DNAs (stDNAs) are highly repeated sequences that constitute large portions of any genome. The evolutionary dynamics of stDNA (e.g. copy number, nucleotide sequence, location) can, therefore, provide an insight into genome organization and evolution. We investigated the evolutionary origin of VSAREP stDNA in 17 monitor lizards (seven Asian, five Australian, and five African) at molecular and cytogenetic level.

RESULTS

Results revealed that VSAREP is conserved in the genome of Asian and Australian varanids, but not in African varanids, suggesting that these sequences are either differentiated or lost in the African varanids. Phylogenetic and arrangement network analyses revealed the existence of at least four VSAREP subfamilies. The similarity of each sequence unit within the same VSAREP subfamily from different species was higher than those of other VSAREP subfamilies belonging to the same species. Additionally, all VSAREP subfamilies isolated from the three Australian species (Varanus rosenbergi, V. gouldii, and V. acanthurus) were co-localized near the centromeric or pericentromeric regions of the macrochromosomes, except for chromosomes 3 and 4 in each Australian varanid. However, their chromosomal arrangements were different among species.

CONCLUSIONS

The VSAREP stDNA family lack homogenized species-specific nucleotide positions in varanid lineage. Most VSAREP sequences were shared among varanids within the four VSAREP subfamilies. This suggests that nucleotide substitutions in each varanid species accumulated more slowly than homogenization rates in each VSAREP subfamily, resulting in non-species-specific evolution of stDNA profiles. Moreover, changes in location of VSAREP stDNA in each Australian varanid suggests a correlation with chromosomal rearrangements, leading to karyotypic differences among these species.

Selection Constrains High Rates of Tandem Repetitive DNA Mutation in Daphnia pulex

A long-standing evolutionary puzzle is that all eukaryotic genomes contain large amounts of tandemly-repeated DNA whose sequence motifs and abundance vary greatly among even closely related species. To elucidate the evolutionary forces governing tandem repeat dynamics, quantification of the rates and patterns of mutations in repeat copy number and tests of its selective neutrality are necessary. Here, we used whole-genome sequences of 28 mutation accumulation (MA) lines of Daphnia pulex, in addition to six isolates from a non-MA population originating from the same progenitor, to both estimate mutation rates of abundances of repeat sequences and evaluate the selective regime acting upon them. We found that mutation rates of individual repeats were both high and highly variable, ranging from additions/deletions of 0.29-105 copies per generation (reflecting changes of 0.12-0.80% per generation). Our results also provide evidence that new repeat sequences are often formed from existing ones. The non-MA population isolates showed a signal of either purifying or stabilizing selection, with 33% lower variation in repeat copy number on average than the MA lines, although the level of selective constraint was not evenly distributed across all repeats. The changes between many pairs of repeats were correlated, and the pattern of correlations was significantly different between the MA lines and the non-MA population. Our study demonstrates that tandem repeats can experience extremely rapid evolution in copy number, which can lead to high levels of divergence in genome-wide repeat composition between closely related species.

Dissecting the Satellite DNA Landscape in Three Cactophilic Drosophila Sequenced Genomes

Eukaryote genomes are replete with repetitive DNAs. This class includes tandemly repeated satellite DNAs (satDNA) which are among the most abundant, fast evolving (yet poorly studied) genomic components. Here, we used high-throughput sequencing data from three cactophilic Drosophila species, D. buzzatii, D. seriema, and D. mojavensis, to access and study their whole satDNA landscape. In total, the RepeatExplorer software identified five satDNAs, three previously described (pBuM, DBC-150 and CDSTR198) and two novel ones (CDSTR138 and CDSTR130). Only pBuM is shared among all three species. The satDNA repeat length falls within only two classes, between 130 and 200 bp or between 340 and 390 bp. FISH on metaphase and polytene chromosomes revealed the presence of satDNA arrays in at least one of the following genomic compartments: centromeric, telomeric, subtelomeric, or dispersed along euchromatin. The chromosomal distribution ranges from a single chromosome to almost all chromosomes of the complement. Fiber-FISH and sequence analysis of contigs revealed interspersion between pBuM and CDSTR130 in the microchromosomes of D. mojavensis Phylogenetic analyses showed that the pBuM satDNA underwent concerted evolution at both interspecific and intraspecific levels. Based on RNA-seq data, we found transcription activity for pBuM (in D. mojavensis) and CDSTR198 (in D. buzzatii) in all five analyzed developmental stages, most notably in pupae and adult males. Our data revealed that cactophilic Drosophila present the lowest amount of satDNAs (1.9-2.9%) within the Drosophila genus reported so far. We discuss how our findings on the satDNA location, abundance, organization, and transcription activity may be related to functional aspects.

High-throughput analysis of the satellitome revealed enormous diversity of satellite DNAs in the neo-Y chromosome of the cricket Eneoptera surinamensis

Satellite DNAs (satDNAs) constitute large portion of eukaryote genomes, comprising non-protein-coding sequences tandemly repeated. They are mostly found in heterochromatic regions of chromosomes such as around centromere or near telomeres, in intercalary heterochromatin, and often in non-recombining segments of sex chromosomes. We examined the satellitome in the cricket Eneoptera surinamensis (2n = 9, neo-X1X2Y, males) to characterize the molecular evolution of its neo-sex chromosomes. To achieve this, we analyzed illumina reads using graph-based clustering and complementary analyses. We found an unusually high number of 45 families of satDNAs, ranging from 4 bp to 517 bp, accounting for about 14% of the genome and showing different modular structures and high diversity of arrays. FISH mapping revealed that satDNAs are located mostly in C-positive pericentromeric regions of the chromosomes. SatDNAs enrichment was also observed in the neo-sex chromosomes in comparison to autosomes. Especially astonishing accumulation of satDNAs loci was found in the highly differentiated neo-Y, including 39 satDNAs over-represented in this chromosome, which is the greatest satDNAs diversity yet reported for sex chromosomes. Our results suggest possible involvement of satDNAs in genome increasing and in molecular differentiation of the neo-sex chromosomes in this species, contributing to the understanding of sex chromosome composition and evolution in Orthoptera.

Structurally Complex Organization of Repetitive DNAs in the Genome of Cobia (Rachycentron canadum)

Repetitive DNAs comprise the largest fraction of the eukaryotic genome. They include microsatellites or simple sequence repeats (SSRs), which play an important role in the chromosome differentiation among fishes. Rachycentron canadum is the only representative of the family Rachycentridae. This species has been focused on several multidisciplinary studies in view of its important potential for marine fish farming. In the present study, distinct classes of repetitive DNAs, with emphasis on SSRs, were mapped in the chromosomes of this species to improve the knowledge of its genome organization. Microsatellites exhibited a diversified distribution, both dispersed in euchromatin and clustered in the heterochromatin. The multilocus location of SSRs strengthened the heterochromatin heterogeneity in this species, as suggested by some previous studies. The colocalization of SSRs with retrotransposons and transposons pointed to a close evolutionary relationship between these repetitive sequences. A number of heterochromatic regions highlighted a greater complex organization than previously supposed, harboring a diversity of repetitive elements. In this sense, there was also evidence of colocalization of active genetic regions and different classes of repetitive DNAs in a common heterochromatic region, which offers a potential opportunity for further researches regarding the interaction of these distinct fractions in fish genomes.

Correlated variation and population differentiation in satellite DNA abundance among lines of Drosophila melanogaster

Tandemly repeating satellite DNA elements in heterochromatin occupy a substantial portion of many eukaryotic genomes. Although often characterized as genomic parasites deleterious to the host, they also can be crucial for essential processes such as chromosome segregation. Adding to their interest, satellite DNA elements evolve at high rates; among Drosophila, closely related species often differ drastically in both the types and abundances of satellite repeats. However, due to technical challenges, the evolutionary mechanisms driving this rapid turnover remain unclear. Here we characterize natural variation in simple-sequence repeats of 2-10 bp from inbred Drosophila melanogaster lines derived from multiple populations, using a method we developed called k-Seek that analyzes unassembled Illumina sequence reads. In addition to quantifying all previously described satellite repeats, we identified many novel repeats of low to medium abundance. Many of the repeats show population differentiation, including two that are present in only some populations. Interestingly, the population structure inferred from overall satellite quantities does not recapitulate the expected population relationships based on the demographic history of D. melanogaster. We also find that some satellites of similar sequence composition are correlated across lines, revealing concerted evolution. Moreover, correlated satellites tend to be interspersed with each other, further suggesting that concerted change is partially driven by higher order structure. Surprisingly, we identified negative correlations among some satellites, suggesting antagonistic interactions. Our study demonstrates that current genome assemblies vastly underestimate the complexity, abundance, and variation of highly repetitive satellite DNA and presents approaches to understand their rapid evolutionary divergence.

Differential amplification of satellite PaB6 in chromosomally hypervariable Prospero autumnale complex (Hyacinthaceae)

BACKGROUND AND AIMS

Chromosomal evolution, including numerical and structural changes, is a major force in plant diversification and speciation. This study addresses genomic changes associated with the extensive chromosomal variation of the Mediterranean Prospero autumnale complex (Hyacinthaceae), which includes four diploid cytotypes each with a unique combination of chromosome number (x = 5, 6, 7), rDNA loci and genome size.

METHODS

A new satellite repeat PaB6 has previously been identified, and monomers were reconstructed from next-generation sequencing (NGS) data of P. autumnale cytotype B(6)B(6) (2n = 12). Monomers of all other Prospero cytotypes and species were sequenced to check for lineage-specific mutations. Copy number, restriction patterns and methylation levels of PaB6 were analysed using Southern blotting. PaB6 was localized on chromosomes using fluorescence in situ hybridization (FISH).

KEY RESULTS

The monomer of PaB6 is 249 bp long, contains several intact and truncated vertebrate-type telomeric repeats and is highly methylated. PaB6 is exceptional because of its high copy number and unprecedented variation among diploid cytotypes, ranging from 10(4) to 10(6) copies per 1C. PaB6 is always located in pericentromeric regions of several to all chromosomes. Additionally, two lineages of cytotype B(7)B(7) (x = 7), possessing either a single or duplicated 5S rDNA locus, differ in PaB6 copy number; the ancestral condition of a single locus is associated with higher PaB6 copy numbers.

CONCLUSIONS

Although present in all Prospero species, PaB6 has undergone differential amplification only in chromosomally variable P. autumnale, particularly in cytotypes B(6)B(6) and B(5)B(5). These arose via independent chromosomal fusions from x = 7 to x = 6 and 5, respectively, accompanied by genome size increases. The copy numbers of satellite DNA PaB6 are among the highest in angiosperms, and changes of PaB6 are exceptionally dynamic in this group of closely related cytotypes of a single species. The evolution of the PaB6 copy numbers is discussed, and it is suggested that PaB6 represents a recent and highly dynamic system originating from a small pool of ancestral repeats.

Transposable elements as stress adaptive capacitors induce genomic instability in fungal pathogen Magnaporthe oryzae

A fundamental problem in fungal pathogenesis is to elucidate the evolutionary forces responsible for genomic rearrangements leading to races with fitter genotypes. Understanding the adaptive evolutionary mechanisms requires identification of genomic components and environmental factors reshaping the genome of fungal pathogens to adapt. Herein, Magnaporthe oryzae, a model fungal plant pathogen is used to demonstrate the impact of environmental cues on transposable elements (TE) based genome dynamics. For heat shock and copper stress exposed samples, eight TEs belonging to class I and II family were employed to obtain DNA profiles. Stress induced mutant bands showed a positive correlation with dose/duration of stress and provided evidences of TEs role in stress adaptiveness. Further, we demonstrate that genome dynamics differ for the type/family of TEs upon stress exposition and previous reports of stress induced MAGGY transposition has underestimated the role of TEs in M. oryzae. Here, we identified Pyret, MAGGY, Pot3, MINE, Mg-SINE, Grasshopper and MGLR3 as contributors of high genomic instability in M. oryzae in respective order. Sequencing of mutated bands led to the identification of LTR-retrotransposon sequences within regulatory regions of psuedogenes. DNA transposon Pot3 was identified in the coding regions of chromatin remodelling protein containing tyrosinase copper-binding and PWWP domains. LTR-retrotransposons Pyret and MAGGY are identified as key components responsible for the high genomic instability and perhaps these TEs are utilized by M. oryzae for its acclimatization to adverse environmental conditions. Our results demonstrate how common field stresses change genome dynamics of pathogen and provide perspective to explore the role of TEs in genome adaptability, signalling network and its impact on the virulence of fungal pathogens.

Isolation and characterization of two satellite DNAs in some Iberian rock lizards (Squamata, Lacertidae)

Satellite DNAs represent a large portion of all high eukaryotic genomes. They consist of numerous very similar repeated sequences, tandemly arranged in large clusters up to 100 million base pairs in length, usually located in the heterochromatic parts of chromosomes. The biological significance of satDNAs is still under discussion, but most of their proposed functions are related to heterochromatin and/or centromere formation and function. Because information about the structure of reptilian satDNA is far from exhaustive, we present a molecular and cytogenetic characterization of two satDNA families in four lacertid species. Two families of tandemly repeated DNAs, namely TaqI and HindIII satDNAs, have been cloned and sequenced from four species belonging to the genus Iberolacerta. These satDNAs are characterized by a monomer length of 171-188 and 170-172 bp, and by an AT content of 60.5% and 58.1%, respectively. FISH experiments with TaqI satDNA probe produced bright signals in pericentromeric regions of a subset of chromosomes whereas all the centromeres were marked by HindIII probe. The results obtained in this study suggest that chromosome location and abundance of satDNAs influence the evolution of these elements, with centromeric families evolving tenfold faster than interstitial/pericentromeric ones. Such different rates render different satellites useful for phylogenetic investigation at different taxonomic ranks.

Screening and characterization of sex-specific DNA fragments in the freshwater fish matrinchã, Brycon amazonicus (Teleostei: Characiformes: Characidae)

The matrinchã Brycon amazonicus, a commercially important freshwater fish resource, has no heteromorphic sex chromosomes so far described. In the present study, we performed a screening of sex-associated DNA markers in this species, through the use of a random amplified polymorphic DNA (RAPD) assay and a genomic DNA restriction digestion analysis. DNA digestions evidenced no differences between sexes. Sixty-six random primers were used in pooled and individual DNA samples of males and females, and the analysis of the RAPD fingerprints revealed one female sex-associated band. Cloning and sequencing of this band led to the identification of two distinct DNA segments. While one of the isolated fragments showed a significant identity with a described protein gene (phosphatidylinositol glycan anchor biosynthesis, class W), the other fragment, composed of 535 bp, corresponds to a novel DNA marker. Further experiments were performed with this second DNA fragment in order to verify its sex-specificity. Data on dot blot hybridization, using total DNA of both sexes, confirmed its female-specificity in B. amazonicus. A primer set was designed based on its sequence data and used in PCR with DNA samples of this species, leading to diagnose the animals' sexes with a 100 % overall accuracy through a sequence characterized amplified region approach. No amplification results were found for two other species of the genus--B. orbignyanus and B. lundii. The obtained data can lead to the hypothesis that B. amazonicus may present heteromorphic sex chromosomes that should be in an early phase of differentiation.

Genetic diversity and mating type distribution of Tuber melanosporum and their significance to truffle cultivation in artificially planted truffieres in Australia

Tuber melanosporum is a truffle native to Europe and is cultivated in countries such as Australia for the gastronomic market, where production yields are often lower than expected. We assessed the genetic diversity of T. melanosporum with six microsatellite loci to assess the effect of genetic drift on truffle yield in Australia. Genetic diversity as assessed on 210 ascocarps revealed a higher allelic diversity compared to previous studies from Europe, suggesting a possible genetic expansion and/or multiple and diverse source populations for inoculum. The results also suggest that the single sequence repeat diversity of locus ME2 is adaptive and that, for example, the probability of replication errors is increased for this locus. Loss of genetic diversity in Australian populations is therefore not a likely factor in limiting ascocarp production. A survey of nursery seedlings and trees inoculated with T. melanosporum revealed that <70% of seedlings and host trees were colonized with T. melanosporum and that some trees had been contaminated by Tuber brumale, presumably during the inoculation process. Mating type (MAT1-1-1 and MAT1-2-1) analyses on seedling and four- to ten-year-old host trees found that 100% of seedlings but only approximately half of host trees had both mating types present. Furthermore, MAT1-1-1 was detected significantly more commonly than MAT1-2-1 in established trees, suggesting a competitive advantage for MAT1-1-1 strains. This study clearly shows that there are more factors involved in ascocarp production than just the presence of both mating types on host trees.

Structure and population dynamics of the major satellite DNA in the red flour beetle Tribolium castaneum

In the beetle genus Tribolium, satellite DNAs comprise a significant amount of pericentromeric heterochromatin and are characterized by rapid turnover resulting in species specific profiles. In the present work we characterize the major pericentromeric satellite DNA TCAST of the beetle T. castaneum and analyse its population dynamics. Using direct sequencing of genomic PCR products we show that the TCAST satellite exists in the form of two related subfamilies: Tcast1a and Tcast1b that make up 20 and 15% of the genome, respectively. Tcast1a and Tcast1b have consensus sequences of 377 and 362 bp respectively, share an average similarity of 79% and are characterized by a divergent, subfamily specific region of approximately 100 bp. The two subfamilies are prevalently organized in the interspersed form, although a portion exists in the form of homogenous tandem arrays composed of only Tcast1a or Tcast1b. The pattern of restriction enzyme digestion indicates that Tcast1a and Tcast1b are organized in composite higher order repeats. Comparison of sequence variability of Tcast1a and Tcast1b among ten strains reveals a difference in the frequency of particular mutations present at some positions. However, no difference in the organization and in the amount of subfamilies was detected among strains. The results show that direct genomic sequencing can be a useful method for the detection of population specific features of satellite DNA. In the case of TCAST satellite DNA, changes in the mutational profiles seem to represent the first step in the genesis of a population specific satellite profile.

Population genetic structure of Culex quinquefasciatus in India by ISSR marker

OBJECTIVE

To characterize the genetic structure of various populations of Culex quinquefasciatus (Cx. quinquefasciatus) from India representing different geoclimatic locations.

METHODS

Inter simple sequence repeat (ISSR) markers were used. A set of 20 primers were screened with the laboratory populations of mosquito species. Finally the IS 40 primer was chosen based on the scorable banding pattern showing 100 percent polymorphism among the various populations. The statistical analysis was done using POPGENE 1.31 software. The consensus tree was generated based on UPGMA modified from NEIGHBOR procedure of PHYLIP Version 3.5.

RESULTS

The cluster analysis shows the main cluster which is divided into two sub cluster representing all the populations separated as per their phylogeographic and geoclimatic condition.

CONCLUSIONS

The findings will be helpful in understanding the population variation under different ecological conditions and development of effective vector management strategies.

Genome nucleotide composition shapes variation in simple sequence repeats

Simple sequence repeats (SSRs) or microsatellites are a common component of genomes but vary greatly across species in their abundance. We tested the hypothesis that this variation is due in part to AT/GC content of genomes, with genomes biased toward either high AT or high CG generating more short random repeats that are long enough to enhance expansion through slippage during replication. To test this hypothesis, we identified repeats with perfect tandem iterations of 1-6 bp from 25 protists with complete or near-complete genome sequences. As expected, the density and the frequency are highly related to genome AT content, with excellent fits to quadratic regressions with minima near a 50% AT content and rising toward both extremes. Within species, the same trends hold, except the limited variation in AT content within each species places each mainly on the descending (GC rich), middle, or ascending (AT rich) part of the curve. The base usages of repeat motifs are also significantly correlated with genome nucleotide compositions: Percentages of AT-rich motifs rise with the increase of genome AT content but vice versa for GC-rich subgroups. Amino acid homopolymer repeats also show the expected quadratic relationship, with higher abundance in species with AT content biased in either direction. Our results show that genome nucleotide composition explains up to half of the variance in the abundance and motif constitution of SSRs.

The non-regular orbit: three satellite DNAs in Drosophila martensis (buzzatii complex, repleta group) followed three different evolutionary pathways

The genome of species from the buzzatii cluster (buzzatii complex, repleta group) is hosted by a number of satellite DNAs (satDNAs) showing contrasting structural characteristics, genomic organization and evolution, such as pBuM-alpha (~190 bp repeats), pBuM-alpha/beta (~370 bp repeats) and the DBC-150 (~150 bp repeats). In the present study, we aimed to investigate the evolution of these three satDNAs by looking for homologous sequences in the genome of the closest outgroup species: Drosophila martensis (buzzatii complex). After PCR, we isolated and sequenced 9 alpha, 8 alpha/beta and 11 DBC-150 sequences from this species. The results were compared to all pBuM and DBC-150 sequences available in literature. After D. martensis split from the buzzatii cluster some 6 Mya, the three satDNAs evolved differently in the genome of D. martensis by: (1) maintenance of a collection of major types of ancestral repeats in the genome (alpha); (2) fixation for a single major type of ancestral repeats (alpha/beta) or (3) fixation for new divergent species-specific repeat types (DBC-150). Curiously, D. seriema and D. martensis, although belonging to different and allopatric clusters, became independently fixed for the same major type of alpha/beta ancestral repeats, illustrating a rare case of parallelism in satDNA evolution. The contrasting pictures illustrate the diversity of evolutionary pathways a satDNA can follow, defining a "non-regular orbit" with outcomes difficult to predict.

Phylogenetic footprinting of non-coding RNA: hammerhead ribozyme sequences in a satellite DNA family of Dolichopoda cave crickets (Orthoptera, Rhaphidophoridae)

BACKGROUND

The great variety in sequence, length, complexity, and abundance of satellite DNA has made it difficult to ascribe any function to this genome component. Recent studies have shown that satellite DNA can be transcribed and be involved in regulation of chromatin structure and gene expression. Some satellite DNAs, such as the pDo500 sequence family in Dolichopoda cave crickets, have a catalytic hammerhead (HH) ribozyme structure and activity embedded within each repeat.

RESULTS

We assessed the phylogenetic footprints of the HH ribozyme within the pDo500 sequences from 38 different populations representing 12 species of Dolichopoda. The HH region was significantly more conserved than the non-hammerhead (NHH) region of the pDo500 repeat. In addition, stems were more conserved than loops. In stems, several compensatory mutations were detected that maintain base pairing. The core region of the HH ribozyme was affected by very few nucleotide substitutions and the cleavage position was altered only once among 198 sequences. RNA folding of the HH sequences revealed that a potentially active HH ribozyme can be found in most of the Dolichopoda populations and species.

CONCLUSIONS

The phylogenetic footprints suggest that the HH region of the pDo500 sequence family is selected for function in Dolichopoda cave crickets. However, the functional role of HH ribozymes in eukaryotic organisms is unclear. The possible functions have been related to trans cleavage of an RNA target by a ribonucleoprotein and regulation of gene expression. Whether the HH ribozyme in Dolichopoda is involved in similar functions remains to be investigated. Future studies need to demonstrate how the observed nucleotide changes and evolutionary constraint have affected the catalytic efficiency of the hammerhead.

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.

A surrogate approach to study the evolution of noncoding DNA elements that organize eukaryotic genomes

Comparative genomics provides a facile way to address issues of evolutionary constraint acting on different elements of the genome. However, several important DNA elements have not reaped the benefits of this new approach. Some have proved intractable to current day sequencing technology. These include centromeric and heterochromatic DNA, which are essential for chromosome segregation as well as gene regulation, but the highly repetitive nature of the DNA sequences in these regions make them difficult to assemble into longer contigs. Other sequences, like dosage compensation X chromosomal sites, origins of DNA replication, or heterochromatic sequences that encode piwi-associated RNAs, have proved difficult to study because they do not have recognizable DNA features that allow them to be described functionally or computationally. We have employed an alternate approach to the direct study of these DNA elements. By using proteins that specifically bind these noncoding DNAs as surrogates, we can indirectly assay the evolutionary constraints acting on these important DNA elements. We review the impact that such "surrogate strategies" have had on our understanding of the evolutionary constraints shaping centromeres, origins of DNA replication, and dosage compensation X chromosomal sites. These have begun to reveal that in contrast to the view that such structural DNA elements are either highly constrained (under purifying selection) or free to drift (under neutral evolution), some of them may instead be shaped by adaptive evolution and genetic conflicts (these are not mutually exclusive). These insights also help to explain why the same elements (e.g., centromeres and replication origins), which are so complex in some eukaryotic genomes, can be simple and well defined in other where similar conflicts do not exist.

The centromere-drive hypothesis: a simple basis for centromere complexity

Centromeres are far more complex and evolutionarily labile than expected based on their conserved, essential function. The rapid evolution of both centromeric DNA and proteins strongly argue that centromeres are locked in an evolutionary conflict to increase their odds of transmission during asymmetric (female) meiosis. Evolutionary success for "cheating" centromeres can result in highly deleterious consequences for the species, either in terms of skewed sex ratios or male sterility. Centromeric proteins evolve rapidly to suppress the deleterious effects of "centromere-drive." This chapter summarizes the mounting evidence in favor of the centromere-drive model, and its implications for centromere evolution in taxa with variations in meiosis.

Survey of extrachromosomal circular DNA derived from plant satellite repeats

BACKGROUND

Satellite repeats represent one of the most dynamic components of higher plant genomes, undergoing rapid evolutionary changes of their nucleotide sequences and abundance in a genome. However, the exact molecular mechanisms driving these changes and their eventual regulation are mostly unknown. It has been proposed that amplification and homogenization of satellite DNA could be facilitated by extrachromosomal circular DNA (eccDNA) molecules originated by recombination-based excision from satellite repeat arrays. While the models including eccDNA are attractive for their potential to explain rapid turnover of satellite DNA, the existence of satellite repeat-derived eccDNA has not yet been systematically studied in a wider range of plant genomes.

RESULTS

We performed a survey of eccDNA corresponding to nine different families and three subfamilies of satellite repeats in ten species from various genera of higher plants (Arabidopsis, Oryza, Pisum, Secale, Triticum and Vicia). The repeats selected for this study differed in their monomer length, abundance, and chromosomal localization in individual species. Using two-dimensional agarose gel electrophoresis followed by Southern blotting, eccDNA molecules corresponding to all examined satellites were detected. EccDNA occurred in the form of nicked circles ranging from hundreds to over eight thousand nucleotides in size. Within this range the circular molecules occurred preferentially in discrete size intervals corresponding to multiples of monomer or higher-order repeat lengths.

CONCLUSION

This work demonstrated that satellite repeat-derived eccDNA is common in plant genomes and thus it can be seriously considered as a potential intermediate in processes driving satellite repeat evolution. The observed size distribution of circular molecules suggests that they are most likely generated by molecular mechanisms based on homologous recombination requiring long stretches of sequence similarity.

Microsatellite accumulation on the Y chromosome in Silene latifolia

The dioecious plant Silene latifolia possesses evolutionarily young sex chromosomes, and so serves as a model system to study the early stages of sex chromosome evolution. Sex chromosomes often differ distinctly from autosomes in both their structure and their patterns of evolution. The S. latifolia Y chromosome is particularly unique owing to its large size, which contrasts with the size of smaller, degenerate mammalian Y chromosomes. It is thought that the suppression of recombination on the S. latifolia Y chromosome could have resulted in the accumulation of repetitive sequences that account for its large size. Here we used fluorescence in situ hybridization (FISH) to study the chromosomal distribution of various microsatellites in S. latifolia including all possible mono-, di-, and tri-nucleotides. Our results demonstrate that a majority of microsatellites are accumulated on the q arm of the Y chromosome, which stopped recombining relatively recently and has had less time to accumulate repetitive DNA sequences compared with the p arm. Based on these results we can speculate that microsatellites have accumulated in regions that predate the genome expansion, supporting the view that the accumulation of repetitive DNA sequences occurred prior to, not because of, the degeneration of genes.

Comparative analyses of human single- and multilocus tandem repeats

Using the compiled human genome sequence, we systematically cataloged all tandem repeats with periods between 20 and 2000 bp and defined two subsets whose consensus sequences were found at either single-locus tandem repeats (slTRs) or multilocus tandem repeats (mlTRs). Parameters compiled for these subsets provide insights into mechanisms underlying the creation and evolution of tandem repeats. Both subsets of tandem repeats are nonrandomly distributed in the genome, being found at higher frequency at many but not all chromosome ends and internal clusters of mlTRs were also observed. Despite the integral role of recombination in the biology of tandem repeats, recombination hotspots colocalized only with shorter microsatellites and not the longer repeats examined here. An increased frequency of slTRs was observed near imprinted genes, consistent with a functional role, while both slTRs and mlTRs were found more frequently near genes implicated in triplet expansion diseases, suggesting a general instability of these regions. Using our collated parameters, we identified 2230 slTRs as candidates for highly informative molecular markers.

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.

Analysis of Drosophila species genome size and satellite DNA content reveals significant differences among strains as well as between species

The size of eukaryotic genomes can vary by several orders of magnitude, yet genome size does not correlate with the number of genes nor with the size or complexity of the organism. Although "whole"-genome sequences, such as those now available for 12 Drosophila species, provide information about euchromatic DNA content, they cannot give an accurate estimate of genome sizes that include heterochromatin or repetitive DNA content. Moreover, genome sequences typically represent only one strain or isolate of a single species that does not reflect intraspecies variation. To more accurately estimate whole-genome DNA content and compare these estimates to newly assembled genomes, we used flow cytometry to measure the 2C genome values, relative to Drosophila melanogaster. We estimated genome sizes for the 12 sequenced Drosophila species as well as 91 different strains of 38 species of Drosophilidae. Significant differences in intra- and interspecific 2C genome values exist within the Drosophilidae. Furthermore, by measuring polyploid 16C ovarian follicle cell underreplication we estimated the amount of satellite DNA in each of these species. We found a strong correlation between genome size and amount of satellite underreplication. Addition and loss of heterochromatin satellite repeat elements appear to have made major contributions to the large differences in genome size observed in the Drosophilidae.

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.

Low rates of homogenization of the DBC-150 satellite DNA family restricted to a single pair of microchromosomes in species from the Drosophila buzzatii cluster

A satellite DNA family, termed DBC-150, comprises slightly GC-rich repeat units of approximately 150 bp that were isolated (by DNA digestions or PCR) from the genome of all seven Drosophila species from the buzzatii cluster (repleta group). The presence of subrepeats suggests that part of the extant DBC-150 monomer originated by the duplication of small sequence motifs. The DBC-150 family is compared to the previously described pBuM satDNA family, an abundant component of the genome of five species of the cluster. The two families are different in several aspects, including primary structure, A + T content, intraspecific and interspecific variability and rates of homogenization (or nucleotide spread). The data indicate a lower rate of homogenization (and absence of complete concerted evolution) of the DBC-150 compared to the pBuM family. FISH on metaphase chromosomes revealed that the DBC-150 family is located exclusively in the microchromosomes. To our knowledge this is the first record of a complex Drosophila satDNA restricted to a single pair of microchromosomes. The observed low rates of homogenization of the DBC-150 family might be related to a presumed reduction or suppression of meiotic recombination in the microchromosomes.

Informativeness of dinucleotide repeat-based primers in fungal pathogen of rice Magnaporthe grisea

Microsatellites or simple sequence repeats (SSRs) occur ubiquitously and show complex patterns in length, motif size and sequence. Among SSRs, dinucleotide repeats occur in high abundance in fungi with shorter length as compared to other organisms. In this study, multilocus profiles obtained in Magnaporthe grisea, a model plant pathogen were evaluated. The results showed lower rate of polymorphism by (GT)(n)/(TG)(n) repeat-based primers and suggested occurrence of (GA)(n)/(AG)(n) repeats as integral repeats and (TC)(n)/(CT)(n) and (AC)(n)/(CA)(n) as non-integral repeats. Low repeat length variation was found to be correlated with less number of repeat motifs. The study provides an insight into the possibility of molecular coevolution of mobile elements and dinucleotide repeats in fungi. The study could be applied to other species for wider applications including evolutionary and population genetics.

The rise, fall and renaissance of microsatellites in eukaryotic genomes

Microsatellites are among the most versatile of genetic markers, being used in an impressive number of biological applications. However, the evolutionary dynamics of these markers remain a source of contention. Almost 20 years after the discovery of these ubiquitous simple sequences, new genomic data are clarifying our understanding of the structure, distribution and variability of microsatellites in genomes, especially for the eukaryotes. While these new data provide a great deal of descriptive information about the nature and abundance of microsatellite sequences within eukaryotic genomes, there have been few attempts to synthesise this information to develop a global concept of evolution. This review provides an up-to-date account of the mutational processes, biases and constraints believed to be involved in the evolution of microsatellites, particularly with respect to the creation and degeneration of microsatellites, which we assert may be broadly viewed as a life cycle. In addition, we identify areas of contention that require further research and propose some possible directions for future investigation.

Characterization of a novel satellite DNA sequence from Flying Dragon (Poncirus trifoliata)

Repetitive sequences constitute a significant component of most eukaryotic genomes, and the isolation and characterization of repetitive DNA sequences provide an insight into the organization of the genome of interest. Here, we report the isolation and the molecular analysis and methylation status of a novel tandemly organized repetitive DNA sequence from the genome of Poncirus trifoliata. Digestion of P. trifoliata DNA with Afa I produced a prominent fragment of approximately 400 bp. Southern blotting analysis of genomic DNA digested with the same enzyme revealed a ladder composed of DNA fragments that are multimers of the 400-bp Afa I band, indicating that the repetitive DNA is arrayed in tandem. This suggests that Afa I isolated a novel satellite that we have called Poncirus trifoliata satellite DNA 400 (PN400). This satellite composes 25% of the genome and it is also present in lemon, sour orange and kumquat. Analysis of the methylation status demonstrated that the cytosines in CCGG sequences in this satellite were methylated.

Genetic conflicts during meiosis and the evolutionary origins of centromere complexity

Centromeric DNA evolves rapidly, ranging in size and complexity over several orders of magnitude. Traditional attempts at studying centromeres have left unexplained the causes underlying this complexity and rapid evolution. Instead of directly studying centromeric DNA sequence, our approach has been to study the proteins that epigenetically determine centromere identity. We have discovered that centromeric histones (CenH3s) have evolved under positive selection in multiple lineages, suggesting an involvement in recurrent genetic conflict. Our hypothesis is that 'centromere-drive' is the source of this conflict. Under this model, centromeres compete via microtubule attachments for preferential transmission in female meioses occurring in animals and plants. Since only one of four meiotic products will become the egg, this competition confers a selfish advantage to chromosomes that can make more microtubule attachments, resulting in runaway expansions of centromeric satellites. While beneficial to the 'driving' chromosome, these expansions can have deleterious effects on the fitness of an organism and of the species. CenH3s as well as other heterochromatin proteins have evolved under positive selection to suppress the deleterious consequences of 'centromere-drive' by restoring meiotic parity.

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.