Multiple levels of single-strand slippage at cetacean tri- and tetranucleotide repeat microsatellite loci

STRs: Ancient Architectures of the Genome beyond the Sequence

Short tandem repeats (STRs) are commonly defined as short runs of repetitive nucleotides, consisting of tandemly repeating 2-6- bp motif units, which are ubiquitously distributed throughout genomes. Functional STRs are polymorphic in the population, and their variations influence gene expression, which subsequently may result in pathogenic phenotypes. To understand STR phenotypic effects and their functional roles, we describe four different mutational mechanisms including the unequal crossing-over model, gene conversion, retrotransposition mechanism and replication slippage. Due to the multi-allelic nature, small length, abundance, high variability, codominant inheritance, nearly neutral evolution, extensive genome coverage and simple assaying of STRs, these markers are widely used in various types of biological research, including population genetics studies, genome mapping, molecular epidemiology, paternity analysis and gene flow studies. In this review, we focus on the current knowledge regarding STR genomic distribution, function, mutation and applications.

Long-term isolation at a low effective population size greatly reduced genetic diversity in Gulf of California fin whales

The Gulf of California, Mexico is home to many cetacean species, including a presumed resident population of fin whales, Balaenoptera physalus. Past studies reported very low levels of genetic diversity among Gulf of California fin whales and a significant level of genetic differentiation from con-specifics in the eastern North Pacific. The aim of the present study was to assess the degree and timing of the isolation of Gulf of California fin whales in a population genetic analysis of 18 nuclear microsatellite genotypes from 402 samples and 565 mitochondrial control region DNA sequences (including mitochondrial sequences retrieved from NCBI). The analyses revealed that the Gulf of California fin whale population was founded ~2.3 thousand years ago and has since remained at a low effective population size (~360) and isolated from the eastern North Pacific (Nem between 0.89-1.4). The low effective population size and high degree of isolation implied that Gulf of California fin whales are vulnerable to the negative effects of genetic drift, human-caused mortality and habitat change.

Assessing Genetic Structure in Common but Ecologically Distinct Carnivores: The Stone Marten and Red Fox

The identification of populations and spatial genetic patterns is important for ecological and conservation research, and spatially explicit individual-based methods have been recognised as powerful tools in this context. Mammalian carnivores are intrinsically vulnerable to habitat fragmentation but not much is known about the genetic consequences of fragmentation in common species. Stone martens (Martes foina) and red foxes (Vulpes vulpes) share a widespread Palearctic distribution and are considered habitat generalists, but in the Iberian Peninsula stone martens tend to occur in higher quality habitats. We compared their genetic structure in Portugal to see if they are consistent with their differences in ecological plasticity, and also to illustrate an approach to explicitly delineate the spatial boundaries of consistently identified genetic units. We analysed microsatellite data using spatial Bayesian clustering methods (implemented in the software BAPS, GENELAND and TESS), a progressive partitioning approach and a multivariate technique (Spatial Principal Components Analysis-sPCA). Three consensus Bayesian clusters were identified for the stone marten. No consensus was achieved for the red fox, but one cluster was the most probable clustering solution. Progressive partitioning and sPCA suggested additional clusters in the stone marten but they were not consistent among methods and were geographically incoherent. The contrasting results between the two species are consistent with the literature reporting stricter ecological requirements of the stone marten in the Iberian Peninsula. The observed genetic structure in the stone marten may have been influenced by landscape features, particularly rivers, and fragmentation. We suggest that an approach based on a consensus clustering solution of multiple different algorithms may provide an objective and effective means to delineate potential boundaries of inferred subpopulations. sPCA and progressive partitioning offer further verification of possible population structure and may be useful for revealing cryptic spatial genetic patterns worth further investigation.

Fine-scale population structure of blue whale wintering aggregations in the Gulf of California

Population differentiation in environments without well-defined geographical barriers represents a challenge for wildlife management. Based on a comprehensive database of individual sighting records (1988–2009) of blue whales from the winter/calving Gulf of California, we assessed the fine-scale genetic and spatial structure of the population using individual-based approaches. Skin samples of 187 individuals were analyzed for nine microsatellite loci. A single population with no divergence among years and months and no isolation by distance (Rxy = 0.1–0.001, p>0.05) were found. We ran two Bayesian clustering methods using Structure and Geneland softwares in two different ways: 1) a general analysis including all individuals in which a single cluster was identified with both softwares; 2) a specific analysis of females only in which two main clusters (Loreto Bay and northern areas, and San Jose-La Paz Bay area) were revealed by Geneland program. This study provides information indicating that blue whales wintering in the Gulf of California are part of a single population unit and showed a fine-scale structure among females, possibly associated with their high site fidelity, particularly when attending calves. It is likely that the loss of genetic variation is minimized by male mediated gene flow, which may reduce the genetic drift effect. Opportunities for kin selection may also influence calf survival and, in consequence, have a positive impact on population demography in this small and endangered population.

