Phase-type distributions in mathematical population genetics: An emerging framework.

A phase-type distribution is the time to absorption in a continuous- or discrete-time Markov chain. Phase-type distributions can be used as a general framework to calculate key properties of the standard coalescent model and many of its extensions. Here, the 'phases' in the phase-type distribution correspond to states in the ancestral process. For example, the time to the most recent common ancestor and the total branch length are phase-type distributed. Furthermore, the site frequency spectrum follows a multivariate discrete phase-type distribution and the joint distribution of total branch lengths in the two-locus coalescent-with-recombination model is multivariate phase-type distributed. In general, phase-type distributions provide a powerful mathematical framework for coalescent theory because they are analytically tractable using matrix manipulations. The purpose of this review is to explain the phase-type theory and demonstrate how the theory can be applied to derive basic properties of coalescent models. These properties can then be used to obtain insight into the ancestral process, or they can be applied for statistical inference. In particular, we show the relation between classical first-step analysis of coalescent models and phase-type calculations. We also show how reward transformations in phase-type theory lead to easy calculation of covariances and correlation coefficients between e.g. tree height, tree length, external branch length, and internal branch length. Furthermore, we discuss how these quantities can be used for statistical inference based on estimating equations. Providing an alternative to previous work based on the Laplace transform, we derive likelihoods for small-size coalescent trees based on phase-type theory. Overall, our main aim is to demonstrate that phase-type distributions provide a convenient general set of tools to understand aspects of coalescent models that are otherwise difficult to derive. Throughout the review, we emphasize the versatility of the phase-type framework, which is also illustrated by our accompanying R-code. All our analyses and figures can be reproduced from code available on GitHub.

Haplotype based testing for a better understanding of the selective architecture.

BACKGROUND: The identification of genomic regions affected by selection is one of the most important goals in population genetics. If temporal data are available, allele frequency changes at SNP positions are often used for this purpose. Here we provide a new testing approach that uses haplotype frequencies instead of allele frequencies. RESULTS: Using simulated data, we show that compared to SNP based test, our approach has higher power, especially when the number of candidate haplotypes is small or moderate. To improve power when the number of haplotypes is large, we investigate methods to combine them with a moderate number of haplotype subsets. Haplotype frequencies can often be recovered with less noise than SNP frequencies, especially under pool sequencing, giving our test an additional advantage. Furthermore, spurious outlier SNPs may lead to false positives, a problem usually not encountered when working with haplotypes. Post hoc tests for the number of selected haplotypes and for differences between their selection coefficients are also provided for a better understanding of the underlying selection dynamics. An application on a real data set further illustrates the performance benefits. CONCLUSIONS: Due to less multiple testing correction and noise reduction, haplotype based testing is able to outperform SNP based tests in terms of power in most scenarios.

Quantifying Selection with Pool-Seq Time Series Data.

Allele frequency time series data constitute a powerful resource for unraveling mechanisms of adaptation, because the temporal dimension captures important information about evolutionary forces. In particular, Evolve and Resequence (E&R), the whole-genome sequencing of replicated experimentally evolving populations, is becoming increasingly popular. Based on computer simulations several studies proposed experimental parameters to optimize the identification of the selection targets. No such recommendations are available for the underlying parameters selection strength and dominance. Here, we introduce a highly accurate method to estimate selection parameters from replicated time series data, which is fast enough to be applied on a genome scale. Using this new method, we evaluate how experimental parameters can be optimized to obtain the most reliable estimates for selection parameters. We show that the effective population size (Ne) and the number of replicates have the largest impact. Because the number of time points and sequencing coverage had only a minor effect, we suggest that time series analysis is feasible without major increase in sequencing costs. We anticipate that time series analysis will become routine in E&R studies.

Approximate maximum likelihood estimation for population genetic inference.

In many population genetic problems, parameter estimation is obstructed by an intractable likelihood function. Therefore, approximate estimation methods have been developed, and with growing computational power, sampling-based methods became popular. However, these methods such as Approximate Bayesian Computation (ABC) can be inefficient in high-dimensional problems. This led to the development of more sophisticated iterative estimation methods like particle filters. Here, we propose an alternative approach that is based on stochastic approximation. By moving along a simulated gradient or ascent direction, the algorithm produces a sequence of estimates that eventually converges to the maximum likelihood estimate, given a set of observed summary statistics. This strategy does not sample much from low-likelihood regions of the parameter space, and is fast, even when many summary statistics are involved. We put considerable efforts into providing tuning guidelines that improve the robustness and lead to good performance on problems with high-dimensional summary statistics and a low signal-to-noise ratio. We then investigate the performance of our resulting approach and study its properties in simulations. Finally, we re-estimate parameters describing the demographic history of Bornean and Sumatran orang-utans.

The long zinc finger domain of PRDM9 forms a highly stable and long-lived complex with its DNA recognition sequence.

PR domain containing protein 9 (PRDM9) is a meiosis-specific, multi-domain protein that regulates the location of recombination hotspots by targeting its DNA recognition sequence for double-strand breaks (DSBs). PRDM9 specifically recognizes DNA via its tandem array of zinc fingers (ZnFs), epigenetically marks the local chromatin by its histone methyltransferase activity, and is an important tether that brings the DNA into contact with the recombination initiation machinery. A strong correlation between PRDM9-ZnF variants and specific DNA motifs at recombination hotspots has been reported; however, the binding specificity and kinetics of the ZnF domain are still obscure. Using two in vitro methods, gel mobility shift assays and switchSENSE, a quantitative biophysical approach that measures binding rates in real time, we determined that the PRDM9-ZnF domain forms a highly stable and long-lived complex with its recognition sequence, with a dissociation halftime of many hours. The ZnF domain exhibits an equilibrium dissociation constant (K D) in the nanomolar (nM) range, with polymorphisms in the recognition sequence directly affecting the binding affinity. We also determined that alternative sequences (15-16 nucleotides in length) can be specifically bound by different subsets of the ZnF domain, explaining the binding plasticity of PRDM9 for different sequences. Finally, longer binding targets are preferred than predicted from the numbers of ZnFs contacting the DNA. Functionally, a long-lived complex translates into an enzymatically active PRDM9 at specific DNA-binding sites throughout meiotic prophase I that might be relevant in stabilizing the components of the recombination machinery to a specific DNA target until DSBs are initiated by Spo11.

A genome-wide, fine-scale map of natural pigmentation variation in Drosophila melanogaster.

Various approaches can be applied to uncover the genetic basis of natural phenotypic variation, each with their specific strengths and limitations. Here, we use a replicated genome-wide association approach (Pool-GWAS) to fine-scale map genomic regions contributing to natural variation in female abdominal pigmentation in Drosophila melanogaster, a trait that is highly variable in natural populations and highly heritable in the laboratory. We examined abdominal pigmentation phenotypes in approximately 8000 female European D. melanogaster, isolating 1000 individuals with extreme phenotypes. We then used whole-genome Illumina sequencing to identify single nucleotide polymorphisms (SNPs) segregating in our sample, and tested these for associations with pigmentation by contrasting allele frequencies between replicate pools of light and dark individuals. We identify two small regions near the pigmentation genes tan and bric-à-brac 1, both corresponding to known cis-regulatory regions, which contain SNPs showing significant associations with pigmentation variation. While the Pool-GWAS approach suffers some limitations, its cost advantage facilitates replication and it can be applied to any non-model system with an available reference genome.

Multiscale DNA partitioning: statistical evidence for segments.

MOTIVATION: DNA segmentation, i.e. the partitioning of DNA in compositionally homogeneous segments, is a basic task in bioinformatics. Different algorithms have been proposed for various partitioning criteria such as Guanine/Cytosine (GC) content, local ancestry in population genetics or copy number variation. A critical component of any such method is the choice of an appropriate number of segments. Some methods use model selection criteria and do not provide a suitable error control. Other methods that are based on simulating a statistic under a null model provide suitable error control only if the correct null model is chosen. RESULTS: Here, we focus on partitioning with respect to GC content and propose a new approach that provides statistical error control: as in statistical hypothesis testing, it guarantees with a user-specified probability [Formula: see text] that the number of identified segments does not exceed the number of actually present segments. The method is based on a statistical multiscale criterion, rendering this as a segmentation method that searches segments of any length (on all scales) simultaneously. It is also accurate in localizing segments: under benchmark scenarios, our approach leads to a segmentation that is more accurate than the approaches discussed in the comparative review of Elhaik et al. In our real data examples, we find segments that often correspond well to features taken from standard University of California at Santa Cruz (UCSC) genome annotation tracks. AVAILABILITY AND IMPLEMENTATION: Our method is implemented in function smuceR of the R-package stepR available at http://www.stochastik.math.uni-goettingen.de/smuce.

