Directional evolution in germline microsatellite mutations

Orthologous microsatellites, transposable elements, and DNA deletions correlate with generation time and body mass in neoavian birds

The rate of mutation accumulation in germline cells can be affected by cell replication and/or DNA damage, which are further related to life history traits such as generation time and body mass. Leveraging the existing datasets of 233 neoavian bird species, here, we investigated whether generation time and body mass contribute to the interspecific variation of orthologous microsatellite length, transposable element (TE) length, and deletion length and how these genomic attributes affect genome sizes. In nonpasserines, we found that generation time is correlated to both orthologous microsatellite length and TE length, and body mass is negatively correlated to DNA deletions. These patterns are less pronounced in passerines. In all species, we found that DNA deletions relate to genome size similarly as TE length, suggesting a role of body mass dynamics in genome evolution. Our results indicate that generation time and body mass shape the evolution of genomic attributes in neoavian birds.

Genome-wide characterization of microsatellite DNA in fishes: survey and analysis of their abundance and frequency in genome-specific regions

BACKGROUND

Microsatellite repeats are ubiquitous in organism genomes and play an important role in the chromatin organization, regulation of gene activity, recombination and DNA replication. Although microsatellite distribution patterns have been studied in most phylogenetic lineages, they are unclear in fish species.

RESULTS

Here, we present the first systematic examination of microsatellite distribution in coding and non-coding regions of 14 fish genomes. Our study showed that the number and type of microsatellites displayed nonrandom distribution for both intragenic and intergenic regions, suggesting that they have potential roles in transcriptional or translational regulation and DNA replication slippage theories alone were insufficient to explain the distribution patterns. Our results showed that microsatellites are dominant in non-coding regions. The total number of microsatellites ranged from 78,378 to 1,012,084, and the relative density varied from 4925.76 bp/Mb to 25,401.97 bp/Mb. Overall, (A + T)-rich repeats were dominant. The dependence of repeat abundance on the length of the repeated unit (1-6 nt) showed a great similarity decrease, whereas more tri-nucleotide repeats were found in exonic regions than tetra-nucleotide repeats of most species. Moreover, the incidence of different repeated types appeared species- and genomic-specific. These results highlight potential mechanisms for maintaining microsatellite distribution, such as selective forces and mismatch repair systems.

CONCLUSIONS

Our data could be beneficial for the studies of genome evolution and microsatellite DNA evolutionary dynamics, and facilitate the exploration of microsatellites structural, function, composition mode and molecular markers development in these species.

Twenty Years Later: A Comprehensive Review of the X Chromosome Use in Forensic Genetics

The unique structure of the X chromosome shaped by evolution has led to the present gender-specific genetic differences, which are not shared by its counterpart, the Y chromosome, and neither by the autosomes. In males, recombination between the X and Y chromosomes is limited to the pseudoautosomal regions, PAR1 and PAR2; therefore, in males, the X chromosome is (almost) entirely transmitted to female offspring. On the other hand, the X chromosome is present in females with two copies that recombine along the whole chromosome during female meiosis and that is transmitted to both female and male descendants. These transmission characteristics, besides the obvious clinical impact (sex chromosome aneuploidies are extremely frequent), make the X chromosome an irreplaceable genetic tool for population genetic-based studies as well as for kinship and forensic investigations. In the early 2000s, the number of publications using X-chromosomal polymorphisms in forensic and population genetic applications increased steadily. However, nearly 20 years later, we observe a conspicuous decrease in the rate of these publications. In light of this observation, the main aim of this article is to provide a comprehensive review of the advances and applications of X-chromosomal markers in population and forensic genetics over the last two decades. The foremost relevant topics are addressed as: (i) developments concerning the number and types of markers available, with special emphasis on short tandem repeat (STR) polymorphisms (STR nomenclatures and practical concerns); (ii) overview of worldwide population (frequency) data; (iii) the use of X-chromosomal markers in (complex) kinship testing and the forensic statistical evaluation of evidence; (iv) segregation and mutation studies; and (v) current weaknesses and future prospects.

Infestation Pattern and Population Dynamics of the Tropical Bed Bug, Cimex hemipterus (F.) (Hemiptera: Cimicidae) Based on Novel Microsatellites and mtDNA Markers

The tropical bed bug, Cimex hemipterus (F.), has now emerged as an important public health pest in the tropics. Despite its alarming infestation rate, the information on its population genetics remains scarce. Here, we described the infestation structure and population dynamics of C. hemipterus in the tropics, especially Malaysia and Singapore, based on eight novel microsatellites and two mtDNA markers, including cytochrome c oxidase I (COI) and 16S rRNA genes. Across populations, microsatellite data revealed high genetic diversity with significant genetic differentiation and restricted gene flow. Analysis within populations revealed evidence of a recent bottleneck. Nonetheless, elevated genetic diversity in nearly all populations suggests that the propagule in C. hemipterus populations were much diverse, distantly related (mean r = 0.373), and not significantly inbred (mean F_IS = 0.24) than that observed in Cimex lectularius from previous studies. We observed seven mtDNA haplotypes across the 18 populations studied (Hd = 0.593) and several populations displayed more than one matrilineal descent. The two markers were generally congruent in suggesting a common, genetically diverse (especially at the nuclear region) source population with possibilities of multiple introductions for the bed bug populations in the present study.

Molecular Population Genetics of the Northern Elephant Seal Mirounga angustirostris

The northern elephant seal, Mirounga angustirostris, was heavily hunted and declared extinct in the 19th century. However, a colony remained on remote Guadalupe Island, Mexico and the species has since repopulated most of its historical distribution. Here, we present a comprehensive evaluation of genetic variation in the species. First, we assess the effect of the demographic bottleneck on microsatellite variability and compare it with that found in other pinnipeds, demonstrating levels of variation similar to that in species that continue to be threatened with extinction. Next, we use sequence data from these markers to demonstrate that some of the limited polymorphism predates the bottleneck. However, most contemporary variation appears to have arisen recently and persisted due to exponential growth. We also describe how we use the range in allele size of microsatellites to estimate ancestral effective population size before the bottleneck, demonstrating a large reduction in effective size. We then employ a classical method for bacteria to estimate the microsatellite mutation rate in the species, deriving an estimate that is extremely similar to that estimated for a similar set of loci in humans, indicating consistency of microsatellite mutation rates in mammals. Finally, we find slight significant structure between some geographically separated colonies, although its biological significance is unclear. This work demonstrates that genetic analysis can be useful for evaluating the population biology of the northern elephant seal, in spite of the bottleneck that removed most genetic variation from the species.

