Population genetics of Australian white sharks reveals fine-scale spatial structure, transoceanic dispersal events and low effective population sizes

Despite international protection of white sharks Carcharodon carcharias, important conservation parameters such as abundance, population structure and genetic diversity are largely unknown. The tissue of 97 predominately juvenile white sharks sampled from spatially distant eastern and southwestern Australian coastlines was sequenced for the mitochondrial DNA (mtDNA) control region and genotyped with 6 nuclear-encoded microsatellite loci. MtDNA population structure was found between the eastern and southwestern coasts (F-ST = 0.142, p < 0.0001), implying female reproductive philopatry. This concurs with recent satellite and acoustic tracking findings which suggest the sustained presence of discrete east coast nursery areas. Furthermore, population subdivision was found between the same regions with biparentally inherited micro satellite markers (F-ST = 0.009, p < 0.05), suggesting that males may also exhibit some degree of reproductive philopatry; 5 sharks captured along the east coast had mtDNA haplotypes that resembled western Indian Ocean sharks more closely than Australian/New Zealand sharks, suggesting that transoceanic dispersal, or migration resulting in breeding, may occur sporadically. Our most robust estimate of contemporary genetic effective population size was low and close to thresholds at which adaptive potential may be lost. For a variety of reasons, these contemporary estimates were at least 1, possibly 2, orders of magnitude below our historical effective size estimates. Population decline could expose these genetically isolated populations to detrimental genetic effects. Regional Australian white shark conservation management units should be implemented until genetic population structure, size and diversity can be investigated in more detail.

Maximum likelihood estimation of a migration matrix and effective population sizes in n subpopulations by using a coalescent approach

A maximum likelihood estimator based on the coalescent for unequal migration rates and different subpopulation sizes is developed. The method uses a Markov chain Monte Carlo approach to investigate possible genealogies with branch lengths and with migration events. Properties of the new method are shown by using simulated data from a four-population n-island model and a source–sink population model. Our estimation method as coded in migrate is tested against genetree; both programs deliver a very similar likelihood surface. The algorithm converges to the estimates fairly quickly, even when the Markov chain is started from unfavorable parameters. The method was used to estimate gene flow in the Nile valley by using mtDNA data from three human populations.

Spatial population genetic structure reveals strong natal site fidelity in Echinocladius martini (Diptera: Chironomidae) in northeast Queensland, Australia

1. A diverse array of patterns has been reported regarding the spatial extent of population genetic structure and effective dispersal in freshwater macroinvertebrates. In river systems, the movements of many taxa can be restricted to varying degrees by the natural stream channel hierarchy. 2. In this study, we sampled populations of the non-biting freshwater midge Echinocladius martini in the Paluma bioregion of tropical northeast Queensland to investigate fine scale patterns of within- and among-stream dispersal and gene flow within a purported historical refuge. We amplified a 639 bp fragment of mitochondrial COI and analysed genetic structure using pairwise ΦST, hierarchical AMOVA, Mantel tests and a parsimony network. Genetic variation was partitioned among stream sections using Streamtree to investigate the effect of potential instream dispersal barriers. 3. The data revealed strong natal site fidelity and significant differentiation among neighbouring, geographically proximate streams. We found evidence for only episodic adult flight among sites on separate stream reaches. Overall, however, our data suggested that both larval and adult dispersal was largely limited to within a stream channel. 4. This may arise from a combination of the high density of riparian vegetation physically restricting dispersal and from the joint effects of habitat stability and large population sizes. Together these may mitigate the requirement for movement among streams to avoid inbreeding and local extinction due to habitat change and may thus enable persistence of upstream populations in the absence of regular compensatory upstream flight. Taken together, these data suggest that dispersal of E. martini is highly restricted, to the scale of only a few kilometres, and hence occurs predominantly within the natal stream.

