Skip to main content

Glyceraldehyde 3-phosphate dehydrogenase negatively regulates human immunodeficiency virus type 1 infection

Abstract

Background

Host proteins are incorporated inside human immunodeficiency virus type 1 (HIV-1) virions during assembly and can either positively or negatively regulate HIV-1 infection. Although the identification efficiency of host proteins is improved by mass spectrometry, how those host proteins affect HIV-1 replication has not yet been fully clarified.

Results

In this study, we show that virion-associated glyceraldehyde 3-phosphate dehydrogenase (GAPDH) does not allosterically inactivate HIV-1 reverse transcriptase (RT) but decreases the efficiency of reverse transcription reactions by decreasing the packaging efficiency of lysyl-tRNA synthetase (LysRS) and tRNALys3 into HIV-1 virions. Two-dimensional (2D) gel electrophoresis demonstrated that some isozymes of GAPDH with different isoelectric points were expressed in HIV-1-producing CEM/LAV-1 cells, and a proportion of GAPDH was selectively incorporated into the virions. Suppression of GAPDH expression by RNA interference in CEM/LAV-1 cells resulted in decreased GAPDH packaging inside the virions, and the GAPDH-packaging-defective virus maintained at least control levels of viral production but increased the infectivity. Quantitative analysis of reverse transcription products indicated that the levels of early cDNA products of the GAPDH-packaging-defective virus were higher than those of the control virus owing to the higher packaging efficiencies of LysRS and tRNALys3 into the virions rather than the GAPDH-dependent negative allosteric modulation for RT. Furthermore, immunoprecipitation assay using an anti-GAPDH antibody showed that GAPDH directly interacted with Pr55gag and p160gag-pol and the overexpression of LysRS in HIV-1-producing cells resulted in a decrease in the efficiency of GAPDH packaging in HIV particles. In contrast, the viruses produced from cells expressing a high level of GAPDH showed decreased infectivity in TZM-bl cells and reverse transcription efficiency in TZM-bl cells and peripheral blood mononuclear cells (PBMCs).

Conclusions

These findings indicate that GAPDH negatively regulates HIV-1 infection and provide insights into a novel function of GAPDH in the HIV-1 life cycle and a new host defense mechanism against HIV-1 infection.

Background

Because the HIV-1 genome only encodes a limited number of viral proteins, HIV-1 must take advantage of multiple functions of host proteins in order to successfully replicate. Several studies of purified HIV-1 virions have shown that in addition to proteins encoded by the virus, host proteins are found in the virions [1]. Some of these proteins may be taken into the virions simply because of their proximity to the viral assembly and budding sites, while other host proteins, such as cyclophilin A and LysRS, are included in HIV-1 particles as a result of their interaction with gag or gag pol proteins during assembly [27]. These host proteins play an important role in facilitating the process of gag protein folding and tRNALys3 packaging. Therefore, one way to elucidate the viral replication capacity gained by the packaging of host proteins is to directly analyze the host proteins inside the virions. A purified HIV-1LAV-1 preparation was analyzed by 2D gel electrophoresis and matrix-assisted laser desorption/ionization-time-of-flight mass spectrometry (MALDI-TOF MS). Proteome analysis demonstrates that GAPDH is inside the virions.

GAPDH is a prototype “moonlighting” protein that is involved in glycolysis, the carbon reduction cycle, the exportation of nuclear RNA, DNA repair, the bundling of microtubules, and apoptosis [814]. Furthermore, GAPDH also regulates viral replication. Interestingly, GAPDH phosphorylates the hepatitis B virus core protein [15] and binds with the cis-acting RNAs of several viruses, such as the 5′ UTR of the hepatitis A virus [16], the 3′ UTR sequences of the human parainfluenza virus [17], and the hepatitis C virus [18]. These findings indicated that GAPDH might also play a role in regulating HIV-1 replication.

In this study, we show that GAPDH is incorporated into the virions and the suppression of GAPDH packaging inside the virions enhances viral infection owing to the high reverse transcription efficiency. These results elucidate the viral replication capacity gained by GAPDH packaging and reveal a novel regulation step of HIV-1 infection.

Results

Multiple GAPDH isozymes were incorporated into virions

To examine the host proteins found in virions produced by an HIV-infected cell line, CEM/LAV-1 cells, a highly purified HIV-1LAV-1 preparation was analyzed by 2D gel electrophoresis, followed by MALDI-TOF MS/MS (Figure 1A). The MS/MS ion search predicted the identity of the 36-kDa proteins to be GAPDH. To confirm whether GAPDH is incorporated into virions, contaminating particles such as microvesicles and exosomes were removed according to previously described methods [19]. Western immunoblotting for CD45, GAPDH, and the capsid (CA) protein was carried out. In the starting sample, the CD45 signal was detected (Figure 1B left panel). On the other hand, no CD45 signal was detected in the lane of the CD45-depleted fraction (flow-through fraction), indicating that most of the CD45-positive contaminating particles such as microvesicles and exosomes in the starting material were removed (Figure 1B left panel). Therefore, as shown in Figure 1B right panel, western immunoblotting for GAPDH, and the CA protein was carried out. The data showed that there were evident GAPDH and CA signals in the CD45-depleted fraction, while there is hardly any GAPDH signal in the CD45-containing fraction, indicating that GAPDH from CEM/LAV-1 cells may exclusively be incorporated into the virions (Figure 1B right panel). Therefore, unless otherwise noted, HIV-1 preparations were systematically prepared without the CD45 depletion procedure in the following study. To gain further insight into the incorporation of GAPDH into virons, we analyzed GAPDH incorporated into virions and present in CEM/LAV-1 cells and the HIV-1-noninfected cell line CEM cells using 2D gel electrophoresis and western immunoblotting. We found that five isozymes of GAPDH (spot i (pI 7.61), spot ii (pI 7.63), spot iii (pI 8.03), spot iv (pI 8.50), and spot v (pI 9.02)) were inside the virions (Figure 1C). Surprisingly, regardless of HIV-1 infection, in the cellular pool, GAPDH was identified as six isozymes that contain an additional isozyme (spot vi (pI 8.68)) with five isozymes packaged inside the virions (Figure 1D). These results suggest that the GAPDH isozymes (pI 7.61, 7.63, 8.03, 8.50, and 9.02) might be selectively packaged inside the virions. Therefore, we focused primarily on the characterization of GAPDH present within the virions.

Figure 1
figure 1

Confirmation of GAPDH incorporation into HIV-1 virions. (A) SYPRO® Ruby-stained 2D gel image of proteins from highly purified HIV-1LAV-1 virions within the pI range of 6.85-9.30 and the molecular weight range of 30–46 kDa. (B) Western immunoblotting of GAPDH in CD45-depleted fraction and CD45-containing fraction. CD45 and CA levels were used as the markers of microvesicles and viral particles, respectively. Samples from the CD45-immunoaffinity-depleted fraction (CD45-depleted fraction) and anti-CD45 beads fraction (CD45-containing fraction) are identified above their respective lanes. The antibodies used are indicated on the left side of each blot. (C) Detection of virion-associated GAPDH isozymes. The upper panel shows the results of western immunoblotting using an anti-GAPDH antibody and the lower panel indicates the corresponding SYPRO® Ruby-stained gel. (D) Comparison of pI values of GAPDH isozymes among HIV-1LAV-1 (upper panel), CEM/LAV-1 (middle panel), and CEM cells (lower panel). GAPDH isozymes were detected by western immunoblotting using the anti-GAPDH antibody. Spot pI values are shown in the figure.

GAPDH-packaging-defective virus exhibits higher infectivity than control viruses

To determine whether a decrease in the packaging level of virion-associated GAPDH affects HIV-1 infectivity, we transfected CEM/LAV-1 cells with GAPDH siRNA or control siRNA, and observed no toxic effects on cell viability and proliferation (Figure 2A). As shown in Figures 2Band C, GAPDH siRNA specifically decreased GAPDH level not only in CEM/LAV-1 cells without altering the level of actin but also in the virions without altering the levels of gp120, HIV-1 RT, Pr55gag, and CA protein. Furthermore, p24 ELISA indicated that the GAPDH-packaging-defective virus did not show a defect in viral release (Figure 2D) and the incorporation of HIV genomic RNA into viral particles (Figure 2E). However, the GAPDH-packaging-defective virus showed an increased infectivity (Figure 2F, p< 0.01) compared with the control virus. To further clarify the relationship between the increased infectivity and the reverse transcription efficiency of the GAPDH-packaging-defective virus in TZM-bl cells, we carried out quantitative real-time PCR analysis to quantify the early strong-stop DNA of reverse transcription. The GAPDH-packaging-defective virus showed a significant increase in the level of the early strong-stop form of the viral cDNA product (Figure 2G, p<0.01). These results suggest that the enhanced infectivity of the GAPDH-packaging-defective virus is linked to a higher efficiency of the early reverse transcription process.

Figure 2
figure 2

Effects of GAPDH packaging defect inside virions on HIV-1 replication. (A) Effects of GAPDH siRNA treatment on CEM/LAV-1 cells. Viable cells (left) and dead cells (right) were assessed by trypan blue staining. (B) GAPDH siRNA knockdown efficiency is confirmed by western immunoblotting. Expression of GAPDH was analyzed in cell lysates from CEM/LAV-1 cells transfected with GAPDH or control siRNA 72 h after transfection. (C) Effects of GAPDH siRNA on incorporation of GAPDH and HIV-1 proteins inside virions. The incorporation of GAPDH, gp120, RT, Pr55gag, or CA was analyzed by western immunoblotting of lysates from viruses produced from CEM/LAV-1 cells transfected with GAPDH or control siRNA. The anti-GAPDH antibody and HIV-1-positive plasma were used for western immunoblotting. (D) Effects of GAPDH siRNA on virus release in CEM/LAV-1 cells transfected with GAPDH or control siRNA. The virus release in culture supernatant was directly determined by p24 ELISA. The mean values of at least three independent experiments are shown. (E) Effects of GAPDH siRNA on the incorporation of HIV genomic RNA into viral particles. The level of genomic RNA in the control virus (normalized to RT activities) was set as 100%. (F) Infectivity of GAPDH-packaging-defective virus. Infectivity was evaluated on the basis of the luciferase activity in lysates of TZM-bl cells. The value in the control experiment was set as 100%. The mean values of at least three independent experiments are shown. (G) Effect of defect in GAPDH packaging on reverse transcription in TZM-bl cells. The early strong-stop DNA products were determined by quantitative real-time PCR analysis as described in “Methods”. The significance of difference (Student’s t-test) is indicated as follows: *, p<0.01; n.s., not significant. The mean values of at least three independent experiments are shown. The error bars denote the standard deviation.

A suppressed packaging of GAPDH is accompanied by an enhanced packaging of LysRS and tRNALys3in the virions

To investigate why the GAPDH-packaging-defective virus showed a significant increase in the level of the early (R/U5) form of the viral cDNA product, in vitro HIV-1 RT activity was quantified in the presence of GAPDH. In RT assay, the ability of HIV-1 RT to synthesize DNA is utilized using hybrid poly(A)-oligo (dT)15 as a template and the primer, in the presence of recombinant GAPDH. The assays showed that GAPDH did not allosterically inhibit the activity of HIV-1 RT (Figure 3A). Gabor et al. [4] reported that a greater packaging of LysRS into virions is accompanied by increased tRNALys3 packaging, initiation of reverse transcription, and increased infectivity of the viral population. Therefore, we next investigated whether a decrease in virion-associated GAPDH level results in an enhanced packaging of LysRS and tRNALys3 inside the virions. Figures 3Band C show the effect of GAPDH siRNA treatment on the packaging of LysRS and tRNALys3 inside the virions. Western immunoblotting and quantitative real-time PCR analysis demonstrated that the reduction in virion-associated GAPDH level correlates with an about 1.5-fold increase in viral LysRS packaging level (Figure 3B) and an about 4-fold increase in viral tRNALys3 packaging level (Figure 3C, p< 0.01). Because Javanbakht et al. [20] and Mak et al. [21] further reported that LysRS packaging requires the viral Pr55gag and tRNALys3 packaging additionally requires the p160gag-pol precursor, we examined whether GAPDH could interact directly with Pr55gag and p160gag-pol. Immunoprecipitation assay using an anti-GAPDH antibody showed that GAPDH directly interacts with Pr55gag and p160gag-pol (Figure 3D). Furthermore, we concluded that GAPDH does not directly interact with LysRS because no sufficient signal was detected (data not shown). However, the overexpression of LysRS in HIV-1-producing cells results in an increase in LysRS packaging level in HIV particles but a decrease in GAPDH packaging level (Figure 3E). These results suggest that the incorporation of GAPDH inside virions may contribute to the suppression of initiation of reverse transcription owing to the suppression of the packaging of LysRS and tRNALys3 inside the virions.

