Fungal community composition in the gut of rove beetles (Coleoptera: Staphylinidae) from the Canadian boreal forest reveals possible endosymbiotic interactions for dietary needs

The goal of this study was to investigate the fungal community composition in the gut of Staphylinidae from boreal forest in order to better understand the diversity and the complexity of fungus-insect relationships. DNA gut content analyses of nine abundant rove beetle species (Coleoptera, Staphylinidae) living in the boreal balsam fir forest ecosystem (Montmorency Forest, Quebec, Canada) were performed to identify the fungal taxa present either as endosymbiotic taxa or as a source of nutrition. A total of 42 fungal operational taxonomic units (OTUs) were recorded from the analysis of 441 fungal ITS rDNA sequences recovered from gut extracts. The OTU richness per species ranged between four in Tachinus quebecensis and 16 in Atheta ventricosa. The fungal mycobiota in posterior gut extracts was dominated by Saccharomycetales (12 OTUs), followed by Sordariomycetes (nine OTUs). No significant difference was observed between the OTU richness recorded within each of the three subfamilies of rove beetles investigated. The core mycobiome of the posterior gut extracts was dominated by three OTUs related to yeasts, with ITS sequences having pairwise similarities equal to or greater than 99% with Candida mesenterica, Debaryomyces spp. and Ophiostoma pluriannulatum. These results provide some evidence of the consumer-resource relationships of these beetles. Predominance of yeast and fungal spores in the posterior gut of rove beetles suggests that they may play an important role in their dietary requirements and as endosymbionts.     

Comparison of the diversity,composition, and host recurrence of xylariaceous endophytes in subtropical,cool temperate,and subboreal regions in Japan

The diversity, composition, and host recurrence of endophytic fungi in the Xylariaceae were compared in subtropical (ST), cool temperate (CT), and subboreal forests (SB) in Japan based on the 28S ribosomal DNA sequences from fungal isolates. A total of 610 isolates were obtained from the leaves of 167 tree species in three sites, which were classified into 42 operational taxonomic units (OTUs) at the 99 % similarity level of the 28S rDNA sequence. ST, CT, and SB yielded 31, 13, and three OTUs, respectively. The OTU richness, diversity, and evenness of fungal communities were in the order: ST > CT > SB. The 42 OTUs were assigned to nine genera in the Xylariaceae: Xylaria, Annulohypoxylon, Anthostomella, Biscogniauxia, Nemania, Hypoxylon, Muscodor, Daldinia, and Rosellinia. Xylarioid isolates in the subfamily Xylarioideae outnumbered Hypoxyloid isolates in the subfamily Hypoxyloideae in ST and CT, whereas the opposite was found in SB. Sørensen’s quotient of similarity was generally low between the three sites. Host recurrence of fungal OTUs was evaluated with the degree of specialization of interaction network between xylariaceous endophytes and plant species and compared between the three sites. We found that the networks in the three sites showed a significantly higher degree of specialization than simulated networks, where partners were associated randomly. Permutational multivariate analyses of variance indicated that plant family and leaf trait significantly affected the OTU composition in ST, which can account for the specialization of interaction network and host recurrence of xylariaceous endophytes.     

Mucosa-Associated Faecalibacterium prausnitzii Phylotype Richness Is Reduced in Patients with Inflammatory Bowel Disease

