Synteny on genomes.jbrowse.org
TL;DR: genomes.jbrowse.org already carries UCSC's pairwise liftOver alignments for every genome, so you can turn one on in a linear genome view and launch a two-panel linear synteny view from any chain block. JBrowse resolves the mate genome on demand, so the second assembly needs no setup.
Prerequisites
- nothing to install: this is a click-path through a hosted site, and no data, config or second assembly is prepared by hand
Where the data comes from
genomes.jbrowse.org's own hosted configs for the three genomes this page compares, each already carrying UCSC's pairwise liftOver chains to the others.
- hg38: https://jbrowse.org/ucsc/hg38/config.json
- Human (hs1, T2T-CHM13): https://jbrowse.org/ucsc/hs1/config.json
- Chimp (panTro6), for Trying other pairs: https://jbrowse.org/ucsc/panTro6/config.json
The hosted configs
genomes.jbrowse.org hosts a ready-made JBrowse config for every UCSC genome, and each one already carries UCSC's pairwise liftOver alignments to the other genomes, so any pair UCSC has a chain file for can be compared in a couple of clicks.
This tutorial compares hg38 against T2T-CHM13 (hs1) at TNNT3, a locus that GRCh38 and T2T-CHM13 lay out differently.
Opening a liftOver track
Open hg38 on genomes.jbrowse.org, then in the
track selector find Pairwise alignments → liftOver and turn on hg38 to
Human (hs1) liftOver. Type TNNT3 into the location box to navigate to the
gene; the hosted config ships a name index, so gene symbols work without any
setup.
The hg38 config declares only hg38. When a track references an assembly JBrowse
does not know, it asks its plugins to resolve the name (the
Core-handleUnrecognizedAssembly extension point) and the site's hub plugin
supplies the config, which is how every liftOver track on the site reaches its
mate genome.
The liftOver track is a synteny track, and in a plain linear genome view it draws the way an alignments track does: one feature per chain block, laid out in rows.
Launching a synteny view
Right-click any chain block and choose Launch synteny view for this position. The dialog that opens controls how the second panel is framed, and its defaults suit the block you clicked. Use CIGAR to map the current visible region to the target is the option to note: JBrowse walks the alignment to find the interval matching what is in view. The chain running through TNNT3 spans the chromosome, so unticked it frames both panels on the whole of chromosome 11. A reverse-strand block adds another checkbox, Horizontally flip inverted targets, ticked by default so the target panel reads left to right.
The dialog then offers two ways out, both building the same view. Open in new view appends it below the linear view you launched from; Replace current view puts it in that view's place, which is usually what you want here, since the synteny view opens on the locus the linear view is already showing.
The panel for the assembly you launched from opens with the tracks that view had on, and Copy this view's tracks into its panel in the dialog turns that off. The hs1 panel opens empty, and its tracks come from its own track selector, whose button in the view header lists one per panel.
A locus that no single chain block covers (a gene that several blocks tile, or one you have only navigated to) takes a second route: drag-select it on the scale bar and pick Launch → Linear synteny view, whose dialog picks the synteny dataset to read it from. It starts from the region, and offers every assembly the session's synteny datasets align to it, so use it when more than one liftOver track is on. See the linear synteny view guide.
Ribbon display settings
Three settings change how the ribbons read. Two are in the synteny view's settings menu, the sliders button in its header:
- Curved lines, ticked, draws each ribbon as a curve, so a block that lands far from where it started is easier to follow across the gap.
- CIGAR indels → Transparent indels leaves the insertions and deletions inside each block as see-through gaps, so with strand coloring on the color means one thing.
The third is the palette button further along the same header, whose menu sets what every ribbon is colored by. Strand paints each block by the orientation it landed in, and is the setting the TNNT3 figures below use.
The TNNT3 rearrangement
TNNT3 is the locus from Fig 5C of the T2T human variation paper. Called against GRCh38 the region reads as a 24 kb inversion plus a 22 kb deletion that ablates LINC01150 in every individual; against T2T-CHM13 that same 22 kb is intact, just sitting on the other side of TNNT3 in the opposite orientation. Colored by strand, that segment is the one off-color ribbon in the view.
Trying other pairs
The same click-path works for any liftOver track under Pairwise alignments → liftOver. How much of the chain survives changes with the pair: a close comparison like hs1 or panTro6 gives long collinear blocks, a distant one short scattered ones, and the CIGAR option matters more the more diverged the pair is. That liftOver category is itself the index of which pairs exist, with one track in it per chain file UCSC publishes against the genome you are in.
The figure below is that route on hg38 to Chimp (panTro6) liftOver, across an intron of FTO, and its clicks are the ones in the clip above.
Its last frame switches the ribbon to curves and to Transparent indels, both covered above, which turns the one gap in it into a hole that lines up against the RepeatMasker track over it. The element under it is an L1HS, the youngest human LINE-1 subfamily, and the chimp panel keeps every other repeat in the window but not that one.
The chimp panel's own track selector offers NCBI RefSeq - RefSeq All and
RepeatMasker, because naming panTro6 also brought in the panTro6 hub's gene,
repeat and gap tracks. The rest of that hub loads from File → Open
connection as a JBrowse 2 hub at
https://jbrowse.org/ucsc/panTro6/config.json, which is where its conservation
and expression tracks are.
To start from a gene, the site's ortholog search takes a gene symbol and lists its NCBI orthologs among the hosted genomes, with a synteny view per row wherever the two assemblies have an alignment.
See also
- Basic usage of genomes.jbrowse.org
- Synteny visualization (pairwise minimap2)
- Synteny from an ortholog table (grape, peach, cacao)
- Synteny visualization (all-vs-all minimap2)
- Linear synteny view
- Dotplot view
- The T2T human variation paper
Feedback on this tutorial is welcome: contact us.