Evolution of the proto-MHC ancestral region: more evidence for the plesiomorphic organisation of human chromosome 9q34 region

The present day structure of the vertebrate major histocompatibility complex (MHC) and its three paralogous regions has always been a focus of interest. In a recent study, nine human anchor genes located in the MHC region were cloned from a Branchiostoma floridae (amphioxus) cosmid library. The identification and analysis of 31 surrounding genes led to the most probable model of two rounds of en bloc duplication giving rise to these regions. These events were estimated to have occurred after the cephalochordata-craniata divergence [approximately 766 million years ago (Mya)] and before the Gnathostomata radiation (approximately 528 Mya). Furthermore, it was also shown that after this large-scale duplication one of these regions, corresponding to the human 9q33-q34, had retained an ancestral organisation. In the present study, four new cosmids in the amphioxus proto-MHC region were identified by the chromosomal walking technique. These cosmids were sequenced, and their structural annotation was performed, leading to the prediction of eleven genes. Their phylogenetic relationships among species corroborate the results obtained previously and provide more evidence for the plesiomorphic state of the human chromosome 9q33-34 MHC paralogous region.An erratum to this article can be found at     

SNP discovery by amplicon sequencing and multiplex SNP genotyping in the allopolyploid species Brassica napus

Oilseed rape (Brassica napus) is an allotetraploid species consisting of two genomes, derived from B. rapa (A genome) and B. oleracea (C genome). The presence of these two genomes makes single nucleotide polymorphism (SNP) marker identification and SNP analysis more challenging than in diploid species, as for a given locus usually two versions of a DNA sequence (based on the two ancestral genomes) have to be analyzed simultaneously during SNP identification and analysis. One hundred amplicons derived from expressed sequence tag (ESTs) were analyzed to identify SNPs in a panel of oilseed rape varieties and within two sister species representing the ancestral genomes. A total of 604 SNPs were identified, averaging one SNP in every 42 bp. It was possible to clearly discriminate SNPs that are polymorphic between different plant varieties from SNPs differentiating the two ancestral genomes. To validate the identified SNPs for their use in genetic analysis, we have developed Illumina GoldenGate assays for some of the identified SNPs. Through the analysis of a number of oilseed rape varieties and mapping populations with GoldenGate assays, we were able to identify a number of different segregation patterns in allotetraploid oilseed rape. The majority of the identified SNP markers can be readily used for genetic mapping, showing that amplicon sequencing and Illumina GoldenGate assays can be used to reliably identify SNP markers in tetraploid oilseed rape and to convert them into successful SNP assays that can be used for genetic analysis.     

Biased gene conversion is not occurring among rDNA repeats in the Brassica triangle.

Hybridization is a common phenomenon that results in complex genomes. How ancestral genomes interact in hybrids has long been of great interest. Recombination among ancestral genomes may increase or decrease genetic variation. This study examines rDNA from members of the Brassica triangle for evidence of gene conversion across ancestral genomes. Gene conversion is a powerful force in the evolution of multigene families. It has previously been shown that biased gene conversion can act to homogenize rDNA repeats within hybrid genomes. Here, we find no evidence for biased gene conversion or unequal crossing over across ancestral genomes in allotetraploid Brassica species. We suggest that, while basic genomic processes are shared by all organisms, the relative frequency of these processes and their evolutionary importance may differ among lineages. Key words : Brassica, rDNA, gene conversion, allotetraploids.     

Proteome analysis of diploid,tetraploid and hexaploid wheat: towards understanding genome interaction in protein expression

Hexaploid wheat (Triticum aestivum L.) is derived from a complex hybridization procedure involving three diploid species carrying the A, B and D genomes. The proteome patterns of diploid, tetraploid and hexaploid wheat were analyzed to explore the genome interaction in protein expression. At least two species from each of the diploid and tetraploid were used to compare their proteome maps with a hexaploid wheat cv. Chinese Spring. The ancestral cultivars were selected based on their history of closeness with the cultivated wheat. Proteins were extracted from seed flour and separated by two-dimensional electrophoresis (2-DE) with isoelectric focusing of pH range from 4-10. 2-DE maps of cultivated and ancestral species were analyzed by computer assisted image analyzer. The region of high molecular weight glutenin subunits of hexaploid wheat showed similarity with those of the diploid donors, BB and DD genomes. The omega gliadin, which is controlled by B genome in common wheat, was assumed to have evolved as a result of interaction between AA and BB genomes. The low molecular weight glutenins and alpha and beta gliadin regions were contributed by the three genomes. This result suggests that the function of donor genomes particularly in the expression of proteins in hexaploid wheat is not totally independent; rather it is the product of interactions among the diploid genomes in the hexaploid nuclear constitutions. The expression of nonstorage proteins was affected substantially due to the removal of the D genome from hexaploid constitution. Location of the structural gene controlling one of the alpha amylase inhibitor proteins in the nonstorage protein region was identified in the short arm of chromosome 3D.     

Ancestral state reconstructions for genomes

The recent expansion of phylogenetic analysis from the traditional field of molecular evolution, analyzing histories of genes, to the nascent field of "genomic evolution", analyzing histories of entire genomes, enables the construction of trees based on genome information, the quantification of the key processes that shape genome content and, ultimately, plausible parsimony reconstructions of ancestral genomes. Thus, when genomes are considered as phylogenetic characters, it is possible to reconstruct not only the history of species but also the ancestral states in terms of genome structure or function. In the future, we might be able to accurately reconstruct--or retrodict--a chain of events that led to the emergence of a specific genome sequence and, ultimately, to synthesize ancestral genomes at will, creating a "Jurassic database" of genomes.     

Inferring Population Histories for Ancient Genomes Using Genome-Wide Genealogies

Ancient genomes anchor genealogies in directly observed historical genetic variation and contextualize ancestral lineages with archaeological insights into their geography and cultural associations. However, the majority of ancient genomes are of lower coverage and cannot be directly built into genealogies. Here, we present a fast and scalable method, Colate, the first approach for inferring ancestral relationships through time between low-coverage genomes without requiring phasing or imputation. Our approach leverages sharing patterns of mutations dated using a genealogy to infer coalescence rates. For deeply sequenced ancient genomes, we additionally introduce an extension of the Relate algorithm for joint inference of genealogies incorporating such genomes. Application to 278 present-day and 430 ancient DNA samples of >0.5x mean coverage allows us to identify dynamic population structure and directional gene flow between early farmer and European hunter-gatherer groups. We further show that the previously reported, but still unexplained, increase in the TCC/TTC mutation rate, which is strongest in West Eurasia today, was already present at similar strength and widespread in the Late Glacial Period ~10k−15k years ago, but is not observed in samples >30k years old. It is strongest in Neolithic farmers, and highly correlated with recent coalescence rates between other genomes and a 10,000-year-old Anatolian hunter-gatherer. This suggests gene-flow among ancient peoples postdating the last glacial maximum as widespread and localizes the driver of this mutational signal in both time and geography in that region. Our approach should be widely applicable in future for addressing other evolutionary questions, and in other species.     

Properties of the new D3-like Pseudomonas aeruginosa bacteriophage phiPMG1: Genome structure and prospects for the use in phage therapy

Results of studying the novel virulent phage phiPMG1 active on Pseudomonas aeruginosa are presented. It is shown that phiPMG1 exhibits significant homology and the similarity in the overall structure with the genome of a temperate phage converts D3. Phage phiPMG1 differs from D3 in that it fails to stably lysogenize bacteria and can grow on strains carrying plasmids that cause growth inhibition of phage D3 and some other phages. This significantly diminishes the probability of horizontal gene transfer with phage phiPMG1 and suggests the possible employment of this phage in phage therapy. A comparison of phages phiPMG1 and D3 structures of genomes in demonstrated not only high homology of 65 genes, but also the presence of 16 genes in the phiPMG1 genome that were not included in the in NCBI database. Apparently, the evolution of genomes in phages of this species is mostly associated with migrations into other species of bacteria, and recombinations with phages of other species (for example, F116). A detailed analysis of structure of one region genomes, which significant nonhomology for the three D3-like phages (D3, phiPMG1 and PAJU2), revealed that the phiPMG1 genome possible closest to a hypothetical genome of ancestral phage of this species.     

Assessing the level of collinearity between Arabidopsis thaliana and Brassica napus for A. thaliana chromosome 5.

