Constrained pattern of viral evolution in acute and early HCV infection limits viral plasticity

Cellular immune responses during acute Hepatitis C virus (HCV) and HIV infection are a known correlate of infection outcome. Viral adaptation to these responses via mutation(s) within CD8+ T-cell epitopes allows these viruses to subvert host immune control. This study examined HCV evolution in 21 HCV genotype 1-infected subjects to characterise the level of viral adaptation during acute and early HCV infection. Of the total mutations observed 25% were within described CD8+ T-cell epitopes or at viral adaptation sites. Most mutations were maintained into the chronic phase of HCV infection (75%). The lack of reversion of adaptations and high proportion of silent substitutions suggests that HCV has structural and functional limitations that constrain evolution. These results were compared to the pattern of viral evolution observed in 98 subjects during a similar phase in HIV infection from a previous study. In contrast to HCV, evolution during acute HIV infection is marked by high levels of amino acid change relative to silent substitutions, including a higher proportion of adaptations, likely reflecting strong and continued CD8+ T-cell pressure combined with greater plasticity of the virus. Understanding viral escape dynamics for these two viruses is important for effective T cell vaccine design.

HLA-B⁎27 subtype specificity determines targeting and viral evolution of a hepatitis C virus-specific CD8+ T cell epitope

Background & Aims: HLA-B⁄27 is associated with spontaneous HCV genotype 1 clearance. HLA-B⁄27-restricted CD8+ T cells target three NS5B epitopes. Two of these epitopes are dominantly targeted in the majority of HLA-B⁄27+ patients. In chronic infection, viral escape occurs consistently in these two epitopes. The third epitope (NS5B2820) was dominantly targeted in an acutely infected patient. This was in contrast, however, to the lack of recognition and viral escape in the large majority of HLA-B⁄27+ patients. Here, we set out to determine the host factors contributing to selective targeting of this epitope. Methods: Four-digit HLA class I typing and viral sequence analyses were performed in 78 HLA-B⁄27+ patients with chronic HCV genotype 1 infection. CD8+ T cell analyses were performed in a subset of patients. In addition, HLA/peptide affinity was compared for HLA-B⁄27:02 and 05. Results: The NS5B2820 epitope is only restricted by the HLA-B⁄27 subtype HLA-B⁄27:02 (that is frequent in Mediterranean populations), but not by the prototype HLA-B⁄27 subtype B⁄27:05. Indeed, the epitope is very dominant in HLA-B⁄27:02+ patients and is associated with viral escape mutations at the anchor position for HLA-binding in 12 out of 13 HLA-B⁄27:02+ chronically infected patients. Conclusions: The NS5B2820 epitope is immunodominant in the context of HLA-B⁄27:02, but is not restricted by other HLA-B⁄27 subtypes. This finding suggests an important role of HLA subtypes in the restriction of HCV-specific CD8+ responses. With minor HLA subtypes covering up to 39% of specific populations, these findings may have important implications for the selection of epitopes for global vaccines.

Asian genotypes of JC virus in Native Americans and in a Pacific Island population: Markers of viral evolution and human migration

The human polyomavirus JC (JCV) causes the central nervous system demyelinating disease progressive multifocal leukoencephalopathy. Previously, we showed that 40% of Caucasians in the United States excrete JCV in the urine as detected by PCR. We have now studied 68 Navaho from New Mexico, 25 Flathead from Montana, and 29 Chamorro from Guam. By using PCR amplification of a fragment of the VP1 gene, JCV DNA was detected in the urine of 45 (66%) Navaho, 14 (56%) Flathead, and 20 (69%) Chamorro. Genotyping of viral DNAs in these cohorts by cycle sequencing showed predominantly type 2 (Asian), rather than type 1 (European). Type 1 is the major type in the United States and Hungary. Type 2 can be further subdivided into 2A, 2B, and 2C. Type 2A is found in China and Japan. Type 2B is a subtype related to the East Asian type, and is now found in Europe and the United States. The large majority (56–89%) of strains excreted by Native Americans and Pacific Islanders were the type 2A subtype, consistent with the origin of these strains in Asia. These findings indicate that JCV infection of Native Americans predates contact with Europeans, and likely predates migration of Amerind ancestors across the Bering land bridge around 12,000–30,000 years ago. If JCV had already differentiated into stable modern genotypes and subtypes prior to first settlement, the origin of JCV in humans may date from 50,000 to 100,000 years ago or more. We conclude that JCV may have coevolved with the human species, and that it provides a convenient marker for human migrations in both prehistoric and modern times.