Development of chickpea EST-SSR markers and analysis of allelic variation across related species

Despite chickpea being the third important grain legume, there is a limited availability of genomic resources, especially of the expressed sequence tag (EST)-based markers. In this study, we generated 822 chickpea ESTs from immature seeds as well as exploited 1,309 ESTs from the chickpea database, thus utilizing a total of 2,131 EST sequences for development of functional EST-SSR markers. Two hundred and forty-six simple sequence repeat (SSR) motifs were identified from which 183 primer pairs were designed and 60 validated as functional markers. Genetic diversity analysis across 30 chickpea accessions revealed ten markers to be polymorphic producing a total of 29 alleles and an observed heterozygosity average of 0.16 thereby exhibiting low levels of intra-specific polymorphism. However, the markers exhibited high cross-species transferability ranging from 68.3 to 96.6% across the six annual Cicer species and from 29.4 to 61.7% across the seven legume genera. Sequence analysis of size variant amplicons from various species revealed that size polymorphism was due to multiple events such as copy number variation, point mutations and insertions/deletions in the microsatellite repeat as well as in the flanking regions. Interestingly, a wide prevalence of crossability-group-specific sequence variations were observed among Cicer species that were phylogenetically informative. The neighbor joining dendrogram clearly separated the chickpea cultivars from the wild Cicer and validated the proximity of C. judaicum with C. pinnatifidum. Hence, this study for the first time provides an insight into the distribution of SSRs in the chickpea transcribed regions and also demonstrates the development and utilization of genic-SSRs. In addition to proving their suitability for genetic diversity analysis, their high rates of transferability also proved their potential for comparative genomic studies and for following gene introgressions and evolution in wild species, which constitute the valuable secondary genepool in chickpea.

New polymorphic tetranucleotide microsatellites improve scoring accuracy in the bottlenose dolphin Tursiops aduncus

We isolated and characterized 19 novel tetranucleotide microsatellite markers in the Indo-Pacific bottlenose dolphin (Tursiops aduncus) in order to improve genotyping accuracy in applications like large-scale population-wide paternity and relatedness assessments. One hundred T. aduncus from Shark Bay, Western Australia, were screened for polymorphism. Cross-amplification was tested on four other small odontocete species. The new tetranucleotide microsatellite loci showed a more than fourfold higher scoring accuracy and significantly lower stutter formation compared to eight dinucleotide loci, although overall allelic diversity was significantly reduced.

Sequence-dependent effect of interruptions on microsatellite mutation rate in mismatch repair-deficient human cells

Although microsatellite mutation rates generally increase with increasing length of the repeat tract, interruptions in a microsatellite may stabilize it. We have performed a direct analysis of the effect of microsatellite interruptions on mutation rate and spectrum in cultured mammalian cells. Two mononucleotide sequences (G(17) and A(17)) and a dinucleotide [(CA)(17)] were compared with interrupted repeats of the same size and with sequences of 8 repeat units. MMR-deficient (MMR(-)) cells were used for these studies to eliminate effects of this repair process. Mutation rates were determined by fluctuation analysis on cells containing a microsatellite sequence at the 5' end of an antibiotic-resistance gene; the vector carrying this sequence was integrated in the genome of the cells. In general, interrupted sequences had lower mutation rates than perfect ones of the same size, but the magnitude of the difference was dependent upon the sequence of the interrupting base(s). Some interrupted repeats had mutation rates that were lower than those of perfect sequences of the same length but similar to those of half the length. This suggests that interrupting bases effectively divide microsatellites into smaller repeat runs with mutational characteristics different from those of the corresponding full-length microsatellite. We conclude that interruptions decrease microsatellite mutation rate and influence the spectrum of frameshift mutations. The sequence of the interrupting base(s) determines the magnitude of the effect on mutation rate.

Contrasting patterns of mitochondrial and microsatellite genetic structure among Western European populations of tawny owls (Strix aluco)

