published_or_final_version / Biological Sciences / Doctoral / Doctor of Philosophy

Soybean mosaic virus (SMV; Genus Potyvirus; Family Potyviridae) is one of the most widespread viruses in soybean (Glycine max [L.] Merr.). Hutcheson, a cultivar developed in Virginia, is resistant to the common strains of SMV. However, new resistance-breaking (RB) isolates of SMV have emerged in natural infections to break the resistance of Hutcheson containing the Rsv1y allele. These RB isolates are SMV-G5 and G6-like based on the differential reactions on soybean cultivars with the Rsv1 locus, and are more G6-like based on the amino acid sequence of the coat protein (CP). The CP of the RB isolates is diverse at the amino and carboxy termini and highly conserved in the core region. RB isolates reduce the yield of susceptible cultivars and cause mottling of the seed coat. Dual infection of soybeans with SMV and BPMV increased the severity of symptoms, including plant stunting and SMV titer in comparison to single SMV inoculations. The reactions of Hutcheson and herbicide-tolerant Hutcheson RR were similar with or without herbicide application. Resistance to SMV is controlled by single dominant genes at three distinct loci, Rsv1, Rsv3 and Rsv4. The mechanisms of resistance at the Rsv3 and Rsv4 loci were investigated by tracking virus accumulation and movement over time using leaf immunoprints. The mechanisms of Rsv3 resistance include extreme resistance, hypersensitive response, or restriction to virus replication and movement, which are strain specific. The Rsv4 gene was found to function in a non-strain specific and non-necrotic manner. The mechanisms of Rsv4 resistance involve restricting both cell-to-cell and long distance movement of SMV. The Rsv1, Rsv3 and Rsv4 resistance genes exhibit a continuum of SMV-soybean interactions, and include complete susceptibility, local and systemic necrosis, restriction of virus movement (both cell-to-cell and long distance), reduction in virus accumulation, and extreme resistance with no detectable virus. Cultivars containing two genes for resistance, Rsv1 and Rsv3 or Rsv1 and Rsv4, were resistant to multiple strains of SMV tested and show great potential for gene pyramiding efforts to ensure a wider and more durable resistance to SMV in soybeans. / Ph. D.

Endogenous viral elements, or viral genomic fossils, have proven extremely valuable in the study of the macroevolution of viruses, providing important, and otherwise unobtainable, insights into the ancient origin of viruses, and how their ancestors might have co-evolved with their hosts in the distant past. This type of investigation falls within the realm of paleovirology—the study of ancient viruses. Investigations of extant viruses and paleovirological analyses, however, often give conflicting results, especially those concerning viral evolutionary rates and timescales. Reconciling these two types of analyses is a necessary step towards a better understanding of the overall long-term evolutionary dynamics of viruses. The main study system of this thesis is foamy viruses (FVs). FVs are characterised by their stable co-speciation history with their hosts, allowing their evolutionary dynamics to be modelled and investigated over various timescales. This unique evolutionary feature makes FVs one of the best subjects for connecting recent and ancient viral evolution. The work here reports the discovery of several endogenous mammalian FVs, and examines how mammalian FVs co-evolve with their hosts. Analyses reveal a co-diversifying history of the two that could be dated back to the basal radiation of eutherians more than 100 million years ago. However, a small number of ancient FV cross-species transmissions could still be found, mostly involving New World monkey FVs. Based on this extended FV-mammal co-speciation pattern, this thesis investigates the long-term evolutionary rate dynamics of FVs, and shows that the rate estimates of FV evolution appear to decrease continuously as the rate measurement timescale increases, following a power-law decay function. The work presented here also shows that this so-called 'time-dependent rate phenomenon' is in fact a pervasive evolutionary feature of all viruses, and surprisingly, the rate estimates of evolution of all viruses seem to decay at the same speed, decreasing by approximately half for every 3-fold increase in the measurement timescale. Based on this power-law rate-decay pattern, we could infer evolutionary timescales of modern-day lentiviruses that are consistent with paleovirological analyses for the first time. Finally, this thesis reports the discovery of basal FV-like endogenous retroviruses (FLERVs) in amphibian and fish genomes. Phylogenetic analyses reveal that the progenitors of ray-finned fish FLERVs co-diversify broadly with their fish hosts, but also suggest that there might have been several ancient viral cross-class transmissions, involving lobe-finned fish, shark, and frog FLERVs. Again, by using the power-law rate-decay model, analyses in this thesis suggest that this major retroviral clade has an ancient Ordovician marine origin, originating together with their jawed vertebrate hosts more than 450 million years ago. This finding implies that the origin of retroviruses as a whole must be in the early Paleozoic Era, if not earlier. The results presented here bridge ancient and recent viral evolution.

