Genomic Evidence for the Complex Evolutionary History of Macaques (Genus Macaca)

The genus Macaca is widely distributed, occupies a variety of habitats, shows diverse phenotypic characteristics, and is one of the best-studied genera of nonhuman primates. Here, we reported five re-sequencing Macaca genomes, including one M. cyclopis, one M. fuscata, one M. thibetana, one M. silenus, and one M. sylvanus. Together with published genomes of other macaque species, we combined 20 genome sequences of 10 macaque species to investigate the gene introgression and genetic differences among the species. The network analysis of the SNV-fragment trees indicates a reticular phylogeny of macaque species. Combining the results from various analytical methods, we identified extensive ancient introgression events among macaque species. The multiple introgression signals between different species groups were also observed, such as between fascicularis group species and silenus group species. However, gene flow signals between fascicularis and sinica group were not as strong as those between fascicularis group and silenus group. On the other hand, the unidirect gene flow in M. arctoides probably occurred between the progenitor of M. arctoides and the common ancestor of fascicularis group. Our study also shows that the genetic backgrounds and genetic diversity of different macaques vary dramatically among species, even among populations of the same species. In conclusion, using whole genome sequences and multiple methods, we have studied the evolutionary history of the genus Macaca and provided evidence for extensive introgression among the species.

Allele age estimators designed for whole genome datasets show only a modest decrease in accuracy when applied to whole exome datasets

Personalized genomics in the healthcare system is becoming increasingly accessible as the costs of sequencing decreases. With the increase in number of genomes, larger numbers of rare variants are being discovered and much work is being done to identify their functional impacts in relation to disease phenotypes. One way to characterize these variants is to estimate the time the mutation entered the population. However, allele age estimators such as Relate, Genealogical Estimator of Variant Age, and time of coalescence, were developed based on the assumption that datasets include the entire genome. We examined the performance of each of these estimators on simulated exome data under a neutral constant population size model and found that each provides usable estimates of allele age from whole-exome datasets. To test the robustness of these methods, analyses were undertaken to simulate data under a population expansion model and background selection. Relate performs the best amongst all three estimators with Pearson coefficients of 0.64 and 0.68 (neutral constant and expansion population model) with a 17 percent and 15 percent drop in accuracy between whole genome and whole exome estimations. Of the three estimators, Relate is best able to parallelize to yield quick results with little resources, however even Relate is only able to scale to thousands of samples making it unable to match the hundreds of thousands of samples being currently released. While more work is needed to expand the capabilities of current methods of estimating allele age, these methods estimate the age of mutations with a modest decrease in performance.

Analyses of allele age and fitness impact reveal human beneficial alleles to be older than neutral controls

A classic population genetic prediction is that alleles experiencing directional selection should swiftly traverse allele frequency space, leaving detectable reductions in genetic variation in linked regions. However, despite this expectation, identifying clear footprints of beneficial allele passage has proven to be surprisingly challenging. We addressed the basic premise underlying this expectation by estimating the ages of large numbers of beneficial and deleterious alleles in a human population genomic data set. Deleterious alleles were found to be young, on average, given their allele frequency. However, beneficial alleles were older on average than non-coding, non-regulatory alleles of the same frequency. This finding is not consistent with directional selection and instead indicates some type of balancing selection. Among derived beneficial alleles, those fixed in the population show higher local recombination rates than those still segregating, consistent with a model in which new beneficial alleles experience an initial period of balancing selection due to linkage disequilibrium with deleterious recessive alleles. Alleles that ultimately fix following a period of balancing selection will leave a modest 'soft' sweep impact on the local variation, consistent with the overall paucity of species-wide 'hard' sweeps in human genomes.

Divergence time estimation using ddRAD data and an isolation-with-migration model applied to water vole populations of Arvicola

Molecular dating methods of population splits are crucial in evolutionary biology, but they present important difficulties due to the complexity of the genealogical relationships of genes and past migrations between populations. Using the double digest restriction-site associated DNA (ddRAD) technique and an isolation-with-migration (IM) model, we studied the evolutionary history of water vole populations of the genus Arvicola, a group of complex evolution with fossorial and semi-aquatic ecotypes. To do this, we first estimated mutation rates of ddRAD loci using a phylogenetic approach. An IM model was then used to estimate split times and other relevant demographic parameters. A set of 300 ddRAD loci that included 85 calibrated loci resulted in good mixing and model convergence. The results showed that the two populations of A. scherman present in the Iberian Peninsula split 34 thousand years ago, during the last glaciation. In addition, the much greater divergence from its sister species, A. amphibius, may help to clarify the controversial taxonomy of the genus. We conclude that this approach, based on ddRAD data and an IM model, is highly useful for analyzing the origin of populations and species.

10-Year Clinical Comparative Study of Ceramic and Composite Veneered Metal Crowns

The rejection of composite veneerings in fixed partial dentures is primarily caused by the inadequacy of the bonds between the metal framework and the composite veneering. The development of improved veneering composites necessitates an investigation of their clinical suitability compared with ceramic materials. Nineteen patients with at least two suitable, adjacent natural teeth for crowning were treated with 64 high noble alloy crowns. The adjacent crowns were veneered with ceramic (IPS inline) and composite materials (SR Adoro). Seven follow-up examinations were carried out over a period of 10 years. The crowns were investigated for mechanical defects, periodontal parameters, and discolorations. The survival rates of the ceramic veneered metal crowns (CeMCs) and composite veneered metal crowns (CoMCs) at the 10-year follow-up were 87.1 and 87.9%, respectively. The success rates of the crowns after 10 years were 83.9% for CeMCs and 51.5% for CoMCs (log-rank test, p = 0.009). No significant differences between the groups were found in the periodontal parameters (Kruskal-Wallis one-way analysis of variance, ANOVA p ⟩ 0.05). After 10 years, discoloration patterns of the two materials differed significantly (Mann-Whitney-U-test, p = 0.017). Thus, despite the improvements associated with CoMCs, CeMCs remain the gold standard for veneered metal crowns.

The Pop-Gen Pipeline Platform: A Software Platform for Population Genomic Analyses

The Pop-Gen Pipeline Platform (PPP) is a software platform for population genomic analyses. The PPP was designed as a collection of scripts that facilitate common population genomic workflows in a consistent and standardized Python environment. Functions were developed to encompass entire workflows, including input preparation, file format conversion, various population genomic analyses, and output generation. The platform has also been developed with reproducibility and extensibility of analyses in mind. The PPP is an open-source package that is available for download and use at https://ppp.readthedocs.io/en/latest/PPP_pages/install.html.

Further resolution of the house mouse (Mus musculus) phylogeny by integration over isolation-with-migration histories

BACKGROUND

The three main subspecies of house mice, Mus musculus castaneus, Mus musculus domesticus, and Mus musculus musculus, are estimated to have diverged ~ 350-500KYA. Resolution of the details of their evolutionary history is complicated by their relatively recent divergence, ongoing gene flow among the subspecies, and complex demographic histories. Previous studies have been limited to some extent by the number of loci surveyed and/or by the scope of the method used. Here, we apply a method (IMa3) that provides an estimate of a population phylogeny while allowing for complex histories of gene exchange.

RESULTS

Results strongly support a topology with M. m. domesticus as sister to M. m. castaneus and M. m. musculus. In addition, we find evidence of gene flow between all pairs of subspecies, but that gene flow is most restricted from M. m. musculus into M. m. domesticus. Estimates of other key parameters are dependent on assumptions regarding generation time and mutation rate in house mice. Nevertheless, our results support previous findings that the effective population size, Ne, of M. m. castaneus is larger than that of the other two subspecies, that the three subspecies began diverging ~ 130 - 420KYA, and that the time between divergence events was short.

CONCLUSIONS

Joint demographic and phylogenetic analyses of genomic data provide a clearer picture of the history of divergence in house mice.

Emerging Frontiers in the Study of Molecular Evolution

A collection of the editors of Journal of Molecular Evolution have gotten together to pose a set of key challenges and future directions for the field of molecular evolution. Topics include challenges and new directions in prebiotic chemistry and the RNA world, reconstruction of early cellular genomes and proteins, macromolecular and functional evolution, evolutionary cell biology, genome evolution, molecular evolutionary ecology, viral phylodynamics, theoretical population genomics, somatic cell molecular evolution, and directed evolution. While our effort is not meant to be exhaustive, it reflects research questions and problems in the field of molecular evolution that are exciting to our editors.

The effect of undetected recombination on genealogy sampling and inference under an isolation-with-migration model

Many methods for fitting demographic models to data sets of aligned sequences rely upon an assumption that the data have a branching coalescent history without recombination within regions or loci. To mitigate the effects of the failure of this assumption, a common approach is to filter data and sample regions that pass the four-gamete criterion for recombination, an approach that allows data to run, but that is expected to detect only a minority of recombination events. A series of empirical tests of this approach were conducted using computer simulations with and without recombination for a variety of isolation-with-migration (IM) model for two and three populations. Only the IMa3 program was used, but the general results should apply to related genealogy-sampling-based methods for IM models or subsets of IM models. It was found that the details of sampling intervals that pass a four-gamete filter have a moderate effect, and that schemes that use the longest intervals, or that use overlapping intervals, gave poorer results. A simple approach of using a random nonoverlapping interval returned the smallest difference between results with and without recombination, with the mean difference between parameter estimates usually less than 20% of the true value (usually much less). However, the posterior probability distributions for migration rates were flatter with recombination, suggesting that filtering based on the four-gamete criterion, while necessary for methods like these, leads to reduced resolution on migration. A distinct, alternative approach, of using a finite sites mutation model and not filtering the data, performed quite poorly.

An estimator of first coalescent time reveals selection on young variants and large heterogeneity in rare allele ages among human populations

Allele age has long been a focus of population genetic research, primarily because it can be an important clue to the fitness effects of an allele. By virtue of their effects on fitness, alleles under directional selection are expected to be younger than neutral alleles of the same frequency. We developed a new coalescent-based estimator of a close proxy for allele age, the time when a copy of an allele first shares common ancestry with other chromosomes in a sample not carrying that allele. The estimator performs well, including for the very rarest of alleles that occur just once in a sample, with a bias that is typically negative. The estimator is mostly insensitive to population demography and to factors that can arise in population genomic pipelines, including the statistical phasing of chromosomes. Applications to 1000 Genomes Data and UK10K genome data confirm predictions that singleton alleles that alter proteins are significantly younger than those that do not, with a greater difference in the larger UK10K dataset, as expected. The 1000 Genomes populations varied markedly in their distributions for singleton allele ages, suggesting that these distributions can be used to inform models of demographic history, including recent events that are only revealed by their impacts on the ages of very rare alleles.

A population genetic assessment of taxonomic species: The case of Lake Malawi cichlid fishes

Organisms sampled for population‐level research are typically assigned to species by morphological criteria. However, if those criteria are limited to one sex or life stage, or the organisms come from a complex of closely related forms, the species assignments may misdirect analyses. The impact of such sampling can be assessed from the correspondence of genetic clusters, identified only from patterns of genetic variation, to the species identified using only phenotypic criteria. We undertook this protocol with the rock‐dwelling mbuna cichlids of Lake Malawi, for which species within genera are usually identified using adult male coloration patterns. Given high local endemism of male colour patterns, and considerable allele sharing among species, there persists considerable taxonomic uncertainty in these fishes. Over 700 individuals from a single transect were photographed, genotyped and separately assigned: (a) to morphospecies using photographs; and (b) to genetic clusters using five widely used methods. Overall, the correspondence between clustering methods was strong for larger clusters, but methods varied widely in estimated number of clusters. The correspondence between morphospecies and genetic clusters was also strong for larger clusters, as well as some smaller clusters for some methods. These analyses generally affirm (a) adult male‐limited sampling and (b) the taxonomic status of Lake Malawi mbuna, as the species in our study largely appear to be well‐demarcated genetic entities. More generally, our analyses highlight the challenges for clustering methods when the number of populations is unknown, especially in cases of highly uneven sample sizes.

