Absence of population genetic differentiation in the New Zealand greenshell mussel Perna canaliculus (Gmelin 1791) as assessed by allozyme variation

Genetic variation in the endemic New Zealand greenshell mussel, Perna canaliculus (Gmelin 1791), was examined using starch-gel electrophoresis at seven protein-coding loci (Idh; Acon-1; Acon-2; Gpd; Pgi; Pgm; Pgd) in 35 populations (N=1038 mussels). For all loci and all populations, Fisher's exact tests indicated highly significant departures from Hardy-Weinberg equilibrium (HWE), but this overall result was caused by significant heterozygote deficiencies at only two loci (Pgm and Pgd), and in only three northern populations (Kuaotunu, Te Haumi and Days Bay). Allelic and genotypic differentiation between population pairs at individual loci and across all loci were nonsignificant, and genotypic disequilibrium at each locus pair was also nonsignificant for all populations. Genetic variation in all populations was high (mean heterozygosity, 0.210+/-0.113), while Nei's D among populations was very low (0.002+/-0.002). Low population subdivision (θ=-0.001+/-0.002) and high levels of gene flow (Nm(p)=10.18; Nm(θ)=infinity) also indicated that the single panmictic unit model best explains population genetic homogeneity in P. canaliculus over a north-south distance >2000 km. Lack of genetic subdivision in this species is discussed in light of two previous allozyme studies, with differing results: one suggested that a north-south division exists between greenshell mussel stocks, and the other suggested that population structure in this species can be explained through isolation by distance model modified by local hydrology.     

The influence of age structure and fecundity on effective population size.

Simple formulae are developed which define the effective size (Ne) of populations with overlapping generations, and their use is illustrated using data from a squirrel population. Two mating systems are considered, the random union of gametes and monogamy, in combination with age-independent fecundity. In the simplest case of age-independent (type 2) survivorship in a population of N adults, Ne = N/(2-T-1) where T is the generation time. As T increases, Ne declines asymptomatically to N/2. A generalization of this result (Ne = N/[1 + k-1-T-1], where k influences survivorship) shows that given type 1 survivorship (k greater than 1) this decline in Ne is less severe. A biased sex ratio results in Ne differing between the two mating systems; however, in both systems, a sex ratio bias resulting from survival differences has much less influence on Ne than a sex ratio bias resulting from recruitment differences. Low fecundity can increase Ne, but realistic levels of variation among breeding individuals (Poisson or greater) negate the effect. The effect on Ne of variation resulting from the presence of non-breeders is also considered.     

Estimating effective population size or mutation rate with microsatellites

Microsatellites are short tandem repeats that are widely dispersed among eukaryotic genomes. Many of them are highly polymorphic; they have been used widely in genetic studies. Statistical properties of all measures of genetic variation at microsatellites critically depend upon the composite parameter theta = 4Nmicro, where N is the effective population size and micro is mutation rate per locus per generation. Since mutation leads to expansion or contraction of a repeat number in a stepwise fashion, the stepwise mutation model has been widely used to study the dynamics of these loci. We developed an estimator of theta, theta; (F), on the basis of sample homozygosity under the single-step stepwise mutation model. The estimator is unbiased and is much more efficient than the variance-based estimator under the single-step stepwise mutation model. It also has smaller bias and mean square error (MSE) than the variance-based estimator when the mutation follows the multistep generalized stepwise mutation model. Compared with the maximum-likelihood estimator theta; (L) by, theta; (F) has less bias and smaller MSE in general. theta; (L) has a slight advantage when theta is small, but in such a situation the bias in theta; (L) may be more of a concern.     

Weak population regulation in ecological time series

How strongly natural populations are regulated has a long history of debate in ecology. Here, we discuss concepts of population regulation appropriate for stochastic population dynamics. We then analyse two large collections of data sets with autoregressive-moving average (ARMA) models, using model selection techniques to find best-fitting models. We estimated two metrics of population regulation: the characteristic return rate of populations to stationarity and the variability of the stationary distribution (the long-term distribution of population abundance). Empirically, longer time series were more likely to show weakly regulated population dynamics. For data sets of length ≥ 20, more than 35% had characteristic return times > 6 years, and more than 29% had stationary distributions whose coefficients of variation were more than two times greater than would be the case if they were maximally regulated. These results suggest that many natural populations are weakly regulated.
Ecology Letters (2010) 13: 21–31     

