Characterization of the Candiru antigenic complex (Bunyaviridae: Phlebovirus), a highly diverse and reassorting group of viruses affecting humans in tropical America

The genus Phlebovirus of the family Bunyaviridae consists of approximately 70 named viruses, currently assigned to nine serocomplexes (species) based on antigenic similarities. Sixteen other named viruses that show little serologic relationship to the nine recognized groups are also classified as tentative species in the genus. In an effort to develop a more precise classification system for phleboviruses, we are attempting to sequence most of the named viruses in the genus with the goal of clarifying their phylogenetic relationships. In this report, we describe the serologic and phylogenetic relationships of 13 viruses that were found to be members of the Candiru serocomplex; 6 of them cause disease in humans. Analysis of full genome sequences revealed branching inconsistencies that suggest five reassortment events, all involving the M segment, and thus appear to be natural reassortants. This high rate of reassortment illustrates the inaccuracy of a classification system based solely on antigenic relationships.     

Phylogenetic comparison of the S3 gene of United States prototype strains of bluetongue virus with that of field isolates from California.

To better define the molecular epidemiology of bluetongue virus (BTV) infection, the genetic characteristics and phylogenetic relationships of the S3 genes of the five U.S. prototype strains of BTV, the commercially available serotype 10 modified live virus vaccine, and 18 field isolates of BTV serotypes 10, 11, 13, and 17 obtained in California during 1980, 1981, 1989, and 1990 were determined. With the exception of the S3 gene of the U.S. prototype strain of BTV serotype 2 (BTV 2), these viruses had an overall sequence homology of between 95 and 100%. Phylogenetic analyses segregated the prototype U.S. BTV 2 strain to a unique branch (100% bootstrap value), whereas the rest of the viruses clustered in two main monophyletic groups that were not correlated with their serotype, year of isolation, or geographical origin. The lack of consistent association between S3 gene sequence and virus serotype likely is a consequence of reassortment of BTV gene segments during natural mixed infections of vertebrate and invertebrate hosts. The prototype strain of BTV 13, which is considered an introduction to the U.S. like BTV 2, presents an S3 gene which is highly homologous to those of some isolates of BTV 10 and especially to that of the vaccine strain. This finding strongly suggests that the U.S. prototype strain of BTV 13 is a natural reassortant. The different topologies of the phylogenetic trees of the L2 and S3 genes of the various viruses indicate that these two genome segments evolve independently. We conclude that the S3 gene segment of populations of BTV in California is formed by different consensus sequences which cocirculate and which cannot be grouped by serotype.     

Infectious salmon anaemia virus (ISAV) isolated from the ISA disease outbreaks in Chile diverged from ISAV isolates from Norway around 1996 and was disseminated around 2005, based on surface glycoprotein gene sequences

Background

All viruses in the family Bunyaviridae possess a tripartite genome, consisting of a small, a medium, and a large RNA segment. Bunyaviruses therefore possess considerable evolutionary potential, attributable to both intramolecular changes and to genome segment reassortment. Hantaviruses (family Bunyaviridae, genus Hantavirus) are known to cause human hemorrhagic fever with renal syndrome or hantavirus pulmonary syndrome. The primary reservoir host of Sin Nombre virus is the deer mouse (Peromyscus maniculatus), which is widely distributed in North America. We investigated the prevalence of intramolecular changes and of genomic reassortment among Sin Nombre viruses detected in deer mice in three western states.

Methods

Portions of the Sin Nombre virus small (S) and medium (M) RNA segments were amplified by RT-PCR from kidney, lung, liver and spleen of seropositive peromyscine rodents, principally deer mice, collected in Colorado, New Mexico and Montana from 1995 to 2007. Both a 142 nucleotide (nt) amplicon of the M segment, encoding a portion of the G2 transmembrane glycoprotein, and a 751 nt amplicon of the S segment, encoding part of the nucleocapsid protein, were cloned and sequenced from 19 deer mice and from one brush mouse (P. boylii), S RNA but not M RNA from one deer mouse, and M RNA but not S RNA from another deer mouse.

Results

Two of 20 viruses were found to be reassortants. Within virus sequences from different rodents, the average rate of synonymous substitutions among all pair-wise comparisons (π_s) was 0.378 in the M segment and 0.312 in the S segment sequences. The replacement substitution rate (π_a) was 7.0 × 10^-4 in the M segment and 17.3 × 10^-4 in the S segment sequences. The low π_a relative to π_s suggests strong purifying selection and this was confirmed by a Fu and Li analysis. The absolute rate of molecular evolution of the M segment was 6.76 × 10^-3 substitutions/site/year. The absolute age of the M segment tree was estimated to be 37 years. In the S segment the rate of molecular evolution was 1.93 × 10^-3 substitutions/site/year and the absolute age of the tree was 106 years. Assuming that mice were infected with a single Sin Nombre virus genotype, phylogenetic analyses revealed that 10% (2/20) of viruses were reassortants, similar to the 14% (6/43) found in a previous report.

Conclusion

Age estimates from both segments suggest that Sin Nombre virus has evolved within the past 37–106 years. The rates of evolutionary changes reported here suggest that Sin Nombre virus M and S segment reassortment occurs frequently in nature.     

Heavily glycosylated, highly fit SIVMne variants continue to diversify and undergo selection after transmission to a new host and they elicit early antibody dependent cellular responses but delayed neutralizing antibody responses

The evolutionary success of La Crosse virus (LACV, family Bunyaviridae) is due to its ability to adapt to changing conditions through intramolecular genetic changes and segment reassortment. Vertical transmission of LACV in mosquitoes increases the potential for segment reassortment. Studies were conducted to determine if segment reassortment was occurring in naturally infected Aedes triseriatus from Wisconsin and Minnesota in 2000, 2004, 2006 and 2007. Mosquito eggs were collected from various sites in Wisconsin and Minnesota. They were reared in the laboratory and adults were tested for LACV antigen by immunofluorescence assay. RNA was isolated from the abdomen of infected mosquitoes and portions of the small (S), medium (M) and large (L) viral genome segments were amplified by RT-PCR and sequenced. Overall, the viral sequences from 40 infected mosquitoes and 5 virus isolates were analyzed. Phylogenetic and linkage disequilibrium analyses revealed that approximately 25% of infected mosquitoes and viruses contained reassorted genome segments, suggesting that LACV segment reassortment is frequent in nature.     

Prospective of Genomics in Revealing Transmission, Reassortment and Evolution of Wildlife-Borne Avian Influenza A (H5N1) Viruses

The outbreak of highly pathogenic avian influenza (HPAI) H5N1 disease has led to significant loss of poultry and wild life and case fatality rates in humans of 60%. Wild birds are natural hosts for all avian influenza virus subtypes and over120 bird species have been reported with evidence of H5N1 infection. Influenza A viruses possess a segmented RNA genome and are characterized by frequently occurring genetic reassortment events, which play a very important role in virus evolution and the spread of novel gene constellations in immunologically naïve human and animal populations. Phylogenetic analysis of whole genome or sub-genomic sequences is a standard means for delineating genetic variation, novel reassortment events, and surveillance to trace the global transmission pathways. In this paper, special emphasis is given to the transmission and circulation of H5N1 among wild life populations, and to the reassortment events that are associated with inter-host transmission of the H5N1 viruses when they infect different hosts, such as birds, pigs and humans. In addition, we review the inter-subtype reassortment of the viral segments encoding inner proteins between the H5N1 viruses and viruses of other subtypes, such as H9N2 and H6N1. Finally, we highlight the usefulness of genomic sequences in molecular epidemiological analysis of HPAI H5N1 and the technical limitations in existing analytical methods that hinder them from playing a greater role in virological research.     

Nucleotide sequence of human influenza A/PR/8/34 segment 2.

The nucleotide sequence of RNA segment 2 of human influenza strain A/PR/8/34 has been determined. Segment 2 in 2341 nucleotides long and encodes a protein of 757 amino acids (86,500 daltons molecular weight) which is involved in RNA synthesis. Although segment 2 is identical in size to segment 1, which encodes a protein of related function, neither the nucleotide sequences of these two RNA segments nor the amino acid sequences of the encoded proteins appear to be homologous. The sequence of segment 2 completes the sequence of the virus (total 13,588 nucleotides).     

Transfection-mediated recombination of influenza A virus.

