Search

The leucine rich repeat gene APL1 is a key component of immunity to Plasmodium and other microbial pathogens in Anopheles mosquitoes. In the malaria vector Anopheles funestus the APL1 gene has four paralogues which occur along the same chromosome arm. We show that APL1 has exceptional levels of non-synonymous polymorphism across the range of An. funestus with an average πn of 0.027 versus a genome-wide average of 0.002, and πn (and πs) is consistently high in populations across Africa. The pattern of APL1 diversity was consistent between independent pooled-template and target-enrichment datasets, however no link between APL1 diversity and insecticide-resistance was observed with the phenotyped target-enrichment dataset. Two further innate immunity genes of the gambicin anti-microbial peptide family had πn/πs ratios greater than one, possibly driven by either positive or balancing selection. Cecropin antimicrobial peptides were expressed much more highly than other anti-microbial peptide genes, an observation discordant with current models of anti-microbial peptide activity. The observed APL1 diversity likely results from gene conversion between paralogs, as evidenced by shared polymorphisms, overlapping read mappings, and recombination events among paralogues. Gene conversion at APL1 versus alternative explanations is concordant with similarly elevated diversity in APL1 and TEP1 loci in An. gambiae. In contrast, the more closely related An. stephensi which also encodes a single-copy of APL1 does not show this elevated diversity. We hypothesise that a more open chromatin formation at the APL1 locus due to higher gene expression than its paralogues enhances gene conversion, and therefore increased polymorphism, at APL1.

Only female insects transmit diseases such as malaria, dengue and Zika; therefore, control methods that bias the sex ratio of insect offspring have long been sought. Genetic elements such as sex-chromosome drives can distort sex ratios to produce unisex populations that eventually collapse, but the underlying molecular mechanisms are unknown. We report a male-biased sex-distorter gene drive (SDGD) in the human malaria vector Anopheles gambiae. We induced super-Mendelian inheritance of the X-chromosome-shredding I-PpoI nuclease by coupling this to a CRISPR-based gene drive inserted into a conserved sequence of the doublesex (dsx) gene. In modeling of invasion dynamics, SDGD was predicted to have a quicker impact on female mosquito populations than previously developed gene drives targeting female fertility. The SDGD at the dsx locus led to a male-only population from a 2.5% starting allelic frequency in 10–14 generations, with population collapse and no selection for resistance. Our results support the use of SDGD for malaria vector control., A sex-distorter gene drive causes population collapse in the malaria mosquito.

Rationale: Longitudinal modeling of current wheezing identified similar phenotypes, but their characteristics often differ between studies. Objectives: We propose that a more comprehensive description of wheeze may better describe trajectories than binary information on the presence/absence of wheezing. Methods: We derived six multidimensional variables of wheezing spells from birth to adolescence (including duration, temporal sequencing, and the extent of persistence/recurrence). We applied partition-around-medoids clustering on these variables to derive phenotypes in five birth cohorts. We investigated within- and between-phenotype differences compared with binary latent class analysis models and ascertained associations of these phenotypes with asthma and lung function and with polymorphisms in asthma loci 17q12-21 and CDHR3 (cadherin-related family member 3). Measurements and Main Results: Analysis among 7,719 participants with complete data identified five spell-based wheeze phenotypes with a high degree of certainty: never (54.1%), early-transient (ETW) (23.7%), late-onset (LOW) (6.9%), persistent (PEW) (8.3%), and a novel phenotype, intermittent wheeze (INT) (6.9%). FEV 1/FVC was lower in PEW and INT compared with ETW and LOW and declined from age 8 years to adulthood in INT. 17q12-21 and CDHR3 polymorphisms were associated with higher odds of PEW and INT, but not ETW or LOW. Latent class analysis- and spell-based phenotypes appeared similar, but within-phenotype individual trajectories and phenotype allocation differed substantially. The spell-based approach was much more robust in dealing with missing data, and the derived clusters were more stable and internally homogeneous. Conclusions: Modeling of spell variables identified a novel intermittent wheeze phenotype associated with lung function decline to early adulthood. Using multidimensional spell variables may better capture wheeze development and provide a more robust input for phenotype derivation.

Kartagener’s syndrome is a subgroup of primary ciliary dyskinesia (PCD), a genetically heterogeneous condition characterised by sinusitis, bronchiectasis, and situs in versus. Genetic testing has importance for their diagnosis. Here, we report a chinese patient with Kartagener’s syndrome. Transthoracic echocardiography showed severely elevated right ventricular systolic pressure. Right heart catheterisation demonstrated a pre-capillary pulmonary hypertension. Whole-exome sequencing indicated that she had a novel homozygous nonsense mutation, c.2845C > T, p.Gln949*, in exon 18 of CCDC40 and a heterozygotic mutation, c.73G > A, p.Ala25Thr, in exon 1 of DNAH11. She was diagnosed as Kartagener’s syndrome with pulmonary hypertension. Her symptoms improved significantly by treatment of antibiotics, expectorant drugs, bronchodilators, and oxygen therapy treatment. Our findings extend the mutation spectrum of CCDC40 gene related Kartagener’s syndrome, which is very important for gene diagnosis of the disease.

Background Cameroon is considering the implementation of indoor residual spraying (IRS) as a complementary measure to control malaria in the context of high pyrethroid resistance in major malaria vectors. Non-pyrethroid insecticide classes such as organophosphates and carbamates may be utilized in IRS due to widespread pyrethroid resistance. However, the success of this strategy depends on good knowledge of the resistance status of malaria vectors to carbamates and organophosphates. Here, we assessed the susceptibility profile of Anopheles gambiae sensu lato with respect to carbamates and organophosphate and the distribution of the molecular mechanism underlying resistance to these insecticides. Methods Anopheles gambiae s.l. mosquitoes were collected from nine settings across the country and bio-assayed with bendiocarb, propoxur and pirimiphos-methyl. The Ace-1 target-site G119S mutation was genotyped using a TaqMan assay. To investigate the polymorphism in the Ace-1 gene, a region of 924 base pairs in a sequence of the gene was amplified from both live and dead females of An. gambiae exposed to bendiocarb. Results Pirimiphos-methyl induced full mortality in An. gambiae s.l. from all study sites, whereas for carbamates, resistance was observed in four localities, with the lowest mortality rate recorded in Mangoum (17.78 ± 5.02% for bendiocarb and 18.61 ± 3.86% for propoxur) in the southern part of Cameroon. Anopheles coluzzii was found to be the predominant species in the northern tropical part of the country where it is sympatric with Anopheles arabiensis. In the localities situated in southern equatorial regions, this species was predominant in urban settings, while An. gambiae was the most abundant species in rural areas. The G119S Ace-1 target-site mutation was detected only in An. gambiae and only in the sites located in southern Cameroon. Phylogenetic analyses showed a clustering according to the phenotype. Conclusion The occurrence of the Ace-1 target-site substitution G119S in An. gambiae s.l. populations highlights the challenge associated with the impending deployment of IRS in Cameroon using carbamates or organophosphates. It is therefore important to think about a resistance management plan including the use of other insecticide classes such as neonicotinoids or pyrrole to guarantee the implementation of IRS in Cameroon. Graphical Abstract

Mobile genetic elements (MGEs) are often associated with antimicrobial resistance genes (ARGs). They are responsible for intracellular transposition between different replicons and intercellular conjugation, and are therefore important agents of ARG dissemination. Detection and characterisation of functional MGEs, especially in clinical isolates, would increase our understanding of the underlying pathways of transposition and recombination, and allow us to determine interventional strategies to interrupt this process.\ud Entrapment vectors can be used to capture active MGEs as they contain a positive selection genetic system conferring a selectable phenotype upon insertion of an MGE within certain regions of that system. Previously, we developed the pBACpAK entrapment vector that results in a tetracycline-resistant phenotype when MGEs translocate and disrupt the cI repressor gene. We have previously used pBACpAK to capture MGEs in clinical Escherichia coli isolates following transformation with pBACpAK.\ud In this study, we aimed to extend the utilisation of pBACpAK to other bacterial species. We utilised an MGE free recipient E. coli strain containing pBACpAK to capture MGEs on conjugative, ARG containing plasmids following conjugation from clinical Enterobacteriaceae donors.\ud Following the conjugative transfer of multiple conjugative plasmids, and screening for tetracycline resistance in these transconjugants, we captured several IS elements and novel transposons (Tn7350 and Tn7351), we detected the de novo formation of novel putative composite transposons where the pBACpAK located tet(A) is flanked by ISKpn25 from the transferred conjugative plasmid and we also detected the ISKpn14 mediated integration of an entire 119kb, blaNDM-1 containing conjugative plasmid from Klebsiella pneumoniae.\ud Importance\ud By analysing transposition activity within our MGE free recipient, we can gain insights into the interaction and evolution of multidrug resistance conferring MGEs following conjugation, including the movement of multiple ISs, the formation of composite transposons and the cointegration and, or recombination between different replicons in the same cell. This combination of recipient and entrapment vector will allow for fine-scale experimental studies into factors affecting intracellular transposition and MGE formation in, and from, ARG encoding MGEs from multiple species of clinically relevant Enterobacteriaceae.

