THE EFFECT OF ECOLOGICAL DIFFERENTIATION ON GENETIC RECOMBINATION IN THE ENTEROBACTERIA

The existence of distinct species of life is generally explained by the genetic process of reproduction without recombination between populations and/or the ecological process of adaptation to different environments. Both processes affect prokaryotes, and have shaped existing genomes. Here, we use comparative genomic techniques to evaluate the dynamics of divergence among species of the Enterobacteriaceae. Bacteria such as Escherichia coli preferentially acquire allelic variants from closely related organisms (i.e. other E. coli) rather than from more diverged bacteria. Ecological differences between donor and recipient affect the probability of allelic variants becoming fixed across the recombining population. We examine the history of recombination among groups of genomes that no longer recombine with each other, but retain sufficient conservation of ancestral nucleotide sequences to allow recombination to be inferred. From these analyses, we conclude that substantial levels of recombination occurred between E. coli and diverging lineages even after some regions of the genomes had acquired many nucleotide differences. We identify two evolutionary radiations leading to E. coli where the disparity among loci confounds the phylogenetic relationships among species, as evidenced by topological incongruence among gene trees. The forces affecting recombination, reflected in both pairwise divergence and topologically informative sites, vary across regions of the genome measuring tens of kilobases. To examine the relationship between ecological differentiation and genetic recombination, we characterize differences that could be responsible for ecological differentiation among these species. Some of the loci with the most apparent functional differences (i.e. the gain and loss of genes) are associated with the greatest levels of sequence divergence between species, consistent with the hypothesis that ecological divergence interferes with homologous recombination, and therefore drives sequence divergence and genetic isolation. To investigate the role of more subtle ecological differentiation, we develop a statistical framework to evaluate codon usage bias of each protein-coding gene, taking into account the stochastic balance between codon selection, which is driven by the need for high expression, and mutational biases. This tool will be useful in future studies examining codon selection as contribution to diversification among the ecologically diverse species of Enterobacteriaceae

Draft Whole-Genome Sequences of Xylella fastidiosa subsp. fastidiosa Strains TPD3 and TPD4, Isolated from Grapevines in Hou-li, Taiwan.

We report the draft assemblies of TPD3 and TPD4, two Xylella fastidiosa subsp. fastidiosa isolates infecting grapevines in Hou-li, Taiwan. TPD3 and TPD4 showed similar characteristics regarding genome size (2,483,503 bp and 2,491,539 bp, respectively), GC content (51.49% and 51.47%, respectively), and number of protein-coding sequences (2,394 and 2,413, respectively)

Heteroresistance to the model antimicrobial peptide polymyxin B in the emerging Neisseria meningitidis lineage 11.2 urethritis clade: mutations in the pilMNOPQ operon

Clusters of Neisseria meningitidis (Nm) urethritis among primarily heterosexual males in multiple US cities have been attributed to a unique non‐encapsulated meningococcal clade (the US Nm urethritis clade, US_NmUC) within the hypervirulent clonal complex 11. Resistance to antimicrobial peptides (AMPs) is a key feature of urogenital pathogenesis of the closely related species, Neisseria gonorrhoeae. The US_NmUC isolates were found to be highly resistant to the model AMP, polymyxin B (PmB, MICs 64–256 µg ml–1). The isolates also demonstrated stable subpopulations of heteroresistant colonies that showed near total resistant to PmB (MICs 384–1024 µg ml–1) and colistin (MIC 256 µg ml–1) as well as enhanced LL‐37 resistance. This is the first observation of heteroresistance in N. meningitidis. Consistent with previous findings, overall PmB resistance in US_NmUC isolates was due to active Mtr efflux and LptA‐mediated lipid A modification. However, whole genome sequencing, variant analyses and directed mutagenesis revealed that the heteroresistance phenotypes and very high‐level AMP resistance were the result of point mutations and IS1655 element movement in the pilMNOPQ operon, encoding the type IV pilin biogenesis apparatus. Cross‐resistance to other classes of antibiotics was also observed in the heteroresistant colonies. High‐level resistance to AMPs may contribute to the pathogenesis of US_NmUC

Transmission of SARS-CoV-2 in free-ranging white-tailed deer in the United States

SARS-CoV-2 is a zoonotic virus with documented bi-directional transmission between people and animals. Transmission of SARS-CoV-2 from humans to free-ranging white-tailed deer (Odocoileus virginianus) poses a unique public health risk due to the potential for reservoir establishment where variantsmay persist and evolve. We collected 8,830 respiratory samples from free-ranging white-tailed deer across Washington, D.C. and 26 states in the United States between November 2021 and April 2022. We obtained 391 sequences and identified 34 Pango lineages including the Alpha, Gamma, Delta, and Omicron variants. Evolutionary analyses showed these white-tailed deer viruses originated fromat least 109 independent spillovers fromhumans,which resulted in 39 cases of subsequent local deer-to-deer transmission and three cases of potential spillover from white-tailed deer back to humans. Viruses repeatedly adapted to white-tailed deer with recurring amino acid substitutions across spike and other proteins. Overall, our findings suggest that multiple SARS-CoV- 2 lineages were introduced, became enzootic, and co-circulated in whitetailed deer

Household Transmission and Symptomology of Severe Acute Respiratory Syndrome Coronavirus 2 Alpha Variant among Children-California and Colorado, 2021.

OBJECTIVE: To assess the household secondary infection risk (SIR) of B.1.1.7 (Alpha) and non-Alpha lineages of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) among children. STUDY DESIGN: During January to April 2021, we prospectively followed households with a SARS-CoV-2 infection. We collected questionnaires, serial nasopharyngeal swabs for reverse transcription polymerase chain reaction testing and whole genome sequencing, and serial blood samples for serology testing. We calculated SIRs by primary case age (pediatric vs adult), household contact age, and viral lineage. We evaluated risk factors associated with transmission and described symptom profiles among children. RESULTS: Among 36 households with pediatric primary cases, 21 (58%) had secondary infections. Among 91 households with adult primary cases, 51 (56%) had secondary infections. SIRs among pediatric and adult primary cases were 45% and 54%, respectively (OR, 0.79; 95% CI, 0.41-1.54). SIRs among pediatric primary cases with Alpha and non-Alpha lineage were 55% and 46%, respectively (OR, 1.52; 95% CI, 0.51-4.53). SIRs among pediatric and adult household contacts were 55% and 49%, respectively (OR, 1.01; 95% CI, 0.68-1.50). Among pediatric contacts, no significant differences in the odds of acquiring infection by demographic or household characteristics were observed. CONCLUSIONS: Household transmission of SARS-CoV-2 from children and adult primary cases to household members was frequent. The risk of secondary infection was similar among child and adult household contacts. Among children, household transmission of SARS-CoV-2 and the risk of secondary infection was not influenced by lineage. Continued mitigation strategies (eg, masking, physical distancing, vaccination) are needed to protect at-risk groups regardless of virus lineage circulating in communities

Household Transmission of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Alpha Variant-United States, 2021.

BACKGROUND: In Spring 2021, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) B.1.1.7 (Alpha) became the predominant variant in the United States. Research suggests that Alpha has increased transmissibility compared with non-Alpha lineages. We estimated household secondary infection risk (SIR), assessed characteristics associated with transmission, and compared symptoms of persons with Alpha and non-Alpha infections. METHODS: We followed households with SARS-CoV-2 infection for 2 weeks in San Diego County and metropolitan Denver, January to April 2021. We collected epidemiologic information and biospecimens for serology, reverse transcription-polymerase chain reaction (RT-PCR), and whole-genome sequencing. We stratified SIR and symptoms by lineage and identified characteristics associated with transmission using generalized estimating equations. RESULTS: We investigated 127 households with 322 household contacts; 72 households (56.7%) had member(s) with secondary infections. SIRs were not significantly higher for Alpha (61.0% [95% confidence interval, 52.4-69.0%]) than non-Alpha (55.6% [44.7-65.9%], P = .49). In households with Alpha, persons who identified as Asian or Hispanic/Latino had significantly higher SIRs than those who identified as White (P = .01 and .03, respectively). Close contact (eg, kissing, hugging) with primary cases was associated with increased transmission for all lineages. Persons with Alpha infection were more likely to report constitutional symptoms than persons with non-Alpha (86.9% vs 76.8%, P = .05). CONCLUSIONS: Household SIRs were similar for Alpha and non-Alpha. Comparable SIRs may be due to saturation of transmission risk in households due to extensive close contact, or true lack of difference in transmission rates. Avoiding close contact within households may reduce SARS-CoV-2 transmission for all lineages among household members

Complex Evolutionary History of the Aeromonas veronii Group Revealed by Host Interaction and DNA Sequence Data

