Hervé Bourhy

Ecology and evolution of rabies virus in Europe.

The evolution of rabies viruses of predominantly European origin was studied by comparing nucleotide sequences of the nucleoprotein and glycoprotein genes, and by typing isolates using RFLP. Phylogenetic analysis of the gene sequence data revealed a number of distinct groups, each associated with a particular geographical area. Such a pattern suggests that rabies virus has spread westwards and southwards across Europe during this century, but that physical barriers such as the Vistula river in Poland have enabled localized evolution. During this dispersal process, two species jumps took place - one into red foxes and another into raccoon dogs, although it is unclear whether virus strains are preferentially adapted to particular animal species or whether ecological forces explain the occurrence of the phylogenetic groups.

PCR technique as an alternative method for diagnosis and molecular epidemiology of rabies virus

We have investigated the PCR amplification technique of viral nucleic acids as an alternative protocol for diagnosis and epidemiological studies of rabies virus. A primer set mapping in the nucleoprotein cistron allowed a specific and sensitive amplification of infected brain material, fulfilling the diagnosis requirements. One hundred samples checked by Southern or dot-blot analysis using both radioactive and non-radioactive probes showed identical results in parallel with routine techniques. For molecular epidemiological studies we selected another set of conserved primers flanking the highly evolutive pseudogene (psi gene) region. This set was found to be efficient for all tested fixed rabies virus strains or wild rabies virus isolates as well as the rabies-related Mokola virus. We describe a progressive characterization of the strain that could be extended from rapid typing by a limited panel of restriction enzymes, to the ultimate identification of the nucleotide sequence by an original direct sequencing technique of amplified segments.

Evolution of European bat lyssaviruses.

Forty-seven European bat lyssaviruses (EBL) and two African insectivorous bat lyssaviruses (Duvenhage viruses) were selected for a comparison to be made of their evolutionary relationships. Studies were based on direct sequencing of the PCR-amplified products of the 400 nucleotides coding for the amino terminus of the nucleoprotein. Phylogenetic relationships were analysed after bootstrap resampling using the maximum parsimony and the neighbour-joining methods. Analyses of both the nucleotide and amino acid sequences placed these viruses in three separate clusters, namely genotype 4 (Duvenhage), genotype 5 (EBL1) and genotype 6 (EBL2). Evolutionary analysis of the nucleoprotein gene of EBL1 and EBL2 indicated low intrinsic heterogeneity mainly due to synonymous substitutions. In addition, both EBL1 and EBL2 evolved into at least two genetically distinguishable lineages (a and b) following geographical drifting. We can speculate that subsequently the lineages EBL1a and EBL1b were introduced into parts of northern Europe from two different geographical directions; EBL1b was probably introduced most recently and was from North Africa. Eptesicus serotinus appears to be the principal reservoir for EBL1 and Myotis dasycneme and M. daubentonii the reservoirs for EBL2.

Evaluation of High-Throughput Sequencing for Identifying Known and Unknown Viruses in Biological Samples

ABSTRACT High-throughput sequencing furnishes a large number of short sequence reads from uncloned DNA and has rapidly become a major tool for identifying viruses in biological samples, and in particular when the target sequence is undefined. In this study, we assessed the analytical sensitivity of a pipeline for detection of viruses in biological samples based on either the Roche-454 genome sequencer or Illumina genome analyzer platforms. We sequenced biological samples artificially spiked with a wide range of viruses with genomes composed of single or double-stranded DNA or RNA, including linear or circular single-stranded DNA. Viruses were added at a very low concentration most often corresponding to 3 or 0.8 times the validated level of detection of quantitative reverse transcriptase PCRs (RT-PCRs). For the viruses represented, or resembling those represented, in public nucleotide sequence databases, we show that the higher output of Illumina is associated with a much greater sensitivity, approaching that of optimized quantitative (RT-)PCRs. In this blind study, identification of viruses was achieved without incorrect identification. Nevertheless, at these low concentrations, the number of reads generated by the Illumina platform was too small to facilitate assembly of contigs without the use of a reference sequence, thus precluding detection of unknown viruses. When the virus load was sufficiently high, de novo assembly permitted the generation of long contigs corresponding to nearly full-length genomes and thus should facilitate the identification of novel viruses.