Influenza causes annual global epidemics and severe morbidity and mortality. The influenza virus evolves to escape from immune system antibodies that bind to it. The immune system produces influenza virus specific antibodies by VDJ recombination and somatic hypermutation. In this dissertation, we analyze the mechanism of influenza virus evolution and immune system dynamics using theoretical modeling and computational simulation. The first half of this thesis discusses influenza virus evolution. The epidemiological data inspires a novel sequence-based antigenic distance measure for subtypes H1N1 and H3N2 virus, which are superior to the conventional measure using hemagglutination inhibition assay. Historical influenza sequences show that the selective pressure increases charge in immunodominant epitopes of the H3 hemagglutinin influenza protein. Statistical mechanics and high-performance computing technology predict fixation tendencies of the H3N2 influenza virus by free energy calculation. We introduce the notion of entropy from physics and informatics to identify the epitope regions of H1-subtype influenza A with application to vaccine efficacy. We also use entropy to quantify selection and diversity in viruses with application to the hemagglutinin of H3N2 influenza. Using the bacterial E. coli as a model, we show the evidence for recombination contributing to the evolution of extended spectrum β-lactamases (ES-BLs) in clinical isolates. A guinea pig experiment supports the discussion on influenza virus evolution. The second half of the thesis discusses immune system dynamics. We design a two-scale model to describe correlation in B cell VDJ usage of zebrafish. We also introduce a dynamical system to model original antigenic sin in influenza. This dissertation aims to help researchers understand the interaction between influenza virus and the immune system with a quantitative approach.

Les invasions biologiques sont un fardeau économique important si elles affectent des ressources critiques pour l’alimentation, la sante humaine ou les productions agricoles. Les ravageurs de la pomme de terre sont un challenge économique important tant ce tubercule est un aliment clé dans les pays andins. Il est possible de suivre la dispersion récente de la teigne du Guatemala, T. solanivora en Amérique du Sud depuis son introduction au Vénézuela à sa propagation progressive vers le sud. Par ailleurs, les invasions récentes fournissent un modèle unique pour analyser les processus d’adaptation de tout l’écosystème receveur au nouveau venu. Cette introduction de T. solanivora et sa coexistence avec la teigne endémique Phthorimaea operculella, nous offre la possibilité d’étudier l’adaptation de populations virales inféodées à P. operculella au nouvel hôte T. solanivora. Une étude de terrain a été engagée dans les régions productrices de pomme de terre en Colombie. A partir des larves de T. solanivora collectées sur 5 sites distincts, des infections à granulovirus ont été détectées. Tous les isolats viraux sont apparentés au Phthorimaea operculella granulovirus (PhopGV) précédemment décrit. Des différences de pathogénicité envers les deux hôtes ont été observées. Une variabilité a été détectée pour certains isolats au niveau de deux marqueurs génétiques. Les populations présentant une diversité génétique s’avèrent plus pathogènes sur les deux hôtes que des populations génétiquement homogènes. Elles offrent une opportunité pour le contrôle biologique de ces ravageurs. Des populations artificielles ont été construites pour mimer des populations naturelles mélangées. Elles se comportent de la même manière d’un point de vue biologique, mais l’évolution de la fréquence des marqueurs n’est pas liée à l’efficacité biologique, ce qui suggère que des différences non détectées dans le génome pourraient être responsables de l’adaptation de l’hôte. La productivité des infections dans les deux hôtes a été étudiée car elle est la clé de voute du développement d’un agent de contrôle biologique. Les productivités sur P. operculella (1,36 à 2,69 × 108 OBs/ mg) et T. solanivora (0,48 à 3,64 × 108 OBs/mg) ne sont pas très différentes. Les populations génétiquement mélangées ne se distinguent pas des populations homogènes par leur production totale dans l’un ou l’autre des deux hôtes, cependant, les rendements (production virale/inoculum) montrent des différences claires, les populations mélangées (naturelles ou artificielles) sont plus performantes sur les deux hôtes. Aucune réduction de la pathogénicité sur l’hôte d’origine n’a été observée après plusieurs cycles de réplication de la population virale sur l’hôte alternatif. Les populations virales originellement adaptées à P. operculella ont évolué pour infecter T. solanivora. Dans les régions où les deux hôtes sont présents, les populations virales développent une stratégie pour être efficaces sur les deux hôtes. / Biological invasions constitute an important economical burden when they affect key resources for human alimentation, health or agronomic productions. Potato pests are important as this tuber is a key food source in Andean countries. The recent dispersion of the Guatemalan potato tuber moth, T. solanivora in South America can be traced back to the introduction in Venezuela, with progressive dispersion towards the South. Recent invasions provide, in addition, a unique model to analyse the process of adaptation of the whole receiving ecosystem to the new comers. This introduction of T. solanivora and its coexistence with the endemic potato tuber moth, Phthorimaea operculella offer us the possibility of studying the adaptation to T. solanivora of virus populations infeodated to the later. A survey has been carried out in the potato-producing regions of Colombia. From the T. solanivora larvae collected, granulovirus infections were detected in five different locations. All virus isolates are related to the previously described Phthorimaea operculella granulovirus (PhopGV). Differences in the pathogenicity against the two hosts were observed. Variability was detected in some isolates at two genetic markers. Genetically diverse populations appear to be more pathogenic for both hosts than genetically homogeneous populations. They provide a possible solution for the biological control of these insect pests. Artificial populations were constructed to mimic the mixed natural populations. They behave similarly from a biological point of view, but the evolution of the markers frequencies is not related to the biological efficacy, suggesting that undetected differences in the genomes could be responsible of this host adaptation. The productivity of the infections in both hosts has been studied as it constitutes a key point for the development of a biocontrol agent. The productivity in P. operculella (1.36 to 2.69 × 108 OBs/ mg) and T. solanivora (0.48 to 3.64 × 108 OBs/mg) are not very different. Genotypically mixed populations cannot be differentiated from homogeneous populations by their total production in one or the other host, however, the yields (virus output/doses to infect) show clear differences, mixed populations (natural or artificial) perform better in both hosts. No reduction in the pathogenicity for one host was observed after few cycles of replication of the virus population in the second host. Virus populations originally adapted to P. operculella had evolved to infect T. solanivora. In regions where both host are present, the populations developed a strategy to be efficient on both hosts.

