PrimerEvalPy: a tool for in-silico evaluation of primers for targeting the microbiome.

BACKGROUND: The selection of primer pairs in sequencing-based research can greatly influence the results, highlighting the need for a tool capable of analysing their performance in-silico prior to the sequencing process. We therefore propose PrimerEvalPy, a Python-based package designed to test the performance of any primer or primer pair against any sequencing database. The package calculates a coverage metric and returns the amplicon sequences found, along with information such as their average start and end positions. It also allows the analysis of coverage for different taxonomic levels. RESULTS: As a case study, PrimerEvalPy was used to test the most commonly used primers in the literature against two oral 16S rRNA gene databases containing bacteria and archaea. The results showed that the most commonly used primer pairs in the oral cavity did not match those with the highest coverage. The best performing primer pairs were found for the detection of oral bacteria and archaea. CONCLUSIONS: This demonstrates the importance of a coverage analysis tool such as PrimerEvalPy to find the best primer pairs for specific niches. The software is available under the MIT licence at https://gitlab.citius.usc.es/lara.vazquez/PrimerEvalPy .

Critical review of 16S rRNA gene sequencing workflow in microbiome studies: From primer selection to advanced data analysis.

The multi-batch reanalysis approach of jointly reevaluating gene/genome sequences from different works has gained particular relevance in the literature in recent years. The large amount of 16S ribosomal ribonucleic acid (rRNA) gene sequence data stored in public repositories and information in taxonomic databases of the same gene far exceeds that related to complete genomes. This review is intended to guide researchers new to studying microbiota, particularly the oral microbiota, using 16S rRNA gene sequencing and those who want to expand and update their knowledge to optimise their decision-making and improve their research results. First, we describe the advantages and disadvantages of using the 16S rRNA gene as a phylogenetic marker and the latest findings on the impact of primer pair selection on diversity and taxonomic assignment outcomes in oral microbiome studies. Strategies for primer selection based on these results are introduced. Second, we identified the key factors to consider in selecting the sequencing technology and platform. The process and particularities of the main steps for processing 16S rRNA gene-derived data are described in detail to enable researchers to choose the most appropriate bioinformatics pipeline and analysis methods based on the available evidence. We then produce an overview of the different types of advanced analyses, both the most widely used in the literature and the most recent approaches. Several indices, metrics and software for studying microbial communities are included, highlighting their advantages and disadvantages. Considering the principles of clinical metagenomics, we conclude that future research should focus on rigorous analytical approaches, such as developing predictive models to identify microbiome-based biomarkers to classify health and disease states. Finally, we address the batch effect concept and the microbiome-specific methods for accounting for or correcting them.

In silico evaluation and selection of the best 16S rRNA gene primers for use in next-generation sequencing to detect oral bacteria and archaea.

BACKGROUND: Sequencing has been widely used to study the composition of the oral microbiome present in various health conditions. The extent of the coverage of the 16S rRNA gene primers employed for this purpose has not, however, been evaluated in silico using oral-specific databases. This paper analyses these primers using two databases containing 16S rRNA sequences from bacteria and archaea found in the human mouth and describes some of the best primers for each domain. RESULTS: A total of 369 distinct individual primers were identified from sequencing studies of the oral microbiome and other ecosystems. These were evaluated against a database reported in the literature of 16S rRNA sequences obtained from oral bacteria, which was modified by our group, and a self-created oral archaea database. Both databases contained the genomic variants detected for each included species. Primers were evaluated at the variant and species levels, and those with a species coverage (SC) ≥75.00% were selected for the pair analyses. All possible combinations of the forward and reverse primers were identified, with the resulting 4638 primer pairs also evaluated using the two databases. The best bacteria-specific pairs targeted the 3-4, 4-7, and 3-7 16S rRNA gene regions, with SC levels of 98.83-97.14%; meanwhile, the optimum archaea-specific primer pairs amplified regions 5-6, 3-6, and 3-6, with SC estimates of 95.88%. Finally, the best pairs for detecting both domains targeted regions 4-5, 3-5, and 5-9, and produced SC values of 95.71-94.54% and 99.48-96.91% for bacteria and archaea, respectively. CONCLUSIONS: Given the three amplicon length categories (100-300, 301-600, and >600 base pairs), the primer pairs with the best coverage values for detecting oral bacteria were as follows: KP_F048-OP_R043 (region 3-4; primer pair position for Escherichia coli J01859.1: 342-529), KP_F051-OP_R030 (4-7; 514-1079), and KP_F048-OP_R030 (3-7; 342-1079). For detecting oral archaea, these were as follows: OP_F066-KP_R013 (5-6; 784-undefined), KP_F020-KP_R013 (3-6; 518-undefined), and OP_F114-KP_R013 (3-6; 340-undefined). Lastly, for detecting both domains jointly they were KP_F020-KP_R032 (4-5; 518-801), OP_F114-KP_R031 (3-5; 340-801), and OP_F066-OP_R121 (5-9; 784-1405). The primer pairs with the best coverage identified herein are not among those described most widely in the oral microbiome literature. Video Abstract.

Impact of 16S rRNA Gene Redundancy and Primer Pair Selection on the Quantification and Classification of Oral Microbiota in Next-Generation Sequencing.

This study aimed to evaluate the number of 16S rRNA genes in the complete genomes of the bacterial and archaeal species inhabiting the human mouth and to assess how the use of different primer pairs would affect the detection and classification of redundant amplicons and matching amplicons (MAs) from different taxa. A total of 518 oral-bacterial and 191 oral-archaeal complete genomes were downloaded from the NCBI database, and their complete 16S rRNA genes were extracted. The numbers of genes and variants per genome were calculated. Next, 39 primer pairs were used to search for matches in the genomes and obtain amplicons. For each primer, we calculated the number of gene amplicons, variants, genomes, and species detected and the percentage of coverage at the species level with no MAs (SC-NMA). The results showed that 94.09% of oral bacteria and 52.59% of oral archaea had more than one intragenomic 16S rRNA gene. From 1.29% to 46.70% of bacterial species and from 4.65% to 38.89% of archaea detected by the primers had MAs. The best primers were the following (SC-NMA; region; position for Escherichia coli [GenBank version no. J01859.1]): KP_F048-OP_R030 for bacteria (93.55%; V3 to V7; 342 to 1079), KP_F018-KP_R063 for archaea (89.63%; V3 to V9; undefined to 1506), and OP_F114-OP_R121 for both domains (92.52%; V3 to V9; 340 to 1405). In addition to 16S rRNA gene redundancy, the presence of MAs must be controlled to ensure an accurate interpretation of microbial diversity data. The SC-NMA is a more useful parameter than the conventional coverage percentage for selecting the best primer pairs. The pairs used the most in the oral microbiome literature were not among the best performers. IMPORTANCE Hundreds of publications have studied the oral microbiome through 16S rRNA gene sequencing. However, none have assessed the number of 16S rRNA genes in the genomes of oral microbes, or how the use of primer pairs targeting different regions affects the detection of MAs from different taxa. Here, we found that almost all oral bacteria and more than half of oral archaea have more than one intragenomic 16S rRNA gene. The performance of the primer pairs in not detecting MAs increases as the length of the amplicon augments. As none of those most employed in the oral literature were among the best performers, we selected a series of primers to detect bacteria and/or archaea based on their percentage of species detected without MAs. The intragenomic 16S rRNA gene redundancy and the presence of MAs between distinct taxa need to be considered to ensure an accurate interpretation of microbial diversity data.

Update on the Role of Cytokines as Oral Biomarkers in the Diagnosis of Periodontitis.

Periodontitis is one of the world's most common chronic human diseases and has a significant impact on oral health. Recent evidence has revealed a link between periodontitis and certain severe systemic conditions. Moreover, periodontal patients remain so for life, even following successful therapy, requiring ongoing supportive care to prevent the disease's recurrence. The first challenge in treating the condition is ensuring a timely and accurate diagnosis since the loss of periodontal bone and soft tissue is progressive and largely irreversible. Although current clinical and radiographic parameters are the best available for identifying and monitoring the disease, the scientific community has a particular interest in finding quantifiable biomarkers in oral fluids that can improve early detection rates of periodontitis and evaluations of its severity. It is widely accepted that periodontitis is associated with polymicrobial dysbiosis and a chronic inflammatory immune response in the host. This response causes the generation of mediators like cytokines. Higher concentrations of cytokines are involved in inflammation and disease progression, acting as a network of biological redundancy. Most of the cytokines investigated concerning the periodontitis pathogenesis are proinflammatory. Of all of them, interleukin (IL) 1beta has been studied the most, followed by tumor necrosis factor (TNF) alpha and IL6. In contrast, only a few papers have evaluated antiinflammatory cytokines, with the most researched being IL4 and IL10. Several systemic reviews have concluded that the specific cytokines present in patients with periodontitis have a distinctive profile, which may indicate their possible discriminatory potential. In this chapter, the focus is on analyzing studies that investigate the accuracy of diagnoses of periodontitis based on the cytokines present in gingival crevicular fluid and saliva. The findings of our research group are also described.

