Inferring the heritability of bacterial traits in the era of machine learning

Quantification of heritability is a fundamental desideratum in genetics, which allows an assessment of the contribution of additive genetic variation to the variability of a trait of interest. The traditional computational approaches for assessing the heritability of a trait have been developed in the field of quantitative genetics. However, the rise of modern population genomics with large sample sizes has led to the development of several new machine learning based approaches to inferring heritability. In this paper, we systematically summarize recent advances in machine learning which can be used to infer heritability. We focus on an application of these methods to bacterial genomes, where heritability plays a key role in understanding phenotypes such as antibiotic resistance and virulence, which are particularly important due to the rising frequency of antimicrobial resistance. By designing a heritability model incorporating realistic patterns of genome-wide linkage disequilibrium for a frequently recombining bacterial pathogen, we test the performance of a wide spectrum of different inference methods, including also GCTA. In addition to the synthetic data benchmark, we present a comparison of the methods for antibiotic resistance traits for multiple bacterial pathogens. Insights from the benchmarking and real data analyses indicate a highly variable performance of the different methods and suggest that heritability inference would likely benefit from tailoring of the methods to the specific genetic architecture of the target organism

Quantifying macrodetritus fluxes from a small temperate estuary

Empirical measurements of estuary-to-coast material fluxes usually exclude the fraction of primary production that is exported as macrodetritus (marine plant litter), potentially leaving a gap in our understanding of the role of estuaries as outwelling systems. To address this gap, we sampled water and suspended material seasonally from the mouth of Pepe Inlet, Tairua Estuary, New Zealand. From samples collected hourly over 24 h, we calculated the lateral tidal fluxes (import, export, net flux) of macrodetritus, particulate and dissolved forms of nitrogen (N) and phosphorus (P). Annually, the inlet was a net exporter of N and P (5145 and 362 kg respectively). However, macrodetritus accounted for 87%). Nevertheless, seasonal pulses in the source and supply of macrodetritus may have consequences for the temporal scales over which this resource subsidy affects receiving ecosystems (e.g. intertidal sandflats). These mensurative investigations are useful to inform estuarine nutrient budgets that quantify the ecosystem services provided by temperate estuaries (e.g. contribution to fisheries food webs)

Robust analysis of prokaryotic pangenome gene gain and loss rates with Panstripe

Horizontal gene transfer (HGT) plays a critical role in the evolution and diversification of many microbial species. The resulting dynamics of gene gain and loss can have important implications for the development of antibiotic resistance and the design of vaccine and drug interventions. Methods for the analysis of gene presence/absence patterns typically do not account for errors introduced in the automated annotation and clustering of gene sequences. In particular, methods adapted from ecological studies, including the pangenome gene accumulation curve, can be misleading as they may reflect the underlying diversity in the temporal sampling of genomes rather than a difference in the dynamics of HGT. Here, we introduce Panstripe, a method based on generalized linear regression that is robust to population structure, sampling bias, and errors in the predicted presence/absence of genes. We show using simulations that Panstripe can effectively identify differences in the rate and number of genes involved in HGT events, and illustrate its capability by analyzing several diverse bacterial genome data sets representing major human pathogens.Peer reviewe

Clonal expansion within pneumococcal serotype 6C after use of seven-valent vaccine

Streptococcus pneumoniae causes invasive infections, primarily at the extremes of life. A seven-valent conjugate vaccine (PCV7) is used to protect against invasive pneumococcal disease in children. Within three years of PCV7 introduction, we observed a fourfold increase in serotype 6C carriage, predominantly due to a single clone. We determined the whole-genome sequences of nineteen S. pneumoniae serotype 6C isolates, from both carriage (n = 15) and disease (n = 4) states, to investigate the emergence of serotype 6C in our population, focusing on a single multi-locus sequence type (MLST) clonal complex 395 (CC395). A phylogenetic network was constructed to identify different lineages, followed by analysis of variability in gene sets and sequences. Serotype 6C isolates from this single geographical site fell into four broad phylogenetically distinct lineages. Variation was seen in the 6C capsular locus and in sequences of genes encoding surface proteins. The largest clonal complex was characterised by the presence of lantibiotic synthesis locus. In our population, the 6C capsular locus has been introduced into multiple lineages by independent capsular switching events. However, rapid clonal expansion has occurred within a single MLST clonal complex. Worryingly, plasticity exists within current and potential vaccine-associated loci, a consideration for future vaccine use, target selection and design

A high-throughput multiplexing and selection strategy to complete bacterial genomes

Background: Bacterial whole-genome sequencing based on short-read technologies often results in a draft assembly formed by contiguous sequences. The introduction of long-read sequencing technologies permits those contiguous sequences to be unambiguously bridged into complete genomes. However, the elevated costs associated with long-read sequencing frequently limit the number of bacterial isolates that can be long-read sequenced. Here we evaluated the recently released 96 barcoding kit from Oxford Nanopore Technologies (ONT) to generate complete genomes on a high-throughput basis. In addition, we propose an isolate selection strategy that optimizes a representative selection of isolates for long-read sequencing considering as input large-scale bacterial collections. Results: Despite an uneven distribution of long reads per barcode, near-complete chromosomal sequences (assembly contiguity = 0.89) were generated for 96 Escherichia coli isolates with associated short-read sequencing data. The assembly contiguity of the plasmid replicons was even higher (0.98), which indicated the suitability of the multiplexing strategy for studies focused on resolving plasmid sequences. We benchmarked hybrid and ONT-only assemblies and showed that the combination of ONT sequencing data with short-read sequencing data is still highly desirable (i) to perform an unbiased selection of isolates for long-read sequencing, (ii) to achieve an optimal genome accuracy and completeness, and (iii) to include small plasmids underrepresented in the ONT library. Conclusions: The proposed long-read isolate selection ensures the completion of bacterial genomes that span the genome diversity inherent in large collections of bacterial isolates. We show the potential of using this multiplexing approach to close bacterial genomes on a high-throughput basis.Peer reviewe

Pre-vaccine serotype composition within a lineage signposts its serotype replacement – a carriage study over 7 years following pneumococcal conjugate vaccine use in the UK

Serotype replacement has been reported in carriage and disease after pneumococcal conjugate vaccine (PCV) introductions in the UK and globally. We previously described concurrent expansion and decline of sequence types associated with serotype replacement over 5 years following PCV introductions in the UK. Here we use whole-genome sequencing to fully characterise the population structure of pneumococcal isolates collected over seven winters encompassing PCV7 and PCV13 introductions in the UK, investigating the importance of lineages in serotype replacement. We analysed 672 pneumococcal genomes from colonised children of 4 years old or less. The temporal prevalence of 20 lineages, defined by hierarchical Bayesian analysis of population structure (BAPS), was assessed in the context of serotype replacement. Multiple serotypes were detected in the primary winter of sampling within three vaccine-type (VT) lineages BAPS4, BAPS10 and BAPS11, in which serotype replacement were observed. In contrast, serotype replacement was not seen in the remaining three VT lineages (BAPS1, BAPS13 and BAPS14), that expressed a single serotype (6B, 6A and 3, respectively) in the primary winter. One lineage, BAPS1 serotype 6B was undetectable in the population towards the end of the study period. The dynamics of serotype replacement, in this UK population, was preceded by the presence or absence of multiple serotypes within VT lineages, in the pre-PCV population. This observation could help predict which non-vaccine types (NVTs) may be involved in replacement in future PCV introductions here and elsewhere. It could further indicate whether any antibiotic resistance associated with the lineages is likely to be affected by replacement.</p

Comparative Genomics of Carriage and Disease Isolates of Streptococcus pneumoniae Serotype 22F Reveals Lineage-Specific Divergence and Niche Adaptation

Streptococcus pneumoniae is a major cause of meningitis, sepsis, and pneumonia worldwide. Pneumococcal conjugate vaccines have been part of the United Kingdom’s childhood immunization program since 2006 and have significantly reduced the incidence of disease due to vaccine efficacy in reducing carriage in the population. Here we isolated two clones of 22F (an emerging serotype of clinical concern, multilocus sequence types 433 and 698) and conducted comparative genomic analysis on four isolates, paired by Sequence Type (ST) with one of each pair being derived from carriage and the other disease (sepsis). The most compelling observation was of nonsynonymous mutations in pgdA, encoding peptidoglycan N-acetylglucosamine deacetylase A, which was found in the carriage isolates of both ST433 and 698. Deacetylation of pneumococcal peptidoglycan is known to enable resistance to lysozyme upon invasion. Althought no other clear genotypic signatures related to disease or carriage could be determined, additional intriguing comparisons between the two STs were possible. These include the presence of an intact prophage, in addition to numerous additional phage insertions, within the carriage isolate of ST433. Contrasting gene repertoires related to virulence and colonization, including bacteriocins, lantibiotics, and toxin-–antitoxin systems, were also observed

Effect on nasopharyngeal pneumococcal carriage of replacing PCV7 with PCV13 in the Expanded Programme of Immunization in The Gambia.