Faecalibacterium prausnitzii depletion in intestinal diseases has been extensively reported, but little is known about intraspecies variability. This work aims to determine if subjects with gastrointestinal disease host mucosa-associated F. prausnitzii populations different from those hosted by healthy individuals. A new species-specific PCR-denaturing gradient gel electrophoresis (PCR-DGGE) method targeting the 16S rRNA gene was developed to fingerprint F. prausnitzii populations in biopsy specimens from 31 healthy control (H) subjects and 36 Crohn''s disease (CD), 23 ulcerative colitis (UC), 6 irritable bowel syndrome (IBS), and 22 colorectal cancer (CRC) patients. The richness of F. prausnitzii subtypes was lower in inflammatory bowel disease (IBD) patients than in H subjects. The most prevalent operational taxonomic units (OTUs) consisted of four phylotypes (OTUs with a 99% 16S rRNA gene sequence similarity [OTU99]), which were shared by all groups of patients. Their distribution and the presence of some disease-specific F. prausnitzii phylotypes allowed us to differentiate the populations in IBD and CRC patients from that in H subjects. At the level of a minimum similarity of 97% (OTU97), two phylogroups accounted for 98% of the sequences. Phylogroup I was found in 87% of H subjects but in under 50% of IBD patients (P = 0.003). In contrast, phylogroup II was detected in >75% of IBD patients and in only 52% of H subjects (P = 0.005). This study reveals that even though the main members of the F. prausnitzii population are present in both H subjects and individuals with gut diseases, richness is reduced in the latter and an altered phylotype distribution exists between diseases. This approach may serve as a basis for addressing the suitability of F. prausnitzii phylotypes to be quantified as a putative biomarker of disease and depicting the importance of the loss of these subtypes in disease pathogenesis.     

Deciphering the Bacterial Microbiome of Citrus Plants in Response to ‘Candidatus Liberibacter asiaticus’-Infection and Antibiotic Treatments

The bacterial microbiomes of citrus plants were characterized in response to ‘Candidatus Liberibacter asiaticus’ (Las)-infection and treatments with ampicillin (Amp) and gentamicin (Gm) by Phylochip-based metagenomics. The results revealed that 7,407 of over 50,000 known Operational Taxonomic Units (OTUs) in 53 phyla were detected in citrus leaf midribs using the PhyloChip™ G3 array, of which five phyla were dominant, Proteobacteria (38.7%), Firmicutes (29.0%), Actinobacteria (16.1%), Bacteroidetes (6.2%) and Cyanobacteria (2.3%). The {"type":"entrez-protein","attrs":{"text":"OTU62806","term_id":"1193547623","term_text":"OTU62806"}}OTU62806, representing ‘Candidatus Liberibacter’, was present with a high titer in the plants graft-inoculated with Las-infected scions treated with Gm at 100 mg/L and in the water-treated control (CK₁). However, the Las bacterium was not detected in the plants graft-inoculated with Las-infected scions treated with Amp at 1.0 g/L or in plants graft-inoculated with Las-free scions (CK₂). The PhyloChip array demonstrated that more OTUs, at a higher abundance, were detected in the Gm-treated plants than in the other treatment and the controls. Pairwise comparisons indicated that 23 OTUs from the Achromobacter spp. and 12 OTUs from the Methylobacterium spp. were more abundant in CK₂ and CK₁, respectively. Ten abundant OTUs from the Stenotrophomonas spp. were detected only in the Amp-treatment. These results provide new insights into microbial communities that may be associated with the progression of citrus huanglongbing (HLB) and the potential effects of antibiotics on the disease and microbial ecology.     

Introducing DOTUR, a Computer Program for Defining Operational Taxonomic Units and Estimating Species Richness

Although copious qualitative information describes the members of the diverse microbial communities on Earth, statistical approaches for quantifying and comparing the numbers and compositions of lineages in communities are lacking. We present a method that addresses the challenge of assigning sequences to operational taxonomic units (OTUs) based on the genetic distances between sequences. We developed a computer program, DOTUR, which assigns sequences to OTUs by using either the furthest, average, or nearest neighbor algorithm for each distance level. DOTUR uses the frequency at which each OTU is observed to construct rarefaction and collector's curves for various measures of richness and diversity. We analyzed 16S rRNA gene libraries derived from Scottish and Amazonian soils and the Sargasso Sea with DOTUR, which assigned sequences to OTUs rapidly and reliably based on the genetic distances between sequences and identified previous inconsistencies and errors in assigning sequences to OTUs. An analysis of the two 16S rRNA gene libraries from soil demonstrated that they do not contain enough sequences to support a claim that they contain different numbers of bacterial lineages with statistical confidence (P > 0.05), nor do they contain enough sequences to provide a robust estimate of species richness when an OTU is defined as containing sequences that are no more than 3% different from each other. In contrast, the richness of OTUs at the 3% level in the Sargasso Sea collection began to plateau after the sampling of 690 sequences. We anticipate that an equivalent extent of sampling for soil would require sampling more than 10,000 sequences, almost 100 times the size of typical sequence collections obtained from soil.     

