Search

Background Whole genome sequencing (WGS) is increasingly used in diagnostics and surveillance of infectious diseases. A major application for WGS is to use the data for identifying outbreak clusters, and there is therefore a need for methods that can accurately and efficiently infer phylogenies from sequencing reads. In the present study we describe a new dataset that we have created for the purpose of benchmarking such WGS-based methods for epidemiological data, and also present an analysis where we use the data to compare the performance of some current methods. Results Our aim was to create a benchmark data set that mimics sequencing data of the sort that might be collected during an outbreak of an infectious disease. This was achieved by letting an E. coli hypermutator strain grow in the lab for 8 consecutive days, each day splitting the culture in two while also collecting samples for sequencing. The result is a data set consisting of 101 whole genome sequences with known phylogenetic relationship. Among the sequenced samples 51 correspond to internal nodes in the phylogeny because they are ancestral, while the remaining 50 correspond to leaves.We also used the newly created data set to compare three different online available methods that infer phylogenies from whole-genome sequencing reads: NDtree, CSI Phylogeny and REALPHY. One complication when comparing the output of these methods with the known phylogeny is that phylogenetic methods typically build trees where all observed sequences are placed as leafs, even though some of them are in fact ancestral. We therefore devised a method for post processing the inferred trees by collapsing short branches (thus relocating some leafs to internal nodes), and also present two new measures of tree similarity that takes into account the identity of both internal and leaf nodes. Conclusions Based on this analysis we find that, among the investigated methods, CSI Phylogeny had the best performance, co

Background: Reduced intake of n-3 long-chain polyunsaturated fatty acids (LCPUFAs) may be a contributing factor to the increasing prevalence of wheezing disorders. We assessed the effect of supplementation with n-3 LCPUFAs in pregnant women on the risk of persistent wheeze and asthma in their offspring. Methods: We randomly assigned 736 pregnant women at 24 weeks of gestation to receive 2.4 g of n-3 LCPUFA (fish oil) or placebo (olive oil) per day. Their children formed the Copenhagen Prospective Studies on Asthma in Childhood2010 (COPSAC2010) cohort and were followed prospectively with extensive clinical phenotyping. Neither the investigators nor the participants were aware of group assignments during follow-up for the first 3 years of the children's lives, after which there was a 2-year follow-up period during which only the investigators were unaware of group assignments. The primary end point was persistent wheeze or asthma, and the secondary end points included lower respiratory tract infections, asthma exacerbations, eczema, and allergic sensitization. Results: A total of 695 children were included in the trial, and 95.5% completed the 3-year, double-blind follow-up period. The risk of persistent wheeze or asthma in the treatment group was 16.9%, versus 23.7% in the control group (hazard ratio, 0.69; 95% confidence interval [CI], 0.49 to 0.97; P=0.035), corresponding to a relative reduction of 30.7%. Prespecified subgroup analyses suggested that the effect was strongest in the children of women whose blood levels of eicosapentaenoic acid and docosahexaenoic acid were in the lowest third of the trial population at randomization: 17.5% versus 34.1% (hazard ratio, 0.46; 95% CI, 0.25 to 0.83; P=0.011). Analyses of secondary end points showed that supplementation with n-3 LCPUFA was associated with a reduced risk of infections of the lower respiratory tract (31.7% vs. 39.1%; hazard ratio, 0.75; 95% CI, 0.58 to 0.98; P=0.033), but there was no statistically sig

