Skip to main content

An exploration of the gut and environmental resistome in a community in northern Vietnam in relation to antibiotic use



Antibiotic resistance is a major global public health threat. Antibiotic use can directly impact the antibiotic resistant genes (ARGs) profile of the human intestinal microbiome and consequently the environment through shedding.


We determined the resistome of human feces, animal stools, human food and environmental (rain, well, and irrigative water) samples (n = 304) in 40 households within a community cohort and related the data to antibiotic consumption. Metagenomic DNA was isolated and qPCR was used to determine presence of mobile colistin resistance (mcr) genes, genes encoding extended-spectrum β-lactamases (ESBL), carbapenemases and quinolone resistance genes.


Nearly 40 % (39.5%, 120/304) of samples contained ESBL genes (most frequent were CTX-M-9 (23.7% [72/304]), CTX-M-1 (18.8% [57/304]). Quinolone resistance genes (qnrS) were detected in all human and 91% (41/45) of animal stool samples. Mcr-1 and mcr-3 were predominantly detected in human feces at 88% (82/93) and 55% (51/93) and animal feces at 93% (42/45) and 51% (23/45), respectively. Mcr-2, mrc-4 and mcr-5 were not detected in human feces, and only sporadically (< 6%) in other samples. Carbapenemase-encoding genes were most common in water (15% [14/91]) and cooked food (13% [10/75]) samples, while their prevalence in human and animal stools was lower at 4% in both human (4/93) and animal (2/45) samples. We did not find an association between recent antibiotic consumption and ARGs in human stools. Principal component analysis showed that the resistome differs between ecosystems with a strong separation of ARGs profiles of human and animal stools on the one hand versus cooked food and water samples on the other.


Our study indicated that ARGs were abundant in human and animal stools in a rural Vietnamese community, including ARGs targeting last resort antibiotics. The resistomes of animal and human stools were similar as opposed to the resistomes from water and food sources. No association between antibiotic use and ARG profiles was found in a setting of high background rates of AMR.


Antimicrobial resistance (AMR) has become omnipresent and is still increasing globally, making treatment of infectious diseases more and more challenging [1]. Current research includes assessing and analyzing the resistome: the combined collection of antibiotic resistance conferring genes and their precursors in both pathogenic and non-pathogenic bacteria [2].

AMR is recognized as a public health issue by the Vietnamese Ministry of Health [3]. Despite the introduction of a prescription law in Vietnam in 2005, 88 to 91% of all antibiotics were still sold without prescription in rural and urban private pharmacies in 2014 [4]. Use of antibiotics without prescription is usually inappropriate and considered an important driver of AMR [5]. In Vietnam, resistance to first-line drugs of treatment is high among common bacterial pathogens in both hospital and community settings [4, 6].

Besides human antibiotic consumption, antibiotics are widely used in livestock breeding and aquaculture for disease treatment, prevention and growth promotion [7]. Vietnam’s agricultural sector contributes 30% to the national Gross Domestic Product (GDP) and Vietnam is the 4th largest aquaculture producer globally [8]. Nearly half (43.7%) of commercial feed products contain at least one antibiotic. 66.7% of industrial-scale farms and 20% of household farms used antibiotics as growth promoters for pigs [9]. Since 2017 the use of antibiotics for growth promotion is illegal but the effect of this ban has yet to be determined [10]. This situation is largely responsible for the high proportion of antibiotic-resistant bacteria in food products and aquatic environments [7, 9, 11].

Antibiotic use has an important impact on the intestinal microbiota composition of humans and animals as well as on their surrounding environment. Recent studies considered the gut microbiome as an important antimicrobial resistance gene (ARG) reservoir in which ARGs between native and transient intestinal bacteria can be transferred [12, 13]. In a recent study, a total of 1093 ARGs were identified in 162 fecal metagenomes from human individuals from China, Denmark and Spain [14]. Another study on healthy individuals from seven countries detected ARGs for 50 of 68 classes of antibiotics in 252 fecal metagenomes, at an average of 21 ARGs per sample [15]. Exposure to antibiotic-resistant bacteria or genes present in animal stools, food products and water in combination with antibiotic use may lead to acquisition of ARGs in the human intestinal microbiome [16]. Here, we present longitudinally collected human antibiotic usage data and single time point ARG profiles of human, animal, food and environmental samples in a rural Vietnamese community.


Study population & design

From November 2014 to June 2015, we conducted a longitudinal study examining the impact of antibiotic use on the microbiome and resistome of humans, their food, animals and the environment in Vietnam. The population of this study was selected from an existing prospective, household-based community cohort established to investigate influenza virus transmission in the commune of Thanh Ha in Ha Nam province, Vietnam in 2007 [17]. From this community, we selected 80 households: 40 households with children under 11 years of age and 40 households without children under 11 years of age. At baseline, a semi-structured questionnaire was used to record demographic characteristics, social economic factors, living conditions, health status. Current antibiotic use data was collected weekly. In case of antibiotic use, we collected the information of drug name, dose, supplier, and indication for drug use.

Household members were included if they were willing to participate and were available for follow up. Household members were excluded if they were unable to collect fecal specimens, were known to be on antibiotic maintenance therapy or had a chronic intestinal condition.

Participants received a clean container with spoon attached to the lid (VWR corporate Headquaters, PA, USA), gloves, and a biohazard bag and study staff demonstrated the for home collection of the stool sample [18]. Processed food and water samples were collected following the sampling protocol for low-biomass samples as described elsewhere [19]. Fresh stool, processed food and water samples with the collection date and time were collected by study staff and sent in cold condition (4–8 °C) to the National Institute of Hygiene and Epidemiology (NIHE) and Oxford University Clinical Research Unit for sample processing within 4 h. Samples were stored at − 80 °C until they were used for microbiome and resistome analysis in Maastricht, Netherlands.

For the purpose of analyzing the resistome in this study, a subset of 40 out of the 80 households was selected: 20 households in which at least one member used any beta-lactam antibiotic and 20 households in which no antibiotic use was reported prior to sampling. At 4 months, we collected domestic animal feces, cooked food products and water samples from these households (Fig. 1). Three hundred four samples were collected including 93 human fecal samples, 45 animal fecal samples, 75 cooked food samples and 91 water samples (Fig. 1).

Fig. 1
figure 1

Sample population of study and sample selected of ARG explorative study

The research was approved by the Oxford University Tropical Research Ethics Committee (OxTREC, 49–14) and the National Institute of Hygiene and Epidemiology, Vietnam (NIHE) institutional review board. Written informed consent was obtained from all participants in the study.

DNA extraction

Repeated Bead-Beating (RBB) combined with column-based purification was used to extract DNA from human and animal fecal samples according to protocol Q (IHMS_SOP 06 V2 - of the International Human Microbiome Standards consortium [20]. Bead-beating was done using the FastPrep™ Instrument (MP Biomedicals, Santa Ana (CA), USA) with 0.1 mm zirconium-silica beads (BioSpec Products, Bartlesville (OK), USA) to homogenize feces. DNA was finally purified by adapting to QIAamp DNA Stool Mini kit columns (Qiagen, Hilden, Germany).

Regarding the isolation of metagenomic DNA from processed food, 200 mg of sample was homogenized in 0.75 ml PBS (pH 7.2), and centrifuged for 3 min at 13,000 x g. Metagenomic DNA was subsequently isolated by the QIAGEN DNeasy Power Food Kit according to the manufacturer’s instructions for DNA isolation from solid food. For isolation of microbial DNA from water samples, 100 ml of collected water was filtered through 0.22 μm mixed cellulose esters membrane filters (Sartorius, Göttingen, Germany) to capture bacteria. One quarter of the filters were used for metagenomic DNA extraction using the QIAGEN DNeasy Power Water kit according to the manufacturer’s protocol.

Upon isolation, DNA concentrations were determined using the Quan-iT PicoGreen dsDNA assay (Invitrogen, Carlsbad (CA), USA).

Quantitative detection of ARG identification and quantification by qPCR

qPCR amplification was performed to quantify the total number of 16S ribosomal RNA gene copies (16S rDNA) using MyiQ single – color real time PCR detection system (Bio-Rad, Hercules (CA), USA) as previously described [21] (Additonal file 1: Supplement 1).

The selected samples were also assessed for the presence and abundance of ARGs by quantitative real-time PCRs using a 7900HT Fast Real-time PCR System (Applied Biosystem Inc., Foster City (CA), USA) (Additonal file 1: Supplement 2).

The qPCRs to detect and to quantify the ESBL, carbapenemases and quinolone resistance genes were performed using primers and probes from previous studies [22,23,24,25] and subsequently validated protocols for their application on metagenomics DNA [26,27,28] (Additonal file 1: Supplement 1).

In order to identify and quantify plasmid-mediated colistin resistance (mcr) genes (mcr-1, mcr-2, mcr-3, mcr-4, mcr-5), we used in-house designed procedures, specifically for metagenomics DNA (Additional file 1: Supplement 1), from the Department of Medical Microbiology at the Maastricht University Medical Center. For all mcr-genes, products from positive samples were sequenced by using the PCR primers and an ABI BigDye Termination v1.1 Cycle Sequencing Kit (Applied Biosystems, Foster City, CA, USA). Sequencing data were analyzed by using BLAST ( Link).

Quantification of the abundance of AMR genes in samples was achieved by comparing the cycle threshold (Ct) values to a standard curve. The standard curses were constructed by subjecting serial 10-fold dilutions of positive plasmid constructs containing the AMR target sequences to the same qPCR assays.

The load of AMR genes in samples was determined using standard curves generated by serial dilutions of plasmid constructs containing the AMR gene of interest. Subsequently, we calculated the ratio between AMR gene copies and 16S rRNA gene copies to adjust for differences in bacterial load between samples. Samples that were negative for a particular ARG, were assigned a ARG copy number of zero.

Statistical analysis

Pearson’s chi-square test, or Fisher’s exact test in case of expected cell numbers < 5, was used when comparing the prevalence of ARGs or other dichotomous data between two groups. The statistical significance level was set at 5%.

To examine the association between the relative abundance of resistance genes in fecal samples and human antibiotic consumption, the ratio of ARGs to 16S rRNA gene copy numbers was calculated. The differential abundance of ARGs among the two antibiotic exposure groups was compared by Wilcoxon test (IBM SPSS Statistics v23.0, IBM, Armonk (NY), USA).

We applied principal component analysis (PCA) from the stats package and used ggbiplot package in R 3.5.3 to visualize the ordination of the ARG profiles generated from human and animal fecal samples, cooked food products and water. Before running PCA, we standardized the data to equally scale all variables and take into account zeros (samples without a specific AMR gene) in the data set.


Demographics of the complete cohort

At baseline, demographic information from 359 individuals from the selected 80 households in the ongoing Ha Nam household cohort study on influenza was available: 189 (52.6%) individuals were female, 192 (53.5%) individuals belonged to the working age group (19 to 55 years of age), and 61 (17.0%) were children under 11 years of age. With 154 (42.9%) farmers, farming was the most common occupation. Thirty-one (8.6%) of the participants had an underlying chronic illness. Two hundred sixty-fivepeople were willing to participate in this study, 94 of 359 people were not willing or unable to participate in this study. The distribution of demographic characteristics of participants was comparable to that of the total population (Table 1a).

Table 1 Study population characteristics (1.a) Age distribution occupation and antibiotic use, (1.b) direct environmental characteristics of households

All 80 households had at least one animal, with dogs in 63 (79%) of the households and chickens in 56 (70%) being most common. In 62 (70%) households, at least one animal species other than dogs or chicken was present, mostly pigs and cows. Rainwater was the main water supply for 75 (94%) households for cooking and well water was the main water source in 74 (93%) households for daily activities such as washing clothes, vegetable, raw meat. We observed that 38 (48%) households used fertilizer in vegetable culture, of which 14 households (37%) used fertilizer of human or animal fecal origin (Table 1b).

Antibiotic consumption in the complete cohort during 6 months of the study follow-up

During the 6 months follow-up period, 97 (27.0%) participants belonging to 37 households used antibiotics at least once. 55 (56.7%) of them reported use of a single course of antibiotics, 23 (23.7%) used two courses of antibiotics and 18 (18.6%) individuals used more than two courses during follow-up. Antibiotics were obtained from various sources in the drug supply system: national level referral hospitals (3%, n = 4), provincial level hospital (8%, n = 15), community level healthcare center (42%, n = 74), private pharmacy (40%, n = 71) or through antibiotic-sharing (8%, n = 14). Beta-lactam antibiotics were most commonly used: 79 (82%) of 96 individuals who consumed antibiotic, had at least one course of beta-lactam antibiotics (1st generation cephalosporin [43%], amoxicillin [22%], 2nd or 3rd generation cephalosporin [17%], penicillin [9%]).

Resistome of a subset of cohort at the fourth month follow-up time-point

A total of 304 human, animal, food and environmental samples were collected at month 4 from 40 households (see above). Samples were subjected to qPCR assays to assess the presence and abundance of ARGs. A total of 14 unique ARGs conferring resistance to different classes of antimicrobials were targeted. The number of detected resistance genes per individual sample ranged between 1 and 8. QnrS dominated in human, animal and environmental samples (220/304, 72.4%), followed by mcr-1 (134/304, 44.1%), CTX-M (127/304, 41.8%) and mcr-3 (91/304, 29.9%). Other commonly detected genes were: qnrA (39/304, 12.8%), blaNDM (25/304, 8.2%), mcr-2, mcr-4, mcr-5 and blaOXA-48 were found only sporadically in the various sample types (Table 2). Proportions of ARG positive human and animal samples were significantly higher than of cooked food and water. In particular, qnrS was found in all (93/93) of the human fecal samples and in 91% (41/45) of animal stools whereas 75% (57/75) and 32% (29/91) of the cooked food and water samples were positive, respectively. In contrast, the prevalence of qnrA was low in all sample types except animal faeces (22/45, 49%). Mcr-1 was detected at a high frequency in animal (42/45, 93%) and human stools (82/93, 88%), and at lower rates in cooked food (9/75, 12%) and water (1/91, 1%).

Table 2 The prevalence of AGRs among human, animal, food and water obtained from a subset 304 samples from 40 households

Mcr-3 was found in 55% (51/93) of human stools followed by 51% (23/45) in animal stools, 28% (21/75) in cooked food and 7% (6/91) in water samples. Human stools also frequently contained the following ESBL-encoding genes: CTX-M-1 (38%, 35/93), CTX-M-2(18%, 17/93) and CTX-M-9 (61%, 57/93). The proportion of CTX-M-1 and CTX-M-2 (13%, 6/45) positive animal stools was similar to the positive proportion of human samples, whereas the proportion positive of CTX-M-9 was much lower (p < 0.001) in animal stools (29%, 13/45). The proportion of CTX-M-1, CTX-M-2 and CTX-M-9 groups in cooked food and water were much lower (p < 0.001) than in human and animal stools (Table 2).

As opposed to the other ARGs, the proportion positive for carbapenemase encoding genes was highest in food and water samples. The carbapenemase encoding gene blaNDM was detected in 13% (10/75) of food and 15% (14/91) of water samples, but only in 4% (4/93) of human and 4% (2/45) of animal stools. Moreover, while blaOXA-48 was absent in human and animal samples, 8% (6/75) of food and 2% (2/91) of water samples were positive for this carbapenemase gene. BlaKPC and blaVIM were not found in any of the samples (Table 2).

ARG in relation to antibiotic consumption in humans

During the 4 months prior to sampling, 44 (47%) of 93 participants included for the resistome analysis took antibiotics. Table 3 shows the proportion of fecal samples from participants who did and did not use antibiotics positive for specific ARGs. The proportion of qnrA was significantly lower (1/49 versus 7/44) among individuals that consumed antibiotics in the 4 months prior to sample collection (Fisher exact, p –value < 0.05). The proportion of the remaining ARG did not differ significantly between the two groups (Table 3).

Table 3 The proportion of ARGs in the comparison between humans who used antibiotic and those who did not use antibiotics

In terms of ARG abundance in fecal samples from participants who used versus those who did not use antibiotics, we only observed a statistically significant higher relative abundance of the CTX-M-1 group among participants who used antibiotics (p < 0.05). For the remaining ARGs no statistically significant differences were found (Fig. 2).

Fig. 2
figure 2

Comparison of the relative abundance of ARGs between people who used antibiotics versus those who did not use antibiotics in the 4 months preceding sample collection. *indicates p < 0.05 as determined by Wilcoxon test

We subsequently conducted a principal component analysis (PCA) to visualize the ordination of ARGs abundance between sample types. ARGs profiles of human and animal stools appeared to be most strongly separated from the other sample types in both the first two principal components (PC1 and PC2) which accounted for the first and the second largest possible variance in our data set (Fig. 3a). Testing for potential differences confirmed that scores on PC1 significantly differed between all sample types, except for human and animal samples. Also, the scores on PC2 are significantly different between most of the sample types, except for water and vegetables, indicating that the resistome differs between sample types (Fig. 3c). In order to identify how the sample types differed in resistome composition, we next plotted the factor loadings (the correlation coefficients between the ARGs abundance and sample types) of the PC1 and PC2. The plot shows that the abundance of mcr-1, mcr-3 and CTX-M-1 genes are highly loaded on PC1 and abundance of blaNDM, blaOXA-48, and mcr-5 genes dominated on PC2 (Fig. 3b). This revealed that mcr-1, mcr-3, and CTX-M-1 variations are mainly driving the separation of the different sample types on PC1, whereas blaNDM blaOXA-48, and mcr-5 are contributing to the separation on PC2 (Fig. 3b).

Fig. 3
figure 3

Principal Component Analysis (PCA) plots (a) showing the ordination of samples according to their ARGs profiles. Samples are colored according to their origin: animal feces (red), human feces (green), meat (dark blue), vegetable (light blue) or water (purple). The ellipses represent a 95% confidence interval of samples of the same origin. Loading scores (b) of the ARGs on the first and second principal component. Boxplots (c) showing the scores on the first and second principle component of samples of the same origin


We used a targeted metagenomics approach to explore the resistome of humans and animals and their environment in northern Vietnam. Our study revealed highly variable detection rates between sample types, indicating that resistome profiles are habitat-dependent.

The detection of mcr-1 in both human (88%, 82/93) and animal (93%, 42/45) stools was higher than in other parts of the world. Using similar metagenomic screening, the acquisition rate of mcr-1 was approximately 4.9% (6/122) among Dutch intercontinental travelers, while none of the participants carried mcr-1 prior to travel departure [26]. A culture-based screening failed to detect mcr-1 positive Enterobacteriaceae in healthy Swiss individuals [4]. Although previous studies in Asian populations also showed high mcr-1 detection rates, proportions were still much lower than in our study (ranging from 0.4 to 2%). A recent study conducted in southern Vietnam reported mcr-1 carriage among 21% of humans and 59% of chickens [29]. The latter study reported highest human mcr-1 colonization rates among farmers exposed to mcr-1 positive chickens, suggesting zoonotic transmission. A study by Chen et al., that investigated the presence of mcr-1 in the ecosystem in China from 2015 to 2016 showed 51% (234/480) of animal faeces were mcr-1 positive [30]. Interestingly, this study reported much higher mcr-1 rates in water (71%, 24/34) and food (36%, 486/1371) than our findings [30].

The high mcr-1 detection rates in our study when compared to previous studies conducted in Vietnam and China is most likely due to widespread use of colistin in animal feed but could also partly be attributed to the culture-independent targeted metagenomic approach, which is known to have a higher sensitivity when compared to culture-based approaches, and also detects ARGs from non-viable bacteria and extracellular DNA. Using our approach provides better insight on colistin resistance gene abundance in a sample type rather than a culture based approach of indicator bacteria.

The mcr-3 gene is evolutionarily distinct from mcr-1 [31] and both have successfully spread among Enterobacteriaceae in Vietnam, as previous studies indicated [29, 32]. Studies about human fecal shedding of mcr-3 in the community are limited. Here, the detection rate of mcr-3 (55%, 51/93) in human feces was much higher than reported previously (5% [8/152] from stool samples from outpatients in China using culture based approaches) [33].

Examining the overall resistome profiles indicated that the separation of human and animal fecal samples from the environmental samples was mainly driven by mcr-1 and mcr-3. There is limited evidence of ESBLs transmission from animals and meat to humans [34,35,36], but our results suggest its likely that mcr ARGs are transferred from animals to humans as colistin is only used in animal husbandry [7, 37]. Altogether our results advocate a comprehensive assessment of the acquisition and subsequent transmission of mcr genes between animals and humans.

The detection rate of blaNDM in human feces was relatively low, and comparable to the prevalence of Gram-negative bacilli carrying blaNDM cultured from throat swab in humans from rural northern Vietnam [38]. However, it is interesting to note that blaNDM and blaOXA-48 were more frequently detected in well and irrigative water and food than in human and animal feces. Possibly, carbapenemase genes in environmental bacteria as Acinetobacter spp. are much more common than in human and animal intestinal bacteria such as E. coli and this may be a result of higher abundance of these bacteria in aquatic environments. Possibly, residues of antibiotics and other resistance gene driving pollutants (e.g. heavy metals) may be present in this environment selecting for resistance [39].

Antibiotic use in humans and animals in Vietnam is largely uncontrolled [40, 41] as reflected by the findings from our cohort where beta-lactams (mostly cephalosporins) and other antibiotics were very commonly used. However, we only found a slightly higher richness of CTX-M-1 genes, and no differences for other genes, among individuals with recent use of beta-lactam antibiotics. The long-term frequent exposure to antibiotics may have led to a ‘stable’ omnipresence of ARGs. Alternatively, increases in the number of ESBL producing E.coli and Salmonella spp also occur without any prior use of cephalosporins [42], potentially explaining the lack of association between recent antibiotic use and the presence and abundance of most ARGs in our study.

The present study was mainly explorative and has several limitations. Because of the cross-sectional nature of the analyzed sample collection it is impossible to study the transmission dynamics of ARGs, as well as the direct impact of antibiotic use on the ARG profile. Moreover, compared with culturing methods, a metagenomic approach has the advantage of being able to detect resistance in a much wider array of species; however, we did not assess in which organisms the ARG detected in our study were present, nor if they were being expressed.

The targeted metagenomics approach also provides higher sensitivity to detect ARGs than culture-based methods and sequence-based metagenomics [43, 44]. Using sequence-based metagenomics, a recent study monitored antimicrobial resistance in urban sewage from 60 countries [45]. Although this study showed that mcr-1 was detected in only a few countries, including Vietnam, the low relative abundance of mcr-1 and absence of mcr-3 in Vietnamese sewage underestimated the wide dissemination of mcr genes in the Vietnamese community as we describe here. This finding indicates that targeted metagenomics can complement sequence-based metagenomics screening to monitor antimicrobial resistance within local communities and across the globe.

In order to acquire a more detailed insight into the impact of antibiotic use on resistome profiles and to understand the dynamics of the resistome in various microbial ecosystems, longitudinal analyses are needed. These analyses are underway using samples from our cohort collected at different time points. Moreover, additional profiling of the microbiome composition will enable more accurate analysis of the correlation between resistome and taxonomic changes in response to antibiotics in a high ARG and usage background.


Our study provides a comprehensive view of ARGs in humans, animals, and their environment in rural northern Vietnam and reveals a high level of circulating ARGs with different resistome profiles per sample type. Mobilized colistin resistance genes (mcr-1 and mcr-3) were detected in human feces at 88% (82/93) and 55% (51/93) and animal feces at 93% (42/45) and 51% (23/45), respectively. Carbapenemase genes in the environment are much more common than in human and animal feces. Although beta-lactam antibiotics were used most frequently in our study community, their impact on the presence and abundance of ARGs was limited and needs to be examined in more detail in a longitudinal study.

Availability of data and materials

Please contact to Corresponding Author. *email:



Antimicrobial resistance


Antibiotic resistant genes


Extended spectrum beta-lactamase


Gross Domestic Product


Mobilized colistin resistance


Principal component 1


Principal component 2


Principal component analysis


  1. Spellberg B, Gilbert DN. The future of antibiotics and resistance: a tribute to a career of leadership by John Bartlett. Clin Infect Dis. 2014;59 Suppl 2(Suppl 2):S71–5.

    Article  PubMed  CAS  Google Scholar 

  2. Wright GD. The antibiotic resistome: the nexus of chemical and genetic diversity. Nat Rev Microbiol. 2007;5:175.

    Article  CAS  PubMed  Google Scholar 

  3. Bordier M, et al. Antibiotic resistance in Vietnam: moving towards a one health surveillance system. BMC Public Health. 2018;18(1):1136.

    Article  PubMed  PubMed Central  Google Scholar 

  4. CDDEP (The Center for Disease Dynamics, E.P.). State of the World's antibiotics, 2015. Washington, D.C: CDDEP; 2015.

    Google Scholar 

  5. Ventola CL. The antibiotic resistance crisis: part 1: causes and threats. P & T Peer-Rev J Formul Manage. 2015;40(4):277–83.

    Google Scholar 

  6. Bui TMH, et al. Carriage of Escherichia coli producing CTX-M-type extended-Spectrum β-lactamase in healthy Vietnamese individuals. Antimicrob Agents Chemother. 2015;59(10):6611–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Van Cuong N, et al. Antimicrobial consumption in medicated feeds in Vietnamese pig and poultry production. Ecohealth. 2016;13(3):490–8.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Vietnam, G.s.o.o. Statistical yearbook of Viet Nam 2016. Hanoi: Statistical publishing house; 2017.

    Google Scholar 

  9. Dang Pham Kim CS, Douny C, Dinh TV, Xuan BH, Vu BD, Hong NP, et al. First survey on the use of antibiotics in pig and poultry production in the red River Delta region of Vietnam. Food Public Health. 2013;3(5):247–5.

    Google Scholar 

  10. Hanoi, V.A.o.S.E.a.P.S.f., Decree 39/2017/ND-CP. Providing the regulatory framework for animal feeds and aqua feeds.Vietnam: M.o.A.a.R. Development; 2017.

  11. Hoa PT, et al. Antibiotic contamination and occurrence of antibiotic-resistant bacteria in aquatic environments of northern Vietnam. Sci Total Environ. 2011;409(15):2894–901.

    Article  CAS  PubMed  Google Scholar 

  12. Hu Y, et al. The abundance of antibiotic resistance genes in human guts has correlation to the consumption of antibiotics in animal. Gut Microbes. 2014;5(2):245–9.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Sommer MO, Church GM, Dantas G. The human microbiome harbors a diverse reservoir of antibiotic resistance genes. Virulence. 2010;1(4):299–303.

    Article  PubMed  Google Scholar 

  14. Hu Y, et al. Metagenome-wide analysis of antibiotic resistance genes in a large cohort of human gut microbiota. Nat Commun. 2013;4:2151.

    Article  PubMed  CAS  Google Scholar 

  15. Forslund K, et al. Country-specific antibiotic use practices impact the human gut resistome. Genome Res. 2013;23(7):1163–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Authority, E.F.S., E.C.f.D. Prevention, and control, The European Union summary report on antimicrobial resistance in zoonotic and indicator bacteria from humans, animals and food in 2014. EFSA J. 2016;14(2):4380.

    Google Scholar 

  17. Horby P, et al. The epidemiology of interpandemic and pandemic influenza in Vietnam, 2007-2010: the Ha Nam household cohort study I. Am J Epidemiol. 2012;175(10):1062–74.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Tedjo DI, et al. The effect of sampling and storage on the fecal microbiota composition in healthy and diseased subjects. PLoS One. 2015;10(5):e0126685.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  19. Luke Thompson, G.A., Greg Humphrey, Jack Gilbert, Janet Jansson, Rob Knight. EMP Sample Submission Guide.

  20. Costea PI, Zeller G, Sunagawa S. Towards standards for human fecal sample processing in metagenomic studies. Nat Biotechnol. 2017;35(11):1069–76.

    Article  CAS  PubMed  Google Scholar 

  21. Klindworth A, et al. Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies. Nucleic Acids Res. 2013;41(1):e1.

    Article  CAS  PubMed  Google Scholar 

  22. Birkett CI, et al. Real-time TaqMan PCR for rapid detection and typing of genes encoding CTX-M extended-spectrum beta-lactamases. J Med Microbiol. 2007;56(Pt 1):52–5.

    Article  CAS  PubMed  Google Scholar 

  23. Naas T, et al. Real-time PCR for detection of blaOXA-48 genes from stools. J Antimicrob Chemother. 2013;68(1):101–4.

    Article  CAS  PubMed  Google Scholar 

  24. Naas T, et al. Real-time PCR for detection of NDM-1 Carbapenemase genes from spiked stool samples. Antimicrob Agents Chemother. 2011;55(9):4038–43.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Vien le TM, et al. The co-selection of fluoroquinolone resistance genes in the gut flora of Vietnamese children. PLoS One. 2012;7(8):e42919.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  26. von Wintersdorff CJH, et al. Detection of the plasmid-mediated colistin-resistance gene mcr-1 in faecal metagenomes of Dutch travellers. J Antimicrob Chemother. 2016;71(12):3416–9.

    Article  CAS  Google Scholar 

  27. Wintersdorff CJ, et al. The gut resistome is highly dynamic during the first months of life. Future Microbiol. 2016;11(4):501–10.

    Article  CAS  Google Scholar 

  28. von Wintersdorff CJ, et al. High rates of antimicrobial drug resistance gene acquisition after international travel, the Netherlands. Emerg Infect Dis. 2014;20(4):649–57.

    Article  CAS  Google Scholar 

  29. Trung N, Matamoros S, Carrique-Mas JJ, Nghia N, Nhung N, Chieu T, et al. Zoonotic transmission of mcr-1 Colistin resistance gene from small-scale poultry farms, Vietnam. Emerg Infect Dis. 2017;23(3):529–32.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Chen K, et al. Widespread distribution of mcr-1-bearing bacteria in the ecosystem, 2015 to 2016. Eurosurveillance. 2017;22(39):17–00206.

    Article  PubMed Central  Google Scholar 

  31. Wang X, et al. Presence of an <em>mcr-3</em> variant in <span class="named-content genus-species" id="named-content-1">Aeromonas caviae</span>, <span class="named-content genus-species" id="named-content-2">Proteus mirabilis</span>, and <span class="named-content genus-species" id="named-content-3">Escherichia coli</span> from one domestic duck. Antimicrob Agents Chemother. 2018;62(2):e02106–17.

    PubMed  PubMed Central  Google Scholar 

  32. Snesrud E, et al. A model for transposition of the Colistin resistance gene mcr-1 by ISApl1. Antimicrob Agents Chemother. 2016;60(11):6973–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Sun Q, et al. Alkaline peptone water-based enrichment method for mcr-3 from acute diarrheic outpatient gut samples. Front Med (Lausanne). 2018;5:99.

    Article  Google Scholar 

  34. Ljungquist O, et al. Evidence of household transfer of ESBL−/pAmpC-producing Enterobacteriaceae between humans and dogs - a pilot study. Infect Ecol Epidemiol. 2016;6:31514.

    PubMed  Google Scholar 

  35. Mevius D, et al. ESBL-attribution analysis (ESBLAT): looking for sources of antibiotica resistance in humansIn Dutch: Rapport ESBL-Attributieanalyse (ESBLAT) : Op zoek naar de bronnen van antibioticaresistentie bij de mens; 2018.

    Google Scholar 

  36. Trung Nguyen V, et al. Limited contribution of non-intensive chicken farming to ESBL-producing Escherichia coli colonization in humans in Vietnam: an epidemiological and genomic analysis. J Antimicrob Chemother. 2019;74:561–70.

    Article  CAS  Google Scholar 

  37. Nguyen NT, et al. Use of Colistin and other critical antimicrobials on pig and chicken farms in southern Vietnam and its association with resistance in commensal Escherichia coli bacteria. Appl Environ Microbiol. 2016;82(13):3727–35.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Dao TT, et al. Klebsiella pneumoniae oropharyngeal carriage in rural and urban Vietnam and the effect of alcohol consumption. PLoS One. 2014;9(3):e91999.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  39. Binh VN, et al. Antibiotics in the aquatic environment of Vietnam: sources, concentrations, risk and control strategy. Chemosphere. 2018;197:438–50.

    Article  CAS  PubMed  Google Scholar 

  40. Van Kinh N, et al. Situation analysis: antibiotic use and resistance in Vietnam. In: The center for disease dynamics, Economics & Policy; 2010.

    Google Scholar 

  41. VNS. End misuse of antibiotics, combat resistance, experts say, in Viet Nam News; 2017. Viet Nam

    Google Scholar 

  42. Ewers C, et al. Extended-spectrum beta-lactamase-producing and AmpC-producing Escherichia coli from livestock and companion animals, and their putative impact on public health: a global perspective. Clin Microbiol Infect. 2012;18(7):646–55.

    Article  CAS  PubMed  Google Scholar 

  43. Melendez JH, et al. Real-time PCR assays compared to culture-based approaches for identification of aerobic bacteria in chronic wounds. Clin Microbiol Infect. 2010;16(12):1762–9.

    Article  CAS  PubMed  Google Scholar 

  44. Koskinen M, et al. Analytical specificity and sensitivity of a real-time polymerase chain reaction assay for identification of bovine mastitis pathogens. J Dairy Sci. 2009;92:952–9.

    Article  CAS  PubMed  Google Scholar 

  45. Hendriksen RS, et al. Global monitoring of antimicrobial resistance based on metagenomics analyses of urban sewage. Nat Commun. 2019;10(1):1124.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

Download references


This work was supported by Radboudumc Revolving Research Fund (R3Fund) grant, and VIDI grant to J. Penders from the Netherlands Organization for Health Research and Development (ZonMw, grant number 016.156.427). We also would like to express our sincere thanks to the volunteers who participated in this research.


The Radboudumc Revolving Research Fund (R3Fund) and Wellcome Trust Major Asia Programme, UK.

Author information

Authors and Affiliations



HFLW was responsible for conception, study design, and funding application. All authors contributed to the study protocol development. BVTN led day- to- day management of the study implementation supervised by HFLW, THH, JP. HFLW, THH and BVTN took part in getting ethics approval and training for health-care centers. BVTN, PDT, TTVP, MO, CD and EB were responsible for sample collection and other laboratory works. LVT and BVTN were responsible for statistical data analysis supervised by JP, HRVD. Drafting of manuscript was done by BVTN, HFLW, JP and HRVD. All authors contributed to the final revision and approved the submission.

Authors’ information

Not applicable.

Corresponding author

Correspondence to Vu Thi Ngoc Bich.

Ethics declarations

Ethics approval and consent to participate

The research was approved by the Oxford University Tropical Research Ethics Committee (OxTREC, 49–14) and the National Institute of Hygiene and Epidemiology, Vietnam (NIHE) institutional review board. Written informed consent was obtained from all participants in the study.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Additional file

Additional file 1: Supplement 1.

List of primers and probes used to identify and quantify the ARGs by qPCR. Supplement 2. qPCR's results of ARGs identification and quantification.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Bich, V.T.N., Thanh, L.V., Thai, P.D. et al. An exploration of the gut and environmental resistome in a community in northern Vietnam in relation to antibiotic use. Antimicrob Resist Infect Control 8, 194 (2019).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: