Evolution and Phylogenetic Analysis of Full-Length VP3 Genes of Eastern Mediterranean Bluetongue Virus Isolates

Bluetongue virus (BTV) is the ‘type’ species of the genus Orbivirus within the family Reoviridae. The BTV genome is composed of ten linear segments of double-stranded RNA (dsRNA), each of which codes for one of ten distinct viral proteins. Previous phylogenetic comparisons have evaluated variations in genome segment 3 (Seg-3) nucleotide sequence as way to identify the geographical origin (different topotypes) of BTV isolates. The full-length nucleotide sequence of genome Seg-3 was determined for thirty BTV isolates recovered in the eastern Mediterranean region, the Balkans and other geographic areas (Spain, India, Malaysia and Africa). These data were compared, based on molecular variability, positive-selection-analysis and maximum-likelihood phylogenetic reconstructions (using appropriate substitution models) to 24 previously published sequences, revealing their evolutionary relationships. These analyses indicate that negative selection is a major force in the evolution of BTV, restricting nucleotide variability, reducing the evolutionary rate of Seg-3 and potentially of other regions of the BTV genome. Phylogenetic analysis of the BTV-4 strains isolated over a relatively long time interval (1979–2000), in a single geographic area (Greece), showed a low level of nucleotide diversity, indicating that the virus can circulate almost unchanged for many years. These analyses also show that the recent incursions into south-eastern Europe were caused by BTV strains belonging to two different major-lineages: representing an ‘eastern’ (BTV-9, -16 and -1) and a ‘western’ (BTV-4) group/topotype. Epidemiological and phylogenetic analyses indicate that these viruses originated from a geographic area to the east and southeast of Greece (including Cyprus and the Middle East), which appears to represent an important ecological niche for the virus that is likely to represent a continuing source of future BTV incursions into Europe.     

Disruption of Specific RNA-RNA Interactions in a Double-Stranded RNA Virus Inhibits Genome Packaging and Virus Infectivity

Bluetongue virus (BTV) causes hemorrhagic disease in economically important livestock. The BTV genome is organized into ten discrete double-stranded RNA molecules (S1-S10) which have been suggested to follow a sequential packaging pathway from smallest to largest segment during virus capsid assembly. To substantiate and extend these studies, we have investigated the RNA sorting and packaging mechanisms with a new experimental approach using inhibitory oligonucleotides. Putative packaging signals present in the 3’untranslated regions of BTV segments were targeted by a number of nuclease resistant oligoribonucleotides (ORNs) and their effects on virus replication in cell culture were assessed. ORNs complementary to the 3’ UTR of BTV RNAs significantly inhibited virus replication without affecting protein synthesis. Same ORNs were found to inhibit complex formation when added to a novel RNA-RNA interaction assay which measured the formation of supramolecular complexes between and among different RNA segments. ORNs targeting the 3’UTR of BTV segment 10, the smallest RNA segment, were shown to be the most potent and deletions or substitution mutations of the targeted sequences diminished the RNA complexes and abolished the recovery of viable viruses using reverse genetics. Cell-free capsid assembly/RNA packaging assay also confirmed that the inhibitory ORNs could interfere with RNA packaging and further substitution mutations within the putative RNA packaging sequence have identified the recognition sequence concerned. Exchange of 3’UTR between segments have further demonstrated that RNA recognition was segment specific, most likely acting as part of the secondary structure of the entire genomic segment. Our data confirm that genome packaging in this segmented dsRNA virus occurs via the formation of supramolecular complexes formed by the interaction of specific sequences located in the 3’ UTRs. Additionally, the inhibition of packaging in-trans with inhibitory ORNs suggests this that interaction is a bona fide target for the design of compounds with antiviral activity.     

Widespread Reassortment Shapes the Evolution and Epidemiology of Bluetongue Virus following European Invasion

Genetic exchange by a process of genome-segment ‘reassortment’ represents an important mechanism for evolutionary change in all viruses with segmented genomes, yet in many cases a detailed understanding of its frequency and biological consequences is lacking. We provide a comprehensive assessment of reassortment in bluetongue virus (BTV), a globally important insect-borne pathogen of livestock, during recent outbreaks in Europe. Full-genome sequences were generated and analysed for over 150 isolates belonging to the different BTV serotypes that have emerged in the region over the last 5 decades. Based on this novel dataset we confirm that reassortment is a frequent process that plays an important and on-going role in evolution of the virus. We found evidence for reassortment in all ten segments without a significant bias towards any particular segment. However, we observed biases in the relative frequency at which particular segments were associated with each other during reassortment. This points to selective constraints possibly caused by functional relationships between individual proteins or genome segments and genome-wide epistatic interactions. Sites under positive selection were more likely to undergo amino acid changes in newly reassorted viruses, providing additional evidence for adaptive dynamics as a consequence of reassortment. We show that the live attenuated vaccines recently used in Europe have repeatedly reassorted with field strains, contributing to their genotypic, and potentially phenotypic, variability. The high degree of plasticity seen in the BTV genome in terms of segment origin suggests that current classification schemes that are based primarily on serotype, which is determined by only a single genome segment, are inadequate. Our work highlights the need for a better understanding of the mechanisms and epidemiological consequences of reassortment in BTV, as well as other segmented RNA viruses.     

Phylogenetic Analysis Reveals Rapid Evolutionary Dynamics in the Plant RNA Virus Genus Tobamovirus

Early studies on the evolutionary dynamics of plant RNA viruses suggested that they may evolve more slowly than their animal counterparts, sometimes dramatically so. However, these estimates were often based on an assumption of virus–host codivergence over time-scales of many millions of years that is difficult to verify. An important example are viruses of the genus Tobamovirus, where the assumption of host–virus codivergence over 100 million years has led to rate estimates in the range of ~1 × 10⁻⁸ nucleotide substitutions per site, per year. Such a low evolutionary rate is in apparent contradiction with the ability of some tobamoviruses to quickly overcome inbred genetic resistance. To resolve how rapidly molecular evolution proceeds in the tobomaviruses, we estimated rates of nucleotide substitution, times to common ancestry, and the extent of congruence between virus and host phylogenies. Using Bayesian coalescent methods applied to time-stamped sequences, we estimated mean evolutionary rates at the nucleotide and amino acid levels of between 1 × 10⁻⁵ and 1.3 × 10⁻³ substitutions per site, per year, and hence similar to those seen in a broad range of animal and plant RNA viruses. Under these rates, a conservative estimate for the time of origin of the sampled tobamoviruses is within the last 100,000 years, and hence a far more recently than proposed assuming codivergence. This is supported by our cophylogeny analysis which revealed significantly discordant evolutionary histories between the tobamoviruses and the plant families they infect.     

Sequential packaging of RNA genomic segments during the assembly of Bluetongue virus

Bluetongue virus (BTV), a member of the Orbivirus genus within the Reoviridae family, has a genome of 10 double-stranded RNA segments, with three distinct size classes. Although the packaging of the viral genome is evidently highly specific such that every virus particle contains a set of 10 RNA segments, the order and mechanism of packaging are not understood. In this study we have combined the use of a cell-free in vitro assembly system with a novel RNA–RNA interaction assay to investigate the mechanism of single-stranded (ss) RNAs packaging during nascent capsid assembly. Exclusion of single or multiple ssRNA segments in the packaging reaction or their addition in different order significantly altered the outcome and suggested a particular role for the smallest segment, S10. Our data suggests that genome packaging probably initiates with the smallest segment which triggers RNA–RNA interaction with other smaller segments forming a complex network. Subsequently, the medium to larger size ssRNAs are recruited until the complete genome is packaging into the capsid. The untranslated regions of the smallest RNA segment, S10, is critical for the instigation of this process. We suggest that the selective packaging observed in BTV may also apply to other members of the Reoviridae family.     

Molecular evolution and phylogeny of the human AIDS viruses LAV,HTLV-III,and ARV

Summary A phylogenetic tree for the human lymphadenopathy-associated virus (LAV), the human T-cell lymphotrophic virus type III (HTLV-III), and the acquired immune deficiency syndrome (AIDS)-associated retrovirus (ARV) has been constructed from comparisons of the amino acid sequences of their gag proteins. A method is proposed for estimating the divergence times among these AIDS viruses and the rates of nucleotide substitution for their RNA genomes. The analysis indicates that the LAV and HTLV-III strains diverged from one another after 1977 and that their common ancestor diverged from the ARV virus no more than 10 years earlier. Hence, the evolutionary diversity among strains of the AIDS viruses apparently has been generated within the last 20 years. It is estimated that the genome of the AIDS virus has a nucleotide substitution rate on the order of 10^–3 per site per year, with the rate in the second half of the genome being double that in the first half.     