A new stochastic model of microsatellite evolution

We introduce a continuous-time Markov chain model for the evolution of microsatellites, simple sequence repeats in DNA. We prove the existence of a unique stationary distribution for our model, and fit the model to data from approximately 106base pairs of DNA from fruit flies, mice, and humans. The slippage rates from the best fit for our model are consistent with experimental findings.

A new stochastic model of microsatellite evolution

We introduce a continuous-time Markov chain model for the evolution of microsatellites, simple sequence repeats in DNA. We prove the existence of a unique stationary distribution for our model, and fit the model to data from approximately 106 base pairs of DNA from fruit flies, mice, and humans. The slippage rates from the best fit for our model are consistent with experimental findings.

Genetic Variability and Geographic Diversity of the Common Bed Bug (Hemiptera: Cimicidae) Populations from the Midwest Using Microsatellite Markers

With the recent global resurgence of the bed bugs (Cimex lectularius L.), there is a need to better understand its biology, ecology, and ability to establish populations. Bed bugs are domestic pests that feed mainly on mammalian blood. Although bed bugs have not been implicated as vectors of pathogens, their biting activity inflicts severe insomnia and allergic reactions. Moreover, they have recently developed resistance to various insecticides, which requires further molecular research to determine genetic variation and appropriate interventions. Population dynamics, including genetic differentiation and genetic distance of 10 populations from the Midwest were analyzed in this study. The bed bug samples collected by pest control companies were genotyped using eight species-specific microsatellite markers. Results showed all eight markers were polymorphic, with 8-16 alleles per locus, suggesting high genetic diversity. The FST values were >0.25, signifying pronounced genetic differentiation. The G-test results also indicated high genetic differentiation among populations. The frequency of the most common allele across all eight loci was 0.42. The coefficient of relatedness between each of the populations was >0.5, indicative of sibling or parent-offspring relationships, while the FIS and its confidence interval values were statistically insignificant within the populations tested. The populations departed from Hardy-Weinberg equilibrium, possibly because of high heterozygosity. The genetic distance analysis using a neighbor-joining tree showed that the populations from Kansas City, MO, were genetically separate from most of those from Nebraska, indicating a geographic pattern of genetic structure. Our study demonstrated the effectiveness of using microsatellite markers to study bed bugs population structure, thereby improving our understanding of bed bug population dynamics in the Midwest. Overall, this study showed a high genetic diversity and identified several new alleles in the bed bug populations in the Midwest.

Single-molecule PCR analysis of an unstable microsatellite for detecting mutations in sperm of mice exposed to chemical mutagens

Single-molecule PCR (SM-PCR) analysis of long and repetitive DNA sequences, known as expanded simple tandem repeats (ESTRs), has been the most efficient method for studying germline mutation induction in endogenous sequences to date. However, the long length of these sequences makes mutation detection imprecise and laborious, and they have been characterized only in mice. Here, we explore the use of unstable microsatellite sequences that can be typed with high precision by capillary electrophoresis as alternative loci for detecting germline mutations. We screened 24 microsatellite loci across inbred mouse strains and identified Mm2.2.1 as the most polymorphic microsatellite locus. We then optimized SM-PCR of Mm2.2.1 to detect mutations in sperm. SM-PCR analysis of sperm from untreated B6C3F1 and Muta(™)Mouse samples revealed mutation frequencies that are consistent with rates derived from family pedigree analysis (∼ 5 × 10(-3)). To determine whether this locus can be used to detect chemically induced germline mutations, Muta(™)Mouse males were exposed by oral gavage to a single dose of 100mg/kg of N-ethyl-N-nitrosourea (ENU) or to 100mg/kg of benzo(a)pyrene (BaP) for 28 days alongside vehicle treated controls. Sperm were collected 10 weeks post-ENU exposure to sample sperm exposed as spermatogonial stem cells and 6 weeks post-BaP exposure to sample sperm that were dividing spermatogonia when the exposure was terminated. Both treatments resulted in a significant (approximately 2-fold) increase in mutation frequency in sperm compared to the control animals. The work establishes the utility of this microsatellite for studying mutation induction in the germ cells of mice. Because microsatellites are found in virtually every species, this approach holds promise for other organisms, including humans.

Recent insertion/deletion (reINDEL) mutations: increasing awareness to boost molecular-based research in ecology and evolution

Today, the comparative analysis of DNA molecules mainly uses information inferred from nucleotide substitutions. Insertion/deletion (INDEL) mutations, in contrast, are largely considered uninformative and discarded, due to our lacking knowledge on their evolution. However, including rather than discarding INDELs would be relevant to any research area in ecology and evolution that uses molecular data. As a practical approach to better understanding INDEL evolution in general, we propose the study of recent INDEL (reINDEL) mutations – mutations where both ancestral and derived state are seen in the sample. The precondition for reINDEL identification is knowledge about the pedigree of the individuals sampled. Sound reINDEL knowledge will allow the improved modeling needed for including INDELs in the downstream analysis of molecular data. Both microsatellites, currently still the predominant marker system in the analysis of populations, and sequences generated by next-generation sequencing, a promising and rapidly developing range of technologies, offer the opportunity for reINDEL identification. However, a 2013 sample of animal microsatellite studies contained unexpectedly few reINDELs identified. As most likely explanation, we hypothesize that reINDELs are underreported rather than absent and that this underreporting stems from common reINDEL unawareness. If our hypothesis applies, increased reINDEL awareness should allow gathering data rapidly. We recommend the routine reporting of either the absence or presence of reINDELs together with standardized key information on the nature of mutations when they are detected and the use of the keyword “reINDEL” to increase visibility in both instances of successful and unsuccessful search.

