Hierarchical Modelling of Haplotype Effects on a Phylogeny

Abstract

We introduce a hierarchical model to estimate haplotype effects based on phylogenetic relationships between haplotypes and their association with observed phenotypes. In a population there are many, but not all possible, distinct haplotypes and few observations per haplotype. Further, haplotype frequencies tend to vary substantially. Such data structure challenge estimation of haplotype effects. However, haplotypes often differ only due to few mutations, and leveraging similarities can improve the estimation of effects. We build on extensive literature and develop an autoregressive model of order one that models haplotype effects by leveraging phylogenetic relationships described with a directed acyclic graph. The phylogenetic relationships can be either in a form of a tree or a network, and we refer to the model as the haplotype network model. The model can be included as a component in a phenotype model to estimate associations between haplotypes and phenotypes. Our key contribution is that we obtain a sparse model, and by using hierarchical autoregression, the flow of information between similar haplotypes is estimated from the data. A simulation study shows that the hierarchical model can improve estimates of haplotype effects compared to an independent haplotype model, especially with few observations for a specific haplotype. We also compared it to a mutation model and observed comparable performance, though the haplotype model has the potential to capture background specific effects. We demonstrate the model with a study of mitochondrial haplotype effects on milk yield in cattle. We provide R code to fit the model with the INLA package.

Genetic characterization of minor carp (Labeo gonius) from Indian rivers revealed through mitochondrial ATPase 6/8 and D-loop region analysis: implications for conservation and management

“…Haplotype phylogeny may differ along the genome regions due to recombination and the process of swapping genome regions between haplotypes during the process of meiosis (Selle et al, 2021). Our aim was to find out haplotype effects by leveraging phylogenetic relationships between haplotypes.…”

The Consequences of Mitochondrial T10432C Mutation in Cika Cattle: A “Potential” Model for Leber’s Hereditary Optic Neuropathy

“…Possible direct mitogenome effects or interactions between the mitogenome and additive nuclear effects were therefore not considered but could influence the results. Possibly, more complex analyses targeting mitogenome effects [ 20 ] or SNP arrays targeting nuclear genes interacting with the mitogenome [ 21 ] would lead to more accurate conclusions. On the other hand, we accounted for possible confounding by age structure in our analysis because 44 “mutant” cows were born equally between 1997 and 2020, which adequately reflects the structure of the sample population (“normal” animals).…”

The big challenge for livestock genomics is to make sequence data pay

“…To solve this problem, one must find representations of relationships between haplotypes. Several attempts have been made using similarities between haplotypes (Hickey et al, 2013), grouping consecutive markers based on linkage disequilibrium (Cuyabano et al, 2015(Cuyabano et al, , 2014, local convolutional neural networks that represent regions of the genome as part of neural network structure (Pook et al, 2020), and by phylogenetic analysis of haplotypes (Edriss et al, 2013;Selle et al, 2021). Selle et al (2021), who developed a model for prediction based on phylogenetic relationships between non-recombining haplotypes, propose that recently developed methods for inferring genealogy along the genome in the presence of recombination and representing it as so-called tree sequences (Kelleher et al, 2019) may be useful.…”