Polymorphism in the Human Leucocyte Antigen (HLA) region of chromosome is the major source of host genetic variability in outcome of HIV-1 infection. However, there is limited understanding of the mechanisms underlying the beneficial effect of protective class I alleles such as HLA-B57, B27 and B51. Taking advantage of a unique cohort (SM cohort) infected with clade B’ HIV-1 through contaminated blood, in which many variables, such as the length of infection, the infecting viral strain and host genetic background are controlled, we performed a comprehensive study in order to understand HLA-B51 associated HIV-1 control. We first focused on the T cell responses against three dominant HLA-B51 restricted epitopes: GagNI9 (NANPDCKTI), Pol TV8 (TAFTIPSV) and Pol LI9 (LPPVVAKEI), and HLA-B51 associated escape mutations in these three epitopes. A sequential selection of epitope mutations (i.e., epitope Pol LI9, Pol TV8 and Gag NI9) was observed. Good control of viral load and higher CD4+ counts were significantly associated with at least one detectable T cell response to un-mutated epitopes. HLA-B51 restricted CD8+ T-cell clones, generated from the patients, could effectively inhibit HIV-1 replication when wild type epitopes are properly processed and presented. We then assessed the evolution of escape mutations under the selecting pressure of HLA-B51 CTLs in vitro by co-culturing HLA-B51 CTL clones with HIV-1 infected target cells (Virus Evolution Assay). Our data showed that three dominant HLA-B51 restricted CTL responses have driven the sequential escape mutations within the epitopes, leading to the loss of viral control, which confirmed our in vivo findings. Furthermore, applying Virus Evolution Assay, we assessed the impact of antigen sensitivity and TCR usage as well as founder virus effect on HIV-1 evolution and control. Our data suggested that antigen sensitivity plays an important role in anti-viral efficacy of CTLs; the TCR usage of CTLs has stronger effect on virus evolution. More importantly, our study highlighted the major impact of the founder virus sequence on viral control. It has been shown that HIV-1 has adapted to the T-cell responses to epitope Pol TI8 in other HLA-B51+ patient cohorts. However, in our cohort, T-cell responses targeting this epitope, with Valine at position 8 (Pol TV8), provide the hosts with a long-term protection against HIV-1 infection, because of a fine balance of efficient viral control, lower level of immune pressure and the slower rate of development of escape mutations. In addition, we assessed the ex vivo phenotypic characteristics of HLA-B51 restricted dominant T cell responses and our preliminary data indicated that the early differentiated and less senescence phenotype of CD8+ T cell responses in HIV-1 chronic infection is likely to be a result of low viral antigen exposure due to T cell driven escape. In conclusion, immune-dominant T-cell responses targeting three HLA-B51 restricted epitopes (Pol LI9, Pol TV8 and Gag NI9) could be advantageous for the host. In particular, the responses against epitope Gag NI9 with slow development of escape mutations or epitope Pol TV8 with a fine balance of moderate immune pressure and delayed escape mutations, are beneficial for long-term control of HIV-1 infection.

