Design and validation of Dolosigranulum pigrum specific PCR primers using the bacterial core genome

Scientific Reports volume 13, Article number: 6110 (2023) Cite this article

651 Accesses

2 Altmetric

Metrics details

Dolosigranulum pigrum—a lactic acid bacterium that is increasingly recognized as an important member of the nasal microbiome. Currently, there are limited rapid and low-cost options for confirming D. pigrum isolates and detecting D. pigrum in clinical specimens. Here we describe the design and validation of a novel PCR assay targeting D. pigrum that is both sensitive and specific. We designed a PCR assay targeting murJ, a single-copy core species gene identified through the analysis of 21 D. pigrum whole genome sequences. The assay achieved 100% sensitivity and 100% specificity against D. pigrum and diverse bacterial isolates and an overall 91.1% sensitivity and 100% specificity using nasal swabs, detecting D. pigrum at a threshold of 1.0 × 104 D. pigrum 16S rRNA gene copies per swab. This assay adds a reliable and rapid D. pigrum detection tool to the microbiome researcher toolkit investigating the role of generalist and specialist bacteria in the nasal environment.

Dolosigranulum pigrum is a gram-positive, non-spore forming bacterium from the family Carnobacteriaceae1 commonly found in the human nasal cavity2,3. First described in 1993 as small, white colonies that displayed beta-hemolysis1, D. pigrum remains poorly understood and the only species of Dolosigranulum known to date. Epidemiologically, D. pigrum has been associated with the healthy state of the nasal microbiome3,4. Specifically, upper airway colonization by D. pigrum is negatively associated with Staphylococcus aureus carriage5,6,7,8,9. More recently, D. pigrum was found to be in higher abundance in the nasopharynx of patients with asymptomatic SARS-CoV-2 infections than patients with more severe symptoms10.

Rapid and cost-effective methods for the identification of D. pigrum are needed to facilitate future clinical and in vitro studies. Standard biochemical methods are expensive and time-consuming, as are sequencing-based methods. Matrix-assisted laser desorption/ionization-time of flight (MALDI-TOF) analysis is cost-effective, but cannot be used to detect D. pigrum directly from clinical samples. Using a core-genome based approach, we designed and validated a PCR-based assay that can be used to confirm D. pigrum isolates and detect the presence of D. pigrum directly from clinical samples.

We first analyzed the genetic diversity of 21 D. pigrum whole genome sequences available (n = 7 from NCBI and n = 14 from in-house D. pigrum genomes, Table S1). We extracted 87,993 SNPs from non-recombined regions of the core genome and examined the genetic diversity based on maximum likelihood phylogeny (Figure S1). This showed multiple distinct D. pigrum lineages, which indicates both the non-clonal nature of D. pigrum and the robustness of the genome collection. We then generated and analyzed the D. pigrum pan-genome to identify 1291 core and 357 accessory genes. For potential assay targets, we focused on the 1291 core genes.

Assay target genes discovery was a multi-step process (Fig. 1). We removed ribosomal genes (n = 71) and genes with homologs in other genera (n = 345). Manual filtering of randomly-chosen assay target genes from the 843 single-copy core species genes was performed, requiring that the assay target gene: (a) must be present in all 21 D. pigrum genomes, (b) must have less than 70% similarity identity and coverage against sequences from non-Dolosigranulum taxa by BLAST, (c) contain forward and reverse primer sequences meeting Primer3 design criteria and that have less than 50% similarity identity and cover against sequences from non-Dolosigranulum taxa by BLAST.

Core genome-based approach for assay design. Schematic representation of the approach taken to mine the pan-genome for assay targets. Each succeeding step in the pangenome analysis workflow illustrates how genes were filtered to finally retain a unique core genome for the organism of interest.