On the role of RNA silencing in the pathogenicity and evolution of viroids and viral satellites

Viroids and most viral satellites have small, noncoding, and highly structured RNA genomes. How they cause disease symptoms without encoding proteins and why they have characteristic secondary structures are two longstanding questions. Recent studies have shown that both viroids and satellites are capable of inducing RNA silencing, suggesting a possible role of this mechanism in the pathology and evolution of these subviral RNAs. Here we show that preventing RNA silencing in tobacco, using a silencing suppressor, greatly reduces the symptoms caused by the Y satellite of cucumber mosaic virus. Furthermore, tomato plants expressing hairpin RNA, derived from potato spindle tuber viroid, developed symptoms similar to those of potato spindle tuber viroid infection. These results provide evidence suggesting that viroids and satellites cause disease symptoms by directing RNA silencing against physiologically important host genes. We also show that viroid and satellite RNAs are significantly resistant to RNA silencing-mediated degradation, suggesting that RNA silencing is an important selection pressure shaping the evolution of the secondary structures of these pathogens.

Drosophila Adaptation to Viral Infection through Defensive Symbiont Evolution

Microbial symbionts can modulate host interactions with biotic and abiotic factors. Such interactions may affect the evolutionary trajectories of both host and symbiont. Wolbachia protects Drosophila melanogaster against several viral infections and the strength of the protection varies between variants of this endosymbiont. Since Wolbachia is maternally transmitted, its fitness depends on the fitness of its host. Therefore, Wolbachia populations may be under selection when Drosophila is subjected to viral infection. Here we show that in D. melanogaster populations selected for increased survival upon infection with Drosophila C virus there is a strong selection coefficient for specific Wolbachia variants, leading to their fixation. Flies carrying these selected Wolbachia variants have higher survival and fertility upon viral infection when compared to flies with the other variants. These findings demonstrate how the interaction of a host with pathogens shapes the genetic composition of symbiont populations. Furthermore, host adaptation can result from the evolution of its symbionts, with host and symbiont functioning as a single evolutionary unit.

Patterns of HIV-1 evolution in individuals with differing rates of CD4 T cell decline

Evolution of HIV-1 env sequences was studied in 15 seroconverting injection drug users selected for differences in the extent of CD4 T cell decline. The rates of increase of either sequence diversity at a given visit or divergence from the first seropositive visit were both higher in progressors than in nonprogressors. Viral evolution in individuals with rapid or moderate disease progression showed selection favoring nonsynonymous mutations, while nonprogressors with low viral loads selected against the nonsynonymous mutations that might have resulted in viruses with higher levels of replication. For 10 of the 15 subjects no single variant predominated over time. Evolution away from a dominant variant was followed frequently at a later time point by return to dominance of strains closely related to that variant. The observed evolutionary pattern is consistent with either selection against only the predominant virus or independent evolution occurring in different environments within the host. Differences in the level to which CD4 T cells fall in a given time period reflect not only quantitative differences in accumulation of mutations, but differences in the types of mutations that provide the best adaptation to the host environment.

Phylogenetic analysis of within -host serially-sampled viral data

The primary goal of this dissertation is the study of patterns of viral evolution inferred from serially-sampled sequence data, i.e., sequence data obtained from strains isolated at consecutive time points from a single patient or host. RNA viral populations have an extremely high genetic variability, largely due to their astronomical population sizes within host systems, high replication rate, and short generation time. It is this aspect of their evolution that demands special attention and a different approach when studying the evolutionary relationships of serially-sampled sequence data. New methods that analyze serially-sampled data were developed shortly after a groundbreaking HIV-1 study of several patients from which viruses were isolated at recurring intervals over a period of 10 or more years. These methods assume a tree-like evolutionary model, while many RNA viruses have the capacity to exchange genetic material with one another using a process called recombination. ^ A genealogy involving recombination is best described by a network structure. A more general approach was implemented in a new computational tool, Sliding MinPD, one that is mindful of the sampling times of the input sequences and that reconstructs the viral evolutionary relationships in the form of a network structure with implicit representations of recombination events. The underlying network organization reveals unique patterns of viral evolution and could help explain the emergence of disease-associated mutants and drug-resistant strains, with implications for patient prognosis and treatment strategies. In order to comprehensively test the developed methods and to carry out comparison studies with other methods, synthetic data sets are critical. Therefore, appropriate sequence generators were also developed to simulate the evolution of serially-sampled recombinant viruses, new and more through evaluation criteria for recombination detection methods were established, and three major comparison studies were performed. The newly developed tools were also applied to "real" HIV-1 sequence data and it was shown that the results represented within an evolutionary network structure can be interpreted in biologically meaningful ways. ^

