Dataset Information

HALPER facilitates the identification of regulatory element orthologs across species.

ABSTRACT:

Summary

Diverse traits have evolved through cis-regulatory changes in genome sequence that influence the magnitude, timing and cell type-specificity of gene expression. Advances in high-throughput sequencing and regulatory genomics have led to the identification of regulatory elements in individual species, but these genomic regions remain difficult to align across taxonomic orders due to their lack of sequence conservation relative to protein coding genes. The groundwork for tracing the evolution of regulatory elements is provided by the recent assembly of hundreds of genomes, the generation of reference-free Cactus multiple sequence alignments of these genomes, and the development of the halLiftover tool for mapping regions across these alignments. We present halLiftover Post-processing for the Evolution of Regulatory Elements (HALPER), a tool for constructing contiguous regulatory element orthologs from the outputs of halLiftover. We anticipate that this tool will enable users to efficiently identify orthologs of regulatory elements across hundreds of species, providing novel insights into the evolution of traits that have evolved through gene expression.

Availability and implementation

HALPER is implemented in python and available on github: https://github.com/pfenninglab/halLiftover-postprocessing.

Supplementary information

Supplementary data are available at Bioinformatics online.

SUBMITTER: Zhang X

PROVIDER: S-EPMC7520040 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Zebrafish (Danio rerio) is an excellent model organism for the study of vertebrate development including skeletogenesis. Studies of mammalian cartilage formation were greatly advanced through the use of a cartilage specific regulatory element of the Collagen type II alpha 1 (Col2a1) gene. In an effort to isolate such an element in zebrafish, we compared the expression of two col2a1 homologues and found that expression of col2a1b, a previously uncharacterized zebrafish homologue, only partially overlaps with col2a1a. We focused our analysis on col2a1a, as it is expressed in both the stacked chondrocytes and the perichondrium. By comparing the genomic sequence surrounding the predicted transcriptional start site of col2a1a among several species of teleosts we identified a small highly conserved sequence (R2) located 1.7 kb upstream of the presumptive transcriptional initiation site. Interestingly, neither the sequence nor location of this element is conserved between teleost and mammalian Col2a1. We generated transient and stable transgenic lines with just the R2 element or the entire 1.7 kb fragment 5' of the transcriptional initiation site. The identified regulatory elements enable the tracking of cellular development in various tissues by driving robust reporter expression in craniofacial cartilage, ear, notochord, floor plate, hypochord and fins in a pattern similar to the expression of endogenous col2a1a. Using a reporter gene driven by the R2 regulatory element, we analyzed the morphogenesis of the notochord sheath cells as they withdraw from the stack of initially uniform cells and encase the inflating vacuolated notochord cells. Finally, we show that like endogenous col2a1a, craniofacial expression of these reporter constructs depends on Sox9a transcription factor activity. At the same time, notochord expression is maintained after Sox9a knockdown, suggesting that other factors can activate expression through the identified regulatory element in this tissue.

Project description:Vertebrate genomes share numerous conserved non-coding elements, many of which function as enhancer elements and are hypothesised to be under evolutionary constraint due to a need to be bound by combinations of sequence-specific transcription factors. In contrast, few such conserved elements can be detected between vertebrates and their closest invertebrate relatives. Despite this lack of sequence identity, cross-species transgenesis has identified some cases where non-coding DNA from invertebrates drives reporter gene expression in transgenic vertebrates in patterns reminiscent of the expression of vertebrate orthologues. Such instances are presumed to reflect the presence of conserved suites of binding sites in the regulatory regions of invertebrate and vertebrate orthologues, such that both regulatory elements can correctly interpret the trans-activating environment. Shuffling of binding sites has been suggested to lie behind loss of sequence conservation; however this has not been experimentally tested. Here we examine the underlying basis of enhancer activity for the Ciona intestinalis ??-crystallin gene, which drives expression in the lens of transgenic vertebrates despite the Ciona lineage predating the evolution of the lens. We construct an interactive gene regulatory network (GRN) for vertebrate lens development, allowing network interactions to be robustly catalogued and conserved network components and features to be identified. We show that a small number of binding motifs are necessary for Ciona ??-crystallin expression, and narrow down the likely factors that bind to these motifs. Several of these overlap with the conserved core of the vertebrate lens GRN, implicating these sites in cross species function. However when we test these motifs in a transgenic vertebrate they prove to be dispensable for reporter expression in the lens. These results show that current models depicting cross species enhancer function as dependent on conserved binding sites can be overly simplistic, with sound evolutionary inference requiring detailed dissection of underlying mechanisms.