BackgroundMonitoring of drug resistance in Plasmodium populations is crucial for malaria control. This has primarily been performed in humans and rarely in mosquitoes where parasites genetic recombination occurs. Here, we characterized the Plasmodium spp populations in wild Anopheles vectors by analyzing the genetic diversity of the P. falciparum kelch13 and mdr1 gene fragments implicated in artemisinin and partner drug resistance across Cameroon in three major malaria vectors.MethodsAnopheles mosquitoes were collected across nine localities in Cameroon and dissected into the head/thorax (H/T) and abdomen (Abd) after species identification. A TaqMan assay was performed to detect Plasmodium infection. Fragments of the Kelch 13 and mdr1 genes were amplified in P. falciparum positive samples and directly sequenced to assess their drug resistance polymorphisms and genetic diversity profile.ResultsThe study revealed a high Plasmodium infection rate in the major Anopheles vectors across Cameroon. Notably, An. funestus vector recorded the highest sporozoite (8.02%) and oocyst (14.41%) infection rates. A high P. falciparum sporozoite rate (80.08%) alongside epidemiological signatures of significant P. malariae (15.94%) circulation were recorded in these vectors. Low genetic diversity with six (A578S, R575I, G450R, L663L, G453D, N458D) and eight (H53H, V62L, V77E, N86Y, G102G, L132I, H143H, Y184F) point mutations were observed in the k13 and mdr1 backbones respectively. Remarkably, the R575I (4.44%) k13 and Y184F (64.2%) mdr1 mutations were the predominant variants in the P. falciparum populations.ConclusionThe emerging signal of the R575I polymorphism in the Pfk13 propeller backbone entails the regular surveillance of molecular markers to inform evidence-based policy decisions. Moreover, the high frequency of the 86N184F allele highlights concerns on the plausible decline in efficacy of artemisinin-combination therapies (ACTs); further implying that parasite genotyping from mosquitoes can provide a more relevant scale for quantifying resistance epidemiology in the field.

Bloodstream infections (BSIs) caused by carbapenem-resistant Klebsiella pneumoniae (CRKP) are potentially life-threatening and an urgent threat to public health. The present study aims to clarify the characteristics of carbapenemase-encoding and virulent plasmids, and their interactions with the host bacterium. A total of 425 Kp isolates were collected from the blood of BSI patients from nine Chinese hospitals, between 2005 and 2019. Integrated epidemiological and genomic data showed that ST11 and ST307 Kp isolates were associated with nosocomial outbreak and transmission. Comparative analysis of 147 Kp genomes and 39 completely assembled chromosomes revealed extensive interruption of acrR by ISKpn26 in all Kp carbapenemase-2 (KPC-2)-producing ST11 Kp isolates, leading to activation of the AcrAB-Tolc multidrug efflux pump and a subsequent reduction in susceptibility to the last-resort antibiotic tigecycline and six other antibiotics. We described 29 KPC-2 plasmids showing diverse structures, two virulence plasmids in two KPC-2-producing Kp, and two novel multidrug-resistant (MDR)-virulent plasmids. This study revealed a multifactorial impact of KPC-2 plasmid on Kp, which may be associated with nosocomial dissemination of MDR isolates.

Genomic epidemiology is a tool for tracing transmission of pathogens based on whole-genome sequencing. We introduce the mGEMS pipeline for genomic epidemiology with plate sweeps representing mixed samples of a target pathogen, opening the possibility to sequence all colonies on selective plates with a single DNA extraction and sequencing step. The pipeline includes the novel mGEMS read binner for probabilistic assignments of sequencing reads, and the scalable pseudoaligner Themisto. We demonstrate the effectiveness of our approach using closely related samples in a nosocomial setting, obtaining results that are comparable to those based on single-colony picks. Our results lend firm support to more widespread consideration of genomic epidemiology with mixed infection samples.

Emergent hypervirulent Klebsiella pneumoniae has been responsible for severe diseases, representing a serious threat to public health. We report the whole-genome sequencing of a novel ST3994-K2 clone, a single locus variant of ST86 K2, which is considered a worrying hypervirulent clone that emerged in several parts of the world. The strain K. pneumonia Kpi144 was isolated in 2013 from a blood culture of a 69-year-old male patient admitted to a tertiary hospital in Teresina, state of Piauí, northeastern Brazil. The strain was susceptible to 41 antibiotics tested, presented hypermucoviscous phenotype, and a virulent behavior was observed in the Galleria mellonella infection model. Moreover, the virulome showed several virulence genes. To the best of our knowledge, this is the first worldwide report of a novel ST3994-K2 K. pneumoniae clone, an SLV of ST86 K2, which is considered a worrying virulent clone that has emerged in several parts of the world, including South America and Brazil.

The sequence diversity of natural and laboratory populations of Brugia pahangi and Brugia malayi was assessed with Illumina resequencing followed by mapping in order to identify single nucleotide variants and insertions/deletions. In natural and laboratory Brugia populations, there is a lack of sequence diversity on chromosome X relative to the autosomes (πX/πA = 0.2), which is lower than the expected (πX/πA = 0.75). A reduction in diversity is also observed in other filarial nematodes with neo-X chromosome fusions in the genera Onchocerca and Wuchereria, but not those without neo-X chromosome fusions in the genera Loa and Dirofilaria. In the species with neo-X chromosome fusions, chromosome X is abnormally large, containing a third of the genetic material such that a sizable portion of the genome is lacking sequence diversity. Such profound differences in genetic diversity can be consequential, having been associated with drug resistance and adaptability, with the potential to affect filarial eradication., Author summary Almost a billion people receive >7.7 billion doses of treatment aimed at eliminating lymphatic filariasis, which is caused by three filarial nematodes: Wuchereria bancrofti, Brugia malayi, and Brugia timori. Drug resistance and adaptation are both associated with pathogen success as well as higher levels of genetic diversity. In an examination of genetic diversity in Brugia malayi and Brugia pahangi, we observed a lack of genetic diversity over a third of the genome that is found on chromosome X. These species have neo-X chromosomes where a chromosome X fused with an autosome. Using publicly-available published data, the other filarial nematodes of greatest human significance are also found to have a similar lack of genetic diversity on their neo-X chromosomes. The two filarial nematodes with publicly-available data that lack neo-X chromosomes did not have this lack of genetic diversity. This lack of sequence diversity in B. malayi, W. bancrofti, and O. volvulus could have profound effects on all traits encoded on chromosome X.

CRISPR-based homing gene drives can be designed to disrupt essential genes whilst biasing their own inheritance, leading to suppression of mosquito populations in the laboratory. This class of gene drives relies on CRISPR-Cas9 cleavage of a target sequence and copying (‘homing’) therein of the gene drive element from the homologous chromosome. However, target site mutations that are resistant to cleavage yet maintain the function of the essential gene are expected to be strongly selected for. Targeting functionally constrained regions where mutations are not easily tolerated should lower the probability of resistance. Evolutionary conservation at the sequence level is often a reliable indicator of functional constraint, though the actual level of underlying constraint between one conserved sequence and another can vary widely. Here we generated a novel adult lethal gene drive (ALGD) in the malaria vector Anopheles gambiae, targeting an ultra-conserved target site in a haplosufficient essential gene (AGAP029113) required during mosquito development, which fulfils many of the criteria for the target of a population suppression gene drive. We then designed a selection regime to experimentally assess the likelihood of generation and subsequent selection of gene drive resistant mutations at its target site. We simulated, in a caged population, a scenario where the gene drive was approaching fixation, where selection for resistance is expected to be strongest. Continuous sampling of the target locus revealed that a single, restorative, in-frame nucleotide substitution was selected. Our findings show that ultra-conservation alone need not be predictive of a site that is refractory to target site resistance. Our strategy to evaluate resistance in vivo could help to validate candidate gene drive targets for their resilience to resistance and help to improve predictions of the invasion dynamics of gene drives in field populations., Author summary Gene drives have the potential to be applied as a novel control strategy of disease-transmitting mosquitoes, by spreading genetic traits that suppress or modify the target population. Many gene drive elements work by recognising and cutting a specific target sequence in the mosquito genome and copying themselves into that target sequence allowing the gene drive to increase in frequency in the population. Like other mosquito control interventions, efficacy will greatly depend on minimising the development of resistance to the gene drive mechanism—most likely via a change in the target sequence that prevents further cutting. One strategy to reduce resistance is to target sequences that are highly conserved, which implies that changes cannot easily be tolerated. We developed a strategy that simulates high selection pressure, under which resistance is most likely to emerge, and therefore provides a stringent test of its propensity to arise. Unlike previous results with another gene drive, we recovered a resistant allele within a few generations of gene drive exposure and at high frequency. Our results show that conserved sequences can vary hugely in ability to tolerate mutations and highlights the need to functionally validate future candidate gene drive target sites for their robustness to resistance.

