A new method for non-parametric multivariate analysis of variance

Biodiversity and ecology of flower-associated actinomycetes in different flowering stages of Protea repens

Actinomycete bacteria have previously been reported from reproductive structures (infructescences) of Protea (sugarbush/suikerbos) species, a niche dominated by fungi in the genera Knoxdaviesia and Sporothrix. It is probable that these taxa have symbiotic interactions, but a lack of knowledge regarding their diversity and general ecology precludes their study. We determined the diversity of actinomycetes within Protea repens inflorescence buds, open inflorescences, young and mature infructescences, and leaf litter surrounding these trees. Since the P. repens habitat is fire-prone, we also considered the potential of these bacteria to recolonise infructescences after fire. Actinomycetes were largely absent from flower buds and inflorescences but were consistently present in young and mature infructescences. Two Streptomyces spp. were the most consistent taxa recovered, one of which was also routinely isolated from leaf litter. Lower colonisation rates were evident in samples from a recently burnt site. One of the most consistent taxa isolated from older trees in the unburnt site was absent from this site. Our findings show that P. repens has a distinct community of actinomycetes dominated by a few species. These communities change over time and infructescence developmental stage, season and the age of the host population. Mature infructescences appear to be important sources of inoculum for some of the actinomycetes, seemingly disrupted by fire. Increased fire frequency limiting maturation of P. repens infructescences could thus impact future actinomycete colonisation in the landscape. Streptomyces spp. are likely to share this niche with the ophiostomatoid fungi, which merits further study regarding their interactions and mode of transfer.

Maize Phyllosphere Microbial Community Niche Development Across Stages of Host Leaf Growth

Background: The phyllosphere hosts a variety of microorganisms, including bacteria, which can play a positive role in the success of the host plant. Bacterial communities in the phylloplane are influenced by both biotic and abiotic factors, including host plant surface topography and chemistry, which change in concert with microbial communities as the plant leaves develop and age.

Methods: We examined how the Zea mays L. leaf microbial community structure changed with plant age. Ribosomal spacer length and scanning electron microscopic imaging strategies were used to assess microbial community composition across maize plant ages, using a novel staggered experimental design.

Results: Significant changes in community composition were observed for both molecular and imaging analyses, and the two analysis methods provided complementary information about bacterial community structure within each leaf developmental stage.

Conclusions: Both taxonomic and cell-size trait patterns provided evidence for niche-based contributions to microbial community development on leaves.

Seasonality Affects the Diversity and Composition of Bacterioplankton Communities in Dongjiang River, a Drinking Water Source of Hong Kong

Water quality ranks the most vital criterion for rivers serving as drinking water sources, which periodically changes over seasons. Such fluctuation is believed associated with the state shifts of bacterial community within. To date, seasonality effects on bacterioplankton community patterns in large rivers serving as drinking water sources however, are still poorly understood. Here we investigated the intra-annual bacterial community structure in the Dongjiang River, a drinking water source of Hong Kong, using high-throughput pyrosequencing in concert with geochemical property measurements during dry, and wet seasons. Our results showed that Proteobacteria, Actinobacteria, and Bacteroidetes were the dominant phyla of bacterioplankton communities, which varied in composition, and distribution from dry to wet seasons, and exhibited profound seasonal changes. Actinobacteria, Bacteroidetes, and Cyanobacteria seemed to be more associated with seasonality that the relative abundances of Actinobacteria, and Bacteroidetes were significantly higher in the dry season than those in the wet season (p < 0.01), while the relative abundance of Cyanobacteria was about 10-fold higher in the wet season than in the dry season. Temperature and NO3--N concentration represented key contributing factors to the observed seasonal variations. These findings help understand the roles of various bacterioplankton and their interactions with the biogeochemical processes in the river ecosystem.

Microbial communities in marine sediments modify success of an invasive macrophyte

Invasive plants have extensive impacts on ecosystem function and biodiversity globally. Our inability to manage invasive species stems in part from a lack of understanding of the processes that control their successful establishment and spread. To date, studies have largely considered how above-ground processes control native/invasive plant interactions. Emerging research from terrestrial and wetland ecosystems demonstrates that below-ground processes under microbial control can determine the outcome of interactions between native and invasive plants. Whether sediment microbes modify the success of invasive macrophytes in marine ecosystems is untested, despite marine sediment microbes controlling many ecological processes (e.g. nutrient cycling) comparable to those in terrestrial ecosystems. We first show that sediment bacterial communities differ between the native seagrass Zostera capricorni and the invasive alga Caulerpa taxifolia and that those differences relate to functional changes in sulfur cycling between the macrophytes. Second, by experimentally manipulating the microbial communities we show that intact microbial communities in Z. capricorni sediments provide biotic resistance by reducing C. taxifolia fragment growth 119% compared to when they are inactive, and intact microbial communities in C. taxifolia sediments have positive feedbacks by increasing fragment growth 200%. Thus, similar to terrestrial ecosystems, microorganisms appear to indirectly control the success of invasive macrophytes in marine ecosystems.

Skin bacterial microbiome of a generalist Puerto Rican frog varies along elevation and land use gradients

Host-associated microbial communities are ubiquitous among animals, and serve important functions. For example, the bacterial skin microbiome of amphibians can play a role in preventing or reducing infection by the amphibian chytrid fungus, Batrachochytrium dendrobatidis. Evidence suggests that environmental bacteria likely serve as a source pool for at least some of the members of the amphibian skin bacterial community, underscoring the potential for local environmental changes to disrupt microbial community source pools that could be critical to the health of host organisms. However, few studies have assessed variation in the amphibian skin microbiome along clear environmental gradients, and so we know relatively little about how local environmental conditions influence microbiome diversity. We sampled the skin bacterial communities of Coqui frogs, Eleutherodactylus coqui (N = 77), along an elevational gradient in eastern Puerto Rico (0-875 m), with transects in two land use types: intact forest (N = 4 sites) and disturbed (N = 3 sites) forest. We found that alpha diversity (as assessed by Shannon, Simpson, and Phylogenetic Diversity indices) varied across sites, but this variation was not correlated with elevation or land use. Beta diversity (community structure), on the other hand, varied with site, elevation and land use, primarily due to changes in the relative abundance of certain bacterial OTUs (∼species) within these communities. Importantly, although microbiome diversity varied, E. coqui maintained a common core microbiota across all sites. Thus, our findings suggest that environmental conditions can influence the composition of the skin microbiome of terrestrial amphibians, but that some aspects of the microbiome remain consistent despite environmental variation.

Culture Media and Individual Hosts Affect the Recovery of Culturable Bacterial Diversity from Amphibian Skin

One current challenge in microbial ecology is elucidating the functional roles of the large diversity of free-living and host-associated bacteria identified by culture-independent molecular methods. Importantly, the characterization of this immense bacterial diversity will likely require merging data from culture-independent approaches with work on bacterial isolates in culture. Amphibian skin bacterial communities have become a recent focus of work in host-associated microbial systems due to the potential role of these skin bacteria in host defense against the pathogenic fungus Batrachochytrium dendrobatidis (Bd), which is associated with global amphibian population declines and extinctions. As there is evidence that some skin bacteria may inhibit growth of Bd and prevent infection in some cases, there is interest in using these bacteria as probiotic therapy for conservation of at-risk amphibians. In this study, we used skin swabs from American toads (Anaxyrus americanus) to: (1) assess the diversity and community structure of culturable amphibian skin bacteria grown on high and low nutrient culture media, (2) determine which culture media recover the highest proportion of the total skin bacterial community of individual toads relative to culture-independent data, and (3) assess whether the plated communities from the distinct media types vary in their ability to inhibit Bd growth in in-vitro assays. Overall, we found that culture media with low nutrient concentrations facilitated the growth of more diverse bacterial taxa and grew distinct communities relative to media with higher nutrient concentrations. Use of low nutrient media also resulted in culturing proportionally more of the bacterial diversity on individual toads relative to the overall community defined using culture-independent methods. However, while there were differences in diversity among media types, the variation among individual hosts was greater than variation among media types, suggesting that swabbing more individuals in a population is the best way to maximize culture collections, regardless of media type. Lastly, the function of the plated communities against Bd did not vary across culture media type or between high and low nutrient media. These results inform current efforts for developing a probiotic-based approach for amphibian conservation and help to ensure that culture collections are capturing the majority of the important diversity in these systems.

Mutation supply and the repeatability of selection for antibiotic resistance

Whether evolution can be predicted is a key question in evolutionary biology. Here we set out to better understand the repeatability of evolution, which is a necessary condition for predictability. We explored experimentally the effect of mutation supply and the strength of selective pressure on the repeatability of selection from standing genetic variation. Different sizes of mutant libraries of antibiotic resistance gene TEM-1 β-lactamase in Escherichia coli, generated by error-prone PCR, were subjected to different antibiotic concentrations. We determined whether populations went extinct or survived, and sequenced the TEM gene of the surviving populations. The distribution of mutations per allele in our mutant libraries followed a Poisson distribution. Extinction patterns could be explained by a simple stochastic model that assumed the sampling of beneficial mutations was key for survival. In most surviving populations, alleles containing at least one known large-effect beneficial mutation were present. These genotype data also support a model which only invokes sampling effects to describe the occurrence of alleles containing large-effect driver mutations. Hence, evolution is largely predictable given cursory knowledge of mutational fitness effects, the mutation rate and population size. There were no clear trends in the repeatability of selected mutants when we considered all mutations present. However, when only known large-effect mutations were considered, the outcome of selection is less repeatable for large libraries, in contrast to expectations. We show experimentally that alleles carrying multiple mutations selected from large libraries confer higher resistance levels relative to alleles with only a known large-effect mutation, suggesting that the scarcity of high-resistance alleles carrying multiple mutations may contribute to the decrease in repeatability at large library sizes.