Phylogeny Estimation by Integration over Isolation with Migration Models

Phylogeny estimation is difficult for closely related populations and species, especially if they have been exchanging genes. We present a hierarchical Bayesian, Markov-chain Monte Carlo method with a state space that includes all possible phylogenies in a full Isolation-with-Migration model framework. The method is based on a new type of genealogy augmentation called a "hidden genealogy" that enables efficient updating of the phylogeny. This is the first likelihood-based method to fully incorporate directional gene flow and genetic drift for estimation of a species or population phylogeny. Application to human hunter-gatherer populations from Africa revealed a clear phylogenetic history, with strong support for gene exchange with an unsampled ghost population, and relatively ancient divergence between a ghost population and modern human populations, consistent with human/archaic divergence. In contrast, a study of five chimpanzee populations reveals a clear phylogeny with several pairs of populations having exchanged DNA, but does not support a history with an unsampled ghost population.

Implementing a new curriculum for computer-assisted restorations in prosthetic dentistry

INTRODUCTION

Computer-aided design/computer-aided manufacturing (CAD/CAM) of fixed prosthetic restorations has gained popularity in the last decade. However, this field of dentistry has not been integrated in the dental curriculum at most universities.

MATERIAL AND METHODS

According to the method of Kern, a curriculum was designed and established on a voluntary basis in the prosthetic education of a German dental school. The success of the implementation was measured by evaluation carried out by the participants on a visual analogue scale. Furthermore, the clinical performance of the fabricated restorations was evaluated.

RESULTS

Ninety-four percent of all students participated in the CAD/CAM curriculum indicating considerable interest. Nearly half of all students used the acquired knowledge to design crowns for their patients. All restorations fabricated by participants of the new CAD/CAM programme showed good clinical performance.

DISCUSSION

By phasing-in the CAD/CAM training programme, independent CAD/CAM-based fabrication of all-ceramic crowns increased student's self-confidence in tooth preparation. A tendency was found that students using CAD/CAM technology prepared more teeth than their fellow students who did not use CAD/CAM technology. Further studies are required to investigate the influence of independent CAD/CAM-based single-crown fabrication on the quality of the preparation.

Bayesian Analysis of Evolutionary Divergence with Genomic Data under Diverse Demographic Models

We present a new Bayesian method for estimating demographic and phylogenetic history using population genomic data. Several key innovations are introduced that allow the study of diverse models within an Isolation-with-Migration framework. The new method implements a 2-step analysis, with an initial Markov chain Monte Carlo (MCMC) phase that samples simple coalescent trees, followed by the calculation of the joint posterior density for the parameters of a demographic model. In step 1, the MCMC sampling phase, the method uses a reduced state space, consisting of coalescent trees without migration paths, and a simple importance sampling distribution without the demography of interest. Once obtained, a single sample of trees can be used in step 2 to calculate the joint posterior density for model parameters under multiple diverse demographic models, without having to repeat MCMC runs. Because migration paths are not included in the state space of the MCMC phase, but rather are handled by analytic integration in step 2 of the analysis, the method is scalable to a large number of loci with excellent MCMC mixing properties. With an implementation of the new method in the computer program MIST, we demonstrate the method's accuracy, scalability, and other advantages using simulated data and DNA sequences of two common chimpanzee subspecies: Pan troglodytes (P. t.) troglodytes and P. t. verus.

USING HITCHHIKING GENES TO STUDY ADAPTATION AND DIVERGENCE DURING SPECIATION WITHIN THE DROSOPHILA MELANOGASTER SPECIES COMPLEX

Several studies of intraspecific and interspecific DNA sequence variation from Drosophila loci have revealed a pattern of low intraspecific variation from genomic regions of low recombination. The mechanisms consistently invoked to explain these patterns are the selective sweep of advantageous mutations together with genetic hitchhiking of linked loci. To examine the effect of selective sweeps on genetic divergence during speciation, we studied two loci in different genomic regions thought to be subject to selective sweeps. We obtained DNA sequences from 1.1kb pair portions of the fourth chromosome locus cubitus interruptus Dominant (ci^D ) and from the asense locus near the telomere of the X chromosome. At ci^D , we found very low variation among multiple lines of Drosophila mauritiana and D. sechellia. This finding is consistent with an earlier report of very low variation in D. melanogaster and D. simulans at ci^D and supports the conclusion of selective sweeps and genetic hitchhiking on the nonrecombining fourth chromosome. The pattern of variation found at asense suggests that a selective sweep has occurred recently at the tip of the X chromosome in D. simulans, but not in D. melanogaster or D. mauritiana. The data from ci^D and asense are compared with data from three X chromosome loci (period, zeste, and yolk protein 2) that experience normal levels of recombination. By examining estimated genealogies and the rates at which different classes of mutations have accumulated, we conclude that selective sweeps are common occurrences on the fourth chromosome but less common near the tip of the X chromosome. An interesting pattern of low variation at ci^D among D. simulans, D. mauritiana, and D. sechellia suggests that a selective sweep may have occurred among these forms even after divergence into separate species had begun.

IMGui-A Desktop GUI Application for Isolation with Migration Analyses

The Isolation with Migration (IM) programs (e.g., IMa2) have been utilized extensively by evolutionary biologists for model-based inference of demographic parameters including effective population sizes, migration rates, and divergence times. Here, we describe a graphical user interface for the latest IM program. IMGui provides a comprehensive set of tools for performing demographic analyses, tracking progress of runs, and visualizing results. Developed using node. js and the Electron framework, IMGui is an application that runs on any desktop operating system, and is available for download at https://github.com/jaredgk/IMgui-electron-packages.

On the occurrence of false positives in tests of migration under an isolation-with-migration model