Influence of antimicrobial feed additives on broiler commensal posthatch gut microbiota development and performance

The effects of avilamycin, zinc bacitracin, and flavophospholipol on broiler gut microbial community colonization and bird performance in the first 17 days posthatch were investigated. Significant differences in gut microbiota associated with gut section, dietary treatment, and age were identified by terminal restriction fragment length polymorphism (T-RFLP), although no performance-related differences between dietary treatments were detected. Similar age-related shifts in the gut microbiota were identified regardless of diet but varied between the ilea and ceca. Interbird variabilities in ileal bacterial communities were reduced (3 to 7 days posthatch) in chicks fed with feed containing antimicrobial agents. Avilamycin and flavophospholipol had the most consistent effect on gut microbial communities. Operational taxonomic units (OTU) linked to changes in gut microbiota in birds on antimicrobial-supplemented diets were characterized and identified. Some OTUs could be identified to the species level; however, the majority could be only tentatively classified to the genus, family, order, or domain level. OTUs 140 to 146 (Lachnospiraceae), OTU 186/188 (Lactobacillus johnsonii), OTU 220 (Lachnospiraceae), OTUs 284 to 288 (unclassified bacterial spp. or Ruminococcaceae), OTU 296/298 (unclassified bacterium or Clostridiales), and OTU 480/482 (Oxalobacteraceae) were less prevalent in the guts of chicks fed antimicrobial-supplemented diets. OTU 178/180 (Lactobacillus crispatus), OTU 152 (Lactobacillus reuteri or unclassified Clostridiales), OTU 198/200 (Subdoligranulum spp.), and OTU 490/492 (unclassified bacterium or Enterobacteriaceae) were less prevalent in the gut of chicks raised on the antimicrobial-free diet. The identification of key bacterial species influenced by antimicrobial-supplemented feed immediately posthatch may assist in the formulation of diets that facilitate beneficial gut microbial colonization and, hence, the development of alternatives to current antimicrobial agents in feed for sustainable poultry production.     

Biochar as a feed supplement for nutrient digestibility and growth performance of Catla catla fingerlings

The purpose of the current research was to determine the impact of various biochar supplements on nutrient digestibility and growth performance of Catla catla fingerlings fed on Moringa oleifera seed meal (MOSM) based diet. An experiment with 90 days of feeding was conducted to investigate the efficacy of biochar obtained from different sources such as parthenium, farmyard manure, poultry waste, vegetable waste, and corncob waste at 2 mg/kg. There were15 fingerlings in each tank in the triplicate set of tanks. Fingerlings were fed at the rate of 5 % of their live wet weight. From each tank, feces were collected twice daily and stored to determine nutrient digestibility. Results showed that poultry waste biochar (test diet-IV) at 2 mg/kg was the best source to significantly (P < 0.05) improve weight gain % (256.58 %) and feed conversion ratio (1.19) than the other biochar sources and control diet. Furthermore, in terms of nutrient digestibility, optimum fat (81.90 %), protein (75.92 %), and gross energy (74.84 kcalg^?1) values were obtained by using the same type of biochar (poultry waste). Conclusively, among all the five biochar sources, poultry waste biochar proved to be the best one, improving fish body performance.     

Characterization and identification of productivity-associated rhizobacteria in wheat

The rhizosphere is populated by a numerous and diverse array of rhizobacteria, and many impact productivity in largely unknown ways. Here we characterize the rhizobacterial community in a wheat variety categorized according to shoot biomass using 16S rRNA pyrosequencing abundance data. Plants were grown in homogenized field soil under greenhouse conditions, and DNA was extracted and pyrosequenced, resulting in 29,007 quality sequences. Operational taxonomic units (OTUs) that were significantly associated with biomass productivity were identified using an exact test adjusted for the false-discovery rate. The productivity deviation expressed as a percentage of the total mean square for regression (PMSR) was determined for each OTU. Out of 719 OTUs, 42 showed significant positive associations and 39 showed significant negative associations (q value, ≤0.05). OTUs with the greatest net positive associations, by genus, were as follows: Duganella, OTU 43 and OTU 3; Janthinobacterium, OTU 278; Pseudomonas, OTU 588; and Cellvibrio, OTU 1847. Those with negative associations were as follows: Bacteria, OTU 273; Chryseobacterium, OTU 508; Proteobacteria, OTU 249; and Enterobacter, OTU 357. Shoot biomass productivity was strongly correlated with the balance between the overall abundances of positive- and negative-productivity-associated OTUs. High-productivity rhizospheres contained 9.2 significant positives for every negatively associated rhizobacterium, while low-productivity rhizospheres showed 2.3 significant negatives for every positively associated rhizobacterium. Overall rhizobacterial community diversity as measured by the Chao1, Shannon, and Simpson indexes was nonlinearly related to productivity, closely fitting a wavelike cubic equation. We conclude that shoot biomass productivity is strongly related to the ratio of positive- to negative-productivity-associated rhizobacteria in the rhizosphere. This study identifies significant OTUs composing the productive and unproductive rhizobacterial communities.     

Microbial diversity in ostrich ceca as revealed by 16S ribosomal RNA gene clone library and detection of novel Fibrobacter species

The ostrich (Struthio camelus) is a herbivorous bird and although the hindgut is known as the site for fiber digestion, little is known about the microbial diversity in the ostrich hindgut. Our aim was to analyze the microbial diversity in ostrich ceca using a 16S ribosomal RNA gene (rDNA) clone library approach. A total of 310 clones were sequenced and phylogenetically analyzed and were classified into 110 operational taxonomy units (OTUs) based on a 98% similarity criterion. The similarity of the sequences ranged from 86 to 99% and 95 OTUs had less than 98% similarity to the sequences in the public databases. Coverage and the Shannon–Wiener index (H′) of the library were 83.9% and 4.29, respectively. The sequences were assigned to the following 6 phyla: Firmicutes (50.9% of the total number of sequences), Bacteroidetes (39.4%), Fibrobacteres (6.5%), Euryarchaeota (1.9%), Spirochaetes (1.0%), and Verrucomicrobia (0.3%); approximately 90% of the sequences were affiliated with Firmicutes and Bacteroidetes. The only OTU of Fibrobacteres (OTU 107), had 93 and 90% similarity to Fibrobacter succinogenes and F. intestinalis, respectively, suggesting a new species of Fibrobacter in ostrich ceca. Clostridium coccoides and C. leptum formed major groups within the Firmicutes. There was no OTU with high similarity (≥98%) to the 16S rDNA of cultivated fibrolytic bacteria in our library. Although two OTUs were affiliated with Euryarchaeota, no sequence was affiliated with methanogenic Archaea. This study presents the very complex ostrich cecal microbial community, in which the majority of the bacterial species have not yet been cultivated.     

Microbial Contaminants of Cord Blood Units Identified by 16S rRNA Sequencing and by API Test System,and Antibiotic Sensitivity Profiling