Amplicon pyrosequencing targets a known genetic region and thus inherently produces reads highly anticipated to have certain features, such as conserved nucleotide sequence, and in the case of protein coding DNA, an open reading frame. Pyrosequencing errors, consisting mainly of nucleotide insertions and deletions, are on the other hand likely to disrupt open reading frames. Such an inverse relationship between errors and expectation based on prior knowledge can be used advantageously to guide the process known as basecalling, i.e. the inference of nucleotide sequence from raw sequencing data. The new basecalling method described here, named Multipass, implements a probabilistic framework for working with the raw flowgrams obtained by pyrosequencing. For each sequence variant Multipass calculates the likelihood and nucleotide sequence of several most likely sequences given the flowgram data. This probabilistic approach enables integration of basecalling into a larger model where other parameters can be incorporated, such as the likelihood for observing a full-length open reading frame at the targeted region. We apply the method to 454 amplicon pyrosequencing data obtained from a malaria virulence gene family, where Multipass generates 20 % more error-free sequences than current state of the art methods, and provides sequence characteristics that allow generation of a set of high confidence error-free sequences. This novel method can be used to increase accuracy of existing and future amplicon sequencing data, particularly where extensive prior knowledge is available about the obtained sequences, for example in analysis of the immunoglobulin VDJ region where Multipass can be combined with a model for the known recombining germline genes. Multipass is available for Roche 454 data at http://www.cbs.dtu.dk/services/MultiPass-1.0 , and the concept can potentially be implemented for other sequencing technologies as well.

Aleutian Mink Disease Virus (AMDV) is a frequently encountered pathogen associated with commercial mink breeding. AMDV infection leads to increased mortality and compromised animal health and welfare. Currently little is known about the molecular evolution of the virus, and the few existing studies have focused on limited regions of the viral genome. This paper describes a robust, reliable, and fast protocol for amplification of the full AMDV genome using long-range PCR. The method was used to generate next generation sequencing data for the non-virulent cell-culture adapted AMDV-G strain as well as for the virulent AMDV-Utah strain. Comparisons at nucleotide- and amino acid level showed that, in agreement with existing literature, the highest variability between the two virus strains was found in the left open reading frame, which encodes the non-structural (NS1–3) genes. This paper also reports a number of differences that potentially can be linked to virulence and host range. To the authors’ knowledge, this is the first study to apply next generation sequencing on the entire AMDV genome. The results from the study will facilitate the development of new diagnostic tools and can form the basis for more detailed molecular epidemiological analyses of the virus.

RNA viruses have the highest known mutation rates. Consequently it is likely that a high proportion of individual RNA virus genomes, isolated from an infected host, will contain lethal mutations and be non-functional. This is problematic if the aim is to clone and investigate high-fitness, functional cDNAs and may also pose problems for sequence-based analysis of viral evolution. To address these challenges we have performed a study of the evolution of classical swine fever virus (CSFV) using deep sequencing and analysis of 84 full-length cDNA clones, each representing individual genomes from a moderately virulent isolate. In addition to here being used as a model for RNA viruses generally, CSFV has high socioeconomic importance and remains a threat to animal welfare and pig production. We find that the majority of the investigated genomes are non-functional and only 12% produced infectious RNA transcripts. Full length sequencing of cDNA clones and deep sequencing of the parental population identified substitutions important for the observed phenotypes. The investigated cDNA clones were furthermore used as the basis for inferring the sequence of functional viruses. Since each unique clone must necessarily be the descendant of a functional ancestor, we hypothesized that it should be possible to produce functional clones by reconstructing ancestral sequences. To test this we used phylogenetic methods to infer two ancestral sequences, which were then reconstructed as cDNA clones. Viruses rescued from the reconstructed cDNAs were tested in cell culture and pigs. Both reconstructed ancestral genomes proved functional, and displayed distinct phenotypes in vitro and in vivo. We suggest that reconstruction of ancestral viruses is a useful tool for experimental and computational investigations of virulence and viral evolution. Importantly, ancestral reconstruction can be done even on the basis of a set of sequences that all correspond to non-functional variants.

