Toll-like receptor variation in the bottlenecked population of the Seychelles warbler: computer simulations see the 'ghost of selection past' and quantify the 'drift debt'.

Balancing selection can maintain immunogenetic variation within host populations, but detecting its signal in a postbottlenecked population is challenging due to the potentially overriding effects of drift. Toll-like receptor genes (TLRs) play a fundamental role in vertebrate immune defence and are predicted to be under balancing selection. We previously characterized variation at TLR loci in the Seychelles warbler (Acrocephalus sechellensis), an endemic passerine that has undergone a historical bottleneck. Five of seven TLR loci were polymorphic, which is in sharp contrast to the low genomewide variation observed. However, standard population genetic statistical methods failed to detect a contemporary signature of selection at any TLR locus. We examined whether the observed TLR polymorphism could be explained by neutral evolution, simulating the population's demography in the software DIYABC. This showed that the posterior distributions of mutation rates had to be unrealistically high to explain the observed genetic variation. We then conducted simulations with an agent-based model using typical values for the mutation rate, which indicated that weak balancing selection has acted on the three TLR genes. The model was able to detect evidence of past selection elevating TLR polymorphism in the prebottleneck populations, but was unable to discern any effects of balancing selection in the contemporary population. Our results show drift is the overriding evolutionary force that has shaped TLR variation in the contemporary Seychelles warbler population, and the observed TLR polymorphisms might be merely the 'ghost of selection past'. Forecast models predict immunogenetic variation in this species will continue to be eroded in the absence of contemporary balancing selection. Such 'drift debt' occurs when a gene pool has not yet reached its new equilibrium level of polymorphism, and this loss could be an important threat to many recently bottlenecked populations.

Pace of life, predators and parasites: predator-induced life-history evolution in Trinidadian guppies predicts decrease in parasite tolerance.

A common evolutionary response to predation pressure is increased investment in reproduction, ultimately resulting in a fast life history. Theory and comparative studies suggest that short-lived organisms invest less in defence against parasites than those that are longer lived (the pace of life hypothesis). Combining these tenets of evolutionary theory leads to the specific, untested prediction that within species, populations experiencing higher predation pressure invest less in defence against parasites. The Trinidadian guppy, Poecilia reticulata, presents an excellent opportunity to test this prediction: guppy populations in lower courses of rivers experience higher predation pressure, and as a consequence have evolved faster life histories, than those in upper courses. Data from a large-scale field survey showed that fish infected with Gyrodactylus parasites were of a lower body condition (quantified using the scaled mass index) than uninfected fish, but only in lower course populations. Although the evidence we present is correlational, it suggests that upper course guppies sustain lower fitness costs of infection, i.e. are more tolerant, than lower course guppies. The data are therefore consistent with the pace of life hypothesis of parasite defence allocation, and suggest that life-history traits mediate the indirect effect of predators on the parasites of their prey.

Divergent selection for opsin gene variation in guppy (Poecilia reticulata) populations of Trinidad and Tobago.

The guppy is known to exhibit remarkable interindividual variations in spectral sensitivity of middle to long wavelength-sensitive (M/LWS) cone photoreceptor cells. The guppy has four M/LWS-type opsin genes (LWS-1, LWS-2, LWS-3 and LWS-4) that are considered to be responsible for this sensory variation. However, the allelic variation of the opsin genes, particularly in terms of their absorption spectrum, has not been explored in wild populations. Thus, we examined nucleotide variations in the four M/LWS opsin genes as well as blue-sensitive SWS2-B and ultraviolet-sensitive SWS1 opsin genes for comparison and seven non-opsin nuclear loci as reference genes in 10 guppy populations from various light environments in Trinidad and Tobago. For the first time, we discovered a potential spectral variation (180 Ser/Ala) in LWS-1 that differed at an amino acid site known to affect the absorption spectra of opsins. Based on a coalescent simulation of the nucleotide variation of the reference genes, we showed that the interpopulation genetic differentiation of two opsin genes was significantly larger than the neutral expectation. Furthermore, this genetic differentiation was significantly related to differences in dissolved oxygen (DO) level, and it was not explained by the spatial distance between populations. The DO levels are correlated with eutrophication that possibly affects the color of aquatic environments. These results suggest that the population diversity of opsin genes is significantly driven by natural selection and that the guppy could adapt to various light environments through color vision changes.

Interactions between males guppies facilitates the transmission of the monogenean ectoparasite Gyrodactylus turnbulli.

In a previous study we found that female guppies shoaled more than males and that there was greater transmission of the ectoparasite Gyrodactylus turnbulli between females. Here, to test for a possible sex bias in parasite transmission, we conducted a similar experiment on single sex shoals of male and female guppies, observing host behaviour before and after the introduction of an infected shoal mate. The initial parasite burden was considerably lower in the present experiment (30 worms versus >100 worms previously) and we used a different stock of ornamental guppies (Green Cobra variety versus a Tuxedo hybrid previously). Contrary to our previous finding, males aggregated significantly more than females. Males performed 'sigmoid' displays towards each other, a courtship behaviour that is more generally directed towards females. Due to the high rate of male-male interactions, parasite transmission was 10 times higher between males than between females. Furthermore, shoaling intensity was highest for the most parasitised fish indicating that these infected fish were not avoided by non-parasitised conspecifics. These studies show that certain social behaviours including shoaling and courtship displays, appear to facilitate the transmission of gyrodactylid parasites.

Cryptic MHC polymorphism revealed but not explained by selection on the class IIb peptide-binding region.

The immune genes of the major histocompatibility complex (MHC) are characterized by extraordinarily high levels of nucleotide and haplotype diversity. This variation is maintained by pathogen-mediated balancing selection that is operating on the peptide-binding region (PBR). Several recent studies have found, however, that some populations possess large clusters of alleles that are translated into virtually identical proteins. Here, we address the question of how this nucleotide polymorphism is maintained with little or no functional variation for selection to operate on. We investigate circa 750-850 bp of MHC class II DAB genes in four wild populations of the guppy Poecilia reticulata. By sequencing an extended region, we uncovered 40.9% more sequences (alleles), which would have been missed if we had amplified the exon 2 alone. We found evidence of several gene conversion events that may have homogenized sequence variation. This reduces the visible copy number variation (CNV) and can result in a systematic underestimation of the CNV in studies of the MHC and perhaps other multigene families. We then focus on a single cluster, which comprises 27 (of a total of 66) sequences. These sequences are virtually identical and show no signal of selection. We use microsatellites to reconstruct the populations' demography and employ simulations to examine whether so many similar nucleotide sequences can be maintained in the populations. Simulations show that this variation does not behave neutrally. We propose that selection operates outside the PBR, for example, on linked immune genes or on the "sheltered load" that is thought to be associated to the MHC. Future studies on the MHC would benefit from extending the amplicon size to include polymorphisms outside the exon with the PBR. This may capture otherwise cryptic haplotype variation and CNV, and it may help detect other regions in the MHC that are under selection.

Rapid loss of MHC class II variation in a bottlenecked population is explained by drift and loss of copy number variation.

Population bottlenecks may reduce genetic variation and potentially increase the risk of extinction. Here, we present the first study to use historic samples to analyse loss of variation at the major histocompatibility complex (MHC), which plays a central role in vertebrate disease resistance. Balancing selection acts on the MHC and could moderate the loss of variation expected from drift; however, in a Wisconsin population of greater prairie-chickens (Tympanuchus cupido), the number of MHC class II B alleles per individual declined by 44% following a population bottleneck, compared to a loss of only 8% at microsatellites. Simulations indicate that drift likely reduced MHC variation at the population level, as well as within individuals by reducing the number of gene copies per individual or by fixing the same alleles across multiple loci. These multiple effects of genetic drift on MHC variation could have important implications for immunity and fitness.

Population genetic analysis of microsatellite variation of guppies (Poecilia reticulata) in Trinidad and Tobago: evidence for a dynamic source-sink metapopulation structure, founder events and population bottlenecks.

Riverine fish populations are traditionally considered to be highly structured and subject to strong genetic drift. Here, we use microsatellites to analyse the population structure of the guppy (Poecilia reticulata), focussing on the headwater floodplain area of the Caroni drainage in Trinidad. We also analyse the population genetics of guppies in the Northern Drainage in Trinidad, a habitat characterized by rivers flowing directly into the sea, and a small isolated population in Tobago. Upland Caroni populations are highly differentiated and display low levels of genetic diversity. However, we found no evidence to suggest that these upland populations experienced recent population crashes and the populations appear to approach mutation-drift equilibrium. Dominant downstream migration over both short- and long-time frames has a strong impact on the population genetics of lowland Caroni populations. This drainage system could be considered a source-sink metapopulation, with the tributary furthest downstream representing a 'super sink', receiving immigrants from rivers upstream in the drainage. Moreover, the effective population size in the lowlands is surprisingly low in comparison with the apparently large census population sizes.

Gyro-scope: an individual-based computer model to forecast gyrodactylid infections on fish hosts.