A recent study of mitochondrial phylogeography of tawny owls (Strix aluco) in western Europe suggested that this species survived the Pleistocene glaciations in three allopatric refugia located in Iberia, Italy, and the Balkans, and the latter was likely the predominant source of postglacial colonization of northern Europe. New data from seven microsatellite loci from 184 individual owls distributed among 14 populations were used to assess the genetic congruence between nuclear and mitochondrial DNA (mtDNA) markers. Microsatellites corroborated the major phylogeographical conclusions reached on the basis of the mtDNA sequences, but also showed important differences leading to novel inferences. Microsatellites corroborated the three major refugia and supported the Balkan origin of northern populations. When corrected for differences in effective population size, microsatellites and mtDNA yielded generally congruent overall estimates of population structure (N*ST=0.12 vs. RST=0.16); however, there was substantial heterogeneity in the RST among the seven nuclear loci that was not correlated with heterozygosity. Populations representing the Balkans postglacial expansion interact with populations from the other two refugia forming two clines near the Alps and the Pyrenees. In both cases, the apparent position of the contact zones differed substantially between markers due to the genetic composition of populations sampled in northern Italy and Madrid. Microsatellite data did not corroborate the lower genetic diversity of northern, recently populated regions as was found with mtDNA; this discrepancy was taken as evidence for a recent bottleneck recovery. Finally, this study suggests that congruence among genetic markers should be more likely in cases of range expansion into new areas than when populations interact across contact zones.

Molecular characterization of allelic variants of (GATA)n microsatellite loci in parthenogenetic lizards Darevskia unisexualis (Lacertidae)

Populations of parthenogenetic lizards of the genus Darevskia consist of genetically identical animals, and represent a unique model for studying the molecular mechanisms underlying the variability and evolution of hypervariable DNA repeats. As unisexual lineages, parthenogenetic lizards are characterized by some level of genetic diversity at microsatellite loci. We cloned and sequenced a number of (GATA)n microsatellite loci of Darevskia unisexualis. PCR products from these loci were also sequenced and the degree of intraspecific polymorphism was assessed. Among the five (GATA)n loci analysed, two (Du215 and Du281) were polymorphic. Cross-species analysis of Du215 and Du281 indicate that the priming sites at the D. unisexualis loci are conserved in the bisexual parental species, D. raddei and D. valentini. Sequencing the PCR products amplified from Du215 and Du281 and from monomorphic Du323 showed that allelic differences at the polymorphic loci are caused by microsatellite mutations and by point mutations in the flanking regions. The haplotypes identified among the allelic variants of Du281 and among its orthologues in the parental species provide new evidence of the cross-species origin of D. unisexualis. To our knowledge, these data are the first to characterize the nucleotide sequences of allelic variants at microsatellite loci within parthenogenetic vertebrate animals.

Mutation patterns at dinucleotide microsatellite loci in humans

Microsatellites are a major type of molecular markers in genetics studies. Their mutational dynamics are not clear. We investigated the patterns and characteristics of 97 mutation events unambiguously identified, from 53 multigenerational pedigrees with 630 subjects, at 362 autosomal dinucleotide microsatellite loci. A size-dependent mutation bias (in which long alleles are biased toward contraction, whereas short alleles are biased toward expansion) is observed. There is a statistically significant negative relationship between the magnitude (repeat numbers changed during mutation) and direction (contraction or expansion) of mutations and standardized allele size. Contrasting with earlier findings in humans, most mutation events (63%) in our study are multistep events that involve changes of more than one repeat unit. There was no correlation between mutation rate and recombination rate. Our data indicate that mutational dynamics at microsatellite loci are more complicated than the generalized stepwise mutation models.

Dinucleotide microsatellites from Eucalyptus sieberi: inheritance, diversity, and improved scoring of single-base differences

Eight dinucleotide microsatellites were developed in Eucalyptus sieberi L. Johnson (silvertop ash), a member of the subgenus Eucalyptus. Transfer of six of these to the subgenus Symphyomyrtus and their Mendelian inheritance are demonstrated using a full-sib cross in Eucalyptus nitens. Genetic diversity parameters are presented for the eight loci based on a sample of 100 old-growth E. sieberi trees from a single natural stand. One locus, Es266, had an atypically high fixation index, and significantly deviated from Hardy-Weinberg equilibrium genotypic proportions, indicating the likely presence of null alleles. Two of the loci, Es076 and Es140, had many alleles that differed in size by only a single base pair, possibly because of short poly(A) or poly(T) stretches in their flanking regions. These two loci were by far the most polymorphic, but were difficult to score reliably on a capillary DNA sequencer. Reliability of scoring of these two one-base microsatellite loci was markedly improved by the incorporation of internal reference alleles into each sample analysed.Key words: SSRs, single base pair alleles, null alleles, internal reference alleles.

