Spatial and temporal variation in population genetic structure of wild Nile tilapia (Oreochromis niloticus) across Africa

Background: Reconstructing the evolutionary history of a species is challenging. It often depends not only on the past biogeographic and climatic events but also the contemporary and ecological factors, such as current connectivity and habitat heterogeneity. In fact, these factors might interact with each other and shape the current species distribution. However, to what extent the current population genetic structure reflects the past and the contemporary factors is largely unknown. Here we investigated spatio-temporal genetic structures of Nile tilapia (Oreochromis niloticus) populations, across their natural distribution in Africa. While its large biogeographic distribution can cause genetic differentiation at the paleo-biogeographic scales, its restricted dispersal capacity might induce a strong genetic structure at micro-geographic scales. Results: Using nine microsatellite loci and 350 samples from ten natural populations, we found the highest genetic differentiation among the three ichthyofaunal provinces and regions (Ethiopian, Nilotic and Sudano-Sahelian) (R(ST) = 0.38 - 0.69). This result suggests the predominant effect of paleo-geographic events at macro-geographic scale. In addition, intermediate divergences were found between rivers and lakes within the regions, presumably reflecting relatively recent interruptions of gene flow between hydrographic basins (R(ST) = 0.24 - 0.32). The lowest differentiations were observed among connected populations within a basin (R(ST) = 0.015 in the Volta basin). Comparison of temporal sample series revealed subtle changes in the gene pools in a few generations (F = 0 - 0.053). The estimated effective population sizes were 23 - 143 and the estimated migration rate was moderate (m similar to 0.094 - 0.097) in the Volta populations. Conclusions: This study revealed clear hierarchical patterns of the population genetic structuring of O. niloticus in Africa. The effects of paleo-geographic and climatic events were predominant at macro-geographic scale, and the significant effect of geographic connectivity was detected at micro-geographic scale. The estimated effective population size, the moderate level of dispersal and the rapid temporal change in genetic composition might reflect a potential effect of life history strategy on population dynamics. This hypothesis deserves further investigation. The dynamic pattern revealed at micro-geographic and temporal scales appears important from a genetic resource management as well as from a biodiversity conservation point of view.

Population genomic tests of models of adaptive radiation in Lake Victoria Region cichlid fish

Adaptive radiation is usually thought to be associated with speciation, but the evolution of intraspecific polymorphisms without speciation is also possible. The radiation of cichlid fish in Lake Victoria (LV) is perhaps the most impressive example of a recent rapid adaptive radiation, with 600+ very young species. Key questions about its origin remain poorly characterized, such as the importance of speciation versus polymorphism, whether species persist on evolutionary time scales, and if speciation happens more commonly in small isolated or in large connected populations. We used 320 individuals from 105 putative species from Lakes Victoria, Edward, Kivu, Albert, Nabugabo and Saka, in a radiation-wide amplified fragment length polymorphism (AFLP) genome scan to address some of these questions. We demonstrate pervasive signatures of speciation supporting the classical model of adaptive radiation associated with speciation. A positive relationship between the age of lakes and the average genomic differentiation of their species, and a significant fraction of molecular variance explained by above-species level taxonomy suggest the persistence of species on evolutionary time scales, with radiation through sequential speciation rather than a single starburst. Finally the large gene diversity retained from colonization to individual species in every radiation suggests large effective population sizes and makes speciation in small geographical isolates unlikely.

Do males matter? The role of males in population dynamics

Models of population dynamics generally neglect the presence of males. While this assumption holds under many circumstances, behavioural ecology increasingly tells us that the presence (or absence) of males may have an impact on female fitness, and hence population sizes. Here we ask the question of whether males matter to population dynamics, operationally defined as a dependency of population growth on the relative density of males. We provide simple models, and evaluate the current empirical evidence for them, that illustrate various mechanisms of such male influence: mate searching behavior, male resource use (including effects of sexual dimorphism), sexual harassment and sexual segregation. In each case, theory predicts that males can have an effect on population densities, and in some extreme cases a positive feedback between an increasingly male-biased sex ratio and the effects on female harassment may theoretically even bring about population extinction. The results of this study, and the literature reviewed, show that the males can have a substantial effect on population dynamics, particularly so when human influences result in biased sex ratios.

ABCtoolbox: a versatile toolkit for approximate Bayesian computations