Taking Multiple Infections of Cells and Recombination into Account Leads to Small Within-Host Effective-Population-Size Estimates of HIV-1

Whether HIV-1 evolution in infected individuals is dominated by deterministic or stochastic effects remains unclear because current estimates of the effective population size of HIV-1 in vivo, N-e, are widely varying. Models assuming HIV-1 evolution to be neutral estimate N-e similar to 10(2)-10(4), smaller than the inverse mutation rate of HIV-1 (similar to 10(5)), implying the predominance of stochastic forces. In contrast, a model that includes selection estimates N-e>10(5), suggesting that deterministic forces would hold sway. The consequent uncertainty in the nature of HIV-1 evolution compromises our ability to describe disease progression and outcomes of therapy. We perform detailed bit-string simulations of viral evolution that consider large genome lengths and incorporate the key evolutionary processes underlying the genomic diversification of HIV-1 in infected individuals, namely, mutation, multiple infections of cells, recombination, selection, and epistatic interactions between multiple loci. Our simulations describe quantitatively the evolution of HIV-1 diversity and divergence in patients. From comparisons of our simulations with patient data, we estimate N-e similar to 10(3)-10(4), implying predominantly stochastic evolution. Interestingly, we find that N-e and the viral generation time are correlated with the disease progression time, presenting a route to a priori prediction of disease progression in patients. Further, we show that the previous estimate of N-e>10(5) reduces as the frequencies of multiple infections of cells and recombination assumed increase. Our simulations with N-e similar to 10(3)-10(4) may be employed to estimate markers of disease progression and outcomes of therapy that depend on the evolution of viral diversity and divergence.

Influence of recombination on acquisition and reversion of immune escape and compensatory mutations in HIV-1

Following transmission, HIV-1 adapts in the new host by acquiring mutations that allow it to escape from the host immune response at multiple epitopes. It also reverts mutations associated with epitopes targeted in the transmitting host but not in the new host. Moreover, escape mutations are often associated with additional compensatory mutations that partially recover fitness costs. It is unclear whether recombination expedites this process of multi-locus adaptation. To elucidate the role of recombination, we constructed a detailed population dynamics model that integrates viral dynamics, host immune response at multiple epitopes through cytotoxic T lymphocytes, and viral evolution driven by mutation, recombination, and selection. Using this model, we compute the expected waiting time until the emergence of the strain that has gained escape and compensatory mutations against the new host's immune response, and reverted these mutations at epitopes no longer targeted. We find that depending on the underlying fitness landscape, shaped by both costs and benefits of mutations, adaptation proceeds via distinct dominant pathways with different effects of recombination, in particular distinguishing escape and reversion. When adaptation at a single epitope is involved, recombination can substantially accelerate immune escape but minimally affects reversion. When multiple epitopes are involved, recombination can accelerate or inhibit adaptation depending on the fitness landscape. Specifically, recombination tends to delay adaptation when a purely uphill fitness landscape is accessible at each epitope, and accelerate it when a fitness valley is associated with each epitope. Our study points to the importance of recombination in shaping the adaptation of HIV-1 following its transmission to new hosts, a process central to T cell-based vaccine strategies. (C) 2015 The Authors. Published by Elsevier B.V. This is an open access article under the CC BY license.

Transmission of single HIV-1 genomes and dynamics of early immune escape revealed by ultra-deep sequencing.

We used ultra-deep sequencing to obtain tens of thousands of HIV-1 sequences from regions targeted by CD8+ T lymphocytes from longitudinal samples from three acutely infected subjects, and modeled viral evolution during the critical first weeks of infection. Previous studies suggested that a single virus established productive infection, but these conclusions were tempered because of limited sampling; now, we have greatly increased our confidence in this observation through modeling the observed earliest sample diversity based on vastly more extensive sampling. Conventional sequencing of HIV-1 from acute/early infection has shown different patterns of escape at different epitopes; we investigated the earliest escapes in exquisite detail. Over 3-6 weeks, ultradeep sequencing revealed that the virus explored an extraordinary array of potential escape routes in the process of evading the earliest CD8 T-lymphocyte responses--using 454 sequencing, we identified over 50 variant forms of each targeted epitope during early immune escape, while only 2-7 variants were detected in the same samples via conventional sequencing. In contrast to the diversity seen within epitopes, non-epitope regions, including the Envelope V3 region, which was sequenced as a control in each subject, displayed very low levels of variation. In early infection, in the regions sequenced, the consensus forms did not have a fitness advantage large enough to trigger reversion to consensus amino acids in the absence of immune pressure. In one subject, a genetic bottleneck was observed, with extensive diversity at the second time point narrowing to two dominant escape forms by the third time point, all within two months of infection. Traces of immune escape were observed in the earliest samples, suggesting that immune pressure is present and effective earlier than previously reported; quantifying the loss rate of the founder virus suggests a direct role for CD8 T-lymphocyte responses in viral containment after peak viremia. Dramatic shifts in the frequencies of epitope variants during the first weeks of infection revealed a complex interplay between viral fitness and immune escape.

A dimensionless number for understanding the evolutionary dynamics of antigenically variable RNA viruses.

Antigenically variable RNA viruses are significant contributors to the burden of infectious disease worldwide. One reason for their ubiquity is their ability to escape herd immunity through rapid antigenic evolution and thereby to reinfect previously infected hosts. However, the ways in which these viruses evolve antigenically are highly diverse. Some have only limited diversity in the long-run, with every emergence of a new antigenic variant coupled with a replacement of the older variant. Other viruses rapidly accumulate antigenic diversity over time. Others still exhibit dynamics that can be considered evolutionary intermediates between these two extremes. Here, we present a theoretical framework that aims to understand these differences in evolutionary patterns by considering a virus's epidemiological dynamics in a given host population. Our framework, based on a dimensionless number, probabilistically anticipates patterns of viral antigenic diversification and thereby quantifies a virus's evolutionary potential. It is therefore similar in spirit to the basic reproduction number, the well-known dimensionless number which quantifies a pathogen's reproductive potential. We further outline how our theoretical framework can be applied to empirical viral systems, using influenza A/H3N2 as a case study. We end with predictions of our framework and work that remains to be done to further integrate viral evolutionary dynamics with disease ecology.

Effet de l'initiation du traitement antirétroviral sur la diversité virale du VIH

L’épidémie du VIH-1 dure maintenant depuis plus de 25 ans. La grande diversité génétique de ce virus est un obstacle majeur en vue de l’éradication de cette pandémie. Au cours des années, le VIH-1 a évolué en plus de cinquante sous-types ou formes recombinantes. Cette diversité génétique est influencée par diverses pressions de sélection, incluant les pressions du système immunitaire de l’hôte et les agents antirétroviraux (ARV). En effet, bien que les ARV aient considérablement réduit les taux de morbidité et de mortalité, en plus d’améliorer la qualité et l’espérance de vie des personnes atteintes du VIH-1, ces traitements sont complexes, dispendieux et amènent leur lot de toxicité pouvant mener à des concentrations plasmatiques sous-optimales pour contrôler la réplication virale. Ceci va permettre l’émergence de variantes virales portant des mutations de résistance aux ARV. Ce phénomène est encore plus complexe lorsque l’on prend en considération l’immense diversité génétique des différents sous-types. De plus, le virus du VIH est capable de persister sous forme latente dans diverses populations cellulaires, rendant ainsi son éradication extrêmement difficile. Des stratégies pouvant restreindre la diversité virale ont donc été préconisées dans le but de favoriser les réponses immunes de l’hôte pour le contrôle de l’infection et d’identifier des variantes virales offrant une meilleure cible pour des stratégies vaccinales ou immunothérapeutiques. Dans cet esprit, nous avons donc étudié, chez des sujets infectés récemment par le VIH-1, l’effet du traitement ARV précoce sur la diversité virale de la région C2V5 du gène enveloppe ainsi que sur la taille des réservoirs. En deuxième lieu, nous avons caractérisé la pression de sélection des ARV sur des souches virales de sous types variés non-B, chez des patients du Mali et du Burkina Faso afin d’évaluer les voies d’échappement viral dans un fond génétique différent du sous-type B largement prévalent en Amérique du Nord. Notre étude a démontré la présence d’une population virale très homogène et peu diversifiée dans les premières semaines suivant l’infection, qui évolue pour atteindre une diversification de +0,23% à la fin de la première année. Cette diversification est plus importante chez les sujets n’ayant pas initié de traitement. De plus, ceci s’accompagne d’un plus grand nombre de particules virales infectieuses dans les réservoirs viraux des cellules mononucléées du sang périphérique (PBMC) chez ces sujets. Ces résultats suggèrent que l’initiation précoce du traitement pourrait avoir un effet bénéfique en retardant l’évolution virale ainsi que la taille des réservoirs, ce qui pourrait supporter une réponse immune mieux ciblée et potentiellement des stratégies immunothérapeutiques permettant d’éradiquer le virus. Nous avons également suivi 801 sujets infectés par des sous-types non-B sur le point de débuter un traitement antirétroviral. Bien que la majorité des sujets ait été à un stade avancé de la maladie, plus de 75% des individus ont obtenu une charge virale indétectable après 6 mois d’ARV, témoignant de l’efficacité comparable des ARV sur les sous-types non-B et B. Toutefois, contrairement aux virus de sous-type B, nous avons observé différentes voies moléculaires de résistance chez les sous type non-B, particulièrement chez les sous-types AGK/AK/K pour lesquels les voies de résistances étaient associées de façon prédominante aux TAM2. De plus, bien que la divergence entre les virus retrouvés chez les patients d’une même région soit faible, nos analyses phylogénétiques ont permis de conclure que ces mutations de résistance se sont produites de novo et non à partir d’un ancêtre commun porteur de résistance. Cependant, notre dernière étude au Mali nous a permis d’évaluer la résistance primaire à près de 10% et des études phylogénétiques seront effectuées afin d’évaluer la circulation de ces souches résistantes dans la population. Ces études suggèrent qu’un contrôle de la réplication virale par les ARV peut freiner la diversité du VIH et ainsi ouvrir la voie à un contrôle immunologique ciblé, utilisant de nouvelles stratégies vaccinales ou immunothérapeutiques. Toutefois, une thérapie antirétrovirale sous-optimale (adhérence, toxicité) peut conduire à l’échappement virologique en favorisant l’émergence et la dissémination de souches résistantes.

