Abstract
Ebolavirus (EBOV) life cycle involves interactions with numerous host factors, but it remains poorly understood, as does pathogenesis. Herein, we synthesized 65 siRNAs targeting host genes mostly connected with aspects of the negative-sense RNA virus life cycle (including viral entry, uncoating, fusion, replication, assembly, and budding). We produced EBOV transcription- and replication-competent virus-like particles (trVLPs) to mimic the EBOV life cycle. After screening host factors associated with the trVLP life cycle, we assessed interactions of host proteins with trVLP glycoprotein (GP), VP40, and RNA by co-immunoprecipitation (Co-IP) and chromatin immunoprecipitation (ChIP). The results demonstrate that RNAi silencing with 11 siRNAs (ANXA5, ARFGAP1, FLT4, GRP78, HSPA1A, HSP90AB1, HSPA8, MAPK11, MEK2, NTRK1, and YWHAZ) decreased the replication efficiency of trVLPs. Co-IP revealed nine candidate host proteins (FLT4, GRP78, HSPA1A, HSP90AB1, HSPA8, MAPK11, MEK2, NTRK1, and YWHAZ) potentially interacting with trVLP GP, and four (ANXA5, GRP78, HSPA1A, and HSP90AB1) potentially interacting with trVLP VP40. Ch-IP identified nine candidate host proteins (ANXA5, ARFGAP1, FLT4, GRP78, HSPA1A, HSP90AB1, MAPK11, MEK2, and NTRK1) interacting with trVLP RNA. This study was based on trVLP and could not replace live ebolavirus entirely; in particular, the interaction between trVLP RNA and host proteins cannot be assumed to be identical in live virus. However, the results provide valuable information for further studies and deepen our understanding of essential host factors involved in the EBOV life cycle.
Introduction
Ebolavirus (EBOV) is a single-stranded, negative-sense RNA virus with a heterogeneous filamentous structure (Jun et al., 2015) that was first reported in 1976 as the cause of Ebola viral disease (EVD) in humans and other primates (Johnson et al., 1977). Five different EBOVs have been defined: Ebola virus (EBOV, previously known as Zaire ebolavirus), Sudan virus (SUDV), Bundibugyo virus (BDBV), Taï Forest virus (TAFV) and Reston virus (RESTV) (Kuhn, 2017). The EVD epidemic from December 2013 to March 2016 in Western Africa was caused by EBOV, and was associated with >28,000 cases and >11,000 deaths in 11 countries (Ajisegiri et al., 2018).
Although numerous studies have been devoted to EVD therapeutics, such as immune-based treatments and small molecule inhibitors, and considerable progress has been made (Fabozzi et al., 2011; Oestereich et al., 2014; Wolf et al., 2015; Sissoko et al., 2016), the molecular basis of the EBOV life cycle and its relevance to pathogenesis remain poorly understood. For example, the mechanisms, proteins and receptors mediating viral entry, the host enzymes triggering and accelerating uncoating and virus fusion, and the host factors facilitating VP40 oligomerisation and association of VP40 with the inner leaflet remain largely unknown (Johnson et al., 2016; Gc et al., 2017; Younan et al., 2017; Kurosaki et al., 2018).
To model the EBOV life cycle using a safe method, a series of trVLPs have been developed (Hoenen et al., 2014). These trVLPs could express the EBOV proteins required for genome replication and transcription to model EBOV life cycles under biosafety level 2 conditions. Based on the EBOV life cycle in host cells, and on previous RNAi screening that identified filovirion-associated and secretory pathway-related host factors (Spurgers et al., 2010; Simpson et al., 2012; Sakurai et al., 2015), we herein synthesized siRNAs targeting host genes connected to membrane traffic machinery, endoplasmic reticulum-Golgi recycling, inhibition of secretion, and lipid droplet formation, and assessed their role in the EBOV life cycle. These genes are mostly connected with virus life cycle, specifically viral entry, uncoating, fusion, replication, assembly and budding.
In the EBOV life cycle, GP is critical for EBOV entry, uncoating and fusion processes (Nanbo et al., 2010; Takada, 2012), while VP40 also plays a crucial role in EBOV assembly and budding (Dessen et al., 2000; Balmith and Soliman, 2017), and both proteins engage in important interactions with host cellular and plasma membranes. Thus, host factors that exert a significant influence on the EBOV life cycle were identified by RNAi silencing, and interactions between EBOV GP, VP40, and nucleic acids were evaluated by Co-IP or ChIP. This is the first systematic survey of host factors that affect the EBOV life cycle, demonstrating interactions between host proteins and EBOV GP, VP40 and RNA.
Materials and Methods
Cell Line and Plasmids
Human embryonic kidney (HEK) 293T cells were cultured in DMEM (Thermo Fisher, Waltham, MA, United States; Cat# 10566016) containing 10% FBS (Gibco, Waltham, MA, United States; Cat# 10099141), 2 mM L-glutamine (Life Technologies, Waltham, MA, United States; Cat# 25030081), and 1% penicillin- streptomycin (Life Technologies; Cat# 10378016) at 37°C with 5% CO2. Plasmids pCAGGS-VP30, pCAGGS-VP35, pCAGGS-NP, pCAGGS-L, p4cis-vRNA-RLuc, pCAGGS-Tim1, and pCAGGS-T7 were kindly provided by Drs. Heinz Feldmann and Thomas Hoenen, Rocky Mountain Laboratories, National Institute of Health (NIH). Plasmid p4cis-vRNA-RLuc containing EBOV non-coding regions, a reporter gene, and three genes (VP40, GP1, 2 and VP24) involved in morphogenesis, budding and cell entry, was used to produce trVLPs. In the following experiments, all tests were carried out with EBOV trVLPs.
Production of trVLPs
Producer 293T cells (p0) in 6-well plates were transfected with plasmids encoding each EBOV structural protein (75 ng pCAGGS-VP30, 125 ng pCAGGS-VP35, 125 ng pCAGGS-NP, 1,000 ng pCAGGS-L, 250 ng p4cis-vRNA-RLuc), and pCAGGS-T7 (250 ng) encoding a Renilla luciferase reporter (Hoenen et al., 2014). Cell supernatants (200 ml) from 10 6-well plates containing released trVLPs were harvested at 72 h post-transfection, and cells and cellular debris were pelleted by gentle centrifugation at 175 × g. The supernatant was then used to infect target 293T cells (p1) previously transfected with RNP components (125 ng pCAGGS-NP, 125 ng pCAGGS-VP35, 75 ng pCAGGS-VP30, 1,000 ng pCAGGS-L, 250 ng pCAGGS-Tim1) (Hoenen et al., 2014). Target 293T cells (p2–p5) were treated in a similar manner. Cleared supernatants (33 ml) in 2 ml of 20% sucrose from the bottom of each tube were centrifuged at 125,000 × g in a SW-28 rotor for 3 h at 4°C. The resulting pellets were resuspended in 100 μl ice-cold NTE buffer (10 mM Tris pH 7.5, 100 mM NaCl, 1 mM EDTA) by tapping gently about 100 times, and trVLPs were stored on ice or in a refrigerator at 4°C until use.
Determining the 50% Tissue Culture Infective Dose (TCID50) of trVLPs
To determine the TCID50 values of trVLPs, CPEs were observed using microscopy, and the results from each dilution were used to calculate TCID50 values using the Reed–Muench method (Ramakrishnan, 2016). After harvest, trVLPs were 10-fold serially diluted with DMEM at concentrations ranging from 10-1 to 10-10. Attenuated trVLPs (100 μl) were added to eight wells in each row of a 96-well plate, and 293T cell suspension (100 μl) was added to each well to a final cell density of 2 × 105 cells/ml. Additionally, 293T cells not infected with trVLPs were included as controls. CPEs were observed and recorded each day for 7 days, and TCID50 values were calculated according the Reed–Muench method, giving a value of 103.76/0.1 ml (detailed data are listed in Supplementary Table S1).
Electron Microscopy (EM) Analysis of trVLPs
Purified trVLPs were pipetted onto a 300-mesh copper grid coated with carbon film, incubated for 5 min at room temperature, and grids were washed twice with distilled water and negatively stained for 15 s with 1% uranyl acetate. Excess liquid was removed with a filter paper and trVLPs were visualized under a Hitachi H7000 transmission electron microscope.
Imaging of DiI-Labeled trVLP Internalization in Live 293T Cells
Purified trVLPs (∼100 μg viral protein) were resuspended in NTE buffer a final volume of 1 ml and 1,1′-dioctadecyl-3,3,3′,3′-tetramethylindocarbocyanine perchlorate (DiI, dissolved in ethanol to a final concentration of 10 μM; Invitrogen, Waltham, MA, United States; Cat# D282) was added and mixed thoroughly. The mixture was shaken gently for 1 h at room temperature then passed through a 0.22 μm filter (Millipore, Darmstadt, Germany; Cat# SLGP033RB). HEK 293T cells cultured in 8-well chamber slides with removable wells were mixed with 4′,6-diamidino-2-phenylindole (DAPI; Invitrogen; Cat# P36931) at 10 μg/ml for 10 min to stain nuclei, and 293T cells were then incubated with DiI-labeled trVLPs at 4°C for 10 min and washed with PBS. The internalization of DiI-labeled trVLPs in 293T cells was imaged at different times from 10 to 60 min at 37°C on the stage of a fluorescence microscope.
Selection of Candidate Genes
To select host factors that are potentially related to various stages of the EBOV life cycle (including viral entry, uncoating, fusion, replication, assembly and budding), we chose host genes that are associated with membrane traffic machinery, endoplasmic reticulum-Golgi recycling, inhibition of secretion, and lipid droplet formation. For example, the knockdown of FLT4, RAB11B, and PLCH1 resulted in high-level secretion inhibition; GOSR1, LTK, and CNKSR1 are linked to endoplasmic reticulum-Golgi recycling; VAMP2, RAB11B, and COPG are membrane traffic machinery-related molecules (Simpson et al., 2012); HSPA8, HSP90AB1, and ANXA5 are believed to be related to filamentous virus (Spurgers et al., 2010). A total of 65 genes probably associated with the EBOV life cycle or virus life cycle-related functions (Spurgers et al., 2010; Simpson et al., 2012; Sakurai et al., 2015) were selected according to preliminary screening. Details and references for these genes are listed in Supplementary Table S3.
SiRNA Screening and trVLP Infection
SiRNAs targeting the 65 selected genes were designed and three sequences were synthesized for each gene. To choose the most efficient sequence for RNA interfere (RNAi), the gene knockdown efficiency of each sequence was evaluated by quantitative PCR, with the 18sRNA gene as a control. Values were calculated using the 2-ΔΔCt method, and the smaller the value, the higher the knockdown efficiency. The most efficient siRNAs were chosen for RNAi screening tests (data listed in Supplementary Table S2), and details of each chosen siRNA are included in Supplementary Table S3. In 24-well plates, 100 μl of opti-MEM medium (Invitrogen; Cat# 31985070) containing 1.4 μl siRNA and 4.5 μl HiPerFect transfection reagent (Qiagen, Dusseldorf, Germany; Cat# 301705) was added to each well and plates were shaken gently for 1 min. After 10 min incubation at room temperature, a cell suspension (400 μl) containing 1 × 105 cells was added to give a final siRNA concentration of 75 nM. Cells were incubated at 37°C and 5% CO2 for 48 h, washed with PBS, and infected with trVLPs (TCID50 = 102.76/0.01 ml) for 48 h. Negative control siRNA (Qiagen; Cat# 1027310) served as a control.
RNA Extraction and Quantitative Real-Time Reverse Transcription PCR (qRT-PCR)
After trVLP infection, 293T cells in 24-well plates were washed with PBS, and total RNA was extracted from cells using TRIzol (Invitrogen; Cat# 15596018) according to the manufacturer’s instructions. qRT-PCR was performed with a EBOV nucleic acid test kit (Zhijiang Bio-tech, Shanghai, China; Cat# QR-0220-02) on an ABI 7500 qPCR system (45°C for 10 min and 95°C for 15 min, followed by 45 cycles at 95°C for 15 s and 60°C for 30 s). The positive control, provided in the kit, was serially diluted over three orders of magnitude, and the absolute quantity of EBOV RNA (copy/ml) was calculated using the comparative Ct method and a standard curve. All qRT-PCR experiments were performed in triplicate and repeated independently three times.
Antibodies
Antibodies recognizing HSP70 (HSPA1A; Cat# 4873S), HSP90β (HSP90AB1; Cat# 5087S), MEK1/2 (Cat# 9122S), ARFGAP1 (Cat# 14522S) and Trk (NTRK1; Cat# 92991S) were purchased from Cell Signaling Technology (Danvers, MA, United States). Antibodies recognizing MAPK11 (Cat# ab183208), VEGF Receptor 3 (FLT4; Cat# ab10284), Annexin V (ANXA5; Cat# ab54775), GRP78 (HSPA5; Cat# ab21685), GAPDH (Cat# ab8245), and Ebola virus GP (Cat# ab1918) were purchased from Abcam (Cambridge, United Kingdom). Antibodies recognizing 14-3-3 zeta/beta (YWHAZ; Cat# NB100-1964) and Hsc70 (HSPA8; Cat# NBP2-12880) were purchased from Novus Biologicals (Littleton, CO, United States). Antibody recognizing Ebola virus VP40 (Cat# sc-51872) was purchased from Santa Cruz Biotechnology (Dallas, TX, United States). Antibodies against HSPA1A, HSP90AB1, MEK1/2, ARFGAP1, NTRK1, MAPK11, FLT4, and GRP78 were rabbit-derived IgGs. Antibodies against ANXA5, YWHAZ, and HSPA8 were mouse-derived IgGs. Normal mouse IgG (Cat# SC-2025) and normal rabbit IgG (Cat# SC-2027) were obtained from Santa Cruz Biotechnology and used as isotype controls in immunoprecipitation tests.
Co-IP and Immunoblot Analysis
HEK 293T cells in 12-well plates were infected with trVLPs (TCID50 = 103.76/0.1 ml) for 48 h and whole-cell extracts were prepared followed by incubation overnight with appropriate antibodies recognizing host proteins plus Protein G beads (GE Healthcare, Uppsala, Sweden; Cat# 28-9440-08). Antibodies mixed with Protein G beads and normal cell/extracts without trVLPs, and Protein G beads mixed with trVLP-treated cell extracts without antibodies, were included as control groups. Beads were then washed five times with low-salt lysis buffer, and immunoprecipitates were eluted with 2× SDS Loading Buffer and resolved by SDS-PAGE. Proteins were transferred to PVDF membranes and further incubated with anti-GP or VP40 antibodies. Proteins of interest were detected with the Super Signal West Pico chemiluminescent substrate (Thermo Fisher; Cat# 34580). To confirm that siRNAs interfered with the expression of target proteins in a specific manner, immunoblot analysis of target host proteins was performed after siRNA transfection. Meanwhile, to confirm the specificity of the antibodies, target host proteins, GP, and VP40 proteins from the above extracts were tested before and after Co-IP separately (normal mouse or rabbit IgG served as isotype controls). The housekeeping protein GAPDH was used as an internal control for immunoblot analyses.
Chromatin Immunoprecipitation (ChIP)
After infection of trVLPs (TCID50 = 102.76/0.01 ml) for 48 h, HEK 293T cells in 12-well plates were fixed with 1% formaldehyde for 10 min, and glycine (Thermo Fisher; Cat# 28363) was added to terminate the cross-linking reaction. Whole-cell extracts were then prepared using a Simple ChIP Enzymatic Chromatin IP Kit (Cell Signaling Technology; Cat# 9003) following the manufacturer’s instructions with appropriate antibodies recognizing host proteins, and normal mouse or rabbit IgGs served as isotype controls. Beads were then washed and cross-linking substances were unlocked according the ChIP kit instructions. Purified nucleic acids were then amplified by qPCR with a qPCR Probe Kit (Vazyme Biotech, Nanjing, China; Cat# Q2223-01). Primers and probe were designed according to the sequence of Zaire ebolavirus isolate Homo sapiens-wt/GIN/2014/Gueckedou-C05 (GenBank No. KJ660348.2) as follows: EN-Z-F, ATGATGGARGCTACGGCG; EN-Z-R, TGCCCTCTGTATGCTGGCCCT; EN-Z-P, CCARAGTTACTCGGAAAACGGCATG.
Sequencing Analysis
The PCR products amplified from ChIP samples were purified by agarose gel electrophoresis, recycled, and sequenced by the Sanger chain termination method. Genotypes were attributed using BioEdit software (Ibis Biosciences, Carlsbad, CA, United States) for allelic discrimination.
Statistical Analysis
Data in normal distributions are presented as the mean ± standard deviation (Kim et al., 2018). Student’s t-tests (two-tailed, paired or unpaired t-tests with Welch’s correction), and multiple comparisons were performed where appropriate. All analysis was performed with Graphpad Prism 5 (GraphPad Software, United Kingdom) and differences between group means with p < 0.05 were considered statistically significant.
Results
EBOV-trVLPs Successfully Simulate Authentic Ebola Virions
EBOV-trVLPs synthesized as described above displayed a filamentous-like morphology with a diameter of 100–500 μm as visualized by EM (Figure 1A). Meanwhile, internalization and tracking of DiI-labeled trVLPs in live 293T cells visualized by fluorescent microscopy indicated that trVLPs gradually translocated into the cytoplasm over time (Figure 1B). It has been reported that filamentous viruses, including authentic EBOV, are internalized into cells following stimulation of the macropinocytosis pathway, while spherical pseudotype virions are internalized via the endocytosis pathway (Nanbo et al., 2010). In our present study, trVLPs are filamentous virions, hence we presumed that trVLPs would be internalized via the macropinocytosis pathway.
FIGURE 1
SiRNA Knockdown of Host Factors Decreases Replication of trVLPs
To identify factors associated with the life cycle of trVLPs, siRNAs targeting 65 selected genes were tested in 293T cells infected with trVLPs in triplicate by qRT-PCR. The results showed that siRNA knockdown of ANXA5, ARFGAP1, FLT4, GRP78, HSPA1A, HSP90AB1, HSPA8, MAPK11, MEK2, NTRK1, and YWHAZ efficiently inhibited trVLP replication (p < 0.05; Figure 2A). To confirm that these siRNAs interfered with the expression of target proteins in a specific manner, we measured the levels of target proteins by western blot analysis, and this demonstrated significant interference for each target protein as expected (Figure 2B). The results suggest that these host proteins are highly likely to participate in the trVLP life cycle, and may perform critical roles. In addition, the results from the siRNA knockdowns of other host proteins suggest that they do not interfere with trVLP replication (Supplementary Figure S1).
FIGURE 2
Among these selected candidate host factors, four are chaperones (HSPA1A, HSP90AB1, HSPA8, and GRP78), four are involved in plasma and/or endosome membrane trafficking (ANXA5, ARFGAP1, FLT4, and NTRK1), and three play important roles in signal transduction pathways (MAPK11, MEK2, and YWHAZ). Further bioinformatic analysis by GeneMANIA revealed that these candidate host factors were markedly enriched in gene categories associated with 14-3-3 protein families [such as YWHAZ, YWHAE and YWHAQ)], annexin families (including ANXA1, ANXA4, and ANXA5), heat shock protein (HSP) families (e.g., HSPA1A, HSP90AB1, GRP78, and HSPA8), and MAP kinase families (including MAPK11 and MAP2K2; Figure 2C). The analysis also indicated that candidate host factors are associated with cytoplasmic vesicle membranes, mitochondrial membranes, mitochondria, the COP9 signalosome, and phospholipid binding, and strongly connected with antigen binding, viral endocytosis, transportation, anti-virus responses, and viral budding.
Since GP and VP40 proteins perform critical functions in the EBOV life cycle and since both interact with host cell membranes, we performed Co-IP and Ch-IP assays to validate host protein-GP/VP40 and host protein-trVLP nucleic acid interactions.
Candidate Host Proteins Interact With trVLP GP
Antibodies recognizing the 11 host proteins were incubated overnight with Protein G beads and cell extracts infected with trVLPs described above. Immunoblot analysis was then carried out with anti-GP antibody, and the results demonstrated that nine host proteins (FLT4, GRP78, HSPA1A, HSP90AB1, HSPA8, MAPK11, MEK2, NTRK1, and YWHAZ) interacted with GP, whereas ARFGAP1 and ANXA5 did not appear to interact with GP (normal mouse or rabbit IgG served as isotype controls; Figure 3A). Meanwhile, immunoblot analyses of the 11 target host proteins before and after IP indicated a high degree of specificity for those antibodies, further verifying the interactions between GP and candidate host proteins (Figures 3B,C).
FIGURE 3
Candidate Host Proteins Interact With trVLP VP40
Next, interactions of the 11 candidate host proteins with VP40 were tested by incubating antibodies recognizing the 11 host proteins overnight with Protein G beads and trVLPs-treated cell extracts. Western blotting with anti-VP40 antibody indicated that ANXA5, GRP78, HSPA1A, and HSP90AB1 interacted with VP40, whereas ARFGAP1, HSPA8, MAPK11, MEK2, FLT4, NTRK1, and YWHAZ did not interact with VP40 (normal mouse or rabbit IgG served as isotype controls; Figure 4A). Immunoblot analyses of the 11 candidate host proteins before and after IP confirmed the specificity of those antibodies, as well as interactions between VP40 and candidate host proteins (Figures 4B,C).
FIGURE 4
Candidate Host Proteins Interact With trVLP RNA
Since interactions between virus RNA and host proteins are likely to be essential for viral replication, we next investigated whether trVLP RNA interacts with the candidate host proteins using Ch-IP assays and purified nucleic acids by qPCR and DNA sequencing. BLAST analysis of the sequencing results indicated that sequences from nine antibodies recognizing candidate host proteins (HSP90AB1, ANXA5, ARFGAP1, FLT4, GRP78, HSPA1A, MEK2, NTRK1, and MAPK11) were consistent with Zaire ebolavirus (GenBank No. KJ660348.2; Figure 5). HSPA8, YWHAZ, and isotype controls did not recognize trVLP RNA.
FIGURE 5
Discussion
A better understanding of the mechanisms of the EBOV life cycle and host–virus interactions (including protein–protein and protein–RNA interactions) would greatly assist the elucidation of EBOV pathogenesis and therapeutic development. Herein, we identified 11 candidate host proteins by RNAi silencing that appear to function in the EBOV-trVLP life cycle among 65 host factors presumed to be associated with the virus life cycle and secretary pathways. Meanwhile, we measured candidate host protein levels after RNAi by western blot analysis, and the results revealed significant interference effects and confirmed the specificity of RNAi silencing. Bioinformatic analysis by GeneMANIA indicated that these candidate host proteins belong to the following families: (1) HSPs that participate in a wide variety of cellular processes, including activation of proteolysis, protecting the proteome against stress, infection, and transport of newly synthesized polypeptides (Alvarez-Sanchez et al., 2015; Chaudhary et al., 2016; Heck et al., 2017); (2) 14-3-3 protein families that bind a number of signaling proteins including kinases, phosphatases and transmembrane receptors (Radhakrishnan et al., 2012; Obsilova et al., 2014; Aghazadeh and Papadopoulos, 2016); (3) Annexin families that are involved in trafficking and organization of endocytosis, vesicles, exocytosis, and calcium ion channel formation (Grieve et al., 2012; Kuehnl et al., 2016); (4) MAP kinase families, which play critical roles in signal transduction, phosphorylation and activation of MAPK1/ERK2 and MAPK3/ERK1 pathways (Vomastek et al., 2008; Robinson and Pitcher, 2013; Cook et al., 2017). Further functional analysis confirmed that the candidate host proteins were strongly associated with cytoplasmic vesicle membranes, mitochondrial membranes, mitochondria, the COP9 signalosome, and phospholipid binding, and connected with antigen binding, viral endocytosis, transportation, anti-virus responses, and viral budding (Bech-Otschir et al., 2001; Horner et al., 2011; van Zuylen et al., 2012; Richard et al., 2015).
Heat shock protein family proteins are known to be associated with many viruses. For example, the Hsp70 chaperone network mediates the dengue virus life cycle, GRP78 is a potential defensive molecule against hepatitis A virus replication, and Hsp90 controls reactivation of human immunodeficiency virus type 1 from latency (Anderson et al., 2014; Taguwa et al., 2015; Vashist et al., 2015; Jiang et al., 2017). Meanwhile, filovirion studies showed that HSP70 family proteins are associated with Ebola virus (Nelson et al., 2017), consistent with our current results. Members of 14-3-3 protein families are also connected with viral infection, since they mediate signaling pathways and interactions during viral dsRNA stimulation and advanced antiviral innate immune responses (Ohman et al., 2010; Ohman et al., 2014). However, interactions of 14-3-3 proteins with EBOV have not been reported previously. Similarly, annexin family proteins also interact closely with RNA viruses such as measles virus and H5N1 avian influenza virus (Ma et al., 2017; Koga et al., 2018), but interactions with EBOV have not been reported previously. MAP kinases are essential in signal transduction and involved in numerous virus life cycle events, including in EBOV (Strong et al., 2008; Johnson et al., 2014), but the detailed mechanisms of interactions between MAP kinases and EBOV remain ambiguous.
To confirm the interactions of candidate host proteins with EBOV, Co-IP, and Ch-IP were performed. The results revealed that HSPA1A, HSP90AB1, HSPA8, MAPK11, NTRK1, MEK2, FLT4, GRP78, and YWHAZ interacted with EBOV GP. Meanwhile, immunoblot analyses of the 11 target host proteins before and after IP confirmed the high degree of specificity for those antibodies, and specific interactions between GP and candidate host proteins. Four of these proteins belong to HSP families, which are involved in a variety of cellular functions, and some of the others function in plasma and/or endosome membrane trafficking and signal transduction pathways. These results are consistent with the functions of EBOV GP, which is critical for EBOV cell entry, and uncoating and fusion processes occurring at the cell membrane through the endosomal pathway. Protein interactions can involve various mechanisms including subunit polymerization, cross polymerization, and molecular recognition. However, the detailed mechanisms of the interactions between GP and host proteins remain unclear. Meanwhile, it is worth noting that additional factors likely act as a bridge between the candidate host proteins and GP.
Co-immunoprecipitation and immunoblot analyses also revealed that ANXA5, GRP78, HSPA1A, and HSP90AB1 interact with VP40, and bioinformatic analysis indicated that ANXA5 engages in phospholipid binding, while GRP78 and HSP90AB1 are associated with the COP9 signalosome, which is believed to act as a scaffold to facilitate spatial sequestration for multiple interacting molecules that regulate protein stability, transcription, protein phosphorylation and intracellular distribution (Lee et al., 2013; Mosadeghi et al., 2016). The functions of these candidate host proteins are consistent with the role of VP40 in the EBOV life cycle, which presides over association with plasma membrane phosphatidylserine, the inner leaflet, and translocation of virions through the plasma membrane (Dessen et al., 2000; Balmith and Soliman, 2017). Regardless of whether the modes of interaction between the candidate host proteins and VP40 are direct, indirect, or complex, these proteins appear to be linked to the EBOV life cycle.
Interactions of host proteins with trVLP RNA were assessed, and the results indicated that HSP90AB1, HSPA1A, GRP78, ANXA5, ARFGAP1, FLT4, MAPK11, MEK2, and NTRK1 interact with trVLP RNA. Protein-RNA interactions are important for various cellular process including transcriptional and post-transcriptional regulation, and activating immune responses against virus infection (Bidet and Garcia-Blanco, 2014; Cheng et al., 2015; Borodavka et al., 2017; Cipriano and Ballarino, 2018). In the EBOV life cycle, viral nucleic acids are exposed after uncoating in the late endosome, and negative-sense RNA then acts as a template for the synthesis of positive-sense RNA, which guides transcription of viral mRNA. Host protein-virus RNA interactions are presumably crucial during the EBOV life cycle, as indicated by our results.
There were some limitations to our study. Firstly, all the experiments were based on EBOV trVLPs and HEK 293T cells, the interactions may not reflect the physiology of live EBOV. Further experiments using EBOV are therefore required to confirm these current findings. Secondly, the interaction patterns are unclear, which may work either through direct binding or through the indirect action of one or more intermediate molecules, or through the formation of protein complexes by bridge proteins. Furthermore, it is unclear which region or segment of the GP or VP40 protein interacts with the host proteins. Finally, as ebolavirus probably infects different cell types via different mechanisms, the interactions we found in 293T cells may not be relevant in other host cells, such as monocytes and hepatocytes (Kindrachuk et al., 2014; Rogers and Maury, 2018). Subsequent experiments will attempt to define the patterns of interactions, clarify the specificity and binding sites.
Conclusion
Analysis of EBOV life cycle-associated host factors by RNAi silencing identified four chaperones (HSPA1A, HSP90AB1, HSPA8, and GRP78), four plasma and endosome membrane trafficking molecules (ANXA5, ARFGAP1, FLT4, and NTRK1) and three signal transduction elements (MAPK11, MEK2, and YWHAZ) that may be essential for EBOV replication. Among these host proteins, nine (FLT4, GRP78, HSPA1A, HSP90AB1, HSPA8, MAPK11, MEK2, NTRK1, and YWHAZ) potentially interact with trVLP GP, four (ANXA5, GRP78, HSPA1A, and HSP90AB1) potentially interact with trVLP VP40, and nine (ANXA5, ARFGAP1, FLT4, GRP78, HSPA1A, HSP90AB1, MAPK11, MEK2, and NTRK1) participate in host protein-trVLP RNA interactions. Although EBOV trVLP cannot replace live ebolavirus entirely, and these interactions cannot be assumed to be identical in live virus, the study provides valuable foundations for further research.
Statements
Author contributions
D-SY, T-HW, C-YH, Z-GW, Y-HL, and L-FC performed the experiments. H-PY and D-SY performed the statistical analysis. H-PY, N-PW, and L-JL designed the study and drafted the manuscript. All authors participated in writing the manuscript.
Funding
This work was supported by the Major Program of National Natural Science Foundation of China [Grant Number 81590763].
Acknowledgments
We greatly appreciate the Elixigen Corporation which is a professional scientific editing company for proofreading the manuscript.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2018.02724/full#supplementary-material
FIGURE S1RNAi silencing analysis of host factors that failed to inhibit ebolavirus transcription- and replication-competent virus-like particle (EBOV-trVLP) replication. Samples (100 μl) of opti-MEM medium containing 1.4 μl siRNA and 4.5 μl HiPerFect were placed in 24-well plates, and a cell suspension (400 μl) containing 1 × 105 293T cells was added to give a final siRNA concentration of 75 nM. After incubation for 48 h, cells were infected with trVLPs for another 48 h, then total RNA was extracted, and the absolute quantity of EBOV RNA was measured using an EBOV nucleic acid test kit. All qRT-PCR experiments were performed in triplicate and repeated three times independently. Cells not transfected with siRNA but infected with trVLPs served as a blank control, and cells transfected with isotype siRNA and infected with trVLPs served as a negative control.
TABLE S1TCID50 of trVLPs calculated by Reed–Muench method.
TABLE S2The gene knock-down efficiency of the most efficient siRNAs chosen for RNAi screening tests.
TABLE S3Information of each gene chosen for RNA interference screening tests.
Abbreviations
- ChIP
chromatin immunoprecipitation
- Co-IP
co-immunoprecipitation
- CPE
cytopathic effect
- DAPI
4′,6-diamidino-2-phenylindole
- DiI
1,1′-Dioctadecyl-3,3,3′,3′-Tetramethylindocarbocyanine Perchlorate
- DMEM
Dulbecco’s modified Eagle’s medium
- EM
electron microscopy
- EVD
Ebola virus disease
- FBS
fetal calf serum
- GP
glycoprotein
- HEK
human embryonic kidney
- PAGE
polyacrylamide gel electrophoresis
- PBS
phosphate-buffered saline
- PVDF
polyvinylidene fluoride
- qRT-PCR
quantitative real-time reverse transcription polymerase chain reaction
- RNAi
RNA interference
- RNP
ribonucleoprotein
- SDS
sodium dodecyl sulfate
- TCID50
50% tissue culture infective dose
- trVLPs
transcription- and replication-competent virus-like particles
- VP40
matrix protein 40
- Zaire EBOV
Zaire Ebola virus
References
1
AghazadehY.PapadopoulosV. (2016). The role of the 14-3-3 protein family in health, disease, and drug development.Drug Discov. Today21278–287. 10.1016/j.drudis.2015.09.012
2
AjisegiriW. S.ChughtaiA. A.MacIntyreC. R. (2018). A risk analysis approach to prioritizing epidemics: ebola virus disease in West Africa as a case study.Risk Anal.38429–441. 10.1111/risa.12876
3
Alvarez-SanchezC.Cancio-LonchesC.Mora-HerediaJ. E.Santos-ValenciaJ. C.Barrera-VazquezO. S.Yocupicio-MonroyM.et al (2015). Negative effect of heat shock on feline calicivirus release from infected cells is associated with the control of apoptosis.Virus Res.19844–52. 10.1016/j.virusres.2015.01.003
4
AndersonI.LowJ. S.WestonS.WeinbergerM.ZhyvoloupA.LabokhaA. A.et al (2014). Heat shock protein 90 controls HIV-1 reactivation from latency.Proc. Natl. Acad. Sci. U.S.A.111E1528–E1537. 10.1073/pnas.1320178111
5
BalmithM.SolimanM. E. (2017). VP40 of the ebola virus as a target for EboV therapy: comprehensive conformational and inhibitor binding landscape from accelerated molecular dynamics.Cell Biochem. Biophys.7565–78. 10.1007/s12013-017-0783-8
6
Bech-OtschirD.KraftR.HuangX.HenkleinP.KapelariB.PollmannC.et al (2001). COP9 signalosome-specific phosphorylation targets p53 to degradation by the ubiquitin system.EMBO J.201630–1639. 10.1093/emboj/20.7.1630
7
BidetK.Garcia-BlancoM. A. (2014). Flaviviral RNAs: weapons and targets in the war between virus and host.Biochem. J.462215–230. 10.1042/bj20140456
8
BorodavkaA.DykemanE. C.SchrimpfW.LambD. C. (2017). Protein-mediated RNA folding governs sequence-specific interactions between rotavirus genome segments.eLife6:e27453. 10.7554/eLife.27453
9
ChaudharyP.KhanS. Z.RawatP.AugustineT.RaynesD. A.GuerrieroV.et al (2016). HSP70 binding protein 1 (HspBP1) suppresses HIV-1 replication by inhibiting NF-kappaB mediated activation of viral gene expression.Nucleic Acids Res.441613–1629. 10.1093/nar/gkv1151
10
ChengZ.ZhouS.GuanJ. (2015). Computationally predicting protein-RNA interactions using only positive and unlabeled examples.J. Bioinform. Comput. Biol.13:1541005. 10.1142/s021972001541005x
11
CiprianoA.BallarinoM. (2018). The ever-evolving concept of the gene: the use of RNA/Protein experimental techniques to understand genome functions.Front. Mol. Biosci.5:20. 10.3389/fmolb.2018.00020
12
CookS. J.StuartK.GilleyR.SaleM. J. (2017). Control of cell death and mitochondrial fission by ERK1/2 MAP kinase signalling.FEBS J.2844177–4195. 10.1111/febs.14122
13
DessenA.VolchkovV.DolnikO.KlenkH. D.WeissenhornW. (2000). Crystal structure of the matrix protein VP40 from Ebola virus.EMBO J.194228–4236. 10.1093/emboj/19.16.4228
14
FabozziG.NabelC. S.DolanM. A.SullivanN. J. (2011). Ebolavirus proteins suppress the effects of small interfering RNA by direct interaction with the mammalian RNA interference pathway.J. Virol.852512–2523. 10.1128/jvi.01160-10
15
GcJ. B.GerstmanB. S.ChapagainP. P. (2017). Membrane association and localization dynamics of the Ebola virus matrix protein VP40.Biochim. Biophys. Acta18592012–2020. 10.1016/j.bbamem.2017.07.007
16
GrieveA. G.MossS. E.HayesM. J. (2012). Annexin A2 at the interface of actin and membrane dynamics: a focus on its roles in endocytosis and cell polarization.Int. J. Cell Biol.2012:852430. 10.1155/2012/852430
17
HeckT. G.ScomazzonS. P.NunesP. R.ScholerC. M.da SilvaG. S.BittencourtA.et al (2017). Acute exercise boosts cell proliferation and the heat shock response in lymphocytes: correlation with cytokine production and extracellular-to-intracellular HSP70 ratio.Cell Stress Chaperones22271–291. 10.1007/s12192-017-0771-3
18
HoenenT.WattA.MoraA.FeldmannH. (2014). Modeling the lifecycle of Ebola virus under biosafety level 2 conditions with virus-like particles containing tetracistronic minigenomes.J. Vis. Exp.91:52381. 10.3791/52381
19
HornerS. M.LiuH. M.ParkH. S.BrileyJ.GaleM.Jr. (2011). Mitochondrial-associated endoplasmic reticulum membranes (MAM) form innate immune synapses and are targeted by hepatitis C virus.Proc. Natl. Acad. Sci. U.S.A.10814590–14595. 10.1073/pnas.1110133108
20
JiangX.KandaT.HagaY.SasakiR.NakamuraM.WuS.et al (2017). Glucose-regulated protein 78 is an antiviral against hepatitis A virus replication.Exp. Ther. Med.133305–3308. 10.3892/etm.2017.4407
21
JohnsonJ. C.MartinezO.HonkoA. N.HensleyL. E.OlingerG. G.BaslerC. F. (2014). Pyridinyl imidazole inhibitors of p38 MAP kinase impair viral entry and reduce cytokine induction by Zaire ebolavirus in human dendritic cells.Antiviral Res.107102–109. 10.1016/j.antiviral.2014.04.014
22
JohnsonK. A.TaghonG. J.ScottJ. L.StahelinR. V. (2016). The Ebola Virus matrix protein, VP40, requires phosphatidylinositol 4,5-bisphosphate (PI(4,5)P2) for extensive oligomerization at the plasma membrane and viral egress.Sci. Rep.6:19125. 10.1038/srep19125
23
JohnsonK. M.LangeJ. V.WebbP. A.MurphyF. A. (1977). Isolation and partial characterisation of a new virus causing acute haemorrhagic fever in Zaire.Lancet1569–571.
24
JunS. R.LeuzeM. R.NookaewI.UberbacherE. C.LandM.ZhangQ.et al (2015). Ebolavirus comparative genomics.FEMS Microbiol. Rev.39764–778. 10.1093/femsre/fuv031
25
KimY. H.LeeJ.KimY. E.ChongC. K.PinchemelY.ReisdorferF.et al (2018). Development of a rapid diagnostic test kit to detect IgG/IgM antibody against zika virus using monoclonal antibodies to the envelope and non-structural protein 1 of the virus.Korean J. Parasitol.5661–70. 10.3347/kjp.2018.56.1.61
26
KindrachukJ.Wahl-JensenV.SafronetzD.TrostB.HoenenT.ArsenaultR.et al (2014). Ebola virus modulates transforming growth factor beta signaling and cellular markers of mesenchyme-like transition in hepatocytes.J. Virol.889877–9892. 10.1128/jvi.01410-14
27
KogaR.KubotaM.HashiguchiT.YanagiY.OhnoS. (2018). Annexin A2 mediates the localization of measles virus matrix protein at the plasma membrane.J. Virol.92:e00181-18. 10.1128/jvi.00181-18
28
KuehnlA.MusiolA.RaabeC. A.RescherU. (2016). Emerging functions as host cell factors - an encyclopedia of annexin-pathogen interactions.Biol. Chem.397949–959. 10.1515/hsz-2016-0183
29
KuhnJ. H. (2017). Guide to the correct use of filoviral nomenclature.Curr. Top. Microbiol. Immunol.411447–460. 10.1007/82_2017_7
30
KurosakiY.UedaM. T.NakanoY.YasudaJ.KoyanagiY.SatoK.et al (2018). Different effects of two mutations on the infectivity of Ebola virus glycoprotein in nine mammalian species.J. Gen. Virol.99181–186. 10.1099/jgv.0.000999
31
LeeJ. H.YiL.LiJ.SchweitzerK.BorgmannM.NaumannM.et al (2013). Crystal structure and versatile functional roles of the COP9 signalosome subunit 1.Proc. Natl. Acad. Sci. U.S.A.11011845–11850. 10.1073/pnas.1302418110
32
MaY.SunJ.GuL.BaoH.ZhaoY.ShiL.et al (2017). Annexin A2 (ANXA2) interacts with nonstructural protein 1 and promotes the replication of highly pathogenic H5N1 avian influenza virus.BMC Microbiol.17:191. 10.1186/s12866-017-1097-0
33
MosadeghiR.ReichermeierK. M.WinklerM.SchreiberA.ReitsmaJ. M.ZhangY.et al (2016). Structural and kinetic analysis of the COP9-signalosome activation and the cullin-RING ubiquitin ligase deneddylation cycle.eLife5:e12102. 10.7554/eLife.12102
34
NanboA.ImaiM.WatanabeS.NodaT.TakahashiK.NeumannG.et al (2010). Ebolavirus is internalized into host cells via macropinocytosis in a viral glycoprotein-dependent manner.PLoS Pathog.6:e1001121. 10.1371/journal.ppat.1001121
35
NelsonE. V.PachecoJ. R.HumeA. J.CresseyT. N.DeflubeL. R.RuedasJ. B.et al (2017). An RNA polymerase II-driven Ebola virus minigenome system as an advanced tool for antiviral drug screening.Antiviral Res.14621–27. 10.1016/j.antiviral.2017.08.005
36
ObsilovaV.KopeckaM.KosekD.KacirovaM.KylarovaS.RezabkovaL.et al (2014). Mechanisms of the 14-3-3 protein function: regulation of protein function through conformational modulation.Physiol. Res.63(Suppl. 1), S155–S164.
37
OestereichL.RiegerT.NeumannM.BernreutherC.LehmannM.KrasemannS.et al (2014). Evaluation of antiviral efficacy of ribavirin, arbidol, and T-705 (favipiravir) in a mouse model for crimean-congo hemorrhagic fever.PLoS Negl. Trop. Dis.8:e2804. 10.1371/journal.pntd.0002804
38
OhmanT.LietzenN.ValimakiE.MelchjorsenJ.MatikainenS.NymanT. A. (2010). Cytosolic RNA recognition pathway activates 14-3-3 protein mediated signaling and caspase-dependent disruption of cytokeratin network in human keratinocytes.J. Proteome Res.91549–1564. 10.1021/pr901040u
39
OhmanT.SoderholmS.HintsanenP.ValimakiE.LietzenN.MacKintoshC.et al (2014). Phosphoproteomics combined with quantitative 14-3-3-affinity capture identifies SIRT1 and RAI as novel regulators of cytosolic double-stranded RNA recognition pathway.Mol. Cell. Proteomics132604–2617. 10.1074/mcp.M114.038968
40
RadhakrishnanV. M.PutnamC. W.MartinezJ. D. (2012). Activation of phosphatidylinositol 3-kinase (PI3K) and mitogen-activated protein kinase (MAPK) signaling and the consequent induction of transformation by overexpressed 14-3-3gamma protein require specific amino acids within 14-3-3gamma N-terminal variable region II.J. Biol. Chem.28743300–43311. 10.1074/jbc.M112.397877
41
RamakrishnanM. A. (2016). Determination of 50% endpoint titer using a simple formula.World J. Virol.585–86. 10.5501/wjv.v5.i2.85
42
RichardA. S.ZhangA.ParkS. J.FarzanM.ZongM.ChoeH. (2015). Virion-associated phosphatidylethanolamine promotes TIM1-mediated infection by Ebola, dengue, and West Nile viruses.Proc. Natl. Acad. Sci. U.S.A.11214682–14687. 10.1073/pnas.1508095112
43
RobinsonJ. D.PitcherJ. A. (2013). G protein-coupled receptor kinase 2 (GRK2) is a Rho-activated scaffold protein for the ERK MAP kinase cascade.Cell. Signal.252831–2839. 10.1016/j.cellsig.2013.08.031
44
RogersK. J.MauryW. (2018). The role of mononuclear phagocytes in Ebola virus infection.J. Leukoc. Biol.104717–727. 10.1002/jlb.4ri0518-183r
45
SakuraiY.KolokoltsovA. A.ChenC. C.TidwellM. W.BautaW. E.KlugbauerN.et al (2015). Ebola virus. Two-pore channels control Ebola virus host cell entry and are drug targets for disease treatment.Science347995–998. 10.1126/science.1258758
46
SimpsonJ. C.JoggerstB.LaketaV.VerissimoF.CetinC.ErfleH.et al (2012). Genome-wide RNAi screening identifies human proteins with a regulatory function in the early secretory pathway.Nat. Cell Biol.14764–774. 10.1038/ncb2510
47
SissokoD.LaouenanC.FolkessonE.M’LebingA. B.BeavoguiA. H.BaizeS.et al (2016). Experimental treatment with favipiravir for ebola virus disease (the JIKI Trial): a historically controlled, single-arm proof-of-concept trial in Guinea.PLoS Med.13:e1001967. 10.1371/journal.pmed.1001967
48
SpurgersK. B.AlefantisT.PeyserB. D.RuthelG. T.BergeronA. A.CostantinoJ. A.et al (2010). Identification of essential filovirion-associated host factors by serial proteomic analysis and RNAi screen.Mol. Cell. Proteomics92690–2703. 10.1074/mcp.M110.003418
49
StrongJ. E.WongG.JonesS. E.GrollaA.TheriaultS.KobingerG. P.et al (2008). Stimulation of Ebola virus production from persistent infection through activation of the Ras/MAPK pathway.Proc. Natl. Acad. Sci. U.S.A.10517982–17987. 10.1073/pnas.0809698105
50
TaguwaS.MaringerK.LiX.Bernal-RubioD.RauchJ. N.GestwickiJ. E.et al (2015). Defining Hsp70 subnetworks in dengue virus replication reveals key vulnerability in flavivirus infection.Cell1631108–1123. 10.1016/j.cell.2015.10.046
51
TakadaA. (2012). Filovirus tropism: cellular molecules for viral entry.Front. Microbiol.3:34. 10.3389/fmicb.2012.00034
52
van ZuylenW. J.DoyonP.ClementJ. F.KhanK. A.D’AmbrosioL. M.DoF.et al (2012). Proteomic profiling of the TRAF3 interactome network reveals a new role for the ER-to-Golgi transport compartments in innate immunity.PLoS Pathog.8:e1002747. 10.1371/journal.ppat.1002747
53
VashistS.UrenaL.Gonzalez-HernandezM. B.ChoiJ.de RougemontA.Rocha-PereiraJ.et al (2015). Molecular chaperone Hsp90 is a therapeutic target for noroviruses.J. Virol.896352–6363. 10.1128/jvi.00315-15
54
VomastekT.IwanickiM. P.BurackW. R.TiwariD.KumarD.ParsonsJ. T.et al (2008). Extracellular signal-regulated kinase 2 (ERK2) phosphorylation sites and docking domain on the nuclear pore complex protein Tpr cooperatively regulate ERK2-Tpr interaction.Mol. Cell. Biol.286954–6966. 10.1128/mcb.00925-08
55
WolfT.KannG.BeckerS.StephanC.BrodtH. R.de LeuwP.et al (2015). Severe ebola virus disease with vascular leakage and multiorgan failure: treatment of a patient in intensive care.Lancet3851428–1435. 10.1016/s0140-6736(14)62384-9
56
YounanP.IampietroM.NishidaA.RamanathanP.SantosR. I.DuttaM.et al (2017). Ebola virus binding to Tim-1 on T lymphocytes induces a cytokine storm.mBio8:e00845-17. 10.1128/mBio.00845-17
Summary
Keywords
Ebola virus life cycle, trVLPs, RNAi screening, glycoprotein, protein 40, immunoprecipitation
Citation
Yu D-S, Weng T-H, Hu C-Y, Wu Z-G, Li Y-H, Cheng L-F, Wu N-P, Li L-J and Yao H-P (2018) Chaperones, Membrane Trafficking and Signal Transduction Proteins Regulate Zaire Ebola Virus trVLPs and Interact With trVLP Elements. Front. Microbiol. 9:2724. doi: 10.3389/fmicb.2018.02724
Received
21 June 2018
Accepted
24 October 2018
Published
12 November 2018
Volume
9 - 2018
Edited by
Ping An, Frederick National Laboratory for Cancer Research (NIH), United States
Reviewed by
Junji Xing, Houston Methodist Research Institute, United States; Fatah Kashanchi, George Mason University, United States; Axel T. Lehrer, University of Hawai‘i at Mānoa, United States
Updates
Copyright
© 2018 Yu, Weng, Hu, Wu, Li, Cheng, Wu, Li and Yao.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Lan-Juan Li, ljli@zju.edu.cn Hang-Ping Yao, yaohangping@zju.edu.cn
†These authors have contributed equally to this work
This article was submitted to Virology, a section of the journal Frontiers in Microbiology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.