Figure 3
figure 3

Packaging of LysRS and tRNALys3is regulated by GAPDH. (A) Effect of GAPDH on enzymatic activity of HIV-1 RT. RT activity assay was performed as described in “Methods”. The value in the control experiment was set as 100%. The activity in the presence of GAPDH (RT:GAPDH ratio= 1:10 or 1:100) is shown as the activity relative to that of the control. The mean values of at least three independent experiments are shown. Packaging of (B) LysRS and (C) tRNALys3 in GAPDH-packaging-defective virus. (B) LysRS and GAPDH were detected by western immunoblotting using anti-LysRS and anti-GAPDH antibodies in lysates from viral particles produced from CEM/LAV-1 cells transfected with GAPDH or control siRNA. (C) Incorporated tRNALys3 level was determined by reverse transcription quantitative PCR analysis as described in “Methods” and normalized by viral genomic RNA level. The amount of tRNALys3 in the control virus was set as 100%. The mean values of at least three independent experiments are shown. (D) Interaction of GAPDH with Pr55gag and p160gag-pol. GAPDH was immunoprecipitated from the clarified lysate from CEM/LAV-1 cells with the anti-GAPDH antibody. The precipitated proteins were analyzed by western immunoblotting using the indicated antibodies to the proteins shown (anti-p24 antibody for Pr55gag and anti-RT antibody for p160gag-pol), and were visualized by enhanced chemiluminescence analysis. The significance of difference (Student’s t-test) is indicated as follows: *, p<0.01; n.s., not significant. The error bars denote the standard deviation. (E) The packaging of GAPDH in the enhanced-LysRS-packaging virus. The packaging efficiencies of LysRS and GAPDH were analyzed in lysates from viral particles produced from HEK293 cells cotransfected with pNL-CH and either the LysRS expression or control vector.

An enhanced GAPDH packaging decreased viral infectivity

We next investigated whether the virus produced from HEK293 cells with high GAPDH expression levels showed decreased infectivity. We cotransfected HEK293 cells with the HIV-1 expression plasmid pNL-CH and either a GAPDH expression vector or a control vector. As shown in Figure 4A, the GAPDH expression vector increased the GAPDH level inside the virions. As expected, the enhanced-GAPDH-packaging virus decreased the viral infectivity in TZM-bl cells because the enhanced GAPDH packaging decreased tRNALys3 packaging (Figures 4Band C). Furthermore, we examined whether the enhanced packaging of GAPDH suppressed the viral replication in primary PBMCs. Interestingly, the enhanced-GAPDH-packaging virus showed a significant decrease in the level of the R/gag DNA product in PBMCs (Figure 4D). These results suggest that an enhanced GAPDH packaging decreased viral infectivity.

Figure 4
figure 4

Effects of enhanced GAPDH packaging inside virions on HIV-1 replication. (A) Enhanced GAPDH packaging inside virions. The packaging efficiency of GAPDH was analyzed in lysates from viral particles produced from HEK293 cells cotransfected with pNL-CH and either the GAPDH expression or control vector. (B) Infectivity of enhanced-GAPDH-packaging virus. TZM-bl cells were incubated with an equal amount (2 ng of p24 antigen) of the enhanced-GAPDH-packaging virus or the control virus. Infectivity was evaluated as described in “Methods”. The mean values of at least three independent experiments are shown. (C) Quantification of tRNALys3 in enhanced-GAPDH-packaging virus. Incorporated tRNALys3 level was determined by reverse transcription quantitative PCR analysis as described in “Methods” and normalized by viral genomic RNA level. The amount of tRNALys3 in the control virus was set as 100%. The mean values of at least three independent experiments are shown. (D) Enhanced GAPDH packaging inside virions decreased reverse transcription efficiency in PMBCs. PBMCs were incubated with an equal amount (20 ng of p24 antigen) of the enhanced-GAPDH-packaging virus or the control virus. Reverse transcription products (R/gag DNA products) were determined as described in “Methods”. The significance of difference (Student’s t-test) is indicated as follows: **, p<0.05. The mean values of at least three independent experiments are shown. The error bars denote the standard deviation.

Discussion

Current drugs against HIV-1 inhibit the function of viral enzymes, namely, RT, integrase, and protease. While these are especially proving useful in combination therapy, drug resistance could still be generated and appears to result from viral replication that allows the error-prone reverse transcription step. Therefore, future directions towards a more effective therapy for HIV-1 infection will rely on the development of novel therapeutic strategies rather than conventional strategies, which target only viral proteins. HIV-1 exploits multiple host proteins during infection, suggesting that the possibility of interrupting virus-host interactions may be an important pathway for the development of antiviral therapies. In concept, this type of antiviral agent developed through this pathway could minimize the generation of drug-resistant mutants because the virus must maintain the ability to interact with the relatively immutable host proteins. Thus, the interaction between the host and viral proteins may offer some potential applications for therapies against HIV-1.

