PubMed Health⌕ Search

Biomedical subjects

Philippe Lemey

Publications and source records attributed to Philippe Lemey.

At least 37 records · Page 2Linked to original sources
PubMed2005

Isolation and cloning of the raccoon (Procyon lotor) papillomavirus type 1 by using degenerate papillomavirus-specific primers.

Partial sequences of a novel papillomavirus were amplified from a cutaneous lesion biopsy of a raccoon (Procyon lotor), by using PCR with degenerate papillomavirus-specific primers. The Procyon lotor papillomavirus type 1 (PlPV-1) DNA was amplified with long template PCR in two overlapping fragments, together encompassing the entire genome, and the complete PlPV-1 genomic sequence was determined. The PlPV-1 genome consists of 8170 bp, and contains the typical papillomaviral open reading frames, encoding five early proteins and two late capsid proteins. Besides the classical non-coding region (NCR1) between the end of L1 and the start of E6, PlPV-1 contains an additional non-coding region (NCR2) of 1065 bp between the early and late protein region, which has previously also been described for the canine oral papillomavirus (COPV) and the Felis domesticus papillomavirus (FdPV-1). Phylogenetic analysis places PlPV-1 together with COPV and FdPV-1 in a monophyletic branch which encompasses the Lambda papillomavirus genus.

Animals↗
PubMed2005

Genetic variability and molecular evolution of the human respiratory syncytial virus subgroup B attachment G protein.

Human respiratory syncytial virus (HRSV) is the most important cause of acute respiratory disease in infants. Two major subgroups (A and B) have been identified based on antigenic differences in the attachment G protein. Antigenic variation between and within the subgroups may contribute to reinfections with these viruses by evading the host immune responses. To investigate the circulation patterns and mechanisms by which HRSV-B viruses evolve, we analyzed the G protein genetic variability of subgroup B sequences isolated over a 45-year period, including 196 Belgian strains obtained over 22 epidemic seasons (1982 to 2004). Our study revealed that the HRSV-B evolutionary rate (1.95 x 10(-3) nucleotide substitutions/site/year) is similar to that previously estimated for HRSV-A (1.83 x 10(-3) nucleotide substitutions/site/year). However, natural HRSV-B isolates appear to accommodate more drastic changes in their attachment G proteins. The most recent common ancestor of the currently circulating subgroup B strains was estimated to date back to around the year 1949. The divergence between the two major subgroups was calculated to have occurred approximately 350 years ago. Furthermore, we have identified 12 positively selected sites in the G protein ectodomain, suggesting that immune-driven selective pressure operates in certain codon positions. HRSV-A and -B strains have similar phylodynamic patterns: both subgroups are characterized by global spatiotemporal strain dynamics, where the high infectiousness of HRSV permits the rapid geographic spread of novel strain variants.

Amino Acid Sequence↗
PubMed2005

Molecular footprint of drug-selective pressure in a human immunodeficiency virus transmission chain.

Known human immunodeficiency virus (HIV) transmission histories are invaluable models for investigating the evolutionary and transmission dynamics of the virus and to assess the accuracy of phylogenetic reconstructions. Here we have characterized an HIV-1 transmission chain consisting of nine infected patients, almost all of whom were treated with antiviral drugs at later stages of infection. Partial pol and env gp41 regions of the HIV genome were directly sequenced from plasma viral RNA for at least one sample from each patient. Phylogenetic analyses in pol using likelihood methods inferred an evolutionary history not fully compatible with the known transmission history. This could be attributed to parallel evolution of drug resistance mutations resulting in the incorrect clustering of multidrug-resistant virus. On the other hand, a fully compatible phylogenetic tree was reconstructed from the env sequences. We were able to identify and quantify the molecular footprint of drug-selective pressure in pol using maximum likelihood inference under different codon substitution models. An increased fixation rate of mutations in the HIV population of the multidrug-resistant patient was demonstrated using molecular clock modeling. We show that molecular evolutionary analyses, guided by a known transmission history, can reveal the presence of confounding factors like natural selection and caution should be taken when accurate descriptions of HIV evolution are required.

Anti-HIV Agents↗
PubMed2005

Complete genomic sequence of human coronavirus OC43: molecular clock analysis suggests a relatively recent zoonotic coronavirus transmission event.

