Genetic diversity, population structure, effective population size and demographic history of the Finnish wolf population

The Finnish wolf population (Canis lupus) was sampled during three different periods (1996-1998, 1999-2001 and 2002-2004), and 118 individuals were genotyped with 10 microsatellite markers. Large genetic variation was found in the population despite a recent demographic bottleneck. No spatial population subdivision was found even though a significant negative relationship between genetic relatedness and geographic distance suggested isolation by distance. Very few individuals did not belong to the local wolf population as determined by assignment analyses, suggesting a low level of immigration in the population. We used the temporal approach and several statistical methods to estimate the variance effective size of the population. All methods gave similar estimates of effective population size, approximately 40 wolves. These estimates were slightly larger than the estimated census size of breeding individuals. A Bayesian model based on Markov chain Monte Carlo simulations indicated strong evidence for a long-term population decline. These results suggest that the contemporary wolf population size is roughly 8% of its historical size, and that the population decline dates back to late 19th century or early 20th century. Despite an increase of over 50% in the census size of the population during the whole study period, there was only weak evidence that the effective population size during the last period was higher than during the first. This may be caused by increased inbreeding, diminished dispersal within the population, and decreased immigration to the population during the last study period.     

Molecular detection of Marteilia sydneyi, pathogen of Sydney rock oysters

The life cycle of Marteilia sydneyi, the aetiological agent of QX disease in the Sydney rock oyster Saccostrea commercialis, is not known. We have developed and optimised 2 diagnostic assays, the polymerase chain reaction (PCR) and in situ hybridisation, for use in investigating the role of possible alternative hosts in the life cycle of this pathogen. PCR primers, designed within the ITS1 rDNA of M. sydneyi, amplified a 195 bp fragment. Sensitivity of the PCR assay was assessed using DNA extracted from known numbers of sporonts purified from infected oyster digestive gland. DNA equivalent to 0.01 sporonts was detectable following agarose gel electrophoresis. The potential inhibitory effect of the presence of host DNA on the PCR assay was tested by the addition of oyster genomic DNA during amplification. Concentrations of host DNA in excess of 50 ng per 20 microliters reaction reduced the sensitivity of the test. Environmental validation of the PCR assay was demonstrated by the amplification of M. sydneyi DNA from 50 ng of genomic DNA extracted from QX-infected oysters. A DNA probe was constructed using the M. sydneyi unique primers and was able to detect 10 pg of M. sydneyi PCR amplified DNA in dot-blot hybridisations. The probe hybridised with presporulating and sporulating M. sydneyi stages in paraffin sections of oyster digestive gland. No non-specific binding was observed. Hybridisation consistency and signal intensity decreased as sporonts matured. While the high sensitivity and specificity of the PCR test will allow rapid screening of large numbers of potential alternative hosts for the presence of parasite DNA, it does not actually identify infective stages. In situ hybridisation conducted on paraffin sections will determine the location of the parasite within the host for morphological characterisation.     

Characterization of phenoloxidase activity in Sydney rock oysters (Saccostrea glomerata)

Phenoloxidase (PO) activity was studied in Sydney rock oysters (Saccostrea glomerata). As in other molluscs, PO was found to exist as a pro-enzyme (proPO) in hemocytes. ProPO could be activated to PO by exogenous proteases (trypsin and chymotrypsin), exposure of hemocytes to pathogen-associated molecular patterns (PAMPs) and by the detergents, Triton X-100 and sodium dodecyl sulphate (SDS). Inhibition studies confirmed the proPO activating system of Sydney rock oysters is a proteinase cascade in which Ca²⁺ dependent serine proteinases proteolytically convert proPO into active PO. Activated PO was found to be a tyrosinase-like enzyme that is responsible for both monophenolase and diphenolase activity. The bifunctional PO had higher affinity for the monophenol, hydroquinine monomethyl ether (4HA) (K_m = 4.45 ± 1.46 mM) than for the diphenol, l-DOPA (K_m = 10.27 ± 1.33 mM). Maximum enzyme activity was evident at 37 °C, pH 8 and at salinities of between 30 and 37 ppt. Melanogenesis catalysed by the active enzyme is a composite of eumelanin and the product of a sclerotin pathway combining DOPA decarboxylase with PO activity.     

Genetic variation and effective population size in isolated populations of coastal cutthroat trout

