Secondary Logo

Journal Logo

BASIC SCIENCE: CONCISE COMMUNICATION

Comparing genotoxic signatures in cord blood cells from neonates exposed in utero to zidovudine or tenofovir

Vivanti, Alexandrea; Soheili, Tayebeh S.a; Cuccuini, Wendyb,c; Luce, Soniaa; Mandelbrot, Laurentd,e; Lechenadec, Jeromee,f; Cordier, Anne-Gaelg; Azria, Elieh; Soulier, Jeanb,c,i,j; Cavazzana, Marinaa,k; Blanche, Stéphanel,m,*; André-Schmutz, Isabellea,*

Author Information
doi: 10.1097/QAD.0000000000000564

Abstract

Introduction

Guidelines on treatment during pregnancy have gradually expanded from zidovudine monotherapy to various combinations of nucleos(t)ide analogs and either a protease inhibitor or a non-nucleoside reverse transcriptase inhibitor. The tenofovir–emtricitabine nucleos(t)ide analog combination is now widely prescribed (either before conception or during pregnancy) and has recently been recommended by the WHO [1,2].

Potential genotoxicity is a key issue for all currently marketed antiretroviral nucleos(t)ide analogs, since they can all integrate into human nuclear DNA and act as terminators of DNA replication [3–6]. However, the affinity of these nucleos(t)ide analogs for human nuclear and mitochondrial DNA varies from one molecule to another, and possibly from one cell type to another [4]. All of these molecules cross the placental barrier freely and some of them concentrate in the amniotic fluid [7,8]. Thus, exposure in utero might give rise to greater genotoxicity than that observed in adults or older children. This risk merits in-depth evaluation.

We and others identified biomarkers of genotoxicity in neonates exposed to zidovudine–lamivudine [9–12], including a high proportion of aneuploid mononuclear cells in cord blood and the abnormal expression of many DNA repair genes in hematopoietic stem and progenitor cells [12]. To date, few studies on the genotoxicity of tenofovir have been published. Some preclinical data suggest that tenofovir has a weak genotoxicity profile [13,14]. However, in two recent studies, tenofovir–emtricitabine combination was found to be more cytotoxic than zidovudine–lamivudine in vitro, and tenofovir was the most potent inhibitor of telomerase activity [15,16].

Here, we compared the aneuploid rate and the gene-expression profiles in newborns exposed in utero to either fixed-dose combinations of zidovudine–lamivudine or tenofovir–emtricitabine with a ritonavir-boosted protease inhibitor.

Patients, materials and methods

Umbilical cord blood was collected after written consent from non-HIV-1-infected and HIV-1-infected pregnant women having been treated with either a tenofovir or a zidovudine-based combination for at least 4 weeks during pregnancy. The noninfected status of infants born to HIV-1-infected women was checked by PCR. CD34+ hematopoietic stem and progenitor cell (HSPC) and CD34/CD3+ T-lymphocyte were sorted as previously described [12].

CD3+ T cells were karyotyped as previously described [12] according to the 2009 International System for Human Cytogenetic Nomenclature guidelines [17].

The methods used in the transcriptome analysis of CD34+ HSPCs are described elsewhere [12]. For CD3+ T cells, hybridization was performed on HumanHT-12 v4.0 BeadChips (Illumina, Inc., San Diego, California, USA). Analyses were performed as previously described, with a few minor changes [18,19].

Gene set enrichment analysis was performed with GSEA software (http://www.broadinstitute.org/gsea/index.jsp) and reactome pathways derived from the Molecular Signatures Database (http://www.broadinstitute.org/gsea/msigdb/index.jsp).

Quantitative PCR reaction was performed in triplicate with SYBR Green, on 15 and 19 genes selected for validation in CD34+ HSPCs in CD3+ T cells, respectively. The fold-change [2-ΔΔ cycle threshold (CT)] was calculated by normalizing the CT values against the mean values of two housekeeping genes (RPS13 and SEC11a) with stable expression levels according to the microarray data.

Statistical analysis (other than transcriptional profiling analysis)

All statistical analyses were performed with SAS software (version 9.2, SAS Institute Inc., Cary, North Carolina, USA).

Results

Characteristics of patients and controls

Umbilical cord blood was collected for eight HIV-1-uninfected newborns exposed to a zidovudine-based regimen, nine exposed to a tenofovir-based regimen and eight born to uninfected women. The three groups did not differ significantly in terms of maternal age, term of pregnancy or birth weight (data not shown). Infected mothers had no distinctive clinical or biological features relative to the overall data from the French national perinatal cohort. There was no significant difference between the zidovudine and the tenofovir groups in terms of the maternal HIV load. However, the CD4+ cell count was higher in the zidovudine group because the duration of exposure to zidovudine in utero was shorter than that of tenofovir (tenofovir was more frequently initiated before pregnancy than zidovudine).

Cytogenetic studies

The proportion of aneuploid CD3+ T cells in the zidovudine-exposed group (n = 7, mean ± SD: 26.3% ± 9.2) was twice that observed in the tenofovir-exposed group (n = 8, 14.2% ± 6.7; P < 0.05) and in the control group (n = 6, 13.3% ± 7.2; P < 0.05) (Fig. 1a). All chromosomes were involved, and the alterations were randomly distributed (Fig. 1b).

F1-6
Fig. 1:
The frequency of aneuploidy in CD3+ cells stimulated with phytohemagglutinin for 72 h.(a) The aneuploidy rate was defined as the percentage of all scored cells that were aneuploid (i.e. with more or less than 46 chromosomes per metaphase) for each of the antiretroviral-exposed and control samples. Bars represent the mean values. TDF, tenofovir; ZDV, zidovudine. (*) P < 0.05. (b) An example of severe hypoploidy (37X, -Y, -1, -2,-6; -8, -11, -14, -16, -20).

Gene expression profiling of CD3+ T cells and CD34+ hematopoietic stem and progenitor cells

We used microarrays to examine the gene-expression profile of CD3+ cells (n = 7 each for zidovudine, tenofovir and controls). The overall analysis of transcriptome profiles (as illustrated by heat maps; fold-change >1.5; P < 0.05) revealed a characteristic molecular pattern for each group, but with many more similarities between zidovudine and tenofovir samples than between zidovudine or tenofovir and control samples (Fig. 2a). Hierarchical clustering confirmed that the neonates could be segregated according to their in-utero exposure (except for two zidovudine and two tenofovir samples that segregated with control samples; data not shown); however, the number of genes that were differentially expressed when comparing the zidovudine and the tenofovir-exposed groups (3) was much lower than the corresponding values for the zidovudine-exposed vs. control samples (837) and for the tenofovir-exposed vs. control samples (852).

F2-6
Fig. 2:
Gene-expression profiles in cord blood cells.Affymetrix analysis of the gene-expression profiles of CD3+ T cells from zidovudine (azidothymidine)-exposed (ZDV), tenofovir disoproxil fumarate-exposed (TDF) and control groups. Significant differences are based on a 1.5-fold difference in expression and P less than 0.05. (a) The heat map. (b) Enrichment plots and heat maps for telomere maintenance and DNA repair gene sets in GSEA are shown for CD34+ HSPCs. Top left panel: the running enrichment score for the gene set as the analysis walks along the ranked list. Middle left panel: location of the genes from the telomere maintenance reactome within the ranked list. Bottom left panel: a plot of the ranked list of genes. Right panel: a heat map of the core enrichment genes (genes that appear in the ranked list before or at the peak in the enrichment score). The range of colors (red to blue) shows the range of expression values (high to low). (c) Validation of mRNA profiles by qPCR in CD34+ HSPCs (left graph) and CD3+ T cells (right graph) is shown for genes that are differentially expressed in the TDF and ZDV groups. The relative expression (mean ± SD) of genes differentially expressed in ZDV-exposed (dotted bars) and TDF-exposed (hatched bars) groups vs. control group (white bars) is shown. (*) P < 0.05; (**) P < 0.01. GSEA, gene set enrichment analysis; HSPC, hematopoietic stem and progenitor cell; qPCR, quantitative PCR.

Gene set enrichment analysis confirmed that very similar sets of genes were affected by tenofovir and zidovudine in CD3+ T cells (Suppl Table 1, https://links.lww.com/QAD/A668). The expression of genes involved in DNA/RNA synthesis, nucleotide metabolism, cell cycle, DNA repair and telomere maintenance in both the zidovudine and the tenofovir-exposed groups differed significantly from that observed in the control group.

The gene expression profiles of CD34+ HSPCs were similar to those obtained for CD3+ T cells. As illustrated for genes involved in telomere maintenance and DNA repair, further analysis of enrichment plots and heat maps gave similar results when comparing either the tenofovir or the zidovudine-exposed groups with controls (Fig. 2b).

These results suggest that zidovudine and tenofovir have an impact on both mature and immature hematopoietic cells.

Validation of the transcriptome profile with quantitative PCRs

To validate the microarray data, we used quantitative PCRs (qPCRs) to examine the expression of genes selected on the basis of their involvement in DNA repair, chromosome maintenance and the cell cycle; and their significant transcriptional deregulation in both zidovudine and tenofovir-exposed samples.

We found that 43% of the genes analyzed in CD34+ cells and 72% of those analyzed in CD3+ cells were deregulated in zidovudine and tenofovir-exposed samples (P < 0.05), confirming the microarray experiments. Furthermore, qPCR showed that most of the genes involved in the DNA damage response (nucleotide excision repair, mismatch repair, nonhomologous end-joining and homologous recombination) and chromosome maintenance were more strongly expressed in both zidovudine and tenofovir-exposed samples than in control samples (data not shown). When we compared zidovudine and tenofovir-exposed groups, some genes appeared to be expressed differentially (Fig. 2c); for example, CCNE2 was more strongly expressed in CD34+ cells from the zidovudine-exposed samples; and BLM, TP53 and SMC1A were strongly expressed in CD3+ cells from the tenofovir-exposed samples.

Discussion

Our results for the proportion of aneuploid cells in cord blood suggest that tenofovir-based combinations are less genotoxic for the fetus than zidovudine-based combinations. Preclinical data suggest that tenofovir has less clastogenic activity than zidovudine [13,14], which was confirmed here in vivo in newborns exposed in utero. The present results also independently confirmed our previous report of the clear-cut increase in aneuploidy induced by in-utero exposure to zidovudine–lamivudine [12]. Aneuploidy is primarily caused by the occurrence of centrosome and spindle abnormalities during mitosis, and is considered to be a predisposing factor for cancer. Although most aneuploid cells are probably eliminated by apoptosis, one cannot rule out the survival of some of them (notably stem cells) and thus the occurrence of a potentially oncogenic first event [20].

The reassuring aneuploidy data on the tenofovir-based combination are tempered by the results of the transcriptome analyses. We found that tenofovir exposure and zidovudine exposure are associated with altered gene expression in DNA repair and telomere maintenance pathways in CD34+ HSPCs and CD3+ T cells. The expression profiles of some key cell cycle genes also displayed distinctive features. In CD34+ cells, zidovudine exposure was associated with overexpression of CCNE2 (essential for cell cycle control in the late G1 phase and early S phase). It is noteworthy that overexpression of CCNE2 has been reported in mammary epithelial cells treated in vitro with zidovudine [21]. CCNE2 is known to be up-regulated in many tumors [22–24], which may contribute to chromosome instability and even tumorigenesis. SMC1A was down-regulated in CD3+ T cells in the zidovudine-exposed group. The protein encoded by SMC1A is an important part of the kinetochore and is required for cohesion between sister chromatids [25]. It has been established that the frequency of spontaneous chromosome aberrations is significantly higher in SMC1A-mutated cell lines than in control cell lines [26]. We also found that exposure to tenofovir in utero was associated with dysregulated expression of key DNA repair genes, such as BLM and TP53. Despite these specific differences between the zidovudine and the tenofovir-exposed groups, a detailed reactome analysis (using enrichment plots) showed that the two groups exhibited very similar patterns of gene expression dysregulation.

One of the main unresolved questions in our hypothesis-generating study is the persistence over time of this molecular signature and its potential long-term clinical consequences. Olivero et al.[27] demonstrated the presence of centrosome amplification and micronuclei in mesenchymal cells up to 3 years after the in-utero exposure of monkeys to zidovudine. The few studies including children exposed to nucleos(t)ide analogs in utero have small sample sizes and short follow-up periods [28]. A continuous, long-term evaluation via international registries [29] is thus required. The use of drugs (of whatever class) capable of preventing mother-to-child HIV transmission without interfering with the fetal cells’ transcription profile should be preferred.

Acknowledgements

We thank all the families who participated in the study. We thank Fulvio Mavilio, Patrick Revy, Sven Kracker and Josiane Warszawski for fruitful discussions and help with the analyses. We also thank ProfileXpert for technical expertise in transcriptome analysis.

Author contributions: A.V. performed the experiments, analyzed the data and wrote the initial draft of the manuscript. T.S.S., W.C., S.L. and J.L. performed experiments and/or participated in analysis of the data. L.M. helped to design the study, monitored the patients and collected clinical and biological data. A.G.C. and E.A. monitored the patients and collected clinical and biological data. J.S. and M.C. conceived the study and analyzed the data. S.B. and I.A.S. conceived and conducted the study, analyzed the data and wrote the manuscript. All authors have seen and approved the final version.

Funding: This work was funded by the Agence Nationale de Recherche sur le SIDA et les Hepatites. A.V. received a fellowship from the Fondation pour la Recherche Medicale.

Conflicts of interest

There are no conflicts of interest.

References

1. Panel on Treatment of HIV-infected Pregnant Women and Prevention of Perinatal Transmission. Recommendations for use of antiretroviral drugs in pregnant HIV-1-infected women for maternal health and interventions to reduce perinatal HIV transmission in the United States. 2014. http://aidsinfo.nih.gov/contentfiles/lvguidelines/PerinatalGL.pdf.
2. WHO consolidated guidelines on the use of antiretroviral drugs for treating and preventing HIV infection; WHO. http://www.who.int/hiv/pub/guidelines/arv2013/en/. [Accessed 10 October 2014]
3. Olivero OA. Mechanisms of genotoxicity of nucleoside reverse transcriptase inhibitors. Environ Mol Mutagen 2007; 48:215–223.
4. Wutzler P, Thust R. Genetic risks of antiviral nucleoside analogues: a survey. Antiviral Res 2001; 49:55–74.
5. Seier T, Zilberberg G, Zeiger DM, Lovett ST. Azidothymidine and other chain terminators are mutagenic for template-switch-generated genetic mutations. Proc Natl Acad Sci U S A 2012; 109:6171–6174.
6. Olivero OA, Shearer GM, Chougnet CA, Kovacs AA, Baker R, Stek AM, et al. Incorporation of zidovudine into cord blood DNA of infants and peripheral blood DNA of their HIV-1-positive mothers. Ann N Y Acad Sci 2000; 918:262–268.
7. Yeh RF, Rezk NL, Kashuba ADM, Dumond JB, Tappouni HL, Tien H-C, et al. Genital tract, cord blood, and amniotic fluid exposures of seven antiretroviral drugs during and after pregnancy in human immunodeficiency virus type 1-infected women. Antimicrob Agents Chemother 2009; 53:2367–2374.
8. Else LJ, Taylor S, Back DJ, Khoo SH. Pharmacokinetics of antiretroviral drugs in anatomical sanctuary sites: the male and female genital tract. Antivir Ther 2011; 16:1149–1167.
9. Olivero OA, Shearer GM, Chougnet CA, Kovacs AA, Landay AL, Baker R, et al. Incorporation of zidovudine into leukocyte DNA from HIV-1-positive adults and pregnant women, and cord blood from infants exposed in utero. AIDS 1999; 13:919–925.
10. Escobar PA, Olivero OA, Wade NA, Abrams EJ, Nesel CJ, Ness RB, et al. Genotoxicity assessed by the comet and GPA assays following in vitro exposure of human lymphoblastoid cells (H9) or perinatal exposure of mother-child pairs to AZT or AZT-3TC. Environ Mol Mutagen 2007; 48:330–343.
11. Witt KL, Cunningham CK, Patterson KB, Kissling GE, Dertinger SD, Livingston E, et al. Elevated frequencies of micronucleated erythrocytes in infants exposed to zidovudine in utero and postpartum to prevent mother-to-child transmission of HIV. Environ Mol Mutagen 2007; 48:322–329.
12. André-Schmutz I, Dal-Cortivo L, Six E, Kaltenbach S, Cocchiarella F, Le Chenadec J, et al. Genotoxic signature in cord blood cells of newborns exposed in utero to a Zidovudine-based antiretroviral combination. J Infect Dis 2013; 208:235–243.
13. EMA 2005 Scientific Discussion Viread. http://www.ema.europa.eu/docs/en_GB/document_library/EPAR_-_Scientific_Discussion/human/000419/WC500051732.pdf.
14. Brambilla G, Mattioli F, Robbiano L, Martelli A. Studies on genotoxicity and carcinogenicity of antibacterial, antiviral, antimalarial and antifungal drugs. Mutagenesis 2012; 27:387–413.
15. Brüning A, Burger P, Gingelmaier A, Mylonas I. The HIV reverse transcriptase inhibitor tenofovir induces cell cycle arrest in human cancer cells. Invest New Drugs 2012; 30:1389–1395.
16. Leeansyah E, Cameron PU, Solomon A, Tennakoon S, Velayudham P, Gouillou M, et al. Inhibition of telomerase activity by human immunodeficiency virus (HIV) nucleos(t)ide reverse transcriptase inhibitors: a potential factor contributing to HIV-associated accelerated aging. J Infect Dis 2013; 207:1157–1165.
17. Brothman AR, Persons DL, Shaffer LG. Nomenclature evolution: changes in the ISCN from the 2005 to the 2009 edition. Cytogenet Genome Res 2009; 127:1–4.
18. Dai M, Wang P, Boyd AD, Kostov G, Athey B, Jones EG, et al. Evolving gene/transcript definitions significantly alter the interpretation of GeneChip data. Nucleic Acids Res 2005; 33:e175.
19. Ferrari F, Bortoluzzi S, Coppe A, Sirota A, Safran M, Shmoish M, et al. Novel definition files for human GeneChips based on GeneAnnot. BMC Bioinformatics 2007; 8:446–451.
20. Gordon DJ, Resio B, Pellman D. Causes and consequences of aneuploidy in cancer. Nat Rev Genet 2012; 13:189–203.
21. Olivero OA, Tejera AM, Fernandez JJ, Taylor BJ, Das S, Divi RL, et al. Zidovudine induces S-phase arrest and cell cycle gene expression changes in human cells. Mutagenesis 2005; 20:139–146.
22. Payton M, Scully S, Chung G, Coats S. Deregulation of cyclin E2 expression and associated kinase activity in primary breast tumors. Oncogene 2002; 21:8529–8534.
23. Caldon CE, Sergio CM, Burgess A, Deans AJ, Sutherland RL, Musgrove EA. Cyclin E2 induces genomic instability by mechanisms distinct from cyclin E1. Cell Cycle Georget Tex 2013; 12:606–617.
24. Deng J, He M, Chen L, Chen C, Zheng J, Cai Z. The loss of miR-26a-mediated posttranscriptional regulation of cyclin E2 in pancreatic cancer cell proliferation and decreased patient survival. PloS One 2013; 8:e76450.
25. Gruber S, Haering CH, Nasmyth K. Chromosomal cohesion forms a ring. Cell 2003; 112:765–777.
26. Revenkova E, Focarelli ML, Susani L, Paulis M, Bassi MT, Mannini L, et al. Cornelia de Lange syndrome mutations in SMC1A or SMC3 affect binding to DNA. Hum Mol Genet 2009; 18:418–427.
27. Olivero OA, Torres LR, Gorjifard S, Momot D, Marrogi E, Divi RL, et al. Perinatal exposure of patas monkeys to antiretroviral nucleoside reverse-transcriptase inhibitors induces genotoxicity persistent for up to 3 years of age. J Infect Dis 2013; 208:244–248.
28. Benhammou V, Warszawski J, Bellec S, Doz F, André N, Lacour B, et al. Incidence of cancer in children perinatally exposed to nucleoside reverse transcriptase inhibitors. AIDS 2008; 22:2165–2177.
29. The Antiretroviral Pregnancy Registry. http://www.apregistry.com. [Accessed 20 January 2015]
Keywords:

antiretroviral; genotoxicity; HIV; newborn; nucleoside reverse transcriptase inhibitor; pregnancy; telomere; tenofovir; zidovudine

Supplemental Digital Content

Copyright © 2015 Wolters Kluwer Health, Inc.