Individual-based computer models (IBM) feature prominently in current theoretical ecology but have only been applied in a small number of parasitological studies. Here we designed an IBM to simulate the infection dynamics of gyrodactylid parasites and immune defence of naïve hosts (i.e. fish previously not exposed to these parasites). We compared the results of the model with empirical data from guppies (Poecilia reticulata) infected with Gyrodactylus parasites. The laboratory experiments on guppies showed that larger fish acquired a heavier parasite load at the peak of the infection. The survival probability declined with increased body size and no fish survived a parasite load of 80 or more worms in this experiment (i.e. lethal load). The model was a good predictor of the Gyrodactylus infection dynamics of guppies and the model output was congruent with previously published data on Gyrodactylus salaris infections of salmon (Salmo salar). Computer simulations indicated that the infections persisted longer on larger hosts and that the parasite load increased exponentially with the body size of the host. Simulations furthermore predicted that the parasite load of fish with a standard length in excess of 17mm (i.e. the size of adult guppies) reached a lethal load. This suggests that in the conditions of the experiment, the immune defence of naïve guppies can offer moderate protection against gyrodactylid infections to juveniles, but not to naïve adult guppies. The model is a useful tool to forecast the development of gyrodactylid infections on single hosts and make predictions about optimal life history strategies of parasites.

The impact of parasites on the life history evolution of guppies (Poecilia reticulata): the effects of host size on parasite virulence.

There is large spatial and temporal variation in the Gyrodactylus parasite fauna across natural guppy (Poecilia reticulata) populations in Trinidad. The life history evolution of these fish could be affected differently in the various habitats depending on the local parasite selection pressure. Here, we experimentally infected three guppy populations with three gyrodactylid strains in the laboratory and monitored the infection by recording the number of parasites and host mortality in a full factorial design. The origin of the guppy population and parasite strain, and the size of the hosts explained significant variation in the survival of hosts. Larger fish carried the highest parasite loads and experienced the highest mortality rates, which suggests that parasite-mediated selection may favour smaller phenotypes, possibly counter-balancing selection pressures by gape-limited predators, mate choice and female fecundity. We observed significant variation in virulence between parasite strains with the captive-bred experimental strain (Gt3) causing the highest mortality of hosts whilst reaching only relatively low maximum burdens. This suggests that adaptations to the captive environment and/or inbreeding depression may alter the virulence of such captive-bred parasites. There were significant differences in survival rate between guppy populations, with infected guppies from the large population of the Lower Aripo River showing a higher survival rate than the fish from the small and genetically less diverse Upper Aripo River population.

Selection by parasites in spate conditions in wild Trinidadian guppies (Poecilia reticulata).

Gyrodactylids are ubiquitous fish parasites and yet, with the notable exception of Gyrodactylus salaris, few studies have reported the effect of these parasites on host survival in natural populations. Here, we assess the impact of the parasite load of gyrodactylids (G. turnbulli and G. bullatarudis) on the survival and migration of guppies (Poecilia reticulata) in their natural habitat of the Aripo River in Trinidad. The recapture rate of males declined by 19% with every additional parasite, a remarkably high figure given that the parasite load in this study ranged from zero to 20 worms. In addition, with an increased number of parasites, males were more prone to be recovered downstream. In contrast, no effect of parasitism was observed in females. The mean parasite load sharply declined after a series of flushing events during heavy seasonal downpours. The parasite load varied significantly between fish depending on their location in the river, and the size of the fish explained variation in parasite load between individuals. The present study indicates that tropical gyrodactylid parasites can play an important role in the ecology of natural fish populations, causing intense bouts of natural selection in guppies during heavy rains in the wet season.

Evolution of MHC class IIB in the genome of wild and ornamental guppies, Poecilia reticulata.

This is the first study to quantify genomic sequence variation of the major histocompatibility complex (MHC) in wild and ornamental guppies, Poecilia reticulata. We sequenced 196-219 bp of exon 2 MHC class IIB (DAB) in 56 wild Trinidadian guppies and 14 ornamental strain guppies. Each of two natural populations possessed high allelic richness (15-16 alleles), whereas only three or fewer DAB alleles were amplified from ornamental guppies. The disparity in allelic richness between wild and ornamental fish cannot be fully explained by fixation of alleles by inbreeding, nor by the presence of non-amplified sequences (ie null alleles). Rather, we suggest that the same allele is fixed at duplicated MHC DAB loci owing to gene conversion. Alternatively, the number of loci in the ornamental strains has contracted during >100 generations in captivity, a hypothesis consistent with the accordion model of MHC evolution. We furthermore analysed the substitution patterns by making pairwise comparisons of sequence variation at the putative peptide binding region (PBR). The rate of non-synonymous substitutions (dN) only marginally exceeded synonymous substitutions (dS) in PBR codons. Highly diverged sequences showed no evidence for diversifying selection, possibly because synonymous substitutions have accumulated since their divergence. Also, the substitution pattern of similar alleles did not show evidence for diversifying selection, plausibly because advantageous non-synonymous substitutions have not yet accumulated. Intermediately diverged sequences showed the highest relative rate of non-synonymous substitutions, with dN/dS>14 in some pairwise comparisons. Consequently, a curvilinear relationship was observed between the dN/dS ratio and the level of sequence divergence.