Complete genome sequence analysis of a reference strain of bluetongue virus serotype 16

The entire genome of the reference strain of bluetongue virus (BTV) serotype 16 (strain RSArrrr/16) was sequenced (a total of 23,518 base pairs). The virus was obtained from the Orbivirus Reference Collection (ORC) at IAH, Pirbright, United Kingdom. The virus strain, which was previously provided by the Onderstepoort Veterinary Research Institute in South Africa, was originally isolated from the Indian subcontinent (Hazara, West Pakistan) in 1960. Previous phylogenetic comparisons show that BTV RNA sequences cluster according to the geographic origins of the virus isolate/lineage, identifying distinct BTV topotypes. Sequence comparisons of segments Seg-1 to Seg-10 show that RSArrrr/16 belongs to the major eastern topotype of BTV (BTV-16e) and can be regarded as a reference strain of BTV-16e for phylogenetic and molecular epidemiology studies. All 10 genome segments of RSArrrr/16 group closely with the vaccine strain of BTV-16 (RSAvvvv/16) that was derived from it, as well as those recently published for a Chinese isolate of BTV-16 (>99% nucleotide identity), suggesting a very recent common ancestry for all three viruses.     

Cell Tropism Predicts Long-term Nucleotide Substitution Rates of Mammalian RNA Viruses

The high rates of RNA virus evolution are generally attributed to replication with error-prone RNA-dependent RNA polymerases. However, these long-term nucleotide substitution rates span three orders of magnitude and do not correlate well with mutation rates or selection pressures. This substitution rate variation may be explained by differences in virus ecology or intrinsic genomic properties. We generated nucleotide substitution rate estimates for mammalian RNA viruses and compiled comparable published rates, yielding a dataset of 118 substitution rates of structural genes from 51 different species, as well as 40 rates of non-structural genes from 28 species. Through ANCOVA analyses, we evaluated the relationships between these rates and four ecological factors: target cell, transmission route, host range, infection duration; and three genomic properties: genome length, genome sense, genome segmentation. Of these seven factors, we found target cells to be the only significant predictors of viral substitution rates, with tropisms for epithelial cells or neurons (P<0.0001) as the most significant predictors. Further, one-tailed t-tests showed that viruses primarily infecting epithelial cells evolve significantly faster than neurotropic viruses (P<0.0001 and P<0.001 for the structural genes and non-structural genes, respectively). These results provide strong evidence that the fastest evolving mammalian RNA viruses infect cells with the highest turnover rates: the highly proliferative epithelial cells. Estimated viral generation times suggest that epithelial-infecting viruses replicate more quickly than viruses with different cell tropisms. Our results indicate that cell tropism is a key factor in viral evolvability.     

Genomic sequences of Australian bluetongue virus prototype serotypes reveal global relationships and possible routes of entry into Australia

Bluetongue virus (BTV) is transmitted by biting midges (Culicoides spp.). It causes disease mainly in sheep and occasionally in cattle and other species. BTV has spread into northern Europe, causing disease in sheep and cattle. The introduction of new serotypes, changes in vector species, and climate change have contributed to these changes. Ten BTV serotypes have been isolated in Australia without apparent associated disease. Simplified methods for preferential isolation of double-stranded RNA (dsRNA) and template preparation enabled high-throughput sequencing of the 10 genome segments of all Australian BTV prototype serotypes. Phylogenetic analysis reinforced the Western and Eastern topotypes previously characterized but revealed unique features of several Australian BTVs. Many of the Australian BTV genome segments (Seg-) were closely related, clustering together within the Eastern topotypes. A novel Australian topotype for Seg-5 (NS1) was identified, with taxa spread across several serotypes and over time. Seg-1, -2, -3, -4, -6, -7, -9, and -10 of BTV_2_AUS_2008 were most closely related to the cognate segments of viruses from Taiwan and Asia and not other Australian viruses, supporting the conclusion that BTV_2 entered Australia recently. The Australian BTV_15_AUS_1982 prototype was revealed to be unusual among the Australian BTV isolates, with Seg-3 and -8 distantly related to other BTV sequences from all serotypes.     