Introduction: RNA viruses are economically important pathogens of fish, and among these viruses, infectious pancreatic necrosis virus (IPNV) is of particular concern for the aquaculture industry, especially for farmed rainbow trout (Oncorhynchus mykiss) and Atlantic salmon (Salmo salar). This non-enveloped aquatic virus, which was first isolated in the UK in 1971, belongs to the family of Birnaviridae and has a bi-segmented dsRNA genome of about 6kb. IPNV is classified in 6 genogroups with correspondence to 10 known serotypes and an additional proposed genogroup of marine aquabirnaviruses (MABV). IPNV causes high mortality in fry and a reduced mortality in adult fish, respectively. Fish, which survive, can become carriers and this can lead to a clinical outbreak by releasing infective material into water or by vertical transmission via oocytes, milt and seminal fluids. Methods: This project aimed at determining the phylogeny and genomic changes of IPNV in Scotland by whole genome sequence analysis of IPNV isolates (diagnostic TCID50 supernatants) spanning 3 decades since 1982, using next generation sequencing technology. Viral RNA of IPNV culture supernatant (CHSE-214 and TO cell culture) was processed for next generation sequencing on an Illumina MiSeq platform. Library preparation was performed using the Nextera XT DNA Library Kit, prior to sequencing according to the manufacturer's MiSeq Reagent Kit v3 (150cycles) protocol. To optimize whole genome next generation sequencing for IPNV, we compared two RNA processing protocols, the Glasgow (GLAP) and the Goettingen protocol (GOEP) with focus on missing terminal nucleotides after a de novo genome assembly. Sequences were used to determine the phylogeny and selection pressure on the genome as well as a possible virus-host adaptation. Results: The results showed that both protocols were able to give full length genomes as well as genomes with missing terminal nucleotides. The phylogenetic analysis of 57 sequenced IPVN isolates shows that 78.95 % of the isolates group within genogroup V, which includes serogroup Sp and 5.26 % within genogroup I which includes serogroup Ja. Segment A of 15.79 % of the isolate grouped within genogroup III, which includes serotype Ca1 and Te but only 7.02 % of the segment B isolates grouped in the genogroup III. The remaining 8.77 % of segment B groups within genogroup II, containing the Ab serotype. Previous research has shown that residue substitutions at positions 217 and 221 in the major capsid protein VP2 have an impact on the virulence of the virus, leading to different virulence types: virulent (T217, A221), low virulence (P217, A221), avirulent (T217, T221) and persistent (P217, T221). Whole genome sequence results show that 58.93 % of the sequenced isolates belong to the persistent, 32.14 % to the low virulent type, only one isolate was of a virulent type and 7.15 % had not virulence assigned amino acid compositions in positions 217 and 221. The selection pressure analysis showed that especially VP2 is experiencing selection pressure in the variable region. In the VP1 protein we see two sites under positive selection pressure within specific motifs. VP5 showed positive selected sites mostly within the truncated region of the protein. Other proteins showed no particular interesting sites of selection. The codon adaptation analysis showed highest adaptation index for VP2. Besides VP5, which had an CAI index below one, therefore showing negative adaptation, other IPNV proteins had an CAI of barely above the value of 1. The dinucleotide abundance, focussing on CpG, showed that CpG is underrepresented in segment A and B. Discussion Phylogenetic analysis of the sequenced IPNV strains shows separate clustering of different genogroups. Genetic reassortment is observed in segment B showing a grouping within genogroup III and II although the segment A of these isolates was grouping exclusively within III. We found that over 50 % of the isolates belong to the persistent and over 30 % to the low virulent type, assuming that due to not sterilising vaccination these types were selected in the vaccinated population. The results from the CAI calculations indicate an adaptation of IPNV to its host. Together with the findings that CpG is underrepresented in IPNV it suggests that this leads to an immune escape. Especially since the selection pressure analysis showed positive selection in VP2 within the virulence determination sites of the protein, indicating that IPNV "tries" to downregulate immune recognition. The prevalence of mostly persistent type of isolates indicates together with the assumption of adaptation and immune escape that IPNV is evolving with the host in order to ensure survival.