In-Silico Detection of Oral Prokaryotic Species With Highly Similar 16S rRNA Sequence Segments Using Different Primer Pairs.

Although clustering by operational taxonomic units (OTUs) is widely used in the oral microbial literature, no research has specifically evaluated the extent of the limitations of this sequence clustering-based method in the oral microbiome. Consequently, our objectives were to: 1) evaluate in-silico the coverage of a set of previously selected primer pairs to detect oral species having 16S rRNA sequence segments with ≥97% similarity; 2) describe oral species with highly similar sequence segments and determine whether they belong to distinct genera or other higher taxonomic ranks. Thirty-nine primer pairs were employed to obtain the in-silico amplicons from the complete genomes of 186 bacterial and 135 archaeal species. Each fasta file for the same primer pair was inserted as subject and query in BLASTN for obtaining the similarity percentage between amplicons belonging to different oral species. Amplicons with 100% alignment coverage of the query sequences and with an amplicon similarity value ≥97% (ASI97) were selected. For each primer, the species coverage with no ASI97 (SC-NASI97) was calculated. Based on the SC-NASI97 parameter, the best primer pairs were OP_F053-KP_R020 for bacteria (region V1-V3; primer pair position for Escherichia coli J01859.1: 9-356); KP_F018-KP_R002 for archaea (V4; undefined-532); and OP_F114-KP_R031 for both (V3-V5; 340-801). Around 80% of the oral-bacteria and oral-archaea species analyzed had an ASI97 with at least one other species. These very similar species play different roles in the oral microbiota and belong to bacterial genera such as Campylobacter, Rothia, Streptococcus and Tannerella, and archaeal genera such as Halovivax, Methanosarcina and Methanosalsum. Moreover, ~20% and ~30% of these two-by-two similarity relationships were established between species from different bacterial and archaeal genera, respectively. Even taxa from distinct families, orders, and classes could be grouped in the same possible OTU. Consequently, regardless of the primer pair used, sequence clustering with a 97% similarity provides an inaccurate description of oral-bacterial and oral-archaeal species, which can greatly affect microbial diversity parameters. As a result, OTU clustering conditions the credibility of associations between some oral species and certain health and disease conditions. This significantly limits the comparability of the microbial diversity findings reported in oral microbiome literature.

Diagnostic accuracy of IL1ß in saliva: The development of predictive models for estimating the probability of the occurrence of periodontitis in non-smokers and smokers.

AIM: To obtain salivary interleukin (IL) 1ß-based models to predict the probability of the occurrence of periodontitis, differentiating by smoking habit. MATERIALS/METHODS: A total of 141 participants were recruited, 62 periodontally healthy controls and 79 subjects affected by periodontitis. Fifty of the diseased patients were given non-surgical periodontal treatment and showed significant clinical improvement in 2 months. IL1ß was measured in the salivary samples using the Luminex instrument. Binary logistic regression models were obtained to differentiate untreated periodontitis from periodontal health (first modelling) and untreated periodontitis from treated periodontitis (second modelling), distinguishing between non-smokers and smokers. The area under the curve (AUC) and classification measures were calculated. RESULTS: In the first modelling, IL1ß presented AUC values of 0.830 for non-smokers and 0.689 for smokers (accuracy = 77.6% and 70.7%, respectively). In the second, the predictive models revealed AUC values of 0.671 for non-smokers and 0.708 for smokers (accuracy = 70.0% and 75.0%, respectively). CONCLUSION: Salivary IL1ß has an excellent diagnostic capability when it comes to distinguishing systemically healthy patients with untreated periodontitis from those who are periodontally healthy, although this discriminatory potential is reduced in smokers. The diagnostic capacity of salivary IL1ß remains acceptable for differentiating between untreated and treated periodontitis.

Accuracy of single molecular biomarkers in saliva for the diagnosis of periodontitis: A systematic review and meta-analysis.

