Identification of Hub Genes and Pathways in Zika Virus Infection Using RNA-Seq Data: A Network-Based Computational Approach

Key Residue in the Precursor Region of M Protein Contributes to the Neurovirulence and Neuroinvasiveness of the African Lineage of Zika Virus

The Zika virus (ZIKV) represents an important global health threat due to its unusual association with congenital Zika syndrome. ZIKV strains are phylogenetically grouped into the African and Asian lineages. However, the viral determinants underlying the phenotypic differences between the lineages remain unknown. Here, multiple sequence alignment revealed a highly conserved residue at position 21 of the premembrane (prM) protein, which is glutamic acid and lysine in the Asian and African lineages, respectively. Using reverse genetics, we generated a recombinant virus carrying an E21K mutation based on the genomic backbone of the Asian lineage strain FSS13025 (termed E21K). The E21K mutation significantly increased viral replication in multiple neural cell lines with a higher ratio of M to prM production. Animal studies showed E21K exhibited increased neurovirulence in suckling mice, leading to more severe defects in mouse brains by causing more neural cell death and destruction of hippocampus integrity. Moreover, the E21K substitution enhanced neuroinvasiveness in interferon alpha/beta (IFN-α/β) receptor knockout mice, as indicated by the increased mortality, and enhanced replication in mouse brains. The global transcriptional analysis showed E21K infection profoundly altered neuron development networks and induced stronger antiviral immune response than wild type (WT) in both neural cells and mouse brains. More importantly, the reverse K21E mutation based on the genomic backbone of the African strain MR766 caused less mouse neurovirulence. Overall, our findings support the 21st residue of prM functions as a determinant for neurovirulence and neuroinvasiveness of the African lineage of ZIKV. IMPORTANCE The suspected link of Zika virus (ZIKV) to birth defects led the World Health Organization to declare ZIKV a Public Health Emergency of International Concern. ZIKV has been identified to have two dominant phylogenetic lineages, African and Asian. Significant differences exist between the two lineages in terms of neurovirulence and neuroinvasiveness in mice. However, the viral determinants underlying the phenotypic differences are still unknown. Here, combining reverse genetics, animal studies, and global transcriptional analysis, we provide evidence that a single E21K mutation of prM confers to the Asian lineage strain FSS130125 significantly enhanced replication in neural cell lines and more neurovirulent and neuroinvasiveness phenotypes in mice. Our findings support that the highly conserved residue at position 21 of prM functions as a determinant of neurovirulence and neuroinvasiveness of the African lineage of ZIKV in mice.

Comparison of Transcriptomic Signatures between Monkeypox-Infected Monkey and Human Cell Lines

Monkeypox virus (MPV) is a smallpox-like virus belonging to the genus Orthopoxvirus of the family Poxviridae. Unlike smallpox with no animal reservoir identified and patients suffering from milder symptoms with less mortality, several animals were confirmed to serve as natural hosts of MPV. The reemergence of a recently reported monkeypox epidemic outbreak in nonendemic countries has raised concerns about a global outburst. Since the underlying mechanism of animal-to-human transmission remains largely unknown, comprehensive analyses to discover principal differences in gene signatures during disease progression have become ever more critical. In this study, two MPV-infected in vitro models, including human immortal epithelial cancer (HeLa) cells and rhesus monkey (Macaca mulatta) kidney epithelial (MK2) cells, were chosen as the two subjects to identify alterations in gene expression profiles, together with co-regulated genes and pathways that are affected during monkeypox disease progression. Using Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and MetaCore analyses, we discovered that elevated expression of genes associated with interleukins (ILs), G protein-coupled receptors (GPCRs), heat shock proteins (HSPs), Toll-like receptors (TLRs), and metabolic-related pathways play major roles in disease progression of both monkeypox-infected monkey MK2 and human HeLa cell lines. Interestingly, our analytical results also revealed that a cluster of differentiation 40 (CD40), plasmin, and histamine served as major regulators in the monkeypox-infected monkey MK2 cell line model, while interferons (IFNs), macrophages, and neutrophil-related signaling pathways dominated the monkeypox-infected human HeLa cell line model. Among immune pathways of interest, apart from traditional monkeypox-regulated signaling pathways such as nuclear factor- (NF-κB), mitogen-activated protein kinases (MAPKs), and tumor necrosis factors (TNFs), we also identified highly significantly expressed genes in both monkey and human models that played pivotal roles during the progression of monkeypox infection, including CXCL1, TNFAIP3, BIRC3, IL6, CCL2, ZC3H12A, IL11, CSF2, LIF, PTX3, IER3, EGR1, ADORA2A, and DUOX1, together with several epigenetic regulators, such as histone cluster family gene members, HIST1H3D, HIST1H2BJ, etc. These findings might contribute to specific underlying mechanisms related to the pathophysiology and provide suggestions regarding modes of transmission, post-infectious sequelae, and vaccine development for monkeypox in the future.