Microsatellite frequencies vary with body mass and body temperature in mammals, suggesting correlated variation in mutation rate

Substitution rate is often found to correlate with life history traits such as body mass, a predictor of population size and longevity, and body temperature. The underlying mechanism is unclear but most models invoke either natural selection or factors such as generation length that change the number of mutation opportunities per unit time. Here we use published genome sequences from 69 mammals to ask whether life history traits impact another form of genetic mutation, the high rates of predominantly neutral slippage in microsatellites. We find that the length-frequency distributions of three common dinucleotide motifs differ greatly between even closely related species. These frequency differences correlate with body mass and body temperature and can be used to predict the phenotype of an unknown species. Importantly, different length microsatellites show complicated patterns of excess and deficit that cannot be explained by a simple model where species with short generation lengths have experienced more mutations. Instead, the patterns probably require changes in mutation rate that impact alleles of different length to different extents. Body temperature plausibly influences mutation rate by modulating the propensity for slippage. Existing hypotheses struggle to account for a link between body mass and mutation rate. However, body mass correlates inversely with population size, which in turn predicts heterozygosity. We suggest that heterozygote instability, HI, the idea that heterozygous sites show increased mutability, could provide a plausible link between body mass and mutation rate.

A call to arms: on refining Plasmodium vivax microsatellite marker panels for comparing global diversity

BACKGROUND

Microsatellite (MS) markers have become an important tool for studying the population diversity, evolutionary history and multiplicity of infection (MOI) of malaria parasite infections. MS are typically selected on the basis of being highly polymorphic. However, it is known that the polymorphic potential (mutability) of each marker can vary as much as two orders of magnitude, which radically changes how diversity is represented in the genome from one marker to the next. Over the past decade, approximately 240 Plasmodium vivax MS have been published, comprising nine major panels of markers. Inconsistent usage of each panel has resulted in a surfeit of descriptive genetic diversity data that are largely incomparable between populations. The objective of this study was to statistically evaluate the quality of individual MS markers in order to validate a refined panel of markers that will provide a balanced picture of P. vivax population diversity.

METHODS

All previously published data, including genetic diversity indices, MS parameters, and population parameters, were assembled from 18 different global studies into a flat file to facilitate statistical analysis and modelling using JMP® Genomics 6.0 (SAS Institute Inc, Cary, NC, USA). Statistical modeling was employed to down-select markers with extreme variation among the mean number of alleles, expected heterozygosity, maximum repeat length and/or chromosomal location of the repeat. Individual MS were analysed by step-down whole model linear regression and standard least squares fit models, both stratified by annual parasite incidence to identify MS markers with values significantly different from the mean.

RESULTS

Of the 42 MS under evaluation in this study, 18 (nine high priority) were identified as ideal candidates for measuring population diversity between global regions, while five (two high priority) additional markers were identified as candidates for MOI studies.

CONCLUSIONS

MS diversity was found to be a function of endemicity and motif structure. Evaluation of individual MS permitted the assembly of a refined panel of markers that can be reliably utilized in the field to compare population structures between global regions.

Evolution and homoplasy at the Bem6 microsatellite locus in three sweetpotato whitefly (Bemisia tabaci) cryptic species

BACKGROUND

The evolution of individual microsatellite loci is often complex and homoplasy is common but often goes undetected. Sequencing alleles at a microsatellite locus can provide a more complete picture of the common evolutionary mechanisms occurring at that locus and can reveal cases of homoplasy. Within species homoplasy can lead to an underestimate of differentiation among populations and among species homoplasy can produce a misleading interpretation regarding shared alleles and hybridization. This is especially problematic with cryptic species.

RESULTS

By sequencing alleles from three cryptic species of the sweetpotato whitefly (Bemisia tabaci), designated MEAM1, MED, and NW, the evolution of the putatively dinucleotide Bem6 (CA₈)imp microsatellite locus is inferred as one of primarily stepwise mutation occurring at four distinct heptaucleotide tandem repeats. In two of the species this pattern yields a compound tandem repeat. Homoplasy was detected both among species and within species.

CONCLUSIONS

In the absence of sequencing, size homoplasious alleles at the Bem6 locus lead to an overestimate of alleles shared and hybridization among cryptic species of Bemisia tabaci. Furthermore, the compound heptanucleotide motif structure of a putative dinucleotide microsatellite has implications for the nomenclature of heptanucleotide tandem repeats with step-wise evolution.

Tandem-repeat patterns and mutation rates in microsatellites of the nematode model organism Pristionchus pacificus

Modern evolutionary biology requires integrative approaches that combine life history, population structure, ecology, and development. The nematode Pristionchus pacificus has been established as a model system in which these aspects can be studied in one organism. P. pacificus has well-developed genetic, genomic, and transgenic tools and its ecologic association with scarab beetles is well described. A recent study provided first mutation rate estimates based on mitochondrial genome sequencing and mutation accumulation line experiments that help resolve rather ancient evolutionary branches. Here, we analyzed the tandem-repeat pattern and studied spontaneous mutation rates for microsatellite markers by using the previously generated mutation accumulation lines. We found that 0.59%–3.83% of the genome is composed of short tandem repeats. We developed 41 microsatellite markers, randomly chosen throughout the genome and analyzed them in 82 mutation accumulation lines after 142 generations. A total of 31 mutations were identified in these lines. There was a strong correlation between allele size and mutation rate in P. pacificus, similar to Caenorhabditis elegans. In contrast to C. elegans, however, there is no evidence for a bias toward multistep mutations. The mutation spectrum of microsatellite loci in P. pacificus shows more insertions than deletions, indicating a tendency toward lengthening, a process that might have contributed to the increase in genome size. The mutation rates obtained for individual microsatellite markers provide guidelines for divergence time estimates that can be applied in P. pacificus next-generation sequencing approaches of wild isolates.

Evolution of microsatellite Loci of tropical and temperate anguilla eels

Anguilla eels are divided into temperate and tropical eels, based on their major distributions. The present study collected two temperate eels, Anguilla japonica and Anguilla anguilla, and two tropical eels, Anguilla marmorata and Anguilla bicolor pacifica, to examine two questions: do temperate and tropical Anguilla eels have different genetic polymorphic patterns?; and do temperate Anguilla japonica and Anguilla anguilla have a closer relationship to each other than to tropical eels? In total, 274 sequences were cloned and sequenced from six conserved microsatellite loci to examine polymorphic patterns of these four catadromous eels. Different mutational events, including substitutions, and repeat-unit deletions and insertions, appeared in major regions, while different point mutations were observed in flanking regions. The results implied that parallel patterns of microsatellite sequences occurred within both tropical and temperate freshwater eels. Consensus flanking sequences of six homologous loci from each of the four species were constructed. Genetic distances ranged from 0.044 (Anguilla bicolor pacifica vs. Anguilla marmorata) to 0.061 (Anguilla marmorata vs. Anguilla anguilla). The tree topology suggests the hypothesis of Anguilla japonica and Anguilla anguilla being a sister group must be rejected.

Molecular markers reveal infestation dynamics of the bed bug (Hemiptera: Cimicidae) within apartment buildings

The bed bug, Cimex lectularius L. (Hemiptera: Cimicidae), has experienced an extraordinary global resurgence in recent years, the reasons for which remain poorly understood. Once considered a pest of lower socioeconomic classes, bed bugs are now found extensively across all residential settings, with widespread infestations established in multiapartment buildings. Within such buildings, understanding the population genetic structure and patterns of dispersal may prove critical to the development of effective control strategies. Here, we describe the development of 24 high-resolution microsatellite markers through next generation 454 pyrosequencing and their application to elucidate infestation dynamics within three multistory apartment buildings in the United States. Results reveal contrasting characteristics potentially representative of geographic or locale differences. In Raleigh, NC, an infestation within an apartment building seemed to have started from a single introduction followed by extensive spread. In Jersey City, NJ, two or more introductions followed by spread are evident in two buildings. Populations within single apartments in all buildings were characterized by high levels of relatedness and low levels of diversity, indicative of foundation from small, genetically depauperate propagules. Regardless of the number of unique introductions, genetic data indicate that spread within buildings is extensive, supporting both active and human-mediated dispersal within and between adjacent rooms or apartments spanning multiple floors.

A study on mutational dynamics of simple sequence repeats in relation to mismatch repair system in prokaryotic genomes

Mutational bias toward expansion or contraction of simple sequence repeats (SSRs) is referred to as directionality of SSR evolution. In this communication, we report the mutational bias exhibited by mononucleotide SSRs occurring in the non-coding regions of several prokaryotic genomes. Our investigations revealed that the strains or species lacking mismatch repair (MMR) system generally show higher number of polymorphic SSRs than those species/strains having MMR system. An exception to this observation was seen in the mycobacterial genomes that are MMR deficient where only a few SSR tracts were seen with mutations. This low incidence of SSR mutations even in the MMR-deficient background could be attributed to the high fidelity of the DNA polymerases as a consequence of high generation time of the mycobacteria. MMR system-deficient species generally did not show any bias toward mononucleotide SSR expansions or contractions indicating a neutral evolution of SSRs in these species. The MMR-proficient species in which the observed mutations correspond to secondary mutations showed bias toward contraction of polymononucleotide tracts, perhaps, indicating low efficiency of MMR system to repair SSR-induced slippage errors on template strands. This bias toward deletion in the mononucleotide SSR tracts might be a probable reason behind scarcity for long poly A|T and G|C tracts in prokaryotic systems which are mostly MMR proficient. In conclusion, our study clearly demonstrates mutational dynamics of SSRs in relation to the presence/absence of MMR system in the prokaryotic system.

Isolating microsatellite loci: looking back, looking ahead

Microsatellite DNA loci are tandemly repeated simple sequence repeats (SSRs) that are ubiquitous in eukaryotic genomes. When flanked by unique sequences, length variation (driven by high rates of strand slippage during DNA replication) at a given repeat locus can be assayed by PCR and electrophoretic separation of the resulting DNA fragments (representing alleles defined by fragment size or repeat number at that locus). In nonmodel organisms that do not have sequence information at SSR loci (or at SSRs in a closely related taxon), microsatellites must be isolated and sequenced de novo. Traditionally, this has been accomplished with cloning of genomic DNA fragments enriched for SSRs, a protocol described in detail here. PCR primers flanking microsatellite repeats can be used to assay repeat length variation among individuals (typically through fluorescent labeling of one strand and capillary electrophoresis), useful for questions related to population variation, individual assignment, mating studies, selection scans, mapping, and phenotypic traits. High-throughput next-generation sequencing will likely supplant traditional cloning methods for the discovery of microsatellite loci.

Mutation biases and mutation rate variation around very short human microsatellites revealed by human-chimpanzee-orangutan genomic sequence alignments

I have studied mutation patterns around very short microsatellites, focusing mainly on sequences carrying only two repeat units. By using human-chimpanzee-orangutan alignments, inferences can be made about both the relative rates of mutations and which bases have mutated. I find remarkable non-randomness, with mutation rate depending on a base's position relative to the microsatellite, the identity of the base itself and the motif in the microsatellite. Comparing the patterns around AC2 with those around other four-base combinations reveals that AC2 does not stand out as being special in the sense that non-repetitive tetramers also generate strong mutation biases. However, comparing AC2 and AC3 with AC4 reveals a step change in both the rate and nature of mutations occurring, suggesting a transition state, AC4 exhibiting an alternating high-low mutation rate pattern consistent with the sequence patterning seen around longer microsatellites. Surprisingly, most changes in repeat number occur through base substitutions rather than slippage, and the relative probability of gaining versus losing a repeat in this way varies greatly with repeat number. Slippage mutations reveal rather similar patterns of mutability compared with point mutations, being rare at two repeats where most cause the loss of a repeat, with both mutation rate and the proportion of expansion mutations increasing up to 6-8 repeats. Inferences about longer repeat tracts are hampered by uncertainties about the proportion of multi-species alignments that fail due to multi-repeat mutations and other rearrangements.

