Adaptive approximate Bayesian computation

Sequential techniques can enhance the efficiency of the approximate Bayesian computation algorithm, as in Sisson et al.'s (2007) partial rejection control version. While this method is based upon the theoretical works of Del Moral et al. (2006), the application to approximate Bayesian computation results in a bias in the approximation to the posterior. An alternative version based on genuine importance sampling arguments bypasses this difficulty, in connection with the population Monte Carlo method of Cappe et al. (2004), and it includes an automatic scaling of the forward kernel. When applied to a population genetics example, it compares favourably with two other versions of the approximate algorithm.

Recent developments in genetic data analysis: what can they tell us about human demographic history?

Over the last decade, a number of new methods of population genetic analysis based on likelihood have been introduced. This review describes and explains the general statistical techniques that have recently been used, and discusses the underlying population genetic models. Experimental papers that use these methods to infer human demographic and phylogeographic history are reviewed. It appears that the use of likelihood has hitherto had little impact in the field of human population genetics, which is still primarily driven by more traditional approaches. However, with the current uncertainty about the effects of natural selection, population structure and ascertainment of single-nucleotide polymorphism markers, it is suggested that likelihood-based methods may have a greater impact in the future.

Identifying adaptive genetic divergence among populations from genome scans

The identification of signatures of natural selection in genomic surveys has become an area of intense research, stimulated by the increasing ease with which genetic markers can be typed. Loci identified as subject to selection may be functionally important, and hence (weak) candidates for involvement in disease causation. They can also be useful in determining the adaptive differentiation of populations, and exploring hypotheses about speciation. Adaptive differentiation has traditionally been identified from differences in allele frequencies among different populations, summarised by an estimate of F-ST. Low outliers relative to an appropriate neutral population-genetics model indicate loci subject to balancing selection, whereas high outliers suggest adaptive (directional) selection. However, the problem of identifying statistically significant departures from neutrality is complicated by confounding effects on the distribution of F-ST estimates, and current methods have not yet been tested in large-scale simulation experiments. Here, we simulate data from a structured population at many unlinked, diallelic loci that are predominantly neutral but with some loci subject to adaptive or balancing selection. We develop a hierarchical-Bayesian method, implemented via Markov chain Monte Carlo (MCMC), and assess its performance in distinguishing the loci simulated under selection from the neutral loci. We also compare this performance with that of a frequentist method, based on moment-based estimates of F-ST. We find that both methods can identify loci subject to adaptive selection when the selection coefficient is at least five times the migration rate. Neither method could reliably distinguish loci under balancing selection in our simulations, even when the selection coefficient is twenty times the migration rate.

Evidence for a residual postglacial founder effect in a highly dispersive freshwater invertebrate

In the northern hemisphere, glacial periods have had profound and lasting effects on the population genetics of numerous species, with founder effects often persisting for many generations in formerly glaciated regions. We found an unusual example of this in the freshwater bryozoan Cristatella mucedo, which showed regional differences in haplotype diversity with relatively low levels of haplotype diversity in northern Europe compared to central/southern Europe despite previous evidence for frequent dispersal between C. mucedo populations. Such contradictions between high dispersal and low gene flow have now been reported in several other freshwater taxa and may be attributed to persistent founder effects following colonization of sites by a few individuals whose efficient reproduction leads to rapid population growth. Alternatively, selection may determine which genotypes can thrive in northerly locations, or it may be that C. mucedo has undergone cryptic speciation. Future work on adaptive genomic regions is required before we can understand how gene flow, local adaptation, and speciation influence the current distribution patterns of bryozoans and other freshwater invertebrates.

Pseudogene rescue: an adaptive mechanism of codon reassignment

Alterations to the genetic code – codon reassignments – have occurred many times in life’s history, despite the fact that genomes are coadapted to their genetic codes and therefore alterations are likely to be maladaptive. A potential mechanism for adaptive codon reassignment, which could trigger either a temporary period of codon ambiguity or a permanent genetic code change, is the reactivation of a pseudogene by a nonsense suppressor mutant transfer RNA. I examine the population genetics of each stage of this process and find that pseudogene rescue is plausible and also readily explains some features of extant variability in genetic codes.