Coronaviruses are enveloped, positive-stranded RNA viruses with a genome of approximately 30 kb. Based on genetic similarities, coronaviruses are classified into three groups. Two group 2 coronaviruses, human coronavirus OC43 (HCoV-OC43) and bovine coronavirus (BCoV), show remarkable antigenic and genetic similarities. In this study, we report the first complete genome sequence (30,738 nucleotides) of the prototype HCoV-OC43 strain (ATCC VR759). Complete genome and open reading frame (ORF) analyses were performed in comparison to the BCoV genome. In the region between the spike and membrane protein genes, a 290-nucleotide deletion is present, corresponding to the absence of BCoV ORFs ns4.9 and ns4.8. Nucleotide and amino acid similarity percentages were determined for the major HCoV-OC43 ORFs and for those of other group 2 coronaviruses. The highest degree of similarity is demonstrated between HCoV-OC43 and BCoV in all ORFs with the exception of the E gene. Molecular clock analysis of the spike gene sequences of BCoV and HCoV-OC43 suggests a relatively recent zoonotic transmission event and dates their most recent common ancestor to around 1890. An evolutionary rate in the order of 4 x 10(-4) nucleotide changes per site per year was estimated. This is the first animal-human zoonotic pair of coronaviruses that can be analyzed in order to gain insights into the processes of adaptation of a nonhuman coronavirus to a human host, which is important for understanding the interspecies transmission events that led to the origin of the severe acute respiratory syndrome outbreak.

Animals↗
PubMed2004

A Bayesian statistical analysis of human T-cell lymphotropic virus evolutionary rates.

HTLV is a genetically-stable retrovirus that is considered to have evolved partly in concert with human migrations. Its rate of evolution is low and therefore, difficult to estimate reliably. In the first part of this study, we provide an improved estimate of HTLV evolutionary rate using anthropological calibration of phylogenetic nodes. We investigate two different anthropological calibrations using a Bayesian method that implements a relaxed molecular clock model and can combine data from multiple genes. The analysis shows that the two calibrations are compatible. In the second part, we develop a Bayesian statistical model to combine and compare the anthropology-based estimates of evolutionary rate with a rate recently calculated using pedigree data from vertically HTLV-infected families. We compare the statistical power of the two estimates and show that the current pedigree estimate, although resulting in considerably higher evolutionary rates, is too statistically weak to warrant a re-examination of the commonly used anthropology-based estimates. Statistical uncertainty burdens HTLV rate estimates based on both anthropological calibrations and on pedigree data; the former method rests on an untested assumption, whilst that latter is affected by small sample sizes.

Bayes Theorem↗
PubMed2004

HIV-1 compartmentalization in diverse leukocyte populations during antiretroviral therapy.

CD4+ T lymphocytes are the primary target of human immunodeficiency virus type 1 (HIV-1), but there is increasing evidence that other immune cells in the blood, including CD8+ T lymphocytes and monocytes, are also productively infected. The extent to which these additional cellular reservoirs contribute to ongoing immunodeficiency and viral persistence during therapy remains unclear. In this study, we conducted a detailed investigation of HIV-1 diversity and genetic structure in CD4+ T cells, CD8+ T cells, and monocytes of 13 patients receiving highly active antiretroviral therapy (HAART). Analysis of molecular variance and nonparametric tests performed on HIV-1 envelope sequences provided statistically significant evidence of viral compartmentalization in different leukocyte populations. Signature pattern analysis and predictions of coreceptor use provided no evidence that selection arising from viral tropism was responsible for the genetic structure observed. Analysis of viral genetic variation in different leukocyte populations demonstrated the action of founder effects as well as significant variation in the extent of genetic differentiation between subpopulations among patients. In the absence of evidence for leukocyte-specific selection, these features were supportive of a metapopulation model of HIV-1 replication as described previously among HIV-1 populations in the spleen. Compartmentalization of the virus in different leukocytes may have significant implications for current models of HIV-1 population genetics and contribute to the highly variable way in which drug resistance evolves in different individuals during HAART.

Anti-Retroviral Agents↗
PubMed2004

Plasma exchange-associated immunoglobulin m-negative hantavirus disease after a camping holiday in southern france.

A 37-year-old Belgian patient presented with acute nephropathia epidemica (NE) shortly after a camping holiday in southern France. Unusual symptoms were initial noncardiogenic lung involvement, followed by severe acute renal failure, acute acalculous cholecystitis, presence of immunoblasts in the bone marrow, and hemolytic anemia, presenting as hemolytic uremic syndrome. Positive immunoglobulin (Ig) A and rising IgG titers against Puumala hantavirus (PUUV) were detected, but IgM remained negative on days 8 and 20. The results of reverse-transcriptase-polymerase chain reaction performed on day 8 were positive for PUUV. This is the first report of an iatrogenically IgM-negative hantavirus case due to the selective removal of heavy-weight molecules during plasma exchange via the centrifugation technique. This is also the first report of proven NE from the Mediterranean part of France.