Over a period of ten months a total of 5618 cord blood units (CBU) were screened for microbial contamination under routine conditions. The antibiotic resistance profile for all isolates was also examined using ATB strips. The detection rate for culture positive units was 7.5%, corresponding to 422 samples.16S rRNA sequence analysis and identification with API test system were used to identify the culturable aerobic, microaerophilic and anaerobic bacteria from CBUs. From these samples we recovered 485 isolates (84 operational taxonomic units, OTUs) assigned to the classes Bacteroidia, Actinobacteria, Clostridia, Bacilli, Betaproteobacteria and primarily to the Gammaproteobacteria. Sixty-nine OTUs, corresponding to 447 isolates, showed 16S rRNA sequence similarities above 99.0% with known cultured bacteria. However, 14 OTUs had 16S rRNA sequence similarities between 95 and 99% in support of genus level identification and one OTU with 16S rRNA sequence similarity of 90.3% supporting a family level identification only. The phenotypic identification formed 29 OTUs that could be identified to the species level and 9 OTUs that could be identified to the genus level by API test system. We failed to obtain identification for 14 OTUs, while 32 OTUs comprised organisms producing mixed identifications. Forty-two OTUs covered species not included in the API system databases. The API test system Rapid ID 32 Strep and Rapid ID 32 E showed the highest proportion of identifications to the species level, the lowest ratio of unidentified results and the highest agreement to the results of 16S rRNA assignments. Isolates affiliated to the Bacilli and Bacteroidia showed the highest antibiotic multi-resistance indices and microorganisms of the Clostridia displayed the most antibiotic sensitive phenotypes.     

Fungal disease and temperature alter skin microbiome structure in an experimental salamander system

Pathogens compete with host microbiomes for space and resources. Their shared environment impacts pathogen–microbiome–host interactions, which can lead to variation in disease outcome. The skin microbiome of red‐backed salamanders (Plethodon cinereus) can reduce infection by the pathogen Batrachochytrium dendrobatidis (Bd) at moderate infection loads, with high species richness and high abundance of competitors as putative mechanisms. However, it is unclear if the skin microbiome can reduce epizootic Bd loads across temperatures. We conducted a laboratory experiment to quantify skin microbiome and host responses (P. cinereus: n = 87) to Bd at mimicked epizootic loads across temperatures (13, 17 and 21°C). We quantified skin microbiomes using 16S rRNA gene metabarcoding and identified operational taxonomic units (OTUs) taxonomically similar to culturable bacteria known to kill Bd (anti‐Bd OTUs). Prior to pathogen exposure, temperature changed the microbiome (OTU richness decreased by 12% and the abundance of anti‐Bd OTUs increased by 18% per degree increase in temperature), but these changes were not predictive of disease outcome. After exposure, Bd changed the microbiome (OTU richness decreased by 0.1% and the abundance of anti‐Bd OTUs increased by 0.2% per 1% increase in Bd load) and caused high host mortality across temperatures (35/45: 78%). Temperature indirectly impacted microbiome change and mortality through its direct effect on pathogen load. We did not find support for the microbiome impacting Bd load or host survival. Our research reveals complex host, pathogen, microbiome and environmental interactions to demonstrate that during epizootic events the microbiome will be unlikely to reduce pathogen invasion, even for putatively Bd‐resistant species.     

Polymerase matters: non-proofreading enzymes inflate fungal community richness estimates by up to 15 %

Rare taxa overwhelm metabarcoding data generated using next-generation sequencing (NGS). Low frequency Operational Taxonomic Units (OTUs) may be artifacts generated by PCR-amplification errors resulting from polymerase mispairing. We analyzed two Internal Transcribed Spacer 2 (ITS2) MiSeq libraries generated with proofreading (ThermoScientific Phusion^®) and non-proofreading (ThermoScientific Phire^®) polymerases from the same MiSeq reaction, the same samples, using the same DNA tags, and with two different clustering methods to evaluate the effect of polymerase and clustering tool choices on the estimates of richness, diversity and community composition. Our data show that, while the overall communities are comparable, OTU richness is exaggerated by the use of the non-proofreading polymerase–up to 15 % depending on the clustering method, and on the threshold of low frequency OTU removal. The overestimation of richness also consistently led to underestimation of community evenness, a result of increase in the low frequency OTUs. Stringent thresholds of eliminating the rare reads remedy this issue; exclusion of reads that occurred ≤10 times reduced overestimated OTU numbers to <0.3 %. As a result of these findings, we strongly recommend the use of proofreading polymerases to improve the data integrity as well as the use of stringent culling thresholds for rare sequences to minimize overestimation of community richness.     

Next Generation Sequencing to Define Prokaryotic and Fungal Diversity in the Bovine Rumen

A combination of Sanger and 454 sequences of small subunit rRNA loci were used to interrogate microbial diversity in the bovine rumen of 12 cows consuming a forage diet. Observed bacterial species richness, based on the V1–V3 region of the 16S rRNA gene, was between 1,903 to 2,432 species-level operational taxonomic units (OTUs) when 5,520 reads were sampled per animal. Eighty percent of species-level OTUs were dominated by members of the order Clostridiales, Bacteroidales, Erysipelotrichales and unclassified TM7. Abundance of Prevotella species varied widely among the 12 animals. Archaeal species richness, also based on 16S rRNA, was between 8 and 13 OTUs, representing 5 genera. The majority of archaeal OTUs (84%) found in this study were previously observed in public databases with only two new OTUs discovered. Observed rumen fungal species richness, based on the 18S rRNA gene, was between 21 and 40 OTUs with 98.4–99.9% of OTUs represented by more than one read, using Good’s coverage. Examination of the fungal community identified numerous novel groups. Prevotella and Tannerella were overrepresented in the liquid fraction of the rumen while Butyrivibrio and Blautia were significantly overrepresented in the solid fraction of the rumen. No statistical difference was observed between the liquid and solid fractions in biodiversity of archaea and fungi. The survey of microbial communities and analysis of cross-domain correlations suggested there is a far greater extent of microbial diversity in the bovine rumen than previously appreciated, and that next generation sequencing technologies promise to reveal novel species, interactions and pathways that can be studied further in order to better understand how rumen microbial community structure and function affects ruminant feed efficiency, biofuel production, and environmental impact.     

Statistical approaches for estimating actinobacterial diversity in marine sediments

Bacterial diversity in a deep-sea sediment was investigated by constructing actinobacterium-specific 16S ribosomal DNA (rDNA) clone libraries from sediment sections taken 5 to 12, 15 to 18, and 43 to 46 cm below the sea floor at a depth of 3,814 m. Clones were placed into operational taxonomic unit (OTU) groups with >/= 99% 16S rDNA sequence similarity; the cutoff value for an OTU was derived by comparing 16S rRNA homology with DNA-DNA reassociation values for members of the class Actinobacteria. Diversity statistics were used to determine how the level of dominance, species richness, and genetic diversity varied with sediment depth. The reciprocal of Simpson's index (1/D) indicated that the pattern of diversity shifted toward dominance from uniformity with increasing sediment depth. Nonparametric estimation of the species richness in the 5- to 12-, 15- to 18-, and 43- to 46-cm sediment sections revealed a trend of decreasing species number with depth, 1,406, 308, and 212 OTUs, respectively. Application of the LIBSHUFF program indicated that the 5- to 12-cm clone library was composed of OTUs significantly (P = 0.001) different from those of the 15- to 18- and 43- to 46-cm libraries. F(ST) and phylogenetic grouping of taxa (P tests) were both significant (P < 0.00001 and P < 0.001, respectively), indicating that genetic diversity decreased with sediment depth and that each sediment community harbored unique phylogenetic lineages. It was also shown that even nonconservative OTU definitions result in severe underestimation of species richness; unique phylogenetic clades detected in one OTU group suggest that OTUs do not correspond to real ecological groups sensu Palys (T. Palys, L. K. Nakamura, and F. M. Cohan, Int. J. Syst. Bacteriol. 47:1145-1156, 1997). Mechanisms responsible for diversity and their implications are discussed.     

