microorganisms - MDPI

6 downloads 0 Views 6MB Size Report
Jul 28, 2018 - Ushida, K.; Hoshi, S.; Ajisaka, K. 13C-NMR studies on lactate metabolism in a porcine gut microbial ecosystem. Microb. Ecol. Health Dis. 2002 ...
microorganisms Article

Cecal Microbiome Analyses on Wild Japanese Rock Ptarmigans (Lagopus muta japonica) Reveals High Level of Coexistence of Lactic Acid Bacteria and Lactate-Utilizing Bacteria Atsushi Ueda 1 , Atsushi Kobayashi 2 , Sayaka Tsuchida 3,4 ID , Takuji Yamada 1,5, *, Koichi Murata 6 , Hiroshi Nakamura 7 and Kazunari Ushida 3,4, * 1 2 3 4 5 6 7

*

Department of Life Science and Technology, School of Life Science and Technology, Tokyo Institute of Technology, Tokyo 152-8550, Japan; [email protected] Faculty of Science, Toho University, Tokyo 143-8540, Japan; [email protected] Graduate School of Life and Environmental Sciences, Kyoto Prefectural University, Kyoto, Kyoto Prefecture 606-8522, Japan; [email protected] Academy of Emerging Sciences, Chubu University, Kasugai, Aichi Prefecture 487-0027, Japan PRESTO, Japan Science and Technology Agency, 4-1-8 Honcho Kawaguchi, Saitama 332-0012, Japan Department of Bioresource Sciences, Nihon University, Kanagawa 252-0880, Japan; [email protected] General Foundation Hiroshi Nakamura International Institute for Ornithology, Nakagosho, Nagano 380-0934, Japan; [email protected] Correspondence: [email protected] (T.Y.); [email protected] (K.U.); Tel.: +81-3-5734-3629 (T.Y.); Tel.: +81-568-51-9520 (K.U.)

Received: 26 June 2018; Accepted: 25 July 2018; Published: 28 July 2018

 

Abstract: Preservation of indigenous gastrointestinal microbiota is critical for successful captive breeding of endangered wild animals, yet its biology is poorly understood. Here, we compared the cecal microbial composition of wild living Japanese rock ptarmigans (Lagopus muta japonica) in different locations of Japanese mountains, and the dominant cecal microbial structure of wild Japanese rock ptarmigans is elucidated. Coriobacteraceae and Lachnospraceae were the two dominant bacterial families in all samples analyzed. At the genus level, 10 genera Olsenella, Actinomyces, Megasphaera, Slackia, Cloacibacillus, Bifidobacterium, Escherichia, Dialister, Megamonas, and Bilophila were dominant. These results reveal the high level of coexistence of lactic acid bacteria (Olsenella and Bifidobacterium) and lactate-utilizing bacteria (Megasphaera). This coexistence should be taken into account for the successful breeding of captive Japanese rock ptarmigans in the national conservation program. Keywords: Japanese rock ptarmigan; Lagopus muta japonica; cecal microbiome; Olsenella; Bifidobacterium; Megasphaera

1. Introduction Japanese rock ptarmigans (Lagopus muta japonica) are birds that only inhabit alpine areas of Japan’s main island. The birds were registered as one of the National Treasures of Japan, and recognized as endangered species according to recent population studies [1]. The number of the birds declined to about 1700 in the early 2000s, from as much as 3000 recorded in the 1980s. In the Southern Japanese Alps at Mt. Kita there were 33 recognized territories of the birds in 1981, which were reduced to four in 2004 [2]. Therefore, the national conservation program has been endorsed so far [3].

Microorganisms 2018, 6, 77; doi:10.3390/microorganisms6030077

www.mdpi.com/journal/microorganisms

Microorganisms 2018, 6, 77

2 of 11

The birds have a pair of big cecum, in which dense bacterial populations efficiently ferment the food grazed by the host. Therefore, the cecal microbiota are essential to the birds. We have previously compared the composition of cecal microbiota between wild Japanese rock ptarmigans living in the Murodo area of Mt. Tateyama and captive-bred Svalbard rock ptarmigans (Lagopus muta hyperbrea) in a Japanese zoo [4]. In that study, huge differences in cecal microbial composition between wild and captive birds were observed. The differences in microbial communities between these birds might be affected by genetic factors, but we speculated that the majority of difference comes from differences between wild and captive habitats. In fact, the daily dose of tetracycline for artificially raised captive chicks must have had a tremendous effect on the development of gut microbiota. In addition, predominant bacterial taxa were revealed in Murodo individuals, which may involve the efficient degradation of plant materials grazed by the host. Indeed, we isolated several bacteria which had an ability to detoxify chemical compounds contained in natural foods that the birds adapt to [5,6]. Previous articles related to other birds showed that these cecal microbial compositions are linked to geographical differences [7]. Regarding the wild Japanese rock ptarmigan, the previous molecular research reported that there was a genetic difference between groups living in different mountain areas, which was probably due to the geographical isolation of habitats [8]. In addition, food plants to which they adapt in each habitat were not necessarily the same [9]. This evidence may give rise to geographical differences of microbial composition. However, the difference in cecal microbial composition of wild individuals living in different locations is still poorly understood. Therefore, in the present study, fresh cecal feces of wild adult Japanese rock ptarmigans in different locations of Northern and Southern Japanese Alps mountain regions were collected, and comparison analyses of cecal microbial composition were performed. We showed a high level of coexistence of lactic acid bacteria and lactate-utilizing bacteria in cecal microbiota of the birds, and this observed result should be taken into account for the successful breeding of captive Japanese rock ptarmigans in the national conservation program. 2. Materials and Methods 2.1. Sample Collections Fresh cecal feces of adult birds in different locations were collected aseptically in a conservation buffer as described previously [4]. We tried to obtain the cecal feces just after defecation, and so we watched rock ptarmigans at a short distance (within 10 m) for hours whenever possible. If it was difficult for us to approach, we watched rock ptarmigans at a relatively long distance (up to 50 m). Although defecation was not necessarily confirmed, feces were certainly collected 10 to 15 min after defecation. This is because we quickly and continuously followed the birds and collected cecal feces from the trails of the birds. Feces were collected from different territories to avoid redundant sampling. In Northern Japanese Alps mountain regions, rock ptarmigans in the following areas were examined (each number of samples is in parenthesis): Mt. On [35◦ 540 22” N, 137◦ 290 00” E] (3), Mt. Norikura [36◦ 060 34” N, 137◦ 330 19” E] (3), Mt. Yake [36◦ 130 30” N, 137◦ 350 19” E] (1), Mt. Jonen [36◦ 200 03” N, 137◦ 430 38” E] (1), Mt. Otensho [36◦ 210 47” N, 137◦ 420 12” E] (3) and Mt. Tateyama [36◦ 340 54” N, 137◦ 350 32” E] (8) (Figure 1). In Southern Japanese Alps mountain regions, samples were collected from three adult females under a cage protection program at Mt. Kita [35◦ 390 49” N, 138◦ 130 54” E] (3) in July 2016. Further details of sample collections are shown in a supplemental Table S1 and Figure 2; samples were collected during the non-snowing season from May to September 2016, and most of samples were collected from the rock or dry podzolic soil. The samples from Mt. Tateyama were collected from snow surfaces, even in May, because there was still unmelted snow in the sampling area. We also collected feces on the snow in Mt. Tateyama early in November 2016. Most of the samples were obtained from male individuals, except for those from Mt. Kita. In the case of adult females in Mt, Kita, three families (hen and her chicks) were individually housed in cages to avoid eventual predation during the night and sometimes during the daytime in severe bad weather (strong

Microorganisms 2018, 6, 77 Microorganisms 2018, 6, x FOR PEER REVIEW Microorganisms 2018, 6, x FOR PEER REVIEW

312 of 11 33ofof12

eventualpredation predationduring duringthe thenight nightand andsometimes sometimesduring duringthe thedaytime daytimein insevere severebad badweather weather(strong (strong eventual

wind and rain). thefamilies familieswere were allowed to graze the natural food outside of their wind and rain).Otherwise, Otherwise,the the allowed tograze graze thenatural natural foodoutside outside oftheir theircages. cages. wind and rain). Otherwise, families were allowed to the food of cages. In cages, the birds were offered foods such as Vaccinium ovalifolium, Oxytropis japonica var. In cages, cages, the the birds birds were were offered offered foods foods such such as as Vaccinium Vaccinium ovalifolium, ovalifolium, Oxytropis Oxytropis japonica japonica var. var. japonica, japonica, In japonica, Polygonum viviparum, and Stellaria nipponica, which collected in the surrounding area. Polygonum viviparum, and Stellaria Stellaria nipponica, which werewere collected in the the surrounding area. Polygonum viviparum, and nipponica, which were collected in surrounding area. Previously collected fruits of Vaccinium vitis-idaea andand commercially available warm (Tenebrio sp.) sp.) were Previously collected fruits of Vaccinium vitis-idaea commercially available warm (Tenebrio Previously collected fruits of Vaccinium vitis-idaea and commercially available warm (Tenebrio sp.) also supplied as supplemental food. food. were alsosupplied supplied assupplemental supplemental food. were also as

Figure 1. Locations ofsample samplecollections. collections. Figure1. 1.Locations Locations of of Figure sample collections.

Figure 2. Japanese rock Ptarmigan (Left: A male Japanese rock Ptarmigan on the unmelted snow at Mt. Tateyama in spring. Right: A female Japanese rock ptarmigan with her chicks at Mt. Kita in summer.).

Microorganisms 2018, 6, 77

4 of 11

2.2. DNA Extractions and High-Throughput Sequencing of 16s rRNA Gene Bacterial DNA was extracted from preserved feces as described previously using the Quick Gene DNA Tissue kit (Kurabo, Tokyo, Japan) [4]. The concentration of DNA was determined using both Nano-Drop 1000 and Quant-iT dsDNA HS assay kits with a Qubit fluorometer (Invitrogen, Carlsbad, CA, USA). Microbial community structure was analyzed by 16S rRNA gene amplicon sequencing using an Illumina Miseq with the MiSeq Reagent Kits v3 MS-102-3003, together with PhiX Control v3 at BGI Japan (Kobe, Japan), and low quality reads were removed at BGI before assembling, as indicated previously [10]. 2.3. Phylogenetic Analyses and Community Comparisons After merging the forward and reverse reads, we obtained high-quality reads, and operational taxonomic unit (OTU) clustering and taxonomic assignment were performed as previously described [11,12]. In brief, each forward and reverse high-quality read for the paired-end library was merged using an USEARCH (version 8.0.1517). We discarded merged reads that (1) contained ambiguous nucleotides and (2) were mapped to the PhiX genome sequence by a Bowtie 2 (version 2.2.3). We obtained high-quality reads without forward and reverse primer sequences after removal of the reads that contained 450 nt and were associated with an average Phred-like quality score of less than 25, as calculated by the Illumina MiSeq sequencer. Sequence clustering of the high-quality reads was conducted using UCLUST (version 8.0.1517) with identity >97%, and query and reference coverage >80%. Chimeric OTUs were detected and removed if the OTUs were assigned to chimaera in both of the following methods: (1) UCHIME (version 6.0.307) reference mode search against the reference gold database (http://drive5.com/uchime/gold.fa) and (2) UCHIME de novo mode search. Taxonomic assignment of the high-quality reads was performed by RDP Classifier (version 2.12). Clustering of data were made with R (3.2.2) add-in commands (pvclust, prcomp, etc.). 2.4. Ethics All the sampling and access to the Japanese rock ptarmigans in this study was permitted by the Ministry of the Environment as the approved study for the Environment Research and Technology Development Fund (4-1604), the Agency for Cultural Affairs, and the Ministry of Education, Science, Sport and Culture (27-4-365), and sampling in a protected area was permitted by Chubu Regional Forest Office, Ministry of Agriculture, Forestry and Fishery (28-Toyama-113). 3. Result and Discussions A total of 645,636 clean reads were obtained, and each sample had 25,825 ± 12,593 reads on average. In the family level taxonomic assignment, 168 families were detected, but we filtered families of which abundance was below 0.01% in each sample. After that, we filtered families which were detected in a sole sample or only in two samples, and we finally obtained 26 families. In the genus level taxonomic assignment, 476 genera were detected, but we filtered genera of which abundance was below 0.01% in each sample. After that, we filtered genera which were detected in a sole sample or only in two samples, and we finally obtained 66 genera. Community structure at the family and the genus levels are shown in Figure 3a,b. At the family level, 26 unique bacterial families were analyzed. The top 10 abundant families include Coriobacteriaceae, Lachnospiraceae, Veillonellaceae, Actinomycetaceae, Ruminococcaceae, Synergistaceae, Bifidobacteriaceae, Enterobacteriaceae, Prevotellaceae, and Erysipelotrichaceae (Figure 3a). These families occupied about 97% of total reads classified at the family level in any sample. Coriobacteraceae and Lachnospraceae were the two dominant bacterial families in all samples analyzed. These two families were followed by Veillonelaceae and Actinomycetales in most cases. The other abundant families were detected in varied proportions, and prevalence of Enterobacteriaceae in three samples of Mt. Otensho was higher than others (Figure 3a).

Microorganisms 2018, 6, 77

5 of 11

At the genus level analysis, 66 unique genera were analyzed. Among them, 10 genera, Olsenella, Actinomyces, Megasphaera, Slackia, Cloacibacillus, Bifidobacterium, Escherichia, Dialister, Megamonas, and Bilophila were dominant (Figure 3b). The sum of reads from these genera explains 85% of total reads identified at the genus level in any sample. These 10 abundant genera were followed by 56 other genera, and the sum of the top 66 genera occupied 99% of total reads classified at the genus level. However, distribution of these abundant genera was not necessarily uniform in all samples collected. Anaerobiospirillum and Gardnerella were detected only in three samples, while the former was detected only in samples from Tateyama (Tate11_2, Tate11_3, Tate2) and the latter was detected only in samples Microorganisms 2018, 6, x FOR PEER REVIEW 6 of 12 from Mt. Norikura (Norikura2. Norikura5) and Mt. Kita (Kita2). 100%

90%

80%

70%

60%

50%

40%

30%

20%

10%

Jo n

en 1 Ki ta 1 Ki ta 2 K No ita3 rik No ura rik 1 No ura rik 2 No ura3 rik No ura4 rik ur a5 On 1 On 2 On 3 Ot On4 en s Ot ho1 en s Ot ho2 en sh o3 Ta Ta te1 te 1 Ta 1_1 te 1 Ta 1_2 te 11 _3 Ta te 2 Ta te 3 Ta te 4 Ta te Ya 5 ke N1

0%

Coriobacteriaceae

UnClassified

Lachnospiraceae

Veillonellaceae

Ac nomycetaceae

Ruminococcaceae

Synergistaceae

Bifidobacteriaceae

Enterobacteriaceae

Desulfovibrionaceae

Prevotellaceae

Erysipelotrichaceae

Dermabacteraceae

Porphyromonadaceae

Bacteroidaceae

Spirochaetaceae

Succinivibrionaceae

Anaeroplasmataceae

Helicobacteraceae

Clostridiaceae 4

Bdellovibrionaceae

Bacteroidales_incertae_sedis

Su erellaceae

Clostridiales_Incertae Sedis XIII

Streptococcaceae

Oxalobacteraceae

(a) Figure 3. Cont.

Microorganisms 2018, 6, 77 Microorganisms 2018, 6, x FOR PEER REVIEW

6 of 11

7 of 12

100% 90% 80% 70% 60% 50% 40% 30% 20% 10%

O Ot n4 en sh Ot o1 en sh Ot o2 en sh o3 Ta t Ta e1 te 11 Ta _1 te 11 Ta _2 te 11 _3 Ta te 2 Ta te 3 Ta te 4 Ta te 5 Ya ke N1

On 3

On 2

On 1

No 3 rik u No r a 1 rik ur No a 2 rik u No r a 3 rik u No ra4 rik ur a5

1

2

Ki ta

Ki ta

Ki ta

Jo

ne n1

0%

UnClassified

Olsenella

Ac nomyces

Megasphaera

Slackia

Cloacibacillus

Bifidobacterium

Escherichia/Shigella

Dialister

Megamonas

Bilophila

Ruminococcus2

Shu leworthia

Pseudoflavonifractor

Erysipelotrichaceae_incertae_sedis

Lachnospiracea_incertae_sedis

Flavonifractor

Parvibacter

Anaeroglobus

Alloprevotella

Dermabacter

Bacteroides

Intes nimonas

Asteroleplasma

Helicobacter

Allisonella

Paraprevotella

Treponema

Clostridium IV

Succina monas

Gardnerella

Anaerobiospirillum

Sporobacter

Clostridium XlVa

Coprococcus

Robinsoniella

Prevotella

Aceta factor

Atopobium

Clostridium XlVb

Vampirovibrio

Jonquetella

Acetobacteroides

Sphaerochaeta

Coprobacter

Syntrophococcus

Phocaeicola

Pseudocitrobacter

Anaeroplasma

Papillibacter

Clostridium XVIII

Collinsella

Oscillibacter

Anaerosporobacter

Blau a

Ace tomaculum

Wolinella

Mogibacterium

Anaerotruncus

Synergistes

Butyricicoccus

Dorea

Tannerella

Lactococcus

Oxalobacter

Eggerthella

(b)

Figure3.3. (a) (a) Community Communitystructure structure family level of cecal microbiota of wild Japanese Figure at at thethe family level of cecal microbiota of wild Japanese rock ptarmigans at different locations in Japanesein mountains. Community structure at the genus level of rock ptarmigans at different locations Japanese(b) mountains. (b) Community structure at cecal microbiota Japanese rock ptarmigans at different locations in Japanese mountains. the genus levelofofwild cecal microbiota of wild Japanese rock ptarmigans at different locations in Japanese mountains. At the OTU level analysis, the top 10 abundant OTUs occupied more than 50% of total reads and At the OTU level analysis, the more top 10than abundant moresample. than 50% of total and top 100 abundant OTUs occupied 75% ofOTUs totaloccupied reads in any These topreads 10 OTUs top 100 abundant OTUs occupied more than 75% of total reads in any sample. These top 10 OTUs were identified as Olsenella sp., Shuttleworthia sp., Actinomyces sp., Bifidobacterium sp., Megasphaera were identified as Alkalibaculum Olsenella sp., sp., Shuttleworthia Actinomyces sp.,(these Megasphaera sp., two species of Slackia sp.,sp., Robinsoniella sp.,sp., andBifidobacterium Paraprevotella sp. results sp., two species of Alkalibaculum sp., Slackia sp., Robinsoniella sp., and Paraprevotella sp. (these results were also confirmed by BLAST). Identification levels of these OTUs were not always high enough. were also confirmed by BLAST). Identification levels of these OTUs were not always high enough. For example, Confidence values calculated by RDP Classifer for OTUs identified as Shuttleworthia sp., For example, Confidence valuessp., calculated RDP Classifer OTUs identified assp. Shuttleworthia sp., Alkalibaculum sp., Robinsoniella anotherby Alkalibaculum sp.for and Parabacteroides were 0.38, 0.18, Alkalibaculum sp., Robinsoniella sp., another Alkalibaculum sp. and Parabacteroides sp. were 0.38, 0.18, 0.23, 0.21, and 0.32, respectively, while the other OTUs showed high confidence values between 0.82 0.23,0.96 0.21,(Table and 0.32, and 1). respectively, while the other OTUs showed high confidence values between 0.82 and 0.96 (Table 1).

Microorganisms 2018, 6, 77

7 of 11

Table 1. Top 10 abundant OTUs in cecal microbiota of wild Japanese rock ptarmigans at different locations in Japanese mountains. OTU Code

Domain

Phylum

Class

Order

Family

Genus

FCAYKGM7

Bacteria

1

Actinobacteria

1

Actinobacteria

1

Coriobacteriales

1

Coriobacteriaceae

1

Olsenella

0.96

FCAYKGM14

Bacteria

1

Firmicutes

1

Clostridia

1

Clostridiales

1

Lachnospiraceae

1

Shuttleworthia

0.38

FCAYKGM12

Bacteria

1

Actinobacteria

1

Actinobacteria

1

Actinomycetales

1

Actinomycetaceae

0.92

Actinomyces

0.9

FCAYKGM22

Bacteria

1

Firmicutes

1

Negativicutes

1

Selenomonadales

1

Veillonellaceae

1

Megasphaera

0.82

FCAYKGM20

Bacteria

1

Firmicutes

0.96

Clostridia

0.87

Clostridiales

0.86

Eubacteriaceae

0.2

Alkalibaculum

0.18

FCAYKGM29

Bacteria

1

Actinobacteria

0.96

Actinobacteria

0.96

Coriobacteriales

0.94

Coriobacteriaceae

0.94

Slackia

0.67

FCAYKGM30

Bacteria

1

Firmicutes

0.99

Clostridia

0.98

Clostridiales

0.98

Lachnospiraceae

0.44

Robinsoniella

0.23

FCAYKGM50

Bacteria

1

Actinobacteria

1

Actinobacteria

1

Bifidobacteriales

1

Bifidobacteriaceae

1

Bifidobacterium

0.91

FCAYKGM90

Bacteria

1

Firmicutes

0.91

Clostridia

0.87

Clostridiales

0.86

Eubacteriaceae

0.22

Alkalibaculum

0.21

FCAYKGM67

Bacteria

1

Bacteroidetes

1

Bacteroidia

1

Bacteroidales

1

Prevotellaceae

0.39

Paraprevotella

0.32

Microorganisms 2018, 6, 77

8 of 11

A cluster dendrogram at the genus level shows the separation of “Norikura3, YakeN1, Norikura2, Kita1, Norikura5, Kita3, Otensho1, and Tate11_1” from others. “Otensho2, and Otensho3” are further separated from the other two clusters, each composed of “On4, Kita2, Norikura1, Norikura4, On3, and, Tate4”, or “Jonen1, Tate11_3, Tate1, On1, Tate11_2, Tate2, Tate3, and Tate5”. The separation of these clusters is clearly shown in a high AU (Approximately Unbiased) p-value >90% (Figure 4a). The samples collected from same area show the relative similarities; for example, “Otensho2 and Otensho3”, “Tate2, Tate3, and Tate5”, “Norikura1 and 4”. However, several samples from geographically distant areas were close to each other, such as “Norikura2, Otensho1, Kita1, and Norikura5”. The cluster dendrogram at the family level shows the same tendencies as the genus level analysis (Figure 4b). We thought, prior to this study, that there would be the differences in the cecal bacterial Microorganisms 2018, 6, x FOR PEER REVIEW 9 of 12 communities of wild individuals located in different geographical areas, due to the differences in plant coverage. However, it is still hard find the systematic difference in cecal microbial composition of A cluster dendrogram at tothe genus level shows the separation of “Norikura3, YakeN1, wild individuals living in different locations, because major bacterial components were mostly present Norikura2, Kita1, Norikura5, Kita3, Otensho1, and Tate11_1” from others. “Otensho2, and Otensho3” in are all samples hierarchical clustering showedeach no systematic area. The main further analyzed, separatedand from the other two clusters, composedseparation of “On4,byKita2, Norikura1, components of cecal microbiota were the same in all birds tested. This means that the essential bacteria Norikura4, On3, and, Tate4”, or “Jonen1, Tate11_3, Tate1, On1, Tate11_2, Tate2, Tate3, and Tate5”. forThe the separation life of wild of Japanese rock ptarmigans the same. these clusters is clearly are shown in a high AU (Approximately Unbiased) p-value Among the top 10 abundant OTUs, which were identified as Olsenella sp., Shuttleworthia sp., >90% (Figure 4a). Actinomyces sp., Bifidobacterium sp., Megasphaera two species of Alkalibaculum sp., Slackia sp., The samples collected from same area showsp., the relative similarities; for example, “Otensho2 and Robinsoniella sp. and Parabacteroides sp., four OTUs corresponding to two Alkalibaculum sp., Otensho3”, “Tate2, Tate3, and Tate5”, “Norikura1 and 4”. However, several samples from one Robinsoniella sp., andareas one Parabactroides were notsuch identified with highOtensho1, confidenceKita1, values, geographically distant were close tosp., each other, as “Norikura2, and even at the family is therefore at suggested that these OTUs betendencies in a previously Norikura5”. The level. clusterItdendrogram the family level shows theare same as theunknown genus level family, in the order4b). Clostridiales in the case of former three OTUs, and in the order Bacteroidales in the analysis (Figure case of latter one OTU. Cluster dendrogram with AU/BP values (%)

0.3

0.4

au bp edge #

0.2

92 16 23

On2 Distance: correlation Cluster method: complete

Figure 4. Cont.

Tate5

Tate2

98 37 14 88 45 1

Tate3

Tate11_2

On1

Tate1

Tate11_3

0 0 18 47 10 72 17 13 48 10 6 84 27 5 3

Jonen1

On3

87 14 12

Tate4

Norikura4

Kita2

36 15 80 16 11 10 99 33 8

Norikura1

Otensho3

89 22 19

0 0 17

87 68 16

Otensho2

Norikura5

Kita1

Norikura2

Otensho1

Tate11_1

YakeN1

97 13 15 3 23 1 19 75 44 9 7 4 89 35 2

Kita3

0.0

Norikura3

0.1

55 24 21 49 10 20

On4

Height

0 0 22

(a)

Microorganisms 2018, 6, 776, x FOR PEER REVIEW Microorganisms 2018,

9 of 1011 of 12

Cluster dendrogram with AU/BP values (%)

0.6

au bp edge #

0.4

0.5

45 10 23

0.3

98 12 20

81 15 17 93 18

Tate5

On4

98 24 8 100 67 1

Tate3

Norikura5

Tate1

Kita1

12 100 25 7 46 40 4 95 54 2

Tate4

Tate2

97 40 11

Tate11_3

Otensho3

99 48 15

Otensho2

97 40 14

Otensho1

On3

Norikura4

On1

Jonen1

Kita2

Tate11_2

93 23 73 18 10 9 73 35 6 94 32 3

100 54 5

Norikura1

Norikura3

87 45 13

Kita3

90 15 16

On2

50 10 18

Norikura2

YakeN1

3 21 19

Tate11_1

0.2 0.0

0.1

Height

51 10 22

45 10 21

(b) Distance: correlation Cluster method: complete

Figure 4. (a) Cluster dendrogram at the genus level of cecal microbiota of wild Japanese

Figure 4. (a) Cluster dendrogram at the genus level of cecal microbiota of wild Japanese rock ptarmigans rock ptarmigans at Japanese differentmountains. locations in(b)Japanese mountains. at (b)the Cluster at different locations in Cluster dendrogram family dendrogram level of cecal at the family level of cecal microbiota of wild Japanese rock ptarmigans at different locations microbiota of wild Japanese rock ptarmigans at different locations in Japanese mountains.

in Japanese mountains.

The other four OTUs, corresponding to Olsenella sp., Actinomyces sp., Megasphaera sp. and We thought, prior to this study, that there would be the differences in the cecal bacterial Bifidobacterium sp. showed high confidence values, and so these OTUs are probably previously communities of wild individuals located in different geographical areas, due to the differences in known species, and their functionality, at least in fermentation, is highly plausible. Olsenella sp. and plant coverage. However, it is still hard to find the systematic difference in cecal microbial Bifidobacterium sp. can both be broadly regarded as health-promoting lactic acid bacteria [13,14], composition of wild individuals living in different locations, because major bacterial components and the high abundance of these species may be related to the high abundance of Megasphaera sp., were mostly present in all samples analyzed, and hierarchical clustering showed no systematic which is previously known as a lactate and acetate utilizing butyrate producer [15]. In our previous separation by area. The main components of cecal microbiota were the same in all birds tested. This study, Megasphaera sp. is metabolically associated with lactate producers in the porcine and murine means that the essential bacteria for the life of wild Japanese rock ptarmigans are the same. ceca [16,17]. Therefore, the abundance of these three bacterial species suggests the existence of one of Among the top 10 abundant OTUs, which were identified as Olsenella sp., Shuttleworthia sp., the major pathways of cecal fermentation in wild Japanese rock ptarmigans, as shown in a porcine or Actinomyces sp., Bifidobacterium sp., Megasphaera sp., two species of Alkalibaculum sp., Slackia sp., a murine cecal model [16,17]. The conversion of lactate occurs rapidly [18]. This conversion seems Robinsoniella sp. and Parabacteroides sp., four OTUs corresponding to two Alkalibaculum sp., one to be important in terms of energy metabolism in the Japanese rock ptarmigans, because lactate Robinsoniella sp., and one Parabactroides sp., were not identified with high confidence values, even at cannot be absorbed rapidly through gut epithelia compared with acetate, propionate and butyrate [19]. the family level. It is therefore suggested that these OTUs are be in a previously unknown family, in If lactate is accumulated in the ceca because of the lack of the conversion of lactate to SCFA (short the order Clostridiales in the case of former three OTUs, and in the order Bacteroidales in the case of chain fatty acid), Japanese rock ptarmigans cannot retrieve the energy effectively from their food. latter one OTU. In addition, these three bacterial genera were not found as major bacterial taxa in captive Svarbard rock The other four OTUs, corresponding to Olsenella sp., Actinomyces sp., Megasphaera sp. and ptarmigans [4]. The differences in microbial communities between captive Svarbard rock ptarmigans Bifidobacterium sp. showed high confidence values, and so these OTUs are probably previously and Japanese rock ptarmigans might be affected by genetic factors, but we speculated that the majority known species, and their functionality, at least in fermentation, is highly plausible. Olsenella sp. and of the differences come from differences between wild and captive habitats, lifestyles and diets. Bifidobacterium sp. can both be broadly regarded as health-promoting lactic acid bacteria [13,14], and These potentially essential bacteria for wild rock ptarmigans may have been removed or replaced by the high abundance of these species may be related to the high abundance of Megasphaera sp., which other bacterial taxa under captive conditions that would be the same as the ex-situ protection program is previously known as a lactate and acetate utilizing butyrate producer [15]. In our previous study, of Japanese rock ptarmigans. Therefore, the reconstitution of this high level of coexistence between Megasphaera sp. is metabolically associated with lactate producers in the porcine and murine ceca lactic acid bacteria and lactate-utilizing bacteria in cecal microbiota should be taken into account for [16,17]. Therefore, the abundance of these three bacterial species suggests the existence of one of the the successful breeding of captive Japanese rock ptarmigans in the national conservation program. major pathways of cecal fermentation in wild Japanese rock ptarmigans, as shown in a porcine or a murine cecal model [16,17]. The conversion of lactate occurs rapidly [18]. This conversion seems to be important in terms of energy metabolism in the Japanese rock ptarmigans, because lactate cannot

Microorganisms 2018, 6, 77

10 of 11

Supplementary Materials: The following are available online at http://www.mdpi.com/2076-2607/6/3/77/s1. Author Contributions: Conceptualization, Koichi Murata (K.M.) and Kazunari Ushida (K.U.); Data curation, Atsushi Ueda (A.U.), Takuji Yamada (T.Y.) and K.U.; Formal analysis, A.U., T.Y. and K.U.; Funding acquisition, K.U.; Investigation, A.U., Atsushi Kobayashi (A.K.), Sayaka Tsuchida (S.T.) and K.U.; Methodology, A.U., T.Y., A.K., S.T. and K.U.; Resources, A.K., S.T., K.M., Hiroshi Nakamura (H.N.) and K.U.; Software, A.U. and T.Y.; Supervision, T.Y., K.M., H.N. and K.U.; Validation, T.Y., K.M. and H.N.; Writing – original draft, A.U.; Writing – review & editing, A.U., A.K., S.T., T.Y., K.M., H.N. and K.U. Funding: This research was funded by the Ministry of Environment as the Environment Research and Technology Development Fund (ERTDF) 4-1604. Acknowledgments: The support from Tsutomu Matsuda, executive director of Toyama Raicho Kenkyukai (Toyama Ptarmigan Association), Raityo-So lodge and Kitadake Sanso Lodge are acknowledged. Conflicts of Interest: The authors declare no conflicts of interest.

References 1. 2. 3.

4. 5.

6.

7.

8. 9. 10.

11.

12.

13. 14.

15.

Suzuki, A.; Kobayashi, A.; Nakamura, H.; Takasu, F. Population viability analysis of the Japanese rock ptarmigan Lagopus muta japonica in Japan. Wildl. Biol. 2013, 19, 339–346. [CrossRef] Nakamura, H. Rock Ptarmigan Lagopus mutus japonicus. Jpn. J. Ornithol. 2007, 56, 93–114. [CrossRef] Ministry of the Environment, Government of Japan 2012. National Conservation Program for Japanese Rock Ptarmigans. Available online: https://www.env.go.jp/nature/kisho/hogozoushoku/pdf/jigyoukeikaku/ raicho.pdf (accessed on 25 July 2018). Ushida, K.; Segawa, T.; Tsuchida, S.; Murata, K. Cecal bacterial communities in wild Japanese rock ptarmigans and captive Svalbard rock ptarmigans. J. Vet. Med. Sci. 2016, 78, 251–257. [CrossRef] [PubMed] Tsuchida, S.; Murata, K.; Ohkuma, M.; Ushida, K. Isolation of Streptococcus gallolyticus with very high degradability of condensed tannins from feces of the wild Japanese rock ptarmigans on Mt. Tateyama. J. Gen. Appl. Microbiol. 2017, 63, 195–198. [CrossRef] [PubMed] Tuschida, S.; Ohara, Y.; Kuramochi, K.; Murata, K.; Ushida, K. Effective degradation of phenolic glycoside rhododendrin and its aglycone ghododendrol by cecal feces of wild Japanese rock ptarmigans. Jpn. J. Zoo Wildl. Med. 2017, 22, 41–45. Hird, S.M.; Carstens, B.C.; Cardiff, S.W.; Dittmann, D.L.; Brumfield, R.T. Sampling locality is more detectable than taxonomy or ecology in the gut microbiota of the brood-parasitic Brown-headed Cowbird (Molothrus ater). PeerJ 2014, 2, e321. [CrossRef] [PubMed] Baba, Y.; Fujimaki, Y.; Yoshii, R.; Koike, H. Genetic variability in the mitochondrial control region of the Japanese Rock Ptarmigan Lagopus mutus japonicas. Jpn. J. Ornithol. 2001, 50, 53–64. [CrossRef] Kobayashi, A.; Nakamura, H. Seasonal change of food items of the Japanese rock ptarmigan. Jpn. J. Ornithol. 2011, 60, 200–215. [CrossRef] Li, T.H.; Qin, Y.; Sham, P.C.; Lau, K.S.; Chu, K.M.; Leung, W.K. Alterations in Gastric Microbiota After H. Pylori Eradication and in Different Histological Stages of Gastric Carcinogenesis. Sci. Rep. 2017, 7, 44935. [CrossRef] [PubMed] Okai, S.; Usui, F.; Yokota, S.; Hori-I, Y.; Hasegawa, M.; Nakamura, T.; Kurosawa, M.; Okada, S.; Yamamoto, K.; Nishiyama, E.; et al. High-affinity monoclonal IgA regulates gut microbiota and prevents colitis in mice. Nat. Microbiol. 2016, 1, 16103. [CrossRef] [PubMed] Watanabe, H.; Goto, S.; Mori, H.; Higashi, K.; Hosomichi, K.; Aizawa, N.; Takahashi, N.; Tsuchida, M.; Suzuki, Y.; Yamada, T.; et al. Comprehensive microbiome analysis of tonsillar crypts in IgA nephropathy. Nephrol. Dial. Transplant. 2016, 32, 2072–2079. [CrossRef] [PubMed] Orla-Jensen, S. La classification des bactéries lactiques. Le Lait 1924, 4, 468–474. [CrossRef] Dewhirst, F.E.; Paster, B.J.; Tzellas, N.; Coleman, B.; Downes, J.; Spratt, D.A.; Wade, W.G. Characterization of novel human oral isolates and cloned 16S rDNA sequences that fall in the family Coriobacteriaceae: Description of Olsenella gen. nov., reclassification of Lactobacillus uli as Olsenella uli comb. nov. and description of Olsenella profusa sp. nov. Int. J. Syst. Evol. Microbiol. 2001, 51, 1797–1804. [PubMed] Gutierrez, J.; Davis, R.E.; Lindahl, I.L.; Warwick, E.J. Bacterial Changes in the Rumen During the Onset of Feed-lot Bloat of Cattle and Characteristics of Peptostreptococcus elsdenii n. sp. Appl. Microbiol. 1959, 7, 16–22. [PubMed]

Microorganisms 2018, 6, 77

16.

17.

18. 19.

11 of 11

Hashizume, K.; Tsukahara, T.; Yamada, K.; Koyama, H.; Ushida, K. Megasphaera elsdenii JCM1772T normalizes hyperlactate production in the large intestine of fructooligosaccharide-fed rats by stimulating butyrate production. J. Nutr. 2003, 133, 3187–3190. [CrossRef] [PubMed] Tsukahara, T.; Koyama, H.; Okada, M.; Ushida, K. Stimulation of butyrate production by gluconic acid in batch culture of pig cecal digesta and identification of butyrate-producing bacteria. J. Nutr. 2002, 132, 2229–2234. [CrossRef] [PubMed] Ushida, K.; Hoshi, S.; Ajisaka, K. 13 C-NMR studies on lactate metabolism in a porcine gut microbial ecosystem. Microb. Ecol. Health Dis. 2002, 14, 241–246. Umesaki, Y.; Yajima, T.; Yokokura, T.; Mutai, M. Effect of organic acid absorption on bicarbonate transport in rat colon. Pflüg. Arch. 1979, 379, 43–47. [CrossRef] © 2018 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).