ACE2 expression in saliva of patients with COVID-19 and its association with Candida albicans and Aggregatibacter actinomycetemcomitans

Endang W Bachtiar; Boy M Bachtiar; Ardiana Kusumaningrum; Hari Sunarto; Yuniarti Soeroso; Benso Sulijaya; Efa Apriyanti; Citra Fragrantia Theodorea; Irandi Putra Pratomo; Yudhistira .; Defi Efendi; Fathilah Abdul Razak

doi:10.12688/f1000research.111965.2

Home Browse ACE2 expression in saliva of patients with COVID-19 and its association...

ALL Metrics

Views

Downloads

Get PDF

Get XML

Export

▬

✚

Research Article

Revised

ACE2 expression in saliva of patients with COVID-19 and its association with Candida albicans and Aggregatibacter actinomycetemcomitans

[version 2; peer review: 2 approved]

Endang W Bachtiar ¹, Boy M Bachtiar¹, Ardiana Kusumaningrum², [...] Hari Sunarto³, Yuniarti Soeroso⁴, Benso Sulijaya⁴, Efa Apriyanti⁵, Citra Fragrantia Theodorea¹, Irandi Putra Pratomo⁶, Yudhistira .⁷, Defi Efendi⁸, Fathilah Abdul Razak⁹

Endang W Bachtiar ¹, Boy M Bachtiar¹, [...] Ardiana Kusumaningrum², Hari Sunarto³, Yuniarti Soeroso⁴, Benso Sulijaya⁴, Efa Apriyanti⁵, Citra Fragrantia Theodorea¹, Irandi Putra Pratomo⁶, Yudhistira .⁷, Defi Efendi⁸, Fathilah Abdul Razak⁹

PUBLISHED 12 Sep 2022

Author details Author details

¹ Department of Oral Biology, Faculty of Dentistry, Universitas Indonesia - Oral Science Research Center, Jakarta, DKI, 10430, Indonesia
² Department of Microbiology, Faculty of Medicine, Universitas Indonesia - Clinical Microbiology Medicine Staff Group, Universitas Indonesia Hospital, Depok, West Java, 16424, Indonesia
³ Dental Center, Universitas Indonesia Hospital, Depok, West Java, 16424, Indonesia
⁴ Department of Periodontology, Faculty of Dentistry, Universitas Indonesia, Jakarta, DKI, 10430, Indonesia
⁵ Department of Pediatric Nursing, Faculty of Nursing, Universitas Indonesia and Paediatric Intensive Care Unit, Universitas Indonesia Hospital, Depok, West Java, 16424, Indonesia
⁶ Department of Pulmonology and Respiratory Medicine, Faculty of Medicine, Universitas Indonesia, Jakarta, DKI, 10430, Indonesia
⁷ Clinical Pathology Medicine Staff Group, Universitas Indonesia Hospital, Depok, West Java, 16424, Indonesia
⁸ Department of Pediatric Nursing, Faculty of Nursing, Universitas Indonesia, and Neonatal Intensive Care Unit, Universitas Indonesia Hospital, Depok, West Java, 16424, Indonesia
⁹ Department of Oral & Craniofacial Sciences, Faculty of Dentistry, University of Malaya, Kuala Lumpur, 50603, Malaysia

Endang W Bachtiar
Roles: Conceptualization, Funding Acquisition, Methodology, Supervision, Writing – Review & Editing

Boy M Bachtiar
Roles: Conceptualization, Formal Analysis, Supervision, Writing – Original Draft Preparation, Writing – Review & Editing

Ardiana Kusumaningrum
Roles: Data Curation, Formal Analysis, Writing – Review & Editing

Hari Sunarto
Roles: Resources, Supervision, Writing – Review & Editing

Yuniarti Soeroso
Roles: Validation, Visualization, Writing – Review & Editing

Benso Sulijaya
Roles: Investigation, Resources, Writing – Review & Editing

Efa Apriyanti
Roles: Data Curation, Resources, Writing – Review & Editing

Citra Fragrantia Theodorea
Roles: Project Administration, Visualization, Writing – Review & Editing

Irandi Putra Pratomo
Roles: Conceptualization, Methodology, Writing – Review & Editing

Yudhistira .
Roles: Data Curation, Resources, Software, Writing – Review & Editing

Defi Efendi
Roles: Investigation, Resources, Writing – Review & Editing

Fathilah Abdul Razak
Roles: Conceptualization, Methodology, Writing – Review & Editing

OPEN PEER REVIEW

REVIEWER STATUS

This article is included in the Pathogens gateway.

This article is included in the Emerging Diseases and Outbreaks gateway.

This article is included in the Coronavirus (COVID-19) collection.

Abstract

Background: A relationship between oral microbiota and susceptibility to SARS-CoV-2 infection has been extensively studied. However, the relationship between oral commensal flora and expression of angiotensin-converting enzyme 2 (ACE2) remains to be established. In this observational study, we collected saliva from patients with COVID-19 and evaluated the relationship between ACE2 expression and Candida albicans as well as with selected gram-negative bacteria (Aggregatibacter actinomycetemcomitans, Fusobacterium nucleatum, and Veillonella parvula). We investigated how this may be directly or indirectly involved in oral dysbiosis in patients with COVID-19.
Methods: We included 23 hospitalized patients admitted to Universitas Indonesia Hospital with PCR-confirmed COVID-19, with six healthy participants serving as controls. Saliva and tongue surface swabs were collected from patients with diabetes (DG) and without diabetes (NDG) and subject controls. Using quantitative PCR (qPCR) we assessed the mRNA expression of ACE2, the abundance of C. albicans, and the transcription levels of its biofilm-associated genes, agglutinin-like protein 3 (ALS3), hyphal wall protein 1 (HWP1), and yeast-form wall protein 1 (YWP1). We also counted the relative proportion of the three selected gram-negative oral bacteria in saliva. All analyses were performed to determine the relationship between ACE2 expression and C. albicans and gram-negative bacteria.
Results: ACE2 mRNA expression was significantly higher in tongue swab samples than in saliva. However, no significant difference was observed between the patient groups. Conversely, DG patients had a significantly higher abundance of C. albicans in saliva compared to NDG patients and control group patients. The correlation and sensitivity/specificity relationship between ACE2 expression and C. albicans or the selected oral bacteria were also observed.
Conclusions: The data show that ACE2 expression can be detected in saliva of patients with COVID-19 and its association with C. albicans and gram-negative oral bacteria might contribute toward developing an oral dysbiosis based predictor for prognosis of COVID-19 severity.

Keywords

COVID-19, ACE2, Candida albicans, Aggregatibacter actinomycetemcomitans, Fusobacterium nucleatum, Veillonella parvula

Corresponding author: Endang W Bachtiar

Competing interests: No competing interests were disclosed.

Grant information: This study was supported by a grant from Universitas Indonesia (NKB-638/UN2.RST/HKP.05.00/2021).
The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Copyright: © 2022 Bachtiar EW et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite: Bachtiar EW, Bachtiar BM, Kusumaningrum A et al. ACE2 expression in saliva of patients with COVID-19 and its association with Candida albicans and Aggregatibacter actinomycetemcomitans [version 2; peer review: 2 approved]. F1000Research 2022, 11:557 (https://doi.org/10.12688/f1000research.111965.2) First published: 23 May 2022, 11:557 (https://doi.org/10.12688/f1000research.111965.1) Latest published: 12 Sep 2022, 11:557 (https://doi.org/10.12688/f1000research.111965.2)

Revised Amendments from Version 1

We modify the manuscript according to the reviewer's suggestion in Introduction, Figure 1 legend, discussion section and the limitation of the project.

See the authors' detailed response to the review by Juni Handajani
See the authors' detailed response to the review by Jukun Song

Introduction

The oral microbiota may be involved in the pathogenesis of SARS-CoV-2 infection, the causative agent of COVID-19. In addition to other oral bacteria, Candida spp., particularly C. albicans, is a keystone commensal in the human oral cavity¹ that may be involved in dysbiotic events. Indeed, C. albicans is commonly reported for its relationships with known constituents of the oral biofilm in individuals with and without oral disease.²^–⁴ As dysbiosis of the oral microbiome has been associated with inflammatory conditions in the oral habitat of COVID-19 patients with comorbidity, as reported by Bachtiar et al. (preprint),⁵ we assumed that fungal-bacterial interactions might also favor the establishment of SARS-CoV-2 infection. Therefore, our objective was to investigate the level of C. albicans, its pathogenicity, and to evaluate its antagonistic relationship with Aggregatibacter actinomycetemcomitans, which is a keystone pathogen associated with periodontitis in adolescents⁶ in the saliva of COVID-19 patients with and without diabetes. We included Fusobacterium nucleatum and Veillonella parvula as their relationship has been previously reported.⁷

Methods

Study design, patients, and specimens

The study was conducted at the Universitas Indonesia Hospital (RSUI), Depok, Indonesia. The eligible patients were recruited consecutively (up to 23), from August 2021 to September 2021. According to medical records, the patients had mild to moderate symptoms with clinically and laboratory confirmed COVID-19 infections at RSUI. Six subjects who visited the RSUI periodontal clinic served as a control. The average age of the participants was 45.1 ± 15.37 years old, and 10 patients had diabetes.

According to guidance provided by the Ethics Committee, written and oral information was given, after which written informed consent was obtained from all participants before enrolment in this study. The study protocol was approved by the ethics committee of Universitas Indonesia Hospital (protocol number: 2021/04/052). The protocol conformed to the criteria of the Declaration of Helsinki and the good clinical practical guidelines of the International Council on Harmonization, and this study was carried out in accordance with the guidelines provided by the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) statement.

Unstimulated saliva (2 ml) was collected by spitting into a sterile Falcon tube. Tongue samples were taken by swabbing the middle third of the tongue dorsum with a sterile cotton swab for a few minutes.⁸ The obtained samples were then put into a microcentrifuge tube. All collected samples (saliva and tongue swab) were delivered promptly to the laboratory for further processing.

Quantification of C. albicans by quantitative PCR

Fungal genome extraction was performed using GENEzol™ reagents, (phenol, guanidine isothiocyanate solution) (Geneaid Biotech Ltd, New Taipei City, Taiwan), accordance with the protocol provided by the company. The concentration and quality of the obtained DNA were determined using Qubit assay reagents (Thermo Fisher Scientific, Waltham, MA, United States). To amplify DNA in saliva samples, we used quantitative PCR (qPCR) with specific primers for C. albicans as follows: Forward: 5′-CACGACGGAGTTTCACAAGA-3′ and Reverse: CGATGGAAGTTTGAGGCAAT-3′.⁹ Further, the fungus amount was calculated by plotting the cycle threshold (Ct) value against the log of a standard curve shown in Figure 1A. The standard curve was constructed using a 10-fold serial dilution of DNA extracted from C. albicans (ATCC 10231).⁸ The amplicon melting curves was set at 95°C for 15 seconds, 60°C for 60 secods, and 95°C for 15 seconds.