The first single-copy core species genes (SCSG) that met our selection criteria as a target gene candidate with conserved regions for primer design was murJ (Pfam ID: PF01943), a gene with a length of 1665 bp encoding a peptidoglycan lipid II flippase protein. The average uncorrected distance between the isolates for the murJ alignment was 35.84 bp (SD = 13.67 bp) (Fig. 2a). After iterations of primer design and in silico analysis, we identified a pair of forward and reverse PCR primers (Table 1, Supplementary Table S4a–d) targeting the murJ gene that produces a 223 bp PCR product. On average the amplicon varied by 2.14 bp (SD = 1.69 bp) between the isolates (Fig. 2b, Supplementary Table S3a,b, Supplementary File S1).

Dolosigranulum pigrum murJ phylogeny and sequence alignment. (a) Neighbor joining tree constructed using full length murJ gene sequences from 21 D. pigrum isolates using Jalview 2.11 37 and ordered by branch lengths, highlighting that murJ is part of the conserved core genome but is also phylogenetically informative; (b) multiple sequence alignment of murJ amplicon region, where the forward primer is located at 1234–1255 bp and the reverse primer is located at 1436–1457 bp.

The murJ assay was highly sensitive and specific in laboratory analysis of DNA from bacterial isolates and from human nasal swabs. We first evaluated the assay using well-characterized D. pigrum isolates (N = 12) and against five common nasal bacterial species namely Moraxella catarrhalis, Staphylococcus aureus, Staphylococcus epidermidis, Corynebacterium pseudodiphtheriticum, Corynebacterium propinquum, Corynebacterium accolens, which showed 100% sensitivity and specificity (Figure S3).

We further evaluated the assay using DNA extracted from human nasal swabs (n = 110) characterized using 16S rRNA V3–V4 gene-based sequencing, including 54 samples that were positive for D. pigrum and 56 samples that were negative for D. pigrum. This showed that the murJ assay was not able to detect D. pigrum in samples (n = 9) with fewer than ten D. pigrum 16S rRNA gene copies per uL of swab eluent, or 1.0 × 104 D. pigrum 16S rRNA gene copies per swab. However, among the 45 D. pigrum-positive samples with more than 1.0 × 104 D. pigrum 16S rRNA gene copies per swab, the murJ PCR assay was able to detect D. pigrum in 41 (91%) samples (Table 2, Figs. S4, S5). There were no false positives in the 56 D. pigrum-negative samples.

By identifying potential assay targets using the D. pigrum core genome, we designed a novel PCR assay that is both sensitive and specific for D. pigrum. In contrast to other commonly used methods for species confirmation, such as biochemical testing, DNA sequencing, or MALDI-TOF, PCR-based assays are rapid and cost-effective and do not require expensive equipment. This method provides a simpler option for D. pigrum detection and avoids the restriction digestion and analysis challenges of T-RFLP11 that has been used previously for detecting microbial communities in anterior nares12. We demonstrated the utility of the core genome mining techniques to develop species confirmation assays. The resultant murJ assay was able to identify D. pigrum and diverse bacterial isolates with a 100% sensitivity and specificity. Our assay was also highly sensitive and specific for detecting D. pigrum in clinical samples.

Dolosigranulum pigrum is gaining interest as a member of the upper respiratory tract microbial community that is potentially beneficial for the host5,6,8,13,14,15,16,17,18,19. Efforts are being made to better understand its metabolic models and defense mechanisms20. There is a critical need to screen samples to detect the presence of D. pigrum or to verify the identity of the organism isolated through culture-based methods. Our single step gel-based PCR method for the species verification of D. pigrum in clinical samples as well as pure isolates provides a useful tool for epidemiological and clinical studies.