Diversity and dynamics of fungal endophytes in leaves, stems and roots of Stellera chamaejasme L. in northwestern China

This study was conducted to explore fungal endophyte communities inhabiting a toxic weed (Stellera chamaejasme L.) from meadows of northwestern China. The effects of plant tissue and growth stage on endophyte assemblages were characterized. Endophytes were recovered from 50 % of the samples, with a total of 714 isolates. 41 operational taxonomical units (OTUs) were identified, consisting of 40 OTUs belonging primarily to Ascomycota and 1 OTU belonging to Basidiomycota. Pleosporales and Hypocreales were the orders contributing the most species to the endophytic assemblages. The total colonization frequency and species richness of endophytic fungi were higher in roots than in leaves and stems. In addition, for the plant tissues, the structure of fungal communities differed significantly by growth stages of leaf emergence and dormancy; for the plant growth stages, the structure of fungal communities differed significantly by plant tissues. This study demonstrates that S. chamaejasme serves as a reservoir for a wide variety of fungal endophytes that can be isolated from various plant tissues.     

Fungicide Effects on Fungal Community Composition in the Wheat Phyllosphere

The fungicides used to control diseases in cereal production can have adverse effects on non-target fungi, with possible consequences for plant health and productivity. This study examined fungicide effects on fungal communities on winter wheat leaves in two areas of Sweden. High-throughput 454 sequencing of the fungal ITS2 region yielded 235 operational taxonomic units (OTUs) at the species level from the 18 fields studied. It was found that commonly used fungicides had moderate but significant effect on fungal community composition in the wheat phyllosphere. The relative abundance of several saprotrophs was altered by fungicide use, while the effect on common wheat pathogens was mixed. The fungal community on wheat leaves consisted mainly of basidiomycete yeasts, saprotrophic ascomycetes and plant pathogens. A core set of six fungal OTUs representing saprotrophic species was identified. These were present across all fields, although overall the difference in OTU richness was large between the two areas studied.     

Methanogen diversity and community composition in peatlands of the central to northern Appalachian Mountain region, North America

Methanogenic archaea are ubiquitous in peat soils; however, their diversity and distributions within and among peatland ecosystems are not well known. We used comprehensive clone libraries of 16S rRNA gene sequences to investigate spatial patterns in diversity (richness, evenness of taxa) and composition (taxonomic, phylogenetic) of the methanogenic community in six peatlands arrayed 775?km from eastern Ontario, Canada to West Virginia, USA. Five sites were Sphagnum (moss) and shrub dominated; one site was sedge dominated; and, potential rates of methane (CH₄) production ranged from 15 to 450?nmol/g?day. The gradient allowed us to examine influences of site conditions, site history, and climate on community composition. The region had representatives of methanogens from four taxonomic orders. We observed 29 operationally defined units (OTUs) based on >97% sequence identity. One OTU accounted for 43% of all clones, whereas 15 OTUs were rare with <1% of the total number of clones. The number of OTUs per site ranged from 4 to 21, and statistical analysis suggested diversity of 4–43 per site. Eighteen of the OTUs were endemic to one site; albeit, most endemics occurred in the sedge dominated site. One OTU was cosmopolitan, occurring in all six sites. We found a positive relationship between methanogen diversity and rates of CH₄ production per site (Pearson r?=?0.93). Turnover in community composition between sites was weakly related to geographic distance between sites, whereas variation in soil pH and annual temperature played larger roles. About 50% of the variation in community composition was unexplained by distance, pH, mean climate, and site age. We conclude that methanogen diversity in peatlands of the central Appalachian region is shaped by present-day environmental conditions, suggesting an influence of impending climatic and environmental changes.     

Plant identity strongly structures the root-associated fungal community in a diverse subtropical forest