Vertical transmission of hepatitis C virus: A tale of multiple outcomes

Globally, hepatitis C virus (HCV) infection affects approximately 130 million people and 3 million new infections occur annually. HCV is also recognized as an important cause of chronic liver disease in children. The absence of proofreading properties of the HCV RNA polymerase leads to a highly error prone replication process, allowing HCV to escape host immune response. The adaptive nature of HCV evolution dictates the outcome of the disease in many ways. Here, we investigated the molecular evolution of HCV in three unrelated children who acquired chronic HCV infection as a result of mother-to-child transmission, two of whom were also coinfected with HIV-1. The persistence of discrete HCV variants and their population structure were assessed using median joining network and Bayesian approaches. While patterns of viral evolution clearly differed between subjects, immune system dysfunction related to HIV coinfection or persistent HCV seronegativity stand as potential mechanisms to explain the lack of molecular evolution observed in these three cases. In contrast, treatment of HCV infection with PegIFN, which did not lead to sustained virologic responses in all 3 cases, was not associated with commensurate variations in the complexity of the variant spectrum. Finally, the differences in the degree of divergence suggest that the mode of transmission of the virus was not the main factor driving viral evolution. (C) 2013 Elsevier B. V. All rights reserved.

Analysis of energetically biased transcripts of viruses and transposable elements

RNA interference (RNAi) is a natural endogenous process by which double-stranded RNA molecules trigger potent and specific gene silencing in eukaryotic cells and is characterized by target RNA cleavage. In mammals, small interfering RNAs (siRNAs) are the trigger molecules of choice and constitute a new class of RNA-based antiviral agents. In an efficient RNAi response, the antisense strand of siRNAs must enter the RNA-induced silencing complex (RISC) in a process mediated by thermodynamic features. In this report, we hypothesize that silent mutations capable of inverting thermodynamic properties can promote resistance to siRNAs. Extensive computational analyses were used to assess whether continuous selective pressure that promotes such mutations could lead to the emergence of viral strains completely resistant to RNAi (i.e., prone to transfer only the sense strands to RISC). Based on our findings, we propose that, although synonymous mutations may produce functional resistance, this strategy cannot be systematically adopted by viruses since the longest RNAi-refractory sequence is only 10 nt long. This finding also suggests that all mRNAs display fluctuating thermodynamic landscapes and that, in terms of thermodynamic features, RNAi is a very efficient antiviral system since there will always be sites susceptible to siRNAs.

Meccanismi evolutivi dei parapoxvirus: caratterizzazione genomica di pseudocowpoxvirus e messa a punto di sistemi per lo studio delle ricombinazioni