We curated a local D. pigrum genome database by downloading publicly available genomes from NCBI RefSeq and adding in-house sequenced and assembled D. pigrum genomes (Table S1). DNA from the inhouse D. pigrum isolates was extracted using a DNeasy Blood and Tissue kit (Qiagen) or MagNA Pure LC DNA Isolation Kit (Roche) and libraries were generated with a Nextera XT DNA Library kit (Illumina) according to manufacturer’s instructions for paired-end sequencing on an Illumina NextSeq 500 (Illumina, Inc., San Diego, CA) with a read length of 150 bp. We assembled Illumina short read sequences from inhouse D. pigrum isolates into contigs using the SPADES assembler (v.3.5)21. Quality of the assembly was assessed using metrics generated by QUAST (v.2.3)22 and all genomes were annotated with Prokka (v. 1.13)23. To maximize assay sensitivity for D. pigrum detection we focused on the core genome. The GFF files from the Prokka annotation step were used as input for the pan-genome analysis with Roary (v.3.12.0)24 [blastp v.2.9.0 identity = 90%, gene presence in isolates to be core = 99%]. We generated a maximum likelihood tree from core genome SNPs to assess relatedness of the D. pigrum isolates using previously described methods25,26. Briefly, Illumina short reads from inhouse D. pigrum isolates were mapped to the chromosome of the published D. pigrum reference genome (strain 83VPs-KB5; GenBank accession no. CP041626.1) using the NASP pipeline that uses BWA-MEM (v.0.7.12)27 to align and GATK (v.3.5)28 to call SNPs. Publicly available genomes downloaded from NCBI RefSeq were aligned to the reference using MUMMER and SNPs were identified. The resultant SNP matrix was processed with Gubbins29 to remove recombinant regions. A Phylogenetic tree was constructed from the core SNPs in PhyML with Smart Model selection (v.3.0)30. The maximum likelihood phylogeny was visualized alongside the pangenome using PHANDANGO31 (Figure S2). Uniprot IDs of the core genes wherever available, were extracted from the GFF files using an inhouse script and were used to retrieve Gene Ontology terms from UniProt database32 (Table S2). The GO terms were analyzed and summarized using GAOTools33.

The core genome was filtered and only SCSG were retained. An in-silico search for homology against non-D. pigrum species was performed using blastn v.2.9.034 using a local copy of the NT database (updated: 2019-03-31). Gene targets with 70% similarity to non-D. pigrum species were removed. A final set of homologous single-copy core genes was used as the candidate pool for targets to design D. pigrum specific assay.

We used Primer335 with default settings to identify candidate forward and reverse primers which were first compared to the D. pigrum gene alignment file then checked for similarity against other nasal bacteria, including Staphylococcus aureus, Staphylococcus epidermidis, Corynebacterium spp., Cutibacterium spp., Moraxella spp., Escherichia coli, Klebsiella spp., Citrobacter spp., Proteus spp., and Alloiococcus spp. Primers were excluded if 5 or more matching bases were found at the 3′-end of the primer.

To assess the sensitivity of our primers, we tested the murJ assay against 12 D. pigrum isolates. These isolates had been previously verified to be D. pigrum by MALDI-TOF and their genomes were sequenced using Illumina HiSeq system (Illumina, San Diego, CA). Furthermore, we screened murJ primers against 110 clinical samples characterized by 16S rRNA gene-based sequencing as described previously36. A non-D. pigrum control collection that included Moraxella catarrhalis, Staphylococcus aureus, Staphylococcus epidermidis, Corynebacterium pseudodiphtheriticum, Corynebacterium propinquum, Corynebacterium accolens species was used to evaluate specificity of our primers.

Ethical approval for this study was granted by the George Washington University Institutional Review Board and The Office of Human Research. The study and its protocols were implemented according to the approved guidelines outlined in the Declaration of Helsinki. Informed consent was obtained from all participants prior to enrollment in the study.

The first study included 16 healthy community-dwelling adults in Washington, DC (IRB#: NCR191444) were included. At enrollment, nasal specimens were self-collected by participants under staff guidance using Puritan HydraFlock swabs (Puritan Medical Products, Guilford, ME) with staff instructions. Samples were placed immediately into Amies transport media and stored at 4 °C until processing. Samples were processed within 4 h then transferred in 100 μL aliquots into labeled 2 mL cryovials and stored at − 80 °C. The second study included 94 healthy community-dwelling adults in Copenhagen, Denmark (IRB # 041631), which were collected by study personnel and collected into DNA/RNA shield (Zymo R1100-250) and stored at − 80 °C until processing.

DNA from human nasal swabs were extracted using MagMax DNA Ultra 2.0 Kit with enzyme and chemical lysis as previously described 5. DNA from bacterial isolates were extracted through heat soak (D. pigrum, S. aureus, and S. epidermidis) or using the DNeasy Blood & Tissue Kit (Qiagen, Valencia, CA) (C. propinquum and C. pseudodiphtheriticum) according to manufacturer instructions.

Each murJ PCR was performed in a 20 μL reaction volume containing 1 μL of template DNA added to 19 μL of PCR reaction mix containing 0.4 μM of forward (5′-CAACAGCGTCCAGCAATCTA-3′) and reverse (5′-ATCGCTGTAATCCCGATGAG-3′) primer, 1× Phusion High-Fidelity PCR Master Mix (ThermoFisher), and molecular-grade water. Amplification was performed on a C1000 Touch Thermocycler (Bio-Rad, Hercules, CA) using the following conditions: 98 °C for 30 s for denaturing, 54 °C for 30 s for annealing, and 72 °C for 1 min for extension × 35 cycles. Amplified DNA was run on a 2% agarose E-gel (ThermoFisher) to assess amplification of D. pigrum DNA. Gels were imaged using a ChemiDoc-It2 (Analytik Jena US, Upland, CA). Presence of a visible band at the 223 bp size indicated successful amplification.

Raw Reads generated from the whole genome sequencing performed for this study were deposited at NCBI SRA (Accession ID: PRJNA770953).

Aguirre, M., Morrison, D., Cookson, B. D., Gay, F. W. & Collins, M. D. Phenotypic and phylogenetic characterization of some Gemella-like organisms from human infections: Description of Dolosigranulum pigrum gen. nov., sp. nov.. J. Appl. Bacteriol. 75, 608–612 (1993).

Article CAS PubMed Google Scholar

Toivonen, L. et al. Early nasal microbiota and acute respiratory infections during the first years of life. Thorax 74, 592–599 (2019).

Article PubMed Google Scholar

Kaspar, U. et al. The culturome of the human nose habitats reveals individual bacterial fingerprint patterns. Environ. Microbiol. 18, 2130–2142 (2016).

Article CAS PubMed Google Scholar

Gan, W. et al. The difference in nasal bacterial microbiome diversity between chronic rhinosinusitis patients with polyps and a control population. Int. Forum Allergy Rhinol. 9, 582–592 (2019).

Article PubMed Google Scholar

Liu, C. M. et al. Staphylococcus aureus and the ecology of the nasal microbiome. Sci. Adv. 1, e1400216 (2015).

Article ADS PubMed PubMed Central Google Scholar

Accorsi, E. K. et al. Determinants of Staphylococcus aureus carriage in the developing infant nasal microbiome. Genome Biol. 21, 301 (2020).

Article CAS PubMed PubMed Central Google Scholar

Biesbroek, G. et al. The impact of breastfeeding on nasopharyngeal microbial communities in infants. Am. J. Respir. Crit. Care Med. 190, 298–308 (2014).

Article PubMed Google Scholar

Biesbroek, G. et al. Early respiratory microbiota composition determines bacterial succession patterns and respiratory health in children. Am. J. Respir. Crit. Care Med. 190, 1283–1292 (2014).

Article PubMed Google Scholar

Teo, S. M. et al. The infant nasopharyngeal microbiome impacts severity of lower respiratory infection and risk of asthma development. Cell Host Microbe 17, 704–715 (2015).

Article CAS PubMed PubMed Central Google Scholar

Hurst, J. H. et al. Age-related changes in the nasopharyngeal microbiome are associated with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection and symptoms among children, adolescents, and young adults. Clin. Infect. Dis. 75, e928–e937 (2022).

Article PubMed PubMed Central Google Scholar

Prakash, O., Pandey, P. K., Kulkarni, G. J., Mahale, K. N. & Shouche, Y. S. Technicalities and glitches of terminal restriction fragment length polymorphism (T-RFLP). Indian J. Microbiol. 54, 255–261 (2014).

Article CAS PubMed PubMed Central Google Scholar

Camarinha-Silva, A., Wos-Oxley, M. L., Jauregui, R., Becker, K. & Pieper, D. H. Validating T-RFLP as a sensitive and high-throughput approach to assess bacterial diversity patterns in human anterior nares. FEMS Microbiol. Ecol. 79, 98–108 (2012).

Article CAS PubMed Google Scholar

RayaTonetti, F. et al. The respiratory commensal bacterium Dolosigranulum pigrum 040417 improves the innate immune response to Streptococcus pneumoniae. Microorganisms 9, 25 (2021).

Google Scholar

De Boeck, I. et al. Anterior nares diversity and pathobionts represent sinus microbiome in chronic rhinosinusitis. mSphere 4, 25 (2019).

Article Google Scholar

Coleman, A. et al. Upper respiratory microbiota in relation to ear and nose health among Australian aboriginal and torres strait islander children. J. Pediatr. Infect. Dis. Soc. 10, 468–476 (2021).

Article Google Scholar

Hurst, J. H. et al. Age-related changes in the upper respiratory microbiome are associated with SARS-CoV-2 susceptibility and illness severity. medRxiv 20, 25 (2021).

Google Scholar

Laufer, A. S. et al. Microbial communities of the upper respiratory tract and otitis media in children. MBio 2, e00245-e1210 (2011).

Article PubMed PubMed Central Google Scholar

Ortiz Moyano, R. et al. The ability of respiratory commensal bacteria to beneficially modulate the lung innate immune response is a strain dependent characteristic. Microorganisms 8, 25 (2020).

Article Google Scholar

Bosch, A. et al. Maturation of the infant respiratory microbiota, environmental drivers, and health consequences. A prospective cohort study. Am. J. Respir. Crit. Care Med. 196, 1582–1590 (2017).

Article PubMed Google Scholar

Flores Ramos, S. et al. Genomic stability and genetic defense systems in Dolosigranulum pigrum, a candidate beneficial bacterium from the human microbiome. MSystems 20, e0042521 (2021).

Article Google Scholar

Nurk, S. et al. Assembling single-cell genomes and mini-metagenomes from chimeric MDA products. J. Comput. Biol. 20, 714–737 (2013).

Article MathSciNet CAS PubMed PubMed Central Google Scholar

Gurevich, A., Saveliev, V., Vyahhi, N. & Tesler, G. QUAST: Quality assessment tool for genome assemblies. Bioinformatics 29, 1072–1075 (2013).

Article CAS PubMed PubMed Central Google Scholar

Seemann, T. Prokka: Rapid prokaryotic genome annotation. Bioinformatics 30, 2068–2069 (2014).

Article CAS PubMed Google Scholar

Page, A. J. et al. Roary: Rapid large-scale prokaryote pan genome analysis. Bioinformatics 31, 3691–3693 (2015).

Article CAS PubMed PubMed Central Google Scholar

Price, L. B. et al. Staphylococcus aureus CC398: Host adaptation and emergence of methicillin resistance in livestock. MBio 3, 15 (2012).

Article Google Scholar

Reid, C. J., McKinnon, J. & Djordjevic, S. P. Clonal ST131-H22 Escherichia coli strains from a healthy pig and a human urinary tract infection carry highly similar resistance and virulence plasmids. Microb. Genom 5, 25 (2019).

Google Scholar

Li, H. & Durbin, R. Fast and accurate short read alignment with Burrows–Wheeler transform. Bioinformatics 25, 1754–1760 (2009).

Article CAS PubMed PubMed Central Google Scholar

McKenna, A. et al. The Genome Analysis Toolkit: A MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 20, 1297–1303 (2010).

Article CAS PubMed PubMed Central Google Scholar

Croucher, N. J. et al. Rapid phylogenetic analysis of large samples of recombinant bacterial whole genome sequences using Gubbins. Nucleic Acids Res. 43, e15 (2015).

Article PubMed Google Scholar

Guindon, S. et al. New algorithms and methods to estimate maximum-likelihood phylogenies: Assessing the performance of PhyML 3.0. Syst. Biol. 59, 307–321 (2010).

Article CAS PubMed Google Scholar

Hadfield, J. et al. Phandango: An interactive viewer for bacterial population genomics. Bioinformatics 34, 292–293 (2018).

Article CAS PubMed Google Scholar

UniProt, C. UniProt: The universal protein knowledgebase in 2021. Nucleic Acids Res. 49, D480–D489 (2021).

Article Google Scholar

Klopfenstein, D. V. et al. GOATOOLS: A Python library for Gene Ontology analyses. Sci. Rep. 8, 10872 (2018).

Article ADS CAS PubMed PubMed Central Google Scholar

Camacho, C. et al. BLAST+: Architecture and applications. BMC Bioinform. 10, 421 (2009).

Article Google Scholar

Untergasser, A. et al. Primer3—new capabilities and interfaces. Nucleic Acids Res. 40, e115 (2012).

Article CAS PubMed PubMed Central Google Scholar

Liu, C. M. et al. Male circumcision significantly reduces prevalence and load of genital anaerobic bacteria. MBio 4, e00076 (2013).

Article PubMed PubMed Central Google Scholar

Waterhouse, A. M., Procter, J. B., Martin, D. M., Clamp, M. & Barton, G. J. Jalview Version 2—a multiple sequence alignment editor and analysis workbench. Bioinformatics 25, 1189–1191 (2009).

Article CAS PubMed PubMed Central Google Scholar

Download references

This work was supported by research Grants from National Institutes of Health (R01AI125562) to LBP and by a training Grant from the Milken Institute School of Public Health to AP. The contents of this publication are solely the responsibility of the authors and do not necessarily represent the official views of the funding agencies.

These authors contributed equally: Maliha Aziz and Amber Palmer.

Antibiotic Resistance Action Center, Department of Environmental and Occupational Health, Milken Institute School of Public Health, George Washington University, 800 22nd Street NW, Washington, DC, 20052, USA

Maliha Aziz, Amber Palmer, Juan E. Salazar, Tony Pham, Kelsey Roach, Lance B. Price & Cindy M. Liu

Department of Bacteria, Parasites, and Fungi, Statens Serum Institut, Artillerivej 5, 2300, Copenhagen, Denmark

Søren Iversen, Sharmin Baig, Marc Stegger & Paal Skytt Andersen

Friedrich Loeffler-Institute of Medical Microbiology, University Medicine Greifswald, Greifswald, Germany

Karsten Becker

Institute of Medical Microbiology, University Hospital Münster, Münster, Germany

Ursula Kaspar

You can also search for this author in PubMed Google Scholar

M.A. and A.P. contributed equally to the preparation of this manuscript. C.L. conceived the study. M.A. performed data mining and bioinformatics analysis. M.A. and A.P. analyzed results. K.B., U.K., P.A., and M.S. provided study and laboratory materials. S.I. and S.B. generated the D. pigrum genome reference. A.P., K.R. T.P. and J.S. performed assay validation experiments. All authors contributed to and reviewed the manuscript.

Correspondence to Cindy M. Liu.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Aziz, M., Palmer, A., Iversen, S. et al. Design and validation of Dolosigranulum pigrum specific PCR primers using the bacterial core genome. Sci Rep 13, 6110 (2023). https://doi.org/10.1038/s41598-023-32709-y

Download citation

Received: 15 December 2022

Accepted: 31 March 2023

Published: 14 April 2023

DOI: https://doi.org/10.1038/s41598-023-32709-y

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.