Health and saliva microbiomes of a semi-urbanized indigenous tribe in Peninsular Malaysia

Li-Fang Yeo; Farhang F. Aghakhanian; James S. Y. Tan; Han Ming Gan; Maude E. Phipps

doi:10.12688/f1000research.17706.3

Home Browse Health and saliva microbiomes of a semi-urbanized indigenous tribe...

ALL Metrics

Views

Downloads

Get PDF

Get XML

Export

▬

✚

Research Article

Revised

Health and saliva microbiomes of a semi-urbanized indigenous tribe in Peninsular Malaysia

[version 3; peer review: 2 approved]

Li-Fang Yeo^1,2, Farhang F. Aghakhanian^1,2, James S. Y. Tan³, Han Ming Gan⁴, Maude E. Phipps¹

Li-Fang Yeo^1,2, Farhang F. Aghakhanian^1,2, [...] James S. Y. Tan³, Han Ming Gan⁴, Maude E. Phipps¹

PUBLISHED 28 May 2019

Author details Author details

¹ Tropical Medicine and Biology Platform, Monash University, Bandar Sunway, Selangor, 46150, Malaysia
² Jeffrey Cheah School of Medicine and Health Sciences, Monash University, Bandar Sunway, Selangor, 46150, Malaysia
³ Imperial College London, South Kensington Campus, London, SW7 2AZ, UK
⁴ School of Life & Env Sciences, Deakin University, Geelong Waurn Ponds Campus, Australia

Li-Fang Yeo
Roles: Formal Analysis, Investigation, Visualization, Writing – Original Draft Preparation, Writing – Review & Editing

Farhang F. Aghakhanian
Roles: Formal Analysis, Investigation, Resources, Validation, Writing – Review & Editing

James S. Y. Tan
Roles: Data Curation, Investigation, Validation, Writing – Review & Editing

Han Ming Gan
Roles: Funding Acquisition, Investigation, Methodology, Resources, Supervision, Writing – Review & Editing

Maude E. Phipps
Roles: Conceptualization, Funding Acquisition, Methodology, Project Administration, Resources, Supervision, Validation, Writing – Review & Editing

OPEN PEER REVIEW

REVIEWER STATUS

Abstract

Background: The indigenous people of Peninsular Malaysia, also known as Orang Asli, have gradually been urbanized. A shift towards non-communicable diseases commonly associated with sedentary lifestyles have been reported in many tribes. This study engaged with a semi-urbanized Temiar tribe from Kampong Pos Piah, Perak, who are experiencing an epidemiological transition.
Methods: Weight, height, waist circumference, blood pressure, HbA1C and lipid levels were measured as indicators of cardio-metabolic health. DNA was extracted from saliva using salting-out method followed by PCR amplification of the V3-V4 region of the 16S rRNA gene and sequencing on Illumina MiSeq. Microbiome analysis was conducted on Qiime v1.9. Statistical analysis was conducted using Qiime v1.9 and R.
Results: The study revealed that 60.4% of the Temiar community were overweight/obese, with a higher prevalence among women. HbA1C levels showed that 45% of Temiar had pre-diabetes. Insulin resistance was identified in 21% of Temiar by using a surrogate marker, TG/HDL. In total, 56.5% of Temiar were pre-hypertensive, and the condition was prevalent across all age-groups. The saliva microbiome profiles of Temiar revealed significant differences by gender, BMI, abdominal obesity as well as smoking status. The relative abundance of the genus Bifidobacterium was increased in men whereas the genera Prevotella, Capnocytophaga, Leptotrichia, Neisseria and Streptococcus were increased in women. Proteobacteria was significantly depleted in smokers.
Conclusions: Temiar from Pos Piah had a high prevalence of cardio-metabolic risks, including general and abdominal obesity, pre-diabetes, prehypertension and hypertension. This phenomenon has not been previously reported in this tribe. The saliva microbiome profiles were significantly different for individuals of different gender, BMI, abdominal obesity and smoking status.

Keywords

Orang Asli, saliva microbiome, anthropometrics, cardio-metabolic health, indigenous people

Corresponding author: Maude E. Phipps

Competing interests: No competing interests were disclosed.

Grant information: This study was funded by Tropical Medicine and Biology Platform (TMB), Monash University Malaysia seed grant TMB -2016-HG3185140920-YLF/MP awarded to LFY and MEP. Ms Li-Fang Yeo is supported by a scholarship awarded by Jeffrey Cheah School of Medicine and Health Sciences, Monash University Malaysia.
The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Copyright: © 2019 Yeo LF 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: Yeo LF, Aghakhanian FF, Tan JSY et al. Health and saliva microbiomes of a semi-urbanized indigenous tribe in Peninsular Malaysia [version 3; peer review: 2 approved]. F1000Research 2019, 8:175 (https://doi.org/10.12688/f1000research.17706.3) First published: 11 Feb 2019, 8:175 (https://doi.org/10.12688/f1000research.17706.1) Latest published: 28 May 2019, 8:175 (https://doi.org/10.12688/f1000research.17706.3)

Revised Amendments from Version 2

The figures and table titles have been edited to be more informative. R2-value and pseudo-F statistic have been added to the Results section. The new supplementary material has already been uploaded to Figshare.

See the authors' detailed response to the review by Andres Gomez

Introduction

The Orang Asli (OA), which means “original people” in the Malay language, comprise approximately 0.5% (150,000) of the total Malaysian population¹. They are categorized into three main groups, namely Negrito, Senoi and Proto Malay. OAs are widely spread across the Peninsular and range from semi-nomadic deep forest hunter-gatherers such as the Jahai to resettled communities such as Mah Meri to urbanized city-fringe dwellers such as Orang Seletar¹. This study focused on the Temiar who are a subtribe of Senoi and are believed to be descendants of the first Neolithic farmers who migrated to the Malay Peninsula².

In recent years, many OA communities were resettled by the government in effort to improve their lives. As the OA became more urbanized and by large left their ancestral habitats and practices, they led more sedentary lifestyles. These factors, coupled with loss of access to forest resources and increasing pressures to turn to store-bought food, may largely explain the rise in cardio-metabolic diseases such as hypertension, diabetes and obesity shown in recent studies^1,3,4.

The launch of the Human Microbiome Project heralded the unprecedented investigations of various microbiomes. Of these, oral microbiomes had been widely studied in human health and diseases^5,6. Studies implied an oral origin to systemic diseases such as cardiovascular diseases as the oral cavity is a major gateway into the body^7,8. Studies also investigated associations between oral microbiome and diabetes⁹, and obesity¹⁰, with mixed results. There were suggestions that obese people may have a different salivary bacterial composition perhaps akin to inflammation, which contributed to periodontal diseases and caries¹¹.

Little is known about the microbiomes of indigenous communities in Asia. To our knowledge, this was a pioneering investigation of their saliva microbiomes. Furthermore, biomedical studies of Temiar were sparse and outdated, despite them being a very large community. With this impetus, our study aimed to address the gap in knowledge by reporting on the anthropometrics and cardio-metabolic health of a resettled Temiar community and investigated their saliva microbiome in association with their health.

Results

Anthropometrics and cardio-metabolic health

A total of 72 Temiars, 33 men and 39 women, participated in the study. The median age was 34 years old. General and abdominal obesity had higher prevalence among Temiar women (Table 1). Notably, 71.4% (n=25) of women and 28.1% (n=9) of men displayed abdominal obesity.

Table 1. Anthropometrics and cardio-metabolic risk factors among Temiar.

Variable	Whole, n (%)	Men, n (%)	Women, (%)
General obesity (BMI)
Underweight (<18.5)	7 (10.3%)	4 (12.5%)	3 (8.3%)
Normal (18.5-22.9)	20 (29.4%)	12 (37.5%)	8 (22.2%)
Overweight (23-24.9)	15 (22.1%)	6 (18.8%)	9 (25%)
Pre-obese (25-29.9)	22 (32.4%)	10 (31.3%)	12 (33.3%)
Obese (≥30)	4 (5.9%)	0 (0%)	4 (11.1%)
Abdominal obesity (waist circumference)
Healthy (M: <90cm; W: <80cm)	-	23 (71.9%)	10 (28.6%)
Risk (M: ≥90cm; W: ≥80cm)	-	9 (28.1%)	25 (71.4%)
HbA1C level
Normal (4.0-5.6%)	36 (52.2%)	15 (46.9%)	21 (56.8%)
Pre-diabetes (5.7-6.4%)	31 (44.9%)	17 (53.1%)	14 (37.8%)
Diabetes	2 (2.9%)	0 (0%)	2 (5.4%)
Insulin resistance (TG/HDL 0.9-1.7)
Normal	51 (78%)	20 (75%)	32 (94%)
Risk	14 (22%)	10 (15%)	4 (6%)
Blood pressure
Normotensive	14 (20.3%)	7 (21.9%)	7 (18.9%)
Pre-hypertensive	39 (56.5%)	14 (43.8%)	24 (64.9%)
Stage 1	12 (17.4%)	9 (28.1%)	3 (8.1%)
Stage 2	4 (5.8%)	2 (6.3%)	3 (8.1%)

HbA1C levels indicated 44.9% of Temiar to be pre-diabetic, with a higher prevalence in men. The high prevalence of prediabetes is worrying because it indicates a rise in non-communicable diseases that was previously of low prevalence in rural communities. Using a TG/HDL as a surrogate marker for insulin resistance, 22% of Temiar were at risk of IR, mostly affecting men.

Blood pressure measurements showed that 56.5% (n=39) had pre-hypertension, which was more prevalent among women and was prevalent across all age groups. Stage 1 hypertension prevalence rate was 17.4% (n=12) and was found more prevalent among men. Raw data for these measurements are available in figshare¹².

Saliva microbiome analysis

To analyse the saliva microbiota diversity, the V3-V4 hypervariable region on the 16S rRNA gene was amplified and sequenced on Illumina MiSeq. After data quality control (QC), a total of 991,006 reads with mean 14,362±78 reads per individual remained.

To investigate whether the samples were sequenced to a sufficient depth, a rarefaction curve was plotted using the alpha diversity metric, Shannon index. Each colour represents a sample. The rarefaction curve indicated that all 69 samples were sequenced to a sufficient depth (Figure 1). Reads were aligned to Greengenes database V13. The major OTUs (Operational Taxonomic Unit) at the phyla level observed include Actinobacteria, Bacteroidetes, Firmicutes, Fusobacteria and Proteobacteria (Figure 2). These are the features of common oral microbiomes¹³. OTUs shall henceforth be referred to as bacterial species for the ease of reading.

Figure 1. Rarefaction curve plotted using alpha diversity metric, Shannon index against number of sequences per sample. Each coloured line represents one biological sample.

Figure 2. Relative abundance of bacteria found in the oral microbiome of Temiar at the phylum level.

PERMANOVA was used to investigate the saliva microbiome compositions using UniFrac¹⁴ distance matrix and were found to be associated with gender, obesity, waist circumference and smoking habits. To determine whether the relative abundance of individual bacterial species was differently represented in association with the factors investigated, we used the Kruskal-Wallis test with the OTU table as input.

Weighted UniFrac which takes into consideration the abundance of bacteria species revealed that the salivary bacteria were not significantly different between the two genders using PERMANOVA (p-value = 0.165, pseudo-F = 1.546, r²=0.02). However, we found that the saliva microbiomes differed significantly between men and women for unweighted UniFrac (p-value = 0.028, pseudo-F = 1.824, r²=0.02; Figure 3a). Unweighted UniFrac is a qualitative distance matrix that considers only the presence/absence of bacteria species. Kruskal-Wallis test revealed that the relative abundance of the genera Prevotella, Capnocytophaga, Leptotrichia, Neisseria and Streptococcus were significantly increased in women’s saliva microbiomes (Table 2). These commensal oral bacteria may become opportunistic pathogens in immuno-compromised states⁸. The relative abundance of the genus Bifidobacterium, was found to be highly elevated in men.

Figure 3. PCoA plots with the larger shapes representing the mean UniFrac distance between samples segregated by (a) gender – unweighted UniFrac; (b) body mass index (BMI) – unweighted UniFrac; (c) smoking status – weighted UniFrac.

Table 2. Relative abundance of bacteria (genus level) differently represented in Temiar using Kruskal-Wallis test.

Bacteria (genus level)	FDR (q-value)	Abundance
Prevotella	0.006	More abundant in women
Capnocytophaga	0.021
Leptotrichia	0.026
Neisseria	0.026
Streptococcus	0.036
Bifidobacterium	0.043	More abundant in men

The saliva microbiome profiles also differed significantly with BMI (weighted UniFrac, p-value = 0.015, pseudo-F = 2.089, r2= 0.065; unweighted UniFrac, p-value = 0.029, pseudo-F = 1.989, r²= 0.07; Figure 3b). A post-hoc Dunn’s test was conducted using the OTU table to test for differences between groups as PERMANOVA does not conduct pairwise comparison. The oral microbiome profiles of underweight individuals differed significantly from both overweight and normal individuals (underweight vs normal group, p-value = 0.0179; underweight vs overweight group, p-value = 0.0007). There was no significant difference between the saliva microbiome profiles of normal and overweight individuals (p-value = 0.0819). Differential abundance testing using Kruskal-Wallis test revealed that none of the bacterial taxa were significantly different.

There was a significant difference in the saliva microbiome of Temiar who had a healthy waist circumference compared to those with abdominal obesity (weighted UniFrac, p-value = 0.022, pseudo-F = 2.289, r² = 0.05). However, there was no significant difference in unweighted UniFrac (p-value = 0.286, pseudo-F = 1.099, r² = 0.021) as well as in individual bacterial taxa (p-value>0.05) among healthy individuals and those with abdominal obesity. The saliva microbiome composition of non-diabetic, pre-diabetic and diabetic individuals suggested some differences, but they were not significant (unweighted p-value = 0.069, pseudo-F = 1.502, r² = 0.043; weighted p-value = 0.122, pseudo-F = 1.579, r² = 0.045). The saliva microbiome composition and relative abundance of specific bacterial species were not statistically different when categorised by age group, lipid levels nor blood pressure levels.

There was a perceptible difference in the saliva microbiomes and smoking habits using weighted UniFrac (p-value = 0.016, pseudo-F = 2.498, r² = 0.07; Figure 3c) but no difference was detected when using unweighted UniFrac (p-value = 0.059, pseudo-F = 1.475, r² = 0.04). Further testing with Kruskal-Wallis showed the relative abundance of Proteobacteria and Firmicutes (phylum level) were significantly different among smokers and non-smokers. Differential abundance testing at the genus level revealed the relative abundance of Neisseria and Aggregatibacter was decreased in smokers compared to never-smokers and former smokers. Current smokers had a lower abundance of the genus Neisseria and Aggregatibacter than former smokers, but the difference was not statistically significant (p-value>0.05). The relative abundance of the genus Campylobacter and the class Clostridia, were higher in both current and former smokers compared to never-smokers (Table 3).

Table 3. Relative abundance of bacteria found among smokers vs former and never-smokers using Kruskal-Wallis test.

Bacteria species		FDR(q-value)	Relative abundance of bacteria
Genus	Neisseria	0.025	Depleted in smokers and former smokers
	Aggregatibacter	0.035	Depleted in smokers and former smokers
	Campylobacter	0.037	Increased in smokers and former smokers
Class	Clostridia	0.037	Increased in smokers and former smokers

Overall, the relative abundance of the genera Neisseria and Aggregatibacter was decreased in current and former smokers, whereas the relative abundance of the genus Campylobacter and the class Clostridia was greater in smokers. The saliva microbiome showed no significant difference between former smokers and never-smokers

Discussion

We reported a high prevalence of cardio-metabolic diseases such as obesity, pre-diabetes, insulin resistance and pre-hypertension among Temiar. These non-communicable diseases were previously not reported in OA^15,16, but recent studies have indicated their high prevalence, especially in OA tribes living near cities^1,3,17. Increased cardio-metabolic risks were also reported in aboriginal Torres Straits Islanders from Australia¹⁸, the Jaguapiru indigenous community in Brazil¹⁹ and the Rang Bothiya tribe in India²⁰.

Many OA tribes lead relatively more sedentary lifestyles compared to their hunter-gatherer ancestors. They can no longer and perhaps have no need to rely entirely on the forest and its resources for survival. Rapid development had given them easier access to high-calorie processed foods, which may have contributed to obesity and cardio-metabolic diseases.

National Health and Morbidity Survey 2015 reported obesity prevalence among the major ethnic groups in Malaysia to be 30.6%, which is comparable to the findings in this study (Temiar obesity = 32.4%); the national prevalence rate of diabetes was reported to be 22.9%²¹, which is very much higher than that of the Temiar (2.9%) in this study. Generally, reported prevalence rates of obesity among OA are still low. Nonetheless, the high prevalence rate of pre-diabetes reported among the Temiar (44.9%) indicates that rural OA communities are in dire need of awareness education and medical intervention.

Saliva microbiome analysis revealed significant difference in microbial composition among men and women, where the genus Prevotella was significantly higher in women. Hitherto, no studies have reported increased prevalence in women’s saliva microbiome. There are studies that suggest that hormones may play a role in the type of bacteria that colonise women’s oral cavity^22–26 Temiar women had a lower prevalence of pre-diabetes and insulin resistance, despite the majority of them presenting with either general or abdominal obesity. Perhaps the women consumed a traditional, indigenous diet, which is richer in plant-fibre and less meat compared to men. Although Prevotella is a naturally occurring member of the oral microbiota²⁷, it is also associated with inflammatory conditions such as rheumatoid arthritis and periodontal infections²⁸.

The relative abundance of the genus Bifidobacterium, was shown to be elevated in Temiar men. While it is uncertain whether Temiar men were exposed to more dairy products than women, food taboos practiced among Temiar may contribute to the differences observed among gender²⁹. Several studies that investigated oral microbiomes of urbanized cohorts in association with gender have reported no differences in oral microbiome profiles^23–25. This may be explained by the relatively homogenous environment that urbanized cohorts were exposed to. Studies have shown that salivary microbiomes are most affected by environmental factors, as the oral microbiome of twins which were similar became highly dissimilar when they lived apart²³.

Temiar, on the other hand, lived in a traditional environment where men and women had different social standings. Men went out to the forest to hunt or forage while women stayed in the village with the children. They also observed certain food taboos, where the bush meat consumption of animals such as river terrapin, gibbons and porcupine were reserved only for men²⁹.

This preliminary investigation suggested links between saliva microbiomes and gender where differences may be attributed to cultural, dietary and environmental factors. Even though the bacteria driving the differences in obesity and gender were of different species, it should be noted that most of the women were overweight/obese, which could be a confounding factor in gender-driven disparities in the saliva microbiome.

Studies have suggested an association between obesity and altered oral microbiome^26,30,31, concurring with the findings our study. However, a significant difference was noted only when comparing overweight and underweight individuals. Both states are considered to be ‘unhealthy’ and thus assumed to be at dysbiosis.

The oral hygiene practices and oral health among Temiar were unknown, although due to their geographical isolation, it was highly unlikely they have regular access to dental health care. Our study revealed that the relative abundance of the phylum Proteobacteria, and the genera Neisseria and Aggregatibacter, were decreased in smokers, compared to non-smokers. This was in line with the findings of Wu et al. (2016)³², although our study did not detect differences in the several other genera reported by them. This was probably due to the differences in sample size of both studies.

Wu et al. (2016) reported that Proteobacteria were associated with the breakdown of toxic hydrocarbons found in cigarette smoke, hence a depletion of the bacteria genus in smokers may prove detrimental to oral health³². An interesting difference noted was an increase in pathogenic anaerobe Streptococcus observed in their study³², whereas Hernandez et al. (2017)³³ reported depletion of Streptococcus among betel-nut chewers, even after controlling for smoking. Our investigations revealed no significant difference, even though many of the Temiar smokers were also self-reported betel-nut chewers, a practice frequently associated with oral cancer³⁴. Further investigations may be required to distinguish the effects of betel-nut chewing and smoking on the saliva microbiome.

Conclusion

Our study revealed a high prevalence of cardio-metabolic diseases among the Temiar, including general and abdominal obesity, pre-diabetes and insulin resistance. Pre-hypertension was found highly prevalent across all age groups.

Interestingly, the saliva microbiome profiles were significantly different for gender where the relative abundance of the genera Prevotella, Capnocytophaga, Leptotrichia, Neisseria, Streptococcus and Bifidobacterium were concerned. Our study also noted a significant difference between the saliva microbiome compositions of underweight vs overweight and normal individuals. The oral microbiome was not significantly different among non-diabetic, pre-diabetic and diabetic individuals. The microbiome profiles differed significantly among smokers and non-smokers where further investigation showed that the phylum Proteobacteria were significantly decreased in smokers. Investigation towards other health parameters such as pre-diabetes were inconclusive.

The sample size was small but sufficient to illustrate the health problems that plague the Temiar, especially the high prevalence of obesity, pre-diabetes and pre-hypertension. It would be very useful if future studies could determine the factors that are contributing to these health problems. Could it be a lack of awareness, overly sedentary lifestyle or perhaps a genetic predisposition? Regardless, studies have shown time and again that the OA are facing health problems that are reversible if they are made aware of them and the taught the correct methods to improve their lifestyle.

Methods

Data and sample collection

The study was approved by Ministry of Health Malaysia under National Medical Research Registry, MNDR ID #09—23-3913, Department of Orang Asli Development, Malaysia (JAKOA) and Monash University Human Research Ethics Committee (MUHREC). Before sampling, a courtesy visit to the Temiar elders in Kampong Pos Piah, Perak was conducted to explain the rationale of the study. Upon agreement to participate in our study, a medical team returned to the village on an agreed date and conducted health screening and sampling. Participants who were over 18 years old with no visible health ailments and able to provide informed consent were recruited for the study through convenience sampling, that is whoever who turned up and was eligible. Participants who were pregnant, with a history of alcohol/drug abuse, or with chronic illness (e.g. kidney failure, cancer, heart disease) were excluded from the study.

The consent form was read aloud by interviewers and queries were addressed before either a signature or thumbprint was provided as a sign of consent. A total of 72 Temiar provided informed consent to participate. Interviews were conducted in Bahasa Malaysia using a questionnaire¹² to collect information about their socio-demography, medical history and diet. Height, weight, waist circumference and blood pressure were measured¹. Participants were also examined by clinicians. Acanthosis negricans, which is darkening of the skin around the neck and creases of elbows indicative of insulin resistance, was noted.

Saliva samples were collected in sterile 50ml polypropylene Falcon tube. Participants were requested to rinse their mouths with water thoroughly 30 minutes prior to collecting saliva. Venous blood samples were taken for biochemical analyses.

Anthropometrics and Biochemical analysis

BMI, waist circumference and blood pressure cut-off values were in accordance to WHO recommendations³⁵. We measured their HbA1C and blood lipid levels (cholesterol, HDL, LDL, Triglyceride). We used TG/HDL ratio as a surrogate marker for insulin resistance with a cut-off value of 0.9-1.7³⁶.

DNA extraction and PCR

DNA was extracted from saliva using a modified high salt-method³⁷. The V3-V4 region of the 16S rRNA gene were targeted, resulting in a PCR product of approximately 550 bp³⁸.

Sequencing on Illumina MiSeq

DNA sequencing was done by Genomics Facility in Monash University Malaysia on Illumina MiSeq to produce paired end reads of approximately 230 bp each.

Data analysis

Microbiome analysis was conducted on QIIME 1.9³⁹. Chimeras were filtered using UCHIME 1.39.3⁴⁰ before being aligned to Greengenes database V13.8⁴¹. The reads were then clustered into operational taxonomic units (OTUs) with open-reference method at 97% similarity level using UCLUST⁴² in the QIIME pipeline. OTU clusters were assigned taxonomy with RDP classifier⁴³. The reads were normalized and OTUs that were present at less than 0.05% were filtered off.

Alpha diversity and beta diversity of the samples were reported using phylogenetic distance (PD) and UniFrac¹⁴, respectively. PCoA plots were generated to visualize beta diversity of the samples.

Statistical analysis

Statistical analyses were completed on QIIME and R 3.4.4. Information taken from the mapping file included gender, BMI and smoking status. PERMANOVA, a non-parametric test was used to test for differences in median among the groups using weighted and unweighted UniFrac distance matrix Usi//ng the R packages vegan (v2.4-2), readr (v1.1.0) and dplyr (v0.5.0). Kruskal-Wallis test was used to test for differences in the relative abundance of OTUs among the different groups. A post-hoc test, Dunn’s test was done for pairwise comparison when testing factors like BMI and smoking, as they had more than two groups. False discovery rate, reported as q-value, was used to control for multiple hypothesis testing and was statistically significant at 5%.

Data availability

Underlying data

Saliva microbiomes of the individuals in this study are available from the Sequence Read Archive, BioProject accession number PRJNA515166; https://identifiers.org/bioproject/PRJNA515166.

Anthropometric data, along with the other variables measured, as well as supplementary data are available on figshare. DOI: https://doi.org/10.26180/5c453ff435883¹².

Extended data

The questionnaire used in this study is available on figshare. DOI: https://doi.org/10.26180/5c453ff435883¹².

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

Grant information

This study was funded by Tropical Medicine and Biology Platform (TMB), Monash University Malaysia seed grant TMB -2016-HG3185140920-YLF/MP awarded to LFY and MEP. Ms Li-Fang Yeo is supported by a scholarship awarded by Jeffrey Cheah School of Medicine and Health Sciences, Monash University Malaysia.

The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Acknowledgements

We thank the Temiar community for their participation and the Department of Orang Asli Development (JAKOA) of Malaysia. We are especially grateful to Prof. Sunil K. Lal for helpful suggestions, Ms Yin Peng Lee for technical advice, and Dr Christina Yap for conducting biochemical analysis, Prof Daniel Reidpath for his assistance in statistical analysis on R, Dr Amreeta Dhanoa for editing the manuscript, Prof Yvonne Lim for advice and Dr Eustacia Lee and the University of Malaya team for participating in sample collection. Furthermore, our special thanks to Dr Badariah Ahmad and Dr Siti Noraida for participating in fieldwork.

Faculty Opinions recommended

References

1. Phipps ME, Chan KK, Naidu R, et al.: Cardio-metabolic health risks in indigenous populations of Southeast Asia and the influence of urbanization. BMC Public Health. 2015; 15: 47. PubMed Abstract | Publisher Full Text | Free Full Text
2. McColl H, Racimo F, Vinner L, et al.: The prehistoric peopling of Southeast Asia. Science. 2018; 361(6397): 88–92. PubMed Abstract | Publisher Full Text
3. Wong CY, Zalilah MS, Chua EY, et al.: Double-burden of malnutrition among the indigenous peoples (Orang Asli) of Peninsular Malaysia. BMC Public Health. 2015; 15: 680. PubMed Abstract | Publisher Full Text | Free Full Text
4. Tuan Abdul Aziz TA, Teh LK, Md Idris MH, et al.: Increased risks of cardiovascular diseases and insulin resistance among the Orang Asli in Peninsular Malaysia. BMC Public Health. 2016; 16: 284. PubMed Abstract | Publisher Full Text | Free Full Text
5. Olsen I: The oral microbiome in health and disease. Infections and General Health. 2016; 97–114. Publisher Full Text
6. Yamashita Y, Takeshita T: The oral microbiome and human health. J Oral Sci. 2017; 59(2): 201–206. PubMed Abstract | Publisher Full Text
7. Sudhakara P, Gupta A, Bhardwaj A, et al.: Oral Dysbiotic Communities and Their Implications in Systemic Diseases. Dent J (Basel). 2018; 6(2): pii: E10. PubMed Abstract | Publisher Full Text | Free Full Text
8. Idris A, Hasnain SZ, Huat LZ, et al.: Human diseases, immunity and the oral microbiota—Insights gained from metagenomic studies. Oral Science Int. 2017; 14(2): 27–32. Publisher Full Text
9. Long J, Cai Q, Steinwandel M, et al.: Association of oral microbiome with type 2 diabetes risk. J Periodontal Res. 2017; 52(3): 636–643. PubMed Abstract | Publisher Full Text | Free Full Text
10. Zeigler CC, Persson GR, Wondimu B, et al.: Microbiota in the oral subgingival biofilm is associated with obesity in adolescence. Obesity (Silver Spring). 2012; 20(1): 157–164. PubMed Abstract | Publisher Full Text
11. Choromańska K, Choromańska B, Dąbrowska E, et al.: Saliva of obese patients - is it different? Postepy Hig Med Dosw (Online). 2015; 69: 1190–1195. PubMed Abstract | Publisher Full Text
12. Yeo LF, Phipps M, Aghakhanian FF, et al.: Health and Saliva Microbiomes of a Semi-urbanized Indigenous Tribe in Peninsular Malaysia. figshare. Fileset. 2019. http://www.doi.org/10.26180/5c453ff435883
13. Dewhirst FE, Chen T, Izard J, et al.: The human oral microbiome. J Bacteriol. 2010; 192(19): 5002–5017. PubMed Abstract | Publisher Full Text | Free Full Text
14. Lozupone C, Knight R: UniFrac: a new phylogenetic method for comparing microbial communities. Appl Environ Microbiol. 2005; 71(12): 8228–8235. PubMed Abstract | Publisher Full Text | Free Full Text
15. Burns-Cox CJ: Splenomegaly and blood pressure in an Orang Asli community in West Malaysia. Am Heart J. 1970; 80(5): 718–719. PubMed Abstract | Publisher Full Text
16. Burns-Cox CJ, Chong YH, Gillman R: Risk factors and the absence of coronary heart disease in aborigines in West Malaysia. Br Heart J. 1972; 34(9): 953–8. PubMed Abstract | Publisher Full Text | Free Full Text
17. Pell C, Allotey P, Evans N, et al.: Coming of age, becoming obese: a cross-sectional analysis of obesity among adolescents and young adults in Malaysia. BMC Public Health. 2016; 16(1): 1082. PubMed Abstract | Publisher Full Text | Free Full Text
18. Burns J, Thomson N: Summary of overweight and obesity among Indigenous peoples. Retrieved [03 November 2015]; 2006.
19. Oliveira GF, Oliveira TR, Ikejiri AT, et al.: Prevalence of Obesity and Overweight in an Indigenous Population in Central Brazil: A Population-Based Cross-Sectional Study. Obes Facts. 2015; 8(5): 302–310. PubMed Abstract | Publisher Full Text | Free Full Text
20. Kandpal V, Sachdeva MP, Saraswathy KN: An assessment study of CVD related risk factors in a tribal population of India. BMC Public Health. 2016; 16(1): 434. PubMed Abstract | Publisher Full Text | Free Full Text
21. Institute of Public Health: National Health and Morbidity Survey. Institute of Public Health; 2015. Reference Source
22. Kumar PS: Sex and the subgingival microbiome: do female sex steroids affect periodontal bacteria? Periodontol 2000. 2013; 61(1): 103–124. PubMed Abstract | Publisher Full Text
23. Stahringer SS, Clemente JC, Corley RP, et al.: Nurture trumps nature in a longitudinal survey of salivary bacterial communities in twins from early adolescence to early adulthood. Genome Res. 2012; 22(11): 2146–2152. PubMed Abstract | Publisher Full Text | Free Full Text
24. Chen H, Jiang W: Application of high-throughput sequencing in understanding human oral microbiome related with health and disease. Front Microbiol. 2014; 5: 508. PubMed Abstract | Publisher Full Text | Free Full Text
25. Li J, Quinque D, Horz HP, et al.: Comparative analysis of the human saliva microbiome from different climate zones: Alaska, Germany, and Africa. BMC Microbiol. 2014; 14(1): 316. PubMed Abstract | Publisher Full Text | Free Full Text
26. Takeshita T, Kageyama S, Furuta M, et al.: Bacterial diversity in saliva and oral health-related conditions: the Hisayama Study. Sci Rep. 2016; 6(1): 22164. PubMed Abstract | Publisher Full Text | Free Full Text
27. Palmer RJ Jr: Composition and development of oral bacterial communities. Periodontol 2000. 2014; 64(1): 20–39. PubMed Abstract | Publisher Full Text | Free Full Text
28. Larsen JM: The immune response to Prevotella bacteria in chronic inflammatory disease. Immunology. 2017; 151(4): 363–374. PubMed Abstract | Publisher Full Text | Free Full Text
29. Benjamin G: Temiar Religion, 1964-2012: Enchantment, Disenchantment and Re-enchantment in Malaysia’s Uplands. 2014. Reference Source
30. Goodson JM, Hartman ML, Shi P, et al.: The salivary microbiome is altered in the presence of a high salivary glucose concentration. PLoS One. 2017; 12(3): e0170437. PubMed Abstract | Publisher Full Text | Free Full Text
31. Janem WF, Scannapieco FA, Sabharwal A, et al.: Salivary inflammatory markers and microbiome in normoglycemic lean and obese children compared to obese children with type 2 diabetes. PLoS One. 2017; 12(3): e0172647. PubMed Abstract | Publisher Full Text | Free Full Text
32. Wu J, Peters BA, Dominianni C, et al.: Cigarette smoking and the oral microbiome in a large study of American adults. ISME J. 2016; 10(10): 2435–2446. PubMed Abstract | Publisher Full Text | Free Full Text
33. Hernandez BY, Zhu X, Goodman MT, et al.: Betel nut chewing, oral premalignant lesions, and the oral microbiome. PLoS One. 2017; 12(2): e0172196. PubMed Abstract | Publisher Full Text | Free Full Text
34. Guha N, Warnakulasuriya S, Vlaanderen J, et al.: Betel quid chewing and the risk of oral and oropharyngeal cancers: a meta-analysis with implications for cancer control. Int J Cancer. 2014; 135(6): 1433–43. PubMed Abstract | Publisher Full Text
35. World Health Organization: Use of glycated haemoglobin (HbA1c) in the diagnosis of diabetes mellitus. Abbreviated report of a WHO consultation 2011. Geneva: World Health Organization. 2013. Reference Source
36. Mostafa SA, Davies MJ, Morris DH, et al.: The association of the triglyceride-to-HDL cholesterol ratio with insulin resistance in White European and South Asian men and women. PLoS One. 2012; 7(12): e50931. PubMed Abstract | Publisher Full Text | Free Full Text
37. Quinque D, Kittler R, Kayser M, et al.: Evaluation of saliva as a source of human DNA for population and association studies. Anal Biochem. 2006; 353(2): 272–277. PubMed Abstract | Publisher Full Text
38. PCR I: 16S Metagenomics Sequencing Library Preparation. Illumina. 2013. Reference Source
39. Caporaso JG, Kuczynski J, Stombaugh J, et al.: QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010; 7(5): 335–336. PubMed Abstract | Publisher Full Text | Free Full Text
40. Edgar RC, Haas BJ, Clemente JC, et al.: UCHIME improves sensitivity and speed of chimera detection. Bioinformatics. 2011; 27(16): 2194–2200. PubMed Abstract | Publisher Full Text | Free Full Text
41. DeSantis TZ, Hugenholtz P, Larsen N, et al.: Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl Environ Microbiol. 2006; 72(7): 5069–5072. PubMed Abstract | Publisher Full Text | Free Full Text
42. Edgar RC: Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010; 26(19): 2460–2461. PubMed Abstract | Publisher Full Text
43. Wang Q, Garrity GM, Tiedje JM, et al.: Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol. 2007; 73(16): 5261–5267. PubMed Abstract | Publisher Full Text | Free Full Text

Comments on this article Comments (2)

Version 3

VERSION 3 PUBLISHED 28 May 2019

Revised

Comment

Version 1

VERSION 1 PUBLISHED 11 Feb 2019

Discussion is closed on this version, please comment on the latest version above.

Author Response 09 May 2019

Li-Fang Yeo, Tropical Medicine and Biology Platform, Monash University, Bandar Sunway, 46150, Malaysia

09 May 2019

Author Response

Thank you Andres and Siti Nursheena for your comments. We have addressed all your concerns, hopefully to your satisfaction, with one exception regarding Andres's suggestion to use boxplots to show ... Continue reading Thank you Andres and Siti Nursheena for your comments. We have addressed all your concerns, hopefully to your satisfaction, with one exception regarding Andres's suggestion to use boxplots to show median and variation of data. For this publication, we have decided to maintain the results in table-form showing q-values due to time constraint and technical limitations on our side. We will definitely take this into consideration to improve our future work.
Thank you Andres and Siti Nursheena for your comments. We have addressed all your concerns, hopefully to your satisfaction, with one exception regarding Andres's suggestion to use boxplots to show median and variation of data. For this publication, we have decided to maintain the results in table-form showing q-values due to time constraint and technical limitations on our side. We will definitely take this into consideration to improve our future work.
Competing Interests: No competing interests. Close
Report a concern
Author Response 24 Apr 2019

Maude Phipps, Jeffrey Cheah School of Medicine and Health Sciences, Monash University, Bandar Sunway, 46150, Malaysia

24 Apr 2019

Author Response

Thanks Andres and Siti Nursheena for your constructive comments. We endeavor to improve this paper and hope it will be confirmed and citable soon.
Competing Interests: No competing interests were disclosed.
Thanks Andres and Siti Nursheena for your constructive comments. We endeavor to improve this paper and hope it will be confirmed and citable soon.
Thanks Andres and Siti Nursheena for your constructive comments. We endeavor to improve this paper and hope it will be confirmed and citable soon.
Competing Interests: No competing interests were disclosed. Close
Report a concern
Discussion is closed on this version, please comment on the latest version above.

Author details Author details

Li-Fang Yeo
Roles: Formal Analysis, Investigation, Visualization, Writing – Original Draft Preparation, Writing – Review & Editing

Farhang F. Aghakhanian
Roles: Formal Analysis, Investigation, Resources, Validation, Writing – Review & Editing

James S. Y. Tan
Roles: Data Curation, Investigation, Validation, Writing – Review & Editing

Han Ming Gan
Roles: Funding Acquisition, Investigation, Methodology, Resources, Supervision, Writing – Review & Editing

Maude E. Phipps
Roles: Conceptualization, Funding Acquisition, Methodology, Project Administration, Resources, Supervision, Validation, Writing – Review & Editing

Competing interests

No competing interests were disclosed.

Grant information

Article Versions (3)

version 3

Revised

Published: 28 May 2019, 8:175

https://doi.org/10.12688/f1000research.17706.3

version 2

Revised

Published: 09 May 2019, 8:175

https://doi.org/10.12688/f1000research.17706.2

version 1

Published: 11 Feb 2019, 8:175

https://doi.org/10.12688/f1000research.17706.1

© 2019 Yeo LF 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.

Download

Export To

metrics

	Views	Downloads
F1000Research	-	-
PubMed Central Data from PMC are received and updated monthly.	-	-

Citations

SEE MORE DETAILS

CITE

how to cite this article

Yeo LF, Aghakhanian FF, Tan JSY et al. Health and saliva microbiomes of a semi-urbanized indigenous tribe in Peninsular Malaysia [version 3; peer review: 2 approved] F1000Research 2019, 8:175 (https://doi.org/10.12688/f1000research.17706.3)

NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.

track

receive updates on this article

Track an article to receive email alerts on any updates to this article.

Open Peer Review

Current Reviewer Status: ?

Key to Reviewer Statuses VIEW HIDE

ApprovedThe paper is scientifically sound in its current form and only minor, if any, improvements are suggested

Approved with reservations A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit.

Not approvedFundamental flaws in the paper seriously undermine the findings and conclusions

Version 3

VERSION 3

PUBLISHED 28 May 2019

Revised

Views

Reviewer Report 29 May 2019

Andres Gomez, Department of Animal Science, University of Minnesota, Saint Paul, MN, USA

Approved

https://doi.org/10.5256/f1000research.21321.r49127

Thanks to the authors for ... Continue reading

CITE

Report a concern

Author Response 31 May 2019

Li-Fang Yeo, Tropical Medicine and Biology Platform, Monash University, Bandar Sunway, 46150, Malaysia

31 May 2019

Author Response

The authors thank the reviewer very much for his time and all comments were very much appreciated.
Competing Interests: No competing interest.
The authors thank the reviewer very much for his time and all comments were very much appreciated.
The authors thank the reviewer very much for his time and all comments were very much appreciated.
Competing Interests: No competing interest. Close
Report a concern
Respond or Comment

COMMENTS ON THIS REPORT

Author Response 31 May 2019

Li-Fang Yeo, Tropical Medicine and Biology Platform, Monash University, Bandar Sunway, 46150, Malaysia

31 May 2019

Author Response

The authors thank the reviewer very much for his time and all comments were very much appreciated.
Competing Interests: No competing interest.
The authors thank the reviewer very much for his time and all comments were very much appreciated.
The authors thank the reviewer very much for his time and all comments were very much appreciated.
Competing Interests: No competing interest. Close
Report a concern

Version 2

VERSION 2

PUBLISHED 09 May 2019

Revised

Views

Reviewer Report 22 May 2019

Andres Gomez, Department of Animal Science, University of Minnesota, Saint Paul, MN, USA

Approved with Reservations

https://doi.org/10.5256/f1000research.21016.r48249

Figure 1. Please specify what each color on the rarefaction curve represents.

Figure 3 and distance based analyses: Authors still have not followed through with adding R2 and pseudo-F values, besides P-values, in the text. I know it is referenced in the supplement (which by the way should be referenced throughout the text with specific S-table or S-figure number), but these three values should be referenced in the main text. Also, authors should specify in the text that the P-values from the distance-based differences calculated (also R2 and pseudo-F) were calculated via PERMANOVA.

Also on figure 3, please omit "...representing the mean of distance matrix for.." It should read "representing the mean UniFrac distance between samples"

PERMANOVA analyses: Authors have a confusion about the difference between then terms ADONIS and PERMANOVA- basically ADONIS, is a function within the R package vegan to perform PERMANOVA- the R2 is not a value from adonis, it's a value generated in PERMANOVA (all, P-values, R2 and pseudo-F).

This reviewer thanks the authors for improving the manuscript. However, before it can finally be endorsed, authors should address the following minor comments:

Table 3: This is a comment made in the previous version: the title "Relative abundance of bacteria found among smokers vs former and never-smokers generated by QIIME" is incorrect - QIIME is am interface with several plugins, authors should specify which of these specific plugins allowed them to detect discriminant taxa in table 3 (e.g. LeFSE, binomial distributions?). Authors should specify I the text that discriminant taxa analyses were made at the genus level. Last, if authors mix the term OTU with several taxonomic levels (e.g. OTUS at the genius level, or OTUS at the phyla level, readers may get confused- typically, we only use OTUs as a proxy for species (>97% 16S rRNA sequence similarity).

Competing Interests: No competing interests were disclosed.

I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard, however I have significant reservations, as outlined above.

CITE

Report a concern

Respond or Comment

Views

Reviewer Report 17 May 2019

Siti Nursheena Mohd Zain, Institute of Biological Science, Faculty of Science, University of Malaya, Kuala Lumpur, Malaysia

Approved

https://doi.org/10.5256/f1000research.21016.r48248

The article reads well ... Continue reading

CITE

Report a concern

Respond or Comment

Version 1

VERSION 1

PUBLISHED 11 Feb 2019

Views

Reviewer Report 10 Apr 2019

Siti Nursheena Mohd Zain, Institute of Biological Science, Faculty of Science, University of Malaya, Kuala Lumpur, Malaysia

Approved with Reservations

https://doi.org/10.5256/f1000research.19359.r46035

An interesting work looking at non-communicable disease risks link to microbiome in saliva as an indicator. The discussion could be written better to pull the results together. What is missing is a comparative study with another sedentary population. I would like if the conclusion is added with recommendations to improve the issue facing the Temiar tribe. Noted that sample size is small to make small but still warrants some finding to improve their health status.

Is the work clearly and accurately presented and does it cite the current literature?

Yes
Is the study design appropriate and is the work technically sound?

Yes
Are sufficient details of methods and analysis provided to allow replication by others?

Yes
If applicable, is the statistical analysis and its interpretation appropriate?

I cannot comment. A qualified statistician is required.
Are all the source data underlying the results available to ensure full reproducibility?

Yes
Are the conclusions drawn adequately supported by the results?

Yes

Competing Interests: No competing interests were disclosed.

Reviewer Expertise: Parasitology and molecular parasitology

CITE

Report a concern

Respond or Comment

Views

Reviewer Report 13 Mar 2019

Andres Gomez, Department of Animal Science, University of Minnesota, Saint Paul, MN, USA

Approved with Reservations

https://doi.org/10.5256/f1000research.19359.r44294

Results
The QIIME microbiome analyses should be briefly described as far as the bioinformatics pipeline used and its details: clustering of OTUS? OTU picking? amplicon sequence variants (DADA2-ASV)? These should be specifically detailed.

It does not seem that the rarefaction curve shows compete plateau, also what do the colors mean?

PERMANOVA is used to mine for distinctions in microbiome composition.

It is unclear how kruskal-wallis test was used in conjunction with Unifrac - KW is used for pairwise comparisons in multivariate space, not for multivariate data.

PERMANOVA results (P values) should be accompanied by R2 and, Pseudo F values also. PERMANOVA models used should also be reported.

"Kruskal-Wallis test revealed that the relative abundance of Prevotella, Capnocytophaga, Leptotrichia, Neisseria and Streptococcus"

Were these discriminant analyses made at genus or OTU/ASV levels? Also, it is recommended to show these results in box plots showing medians and variation of data, not in tables

Authors make several assertions as to the character and typically reported nature of the markers discovered, in the context of "good/bad bacteria" - but many of the claims made are in the context of gut microbiomes not oral communities.
For instance:
"Prevotella, associated with consumption of plant-rich diets"
This claim may be true only for the gut microbiome not for oral communities. Also: "The relative abundance of Bifidobacterium, the ‘good bacteria’ commonly found in dairy products, was found to be highly elevated in men" - this may be associated with effect of dairy products in the gut not in the oral cavity.

"These commensal oral bacteria may become opportunistic pathogens in immune-compromised states." please cite

Fig 3- legends should indicate that a-b were based on Unw-UniFrac, and that C on W-UniFrac - Please depict % variation on axes

It is not clear how the post-hoc Dun's test were conducted-0 is this still based on PERMANOVA?

"The Kruskal-Wallis test generated by QIIME revealed no significant difference." This is not clear - compared to what?

"There was a significant difference between waist circumference and saliva microbiome composition (weighted UniFrac, p-value = 0.022). " This is not clear - you mean differences between the microbiomes of subjects with different waist circumference?

QIIME is just a software with plugins to perform statistical analyses - it is not clear when authors are using qiime or R and how these platforms would yield different results. Please avoid jargon as "generated by QIIME"

"The saliva microbiome showed no significant difference between former smokers and never-smokers, suggesting perhaps the saliva microbiome may recover partially to an environment prior to smoking". This is unclear but also speculative.

The discussion is highly speculative and should be toned down, mainly as far as the hypotheses behind the differences reported.

Authors also seem to interchangeably extrapolate hypothesis on oral microbiome changes with gut microbiome changes in the context of lifestyle and diet (e.g. Prevotella) although diet may also be linked to oral microbiome modulation, this relationship is less clear, and taxa such as Prevotella, abundant in both gut and oral cavity, cannot be assumed to be modulated by diet in both sites

Is the work clearly and accurately presented and does it cite the current literature?

Partly
Is the study design appropriate and is the work technically sound?

Partly
Are sufficient details of methods and analysis provided to allow replication by others?

Partly
If applicable, is the statistical analysis and its interpretation appropriate?

Partly
Are all the source data underlying the results available to ensure full reproducibility?

Partly
Are the conclusions drawn adequately supported by the results?

Partly

Competing Interests: No competing interests were disclosed.

Reviewer Expertise: Microbial ecology

CITE

Report a concern

Author Response 09 May 2019

Li-Fang Yeo, Tropical Medicine and Biology Platform, Monash University, Bandar Sunway, 46150, Malaysia

09 May 2019

Author Response

Details on the bioinformatics pipeline will be addressed in version 2. The different colours on the rarefaction curve represent different individuals. Other metrics such as Shannon and Chao1 shows an ... Continue reading Details on the bioinformatics pipeline will be addressed in version 2. The different colours on the rarefaction curve represent different individuals. Other metrics such as Shannon and Chao1 shows an obvious plateau. We felt that PD was more informative because it included phylogenetic distances.

Our statistics were done mostly with reference to this paper:
Wu, J., Peters, B.A., Dominianni, C., Zhang, Y., Pei, Z., Yang, L., Ma, Y., Purdue, M.P., Jacobs, E.J., Gapstur, S.M. and Li, H., 2016. “Cigarette smoking and the oral microbiome in a large study of American adults”. The ISME journal. https://www.nature.com/articles/ismej201637

However, more recent methods such as using balances and log ratio shall be considered in a larger, extended study. Test scores for PERMANOVA results shall be addressed in version 2 as will several questions that were unclear in version 1.

To date, we have none found literature that was comparable (ie. Having Bifidobacterium elevated in men, Prevotella elevated in women). Hence we attempted to suggest possible explanations to this observation. The discussion will be re-worded to take into account reviewer’s comments.

Details on the bioinformatics pipeline will be addressed in version 2. The different colours on the rarefaction curve represent different individuals. Other metrics such as Shannon and Chao1 shows an obvious plateau. We felt that PD was more informative because it included phylogenetic distances.

Our statistics were done mostly with reference to this paper:
Wu, J., Peters, B.A., Dominianni, C., Zhang, Y., Pei, Z., Yang, L., Ma, Y., Purdue, M.P., Jacobs, E.J., Gapstur, S.M. and Li, H., 2016. “Cigarette smoking and the oral microbiome in a large study of American adults”. The ISME journal. https://www.nature.com/articles/ismej201637

However, more recent methods such as using balances and log ratio shall be considered in a larger, extended study. Test scores for PERMANOVA results shall be addressed in version 2 as will several questions that were unclear in version 1.

To date, we have none found literature that was comparable (ie. Having Bifidobacterium elevated in men, Prevotella elevated in women). Hence we attempted to suggest possible explanations to this observation. The discussion will be re-worded to take into account reviewer’s comments.
Competing Interests: No competing interest. Close
Report a concern
Reviewer Response 20 May 2019

Andres Gomez, Department of Animal Science, University of Minnesota, Saint Paul, USA

20 May 2019

Reviewer Response

Dear authors,

Could you please provide point by point responses to the comments raised in the previous version, along with where changes were made in the manuscript?
That would facilitate reviewing of ... Continue reading Dear authors,

Could you please provide point by point responses to the comments raised in the previous version, along with where changes were made in the manuscript?
That would facilitate reviewing of the new manuscript.
Dear authors,

Could you please provide point by point responses to the comments raised in the previous version, along with where changes were made in the manuscript?
That would facilitate reviewing of the new manuscript.
Competing Interests: No competing interests were disclosed. Close
Report a concern
Author Response 20 May 2019

Li-Fang Yeo, Tropical Medicine and Biology Platform, Monash University, Bandar Sunway, 46150, Malaysia

20 May 2019

Author Response

Dear reviewer,

As per request:

1. The QIIME microbiome analyses should be briefly described as far as the bioinformatics pipeline used and its details: clustering of OTUS? OTU picking? ... Continue reading Dear reviewer,

As per request:

1. The QIIME microbiome analyses should be briefly described as far as the bioinformatics pipeline used and its details: clustering of OTUS? OTU picking? amplicon sequence variants (DADA2-ASV)? These should be specifically detailed.

Microbiome analysis pipeline edited and detailed in Methods section.

2. It does not seem that the rarefaction curve shows compete plateau, also what do the colors mean?

Rarefaction curve now shows Shannon index instead of phylogenetic diversity for better indicator of plateau. Each colour is a biological sample.

3. It is unclear how kruskal-wallis test was used in conjunction with Unifrac - KW is used for pairwise comparisons in multivariate space, not for multivariate data.

KW test was conducted on an OTU table, not UniFrac. Refer to Ying & Sun (2017) https://www.sciencedirect.com/science/article/pii/S2352304217300351
Edits made in Results section.

4. PERMANOVA results (P values) should be accompanied by R2 and, Pseudo F values also. PERMANOVA models used should also be reported.

Pseudo F values are reported. R2 value not available for PERMANOVA but available for adonis. All available values are reported in supplementary material.

5. "Kruskal-Wallis test revealed that the relative abundance of Prevotella, Capnocytophaga, Leptotrichia, Neisseria and Streptococcus"

Were these discriminant analyses made at genus or OTU/ASV levels? Also, it is recommended to show these results in box plots showing medians and variation of data, not in tables

These analyses were made at OTU level 6, equivalent to genus level. Addressed in Results section. Also, due to time constraint and admittedly technical limitations on our part, the results will be shown in tables for this paper. Reviewer's advice are definitely well-received and will be considered for our next paper.

6. Authors make several assertions as to the character and typically reported nature of the markers discovered, in the context of "good/bad bacteria" - but many of the claims made are in the context of gut microbiomes not oral communities.

We have acknowledged the reviewer's comments and rephrased the discussion to be more cautious and less speculative.

7. "These commensal oral bacteria may become opportunistic pathogens in immune-compromised states." please cite

Idris A, Hasnain SZ, Huat LZ, Koh D. Human diseases, immunity and the oral microbiota—Insights gained from metagenomic studies. Oral Science International. 2017 Jul;14(2):27–32.

Cited in-text as well.

8. Fig 3- legends should indicate that a-b were based on Unw-UniFrac, and that C on W-UniFrac - Please depict % variation on axes

Done.

9. It is not clear how the post-hoc Dun's test were conducted-0 is this still based on PERMANOVA?

Post-hoc Dunn's test was conducted using an OTU table. Permanova does not conduct pairwise comparison hence Dunn's test was chosen to test for differences between people with different BMI status (underweight, normal, overweight). And the test revealed that underweight individuals had microbiomes that was significantly different from both normal and overweight individuals.

10. "The Kruskal-Wallis test generated by QIIME revealed no significant difference." This is not clear - compared to what?

12. QIIME is just a software with plugins to perform statistical analyses - it is not clear when authors are using qiime or R and how these platforms would yield different results. Please avoid jargon as "generated by QIIME"

For 10 & 12, previously qiime and R yielded different results (one not significant, one significantly different) for oral microbiomes and BMI. Upon closer inspection we realized this happened because of mismatching samples and OTU table upon input to R. We checked the other sets of data and the qiime results matches with R. Changes addressed in Results section.

11. "There was a significant difference between waist circumference and saliva microbiome composition (weighted UniFrac, p-value = 0.022). " This is not clear - you mean differences between the microbiomes of subjects with different waist circumference?

The above sentence has been rephrased to "There was a significant difference in the saliva microbiome of Temiar who had a healthy waist circumference compared to those with abdominal obesity."

13. "The saliva microbiome showed no significant difference between former smokers and never-smokers, suggesting perhaps the saliva microbiome may recover partially to an environment prior to smoking". This is unclear but also speculative.

The discussion is highly speculative and should be toned down, mainly as far as the hypotheses behind the differences reported.

Authors also seem to interchangeably extrapolate hypothesis on oral microbiome changes with gut microbiome changes in the context of lifestyle and diet (e.g. Prevotella) although diet may also be linked to oral microbiome modulation, this relationship is less clear, and taxa such as Prevotella, abundant in both gut and oral cavity, cannot be assumed to be modulated by diet in both sites

This comment has been well-received and edits have been done to Discussion section.

TQ.
Dear reviewer,

As per request:

1. The QIIME microbiome analyses should be briefly described as far as the bioinformatics pipeline used and its details: clustering of OTUS? OTU picking? amplicon sequence variants (DADA2-ASV)? These should be specifically detailed.

Microbiome analysis pipeline edited and detailed in Methods section.

2. It does not seem that the rarefaction curve shows compete plateau, also what do the colors mean?

Rarefaction curve now shows Shannon index instead of phylogenetic diversity for better indicator of plateau. Each colour is a biological sample.

3. It is unclear how kruskal-wallis test was used in conjunction with Unifrac - KW is used for pairwise comparisons in multivariate space, not for multivariate data.

KW test was conducted on an OTU table, not UniFrac. Refer to Ying & Sun (2017) https://www.sciencedirect.com/science/article/pii/S2352304217300351
Edits made in Results section.

4. PERMANOVA results (P values) should be accompanied by R2 and, Pseudo F values also. PERMANOVA models used should also be reported.

Pseudo F values are reported. R2 value not available for PERMANOVA but available for adonis. All available values are reported in supplementary material.

5. "Kruskal-Wallis test revealed that the relative abundance of Prevotella, Capnocytophaga, Leptotrichia, Neisseria and Streptococcus"

Were these discriminant analyses made at genus or OTU/ASV levels? Also, it is recommended to show these results in box plots showing medians and variation of data, not in tables

These analyses were made at OTU level 6, equivalent to genus level. Addressed in Results section. Also, due to time constraint and admittedly technical limitations on our part, the results will be shown in tables for this paper. Reviewer's advice are definitely well-received and will be considered for our next paper.

6. Authors make several assertions as to the character and typically reported nature of the markers discovered, in the context of "good/bad bacteria" - but many of the claims made are in the context of gut microbiomes not oral communities.

We have acknowledged the reviewer's comments and rephrased the discussion to be more cautious and less speculative.

7. "These commensal oral bacteria may become opportunistic pathogens in immune-compromised states." please cite

Idris A, Hasnain SZ, Huat LZ, Koh D. Human diseases, immunity and the oral microbiota—Insights gained from metagenomic studies. Oral Science International. 2017 Jul;14(2):27–32.

Cited in-text as well.

8. Fig 3- legends should indicate that a-b were based on Unw-UniFrac, and that C on W-UniFrac - Please depict % variation on axes

Done.

9. It is not clear how the post-hoc Dun's test were conducted-0 is this still based on PERMANOVA?

Post-hoc Dunn's test was conducted using an OTU table. Permanova does not conduct pairwise comparison hence Dunn's test was chosen to test for differences between people with different BMI status (underweight, normal, overweight). And the test revealed that underweight individuals had microbiomes that was significantly different from both normal and overweight individuals.

10. "The Kruskal-Wallis test generated by QIIME revealed no significant difference." This is not clear - compared to what?

12. QIIME is just a software with plugins to perform statistical analyses - it is not clear when authors are using qiime or R and how these platforms would yield different results. Please avoid jargon as "generated by QIIME"

For 10 & 12, previously qiime and R yielded different results (one not significant, one significantly different) for oral microbiomes and BMI. Upon closer inspection we realized this happened because of mismatching samples and OTU table upon input to R. We checked the other sets of data and the qiime results matches with R. Changes addressed in Results section.

11. "There was a significant difference between waist circumference and saliva microbiome composition (weighted UniFrac, p-value = 0.022). " This is not clear - you mean differences between the microbiomes of subjects with different waist circumference?

The above sentence has been rephrased to "There was a significant difference in the saliva microbiome of Temiar who had a healthy waist circumference compared to those with abdominal obesity."

13. "The saliva microbiome showed no significant difference between former smokers and never-smokers, suggesting perhaps the saliva microbiome may recover partially to an environment prior to smoking". This is unclear but also speculative.

The discussion is highly speculative and should be toned down, mainly as far as the hypotheses behind the differences reported.

Authors also seem to interchangeably extrapolate hypothesis on oral microbiome changes with gut microbiome changes in the context of lifestyle and diet (e.g. Prevotella) although diet may also be linked to oral microbiome modulation, this relationship is less clear, and taxa such as Prevotella, abundant in both gut and oral cavity, cannot be assumed to be modulated by diet in both sites

This comment has been well-received and edits have been done to Discussion section.

TQ.
Competing Interests: No competing interests. Close
Report a concern
Respond or Comment

COMMENTS ON THIS REPORT

Author Response 09 May 2019

Li-Fang Yeo, Tropical Medicine and Biology Platform, Monash University, Bandar Sunway, 46150, Malaysia

09 May 2019

Author Response

Details on the bioinformatics pipeline will be addressed in version 2. The different colours on the rarefaction curve represent different individuals. Other metrics such as Shannon and Chao1 shows an ... Continue reading Details on the bioinformatics pipeline will be addressed in version 2. The different colours on the rarefaction curve represent different individuals. Other metrics such as Shannon and Chao1 shows an obvious plateau. We felt that PD was more informative because it included phylogenetic distances.

Our statistics were done mostly with reference to this paper:
Wu, J., Peters, B.A., Dominianni, C., Zhang, Y., Pei, Z., Yang, L., Ma, Y., Purdue, M.P., Jacobs, E.J., Gapstur, S.M. and Li, H., 2016. “Cigarette smoking and the oral microbiome in a large study of American adults”. The ISME journal. https://www.nature.com/articles/ismej201637

However, more recent methods such as using balances and log ratio shall be considered in a larger, extended study. Test scores for PERMANOVA results shall be addressed in version 2 as will several questions that were unclear in version 1.

To date, we have none found literature that was comparable (ie. Having Bifidobacterium elevated in men, Prevotella elevated in women). Hence we attempted to suggest possible explanations to this observation. The discussion will be re-worded to take into account reviewer’s comments.

Details on the bioinformatics pipeline will be addressed in version 2. The different colours on the rarefaction curve represent different individuals. Other metrics such as Shannon and Chao1 shows an obvious plateau. We felt that PD was more informative because it included phylogenetic distances.

Our statistics were done mostly with reference to this paper:
Wu, J., Peters, B.A., Dominianni, C., Zhang, Y., Pei, Z., Yang, L., Ma, Y., Purdue, M.P., Jacobs, E.J., Gapstur, S.M. and Li, H., 2016. “Cigarette smoking and the oral microbiome in a large study of American adults”. The ISME journal. https://www.nature.com/articles/ismej201637

However, more recent methods such as using balances and log ratio shall be considered in a larger, extended study. Test scores for PERMANOVA results shall be addressed in version 2 as will several questions that were unclear in version 1.

To date, we have none found literature that was comparable (ie. Having Bifidobacterium elevated in men, Prevotella elevated in women). Hence we attempted to suggest possible explanations to this observation. The discussion will be re-worded to take into account reviewer’s comments.
Competing Interests: No competing interest. Close
Report a concern
Reviewer Response 20 May 2019

Andres Gomez, Department of Animal Science, University of Minnesota, Saint Paul, USA

20 May 2019

Reviewer Response

Dear authors,

Could you please provide point by point responses to the comments raised in the previous version, along with where changes were made in the manuscript?
That would facilitate reviewing of ... Continue reading Dear authors,

Could you please provide point by point responses to the comments raised in the previous version, along with where changes were made in the manuscript?
That would facilitate reviewing of the new manuscript.
Dear authors,

Could you please provide point by point responses to the comments raised in the previous version, along with where changes were made in the manuscript?
That would facilitate reviewing of the new manuscript.
Competing Interests: No competing interests were disclosed. Close
Report a concern
Author Response 20 May 2019

Li-Fang Yeo, Tropical Medicine and Biology Platform, Monash University, Bandar Sunway, 46150, Malaysia

20 May 2019

Author Response

Dear reviewer,

As per request:

1. The QIIME microbiome analyses should be briefly described as far as the bioinformatics pipeline used and its details: clustering of OTUS? OTU picking? ... Continue reading Dear reviewer,

As per request:

1. The QIIME microbiome analyses should be briefly described as far as the bioinformatics pipeline used and its details: clustering of OTUS? OTU picking? amplicon sequence variants (DADA2-ASV)? These should be specifically detailed.

Microbiome analysis pipeline edited and detailed in Methods section.

2. It does not seem that the rarefaction curve shows compete plateau, also what do the colors mean?

Rarefaction curve now shows Shannon index instead of phylogenetic diversity for better indicator of plateau. Each colour is a biological sample.

3. It is unclear how kruskal-wallis test was used in conjunction with Unifrac - KW is used for pairwise comparisons in multivariate space, not for multivariate data.

KW test was conducted on an OTU table, not UniFrac. Refer to Ying & Sun (2017) https://www.sciencedirect.com/science/article/pii/S2352304217300351
Edits made in Results section.

4. PERMANOVA results (P values) should be accompanied by R2 and, Pseudo F values also. PERMANOVA models used should also be reported.

Pseudo F values are reported. R2 value not available for PERMANOVA but available for adonis. All available values are reported in supplementary material.

5. "Kruskal-Wallis test revealed that the relative abundance of Prevotella, Capnocytophaga, Leptotrichia, Neisseria and Streptococcus"

Were these discriminant analyses made at genus or OTU/ASV levels? Also, it is recommended to show these results in box plots showing medians and variation of data, not in tables

These analyses were made at OTU level 6, equivalent to genus level. Addressed in Results section. Also, due to time constraint and admittedly technical limitations on our part, the results will be shown in tables for this paper. Reviewer's advice are definitely well-received and will be considered for our next paper.

6. Authors make several assertions as to the character and typically reported nature of the markers discovered, in the context of "good/bad bacteria" - but many of the claims made are in the context of gut microbiomes not oral communities.

We have acknowledged the reviewer's comments and rephrased the discussion to be more cautious and less speculative.

7. "These commensal oral bacteria may become opportunistic pathogens in immune-compromised states." please cite

Idris A, Hasnain SZ, Huat LZ, Koh D. Human diseases, immunity and the oral microbiota—Insights gained from metagenomic studies. Oral Science International. 2017 Jul;14(2):27–32.

Cited in-text as well.

8. Fig 3- legends should indicate that a-b were based on Unw-UniFrac, and that C on W-UniFrac - Please depict % variation on axes

Done.

9. It is not clear how the post-hoc Dun's test were conducted-0 is this still based on PERMANOVA?

Post-hoc Dunn's test was conducted using an OTU table. Permanova does not conduct pairwise comparison hence Dunn's test was chosen to test for differences between people with different BMI status (underweight, normal, overweight). And the test revealed that underweight individuals had microbiomes that was significantly different from both normal and overweight individuals.

10. "The Kruskal-Wallis test generated by QIIME revealed no significant difference." This is not clear - compared to what?

12. QIIME is just a software with plugins to perform statistical analyses - it is not clear when authors are using qiime or R and how these platforms would yield different results. Please avoid jargon as "generated by QIIME"

For 10 & 12, previously qiime and R yielded different results (one not significant, one significantly different) for oral microbiomes and BMI. Upon closer inspection we realized this happened because of mismatching samples and OTU table upon input to R. We checked the other sets of data and the qiime results matches with R. Changes addressed in Results section.

11. "There was a significant difference between waist circumference and saliva microbiome composition (weighted UniFrac, p-value = 0.022). " This is not clear - you mean differences between the microbiomes of subjects with different waist circumference?

The above sentence has been rephrased to "There was a significant difference in the saliva microbiome of Temiar who had a healthy waist circumference compared to those with abdominal obesity."

13. "The saliva microbiome showed no significant difference between former smokers and never-smokers, suggesting perhaps the saliva microbiome may recover partially to an environment prior to smoking". This is unclear but also speculative.

The discussion is highly speculative and should be toned down, mainly as far as the hypotheses behind the differences reported.

Authors also seem to interchangeably extrapolate hypothesis on oral microbiome changes with gut microbiome changes in the context of lifestyle and diet (e.g. Prevotella) although diet may also be linked to oral microbiome modulation, this relationship is less clear, and taxa such as Prevotella, abundant in both gut and oral cavity, cannot be assumed to be modulated by diet in both sites

This comment has been well-received and edits have been done to Discussion section.

TQ.
Dear reviewer,

As per request:

1. The QIIME microbiome analyses should be briefly described as far as the bioinformatics pipeline used and its details: clustering of OTUS? OTU picking? amplicon sequence variants (DADA2-ASV)? These should be specifically detailed.

Microbiome analysis pipeline edited and detailed in Methods section.

2. It does not seem that the rarefaction curve shows compete plateau, also what do the colors mean?

Rarefaction curve now shows Shannon index instead of phylogenetic diversity for better indicator of plateau. Each colour is a biological sample.

3. It is unclear how kruskal-wallis test was used in conjunction with Unifrac - KW is used for pairwise comparisons in multivariate space, not for multivariate data.

KW test was conducted on an OTU table, not UniFrac. Refer to Ying & Sun (2017) https://www.sciencedirect.com/science/article/pii/S2352304217300351
Edits made in Results section.

4. PERMANOVA results (P values) should be accompanied by R2 and, Pseudo F values also. PERMANOVA models used should also be reported.

Pseudo F values are reported. R2 value not available for PERMANOVA but available for adonis. All available values are reported in supplementary material.

5. "Kruskal-Wallis test revealed that the relative abundance of Prevotella, Capnocytophaga, Leptotrichia, Neisseria and Streptococcus"

Were these discriminant analyses made at genus or OTU/ASV levels? Also, it is recommended to show these results in box plots showing medians and variation of data, not in tables

These analyses were made at OTU level 6, equivalent to genus level. Addressed in Results section. Also, due to time constraint and admittedly technical limitations on our part, the results will be shown in tables for this paper. Reviewer's advice are definitely well-received and will be considered for our next paper.

6. Authors make several assertions as to the character and typically reported nature of the markers discovered, in the context of "good/bad bacteria" - but many of the claims made are in the context of gut microbiomes not oral communities.

We have acknowledged the reviewer's comments and rephrased the discussion to be more cautious and less speculative.

7. "These commensal oral bacteria may become opportunistic pathogens in immune-compromised states." please cite

Idris A, Hasnain SZ, Huat LZ, Koh D. Human diseases, immunity and the oral microbiota—Insights gained from metagenomic studies. Oral Science International. 2017 Jul;14(2):27–32.

Cited in-text as well.

8. Fig 3- legends should indicate that a-b were based on Unw-UniFrac, and that C on W-UniFrac - Please depict % variation on axes

Done.

9. It is not clear how the post-hoc Dun's test were conducted-0 is this still based on PERMANOVA?

Post-hoc Dunn's test was conducted using an OTU table. Permanova does not conduct pairwise comparison hence Dunn's test was chosen to test for differences between people with different BMI status (underweight, normal, overweight). And the test revealed that underweight individuals had microbiomes that was significantly different from both normal and overweight individuals.

10. "The Kruskal-Wallis test generated by QIIME revealed no significant difference." This is not clear - compared to what?

12. QIIME is just a software with plugins to perform statistical analyses - it is not clear when authors are using qiime or R and how these platforms would yield different results. Please avoid jargon as "generated by QIIME"

For 10 & 12, previously qiime and R yielded different results (one not significant, one significantly different) for oral microbiomes and BMI. Upon closer inspection we realized this happened because of mismatching samples and OTU table upon input to R. We checked the other sets of data and the qiime results matches with R. Changes addressed in Results section.

11. "There was a significant difference between waist circumference and saliva microbiome composition (weighted UniFrac, p-value = 0.022). " This is not clear - you mean differences between the microbiomes of subjects with different waist circumference?

The above sentence has been rephrased to "There was a significant difference in the saliva microbiome of Temiar who had a healthy waist circumference compared to those with abdominal obesity."

13. "The saliva microbiome showed no significant difference between former smokers and never-smokers, suggesting perhaps the saliva microbiome may recover partially to an environment prior to smoking". This is unclear but also speculative.

The discussion is highly speculative and should be toned down, mainly as far as the hypotheses behind the differences reported.

Authors also seem to interchangeably extrapolate hypothesis on oral microbiome changes with gut microbiome changes in the context of lifestyle and diet (e.g. Prevotella) although diet may also be linked to oral microbiome modulation, this relationship is less clear, and taxa such as Prevotella, abundant in both gut and oral cavity, cannot be assumed to be modulated by diet in both sites

This comment has been well-received and edits have been done to Discussion section.

TQ.
Competing Interests: No competing interests. Close
Report a concern

Comments on this article Comments (2)

Version 3

VERSION 3 PUBLISHED 28 May 2019

Revised

Comment

Version 1

VERSION 1 PUBLISHED 11 Feb 2019

Discussion is closed on this version, please comment on the latest version above.

Author Response 09 May 2019

Li-Fang Yeo, Tropical Medicine and Biology Platform, Monash University, Bandar Sunway, 46150, Malaysia

09 May 2019

Author Response

Thank you Andres and Siti Nursheena for your comments. We have addressed all your concerns, hopefully to your satisfaction, with one exception regarding Andres's suggestion to use boxplots to show ... Continue reading Thank you Andres and Siti Nursheena for your comments. We have addressed all your concerns, hopefully to your satisfaction, with one exception regarding Andres's suggestion to use boxplots to show median and variation of data. For this publication, we have decided to maintain the results in table-form showing q-values due to time constraint and technical limitations on our side. We will definitely take this into consideration to improve our future work.
Thank you Andres and Siti Nursheena for your comments. We have addressed all your concerns, hopefully to your satisfaction, with one exception regarding Andres's suggestion to use boxplots to show median and variation of data. For this publication, we have decided to maintain the results in table-form showing q-values due to time constraint and technical limitations on our side. We will definitely take this into consideration to improve our future work.
Competing Interests: No competing interests. Close
Report a concern
Author Response 24 Apr 2019

Maude Phipps, Jeffrey Cheah School of Medicine and Health Sciences, Monash University, Bandar Sunway, 46150, Malaysia

24 Apr 2019

Author Response

Thanks Andres and Siti Nursheena for your constructive comments. We endeavor to improve this paper and hope it will be confirmed and citable soon.
Competing Interests: No competing interests were disclosed.
Thanks Andres and Siti Nursheena for your constructive comments. We endeavor to improve this paper and hope it will be confirmed and citable soon.
Thanks Andres and Siti Nursheena for your constructive comments. We endeavor to improve this paper and hope it will be confirmed and citable soon.
Competing Interests: No competing interests were disclosed. Close
Report a concern
Discussion is closed on this version, please comment on the latest version above.

Open Peer Review

Reviewer Status

Reviewer Reports

	Invited Reviewers
	1	2
Version 3 (revision) 28 May 19	read
Version 2 (revision) 09 May 19	read	read
Version 1 11 Feb 19	read	read

Andres Gomez, University of Minnesota, Saint Paul, USA
Siti Nursheena Mohd Zain, University of Malaya, Kuala Lumpur, Malaysia

Comments on this article

All Comments(2)

Add a comment

Browse by related subjects

Back to all reports

Reviewer Report

13 Views

29 May 2019 | for Version 3

Andres Gomez, Department of Animal Science, University of Minnesota, Saint Paul, MN, USA

13 Views Cite this report Responses(1)

Approved

Thanks to the authors for addressing all the pending issues.

Competing Interests

No competing interests were disclosed.

I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard.

Respond to this report

Responses (1)

Author Response

31 May 2019

Li-Fang Yeo, Tropical Medicine and Biology Platform, Monash University, Bandar Sunway, 46150, Malaysia

The authors thank the reviewer very much for his time and all comments were very much appreciated.

View more View less

Competing Interests

No competing interest.

Back to all reports

Reviewer Report

9 Views

22 May 2019 | for Version 2

Andres Gomez, Department of Animal Science, University of Minnesota, Saint Paul, MN, USA

9 Views Cite this report Responses(0)

Approved With Reservations

Competing Interests

No competing interests were disclosed.

Respond to this report

Responses (0)

Back to all reports

Reviewer Report

7 Views

17 May 2019 | for Version 2

Siti Nursheena Mohd Zain, Institute of Biological Science, Faculty of Science, University of Malaya, Kuala Lumpur, Malaysia

7 Views Cite this report Responses(0)

Approved

The article reads well and all issues resolved.

Competing Interests

No competing interests were disclosed.

Reviewer Expertise

Parasitology and molecular parasitology

I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard.

Respond to this report

Responses (0)

Back to all reports

Reviewer Report

12 Views

10 Apr 2019 | for Version 1

Siti Nursheena Mohd Zain, Institute of Biological Science, Faculty of Science, University of Malaya, Kuala Lumpur, Malaysia

12 Views Cite this report Responses(0)

Approved With Reservations

Is the work clearly and accurately presented and does it cite the current literature?

Yes
Is the study design appropriate and is the work technically sound?

Yes
Are sufficient details of methods and analysis provided to allow replication by others?

Yes
If applicable, is the statistical analysis and its interpretation appropriate?

I cannot comment. A qualified statistician is required.
Are all the source data underlying the results available to ensure full reproducibility?

Yes
Are the conclusions drawn adequately supported by the results?

Yes

Competing Interests

No competing interests were disclosed.

Reviewer Expertise

Parasitology and molecular parasitology

Respond to this report

Responses (0)

Back to all reports

Reviewer Report

28 Views

13 Mar 2019 | for Version 1

Andres Gomez, Department of Animal Science, University of Minnesota, Saint Paul, MN, USA

28 Views Cite this report Responses(3)

Approved With Reservations

Is the work clearly and accurately presented and does it cite the current literature?

Partly
Is the study design appropriate and is the work technically sound?

Partly
Are sufficient details of methods and analysis provided to allow replication by others?

Partly
If applicable, is the statistical analysis and its interpretation appropriate?

Partly
Are all the source data underlying the results available to ensure full reproducibility?

Partly
Are the conclusions drawn adequately supported by the results?

Partly

Competing Interests

No competing interests were disclosed.

Reviewer Expertise

Microbial ecology

Respond to this report

Responses (3)

Author Response

09 May 2019

Li-Fang Yeo, Tropical Medicine and Biology Platform, Monash University, Bandar Sunway, 46150, Malaysia

Details on the bioinformatics pipeline will be addressed in version 2. The different colours on the rarefaction curve represent different individuals. Other metrics such as Shannon and Chao1 shows an obvious plateau. We felt that PD was more informative because it included phylogenetic distances.

Our statistics were done mostly with reference to this paper:
Wu, J., Peters, B.A., Dominianni, C., Zhang, Y., Pei, Z., Yang, L., Ma, Y., Purdue, M.P., Jacobs, E.J., Gapstur, S.M. and Li, H., 2016. “Cigarette smoking and the oral microbiome in a large study of American adults”. The ISME journal. https://www.nature.com/articles/ismej201637

However, more recent methods such as using balances and log ratio shall be considered in a larger, extended study. Test scores for PERMANOVA results shall be addressed in version 2 as will several questions that were unclear in version 1.

To date, we have none found literature that was comparable (ie. Having Bifidobacterium elevated in men, Prevotella elevated in women). Hence we attempted to suggest possible explanations to this observation. The discussion will be re-worded to take into account reviewer’s comments.

View more View less

Competing Interests

No competing interest.

Reviewer Response

20 May 2019

Andres Gomez, Department of Animal Science, University of Minnesota, Saint Paul, USA

Dear authors,

Could you please provide point by point responses to the comments raised in the previous version, along with where changes were made in the manuscript?
That would facilitate reviewing of the new manuscript.

View more View less

Competing Interests

No competing interests were disclosed.

Author Response

20 May 2019

Li-Fang Yeo, Tropical Medicine and Biology Platform, Monash University, Bandar Sunway, 46150, Malaysia

Dear reviewer,

As per request:

1. The QIIME microbiome analyses should be briefly described as far as the bioinformatics pipeline used and its details: clustering of OTUS? OTU picking? amplicon sequence variants (DADA2-ASV)? These should be specifically detailed.

Microbiome analysis pipeline edited and detailed in Methods section.

2. It does not seem that the rarefaction curve shows compete plateau, also what do the colors mean?

Rarefaction curve now shows Shannon index instead of phylogenetic diversity for better indicator of plateau. Each colour is a biological sample.

3. It is unclear how kruskal-wallis test was used in conjunction with Unifrac - KW is used for pairwise comparisons in multivariate space, not for multivariate data.

KW test was conducted on an OTU table, not UniFrac. Refer to Ying & Sun (2017) https://www.sciencedirect.com/science/article/pii/S2352304217300351
Edits made in Results section.

4. PERMANOVA results (P values) should be accompanied by R2 and, Pseudo F values also. PERMANOVA models used should also be reported.

Pseudo F values are reported. R2 value not available for PERMANOVA but available for adonis. All available values are reported in supplementary material.

5. "Kruskal-Wallis test revealed that the relative abundance of Prevotella, Capnocytophaga, Leptotrichia, Neisseria and Streptococcus"

Were these discriminant analyses made at genus or OTU/ASV levels? Also, it is recommended to show these results in box plots showing medians and variation of data, not in tables

These analyses were made at OTU level 6, equivalent to genus level. Addressed in Results section. Also, due to time constraint and admittedly technical limitations on our part, the results will be shown in tables for this paper. Reviewer's advice are definitely well-received and will be considered for our next paper.

6. Authors make several assertions as to the character and typically reported nature of the markers discovered, in the context of "good/bad bacteria" - but many of the claims made are in the context of gut microbiomes not oral communities.

We have acknowledged the reviewer's comments and rephrased the discussion to be more cautious and less speculative.

7. "These commensal oral bacteria may become opportunistic pathogens in immune-compromised states." please cite

Idris A, Hasnain SZ, Huat LZ, Koh D. Human diseases, immunity and the oral microbiota—Insights gained from metagenomic studies. Oral Science International. 2017 Jul;14(2):27–32.

Cited in-text as well.

8. Fig 3- legends should indicate that a-b were based on Unw-UniFrac, and that C on W-UniFrac - Please depict % variation on axes

Done.

9. It is not clear how the post-hoc Dun's test were conducted-0 is this still based on PERMANOVA?

Post-hoc Dunn's test was conducted using an OTU table. Permanova does not conduct pairwise comparison hence Dunn's test was chosen to test for differences between people with different BMI status (underweight, normal, overweight). And the test revealed that underweight individuals had microbiomes that was significantly different from both normal and overweight individuals.

10. "The Kruskal-Wallis test generated by QIIME revealed no significant difference." This is not clear - compared to what?

12. QIIME is just a software with plugins to perform statistical analyses - it is not clear when authors are using qiime or R and how these platforms would yield different results. Please avoid jargon as "generated by QIIME"

For 10 & 12, previously qiime and R yielded different results (one not significant, one significantly different) for oral microbiomes and BMI. Upon closer inspection we realized this happened because of mismatching samples and OTU table upon input to R. We checked the other sets of data and the qiime results matches with R. Changes addressed in Results section.

11. "There was a significant difference between waist circumference and saliva microbiome composition (weighted UniFrac, p-value = 0.022). " This is not clear - you mean differences between the microbiomes of subjects with different waist circumference?

The above sentence has been rephrased to "There was a significant difference in the saliva microbiome of Temiar who had a healthy waist circumference compared to those with abdominal obesity."

13. "The saliva microbiome showed no significant difference between former smokers and never-smokers, suggesting perhaps the saliva microbiome may recover partially to an environment prior to smoking". This is unclear but also speculative.

The discussion is highly speculative and should be toned down, mainly as far as the hypotheses behind the differences reported.

Authors also seem to interchangeably extrapolate hypothesis on oral microbiome changes with gut microbiome changes in the context of lifestyle and diet (e.g. Prevotella) although diet may also be linked to oral microbiome modulation, this relationship is less clear, and taxa such as Prevotella, abundant in both gut and oral cavity, cannot be assumed to be modulated by diet in both sites

This comment has been well-received and edits have been done to Discussion section.

TQ.

View more View less

Competing Interests

No competing interests.

Alongside their report, reviewers assign a status to the article:

Approved - the paper is scientifically sound in its current form and only minor, if any, improvements are suggested

Approved with reservations - A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit.

Not approved - fundamental flaws in the paper seriously undermine the findings and conclusions

[1] 1. Phipps ME, Chan KK, Naidu R, et al.: Cardio-metabolic health risks in indigenous populations of Southeast Asia and the influence of urbanization. BMC Public Health. 2015; 15: 47. PubMed Abstract | Publisher Full Text | Free Full Text

[2] 2. McColl H, Racimo F, Vinner L, et al.: The prehistoric peopling of Southeast Asia. Science. 2018; 361(6397): 88–92. PubMed Abstract | Publisher Full Text

[3] 3. Wong CY, Zalilah MS, Chua EY, et al.: Double-burden of malnutrition among the indigenous peoples (Orang Asli) of Peninsular Malaysia. BMC Public Health. 2015; 15: 680. PubMed Abstract | Publisher Full Text | Free Full Text

[4] 4. Tuan Abdul Aziz TA, Teh LK, Md Idris MH, et al.: Increased risks of cardiovascular diseases and insulin resistance among the Orang Asli in Peninsular Malaysia. BMC Public Health. 2016; 16: 284. PubMed Abstract | Publisher Full Text | Free Full Text

[5] 5. Olsen I: The oral microbiome in health and disease. Infections and General Health. 2016; 97–114. Publisher Full Text

[6] 6. Yamashita Y, Takeshita T: The oral microbiome and human health. J Oral Sci. 2017; 59(2): 201–206. PubMed Abstract | Publisher Full Text

[7] 7. Sudhakara P, Gupta A, Bhardwaj A, et al.: Oral Dysbiotic Communities and Their Implications in Systemic Diseases. Dent J (Basel). 2018; 6(2): pii: E10. PubMed Abstract | Publisher Full Text | Free Full Text

[8] 8. Idris A, Hasnain SZ, Huat LZ, et al.: Human diseases, immunity and the oral microbiota—Insights gained from metagenomic studies. Oral Science Int. 2017; 14(2): 27–32. Publisher Full Text

[9] 9. Long J, Cai Q, Steinwandel M, et al.: Association of oral microbiome with type 2 diabetes risk. J Periodontal Res. 2017; 52(3): 636–643. PubMed Abstract | Publisher Full Text | Free Full Text

[10] 10. Zeigler CC, Persson GR, Wondimu B, et al.: Microbiota in the oral subgingival biofilm is associated with obesity in adolescence. Obesity (Silver Spring). 2012; 20(1): 157–164. PubMed Abstract | Publisher Full Text

[11] 11. Choromańska K, Choromańska B, Dąbrowska E, et al.: Saliva of obese patients - is it different? Postepy Hig Med Dosw (Online). 2015; 69: 1190–1195. PubMed Abstract | Publisher Full Text

[12] 12. Yeo LF, Phipps M, Aghakhanian FF, et al.: Health and Saliva Microbiomes of a Semi-urbanized Indigenous Tribe in Peninsular Malaysia. figshare. Fileset. 2019. http://www.doi.org/10.26180/5c453ff435883

[13] 13. Dewhirst FE, Chen T, Izard J, et al.: The human oral microbiome. J Bacteriol. 2010; 192(19): 5002–5017. PubMed Abstract | Publisher Full Text | Free Full Text

[14] 14. Lozupone C, Knight R: UniFrac: a new phylogenetic method for comparing microbial communities. Appl Environ Microbiol. 2005; 71(12): 8228–8235. PubMed Abstract | Publisher Full Text | Free Full Text

[15] 15. Burns-Cox CJ: Splenomegaly and blood pressure in an Orang Asli community in West Malaysia. Am Heart J. 1970; 80(5): 718–719. PubMed Abstract | Publisher Full Text

[16] 16. Burns-Cox CJ, Chong YH, Gillman R: Risk factors and the absence of coronary heart disease in aborigines in West Malaysia. Br Heart J. 1972; 34(9): 953–8. PubMed Abstract | Publisher Full Text | Free Full Text

[17] 17. Pell C, Allotey P, Evans N, et al.: Coming of age, becoming obese: a cross-sectional analysis of obesity among adolescents and young adults in Malaysia. BMC Public Health. 2016; 16(1): 1082. PubMed Abstract | Publisher Full Text | Free Full Text

[18] 18. Burns J, Thomson N: Summary of overweight and obesity among Indigenous peoples. Retrieved [03 November 2015]; 2006.

[19] 19. Oliveira GF, Oliveira TR, Ikejiri AT, et al.: Prevalence of Obesity and Overweight in an Indigenous Population in Central Brazil: A Population-Based Cross-Sectional Study. Obes Facts. 2015; 8(5): 302–310. PubMed Abstract | Publisher Full Text | Free Full Text

[20] 20. Kandpal V, Sachdeva MP, Saraswathy KN: An assessment study of CVD related risk factors in a tribal population of India. BMC Public Health. 2016; 16(1): 434. PubMed Abstract | Publisher Full Text | Free Full Text

[21] 21. Institute of Public Health: National Health and Morbidity Survey. Institute of Public Health; 2015. Reference Source

[22] 22. Kumar PS: Sex and the subgingival microbiome: do female sex steroids affect periodontal bacteria? Periodontol 2000. 2013; 61(1): 103–124. PubMed Abstract | Publisher Full Text

[23] 23. Stahringer SS, Clemente JC, Corley RP, et al.: Nurture trumps nature in a longitudinal survey of salivary bacterial communities in twins from early adolescence to early adulthood. Genome Res. 2012; 22(11): 2146–2152. PubMed Abstract | Publisher Full Text | Free Full Text

[24] 24. Chen H, Jiang W: Application of high-throughput sequencing in understanding human oral microbiome related with health and disease. Front Microbiol. 2014; 5: 508. PubMed Abstract | Publisher Full Text | Free Full Text

[25] 25. Li J, Quinque D, Horz HP, et al.: Comparative analysis of the human saliva microbiome from different climate zones: Alaska, Germany, and Africa. BMC Microbiol. 2014; 14(1): 316. PubMed Abstract | Publisher Full Text | Free Full Text

[26] 26. Takeshita T, Kageyama S, Furuta M, et al.: Bacterial diversity in saliva and oral health-related conditions: the Hisayama Study. Sci Rep. 2016; 6(1): 22164. PubMed Abstract | Publisher Full Text | Free Full Text

[27] 27. Palmer RJ Jr: Composition and development of oral bacterial communities. Periodontol 2000. 2014; 64(1): 20–39. PubMed Abstract | Publisher Full Text | Free Full Text

[28] 28. Larsen JM: The immune response to Prevotella bacteria in chronic inflammatory disease. Immunology. 2017; 151(4): 363–374. PubMed Abstract | Publisher Full Text | Free Full Text

[29] 29. Benjamin G: Temiar Religion, 1964-2012: Enchantment, Disenchantment and Re-enchantment in Malaysia’s Uplands. 2014. Reference Source

[30] 30. Goodson JM, Hartman ML, Shi P, et al.: The salivary microbiome is altered in the presence of a high salivary glucose concentration. PLoS One. 2017; 12(3): e0170437. PubMed Abstract | Publisher Full Text | Free Full Text

[31] 31. Janem WF, Scannapieco FA, Sabharwal A, et al.: Salivary inflammatory markers and microbiome in normoglycemic lean and obese children compared to obese children with type 2 diabetes. PLoS One. 2017; 12(3): e0172647. PubMed Abstract | Publisher Full Text | Free Full Text

[32] 32. Wu J, Peters BA, Dominianni C, et al.: Cigarette smoking and the oral microbiome in a large study of American adults. ISME J. 2016; 10(10): 2435–2446. PubMed Abstract | Publisher Full Text | Free Full Text

[33] 33. Hernandez BY, Zhu X, Goodman MT, et al.: Betel nut chewing, oral premalignant lesions, and the oral microbiome. PLoS One. 2017; 12(2): e0172196. PubMed Abstract | Publisher Full Text | Free Full Text

[34] 34. Guha N, Warnakulasuriya S, Vlaanderen J, et al.: Betel quid chewing and the risk of oral and oropharyngeal cancers: a meta-analysis with implications for cancer control. Int J Cancer. 2014; 135(6): 1433–43. PubMed Abstract | Publisher Full Text

[35] 35. World Health Organization: Use of glycated haemoglobin (HbA1c) in the diagnosis of diabetes mellitus. Abbreviated report of a WHO consultation 2011. Geneva: World Health Organization. 2013. Reference Source

[36] 36. Mostafa SA, Davies MJ, Morris DH, et al.: The association of the triglyceride-to-HDL cholesterol ratio with insulin resistance in White European and South Asian men and women. PLoS One. 2012; 7(12): e50931. PubMed Abstract | Publisher Full Text | Free Full Text

[37] 37. Quinque D, Kittler R, Kayser M, et al.: Evaluation of saliva as a source of human DNA for population and association studies. Anal Biochem. 2006; 353(2): 272–277. PubMed Abstract | Publisher Full Text

[38] 38. PCR I: 16S Metagenomics Sequencing Library Preparation. Illumina. 2013. Reference Source

[39] 39. Caporaso JG, Kuczynski J, Stombaugh J, et al.: QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010; 7(5): 335–336. PubMed Abstract | Publisher Full Text | Free Full Text

[40] 40. Edgar RC, Haas BJ, Clemente JC, et al.: UCHIME improves sensitivity and speed of chimera detection. Bioinformatics. 2011; 27(16): 2194–2200. PubMed Abstract | Publisher Full Text | Free Full Text

[41] 41. DeSantis TZ, Hugenholtz P, Larsen N, et al.: Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl Environ Microbiol. 2006; 72(7): 5069–5072. PubMed Abstract | Publisher Full Text | Free Full Text

[42] 42. Edgar RC: Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010; 26(19): 2460–2461. PubMed Abstract | Publisher Full Text

[43] 43. Wang Q, Garrity GM, Tiedje JM, et al.: Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol. 2007; 73(16): 5261–5267. PubMed Abstract | Publisher Full Text | Free Full Text

Health and saliva microbiomes of a semi-urbanized indigenous tribe in Peninsular Malaysia

Abstract

Keywords

Revised Amendments from Version 2

Introduction

Results

Anthropometrics and cardio-metabolic health

Table 1. Anthropometrics and cardio-metabolic risk factors among Temiar.

Saliva microbiome analysis

Figure 1. Rarefaction curve plotted using alpha diversity metric, Shannon index against number of sequences per sample. Each coloured line represents one biological sample.

Figure 2. Relative abundance of bacteria found in the oral microbiome of Temiar at the phylum level.

Figure 3. PCoA plots with the larger shapes representing the mean UniFrac distance between samples segregated by (a) gender – unweighted UniFrac; (b) body mass index (BMI) – unweighted UniFrac; (c) smoking status – weighted UniFrac.

Table 2. Relative abundance of bacteria (genus level) differently represented in Temiar using Kruskal-Wallis test.

Table 3. Relative abundance of bacteria found among smokers vs former and never-smokers using Kruskal-Wallis test.

Discussion

Conclusion

Methods

Data and sample collection

Anthropometrics and Biochemical analysis

DNA extraction and PCR

Sequencing on Illumina MiSeq

Data analysis

Statistical analysis

Data availability

Underlying data

Extended data

Grant information

Acknowledgements

References

Comments on this article Comments (2)

Open Peer Review

Comments on this article Comments (2)

Open Peer Review

Reviewer Status

Reviewer Reports

Comments on this article

Browse by related subjects

Competing Interests Policy

Stay Updated