Janaina Macedo-da-Silva1,2†
Lívia Rosa-Fernandes1*†
Raquel Hora Barbosa3
Claudia B. Angeli1
Fabiana Rabe Carvalho3,4
Renata Artimos de Oliveira Vianna3,4
Paulo C. Carvalho2Martin R. Larsen5
Claudete Araújo Cardoso3,4
Giuseppe Palmisano1*- 1GlycoProteomics Laboratory, Department of Parasitology, Institute of Biomedical Sciences, University of Sao Paulo, Sao Paulo, Brazil
- 2Laboratory for Structural and Computational Proteomics, Carlos Chagas Institute, Fiocruz, Curitiba, Brazil
- 3Maternal and Child Department, School of Medicine, Universidade Federal Fluminense, Niteroi, Brazil
- 4Multiuser Laboratory for Research in Nephrology and Medical Sciences (LAMAP), School of Medicine, Universidade Federal Fluminense, Niteroi, Brazil
- 5Department of Biochemistry and Molecular Biology, University of Southern Denmark, Odense, Denmark
In 2015, ZIKV infection attracted international attention during an epidemic in the Americas, when neurological disorders were reported in infants who had their mothers exposed to ZIKV during pregnancy. World Health Organization (WHO) epidemiological data show that 5 to 15% of neonates exposed to ZIKV in the uterus have complications included in abnormalities related to Congenital Zika Syndrome (CZS). The risk of complications after birth is not well documented, however, clinical evidence shows that 6% of infants exposed to ZIKV during pregnancy have complications present at birth, and this rate rises to 14% when medical monitoring is performed in all exposed infants, regardless of birth condition. Thus, the evaluation and monitoring of all exposed infants are of foremost importance as the development of late complications has been increasingly supported by clinical evidence. The identification of changes in protein profile of infants exposed to ZIKV without CZS could provide valuable findings to better understand molecular changes in this cohort. Here, we use a shotgun-proteomics approach to investigate alterations in the serum of infants without CZS symptoms but exposed to intrauterine ZIKV (ZIKV) compared to unexposed controls (CTRL). A complex pattern of differentially expressed proteins was identified, highlighting the dysregulation of proteins involved in axon orientation, visual phototransduction, and global protease activity in children exposed to ZIKV without CZS. These data support the importance of monitoring children exposed to ZIKV during gestation and without early CZS symptoms. Our study is the first to assess molecular evidence of possible late disorders in children victims of the ZIKV outbreak in the Americas. We emphasize the importance of medical monitoring of symptomatic and asymptomatic children, as apparently unexplained late neurological and eye disorders may be due to intrauterine ZIKV exposure.
Introduction
Zika is a single-stranded RNA virus belonging to the Flaviviridae family that was first isolated in 1947 from a rhesus monkey in Kampala, Uganda (Dick et al., 1952). The vectors of these viruses are infected mosquitoes, also responsible for the spread of important diseases such as dengue fever, West Nile fever, and yellow fever (Noorbakhsh et al., 2019). Between 1960 and 1980, Zika virus (ZIKV) was associated with mild diseases in Asian and African populations, however, in 2013-2014, during an outbreak in French Polynesia, the possibility of transplacental transmission from mother to fetus has been described (Kindhauser et al., 2016; Noorbakhsh et al., 2019). In 2015, Zika virus infection attracted international attention during an epidemic in the Americas (Venancio et al., 2019), due to its association with increasing cases of microcephaly, congenital malformation and other neurological disorders in newborns who had their mothers infected during pregnancy (Kindhauser et al., 2016; Mlakar et al., 2016). In 2016, WHO declared the infection a public health emergency of international interest (PHEIC) (Lowe et al., 2018), including 11,546 exposed pregnant women between 2016 and 2017 only in Brazil (Heukelbach et al., 2016). In addition, data published by the Brazilian Ministry of Health indicate 14,558 suspected cases of congenital microcephaly and others central nervous system (CNS) malformations between 2015 and 2017 (Pan American Health Organization/World Health Organization, 2017).
The set of abnormalities resulting from ZIKV intrauterine infection is called congenital Zika syndrome (CZS) and is characterized by severe neurological damage and loss of intracranial volume (Wheeler, 2018). Clinical evidences of CZS are divided into structural components, such as changes in cranial morphology, and functional components that include neurological impairment (Moore et al., 2017; Mohr et al., 2018). The latest published epidemiological bulletin by the Brazilian Ministry of Health reported 19,000 suspected cases of CZS from November 2015 to May 2020, being 3,534 (18.6%) confirmed and 2,784 cases under investigation (Pan American Health Organization/World Health Organization, 2017). CZS has symptoms that are common to other congenital infections, however, has characteristics not previously seen, including severe microcephaly, macular scarring, congenital contractures, and hypertonia (Moore et al., 2017). In 2016, the suggested Zika virus seroprevalence peak in Brazil was 63% (Netto et al., 2017), and following reports indicate that ZIKV continues to circulate in the cycle of human transmission in Brazil and the Americas (Lowe et al., 2018).
The risk of complications after birth is not well established, however, clinical evidence show that 6% of infants exposed to Zika virus during pregnancy have complications at birth, and this rate rises to 14% when medical monitoring is performed on all exposed infants, regardless of birth condition (Rice et al., 2018; Musso et al., 2019; Souza et al., 2019). In addition, subcortical calcifications and an enlarged ventricle have been demonstrated in exposed children, but who were born without clinical evidence of CZS. Another interesting clinical data shows that 60% of these children have seizures during development (Rice et al., 2018; Souza et al., 2019). Reports of other congenital infections, such as cytomegalovirus and toxoplasmosis, indicate that hearing loss and eye damage can occur from 33 to 44 months after birth for symptomatic and asymptomatic infants, respectively (Dahle et al., 2000). Thus, stressing the importance of evaluating and monitoring not only symptomatic infants, but all exposed ones as the development of late complications has been increasingly supported by reports of clinical evidence (Pomar et al., 2018; Rice et al., 2018; Soares et al., 2019; Valdes et al., 2019).
Upon in utero zika virus exposure, 12.5% neonates who tested negative for the infection still presented severe, moderate, or mild complications, such as jaundice, hypotonia, hypertonia, hepatomegaly, and elevated liver enzymes (Pomar et al., 2018). Moreover, zika virus exposure without evidence of microcephaly has also been connected to events of asphyxia, hypersalivation, and reflux. In addition, exposed children had lower weight, length, and fat-free masses in the first three months of life compared to unexposed children (Soares et al., 2019). Zika virus-exposed infants without CZS also showed lower receptive language scores in the first year of life (Valdes et al., 2019), besides multiple cerebral and visual abnormalities, which would not have been identified without prolonged medical monitoring (Rice et al., 2018).
In view of the previously reported clinical data, the identification of molecular changes in infants exposed to ZIKV, but who did not have CZS, is important to understand molecular alterations underlying the occurrence of late abnormalities. Here, we use a shotgun-proteomic approach to investigate molecular markers in the serum of infants exposed to intrauterine ZIKV (ZIKV), but without symptoms of CZS, compared to those not exposed to ZIKV (CTRL). A complex pattern of differentially abundant proteins has been identified, highlighting the deregulation of proteins involved in axon guidance, visual phototransduction and global protease activity in children exposed to ZIKV without CZS. These data support the importance of monitoring children exposed to ZIKV during pregnancy and without early symptoms of CZS.
Materials and Methods
Patient Cohort
This study includes 20 infants aged between 3 and 23 months referred to the Pediatrics Service of the Antonio Pedro University Hospital, Universidade Federal Fluminense, Brazil. The cohort was divided in CTRL group (non-exposed to ZIKV and negative maternal qPCR, n = 10) and ZIKV group (positive maternal qPCR, n = 10), which consisted of patients with maternal ZIKV exposure during pregnancy and no clinical evidence of CZS. In both CTRL and ZIKV conditions, the mothers tested negative to other infectious agents (syphilis, toxoplasmosis, rubella, cytomegalovirus, and HIV). This study was approved by the institutional review board and the ethics committee of the Universidade Federal Fluminense (protocol CAAE number 79890517.6.0000.5243) and followed the guidelines of the Declaration of Helsinki. All samples were collected upon informed and written consent from the parents/legal guardians of each participant. All participants were clinically evaluated by a multidisciplinary team and are included in a currently ongoing clinical follow-up program (Vianna et al., 2019). The clinical diagnosis was performed based on the guidelines of the Ministry of Health (Brazil, Ministry of Health, 2017).
Sample Collection
Venous blood samples were collected in Vacutainer blood collection tubes with clot activator (Becton Dickson, USA) and centrifuged (1,210g for 15 min) to obtain the serum. Subsequently, the samples were aliquoted in sterile tubes and frozen at -80°C until further analysis. Clinical data were obtained during the outpatient clinic visit and from the patients’ medical records. Epidemiological and demographic data were retrieved from the questionnaire for investigating suspected cases of microcephaly related to ZIKV infection, made available by the Brazilian Ministry of Health, and applied to all patients.
Sample Preparation for Mass-Spectrometry Based-Proteomics
Serum samples were depleted using the Multiple Affinity Removal System Spin Depletion Cartridge (Agilent Technologies) as per the manufacturer’s instructions. This approach reduces the levels of the 14 most abundant serum proteins (albumin, IgG, antitrypsin, IgA, transferrin, haptoglobin, fibrinogen, alpha2-macroglobulin, alpha1-acid glycoprotein, IgM, apolipoprotein AI, apolipoprotein AII, complement C3, and transthyretin) by approximately 94%. In parallel, non-depleted serum samples were diluted 10× without any pretreatment. Depleted and non-depleted samples were quantified using the Qubit Protein Assay Kit platform (Invitrogen) according to the manufacturer’s instructions. A total of 20 μg proteins were reduced with 10 mM Dithiothreitol (DTT) at 30°C for 45 min and alkylated with 40 mM of iodoacetamide (IAA) for 30 min at room temperature in the dark. The samples were digested with 10% (m/m) trypsin (Promega) during 16 h at 30°C. Following digestion, all reactions were acidified with 1% (v/v) trifluoroacetic acid and tryptic peptides were desalted using C18 in-house stage-tips (3M Empore), dried and suspended in 0.1% formic acid (FA) prior to LC-MS/MS analysis.
Mass Spectrometry Analysis
For analysis of the depleted samples, an UltiMate 3000 Nanoflow LC system (Thermo Scientific) coupled online to a hybrid Quadrupole-Orbitrap mass spectrometer HF-X (Thermo Fisher Scientific) was used. The peptide mixture was loaded on an in-house packed reversed-phase pre-column (4 cm × 100 μm inner diameter, ReproSil-Pur C18-AQ 5 μm particles) and subsequently eluted onto a 20-cm 75-μm inner diameter analytical column containing ReproSil-Pur C18-AQ 3-μm particles. We applied a 66-min gradient using 0 to 35% solvent B in 40 min, 35 to 45% in 15 min, 45 to 99% B in 3 min and 5% B for 2 min (A = 0.1% FA; B = 90% ACN, 0.1% FA). After each run, the trap column and column were equilibrated with mobile phase A. The Quadrupole-Orbitrap HF-X instrument was set to data-dependent acquisition (DDA) and was operated in a positive mode. Survey scans (350–1,600 m/z) were acquired in the Orbitrap system with a resolution of 120,000 at m/z 200. The 20 most intense ions were sequentially isolated and HCD (Higher Energy Collision Dissociation) fragmented using normalized collision energy. The general mass spectrometric conditions were as follows: 2 kV spray voltage, no sheath and auxiliary gas flow, heated capillary temperature of 275°C, predictive automatic gain control (AGC) enabled, and an S-lens RF level of 40%.
The non-depleted samples were analyzed on an LTQ-Orbitrap Velos ETD (Thermo Fisher Scientific) coupled with Easy NanoLC II (Thermo Scientific). The peptide mixture was loaded on a ReproSil-Pur C18-AQ C18 reversed phase column (4 cm × 100 μm inner diameter, 5 μm particles) and subsequently eluted onto a 20 cm 75 inner diameter analytical column containing ReproSil-Pur C18-AQ 3 μm particles. We applied a 130-min gradient using the from 3 to 28% solvent B in 100 min, 28 to 45% in 20 min, 45 - 95% B in 2 min and 8 min at 95% B (A = 0.1% FA; B = 90% ACN, 0.1% FA). After each run, the trap column and column were equilibrated with mobile phase A. The LTQ-Orbitrap Velos instrument was set to data-dependent acquisition (DDA) and was operated in a positive mode. Survey scans (350–1,500 m/z) were acquired in the Orbitrap system with a resolution of 60,000 at m/z 110. The 20 most intense ions, excluding unassigned and 1+ charge state, were sequentially isolated and CID (Collision-induced dissociation) fragmented using normalized collision energy of 35. The general mass spectrometric conditions were as follows: 1.9 kV spray voltage, no sheath and auxiliary gas flow, heated capillary temperature of 280°C, predictive automatic gain control (AGC) enabled, and an S-lens RF level of 65.88%.
Database Search and Statistical Analysis
The raw files corresponding to the depleted and non-depleted serum were searched using Proteome Discoverer v2.3.0.523 and PatternLab for proteomics v4.1.1.17 (Carvalho et al., 2016) (http://patternlabforproteomics.org/), using the SEQUEST search engine. The H. sapiens Swiss-Prot database was downloaded on January 24, 2020; a reversed version of each sequence plus those from 127 common mass spectrometry contaminants was included. Trypsin was used as a cleavage enzyme (fully tryptic and semi-tryptic), allowing a maximum of 2 missed cleavage sites. Cysteine carbamidomethylation and methionine oxidation were considered as a static and dynamic modifications, respectively. In the Proteome Discoverer and PatternLab for Proteomics tools, false Discovery Rate (FDR) was 1% for peptide and protein identification. Label free quantification (LFQ) was performed by applying the extracted ion chromatogram (XIC) area. In PatternLab, the quantitation was performed according to Normalized Ion Abundance Factors (NIAF) as a relative quantitation strategy. NIAF is the equivalent to NSAF (Zybailov et al., 2006), but applied to XIC (Neilson et al., 2011). Statistical analyzes were performed using the Perseus 1.5.3.2 software (Cox and Mann, 2008), Proteome Discoverer, and PatternLab for Proteomics (Carvalho et al., 2016). Differentially regulated proteins and semi-tryptic peptides were determined by applying a t-test with an adjusted p-value <0.1 (Benjamini-Hochberg method).
Bioinformatics Analysis
For gene ontology (GO) analysis, the tool g: Profiler (Raudvere et al., 2019) was employed. The analyses were performed separately for proteins identified upregulated and downregulated between the groups evaluated. A q-value threshold of 0.05 was used, corrected by the Benjamini-Hochberg method. The Gene Enrichment Analysis (GSEA) (Subramanian et al., 2007) was applied to evaluate enriched pathways using the Reactome platform as reference data. Other parameters were used with configurations: permutation method: genes, minimum number of members: 3, maximum number of members: 84, metric for ranking genes: T-test. The Database for Annotation, Visualization and Integrated Discovery (DAVID) v6.8 (Huang et al., 2009) was used to complement bioinformatics analysis. IceLogo tool (Colaert et al., 2009), BRENDA enzyme database (Placzek et al., 2017), Proteasix (http://www.proteasix.org/) and MEROPS - the Peptidase Database (Rawlings et al., 2018) were used to access and analyze peptidases activity.
Analysis of Enzyme Activity
To analyze the enzymatic activity of the samples, a gelatin zymography was performed as previously described (Toth and Fridman, 2001). Proteins were resolved electrophoretically in 12% SDS-PAGE containing 0.1% gelatin. To remove the SDS, the gel was incubated two times of 30 min with wash buffer (2.5% Triton X-100, 50 mM Tris HCl, 5 mM CaCl2, 1 μM ZnCl2, and H2O); followed by a wash in water and incubation for 12 h in incubation buffer (Triton x 100 at 1%, Tris HCl 50 mM pH 7.5, 5 mM CaCl2, 1 µM ZnCl2) at 37°C in a water bath. The gels were stained with 0.5% Coomassie blue and MMPs activity was determined by the intensity of the band using the ImageLab 3.0 software. Statistical significance was assessed by Student t-test using Graphpad Prism 5 software.
Results
The Serum Proteomic Profile of ZIKV Exposed Children Without CZS Is Altered Compared to Controls
The serum of infants without symptoms of CZS, but with intrauterine exposure to ZIKV and unexposed controls was evaluated by a proteomic approach based on mass spectrometry, with and without depletion of the 14 most abundant serum proteins (Figure 1A). The gestational and maternal age was similar between ZIKV and CTRL infants (Supplementary Figures 1A–C and Supplementary File 1). The age of infants varied from 3 to 23 months (CTRL = 14.4 and ZIKV = 13.2) without statistical significance (p-value = 0.3379) (Supplementary Figure 1B and Supplementary File 1). The average head circumference was smaller in the ZIKV group, however, this difference was not statistically significant (p-value = 0.2792) (Supplementary Figure 1D and Supplementary File 1).
Figure 1 Proteins identified in depleted and non-depleted serum. (A) Workflow adopted for sample preparation and data analysis. (B) Venn diagram of proteins identified in the depleted and non-depleted serum. Only proteins identified by the Proteome Discoverer and PatternLab for Proteomics software were considered for posterior analysis. (C, D) Cellular components attributed to the total proteins identified in depleted and non-depleted treatment were determined using the tool g: Profiler. (q-value < 0.05). (E, F) Tissues that express identified proteins were evaluated by the Database for Annotation, Visualization and Integrated Discovery (DAVID) v6.8 platform. Only identifications with a p-value equal to or less than 0.05 and that present at least five proteins were considered.
All mothers presented rash during the second or third trimester of pregnancy. Boys and girls distribution between CTRL and ZIKV groups are six boys and four girls each, and 75% of 20 families were favela residents. One child of the ZIKV exposed group presented developmental delays and was diagnosed with apraxia of speech and attention deficit hyperactivity disorder. This child had abnormal magnetic resonance findings, with the high signal at the periventricular area in T2 and FLAIR, prominence of perivascular spaces with unspecific aspect, and corpus callosum thinning (Supplementary File 1).
Proteomics identified a total of 383 and 449 proteins in the depleted serum using the Proteome Discoverer and PatternLab software, respectively (Supplementary 2). The 383 common proteins between the two software were selected for further statistical analysis. The identifications in the non-depleted serum were 278 and 438 proteins by Proteome Discoverer and PatternLab for Proteomics, respectively. The 278 proteins identified by both search engines were considered for further analysis. A total of 190 proteins were identified in common between the depleted and non-depleted serum: and 193 and 88 exclusive proteins in the non-depleted and depleted treatment, respectively (Figure 1B). The Principal Component Analysis (PCA) of depleted (Supplementary Figure 1E) and non-depleted (Supplementary Figure 1F) serum are different, with separation between two distinct groups only in the depleted serum (Supplementary File 1). GO analysis of each dataset was performed to evaluate the effect of serum depletion on the cellular components (CC) of identified proteins. As expected, the variety of CC identified in the depleted (Figure 1C) and non-depleted (Figure 1D) serum revealed minor variation; however, proteins of the immunoglobulin complex were enriched only in the non-depleted serum. We also evaluated tissues that express the identified proteins (Figures 1E, F). The depleted serum shows greater diversity in the tissue distribution of identified proteins (Figure 1E), with several proteins being expressed in the placenta (Supplementary File 2). In the depleted serum, 84 proteins were regulated, with 36 upregulated and 48 downregulated in the serum of infants with ZIKV intrauterine exposure compared to controls (Figure 2A). A total of 37 proteins showed fold change greater than 1 time (Figure 2B and Table 1) while the non-depleted serum indicated only two regulated proteins, one upregulated and one downregulated (Figure 2C, Supplementary File 3).
Figure 2 Differently abundant proteins identified in the groups of patients CTRL and ZIKV. (A) Heatmap of the proteins regulated in the depleted serum between the evaluated groups. The intensities were normalized using the z-score function. The red and blue colors show upregulated and downregulated proteins, respectively. The blank color indicates that the protein did not show any intensity in the sample. (B, C) Volcano plot of depleted and non-depleted serum proteins, respectively. The fold change is represented by the ratio between ZIKV and CTRL protein intensities; negative values indicate greater abundance in the CTRL group compared to the ZIKV group. In (B) proteins with fold change > 1 were marked, in (C) both regulated proteins were marked.
Table 1 Differently abundant proteins identified in CTRL and ZIKV patient groups.
Disease-Related and Enriched Pathways Analysis
The results of disease-related and enriched protein analysis are broad, due to proteins that are shared between multiple pathological processes (Supplementary File 4). Indeed, evaluation of diseases related to all differentially regulated proteins showed different classes of disorders, which include metabolic, renal, brain, cardiovascular abnormalities, among others (Figure 3A). When evaluation was performed for upregulated proteins, data indicate that proteins associated to macular degeneration, thrombosis, retinopathy of prematurity, and cerebrovascular disease (Figure 3B). These data provide insight of important molecular changes that might be related to exposure to ZIKV in utero.
Figure 3 Disease-related proteins. (A). The “Diseases” module in the Database for Annotation, Visualization and Integrated Discovery (DAVID) tool was used to determine the disease classes related to all regulated proteins. (B) The analysis for proteins related to diseases was performed separately for upregulated and downregulated proteins. Only tissues and diseases that presented p-value <0.05 were represented in the graph.
The GSEA Reactome analysis showed 27 differently regulated pathways (Figure 4A). The colors red and blue represent upregulated and downregulated pathways, respectively; heatmaps indicate proteins and their related pathways. Axon guidance and RHO GTPase signaling were identified upregulated in the ZIKV group, while visual phototransduction and retinoid metabolism are downregulated (Figure 4A). The pathways were represented in relation to p-value; in positive regulation (Figure 4B) and negative regulation (Figure 4C). Complete information on the GSEA analysis is available in Supplementary file 5.
Figure 4 The Gene Enrichment Analysis (GSEA). (A) Results obtained with the GSEA analysis; the reference data used were the Reactome platform pathways. The red and blue dots indicate upregulated and downregulated pathways, respectively. The numbered pathways correspond to: 1: Formation of fibrin clot, 2: Signaling by RHO GTPase, 3: Metabolism of fat-soluble vitamins, 4: Retinoid metabolism and transport, 5: Visual phototransduction, 6: Response to elevated platelet cytosolic Ca2+, 7: Metabolism of lipids, 8: Scavenging of heme from plasma, 9: Plasma lipoprotein remodeling, 10: HDL remodeling, 11: Post-translational protein phosphorylation,12: Regulation of Insulin-like Growth Factor (IGF) transport and uptake by Insulin-like Growth Factor Binding Proteins (IGFBPs), 13: Platelet degranulation, 14: Degradation of the extracellular matrix, 15: Integrin cell surface interactions, 16: Non-integrin membrane-ECM interactions, 17: ECM proteoglycans, 18: Extracellular matrix organization, 19: Signal Transduction, 20: Developmental Biology, 21: Signaling by Interleukins, 22: Immune System, 23: Cytokine signaling in immune system,24: Axon guidance, 25: NCAM1 interactions, 26: Hemostasis, 27: Neutrophil degranulation. Heatmaps represent proteins related to three chosen pathways: Axon guidance, RHO GTPase signaling, and visual phototransduction. (B, C) Upregulated and downregulated pathways, respectively.
GO analysis for upregulated and downregulated proteins was performed for molecular function (MF), biological process (BP), and cellular component (CC). Events related to coagulation, leukocyte aggregation, synapse maturation, fibrinolysis, and platelet degranulation are increased in the ZIKV group (Figure 5A). On the other hand, downregulated proteins are involved with immune responses, tissue homeostasis, and remodeling of the lipid-protein complex. The CC of differently regulated proteins is similar, with emphasis on lipoprotein complexes for downregulated proteins (Figure 5B). Both upregulated and downregulated proteins showed MFs related to the activities of endopeptidases, indicating higher protease activity in the ZIKV exposed serum. The interaction between proteins upregulated (red dots) and downregulated (blue dots) and the respective ontologies are shown in Figure 5C, including TGF-beta signaling and complement activation (Supplementary File 6).
Figure 5 Gene ontology analysis. (A, B) From top to bottom are indicated the molecular function (MF), Biological Process (BP), and Cellular Component (CC) related to upregulated and downregulated proteins, respectively. q-value < 0.05. (C) Interaction between upregulated (red dots) and downregulated (blue dots) proteins, with their respective related ontologies. The gray color show the ontology related to proteins identified upregulated and downregulated. The analysis was performed in the ClueGo app. Only interactions with a p-value less than or equal to 0.05 were considered, with correction by the Benjamini-Hochberg method. The size of the node shows the significance of the interactions, based on the adjusted p-value.
Aberrant Pattern of Protease Activity in the Serum of Infants With ZIKV Intrauterine Exposure Compared to Control
Semi-tryptic peptides are those cleaved at the C-terminus by trypsin (arginine and lysine residues) and at the other terminal (N-terminus) by another endogenous enzyme. The analysis of semi-tryptic peptides can indicate the action of proteases. In depleted serum, we identified a total of 3,258 peptides and among these, 976 were semi-tryptic (Figure 6A). To determine the semi-tryptic peptides differentially regulated between the groups evaluated, the same approach described above for proteins was applied. Our results show 547 semi-tryptic peptides upregulated and 47 downregulated in the ZIKV group (Figure 6B and Supplementary File 7); suggesting a higher proteolytic activity in the serum of ZIKV exposed infants. In total, 87 proteases and protease inhibitors were identified in our dataset (Table 2).
Figure 6 Semi-tryptic peptides. (A) The number of tryptic and semi-tryptic peptides identified in samples of depleted serum. q-value < 0.1. (B) Volcano plot of differently regulated semi-tryptic peptides. A total of 547 upregulated (red dots) and 47 downregulated (blue dots) were identified. (C, D) The most frequent cleavage sites (%) identified by the upregulated and downregulated semi-tryptic peptides, respectively. The numbers 1-6 indicate the position of the amino acids in the peptide sequence and the endopeptidase cleavage site is between positions 1 and 2. In heatmaps, red is the most represented amino acid, while blue is the least represented amino acid.
Table 2 Proteases and protease inhibitors identified in our dataset.
Motif analysis of most frequent cleavage sites for semi-tryptic peptides upregulated in the ZIKV group showed alanine and leucine at the cleaved sites (Figure 6C). On the other hand, the downregulated semi-tryptic peptides showed tyrosine and phenylalanine as the most frequent cleavage sites (Figure 6D). The prediction of proteases involved in the cleavage of the all semi-tryptic peptides identified showed that metalloproteinase family is responsible for the cleavage of 24% and plasminogen of 4%. The individual analysis of upregulated and downregulated semi-tryptic peptides showed that metalloproteinases are responsible for the cleavage of 22.8% and 12.6%, respectively (Table 2 and Supplementary 7).
Evaluation of serum protease activity (Figure 7A) showed an increase in the metalloproteinase 2 (MMP2) (Figure 7B) and metalloproteinase 9 (MMP9) (Figure 7C) in serum from ZIKV group compared to CTRL, which are shown by the intensity of the molecular weight bands of 72 KDa and 92 KDa, respectively (Supplementary File 8). In addition, the formation of MMPs complexes is indicated by the presence of high molecular weight bands (~ 250 KDa) (Roy et al., 2008) (Figure 7D). The intensity of all quantified bands was summed, and the result is shown in (Figure 7E), with p-value <0.0001 (Supplementary File 8).
Figure 7 Gelatin zymography. (A) Gelatinolytic activity of the samples from the CTRL and ZIKV group, with emphasis on the activities of MMP 9 (92 kDa) and MMP 2 (72 KDa). (B–D) The intensities of the molecular weight bands of 72, 92, and 250 kDa, with their respective p-value. (E) Shows the sum of all quantified intensities, with their respective p-value.
Discussion
Exposure to ZIKV Is Related to Coagulation, Thrombotic Events, and Late Cerebrovascular Disorders
Quantitative proteomics approaches have been applied to investigate the molecular mechanisms involved in intrauterine ZIKV infection. Infection of primary human neural stem cells (McGrath et al., 2017), neural progenitor cells (Scaturro et al., 2018), and neurospheres derived from induced pluripotent stem cells are used in modeling the effects of first-trimester infections in vitro (Garcez et al., 2017; Rosa-Fernandes et al., 2019). In this study, we performed large-scale serum quantitative proteomics to identify molecular changes in infants exposed to ZIKV without early clinical symptoms compared to paired controls, in order to explore and predict molecular evidence of possible late abnormalities during the ongoing follow-up of this cohort.
Proteins related to venous thrombosis and blood coagulation, such as Coagulation factor VII (F7) and Fibrinogen alpha chain (FGA), were upregulated in the ZIKV group. The relationship between viral infections and changes in coagulation processes has already been described for ZIKV (Ramacciotti et al., 2019) and other viruses (Kimmel, 1967; Bibas et al., 2011; da Costa et al., 2012; Goeijenbier et al., 2012; Roy et al., 2013; Wang et al., 2015; Wijarnpreecha et al., 2017; Marques et al., 2017; Neppelenbroek et al., 2018; Ngu et al., 2018; Ramacciotti et al., 2019).
After the occurrence of isolated cases of venous thrombosis in patients who had a positive diagnosis for ZIKV, Ramacciotti et al. (2019) assessed blood D-dimer levels, which are usually monitored for the diagnosis of deep venous thrombosis, of 172 patients who had ZIKV or chikungunya, without cross-infection. The results showed an increase of 19.4% and 63.8% in the D-dimer levels of patients with ZIKV and chikungunya, respectively. Moreover, pro-thrombosis effects related to viral infection of the CNS can result in cerebrovascular complications (Tang et al., 2019). Landais et al. (2017) reported a case of stroke in a 10-month-old child who was positive for ZIKV. A second case of fetal cerebral infarction has been reported, in which a child has been exposed to ZIKV and had no brain changes on MRI. However, at 16 days of age, an area of chronic encephalomalacia was found (Mulkey et al., 2018). A third study confirmed the presence of ZIKV RNA in the placenta, umbilical cord and amniotic membrane on the maternal side of the placenta of a pregnant woman who delivered a child without clinical evidence of microcephaly and CZS. The medical follow-up showed that the child had a stroke on the eighth day of life, even without apparent abnormalities in complete blood count, cranial ultrasound and ammonia values (Raymond and Jakus, 2018). In addition, a report of cerebral vasculitis in an adult patient with ZIKV positive PCR was published (Acevedo et al., 2017). Cerebral vasculitis in neonates is characteristic of congenital diseases (Koeppen et al., 1981; Baskin and Hedlund, 2007; Hauer et al., 2019). Although case reports related to ZIKV exposure are punctual, our study agrees with these findings as we identified proteins related to upregulated cerebrovascular disease in children exposed to ZIKV compared to CTRL group. In our cohort, a patient in the ZIKV group showed changes in MRI and developmental delay. This finding corroborates previous studies applied to the larger cohorts, such as Pomar et al. (2018) and Rice et al. (2018), in which 12.5% and 1.6% of the population exposed to ZIKV who were born without microcephaly and CZS developed complications, respectively. Due do that, we emphasize the importance of medical monitoring of symptomatic and asymptomatic children, as apparently unexplained late neurological complications might be due to exposure to ZIKV in the uterus.
Cerebrovascular Disorders Related to Increased Activity of Metalloproteinases
MMPs are proteolytic enzymes that degrade the extracellular matrix and basement membranes and participate in biological homeostatic and pathological processes (Birkedal-Hansen et al., 1993). Among the neuronal pathologies linked to MMPs, neuronal apoptosis and oxidative damage to DNA are apparent (Xie et al., 2017). MMPs are expressed in different tissues, including the CNS, where they perform pathological functions linked to the opening of the blood-brain barrier after cerebral ischemia (Cunningham et al., 2005). Reports show that MMPs play a key role in chronic inflammatory diseases of the CNS and participate in the degradation of myelin components (Walker and Rosenberg, 2010). Our study identified an increased activity of MMPs 2 and 9 in the serum of children exposed to ZIKV, as well as an increase in the levels of the Profilin-1 (PFN-1) protein, which is fundamental for myelination (Montani et al., 2014). Neuronal death after cerebrovascular disorders has been reported, and studies show a possible role for MMPs in this phenomenon. Lee (2004) showed increased activity of MMP2 and 9 located mainly in neurons of rats after transient global ischemia. Yang et al. (2010) demonstrated that MMP2 and 9 cleave the proteins Poly(ADP-Ribose) Polymerase 1 (PARP-1) and DNA repair protein XRCC1 (XRCC1), which facilitates oxidative damage in neurons during early ischemia. These proteins are involved in events of cell survival and apoptosis.
Macular Degeneration and Eye Abnormalities Related Proteins
Ocular complications have already been described in infants diagnosed with CZS (Ventura and Ventura, 2018). Moreover, a recently developed label-free proteomics methodology has also been applied to evaluate the alteration of the ocular protein of children exposed to ZIKV during pregnancy, comparing groups with and without CZS (Rosa-Fernandes et al., 2020). The study showed that the biomolecules involved in neutrophil degranulation, neurodevelopment, cell infiltration, and ocular dysfunction were identified in different abundances between the groups evaluated, being able to discriminate infants exposed to ZIKV during gestation and without early clinical symptoms (Rosa-Fernandes et al., 2020).
In addition, other flaviviruses, such as West Nile fever, transmitted from the mother to the fetus have also been shown to cause eye abnormalities (Alpert et al., 2003). In our analysis, we identified proteins related to premature retinopathy and macular degeneration, such as pigment epithelium-derived factor (SERPINF1), transforming growth factor beta-1 proprotein (TGFB1), fibrinogen alpha chain (FGA), mannan-binding lectin serine protease 1 (MASP1), moesin (MSN), and vitronectin (VTN), upregulated in the ZIKV group, suggesting the possibility of developing these complications. Our analysis of enriched pathways showed that visual phototransduction events are downregulated in this group compared to CTRL. Visual phototransduction is a photochemical and biochemical process that consists of photon absorption by photoreceptor cells, which convert this signal into an electrical cellular response. This electrical response is sent to the brain through action potentials and electrochemical changes. Therefore, the normal activity of this pathway, which is highly conserved in many species, is essential for the proper functioning of vision neurobiology (Mannu, 2014).
Our data together with other reports in the literature, reinforce the possibility of the occurrence of abnormalities related to ZIKV infection, not only including neurological, but also ophthalmological disorders.
The children included in this study are currently monitored by a multidisciplinary clinical team to assess the outcome of exposure to ZIKV. Here, we demonstrate molecular differences between the groups evaluated and we emphasize the activity of MMP2 and MMP9 can emerge as potential biomarkers of exposure to the virus. We have demonstrated the increased activity of MMPs in the cohort exposed to ZIKV, making it possible to apply this finding to monitor MMPs activity in children exposed to ZIKV and who may develop late abnormalities.
Conclusions
Serum provides information on the entire content of circulating proteins and, as the main fluid available for routine clinical evaluations, has great potential in diagnostic and prognostic analyses. We have shown that children exposed to ZIKV during pregnancy, but who were born without any complications, can present molecular evidence for late abnormalities related to CZS. Our study has shown alterations in proteins that participate in processes related to neuronal death and cerebrovascular abnormalities in the ZIKV group, even if these children do not present clinical evidences of CZS at birth. In addition, vision-related proteins have been identified as downregulated, which may indicate ocular and visual impairments, a frequent characteristic in infants that develop CZS. Another key finding was the increased activity of MMP-2 and MMP-9 in all serum samples in the ZIKV group, which could be associated to neuronal death. The children included in this study are currently monitored by multidisciplinary clinical teams to assess the outcome of long-term exposure to the virus. Our study is the first to assess molecular alterations for late disorders in child victims of the ZIKV epidemic in the Americas, demonstrating that medical follow-up should be carried out on all children exposed to the virus, as late complications can occur.
Data Availability Statement
The original contributions presented in the study are publicly available. The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository with the dataset identifier PXD020294.
Ethics Statement
The studies involving human participants were reviewed and approved by Institutional review board and the ethics committee of the Universidade Federal Fluminense (protocol CAAE number 79890517.6.0000.5243). Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin.
Author Contributions
Conceptualization, methodology and project administration : LR-F and GP. Formal analysis: JM-d-S LR-F and GP. Investigation: JM-d-S, LR-F, RB, CA. Resources: RB, FC, RV, PC, ML and CC. Validation: JM-d-S. Visualization: JM-d-S and LR-F. Writing – original draft preparation: JM-d-S, LR-F and GP. Writing – review and editing: all authors. All authors contributed to the article and approved the submitted version.
Funding
RB was supported by FAPERJ (201.779/2017); GP was supported by FAPESP (2014/06863-3, 2018/18257-1, 2018/15549-1) and CNPq (“Bolsa de Produtividade”). This work was also supported by the VILLUM Center for Bioanalytical Sciences at the University of Southern Denmark.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
The authors would like to thank the medical team at the Hospital Universitário Antonio Pedro at Universidade Federal Fluminense for assisting with the sample’s collections, and all mothers and children included in this project. We would like to thank Professor Luigi Palmieri, University of Bari, Italy, for useful commentaries and for reviewing this paper.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2020.577819/full#supplementary-material
Supplementary Figure 1 | Details of the participants. (A) Gestational age (GA) of birth of the study participants. The ZIKV group has nine points, as one participant GA was not noted. (B) The age of blood collection for the study. A t-test was applied, and no statistical significance was observed between the groups (p-value = 0.3379). (C) Age of the mothers included in the study. (D) Circumference of the infants heads in centimeters. The measurements were taken in the first week of the participants life. The ZIKV group has nine points, as one participant head circumference was not noted. (E, F) Principal Component Analysis (PCA) for depleted and non-depleted serum, respectively. The yellow and blue dots represent the CTRL and ZIKV groups, respectively. The “Principal Component Analysis” module, available in the Perseus tool, was used to build the graph. Previously, a filter was applied to select only proteins present in all samples.
Supplementary File 1 | Demographic, clinical, and radiological characteristics of the study population. Data relating to children included in this study, indicating results of maternal PCR, maternal age, child’s age when collecting serum samples, gestational age of maternal rash, head circumference in the first week of life, gestational age of birth, TORCH analysis, gender, place of residence and clinical follow-up data.
Supplementary File 2 | Proteins identified by Proteome Discoverer v2.3.0.523 and PatternLab for proteomics v4.1.1.17 in depleted and non-depleted serum; cellular component analysis, and tissue-related protein. Total signal (NSAF): sum of PatternLab’s normalized label-free quantitation derived from the extracted ion chromatograms (XIC). Adjusted p-value refers to the application of the correction by the Benjamini-Hochberg method. Description: protein description from the Swiss-Prot database. Name: identifier from SwissProt/UniProt database.
Supplementary File 3 | Complete information for protein differently regulated between the CTRL and ZIKV groups. Protein ID: Swiss-Prot protein identifier. Fold change: ratio between ZIKV and CTRL protein intensities; negative values indicate greater abundance in the CTRL group compared to the ZIKV group. Protein description: according to the Swiss-Prot database. All proteins satisfy a q-value < 0.1, corrected by Benjamini-Hochberg method. In the analysis of the heat map, the intensities were normalized using the z-score function. Positive values indicate that the data is above average and, when negative, it means that the data is below average.
Supplementary File 4 | Disease-related proteins. Complete information from DAVID analysis for all differently regulated, upregulated, and downregulated proteins. GAD_Disease: Gene-Disease Associations.
Supplementary File 5 | Full list of pathways disclosed by the GSEA analysis. The results include the pathways identifier, statistical results, and name from Reactome.
Supplementary File 6 | Complete information of the DAVID and ClueGO Gene Ontology analysis. CTRL and ZIKV GO information. BP: biological process; CC: cellular component; MF: molecular function.
Supplementary File 7 | Analysis of semi-tryptic peptides. Result of the search performed by Proteome Discoverer software. Trypsin was used as a cleavage enzyme (semi-tryptic), allowing a maximum of 2 missed cleavage sites. Description: protein description from the Swiss-Prot database. Name: identifier from SwissProt/UniProt database.
Supplementary File 8 | Gelatin zymography. The intensity of the bands corresponds to the results of the ImageLab 3.0 software.
References
Acevedo N., Waggoner J., Rodriguez M., Rivera L., Landivar J., Pinsky B., et al. (2017). Zika Virus, Chikungunya Virus, and Dengue Virus in Cerebrospinal Fluid from Adults with Neurological Manifestations, Guayaquil, Ecuador. Front. Microbiol. 8, 42. doi: 10.3389/fmicb.2017.00042
Alpert S. G., Fergerson J., Noël L.-P. (2003). Intrauterine West Nile virus: ocular and systemic findings. Am. J. Ophthalmol. 136, 733–735. doi: 10.1016/S0002-9394(03)00452-5
Baskin H. J., Hedlund G. (2007). Neuroimaging of herpesvirus infections in children. Pediatr. Radiol. 37, 949–963. doi: 10.1007/s00247-007-0506-1
Bibas M., Biava G., Antinori A. (2011). HIV-ASSOCIATED VENOUS THROMBOEMBOLISM. Mediterr. J. Hematol. Infect. Dis. 3, e2011030. doi: 10.4084/mjhid.2011.030
Birkedal-Hansen H., Moore W. G., Bodden M. K., Windsor L. J., Birkedal-Hansen B., DeCarlo A., et al. (1993). Matrix metalloproteinases: a review. Crit. Rev. Oral. Biol. Med. 4, 197–250. doi: 10.1177/10454411930040020401
Carvalho P. C., Lima D. B., Leprevost F. V., Santos M. D. M., Fischer J. S. G., Aquino P. F., et al. (2016). Integrated analysis of shotgun proteomic data with PatternLab for proteomics 4.0. Nat. Protoc. 11, 102–117. doi: 10.1038/nprot.2015.133
Colaert N., Helsens K., Martens L., Vandekerckhove J., Gevaert K. (2009). Improved visualization of protein consensus sequences by iceLogo. Nat. Methods 6, 786–787. doi: 10.1038/nmeth1109-786
Cox J., Mann M. (2008). MaxQuant enables high peptide identification rates, individualized p.p.b.-range mass accuracies and proteome-wide protein quantification. Nat. Biotechnol. 26, 1367–1372. doi: 10.1038/nbt.1511
Cunningham L. A., Wetzel M., Rosenberg G. A. (2005). Multiple roles for MMPs and TIMPs in cerebral ischemia. Glia 50, 329–339. doi: 10.1002/glia.20169
da Costa P. S. G., Ribeiro G. M., Junior C. S., da Costa Campos L. (2012). Severe thrombotic events associated with dengue fever, Brazil. Am. J. Trop. Med. Hyg. 87, 741–742. doi: 10.4269/ajtmh.2012.11-0692
Dahle A. J., Fowler K. B., Wright J. D., Boppana S. B., Britt W. J., Pass R. F. (2000). Longitudinal investigation of hearing disorders in children with congenital cytomegalovirus. J. Am. Acad. Audiol. 11 (5), 283–290.
Dick G. W. A., Kitchen S. F., Haddow A. J. (1952). Zika virusand serological specificity. Trans. R. Soc Trop. Med. Hyg. 46, 509–520. doi: 10.1016/0035-9203(52)90042-4
Garcez P. P., Nascimento J. M., de Vasconcelos J. M., Madeiro da Costa R., Delvecchio R., Trindade P., et al. (2017). Zika virus disrupts molecular fingerprinting of human neurospheres. Sci. Rep. 7:40780. doi: 10.1038/srep40780
Goeijenbier M., van Wissen M., van de Weg C., Jong E., Gerdes V. E. A., Meijers J. C. M., et al. (2012). Review: Viral infections and mechanisms of thrombosis and bleeding. J. Med. Virol. 84, 1680–1696. doi: 10.1002/jmv.23354
Hauer L., Pikija S., Schulte E. C., Sztriha L. K., Nardone R., Sellner J. (2019). Cerebrovascular manifestations of herpes simplex virus infection of the central nervous system: a systematic review. J. Neuroinflamm. 16, 19. doi: 10.1186/s12974-019-1409-4
Heukelbach J., Alencar C. H., Kelvin A. A., de Oliveira W. K., Pamplona de Góes Cavalcanti L. (2016). Zika virus outbreak in Brazil. J. Infect. Dev. Ctries 10, 116–120. doi: 10.3855/jidc.8217
Huang D. W., Sherman B. T., Lempicki R. A. (2009). Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc. 4, 44–57. doi: 10.1038/nprot.2008.211
Kimmel C. B. (1967). Lysosomes in the spleen of the chick embryo. II. Changes during the graft-versus-host reaction. J. Exp. Zool. 166, 447–458. doi: 10.1002/jez.1401660316
Kindhauser M. K., Allen T., Frank V., Santhana R. S., Dye C. (2016). Zika: the origin and spread of a mosquito-borne virus. Bull. World Health Organ. 94, 675–686C. doi: 10.2471/BLT.16.171082
Koeppen A. H., Lansing L. S., Peng S. K., Smith R. S. (1981). Central nervous system vasculitis in cytomegalovirus infection. J. Neurol. Sci. 51, 395–410. doi: 10.1016/0022-510x(81)90117-9
Landais A., Césaire A., Fernandez M., Breurec S., Herrmann C., Delion F., et al. (2017). ZIKA vasculitis: A new cause of stroke in children? J. Neurol. Sci. 383, 211–213. doi: 10.1016/j.jns.2017.10.045
Lee S.-R. (2004). Role of Matrix Metalloproteinases in Delayed Neuronal Damage after Transient Global Cerebral Ischemia. J. Neurosci. 24, 671–678. doi: 10.1523/JNEUROSCI.4243-03.2004
Lowe R., Barcellos C., Brasil P., Cruz O. G., Honório N. A., Kuper H., et al. (2018). The Zika Virus Epidemic in Brazil: From Discovery to Future Implications. Int. J. Environ. Res. Public Health 15, 96. doi: 10.3390/ijerph15010096
Marques M. A., Adami de Sá F. P., Lupi O., Brasil P., von Ristow A. (2017). Trombose venosa profunda e vírus chicungunha. J. Vasc. Bras. 16, 60–62. doi: 10.1590/1677-5449.009616
McGrath E. L., Rossi S. L., Gao J., Widen S. G., Grant A. C., Dunn T. J., et al. (2017). Differential Responses of Human Fetal Brain Neural Stem Cells to Zika Virus Infection. Stem Cell Rep. 8, 715–727. doi: 10.1016/j.stemcr.2017.01.008
Mlakar J., Korva M., Tul N., Popović M., Poljšak-Prijatelj M., Mraz J., et al. (2016). Zika Virus Associated with Microcephaly. N. Engl. J. Med. 374, 951–958. doi: 10.1056/NEJMoa1600651
Mohr E. L., Block L. N., Newman C. M., Stewart L. M., Koenig M., Semler M., et al. (2018). Ocular and uteroplacental pathology in a macaque pregnancy with congenital Zika virus infection. PloS One 13, e0190617. doi: 10.1371/journal.pone.0190617
Montani L., Buerki-Thurnherr T., de Faria J. P., Pereira J. A., Dias N. G., Fernandes R., et al. (2014). Profilin 1 is required for peripheral nervous system myelination. Development 141, 1553–1561. doi: 10.1242/dev.101840
Moore C. A., Staples J. E., Dobyns W. B., Pessoa A., Ventura C. V., Fonseca E. B., et al. (2017). Characterizing the Pattern of Anomalies in Congenital Zika Syndrome for Pediatric Clinicians. JAMA Pediatr. 171, 288–295. doi: 10.1001/jamapediatrics.2016.3982
Mulkey S. B., DeBiasi R. L., du Plessis A. J. (2018). Cerebral infarction due to Zika virus. J. Neurol. Sci. 387, 109–110. doi: 10.1016/j.jns.2018.01.032
Musso D., Ko A. II, Baud D. (2019). Zika Virus Infection — After the Pandemic. N. Engl. J. Med. 381, 1444–1457. doi: 10.1056/NEJMra1808246
Neilson K. A., Ali N. A., Muralidharan S., Mirzaei M., Mariani M., Assadourian G., et al. (2011). Less label, more free: Approaches in label-free quantitative mass spectrometry. PROTEOMICS 11, 535–553. doi: 10.1002/pmic.201000553
Neppelenbroek S. II, Rootjes P. A., Boxhoorn L., Wagenaar J. F. P., Simsek S., Stam F. (2018). Cytomegalovirus-associated thrombosis. Neth J. Med. 76 (5), 251–254.
Netto E. M., Moreira-Soto A., Pedroso C., Höser C., Funk S., Kucharski A. J., et al. (2017). High Zika Virus Seroprevalence in Salvador, Northeastern Brazil Limits the Potential for Further Outbreaks. mBio 8, e01390-17. doi: 10.1128/mBio.01390-17
Ngu S., Narula N., Jilani T. N., Bershadskiy A. (2018). Venous Thrombosis Secondary to Acute Cytomegalovirus Infection in an Immunocompetent Host: Consideration for New Screening Guidelines. Cureus 10 (6), e2742. doi: 10.7759/cureus.2742
Noorbakhsh F., Abdolmohammadi K., Fatahi Y., Dalili H., Rasoolinejad M., Rezaei F., et al. (2019). Zika Virus Infection, Basic and Clinical Aspects: A Review Article. Iran. J. Public Health 48, 20–31. doi: 10.18502/ijph.v48i1.779
Pan American Health Organization / World Health Organization. (2017). Zika - Epidemiological Report Brazil. (Washington, D.C.: PAHO/WHO). Available at: https://www.paho.org/hq/dmdocuments/2017/2017-phe-zika-situation-report-bra.pdf. Accessed on 08 November 2020.
Placzek S., Schomburg I., Chang A., Jeske L., Ulbrich M., Tillack J., et al. (2017). BRENDA in 2017: new perspectives and new tools in BRENDA. Nucleic Acids Res. 45, D380–D388. doi: 10.1093/nar/gkw952
Pomar L., Vouga M., Lambert V., Pomar C., Hcini N., Jolivet A., et al. (2018). Maternal-fetal transmission and adverse perinatal outcomes in pregnant women infected with Zika virus: prospective cohort study in French Guiana. BMJ 363, k4431. doi: 10.1136/bmj.k4431
Ramacciotti E., Agati L. B., Aguiar V. C. R., Wolosker N., Guerra J. C., de Almeida R. P., et al. (2019). Zika and Chikungunya Virus and Risk for Venous Thromboembolism. Clin. Appl. Thromb. Hemost. 25, 1076029618821184. doi: 10.1177/1076029618821184
Raudvere U., Kolberg L., Kuzmin I., Arak T., Adler P., Peterson H., et al. (2019). g:Profiler: a web server for functional enrichment analysis and conversions of gene listupdate). Nucleic Acids Res. 47, W191–W198. doi: 10.1093/nar/gkz369
Rawlings N. D., Barrett A. J., Thomas P. D., Huang X., Bateman A., Finn R. D. (2018). The MEROPS database of proteolytic enzymes, their substrates and inhibitors in 2017 and a comparison with peptidases in the PANTHER database. Nucleic Acids Res. 46, D624–D632. doi: 10.1093/nar/gkx1134
Raymond A., Jakus J. (2018). Cerebral Infarction and Refractory Seizures in a Neonate with Suspected Zika Virus Infection. Pediatr. Infect. Dis. J. 37, e112–e114. doi: 10.1097/INF.0000000000001832
Rice M. E., Galang R. R., Roth N. M., Ellington S. R., Moore C. A., Valencia-Prado M., et al. (2018). Vital Signs: Zika-Associated Birth Defects and Neurodevelopmental Abnormalities Possibly Associated with Congenital Zika Virus Infection - U.S. Territories and Freely Associated State. MMWR Morb. Mortal. Wkly. Rep. 67, 858–867. doi: 10.15585/mmwr.mm6731e1
Rosa-Fernandes L., Barbosa R. H., dos Santos M. L. B., Angeli C. B., Silva T. P., Melo R. C. N., et al. (2020). Cellular Imprinting Proteomics Assay: A Novel Method for Detection of Neural and Ocular Disorders Applied to Congenital Zika Virus Syndrome. J. Proteome Res. 19 (11), 4496–4515. doi: 10.1021/acs.jproteome.0c00320
Rosa-Fernandes L., Cugola F. R., Russo F. B., Kawahara R., de Melo Freire C. C., Leite P. E. C., et al. (2019). Zika Virus Impairs Neurogenesis and Synaptogenesis Pathways in Human Neural Stem Cells and Neurons. Front. Cell Neurosci. 13, 64. doi: 10.3389/fncel.2019.00064
Roy R., Louis G., Loughlin K. R., Wiederschain D., Kilroy S. M., Lamb C. C., et al. (2008). Tumor-specific urinary matrix metalloproteinase fingerprinting: identification of high molecular weight urinary matrix metalloproteinase species. Clin. Cancer Res. 14, 6610–6617. doi: 10.1158/1078-0432.CCR-08-1136
Roy A., Chaudhuri J., Chakraborty S., Jasani B., Nanavati R., Jain V., et al. (2013). Case reports. Indian Pediatr. 50, 1053–1057. doi: 10.1007/s13312-013-0269-8
Scaturro P., Stukalov A., Haas D. A., Cortese M., Draganova K., Płaszczyca A., et al. (2018). An orthogonal proteomic survey uncovers novel Zika virus host factors. Nature 561, 253–257. doi: 10.1038/s41586-018-0484-5
Soares F., Abranches A. D., Villela L., Lara S., Araújo D., Nehab S., et al. (2019). Zika virus infection in pregnancy and infant growth, body composition in the first three months of life: a cohort study. Sci. Rep. 9, 19198. doi: 10.1038/s41598-019-55598-6
Souza I. N. O., Barros-Aragão F. G. Q., Frost P. S., Figueiredo C. P., Clarke J. R. (2019). Late Neurological Consequences of Zika Virus Infection: Risk Factors and Pharmaceutical Approaches. Pharmaceut. (Basel) 12, 60. doi: 10.3390/ph12020060
Subramanian A., Kuehn H., Gould J., Tamayo P., Mesirov J. P. (2007). GSEA-P: a desktop application for Gene Set Enrichment Analysis. Bioinformatics 23, 3251–3253. doi: 10.1093/bioinformatics/btm369
Tang C., Patel N. K., Miller R., Beale T., Hyare H. (2019). Spectrum of neurovascular complications from central nervous system infections (viral, bacterial and fungal). BJR|Open 1:bjro.20180024. doi: 10.1259/bjro.20180024
Toth M., Fridman R. (2001). Assessment of Gelatinases (MMP-2 and MMP-9 by Gelatin Zymography. Methods Mol. Med. 57, 163–174. doi: 10.1385/1-59259-136-1:163
Valdes V., Zorrilla C. D., Gabard-Durnam L., Muler-Mendez N., Rahman Z. II, Rivera D., et al. (2019). Cognitive Development of Infants Exposed to the Zika Virus in Puerto Rico. JAMA Netw. Open 2, e1914061. doi: 10.1001/jamanetworkopen.2019.14061
Venancio F. A., Bernal M. E. Q., Ramos M., da C. de B. V., Chaves N. R., Hendges M. V., et al. (2019). Congenital Zika Syndrome in a Brazil-Paraguay-Bolivia border region: Clinical features of cases diagnosed between 2015 and 2018. PloS One 14, e0223408. doi: 10.1371/journal.pone.0223408
Ventura C. V., Ventura L. O. (2018). Ophthalmologic Manifestations Associated With Zika Virus Infection. Pediatrics 141, S161–S166. doi: 10.1542/peds.2017-2038E
Vianna R. A., de O., Lovero K. L., Oliveira S. A., Fernandes A. R., Santos T. C. S. D., et al. (2019). Children Born to Mothers with Rash During Zika Virus Epidemic in Brazil: First 18 Months of Life. J. Trop. Pediatr. 65, 592–602. doi: 10.1093/tropej/fmz019
Walker E. J., Rosenberg G. A. (2010). Divergent role for MMP-2 in myelin breakdown and oligodendrocyte death following transient global ischemia. J. Neurosci. Res. 88, 764–773. doi: 10.1002/jnr.22257
Wang C.-C., Chang C.-T., Lin C.-L., Lin I.-C., Kao C.-H. (2015). Hepatitis C Virus Infection Associated With an Increased Risk of Deep Vein Thrombosis: A Population-Based Cohort Study. Med. (Baltimore) 94, e1585. doi: 10.1097/MD.0000000000001585
Wheeler A. C. (2018). Development of Infants With Congenital Zika Syndrome: What Do We Know and What Can We Expect? Pediatrics 141, S154–S160. doi: 10.1542/peds.2017-2038D
Wijarnpreecha K., Thongprayoon C., Panjawatanan P., Ungprasert P. (2017). Hepatitis C Virus Infection and Risk of Venous Thromboembolism: A Systematic Review and Meta-Analysis. Ann. Hepatol. 16, 514–520. doi: 10.5604/01.3001.0010.0279
Xie Y., Mustafa A., Yerzhan A., Merzhakupova D., Yerlan P., N Orakov A., et al. (2017). Nuclear matrix metalloproteinases: functions resemble the evolution from the intracellular to the extracellular compartment. Cell Death Discovery 3, 17036. doi: 10.1038/cddiscovery.2017.36
Yang Y., Candelario-Jalil E., Thompson J. F., Cuadrado E., Estrada E. Y., Rosell A., et al. (2010). Increased intranuclear matrix metalloproteinase activity in neurons interferes with oxidative DNA repair in focal cerebral ischemia. J. Neurochem. 112, 134–149. doi: 10.1111/j.1471-4159.2009.06433.x
Keywords: Zika virus, late abnormalities, serum proteomics, mass spectrometry, biomarker, Congenital Zika Syndrome
Citation: Macedo-da-Silva J, Rosa-Fernandes L, Barbosa RH, Angeli CB, Carvalho FR, de Oliveira Vianna RA, Carvalho PC, Larsen MR, Cardoso CA and Palmisano G (2020) Serum Proteomics Reveals Alterations in Protease Activity, Axon Guidance, and Visual Phototransduction Pathways in Infants With In Utero Exposure to Zika Virus Without Congenital Zika Syndrome. Front. Cell. Infect. Microbiol. 10:577819. doi: 10.3389/fcimb.2020.577819
Received: 30 June 2020; Accepted: 20 October 2020;
Published: 18 November 2020.
Edited by:
Sigrun Lange, University of Westminster, United KingdomReviewed by:
Emma Mohr, University of Wisconsin-Madison, United StatesMichael F. Criscitiello, Texas A&M University, United States
Copyright © 2020 Macedo-da-Silva, Rosa-Fernandes, Barbosa, Angeli, Carvalho, de Oliveira Vianna, Carvalho, Larsen, Cardoso and Palmisano. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Giuseppe Palmisano, cGFsbWlzYW5vLmdwQGdtYWlsLmNvbQ==; Lívia Rosa-Fernandes, bGl2aWFyb3NhLmZAZ21haWwuY29t
†These authors have contributed equally to this work