Next Article in Journal
Identification and Functional Analysis of NLP-Encoding Genes from the Postharvest Pathogen Penicillium expansum
Previous Article in Journal
The Revival of the Battle between David and Goliath in the Enteric Viruses and Microbiota Struggle: Potential Implication for Celiac Disease
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Microorganisms
Volume 7
Issue 6
10.3390/microorganisms7060174
microorganisms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Influence of Modernization and Disease on the Gastric Microbiome of Orang Asli, Myanmars and Modern Malaysians

by
Eng Guan Chua
1,*,
Mun Fai Loke
2,3,†,
Selva Perumal Gunaletchumy
2,†,
Han Ming Gan
4,5,
Kavitha Thevakumar
2,
Chin Yen Tay
1,
Sylvia Young
6,
Than Than Aye
7,
Win Win Maw
8,
Mya Mya Aye
9,
Alex Hwong-Ruey Leow
10,
Ahmad Najib Azmi
10,
Sri Ganesh Kalimuthu
11,
Haji Mohd Akmal Bin Dato Dahlan
12,
Khean Lee Goh
10 and
Jamuna Vadivelu
2,†
1
The Marshall Centre for Infectious Diseases Research and Training, University of Western Australia, Perth, WA 6009, Australia
2
Department of Medical Microbiology, University of Malaya, Kuala Lumpur 50603, Malaysia
3
School of Life Sciences and Chemical Technology, Ngee Ann Polytechnic, Singapore 599489, Singapore
4
School of Science, Monash University Malaysia, Selangor 47500, Malaysia
5
Centre for Integrative Ecology, School of Life and Environmental Sciences, Deakin University, Melbourne, VIC 3216, Australia
6
Harry Perkins Institute of Medical Research, The Centre for Medical Research at the University of Western Australia, Perth, WA 6009, Australia
7
Department of Gastroenterology, Thingangyun Sanpya General Hospital, Yangon 11071, Myanmar
8
Microbiology Department, University of Medicine 2, Yangon 11031, Myanmar
9
University of Medical Technology, Mandalay 05071, Myanmar
10
Department of Medicine, University of Malaya, Kuala Lumpur 50603, Malaysia
11
CliniPath Malaysia Sdn Bhd, Kuala Lumpur 55100, Malaysia
12
Hospital Orang Asli Gombak, Selangor 53100, Malaysia
 
add Show full affiliation list
 
remove Hide full affiliation list
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Microorganisms 2019, 7(6), 174; https://doi.org/10.3390/microorganisms7060174
Submission received: 6 May 2019 / Revised: 5 June 2019 / Accepted: 12 June 2019 / Published: 14 June 2019
Browse Figures
Versions Notes

Abstract

:
The present study explored the differences in gastric microbiome between three distinct populations of Southeast Asia. These include the isolated Orang Asli population and modern Malaysians, as well as patients from Myanmar, the least developed country in the region. All 79 subjects recruited in this study had Helicobacter pylori infection. Based on alpha diversity analysis, Orang Asli had the richest and most diverse gastric microbiome, followed by Myanmar and modern Malaysian groups. Beta diversity analysis revealed significant separation of samples between different populations. These observations are likely to be associated with the level of modernization of each population. Our data further suggested increased bacterial species richness and diversity of the gastric microbiome in individuals who were less modernized, particularly in the Orang Asli group, could suppress the growth of H. pylori. In addition, there were significant variations in the gastric microbiome between modern Malaysians with different types of gastric diseases. Notably, Cutibacterium acnes was present at significantly greater abundance level in patients with non-ulcerative dyspepsia than those with peptic-ulcer diagnosis. This suggests that C. acnes may also play a role in gastritis besides H. pylori, which merits further investigation.

1. Introduction

The human body does not contain only our own cells but is also a host to numerous microorganisms. To illustrate this scenario, Joshua Lederberg coined the term “microbiome” to signify “the ecological community of commensal, symbiotic, and pathogenic microorganisms that literally share our body space” [1]. Up to 60% of human-associated microbes could not be cultured in vitro [2,3,4], leading to substantial underestimation of the actual microbiota diversity in different parts of human body. For years, the study of the gastric microbial composition was stagnated as the human stomach was considered a sterile organ due to its hostile acidic environment [5]. However, the discovery of Helicobacter pylori in 1982 relinquished this central dogma [6].
A recent study conducted on two populations in Colombia found striking differences in microbiome composition between the two populations, which could be related to dietary differences and gastric cancer risk [7]. This supports the need for population-based microbiome studies to better understand how changes related to lifestyle, diet, socioeconomic level and disease status could alter the microbiota, and if possible, to revert or mimic such changes for better human health. To date, while there have been studies comparing the compositions of gut microbiota in populations of different ethnic and geographic origins [8,9,10,11], such studies are generally lacking in Southeast Asia.
Malaysia and Myanmar are both multiracial and multicultural countries located in Southeast Asia. In Malaysia, the major populations are Malay, Chinese and Indian, while Myanmar’s primary residents consist of Burman, Shan, Karen, Rakhine, Chinese, Indian and Mon. Historically, Myanmar and Malaysia were located along the Maritime Silk Road and were once colonized by the British until 1948 and 1957, respectively. Since then, while Malaysia continues to advance, Myanmar civilization, under the isolationistic policies of its military government from 1962 to 2011, had remained stagnant [12]. The level of modernization in these two major populations can be reflected by their access to clean drinking water and sanitary facilities, which are both vital development indicators. Until 2015, only 68% and 64.7% of households in the less-developed Myanmar nation had access to clean drinking water and basic sanitation service [13]. In contrast, 96% and 99.6% of Malaysian households were provided with clean water supply and adequate sanitation, respectively.
Orang Asli are the indigenous minorities residing in Peninsular Malaysia. More importantly, they are the oldest inhabitants of Peninsular Malaysia. The Orang Asli population is a collective of 18 ethnic subgroups further classified under Negrito, Senoi and Proto-Malay for administrative purposes [14]. While the Negrito comprises the Kensiu, Kintaq, Jahai, Lanoh, Mendriq and Bateq subethnic groups, classified under the Senoi are the Semai, Temiar, Jah Hut, Semaq Beri, Che Wong and Mah Meri subgroups. Proto-Malay, on the other hand, contains the Temuan, Temoq, Semelai, Jakun, Orang Kuala, Orang Kanaq and Orang Seletar subpopulations [15]. Orang Asli mostly live in remote mountainous areas, away from urbanization with no access to clean water supply and proper sanitary service [16]. Each subgroup can be found within a specific geographical space usually next to a river valley and isolated from the others.
This study highlighted the differences in gastric microbiome of three distinct populations from the Southeast Asia region (the isolated Orang Asli of Malaysia, Myanmar residents and modern Malaysians), which is likely associated with the level of modernization. We further demonstrated the influence of gastric diseases on human gastric microbiome, albeit interestingly, only in modern Malaysian patients.

2. Materials and Methods

2.1. Ethics Statement

This study was approved by the National Institutes of Health of Malaysia Clinical Research Centre and Medical Research & Ethics Committee (NMRR-12-671-11909), University of Malaya Medical Centre (UMMC) Medical Ethics Committee (MEC Ref. No. 1023.3) and the Ethics and Research Committee of the University of Medicine 2, Myanmar. Gastric biopsy samples were obtained with informed and written consent from subjects who presented for endoscopy at UMMC and Thingangyun Sanpya General Hospital.

2.2. Screening of H. pylori Infections in Orang Asli

Between 28 October and 12 December 2013, 177 asymptomatic volunteers from five different locations: (1) Kampung Sungai Padi, Lipis district, Pahang state; (2) Kampung Donglai, Hulu Langat district, Selangor state; (3) Cameron Highlands district, Pahang state; (4) Gua Musang district, Kelantan state and; (5) Gombak district and Selangor state, were tested for H. pylori infection. Among these 177 volunteers; 69, 45, 4 and 1 were from the Temiar, Semai, Mah Meri and Jah Hut subethnic groups classified under the Senoi group, respectively; 56 and 1 were from the Temuan and Jakun subgroups derived from the Proto-Malay group, respectively; and a remaining volunteer who was of Iban origin that does not belong to any of the three primary Orang Asli tribal groups. The test was conducted using the PYtest 14C-urea breath test (UBT) kit (Tri-Med, West Leederville, WA, Australia). Breath samples were analyzed using the microCOUNT Lite Liquid Scintillation Counter. Results were reported as disintegrations per minute (DPM). A DPM count of less than 50 was interpreted as a negative result while DPM counts in the range of 50–199 were classified as indeterminate. A DPM count of 200 and above was considered positive for H. pylori infection.

2.3. Collection of Gastric Biopsies

Of the 79 gastric tissue biopsy samples collected for this study, from September 2011 to June 2014, 10 and 32 were from UBT-positive; Orang Asli volunteers and non-Orang Asli Malaysian patients who were symptomatic and previously tested positive for H. pylori infection, respectively. These individuals underwent endoscopic examination followed by the collection of gastric biopsies at UMMC’s Endoscopy Unit. The remaining 37 samples were collected from Myanmar patients at the Department of Gastroenterology of Thingangyun Sanpya General Hospital, Yangon, Myanmar. All subjects were further diagnosed with gastric cancer (GC), peptic ulcer disease (PUD) or non-ulcerative dyspepsia (NUD) based on endoscopy evaluation.
Biopsy specimens were obtained from the gastric antrum and fundus for metagenomics purposes, CLOtest rapid urease testing, H. pylori culturing, and histopathological examination. Culturing of H. pylori isolate was not performed on Myanmar patient samples in this study. Histopathological examination was performed independently by Dr. Sri Ganesh Kalimuthu of the CliniPath Malaysia Sdn Bhd.

2.4. Illumina 16S rRNA Gene Sequencing

Genomic DNA was extracted from gastric biopsy samples using MasterPureTM DNA purification kit (Epicentre, Madison, WI, USA) according to manufacturer’s instructions. The V3–V4 region of 16S rRNA gene was amplified using S-D-Bact-0341-b-S-17 and S-D-Bact-0785-a-A-21 primers designed to include the Illumina-compatible adaptors [17,18]. The 16S amplicon libraries were prepared according to Illumina 16S library preparation protocol [19]. Initial PCR amplification was performed using NEBNext High-Fidelity Master Mix (New England Biolabs, Ipswich, MA, USA) with the following conditions: An initial denaturation at 98 °C for 30 s, followed by 30 cycles consisting of denaturation (98 °C for 10 s), annealing (60 °C for 2 min) and extension (72 °C for 20 s), and a final extension step at 72 °C for 1 min. Automated cluster generation and a 2 × 250 bp paired-end sequencing (MiSeq 500-cycle reagent kit V2) was carried out on the MiSeq platform (Illumina, San Diego, CA, USA) at the Monash University Malaysia Genomics Facility. Metagenomics data reported in this paper are available for public access through the MG-RAST server. The accession number and demographic characteristics of each sample are listed in Table S1.

2.5. Data Analysis

The 16S sequencing data were quality-trimmed using Sickle (https://github.com/najoshi/sickle) using the following parameters: –q 20 –l 200. Merging of overlapping paired-end sequences was performed using MeFit with default parameters [20]. Filtering of chimeric sequences, de novo greedy clustering of 16S rRNA sequences into Operational Taxonomic Units (OTUs) at 97% similarity threshold, and removal of singleton and chimeric OTUs were conducted using Micca (version 1.7.0) [21]. Taxonomic assignment of the representative OTUs was performed using the Bayesian LCA-based taxonomic classification method with a 1−100 cut-off e-value and 100 bootstrap replications, against NCBI 16S microbial database [22]. Taxonomic assignment at each level was accepted with a minimum confidence score of 80. The OTU table including taxonomic information is available in Table S2. Unclassified OTUs were retained and included in all following analyses.
Microbial diversity analysis was performed on rarefied OTU abundance matrix (depth value of 17358) using QIIME (version 1.9.1) [23]. Alpha diversity was evaluated based on the following metrics: Observed species and Shannon diversity indices. A non-parametric two-sample t-test with Bonferroni correction was used to compare the alpha diversity metrics between different groups. Principle coordinate analysis (PcoA) using Bray–Curtis and Jaccard distance metrics was performed to visualize separation of samples. Non-parametric statistical analysis of the distance metric was performed using ANOSIM with 1000 permutations. Plots were generated using PhyloToAST software [24]. Prediction of KEGG functional pathways was conducted using Piphillin and the results are available in Table S3 [25,26]. Pearson correlations between KEGG pathways and bacterial genera were examined using microbiomeSeq R package (https://github.com/umerijaz/microbiomeSeq.git). Only negative and positive correlations of ≤−0.8 and ≥0.8, respectively, and with p-values less than 0.01, were considered.

2.6. Statistical Analysis

Bacterial genera and KEGG functional pathways with minimum average relative abundances of 0.1% were compared between groups by Wilcoxon rank-sum testing, followed by Bonferroni correction. Differences with adjusted p-values of less than or equal to 0.01 were considered significant.

3. Results

3.1. H. pylori Prevalence among Orang Asli

Among 177 Orang Asli subjects who were screened for H. pylori infection, only 18 (10.2%) and 53 (30.3%) were tested positive and negative, respectively. The rest, interestingly, had indeterminate UBT readings. The demographics and UBT results of the Orang Asli subjects in this study are summarized in Table 1. Ten UBT-positive Orang Asli individuals consented for further endoscopy examination and provided gastric biopsies for this study.

3.2. The Gastric Microbial Community of H. pylori-Positive Orang Asli, Myanmar and Modern Malaysian Subjects

The general characteristics and clinical diagnosis of study subjects who underwent endoscopic examination are summarized in Table 2. A total of 22,526,214 reads were generated for all 79 samples. Following quality trimming and merging of overlapping paired-end sequences, 9,912,950 sequences with an average length of 454.3 ± 8.1 bp were retained. The number of sequences per sample ranged from 20,042 to 233,069 (median 115,971). Of 2233 OTUs acquired via de novo clustering algorithm, 1466 and 1105 were taxonomically assigned down to and phylum and genus level, respectively, with an 80% confidence threshold, which were subsequently aggregated to 20 bacterial phyla and 443 bacterial genera, respectively, for further analysis.
As our study subjects differed not only by population but also by disease status, a comparison of those with only NUD was therefore performed to exclude the influence of PUD and GC when examining the microbial differences between three major populations including Orang Asli, modern Malaysians comprising Chinese, Indian and Malay ethnic groups, and Myanmar subjects. Overall, Proteobacteria, Firmicute and Actinobacteria constituted the three most predominant phyla, at 71.9%, 10.4% and 9.1% of relative abundances (Figure 1A), respectively, in the gastric microbiome of modern Malaysian patients. Similarly, in the Orang Asli samples, while having both Proteobacteria and Actinobacteria as the first and third most abundant phyla, at 57.3% and 5.6%, respectively, Bacteroidetes, not Firmicute, was the second most common bacterial phylum (13.4%) identified. In Myanmar patients, the top three most abundant phyla were Firmicutes, Proteobacteria and Bacteroidetes, at 39%, 28.1% and 15.5%, respectively.
At the genus level, Helicobacter was the most abundant bacterial genus identified in the modern Malaysians, with a relative abundance of 61% (Figure 1B). In both Myanmar and Orang Asli patients, however, the Helicobacter bacterial loads were only 9.8% and 0.25%, respectively. While a high abundance of Lactobacillus (7.2%) was detected in Myanmar patients, Orang Asli group contained mostly environmental bacterial genera including Flavobacterium, Sulfuritalea, Desulfatiglans, Brevibacillus and Massilia.
We next analyzed the diversity of bacterial populations amongst Orang Asli, Myanmar and modern Malaysian samples. Alpha diversity based on species richness and Shannon diversity indices was estimated. Both modern Malaysian and Myanmar groups had significantly reduced microbial richness (Figure 2A) and diversity (Figure 2B) relative to Orang Asli. It is worth mentioning that the gastric microbial community of modern Malaysians was the least diverse among all. To assess bacterial community differences between groups, principal coordinate analysis (PCoA) based on Jaccard metric was performed. The analysis revealed significant dissimilarities in microbial composition among the groups (Figure 3, R = 0.75, p < 0.001).
To identify distinguishing bacterial genera between Orang Asli, Myanmar and modern Malaysian NUD samples, 76 genera with minimum average abundances of 0.1% were subjected to Wilcoxon rank-sum testing followed by Bonferroni correction. Of 28 significant genera identified, as summarized in Table 3, 17 were profoundly enriched in the Orang Asli gastric microbiome compared to that of Myanmar and modern Malaysian samples. While Nocardioides was nearly absent in the Myanmar samples, modern Malaysians had the highest Helicobacter bacterial load, significantly greater than that of Orang Asli and Myanmar cohorts.
Functional prediction using Piphillin revealed 306 KEGG pathways (Table S3). Following rarefaction and filtering non-NUD samples, we then compared between different populations the pathways with total average abundances of at least 0.1% and considered only those with a minimum relative fold change of 1.25 and a Bonferroni-adjusted p-value of less than or equal to 0.01 as statistically significant. Seventy-six pathways were significant (Table S4). Notably more than one third of pathways (28/76) were significantly up- or down-regulated in Orang Asli samples in comparison to both Myanmar and modern Malaysian cohorts. Of 28 significantly differential pathways identified in Orang Asli, half were involved in metabolism of amino acids and organic compounds like ascorbate and aldarate, histidine, tryptophan, inositol phosphate, arginine and proline, phenylalanine, tyrosine and beta-alanine, which had positive significant association with Bosea, Brevundimonas, Emticicia, Sphingorhabdus, Tardiphaga and Variovorax; six environmental bacterial genera that were found present almost exclusively in the Orang Asli population (Table S5). Compared to both Orang Asli and Myanmar cohorts, significant up-regulation of FoxO signaling was observed in the modern Malaysian samples and was shown to be positively and negatively associated with Helicobacter and Streptococcus, respectively.

3.3. Gastric Diseases Altered the Microbial Community Structure in Modern Malaysian Patients

We attempted additional diversity analysis to determine whether different gastric disease outcomes would affect the structural composition of gastric microbial community. The analysis was performed separately for each population to eliminate the confounding effect of population factor on disease status. Orang Asli group was excluded as NUD was the only gastric disorder diagnosed in this population. While clear separations between NUD, PUD and GC samples were observed in modern Malaysian samples (Figure 4, Jaccard index, R = 0.4, p = 0.002), interestingly, no evident grouping of Myanmar samples according to disease types could be observed.
Both NUD and PUD samples shared similar major bacterial phyla compositions with Proteobacteria and Firmicutes being the two most predominant phyla, at 72.1% and 10.1% in the former and 85.6% and 7.4% in the latter (Figure 5). The GC samples, however, had comparable abundances of both Proteobacteria and Firmicutes, at 33.2% and 36.3%, respectively. At genus level, while Helicobacter bacteria constituted 60.5% and 74.4% of gastric microbiome in both NUD and PUD groups, respectively, in GC group the total abundance of top five most predominant genera including Helicobacter was 37.8%, indicating that GC samples had a more diverse and complex gastric bacterial population.
Due to small GC sample size (n = 4), comparison of bacterial genera abundances by Wilcoxon ranked-sum testing was performed only between NUD and PUD groups. Following Bonferroni correction, of 55 genera tested, only one appeared significant, which is Cutibacterium (adjusted p = 0.005, median abundances of 0.159% and 0.002% for NUD and PUD, respectively).

4. Discussion

In this study, the gastric microbial compositions of three distinct populations, Orang Asli, Myanmars and modern Malaysians, who all had H. pylori infection based on positive UBT result, were investigated. It is important to note that the median abundance of Helicobacter bacteria in Orang Asli was merely 0.25%, which is several orders of magnitude less than that of modern Malaysian and Myanmar individuals. This is likely due to the significantly more diverse bacterial makeup in the gastric environment of Orang Asli compared to that of modern Malaysians and Myanmars. Although the human stomach serves as the primary reservoir for H. pylori colonization, the findings above indicate that the degree of colonization can be affected by the changes in the gastric microenvironment. It is evident that the presence of environmental bacterial species in the stomach of Orang Asli, as well as in the Myanmar individuals to a certain extent but less, could suppress the growth of H. pylori, while the lack of these foreign environmental isolates in modern Malaysians allows H. pylori to flourish. This finding could possibly explain why most Orang Asli individuals had a borderline UBT outcome, presuming most of them are likely to be H. pylori-positive.
In both Orang Asli and modern Malaysian groups, Proteobacteria was found to be the predominant bacterial phylum. While Proteobacteria in modern Malaysians contained primarily Helicobacter bacteria belonging to the class Epsilon Proteobacteria, Orang Asli harbored mostly environmental microbes of Alpha-, Beta-, Delta- and Gamma-proteobacteria classes. Notably, the Orang Asli gastric microbiome had a Firmicutes/Bacteroidetes ratio of 0.41, which is at least 6-fold less than that of Myanmar and modern Malaysian groups, at 2.51 and 2.65, respectively. It was previously shown that a high-fat/high-sugar diet promotes the growth of Firmicutes while reducing the abundance of Bacteroidetes in the distal gut of rodents [27]. Additionally, when an obese individual shifted to a fat-restricted diet or a carbohydrate-restricted low-calorie diet, the abundance levels of Bacteroidetes and Firmicutes were increased and decreased, respectively [28]. Together, these observations suggest that both modern Malaysian and Myanmar populations are likely to consume a diet high in carbohydrate and/or fat content possibly attributed to influx of Western culture as part of modernization process, while Orang Asli are still maintaining traditional plant-based diets low in fat and carbohydrate corresponding to their hunting/planting lifestyle [29].
Of the several bacterial genera that were significantly enriched in Orang Asli, Bosea, Sphingorhabdus, and Variovorax, they also displayed significant positive correlations to arginine and proline metabolism, and β-alanine metabolism, respectively. While a diet supplemented with arginine and proline was shown to improve the wound healing process in diabetic rats [30], β-alanine supplementation as reported in a recent meta-analysis was able to improve exercise performance [31]. It is thought that these bacteria may facilitate the metabolism and uptake of these amino acids that are widely available in vegetables and grains consumed by Orang Asli, which is then likely to increase the wound healing and physical capacities in Orang Asli to enhance their survival fitness in an outback setting.
While bacteria of the genus Nocardioides was found in nearly all modern Malaysian (95%, 19/20) and Orang Asli (100%, 10/10) samples, it was detected in only 23.1% (3/13) of Myanmar samples. Members of the genus Nocardioides are mostly environmental bacteria found in soil or plants and only one Nocardioides species has thus far been isolated from a human faecal sample [32]. The high prevalence of Nocardioides species could be a distinct feature of gastric microbiome in peninsular Malaysians including Orang Asli, albeit the influence of this Nocardioides species on human health is uncertain.
Contrary to modern Malaysian samples, the Myanmar population did not show distinct clustering of samples by gastric disease conditions. It is thought that lifestyle and dietary habits, geographical and socioeconomic factors, and/or host genetics may have greater influence on the gastric microbiome than disease conditions in the Myanmar population. In modern Malaysians, PUD samples harbored the most H. pylori, which is expected as this bacterium is an important driving factor for peptic ulcer disease outcome [33]. Notably, NUD samples had significantly more Cutibacterium bacteria than in the PUD samples. At the species level, it was further identified as C. acnes, formerly known as Propionibacterium acnes [34]. Interestingly, while H. pylori has been known for its dominant role in causing different gastric disorders, C. acnes was identified as a possible causative agent for lymphocytic gastritis [35]. It is important to clarify that NUD, also known as functional dyspepsia, is the preferred clinical diagnosis term assigned to patients with gastritis, on the basis of endoscopic and histological results [36]. To examine the potential role of C. acnes in the development of gastric diseases alongside H. pylori, further investigation is required. This would help to expand our current knowledge of bacterium-driven gastric diseases.

5. Conclusions

This study highlighted the taxonomic differences between the gastric microbiome of three distinct populations from Southeast Asia including Orang Asli, Myanmars and modern Malaysia, and such differences are possibly due to different levels of modernization. Distinguishing gastric microbial community structure could also be seen between different gastric disease conditions in modern Malaysians, albeit not in Myanmar patients. Importantly, the gastric microbiome of Orang Asli, who lived in isolation from urbanization, and a near primitive lifestyle with no use of any modern antibiotics and few modern medicines, may represent one of the most “original” snapshots of human gastric microbiome to-date.

Supplementary Materials

The following are available online at https://www.mdpi.com/2076-2607/7/6/174/s1, Table S1, Demographics data and accession numbers; Table S2, OTU table with taxonomic information; Table S3, KEGG pathway abundances predicted based on the gastric microbiome of Orang Asli, Myanmar and modern Malaysian groups; Table S4, Distinguishing KEGG pathways between Orang Asli, Myanmar and modern Malaysians groups; Table S5, Pearson correlation analysis between KEGG pathways and bacterial genera of Orang Asli, Myanmar and modern Malaysian groups.

Author Contributions

Conceptualization, M.F.L., H.M.A.B.D.D. and J.V.; formal analysis, E.G.C. and M.F.L.; funding acquisition, M.F.L., C.Y.T. and J.V.; investigation, S.P.G., H.M.G., K.T., T.T.A., W.W.M., M.M.A., A.H.-R.L., A.N.A., S.G.K. and K.L.G.; validation, E.G.C. and S.Y.; writing—original draft, S.P.G.; writing—review & editing, E.G.C. and M.F.L.

Funding

This study was supported by the University of Malaya Ministry of Education under the High Impact Research (HIR) Grant UM.C/625/1/HIR/MoE/CHAN/13/3 (No. H-50001-A000030) awarded to MFL and JV, and Shenzhen’s Sanming Project (Grant No. SZSM201510050) awarded to Chin Yen Tay.

Acknowledgments

We sincerely thank the volunteers and patients from Malaysia and Myanmar for their participation in this study. We would also like to thank everyone from the Endoscopy Unit at the University of Malaya Medical Center (Malaysia), Orang Asli Hospital Gombak (Malaysia) and Department of Gastroenterology at the Thingangyun Sanpya General Hospital (Myanmar) for their assistance. In addition, we would like to extend our appreciation to Prof. Tint Swe Latt (former Rector of the University of Medicine 2, Yangon) for his support of this collaborative project.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

References

  1. Lederberg, J.; McCray, A.T. ‘Ome Sweet ‘Omics-A Genealogical Treasury of Words. Scientist 2001, 15, 8. [Google Scholar]
  2. Pei, Z.; Bini, E.J.; Yang, L.; Zhou, M.; Francois, F.; Blaser, M.J. Bacterial biota in the human distal esophagus. Proc. Natl. Acad. Sci. USA 2004, 101, 4250–4255. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Aas, J.A.; Paster, B.J.; Stokes, L.N.; Olsen, I.; Dewhirst, F.E. Defining the normal bacterial flora of the oral cavity. J. Clin. Microbiol. 2005, 43, 5721–5732. [Google Scholar] [CrossRef] [PubMed]
  4. Verhelst, R.; Verstraelen, H.; Claeys, G.; Verschraegen, G.; Delanghe, J.; Van Simaey, L.; De Ganck, C.; Temmerman, M.; Vaneechoutte, M. Cloning of 16S rRNA genes amplified from normal and disturbed vaginal microflora suggests a strong association between Atopobium vaginae, Gardnerella vaginalis and bacterial vaginosis. BMC Microbiol. 2004, 4, 16. [Google Scholar] [CrossRef] [PubMed]
  5. Nardone, G.; Compare, D. The human gastric microbiota: Is it time to rethink the pathogenesis of stomach diseases? United Eur. Gastroenterol. J. 2015, 3, 255–260. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Marshall, B.J.; Warren, J.R. Unidentified curved bacilli in the stomach of patients with gastritis and peptic ulceration. Lancet 1984, 1, 1311–1315. [Google Scholar] [CrossRef]
  7. Yang, I.; Woltemate, S.; Piazuelo, M.B.; Bravo, L.E.; Yepez, M.C.; Romero-Gallo, J.; Delgado, A.G.; Wilson, K.T.; Peek, R.M.; Correa, P.; et al. Different gastric microbiota compositions in two human populations with high and low gastric cancer risk in Colombia. Sci. Rep. 2016, 6, 18594. [Google Scholar] [CrossRef]
  8. Yatsunenko, T.; Rey, F.E.; Manary, M.J.; Trehan, I.; Dominguez-Bello, M.G.; Contreras, M.; Magris, M.; Hidalgo, G.; Baldassano, R.N.; Anokhin, A.P.; et al. Human gut microbiome viewed across age and geography. Nature 2012, 486, 222–227. [Google Scholar] [CrossRef]
  9. De Filippo, C.; Cavalieri, D.; Di Paola, M.; Ramazzotti, M.; Poullet, J.B.; Massart, S.; Collini, S.; Pieraccini, G.; Lionetti, P. Impact of diet in shaping gut microbiota revealed by a comparative study in children from Europe and rural Africa. Proc. Natl. Acad. Sci. USA 2010, 107, 14691–14696. [Google Scholar] [CrossRef] [Green Version]
  10. Schnorr, S.L.; Candela, M.; Rampelli, S.; Centanni, M.; Consolandi, C.; Basaglia, G.; Turroni, S.; Biagi, E.; Peano, C.; Severgnini, M.; et al. Gut microbiome of the Hadza hunter-gatherers. Nat. Commun. 2014, 5, 3654. [Google Scholar] [CrossRef]
  11. Suzuki, T.A.; Worobey, M. Geographical variation of human gut microbial composition. Biol. Lett. 2014, 10, 20131037. [Google Scholar] [CrossRef] [PubMed]
  12. The World Factbook. Available online: https://www.cia.gov/library/publications/the-world-factbook/ (accessed on 31 January 2019).
  13. UNICEF. Drinking Water, Sanitation and Hygiene Database; The WHO/UNICEF: Geneva, Switzerland, 2017. [Google Scholar]
  14. Nicholas, C.; International Work Group for Indigenous Affairs (IWGIA); Center for Orang Asli Concerns (COAC). The Orang Asli and the Contest for Resources: Indigenous Politics, Development, and Identity in Peninsular Malaysia; International Work Group for Indigenous Affairs: Copenhagen, Denmark; Center for Orang Asli Concerns: Subang Jaya, Malaysia, 2000. [Google Scholar]
  15. Ching, C.C.S.; Ross, V. Preliminary Exploration of a Semai Musician’s Transmission of Indigenous Musical Traditions in Peninsular Malaysia. Pertanika J. Soc. Sci. Humanit. 2015, 23, 20. [Google Scholar]
  16. Masron, T.; Masami, F.; Ismail, N. Orang Asli in Peninsular Malaysia: Population, Spatial Distribution and Socio-Economic Condition. J. Ritsumeikan Soc. Sci. Hum. 2013, 6, 75–115. [Google Scholar]
  17. Klindworth, A.; Pruesse, E.; Schweer, T.; Peplies, J.; Quast, C.; Horn, M.; Glockner, F.O. Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies. Nucleic Acids Res. 2013, 41, e1. [Google Scholar] [CrossRef] [PubMed]
  18. Bartram, A.K.; Lynch, M.D.; Stearns, J.C.; Moreno-Hagelsieb, G.; Neufeld, J.D. Generation of multimillion-sequence 16S rRNA gene libraries from complex microbial communities by assembling paired-end illumina reads. Appl. Environ. Microbiol. 2011, 77, 3846–3852. [Google Scholar] [CrossRef] [PubMed]
  19. 16S Metagenomic Sequencing Library Preparation. Available online: https://support.illumina.com/documents/documentation/chemistry_documentation/16s/16s-metagenomic-library-prep-guide-15044223-b.pdf (accessed on 31 January 2019).
  20. Parikh, H.I.; Koparde, V.N.; Bradley, S.P.; Buck, G.A.; Sheth, N.U. MeFiT: Merging and filtering tool for illumina paired-end reads for 16S rRNA amplicon sequencing. BMC Bioinform. 2016, 17, 491. [Google Scholar] [CrossRef] [PubMed]
  21. Albanese, D.; Fontana, P.; De Filippo, C.; Cavalieri, D.; Donati, C. MICCA: A complete and accurate software for taxonomic profiling of metagenomic data. Sci Rep. 2015, 5, 9743. [Google Scholar] [CrossRef] [PubMed]
  22. Gao, X.; Lin, H.; Revanna, K.; Dong, Q. A Bayesian taxonomic classification method for 16S rRNA gene sequences with improved species-level accuracy. BMC Bioinform. 2017, 18, 247. [Google Scholar] [CrossRef]
  23. Caporaso, J.G.; Kuczynski, J.; Stombaugh, J.; Bittinger, K.; Bushman, F.D.; Costello, E.K.; Fierer, N.; Pena, A.G.; Goodrich, J.K.; Gordon, J.I.; et al. QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 2010, 7, 335–336. [Google Scholar] [CrossRef] [Green Version]
  24. Dabdoub, S.M.; Fellows, M.L.; Paropkari, A.D.; Mason, M.R.; Huja, S.S.; Tsigarida, A.A.; Kumar, P.S. PhyloToAST: Bioinformatics tools for species-level analysis and visualization of complex microbial datasets. Sci Rep. 2016, 6, 29123. [Google Scholar] [CrossRef] [Green Version]
  25. Iwai, S.; Weinmaier, T.; Schmidt, B.L.; Albertson, D.G.; Poloso, N.J.; Dabbagh, K.; DeSantis, T.Z. Piphillin: Improved Prediction of Metagenomic Content by Direct Inference from Human Microbiomes. PLoS ONE 2016, 11, e0166104. [Google Scholar] [CrossRef] [PubMed]
  26. Kanehisa, M.; Goto, S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000, 28, 27–30. [Google Scholar] [CrossRef] [PubMed]
  27. Turnbaugh, P.J.; Backhed, F.; Fulton, L.; Gordon, J.I. Diet-induced obesity is linked to marked but reversible alterations in the mouse distal gut microbiome. Cell Host Microbe 2008, 3, 213–223. [Google Scholar] [CrossRef] [PubMed]
  28. Ley, R.E.; Turnbaugh, P.J.; Klein, S.; Gordon, J.I. Microbial ecology: Human gut microbes associated with obesity. Nature 2006, 444, 1022–1023. [Google Scholar] [CrossRef] [PubMed]
  29. Mirnalini, K., Jr.; Zalilah, M.S.; Safiah, M.Y.; Tahir, A.; Siti Haslinda, M.D.; Siti Rohana, D.; Khairul Zarina, M.Y.; Mohd Hasyami, S.; Normah, H. Energy and Nutrient Intakes: Findings from the Malaysian Adult Nutrition Survey (MANS). Malays. J. Nutr. 2008, 14, 1–24. [Google Scholar] [PubMed]
  30. Raynaud-Simon, A.; Belabed, L.; Le Naour, G.; Marc, J.; Capron, F.; Cynober, L.; Darquy, S. Arginine plus proline supplementation elicits metabolic adaptation that favors wound healing in diabetic rats. Am. J. Physiol. Regul. Integr. Comp. Physiol. 2012, 303, R1053–R1061. [Google Scholar] [CrossRef]
  31. Hobson, R.M.; Saunders, B.; Ball, G.; Harris, R.C.; Sale, C. Effects of beta-alanine supplementation on exercise performance: A meta-analysis. Amino Acids 2012, 43, 25–37. [Google Scholar] [CrossRef]
  32. Dubourg, G.; Sankar, S.A.; Rathored, J.; Lagier, J.C.; Robert, C.; Couderc, C.; Papazian, L.; Raoult, D.; Fournier, P.E. Noncontiguous finished genome sequence and description of Nocardioides massiliensis sp. nov. GD13(T). New Microbes New Infect. 2016, 10, 47–57. [Google Scholar] [CrossRef]
  33. Covacci, A.; Telford, J.L.; Del Giudice, G.; Parsonnet, J.; Rappuoli, R. Helicobacter pylori virulence and genetic geography. Science 1999, 284, 1328–1333. [Google Scholar] [CrossRef]
  34. Scholz, C.F.; Kilian, M. The natural history of cutaneous propionibacteria, and reclassification of selected species within the genus Propionibacterium to the proposed novel genera Acidipropionibacterium gen. nov., Cutibacterium gen. nov. and Pseudopropionibacterium gen. nov. Int J. Syst. Evol. Microbiol. 2016, 66, 4422–4432. [Google Scholar] [CrossRef]
  35. Montalban-Arques, A.; Wurm, P.; Trajanoski, S.; Schauer, S.; Kienesberger, S.; Halwachs, B.; Hogenauer, C.; Langner, C.; Gorkiewicz, G. Propionibacterium acnes overabundance and natural killer group 2 member D system activation in corpus-dominant lymphocytic gastritis. J. Pathol. 2016, 240, 425–436. [Google Scholar] [CrossRef] [PubMed]
  36. Madisch, A.; Andresen, V.; Enck, P.; Labenz, J.; Frieling, T.; Schemann, M. The Diagnosis and Treatment of Functional Dyspepsia. Dtsch Arztebl. Int. 2018, 115, 222–232. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Microorganisms 07 00174 g001 550
Figure 1. Bacterial phyla and genera compositions in Orang Asli, Myanmar and modern Malaysian non-ulcerative dyspepsia (NUD) patients. (A) Mean relative abundances of the five most predominant bacterial phyla in each population. (B) Mean relative abundances of the five most predominant bacterial genera in each cohort.
Figure 1. Bacterial phyla and genera compositions in Orang Asli, Myanmar and modern Malaysian non-ulcerative dyspepsia (NUD) patients. (A) Mean relative abundances of the five most predominant bacterial phyla in each population. (B) Mean relative abundances of the five most predominant bacterial genera in each cohort.
Microorganisms 07 00174 g001
Microorganisms 07 00174 g002 550
Figure 2. Alpha diversity analysis between Orang Asli, Myanmar and modern Malaysian groups diagnosed with NUD. (A) Alpha diversity estimation based on species richness index. (B) Alpha diversity estimation based on Shannon diversity index.
Figure 2. Alpha diversity analysis between Orang Asli, Myanmar and modern Malaysian groups diagnosed with NUD. (A) Alpha diversity estimation based on species richness index. (B) Alpha diversity estimation based on Shannon diversity index.
Microorganisms 07 00174 g002
Microorganisms 07 00174 g003 550
Figure 3. Beta diversity analysis between Orang Asli, Myanmar and modern Malaysian groups diagnosed with NUD. Principal coordinate analysis (PCoA) plot based on Jaccard index (R = 0.75, p < 0.001).
Figure 3. Beta diversity analysis between Orang Asli, Myanmar and modern Malaysian groups diagnosed with NUD. Principal coordinate analysis (PCoA) plot based on Jaccard index (R = 0.75, p < 0.001).
Microorganisms 07 00174 g003
Microorganisms 07 00174 g004 550
Figure 4. Beta diversity analysis between modern Malaysians with gastric cancer, peptic ulcer disease and non-ulcerative dyspepsia. PCoA plot based on Jaccard index (R = 0.4, p = 0.002).
Figure 4. Beta diversity analysis between modern Malaysians with gastric cancer, peptic ulcer disease and non-ulcerative dyspepsia. PCoA plot based on Jaccard index (R = 0.4, p = 0.002).
Microorganisms 07 00174 g004
Microorganisms 07 00174 g005 550
Figure 5. Bacterial phyla and genera compositions in modern Malaysians with gastric cancer, peptic ulcer disease and non-ulcerative dyspepsia. (A) Mean relative abundances of the five most dominating bacterial phyla in each disease group. (B) Mean relative abundances of the five most dominating bacterial genera in each disease group.
Figure 5. Bacterial phyla and genera compositions in modern Malaysians with gastric cancer, peptic ulcer disease and non-ulcerative dyspepsia. (A) Mean relative abundances of the five most dominating bacterial phyla in each disease group. (B) Mean relative abundances of the five most dominating bacterial genera in each disease group.
Microorganisms 07 00174 g005
Table
Table 1. Urea breath test results and demographics of Orang Asli subjects.
Table 1. Urea breath test results and demographics of Orang Asli subjects.
OriginNPositiveIndeterminateNegativeMean Age (years)Married (%)Female (%)
Temiar693491733 (16–80)85%67%
Temuan563322129 (15–56)69%51%
Semai4511191532 (14–60)83%77%
Others (Mah Meri, Iban, Jah Hut and Jakun)716039 (27–53)43%57%
Table
Table 2. Demographics and clinical findings of study subjects who underwent endoscopic examination.
Table 2. Demographics and clinical findings of study subjects who underwent endoscopic examination.
CountryEthnic GroupNMean Age (years)NUDPUDGC
MalaysiaOrang Asli (Semai)333.0300
Orang Asli (Temiar)732.3700
Chinese1647.9763
Indian1034.6811
Malay637.5510
Overall4239.63084
MyanmarBurmese2248.35116
Rakhine547.4212
Karen336.7201
Chin241.5200
Mon231.5200
Others (Kachin, Paoh)361.7021
Overall3747.1131410
GC: Gastric cancer; NUD: Non-ulcerative dyspepsia; PUD: Peptic ulcer disease.
Table
Table 3. List of bacterial genera abundances that significantly differed between Orang Asli, Myanmar and modern Malaysian groups diagnosed with NUD.
Table 3. List of bacterial genera abundances that significantly differed between Orang Asli, Myanmar and modern Malaysian groups diagnosed with NUD.
GenusMedian Relative Abundance (%)Adjusted p-Values
OA
(n = 10)		MYA
(n = 13)		MM
(n = 20)		OA vs. MYAOA vs. MMMYA vs. MM
Aquabacterium0.69000.005<0.001N.S
Asticcacaulis0.77000.004<0.001N.S
Bosea1.24000.004<0.001N.S
Brevibacillus2.09000.004<0.001N.S
Brevundimonas1.860.0100.0040.001N.S
Caulobacter0.72000.0040.001N.S
Emticicia3.11000.004<0.001N.S
Flavobacterium7.7000.004<0.001N.S
Ideonella2.07000.004<0.001N.S
Legionella0.84000.004<0.001N.S
Massilia3.100.050.004<0.001N.S
Methyloversatilis0.63000.0040.001N.S
Phenylobacterium2.96000.004<0.001N.S
Sphingorhabdus2.18000.004<0.001N.S
Sulfuritalea3.96000.004<0.001N.S
Tardiphaga0.55000.004<0.001N.S
Variovorax1.49000.004<0.001N.S
Helicobacter0.010.4971.69N.S0.0010.01
Nocardioides0.8800.420.005N.S0.001
Agrobacterium0.62000.005N.SN.S
Streptococcus0.3611.061.890.007N.SN.S
Anaerococcus0.010.010.59N.S0.004N.S
Corynebacterium0.10.411.6N.S0.001N.S
Dehalogenimonas0.2900N.S0.01N.S
Desulfatiglans1.8500N.S<0.001N.S
Paludibaculum0.7100N.S0.006N.S
Rothia00.440.33N.S0.006N.S
Lactobacillus00.590N.SN.S0.004
OA: Orang Asli; MYA: Myanmars; MM: Modern Malaysians; N.S: Not significant.

Share and Cite

MDPI and ACS Style

Chua, E.G.; Loke, M.F.; Gunaletchumy, S.P.; Gan, H.M.; Thevakumar, K.; Tay, C.Y.; Young, S.; Aye, T.T.; Maw, W.W.; Aye, M.M.; et al. The Influence of Modernization and Disease on the Gastric Microbiome of Orang Asli, Myanmars and Modern Malaysians. Microorganisms 2019, 7, 174. https://doi.org/10.3390/microorganisms7060174

AMA Style

Chua EG, Loke MF, Gunaletchumy SP, Gan HM, Thevakumar K, Tay CY, Young S, Aye TT, Maw WW, Aye MM, et al. The Influence of Modernization and Disease on the Gastric Microbiome of Orang Asli, Myanmars and Modern Malaysians. Microorganisms. 2019; 7(6):174. https://doi.org/10.3390/microorganisms7060174

Chicago/Turabian Style

Chua, Eng Guan, Mun Fai Loke, Selva Perumal Gunaletchumy, Han Ming Gan, Kavitha Thevakumar, Chin Yen Tay, Sylvia Young, Than Than Aye, Win Win Maw, Mya Mya Aye, and et al. 2019. "The Influence of Modernization and Disease on the Gastric Microbiome of Orang Asli, Myanmars and Modern Malaysians" Microorganisms 7, no. 6: 174. https://doi.org/10.3390/microorganisms7060174

APA Style

Chua, E. G., Loke, M. F., Gunaletchumy, S. P., Gan, H. M., Thevakumar, K., Tay, C. Y., Young, S., Aye, T. T., Maw, W. W., Aye, M. M., Leow, A. H. -R., Azmi, A. N., Kalimuthu, S. G., Dahlan, H. M. A. B. D., Goh, K. L., & Vadivelu, J. (2019). The Influence of Modernization and Disease on the Gastric Microbiome of Orang Asli, Myanmars and Modern Malaysians. Microorganisms, 7(6), 174. https://doi.org/10.3390/microorganisms7060174

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop