Search

International audience; The formation of molecular complexes is a key feature of intracellular signaling pathways which governs to the initiation and execution of dedicated cellular processes. Tumor Necrosis Factor (TNF) and Reactive Oxygen Species (ROS) function as signaling molecules and are both involved in balancing cell fate decision between cell survival or cell demise. As master regulators of cell signaling, they are also instrumental in controlling various cellular processes towards tissue homeostasis, innate immunity and inflammation. Interestingly, TNF and ROS are interlinked and involved in regulating each other’s production via the engagement of molecular signaling complexes. This relationship calls for detailed reviewing of both TNF-induced and ROS-producing molecular complexes in the context of regulated cell death (RCD) modalities. Here, we outline biotechnological approaches that were used to investigate the TNF- and, concerning ROS, the NADPH oxidase-related molecular complexes with an emphasis on different regulated cell death modalities. This systematic review highlights how the cell death field has benefited from both biochemical and live-cell fluorescence imaging approaches. This knowledge and established workflows are highly generalizable, can be of a broader use for any protein-complex studies, and well suited for addressing new challenges in signaling dynamics. These will help understand molecular signaling complexes as ensembles organized into signaling platforms, most likely the key sites of signaling dynamics integration toward cell fate regulation.

Protein activities depend heavily on protein complex formation and dynamic post-translational modifications, such as phosphorylation. Their dynamic nature is notoriously difficult to monitor in planta at cellular resolution, often requiring extensive optimization and high-end microscopy. Here, we generated and exploited the SYnthetic Multivalency in PLants (SYMPL)-vector set to study protein-protein interactions (PPIs) and kinase activities in planta based on phase separation. This technology enabled easy detection of inducible, binary and ternary protein-protein interactions among cytoplasmic, nuclear and plasma membrane proteins in plant cells via a robust image-based readout. Moreover, we applied the SYMPL toolbox to develop an in vivo reporter for SnRK1 kinase activity, allowing us to visualize tissue-specific, dynamic SnRK1 activation upon energy deprivation in stable transgenic Arabidopsis plants. The applications of the SYMPL cloning toolbox lay the foundation for the exploration of PPIs, phosphorylation and other post-translational modifications with unprecedented ease and sensitivity.

Secondary necrosis has long been perceived as an uncontrolled process resulting in total lysis of the apoptotic cell. Recently, it was shown that progression of apoptosis to secondary necrosis is regulated by Gasdermin E (GSDME), which requires activation by caspase-3. Although the contribution of GSDME in this context has been attributed to its pore-forming capacity, little is known about the kinetics and size characteristics of this. Here we report on the membrane permeabilizing features of GSDME by monitoring the influx and efflux of dextrans of different sizes into/from anti-Fas-treated L929sAhFas cells undergoing apoptosis-driven secondary necrosis. We found that GSDME accelerates cell lysis measured by SYTOX Blue staining but does not affect the exposure of phosphatidylserine on the plasma membrane. Furthermore, loss of GSDME expression clearly hampered the influx of fluorescently labeled dextrans while the efflux happened independently of the presence or absence of GSDME expression. Importantly, both in- and efflux of dextrans were dependent on their molecular weight. Altogether, our results demonstrate that GSDME regulates the passage of compounds together with other plasma membrane destabilizing subroutines. Supplementary Information The online version contains supplementary material available at 10.1007/s00018-021-04078-0.

Secondary necrosis has long been perceived as an uncontrolled process resulting in total lysis of the apoptotic cell. Recently, it was shown that progression of apoptosis to secondary necrosis is regulated by Gasdermin E (GSDME), which requires activation by caspase-3. Although the contribution of GSDME in this context has been attributed to its pore-forming capacity, little is known about the kinetics and size characteristics of this. Here we report on the membrane permeabilizing features of GSDME by monitoring the influx and efflux of dextrans of different sizes into/from anti-Fas stimulated L929sAhFas cells undergoing apoptosis-driven secondary necrosis. We found that GSDME accelerates cell lysis measured by SYTOX Blue staining but does not affect the exposure of phosphatidylserine on the plasma membrane. Furthermore, loss of GSDME expression clearly hampered the influx of fluorescently labeled dextrans while the efflux happened independent of the presence or absence of GSDME expression. Importantly, both in- and efflux of dextrans was dependent on their molecular weight. Altogether, our results demonstrate that GSDME regulates the passage of compounds together with other plasma membrane destabilizing subroutines.