Figure 1. Standard curves and melt curves. Standard curve of C. albicans ATCC 10231 (A) and representative melt curve graph for C. albicans count in saliva (B) and tongue surface (C) all participants (DG, NDG, and control/ n = 29).

The standard curve was linear as indicated by a linear correlation (r²) of 0.9974 between the cycle threshold (Ct) value and template concentration, and a slope of -3.705.

The PCR reactions were performed at a 10-minute initial denaturation at 95°C, followed by 40 cycles of denaturation at 95°C for 15-seconds annealing at 60°C for 60 second and elongation at 95°C for 15 seconds. As shown in Figure 1B and C, the qPCR products were visualized as a melting curve, that was set at 95°C for 15 seconds, 60°C for 60 seconds, and 95°C for 15 seconds.

Proportion of A. actinomycetemcomitans, F. nucleatum, and V. parvula in saliva microbiota

Bacterial DNA was extracted from the saliva samples with a similar procedure described above. The levels of A. actinomycetemcomitans, F. nucleatum, and V. parvula colonization and the total amount of bacteria in saliva were determined using qPCR with the specific primers,⁵ except for A. actinomycetemcomitans, for which we used oligonucleotides as reported elsewhere,¹⁰ as follows: GTGGGGAGCAAACAGGATTAG (forward) and CCTAAGGCACAAACCCATCTC (reverse).

