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 config for every UCSC genome, each carrying UCSC's pairwise liftOver alignments to the others. This page compares hg38 against T2T-CHM13 (hs1) at TNNT3, a locus the two lay out differently.
Opening a liftOver track
Open hg38 on genomes.jbrowse.org, find Pairwise
alignments → liftOver in the track selector and turn on hg38 to Human
(hs1) liftOver. Type TNNT3 into the location box; the hosted config ships a
name index.
The hg38 config declares only hg38. A track referencing an assembly JBrowse does
not know is resolved through the Core-handleUnrecognizedAssembly extension
point, where the site's hub plugin supplies the mate genome's config.
In a plain linear genome view the liftOver track draws one feature per chain block, laid out in rows.
Launching a synteny view
Right-click any chain block and choose Linear synteny view with ..., which names the other assembly. The dialog frames the second panel. Use CIGAR to map the current visible region to the target walks the alignment to find the interval matching what is in view; the chain through TNNT3 spans the chromosome, so unticked it frames both panels on all of chromosome 11. A reverse-strand block adds Horizontally flip inverted targets, ticked by default.
Open in new view appends the result below the linear view; Replace current view puts it in that view's place.
The panel you launched from keeps the tracks that view had on (Copy this view's tracks into its panel turns that off). The hs1 panel opens empty; the view header's track selector button lists one selector per panel.
A locus no single chain block covers takes a second route: drag-select it on the scale bar and pick Launch → Linear synteny view, whose dialog offers every assembly the session's synteny datasets align to it. See the linear synteny view guide.
Ribbon display settings
Two settings live in the synteny view's settings menu, the sliders button in its header:
- Curved lines draws each ribbon as a curve, easier to follow across a gap
- CIGAR indels → Transparent indels leaves insertions and deletions inside each block as see-through gaps
The palette button in the same header sets what ribbons are colored by. Strand paints each block by its orientation, and is what the figures below use.
The TNNT3 rearrangement
TNNT3 is the locus from Fig 5C of the T2T human variation paper. Against GRCh38 the region reads as an inversion plus a deletion that ablates LINC01150 in every individual; against T2T-CHM13 that segment is intact on the other side of TNNT3 in the opposite orientation. Colored by strand, it is the one off-color ribbon.
Trying other pairs
The same click-path works for any track under Pairwise alignments → liftOver, one per chain file UCSC publishes against the genome you are in. A close pair gives long collinear blocks, a distant one short scattered ones.
The figure below is that route on hg38 to Chimp (panTro6) liftOver, across an intron of FTO.
Its last frame switches to curves and Transparent indels (above), which turns the one gap into a hole lining up against the RepeatMasker track. The element under it is an L1HS, the youngest human LINE-1 subfamily, and the chimp panel has every other repeat in the window but not that one.
The chimp panel's track selector offers NCBI RefSeq - RefSeq All and
RepeatMasker, brought in with the panTro6 hub. The rest of that hub loads
from File → Open connection as a JBrowse 2 hub at
https://jbrowse.org/ucsc/panTro6/config.json.
To start from a gene, the site's ortholog search lists a symbol's NCBI orthologs among the hosted genomes, with a synteny view per row where an alignment exists.
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.