Integrated Analyses of Microbiome and Longitudinal Metabolome Data Reveal Microbial-Host Interactions on Sulfur Metabolism in Parkinson's Disease.

Parkinson's disease (PD) exhibits systemic effects on the human metabolism, with emerging roles for the gut microbiome. Here, we integrate longitudinal metabolome data from 30 drug-naive, de novo PD patients and 30 matched controls with constraint-based modeling of gut microbial communities derived from an independent, drug-naive PD cohort, and prospective data from the general population. Our key results are (1) longitudinal trajectory of metabolites associated with the interconversion of methionine and cysteine via cystathionine differed between PD patients and controls; (2) dopaminergic medication showed strong lipidomic signatures; (3) taurine-conjugated bile acids correlated with the severity of motor symptoms, while low levels of sulfated taurolithocholate were associated with PD incidence in the general population; and (4) computational modeling predicted changes in sulfur metabolism, driven by A. muciniphila and B. wadsworthia, which is consistent with the changed metabolome. The multi-omics integration reveals PD-specific patterns in microbial-host sulfur co-metabolism that may contribute to PD severity.

Sweepstake evolution revealed by population-genetic analysis of copy-number alterations in single genomes of breast cancer.

Single-cell sequencing is a promising technology that can address cancer cell evolution by identifying genetic alterations in individual cells. In a recent study, genome-wide DNA copy numbers of single cells were accurately quantified by single-cell sequencing in breast cancers. Phylogenetic-tree analysis revealed genetically distinct populations, each consisting of homogeneous cells. Bioinformatics methods based on population genetics should be further developed to quantitatively analyse the single-cell sequencing data. We developed a bioinformatics framework that was combined with molecular-evolution theories to analyse copy-number losses. This analysis revealed that most deletions in the breast cancers at the single-cell level were generated by simple stochastic processes. A non-standard type of coalescent theory, the multiple-merger coalescent model, aided by approximate Bayesian computation fit well with the data, allowing us to estimate the population-genetic parameters in addition to false-positive and false-negative rates. The estimated parameters suggest that the cancer cells underwent sweepstake evolution, where only one or very few parental cells produced a descendent cell population. We conclude that breast cancer cells successively substitute in a tumour mass, and the high reproduction of only a portion of cancer cells may confer high adaptability to this cancer.

Genome wide association study of SNP-, gene-, and pathway-based approaches to identify genes influencing susceptibility to Staphylococcus aureus infections.

BACKGROUND: We conducted a genome-wide association study (GWAS) to identify specific genetic variants that underlie susceptibility to diseases caused by Staphylococcus aureus in humans. METHODS: Cases (n = 309) and controls (n = 2925) were genotyped at 508,921 single nucleotide polymorphisms (SNPs). Cases had at least one laboratory and clinician confirmed disease caused by S. aureus whereas controls did not. R-package (for SNP association), EIGENSOFT (to estimate and adjust for population stratification) and gene- (VEGAS) and pathway-based (DAVID, PANTHER, and Ingenuity Pathway Analysis) analyses were performed. RESULTS: No SNP reached genome-wide significance. Four SNPs exceeded the p < 10(-5) threshold including two (rs2455012 and rs7152530) reaching a p-value < 10(-7). The nearby genes were PDE4B (rs2455012), TXNRD2 (rs3804047), VRK1 and BCL11B (rs7152530), and PNPLA5 (rs470093). The top two findings from the gene-based analysis were NMRK2 (p gene = 1.20E-05), which codes an integrin binding molecule (focal adhesion), and DAPK3 (p gene = 5.10E-05), a serine/threonine kinase (apoptosis and cytokinesis). The pathway analyses identified epithelial cell responses to mechanical and non-mechanical stress. CONCLUSION: We identified potential susceptibility genes for S. aureus diseases in this preliminary study but confirmation by other studies is needed. The observed associations could be relevant given the complexity of S. aureus as a pathogen and its ability to exploit multiple biological pathways to cause infections in humans.

A novel analytical method, Birth Date Selection Mapping, detects response of the Angus (Bos taurus) genome to selection on complex traits.

BACKGROUND: Several methods have recently been developed to identify regions of the genome that have been exposed to strong selection. However, recent theoretical and empirical work suggests that polygenic models are required to identify the genomic regions that are more moderately responding to ongoing selection on complex traits. We examine the effects of multi-trait selection on the genome of a population of US registered Angus beef cattle born over a 50-year period representing approximately 10 generations of selection. We present results from the application of a quantitative genetic model, called Birth Date Selection Mapping, to identify signatures of recent ongoing selection. RESULTS: We show that US Angus cattle have been systematically selected to alter their mean additive genetic merit for most of the 16 production traits routinely recorded by breeders. Using Birth Date Selection Mapping, we estimate the time-dependency of allele frequency for 44,817 SNP loci using genomic best linear unbiased prediction, generalized least squares, and BayesCπ analyses. Finally, we reconstruct the primary phenotypes that have historically been exposed to selection from a genome-wide analysis of the 16 production traits and gene ontology enrichment analysis. CONCLUSIONS: We demonstrate that Birth Date Selection Mapping utilizing mixed models corrects for time-dependent pedigree sampling effects that lead to spurious SNP associations and reveals genomic signatures of ongoing selection on complex traits. Because multiple traits have historically been selected in concert and most quantitative trait loci have small effects, selection has incrementally altered allele frequencies throughout the genome. Two quantitative trait loci of large effect were not the most strongly selected of the loci due to their antagonistic pleiotropic effects on strongly selected phenotypes. Birth Date Selection Mapping may readily be extended to temporally-stratified human or model organism populations.

Inference for ecological dynamical systems: a case study of two endemic diseases.

A Bayesian Markov chain Monte Carlo method is used to infer parameters for an open stochastic epidemiological modEL: the Markovian susceptible-infected-recovered (SIR) model, which is suitable for modeling and simulating recurrent epidemics. This allows exploring two major problems of inference appearing in many mechanistic population models. First, trajectories of these processes are often only partly observed. For example, during an epidemic the transmission process is only partly observable: one cannot record infection times. Therefore, one only records cases (infections) as the observations. As a result some means of imputing or reconstructing individuals in the susceptible cases class must be accomplished. Second, the official reporting of observations (cases in epidemiology) is typically done not as they are actually recorded but at some temporal interval over which they have been aggregated. To address these issues, this paper investigates the following problems. Parameter inference for a perfectly sampled open Markovian SIR is first considered. Next inference for an imperfectly observed sample path of the system is studied. Although this second problem has been solved for the case of closed epidemics, it has proven quite difficult for the case of open recurrent epidemics. Lastly, application of the statistical theory is made to measles and pertussis epidemic time series data from 60 UK cities.

A fast and reliable computational method for estimating population genetic parameters.

The estimation of ancestral and current effective population sizes in expanding populations is a fundamental problem in population genetics. Recently it has become possible to scan entire genomes of several individuals within a population. These genomic data sets can be used to estimate basic population parameters such as the effective population size and population growth rate. Full-data-likelihood methods potentially offer a powerful statistical framework for inferring population genetic parameters. However, for large data sets, computationally intensive methods based upon full-likelihood estimates may encounter difficulties. First, the computational method may be prohibitively slow or difficult to implement for large data. Second, estimation bias may markedly affect the accuracy and reliability of parameter estimates, as suggested from past work on coalescent methods. To address these problems, a fast and computationally efficient least-squares method for estimating population parameters from genomic data is presented here. Instead of modeling genomic data using a full likelihood, this new approach uses an analogous function, in which the full data are replaced with a vector of summary statistics. Furthermore, these least-squares estimators may show significantly less estimation bias for growth rate and genetic diversity than a corresponding maximum-likelihood estimator for the same coalescent process. The least-squares statistics also scale up to genome-sized data sets with many nucleotides and loci. These results demonstrate that least-squares statistics will likely prove useful for nonlinear parameter estimation when the underlying population genomic processes have complex evolutionary dynamics involving interactions between mutation, selection, demography, and recombination.

Tracking the dynamics of pathogen interactions: modeling ecological and immune-mediated processes in a two-pathogen single-host system.

Traditionally, epidemiological studies have focused on understanding the dynamics of a single pathogen, assuming no interactions with other pathogens. Recently, a large body of work has begun to explore the effects of immune-mediated interactions, arising from cross-immunity and antibody-dependent enhancement, between related pathogen strains. In addition, ecological processes such as a temporary period of convalescence and pathogen-induced mortality have led to the concept of ecological interference between unrelated diseases. There remains, however, the need for a systematic study of both immunological and ecological processes within a single framework. In this paper, we develop a general two-pathogen single-host model of pathogen interactions that simultaneously incorporates these mechanisms. We are then able to mechanistically explore how immunoecological processes mediate interactions between diseases for a pool of susceptible individuals. We show that the precise nature of the interaction can induce either competitive or cooperative associations between pathogens. Understanding the dynamic implications of multi-pathogen associations has potentially important public health consequences. Such a framework may be especially helpful in disentangling the effects of partially cross-immunizing infections that affect populations with a pre-disposition towards immunosuppression such as children and the elderly.