A Reliable Diagnosis of Human Rabies Based on Analysis of Skin Biopsy Specimens

BACKGROUND The number of human deaths due to rabies is currently underestimated to be 55,000 deaths per year. Biological diagnostic methods for confirmation of rabies remain limited, because testing on postmortem cerebral samples is the reference method, and in many countries, sampling brain tissue is rarely practiced. There is a need for a reliable method based on a simple collection of nonneural specimens. METHODS A new reverse-transcription, heminested polymerase chain reaction (RT-hnPCR) protocol was standardized at 3 participating centers in Cambodia, Madagascar, and France. Fifty-one patients from Cambodia, Madagascar, Senegal, and France were prospectively enrolled in the study; 43 (84%) were ultimately confirmed as having rabies. A total of 425 samples were collected from these patients during hospitalization. We studied the accuracy of the diagnosis by comparing the results obtained with use of biological fluid specimens (saliva and urine) and skin biopsy specimens with the results obtained with use of the standard rabies diagnostic procedure performed with a postmortem brain biopsy specimen. RESULTS The data obtained indicate a high specificity (100%) of RT-hnPCR and a higher sensitivity (>/=98%) when the RT-hnPCR was performed with skin biopsy specimens than when the test was performed with fluid specimens, irrespective of the time of collection (i.e., 1 day after the onset of symptoms or just after death). Also, a sensitivity of 100% was obtained with the saliva sample when we analyzed at least 3 successive samples per patient. CONCLUSIONS Skin biopsy specimens should be systematically collected in cases of encephalitis of unknown origin. These samples should be tested by RT-hnPCR immediately to confirm rabies; if the technique is not readily available locally, the samples should be tested retrospectively for epidemiological purposes.

European Bat Lyssavirus Infection in Spanish Bat Populations

From 1992 to 2000, 976 sera, 27 blood pellets, and 91 brains were obtained from 14 bat species in 37 localities in Spain. Specific anti-European bat lyssavirus 1 (EBL1)-neutralizing antibodies have been detected in Myotis myotis, Miniopterus schreibersii, Tadarida teniotis, and Rhinolophus ferrumequinum in the region of Aragon and the Balearic Islands. Positive results were also obtained by nested reverse transcription-polymerase chain reaction on brain, blood pellet, lung, heart, tongue, and esophagus-larynx-pharynx of M. myotis, Myotis nattereri, R. ferrumequinum, and M. schreibersii. Determination of nucleotide sequence confirmed the presence of EBL1 RNA in the different tissues. In one colony, the prevalence of seropositive bats over time corresponded to an asymmetrical curve, with a sudden initial increase peaking at 60% of the bats, followed by a gradual decline. Banded seropositive bats were recovered during several years, indicating that EBL1 infection in these bats was nonlethal. At least one of this species (M. schreibersii) is migratory and thus could be partially responsible for the dissemination of EBL1 on both shores of the Mediterranean Sea.

Phylodynamics and human-mediated dispersal of a zoonotic virus

Understanding the role of humans in the dispersal of predominately animal pathogens is essential for their control. We used newly developed Bayesian phylogeographic methods to unravel the dynamics and determinants of the spread of dog rabies virus (RABV) in North Africa. Each of the countries studied exhibited largely disconnected spatial dynamics with major geo-political boundaries acting as barriers to gene flow. Road distances proved to be better predictors of the movement of dog RABV than accessibility or raw geographical distance, with occasional long distance and rapid spread within each of these countries. Using simulations that bridge phylodynamics and spatial epidemiology, we demonstrate that the contemporary viral distribution extends beyond that expected for RABV transmission in African dog populations. These results are strongly supportive of human-mediated dispersal, and demonstrate how an integrated phylogeographic approach will turn viral genetic data into a powerful asset for characterizing, predicting, and potentially controlling the spatial spread of pathogens.

Dynamics of rabies virus quasispecies during serial passages in heterologous hosts