Many bacterial, viral and parasitic pathogens undergo antigenic variation to counter host immune defense mechanisms. In Plasmodium falciparum, the most lethal of human malaria parasites, switching of var gene expression results in alternating expression of the adhesion proteins of the Plasmodium falciparum-erythrocyte membrane protein 1 class on the infected erythrocyte surface. Recombination clearly generates var diversity, but the nature and control of the genetic exchanges involved remain unclear. By experimental and bioinformatic identification of recombination events and genome-wide recombination hotspots in var genes, we show that during the parasite's sexual stages, ectopic recombination between isogenous var paralogs occurs near low folding free energy DNA 50-mers and that these sequences are heavily concentrated at the boundaries of regions encoding individual Plasmodium falciparum-erythrocyte membrane protein 1 structural domains. The recombinogenic potential of these 50-mers is not parasite-specific because these sequences also induce recombination when transferred to the yeast Saccharomyces cerevisiae. Genetic cross data suggest that DNA secondary structures (DSS) act as inducers of recombination during DNA replication in P. falciparum sexual stages, and that these DSS-regulated genetic exchanges generate functional and diverse P. falciparum adhesion antigens. DSS-induced recombination may represent a common mechanism for optimizing the evolvability of virulence gene families in pathogens.

Classical swine fever virus (CSFV) may be highly virulent in pigs with a mortality rate close to 100%. The CSFV “Koslov strain” is known to be one of the most virulent CSFV, but so far a functional cloned cDNA of this strain has not been described. We suggest that this may be due to the error-prone nature of the RNA-dependent RNA polymerase resulting in the majority of circulating forms being non-functional. However, since any infectious virus particle should necessarily be the offspring of a functional virus, we hypothesized that it should be possible to synthesize a highly virulent form by reconstructing ancestral sequences. To test this hypothesis, we inferred sequences that correspond to ancestral nodes in a phylogenetic tree built from full-length nucleotide sequences of non-functional Koslov cDNAs and then proceeded to test the reconstructions. Specifically, we altered a non-functional cDNA by site directed mutagenesis, removing non-synonymous mutations step by step. In vitro testing of modified constructs did indeed lead to fully functional viruses with similar growth kinetics as the wild-type strain. Moreover, viruses rescued from the construct had the ancestral amino acid sequence and, when tested in pigs, were at least as virulent as the Koslov strain. The ancestral reconstruction therefore proved to give rise to a functional cDNA of the highly virulent Koslov strain. In vivo studies confirmed our methods and enabled us to identify nucleotide positions within the viral genome important for virulence.

Next Generation Sequencing (NGS) is becoming more adopted into viral research and will be the preferred technology in the years to come. We have recently sequenced several strains of Classical Swine Fever Virus (CSFV) by NGS on both Genome Sequencer FLX (GS FLX) and Iontorrent PGM platforms. In this study, we analyzed NGS data of virus rescued from a CSFV C-strain vaccine strain cDNA clone. The virus analyzed was obtained from a 4th and a 12th passage of rescued virus in SFT cell culture, which had shown a difference in growth kinetics between the passages, and NGS analysis was chosen in order to look for molecular differences. Identical RT-PCR products were run on both GS FLX and an Iontorrent PGM platform for comparison. The NGS data was compared by quality and the percentage mapped reads. Results showed good quality of reads for both platforms and a close to 100% of the reads mapped to the consensus sequence. Additionally, we got an average sequence depth for the genome of 4000 for the Iontorrent PGM and 400 for the FLX platform making the mapping suitable for single nucleotide variant (SNV) detection. The analysis revealed a single non-silent SNV A10665G leading to the amino acid change D3431G in the RNAdependent RNA polymerase NS5B. This SNV was present at 100% frequency in the 12th passage and only at 55% in the 4th passage, which could explain the difference in growth kinetics between the passages.