Following glacial recession in southeast Alaska, waterfalls created by isostatic rebound have isolated numerous replicate populations of coastal cutthroat trout (Oncorhynchus clarkii clarkii) in short coastal streams. These replicate isolated populations offer an unusual opportunity to examine factors associated with the maintenance of genetic diversity. We used eight microsatellites to examine genetic variation within and differentiation among 12 population pairs sampled from above and below these natural migration barriers. Geological evidence indicated that the above-barrier populations have been isolated for 8,000–12,500 years. Genetic differentiation among below-barrier populations (F _ST = 0.10, 95% C.I. 0.08–0.12) was similar to a previous study of more southern populations of this species. Above-barrier populations were highly differentiated from adjacent below-barrier populations (mean pairwise F _ST = 0.28; SD 0.18) and multiple lines of evidence were consistent with asymmetric downstream gene flow that varied among streams. Each above-barrier population had reduced within-population genetic variation when compared to the adjacent below-barrier population. Within-population genetic diversity was significantly correlated with the amount of available habitat in above-barrier sites. Increased genetic differentiation of above-barrier populations with lower genetic diversity suggests that genetic drift has been the primary cause of genetic divergence. Long-term estimates of N _e based on loss of heterozygosity over the time since isolation were large (3,170; range 1,077–7,606) and established an upper limit for N _e if drift were the only evolutionary process responsible for loss of genetic diversity. However, it is likely that a combination of mutation, selection, and gene flow have also contributed to the genetic diversity of above-barrier populations. Contemporary above-barrier N _e estimates were much smaller than long-term N _e estimates, not correlated with within-population genetic diversity, and not consistent with the amount of genetic variation retained, given the approximate 10,000-year period of isolation. The populations isolated by waterfalls in this study that occur in larger stream networks have retained substantial genetic variation, which suggests that the amount of habitat in headwater streams is an important consideration for maintaining the evolutionary potential of isolated populations.     

Genetic variation, genetic structure and effective population size in the tropical holoparasitic endophyteBdallophyton bambusarum (Rafflesiaceae)