Estimating Haplotype Structure and Frequencies: A Bayesian Approach to Unknown Design in Pooled Genomic Data.

The estimation of haplotype structure and frequencies provides crucial information about the composition of genomes. Techniques, such as single-individual haplotyping, aim to reconstruct individual haplotypes from diploid genome sequencing data. However, our focus is distinct. We address the challenge of reconstructing haplotype structure and frequencies from pooled sequencing samples where multiple individuals are sequenced simultaneously. A frequentist method to address this issue has recently been proposed. In contrast to this and other methods that compute point estimates, our proposed Bayesian hierarchical model delivers a posterior that permits us to also quantify uncertainty. Since matching permutations in both haplotype structure and corresponding frequency matrix lead to the same reconstruction of their product, we introduce an order-preserving shrinkage prior that ensures identifiability with respect to permutations. For inference, we introduce a blocked Gibbs sampler that enforces the required constraints. In a simulation study, we assessed the performance of our method. Furthermore, by using our approach on two distinct sets of real data, we demonstrate that our Bayesian approach can reconstruct the dominant haplotypes in a challenging, high-dimensional set-up.

Detecting selective sweeps from pooled next-generation sequencing samples.

Due to its cost effectiveness, next-generation sequencing of pools of individuals (Pool-Seq) is becoming a popular strategy for characterizing variation in population samples. Because Pool-Seq provides genome-wide SNP frequency data, it is possible to use them for demographic inference and/or the identification of selective sweeps. Here, we introduce a statistical method that is designed to detect selective sweeps from pooled data by accounting for statistical challenges associated with Pool-Seq, namely sequencing errors and random sampling among chromosomes. This allows for an efficient use of the information: all base calls are included in the analysis, but the higher credibility of regions with higher coverage and base calls with better quality scores is accounted for. Computer simulations show that our method efficiently detects sweeps even at very low coverage (0.5× per chromosome). Indeed, the power of detecting sweeps is similar to what we could expect from sequences of individual chromosomes. Since the inference of selective sweeps is based on the allele frequency spectrum (AFS), we also provide a method to accurately estimate the AFS provided that the quality scores for the sequence reads are reliable. Applying our approach to Pool-Seq data from Drosophila melanogaster, we identify several selective sweep signatures on chromosome X that include some previously well-characterized sweeps like the wapl region.

DNA pooling and statistical tests for the detection of single nucleotide polymorphisms.

The development of next generation genome sequencers gives the opportunity of learning more about the genetic make-up of human and other populations. One important question involves the location of sites at which variation occurs within a population. Our focus will be on the detection of rare variants. Such variants will often not be present in smaller samples and are hard to distinguish from sequencing errors in larger samples. This is particularly true for pooled samples which are often used as part of a cost saving strategy. The focus of this article is on experiments that involve DNA pooling. We derive experimental designs that optimize the power of statistical tests for detecting single nucleotide polymorphisms (SNPs, sites at which there is variation within a population). We also present a new simple test that calls a SNP, if the maximum number of reads of a prospective variant across lanes exceeds a certain threshold. The value of this threshold is defined according to the number of available lanes, the parameters of the genome sequencer and a specified probability of accepting that there is variation at a site when no variation is present. On the basis of this test, we derive pool sizes which are optimal for the detection of rare variants. This test is compared with a likelihood ratio test, which takes into account the number of reads of a prospective variant from all the lanes. It is shown that the threshold based rule achieves a comparable power to this likelihood ratio test and may well be a useful tool in determining near optimal pool sizes for the detection of rare alleles in practical applications.

Development and evaluation of an animal health and welfare monitoring system for veterinary supervision of pullet farms.