Effects of Mutation on Selection Limits in Finite Populations with Multiple Alleles

The ultimate response to directional selection (i.e., the selection limit) under recurrent mutation is analyzed by a diffusion approximation for a population in which there are k possible alleles at a locus. The limit mainly depends on two scaled parameters S (= 4Ns sigma a) and theta (= 4Nu) and k, the number of alleles, where N is the effective population size, u is the mutation rate, s is the selection coefficient, and sigma 2a is the variance of allelic effects. When the selection pressure is weak (S less than or equal to 0.5), the limit is given approximately by 2S sigma a[1 - (1 + c2)/k]/(theta + 1) for additive effects of alleles, where c is the coefficient of variation of the mutation rates among alleles. For strong selection, other approximations are devised to analyze the limit in different parameter regions. The effect of mutation on selection limits largely relies on the potential of mutation to introduce new and better alleles into the population. This effect is, however, bounded under the present model. Unequal mutation rates among alleles tend to reduce the selection limit, and can have a substantial effect only for small numbers of alleles and weak selection. The selection limit decreases as the mutation rate increases.     

Histone Variant HTZ1 Shows Extensive Epistasis with,but Does Not Increase Robustness to,New Mutations

Biological systems produce phenotypes that appear to be robust to perturbation by mutations and environmental variation. Prior studies identified genes that, when impaired, reveal previously cryptic genetic variation. This result is typically interpreted as evidence that the disrupted gene normally increases robustness to mutations, as such robustness would allow cryptic variants to accumulate. However, revelation of cryptic genetic variation is not necessarily evidence that a mutationally robust state has been made less robust. Demonstrating a difference in robustness requires comparing the ability of each state (with the gene perturbed or intact) to suppress the effects of new mutations. Previous studies used strains in which the existing genetic variation had been filtered by selection. Here, we use mutation accumulation (MA) lines that have experienced minimal selection, to test the ability of histone H2A.Z (HTZ1) to increase robustness to mutations in the yeast Saccharomyces cerevisiae. HTZ1, a regulator of chromatin structure and gene expression, represents a class of genes implicated in mutational robustness. It had previously been shown to increase robustness of yeast cell morphology to fluctuations in the external or internal microenvironment. We measured morphological variation within and among 79 MA lines with and without HTZ1. Analysis of within-line variation confirms that HTZ1 increases microenvironmental robustness. Analysis of between-line variation shows the morphological effects of eliminating HTZ1 to be highly dependent on the line, which implies that HTZ1 interacts with mutations that have accumulated in the lines. However, lines without HTZ1 are, as a group, not more phenotypically diverse than lines with HTZ1 present. The presence of HTZ1, therefore, does not confer greater robustness to mutations than its absence. Our results provide experimental evidence that revelation of cryptic genetic variation cannot be assumed to be caused by loss of robustness, and therefore force reevaluation of prior claims based on that assumption.     

Maintenance of genetic differentiation across a transition zone in the sea: discordance between nuclear and cytoplasmic markers

To investigate the origin and maintenance of the genetic discontinuity between Atlantic and Mediterranean populations of the common sea bass (Dicentrarchus labrax) we analysed the genetic variation at a fragment of mitochondrial cytochrome b sequence for 18 population samples. The result were also compared with new or previously published microsatellite data. Seven mitochondrial haplotypes and an average nucleotidic divergence of 0.02 between Atlantic and Mediterranean populations that matches a Pleistocene allopatric isolation were found. The frequency variation at the cytochrome b locus was many times greater between Atlantic and Mediterranean populations (theta(C) = 0.67) than at microsatellite loci (theta(N)= 0.02). The examination of the different evolutionary forces at play suggests that a sex-biased hybrid breakdown is a likely explanation for part of the observed discrepancy between mitochondrial and nuclear loci. In addition, an analysis is made of the correlation between microsatellite loci points towards the possible existence of a hybrid zone in samples from the Alboran Sea.     

Population Growth Inflates the Per-Individual Number of Deleterious Mutations and Reduces Their Mean Effect

