Mapping similarities in mTOR pathway perturbations in mouse lupus nephritis models and human lupus nephritis
- Padmalatha S Reddy†1,
- Holly M Legault†1,
- Joseph P Sypek2,
- Mark J Collins3,
- Elizabeth Goad4,
- Samuel J Goldman2,
- Wei Liu1,
- Stuart Murray1,
- Andrew J Dorner5 and
- Margot O'Toole1Email author
© Reddy et al.; licensee BioMed Central Ltd. 2008
Received: 20 June 2008
Accepted: 03 November 2008
Published: 03 November 2008
The Erratum to this article has been published in Arthritis Research & Therapy 2009 11:407
Treatment with sirolimus, a mammalian target of rapamycin (mTOR) inhibitor, has been shown to be efficacious in the MRL/lpr and NZB × NZW F1 mouse models of lupus nephritis, indicating a critical role for the mTOR pathway in both models. This type of demonstration of efficacy in animal models is usually a pre-requisite for advancement into clinical development. However, efficacy in an animal model often has not translated to the desired activity in the clinic. Therefore, a more profound understanding of the mechanistic similarities and differences between various animal models and human diseases is highly desirable.
Transcriptional profiling was performed on kidneys from mice with lupus nephritis; from mice who had efficacious drug treatment; and from mice before they developed nephritis. Analysis of variance with false discovery rate adjusted to p < 0.05 and an average fold change of two or more was used to identify transcripts significantly associated with disease and response to therapy. Pathway analyses (using various bioinformatics tools) were carried out to understand the basis for drug efficacy in the mouse model. The relevance in human lupus of the pathways identified in the mouse model was explored using information from several databases derived from the published literature.
We identified a set of nephritis-associated genes in mouse kidney. Expression of the majority of these returned to asymptomatic levels on sirolimus treatment, confirming the correlation between expression levels and symptoms of nephritis. Network analysis showed that many of these nephritis genes are known to interact with the mTOR pathway. This led us to ask what human diseases are linked to the mTOR pathway. We constructed the mTOR pathway interactome consisting of proteins that interact with members of the mTOR pathway and identified a strong association between mTOR pathway genes and genes reported in the literature as being involved in human lupus.
Our findings implicate the mTOR pathway as a critical contributor to human lupus. This broad pathway-based approach to understanding the similarities in, and differences between, animal models and human diseases may have broader utility.
Clinical development of therapies is heavily dependent on demonstrated efficacy in animal model(s), but efficacy in animal models often does not translate into clinical success. A number of factors have been proposed as contributing to this lack of concordance between efficacy in animal and clinical studies [1–5]. One clear limitation of relying on disease models in inbred strains is that the genes that produce the disease phenotype in a given model may represent only a subset of the genes that can cause the phenotype in complex human diseases such as lupus. Using our own animal model transcriptomics, the vast and rapidly accumulating literature on genes linked to human disease and pathway tools, we have taken a broad analytical approach to identifying similarities between the mouse and human 'lupus phenotype' at the level of biological pathway perturbations. The potential advantage of this approach is that, by linking the human disease phenotype to a pathway, drug development efforts can be targeted to the pathway. Animal models with involvement of the same pathway can then be chosen and/or derived.
Systemic lupus erythematosus (SLE) is a chronic inflammatory autoimmune disease [6–8]. The pathophysiology of disease is manifested by the production of autoantibodies directed against multiple self-antigens. This dysregulation of the immune system resulting in the loss of tolerance appears to be mediated by both T cells and B cells. Many organs including the kidney can be affected . Direct action of autoantibodies, deposition of immune complexes and pro-inflammatory cytokines, particularly interferon (IFN) γ, have all been implicated in disease pathophysiology [10–13].
There are at least four mouse models of lupus nephritis . Both NZB × NZW F1  and MRL/lpr mouse [16, 17] strains spontaneously develop autoimmune lupus nephritis. Female mice from the NZB × NZW F1 cross (NZB/W) develop proteinuria and only a small number (< 20%) survive to 52 weeks. In MLR/lpr mice, the disease develops in both males and females and is associated with the fas lpr mutation on the MLR background. Mice develop significant proteinuria at 16 weeks and show significant mortality rates (about 50%) by 20 weeks [18, 19]. Despite their independent derivation, lupus nephritis in both MLR/lpr and NZB/W mouse models shows a remarkably efficacious response to sirolimus treatment [20–22].
Sirolimus (rapamycin) is an immunosuppressive drug that binds to mTOR (mammalian target of rapamycin), a serine/threonine kinase that regulates cellular proliferation and metabolism and blocks G1 to S phase cell cycle progression, interfering with T and B cell activation [23–25]. Sirolimus is approved for the prevention of transplant rejection . We used our own data (presented here) and previously published data on the efficacy of mTOR inhibitors in two mouse models of lupus nephritis to infer that perturbations of the mTOR pathway are critical to the development of lupus nephritis in both these models. In order to assess the likelihood of mTOR pathway involvement in human lupus, we examined the concordance between the mTOR pathway interactome and genes linked to human lupus and report the results of this analysis here.
Materials and methods
NZB/W (NZB × NZW F1 cross) females were purchased from the Jackson Laboratory (Bar Harbor, ME). These mice were maintained and studied under pathogen free conditions in accordance with guidelines from the American Association for the Accreditation of Laboratory Animal Care and the Institutional Animal Care and Use Committee of Wyeth Research.
Beginning at 20 weeks of age, disease progression was monitored weekly by assessing proteinuria. A cohort of mice, selected at 20 weeks of age, served as the 'asymptomatic' group (n = 4). Once fixed proteinuria of 30 to 100 mg/dL had appeared on two consecutive occasions, (generally at 25 to 29 weeks of age), the diseased mice were randomly assigned to either the sirolimus-treated group (n = 6 mice) or the untreated (disease) group (n = 6). Sirolimus, dissolved in carboxymethylcellulose (vehicle), was subcutaneously administered three times weekly in single doses of 1 mg/kg or 5 mg/kg for eight weeks. Mice were monitored weekly until 52 weeks of age.
Assessment of proteinuria
Urine was manually expressed from each mouse on a weekly basis, collected into a sterile container and assayed for the presence of protein (specifically albumin) using a colorimetric method (Albustix Reagent Strips, Bayer Corporation, Elkhart, IN).
Proteinuria evaluations were scored as follows: grade 0.5 = 'trace' proteinuria; grade 1 = about 30 mg/dL; grade 2 = about 100 mg/dL; grade 3 = about 300 mg/dL; and grade 4 = more than 2000 mg/dL. If mice achieved a grade 4 reading on two consecutive days they were euthanised.
Assessment of renal pathology
Kidneys were harvested from mice one to four months after an eight-week course of treatment. Three to five mice were examined in each group. One-half of a kidney was fixed by overnight immersion in 10% formaldehyde and paraffin-embedded. The other half was snap frozen for RNA preparation. To determine the extent of renal damage, sections were stained with H & E and periodic acid-Schiff (PAS) and scored for pathological changes. In addition, glomerulopathy was scored on a 0 to 5 scale. Severity grades were as follows: 0 = normal or within normal limits; 1 = minimal or slight; 2 = mild; 3 = moderate; 4 = marked; 5 = severe.
RNA purification and microarray hybridisation
Snap frozen murine kidney tissue was homogenised in RLT buffer containing 1% beta-mercaptoethanol using the polytron and RNA purified by Qiagen RNeasy columns (Qiagen, Valencia, CA). Eluted RNA was quantified using a Spectramax 96-well plate UV reader (Molecular Devices, Sunnyvale, CA, USA) monitoring A260/280 OD values. The quality of each RNA sample was assessed by capillary electrophoresis alongside an RNA molecular weight ladder on the Agilent 2100 bioanalyser (Agilent Technologies, Palo Alto, CA, USA).
Five micrograms of total kidney RNA was prepared from individual NZB/W F1 mice of the following groups; untreated mice at 12 weeks (n = 4); untreated F1 mice at 36 and 42 weeks combined (n = 6); and sirolimus-treated mice at 36 and 42 weeks combined (n = 6). The animals selected for expression analysis reported here were representative of multiple studies performed that confirmed the data on proteinuria, mortality rates and histopathology here.
Biotin-labelled cRNA was prepared using an oligo T7-primed reverse transcription reaction followed by in vitro transcription reaction with biotin-labelled UTP and CTP. cRNA 15 μg was fragmented and hybridised to Mu11KsubA and Mu11KsubB arrays (Affymetrix, Santa Clara, CA). Hybridised arrays were washed and stained with Streptavidin R-phycoerythrin (Molecular Probes Inc., Eugene, OR) using the GeneChip Fluidics Station 400 (Affymetrix, Santa Clara, CA) and scanned with a Hewlett Packard GeneArray Scanner (Hewlett Packard Company, Palo Alto, CA) according to the manufacturer's protocols. All array images were visually inspected for defects and quality. Arrays with excessive background, low signal intensity or major defects within the array were eliminated from further analysis. GeneChip MAS 5.0 software was used to evaluate the hybridisation intensity, compute the signal value for each probe set and make an absent/present call. GeneChip signal data for the samples analysed in this study are available under Accession Number [GEO:GSE12924] .
Data normalisation and filtering
GeneChips were required to pass standardised quality control criteria. RNA quality was monitored by the ratio of frequencies measured by independent probe sets representing 5' and 3' regions of glyceraldehyde 3-phosphate dehydrogenase. This ratio must be more than 0.4. Filtering criteria for individual probe sets required that a probe set was called 'present' or a signal of 50 or more in at least one of the samples. All filtering criteria were passed by 6384 probe sets and were subject to the statistical analysis described below, and probe sets that did not meet these criteria were not included in subsequent analyses.
For hierarchical clustering of probe sets and arrays, the Log-2 scale MAS5 expression values from each probe set were first z-normalised so each probe set had a mean expression level of zero and a standard deviation (SD) of one across all samples. Then these normalised profiles were clustered hierarchically using an unweighted paired group method with arithmetic mean, and the Euclidean distance measure.
Identification of genes associated with lupus nephritis and response to sirolimus therapy
The disease-related fold change differences were calculated by determining the difference in the log-2 signal of the 12-week-old asymptomatic mice and the combined 36 and 42-week-old diseased mice. Analysis of variance (ANOVA) was performed using this metric to identify disease-related differences. Raw p values were adjusted for multiplicity of testing using the false discovery rate (FDR) procedure of Reiner and colleagues  using Spotfire (Somerville, MA). Genes with a FDR p < 0.05 and an absolute fold change of two of more in the comparison between disease and asymptomatic groups were identified as being significantly associated with lupus nephritis. Lupus nephritis genes were identified as being significantly associated with response to sirolimus treatment if they met an FDR p < 0.05, in comparison between sirolimus-treated and disease groups. Those that failed to meet the FDR p < 0.05 in the comparison of the sirolimus-treated group to the disease group, did have a significant difference (FDR p < 0.05) in the comparison between sirolimus-treated and asymptomatic groups, confirming a resistance to sirolimus therapy.
Pathway analysis was performed using Ingenuity Pathways Analysis (IPA) (Ingenuity Systems, Redwood City, CA), MetaCore (GeneGo Inc., St. Joseph, MI) and an in-house implementation of the sigPathway algorithm . SigPathway is an algorithm that identifies differentially expressed gene sets. An FDR p < 0.01 was used to identify significantly changing gene sets. Networks were built using genes that met FDR p < 0.05 within a gene set. Networks, canonical pathways and functional processes for genes passing either the FDR p < 0.05 and an absolute fold change of two or more criteria and/or SigPathway filter FDR p < 0.05 were analysed using IPA. Human, mouse and rat lupus-associated genes were identified using the search tools within IPA and MetaCore. Rapalog-mTOR pathway and its connectivity to the lupus nephritis disease data and the lupus disease genes were explored using the pathway building tools in IPA. The rapalog-mTOR pathway interactome was built using IPA and all proteins were exported to MetaCore to explore the human disease representation on the mTOR pathway interactome. In MetaCore, a gene is considered to be associated with a condition (disease, pathological or toxic process) as a biomarker if this gene or its product (RNA or protein) has different properties in the disease and healthy states. Such properties may include DNA characteristics (eg, mutation, single nucleotide polymorphisms, chromosomal rearrangement); epigenetics (eg, methylation of a promoter region, histone acethylation, nucleosome misplacement); RNA level (eg, disease specific splice variant, higher/lower expression level); or protein level (eg, disease-specific isoform, mutant isoform, subcellular localisation, post-translational modification, abundance). These genes are used to generate human 'disease biomarkers' networks using direct interactions between biomarker genes and proteins from their MetaBase database. I2E (from Linguamatics Ltd., Cambridge, UK) and clinical database findings were used to validate some of the relationships between sirolimus analogues and various human diseases identified in the mTOR pathway interactome.
Short course therapy with sirolimus prevents onset of murine lupus and renal damage
Histopathology evaluation of kidneys from NZB/W F1 hybrid female mice at 42 weeks of age
p value, untreated versus treated
Identification of the disease-associated transcriptome
Identification of nephritis-associated probe sets modulated by sirolimus treatment
Biological annotation of disease and drug-responsive genes
A critical signalling pathway involved in mediating an inflammatory response is the JAK-STAT and MAP kinase pathway. Increased levels of transcripts for pathway members including JAK3, STAT3, SOCS3, PTPN1, CDKN1A, RRAS and MAPK1 were observed in nephritis. After treatment with sirolimus these pathways exhibit transcriptional expression levels similar to asymptomatic levels.
Rapalog-mTOR canonical pathway and links to mouse lupus nephritis genes
Networks were built in an effort to understand the broad ranging beneficial effects of the mTOR pathway inhibitor, sirolimus, in NZB/W lupus nephritis. The first step in this process was to build a rapalog-mTOR pathway. This pathway consisted of the mTORC1 complex (mTOR, GBL, Raptor), the mTORC2 complex (mTOR, GBL, AVO3), the immediate downstream targets of mTOR – RPS6KB1 and RPS6KB2, and the upstream effectors of mTOR – AKT1, AKT2, TSC1, TSC2. In addition, rapalogs, such as sirolimus, temsirolimus and everolimus as well as all members of the immunophilin protein family that bind to the immunosuppressants FK506 and rapamycin, were included in the rapalog-mTOR pathway. Downstream connectivity of the rapalog-mTOR pathway to the 387 lupus nephritis genes was explored using the IPA. Of the 387 genes, 32 can be placed immediately downstream of the rapalog-mTOR pathway. An additional 25 of these are connected to the pathway through various types of functional protein interactions (eg, phosphorylation, binding etc).
Building the mTOR pathway interactome
All proteins known to interact with members of the rapalog-mTOR pathway, upstream and downstream were used to create the mTOR pathway interactome, consisting of 570 proteins (see Additional file 5). The magnitude of this 570 protein interactome is likely to be due in part to some high connectivity proteins such as the AKT gene group, which regulate many pathways. The 570 proteins from this interactome were imported into Metacore to query human diseases significantly associated with the mTOR pathway interactome. Significantly associated human diseases likely to be perturbed by mTOR pathway dysfunction were identified using the functional enrichment category within Metacore that relies on curated human disease specific networks.
Human lupus and genes of the mTOR pathway interactome
Various cancer and non-cancer human diseases were identified as being significantly associated with the mTOR pathway interactome. A few cancers, such as ovarian neoplasm and pancreatic neoplasm, topped the list of associated diseases, and remarkably these were followed by lupus with significance of association p value of about E-8.5. Other non-cancer diseases that showed a less significant association are diabetes, obesity, Alzheimer's disease, multiple sclerosis and arthritis, all having an association of p < 0.01. Notably, the p values for the other pro-inflammation diseases were much less significant than those for lupus. For example, p values for both multiple sclerosis and arthritis are about 1E-3.
The genes in the mTOR pathway interactome that are associated with lupus are, in large part, distinct from the genes involved in the other non-cancer diseases. This is evident from the top scoring Gene Ontology cellular processes for each of these diseases. For example, in Metacore, the top scoring Gene Ontology cellular processes for lupus are related to apoptosis and cell death, while those for diabetes are related to carbohydrate metabolism, and those for obesity are related to fat metabolism. Both diabetes and obesity show a much more significant association with the mTOR interactome than any pro-inflammatory disease other than lupus.
NZB/W mice develop nephritis closely resembling that seen in human patients with lupus nephritis. Here we show that an eight-week course of 5 mg/kg sirolimus delivered three times a week starting at disease onset (about 25 weeks) and continuing for eight weeks markedly reduced symptoms of disease as measured by proteinuria, kidney histopathology and survival. By 52 weeks of age all treated mice were alive despite cessation of treatment many weeks earlier. Both results described here and previously reported studies [20, 22, 30], establish that sirolimus treatment prevents progression of autoimmune nephritis and prolongs survival in NZB/W mice. A comparable effect of sirolimus treatment has been observed in MRL/lpr mice, a second model of lupus and lupus nephritis .
We identified genes expressed at abnormal levels in NZB/W kidneys by comparing RNA levels in asymptomatic young mice and older mice with symptoms of lupus nephritis. These lupus nephritis genes were further shown to not be associated with the normal ageing process based on the observed differences between healthy young and old C57BL6 mice. A broad range of biological functions was represented among the lupus nephritis genes identified in this study. As expected, given the loss of kidney function, the vast majority of genes involved in metabolic pathways are down-regulated in nephritis and, given the inflammatory nature of the disease, many of the signalling pathway genes are up-regulated. Glomerular disease is a significant component in lupus nephritis. A recent study identified a glomerulus-enriched gene set . We used data from this study to determine if the nephritis-associated genes are enriched in the glomerular gene set. We found a highly significant over-representation of the glomerular genes (chi-square value of 49.29 with one degree of freedom) consistent with glomerular involvement.
A recent study by Liu and colleagues reported on 126 nephritis-associated genes in the MRL/lpr model . Of these, 37 genes were present in the nephritis-signature reported here. Commonalities were noted in the nephritis signatures (and this includes the directionality of change in the disease state) between these two models, such as the antigen presentation and complement pathways as well as various IFN-regulated genes and immunoglobulins. A good overlap (21 genes) was also noted between our mouse nephritis gene set and 68 human lupus nephritis genes derived from laser-captured glomeruli from SLE patients . Additional similarities may be present, but probably lie outside the statistical parameters of both datasets.
A profound 'normalisation' of expression levels of lupus nephritis genes was observed in mice treated with sirolimus, both for metabolic as well as signalling pathways. Affected metabolic pathways in lupus nephritis include fatty acid degradation, glycolysis pathways and leucine/valine/isoleucine degradation. Transcripts for BCKDHA and DBT, two enzymes in the branched chain amino acid metabolism pathway required for the catabolism of leucine, valine and isoleucine, are reduced in nephritis, perhaps leading to the accumulation of leucine in diseased tissue. Interestingly, leucine activates the target of sirolimus inhibition, mTOR, leading to increased protein synthesis , and in addition we noted an increase in ribosomal RNA transcripts in the disease state. This physiological connection suggests that mTOR pathway activation may be increased by leucine in disease, providing perhaps an additional mechanism for sirolimus efficacy. Levels of these transcripts were returned to asymptomatic levels in sirolimus-treated mice. Several genes in the mitochondrial electron transport chain are also down-regulated in the disease state, and mitochondrial dysfunction has been implicated in kidney function impairment .
Reflecting the pro-inflammatory functions of nephritis, genes such as JAK3, STAT3 and MAPK1 involved in signalling pathways (JAK/STAT, MAP kinase, antigen presentation, IL1/IL10) are expressed at higher levels in the disease state. Also SOCS3, a negative regulator of JAKs and PTPN1 and CDKN1A, a negative regulator of STATs, are also elevated in the disease state. Although activation of these signalling pathways occurs through phosphorylation-dephosphorylation events of pathway components, it can be noted here that this pathway is also dysregulated at the transcriptional level in lupus nephritis. This complex dysregulation of the JAK/STAT pathway, which drives production of multiple cytokines and other inflammatory mediators, is returned to asymptomatic levels on sirolimus treatment. PTPN1, a negative regulator of STATs, is a notable exception, suggesting a link between the quiescence of this pathway with amelioration of disease. Consistent with the activation of this signalling pathway, genes involved in immune system cascades, such the IFN-regulated genes (IRF1, IRF7, OAS1, IFITM3, IFI27 and IGTP), and signalling by IL2 subfamily of type 1 cytokines (IL2Rγ) were also up-regulated in the disease state and are down-regulated by sirolimus.
Genes of the complement pathway known to be involved in renal damage, such as C3, C4, C1QA, CCL13 and FCGR2a, are also expressed at higher levels than in the untreated group. C3, C4 and C1QA play a role in antigen clearance. Using sigPathway , an algorithm that identifies differentially expressed gene sets, additional components of the complement pathway are transcriptionally elevated in the diseased renal tissue. Our results suggest that the complement components in the early parts of both the classical and alternate pathways are elevated in nephritis, while one component of the membrane attack complex, further downstream in the complement pathway, is down-regulated. C1q and C3, but not C4, were normalised by treatment. Complement pathway components are known to be significant contributors to renal damage. C3 deposition in the kidney has been observed in both human lupus nephritis and in murine models . The elevated levels of C4 during disease amelioration is consistent with the concept that the early members of the classical pathway may be important in reducing disease pathology by clearing immune complexes and apoptotic cells [35, 36]. Our profiling analysis also identified a large number of immunoglobulin transcripts elevated in the kidney tissue consistent with the role of autoantibodies and immune complex deposition in pathology.
To understand the mechanism by which sirolimus normalised such a wide range of biological processes, networks were built around the nephritis genes and the rapalog-mTOR pathway. Using curated findings from the literature, the shortest path for about one sixth of the 387 nephritis genes was defined to be either 0 or one-step downstream of the rapalog-mTOR pathway. This suggests a close functional association of mTOR pathway with disease mechanisms. In the context of the findings reported here, it is worth noting that steroid and cyclophosphamide, known to ameliorate lupus, directly impact some components of the mTOR pathway (data not shown).
In additional to preventing nephritis, sirolimus also had striking effects on the anti-DNA antibody titres in mice with lupus [20–22, 30], so we addressed the connectivity of genes linked to any form of lupus with the mTOR pathway. About 50% of the lupus genes curated as lupus disease genes from human and rodent species in Ingenuity and MetaCore can be linked to the rapalog-mTOR pathway. The connectivity would, no doubt have been higher with the use of automatically extracted relationships from the biomedical literature. However, the algorithms used in automatic extractions cannot approximate human reasoning and return a mixture of true and false positives. Therefore, we relied exclusively on manually curated databases of protein findings and our results should be viewed as a lower estimate of connectivity.
To assess the significance of the association between human lupus genes and the mTOR pathway, we built an mTOR pathway interactome (ie, a network consisting of proteins that interact with the mTOR pathway) using IPA. We then queried which human disease networks (manually curated) in Metacore were best represented in the mTOR pathway interactome. Of the 87 human disease networks represented in Metacore, human lupus was identified as being highly significant, with only two cancers showing more significant associations. Additional cancer and non-cancer diseases were also identified through this process, including Alzheimer's disease and other autoimmune diseases such as multiple sclerosis and arthritis. Indeed recent work has uncovered a strong link between the mTOR pathway, Treg function and autoimmunity. Rapamycin was shown to inhibit AKT-mediated repression of FOXP3 (and AKT is upstream of mTOR). FOXP3 is a critical player in Treg cell differentiation and maintenance and deficiency of FOXP3 in both humans and mice is associated with multi-organ autoimmunity and lymphoproliferative disorders [37–40].
Having investigated the human disease-mTOR pathway connectivity, we then widened our analysis by exploring the validity of the claim of connectivity by searching the literature for data showing the effects on rapalogs on these human diseases. By conducting these analyses independently of Metacore, we confirmed the relationship between the mTOR pathway and some human diseases, such as multiple sclerosis [41, 42], diabetes [43–45], arthritis  and some cancers [47, 48]. A search of the clinical trial database  reports ongoing clinical studies with rapalogs in a number of these diseases, and the analyses we present here support such studies. Indeed early clinical results on the effects of sirolimus treatment of lupus patients show promise. Nine SLE patients that had been treated unsuccessfully with other immunosuppressive drugs had significantly improved disease scores after sirolimus treatment (BILAG p = 0.0218, SELDAI p = 0.00002) , and another clinical study is in progress .
Our analyses indicate that the coverage of protein-protein interactions in curated databases such as Ingenuity and MetaCore is comparable with up-to-date text-mined content derived using MedScan, a data mining/natural language processing tool (Ariadne Inc., Rockville, MD). For example, Ingenuity has 80 and MetaCore has 65 proteins/complex/groups that interact with the mTOR protein and MedScan identifies 115 proteins that interact with the mTOR protein. This level of overlap indicates a comprehensive coverage in the databases used for these analyses.
Given our results and the results of others showing that inhibition of the mTOR pathway prevents progression of lupus nephritis in various mice models, we reasoned that perturbations of the mTOR pathway can lead to the phenotype of lupus nephritis. We also assessed the involvement of the mTOR pathway in human lupus by building the mTOR pathway interactome (genes connected to the mTOR pathway) and using bioinformatic algorithms to determine the significance of the overlap between the mTOR interactome and the published findings on genes involved in human lupus. We found a highly significant overlap. We suggest a similar approach of assessing significance of overlap between genes linked to human diseases and networks controlling animal model perturbations can be useful in understanding the relevance of animal models and the exploration of new indications for established therapies.
analysis of variance
false discovery rate
- H & E:
haematoxylin & eosin
Ingenuity Pathway Analysis
mammalian target of rapamycin
systemic lupus erythromatous.
We are grateful for Andrew Hill's help in suggesting appropriate reference datasets for assessing over-representation of glomerular genes in the nephritis gene set, and for Lawrence Mason's processing and staining of tissues used for the histological assessments.
- Schroder NW, Maurer M: The role of innate immunity in asthma: leads and lessons from mouse models. Allergy. 2007, 62: 579-590. 10.1111/j.1398-9995.2007.01386.x.View ArticlePubMedGoogle Scholar
- Sena E, Worp van der HB, Howells D, Macleod M: How can we improve the pre-clinical development of drugs for stroke?. Trends Neurosci. 2007, 30: 433-439. 10.1016/j.tins.2007.06.009.View ArticlePubMedGoogle Scholar
- Sena E, Wheble P, Sandercock P, Macleod M: Systematic review and meta-analysis of the efficacy of tirilazad in experimental stroke. Stroke. 2007, 38: 388-394. 10.1161/01.STR.0000254462.75851.22.View ArticlePubMedGoogle Scholar
- Zosky GR, Larcombe AN, White OJ, Burchell JT, Janosi TZ, Hantos Z, Holt PG, Sly PD, Turner DJ: Ovalbumin-sensitized mice are good models for airway hyperresponsiveness but not acute physiological responses to allergen inhalation. Clin Exp Allergy. 2008, 38: 829-838. 10.1111/j.1365-2222.2007.02884.x.View ArticlePubMedGoogle Scholar
- Zosky GR, Sly PD: Animal models of asthma. Clin Exp Allergy. 2007, 37: 973-988. 10.1111/j.1365-2222.2007.02740.x.View ArticlePubMedGoogle Scholar
- Bertsias G, Ioannidis JP, Boletis J, Bombardieri S, Cervera R, Dostal C, Font J, Gilboe IM, Houssiau F, Huizinga T, Isenberg D, Kallenberg CG, Khamashta M, Piette JC, Schneider M, Smolen J, Sturfelt G, Tincani A, van Vollenhoven R, Gordon C, Boumpas DT, Task Force of the EULAR Standing Committee for International Clinical Studies Including Therapeutics: EULAR recommendations for the management of systemic lupus erythematosus. Report of a Task Force of the EULAR Standing Committee for International Clinical Studies Including Therapeutics. Ann Rheum Dis. 2008, 67: 195-205. 10.1136/ard.2007.070367.View ArticlePubMedGoogle Scholar
- Rahman A, Isenberg DA: Systemic lupus erythematosus. N Engl J Med. 2008, 358: 929-939. 10.1056/NEJMra071297.View ArticlePubMedGoogle Scholar
- Mason L, Berden J: Pathogenic factors for the development of lupus nephritis. Lupus. 2008, 17: 251-255. 10.1177/0961203307088250.View ArticlePubMedGoogle Scholar
- Boumpas DT: Lupus nephritis. Nephrol Dial Transplant. 2001, 16 (Suppl 6): 55-57.View ArticlePubMedGoogle Scholar
- Aringer M, Smolen JS: Cytokine expression in lupus kidneys. Lupus. 2005, 14: 13-18. 10.1191/0961203305lu2053oa.View ArticlePubMedGoogle Scholar
- Baudino L, Azeredo da Silveira S, Nakata M, Izui S: Molecular and cellular basis for pathogenicity of autoantibodies: lessons from murine monoclonal autoantibodies. Springer Semin Immunopathol. 2006, 28: 175-184. 10.1007/s00281-006-0037-0.View ArticlePubMedGoogle Scholar
- Haas C, Ryffel B, Le Hir M: IFN-gamma is essential for the development of autoimmune glomerulonephritis in MRL/Ipr mice. J Immunol. 1997, 158: 5484-5491.PubMedGoogle Scholar
- Miyake K, Nakashima H, Akahoshi M, Inoue Y, Nagano S, Tanaka Y, Masutani K, Hirakata H, Gondo H, Otsuka T, Harada M: Genetically determined interferon-gamma production influences the histological phenotype of lupus nephritis. Rheumatology (Oxford). 2002, 41: 518-524. 10.1093/rheumatology/41.5.518.View ArticleGoogle Scholar
- Santiago-Raber ML, Laporte C, Reininger L, Izui S: Genetic basis of murine lupus. Autoimmun Rev. 2004, 3: 33-39. 10.1016/S1568-9972(03)00062-4.View ArticlePubMedGoogle Scholar
- Foster MH: Relevance of systemic lupus erythematosus nephritis animal models to human disease. Semin Nephrol. 1999, 19: 12-24.PubMedGoogle Scholar
- Reilly CM, Gilkeson GS: Use of genetic knockouts to modulate disease expression in a murine model of lupus, MRL/lpr mice. Immunol Res. 2002, 25: 143-153. 10.1385/IR:25:2:143.View ArticlePubMedGoogle Scholar
- Liu J, Karypis G, Hippen KL, Vegoe AL, Ruiz P, Gilkeson GS, Behrens TW: Genomic view of systemic autoimmunity in MRLlpr mice. Genes Immun. 2006, 7: 156-168. 10.1038/sj.gene.6364286.View ArticlePubMedGoogle Scholar
- Shlomchik MJ, Craft JE, Mamula MJ: From T to B and back again: positive feedback in systemic autoimmune disease. Nat Rev Immunol. 2001, 1: 147-153. 10.1038/35100573.View ArticlePubMedGoogle Scholar
- Hutloff A, Buchner K, Reiter K, Baelde HJ, Odendahl M, Jacobi A, Dorner T, Kroczek RA: Involvement of inducible costimulator in the exaggerated memory B cell and plasma cell generation in systemic lupus erythematosus. Arthritis Rheum. 2004, 50: 3211-3220. 10.1002/art.20519.View ArticlePubMedGoogle Scholar
- Ramos-Barron A, Pinera-Haces C, Gomez-Alamillo C, Santiuste-Torcida I, Ruiz JC, Buelta-Carrillo L, Merino R, de Francisco AL, Arias M: Prevention of murine lupus disease in (NZBxNZW)F1 mice by sirolimus treatment. Lupus. 2007, 16: 775-781. 10.1177/0961203307081401.View ArticlePubMedGoogle Scholar
- Warner LM, Adams LM, Sehgal SN: Rapamycin prolongs survival and arrests pathophysiologic changes in murine systemic lupus erythematosus. Arthritis Rheum. 1994, 37: 289-297. 10.1002/art.1780370219.View ArticlePubMedGoogle Scholar
- Lui S, Yung S, Tsang R, Zhang F, Chan K, Tam S, Chan T: Rapamycin prevents the development of nephritis in lupus-prone NZB/W F1 mice. Lupus. 2008, 17: 305-313. 10.1177/0961203307088289.View ArticlePubMedGoogle Scholar
- Chiang GG, Abraham RT: Targeting the mTOR signaling network in cancer. Trends Mol Med. 2007, 13: 433-442. 10.1016/j.molmed.2007.08.001.View ArticlePubMedGoogle Scholar
- Brunet M, Campistol JM, Diekmann F, Guillen D, Millan O: T-cell function monitoring in stable renal transplant patients treated with sirolimus monotherapy. Mol Diagn Ther. 2007, 11: 247-256.View ArticlePubMedGoogle Scholar
- Sarbassov DD, Ali SM, Sabatini DM: Growing roles for the mTOR pathway. Curr Opin Cell Biol. 2005, 17: 596-603. 10.1016/j.ceb.2005.09.009.View ArticlePubMedGoogle Scholar
- MacDonald AS: A worldwide, phase III, randomized, controlled, safety and efficacy study of a sirolimus/cyclosporine regimen for prevention of acute rejection in recipients of primary mismatched renal allografts. Transplantation. 2001, 71: 271-280. 10.1097/00007890-200101270-00019.View ArticlePubMedGoogle Scholar
- GEO – Gene Expression Omnibus. [http://www.ncbi.nlm.nih.gov/geo]
- Reiner A, Yekutieli D, Benjamini Y: Identifying differentially expressed genes using false discovery rate controlling procedures. Bioinformatics. 2003, 19: 368-375. 10.1093/bioinformatics/btf877.View ArticlePubMedGoogle Scholar
- Tian L, Greenberg SA, Kong SW, Altschuler J, Kohane IS, Park PJ: Discovering statistically significant pathways in expression profiling studies. Proc Natl Acad Sci USA. 2005, 102: 13544-13549. 10.1073/pnas.0506577102.PubMed CentralView ArticlePubMedGoogle Scholar
- Alperovich G, Rama I, Lloberas N, Franquesa M, Poveda R, Gomà M, Herrero-Fresneda I, Cruzado JM, Bolaños N, Carrera M, Grinyó JM, Torras J: New immunosuppresor strategies in the treatment of murine lupus nephritis. Lupus. 2007, 16: 18-24. 10.1177/0961203306073136.View ArticlePubMedGoogle Scholar
- He L, Sun Y, Takemoto M, Norlin J, Tryggvason K, Samuelsson T, Betsholtz C: The glomerular transcriptome and a predicted protein-protein interaction network. J Am Soc Nephrol. 2008, 19: 260-268. 10.1681/ASN.2007050588.PubMed CentralView ArticlePubMedGoogle Scholar
- Peterson KS, Huang JF, Zhu J, D'Agati V, Liu X, Miller N, Erlander MG, Jackson MR, Winchester RJ: Characterization of heterogeneity in the molecular pathogenesis of lupus nephritis from transcriptional profiles of laser-captured glomeruli. J Clin Invest. 2004, 113: 1722-1733.PubMed CentralView ArticlePubMedGoogle Scholar
- Raught B, Gingras AC, Sonenberg N: The target of rapamycin (TOR) proteins. Proc Natl Acad Sci USA. 2001, 98: 7037-7044. 10.1073/pnas.121145898.PubMed CentralView ArticlePubMedGoogle Scholar
- Rotig A, Munnich A: Genetic features of mitochondrial respiratory chain disorders. J Am Soc Nephrol. 2003, 14: 2995-3007. 10.1097/01.ASN.0000095481.24091.C9.View ArticlePubMedGoogle Scholar
- Berger SP, Daha MR: Complement in glomerular injury. Semin Immunopathol. 2007, 29: 375-384. 10.1007/s00281-007-0090-3.PubMed CentralView ArticlePubMedGoogle Scholar
- Berger SP, Daha MR: Complement in glomerular nephritis. Semin Immunopathol. 2007, 29: 375-384. 10.1007/s00281-007-0090-3.PubMed CentralView ArticlePubMedGoogle Scholar
- Haxhinasto S, Mathis D, Benoist C: The AKT-mTOR axis regulates de novo differentiation of CD4+Foxp3+ cells. J Exp Med. 2008, 205: 565-574. 10.1084/jem.20071477.PubMed CentralView ArticlePubMedGoogle Scholar
- Qu Y, Zhang B, Zhao L, Liu G, Ma H, Rao E, Zeng C, Zhao Y: The effect of immunosuppressive drug rapamycin on regulatory CD4+CD25+Foxp3+T cells in mice. Transpl Immunol. 2007, 17: 153-161. 10.1016/j.trim.2007.01.002.View ArticlePubMedGoogle Scholar
- Strauss L, Whiteside TL, Knights A, Bergmann C, Knuth A, Zippelius A: Selective survival of naturally occurring human CD4+CD25+Foxp3+ regulatory T cells cultured with rapamycin. J Immunol. 2007, 178: 320-329.View ArticlePubMedGoogle Scholar
- Keever-Taylor CA, Browning MB, Johnson BD, Truitt RL, Bredeson CN, Behn B, Tsao A: Rapamycin enriches for CD4(+) CD25(+) CD27(+) Foxp3(+) regulatory T cells in ex vivo-expanded CD25-enriched products from healthy donors and patients with multiple sclerosis. Cytotherapy. 2007, 9: 144-157. 10.1080/14653240601145223.View ArticlePubMedGoogle Scholar
- Kappos L, Barkhof F, Desmet A: The effect of oral temsirolimus on new magnetic resonance imaging scan lesions, brain atrophy, and the number of relapses in multiple sclerosis; results from a randomized, controlled clinical trial. J Neurol. 2005, 252 ((Supplement 2);S46 (abstract)): II/46-Google Scholar
- Neuhaus O, Kieseier BC, Hartung HP: Immunosuppressive agents in multiple sclerosis. Neurotherapeutics. 2007, 4: 654-660. 10.1016/j.nurt.2007.08.003.View ArticlePubMedGoogle Scholar
- Dowling RJ, Zakikhani M, Fantus IG, Pollak M, Sonenberg N: Metformin inhibits mammalian target of rapamycin-dependent translation initiation in breast cancer cells. Cancer Res. 2007, 67: 10804-10812. 10.1158/0008-5472.CAN-07-2310.View ArticlePubMedGoogle Scholar
- Penfornis A, Kury-Paulin S: Immunosuppressive drug-induced diabetes. Diabetes Metab. 2006, 32: 539-546.View ArticlePubMedGoogle Scholar
- Battaglia M, Stabilini A, Migliavacca B, Horejs-Hoeck J, Kaupper T, Roncarolo MG: Rapamycin promotes expansion of functional CD4+CD25+FOXP3+ regulatory T cells of both healthy subjects and type 1 diabetic patients. J Immunol. 2006, 177: 8338-8347.View ArticlePubMedGoogle Scholar
- Foroncewicz B, Mucha K, Paczek L, Chmura A, Rowinski W: Efficacy of rapamycin in patient with juvenile rheumatoid arthritis. Transpl Int. 2005, 18: 366-368. 10.1111/j.1432-2277.2004.00070.x.View ArticlePubMedGoogle Scholar
- Tsang CK, Qi H, Liu LF, Zheng XF: Targeting mammalian target of rapamycin (mTOR) for health and diseases. Drug Discov Today. 2007, 12: 112-124. 10.1016/j.drudis.2006.12.008.View ArticlePubMedGoogle Scholar
- Faivre S, Kroemer G, Raymond E: Current development of mTOR inhibitors as anticancer agents. Nat Rev Drug Discov. 2006, 5: 671-688. 10.1038/nrd2062.View ArticlePubMedGoogle Scholar
- ClinicalTrial.gov, a service of the National Insitutes of Health. [http://clinicaltrials.gov]
- Fernandez D, Bonilla E, Mirza N, Niland B, Perl A: Rapamycin reduces disease activity and normalizes T cell activation-induced calcium fluxing in patients with systemic lupus erythematosus. Arthritis Rheum. 2006, 54: 2983-2988. 10.1002/art.22085.PubMed CentralView ArticlePubMedGoogle Scholar
- Sirolimus therapy for idopathic and lupus membranous nephropathy. [http://clinicaltrials.gov/ct2/show/NCT00050713]
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.