The var gene encoded hyper-variable Plasmodium falciparum erythrocyte membrane protein 1 (PfEMP1) family mediates cytoadhesion of infected erythrocytes to human endothelium. Antibodies blocking cytoadhesion are important mediators of malaria immunity acquired by endemic populations. The development of a PfEMP1 based vaccine mimicking natural acquired immunity depends on a thorough understanding of the evolved PfEMP1 diversity, balancing antigenic variation against conserved receptor binding affinities. This study redefines and reclassifies the domains of PfEMP1 from seven genomes. Analysis of domains in 399 different PfEMP1 sequences allowed identification of several novel domain classes, and a high degree of PfEMP1 domain compositional order, including conserved domain cassettes not always associated with the established group A–E division of PfEMP1. A novel iterative homology block (HB) detection method was applied, allowing identification of 628 conserved minimal PfEMP1 building blocks, describing on average 83% of a PfEMP1 sequence. Using the HBs, similarities between domain classes were determined, and Duffy binding-like (DBL) domain subclasses were found in many cases to be hybrids of major domain classes. Related to this, a recombination hotspot was uncovered between DBL subdomains S2 and S3. The VarDom server is introduced, from which information on domain classes and homology blocks can be retrieved, and new sequences can be classified. Several conserved sequence elements were found, including: (1) residues conserved in all DBL domains predicted to interact and hold together the three DBL subdomains, (2) potential integrin binding sites in DBLα domains, (3) an acylation motif conserved in group A var genes suggesting N-terminal N-myristoylation, (4) PfEMP1 inter-domain regions proposed to be elastic disordered structures, and (5) several conserved predicted phosphorylation sites. Ideally, this comprehensive categorization of PfEMP1 will provide a platform f

BACKGROUND: The presence of gaps in an alignment of nucleotide or protein sequences is often an inconvenience for bioinformatical studies. In phylogenetic and other analyses, for instance, gapped columns are often discarded entirely from the alignment. RESULTS: MaxAlign is a program that optimizes the alignment prior to such analyses. Specifically, it maximizes the number of nucleotide (or amino acid) symbols that are present in gap-free columns - the alignment area - by selecting the optimal subset of sequences to exclude from the alignment. MaxAlign can be used prior to phylogenetic and bioinformatical analyses as well as in other situations where this form of alignment improvement is useful. In this work we test MaxAlign's performance in these tasks and compare the accuracy of phylogenetic estimates including and excluding gapped columns from the analysis, with and without processing with MaxAlign. In this paper we also introduce a new simple measure of tree similarity, Normalized Symmetric Similarity (NSS) that we consider useful for comparing tree topologies. CONCLUSION: We demonstrate how MaxAlign is helpful in detecting misaligned or defective sequences without requiring manual inspection. We also show that it is not advisable to exclude gapped columns from phylogenetic analyses unless MaxAlign is used first. Finally, we find that the sequences removed by MaxAlign from an alignment tend to be those that would otherwise be associated with low phylogenetic accuracy, and that the presence of gaps in any given sequence does not seem to disturb the phylogenetic estimates of other sequences. The MaxAlign web-server is freely available online at http://www.cbs.dtu.dk/services/MaxAlign where supplementary information can also be found. The program is also freely available as a Perl stand-alone package.

The Mannheimia subclades belong to the same bacterial genus but have taken divergent paths toward their distinct lifestyles. M. haemolytica + M. glucosida are potential pathogens of the respiratory tract in the mammalian suborder Ruminantia, whereas M. ruminalis, the supposed sister group, lives as a commensal in the ovine rumen. We have tested the hypothesis that horizontal gene transfer of the leukotoxin operon has catalyzed pathogenic adaptation and speciation of M. haemolytica + M. glucosida, or other major subclades, by using a strategy that combines compositional and phylogenetic methods. We show that it has been vertically inherited from the last common ancestor of the diverging Mannheimia subclades, although several strains belonging to M. ruminalis have lost the operon. Our analyses support that divergence within M. ruminalis following colonization of the ovine rumen was very rapid and that functional decay of most of the leukotoxin operons occurred early when the adaptation to the rumen was fastest, suggesting that antagonistic pleiotropy was the main contributor to losses in the radiating lineages of M. ruminalis. To sum up, the scenario derived from these analyses reflects two aspects. On one hand, it opposes the hypothesis of horizontal gene transfer as a catalyst of pathogenic adaptation and speciation. On the other hand, it indicates that losses of the leukotoxin operons in the radiating lineages of M. ruminalis have catalyzed their adaptation to a commensal environment and reproductive isolation (speciation).

