This is an Open Access article distributed under the terms of the Creative Commons Attribution License (
Altered gene expression is an important feature of ischemic cerebral injury and affects proteins of many functional classes. We have used microarrays to investigate the changes in gene expression at various times after middle cerebral artery occlusion in human and rat brain.
Our results demonstrated a significant difference in the number of genes affected and the time-course of expression between the two cases. The total number of deregulated genes in the rat was 335 versus 126 in the human, while, of 393 overlapping genes between the two array sets, 184 were changed only in the rat and 36 in the human with a total of 41 genes deregulated in both cases. Interestingly, the mean fold changes were much higher in the human. The expression of novel genes, including p21-activated kinase 1 (PAK1), matrix metalloproteinase 11 (MMP11) and integrase interactor 1, was further analyzed by RT-PCR, Western blotting and immunohistochemistry. Strong neuronal staining was seen for PAK1 and MMP11.
Our findings confirmed previous studies reporting that gene expression screening can detect known and unknown transcriptional features of stroke and highlight the importance of research using human brain tissue in the search for novel therapeutic agents.
Ischaemic stroke results from obstruction of blood flow in a major cerebral vessel and leads to deregulation of genes whose expression promotes ischemic neuronal death and subsequent neurological dysfunction [
The precise molecular mechanisms involved in ischemia-induced brain injury remain poorly understood. Limited knowledge of the molecular mechanisms involved in tissue regeneration has been gained from animal experiments using the middle cerebral artery occlusion (MCAO) model which replicates, in many aspects, the neuropathological changes following stroke in humans [
Analysis of ischemic brain tissue with techniques capable of studying multiple transcripts simultaneously can identify gene expression changes previously not known to be implicated in ischemic pathophysiology and may lead to development of new targets for stroke therapy [
Previous studies employing microarray approaches to study stroke
| Soriano et al. 2000 | Jin et al. 2001 | Kim et al. 2002 | Rao et al. 2002 | Schmidt-Kastner et al. 2002 | Tang et al. 2002 | Roth et al. 2003 | Kim et al. 2004 | Lu et al. 2004 | Moore et al. 2005 | Ford et al. 2006 | Tang et al. 2006 | Vikman and Edvinsson 2006 | Our data | |
| Material used | Rat brain tissue | Rat brain tissue | Rat brain tissue | Rat brain tissue | Rat brain tissue | Rat brain tissue | Rat brain tissue | Rat brain tissue | Rat brain tissue | Rat brain tissue | ||||
|
|
|
|||||||||||||
| Model of ischemia | Permanent focal MCAO | Transient global MCAO | Permanent focal MCAO | Transient focal MCAO | Transient focal MCAO | Permanent focal MCAO | Permanent focal MCAO | Transient focal MCAO | Transient focal MCAO | Blood from ischemic stroke patients | Permanent and transient focal MCAO | Blood from ischemic stroke patients | Post-mortem brain tissue from 11 stroke patients | Post-mortem brain tissue from 12 stroke patients and permanent focal rat MCAO |
|
|
||||||||||||||
| No of genes | 750 | 374 | 1176 | 1263 | 9044 | ~8,000 | ~13,000 | 5,000 | 1,322 | ~19,000 | 8784 | ~39,000 | 7458 | 1176 |
|
|
||||||||||||||
| Time after ischemia | 3 hours | 4 hours |
6 hours | 6 hours |
5 hours | 24 hours | 1 hours |
3 hours |
30 min |
As soon as possible after hospitalization | 24 hours | 3 hours |
7–10 days (obtained 2–3 days post-mortem) | 1 hour-21 days (rat) and 2–37 days (human, obtained by 6 hours post-mortem) |
|
|
||||||||||||||
| Cut-off values | 2.0-fold | 1.7-fold | 2.0-fold | 2.5-fold | 1.7-fold | 2.0-fold | 3.0-fold | 2.0-fold | 2.0-fold | - | 2.0-fold | 1.5-fold | - | 2.0-fold |
|
|
||||||||||||||
| Confirmation of results | In situ hybridization, western blotting | Western blotting, immuno-histochemistry | RT-PCR | Real-time PCR, antisense knockdown, western blotting, immuno-histochemistry | Microarray analysis only | Real-time RT-PCR | Cell culture, in situ hybridization, western blotting, immuno-fluorescence | Cell culture, northern blotting, RT-PCR, western blotting, immuno-histochemistry | Real-time RT-PCR | Real-time RT-PCR | Microarray analysis only | Microarray analysis only | Real-time PCR, immuno-histochemistry | Cell culture, RT-PCR, western blotting, immuno-histochemistry, immuno-fluorescence |
|
|
|
|
||||||||||||
| Selected molecules | NGFI-C |
GRB2 |
IFN-IP |
SOCS-3 | NARP |
PC4 | FAK | Synaptic proteins | CD14 |
LY64 |
PAK1 |
|||
Until recently, gene expression profiling had not been applied to patients dying of ischemic stroke, in part because human brain autopsies are not regularly obtained. Although tissue obtained from brain autopsies is generally of lower quality than that of brain biopsies obtained from living patients, the majority of RNA transcripts and proteins in the human brain are reasonably stable (compared to other tissues such as blood and kidney) and degrade to only a minor degree following death, thus making autopsy tissue a useful source for the isolation of nucleic acids and proteins [
To identify the genes whose expression was changed in the human brain following ischaemia, we investigated the dynamic changes in gene expression in brain samples (collected within 6 h of death) from patients with various times of survival (2–37 days; Table
Clinical Details of Patients
| Patient no. | Age/sex | Survival after stroke | NIHSS on admission | Hypertensiona | Coronary artery disease | Atrial fibrillation | History of TIA/previous stroke | Hypercholesterolemiab | Smoking | Obesityc | Cause of death | Antiplatelets | Statinsd | RSA-be |
| 1 | 63/F | 2 days | 26 | Yes | No | No | No | No | No | No | Large ischemic stroke | No | No | Yes |
| 2 | 84/M | 3 days | 21 | Yes | Yes | No | No | No | Yes | No | Malignant stroke | No | No | No |
| 3 | 68/M | 3 days | 24 | Yes | Yes | No | No | Yes | No | No | Brain oedema | No | No | Yes |
| 4 | 84/M | 6 days | 22 | Yes | Yes | No | No | No | No | No | Cardiac failure | Yes | Yes | No |
| 5 | 51/M | 9 days | 25 | Yes | No | No | Yes | No | No | No | Respiratory infection | Yes | No | No |
| 6 | 74/M | 15 days | 22 | Yes | Yes | No | No | No | No | Yes | Heart attack | Yes | No | Yes |
| 7 | 86/M | 15 days | 14 | Yes | Yes | No | No | Yes | No | No | Urinary infection | Yes | No | No |
| 8 | 58/M | 17 days | 16 | Yes | Yes | No | No | Yes | Yes | No | Cardiac infarction | Yes | No | No |
| 9 | 74/M | 20 days | 12 | Yes | Yes | No | No | No | No | No | Bronchial aspiration | Yes | No | No |
| 10 | 73/M | 26 days | 14 | Yes | Yes | No | No | Yes | No | Yes | Respiratory infection | Yes | No | Yes |
| 11 | 75/M | 29 days | 20 | No | Yes | Yes | No | No | Yes | No | Septic shock | Yes | No | Yes |
| 12 | 60/F | 37 days | 18 | Yes | Yes | No | No | Yes | Yes | No | Pulmonary embolism | Yes | No | No |
a Blood pressure greater than 135/85 mmHg.
b Serum total cholesterol levels greater than 5.2 mmol.
c Body mass index greater than 30.
d Patients who were on statins before the ischemic stroke.
e Patients taking either angiotensin converting enzyme inhibitors or angiotensin type I receptor antagonists.
M = male; F = female; NIHSS = NIH Stroke Scale; TIA = Transient Ischaemic Attack.
The expression of ischemia-related genes was determined by comparing the infarct-induced expression (combined samples from infarcted and peri-infarcted areas) to that in the contralateral hemisphere: 77, 92 and 15 genes were de-regulated in stroke-affected regions in the 3 patient survival groups respectively, while 9, 51, 48, 166, 253, 117 and 261 genes were altered at the 7 different time-points in the animal model compared to the controls (Figure
In total, 126 genes were deregulated after stroke in humans and 335 in the rat MCAO model. However, these data are not directly comparable since many transcripts in the human array were not present in the rat array and
Genes deregulated in both human and animal stroke microarrays
| Human | Rat | Human | Rat | |||||
|
|
||||||||
| Gene name | GenBank | SwissProt | GenBank | SwissProt | Max/min | Days | Max/min | Time |
| c-jun proto-oncogene |
|
|
|
|
4.4-fold | 9 – 20 | 3.5-fold | 1 h – 24 h |
| Matrix metalloproteinase 11 |
|
|
|
|
3.2-fold | 2 – 20 | 2.6-fold | 3 days |
| Calcium/calmodulin-dependent kinase (CAMK1) |
|
|
|
|
17.2-fold | 2 – 20 | 0.05-fold | 21 days |
|
|
|
|||||||
| LIM domain kinase 1 |
|
|
|
|
3.6-fold | 2 – 20 | 2.4-fold | 3 days |
| 0.4-fold | 21 days | |||||||
| T-Lymphocyte maturation-associated protein |
|
|
|
|
1.7-fold | 2 – 6 | 0.2-fold | 21 days |
| Retinoic Acid Receptor beta |
|
|
|
|
2.0-fold | 2 – 6 | 0.1-fold | 21 days |
| S54072 |
|
|||||||
| Tyrosine Phosphatase 1B |
|
|
|
|
3.4-fold | 2 – 6 | 0.2-fold | 21 days |
| Adenosine A1 Receptor |
|
|
|
|
2.6-fold | 2 – 6 | 5.2-fold | 4 hrs |
| Growth arrest & DNA damage-inducible 153 |
|
|
|
|
2.4-fold | 2 – 6 | 2.1-fold | 3 days |
|
|
||||||||
| Glutamate Decarboxylase 67 |
|
|
|
|
5.6-fold | 2 – 6 | 2.5-fold | 21 days |
| Glutamate Decarboxylase 65 |
|
|
|
|
22.7-fold | 2 – 20 | 2.2-fold | 3 days |
| Neurotrophin 3 |
|
|
|
|
5.1-fold | 2 – 37 | 2.2-fold | 12 hrs |
| Inhibitor of DNA binding 2 |
|
|
|
|
5.6-fold | 2 – 20 | 0.4-fold | 21 days |
| Neuropeptide Y |
|
|
|
|
8.8-fold | 2 – 20 | 0.04-fold | 21 days |
| Glia Maturation Factor beta |
|
|
|
|
7.6-fold | 2 – 6 | 0.04-fold | 21 days |
| High Mobility Group Protein 1 |
|
|
|
|
4.3-fold | 2 – 37 | 3-fold | 4 h – 3 d |
|
|
0.3-fold | 21 days | ||||||
| Early Growth Response Protein 1 |
|
|
|
|
4.4-fold | 2 – 20 | 3.9-fold | 1 h – 12 h |
|
|
|
0.2-fold | 21 days | |||||
| TAT-Binding Protein 1 |
|
|
|
|
3.8-fold | 2 – 20 | 0.4-fold | 21 days |
| Glutathione S-Transferase 1 |
|
|
|
|
17.5-fold | 2 – 20 | 10.8-fold | 24 h – 21 d |
| Fibroblast Growth Factor Receptor 1 |
|
|
|
|
10.1-fold | 2 – 20 | 4-fold | 4 h – 24 h |
|
|
|
|||||||
|
|
||||||||
| Interleukin 10 |
|
|
|
|
2.4-fold | 2 – 20 26 – 37 | 6.4-fold | 21 days |
|
|
0.2-fold | |||||||
| Heat Shock Protein 27 |
|
|
|
|
0.6-fold | 2 – 20 | 15.2-fold | 4 h – 24 h |
| Heat Shock Protein 70 |
|
|
|
|
0.6-fold | 2 – 6 | 9.4-fold | 1 h – 24 h |
|
|
||||||||
| Thioredoxin Peroxidase 1 |
|
|
|
|
4.9-fold | 2 – 20 | 3.9-fold | 21 days |
|
|
|
|||||||
| Platelet-Derived Growth Factor A |
|
|
|
|
1.6-fold | 2 – 6 | 0.5-fold | 21 days |
| Matrix Metalloproteinase 14 |
|
|
|
|
6.9-fold | 2 – 6 | 3.3-fold | 24 h – 3 d |
| Kinase receptor TYRO3 Sky proto-oncogene |
|
|
|
|
3.1-fold | 9 – 20 | 4.0-fold | 24 h – 3 d |
| 0.4-fold | 21 days | |||||||
| CSF-1-Receptor |
|
|
|
|
89.2-fold | 9 – 20 | 2.8-fold | 3 days |
| Insulin-like Growth Factor Binding Protein 2 |
|
|
|
|
79.7-fold | 9 – 20 | 2.1-fold | 3 days |
| Mitogen activated kinase 1/2 |
|
|
|
|
48.6-fold | 9 – 20 | 0.3-fold | 21 days |
| Aquaporin 4 |
|
|
|
|
18-fold | 9 – 20 | 3.2-fold | 3 days |
| erbB2 proto-oncogene Neu proto-oncogene |
|
|
|
|
11.2-fold | 9 – 20 | 2.7-fold | 12 hrs |
|
|
|
|||||||
| L-type calcium channel β3 |
|
|
|
|
10.3-fold | 9 – 20 | 8.6-fold | 21 days |
| Ras-related protein RAB3A |
|
|
|
|
13.9-fold | 9 – 20 | 0.3-fold | 21 days |
| CAMK-II beta |
|
|
|
|
1.8-fold | 9 – 20 | 0.3-fold | 21 days |
| Growth Factor Receptor-Bound 2 |
|
|
|
|
19.9-fold | 9 – 20 | 2.7-fold | 3 days |
|
|
|
|
||||||
| Signal Transducer & Activator of Transcription 3 |
|
|
|
|
0.4-fold | 9 – 20 | 6.6-fold | 4 h – 3 d |
| 0.05-fold | 21 days | |||||||
| Neuronatin |
|
|
|
|
11.1-fold | 9 – 20 | 0.4-fold | 21 days |
|
|
||||||||
| Glutathione S-Transferase P |
|
|
|
|
3.1-fold | 9 – 20 | 0.1-fold | 1 hr |
| Glucocorticoid-regulated serine/threonine kinase GSK |
|
|
|
|
0.6-fold | 26 – 37 | 2.4-fold | 3 days |
| 0.05-fold | 21 days | |||||||
| Glucose Transporter 1 |
|
|
|
|
0.6-fold | 26 – 37 | 11.6-fold | 4 h – 21 d |
Amongst these genes we examined in more detail a small subset with no prior report of a role in stroke (PAK1, MMP11 and INI1). PAK1 was only induced in the human study although present in both microarray sets, MMP11 was induced in both cases, while INI1 was induced in the human but was not present in the rat microarray set. To confirm the microarray data, RT-PCR was carried out on selected deregulated genes. The temporal expression patterns of these genes following RT-PCR showed good agreement with the corresponding expression profiles obtained from the microarray analysis, supporting the validity of the data obtained from the microarrays. Using Western blotting and immunohistochemistry, PAK1, INI1 and MMP11 protein expression and localization was determined in the contralateral and ipsilateral brain areas of individual stroke patients and rats subjected to MCAO, and in HBMEC and HFN exposed to OGD and reperfusion.
In agreement with the microarray data, RT-PCR demonstrated an increase in
Protein expression in infarcted (I) and peri-infarcted (P) areas (Fold increase compared to contralateral hemisphere)
| PAK1 | INI1 | MMP11 | |||||
|
|
|||||||
| Patient no. | Survival (days) | P | I | P | I | P | I |
| 1 | 2 | 2.2 | 1.0 | 1.5 | 1.5 | 1.5 | 1.5 |
| 2 | 3 | 3.3 | 4.0 | 0.2 | 0.4 | 1.0 | 1.0 |
| 3 | 3 | 1.0 | 1.0 | 4.2 | 4.3 | 0.7 | 0.7 |
| 4 | 6 | 1.0 | 1.0 | 4.3 | 5.8 | 1.6 | 1.5 |
| 5 | 9 | 1.5 | 0.4 | 3.2 | 3.3 | ND | ND |
| 6 | 15 | 2.3 | 1.5 | 2.8 | 1.0 | 1.7 | 1.6 |
| 7 | 15 | 3.0 | 3.2 | 1.7 | 2.0 | 1.0 | 1.0 |
| 8 | 17 | 1.0 | 1.0 | 1.0 | 1.0 | 1.0 | 1.0 |
| 9 | 20 | 1.0 | 1.0 | 1.0 | 1.0 | 1.0 | 1.0 |
| 10 | 26 | 1.5 | 1.5 | 1.0 | 1.6 | 5.1 | 2.2 |
| 11 | 29 | 1.0 | 1.5 | 2.2 | 2.8 | 1.8 | 3.5 |
| 12 | 37 | 1.5 | 1.0 | 1.7 | 1.7 | 1.0 | 1.5 |
| Total | Upregulated | 7 | 5 | 8 | 8 | 5 | 6 |
| Downregulated | 0 | 1 | 1 | 1 | 1 | 1 | |
| No change | 5 | 6 | 3 | 3 | 5 | 4 | |
| No detection | 0 | 0 | 0 | 0 | 1 | 1 | |
For MMP11, RT-PCR data agreed with the findings from the microarray study, showing increased mRNA levels in infarcted and peri-infarcted tissue from patients surviving 2–20 days following stroke (Figure
RT-PCR confirmed the upregulation of
In the human brain, many differentially expressed genes were observed from 2 to 6 days and from 9 to 20 days after stroke, with the majority being upregulated. The number of deregulated genes declined during 26 to 37 days after stroke, indicating that dynamic changes in gene expression occur during the first days to few weeks in the human postischaemic brain. In the rat brain, few differences were observed at 1 hour, while the number of differentially expressed genes steadily increased with time after MCAO, with a peak after 3 days, supporting the concept of active mechanisms initiated during the acute phase after experimental stroke and lasting for several days. The number of upregulated genes gradually increased, peaking at 3 days, while downregulated genes were detected 24 h after MCAO and increased dramatically until the final measured time-point at 21 days (Figure
The limitations of post-mortem brain samples in cDNA microarray analysis concern the small sample size and potential low quality and the genetic heterogeneity and diversity in terms of age, sex and previous medical history within a group of patients [
The present gene expression profile study is the first large-scale microarray report showing altered expression of several genes following human stroke. These included genes participating in transcription, apoptosis, inflammation and neuroprotection. Many genes/proteins previously shown to be deregulated following stroke were reported in our study too e.g. IL-10 [
PAK1 is a downstream Rac effector and a major cyclin-dependent kinase 5 (Cdk5) substrate and target that co-localizes with p35/Cdk5 at neuronal peripheries. P35/Cdk5 causes PAK1 hyperphosphorylation, which results in PAK1 down-regulation and is likely to have an impact on the dynamics of the reorganization of the actin cytoskeleton in neurons during dendrite development [
MMP11 or stromelysin-3 (ST3), first isolated as a breast cancer-associated protease, is not expressed in the majority of normal adult organs but is expressed during a number of pathological processes, including wound healing and atherosclerotic lesions [
INI1 is a tumour suppressor gene, thought to exert its tumour suppressor function by mediating cell cycle arrest [
Many experimental trials of stroke therapies have failed to translate to human clinical trials and one possible way to improve the success rate can be through comparative genomics. As it has been recently commented, it is very surprising that the exciting developments observed in basic and clinical stroke research over the past two decades have occurred in parallel, with too little direct translation between bench and bedside [
Human brain tissue samples were obtained from 12 patients who died from acute ischaemic stroke, with the approval of the local Ethics Committee and Brain Bank at the Department of Neuropathology, Collegium Medicum, Jagiellonian University, Krakow, Poland. All patients were admitted with large middle cerebral artery strokes confirmed by CT-scan or MRI. The patients, 10 male and 2 female, were aged between 51 and 86 years and had survived between 2–37 days following ischaemic stroke (Table
Stroke experiments were performed on female Sprague-Dawley rats (weight: 230–270 g) as they suffer less than male during ischaemia. Cerebral ischaemia was produced using a modified method of Baron [
Human brain microvascular endothelial cells (HBMEC) were obtained from TCS CellWorks (Buckingham, UK) and cultured according to the supplier's instructions. Human foetal (cerebral cortical) neurons (HFN) were extracted and cultured with permission from the Local Ethics Committee. Brain tissue from foetus specimens of 14–19 weeks gestational age, legally aborted and with the appropriate written consent, were collected in cold preservation medium and cells were isolated and cultured as described elsewhere [
We established mRNA expression profiles of the damaged brain tissues between 2 to 6 days, 9 to 20 days, and 26 to 37 days after stroke in human patients and 1, 4, 12, 24 hours and 3, 7 and 21 days after the ischaemic insult in rats. The corresponding samples from the non-ischemic control hemisphere were used to measure the normal mRNA abundance of the modulated genes in each tissue at each time point. RNA from three stroke patients was pooled for each patient survival group while RNA from three MCAO rats was also pooled at each time-point to improve yields in preparation of poly A+ RNA. Although pooling was previously thought to affect data quality, Kendziorski
RNA was extracted according to the manufacturer's protocols (BD Biosciences, Oxfordshire, UK) and its quality was measured spectrophotometically. The protocol recommended by Clontech in their Atlas 1.2 microarray kit was used without any modification. Briefly, RNA was reverse-transcribed to cDNA, 32P-labelled and applied to the array for overnight hybridisation at 68°C. Following washing, the array was exposed to a phosphorimaging plate for 12–72 hours and data analysis was performed using the AtlasImage 1.5 software. The results were normalized using two housekeeping genes, ubiquitin and glyceraldehyde 3-phosphate dehydrogenase (GAPDH). As in the majority of microarray studies mentioned before, only those genes upregulated > 2-fold or downregulated < 0.5-fold were counted as deregulated and taken into consideration. The microarray data are available in Gene Expression Omnibus under the accession number GSE9391.
Gene expression was examined by semi-quantitative RT-PCR with standard reaction conditions of a 10 min denaturation at 94°C, followed by 35 cycles of 1 min at 94°C, 1 min at primer-specific annealing temperatures (Table
Primer sequences
|
|
|
|
|
|
|
Human | 5'-TAAAGGTATGGAGCGATGTGAC-3' (forward) | 58°C |
|
|
5'-TGGGTAGCGAAAGGTGTAGAAG-3' (reverse) | ||
|
|
Rat | 5'-GATGGAGGCCAGCTAGTCAG-3' (forward) | 60°C |
|
|
5'-ATGGTACATGACCACGCAGA-3' (reverse) | ||
|
|
Human | 5'-ACCCTGTCCAACAGCTCCCA-3' (forward) | 64°C |
|
|
5'-GGCCCAATCTTCTGAGATGC-3' (reverse) | ||
|
|
Rat | 5'-CCTGGGGCTCCTATACAAAA-3' (forward) | 60°C |
|
|
5'-CCATGACCGAGCAAATGAC-3' (reverse) | ||
|
|
Human | 5'-GCTGTTCTGGATGTGTTGGA-3' (forward) | 60°C |
|
|
5'-TCTGCTCTGGGGTTATCTGTG-3' (reverse) | ||
|
|
Rat | 5'-AGCAAAAGAGGCAACCAAGA-3' (forward) | 60°C |
|
|
5'-GGGTAAGGAATGGGATGGTT-3' (reverse) | ||
|
|
Human | 5'-ATGATCTTGAGGCTGTTG-3' (forward) | 58°C |
|
|
5'-CTCAGACACCATGGGGAA-3' (reverse) |
Proteins were extracted from tissues and the protein concentration of each sample was determined using the BioRad assay. For Western blotting, 10 μg of protein were separated by SDS-PAGE (13% w/v) and the proteins were electro-blotted onto nitrocellulose filters as described previously [
Paraffin-embedded tissue samples were processed and serial 5 μm sections were cut. The Avidin-Biotin Peroxidase (ABC Vectastain kit, Vector Laboratories, Peterborough, UK) method was used and antibodies to MMP11, PAK1 and INI1 were used at a dilution of 1:50. Paraffin-embedded sections were deparaffinized, rehydrated and boiled for 10 min in an antigen unmasking solution of concentrated citric acid pH 6.0 as described elsewhere [
GAPDH: glyceraldehyde 3-phosphate dehydrogenase; HBMEC: human brain microvascular endothelial cells; HFN: human foetal neurons; INI1: integrase interactor 1; MCAO: middle cerebral artery occlusion; MMP11: matrix metalloproteinase; OGD: oxygen-glucose deprivation; PAK1: p21-activated kinase 1.
NM carried out the human studies and drafted the manuscript. MOS carried out the experimental work in rats. JK and FR provided the material for the study. RP and CS carried out the
This work was supported by the Higher Education Funding Council for England (HEFCE) and the Research Institute for Health and Social Change (RIHSC). We also thank Mick Hoult for his assistance in the preparation of this manuscript.