Adult↗
PubMed2004

Phylogenetic reconstruction of a known HIV-1 CRF04_cpx transmission network using maximum likelihood and Bayesian methods.

The CRF04_cpx strains of HIV-1 accounts for approximately 2-10% of the infected population in Greece, across different transmission risk groups. CRF04_cpx was the lineage documented in an HIV-1 transmission network in Thessalonica, northern Greece. Most of the transmissions occurred through unprotected heterosexual contacts between 1989 and 1993. Blood samples were available for six patients, obtained 6-10 years later, except for one patient sampled in 1991. Our objective was to examine whether the transmission history is compatible with the evolutionary tree of the virus, in partial gag, partial env, and partial gag+env. The inferred phylogenetic tree obtained using maximum likelihood and Bayesian methods in partial gag+env was much closer to the transmission tree than that using either env or gag separately. Our findings suggest that the epidemiological relationships among patients who have been infected by a common source correspond almost exactly to the evolutionary trees of the virus, given that enough phylogenetic signal is present in the alignment. Moreover, we found evidence that recombination is not the most parsimonious explanation for the phylogenetic incongruence between gag and env. For patients with known infection dates, the estimated dates of the coalescent events obtained using molecular clock calculations based on a newly developed Bayesian method in gag + env were in agreement with the actual infection dates.

Bayes Theorem↗
PubMed2004

Analysis of the serotype and genotype correlation of VP1 and the 5' noncoding region in an epidemiological survey of the human enterovirus B species.

The sequence identity of the enterovirus VP1 gene has been shown to correlate with the serotype concept. Enterovirus molecular typing methods are therefore often based on sequencing of the VP1 genomic region and monophyletic clustering of VP1 sequences of a homologous serotype. For epidemiological surveillance, 342 enterovirus samples obtained from patients with aseptic meningitis in Belgium from 1999 to 2002 were first diagnosed as being enterovirus positive by amplification of the 5' noncoding region (5'NCR) by reverse transcription (RT)-PCR. Subsequently, samples were molecularly typed by RT-nested PCR amplification and sequencing of a portion of the VP1 gene. Phylogenetic analyses were performed to investigate enteroviral evolution and to examine the serotype and genotype correlation of the two genomic regions. Our typing results demonstrated echovirus 30, echovirus 13, echovirus 18, and echovirus 6 to be the most predominant types. Echoviruses 13 and 18 were considered to be emerging human serotypes since 2000 and 2001, respectively, as they had been rarely reported before. Several serotypes existed as multiple genotypes (subtypes) from 1999 to 2002, but genomic differences mainly resided at synonymous sites; these results strongly suggest that the subtypes exhibit similar antigenic properties. Phylogenetic analyses confirmed that VP1 is an adequate region for molecular typing. Serotype-specific clusters are not observed commonly in phylogenetic trees based on the 5'NCR, and the phylogenetic signal in the 5'NCR was found to be particularly low. However, some substructure in the 5'NCR tree made a tentative prediction of the enterovirus type possible and was therefore helpful in PCR strategies for VP1 (e.g., primer choice), provided some background knowledge on the local spectrum of enteroviruses already exists.

Base Sequence↗
PubMed2004

Molecular evolution and circulation patterns of human respiratory syncytial virus subgroup a: positively selected sites in the attachment g glycoprotein.

Human respiratory syncytial virus (HRSV) is the most common etiological agent of acute lower respiratory tract disease in infants and can cause repeated infections throughout life. In this study, we have analyzed nucleotide sequences encompassing 629 bp at the carboxy terminus of the G glycoprotein gene for HRSV subgroup A strains isolated over 47 years, including 112 Belgian strains isolated over 19 consecutive years (1984 to 2002). By using a maximum likelihood method, we have tested the presence of diversifying selection and identified 13 positively selected sites with a posterior probability above 0.5. The sites under positive selection correspond to sites of O glycosylation or to amino acids that were previously described as monoclonal antibody-induced in vitro escape mutants. Our findings suggest that the evolution of subgroup A HRSV G glycoprotein is driven by immune pressure operating in certain codon positions located mainly in the second hypervariable region of the ectodomain. Phylogenetic analysis revealed the prolonged cocirculation of two subgroup A lineages among the Belgian population and the possible extinction of three other lineages. The evolutionary rate of HRSV subgroup A isolates was estimated to be 1.83 x 10(-3) nucleotide substitutions/site/year, projecting the most recent common ancestor back to the early 1940s.

Belgium↗
PubMed2004

The molecular population genetics of HIV-1 group O.

HIV-1 group O originated through cross-species transmission of SIV from chimpanzees to humans and has established a relatively low prevalence in Central Africa. Here, we infer the population genetics and epidemic history of HIV-1 group O from viral gene sequence data and evaluate the effect of variable evolutionary rates and recombination on our estimates. First, model selection tools were used to specify suitable evolutionary and coalescent models for HIV group O. Second, divergence times and population genetic parameters were estimated in a Bayesian framework using Markov chain Monte Carlo sampling, under both strict and relaxed molecular clock methods. Our results date the origin of the group O radiation to around 1920 (1890-1940), a time frame similar to that estimated for HIV-1 group M. However, group O infections, which remain almost wholly restricted to Cameroon, show a slower rate of exponential growth during the twentieth century, explaining their lower current prevalence. To explore the effect of recombination, the Bayesian framework is extended to incorporate multiple unlinked loci. Although recombination can bias estimates of the time to the most recent common ancestor, this effect does not appear to be important for HIV-1 group O. In addition, we show that evolutionary rate estimates for different HIV genes accurately reflect differential selective constraints along the HIV genome.

Bayes Theorem↗
PubMed2003

Site stripping based on likelihood ratio reduction is a useful tool to evaluate the impact of non-clock-like behavior on viral phylogenetic reconstructions.

The site stripping for clock detection procedure was implemented in the recently developed maximum likelihood framework for estimating evolutionary rates and divergence times in measurably evolving populations. The method was used to investigate the effect of rate variability on estimating divergence times in non-clock-like trees for human immunodeficiency viruses and hepatitis C viruses. We validate our approach by comparing dated coalescent nodes in molecular phylogenies with known dates of transmission. Our method was able to rapidly recover clock-like behavior and to indicate the presence and direction of a bias when estimates of divergence times using the unstripped data were flawed.

Algorithms↗
PubMed2003

Tracing the origin and history of the HIV-2 epidemic.

In this study we date the introduction of HIV-2 into the human population and estimate the epidemic history of HIV-2 subtype A in Guinea-Bissau, the putative geographic origin of HIV-2. The evolutionary history of the simian immunodeficiency virussooty mangabey/HIV-2 lineage was reconstructed by using available database sequences with known sampling dates, and a timescale for this history was calculated by using maximum likelihood methods. The date of the most recent common ancestor of HIV-2 subtype A strains was estimated to be 1940 +/- 16 and that of B strains was estimated to be 1945 +/- 14. In addition we used coalescent theory to estimate the past population dynamics of HIV-2 subtype A in a rural population of Guinea-Bissau. Parametric and nonparametric estimates of the effective number of infections through time were obtained for an equal sample of gag, pol, and env sequences. Our estimates of the epidemic history of HIV-2 subtype A in Guinea-Bissau show a transition from constant size to rapid exponential growth around 1955-1970. Our analysis provides evidence for a zoonotic transfer of HIV-2 during the first half of the 20th century and an epidemic initiation in Guinea-Bissau that coincides with the independence war (1963-1974), suggesting that war-related changes in sociocultural patterns had a major impact on the HIV-2 epidemic.

HIV Infections↗
PubMed2003

Molecular typing and epidemiology of enteroviruses identified from an outbreak of aseptic meningitis in Belgium during the summer of 2000.

Non-polio enteroviruses are the most common cause of aseptic meningitis worldwide. From May to September 2000, a major outbreak of aseptic meningitis occurred in Belgium. Cerebrospinal fluid samples (CSF) of 122 patients were found to contain enterovirus RNA using diagnostic RT-PCR that targeted a 231-bp gene fragment in the 5' noncoding region. In addition, a molecular typing method was developed based on RT-nested PCR and sequencing directly from CSF(a) 358-bp fragment in the aminoterminal part of the VP1 capsid protein. To identify the enterovirus type, nucleotide sequences of the VP1 amplicons were compared to all the enterovirus VP1 sequences available in GenBank. Echovirus 30 (31.2%), echovirus 13 (23.8%), and echovirus 6 (20.5%) were identified most frequently during the epidemic. Coxsackievirus B5 was present in 15.6% of the samples, and could be subdivided in two distinct epidemic clusters, coxsackievirus B5a (10.7%) and B5b (4.9%). Other enteroviruses encountered were echovirus 16 (5.7%), echovirus 18 (1.6%), coxsackievirus B4 (0.8%) and echovirus 7 (0.8%). The high prevalence of echovirus 13, considered previously a rare serotype, indicates it is an emerging epidemic type. To verify the typing results and to explore further the intratypical genetic variation, phylogenetic analysis was carried out. Geographical clustering of most of the strains within each type and subtype could be observed. The RT-nested PCR strategy, carried out directly on clinical samples, is a simple and rapid method for adequate molecular typing of the Group B enteroviruses causing aseptic meningitis.

