Modeling DNA methylation dynamics with approaches from phylogenetics

Modeling DNA methylation dynamics with approaches from phylogenetics

John A. Capra, Dennis Kostka
(Submitted on 11 Apr 2014)

Methylation of CpG dinucleotides is a prevalent epigenetic modification that is required for proper development in vertebrates, and changes in CpG methylation are essential to cellular differentiation. Genome-wide DNA methylation assays have become increasingly common, and recently distinct stages across differentiating cellular lineages have been assayed. How- ever, current methods for modeling methylation dynamics do not account for the dependency structure between precursor and dependent cell types. We developed a continuous-time Markov chain approach, based on the observation that changes in methylation state over tissue differentiation can be modeled similarly to DNA nucleotide changes over evolutionary time. This model explicitly takes precursor to descendant relationships into account and enables inference of CpG methylation dynamics. To illustrate our method, we analyzed a high-resolution methylation map of the differentiation of mouse stem cells into several blood cell types. Our model can successfully infer unobserved CpG methylation states from observations at the same sites in related cell types (90% correct), and this approach more accurately reconstructs missing data than imputation based on neighboring CpGs (84% correct). Additionally, the single CpG resolution of our methylation dynamics estimates enabled us to show that DNA sequence context of CpG sites is informative about methylation dynamics across tissue differentiation. Finally, we identified genomic regions with clusters of highly dynamic CpGs and present a likely functional example. Our work establishes a framework for inference and modeling that is well-suited to DNA methylation data, and our success suggests that other methods for analyzing DNA nucleotide substitutions will also translate to the modeling of epigenetic phenomena.

Flexible methods for estimating genetic distances from nucleotide data

Flexible methods for estimating genetic distances from nucleotide data

Simon Joly, David J Bryant, Peter J Lockhart

With the increasing use of massively parallel sequencing approaches in evolutionary biology, the need for fast and accurate methods suitable to investigate genetic structure and evolutionary history are more important than ever. We propose new distance measures for estimating genetic distances between individuals when allelic variation, gene dosage and recombination could compromise standard approaches. We present four distance measures based on single nucleotide polymorphisms (SNP) and evaluate them against previously published measures using coalescent-based simulations. Simulations were used to test (i) whether the measures give unbiased and accurate distance estimates, (ii) if they can accurately identify the genomic mixture of hybrid individuals and (iii) if they give precise (low variance) estimates. The results showed that the SNP-based genpofad distance we propose appears to work well in the widest circumstances. It was the most accurate method for estimating genetic distances and is also relatively good at estimating the genomic mixture of hybrid individuals. Our simulations provide benchmarks to compare the performance of different distance measures in specific situations.

Genetic Influences on Brain Gene Expression in Rats Selected for Tameness and Aggression

Genetic Influences on Brain Gene Expression in Rats Selected for Tameness and Aggression

Henrike O. Heyne, Susann Lautenschläger, Ronald Nelson, François Besnier, Maxime Rotival, Alexander Cagan, Rimma Kozhemyakina, Irina Z. Plyusnina, Lyudmila Trut, Örjan Carlborg, Enrico Petretto, Leonid Kruglyak, Svante Pääbo, Torsten Schöneberg, Frank W. Albert
(Submitted on 14 Apr 2014)

Inter-individual differences in many behaviors are partly due to genetic differences, but the identification of the genes and variants that influence behavior remains challenging. Here, we studied an F2 intercross of two outbred lines of rats selected for tame and aggressive behavior towards humans for more than 64 generations. By using a mapping approach that is able to identify genetic loci segregating within the lines, we identified four times more loci influencing tameness and aggression than by an approach that assumes fixation of causative alleles, suggesting that many causative loci were not driven to fixation by the selection. We used RNA sequencing in 150 F2 animals to identify hundreds of loci that influence brain gene expression. Several of these loci colocalize with tameness loci and may reflect the same genetic variants. Through analyses of correlations between allele effects on behavior and gene expression, differential expression between the tame and aggressive rat selection lines, and correlations between gene expression and tameness in F2 animals, we identify the genes Gltscr2, Lgi4, Zfp40 and Slc17a7 as candidate contributors to the strikingly different behavior of the tame and aggressive animals.

