Operational Taxonomic Units Sequences Research Articles

Intratumor microbiome features reveal antitumor potentials of intrahepatic cholangiocarcinoma

ABSTRACT Intrahepatic cholangiocarcinoma (ICC) is a rare malignancy with a high prevalence in China. This study aimed to characterize the ICC tissues’ bacterial metagenomics signature and explore its antitumor potential for cancer. In this study, 16S rRNA sequencing was carried out on 99 tissues to characterize the features of intratumoral microbiota, followed by single-cell RNA sequencing (scRNA-seq) and multilevel validation. The presence of microbial DNA in tissues was determined using staining, fluorescence in situ hybridization (FISH), and transmission electron microscopy (TEM). A Gram-positive aerobic bacterium, identified as Staphylococcus capitis, was cultured from fresh tissues. Meanwhile, scRNA-seq showed that intratumoral bacteria could be present in multiple cell types. Using 16S rRNA sequencing, we identified a total of 2,320,287 high-quality reads corresponding to 4,594 OTU (operational taxonomic units) sequences. The most abundant bacterial orders include Burkholderiales, Pseudomonadales, Xanthomonadales, Bacillales and Clostridiales. Alpha and Beta diversity analysis revealed specific features in different tissues. In addition, the content of Paraburkholderia fungorum was significantly higher in the paracancerous tissues and negatively correlated with CA199 (Carbohydrate antigen199) levels. The results of in vitro and in vivo experiments suggest that P. fungorum possesses an antitumor activity against tumors. Metabolomics and transcriptomics showed that P. fungorum could inhibit tumor growth through alanine, aspartate and glutamate metabolism. We determined the characteristic profile of the intratumoral microbiota and the antitumor effect of P. fungorum in ICC.

Intraspecies characterization of bacteria via evolutionary modeling of protein domains

The ability to detect and characterize bacteria within a biological sample is crucial for the monitoring of infections and epidemics, as well as for the study of human health and its relationship with commensal microorganisms. To this aim, a commonly used technique is the 16S rRNA gene targeted sequencing. PCR-amplified 16S sequences derived from the sample of interest are usually clustered into the so-called Operational Taxonomic Units (OTUs) based on pairwise similarities. Then, representative OTU sequences are compared with reference (human-made) databases to derive their phylogeny and taxonomic classification. Here, we propose a new reference-free approach to define the phylogenetic distance between bacteria based on protein domains, which are the evolving units of proteins. We extract the protein domain profiles of 3368 bacterial genomes and we use an ecological approach to model their Relative Species Abundance distribution. Based on the model parameters, we then derive a new measurement of phylogenetic distance. Finally, we show that such model-based distance is capable of detecting differences between bacteria in cases in which the 16S rRNA-based method fails, providing a possibly complementary approach , which is particularly promising for the analysis of bacterial populations measured by shotgun sequencing.

Bacterial composition and colony structure of the lower respiratory tract in infants and children with recurrent wheezing: a case–control study

BackgroundThe bacterial load of the human lower respiratory tract is at least several times lower than that of the other parts of the body. This study aimed to identify the bacterial composition and colony structure of the lower respiratory tract in infants and children with recurrent wheezing compared with those of children with a bronchial foreign body and clarify whether the length of wheezing in infants can contribute to differences in the lower respiratory tract’s bacterial colony structure.MethodsWe collected specimens of alveolar lavage fluid from 48 infants and children who underwent fiberoptic bronchoscopy and were divided into groups: A1 (multiple wheezing: wheezing more than three times in < 1 month), A2 (persistent wheezing: > 1 month), and B (bronchial foreign body; control group). We analyzed the bacterial community structure of alveolar lavage fluid using high-throughput sequencing. The richness and diversity of the microbial communities were assessed by α and β diversity analyses.ResultsA total of 6,644 operational taxonomic units (OTUs) were obtained based on the Illumina Nova sequencing platform and clustered according to those that met the 97% identity threshold, followed by species annotation of the OTU sequences. In the annotation results, 2,608 (39.25%) OTUs were annotated at the genus level. At the genus level, Sphingomonas and Phyllobacterium were significantly higher in group A1 than in group B. There were significantly more Phyllobacterium in group A2 than in group B. Prevotella, Neisseria, and Haemophilus were higher in group B than in groups A1 and A2. The differences in the between-group α and β diversity analyses were statistically significant. The microbial diversity in groups A1 and A2 was significantly less than that in group B, but there was no statistical difference in bacterial community diversity between groups A1 and A2.ConclusionRecurrent wheezing in infants and children is more likely due to alterations in the overall bacterial microecology and disruption of host respiration and immune homeostasis than the effects of a single bacterium.