This study addresses the question of how purifying selection operates during recent rapid population growth such as has been experienced by human populations. This is not a straightforward problem because the human population is not at equilibrium: population genetics predicts that, on the one hand, the efficacy of natural selection increases as population size increases, eliminating ever more weakly deleterious variants; on the other hand, a larger number of deleterious mutations will be introduced into the population and will be more likely to increase in their number of copies as the population grows. To understand how patterns of human genetic variation have been shaped by the interaction of natural selection and population growth, we examined the trajectories of mutations with varying selection coefficients, using computer simulations. We observed that while population growth dramatically increases the number of deleterious segregating sites in the population, it only mildly increases the number carried by each individual. Our simulations also show an increased efficacy of natural selection, reflected in a higher fraction of deleterious mutations eliminated at each generation and a more efficient elimination of the most deleterious ones. As a consequence, while each individual carries a larger number of deleterious alleles than expected in the absence of growth, the average selection coefficient of each segregating allele is less deleterious. Combined, our results suggest that the genetic risk of complex diseases in growing populations might be distributed across a larger number of more weakly deleterious rare variants.     

The loss of adaptive plasticity during long periods of environmental stasis

Adaptive plasticity allows populations to adjust rapidly to environmental change. If this is useful only rarely, plasticity may undergo mutational degradation and be lost from a population. We consider a population of constant size N undergoing loss of plasticity at functional mutation rate m and with selective advantage s associated with loss. Environmental change events occur at rate theta per generation, killing all individuals that lack plasticity. The expected time until loss of plasticity in a fluctuating environment is always at least tau, the expected time until loss of plasticity in a static environment. When mN > 1 and N theta > 1, we find that plasticity will be maintained for an average of at least 10(8) generations in a single population, provided tau > 18/theta. In a metapopulation, plasticity is retained under the more lenient condition tau > 1.3/theta, irrespective of mN, for a modest number of demes. We calculate both exact and approximate solutions for tau and find that it is linearly dependent only on the logarithm of N, and so, surprisingly, both the population size and the number of demes in the metapopulation make little difference to the retention of plasticity. Instead, tau is dominated by the term 1/(m+s/2).     

Reproduction numbers and thresholds in stochastic epidemic models. I. Homogeneous populations.

We compare threshold results for the deterministic and stochastic versions of the homogeneous SI model with recruitment, death due to the disease, a background death rate, and transmission rate beta cXY/N. If an infective is introduced into a population of susceptibles, the basic reproduction number, R0, plays a fundamental role for both, though the threshold results differ somewhat. For the deterministic model, no epidemic can occur if R0 less than or equal to 1 and an epidemic occurs if R0 greater than 1. For the stochastic model we find that on average, no epidemic will occur if R0 less than or equal to 1. If R0 greater than 1, there is a finite probability, but less than 1, that an epidemic will develop and eventuate in an endemic quasi-equilibrium. However, there is also a finite probability of extinction of the infection, and the probability of extinction decreases as R0 increases above 1.     

Variance of Neutral Genetic Variances within and between Populations for a Quantitative Character

The variances of genetic variances within and between finite populations were systematically studied using a general multiple allele model with mutation in terms of identity by descent measures. We partitioned the genetic variances into components corresponding to genetic variances and covariances within and between loci. We also analyzed the sampling variance. Both transient and equilibrium results were derived exactly and the results can be used in diverse applications. For the genetic variance within populations, sigma 2 omega, the coefficient of variation can be very well approximated as [formula: see text] for a normal distribution of allelic effects, ignoring recurrent mutation in the absence of linkage, where m is the number of loci, N is the effective population size, theta 1(0) is the initial identity by descent measure of two genes within populations and t is the generation number. The first term is due to genic variance, the second due to linkage disequilibrium, and third due to sampling. In the short term, the variation is predominantly due to linkage disequilibrium and sampling; but in the long term it can be largely due to genic variance. At equilibrium with mutation [formula: see text] where u is the mutation rate. The genetic variance between populations is a parameter. Variance arises only among sample estimates due to finite sampling of populations and individuals. The coefficient of variation for sample gentic variance between populations, sigma 2b, can be generally approximated as [formula: see text] when the number of loci is large where S is the number of sampling populations.     

Estimating sex-specific processes in human populations: Are XY-homologous markers an effective tool?