Eight atypical Mannheimia isolates were isolated from lesions in roe deer (Capreolus capreolus). Traditional classification based on morphologic and physiologic traits showed that they belong to a distinct biogroup (taxon) within genus Mannheimia. Extensive phenotypic characterization suggested that the isolates should be classified as M. granulomatis, although the presence of distinct traits justified their classification into a separate biogroup within this species. Phylogenetic analyses based on 16S rRNA sequences from two roe deer isolates and 41 other Mannheimia strains supported that the roe deer isolates form a monophyletic group within M. granulomatis. The lktA genotype was present in all roe deer isolates based on Southern blot analysis, whereas the corresponding beta-hemolytic phenotype was absent in one of these isolates.

BACKGROUND: Some amino acid residues functionally interact with each other. This interaction will result in an evolutionary co-variation between these residues - coevolution. Our goal is to find these coevolving residues. RESULTS: We present six new methods for detecting coevolving residues. Among other things, we suggest measures that are variants of Mutual Information, and measures that use a multidimensional representation of each residue in order to capture the physico-chemical similarities between amino acids. We created a benchmarking system, in silico, able to evaluate these methods through a wide range of realistic conditions. Finally, we use the combination of different methods as a way of improving performance. CONCLUSION: Our best method (Row and Column Weighed Mutual Information) has an estimated accuracy increase of 63% over Mutual Information. Furthermore, we show that the combination of different methods is efficient, and that the methods are quite sensitive to the different conditions tested.

In order to shed light on the nature of the persistent reservoir of human immunodeficiency virus type 1 (HIV-1), we investigated signs of recent evolution in the pool of proviral DNA in patients on successful HAART. Pro-viral DNA, corresponding to the C2-V3-C3 region of the HIV-1 env gene, was collected from PBMCs isolated from 57 patients. Both "consensus" (57 patients) and clonal (7 patients) sequences were obtained from five time points spanning a 24-month period. The main computational strategy was to use maximum likelihood to fit a set of alternative phylogenetic models to the clonal data, and then determine the support for models that imply evolution between time points. Model fit and model-selection uncertainty was assessed using the Akaike information criterion (AIC) and Akaike weights. The consensus sequence data was also analyzed using a range of phylogenetic techniques to determine whether there were temporal trends indicating ongoing replication and evolution. In summary, it was not possible to detect definitive signs of ongoing evolution in either the bulk-sequenced or the clonal data with the methods employed here, but our results could be consistent with localized expression of archival HIV genomes in some patients. Interestingly, stop-codons were present at the same two positions in several independent clones and across patients. Simulation studies indicated that this phenomenon could be explained as the result of parallel evolution and that some sites were inherently more likely to evolve into stop codons.

We performed a systematic, large-scale analysis of human protein complexes comprising gene products implicated in many different categories of human disease to create a phenome-interactome network. This was done by integrating quality-controlled interactions of human proteins with a validated, computationally derived phenotype similarity score, permitting identification of previously unknown complexes likely to be associated with disease. Using a phenomic ranking of protein complexes linked to human disease, we developed a Bayesian predictor that in 298 of 669 linkage intervals correctly ranks the known disease-causing protein as the top candidate, and in 870 intervals with no identified disease-causing gene, provides novel candidates implicated in disorders such as retinitis pigmentosa, epithelial ovarian cancer, inflammatory bowel disease, amyotrophic lateral sclerosis, Alzheimer disease, type 2 diabetes and coronary heart disease. Our publicly available draft of protein complexes associated with pathology comprises 506 complexes, which reveal functional relationships between disease-promoting genes that will inform future experimentation.