Background. Regulatory clinical trials are required to ensure the continued supply and deployment of effective antimalarial drugs. Patient follow-up in such trials typically lasts several weeks as the drugs have long half-lives and new infections often occur during this period. "Molecular correction" is therefore used to distinguish drug failures from new infections. The current WHO-recommend method for molecular correction uses length-polymorphic alleles at highly diverse loci but is inherently poor at detecting low density clones in polyclonal infections. This likely leads to substantial underestimates of failure rates, delaying the replacement of failing drugs with potentially lethal consequences. Deep sequenced amplicons (AmpSeq) substantially increase the detectability of low-density clones and may offer a new "gold standard" for molecular correction. Methods. Pharmacological simulation of clinical trials was used to evaluate the suitability of AmpSeq for molecular correction. We investigated the impact of factors such as the number of amplicon loci analysed, the informatics criteria used to distinguish genotyping 'noise' from real low density signals, the local epidemiology of malaria transmission, and the potential impact of genetic signals from gametocytes. Results. AmpSeq greatly improved molecular correction and provided accurate drug failure rate estimates. The use of 3 to 5 amplicons was sufficient, and simple, non-statistical, criteria could be used to classify recurrent infections as drug failures or new infections. Conclusions. These results suggest AmpSeq is strongly placed to become the new standard for molecular correction in regulatory trials, with its potential extension into routine surveillance once the requisite technical support becomes established.

Transposable elements in prokaryotes are found in many forms and therefore a robust nomenclature system is needed in order to allow researchers to describe and search for them in publications and databases. Here we provide an update on The Transposon Registry which allocates numbers to any prokaryotic transposable element. Additionally, we present the completion of registry records for all transposons assigned Tn numbers from Tn1 onwards where sequence data or publications exist.

The Anopheles gambiae complex consists of multiple morphologically indistinguishable mosquito species including the most important vectors of the malaria parasite Plasmodium falciparum in sub-Saharan Africa. Nine cryptic species have been described so far within the complex. The ecological, immunological and reproductive differences among these species will critically impact population responses to disease control strategies and environmental changes. Here, we examine whole-genome sequencing data from a longitudinal study of putative A. coluzzii in western Burkina Faso. Surprisingly, many specimens are genetically divergent from A. coluzzii and all other Anopheles species and represent a new taxon, here designated Anopheles TENGRELA (AT). Population genetic analysis suggests that the cryptic GOUNDRY subgroup, previously collected as larvae in central Burkina Faso, represents an admixed population descended from both A. coluzzii and AT. AT harbours low nucleotide diversity except for the 2La inversion polymorphism which is maintained by overdominance. It shows numerous fixed differences with A. coluzzii concentrated in several regions reflecting selective sweeps, but the two taxa are identical at standard diagnostic loci used for taxon identification, and thus, AT may often go unnoticed. We present an amplicon-based genotyping assay for identifying AT which could be usefully applied to numerous existing samples. Misidentified cryptic taxa could seriously confound ongoing studies of Anopheles ecology and evolution in western Africa, including phenotypic and genotypic surveys of insecticide resistance. Reproductive barriers between cryptic species may also complicate novel vector control efforts, for example gene drives, and hinder predictions about evolutionary dynamics of Anopheles and Plasmodium. [Abstract copyright: © 2020 John Wiley & Sons Ltd.]

Homing-based gene drives use a germline source of nuclease to copy themselves at specific target sites in a genome and bias their inheritance. Such gene drives can be designed to spread and deliberately suppress populations of malaria mosquitoes by impairing female fertility. However, strong unintended fitness costs of the drive and a propensity to generate resistant mutations can limit a gene drive’s potential to spread. Alternative germline regulatory sequences in the drive element confer improved fecundity of carrier individuals and reduced propensity for target site resistance. This is explained by reduced rates of end-joining repair of DNA breaks from parentally deposited nuclease in the embryo, which can produce heritable mutations that reduce gene drive penetrance. We tracked the generation and selection of resistant mutations over the course of a gene drive invasion of a population. Improved gene drives show faster invasion dynamics, increased suppressive effect and later onset of target site resistance. Our results show that regulation of nuclease expression is as important as the choice of target site when developing a robust homing-based gene drive for population suppression., Author summary Gene drives are selfish genetic elements that are able to drastically bias their own inheritance. They can rapidly invade populations, even starting from a very low frequency. Recent advances have allowed the engineering of gene drives deliberately designed to spread genetic traits of choice into populations of malaria-transmitting mosquito species–for example traits that impair a mosquito’s ability to reproduce or its ability to transmit parasites. The class of gene drive in question uses a very precise cutting and copying mechanism, termed ‘homing’, that allows it to increase its numbers in the cells that go on to form sperm or eggs, thereby increasing the chances that a copy of the gene drive is transmitted to offspring. However, while this type of gene drive can rapidly invade a mosquito population, mosquitoes can also eventually become resistant to the gene drive in some cases. Here we show that restricting the cutting activity of the gene drive to the germline tissue is crucial to maintaining its potency and we illustrate how failure to restrict this activity can lead to the generation of mutations that can make mosquitoes resistant to the gene drive.

© The Author(s) 2020., The causes and consequences of genome reduction in animals are unclear because our understanding of this process mostly relies on lineages with often exceptionally high rates of evolution. Here, we decode the compact 73.8-megabase genome of Dimorphilus gyrociliatus, a meiobenthic segmented worm. The D. gyrociliatus genome retains traits classically associated with larger and slower-evolving genomes, such as an ordered, intact Hox cluster, a generally conserved developmental toolkit and traces of ancestral bilaterian linkage. Unlike some other animals with small genomes, the analysis of the D. gyrociliatus epigenome revealed canonical features of genome regulation, excluding the presence of operons and trans-splicing. Instead, the gene-dense D. gyrociliatus genome presents a divergent Myc pathway, a key physiological regulator of growth, proliferation and genome stability in animals. Altogether, our results uncover a conservative route to genome compaction in annelids, reminiscent of that observed in the vertebrate Takifugu rubripes., This study was supported by Sars Centre core budget and the European Research Council (ERC) grant agreement no. 648861 to A.H. J.M.M.-D. was additionally supported by the ERC grant agreement no. 801669, and B.C.V. by an EMBO Long-Term Fellowship (ALTF 74-2018). J.L.G.-S. received funding from the ERC (grant agreement no. 740041) and the Spanish Ministerio de Economía y Competitividad (grant no. BFU2016-74961-P) and the institutional grant Unidad de Excelencia María de Maeztu (MDM-2016-0687). We dedicate this manuscript to the loving memory of José Luis Gómez-Skarmeta.

MalariaGEN is a data-sharing network that enables groups around the world to work together on the genomic epidemiology of malaria. Here we describe a new release of curated genome variation data on 7,000 Plasmodium falciparum samples from MalariaGEN partner studies in 28 malaria-endemic countries. High-quality genotype calls on 3 million single nucleotide polymorphisms (SNPs) and short indels were produced using a standardised analysis pipeline. Copy number variants associated with drug resistance and structural variants that cause failure of rapid diagnostic tests were also analysed. Almost all samples showed genetic evidence of resistance to at least one antimalarial drug, and some samples from Southeast Asia carried markers of resistance to six commonly-used drugs. Genes expressed during the mosquito stage of the parasite life-cycle are prominent among loci that show strong geographic differentiation. By continuing to enlarge this open data resource we aim to facilitate research into the evolutionary processes affecting malaria control and to accelerate development of the surveillance toolkit required for malaria elimination.