A two-locus forensic match probability for subdivided populations

A two-locus match probability is presented that incorporates the effects of within-subpopulation inbreeding (consanguinity) in addition to population subdivision. The usual practice of calculating multi-locus match probabilities as the product of single-locus probabilities assumes independence between loci. There are a number of population genetics phenomena that can violate this assumption: in addition to consanguinity, which increases homozygosity at all loci simultaneously, gametic disequilibrium will introduce dependence into DNA profiles. However, in forensics the latter problem is usually addressed in part by the careful choice of unlinked loci. Hence, as is conventional, we assume gametic equilibrium here, and focus instead on between-locus dependence due to consanguinity. The resulting match probability formulae are an extension of existing methods in the literature, and are shown to be more conservative than these methods in the case of double homozygote matches. For two-locus profiles involving one or more heterozygous genotypes, results are similar to, or smaller than, the existing approaches.

Selfishness and altruism can coexist when help is subject to diminishing returns

Altruism and selfishness are 30–50% heritable in man in both Western and non-Western populations. This genetically based variation in altruism and selfishness requires explanation. In non-human animals, altruism is generally directed towards relatives, and satisfies the condition known as Hamilton's rule. This nepotistic altruism evolves under natural selection only if the ratio of the benefit of receiving help to the cost of giving it exceeds a value that depends on the relatedness of the individuals involved. Standard analyses assume that the benefit provided by each individual is the same but it is plausible in some cases that as more individuals contribute, help is subject to diminishing returns. We analyse this situation using a single-locus two-allele model of selection in a diploid population with the altruistic allele dominant to the selfish allele. The analysis requires calculation of the relationship between the fitnesses of the genotypes and the frequencies of the genes. The fitnesses vary not only with the genotype of the individual but also with the distribution of phenotypes amongst the sibs of the individual and this depends on the genotypes of his parents. These calculations are not possible by direct fitness or ESS methods but are possible using population genetics. Our analysis shows that diminishing returns change the operation of natural selection and the outcome can now be a stable equilibrium between altruistic and selfish alleles rather than the elimination of one allele or the other. We thus provide a plausible genetic model of kin selection that leads to the stable coexistence in the same population of both altruistic and selfish individuals. This may explain reported genetic variation in altruism in man.

Likelihood-free estimation of model evidence

Statistical methods of inference typically require the likelihood function to be computable in a reasonable amount of time. The class of “likelihood-free” methods termed Approximate Bayesian Computation (ABC) is able to eliminate this requirement, replacing the evaluation of the likelihood with simulation from it. Likelihood-free methods have gained in efficiency and popularity in the past few years, following their integration with Markov Chain Monte Carlo (MCMC) and Sequential Monte Carlo (SMC) in order to better explore the parameter space. They have been applied primarily to estimating the parameters of a given model, but can also be used to compare models. Here we present novel likelihood-free approaches to model comparison, based upon the independent estimation of the evidence of each model under study. Key advantages of these approaches over previous techniques are that they allow the exploitation of MCMC or SMC algorithms for exploring the parameter space, and that they do not require a sampler able to mix between models. We validate the proposed methods using a simple exponential family problem before providing a realistic problem from human population genetics: the comparison of different demographic models based upon genetic data from the Y chromosome.

Evolution of discrimination in populations at equilibrium between selfishness and altruism

Where there is genetically based variation in selfishness and altruism, as in man, altruists with an innate ability to recognise and thereby only help their altruistic relatives may evolve. Here we use diploid population genetic models to chart the evolution of genetically-based discrimination in populations initially in stable equilibrium between altruism and selfishness. The initial stable equilibria occur because help is assumed subject to diminishing returns. Similar results were obtained whether we used a model with two independently inherited loci, one controlling altruism the other discrimination, or a one locus model with three alleles. The latter is the opposite extreme to the first model, and can be thought of as involving complete linkage between two loci on the same chromosome. The introduction of discrimination reduced the benefits obtained by selfish individuals, more so as the number of discriminators increased, and selfishness was eventually eliminated in some cases. In others selfishness persisted and the evolutionary outcome was a stable equilibrium involving selfish individuals and both discriminating and non-discriminating altruists. Heritable variation in selfishness, altruism and discrimination is predicted to be particularly evident among full sibs. The suggested coexistence of these three genetic dispositions could explain widespread interest within human social groups as to who will and who will not help others. These predictions merit experimental and observational investigation by primatologists, anthropologists and psychologists. Keywords: Population genetics, Diploid, Heritability, Prosocial, Behaviour genetics

