Zika virus targets the human thymic epithelium

Products Related to Zika, WestNile, Dengue, Malaria, T.B, Chikungunya, HIV, SARS

Product# 10003-10: Recombinant Zika NS1 Protein (E.coli)

Product# 10003: Recombinant Zika NS1 Protein (E.coli)

Product#  10001: Recombinant Zika Envelope Protein (E.coli)

 

Abstract

 

Introduction

Zika virus (ZIKV) epidemics in 2015-2016 resulted in devastating effects, causing microcephaly, other related congenital defects at birth and neurodevelopmental delay after two years in children born from mothers infected by the virus during pregnancy¹–⁴. Additionally, ZIKV infection in adults correlated with a rise in the frequency of cases of Guillain-Barré syndrome⁵,⁶.

Although the knowledge on the cellular and molecular alterations caused by the ZIKV in the nervous tissue largely increased in the last few years²,⁷,⁸, the effects of this virus upon hematopoietic tissues are much less defined.

In terms of secondary lymphoid organs, ZIKV RNA and protein have been described in lymph nodes of Rhesus monkeys in both paracortex and germinal centers⁹,¹⁰. The virus was also found in macrophages, dendritic cells, and B-cells, in both spleen and axillary lymph nodes of this non-human primate. However, in the same study, the presence of ZIKV RNA in T-cells was observed only in axillary lymph nodes from one animal¹⁰.

Much less is known on the putative infection of primary lymphoid organs, and more particularly in the thymus. This central lymphoid organ is responsible for the generation of T lymphocytes under the control of the thymic microenvironment, a three-dimensional cellular network mainly composed by thymic epithelial cells – TEC¹¹. In this respect, it is interesting to note the thymus as a target organ for other RNA viruses, such as HIV¹² and HTLV-1¹³,¹⁴. In fact, we showed that cultured human TEC can be infected by HTLV-1 and convey virus particles to lymphoblastic T cells¹⁵,¹⁶. We thus hypothesized that the thymic epithelium might also be infected by the Zika virus. Herein we provide compelling evidence that ZIKV targets human TEC both in vivo and in vitro. Data are provided showing that ZIKV-infection enhances keratinization of TEC, with a decrease in proliferation and increase in cell death. Moreover, in vitro data revealed that the virus could modulate a high amount of coding and non-coding genes (including miRNAs, circRNAs, and lncRNAs), with upregulation of genes encoding various cell adhesion and cell migration-related proteins. Accordingly, infected TEC enhanced adhesion of lymphoblastic lymphoma-derived T cells and conveyed the virus to this same cell type. Lastly, alterations in the thymic microenvironment from babies prenatally infected by ZIKV were observed, as well as the presence of the viral envelope protein in the thymus. Taken together, our data reveal that the thymus is a target for the ZIKV and may function as a reservoir of the virus during congenital infection.

Go to:

Results

Cultured human TEC can be infected by the Zika virus

We first investigated the infectivity and growth capacity in ZIKV-infected human TEC, evaluating the virus yield in cell monolayers. The human postnatal TEC line used was obtained by explant technique and limiting dilution cloning, being derived from fragments of a postnatal thymus from a child undergoing cardiac surgery¹⁷. Functionally these cells are able to produce cytokines, chemokines, and extracellular matrix proteins, and can adhere to freshly-isolated thymocytes, as well as acute lymphoblastic leukemia derived T cells¹⁶,¹⁸. In the first group of experiments, we applied infection doses of 0.1 and 1.0 MOI, with the cells being then harvested 24, 48 and 72 hours post-infection, and subjected to cytofluorometry or immunofluorescence with the 4G2 antibody for intracellular detection of the viral envelope protein. The mouse monoclonal 4G2 antibody recognizes an epitope on the envelope protein conserved in the flavivirus family, including Dengue virus, West Nile virus, Japanese Encephalitis virus and Zika virus¹⁹. As shown in Fig. 1a,b, we found a progressive increase in the relative numbers of infected cells, as ascertained by the two applied viral doses. Yet, the use of 1.0 MOI was much more efficient in promoting infection, so that around 90% of the cells were 4G2-positive after 72 hours. Accordingly, all further experiments were done using this infection dose at 72 hours post-infection and MOCK cells were used as control, as the 4G2 antibody present a slight background in MOCK cells when compared to the isotype control. The same percentage of infection was detected by immunofluorescence on adhered cells (Fig. 1c,d), and the presence of the ZIKV envelope protein could be detected in the cytoplasm - especially around the nucleus (Fig. 1e).

 

 

Open in a separate window

Figure 1

Human TEC can be infected in vitro by the Zika virus. (a) TEC were infected with ZIKV (MOI = 0.1 or 1) and the relative expression of 4G2⁺ cells was detected 24, 28 and 72 hours post-infection (hpi) by flow cytometry. Each time point represents the mean ± standard error. Asterisks represent statistical significance between MOCK and MOI = 1 (24, 48 and 72 h) and between MOCK and MOI = 0.1 (72 h). Hash marks represent statistical significance between MOI = 0.1 and MOI = 1. (b) Representative histograms of the relative expression of 4G2⁺ cells 24 and 72 hpi with MOI = 1. Orange curves represent isotype controls, blue curves represent MOCK and red curves represent ZIKV (n = 3). TEC were infected with ZIKV (MOI = 1) and (c) the percentage of 4G2⁺ cells were detected 72 hpi by immunofluorescence. Results are represented by the mean percentage of infection of the three replicates of each independent experiment. Representative images of immunofluorescence for 4G2 (viral protein, in red), cytokeratin (CK, in green) and DAPI (nuclei, in blue) (d) in lower and (e) higher magnifications (n = 3, in triplicates). (f) Representative images of TEC culture 72 hpi (n = 3). (g) ZIKV growth curve in TEC culture. Supernatant of ZIKV-infected TEC culture was harvested, and the presence of infective viral particles was verified in Vero cells. Results are shown as plaque-forming units and each time point represents the mean ± SD, of at least 5 wells (n ≥ 5). (h) Relative expression of 4G2 ⁺ cells detected by flow cytometry in TEC culture 72 hpi with 1 mL of supernatant from MOCK or ZIKV-infected TEC cultures (n = 3). (i) Relative expression of AXL, TIM1, TIM4 e TYRO3 receptors detected by flow cytometry in TEC (n = 3). (j) Relative expression and (k) Mean Fluorescence Intensity (MFI) of 4G2⁺ cells detected by flow cytometry 72 hpi (MOI = 1) after treatment with anti-AXL polyclonal antibody (10 µg/mL) (n = 3). *or p < 0.05, **or p < 0.01, ***p < 0.001, ****or p < 0.0001, (unpaired t-test or one-way ANOVA, followed by Tukey multiple comparison test).

ZIKV infection induced cytopathic effects in TEC with the presence of higher amounts of floating cells in the culture and vacuole-like structures around the nucleus (Fig. 1f).

We also showed that infection of growing TEC resulted in the generation and release of viral infectious particles in the culture supernatants. Following 24, 48 and 72 hours post-infection, the obtained viral titers from TEC cultures were 4.29 ± 0.16, 6.03 ± 0.25, and 6.96 ± 0.11 log10 PFU/mL, as ascertained in Vero cell cultures (Fig. 1g). Furthermore, infected TEC-derived supernatants were able to convey ZIKV to uninfected TEC cultures (Fig. 1h).

Since human TEC could be infected by ZIKV, we searched for the membrane expression of well-known ZIKV receptors, namely TIM1, TIM4, TYRO3, and AXL. Actually, cultured human TEC constitutively expresses AXL, but not TIM1, TIM4, TYRO3 receptors (Fig. 1i). Blocking the AXL receptor with the corresponding neutralizing antibody completely prevented TEC infection by ZIKV (Fig. 1j,k), showing that ZIKV uses the AXL receptor to invade human TEC.

Zika virus infection disturbs cell growth in cultured human TEC

Considering that human TEC were sensitive to ZIKV infection, we wondered whether such an infection could also modulate the growth of these cells in vitro. The TEC counting in infected versus uninfected cultures showed that ZIKV infection induced a progressive reduction in the numbers of growing TEC (Fig. 2a). Experiments aiming at defining whether such a diminution was due to a decrease in cell proliferation or increase in cell death, revealed a diminution in the proportion of dividing cells (defined by the high expression of the Ki67 marker) and an increase in the relative numbers of apoptotic (but not necrotic) cells (Fig. 2b,c), as ascertained by PI and annexin-V double labeling. Differences in the proportion of apoptotic and dividing cells begin to appear 48 hours post-infection, although significant differences are only detected 72 hours post-infection. Of note, in addition to the small, but significant increase in the relative numbers of apoptotic cells, we observed increased amounts of cellular debris in ZIKV-infected TEC cultures when compared to controls (data not shown). It is interesting to note that such modulation of cell growth was not seen when we infected the same human TEC line by HTLV-1 neither in cell proliferation nor cell death¹⁶.

 

 