Revealing the relationship between plants and root-associated fungi is very important in understanding diversity maintenance and community assembly in ecosystems. However, the community assembly of root-associated fungi of focal plant species along a subtropical plant species diversity gradient is less documented. Here, we examined root-associated fungal communities associated with five ectomycorrhizal (EM) plant species (Betula luminifera, Castanea henryi, Castanopsis fargesii, C. sclerophylla, and Quercus serrate) in a Chinese subtropical woody plant species diversity (1, 2, 4, 8, 16 and 24 species) experiment, using paired-end Illumina MiSeq sequencing of the ITS2 region. In total, we detected 1933 root-associated fungal operational taxonomic units (OTUs) at a 97% sequence similarity level. Plant identity had a significant effect on total and saprotrophic fungal OTU richness, but plant species diversity level had a significant effect on saprotrophic and pathogenic fungal OTU richness. The community composition of total, saprotrophic and EM fungi was structured by plant identity and plant species diversity level. However, the community composition of pathogenic fungi was only shaped by plant identity. This study highlights that plant identity has a stronger effect on the root-associated fungal community than plant species diversity level in a diverse subtropical forest ecosystem.     

Global patterns of diversity and community structure in marine bacterioplankton

Because of their small size, great abundance and easy dispersal, it is often assumed that marine planktonic microorganisms have a ubiquitous distribution that prevents any structured assembly into local communities. To challenge this view, marine bacterioplankton communities from coastal waters at nine locations distributed world-wide were examined through the use of comprehensive clone libraries of 16S ribosomal RNA genes, used as operational taxonomic units (OTU). Our survey and analyses show that there were marked differences in the composition and richness of OTUs between locations. Remarkably, the global marine bacterioplankton community showed a high degree of endemism, and conversely included few cosmopolitan OTUs. Our data were consistent with a latitudinal gradient of OTU richness. We observed a positive relationship between the relative OTU abundances and their range of occupation, i.e. cosmopolitans had the largest population sizes. Although OTU richness differed among locations, the distributions of the major taxonomic groups represented in the communities were analogous, and all local communities were similarly structured and dominated by a few OTUs showing variable taxonomic affiliations. The observed patterns of OTU richness indicate that similar evolutionary and ecological processes structured the communities. We conclude that marine bacterioplankton share many of the biogeographical and macroecological features of macroscopic organisms. The general processes behind those patterns are likely to be comparable across taxa and major global biomes.     

Divergence thresholds and divergent biodiversity estimates: can metabarcoding reliably describe zooplankton communities?

DNA metabarcoding is a promising method for describing communities and estimating biodiversity. This approach uses high‐throughput sequencing of targeted markers to identify species in a complex sample. By convention, sequences are clustered at a predefined sequence divergence threshold (often 3%) into operational taxonomic units (OTUs) that serve as a proxy for species. However, variable levels of interspecific marker variation across taxonomic groups make clustering sequences from a phylogenetically diverse dataset into OTUs at a uniform threshold problematic. In this study, we use mock zooplankton communities to evaluate the accuracy of species richness estimates when following conventional protocols to cluster hypervariable sequences of the V4 region of the small subunit ribosomal RNA gene (18S) into OTUs. By including individually tagged single specimens and “populations” of various species in our communities, we examine the impact of intra‐ and interspecific diversity on OTU clustering. Communities consisting of single individuals per species generated a correspondence of 59–84% between OTU number and species richness at a 3% divergence threshold. However, when multiple individuals per species were included, the correspondence between OTU number and species richness dropped to 31–63%. Our results suggest that intraspecific variation in this marker can often exceed 3%, such that a single species does not always correspond to one OTU. We advocate the need to apply group‐specific divergence thresholds when analyzing complex and taxonomically diverse communities, but also encourage the development of additional filtering steps that allow identification of artifactual rRNA gene sequences or pseudogenes that may generate spurious OTUs.