This study describes a comprehensive comparison of chromosome 5 of the model crucifer Arabidopsis with the genome of its amphidiploid crop relative Brassica napus and introduces the use of in silico sequence homology to identify conserved loci between the two species. A region of chromosome 5, spanning 8 Mb, was found in six highly conserved copies in the B. napus genome. A single inversion appeared to be the predominant rearrangement that had separated the two lineages leading to the formation of Arabidopsis chromosome 5 and its homologues in B. napus. The observed results could be explained by the fusion of three ancestral genomes with strong similarities to modern-day Arabidopsis to generate the constituent diploid genomes of B. napus. This supports the hypothesis that the diploid Brassica genomes evolved from a common hexaploid ancestor. Alignment of the genetic linkage map of B. napus with the genomic sequence of Arabidopsis indicated that for specific regions a genetic distance of 1 cM in B. napus was equivalent to 285 Kb of Arabidopsis DNA sequence. This analysis strongly supports the application of Arabidopsis as a tool in marker development, map-based gene cloning, and candidate gene identification for the larger genomes of Brassica crop species.     

The complete mitochondrial genomes of two ghost moths, Thitarodes renzhiensis and Thitarodes yunnanensis: the ancestral gene arrangement in Lepidoptera

ABSTRACT: BACKGROUND: Lepidoptera encompasses more than 160,000 described species that have been classified into 45-48 superfamilies. The previously determined Lepidoptera mitochondrial genomes (mitogenomes) are limited to six superfamilies of the most derived lepidopteran lineage Ditrysia. Compared with the ancestral insect gene order, these mitogenomes all contain a tRNA rearrangement. To gain new insights into Lepidoptera mitogenome evolution, we sequenced the mitogenomes of two ghost moths that belong to primitive lepidopteran lineages and conducted a comparative mitogenomic analysis across Lepidoptera. RESULTS: The mitogenomes of Thitarodes renzhiensis and T. yunnanensis are 16,173 bp and 15,814 bp long with an A+T content of 81.28% and 82.33%, respectively. Different tandem repeats in the A+T-rich region mainly account for the size difference between the two mitogenomes. Both mitogenomes include 13 protein-coding genes, 22 transfer RNA genes, and 2 ribosomal RNA genes. The 1,584-bp sequence from rrnS to nad2 was also determined for Thitarodes sp.QL, which has no repetitive sequence in the A+T-rich region. All three Thitarodes species possess the ancestral gene order with trnI-trnQ-trnM located between the A+T-rich region and nad2, which is different from the gene order trnM-trnI-trnQ in all previously sequenced Lepidoptera species. The formerly identified conserved elements of Lepidoptera mitogenomes (i.e. the motif 'ATAGA' and poly-T stretch in the A+T-rich region and the long intergenic spacer upstream of nad2) are absent in the Thitarodes mitogenomes. The phylogenetic analysis supports that Hepialoidea, represented by T. renzhiensis and T. yunnanensis, occupies a basal position in the currently sampled seven superfamilies. The relationships of the other six superfamilies are (((((Bombycoidea + Geometroidea) + Noctuoidea) + Pyraloidea) + Papilionoidea) + Tortricoidea). CONCLUSION: The mitogenomes of T. renzhiensis and T. yunnanensis exhibit unusual features compared with the previously determined Lepidoptera mitogenomes. Their ancestral gene order indicates that the tRNA rearrangement event occurred after Lepidoptera diverged from other holometabolous insect orders. Phylogenetic analysis based on mitogenome sequences is a power tool for addressing phylogenetic relationships among major Lepidoptera superfamilies. Characterization of the two ghost moth mitogenomes has enriched our knowledge of Lepidoptera mitogenomes and contributed to our understanding of the mechanisms underlying mitogenome evolution, especially gene rearrangements.     

Sequence and phylogenetic analysis of the complete mitochondrial genome of the flour beetle Tribolium castanaeum

We describe the first complete mitochondrial genome sequence from a representative of the insect order Coleoptera, the flour beetle Tribolium castaneum. The 15,881 bp long Tribolium mitochondrial genome encodes 13 putative proteins, two ribosomal RNAs and 22 tRNAs canonical for animal mitochondrial genomes. Their arrangement is identical to that in Drosophila melanogaster, which is considered ancestral for insects and crustaceans (Boore et al., 1998; Hwang, et al., 2001a). Nucleotide composition, amino acid composition, and codon usage fall within the range of values observed in other insect mitochondrial genomes. Most notable features are the use of TCT as tRNA(Ser(AGN)) anticodon instead of GCT, which is used in most other arthropod species, and the relative scarcity of special sequence motifs in the 1431 bp long control region. Phylogenetic analysis confirmed resolving power in the conserved regions of the mitochondrial proteome regarding diversification events, which predate the emergence of pterygote insects, while little resolution was obtained at the level of basal perygote diversification. The partition of faster evolving amino acid sites harbored strong support for joining Lepidoptera with Diptera, which is consistent with a monophyletic Mecopterida.     