Regular welfare monitoring throughout rearing of pullets may help to identify problems early and take counteractions timely, which helps in guaranteeing good welfare. The aims of our observational study were (i) to establish and test a welfare monitoring system that can be used during (short) routine veterinary and technical staff visits for pullet flocks, (ii) to use the monitoring system to investigate variability between flocks and (iii) to analyse factors that potentially affect pullets' body weight, uniformity in body weight and mortality. The developed monitoring system tries to minimise the time required while not losing important information. Age-specific recording sheets comprise animal-based indicators of welfare and relevant environmental factors (housing, management, care) to allow for identifying causes of problems and targeted action. Finally, the system was implemented in a cross-sectional study and data collected in 100 flocks (67 organic, 33 conventional) on 28 rearing farms in Austria. Linear mixed models were used to identify factors influencing body weight, uniformity and mortality, both including all flocks (A) and only organic flocks (O) and a linear regression model with all flocks to investigate associations within animal-based indicators. High variability was found between flocks in animal-based indicators. Body weight was higher when the pre-rearing period was shorter (p ≤ 0.001, A&O), with higher intensities of light (p = 0.012, O), with only one compared to more stockpersons (p ≤ 0.007, A&O), with a higher number of flock visits per day (p ≤ 0.018, A&O), and a lower avoidance distance (p = 0.034, A). Body weight uniformity increased, with age and decreased with the duration of the light period (p = 0.046, A), and, amongst others, was higher on organic farms (farming type; p = 0.041). The latter may reflect a more uniform level of welfare due to a lower stocking density and lowered effects of social competition. Within organic flocks mortality was lower if pullets had access to a covered veranda (p = 0.025) resulting in an overall lower stocking density inside the barn, while in the model including all farms mortality was higher in cases where a disease had been diagnosed. We conclude that our monitoring system can easily be implemented in regular veterinary and technical staff visits, but could also be used by the farmers'. Several easy-to-record animal-based indicators of animal welfare could be analysed more frequently to increase early detection of problems. Implementation of such a routine-based monitoring system with easy-to-assess animal-based parameters and input measures can contribute to better animal health and welfare in pullets.

Using fecal immmunochemical cartridges for gut microbiome analysis within a colorectal cancer screening program.

The colorectal cancer (CRC) screening program B-PREDICT is an invited two-stage screening project using a fecal immunochemical test (FIT) for initial screening followed by a colonoscopy for those with a positive FIT. Since the gut microbiome likely plays a role in the etiology of CRC, microbiome-based biomarkers in combination with FIT could be a promising tool for optimizing CRC screening. Therefore, we evaluated the usability of FIT cartridges for microbiome analysis and compared it to Stool Collection and Preservation Tubes. Corresponding FIT cartridges as well as Stool Collection and Preservation Tubes were collected from participants of the B-PREDICT screening program to perform 16S rRNA gene sequencing. We calculated intraclass correlation coefficients (ICCs) based on center log ratio transformed abundances and used ALDEx2 to test for significantly differential abundant taxa between the two sample types. Additionally, FIT and Stool Collection and Preservation Tube triplicate samples were obtained from volunteers to estimate variance components of microbial abundances. FIT and Preservation Tube samples produce highly similar microbiome profiles which cluster according to subject. Significant differences between the two sample types can be found for abundances of some bacterial taxa (e.g. 33 genera) but are minor compared to the differences between the subjects. Analysis of triplicate samples revealed slightly worse repeatability of results for FIT than for Preservation Tube samples. Our findings indicate that FIT cartridges are appropriate for gut microbiome analysis nested within CRC screening programs.

Waves Out of the Korean Peninsula and Inter- and Intra-Species Replacements in Freshwater Fishes in Japan.

