Altered gene expression in asymptomatic SHIV-infected rhesus macaques (Macacca mulatta)
- Erica E Carroll†1,
- Rasha Hammamieh†1,
- Nabarun Chakraborty1,
- Aaron T Phillips1,
- Stacy-Ann M Miller1 and
- Marti Jett1Email author
© Carroll et al; licensee BioMed Central Ltd. 2006
Received: 06 July 2006
Accepted: 06 September 2006
Published: 06 September 2006
Simian-Human immunodeficiency virus is a chimeric virus which, in rhesus macaques (Macacca mulatta) closely imitates immunodeficiency virus infection in human (HIV). A relatively new way to study pathogenesis of viral infection is to study alterations in host gene expression induced by the virus. SHIV infection with certain strains does not result in clinical signs. We hypothesized that alterations in gene expression relating to the immune system would be present in SHIV-infected animals despite the lack of clinical signs. Splenic tissue from four adult male Indian-origin Rhesus monkeys serologically positive for non-pathogenic SHIV 89.6 was processed by cDNA microarray analysis. Results were compared with the corresponding outcome using splenic tissues from four unexposed adult male Rhesus monkeys. Subsequent gene analysis confirmed statistically significant variations between control and infected samples. Interestingly, SHIV-infected monkeys exhibited altered expression in genes related to apoptosis, signal transduction, T and B lymphocyte activation and importantly, to immune regulation. Although infected animals appeared asymptomatic, our study demonstrated that SHIV-infected monkeys cannot reliably be used in studies of other infectious agents as their baseline gene expression differs from that of normal Rhesus monkeys. The gene expression differences in SHIV-infected animals relative to uninfected animals offer additional clues to the pathogenesis of altered immune function in response to secondary infection.
Simian immunodeficiency virus (SIV) infection of rhesus macaques exhibits many similarities to human immunodeficiency viral (HIV) infection of humans. Most pathogenesis and vaccine studies for HIV-1 have been undertaken in either SIV-macaque or a chimeric simian-human immunodeficiency (SHIV)-macaque model . SHIV strains have the viral envelope of HIV but the gag/pol genes of SIV. Pathogenesis is similar with respect to macrophage and T lymphocyte cell tropism, histopathologic changes, CD4-cell depletion and clinical signs of auto-immune deficiency syndrome (AIDS) in virulent strains. HIV and SIV additionally cause cognitive and motor impairments in infected patients and monkeys, respectively . Host factors may play a role in degree of pathogenesis between varying SHIV constructs, as one study reported observing similar viral loads in rhesus monkeys infected with pathogenic and non-pathogenic SHIV constructs .
Gene expression studies have grown increasingly popular as a tool to mine large amounts of data from treated and control populations. Such data can be used to examine host factors involved in SHIV, and thereby HIV, pathogenesis. To our knowledge, microarray data from SHIV-infected Rhesus macaques have not yet been examined for genes affecting immune response and inflammation. Gene expression data have the potential to greatly expand the understanding of SHIV-host interaction beyond the limited number of cell types or cytokines generally examined.
In animals free of clinical signs of SHIV, altered baseline gene expression data may give clues to the pathogenesis of altered immune response to secondary infections. Studies involving HIV-infected humans demonstrated suppression of IL-2 in response to select antigens and increase in TNF-α even prior to the onset of CD4+ T-cell depletion [3, 4]. Gene expression data collected in this study from SHIV 89.6-infected monkeys demonstrate that these animals are not genetically 'normal' and cannot ethically be used for studies involving other infectious agents, if at all, without an explicit caveat listing their SHIV status. Comparison of gene expression patterns collected from SHIV-infected and uninfected animals to that of the matched animals exposed to select bacterial and viral agents would provide a more complete understanding of SHIV effect on immune response to particular infectious agents. Extrapolation to the HIV-patient response to secondary agents may then be attempted. Gene expression data may also provide clues to pathogenesis of cognitive and related ailments arising with HIV infection.
An overview of the Rhesus macaques used in SHIV gene expression study
Herpes B Status
SHIV 89.6 status
Clinical pathology of SHIV-positive rhesus macaques
Abnormal findings in complete blood count and serum chemistry analysis.
Sodium 144 mg/dl (reference range 147–158)
Chloride 108 mg/dl (range 110–120)
Lymphocytes 65.4% (reference range 14–64%)
Sodium 146 mg/dl (range 147–158)
Carbon dioxide 29 mmol/L (range 19–29)
Total protein 6.4 g/dl (range 6.7–8.0)
ALT 113 U/L (range 20–91)
LDH 538 U/L (range 638–3012)
Sodium 145 mg/dl (range 147–158)
Chloride 110 mg/dl (range 110–120)
AST 29 U/L (range 29–64)
JGH (Herpes B+)
Sodium 147 mg/dl (range 147–158)
Chloride 109 mg/dl (range 110–120)
Carbon dioxide 30 mmol/L (range 19–29)
Triglycerides 18 mg/dl (range 35–137)
Total protein 6.6 g/dl (range 6.7–8.0)
AST 26 U/L (range 29–64)
The sequences of the primers used in the present project
Gene Bank ID
5'-TCC AAG TCC TGT GTC TCC TC-3'
5'-TCT TGA GTT CAG CCC TCT CC-3'
5'-TGG CTG TTG GAG TTG CTT G-3'
5'-TGT GTG TCT TGC TCG TCT TG-3'
chloride channel. calcium activated, fam
5'-CAA CCA AGA AGC ACC AA CC-3'
5'-CAT CCA GCA CTA AAC AGA CCA C-3'
postmeiotic segregation increased 2-like
5'-GTT TCA GGC AAT GGA TGT GG-3'
5'-CAT GGC AGG TAG AAA TGG TG-3'
collagen, type XV, alpha 1
5'-CCA CCT ACC GAG CAT TCT TAT C-3'
5'-CAA TAC GTC TCG ACC ATC AAA G-3'
interleukin 2 receptor, alpha
5'-CTG AGA GCA TCT GCA AAA TGA C-3'
5'-GGC CAC TGC TAC TTG GTA CTC T-3'
programmed cell death 4
5'-CCG GTG ATG AAG AAA ATG CT-3'
5'-TGG TTG GCA CAG TTA ATC CA-3'
adenosine A2a receptor
5'-TCA ACA GCA ACC TGC AGA AC-3'
5'-ATG GCA ATG TAG CGG TCA AT-3'
RNA binding motif protein 9
5'-AAC TCC TGA CTC AAT GGT TC-3'
5'-CAT TTT GTG TGC TGG GTG AG-3'
mitogen-activated protein kinase kinase
5'-ACC AGG CAG AAA TCA ACG AC-3'
5'-GAT GAA CGT CCC AAA GCA CT-3'
collagen, tykpe VII, alpha 1 (epidermolysin)
5'-AGC CCA GAT GTT TCC ACT CA-3'
5'-ACA AGA GGC AAT CCT TGG AGA-3'
Micro-array analysis of SHIV-infected versus uninfected
The list of some of the genes of interest.
Cellular defense immunity:
golgi autoantigen, subfamily a2
major histocompatibility complex, class II, DR beta 1
leukocyte immunoglobulin-like receptor, subfamily B, member 3
succinate dehydrogenase complex
3-hydroxyisobutyryl-Coenzyme A hydrolase
hepatic transcription factor 2
transition protein 2
aldehyde dehydrogenase 5 family
peter pan homolog
phorbol-12-myristate-13-acetate-induced protein 1
apoptosis inhibitor 5
nucleosome assembly protein 1-like 1
proline rich 2
histidine-rich calcium-binding protein
GTP-binding protein, member RAS oncogene family
alpha tubulin 2
15 kDa selenoprotein
collagen, type VI, alpha 2
postmeiotic segregation increased 2-like 5
cyclin-dependent kinase 2-interacting protein
nuclear fragile X mental retardation protein interacting protein 1
IQ motif containing GTPase activating protein
T-cell receptor active alpha-chain
calcium activated chloride channel
inositol 1,4,5-trisphosphate 3-kinase C
Confirmation of gene expression changes by Real-Time PCR analysis
Simian immunodeficiency virus (SIV), previously referred to as simian T-cell lymphotropic virus type III (STLV-III), induces an AIDS-like disease in its natural host, rhesus macaques. HIV and SIV, members of the lentivirus subfamily of retroviruses, not only resemble each other by their antigenicity, but also bear remarkable similarity in their biological properties, such as cytopathic effect and tropism for CD4-bearing cells. These criteria render the chimeric SHIV the best animal model currently available for HIV study.
In this study, we examined gene expression in SHIV- infected male rhesus macaques of Indian origin using a genomic perspective and compared the results to uninfected age, gender and Herpes B-status-matched controls. Although infected animals were without clinical signs related to SHIV infection, a significant number of genes exhibited significantly altered expression concurrent with SIV infection.
Ontological research revealed that several genes, namely FOS-like antigen 1 (FOSL1, ID: H96643), golgi autoantigen (GOLGA2, ID: AA424786), major histocompatibility complex (MHC), class II, DR beta 1 (HLA-DRB1, ID: AA664195) and leukocyte immunoglobulin-like receptor (LILRB3, ID: AI815229) are associated with human immune defense. LILRB3 is a leukocyte inhibitory receptor which, upon binding to MHC Class I molecules, transmits inhibitory signals to the nucleus. HLA-DRB1, down regulated by SIV infection, is a cell-surface-associated immunoregulatory protein. Interestingly, this human leukocyte antigen (HLA)-associated gene has been correlated with non-responsiveness to recombinant hepatitis B virus (HBV) vaccine but does not alter susceptibility to viral persistence . Another MHC protein binding unit, T cell receptor alpha locus (TRAC, ID: AA427491) is ontologically related to signal transduction.
Gene ontology investigation classified a significant subset of the genome of interest as a regulator of cell growth and apoptosis. SIV infection results in down-regulation of apoptosis inhibitor 5 (API5, ID: AI972925) and up-regulation of pro-apoptotic protein phorbol-12-myristate-13-acetate-induced protein 1 (NOXA, ID: AA458838) . These alterations in gene expression might instigate opportunistic infections by inducing apoptosis among T-helper lymphocytes. Likewise, SIV infection alters several metabolism and cell growth regulating factors. For example, SIV-infected genome contains upregulated aldehyde dehydrogenase 5 family member A1 (ALDH5A1, ID: H06676); and concurrent down regulated succinate dehydrogenase complex, subunit D (SDHD, ID: AA035384) and nephropathic cystinosis (CTNS, ID: W94331). Reports suggest that overexpressed ALDH5A1 changes the concentration of gamma-aminobutyric acid (GABA) and glutamate, commencing henceforth excitotoxic damage, a well-established clinical marker of HIV activity . Underexpressed SDHD and CTNS are associated with immunodeficiency through curbed monocyte and CD4+ T cell -induced immunoregulation , respectively.
Several entries of the present genome are functionally related to cellular and molecular transportation and binding. Interestingly, five actin-binding genes appeared in the list; namely: anillin (ANLN, ID: R16712), destrin (DSTN, ID: AA424824), utrophin (UTRN, ID: AA676840), cyclin-dependent kinase 2-interacting protein (CINP, ID: AI364103) and IQ motif containing GTPase activating protein 2 (IQGAP2, ID: W32272). Actin, the ubiquitously present cellular protein, has been reported to guide the direct cell-to-cell HIV-1 propagation by making of a stable adhesive junction at the target-effector cell interface . Table 4 displays the down regulation of another molecular binding protein, 15 kDa selenoprotein (SEP15, ID: AA521350). Reduced level of selenoprotein in cells is a known marker of in vitro infection of SHIV . Our data also supports the fact that immunodeficiency is correlated with altered calcium ion binding (UTRN, ID: AA676840; CDH6, ID: AA421819, CASQ2, ID: AA055163) and also is influenced by calcium- activated chloride channels (CLCA2, ID: AI675394) of host cells. Those are well established pathoregulating markers of activ HIV-1 negative factor (Nef) [13–15].
In summary, in this small sample of SHIV-infected Rhesus macaques, expression was consistently altered in specific groups of genes which regulate a broad range of biochemical functions. A few important members of the genome of interest are discussed here. The present study, along with correlating some genes with SHIV and HIV model, identifies several novel genes as potential therapeutic markers for immune deficiency studies. Furthermore, results of this study suggest that SHIV infection of rhesus macaques may influence immune response to a second agent, even if baseline levels of clinical measurements appear normal. This study substantiates and validates the concern that an infected (i.e., antibody-producing) but asymptomatic animal should not be used in any other study involving infectious agents unless the pattern of gene expression to that agent is compared to normal animals' pattern, one agent at a time.
Note: microarray data have been submitted to the Gene Expression Omnibus (GEO) and can be searched using the Platform ID: GPL3395.
Materials and methods
Animals and virus
Four adult (7–8 years old) male Rhesus macaques (one Herpes B-positive and three Herpes B-negative) that were previously exposed to SHIV 89.6 strain (Animal identifications: FFG, JGH, PHB and TTH) were euthanized due to being declared 'excess' and no longer usable due to their serologically positive SHIV status. Splenic tissue was collected from each animal upon euthanasia and immersed in RNA Later® for 30–60 minutes before freezing at -80C.
SHIV 89.6, like all SHIV strains, has the env gene from the HIV-1 strain. All four animals had been challenged with 1.0 ml intravenous SHIV 89.6, a non-pathogenic strain, and became seropositive. Previous studies by the same researchers showed that seropositive animals were PCR positive as well (WRAIR Protocol TO03-98). All animals remained free of clinical signs. Complete blood counts and serum chemistry profiles were performed on the SHIV-positive animals and were within or very close to normal limits. The negative control animals were Indian-origin adult male Herpes B-negative Rhesus macaques. Splenic tissues were kindly provided by Scripps Institute, the National Institute of Health, and the Oregon National Primate Research Center. Tissue from a SHIV-negative animal (DB-87, provided by the Tulane National Regional Primate Research Center) was Herpes B-positive to control for the Herpes B-positive status of one SHIV-infected animal Table 1 represents an overview of the Rhesus macaques used in this study. Table 2 shows the clinical results of the SHIV-positive rhesus macaques
Splenic tissue samples stored in RNALater® (Ambion, TX) at -80C were thawed in 1.5 mL tubes on ice. Tissue was submerged in Trizol ™ (Invitrogen, CA) solution and RNA isolation was carried out paccording to the Trizol ™ Reagent manufacturer's recommended instructions. RNA was ethanol-precipitated, air-dried and re-suspended in 20 ul/sample of nuclease-free water. RNA quantity was measured via spectrophotometry followed by analysis with a Bioanalyzer 2100 (Agilent Technologies, CA)
Custom made cDNA microarray SlidePreparation and hybridization
The gene library for the present project was commercially obtained from Research Genetics (Invitrogen, CA), containing 7489 genes, including 7019 known genes, 249 unknown genes and 110 expressed sequence tagged genes (ESTs). Superamine coated Telechem slides (Telechem Inc., OR) were used for printing the cDNA clones using 12 × 4 pin format, on a Virtek chip writer professional microarrayer in KemTek, Inc, MD. The printed slides underwent UV cross-linking, followed by post-processed by succinic anhydride treatment. The Micromax™ Tyramide Signal Amplification (TSA)™ Labeling and Detection Kit (PerkinElmer, Inc., MA) was used as directed by the manufacturer to determine relative gene expression of the collected samples. Custom-made reference RNA was prepared by combining aliquots of RNA from 33 normal Rhesus tissues and was used on every slide as the array controller, to check overall sensitivity of array printing, and to monitor reverse transcription, labeling and hybridization efficiency. Sample hybridization was carried out at 55°C for sixteen hours. A laser detection system was used (GenePix 4000b, Axon Instruments, CA) to scan the finished slides. Intensity of the scanned images was digitalized through Genepix 4.0 software (Axon Inc., CA).
Data cleansing and statistical analysis was carried out using Genespring® 7.0 (Agilent Tech., CA). Local background was subtracted from individual spot intensity. Genes that failed this 'background check' in any of the eight given experiments were eliminated from further analysis. Each chip was next subjected to intra-chip normalization (LOWESS). The genes that varied most between control and treated sample sets were selected via t-test analysis. The p-value cutoff was set at 0.05. Four hundreds and thirty two genes were differentially expressed between SHIV -infected and control uninfected animals with p < 0.05.
The pattern of gene expression variability of the experimental set having reduced dimension was evaluated using principal component analysis (PCA) classifying SHIV positive and negative samples as the two variable classes .
Real Time PCR
The t-test result was corroborated through real time polymerized chain reaction (Real-time PCR). A web-based primer designing tool was used to design the primers for the selected genes . The specificity of each primer sequence was further confirmed by running a blast search. Reverse transcription and Real-time PCR reactions were carried out using reverse transcription kit (Invitrogen, CA) and Real-time PCR kit (Roche, IN), respectively. Each reaction with five technical duplicates was run in I-Cycler machine (Bio-Rad, CA). Each sample was also amplified against the house-keeping probe of the experiment: glyceraldehyde 3 phosphate dehydrogenase (GAPDH). The resultant cycle threshold data from each real-time-PCR 'run' was converted to fold-change using an established algorithm .
Quantitative and qualitative verification of the PCR product was accomplished by performing 1% agarose gel electrophoresis using SYBR Green I (Kemtek, Rockville, MD). Gel images were captured using PharosFX Molecular Imager system (Bio-Rad, CA) scanner and analyzed using Quantity One software (Bio-Rad, CA).
EEC wants to extend thanks to LTC Gary D. Coleman and LTC Keith E. Steele for giving her the time to devote to this project in the face of other equally pressing mission requirements.
- Nath BM, Schumann KE, Boyer JD: The chimpanzee and other non-human-primate models in HIV-1 vaccine research. Trends Microbiol 2000,8(9):426-31. 10.1016/S0966-842X(00)01816-3View ArticlePubMedGoogle Scholar
- Murray EA, Rausch DM, Lendvay J, Sharer LR, Eiden LE: Cognitive and motor impairments associated with SIV infection in rhesus monkeys. Science 1992,255(5049):1246-9.View ArticlePubMedGoogle Scholar
- Blackburn R, Clerici M, Mann D, Lucey DR, Goedert J, Golding B, Shearer GM, Golding H: Common sequence in HIV 1 GP41 and HLA class II beta chains can generate crossreactive autoantibodies with immunosuppressive potential early in the course of HIV 1 infection. Adv Exp Med Biol 1991, 303: 63-9.View ArticlePubMedGoogle Scholar
- Lane BR, Markovitz DM, Woodford NL, Rochford R, Strieter RM, Coffey MJ: TNF-alpha inhibits HIV-1 replication in peripheral blood monocytes and alveolar macrophages by inducing the production of RANTES and decreasing C-C chemokine receptor 5 (CCR5) expression. J Immunol 1999,163(7):3653-61.PubMedGoogle Scholar
- Hammamieh R, Chakraborty N, Das R, Jett M: Molecular impacts of antisense complementary to the liver fatty acid binding protein (FABP) mRNA in DU 145 prostate cancer cells in vitro. J Exp Ther Oncol 2004,4(3):195-202.PubMedGoogle Scholar
- Hammamieh R, Chakraborty N, Laing M, Liu Z, Mulligan J, Wang Y, Jett M: GeneCite: tool for high throughput literature and pathway mining. , in press.Google Scholar
- Wang C, Tang J, Song W, Lobashevsky E, Wilson CM, Kaslow RA: HLA and cytokine gene polymorphisms are independently associated with responses to hepatitis B vaccination. Hepatology 2004,39(4):978-88. 10.1002/hep.20142View ArticlePubMedGoogle Scholar
- Leal DB, Streher CA, Bertoncheli Cde M, Carli LF, Leal CA, da Silva JE, Morsch VM, Schetinger MR: HIV infection is associated with increased NTPDase activity that correlates with CD39-positive lymphocytes. Biochim Biophys Acta 2005,1746(2):129-34. 10.1016/j.bbamcr.2005.10.009View ArticlePubMedGoogle Scholar
- Koutsilieri E, Sopper S, Heinemann T, Scheller C, Lan J, Stahl-Hennig C, ter Meulen V, Riederer P, Gerlach M: Involvement of microglia in cerebrospinal fluid glutamate increase in SIV-infected rhesus monkeys (Macaca mulatta). AIDS Res Hum Retroviruses 1999,15(5):471-7. 10.1089/088922299311222View ArticlePubMedGoogle Scholar
- Flo RW, Naess A, Nilsen A, Harthug S, Solberg CO: A longitudinal study of phagocyte function in HIV-infected patients. Aids 1994,8(6):771-7.View ArticlePubMedGoogle Scholar
- Jolly C, Kashefi K, Hollinshead M, Sattentau QJ: HIV-1 cell to cell transfer across an Env-induced, actin-dependent synapse. J Exp Med 2004,199(2):283-93. 10.1084/jem.20030648PubMed CentralView ArticlePubMedGoogle Scholar
- Torrealba J: Selenium binding protein 1: passive or active role in disease? Am J Transplant 2005,5(10):2593. 10.1111/j.1600-6143.2005.01059.xView ArticleGoogle Scholar
- Shoeman RL, Kesselmier C, Mothes E, Honer B, Traub P: Non-viral cellular substrates for human immunodeficiency virus type 1 protease. FEBS Lett 1991,278(2):199-203. 10.1016/0014-5793(91)80116-KView ArticlePubMedGoogle Scholar
- Matsubara M, Jing T, Kawamura K, Shimojo N, Titani K, Hashimoto K, Hayashi N: Myristoyl moiety of HIV Nef is involved in regulation of the interaction with calmodulin in vivo. Protein Sci 2005,14(2):494-503. 10.1110/ps.04969605PubMed CentralView ArticlePubMedGoogle Scholar
- Liu QH, Williams DA, McManus C, Baribaud F, Doms RW, Schols D, De Clercq E, Kotlikoff MI, Collman RG, Freedman BD: HIV-1 gp120 and chemokines activate ion channels in primary macrophages through CCR5 and CXCR4 stimulation. Proc Natl Acad Sci USA 2000,97(9):4832-7. 10.1073/pnas.090521697PubMed CentralView ArticlePubMedGoogle Scholar
- Raw data[http://www.ncbi.nlm.nih.gov/geo/]
- PCR design[http://frodo.wi.mit.edu/cgi-bin/primer3]
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.