Late Pleistocene climate change and the global expansion of anatomically modern humans

The extent to which past climate change has dictated the pattern and timing of the out-of-Africa expansion by anatomically modern humans is currently unclear [Stewart JR, Stringer CB (2012) Science 335:1317–1321]. In particular, the incompleteness of the fossil record makes it difficult to quantify the effect of climate. Here, we take a different approach to this problem; rather than relying on the appearance of fossils or archaeological evidence to determine arrival times in different parts of the world, we use patterns of genetic variation in modern human populations to determine the plausibility of past demographic parameters. We develop a spatially explicit model of the expansion of anatomically modern humans and use climate reconstructions over the past 120 ky based on the Hadley Centre global climate model HadCM3 to quantify the possible effects of climate on human demography. The combinations of demographic parameters compatible with the current genetic makeup of worldwide populations indicate a clear effect of climate on past population densities. Our estimates of this effect, based on population genetics, capture the observed relationship between current climate and population density in modern hunter–gatherers worldwide, providing supporting evidence for the realism of our approach. Furthermore, although we did not use any archaeological and anthropological data to inform the model, the arrival times in different continents predicted by our model are also broadly consistent with the fossil and archaeological records. Our framework provides the most accurate spatiotemporal reconstruction of human demographic history available at present and will allow for a greater integration of genetic and archaeological evidence.

Isolation of novel multipotent neural crest-derived stem cells from adult human inferior turbinate

Adult human neural crest-derived stem cells (NCSCs) are of extraordinary high plasticity and promising candidates for the use in regenerative medicine. Here we describe for the first time a novel neural crest-derived stem cell population within the respiratory epithelium of human adult inferior turbinate. In contrast to superior and middle turbinates, high amounts of source material could be isolated from human inferior turbinates. Using minimally-invasive surgery methods isolation is efficient even in older patients. Within their endogenous niche, inferior turbinate stem cells (ITSCs) expressed high levels of nestin, p75(NTR), and S100. Immunoelectron microscopy using anti-p75 antibodies displayed that ITSCs are of glial origin and closely related to nonmyelinating Schwann cells. Cultivated ITSCs were positive for nestin and S100 and the neural crest markers Slug and SOX10. Whole genome microarray analysis showed pronounced differences to human ES cells in respect to pluripotency markers OCT4, SOX2, LIN28, and NANOG, whereas expression of WDR5, KLF4, and c-MYC was nearly similar. ITSCs were able to differentiate into cells with neuro-ectodermal and mesodermal phenotype. Additionally ITSCs are able to survive and perform neural crest typical chain migration in vivo when transplanted into chicken embryos. However ITSCs do not form teratomas in severe combined immunodeficient mice. Finally, we developed a separation strategy based on magnetic cell sorting of p75(NTR) positive ITSCs that formed larger neurospheres and proliferated faster than p75(NTR) negative ITSCs. Taken together our study describes a novel, readily accessible source of multipotent human NCSCs for potential cell-replacement therapy.

Tales of the unexpected: research discoveries that occurred by accident

Data from four experimental research projects are presented which have in common that unexpected results caused a change in direction of the research. A plant growth accelerator caused the appearance of white black bean aphids, a synthetic pyrethroid suspected of enhancing aphid reproduction proved to enhance plant growth, a chance conversation with a colleague initiated a search for fungal DNA in aphids, and the accidental invasion of aphid cultures by a parasitoid reversed the aphid population ranking of two Brussels sprout cultivars. This last result led to a whole series of studies on the plant odour preferences of emerging parasitoids which in turn revealed the unexpected phenomenon that chemical cues to the maternal host plant are left with the eggs at oviposition. It is pointed out that, too often, researchers fail to follow up unexpected results because they resist accepting flaws in their hypotheses; also that current application criteria for research funding make it hard to accommodate unexpected findings.