AIM: To analyse, using a meta-analytical approach, the diagnostic accuracy of single molecular biomarkers in saliva for the detection of periodontitis in systemically healthy subjects. MATERIALS AND METHODS: Articles on molecular biomarkers in saliva providing a binary contingency table (or sensitivity and specificity values and group sample sizes) in individuals with clinically diagnosed periodontitis were considered eligible. Searches for candidate articles were conducted in six electronic databases. The methodological quality was assessed through the tool Quality Assessment of Diagnostic Studies. Meta-analyses were performed using the Hierarchical Summary Receiver Operating Characteristic model. RESULTS: Meta-analysis was possible for 5 of the 32 biomarkers studied. The highest values of sensitivity for the diagnosis of periodontitis were obtained for IL1beta (78.7%), followed by MMP8 (72.5%), IL6 and haemoglobin (72.0% for both molecules); the lowest sensitivity value was for MMP9 (70.3%). In terms of specificity estimates, MMP9 had the best result (81.5%), followed by IL1beta (78.0%) and haemoglobin (75.2%); MMP8 had the lowest specificity (70.5%). CONCLUSIONS: MMP8, MMP9, IL1beta, IL6 and Hb were salivary biomarkers with good capability to detect periodontitis in systemically healthy subjects. MMP8 and IL1beta are the most researched biomarkers in the field, both showing clinically fair effectiveness for the diagnosis of periodontitis.

Accuracy of single molecular biomarkers in gingival crevicular fluid for the diagnosis of periodontitis: A systematic review and meta-analysis.

AIM: To analyse, by means of a meta-analytical approach, the diagnostic accuracy of molecular biomarkers in gingival crevicular fluid (GCF) for the detection of periodontitis in systemically healthy subjects. MATERIAL AND METHODS: Studies on GCF molecular biomarkers providing a binary classification table (or sensitivity and specificity values and group sample sizes) in individuals with clinically diagnosed periodontitis were considered eligible. The search was performed using six electronic databases. The methodological quality of studies was assessed through the tool Quality Assessment of Diagnostic Studies. Meta-analyses were performed using the Hierarchical Summary Receiver Operating Characteristic, which adjusts classification data using random effects logistic regression. RESULTS: The included papers identified 36 potential biomarkers for the detection of periodontitis and for four of them meta-analyses were performed. The median sensitivity and specificity were for MMP8, 76.7% and 92.0%; for elastase, 74.6% and 81.1%; for cathepsin, 72.8% and 67.3%, respectively. The worst estimates of sensitivity and specificity were for trypsin (71.3% and 66.1%, respectively). CONCLUSIONS: MMP8 showed good sensitivity and excellent specificity, which resulted in this biomarker being clinically the most useful or effective for the diagnosis of periodontitis in systemically healthy subjects, regardless of smoking condition.

Quantification by qPCR of Pathobionts in Chronic Periodontitis: Development of Predictive Models of Disease Severity at Site-Specific Level.

Currently, there is little evidence available on the development of predictive models for the diagnosis or prognosis of chronic periodontitis based on the qPCR quantification of subgingival pathobionts. Our objectives were to: (1) analyze and internally validate pathobiont-based models that could be used to distinguish different periodontal conditions at site-specific level within the same patient with chronic periodontitis; (2) develop nomograms derived from predictive models. Subgingival plaque samples were obtained from control and periodontal sites (probing pocket depth and clinical attachment loss <4 mm and >4 mm, respectively) from 40 patients with moderate-severe generalized chronic periodontitis. The samples were analyzed by qPCR using TaqMan probes and specific primers to determine the concentrations of Actinobacillus actinomycetemcomitans (Aa), Fusobacterium nucleatum (Fn), Parvimonas micra (Pm), Porphyromonas gingivalis (Pg), Prevotella intermedia (Pi), Tannerella forsythia (Tf), and Treponema denticola (Td). The pathobiont-based models were obtained using multivariate binary logistic regression. The best models were selected according to specified criteria. The discrimination was assessed using receiver operating characteristic curves and numerous classification measures were thus obtained. The nomograms were built based on the best predictive models. Eight bacterial cluster-based models showed an area under the curve (AUC) ≥0.760 and a sensitivity and specificity ≥75.0%. The PiTfFn cluster showed an AUC of 0.773 (sensitivity and specificity = 75.0%). When Pm and AaPm were incorporated in the TdPiTfFn cluster, we detected the two best predictive models with an AUC of 0.788 and 0.789, respectively (sensitivity and specificity = 77.5%). The TdPiTfAa cluster had an AUC of 0.785 (sensitivity and specificity = 75.0%). When Pm was incorporated in this cluster, a new predictive model appeared with better AUC and specificity values (0.787 and 80.0%, respectively). Distinct clusters formed by species with different etiopathogenic role (belonging to different Socransky's complexes) had a good predictive accuracy for distinguishing a site with periodontal destruction in a periodontal patient. The predictive clusters with the lowest number of bacteria were PiTfFn and TdPiTfAa, while TdPiTfAaFnPm had the highest number. In all the developed nomograms, high concentrations of these clusters were associated with an increased probability of having a periodontal site in a patient with chronic periodontitis.