Population dynamics of Neisseria gonorrhoeae in Shanghai, China: a comparative study

Investigation of cross-regional spread and evolution of equine influenza H3N8 at US and global scales using Bayesian phylogeography based on balanced subsampling

Equine influenza virus (EIV) is a highly contagious pathogen of equids, and a well-known burden in global equine health. EIV H3N8 variants seasonally emerged and resulted in EIV outbreaks in the United States (US) and worldwide. The present study evaluated the pattern of cross-regional EIV H3N8 spread and evolutionary characteristics at US and global scales using Bayesian phylogeography with balanced subsampling based on regional horse population size. A total of 297 Haemagglutinin (HA) sequences of global EIV H3N8 were collected from 1963 to 2019 and subsampled to global subset (n = 67), raw US sequences (n = 100) and US subset (n = 44) datasets. Discrete trait phylogeography analysis was used to estimate the transmission history of EIV using four global and US genome datasets. The North American lineage was the major source of globally dominant EIV variants and spread to other global regions. The US EIV strains generally spread from the southern and midwestern regions to other regions. The EIV H3N8 accumulated approximately three nucleotide substitutions per year in the HA gene under heterogenous local positive selection. Our findings will guide better decision making of target intervention strategies of EIV H3N8 infection and provide the better scheme of genomic surveillance in the US and global equine health. This article is protected by copyright. All rights reserved.

Microbial sequence typing in the genomic era

Next-generation sequencing (NGS), also known as high-throughput sequencing, is changing the field of microbial genomics research. NGS allows for a more comprehensive analysis of the diversity, structure and composition of microbial genes and genomes compared to the traditional automated Sanger capillary sequencing at a lower cost. NGS strategies have expanded the versatility of standard and widely used typing approaches based on nucleotide variation in several hundred DNA sequences and a few gene fragments (MLST, MLVA, rMLST and cgMLST). NGS can now accommodate variation in thousands or millions of sequences from selected amplicons to full genomes (WGS, NGMLST and HiMLST). To extract signals from high-dimensional NGS data and make valid statistical inferences, novel analytic and statistical techniques are needed. In this review, we describe standard and new approaches for microbial sequence typing at gene and genome levels and guidelines for subsequent analysis, including methods and computational frameworks. We also present several applications of these approaches to some disciplines, namely genotyping, phylogenetics and molecular epidemiology.

The Impact of Selection, Gene Conversion, and Biased Sampling on the Assessment of Microbial Demography

Recent studies have linked demographic changes and epidemiological patterns in bacterial populations using coalescent-based approaches. We identified 26 studies using skyline plots and found that 21 inferred overall population expansion. This surprising result led us to analyze the impact of natural selection, recombination (gene conversion), and sampling biases on demographic inference using skyline plots and site frequency spectra (SFS). Forward simulations based on biologically relevant parameters from Escherichia coli populations showed that theoretical arguments on the detrimental impact of recombination and especially natural selection on the reconstructed genealogies cannot be ignored in practice. In fact, both processes systematically lead to spurious interpretations of population expansion in skyline plots (and in SFS for selection). Weak purifying selection, and especially positive selection, had important effects on skyline plots, showing patterns akin to those of population expansions. State-of-the-art techniques to remove recombination further amplified these biases. We simulated three common sampling biases in microbiological research: uniform, clustered, and mixed sampling. Alone, or together with recombination and selection, they further mislead demographic inferences producing almost any possible skyline shape or SFS. Interestingly, sampling sub-populations also affected skyline plots and SFS, because the coalescent rates of populations and their sub-populations had different distributions. This study suggests that extreme caution is needed to infer demographic changes solely based on reconstructed genealogies. We suggest that the development of novel sampling strategies and the joint analyzes of diverse population genetic methods are strictly necessary to estimate demographic changes in populations where selection, recombination, and biased sampling are present.

Bacteriocins and other bioactive substances of probiotic lactobacilli as biological weapons against Neisseria gonorrhoeae

