We aim to understand the microbial ecology of noma (cancrum oris), a devastating ancient illness which causes severe facial disfigurement in>140,000 malnourished children every year. The cause of noma is still elusive. A chaotic mix of microbial infection, oral hygiene and weakened immune system likely contribute to the development of oral lesions. These lesions are a plausible entry point for unidentified microorganisms that trigger gangrenous facial infections. To catalog bacteria present in noma lesions and identify candidate noma-triggering organisms, we performed a cross-sectional sequencing study of 16S rRNA gene amplicons from sixty samples of gingival fluid from twelve healthy children, twelve children suffering from noma (lesion and healthy sites), and twelve children suffering from Acute Necrotizing Gingivitis (ANG) (lesion and healthy sites). Relative to healthy individuals, samples taken from lesions in diseased mouths were enriched with Spirochaetes and depleted for Proteobacteria. Samples taken from healthy sites of diseased mouths had proportions of Spirochaetes and Proteobacteria that were similar to healthy control individuals. Samples from noma mouths did not have a higher abundance of Fusobacterium, casting doubt on its role as a causative agent of noma. Microbial communities sampled from noma and ANG lesions were dominated by the same Prevotella intermedia OTU, which was much less abundant in healthy sites sampled from the same mouths. Multivariate analysis confirmed that bacterial communities in healthy and lesion sites were significantly different. Several OTUs in the Orders Erysipelotrichales, Clostridiales, Bacteroidales, and Spirochaetales were identified as indicators of noma, suggesting that one or more microbes within these Orders is associated with the development of noma lesions. Future studies should include longitudinal sampling of viral and microbial components of this community, before and early in noma lesion development.
Noma, or cancrum oris, has caused devastating gangrenous oral-lesions throughout human history. It is a rapid and destructive infection that is thought to have a microbial origin in the context of poor immune function. Within 24–48 hours of the emergence of swelling in the cheek, often accompanied by gingivitis, gangrenous infection overtakes a region of the cheek, gums and facial bones. The ability of the infection to overcome structural barriers like muscle and bone is unusual. The term noma is derived from the Greek verb numein (to devour). It is a complex disease which arises most often in young children who are malnourished, exposed to infectious diseases, and lack access to clean water and dental hygiene. Noma has not been prevalent in developed countries since the 20th century, with the exception of concentration camps in WWII, and isolated instances in Europe and the United States in people who are immunocompromised [1]. The WHO estimated a world-wide incidence of 140,000 cases per year in 1998, with 100,000 occurring in children who are 1–7 years old in sub-Saharan Africa [2]. The mortality rate is estimated to be 80–90%, though 770,000 people are thought to have survived the acute phase of infection. This may be an underestimate because noma occurs in remote parts of the world and is not often reported because of the great deal of shame associated with the disease. Current treatment includes hydration and nutrition, local disinfection of the wound, antibiotic therapy, and reconstructive surgery after the wounds have healed. Noma is not considered an infectious disease. There have not been outbreaks reported in villages, and recurrence is rare. In one family, twins who both developed noma were affected in different years [3]. Attempts to infect mammals by inoculation of sub-gingival microbiota were only successful with the addition of steroid injections [4]. The exact etiology of noma is unknown. Several risk factors are thought to work together to increase the risk of developing the disease, including: 1) malnutrition 2) poor oral hygiene, potentially including oral lesions from gingivitis 3) a compromised immune system and 4) an unidentified triggering micro-organism, previously hypothesized to be a member of Fusobacterium or Prevotella [5]–[8]. The type of necrosis and the odor of the lesions are the main evidence supporting the long-held theory that a bacterium triggers the infection [1], [3], [6], [9]. Several limitations have hindered identification of the microbiological cause of noma. Because it develops rapidly in remote regions, obtaining reliable samples during early development of the disease is difficult. Earlier studies of noma relied on microbiological culturing of samples, which biases the results towards organisms that can grow in the culture conditions and may not reveal important members of the microbial community. Also, without knowing the composition of oral bacterial populations in healthy children, or children in similar circumstances but who do not have noma, it is difficult to separate disease-causing organisms from those that are common in the local population. The recent revolution in DNA-based taxonomy methods enabled the discovery of the 99% of microorganisms that are difficult to grow in culture, but thrive in various environments, including the human body [1], [10]. Increased knowledge of these commensal microbial communities has refined our thinking about how microbes cause disease. For more than a century, medical doctors have focused on individual pathogens when studying infectious disease. In reality many infections may originate from, or be influenced by, the thriving communities of understudied commensal microbiota. In many cases, disease may arise from a disturbed microbial community, rather than a single pathogen [11], [12]. The obscure microbial origin of noma infection is well-suited for an investigation of how a shifted microbial community can result in disease. A single microbial trigger may not be revealed, but rather a combination of altered microbial abundances, which lead to disease in the context of malnutrition and a compromised immune system. The interplay with a background of malnutrition may mean that a commensal or opportunistic pathogen that is also present in healthy people triggers the disease. Ideally, microbial communities would be sampled before, during and after the appearance of noma wounds. This is still an elusive goal, as obtaining early samples is challenging under any circumstance, and becomes nearly impossible in rural areas where modern medicine is not accessible. The GESNOMA group has led DNA-based efforts to identify the etiology of noma, having obtained and studied the bacterial communities in 413 samples over the last 10 years [9], [13], [14]. Sampled children in Niger included those deemed healthy, diagnosed with noma, and diagnosed with Acute Necrotizing Gingivitis (ANG), a related periodontal disease characterized by painful infections of the gums (also known as trench mouth) [15]. The ages of healthy controls were matched as closely as possible to those suffering from noma. Bacterial communities in these samples were characterized using DNA-based methods to learn about the epidemiology and microbiology of Noma [1], [9], [13], [14]. These recent GESNOMA studies suggest that the same (or very similar) bacteria are present in noma and ANG lesions. For example, both types of lesion had members of the genera Prevotella, Peptostreptococcus, and Treponema present [13], [14]. The healthy controls used in these studies were crucial for examining the role of Fusobacterium necrophorum, which was considered to be a culprit in the development of noma [5]. GESNOMA studies based on both low-throughput sequencing and microarray hybridization identified Fusobacterium to be as common in healthy controls as in lesions, suggesting its mere presence does not trigger noma [13], [14]. GESNOMA' s low-throughput sequencing and microarray approaches led to discordant estimations of bacterial diversity between healthy and diseased sites: low-throughput sequencing suggested these communities were equally diverse [14], whereas hybridization of community DNA to microarrays indicated that noma had less diversity than healthy controls [13]. Though both of these technologies were powerful in their abilities to detect bacteria without imposing biases associated with culturing, they both lacked the ability to detect microbes of lower abundance, and the microarrays could not detect nucleic acids that were not included in the microarray design. Our current study builds on these previous reports by using high-throughput sequencing to examine the bacterial communities in sixty samples from children in Niger. Improving the sequence coverage of these communities enabled the detection and quantification of bacteria that were less abundant and whose identities were not known in advance. We investigated whether the bacterial communities in lesions from noma differed from those in ANG lesions, compared to healthy gingival samples from the same children and gingival samples from healthy controls. Our goals were to begin cataloging the species present in individuals from the same villages who were healthy, suffering from noma, or suffering from ANG. This catalog will eventually serve as a powerful comparative tool for future longitudinal studies. We used this information to identify candidate bacterial species that are characteristic of acute noma lesions.
As part of the GESNOMA study, subgingival fluid samples were collected with cotton points from noma, acute necrotizing gingivitis and control volunteers by nurses based in Zinder, Niger as described in [13], [14]. Samples were subsequently stored in guanidinium isothiocyanate medium (RLT buffer, Qiagen) at -80°C until processing. This study is part of the GESNOMA project, an epidemiological data and sampling collection effort conducted between September 2001 and October 2006 in the Zinder area of Niger, Africa. Gingival fluids were collected by local nurses from children affected by ANG or noma along with gender matched controls. The children were between 6 months and 12 years of age at the time of sample collection. The study was carried out with the prior approval of the Niger Ministry of Health, the WHO noma programme in Niger, and the Swiss non-governmental organization “Sentinelles”, and complied with the Declaration of Helsinki. A brochure with an explanation of informed consent was prepared in French and translated orally by local staff to the children' s parents or guardians. The parent or guardian of each child signed the consent form, usually by a cross and/or a fingerprint. The consent forms were all sent to the GESNOMA office in Geneva, Switzerland. The Republic of Niger Ministry of Public Health approved the protocol in 2001. The Head of the Department of Studies and Programming of the Republic of Niger Ministry of Public Health approved the use of fingerprints as a signature in the case of illiteracy in advance. Additional information about the children in the study can be found in Table S1 and GESNOMA' s epidemiological study [9]. DNA was extracted from stored samples using the DNeasy Blood & Tissue Kit (Qiagen). To facilitate cell lysis, samples were vortexed for 30 seconds with acid washed glass beads (Sigma #G1277,212–300 µm). Then, 200 µl buffer AL was added and the protocol ‘Purification of Total DNA from Animal Blood’ from the’ DNeasy Blood & Tissue Handbook’ was followed. DNA was eluted in 30 µL of sterile water. To amplify the V1–3 region of the bacterial 16S rRNA gene, 1 µL of DNA was combined with 0. 5 µM of the forward primer 5′–gccttgccagcccgctcagacGAGTTTGATCMTGGCTCAG–3′, 0. 05 µM of an additional forward primer 5′–gccttgccagcccgctcagacAGGGTTCGATTCTGGCTCAG–3′ and 0. 5 µM barcoded reverse primer 5′–gcctccctcgcgccatcagNNNNNNNNatCCGCGRCTGCTGGCAC–3′ in 50 µl PrimeStar HS Premix (Takara). The structure of composite PCR primers is described in [16]. Amplification proceeded for 30 cycles of 98°C for 10 s, 56°C for 15 s, and 72°C for 1 min. To amplify a segment of archaeal 16S rRNA gene, a two-step PCR was carried out to obtain sufficient material for sequencing. Initially, 2 µL of lysate was combined with 0. 3 µM each of the forward 5′–CCGACGGTGAGRGRYGAA–3′ and reverse 5′–ACGGGCGGTGWGTRCAA–3′ primers in 25 µl PrimeStar HS Premix (Takara). Amplification proceeded for 33 cycles of 98°C for 10 s, 69°C for 15 s, and 72°C for 1 min. The second step used a 0. 5 µl aliquot from the first reaction as template in a new PCR with 50 µL PrimeStar HS Premix (Takara) containing 0. 5 µM each barcoded forward primer 5′–gcctccctcgcgccatcagNNNNNNNNtaCCGACGGTGAGRGRYGAA–3′ and reverse primer 5′–gccttgccagcccgctcagacACGGGCGGTGWGTRCAA–3′. Amplification proceeded for 28 cycles of 98°C for 10 s, 58°C for 15 s, and 72°C for 45 s. For each sample, two identical PCRs were pooled and purified using the MinElute PCR purification kit (Qiagen), using 25 uL water for elution. Amplicon sizes and concentrations were estimated by running 1 µl of each amplicon on the Agilent 2100 Bioanalyzer using a DNA1000 lab chip. 100 ng of each purified amplicon was pooled and sequenced on a Genome Sequencer FLX system (Roche). The 91,114 sequence reads were trimmed using Qiime 1. 2. 0 [17]. Sequences were removed if they: contained any ambiguous bases (e. g. , ‘N’), did not match the reverse primer or barcode, or had homopolymer stretches greater than 6 nucleotides. A large percentage of sequences were found to contain ambiguities at the ‘A’ in the second to last position at the 3′ end of the reverse primer; we allowed for only these ambiguities in order to retain more data. After trimming sequences to 200 bp, 71,715 sequences remained. These were binned by the 60 unique barcodes, representing the 60 samples that included 12 samples from each of the 5 health status groups. The sequences were clustered, using uclust implemented in Qiime [17], into OTUs using three similarity cut-offs: 90%, 97% and 99%. Chimeras were identified and removed using Chimera Slayer in Mothur [18]. Taxonomy was assigned to a representative sequence from each OTU using RDP in Qiime, and the output was used to create an OTU abundance table. Prior to OTU clustering, each read was also separately classified using RDP in Qiime. Results from these non-OTU based classifications are shown in Figure 1. A Bray-Curtis dissimilarity matrix based on the OTU abundance table was calculated using Primer 6. 0 [19], and used as input for PERMANOVA [19] and Multi-Dimensional Scaling (MDS). A cytoscape network was built with the make_cytoscape_network. py command in Qiime; network edges were defined by the number of sequences in a particular OTU. The network was visualized with the organic layout in Cytoscape, so that the distance between nodes and edges in the network minimizes the force between them. Indicator Species Analysis (ISA) [20] was carried out using the labdsv package in R [21], [22]. The ISA score is calculated based on the relative abundance and frequency of each species by group [20]. Monte Carlo randomization procedures were then used to test the statistical significance of the highest indicator values. FDR values were estimated using the p. adjust function in R [22].
Sixty samples collected from 36 children living in villages within a few hundred kilometers of Zinder, Niger were chosen for inclusion in the current study. Children were examined and sampled before treatment at the Sentinelles clinic by members of the GESNOMA study group [9], [14]. Gingival fluid from healthy and lesion sites of twelve children with noma and twelve children with ANG were included, along with twelve healthy controls [9], [14]. On average, noma lesions were sampled 11 days after their appearance. These sixty samples form five groups of twelve samples for comparison in this study: Noma lesions (N), healthy gingival samples from a noma affected mouth (NH), Acute Necrotizing Gingivitis lesions (ANG), healthy gingival samples from an Acute Necrotizing Gingivitis affected mouth (ANGH), and unaffected mouths from age, location and gender matched-controls (C). The V1-3 segment of the bacterial 16S rRNA gene was amplified and sequenced (unidirectionally) using 454 pyrosequencing technology. We obtained 91,114 reads; after quality filtering 71,715 reads remained (see Methods), of which 11,438 were unique singletons. The distribution of sequence reads among the five health status categories was quite even, ranging from 13,520 to 15,993 reads per sample category (Table 1). To examine bacterial communities at different taxonomic levels, OTUs were calculated three separate times using 90%, 97%, and 99% sequence identity cut-offs. Here we present results from the 97% cutoff; results from the 90% and 99% cutoffs were similar (Table S4 and S5). OTUs shared between different health statuses varied from 262 to 358 (out of 1352 total; Table S4). OTUs from all three identity cutoffs were used to calculate the dissimilarities between communities using the Bray-Curtis metric. In all three cases, samples from sites similar in their health statuses harbored similar communities (i. e. , lower dissimilarity values, in Table S5). Bacterial diversity was measured by calculating Chao1 and Shannon diversity values for each health status group. Chao1 estimates diversity using only the number of OTUs present, whereas Shannon also considers the relative abundance of each OTU. Both diversity measures indicated that noma lesions support a greater diversity of bacteria than do healthy sites in diseased mouths; however, neither Chao1 or Shannon diversity measures differed significantly between wound and healthy sites (Table 1). This suggests that noma lesions are not dominated by a single pathogen that takes over the community; rather these lesions may develop due to a shift in the abundances of community members. The distribution of phyla among the five sample categories revealed pronounced differences in their average abundances (Figure 1); we briefly summarize these differences here and present results from statistical analyses below. In healthy individuals, over 35% of sequences were assigned to Proteobacteria; in individuals with noma or ANG this dropped to fewer than 30%, regardless of whether a lesion or healthy site was sampled. Conversely, in healthy individuals Firmicutes accounted for less than 15% of reads, whereas in samples from both noma and ANG lesions this phylum accounted for over 20% of reads. Reads mapping to members of the phylum Bacteroidetes were more numerous in noma and ANG lesions, where they accounted for greater than 40% of the total, compared to less than 40% for healthy sites. Spirochaetes accounted for well over 5% of reads in lesion sites from both noma and ANG, which was noticeably higher than its average of ∼1% in healthy sites. The top ten most abundant 97% OTUs shown in Figure 2 reveal important differences in community structure between the lesion sites (noma and ANG) and healthy sites (from unaffected gingival fluid in noma and ANG volunteers, and the healthy volunteers). OTU 274, assigned to Prevotella intermedia, dominated wound sites, and had median values at least 3-fold higher than the next most abundant OTU. OTU 274 was also present in samples from healthy sites, though it was similar in abundance to other OTUs. Multi-dimensional scaling (MDS) of the Bray-Curtis dissimilarity matrix was used to visualize differences in community composition for each sample (Figure 3). The healthy controls had more similar microbial community composition, as reflected by their proximity to each other in the MDS plot. Samples from healthy sites in noma and ANG mouths spanned multivariate space between healthy individuals and samples from noma and ANG lesions. Lesions from noma and ANG affected mouths were variable in their community similarities, and although some samples were similar to healthy sites/controls, most of the lesion samples have unique community composition, as reflected by their distance from the other samples in the MDS projection. Samples from the same mouth did not typically cluster together, rather each sample clustered with samples from other mouths that had similar health statuses. To test whether groups significantly differed in their microbiota, we used PERMANOVA [23], a non-parametric multivariate analysis of variance, to test the null hypothesis that there were no differences between the health status groups (Table 2). Overall, PERMANOVA showed that the sample groups were significantly different from each other (F999 = 3. 56, P = 0. 001). Pairwise comparisons with p-values greater than 0. 05 indicate that those two health status groups were not significantly different in the communities they support. We found that the communities from the lesion sites in noma and ANG mouths were more similar to each other (F999 = 1. 22, P = 0. 08) than they were to healthy sites in the same mouths or healthy mouths. This was also true for the healthy parts of the ANG and noma affected mouths (F999 = 1. 08, P = 0. 245). The control samples were most different from the ANG and N samples (Table 2). Similar results were obtained using Anosim [19] (Table S3). Attempts to understand the origin of noma have identified several bacterial species as candidate instigators of noma lesions. We report their abundances in each of the 5 health status categories as numbers of 97% OTUs in Table 3. Although a species in the genus Fusobacterium has been hypothesized to trigger noma, we found this phylum to be more abundant in the controls (Figure 4B). This pattern was also seen for the following hypothesized noma-triggering bacterial genera: Porphyromonas, Capnocytophaga, Aggregatibacter and Streptococcus. The opposite pattern, where abundances were higher in lesions than healthy controls, was observed for the phylum Spirochaetes, and the genera Prevotella and Peptostreptococcus (Figure 4A). Several organisms identified in noma infections by past studies were not found to be more abundant in noma in our dataset, including Corynebacterium, Bacteroides, and Bacillus (Table 3, [1]). To determine which taxa were significantly associated with noma lesion sites, we used indicator species analysis (ISA) [20]. The 12 noma lesion samples were compared against all other samples, including the healthy site samples from the same 12 individuals, and the ANG lesion site samples which had similar composition. ISA calculates the relative abundance and relative frequency of each species in each group, and compares this to a null distribution obtained by randomizing the taxon-sample relationships. These calculations were carried out separately with all three OTU sets (i. e. , 90%, 97%, and 99%). We have reported results from the 97% OTUs here (Table 4); similar results were obtained with the 90% and 99% OTUs. The ISA identified several families that were significantly associated with noma lesions, which overall are consistent with the abundance plots shown in Figure 1. Sharpea (Family: Coprobacillaceae) was the taxon most discerning of noma, and its taxonomic Order Erysipelotrichales was consistently present as a significant indicator no matter which OTU cutoff was employed. Also associated with noma were several OTUs within the genera Veillonella, Peptostreptococcus, Prevotella, and Treponema. These genera are members of the Firmicutes phylum, except for the latter two, which are members of Bacteroidetes and Spirochaetes, respectively. Genera with OTUs indicative of non-noma sites included Capnocytophaga and Aggregatibacter, of the phyla Flavobacteria and Proteobacteria, respectively. Separate OTUs in the Veillonella and Prevotella genera were present in both noma and non-noma sites, indicating that species or strains in these genera differ in their ability to affect lesion development. We specifically sought archaeal sequences in these samples, by repeating the sequencing using primers appropriate for the archaeal 16S rRNA gene. These sequences were analyzed with the same approach as sequences obtained using bacterial primers. Although all individuals produced amplicons with the archaeal primers, the vast majority of sequence reads mapped to bacterial taxa. Presumably this was due to little or no archaeal DNA being present, which then led to non-specific priming and amplification of bacterial templates. However, 15 individuals did produce some archaeal reads, ranging from 1 to 174 reads. Four individuals had 10 or more reads that corresponded to the archaeal phylum Euryarchaeota, with the highest read counts corresponding to the genus Methanobrevibacter. To obtain a global view of the distribution of taxa across the health status groups, and at different levels of similarity, 99% OTUs and the top 35 OTUs with 90% identity were separately visualized in a Cytoscape network using the Qiime network building algorithm (Figure 5) [17], [24]. Individual samples and OTUs are displayed in the network, with lines connecting OTUs to the samples where they are present; shorter distances between the sample and OTU nodes reflect larger numbers of sequences from that OTU found in the connected sample. The OTUs that occur in the center of the network are more evenly distributed between all the samples, while those on the edges occur preferentially in certain samples. Spirochaetaceae and Prevotellaceae are both examples of taxa that were more prevalent in the lesion site samples, and they appeared tightly connected to disease site samples (Figure 5B). Overall, the courser grouping and inclusion of only the top 35 most abundant taxa led to a network where each node is connected to many more neighbors (38. 5 neighbors in the 90% OTU network vs 6. 5 in the 99% network). However strong clustering of sample site by color supports the statistically significant differences in community structure found using PERMANOVA (Table 2).
The 16S rRNA gene amplicon sequence data are available at the MG-RAST website (http: //metagenomics. anl. gov/) with ID# 4585252. 3. The fasta sequence file available there contains output from the ‘split. libraries. py’ command in Qiime; each fasta header includes information identifying which sample the sequence was derived from. The samples are further described, using the same identifiers found in the fasta file headers, in Table S1.
Noma is a traumatic disease characterized by oral-facial lesions that often lead to severe disfigurement and ultimately shame and isolation from the community. Because the causes of noma are likely to be numerous, and reaching those who suffer from this illness is challenging, the etiology of noma remains ill-defined. Although it is known that oral hygiene and nutrition influence the development of noma, evidence suggests that one or more microbes play a crucial role in development of noma lesions. Previous studies have examined the DNA of microbes in lesions to determine which species are present and how their abundances differ between healthy mouth sites and noma lesions. These studies used techniques that were state-of-the-art at the time, though we know they likely only scratched the surface of the resident microbial diversity. Here we extend these studies by digging deeper to characterize a larger diversity of microbial species in noma and control samples, with the goal of better identifying which microbes are uniquely present or have altered abundances in noma lesions.
infectious diseases
medicine and health sciences
nutrition
ecology
oral medicine
biology and life sciences
tropical diseases
microbiology
molecular biology
health care
Noma Affected Children from Niger Have Distinct Oral Microbial Communities Based on High-Throughput Sequencing of 16S rRNA Gene Fragments