- Lead Trainer
- Ben Moore
- Associate Trainers
- Event Date
- Virtual: Noblekinmat
- Work with the Ensembl Outreach team to get to grips with the Ensembl browser and accessing genomics data.
- Ensembl Browser Workshop: Noblekinmat Feedback Survey
Demos and exercises
The front page of Ensembl is found at ensembl.org. It contains lots of information and links to help you navigate Ensembl:
At the top left you can see the current release number and what has come out in this release. To access old releases, scroll to the bottom of the page and click on View in archive site.
Click on the links to go to the archives. Alternatively, you can jump quickly to the correct release by putting it into the URL, for example e98.ensembl.org jumps to release 98.
Click on View full list of all species.
Click on the common name of your species of interest to go to the species homepage. We’ll click on Human.
Here you can see links to example pages and to download flatfiles. To find out more about the genome assembly and genebuild, click on More information and statistics.
Here you’ll find a detailed description of how to the genome was produced and links to the original source. You will also see details of how the genes were annotated.
The current genome assembly for human is GRCh38. If you want to see the previous assembly, GRCh37, visit our dedicated site grch37.ensembl.org.
Region in detail view
Start at the Ensembl front page, ensembl.org. You can search for a region by typing it into a search box, but you have to specify the species.
To bypass the text search, you need to input your region coordinates in the correct format, which is
chromosome, colon, start coordinate, dash, end coordinate, with no spaces for example:
human 4:122868000-122946000. Type (or copy and paste) these coordinates into either search box.
Press Enter or click Go to jump directly to the Region in detail Page.
Click on the button to view page-specific help. The help pages provide text, labelled images and, in some cases, help videos to describe what you can see on the page and how to interact with it.
The Region in detail page is made up of three images, let’s look at each one in detail.
- The first image shows the chromosome:
The region we’re looking at is highlighted on the chromosome. You can jump to a different region by dragging out a box in this image. Drag out a box on the chromosome, a pop-up menu will appear.
If you wanted to move to the region, you could click on Jump to region (### bp). If you wanted to highlight it, click on Mark region (###bp). For now, we’ll close the pop-up by clicking on the X on the corner.
- The second image shows a 1Mb region around our selected region. This is always 1Mb in human, but the fixed size of this view varies between species. This view allows you to scroll back and forth along the chromosome.
You can also drag out and jump to or mark a region.
Click on the X to close the pop-up menu.
Click on the Drag/Select button to change the action of your mouse click. Now you can scroll along the chromosome by clicking and dragging within the image. As you do this you’ll see the image below grey out and two blue buttons appear. Clicking on Update this image would jump the lower image to the region central to the scrollable image. We want to go back to where we started, so we’ll click on Reset scrollable image.
- The third image is a detailed, configurable view of the region.
Here you can see various tracks, which is what we call a data type that you can plot against the genome. Some tracks, such as the transcripts, can be on the forward or reverse strand. Forward stranded features are shown above the blue contig track that runs across the middle of the image, with reverse stranded features below the contig. Other tracks, such as variants, regulatory features or conserved regions, refer to both strands of the genome, and these are shown by default at the very top or very bottom of the view.
You can use click and drag to either navigate around the region or highlight regions of interest, Click on the Drag/Select option at the top or bottom right to switch mouse action. On Drag, you can click and drag left or right to move along the genome, the page will reload when you drop the mouse button. On Select you can drag out a box to highlight or zoom in on a region of interest.
With the tool set to Select, drag out a box around an exon and choose Mark region.
The highlight will remain in place if you zoom in and out or move around the region. This allows you to keep track of regions or features of interest.
We can edit what we see on this page by clicking on the blue Configure this page menu at the left.
This will open a menu that allows you to change the image.
There are thousands of possible tracks that you can add. When you launch the view, you will see all the tracks that are currently turned on with their names on the left and an info icon on the right, which you can click on to expand the description of the track. Turn them on or off, or change the track style by clicking on the box next to the name. More details about the different track styles are in this FAQ: http://www.ensembl.org/Help/Faq?id=335.
You can find more tracks to add by either exploring the categories on the left, or using the Find a track option at the top left. Type in a word or phrase to find tracks with it in the track name or description.
Let’s add some tracks to this image. Add:
- Proteins (mammal) from UniProt – Labels
- 1000 Genomes - All - short variants (SNPs and indels) – Normal
Now click on the tick in the top left hand to save and close the menu. Alternatively, click anywhere outside of the menu. We can now see the tracks in the image. The proteins track is stranded, so you will see two tracks, one above and one below the contig, representing the proteins mapped to the forward and reverse strands respectively. The variants track is not stranded, so is found near the bottom of the image.
If the track is not giving you can information you need, you can easily change the way the tracks appear by hovering over the track name then the cog wheel to open a menu. To make it easier to compare information between tracks, such as spotting overlaps, you can move tracks around by clicking and dragging on the bar to the left of the track name.
Now that you’ve got the view how you want it, you might like to show something you’ve found to a colleague or collaborator. Click on the Share this page button to generate a link. Email the link to someone else, so that they can see the same view as you, including all the tracks you’ve added. These links contain the Ensembl release number, so if a new release or even assembly comes out, your link will just take you to the archive site for the release it was made on.
To return this to the default view, go to Configure this page and select Reset configuration at the bottom of the menu.
Genes and transcripts
You can find out lots of information about Ensembl genes and transcripts using the browser. If you’re already looking at a region view, you can click on any transcript and a pop-up menu will appear, allowing you to jump directly to that gene or transcript.
Alternatively, you can find a gene by searching for it. You can search for gene names or identifiers, and also phenotypes or functions that might be associated with the genes.
We’re going to look at the human UQCRQ gene. From ensembl.org, type UQCRQ into the search bar and click the Go button. You will get a list of hits with the human gene at the top.
Where you search for something without specifying the species, or where the ID is not restricted to a single species, the most popular species will appear first, in this case, human, mouse and zebrafish appear first. You can restrict your query to species or features of interest using the options on the left.
The gene tab
Click on the gene name or Ensembl ID. The Gene tab should open:
This page summarises the gene, including its location, name and equivalents in other databases. At the bottom of the page, a graphic shows a region view with the transcripts. We can see exons shown as blocks with introns as lines linking them together. Coding exons are filled, whereas non-coding exons are empty. We can also see the overlapping and neighbouring genes and other genomic features.
There are different tabs for different types of features, such as genes, transcripts or variants. These appear side-by-side across the blue bar, allowing you to jump back and forth between features of interest. Each tab has its own navigation column down the left hand side of the page, listing all the things you can see for this feature.
Let’s walk through this menu for the gene tab. How can we view the genomic sequence? Click Sequence at the left of the page.
The sequence is shown in FASTA format. The FASTA header contains the genome assembly, chromosome, coordinates and strand (1 or -1) – this gene is on the positive strand.
Exons are highlighted within the genomic sequence, both exons of our gene of interest and any neighbouring or overlapping gene. By default, 600 bases are shown up and downstream of the gene. We can make changes to how this sequence appears with the blue Configure this page button found at the left. This allows us to change the flanking regions, add variants, add line numbering and more. Click on it now.
Once you have selected changes (in this example, Show variants, 1000 Genomes variants and Line numbering) click at the top right.
You can download this sequence by clicking in the Download sequence button above the sequence:
This will open a dialogue box that allows you to pick between plain FASTA sequence, or sequence in RTF, which includes all the coloured annotations and can be opened in a word processor. If you want run a sequence analysis tool, download as FASTA sequence, whereas if you want to analyse the sequence visually, RTF is best for this. This button is available for all sequence views.
To find out what the protein does, have a look at GO terms from the Gene Ontology consortium. There are three pages of GO terms, representing the three divisions in GO: Biological process (what the protein does), Cellular component (where the protein is) and Molecular function (how it does it). Click on GO: Biological process to see an example of the GO pages.
Here you can see the functions that have been associated with the gene. There are three-letter codes that indicate how the association was made, as well as links to the specific transcript they are linked to.
You can also see the phenotypes associated with a gene. Click on Phenotype in the left hand menu.
To view all the sequence variants in table form, click the Variant table link at the left of the gene tab.
We also have links out to other databases which have information about our genes and may focus on other topics that we don’t cover, like Gene Expression Atlas or OMIM. Go up the left-hand menu to External references:
The transcript tab
We’re now going to explore the different transcripts of UQCRQ. Click on Show transcript table at the top.
Here we can see a list of all the transcripts of UQCRQ with their identifiers, lengths, biotypes and flags to help you decide which ones to look at.
If we were to only choose one transcript to analyse, we would choose UQCRQ-203 because it is the MANE Select and Ensembl Canonical. This means it is both 100% identical to the RefSeq transcript NM_014402.5 and both Ensembl and NCBI agree that it is the most biologically important transcript.
Click on the ID, ENST00000378670.8.
You are now in the Transcript tab for UQCRQ-203. We can still see the gene tab so we can easily jump back. The left hand navigation column provides several options for the transcript UQCRQ-203 - many of these are similar to the options you see in the gene tab, but not all of them. If you can’t find the thing you’re looking for, often the solution is to switch tabs.
Click on the Exons link. This page is useful for designing RT-PCR primers because you can see the sequences of the different exons and their lengths.
You may want to change the display (for example, to show more flanking sequence, or to show full introns). In order to do so click on Configure this page and change the display options accordingly.
Now click on the cDNA link to see the spliced transcript sequence with the amino acid sequence. This page is useful for mapping between the RNA and protein sequences, particularly genetic variants.
UnTranslated Regions (UTRs) are highlighted in dark yellow, codons are highlighted in light yellow, and exon sequence is shown in black or blue letters to show exon divides. Sequence variants are represented by highlighted nucleotides and clickable IUPAC codes are above the sequence.
Next, follow the General identifiers link at the left. Just like the External References page in the gene tab, this page shows links out to other databases such as RefSeq, UniProtKB, PDBe and others, this time linked to the transcript or protein product, rather than the gene.
Exploring the human MYH9 gene
(a) Find the human MYH9 (myosin, heavy chain 9, non-muscle) gene, and go to the Gene tab.
- On which chromosome and which strand of the genome is this gene located?
- How many transcripts (splice variants) are there and how many are protein coding?
- What is the longest transcript, and how long is the protein it encodes?
- Which transcript would you take forward for further study?
(b) Click on Phenotypes at the left side of the page. Are there any diseases associated with this gene, according to MIM (Mendelian Inheritance in Man)?
(c) What are some functions of MYH9 according to the Gene Ontology consortium? Have a look at the GO pages for this gene.
(d) In the transcript table, click on the transcript ID for MYH9-201, and go to the Transcript tab.
- How many exons does it have?
- Are any of the exons completely or partially untranslated?
- Is there an associated sequence in UniProtKB/Swiss-Prot? Have a look at the General identifiers for this transcript.
(e) Are there microarray (oligo) probes that can be used to monitor ENST00000216181 expression?
(a) Go to the Ensembl homepage (http://www.ensembl.org).
Select Search: Human and type MYH9. Click Go.
Click on either the Ensembl ID ENSG00000100345 or the HGNC official gene name MYH9.
- Chromosome 22 on the reverse strand.
- Ensembl has 23 transcripts annotated for this gene, of which six are protein coding.
- The longest transcript is MYH9-215 and it codes for a protein of 1,981 amino acids
- MYH9-201 is the transcript I would take forward for further study, as it is the MANE Select.
(b) Click on Phenotypes at the left to see the associated phenotypes. There is a large table of phenotypes. To see only the ones from MIM, type mim into the filter box at the top right of the table.
These are some of the phenotypes associated with MYH9 according to MIM: autosomal dominant deafness and Macrothrombocytopenia and granulocyte inclusions with or without nephritis or sensorineural hearing loss. Click on the records for more information.
(c) > The Gene Ontology project (http://www.geneontology.org/) maps terms to a protein in three classes: biological process, cellular component, and molecular function. Meiotic spindle organisation, cell morphogenesis, and cytokinesis are some of the roles associated with MYH9.
(d) Click on ENST00000216181.11
- It has 41 exons, shown in the Transcript summary.
Click on the Exons link in this side menu.
- Exon 1 is completely untranslated, and exons 2 and 41 are partially untranslated (UTR sequence is shown in orange). You can also see this in the cDNA view if you click on the cDNA link in the left side menu.
- P35579-1 from UniProt/Swiss-Prot matches the translation of the Ensembl transcript. Click on P35579-1 to go to UniProtKB, or click align for the alignment.
(e) Click on Oligo probes in the side menu.
Probesets from Affymetrix, Agilent, Codelink, Illumina, and Phalanx match to this transcript sequence. Expression analysis with any of these probesets would reveal information about the transcript. Hint: this information can sometimes be found in the ArrayExpress Atlas: www.ebi.ac.uk/arrayexpress/
Finding a gene associated with a phenotype
Phenylketonuria is a genetic disorder caused by an inability to metabolise phenylalanine in any body tissue. This results in an accumulation of phenylalanine causing seizures and intellectual disability.
(a) Search for phenylketonuria from the Ensembl homepage and narrow down your search to only genes. What gene is associated with this disorder?
(b) How many protein coding transcripts does this gene have? View all of these in the transcript comparison view.
(c) What is the MIM gene identifier for this gene?
(d) Go to the MANE Select transcript and look at its 3D structure. In the model 2pah, how many protein molecules can you see?
(a) Start at the Ensembl homepage (http://www.ensembl.org).
Type phenylketonuria into the search box then click Go. Choose Gene from the left hand menu.
The gene associated with this disorder is PAH, phenylalanine hydroxylase, ENSG00000171759.
(b) If the transcript table is hidden, click on Show transcript table to see it.
There are six protein coding transcripts.
Click on Transcript comparison in the left hand menu. Click on Select transcripts. Either select all the transcripts labelled protein coding one-by-one, or click on the drop down and select Protein coding. Close the menu.
(c) Click on External references.
The MIM gene ID is 612349.
(d) Open the transcript table and click on the ID for the MANE Select: ENST00000553106.6. Go to PDB 3D protein model in the left-hand menu.
The model 2pah is shown by default. It has two protein molecules in it. You may need to rotate the model to see this clearly.
Exploring a genomic region in human
(a) Go to the region from 31,937,000 to 32,633,000 bp on human chromosome 13. On which cytogenetic band is this region located? How many contigs make up this portion of the assembly (contigs are contiguous stretches of DNA sequence that have been assembled solely based on direct sequencing information)?
(b) Zoom in on the BRCA2 gene.
(c) Configure this page to turn on the LTR (repeat) track in this view. What tool was used to annotate the LTRs according to the track information? How many LTRs can you see within the BRCA2 gene? Do any overlap exons?
(d) Create a Share link for this display. Email it to your neighbour. Open the link they sent you and compare. If there are differences, can you work out why?
(e) Export the genomic sequence of the region you are looking at in FASTA format.
(f) Turn off all tracks you added to the Region in detail page.
(a) Go to the Ensembl homepage.
Select Search: Human and type 13:31937000-32633000 in the text box (or alternatively leave the Search drop-down list like it is and type human 13:31937000-32633000 in the text box). Click Go.
This genomic region is located on cytogenetic band q13.1. It is made up of eight contigs, indicated by the alternating light and dark blue coloured bars in the Contigs track. Note that KF455761.1 is a tiny contig that splits AL137143.8 in two.
(b) Draw with your mouse a box encompassing the BRCA2 transcripts. Click on Jump to region in the pop-up menu.
(c) Click Configure this page in the side menu (or on the cog wheel icon in the top left hand side of the bottom image).
Go into Repeats in the left-hand menu then select LTR. Click on the (i) button to find out more
Repeat Masker was used to annotate LTRs onto the genome.
Save and close the new configuration by clicking on ✓ (or anywhere outside the pop-up window).
There are seven LTRs overlapping BRCA2, none of them overlap exons.
(d) Click Share this page in the side menu. Select the link and copy. Get your neighbour’s email address and compose an email to them, paste the link in and send the message.
When you receive the link from them, open the email and click on your link. You should be able to view the page with the new configuration and data tracks they have added to in the Location tab. You might see differences where they specified a slightly different region to you, or where they have added different tracks.
Here is the Share link from the video answer: http://may2021.archive.ensembl.org/Homo_sapiens/Share/71a173bba78f0dbe03e48d3240424943?redirect=no;mobileredirect=no
(e) Click Export data in the side menu. Leave the default parameters as they are (FASTA sequence should already be selected). Click Next>. Click on Text.
Note that the sequence has a header that provides information about the genome assembly (GRCh38), the chromosome, the start and end coordinates and the strand. For example:
>13 dna:chromosome chromosome:GRCh38:13:32311910:32405865:1
(f) Click Configure this page in the side menu. Click Reset configuration. Click ✓.
We’re going to have a look for regulatory features in the region of a gene and investigate their activity in different cell types. We’ll start by searching for the gene KPNA2 and jumping to the location tab. Zoom out a little to see the gene plus some flanking regions.
The Regulatory Build is shown by default.
In this region we can see a number of regulatory features, including a red promoter with pink promoter flanks, cyan CTCF binding sites, yellow enhancers and lilac TF binding sites (don’t worry if you have zoomed out further or not as far and can see more/less). Refer to the legend at the bottom to see what the colours mean.
You can also click on the regulatory features to learn more. Click on the red promoter to get a pop-up.
Click on the stable ID, ENSR00000097453, to jump to the Regulation tab.
Black and grey lines indicate the positions of transcription factor binding motifs (TFBMs). Black lines are verified, which means there is evidence of the TF binding at this locus in at least one studied cell type, whereas grey ones have not been observed. Click on a TFBM then the matrix ID to see the matrix as a pop-up window.
Close the motif pop-up and scroll down the page. We can see that this promoter is active in one out of the 118 cell types currently in Ensembl.
We can explore more detailed data in Details by cell type – click on the button at the top.
At the moment the page is not displaying any data as we haven’t chosen any cell types. Click on Configure Cell/Tissue to add more to the view.
We can add cells by clicking on them. Find them using the search or the alphabet ribbon. Let’s add a cell type where the promoter is inactive, aorta, and one where it’s active astrocytes. Once you’ve selected the cells, they will appear in the menu on the right, where you can easily view the list by clicking on the + and deselect them.
To choose the experiments to see data on, click on the Experiments tab at the top of the menu. You can navigate this the same as the Cell/Tissue tab, except that you have to choose between Histone modifications, Open Chromatin and Transcription Factors. Let’s Select all in all categories.
When you’ve chosen your experiments and cells, you can click on the green Configure track display button.
Now we can see the active feature in astrocytes compared to the inactive feature in aorta.
Gene regulation: Human INSIG1 and BLACE
(a) Find the Location tab (Region in detail page) for the region between the genes INSIG1 and BLACE. Are there any predicted enhancers in this region?
(b) Go to the Regulation tab for the enhancer ENSR00001133586. How many cell types is this enhancer active in? Are there any cell types where its activity is repressed?
(c) Go to the Details by cell type page. Take a look at the histone modifications across this enhancer in neutro myelocyte cells, where this enhancer is active, compared to neutrophil (CB) cells, where it is poised. What differences can you observe?
(d) Go back to the Summary page. Are there any verified transcription factor binding motifs in this enhancer? In what cells?
(a) Search for human INSIG1 from the home page. Click on INSIG1 genomic coordinates 7:155297776-155310235:1 in the search results to open the Location tab.
In the region overview, drag out a box to encompass the neighbouring BLACE gene. Have a look at the Regulatory Build track.
There is one yellow enhancer in the region between these genes.
(b) Click on the enhancer to get a pop-up with its ID and other information. Click on the ID ENSR00001133586 to open the Regulation tab.
Scroll down to see the summary of regulatory activity across different cell types.
ENSR00001133586 is active in neutro myelocyte cells. It is repressed in 34 cell types.
(c) Click on Details by cell type at the top of the page or in the left-hand menu.
Choose cells by clicking on Configure Cell/Tissue blue button, then selecting neutro myelocyte in which this enhancer is active and neutrophil (CB) in which it is poised.
Add experiment tracks by clicking on Experiments tab and Select all under Histone. Click Configure track display, then View tracks to load the page.
Both cell types have H3K27me3, H3K4me1 and H3K9me3 modifications at this locus, while neutro myelocyte cells also have H3K27ac and H3K36me3 modifications. Different clusters of peaks indicate different epigenetic profiles, which might explain the difference in the enhancer activity between these two cell types.
(d) Go back to the Summary view. In the Motif features track there are two black markers indicating verified TF motifs.
Click on them to tell which motifs and which cells.
The two motifs are both verified in K562 cells and bind a number of different TFs. One binds ELF1, ELF2, ELK1, FLI1, ERG, ETS1, ETV6, FOXO1::ELK3, FOXO1::ETV1, ETV1, ETV2, ERF, ELK3, ETV3, GABPA, ETS2, ELK4, FEV, ETV5 and ETV4. The other binds ETV7, ETS1 and ELK1::SPDEF.
Regulatory features in human
(a) Search for the regulatory feature ENSR00000262400. What type of feature is this? What is its genomic location?
(b) Which cell types is this feature inactive or repressed in? Look at the Details by cell type, why do you think it has been called differently in these cell types when compared to all other cells?
(c) Why do so many cells have this feature listed as NA on the summary page?
(a) Search for ENSR00000262400 from the home page. Click on the search result to open the Regulation tab.
ENSR00000262400 is a CTCF binding site found at Chromosome 11: 1,998,001-2,001,400, which can be found at the top of the Summary page.
(b) Scroll down to see the summary of regulatory activity across different cell types.
The CTCF binding site is inactive in H1-hESC_3 and HepG2 cells. It is repressed in A673.
Click on Details by cell type at the top of the page or in the left-hand menu. Choose cells by clicking on Configure Cell/Tissue blue button, then Selecting all 148 cell types. Add the CTCF binding track by clicking on Experiments tab, Transcription factors and CTCF. Click Configure track display, then View tracks to load the page.
CTCF-binding has been observed in cells where it is active and in some cells where it is poised or inactive. It cannot be observed in cells where it is repressed.
(c) Note that many cell types have this feature represented as white block with no corresponding CTCF signal or peaks.
Cells which do not have CTCF ChIP-seq data cannot have an activity listed for this feature.
Demo: Upload small files
We have some patients that present with microcephaly and developmental delay. They all have large scale deletions on chromosome five:
We can turn them into a BED file and view them in the genome browser:
chr5 36821632 37091234 P1
chr5 36731476 36978306 P2
chr5 36908552 37108671 P3
You can add data from a Region in Detail page by clicking on the Custom tracks button at the left. Alternatively, go to a species homepage and click on Display your data in Ensembl.
A menu will appear:
The interface detects file types if you upload or attach a file. When you paste in your data, it can’t do this so we have to tell it what our file type is. It will give you an option where you can select BED.
Click Add data.
You should get to a dialogue box telling you your upload has been successful.
Click on the genomic coordinates link to go to the nearest region with data.
To have a look at the file, click on Custom tracks.
If you’ve got an Ensembl account, you can save this data to your account. Accounts are free to set up and allow you to save configurations and data, and share with groups.
Demo: Attach URLs of large files
Larger files, such as BAM files generated by NGS, need to be attached by URL. I’ve put a BAM file of human chromosome 20 RNASeq data online at: http://ftp.ebi.ac.uk/pub/databases/ensembl/training/emily_BAM/
Let’s take a look at the folder.
Here you can see a number of BAM files (
.bam) with corresponding index files (
.bam.bai). We’re interested in the files GRCh38.20.illumina.merged.1.bam and GRCh38.20.illumina.merged.1.bam.bai. These files are the BAM file and the index file respectively. When attaching a BAM file to Ensembl, there must be an index file in the same folder.
To attach the file, click on Custom tracks, then click on Add more data to add a new track.
We get to the same dialogue box as before. This time we’ll name our data Illumina reads.
Paste in the URL of the BAM file itself (http://ftp.ebi.ac.uk/pub/databases/ensembl/training/emily_BAM/GRCh38.20.illumina.merged.1.bam).
Since this is a file, the interface is able to detect the “.BAM” file extension, so automatically labels the format as BAM. Click on Add data and close the menu.
To see this data, jump to a region on chromosome 20. Let’s go to the region of the CDH22 gene. Search for the gene and click on the location.
We can zoom in to see the sequence itself. Drag out boxes in the view to zoom in, until you see a view like this. Alternatively, jump to location 20:46241000-46241030.
Demo: Track hub registry
Our regulatory data incorporates data from sources such as ENCODE, Blueprint, and Roadmap Epigenomics. To see the data directly from these sources, you can add track hubs.
You can search for track hubs to add in different ways:
- Search for track hubs in the Track Hub Registry and choose to add them to your genome browser of choice.
- Search the track hub registry using the Track Hub Registry interface in Ensembl (there is a link from the homepage).
We will now add the track hub containing data from the Blueprint project.
You can add track hubs to view in Ensembl directly via the Track Hub Registry. Go to the Track Hub Registry homepage and search for blueprint.
There are two results for the Blueprint Hub, one for adding the track hub to GRCh37 and one for adding it to GRCh38, plus one RNA-seq alignment hub.
Alternatively, you can add track hubs by searching the Track Hub Registry through Ensembl. Click the Custom tracks -> Track Hub Registry Search in any region view within Ensembl.
You can only find track hubs for the selected species and assembly denoted in the search box.
Search for blueprint.
Click Attach this hub in the search results page.
Track Hubs often contain vast amounts of data, which can slow Ensembl down, so only add them if you need them, and trash them when you are finished with them.
Go to Configure this Page to see that a new category has been added to your menu.
You can add tracks to the Region in Detail view using the matrix.
Attaching the ENCODE track hub in human
(a) Add the ENCODE Track Hub to the Region in Detail View for the genomic region surrounding the BRCA2 gene. Hint: You will need to add and view this Track Hub to the human GRCh37 genome assembly.
(b) Turn on all the available tracks relating to Histone Modification Peaks and Transcription Factor Peaks in HeLa-S3 cells.
(c) Which Transcription Factors and Histone Modifications have features in this region?
(d) Add the Tracks showing Signals for the Histone Modifications and Transcription Factors that have peaks in this region. Compare the signal intensity to the location of annotated peaks.
(e) Remove the ENCODE Track Hub from your list of custom tracks.
(a) There are two ways to add the ENCODE Track Hub to view in Ensembl.
- Search for Encode from the Track Hub Registry home page. From the search results, find the Encode Analysis Hub and select Ensembl in the View in Genome Browser list.
- Search for the BRCA2 gene in the GRCh37 human genome assembly. Switch to the Location Tab. Click on Custom tracks in the side panel. Click Track Hub Registry Search and Search for Encode. Click Attach this hub in the ENCODE Analysis Hub option.
(b) Go to Configure this page and click on ENCODE Histone Modification Peaks under ENCODE Analysis Hub. Turn on all available tracks for HeLaS3 cells by selecting HeLaS3 from the Cell Line tab and clicking Select all in the Factor tab. Go to the second step 2. Refine selection and make sure all tracks are on (blue colour). You can go to the third step to change the track display (3. Configure track display).
Click on ENCODE Transcription Factor Peaks now. Turn on all available tracks for HeLaS3 cells by selecting HeLaS3 from the Cell Line tab and clicking Select all in the Factor tab. Go to the second step 2. Refine selection and make sure all tracks are on (blue colour). You can go to the third step to change the track display (3. Configure track display) or proceed straight to View tracks.
(c) There are features for a number of different Transcription Factors and Histone modifications, mainly surrounding the BRCA2 5’ region.
Transcription Factors: CTCF, USF2, TBP, STAT1, MAX, POL2S2, POLR2A, POL2, POL2B, MXI1, INI1, E2F1, E2F4, E2F6, HAE2F1, ELK4, MYC, CMYC, CEBPB, TAF1, TFAP2A, TFAP2C, AP2alpha and AP2gamma.
Histone Modifications: H3K4me3, H3K9ac, H3K79me2, H3K4me3, H3K4me2, H3K36me3, H3K27ac
(d) Go to Configure this page and click on ENCODE Histone Modifications Signal. Turn on all available tracks for HeLaS3 cells by selecting HeLaS3 from the Cell Line tab and clicking Select all in the Factor tab. Go to the second step 2. Refine selection and make sure all tracks are on (blue colour). You can go to the third step to change the track display (3. Configure track display).
Go to Configure this page and click on ENCODE Transcription Factor Signal. Turn on all available tracks for HeLaS3 cells by selecting HeLaS3 from the Cell Line tab and clicking Select all in the Factor tab. Go to the second step 2. Refine selection and make sure all tracks are on (blue colour). You can go to the third step to change the track display (3. Configure track display) or proceed straight to View tracks.
By comparing the signal intensity and annotated peaks for each of the Histone Modifications and Transcription factors, you can see that the increased signal intensity corresponds to the regions where a peak has been annotated.
(e) Go to Custom tracks and click the Trash icon from the Actions section of the ENCODE Track Hub.
Adding Wiggle files to Ensembl Bacteria
Upload the GD_wiggle.wig file to the Gluconacetobacter diazotrophicus PA1 5 (GCA_000021325) genome in Ensembl Bacteria. View this track across the region Chromosome:2884000-2898000. What is the highest score in this region?
Go to Ensembl Bacteria and put Gluconacetobacter diazotrophicus PA1 5 into the Search for a genome box. Select Gluconacetobacter diazotrophicus PA1 5 (GCA_000021325) to go to the species homepage.
Select Display your data in Ensembl Bacteria to get to the custom track menu. Select Choose file and select the file location. The file type should be automatically selected. Click Add data.
Click on the Nearest region with data in the results page. From the region page you reach, put the coordinates Chromosome:2884000-2898000 into the Location box to jump to the region.
The highest score is 99 and it overlaps the ACI52364 transcript.
Adding BAM files to Ensembl Fungi
Attach the file Spom_all_61G9EAAXX_and_61G9UAAXX.+.sorted.bam, which can be found here, to view in Schizoaccharomyces pombe. Go to the region I:490843-490870, can you see a mismatch between a read and the reference assembly?
Go to the provided public directory and right click on the file Spom_all_61G9EAAXX_and_61G9UAAXX.+.sorted.bam to copy its URL, which is: http://ftp.ebi.ac.uk/ensemblgenomes/pub/misc_data/bam/fungi/Spom/Spom_all_61G9EAAXX_and_61G9UAAXX.%2B.sorted.bam Go to Ensembl Fungi and select Schizosaccharomyces pombe to go to the species homepage.
Select Display your data in Ensembl Fungi to get to the custom track menu. Paste the file URL into the Data box. The file type should be automatically selected. Click Add data.
Close menu and then put I:490843-490870 into the search box.
In the reads there are three red bases which do not match the reference assembly: A on the forward strand, and G and T on the reverse.