Rates and patterns of microsatellite mutations in pink salmon

A 13-plex of tetra- and penta-STRs to identify zebrafish

The zebrafish species Danio rerio has become one of the major vertebrate model organisms used in biomedical research. However, there are aspects of the model that need to be improved. One of these is the ability to identify individual fish and fish lines by DNA profiling. Although many dinucleotide short tandem repeat (diSTR) markers are available for this and similar purposes, they have certain disadvantages such as an excessive polymerase slippage ("stutter") that causes difficulties in automated genotyping and cross-laboratory comparisons. Here we report on the development of a 13-plex of tetranucleotide and pentanucleotide STRs (tetraSTRs and pentaSTRs, respectively) that have low stutter. The system uses an inexpensive universal primer labelling system, which can easily be converted to a direct labeling system if desired. This 13-plex was examined in three zebrafish lines (NHGRI-1, kca33Tg, and kca66Tg, originally obtained from ZIRC). The average observed heterozygosity (Ho) and expected heterozygosity (He) in these highly inbred lines were 0.291 and 0.359, respectively, which is very similar to what has been found with diSTRs. The probability of identity (PI) for all fish tested was 2.1 × 10^-5 and the PI for siblings (PIsib) was 6.4 × 10^-3, as calculated by the Genalex package. Ninety percent of the fish tested were correctly identified with their respective strains. It is also demonstrated that this panel can be used to confirm doubled-haploid cell lines. This multiplex should find multiple uses for improving the accuracy and reproducibility of studies using the zebrafish model.

Phylogeography and Population Genetic Analyses in the Iberian Toothcarp (Aphanius iberus Valenciennes, 1846) at Different Time Scales

Secondary freshwater fish species inhabiting fluctuating and extreme environments are susceptible to changes in dispersion, effective population size, and genetic structure. The Iberian toothcarp Aphanius iberus is an endemic cyprinodontid of the Iberian Peninsula restricted to brackish water of salt marshes and coastal lagoons on the eastern Spanish Mediterranean coast. In this study, we analyzed mitochondrial cytochrome b (cyt b) DNA and microsatellite variation to evaluate ways in which the processes of extinction, dispersal, and colonization of A. iberus across its geographic distribution have affected its population genetic structure over time and space. The A. iberus network reconstruction indicated subtle levels of phylogeographic structuring. This, combined with substantial mitochondrial DNA (mtDNA) genetic diversity, suggests that Pleistocene glaciations had a lesser effect on the demographic structure of its populations than was the case for Iberian freshwater species with a similar distribution. Haplotype network, hierarchical analysis of molecular variance, and pairwise ΦST comparisons involving some Levantine samples showed a relatively high degree of mtDNA differentiation, which could be explained by historical isolation of the Villena Lagoon population. Conversely, significant genetic differentiation that follows an isolation-by-distance pattern, and a reduction in Ne though time was detected with microsatellites, suggesting extensive habitat fragmentation on the Mediterranean coast of the Iberian Peninsula over the past hundreds of years. At a smaller geographical scale (Mar Menor Lagoon), habitat fragmentation, probably due to human activity, appears to have resulted in substantially reduced migration and increased genetic drift, as shown by expanded genetic differentiation of populations.

Microsatellite Mutation Rate in Atlantic Sturgeon (Acipenser oxyrinchus)

Understanding mutation rates can greatly extend the utility of population and conservation genetic analyses. Herein, we present an estimate of genome-wide microsatellite mutation rate in Atlantic sturgeon (Acipenser oxyrinchus) based on parent-offspring transmission patterns. We screened 307 individuals for parentage and mutation-rate analysis applying 43 variable markers. Out of 13228 allele transfers, 11 mutations were detected, producing a mutation rate of 8.3 × 10-4 per locus per generation (95% confidence interval: 1.48 × 10-3, 4.15 × 10-4). Single-step mutations predominated and there were trends toward mutations in loci with greater polymorphism and allele length. Two of the detected mutations were most probably cluster mutations, being identified in 12 and 28 sibs, respectively. Finally, we observed evidences of polyploidy based on the sporadic presence of 3 or 4 alleles per locus in the genotyped individuals, supporting previous reports of incomplete diploidization in Atlantic sturgeon.

Reconstructing the demographic history of divergence between European river and brook lampreys using approximate Bayesian computations

Inferring the history of isolation and gene flow during species divergence is a central question in evolutionary biology. The European river lamprey (Lampetra fluviatilis) and brook lamprey (L. planeri) show a low reproductive isolation but have highly distinct life histories, the former being parasitic-anadromous and the latter non-parasitic and freshwater resident. Here we used microsatellite data from six replicated population pairs to reconstruct their history of divergence using an approximate Bayesian computation framework combined with a random forest model. In most population pairs, scenarios of divergence with recent isolation were outcompeted by scenarios proposing ongoing gene flow, namely the Secondary Contact (SC) and Isolation with Migration (IM) models. The estimation of demographic parameters under the SC model indicated a time of secondary contact close to the time of speciation, explaining why SC and IM models could not be discriminated. In case of an ancient secondary contact, the historical signal of divergence is lost and neutral markers converge to the same equilibrium as under the less parameterized model allowing ongoing gene flow. Our results imply that models of secondary contacts should be systematically compared to models of divergence with gene flow; given the difficulty to discriminate among these models, we suggest that genome-wide data are needed to adequately reconstruct divergence history.

Microsatellite evolutionary rate and pattern in Schistocerca gregaria inferred from direct observation of germline mutations

Unravelling variation among taxonomic orders regarding the rate of evolution in microsatellites is crucial for evolutionary biology and population genetics research. The mean mutation rate of microsatellites tends to be lower in arthropods than in vertebrates, but data are scarce and mostly concern accumulation of mutations in model species. Based on parent-offspring segregations and a hierarchical Bayesian model, the mean rate of mutation in the orthopteran insect Schistocerca gregaria was estimated at 2.1e(-4) per generation per untranscribed dinucleotide locus. This is close to vertebrate estimates and one order of magnitude higher than estimates from species of other arthropod orders, such as Drosophila melanogaster and Daphnia pulex. We also found evidence of a directional bias towards expansions even for long alleles and exceptionally large ranges of allele sizes. Finally, at transcribed microsatellites, the mean rate of mutation was half the rate found at untranscribed loci and the mutational model deviated from that usually considered, with most mutations involving multistep changes that avoid disrupting the reading frame. Our direct estimates of mutation rate were discussed in the light of peculiar biological and genomic features of S. gregaria, including specificities in mismatch repair and the dependence of its activity to allele length. Shedding new light on the mutational dynamics of grasshopper microsatellites is of critical importance for a number of research fields. As an illustration, we showed how our findings improve microsatellite application in population genetics, by obtaining a more precise estimation of S. gregaria effective population size from a published data set based on the same microsatellites.

Fitness of hatchery-reared salmonids in the wild

Accumulating data indicate that hatchery fish have lower fitness in natural environments than wild fish. This fitness decline can occur very quickly, sometimes following only one or two generations of captive rearing. In this review, we summarize existing data on the fitness of hatchery fish in the wild, and we investigate the conditions under which rapid fitness declines can occur. The summary of studies to date suggests: nonlocal hatchery stocks consistently reproduce very poorly in the wild; hatchery stocks that use wild, local fish for captive propagation generally perform better than nonlocal stocks, but often worse than wild fish. However, the data above are from a limited number of studies and species, and more studies are needed before one can generalize further. We used a simple quantitative genetic model to evaluate whether domestication selection is a sufficient explanation for some observed rapid fitness declines. We show that if selection acts on a single trait, such rapid effects can be explained only when selection is very strong, both in captivity and in the wild, and when the heritability of the trait under selection is high. If selection acts on multiple traits throughout the life cycle, rapid fitness declines are plausible.

Density drives polyandry and relatedness influences paternal success in the Pacific gooseneck barnacle, Pollicipes elegans

BACKGROUND

Polyandry is a common mating strategy in animals, increasing female fitness through direct (material) and indirect (genetic) benefits. Most theories about the benefits of polyandry come from studies of terrestrial animals, which have relatively complex mating systems and behaviors; less is known about the potential benefits of polyandry in sessile marine animals, for which potential mates may be scarce and females have less control over pre-copulatory mate choice. Here, we used microsatellite markers to examine multiple paternity in natural aggregations of the Pacific gooseneck barnacle Pollicipes elegans, testing the effect of density on paternity and mate relatedness on male reproductive success.

RESULTS

We found that multiple paternity was very common (79% of broods), with up to five fathers contributing to a brood, though power was relatively low to detect more than four fathers. Density had a significant and positive linear effect on the number of fathers siring a brood, though this relationship leveled off at high numbers of fathers, which may reflect a lack of power and/or an upper limit to polyandry in this species. Significant skew in male reproductive contribution in multiply-sired broods was observed and we found a positive and significant relationship between the proportion of offspring sired and the genetic similarity between mates, suggesting that genetic compatibility may influence reproductive success in this species.

