To analyze clustering in the occurrence of ''extreme expression'' status we coded extreme (1) vs non-extreme (0) expression as a binary variable and carried out retrospective spatial, temporal and spatio-temporal scans [29] with the software SaTScan TM v9.1.0, specifying a Bernoulli model. Following identification of a most significant cluster, subsequent clusters were identified by adjusting for the presence of more significant clusters.