Figure 2

Zika virus infection down modulates the growth of cultured human thymic epithelial cells. TECs were infected with ZIKV (MOI = 1) and the (a) numbers of cells in culture were detected 24, 48 and 72 hours post-infection (hpi). Each time point represents the mean ± SEM (n = 3). (b) Percentage of apoptotic cells (annexin-V⁺PI⁻) detected by flow cytometry 24, 48 and 72 hpi. Each time point represents the mean ± SEM (n = 3) (c) Percentage of Ki67 high⁺ cells detected by flow cytometry 24, 48 and 72 hpi. Values represent the mean ± SEM (n = 3). **p < 0.01, ***p < 0.001 (unpaired t-test).

Zika infection induces major changes in human TEC gene expression profile

The transcriptome analysis of ZIKV-infected TEC revealed 2,934 differentially expressed genes (DEGs) with p-adjusted < 0.01 and Log₂ Fold Change > |1|, distributed in protein-coding genes, long non-coding RNAs, circRNAs and other RNAs classes (Table 1, Supplemental Table ¹).

Table 1

Summary of the differentially Expressed Genes and miRNAs identified by RNAseq and small RNAseq in ZIKV-infected TEC.

Molecule	Expressed Genes	Upregulated*	Downregulated*	Total*
Protein Coding	15935	1471	106	1577
Long Non-Coding RNA	4113	629	26	655
circRNAs	13342	562	91	653
other RNA Classes**	321	3	46	49
Total DEGs	33711	2665	269	2934
miRNAs	591	42	25	67

*Statistical significance of DESEq. 2 with p-value adjusted < 0.01 and Log2 Fold Change > |1.0|; **miscRNA, scaRNA, scRNA, snoRNA and snRNA were included in this category.

The enrichment analysis of DEGs in ZIKV-infected TEC allowed to identify upregulated pathways involved with interleukin and interferon signaling, extracellular matrix organization, integrin cell surface interactions, NCAM1 interactions, NCAM signaling and ECM proteoglycans in Reactome enrichment (Fig. 3a,b, Supplemental Table ²). The GO enrichment analysis revealed categories also associated with extracellular matrix, response to virus and bacteria, cytokine production, defense response to virus (Supplemental Table ²). In addition, we found genes with enhanced expression by ZIKV-infected TEC such as CCL4, CCL5, CCL20, CXCL10, GBP2, GBP4, GBP5, GBP6, IFNB1, IFNL1, 2 and 3, IL12B, IL18RAP, MX1, MX2, NOS2, NOD2, OAS2, OASL, RSAD2, and TNF. Other upregulated antiviral genes were found, including IFIH1 and DDX58, which act as cytoplasmic sensors of viral nucleic acids, and DHX58, which mediates antiviral signaling and induction of type I interferons and proinflammatory cytokines (Supplemental Tables ¹ and ²).

 

Open in a separate window

Figure 3

Enrichment Analyses of upregulated and downregulated genes from transcriptome of Zika virus-infected TEC. (a) Top 10 significant terms from the Reactome Enrichment Analysis of the upregulated genes. Horizontal red bars indicate the p-adjusted values (−Log₁₀ scale) from the Reactome Enrichment Analysis. Numbers in the top of the red bars indicate the gene counting of each Reactome Term. (b) Cnetplot showing both upregulated genes and enriched Reactome terms. Pale yellow circles indicate the size of Reactome terms. The Log₂ Fold Changes are indicated by the dot colors. Warmer colors indicate those most upregulated genes. (c) Top 10 significant terms from the Reactome Enrichment Analysis of the downregulated genes. Horizontal blue bars indicate the p-adjusted values (−Log₁₀ scale) from the Reactome Enrichment Analysis. Numbers in the top of the blue bars indicates the gene counting of each Reactome Term. (d) Cnetplot showing both downregulated genes and enriched Reactome terms. Pale yellow circles indicate the size of the Reactome terms. The Log₂ Fold Changes are indicated by the dot colors. Greener colors indicate the most downregulated genes.

The Enrichment Analysis of Reactome also revealed upregulated genes included in extracellular matrix organization, and recognized as adhesion molecules, such as ADAMTS4, ACAM, ICAM1, ICAM4, VCAM1, ITGAM, PECAM1, CEACAM1, JAM2, MMP8, ITGA9, LAMA2 LAMA3, LAMA4, LAMC3, COL2A1, COL20A1, and COL26A1 (Fig. 3b). Interestingly, the GO enrichment revealed the category named as collagen fibril organization enriched, containing upregulated collagen fibril forming genes (COL2A1, COL3A1, COL5A1, COL5A2 COL5A3, COL11A1, and COL11A2), as shown Fig. 3 and Supplemental Table ².

We also found upregulated genes that integrate extracellular matrix organization, interferon signaling, interleukin signaling, and integrin cell surface interaction, including ITGAM, VCAM1, ICAM1 in Reactome enrichment (Fig. 3). The transcriptome results were coherent with the cell growth assay that results in progressive reduction in the numbers of growing TEC, reduced cell proliferation and cell division, with a slight increase of apoptosis after ZIKV infection (Supplemental Table ² and Fig. 2).

The Enrichment Analysis also revealed downregulated protein-coding genes related to histone 1 and 2 cluster, ubiquitin, cyclin, E2F Transcription Factor 2 (Fig. 3c, Supplemental Table ²). The downregulation of FBXO43, CKS2, AURKA, CCNB1, CDC20, UBE2C, CCNF, and CCNA2 genes could also contribute to the decreasing of cell growth (Fig. 2, Supplemental Tables ¹ and ²) corroborating the reduction of cell proliferation in vitro. Moreover, we found downregulated genes (CCNB1, AURKA, PMAIP, and PLK3), which are also associated with the response to DNA damage.

Long non-coding RNAs involved with an antiviral response in ZIKV-infected TEC

The DGE analysis also identified 655 differentially expressed long non-coding RNAs (lncRNAs), with 629 upregulated and 26 downregulated in ZIKV-infected TEC (Table 1). Comprising among others, three upregulated lncRNAs, previously related with antiviral response (LUARIS, EGOT, TNFRSF14-AS1), two associated with regulation of apoptosis (PANDAR, MBNL1-AS1), two associated with both apoptosis and proliferation (UCA1, ZEB2-AS1), four with cell proliferation (MDC1-AS1, HOXB-AS3, TFAP2A-AS1, PICART1), one with cell cycle arrest (BDNF-AS) and one with cell adhesion (SENCR). Those lncRNAs showed log2 Fold Change ranging from 1.05 to 5.40 (Fig. 4a, Supplemental Table ¹).

 

 

Figure 4

Non-coding genes response in human TEC infected with the Zika virus. (a) Expression of Long non-coding genes related to Immune System, Apoptosis, Apoptosis and Proliferation, Proliferation, Cell Cycle Arrest and Cellular Adhesion. (b) Network formed by MicroRNA and respective target transcripts, showing upregulated microRNAs and down-regulated targets. (c) Network formed by MicroRNA and respective target transcripts, showing down regulated microRNAs and upregulated targets. The Log₂ Fold Changes were indicated by the continuous scaled colors. Warmer colors indicate upregulation, cooler colors indicate down regulation. (d) Expression analysis of Hsa-miR-542-3p and its putative circRNAs targets. (e) Expression analysis of Hsa-miR-92a-1-5p and its putative circRNAs targets. The data are expressed in Log2 Fold Change (Log2FC). Red and blue bars indicate up and downregulated genes, respectively.

miRNAs and circRNAs expression patterns reveal a fine control of transcription in ZIKV-infected TEC

Small RNA-seq analysis revealed 67 miRNAs differentially expressed in ZIKV-infected TEC, resulting in 42 upregulated and 25 downregulated microRNAs (Table 1 and Supplemental Table ¹). The hsa-miR-7-3-3p was the most upregulated miRNA (log2 FC = 5.16) after ZIKV infection. The majority of miRNAs presented a slightly upregulation or downregulation with log2 FC > |1| including hsa-miR-155-5p (1.55) and hsa-miR-203a-3p (1.19).

The microRNA-target network (Fig. 4b,c) revealed 13 upregulated microRNAs that can modulate the expression of 12 downregulated target transcripts, including the hsa-miR-3065-5p regulating HIST2H2BE and MYB, hsa-miR-155-5p regulating MYB and E2F2, hsa-miR-449a-5p regulating HIST1H4H and hsa-miR-215-5p regulating CDCA4 (Fig. 4b). This analysis also showed 20 downregulated microRNAs, which target upregulated genes involved with immune system, viral response and extracellular matrix, as shown in hsa-miR-4458-5p that regulates the IL-6, CXCL8, C2, COL3A1, COL5A2, ITGB8, TNFSF9 genes, hsa-miR19b-5p targeting ITGB8 and CXCL8, hsa-miR-100-3p targeting CXCL8 and hsa-miR-222-5p regulating HLA-B and DDX58 genes (Fig. 4c).

The transcriptome analysis also revealed 562 upregulated and 91 downregulated circRNAs in ZIKV-infected TEC (Table 1). In order to identify the possible role of microRNA sponge by circRNAs, we conducted a microRNA target prediction with miRanda tool, using the circRNAs as targets. This strategy revealed the upregulated hsa-miR-542-3p targeting 3 circRNAs (Fig. 4d) and the downregulated hsa-miR-92a-1-5p, possibly regulated by 12 circRNAs (Fig. 4e).

Taken together, the transcriptome and small RNA-seq analyses revealed complex molecular responses of TEC in response to ZIKV infection, showing an intricate regulatory network, including different layers of transcription regulation. In this context, ICAM-1 and VCAM appeared as a main hub connecting distinct clusters of genes whose expression was altered secondary to ZIKV infection.

Co-cultures of ZIKV-infected TEC with lymphoblastic T cells enhance heterocellular cell adhesion and virus transfer