For both C. albicans and bacterial abundance, the PCR cycling process was performed in a total volume of 10 μl (comprising 5 μl of SYBR1 Selected Master Mix (Thermo Fisher Scientific, Waltham, MA, USA), DNA template (2 μl), and primer pair solution (1 μl, 300 nM/reaction). The abundance of each bacterium was determined by using the 2^-ΔΔCt method.¹¹ ΔCt was the difference between the Ct value using the primers for each bacterium and the Ct value obtained by using the primers for total bacteria in saliva. ΔΔCt was the difference between the ΔCt of the patient and control subjects, where the value of 2^-ΔΔCt shows the changes in bacterial proportion in the sample of patients relative to those of the control subjects.

Expression of mRNA ALS3, HWP1, and YWP1 by qPCR

For extracting total RNA, we used GENEzolTM reagent (Geneaid; Biotech Ltd, New Taipei City, Taiwan), followed by a reverse transcription kit (High-Capacity cDNA Reverse Transcription Kit, Applied Biosystems^TM). We followed all instructions provided by the kits. The resulting cDNA was amplified by qPCR with specific primers, as follows: YWP1; F: 5′-GCTACTGCTACTGGTGCTA-3′, R: 5′-AACGGTGGTTTCTTGAC-3′, HWP1; F: 5′-GCTCCTGCTCCTGAAATGAC-3′, R: 5′-CTGGAGCAATTGGTGAGGTT-3′, and ALS1; F: 5′-CAACTTGGGTTATTGAAACAAAAACA-3′, R: 5′-AGAAACAGAAACCCAAGAACAACC-3′.⁹

Quantitative PCR analysis was performed in triplicate on an ABI StepOnePlus Real-Time PCR System with SYBR Green PCR Master Mix (Applied Biosystems). The qPCR cycling conditions consisted of a 10-minute initial denaturation at 95°C followed by 40 PCR cycles of 15 seconds at 95°C and 1 minute at 60°C. The formula of fold change 2-ΔΔCt was used to calculate the relative mRNA expression, which was compared with that of the housekeeping gene, ACT1 with primers: F: 5′-TTTCATCTTCTGTATCAGAGGAACTTATTT-3′, R: 5′-ATGGGATGAATCATCAAACAAGAG -3′.¹² The formula of fold change 2^-∆∆Ct was used to calculate the relative mRNA expression of genes (ALS3, HWP1, and YWP1), which was normalized to that of the housekeeping gene ACT1.¹³ All values obtained from the tested patient groups were standardized and compared to the values obtained from the control subjects.

Data analysis

In this study, we compared the amount of C. albicans, relative abundance of bacteria, and mRNA transcription levels of the targeted genes in two groups: patients with COVID-19 with diabetes (DG) and patients without diabetes (NDG). Statistical analyses were conducted using GraphPad Prism 9.0 (GraphPad Software, San Diego, CA, USA) (RRID:SCR_002798) (An open-access alternative is the R Stats Package 4.3.0). One way ANOVA and an unpaired Student’s t-test were used to determine the p-values between and within groups, respectively. Data are presented as the mean ± standard error (SE), and p < 0.05 was considered significant. Spearman’s correlation coefficient (r) with two-tailed p-values was used to measure the degree of association between two variables tested. The line of best fit (95% confidence interval) was shown by using linear regression. The receiver operating characteristic (ROC) method was also used to determine the sensitivity and specificity of the relationship between C. albicans and A. actinomycetemcomitans as predictors of oral dysbiosis in patients with COVID-19.

Results

The mRNA level of ACE2 and the abundance of C. albicans

As shown in Figure 2A. ACE2 mRNA expression was found in all saliva samples collected from either subject tested (DG, NDG, and control). This finding was confirmed by comparing ACE2 transcription levels on tongue surface sample (TS), where ACE2 transcription level is highly expressed in this niche.¹⁴ We observed, that mean transcription level of ACE2 detected in saliva was lower than on the TS (p < 0.05). When comparing the two group, we found that in both saliva and TS, the transcription levels of ACE2 was higher in NDG than in DG, but the difference was not significant (p > 0.05). We further determined that C. albicans was also present in the saliva and TS samples collected from all subjects. In general, the count of C. albicans in saliva samples was higher than that of TS samples (p < 0.05). Additionally, the amount (log DNA copies) of C. albicans in the saliva of DG was significantly higher than that in NDG (p < 0.05). In contrast, the different number of C. albicans on tongue surface found in either group was not significant (p > 0.05) (Figure 2B).

Figure 2. Evaluation of ACE2 expression, the number of C. albicans/the proportion of selected gram-negative oral bacteria [Aggregatibacter actinomycetemcomitans (Aa), Fusobacterium nucleatum (Fn), and Veillonella parvula (Vp)], and mRNA expression of selected biofilm-associated genes (ALS3, HWP1, and YWP1), in different patient groups (DG and NDG) of COVID-19.

In general, the transcription level of ACE2 mRNA is higher in tongue surface (Ts) than those in saliva (Sa). Between groups, the different level of ACE2 expression is not significant (A). The number of C. albicans are significantly higher (p < 0.05) in saliva than in tongue surface samples, and C. albicans abundance is significantly higher in DG than in NDG (B). The proportion of Vp is found to be the highest compared to the other two species (Aa and Fn). Between groups, the abundance of Aa and Vp is significantly higher in DG than in NDG (C). Fold changes in gene expression detected in saliva (D) and tongue surface (E) were each compared between and within patient groups (DG and NDG). All quantitative polymerase chain reaction values were normalized according to the expression of housekeeping gene ACT1. All data are expressed as mean ± SE. *p < 0.05, **p < 0.001.

Relative abundance of A. actinomycetemcomitans, F. nucleatum, and V. parvula

The qPCR results showed that the proportion of each bacterium in all subjects tested was lower in DG than in NDG. We found that the proportion of A. actinomycetemcomitans and V. parvula in NDG was >20% higher than that in DG (<10%). Subsequently, in either group tested, the abundance of F. nucleatum was found to be the lowest (<5%) compared to the proportion of the other two species (Figure 2C).

Transcription levels of C. albicans biofilm-associated genes

The qPCR results showed that the transcription levels of ALS3 and HWP1 in saliva were significantly lower than those in YWP1 (p < 0.05). We found that the transcription of both hypha-associated genes (ALS3 and HWP1) was significantly higher in NDG than in DG (p < 0.05), whereas no difference was found in the expression of YWP1 mRNA (Figure 2D). Furthermore, we analyzed the relative expression levels of each gene on the tongue surface (TS). As expected, in both groups tested, the hypha-related genes (ALS3 and HWP1) showed upregulation at a similar level. Conversely, a higher level of YWP1 mRNA expression was significantly detected on the TS of NDG subjects than in DG subjects (Figure 2E).

Moreover, as shown in Figures 3A and B, a strong negative linear correlation was observed between the abundance of C. albicans and the relative proportion of A. actinomycetemcomitans in DG (r = -0.79, p = 0.02), whereas in NDG, a low negative, non-significant correlation was observed between the two oral microorganisms (r = 0.08, p = 0.79). We noted, that in DG there was a positive but not significant correlation between the proportion of C. albicans and F. nucleatum (r = 0.05; p = 0.9) as well as with V. parvula (r = 0.24; p = 0.29) Conversely, a negative non-significant correlation was observed between C albicans and F. nucleatum/V. parvula. The correlation coefficients were r = -0.5, p = 0.08, and r = -0.41, p = 0.15, respectively (Figure 3C–F).

Figure 3. Relationship between the number (colony forming units/mL) of C. albicans (Ca) and the three gram-negative oral bacteria A. actinomycetemcomitans (Aa), F. nucleatum (Fn), and V. parvula (Vp).

In either patient group (DG or NDG), a strong negative correlation is consistently seen between Ca and Aa (A and B), while the other two bacteria show a weak positive in DG (C and E) and negative correlations in NDG, respectively (D and F). Spearman correlation coefficient (r²) and exact p-values are given.

As SARS-CoV-2 has been consistently detected in the saliva of infected patients,¹⁵ we further examined the relationship between ACE2 expression and the relative abundance of C. albicans and A. actinomycetemcomitans in the saliva of patients with COVID-19. As shown in Figure 4A–D, a strong and significant positive correlation between ACE2 mRNA transcription and the abundance of C. albicans/proportion of A. actinomycetemcomitans was observed in DG (r = 0.81, p = 0. 01 and r = 0.75, p = 0.03, respectively). In NDG, the correlation was positive, but not statistically significant (r = 0.53, p = 0.06 and r = 0.06, p = 0.83, respectively). Based on these results, we evaluated the accuracy of the combination of ACE2 and C. albicans/A. actinomycetemcomitans relationship analyses. We revealed that the area under the curve (AUC) of the ACE2/C. albicans association, in DG was 1 (95% CI: 1 to 1, p < 0.0002; Figure 5A), and in NDG, the AUC was 0.76 (95% CI: 0.54 to 0.99, p < 0.019; Figure 5B). For the relationship between ACE2 and A. actinomycetemcomitans, in DG the AUC was 0.80 (95% CI: 0.56 to 1, p < 0.02; Figure 5C). In NDG, the AUC was 0.82 (95% CI: 0.66 to 0.98, p < 0.005; Figure 5D).

Figure 4. Scatter diagram illustrating the correlation between mRNA expression of ACE2 and C. albicans/A. actinomycetemcomitans in patients with COVID-19 with (DG) and without (NDG) diabetes.

These observations indicate that in DG, the correlation between ACE2 mRNA expression and the abundance of C. albicans/A. actinomycetemcomitans is strongly positive (A and C), while in NDG the correlation is weakly positive (B and D). Spearman correlation coefficient (r2) and exact p-values are given.

Figure 5. Receiver operating characteristic curve (ROC) showing the plot and the best cut-off point of the relationship between the mRNA expression of ACE2 and C. albicans or A. actinomycetemcomitans in DG (A and C) and NDG (B and D), respectively.

Discussion

This study found that ACE2 expression was detected at a lower level in saliva than on the tongue surface, indicating that saliva from patients with COVID-19 might harbor epithelial cells containing SARS-CoV-2. However, the epithelium of the tongue is likely the primary target of SARS-COV-2 in the oral cavity, which is consistent with previous findings.¹⁶^,¹⁷ Therefore, although the oral environment is not the main target of SARS-CoV-2, the oral cavity could constitute both a portal of entry and a reservoir for the virus. This finding also indicated that ACE2 may be crucial for the progression and prognosis of COVID-19. Hence, exploring ACE2 expression under different physiological conditions may help predict the susceptibility of SARS-CoV-2 in different cohorts, such as COVID-19 patients with and without comorbid diabetes. The comorbidity observed in our subjects appeared to accelerate the expression of ACE2 mRNA in oral niches. Since the susceptibility to SARS-CoV-2-driven infection correlates with ACE2 expression,¹⁸ it is possible that in our subjects, either the tongue or epithelia-containing saliva have been exposed to SARS-CoV-2 infection. We assumed that this result could be linked to the clinical status of our COVID-19 patients, which, according to the patient’s medical history, only developed mild to moderate illness.¹⁹

In this study, we aimed to compare the interactions observed between C. albicans and gram-negative oral bacteria. As an important constituent of oral commensal flora, C. albicans shows a diverse inter-kingdom relationship under certain conditions, ranging from synergistic to antagonistic.²⁰^–²² The results of this study showed that the number of C. albicans in the saliva of both patient groups was increased compared to the fungal count on the tongue surface, but a significant increase was only observed in DG. This finding indicates that although saliva provides antimicrobial activity against the opportunistic oral fungal pathogen,²² this critical function is less effective in patients with COVID-19 with diabetes. Therefore, it is possible that the protective functions normally observed in saliva, such as cleansing, lubrication, and antibacterial activity, had changed in our patients with COVID-19.

Additionally, studies have shown that most patients with COVID-19 have one or more systemic (e.g., use of broad-spectrum antibiotics, use of corticosteroids, immunosuppression) or local (e.g., use of dental prostheses, reduced salivary flow due to use of medication) risk factors that favor Candida proliferation.²³^–²⁵

Additionally, most participants had poor oral hygiene (not shown), and we did not measure the local risk factors. However, the high count of salivary C. albicans observed may indicate that low salivary flow rate influences salivary colonization by this fungus.²⁶

According to the medical records, all patients had received a non-invasive supplemental device. However, we did not have any data regarding how long the patients had received the device at the time the samples were collected.

Previous studies have shown that Candida airway colonization is associated with prolonged use of mechanical ventilation and length of hospital stay.²⁶^,²⁷ According to the medical records, all patients had received a noninvasive supplemental device. However, we did not have any data regarding how long the patients had received the device at the time the samples were collected. Additionally, although the local risk factors had not been measured, the high count of salivary C. albicans found in this study may suggest that low salivary flow rate facilitates salivary colonization by this fungus,²⁸ as reported previously in elderly populations.⁸ Thus, we assumed that the physiological effects of COVID-19 on salivary gland secretion²⁹^,³⁰ might affect the salivary flow rate of patients with COVID-19.

Additionally, we analyzed the pathogenicity of C. albicans in saliva and on the tongue surface by comparing the expression of selected biofilm-associated genes (ALS3, HWP1, and YWP1), since these genes alters the morphology of the fungus, from the yeast to hyphal form.³¹

Our data showed that in both patient groups (DG and NDG), the transcription levels of ALS3 and HWP1 were downregulated in saliva but conversely upregulated in tongue sample swabs. This finding suggests that the tongue is an ideal biotic surface in oral habitats for C. albicans attachment and growth. This study also revealed no difference in the pathogenicity of C. albicans (evaluated by qPCR) between DG and NDG patients with COVID-19. Hence, the results showed that the tongue surface is a better oral sample for detecting C. albicans pathogenicity in patients with COVID-19, irrespective of the diabetic condition.

ALS3 (adhesion-related gene) and HWP1 (hypha-specific gene)³² are produced predominantly during biofilm formation, while YWP1 is a yeast-associated gene involved in the anti-adhesive activity of C. albicans.³¹ Therefore, the higher expression of ALS3 and HWP1 mRNAs on the tongue surface indicate that this is an important reservoir of C. albicans colonization. From here the fungus may dislodge into the saliva, where the expression of YWP1 was found to be higher.

An additional phenomenon revealed in this study was the relative proportion of gram-negative bacteria in patients with COVID-19. Our data indicated that the selected oral bacteria were found in the saliva of all participants recruited in the current investigation (DG, NDG, and control). This finding suggests that the bacteria exist as normal microflora in the oral cavity, and they may have been involved in disease processes observed in our patients with COVID-19.⁵ The current study highlighted that in the presence of C. albicans, the most abundant gram-negative bacteria in both diseased groups (DG and NDG) were V. parvula, followed by A. actinomycetemcomitans, while F. nucleatum was the least abundant. Moreover, all species can be detected in periodontitis patients with diabetes,³³ but only F. nucleatum and V. parvula have been reported to be positively associated with COVID-19-associated events.³⁴^,³⁵ Therefore, in this study, it was deemed pertinent to assess the relationship between the proportions of these three gram-negative bacteria and the abundance of C. albicans in the COVID-19 related oral environment. Analysis of saliva samples showed that in NDG (patients without diabetes), the proportion of all selected gram-negative bacteria had a significant negative correlation with the increasing load of C. albicans DNA. Interestingly, in DG (patients with diabetes), only the proportion of A. actinomyces was consistently and significantly negatively correlated with a higher count of C. albicans. This finding supports our previous work regarding the ability of A. actinomycetemcomitans to reduce biofilm formation by C. albicans when grown in mono³⁶ or dual species with Streptococcus mutans.⁶ The current findings provide additional information in vivo, whereby the salivary component in patients with COVID-19 modulates the relationship pattern shown by C. albicans when interacting with A. actinomycetemcomitans. In the present study, a negative association between the fungus and periodontal pathogen was observed by counting both microflora in saliva samples collected from patients with COVID-19. The relationship pattern may be relevant, since it not only translates closer to the real inter-kingdom relationships in vivo (oral cavity), but also demonstrates the relationships between C. albicans and other selected gram-negative bacteria in the salivary environment. From the results here, the negative association between C. albicans and A. actinomycetemcomitans observed in salivary microbiota could be multifactorial and not be simply caused by the effect of SARS-CoV-2 in the oral cavity. However, this may explain why the presence of SARS-CoV-2 in the oral cavity favors the emergence and persistence of dysbiosis in another oral niche, including the periodontal microenvironment.⁵ Our observations provide additional information, in which there is a strong negative correlation between the “key stone commensal” oral microflora, C. albicans¹ and periodontopathogen, A. actinomycetemcomitans.³⁷ Interestingly, we found that both C. albicans and A. actinomycetemcomitans showed a positive correlation with ACE2 expression, and a strong correlation was observed in patients with COVID-19 accompanied by diabetes. Therefore, it was deemed relevant to assess these relationships by determining the sensitivity and specificity of the association. Using ROC curve analyses, we found that the antagonistic relationship between C. albicans and A. actinomycetemcomitans and its respective correlation with ACE2 expression had an obvious effect in distinguishing between DG and NDG and could be used as a biomarker with a certain degree of accuracy. Thus, the correlation of ACE2 expression with C. albicans/A. actinomycetemcomitans had a higher predictive value for oral dysbiosis in patients with COVID-19 with diabetes.

The literature shows that C. albicans is implicated in oral diseases, including caries,² periodontitis,³⁸ denture stomatitis,³⁹ and endodontic lesions.⁴⁰ We suggest that there is a synergistic relationship between the receptor for SARS-CoV-2 entry (ACE2) and C. albicans/A. actinomycetemcomitans as the underlying mechanism of oral dysbiosis in COVID-19 patients with diabetes. These relationships may be crucial to the persistence of C. albicans and A. actinomycetemcomitans as part of the oral commensal flora and may potentially contribute to the progression of polymicrobial infection-associated dysbiosis under COVID-19 conditions.

This study has limitations. It was not possible to include all potential confounding variables. Unlike the abundance of C. albicans, we used the relative abundance of each targeted bacterium species as a proportion¹¹ rather than the actual levels. Moreover, due to the sample size is too small, the statistical power might lack of power.

Conclusions

This study revealed that the saliva of patients with COVID-19 with diabetes retained a special relationship between SARS-CoV-2 host entry and the oral dysbiotic atmosphere represented by a unique pattern of C. albicans and A. actinomycetemcomitans. This relationship could be associated with the existence of SARS-CoV-2, but it is necessary to consider the complicity of diabetes.

Data availability

Underlying data

Open Science Framework: Underlying data for ‘ACE2 expression in saliva of patients with COVID-19 and its association with Candida albicans and Aggregatibacter actinomycetemcomitans’. https://doi.org/10.17605/OSF.IO/ENFY3¹⁹

This project contains the following underlying data:

• Data file 1: Covid-19.xlsx
• Data file 2: Ct-values.xlsx

Data are available under the terms of the Creative Commons Zero “No rights reserved” data waiver (CC0 1.0 Public domain dedication).

Consent

Written informed consent for publication of the patients’ details was obtained from the patients.

Acknowledgements

All the authors would like to acknowledge all the study participants, Director of RSUI and the Ethical Clearance Committee for providing permission and ethical clearance to conduct this research project. We acknowledge the contribution provided by NR. Utami, A. Deviana, AD. Andiryani, for technical lab., and gathering data.

References

1. Janus MM, Willems HM, Krom BP: Candida albicans in Multispecies Oral Communities; A Keystone Commensal?. Adv. Exp. Med. Biol. 2016; 931: 13–20. PubMed Abstract | Publisher Full Text
2. Bachtiar EW, Bachtiar BM: Relationship between Candida albicans and Streptococcus mutans in early childhood caries, evaluated by quantitative PCR. F1000Res. 2018; 7: 1645. Epub 2019/01/11. PubMed Abstract | Publisher Full Text | Free Full Text
3. Thurnheer T, Karygianni L, Flury M, et al.: Fusobacterium Species and Subspecies Differentially Affect the Composition and Architecture of Supra- and Subgingival Biofilms Models. Front. Microbiol. 2019; 10: 1716. Epub 2019/08/17. PubMed Abstract | Publisher Full Text | Free Full Text
4. Sztukowska MN, Dutton LC, Delaney C, et al.: Community Development between Porphyromonas gingivalis and Candida albicans Mediated by InlJ and Als3. MBio. 2018; 9(2). Epub 2018/04/25. PubMed Abstract | Publisher Full Text | Free Full Text
5. Bachtiar BB, Sunarto EW, Soeroso H, et al.: ACE2 gene expression and inflammatory conditions in periodontal microenvironment of COVID-19 patients with and without diabetes evaluated by qPCR. MedRxiv [Preprint]. 2022; March 14, 2022. [accessed 2022 April 14]. Publisher Full Text
6. Bachtiar EW, Bachtiar BM: Effect of cell-free spent media prepared from Aggregatibacter actinomycetemcomitans on the growth of Candida albicans and Streptococcus mutans in co-species biofilms. Eur. J. Oral Sci. 2020; 128(5): 395–404. PubMed Abstract | Publisher Full Text
7. Periasamy S, Kolenbrander PE: Aggregatibacter actinomycetemcomitans builds mutualistic biofilm communities with Fusobacterium nucleatum and Veillonella species in saliva. Infect. Immun. 2009; 77(9): 3542–3551. PubMed Abstract | Publisher Full Text | Free Full Text
8. Bachtiar BM, Fath T, Widowati R, et al.: Quantification and Pathogenicity of Candida albicans in Denture-Wearing and Nondenture-Wearing Elderly. Eur J Dent. 2020; 14(3): 423–428. PubMed Abstract | Publisher Full Text | Free Full Text
9. Feldman M, Ginsburg I, Al-Quntar A, et al.: Thiazolidinedione-8 Alters Symbiotic Relationship in C. albicans-S. mutans Dual Species Biofilm. Front. Microbiol. 2016; 7: 140. PubMed Abstract | Publisher Full Text | Free Full Text
10. Belibasakis GN, Ozturk VO, Emingil G, et al.: Soluble triggering receptor expressed on myeloid cells 1 (sTREM-1) in gingival crevicular fluid: association with clinical and microbiologic parameters. J. Periodontol. 2014; 85(1): 204–210. PubMed Abstract | Publisher Full Text
11. Navidshad B, Liang JB, Jahromi MF: Correlation coefficients between different methods of expressing bacterial quantification using real time PCR. Int. J. Mol. Sci. 2012; 13(2): 2119–2132. PubMed Abstract | Publisher Full Text | Free Full Text
12. Alonso GC, Pavarina AC, Sousa TV, et al.: A quest to find good primers for gene expression analysis of Candida albicans from clinical samples. J. Microbiol. Methods. 2018; 147: 1–13. PubMed Abstract | Publisher Full Text
13. Nailis H, Kucharikova S, Ricicova M, et al.: Real-time PCR expression profiling of genes encoding potential virulence factors in Candida albicans biofilms: identification of model-dependent and -independent gene expression. BMC Microbiol. 2010; 10: 114. Epub 2010/04/20. PubMed Abstract | Publisher Full Text | Free Full Text
14. Xu H, Zhong L, Deng J, et al.: High expression of ACE2 receptor of 2019-nCoV on the epithelial cells of oral mucosa. Int. J. Oral Sci. 2020; 12(1): 8. Epub 2020/02/26. PubMed Abstract | Publisher Full Text | Free Full Text
15. Azzi L, Carcano G, Gianfagna F, et al.: Saliva is a reliable tool to detect SARS-CoV-2. J. Infect. 2020; 81(1): e45–e50. PubMed Abstract | Publisher Full Text | Free Full Text
16. Huang N, Perez P, Kato T, et al.: SARS-CoV-2 infection of the oral cavity and saliva. Nat. Med. 2021; 27(5): 892–903. PubMed Abstract | Publisher Full Text | Free Full Text
17. Wang Z, Zhou J, Marshall B, et al.: SARS-CoV-2 Receptor ACE2 Is Enriched in a Subpopulation of Mouse Tongue Epithelial Cells in Nongustatory Papillae but Not in Taste Buds or Embryonic Oral Epithelium. ACS Pharmacol. Transl. Sci. 2020; 3(4): 749–758. PubMed Abstract | Publisher Full Text | Free Full Text
18. Mazucanti CH, Egan JM: SARS-CoV-2 disease severity and diabetes: why the connection and what is to be done?. Immun. Ageing. 2020; 17: 21. Epub 2020/07/03. PubMed Abstract | Publisher Full Text | Free Full Text
19. Bachtiar EW: ACE2 expression in saliva of patients with COVID-19 and its association with Candida albicans and Aggregatibacter actinomycetemcomitans. Open Science Framework. 2022. Publisher Full Text
20. Bachtiar EW, Dewiyani S, Surono Akbar SM, et al.: Inhibition of Candida albicans biofilm development by unencapsulated Enterococcus faecalis cps2. J Dent Sci. 2016; 11(3): 323–330. PubMed Abstract | Publisher Full Text | Free Full Text
21. Ricker A, Vickerman M, Dongari-Bagtzoglou A: Streptococcus gordonii glucosyltransferase promotes biofilm interactions with Candida albicans. J. Oral Microbiol. 2014; 6: 6. Epub 2014/02/04. PubMed Abstract | Publisher Full Text | Free Full Text
22. Xu H, Sobue T, Thompson A, et al.: Streptococcal co-infection augments Candida pathogenicity by amplifying the mucosal inflammatory response. Cell. Microbiol. 2014; 16(2): 214–231. PubMed Abstract | Publisher Full Text | Free Full Text
23. Serrano J, Lopez-Pintor RM, Ramirez L, et al.: Risk factors related to oral candidiasis in patients with primary Sjogren's syndrome. Med. Oral Patol. Oral Cir. Bucal. 2020; 25(5): e700–e705. PubMed Abstract | Publisher Full Text | Free Full Text
24. Sugio CYC, Garcia A, Albach T, et al.: Candida-Associated Denture Stomatitis and Murine Models: What Is the Importance and Scientific Evidence?. J Fungi (Basel). 2020; 6(2). Epub 2020/05/28. PubMed Abstract | Publisher Full Text | Free Full Text
25. Segrelles-Calvo G d SAGR, Llopis-Pastor E, Carrillo J, et al.: Candida spp. co-infection in COVID-19 patients with severe pneumonia: Prevalence study and associated risk factors. Respir. Med. 2021; 188: 106619. Epub 2021/09/24. PubMed Abstract | Publisher Full Text | Free Full Text
26. Delisle MS, Williamson DR, Perreault MM, et al.: The clinical significance of Candida colonization of respiratory tract secretions in critically ill patients. J. Crit. Care. 2008; 23(1): 11–17. PubMed Abstract | Publisher Full Text
27. Heyland D, Jiang X, Day AG, et al.: Serum beta-d-glucan of critically ill patients with suspected ventilator-associated pneumonia: preliminary observations. J. Crit. Care. 2011; 26(5): 536.e1–536.e9. PubMed Abstract | Publisher Full Text
28. Thaweboon S, Thaweboon B, Nakornchai S, et al.: Salivary secretory IgA, pH, flow rates, mutans streptococci and Candida in children with rampant caries. Southeast Asian J. Trop. Med. Public Health. 2008; 39(5): 893–899. PubMed Abstract
29. Dong L, Bouey J: Public Mental Health Crisis during COVID-19 Pandemic, China. Emerg. Infect. Dis. 2020; 26(7): 1616–1618. PubMed Abstract | Publisher Full Text | Free Full Text
30. Duan L, Zhu G: Psychological interventions for people affected by the COVID-19 epidemic. Lancet Psychiatry. 2020; 7(4): 300–302. PubMed Abstract | Publisher Full Text | Free Full Text
31. Granger BL: Insight into the antiadhesive effect of yeast wall protein 1 of Candida albicans. Eukaryot. Cell. 2012; 11(6): 795–805. PubMed Abstract | Publisher Full Text | Free Full Text
32. Hoyer LL, Green CB, Oh SH, et al.: Discovering the secrets of the Candida albicans agglutinin-like sequence (ALS) gene family--a sticky pursuit. Med. Mycol. 2008; 46(1): 1–15. Epub 2007/09/14. PubMed Abstract | Publisher Full Text | Free Full Text
33. Bachtiar BM, Theodorea CF, Tahapary DL, et al.: A pilot study of red complex and three genera subgingival microbiome in periodontitis subjects with and without diabetes, evaluated by MinION platform. F1000Res. 2021; 10: 79. Epub 2021/07/28. PubMed Abstract | Publisher Full Text | Free Full Text
34. Gu X, Song LJ, Li LX, et al.: Fusobacterium nucleatum Causes Microbial Dysbiosis and Exacerbates Visceral Hypersensitivity in a Colonization-Independent Manner. Front. Microbiol. 2020; 11: 1281. Epub 2020/08/01. PubMed Abstract | Publisher Full Text | Free Full Text
35. Takahashi Y, Watanabe N, Kamio N, et al.: Expression of the SARS-CoV-2 Receptor ACE2 and Proinflammatory Cytokines Induced by the Periodontopathic Bacterium Fusobacterium nucleatum in Human Respiratory Epithelial Cells. Int. J. Mol. Sci. 2021; 22(3). Epub 2021/02/13. PubMed Abstract | Publisher Full Text | Free Full Text
36. Bachtiar EW, Bachtiar BM, Jarosz LM, et al.: AI-2 of Aggregatibacter actinomycetemcomitans inhibits Candida albicans biofilm formation. Front. Cell. Infect. Microbiol. 2014; 4: 94. Epub 2014/08/08. PubMed Abstract | Publisher Full Text | Free Full Text
37. Henderson B, Ward JM, Ready D: Aggregatibacter (Actinobacillus) actinomycetemcomitans: a triple A* periodontopathogen?. Periodontol. 2010; 54(1): 78–105. PubMed Abstract | Publisher Full Text
38. Reynaud AH, Nygaard-Ostby B, Boygard GK, et al.: Yeasts in periodontal pockets. J. Clin. Periodontol. 2001; 28(9): 860–864. PubMed Abstract | Publisher Full Text
39. Altarawneh S, Bencharit S, Mendoza L, et al.: Clinical and histological findings of denture stomatitis as related to intraoral colonization patterns of Candida albicans, salivary flow, and dry mouth. J. Prosthodont. 2013; 22(1): 13–22. PubMed Abstract | Publisher Full Text | Free Full Text
40. Baumgartner JC, Watts CM, Xia T: Occurrence of Candida albicans in infections of endodontic origin. J. Endod. 2000; 26(12): 695–698. PubMed Abstract | Publisher Full Text

Comments on this article Comments (0)

Version 2

VERSION 2 PUBLISHED 23 May 2022