Metabarcoding of insect-associated fungal communities: a comparison of internal transcribed spacer (ITS) and large-subunit (LSU) rRNA markers.

Full taxonomic characterisation of fungal communities is necessary for establishing ecological associations and early detection of pathogens and invasive species. Complex communities of fungi are regularly characterised by metabarcoding using the Internal Transcribed Spacer (ITS) and the Large-Subunit (LSU) gene of the rRNA locus, but reliance on a single short sequence fragment limits the confidence of identification. Here we link metabarcoding from the ITS2 and LSU D1-D2 regions to characterise fungal communities associated with bark beetles (Scolytinae), the likely vectors of several tree pathogens. Both markers revealed similar patterns of overall species richness and response to key variables (beetle species, forest type), but identification against the respective reference databases using various taxonomic classifiers revealed poor resolution towards lower taxonomic levels, especially the species level. Thus, Operational Taxonomic Units (OTUs) could not be linked via taxonomic classifiers across ITS and LSU fragments. However, using phylogenetic trees (focused on the epidemiologically important Sordariomycetes) we placed OTUs obtained with either marker relative to reference sequences of the entire rRNA cistron that includes both loci and demonstrated the largely similar phylogenetic distribution of ITS and LSU-derived OTUs. Sensitivity analysis of congruence in both markers suggested the biologically most defensible threshold values for OTU delimitation in Sordariomycetes to be 98% for ITS2 and 99% for LSU D1-D2. Studies of fungal communities using the canonical ITS barcode require corroboration across additional loci. Phylogenetic analysis of OTU sequences aligned to the full rRNA cistron shows higher success rate and greater accuracy of species identification compared to probabilistic taxonomic classifiers.

Revealing microbial species diversity using sequence capture by hybridization.

Targeting small parts of the 16S rDNA phylogenetic marker by metabarcoding reveals microorganisms of interest but cannot achieve a taxonomic resolution at the species level, precluding further precise characterizations. To identify species behind operational taxonomic units (OTUs) of interest, even in the rare biosphere, we developed an innovative strategy using gene capture by hybridization. From three OTU sequences detected upon polyphenol supplementation and belonging to the rare biosphere of the human gut microbiota, we revealed 59 nearly full-length 16S rRNA genes, highlighting high bacterial diversity hidden behind OTUs while evidencing novel taxa. Inside each OTU, revealed 16S rDNA sequences could be highly distant from each other with similarities down to 85 %. We identified one new family belonging to the order Clostridiales , 39 new genera and 52 novel species. Related bacteria potentially involved in polyphenol degradation have also been identified through genome mining and our results suggest that the human gut microbiota could be much more diverse than previously thought.

Different analysis strategies of 16S rRNA gene data from rodent studies generate contrasting views of gut bacterial communities associated with diet, health and obesity.

BackgroundOne of the main functions of diet is to nurture the gut microbiota and this relationship affects the health of the host. However, different analysis strategies can generate different views on the relative abundance of each microbial taxon, which can affect our conclusions about the significance of diet to gut health in lean and obese subjects. Here we explored the impact of using different analysis strategies to study the gut microbiota in a context of diet, health and obesity.MethodsOver 15 million 16S rRNA gene sequences from published studies involving dietary interventions in obese laboratory rodents were analyzed. Three strategies were used to assign the 16S sequences to Operational Taxonomic Units (OTUs) based on the GreenGenes reference OTU sequence files clustered at 97% and 99% similarity.ResultsDifferent strategies to select OTUs influenced the relative abundance of all bacterial taxa, but the magnitude of this phenomenon showed a strong study effect. Different taxa showed up to 20% difference in relative abundance within the same study, depending on the analysis strategy. Very few OTUs were shared among the samples. ANOSIM test on unweighted UniFrac distances showed that study, sequencing technique, animal model, and dietary treatment (in that order) were the most important factors explaining the differences in bacterial communities. Except for obesity status, the contribution of diet and other factors to explain the variability in bacterial communities was lower when using weighted UniFrac distances. Predicted functional profile and high-level phenotypes of the microbiota showed that each study was associated with unique features and patterns.ConclusionsThe results confirm previous findings showing a strong study effect on gut microbial composition and raise concerns about the impact of analysis strategies on the membership and composition of the gut microbiota. This study may be helpful to guide future research aiming to investigate the relationship between diet, health, and the gut microbiota.

Association of Tongue Bacterial Flora and Subtypes of Liver-Fire Hyperactivity Syndrome in Hypertensive Patients.

Structural changes in symbiotic human microorganisms can affect host phenotype. Liver-fire hyperactivity syndrome (LFHS) presents as bitter taste, halitosis, xerostomia, odontalgia, and other oral symptoms. LFHS is associated with hypertension (EH). In this study, tongue flora was analyzed to further understand the intrinsic relationship between tongue flora and LFHS. Samples of tongue coating, from 16 patients with EH-LFHS, 16 with EH-non-LFHS, and 16 controls, were obtained; then, 16S rRNA variable (V3-V4) regions were amplified and sequenced by MiSeq PE300 Sequencing. Tag clustering and Operational Taxonomic Units (OTUs) abundance analysis were used to compare the OTU sequence with the 16S database. The species were classified, and diversity and structure of the bacterial flora were compared between the three groups. Alpha diversity analysis, including Observed Species index and Chao index, indicated significantly higher richness of species in patients with EH-LFHS (p < 0.05). Higher phylogenetic diversity, in patients with EH-non-LFHS, indicates greater differences in evolutionary history than in patients with EH-LFHS. Streptococcus, Rothia, Neisseria, and Sphingomonas were the most prevalent in patients with EH-LFHS, differed from the other two groups. This indicates that richer bacterial diversity, and structure associated with EH-LFHS, may affect the occurrence, development, and outcome of hypertension and syndrome subtypes recognized by TCM.

The Coal Seam Microbiome (CSMB) reference set, a lingua franca for the microbial coal-to-methane community

Microbial generation of methane from organic matter in coal represents an opportunity to improve environmental and economic outcomes of coal seam gas production. This has prompted research interest in microbial communities responsible for its generation. While advances have been made in understanding the overall structure of these communities and in methods to stimulate methane generation through nutrient additions, many questions still remain as to the metabolic roles played by the majority of the taxa in coal seams. One confounding issue faced by the coal seam microbiology community is an inability to compare OTUs (Operational Taxonomic Units) from 16S rDNA community surveys between studies. The aim of this study is to provide a common language for researchers studying coal seam microbes so that any given OTU from a coal seam community can be readily referenced to OTUs from other studies. To meet this aim a new Coal Seam MicroBiome (CSMB) reference set of OTU sequences was constructed from both new data obtained through community surveys of 28 coal seam formation water samples across three eastern Australian coal basins, along with all compatible, publicly available community surveys. To illustrate the power of the CSMB set for enabling comparisons the set was used to compare the co-occurrence of OTUs from publicly available compatible coal seam community surveys. Additionally, 16S rDNA sequence data from all compatible coal seam isolates and/or reconstructed genomes from published studies were mapped to the CSMB set. Resultant analysis revealed that a large number of published surveys have shared CSMB OTUs and a number of isolate and reconstructed genomes matched CSMB OTUs observed in surveys of the coal seam formation waters across the USA, Canada, China, India, Japan and Australia. This report details these observations of co-occurrence, as well as providing information on downloading, mapping and adding new sequences to the CSMB set.

PrimerMiner: an r package for development and in silico validation of DNA metabarcoding primers

