A Guide to Dietary Pattern-Microbiome Data Integration

The human gut microbiome is linked to metabolic and cardiovascular disease risk. Dietary modulation of the human gut microbiome offers an attractive pathway to manipulate the microbiome to prevent microbiome-related disease. However, this promise has not been realized. The complex system of diet and microbiome interactions is poorly understood. Integrating observational human diet and microbiome data can help researchers and clinicians untangle the complex systems of interactions that predict how the microbiome will change in response to foods. The use of dietary patterns to assess diet-microbiome relations holds promise to identify interesting associations and result in findings that can directly translate into actionable dietary intake recommendations and eating plans. In this article, we first highlight the complexity inherent in both dietary and microbiome data and introduce the approaches generally used to explore diet and microbiome simultaneously in observational studies. Second, we review the food group and dietary pattern-microbiome literature focusing on dietary complexity-moving beyond nutrients. Our review identified a substantial and growing body of literature that explores links between the microbiome and dietary patterns. However, there was very little standardization of dietary collection and assessment methods across studies. The 54 studies identified in this review used ≥7 different methods to assess diet. Coupled with the variation in final dietary parameters calculated from dietary data (e.g., dietary indices, dietary patterns, food groups, etc.), few studies with shared methods and assessment techniques were available for comparison. Third, we highlight the similarities between dietary and microbiome data structures and present the possibility that multivariate and compositional methods, developed initially for microbiome data, could have utility when applied to dietary data. Finally, we summarize the current state of the art for diet-microbiome data integration and highlight ways dietary data could be paired with microbiome data in future studies to improve the detection of diet-microbiome signals.

Ulcerative Colitis Seems to Imply Oral Microbiome Dysbiosis

Ulcerative colitis (UC) is a recurrent pathology of complex etiology that has been occasionally associated with oral lesions, but the overall composition of the oral microbiome in UC patients and its role in the pathogenesis of the disease are still poorly understood. In this study, the oral microbiome of UC patients and healthy individuals was compared to ascertain the possible changes in the oral microbial communities associated with UC. For this, the salivary microbiota of 10 patients diagnosed with an active phase of UC and 11 healthy controls was analyzed by 16S rRNA gene sequencing (trial ref. ISRCTN39987). Metataxonomic analysis revealed a decrease in the alpha diversity and an imbalance in the relative proportions of some key members of the oral core microbiome in UC patients. Additionally, Staphylococcus members and four differential species or phylotypes were only present in UC patients, not being detected in healthy subjects. This study provides a global snapshot of the existence of oral dysbiosis associated with UC, and the possible presence of potential oral biomarkers.

Parasite communities in English Sole (Parophrys vetulus) have changed in composition but not richness in the Salish Sea, Washington, USA since 1930

Earth is rapidly losing free-living species. Is the same true for parasitic species? To reveal temporal trends in biodiversity, historical data are needed, but often such data do not exist for parasites. Here, parasite communities of the past were reconstructed by identifying parasites in fluid-preserved specimens held in natural history collections. Approximately 2500 macroparasites were counted from 109 English Sole (Parophrys vetulus ) collected between 1930 and 2019 in the Salish Sea, Washington, USA. Alpha and beta diversity were measured to determine if and how diversity changed over time. Species richness of parasite infracommunities and community dispersion did not vary over time, but community composition of decadal component communities varied significantly over the study period. Community dissimilarity also varied: prior to the mid-20th century, parasites shifted in abundance in a seemingly stochastic manner and, after this time period, a canalization of community change was observed, where species' abundances began to shift in consistent directions. Further work is needed to elucidate potential drivers of these changes and to determine if these patterns are present in the parasite communities of other fishes of the Salish Sea.

Gut Microbiota in Heart Failure Patients With Preserved Ejection Fraction (GUMPTION Study)

Introduction: Heart failure with preserved ejection fraction (HFpEF) is associated with disrupted intestinal epithelial function, resulting from intestinal congestion. Intestinal congestion changes the morphology and permeability of the intestinal wall, and it becomes easy for the gut microbiota to change and transfer. Intervention on gut microbiota may become a new target for HFpEF treatment. However, the characteristics of gut microbiota in patients with HFpEF remain unknown. This preliminary report aims to detect the structure of gut microbiota in HFpEF patients so as to explore their characteristic changes, thereby providing a theoretical basis for future research. Methods: This research recruited 30 patients diagnosed with HFpEF and 30 healthy individuals. Stool specimens of research subjects were collected separately, and the microarray analyses of gut microbiota were conducted by Illumina high-throughput DNA sequencing. The differences in gut microbiota composition, alpha diversity, and beta diversity between the two groups were finally obtained. Results: The composition of gut microbiota was significantly different between the two groups. At the phylum classification level, the abundance of Synergistetes tended to be higher in the HFpEF group (P = 0.012). At genus classification level, the abundance of Butyricicoccus (P < 0.001), Sutterella (P = 0.004), Lachnospira (P = 0.003), and Ruminiclostridium (P = 0.009) in the HFpEF group were lower, while the abundance of Enterococcus (P < 0.001) and Lactobacillus (P = 0.005) were higher. According to the Chao index of alpha diversity analysis, HFpEF patients showed a nominally significant lower species richness when compared with controls (P = 0.046). However, there was no statistical difference in the Shannon index (P = 0.159) and Simpson index (P = 0.495), indicating that there was no difference in species diversity between the two groups. Beta diversity analysis revealed a highly significant separation of HFpEF patients and controls. Conclusions: An imbalance in the gut microbiota of HFpEF patients was observed. Patients with HFpEF have an increased abundance of microbiota associated with inflammation and a decreased abundance of microbiota associated with anti-inflammatory effects in the gut environment. In line with that, the species richness of gut microbiota in HFpEF patients tended to be lower.

Contrasting Community Assembly Mechanisms Underlie Similar Biogeographic Patterns of Surface Microbiota in the Tropical North Pacific Ocean

Marine microbiota are critical components of global biogeochemical cycles. However, the biogeographic patterns and ecological processes that structure them remain poorly understood, especially in the oligotrophic ocean. In this study, we used high-throughput sequencing of 16S and 18S rRNA genes to investigate the distribution patterns of bacterial and microeukaryotic communities and their assembly mechanisms in the surface waters of the tropical North Pacific Ocean. The fact that both the bacterial and the microeukaryotic communities showed similar distribution patterns (i.e., similar distance-decay patterns) and were clustered according to their geographic origin (i.e., the western tropical North Pacific and central tropical North Pacific) suggested that there was a significant biogeographic pattern of microbiota in the North Pacific Ocean. Indices of alpha diversity such as species richness, phylogenetic diversity, and the Shannon diversity index also differed significantly between regions. The correlations were generally similar between spatial and environmental variables and the alpha and beta diversities of bacteria and microeukaryotes across the entire region. The relative importance of ecological processes differed between bacteria and microeukaryotes: ecological drift was the principal mechanism that accounted for the structure of bacterial communities; heterogeneous selection, dispersal limitation, and ecological drift collectively explained much of the turnover of the microeukaryote communities.

IMPORTANCE Bacteria and microeukaryotes are extremely diverse groups in the ocean, where they regulate elemental cycling and energy flow. Studies of marine microbial ecology have benefited greatly from the rapid progress that has been made in genomic sequencing and theoretical microbial ecology. However, the spatial distribution of marine bacteria and microeukaryotes and the nature of the assembly mechanisms that determine their distribution patterns in oligotrophic marine waters are poorly understood. In this study, we used high-throughput sequencing methods to identify the distribution patterns and ecological processes of bacteria and microeukaryotes in an oligotrophic, tropical ocean. Our study showed that contrasting community assembly mechanisms underlaid similar biogeographic patterns of surface bacterial and microeukaryotic communities in the tropical North Pacific Ocean.

Multigene Phylogeny, Diversity and Antimicrobial Potential of Endophytic Sordariomycetes From Rosa roxburghii

Rosa roxburghii Tratt. is widely applied in food, cosmetics, and traditional medicine, and has been demonstrated to possess diverse bioactivities. Plant endophytic fungi are important microbial resources with great potential for application in many fields. They not only establish mutualistic symbiosis with host plants but also produce a variety of bioactive compounds. Therefore, in the present study, endophytic fungi were isolated from R. roxburghii, the diversity and antimicrobial activities were evaluated. As a result, 242 strains of endophytic Sordariomycetes were successfully isolated. Multigene phylogenetic analyses showed that these isolates included eight orders, 19 families, 33 genera. The dominant genera were Diaporthe (31.4%), Fusarium (14.4%), Chaetomium (7.9%), Dactylonectria (7.0%), Graphium (4.5%), Colletotrichum (4.1%), and Clonostachys (4.1%). For different tissues of R. roxburghii, alpha diversity analysis revealed that the diversity of fungal communities decreased in the order of root, fruit, stem, flower, leaf, and seed, and Clonostachys and Dactylonectria exhibited obvious tissue specificity. Meanwhile, functional annotation of 33 genera indicated that some fungi have multitrophic lifestyles combining endophytic, pathogenic, and saprophytic behavior. Additionally, antimicrobial activities of endophytic Sordariomycetes against Lasiodiplodia theobromae, Botryosphaeria dothidea, Colletotrichum capsici, Pyricularia oryzae, Rhizoctonia solani, Fusarium oxysporum, Pseudomonas syringae, Pantoea agglomerans, Staphylococcus aureus, Bacillus subtilis, Escherichia coli, and Pseudomonas aeruginosa were screened. Dual culture test assays showed that there were 40 different endophytic species with strong inhibition of at least one or moderate inhibition of two or more against the 12 tested strains. The results from the filter paper diffusion method suggested that extracellular metabolites may be more advantageous than intracellular metabolites in the development of antimicrobial agents. Eleven isolates with good activities were screened. In particular, Hypomontagnella monticulosa HGUP194009 and Nigrospora sphaerica HGUP191020 have shown promise in both broad-spectrum and intensity. Finally, some fungi that commonly cause disease have been observed to have beneficial biological activities as endophytic fungi. In conclusion, this study showed the species composition, alpha diversity, and lifestyle diversity of endophytic Sordariomycetes from R. roxburghii and demonstrated these isolates are potential sources for exploring antimicrobial agents.

Application of Clustering Method to Explore the Correlation Between Dominant Flora and the Autism Spectrum Disorder Clinical Phenotype in Chinese Children