Considering that the expression of cell adhesion related-genes was enhanced in ZIKV-infected TEC, and the fact that ICAM-1 as well as VCAM are adhesion molecules known to play a role in TEC thymocyte interaction²⁰, we investigated whether these proteins were upregulated in the membrane of ZIKV-infected TEC. We found enhanced numbers of TEC expressing ICAM-1 (CD54 – data not shown), as well as enhanced ICAM-1 in infected cells (Fig. 5a). This was more clear-cut when we analyzed the cells expressing high amounts of ICAM-1 (5.2% (±1.2) of ZIKV infected cells and 2.6% (±0.5) for controls (p = 0.05). As regards the expression of VCAM-1 (CD106), we did find higher membrane expression levels in infected cultures, although the differences were not statistically significant (p-value = 0.1).

 

 

 

Figure 5

Zika virus infection in human cultured epithelial cells enhances cell adhesion and conveys the virus to lymphoblastic T cells. (a) Human TEC were infected with ZIKV (MOI = 1) and 72 hpi the expression of ICAM-1 (CD54) and VCAM-1 (CD106) were analyzed by flow cytometry. Graphs show the Mean Fluorescence Intensity (MFI) of these molecules. MFI values were normalized by the mean of the control group (n = 3). Human TEC were infected with ZIKV (MOI = 1) and 72 hpi uninfected CEM acute T-cell leukemia cell line was added to the culture. (b) Association index between CEM and TEC after 1 h of co-culture (adhesion). Values in the upper panel represent mean ± SEM. Representative microscopic fields after Giemsa staining are shown in the lower panel (n = 4). Relative expression of 4G2⁺ cells in (c) TEC and (d) CEM subpopulations were detected by flow cytometry after 24 h of co-culture. Values represent the mean ± SEM (n = 3, two experiments in duplicate). (e) CEM cells were not directly infected by ZIKV (MOI = 1). CEM cells were infected with ZIKV (MOI = 1) for 2 h and 72 hpi 4G2⁺ cells were analyzed by flow cytometry. Values represent the mean ± SEM (n = 3). *p < 0.05, **p < 0.01, ****p < 0.0001 (as determined by the unpaired t-test).

We searched for a functional change on the ability of infected TEC to bind to lymphoblastic T cells. For that, the CEM acute T-cell leukemia cell line, previously known to adhere to human TEC¹⁶, was led to adhere to growing TEC 72 hours post-infection. Figure 5b shows that adhesion was significantly increased onto ZIKV-infected TEC cultures, as compared to controls. Importantly, we showed that 72 hours after TEC infection plus 24 hours of co-culture with CEM, a proportion (10.76% ± 1.98) of the T lymphoblasts were infected as revealed by immunolabeling with the 4G2 antibody (Fig. 5d). At this time point (96 hours post-infection), 86% of the TEC present in the co-culture were infected (Fig. 5c). Interestingly enough, we could not directly infect CEM cells with ZIKV (Fig. 5e) or with Zika virus-infected TEC supernatant (data not shown).

The human thymic epithelium is altered after congenital in vivo ZIKV infection

We also studied sections of six thymuses from children that died in postnatal life born from mothers infected by ZIKV during pregnancy. The histological profile of these specimens varied from virtually normal cortico-medullary pattern to a severe disruption of the cortex in one thymus, with clear-cut decrease in thymocytes (Fig. 6a). Yet, in all cases, we observed an increase in the numbers of the Hassall’s corpuscles, indicating an over differentiation of the TEC network. This feature was further confirmed by immunohistochemistry for cytokeratin detection (applied herein as a pan-marker of the whole TEC network) in TEC cultures, as seen in Supplemental Fig. ¹. In fact, cytokeratin immunostaining revealed, in some thymus samples, a complete loss of the cortical/medullary distinction (Fig. 6b, case 2).

 

Open in a separate window

Figure 6

Congenital Zika virus infection induces major changes in the human thymic epithelial network. Representative images of (a) Hematoxylin-eosin and (b) cytokeratin staining of paraffin-embedded human thymus sections. A control case, case 6 with intermediate alterations and case 2 with severe alterations in the human thymic epithelial network are represented (n = 2 for controls and n = 6 for congenital ZIKV infection). Arrows show some of the Hassal’s corpuscles. C: cortical region; M: medullary region.

Considering the transcriptome data, we also performed histological evaluations for ECM-containing fibers and fibrils. In the thymuses from ZIKV-infected newborns, there was an important increase in the reticulin network, indicating thickness of basement membranes in the thymic microenvironment. Additionally, Masson’s trichrome staining revealed a higher density of type I collagen-containing fibers in the thymic septae. These data are summarized in Fig. 7.

 

 

 

Open in a separate window

Figure 7

Congenital Zika virus infection induces changes in the human thymic extracellular matrix. Representative images of (a) Gomori’s reticulin (dark blue fibers) and (b) Masson’s trichrome staining (blue fibers) of paraffin-embedded human thymus sections. A control case, case 6 (with intermediate alterations) and case 2 (with severe alterations) in the human thymic epithelial network are represented. Arrows show some collagen fibers present in the images. (n = 2 for controls and n = 6 for congenital ZIKV infection).

We also searched for the presence of the viral envelope protein in the same thymuses, using the 4G2 antibody, which recognizes any Flavivirus family envelope protein. As compared to normal thymuses (that were completely negative for 4G2 labeling), we did find specific immunostaining in thymus sections from 5 out of 6 thymuses evaluated. The presence of the envelope viral protein was seen in scattered TEC, both in the cortex and medulla of the thymic lobules, as well as in Hassall’s corpuscles (Fig. 8). It is important to point out that one thymus was positive for the presence of ZIKV ascertained by the RT-PCR (case 4 - Supplemental Table ⁴), although we were not able to detect the viral envelope protein in this case (case 4 – Fig. 8 and Supplemental Fig. ^2a). None of the thymus preparations were positive for the presence of Dengue as well as Chikungunya viruses (Supplemental Table ⁴). Secondary antibodies were used as negative controls and did not generate any significant labeling (Supplemental Fig. ^2b).

 

Open in a separate window

Figure 8

ZIKV envelope protein is found in human thymus of babies with congenital Zika virus infection. Representative images of immunofluorescence staining of 4G2 (viral protein, in red), cytokeratin (in green) and DAPI (nuclei, in blue) of paraffin-embedded human thymus sections (n = 2 for controls and n = 6 for congenital ZIKV infection). First column shows the merge of the three stainings. The last column shows the colocalization analysis for 4G2 and cytokeratin stainings; white points represent overlap between these two immunostainings. We could identify viral protein in five out of the six thymuses studied; most of cells labeled being epithelial. Arrows show some of the cells that were stained for both 4G2 and cytokeratin. Differences in bar sizes are due to the digital zoom used during image acquisition.

Go to:

Discussion

The behavior of the human immune system secondary to a congenital syndrome caused by the ZIKV remains largely unknown, although there is evidence that innate and acquired immune cells can be infected by the virus, as ascertained by studies on peripheral blood leukocytes²¹,²². Nevertheless, whether primary lymphoid organs are targeted by the ZIKV, as well as putative consequences upon interactions with B- and T-cell precursors and their corresponding microenvironments, remain to be demonstrated. Herein we tackled this issue by investigating whether or not the human thymic epithelium could be infected by ZIKV.

We first demonstrated that growing human TEC could indeed be infected, and that such infection occurred through the AXL receptor. Moreover, viral particles produced within TEC were infective. Interestingly, the human TEC line could also be infected with other RNA viruses, namely the retrovirus HTLV-1¹⁵. Moreover, pioneer work has shown that primary cultures of human TEC can also be infected with HIV¹². This shows that human TEC is permissive to different types of RNA viruses.

Importantly, ZIKV did trigger an antiviral machinery response in human TEC, as ascertained by the upregulation in the expression of several genes related to antiviral response, including, among others, those related to the type 1 interferon signaling.

Major changes in TEC growth were also found after in vitro ZIKV infection, including a decrease in proliferation and increase in cell death, resulting in a reduction of total cell numbers after 3-day TEC cultures. RNA viruses such as ZIKV use strategies to modulate cell-cycle control and cell proliferation from the host and obtain cellular conditions favorable for their replication²³,²⁴. Accordingly, the data from RNA sequencing revealed few downregulated genes associated with the regulation of cell proliferation and apoptosis. These results are consistent with previous studies on ZIKV-infected human cell types, such as Sertoli cells in the seminiferous tubules, retinal pigment epithelium, neural progenitors, and astrocytes, which reported suppression of cell proliferation, as well as induction of antiviral defense pathways²⁴–²⁷.

It is known that adhesion of developing T-cells to the thymic epithelium is controlled by various cell-cell and cell-extracellular matrix interactions, including ICAM-1, VCAM-1, fibronectin, laminins, and their corresponding integrin receptors²⁸–³². We found a large number of upregulated genes related to extracellular matrix, cell adhesion, and cell migration. Significant increase in the gene expression of various laminin alpha-chains as well as fibril forming collagen chains was seen in ZIKV-infected TEC, as well as higher levels of VCAM1 and ICAM-1 expression in TEC after ZIKV infection.

Enhanced expression of some matrix metalloproteinase genes, such as MMP8 and MMP19 (both having collagenolytic activity) was also observed, suggesting that the ECM-related metabolism of TEC is altered by ZIKV. Moreover, ZIKV augmented the gene expression of various chemokines; some of which (CCL4, CCL5, CCL20, and CXCL10, for example) being involved in intrathymic T-cell migration³³–³⁶.

In keeping with the transcriptome data related to extracellular matrix, cell adhesion, and cell migration, we demonstrated that infected TEC expressed higher levels of the ICAM-1 protein on their membranes. In keeping with the literature showing that ICAM-1 plays a major role in TEC adhesion to developing T-cells³⁷, we found that lymphoblastic T cells adhered more to ZIKV-infected TEC than to uninfected TEC cultures. We can thus hypothesize that enhanced expression of adhesion molecules and chemokines would facilitate the encounter of developing T-cells with ZIKV-infected TEC, favoring in one hand the transfer of the virus to lymphocytes, and on the other hand, the establishment of an immunological synapse between the two cell types. This should also be placed in the context of the control of HLA-B gene expression, which was upregulated (Fig. 3b) in ZIKV-infected human TEC.

The control of gene expression in ZIKV-infected human TEC was also evaluated through the analysis of different non-coding RNAs. In fact, our robust transcriptome and microRNA expression analyses in ZIKV-infected TEC highlighted the molecular response of protein-coding genes and non-coding genes (lncRNAs, circRNAs and miRNAs). We observed miRNAs acting as probable regulators of genes associated with immune response, such as hsa-miR-222-5p (DDX58 and HLA-B), hsa-miR-1303-5p (HLA-DOA), hsa-miR-193b-5p (IFIT2), hsa-miR-100-3p (IFIT1B). Some miRNAs also have shown response to Cytokine Signaling in Immune system such as hsa-miR-503-3p (TNFRSF9) and hsa-miR-4458-5p (TNFSF9). The chemokines CXCL8 and CXCL11 could also be regulated by miRNAs in ZIKV-infected TEC. Taken together, these analyses indicate that miRNAs may play a relevant role in the intracellular control of antiviral immunity against ZIKV-infected TEC.

We also identified several cirRNAs which can serve as miRNA sponges and are crucial regulators of gene expression³⁸. Some circRNAs present antiviral activity, interacting directly with the antiviral genes DDX58 (RIG-I), IFIH1 (MDA5) and TLR3³⁹. The major class of non-coding lncRNAs also shows antiviral activity and involvement in cell proliferation and apoptosis. The LUARIS lncRNA expression level was decreased, indicating that it may be suppressed by IFNB1, which, in turn, was highly induced in ZIKV-infected human TEC. The expression of some interferon stimulated genes (ISGs) and chemokines showed altered levels of transcripts, indicating that LUARIS could act in antiviral immunity. The same regulation pattern of ISGs and chemokines was found in infection by other human viruses such as HBV, HCV, and EMCV⁴⁰. Following ZIKV infection in TEC, we also detected an upregulation of the lncRNA EGOT, which is also induced after HCV infection, increasing viral replication⁴¹.

Studies involving lncRNA associated with cell proliferation and apoptosis and cell adhesion have been already described. We hypothesize that upregulation of PICART1 may also suppress proliferation and induce apoptosis in ZIKV-infected human TEC, as shown by other authors that this lncRNA impairing proliferation⁴² and inducing apoptosis⁴³ in lung cancer cells. Another lncRNA, SENCR can stabilize the vascular endothelial cell adherent junctions⁴⁴, which mediate cell adhesion and intracellular signals⁴⁵. Accordingly, the upregulation of SENCR may be involved in regulating the enhancement of cell adhesion seen in ZIKV-infected TEC. The transcriptome and miRNA sequencing also revealed that genes associated with adhesion and ECM molecules were activated by infection by ZIKV, including coding and non-coding genes. Overall, our findings highlight a high complexity in the intracellular circuitry controlling ZIKV replication within human TEC.

In addition to extended in vitro studies showing that human TEC is infected and respond to ZIKV, we evaluated the presence of the virus as well as the general thymus microarchitecture in histological sections obtained from children who died early in life, whose mothers were infected with ZIKV during pregnancy. Although the number of specimens was low, thymuses were widely affected, varying from massive loss of cortical thymocytes, to preservation of the cortex in some cases. Within the thymic lobules and septae, we could see an increase in ECM deposition, as revealed by histological staining of collagen-containing fibrils. Most importantly, there was an alteration in the thymic epithelial network with a clear-cut increase in the formation of Hassall’s corpuscles, denoting a higher degree of keratinization (corresponding to terminal epithelial cell differentiation). In the same vein, we observed an increase in keratin contents in ZIKV-infected TEC cultures. It is important to point out that such enhancement of Hassall’s corpuscles is also found in cases of glucocorticoid treatment of patients with primary immunodeficiency⁴⁶, as well as in mice treated with dexamethasone⁴⁷.

In 5 out of the 6 thymuses analyzed, we detected the presence of viral protein colocalized with TEC, as ascertained by immunofluorescence using the 4G2 anti-flavivirus and anti-cytokeratin antibodies. Not only we detected the viral protein in Hassall’s corpuscles, as briefly reported by Valdespino-Vázquez and co-workers⁴⁸, but also in TEC outside the corpuscles and that are located in close contact with developing thymocytes.

Rather surprisingly however, the virus genome was detected, by using qRT-PCR, only in one case. This can be explained by the fact the material applied was the bulk of thymic tissues whereas the immunohistochemical data revealed detectable amounts of the virus only in some microenvironmental cells (most of them epithelial). Moreover, in this particular case we did not see 4G2⁺ cells in the thymus sections we examined. One conceivable explanation for the fact that the only thymic extract positive for qRT-PCR was negative for immunolabeling with the 4G2 monoclonal antibody is that in this particular case, the presence of the virus was restricted to one part of the thymus that was not evaluated by immunohistochemistry.

Finally, from a conceptual point of view, it is conceivable that the human thymic epithelium is a target for ZIKV infection, with general but also specific consequences, as for example the modulation of adhesion molecules, extracellular matrix and cytokines and chemokines known to modulate intrathymic T-cell development.

Go to:

Methods and Patients

Human TEC line and Zika virus infection

The human postnatal TEC line¹⁷ was maintained in RPMI-1640 medium - supplemented with 10% fetal calf serum, 2 g/L sodium bicarbonate, 2 g/L HEPES, 1X antimycotic solution (Sigma-Aldrich) and 5 µg/mL of plasmocin prophylactic (InvivoGen) at 37 °C, in an atmosphere containing 5% CO₂.

The ZIKV strain used in this study was RIO-U1, which was isolated in 2016 from urine of a patient from Rio de Janeiro state, Brazil (GenBank Accession number: KU926309), as described elsewhere⁴⁹. Viral stocks of ZIKV Rio-U1 were prepared by infecting Vero cell monolayers. The titer (PFU/ml) and the genome integrity were determined as described elsewhere⁴⁹. Samples of the stock were used in the TEC infection assays. Vero cell cultures were also applied to test the infectivity of virus particles derived from infected TEC.

Unless described elsewhere, growing human TEC was plated in T-25 flasks (10⁵ cells) and infected 24 h later with 0.1 or 1.0 multiplicity of infection (MOI) or with 1 mL of infected TEC-derived supernatants for 1 h at 37°C in an atmosphere containing 5% CO₂. After infection, culture medium was changed and cells were maintained in culture for 72 h. Controls (MOCK group) corresponded to cell grown under the same culture conditions but not infected by ZIKV.

CEM cell line and ZIKV infection

CEM is an acute T-cell lymphoblastic leukemia cell line that was maintained in RPMI-1640 medium supplemented as described above. ZIKV Infection was performed with 1.0 MOI, as described previously for peripheral blood mononuclear cells⁵⁰.

Viral Multiplication in TEC line

An aliquot of supernatant was collected every 24 h after infection for three days and frozen at −80 °C. Briefly, Vero cells were seeded in a 24-well plate (5 × 10⁴ cells/cm²), 24 h before inoculation. Serial dilutions of the supernatant were used to infect cell monolayers. After 1 h incubation at 37 °C, the supernatant was replaced by 2.4% CMC (carboxymethyl cellulose) in Earle’s 199 complete medium supplemented with 5% FBS, followed by incubation for 7 days at 37 °C. Cells were fixed with 10% formaldehyde, washed, and dyed with 0.4% crystal violet. Viral titer was determined from the numbers of plaques visualized.

Inhibition of infection

TEC were plated and 24 h later incubated for 30 min at 37 °C in an atmosphere containing 5% CO₂ with culture media containing 10 µg/mL of goat anti-AXL polyclonal antibody (R&D Systems). Identical concentrations of IgG from normal goat serum were used as control. IgG amounts present in normal goat serum were quantified in Nanodrop™ 2000/2000c Spectrophotometers (Thermo Fisher Scientific). Cells were then infected as described previously in the presence of anti-AXL. Our experiments on the role of the AXL receptor in ZIKV entry in TECs were based on the paper of Hamel and co-workers, showing that exposure of the human skin fibroblast cell line HFF1 to ZIKV (and to Dengue virus, used as a positive control) resulted in comparable numbers of infected cells that were inhibited 70% and 50%, respectively, in the presence of the same neutralizing anti-AXL antibody⁵¹. Infection was quantified by flow cytometry as indicated below.

Antibodies

See Supplementary Table ³.

Flow cytometry

Intracellular staining of viral protein was performed with True-Nuclear™ transcriptional factor buffer set (Biolegend). Cells were incubated with 1x Fix Concentrate buffer for 40 min at 4 °C, followed by incubation with normal human serum for 20 min at 4 °C. Then, cells were stained with 4G2 antibody for 1 h at 4 °C and stained with secondary antibody for 30 min at 4 °C. After, cells were washed with 1x Perm buffer, as well as between each incubation. Subsequently, cells were fixed with formaldehyde 2% and evaluated by flow cytometry.

Cell proliferation status was performed by intracellular staining of Ki67 protein as indicated above for detection of viral protein, with the exception of that a single antibody incubation was realized for 30 min. We further evaluated those cells which are Ki67 high expressers, which correspond to those leaving cells in the mitosis phase of the cell cycle.

Cell death was analyzed by staining with Annexin-V-APC (Immunotools) and PI (Propidium iodide– Sigma-Aldrich). Cells were resuspended in 100 µL of 1x Annexin-V binding buffer (BD Biosciences) and incubated with Annexin-V-APC (1 µL) for 10 min at room temperature. PI (10 µg/mL) was added to the suspension prior acquisition on the flow cytometer.

To analyze ICAM-1 (CD54) and VCAM-1 (CD106) expression cells were incubated with normal human serum for 20 min at 4 °C and then incubated with CD54-PECY5 or CD106-APC antibodies for 30 min at 4 °C. After incubations cells were washed with PBS. Then, the cells were fixed with formaldehyde 2% and evaluated by flow cytometry.

The fluorescence intensity of labeled TEC was determined with FACSCanto II flow cytometer (BD Biosciences) Analyses were carried out using the BD FACSDiva 6.1.3 software (BD Biosciences).

Immunofluorescence

TEC were seeded in 8 well (10³ cells/well) Lab-Tek^® Chamber Slide™ System and infected as described previously. After infection, the slides were washed twice with PBS and fixed in 4% paraformaldehyde for 15 min. Slides were washed again and treated with permeabilization buffer (0.1% saponin and 1% BSA – bovine albumin serum) for 15 min at 37 °C. Thereafter, samples were incubated for 40 min with blocking solution (0.2% normal goat serum, 1% FBS, 1% BSA, 0.3% X-100 triton) and submitted to staining with primary antibody for 1 h and with secondary antibody for 30 min. After staining, the slides were washed twice with permeabilization buffer, stained with DAPI for 10 min, washed and mounted with Dako Fluorescent Mounting Medium (DAKO).

Immunostained samples were analyzed by an AxioImager A2 device using the AxioVision Rel 4.8 software (Zeiss) or by Leica TSC SP8 Confocal device using LAS-X Software (Leica Microsystems). Negative controls, in which the secondary antibody was used alone, did not generate any significant labeling. The percentage of infection was achieved by the number of 4G2⁺ cells divided by the total number of cells per field. We analyzed 5 fields from 3 independent experiments and each experiment was performed in triplicate. The quantitative fluorescence analysis of cytokeratin expression was performed by transforming the specific staining into an eight-bit grey image for analysis of the corrected total cell fluorescence (CTCF) through the ImageJ software (NIH). At least 50 cells per condition were analyzed CTCF is calculated by the integrated density minus the area of a given cell multiplies by the mean fluorescence of background readings. When necessary, brightness and contrast adjustments were performed in the entire images.

TEC adhesion to lymphoblastic T cells

Human TEC were seeded (10⁴ or 5 × 10⁴ cells) and 24 h later infected as described previously. After 72 h in culture, CEM cells (human T-cell acute lymphoblastic leukemia cell line) were added to the infected culture in a proportion of 50 CEM cells per TEC. We led these cells to adhere onto ZIKV-infected or non-infected-TEC during 1 h in serum-free RPMI medium at 37 °C in an atmosphere containing 5% CO₂. Non-adherent CEM cells were gentle washed out with PBS at 37 °C. Culture flasks were fixed with methanol for 7 min and stained with Giemsa. The number of adherent CEM cells per at least 500 TEC were counted and the association index (AI) was calculated by the following formula, previously validated for this type of analysis⁵²,⁵³.

AI=numberofTECwithCEMcells×numberofCEMcellsboundtoTECtotalTECnumber×100

Thymic epithelial and lymphoblastic T cell co-cultures

TEC were seeded in T-25 flasks (5 × 10⁴ cells) and 24 h later infected as described previously. Seventy-two hours post-infection, the medium was removed and CEM cells were added to the culture in a proportion of 25 CEM cells per TEC with fresh RPMI medium supplemented with 15% FBS. After 24 h of co-culture, non-adherent cells were removed by washing the flasks with PBS. Remaining adherent cells were recovered after addition of trypsin to the flasks. Both adherent and non-adherent cells were stained with 4G2 and anti-cytokeratin antibody by flow cytometry, as described previously. Anti-cytokeratin antibody was used for distinguishing CEM cells and TEC in adherent and non-adherent cells suspensions.

RNA extraction

TEC were seeded in T-25 flasks (10⁵ cells) and 24 h later infected as described previously. Seventy-two hours post-infection, RNA extraction was performed with MOCK and ZIKV-infected TEC using RNeasy Plus Mini Kit (Qiagen, Hilden, Germany), following the manufacturer’s instructions for mRNA and small RNA extraction. Total RNA concentration and purity were determined by Nanodrop spectrophotometer. RNA integrity was verified through Agilent 2100 Bioanalyzer (Agilent Technologies).

Transcriptome and small RNA library preparation and sequencing

The 12 transcriptome libraries were prepared from 500 ng of total RNA using the TruSeq Stranded Total RNA Sample Preparation Kit with Ribo-Zero Gold (Illumina). Additionally, the 12 small RNA libraries were prepared from 1 μg from the same total RNA of transcriptome analysis, using the TruSeq Small RNA Library Prep Kit (Illumina). All the steps of the protocols were carried out without any modifications, with the exception of the cDNA construction purification step on the small RNA protocol, which was performed using 3% Agarose Gel Cassette for targets between 100 bp − 250 bp in a BluePippin system (Sage Science). The quality control of each library was performed using the 2100 Bioanalyzer System with the Agilent High Sensitivity DNA Kit (Agilent). The libraries were quantified via qPCR using a KAPA Library Quantification Kits for Illumina platforms (KAPA Biosystems). The Libraries were sequenced on a NextSeq. 500 sequencing system (Illumina), for the RNA-seq paired-end reads (2 × 75 bp) were obtained using NextSeq. 500/550 High Output v2 kit (150 cycles) and for the small RNA, single-end reads (1 × 50 bp) were obtained using a NextSeq. 500/550 High Output v2 kit (75 cycles).

Transcriptome analysis

The 12 transcriptome and 12 small RNA libraries were quality checked with FASTQC tool (https://www.bioinformatics.babraham.ac.uk/projects/fastqc/) and further trimmed with Trimmomatic tool⁵⁴. For small RNAse libraries, only reads with length between 16-28 nucleotides were kept for posterior analysis. The high-quality transcriptome libraries were mapped separately on the Human genome (GRCH38.p12) from ENSEMBL databank (http://www.ensembl.org/info/data/ftp/index.html) using the STAR mapper (Version 2.6)⁵⁵, with default parameters. The small RNA reads were mapped in the human microRNA precursors retrieved from miRBase (Release 22) (http://www.mirbase.org/) with Bowtie2 tool (version 2.2.1)⁵⁶ running with default parameters.

The Differential Gene Expression (DGE) analysis was performed with DESeq. 2 package⁵⁷ from R. The genes (protein-coding, lncRNAs, circRNAs, and other transcripts) and microRNAs with adjusted p-value < 0.01 and Log2 Fold Change > |1.0| were considered for further analyses. Both transcriptome and microRNA expression analysis was carried out in six biological replicates.

In order to identify microRNAs and their real targets, we retrieved the experimentally validated microRNA-target interactions from miRTarBase (http://mirtarbase.mbc.nctu.edu.tw/php/index.php) plus the predicted ones from SpidermiR package⁵⁸, which includes DIANA⁵⁹, miRanda⁶⁰, PicTar⁶¹ and TargetScan⁶² microRNA target datasets. We compared with differentially expressed microRNAs and target genes from both transcriptome and small RNA seq experiments.

Circular RNA (circRNAs) Analysis

Using the 12 transcriptome libraries, we performed a DGE analysis with DESEQ. 2 library⁵⁷ of the circular RNAs (circRNAs) from human genome, based on circBase (http://www.circbase.org/). The circRNAs with adjusted p-value lower than 0.01 and Log₂ Fold Change > |1.0| were used to characterize the interaction of microRNAs and circRNAs with miRanda tool⁶⁰, using strict 5′ seed pairing, minimum score threshold, alignment length and identity of 140, 17 and 80%, respectively. The interactions were visualized with Cytoscape tool⁶³.

Enrichment analysis of DEGs

The Gene Ontology (GO)⁶⁴ and ReactomePA⁶⁵ were used to determine enrichment of biological processes and pathways, respectively on DEGs. We used the ClusterProfiler package from R⁶⁶ for ReactomePA generation and GO enrichment by over-representation test. The enrichment was generated with Benjamini & Hochberg (BH) adjusted p-value method with a p-value cutoff of 0.01.

Human thymus histology and immunostainings

We evaluated human thymus fragments from children who died up to 6 months after birth, and whose mothers were infected with ZIKV during pregnancy. The clinical features of these newborns are summarized in Supplemental Table ⁴. The use of organ histological sections for research has been approved by the Fernandes Figueira National Institute for the Health of Mother, Child, and Adolescent, Oswaldo Cruz Foundation review board under the number CAAE: 52675616.0.000.5269. All participants provided written informed consent. All methods were performed in accordance with the relevant guidelines and regulations. General histological profiles were evaluated through haematoxin-eosin staining. Gomori’s reticulin and Masson’s trichrome staining procedures were applied for revealing extracellular matrix containing compartments.

In brief, paraffin-embedded thymus sections (4 µm) were maintained at 65 °C overnight. The slides were deparaffinized in two passages in xylene and three passages in absolute alcohol for two min each. Samples were hydrated in distilled water and the peroxidase blockade was performed. The slides were washed with PBS and antigenic recovery was performed in water bath at 96 °C with 10% Trilogy solution for 40 min. After cooling for 20 min, the slides were washed with PBS and blocked with 2% BSA for 40 min. Thereafter, samples were incubated overnight at 4 °C with anti-cytokeratin and 4G2 antibody. The slides were washed twice with permeabilization buffer and incubated for 30 min with secondary antibodies for 30 min. After being washed twice with permeabilization buffer, samples were incubated for 10 min with DAPI and washed again. Immunostained samples were analyzed by the Leica TSC SP8 Confocal device using LAS-X Software.

Two thymuses derived from 2 infants were applied for controls. Although in both cases there was a pathological situation associated (pulmonary hypoplasia), their normal histological profile was confirmed by hematoxylin-eosin staining. Sections of these specimens were subjected to the same histological and immunofluorescence procedures.

Statistical analysis

Unless described elsewhere, the statistical significance was approached using unpaired one- or two-tailed t-test or One-way ANOVA followed by Tukey’s multiple comparison test.

Go to:

Supplementary information

Supplemental material.^{(1.0M, docx)}

Supplementary Table 1.^{(11M, pdf)}

Supplementary Table 2.^{(808K, pdf)}

Go to:

Acknowledgements

The human postnatal TEC line was kindly provided by Dr. Maria Luiza Toríbio (Universidad Autonoma de Madrid, Madrid, Spain). The mouse monoclonal 4G2 antibody (anti-flavivirus) was gently provided by the Laboratory of Monoclonal Antibody Technology from Bio-Manguinhos (Oswaldo Cruz Foundation, Rio de Janeiro, Brazil). The authors thank Raquel de Souza Martins for confocal analysis assistance. This work was funded by Brazilian Funding Agency for Studies and Projects (FINEP) #01.16.0078.00; Brazilian National Council for Scientific and Technological Development (CNPq) ATRV supported by #303170/2017-4, WS by #306699/2018-4, DAMC by #307720/2016-0; Carlos Chagas Filho Foundation for Funding Research in the State of Rio de Janeiro (Faperj) ATRV supported by #E-26/202.826/2018; the Oswaldo Cruz Foundation (Fiocruz) and National Laboratory for Scientific Computation; Coordination of Superior Level Staff Improvement (CAPES) - Brazil - Finance Code 001 and GLM, JBC and LWPA supported #23038.010041/2013-13; Mercosur Fund for Structural Convergence (FOCEM, Mercosur #03/11) and the Program of Preventing and Combating the Zika Virus, MCTIC/FNDCT - CNPq/MEC-CAPES/MS-Decit (Grants 426767/2018-7 and 88881.130684/2016-01 for MCB). This work was developed in the framework of the Brazilian National Institute of Science and Technology on Neuroimmunomodulation (INCT-NIM) WS supported by #465489/2014-1.

Go to:

Author contributions

C.V.M., G.L.M. and J.B.C.: conceived and designed the experiments, performed the experiments, analyzed data and wrote the manuscript. M.N.G., D.P.C., D.A.F.O., A.L.G. and L.M.S.R.: conceived and designed the experiments, performed the experiments, analyzed data. Z.V., E.A.P., M.E.L.M., M.C.B., I.R. and V.C.A.: conceived and designed the experiments; L.W.P.A.: analyzed data; N.F: performed the experiments; A.T.R.V., D.A.M.C. and W.S.: conceived and designed the experiments, wrote the manuscript and coordinated the research.

Go to:

Data availability

The RNA-seq and small-RNA-seq data that support the findings of this study have been deposited in SRA database from NCBI with the accession codes PRJNA560342 and PRJNA561078. All other data are available in the article and its supplementary information files or from the corresponding authors upon request.

Go to:

Code availability

We use already published algorithms. See Methods section for full description of each analysis including input data, library, and algorithm version used, for which we provide corresponding web links or references.

Go to:

Competing interests

The authors declare no competing interests.

Go to:

Footnotes

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

These authors contributed equally: Carolina V. Messias, Guilherme Loss Morais and Joseane Biso de Carvalho.

These authors jointly supervised this work: Ana T. R. Vasconcelos, Daniella A. Mendes-da-Cruz and Wilson Savino.

Go to:

Contributor Information

Ana T. R. Vasconcelos, Email: rb.ccnl@vrta.

Daniella A. Mendes-da-Cruz, Email: rb.zurcoif.coi@alleinad.

Wilson Savino, Email: rb.zurcoif@onivas.nosliw.

Go to:

Supplementary information

is available for this paper at 10.1038/s41598-020-58135-y.

Go to:

References

1. Martines RB, et al. Pathology of congenital Zika syndrome in Brazil: a case series. Lancet. 2016;388:898–904. doi: 10.1016/S0140-6736(16)30883-2. [PubMed] [CrossRef] [Google Scholar]
2. Carod-Artal FJ. Neurological complications of Zika virus infection. Expert Rev. Anti. Infect. Ther. 2018;16:399–410. doi: 10.1080/14787210.2018.1466702. [PubMed] [CrossRef] [Google Scholar]
3. Lopes Moreira ME, et al. Neurodevelopment in Infants Exposed to Zika Virus In Utero. N. Engl. J. Med. 2018;379:2377–2379. doi: 10.1056/NEJMc1800098. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
4. Nielsen-Saines K, et al. Delayed childhood neurodevelopment and neurosensory alterations in the second year of life in a prospective cohort of ZIKV-exposed children. Nat. Med. 2019;25:1213–1217. doi: 10.1038/s41591-019-0496-1. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
5. Da Silva IRF, Frontera JA, De Filippis AMB, Do Nascimento OJM. Neurologic complications associated with the Zika virus in Brazilian adults. JAMA Neurol. 2017;74:1190–1198. doi: 10.1001/jamaneurol.2017.1703. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
6. Savino W, et al. Zika Virus Infection in the Elderly: Possible Relationship with Guillain-Barré Syndrome. Gerontology. 2017;63:210–215. doi: 10.1159/000453579. [PubMed] [CrossRef] [Google Scholar]
7. Garcez PP, et al. Zika virus impairs growth in human neurospheres and brain organoids. Science (80-.). 2016;352:816–818. doi: 10.1126/science.aaf6116. [PubMed] [CrossRef] [Google Scholar]
8. Azouz F, Arora K, Krause K, Nerurkar VR, Kumar M. Integrated MicroRNA and mRNA profiling in zika virus-infected neurons. Viruses. 2019;11:1–22. doi: 10.3390/v11020162. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
9. Aid M, et al. Zika Virus Persistence in the Central Nervous System and Lymph Nodes of Rhesus Monkeys. Cell. 2017;169:610–620. doi: 10.1016/j.cell.2017.04.008. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
10. Hirsch AJ, et al. Zika Virus infection of rhesus macaques leads to viral persistence in multiple tissues. PLoS Pathog. 2017;13:e1006219. doi: 10.1371/journal.ppat.1006219. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
11. Savino W, Mendes-da-Cruz DA, Lepletier A, Dardenne M. Hormonal control of T-cell development in health and disease. Nat. Rev. Endocrinol. 2016;12:77–89. doi: 10.1038/nrendo.2015.168. [PubMed] [CrossRef] [Google Scholar]
12. Braun J, et al. Productive and Persistent Infection of Human Thymic Epithelial Cells in Vitrowith HIV-1. Virology. 1996;225:413–418. doi: 10.1006/viro.1996.0617. [PubMed] [CrossRef] [Google Scholar]
13. Maguer-Satta V, Gazzolo L, Dodon MD. Human Immature Thymocytes as Target Cells of the Leukemogenic Activity of Human T-cell Leukemia Virus Type I. Blood. 1995;86:1444–1452. doi: 10.1182/blood.V86.4.1444.bloodjournal8641444. [PubMed] [CrossRef] [Google Scholar]
14. Manca N, et al. Detection of HTLV-I tax-rex and pol gene sequences of thymus gland in a large group of patients with myasthenia gravis. J. Acquir. immune Defic. Syndr. 2002;29:300–306. doi: 10.1097/00126334-200203010-00012. [PubMed] [CrossRef] [Google Scholar]
15. Moreira-Ramos K, Castro FMM, De, Linhares-Lacerda L, Savino W. Can thymic epithelial cells be infected by human T-lymphotropic virus type 1? Mem. Inst. Oswaldo Cruz. 2011;106:759–762. doi: 10.1590/S0074-02762011000600018. [PubMed] [CrossRef] [Google Scholar]
16. Carvalho Barros LR, et al. HTLV-1-infected thymic epithelial cells convey the virus to CD4 + T lymphocytes. Immunobiology. 2017;222:1053–1063. doi: 10.1016/j.imbio.2017.08.001. [PubMed] [CrossRef] [Google Scholar]
17. Fernández E, et al. Establishment and characterization of cloned human thymic epithelial cell lines. Analysis of adhesion molecule expression and cytokine production. Blood. 1994;83:3245–3254. [PubMed] [Google Scholar]
18. Ribeiro-Carvalho MM, et al. Triiodothyronine modulates thymocyte migration. Scand. J. Immunol. 2007;66:17–25. doi: 10.1111/j.1365-3083.2007.01928.x. [PubMed] [CrossRef] [Google Scholar]
19. Henchal E, Gentry MK, CCown JM, Brandt WE. Dengue Virus-Specific and Flavivirus Group Determinants Identified with Monoclonal Antibodies by Indirect Immunofluorescence. Am. J. Trop. Med. Hyg. 1982;31:830–836. doi: 10.4269/ajtmh.1982.31.830. [PubMed] [CrossRef] [Google Scholar]
20. Paessens LC, Singh SK, Fernandes RJ, van Kooyk Y. Vascular cell adhesion molecule-1 (VCAM-1) and intercellular adhesion molecule-1 (ICAM-1) provide co-stimulation in positive selection along with survival of selected thymocytes. Mol. Immunol. 2008;45:42–8. doi: 10.1016/j.molimm.2007.05.016. [PubMed] [CrossRef] [Google Scholar]
21. Foo SS, et al. Asian Zika virus strains target CD14+blood monocytes and induce M2-skewed immunosuppression during pregnancy. Nat. Microbiol. 2017;2:1558–1570. doi: 10.1038/s41564-017-0016-3. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
22. Michlmayr D, Andrade P, Gonzalez K, Balmaseda A, Harris E. CD14+CD16+monocytes are the main target of Zika virus infection in peripheral blood mononuclear cells in a paediatric study in Nicaragua. Nat. Microbiol. 2017;2:1462–1470. doi: 10.1038/s41564-017-0035-0. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
23. Kannan RP, Hensley LL, Evers LE, Lemon SM, McGivern DR. Hepatitis C Virus Infection Causes Cell Cycle Arrest at the Level of Initiation of Mitosis. J. Virol. 2011;85:7989–8001. doi: 10.1128/JVI.00280-11. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
24. Strange DP, Green R, Siemann DN, Gale M, Verma S. Immunoprofiles of human Sertoli cells infected with Zika virus reveals unique insights into host-pathogen crosstalk. Sci. Rep. 2018;8:1–15. doi: 10.1038/s41598-018-27027-7. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
25. Limonta D, et al. Human fetal astrocytes infected with zika virus exhibit delayed apoptosis and resistance to interferon: Implications for persistence. Viruses. 2018;10:1–19. doi: 10.3390/v10110646. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
26. Singh PK, et al. Determination of system level alterations in host transcriptome due to Zika virus (ZIKV) Infection in retinal pigment epithelium. Sci. Rep. 2018;8:1–16. doi: 10.1038/s41598-017-17765-5. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
27. Liu Ling, Chen Zhenyu, Zhang Xin, Li Shun, Hui Yi, Feng Hexi, Du Yanhua, Jin Guohua, Zhou Xiaohui, Zhang Xiaoqing. Protection of ZIKV infection-induced neuropathy by abrogation of acute antiviral response in human neural progenitors. Cell Death & Differentiation. 2019;26(12):2607–2621. doi: 10.1038/s41418-019-0324-7. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
28. Bertoni A, Alabiso O, Galetto AS, Baldanzi G. Integrins in T cell physiology. Int. J. Mol. Sci. 2018;19:E485. doi: 10.3390/ijms19020485. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
29. Lannes-Vieira J, et al. Extracellular matrix components of the mouse thymic microenvironment. III. Thymic epithelial cells express the VLA6 complex that is involved in laminin-mediated interactions with thymocytes. Int. Immunol. 1993;5:1421–1430. doi: 10.1093/intimm/5.11.1421. [PubMed] [CrossRef] [Google Scholar]
30. Le PT, Singer KH. Human thymic epithelial cells: adhesion molecules and cytokine production. Int. J. Clin. Lab. Res. 1993;23:56–60. doi: 10.1007/BF02592284. [PubMed] [CrossRef] [Google Scholar]
31. Ocampo JSP, et al. Laminin-211 controls thymocyte-thymic epithelial cell interactions. Cell. Immunol. 2008;254:1–9. doi: 10.1016/j.cellimm.2008.06.005. [PubMed] [CrossRef] [Google Scholar]
32. Savino W, Mendes-da-Cruz DA, Golbert DCF, Riederer I, Cotta-de-Almeida V. Laminin-Mediated Interactions in Thymocyte Migration and Development. Front. Immunol. 2015;6:579. doi: 10.3389/fimmu.2015.00579. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
33. Dairaghi DJ, et al. Macrophage Inflammatory Protein-1 β Induces Migration and Activation of Human Thymocytes. Blood. 1998;91:2905–2913. doi: 10.1182/blood.V91.8.2905.2905_2905_2913. [PubMed] [CrossRef] [Google Scholar]
34. Savino W, Mendes-da-Cruz DA, Silva JS, Dardenne M, Cotta-de-Almeida V. Intrathymic T-cell migration: a combinatorial interplay of extracellular matrix and chemokines? Trends Immunol. 2002;23:305–313. doi: 10.1016/S1471-4906(02)02224-X. [PubMed] [CrossRef] [Google Scholar]
35. Zamarchi R, et al. Expression and functional activity of CXCR-4 and CCR-5 chemokine receptors in human thymocytes. Clin. Exp. Immunol. 2002;127:321–330. doi: 10.1046/j.1365-2249.2002.01775.x. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
36. Bunting MD, Comerford I, Kara EE, Korner H, McColl SR. CCR6 supports migration and differentiation of a subset of DN1 early thymocyte progenitors but is not required for thymic nTreg development. Immunol. Cell Biol. 2014;92:489–498. doi: 10.1038/icb.2014.14. [PubMed] [CrossRef] [Google Scholar]
37. Singer KH, Denning SM, Whichard LP, Haynes BF. Thymocyte LFA-1 and thymic epithelial cell ICAM-1 molecules mediate binding of activated human thymocytes to thymic epithelial cells. J. Immunol. 1990;144:2931–9. [PubMed] [Google Scholar]
38. Wang M, et al. Circular RNAs: A novel type of non-coding RNA and their potential implications in antiviral immunity. Int. J. Biol. Sci. 2017;13:1497–1506. doi: 10.7150/ijbs.22531. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
39. Cadena C, Hur S. Antiviral Immunity and Circular RNA: No End in Sight. Mol. Cell. 2017;67:163–164. doi: 10.1016/j.molcel.2017.07.005. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
40. Nishitsuji H, et al. Long noncoding RNA #32 contributes to antiviral responses by controlling interferon-stimulated gene expression. Proc. Natl. Acad. Sci. 2016;113:10388–10393. doi: 10.1073/pnas.1525022113. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
41. Carnero E, et al. Long noncoding RNA EGOT negatively affects the antiviral response and favors HCV replication. EMBO Rep. 2016;17:1013–1028. doi: 10.15252/embr.201541763. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
42. Zhao JM, et al. Long non-coding RNA PICART1 suppresses proliferation and promotes apoptosis in lung cancer cells by inhibiting JAK2/STAT3 signaling. Neoplasma. 2018;65:779–789. doi: 10.4149/neo_2018_171130N778. [PubMed] [CrossRef] [Google Scholar]
43. Zhang C, et al. LncRNA PICART1 suppressed non-small cell lung cancer cells proliferation and invasion by targeting AKT1 signaling pathway. Am. J. Clin. Res. 2018;10:4193–4201. [PMC free article] [PubMed] [Google Scholar]
44. Lyu Q, et al. SENCR stabilizes vascular endothelial cell adherens junctions through interaction with CKAP4. Proc. Natl. Acad. Sci. USA. 2019;116:546–555. doi: 10.1073/pnas.1810729116. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
45. Dejana E, Orsenigo F. Endothelial adherens junctions at a glance. J. Cell Sci. 2013;126:2545–2549. doi: 10.1242/jcs.124529. [PubMed] [CrossRef] [Google Scholar]
46. Hale LP, Markert ML. Corticosteroids Regulate Epithelial Cell Differentiation and Hassall Body Formation in the Human Thymus. J. Immunol. 2004;172:617–624. doi: 10.4049/jimmunol.172.1.617. [PubMed] [CrossRef] [Google Scholar]
47. Savino W, et al. Hydrocortisone increases the numbers of KL1+ cells, a discrete thymic epithelial cell subset characterized by high molecular weight cytokeratin expression. Endocrinology. 1988;123:2557–64. doi: 10.1210/endo-123-5-2557. [PubMed] [CrossRef] [Google Scholar]
48. Valdespino-Vázquez MY, et al. Congenital Zika Syndrome and Extra-Central Nervous System Detection of Zika Virus in a Pre-term Newborn in Mexico. Clin. Infect. Dis. 2019;68:903–912. doi: 10.1093/cid/ciy616. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
49. Bonaldo MC, et al. Isolation of Infective Zika Virus from Urine and Saliva of Patients in Brazil. PLoS Negl. Trop. Dis. 2016;10:e0004816. doi: 10.1371/journal.pntd.0004816. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
50. Messias CV, et al. Zika virus infects human blood mononuclear cells. BMC Infect. Dis. 2019;19:986. doi: 10.1186/s12879-019-4622-y. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
51. Hamel R, et al. Biology of Zika Virus Infection in Human Skin Cells. J. Virol. 2015;89:8880–8896. doi: 10.1128/JVI.00354-15. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
52. Ribeiro-Carvalho MM, Farias-de-Oliveira DA, Villa-Verde DMS, Savino W. Triiodothyronine modulates extracellular matrix-mediated interactions between thymocytes and thymic microenvironmental cells. Neuroimmunomodulation. 2002;10:142–152. doi: 10.1159/000067175. [PubMed] [CrossRef] [Google Scholar]
53. Linhares-Lacerda L, et al. RNA interference-mediated knockdown of CD49e (α5 integrin chain) in human thymic epithelial cells modulates the expression of multiple genes and decreases thymocyte adhesion. BMC Genomics. 2010;11:S2. doi: 10.1186/1471-2164-11-S5-S2. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
54. Bolger AM, Lohse M, Usadel B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–2120. doi: 10.1093/bioinformatics/btu170. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
55. Dobin A, et al. STAR: Ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29:15–21. doi: 10.1093/bioinformatics/bts635. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
56. Langmead B, Slazberg SL. Fast gapped-read alignmnet with Bowtie 2. Nat. Methods. 2012;9:357–359. doi: 10.1038/nmeth.1923. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
57. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq. 2. Genome Biol. 2014;15:550. doi: 10.1186/s13059-014-0550-8. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
58. Cava C, et al. SpidermiR: An R/bioconductor package for integrative analysis with miRNA data. Int. J. Mol. Sci. 2017;18:E274. doi: 10.3390/ijms18020274. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
59. Paraskevopoulou MD, et al. DIANA-microT web server v5.0: service integration into miRNA functional analysis workflows. Nucleic Acids Res. 2013;41:169–173. doi: 10.1093/nar/gkt393. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
60. John B, et al. Human microRNA targets. PLoS Biol. 2004;2:e363. doi: 10.1371/journal.pbio.0020363. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
61. Krek A, et al. Combinatorial microRNA target predictions. Nat. Genet. 2005;37:495–500. doi: 10.1038/ng1536. [PubMed] [CrossRef] [Google Scholar]
62. Agarwal V, Bell GW, Nam J-W, Bartel DP. Predicting effective microRNA target sites in mammalian mRNAs. Elife. 2015;4:1–38. doi: 10.7554/eLife.05005. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
63. Shannon P, et al. Cytoscape: A Software Environment for Integrated Models of Biomolecular Interaction Networks. Genome Res. 2003;13:2498–2504. doi: 10.1101/gr.1239303. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
64. Ashburner M, et al. Gene Ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat. Genet. 2000;25:25–29. doi: 10.1038/75556. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
65. Fabregat A, et al. The Reactome Pathway Knowledgebase. Nucleic Acids Res. 2018;46:D649–D655. doi: 10.1093/nar/gkx1132. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
66. Yu G, Wang L-G, Han Y, He Q-Y. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16:284–287. doi: 10.1089/omi.2011.0118. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

Articles from Scientific Reports are provided here courtesy of Nature Publishing Group