In the search of new antimicrobial agents against Neisseria gonorrhoeae, the bacteriocins-producing probiotic lactobacilli deserve special attention. The inhibitory effects of biosubstances such as organic acids, hydrogen peroxide and each bacteriocin-like inhibitory substance (BLIS) L23 and L60 on the growth of different gonococcal strains were investigated. Different non-treated and treated cell-free supernatants of two probiotic lactobacilli containing these metabolites were used. The aims of this work were (i) to evaluate the antimicrobial activity of the biosubstances produced by two probiotic lactobacilli, estimating the proportion in which each of them is responsible for the inhibitory effect, (ii) to define their minimum inhibitory concentrations (MICs) and, (iii) to determine the potential interactions between these biosubstances against N. gonorrhoeae. The main antimicrobial metabolites were the BLIS-es L23 and L60 in comparison with other biosubstances. Proportionally, their contributions to the inhibition on the gonococcal growth were 87.28% and 80.66%, respectively. The MIC values of bacteriocins were promising since these substances, when diluted, showed considerable inhibitory activity for all gonococci. In the interaction between bacteriocins, 100% of synergism was found on the gonococcal growth. In summary, this study indicates that both L23 and L60 could potentially serve to design new bioproducts against N. gonorrhoeae.

Population expansions shared among coexisting bacterial lineages are revealed by genetic evidence

Comparative population studies can help elucidate the influence of historical events upon current patterns of biodiversity among taxa that coexist in a given geographic area. In particular, comparative assessments derived from population genetics and coalescent theory have been used to investigate population dynamics of bacterial pathogens in order to understand disease epidemics. In contrast, and despite the ecological relevance of non-host associated and naturally occurring bacteria, there is little understanding of the processes determining their diversity. Here we analyzed the patterns of genetic diversity in coexisting populations of three genera of bacteria (Bacillus, Exiguobacterium, and Pseudomonas) that are abundant in the aquatic systems of the Cuatro Cienegas Basin, Mexico. We tested the hypothesis that a common habitat leaves a signature upon the genetic variation present in bacterial populations, independent of phylogenetic relationships. We used multilocus markers to assess genetic diversity and (1) performed comparative phylogenetic analyses, (2) described the genetic structure of bacterial populations, (3) calculated descriptive parameters of genetic diversity, (4) performed neutrality tests, and (5) conducted coalescent-based historical reconstructions. Our results show a trend of synchronic expansions across most populations independent of both lineage and sampling site. Thus, we provide empirical evidence supporting the analysis of coexisting bacterial lineages in natural environments to advance our understanding of bacterial evolution beyond medical or health-related microbes.

Patterns of nucleotide diversity of the ldpA circadian gene in closely related species of cyanobacteria from extreme cold deserts

In the circadian system of cyanobacteria, the ldpA gene is a component of the input to the clock. We comparatively analyzed nucleotide polymorphism of this gene in populations of two closely related species of cyanobacteria (denoted as Synechococcus species S1 and S2, respectively) from extreme cold deserts in Antarctica, the Canadian Arctic, and Tibet. Although both species manifested similarly high haplotype diversities (0.990 and 0.809, respectively), the nucleotide diversity differed significantly (0.0091 in S1 and 0.0037 in S2). The populations of species S2 were more differentiated (F(ST) = 0.2242) compared to those of species S1 (F(ST) between 0.0296 and 0.1188). An analysis of positive selection with several tests yielded highly significant values (P < 0.01) for both species. On the other hand, these results may be somewhat compromised by fluctuating population sizes of the species. The apparent selection pressure coupled with the pronounced demographic factors, such as population expansion, small effective population size, and genetic drift, may thus result in the observed significant interpopulation differentiation and subsequent speciation of cyanobacteria.

Pathogen typing in the genomics era: MLST and the future of molecular epidemiology

Multi-locus sequence typing (MLST) is a high-resolution genetic typing approach to identify species and strains of pathogens impacting human health, agriculture (animals and plants), and biosafety. In this review, we outline the general concepts behind MLST, molecular approaches for obtaining MLST data, analytical approaches for MLST data, and the contributions MLST studies have made in a wide variety of areas. We then look at the future of MLST and their relative strengths and weaknesses with respect to whole genome sequence typing approaches that are moving into the research arena at an ever-increasing pace. Throughout the paper, we provide exemplar references of these various aspects of MLST. The literature is simply too vast to make this review comprehensive, nevertheless, we have attempted to include enough references in a variety of key areas to introduce the reader to the broad applications and complications of MLST data.

Longitudinal analysis of the evolution and dissemination of Neisseria gonorrhoeae strains (Saskatchewan, Canada, 2005 to 2008) reveals three major circulating strains and convergent evolution of ciprofloxacin and azithromycin resistance

