Rapid adaptation to elevated extracellular potassium in the pyloric circuit of the crab, Cancer borealis.

Elevated potassium concentration ([K+]) is often used to alter excitability in neurons and networks by shifting the potassium equilibrium potential (EK) and, consequently, the resting membrane potential. We studied the effects of increased extracellular [K+] on the well-described pyloric circuit of the crab Cancer borealis. A 2.5-fold increase in extracellular [K+] (2.5×[K+]) depolarized pyloric dilator (PD) neurons and resulted in short-term loss of their normal bursting activity. This period of silence was followed within 5-10 min by the recovery of spiking and/or bursting activity during continued superfusion of 2.5×[K+] saline. In contrast, when PD neurons were pharmacologically isolated from pyloric presynaptic inputs, they exhibited no transient loss of spiking activity in 2.5×[K+], suggesting the presence of an acute inhibitory effect mediated by circuit interactions. Action potential threshold in PD neurons hyperpolarized during an hour-long exposure to 2.5×[K+] concurrent with the recovery of spiking and/or bursting activity. Thus the initial loss of activity appears to be mediated by synaptic interactions within the network, but the secondary adaptation depends on changes in the intrinsic excitability of the pacemaker neurons. The complex sequence of events in the responses of pyloric neurons to elevated [K+] demonstrates that electrophysiological recordings are necessary to determine both the transient and longer term effects of even modest alterations of K+ concentrations on neuronal activity.NEW & NOTEWORTHY Solutions with elevated extracellular potassium are commonly used as a depolarizing stimulus. We studied the effects of high potassium concentration ([K+]) on the pyloric circuit of the crab stomatogastric ganglion. A 2.5-fold increase in extracellular [K+] caused a transient loss of activity that was not due to depolarization block, followed by a rapid increase in excitability and recovery of spiking within minutes. This suggests that changing extracellular potassium can have complex and nonstationary effects on neuronal circuits.

Positional Correlation Natural Vector: A Novel Method for Genome Comparison.

Advances in sequencing technology have made large amounts of biological data available. Evolutionary analysis of data such as DNA sequences is highly important in biological studies. As alignment methods are ineffective for analyzing large-scale data due to their inherently high costs, alignment-free methods have recently attracted attention in the field of bioinformatics. In this paper, we introduce a new positional correlation natural vector (PCNV) method that involves converting a DNA sequence into an 18-dimensional numerical feature vector. Using frequency and position correlation to represent the nucleotide distribution, it is possible to obtain a PCNV for a DNA sequence. This new numerical vector design uses six suitable features to characterize the correlation among nucleotide positions in sequences. PCNV is also very easy to compute and can be used for rapid genome comparison. To test our novel method, we performed phylogenetic analysis with several viral and bacterial genome datasets with PCNV. For comparison, an alignment-based method, Bayesian inference, and two alignment-free methods, feature frequency profile and natural vector, were performed using the same datasets. We found that the PCNV technique is fast and accurate when used for phylogenetic analysis and classification of viruses and bacteria.

A novel alignment-free vector method to cluster protein sequences.

Classification of protein are crucial topics in biology. The number of protein sequences stored in databases increases sharply in the past decade. Traditionally, comparison of protein sequences is usually carried out through multiple sequence alignment methods. However, these methods may be unsuitable for clustering of protein sequences when gene rearrangements occur such as in viral genomes. The computation is also very time-consuming for large datasets with long genomes. In this paper, based on three important biochemical properties of amino acids: the hydropathy index, polar requirement and chemical composition of the side chain, we propose a 24 dimensional feature vector describing the composition of amino acids in protein sequences. Our method not only utilizes the chemical properties of amino acids but also counts on their numbers and positions. The results on beta-globin, mammals, and three virus datasets show that this new tool is fast and accurate for classifying proteins and inferring the phylogeny of organisms.

Effect of prostaglandin E2 on multidrug resistance transporters in human placental cells.