Comparing Evolutionary Rates Using An Exact Test for 2×2 Tables with Continuous Cell Entries

Comparing Evolutionary Rates Using An Exact Test for 2×2 Tables with Continuous Cell Entries

A. Morgan Thompson, M. Cyrus Maher, Lawrence H. Uricchio, Zachary A. Szpiech, Ryan D. Hernandez
(Submitted on 11 Apr 2014)

Assessing the statistical significance of an observed 2×2 contingency table can easily be accomplished using Fisher’s exact test (FET). However, if the cell entries are continuous or represent values inferred from a continuous parametric model, then FET cannot be applied. Such tables arise frequently in areas of biostatistical research including population genetics and evolutionary genomics, where cell entries are estimated by computational methods and result in cell entries drawn from the non-negative real line R+. Simply rounding cell entries to conform to the assumptions of FET is an ill-suited approach that we show creates problems related to both type-I and type-II errors. Pearson’s chi^2 test for independence, while technically applicable, is not often effective for these tables, as the test has several limiting assumptions that make application of this method inadvisable in many common instances (particularly with small cell entries). Here we develop a novel method for tables with continuous entries, which we term continuous Fisher’s Exact Test (cFET). Through simulations, we show that cFET has a close-to-uniform distribution of p-values under the null hypothesis of independence, and more power when applied to tables where the null hypothesis is false (compared to FET applied to rounded cell entries). We apply cFET to an example from comparative genomics to confirm an overall increased evolutionary rate among primates compared to rodents, and identify several genes that show particularly elevated evolutionary rates in primates. Some of these genes exhibit signatures of continued positive selection along the human lineage since our divergence with chimpanzee 5-7 million years ago, as well as ongoing selection in modern humans.

Selection signatures in worldwide Sheep populations

Selection signatures in worldwide Sheep populations

Maria-Ines Fariello, Bertrand Servin, Gwenola Tosser-Klopp, Rachelle Rupp, Carole Moreno, International Sheep Genomics Consortium n.a., Magali San Cristobal, simon boitard

The diversity of populations in domestic species offers great opportunities to study genome response to selection. The recently published Sheep HapMap dataset is a great example of characterization of the world wide genetic diversity in sheep. In this study, we re-analyzed the Sheep HapMap dataset to identify selection signatures in worldwide sheep populations. Compared to previous analyses, we made use of statistical methods that (i) take account of the hierarchical structure of sheep populations, (ii) make use of linkage disequilibrium information and (iii) focus specifically on either recent or older selection signatures. We show that this allows pinpointing several new selection signatures in the sheep genome and distinguishing those related to modern breeding objectives and to earlier post-domestication constraints. The newly identified regions, together with the ones previously identified, reveal the extensive genome response to selection on morphology, color and adaptation to new environments.

Natural CMT2 variation is associated with genome-wide methylation changes and temperature adaptation

Natural CMT2 variation is associated with genome-wide methylation changes and temperature adaptation

Xia Shen, Jennifer De Jonge, Simon Forsberg, Mats Pettersson, Zheya Sheng, Lars Hennig, Örjan Carlborg

As Arabidopsis thaliana has colonized a wide range of habitats across the world it is an attractive model for studying the genetic mechanisms underlying environmental adaptation. Here, we used public data from two collections of A. thaliana accessions to associate genetic variability at individual loci with differences in climates at the sampling sites. We use a novel method to screen the genome for plastic alleles that tolerate a broader climate range than the major allele. This approach reduces confounding with population structure and increases power compared to standard genome-wide association methods. Sixteen novel loci were found, including an association between Chromomethylase 2 (CMT2) and variability in seasonal temperatures where the plastic allele had reduced genome-wide CHH methylation. Cmt2 mutants were more tolerant to heat-stress, suggesting genetic regulation of epigenetic modifications as a likely mechanism underlying natural adaptation to variable temperatures, potentially through differential allelic plasticity to temperature- stress.

Majority rule has transition ratio 4 on Yule trees under a 2-state symmetric model

Majority rule has transition ratio 4 on Yule trees under a 2-state symmetric model

