Immunoglobulin G N-glycan markers of accelerated biological aging during chronic HIV infection.

People living with HIV (PLWH) experience increased vulnerability to premature aging and inflammation-associated comorbidities, even when HIV replication is suppressed by antiretroviral therapy (ART). However, the factors associated with this vulnerability remain uncertain. In the general population, alterations in the N-glycans on IgGs trigger inflammation and precede the onset of aging-associated diseases. Here, we investigate the IgG N-glycans in cross-sectional and longitudinal samples from 1214 women and men, living with and without HIV. PLWH exhibit an accelerated accumulation of pro-aging-associated glycan alterations and heightened expression of senescence-associated glycan-degrading enzymes compared to controls. These alterations correlate with elevated markers of inflammation and the severity of comorbidities, potentially preceding the development of such comorbidities. Mechanistically, HIV-specific antibodies glycoengineered with these alterations exhibit a reduced ability to elicit anti-HIV Fc-mediated immune activities. These findings hold potential for the development of biomarkers and tools to identify and prevent premature aging and comorbidities in PLWH.

Release of P-TEFb from the Super Elongation Complex promotes HIV-1 latency reversal.

The persistence of HIV-1 in long-lived latent reservoirs during suppressive antiretroviral therapy (ART) remains one of the principal barriers to a functional cure. Blocks to transcriptional elongation play a central role in maintaining the latent state, and several latency reversal strategies focus on the release of positive transcription elongation factor b (P-TEFb) from sequestration by negative regulatory complexes, such as the 7SK complex and BRD4. Another major cellular reservoir of P-TEFb is in Super Elongation Complexes (SECs), which play broad regulatory roles in host gene expression. Still, it is unknown if the release of P-TEFb from SECs is a viable latency reversal strategy. Here, we demonstrate that the SEC is not required for HIV-1 replication in primary CD4+ T cells and that a small molecular inhibitor of the P-TEFb/SEC interaction (termed KL-2) increases viral transcription. KL-2 acts synergistically with other latency reversing agents (LRAs) to reactivate viral transcription in several cell line models of latency in a manner that is, at least in part, dependent on the viral Tat protein. Finally, we demonstrate that KL-2 enhances viral reactivation in peripheral blood mononuclear cells (PBMCs) from people living with HIV on suppressive ART, most notably in combination with inhibitor of apoptosis protein antagonists (IAPi). Taken together, these results suggest that the release of P-TEFb from cellular SECs may be a novel route for HIV-1 latency reactivation.

Impact of HLA class I functional divergence on HIV control.

Heterozygosity of Human leukocyte antigen (HLA) class I genes is linked to beneficial outcomes after HIV infection, presumably through greater breadth of HIV epitope presentation and cytotoxic T cell response. Distinct allotype pairs, however, differ in the extent to which they bind shared sets of peptides. We developed a functional divergence metric that measures pairwise complementarity of allotype-associated peptide binding profiles. Greater functional divergence for pairs of HLA-A and/or HLA-B allotypes was associated with slower AIDS progression and independently with enhanced viral load control. The metric predicts immune breadth at the peptide level rather than gene level and redefines HLA heterozygosity as a continuum differentially affecting disease outcome. Functional divergence may affect response to additional infections, vaccination, immunotherapy, and other diseases where HLA heterozygote advantage occurs.

Association of urine biomarkers of kidney health with subclinical cardiovascular disease among men with and without HIV.

OBJECTIVE: The aim of this study was to determine whether urine biomarkers of kidney health are associated with subclinical cardiovascular disease among men with and without HIV. DESIGN: A cross-sectional study within the Multicenter AIDS Cohort Study (MACS) among 504 men with and without HIV infection who underwent cardiac computed tomography scans and had urine biomarkers measured within the preceding 2 years. METHODS: Our primary predictors were four urine biomarkers of endothelial (albuminuria), proximal tubule dysfunction (alpha-1-microglobulin [A1 M] and injury (kidney injury molecule-1 [KIM-1]) and tubulointerstitial fibrosis (pro-collagen-III N-terminal peptide [PIIINP]). These were evaluated for association with coronary artery calcium (CAC) prevalence, CAC extent, total plaque score, and total segment stenosis using multivariable regression. RESULTS: Of the 504 participants, 384 were men with HIV (MWH) and 120 were men without HIV. In models adjusted for sociodemographic factors, cardiovascular disease risk factors, eGFR, and HIV-related factors, each two-fold higher concentration of albuminuria was associated with a greater extent of CAC (1.35-fold higher, 95% confidence interval 1.11-1.65), and segment stenosis (1.08-fold greater, 95% confidence interval 1.01-1.16). Associations were similar between MWH and men without HIV in stratified analyses. The third quartile of A1 M showed an association with greater CAC extent, total plaque score, and total segment stenosis, compared with the lowest quartile. CONCLUSION: Worse endothelial and proximal tubule dysfunction, as reflected by higher urine albumin and A1 M, were associated with greater CAC extent and coronary artery stenosis.

Proteomic and genetic analyses of influenza A viruses identify pan-viral host targets.

Influenza A Virus (IAV) is a recurring respiratory virus with limited availability of antiviral therapies. Understanding host proteins essential for IAV infection can identify targets for alternative host-directed therapies (HDTs). Using affinity purification-mass spectrometry and global phosphoproteomic and protein abundance analyses using three IAV strains (pH1N1, H3N2, H5N1) in three human cell types (A549, NHBE, THP-1), we map 332 IAV-human protein-protein interactions and identify 13 IAV-modulated kinases. Whole exome sequencing of patients who experienced severe influenza reveals several genes, including scaffold protein AHNAK, with predicted loss-of-function variants that are also identified in our proteomic analyses. Of our identified host factors, 54 significantly alter IAV infection upon siRNA knockdown, and two factors, AHNAK and coatomer subunit COPB1, are also essential for productive infection by SARS-CoV-2. Finally, 16 compounds targeting our identified host factors suppress IAV replication, with two targeting CDK2 and FLT3 showing pan-antiviral activity across influenza and coronavirus families. This study provides a comprehensive network model of IAV infection in human cells, identifying functional host targets for pan-viral HDT.

Plasma Glycomic Markers of Accelerated Biological Aging During Chronic HIV Infection.

People with HIV (PWH) experience an increased vulnerability to premature aging and inflammation-associated comorbidities, even when HIV replication is suppressed by antiretroviral therapy (ART). However, the factors that contribute to or are associated with this vulnerability remain uncertain. In the general population, alterations in the glycomes of circulating IgGs trigger inflammation and precede the onset of aging-associated diseases. Here, we investigate the IgG glycomes of cross-sectional and longitudinal samples from 1,216 women and men, both living with virally suppressed HIV and those without HIV. Our glycan-based machine learning models indicate that living with chronic HIV significantly accelerates the accumulation of pro-aging-associated glycomic alterations. Consistently, PWH exhibit heightened expression of senescence-associated glycan-degrading enzymes compared to their controls. These glycomic alterations correlate with elevated markers of inflammatory aging and the severity of comorbidities, potentially preceding the development of such comorbidities. Mechanistically, HIV-specific antibodies glycoengineered with these alterations exhibit reduced anti-HIV IgG-mediated innate immune functions. These findings hold significant potential for the development of glycomic-based biomarkers and tools to identify and prevent premature aging and comorbidities in people living with chronic viral infections.

Erratum: Accelerated aging with HIV begins at the time of initial HIV infection.

[This corrects the article DOI: 10.1016/j.isci.2022.104488.].

Enhancing HIV-1 latency reversal through regulating the elongating RNA Pol II pause-release by a small-molecule disruptor of PAF1C.

The polymerase-associated factor 1 complex (PAF1C) is a key, post-initiation transcriptional regulator of both promoter-proximal pausing and productive elongation catalyzed by RNA Pol II and is also involved in transcriptional repression of viral gene expression during human immunodeficiency virus-1 (HIV-1) latency. Using a molecular docking-based compound screen in silico and global sequencing-based candidate evaluation in vivo, we identified a first-in-class, small-molecule inhibitor of PAF1C (iPAF1C) that disrupts PAF1 chromatin occupancy and induces global release of promoter-proximal paused RNA Pol II into gene bodies. Transcriptomic analysis revealed that iPAF1C treatment mimics acute PAF1 subunit depletion and impairs RNA Pol II pausing at heat shock-down-regulated genes. Furthermore, iPAF1C enhances the activity of diverse HIV-1 latency reversal agents both in cell line latency models and in primary cells from persons living with HIV-1. In sum, this study demonstrates that efficient disruption of PAF1C by a first-in-class, small-molecule inhibitor may have therapeutic potential for improving current HIV-1 latency reversal strategies.

Mosaic chromosomal alterations detected in men living with HIV and the relationship to non-Hodgkin lymphoma.

OBJECTIVES: People with HIV (PWH) have an elevated risk of non-Hodgkin lymphoma (NHL) and other diseases. Studying clonal hematopoiesis (CH), the clonal expansion of mutated hematopoietic stem cells, could provide insights regarding elevated NHL risk. DESIGN: Cohort analysis of participants in the Multicenter AIDS Cohort Study ( N  = 5979). METHODS: Mosaic chromosomal alterations (mCAs), a type of CH, were detected from genotyping array data using MoChA. We compared CH prevalence in men with HIV (MWH) to HIV-uninfected men using logistic regression, and among MWH, assessed the associations of CH with NHL incidence and overall mortality using Poisson regression. RESULTS: Comparing MWH to HIV-uninfected men, we observed no difference in the frequency of autosomal mCAs (3.9% vs. 3.6%, P -value = 0.09) or mosaic loss of the Y chromosome (mLOY) (1.4% vs. 2.9%, P -value = 0.13). Autosomal mCAs involving copy-neutral loss of heterozygosity (CN-LOH) of chromosome 14q were more common in MWH. Among MWH, mCAs were not associated with subsequent NHL incidence (autosomal mCA P -value = 0.65, mLOY P -value = 0.48). However, two MWH with diffuse large B-cell lymphoma had overlapping CN-LOH mCAs on chromosome 19 spanning U2AF2 (involved in RNA splicing), and one MWH with Burkitt lymphoma had high-frequency mCAs involving chromosome 1 gain and chromosome 17 CN-LOH (cell fractions 22.1% and 25.0%, respectively). mCAs were not associated with mortality among MWH (autosomal mCA P -value = 0.52, mLOY P -value = 0.93). CONCLUSIONS: We found limited evidence for a relationship between HIV infection and mCAs. Although mCAs were not significantly associated with NHL, mCAs detected in several NHL cases indicate a need for further investigation.

No Association of IFNL4 Genotype With Opportunistic Infections and Cancers Among Men With Human Immunodeficiency Virus 1 Infection.

BACKGROUND: IFNL4 genetic variants that are strongly associated with clearance of hepatitis C virus have been linked to risk of certain opportunistic infections (OIs) and cancers, including Kaposi sarcoma, cytomegalovirus infection, and herpes simplex virus infection. As the interferon (IFN) λ family plays a role in response to viral, bacterial, and fungal infections, IFNL4 genotype might affect risk for a wide range of OIs/cancers. METHODS: We examined associations between genotype for the functional IFNL4 rs368234815 polymorphism and incidence of 16 OIs/cancers among 2310 men with human immunodeficiency virus (2038 white; 272 black) enrolled in the Multicenter AIDS Cohort Study during 1984-1990. Our primary analyses used Cox proportional hazards models adjusted for self-reported racial ancestry to estimate hazard ratios with 95% confidence intervals, comparing participants with the genotypes that generate IFN-λ4 and those with the genotype that abrogates IFN-λ4. We censored follow-up at the introduction of highly effective antiretroviral therapies. RESULTS: We found no statistically significant association between IFNL4 genotype and the incidence of Kaposi sarcoma (hazard ratio, 0.92 [95% confidence interval, .76-1.11]), cytomegalovirus infection (0.94 [.71-1.24]), herpes simplex virus infection (1.37 [.68-2.93]), or any other OI/cancer. We observed consistent results using additive genetic models and after controlling for CD4 cell count through time-dependent adjustment or restriction to participants with a low CD4 cell count. CONCLUSIONS: The absence of associations between IFNL4 genotype and these OIs/cancers provides evidence that this gene does not affect the risk of disease from opportunistic pathogens.

Accelerated aging with HIV begins at the time of initial HIV infection.

Living with HIV infection is associated with early onset of aging-related chronic conditions, sometimes described as accelerated aging. Epigenetic DNA methylation patterns can evaluate acceleration of biological age relative to chronological age. The impact of initial HIV infection on five epigenetic measures of aging was examined before and approximately 3 years after HIV infection in the same individuals (n=102). Significant epigenetic age acceleration (median 1.9-4.8 years) and estimated telomere length shortening (all p≤ 0.001) were observed from pre-to post-HIV infection, and remained significant in three epigenetic measures after controlling for T cell changes. No acceleration was seen in age- and time interval-matched HIV-uninfected controls. Changes in genome-wide co-methylation clusters were also significantly associated with initial HIV infection (p≤ 2.0 × 10-4). These longitudinal observations clearly demonstrate an early and substantial impact of HIV infection on the epigenetic aging process, and suggest a role for HIV itself in the earlier onset of clinical aging.

Gut Microbiota, Plasma Metabolomic Profiles, and Carotid Artery Atherosclerosis in HIV Infection.

BACKGROUND: Alterations in gut microbiota and blood metabolomic profiles have been implicated in HIV infection and cardiovascular disease. However, it remains unclear whether alterations in gut microbiota may contribute to disrupted host blood metabolomic profiles in relation to atherosclerosis, especially in the context of HIV infection. METHODS: We analyzed cross-sectional associations between gut microbiota features and carotid artery plaque in 361 women with or at high risk of HIV (67% HIV+), and further integrated plaque-associated microbial features with plasma lipidomic/metabolomic profiles. Furthermore, in 737 women and men, we examined prospective associations of baseline gut bacteria-associated lipidomic and metabolomic profiles with incident carotid artery plaque over 7-year follow-up. RESULTS: We found 2 potentially pathogenic bacteria, Fusobacterium and Proteus, were associated with carotid artery plaque; while the beneficial butyrate producer Odoribacter was inversely associated with plaque. Fusobacterium and Proteus were associated with multiple lipids/metabolites which were clustered into 8 modules in network. A module comprised of 9 lysophosphatidylcholines and lysophosphatidylethanolamines and a module comprised of 9 diglycerides were associated with increased risk of carotid artery plaque (risk ratio [95% CI], 1.34 [1.09-1.64] and 1.24 [1.02-1.51] per SD increment, respectively). Functional analyses identified bacterial enzymes in lipid metabolism associated with these plasma lipids. In particular, phospholipase A1 and A2 are the key enzymes in the reactions producing lysophosphatidylcholines and lysophosphatidylethanolamines. CONCLUSIONS: Among individuals with or at high risk of HIV infection, we identified altered gut microbiota and related functional capacities in the lipid metabolism associated with disrupted plasma lipidomic profiles and carotid artery atherosclerosis.

Transcriptome-wide changes in gene expression, splicing, and lncRNAs in response to a live attenuated dengue virus vaccine.

The tetravalent dengue vaccine candidate, TAK-003, induces a functional antibody response, but the titers of antibodies against the four serotypes of the dengue virus (DENV) can vary. Here, through a transcriptomic analysis on whole blood collected from recipients of a two-dose schedule of TAK-003, we examine gene expression, splicing, and transcript isoform-level changes for both protein-coding and noncoding genes to broaden our understanding of the immune response. Our analysis reveals a dynamic pattern of vaccine-associated regulation of long noncoding RNAs (lncRNAs), differential splicing of interferon-stimulated gene exons, and gene expression changes related to multiple signaling pathways that detect viral infection. Co-expression networks isolate immune cell-type-related and interferon-response modules that represent specific biological processes that correlate with more robust antibody responses. These data provide insights into the early determinants of the variable immune response to the vaccine, highlighting the significance of splicing and isoform-level gene regulatory mechanisms in defining vaccine immunogenicity.

Acrolein and other toxicant exposures in relation to cardiovascular disease among marijuana and tobacco smokers in a longitudinal cohort of HIV-positive and negative adults.

BACKGROUND: Marijuana smoke contains some of the same toxicants present in tobacco smoke. Marijuana smoking is prevalent among HIV+ individuals, but few studies have characterized smoke-related toxicants or associated health outcomes in exclusive marijuana users. METHODS: This longitudinal study included 245 participants over age 40 (76% HIV+). 33 plasma and 28 urine metabolites of nicotine, ∆-9-trans-tetrahydrocannabinol, polycyclic aromatic hydrocarbons, and volatile organic compounds were assayed by liquid or gas chromatography/mass spectrometry. Exposures and health outcomes were assessed from surveys and medical records. FINDINGS: At baseline, 18% of participants were marijuana-only smokers, 20% tobacco-only smokers, and 24% dual marijuana-tobacco smokers (median (IQR) age 53 (47-60) years, 78% male, 54% white race). Marijuana smoking was independently associated with elevated plasma naphthalenes, 2-hydroxyfluorene sulfate, 4-vinylphenol sulfate, and o-cresol sulfate (p<0·05) and urine acrylonitrile and acrylamide metabolites (p<0·05), but levels were lower than those associated with tobacco smoking. Acrolein metabolite N-Acetyl-S-(3-hydroxypropyl)-l-cysteine (3HPMA) was significantly elevated in plasma and urine in tobacco-only and dual but not marijuana-only smokers, and correlated with nicotine metabolites (p<0·05). The highest tertile of 3HPMA was associated with increased cardiovascular disease diagnoses independent of tobacco smoking, traditional risk factors, and HIV status (odds ratio [95% CI] 3·34 [1·31-8·57]; p = 0·012). INTERPRETATION: Smoke-related toxicants, including acrylonitrile and acrylamide metabolites, are detectable in exclusive marijuana smokers, but exposures are lower compared with tobacco or dual smokers. Acrolein exposure is increased by tobacco smoking but not exclusive marijuana smoking in HIV+ and HIV- adults, and contributes to cardiovascular disease in tobacco smokers. FUNDING: U.S. NIH.

Outcomes of acute hepatitis B virus (HBV) in HIV infection with and without HBV-active antiretroviral therapy.

Men with acute hepatitis B virus (HBV) infection in the Multicenter AIDS Cohort Study from 1985 to 2013 had serological testing to determine proportions with HBV recovery or chronic hepatitis B (CHB). A similar proportion of men without human immunodeficiency virus (HIV) and men with HIV receiving HBV-active antiretroviral therapy (ART) developed CHB [8.2%, 95% confidence interval (CI) 3.8-15.0% vs. 7.7%, 95% CI 2.00-36.0%]. In contrast, 17.5% (95% CI 8.7-29.9%) of men living with HIV, not on HBV-active ART developed CHB. HBV-active ART protects against developing CHB.

Increased Rate of Epigenetic Aging in Men Living With HIV Prior to Treatment.

Background: Epigenetic aging is accelerated in tissues of persons living with HIV (PLWH) and may underlie the early onset of age-related illnesses. This study examines the rate-of-change in epigenetic age in PLWH following HIV infection but before HAART, using archived longitudinal samples from the Multicenter AIDS Cohort Study. Methods: DNA was isolated from cryopreserved peripheral blood mononuclear cells from 101 men living with HIV, with baseline visit <2.5 years after HIV seroconversion (Visit 1) and follow-up visit <1.5 years before the initiation of HAART (Visit 2), and 100 HIV-uninfected men matched on age and visits with comparable time intervals. DNA methylation (DNAm) age was estimated for five clocks (Pan-tissue, Extrinsic, Phenotypic, Grim, and Skin & Blood age), and a DNAm-based estimate of telomere length (DNAmTL). Multivariate linear regression models were used to examine baseline factors associated with rate-of-aging, defined as (DNAm age visit 2-DNAm age visit 1)/(age visit 2-age visit 1). Results: Epigenetic age increased approximately twice as fast in PLWH as uninfected controls (Pan-tissue, Extrinsic, and Phenotypic clocks). Shortening of DNAmTL was nearly 3-fold faster in PLWH than controls. Faster rate-of-aging was associated with HIV status (Pan-Tissue, Extrinsic, Phenotypic, and DNAmTL), white race (Extrinsic, DNAmTL), higher cumulative HIV viral load (Grim), and lower baseline DNAm age (Phenotypic, Skin & Blood). Conclusion: Epigenetic rates-of-aging were significantly faster for untreated PLWH. Our findings expand on the important impact of HIV infection on biologic aging, both in elevating epigenetic age and increasing the rate-of-aging in the years following infection.

Small RNA sequencing of extracellular vesicles identifies circulating miRNAs related to inflammation and oxidative stress in HIV patients.

BACKGROUND: Extracellular vesicles (EVs) are nano-sized particles secreted by most cells. EVs carry nucleic acids that hold promise as potential biomarkers in various diseases. Human immunodeficiency virus type 1 (HIV) infects CD4+ T cells and induces immune dysfunction, inflammation, and EV secretion, but little is known about EV small RNA cargo in relation to immune dysregulation in HIV-infected individuals. Here, we characterize small RNA carried by circulating EVs in HIV-positive subjects on antiretroviral therapy (ART) relative to uninfected controls by next-generation RNA sequencing. RESULTS: Plasma EVs isolated from HIV-positive and HIV-negative subjects in test (n = 24) and validation (n = 16) cohorts were characterized by electron microscopy, nanoparticle tracking analysis, and immunoblotting for exosome markers. EVs were more abundant in plasma from HIV-positive compared to HIV-negative subjects. Small RNA sequencing of plasma EVs in the test cohort identified diverse small RNA species including miRNA, piRNA, snRNA, snoRNA, tRNA, and rRNA, with miRNA being the most abundant. A total of 351 different miRNAs were detected in plasma EVs, with the top 50 miRNAs accounting for 90% of all miRNA reads. miR-26a-5p was the most abundant miRNA, followed by miR-21-5p and miR-148-3p. qRT-PCR analysis showed that six miRNAs (miR-10a-5p, - 21-5p, -27b-3p, - 122-5p, -146a-5p, - 423-5p) were significantly increased in plasma EVs from HIV-positive compared to HIV-negative subjects in the validation cohort. Furthermore, miR-21-5p, -27b-3p, -146a-5p, and - 423-5p correlated positively with metabolite markers of oxidative stress and negatively with anti-inflammatory polyunsaturated fatty acids. Over-representation and pathway enrichment analyses of miRNAs and their target genes predicted functional association with oxidative stress responses, interferon gamma signaling, Toll-like receptor signaling, TGF beta signaling, and Notch signaling. CONCLUSIONS: HIV-positive individuals on ART have increased abundance of circulating EVs carrying diverse small RNAs, with miRNAs being the most abundant. Several miRNAs associated with inflammation and oxidative stress are increased in circulating EVs of HIV-positive individuals, representing potential biomarkers of targetable pathways that contribute to disease pathogenesis.

HLA tapasin independence: broader peptide repertoire and HIV control.

Human leukocyte antigen (HLA) class I allotypes vary in their ability to present peptides in the absence of tapasin, an essential component of the peptide loading complex. We quantified tapasin dependence of all allotypes that are common in European and African Americans (n = 97), which revealed a broad continuum of values. Ex vivo examination of cytotoxic T cell responses to the entire HIV-1 proteome from infected subjects indicates that tapasin-dependent allotypes present a more limited set of distinct peptides than do tapasin-independent allotypes, data supported by computational predictions. This suggests that variation in tapasin dependence may impact the strength of the immune responses by altering peptide repertoire size. In support of this model, we observed that individuals carrying HLA class I genotypes characterized by greater tapasin independence progress more slowly to AIDS and maintain lower viral loads, presumably due to increased breadth of peptide presentation. Thus, tapasin dependence level, like HLA zygosity, may serve as a means to restrict or expand breadth of the HLA-I peptide repertoire across humans, ultimately influencing immune responses to pathogens and vaccines.

Comprehensive Immunoprofiling of Pediatric Zika Reveals Key Role for Monocytes in the Acute Phase and No Effect of Prior Dengue Virus Infection.

Zika virus (ZIKV) is an emerging, mosquito-borne flavivirus responsible for recent epidemics across the Americas, and it is closely related to dengue virus (DENV). Here, we study samples from 46 DENV-naive and 43 DENV-immune patients with RT-PCR-confirmed ZIKV infection at early-acute, late-acute, and convalescent time points from our pediatric cohort study in Nicaragua. We analyze the samples via RNA sequencing (RNA-seq), CyTOF, and multiplex cytokine/chemokine Luminex to generate a comprehensive, innate immune profile during ZIKV infection. Immunophenotyping and analysis of cytokines/chemokines reveal that CD14+ monocytes play a key role during ZIKV infection. Further, we identify CD169 (Siglec-1) on CD14+ monocytes as a potential biomarker of acute ZIKV infection. Strikingly distinct transcriptomic and immunophenotypic signatures are observed at all three time points. Interestingly, pre-existing dengue immunity has minimal impact on the innate immune response to Zika. Finally, this comprehensive immune profiling and network analysis of ZIKV infection in children serves as a valuable resource.