Background The estimation of demographic parameters from genetic data often requires the computation of likelihoods. However, the likelihood function is computationally intractable for many realistic evolutionary models, and the use of Bayesian inference has therefore been limited to very simple models. The situation changed recently with the advent of Approximate Bayesian Computation (ABC) algorithms allowing one to obtain parameter posterior distributions based on simulations not requiring likelihood computations. Results Here we present ABCtoolbox, a series of open source programs to perform Approximate Bayesian Computations (ABC). It implements various ABC algorithms including rejection sampling, MCMC without likelihood, a Particle-based sampler and ABC-GLM. ABCtoolbox is bundled with, but not limited to, a program that allows parameter inference in a population genetics context and the simultaneous use of different types of markers with different ploidy levels. In addition, ABCtoolbox can also interact with most simulation and summary statistics computation programs. The usability of the ABCtoolbox is demonstrated by inferring the evolutionary history of two evolutionary lineages of Microtus arvalis. Using nuclear microsatellites and mitochondrial sequence data in the same estimation procedure enabled us to infer sex-specific population sizes and migration rates and to find that males show smaller population sizes but much higher levels of migration than females. Conclusion ABCtoolbox allows a user to perform all the necessary steps of a full ABC analysis, from parameter sampling from prior distributions, data simulations, computation of summary statistics, estimation of posterior distributions, model choice, validation of the estimation procedure, and visualization of the results.

Monitoring Butterfly Abundance: Beyond Pollard Walks

Most butterfly monitoring protocols rely on counts along transects (Pollard walks) to generate species abundance indices and track population trends. It is still too often ignored that a population count results from two processes: the biological process (true abundance) and the statistical process (our ability to properly quantify abundance). Because individual detectability tends to vary in space (e.g., among sites) and time (e.g., among years), it remains unclear whether index counts truly reflect population sizes and trends. This study compares capture-mark-recapture (absolute abundance) and count-index (relative abundance) monitoring methods in three species (Maculinea nausithous and Iolana iolas: Lycaenidae; Minois dryas: Satyridae) in contrasted habitat types. We demonstrate that intraspecific variability in individual detectability under standard monitoring conditions is probably the rule rather than the exception, which questions the reliability of count-based indices to estimate and compare specific population abundance. Our results suggest that the accuracy of count-based methods depends heavily on the ecology and behavior of the target species, as well as on the type of habitat in which surveys take place. Monitoring programs designed to assess the abundance and trends in butterfly populations should incorporate a measure of detectability. We discuss the relative advantages and inconveniences of current monitoring methods and analytical approaches with respect to the characteristics of the species under scrutiny and resources availability.

The Oxytricha trifallax macronuclear genome: a complex eukaryotic genome with 16,000 tiny chromosomes.

The macronuclear genome of the ciliate Oxytricha trifallax displays an extreme and unique eukaryotic genome architecture with extensive genomic variation. During sexual genome development, the expressed, somatic macronuclear genome is whittled down to the genic portion of a small fraction (∼5%) of its precursor "silent" germline micronuclear genome by a process of "unscrambling" and fragmentation. The tiny macronuclear "nanochromosomes" typically encode single, protein-coding genes (a small portion, 10%, encode 2-8 genes), have minimal noncoding regions, and are differentially amplified to an average of ∼2,000 copies. We report the high-quality genome assembly of ∼16,000 complete nanochromosomes (∼50 Mb haploid genome size) that vary from 469 bp to 66 kb long (mean ∼3.2 kb) and encode ∼18,500 genes. Alternative DNA fragmentation processes ∼10% of the nanochromosomes into multiple isoforms that usually encode complete genes. Nucleotide diversity in the macronucleus is very high (SNP heterozygosity is ∼4.0%), suggesting that Oxytricha trifallax may have one of the largest known effective population sizes of eukaryotes. Comparison to other ciliates with nonscrambled genomes and long macronuclear chromosomes (on the order of 100 kb) suggests several candidate proteins that could be involved in genome rearrangement, including domesticated MULE and IS1595-like DDE transposases. The assembly of the highly fragmented Oxytricha macronuclear genome is the first completed genome with such an unusual architecture. This genome sequence provides tantalizing glimpses into novel molecular biology and evolution. For example, Oxytricha maintains tens of millions of telomeres per cell and has also evolved an intriguing expansion of telomere end-binding proteins. In conjunction with the micronuclear genome in progress, the O. trifallax macronuclear genome will provide an invaluable resource for investigating programmed genome rearrangements, complementing studies of rearrangements arising during evolution and disease.

Minimizing the Threat of a Positive Majority Deficit in Two-tier Voting Systems with Equipopulous Units

The mean majority deficit in a two-tier voting system is a function of the partition of the population. We derive a new square-root rule: For odd-numbered population sizes and equipopulous units the mean majority deficit is maximal when the member size of the units in the partition is close to the square root of the population size. Furthermore, within the partitions into roughly equipopulous units, partitions with small even numbers of units or small even-sized units yield high mean majority deficits. We discuss the implications for the winner-takes-all system in the US Electoral College.

Living in Heterogeneous Woodlands – Are Habitat Continuity or Quality Drivers of Genetic Variability in a Flightless Ground Beetle?