The genetic population structure inBdallophyton bambusarum, an endoparasite, was studied in ten subpopulations from a subdeciduous tropical forest in Veracruz Mexico. The sample was analyzed using seven polymorphic loci in cellulose acetate electrophoresis. Isozyme data indicated that the subpopulations ofB. bambusarum contained high genetic variability (H_ep = 0.452 ± 0.045, S.E.). Our analysis suggests that almost each inflorescence ofB. bambusarum is an individual. The subpopulations studied were genetically similar (average Nei's genetic identity 0.941 ± 0.051 and F _st values 0.097 ± 0.026), suggesting that genetic differentiation among subpopulations was small. Direct estimates of effective population size was derived from observations of three fluorescent dyes, and from the genetic neighborhood area derived from these data. The neighborhood area, multiplied by the total density of individuals, gave an N_e = 124.84 plants, and when corrected to consider the proportion of males and females gave an N_e = 118.59 individuals. An indirect estimate of Nm was obtained from the F _st values (mean Nm=2.037), giving an indirect estimate of the effective population size N_b = 12.8 individuals. Both values are relatively high when compared to other plant studies. The gene flow and/or effective populations size of the studied subpopulations ofB. bambusarum are believed to be large enough to prevent differentiation among subpopulations due to genetic drift.     

The effective population size of Anopheles gambiae in Kenya: implications for population structure

We estimated current and long-term effective population size (Ne) of two Anopheles gambiae (savanna cytotype) populations in Kenya. Temporal variation at nine microsatellite loci in each population sampled 7 and 9 years apart and genetic diversity in each sample were analyzed to answer the following questions. (1) Do bottlenecks occur in Kenyan populations of A. gambiae? (2) How variable are different populations with respect to their current and long-term Ne values? (3) What are the implications of these results on population structure and history? The estimates of Ne of Asembo and Jego were 6,359 and 4,258, respectively, and the lower 95% limits were 2,455 and 1,669, respectively. Thus, despite the typical observation of low density at the village level during the dry season, large populations are maintained annually. Large current Ne is consistent with previous studies showing low differentiation across the continent, especially under Wright's isolation-by-distance model. Current Ne in Asembo was 1.5-fold higher than in Jego, but this difference was not significant. Long-term Ne in Asembo (22,667) was 2.9-fold higher than that in Jego (7,855) based on the stepwise mutation model. The difference between populations was significant at both time points regardless of whether long-term Ne values were calculated based on the stepwise mutation model or the infinite-alleles model. Heterozygosity in Jego declined significantly between 1987 (59%) and 1996 (54%), whereas heterozygosity in Asembo was stable (66%-65%). Despite the relatively high and significant differentiation between Asembo and Jego (FST = 0.072-0.10, RST = 0.037- 0.038), all alleles in Jego were found in Asembo but not vice versa. All of these findings suggest that lower Ne in Jego magnifies differentiation between the two populations. The long-term Ne was biased downward, because its calculation was based on an upper bound estimate of microsatellite mutation rate. Ne values based on mtDNA and allozymes were an order of magnitude higher. Long-term Ne therefore, is probably measured in hundreds of thousands and hence does not support a recent expansion of this species from a small population.      

Genetic population structure and neighbourhood population size estimates of the caddisfly Plectrocnemia conspersa

1. We used both genetic and ecological methods to evaluate the role of history and the scale of colonisation in structuring populations of the caddisfly Plectrocnemia conspersa. There was no genetic differentiation between sites up to 20 km apart, despite population sizes suggesting that genetic drift could create substantial variation at this scale. 2. Genetic differentiation between populations separated by more than 20 km was greater than expected given the contrasting short‐range trend, and implied a neighbourhood population size that is implausibly small. Therefore, the evolutionary processes that affect the short‐range trend do not explain differentiation over greater distances. 3. At small scales (<20 km), relatively short flights by winged adults spread over a number of generations could account for the spread of genes. For instance, dispersing individuals could found small (often temporary) populations, which may then grow and exchange genes with larger and more permanent local populations, amplifying the effects of the initial gene flow. 4. Over larger scales (20–500 km), substantial gaps between regions containing suitable habitat patches could reduce the number of colonisation events. Genetic patterns at this scale may date from the time they were last colonised. Previous ecological studies have rarely examined the dynamics of aquatic insect populations over these larger geographical scales, yet these processes may be central to their persistence and spread.     

Induction of phenoloxidase and other immunological activities in Sydney rock oysters challenged with microbial pathogen-associate molecular patterns

This study investigates the effects of two pathogen-associated molecular patterns (PAMPs), LPS and zymosan, on the Sydney rock oyster (Saccostrea glomerata) immune system. Phenoloxidase and phagocytic activities, total and differential haemocyte frequencies, as well as peroxide and superoxide concentrations were measured after the injection of lipopolysaccharide and zymosan. All of the immunological parameters were induced by both PAMPs. Phenoloxidase (monophenolase and diphenolase) and phagocytic activities, as well as the frequencies of phenoloxidase-positive haemocytes, hyalinocytes and granulocytes in the haemolymph, increased within 24 h of PAMP injection. Values for all of these parameters peaked within 48 h of challenge and began to decrease to levels that were indistinguishable from those of controls within 96h. The only exception to this pattern was diphenolase activity, which remained elevated for at least 96 h. Control saline injections that lacked PAMPs also induced responses in most of the parameters measured. However, reactions to saline injections were of far lower magnitude compared to those induced by PAMPs. All of the data suggest that the phenoloxidase and phagocytic systems of oysters are inducible components of the Sydney rock oyster immune system, and that induction is primarily due to increased frequencies of specialised haemocytes in the haemolymph.     

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.     

Genetic diversity,effective population size,and structure among black bear populations in the Lower Mississippi Alluvial Valley,USA

Multiple small populations of American black bears Ursus americanus, including the recently delisted Louisiana black bear subspecies U. a. luteolus, occupy a fragmented landscape in the Lower Mississippi Alluvial Valley, USA (LMAV). Populations include bears native to the LMAV, bears translocated from Minnesota during the 1960s, and recently reintroduced and colonizing populations sourced from within the LMAV. We estimated population structure, gene flow, and genetic parameters important to conservation of small populations using genotypes at 23 microsatellite markers for 265 bears from seven populations. We inferred five genetic clusters corresponding to the following populations: White River and western Mississippi, Tensas River and Three Rivers, Upper Atchafalaya, Lower Atchafalaya, and Minnesota. Upper Atchafalaya was suggested as the product of Minnesota-sourced translocations, but those populations have since diverged, likely because of a founder effect followed by genetic drift and isolation. An admixture zone recently developed in northeastern Louisiana and western Mississippi between migrants from White River and Tensas River, resulting in a Wahlund effect. However, gene flow among most populations has been limited and considerable genetic differentiation accumulated (global F_ST?=?0.22), particularly among the three Louisiana black bear populations that existed when federal listing occurred. Consistent with previous bottlenecks, founder effects, and persisting isolation, all LMAV bear populations had low genetic diversity (A_R?=?2.08–4.81; H_E?=?0.36–0.63) or small effective population size (N_E?=?3–49). Translocating bears among populations as part of a regional genetic restoration program may help improve genetic diversity and increase effective population sizes.     

Molecular characterisation of a haplosporidian parasite infecting rock oysters Saccostrea cuccullata in north Western Australia

A haplosporidian parasite was identified in rock oysters (Saccostrea cuccullata Born, 1778) from the Montebello Islands (latitude -20.4'S longitude 115.53'E) off the northern coast of Western Australia by histopathological examination, PCR amplification and DNA sequencing of a segment of the SSU region of the parasite's rRNA gene. An oligonucleotide probe was constructed from the parasite's SSU rRNA gene in order to confirm its presence by in situ hybridisation. The parasite was disseminated throughout the gonad follicles of the host and to a lesser extent in the gills. The only parasite life stages thus far observed in this study were a uninucleate naked cell assumed to be a precursor to multinucleate plasmodial stages and a binucleate plasmodial stage. Whilst no parasite spores were detected in affected rock oysters, a phylogenetic analysis of the SSU region of the parasite's rRNA gene indicates the parasite belongs to the genus Minchinia. A PCR and in situ hybridisation assay for the Minchinia sp. was used to identify haplosporidians described by Hine and Thorne [Hine, P.M.., Thorne, T., 2002. Haplosporidium sp. (Haplosporidia: Haplosporidiidae) associated with mortalities among rock oysters Saccostrea cuccullata in north Western Australia. Dis. Aquat. Organ. 51, 123-13], in archived rock oyster tissues from the same coastline.     

Coalescent process with fluctuating population size and its effective size

We consider a Wright-Fisher model whose population size is a finite Markov chain. We introduce a sequence of two-dimensional discrete time Markov chains whose components describe the coalescent process and the fluctuation of population size. For the limiting process of the sequence of Markov chains, the relationship of the expectation of coalescence time to the harmonic and the arithmetic means of population sizes is shown, and the Laplace transform of the distribution of coalescence time is calculated. We define the coalescence effective population size (cEPS) by the expectation of coalescence time. We show that cEPS is strictly larger (resp. smaller) than the harmonic (resp. arithmetic) mean. As the population size fluctuates more quickly (resp. slowly), cEPS is closer to the harmonic (resp. arithmetic) mean. For the case of a two-valued Markov chain, we show the explicit expression of cEPS and its dependency on the sample size.     

Testing demographic models of effective population size

Due to practical difficulties in obtaining direct genetic estimates of effective sizes, conservation biologists have to rely on so-called 'demographic models' which combine life-history and mating-system parameters with F-statistics in order to produce indirect estimates of effective sizes. However, for the same practical reasons that prevent direct genetic estimates, the accuracy of demographic models is difficult to evaluate. Here we use individual-based, genetically explicit computer simulations in order to investigate the accuracy of two such demographic models aimed at investigating the hierarchical structure of populations. We show that, by and large, these models provide good estimates under a wide range of mating systems and dispersal patterns. However, one of the models should be avoided whenever the focal species' breeding system approaches monogamy with no sex bias in dispersal or when a substructure within social groups is suspected because effective sizes may then be strongly overestimated. The timing during the life cycle at which F-statistics are evaluated is also of crucial importance and attention should be paid to it when designing field sampling since different demographic models assume different timings. Our study shows that individual-based, genetically explicit models provide a promising way of evaluating the accuracy of demographic models of effective size and delineate their field of applicability.     

Generation time and effective population size in Polar Eskimos

North Greenland Polar Eskimos are the only hunter-gatherer population, to our knowledge, who can offer precise genealogical records spanning several generations. This is the first report from Eskimos on two key parameters in population genetics, namely, generation time (T) and effective population size (Ne). The average mother-daughter and father-son intervals were 27 and 32 years, respectively, roughly similar to the previously published generation times obtained from recent agricultural societies across the world. To gain an insight for the generation time in our distant ancestors, we calculated maternal generation time for two wild chimpanzee populations. We also provide the first comparison among three distinct approaches (genealogy, variance and life table methods) for calculating Ne, which resulted in slightly differing values for the Eskimos. The ratio of the effective to the census population size is estimated as 0.6-0.7 for autosomal and X-chromosomal DNA, 0.7-0.9 for mitochondrial DNA and 0.5 for Y-chromosomal DNA. A simulation of alleles along the genealogy suggested that Y-chromosomal DNA may drift a little faster than mitochondrial DNA in this population, in contrast to agricultural Icelanders. Our values will be useful not only in prehistoric population inference but also in understanding the shaping of our genome today.     

Genetic variation in a semi-natural Drosophila population after a bottleneck I. Lethals,their allelism and effective population size

A semi-natural Drosophila melanogaster population was twice forced through a genetic bottleneck and allowed to recover naturally. In one case additional variation was introduced to the recovering population. The percentage of lethal chromosomes, the level of allelism between these lethals, and the effective population size calculated from the allelism of these lethals all rose sharply in the few generations following each bottleneck, though this was not the case in the very first generation. Thereafter this rise decelerated rapidly and never returned to pre-bottleneck levels. Additional introduced variation had little effect. The reasons for and implications of this pattern have been considered.     

Changing effective population size and the McDonald-Kreitman test

Artifactual evidence of adaptive amino acid substitution can be generated within a McDonald-Kreitman test if some amino acid mutations are slightly deleterious and there has been an increase in effective population size. Here I investigate the conditions under which this occurs. I show that fairly small increases in effective population size can generate artifactual evidence of positive selection if there is no selection upon synonymous codon use. This problem is exacerbated by the removal of low-frequency polymorphisms. However, selection on synonymous codon use restricts the conditions under which artifactual evidence of adaptive evolution is produced.