Prostaglandin (PG) E2, a major product of cyclooxygenase (COX)-2, acts as an immunomodulator at the maternal-fetal interface during pregnancy. It exerts biologic function through interaction with E-prostanoid (EP) receptors localized to the placenta. The activation of the COX-2/PGE2/EP signal pathway can alter the expression of the ATP-binding cassette (ABC) transporters, multidrug resistance protein 1 [P-glycoprotein (Pgp); gene: ABCB1], and breast cancer resistance protein (BCRP; gene: ABCG2), which function to extrude drugs and xenobiotics from cells. In the placenta, PGE2-mediated changes in ABC transporter expression could impact fetal drug exposure. Furthermore, understanding the signaling cascades involved could lead to strategies for the control of Pgp and BCRP expression levels. We sought to determine the impact of PGE2 signaling mechanisms on Pgp and BCRP in human placental cells. The treatment of placental cells with PGE2 up-regulated BCRP expression and resulted in decreased cellular accumulation of the fluorescent substrate Hoechst 33342. Inhibiting the EP1 and EP3 receptors with specific antagonists attenuated the increase in BCRP. EP receptor signaling results in activation of transcription factors, which can affect BCRP expression. Although PGE2 decreased nuclear factor κ-light chain-enhancer of activated B activation and increased activator protein 1, chemical inhibition of these inflammatory transcription factors did not blunt BCRP up-regulation by PGE2. Though PGE2 decreased Pgp mRNA, Pgp expression and function were not significantly altered. Overall, these findings suggest a possible role for PGE2 in the up-regulation of placental BCRP expression via EP1 and EP3 receptor signaling cascades.

Small Heterodimer Partner Modulates Macrophage Differentiation during Innate Immune Response through the Regulation of Peroxisome Proliferator Activated Receptor Gamma, Mitogen-Activated Protein Kinase, and Nuclear Factor Kappa B Pathways.

Hepatic macrophages act as the liver's first line of defense against injury. Their differentiation into proinflammatory or anti-inflammatory subpopulations is a critical event that maintains a delicate balance between liver injury and repair. In our investigation, we explored the influence of the small heterodimer partner (SHP), a nuclear receptor primarily associated with metabolism, on macrophage differentiation during the innate immune response. During macrophage differentiation, we observed significant alterations in Shp mRNA expression. Deletion of Shp promoted M1 differentiation while interfering with M2 polarization. Conversely, overexpression of SHP resulted in increased expression of peroxisome proliferator activated receptor gamma (Pparg), a master regulator of anti-inflammatory macrophage differentiation, thereby inhibiting M1 differentiation. Upon lipopolysaccharide (LPS) injection, there was a notable increase in the proinflammatory M1-like macrophages, accompanied by exacerbated infiltration of monocyte-derived macrophages (MDMs) into the livers of Shp myeloid cell specific knockout (Shp-MKO). Concurrently, we observed significant induction of tumor necrosis factor alpha (Tnfa) and chemokine (C-C motif) ligand 2 (Ccl2) expression in LPS-treated Shp-MKO livers. Additionally, the mitogen-activated protein kinase (MAPK) and nuclear factor kappa B (NF-κB) pathways were activated in LPS-treated Shp-MKO livers. Consistently, both pathways were hindered in SHP overexpression macrophages. Finally, we demonstrated that SHP interacts with p65, thereby influencing macrophage immune repones. In summary, our study uncovered a previously unrecognized role of SHP in promoting anti-inflammatory macrophage differentiation during the innate immune response. This was achieved by SHP acting as a regulator for the Pparg, MAPK, and NF-κB pathways.

Functional inhibition of the RNA-binding protein HuR sensitizes triple-negative breast cancer to chemotherapy.

Chemotherapy remains the standard treatment for triple-negative breast cancer (TNBC); however, chemoresistance compromises its efficacy. The RNA-binding protein Hu antigen R (HuR) could be a potential therapeutic target to enhance the chemotherapy efficacy. HuR is known to mainly stabilize its target mRNAs, and/or promote the translation of encoded proteins, which are implicated in multiple cancer hallmarks, including chemoresistance. In this study, a docetaxel-resistant cell subline (231-TR) was established from the human TNBC cell line MDA-MB-231. Both the parental and resistant cell lines exhibited similar sensitivity to the small molecule functional inhibitor of HuR, KH-3. Docetaxel and KH-3 combination therapy synergistically inhibited cell proliferation in TNBC cells and tumor growth in three animal models. KH-3 downregulated the expression levels of HuR targets (e.g., ß-Catenin and BCL2) in a time- and dose-dependent manner. Moreover, KH-3 restored docetaxel's effects on activating Caspase-3 and cleaving PARP in 231-TR cells, induced apoptotic cell death, and caused S-phase cell cycle arrest. Together, our findings suggest that HuR is a critical mediator of docetaxel resistance and provide a rationale for combining HuR inhibitors and chemotherapeutic agents to enhance chemotherapy efficacy.

