Publisher Correction: Genomic insights into the 2016-2017 cholera epidemic in Yemen.

In the HTML version of this Letter, the affiliations for authors Andrew S. Azman, Dhirendra Kumar and Thandavarayan Ramamurthy were inverted (the PDF and print versions of the Letter were correct); the affiliations have been corrected online.

Genomic insights into the 2016-2017 cholera epidemic in Yemen.

Yemen is currently experiencing, to our knowledge, the largest cholera epidemic in recent history. The first cases were declared in September 2016, and over 1.1 million cases and 2,300 deaths have since been reported1. Here we investigate the phylogenetic relationships, pathogenesis and determinants of antimicrobial resistance by sequencing the genomes of Vibrio cholerae isolates from the epidemic in Yemen and recent isolates from neighbouring regions. These 116 genomic sequences were placed within the phylogenetic context of a global collection of 1,087 isolates of the seventh pandemic V. cholerae serogroups O1 and O139 biotype El Tor2-4. We show that the isolates from Yemen that were collected during the two epidemiological waves of the epidemic1-the first between 28 September 2016 and 23 April 2017 (25,839 suspected cases) and the second beginning on 24 April 2017 (more than 1 million suspected cases)-are V. cholerae serotype Ogawa isolates from a single sublineage of the seventh pandemic V. cholerae O1 El Tor (7PET) lineage. Using genomic approaches, we link the epidemic in Yemen to global radiations of pandemic V. cholerae and show that this sublineage originated from South Asia and that it caused outbreaks in East Africa before appearing in Yemen. Furthermore, we show that the isolates from Yemen are susceptible to several antibiotics that are commonly used to treat cholera and to polymyxin B, resistance to which is used as a marker of the El Tor biotype.

Neisseria leonii sp. nov., isolated from the nose, lung, and liver of rabbits.

A taxogenomic study of three strains (3986T, 51.81, and JF 2415) isolated from rabbits between 1972 and 2000 led to the description of a new Neisseria species. The highest sequence similarity of the 16S rRNA gene was found to Neisseria animalis NCTC 10212T (96.7 %). The 16S rRNA gene similarity above 99 % and average nucleotide identity (ANI) values above 96 % among the strains, indicated that they belong to the same species. At the same time, the strains shared ANI values below 81 % and dDDH values below 24 % with all described Neisseria species. In the bac120 gene phylogenetic tree, the three strains clustered near Neisseria elongata and Neisseria bacilliformis in the Neisseria clade. However, the Neisseria clade is not monophyletic, and includes the type strains of Morococcus cerebrosus, Bergeriella denitrificans, Kingella potus, Uruburuella suis, and Uruburuella testudinis. Neisseria shayeganii clustered outside the clade with members of the genus Eikenella. Amino acid identity (AAI) values were calculated, and a threshold of 71 % was used to circumscribe the genus Neisseria. According to this proposed AAI threshold, strains 3986T, 51.81, and JF 2415 were placed within the genus Neisseria. The cells of the three strains were Gram-stain-negative diplococcobacilli and non-motile. Optimal growth on trypticase soy agar occurred at 37 °C and pH 8.5 in aerobic conditions. Notably, all strains exhibited indole production in the API-NH test, which is atypical for Neisseria and the family Neisseriaceae. The strains exhibited a common set of 68 peaks in their MALDI-TOF MS profiles, facilitating the swift and accurate identification of this species. Based on genotypic and phenotypic data, it is proposed that strains 3986T, 51.81, and JF 2415 represent a novel species within the genus Neisseria, for which the name Neisseria leonii sp. nov. is proposed (type strain 3986T=R726T=CIP 109994T=LMG 32907T).

Comparison of intra- and inter-host genetic diversity in rabies virus during experimental cross-species transmission.

The development of high-throughput genome sequencing enables accurate measurements of levels of sub-consensus intra-host virus genetic diversity and analysis of the role played by natural selection during cross-species transmission. We analysed the natural and experimental evolution of rabies virus (RABV), an important example of a virus that is able to make multiple host jumps. In particular, we (i) analyzed RABV evolution during experimental host switching with the goal of identifying possible genetic markers of host adaptation, (ii) compared the mutational changes observed during passage with those observed in natura, and (iii) determined whether the colonization of new hosts or tissues requires adaptive evolution in the virus. To address these aims, animal infection models (dog and fox) and primary cell culture models (embryo brain cells of dog and fox) were developed and viral variation was studied in detail through deep genome sequencing. Our analysis revealed a strong unidirectional host evolutionary effect, as dog-adapted rabies virus was able to replicate in fox and fox cells relatively easily, while dogs or neuronal dog cells were not easily susceptible to fox adapted-RABV. This suggests that dog RABV may be able to adapt to some hosts more easily than other host variants, or that when RABV switched from dogs to red foxes it lost its ability to adapt easily to other species. Although no difference in patterns of mutation variation between different host organs was observed, mutations were common following both in vitro and in vivo passage. However, only a small number of these mutations also appeared in natura, suggesting that adaptation during successful cross-species virus transmission is a complex, multifactorial evolutionary process.

Distinct lineages of Ebola virus in Guinea during the 2014 West African epidemic.

An epidemic of Ebola virus disease of unprecedented scale has been ongoing for more than a year in West Africa. As of 29 April 2015, there have been 26,277 reported total cases (of which 14,895 have been laboratory confirmed) resulting in 10,899 deaths. The source of the outbreak was traced to the prefecture of Guéckédou in the forested region of southeastern Guinea. The virus later spread to the capital, Conakry, and to the neighbouring countries of Sierra Leone, Liberia, Nigeria, Senegal and Mali. In March 2014, when the first cases were detected in Conakry, the Institut Pasteur of Dakar, Senegal, deployed a mobile laboratory in Donka hospital to provide diagnostic services to the greater Conakry urban area and other regions of Guinea. Through this process we sampled 85 Ebola viruses (EBOV) from patients infected from July to November 2014, and report their full genome sequences here. Phylogenetic analysis reveals the sustained transmission of three distinct viral lineages co-circulating in Guinea, including the urban setting of Conakry and its surroundings. One lineage is unique to Guinea and closely related to the earliest sampled viruses of the epidemic. A second lineage contains viruses probably reintroduced from neighbouring Sierra Leone on multiple occasions, while a third lineage later spread from Guinea to Mali. Each lineage is defined by multiple mutations, including non-synonymous changes in the virion protein 35 (VP35), glycoprotein (GP) and RNA-dependent RNA polymerase (L) proteins. The viral GP is characterized by a glycosylation site modification and mutations in the mucin-like domain that could modify the outer shape of the virion. These data illustrate the ongoing ability of EBOV to develop lineage-specific and potentially phenotypically important variation.

A molecular marker of artemisinin-resistant Plasmodium falciparum malaria.

Plasmodium falciparum resistance to artemisinin derivatives in southeast Asia threatens malaria control and elimination activities worldwide. To monitor the spread of artemisinin resistance, a molecular marker is urgently needed. Here, using whole-genome sequencing of an artemisinin-resistant parasite line from Africa and clinical parasite isolates from Cambodia, we associate mutations in the PF3D7_1343700 kelch propeller domain ('K13-propeller') with artemisinin resistance in vitro and in vivo. Mutant K13-propeller alleles cluster in Cambodian provinces where resistance is prevalent, and the increasing frequency of a dominant mutant K13-propeller allele correlates with the recent spread of resistance in western Cambodia. Strong correlations between the presence of a mutant allele, in vitro parasite survival rates and in vivo parasite clearance rates indicate that K13-propeller mutations are important determinants of artemisinin resistance. K13-propeller polymorphism constitutes a useful molecular marker for large-scale surveillance efforts to contain artemisinin resistance in the Greater Mekong Subregion and prevent its global spread.

Generating genomic platforms to study Candida albicans pathogenesis.

The advent of the genomic era has made elucidating gene function on a large scale a pressing challenge. ORFeome collections, whereby almost all ORFs of a given species are cloned and can be subsequently leveraged in multiple functional genomic approaches, represent valuable resources toward this endeavor. Here we provide novel, genome-scale tools for the study of Candida albicans, a commensal yeast that is also responsible for frequent superficial and disseminated infections in humans. We have generated an ORFeome collection composed of 5099 ORFs cloned in a Gateway™ donor vector, representing 83% of the currently annotated coding sequences of C. albicans. Sequencing data of the cloned ORFs are available in the CandidaOrfDB database at http://candidaorfeome.eu. We also engineered 49 expression vectors with a choice of promoters, tags and selection markers and demonstrated their applicability to the study of target ORFs transferred from the C. albicans ORFeome. In addition, the use of the ORFeome in the detection of protein-protein interaction was demonstrated. Mating-compatible strains as well as Gateway™-compatible two-hybrid vectors were engineered, validated and used in a proof of concept experiment. These unique and valuable resources should greatly facilitate future functional studies in C. albicans and the elucidation of mechanisms that underlie its pathogenicity.

Multiple major disease-associated clones of Legionella pneumophila have emerged recently and independently.

Legionella pneumophila is an environmental bacterium and the leading cause of Legionnaires' disease. Just five sequence types (ST), from more than 2000 currently described, cause nearly half of disease cases in northwest Europe. Here, we report the sequence and analyses of 364 L. pneumophila genomes, including 337 from the five disease-associated STs and 27 representative of the species diversity. Phylogenetic analyses revealed that the five STs have independent origins within a highly diverse species. The number of de novo mutations is extremely low with maximum pairwise single-nucleotide polymorphisms (SNPs) ranging from 19 (ST47) to 127 (ST1), which suggests emergences within the last century. Isolates sampled geographically far apart differ by only a few SNPs, demonstrating rapid dissemination. These five STs have been recombining recently, leading to a shared pool of allelic variants potentially contributing to their increased disease propensity. The oldest clone, ST1, has spread globally; between 1940 and 2000, four new clones have emerged in Europe, which show long-distance, rapid dispersal. That a large proportion of clinical cases is caused by recently emerged and internationally dispersed clones, linked by convergent evolution, is surprising for an environmental bacterium traditionally considered to be an opportunistic pathogen. To simultaneously explain recent emergence, rapid spread and increased disease association, we hypothesize that these STs have adapted to new man-made environmental niches, which may be linked by human infection and transmission.

Sequanix: a dynamic graphical interface for Snakemake workflows.

Summary: We designed a PyQt graphical user interface-Sequanix-aimed at democratizing the use of Snakemake pipelines in the NGS space and beyond. By default, Sequanix includes Sequana NGS pipelines (Snakemake format) (http://sequana.readthedocs.io), and is also capable of loading any external Snakemake pipeline. New users can easily, visually, edit configuration files of expert-validated pipelines and can interactively execute these production-ready workflows. Sequanix will be useful to both Snakemake developers in exposing their pipelines and to a wide audience of users. Availability and implementation: Source on http://github.com/sequana/sequana, bio-containers on http://bioconda.github.io and Singularity hub (http://singularity-hub.org). Contact: dimitri.desvillechabrol@pasteur.fr or thomas.cokelaer@pasteur.fr. Supplementary information: Supplementary data are available at Bioinformatics online.

Key experimental evidence of chromosomal DNA transfer among selected tuberculosis-causing mycobacteria.

Horizontal gene transfer (HGT) is a major driving force of bacterial diversification and evolution. For tuberculosis-causing mycobacteria, the impact of HGT in the emergence and distribution of dominant lineages remains a matter of debate. Here, by using fluorescence-assisted mating assays and whole genome sequencing, we present unique experimental evidence of chromosomal DNA transfer between tubercle bacilli of the early-branching Mycobacterium canettii clade. We found that the obtained recombinants had received multiple donor-derived DNA fragments in the size range of 100 bp to 118 kbp, fragments large enough to contain whole operons. Although the transfer frequency between M. canettii strains was low and no transfer could be observed among classical Mycobacterium tuberculosis complex (MTBC) strains, our study provides the proof of concept for genetic exchange in tubercle bacilli. This outstanding, now experimentally validated phenomenon presumably played a key role in the early evolution of the MTBC toward pathogenicity. Moreover, our findings also provide important information for the risk evaluation of potential transfer of drug resistance and fitness mutations among clinically relevant mycobacterial strains.

Large-Scale Phylogenomic Analysis Reveals the Complex Evolutionary History of Rabies Virus in Multiple Carnivore Hosts.

The natural evolution of rabies virus (RABV) provides a potent example of multiple host shifts and an important opportunity to determine the mechanisms that underpin viral emergence. Using 321 genome sequences spanning an unprecedented diversity of RABV, we compared evolutionary rates and selection pressures in viruses sampled from multiple primary host shifts that occurred on various continents. Two major phylogenetic groups, bat-related RABV and dog-related RABV, experiencing markedly different evolutionary dynamics were identified. While no correlation between time and genetic divergence was found in bat-related RABV, the evolution of dog-related RABV followed a generally clock-like structure, although with a relatively low evolutionary rate. Subsequent molecular clock dating indicated that dog-related RABV likely underwent a rapid global spread following the intensification of intercontinental trade starting in the 15th century. Strikingly, although dog RABV has jumped to various wildlife species from the order Carnivora, we found no clear evidence that these host-jumping events involved adaptive evolution, with RABV instead characterized by strong purifying selection, suggesting that ecological processes also play an important role in shaping patterns of emergence. However, specific amino acid changes were associated with the parallel emergence of RABV in ferret-badgers in Asia, and some host shifts were associated with increases in evolutionary rate, particularly in the ferret-badger and mongoose, implying that changes in host species can have important impacts on evolutionary dynamics.

DNA minicircles clarify the specific role of DNA structure on retroviral integration.

Chromatin regulates the selectivity of retroviral integration into the genome of infected cells. At the nucleosome level, both histones and DNA structure are involved in this regulation. We propose a strategy that allows to specifically study a single factor: the DNA distortion induced by the nucleosome. This strategy relies on mimicking this distortion using DNA minicircles (MCs) having a fixed rotational orientation of DNA curvature, coupled with atomic-resolution modeling. Contrasting MCs with linear DNA fragments having identical sequences enabled us to analyze the impact of DNA distortion on the efficiency and selectivity of integration. We observed a global enhancement of HIV-1 integration in MCs and an enrichment of integration sites in the outward-facing DNA major grooves. Both of these changes are favored by LEDGF/p75, revealing a new, histone-independent role of this integration cofactor. PFV integration is also enhanced in MCs, but is not associated with a periodic redistribution of integration sites, thus highlighting its distinct catalytic properties. MCs help to separate the roles of target DNA structure, histone modifications and integrase (IN) cofactors during retroviral integration and to reveal IN-specific regulation mechanisms.

Insights into the smooth-to-rough transitioning in Mycobacterium bolletii unravels a functional Tyr residue conserved in all mycobacterial MmpL family members.

In mycobacteria, MmpL proteins represent key components that participate in the biosynthesis of the complex cell envelope. Whole genome analysis of a spontaneous rough morphotype variant of Mycobacterium abscessus subsp. bolletii identified a conserved tyrosine that is crucial for the function of MmpL family proteins. Isogenic smooth (S) and rough (R) variants differed by a single mutation linked to a Y842H substitution in MmpL4a. This mutation caused a deficiency in glycopeptidolipid production/transport in the R variant and a gain in the capacity to produce cords in vitro. In zebrafish, increased virulence of the M. bolletii R variant over the parental S strain was found, involving massive production of serpentine cords, abscess formation and rapid larval death. Importantly, this finding allowed us to demonstrate an essential role of Tyr842 in several different MmpL proteins, including Mycobacterium tuberculosis MmpL3. Structural homology models of MmpL4a and MmpL3 identified two additional critical residues located in the transmembrane regions TM10 and TM4 that are facing each other. We propose that these central residues are part of the proton-motive force that supplies the energy for substrate transport. Hence, we provide important insights into mechanistic/structural aspects of MmpL proteins as lipid transporters and virulence determinants in mycobacteria.

A new piperidinol derivative targeting mycolic acid transport in Mycobacterium abscessus.

The natural resistance of Mycobacterium abscessus to most commonly available antibiotics seriously limits chemotherapeutic treatment options, which is particularly challenging for cystic fibrosis patients infected with this rapid-growing mycobacterium. New drugs with novel molecular targets are urgently needed against this emerging pathogen. However, the discovery of such new chemotypes has not been appropriately performed. Here, we demonstrate the utility of a phenotypic screen for bactericidal compounds against M. abscessus using a library of compounds previously validated for activity against M. tuberculosis. We identified a new piperidinol-based molecule, PIPD1, exhibiting potent activity against clinical M. abscessus strains in vitro and in infected macrophages. Treatment of infected zebrafish with PIPD1 correlated with increased embryo survival and decreased bacterial burden. Whole genome analysis of M. abscessus strains resistant to PIPD1 identified several mutations in MAB_4508, encoding a protein homologous to MmpL3. Biochemical analyses demonstrated that while de novo mycolic acid synthesis was unaffected, PIPD1 strongly inhibited the transport of trehalose monomycolate, thereby abrogating mycolylation of arabinogalactan. Mapping the mutations conferring resistance to PIPD1 on a MAB_4508 tridimensional homology model defined a potential PIPD1-binding pocket. Our data emphasize a yet unexploited chemical structure class against M. abscessus infections with promising translational development possibilities.

Resistance to Thiacetazone Derivatives Active against Mycobacterium abscessus Involves Mutations in the MmpL5 Transcriptional Repressor MAB_4384.

Available chemotherapeutic options are very limited against Mycobacterium abscessus, which imparts a particular challenge in the treatment of cystic fibrosis (CF) patients infected with this rapidly growing mycobacterium. New drugs are urgently needed against this emerging pathogen, but the discovery of active chemotypes has not been performed intensively. Interestingly, however, the repurposing of thiacetazone (TAC), a drug once used to treat tuberculosis, has increased following the deciphering of its mechanism of action and the detection of significantly more potent analogues. We therefore report studies performed on a library of 38 TAC-related derivatives previously evaluated for their antitubercular activity. Several compounds, including D6, D15, and D17, were found to exhibit potent activity in vitro against M. abscessus, Mycobacterium massiliense, and Mycobacterium bolletii clinical isolates from CF and non-CF patients. Similar to TAC in Mycobacterium tuberculosis, the three analogues act as prodrugs in M. abscessus, requiring bioactivation by the EthA enzyme, MAB_0985. Importantly, mutations in the transcriptional TetR repressor MAB_4384, with concomitant upregulation of the divergently oriented adjacent genes encoding an MmpS5/MmpL5 efflux pump system, accounted for high cross-resistance levels among all three compounds. Overall, this study uncovered a new mechanism of drug resistance in M. abscessus and demonstrated that simple structural optimization of the TAC scaffold can lead to the development of new drug candidates against M. abscessus infections.

Molecular bases for multidrug resistance in Yersinia pseudotuberculosis.

The enteropathogen Yersinia pseudotuberculosis causes gastrointestinal infections in humans. Although this species is usually susceptible to antibiotics active against Gram-negative bacteria, we identified three multidrug resistant (MDR) strains of Y. pseudotuberculosis that were isolated from the environment in Russia and from a patient in France. The resistance traits of the two Russian isolates were transferable at high frequencies (≈2×10-1/donor CFU) to Y. pseudotuberculosis. In contrast no transfer of the antibiotic resistances carried by the French strain was observed. Sequencing of the plasmid extracts of the Y. pseudotuberculosis transconjugants for the Russian isolates revealed the presence of conjugative plasmids of the IncN group that carried genes conferring resistance to four to six classes of antibiotics. The French strain harbored a large MDR plasmid of the IncHI2 group that carried resistance genes to six families of antibiotics, and contained a truncated set of transfer genes, accounting for the lack of plasmid transfer. All three Y. pseudotuberculosis plasmids were homologous to MDR plasmids found in various enterobacteria. A phylogenetic analysis showed that the two Russian strain plasmids were closely related to each other and were more distant from the French plasmid. To the best of our knowledge, this is the first molecular characterization of MDR plasmids in Y. pseudotuberculosis. Due to the propensity of this species to acquire exogenous plasmids, the risk of emergence of new MDR Y. pseudotuberculosis isolates should be seriously taken into consideration.

Parallel independent evolution of pathogenicity within the genus Yersinia.

The genus Yersinia has been used as a model system to study pathogen evolution. Using whole-genome sequencing of all Yersinia species, we delineate the gene complement of the whole genus and define patterns of virulence evolution. Multiple distinct ecological specializations appear to have split pathogenic strains from environmental, nonpathogenic lineages. This split demonstrates that contrary to hypotheses that all pathogenic Yersinia species share a recent common pathogenic ancestor, they have evolved independently but followed parallel evolutionary paths in acquiring the same virulence determinants as well as becoming progressively more limited metabolically. Shared virulence determinants are limited to the virulence plasmid pYV and the attachment invasion locus ail. These acquisitions, together with genomic variations in metabolic pathways, have resulted in the parallel emergence of related pathogens displaying an increasingly specialized lifestyle with a spectrum of virulence potential, an emerging theme in the evolution of other important human pathogens.

Genome-wide mosaicism within Mycobacterium abscessus: evolutionary and epidemiological implications.

BACKGROUND: In mycobacteria, conjugation differs from the canonical Hfr model, but is still poorly understood. Here, we quantified this evolutionary processe in a natural mycobacterial population, taking advantage of a large clinical strain collection of the emerging pathogen Mycobacterium abscessus (MAB). RESULTS: Multilocus sequence typing confirmed the existence of three M. abscessus subspecies, and unravelled extensive allelic exchange between them. Furthermore, an asymmetrical gene flow occurring between these main lineages was detected, resulting in highly admixed strains. Intriguingly, these mosaic strains were significantly associated with cystic fibrosis patients with lung infections or chronic colonization. Genome sequencing of those hybrid strains confirmed that half of their genomic content was remodelled in large genomic blocks, leading to original tri-modal 'patchwork' architecture. One of these hybrid strains acquired a locus conferring inducible macrolide resistance, and a large genomic insertion from a slowly growing pathogenic mycobacteria, suggesting an adaptive gene transfer. This atypical genomic architecture of the highly recombinogenic strains is consistent with the distributive conjugal transfer (DCT) observed in M. smegmatis. Intriguingly, no known DCT function was found in M. abscessus chromosome, however, a p-RAW-like genetic element was detected in one of the highly admixed strains. CONCLUSION: Taken together, our results strongly suggest that MAB evolution is sporadically punctuated by dramatic genome wide remodelling events. These findings might have far reaching epidemiological consequences for emerging mycobacterial pathogens survey in the context of increasing numbers of rapidly growing mycobacteria and M. tuberculosis co-infections.

Evidence of lymphocytic choriomeningitis virus (LCMV) in domestic mice in Gabon: risk of emergence of LCMV encephalitis in Central Africa.

Lymphocytic choriomeningitis virus (LCMV) can cause acute fatal disease on all continents but was never detected in Africa. We report the first detection of LCMV RNA in a common European house mouse (Mus musculus domesticus) in Africa. Phylogenetic analyses show a close relationship with North American strains. These findings suggest that there is a risk of the appearance of LCMV acute encephalitis cases. This is a perfect example of virus dissemination by its natural host that may have dramatic public health consequences.

Molecular characterization of Orf virus in goats in Gabon, Central Africa.

BACKGROUND: Orf or contagious ecthyma is a zoonotic viral infection with a potential serious health threat for the small ruminants industry as well as humans. It is currently emerging in new territories. RESULTS: Eight suspected clinical cases of pustular dermatitis in goats occurred in the rural area of Tebe, in south-eastern Gabon, in January 2013. The orf virus (ORFV) was detected by high-throughput sequencing on sera, buccal swabs and scab pool samples. It was confirmed in six out of eight sick goats by using specific PCR targeting the major envelope protein (B2L) and the orf virus interferon resistance (VIR) genes. Phylogenetic analysis revealed that the Gabonese strain and South Korean strains evolved from a common ancestor, suggesting an Asian origin of the ORFV' Gabonese strain. CONCLUSIONS: This study provides the molecular detection of the ORFV strain involved in the cases of pustular dermatitis in goats and highlights its circulation in Gabon.