Gene expression profiling of corona virus microarray datasets to identify crucial targets in COVID-19 patients

The current outbreak of coronavirus disease (COVID-19) has been affecting millions of people and has caused devastating mortality worldwide. Moreover, it is to be noted that cytokine storm has become an important cause for the rising mortality. However, the efforts for the development of drugs, vaccines and treatment has also been intervened due to poor understanding of host's defense mechanism and also due to the development of cytokine storm against this viral infection. Thus, a deeper understanding of the mechanism behind the immune dysregulation and cytokine storm development might give us clues for the clinical management of the severe cases. Hence, we have implemented differential gene expression analysis together with protein-protein interaction and Gene Ontology (GO) studies with the help of Severe Acute respiratory syndrome coronavirus (SARS-CoV) data sets such as GSE1739 and GSE33267 to give us more knowledge on the host immune response for the pathogenic coronavirus which in turn reduces the mortality. A total of 79 differentially-expressed genes (DEGs) were identified in our data set using the filters such as P-value and log2 fold change values of less than 0.05 and 1.5 respectively. Further, network analysis and GO studies showed that differential expression of two hub genes namely ELANE and LTF which could induce higher levels of pro-inflammatory cytokines in the lungs. We are certain that differential expression of ELANE and LTF results in an excessive inflammatory reaction known as the cytokine storm and ultimately leading to death. Therefore, targeting these key drivers of cytokine storm genes appears to be the potential therapeutic targets for combating the Severe Acute respiratory syndrome coronavirus - 2 (SARS-CoV-2) infection ultimately resulting in reduced mortality. Indeed, this predictive view may open new insights for designing an immune intervention for COVID-19 in the near future resulting in the mitigation of mortality rate.

The Transcriptional and Protein Profile From Human Infected Neuroprogenitor Cells Is Strongly Correlated to Zika Virus Microcephaly Cytokines Phenotype Evidencing a Persistent Inflammation in the CNS

Zika virus (ZIKV) infection during pregnancy is associated with microcephaly, a congenital malformation resulting from neuroinflammation and direct effects of virus replication on the developing central nervous system (CNS). However, the exact changes in the affected CNS remain unknown. Here, we show by transcriptome analysis (at 48 h post-infection) and multiplex immune profiling that human induced-neuroprogenitor stem cells (hiNPCs) respond to ZIKV infection with a strong induction of type-I interferons (IFNs) and several type-I IFNs stimulated genes (ISGs), notably cytokines and the pro-apoptotic chemokines CXCL9 and CXCL10. By comparing the inflammatory profile induced by a ZIKV Brazilian strain with an ancestral strain isolated from Cambodia in 2010, we observed that the response magnitude differs among them. Compared to ZIKV/Cambodia, the experimental infection of hiNPCs with ZIKV/Brazil resulted in a diminished induction of ISGs and lower induction of several cytokines (IFN-α, IL-1α/β, IL-6, IL-8, and IL-15), consequently favoring virus replication. From ZIKV-confirmed infant microcephaly cases, we detected a similar profile characterized by the presence of IFN-α, CXCL10, and CXCL9 in cerebrospinal fluid (CSF) samples collected after birth, evidencing a sustained CNS inflammation. Altogether, our data suggest that the CNS may be directly affected due to an unbalanced and chronic local inflammatory response, elicited by ZIKV infection, which contributes to damage to the fetal brain.