Deciphering the evolutionary interplay between subgenomes following polyploidy: A paleogenomics approach in grasses

How did plant species emerge from their most recent common ancestors (MRCAs) 250 million years ago? Modern plant genomes help to address such key questions in unveiling precise species genealogies. The field of paleogenomics is undergoing a paradigm shift for investigating species evolution from the study of ancestral genomes from extinct species to deciphering the evolutionary forces (in terms of duplication, fusion, fission, deletion, and translocation) that drove present‐day plant diversity (in terms of chromosome/gene number and genome size). In this review, inferred ancestral karyotype genomes are shown to be powerful tools to (1) unravel the past history of extant species by recovering the variations of ancestral genomic compartments and (2) accelerate translational research by facilitating the transfer of genomic information from model systems to species of agronomic interest.     

A hotspot of gene order rearrangement by tandem duplication and random loss in the vertebrate mitochondrial genome

Most reported examples of change in vertebrate mitochondrial (mt) gene order could be explained by a tandem duplication followed by random loss of redundant genes (tandem duplication-random loss [TDRL] model). Under this model of evolution, independent loss of genes arising from a single duplication in an ancestral species are predicted, and remnant pseudogenes expected, intermediate states that may remain in rearranged genomes. However, evidence for this is rare and largely scattered across vertebrate lineages. Here, we report new derived mt gene orders in the vertebrate "WANCY" region of four closely related caecilian amphibians. The novel arrangements found in this genomic region (one of them is convergent with the derived arrangement of marsupials), presence of pseudogenes, and positions of intergenic spacers fully satisfy predictions from the TDRL model. Our results, together with comparative data for the available vertebrate complete mt genomes, provide further evidence that the WANCY genomic region is a hotspot for gene order rearrangements and support the view that TDRL is the dominant mechanism of gene order rearrangement in vertebrate mt genomes. Convergent gene rearrangements are not unlikely in hotspots of gene order rearrangement by TDRL.     

Genome size evolution of the extant lycophytes and ferns

Ferns and lycophytes have remarkably large genomes. However, little is known about how their genome size evolved in fern lineages. To explore the origins and evolution of chromosome numbers and genome size in ferns, we used flow cytometry to measure the genomes of 240 species (255 samples) of extant ferns and lycophytes comprising 27 families and 72 genera, of which 228 species (242 samples) represent new reports. We analyzed correlations among genome size, spore size, chromosomal features, phylogeny, and habitat type preference within a phylogenetic framework. We also applied ANOVA and multinomial logistic regression analysis to preference of habitat type and genome size. Using the phylogeny, we conducted ancestral character reconstruction for habitat types and tested whether genome size changes simultaneously with shifts in habitat preference. We found that 2C values had weak phylogenetic signal, whereas the base number of chromosomes (x) had a strong phylogenetic signal. Furthermore, our analyses revealed a positive correlation between genome size and chromosome traits, indicating that the base number of chromosomes (x), chromosome size, and polyploidization may be primary contributors to genome expansion in ferns and lycophytes. Genome sizes in different habitat types varied significantly and were significantly correlated with habitat types; specifically, multinomial logistic regression indicated that species with larger 2C values were more likely to be epiphytes. Terrestrial habitat is inferred to be ancestral for both extant ferns and lycophytes, whereas transitions to other habitat types occurred as the major clades emerged. Shifts in habitat types appear be followed by periods of genomic stability. Based on these results, we inferred that habitat type changes and multiple whole-genome duplications have contributed to the formation of large genomes of ferns and their allies during their evolutionary history.     

Ancient linkage of a POU class 6 and an anterior Hox-like gene in cnidaria: implications for the evolution of homeobox genes

