Session 0
Session 0
VISUALIZATION GUIDELINES
MARTIN KRZYWINSKI
Genome Sciences Center
BC Cancer Agency
Vancouver, Canada
/BGA2017/Course_material/June08
including
—slides
http://www.circos.ca/documentation/course
what is circos?
are there any useful general rules for data visualization? (yes)
practical sessions
comparing genomes
to solve problems—
Thomson, N.R., et al., Comparative genome analysis of Salmonella Thomson, N.R., et al., Comparative genome analysis of Salmonella
Enteritidis PT4 and Salmonella Gallinarum 287/91 provides insights into Enteritidis PT4 and Salmonella Gallinarum 287/91 provides insights into
evolutionary and host adaptation pathways. Genome Res, 2008. 18(10): p. evolutionary and host adaptation pathways. Genome Res, 2008. 18(10): p.
1624-37. 1624-37.
Nielsen C, Wong B (2012) Representing genomic structural variation. Nat Methods 9:631.
Hillmer AM, Yao F, Inaki K et al. 2011 Comprehensive long-span paired-end-tag mapping reveals characteristic patterns of structural variations in
epithelial cancer genomes. Genome research 21:665-675.
The town of Caceres, Spain, a UNESCO World Heritage Site, used Circos to illustrate the relationships between businesses in their urban
planning strategy.
http://sperlingmosaics.com/
http://www.bl.uk/whatson/exhibitions/beautiful-science/
Circos addresses the need Driven entirely by plain-text Large number of data
to visualize differences in configuration files. tracks, which can be
disease genomes and stacked and layered.
assess variation in genomic Data agnostic.
content across many Format of everything in the
samples. Simple format for data input. figure can be dynamically
adjusted based on rules that
Dynamic rules provide a Highly automatable. react to data values.
way to adjust the format of
figure elements based on Fits naturally into any data Utility tools assist with
data values. pipeline. manipulating data files (e.g.
binning links and ordering
SVG output is designed for Extended longevity: ideograms to optimize
publication-quality performs only visualization, layout).
visualizations. not analysis.
(A) histogram (B) ideograms (C) histogram (D) heat map (E) links (F) highlights (G) grid (H) ticks. Format of data in tracks A, C, D, E is adjusted by
rules based on data values.
(A) Synteny between dog and human genomes. Each image represents the comparison of a single dog chromosome (bottom half of circle) with the
entire human genome. Krzywinski, M. et al. (2009). "Circos: an information aesthetic for comparative genomics." Genome Res 19(9): 1639-1645. (B)
Each image contains the entire dog and human genomes (bottom and top half of circle, respectively). Links shown are based on the same data as in
(A), but limited to a single chromosome (dog or human) for each image in the panel. http://mkweb.bcgsc.ca/circos/presentations/articles/amsci_cover
Ideogram order, scale and axis breaks help arrange ideograms to suit the data. (A) four chromosomes 5, 10, 15 and 20Mb in size. (B) ideograms can
be reordered. (C) ideogram scale can be reversed. (D) global scale of ideograms can be changed – here, chr2 2.5x and chr4 0.5x
Cropped ideogram regions are treated as individual ideograms. (E) axis breaks are used to remove the following ideogram regions chr1:2.5-3.5Mb,
chr1:4.5-5, chr2:0-2.5Mb, chr2:5.5-6.5Mb, chr2:9.5-10Mb, chr3:5.5-9.5Mb. (F) order of ideogram regions can be changed
The most frequent complex rearrangements involving MLL and (A) AFF1/AF4. Localization of chromosomal breakpoints and UPN of individual patients
are indicated. Colored lines indicate in-frame fusions (green), out-of-frame fusions (red), no partner gene present at the recombination site (blue).
Meyer, C., E. Kowarz, et al. (2009). "New insights to the MLL recombinome of acute leukemias." Leukemia 23(8): 1490-1499. Figure by M Krzywinski.
(A) each data track confined to an annulus bounded by radii r0 and r1. (B) any number of tracks can be placed on the figure. (C) tracks can be placed
at any radial position, including inside/outside ideogram circle and inside/outside ticks. (D) tracks can be made to overlap and can be drawn in any
order.
By defining three histogram tracks within the same radial region, and drawing the data in a specific order, a compound track can be created.
Various types of data tracks can be stacked. Five instances of a compound track each represent copy number information from a different sample.
Two histograms, a line plot and a scatter plot are used to form a compound track. Using links and highlights, attention is drawn to the progression of
scale increase within chr17:53-63Mb. This region is magnified at 5x and smaller subregions are further magnified to 40x.
Krzywinski, M., J. Schein, et al. (2009). "Circos: an information aesthetic for comparative genomics." Genome Res 19(9): 1639-1645.
Data sets which do not sample the genome uniformly (A) can be effectively shown by using a connector track (B) to show the remapping onto an index
scale (C). Shown in the figure are methylation values (A) for 7 tissues are summarized using stacked histograms (C), whose bins represent statistics for
remapped methylation probe positions. Zimmer, C. (2008). Now: The Rest of the Genome. New York Times. Figure by M Krzywinski.
The same data set is shown in all panels. (A) each link represents one of a subset of 2,500 segmental duplications within the human genome. (B) rules
are used to change link color and thickness. (C) rules are used to show only links to chrY. (D) in addition to rules in (C), other rules add a second layer
of links from chr8.
(E,F) adjacent links are grouped into thicker links (bundles) to reduce the complexity of the figure.
Regions of similarity between human and dog genomes. (A) human genome. (B) human ideograms. (C) dog genome. (D) dog ideograms, coded by
most similar human chromosome. (E,F) link bundles connect similar regions. (F1) rules are used to color bundles by size. (F2) bundles twist when
similarity involves opposite strands. American Scientist, Sept-Oct 2007. Cover figure by M Krzywinski.
Link color, thickness and geometry can be dynamically adjusted based on position and size of the link. Here, blue links are made to point outward and
show intrachromosomal connections. Red links emphasize connections from a short region on chr2.
Resequencing with naïve and log pool designs. Prabhu, S. and I. Pe'er, The size and outline of each scatter plot glyph is influenced by the data
Overlapping pools for high-throughput targeted resequencing. Genome value. The data value itself can be altered, as see in the two outermost
Res, 2009. 19(7): p. 1254-61. collapsed scatter plots, where the value for each point has been set to 0 to
display the glyphs at the same radius.
#
species.26.txt
gh
0
235
id=gh3
gh
235
454
id=gh4
gh
454
534
id=gh7
...
<rules>
<rule>
#
data
point
must
match
the
condition
#
for
the
rule
to
apply
condition
=
var(id)
eq
"gh1"
fill_color
=
spectral-‐10-‐div-‐1
</rule>
<rule>
condition
=
var(id)
eq
"gh2"
fill_color
=
spectral-‐10-‐div-‐2
</rule>
...
.
Each track is associated with several internal counters. The value of the By referencing the template multiple times, new tracks can be created
counters are different for each track and can be used to drive track automatically, without having change the template.
generation from a single template.
Central configuration file defines data track information and imports other configuration files that store parameters that change less frequently. Each
data file can be used for multiple tracks. PNG image output is used for immediate viewing, web-based reporting or presentation. SVG output is ideal
for high-res publication and post-processing individual elements.
circos does not perform any analysis, several tools for this are included in tools/
beautiful visualizations—yes
ugly visualizations—yes
K3
K3 Cancer CANCER
Krzywinski, M., Elements of visual style. Nat Methods, 2013. 10(5): p. 371.
top-down bottom-up
conciseness color
clarity S AT I S F I E S typeface
FORWARD REVERSE
RUN RUN
REVERSE FORWARD
SLOWING
RUN RUN
DORSAL VENTRAL
TURN TURN
VENTRAL TURN
FORWARD REVERSE
RUN RUN
REVERSE FORWARD
SLOWING
RUN RUN
DORSAL VENTRAL
TURN TURN
VENTRAL TURN
AatII NcoI
HA CB1**
genomic integration
+ removal of Neor
FLP
.
421 QPLDNSMGDSDCLHKHANN Wt
QPLDNAMGDADCLHKHANN Mutant
AatII NcoI
HA CB1**
genomic integration
+ removal of Neor
FLP
.
J. Neurosci., May 13, 2015 • 35(19):7349 –7364
ExF region C-terminus
420 430 500
RESIDUE VARIATION 3 3 2 0 3 6 4 3 3 4 2 0 0 4 2 2 0 3 4 7 7 8 6 2 0 2 0 3 2 0 3 0 3
splgrrdclvklecfrflppedt gdglddeiav
0
tasmanian devil 0
1 s
1 t
2 a n
2 s t
6 ssa a g r
7 i s y nsm–
7 i k y d–sg
8 i s y p gaa–
7 n yh s–sg
8 i n yh dssg
9 s r smg– e v l
10 mdn i h h rd kr
11 dqc l gppg p c g
12 dqc l gppg p c g i
12 dqc l r gppg p c g
12 nh dqc l gppt p c g
19 av ––rsg r str ep p c ge m v
20 20
% cells
10 10
0 0
Iso Iso
Obs Obs
Sac Sac
Ferrari-Toniolo, S., et al., Posterior Parietal Cortex Encoding of Dynamic Hand Force Underlying Hand-Object Interaction. Journal of
Neuroscience, 2015. 35(31): p. 10899-10910.
1000 1000
500 500
0 0
CitExp CitExp
CitSta CitSta
GluExp GluExp
GluSta GluSta
UpSet encoding
J. Neurosci., August 19, 2015 • 35(33):11729 –11742
17.4 28.1 neuromuscular disease
16.9 behaviour
14.2 23.0 neurological signs
19.0 neuronal cell death
6.4 12.4 apoptosis of neurons
8.9 learning
8.8 migration of neurons
8.1 neurotransmission
7.1 progressive motor neuropathy
6.2 differentiation of neurons
4.7 memory
3.6 development of hippocampus
1.8 demylinating peripheral neuropathy
1.5 cell movement of cortical neurons
1.2 development of hippocampal neurons
0.6 stress response of neurons
.
typical network analysis pipeline
FOXP2
CDH8
CDH5
typical network analysis pipeline
pipeline
identify
annotation with clustering candidate
seed list genes genes
FOXP2
FOXP2
CDH8
CDH8
BIOLOGICAL NETWORK CDH5
CDH5
typical network analysis pipeline
pipeline
identify
annotation with clustering candidate
seed list genes genes
FOXP2
FOXP2
CDH8
CDH8
BIOLOGICAL NETWORK CDH5
CDH5
identify
annotation with clustering candidate
seed list genes genes
FOXP2
CDH8
BIOLOGICAL NETWORK CDH5
Figure 1 (A) Damaged cell plasma membrane (Blue) and protein generated in Endoplasmic
EXPERIMENTAL PROCEDURE Figure 5 Raman spectra and (insert) white
light image of T cell. Each peak correspond
with unique molecular composition.
reticulum (ER) (yellow) (B) repaired plasma membrane with patch (yellow), which is synthesized
from these protein.
CONCLUSION
TECHNICAL CHALLENGES
Multi-modality system combined two photon laser wound-
Imaging speed — high enough to monitor membrane
ing technology with LTRS, RCM imaging and TPF imaging
repair process in 30 seconds.
1 3 5 for repairing monitoring has been demonstrated.
Resolution — determines the image quality and how Culture T lymphoma Trap single cell with Observe membrane
cell line LTRS repair process
accurate the laser would induce the localized wounds.
Combined multi-modality system has potential to analyze
Imaging cell stability — laser tweezer technique stabilize 2 4
Membrane staining Induce membrane injury cellular dynamics.
the cell avoiding any movement. with two photon laser
wounding technique
Declaration of conflict of interest: The authors declare no conflict of interest Contact Giselle (Yunxian Tian) [email protected] for further inquiry
Real time imaging of live cell membrane using laser
trapping, reflectance confocal microscopy, and multiphoton
fluorescence microscopy
Yunxian Tian Shangyuan Feng Yimei Huang Jianhua Zhao Eddie Shen Wenbo Wang Caigan Dud Haishan Zenga
Imaging Unit-Integrative Oncology Department, British Columbia Cancer Research Centre
Photomedicine Institute – Department of Dermatology and Skin Science, University of British Columbia
Department of Biomedical Engineering, University of British Columbia dDepartment of Urologic Sciences, University of British Columbia
Figure 1 (A) Damaged cell plasma membrane (Blue) and protein generated in Endoplasmic
EXPERIMENTAL PROCEDURE Figure 5 Raman spectra and (insert) white
light image of T cell. Each peak correspond
with unique molecular composition.
reticulum (ER) (yellow) (B) repaired plasma membrane with patch (yellow), which is synthesized
from these protein.
CONCLUSION
TECHNICAL CHALLENGES
Multi-modality system combined two photon laser wound-
Imaging speed — high enough to monitor membrane
ing technology with LTRS, RCM imaging and TPF imaging
repair process in 30 seconds.
1 3 5 for repairing monitoring has been demonstrated.
Resolution — determines the image quality and how Culture T lymphoma Trap single cell with Observe membrane
cell line LTRS repair process
accurate the laser would induce the localized wounds.
Combined multi-modality system has potential to analyze
Imaging cell stability — laser tweezer technique stabilize 2 4
Membrane staining Induce membrane injury cellular dynamics.
the cell avoiding any movement. with two photon laser
wounding technique
Declaration of conflict of interest: The authors declare no conflict of interest Contact Giselle (Yunxian Tian) [email protected] for further inquiry
Real time imaging of live cell membrane using laser
trapping, reflectance confocal microscopy, and multiphoton
fluorescence microscopy
Yunxian Tian Shangyuan Feng Yimei Huang Jianhua Zhao Eddie Shen Wenbo Wang Caigan Dud Haishan Zenga
Imaging Unit-Integrative Oncology Department, British Columbia Cancer Research Centre
Photomedicine Institute – Department of Dermatology and Skin Science, University of British Columbia
Department of Biomedical Engineering, University of British Columbia dDepartment of Urologic Sciences, University of British Columbia
Figure 5 Raman
spectra and (insert)
white light image of
Figure 1 (A) Damaged cell plasma membrane (Blue) and protein generated in Endoplasmic
EXPERIMENTAL PROCEDURE T cell. Each peak
correspond with
reticulum (ER) (yellow) (B) repaired plasma membrane with patch (yellow), which is synthesized unique molecular
from these protein. composition.
CONCLUSION
TECHNICAL CHALLENGES
Multi-modality system combined two photon laser wound-
Imaging speed — high enough to monitor membrane
ing technology with LTRS, RCM imaging and TPF imaging
repair process in 30 seconds.
1 3 5 for repairing monitoring has been demonstrated.
Resolution — determines the image quality and how Culture T lymphoma Trap single cell with Observe membrane
cell line LTRS repair process
accurate the laser would induce the localized wounds.
Combined multi-modality system has potential to analyze
Imaging cell stability — laser tweezer technique stabilize 2 4
Membrane staining Induce membrane injury cellular dynamics.
the cell avoiding any movement. with two photon laser
wounding technique
Declaration of conflict of interest: The authors declare no conflict of interest Contact Giselle (Yunxian Tian) [email protected] for further inquiry
Real time imaging of live cell membrane using laser
trapping, reflectance confocal microscopy, and multiphoton
fluorescence microscopy
Yunxian Tian Shangyuan Feng Yimei Huang Jianhua Zhao Eddie Shen Wenbo Wang Caigan Dud Haishan Zenga
Imaging Unit-Integrative Oncology Department, British Columbia Cancer Research Centre
Photomedicine Institute – Department of Dermatology and Skin Science, University of British Columbia
Department of Biomedical Engineering, University of British Columbia dDepartment of Urologic Sciences, University of British Columbia
Laser scanning microscopy provides reflectance Raman spectra provide molecular composition in vivo
confocal (RCM) image and two photon fluorescence
(TPF) image simultaneously at 15 frames per second.
Figure 3 (A) RCM imaging; (B) TPF
imaging of nucleotide distribution; (C)
TPF imaging of cell membrane; (D)
Overlay of RCM+TPF images from cell
nucleotide and membrane.
Figure 5 Raman
spectra and (insert)
CONCLUSION
TECHNICAL CHALLENGES
Multi-modality system combined two photon laser wound-
Imaging speed — high enough to monitor membrane
1 3 5 ing technology with LTRS, RCM imaging and TPF imaging
repair process in 30 seconds. Culture T lymphoma Trap single cell with Observe membrane
cell line LTRS repair process for repairing monitoring has been demonstrated.
Resolution — determines the image quality and how
accurate the laser would induce the localized wounds. 2 4
Membrane staining Induce membrane injury Combined multi-modality system has potential to analyze
Imaging cell stability — laser tweezer technique stabilize with two photon laser
wounding technique cellular dynamics.
the cell avoiding any movement.
Declaration of conflict of interest: The authors declare no conflict of interest Contact Giselle (Yunxian Tian) [email protected] for further inquiry
AIBL Active-Physical Activity
& Alzheimer’s Disease
Bernd Merkel
Department of Radiology
Royal Melbourne Hospital
postgraduate masterclass
VLSCI 10.10.2014
Click to add title
Population growth & urbanization
2
Challenge 1: water pollution
3
Challenge 2: water sanitation
4
phylogenomics of land snails
VISUALIZATION GUIDELINES
What are the major questions that the figure should help the reader answer?
What are you trying to communicate? Does the figure communicate it clearly?
Have you selected the simplest visual representation sufficient for your purpose?
Are there extraneous or ornamental elements? What can you safely remove?