A longitudinal study combining multilocus sequence typing with molecular evolutionary analysis determined the distribution, population structure, and evolution of antibiotic resistance in Neisseria gonorrhoeae isolates in Saskatchewan that were collected between 2005 and 2008. Of 195 gonococcal isolates examined, 29 sequence types (STs) were identified with 3 major circulating strains (ST-1 through ST-3) comprising 52% of all gonococcal isolates studied. The prevalences, persistence, distribution patterns, and clonalities of these isolates strongly suggest that gonorrhea endemicity within this broad geographic region was driven by these 3 circulating strains. ST-1 exhibited a significantly (P = 0.001) higher prevalence throughout the study than did the others, accounting for ∼25% of the tested isolates each year. The spatial distributions of the gonococcal strains indicated that ST-1 in 2007 entered a linear component of the sexual network, reaching the remote north and resulting in the further spread and maintenance of infection. Ciprofloxacin and azithromycin resistances were observed in distantly related gonococcal lineages, clearly indicating the convergent acquisition of these antibiotic-resistant phenotypes. In addition, all ciprofloxacin- and azithromycin-resistant lineages were found at the edges of the minimum spanning tree, far from the major lineages, suggesting that these antibiotic phenotypes were most likely introduced into the province. In contrast, resistance to penicillin was found mostly in the endemic gonococcal lineages, suggesting that penicillin resistance was probably acquired in Saskatchewan as a result of spontaneous mutations in already-established lineages. Tetracycline resistance was present in all STs except one, indicating its ubiquitous nature in the gonococcal population studied.

Global fluoroquinolone resistance epidemiology and implictions for clinical use

This paper on the fluoroquinolone resistance epidemiology stratifies the data according to the different prescription patterns by either primary or tertiary caregivers and by indication. Global surveillance studies demonstrate that fluoroquinolone resistance rates increased in the past years in almost all bacterial species except S. pneumoniae and H. influenzae, causing community-acquired respiratory tract infections. However, 10 to 30% of these isolates harbored first-step mutations conferring low level fluoroquinolone resistance. Fluoroquinolone resistance increased in Enterobacteriaceae causing community acquired or healthcare associated urinary tract infections and intraabdominal infections, exceeding 50% in some parts of the world, particularly in Asia. One to two-thirds of Enterobacteriaceae producing extended spectrum β-lactamases were fluoroquinolone resistant too. Furthermore, fluoroquinolones select for methicillin resistance in Staphylococci. Neisseria gonorrhoeae acquired fluoroquinolone resistance rapidly; actual resistance rates are highly variable and can be as high as almost 100%, particularly in Asia, whereas resistance rates in Europe and North America range from <10% in rural areas to >30% in established sexual networks. In general, the continued increase in fluoroquinolone resistance affects patient management and necessitates changes in some guidelines, for example, treatment of urinary tract, intra-abdominal, skin and skin structure infections, and traveller's diarrhea, or even precludes the use in indications like sexually transmitted diseases and enteric fever.

TaqMan real-time quantitative PCR assay for detection of fluoroquinolone-resistant Neisseria gonorrhoeae

It is noted that more than 99 % of fluoroquinolone resistance in Neisseria gonorrhoeae (QRNG) specimens have been shown to have the mutation of Ser91/Phe in the gyrA gene. In order to detect QRNG isolates as quickly as possible, the real-time TaqMan quantitative PCR assay was established for detection of the point mutation of Ser91/Phe in gyrA gene. The standard curve was generated automatically on ABI Prism PE7500. The correlation coefficient (r) of the standard curve was -0.9984 (R(2) = 0.9968), indicating a quietly precise log-linear relationship between the concentration of target DNA and the Ct value. Presently, correlated, cultured antimicrobial susceptibility testing of N. gonorrhoeae isolates continues to be the gold standard method for the detection of antimicrobial resistance. Comparison to the correlated, cultured antimicrobial susceptibility testing, the sensitivity and specificity of the established TaqMan assay for the detection of the QRNG specimens were 100 and 99 %, respectively. The TaqMan assay also allows for rapid detection of QRNG isolates without complex laboratory techniques. Therefore, real-time TaqMan quantitative PCR assay is a rapid, simple, highly sensitive, highly specific, and easy-to-perform method for the detection of the QRNG specimens. It can be applied as a quick screening method for QRNG isolates to help clinical determination of optimal treatment prescription.

Population genomics and phylogeography of an Australian dairy factory derived lytic bacteriophage