Proteomic analysis of HIV-1 particles is a powerful tool to identify not only novel host proteins that are packaged during viral budding but also post-translational modifications of host and viral proteins in the virions. In our previous study, three isoforms of cyclophilin A associated with HIV-1 and the post-translational modifications such as N-acetylation of cyclophilin A and formylation of CA protein were identified [22, 23]. In this study, GAPDH isozymes with different isoelectric points were detected in the virions. Although Chertova et al. [24] and Ott et al. [25] also reported the GAPDH incorporation into HIV-1 virions produced from lymphocytes and macrophages, the virological significance of GAPDH has remained unclear. To gain insight into the role of GAPDH in viral replication, we prepared the GAPDH-packaging-defective virus by suppressing GAPDH expression by RNA interference in CEM/LAV-1 cells, and investigated the phenotypic and functional property changes of the virus. Here, we demonstrated that the suppression of GAPDH packaging increased the efficiencies of LysRS and tRNALys3 packaging into the virions, resulting in the increase in infectivity, and that GAPDH directly interacted with Pr55gag and p160gag-pol, which are required for the packaging of the LysRS and tRNALys3 complex. Javanbakht et al. [20, 26] reported that the interaction between the gag protein and LysRS is dependent on the last 54 amino acids of the CA C-terminal domain of the gag protein and amino acids 208–259 of LysRS and that the packaging of tRNALys3 requires interaction with LysRS. Furthermore, Khorchid et al. [27] reported that the interaction between p160gag-pol and tRNALys3 involves the thumb domain (TH) sequences in RT. From these findings, Kleiman et al. postulated a conventional model for the formation of a tRNALys3 packaging complex, in which a Pr55gag/p160gag-pol/viral genomic RNA complex interacts with a tRNALys3/LysRS complex. On the other hand, our findings suggest that GAPDH and tRNALys3/LysRS complex may compete with each other for the interaction to the binding domain within the Pr55gag/p160gag-pol/viral genomic RNA complex (Figure 5). Taken together, these results indicate that GAPDH negatively regulates HIV-1 infection, and small molecules that reconstitute the binding mode of GAPDH to Pr55gag and p160gag-pol may interrupt Pr55gag-LysRS or p160gag-pol-LysRS interactions.

Figure 5
figure 5

Schematic illustration showing the mechanism by which GAPDH negatively regulates HIV- 1 infection. A single copy of each molecule is shown for simplicity. As described in detail elsewhere [20, 26, 27], a Pr55gag/p160gag-pol/viral RNA complex interacts with a tRNALys3/LysRS complex. Although the actual conformation of the Pr55gag/p160gag-pol/viral RNA complex is still unclear, findings of this study suggest that GAPDH and tRNALys3/LysRS complex may compete with each other for interaction with the binding site within Pr55gag and p160gag-pol.

Conclusions

We conclude that GAPDH is incorporated into virions and that this negatively regulates HIV-1 infection by reducing LysRS and tRNALys3 packaging efficiency. Elucidating the role of host proteins in HIV-1 infection is a critical issue and may reveal novel mechanisms of pathogenesis that may lead to the discovery of new antiviral targets.

Methods

Cell culture

A chronically HIV-1LAV-1-infected T-cell line (CEM/LAV-1) and HEK293 cells were maintained at 37°C in RPMI-1640 and DMEM supplemented with 10% fetal calf serum (FCS) containing 100 IU/ml penicillin and 100 μg/ml streptomycin in 5% CO2. TZM-bl cells were obtained from the NIH AIDS Research and Reference Reagent Program. This is a HeLa cell clone that was engineered to express CD4 and CCR5 and contains integrated reporter genes for firefly luciferase and E. coli β-galactosidase under the control of an HIV-1 LTR, permitting sensitive and accurate measurements of infection. PBMCs were isolated using Ficoll-Paque (Pharmacia Biotech, Piscataway, NJ) density gradient centrifugation. For concanavalin A (ConA) stimulation of PBMCs, 2 × 106 cells were incubated for 24 h in RPMI-1640 containing 10% FCS, 4 mM l-glutamine, 2 mM sodium pyruvate, 50 μM 2-ME, and 25 μg/mL ConA (Sigma, St. Louis, MO). After stimulation with ConA, PBMCs were further incubated for 72 h in RPMI-1640 containing 10% FCS, 4 mM l-glutamine, 2 mM sodium pyruvate, 50 μM 2-ME, and 100 U/ml IL-2.

Plasmids

The coding regions of human GAPDH cDNAs (GenBank™ accession number M33197.1) were amplified by PCR using the primers GF (5′-AGGATCCGCCATGGGGAAGGTGAAGGTCGG-3′) and GR (5′-GCCCACATGGCCTCCAAGGAGGATATCA-3′). The regions were cloned into the EcoRV and BamHI sites of the pcDNA™3.1D/V5-His-TOPO® vector (Life Technologies Corporation). The coding regions of the human LysRS cDNAs (GenBank™ accession number NM_001130089) were amplified by PCR using the primers LF (5′-AGGTACCGCCATGTTGACGCAAGCTGCTGTAAGG-3′) and LR (5′-ATTCGAAGACAGAAGTGCCAACTGTTG-3′). The regions were cloned into the KpnI and BstBI sites of the pcDNA™3.1D/V5-His-TOPO® vector (Life Technologies Corporation).

Viruses