Nucleotide sequence analysis of VP7 gene of Indian isolates of bluetongue virus vis-à-vis other serotypes from different parts of the world.

Bluetongue virus (BTV), a member of genus Orbivirus, a family Reoviridae, is a non-enveloped with double shelled structure and ten segmented double stranded (ds) RNA genome. The RNA segment S7 encodes an inner capsid serogroup specific viral protein VP7. To amplify coding region of VP7 gene of BTV, new primers, forward primer (18-38 bp) and reverse primer (1156-1136 bp), were designed using VP7 gene sequences available in GenBank. This primer pair successfully amplified cell culture adapted Indian isolates of BTV belonging to two different serotypes 1 and 18. The coding sequences of two Indian isolates of BTV (BTV-1H and BTV-18B) were cloned into pPCR Script-Amp SK (+) plasmid vector and transformed into XL10-Gold Kan ultracompetent E. coli cells. The positive clones selected by blue-white screening and colony touch PCR were sequenced. The sequence analysis revealed that there was 93-97% nucleotide sequence identity in VP7 gene of three different Indian serotypes of BTV. The VP7 gene sequences of Indian isolates have comparatively less sequence homology (< 80%) with American (US), and French isolates compared to South African (SA), Australian (AUS) and Chinese (PRC) isolates. In silico restriction enzyme profile analysis of VP7 gene sequences revealed that Indian isolates of BTV-1 can be differentiated from other BTV-1 isolates reported from SA, AUS and PRC using TaqI. Similarly the Indian isolates of BTV belonging to three different serotypes can be differentiated using EcoRI, Hae III and TaqI restriction enzymes.     

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.     

Bluetongue virus non-structural protein 1 is a positive regulator of viral protein synthesis

ABSTRACT: BACKGROUND: Bluetongue virus (BTV) is a double-stranded RNA (dsRNA) virus of the Reoviridae family, which encodes its genes in ten linear dsRNA segments. BTV mRNAs are synthesised by the viral RNA-dependent RNA polymerase (RdRp) as exact plus sense copies of the genome segments. Infection of mammalian cells with BTV rapidly replaces cellular protein synthesis with viral protein synthesis, but the regulation of viral gene expression in the Orbivirus genus has not been investigated. RESULTS: Using an mRNA reporter system based on genome segment 10 of BTV fused with GFP we identify the protein characteristic of this genus, non-structural protein 1 (NS1) as sufficient to upregulate translation. The wider applicability of this phenomenon among the viral genes is demonstrated using the untranslated regions (UTRs) of BTV genome segments flanking the quantifiable Renilla luciferase ORF in chimeric mRNAs. The UTRs of viral mRNAs are shown to be determinants of the amount of protein synthesised, with the pre-expression of NS1 increasing the quantity in each case. The increased expression induced by pre-expression of NS1 is confirmed in virus infected cells by generating a replicating virus which expresses the reporter fused with genome segment 10, using reverse genetics. Moreover, NS1-mediated upregulation of expression is restricted to mRNAs which lack the cellular 3[PRIME] poly(A) sequence identifying the 3[PRIME] end as a necessary determinant in specifically increasing the translation of viral mRNA in the presence of cellular mRNA. CONCLUSIONS: NS1 is identified as a positive regulator of viral protein synthesis. We propose a model of translational regulation where NS1 upregulates the synthesis of viral proteins, including itself, and creates a positive feedback loop of NS1 expression, which rapidly increases the expression of all the viral proteins. The efficient translation of viral reporter mRNAs among cellular mRNAs can account for the observed replacement of cellular protein synthesis with viral protein synthesis during infection.     

Functional Mapping of Bluetongue Virus Proteins and Their Interactions with Host Proteins During Virus Replication