The population genetic study of divergence is often carried out using a Bayesian genealogy sampler, like those implemented in ima2 and related programs, and these analyses frequently include a likelihood ratio test of the null hypothesis of no migration between populations. Cruickshank and Hahn (2014, Molecular Ecology, 23, 3133-3157) recently reported a high rate of false-positive test results with ima2 for data simulated with small numbers of loci under models with no migration and recent splitting times. We confirm these findings and discover that they are caused by a failure of the assumptions underlying likelihood ratio tests that arises when using marginal likelihoods for a subset of model parameters. We also show that for small data sets, with little divergence between samples from two populations, an excellent fit can often be found by a model with a low migration rate and recent splitting time and a model with a high migration rate and a deep splitting time.

Inferring Very Recent Population Growth Rate from Population-Scale Sequencing Data: Using a Large-Sample Coalescent Estimator

Large-sample or population-level sequencing data provide unprecedented opportunities for inferring detailed population histories, especially recent demographic histories. On the other hand, it challenges most existing population genetic methods: Simulation-based approaches require intensive computation, and analytical approaches are often numerically intractable when the sample size is large. We propose a computationally efficient method for simultaneous estimation of population size, the rate, and onset time of population growth in the very recent history, using the pattern of the total number of segregating sites as a function of sample size. Coalescent simulation shows that it can accurately and efficiently estimate the parameters of recent population growth from large-scale data. This approach has the flexibility to model population history with multiple growth stages or other epochs, and it is robust when the sample size is very large or at the population scale, for which the Kingman's coalescent assumption is not valid. This approach is applied to recently published data and estimates the recent population growth rate in the European population to be 1.49% with the onset time 7.26 ka, and the rate in the African population to be 0.735% with the onset time 10.01 ka.

IMa2p--parallel MCMC and inference of ancient demography under the Isolation with migration (IM) model

IMa2 and related programs are used to study the divergence of closely related species and of populations within species. These methods are based on the sampling of genealogies using MCMC, and they can proceed quite slowly for larger data sets. We describe a parallel implementation, called IMa2p, that provides a nearly linear increase in genealogy sampling rate with the number of processors in use. IMa2p is written in OpenMPI and C++, and scales well for demographic analyses of a large number of loci and populations, which are difficult to study using the serial version of the program.

A hidden Markov model for investigating recent positive selection through haplotype structure

Recent positive selection can increase the frequency of an advantageous mutant rapidly enough that a relatively long ancestral haplotype will be remained intact around it. We present a hidden Markov model (HMM) to identify such haplotype structures. With HMM identified haplotype structures, a population genetic model for the extent of ancestral haplotypes is then adopted for parameter inference of the selection intensity and the allele age. Simulations show that this method can detect selection under a wide range of conditions and has higher power than the existing frequency spectrum-based method. In addition, it provides good estimate of the selection coefficients and allele ages for strong selection. The method analyzes large data sets in a reasonable amount of running time. This method is applied to HapMap III data for a genome scan, and identifies a list of candidate regions putatively under recent positive selection. It is also applied to several genes known to be under recent positive selection, including the LCT, KITLG and TYRP1 genes in Northern Europeans, and OCA2 in East Asians, to estimate their allele ages and selection coefficients.

Apparent variation in Neanderthal admixture among African populations is consistent with gene flow from Non-African populations

Recent studies have found evidence of introgression from Neanderthals into modern humans outside of sub-Saharan Africa. Given the geographic range of Neanderthals, the findings have been interpreted as evidence of gene exchange between Neanderthals and modern humans descended from the Out-of-Africa (OOA) migration. Here, we examine an alternative interpretation in which the introgression occurred earlier within Africa, between ancestors or relatives of Neanderthals and a subset of African modern humans who were the ancestors of those involved in the OOA migration. Under the alternative model, if the population structure among present-day Africans predates the OOA migration, we might find some African populations show a signal of Neanderthal introgression whereas others do not. To test this alternative model, we compiled a whole-genome data set including 38 sub-Saharan Africans from eight populations and 25 non-African individuals from five populations. We assessed differences in the amount of Neanderthal-like single-nucleotide polymorphism alleles among these populations and observed up to 1.5% difference in the number of Neanderthal-like alleles among African populations. Further analyses suggest that these differences are likely due to recent non-African admixture in these populations. After accounting for recent non-African admixture, our results do not support the alternative model of older (e.g., >100 kya) admixture between modern humans and Neanderthal-like hominids within Africa.

Identification of cichlid fishes from Lake Malawi using computer vision

BACKGROUND

The explosively radiating evolution of cichlid fishes of Lake Malawi has yielded an amazing number of haplochromine species estimated as many as 500 to 800 with a surprising degree of diversity not only in color and stripe pattern but also in the shape of jaw and body among them. As these morphological diversities have been a central subject of adaptive speciation and taxonomic classification, such high diversity could serve as a foundation for automation of species identification of cichlids.

METHODOLOGY/PRINCIPAL FINDING

Here we demonstrate a method for automatic classification of the Lake Malawi cichlids based on computer vision and geometric morphometrics. For this end we developed a pipeline that integrates multiple image processing tools to automatically extract informative features of color and stripe patterns from a large set of photographic images of wild cichlids. The extracted information was evaluated by statistical classifiers Support Vector Machine and Random Forests. Both classifiers performed better when body shape information was added to the feature of color and stripe. Besides the coloration and stripe pattern, body shape variables boosted the accuracy of classification by about 10%. The programs were able to classify 594 live cichlid individuals belonging to 12 different classes (species and sexes) with an average accuracy of 78%, contrasting to a mere 42% success rate by human eyes. The variables that contributed most to the accuracy were body height and the hue of the most frequent color.

CONCLUSIONS

Computer vision showed a notable performance in extracting information from the color and stripe patterns of Lake Malawi cichlids although the information was not enough for errorless species identification. Our results indicate that there appears an unavoidable difficulty in automatic species identification of cichlid fishes, which may arise from short divergence times and gene flow between closely related species.

Performance of chromID® CARBA medium for carbapenemases-producing Enterobacteriaceae detection during rectal screening