Diruthenium complexes with bridging diethynyl polyaromatic ligands: synthesis, spectroelectrochemistry, and theoretical calculations

This work describes syntheses and electrochemical, spectroscopic, and bonding properties in a new series of dinuclear ruthenium(II) complexes bridged by polyaromatic (biphenyl, fluorene, phenanthrene, and pyrene) alkynyl ligands. Longitudinal expansion of the π-conjugated polyaromatic core of the bridging ligands caused a reduced potential difference between the anodic steps and reinforced their bridge-localized nature, as evidenced by UV/vis/near-IR and IR spectroelectrochemical data combined with DFT and TDDFT calculations. Importantly, the intricate multiple IR ν(CC) absorption bands for the singly oxidized states imply a thermal population of a range of conformers (rotamers) with distinct electronic character. This behavior was demonstrated with more accurate DFT calculations of selected nontruncated 1e− oxidized complexes in three different conformations. The combined experimental and theoretical data reveal that thermally populated rotamers featuring various mutual orientations of the ligated metal termini and the bridging diethynyl polyaromatic moieties have a significant impact on the electronic absorption and ν(CC) wavenumbers of the singly oxidized systems.

Concordant Phylogeographies of 2 Malaria Vectors Attest to Common Spatial and Demographic Histories

The phylogeography of South American lineages is a topic of heated debate. Although a single process is unlikely to describe entire ecosystems, related species, which incur similar habitat limitations, can inform the history for a subsection of assemblages. We compared the phylogeographic patterns of the cytochrome oxidase I marker from Anopheles triannulatus (N = 72) and previous results for A. darlingi (N = 126) in a broad portion of their South American distributions. Both species share similar population subdivisions, with aggregations northeast of the Amazon River, in southern coastal Brazil and 2 regions in central Brazil. The average (ST) between these groups was 0.39 for A. triannulatus. Populations northeast of the Amazon and in southeastern Brazil are generally reciprocally monophyletic to the remaining groups. Based on these initial analyses, we constructed the a priori hypothesis that the Amazon and regions of high declivity pose geographic barriers to dispersal in these taxa. Mantel tests confirmed that these areas block gene flow for more than 1000 km for both species. The efficacy of these impediments was tested using landscape genetics, which could not reject our a priori hypothesis but did reject simpler scenarios. Results form summary statistics and phylogenetics suggest that both lineages originated in central Amazonia (south of the Amazon River) during the late Pleistocene (579 000 years ago) and that they followed the same paths of expansion into their contemporary distributions. These results may have implications for other species sharing similar ecological limitations but probably are not applicable as a general paradigm of Neotropical biogeography.

Forest-obligate Sabethes mosquitoes suggest palaeoecological perturbations

The origin of tropical forest diversity has been hotly debated for decades. Although specific mechanisms vary, many such explanations propose some vicariance in the distribution of species during glacial cycles and several have been supported by genetic evidence in Neotropical taxa. However, no consensus exists with regard to the extent or time frame of the vicariance events. Here, we analyse the cytochrome oxidase II mitochondrial gene of 250 Sabethes albiprivus B mosquitoes sampled from western Sao Paulo in Brazil. There was very low population structuring among collection sites (Phi(ST) = 0.03, P = 0.04). Historic demographic analyses and the contemporary geographic distribution of genetic diversity suggest that the populations sampled are not at demographic equilibrium. Three distinct mitochondrial clades were observed in the samples, one of which differed significantly in its geographic distribution relative to the other two within a small sampling area (similar to 70 x 35 km). This fact, supported by the inability of maximum likelihood analyses to achieve adequate fits to simple models for the population demography of the species, suggests a more complex history, possibly involving disjunct forest refugia. This hypothesis is supported by a genetic signal of recent population growth, which is expected if population sizes of this forest-obligate insect increased during the forest expansions that followed glacial periods. Although a time frame cannot be reliably inferred for the vicariance event leading to the three genetic clades, molecular clock estimates place this at similar to 1 Myr before present.