Search

High-risk human papillomavirus (HPV) infections are responsible for cervical cancer. However, little is known about the differences between HPV types and risk categories regarding their genetic diversity and particularly APOBEC3-induced mutations - which contribute to the innate immune response to HPV. Using a capture-based next-generation sequencing, 156 HPV whole genome sequences covering 43 HPV types were generated from paired cervical and anal swabs of 30 Togolese female sex workers (FSWs) sampled in 2017. Genetic diversity and APOBEC3-induced mutations were assessed at the viral whole genome and gene levels. Thirty-four pairwise sequence comparisons covering 24 HPV types in cervical and anal swabs revealed identical infections in the two anatomical sites. Differences in genetic diversity among HPV types was observed between patients. The E6 gene was significantly less conserved in low-risk HPVs (lrHPVs) compared to high-risk HPVs (hrHPVs) (p = 0.009). APOBEC3-induced mutations were found to be more common in lrHPVs than in hrHPVs (p = 0.005), supported by our data and by using large HPV sequence collections from the GenBank database. Focusing on the most common lrHPVs 6 and 11 and hrHPVs 16 and 18, APOBEC3-induced mutations were predominantly found in the E4 and E6 genes in lrHPVs, but were almost absent in these genes in hrHPVs. The variable APOBEC3 mutational signatures could contribute to the different oncogenic potentials between HPVs. Further studies are needed to conclusively determine whether APOBEC3 editing levels are associated to the carcinogenic potential of HPVs at the type and sublineage scales., (© 2024 The Author(s). Journal of Medical Virology published by Wiley Periodicals LLC.)

Artemisinin-based combination therapies (ACTs) were introduced as the standard of care for uncomplicated malaria in Africa almost two decades ago. Recent studies in East Africa have reported a gradual increase in kelch13 ( k13 ) mutant parasites associated with reduced artesunate efficacy. As part of the Community Access to Rectal Artesunate for Malaria project, we collected blood samples from 697 children with signs of severe malaria in northern Uganda between 2018 and 2020, before and after the introduction of rectal artesunate (RAS) in 2019. K13 polymorphisms were assessed, and parasite editing and phenotyping were performed to assess the impact of mutations on parasite resistance. Whole-genome sequencing was performed, and haplotype networks were constructed to determine the geographic origin of k13 mutations. Of the 697 children, 540 were positive for Plasmodium falciparum malaria by PCR and were treated with either RAS or injectable artesunate monotherapy followed in most cases by ACT. The most common k13 mutation was C469Y (6.7%), which was detected more frequently in samples collected after RAS introduction. Genome editing confirmed reduced in vitro susceptibility to artemisinin in C469Y-harboring parasites compared to wild-type controls ( P < 0.001). The haplotypic network showed that flanking regions of the C469Y mutation shared the same African genetic background, suggesting a single and indigenous origin of the mutation. Our data provide evidence of selection for the artemisinin-resistant C469Y mutation. The realistic threat of multiresistant parasites emerging in Africa should encourage careful monitoring of the efficacy of artemisinin derivatives and strict adherence to ACT treatment regimens., Competing Interests: The authors declare no conflict of interest.

Background: Given the altered responses to both artemisinins and lumefantrine in Eastern Africa, monitoring antimalarial drug resistance in all African countries is paramount. Methods: We measured the susceptibility to six antimalarials using ex vivo growth inhibition assays (IC 50 ) for a total of 805 Plasmodium falciparum isolates obtained from travelers returning to France (2016-2023), mainly from West and Central Africa. Isolates were sequenced using molecular inversion probes (MIPs) targeting fourteen drug resistance genes across the parasite genome. Findings: Ex vivo susceptibility to several drugs has significantly decreased in 2019-2023 versus 2016-2018 parasite samples: lumefantrine (median IC 50 : 23·0 nM [IQR: 14·4-35·1] in 2019-2023 versus 13·9 nM [8·42-21·7] in 2016-2018, p<0·0001), monodesethylamodiaquine (35·4 [21·2-51·1] versus 20·3 nM [15·4-33·1], p<0·0001), and marginally piperaquine (20·5 [16·5-26·2] versus 18.0 [14·2-22·4] nM, p<0·0001). Only four isolates carried a validated pfkelch13 mutation. Multiple mutations in pfcrt and one in pfmdr1 (N86Y) were significantly associated with altered susceptibility to multiple drugs. The susceptibility to lumefantrine was altered by pfcrt and pfmdr1 mutations in an additive manner, with the wild-type haplotype ( pfcrt K76- pfmdr1 N86) exhibiting the least susceptibility. Interpretation: Our study on P. falciparum isolates from West and Central Africa indicates a low prevalence of molecular markers of artemisinin resistance but a significant decrease in susceptibility to the partner drugs that have been the most widely used since a decade -lumefantrine and amodiaquine. These phenotypic changes likely mark parasite adaptation to sustained drug pressure and call for intensifying the monitoring of antimalarial drug resistance in Africa.

Background: Given the altered responses to both artemisinins and lumefantrine in Eastern Africa, monitoring antimalarial drug resistance in all African countries is paramount., Methods: We measured the susceptibility to six antimalarials using ex vivo growth inhibition assays (IC 50 ) for a total of 805 Plasmodium falciparum isolates obtained from travelers returning to France (2016-2023), mainly from West and Central Africa. Isolates were sequenced using molecular inversion probes (MIPs) targeting fourteen drug resistance genes across the parasite genome., Findings: Ex vivo susceptibility to several drugs has significantly decreased in 2019-2023 versus 2016-2018 parasite samples: lumefantrine (median IC 50 : 23·0 nM [IQR: 14·4-35·1] in 2019-2023 versus 13·9 nM [8·42-21·7] in 2016-2018, p<0·0001), monodesethylamodiaquine (35·4 [21·2-51·1] versus 20·3 nM [15·4-33·1], p<0·0001), and marginally piperaquine (20·5 [16·5-26·2] versus 18.0 [14·2-22·4] nM, p<0·0001). Only four isolates carried a validated pfkelch13 mutation. Multiple mutations in pfcrt and one in pfmdr1 (N86Y) were significantly associated with altered susceptibility to multiple drugs. The susceptibility to lumefantrine was altered by pfcrt and pfmdr1 mutations in an additive manner, with the wild-type haplotype ( pfcrt K76- pfmdr1 N86) exhibiting the least susceptibility., Interpretation: Our study on P. falciparum isolates from West and Central Africa indicates a low prevalence of molecular markers of artemisinin resistance but a significant decrease in susceptibility to the partner drugs that have been the most widely used since a decade -lumefantrine and amodiaquine. These phenotypic changes likely mark parasite adaptation to sustained drug pressure and call for intensifying the monitoring of antimalarial drug resistance in Africa., Funding: This work was supported by the French Ministry of Health (grant to the French National Malaria Reference Center) and by the Agence Nationale de la Recherche (ANR-17-CE15-0013-03 to JC). JAB was supported by NIH R01AI139520. JR postdoctoral fellowship was funded by Institut de Recherche pour le Développement., Competing Interests: Conflict of interest The authors declare no conflicts of interest.

Background: Mutations in the Plasmodium falciparum dhfr gene confer resistance to pyrimethamine, which is widely used for malaria chemoprevention in Africa. We aimed to evaluate the frequency and evolution of dhfr mutations in Plasmodium ovale spp in Africa and their functional consequences, which are incompletely characterised., Methods: We analysed dhfr mutations and their frequencies in P ovale spp isolates collected between Feb 1, 2004, and Aug 31, 2023, from the French National Malaria Reference Centre collection and from field studies in Benin, Gabon, and Kenya. Genetic patterns of positive selection were investigated. Full-length recombinant wild-type and mutant DHFR enzymes from both P ovale curtisi and P ovale wallikeri were expressed in bacteria to test whether the most common mutations reduced pyrimethamine susceptibility., Findings: We included 518 P ovale spp samples (314 P ovale curtisi and 204 P ovale wallikeri). In P ovale curtisi, Ala15Ser-Ser58Arg was the most common dhfr mutation (39%; 124 of 314 samples). In P ovale wallikeri, dhfr mutations were less frequent, with Phe57Leu-Ser58Arg reaching 17% (34 of 204 samples). These two mutants were the most prevalent in central and east Africa and were fixed in Kenyan isolates. We detected six and four other non-synonymous mutations, representing 8% (24 isolates) and 2% (five isolates) of the P ovale curtisi and P ovale wallikeri isolates, respectively. Whole-genome sequencing and microsatellite analyses revealed reduced genetic diversity around the mutant pocdhfr and powdhfr genes. The mutant DHFR proteins showed structural changes at the pyrimethamine binding site in-silico, confirmed by a 4-times increase in pyrimethamine half-maximal inhibitory concentration in an Escherichia coli growth assay for the Phe57Leu-Ser58Arg mutant and 50-times increase for the Ala15Ser-Ser58Arg mutant, compared with the wild-type counterparts., Interpretation: The widespread use of sulfadoxine-pyrimethamine for malaria chemoprevention might have exerted fortuitous selection pressure for dhfr mutations in P ovale spp. This calls for closer monitoring of dhfr and dhps mutations in P ovale spp., Funding: French Ministry of Health, Agence Nationale de la Recherche, and Global Emerging Infections Surveillance branch of the Armed Forces Health Surveillance Division., Competing Interests: Declaration of interests We declare no competing interests., (Copyright © 2024 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

Introduction: In Senegal, molecular diagnosis was widely used for the detection and management of COVID-19 patients. However, genomic surveillance was very limited in the public sector. This study aimed to share the experience of a Senegalese public sector laboratory in response to the COVID-19 pandemic, and to describe the distribution of variants circulating in 2020 and 2021., Methodology: From July 2020 to December 2021, SARS-CoV-2 qRT-PCR was performed on nasopharyngeal samples from travelers and symptomatic patients at the Bacteriology and Virology Laboratory (LBV) of the Aristide le Dantec University Teaching Hospital. Samples with a cycle threshold (Ct) ≤ 30 were selected for whole-genome sequencing (WGS) using the Nanopore technology. In-house scripts were developed to study the spatial and temporal distribution of SARS-CoV-2 variants in Senegal, using our sequences and those retrieved from the GISAID database., Results: Of 8,207 patients or travelers screened for SARS-CoV-2, 970 (11.8%) were positive and 386 had a Ct ≤ 30. WGS was performed on 133 samples. Concomitantly with high-quality sequences deposited in the GISAID database covering nine cities in Senegal in 2020 and 2021 (n = 1,539), we observed a high circulation of the 20A (B.1, B.1.416 and B.1.620) and 20B (B.1.1.420) lineages in 2020, while most of the samples belonged to Delta variants (AY34 and AY.34.1, 22%) in 2021., Conclusions: Despite its late involvement, COVID-19 diagnosis was routinely performed in LBV, but genomic characterization remained challenging. The genomic diversity of SARS-CoV-2 strains in Senegal reflected that observed worldwide during the first waves of the pandemic., Competing Interests: No Conflict of Interest is declared, (Copyright (c) 2024 Adjiratou Aissatou Ba, Romain Coppée, Assane Dieng, Jarra Manneh, Mengue Fall, Khadim Gueye, Pauline Yacine Sene, Samba Ndiour, Dianké Samaté, Pascaline Manga, Oumy Diop Diongue, Abdoulaye Camara, Abdou Padane, Sada Diallo; Sokhna Bousso Gueye; Makhtar Camara, Abdul Karim Sesay, Coumba Toure-Kane, Souleymane Mboup, Charlotte Charpentier, Halimatou Diop-Ndiaye.)

Background: Advances in high-throughput sequencing (HTS) technologies and reductions in sequencing costs have revolutionised the study of genomics and molecular biology by making whole-genome sequencing (WGS) accessible to many laboratories. However, the analysis of WGS data requires significant computational effort, which is the major drawback in implementing WGS as a routine laboratory technique., Objective: Automated pipelines have been developed to overcome this issue, but they do not exist for all organisms. This is the case for human respiratory syncytial virus (RSV), which is a leading cause of lower respiratory tract infections in infants, the elderly, and immunocompromised adults., Results: We present RSV-GenoScan, a fast and easy-to-use pipeline for WGS analysis of RSV generated by HTS on Illumina or Nanopore platforms. RSV-GenoScan automates the WGS analysis steps directly from the raw sequence data. The pipeline filters the sequence data, maps the reads to the RSV reference genomes, generates a consensus sequence, identifies the RSV subgroup, and lists amino acid mutations, insertions and deletions in the F and G viral genes. This enables the rapid identification of mutations in these coding genes that are known to confer resistance to monoclonal antibodies., Availability: RSV-GenoScan is freely available at https://github.com/AlexandreD-bio/RSV-GenoScan., Competing Interests: Declaration of Competing Interest The authors declare that they have no competing interests regarding the current work., (Copyright © 2024 The Authors. Published by Elsevier B.V. All rights reserved.)

Background: Few data are available on the real-world efficacy of receiving tenofovir-lamivudine-dolutegravir (DTG) as HIV treatment, particularly among young people in West Africa. Here, we evaluated pharmaco-virological outcomes and resistance profiles among Togolese children and adolescents., Methods: A cross-sectional study was conducted in Lomé, Togo, enrolling antiretroviral-treated people with HIV aged from 18 months to 24 years. Plasma HIV-1 viral load and antiretroviral concentrations were measured. Next-Generation Sequencing (NGS) of protease, Reverse Transcriptase (RT) and integrase was performed on all samples with viral load >200 c/mL. Drug resistance mutations (DRMs) were identified and interpreted using the ANRS-MIE algorithm., Results: 264 participants were enrolled (median age=17 years), 226 received a DTG-based regimen for a median of 20.5 months. Among them, virological suppression at the 200 c/mL threshold in 80.0% of the participants. Plasma DTG concentrations were adequate (i.e., >640 ng/mL), suboptimal and below the limit of quantification in 74.1%, 6.7% and 19.2% of participants receiving DTG, respectively. Overall, viruses resistant to any of Nucleoside RT Inhibitors, Non-NRTIs, and protease inhibitors were found in 52%, 66% and 1.6% of participants, respectively. A major integrase inhibitor DRM was observed in 9.4% (n=3/32, R263K, E138A-G140A-Q148R, and N155H) of participants with a viral load >200 c/mL., Conclusions: These first findings in such a large series of adolescents in a low-income country, showed a good virological response of 80% and the presence of an integrase DRM in 9.4% of the virological failures, supporting the need to monitor DTG drug resistance to reduce the risk of resistance acquisition., (© The Author(s) 2024. Published by Oxford University Press on behalf of Infectious Diseases Society of America. All rights reserved. For commercial re-use, please contact reprints@oup.com for reprints and translation rights for reprints. All other permissions can be obtained through our RightsLink service via the Permissions link on the article page on our site—for further information please contact journals.permissions@oup.com.)

Background: A severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection lasts longer in immunocompromised hosts than in immunocompetent patients. Prolonged infection is associated with a higher probability of selection for novel SARS-CoV-2 mutations, particularly in the spike protein, a critical target for vaccines and therapeutics., Methods: From December 2020 to September 2022, respiratory samples from 444 immunocompromised patients and 234 health care workers positive for SARS-CoV-2, diagnosed at 2 hospitals in Paris, France, were analyzed using whole-genome sequencing using Nanopore technology. Custom scripts were developed to assess the SARS-CoV-2 genetic diversity between the 2 groups and within the host., Results: Most infections were SARS-CoV-2 Delta or Omicron lineages. Viral genetic diversity was significantly higher in infections of immunocompromised patients than those of controls. Minor mutations were identified in viruses sequenced from immunocompromised individuals, which became signature mutations for newer SARS-CoV-2 variants as the epidemic progressed. Two patients were coinfected with Delta and Omicron variants. The follow-up of immunocompromised patients revealed that the SARS-CoV-2 genome evolution differed in the upper and lower respiratory tracts., Conclusions: This study found that SARS-CoV-2 infection in immunocompromised patients is associated with higher genetic diversity, which could lead to the emergence of new SARS-CoV-2 variants with possible immune evasion or different virulence characteristics., Competing Interests: Potential conflicts of interest. All authors: No reported conflicts. All authors have submitted the ICMJE Form for Disclosure of Potential Conflicts of Interest. Conflicts that the editors consider relevant to the content of the manuscript have been disclosed., (© The Author(s) 2023. Published by Oxford University Press on behalf of Infectious Diseases Society of America. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Background: Efficacy of sulfadoxine-pyrimethamine, the malaria chemoprophylaxis used in pregnant women, and in children when combined with amodiaquine, is threatened by the accumulation of mutations in the Plasmodium falciparum dihydropteroate synthase (pfdhps) and dihydrofolate reductase (pfdhfr) genes. Data on the prevalence of resistant alleles in central Africa and the new pfdhps I431V mutation, particularly associated with other mutations to form the pfdhps vagKgs allele, are scarce. We explored the frequency and geographical distribution of pfdhps and pfdhfr mutations in central Africa in 2014-18, and assessed the evolutionary origin of the vagKgs allele., Methods: Samples were collected at 18 health-care centres in seven countries (Angola, Cameroon, Central African Republic, Democratic Republic of the Congo, Gabon, Nigeria, and Republic of the Congo) from patients who showed possible symptoms of malaria between March 1, 2014, and Oct 31, 2018. Samples that were positive for P falciparum were transported to a laboratory in Toulouse, France, and genotyped. The frequency of pfdhfr and pfdhps mutations was studied in 1749 samples. Microsatellites in pfdhps flanking regions and whole-genome analysis compared with parasite genomes from the data-sharing network MalariaGEN were performed on samples carrying the vagKgs allele., Findings: Mapping of the prevalence of single nucleotide polymorphisms and corresponding alleles of pfdhfr and pfdhps showed a substantial spread of alleles associated with sulfadoxine-pyrimethamine resistance in central Africa during the 2014-18 period, especially an increase going west to east in pfdhps alleles carrying the K540E and A581G mutations. A high prevalence of the pfdhps I431V mutation was observed in Cameroon (exceeding 50% in the northern region) and Nigeria. Genomic analysis showed a recent African emergence and a clonal expansion of the most frequent pfdhps vagKgs allele., Interpretation: Reduced sulfadoxine-pyrimethamine efficacy due to increased resistance is a worrying situation, especially because the malaria transmission level is high in central Africa. Although the resistance phenotype remains to be confirmed, the emergence and spread of the vagKgs allele in west and central Africa could challenge the use of sulfadoxine-pyrimethamine., Funding: Toulouse Institute for Infectious and Inflammatory Diseases., Competing Interests: Declaration of interests We declare no competing interests., (Copyright © 2023 The Author(s). Published by Elsevier Ltd. This is an Open Access article under the CC BY-NC-ND 4.0 license. Published by Elsevier Ltd.. All rights reserved.)

Background: Hypermutated viruses induced by APOBEC3 (apolipoprotein B mRNA-editing, enzyme-catalytic, polypeptide-like 3) proteins comprise some of the defective viruses in the HIV reservoir. Here, we assessed the proportion of APOBEC3-induced defective proviruses in HIV-positive patients before and after receiving dolutegravir + lamivudine dual therapy., Methods: PBMCs of virologically suppressed patients enrolled in the ANRS 167 LAMIDOL trial, evaluating a switch from triple therapy to dolutegravir + lamivudine, were collected 8 weeks before (W-8) and 48 weeks after (W48) dual-therapy initiation. The Vif and RT regions were subject to next-generation sequencing. Bioinformatic algorithms were developed to identify APOBEC3-defective sequences and APOBEC3-related drug resistance mutations (APOMuts). All hypermutated sequences and those containing at least one stop codon were considered as defective., Results: One hundred and four patients were enrolled (median virological suppression duration: 4.2 years; IQR: 2.0-9.1). Proviral defective reads at W-8 and W48 were detected in Vif in 22% and 29% of patients, respectively, and in RT in 38% and 42% of patients, respectively. At least one APOMut was present in proviruses of 27% and 38% of patients at W-8 and W48, respectively. The ratio of APOMuts/number of potential APOMut sites was significantly higher at W48 (16.5%) than at W-8 (9.8%, P = 0.007). The presence of APOBEC3-defective viruses at W-8 was not associated with HIV total DNA level, nor with the third drug class received prior to switching to dolutegravir + lamivudine, nor with the duration of virological suppression., Conclusions: Whereas no significant change in the proportion of patients with APOBEC3-defective proviruses was evidenced after 1 year of dolutegravir + lamivudine maintenance, enrichment in APOMuts was observed. Further longer-term studies are needed to assess the other forms of defective viruses with dual-therapy., (© The Author(s) 2023. Published by Oxford University Press on behalf of British Society for Antimicrobial Chemotherapy. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Like Plasmodium vivax, both Plasmodium ovale curtisi and Plasmodium ovale wallikeri have the ability to cause relapse in humans, defined as recurring asexual parasitemia originating from liver-dormant forms subsequent to a primary infection. Here, we investigated relapse patterns in P ovale wallikeri infections from a cohort of travelers who were exposed to the parasite in sub-Saharan Africa and then experienced relapses after their return to France. Using a novel set of 8 highly polymorphic microsatellite markers, we genotyped 15 P ovale wallikeri relapses. For most relapses, the paired primary and relapse infections were highly genetically related (with 12 being homologous), an observation that was confirmed by whole-genome sequencing for the 4 relapses we further studied. This is, to our knowledge, the first genetic evidence of relapses in P ovale spp., Competing Interests: Potential conflicts of interest. The authors: No reported conflicts of interest. All authors have submitted the ICMJE Form for Disclosure of Potential Conflicts of Interest. Conflicts that the editors consider relevant to the content of the manuscript have been disclosed., (© The Author(s) 2023. Published by Oxford University Press on behalf of Infectious Diseases Society of America. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

In May 2022, several countries reported mpox cases from patients without history of traveling to endemic areas. France was one of the most affected European countries by this outbreak. In this study, the clinical characteristics of mpox cases in France were described, and the genetic diversity of the virus was studied. Patients diagnosed with mpox infection (quantitative polymerase chain reaction c t  < 28) between May 21, and July 4, 2022 and between 16th August and 10th September 2022 were included to this study. Twelve amplicons corresponding to the most polymorphic regions of the mpox genome and covering ~30 000 nucleotides were generated and sequenced using the S5 XL Ion Torrent technology to evaluate the genetic diversity of mpox sequences. One hundred and forty-eight patients were diagnosed with mpox-infection. 95% were men, 5% transgender (M-to-F), 50% were taking human immunodeficiency virus (HIV) pre-exposure prophylaxis, and 25% were HIV seropositive. One hundred and sixty-two samples (some patients had two samples) were sequenced and compared to GenBank sequences. Overall, low genetic diversity of mpox sequences was found compared with pre-epidemic Western-African sequences, with 32 distinct mutational patterns. This study provides a first glance at the mutational landscape of early mpox 2022 circulating strains in Paris (France)., (© 2023 The Authors. Journal of Medical Virology published by Wiley Periodicals LLC.)

Background: Malaria control is highly dependent on the effectiveness of artemisinin-based combination therapy (ACT), the current frontline malaria curative treatment. Unfortunately, the emergence and spread of parasites resistant to artemisinin (ART) derivatives in Southeast Asia and South America, and more recently in Rwanda and Uganda (East Africa), compromise their long-term use in sub-Saharan Africa, where most malaria deaths occur., Methods: Here, ex vivo susceptibility to dihydroartemisinin (DHA) was evaluated from 38 Plasmodium falciparum isolates collected in 2017 in Thiès (Senegal) expressed in the Ring-stage Survival Assay (RSA). Both major and minor variants were explored in the three conserved-encoding domains of the pfkelch13 gene, the main determinant of ART resistance using a targeted-amplicon deep sequencing (TADS) approach., Results: All samples tested in the ex vivo RSA were found to be susceptible to DHA (parasite survival rate < 1%). The non-synonymous mutations K189T and K248R in pfkelch13 were observed each in one isolate, as major (99%) or minor (5%) variants, respectively., Conclusion: The results suggest that ART is still fully effective in the Thiès region of Senegal in 2017. Investigations combining ex vivo RSA and TADS are a useful approach for monitoring ART resistance in Africa., (© 2023. The Author(s).)

Although France was one of the most affected European countries by the COVID-19 pandemic in 2020, the dynamics of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) movement within France, but also involving France in Europe and in the world, remain only partially characterized in this timeframe. Here, we analyzed GISAID deposited sequences from January 1 to December 31, 2020 ( n = 638,706 sequences at the time of writing). To tackle the challenging number of sequences without the bias of analyzing a single subsample of sequences, we produced 100 subsamples of sequences and related phylogenetic trees from the whole dataset for different geographic scales (worldwide, European countries, and French administrative regions) and time periods (from January 1 to July 25, 2020, and from July 26 to December 31, 2020). We applied a maximum likelihood discrete trait phylogeographic method to date exchange events (i.e., a transition from one location to another one), to estimate the geographic spread of SARS-CoV-2 transmissions and lineages into, from and within France, Europe, and the world. The results unraveled two different patterns of exchange events between the first and second half of 2020. Throughout the year, Europe was systematically associated with most of the intercontinental exchanges. SARS-CoV-2 was mainly introduced into France from North America and Europe (mostly by Italy, Spain, the United Kingdom, Belgium, and Germany) during the first European epidemic wave. During the second wave, exchange events were limited to neighboring countries without strong intercontinental movement, but Russia widely exported the virus into Europe during the summer of 2020. France mostly exported B.1 and B.1.160 lineages, respectively, during the first and second European epidemic waves. At the level of French administrative regions, the Paris area was the main exporter during the first wave. But, for the second epidemic wave, it equally contributed to virus spread with Lyon area, the second most populated urban area after Paris in France. The main circulating lineages were similarly distributed among the French regions. To conclude, by enabling the inclusion of tens of thousands of viral sequences, this original phylodynamic method enabled us to robustly describe SARS-CoV-2 geographic spread through France, Europe, and worldwide in 2020., Competing Interests: RC, FB, AJ, VL, VF, AF, LD, CC, SL, KZ, VC, DD, AM, BV, AB No competing interests declared, (© 2023, Coppée et al.)

Background: Exposure of blood to malaria parasites can lead to infection even in the absence of the mosquito vector. During a stay in a healthcare facility, accidental inoculation of the skin with blood from a malaria patient might occur, referred to as nosocomial malaria., Methods: Between 2007 and 2021, we identified 6 autochthonous malaria cases that occurred in different French hospitals, originating from nosocomial transmission and imported malaria cases being the infection source. Four cases were observed during the coronavirus disease 2019 pandemic. The genetic relatedness between source and nosocomial infections was evaluated by genome-wide short tandem repeats (STRs) and single-nucleotide polymorphisms (SNPs)., Results: None of the patients with autochthonous malaria had travel history to an endemic area nor had been transfused. For each case, both the source and recipient patients stayed a few hours in the same ward. After diagnosis, autochthonous cases were treated with antimalarials and all recovered except 1. Genetically, each pair of matched source/nosocomial parasite infections showed <1% of different STRs and <6.9% (<1.5% for monoclonal infections) of different SNPs. Similar levels of genetic differences were obtained for parasite DNA samples that were independently sequenced twice as references of identical infections. Parasite phylogenomics were consistent with travel information reported by the source patients., Conclusions: Our study demonstrates that genomics analyses may resolve nosocomial malaria transmissions, despite the uncertainty regarding the modes of contamination. Nosocomial transmission of potentially life-threatening parasites should be taken into consideration in settings or occasions where compliance with universal precautions is not rigorous., Competing Interests: Potential conflicts of interest . The authors: No reported conflicts of interest. All authors have submitted the ICMJE Form for Disclosure of Potential Conflicts of Interest., (© The Author(s) 2022. Published by Oxford University Press on behalf of Infectious Diseases Society of America. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

No data about antiretroviral (ARV) treatment coverage and virological response are available among key populations (female sex workers [FSW] and Men having Sex with Men [MSM]) in Togo. This study aimed to both describe Human Immunodeficiency Virus (HIV) immunovirological status and evaluate the pertinence of an original algorithm combining pharmacology (PK) and viral load (VL) to identify subjects at risk of ARV drug resistance. A cross-sectional multicentric study was conducted in 2017 in Togo. Our PK-virological algorithm (PK-VA) defines subjects at risk of resistance when exhibiting both detectable plasma drug concentrations and VL > 200 c/mL. Among the 123 FSW and 136 MSM included, 50% and 66% were receiving ARV, with 69% and 80% of them successfully-treated, respectively. Genotypes showed drug-resistance mutation in 58% and 63% of nonvirologically controlled (VL > 200 c/mL) ARV-treated FSW and MSM, respectively. PK-VA would have enabled to save 75% and 72% of genotypic tests, for FSW and MSM, respectively. We reported first data about HIV care cascade among key populations in Togo, highlighting they are tested for HIV but linkage to care remains a concern. Furthermore, 70%-80% of ARV-treated participants experienced virological success. In limited resources settings, where genotyping tests are beyond reach, PK-VA might be an easiest solution to sort out patients needing ARV adaptation due to resistance., (© 2023 The Authors. Journal of Medical Virology published by Wiley Periodicals LLC.)

Plasmodium sporozoites are transmitted to a mammalian host during blood feeding by an infected mosquito and invade hepatocytes for initial replication of the parasite into thousands of erythrocyte-invasive merozoites. Here we report that the B9 protein, a member of the 6-cysteine domain protein family, is secreted from sporozoite micronemes and is required for productive invasion of hepatocytes. The N-terminus of B9 forms a beta-propeller domain structurally related to CyRPA, a cysteine-rich protein forming an essential invasion complex in Plasmodium falciparum merozoites. The beta-propeller domain of B9 is essential for sporozoite infectivity and interacts with the 6-cysteine proteins P36 and P52 in a heterologous expression system. Our results suggest that, despite using distinct sets of parasite and host entry factors, Plasmodium sporozoites and merozoites may share common structural modules to assemble protein complexes for invasion of host cells., Competing Interests: The authors declare no competing interests., (© 2023 The Author(s).)

Human respiratory syncytial virus (RSV) is responsible of lower respiratory tract infections which may be severe in infants, elderly and immunocompromised adults. Europe and North-American countries have observed a massive reduction of RSV incidence during the 2020-2021 winter season. Using a systematic RSV detection coupled to SARS-CoV-2 for all adult patients admitted at the Foch hospital (Suresnes, France) between January and March 2021 (n = 11,324), only eight RSV infections in patients with prolonged RNA shedding were diagnosed. RSV whole-genome sequencing revealed that six and two patients were infected by RSV groups A and B, respectively. RSV carriage lasted from 7 to at least 30 days disregarding of RSV lineage. The most prolonged RSV shedding was observed in an asymptomatic patient. We detected novel patient-specific non-synonymous mutations in the G glycoprotein gene, including a double identical mutation in the repeated region for one patient. No additional mutation occurred in the RSV genome over the course of infection in the four patients tested for. In conclusion, our results suggest that the temporal shift in the RSV epidemic is not likely to be explained by the emergence of a high frequency, unreported variant. Moreover, prolonged RSV carriages in asymptomatic patients could play a role in virus spread., Competing Interests: Declaration of interests The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2022. Published by Elsevier B.V.)

Since 2010, the human-infecting malaria parasite Plasmodium ovale spp. has been divided into two genetically distinct species, P. ovale wallikeri and P. ovale curtisi. In recent years, application of whole-genome sequencing (WGS) to P. ovale spp. allowed to get a better understanding of its evolutionary history and discover some specific genetic patterns. Nevertheless, WGS data from P. ovale spp. are still scarce due to several drawbacks, including a high level of human DNA contamination in blood samples, infections with commonly low parasite density, and the lack of robust in vitro culture. Here, we developed two selective whole-genome amplification (sWGA) protocols that were tested on six P. ovale wallikeri and five P. ovale curtisi mono-infection clinical samples. Blood leukodepletion by a cellulose-based filtration was used as the gold standard for intraspecies comparative genomics with sWGA. We also demonstrated the importance of genomic DNA preincubation with the endonuclease McrBC to optimize P. ovale spp. sWGA. We obtained high-quality WGS data with more than 80% of the genome covered by ≥5 reads for each sample and identified more than 5,000 unique single-nucleotide polymorphisms (SNPs) per species. We also identified some amino acid changes in pocdhfr and powdhfr for which similar mutations in P. falciparum and P. vivax are associated with pyrimethamine or cycloguanil resistance. In conclusion, we developed two sWGA protocols for P. ovale spp. WGS that will help to design much-needed large-scale P. ovale spp. population studies. IMPORTANCE Plasmodium ovale spp. has the ability to cause relapse, defined as recurring asexual parasitemia originating from liver-dormant forms. Whole-genome sequencing (WGS) data are of importance to identify putative molecular markers associated with relapse or other virulence mechanisms. Due to low parasitemia encountered in P. ovale spp. infections and no in vitro culture available, WGS of P. ovale spp. is challenging. Blood leukodepletion by filtration has been used, but no technique exists yet to increase the quantity of parasite DNA over human DNA when starting from genomic DNA extracted from whole blood. Here, we demonstrated that selective whole-genome amplification (sWGA) is an easy-to-use protocol to obtain high-quality WGS data for both P. ovale spp. species from unprocessed blood samples. The new method will facilitate P. ovale spp. population genomic studies.

Monoclonal antibodies (MAbs) targeting the Spike glycoprotein of SARS-CoV-2 is a key strategy to prevent severe COVID-19. Here, the efficacy of two monoclonal antibody bitherapies against SARS-CoV-2 was assessed on 92 patients at high risk of severe COVID-19 between March and October 2021 (Bichat-Claude Bernard Hospital, Paris, France). Nine patients died despite appropriate management. From 14 days following treatment initiation, we observed a slower viral load decay for patients treated with the bitherapy Bamlanivimab/Etsevimab compared to the Casirivimab/Imdevimab association therapy ( P  = 0.045). The emergence of several mutations on the Spike protein known to diminish antiviral efficacy was observed from 1 to 3 weeks after infusion. The Q493R mutation was frequently selected, located in a region of joint structural overlap by Bamlanivimab/Etsevimab antibodies. Despite that this study was done on former SARS-CoV-2 variants (Alpha and Delta), the results provide new insights into resistance mechanisms in SARS-CoV-2 antibodies neutralization escape and should be considered for current and novel variants. IMPORTANCE Monoclonal antibody bitherapies (MAbs) are commonly prescribed to treat severe SARS-CoV-2-positive patients, and the rapid growth of resistance mutation emergence is alarming globally. To explore this issue, we conducted both clinical and genomic analyses of SARS-CoV-2 in a series of patients treated in 2021. We first noticed that the two dual therapies prescribed during the study had different kinetics of viral load decay. Rapidly after initiation of the treatments, resistance mutations emerged in the interface between the MAbs and the target Spike glycoprotein, demonstrating the importance to continuously screen the viral genome during treatment course. Taken together, the results highlight that viral mutations may emerge under selective pressure, conferring a putative competitive advantage, and could rapidly spread, as observed for the Omicron variant.

A returned traveler to Uganda presented with a Plasmodium falciparum kelch13 A675V mutant infection that exhibited delayed clearance under artesunate therapy. Parasites were genetically related to recently reported Ugandan artemisinin-resistant A675V parasites. Adequate malaria prevention measures and clinical and genotypic surveillance are important tools to avoid and track artemisinin resistance., Competing Interests: Potential conflicts of interest. All authors: No reported conflicts. All authors have submitted the ICMJE Form for Disclosure of Potential Conflicts of Interest., (© The Author(s) 2022. Published by Oxford University Press for the Infectious Diseases Society of America. All rights reserved. For permissions, e-mail: journals.permissions@oup.com.)

Background: Whole-genome sequencing (WGS) is becoming increasingly helpful to assist malaria control programmes. A major drawback of this approach is the large amount of human DNA compared to parasite DNA extracted from unprocessed whole blood. As red blood cells (RBCs) have a diameter of about 7-8 µm and exhibit some deformability, it was hypothesized that cheap and commercially available 5 µm filters might retain leukocytes but much less of Plasmodium falciparum-infected RBCs. This study aimed to test the hypothesis that such a filtration method, named 5WBF (for 5 µm Whole Blood Filtration), may provide highly enriched parasite material suitable for P. falciparum WGS., Methods: Whole blood was collected from five patients experiencing a P. falciparum malaria episode (ring-stage parasitaemia range: 0.04-5.5%) and from mock samples obtained by mixing synchronized, ring-stage cultured P. falciparum 3D7 parasites with uninfected human whole blood (final parasitaemia range: 0.02-1.1%). These whole blood samples (50 to 400 µL) were diluted in RPMI 1640 medium or PBS 1× buffer and filtered with a syringe connected to a 5 µm commercial filter. DNA was extracted from 5WBF-treated and unfiltered counterpart blood samples using a commercial kit. The 5WBF method was evaluated on the ratios of parasite:human DNA assessed by qPCR and by sequencing depth and percentages of coverage from WGS data (Illumina NextSeq 500). As a comparison, the popular selective whole-genome amplification (sWGA) method, which does not rely on blood filtration, was applied to the unfiltered counterpart blood samples., Results: After applying 5WBF, qPCR indicated an average of twofold loss in the amount of parasite template DNA (Pf ARN18S gene) and from 4096- to 65,536-fold loss of human template DNA (human β actin gene). WGS analyses revealed that > 95% of the  parasite nuclear and organellar genomes were all covered at ≥ 10× depth for all samples tested. In sWGA counterparts, the organellar genomes were poorly covered and from 47.7 to 82.1% of the nuclear genome was covered at ≥ 10× depth depending on parasitaemia. Sequence reads were homogeneously distributed across gene sequences for 5WBF-treated samples (n = 5460 genes; mean coverage: 91×; median coverage: 93×; 5th percentile: 70×; 95th percentile: 103×), allowing the identification of gene copy number variations such as for gch1. This later analysis was not possible for sWGA-treated samples, as a much more heterogeneous distribution of reads across gene sequences was observed (mean coverage: 80×; median coverage: 51×; 5th percentile: 7×; 95th percentile: 245×)., Conclusions: The novel 5WBF leucodepletion method is simple to implement and based on commercially available, standardized 5 µm filters which cost from 1.0 to 1.7€ per unit depending on suppliers. 5WBF permits extensive genome-wide analysis of P. falciparum ring-stage isolates from minute amounts of whole blood even with parasitaemias as low as 0.02%., (© 2022. The Author(s).)

Partial artemisinin resistance, defined in patients as a delayed parasite clearance following artemisinin-based treatment, is conferred by non-synonymous mutations in the Kelch beta-propeller domain of the Plasmodium falciparum k13 ( pfk13 ) gene. Here, we carried out in vitro selection over a 1-year period on a West African P. falciparum strain isolated from Kolle (Mali) under a dose-escalating artemisinin regimen. After 18 cycles of sequential drug pressure, the selected parasites exhibited enhanced survival to dihydroartemisinin in the ring-stage survival assay (RSA 0-3h = 9.2%). Sanger and whole-genome sequence analyses identified the PfK13 P413A mutation, localized in the BTB/POZ domain, upstream of the propeller domain. This mutation was sufficient to confer in vitro artemisinin resistance when introduced into the PfK13 coding sequence of the parasite strain Dd2 by CRISPR/Cas9 gene editing. These results together with structural studies of the protein demonstrate that the propeller domain is not the sole in vitro mediator of PfK13-mediated artemisinin resistance, and highlight the importance of monitoring for mutations throughout PfK13.

Since its emergence in China at the end of 2019, SARS-CoV-2 has rapidly spread across the world to become a global public health emergency. Since then, the pandemic has evolved with the large worldwide emergence of new variants, such as the Alpha (B.1.1.7 variant), Beta (B.1.351 variant), and Gamma (P.1 variant), and some other under investigation such as the A.27 in France. Many studies are focusing on antibody neutralisation changes according to the spike mutations, but to date, little is known regarding their respective replication capacities. In this work, we demonstrate that the Alpha variant provides an earlier replication in vitro , on Vero E6 and A549 cells, than Beta, Gamma, A.27, and historical lineages. This earlier replication was associated with higher infectious titres in cell-culture supernatants, in line with the higher viral loads observed among Alpha-infected patients. Interestingly, Beta and Gamma variants presented similar kinetic and viral load than the other non-Alpha-tested variants., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Lebourgeois, Chenane, Houhou-Fidouh, Menidjel, Ferré, Collin, Benmalek, Coppée, Larrouy, Yazdanpanah, Timsit, Charpentier, Descamps and Visseaux.)

Sporozoite forms of the Plasmodium parasite, the causative agent of malaria, are transmitted by mosquitoes and first infect the liver for an initial round of replication before parasite proliferation in the blood. The molecular mechanisms involved during sporozoite invasion of hepatocytes remain poorly understood. Two receptors of the Hepatitis C virus (HCV), the tetraspanin CD81 and the scavenger receptor class B type 1 (SR-B1), play an important role during the entry of Plasmodium sporozoites into hepatocytes. In contrast to HCV entry, which requires both CD81 and SR-B1 together with additional host factors, CD81 and SR-B1 operate independently during malaria liver infection. Sporozoites from human-infecting P. falciparum and P. vivax rely respectively on CD81 or SR-B1. Rodent-infecting P. berghei can use SR-B1 to infect host cells as an alternative pathway to CD81, providing a tractable model to investigate the role of SR-B1 during Plasmodium liver infection. Here we show that mouse SR-B1 is less functional as compared to human SR-B1 during P. berghei infection. We took advantage of this functional difference to investigate the structural determinants of SR-B1 required for infection. Using a structure-guided strategy and chimeric mouse/human SR-B1 constructs, we could map the functional region of human SR-B1 within apical loops, suggesting that this region of the protein may play a crucial role for interaction of sporozoite ligands with host cells and thus the very first step of Plasmodium infection.

The continuous spread of antimalarial drug resistance is a threat to current chemotherapy efficacy. Therefore, characterizing the genetic diversity of drug resistance markers is needed to follow treatment effectiveness and further update control strategies. Here, we genotyped Plasmodium falciparum resistance gene markers associated with sulfadoxine-pyrimethamine (SP) and artemisinin-based combination therapy (ACT) in isolates from pregnant women in Ghana. The prevalence of the septuple IRN I- A/FG K GS/T pfdhfr / pfdhps haplotypes, including the pfdhps A581G and A613S/T mutations, was high at delivery among post-SP treatment isolates (18.2%) compared to those of first antenatal care (before initiation of intermittent preventive treatment of malaria in pregnancy with sulfadoxine-pyrimethamine [IPTp-SP]; 6.1%; P =  0.03). Regarding the pfk13 marker gene, two nonsynonymous mutations (N458D and A481C) were detected at positions previously related to artemisinin resistance in isolates from Southeast Asia. These mutations were predicted in silico to alter the stability of the pfk13 propeller-encoding domain. Overall, these findings highlight the need for intensified monitoring and surveillance of additional mutations associated with increased SP resistance as well as emergence of resistance against artemisinin derivatives., (Copyright © 2020 American Society for Microbiology.)

Mutations in the Plasmodium falciparum chloroquine resistance transporter (PfCRT) confer resistance to several antimalarial drugs such as chloroquine (CQ) or piperaquine (PPQ), a partner molecule in current artemisinin-based combination therapies. As a member of the Drug/Metabolite Transporter (DMT) superfamily, the vacuolar transporter PfCRT may translocate substrate molecule(s) across the membrane of the digestive vacuole (DV), a lysosome-like organelle. However, the physiological substrate(s), the transport mechanism and the functional regions of PfCRT remain to be fully characterized. Here, we hypothesized that identification of evolutionary conserved sites in a tertiary structural context could help locate putative functional regions of PfCRT. Hence, site-specific substitution rates were estimated over Plasmodium evolution at each amino acid sites, and the PfCRT tertiary structure was predicted in both inward-facing (open-to-vacuole) and occluded states through homology modeling using DMT template structures sharing <15% sequence identity with PfCRT. We found that the vacuolar-half and membrane-spanning domain (and especially the transmembrane helix 9) of PfCRT were more conserved, supporting that its physiological substrate is expelled out of the parasite DV. In the PfCRT occluded state, some evolutionary conserved sites, including positions related to drug resistance mutations, participate in a putative binding pocket located at the core of the PfCRT membrane-spanning domain. Through structural comparison with experimentally-characterized DMT transporters, we identified several conserved PfCRT amino acid sites located in this pocket as robust candidates for mediating substrate transport. Finally, in silico mutagenesis revealed that drug resistance mutations caused drastic changes in the electrostatic potential of the transporter vacuolar entry and pocket, facilitating the escape of protonated CQ and PPQ from the parasite DV.

Numerous mutations in the Plasmodium falciparum Kelch13 (K13) protein confer resistance to artemisinin derivatives, the current front-line antimalarial drugs. K13 is an essential protein that contains BTB and Kelch-repeat propeller (KREP) domains usually found in E3 ubiquitin ligase complexes that target substrate protein(s) for ubiquitin-dependent degradation. K13 is thought to bind substrate proteins, but its functional/interaction sites and the structural alterations associated with artemisinin resistance mutations remain unknown. Here, we screened for the most evolutionarily conserved sites in the protein structure of K13 as indicators of structural and/or functional constraints. We inferred structure-dependent substitution rates at each amino acid site of the highly conserved K13 protein during the evolution of Apicomplexa parasites. We found two solvent-exposed patches of extraordinarily conserved sites likely involved in protein-protein interactions, one in BTB and the other one in KREP. The conserved patch in K13 KREP overlaps with a shallow pocket that displays a differential electrostatic surface potential, relative to neighboring sites, and that is rich in serine and arginine residues. Comparative structural and evolutionary analyses revealed that these properties were also found in the functionally-validated shallow pocket of other KREPs including that of the cancer-related KEAP1 protein. Finally, molecular dynamics simulations carried out on PfK13 R539T and C580Y artemisinin resistance mutant structures revealed some local structural destabilization of KREP but not in its shallow pocket. These findings open new avenues of research on one of the most enigmatic malaria proteins with the utmost clinical importance.

Background: PfEMP1 is the major protein from parasitic origin involved in the pathophysiology of severe malaria, and PfEMP1 domain subtypes are associated with the infection outcome. In addition, PfEMP1 variability is endless and current publicly available protein repositories do not reflect the high diversity of the sequences of PfEMP1 proteins. The identification of PfEMP1 protein sequences expressed with samples remains challenging. The aim of our study is to identify the different PfEMP1 proteins variants expressed within patient samples, and therefore identify PfEMP1 proteins domains expressed by patients presenting uncomplicated malaria or severe malaria in malaria endemic setting in Cotonou, Benin., Methods: We performed a multi-omic approach to decipher PfEMP1 expression at the patient's level in different clinical settings. Using a combination of whole genome sequencing approach and RNA sequencing, we were able to identify new PfEMP1 sequences and created a new custom protein database. This database was used for protein identification in mass spectrometry analysis., Results: The differential expression analysis of RNAsequencing data shows an increased expression of the var domains transcripts DBLα1.7, DBLα1.1, DBLα2 and DBLβ12 in samples from patients suffering from Cerebral Malaria compared to Uncomplicated Malaria. Our approach allowed us to attribute PfEMP1 sequences to each sample and identify new peptides associated to PfEMP1 proteins in mass spectrometry., Conclusion: We highlighted the diversity of the PfEMP1 sequences from field sample compared to reference sequences repositories and confirmed the validity of our approach. These findings should contribute to further vaccine development strategies based on PfEMP1 proteins., Competing Interests: EG is a salaried employee of Inovarion (Paris, France https://www.inovarion.com/). The 3P5 group (DR, CB, FG) group received funding from Cancéropôle Île-de-France (https://www.canceropole-idf.fr/) to purchase the Orbitrap. There are no patents, products in development or marketed products associated with this research to declare. This does not alter our adherence to PLOS ONE policies on sharing data and materials. All other authors declare that no competing interests exist.

Background: Although trimethoprim-sulfamethoxazole is the more efficient drug for prophylactic and curative treatment of pneumocystosis, atovaquone is considered a second-line prophylactic treatment in immunocompromised patients. Variations in atovaquone absorption and mutant fungi selection after atovaquone exposure have been associated with atovaquone prophylactic failure. We report here a Pneumocystis jirovecii cytochrome b (cyt b) mutation (A144V) associated with such prophylactic failure during a pneumocystosis outbreak among heart transplant recipients., Methods: Analyses of clinical data, serum drug dosage, and molecular modeling of the P. jirovecii Rieske-cyt b complex were performed to investigate these prophylactic failures., Results: The cyt b A144V mutation was detected in all infected, heart transplant recipient patients exposed to atovaquone prophylaxis but in none of 11 other immunocompromised, infected control patients not treated with atovaquone. Serum atovaquone concentrations associated with these prophylactic failures were similar than those found in noninfected exposed control patients under a similar prophylactic regimen. Computational modeling of the P. jirovecii Rieske-cyt b complex and in silico mutagenesis indicated that the cyt b A144V mutation might alter the volume of the atovaquone-binding pocket, which could decrease atovaquone binding., Conclusions: These data suggest that the cyt b A144V mutation confers diminished sensitivity to atovaquone, resulting in spread of Pneumocystis pneumonia among heart transplant recipients submitted to atovaquone prophylaxis. Potential selection and interhuman transmission of resistant P. jirovecii strain during atovaquone prophylactic treatment has to be considered and could limit its extended large-scale use in immucompromised patients.