Escherichia coli is a highly diverse organism that includes a range of commensal and pathogenic variants found across a range of niches and worldwide. In addition to causing severe intestinal and extraintestinal disease, E. coli is considered a priority pathogen due to high levels of observed drug resistance. The diversity in the E. coli population is driven by high genome plasticity and a very large gene pool. All these have made E. coli one of the most well-studied organisms, as well as a commonly used laboratory strain. Today, there are thousands of sequenced E. coli genomes stored in public databases. While data is widely available, accessing the information in order to perform analyses can still be a challenge. Collecting relevant available data requires accessing different sources, where data may be stored in a range of formats, and often requires further manipulation and processing to apply various analyses and extract useful information. In this study, we collated and intensely curated a collection of over 10 000 E. coli and Shigella genomes to provide a single, uniform, high-quality dataset. Shigella were included as they are considered specialized pathovars of E. coli . We provide these data in a number of easily accessible formats that can be used as the foundation for future studies addressing the biological differences between E. coli lineages and the distribution and flow of genes in the E. coli population at a high resolution. The analysis we present emphasizes our lack of understanding of the true diversity of the E. coli species, and the biased nature of our current understanding of the genetic diversity of such a key pathogen.

Insulin resistance (IR) precedes the development of type 2 diabetes (T2D) and increases cardiovascular disease risk. Although genome wide association studies (GWAS) have uncovered new loci associated with T2D, their contribution to explain the mechanisms leading to decreased insulin sensitivity has been very limited. Thus, new approaches are necessary to explore the genetic architecture of insulin resistance. To that end, we generated an iPSC library across the spectrum of insulin sensitivity in humans. RNA-seq based analysis of 310 induced pluripotent stem cell (iPSC) clones derived from 100 individuals allowed us to identify differentially expressed genes between insulin resistant and sensitive iPSC lines. Analysis of the co-expression architecture uncovered several insulin sensitivity-relevant gene sub-networks, and predictive network modeling identified a set of key driver genes that regulate these co-expression modules. Functional validation in human adipocytes and skeletal muscle cells (SKMCs) confirmed the relevance of the key driver candidate genes for insulin responsiveness., Author summary Insulin resistance is characterized by a defective response (“resistance”) to normal insulin concentrations to uptake the glucose present in the blood, and is the underlying condition that leads to type 2 diabetes (T2D) and increases the risk of cardiovascular disease. It is estimated that 25–33% of the US population are insulin resistant enough to be at risk of serious clinical consequences. For more than a decade, large population studies have tried to discover the genes that participate in the development of insulin resistance, but without much success. It is now increasingly clear that the complex genetic nature of insulin resistance requires novel approaches centered in patient specific cellular models. To fill this gap, we have generated an induced pluripotent stem cell (iPSC) library from individuals with accurate measurements of insulin sensitivity, and performed gene expression and key driver analyses. Our work demonstrates that iPSCs can be used as a revolutionary technology to model insulin resistance and to discover key genetic drivers. Moreover, they can develop our basic knowledge of the disease, and are ultimately expected to increase the therapeutic targets to treat insulin resistance and type 2 diabetes.

The known species richness of oak gallwasps in Asia has increased tremendously in the past decade. However, the vast majority of taxa have been described from the east coast of Asia, and knowledge of oak gallwasps from Central Asia is still scant. Here we use molecular and morphological characters to describe a new genus of cynipid oak gallwasp, Heocynips Fang, Nieves-Aldrey, and Melika, new genus (Hymenoptera: Cynipidae: Cynipini), and a new species Heocynips furvoaurantius Tang, Nicholls, and Stone, new species from Emeishan, Sichuan Province, China. We present diagnostic characters for the new genus, including adult morphology and DNA sequence data, and provide biological information on host plant association, gall morphology and developmental phenology for the new species., Workin China by GNS, GM, FS and CTT was supported by the Fundamental Research Funds for the Central Universities of P.R.China, and NERC grants NE/E014453/1and NE/T000120/1 to GNS.

International audience; CRISPR guide RNAs (gRNAs) can be programmed with relative ease to allow the genetic editing of nearly any DNA or RNA sequence. Here, we propose novel molecular architectures to achieve RNA-dependent modulation of CRISPR activity in response to specific RNA molecules. We designed and tested, in both living Escherichia coli cells and cell-free assays for rapid prototyping, cis-repressed RNA-interacting guide RNA (igRNA) that switch to their active state only upon interaction with small RNA fragments or long RNA transcripts, including pathogen-derived mRNAs of medical relevance such as the human immunodeficiency virus infectivity factor. The proposed CRISPR-igRNAs are fully customizable and easily adaptable to the majority if not all the available CRISPR-Cas variants to modulate a variety of genetic functions in response to specific cellular conditions, providing orthogonal activation and increased specificity. We thereby foresee a large scope of application for therapeutic, diagnostic, and biotech applications in both prokaryotic and eukaryotic systems.

Objectives\ud Mobile colistin resistance genes (mcr) encoded on conjugative plasmids, although described only relatively recently, have been reported globally in both humans and livestock. The genes are often associated with ISApl1 which can transpose the genes to novel genetic locations. Since the first report, multiple variants of mcr have been discovered on a variety of genetic locations in Escherichia coli, in plasmids and integrated into the chromosome. \ud Methods\ud Using hybrid assembly of short-read and long-read WGS data, the presence of mcr-1 was confirmed on an IncI1 plasmid in E. coli. In vitro conjugation assays were performed to determine the potential to transfer between strains. Genetic comparison to previously reported IncI1 plasmids was performed.\ud Results\ud The genomic sequence has identified that mcr-1 is present on a complete IncI1 plasmid. Comparison to previously reported ESBL-encoding plasmids from E. coli in the Netherlands from the same time period indicated a distinct lineage for this plasmid.\ud Conclusion\ud The observation of mcr-1¬ on an IncI1 plasmid confirms that the genetic region of this gene is actively transposed between genetic locations. This active transposition has consequences for the study of the epidemiology of mcr in populations.

Carbapenemases inactivate most β-lactam antibiotics, including carbapenems, and have frequently been reported among Enterobacteriaceae , Acinetobacter spp. and Pseudomonas spp. Traditionally, the horizontal gene transfer of carbapenemase-encoding genes (CEGs) has been linked to plasmids. However, given that integrative and conjugative elements (ICEs) are possibly the most abundant conjugative elements among prokaryotes, we conducted an in silico analysis to ascertain the likely role of ICEs in the spread of CEGs among all bacterial genomes (n=182 663). We detected 17 520 CEGs, of which 66 were located within putative ICEs among several bacterial species (including clinically relevant bacteria, such as Pseudomonas aeruginosa , Klebsiella pneumoniae and Escherichia coli ). Most CEGs detected within ICEs belong to the IMP, NDM and SPM metallo-beta-lactamase families, and the serine beta-lactamase KPC and GES families. Different mechanisms were likely responsible for acquisition of these genes. The majority of CEG-bearing ICEs belong to the MPFG, MPFT and MPFF classes and often encode resistance to other antibiotics (e.g. aminoglycosides and fluoroquinolones). This study provides a snapshot of the different CEGs associated with ICEs among available bacterial genomes and sheds light on the underappreciated contribution of ICEs to the spread of carbapenem resistance globally.

Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) are diverse immune systems found in many prokaryotic genomes that target invading foreign DNA such as bacteriophages and plasmids. There are multiple types of CRISPR with arguably the most enigmatic being Type IV. During an investigation of CRISPR carriage in clinical, multi-drug resistant, Klebsiella pneumoniae, a Type IV-A3 CRISPR-Cas system was detected on plasmids from two K. pneumoniae isolates from Egypt (isolated in 2002-2003) and a single K. pneumoniae isolate from the United Kingdom (isolated in 2017). Sequence analysis of all other genomes available in GenBank revealed that this CRISPR-Cas system was present on 28 other plasmids from various Enterobacteriaceae hosts and was never found on a bacterial chromosome. This system is exclusively located on IncHI1B/IncFIB plasmids and is associated with multiple putative transposable elements. Expression of the cas loci was confirmed in the available clinical isolates by RT-PCR. In all cases, the CRISPR-Cas system has a single CRISPR array (CRISPR1) upstream of the cas loci which has several, conserved, spacers which, amongst things, match regions within conjugal transfer genes of IncFIIK/IncFIB(K) plasmids. Our results reveal a Type IV-A3 CRISPR-Cas system exclusively located on IncHI1B/IncFIB plasmids in Enterobacteriaceae that is likely to be able to target IncFIIK/IncFIB(K) plasmids presumably facilitating intracellular, inter-plasmid competition.

