- Volume 6, Issue 9, 2020
Volume 6, Issue 9, 2020
- Research Article
-
- Microbial Evolution and Epidemiology
- Population Genomics
-
-
Diversification in immunogenicity genes caused by selective pressures in invasive meningococci
We studied population genomics of 486 Neisseria meningitidis isolates causing meningitis in the Netherlands during the period 1979–2003 and 2006–2013 using whole-genome sequencing to evaluate the impact of a hyperendemic period of serogroup B invasive disease. The majority of serogroup B isolates belonged to ST-41/44 (41 %) and ST-32 complex (16 %). Comparing the time periods, before and after the decline of serogroup B invasive disease, there was a decrease of ST-41/44 complex sequences (P=0.002). We observed the expansion of a sub-lineage within ST-41/44 complex sequences being associated with isolation from the 1979–2003 time period (P=0.014). Isolates belonging to this sub-lineage expansion within ST-41/44 complex were marked by four antigen allele variants. Presence of these allele variants was associated with isolation from the 1979–2003 time period after correction for multiple testing (Wald test, P=0.0043 for FetA 1–5; P=0.0035 for FHbp 14; P=0.012 for PorA 7–2.4 and P=0.0031 for NHBA two peptide allele). These sequences were associated with 4CMenB vaccine coverage (Fisher’s exact test, P<0.001). Outside of the sub-lineage expansion, isolates with markedly lower levels of predicted vaccine coverage clustered in phylogenetic groups showing a trend towards isolation in the 2006–2013 time period (P=0.08). In conclusion, we show the emergence and decline of a sub-lineage expansion within ST-41/44 complex isolates concurrent with a hyperendemic period in meningococcal meningitis. The expansion was marked by specific antigen peptide allele combinations. We observed preliminary evidence for decreasing 4CMenB vaccine coverage in the post-hyperendemic period.
-
- Mechanisms of Evolution
-
-
Virulence genes and previously unexplored gene clusters in four commensal Neisseria spp. isolated from the human throat expand the neisserial gene repertoire
Alan Calder, Chukwuma Jude Menkiti, Aylin Çağdaş, Jefferson Lisboa Santos, Ricarda Streich, Alice Wong, Amir H. Avini, Ebrima Bojang, Karththeepan Yogamanoharan, Nivetha Sivanesan, Besma Ali, Mariam Ashrafi, Abdirizak Issa, Tajinder Kaur, Aisha Latif, Hani A. Sheik Mohamed, Atifa Maqsood, Laxmi Tamang, Emily Swager, Alex J. Stringer and Lori A.S. SnyderCommensal non-pathogenic Neisseria spp. live within the human host alongside the pathogenic Neisseria meningitidis and Neisseria gonorrhoeae and due to natural competence, horizontal gene transfer within the genus is possible and has been observed. Four distinct Neisseria spp. isolates taken from the throats of two human volunteers have been assessed here using a combination of microbiological and bioinformatics techniques. Three of the isolates have been identified as Neisseria subflava biovar perflava and one as Neisseria cinerea . Specific gene clusters have been identified within these commensal isolate genome sequences that are believed to encode a Type VI Secretion System, a newly identified CRISPR system, a Type IV Secretion System unlike that in other Neisseria spp., a hemin transporter, and a haem acquisition and utilization system. This investigation is the first to investigate these systems in either the non-pathogenic or pathogenic Neisseria spp. In addition, the N. subflava biovar perflava possess previously unreported capsule loci and sequences have been identified in all four isolates that are similar to genes seen within the pathogens that are associated with virulence. These data from the four commensal isolates provide further evidence for a Neisseria spp. gene pool and highlight the presence of systems within the commensals with functions still to be explored.
-
- Communicable Disease Genomics
-
-
Epidemiological and molecular characterization of a novel adenovirus of squirrel monkeys after fatal infection during immunosuppression
Adenoviruses are a frequent cause of acute upper respiratory tract infections that can also cause disseminated disease in immunosuppressed patients. We identified a novel adenovirus, squirrel monkey adenovirus 1 (SqMAdV-1), as the cause of fatal infection in an immunocompromised squirrel monkey (Saimiri boliviensis) at the Keeling Center for Comparative Medicine and Research (KCCMR). Sequencing of SqMAdV-1 revealed that it is most closely related (80.4 % pairwise nucleotide identity) to the titi monkey (Plecturocebus cupreus) adenovirus (TMAdV). Although identified in the titi monkey, TMAdV is highly lethal in these monkeys, and they are not thought to be the natural host. While SqMAdV-1 is similar to other primate adenoviruses in size and genomic characteristics, a nucleotide polymorphism at the expected stop codon of the DNA polymerase gene results in a 126 amino acid extension at the carboxy terminus, a feature not previously observed among other primate adenoviruses. PCR testing and partial sequencing of 95 archived faecal samples from other squirrel monkeys (Saimiri boliviensis and Saimiri sciureus) housed at the KCCMR revealed the presence of three distinct, and apparently endemic species of adenoviruses. A grouping of ten squirrel monkey adenovirus variants has high similarity to SqMAdV-1. A single adenovirus variant (designated SqMAdV-3), detected in five monkeys, has similarity to tufted capuchin (Sapajus apella) adenoviruses. The largest group of adenovirus variants detected (designated SqMAdV-2.0–2.16) has very high similarity (93–99 %) to the TMAdV, suggesting that squirrel monkeys may be the natural host of the TMAdV.
-
- Zoonosis/Anthroponosis
-
-
Comparing serotyping with whole-genome sequencing for subtyping of non-typhoidal Salmonella enterica: a large-scale analysis of 37 serotypes with a public health impact in the USA
More LessSerotyping has traditionally been used for subtyping of non-typhoidal Salmonella (NTS) isolates. However, its discriminatory power is limited, which impairs its use for epidemiological investigations of source attribution. Whole-genome sequencing (WGS) analysis allows more accurate subtyping of strains. However, because of the relative newness and cost of routine WGS, large-scale studies involving NTS WGS are still rare. We aimed to revisit the big picture of subtyping NTS with a public health impact by using traditional serotyping (i.e. reaction between antisera and surface antigens) and comparing the results with those obtained using WGS. For this purpose, we analysed 18 282 sequences of isolates belonging to 37 serotypes with a public health impact that were recovered in the USA between 2006 and 2017 from multiple sources, and were available at the National Center for Biotechnology Information (NCBI). Phylogenetic trees were reconstructed for each serotype using the core genome for the identification of genetic subpopulations. We demonstrated that WGS-based subtyping allows better identification of sources potentially linked with human infection and emerging subpopulations, along with providing information on the risk of dissemination of plasmids and acquired antimicrobial resistance genes (AARGs). In addition, by reconstructing a phylogenetic tree with representative isolates from all serotypes (n=370), we demonstrated genetic variability within and between serotypes, which formed monophyletic, polyphyletic and paraphyletic clades. Moreover, we found (in the entire data set) an increased detection rate for AARGs linked to key antimicrobials (such as quinolones and extended-spectrum cephalosporins) over time. The outputs of this large-scale analysis reveal new insights into the genetic diversity within and between serotypes; the polyphyly and paraphyly of certain serotypes may suggest that the subtyping of NTS to serotypes may not be sufficient. Moreover, the results and the methods presented here, leading to differentiation between genetic subpopulations based on their potential risk to public health, as well as narrowing down the possible sources of these infections, may be used as a baseline for subtyping of future NTS infections and help efforts to mitigate and prevent infections in the USA and globally.
-
-
-
Genomic diversity of Escherichia coli isolates from non-human primates in the Gambia
Increasing contact between humans and non-human primates provides an opportunity for the transfer of potential pathogens or antimicrobial resistance between host species. We have investigated genomic diversity and antimicrobial resistance in Escherichia coli isolates from four species of non-human primates in the Gambia: Papio papio (n=22), Chlorocebus sabaeus (n=14), Piliocolobus badius (n=6) and Erythrocebus patas (n=1). We performed Illumina whole-genome sequencing on 101 isolates from 43 stools, followed by nanopore long-read sequencing on 11 isolates. We identified 43 sequence types (STs) by the Achtman scheme (ten of which are novel), spanning five of the eight known phylogroups of E. coli . The majority of simian isolates belong to phylogroup B2 – characterized by strains that cause human extraintestinal infections – and encode factors associated with extraintestinal disease. A subset of the B2 strains (ST73, ST681 and ST127) carry the pks genomic island, which encodes colibactin, a genotoxin associated with colorectal cancer. We found little antimicrobial resistance and only one example of multi-drug resistance among the simian isolates. Hierarchical clustering showed that simian isolates from ST442 and ST349 are closely related to isolates recovered from human clinical cases (differences in 50 and 7 alleles, respectively), suggesting recent exchange between the two host species. Conversely, simian isolates from ST73, ST681 and ST127 were distinct from human isolates, while five simian isolates belong to unique core-genome ST complexes – indicating novel diversity specific to the primate niche. Our results are of planetary health importance, considering the increasing contact between humans and wild non-human primates.
-
- Microbe-Niche Interactions
- Mutualism, Commensalism and Parasitism
-
-
Phylogenetic distribution and evolutionary dynamics of nod and T3SS genes in the genus Bradyrhizobium
More LessBradyrhizobium are abundant soil bacteria and the major symbiont of legumes. The recent availability of Bradyrhizobium genome sequences provides a large source of information for analysis of symbiotic traits. In this study, we investigated the evolutionary dynamics of the nodulation genes (nod) and their relationship with the genes encoding type III secretion systems (T3SS) and their effectors among bradyrhizobia. Based on the comparative analysis of 146 Bradyrhizobium genome sequences, we identified six different types of T3SS gene clusters. The two predominant cluster types are designated RhcIa and RhcIb and both belong to the RhcI-T3SS family previously described in other rhizobia. They are found in 92/146 strains, most of them also containing nod genes. RhcIa and RhcIb gene clusters differ in the genes they carry: while the translocon-encoding gene nopX is systematically found in strains containing RhcIb, the nopE and nopH genes are specifically conserved in strains containing RhcIa, suggesting that these last two genes might functionally substitute nopX and play a role related to effector translocation. Phylogenetic analysis suggests that bradyrhizobia simultaneously gained nod and RhcI-T3SS gene clusters via horizontal transfer or subsequent vertical inheritance of a symbiotic island containing both. Sequence similarity searches for known Nop effector proteins in bradyrhizobial proteomes revealed the absence of a so-called core effectome, i.e. that no effector is conserved among all Bradyrhizobium strains. However, NopM and SUMO proteases were found to be the main effector families, being represented in the majority of the genus. This study indicates that bradyrhizobial T3SSs might play a more significant symbiotic role than previously thought and provides new candidates among T3SS structural proteins and effectors for future functional investigations.
-
-
-
Symbiosis islands of Loteae-nodulating Mesorhizobium comprise three radiating lineages with concordant nod gene complements and nodulation host-range groupings
Mesorhizobium is a genus of soil bacteria, some isolates of which form an endosymbiotic relationship with diverse legumes of the Loteae tribe. The symbiotic genes of these mesorhizobia are generally carried on integrative and conjugative elements termed symbiosis islands (ICESyms). Mesorhizobium strains that nodulate Lotus spp. have been divided into host-range groupings. Group I (GI) strains nodulate L. corniculatus and L. japonicus ecotype Gifu, while group II (GII) strains have a broader host range, which includes L. pedunculatus. To identify the basis of this extended host range, and better understand Mesorhizobium and ICESym genomics, the genomes of eight Mesorhizobium strains were completed using hybrid long- and short-read assembly. Bioinformatic comparison with previously sequenced mesorhizobia genomes indicated host range was not predicted by Mesorhizobium genospecies but rather by the evolutionary relationship between ICESym symbiotic regions. Three radiating lineages of Loteae ICESyms were identified on this basis, which correlate with Lotus spp. host-range grouping and have lineage-specific nod gene complements. Pangenomic analysis of the completed GI and GII ICESyms identified 155 core genes (on average 30.1 % of a given ICESym). Individual GI or GII ICESyms carried diverse accessory genes with an average of 34.6 % of genes unique to a given ICESym. Identification and comparative analysis of NodD symbiotic regulatory motifs – nod boxes – identified 21 branches across the NodD regulons. Four of these branches were associated with seven genes unique to the five GII ICESyms. The nod boxes preceding the host-range gene nodZ in GI and GII ICESyms were disparate, suggesting regulation of nodZ may differ between GI and GII ICESyms. The broad host-range determinant(s) of GII ICESyms that confer nodulation of L. pedunculatus are likely present amongst the 53 GII-unique genes identified.
-
- Pathogenesis
-
-
Acquisition of plasmids conferring carbapenem and aminoglycoside resistance and loss of surface-exposed macromolecule structures as strategies for the adaptation of Acinetobacter baumannii CC104O/CC15P strains to the clinical setting
Acinetobacter baumannii (Aba) is an emerging opportunistic pathogen associated to nosocomial infections. The rapid increase in multidrug resistance (MDR) among Aba strains underscores the urgency of understanding how this pathogen evolves in the clinical environment. We conducted here a whole-genome sequence comparative analysis of three phylogenetically and epidemiologically related MDR Aba strains from Argentinean hospitals, assigned to the CC104O/CC15P clonal complex. While the Ab244 strain was carbapenem-susceptible, Ab242 and Ab825, isolated after the introduction of carbapenem therapy, displayed resistance to these last resource β-lactams. We found a high chromosomal synteny among the three strains, but significant differences at their accessory genomes. Most importantly, carbapenem resistance in Ab242 and Ab825 was attributed to the acquisition of a Rep_3 family plasmid carrying a bla OXA-58 gene. Other differences involved a genomic island carrying resistance to toxic compounds and a Tn10 element exclusive to Ab244 and Ab825, respectively. Also remarkably, 44 insertion sequences (ISs) were uncovered in Ab825, in contrast with the 14 and 11 detected in Ab242 and Ab244, respectively. Moreover, Ab825 showed a higher killing capacity as compared to the other two strains in the Galleria mellonella infection model. A search for virulence and persistence determinants indicated the loss or IS-mediated interruption of genes encoding many surface-exposed macromolecules in Ab825, suggesting that these events are responsible for its higher relative virulence. The comparative genomic analyses of the CC104O/CC15P strains conducted here revealed the contribution of acquired mobile genetic elements such as ISs and plasmids to the adaptation of A. baumannii to the clinical setting.
-
- Host Adaptation
-
-
Insertion sequences drive the emergence of a highly adapted human pathogen
Pseudomonas aeruginosa is a highly adaptive opportunistic pathogen that can have serious health consequences in patients with lung disorders. Taxonomic outliers of P. aeruginosa of environmental origin have recently emerged as infectious for humans. Here, we present the first genome-wide analysis of an isolate that caused fatal haemorrhagic pneumonia. In two clones, CLJ1 and CLJ3, sequentially recovered from a patient with chronic pulmonary disease, insertion of a mobile genetic element into the P. aeruginosa chromosome affected major virulence-associated phenotypes and led to increased resistance to the antibiotics used to combat the infection. Comparative genome, proteome and transcriptome analyses revealed that this ISL3-family insertion sequence disrupted the genes for flagellar components, type IV pili, O-specific antigens, translesion polymerase and enzymes producing hydrogen cyanide. Seven-fold more insertions were detected in the later isolate, CLJ3, than in CLJ1, some of which modified strain susceptibility to antibiotics by disrupting the genes for the outer-membrane porin OprD and the regulator of β-lactamase expression AmpD. In the Galleria mellonella larvae model, the two strains displayed different levels of virulence, with CLJ1 being highly pathogenic. This study revealed insertion sequences to be major players in enhancing the pathogenic potential of a P. aeruginosa taxonomic outlier by modulating both its virulence and its resistance to antimicrobials, and explains how this bacterium adapts from the environment to a human host.
-
-
-
Complete genome of a unicellular parasite (Antonospora locustae) and transcriptional interactions with its host locust
Microsporidia are a large group of unicellular parasites that infect insects and mammals. The simpler life cycle of microsporidia in insects provides a model system for understanding their evolution and molecular interactions with their hosts. However, no complete genome is available for insect-parasitic microsporidian species. The complete genome of Antonospora locustae, a microsporidian parasite that obligately infects insects, is reported here. The genome size of A. locustae is 3 170 203 nucleotides, composed of 17 chromosomes onto which a total of 1857 annotated genes have been mapped and detailed. A unique feature of the A. locustae genome is the presence of an ultra-low GC region of approximately 25 kb on 16 of the 17 chromosomes, in which the average GC content is only 20 %. Transcription profiling indicated that the ultra-low GC region of the parasite could be associated with differential regulation of host defences in the fat body to promote the parasite’s survival and propagation. Phylogenetic gene analysis showed that A. locustae, and the microsporidian family in general, is likely at an evolutionarily transitional position between prokaryotes and eukaryotes, and that it evolved independently. Transcriptomic analysis showed that A. locustae can systematically inhibit the locust phenoloxidase PPO, TCA and glyoxylate cycles, and PPAR pathways to escape melanization, and can activate host energy transfer pathways to support its reproduction in the fat body, which is an insect energy-producing organ. Our study provides a platform and model for studies of the molecular mechanisms of microsporidium–host interactions in an energy-producing organ and for understanding the evolution of microsporidia.
-
- Systems Microbiology
- Transcriptomics, Proteomics, Networks
-
-
Dual transcriptome analysis reveals differential gene expression modulation influenced by Leishmania arginase and host genetic background
The outcome of Leishmania infection is strongly influenced by the host’s genetic background. BALB/c mice are susceptible to Leishmania infection, while C57BL/6 mice show discrete resistance. Central to the fate of the infection is the availability of l-arginine and the related metabolic processes in the host and parasite. Depending on l-arginine availability, nitric oxide synthase 2 (NOS2) of the host cell produces nitric oxide (NO) controlling the parasite growth. On the other hand, Leishmania can also use host l-arginine for the production of polyamines through its own arginase activity, thus favouring parasite replication. Considering RNA-seq data, we analysed the dual modulation of host and parasite gene expression of BALB/c or C57BL/6 mouse bone marrow-derived macrophages (BMDMs) after 4 h of infection with Leishmania amazonensis wild-type (La-WT) or L. amazonensis arginase knockout (La-arg-). We identified 12 641 host transcripts and 8282 parasite transcripts by alignment analysis with the respective Mus musculus and L. mexicana genomes. The comparison of BALB/c_La-arg- versus BALB/c_La-WT revealed 233 modulated transcripts, with most related to the immune response and some related to the amino acid transporters and l-arginine metabolism. In contrast, the comparison of C57BL/6_La-arg- vs. C57BL/6_La-WT revealed only 30 modulated transcripts, including some related to the immune response but none related to amino acid transport or l-arginine metabolism. The transcriptome profiles of the intracellular amastigote revealed 94 modulated transcripts in the comparison of La-arg-_BALB/c vs. La-WT_BALB/c and 45 modulated transcripts in the comparison of La-arg-_C57BL/6 vs. La-WT_C57BL/6. Taken together, our data present new insights into the impact of parasite arginase activity on the orchestration of the host gene expression modulation, including in the immune response and amino acid transport and metabolism, mainly in susceptible BALB/c-infected macrophages. Moreover, we show how parasite arginase activity affects parasite gene expression modulation, including amino acid uptake and amastin expression.
-
- Pangenome Analysis, Big-data Approaches
-
-
Comparative genomic insights into Yersinia hibernica – a commonly misidentified Yersinia enterocolitica-like organism
Food-associated outbreaks linked to enteropathogenic Yersinia enterocolitica are of concern to public health. Pigs and their meat are recognized risk factors for transmission of Y. enterocolitica . This study aimed to describe the comparative genomics of Y. enterocolitica along with a number of misclassified Yersinia isolates, now constituting the recently described Yersinia hibernica . The latter was originally cultured from an environmental sample taken at a pig slaughterhouse. Unique features were identified in the genome of Y. hibernica, including a novel integrative conjugative element (ICE), denoted as ICE Yh-1 contained within a 255 kbp region of plasticity. In addition, a zebrafish embryo infection model was adapted and applied to assess the virulence potential among Yersinia isolates including Y. hibernica .
-