Evolutionary History and Attenuation of Myxoma Virus on Two Continents

The attenuation of myxoma virus (MYXV) following its introduction as a biological control into the European rabbit populations of Australia and Europe is the canonical study of the evolution of virulence. However, the evolutionary genetics of this profound change in host-pathogen relationship is unknown. We describe the genome-scale evolution of MYXV covering a range of virulence grades sampled over 49 years from the parallel Australian and European epidemics, including the high-virulence progenitor strains released in the early 1950s. MYXV evolved rapidly over the sampling period, exhibiting one of the highest nucleotide substitution rates ever reported for a double-stranded DNA virus, and indicative of a relatively high mutation rate and/or a continually changing selective environment. Our comparative sequence data reveal that changes in virulence involved multiple genes, likely losses of gene function due to insertion-deletion events, and no mutations common to specific virulence grades. Hence, despite the similarity in selection pressures there are multiple genetic routes to attain either highly virulent or attenuated phenotypes in MYXV, resulting in convergence for phenotype but not genotype.     

Comparative Analysis of the Complete Genome Sequence of the California MSW Strain of Myxoma Virus Reveals Potential Host Adaptations

Myxomatosis is a rapidly lethal disease of European rabbits that is caused by myxoma virus (MYXV). The introduction of a South American strain of MYXV into the European rabbit population of Australia is the classic case of host-pathogen coevolution following cross-species transmission. The most virulent strains of MYXV for European rabbits are the Californian viruses, found in the Pacific states of the United States and the Baja Peninsula, Mexico. The natural host of Californian MYXV is the brush rabbit, Sylvilagus bachmani. We determined the complete sequence of the MSW strain of Californian MYXV and performed a comparative analysis with other MYXV genomes. The MSW genome is larger than that of the South American Lausanne (type) strain of MYXV due to an expansion of the terminal inverted repeats (TIRs) of the genome, with duplication of the M156R, M154L, M153R, M152R, and M151R genes and part of the M150R gene from the right-hand (RH) end of the genome at the left-hand (LH) TIR. Despite the extreme virulence of MSW, no novel genes were identified; five genes were disrupted by multiple indels or mutations to the ATG start codon, including two genes, M008.1L/R and M152R, with major virulence functions in European rabbits, and a sixth gene, M000.5L/R, was absent. The loss of these gene functions suggests that S. bachmani is a relatively recent host for MYXV and that duplication of virulence genes in the TIRs, gene loss, or sequence variation in other genes can compensate for the loss of M008.1L/R and M152R in infections of European rabbits.     

The Recent Establishment of North American H10 Lineage Influenza Viruses in Australian Wild Waterfowl and the Evolution of Australian Avian Influenza Viruses

Influenza A H10N7 virus with a hemagglutinin gene of North American origin was detected in Australian chickens and poultry abattoir workers in New South Wales, Australia, in 2010 and in chickens in Queensland, Australia, on a mixed chicken and domestic duck farm in 2012. We investigated their genomic origins by sequencing full and partial genomes of H10 viruses isolated from wild aquatic birds and poultry in Australia and analyzed them with all available avian influenza virus sequences from Oceania and representative viruses from North America and Eurasia. Our analysis showed that the H10N7 viruses isolated from poultry were similar to those that have been circulating since 2009 in Australian aquatic birds and that their initial transmission into Australia occurred during 2007 and 2008. The H10 viruses that appear to have developed endemicity in Australian wild aquatic birds were derived from several viruses circulating in waterfowl along various flyways. Their hemagglutinin gene was derived from aquatic birds in the western states of the United States, whereas the neuraminidase was closely related to that from viruses previously detected in waterfowl in Japan. The remaining genes were derived from Eurasian avian influenza virus lineages. Our analysis of virological data spanning 40 years in Oceania indicates that the long-term evolutionary dynamics of avian influenza viruses in Australia may be determined by climatic changes. The introduction and long-term persistence of avian influenza virus lineages were observed during periods with increased rainfall, whereas bottlenecks and extinction were observed during phases of widespread decreases in rainfall. These results extend our understanding of factors affecting the dynamics of avian influenza and provide important considerations for surveillance and disease control strategies.     

Australia as a global sink for the genetic diversity of avian influenza A virus

Most of our understanding of the ecology and evolution of avian influenza A virus (AIV) in wild birds is derived from studies conducted in the northern hemisphere on waterfowl, with a substantial bias towards dabbling ducks. However, relevant environmental conditions and patterns of avian migration and reproduction are substantially different in the southern hemisphere. Through the sequencing and analysis of 333 unique AIV genomes collected from wild birds collected over 15 years we show that Australia is a global sink for AIV diversity and not integrally linked with the Eurasian gene pool. Rather, AIV are infrequently introduced to Australia, followed by decades of isolated circulation and eventual extinction. The number of co-circulating viral lineages varies per subtype. AIV haemagglutinin (HA) subtypes that are rarely identified at duck-centric study sites (H8-12) had more detected introductions and contemporary co-circulating lineages in Australia. Combined with a lack of duck migration beyond the Australian-Papuan region, these findings suggest introductions by long-distance migratory shorebirds. In addition, on the available data we found no evidence of directional or consistent patterns in virus movement across the Australian continent. This feature corresponds to patterns of bird movement, whereby waterfowl have nomadic and erratic rainfall-dependant distributions rather than consistent intra-continental migratory routes. Finally, we detected high levels of virus gene segment reassortment, with a high diversity of AIV genome constellations across years and locations. These data, in addition to those from other studies in Africa and South America, clearly show that patterns of AIV dynamics in the Southern Hemisphere are distinct from those in the temperate north.     

Timing of transmission and the evolution of virulence of an insect virus

We used the nuclear polyhedrosis virus of the gypsy moth, Lymantria dispar, to investigate whether the timing of transmission influences the evolution of virulence. In theory, early transmission should favour rapid replication and increase virulence, while late transmission should favour slower replication and reduce virulence. We tested this prediction by subjecting one set of 10 virus lineages to early transmission (Early viruses) and another set to late transmission (Late viruses). Each lineage of virus underwent nine cycles of transmission. Virulence assays on these lineages indicated that viruses transmitted early were significantly more lethal than those transmitted late. Increased exploitation of the host appears to come at a cost, however. While Early viruses initially produced more progeny, Late viruses were ultimately more productive over the entire duration of the infection. These results illustrate fitness trade-offs associated with the evolution of virulence and indicate that milder viruses can obtain a numerical advantage when mild and harmful strains tend to infect separate hosts.     

Divergence of the mRNA targets for the Ssb proteins of bacteriophages T4 and RB69

Background

Mortality rates have differed during distemper outbreaks among free-ranging raccoons (Procyon lotor) living around a large Chicago-area zoo, and appeared higher in year 2001 than in 1998 and 2000. We hypothesized that a more lethal variant of the local Canine distemper virus (CDV) lineage had emerged in 2001, and sought the genetic basis that led to increased virulence. However, a more complex model surfaced during preliminary analyses of CDV genomic sequences in infected tissues and of virus isolated in vitro from the raccoons.

Results

Phylogenetic analyses of subgenomic CDV fusion (F) -, phosphoprotein (P) -, and complete hemagglutinin (H) – gene sequences indicated that distinct American CDV lineages caused the distemper epizootics. The 1998 outbreak was caused by viruses that are likely from an old CDV lineage that includes CDV Snyder Hill and Lederle, which are CDV strains from the early 1950's. The 2000 and 2001 viruses appear to stem from the lineage of CDV A75/17, which was isolated in the mid 1970's. Only the 2001 viruses formed large syncytia in brain and/or lung tissue, and during primary isolation in-vitro in Vero cells, demonstrating at least one phenotypic property by which they differed from the other viruses.

Conclusions

Two different American CDV lineages caused the raccoon distemper outbreaks. The 1998 viruses are genetically distant to the 2000/2001 viruses. Since CDV does not cause persistent infections, the cycling of different CDV lineages within the same locale suggests multiple reintroductions of the virus to area raccoons. Our findings establish a precedent for determining whether the perceived differences in mortality rates are actual and attributable in part to inherent differences between CDV strains arising from different CDV lineages.     

Genetic diversity of human immunodeficiency virus type 2: evidence for distinct sequence subtypes with differences in virus biology.

The virulence properties of human immunodeficiency virus type 2 (HIV-2) are known to vary significantly and to range from relative attenuation in certain individuals to high-level pathogenicity in others. These differences in clinical manifestations may, at least in part, be determined by genetic differences among infecting virus strains. Evaluation of the full spectrum of HIV-2 genetic diversity is thus a necessary first step towards understanding its molecular epidemiology, natural history of infection, and biological diversity. In this study, we have used nested PCR techniques to amplify viral sequences from the DNA of uncultured peripheral blood mononuclear cells from 12 patients with HIV-2 seroreactivity. Sequence analysis of four nonoverlapping genomic regions allowed a comprehensive analysis of HIV-2 phylogeny. The results revealed (i) the existence of five distinct and roughly equidistant evolutionary lineages of HIV-2 which, by analogy with HIV-1, have been termed sequence subtypes A to E; (ii) evidence for a mosaic HIV-2 genome, indicating that coinfection with genetically divergent strains and recombination can occur in HIV-2-infected individuals; and (iii) evidence supporting the conclusion that some of the HIV-2 subtypes may have arisen from independent introductions of genetically diverse sooty mangabey viruses into the human population. Importantly, only a subset of HIV-2 strains replicated in culture: all subtype A viruses grew to high titers, but attempts to isolate representatives of subtypes C, D, and E, as well as the majority of subtype B viruses, remained unsuccessful. Infection with all five viral subtypes was detectable by commercially available serological (Western immunoblot) assays, despite intersubtype sequence differences of up to 25% in the gag, pol, and env regions. These results indicate that the genetic and biological diversity of HIV-2 is far greater than previously appreciated and suggest that there may be subtype-specific differences in virus biology. Systematic natural history studies are needed to determine whether this heterogeneity has clinical relevance and whether the various HIV-2 subtypes differ in their in vivo pathogenicity.     

Phenotypic divergence during the invasion of Phyla canescens in Australia and France: evidence for selection-driven evolution

Rapid adaptive evolution has been advocated as a mechanism that promotes invasion. Demonstrating adaptive evolution in invasive species requires rigorous analysis of phenotypic shifts driven by selection. Here, we document selection-driven evolution of Phyla canescens , an Argentine weed, in two invaded regions (Australia and France). Invasive populations possessed similar or higher diversity than native populations, and displayed mixed lineages from different sources, suggesting that genetic bottlenecks in both countries might have been alleviated by multiple introductions. Compared to native populations, Australian populations displayed more investment in sexual reproduction, whereas French populations possessed enhanced vegetative reproduction and growth. We partitioned evolutionary forces (selection vs. stochastic events) using two independent methods. Results of both analyses suggest that the pattern of molecular and phenotypic variability among regions was consistent with selection-driven evolution, rather than stochastic events. Our findings indicate that selection has shaped the evolution of P . canescens in two different invaded regions.
Ecology Letters (2010) 13: 32–44     

Evolutionary genomics of mycovirus-related dsRNA viruses reveals cross-family horizontal gene transfer and evolution of diverse viral lineages

ABSTRACT: BACKGROUND: Double-stranded (ds) RNA fungal viruses are typically isometric single-shelled particles that are classified into three families, Totiviridae, Partitiviridae and Chrysoviridae, the members of which possess monopartite, bipartite and quadripartite genomes, respectively. Recent findings revealed that mycovirus-related dsRNA viruses are more diverse than previously recognized. Although an increasing number of viral complete genomic sequences have become available, the evolution of these diverse dsRNA viruses remains to be clarified. This is particularly so since there is little evidence for horizontal gene transfer (HGT) among dsRNA viruses. RESULTS: In this study, we report the molecular properties of two novel dsRNA mycoviruses that were isolated from a field strain of Sclerotinia sclerotiorum, Sunf-M: one is a large monopartite virus representing a distinct evolutionary lineage of dsRNA viruses; the other is a new member of the family Partitiviridae. Comprehensive phylogenetic analysis and genome comparison revealed that there are at least ten monopartite, three bipartite, one tripartite and three quadripartite lineages in the known dsRNA mycoviruses and that the multipartite lineages have possibly evolved from different monopartite dsRNA viruses. Moreover, we found that homologs of the S7 Domain, characteristic of members of the genus phytoreovirus in family Reoviridae are widely distributed in diverse dsRNA viral lineages, including chrysoviruses, endornaviruses and some unclassified dsRNA mycoviruses. We further provided evidence that multiple HGT events may have occurred among these dsRNA viruses from different families. CONCLUSIONS: Our study provides an insight into the phylogeny and evolution of mycovirus-related dsRNA viruses and reveals that the occurrence of HGT between different virus species and the development of multipartite genomes during evolution are important macroevolutionary mechanisms in dsRNA viruses.     

Evolution of viral sensing RIG-I-like receptor genes in Leporidae genera Oryctolagus, Sylvilagus, and Lepus

One of the most severe European rabbit (Oryctolagus cuniculus) pathogens is myxoma virus (MYXV), a rabbit-specific leporipoxvirus that causes the highly lethal disease myxomatosis. Other leporid genera, Sylvilagus and Lepus, encompass species with variable susceptibilities to MYXV, but these do not develop the lethal form of the disease. The protective role of the retinoic acid-inducible gene-I (RIG-I/DDX58) in sensing MYXV in nonpermissive human myeloid cells prompted the study of the RIG-I-like receptor (RLR) family evolution in the three leporid genera. This viral-sensor family also includes the melanoma differentiation-associated factor 5 (MDA5/IFIH1), and the laboratory of genetics and physiology 2 (LGP2/DHX58). Considering specifically the MYXV susceptible host (European rabbit) and one of the virus natural long-term hosts (Sylvilagus bachmani, brush rabbit), the amino acid differences of positively selected sites in RIG-I between the two species were located in the protein region responsible for viral RNA recognition and binding, the repressor domain. Such differences might play a determinant role in how MYXV is sensed. When looking for episodic selection on MDA5 and LGP2 of the eastern cottontail (Sylvilagus floridanus), we also uncovered evidence of selective pressures that might be exerted by a species-specific leporipoxvirus, the Shope fibroma virus. Finally, a putative alternative splicing case was identified in Oryctolagus and Lepus MDA5 isoforms, corresponding to the deletion of one specific exon. This study provided the first insights into the evolution of the leporid RLR gene family that helps illuminate the origins of the species-specific innate responses to pathogens and more specifically to MYXV.     

表达兔出血症病毒VP60蛋白的重组兔粘液瘤病毒的构建与初步评价

兔出血症病毒 (Rabbit hemorrhagic disease virus,RHDV) 及兔粘液瘤病毒 (Myxoma virus,MYXV) 分别引起兔出血症 (兔瘟) 和兔粘液瘤病,是两种严重危害家兔养殖业以及导致原产地欧洲野兔-穴兔 (Oryctolagus cuniculus) 种群近濒危的重要病原。VP60为构成RHDV衣壳的主要抗原蛋白。为研制能同时免疫预防该两种疫病的重组二联疫苗,本研究分别以MYXV和其复制非必需基因——胸腺激酶 (Thymidine kinase,TK) 基因为重组载体和同源重组靶基因,构建穿梭载体p7.5-VP60-GFP。将p7.5-VP60-GFP载体转染被MYXV感染的兔肾细胞株RK13,经同源重组后,在荧光显微镜下筛选出表达GFP的重组病毒,并将其命名为rMV-VP60-GFP。通过PCR和Western blotting进行重组病毒vp60基因特异性插入和表达验证结果显示,vp60和gfp基因成功插入MYXV基因组中并且可成功表达,表明成功构建了表达RHDV衣壳蛋白基因vp60的重组MYXV。动物攻毒保护试验表明,制备的重组病毒能保护家兔抵抗MYXV的致死性攻击,这为后续疫苗的研发奠定了基础。     

High adaptive variability and virus-driven selection on major histocompatibility complex (MHC) genes in invasive wild rabbits in Australia

The rabbit haemorrhagic disease virus (RHDV) was imported into Australia in 1995 as a biocontrol agent to manage one of the most successful and devastating invasive species, the European rabbit (Oryctolagus cuniculus cuniculus). During the first disease outbreaks, RHDV caused mortality rates of up to 97% and reduced Australian rabbit numbers to very low levels. However, recently increased genetic resistance to RHDV and strong population growth has been reported. Major histocompatibility complex (MHC) class I immune genes are important for immune responses against viruses, and a high MHC variability is thought to be crucial in adaptive processes under pathogen-driven selection. We asked whether strong population bottlenecks and presumed genetic drift would have led to low MHC variability in wild Australian rabbits, and if the retained MHC variability was enough to explain the increased resistance against RHD. Despite the past bottlenecks we found a relatively high number of MHC class I sequences distributed over 2–4 loci. We identified positive selection on putative antigen-binding sites of the MHC. We detected evidence for RHDV-driven selection as one MHC supertype was negatively associated with RHD survival, fitting expectations of frequency-dependent selection. Gene duplication and pathogen-driven selection are possible (and likely) mechanisms that maintained the adaptive potential of MHC genes in Australian rabbits. Our findings not only contribute to a better understanding of the evolution of invasive species, they are also important in the light of planned future rabbit biocontrol in Australia.     

Watching Every Step of the Way: Junín Virus Attenuation Markers in the Vaccine Lineage

The Arenaviridae family includes several hemorrhagic fever viruses which are important emerging pathogens. Junín virus, a member of this family, is the etiological agent of Argentine Hemorrhagic Fever (AHF). A collaboration between the Governments of Argentina and the USA rendered the attenuated Junín virus vaccine strain Candid#1. Arenaviruses are enveloped viruses with genomes consisting of two single-stranded RNA species (L and S), each carrying two coding regions separated by a stably structured, non-coding intergenic region. Molecular characterization of the vaccine strain and of its more virulent ancestors, XJ13 (prototype) and XJ#44, allows a systematic approach for the discovery of key elements in virulence attenuation. We show comparisons of sequence information for the S RNA of the strains XJ13, XJ#44 and Candid#1 of Junín virus, along with other strains from the vaccine lineage and a set of Junín virus field strains collected at the AHF endemic area. Comparisons of nucleotide and amino acid sequences revealed different point mutations which might be linked to the attenuated phenotype. The majority of changes are consistent with a progressive attenuation of virulence between XJ13, XJ#44 and Candid#1. We propose that changes found in genomic regions with low natural variation frequencies are more likely to be associated with the virulence attenuation process. We partially sequenced field strains to analyze the genomic variability naturally occurring for Junín virus. This information, together with the sequence analysis of strains with intermediate virulence, will serve as a starting point to study the molecular bases for viral attenuation.     

Changes in population dynamics during long-term evolution of sabin type 1 poliovirus in an immunodeficient patient

The evolution of the Sabin strain of type 1 poliovirus in a hypogammaglobulinemia patient for a period of 649 days is described. Twelve poliovirus isolates from sequential stool samples encompassing days 21 to 649 after vaccination with Sabin 1 were characterized in terms of their antigenic properties, virulence in transgenic mice, sensitivity for growth at high temperatures, and differences in nucleotide sequence from the Sabin 1 strain. Poliovirus isolates from the immunodeficient patient evolved gradually toward non-temperature-sensitive and neurovirulent phenotypes, accumulating mutations at key nucleotide positions that correlated with the observed reversion to biological properties typical of wild polioviruses. Analysis of plaque-purified viruses from stool samples revealed complex genetic and evolutionary relationships between the poliovirus strains. The generation of various coevolving genetic lineages incorporating different mutations was observed at early stages of virus excretion. The main driving force for genetic diversity appeared to be the selection of mutations at attenuation sites, particularly in the 5' noncoding region and the VP1 BC loop. Recombination between virus strains from the two main lineages was observed between days 63 and 88. Genetic heterogeneity among plaque-purified viruses at each time point seemed to decrease with time, and only viruses belonging to a unique genotypic lineage were seen from day 105 after vaccination. The relevance of vaccine-derived poliovirus strains for disease surveillance and future polio immunization policies is discussed in the context of the Global Polio Eradication Initiative.     

Myxoma Virus Oncolytic Efficiency Can Be Enhanced Through Chemical or Genetic Disruption of the Actin Cytoskeleton

Myxoma virus (MYXV) is one of many animal viruses that exhibit oncolytic properties in transformed human cells. Compared to orthopoxviruses like vaccinia (VACV), MYXV spreads inefficiently, which could compromise its use in treating tumors and their associated metastases. The VACV F11 protein promotes virus exit and rapid spread by inhibiting Rho signalling, which results in a disruption of cortical actin. We have previously shown that although MYXV lacks an F11 homolog, the F11L gene can be introduced into MYXV promoting the spread of this Leporipoxvirus in natural host cells. Here we show that the F11-encoding (F11L⁺) MYXV strain replicates to higher levels in a number of human cancer cells. We also show that F11L⁺ MYXV induces better tumor control and prolonged survival of mice bearing MDA-MB-231 cancer cells. Furthermore, we show that this virus also spreads more efficiently from the site of growth in one injected tumor, to a second untreated tumor.While we focused mostly on the use of a modified MYXV we were able to show that the effects of F11 on MYXV growth in cancer cells could be mimicked through the use of pharmacological inhibition or siRNA-mediated silencing of key regulators of cortical actin (RhoA, RhoC, mDia1, or LIMK2). These data suggest that it may be possible to increase the oncolytic efficacy of wild-type MYXV using chemical inhibitors of RhoA/C or their downstream targets. Furthermore, since all viruses must overcome barriers to exit posed by structures like cortical actin, these findings suggest that the oncolytic activity of other viruses may be enhanced through similar strategies.     

The Evolution of Virulence in RNA Viruses under a Competition–Colonization Trade-Off

RNA viruses exist in large intra-host populations which display great genotypic and phenotypic diversity. We analyze a model of viral competition between two viruses infecting a constantly replenished cell pool. We assume a trade-off between the ability of the virus to colonize new cells (cell killing rate or virulence) and its local competitiveness (replicative success within coinfected cells). We characterize the conditions that allow for viral spread by means of the basic reproductive number and show that a local coexistence equilibrium exists, which is asymptotically stable. At this equilibrium, the less virulent competitor has a reproductive advantage over the more virulent colonizer reflected by a larger equilibrium population size of the competitor. The equilibria at which one virus outcompetes the other one are unstable, i.e., a second virus is always able to permanently invade. We generalize the two-virus model to multiple viral strains, each displaying a different virulence. To account for the large phenotypic diversity in viral populations, we consider a continuous spectrum of virulences and present a continuum limit of this multiple viral strains model that describes the time evolution of an initial continuous distribution of virulence without mutations. We provide a proof of the existence of solutions of the model equations, analytically assess the properties of stationary solutions, and present numerical approximations of solutions for different initial distributions. Our simulations suggest that initial continuous distributions of virulence evolve toward a distribution that is extremely skewed in favor of competitors. At equilibrium, only the least virulent part of the population survives. The discrepancy of this finding in the continuum limit with the two-virus model is attributed to the skewed equilibrium subpopulation sizes and to the transition to a continuum. Consequently, in viral quasispecies with high virulence diversity, the model predicts collective virulence attenuation. This result may contribute to understanding virulence attenuation, which has been reported in several experimental studies.     

Biotic exchange leaves detectable genomic patterns in the Australian rain forest flora

The movement (or invasion) of plant lineages from Sunda (the Malay Archipelago) into Sahul (mainland Australia) has resulted in a present-day Australian rain forest flora of mixed ancestries. Floristic integration increased during the Quaternary when continental vegetation was subjected to recurrent expansion/contraction cycles. To date, this expansion history has yet to be investigated through multi-species, landscape-level genetic analyses within tropical Northern Australia, presumably the main point of contact for Sunda lineages. Here, we characterize and compare the dynamics of 53 species of Sunda and Sahul ancestry co-distributed across the Australian Tropics and Subtropics. We use whole chloroplast genomic datasets to obtain comparable measures of species-level diversity and estimate community dynamics through time across multiple rain forest sites. Unlike Sahul-derived species, Sunda-derived species show consistently low genomic diversities, with recent accumulation rates for Sunda species being detected across all sites, confirming recent arrival and expansion across eastern Australia. A subset of Sunda-derived species with continental distributions consistently exhibited highest diversity at the most northerly site sampled, suggesting north to south colonization processes. The same species, however, differed in the levels of genomic differentiation between the Tropics and Subtropics, suggesting that continental expansion occurs at different temporal scales, with some species experiencing a northern time lag before a southern expansion along the east coast of Australia.     

Establishment and lineage replacement of H6 influenza viruses in domestic ducks in southern China

Domestic ducks in southern China act as an important reservoir for influenza viruses and have also facilitated the establishment of multiple H6 influenza virus lineages. To understand the continuing evolution of these established lineages, 297 H6 viruses isolated from domestic ducks during 2006 and 2007 were genetically and antigenically analyzed. Phylogenetic analyses showed that group II duck H6 viruses had replaced the previously predominant group I lineage and extended their geographic distribution from coastal to inland regions. Group II H6 virus showed that the genesis and development of multiple types of deletions in the neuraminidase (NA) stalk region could occur in the influenza viruses from domestic ducks. A gradual replacement of the N2 NA subtype with N6 was observed. Significant antigenic changes occurred within group II H6 viruses so that they became antigenically distinguishable from group I and gene pool viruses. Gene exchange between group II H6 viruses and the established H5N1, H9N2, or H6N1 virus lineages in poultry in the region was very limited. These findings suggest that domestic ducks can facilitate significant genetic and antigenic changes in viruses established in this host and highlight gaps in our knowledge of influenza virus ecology and even the evolutionary behavior of this virus family in its aquatic avian reservoirs.     

Molecular epidemiology of Rabbit Haemorrhagic Disease Virus in Australia: when one became many

Rabbit Haemorrhagic Disease Virus (RHDV) was introduced into Australia in 1995 as a biological control agent against the wild European rabbit (Oryctolagus cuniculus). We evaluated its evolution over a 16‐year period (1995–2011) by examining 50 isolates collected throughout Australia, as well as the original inoculum strains. Phylogenetic analysis of capsid protein VP60 sequences of the Australian isolates, compared with those sampled globally, revealed that they form a monophyletic group with the inoculum strains (CAPM V‐351 and RHDV351INOC). Strikingly, despite more than 3000 rereleases of RHDV351INOC since 1995, only a single viral lineage has sustained its transmission in the long‐term, indicative of a major competitive advantage. In addition, we find evidence for widespread viral gene flow, in which multiple lineages entered individual geographic locations, resulting in a marked turnover of viral lineages with time, as well as a continual increase in viral genetic diversity. The rate of RHDV evolution recorded in Australia ?4.0 (3.3–4.7) × 10^?3 nucleotide substitutions per site per year – was higher than previously observed in RHDV, and evidence for adaptive evolution was obtained at two VP60 residues. Finally, more intensive study of a single rabbit population (Turretfield) in South Australia provided no evidence for viral persistence between outbreaks, with genetic diversity instead generated by continual strain importation.