RNA editing is a source of molecular diversity that regulates the functional repertoire of animal transcriptomes. Multiple studies in Drosophila have revealed that conserved editing events can be a 10 source of evolutionary ad aptations, and there is a solid body of evidence linking editing and the fine tuning of neural genes, which are often targeted by insecticides used in vector control. Yet, 12 despite these suggestive connections, genome wide analyses of editing in insect vect ors are conspicuously lacking. Future advances will require complementing the growing wealth of vector 14 genomes with targeted transcriptome analyses. Here, we review recent investigations of the genetic\ud footprints of adaptive RNA editing in insects and prov ide an overview of new methodologies 16 applicable to studies of RNA editing in insect vectors.

Mobile genetic elements (MGEs) are key factors responsible for dissemination of virulence determinants and antimicrobial-resistance genes amongst pathogenic bacteria. Conjugative MGEs are notable for their high gene loads donated per transfer event, broad host ranges and phylogenetic ubiquity amongst prokaryotes, with the subclass of chromosomally inserted integrative and conjugative elements (ICEs) being particularly abundant. The focus on a small number of model systems has biased the study of ICEs towards those conferring readily selectable phenotypes to host cells, whereas the identification and characterization of integrated cryptic elements remains challenging. Even though antimicrobial resistance and horizontally acquired virulence genes are major factors aggravating neisserial infection, conjugative MGEs of Neisseria gonorrhoeae and Neisseria meningitidis remain poorly characterized. Using a phenotype-independent approach based on atypical distributions of DNA uptake sequences (DUSs) in MGEs relative to the chromosomal background, we have identified two groups of chromosomally integrated conjugative elements in Neisseria: one found almost exclusively in pathogenic species possibly deriving from the genus Kingella, the other belonging to a group of Neisseria mucosa-like commensals. The former element appears to enable transfer of traditionally gonococcal-specific loci such as the virulence-associated toxin–antitoxin system fitAB to N. meningitidis chromosomes, whilst the circular form of the latter possesses a unique attachment site (attP) sequence seemingly adapted to exploit DUS motifs as chromosomal integration sites. In addition to validating the use of DUS distributions in Neisseriaceae MGE identification, the >170 identified ICE sequences provide a valuable resource for future studies of ICE evolution and host adaptation.

Mutations are incorporated into the genomes of RNA viruses at an optimal frequency and altering this precise frequency has been proposed as a strategy to create live-attenuated vaccines. However, determining the effect of specific mutations that alter fidelity has been difficult because of the rapid selection of the virus population during replication. By deleting residues of the structural polyprotein PE2 cleavage site, E3Δ56-59, in Venezuelan equine encephalitis virus (VEEV) TC-83 vaccine strain, non-infectious virus particles were used to assess the effect of single mutations on mutation frequency without the interference of selection that results from multiple replication cycles. Next-generation sequencing analysis revealed a significantly lower frequency of transversion mutations and overall mutation frequency for the fidelity mutants compared to VEEV TC-83 E3Δ56-59. We demonstrate that deletion of the PE2 cleavage site halts virus infection while making the virus particles available for downstream sequencing. The conservation of the site will allow the evaluation of suspected fidelity mutants across alphaviruses of medical importance.

Tn916 is a conjugative transposon (CTn) and the first reported and most well characterised of the Tn916/Tn1545 family of CTns. Tn916-like elements have a characteristic modular structure and different members of this family have been identified based on similarities and variations in these modules. In addition to carrying genes encoding proteins required for their conjugation, Tn916-like elements also carry accessory, antimicrobial resistance genes, most commonly the tetracycline resistance gene, tet(M). Our study aimed to identify and characterise tetracycline resistance genes from the human saliva metagenome using a functional metagenomic approach. We identified a tetracycline-resistant clone, TT31, the sequencing of which revealed it to encode both tet(M) and tet(L). Comparison of the TT31 sequence with the accessory, regulation, and recombination modules of other Tn916-like elements indicated that a partial Tn916-like element encoding a truncated orf9 was cloned in TT31. Analysis indicated that a previous insertion within the truncated orf9 created the full length orf9 found in most Tn916-like transposons, demonstrating that orf9 is, in fact, the result of a gene fusion event. Thus, we hypothesise that the Tn916-like element cloned in TT31 likely represents an ancestral Tn916.

Increasing insecticide resistance in malaria-transmitting vectors represents a public health threat, but underlying mechanisms are poorly understood. Here, a data integration approach is used to analyse transcriptomic data from comparisons of insecticide resistant and susceptible Anopheles populations from disparate geographical regions across the African continent. An unbiased, integrated analysis of this data confirms previously described resistance candidates but also identifies multiple novel genes involving alternative resistance mechanisms, including sequestration, and transcription factors regulating multiple downstream effector genes, which are validated by gene silencing. The integrated datasets can be interrogated with a bespoke Shiny R script, deployed as an interactive web-based application, that maps the expression of resistance candidates and identifies co-regulated transcripts that may give clues to the function of novel resistance-associated genes., Increasing insecticide resistance of mosquitoes represents a public health threat, and underlying mechanisms are poorly understood. Here, Ingham et al. identify putative insecticide resistance genes in Anopheles gambiae populations across Africa and develop a web-based application that maps their expression.

Resistance in to members of all 4 major classes (pyrethroids, carbamates, organochlorines, and organophosphates) of public health insecticides limits effective control of malaria transmission in Africa. Increase in expression of detoxifying enzymes has been associated with insecticide resistance, but their direct functional validation in is still lacking. Here, we perform transgenic analysis using the GAL4/UAS system to examine insecticide resistance phenotypes conferred by increased expression of the 3 genes-, , and -most often found up-regulated in resistant We report evidence in that organophosphate and organochlorine resistance is conferred by overexpression of GSTE2 in a broad tissue profile. Pyrethroid and carbamate resistance is bestowed by similar overexpression, and confers only pyrethroid resistance when overexpressed in the same tissues. Conversely, such overexpression increases susceptibility to the organophosphate malathion, presumably due to conversion to the more toxic metabolite, malaoxon. No resistant phenotypes are conferred when either gene overexpression is restricted to the midgut or oenocytes, indicating that neither tissue is involved in insecticide resistance mediated by the candidate P450s examined. Validation of genes conferring resistance provides markers to guide control strategies, and the observed negative cross-resistance due to gives credence to proposed dual-insecticide strategies to overcome pyrethroid resistance. These transgenic -resistant lines are being used to test the "resistance-breaking" efficacy of active compounds early in their development.

In Argentina, NDM metallo-β-lactamase was first reported in 2013. By now, it has disseminated throughout the country in diverse Gram negative bacteria. Here, we report the case of a paediatric patient that underwent a 1-year hospitalisation due to erythrodermic psoriasis in 2014 and received multiple antimicrobial treatments. During his stay, five isolates were obtained from rectal swabs (rs) or blood culture (bc) suspicious of carbapenemase production: a K. quasipneumoniae subsp. quasipneumoniae (rs), Citrobacter freundii (rs), Escherichia coli (bc), Enterobacter cloacae (rs), and a Serratia marcescens (bc). The isolates were studied with broth microdilution, biparental conjugation and plasmid and whole genome sequencing (Illumina). All isolates harboured an 138,998-bp type 1 IncC plasmid that carried blaNDM-1, bleMBL, blaCMY-6, rmtC, aac(6')-Ib, and sul1 resistance genes. Additionally, the blaNDM-plasmids contained ISKpn8 an insertion sequence previously described as associated only to blaKPC. One isolate, a colistin-resistant E. coli, also carried a mcr-1-containing an IncI2 plasmid, which did not harbour additional resistance. The whole genome of K. quasipneumoniae subsp. quasipneumoniae isolate was fully sequenced. This isolate harboured, additionally to blaNDM, three plasmid-mediated quinolone resistance genes: qnrB4, qnrB52 and aac(6')-Ib-cr1. The E. cloacae isolate also harboured qnrA1. These findings alert to the underestimated horizontal dissemination of multidrug-resistant plasmids limiting treatment options with last resort antimicrobials.