ERK involvement in cell death remains unclear, although many studies have demonstrated the importance of ERK dynamics in determining cellular responses. To untangle ERK’s contribution in two cell death programs, we investigated ERK signaling dynamics during hFasL-induced apoptosis and TNF-induced necroptosis in L929sAhFas cells. We observed that ERK inhibition sensitizes cells to apoptosis while delaying necroptosis. By monitoring ERK activity by live-cell imaging using an improved ERK biosensor (EKAR4.0), we reported differential ERK signaling dynamics between cell survival, apoptosis, and necroptosis. We also decrypted a temporally shifted amplitude- and frequency-modulated (AM/FM) ERK activity profile in necroptosis versus apoptosis. ERK inhibition disrupted these ERK signaling dynamics and prevented TNF and IL-6 gene expression increase during TNF-induced necroptosis. Using inducible cell line for activated MLKL, the final executioner of necroptosis, we propose ERK and its distinctive necroptotic ERK activity dynamics to be positioned downstream of MLKL.

The gasdermin (GSDM) family has evolved as six gene clusters (GSDMA-E and Pejvakin, PJVK), and GSDM proteins are characterized by a unique N-terminal domain (N-GSDM). With the exception of PJVK, the N-GSDM domain is capable of executing plasma membrane permeabilization. Depending on the cell death modality, several protease-and kinase-dependent mechanisms directly regulate the activity of GSDME and GSDMD, the two most widely expressed and best-studied GSDMs. We provide an overview of all GSDMs in terms of biological function, tissue expression, activation, regulation, and structure. In-depth phylogenetic analysis reveals that GSDM genes show many gene duplications and deletions, suggesting that strong evolutionary forces and a unique position of the PJVK gene are associated with the occurrence of complex inner-ear development in vertebrates.

Direct targeting of the downstream mitogen-activated protein kinase (MAPK) pathway to suppress extracellular-regulated kinase (ERK) activation in KRAS and BRAF mutant colorectal cancer (CRC) has proven clinically unsuccessful, but promising results have been obtained with combination therapies including epidermal growth factor receptor (EGFR) inhibition. To elucidate the interplay between EGF signalling and ERK activation in tumours, we used patient-derived organoids (PDOs) from KRAS and BRAF mutant CRCs. PDOs resemble in vivo tumours, model treatment response and are compatible with live-cell microscopy. We established real-time, quantitative drug response assessment in PDOs with single-cell resolution, using our improved fluorescence resonance energy transfer (FRET)-based ERK biosensor EKAREN5. We show that oncogene-driven signalling is strikingly limited without EGFR activity and insufficient to sustain full proliferative potential. In PDOs and in vivo, upstream EGFR activity rigorously amplifies signal transduction efficiency in KRAS or BRAF mutant MAPK pathways. Our data provide a mechanistic understanding of the effectivity of EGFR inhibitors within combination therapies against KRAS and BRAF mutant CRC. Ponsioen et al. use a FRET‐based ERK biosensor EKAREN5 in patient‐derived organoids to show that EGFR activity amplifies signal transduction efficiency in KRAS or BRAF mutant MAPK pathways.

International audience; ERK1/2 involvement in cell death remains unclear, although many studies have demonstrated the importance of ERK1/2 dynamics in determining cellular re- sponses. To untangle how ERK1/2 contributes to two cell death programs, we investigated ERK1/2 signaling dynamics during hFasL-induced apoptosis and TNF-induced necroptosis in L929 cells. We observed that ERK1/2 inhibition sensi- tizes cells to apoptosis while delaying necroptosis. By monitoring ERK1/2 activity by live-cell imaging using an improved ERK1/2 biosensor (EKAR4.0), we reported differential ERK1/2 signaling dynamics between cell survival, apoptosis, and nec- roptosis. We also decrypted a temporally shifted amplitude- and frequency-modu- lated (AM/FM) ERK1/2 activity profile in necroptosis versus apoptosis. ERK1/2 inhibition, which disrupted ERK1/2 signaling dynamics, prevented TNF and IL-6 gene expression increase during TNF-induced necroptosis. Using an inducible cell line for activated MLKL, the final executioner of necroptosis, we showed ERK1/2 and its distinctive necroptotic ERK1/2 activity dynamics to be positioned downstream of MLKL.

RIPK3 kinase-mediated phosphorylation of MLKL pseudokinase is the execution event of necroptosis. Two independent reports—in Immunity (Yoon et al., 2017) and Cell (Gong et al., 2017)—reveal that MLKL affects homeostatic membrane trafficking and necroptosis-enhanced bubble formation involving interaction with the ESCRT machinery.

Monitoring of different signalling enzymes in a single assay using multiplex biosensing provides a multidimensional workspace to elucidate biological processes, signalling pathway crosstalk, and determine precise sequence of events at the single living cell level. In this study, we interrogate the complexity in cAMP/PKA-MAPK/ERK1&2 crosstalk by using multi-parameter biosensing experiments to correlate biochemical activities simultaneously in time and space. Using a single excitation wavelength dual colour FLIM method we are able to detect fluorescence lifetime images of two donors to simultaneously measure PKA and ERK1&2 kinase activities in the same cellular localization by using FRET biosensors. To this end, we excite two FRET donors mTFP1 and LSSmOrange with a 440 nm wavelength and we alleviate spectral bleed-through associated limitations with the very dim-fluorescent acceptor ShadowG for mTFP1 and the red-shifted mKate2 for LSSmOrange. The simultaneous recording of PKA and ERK1&2 kinase activities reveals concomitant EGF-mediated activations of both kinases in HeLa cells. Under these conditions the subsequent Forskolin-induced cAMP release reverses the transient increase of EGF-mediated ERK1&2 kinase activity while reinforcing PKA activation. Here we propose a validated methodology for multiparametric kinase biosensing in living cells using FRET-FLIM.

Over the last decade, our picture of cell death signals involved in experimental disease models totally shifted. Indeed, in addition to apoptosis, multiple forms of regulated necrosis have been associated with an increasing number of pathologies such as ischemia-reperfusion injury in brain, heart and kidney, inflammatory diseases, sepsis, retinal disorders, neurodegenerative diseases and infectious disorders. Especially necroptosis is currently attracting the attention of the scientific community. However, the in vivo identification of ongoing necroptosis in experimental disease conditions remains troublesome, mainly due to the lack of specific biomarkers. Initially, Receptor-Interacting Protein Kinase 1 (RIPK1) and RIPK3 kinase activity were uniquely associated with induction of necroptosis, however recent evidence suggests pleiotropic functions in cell death, inflammation and survival, obscuring a clear picture. In this review, we will present the last methodological advances for in vivo necroptosis identification and discuss past and recent data to provide an update of the so-called "necroptosis-associated pathologies".

International audience; Among biosensors, genetically-encoded FRET-based biosensors are widely used to localize and measure enzymatic activities. Kinases activities are of particular interest as their spatiotemporal regulation has become crucial for the deep understanding of cell fate decisions. This is especially the case for ERK, whose activity is a key node in signal transduction pathways and can direct the cell into various processes. There is a constant need for better tools to analyze kinases in vivo, and to detect even the slightest variations of their activities. Here we report the optimization of the previous ERK activity reporters, EKAR and EKAREV. Those tools are constituted by two fluorophores adapted for FRET experiments, which are flanking a specific substrate of ERK, and a domain able to recognize and bind this substrate when phosphorylated. The latter phosphorylation allows a conformational change of the biosensor and thus a FRET signal. We improved those biosensors with modifications of: (i) fluorophores and (ii) linkers between substrate and binding domain, resulting in new versions that exhibit broader dynamic ranges upon EGF stimulation when FRET experiments are carried out by fluorescence lifetime and ratiometric measurements. Herein, we characterize those new biosensors and discuss their observed differences that depend on their fluorescence properties.

Cell death contributes to the maintenance of homeostasis, but mounting evidence has confirmed the involvement of programmed cell death in some diseases. The concept of programmed cell death, which was coined several decades ago to refer to apoptosis, now also encompasses necroptosis, a newly characterized cell death program. Research on programmed cell death has become essential for the development of some new therapies. To study cell death signaling and its molecular mechanisms, new biochemical and fluorogenic approaches have been devised. Here, we first provide an overview of programmed cell death modes and the importance of dynamic cell death studies. Next, we focus on both apoptotic and necroptotic signaling and their mechanisms by providing a systematic review of all the methods and approaches that have been used. We emphasize the contribution of advanced approaches based on fluorescent probes, reporters, and Forster resonance energy transfer (FRET)-based biosensors for studying programmed cell death. Because apoptosis and necroptosis signaling pathways share some effectors molecules, we discuss how these new tools could be used to discriminate between apoptosis and necroptosis. We also describe how we developed specific FRET-based biosensors for detecting necroptosis. Finally, we touch on how dynamic measurement of biomolecules in living models will play a role in personalized prognosis and therapy.

The mitogen-activated protein kinase (MAPK) cascade in Xenopus oocytes exhibits an all-or-none, ultrasensitive response, which is believed to result from a positive feed-back loop. Here we describe a context where 1,10-phenanthroline slowly, but strongly, activates MAPK while it impairs protein synthesis in a zincdependent manner, abolishing any feed-back loop. The induced-MAPK response was found to be strongly ultrasensitive. We argue that underlying this behaviour is a regulation motif akin to a feed-forward loop acting in vivo.

Biological processes are intrinsically dynamic. Although traditional methods provide valuable insights for the understanding of many biological phenomena, the possibility of measuring, quantifying, and localizing proteins within a cell, a tissue, and even an embryo has revolutionized our train of thoughts and has encouraged scientists to develop molecular tools for the assessment of protein or protein complex dynamics within their physiological context. These ongoing efforts rest on the emergence of biophotonic techniques and the continuous improvement of fluorescent probes, allowing precise and reliable measurements of dynamic cellular functions. The march of the “ in vivo biochemistry” has begun, already yielding breathtaking results.

Biological processes are intrinsically dynamic. Although traditional methods provide valuable insights for the understanding of many biological phenomena, the possibility of measuring, quantifying, and localizing proteins within a cell, a tissue, and even an embryo has revolutionized our train of thoughts and has encouraged scientists to develop molecular tools for the assessment of protein or protein complex dynamics within their physiological context. These ongoing efforts rest on the emergence of biophotonic techniques and the continuous improvement of fluorescent probes, allowing precise and reliable measurements of dynamic cellular functions. The march of the "in vivo biochemistry" has begun, already yielding breathtaking results.

Genetically encoded Forster Resonance Energy Transfert (FRET) biosensors are powerful tools for monitoring spatiotemporal biochemical activities in living samples. By labeling a probe protein with a pair of fluorescent proteins, FRET measurement allows to follow a conformational change of the probe sensor to a specific activity. A very exciting challenge is to follow two FRET biosensors at the same time in the same sample and in the same cellular compartment. But the multiplex approach suffers from two limitations: (i) a spectral bleed-through of the first acceptor in the second donor emission band that depends on the concentration of the two biosensors and (ii) the multiple excitation wavelengths which necessitates sequential acquisition that is not adequate to follow fast signal changes in highly dynamic biochemical activities.Here, we report a method alleviating from both limitations. Taking advantage of the long stoke shift of LSSmOrange, we have used 440 nm single excitation wavelength of the two donor mTFP1 and LSSmOrange and a dual color FLIM to simultaneously measure two genetically encoded FRET biosensors. Moreover, thanks to the non-fluorescent acceptor sREACh for mTFP1 and of red-shifted mKate2 for LSSmOrange, we were able to neglect any spectral bleed-trough. With a dual spectral FLIM system we were able to detect fluorescence lifetime images of mTFP1 and LSSmOrange simultaneously and in the same cellular localization. We validated our approach by applying this methodology to simultaneously determine ERK and PKA activation in the same Hela cell using EKAR2G and AKAR4 biosensors respectively modified with mTFP1/sREACh and LSSmOrange/mKate2 pairs. By activating PKA with forskolin, the ERK pathway is not activated as expected. But surprisingly, by activating ERK with EGF, PKA is also activated denoting a cross talk between these two signaling pathway in Hela cell line.

Frequency-domain fluorescence lifetime imaging microscopy (FLIM) has become a commonly used technique to measure lifetimes in biological systems. However, lifetime measurements are strongly dependent on numerous experimental parameters. Here, we describe a complete calibration and characterization of a FLIM system and suggest parameter optimization for minimizing measurement errors during acquisition. We used standard fluorescent molecules and reference biological samples, exhibiting both single and multiple lifetime components, to calibrate and evaluate our frequency domain FLIM system. We identify several sources of lifetime precision degradation that may occur in FLIM measurements. Following a rigorous calibration of the system and a careful optimization of the acquisition parameters, we demonstrate fluorescence lifetime measurements accuracy and reliability. In addition, we show its potential on living cells by visualizing FRET in CHO cells. The proposed calibration and optimization protocol is suitable for the measurement of multiple lifetime components sample and is applicable to any frequency domain FLIM system. Using this method on our FLIM microscope enabled us to obtain the best fluorescence lifetime precision accessible with such a system. Microsc. Res. Tech., 2009. © 2008 Wiley-Liss, Inc.

Pathogenic circulating vaccine-derived polioviruses (cVDPVs) have become a major obstacle to the successful completion of the global polio eradication program. Most cVDPVs are recombinant between the oral poliovirus vaccine (OPV) and human enterovirus species C (HEV-C). To study the role of HEV-C sequences in the phenotype of cVDPVs, we generated a series of recombinants between a Madagascar cVDPV isolate and its parental OPV type 2 strain. Results indicated that the HEV-C sequences present in this cVDPV contribute to its characteristics, including pathogenicity, suggesting that interspecific recombination contributes to the phenotypic biodiversity of polioviruses and may favor the emergence of cVDPVs.

In combination with two photon excitation, FLIM is currently one of the best techniques to quantitatively study the subcellular localization of protein-protein interactions in living cells. An appropriate analysis procedure is crucial to obtain reliable results. TCSPC is an accurate method to measure FLIM. It is however an indirect process that requires photon decay curve fitting, using an exponential decay equation. Although choosing the number of exponential terms is essential, it is labor-intensive and time consuming. Therefore, a mono-model is usually applied to a whole image. Here we propose an algorithm, named Lichi, allowing pixel by pixel analysis based on the Deltachi(2) value. Lichi was validated using simulated photon decay curves with known lifetimes and proportions. It showed a high robustness for decay curves with more than 10(3) photons. When applied to lifetime images acquired from living cells, it resulted in a more realistic representation of the interaction maps. We developed an easy-to-use procedure for multi-model FLIM analysis, which enables optimized FRET quantification for all interaction texture studies, and is especially suitable to avoid the classical misinterpretation of heterogeneous samples.

Transcriptional mechanisms controlling human UCP3 gene expression in skeletal muscle remain poorly understood. Experiments based on plasmid electrotransfer into tibialis anterior muscle of C57/BL6 male mice were set up in order to functionally analyze the hUCP3 gene promoter. These transfection experiments showed that a 6300 bp region upstream of the transcription initiation site was sufficient to mediate maximal promoter activity. Further analyses with a series of 5(')-deleted constructs demonstrated that the hUCP3 gene minimal promoter was located between nucleotides -284 and -40. Furthermore, an essential region was identified between nucleotides -284 and -224. The analysis of this region revealed a putative response element for PPAR located between nucleotides -281 and -269. Finally, mutations of potential cis-acting elements within the hUCP3 minimal promoter showed the presence of two TATA boxes (-198/-194 and -45/-41) required for constitutive UCP3 gene expression. To our knowledge, this is the first time that molecular characterization of the UCP3 promoter has been achieved using an in vivo experimental model.

Both cell-specific and ubiquitous transcription factors in fibroblasts have been identified as critical for expression of the Col1a1 gene, which encodes the alpha1 chain of type I collagen. Here, we report that Yin Yang 1 (YY1) binds to the Col1a1 promoter immediately upstream of the TATA box, and we examine the functional implications of YY1 binding for regulation of Col1a1 gene expression in BALBc/3T3 fibroblasts. The Col1a1 promoter region spanning base pairs (bp) -56 to -9 bound purified recombinant YY1 and the corresponding binding activity in nuclear extracts was supershifted using a YY1-specific antibody. Mutation of the TATA box to TgTA enhanced YY1 complex formation. Mutation analysis revealed two YY1 core binding sites at -40/-37 bp (YY1A) and, on the reverse strand, at -32/-29 bp (YY1B) immediately adjacent to the TATA box. In transfections using Col1a1-luciferase constructs, mutation of YY1A decreased activity completely (wild-type p350 (p350wt), -222/+113 bp) or partially (p130wt, -84 bp/+13 bp), whereas mutation of YY1B blocked the expression of both promoter constructs. Cotransfection with pCMV-YY1 increased p350wt and p130wt activities by as much as 10-fold, whereas antisense YY1 decreased constitutive expression and blocked the increased activity due to pCMV-YY1 overexpression. The mTgTA constructs were devoid of activity, arguing for a requirement for cognate binding of the TATA box-binding protein (TBP). Electrophoretic mobility shift assays performed under conditions permitting TBP binding showed that recombinant TBP/TFIID and YY1 could bind to the -56/-9 bp fragment and that YY1B was the preferred site for YY1 binding. Our results indicate that YY1 binds to the Col1a1 proximal promoter and functions as a positive regulator of constitutive activity in fibroblasts. Although YY1 is not sufficient for transcriptional initiation, it is a required component of the transcription machinery in this promoter.

Recently, nanodiamond particles (ND) have emerged as a promising tool in the field of nanobiotechnology. However, studies about the impact of ND on living organisms are still limited to raw materials and primarily confined to in vitro studies. In this work, we investigated the cytotoxicity and in vivo toxicity of ND correlated with their chemical surface functionality (-OH, -NH2 or -CO2H). Two model systems have been used, human embryonic kidney 293 (HEK293) cells and Xenopus laevis embryos. Cell viability assays showed that ND were not cytotoxic to HEK293 cells for concentrations below 50 μg mL−1. Our data suggest that the cytotoxicity may be due to the affinity of cationic particles for the negatively charged cell membrane. In parallel, visual monitoring of microinjected early-stage embryos showed a potential embryotoxicity and teratogenicity for carboxylated ND-CO2H. ND seem to have a negative impact on the gastrulation and neurulation stages inducing phenotypical abnormalities and high mortality.

Between October 2001 and April 2002, five cases of acute flaccid paralysis (AFP) associated with type 2 vaccine-derived polioviruses (VDPVs) were reported in the southern province of the Republic of Madagascar. To determine viral factors that favor the emergence of these pathogenic VDPVs, we analyzed in detail their genomic and phenotypic characteristics and compared them with co-circulating enteroviruses. These VDPVs appeared to belong to two independent recombinant lineages with sequences from the type 2 strain of the oral poliovaccine (OPV) in the 5′-half of the genome and sequences derived from unidentified species C enteroviruses (HEV-C) in the 3′-half. VDPV strains showed characteristics similar to those of wild neurovirulent viruses including neurovirulence in poliovirus-receptor transgenic mice. We looked for other VDPVs and for circulating enteroviruses in 316 stools collected from healthy children living in the small area where most of the AFP cases occurred. We found vaccine PVs, two VDPVs similar to those found in AFP cases, some echoviruses, and above all, many serotypes of coxsackie A viruses belonging to HEV-C, with substantial genetic diversity. Several coxsackie viruses A17 and A13 carried nucleotide sequences closely related to the 2C and the 3Dpol coding regions of the VDPVs, respectively. There was also evidence of multiple genetic recombination events among the HEV-C resulting in numerous recombinant genotypes. This indicates that co-circulation of HEV-C and OPV strains is associated with evolution by recombination, resulting in unexpectedly extensive viral diversity in small human populations in some tropical regions. This probably contributed to the emergence of recombinant VDPVs. These findings give further insight into viral ecosystems and the evolutionary processes that shape viral biodiversity., Author Summary Following extensive vaccination campaigns using the attenuated oral polio vaccine, wild polioviruses remain endemic in only a few countries. Nevertheless, several poliomyelitis outbreaks associated with vaccine-derived polioviruses (VDPVs) were reported in different parts of the world in recent years, particularly in Madagascar in 2002. We analyzed the molecular characteristics of Madagascar VDPVs and compared them with those of co-circulating enteroviruses. These VDPVs appear to be recombinant viruses between vaccine polioviruses and human enteroviruses of species C (HEV-C) and to present phenotypic characteristics similar to those of wild polioviruses including pathogenicity. Similar VDPVs and other enteroviruses, including several HEV-C of different types, were found in the stools of healthy children living in neighboring villages to where most of the poliomyelitis cases occurred. Some HEV-Cs showed sequences closely related to those of VDPVs, indicating genetic recombination between these viruses and vaccine polioviruses. There was also evidence of multiple genetic recombination events among other HEV-C isolates resulting in numerous different genotypes. These findings indicate that co-circulation of HEV-C and vaccine polioviruses and their evolution by recombination results in unexpectedly extensive viral diversity, at least in some small human populations, probably contributing to the emergence of recombinant VDPVs. Results of this study give further insight into the world of viruses and their biodiversity.