Autism spectrum disorder (ASD) is characterized by deficits in social interactions and repetitive, stereotypic behaviors. Evidence shows that bidirectional communication of the gut-brain axis plays an important role. Here, we recruited 62 patients with ASD in southern China, and performed a cross-sectional study to test the relationship between repeated behavior, gut microbiome composition, and alpha diversity. We divided all participants into two groups based on the clustering results of their microbial compositions and found Veillonella and Ruminococcus as the seed genera in each group. Repetitive behavior differed between clusters, and cluster 2 had milder repetitive symptoms than Cluster 1. Alpha diversity between clusters was significantly different, indicating that cluster 1 had lower alpha diversity and more severe repetitive, stereotypic behaviors. Repetitive behavior had a negative correlation with alpha diversity. We demonstrated that the difference in intestinal microbiome composition and altered alpha diversity can be associated with repetitive, stereotypic behavior in autism. The role of Ruminococcus and Veillonella in ASD is not yet understood.

Advanced Dental Cleaning is Associated with Reduced Risk of COPD Exacerbations - A Randomized Controlled Trial

Purpose

Infections from the oral microbiome may lead to exacerbations of chronic obstructive pulmonary disease (COPD). We investigated whether advanced dental cleaning could reduce exacerbation frequency. Secondary outcomes were disease-specific health status, lung function, and whether the bacterial load and composition of plaque microbiome at baseline were associated with a difference in outcomes.

Patients and Methods

One-hundred-one primary and secondary care patients with COPD were randomized to intervention with advanced dental cleaning or to dental examination only, repeated after six months. At baseline and at 12 months, data of exacerbations, lung function, COPD Assessment Test (CAT) score, and periodontal status were collected from questionnaires, record review, and periodontal examination. Student’s t-test and Mann–Whitney-U (MWU) test compared changes in outcomes. The primary outcome variable was also assessed using multivariable linear regression with adjustment for potential confounders. Microbiome analyses of plaque samples taken at baseline were performed using Wilcoxon signed ranks tests for calculation of alpha diversity, per mutational multivariate analysis of variance for beta diversity, and receiver operating characteristic curves for prediction of outcomes based on machine learning models.

Results

In the MWU test, the annual exacerbation frequency was significantly reduced in patients previously experiencing frequent exacerbations (p = 0.020) and in those with repeated advanced dental cleaning (p = 0.039) compared with the non-treated control group, but not in the total population including both patients with a single and repeated visits (p = 0.207). The result was confirmed in multivariable linear regression, where the risk of new exacerbations was significantly lower in patients both in the intention to treat analysis (regression coefficient 0.36 (95% CI 0.25–0.52), p < 0.0001) and in the population with repeated dental cleaning (0.16 (0.10–0.27), p < 0.0001). The composition of microbiome at baseline was moderately predictive of an increased risk of worsened health status at 12 months (AUC = 0.723).

Conclusion

Advanced dental cleaning is associated with a reduced frequency of COPD exacerbations. Regular periodontal examination and dental cleaning may be of clinical importance to prevent COPD exacerbations.

Wildfire and topography drive woody plant diversity in a Sky Island mountain range in the Southwest USA

AIM

Drastic changes in fire regimes are altering plant communities, inspiring ecologists to better understand the relationship between fire and plant species diversity. We examined the impact of a 90,000-ha wildfire on woody plant species diversity in an arid mountain range in southern Arizona, USA. We tested recent fire-diversity hypotheses by addressing the impacts on diversity of fire severity, fire variability, historical fire regimes, and topography.

LOCATION

Chiricahua National Monument, Chiricahua Mountains, Arizona, USA, part of the Sky Islands of the US-Mexico borderlands.

TAXON

Woody plant species.

METHODS

We sampled woody plant diversity in 138 plots before (2002-2003) and after (2017-2018) the 2011 Horseshoe Two Fire in three vegetation types and across fire severity and topographic gradients. We calculated gamma, alpha, and beta diversity and examined changes over time in burned versus unburned plots and the shapes of the relationships of diversity with fire severity and topography.

RESULTS

Alpha species richness declined, and beta and gamma diversity increased in burned but not unburned plots. Fire-induced enhancement of gamma diversity was confined to low fire severity plots. Alpha diversity did not exhibit a clear continuous relationship with fire severity. Beta diversity was enhanced by variation in fire severity among plots and increased with fire severity up to very high severity, where it declined slightly.

MAIN CONCLUSIONS

The results reject the intermediate disturbance hypothesis for alpha diversity but weakly support it for gamma diversity. Spatial variation in fire severity promoted variation among plant assemblages, supporting the pyrodiversity hypothesis. Long-term drought probably amplified fire-driven diversity changes. Despite the apparent benign impact of the fire on diversity, the replacement of two large conifer species with a suite of drought-tolerant shrubs signals the potential loss of functional diversity, a pattern that may warrant restoration efforts to retain these important compositional elements.

Revisiting the effect of PCR replication and sequencing depth on biodiversity metrics in environmental DNA metabarcoding

Environmental DNA (eDNA) metabarcoding is an increasingly popular tool for measuring and cataloguing biodiversity. Because the environments and substrates in which DNA is preserved differ considerably, eDNA research often requires bespoke approaches to generating eDNA data. Here, we explore how two experimental choices in eDNA study design-the number of PCR replicates and the depth of sequencing of PCR replicates-influence the composition and consistency of taxa recovered from eDNA extracts. We perform 24 PCR replicates from each of six soil samples using two of the most common metabarcodes for Fungi and Viridiplantae (ITS1 and ITS2), and sequence each replicate to an average depth of ~84,000 reads. We find that PCR replicates are broadly consistent in composition and relative abundance of dominant taxa, but that low abundance taxa are often unique to one or a few PCR replicates. Taxa observed in one out of 24 PCR replicates make up 21-29% of the total taxa detected. We also observe that sequencing depth or rarefaction influences alpha diversity and beta diversity estimates. Read sampling depth influences local contribution to beta diversity, placement in ordinations, and beta dispersion in ordinations. Our results suggest that, because common taxa drive some alpha diversity estimates, few PCR replicates and low read sampling depths may be sufficient for many biological applications of eDNA metabarcoding. However, because rare taxa are recovered stochastically, eDNA metabarcoding may never fully recover the true amplifiable alpha diversity in an eDNA extract. Rare taxa drive PCR replicate outliers of alpha and beta diversity and lead to dispersion differences at different read sampling depths. We conclude that researchers should consider the complexity and unevenness of a community when choosing analytical approaches, read sampling depths, and filtering thresholds to arrive at stable estimates.

Altered Plasma Fatty Acids Associate with Gut Microbial Composition in Common Variable Immunodeficiency

PURPOSE

Fatty acid (FA) abnormalities are found in various inflammatory disorders and have been related to disturbed gut microbiota. Patients with common variable immunodeficiency (CVID) have inflammatory complications associated with altered gut microbial composition. We hypothesized that there is an altered FA profile in CVID patients, related to gut microbial dysbiosis.

METHODS

Plasma FAs were measured in 39 CVID patients and 30 healthy controls. Gut microbial profile, a food frequency questionnaire, and the effect of the oral antibiotic rifaximin were investigated in CVID patients.

RESULTS

The n-3 polyunsaturated fatty acids (PUFAs), eicosapentaenoic acid (EPA) (1.4 [1.0-1.8] vs. 1.9 [1.2-2.5], median (IQR), P < 0.05), and docosahexaenoic acid (DHA) (3.2 [2.4-3.9] vs. 3.5 [2.9-4.3], P < 0.05), all values expressed as weight percent of total plasma FAs, were reduced in CVID compared to controls. Also, n-6 PUFAs (34.3 ± 3.4 vs. 37.1 ± 2.8, mean ± SD, P < 0.001) and linoleic acid (LA) (24.5 ± 3.3 vs. 28.1 ± 2.7, P < 0.0001) and the FA anti-inflammatory index (98.9 [82.1-119.4] vs. 117.0 [88.7-153.1], median (IQR), P < 0.05) were reduced in CVID. The microbial alpha diversity was positively associated with plasma n-6 PUFAs (r = 0.41, P < 0.001) and LA (r = 0.51, P < 0.001), but not n-3 PUFAs (P = 0.78). Moreover, a 2-week course of rifaximin significantly reduced the proportion of n-6 PUFAs (P = 0.04, UNIANOVA). Serum immunoglobulin G (IgG) levels correlated with plasma n-3 PUFAs (rho = 0.36, P = 0.03) and DHA (rho = 0.41, P = 0.009).

CONCLUSION

We found a potentially unfavorable FA profile in CVID, related to low IgG levels. High plasma n-6 PUFAs were related to increased gut microbial diversity and altered by rifaximin therapy.

Moderate-Intensity Physical Exercise Affects the Exercise Performance and Gut Microbiota of Mice

The gut microbiota is closely associated with the health of the host and is affected by many factors, including exercise. In this study, we compared the gut microbial changes and exercise performance over a 14-week period in mice that performed exercise (NE; n = 15) and mice that did not perform exercise (NC; n = 15). Mice were subjected to stool collection and exercise tests one week prior to adaptive training and after 2, 6, 10, and 14 weeks of exercise. Bacteria associated with the stool samples were assessed via Illumina-based 16S rRNA gene sequencing. While there was no significant difference in body weight between the groups (p > 0.05), the NE group had a significantly higher exercise performance from weeks 2-14 (p < 0.01) and lower fat coefficient (p < 0.01) compared with the NC group. The Shannon index of the gut microbiota in the NC group was higher than that in the NE group at weeks 6 and 10, and the Chao1 index was higher than that in the NE group at week 14. Exercise performance positively correlated with the relative abundance of Phascolarctobacterium. Grouped time series data analysis demonstrated that Bifidobacteria, Coprococcus, and one unnamed genus in the Clostridiales order were significantly increased in the NE group, which correlated with increased glucose, flavonoid, arginine, and proline metabolism. In conclusion, moderate-intensity treadmill exercise significantly increased the exercise performance of mice and changed the core bacteria and bacterial metabolic activity. These results provide a reference for studying the effects of exercise intervention and exercise performance on the gut microbiota of mice.

A posteriori dietary patterns better explain variations of the gut microbiome than individual markers in the American Gut Project

BACKGROUND