CONCLUSIONS

To our knowledge, this is the first study to show high levels of multiple paternity in a barnacle, and overall, patterns of paternity in P. elegans appear to be driven primarily by mate availability. Evidence of paternity bias for males with higher relatedness suggests some form of post-copulatory sexual selection is taking place, but more work is needed to determine whether it operates during or post-fertilization. Overall, our results suggest that while polyandry in P. elegans is driven by mate availability, it may also provide a mechanism for females to ensure fertilization by compatible gametes and increase reproductive success in this sessile species.

Evaluation of rockfish conservation area networks in the United States and Canada relative to the dispersal distance for black rockfish (Sebastes melanops)

Marine reserves networks are implemented as a way to mitigate the impact of fishing on marine ecosystems. Theory suggests that a reserve network will function synergistically when connected by dispersal, but the scale of dispersal is often unknown. On the Pacific coast of the United States and Canada, both countries have recently implemented a number of rockfish conservation areas (RCAs) to protect exploited rockfish species, but no study has evaluated the connectivity within networks in each country or between the two countries. We used isolation-by-distance theory to estimate the scale of dispersal from microsatellite data in the black rockfish, Sebastes melanops, and compared this estimate with the distance between RCAs that would protect this species. Within each country, we found that the distance between RCAs was generally within the confidence intervals of mean dispersal per generation. The distance between these two RCA networks, however, was greater than the average dispersal per generation. The data were also consistent with a genetic break between southern Oregon and central Oregon. We discuss whether additional nearshore RCAs in southern Oregon and Washington would help promote connectivity between RCA's for shallow-water rockfishes.

Characterization of unstable microsatellites in mice: no evidence for germline mutation induction following gamma-radiation exposure

Large tandem repeat DNA loci such as expanded simple tandem repeats and minisatellites are efficient markers for detecting germline mutations; however, mutation detection using these loci can be imprecise and difficult to standardize across labs. Short-tandem repeats, such as microsatellites, offer more precise and high-throughput mutation detection, but germline mutation induction at these loci has not yet been studied in model organisms such as mice. In this study, we used microsatellite enrichment and large-scale DNA sequencing of several closely related inbred mouse lines to identify a panel of 19 polymorphic microsatellites with potentially high spontaneous mutation frequencies. We used this panel and four additional loci from other sources to quantify spontaneous mutation frequency in pedigrees of outbred Swiss-Webster mice. In addition, we also examined mutation induction in families in which sires were treated with acute doses of either 0.5 Gy or 1.0 Gy gamma-irradiation to spermatogonial stem cells. Per locus mutation frequencies ranged from 0 to 5.03 × 10(-3). Considering only the 11 loci with mutations, the mutation frequencies were: control 2.78 × 10(-3), 0.5 Gy 4.09 × 10(-3), and 1.0 Gy 1.82 × 10(-3). There were no statistically significant changes in mutation frequencies among treatment groups. Our study provides the first direct quantification of microsatellite mutation frequency in the mouse germline, but shows no evidence for mutation induction at pre-meiotic male germ cells following acute gamma-irradiation. Further work using the panel is needed to examine mutation induction at different doses of radiation, exposure durations, and stages during spermatogenesis.

Reliability of genetic bottleneck tests for detecting recent population declines

The identification of population bottlenecks is critical in conservation because populations that have experienced significant reductions in abundance are subject to a variety of genetic and demographic processes that can hasten extinction. Genetic bottleneck tests constitute an appealing and popular approach for determining if a population decline has occurred because they only require sampling at a single point in time, yet reflect demographic history over multiple generations. However, a review of the published literature indicates that, as typically applied, microsatellite-based bottleneck tests often do not detect bottlenecks in vertebrate populations known to have experienced declines. This observation was supported by simulations that revealed that bottleneck tests can have limited statistical power to detect bottlenecks largely as a result of limited sample sizes typically used in published studies. Moreover, commonly assumed values for mutation model parameters do not appear to encompass variation in microsatellite evolution observed in vertebrates and, on average, the proportion of multi-step mutations is underestimated by a factor of approximately two. As a result, bottleneck tests can have a higher probability of 'detecting' bottlenecks in stable populations than expected based on the nominal significance level. We provide recommendations that could add rigor to inferences drawn from future bottleneck tests and highlight new directions for the characterization of demographic history.

The West Pacific diversity hotspot as a source or sink for new species? Population genetic insights from the Indo-Pacific parrotfish Scarus rubroviolaceus

We used a population genetic approach to quantify major population subdivisions and patterns of migration within a broadly distributed Indo-Pacific parrotfish. We genotyped 15 microsatellite loci in Scarus rubroviolaceus collected from 20 localities between Africa and the Americas. A STRUCTURE model indicates the presence of four major populations: Eastern Pacific, Hawaii, Central-West Pacific and a less well-differentiated Indian Ocean. We used the isolation and migration model to estimate splitting times, population sizes and migration patterns between sister population pairs. To eliminate loci under selection, we used BayeScan to select loci for three isolation and migration models: Eastern Pacific and Central-West Pacific, Hawaii and the Central-West Pacific, and Indian Ocean and the Central-West Pacific. To test the assumption of a stepwise mutation model (SMM), we used likelihood to test the SMM against a two-phase model that allowed mutational complexity. A posteriori, minor departures from SMM were estimated to affect ≤2% of the alleles in the data. The data were informative about the contemporary and ancestral population sizes, migration rates and the splitting time in the eastern Pacific/Central-West Pacific comparison. The model revealed a splitting time ∼17,000 BP, a larger contemporary N(e) in the Central-West Pacific than in the eastern Pacific and a strong bias of east to west migration. These characteristics support the Center of Accumulation model of peripatric diversification in low-diversity peripheral sites and perhaps migration from those sites to the western Pacific diversity hotspot.

Evolutionary mechanisms shaping the genetic population structure of marine fishes; lessons from the European flounder (Platichthys flesusL.)

A number of evolutionary mechanisms have been suggested for generating low but significant genetic structuring among marine fish populations. We used nine microsatellite loci and recently developed methods in landscape genetics and coalescence-based estimation of historical gene flow and effective population sizes to assess temporal and spatial dynamics of the population structure in European flounder (Platichthys flesus L.). We collected 1062 flounders from 13 localities in the northeast Atlantic and Baltic Seas and found temporally stable and highly significant genetic differentiation among samples covering a large part of the species' range (global F(ST) = 0.024, P < 0.0001). In addition to historical processes, a number of contemporary acting evolutionary mechanisms were associated with genetic structuring. Physical forces, such as oceanographic and bathymetric barriers, were most likely related with the extreme isolation of the island population at the Faroe Islands. A sharp genetic break was associated with a change in life history from pelagic to benthic spawners in the Baltic Sea. Partial Mantel tests showed that geographical distance per se was not related with genetic structuring among Atlantic and western Baltic Sea samples. Alternative factors, such as dispersal potential and/or environmental gradients, could be important for generating genetic divergence in this region. The results show that the magnitude and scale of structuring generated by a specific mechanism depend critically on its interplay with other evolutionary mechanisms, highlighting the importance of investigating species with wide geographical and ecological distributions to increase our understanding of evolution in the marine environment.

Mutation rate and pattern of microsatellites in common carp (Cyprinus carpio L.)

Microsatellites are popular molecular markers in genetic and evolutionary studies. Their mutational dynamics have been extensively studied in humans and fruit flies, but few data were available in fish. By genotyping 55 individuals of a F1 pedigree, we investigated the mutation rates and patterns of 49 microsatellites in one of the most important fresh water fish species, the common carp (Cyprinus carpio L.). The overall mutation rate of the 49 loci was 5.56 x 10(-4)/locus/generation (95% confidence interval 1.52 x 10(-4) and 1.63 x 10(-3)). The change of allele size was between +2 to -5 repeat units, assuming that the mutation allele arose from the parental allele most similar in size to the mutant.

The Fundamental Theorem of Neutral Evolution: Rates of Substitution and Mutation Should Factor in Premeiotic Clusters

Mutations do not always arise as single events. Many new mutations actually occur in the cell lineage before germ cell formation or meiosis and are therefore replicated pre-meiotically. The increased likelihood of substitutions caused by these clusters of new mutant alleles can change the fundamental theorem of neutral evolution.

Efficient isothermal expansion of human telomeric and minisatellite repeats by Thermococcus litoralis DNA polymerase