DNA slippage occurs at microsatellite loci without minimal threshold length in humans: a comparative genomic approach

The dynamics of microsatellite, or short tandem repeats (STRs), is well documented for long, polymorphic loci, but much less is known for shorter ones. For example, the issue of a minimum threshold length for DNA slippage remains contentious. Model-fitting methods have generally concluded that slippage only occurs over a threshold length of about eight nucleotides, in contradiction with some direct observations of tandem duplications at shorter repeated sites. Using a comparative analysis of the human and chimpanzee genomes, we examined the mutation patterns at microsatellite loci with lengths as short as one period plus one nucleotide. We found that the rates of tandem insertions and deletions at microsatellite loci strongly deviated from background rates in other parts of the human genome and followed an exponential increase with STR size. More importantly, we detected no lower threshold length for slippage. The rate of tandem duplications at unrepeated sites was higher than expected from random insertions, providing evidence for genome-wide action of indel slippage (an alternative mechanism generating tandem repeats). The rate of point mutations adjacent to STRs did not differ from that estimated elsewhere in the genome, except around dinucleotide loci. Our results suggest that the emergence of STR depends on DNA slippage, indel slippage, and point mutations. We also found that the dynamics of tandem insertions and deletions differed in both rates and size at which these mutations take place. We discuss these results in both evolutionary and mechanistic terms.

Origins and spread of pfdhfr mutant alleles in Plasmodium falciparum

The emergence and spread of Plasmodium falciparum parasite resistant to sulfadoxine and pyrimethamine (SP) poses a serious public health problem. Resistance is caused by point mutations in dihydrofolate reductase (pfdhfr) and dihydropteroate synthase (pfdhps), the two key enzymes in the folate biosynthetic pathway. The use of microsatellite markers flanking pfdhfr has recently shown that the invasion of limited resistant lineages may explain the widespread SP resistance in many endemic regions. In Africa, however, multiple indigenous origins of pfdhfr triple mutants have been demonstrated. More new independent lineages and routes of geographical spread of resistance may be found by further molecular evolutionary analyses using samples from various endemic regions. Here, I review recent studies about the history of SP usage and the evolution and spread of resistant lineages while addressing the technical issue of microsatellite analysis.

Sequence determinants of human microsatellite variability

BACKGROUND

Microsatellite loci are frequently used in genomic studies of DNA sequence repeats and in population studies of genetic variability. To investigate the effect of sequence properties of microsatellites on their level of variability we have analyzed genotypes at 627 microsatellite loci in 1,048 worldwide individuals from the HGDP-CEPH cell line panel together with the DNA sequences of these microsatellites in the human RefSeq database.

RESULTS

Calibrating PCR fragment lengths in individual genotypes by using the RefSeq sequence enabled us to infer repeat number in the HGDP-CEPH dataset and to calculate the mean number of repeats (as opposed to the mean PCR fragment length), under the assumption that differences in PCR fragment length reflect differences in the numbers of repeats in the embedded repeat sequences. We find the mean and maximum numbers of repeats across individuals to be positively correlated with heterozygosity. The size and composition of the repeat unit of a microsatellite are also important factors in predicting heterozygosity, with tetra-nucleotide repeat units high in G/C content leading to higher heterozygosity. Finally, we find that microsatellites containing more separate sets of repeated motifs generally have higher heterozygosity.

CONCLUSIONS

These results suggest that sequence properties of microsatellites have a significant impact in determining the features of human microsatellite variability.

Genetic structure and differentiation of Oryza sativa L. in China revealed by microsatellites

China is one of the largest centers of genetic diversity of Oryza sativa L. in the world. Using a genetically representative primary core collection of 3,024 rice landraces in China, we analyzed the genetic structure and intraspecific differentiation of O. sativa, and the directional evolution of SSR. The genetic structure was investigated by model-based structure analysis and construction of neighbor-joining phylogenetic tree. Comparison between genetic structure and predefined populations according to Ting's taxonomic system revealed a hierarchical genetic structure: two distinct subspecies, each with three ecotypes and different numbers of geo-ecogroups within each ecotype. Two subspecies evidently resulted from adaptation to different environments. The different cropping systems imposed on the subspecies led to further differentiation, but the variation within each subspecies resulted from different causes. Indica, under tropical-like or lowland-like environments, exhibited clear differentiation among seasonal ecotypes, but not among soil-watery ecotypes; and japonica showed clear differences between soil water regime ecotypes, but not among seasonal ecotypes. Chinese cultivated rice took on evident directional evolution in microsatellite allele size at several aspects, such as subspecies and geographical populations. Japonica has smaller allele sizes than indica, and this may partly be the result of their different domestication times. Allele size was also negatively correlated with latitude and altitude, and this may be interpreted by different mutation rates, selection pressures, and population size effects under different environments and cropping systems.

Low genetic differentiation across three major ocean populations of the whale shark, Rhincodon typus

BACKGROUND

Whale sharks are a declining species for which little biological data is available. While these animals are protected in many parts of their range, they are fished legally and illegally in some countries. Baseline biological and ecological data are needed to allow the formulation of an effective conservation plan for whale sharks. It is not known, for example, whether the whale shark is represented by a single worldwide panmictic population or by numerous, reproductively isolated populations. Genetic analysis of population structure is one essential component of the baseline data required for whale shark conservation.

METHODOLOGY/PRINCIPAL FINDINGS

We have identified 8 polymorphic microsatellites in the whale shark and used these markers to assess genetic variation and population structure in a panel of whale sharks covering a broad geographic region. This is the first record of microsatellite loci in the whale shark, which displayed an average of 9 alleles per locus and mean H(o) = 0.66 and H(e) = 0.69. All but one of the eight loci meet the expectations of Hardy-Weinberg equilibrium. Analysis of these loci in whale sharks representing three major portions of their range, the Pacific (P), Caribbean (C), and Indian (I) Oceans, determined that there is little population differentiation between animals sampled in different geographic regions, indicating historical gene flow between populations. F(ST) values for inter-ocean comparisons were low (PxC = 0.0387, CxI = 0.0296 and PxI = -0.0022), and only CxI approached statistical significance (p = 0.0495).