Lan-Lan Zhong,1,2,* Qian Zhou,3,*Cui-Yan Tan,3,*Adam P Roberts,4,5Mohamed Abd El-Gawad El-Sayed Ahmed,1–2,6 Guanping Chen,1,2Min Dai,7 Fan Yang,8 Yong Xia,9 Kang Liao,10 Yingjian Liang,3 Yongqiang Yang,1,11Siyuan Feng,1,2Xiaobin Zheng,3,*Guo-Bao Tian1,2,*1Program in Pathobiology, The Fifth Affiliated Hospital, Zhongshan School of Medicine, Sun Yat-Sen University, Guangdong 519000, People’s Republic of China; 2Ministry of Education, Key Laboratory of Tropical Diseases Control, Sun Yat-sen University, Guangzhou 510080, People’s Republic of China; 3Department of Respiratory Medicine, The Fifth Affiliated Hospital of Sun Yat-sen University, Zhuhai 519000, People’s Republic of China; 4Department of Tropical Disease Biology, Liverpool School of Tropical Medicine, Liverpool, L3 5QA, UK; 5Centre for Drugs and Diagnostics, Liverpool School of Tropical Medicine, Liverpool L3 5QA, UK; 6Department of Microbiology and Immunology, Faculty of Pharmaceutical Sciences and Drug Manufacturing, Misr University for Science and Technology (MUST), Cairo, Egypt; 7School of Laboratory Medicine, Chengdu Medical College, Chengdu 610500, People’s Republic of China; 8Department of Microbiology, School of Basic Medical Science, Xinxiang Medical University, Xinxiang 453003, People’s Republic of China; 9Department of Clinical Laboratory Medicine, The Third Affiliated Hospital of Guangzhou Medical University, Guangzhou 510150, People’s Republic of China; 10Department of Clinical Laboratory, The First Affiliated Hospital of Sun Yat-Sen University, Guangzhou 510080, People’s Republic of China; 11School of Pharmaceutical Sciences (Shenzhen), Sun Yat-sen University, Guangzhou 510006, People’s Republic of China*These authors contributed equally to this workPurpose: The discovery of the plasmid-mediated colistin resistance genes, mcr, revealed a mechanism of transmission of colistin resistance, which is a major, global public health concern especially among individuals infected with carbapenem-resistant Gram-negative bacteria. To monitor the spread and epidemiology of mcr genes, a convenient and reliable method to detect mcr genes in clinical isolates is needed, especially in the primary care institutions. This study aimed to establish a restriction endonuclease-based multiplex loop-mediated isothermal amplification (multi-LAMP) assay to detect mcr genes (mcr-1 to mcr-5) harbored by colistin-resistant bacteria.Methods: A triple-LAMP assay for mcr-1, mcr-3, and mcr-4 and a double-LAMP assay for mcr-2 and mcr-5 were established. The sensitivity and specificity of the LAMP reactions were determined via electrophoresis and visual detection.Results: The sensitivity of the LAMP assay was 10-fold greater than that of PCR, with high specificity among the screened primers. Specific mcr genes were distinguished in accordance with band numbers and the fragment length of the digested LAMP amplification products. Furthermore, the LAMP assay was confirmed as a rapid and reliable diagnostic technique upon application for clinical samples, and the results were consistent with those of conventional PCR assay.Conclusion: The multi-LAMP assay is a potentially promising method to detect mcr genes and will, if implemented, help prevent infections by drug-resistant bacteria in primary-care hospitals due to rapid and reliable surveillance. To our knowledge, this is the first study to report the application of LAMP to detect mcr-2 to mcr-5 genes and the first time that multi-LAMP has been applied to detect mcr genes.Keywords: mcr genes, colistin resistance, multi-LAMP, rapid detection, enzyme digestion

Integrons are genetic elements consisting of a functional platform for recombination and expression of gene cassettes (GCs). GCs usually carry promoter-less open reading frames (ORFs), encoding proteins with various functions including antibiotic resistance. The transcription of GCs relies mainly on a cassette promoter (PC), located upstream of an array of GCs. Some integron GCs, called ORF-less GCs, contain no identifiable ORF with a small number shown to be involved in antisense mRNA mediated gene regulation.In this study, promoter sequences were identified, usingin silicoanalysis, within GCs PCR amplified from the oral metagenome. The promoter activity of ORF-less GCs was verified by cloning them upstream of agusAreporter, proving they can function as a promoter, presumably allowing bacteria to adapt to multiple stresses within the complex physico-chemical environment of the human oral cavity. A bi-directional promoter detection system was also developed allowing direct identification of clones with promoter-containing GCs on agar plates. Novel promoter-containing GCs were identified from the human oral metagenomic DNA using this construct, called pBiDiPD.This is the first demonstration and detection of promoter activity of ORF-less GCs and the development of an agar plate-based detection system will enable similar studies in other environments.

Background Salmonella enterica serovar Enteritidis is a cause of both poultry- and egg-associated enterocolitis globally and bloodstream-invasive nontyphoidal Salmonella (iNTS) disease in sub-Saharan Africa (sSA). Distinct, multi-drug resistant genotypes associated with iNTS disease in sSA have recently been described, often requiring treatment with fluoroquinolone antibiotics. In industrialised countries, antimicrobial use in poultry production has led to frequent fluoroquinolone resistance amongst globally prevalent enterocolitis-associated lineages. Methodology/Principal findings Twenty seven S. Enteritidis isolates from patients with iNTS disease and two poultry isolates, collected between 2007 and 2015 in the Ashanti region of Ghana, were whole-genome sequenced. These isolates, notable for a high rate of diminished ciprofloxacin susceptibility (DCS), were placed in the phyletic context of 1,067 sequences from the Public Health England (PHE) S. Enteritidis genome database to understand whether DCS was associated with African or globally-circulating clades of S. Enteritidis. Analysis showed four of the major S. Enteritidis clades were represented, two global and two African. All thirteen DCS isolates, containing a single gyrA mutation at codon 87, belonged to a global PT4-like clade responsible for epidemics of poultry-associated enterocolitis. Apart from two DCS isolates, which clustered with PHE isolates associated with travel to Spain and Brazil, the remaining DCS isolates, including one poultry isolate, belonged to two monophyletic clusters in which gyrA 87 mutations appear to have developed within the region. Conclusions/Significance Extensive phylogenetic diversity is evident amongst iNTS disease-associated S. Enteritidis in Ghana. Antimicrobial resistance profiles differed by clade, highlighting the challenges of devising empirical sepsis guidelines. The detection of fluoroquinolone resistance in phyletically-related poultry and human isolates is of major concern and surveillance and control measures within the region’s burgeoning poultry industry are required to protect a human population at high risk of iNTS disease., Author summary Salmonella enterica serovar Enteritidis is both a prominent global cause of zoonotic gastroenteritis, in association with the industrial production of eggs and poultry, and of bloodstream-invasive infection in sub-Saharan Africa, a clinical syndrome referred to as invasive nontyphoidal Salmonella (iNTS) disease. African epidemic iNTS strains are frequently multi-drug resistant, leaving fluoroquinolone antibiotics as key agents in reducing iNTS-associated morbidity and mortality. We analysed the genomes of S. Enteritidis collected in Ghana from patients with iNTS disease to investigate the emergence of fluoroquinolone resistance in the region. Extensive phylogenetic diversity was present, however, fluoroquinolone resistance was confined to a single clade causing global epidemics of poultry-associated gastroenteritis. This resistance has predominantly emerged locally, rather than being imported. We found that antimicrobial resistance patterns differed by S. Enteritidis clade, highlighting the challenges of devising empirical sepsis treatment guidelines in the absence of diagnostic microbiology facilities to monitor changes in resistance profiles. Furthermore, we detected fluoroquinolone resistance in closely related poultry and human isolates, suggesting a role of antimicrobial use within the growing local poultry industry in driving the emergence of resistance and a need for surveillance measures within this industry to protect a human population at high risk of iNTS disease.

The goliath frog (Conraua goliath) is an Endangered species exclusively found in Cameroon and Equatorial Guinea. Climate change, deforestation and overhunting are principal causes driving this species to extinction. Therefore, a better understanding of the genetic diversity and population structure of this species is necessary to improve conservation efforts. Here we used two mitochondrial genes (Cytochrome Oxidase subunit 1 (COI) and 16S) extracted from 54 C. goliath individuals from six localities in Cameroon to examine their genetic diversity. The result shows a low DNA substitution between the sequences. There were four 16S and two COI haplotypes in total. Overall, genetic diversity was very low for all the genes with nucleotide diversity of 0.00106 and 0.00007 for 16S and COI respectively. The Tajima D and Fu Fs statistics were negative. The TCS haplotype network revealed a predominant and ancestral haplotype (H1) for these genes which is distributed in the 6 populations. Pairwise genetic differentiation (F) generated between these populations using 16S revealed very high differentiation between populations from Nkam and Mungo Administrative Divisions in Cameroon. In contrast, we observed low differentiation among the geographically clustered Mungo and Nkam populations. Overall, human activities and perhaps climate change can appear to have depleted genetic diversity in the scattered populations that remain of this amphibian. To sustain the Goliath frog, we suggest to the Cameroonian government to implement more effective strategies to conserve and manage remnant populations of this iconic species through more effort against poaching which contribute to reduce the genetic diversity.