In this study, we present the full genomic sequences and evolutionary analyses of a serially sampled population of 28 Lactococcus lactis-infecting phage belonging to the 936-like group in Australia. Genome sizes were consistent with previously available genomes ranging in length from 30.9 to 32.1 Kbp and consisted of 55-65 open reading frames. We analyzed their genetic diversity and found that regions of high diversity are correlated with high recombination rate regions (P value = 0.01). Phylogenetic inference showed two major clades that correlate well with known host range. Using the extended Bayesian Skyline model, we found that population size has remained mostly constant through time. Moreover, the dispersion pattern of these genomes is in agreement with human-driven dispersion as suggested by phylogeographic analysis. In addition, selection analysis found evidence of positive selection on codon positions of the Receptor Binding Protein (RBP). Likewise, positively selected sites in the RBP were located within the neck and head region in the crystal structure, both known determinants of host range. Our study demonstrates the utility of phylogenetic methods applied to whole genome data collected from populations of phage for providing insights into applied microbiology.

The evolution of HIV: inferences using phylogenetics

Molecular phylogenetics has revolutionized the study of not only evolution but also disparate fields such as genomics, bioinformatics, epidemiology, ecology, microbiology, molecular biology and biochemistry. Particularly significant are its achievements in population genetics as a result of the development of coalescent theory, which have contributed to more accurate model-based parameter estimation and explicit hypothesis testing. The study of the evolution of many microorganisms, and HIV in particular, have benefited from these new methodologies. HIV is well suited for such sophisticated population analyses because of its large population sizes, short generation times, high substitution rates and relatively small genomes. All these factors make HIV an ideal and fascinating model to study molecular evolution in real time. Here we review the significant advances made in HIV evolution through the application of phylogenetic approaches. We first examine the relative roles of mutation and recombination on the molecular evolution of HIV and its adaptive response to drug therapy and tissue allocation. We then review some of the fundamental questions in HIV evolution in relation to its origin and diversification and describe some of the insights gained using phylogenies. Finally, we show how phylogenetic analysis has advanced our knowledge of HIV dynamics (i.e., phylodynamics).

Molecular phylodynamics and protein modeling of infectious salmon anemia virus (ISAV)

BACKGROUND

ISAV is a member of the Orthomyxoviridae family that affects salmonids with disastrous results. It was first detected in 1984 in Norway and from then on it has been reported in Canada, United States, Scotland and the Faroe Islands. Recently, an outbreak was recorded in Chile with negative consequences for the local fishing industry. However, few studies have examined available data to test hypotheses associated with the phylogeographic partitioning of the infecting viral population, the population dynamics, or the evolutionary rates and demographic history of ISAV. To explore these issues, we collected relevant sequences of genes coding for both surface proteins from Chile, Canada, and Norway. We addressed questions regarding their phylogenetic relationships, evolutionary rates, and demographic history using modern phylogenetic methods.

RESULTS

A recombination breakpoint was consistently detected in the Hemagglutinin-Esterase (he) gene at either side of the Highly Polymorphic Region (HPR), whereas no recombination breakpoints were detected in Fusion protein (f) gene. Evolutionary relationships of ISAV revealed the 2007 Chilean outbreak group as a monophyletic clade for f that has a sister relationship to the Norwegian isolates. Their tMRCA is consistent with epidemiological data and demographic history was successfully recovered showing a profound bottleneck with further population expansion. Finally, selection analyses detected ongoing diversifying selection in f and he codons associated with protease processing and the HPR region, respectively.

CONCLUSIONS

Our results are consistent with the Norwegian origin hypothesis for the Chilean outbreak clade. In particular, ISAV HPR0 genotype is not the ancestor of all ISAV strains, although SK779/06 (HPR0) shares a common ancestor with the Chilean outbreak clade. Our analyses suggest that ISAV shows hallmarks typical of RNA viruses that can be exploited in epidemiological and surveillance settings. In addition, we hypothesized that genetic diversity of the HPR region is governed by recombination, probably due to template switching and that novel fusion gene proteolytic sites confer a selective advantage for the isolates that carry them. Additionally, protein modeling allowed us to relate the results of phylogenetic studies with the predicted structures. This study demonstrates that phylogenetic methods are important tools to predict future outbreaks of ISAV and other salmon pathogens.

Review and international recommendation of methods for typing neisseria gonorrhoeae isolates and their implications for improved knowledge of gonococcal epidemiology, treatment, and biology