Individual diet components and specific dietary regimens have been shown to impact the gut microbiome.

OBJECTIVES

Here, we explored the contribution of long-term diet by searching for dietary patterns that would best associate with the gut microbiome in a population-based cohort.

METHODS

Using a priori and a posteriori approaches, we constructed dietary patterns from an FFQ completed by 1800 adults in the American Gut Project. Dietary patterns were defined as groups of participants or combinations of food variables (factors) driven by criteria ranging from individual nutrients to overall diet. We associated these patterns with 16S ribosomal RNA-based gut microbiome data for a subset of 744 participants.

RESULTS

Compared to individual features (e.g., fiber and protein), or to factors representing a reduced number of dietary features, 5 a posteriori dietary patterns based on food groups were best associated with gut microbiome beta diversity (P ≤ 0.0002). Two patterns followed Prudent-like diets-Plant-Based and Flexitarian-and exhibited the highest Healthy Eating Index 2010 (HEI-2010) scores. Two other patterns presented Western-like diets with a gradient in HEI-2010 scores. A fifth pattern consisted mostly of participants following an Exclusion diet (e.g., low carbohydrate). Notably, gut microbiome alpha diversity was significantly lower in the most Western pattern compared to the Flexitarian pattern (P ≤ 0.009), and the Exclusion diet pattern was associated with low relative abundance of Bifidobacterium (P ≤ 1.2 × 10-7), which was better explained by diet than health status.

CONCLUSIONS

We demonstrated that global-diet a posteriori patterns were more associated with gut microbiome variations than individual dietary features among adults in the United States. These results confirm that evaluating diet as a whole is important when studying the gut microbiome. It will also facilitate the design of more personalized dietary strategies in general populations.

Effects of Helicobacter pylori Infection on the Oral Microbiota of Reflux Esophagitis Patients

The human oral microbiota plays a vital role in maintaining metabolic homeostasis. To explore the relationship between Helicobacter pylori (Hp) and reflux esophagitis, we collected 86 saliva samples from reflux esophagitis patients (RE group) and 106 saliva samples from healthy people (C group) for a high-throughput sequencing comparison. No difference in alpha diversity was detected between the RE and the C groups, but beta diversity of the RE group was higher than the C group. Bacteroidetes was more abundant in the RE group, whereas Firmicutes was more abundant in the C group. The linear discriminant analysis effect size analysis demonstrated that the biomarkers of the RE group were Prevotella, Veillonella, Leptotrichia, and Actinomyces, and the biomarkers of the C group were Lautropia, Gemella, Rothia, and Streptococcus. The oral microbial network structure of the C group was more complex than that of the RE group. Second, to explore the effect of Hp on the oral microbiota of RE patients, we performed the ¹⁴C-urea breath test on 45 of the 86 RE patients. We compared the oral microbiota of 33 Hp-infected reflux esophagitis patients (REHpp group) and 12 non-Hp-infected reflux esophagitis patients (REHpn group). No difference in alpha diversity was observed between the REHpn and REHpp groups, and beta diversity of the REHpp group was significantly lower than that of the REHpn group. The biomarkers in the REHpp group were Veillonella, Haemophilus, Selenomonas, Megasphaera, Oribacterium, Butyrivibrio, and Campylobacter; and the biomarker in the REHpn group was Stomatobaculum. Megasphaera was positively correlated with Veillonella in the microbial network of the REHpp group. The main finding of this study is that RE disturbs the human oral microbiota, such as increased beta diversity. Hp infection may inhibit this disorderly trend.

The Influence of Habitat on Viral Diversity in Neotropical Rodent Hosts

Rodents are important reservoirs of numerous viruses, some of which have significant impacts on public health. Ecosystem disturbances and decreased host species richness have been associated with the emergence of zoonotic diseases. In this study, we aimed at (a) characterizing the viral diversity in seven neotropical rodent species living in four types of habitats and (b) exploring how the extent of environmental disturbance influences this diversity. Through a metagenomic approach, we identified 77,767 viral sequences from spleen, kidney, and serum samples. These viral sequences were attributed to 27 viral families known to infect vertebrates, invertebrates, plants, and amoeba. Viral diversities were greater in pristine habitats compared with disturbed ones, and lowest in peri-urban areas. High viral richness was observed in savannah areas. Differences in these diversities were explained by rare viruses that were generally more frequent in pristine forest and savannah habitats. Moreover, changes in the ecology and behavior of rodent hosts, in a given habitat, such as modifications to the diet in disturbed vs. pristine forests, are major determinants of viral composition. Lastly, the phylogenetic relationships of four vertebrate-related viral families (Polyomaviridae, Flaviviridae, Togaviridae, and Phenuiviridae) highlighted the wide diversity of these viral families, and in some cases, a potential risk of transmission to humans. All these findings provide significant insights into the diversity of rodent viruses in Amazonia, and emphasize that habitats and the host’s dietary ecology may drive viral diversity. Linking viral richness and abundance to the ecology of their hosts and their responses to habitat disturbance could be the starting point for a better understanding of viral emergence and for future management of ecosystems.

Fish taxonomic, functional, and phylogenetic diversity and their vulnerabilities in the largest river in southeastern China