INTRODUCTION: In 2011, two years after the introduction of 7-valent Pneumococcal conjugate vaccine (PCV7), the Gambian immunization programme replaced PVC7 with PCV13 (13-valent). Our objective was to assess the additional impact of PCV13 on prevalence of pneumococcal nasopharyngeal carriage. METHODS: We recruited healthy Gambian infants who had received three PCV doses. Nasopharyngeal swabs were collected from infants and their mothers during two cross-sectional surveys (CSS) conducted in infants vaccinated with PCV7 (CSS1) and vaccinated with PCV13 (CSS2). Pneumococci were isolated and serotyped following standardized methods. Whole genome sequencing was performed on non-typable pneumococcus isolated in CSS1 and CSS2. RESULTS: 339 and 350 infants and their mothers were recruited in CSS1 and CSS2, respectively. Overall prevalence of pneumococcal carriage was 85.4% in infants. Among infants, prevalence of vaccine type (VT) carriage was lower in CSS2 [9.4% versus 4.9% (p=0.025) for PCV7-VT; 33.3% versus 18.3% (p<0.001) for PCV13-VT and 23.9% versus 13.7% (p=0.001) for the 6 additional serotypes included in PCV13]. At CSS2, there was a decrease of serotypes 6A (from 15.3% to 5.7%, p<0.001) and 19F (from 5.6% to 1.7%, p=0.007), and an increase of non-typable pneumococci (0.3-6.0%, p<0.001), most of which (82.4%) were from typable serotype backgrounds that had lost the ability to express a capsule. Prevalence of overall and VT carriage in mothers was similar in CSS1 and CSS2. CONCLUSIONS: Replacing PCV7 for PCV13 rapidly decreased prevalence of VT carriage among vaccinated Gambian infants. An indirect effect in mothers was not observed yet. Vaccine-driven selection pressure may have been responsible for the increase of non-typable isolates

Strong pathogen competition in neonatal gut colonisation

Opportunistic bacterial pathogen species and their strains that colonise the human gut are generally understood to compete against both each other and the commensal species colonising this ecosystem. Currently we are lacking a population-wide quantification of strain-level colonisation dynamics and the relationship of colonisation potential to prevalence in disease, and how ecological factors might be modulating these. Here, using a combination of latest high-resolution metagenomics and strain-level genomic epidemiology methods we performed a characterisation of the competition and colonisation dynamics for a longitudinal cohort of neonatal gut microbiomes. We found strong inter- and intra-species competition dynamics in the gut colonisation process, but also a number of synergistic relationships among several species belonging to genus Klebsiella, which includes the prominent human pathogen Klebsiella pneumoniae. No evidence of preferential colonisation by hospital-adapted pathogen lineages in either vaginal or caesarean section birth groups was detected. Our analysis further enabled unbiased assessment of strain-level colonisation potential of extra-intestinal pathogenic Escherichia coli (ExPEC) in comparison with their propensity to cause bloodstream infections. Our study highlights the importance of systematic surveillance of bacterial gut pathogens, not only from disease but also from carriage state, to better inform therapies and preventive medicine in the future.Peer reviewe

Global distribution of invasive serotype 35D streptococcus pneumoniae isolates following introduction of 13-valent pneumococcal conjugate vaccine

A newly recognized pneumococcal serotype 35D, which differs from the 35B polysaccharide in structure and serology by not binding to factor serum 35a, was recently reported. The genetic basis for this distinctive serology is due to the presence of an inactivating mutation in wciG, which encodes an O-acetyltransferase responsible for O-acetylation of a galactofuranose. Here, we assessed the genomic data of a worldwide pneumococcal collection to identify serotype 35D isolates and understand their geographical distribution, genetic background and invasiveness potential. Of 21,980 pneumococcal isolates, 444 were originally typed as serotype 35B by PneumoCaT. Analysis of wciGrevealed 23 isolates from carriage (n=4) and disease (n=19) with partial or complete loss-of-funtion mutations, including mutations resulting in pre-mature stop codons (n=22) and an in-frame mutation (n=1). These were selected for further analysis. The putative 35D isolates were geographically widespread and 65.2% (15/23) of them was recovered after PCV13 introduction. Compared with serotype 35B, putative serotype 35D isolates have higher invasive disease potentials based on odds ratio (OR) (11.58; 95% CI, 1.42-94.19 vs 0.61; 95% CI, 0.40-0.92) and a higher prevalence of macrolide resistance mediated by mefA (26.1% vs 7.6%, p=0.009). Using Quellung, 50% (10/20) of viable isolates were serotype 35D, 25% (5/20) serotype 35B, and 25% (5/20) a mixture of 35B/35D. The discrepancy between phenotype and genotype requires further investigation. These findings illustrated a global distribution of an invasive serotype 35D among young children post-PCV13 introduction and underlined the invasive potential conferred by the loss of O-acetylation in the pneumococcal capsule