Correlations of age and growth rate with microbiota composition in Atlantic cod (Gadus morhua) larvae

Little information is available on the link between host development (growth rate and ontogeny) and the composition of the microbiota in fish larvae. This study was carried out to examine potential correlations of microbiota composition with age and growth rate of Atlantic cod larvae. Small and large cod larvae of the same age, representing slow and fast growing individuals, were sampled 10 times during a period of 42 days post hatching (dph), and the composition of the larval microbiota was investigated using a PCR/DGGE (Denaturing Gradient Gel Electrophoresis) strategy. We found significant differences in the intestinal microbiota of small and large larvae of the same age for 4 of the 10 age stages studied. We further found that the variation in the composition of the larval microbiota was more strongly correlated to age than to growth rate for larvae up to 28 dph, whereas for the older larvae growth rate and age was equally correlated to the composition of the microbiota. These results indicate that larval development may structure the microbiota through a change in selection pressure due to host-microbe and microbe-microbe interactions, and that the composition of the microbiota may influence larval development through improved energy gain.

Microbiome succession during ammonification in eelgrass bed sediments

BACKGROUND

Eelgrass (Zostera marina) is a marine angiosperm and foundation species that plays an important ecological role in primary production, food web support, and elemental cycling in coastal ecosystems. As with other plants, the microbial communities living in, on, and near eelgrass are thought to be intimately connected to the ecology and biology of eelgrass. Here we characterized the microbial communities in eelgrass sediments throughout an experiment to quantify the rate of ammonification, the first step in early remineralization of organic matter, also known as diagenesis, from plots at a field site in Bodega Bay, CA.

METHODS

Sediment was collected from 72 plots from a 15 month long field experiment in which eelgrass genotypic richness and relatedness were manipulated. In the laboratory, we placed sediment samples (n = 4 per plot) under a N₂ atmosphere, incubated them at in situ temperatures (15 °C) and sampled them initially and after 4, 7, 13, and 19 days to determine the ammonification rate. Comparative microbiome analysis using high throughput sequencing of 16S rRNA genes was performed on sediment samples taken initially and at seven, 13 and 19 days to characterize changes in the relative abundances of microbial taxa throughout ammonification.

RESULTS

Within-sample diversity of the sediment microbial communities across all plots decreased after the initial timepoint using both richness based (observed number of OTUs, Chao1) and richness and evenness based diversity metrics (Shannon, Inverse Simpson). Additionally, microbial community composition changed across the different timepoints. Many of the observed changes in relative abundance of taxonomic groups between timepoints appeared driven by sulfur cycling with observed decreases in predicted sulfur reducers (Desulfobacterales) and corresponding increases in predicted sulfide oxidizers (Thiotrichales). None of these changes in composition or richness were associated with variation in ammonification rates.

DISCUSSION

Our results showed that the microbiome of sediment from different plots followed similar successional patterns, which we infer to be due to changes related to sulfur metabolism. These large changes likely overwhelmed any potential changes in sediment microbiome related to ammonification rate. We found no relationship between eelgrass presence or genetic composition and the microbiome. This was likely due to our sampling of bulk sediments to measure ammonification rates rather than sampling microbes in sediment directly in contact with the plants and suggests that eelgrass influence on the sediment microbiome may be limited in spatial extent. More in-depth functional studies associated with eelgrass microbiome will be required in order to fully understand the implications of these microbial communities in broader host-plant and ecosystem functions (e.g., elemental cycling and eelgrass-microbe interactions).

Changes in zooplankton community, and seston and zooplankton fatty acid profiles at the freshwater/saltwater interface of the Chowan River, North Carolina

The variability in zooplankton fatty acid composition may be an indicator of larval fish habitat quality as fatty acids are linked to fish larval growth and survival. We sampled an anadromous fish nursery, the Chowan River, during spring of 2013 in order to determine how the seston fatty acid composition varied in comparison with the zooplankton community composition and fatty acid composition during the period of anadromous larval fish residency. The seston fatty acid profiles showed no distinct pattern in relation to sampling time or location. The mesozooplankton community composition varied spatially and the fatty acid profiles were typical of freshwater species in April. The Chowan River experienced a saltwater intrusion event during May, which resulted in brackish water species dominating the zooplankton community and the fatty acid profile showed an increase in polyunsaturated fatty acids (PUFA), in particular eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA). The saltwater intrusion event was followed by an influx of freshwater due to high precipitation levels in June. The zooplankton community composition once again became dominated by freshwater species and the fatty acid profiles shifted to reflect this change; however, EPA levels remained high, particularly in the lower river. We found correlations between the seston, microzooplankton and mesozooplankton fatty acid compositions. Salinity was the main factor correlated to the observed pattern in species composition, and fatty acid changes in the mesozooplankton. These data suggest that anadromous fish nursery habitat likely experiences considerable spatial variability in fatty acid profiles of zooplankton prey and that are correlated to seston community composition and hydrodynamic changes. Our results also suggest that sufficient prey density as well as a diverse fatty acid composition is present in the Chowan River to support larval fish production.

Breast tissue, oral and urinary microbiomes in breast cancer

It has long been proposed that the gut microbiome contributes to breast carcinogenesis by modifying systemic estrogen levels. This is often cited as a possible mechanism linking breast cancer and high-fat, low-fiber diets as well as antibiotic exposure, associations previously identified in population-based studies. More recently, a distinct microbiome has been identified within breast milk and tissue, but few studies have characterized differences in the breast tissue microbiota of patients with and without cancer, and none have investigated distant body-site microbiomes outside of the gut. We hypothesize that cancerous breast tissue is associated with a microbiomic profile distinct from that of benign breast tissue, and that microbiomes of more distant sites, the oral cavity and urinary tract, will reflect dysbiosis as well. Fifty-seven women with invasive breast cancer undergoing mastectomy and 21 healthy women undergoing cosmetic breast surgery were enrolled. The bacterial 16S rRNA gene was amplified from urine, oral rinse and surgically collected breast tissue, sequenced, and processed through a QIIME-based bioinformatics pipeline. Cancer patient breast tissue microbiomes clustered significantly differently from non-cancer patients (p=0.03), largely driven by decreased relative abundance of Methylobacterium in cancer patients (median 0.10 vs. 0.24, p=0.03). There were no significant differences in oral rinse samples. Differences in urinary microbiomes were largely explained by menopausal status, with peri/postmenopausal women showing decreased levels of Lactobacillus. Independent of menopausal status, however, cancer patients had increased levels of gram-positive organisms including Corynebacterium (p<0.01), Staphylococcus (p=0.02), Actinomyces (p<0.01), and Propionibacteriaceae (p<0.01). Our observations suggest that the local breast microbiota differ in patients with and without breast cancer. Cancer patient urinary microbiomes were characterized by increased levels of gram-positive organisms in this study, but need to be further studied in larger cohorts.

Use of multivariate control charts to assess the status of reef fish assemblages in the Northwestern Hawaiian Islands

A distance-based multivariate control chart is a useful tool for ecological monitoring to detect changes in biological community resulting from natural or anthropogenic disturbances at permanent monitoring sites. It is based on a matrix of any distances or dissimilarities among observations obtained from species composition and abundance data, and bootstrapping techniques are used to set upper confidence bounds that trigger an alarm for further investigations. We extended the use of multivariate control charts to stratified random sampling and analyzed reef fish monitoring data collected annually on shallow (≤30 m) reefs across the Northwestern Hawaiian Islands (NWHI), part of the Papahānaumokuākea Marine National Monument. Fish assemblages in the NWHI were mostly stable, with exceptions in the south region (Nihoa, Mokumanamana and French Frigate Shoals) in 2012 and 2015 where changes in the assemblage structure exceeded the upper confidence bounds of multivariate control charts. However, these were due to changes in relative abundances of native species, and potentially related to the small numbers of survey sites and relatively low coral covers at the sites, particularly in 2015. The present study showed that multivariate control charts can be used to evaluate the status of biological communities in a very large protected area. Future monitoring of fish assemblages in the Papahānaumokuākea Marine National Monument should be accompanied by specific habitat or environmental variables that are related to potential threats to its shallow-water ecosystems. This should allow for more detailed investigations into potential causes and mechanisms of changes in fish assemblages when a multivariate control chart triggers an alarm.

Mortality, recovery, and community shifts of scleractinian corals in Puerto Rico one decade after the 2005 regional bleaching event

This work analyzes the mortality, recovery, and shifts in the composition of scleractinian corals from Puerto Rico one decade after the 2005 regional coral bleaching event. Temporal and spatial patterns of coral community structure were examined using a stratified, non-random sampling approach based on five permanent transects per reef at 16 reef stations. A negative correlation between percent coral cover loss and light attenuation coefficient (Kd₄₉₀) was observed, suggesting that light attenuation, as influenced by water turbidity and depth, played a major role in coral protection during the bleaching event (“sunblock effect”). Responses of coral assemblages varied after the bleaching event, including shifts of cover from massive corals (Orbicella spp.) to opportunistic (Porites astreoides) and branching corals (Madracis auretenra, P. porites) and/or turf algae; partial recovery of reef substrate cover by O. annularis complex; and no measurable changes in coral assemblages before and after the event.

The morphology of the inner ear of squamate reptiles and its bearing on the origin of snakes

The inner ear morphology of 80 snake and lizard species, representative of a range of ecologies, is here analysed and compared to that of the fossil stem snake Dinilysia patagonica, using three-dimensional geometric morphometrics. Inner ear morphology is linked to phylogeny (we find here a strong phylogenetic signal in the data that can complicate ecological correlations), but also correlated with ecology, with Dinilysia resembling certain semi-fossorial forms (Xenopeltis and Cylindrophis), consistent with previous reports. We here also find striking resemblances between Dinilysia and some semi-aquatic snakes, such as Myron (Caenophidia, Homalopsidae). Therefore, the inner ear morphology of Dinilysia is consistent with semi-aquatic as well as semi-fossorial habits: the most similar forms are either semi-fossorial burrowers with a strong affinity to water (Xenopeltis and Cylindrophis) or amphibious, intertidal forms which shelter in burrows (Myron). Notably, Dinilysia does not cluster as closely with snakes with exclusively terrestrial or obligate burrowing habits (e.g. scolecophidians and uropeltids). Moreover, despite the above similarities, Dinilysia also occupies a totally unique morphospace, raising issues with linking it with any particular ecological category.

Effect of oil palm on the Plecoptera and Trichoptera (Insecta) assemblages in streams of eastern Amazon

The production of oil palm is expected to increase in the Amazon region. However, expansion of oil palm plantation leads to significant changes in the physical structure of aquatic ecosystems, mainly through the reduction of riparian vegetation that is essential for aquatic biodiversity. Here, we evaluated the effects of oil palm on the physical habitat structure of Amazonian stream environments and assemblages of Plecoptera and Trichoptera (PT), ​both found in these streams. We compared streams sampled in oil palm plantations (n = 13) with natural forest areas ("reference" streams, n = 8), located in the eastern Amazon, Brazil. Our results showed that oil palm streams were more likely to be in close proximity to roads, had higher pH values, and higher amounts of fine substrate deposited in the channel than reference streams. Further, these environmental changes had important effects on the aquatic invertebrate assemblages, reducing the abundance and richness of PT. Nevertheless, the genera composition of the assemblages did not differ between reference and oil palm (PERMANOVA, pseudo-F _(1,19) = 1.891; p = 0.111). We conclude that oil palm production has clear negative impacts on aquatic environments and PT assemblages in Amazonian streams. We recommend that oil palm producers invest more in planning of road networks to avoid the construction of roads near to the riparian vegetation. This planning can minimize impacts of oil palm production on aquatic systems in the Amazon.

High-protein diet improves sensitivity to cholecystokinin and shifts the cecal microbiome without altering brain inflammation in diet-induced obesity in rats

High-protein diet (HPD) curtails obesity and/or fat mass, but it is unknown whether it reverses neuroinflammation or alters glucose levels, CCK sensitivity, and gut microbiome in rats fed a Western diet (WD)-induced obesity (DIO). Male rats fed a WD (high fat and sugar) for 12 wk were switched to a HPD for 6 wk. Body composition, food intake, meal pattern, sensitivity to intraperitoneal CCK-8S, blood glucose, brain signaling, and cecal microbiota were assessed. When compared with a normal diet, WD increased body weight (9.3%) and fat mass (73.4%). CCK-8S (1.8 or 5.2 nmol/kg) did not alter food intake and meal pattern in DIO rats. Switching to a HPD for 6 wk reduced fat mass (15.7%) with a nonsignificantly reduced body weight gain, normalized blood glucose, and decreased feeding after CCK-8S. DIO rats on the WD or switched to a HPD showed comparable microbial diversity. However, in HPD versus WD rats, there was enrichment of 114 operational taxonomic units (OTUs) and depletion of 188 OTUs. Of those, Akkermansia muciniphila (enriched on a HPD), an unclassified Clostridiales, a member of the RF39 order, and a Phascolarctobacterium were significantly associated with fat mass. The WD increased cytokine expression in the hypothalamus and dorsal medulla that was unchanged by switching to HPD. These data indicate that HPD reduces body fat and restores glucose homeostasis and CCK sensitivity, while not modifying brain inflammation. In addition, expansion of cecal Akkermansia muciniphila correlated to fat mass loss may represent a potential peripheral mechanism of HPD beneficial effects.

Fertilization Shapes Bacterial Community Structure by Alteration of Soil pH

Application of chemical fertilizer or manure can affect soil microorganisms directly by supplying nutrients and indirectly by altering soil pH. However, it remains uncertain which effect mostly shapes microbial community structure. We determined soil bacterial diversity and community structure by 454 pyrosequencing the V1-V3 regions of 16S rRNA genes after 7-years (2007-2014) of applying chemical nitrogen, phosphorus and potassium (NPK) fertilizers, composted manure or their combination to acidic (pH 5.8), near-neutral (pH 6.8) or alkaline (pH 8.4) Eutric Regosol soil in a maize-vegetable rotation in southwest China. In alkaline soil, nutrient sources did not affect bacterial Operational Taxonomic Unit (OTU) richness or Shannon diversity index, despite higher available N, P, K, and soil organic carbon in fertilized than in unfertilized soil. In contrast, bacterial OTU richness and Shannon diversity index were significantly lower in acidic and near-neutral soils under NPK than under manure or their combination, which corresponded with changes in soil pH. Permutational multivariate analysis of variance showed that bacterial community structure was significantly affected across these three soils, but the PCoA ordination patterns indicated the effect was less distinct among nutrient sources in alkaline than in acidic and near-neural soils. Distance-based redundancy analysis showed that bacterial community structures were significantly altered by soil pH in acidic and near-neutral soils, but not by any soil chemical properties in alkaline soil. The relative abundance (%) of most bacterial phyla was higher in near-neutral than in acidic or alkaline soils. The most dominant phyla were Proteobacteria (24.6%), Actinobacteria (19.7%), Chloroflexi (15.3%) and Acidobacteria (12.6%); the medium dominant phyla were Bacterioidetes (5.3%), Planctomycetes (4.8%), Gemmatimonadetes (4.5%), Firmicutes (3.4%), Cyanobacteria (2.1%), Nitrospirae (1.8%), and candidate division TM7 (1.0%); the least abundant phyla were Verrucomicrobia (0.7%), Armatimonadetes (0.6%), candidate division WS3 (0.4%) and Fibrobacteres (0.3%). In addition, Cyanobacteria and candidate division TM7 were more abundant in acidic soil, whereas Gemmatimonadetes, Nitrospirae and candidate division WS3 were more abundant in alkaline soil. We conclude that after 7-years of fertilization, soil bacterial diversity and community structure were shaped more by changes in soil pH rather than the direct effect of nutrient addition.

Functional diversity and redundancy across fish gut, sediment and water bacterial communities

This article explores the functional diversity and redundancy in a bacterial metacommunity constituted of three habitats (sediment, water column and fish gut) in a coastal lagoon under anthropogenic pressure. Comprehensive functional gene arrays covering a wide range of ecological processes and stress resistance genes to estimate the functional potential of bacterial communities were used. Then, diversity partitioning was used to characterize functional diversity and redundancy within (α), between (β) and across (γ) habitats. It was showed that all local communities exhibit a highly diversified potential for the realization of key ecological processes and resistance to various environmental conditions, supporting the growing evidence that macro-organisms microbiomes harbour a high functional potential and are integral components of functional gene dynamics in aquatic bacterial metacommunities. Several levels of functional redundancy at different scales of the bacterial metacommunity were observed (within local communities, within habitats and at the metacommunity level). The results suggested a high potential for the realization of spatial ecological insurance within this ecosystem, that is, the functional compensation among microorganisms for the realization and maintenance of key ecological processes, within and across habitats. Finally, the role of macro-organisms as dispersal vectors of microbes and their potential influence on marine metacommunity dynamics were discussed.

Regional differences in the whistles of Australasian humpback dolphins (genus Sousa)

Characteristics of whistles may be used to study differentiation in dolphins to complement morphological and genetic studies. The whistles of four populations of Chinese humpback dolphins (Sousa chinensis chinensis (Osbeck, 1765)), one population of Taiwanese humpback dolphins (Sousa chinensis taiwanensis Wang, Yang, and Hung, 2015), and one population of Australian humpback dolphins (Sousa sahulensis Jefferson and Rosenbaum, 2014) were compared to determine if differences in whistles support current views of population structure and regional and species differentiation in the genus Sousa Gray, 1866. Acoustic features were extracted from whistles captured by broadband recording systems. Permutational MANOVAs were conducted to test for differences between populations, regions, and species. Random forest trees were also used to classify similar whistles. A significant amount of variation in acoustic features was explained by population (pseudo F[5,2742] = 191.66, p < 0.001), regional (pseudo F[3,2741] = 280.62, p < 0.001), and species (pseudo F[1,999] = 3.7, p < 0.05) differences in humpback dolphin whistles. Random forest trees correctly classified whistles into populations from 40% to 67%, regions from 51% to 80%, and species from 74% to 80%. Differences in whistles were consistent with the current ideas of population-, regional-, and species-level differences within the genus Sousa, based on morphological and genetic data, as well as geographic distance and barriers to movement.

Variation in Metabolite Profiles of Amphibian Skin Bacterial Communities Across Elevations in the Neotropics

Both the structure and function of host-associated microbial communities are potentially impacted by environmental conditions, just as the outcomes of many free-living species interactions are context-dependent. Many amphibian populations have declined around the globe due to the fungal skin pathogen, Batrachochytrium dendrobatidis (Bd), but enivronmental conditions may influence disease dynamics. For instance, in Panamá, the most severe Bd outbreaks have occurred at high elevation sites. Some amphibian species harbor bacterial skin communities that can inhibit the growth of Bd, and therefore, there is interest in understanding whether environmental context could also alter these host-associated microbial communities in a way that might ultimately impact Bd dynamics. In a field survey in Panamá, we assessed skin bacterial communities (16S rRNA amplicon sequencing) and metabolite profiles (HPLC-UV/Vis) of Silverstoneia flotator from three high- and three low-elevation populations representing a range of environmental conditions. Across elevations, frogs had similar skin bacterial communities, although one lowland site appeared to differ. Interestingly, we found that bacterial richness decreased from west to east, coinciding with the direction of Bd spread through Panamá. Moreover, metabolite profiles suggested potential functional variation among frog populations and between elevations. While the frogs have similar bacterial community structure, the local environment might shape the metabolite profiles. Ultimately, host-associated community structure and function could be dependent on environmental conditions, which could ultimately influence host disease susceptibility across sites.

The identification of Oligo-Miocene mammalian palaeocommunities from the Riversleigh World Heritage Area, Australia and an appraisal of palaeoecological techniques

Fourteen of the best sampled Oligo-Miocene local faunas from the Riversleigh World Heritage Area, north-western Queensland, Australia are analysed using classification and ordination techniques to identify potential mammalian palaeocommunities and palaeocommunity types. Abundance data for these faunas are used, for the first time, in conjunction with presence/absence data. An early Miocene Faunal Zone B and two middle Miocene Faunal Zone C palaeocommunities are recognised, as well as one palaeocommunity type. Change in palaeocommunity structure, between the early Miocene and middle Miocene, may be the result of significant climate change during the Miocene Carbon Isotope Excursion. The complexes of local faunas identified will allow researchers to use novel palaeocommunities in future analyses of Riversleigh’s fossil faunas. The utility of some palaeoecological multivariate indices and techniques is examined. The Dice index is found to outperform other binary similarity/distance coefficients, while the UPGMA algorithm is more useful than neighbour joining. Evidence is equivocal for the usefulness of presence/absence data compared to abundance.

Functional structure of the bromeliad tank microbiome is strongly shaped by local geochemical conditions

Phytotelmata in tank-forming Bromeliaceae plants are regarded as potential miniature models for aquatic ecology, but detailed investigations of their microbial communities are rare. Hence, the biogeochemistry in bromeliad tanks remains poorly understood. Here we investigate the structure of bacterial and archaeal communities inhabiting the detritus within the tanks of two bromeliad species, Aechmea nudicaulis and Neoregelia cruenta, from a Brazilian sand dune forest. We used metagenomic sequencing for functional community profiling and 16S sequencing for taxonomic profiling. We estimated the correlation between functional groups and various environmental variables, and compared communities between bromeliad species. In all bromeliads, microbial communities spanned a metabolic network adapted to oxygen-limited conditions, including all denitrification steps, ammonification, sulfate respiration, methanogenesis, reductive acetogenesis and anoxygenic phototrophy. Overall, CO2 reducers dominated in abundance over sulfate reducers, and anoxygenic phototrophs largely outnumbered oxygenic photoautotrophs. Functional community structure correlated strongly with environmental variables, between and within a single bromeliad species. Methanogens and reductive acetogens correlated with detrital volume and canopy coverage, and exhibited higher relative abundances in N. cruenta. A comparison of bromeliads to freshwater lake sediments and soil from around the world, revealed stark differences in terms of taxonomic as well as functional microbial community structure.

Temporal Dynamics of Soil Microbial Communities below the Seedbed under Two Contrasting Tillage Regimes

Agricultural productivity relies on a wide range of ecosystem services provided by the soil biota. Plowing is a fundamental component of conventional farming, but long-term detrimental effects such as soil erosion and loss of soil organic matter have been recognized. Moving towards more sustainable management practices such as reduced tillage or crop residue retention can reduce these detrimental effects, but will also influence structure and function of the soil microbiota with direct consequences for the associated ecosystem services. Although there is increasing evidence that different tillage regimes alter the soil microbiome, we have a limited understanding of the temporal dynamics of these effects. Here, we used high-throughput sequencing of bacterial and fungal ribosomal markers to explore changes in soil microbial community structure under two contrasting tillage regimes (conventional and reduced tillage) either with or without crop residue retention. Soil samples were collected over the growing season of two crops (Vicia faba and Triticum aestivum) below the seedbed (15-20 cm). Tillage, crop and growing stage were significant determinants of microbial community structure, but the impact of tillage showed only moderate temporal dependency. Whereas the tillage effect on soil bacteria showed some temporal dependency and became less strong at later growing stages, the tillage effect on soil fungi was more consistent over time. Crop residue retention had only a minor influence on the community. Six years after the conversion from conventional to reduced tillage, soil moisture contents and nutrient levels were significantly lower under reduced than under conventional tillage. These changes in edaphic properties were related to specific shifts in microbial community structure. Notably, bacterial groups featuring copiotrophic lifestyles or potentially carrying the ability to degrade more recalcitrant compounds were favored under conventional tillage, whereas taxa featuring more oligotrophic lifestyles were more abundant under reduced tillage. Our study found that, under the specific edaphic and climatic context of central Belgium, different tillage regimes created different ecological niches that select for different microbial lifestyles with potential consequences for the ecosystem services provided to the plants and their environment.

Limited engraftment of donor microbiome via one-time fecal microbial transplantation in treated HIV-infected individuals

Many HIV-infected individuals on antiretroviral therapy (ART) exhibit persistent systemic inflammation, which predicts morbidity and mortality. ART-treated subjects concurrently exhibit marked compositional alterations in the gut bacterial microbiota and the degree of dysbiosis correlates with systemic inflammation. Whether interventions to modulate the microbiome can affect systemic inflammation is unknown. An open-label fecal microbial transplantation (FMT) was delivered by colonoscopy to asymptomatic HIV-infected ART-suppressed individuals without antibiotic pre-treatment. Stool was assessed before and after FMT for engraftment of donor microbes, and peripheral blood was assayed for immune activation biomarkers. Six participants received FMT and 2 participants served as controls. No serious adverse effects occurred during 24 weeks of follow-up. At baseline, HIV-infected individuals exhibited microbiota profiles distinct from uninfected donors. During the 8 weeks post-FMT, recipients demonstrated partial engraftment of the donor microbiome (P < 0.05). Recipient microbiota remained significantly distant from donors, unlike that observed following FMT for treatment of C. difficile infection. Systemic inflammatory markers showed no significant change post-FMT. FMT was well-tolerated in ART-treated, HIV-infected individuals. Engraftment was detectable but modest, and appeared to be limited to specific bacterial taxa. Whether antibiotic conditioning can enhance engraftment and the capacity of microbiota to modulate inflammation remains to be investigated.

Bacterial Community Associated with the Reef Coral Mussismilia braziliensis's Momentum Boundary Layer over a Diel Cycle

Corals display circadian physiological cycles, changing from autotrophy during the day to heterotrophy during the night. Such physiological transition offers distinct environments to the microbial community associated with corals: an oxygen-rich environment during daylight hours and an oxygen-depleted environment during the night. Most studies of coral reef microbes have been performed on samples taken during the day, representing a bias in the understanding of the composition and function of these communities. We hypothesized that coral circadian physiology alters the composition and function of microbial communities in reef boundary layers. Here, we analyzed microbial communities associated with the momentum boundary layer (MBL) of the Brazilian endemic reef coral Mussismilia braziliensis during a diurnal cycle, and compared them to the water column. We determined microbial abundance and nutrient concentration in samples taken within a few centimeters of the coral's surface every 6 h for 48 h, and sequenced microbial metagenomes from a subset of the samples. We found that dominant taxa and functions in the coral MBL community were stable over the time scale of our sampling, with no significant shifts between night and day samples. Interestingly, the two water column metagenomes sampled 1 m above the corals were also very similar to the MBL metagenomes. When all samples were analyzed together, nutrient concentration significantly explained 40% of the taxonomic dissimilarity among dominant genera in the community. Functional profiles were highly homogenous and not significantly predicted by any environmental variables measured. Our data indicated that water flow may overrule the effects of coral physiology in the MBL bacterial community, at the scale of centimeters, and suggested that sampling resolution at the scale of millimeters may be necessary to address diurnal variation in community composition.