The Japanese archipelago is located at the periphery of the continent of Asia. Rivers in the Japanese archipelago, separated from the continent of Asia by about 17 Ma, have experienced an intermittent exchange of freshwater fish taxa through a narrow land bridge generated by lowered sea level. As the Korean Peninsula and Japanese archipelago were not covered by an ice sheet during glacial periods, phylogeographical analyses in this region can trace the history of biota that were, for a long time, beyond the last glacial maximum. In this study, we analyzed the phylogeography of four freshwater fish taxa, Hemibarbus longirostris, dark chub Nipponocypris temminckii, Tanakia ssp. and Carassius ssp., whose distributions include both the Korean Peninsula and Western Japan. We found for each taxon that a small component of diverse Korean clades of freshwater fishes migrated in waves into the Japanese archipelago to form the current phylogeographic structure of biota. The replacements of indigenous populations by succeeding migrants may have also influenced the phylogeography.

Multiple haplotype reconstruction from allele frequency data.

Because haplotype information is of widespread interest in biomedical applications, effort has been put into their reconstruction. Here, we propose an efficient method, called haploSep, that is able to accurately infer major haplotypes and their frequencies just from multiple samples of allele frequency data. Even the accuracy of experimentally obtained allele frequencies can be improved by re-estimating them from our reconstructed haplotypes. From a methodological point of view, we model our problem as a multivariate regression problem where both the design matrix and the coefficient matrix are unknown. Compared to other methods, haploSep is very fast, with linear computational complexity in the haplotype length. We illustrate our method on simulated and real data focusing on experimental evolution and microbial data.

Detecting selective sweeps: a new approach based on hidden markov models.

Detecting and localizing selective sweeps on the basis of SNP data has recently received considerable attention. Here we introduce the use of hidden Markov models (HMMs) for the detection of selective sweeps in DNA sequences. Like previously published methods, our HMMs use the site frequency spectrum, and the spatial pattern of diversity along the sequence, to identify selection. In contrast to earlier approaches, our HMMs explicitly model the correlation structure between linked sites. The detection power of our methods, and their accuracy for estimating the selected site location, is similar to that of competing methods for constant size populations. In the case of population bottlenecks, however, our methods frequently showed fewer false positives.

Talking to Cows: Reactions to Different Auditory Stimuli During Gentle Human-Animal Interactions.

The quality of the animal-human relationship and, consequently, the welfare of animals can be improved by gentle interactions such as stroking and talking. The perception of different stimuli during these interactions likely plays a key role in their emotional experience, but studies are scarce. During experiments, the standardization of verbal stimuli could be increased by using a recording. However, the use of a playback might influence the perception differently than "live" talking, which is closer to on-farm practice. Thus, we compared heifers' (n = 28) reactions to stroking while an experimenter was talking soothingly ("live") or while a recording of the experimenter talking soothingly was played ("playback"). Each animal was tested three times per condition and each trial comprised three phases: pre-stimulus, stimulus (stroking and talking) and post-stimulus. In both conditions, similar phrases with positive content were spoken calmly, using long low-pitched vowels. All tests were video recorded and analyzed for behaviors associated with different affective states. Effects on the heifers' cardiac parameters were assessed using analysis of heart rate variability. Independently of the auditory stimuli, longer durations of neck stretching occurred during stroking, supporting our hypothesis of a positive perception of stroking. Observation of ear positions revealed longer durations of the "back up" position and less ear flicking and changes of ear positions during stroking. The predicted decrease in HR during stroking was not confirmed; instead we found a slightly increased mean HR during stroking with a subsequent decrease in HR, which was stronger after stroking with live talking. In combination with differences in HRV parameters, our findings suggest that live talking might have been more pleasurable to the animals and had a stronger relaxing effect than "playback." The results regarding the effects of the degree of standardization of the stimulus on the variability of the data were inconclusive. We thus conclude that the use of recorded auditory stimuli to promote positive affective states during human-animal interactions in experimental settings is possible, but not necessarily preferable.

Gentle interactions with restrained and free-moving cows: Effects on the improvement of the animal-human relationship.