CONCLUSIONS/SIGNIFICANCE

We have shown only low levels of genetic differentiation between geographically distinct whale shark populations. Existing satellite tracking data have revealed both regional and long-range migration of whale sharks throughout their range, which supports the finding of gene flow between populations. Whale sharks traverse geographic and political boundaries during their life history and interbreed with animals from distant populations; conservation efforts must therefore target international protection for this species.

Cryptic speciation in a benthic isopod from Patagonian and Falkland Island waters and the impact of glaciations on its population structure

Background

The Falkland Islands and Patagonia are traditionally assigned to the Magellan Biogeographic Province. Most marine species in Falkland waters are also reported from southern Patagonia. It remains unclear if relatively immobile, marine benthic, shallow-water species maintain gene flow, and by what mechanism. Recurrent fluctuations in sea level during glacial cycles are regarded as a possible mechanism that might have allowed genetic exchange between the regions. However, the realized genetic exchange between the Falkland Islands and Patagonia has never been estimated.

Results

This study analyses the genetic structure of three populations of the marine shallow-water isopod Serolis paradoxa (Fabricius, 1775) from the Falkland Islands and southern Patagonia (central Strait of Magellan and the Atlantic opening) applying seven nuclear microsatellites and a fragment of the mitochondrial 16S rRNA gene. Both marker systems report highest genetic diversity for the population from the central Strait of Magellan and lowest for the Falkland Islands. The estimated effective population sizes were large for all populations studied. Significant differentiation was observed among all three populations. The magnitude of differentiation between Patagonia and the Falkland Islands (16S: uncorrected p-distance 2.1%; microsatellites: standardized F'_ST > 0.86) was an order of magnitude higher than between populations from within Patagonia. This indicates that there is currently no effective gene flow for nominal S. paradoxa between these two regions and it has been absent for time exceeding the last glacial maximum. We argue that specimens from the Strait of Magellan and the Falkland Islands very likely represent two distinct species that separated in the mid-Pleistocene (about 1 MY BP).

Conclusion

The results of this study indicate limited gene flow between distant populations of the brooding isopod Serolis paradoxa. The patterns of genetic diversity suggest that the only recently inundated Strait of Magellan was colonized by different source populations, most likely from Atlantic and Pacific coastal waters. Our results demonstrate that more systematic testing of shared faunal inventory and realized genetic exchange between Patagonia and the Falkland Islands is needed before a consensus concerning the position of the Falkland Islands relative to the Magellan zoogeographic province can be reached.

Body temperature predicts maximum microsatellite length in mammals

A long-standing mystery in genome evolution is why short tandem repeats vary so much in length and frequency. Here, we test the hypothesis that body temperature acts to influence the rate and nature of slippage-based mutations. Using the data from both 28 species where genome sequencing is advanced and 76 species from which marker loci have been published, we show that in mammals, maximum repeat number is inversely correlated with body temperature, with warmer-blooded species having shorter 'long' microsatellites. Our results support a model of microsatellite evolution in which maximum length is limited by a temperature-dependent stability threshold.

Amniote phylogenomics: testing evolutionary hypotheses with BAC library scanning and targeted clone analysis of large-scale DNA sequences from reptiles

Phylogenomics research integrating established principles of systematic biology and taking advantage of the wealth of DNA sequences being generated by genome science holds promise for answering long-standing evolutionary questions with orders of magnitude more primary data than in the past. Although it is unrealistic to expect whole-genome initiatives to proceed rapidly for commercially unimportant species such as reptiles, practical approaches utilizing genomic libraries of large-insert clones pave the way for a phylogenomics of species that are nevertheless essential for testing evolutionary hypotheses within a phylogenetic framework. This chapter reviews the case for adopting genome-enabled approaches to evolutionary studies and outlines a program for using bacterial artificial chromosome (BAC) libraries or plasmid libraries as a basis for completing "genome scans" of reptiles. We have used BACs to close a critical gap in the genome database for Reptilia, the sister group of mammals, and present the methodological approaches taken to achieve this as a guideline for designing similar comparative studies. In addition, we provide a detailed step-by-step protocol for BAC-library screening and shotgun sequencing of specific clones containing target genes of evolutionary interest. Taken together, the genome scanning and shotgun sequencing techniques offer complementary diagnostic potential and can substantially increase the scale and power of analyses aimed at testing evolutionary hypotheses for nonmodel species.

Genetic variation and de novo mutations in the parthenogenetic Caucasian rock lizard Darevskia unisexualis

Unisexual all-female lizards of the genus Darevskia that are well adapted to various habitats are known to reproduce normally by true parthenogenesis. Although they consist of unisexual lineages and lack effective genetic recombination, they are characterized by some level of genetic polymorphism. To reveal the mutational contribution to overall genetic variability, the most straightforward and conclusive way is the direct detection of mutation events in pedigree genotyping. Earlier we selected from genomic library of D. unisexualis two polymorphic microsatellite containing loci Du281 and Du215. In this study, these two loci were analyzed to detect possible de novo mutations in 168 parthenogenetic offspring of 49 D. unisexualis mothers and in 147 offspring of 50 D. armeniaca mothers. No mutant alleles were detected in D. armeniaca offspring at both loci, and in D. unisexualis offspring at the Du215 locus. There were a total of seven mutational events in the germ lines of four of the 49 D. unisexualis mothers at the Du281 locus, yielding the mutation rate of 0.1428 events per germ line tissue. Sequencing of the mutant alleles has shown that most mutations occur via deletion or insertion of single microsatellite repeat being identical in all offspring of the family. This indicates that such mutations emerge at the early stages of embryogenesis. In this study we characterized single highly unstable (GATA)(n) containing locus in parthenogenetic lizard species D. unisexualis. Besides, we characterized various types of mutant alleles of this locus found in the D. unisexualis offspring of the first generation. Our data has shown that microsatellite mutations at highly unstable loci can make a significant contribution to population variability of parthenogenetic lizards.

Microsatellite evolution: Mutations, sequence variation, and homoplasy in the hypervariable avian microsatellite locus HrU10

Background

Microsatellites are frequently used genetic markers in a wide range of applications, primarily due to their high length polymorphism levels that can easily be genotyped by fragment length analysis. However, the mode of microsatellite evolution is yet not fully understood, and the role of interrupting motifs for the stability of microsatellites remains to be explored in more detail. Here we present a sequence analysis of mutation events and a description of the structure of repeated regions in the hypervariable, pentanucleotide microsatellite locus HrU10 in barn swallows (Hirundo rustica) and tree swallows (Tachycineta bicolor).

Results

In a large-scale parentage analysis in barn swallows and tree swallows, broods were screened for mutations at the HrU10 locus. In 41 cases in the barn swallows and 15 cases in the tree swallows, mutations corresponding to the loss or gain of one or two repeat units were detected. The parent and mutant offspring alleles were sequenced for 33 of these instances (26 in barn swallows and 7 in tree swallows). Replication slippage was considered the most likely mutational process. We tested the hypothesis that HrU10, a microsatellite with a wide allele size range, has an increased probability of introductions of interruptive motifs (IMs) with increasing length of the repeated region. Indeed, the number and length of the IMs was strongly positively correlated with the total length of the microsatellite. However, there was no significant correlation with the length of the longest stretch of perfectly repeated units, indicating a threshold level for the maximum length of perfectly repeated pentanucleotide motifs in stable HrU10 alleles. The combination of sequence and pedigree data revealed that 15 barn swallow mutations (58%) produced alleles that were size homoplasic to other alleles in the data set.

Conclusion

Our results give further insights into the mode of microsatellite evolution, and support the assumption of increased slippage rate with increased microsatellite length and a stabilizing effect of interrupting motifs for microsatellite regions consisting of perfect repeats. In addition, the observed extent of size homoplasy may impose a general caution against using hypervariable microsatellites in genetic diversity measures when alleles are identified by fragment length analysis only.

An unusually low microsatellite mutation rate in Dictyostelium discoideum, an organism with unusually abundant microsatellites

The genome of the social amoeba Dictyostelium discoideum is known to have a very high density of microsatellite repeats, including thousands of triplet microsatellite repeats in coding regions that apparently code for long runs of single amino acids. We used a mutation accumulation study to see if unusually high microsatellite mutation rates contribute to this pattern. There was a modest bias toward mutations that increase repeat number, but because upward mutations were smaller than downward ones, this did not lead to a net average increase in size. Longer microsatellites had higher mutation rates than shorter ones, but did not show greater directional bias. The most striking finding is that the overall mutation rate is the lowest reported for microsatellites: approximately 1 x 10(-6) for 10 dinucleotide loci and 6 x 10(-6) for 52 trinucleotide loci (which were longer). High microsatellite mutation rates therefore do not explain the high incidence of microsatellites. The causal relation may in fact be reversed, with low mutation rates evolving to protect against deleterious fitness effects of mutation at the numerous microsatellites.

Marker stability throughout 400 days of in vitro hyphal growth in the filamentous ascomycete, Sclerotinia sclerotiorum

The stability of routinely used, population genetic markers through approximately 1 year of continuous laboratory growth was investigated in the common, plant pathogentic ascomycete Sclerotinia sclerotiorum. Given reports of accelerated mutation rates at higher temperatures, both a permissive temperature, 22 degrees C, and a temperature at the high end of tolerance, 30 degrees C, were employed. Because mycelial growth rate was tracked among mitotic lineages established for each strain, a subsidiary objective was addressed, testing the stability of a 30 degrees C-competent phenotype. Twelve laboratory strains of S. sclerotiorum, including the genome sequence isolate, 1980, were propagated serially for up to 400 days at 22 degrees C. Five of these strains were also propagated at 30 degrees C. No mutations were observed in mycelial compatibility groupings (MCGs), DNA fingerprints, alleles at 7 microsatellite loci, or alleles at 56 AFLP loci. All of these markers show variation in field populations, which are likely much larger and influenced by different and more stochastic environmental processes. In S. sclerotiorum, population genetic markers were stable over time through serial transfer and growth of laboratory strains at both 22 degrees C and 30 degrees C. The strain isolated after extended drought and capable of infecting plants at 28 degrees C demonstrated the stability of its high temperature-competent phenotype, in addition to its stable growth rate at 22 degrees C. This observation has implications for modeling pathogen tolerance or adaptation under conditions of environmental stochasticity, including climate warming.

Dynamic mutations as digital genetic modulators of brain development, function and dysfunction

A substantial portion of the human genome has been found to consist of simple sequence repeats, including microsatellites and minisatellites. Microsatellites, tandem repeats of 1-6 nucleotides, form the template for dynamic mutations, which involve heritable changes in the lengths of repeat sequences. In recent years, a large number of human disorders have been found to be caused by dynamic mutations, the most common of which are trinucleotide repeat expansion diseases. Dynamic mutations are common to numerous nervous system disorders, including Huntington's disease, various spinocerebellar ataxias, fragile X syndrome, fragile X tremor/ataxia syndrome, Friedreich ataxia and other neurodegenerative disorders. The involvement of dynamic mutations in brain disorders will be reviewed, with a focus on the large group caused by CAG/glutamine repeat expansions. We will also outline a proposed role of tandem repeat polymorphisms (TRPs), with unique 'digital' genetic distributions, in modulating brain development and normal function, so as to generate additional mutational diversity upon which natural selection may act.

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.

Evolution of hypervariable microsatellites in apomictic polyploid lineages of Ranunculus carpaticola: directional bias at dinucleotide loci

Microsatellites are widely used in genetic and evolutionary analyses, but their own evolution is far from simple. The mechanisms maintaining the mutational patterns of simple repeats and the typical stable allele-frequency distributions are still poorly understood. Asexual lineages may provide particularly informative models for the indirect study of microsatellite evolution, because their genomes act as complete linkage groups, with mutations being the only source of genetic variation. Here, we study the direction of accumulated dinucleotide microsatellite mutations in wild asexual lineages of hexaploid Ranunculus carpaticola. Whereas the overall number of contractions is not significantly different from that of expansions, the within-locus frequency of contractions, but not of expansions, significantly increases with allele length. Moreover, within-locus polymorphism is positively correlated with allele length, but this relationship is due solely to the influence of contraction mutations. Such asymmetries may explain length constraints generally observed with microsatellites and are consistent with stable, bell-shaped allele-frequency distributions. Although apomictic and allohexaploid, the R. carpaticola lineages show mutational patterns resembling the trends observed in a broad range of organisms, including sexuals and diploids, suggesting that, even if not of germline origin, the mutations in these apomicts may be the consequence of similar mechanisms.

Mutation rates at Y chromosome specific microsatellites

A collaborative work was carried out by the Spanish and Portuguese ISFG Working Group (GEP-ISFG) to estimate Y-STR mutation rates. Seventeen Y chromosome STR loci (DYS19, DYS385, DYS389I and II, DYS390, DYS391, DYS392, DYS393, DYS437, DYS438, DYS439, DYS460, DYS461, DYS635 [GATA C4], GATA H4, and GATA A10) were analyzed in a sample of 3,026 father/son pairs. Among 27,029 allele transfers, 54 mutations were observed, with an overall mutation rate across the 17 loci of 1.998 x 10(-3) (95% CI, 1.501 x 10(-3) to 2.606 x 10(-3)). With just one exception, all of the mutations were single-step, and they were observed only once per gametogenesis. Repeat gains were more frequent than losses, longer alleles were found to be more mutable, and the mutation rate seemed to increase with the father's age. Hum Mutat 26(6), 520-528, 2005. (c) 2005 Wiley-Liss, Inc.

The M53I mutation inCDKN2A is a founder mutation that predominates in melanoma patients with Scottish ancestry

Germline mutations in the tumor suppressor gene CDKN2A have been shown to predispose to cutaneous malignant melanoma. The M53I mutation is the most common CDKN2A mutation identified in Scottish melanoma patients and is also found in a small number of families in other countries. The aim of this study was to determine whether the occurrence of this mutation is due to a common ancestor originating from Scotland, and if so, to estimate how long ago the mutation arose. We examined 18 families carrying the M53I mutation: six from Scotland, five from Canada, four from Australia, and three from America. Haplotypes derived from segregation of seven informative microsatellite markers flanking CDKN2A were constructed in each family. Our findings show that 14 of 18 families carry a common ancestral haplotype on which the mutation arose approximately 88 generations ago (1-LOD-unit support interval 44-198 generations). This haplotype is very rare in controls, which supports the idea that it is a common founder mutation haplotype. The four M53I families that do not share the consensus haplotype may in fact have arisen from the same founder, but this is potentially obscured by presumed replication slippage for some of the microsatellite markers tested.

Inference on microsatellite mutation processes in the invasive mite, Varroa destructor, using reversible jump Markov chain Monte Carlo

Varroa destructor is a parasitic mite of the Eastern honeybee Apis cerana. Fifty years ago, two distinct evolutionary lineages (Korean and Japanese) invaded the Western honeybee Apis mellifera. This haplo-diploid parasite species reproduces mainly through brother-sister matings, a system which largely favors the fixation of new mutations. In a worldwide sample of 225 individuals from 21 locations collected on Western honeybees and analyzed at 19 microsatellite loci, a series of de novo mutations was observed. Using historical data concerning the invasion, this original biological system has been exploited to compare three mutation models with allele size constraints for microsatellite markers: stepwise (SMM) and generalized (GSM) mutation models, and a model with mutation rate increasing exponentially with microsatellite length (ESM). Posterior probabilities of the three models have been estimated for each locus individually using reversible jump Markov Chain Monte Carlo. The relative support of each model varies widely among loci, but the GSM is the only model that always receives at least 9% support, whatever the locus. The analysis also provides robust estimates of mutation parameters for each locus and of the divergence time of the two invasive lineages (67,000 generations with a 90% credibility interval of 35,000-174,000). With an average of 10 generations per year, this divergence time fits with the last post-glacial Korea-Japan land separation.

Microsatellite evolution in the mitochondrial genome of Bechstein's bat (Myotis bechsteinii)

Being highly polymorphic, microsatellites are widely used genetic markers. They are abundant throughout the nuclear genomes of eukaryotes but rare in the mitochondrial genomes (mtDNA) of animals. We describe a short but highly polymorphic AT microsatellite in the mtDNA control region of Bechstein's bat and discuss the role of mutation, genetic drift, and selection in maintaining its variability. As heteroplasmy and hence mutation rate were positively correlated with repeat number, a simple mutation model cannot explain the observed frequency distribution of AT copy numbers. Because of the unimodal distribution of repeat numbers found in heteroplasmic individuals, single step mutations are likely to be the predominant mechanism of copy number alternations. Above a certain copy number (seven repeats), deletions of single dinucleotide repeats seem to be more common than additions, which results in a decrease in frequency of long alleles. Heteroplasmy was inherited from mothers to their offspring and no evidence of paternal inheritance of mitochondria was found. Genetic differences accumulated with more distant ancestry, which suggests that microsatellites can be useful genetic markers in population genetics.

Standardizing for microsatellite length in comparisons of genetic diversity

Mutation rates at microsatellites tend to increase with the number of repeats of the motif, leading to higher levels of polymorphism at long microsatellites. To standardize levels of diversity when microsatellites differ in size, we investigate the relationship between tract length and variation and provide a formula to adjust allelic richness to a fixed mean number of repeats in the specific case of chloroplast microsatellites. A comparison between 39 loci from eight species of conifers (where chloroplast DNA is paternally inherited) and 64 loci from 12 species of angiosperms (where chloroplast DNA is generally predominantly maternally inherited) indicates that the greater allelic richness found in conifers remains significant after controlling for number of repeats. The approach stresses the advantage of reporting variation in number of repeats instead of relative fragment sizes.