The Poxviruses are a family of double stranded DNA (dsDNA) viruses that cause disease in many species, both vertebrate and invertebrate. Their genomes range in size from 135 to 365 kbp and show conservation in both organization and content. In particular, the central genomic regions of the chordopoxvirus subfamily (those capable of infecting vertebrates) contain 88 genes which are present in all the virus species characterised to date and which mostly occur in the same order and orientation. In contrast, however, the terminal regions of the genomes frequently contain genes that are species or genera-specific and that are not essential for the growth of the virus in vitro but instead often encode factors with important roles in vivo including modulation of the host immune response to infection and determination of the host range of the virus. The Parapoxviruses (PPV), of which Orf virus is the prototypic species, represent a genus within the chordopoxvirus subfamily of Poxviridae and are characterised by their ability to infect ruminants and humans. The genus currently contains four recognised species of virus, bovine papular stomatitis virus (BPSV) and pseudocowpox virus (PCPV) both of which infect cattle, orf virus (OV) that infects sheep and goats, and parapoxvirus of red deer in New Zealand (PVNZ). The ORFV genome has been fully sequenced, as has that of BPSV, and is ~138 kb in length encoding ~132 genes. The vast majority of these genes allow the virus to replicate in the cytoplasm of the infected host cell and therefore encode proteins involved in replication, transcription and metabolism of nucleic acids. These genes are well conserved between all known genera of poxviruses. There is however another class of genes, located at either end of the linear dsDNA genome, that encode proteins which are non-essential for replication and generally dictate host range and virulence of the virus. The non-essential genes are often the most variable within and between species of virus and therefore are potentially useful for diagnostic purposes. Given their role in subverting the host-immune response to infection they are also targets for novel therapeutics. The function of only a relatively small number of these proteins has been elucidated and there are several genes whose function still remains obscure principally because there is little similarity between them and proteins of known function in current sequence databases. It is thought that by selectively removing some of the virulence genes, or at least neutralising the proteins in some way, current vaccines could be improved. The evolution of poxviruses has been proposed to be an adaptive process involving frequent events of gene gain and loss, such that the virus co-evolves with its specific host. Gene capture or horizontal gene transfer from the host to the virus is considered an important source of new viral genes including those likely to be involved in host range and those enabling the virus to interfere with the host immune response to infection. Given the low rate of nucleotide substitution, recombination can be seen as an essential evolutionary driving force although it is likely underestimated. Recombination in poxviruses is intimately linked to DNA replication with both viral and cellular proteins participate in this recombination-dependent replication. It has been shown, in other poxvirus genera, that recombination between isolates and perhaps even between species does occur, thereby providing another mechanism for the acquisition of new genes and for the rapid evolution of viruses. Such events may result in viruses that have a selective advantage over others, for example in re-infections (a characteristic of the PPV), or in viruses that are able to jump the species barrier and infect new hosts. Sequence data related to viral strains isolated from goats suggest that possible recombination events may have occurred between OV and PCPV (Ueda et al. 2003). The recombination events are frequent during poxvirus replication and comparative genomic analysis of several poxvirus species has revealed that recombinations occur frequently on the right terminal region. Intraspecific recombination can occur between strains of the same PPV species, but also interspecific recombination can happen depending on enough sequence similarity to enable recombination between distinct PPV species. The most important pre-requisite for a successful recombination is the coinfection of the individual host by different virus strains or species. Consequently, the following factors affecting the distribution of different viruses to shared target cells need to be considered: dose of inoculated virus, time interval between inoculation of the first and the second virus, distance between the marker mutations, genetic homology. At present there are no available data on the replication dynamics of PPV in permissive and non permissive hosts and reguarding co-infetions there are no information on the interference mechanisms occurring during the simultaneous replication of viruses of different species. This work has been carried out to set up permissive substrates allowing the replication of different PPV species, in particular keratinocytes monolayers and organotypic skin cultures. Furthermore a method to isolate and expand ovine skin stem cells was has been set up to indeep further aspects of viral cellular tropism during natural infection. The study produced important data to elucidate the replication dynamics of OV and PCPV virus in vitro as well as the mechanisms of interference that can arise during co-infection with different viral species. Moreover, the analysis carried on the genomic right terminal region of PCPV 1303/05 contributed to a better knowledge of the viral genes involved in host interaction and pathogenesis as well as to locate recombination breakpoints and genetic homologies between PPV species. Taken together these data filled several crucial gaps for the study of interspecific recombinations of PPVs which are thought to be important for a better understanding of the viral evolution and to improve the biosafety of antiviral therapy and PPV-based vectors.