An examination of spatial and temporal genetic variation in walleye pollock (Theragra chalcogramma) using allozyme, mitochondrial DNA, and microsatellite data

We used allozyme, microsatellite, and mitochondrial DNA (mtDNA) data to test for spatial and interannual genetic diversity in wall-eye pollock (Theragra chalcogramma) from six spawning aggregations representing three geographic regions: Gulf of Alaska, eastern Bering Sea, and eastern Kamchatka. Interpopulation genetic diversity was evident primarily from the mtDNA and two allozyme loci (SOD-2*, MPI*). Permutation tests ˆindicated that FST values for most allozyme and microsatellite loci were not significantly greater than zero. The microsatellite results suggested that high locus polymorphism may not be a reliable indicator of power for detecting population differentiation in walleye pollock. The fact that mtDNA revealed population structure and most nuclear loci did not suggests that the effective size of most walleye pollock populations is large (genetic drift is weak) and migration is a relatively strong homogenizing force. The allozymes and mtDNA provided mostly concordant estimates of patterns of spatial genetic variation. These data showed significant genetic variation between North American and Asian populations. In addition, two spawning aggregations in the Gulf of Alaska, in Prince William Sound, and off Middleton Island, appeared genetically distinct from walleye pollock spawning in the Shelikof Strait and may merit management as a distinct stock. Finally, we found evidence of interannual genetic variation in two of three North American spawning aggregations, similar in magnitude to the spatial variation among North American walleye pol-lock. We suggest that interannual genetic variation in walleye pollock may be indicative of one or more of the following factors: highly variable reproductive success, adult philopatry, source-sink metapopulation structure, and intraannual variation (days) in spawning timing among genetically distinct but spatially identical spawning aggregates.

Are the native giant tortoises from the Seychelles really extinct? A genetic perspective based on mtDNA and microsatellite data.

The extinction of the giant tortoises of the Seychelles Archipelago has long been suspected but is not beyond doubt. A recent morphological study of the giant tortoises of the western Indian Ocean concluded that specimens of two native Seychelles species survive in captivity today alongside giant tortoises of Aldabra, which are numerous in zoos as well as in the wild. This claim has been controversial because some of the morphological characters used to identify these species, several measures of carapace morphology, are reputed to be quite sensitive to captive conditions. Nonetheless, the potential survival of giant tortoise species previously thought extinct presents an exciting scenario for conservation. We used mitochondrial DNA sequences and nuclear microsatellites to examine the validity of the rediscovered species of Seychelles giant tortoises. Our results indicate that the morphotypes suspected to represent Seychelles species do not show levels of variation and genetic structuring consistent with long periods of reproductive isolation. We found no variation in the mitochondrial control region among 55 individuals examined and no genetic structuring in eight microsatellite loci, pointing to the survival of just a single lineage of Indian Ocean tortoises.

Extensive genome duplications in sturgeons:New evidence from microsatellite data

Summary: Genome duplications and polyploidization events are thought to have played relevant roles in the early stages of vertebrate evolution, in particular near the time of divergence of the lamprey lineage. Additional genome duplications, specifically in ray-finned fish, may have occurred before the divergence of the teleosts. The role of polyploidization in vertebrate genome evolution is a thriving area of research. Sturgeons (order Acipenseriformes) provide a unique model for the investigation of genome duplication, with existing species possessing 120, 250 or 360 chromosomes. In the present study, data from 240 sturgeon specimens representing 11 species were used for analysis of ploidy levels. Allele numbers were assessed at eleven microsatellite loci. The results provide further evidence for functional diploidy, tetraploidy and hexaploidy in species possessing 120, 250 and 360 chromosomes, respectively. The analysis also uncovered novel evidence for functional hexaploidy in the shortnose sturgeon (Acipenser brevirostrum). In conclusion, the process of functional genome reduction is demonstrated to be an on-going process in this fish lineage. © 2013 Blackwell Verlag GmbH.

Rapid radiation in spiny lobsters (Palinurus spp) as revealed by classic and ABC methods using mtDNA and microsatellite data

Background: Molecular tools may help to uncover closely related and still diverging species from a wide variety of taxa and provide insight into the mechanisms, pace and geography of marine speciation. There is a certain controversy on the phylogeography and speciation modes of species-groups with an Eastern Atlantic-Western Indian Ocean distribution, with previous studies suggesting that older events (Miocene) and/or more recent (Pleistocene) oceanographic processes could have influenced the phylogeny of marine taxa. The spiny lobster genus Palinurus allows for testing among speciation hypotheses, since it has a particular distribution with two groups of three species each in the Northeastern Atlantic (P. elephas, P. mauritanicus and P. charlestoni) and Southeastern Atlantic and Southwestern Indian Oceans (P. gilchristi, P. delagoae and P. barbarae). In the present study, we obtain a more complete understanding of the phylogenetic relationships among these species through a combined dataset with both nuclear and mitochondrial markers, by testing alternative hypotheses on both the mutation rate and tree topology under the recently developed approximate Bayesian computation (ABC) methods. Results: Our analyses support a North-to-South speciation pattern in Palinurus with all the South-African species forming a monophyletic clade nested within the Northern Hemisphere species. Coalescent-based ABC methods allowed us to reject the previously proposed hypothesis of a Middle Miocene speciation event related with the closure of the Tethyan Seaway. Instead, divergence times obtained for Palinurus species using the combined mtDNA-microsatellite dataset and standard mutation rates for mtDNA agree with known glaciation-related processes occurring during the last 2 my. Conclusion: The Palinurus speciation pattern is a typical example of a series of rapid speciation events occurring within a group, with very short branches separating different species. Our results support the hypothesis that recent climate change-related oceanographic processes have influenced the phylogeny of marine taxa, with most Palinurus species originating during the last two million years. The present study highlights the value of new coalescent-based statistical methods such as ABC for testing different speciation hypotheses using molecular data.

A note on the accuracy of PAC-likelihood inference with microsatellite data

Stephens and Donnelly have introduced a simple yet powerful importance sampling scheme for computing the likelihood in population genetic models. Fundamental to the method is an approximation to the conditional probability of the allelic type of an additional gene, given those currently in the sample. As noted by Li and Stephens, the product of these conditional probabilities for a sequence of draws that gives the frequency of allelic types in a sample is an approximation to the likelihood, and can be used directly in inference. The aim of this note is to demonstrate the high level of accuracy of "product of approximate conditionals" (PAC) likelihood when used with microsatellite data. Results obtained on simulated microsatellite data show that this strategy leads to a negligible bias over a wide range of the scaled mutation parameter theta. Furthermore, the sampling variance of likelihood estimates as well as the computation time are lower than that obtained with importance sampling on the whole range of theta. It follows that this approach represents an efficient substitute to IS algorithms in computer intensive (e.g. MCMC) inference methods in population genetics. (c) 2006 Elsevier Inc. All rights reserved.

Allelic size standards and reference genotypes to unify international cocoa (Theobroma cacao L.) microsatellite data

Standardisation of microsatellite allele profiles between laboratories is of fundamental importance to the transferability of genetic fingerprint data and the identification of clonal individuals held at multiple sites. Here we describe two methods of standardisation applied to the microsatellite fingerprinting of 429 Theobroma cacao L. trees representing 345 accessions held in the worlds largest Cocoa Intermediate Quarantine facility: the use of a partial allelic ladder through the production of 46 cloned and sequenced allelic standards (AJ748464 to AJ48509), and the use of standard genotypes selected to display a diverse allelic range. Until now a lack of accurate and transferable identification information has impeded efforts to genetically improve the cocoa crop. To address this need, a global initiative to fingerprint all international cocoa germplasm collections using a common set of 15 microsatellite markers is in progress. Data reported here have been deposited with the International Cocoa Germplasm Database and form the basis of a searchable resource for clonal identification. To our knowledge, this is the first quarantine facility to be completely genotyped using microsatellite markers for the purpose of quality control and clonal identification. Implications of the results for retrospective tracking of labelling errors are briefly explored.

Sturnus vulgaris, mitochondrial selection, MATLAB code and microsatellite data

Here we provide MATLAB code used to simulate drift and selection between and within individuals, which has been used to investigate mitochondrial haplotype frequency shifts in Sturnus vulgaris. Also provided is a microsatellite data set used to assess whether empirical allele frequency shifts were likely to be caused by admixture. These files support and upcoming publication, which concludes that within-individuals selection on mitochondrial DNA best explains empirical data.

A new data source for fisheries resource assessment: genetic estimates of the effective number of spawners. Final Report to the Fisheries Research and Development Corporation.

The development of innovative methods of stock assessment is a priority for State and Commonwealth fisheries agencies. It is driven by the need to facilitate sustainable exploitation of naturally occurring fisheries resources for the current and future economic, social and environmental well being of Australia. This project was initiated in this context and took advantage of considerable recent achievements in genomics that are shaping our comprehension of the DNA of humans and animals. The basic idea behind this project was that genetic estimates of effective population size, which can be made from empirical measurements of genetic drift, were equivalent to estimates of the successful number of spawners that is an important parameter in process of fisheries stock assessment. The broad objectives of this study were to 1. Critically evaluate a variety of mathematical methods of calculating effective spawner numbers (Ne) by a. conducting comprehensive computer simulations, and by b. analysis of empirical data collected from the Moreton Bay population of tiger prawns (P. esculentus). 2. Lay the groundwork for the application of the technology in the northern prawn fishery (NPF). 3. Produce software for the calculation of Ne, and to make it widely available. The project pulled together a range of mathematical models for estimating current effective population size from diverse sources. Some of them had been recently implemented with the latest statistical methods (eg. Bayesian framework Berthier, Beaumont et al. 2002), while others had lower profiles (eg. Pudovkin, Zaykin et al. 1996; Rousset and Raymond 1995). Computer code and later software with a user-friendly interface (NeEstimator) was produced to implement the methods. This was used as a basis for simulation experiments to evaluate the performance of the methods with an individual-based model of a prawn population. Following the guidelines suggested by computer simulations, the tiger prawn population in Moreton Bay (south-east Queensland) was sampled for genetic analysis with eight microsatellite loci in three successive spring spawning seasons in 2001, 2002 and 2003. As predicted by the simulations, the estimates had non-infinite upper confidence limits, which is a major achievement for the application of the method to a naturally-occurring, short generation, highly fecund invertebrate species. The genetic estimate of the number of successful spawners was around 1000 individuals in two consecutive years. This contrasts with about 500,000 prawns participating in spawning. It is not possible to distinguish successful from non-successful spawners so we suggest a high level of protection for the entire spawning population. We interpret the difference in numbers between successful and non-successful spawners as a large variation in the number of offspring per family that survive – a large number of families have no surviving offspring, while a few have a large number. We explored various ways in which Ne can be useful in fisheries management. It can be a surrogate for spawning population size, assuming the ratio between Ne and spawning population size has been previously calculated for that species. Alternatively, it can be a surrogate for recruitment, again assuming that the ratio between Ne and recruitment has been previously determined. The number of species that can be analysed in this way, however, is likely to be small because of species-specific life history requirements that need to be satisfied for accuracy. The most universal approach would be to integrate Ne with spawning stock-recruitment models, so that these models are more accurate when applied to fisheries populations. A pathway to achieve this was established in this project, which we predict will significantly improve fisheries sustainability in the future. Regardless of the success of integrating Ne into spawning stock-recruitment models, Ne could be used as a fisheries monitoring tool. Declines in spawning stock size or increases in natural or harvest mortality would be reflected by a decline in Ne. This would be good for data-poor fisheries and provides fishery independent information, however, we suggest a species-by-species approach. Some species may be too numerous or experiencing too much migration for the method to work. During the project two important theoretical studies of the simultaneous estimation of effective population size and migration were published (Vitalis and Couvet 2001b; Wang and Whitlock 2003). These methods, combined with collection of preliminary genetic data from the tiger prawn population in southern Gulf of Carpentaria population and a computer simulation study that evaluated the effect of differing reproductive strategies on genetic estimates, suggest that this technology could make an important contribution to the stock assessment process in the northern prawn fishery (NPF). Advances in the genomics world are rapid and already a cheaper, more reliable substitute for microsatellite loci in this technology is available. Digital data from single nucleotide polymorphisms (SNPs) are likely to super cede ‘analogue’ microsatellite data, making it cheaper and easier to apply the method to species with large population sizes.

Microsatellite standardization and evaluation of genotyping error in a large multi-partner research programme for conservation of Atlantic salmon (Salmo salar L.).

Microsatellite genotyping is a common DNA characterization technique in population, ecological and evolutionary genetics research. Since different alleles are sized relative to internal size-standards, different laboratories must calibrate and standardize allelic designations when exchanging data. This interchange of microsatellite data can often prove problematic. Here, 16 microsatellite loci were calibrated and standardized for the Atlantic salmon, Salmo salar, across 12 laboratories. Although inconsistencies were observed, particularly due to differences between migration of DNA fragments and actual allelic size ('size shifts'), inter-laboratory calibration was successful. Standardization also allowed an assessment of the degree and partitioning of genotyping error. Notably, the global allelic error rate was reduced from 0.05 ± 0.01 prior to calibration to 0.01 ± 0.002 post-calibration. Most errors were found to occur during analysis (i.e. when size-calling alleles; the mean proportion of all errors that were analytical errors across loci was 0.58 after calibration). No evidence was found of an association between the degree of error and allelic size range of a locus, number of alleles, nor repeat type, nor was there evidence that genotyping errors were more prevalent when a laboratory analyzed samples outside of the usual geographic area they encounter. The microsatellite calibration between laboratories presented here will be especially important for genetic assignment of marine-caught Atlantic salmon, enabling analysis of marine mortality, a major factor in the observed declines of this highly valued species.

A novel method of microsatellite genotyping-by-sequencing using individual combinatorial barcoding

This study examines the potential of next-generation sequencing based ‘genotyping-by-sequencing’ (GBS) of microsatellite loci for rapid and cost-effective genotyping in large-scale population genetic studies. The recovery of individual genotypes from large sequence pools was achieved by PCR-incorporated combinatorial barcoding using universal primers. Three experimental conditions were employed to explore the possibility of using this approach with existing and novel multiplex marker panels and weighted amplicon mixture. The GBS approach was validated against microsatellite data generated by capillary electrophoresis. GBS allows access to the underlying nucleotide sequences that can reveal homoplasy, even in large datasets and facilitates cross laboratory transfer. GBS of microsatellites, using individual combinatorial barcoding, is potentially faster and cheaper than current microsatellite approaches and offers better and more data.

Cross-species amplification and phylogenetic reconstruction using fragaria microsatellite primers

Twenty-eight microsatellite primer pairs developed from Fragaria vesca ‘Rügen’ were applied to sixteen accessions representing eight diploid Fragaria species. The number of alleles generated, the power of discrimination and the percentage of accessions where no PCR product could be amplified were calculated for each locus for the thirteen non-F. vesca accessions. A phylogeny was then generated for the species accessions sampled, using the presence or absence of alleles at the polymorphic loci as character states. Despite the problems inherent in phylogeny reconstruction from microsatellite data, the phylogeny showed some congruence with a previously published phylogeny of Fragaria, based on nucleotide sequence data. However, relationships inferred from microsatellite allele data were relatively unresolved and poorly supported. The genetic basis of allelic polymorphisms at specific loci was investigated through direct sequencing of the PCR products amplified by three primer pairs. The potential utility of sequence data generated from microsatellite loci in evolutionary studies of closely related species groups is briefly explored.

Microsatellite variation in the yellowhammer Emberiza citrinella : population structure of a declining farmland bird

In recent years, there has been much concern in the UK about population declines of widespread species in agricultural habitats. Conservation-orientated research on declining birds has focused on vital rates of survival and productivity. However, the environmental factors which may influence movements between populations of widespread species is poorly understood. Population genetic structure is an indirect description of dispersal between groups of individuals. To attempt to develop an understanding of genetic structuring in a widespread, but declining, farmland bird, we therefore investigated the yellowhammer, Emberiza citrinella, population in England and Wales using microsatellite data. Our first aim was to investigate whether there was genetic substructuring in the population. A second aim was to investigate if there was a relationship between genetic distances and various environmental variables. Finally, we analysed the microsatellite data for evidence of loss of genetic variation due to population decline. Our data showed a slight but significant structure within the yellowhammer population. This therefore cannot be considered a panmictic population. Our example from South Cumbria implies that high-altitude barriers may have a slight influence on population structure. However, on the whole, genetic distances between sample sites were not significantly correlated with geographical distances, degrees of population connectivity, high altitudes, or differences in precipitation between sites. Finally, we detected departures from mutation-drift equilibrium (excess heterozygosity), which is indicative of a loss of genetic variation through recent decline.

Genetic analysis reveals two stocks of skipjack tuna (Katsuwonus pelamis) in the north western Indian Ocean

Skipjack (SJT) (Katsuwonus pelamis) is a medium sized, pelagic, highly dispersive tuna species that occurs widely across tropical and subtropical waters. SJT constitute the largest tuna fishery in the Indian Ocean, and are currently managed as a single stock. Patterns of genetic variation in a mtDNA gene and 6 microsatellite loci were examined to test for stock structure in the northwestern Indian Ocean. 324 individuals were sampled from five major fishing grounds around Sri Lanka, and single sites in the Maldive Islands and the Laccadive Islands. Phylogenetic reconstruction of mtDNA revealed two coexisting divergent clades in the region. AMOVA (Analysis of Molecular Variance) of mtDNA data revealed significant genetic differentiation among sites (ΦST = 0.2029, P < 0.0001), also supported by SAMOVA results. AMOVA of microsatellite data also showed significant differentiation among most sampled sites (FST = 0.0256, P<0.001) consistent with the mtDNA pattern. STRUCTURE analysis of the microsatellite data revealed two differentiated stocks. While the both two marker types examined identified two genetic groups, microsatellite analysis indicates that the sampled SJT are likely to represent individuals sourced from discrete breeding grounds that are mixed in feeding grounds in Sri Lankan waters.

Evolutionary genetics and human assisted movement of a globally invasive pest (Russian wheat aphid : Diuraphis noxia)

This PhD study has examined the population genetics of the Russian wheat aphid (RWA, Diuraphis noxia), one of the world’s most invasive agricultural pests, throughout its native and introduced global range. Firstly, this study investigated the geographic distribution of genetic diversity within and among RWA populations in western China. Analysis of mitochondrial data from 18 sites provided evidence for the long-term existence and expansion of RWAs in western China. The results refute the hypothesis that RWA is an exotic species only present in China since 1975. The estimated date of RWA expansion throughout western China coincides with the debut of wheat domestication and cultivation practices in western Asia in the Holocene. It is concluded that western China represents the limit of the far eastern native range of this species. Analysis of microsatellite data indicated high contemporary gene flow among northern populations in western China, while clear geographic isolation between northern and southern populations was identified across the Tianshan mountain range and extensive desert regions. Secondly, this study analyzed the worldwide pathway of invasion using both microsatellite and endosymbiont genetic data. Individual RWAs were obtained from native populations in Central Asia and the Middle East and invasive populations in Africa and the Americas. Results indicated two pathways of RWA invasion from 1) Syria in the Middle East to North Africa and 2) Turkey to South Africa, Mexico and then North and South America. Very little clone diversity was identified among invasive populations suggesting that a limited founder event occurred together with predominantly asexual reproduction and rapid population expansion. The most likely explanation for the rapid spread (within two years) from South Africa to the New World is by human movement, probably as a result of the transfer of wheat breeding material. Furthermore, the mitochondrial data revealed the presence of a universal haplotype and it is proposed that this haplotype is representative of a wheat associated super-clone that has gained dominance worldwide as a result of the widespread planting of domesticated wheat. Finally, this study examined salivary gland gene diversity to determine whether a functional basis for RWA invasiveness could be identified. Peroxidase DNA sequence data were obtained for a selection of worldwide RWA samples. Results demonstrated that most native populations were polymorphic while invasive populations were monomorphic, supporting previous conclusions relating to demographic founder effects in invasive populations. Purifying selection most likely explains the existence of a universal allele present in Middle Eastern populations, while balancing selection was evident in East Asian populations. Selection acting on the peroxidase gene may provide an allele-dependent advantage linked to the successful establishment of RWAs on wheat, and ultimately their invasion potential. In conclusion, this study is the most comprehensive molecular genetic investigation of RWA population genetics undertaken to date and provides significant insights into the source and pathway of global invasion and the potential existence of a wheat-adapted genotype that has colonised major wheat growing countries worldwide except for Australia. This research has major biosecurity implications for Australia’s grain industry.

Pleistocene climate fluctuations influence phylogeographical patterns in Melomys cervinipes across the mesic forests of eastern Australia

Aim Our aim was to clarify the lineage-level relationships for Melomys cervinipes and its close relatives and investigate whether the patterns of divergence observed for these wet-forest-restricted mammals may be associated with recognized biogeographical barriers. Location Mesic closed forest along the east coast of Australia, from north Queensland to mid-eastern New South Wales. Methods To enable rigorous phylogenetic reconstruction, divergence-date estimation and phylogeographical inference, we analysed DNA sequence and microsatellite data from 307 specimens across the complete distribution of M. cervinipes (45 localities). Results Three divergent genetic lineages were found within M. cervinipes, corresponding to geographically delineated northern, central and southern clades. Additionally, a fourth lineage, comprising M. rubicola and M. capensis, was identified and was most closely related to the northern M. cervinipes lineage. Secondary contact of the northern and central lineages was identified at one locality to the north of the Burdekin Gap. Main conclusions Contemporary processes of repeated habitat fragmentation and contraction, local extinction events and subsequent re-expansion across both small and large areas, coupled with the historical influence of the Brisbane Valley Barrier, the St Lawrence Gap and the Burdekin Gap, have contributed to the present phylogeographical structure within M. cervinipes. Our study highlights the need to sample close to the periphery of putative biogeographical barriers or risk missing vital phylogeographical information that may significantly alter the interpretation of biogeographical hypotheses.