Bacterial diversity in the waterholes of the Kruger National Park: an eDNA metabarcoding approach 1

Exploration and comparison of bacterial communities present in bovine faeces, milk and blood using 16S rRNA metagenomic sequencing

Cattle by-products like faeces, milk and blood have many uses among rural communities; aiding to facilitate everyday household activities and occasional rituals. Ecologically, the body sites from which they are derived consist of distinct microbial communities forming a complex ecosystem of niches. We aimed to explore and compare the faecal, milk and blood microbiota of cows through 16S rRNA sequencing. All downstream analyses were performed using applications in R Studio (v3.6.1). Alpha-diversity metrics showed significant differences between faeces and blood; faeces and milk; but non-significant between blood and milk using Kruskal-Wallis test, P < 0,05. The beta-diversity metrics on Principal Coordinate Analysis and Non-Metric Dimensional Scaling significantly clustered samples by type (PERMANOVA test, P < 0,05). The overall analysis revealed a total of 30 phyla, 74 classes, 156 orders, 243 families and 408 genera. Firmicutes, Bacteroidota and Proteobacteria were the most abundant phyla overall. A total of 58 genus-level taxa occurred concurrently between the body sites. The important taxa could be categorized into four potentially pathogenic clusters i.e. arthropod-borne; food-borne and zoonotic; mastitogenic; and metritic and abortigenic. A number of taxa were significantly differentially abundant (DA) between sites based on the Wald test implemented in DESeq2 package. Majority of the DA taxa (i.e. Romboutsia, Paeniclostridium, Monoglobus, Akkermansia, Turicibacter, Bacteroides, Candidatus_Saccharimonas, UCG-005 and Prevotellaceae_UCG-004) were significantly enriched in faeces in comparison to milk and blood, except for Anaplasma which was greatly enriched in blood and was in turn the largest microbial genus in the entire analysis. This study provides insights into the microbial community composition of the sampled body sites and its extent of overlapping. It further highlights the potential risk of disease occurrence and transmission between the animals and the community of Waaihoek in KwaZulu-Natal, Republic of South Africa pertaining to their unsanitary practices associated with the use of cattle by-products.

Phylogenies of the 16S rRNA gene and its hypervariable regions lack concordance with core genome phylogenies

BACKGROUND

The 16S rRNA gene is used extensively in bacterial phylogenetics, in species delineation, and now widely in microbiome studies. However, the gene suffers from intragenomic heterogeneity, and reports of recombination and an unreliable phylogenetic signal are accumulating. Here, we compare core gene phylogenies to phylogenies constructed using core gene concatenations to estimate the strength of signal for the 16S rRNA gene, its hypervariable regions, and all core genes at the intra- and inter-genus levels. Specifically, we perform four intra-genus analyses (Clostridium, n = 65; Legionella, n = 47; Staphylococcus, n = 36; and Campylobacter, n = 17) and one inter-genus analysis [41 core genera of the human gut microbiome (31 families, 17 orders, and 12 classes), n = 82].

RESULTS

At both taxonomic levels, the 16S rRNA gene was recombinant and subject to horizontal gene transfer. At the intra-genus level, the gene showed one of the lowest levels of concordance with the core genome phylogeny (50.7% average). Concordance for hypervariable regions was lower still, with entropy masking providing little to no benefit. A major factor influencing concordance was SNP count, which showed a positive logarithmic association. Using this relationship, we determined that 690 ± 110 SNPs were required for 80% concordance (average 16S rRNA gene SNP count was 254). We also found a wide range in 16S-23S-5S rRNA operon copy number among genomes (1-27). At the inter-genus level, concordance for the whole 16S rRNA gene was markedly higher (73.8% - 10th out of 49 loci); however, the most concordant hypervariable regions (V4, V3-V4, and V1-V2) ranked in the third quartile (62.5 to 60.0%).

CONCLUSIONS

Ramifications of a poor phylogenetic performance for the 16S rRNA gene are far reaching. For example, in addition to incorrect species/strain delineation and phylogenetic inference, it has the potential to confound community diversity metrics if phylogenetic information is incorporated - for example, with popular approaches such as Faith's phylogenetic diversity and UniFrac. Our results highlight the problematic nature of these approaches and their use (along with entropy masking) is discouraged. Lastly, the wide range in 16S rRNA gene copy number among genomes also has a strong potential to confound diversity metrics. Video Abstract.

Diversity and distribution of benthic invertebrates dwelling rivers of the Kruger National Park, South Africa

Meiobenthos (or meiofauna) are microscopic invertebrates that inhabit biofilms and interstitial spaces in rivers. They are diverse and extremely abundant, and they perform essential ecological functions by linking microbial production to higher trophic levels (e.g. macrobenthic invertebrates and fishes). However, meiobenthic communities remain poorly studied in Africa. Here, we sampled meio- and macrobenthic invertebrate communities associated with biofilms and sediments across an upstream–downstream gradient along the Olifants, Sabie and Crocodile rivers flowing through the Kruger National Park (KNP). We expected to link differences in community structure to environmental gradients as those rivers show different degrees of anthropogenic stress as they enter the park. Both meio- and macrobenthic communities differed across rivers and also structured along an upstream–downstream gradient. The upstream sites, which were the closest to the park borders, consistently showed a lower diversity in all three rivers. There, the invasive snail Tarebia granifera strongly dominated (making up 73% – 87% of the macrobenthos), crowding hard substrates, while concomitantly the abundances of biofilm-dwelling meiobenthos like nematodes and rotifers were substantially reduced. Nevertheless, the diversity and evenness of communities then tended to increase as water flowed downstream through the park, suggesting a beneficial effect of protected river reaches on benthic invertebrate diversity. However, for the Crocodile River, which makes up the southern border of the park, this trend was less conspicuous, suggesting that this river may experience the greatest threats. More generally, benthic invertebrate communities were driven by the concentrations of phosphates, sulphates, ammonium and organic matter and by substrate characteristics.Conservation implications: Meiobenthic organisms are very abundant in KNP rivers and react to environmental gradients; thus, they should be more considered for bio-monitoring or conservation of comprehensive assemblages of animals. Interestingly, protected reaches tended to show a reduced dominance of the invasive T. granifera and a higher diversity of benthic invertebrates.