APOBEC3G mRNA expression in exposed seronegative and early stage HIV infected individuals decreases with removal of exposure and with disease progression
© Vázquez-Pérez et al; licensee BioMed Central Ltd. 2009
Received: 23 December 2008
Accepted: 02 March 2009
Published: 02 March 2009
APOBEC3G is an antiretroviral factor that acts by inducing G to A mutations. In this study, we examined the expression of APOBEC3G in uninfected HIV-1 exposed individuals at the time of their partner's diagnosis and one year later. We then compared this expression with that of infected individuals at different disease stages. APOBEC3G mRNA was measured in PBMCs from three groups: healthy controls with no known risk factor to HIV infection (n = 26), exposed uninfected individuals who had unprotected sex with their HIV+ partners for at least 3 months (n = 37), and HIV infected patients at various disease stages (n = 45), including 8 patients with low HIV viral loads < 10,000 copies/mL (LVL) for at least 3 years. Additionally, we obtained sequences from the env, gag, pol, nef, vif and the LTR of the patients' virus.
Exposed uninfected individuals expressed higher APOBEC3G than healthy controls (3.86 vs. 1.69 relative expression units), and their expression significantly decreased after a year from the HIV diagnosis and subsequent treatment of their partners. Infected individuals showed a positive correlation (Rho = 0.57, p = 0.00006) of APOBEC3G expression with CD4+ T cell count, and a negative correlation with HIV viremia (Rho = -0.54, p = 0.00004). The percentage of G to A mutations had a positive correlation (Rho = 0.43, p = 0.0226) with APOBEC3G expression, and it was higher in LVL individuals than in the other patients (IQR 8.27 to 9.64 vs. 7.06 to 8.1, p = 0.0084). Out of 8 LVLs, 3 had hypermutations, and 4 had premature stop codons only in viral vif.
The results suggest that exposure to HIV may trigger APOBEC3G expression in PBMCs, in the absence of infection. Additionally, cessation of exposure or advanced disease is associated with decreased APOBEC3G expression.
Human APOBEC3G (hA3G) is a cellular antiretroviral factor with a potent inhibitory effect on HIV replication. The cytidine deaminase activity of hA3G catalyses the conversion of cytosine to uracil on the negative-strand viral cDNA. In vitro, hA3G inhibits HIV replication by causing G to A hypermutations, preferentially in the GG dinucleotide context. These mutations often alter the amino acid sequences of multiple viral gene products and introduce lethal stop codons, thereby compromising HIV replication [2–4]. In contrast, APOBEC3F acts preferably in the GA context. In addition to G to A hypermutation, other h3AG antiretroviral mechanisms have been proposed, which include interaction with the HIV nucleocapsid with subsequent inhibition of tRNALys3 annealing to viral RNA, and DNA strand transfer during reverse transcription [2, 3]. Recently, studies in murine models suggest a link between hA3G and virus-specific neutralizing antibody responses in Friend Virus (FV) infection. As a countermeasure to hA3G restriction, HIV-1 neutralizes the antiretroviral activity of hA3G by inducing its ubiquitination and degradation through the Vif protein [5–7]. In vivo, the significance of hA3G-induced hypermutation as a protecting agent of clinical pathogenesis and HIV disease progression remains uncertain. One study suggested that hA3G mRNA levels in activated PBMCs correlated negatively with plasma HIV RNA levels and positively with CD4+ T cell counts. In contrast, a different research group did not find any correlation. On the other hand, one report showed significantly increased hA3G mRNA in PBMCs and cervical biopsy cells from HIV-exposed seronegative individuals (ES), and a separate report found that stimulated PBMCs of long term non-progressors (LNTP) had significantly higher hA3G mRNA levels than either uninfected controls or individuals with progressive HIV disease. Finally, Ulenga et al found that patients with a low viral set point had a higher hA3G expression than patients with high viral set point. Taken together, these studies point to an in vivo HIV neutralizing activity of hA3G.
All the above reports were carried out in different patient groups and experimental conditions, which preclude direct comparisons about differential expression levels between groups like LTNP and ES. Additionally, the question remains whether hA3G is constitutively expressed by these protected patients, or whether its expression is the result of exposure to HIV or its gene products. To address this question, we measured hA3G mRNA expression in unstimulated PBMCs from an ES cohort at the time of first diagnosis of their sexual partners, and after one year of antiretroviral treatment of the infected partner, in order to observe if the expression levels varied with the decrease in viremia of their partners, or remained unchanged. The findings would indicate either induced or constitutive hA3G expression, respectively. In parallel, we contrasted these results with hA3G mRNA expression in HIV infected persons at different disease stages and correlated the expression levels with CD4+ counts and HIV viral load. Separately from these patients, we also studied a subgroup of 8 subjects (low viral load group, LVL) that had remained with HIV viral load < 10,000 copies/mL for at least three years, in order to assess if a sustained control of viremia was related to hA3G expression. Additionally, we measured G to A mutation activity on the viral vif, gag, pol, nef, env and LTR sequences from 20 randomly selected patients. We also included twenty six healthy subjects (HC group), with no known HIV risk factors as a control group for this study.
hA3G mRNA expression in ES, HC and HIV infected patients at different stages of disease
hA3G mRNA expression in ES after one year of partner's diagnosis
Detection of G to A mutation and hypermutation in HIV sequences
The role of hA3G in regulating HIV replication in vivo is unclear. In this study, we demonstrated that individuals exposed to HIV, either through unprotected sex with an infected partner or during a non-AIDS phase of HIV infection, have an increased expression of hA3G mRNA, as shown directly by the expression levels in ES individuals, and by the positive correlation between CD4+ cell count and hA3G mRNA expression. Once the virus overcomes the immune system of the patient, hA3G mRNA expression is almost completely suppressed. This can be seen by the negative correlation between hA3G mRNA expression and HIV plasma viral load, and by the lower hA3G mRNA expression seen in the advanced HIV+ patients.
Further supporting the notion that hA3G mRNA expression is related to exposure to HIV gene products, the viremia levels of the ES's partners were significantly reduced a year after HIV diagnosis and subsequent antiretroviral treatment, coinciding with a reduction in hA3G mRNA expression in the ES subjects. The dramatic reduction in plasma HIV viral load most probably reduced exposure of ES subjects to the viral gene products of their partners. An additional reduction of exposure to HIV gene products would result if the ES individuals adopted safer sex practices, and it is reasonable to assume that at least some of the ES initiated these safer sex practices after being informed of their partner's diagnosis. Even though we cannot fully document a lack of further HIV exposure in our ES cohort, this notion appears to be the best explanation for the reduction in hA3G mRNA expression in the ES individuals.
In this study, we found increased levels of hA3G mRNA in the unstimulated PBMCs of two groups in the same cohort: exposed seronegative and LVL individuals. These increased hA3G mRNA levels are in agreement with previous reports that examined CD3 and CD28-stimulated PBMCs in long term non-progressors and in ES individuals. Here we found that ES and LVL groups have similar values of hA3G mRNA expression, which may suggest that ES individuals regulate hA3G in a similar manner as LVL patients.
The negative correlation of hA3G expression with viremia described in the present study (Figure 2) is in agreement with and clarifies other reports. Jin et al have speculated that HIV infection can induce different discrete stages of hA3G expression, which is in concordance with our findings. Additionally, Ulenga et al.  recently reported that patients with a higher viral set point expressed less hA3G and hA3F mRNAs than low viremic set point patients, also suggesting that hA3G expression is suppressed in advanced disease.
Jin et al suggested that the observed discrete stages of hA3G expression arise by naturally selecting from the general population individuals that constitutively express high levels of hA3G. Our data do not support this contention, since the ES individuals decreased their hA3G mRNA expression after a year of their partners' diagnoses and treatment, and return to healthy control levels, which suggested that they had the same constitutive hA3G expression than that of the general population, but that they induced this expression after coming into contact with HIV gene products. However, there is a high degree of variability in the hA3G mRNA expression in the ES individuals, and there are suggestions that expression may be bounded, so that the final level of expression is the result of the contributions from host genetic susceptibility and the amount of HIV exposure. It remains to be studied if the ES individuals have a genetic susceptibility to have a higher or faster induced expression to initial HIV exposure, and this is what conferred to them their apparent immunity.
The present methodology cannot definitively assess if the reduction in hA3G mRNA in PBMCs is the result of a larger proportion of hA3G expressing cells or an increase in hA3G expression at the single cell level. However, the fact that ES subjects had similar CD4+, CD8+ T cells and monocyte PBMC composition at the basal and one year measurements (data not shown) suggests that the different expression levels are reflected at the single cell level. Furthermore, it has been established that ES individuals can generate anti-HIV-1 T cell responses while remaining uninfected, but lose these responses after ceasing exposure to HIV. T cells produce interleukin-2 (IL-2) in response to stimulation by synthetic HIV-1 env-derived peptides, and Stopak et al. have observed that IL-2 and IL-15 activated de novo hA3G gene expression in primary PBMCs. Thus, the secretion of these cytokines could explain the increased levels of hA3G mRNA in ES individuals.
The mutation data obtained from the 20 individuals showed a significant correlation between the percentage of G to A mutation and hA3G mRNA levels and a higher G to A mutation rate in the LVL subjects. Potentially dysfunctional hypermutated regions were found in three sequences, and four sequences had premature stop codons. This incidence of hypermutated sequences did not correlate with hA3G mRNA expression, as was reported elsewhere . This is possibly due to differences in the sample size and the length of the sequences analyzed. Taken together, the data suggest that the expression of hA3G mRNA is associated with the telltale signs of APOBEC3 activity on the viral genome.
Pillai et al. have hypothesized that there is a dose-dependent response between intracellular hA3G concentrations and the degree of viral hypermutation. Also, they assumed that while an effective inhibition of vif may result in mutational extinction, weak inhibition may accelerate evolution of drug resistance, and immune escape. Our data do not shed light on ether a mutational extinction or accelerated evolution.
Chiu et al.  reported a distinction between the active low molecular weight form (LMW) of hA3G and its inactive high molecular weight form (HMW). In this regard Stopak et al.  had demonstrated that activation by cytokine treatment induced a shift of LMM hA3G to its HMM conformation, paralleled by an increase in susceptibility to HIV infection. On the other hand, Vif could also induce conformational changes of hA3G into HMM complexes [20, 21]. With the methodology described here, we cannot quantify how much mRNA will be translated into LMM or HMM forms, which could be important for viral control. However, we could observe that APOBEC3-like action on the viral genomes correlated with mRNA expression. This could be due to LMM forms induced in dendritic cells .
Recently one report showed that HIV hypermutations correlated with CD4+ counts ; however, we did not find this correlation, and we only found an increased G to A mutation rate, without reaching a statistically significant hypermutated score. In conclusion, our study shows that levels of hA3G mRNA are increased in ES individuals and in early stages of HIV infection, while the levels are decreased a year after their partners' diagnoses and treatment, suggesting that hA3G expression is induced by exposure to HIV.
The study sample consisted of 108 subjects. 26 were healthy controls (HC) free of HIV risk factors, screened with a blood donation questionnaire. 37 were HIV-exposed seronegative individuals (ES) and had a history of multiple unprotected sexual episodes within three months before entering the study. The continued seronegative status was verified for at least one year with subsequent serum samples analyzed by ELISA. Additionally, they were tested for the presence of HIV DNA in PBMCs by gag and LTR nested PCR . 45 were infected patients without antiretroviral treatment (HIV+) and were at different disease stages, with a median of 305 CD4+ cells/μL (IQR 53.25 to 535). These included 8 individuals with low viral loads (LVL), < 10,000 HIV-1 RNA copies/mL for at least three years of follow-up and who were naive to antiretroviral therapy. This study was approved by the Bioethics and Research Committee of the National Institute of Respiratory Diseases, Mexico, and all patients and subjects read and signed informed consent letters.
Relative expression of hA3G mRNA
Peripheral venous blood samples were collected in EDTA tubes; PBMCs were separated by Ficoll density gradient and stored at -80°C. DNA and total RNA was extracted from PBMCs by QIAGEN Blood minikit and RNeasy minikit (Qiagen, Valencia, CA, USA), respectively.
To measure hA3G mRNA expression, cDNA products were first synthesized from the total RNA with random primers using a Transcriptor First Strand cDNA Synthesis Kit (Roche Diagnostics, Mannheim, Germany). cDNA quantification for hA3G and G6PDH (LightMix, for the detection of human G6PDH, Roche) was performed by real time PCR using FastStart DNA Master Hybprobes and LightCycler thermal cycler 2.0 (Roche). Primer and hybprobe design were specific for hA3G (NM 021822), hu ApoB 3G F2 (CAATAATGACATACAGTGAATTT), hu ApoB 3G R2 (CAGGTCTCTGCCTTCCTTAGA), huApoB3GFL (GACATCCCTGGTGGTCCACA-FL) and hu ApoB 3G LC (LC640-GGTGTCCCAGCAGTGCTTAAA-PH). For each sample, the number of copies of hA3G mRNA was divided by the number of copies of G6PDH mRNA, in order to normalize for hA3G mRNA expression in different cell samples. Relative expression units were calculated as expression of hA3G/expression of G6PDH. The results are given as median relative expression units of triplicate assays.
All statistical analyses were carried out with R statistical programming environment version 2.7.1. Statistical differences between groups were assessed with Mann-Whitney U test, correlations were calculated with Spearman's rho, and paired before-after comparisons were carried out with Wilcoxon's signed rank test. hA3G mRNA expression is shown in figures as a logarithm, since it was derived from a ratio, and was therefore log-normal. This transformation did not affect the results of the statistical tests, since they relied on the ranks of the data.
Determination of hypermutated sequences
We amplified and sequenced 5 different regions of HIV proviral DNA: P2, P7, P1 and P6 (gag), protease and RT (pol), V3 loop (env), and the entire vif and nef genes. nef and the V3 loop were amplified as previously reported[25, 26]. Protease and RT sequences were amplified and analyzed by a Viroseq kit (Celera Diagnostics, Alameda, CA, USA). P2, P7, P1 and P6 of gag and the entire vif gene were amplified by nested PCR. Initial amplification of gag was performed using primers p2p7p1p6F outer (5'ATTGGATGACAGAAACCTTGTTGG3') and p2p7p1p6R outer (5'CTTCTAATACTGTATCATCTGCTCC3'), vif initial amplification was performed using primers JVPvifF (5'ACAGCAGAGATCCACT3') and JVPvifR2 (5'AGAATTCTTATTATGGCTTCCA 3'). An aliquot (5 μL) of first round PCR product was then used as a template in a second PCR reaction with primers p2p7p1p6F inner (5'GAAGAAATGATGACAGCATGTC3') and p2p7p1p6R inner (5'CATCTGCTCCTGTATCTAATAG3'), and JVPvifF2 (5'TGGAAAGGACCAGCAAAGCT3') and JVPvifR (CTAGGAAAATGTCTAACAGCTT), respectively. We used High Fidelity Polimerase (Platinum Taq DNA Polymerase High Fidelity, Invitrogen, Carlsbad, California, United States) in all PCR reactions.
Nucleotide sequences of PCR products were determined using BigDye Terminator cycle sequencing kit (Applied Biosystems) and ran on an ABI Prism 3100 analyzer sequencer (Applied Biosystems). Nucleotide sequences of PCR products were aligned and analyzed using Clustal X and MEGA 4.0 respectively. Analysis of G to A substitution in proviral sequences was performed using HYPERMUT 2.0 with default settings, available from http://www.hiv.lanl.gov/content/sequence/HYPERMUT/hypermut.html. This tool can distinguish plus-strand hA3G/F hypermutation, by measuring the number of G to A changes in the GG and GA dinucleotide context, respectively, and testing for a significant increase (p < 0.05) from background levels of mutation.
The present study was supported by Comisión de Equidad y Género de la H. Cámara de Diputados, México. We would like to thank all the patients, and specially the persons of the ES cohort, who kindly cooperated with this study. We thank Dr. Enrique Espinosa for careful reading and critical review of this manuscript.
- Sheehy AM, Gaddis NC, Choi JD, Malim MH: Isolation of a human gene that inhibits HIV-1 infection and is suppressed by the viral Vif protein. Nature. 2002, 418: 646-650. 10.1038/nature00939.View ArticlePubMedGoogle Scholar
- Guo F, Cen S, Niu M, Yang Y, Gorelick RJ, Kleiman L: The interaction of APOBEC3G with HIV-1 nucleocapsid inhibits tRNALys3 annealing to viral RNA. J Virol. 2007, 81: 11322-11331. 10.1128/JVI.00162-07.PubMed CentralView ArticlePubMedGoogle Scholar
- Li XY, Guo F, Zhang L, Kleiman L, Cen S: APOBEC3G inhibits DNA strand transfer during HIV-1 reverse transcription. J Biol Chem. 2007, 282 (44): 32065-32074. 10.1074/jbc.M703423200.View ArticlePubMedGoogle Scholar
- Santiago ML, Montano M, Benitez R, Messer RJ, Yonemoto W, Chesebro B, et al: Apobec3 encodes Rfv3, a gene influencing neutralizing antibody control of retrovirus infection. Science. 2008, 321: 1343-1346. 10.1126/science.1161121.PubMed CentralView ArticlePubMedGoogle Scholar
- Sheehy AM, Gaddis NC, Malim MH: The antiretroviral enzyme APOBEC3G is degraded by the proteasome in response to HIV-1 Vif. Nat Med. 2003, 9: 1404-1407. 10.1038/nm945.View ArticlePubMedGoogle Scholar
- Yu X, Yu Y, Liu B, Luo K, Kong W, Mao P, et al: Induction of APOBEC3G ubiquitination and degradation by an HIV-1 Vif-Cul5-SCF complex. Science. 2003, 302: 1056-1060. 10.1126/science.1089591.View ArticlePubMedGoogle Scholar
- Kao S, Khan MA, Miyagi E, Plishka R, Buckler-White A, Strebel K: The human immunodeficiency virus type 1 Vif protein reduces intracellular expression and inhibits packaging of APOBEC3G (CEM15), a cellular inhibitor of virus infectivity. J Virol. 2003, 77: 11398-11407. 10.1128/JVI.77.21.11398-11407.2003.PubMed CentralView ArticlePubMedGoogle Scholar
- Jin X, Brooks A, Chen H, Bennett R, Reichman R, Smith H: APOBEC3G/CEM15 (hA3G) mRNA levels associate inversely with human immunodeficiency virus viremia. J Virol. 2005, 79: 11513-11516. 10.1128/JVI.79.17.11513-11516.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- Cho SJ, Drechsler H, Burke RC, Arens MQ, Powderly W, Davidson NO: APOBEC3F and APOBEC3G mRNA levels do not correlate with human immunodeficiency virus type 1 plasma viremia or CD4+ T-cell count. J Virol. 2006, 80: 2069-2072. 10.1128/JVI.80.4.2069-2072.2006.PubMed CentralView ArticlePubMedGoogle Scholar
- Biasin M, Piacentini L, Lo CS, Kanari Y, Magri G, Trabattoni D, et al: Apolipoprotein B mRNA-editing enzyme, catalytic polypeptide-like 3G: a possible role in the resistance to HIV of HIV-exposed seronegative individuals. J Infect Dis. 2007, 195: 960-964. 10.1086/511988.View ArticlePubMedGoogle Scholar
- Ulenga NK, Sarr AD, Thakore-Meloni S, Sankale JL, Eisen G, Kanki PJ: Relationship between Human Immunodeficiency Type 1 Infection and Expression of Human APOBEC3G and APOBEC3F. J Infect Dis. 2008Google Scholar
- Huthoff H, Malim MH: Identification of amino acid residues in APOBEC3G required for regulation by human immunodeficiency virus type 1 Vif and Virion encapsidation. J Virol. 2007, 81: 3807-3815. 10.1128/JVI.02795-06.PubMed CentralView ArticlePubMedGoogle Scholar
- Russell RA, Pathak VK: Identification of two distinct human immunodeficiency virus type 1 Vif determinants critical for interactions with human APOBEC3G and APOBEC3F. J Virol. 2007, 81: 8201-8210. 10.1128/JVI.00395-07.PubMed CentralView ArticlePubMedGoogle Scholar
- Jin X, Wu H, Smith H: APOBEC3G levels predict rates of progression to AIDS. Retrovirology. 2007, 4: 20-10.1186/1742-4690-4-20.PubMed CentralView ArticlePubMedGoogle Scholar
- Kaul R, Rowland-Jones SL, Kimani J, Fowke K, Dong T, Kiama P, et al: New insights into HIV-1 specific cytotoxic T-lymphocyte responses in exposed, persistently seronegative Kenyan sex workers. Immunol Lett. 2001, 79: 3-13. 10.1016/S0165-2478(01)00260-7.View ArticlePubMedGoogle Scholar
- Stopak KS, Chiu YL, Kropp J, Grant RM, Greene WC: Distinct patterns of cytokine regulation of APOBEC3G expression and activity in primary lymphocytes, macrophages, and dendritic cells. J Biol Chem. 2007, 282: 3539-3546. 10.1074/jbc.M610138200.View ArticlePubMedGoogle Scholar
- Gandhi SK, Siliciano JD, Bailey JR, Siliciano RF, Blankson JN: The role of APOBEC3G/F-mediated hypermutation in the control of HIV-1 in Elite Suppressors. J Virol. 2007, 82 (6): 3125-3130. 10.1128/JVI.01533-07.PubMed CentralView ArticlePubMedGoogle Scholar
- Pillai SK, Wong JK, Barbour JD: Turning up the volume on mutational pressure: is more of a good thing always better? (A case study of HIV-1 Vif and APOBEC3). Retrovirology. 2008, 5: 26-10.1186/1742-4690-5-26.PubMed CentralView ArticlePubMedGoogle Scholar
- Chiu YL, Soros VB, Kreisberg JF, Stopak K, Yonemoto W, Greene WC: Cellular APOBEC3G restricts HIV-1 infection in resting CD4+ T cells. Nature. 2005, 435: 108-114. 10.1038/nature03493.View ArticlePubMedGoogle Scholar
- Goila-Gaur R, Strebel K: HIV-1 Vif, APOBEC, and intrinsic immunity. Retrovirology. 2008, 5: 51-10.1186/1742-4690-5-51.PubMed CentralView ArticlePubMedGoogle Scholar
- Goila-Gaur R, Khan MA, Miyagi E, Kao S, Opi S, Takeuchi H, et al: HIV-1 Vif promotes the formation of high molecular mass APOBEC3G complexes. Virology. 2008, 372: 136-146. 10.1016/j.virol.2007.10.017.PubMed CentralView ArticlePubMedGoogle Scholar
- Land AM, Ball TB, Luo M, Pilon R, Sandstrom P, Embree JE, et al: HIV-1 Proviral Hypermutation Correlates with CD4 Count in HIV Infected Women from Kenya. J Virol. 2008, 18 (16): 8172-8182. 10.1128/JVI.01115-08.View ArticleGoogle Scholar
- Vazquez Perez JA, Basualdo Sigales MC, Reyes-Teran G, Gudino Rosales JC, Soler CC: Human Immunodeficiency Virus type 1 in seronegative infants born to HIV-1-infected mothers. Virol J. 2006, 3: 52-10.1186/1743-422X-3-52.View ArticlePubMedGoogle Scholar
- R Development Core Team 2008: A language and environment for statistical computing. Vienna, Austria. 2008Google Scholar
- Tanaka M, Ueno T, Nakahara T, Sasaki K, Ishimoto A, Sakai H: Downregulation of CD4 is required for maintenance of viral infectivity of HIV-1. Virology. 2003, 311: 316-325. 10.1016/S0042-6822(03)00126-0.View ArticlePubMedGoogle Scholar
- Delwart EL, Shpaer EG, Louwagie J, McCutchan FE, Grez M, Rubsamen-Waigmann H, et al: Genetic relationships determined by a DNA heteroduplex mobility assay: analysis of HIV-1 env genes. Science. 1993, 262: 1257-1261. 10.1126/science.8235655.View ArticlePubMedGoogle Scholar
- Chenna R, Sugawara H, Koike T, Lopez R, Gibson TJ, Higgins DG, et al: Multiple sequence alignment with the Clustal series of programs. Nucleic Acids Res. 2003, 31: 3497-3500. 10.1093/nar/gkg500.PubMed CentralView ArticlePubMedGoogle Scholar
- Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24: 1596-1599. 10.1093/molbev/msm092.View ArticlePubMedGoogle Scholar
- Rose PP, Korber BT: Detecting hypermutations in viral sequences with an emphasis on G --> A hypermutation. Bioinformatics. 2000, 16: 400-401. 10.1093/bioinformatics/16.4.400.View ArticlePubMedGoogle Scholar
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.