Freshwater biodiversity is currently under multiple threats. Conservation of freshwater fish biodiversity needs to be prioritized because natural conservation resources are always limited.Samples were collected at 24 sites in the Min River, the largest basin in southeastern China. Taxonomic, functional, and phylogenetic diversity were analyzed. Biodiversity vulnerability was measured by removing one species each time out of the community with replacement.Results suggested that hotspots for taxonomic and phylogenetic diversity were located at two impounded sites, while for functional diversity were those sites with no upstream dams. Little congruence was observed between taxonomic, functional, and phylogenetic diversity. Fragmentation of river network connectivity caused by dams was a significant factor affecting the biodiversity patterns. Beta turnover was the driving component for beta diversity, indicating that biodiversity dissimilarity along the river was mostly explained by environmental sorting. Fifteen out of 16 species that contributed the most to different facets of biodiversity were mostly endemic, either they had distinctive functional traits or they were the most prevalent species. Sites with the highest diversity vulnerability were characterized by these distinctive species. Functional diversity was more vulnerable to species loss comparing with the other two biodiversity facets.Prioritizing those biodiversity hotspots, sites with extreme functional vulnerability, and those distinctive endemic species which contributed the most to biodiversity vulnerability is suggested in the Min River. The study found evidence that congruence among different facets of biodiversity is hard to achieve, and functional diversity is the most vulnerable in a freshwater system fragmented by intensive dam constructions. This work will help to develop systematic conservation planning from the perspective of different biodiversity facets.

Plastome-based phylogeny improves community phylogenetics of subtropical forests in China

Phylogenetic trees have been extensively used in community ecology. However, how the phylogeny construction affects ecological inferences is poorly understood. In this study, we constructed three different types of phylogenetic trees (a synthetic-tree generated using V.PhyloMaker, a barcode-tree generated using rbcL+matK+trnH-psbA, and a plastome-tree generated from plastid genomes) that represented an increasing level of phylogenetic resolution among 580 woody plant species from six forest dynamic plots in subtropical evergreen broadleaved forests of China. We then evaluated the performance of each phylogeny in estimations of community phylogenetic structure, turnover and phylogenetic signal in functional traits. As expected, the plastome-tree was most resolved and most supported for relationships among species. For local phylogenetic structure, the three trees showed consistent results with Faith's PD and MPD; however, only the synthetic-tree produced significant clustering patterns using MNTD for some plots. For phylogenetic turnover, contrasting results between the molecular trees and the synthetic-tree occurred only with nearest neighbor distance. The barcode-tree agreed more with the plastome-tree than the synthetic-tree for both phylogenetic structure and turnover. For functional traits, both the barcode-tree and plastome-tree detected phylogenetic signal in maximum height, but only the plastome-tree detected signal in leaf width. This is the first study that uses plastid genomes in large-scale community phylogenetics. Our results highlight the improvement of plastome-trees over barcode-trees and synthetic-trees for the analyses studied here. Our results also point to the possibility of type I and II errors in estimation of phylogenetic structure and turnover and detection of phylogenetic signal when using synthetic-trees.

Marine biodiversity refugia in a climate-sensitive subarctic shelf

The subarctic shelf of the Eastern Bering Sea (EBS) is one of the world's most productive marine environments, exposed to drastic climate changes characterized by extreme fluctuations in temperature, sea ice concentration, timing, and duration. These climatic changes elicit profound responses in species distribution, abundance, and community composition. Here, we examined the patterns of alpha and temporal beta diversity of 159 marine taxa (66 vertebrates and 93 invertebrate species) from 29 years (1990-2018) of species observations from the NOAA bottom trawl surveys in the EBS. Based on these data, we identified geographically distinct refugial zones in the northern and southern regions of the middle shelf, defined by high species richness and similarity in community species composition over time. These refugial zones harbor higher frequencies of occurrence for representative taxa relative to the regions outside of refugia. We also explored the primary environmental factors structuring marine biodiversity distributions, which underpinned the importance of the winter sea ice concentration to alpha and temporal beta diversity. The spatial biodiversity distributions between high and low winter sea ice regimes highlighted contrasting signals. In particular, the latter showed elevated species richness compared to the former. Further, the temporal beta diversity between the high and low winter sea ice periods underpinned an overall increase in the compositional similarity of marine communities in the EBS. Despite these spatiotemporal differences in biodiversity distributions, the identified refugia represent safe havens of marine biodiversity in the EBS. Distinguishing these areas can help facilitate conservation and management efforts under accelerated and ongoing climatic changes.

Impact of Dietary Protein Content on Soil Bacterial and Fungal Communities in a Rice-Crab Co-culture System

Although co-culture of paddy fields with aquatic animals is lucrative, the ecological impacts of high-protein content entering the agricultural soil via animal pellet feed and feces have not been well studied. Moreover, the effects of dietary protein on soils and soil microbes remain unclear. To elucidate this, we examined soil bacterial and fungal community composition and temporal changes in paddy fields subjected to different protein-content diets via 16S/18S rRNA gene amplicon sequencing analysis with a high-throughput next-generation sequencer. MiSeq sequencing revealed that protein content significantly impacted fungal community structure. High-protein diets reduced bacterial community diversity and relative abundance in both July and October. The phylum-level bacterial taxonomic composition was not affected by diet treatment, while in fungi, a major phylum-level shift was evident. Hierarchically clustered analysis showed that high-protein diets significantly reduced the relative abundance of Brevundimonas in both July and October. Saprotrophic macrofungal diversity was negatively related to dietary protein content. Considering microbial community structure and environmental factors, ca. 15% protein content is appropriate for the rice-crab co-culture system that we studied.

