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.
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
The development of CRISPR/Cas9 technology has made it relatively straightforward to selectively edit genomes and has revolutionized the way in which we approach biological questions. CRISPR stands for Clustered Regularly Interspaced Short Palindromic Repeats and in simple terms, this technique allows you to direct a nuclease to cut at a specific site in the genome of interest. The…
Maybe you want to examine the entire transcriptome or maybe you want to investigate changes in expression from your favorite gene. You could do whole transcriptome sequencing or mRNA-seq. But which one is right for your project? From budget considerations to sample collection, let’s briefly look at both to see which might be best for your…
Size selection is a critical step in NGS pipelines, but may be most challenging for studies of small RNAs. The concept behind size selection is simple: separate a sheared DNA or cDNA sample by fragment size, and then use the resulting sizes to remove unwanted fragments. This is a tried-and-true way to get rid of…
I’ve never run a sequencing gel in my life, but people around me did, and they spent a lot of time on getting it just right. Although the principle described by Sanger in 1975 sounds straightforward (1), sequencing gels are very long and very thin – less than a millimeter thick! They were easy to…
As we discussed previously, the gaps in our understanding of the human genome make variant classification an extremely difficult job. However, with each passing day our knowledge increases, and the tools to help us become increasingly more efficient. Let’s pick up where we left off in our first article about variants. After checking Ensemble to…
It took scientists a little while to warm up to long-read sequencing, but now you couldn’t pry most of them away from their sequencers with a crowbar. Long reads — we’re talking 10,000 bases and more — provide a level of contiguity and completeness in genome assemblies that simply isn’t possible with short-read sequencers. They…
10 Things Every Molecular Biologist Should Know
The eBook with top tips from our Researcher community.