Population genetic structure and patterns of dispersal in the Giant Long-Armed Prawn, Macrobrachium lar (Fabricius, 1798) (Decapoda : Palaemonidae)

The Giant Long-Armed Prawn, Macrobrachium lar is a freshwater species native to the Indo-Pacific. M. lar has a long-lived, passive, pelagic marine larval stage where larvae need to colonise freshwater within three months to complete their development. Dispersal is likely to be influenced by the extensive distances larvae must transit between small oceanic islands to find suitable freshwater habitat, and by prevailing east to west wind and ocean currents in the southern Pacific Ocean. Thus, both intrinsic and extrinsic factors are likely to influence wild population structure in this species. The present study sought to define the contemporary broad and fine-scale population genetic structure of Macrobrachium lar in the south-western Pacific Ocean. Three polymorphic microsatellite loci were used to assess patterns of genetic variation within and among 19 wild adult sample sites. Statistical procedures that partition variation implied that at both spatial scales, essentially all variation was present within sample sites and differentiation among sites was low. Any differentiation observed also was not correlated with geographical distance. Statistical approaches that measure genetic distance, at the broad-scale, showed that all south-western Pacific Islands were essentially homogeneous, with the exception of a well supported divergent Cook Islands group. These findings are likely the result of some combination of factors that may include the potential for allelic homoplasy, through to the effects of sampling regime. Based on the findings, there is most likely a divergent M. lar Cook Islands clade in the south-western Pacific Ocean, resulting from prevailing ocean currents. Confirmation of this pattern will require a more detailed analysis of nDNA variation using a larger number of loci and, where possible, use of larger population sizes.

Computational experiments involving population size for FPGA-based implementation of a GA for the TSP

The feasibility of using an in-hardware implementation of a genetic algorithm (GA) to solve the computationally expensive travelling salesman problem (TSP) is explored, especially in regard to hardware resource requirements for problem and population sizes. We investigate via numerical experiments whether a small population size might prove sufficient to obtain reasonable quality solutions for the TSP, thereby permitting relatively resource efficient hardware implementation on field programmable gate arrays (FPGAs). Software experiments on two TSP benchmarks involving 48 and 532 cities were used to explore the extent to which population size can be reduced without compromising solution quality, and results show that a GA allowed to run for a large number of generations with a smaller population size can yield solutions of comparable quality to those obtained using a larger population. This finding is then used to investigate feasible problem sizes on a targeted Virtex-7 vx485T-2 FPGA platform via exploration of hardware resource requirements for memory and data flow operations.

A Model to Estimate the Population Contributing to the Wastewater Using Samples Collected on Census Day

An important uncertainty when estimating per capita consumption of, for example, illicit drugs by means of wastewater analysis (sometimes referred to as “sewage epidemiology”) relates to the size and variability of the de facto population in the catchment of interest. In the absence of a day-specific direct population count any indirect surrogate model to estimate population size lacks a standard to assess associated uncertainties. Therefore, the objective of this study was to collect wastewater samples at a unique opportunity, that is, on a census day, as a basis for a model to estimate the number of people contributing to a given wastewater sample. Mass loads for a wide range of pharmaceuticals and personal care products were quantified in influents of ten sewage treatment plants (STP) serving populations ranging from approximately 3500 to 500 000 people. Separate linear models for population size were estimated with the mass loads of the different chemical as the explanatory variable: 14 chemicals showed good, linear relationships, with highest correlations for acesulfame and gabapentin. De facto population was then estimated through Bayesian inference, by updating the population size provided by STP staff (prior knowledge) with measured chemical mass loads. Cross validation showed that large populations can be estimated fairly accurately with a few chemical mass loads quantified from 24-h composite samples. In contrast, the prior knowledge for small population sizes cannot be improved substantially despite the information of multiple chemical mass loads. In the future, observations other than chemical mass loads may improve this deficit, since Bayesian inference allows including any kind of information relating to population size.