Chromogenic chromID® CARBA medium was compared with CDC method and MacConkey agar with imipenem for its performance in detecting carbapenemase-producing Enterobacteriaceae (CPE) during a faecal screening surveillance program. Double rectal swabs were collected from patients hospitalized in the ICU. One swab was inoculated onto the solid media chromID® CARBA and MacConkey agar with imipenem, while the other was tested according to CDC protocol. Suspected colonies from all procedures were identified to species level and tested for their susceptibility to carbapenems by phenotypic tests. All carbapenem non-susceptible isolates were tested by the modified Hodge test (MHT) and synergy tests. Positive results were confirmed by PCR testing for carbapenemase gene detection. Performance of all three procedures applied was statistically analyzed as compared to MHT and PCR results for the presence of carbapenemase-encoding genes. Out of 177 rectal samples tested, 86 samples were found to contain one or more CPE verified by molecular detection of carbapenemase encoding genes among isolated Enterobacteriaceae. Sensitivity of chromID® CARBA was 96.5 % in clinical samples. Specificity was 91.2 % at the reading level and 100.0 % after Gram staining. chromID® CARBA performed with high accuracy among the phenotypic methods applied, giving early results.

Understanding the origin of species with genome-scale data: modelling gene flow

As it becomes easier to sequence multiple genomes from closely related species, evolutionary biologists working on speciation are struggling to get the most out of very large population genomic data sets. Such data hold the potential to resolve long-standing questions in evolutionary biology about the role of gene exchange in species formation. In principle, the new population genomic data can be used to disentangle the conflicting roles of natural selection and gene flow during the divergence process. However, there are great challenges in taking full advantage of such data, especially with regard to including recombination in genetic models of the divergence process. Current data, models, methods and the potential pitfalls in using them will be considered here.

Diverse captive non-human primates with phytanic acid-deficient diets rich in plant products have substantial phytanic acid levels in their red blood cells

BACKGROUND

Humans and rodents with impaired phytanic acid (PA) metabolism can accumulate toxic stores of PA that have deleterious effects on multiple organ systems. Ruminants and certain fish obtain PA from the microbial degradation of dietary chlorophyll and/or through chlorophyll-derived precursors. In contrast, humans cannot derive PA from chlorophyll and instead normally obtain it only from meat, dairy, and fish products.

RESULTS

Captive apes and Old world monkeys had significantly higher red blood cell (RBC) PA levels relative to humans when all subjects were fed PA-deficient diets. Given the adverse health effects resulting from PA over accumulation, we investigated the molecular evolution of thirteen PA metabolism genes in apes, Old world monkeys, and New world monkeys. All non-human primate (NHP) orthologs are predicted to encode full-length proteins with the marmoset Phyh gene containing a rare, but functional, GA splice donor dinucleotide. Acox2, Scp2, and Pecr sequences had amino acid positions with accelerated substitution rates while Amacr had significant variation in evolutionary rates in apes relative to other primates.

CONCLUSIONS

Unlike humans, diverse captive NHPs with PA-deficient diets rich in plant products have substantial RBC PA levels. The favored hypothesis is that NHPs can derive significant amounts of PA from the degradation of ingested chlorophyll through gut fermentation. If correct, this raises the possibility that RBC PA levels could serve as a biomarker for evaluating the digestive health of captive NHPs. Furthermore, the evolutionary rates of the several genes relevant to PA metabolism provide candidate genetic adaptations to NHP diets.

Origins of shared genetic variation in African cichlids

Cichlid fishes have evolved tremendous morphological and behavioral diversity in the waters of East Africa. Within each of the Great Lakes Tanganyika, Malawi, and Victoria, the phenomena of hybridization and retention of ancestral polymorphism explain allele sharing across species. Here, we explore the sharing of single nucleotide polymorphisms (SNPs) between the major East African cichlid assemblages. A set of approximately 200 genic and nongenic SNPs was ascertained in five Lake Malawi species and genotyped in a diverse collection of ∼160 species from across Africa. We observed segregating polymorphism outside of the Malawi lineage for more than 50% of these loci; this holds similarly for genic versus nongenic SNPs, as well as for SNPs at putative CpG versus non-CpG sites. Bayesian and principal component analyses of genetic structure in the data demonstrate that the Lake Malawi endemic flock is not monophyletic and that river species have likely contributed significantly to Malawi genomes. Coalescent simulations support the hypothesis that river cichlids have transported polymorphism between lake assemblages. We observed strong genetic differentiation between Malawi lineages for approximately 8% of loci, with contributions from both genic and nongenic SNPs. Notably, more than half of these outlier loci between Malawi groups are polymorphic outside of the lake. Cichlid fishes have evolved diversity in Lake Malawi as new mutations combined with standing genetic variation shared across East Africa.

Population genetics and objectivity in species diagnosis

Species as evolutionary lineages are expected to show greater evolutionary independence from one another than are populations within species. Two measures of evolutionary independence that stem from the study of isolation-with-migration models, one reflecting the amount of gene exchange and one reflecting the time of separation, were drawn from the literature for a large number of pairs of closely related species and pairs of populations within species. Both measures, for gene flow and time, showed broadly overlapping distributions for pairs of species and for pairs of populations within species. Species on average show more time and less gene flow than populations, but the similarity of the distributions argues against there being a qualitative difference associated with species status, as compared to populations. The two measures of evolutionary independence were similarly correlated with F(ST) estimates, which in turn also showed similar distributions for species comparisons relative to population comparisons. The measures of gene flow and separation time were examined for the capacity to discriminate intraspecific differences from interspecific differences. If used together, the two measures could be used to develop an objective (in the sense of being repeatable) measure for species diagnosis.

Regarding the confusion between the population concept and Mayr's "population thinking"