Gonorrhea, which may become untreatable due to multiple resistance to available antibiotics, remains a public health problem worldwide. Precise methods for typing Neisseria gonorrhoeae, together with epidemiological information, are crucial for an enhanced understanding regarding issues involving epidemiology, test of cure and contact tracing, identifying core groups and risk behaviors, and recommending effective antimicrobial treatment, control, and preventive measures. This review evaluates methods for typing N. gonorrhoeae isolates and recommends various methods for different situations. Phenotypic typing methods, as well as some now-outdated DNA-based methods, have limited usefulness in differentiating between strains of N. gonorrhoeae. Genotypic methods based on DNA sequencing are preferred, and the selection of the appropriate genotypic method should be guided by its performance characteristics and whether short-term epidemiology (microepidemiology) or long-term and/or global epidemiology (macroepidemiology) matters are being investigated. Currently, for microepidemiological questions, the best methods for fast, objective, portable, highly discriminatory, reproducible, typeable, and high-throughput characterization are N. gonorrhoeae multiantigen sequence typing (NG-MAST) or full- or extended-length porB gene sequencing. However, pulsed-field gel electrophoresis (PFGE) and Opa typing can be valuable in specific situations, i.e., extreme microepidemiology, despite their limitations. For macroepidemiological studies and phylogenetic studies, DNA sequencing of chromosomal housekeeping genes, such as multilocus sequence typing (MLST), provides a more nuanced understanding.

Testing spatiotemporal hypothesis of bacterial evolution using methicillin-resistant Staphylococcus aureus ST239 genome-wide data within a bayesian framework

Staphylococcus aureus is a common cause of infections that has undergone rapid global spread over recent decades. Formal phylogeographic methods have not yet been applied to the molecular epidemiology of bacterial pathogens because the limited genetic diversity of data sets based on individual genes usually results in poor phylogenetic resolution. Here, we investigated a whole-genome single nucleotide polymorphism (SNP) data set of health care-associated Methicillin-resistant S. aureus sequence type 239 (HA-MRSA ST239) strains, which we analyzed using Markov spatial models that incorporate geographical sampling distributions. The reconstructed timescale indicated a temporal origin of this strain shortly after the introduction of Methicillin, followed by global pandemic spread. The estimate of the temporal origin was robust to the molecular clock, coalescent prior, full/intergenic/synonymous SNP inclusion, and correction for excluded invariant site patterns. Finally, phylogeographic analyses statistically supported the role of human movement in the global dissemination of HA-MRSA ST239, although it was unable to conclusively resolve the location of the root. This study demonstrates that bacterial genomes can indeed contain sufficient evolutionary information to elucidate the temporal and spatial dynamics of transmission. Future applications of this approach to other bacterial strains may provide valuable epidemiological insights that may justify the cost of genome-wide typing.

Predominant porB1A and porB1B genotypes and correlation of gene mutations with drug resistance in Neisseria gonorrhoeae isolates in Eastern China

BACKGROUND

Variations of porB1A and porB1B genes and their serotypes exist in Neisseria gonorrhoeae isolates from different geographical areas, and some site mutations in the porB1B gene correlate with drug resistance.

METHODS

The β-lactamase production of N. gonorrhoeae isolates was determined by paper acidometric test and nitrocefin discs. The porB1A and porB1B genes of 315 non-penicillinase-producting N. gonorrhoeae (non-PPNG) strains were amplified by PCR for sequencing to determine serotypes and site mutations. A duplex PCR was designed to simultaneously detect both porB1A and porB1B genes. Penicillin and tetracycline resistance was assessed by an in vitro drug sensitivity test.

RESULTS

Of the N. gonorrhoeae isolates, 31.1% tested positive for porB1A and 68.9% for porB1B genes. All the 98 porB1A+ isolates belonging to IA6 serotype with either no mutation at the 120 and 121 sites (88.8%) or a D120G (11.2%) mutation and were no resistance to both penicillin and tetracycline. Among the 217 porB1B+ isolates, 26.7%, 22.6% and 11.5% belonged to IB3, IB3/6 and IB4 serotypes, respectively. Particularly, two novel chimeric serotypes, IB3/6-IB2 and IB2-IB4-IB2, were found in 77 and 8 porB1B+ isolates. Two hundred and twelve (97.7%) of the porB1B+ isolates were presented G120 and/or A121 mutations with 163 (76.9%) at both sites. Interestingly, within the 77 porB1B+ isolates belonging to IB3/6-IB2 serotype, 15 were discovered to possess novel deletions at both A121 and N122 sites. All the replacement mutations at these sites in PorB1B were correlated with resistance and the deletion mutation showed the highest resistance.

CONCLUSION

N. gonorrhoeae isolates circulating in Eastern China include a sole PorB1A serotype (IA6) and five PorB1B serotypes. Multiple mutations in porB1B genes, including novel A121 and N122 deletions, are correlated with high levels of penicillin and tetracycline resistance.