Summary DNA metabarcoding is a powerful tool to assess biodiversity by amplifying and sequencing a standardized gene marker region. Its success is often limited due to variable binding sites that introduce amplification biases. Thus, the development of optimized primers for communities or taxa under study in a certain geographic region and/or ecosystems is of critical importance. However, no tool for obtaining and processing of reference sequence data in bulk that can serve as a backbone for primer design is currently available. We developed the r package PrimerMiner, which batch downloads DNA barcode gene sequences from BOLD and NCBI data bases for specified target taxonomic groups and then applies sequence clustering into operational taxonomic units (OTUs) to reduce biases introduced by the different number of available sequences per species. Additionally, PrimerMiner offers functionalities to evaluate primers in silico, which are in our opinion more realistic than the strategy employed in another available software for that purpose, ecoPCR. We used PrimerMiner to download cytochrome c oxidase subunit I (COI) sequences for 15 important freshwater invertebrate groups, relevant for ecosystem assessment. By processing COI markers from both data bases, we were able to increase the amount of reference data 249‐fold on average, compared to using complete mitochondrial genomes alone. Furthermore, we visualized the generated OTU sequence alignments and describe how to evaluate primers in silico using PrimerMiner. With PrimerMiner, we provide a useful tool to obtain relevant sequence data for targeted primer development and evaluation. The OTU‐based reference alignments generated with PrimerMiner can be used for manual primer design or processed with bioinformatic tools for primer development.

Amplicon-based metabarcoding reveals temporal response of soil microbial community to fumigation-derived products

The use of soil fumigation products to manage soilborne pathogens raises the question of whether it has undesirable effects on the soil ecosystem. Therefore, the aim of this study was to investigate the impacts of alternative fumigation products on soil microbial population structure. We characterised soil bacterial and fungal communities in response to the soil treatment with microencapsulated terpene, Brassica seedmeal (BioFence™) and chloropicrin in a field trial. The effect of soil treatments on the overall microbial population structure and relative abundance of individual microbial OTUs (operational taxonomic units) was assessed using an amplicon-based metabarcoding approach at three time points. Classifying representative OTU sequences into taxonomic groups was more uncertain for bacteria than for fungi. Chloropicrin dramatically altered both bacterial and fungal populations within four weeks of application. The effect on bacterial population structure is short-lived and became non-significant 16 weeks after treatment; however, fungal population structure was more persistently affected by chloropicrin. Neither terpene nor BioFence™ significantly affected soil microbiota. This study highlights the need for reliable algorithms in classifying sequences into taxonomic units and also the importance of identifying microbes into finer taxonomic groups for understanding soil microbiota and their effects on crop production.

The assessment of epiphytic yeast diversity in sugarcane phyllosphere in Thailand by culture-independent method

The diversity of epiphytic yeasts from sugarcane (Saccharum officinarum Linn.) phyllospheres in Thailand was investigated by culture-independent method based on the analysis of the D1/D2 domains of the large subunit rRNA gene sequences. Forty-five samples of sugarcane leaf were collected randomly from ten provinces in Thailand. A total of 1342 clones were obtained from 45 clone libraries. 426 clones (31.7 %) were closely related to yeast strains in the GenBank database, and they were clustered into 31 operational taxonomic units (OTUs) with a similarity threshold of 99 %. All OTU sequences were classified in phylum Basidiomycota which were closely related to 11 yeast species in seven genera including Cryptococcus flavus, Hannaella coprosmaensis, Rhodotorula taiwanensis, Jaminaea angkoreiensis, Malassezia restricta, Pseudozyma antarctica, Pseudozyma aphidis, Pseudozyma hubeiensis, Pseudozyma prolifica, Pseudozyma shanxiensis, and Sporobolomyces vermiculatus. The most predominant yeasts detected belonged to Ustilaginales with 89.4 % relative frequency and the prevalent yeast genus was Pseudozyma. However, the majority were unable to be identified as known yeast species and these sequences may represent the sequences of new yeast taxa. In addition, The OTU that closely related to P. prolifica was commonly detected in sugarcane phyllosphere.

Characterization of bacteria isolated from palaeoproterozoic metasediments for sequestration of carbon dioxide and formation of calcium carbonate.