Udgivelsesdato: 2006-Nov, Pregnancy-associated malaria is a major health problem, which mainly affects primigravidae living in malaria endemic areas. The syndrome is precipitated by accumulation of infected erythrocytes in placental tissue through an interaction between chondroitin sulphate A on syncytiotrophoblasts and a parasite-encoded protein on the surface of infected erythrocytes, believed to be VAR2CSA. VAR2CSA is a polymorphic protein of approximately 3,000 amino acids forming six Duffy-binding-like (DBL) domains. For vaccine development it is important to define the antigenic targets for protective antibodies and to characterize the consequences of sequence variation. In this study, we used a combination of in silico tools, peptide arrays, and structural modeling to show that sequence variation mainly occurs in regions under strong diversifying selection, predicted to form flexible loops. These regions are the main targets of naturally acquired immunoglobulin gamma and accessible for antibodies reacting with native VAR2CSA on infected erythrocytes. Interestingly, surface reactive anti-VAR2CSA antibodies also target a conserved DBL3X region predicted to form an alpha-helix. Finally, we could identify DBL3X sequence motifs that were more likely to occur in parasites isolated from primi- and multigravidae, respectively. These findings strengthen the vaccine candidacy of VAR2CSA and will be important for choosing epitopes and variants of DBL3X to be included in a vaccine protecting women against pregnancy-associated malaria.

The leukotoxin protein of Mannheimia haemolytica belongs to the HlyA-like subfamily of cytotoxic RTX (repeats in toxin) proteins. To test the hypothesis that different lineages of genus Mannheimia gained the leukotoxin operon via horizontal gene transfer we used a strategy that combines compositional and phylogenetic methods: (i) ranking of genes according to their convergence to the average genome signature of M. haemolytica based on the relative 3:1 dinucleotide bias in 56 individual genes; (ii) estimation of ancestral character states on a phylogeny of 43 Mannheimia strains based on 16S rRNA sequences; (iii) phylogeny of 24 leukotoxin gene sequences and 16 homologous genes retrieved from SWISS-PROT by using PSI-BLAST. Our data show no evidence for horizontal gene transfer into this clade. We propose that vertical descent from the common ancestor of genus Mannheimia, with subsequent loss of these genes in many strains can explain the patchy distribution of leukotoxin in this clade. Finally, we argue that antagonistic pleiotropy was the main contributor to gene loss in the radiating M. ruminalis lines, which have adapted to different ecological niches in the rumen offering different conditions of growth.

Comparisons of evolutionary divergent promoter sequences have revealed that gross rearrangements of pre-existing blocks into new combinations play a central role during the evolution of regulatory regions. The noncoding region upstream of the leukotoxin operon was analyzed in six Mannheimia strains (Mhae1, Mhae2, Mglu, Mrum, Mgra, Mvar). The region from +30 to -208 showed sequence conservation in the Mannheimia strains, suggesting that this region represents an ancestral promoter block that has been maintained vertically in genus Mannheimia. Analyses of the 5¢ flanking regions showed the existence of three blocks that are restricted to specific lineages, of which the ancestral block is conserved in Mrum and Mvar. The transcriptional unit artJ has been shuffled into the upstream region in the last common ancestor of Mha1, Mha2 and Mglu. The corresponding non-coding DNA region ranging from -209 to -376 contains at least one positive regulatory element, suggesting early divergence within the promoter family and the emergence of a lineage with increased promoter strength. Tandem duplication of the pattern ACAAAAAACA has occurred in the ancestral Mhae1 promoter after the Mhae1 and Mhae2 ancestors diverged. The build-in A6-tract is inserted in phase with three ancestral tandem repeats containing one A6-tract each, which might have increased gene expression from the Mhae1 promoter. Phylogenetic analyses of the promoter and 16S rRNA sequences showed that the Mgra branch in the promoter tree was significantly longer than the corresponding branch in the 16S rRNA tree (p = 0.004). These data reflect a high rate of nucleotide substitution in the ancestral Mgra promoter and can be interpreted as pseudogenization. Thus, the mosaic structure of the leukotoxin promoter region has continually been shuffled into new combinations through rearrangements of pre-existing blocks. We suggest that some of these changes could

