Skip to content

A Beginners’ Guide to Non-coding Sequence Alignment

Posted in: Genomics and Epigenetics
An image of colors to depict care for your pH meter.

There is no such thing as “junk” DNA

Until recently, vast areas of the genome had been denounced as “junk” DNA, because they do not encode proteins. However, it has become clear that these regions have a large diversity of other functions, from transcriptional and translational regulation to the protection of genes and genome integrity. The ENCODE project reported in 2012 that at least 78% of the genomic sequence (in humans) serve a specific function. Most of the functions are yet unknown, and there is strong interest in developing algorithms that help uncover the logical patterns within non-coding sequences. In this article, we’ll discuss a few different software options that you can use to identify conserved non-coding elements.

Non-coding sequence alignment using MULAN

When aligning protein or mRNA sequences the software usually matches sequences by conservation, since these sequences are assumed to share common origin. However a characteristic of non-coding DNA is that functional elements can rearrange (change position, break up, invert) without losing their functionality, which makes them impossible to align with the same software. The free online software MULAN (MUltiple sequence Local AligNment and visualization tool) uses genes and surrounding regions to look for conservation in the non-coding DNA. The user has to provide sequence data from several species (depending on the depth of conservation you are looking for) for the same gene; for example, a gene plus 5kb of upstream sequence. Additionally, all exons in the area have to be annotated, because naturally they will show up as highly conserved areas. The alignment is performed pairwise, comparing each species with the species that was selected as a reference (see Figure 1: Screenshot of MULAN Output). The software can find patches of conservation that are in a different order or backwards, as is often the case with enhancer elements. Their position on the reference sequence is highlighted and the sequence alignment can be viewed and analyzed.

A Beginners' Guide to Non-coding Sequence Alignment

Figure 1: Screenshot of MULAN output. Here, the 10kb upstream region of a gene was compared between teleost (a large and extremely diverse group of ray-finned fish) and humans, using zebrafish as a reference genome. The red areas show high conservation; the closer related the species are the more non-coding elements can be expected to be conserved. In this example, one upstream and one intronic element are highly conserved from fish to human; a regulatory function of these elements is highly likely. (Tetraodon (Tetraodontidae) and fugu belong to the pufferfish genus; medaka (Oryzias latipes) belongs to the ricefish genus.)

Uncovering synteny using Genomicus

Another way to approach conservation is taking synteny into account. Genes are said to be in synteny if the same genes occur in close proximity to one another across several species. A common feature of syntenic loci is that they also share regulatory elements (the most famous example of this is the hox gene clusters). An online genome browser that searches for syntenic genes is Genomicus. The user selects gene and species and the software calculates a phylogenetic tree based on this gene and shows the surrounding genes if they are in synteny. By ticking “CNE” (conserved non-coding elements) in the view menu, the software will also show areas of non-coding conservation between the syntenic genes.

Exploring synteny and conservation with UCNE base

Finally, conservation and synteny information can also be found conveniently presented in a browse-able database at UCNE base, “a database of ultraconserved non-coding elements and genomic regulatory blocks”. It includes many model species and a lot of data, so knowing exactly what you are looking for is mandatory in this case.

Any questions? Let us know in the comments section!

Check out our related article on identifying conserved elements in genes.

References:
Dimitrieva, S., & Bucher, P. (2012). UCNEbase–a database of ultraconserved non-coding elements and genomic regulatory blocks. Nucleic acids research, 1–9. doi:10.1093/nar/gks1092

Kikuta, H., Laplante, M., Navratilova, P., Komisarczuk, A. Z., Engström, P. G., Fredman, D., Akalin, A., et al. (2007). Genomic regulatory blocks encompass multiple neighboring genes and maintain conserved synteny in vertebrates. Genome research, 17(5), 545–55. doi:10.1101/gr.6086307

Share this to your network:
Image Credit: fras1977

2 Comments

  1. James on May 3, 2016 at 3:49 pm

    I am a beginner in the field of reading genes and dna. I am looking for a beginner course. Do you know of where I should begin?

    • Dr Amanda Welch on May 3, 2016 at 4:28 pm

      What do you want or need to know? You can hop around on Bitesize Bio to learn quite a bit. If you’re looking for a course there might be something on Coursera or other MOOC site (I haven’t checked this, though).

      Anyhow, if you tell me what you’re looking for, I’m happy to help find some resources for you!

Leave a Comment

You must be logged in to post a comment.

This site uses Akismet to reduce spam. Learn how your comment data is processed.

Scroll To Top