Extremely Low Microsatellite Diversity but Distinct Population Structure in a Long-Lived Threatened Species, the Australian Lungfish Neoceratodus forsteri (Dipnoi)

The Australian lungfish is a unique living representative of an ancient dipnoan lineage, listed as ‘vulnerable’ to extinction under Australia’s Environment Protection and Biodiversity Conservation Act 1999. Historical accounts indicate this species occurred naturally in two adjacent river systems in Australia, the Burnett and Mary. Current day populations in other rivers are thought to have arisen by translocation from these source populations. Early genetic work detected very little variation and so had limited power to answer questions relevant for management including how genetic variation is partitioned within and among sub-populations. In this study, we use newly developed microsatellite markers to examine samples from the Burnett and Mary Rivers, as well as from two populations thought to be of translocated origin, Brisbane and North Pine. We test whether there is significant genetic structure among and within river drainages; assign putatively translocated populations to potential source populations; and estimate effective population sizes. Eleven polymorphic microsatellite loci genotyped in 218 individuals gave an average within-population heterozygosity of 0.39 which is low relative to other threatened taxa and for freshwater fishes in general. Based on F_ST values (average over loci = 0.11) and STRUCTURE analyses, we identify three distinct populations in the natural range, one in the Burnett and two distinct populations in the Mary. These analyses also support the hypothesis that the Mary River is the likely source of translocated populations in the Brisbane and North Pine rivers, which agrees with historical published records of a translocation event giving rise to these populations. We were unable to obtain bounded estimates of effective population size, as we have too few genotype combinations, although point estimates were low, ranging from 29 - 129. We recommend that, in order to preserve any local adaptation in the three distinct populations that they be managed separately.

Phylogeography and Population Genetics of the Moor frog (Rana arvalis Nilsson) in Northern Europe

The ongoing climate change along with increasing levels of pollutants, diseases, habitat loss and fragmentation constitute global threats to the persistence of many populations, species and ecosystems. However, for the long-term persistence of local populations, one of the biggest threats is the intrinsic loss of genetic variation. In order to adapt to changes in the environment, organisms must have a sufficient supply of heritable variation in traits important for their fitness. With a loss of genetic variation, the risk of extinction will increase. For conservational practices, one should therefore understand the processes that shape the genetic population structure and also the broader (historical) phylogenetic patterning of the species in focus. In this thesis, microsatellite markers were applied to study genetic diversity and population differentiation of the protected moor frog (Rana arvalis) in Fennoscandia from both historical (evolutionary) and applied (conservation) perspectives. The results demonstrate that R. arvalis populations are highly structured over rather short geographic distances. Moreover, the results suggest that R. arvalis recolonized Fennoscandia from two directions after the last ice age. This has had implications for the genetic structuring and population differentiation, especially in the northernmost parts where the two lineages have met. Compared to more southern populations, the genetic variation decreases and the interpopulation differentiation increases dramatically towards north. This could be an outcome of serial population bottlenecking along the recolonization route. Also, current isolation and small population sizes increase the effect of drift, thus reinforcing the observed pattern. The same pattern can also be seen in island populations. However, though R. arvalis on the island of Gotland has lost most of its neutral genetic variability, our results indicate that the levels of additive genetic variation have remained high. This conforms to the conjecture that though neutral markers are widely used in conservation purposes, they may be quite uninformative about the levels of genetic variation in ecologically important traits. Finally, the evolutionary impact of the typical amphibian mating behaviour on genetic diversity was investigated. Given the short time available for larval development, it is important that mating takes place as early as possible. The genetic data and earlier capture-recapture data suggest that R. arvalis gather at mating grounds they are familiar with. However, by forming leks in random to relatedness, and having multiple paternities in single clutches, the risk of inbreeding may be minimized in this otherwise highly philopatric species.

A finite population model of mlecular evolution: theory and computation