The human immunodeficiency virus (HIV) epidemic is far from over. Antiretroviral therapy (ART) is effective at suppressing viral replication within a patient but it must be taken daily and is life-long. Therefore, the development of a therapy that could induce drug-free remission or constitute a functional cure is a key focus of HIV research. In this thesis I explore three mechanisms which could lead to more individuals being able to control their viraemia in the absence of ART: (1) T-cell immunity, (2) early initiation of ART, and (3) viral evolution. Firstly, a strong HIV-specific T-cell response has been linked to rare cases of spontaneous viral control, but the extent to which this arm of the immune response contributes to viral control is debated. Several types of data are used to answer this question, including the rate at which the virus evolves to escape the CD8+ T-cell response. I study the frequency of incident immune escape in the largest cohort used for this purpose to date. Secondly, some patients, with characteristics dissimilar to spontaneous HIV controllers, are able to control the virus for years after the interruption of ART that was initiated early in infection. I use mathematical models to investigate a new hypothesis for the differing outcomes of early- and late- initiated ART. Thirdly, since HIV is a relatively new infection of humans it is still adapting to its new host. Recent studies suggest that the virus could be evolving towards decreased virulence at the population level. I study whether the widespread administration of ART has the potential to alter the course of virulence evolution and might result in a further attenuated virus. I conclude by discussing the implications of these results for viral control at the individual level and also for population level epidemic control.

Highly pathogenic (HP) avian influenza viruses (AIVs) are naturally restricted to H5 and H7 subtypes with a polybasic cleavage site (CS) in hemagglutinin (HA) and any AIV with an intravenous pathogenicity index (IVPI) ≥ 1.2. Although only a few non-H5/H7 viruses fulfill the criteria of HPAIV; it remains unclear why these viruses did not spread in domestic birds. In 2012, a unique H4N2 virus with a polybasic CS 322PEKRRTR/G329 was isolated from quails in California which, however, was avirulent in chickens. This is the only known non-H5/H7 virus with four basic amino acids in the HACS. Here, we investigated the virulence of this virus in chickens after expansion of the polybasic CS by substitution of T327R (322PEKRRRR/G329) or T327K (322PEKRRKR/G329) with or without reassortment with HPAIV H5N1 and H7N7. The impact of single mutations or reassortment on virus fitness in vitro and in vivo was studied. Efficient cell culture replication of T327R/K carrying H4N2 viruses increased by treatment with trypsin, particularly in MDCK cells, and reassortment with HPAIV H5N1. Replication, virus excretion and bird-to-bird transmission of H4N2 was remarkably compromised by the CS mutations, but restored after reassortment with HPAIV H5N1, although not with HPAIV H7N7. Viruses carrying the H4-HA with or without R327 or K327 mutations and the other seven gene segments from HPAIV H5N1 exhibited high virulence and efficient transmission in chickens. Together, increasing the number of basic amino acids in the H4N2 HACS was detrimental for viral fitness particularly in vivo but compensated by reassortment with HPAIV H5N1. This may explain the absence of non-H5/H7 HPAIV in poultry.

Highly pathogenic (HP) avian influenza viruses (AIVs) are naturally restricted to H5 and H7 subtypes with a polybasic cleavage site (CS) in hemagglutinin (HA) and any AIV with an intravenous pathogenicity index (IVPI) 1.2. Although only a few non-H5/H7 viruses fulfill the criteria of HPAIV; it remains unclear why these viruses did not spread in domestic birds. In 2012, a unique H4N2 virus with a polybasic CS 322PEKRRTR/G329 was isolated from quails in California which, however, was avirulent in chickens. This is the only known non-H5/H7 virus with four basic amino acids in the HACS. Here, we investigated the virulence of this virus in chickens after expansion of the polybasic CS by substitution of T327R (322PEKRRRR/G329) or T327K (322PEKRRKR/G329) with or without reassortment with HPAIV H5N1 and H7N7. The impact of single mutations or reassortment on virus fitness in vitro and in vivo was studied. Ecient cell culture replication of T327R/K carrying H4N2 viruses increased by treatment with trypsin, particularly in MDCK cells, and reassortment with HPAIV H5N1. Replication, virus excretion and bird-to-bird transmission of H4N2 was remarkably compromised by the CS mutations, but restored after reassortment with HPAIV H5N1, although not with HPAIV H7N7. Viruses carrying the H4-HA with or without R327 or K327 mutations and the other seven gene segments from HPAIV H5N1 exhibited high virulence and ecient transmission in chickens. Together, increasing the number of basic amino acids in the H4N2 HACS was detrimental for viral fitness particularly in vivo but compensated by reassortment with HPAIV H5N1. This may explain the absence of non-H5/H7 HPAIV in poultry.