Molecular tuning of sea anemone stinging.

Jellyfish and sea anemones fire single-use, venom-covered barbs to immobilize prey or predators. We previously showed that the anemone Nematostella vectensis uses a specialized voltage-gated calcium (CaV) channel to trigger stinging in response to synergistic prey-derived chemicals and touch (Weir et al., 2020). Here, we use experiments and theory to find that stinging behavior is suited to distinct ecological niches. We find that the burrowing anemone Nematostella uses uniquely strong CaV inactivation for precise control of predatory stinging. In contrast, the related anemone Exaiptasia diaphana inhabits exposed environments to support photosynthetic endosymbionts. Consistent with its niche, Exaiptasia indiscriminately stings for defense and expresses a CaV splice variant that confers weak inactivation. Chimeric analyses reveal that CaVß subunit adaptations regulate inactivation, suggesting an evolutionary tuning mechanism for stinging behavior. These findings demonstrate how functional specialization of ion channel structure contributes to distinct organismal behavior.

Niclosamide improves cancer immunotherapy by modulating RNA-binding protein HuR-mediated PD-L1 signaling.

BACKGROUND: Immune checkpoint blockade (ICB) represents a revolutionary advance in cancer treatment but remains limited success in triple-negative breast cancer (TNBC). Here we aim to explore the mechanism of RNA-binding protein (RBP) HuR in cancer immune evasion by post-transcriptionally regulating PD-L1 and evaluate the potential of HuR inhibition to improve immune response. METHODS: The binding between HuR and PD-L1 mRNA was determined by ribonucleoprotein immunoprecipitation and RNA pull-down assays. The HuR knockout clones were established by CRISPR/Cas9 technology. The protein levels were assessed by Western blot, immunohistochemistry, and immunocytochemistry. The function and molecular mechanism of HuR-PD-L1 were determined by in vitro T cell activation and killing assay and in vivo efficacy assay. RESULTS: We found that HuR directly bound to and stabilized PD-L1 mRNA. Knocking out HuR reduced PD-L1 levels and promoted T cell activation. We discovered that niclosamide reduced PD-L1 by inhibiting HuR cytoplasmic translocation, and diminished glycosylation of PD-L1. Niclosamide enhanced T cell-mediated killing of cancer cells and significantly improved the efficacy of anti-PD-1 immunotherapy in two syngeneic animal tumor models. CONCLUSION: We identified HuR as a novel posttranscriptional regulator of PD-L1, which plays an important role in tumor immune evasion. Niclosamide might be a promising repurposed drug to improve the patient response to immunotherapy by targeting HuR-PD-L1 axis. Our study demonstrates a novel strategy for targeting HuR/PD-L1 and provides the first proof-of-principle for repurposing niclosamide as a HuR inhibitor to overcome cancer immune evasion and improve response to ICB immunotherapy.

Molecular tuning of sea anemone stinging.

Jellyfish and sea anemones fire single-use, venom-covered barbs to immobilize prey or predators. We previously showed that the anemone Nematostella vectensis uses a specialized voltage-gated calcium (CaV) channel to trigger stinging in response to synergistic prey-derived chemicals and touch (Weir et al., 2020). Here we use experiments and theory to find that stinging behavior is suited to distinct ecological niches. We find that the burrowing anemone Nematostella uses uniquely strong CaV inactivation for precise control of predatory stinging. In contrast, the related anemone Exaiptasia diaphana inhabits exposed environments to support photosynthetic endosymbionts. Consistent with its niche, Exaiptasia indiscriminately stings for defense and expresses a CaV splice variant that confers weak inactivation. Chimeric analyses reveal that CaVß subunit adaptations regulate inactivation, suggesting an evolutionary tuning mechanism for stinging behavior. These findings demonstrate how functional specialization of ion channel structure contributes to distinct organismal behavior.

Hepatocyte-specific deletion of small heterodimer partner protects mice against acetaminophen-induced hepatotoxicity via activation of Nrf2.