Fecal Methylmercury Correlates With Gut Microbiota Taxa in Pacific Walruses (Odobenus rosmarus divergens)

OBJECTIVES

Methylmercury metabolism was investigated in Pacific walruses (Odobenus rosmarus divergens) from St. Lawrence Island, Alaska, United States.

METHODS

Total mercury and methylmercury concentrations were measured in fecal samples and paired colon samples (n = 16 walruses). Gut microbiota composition and diversity were determined using 16S rRNA gene sequencing. Associations between fecal and colon mercury and the 24 most prevalent gut microbiota taxa were investigated using linear models.

RESULTS

In fecal samples, the median values for total mercury, methylmercury, and %methylmercury (of total mercury) were 200 ng/g, 4.7 ng/g, and 2.5%, respectively, while in colon samples, the median values for the same parameters were 28 ng/g, 7.8 ng/g, and 26%, respectively. In fecal samples, methylmercury was negatively correlated with one Bacteroides genus, while members of the Oscillospirales order were positively correlated with both methylmercury and %methylmercury (of total mercury). In colon samples, %methylmercury (of total mercury) was negatively correlated with members of two genera, Romboutsia and Paeniclostridium.

CONCLUSIONS

Median %methylmercury (of total mercury) was 10 times higher in the colon compared to the fecal samples, suggesting that methylmercury was able to pass through the colon into systemic circulation. Fecal total mercury and/or methylmercury concentrations in walruses were comparable to some human studies despite differences in seafood consumption rates, suggesting that walruses excreted less mercury. There are no members (at this time) of the Oscillospirales order which are known to contain the genes to methylate mercury, suggesting the source of methylmercury in the gut was from diet and not in vivo methylation.

Exploring Multiple Aspects of Taxonomic and Functional Diversity in Microphytobenthic Communities: Effects of Environmental Gradients and Temporal Changes

Biodiversity has traditionally been quantified using taxonomic information but the importance of also considering its functional characteristics has recently gained an increasing attention among microorganisms. However, studies exploring multiple aspects of taxonomic and functional diversity and their temporal variations are scarce for diatoms, which is one of the most important microbial groups in aquatic ecosystems. Here, our aim was to examine the taxonomic and functional alpha and beta diversities of diatoms in a coastal rock pool system characterized by a naturally high environmental heterogeneity. We also investigated the temporal differences in the diversity patterns and drivers. The relationship between the species richness and functional dispersion was temporally coherent, such that species-poor communities tended to be functionally clustered. The trend between the species richness and taxonomic uniqueness of community composition was temporally inconsistent, changing from negative to non-significant over time. Conductivity or distance to the sea or both were key determinants of species richness, functional dispersion, and uniqueness of community composition. The increase of community dissimilarity with an increasing environmental distance was stronger for the taxonomic than the functional composition. Our results suggest that even minor decreases in the species richness may result in a lowered functional diversity and decreased ecosystem functioning. Species-poor ecosystems may, however, have unique species compositions and high contributions to regional biodiversity. Despite changing the species compositions along the environmental gradients, communities may remain to have a high functional similarity and robustness in the face of environmental changes. Our results highlight the advantage of considering multiple biodiversity metrics and incorporating a temporal component for a deeper understanding of the effects of environmental changes on microbial biodiversity.

Organic Matter Composition at Ocean Station Papa Affects Its Bioavailability, Bacterioplankton Growth Efficiency and the Responding Taxa

The bioavailability of organic matter (OM) to marine heterotrophic bacterioplankton is determined by both the chemical composition of OM and the microbial community composition. In the current study, changes in OM bioavailability were identified at Ocean Station Papa as part of the 2018 Export Processes in the Ocean from Remote Sensing (EXPORTS) field study. Removal rates of carbon (C) in controlled experiments were significantly correlated with the initial composition of total hydrolyzable amino acids, and C removal rates were high when the amino acid degradation index suggested a more labile composition. Carbon remineralization rates averaged 0.19 ± 0.08 μmol C L-1 d-1 over 6-10 days while bacterial growth efficiencies averaged 31 ± 7%. Amino acid composition and tandem mass spectrometry analysis of compound classes also revealed transformations to a more degraded OM composition during experiments. There was a log2-fold increase in the relative abundances of 16S rDNA-resolved bacterioplankton taxa in most experiments by members of the Methylophilaceae family (OM43 genus) and KI89A order. Additionally, when OM was more bioavailable, relative abundances increased by at least threefold for the classes Bacteroidetes (Flavobacteriaceae NS2b genus), Alphaproteobacteria (Rhodobacteraceae Sulfitobacter genus), and Gammaproteobacteria (Alteromonadales and Ectothiorhodospiraceae orders). Our data suggest that a diverse group of bacterioplankton was responsible for removing organic carbon and altering the OM composition to a more degraded state. Elevated community diversity, as inferred from the Shannon-Wiener H index, may have contributed to relatively high growth efficiencies by the bacterioplankton. The data presented here shed light on the interconnections between OM bioavailability and key bacterioplankton taxa for the degradation of marine OM.

Gut Microbiota May Play a Significant Role in the Pathogenesis of Graves' Disease

Background: Gut microbiota are considered to be intrinsic regulators of thyroid autoimmunity. We designed a cross-sectional study to examine the makeup and metabolic function of microbiota in Graves' disease (GD) patients, with the ultimate aim of offering new perspectives on the diagnosis and treatment of GD. Methods: The 16S ribosomal RNA (rRNA) V3-V4 DNA regions of microbiota were obtained from fecal samples collected from 45 GD patients and 59 controls. Microbial differences between the two groups were subsequently analyzed based on high-throughput sequencing. Results: Compared with controls, GD patients had reduced alpha diversity (p < 0.05). At the phylum level, GD patients had a significantly lower proportion of Firmicutes (p = 0.008) and a significantly higher proportion of Bacteroidetes (p = 0.002) compared with the controls. At the genus level, GD patients had greater numbers of Bacteroides and Lactobacillus, although fewer Blautia, [Eubacterium]_hallii_group, Anaerostipes, Collinsella, Dorea, unclassified_f_Peptostreptococcaceae, and [Ruminococcus]_torques_group than controls (all p < 0.05). Subgroup analysis of GD patients revealed that Lactobacillus may play a key role in the pathogenesis of autoimmune thyroid diseases. Nine distinct genera showed significant correlations with certain thyroid function tests. Functional prediction revealed that Blautia may be an important microbe in certain metabolic pathways that occur in the hyperthyroid state. In addition, linear discriminant analysis (LDA) and effect size (LEfSe) analysis showed that there were significant differences in the levels of 18 genera between GD patients and controls (LDA >3.0, all p < 0.05). A diagnostic model using the top nine genera had an area under the curve of 0.8109 [confidence interval: 0.7274-0.8945]. Conclusions: Intestinal microbiota are different in GD patients. The microbiota we identified offer an alternative noninvasive diagnostic methodology for GD. Microbiota may also play a role in thyroid autoimmunity, and future research is needed to further elucidate the role.

Effects of site-selection bias on estimates of biodiversity change

Estimates of biodiversity change are essential for the management and conservation of ecosystems. Accurate estimates rely on selecting representative sites, but monitoring often focuses on sites of special interest. How such site-selection biases influence estimates of biodiversity change is largely unknown. Site-selection bias potentially occurs across four major sources of biodiversity data, decreasing in likelihood from citizen science, museums, national park monitoring, and academic research. We defined site-selection bias as a preference for sites that are either densely populated (i.e., abundance bias) or species rich (i.e., richness bias). We simulated biodiversity change in a virtual landscape and tracked the observed biodiversity at a sampled site. The site was selected either randomly or with a site-selection bias. We used a simple spatially resolved, individual-based model to predict the movement or dispersal of individuals in and out of the chosen sampling site. Site-selection bias exaggerated estimates of biodiversity loss in sites selected with a bias by on average 300-400% compared with randomly selected sites. Based on our simulations, site-selection bias resulted in positive trends being estimated as negative trends: richness increase was estimated as 0.1 in randomly selected sites, whereas sites selected with a bias showed a richness change of -0.1 to -0.2 on average. Thus, site-selection bias may falsely indicate decreases in biodiversity. We varied sampling design and characteristics of the species and found that site-selection biases were strongest in short time series, for small grains, organisms with low dispersal ability, large regional species pools, and strong spatial aggregation. Based on these findings, to minimize site-selection bias, we recommend use of systematic site-selection schemes; maximizing sampling area; calculating biodiversity measures cumulatively across plots; and use of biodiversity measures that are less sensitive to rare species, such as the effective number of species. Awareness of the potential impact of site-selection bias is needed for biodiversity monitoring, the design of new studies on biodiversity change, and the interpretation of existing data.

Early ecological succession patterns of bacterial, fungal and plant communities along a chronosequence in a recently deglaciated area of the Italian Alps

In this study, the early ecological succession patterns of Forni Glacier (Ortles-Cevedale group, Italian Alps) forefield along an 18-year long chronosequence (with a temporal resolution of 1 year) has been reported. Bacterial and fungal community structures were inferred by high-throughput sequencing of 16S rRNA gene and ITS, respectively. In addition, the occurrence of both herbaceous and arboreous plants was also recorded at each plot. A significant decrease of alpha-diversity in more recently deglaciated areas was observed for both bacteria and plants. Time since deglaciation and pH affected the structure of both fungal and bacterial communities. Pioneer plants could be a major source of colonization for both bacterial and fungal communities. Consistently, some of the most abundant bacterial taxa and some of those significantly varying with pH along the chronosequence (Polaromonas, Granulicella, Thiobacillus, Acidiferrobacter) are known to be actively involved in rock-weathering processes due to their chemolithotrophic metabolism, thus suggesting that the early phase of the chronosequence could be mainly shaped by the biologically controlled bioavailability of metals and inorganic compounds. Fungal communities were dominated by ascomycetous filamentous fungi and basidiomycetous yeasts. Their role as cold-adapted organic matter decomposers, due to their heterotrophic metabolism, was suggested.