The potential use of cuticular hydrocarbons and multivariate analysis to age empty puparial cases of Calliphora vicina and Lucilia sericata

Cuticular hydrocarbons (CHC) have been successfully used in the field of forensic entomology for identifying and ageing forensically important blowfly species, primarily in the larval stages. However in older scenes where all other entomological evidence is no longer present, Calliphoridae puparial cases can often be all that remains and therefore being able to establish the age could give an indication of the PMI. This paper examined the CHCs present in the lipid wax layer of insects, to determine the age of the cases over a period of nine months. The two forensically important species examined were Calliphora vicina and Lucilia sericata. The hydrocarbons were chemically extracted and analysed using Gas Chromatography – Mass Spectrometry. Statistical analysis was then applied in the form of non-metric multidimensional scaling analysis (NMDS), permutational multivariate analysis of variance (PERMANOVA) and random forest models. This study was successful in determining age differences within the empty cases, which to date, has not been establish by any other technique.

Predicting the HMA-LMA Status in Marine Sponges by Machine Learning

The dichotomy between high microbial abundance (HMA) and low microbial abundance (LMA) sponges has been observed in sponge-microbe symbiosis, although the extent of this pattern remains poorly unknown. We characterized the differences between the microbiomes of HMA (n = 19) and LMA (n = 17) sponges (575 specimens) present in the Sponge Microbiome Project. HMA sponges were associated with richer and more diverse microbiomes than LMA sponges, as indicated by the comparison of alpha diversity metrics. Microbial community structures differed between HMA and LMA sponges considering Operational Taxonomic Units (OTU) abundances and across microbial taxonomic levels, from phylum to species. The largest proportion of microbiome variation was explained by the host identity. Several phyla, classes, and OTUs were found differentially abundant in either group, which were considered “HMA indicators” and “LMA indicators.” Machine learning algorithms (classifiers) were trained to predict the HMA-LMA status of sponges. Among nine different classifiers, higher performances were achieved by Random Forest trained with phylum and class abundances. Random Forest with optimized parameters predicted the HMA-LMA status of additional 135 sponge species (1,232 specimens) without a priori knowledge. These sponges were grouped in four clusters, from which the largest two were composed of species consistently predicted as HMA (n = 44) and LMA (n = 74). In summary, our analyses shown distinct features of the microbial communities associated with HMA and LMA sponges. The prediction of the HMA-LMA status based on the microbiome profiles of sponges demonstrates the application of machine learning to explore patterns of host-associated microbial communities.

Differences in nutrient content of eggs and larvae as indicators for improvement of broodstock nutrition in walleye (Sander vitreus) production

The aim of this study was to increase our knowledge about the nutritional needs of walleye (Sander vitreus (Mitchill, 1818)) broodstock. Two dietary treatments were tested: the first consisted of frozen Atlantic mackerel (Scomber scombrus L., 1758) pieces with a supplement of experimental dry pellets and the second consisted of frozen Atlantic mackerel pieces with a supplement of commercial trout pellets. Hatching success was highest in eggs from the wild broodstock. Eggs and larvae from wild fish were characterized by the highest levels of linoleic (LOA; 18:2 n–6) and linolenic (LNA; 18:3 n–3) fatty acids in the neutral lipids. They had the lowest level of docosahexaenoic acid (DHA) in both neutral and polar lipids. The results suggest the importance of the amino acids leucine, lysine, methionine, and serine for better egg survival. Vitamin A1 content was similar in eggs from wild fish and the younger broodstock from both treatment groups, but was 10 times lower in eggs from older broodstock. Vitamin A2 was highest in eggs from wild fish. Our results showed that neither of the two diets really resulted in egg or larval quality that were comparable with that of wild fish.

Corexit 9500 Enhances Oil Biodegradation and Changes Active Bacterial Community Structure of Oil-Enriched Microcosms

To better understand the impacts of Corexit 9500 on the structure and activity levels of hydrocarbon-degrading microbial communities, we analyzed next-generation 16S rRNA gene sequencing libraries of hydrocarbon enrichments grown at 5 and 25°C using both DNA and RNA extracts as the sequencing templates. Oil biodegradation patterns in both 5 and 25°C enrichments were consistent with those reported in the literature (i.e., aliphatics were degraded faster than aromatics). Slight increases in biodegradation were observed in the presence of Corexit at both temperatures. Differences in community structure were observed between treatment conditions in the DNA-based libraries. The 25°C consortia were dominated by Vibrio, Idiomarina, Marinobacter, Alcanivorax, and Thalassospira species, while the 5°C consortia were dominated by several species of the genera Flavobacterium, Alcanivorax, and Oleispira Most of these genera have been linked to hydrocarbon degradation and have been observed after oil spills. Colwellia and Cycloclasticus, known aromatic degraders, were also found in these enrichments. The addition of Corexit did not have an effect on the active bacterial community structure of the 5°C consortia, while at 25°C, a decrease in the relative abundance of Marinobacter was observed. At 25°C, Thalassospira, Marinobacter, and Idiomarina were present at higher relative abundances in the RNA than DNA libraries, suggesting that they were active in degradation. Similarly, Oleispira was greatly stimulated by the addition of oil at 5°C.IMPORTANCE While dispersants such as Corexit 9500 can be used to treat oil spills, there is still debate on the effectiveness on enhancing oil biodegradation and its potential toxic effect on oil-degrading microbial communities. The results of this study provide some insights on the microbial dynamics of hydrocarbon-degrading bacterial populations in the presence of Corexit 9500. Operational taxonomic unit (OTU) analyses indicated that several OTUs were inhibited by the addition of Corexit. Conversely, a number of OTUs were stimulated by the addition of the dispersant, many of which were identified as known hydrocarbon-degrading bacteria. The results highlight the value of using RNA-based methods to further understand the impact of dispersant on the overall activity of different hydrocarbon-degrading bacterial groups.

Calcium decline reduces population growth rates of zooplankton in field mesocosms

Regional calcium (Ca) decline, a legacy of acid deposition and logging, is a potential threat to aquatic organisms. Lake surveys and laboratory studies indicate that Ca-rich daphniids are likely most susceptible, allowing for competitive release of other taxa with low Ca demand. Indeed, dramatic shifts in zooplankton community structure have been documented in lakes where Ca has declined, amid multiple other stressors. Given the perceived threat of this large-scale stressor, manipulative studies are needed to evaluate causal relationships between Ca decline and zooplankton community structure. We analysed per capita growth rates of zooplankton from three independent mesocosm experiments where we manipulated aqueous Ca concentrations to reflect current and future Ca concentrations. In two experiments where Ca concentration was reduced to 0.6 or 0.9 mg/L, we observed reduced growth rates for several taxa, including daphniids, bosminids, and copepods. No effect of Ca was detected in the experiment where Ca concentrations ranged from 1.2 to 2.5 mg/L, a gradient representing 68% of lakes in south-central Ontario. These results suggest that future Ca decline in soft-water Canadian Shield lakes may be accompanied by shifts in community structure and overall declines in zooplankton production.

How and where to perform biomonitoring studies: different levels of toxic metal pollution are detected in the Alagados Reservoir in Southern Brazil

The Alagados Reservoir located in the Brazilian city of Ponta Grossa is used to supply water for human consumption but is impacted by toxic metals. The current study combined chemical, biochemical, and multivariate analyses to determine the bioavailability of toxic metals at three sites (AL-A, AL-B, and AL-C) within the Alagados Reservoir. Metal bioaccumulation was analyzed in the liver, gills, and muscle tissue of a native fish species (Geophagus brasiliensis), and neurotoxicity, xenobiotic metabolism, and oxidative stress were evaluated using biochemical biomarkers. Additionally, histopathological studies were performed on the gills and the liver using scanning electron microscopy (SEM) and conventional light microscopy (LM), respectively. Overall, the bioaccumulation of metals, biomarkers of oxidative stress, and melanomacrophage counts indicate that the AL-C and the AL-A sites are the most and least affected by metals, respectively. The AL-B site presented the lowest acetylcholinesterase enzyme activity, a finding which was probably associated with the agricultural activities around this area of the reservoir. The biomarkers clearly revealed that toxic metals negatively affect all three sites studied herein and that human activity is the major source of pollutants. Despite the existence of different pollution levels within the Alagados Reservoir, it is still used as a human water supply.

Ericoid Roots and Mycospheres Govern Plant-Specific Bacterial Communities in Boreal Forest Humus

In this study, the bacterial populations of roots and mycospheres of the boreal pine forest ericoid plants, heather (Calluna vulgaris), bilberry (Vaccinium myrtillus), and lingonberry (Vaccinium vitis-idaea), were studied by qPCR and next-generation sequencing (NGS). All bacterial communities of mycosphere soils differed from soils uncolonized by mycorrhizal mycelia. Colonization by mycorrhizal hyphae increased the total number of bacterial 16S ribosomal DNA (rDNA) gene copies in the humus but decreased the number of different bacterial operational taxonomic units (OTUs). Nevertheless, ericoid roots and mycospheres supported numerous OTUs not present in uncolonized humus. Bacterial communities in bilberry mycospheres were surprisingly similar to those in pine mycospheres but not to bacterial communities in heather and lingonberry mycospheres. In contrast, bacterial communities of ericoid roots were more similar to each other than to those of pine roots. In all sample types, the relative abundances of bacterial sequences belonging to Alphaproteobacteria and Acidobacteria were higher than the sequences belonging to other classes. Soil samples contained more Actinobacteria, Deltaproteobacteria, Opitutae, and Planctomycetia, whereas Armatimonadia, Betaproteobacteria, Gammaproteobacteria, and Sphingobacteriia were more common to roots. All mycosphere soils and roots harbored bacteria unique to that particular habitat. Our study suggests that the habitation by ericoid plants increases the overall bacterial diversity of boreal forest soils.

Effects of Gelling Agent and Extracellular Signaling Molecules on the Culturability of Marine Bacteria

Only 1% of marine bacteria are currently culturable using standard laboratory procedures, and this is a major obstacle for our understanding of the biology of marine microorganisms and for the discovery of novel microbial natural products. Therefore, the purpose of this study was to investigate if improved cultivation conditions, including the use of an alternative gelling agent and supplementation with signaling molecules, improve the culturability of bacteria from seawater. Replacing agar with gellan gum improved viable counts 3- to 40-fold, depending on medium composition and incubation conditions, with a maximum of 6.6% culturability relative to direct cell counts. Through V4 amplicon sequencing we found that culturable diversity was also affected by a change in gelling agent, facilitating the growth of orders not culturable on agar-based substrates. Community analyses showed that communities grown on gellan gum substrates were significantly different from communities grown on agar and that they covered a larger fraction of the seawater community. Other factors, such as incubation temperature and time, had less obvious effects on viable counts and culturable diversity. Supplementation with acylated homoserine lactones (AHLs) did not have a positive effect on total viable counts or a strong effect on culturable diversity. However, low concentrations of AHLs increased the relative abundance of sphingobacteria. Hence, with alternative growth substrates, it is possible to significantly increase the number and diversity of cultured marine bacteria.IMPORTANCE Serious challenges to human health, such as the occurrence and spread of antibiotic resistance and an aging human population in need of bioactive pharmaceuticals, have revitalized the search for natural microbial products. The marine environment, representing the largest ecosystem in the biosphere, harbors an immense and virtually untapped microbial diversity producing unique bioactive compounds. However, we are currently able to cultivate only a minute fraction of this diversity. The lack of cultivated microbes is hampering not only bioprospecting efforts but also our general understanding of marine microbes. In this study, we present a means to increase the number and diversity of cultured bacteria from seawater, showing that relatively simple changes to medium components may facilitate the isolation and growth of hitherto unknown bacterial orders.

Endosymbiont interference and microbial diversity of the Pacific coast tick, Dermacentor occidentalis, in San Diego County, California

The Pacific coast tick, Dermacentor occidentalis Marx, is found throughout California and can harbor agents that cause human diseases such as anaplasmosis, ehrlichiosis, tularemia, Rocky Mountain spotted fever and rickettsiosis 364D. Previous studies have demonstrated that nonpathogenic endosymbiotic bacteria can interfere with Rickettsia co-infections in other tick species. We hypothesized that within D. occidentalis ticks, interference may exist between different nonpathogenic endosymbiotic or nonendosymbiotic bacteria and Spotted Fever group Rickettsia (SFGR). Using PCR amplification and sequencing of the rompA gene and intergenic region we identified a cohort of SFGR-infected and non-infected D. occidentalis ticks collected from San Diego County. We then amplified a partial segment of the 16S rRNA gene and used next-generation sequencing to elucidate the microbiomes and levels of co-infection in the ticks. The SFGR R. philipii str. 364D and R. rhipicephali were detected in 2.3% and 8.2% of the ticks, respectively, via rompA sequencing. Interestingly, next generation sequencing revealed an inverse relationship between the number of Francisella-like endosymbiont (FLE) 16S rRNA sequences and Rickettsia 16S rRNA sequences within individual ticks that is consistent with partial interference between FLE and SFGR infecting ticks. After excluding the Rickettsia and FLE endosymbionts from the analysis, there was a small but significant difference in microbial community diversity and a pattern of geographic isolation by distance between collection locales. In addition, male ticks had a greater diversity of bacteria than female ticks and ticks that weren’t infected with SFGR had similar microbiomes to canine skin microbiomes. Although experimental studies are required for confirmation, our findings are consistent with the hypothesis that FLEs and, to a lesser extent, other bacteria, interfere with the ability of D. occidentalis to be infected with certain SFGR. The results also raise interesting possibilities about the effects of putative vertebrate hosts on the tick microbiome.

Multiple environmental changes drive forest floor vegetation in a temperate mountain forest

Human-induced changes of the environment and their possible impacts on temperate forest understory plant communities have been examined in many studies. However, the relative contribution of individual environmental factors to these changes in the herb layer is still unclear. In this study, we used vegetation survey data covering a time period of 21 years and collected from 143 permanent plots in the Northern Limestone Alps, Austria. Data on soil chemistry (49 plots), light condition (51 plots), soil temperature and moisture (four and six plots), disturbance (all plots), climate (one station in a clearing area), and airborne sulfur (S) and nitrogen (N) deposition (two forest stands) were available for analyses. We used these data together with plot mean Ellenberg indicator values in a path analysis to attribute their relative contributions to observed vegetation changes. Our analysis reveals a strong directional shift of the forest understory plant community. We found strong evidence for a recovery of the ground-layer vegetation from acidification as response to decreased S deposition. We did not observe a community response to atmospheric N deposition, but we found a response to altered climatic conditions (thermophilization and drying). The path analysis revealed that changes in the light regime, which were related to small-scale disturbances, had most influence on herb layer community shifts. Thermophilization and drying were identified as drivers of understory community changes independent of disturbance events.

Microbial Community Functional Potential and Composition Are Shaped by Hydrologic Connectivity in Riverine Floodplain Soils

Riverine floodplains are ecologically and economically valuable ecosystems that are heavily threatened by anthropogenic stressors. Microbial communities in floodplain soils mediate critical biogeochemical processes, yet we understand little about the relationship between these communities and variation in hydrologic connectivity related to land management or topography. Here, we present metagenomic evidence that differences among microbial communities in three floodplain soils correspond to a long-term gradient of hydrologic connectivity. Specifically, all strictly anaerobic taxa and metabolic pathways were positively associated with increased hydrologic connectivity and flooding frequency. In contrast, most aerobic taxa and all strictly aerobic pathways were negatively related to hydrologic connectivity and flooding frequency. Furthermore, the genetic potential to metabolize organic compounds tended to decrease as hydrologic connectivity increased, which may reflect either the observed concomitant decline of soil organic matter or the parallel increase in both anaerobic taxa and pathways. A decline in soil N, accompanied by an increased genetic potential for oligotrophic N acquisition subsystems, suggests that soil nutrients also shape microbial communities in these soils. We conclude that differences among floodplain soil microbial communities can be conceptualized along a gradient of hydrologic connectivity. Additionally, we show that these differences are likely due to connectivity-related variation in flooding frequency, soil organic matter, and soil N. Our findings are particularly relevant to the restoration and management of microbially mediated biogeochemical processes in riverine floodplain wetlands.

A model for paired-multinomial data and its application to analysis of data on a taxonomic tree

In human microbiome studies, sequencing reads data are often summarized as counts of bacterial taxa at various taxonomic levels specified by a taxonomic tree. This article considers the problem of analyzing two repeated measurements of microbiome data from the same subjects. Such data are often collected to assess the change of microbial composition after certain treatment, or the difference in microbial compositions across body sites. Existing models for such count data are limited in modeling the covariance structure of the counts and in handling paired multinomial count data. A new probability distribution is proposed for paired-multinomial count data, which allows flexible covariance structure and can be used to model repeatedly measured multivariate count data. Based on this distribution, a test statistic is developed for testing the difference in compositions based on paired multinomial count data. The proposed test can be applied to the count data observed on a taxonomic tree in order to test difference in microbiome compositions and to identify the subtrees with different subcompositions. Simulation results indicate that proposed test has correct type 1 errors and increased power compared to some commonly used methods. An analysis of an upper respiratory tract microbiome data set is used to illustrate the proposed methods.

Bacterial Community Composition Associated with Pyrogenic Organic Matter (Biochar) Varies with Pyrolysis Temperature and Colonization Environment

Pyrogenic organic matter (PyOM) is widely distributed in soil and fluvial ecosystems and plays an important role in biogeochemical cycling. Many studies have reported changes in soil microbial communities stimulated by PyOM, but very little is known about the microbial communities associated with PyOM. The microbes that colonize PyOMs can participate in the mineralization of PyOM, so changing its structure affects the fate of PyOMs and contributes to soil biogeochemical cycling. This study identified the bacterial community composition associated with PyOMs on the basis of high-throughput sequencing and demonstrated that both PyOM pyrolysis temperature and the colonization environment determined the bacterial community composition. Our work increases our understanding of the dominant phylogenetic taxa associated with PyOMs, demonstrates mechanisms mediating microbial metabolism and growth in PyOMs, and expands a new research area for pyrogenic organic matter. This study identified the bacterial community composition associated with PyOM, which is widely distributed in the environment. Most bacterial OTUs preferentially thrived on PyOM pyrolyzed at low temperature, while some specific OTUs thrived on PyOM pyrolyzed at high temperature.