Aeromonas veronii biovar sobria, Aeromonas veronii biovar veronii, and Aeromonas allosaccharophila are a closely related group of organisms, the Aeromonas veronii Group, that inhabit a wide range of host animals as a symbiont or pathogen. In this study, the ability of various strains to colonize the medicinal leech as a model for beneficial symbiosis and to kill wax worm larvae as a model for virulence was determined. Isolates cultured from the leech out-competed other strains in the leech model, while most strains were virulent in the wax worms. Three housekeeping genes, recA, dnaJ and gyrB, the gene encoding chitinase, chiA, and four loci associated with the type three secretion system, ascV, ascFG, aexT, and aexU were sequenced. The phylogenetic reconstruction failed to produce one consensus tree that was compatible with most of the individual genes. The Approximately Unbiased test and the Genetic Algorithm for Recombination Detection both provided further support for differing evolutionary histories among this group of genes. Two contrasting tests detected recombination within aexU, ascFG, ascV, dnaJ, and gyrB but not in aexT or chiA. Quartet decomposition analysis indicated a complex recent evolutionary history for these strains with a high frequency of horizontal gene transfer between several but not among all strains. In this study we demonstrate that at least for some strains, horizontal gene transfer occurs at a sufficient frequency to blur the signal from vertically inherited genes, despite strains being adapted to distinct niches. Simply increasing the number of genes included in the analysis is unlikely to overcome this challenge in organisms that occupy multiple niches and can exchange DNA between strains specialized to different niches. Instead, the detection of genes critical in the adaptation to specific niches may help to reveal the physiological specialization of these strains

Quantification of codon selection for comparative bacterial genomics

<p>Abstract</p> <p>Background</p> <p>Statistics measuring codon selection seek to compare genes by their sensitivity to selection for translational efficiency, but existing statistics lack a model for testing the significance of differences between genes. Here, we introduce a new statistic for measuring codon selection, the Adaptive Codon Enrichment (ACE).</p> <p>Results</p> <p>This statistic represents codon usage bias in terms of a probabilistic distribution, quantifying the extent that preferred codons are over-represented in the gene of interest relative to the mean and variance that would result from stochastic sampling of codons. Expected codon frequencies are derived from the observed codon usage frequencies of a broad set of genes, such that they are likely to reflect nonselective, genome wide influences on codon usage (<it>e.g</it>. mutational biases). The relative adaptiveness of synonymous codons is deduced from the frequency of codon usage in a pre-selected set of genes relative to the expected frequency. The ACE can predict both transcript abundance during rapid growth and the rate of synonymous substitutions, with accuracy comparable to or greater than existing metrics. We further examine how the composition of reference gene sets affects the accuracy of the statistic, and suggest methods for selecting appropriate reference sets for any genome, including bacteriophages. Finally, we demonstrate that the ACE may naturally be extended to quantify the genome-wide influence of codon selection in a manner that is sensitive to a large fraction of codons in the genome. This reveals substantial variation among genomes, correlated with the tRNA gene number, even among groups of bacteria where previously proposed whole-genome measures show little variation.</p> <p>Conclusions</p> <p>The statistical framework of the ACE allows rigorous comparison of the level of codon selection acting on genes, both within a genome and between genomes.</p

BMScan: using whole genome similarity to rapidly and accurately identify bacterial meningitis causing species

Abstract Background Bacterial meningitis is a life-threatening infection that remains a public health concern. Bacterial meningitis is commonly caused by the following species: Neisseria meningitidis, Streptococcus pneumoniae, Listeria monocytogenes, Haemophilus influenzae and Escherichia coli. Here, we describe BMScan (Bacterial Meningitis Scan), a whole-genome analysis tool for the species identification of bacterial meningitis-causing and closely-related pathogens, an essential step for case management and disease surveillance. BMScan relies on a reference collection that contains genomes for 17 focal species to scan against to identify a given species. We established this reference collection by supplementing publically available genomes from RefSeq with genomes from the isolate collections of the Centers for Disease Control Bacterial Meningitis Laboratory and the Minnesota Department of Health Public Health Laboratory, and then filtered them down to a representative set of genomes which capture the diversity for each species. Using this reference collection, we evaluated two genomic comparison algorithms, Mash and Average Nucleotide Identity, for their ability to accurately and rapidly identify our focal species. Results We found that the results of Mash were strongly correlated with the results of ANI for species identification, while providing a significant reduction in run-time. This drastic difference in run-time enabled the rapid scanning of large reference genome collections, which, when combined with species-specific threshold values, facilitated the development of BMScan. Using a validation set of 15,503 genomes of our species of interest, BMScan accurately identified 99.97% of the species within 16 min 47 s. Conclusions Identification of the bacterial meningitis pathogenic species is a critical step for case confirmation and further strain characterization. BMScan employs species-specific thresholds for previously-validated, genome-wide similarity statistics compiled from a curated reference genome collection to rapidly and accurately identify the species of uncharacterized bacterial meningitis pathogens and closely related pathogens. BMScan will facilitate the transition in public health laboratories from traditional phenotypic detection methods to whole genome sequencing based methods for species identification