Major histocompatibility complex (MHC) proteins are encoded by extremely polymorphic genes and play a crucial role in immunity. However, not all genetically different MHC molecules are functionally different. Sette and Sidney (1999) have defined nine HLA class I supertypes and showed that with only nine main functional binding specificities it is possible to cover the binding properties of almost all known HLA class I molecules. Here we present a comprehensive study of the functional relationship between all HLA molecules with known specificities in a uniform and automated way. We have developed a novel method for clustering sequence motifs. We construct hidden Markov models for HLA class I molecules using a Gibbs sampling procedure and use the similarities among these to define clusters of specificities. These clusters are extensions of the previously suggested ones. We suggest splitting some of the alleles in the A1 supertype into a new A26 supertype, and some of the alleles in the B27 supertype into a new B39 supertype. Furthermore the B8 alleles may define their own supertype. We also use the published specificities for a number of HLA-DR types to define clusters with similar specificities. We report that the previously observed specificities of these class II molecules can be clustered into nine classes, which only partly correspond to the serological classification. We show that classification of HLA molecules may be done in a uniform and automated way. The definition of clusters allows for selection of representative HLA molecules that can cover the HLA specificity space better. This makes it possible to target most of the known HLA alleles with known specificities using only a few peptides, and may be used in construction of vaccines. Supplementary material is available at http://www.cbs.dtu.dk/researchgroups/immunology/supertypes.html.

Among 167 rotaviruis specimens collected from young children in a suburban area of Bissau, Guinea-Bissau, from 1996 to 1998, most identifiable strains belonged to the uncommon P[6], G2 type and approximately 50% remained incompletely typed. In the present study, 76 such strains were further characterized. Due to interprimer interaction during the standard multiplex PCR approach, modifications of this procedure were implemented. The modified analyses revealed a high frequency of G2, G8, and G9 genotypes, often combined with P[4] and/or P[6]. The Guinean G8 and G9 strains were 97 and 98%, respectively, identical to other African G8 and G9 strains. Multiple G and/or P types were identified at a high frequency (59%), including two previously undescribed mixed infections, P[4]P[6], G2G8 and P[4]P[6], G2G9. These mixed infections most likely represent naturally occurring reassortance of rotavirus strains. Detection of such strains among the previously incompletely typed strains indicates a potential underestimation of mixed infections, if only a standard multiplex PCR procedure is followed. Furthermore cross-priming of the G3 primer with the G8 primer binding site and silent mutations at the P[4] and P[6] primer binding sites were detected. These findings highlight the need for regular evaluation of the multiplex primer PCR method and typing primers. The high frequency of uncommon as well as reassortant rotavirus strains in countries where rotavirus is an important cause of child mortality underscores the need for extensive strain surveillance as a basis to develop appropriate rotavirus vaccine candidates.

The simple fact that proteins are built from 20 amino acids while DNA only contains four different bases, means that the 'signal-to-noise ratio' in protein sequence alignments is much better than in alignments of DNA. Besides this information-theoretical advantage, protein alignments also benefit from the information that is implicit in empirical substitution matrices such as BLOSUM-62. Taken together with the generally higher rate of synonymous mutations over non-synonymous ones, this means that the phylogenetic signal disappears much more rapidly from DNA sequences than from the encoded proteins. It is therefore preferable to align coding DNA at the amino acid level and it is for this purpose we have constructed the program RevTrans. RevTrans constructs a multiple DNA alignment by: (i) translating the DNA; (ii) aligning the resulting peptide sequences; and (iii) building a multiple DNA alignment by 'reverse translation' of the aligned protein sequences. In the resulting DNA alignment, gaps occur in groups of three corresponding to entire codons, and analogous codon positions are therefore always lined up. These features are useful when constructing multiple DNA alignments for phylogenetic analysis. RevTrans also accepts user-provided protein alignments for greater control of the alignment process. The RevTrans web server is freely available at http://www.cbs.dtu.dk/services/RevTrans/.