Bacterial community of palaeoproterozoic metasediments was enriched in the chemostat in the presence of different concentrations of NaHCO3. Six bacterial isolates were isolated from the chemostat on nutrient agar plates on the basis of distinct morphology. Denaturing gradient gel electrophoresis (DGGE) proved the presence of six operational taxonomic units (OTUs) at 50 and 100 mM NaHCO3. The OTU was reduced to three and one at enrichment concentration of 150 and 200 mM NaHCO3 respectively. These six isolates were tested for sequestration of carbon dioxide by (14)C metabolic labeling of NaH(14)CO3. Among the six isolates, one of the bacterium showed better potency to fix radiolabeled NaH(14)CO3. The isolate (ISTD04) was identified as Serratia sp. by 16S ribosomal RNA (16S rRNA) sequence analysis and was found to be same as the DGGE OTU sequence at 200-mM NaHCO3 concentration. The bacterium was tested for product formation in form of calcite crystals in presence of 5 % CO2. Scanning electron microscopy (SEM) of product formed by the bacterium revealed defined faceted rhombohedral structure which resembled calcite and vaterite phases of the crystal. Formation of calcium carbonate crystals was further confirmed by Fourier transform infrared (FTIR) spectroscopy as carbonate group showing strong vibration at 1,456 cm(-1). Major calcite phase diffraction peaks were determined by X-ray diffraction (XRD) analysis, and energy-dispersive X-ray (EDX) analysis showed the presence of CaO (72 %) and carbon (18 %). Bacterium use bicarbonate as carbon source for their growth as well as by-product formation in form of calcite shows carbon circulation and storage.

Simultaneous profiling of seed-associated bacteria and fungi reveals antagonistic interactions between microorganisms within a shared epiphytic microbiome on Triticum and Brassica seeds.

In order to address the hypothesis that seeds from ecologically and geographically diverse plants harbor characteristic epiphytic microbiota, we characterized the bacterial and fungal microbiota associated with Triticum and Brassica seed surfaces.The total microbial complement was determined by amplification and sequencing of a fragment of chaperonin 60 (cpn60). Specific microorganisms were quantified by qPCR. Bacteria and fungi corresponding to operational taxonomic units (OTU) that were identified in the sequencing study were isolated and their interactions examined.A total of 5477 OTU were observed from seed washes. Neither total epiphytic bacterial load nor community richness/evenness was significantly different between the seed types; 578 OTU were shared among all samples at a variety of abundances. Hierarchical clustering revealed that 203 were significantly different in abundance on Triticum seeds compared with Brassica. Microorganisms isolated from seeds showed 99–100% identity between the cpn60 sequences of the isolates and the OTU sequences from this shared microbiome. Bacterial strains identified as Pantoea agglomerans had antagonistic properties toward one of the fungal isolates (Alternaria sp.), providing a possible explanation for their reciprocal abundances on both Triticum and Brassica seeds.cpn60 enabled the simultaneous profiling of bacterial and fungal microbiota and revealed a core seed-associated microbiota shared between diverse plant genera.

MPUMA: a computational approach to microbiota analysis by de novo assembly of operational taxonomic units based on protein-coding barcode sequences

