Evolution of gene regulation in ruminants differs between evolutionary breakpoint regions and homologous synteny blocks
- Marta Farré1,
- Jaebum Kim2,
- Anastasia A. Proskuryakova3,
- Yang Zhang4,
- Anastasia I. Kulemzina3,
- Qiye Li5,
- Yang Zhou5,
- Yingqi Xiong5,
- Jennifer L. Johnson6,
- Polina Perelman3,
- Warren E. Johnson7,
- Wesley C. Warren8,
- Anna V. Kukekova6,
- Guojie Zhang5,
- Stephen J. O'Brien9,
- Oliver A. Ryder10,
- Alexander S. Graphodatsky3,
- Jian Ma4,
- Harris A. Lewin11 and
- Denis M. Larkin1,12
- 1 Royal Veterinary College;
- 2 Konkuk University;
- 3 Institute of Molecular and Cellular Biology;
- 4 Carnegie Mellon University;
- 5 BGI-Shenzen;
- 6 University of Illinois at Urbana-Champaign;
- 7 Smithsonian Conservation Biology Institute;
- 8 University of Missouri;
- 9 St. Petersburg State University;
- 10 San Diego Zoo;
- 11 University of California
Abstract
The role of chromosome rearrangements in driving evolution has been a long-standing question of evolutionary biology. Here we focused on ruminants as a model to assess how rearrangements may have contributed to the evolution of gene regulation. Using reconstructed ancestral karyotypes of Cetartiodactyls, Ruminants, Pecorans, and Bovids, we traced patterns of gross chromosome changes. We found that the lineage leading to the ruminant ancestor after the split from other cetartiodactyls, was characterized by mostly intrachromosomal changes while the lineage leading to the pecoran ancestor (including all livestock ruminants) included multiple interchromosomal changes. We observed that the liver cell putative enhancers in the ruminant evolutionary breakpoint regions are highly enriched for DNA sequences under selective constraint acting on lineage-specific transposable elements (TEs) and a set of 25 specific transcription factor (TF) binding motifs associated with recently active TEs. Coupled with gene expression data, we found that genes near ruminant breakpoint regions exhibit more divergent expression profiles among species, particularly in cattle, which is consistent with the phylogenetic origin of these breakpoint regions. Notably, this divergence was significantly greater in genes with enhancers that contain at least one of the 25 specific TF binding motifs and located near bovidae-to-cattle lineage breakpoint regions. Taken together, by combining ancestral karyotype reconstructions with analysis of cis regulatory element and gene expression evolution, our work demonstrated that lineage-specific regulatory elements co-localized with gross chromosome rearrangements may have provided valuable functional modifications that helped to shape ruminant evolution.
- Received June 25, 2018.
- Accepted February 8, 2019.
- Published by Cold Spring Harbor Laboratory Press
This article is distributed exclusively by Cold Spring Harbor Laboratory Press for the first six months after the full-issue publication date (see http://genome.cshlp.org/site/misc/terms.xhtml). After six months, it is available under a Creative Commons License (Attribution-NonCommercial 4.0 International), as described at http://creativecommons.org/licenses/by-nc/4.0/.