Homologous markers on the sex-specific regions of the X- and Y-chromosomes are differentially inherited through males and females, and have similar molecular characteristics. They may therefore be useful as a complement to the comparison of mtDNA and Y-chromosomal haplotypes for estimating sex-specific processes shaping human population structure. To test this idea, we analyzed XY-homologous microsatellite diversity in 33 human populations from Africa, Asia and Europe. Interpopulation comparisons suggest that the generally discordant pattern of genetic variation observed for X- and Y-linked markers could be an outcome of sex-specific migration processes (m(females)/m(males) approximately 3) or sex-specific demographic processes (N(females)/N(males) approximately 11) or a combination of both. However, intrapopulation diversity estimated by the X/Y ratio Watterson estimator (theta(H(Y))/theta(H(X))) suggests that the scenarios required to explain the global genetic variation of XY-homologous markers are many and complex, and that the sex-specific processes (effective population size and migration rate) shaping human population structures are likely to be specific to each population under study. XY-homologous markers provide an insight into the genuine complexity of sex-specific processes, and their further exploitation in human population studies seems worthwhile.     

Selection intensity against deleterious mutations in RNA secondary structures and rate of compensatory nucleotide substitutions

A two-locus model of reversible mutations with compensatory fitness interactions is presented; single mutations are assumed to be deleterious but neutral in appropriate combinations. The expectation of the time of compensatory nucleotide substitutions is calculated analytically for the case of tight linkage between sites. It is shown that selection increases the substitution time dramatically when selection intensity Ns > 1, where N is the diploid population size and s the selection coefficient. Computer simulations demonstrate that recombination increases the substitution time, but the effect of recombination is small when selection is weak. The amount of linkage disequilibrium generated in the process of compensatory substitution is also investigated. It is shown that significant linkage disequilibrium is expected to be rare in natural populations. The model is applied to the mRNA secondary structure of the bicoid 3' untranslated region of Drosophila. It is concluded that average selection intensity Ns against single deleterious mutations is not likely to be much larger than 1.     

Genome-wide selection on codon usage at the population level in the fungal model organism Neurospora crassa

Many organisms exhibit biased codon usage in their genome, including the fungal model organism Neurospora crassa. The preferential use of subset of synonymous codons (optimal codons) at the macroevolutionary level is believed to result from a history of selection to promote translational efficiency. At present, few data are available about selection on optimal codons at the microevolutionary scale, that is, at the population level. Herein, we conducted a large-scale assessment of codon mutations at biallelic sites, spanning more than 5,100 genes, in 2 distinct populations of N. crassa: the Caribbean and Louisiana populations. Based on analysis of the frequency spectra of synonymous codon mutations at biallelic sites, we found that derived (nonancestral) optimal codon mutations segregate at a higher frequency than derived nonoptimal codon mutations in each population; this is consistent with natural selection favoring optimal codons. We also report that optimal codon variants were less frequent in longer genes and that the fixation of optimal codons was reduced in rapidly evolving long genes/proteins, trends suggestive of genetic hitchhiking (Hill-Robertson) altering codon usage variation. Notably, nonsynonymous codon mutations segregated at a lower frequency than synonymous nonoptimal codon mutations (which impair translational efficiency) in each N. crassa population, suggesting that changes in protein composition are more detrimental to fitness than mutations altering translation. Overall, the present data demonstrate that selection, and partly genetic interference, shapes codon variation across the genome in N. crassa populations.     

Incorporating spatial variation in density enhances the stability of simple population dynamics models

Simple discrete time models of population growth admit a wide variety of dynamic behaviors, including population cycles and chaos. Yet studies of natural and laboratory populations typically reveal their dynamics to be relatively stable. Many explanations for the apparent rarity of unstable or chaotic behavior in real populations have been developed, including the possible stabilizing roles of migration, refugia, abrupt density-dependence, and genetic variation in sensitivity to density. We develop a theoretical framework for incorporating random spatial variation in density into simple models of population growth, and apply this approach to two commonly used models in ecology: the Ricker and Hassell maps. We show that the incorporation of spatial density variation into both these models has a strong stabilizing influence on their dynamic behavior, and leads to their exhibiting stable point equilibria or stable limit cycles over a relatively much larger range of parameter values. We suggest that one reason why chaotic population dynamics are less common than the simple models indicate is, these models typically neglect the potentially stabilizing role of spatial variation in density.     

Microsatellite analysis of genetic diversity and population genetic structure of a wild rice (<Emphasis Type="Italic">Oryza rufipogon</Emphasis> Griff.) in China