Acetaminophen (APAP) overdose stands as the primary cause of acute liver failure in the United States. APAP hepatotoxicity involves hepatic glutathione (GSH) depletion and mitochondrial damage. To counteract the toxicity of APAP, the nuclear factor erythroid 2 like 2 (Nrf2) activates the expression of genes responsible for drug detoxification and GSH synthesis. In this study, we present evidence that the elimination of hepatocyte small heterodimer partner, a critical transcriptional repressor for liver metabolism, results in Nrf2 activation and protects mice from APAP-induced acute liver injury. Initial investigations conducted on wildtype (WT) mice revealed a swift downregulation of Shp mRNA within the first 24 h after APAP administration. Subsequent treatment of hepatocyte-specific Shp knockout (ShpHep-/-) mice with 300 mg/kg APAP for 2 h exhibited comparable bioactivation of APAP with that observed in the WT controls. However, a significant reduction in liver injury was observed in ShpHep-/- after APAP treatment for 6 and 24 h. The decreased liver injury correlated with a faster recovery of GSH, attributable to heightened expression of Nrf2 target genes involved in APAP detoxification and GSH synthesis. Moreover, in vitro studies revealed that SHP protein interacted with NRF2 protein, inhibiting the transcription of Nrf2 target genes. These findings hold relevance for humans, as overexpression of SHP hindered APAP-induced NRF2 activation in primary human hepatocytes. In conclusion, our studies have unveiled a novel regulatory axis involving SHP and NRF2 in APAP-induced acute liver injury, emphasizing SHP as a promising therapeutic target in APAP overdose-induced hepatotoxicity.

Hepatic transcriptome profiling reveals early signatures associated with disease transition from non-alcoholic steatosis to steatohepatitis.

Background and aim: Non-alcoholic fatty liver disease (NAFLD) is becoming a leading cause of chronic liver disease worldwide. The molecular events that influence disease progression from non-alcoholic fatty liver (NAFL) to aggressive non-alcoholic steatohepatitis (NASH) remain incompletely understood, leading to lack of mechanism-based targeted treatment options for NASH. This study aims to identify early signatures associated with disease progression from NAFL to NASH in mice and humans. Materials and methods: Male C57BL/6J mice were fed a high-fat, -cholesterol, and - fructose (HFCF) diet for up to 9 months. The extent of steatosis, inflammation, and fibrosis was evaluated in liver tissues. Total RNA sequencing (RNA-seq) was conducted to determine liver transcriptomic changes. Results: After being fed the HFCF diet, mice sequentially developed steatosis, early steatohepatitis, steatohepatitis with fibrosis, and eventually spontaneous liver tumor. Hepatic RNA-seq revealed that the key signatures during steatosis progression to early steatohepatitis were pathways related to extracellular matrix organization and immune responses such as T cell migration, arginine biosynthesis, C-type lectin receptor signaling, and cytokine-cytokine receptor interaction. Genes regulated by transcription factors forkhead box M1 (FOXM1) and negative elongation factor complex member E (NELFE) were significantly altered during disease progression. This phenomenon was also observed in patients with NASH. Conclusions: In summary, we identified early signatures associated with disease progression from NAFL to early NASH in a mouse model that recapitulated key metabolic, histologic, and transcriptomic changes seen in humans. The findings from our study may shed light on the development of novel preventative, diagnostic, and therapeutic strategies for NASH.

Alignment-free sequence comparison for virus genomes based on location correlation coefficient.

Coronaviruses (especially SARS-CoV-2) are characterized by rapid mutation and wide spread. As these characteristics easily lead to global pandemics, studying the evolutionary relationship between viruses is essential for clinical diagnosis. DNA sequencing has played an important role in evolutionary analysis. Recent alignment-free methods can overcome the problems of traditional alignment-based methods, which consume both time and space. This paper proposes a novel alignment-free method called the correlation coefficient feature vector (CCFV), which defines a correlation measure of the L-step delay of a nucleotide location from its location in the original DNA sequence. The numerical feature is a 16×L-dimensional numerical vector describing the distribution characteristics of the nucleotide positions in a DNA sequence. The proposed L-step delay correlation measure is interestingly related to some types of L+1 spaced mers. Unlike traditional gene comparison, our method avoids the computational complexity of multiple sequence alignment, and hence improves the speed of sequence comparison. Our method is applied to evolutionary analysis of the common human viruses including SARS-CoV-2, Dengue virus, Hepatitis B virus, and human rhinovirus and achieves the same or even better results than alignment-based methods. Especially for SARS-CoV-2, our method also confirms that bats are potential intermediate hosts of SARS-CoV-2.