This article is concerned with the evolution of haploid organisms that reproduce asexually. In a seminal piece of work, Eigen and coauthors proposed the quasispecies model in an attempt to understand such an evolutionary process. Their work has impacted antiviral treatment and vaccine design strategies. Yet, predictions of the quasispecies model are at best viewed as a guideline, primarily because it assumes an infinite population size, whereas realistic population sizes can be quite small. In this paper we consider a population genetics-based model aimed at understanding the evolution of such organisms with finite population sizes and present a rigorous study of the convergence and computational issues that arise therein. Our first result is structural and shows that, at any time during the evolution, as the population size tends to infinity, the distribution of genomes predicted by our model converges to that predicted by the quasispecies model. This justifies the continued use of the quasispecies model to derive guidelines for intervention. While the stationary state in the quasispecies model is readily obtained, due to the explosion of the state space in our model, exact computations are prohibitive. Our second set of results are computational in nature and address this issue. We derive conditions on the parameters of evolution under which our stochastic model mixes rapidly. Further, for a class of widely used fitness landscapes we give a fast deterministic algorithm which computes the stationary distribution of our model. These computational tools are expected to serve as a framework for the modeling of strategies for the deployment of mutagenic drugs.

Population structure, long-term connectivity, and effective size of mutton snapper (Lutjanus analis) in the Caribbean Sea and Florida Keys

Genetic structure and average long-term connectivity and effective size of mutton snapper (Lutjanus analis) sampled from offshore localities in the U.S. Caribbean and the Florida Keys were assessed by using nuclear-encoded microsatellites and a fragment of mitochondrial DNA. No significant differences in allele, genotype (microsatellites), or haplotype (mtDNA) distributions were detected; tests of selective neutrality (mtDNA) were nonsignificant after Bonferroni correction. Heuristic estimates of average long-term rate of migration (proportion of migrant individuals/generation) between geographically adjacent localities varied from 0.0033 to 0.0054, indicating that local subpopulations could respond independently of environmental perturbations. Estimates of average longterm effective population sizes varied from 341 to 1066 and differed significantly among several of the localities. These results indicate that over time larval drift and interregional adult movement may not be sufficient to maintain population sustainability across the region and that there may be different demographic stocks at some of the localities studied. The estimate of long-term effective population size at the locality offshore of St. Croix was below the minimum threshold size considered necessary to maintain the equilibrium between the loss of adaptive genetic variance from genetic drift and its replacement by mutation. Genetic variability in mutton snapper likely is maintained at the intraregional level by aggregate spawning and random mating of local populations. This feature is perhaps ironic in that aggregate spawning also renders mutton snapper especially vulnerable to overexploitation.

New Nuclear SNP Markers Unravel the Genetic Structure and Effective Population Size of Albacore Tuna (Thunnus alalunga)

In the present study we have investigated the population genetic structure of albacore (Thunnus alalunga, Bonnaterre 1788) and assessed the loss of genetic diversity, likely due to overfishing, of albacore population in the North Atlantic Ocean. For this purpose, 1,331 individuals from 26 worldwide locations were analyzed by genotyping 75 novel nuclear SNPs. Our results indicated the existence of four genetically homogeneous populations delimited within the Mediterranean Sea, the Atlantic Ocean, the Indian Ocean and the Pacific Ocean. Current definition of stocks allows the sustainable management of albacore since no stock includes more than one genetic entity. In addition, short-and long-term effective population sizes were estimated for the North Atlantic Ocean albacore population, and results showed no historical decline for this population. Therefore, the genetic diversity and, consequently, the adaptive potential of this population have not been significantly affected by overfishing.

Impact of simian immunodeficiency virus infection on chimpanzee population dynamics.

Like human immunodeficiency virus type 1 (HIV-1), simian immunodeficiency virus of chimpanzees (SIVcpz) can cause CD4+ T cell loss and premature death. Here, we used molecular surveillance tools and mathematical modeling to estimate the impact of SIVcpz infection on chimpanzee population dynamics. Habituated (Mitumba and Kasekela) and non-habituated (Kalande) chimpanzees were studied in Gombe National Park, Tanzania. Ape population sizes were determined from demographic records (Mitumba and Kasekela) or individual sightings and genotyping (Kalande), while SIVcpz prevalence rates were monitored using non-invasive methods. Between 2002-2009, the Mitumba and Kasekela communities experienced mean annual growth rates of 1.9% and 2.4%, respectively, while Kalande chimpanzees suffered a significant decline, with a mean growth rate of -6.5% to -7.4%, depending on population estimates. A rapid decline in Kalande was first noted in the 1990s and originally attributed to poaching and reduced food sources. However, between 2002-2009, we found a mean SIVcpz prevalence in Kalande of 46.1%, which was almost four times higher than the prevalence in Mitumba (12.7%) and Kasekela (12.1%). To explore whether SIVcpz contributed to the Kalande decline, we used empirically determined SIVcpz transmission probabilities as well as chimpanzee mortality, mating and migration data to model the effect of viral pathogenicity on chimpanzee population growth. Deterministic calculations indicated that a prevalence of greater than 3.4% would result in negative growth and eventual population extinction, even using conservative mortality estimates. However, stochastic models revealed that in representative populations, SIVcpz, and not its host species, frequently went extinct. High SIVcpz transmission probability and excess mortality reduced population persistence, while intercommunity migration often rescued infected communities, even when immigrating females had a chance of being SIVcpz infected. Together, these results suggest that the decline of the Kalande community was caused, at least in part, by high levels of SIVcpz infection. However, population extinction is not an inevitable consequence of SIVcpz infection, but depends on additional variables, such as migration, that promote survival. These findings are consistent with the uneven distribution of SIVcpz throughout central Africa and explain how chimpanzees in Gombe and elsewhere can be at equipoise with this pathogen.

High dispersal potential has maintained long-term population stability in the North Atlantic copepod Calanus finmarchicus

The cool-water copepod Calanus finmarchicus is a key species in North Atlantic marine ecosystems since it represents an important food resource for the developmental stages of several fish of major economic value. Over the last 40 years, however, data from the Continuous Plankton Recorder survey have highlighted a 70 per cent reduction in C. finmarchicus biomass, coupled with a gradual northward shift in the species's distribution, which have both been linked with climate change. To determine the potential for C. finmarchicus to track changes in habitat availability and maintain stable effective population sizes, we have assessed levels of gene flow and dispersal in current populations, as well as using a coalescent approach together with palaeodistribution modelling to elucidate the historical population demography of the species over previous changes in Earth's climate. Our findings indicate high levels of dispersal and a constant effective population size over the period 359 000–566 000 BP and suggest that C. finmarchicus possesses the capacity to track changes in available habitat, a feature that may be of crucial importance to the species's ability to cope with the current period of global climate change.

High dispersal potential has maintained long-term population stability in the North Atlantic copepod Calanus finmarchicus

The cool-water copepod Calanus finmarchicus is a key species in North Atlantic marine ecosystems since it represents an important food resource for the developmental stages of several fish of major economic value. Over the last 40 years, however, data from the Continuous Plankton Recorder survey have highlighted a 70 per cent reduction in C. finmarchicus biomass, coupled with a gradual northward shift in the species's distribution, which have both been linked with climate change. To determine the potential for C. finmarchicus to track changes in habitat availability and maintain stable effective population sizes, we have assessed levels of gene flow and dispersal in current populations, as well as using a coalescent approach together with palaeodistribution modelling to elucidate the historical population demography of the species over previous changes in Earth's climate. Our findings indicate high levels of dispersal and a constant effective population size over the period 359 000-566 000 BP and suggest that C. finmarchicus possesses the capacity to track changes in available habitat, a feature that may be of crucial importance to the species's ability to cope with the current period of global climate change.