Several mechanisms, including a high mutation rate and reassortment of genes, have been found to be responsible for the variability of influenza A viruses. RNA recombination would be another mechanism leading to genetic variation; however, recombination has only rarely been reported to occur in influenza viruses. During ribonucleoprotein transfection experiments designed to generate viable influenza viruses from in vitro-synthesized RNA, we discovered several viruses which must have originated from recombination events. The ribonucleoprotein transfection system may enhance the formation of viruses which result from jumping of the viral polymerase between RNAs or from ligation of different viral RNAs. Five different recombinant viruses are described. Two of these, REC1 and REC2, contain a neuraminidase (NA) gene whose defective polyadenylation signal has been repaired via intergenic recombination; 124 and 95 nucleotides have been added, respectively. Another virus, REC5, must have originated by multiple recombination events since it contains a mosaic gene with sequences derived from the NA gene of influenza A/WSN/33 virus and the matrix, polymerase protein PB1, and NA genes of influenza A/PR/8/34 virus.     

Conserved Region of Mammalian Retrovirus RNA

The viral RNAs of various mammalian retroviruses contain highly conserved sequences close to their 3' ends. This was demonstrated by interviral molecular hybridization between fractionated viral complementary DNA (cDNA) and RNA. cDNA near the 3' end (cDNA(3')) from a rat virus (RPL strain) was fractionated by size and mixed with mouse virus RNA (Rauscher leukemia virus). No hybridization occurred with total cDNA (cDNA(total)), in agreement with previous results, but a cross-reacting sequence was found with the fractionated cDNA(3'). The sequences between 50 to 400 nucleotides from the 3' terminus of heteropolymeric RNA were most hybridizable. The rat viral cDNA(3') hybridized with mouse virus RNA more extensively than with RNA of remotely related retroviruses. The related viral sequence of the rodent viruses (mouse and rat) showed as much divergence in heteroduplex thermal denaturation profiles as did the unique sequence DNA of these two rodents. This suggests that over a period of time, rodent viruses have preserved a sequence with changes correlated to phylogenetic distance of hosts. The cross-reacting sequence of replication-competent retroviruses was conserved even in the genome of the replication-defective sarcoma virus and was also located in these genomes near the 3' end of 30S RNA. A fraction of RD114 cDNA(3'), corresponding to the conserved region, cross-hybridized extensively with RNA of a baboon endogenous virus (M7). Fractions of similar size prepared from cDNA(3') of MPMV, a primate type D virus, hybridized with M7 RNA to a lesser extent. Hybridization was not observed between Mason-Pfizer monkey virus and M7 if total cDNA's were incubated with viral RNAs. The degree of cross-reaction of the shared sequence appeared to be influenced by viral ancestral relatedness and host cell phylogenetic relationships. Thus, the strikingly high extent of cross-reaction at the conserved region between rodent viruses and simian sarcoma virus and between baboon virus and RD114 virus may reflect ancestral relatedness of the viruses. Slight cross-reaction at the site between type B and C viruses of rodents (mouse mammary tumor virus and RPL virus, 58-2T) or type C and D viruses of primates (M7, RD114, and Mason-Pfizer monkey virus) may have arisen at the conserved region through a mechanism that depends more on the phylogenetic relatedness of the host cells than on the viral type or origin. Determining the sequence of the conserved region may help elucidate this mechanism. The conserved sequences in retroviruses described here may be an important functional unit for the life cycle of many retroviruses.     

The Evolutionary Dynamics of Human Influenza B Virus

Despite their close phylogenetic relationship, type A and B influenza viruses exhibit major epidemiological differences in humans, with the latter both less common and less often associated with severe disease. However, it is unclear what processes determine the evolutionary dynamics of influenza B virus, and how influenza viruses A and B interact at the evolutionary scale. To address these questions we inferred the phylogenetic history of human influenza B virus using complete genome sequences for which the date (day) of isolation was available. By comparing the phylogenetic patterns of all eight viral segments we determined the occurrence of segment reassortment over a 30-year sampling period. An analysis of rates of nucleotide substitution and selection pressures revealed sporadic occurrences of adaptive evolution, most notably in the viral hemagglutinin and compatible with the action of antigenic drift, yet lower rates of overall and nonsynonymous nucleotide substitution compared to influenza A virus. Overall, these results led us to propose a model in which evolutionary changes within and between the antigenically distinct 'Yam88' and 'Vic87' lineages of influenza B virus are the result of changes in herd immunity, with reassortment continuously generating novel genetic variation. Additionally, we suggest that the interaction with influenza A virus may be central in shaping the evolutionary dynamics of influenza B virus, facilitating the shift of dominance between the Vic87 and the Yam88 lineages.     

The evolutionary genetics and emergence of avian influenza viruses in wild birds

We surveyed the genetic diversity among avian influenza virus (AIV) in wild birds, comprising 167 complete viral genomes from 14 bird species sampled in four locations across the United States. These isolates represented 29 type A influenza virus hemagglutinin (HA) and neuraminidase (NA) subtype combinations, with up to 26% of isolates showing evidence of mixed subtype infection. Through a phylogenetic analysis of the largest data set of AIV genomes compiled to date, we were able to document a remarkably high rate of genome reassortment, with no clear pattern of gene segment association and occasional inter-hemisphere gene segment migration and reassortment. From this, we propose that AIV in wild birds forms transient "genome constellations," continually reshuffled by reassortment, in contrast to the spread of a limited number of stable genome constellations that characterizes the evolution of mammalian-adapted influenza A viruses.     

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.     

Variation in nucleotide sequences coding for the N-terminal regions of the matrix and nonstructural proteins of influenza A viruses.

Nucleotide sequences have been determined for complementary DNA transcribed from the 3' ends of RNA segments 7 (matrix gene) and 8 (nonstructural gene) from a number of human influenza A viruses isolated over a period of 43 years and representing H0N1, H1N1, H2N2, and H3N2 subtypes. The pattern of nucleotide variation in both genes suggests that RNA segments 7 and 8 were conserved during the reassortment events which were responsible for the antigenic shifts H1N1 leads to H2N2 and H2N2 leads to H3N2. During the 23-year period between the isolation of A/PR/8/34(H0N1) and A/RI/5-/57(H2N2), substitutions have occurred at 7 of 230 nucleotides in RNA segment 7 and 13 of 220 nucleotides in RNA segment 8, and in 20 years A/RI/5-/57(H2N2) to A/Canberra Grammar/77(H3N2) substitutions have occurred at 5 of 230 nucleotides in RNA segment 7 and 12 of 220 nucleotides in RNA segment 8. These give rise to 2 of 67, 5 of 64, 1 of 67, and 5 of 64 amino acid changes, respectively. The number of nucleotide and amino acid changes observed is of the same order of magnitude as that which occurs over a comparable period of drift in RNA segments 4 and 6, which code for the variable antigenic determinants hemagglutinin and neuraminidase.     

Natural history of human respiratory syncytial virus inferred from phylogenetic analysis of the attachment (G) glycoprotein with a 60-nucleotide duplication

A total of 47 clinical samples were identified during an active surveillance program of respiratory infections in Buenos Aires (BA) (1999 to 2004) that contained sequences of human respiratory syncytial virus (HRSV) with a 60-nucleotide duplication in the attachment (G) protein gene. This duplication was analogous to that previously described for other three viruses also isolated in Buenos Aires in 1999 (A. Trento et al., J. Gen. Virol. 84:3115-3120, 2003). Phylogenetic analysis indicated that BA sequences with that duplication shared a common ancestor (dated about 1998) with other HRSV G sequences reported worldwide after 1999. The duplicated nucleotide sequence was an exact copy of the preceding 60 nucleotides in early viruses, but both copies of the duplicated segment accumulated nucleotide substitutions in more recent viruses at a rate apparently higher than in other regions of the G protein gene. The evolution of the viruses with the duplicated G segment apparently followed the overall evolutionary pattern previously described for HRSV, and this genotype has replaced other prevailing antigenic group B genotypes in Buenos Aires and other places. Thus, the duplicated segment represents a natural tag that can be used to track the dissemination and evolution of HRSV in an unprecedented setting. We have taken advantage of this situation to reexamine the molecular epidemiology of HRSV and to explore the natural history of this important human pathogen.     

GEOGRAPHIC DIFFERENCES IN SEXUAL REASSORTMENT IN RNA PHAGE

The genetic structure of natural bacteriophage populations is poorly understood. Recent metagenomic studies suggest that phage biogeography is characterized by frequent migration. Using virus samples mostly isolated in Southern California, we recently showed that very little population structure exists in segmented RNA phage of the Cystoviridae family due to frequent segment reassortment (sexual genetic mixis) between unrelated virus individuals. Here we use a larger genetic dataset to examine the structure of Cystoviridae phage isolated from three geographic locations in Southern New England. We document extensive natural variation in the physical sizes of RNA genome segments for these viruses. In addition, consistent with earlier findings, our phylogenetic analyses and calculations of linkage disequilibrium (LD) show no evidence of within‐segment recombination in wild populations. However, in contrast to the prior study, our analysis finds that reassortment of segments between individual phage plays a lesser role among cystoviruses sampled in New England, suggesting that the evolutionary importance of genetic mixis in Cystoviridae phage may vary according to geography. We discuss possible explanations for these conflicting results across the studies, such as differing local ecology and its impact on phage growth, and geographic differences in selection against hybrid phage genotypes.     

Phylogenetic analysis based on full genome sequencing of Italian tomato spotted wilt virus isolates identified in “Roggianese” sweet pepper and chilli pepper

A study aimed at defining population structure of Italian tomato spotted wilt virus (TSWV) isolates was performed. Full genome sequencing of six TSWV isolates found in two Italian regions (two from Latium: Lazio 17 and Tarquinia; and four from Calabria: PepCal 10, 12, 22 and 24) were assembled. Identity percentages in nucleotide sequence among these TSWV isolates are here provided. The six full length genome sequences were compared with other two Italian isolates (p105 and p202/3WT) already fully sequenced, as well as full TSWV genomes that could be retrieved from GenBank. Phylogenetic analysis, performed in concatenated sequences and for each gene of each genome segment (L, M and S), confirmed the presence of two clades, namely A-like and D-like. In particular, the phylogenetic tree based on segment L grouped all the newly sequenced TSWV isolates in D-like clade. In the M segment phylogenetic tree, all our TSWV isolates shifted in the A-like clade. Isolates separation was not correlated to their geographical origin in phylogenetic study of distinct ORFs encoded by the RNA S segment. In fact, in nucleocapsid-encoding phylogenetic tree, PepCal 10 and 22 grouped in an A-like clade with p105, PepCal 12 and 24 in a D-like clade with p202/3WT, whereas Lazio 17 and Tarquinia in a third well distinct group. NSs tree displayed only PepCal 10 with p105 in A-like clade, whereas PepCal 12, 22, 24 with p202/3WT in D-like subclade; and isolates from Latium grouped a separated clade adjacent to D-like isolates. Additional analysis on putative reassortment events showed that TSWV Calabrian isolates likely originated from a reassortment event in M RNA and other in S RNA with p105 as major parent. Recombination events were detected in isolates from Latium in L and S RNAs with Chinese isolates as putative major parent.     

Identification of a 30S RNA with properties of a defective type C virus in murine cells.

A novel species of 30S RNA has been detected in a variety of mouse cell lines. The 30S RNA is specifically packaged by helper-independent type C viruses propagated in such cells. Nucleic acid hybridization detects no homology between the 30SRNA and the genomic RNA of helper-independent mouse type C viruses. The properties of the 30S RNA suggest that it is a defective endogenous mouse type C virus and that it is analogous to a previously described class of defective endogenous rat type C virus, which has been shown previously to be the progenitor of Kirsten and Harvey murine sarcoma viruses.     

Evidence for the Circulation and Inter-Hemispheric Movement of the H14 Subtype Influenza A Virus

Three H14 influenza A virus (IAV) isolates recovered in 2010 during routine virus surveillance along the Mississippi Migratory Bird Flyway in Wisconsin, U.S.A. raised questions about the natural history of these rare viruses. These were the first H14 IAV isolates recovered in the Western Hemisphere and the only H14 IAV isolates recovered since the original four isolates in 1982 in Asia. Full length genomic sequencing of the 2010 H14 isolates demonstrated the hemagglutinin (HA) gene from the 1982 and 2010 H14 isolates showed 89.6% nucleotide and 95.6% amino acid similarity and phylogenetic analysis of these viruses placed them with strong support within the H14 subtype lineage. The level of genomic divergence observed between the 1982 and 2010 viruses provides evidence that the H14 HA segment was circulating undetected in hosts and was not maintained in environmental stasis. Further, the evolutionary relationship observed between 1982 H14 and the closely related H4 subtype HA segments were similar to contemporary comparisons suggesting limited adaptive divergence between these sister subtypes. The nonstructural (NS) segment of one 2010 isolate was placed in a NS clade isolated infrequently over the last several decades that includes the NS segment from a previously reported 1982 H14 isolate indicating the existence of an unidentified pool of genomic diversity. An additional neuraminidase reassortment event indicated a recent inter-hemispheric gene flow from Asia into the center of North America. These results demonstrate temporal and spatial gaps in the understanding of IAV natural history. Additionally, the reassortment history of these viruses raises concern for the inter-continental spread of IAVs and the efficacy of current IAV surveillance efforts in detecting genomic diversity of viruses circulating in wild birds.