Geometric construction of viral genome space and its applications.

Understanding the relationships between genomic sequences is essential to the classification and characterization of living beings. The classes and characteristics of an organism can be identified in the corresponding genome space. In the genome space, the natural metric is important to describe the distribution of genomes. Therefore, the similarity of two biological sequences can be measured. Here, we report that all of the viral genomes are in 32-dimensional Euclidean space, in which the natural metric is the weighted summation of Euclidean distance of k-mer natural vectors. The classification of viral genomes in the constructed genome space further proves the convex hull principle of taxonomy, which states that convex hulls of different families are mutually disjoint. This study provides a novel geometric perspective to describe the genome sequences.

A novel alignment-free method for HIV-1 subtype classification.

HIV-1 is the most common and pathogenic strain of human immunodeficiency virus consisting of many subtypes. To study the difference among HIV-1 subtypes in infection, diagnosis and drug design, it is important to identify HIV-1 subtypes from clinical HIV-1 samples. In this work, we propose an effective numeric representation called Subsequence Natural Vector (SNV) to encode HIV-1 sequences. Using the representation, we introduce an improved linear discriminant analysis method to classify HIV-1 viruses correctly. SNV is based on distribution of nucleotides in HIV-1 viral sequences. It not only computes the number of nucleotides, but also describes the position and variance of nucleotides in viruses. To validate our alignment-free method, 6902 complete genomes and 11,668 pol gene sequences of HIV-1 subtypes were collected from the up-to-date Los Alamos HIV database. SNV outperforms the three popular methods, Kameris, Comet and REGA, with almost 100% Sensitivity and Specificity, also with much less time. Our subtyping algorithm especially works better for circulating recombinant forms (CRFs) consisting of a few sequences. Our approach is also powerful to separate unique recombinant forms (URFs) from other subtypes with 100% Sensitivity and Specificity. Moreover, phylogenetic trees based on SNV representation are constructed using full-length HIV-1 genomes and pol genes respectively, where viruses from the same subtype are clustered together correctly.

A Novel Approach to Clustering Genome Sequences Using Inter-nucleotide Covariance.

Classification of DNA sequences is an important issue in the bioinformatics study, yet most existing methods for phylogenetic analysis including Multiple Sequence Alignment (MSA) are time-consuming and computationally expensive. The alignment-free methods are popular nowadays, whereas the manual intervention in those methods usually decreases the accuracy. Also, the interactions among nucleotides are neglected in most methods. Here we propose a new Accumulated Natural Vector (ANV) method which represents each DNA sequence by a point in ℝ18. By calculating the Accumulated Indicator Functions of nucleotides, we can further find an Accumulated Natural Vector for each sequence. This new Accumulated Natural Vector not only can capture the distribution of each nucleotide, but also provide the covariance among nucleotides. Thus global comparison of DNA sequences or genomes can be done easily in ℝ18. The tests of ANV of datasets of different sizes and types have proved the accuracy and time-efficiency of the new proposed ANV method.

Expression, Regulation, and Function of the Calmodulin Accessory Protein PCP4/PEP-19 in Myometrium.

OBJECTIVE: Calmodulin (CaM) plays a key role in the orchestration of Ca2+ signaling events, and its regulation is considered an important component of cellular homeostasis. The control of uterine smooth muscle function is largely dependent on the regulation of Ca2+ and CaM signaling. The objective of this study was to investigate the expression, function, and regulation of CaM regulatory proteins in myometrium during pregnancy. STUDY DESIGN: Myometrium was obtained from nonpregnant women and 4 groups of pregnant women at the time their primary cesarean delivery: (i) preterm not in labor, (ii) preterm in labor with clinical and/or histological diagnosis of chorioamnionitis, (3) term not in labor; and (4) term in labor. The effect of perinatal inflammation on pcp4/pep-19 expression was evaluated in a mouse model of Ureaplasma parvum-induced chorioamnionitis. Human myometrial cells stably expressing wild-type and mutant forms of PCP4/PEP-19 were used in the evaluation of agonist-induced intracellular Ca2+ mobilization. RESULTS: Compared to other CaM regulatory proteins, PCP4/PEP-19 transcripts were more abundant in human myometrium. The expression of PCP4/PEP-19 was lowest in myometrium of women with preterm pregnancy and chorioamnionitis. In the mouse uterus, pcp4/pep-19 expression was lower in late compared to mid-gestation and decreased in mice injected intra-amniotic with Ureaplasma parvum. In myometrial smooth muscle cells, tumor necrosis factor alpha and progesterone decreased and PCP4/PEP-19 promoter activity increased. Finally, the overexpression of PCP4/PEP-19 reduced agonist-induced intracellular Ca2+ levels in myometrial cells. CONCLUSION: The decreased expression of PCP4/PEP-19 in myometrium contributes to a loss of quiescence in response to infection-induced inflammation at preterm pregnancy.