Vector-borne disease control relies on efficient vector surveillance, mostly carried out using traps whose number and locations are often determined by expert opinion rather than a rigorous quantitative sampling design. In this work we propose a framework for ecological sampling design which in its preliminary stages can take into account environmental conditions obtained from open data (i.e. remote sensing and meteorological stations) not necessarily designed for ecological analysis. These environmental data are used to delimit the area into ecologically homogeneous strata. By employing Bayesian statistics within a model-based sampling design, the traps are deployed among the strata using a mixture of random and grid\ud locations which allows balancing predictions and model-fitting accuracies. Sample sizes and the effect of ecological strata on sample sizes are estimated from previous mosquito sampling campaigns open data. Notably, we found that a configuration of 30 locations with four households each (120 samples) will have a similar accuracy in the predictions of mosquito abundance as 200 random samples. In addition, we show that random sampling independently from ecological strata, produces biased estimates of the mosquito\ud abundance. Finally, we propose standardizing reporting of sampling\ud designs to allow transparency and repetition/re-use in subsequent sampling campaigns

BACKGROUND\ud Cryptococcal meningitis (CM) causes an estimated 180,000 deaths annually, predominantly in sub-Saharan Africa, where most patients receive fluconazole (FLC) monotherapy. While relapse after FLC monotherapy with resistant strains is frequently observed, the mechanisms and impact of emergence of FLC resistance in human CM are poorly understood. Heteroresistance (HetR) - a resistant subpopulation within a susceptible strain - is a recently described phenomenon in Cryptococcus neoformans (Cn) and Cryptococcus gattii (Cg), the significance of which has not previously been studied in humans.\ud \ud METHODS\ud A cohort of 20 patients with HIV-associated CM in Tanzania was prospectively observed during therapy with either FLC monotherapy or in combination with flucytosine (5FC). Total and resistant subpopulations of Cryptococcus spp. were quantified directly from patient cerebrospinal fluid (CSF). Stored isolates underwent whole genome sequencing and phenotypic characterization.\ud \ud RESULTS\ud Heteroresistance was detectable in Cryptococcus spp. in the CSF of all patients at baseline (i.e., prior to initiation of therapy). During FLC monotherapy, the proportion of resistant colonies in the CSF increased during the first 2 weeks of treatment. In contrast, no resistant subpopulation was detectable in CSF by day 14 in those receiving a combination of FLC and 5FC. Genomic analysis revealed high rates of aneuploidy in heteroresistant colonies as well as in relapse isolates, with chromosome 1 (Chr1) disomy predominating. This is apparently due to the presence on Chr1 of ERG11, which is the FLC drug target, and AFR1, which encodes a drug efflux pump. In vitro efflux levels positively correlated with the level of heteroresistance.\ud \ud CONCLUSION\ud Our findings demonstrate for what we believe is the first time the presence and emergence of aneuploidy-driven FLC heteroresistance in human CM, association of efflux levels with heteroresistance, and the successful suppression of heteroresistance with 5FC/FLC combination therapy.\ud \ud FUNDING\ud This work was supported by the Wellcome Trust Strategic Award for Medical Mycology and Fungal Immunology 097377/Z/11/Z and the Daniel Turnberg Travel Fellowship.