The animal-human relationship is essential for farm animal welfare and production. Generally, gentle tactile and vocal interactions improve the animal-human relationship in cattle. However, cows that are fearful of humans avoid their close presence and touch; thus, the animal-human relationship first has to be improved to a point where the animals accept stroking before their perception of the interactions and consequently the animal-human relationship can become positive. We tested whether the animal-human relationship of cows fearful of humans is improved more effectively by gentle interactions during restraint, allowing physical contact from the beginning, or if the gentle interactions are offered while the animals are free to move, giving them more control over the situation and thus probably a higher level of agency and a more positive perception of the interactions. Thirty-six dairy cows (median avoidance distance 1.6 m) were assigned to three treatments (each n = 12): gentle vocal and tactile interactions during restraint in the feeding rack (LOCK); gentle vocal and, if possible, tactile interactions while free in the barn (FREE); routine management without additional interactions (CON). Treatments were applied for 3 min per cow on 10 d per fortnight for 6 weeks (i.e., three periods). Avoidance and approach behaviour towards humans was tested before the start of the treatment period, and then at 2-week intervals. The recorded variables were reduced to one score by Principal Component Analysis. The resulting relationship score (higher values implying a better relationship with humans) increased in all groups; the increase was stronger in FREE than in CON, with the increase in LOCK being not significantly different from the other treatment groups. Thus, we recommend that gentle interactions with cows should take place while they are unrestrained, if possible.

Effects of Different Stroking Styles on Behaviour and Cardiac Parameters in Heifers.

Gentle animal-human interactions, such as stroking, can promote positive emotions and thus welfare in cattle. While previous studies showed that stroking at the ventral neck elicited the most positive reactions in cows, intra-specific allogrooming in cattle includes different body regions and is probably guided partly by the receiver. Thus, we compared heifers' (n = 28) reactions to stroking with the experimenter either reactively responding to perceived momentary preferences of the heifers or exclusively stroking the ventral neck. Independently of the stroking style, longer durations of neck stretching and contact occurred during stroking, supporting our hypothesis of a positive perception of stroking. We did not confirm the predicted decrease in heart rate and increase in heart rate variability, but instead found a slightly increased mean heart rate during stroking. The different stroking styles elicited differences in the heifers' ear positions: "reactive" stroking led to longer durations of low ear positions during stroking, while during "ventral neck" stroking, the duration of back up increased. However, no other behaviours differed significantly between different stroking styles, indicating that the exact manner of stroking applied in our treatments seemed to be less important in the promotion of positive affective states in cattle through gentle human-animal interactions.

An omnibus test for the global null hypothesis.

Global hypothesis tests are a useful tool in the context of clinical trials, genetic studies, or meta-analyses, when researchers are not interested in testing individual hypotheses, but in testing whether none of the hypotheses is false. There are several possibilities how to test the global null hypothesis when the individual null hypotheses are independent. If it is assumed that many of the individual null hypotheses are false, combination tests have been recommended to maximize power. If, however, it is assumed that only one or a few null hypotheses are false, global tests based on individual test statistics are more powerful (e.g. Bonferroni or Simes test). However, usually there is no a priori knowledge on the number of false individual null hypotheses. We therefore propose an omnibus test based on cumulative sums of the transformed p-values. We show that this test yields an impressive overall performance. The proposed method is implemented in an R-package called omnibus.

The Effects of Play Behavior, Feeding, and Time of Day on Salivary Concentrations of sIgA in Calves.

The focus of animal welfare science has shifted over the last decades from efforts to avoid negative states to ways of allowing animals the experience of positive emotions. They may influence physiological processes in farmed animals, potentially providing health benefits; in addition, the physiological changes might be used as indicators of emotional states. We investigated calves' salivary secretory immunoglobulin A (sIgA) concentrations with regard to a possible circadian rhythm and two situations that elicit positive emotions. Ten saliva samples of 14 calves were taken on two consecutive days; within the course of a day we observed a significant decline in salivary sIgA concentrations at 14:00 h. Further, we probed the animals before and after milk feeding and, contrarily to our prediction, detected lower sIgA concentrations 5 min after feeding than 15 min before. A probable explanation might be an increase in salivary flow rate caused by milk ingestion. We also took samples before and after we stimulated play behavior in calves. There was no significant difference in sIgA concentrations between samples taken before and after play. Although there was a significant correlation between the change in sIgA concentrations and the amount of play behavior shown, the correlation depended on an unexpected decrease of sIgA in animals that played little, and thus, does not support our hypothesis. In general, the data showed a large variability that might arise from different factors that are difficult to standardize in animals. Thus, the use of salivary sIgA concentrations as a marker of positive emotions in calves is not supported conclusively by the present data.