Infectious HIV-1LAV-1 stocks were prepared from culture supernatants of CEM/LAV-1 cells [28]. Briefly, the supernatant from the culture medium of CEM/LAV-1 cells was filtered through a 0.22-μm-pore size disposable filter and then centrifuged at 43,000 × g for 3 h at 4°C. The obtained pellet was resuspended in PBS(−) (0.02% KH2PO4, 0.29% Na2HPO4·12H2O, 0.8% NaCl, 0.02% KCl) and then centrifuged at 100,000 × g for 1 h at 4°C. The resulting pellet was resuspended in PBS(−) and used as an infectious HIV-1LAV-1 stock. GAPDH-packaging-defective viruses were prepared by GAPDH siRNA treatment of CEM/LAV-1 cells. CEM/LAV-1 cells (5 × 105 cells) were transfected with 50 nM validated commercially available Silencer GAPDH siRNA (Catalog #:AM4605, the siRNA sequence is not released, Life Technologies Corporation) or 50 nM control siRNA. After 24 h, the CEM/LAV-1 cells were washed and cultured for 48 h with antibiotic-free RPMI containing 10% FCS. The virus-containing supernatant was collected 72 h after transfection. On the other hand, the enhanced-GAPDH- or enhanced-LysRS-packaging virus was prepared by cotransfection of HEK293 cells (5 × 105 cells) with pNL-CH [29] and either the pcDNA™3.1D-GAPDH-V5-His-TOPO® (GAPDH expression vector) or pcDNA™3.1D-LysRS-V5-His-TOPO® (LysRS expression vector). pNL-CH, derived from the pNL4-3 clone of HIV-1, contains a silent T-to-C mutation at nucleotide 2600 to introduce an RsrII restriction enzyme site near the 5′ end of pol.

CD45 affinity depletion

HIV-1LAV-1 was collected from CEM/LAV-1 culture medium and prepared according to previously described methods [28], with slight modifications. The supernatant from the culture medium of CEM/LAV-1 cells was filtered through a 0.22-μm-pore size disposable filter and then centrifuged at 43,000 × g for 3 h at 4°C. The obtained pellet was resuspended in PBS(−) and then centrifuged at 100,000 × g for 1 h at 4°C. The resulting pellet was resuspended in PBS(−). To remove CD45-positive contaminating particles such as microvesicles and exosomes from HIV-1 preparations, HIV-1 preparations containing microvesicles and exosomes were incubated with 2 μl of prewashed CD45 microbeads (Miltenyi Biotech, Auburn, CA) per 1 μg of total protein overnight on ice according to previously described methods [19], with slight modifications. In Figure 1B, forty micrograms (total protein) of purified HIV-1LAV-1 produced from CEM/LAV-1 cells (starting sample) was subjected to the CD45 depletion procedure. The mixture of sample-CD45 microbeads was then applied to a washed and magnetized Macs separation column (Miltenyi Biotech). The flow-through fraction was collected to obtain the CD45-depleted fraction. The column was sufficiently washed with PBS(−). To recover the contaminating particles retained on the beads as the CD45-containing fraction, preheated (95°C) SDS gel loading buffer (62.5 mM Tris–HCl, 2% SDS, 10% glycerol (pH 6.8), and 10% 2-mercaptoethanol) was applied on to the column matrix. To compare the CD45 level between the “CD45-depleted” and “CD45-containing” lanes in Figure 1B, equal amounts of samples by volume were loaded on the gel in accordance with the method of Ott et al. (J. Virol. 2003, 77:12699–12709.). Twenty microliter portions of fractions equivalent to the starting sample (by volume) were analyzed by CD45 immunoblotting (Figure 1B).

Two-dimensional electrophoresis analysis

For 2D gel electrophoresis, the viral pellets derived from HIV-1LAV-1 were boiled for 5 min and lysed with 100 μl of lysis buffer (9.5 M urea, 2% (w/v) NP-40, 2% ampholine (pHs 3–10), and 5% 2-mercaptoethanol). Furthermore, CEM/LAV-1 or CEM cells were also lysed with the same lysis buffer. An aliquot of the lysate (1 mg of proteins for in-gel digestion, or 300 μg of proteins for western immunoblotting) was subjected to 2D polyacrylamide gel electrophoresis according to the method of O’Farrell [30].

Protein identification by mass spectrometry

The gels obtained after 2D gel electrophoresis were stained with SYPRO® Ruby (Invitrogen). SYPRO® Ruby-stained gel pieces were excised from the gel and in-gel-digested with trypsin, as previously described [28]. For the analysis in the MS and MS/MS mode, 1 μl of the matrix (alpha-cyano-4-cinnamic acid saturated in 40% acetonitrile, 0.1% trifluoroacetic acid) per microliter of the tryptic digest was deposited on a ground steel target plate. Analysis was conducted using a MALDI-TOF/TOF UltrafleXtreme (Bruker Daltonics Inc.). Database search was performed on the NCBInr database using Mascot software (Matrix Science, London, UK).

One-dimensional electrophoresis and western immunoblotting

SDS-polyacrylamide gel electrophoresis was performed on 5-20% polyacrylamide gels according to the method of Laemmli [31]. Electrophoresis was carried out at 20 mA per gel for 90 min. Proteins were then transferred onto a PVDF membrane at 0.8 mA/cm2 for 70 min. The membrane was saturated in 5% skimmed milk powder in TBS (10 mM Tris–HCl (pH 7.5), 0.5 M NaCl) and incubated for 4 h with each primary antibody. After washing in TBS-T (10 mM Tris–HCl (pH 7.5), 0.5 M NaCl, 0.5% Tween 20), the membrane was further incubated for 1 h with each secondary antibody, and immunoreaction products were visualized using a SuperSignal® West Pico Chemiluminescent substrate (Thermo Fisher Scientific Inc.).

Suppression of GAPDH packaging by siRNA and virus release assay

To examine the cytotoxicity induced by siRNA treatment, trypan blue dye exclusion assay was carried out. GAPDH siRNA knockdown efficiency in CEM/LAV-1 cells, which were transfected with GAPDH or control siRNA, 72 h after transfection was confirmed by western immunoblotting using an anti-GAPDH antibody (SIGMA). The packaging efficiency of GAPDH, LysRS, and viral structural proteins into virions was monitored by western immunoblotting using an anti-GAPDH antibody, an anti-LysRS antibody (Cell Signaling Technology), and HIV-1-positive plasma (a kind gift from Dr. Shuzo Matsusita, AIDS Research Institute, Kumamoto University, Kumamoto, Japan). The release of the GAPDH-packaging-defective virus was directly monitored by p24 ELISA (ZeptoMetrix Corporation) as the amount of the CA protein in culture supernatants of CEM/LAV-1 cells transfected with GAPDH or control siRNA.

Determination of HIV-1 infectivity

The infectivity of the GAPDH-packaging-defective virus or enhanced-GAPDH-packaging virus was quantified using TZM-bl cells, which express a luciferase gene and a beta-galactosidase gene under the control of the HIV-1 LTR promoter [32]. TZM-bl cells (4 × 105 cells) were incubated with the GAPDH-packaging-defective, enhanced-GAPDH-packaging, or control virus with 20 μg/ml DEAE dextran for 2 h at 37°C and then cultured in DMEM supplemented with 10% FCS (200 μl) for 48 h. The cells were fixed, and the HIV-1 infection of TZM-bl cells was determined by measuring the luciferase activity in cell lysates.

Quantitative analysis of HIV-1 reverse transcription during acute infection

De novo-synthesized HIV-1 cDNA was analyzed using the protocol of Ikeda et al. [33]. Briefly, the TZM-bl cells or PBMCs (1 × 106 cells) were infected with either the GAPDH-packaging-defective virus or the enhanced-GAPDH-packaging virus and incubated for 4 h at 37°C. The cells were washed with PBS(−), incubated for 20 h at 37°C, washed with PBS(−), and further incubated for 5 min at 37°C in PBS(−) containing 0.25% trypsin. After trypsinization, the cells were washed twice with PBS(−) and then digested in 200 μl of digestion buffer (10 mM Tris–HCl (pH 8.0), 150 mM NaCl, 10 mM EDTA, 0.1% SDS, 100 μg/ml proteinase K) for 2 h at 50°C. After digestion, proteinase K was heat-inactivated for 10 min at 95°C. To measure the amounts of early reverse transcription products in GAPDH-packaging-defective virus infection, the sample was subjected to quantitative real-time PCR with a primer pair specific for the R/U5 region (M667, 5′-GGCTAACTAGGGAACCCACTG-3′; AA55: 5′-CTGCTAGAGATTTTCCACACTGAC-3′). To further measure the amounts of late reverse transcription products in enhanced-GAPDH-packaging virus infection, a primer pair specific for the R/gag region (M667, 5′-GGCTAACTAGGGAACCCACTG-3′; M661, 5′-CCTGCGTCGAGAGAGCT CCTCTGG-3′) was used. Because the primer pair R/U5 used detects both early and late products, the following computation was used to determine the amount of early strong-stop DNA: the copy number of strong-stop DNA=R/U5 DNA-R/gag DNA copies.

RT activity assay

To investigate whether GAPDH could allosterically reduce RT activity, recombinant GAPDH (Sigma-Aldrich Co., LLC.) and a reverse transcription assay kit (F. Hoffmann-La Roche Ltd.) were used in this assay. Briefly, the solution (46 mM Tris–HCl, 266 mM potassium chloride, 27.5 mM magnesium chloride, 9.2 mM DTT, digoxigenin (DIG)-labeled dUTP, biotin-labeled dUTP, dTTP, and poly(A) × oligo(dT)15 template/primer hybrid) was added to the reaction tube containing the HIV-1 RT standard preincubated for 1 h with recombinant GAPDH at a ratio of 1:10 or 1:100, and then incubated for 1 h at 37°C. After finishing the RT reaction, the reaction mixture was transferred to streptavidin-coated microtitre plates. DIG-labeled DNA was detected with an anti-DIG-POD conjugate, reacted with 2,2-azino-di(3-ethylbenzthiazoline) sulfonic acid, and quantitated by measuring OD at 405/490 nm. The HIV-1-RT inhibition assay was performed as described in the kit protocol.

Quantification of viral genomic RNA and tRNALys3packaging levels in virions

Both viral genomic RNA and tRNALys3 were collected using a QIAamp® Viral RNA Mini kit (Qiagen) or Nucleo Spin® miRNA (Macherey-Nagel). Genomic RNA was reverse-transcribed using a SuperScript® VILO™ cDNA Synthesis kit and quantified using a primer pair specific for the R/gag region (M667, 5′-GGCTAACTAGGGAACCCACTG-3′; M661, 5′-CCTGCGTCGAGAGAGCTCCTCTGG-3′), or the primers SK38 (5′-ATAATCCACCTATCCCAGTAGGAGAAAT-3′) and SK39 (5′-TTTGGTCCTTGTCTTATGTCCAGAATGC-3′). On the other hand, tRNALys3 was reverse-transcribed by a SuperScript™III First-Strand Synthesis System for RT-PCR (Life Technologies Corporation) using a tRNALys3-F-primer (5′-TGGCGCCCGAACAGGGAC-3′) and then quantified by a tRNALys3-F-primer and a tRNALys3-R-primer (5′-GCATCAGACTTTTAATCTGAGGG-3′). Quantitative real-time PCR was carried out with a SsoFAST™ EvaGreen® Supermix (Bio-Rad Laboratories, Inc.); the cycling conditions were 98°C for 2 min, then 98°C for 5 sec, followed by 40 cycles of 15 sec at 60°C.

Coimmunoprecipitation

CEM/LAV-1 cells were lysed in 350 μl of RIPA buffer (50 mM Tris–HCl (pH 7.4), 150 mM NaCl, 1% sodium deoxycholate, 0.1% SDS, 1% Triton X-100). The insoluble material was pelleted, and the supernatant was used for coimmunoprecipitation. The supernatant was precleaned with Protein-G Sepharose® without any antibody, incubated with 10 μl of an anti-GAPDH antibody (Santa Cruz Biotechnology, INC) or an isotype control goat IgG antibody (Southern Biotechnology Associates, Inc.) for 4 h at 4°C, and further incubated with 50 μl of Protein-G Sepharose® resin slurry (50% slurry in RIPA buffer) for 4 h at 4°C. After washing with RIPA buffer, the bound proteins were eluted using the SDS gel loading buffer. The precipitated proteins were detected by western immunoblotting. Pr55gag and p160gag-pol were detected using an anti-p24 antibody (ViroGen) and an anti-RT antibody (Bio Academia), respectively.

Enhancement of GAPDH packaging by GAPDH expression vector

Enhanced-GAPDH-packaging viruses were prepared by cotransfection of HEK293 cells with pNL-CH and the GAPDH expression vector. Enhanced GAPDH packaging efficiency in viral particles is confirmed by western immunoblotting using an anti-GAPDH antibody (SIGMA). The signal of the CA protein was used as the loading control.

Enhancement of LysRS packaging by LysRS expression vector

Enhanced-LysRS-packaging viruses were prepared by cotransfection of HEK293 cells with pNL-CH and the LysRS expression vector. Enhanced LysRS packaging efficiency in viral particles is confirmed by western immunoblotting using an anti-LysRS antibody (Cell Signaling Technology). The signal of the CA protein was used as the loading control.

Abbreviations

GAPDH:

Glyceraldehyde 3-phosphate dehydrogenase

HIV-1:

Human immunodeficiency virus type 1

RT:

Reverse transcriptase

LysRS:

Lysyl-tRNA synthetase

PBMCs:

Peripheral blood mononuclear cells

MALDI-TOF MS:

Matrix-assisted laser desorption/ionization-time-of-flight mass spectrometry

CA:

Capsid

MA:

Matrix.

References

  1. Ott DE: Cellular proteins in HIV virions. Rev Med Virol. 1997, 7: 167-180. 10.1002/(SICI)1099-1654(199709)7:3<167::AID-RMV199>3.0.CO;2-K.

    Article  CAS  PubMed  Google Scholar 

  2. Cen S, Khorchid A, Javanbakht H, Gabor J, Stello T, Shiba K, Musier-Forsyth K, Kleiman L: Incorporation of lysyl-tRNA synthetase into human immunodeficiency virus type 1. J Virol. 2001, 75: 5043-5048. 10.1128/JVI.75.11.5043-5048.2001.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  3. Franke EK, Yuan HE, Luban J: Specific incorporation of cyclophilin A into HIV-1 virions. Nature. 1994, 372: 359-362. 10.1038/372359a0.

    Article  CAS  PubMed  Google Scholar 

  4. Gabor J, Cen S, Javanbakht H, Niu M, Kleiman L: Effect of altering the tRNA(Lys)(3) concentration in human immunodeficiency virus type 1 upon its annealing to viral RNA, GagPol incorporation, and viral infectivity. J Virol. 2002, 76: 9096-9102. 10.1128/JVI.76.18.9096-9102.2002.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  5. Ott DE, Coren LV, Johnson DG, 2nd Sowder RC, Arthur LO, Henderson LE: Analysis and localization of cyclophilin A found in the virions of human immunodeficiency virus type 1 MN strain. AIDS Res Hum Retroviruses. 1995, 11: 1003-1006. 10.1089/aid.1995.11.1003.

    Article  CAS  PubMed  Google Scholar 

  6. Thali M, Bukovsky A, Kondo E, Rosenwirth B, Walsh CT, Sodroski J, Gottlinger HG: Functional association of cyclophilin A with HIV-1 virions. Nature. 1994, 372: 363-365. 10.1038/372363a0.

    Article  CAS  PubMed  Google Scholar 

  7. Kobbi L, Octobre G, Dias J, Comisso M, Mirande M: Association of mitochondrial Lysyl-tRNA synthetase with HIV-1 GagPol involves catalytic domain of the synthetase and transframe and integrase domains of Pol. J Mol Biol. 2011, 410: 875-886. 10.1016/j.jmb.2011.03.005.

    Article  CAS  PubMed  Google Scholar 

  8. Fermani S, Ripamonti A, Sabatino P, Zanotti G, Scagliarini S, Sparla F, Trost P, Pupillo P: Crystal structure of the non-regulatory A(4)isoform of spinach chloroplast glyceraldehyde-3-phosphate dehydrogenase complexed with NADP. J Mol Biol. 2001, 314: 527-542. 10.1006/jmbi.2001.5172.

    Article  CAS  PubMed  Google Scholar 

  9. Glaser PE, Han X, Gross RW: Tubulin is the endogenous inhibitor of the glyceraldehyde 3-phosphate dehydrogenase isoform that catalyzes membrane fusion: Implications for the coordinated regulation of glycolysis and membrane fusion. Proc Natl Acad Sci U S A. 2002, 99: 14104-14109. 10.1073/pnas.222542999.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  10. Ishitani R, Chuang DM: Glyceraldehyde-3-phosphate dehydrogenase antisense oligodeoxynucleotides protect against cytosine arabinonucleoside-induced apoptosis in cultured cerebellar neurons. Proc Natl Acad Sci U S A. 1996, 93: 9937-9941. 10.1073/pnas.93.18.9937.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  11. Sirover MA: New insights into an old protein: the functional diversity of mammalian glyceraldehyde-3-phosphate dehydrogenase. Biochim Biophys Acta. 1999, 1432: 159-184. 10.1016/S0167-4838(99)00119-3.

    Article  CAS  PubMed  Google Scholar 

  12. Sirover MA: On the functional diversity of glyceraldehyde-3-phosphate dehydrogenase: biochemical mechanisms and regulatory control. Biochim Biophys Acta. 2011, 1810: 741-751. 10.1016/j.bbagen.2011.05.010.

    Article  CAS  PubMed  Google Scholar 

  13. Sirover MA: Role of the glycolytic protein, glyceraldehyde-3-phosphate dehydrogenase, in normal cell function and in cell pathology. J Cell Biochem. 1997, 66: 133-140. 10.1002/(SICI)1097-4644(19970801)66:2<133::AID-JCB1>3.0.CO;2-R.

    Article  CAS  PubMed  Google Scholar 

  14. Volker KW, Reinitz CA, Knull HR: Glycolytic enzymes and assembly of microtubule networks. Comp Biochem Physiol B Biochem Mol Biol. 1995, 112: 503-514. 10.1016/0305-0491(95)00096-8.

    Article  CAS  PubMed  Google Scholar 

  15. Duclos-Vallee JC, Capel F, Mabit H, Petit MA: Phosphorylation of the hepatitis B virus core protein by glyceraldehyde-3-phosphate dehydrogenase protein kinase activity. J Gen Virol. 1998, 79 (Pt 7): 1665-1670.

    Article  CAS  PubMed  Google Scholar 

  16. Schultz DE, Hardin CC, Lemon SM: Specific interaction of glyceraldehyde 3-phosphate dehydrogenase with the 5′-nontranslated RNA of hepatitis A virus. J Biol Chem. 1996, 271: 14134-14142. 10.1074/jbc.271.24.14134.

    Article  CAS  PubMed  Google Scholar 

  17. De BP, Gupta S, Zhao H, Drazba JA, Banerjee AK: Specific interaction in vitro and in vivo of glyceraldehyde-3-phosphate dehydrogenase and LA protein with cis-acting RNAs of human parainfluenza virus type 3. J Biol Chem. 1996, 271: 24728-24735. 10.1074/jbc.271.40.24728.

    Article  CAS  PubMed  Google Scholar 

  18. Petrik J, Parker H, Alexander GJ: Human hepatic glyceraldehyde-3-phosphate dehydrogenase binds to the poly(U) tract of the 3′ non-coding region of hepatitis C virus genomic RNA. J Gen Virol. 1999, 80 (Pt 12): 3109-3113.

    Article  CAS  PubMed  Google Scholar 

  19. Trubey CM, Chertova E, Coren LV, Hilburn JM, Hixson CV, Nagashima K, Lifson JD, Ott DE: Quantitation of HLA class II protein incorporated into human immunodeficiency type 1 virions purified by anti-CD45 immunoaffinity depletion of microvesicles. J Virol. 2003, 77: 12699-12709. 10.1128/JVI.77.23.12699-12709.2003.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  20. Javanbakht H, Halwani R, Cen S, Saadatmand J, Musier-Forsyth K, Gottlinger H, Kleiman L: The interaction between HIV-1 Gag and human lysyl-tRNA synthetase during viral assembly. J Biol Chem. 2003, 278: 27644-27651. 10.1074/jbc.M301840200.

    Article  CAS  PubMed  Google Scholar 

  21. Mak J, Jiang M, Wainberg MA, Hammarskjold ML, Rekosh D, Kleiman L: Role of Pr160gag-pol in mediating the selective incorporation of tRNA(Lys) into human immunodeficiency virus type 1 particles. J Virol. 1994, 68: 2065-2072.

    PubMed Central  CAS  PubMed  Google Scholar 

  22. Fuchigami T, Misumi S, Takamune N, Takahashi I, Takama M, Shoji S: Acid-labile formylation of amino terminal proline of human immunodeficiency virus type 1 p24(gag) was found by proteomics using two-dimensional gel electrophoresis and matrix-assisted laser desorption/ionization-time-of-flight mass spectrometry. Biochem Biophys Res Commun. 2002, 293: 1107-1113. 10.1016/S0006-291X(02)00329-7.

    Article  CAS  PubMed  Google Scholar 

  23. Misumi S, Fuchigami T, Takamune N, Takahashi I, Takama M, Shoji S: Three isoforms of cyclophilin A associated with human immunodeficiency virus type 1 were found by proteomics by using two-dimensional gel electrophoresis and matrix-assisted laser desorption ionization-time of flight mass spectrometry. J Virol. 2002, 76: 10000-10008. 10.1128/JVI.76.19.10000-10008.2002.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  24. Chertova E, Chertov O, Coren LV, Roser JD, Trubey CM, Bess JW, Sowder RC, Barsov E, Hood BL, Fisher RJ, Nagashima K, Conrads TP, Veenstra TD, Lifson JD, Ott DE: Proteomic and biochemical analysis of purified human immunodeficiency virus type 1 produced from infected monocyte-derived macrophages. J Virol. 2006, 80: 9039-9052. 10.1128/JVI.01013-06.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  25. Ott DE, Coren LV, Johnson DG, Kane BP, 2nd Sowder RC, Kim YD, Fisher RJ, Zhou XZ, Lu KP, Henderson LE: Actin-binding cellular proteins inside human immunodeficiency virus type 1. Virology. 2000, 266: 42-51. 10.1006/viro.1999.0075.

    Article  CAS  PubMed  Google Scholar 

  26. Javanbakht H, Cen S, Musier-Forsyth K, Kleiman L: Correlation between tRNALys3 aminoacylation and incorporation into HIV-1. J Biol Chem. 2002, 277: 17389-17396. 10.1074/jbc.M112479200.

    Article  CAS  PubMed  Google Scholar 

  27. Khorchid A, Javanbakht H, Parniak MA, Wainberg MA, Kleiman L: Sequences within Pr160gag–pol affecting the selective packaging of tRNALys into HIV-1. J Mol Biol. 2000, 299: 17-26. 10.1006/jmbi.2000.3709.

    Article  CAS  PubMed  Google Scholar 

  28. Misumi S, Inoue M, Dochi T, Kishimoto N, Hasegawa N, Takamune N, Shoji S: Uncoating of human immunodeficiency virus type 1 requires prolyl isomerase Pin1. J Biol Chem. 2010, 285: 25185-25195. 10.1074/jbc.M110.114256.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  29. Lee SK, Harris J, Swanstrom R: A strongly transdominant mutation in the human immunodeficiency virus type 1 gag gene defines an Achilles heel in the virus life cycle. J Virol. 2009, 83: 8536-8543. 10.1128/JVI.00317-09.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  30. O′Farrell PH: High resolution two-dimensional electrophoresis of proteins. J Biol Chem. 1975, 250: 4007-4021.

    PubMed Central  PubMed  Google Scholar 

  31. Laemmli UK: Cleavage of structural proteins during the assembly of the head of bacteriophage T4. Nature. 1970, 227: 680-685. 10.1038/227680a0.

    Article  CAS  PubMed  Google Scholar 

  32. Derdeyn CA, Decker JM, Sfakianos JN, Wu X, O’Brien WA, Ratner L, Kappes JC, Shaw GM, Hunter E: Sensitivity of human immunodeficiency virus type 1 to the fusion inhibitor T-20 is modulated by coreceptor specificity defined by the V3 loop of gp120. J Virol. 2000, 74: 8358-8367. 10.1128/JVI.74.18.8358-8367.2000.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  33. Ikeda T, Nishitsuji H, Zhou X, Nara N, Ohashi T, Kannagi M, Masuda T: Evaluation of the functional involvement of human immunodeficiency virus type 1 integrase in nuclear import of viral cDNA during acute infection. J Virol. 2004, 78: 11563-11573. 10.1128/JVI.78.21.11563-11573.2004.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We thank Dr. R. Swanstrom (Lineberger Comprehensive Cancer Center, University of North Carolina at Chapel Hill) for providing pNL-CH and helpful discussions. We also thank Dr. Shuzo Matsusita (AIDS Research Institute, Kumamoto University, Kumamoto, Japan) for providing the HIV-1-positive plasma. TZM-bl cells were obtained from the AIDS Research and Reference Reagent Program, Division of AIDS, NIAID, NIH. This work was supported in part by a Grant-in-Aid for Scientific Research from the Ministry of Education, Culture, Sports, Science and Technology of Japan and a Health Science Research Grant from the Ministry of Health, Labour, and Welfare of Japan.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Shogo Misumi.

Additional information

Competing interests

The authors have no conflicting financial interests.

Authors’ contributions

SM conceptualized and designed the study, NK, AO, KK, NT, SS, and SM performed the study, and analyzed data; NK and SM wrote and critically read the paper. All the authors reviewed the manuscript and approved the final version.

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Kishimoto, N., Onitsuka, A., Kido, K. et al. Glyceraldehyde 3-phosphate dehydrogenase negatively regulates human immunodeficiency virus type 1 infection. Retrovirology 9, 107 (2012). https://doi.org/10.1186/1742-4690-9-107

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1742-4690-9-107

Keywords