The surveillance of antimicrobial resistant isolates has proven to be one of the most valuable tools to understand the global rise of multidrug-resistant bacterial pathogens. We report the first insights into the current situation in the Caribbean, where a pilot project to monitor antimicrobial resistance through phenotypic resistance measurements combined with whole-genome sequencing was set up in collaboration with the Caribbean Public Health Agency (CARPHA). Our first study focused on Klebsiella pneumoniae, a highly relevant organism amongst the Gram-negative opportunistic pathogens world-wide today causing hospital, as well as community-acquired, infections. Our results show that not only carbapenem resistance, but also hypervirulent strains, are circulating in patients in the Caribbean. Our current data does not allow us to infer their prevalence in the population. We argue for the urgent need to further support antimicrobial resistance surveillance and stewardship in this almost uncharted territory, which can make a significant impact on the reduction of antimicrobial usage.DATA SUMMARY-Raw sequencing data is deposited at the sequence read archive (SRA), and assemblies are deposited at GenBank, accession numbers for all are given in Dataset S1.-The data of measured resistance phenotypes (Vitek) is also provided in Dataset S1.-The tree file and associated metadata can be investigated and downloaded through the free online platform microreact (https://microreact.org/project/S1-a7KAkV).✓ I/We confirm all supporting data, code and protocols have been provided within the article or through supplementary data files. ⊠Significance as a BioResource to the communityThis BioResource contains whole-genome sequence data of 270 Klebsiella pneumoniae isolates, information about encoded resistance genes and the phylogeny of the isolates, their distribution in the global K. pneumoniae population, and their resistance phenotype data as determined by the VITEK 2 compact system. The isolates are recent (2017 through 2018) and represent clinically relevant patient isolates from 15 different sites in 12 Caribbean states. These data will be of interest for researchers working on K. pneumoniae and other opportunistic pathogens as well as those interested in mobile genetic elements carrying antimicrobial resistance cassettes. Our data is the only recent survey of antimicrobial resistant opportunistic pathogens from multiple sites within the Caribbean, and is of high significance for the global surveillance of K. pneumoniae and antimicrobial resistance elements. This BioResource is made available through data tables provided with this article, as well deposition of the raw data in the relevant archives, and an interactive platform (microreact) to enquire and download analyses (phylogenetic tree, metadata).

Decades of unmanaged insecticide use and routine exposure to agrochemicals have left many populations of malaria vectors in the Americas resistant to multiple classes of insecticides, including pyrethroids. The molecular basis of pyrethroid resistance is relatively uncharacterised in American malaria vectors, preventing the design of suitable resistance management strategies. Using whole transcriptome sequencing, we characterized the mechanisms of pyrethroid resistance in Anopheles albimanus from Peru and Guatemala. An. albimanus were phenotyped as either deltamethrin or alpha-cypermethrin resistant. RNA from 1) resistant, 2) unexposed, and 3) a susceptible laboratory strain of An. albimanus was sequenced and analyzed using RNA-Seq. Expression profiles of the three groups were compared based on the current annotation of the An. albimanus reference genome. Several candidate genes associated with pyrethroid resistance in other malaria vectors were found to be overexpressed in resistant An. albimanus. In addition, gene ontology terms related to serine-type endopeptidase activity, extracellular activity and chitin metabolic process were also commonly overexpressed in the field caught resistant and unexposed samples from both Peru and Guatemala when compared to the susceptible strain. The cytochrome P450 CYP9K1 was overexpressed 14x in deltamethrin and 8x in alpha-cypermethrin-resistant samples from Peru and 2x in deltamethrin-resistant samples from Guatemala, relative to the susceptible laboratory strain. CYP6P5 was overexpressed 68x in deltamethrin-resistant samples from Peru but not in deltamethrin-resistant samples from Guatemala. When comparing overexpressed genes between deltamethrin-resistant and alpha-cypermethrin-resistant samples from Peru, a single P450 gene, CYP4C26, was overexpressed 9.8x (p

One of the factors that can affect conjugation of IncI1 plasmids, amongst others, is the genetic region known as the shufflon. This multiple inversion system modifies the pilus tip proteins used during conjugation, thus affecting the affinity for different recipient cells. Although recombination is known to occur in in vitro conditions, little is known about the regulation and the extent of recombination that occurs. To measure the recombination of the shufflon, we have amplified the entire shufflon region and sequenced the amplicons using nanopore long-read sequencing. This method was effective to determine the order of the segments of the shufflon and allow for the analysis of the shufflon variants that are present in a heterogeneous pool of templates. Analysis was performed over different growth phases and after addition of cefotaxime. Furthermore, analysis was performed in different E. coli host cells to determine if recombination is likely to be influenced. Recombination of the shufflon was constantly ongoing in all conditions that were measured, although no differences in the amount of different shufflon variants or the rate at which novel variants were formed could be found. As previously reported, some variants were abundant in the population while others were scarce. This leads to the conclusion that the shufflon is continuously recombining at a constant rate, or that the method used here was not sensitive enough to detect differences in this rate. For one of the plasmids, the host cell appeared to have an effect on the specific shufflon variants that were formed which were not predominant in another host, indicating that host factors may be involved. As previously reported, the pilV-A and pilV-A' ORFs are formed at higher frequencies than other pilV ORFs. These results demonstrate that the recombination that occurs within the shufflon is not random. While any regulation of the shufflon affected by these in vitro conditions could not be revealed, the method of amplifying large regions for long-read sequencing for the analysis of multiple inversion systems proved effective.

Background: Tn916 and Tn1545 are paradigms of a large family of related, broad host range, conjugative transposons that are widely distributed in bacteria and contribute to the spread of antibiotic resistance genes (ARGs). Variation in the copy number (CN) of Tn916-Tn1545 elements and the circularization ratio (CR) may play an important role in propagation of ARGs carried by these elements. Objectives and Design: In this study, the CN and CR of Tn916-Tn1545 elements in oral streptococci were determined using droplet digital PCR (ddPCR). In addition, we investigated the influence of tetracycline on the CR of Tn916-Tn1545 elements. Results: The ddPCR assay designed in this study is a reliable way to rapidly determine CN and CR of Tn916-Tn1545 elements. Conclusions: Our data also suggest that Tn916-Tn1545 elements are generally stable without selective pressure in the clinical oral Streptococcus strains investigated in this study.

Colubridae represents the most phenotypically diverse and speciose family of snakes, yet no well-assembled and annotated genome exists for this lineage. Here we report and analyze the genome of the garter snake, Thamnophis sirtalis, a colubrid snake that is an important model species for research in evolutionary biology, physiology, genomics, behavior and the evolution of toxin resistance. Using the garter snake genome, we show how snakes have evolved numerous adaptations for sensing and securing prey, and identify features of snake genome structure that provide insight into the evolution of amniote genomes. Analyses of the garter snake and other squamate reptile genomes highlight shifts in repeat element abundance and expansion within snakes, uncover evidence of genes under positive selection, and provide revised neutral substitution rate estimates for squamates. Our identification of Z and W sex chromosome-specific scaffolds provides evidence for multiple origins of sex chromosome systems in snakes and demonstrates the value of this genome for studying sex chromosome evolution. Analysis of gene duplication and loss in visual and olfactory gene families supports a dim-light ancestral condition in snakes and indicates that olfactory receptor repertoires underwent an expansion early in snake evolution. Additionally, we provide some of the first links between secreted venom proteins, the genes that encode them, and their evolutionary origins in a rear-fanged colubrid snake, together with new genomic insight into the coevolutionary arms race between garter snakes and highly toxic newt prey that led to toxin resistance in garter snakes.

Streptococcus pneumoniae, a normal commensal of the upper respiratory tract, is a major public health concern, responsible for substantial global morbidity and mortality due to pneumonia, meningitis and sepsis. Why some pneumococci invade the bloodstream or CSF (so-called invasive pneumococcal disease; IPD) is uncertain. In this study we identify genes associated with IPD. We transform whole genome sequence (WGS) data into a sequence typing scheme, while avoiding the caveat of using an arbitrary genome as a reference by substituting it with a constructed pangenome. We then employ a random forest machine-learning algorithm on the transformed data, and find 43 genes consistently associated with IPD across three geographically distinct WGS data sets of pneumococcal carriage isolates. Of the genes we identified as associated with IPD, we find 23 genes previously shown to be directly relevant to IPD, as well as 18 uncharacterized genes. We suggest that these uncharacterized genes identified by us are also likely to be relevant for IPD.

Comparison of the genome of the Gram negative human pathogen Haemophilus quentini MP1 with other Haemophilus species revealed that, although it is more closely related to Haemophilus haemolyticus than Haemophilus influenzae, the pathogen is in fact genetically distinct, a finding confirmed by phylogenetic analysis using the H. influenzae multilocus sequence typing genes. Further comparison with two other H. quentini strains recently identified in Canada revealed that these three genomes are more closely related than any other Haemophilus species, however there is still some sequence variation. There was no evidence of acquired antimicrobial resistance within the H. quentini MP1 genome nor any mutations within the DNA gyrase or topoisomerase IV genes known to confer resistance to fluoroquinolones, which has been previously identified in other H. quentini isolates. We hope by presenting the annotation and genetic comparison of the H. quentini MP1 genome it will aid the future molecular detection of this potentially emerging pathogen via the identification of unique genes that differentiate it from other Haemophilus species.

Clostridium difficile is a Gram-positive and sporulating enteropathogen that is a major cause of healthcare-associated infections. Even though a large number of genomes of this species have been sequenced, only a few plasmids have been described in the literature. Here, we use a combination of in silico analyses and laboratory experiments to show that plasmids are common in C. difficile. We focus on a group of plasmids that share similarity with the plasmid pCD630, from the reference strain 630. The family of pCD630-like plasmids is defined by the presence of a conserved putative helicase that is likely part of the plasmid replicon. This replicon is compatible with at least some other C. difficile replicons, as strains can carry pCD630-like plasmids in addition to other plasmids. We find two distinct sub-groups of pCD630-like plasmids that differ in size and accessory modules. This study is the first to describe a family of plasmids in C. difficile.

The Plasmodium falciparum variant surface antigen PfEMP1 expressed on the surface of infected erythrocytes is thought to play a major role in the pathology of severe malaria. As the sequence pool of the var genes encoding PfEMP1 expands there are opportunities, despite the high degree of sequence diversity demonstrated by this gene family, to reconstruct full-length var genes from small sequence tags generated from patient isolates. To test whether this is possible we have used a set of recently laboratory adapted ICAM-1-binding parasite isolates to generate sequence tags and, from these, to identify the full-length PfEMP1 being expressed by them. In a subset of the strains available we were able to produce validated, full-length var gene sequences and use these to conduct biophysical analyses of the ICAM-1 binding regions.

Background Culex quinquefasciatus collected in Uganda, where no vector control interventions directly targeting this species have been conducted, was used as a model to determine if it is possible to detect heterogeneities in selection pressure driven by insecticide application targeting other insect species. Methodology/Principal findings Population genetic structure was assessed through microsatellite analysis, and the impact of insecticide pressure by genotyping two target-site mutations, Vgsc-1014F of the voltage-gated sodium channel target of pyrethroid and DDT insecticides, and Ace1-119S of the acetylcholinesterase gene, target of carbamate and organophosphate insecticides. No significant differences in genetic diversity were observed among populations by microsatellite markers with HE ranging from 0.597 to 0.612 and low, but significant, genetic differentiation among populations (FST = 0.019, P = 0.001). By contrast, the insecticide-resistance markers display heterogeneous allelic distributions with significant differences detected between Central Ugandan (urban) populations relative to Eastern and Southwestern (rural) populations. In the central region, a frequency of 62% for Vgsc-1014F, and 32% for the Ace1-119S resistant allele were observed. Conversely, in both Eastern and Southwestern regions the Vgsc-1014F alleles were close to fixation, whilst Ace1-119S allele frequency was 12% (although frequencies may be underestimated due to copy number variation at both loci). Conclusions/Significance Taken together, the microsatellite and both insecticide resistance target-site markers provide evidence that in the face of intense gene flow among populations, disjunction in resistance frequencies arise due to intense local selection pressures despite an absence of insecticidal control interventions targeting Culex., Author summary Culex quinquefasciatus is a primary vector of arboviruses such as West Nile virus (WNV) and St. Louis encephalitis virus (SLEV) in temperate regions, and in many tropical/sub-tropical regions is implicated in the transmission of lymphatic filariasis (LF). Insecticide-based approaches are one of the most important interventions to mitigate disease burden; nevertheless, increased resistance of vectors to insecticides poses a challenge for sustainability and effectiveness of both current and future vector control interventions. Here, Cx. quinquefasciatus collected across a transect from Eastern to South Western Uganda were utilized to infer the likely pattern of the evolution of insecticide resistance among populations through a combination of putative selective loci (the target-site mutations Vgsc-1014F and Ace1-119S) and 26 neutral microsatellite markers. Taken together, the results yielded by these markers provide evidence that Ugandan Cx. quinquefasciatus mosquitoes are under heterogeneous selection pressure imposed by insecticides from distinct classes.