Regulation of human placental drug transporters in HCV infection and their influence on direct acting antiviral medications.

INTRODUCTION: The objectives of this study were to determine how HCV infection affects placental drug transporters, and to determine the role of drug transporters on the cellular accumulation of direct-acting antiviral drugs in human trophoblasts. METHODS: Eighty-four ABC and SLC transporter genes were first screened in normal and HCV infected pregnant women using PCR profiler array. The changes in expression were confirmed by qPCR and Western blot. The impact of selected drug transporters on the cellular accumulation of radiolabeled antiviral drugs sofosbuvir, entecavir, and tenofovir was measured in primary human trophoblasts (PHT) and BeWo b30 cells in the presence or absence of transporter-specific inhibitors. PHT were then treated with CL097, ssRNA40, and imquimod to determine the impact of Toll-like receptor (TLR) 7/8 activation on drug transporter expression. RESULTS: The expression of the ABC efflux transporters ABCB1/P-gp and ABCG2/BCRP was increased in placenta of women with HCV, while the nucleoside transporters SLC29A1/ENT1 and SLC29A2/ENT2 remained unchanged. The accumulation of sofosbuvir and tenofovir was unaffected by inhibition of these transporters in trophoblast cells. Entecavir accumulation was decreased by the inhibition of ENT2. P-gp and BCRP inhibition enhanced entecavir accumulation in BeWo b30, but not PHT. Overall, there was little effect of TLR7/8 activation on these drug transporters, and the accumulation of entecavir in PHT. DISCUSSION: The data suggest that expression of placental drug transporters and selection of antiviral drug may impact fetal drug exposure in pregnancies complicated by HCV infections.

A novel fast vector method for genetic sequence comparison.

With sharp increasing in biological sequences, the traditional sequence alignment methods become unsuitable and infeasible. It motivates a surge of fast alignment-free techniques for sequence analysis. Among these methods, many sorts of feature vector methods are established and applied to reconstruction of species phylogeny. The vectors basically consist of some typical numerical features for certain biological problems. The features may come from the primary sequences, secondary or three dimensional structures of macromolecules. In this study, we propose a novel numerical vector based on only primary sequences of organism to build their phylogeny. Three chemical and physical properties of primary sequences: purine, pyrimidine and keto are also incorporated to the vector. Using each property, we convert the nucleotide sequence into a new sequence consisting of only two kinds of letters. Therefore, three sequences are constructed according to the three properties. For each letter of each sequence we calculate the number of the letter, the average position of the letter and the variation of the position of the letter appearing in the sequence. Tested on several datasets related to mammals, viruses and bacteria, this new tool is fast in speed and accurate for inferring the phylogeny of organisms.

Zika and Flaviviruses Phylogeny Based on the Alignment-Free Natural Vector Method.

Zika virus (ZIKV) is a mosquito-borne flavivirus. It was first isolated from Uganda in 1947 and has become an emergent event since 2007. However, because of the inconsistency of alignment methods, the evolution of ZIKV remains poorly understood. In this study, we first use the complete protein and an alignment-free method to build a phylogenetic tree of 87 Zika strains in which Asian, East African, and West African lineages are characterized. We also use the NS5 protein to construct the genetic relationship among 44 Zika strains. For the first time, these strains are divided into two clades: African 1 and African 2. This result suggests that ZIKV originates from Africa, then spread to Asia, Pacific islands, and throughout the Americas. We also perform the phylogeny analysis for 53 viruses in genus Flavivirus to which ZIKV belongs using complete proteins. Our conclusion is consistent with the classification by the hosts and transmission vectors.