Bluetongue virus (BTV) is a double-stranded RNA (dsRNA) virus which is transmitted by blood-feeding gnats to wild and domestic ruminants, causing high morbidity and often high mortality. Partly due to this BTV has been in the forefront of molecular studies for last three decades and now represents one of the best understood viruses at the molecular and structural levels. BTV, like the other members of the Reoviridae family is a complex non-enveloped virus with seven structural proteins and a RNA genome consisting of 10 dsRNA segments of different sizes. In virus infected cells, three other virus encoded nonstructural proteins are synthesized. Significant recent advances have been made in understanding the structure–function relationships of BTV proteins and their interactions during virus assembly. By combining structural and molecular data it has been possible to make progress on the fundamental mechanisms used by the virus to invade, replicate in, and escape from, susceptible host cells. Data obtained from studies over a number of years have defined the key players in BTV entry, replication, assembly and egress. Specifically, it has been possible to determine the complex nature of the virion through three dimensional structure reconstructions; atomic structure of proteins and the internal capsid; the definition of the virus encoded enzymes required for RNA replication; the ordered assembly of the capsid shell and the protein sequestration required for it; and the role of three NS proteins in virus replication, assembly and release. Overall, this review demonstrates that the integration of structural, biochemical and molecular data is necessary to fully understand the assembly and replication of this complex RNA virus.     

Bayesian Coalescent Analysis Reveals a High Rate of Molecular Evolution in GB Virus C

GB virus C/hepatitis G (GBV-C) is an RNA virus of the family Flaviviridae. Despite replicating with an RNA-dependent RNA polymerase, some previous estimates of rates of evolutionary change in GBV-C suggest that it fixes mutations at the anomalously low rate of ∼10⁻⁷ nucleotide substitution per site, per year. However, these estimates were largely based on the assumption that GBV-C and its close relative GBV-A (New World monkey GB viruses) codiverged with their primate hosts over millions of years. Herein, we estimated the substitution rate of GBV-C using the largest set of dated GBV-C isolates compiled to date and a Bayesian coalescent approach that utilizes the year of sampling and so is independent of the assumption of codivergence. This revealed a rate of evolutionary change approximately four orders of magnitude higher than that estimated previously, in the range of 10⁻² to 10⁻³ sub/site/year, and hence in line with those previously determined for RNA viruses in general and the Flaviviridae in particular. In addition, we tested the assumption of host-virus codivergence in GBV-A by performing a reconciliation analysis of host and virus phylogenies. Strikingly, we found no statistical evidence for host-virus codivergence in GBV-A, indicating that substitution rates in the GB viruses should not be estimated from host divergence times.     

Shifts in the evolutionary rate and intensity of purifying selection between two Brassica genomes revealed by analyses of orthologous transposons and relics of a whole genome triplication

Recent sequencing of the Brassica rapa and Brassica oleracea genomes revealed extremely contrasting genomic features such as the abundance and distribution of transposable elements between the two genomes. However, whether and how these structural differentiations may have influenced the evolutionary rates of the two genomes since their split from a common ancestor are unknown. Here, we investigated and compared the rates of nucleotide substitution between two long terminal repeats (LTRs) of individual orthologous LTR‐retrotransposons, the rates of synonymous and non‐synonymous substitution among triplicated genes retained in both genomes from a shared whole genome triplication event, and the rates of genetic recombination estimated/deduced by the comparison of physical and genetic distances along chromosomes and ratios of solo LTRs to intact elements. Overall, LTR sequences and genic sequences showed more rapid nucleotide substitution in B. rapa than in B. oleracea. Synonymous substitution of triplicated genes retained from a shared whole genome triplication was detected at higher rates in B. rapa than in B. oleracea. Interestingly, non‐synonymous substitution was observed at lower rates in the former than in the latter, indicating shifted densities of purifying selection between the two genomes. In addition to evolutionary asymmetry, orthologous genes differentially regulated and/or disrupted by transposable elements between the two genomes were also characterized. Our analyses suggest that local genomic and epigenomic features, such as recombination rates and chromatin dynamics reshaped by independent proliferation of transposable elements and elimination between the two genomes, are perhaps partially the causes and partially the outcomes of the observed inter‐specific asymmetric evolution.     

Rate of divergence of cellular sequences homologous to segments of Moloney sarcoma virus.

The RNA genome of the Moloney isolate of murine sarcoma virus (M-MSV) consists of two parts--a sarcoma-specific region with no homology to known leukemia viral RNAs, and a shared region present also in Moloney murine leukemia virus RNA. Complementary DNA was isolated which was specific for each part of the M-MSV genome. The DNA of a number of mammalian species was examined for the presence of nucleotide sequences homologous with the two M-MSV regions. Both sets of viral sequences had homologous nucleotide sequences present in normal mouse cellular DNA. MSV-specific sequences found in mouse cellular DNA closely matched those nucleotide sequences found in M-MSV as seen by comparisons of thermal denaturation profiles. In all normal mouse cells tested, the cellular set of M-MSV-specific nucleotide sequences was present in DNA as one to a few copies per cell. The rate of base substitution of M-MSV nucleotide sequences was compared with the rate of evolution of both unique sequences and the hemoglobin gene of various species. Conservation of MSV-specific nucleotide sequences among species was similar to that of mouse globin gene(s) and greater than that of average unique cellular sequences. In contrast, cellular nucleotide sequences that are homologous to the M-MSV-murine leukemia virus "common" nucleotide region were present in multiple copies in mouse cells and were less well matched, as seen by reduced melting profiles of the hybrids. The cellular common nucleotide sequences diverged very rapidly during evolution, with a base substitution rate similar to that reported for some primate and avian endogenous virogenes. The observation that two sets of covalently linked viral sequences evolved at very different rates suggests that the origin of M-MSV may be different from endogenous helper viruses and that cellular sequences homologous to MSV-specific nucleotide sequences may be important to survival.     

Rates of Nucleotide Substitution in Angiosperm Mitochondrial DNA Sequences and Dates of Divergence Between Brassica and Other Angiosperm Lineages

We obtained 16 nucleotide sequences (∼1400 bp each) of the first intron of the mitochondrial (mt) gene for NADH subunit 4 (nad4) from 10 species of Brassicaceae. Using these new sequences and five published sequences from GenBank, we constructed a phylogenetic tree of the Brassicaceae species under study and showed that the rate of nucleotide substitution in the first intron of nad4 is very low, about 0.16–0.23 × 10⁻⁹ substitution per site per year, which is about half of the silent rate in exons of nad4. The ratios of substitution rates in this intron, ITS, and IGS are approximately 1:23:73, where ITS is the nuclear intergenic spacer between 18S and 25S rRNA genes and IGS is the intergenic spacer of 5S rRNA genes. A segment (335 bp) in the first intron of nad4 in Brassicaceae species that is absent in wheat was considered as a nonfunctional sequence and used to estimate the neutral rate (the rate of mutation) in mtDNA to be 0.5–0.7 × 10⁻⁹ substitution per site per year, which is about three times higher than the substitution rate in the rest of the first intron of nad4. We estimated that the dates of divergence are 170–235 million years (Myr) for the monocot–dicot split, 112–156 Myr for the Brassicaceae–Lettuce split, 14.5–20.4 Myr for the Brassica–Arabidopsis split, and 14.5–20.4 Myr for the Arabidopsis–Arabideae split. Received: 14 July 1998 / Accepted: 1 October 1998     

Codon conservation in the influenza A virus genome defines RNA packaging signals

Genome segmentation facilitates reassortment and rapid evolution of influenza A virus. However, segmentation complicates particle assembly as virions must contain all eight vRNA species to be infectious. Specific packaging signals exist that extend into the coding regions of most if not all segments, but these RNA motifs are poorly defined. We measured codon variability in a large dataset of sequences to identify areas of low nucleotide sequence variation independent of amino acid conservation in each segment. Most clusters of codons showing very little synonymous variation were located at segment termini, consistent with previous experimental data mapping packaging signals. Certain internal regions of conservation, most notably in the PA gene, may however signify previously unidentified functions in the virus genome. To experimentally test the bioinformatics analysis, we introduced synonymous mutations into conserved codons within known packaging signals and measured incorporation of the mutant segment into virus particles. Surprisingly, in most cases, single nucleotide changes dramatically reduced segment packaging. Thus our analysis identifies cis-acting sequences in the influenza virus genome at the nucleotide level. Furthermore, we propose that strain-specific differences exist in certain packaging signals, most notably the haemagglutinin gene; this finding has major implications for the evolution of pandemic viruses.