Repeating DNA sequences, such as telomeres, centromeres, and micro- and mini-satellites, comprise 50% of the genome and play important roles in regulatory and pathogenic mechanisms. In order to study structures and functions of such repeating sequences, it is important to have simple and efficient methods for making them in vitro. Here, we describe the efficient and convenient expansion of repetitive telomeric and minisatellite DNA sequences starting from small synthetic templates to final product lengths of several hundreds to thousands of nucleotides by the thermostable DNA polymerase from Thermococcus litoralis (Vent DNA polymerase). This enzyme was so far unknown to catalyze repeat expansion. Either single-stranded or double-stranded DNAs could be produced, depending on nucleotides present. Compared to earlier results obtained with other enzymes, the expansion reaction is highly efficient both in its yield and product length, and proceeds without thermal cycling. Moreover, the products are characterized by a narrow length distribution.

A high incidence of clustered microsatellite mutations revealed by parent-offspring analysis in the African freshwater snail, Bulinus forskalii (Gastropoda, Pulmonata)

Genotyping of 11 microsatellites in 432 offspring from 28 families of the hermaphroditic, freshwater snail Bulinusforskalii detected 10 de novo mutant alleles. This gave an estimated mutation rate of 1.1 x 10(-3) per locus per gamete per generation. There was a trend towards repeat length expansion and, unlike most studies, multi-step mutations predominated, suggesting that the microsatellite mutation process does not conform to a strict stepwise mutation model. Interestingly, the ten mutant alleles appear to have arisen from only six independent germline mutation events within the microsatellite array, with seven of them residing in three mutational clusters. Our results extend observations of clustered microsatellite mutations to another taxonomic group and type of mating system, self-fertile gastropods, and provide compelling evidence of premeiotic germline mutations, a phenomenon that could greatly impact upon our understanding of mutation dynamics but which has received little attention.

The impact of microsatellite electromorph size homoplasy on multilocus population structure estimates in a tropical tree (Corythophora alta) and an anadromous fish (Morone saxatilis)

Microsatellite allelic states are determined by electrophoretic sizing of polymerase chain reaction fragments to define electromorphs. Numerous studies have documented that identical microsatellite electromorphs are potentially heterogeneous at the DNA sequence level, a phenomenon called electromorph size homoplasy. Few studies have examined the impact of electromorph size homoplasy on estimates of population genetic parameters. We investigated the frequency of microsatellite electromorph size homoplasy for 12 loci in the tropical tree Corythophora alta and 11 loci in the anadromous fish Morone saxatilis by sequencing 14-23 homozygotes per locus sampled from multiple populations for a total of 453 sequences. Sequencing revealed no homoplasy for M. saxatilis loci. Seven C. alta loci exhibited homoplasy, including single and compound repeat motifs both with and without interruptions. Between 12.5 and 42.9% of electromorphs sampled per locus showed size homoplasy. Two methods of correction for homoplasy in C. alta generally produced little or no change in single-locus estimates of RST, except for two loci in which some additional differentiation among populations was revealed. Twelve-locus estimates of RST (including the seven loci corrected for homoplasy) were slightly greater than estimates from uncorrected data, although the 95% confidence intervals overlapped. The frequency of methodological errors such as clerical mistakes or sample mislabelling per genotype scored was estimated at 5.4 and 7.3% for C. alta and M. saxatilis, respectively. Simulations showed that the increase in RST produced by homoplasy correction was only slightly larger than variation in RST estimates expected to be caused by methodological errors.

High somatic instability of a microsatellite locus in a clonal tree, Robinia pseudoacacia

Robinia pseudoacacia L. is a clonal tree species. To investigate a mutation within eight microsatellite loci of R. pseudoacacia, we analyzed DNA samples obtained from different leaf samples within each ramet, leaves from ramets within the genet, and seeds. Of the eight loci, locus Rops15 (AG motif) displayed hypermutability. The mutation rates of Rops15 within each ramet, among ramets within the genet, and offspring were 6.27% (ranging from 0 to 31.1%), 6.11% (from 0 to 25.0%) and 3.78% (from 0 to 10.9%), respectively. The mutation rate increased with allele size (13-71 repeat units). The mutation patterns observed in Rops15 were distinctive in two ways. First, there was a significant bias toward additions over deletions, and both addition and deletion of single repeats were dominant at alleles with lengths less than 232 bp (63 repeats). Second, for the longest allele of 248 bp (71 repeats), the number of losses was higher than the number of gains. These observations suggest that the mutation patterns of microsatellites in R. pseudoacacia may follow a generalized stepwise mutation model, and that the tendency of long alleles to mutate to shorter lengths acts to prevent infinite growth. Finally, the observation of somatic hypermutability at locus Rops15 highlights the need for caution when using highly polymorphic microsatellites for population genetic structure and paternity analysis in tree species.