Next Article in Journal
Enlisting the Ixodes scapularis Embryonic ISE6 Cell Line to Investigate the Neuronal Basis of Tick—Pathogen Interactions
Next Article in Special Issue
The Impact of Intraspecies and Interspecies Bacterial Interactions on Disease Outcome
Previous Article in Journal
The Remarkable Metabolism of Vickermania ingenoplastis: Genomic Predictions
Previous Article in Special Issue
Altered Salivary Microbiome in the Early Stage of HIV Infections among Young Chinese Men Who Have Sex with Men (MSM)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

In Silico Analysis of Possible Interaction between Host Genomic Transcription Factors (TFs) and Zika Virus (ZikaSPH2015) Strain with Combinatorial Gene Regulation; Virus Versus Host—The Game Reloaded

1
U.O.C. Genetica Medica e di Laboratorio, Ospedale Antonio Cardarelli, 80131 Napoli, Italy
2
Clinica di Malattie Infettive, Dipartimento di Scienze Biomediche ed Oncologia Umana, Università degli Studi “Aldo Moro” di Bari, 70124 Bari, Italy
3
Genetica Medica, Azienda Ospedaliero Universitaria Consorziale Policlinico di Bari, 70124 Bari, Italy
*
Author to whom correspondence should be addressed.
Submission received: 5 December 2020 / Revised: 30 December 2020 / Accepted: 12 January 2021 / Published: 14 January 2021
(This article belongs to the Special Issue Microbial Interactions during Infection)

Abstract

:
In silico analysis is a promising approach for understanding biological events in complex diseases. Herein we report on the innovative computational workflow allowed to highlight new direct interactions between human transcription factors (TFs) and an entire genome of virus ZikaSPH2015 strain in order to identify the occurrence of specific motifs on a genomic Zika Virus sequence that is able to bind and, therefore, sequester host’s TFs. The analysis pipeline was performed using different bioinformatics tools available online (free of charge). According to obtained results of this in silico analysis, it is possible to hypothesize that these TFs binding motifs might be able to explain the complex and heterogeneous phenotype presentation in Zika-virus-affected fetuses/newborns, as well as the less severe condition in adults. Moreover, the proposed in silico protocol identified thirty-three different TFs identical to the distribution of TFBSs (Transcription Factor Binding Sites) on ZikaSPH2015 strain, potentially able to influence genes and pathways with biological functions confirming that this approach could find potential answers on disease pathogenesis.

1. Introduction

The Zika virus (ZIKV) is a positive-sense, single-strand RNA arbovirus belonging to the Flaviviridae family, genus Flavivirus, which includes other well known viruses, such as yellow fever, the dengue virus and the West Nile virus [1]. In 1947, ZIKV was originally isolated from a sentinel monkey, commonly found in the Zika forest in Uganda, east Africa. Shortly thereafter, other epidemiological studies evidenced the broad geographical spread of ZIKV throughout India, Southeast Asia and sub-Saharan Africa [2].
The infection causes a disease known as “Zika virus disease“, which is asymptomatic or paucisintomatic in approximately 80% of cases, and when symptomatic, adult populations experience low-grade fever, often accompanied by malaise, headache, dizziness, anorexia, conjunctivitis with retro-orbital pain, arthritis, transient arthralgia and maculopapular rash, which appear between 3 and 12 days after infection [2,3]. In some cases, neurological and autoimmune complications may occur with manifestations similar to those seen in Guillain-Barré syndrome. However, because fetal complications resulting from ZIKV infection contracted in pregnancy (first trimester) have been demonstrated [3,4], scientific interest has been increasing. In particular, Brazil was the first country where important fetal development issues were described, including microcephaly babies, with a high number of case notifications (2010: n = 153, 2011: n = 139, 2012: n = 175, 2013: n = 167 and 2014: n = 147), although high variability amongst different regions has been recorded [5].
ZIKV is spread by numerous Aedes species including Aedes aegypti, Aedes albopictus and Aedes henselli. In North America and Europe, Aedes albopictus is spreading as a new vector of ZIKV (spread from Gabon), as previously reported. The genus Aedes mosquito is common in warm climates; in the northern hemisphere it is present in North America and in different areas of the Mediterranean. Before 2007, a few cases of infections were reported at different places in Asian and African countries. Changes in the last years, due to human activities and climate changes, likely induced adaptive evolution mutation in the RNA genome of ZIKV [6]. Furthermore, during 2014–2015, ZIKV infection spread in South and Central America, the Caribbean regions and northeastern parts of Brazil. In Italy, the spread of these viruses is monitored by specific programs, as indicated in the letter of the Ministry of Health “National Surveillance Plan and response to arboviruses transmitted by mosquitoes 2018” [7].
Epidemic spread and the finding of autochthonous cases of ZIKV infection in various countries caused an international public health emergency. The health organizations therefore focused great attention on the potential risk of ZIKV transmission by transfusion of labile blood components (i.e., whole blood, packed red cells, fresh-frozen plasma, cryoprecipitate and platelets) [8]. Additionally, the transmission could also occur via breast milk, saliva and sexual intercourse. In fact, genomic ZIKV RNA was found in semen after six months in asymptomatic subjects and, in some cases, for at least ten weeks, in prostate, testicles and vaginal secretions [9].
Considering the impact of the ZIKV infection and its clinical consequences, studies are warranted to understand reasons underlining differences in clinical presentation, particularly with regards to occurrence of fetal abnormalities.
The ZIKV genome is a single-stranded of ~10.500 nt with an RNA 5′ cap structure that codes one polyprotein that contains three structural peptides (capsid protein (C), envelope (E) and membrane (prM/M) glycoproteins) and seven nonstructural proteins (NS1, NS2A, NS2B, NS3, NS4A, NS4B and NS5) [10]. The translation of ZIKV RNA begins in the cytosol immediately after release from the virion. The N-terminal part of the polyprotein contains a signal localization at the level of the endoplasmic reticulum (ER) and promotes the association of ribosomes for the viral RNA translation. This determined a particular distribution of the protein domains: PrM, E, NS1 and some extended traits of NS2A, NS4A and NS4B which are localized intraluminally, while C, NS3 and NS5 are located on the cytoplasmic side of the ER, with different transmembrane sequences pre-sent in NS2A, NS2B and NS4B [10].
Furthermore, recent discoveries revealed how viruses use nascent RNAs to benefit their own gene expression; by adopting a parasitic lifestyle, they exploit parts of the host’s transcription pathways [11]. The viruses co-evolved many pathways to transcribe their own genetic material with consequent alteration of host transcriptional machinery [11,12]. ZIKV and other flaviviruses critically depend on the cellular secretory pathways, which transfer proteins and membranes from the ER through the Golgi and to the plasma membrane, for virion assembly, maturation and release. For this reason, the recruitment of transcription factors (TFs) could take place in cytoplasm of host cells; consequently, normal processes and “trafficking” like translocation into the nucleus could be blocked or sensitively decreased.
Herein we report on the innovative computational workflow allowed to highlight proposed interaction between host TFs and viral RNA (the first complete genome of Zika virus, has been deposited in the GenBank under the accession number KU321639) [13], which, as we hypothesized, could potentially lead to an effect comparable to haploinsufficiency. The phenomenon of haploinsufficiency is defined as reduced expression levels of involved protein as a result of a mutation in one allele of the corresponding TF, leading to a phenotypic effect. Depletion of cytoplasmic TFs, following viral infection, has recently been proposed and applied on the SARS-CoV-2 strain, demonstrating how the depletion of these TFs has generated the perturbation of different metabolic pathways with biological functions, confirming that this approach, though it is an in silico study, could bring useful insights also for ZIKV infection [14].
Furthermore, it is well known that during infection the accumulation of viral proteins is accompanied by a progressive global reduction in the production of host proteins, a phenomenon that has been termed “host shutoff” [11]. The functional annotation of putative in silico target genes of these motifs (TFs) are investigated to explore and understand processes which are under the complex and heterogeneous clinical presentation in ZIKV affected fetuses/newborns, as well as less severe presentation of infection in adults.