Ernst Mayr said that one of Darwin's greatest contributions was to show scholars the way to population thinking, and to help them discard a mindset of typological thinking. Population thinking rejects a focus on a central representative type, and emphasizes the variation among individuals. However, Mayr's choice of terms has led to confusion, particularly among biologists who study natural populations. Both population thinking and the concept of a biological population were inspired by Darwin, and from Darwin the chain for both concepts runs through Francis Galton who introduced the statistical usage of "population" that appears in Mayr's population thinking. It was Galton's "population" that was modified by geneticists and biometricians in the early 20th century to refer to an interbreeding and evolving community of organisms. Under this meaning, a population is a biological entity and so paradoxically population thinking, which emphasizes variation at the expense of dwelling on entities, is usually not about populations. Mayr did not address the potential for misunderstanding but for him the important part of the population concept was that the organisms within a population were variable, and so he probably thought there should not be confusion between population thinking and the concept of a population.

[Radiation-related caries: etiology and possible preventive strategies. What should the radiotherapist know?]

BACKGROUND

After radiation treatment of head-and-neck cancer, the impairment of patient's quality of life still remains an issue. After completion of the treatment course, a substantial number of patients develop so-called radiation caries. In addition, almost 50% of all cases of infectious osteoradionecrosis (iORN) of the jaws are directly associated with radiation caries. This review addresses our current knowledge on the etiology and pathogenesis of radiation caries including possible preventive strategies.

MATERIALS AND METHODS

A PubMed search using the terms "radiation caries" ("radiation related caries", "radiation related damage to dentition") and "radiogenic caries" ("postradiation caries", "dental complications and radiotherapy") was performed. The analysis of its content focused on the etiology, the pathogenesis, and the available knowledge on prophylaxis as well as treatment of radiation caries.

RESULTS

For this review, 60 publications were selected. As main causal factors for radiogenic caries, either indirect impairment, resulting from alterations in the oral environment (e.g., radiation-induced xerostomia) or direct radiation-induced damage in teeth hard tissues are discussed. Radiation caries remains a lifelong threat and, therefore, requires permanent prevention programs.

CONCLUSION

To enable optimal medical care of the patients during the time course of radiotherapy as well as afterwards, close interdisciplinary cooperation between radiotherapists, oral surgeons, otorhinolaryngologists, and dentists is absolutely essential.

Vicariance divergence and gene flow among islet populations of an endemic lizard

Allopatry and allopatric speciation can arise through two different mechanisms: vicariance or colonization through dispersal. Distinguishing between these different allopatric mechanisms is difficult and one of the major challenges in biogeographical research. Here, we address whether allopatric isolation in an endemic island lizard is the result of vicariance or dispersal. We estimated the amount and direction of gene flow during the divergence of isolated islet populations and subspecies of the endemic Skyros wall lizard Podarcis gaigeae, a phenotypically variable species that inhabits a major island and small islets in the Greek archipelago. We applied isolation-with-migration models to estimate population divergence times, population sizes and gene flow between islet-mainland population pairs. Divergence times were significantly correlated with independently estimated geological divergence times. This correlation strongly supports a vicariance scenario where islet populations have sequentially become isolated from the major island. We did not find evidence for significant gene flow within P. g. gaigeae. However, gene-flow estimates from the islet to the mainland populations were positively affected by islet area and negatively by distance between the islet and mainland. We also found evidence for gene flow from one subspecies (P. g. weigandi) into another (P. g. gaigeae), but not in the other direction. Ongoing gene flow between the subspecies suggests that even in this geographically allopatric scenario with the sea posing a strong barrier to dispersal, divergence with some gene flow is still feasible.

Non-equilibrium allele frequency spectra via spectral methods

A major challenge in the analysis of population genomics data consists of isolating signatures of natural selection from background noise caused by random drift and gene flow. Analyses of massive amounts of data from many related populations require high-performance algorithms to determine the likelihood of different demographic scenarios that could have shaped the observed neutral single nucleotide polymorphism (SNP) allele frequency spectrum. In many areas of applied mathematics, Fourier Transforms and Spectral Methods are firmly established tools to analyze spectra of signals and model their dynamics as solutions of certain Partial Differential Equations (PDEs). When spectral methods are applicable, they have excellent error properties and are the fastest possible in high dimension; see Press et al. (2007). In this paper we present an explicit numerical solution, using spectral methods, to the forward Kolmogorov equations for a Wright-Fisher process with migration of K populations, influx of mutations, and multiple population splitting events.

Enrichment of mRNA-like noncoding RNAs in the divergence of Drosophila males

With the advent of transcriptome data, it has become clear that mRNA-like noncoding RNAs (mlncRNAs) are widespread in eukaryotes. Although their functions are poorly understood, these transcripts may play an important role in development and could thus be involved in determining developmental complexity and phenotypic diversification. However, few studies have assessed their potential roles in the divergence of closely related species. Here, we identify and study patterns of sequence and expression divergence in ten novel candidate mlncRNAs from Drosophila pseudoobscura and its close relative D. persimilis. The candidate mlncRNAs were identified by randomly sequencing a group of 734 cDNA clones from a microarray that showed either no difference in expression (187 clones) or differential expression (547 clones) in comparisons between D. pseudoobscura and D. persimilis and between these two species and their F(1) hybrids. Candidate mlncRNAs are overrepresented among differentially expressed transcripts between males of D. pseudoobscura and D. persimilis, and although they have high sequence conservation between these two species, seven of them have no putative homologs in any of the other ten Drosophila species whose genomes have been sequenced. Expression of eight of the ten candidate mlncRNAs was detected either in whole bodies (adults) or testes using a custom-designed oligonucleotide microarray. Three of the ten candidate mlncRNAs are highly expressed (in the top 4% of the male transcriptome), differentially expressed between species, and show extreme levels of sex-bias, with one transcript having the highest level of male bias in the whole transcriptome. Proteomic data from testes show no traces of any predicted peptides from the candidate mlncRNAs. Our results suggest that these mlncRNAs may be important in male-specific processes related to sexual dimorphism and species divergence in this species group.

In defence of model-based inference in phylogeography