To understand the mutations and genetic rearrangements that allow rabies virus infections of new hosts and adaptation in nature, the quasispecies structure of the nucleoprotein and glycoprotein genes as well as two noncoding sequences of a rabies virus genome were determined. Gene sequences were obtained from the brain and from the salivary glands of the original host, a naturally infected European fox, and after serial passages in mice, dogs, cats and cell culture. A relative genetic stasis of the consensus sequences confirmed previous results about the stability of rabies virus. At the quasispecies level, the mutation frequency varies, in the following order: glycoprotein region (21.9 x 10(-4) mutations per bp), noncoding sequence nucleoprotein-phosphoprotein region (7.2-7.9 x 10(-4) mutations per bp) and nucleoprotein gene region (2.9-3.7 x 10(-4) mutations per bp). These frequencies varied according to the number, type of heterologous passages and the genomic region considered. The shape of the quasispecies structure was dramatically modified by passages in mice, in which the mutation frequencies increased by 12-31 x 10(-4) mutations per bp, depending on the region considered. Non-synonymous mutations were preponderant particularly in the glycoprotein gene, stressing the importance of positive selection in the maintenance and fixation of substitutions. Two mechanisms of genomic evolution of the rabies virus quasispecies, while adapting to environmental changes, have been identified: a limited accumulation of mutations with no replacement of the original master sequence and a less frequent but rapid selective overgrowth of favoured variants.

Evolution of Bovine Respiratory Syncytial Virus

ABSTRACT Until now, the analysis of the genetic diversity of bovine respiratory syncytial virus (BRSV) has been based on small numbers of field isolates. In this report, we determined the nucleotide and deduced amino acid sequences of regions of the nucleoprotein (N protein), fusion protein (F protein), and glycoprotein (G protein) of 54 European and North American isolates and compared them with the sequences of 33 isolates of BRSV obtained from the databases, together with those of 2 human respiratory syncytial viruses and 1 ovine respiratory syncytial virus. A clustering of BRSV sequences according to geographical origin was observed. We also set out to show that a continuous evolution of the sequences of the N, G, and F proteins of BRSV has been occurring in isolates since 1967 in countries where vaccination was widely used. The exertion of a strong positive selective pressure on the mucin-like region of the G protein and on particular sites of the N and F proteins is also demonstrated. Furthermore, mutations which are located in the conserved central hydrophobic part of the ectodomain of the G protein and which result in the loss of four Cys residues and in the suppression of two disulfide bridges and an α helix critical to the three-dimensional structure of the G protein have been detected in some recent French BRSV isolates. This conserved central region, which is immunodominant in BRSV G protein, thus has been modified in recent isolates. This work demonstrates that the evolution of BRSV should be taken into account in the rational development of future vaccines.

Phi29 polymerase based random amplification of viral RNA as an alternative to random RT-PCR

BackgroundPhi29 polymerase based amplification methods provides amplified DNA with minimal changes in sequence and relative abundance for many biomedical applications. RNA virus detection using microarrays, however, can present a challenge because phi29 DNA polymerase cannot amplify RNA nor small cDNA fragments (<2000 bases) obtained by reverse transcription of certain viral RNA genomes. Therefore, ligation of cDNA fragments is necessary prior phi29 polymerase based amplification. We adapted the QuantiTect Whole Transcriptome Kit (Qiagen) to our purposes and designated the method as Whole Transcriptome Amplification (WTA).ResultsWTA successfully amplified cDNA from a panel of RNA viruses representing the diversity of ribovirus genome sizes. We amplified a range of genome copy numbers from 15 to 4 × 107 using WTA, which yielded quantities of amplified DNA as high as 1.2 μg/μl or 1010 target copies. The amplification factor varied between 109 and 106. We also demonstrated that co-amplification occurred when viral RNA was mixed with bacterial DNA.ConclusionThis is the first report in the scientific literature showing that a modified WGA (WTA) approach can be successfully applied to viral genomic RNA of all sizes. Amplifying viral RNA by WTA provides considerably better sensitivity and accuracy of detection compared to random RT-PCR.