Microbes that colonize pyrogenic organic matter (PyOM) (also called biochar) play an important role in PyOM mineralization and crucially affect soil biogeochemical cycling, while the microbial community composition associated with PyOM particles is poorly understood. We generated two manure-based PyOMs with different characteristics (PyOM pyrolyzed at the low temperature of 300°C [i.e., PyOM300] and at the high temperature of 700°C [i.e., PyOM700]) and added them to high-carbon (4.15%) and low-C (0.37%) soil for microbial colonization. 16S rRNA gene sequencing showed that Actinobacteria, particularly Actinomycetales, was the dominant taxon in PyOM, regardless of the PyOM pyrolysis temperature and soil type. Bacterial communities associated with PyOM particles from high-C soils were similar to those in non-PyOM-amended soils. PyOM300 had higher total microbial activity and more differential bacterial communities than PyOM700. More bacterial operational taxonomic units (OTUs) preferentially thrived on the low-pyrolysis-temperature PyOM, while some specific OTUs thrived on high-pyrolysis-temperature PyOM. In particular, Chloroflexi species tended to be more prevalent in high-pyrolysis-temperature PyOM in low-C soils. In conclusion, the differences in colonized bacterial community composition between the different PyOMs were strongly influenced by the pyrolysis temperatures of PyOM, i.e., under conditions of easily mineralizable C or fused aromatic C, and by other properties, e.g., pH, surface area, and nutrient content.

IMPORTANCE Pyrogenic organic matter (PyOM) is widely distributed in soil and fluvial ecosystems and plays an important role in biogeochemical cycling. Many studies have reported changes in soil microbial communities stimulated by PyOM, but very little is known about the microbial communities associated with PyOM. The microbes that colonize PyOMs can participate in the mineralization of PyOM, so changing its structure affects the fate of PyOMs and contributes to soil biogeochemical cycling. This study identified the bacterial community composition associated with PyOMs on the basis of high-throughput sequencing and demonstrated that both PyOM pyrolysis temperature and the colonization environment determined the bacterial community composition. Our work increases our understanding of the dominant phylogenetic taxa associated with PyOMs, demonstrates mechanisms mediating microbial metabolism and growth in PyOMs, and expands a new research area for pyrogenic organic matter. This study identified the bacterial community composition associated with PyOM, which is widely distributed in the environment. Most bacterial OTUs preferentially thrived on PyOM pyrolyzed at low temperature, while some specific OTUs thrived on PyOM pyrolyzed at high temperature.

Mineral Ecology: Surface Specific Colonization and Geochemical Drivers of Biofilm Accumulation, Composition, and Phylogeny

This study tests the hypothesis that surface composition influences microbial community structure and growth of biofilms. We used laboratory biofilm reactors (inoculated with a diverse subsurface community) to explore the phylogenetic and taxonomic variability in microbial communities as a function of surface type (carbonate, silicate, aluminosilicate), media pH, and carbon and phosphate availability. Using high-throughput pyrosequencing, we found that surface type significantly controlled ~70–90% of the variance in phylogenetic diversity regardless of environmental pressures. Consistent patterns also emerged in the taxonomy of specific guilds (sulfur-oxidizers/reducers, Gram-positives, acidophiles) due to variations in media chemistry. Media phosphate availability was a key property associated with variation in phylogeny and taxonomy of whole reactors and was negatively correlated with biofilm accumulation and α-diversity (species richness and evenness). However, mineral-bound phosphate limitations were correlated with less biofilm. Carbon added to the media was correlated with a significant increase in biofilm accumulation and overall α-diversity. Additionally, planktonic communities were phylogenetically distant from those in biofilms. All treatments harbored structurally (taxonomically and phylogenetically) distinct microbial communities. Selective advantages within each treatment encouraged growth and revealed the presence of hundreds of additional operational taxonomix units (OTU), representing distinct consortiums of microorganisms. Ultimately, these results provide evidence that mineral/rock composition significantly influences microbial community structure, diversity, membership, phylogenetic variability, and biofilm growth in subsurface communities.

Random sampling causes the low reproducibility of rare eukaryotic OTUs in Illumina COI metabarcoding

DNA metabarcoding, the PCR-based profiling of natural communities, is becoming the method of choice for biodiversity monitoring because it circumvents some of the limitations inherent to traditional ecological surveys. However, potential sources of bias that can affect the reproducibility of this method remain to be quantified. The interpretation of differences in patterns of sequence abundance and the ecological relevance of rare sequences remain particularly uncertain. Here we used one artificial mock community to explore the significance of abundance patterns and disentangle the effects of two potential biases on data reproducibility: indexed PCR primers and random sampling during Illumina MiSeq sequencing. We amplified a short fragment of the mitochondrial Cytochrome c Oxidase Subunit I (COI) for a single mock sample containing equimolar amounts of total genomic DNA from 34 marine invertebrates belonging to six phyla. We used seven indexed broad-range primers and sequenced the resulting library on two consecutive Illumina MiSeq runs. The total number of Operational Taxonomic Units (OTUs) was ∼4 times higher than expected based on the composition of the mock sample. Moreover, the total number of reads for the 34 components of the mock sample differed by up to three orders of magnitude. However, 79 out of 86 of the unexpected OTUs were represented by <10 sequences that did not appear consistently across replicates. Our data suggest that random sampling of rare OTUs (e.g., small associated fauna such as parasites) accounted for most of variation in OTU presence–absence, whereas biases associated with indexed PCRs accounted for a larger amount of variation in relative abundance patterns. These results suggest that random sampling during sequencing leads to the low reproducibility of rare OTUs. We suggest that the strategy for handling rare OTUs should depend on the objectives of the study. Systematic removal of rare OTUs may avoid inflating diversity based on common β descriptors but will exclude positive records of taxa that are functionally important. Our results further reinforce the need for technical replicates (parallel PCR and sequencing from the same sample) in metabarcoding experimental designs. Data reproducibility should be determined empirically as it will depend upon the sequencing depth, the type of sample, the sequence analysis pipeline, and the number of replicates. Moreover, estimating relative biomasses or abundances based on read counts remains elusive at the OTU level.

Diurnal Dynamics of Gaseous and Dissolved Metabolites and Microbiota Composition in the Bovine Rumen

Diurnal patterns of ruminal fermentation metabolites and microbial communities are not commonly assessed when investigating variation in ruminal CH4 production. The aims of this study were to monitor diurnal patterns of: (i) gaseous and dissolved metabolite concentrations in the bovine rumen, (ii) H2 and CH4 emitted, and (iii) the rumen microbiota. Furthermore, the effect of dietary inclusion of linseed oil on these patterns was assessed. Four rumen cannulated multiparous cows were used in a cross-over design with two 17 days periods and two dietary treatments: a control diet and a linseed oil supplemented diet [40% maize silage, 30% grass silage, 30% concentrate on dry matter (DM) basis for both diets; fat contents of 33 vs. 56 g/kg of DM]. On day 11, rumen contents were sampled for 10 h after morning feeding to profile gaseous and dissolved metabolite concentrations and microbiota composition. H2 and CH4 emission (mass per unit of time) was measured in respiration chambers from day 13 to 17. A 100-fold increase in ruminal H2 partial pressure (contribution to the total pressure of rumen headspace gases) was observed at 0.5 h after feeding. This peak was followed by a decline to basal level. Qualitatively similar patterns after feeding were also observed for H2 and CH4 emission, ethanol and lactate concentrations, and propionate molar proportion, although the opposite pattern was seen for acetate molar proportion. Associated with these patterns, a temporal biphasic change in the microbial composition was observed as based on 16S ribosomal RNA with certain taxa specifically associated with each phase. Bacterial concentrations (log10 16S ribosomal RNA gene copies based) were affected by time, and were increased by linseed oil supplementation. Archaeal concentrations (log10 16S ribosomal RNA gene copies based) tended to be affected by time and were not affected by diet, despite linseed oil supplementation decreasing CH4 emission, tending to decrease the partial pressure of CH4, and tending to increase propionate molar proportion. Linseed oil supplementation affected microbiota composition, and was most associated with an uncultivated Bacteroidales taxon. In summary, our findings support the importance of diurnal dynamics for the understanding of VFA, H2, and CH4 production.

An Endohyphal Bacterium (Chitinophaga, Bacteroidetes) Alters Carbon Source Use by Fusarium keratoplasticum (F. solani Species Complex, Nectriaceae)