Recent papers have promoted the view that model-based methods in general, and those based on Approximate Bayesian Computation (ABC) in particular, are flawed in a number of ways, and are therefore inappropriate for the analysis of phylogeographic data. These papers further argue that Nested Clade Phylogeographic Analysis (NCPA) offers the best approach in statistical phylogeography. In order to remove the confusion and misconceptions introduced by these papers, we justify and explain the reasoning behind model-based inference. We argue that ABC is a statistically valid approach, alongside other computational statistical techniques that have been successfully used to infer parameters and compare models in population genetics. We also examine the NCPA method and highlight numerous deficiencies, either when used with single or multiple loci. We further show that the ages of clades are carelessly used to infer ages of demographic events, that these ages are estimated under a simple model of panmixia and population stationarity but are then used under different and unspecified models to test hypotheses, a usage the invalidates these testing procedures. We conclude by encouraging researchers to study and use model-based inference in population genetics.

The divergence of chimpanzee species and subspecies as revealed in multipopulation isolation-with-migration analyses

The divergence of bonobos and three subspecies of the common chimpanzee was examined under a multipopulation isolation-with-migration (IM) model with data from 73 loci drawn from the literature. A benefit of having a full multipopulation model, relative to conducting multiple pairwise analyses between sampled populations, is that a full model can reveal historical gene flow involving ancestral populations. An example of this was found in which gene flow is indicated between the western common chimpanzee subspecies and the ancestor of the central and the eastern common chimpanzee subspecies. The results of a full analysis on all four populations are strongly consistent with analyses on pairs of populations and generally similar to results from previous studies. The basal split between bonobos and common chimpanzees was estimated at 0.93 Ma (0.68-1.54 Ma, 95% highest posterior density interval), with the split among the ancestor of three common chimpanzee populations at 0.46 Ma (0.35-0.65), and the most recent split between central and eastern common chimpanzee populations at 0.093 Ma (0.041-0.157). Population size estimates mostly fell in the range from 5,000 to 10,000 individuals. The exceptions are the size of the ancestor of the common chimpanzee and the bonobo, at 17,000 (8,000-28,000) individuals, and the central common chimpanzee and its immediate ancestor with the eastern common chimpanzee, which have effective size estimates at 27,000 (16,000-44,000) and 32,000 (19,000-54,000) individuals, respectively.

Isolation with migration models for more than two populations

A method for studying the divergence of multiple closely related populations is described and assessed. The approach of Hey and Nielsen (2007, Integration within the Felsenstein equation for improved Markov chain Monte Carlo methods in population genetics. Proc Natl Acad Sci USA. 104:2785-2790) for fitting an isolation-with-migration model was extended to the case of multiple populations with a known phylogeny. Analysis of simulated data sets reveals the kinds of history that are accessible with a multipopulation analysis. Necessarily, processes associated with older time periods in a phylogeny are more difficult to estimate; and histories with high levels of gene flow are particularly difficult with more than two populations. However, for histories with modest levels of gene flow, or for very large data sets, it is possible to study large complex divergence problems that involve multiple closely related populations or species.

Evolution of Population Structure in a Highly Social Top Predator, the Killer Whale

Intraspecific resource partitioning and social affiliations both have the potential to structure populations, though it is rarely possible to directly assess the impact of these mechanisms on genetic diversity and population divergence. Here, we address this for killer whales (Orcinus orca), which specialize on prey species and hunting strategy and have long-term social affiliations involving both males and females. We used genetic markers to assess the structure and demographic history of regional populations and test the hypothesis that known foraging specializations and matrifocal sociality contributed significantly to the evolution of population structure. We find genetic structure in sympatry between populations of foraging specialists (ecotypes) and evidence for isolation by distance within an ecotype. Fitting of an isolation with migration model suggested ongoing, low-level migration between regional populations (within and between ecotypes) and small effective sizes for extant local populations. The founding of local populations by matrifocal social groups was indicated by the pattern of fixed mtDNA haplotypes in regional populations. Simulations indicate that this occurred within the last 20,000 years (after the last glacial maximum). Our data indicate a key role for social and foraging behavior in the evolution of genetic structure among conspecific populations of the killer whale.

Integration within the Felsenstein equation for improved Markov chain Monte Carlo methods in population genetics

In 1988, Felsenstein described a framework for assessing the likelihood of a genetic data set in which all of the possible genealogical histories of the data are considered, each in proportion to their probability. Although not analytically solvable, several approaches, including Markov chain Monte Carlo methods, have been developed to find approximate solutions. Here, we describe an approach in which Markov chain Monte Carlo simulations are used to integrate over the space of genealogies, whereas other parameters are integrated out analytically. The result is an approximation to the full joint posterior density of the model parameters. For many purposes, this function can be treated as a likelihood, thereby permitting likelihood-based analyses, including likelihood ratio tests of nested models. Several examples, including an application to the divergence of chimpanzee subspecies, are provided.

Divergent haplotypes and human history as revealed in a worldwide survey of X-linked DNA sequence variation

The population genetic history of a 10.1-kbp noncoding region of the human X chromosome was studied using the males of the HGDP-CEPH Human Genome Diversity Panel (672 individuals from 52 populations). The geographic distribution of patterns of variation was roughly consistent with previous studies, with the major exception that 1 highly divergent haplotype (haplotype X, hX) was observed at low frequency in widely scattered non-African populations and not at all observed in sub-Saharan African populations. Microsatellite (short tandem repeat) variation within the sequenced region was low among copies of hX, even though the estimated time of ancestry of hX and other sequences was 1.44 Myr. The estimated age of the common ancestor of all hX copies was 5,230 years (95% consistency index: 2,000-75,480 years). To further address the presence of hX in Africa, additional samples from Chad and Tanzania were screened. Five additional copies of hX were observed, consistent with a history in which hX was present in Africa prior to the migration of modern humans out of Africa and with eastern Africa being the source of non-African modern human populations. Taken together, these features of hX-that it is much older than other haplotypes and uncommon and patchily distributed throughout Africa, Europe, and Asia-present a cautionary tale for interpretations of human history.

On the failure of modern species concepts