Gyrodactylus pictae n. sp. (Monogenea: Gyrodactylidae) from the Trinidadian swamp guppy Poecilia picta Regan, with a discussion on species of Gyrodactylus von Nordmann, 1832 and their poeciliid hosts.

Gyrodactylus pictae n. sp. is recorded from Poecilia picta in heterospecific shoals with the guppy P. reticulata in Northern Trinidad. G. pictae is morphologically similar to G. turnbulli Harris, 1986, but the hamuli and marginal hooks are slightly smaller and more gracile. The toe and the point of the marginal hook have a distinctly different shape, providing the best morphological characters for distinguishing the two species. The rDNA ITS1 and ITS2 sequences differ from those of G. turnbulli (the closest relative) by >5, suggesting that these two taxa are not sibling species. The origin of the two species on poeciliids of the subgenus Micropoecilia is discussed, and it is suggested that this may represent a case of host-parasite co-evolution.

On the neutrality of molecular genetic markers: pedigree analysis of genetic variation in fragmented populations.

Many studies employ molecular markers to infer ecological and evolutionary processes, assuming that variation found at genetic loci offers a reliable representation of stochastic events in natural populations. Increasingly, evidence emerges that molecular markers might not always be selectively neutral. However, only a few studies have analysed how deviations from neutrality could affect estimates of genetic variation, using populations with known genealogy. We monitored changes in allozyme variation over eight generations in captive metapopulations of the butterfly Bicyclus anynana. Population demography was recorded by individually marking 35 000 butterflies and constructing pedigrees. We designed a computer program that simulated the inheritance of founder allozyme alleles in butterfly pedigrees. We thus tested whether the observed transmission of allozyme alleles could be explained by random genetic drift alone, or whether there was evidence for positive or negative selection. This analysis showed that in the smallest metapopulations the loss of allozyme variation exceeded the neutral rate. Possibly, linkage disequilibria between deleterious mutations and marker alleles resulted in background selection and a faster erosion of allozyme variation. In larger metapopulations, one locus (MDH) showed a significant heterozygote excess and smaller than expected loss in heterozygosity, observations consistent with (associative) overdominance. This study demonstrates that the neutrality of molecular markers cannot always be assumed, particularly in small populations with a high mutation load.

Inbreeding depression and genetic load of sexually selected traits: how the guppy lost its spots.

To date, few studies have investigated the effects of inbreeding on sexually selected traits, although inbreeding depression on such traits can play an important role in the evolution and ecology of wild populations. Sexually selected traits such as ornamentation and courtship behaviour may not be primary fitness characters, but selection and dominance coefficients of their mutations will resemble those of traits under natural selection. Strong directional selection, for instance, through female mate-choice, purges all but the most recessive deleterious mutations, and the remaining dominance variation will result in inbreeding depression once populations undergo bottlenecks. We analysed the effects of inbreeding on sexually selected traits (colour pattern and courtship behaviour) in the male guppy, Poecilia reticulata, from Trinidad, and found a significant decline in the frequency of mating behaviour and colour spots. Such effects occurred although the genetic basis of these traits, many of which are Y-linked and hemizygous, would be expected to leave relatively little scope for inbreeding depression. Findings suggest that these sexually selected traits could reflect the genetic condition or health of males, and thus may be informative mate-cue characters for female choice as suggested by the 'good genes' model.

Inbreeding depression and genetic load in laboratory metapopulations of the butterfly Bicyclus anynana.

We investigated the effects of inbreeding on various fitness components and their genetic load in laboratory metapopulations of the butterfly Bicyclus anynana. Six metapopulations each consisted of four subpopulations with breeding population sizes of N = 6 or N = 12 and migration rate of m = 0 or m = 0.33. Metapopulations were maintained for seven generations during which coancestries and pedigrees were established. Individual inbreeding coefficients at the F7 were calculated and ranged between 0.01 and 0.51. Even though considerable purging had occurred during inbreeding, the genetic load remained higher than that of many outbreeding species: approximately two lethal equivalents were detected for egg sterility, one for zygote survival, one for juvenile survival, and one for longevity. Severe inbreeding depression occurred after seven generations of inbreeding, which jeopardized the metapopulation survival. This finding suggests that the purging of genetic load by intentional inbreeding cannot be recommended for the genetic conservation of species with a high number of lethal equivalents.