Alzheimer's gene co-expression · Workflow

From hippocampus expression data for 31 people to the candidate genes whose expression pattern follows PLD2 or TPP2 in every disease stage.

Each row is one probe, a single gene measurement on the microarray; each column is one person. The 31 columns run Control (9 people), Incipient (7), Moderate (8) and Severe (7).

Each group’s matrix is standardized with a single mean and standard deviation computed over all its values, rather than one per probe. t-SNE and K-Means then run on each disease stage separately, with the settings below.

Analysis stepSettings
Standardization(value − matrix mean) / matrix standard deviation
t-SNEPerplexity 40 · pca initialization · 300 iterations · seed 42
K-Means35 clusters · 10 initializations · seed 42

Each dot is one probe. For the chosen disease stage, t-SNE arranges all probes so that those with similar expression patterns across that stage's people sit close together, and K-Means cuts the map into clusters. Anchor’s cluster shows the anchor, the candidates and the other probes in the anchor's cluster at this stage. All clusters gives every cluster its own colour; the colours are labels, not an order.

Disease stage
Anchor gene
Colour by
anchor probe candidate: in the anchor's cluster in all three stages in the anchor's cluster at this stage only other probes

The candidates as a table, with their correlations to the anchor
RowProbeGene r · incipientr · moderater · severe
2243202716_atPTPN10.7110.2810.630
2482202954_atUBE2C0.7420.9340.765
2701203174_s_atARFRP10.7620.2060.680
3151203625_x_atSKP20.7560.5330.817
3623204097_s_atRBMX20.8510.7170.731
4272204746_s_atPICK10.441-0.0210.400
4547205021_s_atFOXN30.8560.4280.608
6118206593_s_atMED220.8090.7160.851
6677207152_atNTRK20.458-0.3720.217
7444207929_atGRPR0.6400.0220.803
8581209088_s_atUBN10.824-0.1220.648
10879211471_s_atRAB360.8410.2360.818
11541212156_atVPS390.766-0.1380.826
12170212786_atCLEC16A0.7400.1730.697
12376212994_atTHOC20.770-0.1860.743
12769213389_atZNF5920.7310.7600.697
13112213733_atMYO1F0.8290.5280.755
13700214322_atCAMK2G0.349-0.1730.676
13936214558_atGPR120.6980.7920.682
14739215366_atSNX130.7560.7160.474
15126215754_atSCARB20.7990.7110.631
15130215758_x_atZNF930.8850.6030.557
15984216615_s_atHTR3A0.7680.3480.740
16160216791_atAK0264490.8450.7530.790
16633217269_s_atTMPRSS150.8110.2290.618
16943217579_x_at217579_x_at0.6310.6180.576
17742218378_s_atLOC1006309230.9200.3950.881
18235218872_atTESC0.7160.3510.464
18536219173_atMYO15B0.8700.5970.878
18557219194_atSEMA4G0.8310.7280.700
18721219358_s_atADAP20.9570.7850.935
19260219897_atRNF1220.7610.5840.584
20102220739_s_atCNNM30.7180.0040.786
20435221073_s_atNOD10.9210.6830.579
20797221435_x_atHYI0.8920.3110.668
20936221575_atSCLY0.6760.5780.634
20990221629_x_atHGH10.8010.6590.585
21286221927_s_atABHD110.6680.5050.767
21411222052_atC19orf540.8410.7220.785
21655222297_x_atRPL180.9010.9670.773
2187737408_atMRC20.8140.5960.855
2189137796_atSAP250.4890.8330.615
2189937953_s_atASIC10.1080.0850.775
2198941858_atPGAP20.7790.7300.846
2218565770_atRHOT20.6150.4530.736

K-Means needs the number of clusters, k, to be set in advance. The curve scores every k that the original notebook tried with the Calinski–Harabasz score, which is higher when clusters are tight and well separated. Both the chosen k and the best-scoring k are marked.

The analysis finds 45 candidates for PLD2 and 9 for TPP2: probes that share the anchor’s cluster on the incipient, moderate and severe maps. The control group is not used to select them.

Main findings

  • Each candidate shares the anchor's cluster on the t-SNE map of every disease stage, which reflects a similar expression pattern across that stage's people.
  • The number of clusters is a choice: k = 35 sits where the score's rise slows, and the sweep keeps climbing to its highest value at k = 45.

probe (probe set)
one measurement feature on the microarray. A row of the matrix is one probe; its columns are the people it was measured in.
anchor probe
the probe the search starts from — here the ones for the genes PLD2 and TPP2. A candidate is a probe that shares the anchor's cluster in all three disease stages.
global z-score
as used here: every value in one severity matrix has that whole matrix's mean subtracted and is divided by that whole matrix's standard deviation. One number per matrix, not per row — which is why changing it would move every distance on the map.
t-SNE
a way of placing high-dimensional points on a flat map so that points which were close stay close. As used here, close means similar expression pattern across people. The axes, the compass directions and the distances between separate clumps mean nothing at all.
K-Means, and k
a way of cutting a map into k groups. As used here k is 35, chosen where the quality score starts to level off.
Calinski-Harabasz score
a number that rewards clusters for being tight and far apart. Used here to show where the chosen k sits on the curve.
the anchor’s cluster
the K-Means cluster that contains the anchor on one disease stage’s map. It groups probes by expression pattern.
candidates (the shortlist)
the probes that share the anchor's cluster in all three disease stages. A probe that qualifies in two stages out of three does not appear. The control group is not used.
Pearson correlation, r
how tightly two probes' readings rise and fall together across the people in one stage, from −1 to 1. Computed here on the raw readings, not the standardised ones, and shown per stage because a single number would hide the disagreement between stages.