Genetic diversity and population genetic structure of natural Oryza rufipogon populations in China were studied based on ten microsatellite loci. For a total of 237 individuals of 12 populations collected from four regions, a moderate to high level of genetic diversity was observed at population levels with the number of alleles per locus ( A) ranging from 2 to 18 (average 10.6), and polymorphic loci ( P) from 40.0% to 100% (average 83.3%). The observed heterozygosity ( H(O)) varied from 0.163 to 0.550 with the mean of 0.332, and the expected heterozygosity ( H(E)) from 0.164 to 0.648 with the mean of 0.413. The level of genetic diversity for Guangxi was the highest. These results are in good agreement with previous allozyme and RAPD studies. However, it was unexpected that high genetic differentiation among populations was found ( R(ST) = 0.5199, theta = 0.491), suggesting that about one-half of the genetic variation existed between the populations. Differentiation (pairwise theta) was positively correlated with geographical distance ( r = 0.464), as expected under the isolation by distance model. The habitat destruction and degradation throughout the geographic range of O. rufipogon may be the main factor attributed to high genetic differentiation among populations of O. rufipogon in China.     

Simple method for analyzing the pattern of DNA polymorphism and its application to SNP data of human

In order to analyze the pattern of DNA polymorphism in detail, we have developed a simple method using a new statistic theta(i) which estimates 4Nmu from the number of segregating sites whose allelic nucleotide frequency is i/n among n DNA sequences, where N is the effective population size and mu is the mutation rate per generation per nucleotide site. Under the assumption that mutations are selectively neutral and a population size is constant, the expectation of theta(i) is equal to that of theta, which estimates 4Nmu from the number of segregating sites, so that the distribution of theta(i) is flat. Therefore, the departure of the distribution of theta(i) from the horizontal line, which represents the value of theta, reflects change in population size and natural selection. Results of the coalescent simulation show that the distributions of theta(i) in the populations which experienced expansion and reduction are U-shaped and upside-down U-shaped, respectively. And the distributions of theta(i) in some populations that experienced bottleneck are W-shaped. Furthermore, we have applied this method to the SNP data in the International HapMap Project. Results of data analyses show that the distributions of theta(i) in the CEU (European), CHB and JPT (Asian) populations are different from that in the YRI population (African). From these results of data analyses in nuclear DNA and the pattern of polymorphism in human mitochondrial DNA already known, we infer that the CEU, CHB and JPT populations experienced the bottleneck.     

High genetic diversity in Sarracenia leucophylla(Sarraceniaceae), a carnivorous wetland herb

Eighteen allozyme loci were used to examine genetic diversity in 10 natural populations of Sarracenia leucophylla Raf., a pitcher plant restricted to the southeastern United States. One ex situ population propagated for restoration in Georgia was also analyzed. S. leucophylla is an insect-pollinated, outcrossing perennial wetland herb that is threatened over much of its geographic range. Fifteen loci (83.3%) were polymorphic, with a mean number of alleles of 3.33. Compared to species having similar life-history traits and to previously analyzed Sarracenia species, S. leucophylla displayed unexpectedly high genetic diversity. For example, genetic diversity within the species (Hes) was 0.224 and mean population genetic diversity (Hep) was 0.183. Although small S. leucophylla populations maintained less genetic diversity than larger ones, these differences were not statistically significant. Nonetheless, this suggests that small populations may have lost rare alleles. Statistically significant genetic differentiation among populations was found (theta = 0.192, P < .01), although it was not atypical considering the species' life-history characteristics. A significant correlation (P < .01) between genetic and geographic distance was found, indicating an isolation-by-distance effect. However, the correlation coefficient for this relationship was low (r = 0.46), suggesting that factors other than gene flow play a prominent role in the geographic distribution of genetic diversity within the species. The ex situ population captured most of the allozyme variation found in its source population.     

The rate of antigenic variation in fly-transmitted and syringe-passaged infections of Trypanosoma brucei

Rates of antigenic variation were measured in vivo in several populations of cloned lines of Trypanosoma brucei before and after cyclical transmission through tsetse flies. Two cloned lines were adapted for use in laboratory conditions by extensive syringe passaging and rates of antigenic switching/cell/generation were less than 3×10⁻⁶ and 1×10⁻⁴ in each line. Rates of switching were then determined after fly transmission of the first line and generated per capita rate values of greater than 2×10⁻³ in three of four populations examined. In the fourth population the switch rate was lower: less than 7×10⁻⁵ switches/cell/generation. These data show that rates of antigenic variation are several orders of magnitude lower in syringe-passaged lines, such as those routinely used in the majority of laboratory studies, compared with most recently fly-transmitted lines. They also show that the reduction in switching rate associated with syringe passaging is reversible and is thus unlikely to be caused by mutation.