Strain-level resolution and pneumococcal carriage dynamics by single-molecule real-time (SMRT) sequencing of the plyNCR marker: a longitudinal study in Swiss infants

BACKGROUND: Pneumococcal carriage has often been studied from a serotype perspective; however, little is known about the strain-specific carriage and inter-strain interactions. Here, we examined the strain-level carriage and co-colonization dynamics of Streptococcus pneumoniae in a Swiss birth cohort by PacBio single-molecule real-time (SMRT) sequencing of the plyNCR marker. METHODS: A total of 872 nasal swab (NS) samples were included from 47 healthy infants during the first year of life. Pneumococcal carriage was determined based on the quantitative real-time polymerase chain reaction (qPCR) targeting the lytA gene. The plyNCR marker was amplified from 214 samples having lytA-based carriage for pneumococcal strain resolution. Amplicons were sequenced using SMRT technology, and sequences were analyzed with the DADA2 pipeline. In addition, pneumococcal serotypes were determined using conventional, multiplex PCR (cPCR). RESULTS: PCR-based plyNCR amplification demonstrated a 94.2% sensitivity and 100% specificity for Streptococcus pneumoniae if compared to lytA qPCR. The overall carriage prevalence was 63.8%, and pneumococcal co-colonization (≥ 2 plyNCR amplicon sequence variants (ASVs)) was detected in 38/213 (17.8%) sequenced samples with the relative proportion of the least abundant strain(s) ranging from 1.1 to 48.8% (median, 17.2%; IQR, 5.8–33.4%). The median age to first acquisition was 147 days, and having ≥ 2 siblings increased the risk of acquisition. CONCLUSION: The plyNCR amplicon sequencing is species-specific and enables pneumococcal strain resolution. We therefore recommend its application for longitudinal strain-level carriage studies of Streptococcus pneumoniae. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1186/s40168-022-01344-6.