Bacterial endosymbionts occur in diverse fungi, including members of many lineages of Ascomycota that inhabit living plants. These endosymbiotic bacteria (endohyphal bacteria, EHB) often can be removed from living fungi by antibiotic treatment, providing an opportunity to assess their effects on functional traits of their fungal hosts. We examined the effects of an endohyphal bacterium (Chitinophaga sp., Bacteroidetes) on substrate use by its host, a seed-associated strain of the fungus Fusarium keratoplasticum, by comparing growth between naturally infected and cured fungal strains across 95 carbon sources with a Biolog® phenotypic microarray. Across the majority of substrates (62%), the strain harboring the bacterium significantly outperformed the cured strain as measured by respiration and hyphal density. These substrates included many that are important for plant- and seed-fungus interactions, such as D-trehalose, myo-inositol, and sucrose, highlighting the potential influence of EHB on the breadth and efficiency of substrate use by an important Fusarium species. Cases in which the cured strain outperformed the strain harboring the bacterium were observed in only 5% of substrates. We propose that additive or synergistic substrate use by the fungus-bacterium pair enhances fungal growth in this association. More generally, alteration of the breadth or efficiency of substrate use by dispensable EHB may change fungal niches in short timeframes, potentially shaping fungal ecology and the outcomes of fungal-host interactions.

Hard time to be parents? Sea urchin fishery shifts potential reproductive contribution of population onto the shoulders of the young adults

Background

In Sardinia, as in other regions of the Mediterranean Sea, sustainable fisheries of the sea urchin Paracentrotus lividus have become a necessity. At harvesting sites, the systematic removal of large individuals (diameter ≥ 50 mm) seriously compromises the biological and ecological functions of sea urchin populations. Specifically, in this study, we compared the reproductive potential of the populations from Mediterranean coastal areas which have different levels of sea urchin fishing pressure. The areas were located at Su Pallosu Bay, where pressure is high and Tavolara-Punta Coda Cavallo, a marine protected area where sea urchin harvesting is low.

Methods

Reproductive potential was estimated by calculating the gonadosomatic index (GSI) from June 2013 to May 2014 both for individuals of commercial size (diameter without spines, TD ≥ 50 mm) and the undersized ones with gonads (30 ≤ TD < 40 mm and 40 ≤ TD < 50 mm). Gamete output was calculated for the commercial-size class and the undersized individuals with fertile gonads (40 ≤ TD < 50 mm) in relation to their natural density (gamete output per m²).

Results

The reproductive potential of populations was slightly different at the beginning of the sampling period but it progressed at different rates with an early spring spawning event in the high-pressure zone and two gamete depositions in early and late spring in the low-pressure zone. For each fertile size class, GSI values changed significantly during the year of our study and between the two zones. Although the multiple spawning events determined a two-fold higher total gamete output of population (popTGO) in the low-pressure zone, the population mean gamete output (popMGO) was similar in the two zones. In the high-pressure zone, the commercial-sized individuals represented approximatively 5% of the population, with almost all the individuals smaller than 60 mm producing an amount of gametes nearly three times lower than the undersized ones. Conversely, the high density of the undersized individuals released a similar amount of gametes to the commercial-size class in the low-pressure zone.

Discussion

Overall, the lack of the commercial-size class in the high-pressure zone does not seem to be very alarming for the self-supporting capacity of the population, and the reproductive potential contribution seems to depend more on the total density of fertile sea urchins than on their size. However, since population survival in the high-pressure zone is supported by the high density of undersized sea urchins between 30 and 50 mm, management measures should be addressed to maintain these sizes and to shed light on the source of the larval supply.

Sorghum and wheat differentially affect caecal microbiota and associated performance characteristics of meat chickens

This study compared the effects of wheat- and sorghum-based diets on broiler chickens. The growth performance and caecal microbial community of chickens were measured and correlations between productivity and specific gut microbes were observed. Cobb broilers 15 days of age were individually caged and two dietary treatments were used, one with a wheat-based diet (n = 48) and another one with a sorghum-based diet (n = 48). Growth performance measurements were taken over a 10 day period and samples for microbiota analysis were taken at the end of that period. Caecal microbiota was characterised by sequencing of 16S bacterial rRNA gene amplicons. Overall, the results indicated that a sorghum-based diet produced higher apparent metabolisable energy (AME) and body-weight gain (BWG) values in chickens, compared to a wheat-based diet. Nevertheless, sorghum-fed birds had higher feed conversion ratio (FCR) values than wheat-fed birds, possibly because of some anti-nutritional factors in sorghum. Further analyses showed that caecal microbial community was significantly associated with AME values, but microbiota composition differed between dietary treatments. A number of bacteria were individually correlated with growth performance measurements. Numerous OTUs assigned to strains of Lactobacillus crispatus and Lachnospiraceae, which were prevalent in sorghum-fed chickens, were correlated with high AME and BWG values, respectively. Additionally, a number of OTUs assigned to Clostridiales that were prevalent in wheat-fed chickens were correlated with low FCR values. Overall, these results suggest that between-diet variations in growth performance were partly associated with changes in the caecal microbiota.

Effect of predatory bacteria on the gut bacterial microbiota in rats

Bdellovibrio bacteriovorus and Micavibrio aeruginosavorus are Gram-negative proteobacteria that are obligate predators of other Gram-negative bacteria and are considered potential alternatives to antibiotics. Most studies focusing on predatory bacteria have been performed in vitro, thus the effect of predatory bacteria on a live host, including the impact on the ecology of the native microbiota, has yet to be fully examined. In this study, intrarectal inoculations of Sprague-Dawley rats with predatory bacteria were performed. Additionally, feces were collected for seven days post-inoculation to determine the effect on gut bacterial diversity. Rat colonic tissue exhibited no abnormal histopathological effects due to predatory bacteria. A modest increase in pro-inflammatory cytokines was measured in the colons of rats inoculated with predatory bacteria by 24 and 48 hours, with all but IL-13 returning to baseline by seven days. V4 16S rRNA gene sequencing of fecal DNA demonstrated minimal shifts in taxonomic representation over the week due to predatory bacteria. Changes in bacterial populations due to exposure to B. bacteriovorus are predicted to contribute to health, however, an overgrowth of Prevotella was observed due to exposure to M. aeruginosavorus. This study further addresses safety concerns associated with the potential use of predatory bacteria to treat infections.

Normalization and microbial differential abundance strategies depend upon data characteristics

BACKGROUND

Data from 16S ribosomal RNA (rRNA) amplicon sequencing present challenges to ecological and statistical interpretation. In particular, library sizes often vary over several ranges of magnitude, and the data contains many zeros. Although we are typically interested in comparing relative abundance of taxa in the ecosystem of two or more groups, we can only measure the taxon relative abundance in specimens obtained from the ecosystems. Because the comparison of taxon relative abundance in the specimen is not equivalent to the comparison of taxon relative abundance in the ecosystems, this presents a special challenge. Second, because the relative abundance of taxa in the specimen (as well as in the ecosystem) sum to 1, these are compositional data. Because the compositional data are constrained by the simplex (sum to 1) and are not unconstrained in the Euclidean space, many standard methods of analysis are not applicable. Here, we evaluate how these challenges impact the performance of existing normalization methods and differential abundance analyses.

RESULTS

Effects on normalization: Most normalization methods enable successful clustering of samples according to biological origin when the groups differ substantially in their overall microbial composition. Rarefying more clearly clusters samples according to biological origin than other normalization techniques do for ordination metrics based on presence or absence. Alternate normalization measures are potentially vulnerable to artifacts due to library size. Effects on differential abundance testing: We build on a previous work to evaluate seven proposed statistical methods using rarefied as well as raw data. Our simulation studies suggest that the false discovery rates of many differential abundance-testing methods are not increased by rarefying itself, although of course rarefying results in a loss of sensitivity due to elimination of a portion of available data. For groups with large (~10×) differences in the average library size, rarefying lowers the false discovery rate. DESeq2, without addition of a constant, increased sensitivity on smaller datasets (<20 samples per group) but tends towards a higher false discovery rate with more samples, very uneven (~10×) library sizes, and/or compositional effects. For drawing inferences regarding taxon abundance in the ecosystem, analysis of composition of microbiomes (ANCOM) is not only very sensitive (for >20 samples per group) but also critically the only method tested that has a good control of false discovery rate.

CONCLUSIONS

These findings guide which normalization and differential abundance techniques to use based on the data characteristics of a given study.

Combined magnetic resonance elastography and collagen molecular magnetic resonance imaging accurately stage liver fibrosis in a rat model

Hepatic fibrosis is associated with an overproduction of matrix proteins and a pathological increase of liver stiffness. Noninvasive magnetic resonance (MR) quantification of matrix can be assessed with a collagen-binding molecular MR probe and stiffness by MR elastography, complementary techniques. This study used both imaging techniques to more accurately stage hepatic fibrosis in a rat model. Thirty rats with varying levels of diethylnitrosamine-induced liver fibrosis were imaged before and 45 minutes after injection of collagen-specific probe EP-3533. MR elastography was performed in the same imaging session. Changes in liver relaxation rate post-EP-3533 and liver stiffness were compared to the collagen proportional area determined by histology and to Ishak scoring using receiver operating characteristic analysis. Collagen imaging was most sensitive to early fibrosis, while elastography was more sensitive to advanced fibrosis. This complementary feature enabled the formulation of a composite model using multivariate analysis of variance. This model incorporated the discriminating advantages of both MR techniques, resulting in more accurate staging throughout fibrotic progression.

CONCLUSION

Collagen molecular MR imaging is complementary to MR elastography, and combining the two techniques in a single exam leads to increased diagnostic accuracy for all stages of fibrosis. (Hepatology 2017;65:1015-1025).