Adolescent↗
PubMed2003

U.S. Human immunodeficiency virus type 1 epidemic: date of origin, population history, and characterization of early strains.

Human immunodeficiency virus (HIV) type 1 subtype B sequences (whole envelope and the p17 region of gag) were obtained from peripheral blood mononuclear cell samples collected in 1981 from seven HIV-infected U.S. individuals and in 1982 from one infected Canadian resident. Phylogenetic and nucleotide distance analyses were performed by using database sequences representing North American strains collected from 1978 to 1995. The estimated phylogeny was starlike, with early strains represented on different lineages. When sequences were grouped by years of collection, nucleotide distance comparisons demonstrated an increase in diversity over time and indicated that contemporary strains are more closely related to early epidemic strains than to each other. Using a recently developed likelihood ratio reduction procedure, the date of origin of the U.S. epidemic was estimated to be 1968 +/- 1.4 years. A coalescent approach was also used to estimate the population history of the U.S. subtype B epidemic. Our analyses provide new information that implies an exponential growth rate from the beginning of the U.S. HIV epidemic. The dating results suggest a U.S. introduction date (or date of divergence from the most recent common ancestor) that precedes the date of the earliest known AIDS cases in the late 1970s. Furthermore, the estimated epidemic growth curve shows a period of exponential growth that preceded most of the early documented cases and also indicates a leveling of prevalence rates in the recent past.

Disease Outbreaks↗
PubMed2002

Phylogenetic classification of TT virus groups based on the N22 region is unreliable.

Evolutionary studies on the recently discovered TT virus (TTV) are currently focused on approximately 220 nts of the N22 region, since this is the region for which most sequence information is available. Regarding the extensive sequence heterogeneity in this region, within group classification can be sufficiently reliable, however, between group classification becomes problematic. We observed high divergence at the nucleotide level between distant related strains (TTV groups) preventing unambiguous alignments, saturation in transitions within TTV groups and considerable phylogenetic noise due to conflicting signals within distinct genotypes. Consequently, analysing all TTV groups in one tree, using this 220 nt region provides unreliable results. Also within genotype analysis can produce conflicting results. On the other hand, this region is still suitable to some extent for within TTV group phylogenetic analysis. We suggest that care should be taken in future TTV phylogenetic analysis, in particular, larger and more conserved regions should be sequenced to allow between group comparisons.

Genotype↗
PubMed2002

Complete genome sequence, taxonomic assignment, and comparative analysis of the untranslated regions of the Modoc virus, a flavivirus with no known vector.

We recently developed a model for flavivirus infection in mice and hamsters using the Modoc virus (MODV), a flavivirus with no known vector (P. Leyssen, A. Van Lommel, C. Drosten, H. Schmitz, E. De Clercq, and J. Neyts, 2001, Virology 279, 27-37). We now present the coding and noncoding sequence of MODV. The Modoc virus genome was determined to be 10,600 nucleotides in length with a single open reading frame extending from nucleotides 110 to 10,234, encoding 3374 amino acids. The deduced gene order of the single open reading frame is C-prM-E-NS1-NS2A-NS2B-NS3-NS4A-NS4B-NS5, which is exactly the same as that of the mosquito- and tick-borne flaviviruses. It is flanked by a 5'- and 3'-untranslated region (UTR) of 109 and 366 nucleotides, respectively. Alignment of the MODV amino acid sequence with that of 20 other flaviviruses revealed several regions with high sequence similarity corresponding to functionally important domains (e.g., the serine protease/helicase/NTPase of NS3 and the methyltransferase/RNA-dependent RNA polymerase of NS5) and conserved sites for proteolytic cleavage by viral and cellular proteases. Phylogenetic analysis of the entire coding region confirmed the classification of MODV within the flaviviruses with no known vector, which is in agreement with previous findings based on partial NS5 sequences. A detailed comparative analysis of the putative folding patterns of the 5'- and 3'-UTR of MODV and of the tick- and mosquito-borne viruses was carried out. Structural elements in the 5'- and 3' UTR of MODV that are preserved among vector-borne flaviviruses were noted and so were structural elements distinguishing the MODV UTRs from mosquito-borne and tick-borne flaviviruses. Also the putative secondary structure of circularized MODV RNA is presented.

3' Untranslated Regions↗