Statistical approaches to paternity analysis in natural populations and applications to the North Atlantic humpback whale

We present a new method for paternity analysis in natural populations that is based on genotypic data that can take the sampling fraction of putative parents into account. The method allows paternity assignment to be performed in a decision theoretic framework. Simulations are performed to evaluate the utility and robustness of the method and to assess how many loci are necessary for reliable paternity inference. In addition we present a method for testing hypotheses regarding relative reproductive success of different ecologically or behaviorally defined groups as well as a new method for estimating the current population size of males from genotypic data. This method is an extension of the fractional paternity method to the case where only a proportion of all putative fathers have been sampled. It can also be applied to provide abundance estimates of the number of breeding males from genetic data. Throughout, the methods were applied to genotypic data collected from North Atlantic humpback whales (Megaptera novaeangliae) to test if the males that appear dominant during the mating season have a higher reproductive success than the subdominant males.

Microsatellite mutations in the germline: implications for evolutionary inference

Microsatellite DNA sequences mutate at rates several orders of magnitude higher than that of the bulk of DNA. Such high rates mean that spontaneous mutations that form new-length variants can realistically be seen in pedigree analysis. Data on observed mutation events from various organisms are now accumulating, allowing inferences on DNA sequence evolution to be made through an unusually direct approach. Here I discuss and integrate microsatellite mutation data in an evolutionary context. A striking feature of the mutation process is that it seems highly heterogeneous, with distinct differences between species, repeat types, loci and alleles. Age and sex also affect the mutation rate. Within genomes at equilibrium, the microsatellite-length distribution is a delicate balance between biased mutation processes and point mutations acting towards the decay of repetitive DNA. Indeed, simple repeats do not evolve simply.

Complex mutations in a high proportion of microsatellite loci from the protozoan parasite Plasmodium falciparum

Microsatellite loci are generally assumed to evolve via a stepwise mutational process and a battery of statistical techniques has been developed in recent years based on this or related mutation models. It is therefore important to investigate the appropriateness of these models in a wide variety of taxa. We used two approaches to examine mutation patterns in the malaria parasite Plasmodium falciparum: (i) we examined sequence variation at 12 tri-nucleotide repeat loci; and (ii) we analysed patterns of repeat structure and heterozygosity at 114 loci using data from 12 laboratory parasite lines. The sequencing study revealed complex patterns of mutation in five of the 12 loci studied. Alleles at two loci contain indels of 24 bp and 57 bp in flanking regions, while in the other three loci, blocks of imperfect microsatellites appear to be duplicated or inserted; these loci essentially consist of minisatellite repeats, with each repeat unit containing four to eight microsatellites. The survey of heterozygosity revealed a positive relationship between repeat number and microsatellite variability for both di- and trinucleotides, indicating a higher mutation rate in loci with longer repeat arrays. Comparisons of levels of variation in different repeat types indicate that the mutation rate of dinucleotide-bearing loci is 1.6-2.1 times faster than trinucleotides, consistent with the lower mean number of repeats in trinucleotide-bearing loci. However, despite the evidence that microsatellite arrays themselves are evolving in a manner consistent with stepwise mutation model in P. falciparum, the high frequency of complex mutations precludes the use of analytical tools based on this mutation model for many microsatellite-bearing loci in this protozoan. The results call into question the generality of models based on stepwise mutation for analysing microsatellite data, but also demonstrate the ease with which loci that violate model assumptions can be detected using minimal sequencing effort.

Removal of microsatellite interruptions by DNA replication slippage: phylogenetic evidence from Drosophila

Microsatellites are tandem repetitions of short (1-6 bp) motifs. It is widely assumed that microsatellites degenerate through the accumulation of base substitutions in the repeat array. Using a phylogenetic framework, we studied the evolutionary dynamics of interruptions in three Drosophila microsatellite loci. For all three loci, we show that the interruptions in a microsatellite can be lost, resulting in a longer uninterrupted microsatellite stretch. These results indicate that mutations in the microsatellite array do not necessarily lead to decay but may represent only a transition state during the evolution of a microsatellite. Most likely, this purification of interrupted microsatellites is caused by DNA replication slippage.