Linkage analyses in metazoan genomes suggest two ancestral arrays for the majority of homeobox genes. The related homeobox genes and chromosomal regions that are dispersed in extant species derived possibly from only two single common ancestor regions. One proposed ancestral array, designated as ANTP mega-array, contains most of the ANTP class homeobox genes; the second, named the contraHox super-paralogon, would consist of the classes PRD, POU, LIM, CUT, prospero, TALE and SIX. Here, we report the tight linkage of a POU class 6 gene to an anterior Hox-like gene in the hydrozoan Eleutheria dichotoma and discuss its possible significance for the evolution of homeobox genes. POU class 6 genes also seem to be ancestrally linked to the HoxC and A clusters in vertebrates, despite POU homeobox genes belonging to the contraHox paralogon. Hence, the much tighter linkage of a POU class 6 gene to an anterior Hox-like gene in a cnidarian is possibly the evolutionary echo of an ancestral genomic region from which most metazoan homeobox classes emerged.     

Comparative Mitogenomics of the Assassin Bug Genus Peirates (Hemiptera: Reduviidae: Peiratinae) Reveal Conserved Mitochondrial Genome Organization of P. atromaculatus,P. fulvescens and P. turpis

In this study, we sequenced four new mitochondrial genomes and presented comparative mitogenomic analyses of five species in the genus Peirates (Hemiptera: Reduviidae). Mitochondrial genomes of these five assassin bugs had a typical set of 37 genes and retained the ancestral gene arrangement of insects. The A+T content, AT- and GC-skews were similar to the common base composition biases of insect mtDNA. Genomic size ranges from 15,702 bp to 16,314 bp and most of the size variation was due to length and copy number of the repeat unit in the putative control region. All of the control region sequences included large tandem repeats present in two or more copies. Our result revealed similarity in mitochondrial genomes of P. atromaculatus, P. fulvescens and P. turpis, as well as the highly conserved genomic-level characteristics of these three species, e.g., the same start and stop codons of protein-coding genes, conserved secondary structure of tRNAs, identical location and length of non-coding and overlapping regions, and conservation of structural elements and tandem repeat unit in control region. Phylogenetic analyses also supported a close relationship between P. atromaculatus, P. fulvescens and P. turpis, which might be recently diverged species. The present study indicates that mitochondrial genome has important implications on phylogenetics, population genetics and speciation in the genus Peirates.     

Reconstruction of ancestral karyotype illuminates chromosome evolution in the genus Cucumis

Karyotype dynamics driven by complex chromosome rearrangements constitute a fundamental issue in evolutionary genetics. The evolutionary events underlying karyotype diversity within plant genera, however, have rarely been reconstructed from a computed ancestral progenitor. Here, we developed a method to rapidly and accurately represent extant karyotypes with the genus, Cucumis, using highly customizable comparative oligo-painting (COP) allowing visualization of fine-scale genome structures of eight Cucumis species from both African-origin and Asian-origin clades. Based on COP data, an evolutionary framework containing a genus-level ancestral karyotype was reconstructed, allowing elucidation of the evolutionary events that account for the origin of these diverse genomes within Cucumis. Our results characterize the cryptic rearrangement hotspots on ancestral chromosomes, and demonstrate that the ancestral Cucumis karyotype (n = 12) evolved to extant Cucumis genomes by hybridizations and frequent lineage- and species-specific genome reshuffling. Relative to the African species, the Asian species, including melon (Cucumis melo, n = 12), Cucumis hystrix (n = 12) and cucumber (Cucumis sativus, n = 7), had highly shuffled genomes caused by large-scale inversions, centromere repositioning and chromothripsis-like rearrangement. The deduced reconstructed ancestral karyotype for the genus allowed us to propose evolutionary trajectories and specific events underlying the origin of these Cucumis species. Our findings highlight that the partitioned evolutionary plasticity of Cucumis karyotype is primarily located in the centromere-proximal regions marked by rearrangement hotspots, which can potentially serve as a reservoir for chromosome evolution due to their fragility.     

Evolutionary convergence on highly-conserved 3' intron structures in intron-poor eukaryotes and insights into the ancestral eukaryotic genome