2. Materials and Methods

The analysis pipeline was performed using different bioinformatics tools available online and consists of four main steps (Figure 1) [15].
The entire ZikaSPH2015 strain was analyzed using a MEME-ChIP that performs a comprehensive motif analysis (including motif discovery) on large sets of sequences identified by ChIP-seq or CLIP-seq experiments on Human DNA (http://meme-suite.org/tools/meme-chip) [16].
All motifs identified were used as query for Tomtom (http://meme-suite.org/doc/tomtom.html), another tool of MEME suite that compared the motifs against a database of known motifs (i.e., HOCOMOCO human v11 full). Hocomoco is a complete collection of transcription factor binding models for humans, using a large-scale ChIP-Seq analysis. Tomtom ranked the motifs in the database and produced an alignment for each significant match, searching one or more query motifs against one or more databases of human target motifs (and their reverse complements when applicable). The report for each query was a list of target motifs, ranked by p-value in the order that the queries appear in the input file. The E-value and the q-value for each match were also reported. The q-value is the minimal false discovery rate at which the observed similarity would be considered significant. Tomtom estimated q-values from all the match p-values using the Benjamini and Hochberg method. By default, significance was measured by q-value of the match [17].
For all motifs queries, a list of TFs that contained the common conserved domain was obtained.
The list of TFs was loaded in STITCH (http://stitch.embl.de/), another software tool, to identify the most related correlation within query TFs set using a guilt-by-association approach. The bioinformatics tool used a large database of functional interaction net-works from multiple organisms, and each related TF is traceable to the source network used to make the prediction [18].

3. Results

The results of our in silico analysis allowed the identification of specific motifs on the ZikaSPH2015 strain able to bind and, therefore, according to our hypothesis, to sequester thirty-three different TFs. The distribution of TFBSs (Transcription Factor Binding Sites) on ZikaSPH2015 strain is shown in the Figure 2.
Using STITCH analysis, twenty-three of those TFs were connected in a gene regulatory knowledge network, while ten seemed to be out of this network and had no interactions (Figure 3).
This in silico approach allowed us to predict recruitment of some TFs directly related to different clinical manifestation (some of them syndromic) characterized by growth retardation, dysmorphic features, intellectual disability and others.
In particular, the analysis revealed a possible reduction of BCL11A (BAF Chromatin Remodeling Complex Subunit; OMIM 606557) TFs reported in association with intellectual developmental disorder with persistence of fetal hemoglobin (HbF) including microcephaly, down-slanting palpebral fissures, strabismus and external ear abnormalities [19].
The decrease of POU1F1 (POU Domain, Class 1, Transcription Factor 1; OMIM 173110) is implicated in growth retardation combined with pituitary hormone deficiency [20], reported in different patients affected by ZIKV, same as possible subtraction of ZNF214 (Zinc Finger Protein 214; OMIM 605015) involved in Beckwith-Wiedemann syndrome-BWSCR2 (Beckwith-Wiedemann Syndrome Chromosome Region-2) characterized by mild intellectual disability, obesity and characteristic dysmorphic features [21].
Indications regarding the relationship of PRDM6 (PR domain-containing protein 6; OMIM 616982) and ZFP28 (zinc finger protein 28; OMIM 616798) were put in evidence by our approach. PRDM6 is associated with isolated nonsyndromic patent ductus arteriosus, while ZFP28 plays an essential role in controlling gene expression during cardiac and vascular pathogeneses [22,23]. This approach endorsed an important interest considering involvement of MAF (MAF bZIP transcription factor; OMIM 177075), which is associated with juvenile-onset pulverulent cataract as well as congenital cerulean cataract (Cataract Multiple Types). Beside MAF, also two other regulators of several noncrystalline human-cataract-associated genes MAFG (MAF BZIP Transcription Factor G) and MAFK (MAF BZIP Transcription Factor K) which were also evidenced [24]. These observations could explain some ocular and heart defects manifestations frequently observed in ZIKV-affected babies/patients [22]. Furthermore, MAF haploinsufficiency was associated with brachycephaly and Fine-Lubinsky syndrome characterized by psychomotor delay, brachycephaly with flat face, small nose, microstomia, cleft palate, cataract and hearing loss [25], just as with Ayme-Gripp syndrome characterized by congenital cataracts, sensorineural hearing loss, intellectual disability, seizures, brachycephaly, a distinctive flat facial appearance and reduced growth [26].
Moreover, our analysis highlighted DUX4 (Double Homeobox 4; OMIM 606009), which, when depleted, is related to muscle weakness or atrophy, one of the main characteristics of Guillaume-Barrett syndrome [27]. Indeed, KLF5 (Krüppel-like factor 5; OMIM602903) belongs to the family of Krüppel-like TFs. The deficiency of this TFs is put in relation to unusual embryonic development, cardiovascular remodeling, inflammatory stress responses and intestinal development [28].
Another group of TFs, highlighted by this approach, were those involved in immune processes as they relate to immune specific regulatory elements. In particular, ZNF680 (Zinc Finger Protein 680), IRF4 (Interferon Regulatory Factor 4; OMIM 601900), SPI1 (Spleen Focus Forming Virus Proviral Integration Oncogene; OMIM 165170) and MEIS1 (Meis Homeobox 1; OMIM 601739), which are involved in differentiation and activation of macrophages or B-cells [29,30]. In this specific group of TFs, IRF8 (Interferon Regulatory Factor-8; OMIM 601565), SPIB (SPIB Transcription Factor; OMIM 606802) and TCF7 (Transcription Factor 7; OMIM 189908) have also been noticed. IRF8 is essential for development of monocytes, Plasmacytoid dendritic cells (pDCs) and Type 1 Conventional Dendritic Cells (cDC1s) [31]. SPIB is important in production of large amounts of interferon which block viral replication [32], while TCF7, predominantly expressed in T-cells, plays a critical role in natural killer cell and innate lymphoid cell development [33]. Our workflow allowd us to to take over involvement of AR (Androgen Receptor; OMIM 313700), which, when decreased in AR knockout mice, is associated with reduction of neutrophils function [34].
NFE2L2 (Nuclear Factor Erythroid 2-Like 2; OMIM 600492), ZNF250 (Zinc Finger Protein 250) and ZNF549 (Zinc Finger Protein 549) were activated in response to infection (i.e., Herpes Simplex Virus) and inflammation. The sequestration of these TFs makes the cells more susceptible to the virus infection [35,36].
By this approach, the indications on possible recruitment of different TFs involved in embryonic development was also underlined. In detail, the following have been indicated: LYL1 (LYM-phoblastic Leuke-mia-Derived Sequence 1; OMIM 151440) [37], SNAI1 (Snail Family Transcriptional Repressor 1; OMIM 604238) [38], ZFP28 (Zinc Finger Protein 28; OMIM 616798) [23], FEZF1 (FEZ Family Zinc Finger Protein 1; OMIM 613301) [39], HOXB4 (Homeobox B4; OMIM 142965), HOXB8 (Homeobox B8; OMIM 142963) [40], TCF7L1 (Transcription Factor 7-Like 1; OMIM 604652) [41], as well as those involved in brain development such as SMCA1 (Structural Maintenance Of Chromosomes 1A; OMIM 300040) [42], TBX3 (T-Box 3; OMIM 601621) [43], ZNF323 (Zinc Finger Protein 323; OMIM 610794) [44].
Finally, considering recent emerging evidences regarding the correlation between ZIKV and diabetes mellitus, our approach explains how the reduction of some TFs can be associated to pancreatic agenesis and/or disturbing its normal function. In particular PDX1 (Pancreas/Duodenum Homeobox Protein 1; OMIM 600733) is associated with early-onset of insulin dependent diabetes mellitus [45]; NKX6-1 (NK6, Homolog of Drosophila 1O; OMIM 602563) is implicated in development of insulin producing beta cells in the islets of Langerhans [46]; and TCF7L2 (Transcription Factor 7-Like 2; OMIM 602228) is involved in blood glucose homeostasis [47].
However, the analysis performed did not revealed a direct correlation with microcephaly.

4. Discussion and Conclusions

Based on the performed in silico analysis, it was possible to identify the occurrence of specific motifs in the ZikaSPH2015 strain. These motifs could be able to bind specific host (human) TFs, and, consequently, recruitment/subtraction could be provoked. In this manner, as we hypothesized, a haploinsufficiency-like effect is occurring, with consequences on involved target genes which are predicted by our in silico approach. Even though many studies are concentrated on the molecular mechanism of ZIKV infection, the exact mechanisms that induce the complex symptomatology, including microcephaly and brain anomalies in fetuses, remain unclear. An intricate network of TFs has been suggested by our analysis (Figure 3). Interestingly, the functional annotation of putative target genes revealed evidence that TFs are specifically enriched in biological networks involved in embryonic development or other functional clusters. All of these observations lead to a new insight, which should make possible the explanation of the complex and heterogeneous phenotype expression in ZIKV affected fetus/newborns, as well as less severe clinical picture of infection in adults. Indeed, this in silico analysis showed evidence that different TFs can regulate multiple target genes, and it provides original predictive computational data, paving the way as a testable hypothesis to further studies.
The Zika virus remains ambiguous in its role in development of severe microcephaly (primary microcephaly is characterized by significant microcephaly, below-2 SD for age, usually present at birth and always present before age one year and the absence of other congenital anomalies). In fact, specific TFs directly related to microcephaly were not evidenced by our approach. Our data is suggesting, according to recent observations, that microcephaly could be a result of the combined insufficiency of several embryonic neurological developmental TFs, able to produce a synergic effect on target genes which provoke this condition [48]; or, other factors could be engaged in this phenotypic expression. Namely, ZIKV, at least some strains, seems to have a specific neurotropism, and that is why major anomalies observed in fetuses and neonates include microcephaly, ventriculomegaly, diffuse calcifications, cerebral atrophy, signs of abnormal gyration, cortical development that might lead to severe mental retardation and motor disabilities together with visual and hearing impairments [49]. However, microcephaly might not always be observed, considering that a normal head circumference was observed in 20% of ZIKV congenital syndromes as reported by França et al. [50]. Furthermore, screening for microcephaly at birth is therefore not sufficient to detect congenital syndromes [51], considering that head growth deceleration has been reported in infants born with a normal head circumference, leading to development of microcephaly after birth [52]. The first genome sequence of ZIKV, from an autochthonous transmission in Brazil, was completed in 2016, and the authors concluded that it could be used as a starting point to study genome changes and gave an explanation why microcephaly was detected in Brazil and not in other countries [10]. Furthermore, Teng et al. [53] observed that ZIKV-related proteins might overlap the microcephaly associated proteins at the P53 signaling pathway, where P53 is the hub of the genetic regulatory network, causing viral-infection-induced cell death, thereby leading to a disorder characterized by abnormal brain development. Recently Chesnut et al. [54] summarized the research on in vitro and in silico models to study mosquito-borne flavivirus neuropathogenesis, prevention and treatment. This review article reporting on an in silico genetic screening approach performed by Kumar et al. [55] was utilized to elucidate that ZIKV-elicited congenital microcephaly may act through dysregulation of the retinoic acid response element (RARE). In addition, Pylro et al. employed an in silico analysis to investigate the role of both viral and human microRNA in congenital ZIKV microcephaly [56]. Finally, the observations regarding epigenetic and epitranscriptomic gene regulation has only recently begun to emerge, and it is well established that gene expression is not determined simply by the sequence information, but rather is subject to multiple levels of control [57].
In conclusion, the suggested pipeline is easy to access, and this easily reproducible approach could be used not only to understand the mechanism of infection and therefore the characterization of TFs involved, but also to identify metabolic networks suitable to identify tissue specific biomarkers and possible pathway dysfunction, which are useful for development of a possible vaccine or therapeutic strategy [14].
The application of these integrated analysis pipelines, which can be developed in a single user-friendly interface, could help support infection investigations by providing rapid information on multiple datasets analysis (i.e., genotype and transcriptome data), especially in times of crisis such as those currently in progress with SARS-CoV-2 [14].
Certainly, some of the differences observed in the phenotypic responses may be caused by interindividual genetic variations, but the identification of these regulatory elements should increase diagnostic sensitivity and accuracy of cohort classification, and therefore guide treatment [58].

Author Contributions

Conceptualization, M.C.; Data curation, M.T., L.V. and N.B.; Formal analysis, M.C.; Project administration, N.B.; Resources, L.V.; Supervision, A.S. and N.B.; Visualization, M.T.; Writing—original draft, M.C.; Writing—review & editing, N.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

The authors would like to express deepest gratitude to Teresa Santantonio for suggestions, advices and critical reading of manuscript. The authors are grateful to Urska Gojkovic for language editing and proofreading.

Conflicts of Interest

The authors declare no conflict of interest. The authors declare that they have no competing financial and/or non-financial interests, or other interests that might be perceived to influence the results and/or discussion reported in this paper.

References

  1. Laureti, M.; Narayanan, D.; Rodriguez-Andres, J.; Fazakerley, J.K.; Kedzierski, L. Flavivirus Receptors: Diversity, Identity, and Cell Entry. Front. Immunol. 2018, 9, 2180. [Google Scholar] [CrossRef] [Green Version]
  2. Baud, D.; Gubler, D.J.; Schaub, B.; Lanteri, M.C.; Musso, D. An update on Zika virus infection. Lancet 2017, 390, 2099–2109. [Google Scholar] [CrossRef] [Green Version]
  3. Panchaud, A.; Stojanov, M.; Ammerdorffer, A.; Vouga, M.; Baud, D. Emerging Role of Zika Virus in Adverse Fetal and Neonatal Outcomes. Clin. Microbiol. Rev. 2016, 29, 659–694. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Mittal, R.; Nguyen, D.; Debs, L.H.; Patel, A.P.; Liu, G.; Jhaveri, V.M.; S Kay, S.I.; Mittal, J.; Bandstra, E.S.; Younis, R.T.; et al. Zika Virus: An Emerging Global Health Threat. Front. Cell. Infect. Microbiol. 2017, 7, 486. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. EDAC. Zika Virus Epidemic in the Americas: Potential Association with Microcephaly and Guillain-Barré Syndrome. 2015. Available online: https://www.ecdc.europa.eu/en/publications-data/rapid-risk-assessment-zika-virus-epidemic-americas-potential-association (accessed on 13 January 2021).
  6. Pyzocha, N.J.; Chinchen, S.E.; Maurer, D.M. Zany Over Zika Virus: An Overview of Diagnosis and Treatment Modalities. Curr. Sports. Med. Rep. 2017, 16, 109–113. [Google Scholar] [CrossRef] [PubMed]
  7. Rizzo, C.; Napoli, C.; Venturi, G.; Pupella, S.; Lombardini, L.; Calistri, P.; Monaco, F.; Cagarelli, R.; Angelini, P.; Bellini, R.; et al. West Nile virus transmission: Results from the integrated surveillance system in Italy, 2008 to 2015. Euro Surveill. 2016, 21, 30340. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Prisant, N.; Joguet, G.; Herrmann-Stock, C.; Moriniere, C.; Pavili, L.; Lurel, S.; Bujan, L. Upper and lower genital tract Zika virus screening in a large cohort of reproductive-age women during the Americas epidemic. Reprod. Biomed. Online 2019, 39, 624–632. [Google Scholar] [CrossRef] [PubMed]
  9. Hermanns, K.; Göhner, C.; Kopp, A.; Schmidt, A.; Merz, W.M.; Markert, U.R.; Junglen, S.; Drosten, C. Zika virus infection in human placental tissue explants is enhanced in the presence of dengue virus antibodies in-vitro. Emerg. Microbes Infect. 2018, 7, 198. [Google Scholar] [CrossRef] [Green Version]
  10. Sager, G.; Gabaglio, S.; Sztul, E.; Belov, G.A. Role of Host Cell Secretory Machinery in Zika Virus Life Cycle. Viruses 2018, 10, 559. [Google Scholar] [CrossRef] [Green Version]
  11. Harwing, A.; Landick, R.; Berkhout, B. The battle of RNA synthesis: Virus versus host. Viruses 2017, 9, 309. [Google Scholar] [CrossRef] [Green Version]
  12. Wei, H.; Zhou, M.M. Viral-encoded enzymes that target host chromatin functions. Biochim. Biophys. Acta 2010, 1799, 296–301. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Cunha, M.S.; Esposito, D.L.A.; Rocco, I.M.; Maeda, A.Y.; Vasami, F.G.S.; Nogueira, J.S.; de Souza, R.P.; Suzuki, A.; Addas-Carvalho, M.; Barjas-Castro Mde, L.; et al. First complete genome sequence of Zika virus (Flaviviridae, Flavivirus) from an autochthonous transmission in Brazil. Genome Announc. 2016, 4, e00032-16. [Google Scholar] [CrossRef] [Green Version]
  14. Chetta, M.; Rosati, A.; Marzullo, L.; Tarsitano, M.; Bukvic, N. A SARS-CoV-2 host infection model network based on genomic Human Transcription Factors (TFs) depletion. Heliyon 2020, 6, e05010. [Google Scholar] [CrossRef] [PubMed]
  15. Chetta, M.; Di Pietro, L.; Bukvic, N.; Lattanzi, W. Rising Roles of Small Noncoding RNAs in Cotranscriptional Regulation: In Silico Study of miRNA and piRNA Regulatory Network in Humans. Genes 2020, 11, 482. [Google Scholar] [CrossRef]
  16. Machanick, P.; Bailey, T.L. MEME-ChIP: Motif analysis of large DNA datasets. Bioinformatics 2011, 27, 1696–1697. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Gupta, S.; Stamatoyannopoulos, J.A.; Bailey, T.L.; Noble, W.S. Quantifying similarity between motifs. Genome Biol. 2007, 8, R24. [Google Scholar] [CrossRef] [Green Version]
  18. Kuhn, M.; von Mering, C.; Campillos, M.; Jensen, L.J.; Bork, P. STITCH: Interaction networks of chemicals and proteins. Nucleic Acids Res. 2008, 36, D684–D688. [Google Scholar] [CrossRef]
  19. Soblet, J.; Dimov, I.; Graf von Kalckreuth, C.; Cano-Chervel, J.; Baijot, S.; Pelc, K.; Sottiaux, M.; Vilain, C.; Smits, G.; Deconinck, N. BCL11A frameshift mutation associated with dyspraxia and hypotonia affecting the fine, gross, oral, and speech motor systems. Am. J. Med. Genet. A 2018, 176, 201–208. [Google Scholar] [CrossRef] [Green Version]
  20. Fang, Q.; George, A.S.; Brinkmeier, M.L.; Mortensen, A.H.; Gergics, P.; Cheung, L.Y.; Daly, A.Z.; Ajmal, A.; Pérez Millán, M.I.; Ozel, A.B.; et al. Genetics of Combined Pituitary Hormone Deficiency: Roadmap into the Genome Era. Endocr. Rev. 2016, 37, 636–675. [Google Scholar] [CrossRef] [Green Version]
  21. Sofos, E.; Pescosolido, M.F.; Quintos, J.B.; Abuelo, D.; Gunn, S.; Hovanes, K.; Morrow, E.M.; Shur, N.A. A novel familial 11p15.4 microduplication associated with intellectual disability, dysmorphic features, and obesity with involvement of the ZNF214 gene. Am. J. Med. Genet. A. 2012, 158A, 50–58. [Google Scholar] [CrossRef]
  22. Li, N.; Subrahmanyan, L.; Smith, E.; Yu, X.; Zaidi, S.; Choi, M.; Mane, S.; Nelson-Williams, C.; Behjati, M.; Kazemi, M.; et al. Mutations in the Histone Modifier PRDM6 Are Associated with Isolated Nonsyndromic Patent Ductus Arteriosus. Am. J. Hum. Genet. 2016, 98, 1082–1091. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Zhou, L.; Zhu, C.; Luo, K.; Li, Y.; Pi, H.; Yuan, W.; Wang, Y.; Huang, C.; Liu, M.; Wu, X. Identification and characterization of two novel zinc finger genes, ZNF359 and ZFP28, in human development. Biochem. Biophys. Res. Commun. 2002, 295, 862–868. [Google Scholar] [CrossRef]
  24. de Oliveira Dias, J.R.; Ventura, C.V.; de Paula Freitas, B.; Prazeres, J.; Ventura, L.O.; Bravo-Filho, V.; Aleman, T.; Ko, A.I.; Zin, A.; Belfort, R., Jr.; et al. Zika and the Eye: Pieces of a Puzzle. Prog. Retin. Eye. Res. 2018, 66, 85–106. [Google Scholar] [CrossRef] [PubMed]
  25. Aymé, S.; Philip, N. Fine-Lubinsky syndrome: A fourth patient with brachycephaly, deafness, cataract, microstomia and mental retardation. Clin. Dysmorphol. 1996, 5, 55–60. [Google Scholar] [PubMed]
  26. Niceta, M.; Stellacci, E.; Gripp, K.W.; Zampino, G.; Kousi, M.; Anselmi, M.; Traversa, A.; Ciolfi, A.; Stabley, D.; Bruselles, A.; et al. Mutations impairing GSK3-mediated MAF phosphorylation cause cataract, deafness, intellectual disability, seizures, and a Down syndrome-like facies. Am. J. Hum. Genet. 2015, 96, 816–825. [Google Scholar] [CrossRef] [Green Version]
  27. Lim, K.R.Q.; Nguyen, Q.; Yokota, T. DUX4 Signalling in the Pathogenesis of Facioscapulohumeral Muscular Dystrophy. Int. J. Mol. Sci. 2020, 21, 729. [Google Scholar] [CrossRef] [Green Version]
  28. Alaiti, M.A.; Orasanu, G.; Tugal, D.; Lu, Y.; Jain, M.K. Kruppel-like factors and vascular inflammation: Implications for atherosclerosis. Curr. Atheroscler. Rep. 2012, 14, 438–449. [Google Scholar] [CrossRef] [Green Version]
  29. Gilmour, J.; Assi, S.A.; Jaegle, U.; Kulu, D.; van de Werken, H.; Clarke, D.; Westhead, D.R.; Philipsen, S.; Bonifer, C. A crucial role for the ubiquitously expressed transcription factor Sp1 at early stages of hematopoietic specification. Development 2014, 141, 2391–2401. [Google Scholar] [CrossRef] [Green Version]
  30. Shukla, V.; Lu, R. IRF4 and IRF8: Governing the virtues of B Lymphocytes. Front. Biol. 2014, 9, 269–282. [Google Scholar] [CrossRef] [Green Version]
  31. Sichien, D.; Scott, C.L.; Martens, L.; Vanderkerken, M.; Van Gassen, S.; Plantinga, M.; Joeris, T.; De Prijck, S.; Vanhoutte, L.; Vanheerswynghels, M.; et al. IRF8 Transcription Factor Controls Survival and Function of Terminally Differentiated Conventional and Plasmacytoid Dendritic Cells, Respectively. Immunity 2016, 45, 626–640. [Google Scholar] [CrossRef] [Green Version]
  32. Cisse, B.; Caton, M.L.; Lehner, M.; Maeda, T.; Scheu, S.; Locksley, R.; Holmberg, D.; Zweier, C.; den Hollander, N.S.; Kant, S.G.; et al. Transcription factor E2-2 is an essential and specific regulator of plasmacytoid dendritic cell development. Cell 2008, 135, 37–48. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Harly, C.; Kenney, D.; Wang, Y.; Ding, Y.; Zhao, Y.; Awasthi, P.; Bhandoola, A. A Shared Regulatory Element Controls the Initiation of Tcf7 Expression During Early T Cell and Innate Lymphoid Cell Developments. Front. Immunol. 2020, 11, 470. [Google Scholar] [CrossRef] [PubMed]
  34. Lai, J.J.; Lai, K.P.; Zeng, W.; Chuang, K.H.; Altuwaijri, S.; Chang, C. Androgen receptor influences on body defense system via modulation of innate and adaptive immune systems: Lessons from conditional AR knockout mice. Am. J. Pathol. 2012, 181, 1504–1512. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Nunes Dos Santos, K.; Florentino, R.M.; França, A.; Lima Filho, A.C.M.; Santos, M.L.D.; Missiaggia, D.; Fonseca, M.C.; Brasil Costa, I.; Vidigal, P.V.T.; Nathanson, M.H.; et al. Polymorphism in the Promoter Region of NFE2L2 Gene Is a Genetic Marker of Susceptibility to Cirrhosis Associated with Alcohol Abuse. Int. J. Mol. Sci. 2019, 3589. [Google Scholar] [CrossRef] [Green Version]
  36. Johung, K.; Goodwin, E.C.; Di Maio, D. Human papillomavirus E7 repression in cervical carcinoma cells initiates a transcriptional cascade driven by the retinoblastoma family, resulting in senescence. J. Virol. 2007, 81, 2102–2116. [Google Scholar] [CrossRef] [Green Version]
  37. Capron, C.; Lécluse, Y.; Kaushik, A.L.; Foudi, A.; Lacout, C.; Sekkai, D.; Godin, I.; Albagli, O.; Poullion, I.; Svinartchouk, F.; et al. The SCL relative LYL-1 is required for fetal and adult hematopoietic stem cell function and B-cell differentiation. Blood 2006, 107, 4678–4686. [Google Scholar] [CrossRef] [Green Version]
  38. Goldman, O.; Valdes, V.J.; Ezhkova, E.; Gouon-Evans, V. The mesenchymal transcription factor SNAI-1 instructs human liver specification. Stem. Cell. Res. 2016, 17, 62–68. [Google Scholar] [CrossRef] [Green Version]
  39. Eckler, M.J.; Chen, B. Fez family transcription factors: Controlling neurogenesis and cell fate in the developing mammalian nervous system. Bioessays 2014, 36, 788–797. [Google Scholar] [CrossRef] [Green Version]
  40. Oshima, M.; Endoh, M.; Endo, T.A.; Toyoda, T.; Nakajima-Takagi, Y.; Sugiyama, F.; Koseki, H.; Kyba, M.; Iwama, A.; Osawa, M. Genome-wide analysis of target genes regulated by HoxB4 in hematopoietic stem and progenitor cells developing from embryonic stem cells. Blood 2011, 117, e142–e150. [Google Scholar] [CrossRef] [Green Version]
  41. Gaston-Massuet, C.; McCabe, M.J.; Scagliotti, V.; Young, R.M.; Carreno, G.; Gregory, L.C.; Jayakody, S.A.; Pozzi, S.; Gualtieri, A.; Basu, B.; et al. Transcription factor 7-like 1 is involved in hypothalamo-pituitary axis development in mice and humans. Proc. Natl. Acad. Sci. USA 2016, 113, E548–E557. [Google Scholar] [CrossRef] [Green Version]
  42. Karaca, E.; Harel, T.; Pehlivan, D.; Jhangiani, S.N.; Gambin, T.; Coban Akdemir, Z.; Gonzaga-Jauregui, C.; Erdin, S.; Bayram, Y.; Campbell, I.M.; et al. Genes that Affect Brain Structure and Function Identified by Rare Variant Analyses of Mendelian Neurologic Disease. Neuron 2015, 88, 499–513. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Motahari, Z.; Martinez-De Luna, R.I.; Viczian, A.S.; Zuber, M.E. Tbx3 represses bmp4 expression and, with Pax6, is required and sufficient for retina formation. Development 2016, 143, 3560–3572. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Luo, X.J.; Mattheisen, M.; Li, M.; Huang, L.; Rietschel, M.; Børglum, A.D.; Als, T.D.; van den Oord, E.J.; Aberg, K.A.; Mors, O.; et al. Systematic Integration of Brain eQTL and GWAS Identifies ZNF323 as a Novel Schizophrenia Risk Gene and Suggests Recent Positive Selection Based on Compensatory Advantage on Pulmonary Function. Schizophr. Bull. 2015, 41, 1294–1308. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Nicolino, M.; Claiborn, K.C.; Senée, V.; Boland, A.; Stoffers, D.A.; Julier, C. A novel hypomorphic PDX1 mutation responsible for permanent neonatal diabetes with subclinical exocrine deficiency. Diabetes 2010, 59, 733–740. [Google Scholar] [CrossRef] [Green Version]
  46. Bansal, V.S.; Raja, C.P.; Venkataraman, K.; Vijayalakshmi, M.A. Genes involved in pancreatic islet cell rejuvenation. Indian J Med. Res. 2013, 137, 695–703. [Google Scholar]
  47. Zhuang, Y.; Niu, F.; Liu, D.; Sun, J.; Zhang, X.; Zhang, J.; Guo, S. Associations of TCF7L2 gene polymorphisms with the risk of diabetic nephropathy: A case-control study. Medicine (Baltimore) 2018, 97, e8388. [Google Scholar] [CrossRef]
  48. Lengner, C.J.; Welstead, G.G.; Jaenisch, R. The pluripotency regulator Oct4: A role in somatic stem cells? Cell Cycle 2008, 7, 725–728. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Vouga, M.; Baud, D. Imaging of congenital Zika virus infection: The route to identification of prognostic factors. Prenat. Diagn. 2016, 36, 799–811. [Google Scholar] [CrossRef]
  50. França, G.V.A.; Schuler-Faccini, L.; Oliveira, W.K.; Henriques, C.M.; Carmo, E.H.; Pedi, V.D.; Nunes, M.L.; Castro, M.C.; Serruya, S.; Silveira, M.F.; et al. Congenital Zika virus syndrome in Brazil: A case series of the first 1501 livebirths with complete investigation. Lancet 2016, 388, 891–897. [Google Scholar] [CrossRef] [Green Version]
  51. Eppes, C.; Rac, M.; Dunn, J.; Versalovic, J.; Murray, K.O.; Suter, M.A.; Sanz Cortes, M.; Espinoza, J.; Seferovic, M.D.; Lee, W.; et al. Testing for Zika virus infection in pregnancy: Key concepts to deal with an emerging epidemic. Am. J. Obstet. Gynecol. 2017, 216, 209–225. [Google Scholar] [CrossRef] [Green Version]
  52. van der Linden, V.; Pessoa, A.; Dobyns, W.; Barkovich, A.J.; Júnior, H.V.; Filho, E.L.; Ribeiro, E.M.; Leal, M.C.; Coimbra, P.P.; Aragão, M.F.; et al. Description of 13 infants born during October 2015–January 2016 with congenital Zika virus infection without microcephaly at birth Brazil. Morb. Mortal. Wkly. Rep. 2016, 65, 1343–1348. [Google Scholar] [CrossRef] [PubMed]
  53. Teng, Y.; Liu, S.; Guo, X.; Liu, S.; Jin, Y.; He, T.; Bi, D.; Zhang, P.; Lin, B.; An, X.; et al. An Integrative Analysis Reveals a Central Role of P53 Activation via MDM2 in Zika Virus Infection Induced Cell Death. Front. Cell. Infect. Microbiol. 2017, 7, 327. [Google Scholar] [CrossRef] [PubMed]
  54. Chesnut, M.; Muñoz, L.S.; Harris, G.; Freeman, D.; Gama, L.; Pardo, C.A.; Pamies, D. In vitro and in silico Models to Study Mosquito-Borne Flavivirus Neuropathogenesis. Prevention, and Treatment. Front. Cell. Infect. Microbiol. 2019, 9, 223. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Kumar, A.; Singh, H.N.; Pareek, V.; Raza, K.; Dantham, S.; Kumar, P.; Mochan, S.; Faiq, M.A. A possible mechanism of zika virus associated microcephaly: Imperative role of retinoic acid response element (RARE) consensus sequence repeats in the viral genome. Front. Hum. Neurosci. 2016, 10, 403. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Pylro, V.S.; Oliveira, F.S.; Morais, D.K.; Cuadros-Orellana, S.; Pais, F.S.; Medeiros, J.D.; Geraldo, J.A.; Gilbert, J.; Volpini, A.C.; Fernandes, G.R. ZIKV–CDB: A collaborative database to guide research linking SncRNAs and ZIKA virus disease symptoms. PLoS Negl. Trop. Dis. 2016, 10, e0004817. [Google Scholar] [CrossRef]
  57. Tsai, K.; Cullen, B.R. Epigenetic and epitranscriptomic regulation of viral replication. Nat. Rev. Microbiol. 2020, 18, 559–570. [Google Scholar] [CrossRef]
  58. Kalayci, S.; Selvan, M.E.; Ramos, I.; Cotsapas, C.; Harris, E.; Kim, E.Y.; Montgomery, R.R.; Poland, G.; Pulendran, B.; Tsang, J.; et al. ImmuneRegulation: A web-based tool for identifying human immune regulatory elements. Nucleic Acids Res. 2019, 47, 142–150. [Google Scholar] [CrossRef] [Green Version]
Figure 1. A schematic representation of the pipeline integrating the different analytical software. The details of the analytcal steps are described in the text (Materials and Methods section).
Figure 1. A schematic representation of the pipeline integrating the different analytical software. The details of the analytcal steps are described in the text (Materials and Methods section).
Pathogens 10 00069 g001
Figure 2. The distribution of TFBSs (transcription factor binding sites) on the ZikaSPH2015 strain.
Figure 2. The distribution of TFBSs (transcription factor binding sites) on the ZikaSPH2015 strain.
Pathogens 10 00069 g002
Figure 3. The scheme provides an overview of the TFs’ computational relationship analysis (based on STITCH tool) associated with each functional cluster.
Figure 3. The scheme provides an overview of the TFs’ computational relationship analysis (based on STITCH tool) associated with each functional cluster.
Pathogens 10 00069 g003
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chetta, M.; Tarsitano, M.; Vicari, L.; Saracino, A.; Bukvic, N. In Silico Analysis of Possible Interaction between Host Genomic Transcription Factors (TFs) and Zika Virus (ZikaSPH2015) Strain with Combinatorial Gene Regulation; Virus Versus Host—The Game Reloaded. Pathogens 2021, 10, 69. https://0-doi-org.brum.beds.ac.uk/10.3390/pathogens10010069

AMA Style

Chetta M, Tarsitano M, Vicari L, Saracino A, Bukvic N. In Silico Analysis of Possible Interaction between Host Genomic Transcription Factors (TFs) and Zika Virus (ZikaSPH2015) Strain with Combinatorial Gene Regulation; Virus Versus Host—The Game Reloaded. Pathogens. 2021; 10(1):69. https://0-doi-org.brum.beds.ac.uk/10.3390/pathogens10010069

Chicago/Turabian Style

Chetta, Massimiliano, Marina Tarsitano, Laura Vicari, Annalisa Saracino, and Nenad Bukvic. 2021. "In Silico Analysis of Possible Interaction between Host Genomic Transcription Factors (TFs) and Zika Virus (ZikaSPH2015) Strain with Combinatorial Gene Regulation; Virus Versus Host—The Game Reloaded" Pathogens 10, no. 1: 69. https://0-doi-org.brum.beds.ac.uk/10.3390/pathogens10010069

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop