BACKGROUND: Many arthropods rely on their gut microbiome to digest plant material, which is often low in nitrogen but high in complex polysaccharides. Detritivores, such as millipedes, live on a particularly poor diet, but the identity and nutritional contribution of their microbiome are largely unknown. In this study, the hindgut microbiota of the tropical millipede Epibolus pulchripes (large, methane emitting) and the temperate millipede Glomeris connexa (small, non-methane emitting), fed on an identical diet, were studied using comparative metagenomics and metatranscriptomics. RESULTS: The results showed that the microbial load in E. pulchripes is much higher and more diverse than in G. connexa. The microbial communities of the two species differed significantly, with Bacteroidota dominating the hindguts of E. pulchripes and Proteobacteria (Pseudomonadota) in G. connexa. Despite equal sequencing effort, de novo assembly and binning recovered 282 metagenome-assembled genomes (MAGs) from E. pulchripes and 33 from G. connexa, including 90 novel bacterial taxa (81 in E. pulchripes and 9 in G. connexa). However, despite this taxonomic divergence, most of the functions, including carbohydrate hydrolysis, sulfate reduction, and nitrogen cycling, were common to the two species. Members of the Bacteroidota (Bacteroidetes) were the primary agents of complex carbon degradation in E. pulchripes, while members of Proteobacteria dominated in G. connexa. Members of Desulfobacterota were the potential sulfate-reducing bacteria in E. pulchripes. The capacity for dissimilatory nitrate reduction was found in Actinobacteriota (E. pulchripes) and Proteobacteria (both species), but only Proteobacteria possessed the capacity for denitrification (both species). In contrast, some functions were only found in E. pulchripes. These include reductive acetogenesis, found in members of Desulfobacterota and Firmicutes (Bacillota) in E. pulchripes. Also, diazotrophs were only found in E. pulchripes, with a few members of the Firmicutes and Proteobacteria expressing the nifH gene. Interestingly, fungal-cell-wall-degrading glycoside hydrolases (GHs) were among the most abundant carbohydrate-active enzymes (CAZymes) expressed in both millipede species, suggesting that fungal biomass plays an important role in the millipede diet. CONCLUSIONS: Overall, these results provide detailed insights into the genomic capabilities of the microbial community in the hindgut of millipedes and shed light on the ecophysiology of these essential detritivores. Video Abstract.

• Keywords
• Acetogens, Ecosystem engineers, Glycoside hydrolases, Hindgut microbiota, Millipede holobiont, Nutrient cycling, Polysaccharide degradation, Symbiosis,
• MeSH
• Bacteria MeSH
• Bacteroidetes genetics   MeSH
• Arthropods * genetics   MeSH
• Nitrogen metabolism   MeSH
• Phylogeny MeSH
• Metagenome MeSH
• Metagenomics MeSH
• Proteobacteria genetics   MeSH
• Carbohydrates MeSH
• Sulfates metabolism   MeSH
• Gastrointestinal Microbiome * genetics   MeSH
• Animals MeSH
• Check Tag
• Animals MeSH
• Publication type
• Video-Audio Media MeSH
• Journal Article MeSH
• Names of Substances
• Nitrogen MeSH
• Carbohydrates MeSH
• Sulfates MeSH

Amplified ribosomal DNA restriction analysis (ARDRA) was used to compare the bacterial communities of the food, the gut sections (ceca, anterior and posterior midgut, hindgut) and the excrement of the litter feeding bibionid larvae of Penthetria holosericea. For universal eubacterial primers ARDRA patterns were complex with only minor differences among samples. Taxon specific primers were also applied to characterize the samples. Fragment composition was transformed to presence/absence binary data and further analyzed. Cluster analysis revealed that bacterial communities of gut highly resembled each other with the exception of the ceca. ARDRA patterns of consumed leaves clustered together with the intact leaves but differed from those of the excrement. ARDRA results were compared with microbial community structure based on phospholipid fatty acid (PLFA) fingerprints. The cluster analysis of PLFA (presence/absence binary) data resulted in a pattern similar to the ARDRA data. The PCA analysis of PLFA relative content separated microbial communities into five groups: (1) anterior and posterior midgut, (2) hindgut, (3) ceca, (4) consumed and intact litter, (5) excrement. Both methods indicated that conditions in the larval gut result in formation of a specific microbial community which differs from both the food and excrement ones. Particularly ceca--(blind appendages, harbor very specific microbial community) are divided from the rest of the gut by perithropic membrane.

• MeSH
• Bacteria chemistry   genetics   isolation & purification   MeSH
• Cecum microbiology   MeSH
• Diptera microbiology   MeSH
• DNA Fingerprinting methods   MeSH
• Feces microbiology   MeSH
• Phospholipids analysis   MeSH
• Animal Feed microbiology   MeSH
• Larva microbiology   MeSH
• Fatty Acids analysis   MeSH
• Polymorphism, Restriction Fragment Length MeSH
• Cluster Analysis MeSH
• Intestines microbiology   MeSH
• Animals MeSH
• Check Tag
• Animals MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Comparative Study MeSH
• Names of Substances
• Phospholipids MeSH
• Fatty Acids MeSH