Elchanan Mossel, Mike Steel
(Submitted on 10 Apr 2014)

Inferring the ancestral state at the root of a phylogenetic tree from states observed at the leaves is a problem arising in evolutionary biology. The simplest technique — majority rule — estimates the root state by the most frequently occurring state at the leaves. Alternative methods — such as maximum parsimony – explicitly take the tree structure into account. Since either method can outperform the other on particular trees, it is useful to consider the accuracy of the methods on trees generated under some evolutionary null model, such as a Yule pure-birth model. In this short note, we answer a recently posed question concerning the performance of majority rule on Yule trees under a symmetric 2-state Markovian substitution model of character state change. We show that majority rule is accurate precisely when the ratio of the birth (speciation) rate of the Yule process to the substitution rate exceeds the value 4. By contrast, maximum parsimony has been shown to be accurate only when this ratio is at least 6. Our proof relies on a second moment calculation, coupling, and a novel application of a reflection principle.

The relationships among GC content, nucleosome occupancy, and exon size

The relationships among GC content, nucleosome occupancy, and exon size

Liya Wang, Lincoln Stein, Doreen Ware
(Submitted on 9 Apr 2014)

The average size of internal translated exons, ranging from 120 to 165 nt across metazoans, is approximately the size of the typical mononucleosome (147 nt). Genome-wide study has also shown that nucleosome occupancy is significantly higher in exons than in introns, which might indicate that the evolution of exon size is related to its nucleosome occupancy. By grouping exons by the GC contents of their flanking introns, we show that the average exon size is positively correlated with its GC content. Using the sequencing data from direct mapping of Homo sapiens nucleosomes with limited nuclease digestion, we show that the level of nucleosome occupancy is also positively correlated with the exon GC content in a similar fashion. We then demonstrated that exon size is positively correlated with their nucleosome occupancy. The strong correlation between exon size and the nucleosome occupancy suggests that chromatin organization may be related to the evolution of exon sizes.

Estimating Phylogeny from microRNA Data: A Critical Appraisal

Estimating Phylogeny from microRNA Data: A Critical Appraisal

Robert Thomson, David Plachetzki, Luke Mahler, Brian Moore

As progress toward a highly resolved tree of life continues to expose nodes that resist resolution, interest in new sources of phylogenetic information that are informative for these most difficult relationships continues to increase. One such potential source of information, the presence and absence of microRNA families, has been vigorously promoted as an ideal phylogenetic marker and has been recently deployed to resolve several long-standing phylogenetic questions. Understanding the utility of such markers for phylogenetic inference hinges on developing a better understanding for how such markers behave under suitable evolutionary models, as well as how they perform in real inference scenarios. However, as yet, no study has rigorously characterized the statistical behavior or utility of these markers. Here we examine the behavior and performance of microRNA presence/absence data under a variety of evolutionary models and reexamine datasets from several previous studies. We find that highly heterogeneous rates of microRNA gain and loss, pervasive secondary loss, and sampling error collectively render microRNA-based inference of phylogeny difficult, and fundamentally alter the conclusions for four of the five studies that we re-examine. Our results indicate that miRNA data have far less phylogenetic utility in resolving the tree of life than is currently recognized and we urge ample caution in their interpretation.

Bias and measurement error in comparative analyses: a case study with the Ornstein Uhlenbeck model

Bias and measurement error in comparative analyses: a case study with the Ornstein Uhlenbeck model

Gavin Huw Thomas, Natalie Cooper, Chris Venditti, Andrew Meade, Robert P Freckleton

Phylogenetic comparative methods are increasingly used to give new insight into variation, causes and consequences of trait variation among species. The foundation of these methods is a suite of models that attempt to capture evolutionary patterns by extending the Brownian constant variance model. However, the parameters of these models have been hypothesised to be biased and only asymptotically behave in a statistically predictable way as datasets become large. This does not seem to be widely appreciated. We show that a commonly used model in evolutionary biology (the Ornstein-Uhlenbeck model) is biased over a wide range of conditions. Many studies fitting this model use datasets that are small and prone to substantial biases. Our results suggest that simulating fitted models and comparing with empirical results is critical when fitting OU and other extensions of the Brownian model.