The modern age of species concepts began in 1942, when Ernst Mayr gave concept names to several different approaches to species identification. A long list of species concepts then followed, as well as a complex literature on their merits, motivations and uses. Some of these complexities arose as a consequence of the semantic shift that Mayr introduced, in which procedures for identifying species were elevated to concepts. Much of the debate in recent decades over concepts, and over pluralism versus monism, can be seen as an unnecessary consequence of treating species identification criteria as if they were more fundamental concepts. Recently, biologists have begun to recognize both the shortcomings of a lexicon of multiple species concepts and a common evolutionary idea that underlies them.

Sampling from natural populations with RNAI reveals high outcrossing and population structure in Caenorhabditis elegans

Despite a nearly worldwide distribution in nature, Caenorhabditis elegans exhibits low levels of genetic polymorphism, possibly as an indirect consequence of low levels of outcrossing. In the laboratory, Caenorhabditis elegans males are produced at low rates and are steadily eliminated from cultures, so that reproduction happens largely through self-fertilization in hermaphrodites. C. elegans is increasingly the focus of evolutionary research; however, natural outcrossing rates are difficult to measure because mating tests with laboratory strains are usually required to identify C. elegans. We sampled natural populations of C. elegans with an RNA interference (RNAi) assay. Heterozygosities and polymorphism patterns revealed surprisingly high levels of population structure and outcrossing (approximately 22% of individuals are estimated to be the result of outcrossing and not self-fertilization). The finding of strong local population structure, together with low levels of diversity on local and global scales, suggests a metapopulation model of frequent extinction and recolonization of local populations. The occurrence of substantial outcrossing suggests that the extinction of local populations is probably not driven by the accumulation of harmful mutations.

Nuclear gene variation and molecular dating of the cichlid species flock of Lake Malawi

The cichlid fishes of Lake Malawi are famously diverse. However, evolutionary studies have been difficult because of their recent and uncertain phylogenetic history. Portions of 12 nuclear loci were sequenced in nine rock-dwelling species (mbuna) and three representatives of pelagic nonmbuna species. In contrast to the pattern of variation at mitochondrial genes, which do provide phylogenetic resolution at the level of mbuna versus nonmbuna, and among some genera, the nuclear loci were virtually devoid of phylogenetic signal. Only a small minority of variable positions were phylogenetically informative, and no phylogenetic branches are supported by more than one site. From the nuclear gene perspective the Malawian radiation appears to be a star phylogeny, as if the founding of the lake was accompanied by a partial bottleneck. The pattern is different from that found in Lake Victoria, in which nuclear loci share large amounts of ancestral variation. In the case of nuclear genes of Lake Malawi, the absence of phylogenetically informative variation suggests a relative absence of ancestral variation. Nuclear genes also differed from the mitochondria in having nearly twice the amount of divergence from Oreochromis (tilapia). An approximate splitting time between mbuna and nonmbuna lineages was estimated as 0.7 Myr. Oreochromis is estimated to have diverged from the cichlids in Lake Malawi and Lake Tanganyika about 18 MYA.

On the number of New World founders: a population genetic portrait of the peopling of the Americas

The founding of New World populations by Asian peoples is the focus of considerable archaeological and genetic research, and there persist important questions on when and how these events occurred. Genetic data offer great potential for the study of human population history, but there are significant challenges in discerning distinct demographic processes. A new method for the study of diverging populations was applied to questions on the founding and history of Amerind-speaking Native American populations. The model permits estimation of founding population sizes, changes in population size, time of population formation, and gene flow. Analyses of data from nine loci are consistent with the general portrait that has emerged from archaeological and other kinds of evidence. The estimated effective size of the founding population for the New World is fewer than 80 individuals, approximately 1% of the effective size of the estimated ancestral Asian population. By adding a splitting parameter to population divergence models it becomes possible to develop detailed portraits of human demographic history. Analyses of Asian and New World data support a model of a recent founding of the New World by a population of quite small effective size.

A new population-genetic method for assessing human demographic history reveals that the effective size of the founding population of the New World comprised less than 80 individuals.

On the origin of Lake Malawi cichlid species: a population genetic analysis of divergence

The cichlid fishes of Lake Malawi are famously diverse. However, phylogenetic and population genetic studies of their history have been difficult because of the great amount of genetic variation that is shared between species. We apply a recently developed method for fitting the "isolation with migration" divergence model to a data set of specially designed compound loci to develop portraits of cichlid species divergence. Outgroup sequences from a cichlid from Lake Tanganyika permit model parameter estimates in units of years and effective population sizes. Estimated speciation times range from 1,000 to 17,000 years for species in the genus Tropheops. These exceptionally recent dates suggest that Malawi cichlids as a group experience a very active and dynamic diversification process. Current effective population size estimates range form 2,000 to near 40,000, and to >120,000 for estimates of ancestral population sizes. It appears that very recent speciation and gene flow are among the reasons why it has been difficult to discern the phylogenetic history of Malawi cichlids.

Multilocus methods for estimating population sizes, migration rates and divergence time, with applications to the divergence of Drosophila pseudoobscura and D. persimilis

The genetic study of diverging, closely related populations is required for basic questions on demography and speciation, as well as for biodiversity and conservation research. However, it is often unclear whether divergence is due simply to separation or whether populations have also experienced gene flow. These questions can be addressed with a full model of population separation with gene flow, by applying a Markov chain Monte Carlo method for estimating the posterior probability distribution of model parameters. We have generalized this method and made it applicable to data from multiple unlinked loci. These loci can vary in their modes of inheritance, and inheritance scalars can be implemented either as constants or as parameters to be estimated. By treating inheritance scalars as parameters it is also possible to address variation among loci in the impact via linkage of recurrent selective sweeps or background selection. These methods are applied to a large multilocus data set from Drosophila pseudoobscura and D. persimilis. The species are estimated to have diverged approximately 500,000 years ago. Several loci have nonzero estimates of gene flow since the initial separation of the species, with considerable variation in gene flow estimates among loci, in both directions between the species.