Although genetic diversity is one of the key components of biodiversity, its drivers are still not fully understood. While it is known that genetic diversity is affected both by environmental parameters as well as habitat history, these factors are not often tested together. Therefore, we analyzed 14 microsatellite loci in Abax parallelepipedus, a flightless, forest dwelling ground beetle, from 88 plots in two study regions in Germany. We modeled the effects of historical and environmental variables on allelic richness, and found for one of the regions, the Schorfheide-Chorin, a significant effect of the depth of the litter layer, which is a main component of habitat quality, and of the sampling effort, which serves as an inverse proxy for local population size. For the other region, the Schwäbische Alb, none of the potential drivers showed a significant effect on allelic richness. We conclude that the genetic diversity in our study species is being driven by current local population sizes via environmental variables and not by historical processes in the studied regions. This is also supported by lack of genetic differentiation between local populations sampled from ancient and from recent woodlands. We suggest that the potential effects of former fragmentation and recolonization processes have been mitigated by the large and stable local populations of Abax parallelepipedus in combination with the proximity of the ancient and recent woodlands in the studied landscapes.

Spatiotemporal dynamics of Puumala hantavirus associated with its rodent host, Myodes glareolus

Many viruses significantly impact human and animal health. Understanding the population dynamics of these viruses and their hosts can provide important insights for epidemiology and virus evolution. Puumala virus (PUUV) is a European hantavirus that may cause regional outbreaks of hemorrhagic fever with renal syndrome in humans. Here, we analyzed the spatiotemporal dynamics of PUUV circulating in local populations of its rodent reservoir host, the bank vole (Myodes glareolus) during eight years. Phylogenetic and population genetic analyses of all three genome segments of PUUV showed strong geographical structuring at a very local scale. There was a high temporal turnover of virus strains in the local bank vole populations, but several virus strains persisted through multiple years. Phylodynamic analyses showed no significant changes in the local effective population sizes of PUUV, although vole numbers and virus prevalence fluctuated widely. Microsatellite data demonstrated also a temporally persisting subdivision between local vole populations, but these groups did not correspond to the subdivision in the virus strains. We conclude that restricted transmission between vole populations and genetic drift play important roles in shaping the genetic structure and temporal dynamics of PUUV in its natural host which has several implications for zoonotic risks of the human population.

[Occurrence of Clostridium perfringens type A and type C in piglets of the Swiss swine population]

Necrotizing enteritis (NE) of newborn piglets still represents an economical problem in Swiss pig breeding and production. The aim of our study was to identify risk factors for NE and evaluate the prevalence of C. perfringens with the toxingenes cpb and cpb2 in Swiss pig breeding farms. The prevalence of theses C. perfringens was investigated using fecal swabs followed by bacteriological culturing and genotyping. Close proximity to other breeding farms and large herd sizes were shown to predispose to NE. C. perfringens type C, carrying the genes cpa, cpb and cpb2 were frequently identified in herds with acute outbreaks of NE. Farms not affected by NE or those using prophylactic vaccination against NE were predominantly positive for C. perfringens type A strains with cpb2 and showed much lower prevalence of C. perfringens type C, compared to acutely affected herds. Our results demonstrate that C. perfringens type A strains with cpb2 are not associated with NE. Besides typical necropsy finding, only the identification of cpb can be used for the diagnosis of NE in affected herds.

Comparison of pain-resilient working individuals to population-based case controls with/without momentary low back pain

Background: Few studies have examined the 20% of individuals who never experience an episode of low back pain (LBP). To date, no investigation has been undertaken that examines a group who claim to have never experienced LBP in their lifetime in comparison to two population-based case–control groups with and without momentary LBP. This study investigates whether LBP-resilient workers between 50 and 65 years had better general health, demonstrated more positive health behaviour and were better able to achieve routine activities compared with both case–control groups. Methods: Forty-two LBP-resilient participants completed the same pain assessment questionnaire as a population-based LBP sample from a nationwide, large-scale cross-sectional survey in Switzerland. The LBP-resilient participants were pairwise compared to the propensity score-matched case controls by exploring differences in demographic and work characteristics, and by calculating odds ratios (ORs) and effect sizes. A discriminant analysis explored group differences, while the multiple logistic regression analysis specified single indicators which accounted for group differences. Results: LBP-resilient participants were healthier than the case controls with momentary LBP and achieved routine activities more easily. Compared to controls without momentary LBP, LBP-resilient participants had a higher vitality, a lower workload, a healthier attitude towards health and behaved more healthily by drinking less alcohol. Conclusions: By demonstrating a difference between LBP-resilient participants and controls without momentary LBP, the question that arises is what additional knowledge can be attained. Three underlying traits seem to be relevant about LBP-resilient participants: personality, favourable work conditions and subjective attitudes/attributions towards health. These rationales have to be considered with respect to LBP prevention.