Accelerating Bayesian inference for evolutionary biology models

Genome-wide identification and expression of the peroxisome proliferator-activated receptor gene family in the Tibetan highland fish Gymnocypris przewalskii

Peroxisome proliferator-activated receptor (PPAR) plays an important role in the regulation of lipid metabolism and has been widely identified in diverse species. Gymnocypris przewalskii is a native fish of the Qinghai Tibetan Plateau that survives in a chronically cold environment. In the current study, we conducted genome-wide identification of PPAR genes, revealing the existence of seven PPARs in the G. przewalskii genome. Collinearity was observed between two copies of PPARαb and PPARγ in G. przewalskii, suggesting that the additional copy might be gained through whole genome duplication. Both phylogenetic and multiple sequence alignment analyses indicated that PPARs in G. przewalskii were conserved with teleosts. The cold treatment (10 °C and 4 °C) led to the developmental delay of G. przewalskii embryos. Continuous expression of PPARs was observed during the embryonic development of G. przewalskii under normal and cold conditions, with significantly different transcriptional patterns. These results indicated that PPARs participated in the embryonic development of G. przewalskii, and were involved in the cold response during development. The current study proposed a potential role of PPARs in the cold response in the embryonic development of G. przewalskii, which shed light on understanding cold adaptation in Tibetan highland fish.

Adaptive Tree Proposals for Bayesian Phylogenetic Inference

Bayesian inference of phylogeny with MCMC plays a key role in the study of evolution. Yet, this method still suffers from a practical challenge identified more than two decades ago: designing tree topology proposals that efficiently sample tree spaces. In this article, I introduce the concept of adaptive tree proposals for unrooted topologies, that is tree proposals adapting to the posterior distribution as it is estimated. I use this concept to elaborate two adaptive variants of existing proposals and an adaptive proposal based on a novel design philosophy in which the structure of the proposal is informed by the posterior distribution of trees. I investigate the performance of these proposals by first presenting a metric that captures the performance of each proposal within a mixture of proposals. Using this metric, I compare the performance of the adaptive proposals to the performance of standard and parsimony-guided proposals on 11 empirical datasets. Using adaptive proposals led to consistent performance gains and resulted in up to 18-fold increases in mixing efficiency and 6-fold increases in convergence rate without increasing the computational cost of these analyses.

Simultaneous Bayesian inference of phylogeny and molecular coevolution

Patterns of molecular coevolution can reveal structural and functional constraints within or among organic molecules. These patterns are better understood when considering the underlying evolutionary process, which enables us to disentangle the signal of the dependent evolution of sites (coevolution) from the effects of shared ancestry of genes. Conversely, disregarding the dependent evolution of sites when studying the history of genes negatively impacts the accuracy of the inferred phylogenetic trees. Although molecular coevolution and phylogenetic history are interdependent, analyses of the two processes are conducted separately, a choice dictated by computational convenience, but at the expense of accuracy. We present a Bayesian method and associated software to infer how many and which sites of an alignment evolve according to an independent or a pairwise dependent evolutionary process, and to simultaneously estimate the phylogenetic relationships among sequences. We validate our method on synthetic datasets and challenge our predictions of coevolution on the 16S rRNA molecule by comparing them with its known molecular structure. Finally, we assess the accuracy of phylogenetic trees inferred under the assumption of independence among sites using synthetic datasets, the 16S rRNA molecule and 10 additional alignments of protein-coding genes of eukaryotes. Our results demonstrate that inferring phylogenetic trees while accounting for dependent site evolution significantly impacts the estimates of the phylogeny and the evolutionary process.

Characterization of tree shrew (Tupaia belangeri) interleukin-6 and its expression pattern in response to exogenous challenge

Tree shrews, one of the closest relatives of primates, have attracted increasing attention as a model of human diseases, particularly for viral infections. As the first line of defense against microbial pathogens, the innate immune system is crucial in tree shrews. Interleukin-6 (IL-6) is important in the pathophysiology of infection, inflammation and cancer, where it promotes disease development or sustains immune reactions. The present study aimed to obtain further insight into the tree shrew IL-6 (tsIL-6) system, and the function of tsIL-6 in the antiviral and antibacterial response. In the present study, the mRNA and genomic sequence of the tsIL-6 gene were characterized, and the tissue distribution and expression profile of this gene were analyzed in response to lipopolysaccharide (LPS) and polyinosinic:polycytidylic acid (poly I:C) treatment. The full-length tsIL-6 mRNA consisted of 1,152 bp with an open reading frame of 627 bp encoding 208 amino acids, a 5′-untranslated region (UTR) of 62 bp, and a 3′-UTR of 436 bp. The genome sequence of the tsIL-6 gene was 5,265 bp in length, comprising of five exons and four introns. The predicted tsIL-6 protein contained a 25-amino-acid-long signal peptide and a conserved IL-6 domain. Phylogenetic analysis based on the coding sequences revealed that tsIL-6 was closely related to IL-6 in humans. Residues crucial for receptor binding were completely conserved in the tree shrew protein. Reverse transcription-polymerase chain reaction analysis revealed that tsIL-6 mRNA was expressed in all examined tissues of healthy tree shrews, with high levels in the muscle and spleen. Following poly I:C challenge, the expression levels of tsIL-6 were upregulated in four tissues associated with immune system, the liver, spleen, kidney and intestine. Taken together, the molecular and bioinformatics analyses based on the IL-6 sequence revealed that the tree shrew has a close phylogenetic association with humans. These results provide insight for future investigations on the structure and function of tsIL-6.