Martin Krzywinski / Genome Sciences Center / mkweb.bcgsc.ca Martin Krzywinski / Genome Sciences Center / mkweb.bcgsc.ca - contact me Martin Krzywinski / Genome Sciences Center / mkweb.bcgsc.ca on Twitter Martin Krzywinski / Genome Sciences Center / mkweb.bcgsc.ca - Lumondo Photography Martin Krzywinski / Genome Sciences Center / mkweb.bcgsc.ca - Pi Art Martin Krzywinski / Genome Sciences Center / mkweb.bcgsc.ca - Hilbertonians - Creatures on the Hilbert Curve
Trance opera—Spente le Stellebe dramaticmore quotes

space: exciting



UCD Computational and Molecular Biology Symposium, Dublin, Ireland. 1-2 Dec 2016.


visualization + design

Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
Cover image accompanying our article on mouse vasculature development. Biology turns astrophysical. PNAS 1 May 2012; 109 (18) (zoom, PNAS)

Creating the PNAS Cover

One of my goals in life, which I can now say has been accomplished, is to make biology look like astrophysics. Call it my love for the Torino Impact Hazard Scale.

Recently, I was given an opportunity to attend to this (admittedly vague) goal when Linda Chang from Aly Karsan's group approached me with some microscopy photos of mouse veins. I was asked to do "something" with these images for a cover submission to accompany the manuscript.

When people see my covers, sometimes they ask "How did you do that?" Ok, actually they never ask this. But being a scientist, I'm trained me to produce answers in anticipation of such questions. So, below, I show you how the image was constructed.

The image was published on the cover of PNAS (PNAS 1 May 2012; 109 (18))

Tools

Photoshop CS5, Nik Color Efex Pro 4, Alien Skin Bokeh 2 and a cup of coffee from a Rancilio Silvia.

source images

Below are a few of the images I had the option to work with. These are mouse embryonic blood vessels, with a carotid artery shown in the foreground with endothelial cells in green, vascular smooth muscle cells in red and the nuclei in blue.

Of course, as soon as I saw the images, I realized that there was very little that I needed to do to trigger the viewer's imagination. These photos were great!

Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
Mouse carotid arteries. (zoom)
Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
Mouse carotid arteries. (zoom)
Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
Mouse carotid arteries. (zoom)
Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
Mouse carotid arteries. (zoom)
Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
Mouse carotid arteries. (zoom)
Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
Mouse carotid arteries. (zoom)
Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
Mouse carotid arteries. (zoom)
Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
Mouse carotid arteries. (zoom)
Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
Mouse carotid arteries. (zoom)

memories of star trek

Immediately I thought of two episodes of Star Trek (original series): Doomsday Machine and the Immunity Syndrome, as well as of images from the Hubble Telescope.

Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
Enterprise is about to be consumed by a horror tube: a planet killer! (The Doomsday Machine)
Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
Enterprise heads into a giant amoeba. Who eats whom? I'll let you guess. (The Immunity Syndrome)
Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
Orion nebula (M42) as seen by the Hubble telescope. (zoom)

I though it would be pretty easy to make the artery images look all-outer-spacey. They already looked it.

centerpiece image

And then I saw the image below.

Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
A particularly spectacular image of a mouse carotid artery. I'm thinking 10 on the Torino scale. (zoom)

constructing the cover

background

The background was created from the two images shown here. The second image was sampled three times, at different rotations.

Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
Images used for background. (zoom)
Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
Images used for background. (zoom)
Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
Layer composition for background elements. (zoom)

The channel mixer was used to remove the green channel and leave only red and blue.

Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
Background elements for PNAS cover image. (zoom)

middle ground

The next layer was composed of what looked like ribbons of blue gas. This was created by sampling the oval shapes from the source images. Here the red channel was a great source for cloud shapes, and this was the only channel that was kept. The hue was shifted to blue and a curve adjustment was applied to increase the contrast.

Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
First set of middle ground elements, before adjustments. (zoom)
Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
First set of middle ground elements, after channel adjustments. (zoom)
Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
Second set of middle ground elements. (zoom)
Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
Layer composition for middle ground elements. (zoom)

When the foreground and middle ground elements were combined, the result was already 40 parsecs away.

Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
Background and foreground elements for PNAS cover image. (zoom)

foreground

The foreground was created from the spectacular comet-like image of a mouse artery. Very little had to be done to make this element look good. It already looked good.

I applied a little blur using Alien Skin's Bokeh 2 to narrow the apparent depth of field, masked out elements at the bottom of the image and removed some of the green channel. The entire blue channel was removed altogether (this gave the tail of the comet a mottled, flame-like appearance).

Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
Foreground element, before adjustments. (zoom)
Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
Foreground element, after channel adjustments. (zoom)
Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
Layer composition for foreground element. (zoom)