BackgroundFormation of operational taxonomic units (OTU) is a common approach to data aggregation in microbial ecology studies based on amplification and sequencing of individual gene targets. The de novo assembly of OTU sequences has been recently demonstrated as an alternative to widely used clustering methods, providing robust information from experimental data alone, without any reliance on an external reference database.ResultsHere we introduce mPUMA (microbial Profiling Using Metagenomic Assembly, http://mpuma.sourceforge.net), a software package for identification and analysis of protein-coding barcode sequence data. It was developed originally for Cpn60 universal target sequences (also known as GroEL or Hsp60). Using an unattended process that is independent of external reference sequences, mPUMA forms OTUs by DNA sequence assembly and is capable of tracking OTU abundance. mPUMA processes microbial profiles both in terms of the direct DNA sequence as well as in the translated amino acid sequence for protein coding barcodes. By forming OTUs and calculating abundance through an assembly approach, mPUMA is capable of generating inputs for several popular microbiota analysis tools. Using SFF data from sequencing of a synthetic community of Cpn60 sequences derived from the human vaginal microbiome, we demonstrate that mPUMA can faithfully reconstruct all expected OTU sequences and produce compositional profiles consistent with actual community structure.ConclusionsmPUMA enables analysis of microbial communities while empowering the discovery of novel organisms through OTU assembly.

Differences in the Rumen Methanogen Populations of Lactating Jersey and Holstein Dairy Cows under the Same Diet Regimen

In the dairy cattle industry, Holstein and Jersey are the breeds most commonly used for production. They differ in performance by various traits, such as body size, milk production, and milk composition. With increased concerns about the impact of agriculture on climate change, potential differences in other traits, such as methane emission, also need to be characterized further. Since methane is produced in the rumen by methanogenic archaea, we investigated whether the population structure of methanogen communities would differ between Holsteins and Jerseys. Breed-specific rumen methanogen 16S rRNA gene clone libraries were constructed from pooled PCR products obtained from lactating Holstein and Jersey cows, generating 180 and 185 clones, respectively. The combined 365 sequences were assigned to 55 species-level operational taxonomic units (OTUs). Twenty OTUs, representing 85% of the combined library sequences, were common to both breeds, while 23 OTUs (36 sequences) were found only in the Holstein library and 12 OTUs (18 sequences) were found only in the Jersey library, highlighting increased diversity in the Holstein library. Other differences included the observation that sequences with species-like sequence identity to Methanobrevibacter millerae were represented more highly in the Jersey breed, while Methanosphaera-related sequences and novel uncultured methanogen clones were more frequent in the Holstein library. In contrast, OTU sequences with species-level sequence identity to Methanobrevibacter ruminantium were represented similarly in both libraries. Since the sampled animals were from a single herd consisting of two breeds which were fed the same diet and maintained under the same environmental conditions, the differences we observed may be due to differences in host breed genetics.

Geographic Variation in Bacterial Communities Associated With the Red Turpentine Beetle (Coleoptera: Curculionidae)

Bacterial communities are known to play important roles in insect life histories, yet their consistency or variation across populations is poorly understood. Bacteria associated with the bark beetle Dendroctonus valens LeConte from eight populations, ranging from Wisconsin to Oregon, were evaluated and compared. We used the culture-independent technique of denaturing gradient gel electrophoresis to visualize bacterial diversity, or individual operational taxonomic units (OTUs), from individual beetles. One-way analysis of similarities was used to test for differences of bacterial communities between sites. Analysis of community profiles showed that individual beetles on average contained 10 OTUs, with frequency of association from 2 to 100% of beetles. OTU sequences most closely matched beta- and gamma-proteobacteria, and one each matched Bacilli and Actinobacteria. Several OTUs were particularly abundant, most notably an Actinobacterium from 100% and two Proteobacteria from 60% of beetles sampled. Some OTUs were similar to previously described bacteria with known biochemical capabilities and ecological functions, suggesting that some bacterial associates of D. valens may contribute to its ability to exploit a resource low in nutrients and high in defensive compounds. There were significant differences of bacterial communities between sites. The strength of these differences was positively correlated with distance between sites, although additional unexplained factors also contribute to the variation.

Phylogenetic diversity of Synechococcus in the Chesapeake Bay revealed by Ribulose-1,5-bisphosphate carboxylase-oxygenase (RuBisCO) large subunit gene (rbcL) sequences

In order to understand how Synechococcus in the estuarine environment (Chesapeake Bay) are phylogenetically related to other known marine Synechococcus, partial rbcL gene sequences from 25 strains of Synechococcus spp. isolated from estuarine, coastal and oceanic waters were sequenced. The rbcL gene phylogeny showed that Chesapeake Bay Synechococcus isolates together with other marine Synechococcus spp. formed a monophyletic group which belongs to the form IA RuBisCO. All the Chesapeake Bay Synechococcus were able to grow in a wide range of salinity (0 to 30‰), and most of them belong to the Marine Cluster B (MC-B). Interestingly, several phycoerythrin (PE)-containing Synechococcus isolated from the bay were clustered in the MC-B group, which had previously only contained the non-PE Synechococcus. A set of PCR primers was developed to specifically amplify the rbcL gene from natural marine Synechococcus populations. After screening 232 clones randomly selected from 5 clone libraries (built on 5 estuarine samples respectively), at least 7 different rbcL genotypes or the operational taxonomic units (OTUs) were identified. Despite the great genetic diversity among the OTU sequences, they were all clustered with 13 Chesapeake Bay isolates. The distribution frequency of these OTUs varied dramatically from the upper to lower bay. Our results suggest that the Chesapeake estuary provides an ideal environ- ment for the MC-B type Synechococcus populations to thrive. Marine Synechococcus appear to adapt well to various ecological niches, and a clear boundary solely based on phenotypic features may not exist when more and more novel ecotypes and genotypes are unveiled with molecular tools.