Assembling systems biology, embryo development and teratogenesis: What do we know so far and where to go next?

The recognition of molecular mechanisms of a teratogen can provide insights to understand its embryopathy, and later to plan strategies for the prevention of new exposures. In this context, experimental research is the most invested approach. Despite its relevance, these assays require financial and time investment. Hence, the evaluation of such mechanisms through systems biology rise as an alternative for this conventional methodology. Systems biology is an integrative field that connects experimental and computational analyses, assembling interaction networks between genes, proteins, and even teratogens. It is a valid strategy to generate new hypotheses, that can later be confirmed in experimental assays. Here, we present a literature review of the application of systems biology in embryo development and teratogenesis studies. We provide a glance at the data available in public databases, and evaluate common mechanisms between different teratogens. Finally, we discuss the advantages of using this strategy in future teratogenesis researches.

RNA-seq analysis of blood of valproic acid-responsive and non-responsive pediatric patients with epilepsy

Epilepsy is the most common chronic neurological disorder, affecting ~70 million individuals worldwide. However, approximately one-third of the patients are refractory to epilepsy medication. Of note, 100% of patients with genetic epilepsy who are resistant to the traditional drug, valproic acid (VPA), are also refractory to the other anti-epileptic drugs. The aim of the present study was to compare the transcriptomes in VPA responders and non-responders, to explore the mechanism of action of VPA and identify possible biomarkers to predict VPA resistance. Thus, RNA-seq was employed for transcriptomic analysis, differentially expressed genes (DEGs) were analyzed using Cuffdiff software and the DAVID database was used to infer the functions of the DEGs. A protein-protein interaction network was obtained using STRING and visualized with Cytoscape. A total of 389 DEGs between VPA-responsive and non-responsive pediatric patients were identified. Of these genes, 227 were upregulated and 162 were downregulated. The upregulated DEGs were largely associated with cytokines, chemokines and chemokine receptor-binding factors, whereas the downregulated DEGs were associated with cation channels, iron ion binding proteins, and immunoglobulin E receptors. In the pathway analysis, the toll-like receptor signaling pathway, pathways in cancer, and cytokine-cytokine receptor interaction were mostly enriched by the DEGs. Furthermore, three modules were identified by protein-protein interaction analysis, and the potential hub genes, chemokine (C-C motif) ligand 3 and 4, chemokine (C-X-C motif) ligand 9, tumor necrosis factor-α and interleukin-1β, which are known to be closely associated with epilepsy, were identified. These specific chemokines may participate in processes associated with VPA resistance and may be potential biomarkers for monitoring the efficacy of VPA.

Zika Virus Infection of Human Mesenchymal Stem Cells Promotes Differential Expression of Proteins Linked to Several Neurological Diseases

The recent microcephaly outbreak in Brazil has been associated with Zika virus (ZIKV) infection. The current understanding of damage caused by ZIKV infection is still unclear, since it has been implicated in other neurodegenerative and developmental complications. Here, the differential proteome analysis of human mesenchymal stem cells (hMSC) infected with a Brazilian strain of ZIKV was identified by shotgun proteomics (MudPIT). Our results indicate that ZIKV induces a potential reprogramming of the metabolic machinery in nucleotide metabolism, changes in the energy production via glycolysis and other metabolic pathways, and potentially inhibits autophagy, neurogenesis, and immune response by downregulation of signaling pathways. In addition, proteins previously described in several brain pathologies, such as Alzheimer's disease, autism spectrum disorder, amyotrophic lateral sclerosis, and Parkinson's disease, were found with altered expression due to ZIKV infection in hMSC. This potential link between ZIKV and several neuropathologies beyond microcephaly is being described here for the first time and can be used to guide specific follow-up studies concerning these specific diseases and ZIKV infection.