post processing

Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
Initial composition of background, middle ground and foreground elements. (zoom)
Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
40% localized application of Nik's Tonal Contrast (Color Efex 4 plugin) to increase structure in red channel. (zoom)
Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
50% blend with Nik's Pro Contrast (Color Efex 4 plugin). (zoom)

And here we have the final image.

Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
Final PNAS cover. Spacey! (zoom)
VIEW ALL

news + thoughts

Visualizing Clonal Evolution in Cancer

Thu 02-06-2016

Genomic instability is one of the defining characteristics of cancer and within a tumor, which is an ever-evolving population of cells, there are many genomes. Mutations accumulate and propagate to create subpopulations and these groups of cells, called clones, may respond differently to treatment.

It is now possible to sequence individual cells within a tumor to create a profile of genomes. This profile changes with time, both in the kinds of mutation that are found and in their proportion in the overall population.

Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
Ways to present temporal and phylogenetic evolution of clones in cancer. M Krzywinski (2016) Molecular Cell 62:652-656. (read)

Clone evolution diagrams visualize these data. These diagrams can be qualitative, showing only trends, or quantitative, showing temporal and population changes to scale. In this Molecular Cell forum article I provide guidelines for drawing these diagrams, focusing with how to use color and navigational elements, such as grids, to clarify the relationships between clones.

Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
How to draw clone evolution diagrams better. M Krzywinski (2016) Molecular Cell xxx:xxx-xxx. (read)

I'd like to thank Maia Smith and Cydney Nielsen for assistance in preparing some of the figures in the paper.

Krzywinski, M. (2016) Visualizing Clonal Evolution in Cancer. Mol Cell 62:652-656.

Binning High-Resolution Data

Wed 01-06-2016

Limitations in print resolution and visual acuity impose limits on data density and detail.

Your printer can print at 1,200 or 2,400 dots per inch. At reading distance, your reader can resolve about 200–300 lines per inch. This large gap—how finely we can print and how well we can see—can create problems when we don't take visual acuity into account.

Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
Nature Methods Points of View column: Binning high-resolution data. (read)

The column provides some guidelines—particularly relevant when showing whole-genome data, where the scale of elements of interest such as genes is below the visual acuity limit—for binning data so that they are represented by elements that can be comfortably discerned.

Krzywinski, M. (2016) Points of view: Binning high-resolution data. Nature Methods 13:463.

...more about the Points of View column

Regression diagnostics

Wed 11-05-2016

Residual plots can be used to validate assumptions about the regression model.

Continuing with our series on regression, we look at how you can identify issues in your regression model.

The difference between the observed value and the model's predicted value is the residual, `r = y_i - \hat{y}_i`, a very useful quantity to identify the effects of outliers and trends in the data that might suggest your model is inadequate.

Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
Nature Methods Points of Significance column: Regression diagnostics? (read)

We also discuss normal probability plots (or Q-Q plots) and show how these can be used to check that the residuals are normally distributed, which is one of the assumptions of regression (constant variance being another).

Background reading

Altman, N. & Krzywinski, M. (2016) Points of Significance: Analyzing outliers: Influential or nuisance? Nature Methods 13:281-282.

Altman, N. & Krzywinski, M. (2015) Points of Significance: Multiple Linear Regression Nature Methods 12:1103-1104.

Altman, N. & Krzywinski, M. (2015) Points of significance: Simple Linear Regression Nature Methods 12:999-1000.

...more about the Points of Significance column

Analyzing Outliers: Influential or Nuisance?

Fri 08-04-2016

Some outliers influence the regression fit more than others.

This month our column addresses the effect that outliers have on linear regression.

You may be surprised, but not all outliers have the same influence on the fit (e.g. regression slope) or inference (e.g. confidence or prediction intervals). Outliers with large leverage—points that are far from the sample average—can have a very large effect. On the other hand, if the outlier is close to the sample average, it may not influence the regression slope at all.

Martin Krzywinski @MKrzywinski mkweb.bcgsc.ca
Nature Methods Points of Significance column: Analyzing Outliers: Influential or Nuisance? (read)

Quantities such as Cook's distance and the so-called hat matrix, which defines leverage, are useful in assessing the effect of outliers.

Background reading

Altman, N. & Krzywinski, M. (2015) Points of Significance: Multiple Linear Regression Nature Methods 12:1103-1104.

Altman, N. & Krzywinski, M. (2015) Points of significance: Simple Linear Regression Nature Methods 12:999-1000.

...more about the Points of Significance column