The presence of spliceosomal introns in eukaryotes raises a range of questions about genomic evolution. Along with the fundamental mysteries of introns' initial proliferation and persistence, the evolutionary forces acting on intron sequences remain largely mysterious. Intron number varies across species from a few introns per genome to several introns per gene, and the elements of intron sequences directly implicated in splicing vary from degenerate to strict consensus motifs. We report a 50-species comparative genomic study of intron sequences across most eukaryotic groups. We find two broad and striking patterns. First, we find that some highly intron-poor lineages have undergone evolutionary convergence to strong 3' consensus intron structures. This finding holds for both branch point sequence and distance between the branch point and the 3' splice site. Interestingly, this difference appears to exist within the genomes of green alga of the genus Ostreococcus, which exhibit highly constrained intron sequences through most of the intron-poor genome, but not in one much more intron-dense genomic region. Second, we find evidence that ancestral genomes contained highly variable branch point sequences, similar to more complex modern intron-rich eukaryotic lineages. In addition, ancestral structures are likely to have included polyT tails similar to those in metazoans and plants, which we found in a variety of protist lineages. Intriguingly, intron structure evolution appears to be quite different across lineages experiencing different types of genome reduction: whereas lineages with very few introns tend towards highly regular intronic sequences, lineages with very short introns tend towards highly degenerate sequences. Together, these results attest to the complex nature of ancestral eukaryotic splicing, the qualitatively different evolutionary forces acting on intron structures across modern lineages, and the impressive evolutionary malleability of eukaryotic gene structures.     

Rates of gene rearrangement and nucleotide substitution are correlated in the mitochondrial genomes of insects

A number of studies indicated that lineages of animals with high rates of mitochondrial (mt) gene rearrangement might have high rates of mt nucleotide substitution. We chose the hemipteroid assemblage and the Insecta to test the idea that rates of mt gene rearrangement and mt nucleotide substitution are correlated. For this purpose, we sequenced the mt genome of a lepidopsocid from the Psocoptera, the only order of hemipteroid insects for which an entire mtDNA sequence is not available. The mt genome of this lepidopsocid is circular, 16,924 bp long, and contains 37 genes and a putative control region; seven tRNA genes and a protein-coding gene in this genome have changed positions relative to the ancestral arrangement of mt genes of insects. We then compared the relative rates of nucleotide substitution among species from each of the four orders of hemipteroid insects and among the 20 insects whose mt genomes have been sequenced entirely. All comparisons among the hemipteroid insects showed that species with higher rates of gene rearrangement also had significantly higher rates of nucleotide substitution statistically than did species with lower rates of gene rearrangement. In comparisons among the 20 insects, where the mt genomes of the two species differed by more than five breakpoints, the more rearranged species always had a significantly higher rate of nucleotide substitution than the less rearranged species. However, in comparisons where the mt genomes of two species differed by five or less breakpoints, the more rearranged species did not always have a significantly higher rate of nucleotide substitution than the less rearranged species. We tested the statistical significance of the correlation between the rates of mt gene rearrangement and mt nucleotide substitution with nine pairs of insects that were phylogenetically independent from one another. We found that the correlation was positive and statistically significant (R2 = 0.73, P = 0.01; Rs = 0.67, P < 0.05). We propose that increased rates of nucleotide substitution may lead to increased rates of gene rearrangement in the mt genomes of insects.     

Rapid fixation of a distinctive sequence motif in the 3' noncoding region of the clade of West Nile virus invading North America

Phylogenetic analysis of complete genomes of West Nile virus (WNV) by a variety of methods supported the hypothesis that North American isolates of WNV constitute a monophyletic group, together with an isolate from Israel and one from Hungary. We used ancestral sequence reconstruction in order to obtain evidence for evolutionary changes that might be correlated with increased virulence in this clade (designated the N.A. clade). There was one amino acid change (I-->T at residue 356 of the NS3 protein) that occurred in the ancestor of the N.A. clade and remained conserved in all N.A. clade genomes analyzed. There were four changes in the upstream portion of the 3' noncoding region (the AT-enriched region) that occurred in the ancestor of the N.A. clade and remained conserved in all N.A. clade genomes analyzed, changes predicted to alter RNA secondary structure. The AT-enriched region showed a higher rate of substitution in the branch ancestral to the N.A. clade, relative to polymorphism, than did the remainder of the noncoding regions, synonymous sites in coding regions, or nonsynonymous sites in coding regions. The high rate of occurrence of fixed nucleotide substitutions in this region suggests that positive Darwinian selection may have acted on this portion of the 3'NCR and that these fixed changes, possibly in concert with the amino acid change in NS3, may underlie phenotypic effects associated with increased virulence in North American WNV.