Strain/species-specific probe design for microbial identification microarrays

Synthetic phylogenetically diverse communities promote denitrification and stability

Denitrification is an important process of the global nitrogen cycle as some of its intermediates are environmentally important or related to global warming. However, how the phylogenetic diversity of denitrifying communities affects their denitrification rates and temporal stability remains unclear. Here we selected denitrifiers based on their phylogenetic distance to construct two groups of synthetic denitrifying communities: one closely related (CR) group with all strains from the genus Shewanella and the other distantly related (DR) group with all constituents from different genera. All synthetic denitrifying communities (SDCs) were experimentally evolved for 200 generations. The results showed that high phylogenetic diversity followed by experimental evolution promoted the function and stability of synthetic denitrifying communities. Specifically, the productivity and denitrification rates were significantly (P < 0.05) higher with Paracocus denitrificans as the dominant species (since the 50th generation) in the DR community than those in the CR community. The DR community also showed significantly (t = 7.119, df = 10, P < 0.001) higher stability through overyielding and asynchrony of species fluctuations, and showed more complementarity than the CR group during the experimental evolution. This study has important implications for applying synthetic communities to remediate environmental problems and mitigate greenhouse gas emissions.

Synthetic Denitrifying Communities Reveal a Positive and Dynamic Biodiversity-Ecosystem Functioning Relationship during Experimental Evolution

Biodiversity is vital for ecosystem functions and services, and many studies have reported positive, negative, or neutral biodiversity-ecosystem functioning (BEF) relationships in plant and animal systems. However, if the BEF relationship exists and how it evolves remains elusive in microbial systems. Here, we selected 12 Shewanella denitrifiers to construct synthetic denitrifying communities (SDCs) with a richness gradient spanning 1 to 12 species, which were subjected to approximately 180 days (with 60 transfers) of experimental evolution with generational changes in community functions continuously tracked. A significant positive correlation was observed between community richness and functions, represented by productivity (biomass) and denitrification rate, however, such a positive correlation was transient, only significant in earlier days (0 to 60) during the evolution experiment (180 days). Also, we found that community functions generally increased throughout the evolution experiment. Furthermore, microbial community functions with lower richness exhibited greater increases than those with higher richness. Biodiversity effect analysis revealed positive BEF relationships largely attributable to complementary effects, which were more pronounced in communities with lower richness than those with higher richness. This study is one of the first studies that advances our understanding of BEF relationships and their evolutionary mechanisms in microbial systems, highlighting the crucial role of evolution in predicting the BEF relationship in microbial systems. IMPORTANCE Despite the consensus that biodiversity supports ecosystem functioning, not all experimental models of macro-organisms support this notion with positive, negative, or neutral biodiversity-ecosystem functioning (BEF) relationships reported. The fast-growing, metabolically versatile, and easy manipulation nature of microbial communities allows us to explore well the BEF relationship and further interrogate if the BEF relationship remains constant during long-term community evolution. Here, we constructed multiple synthetic denitrifying communities (SDCs) by randomly selecting species from a candidate pool of 12 Shewanella denitrifiers. These SDCs differ in species richness, spanning 1 to 12 species, and were monitored continuously for community functional shifts during approximately 180-day parallel cultivation. We demonstrated that the BEF relationship was dynamic with initially (day 0 to 60) greater productivity and denitrification among SDCs of higher richness. However, such pattern was reversed thereafter with greater productivity and denitrification increments in lower-richness SDCs, likely due to a greater accumulation of beneficial mutations during the experimental evolution.

Niche differentiation among comammox (Nitrospira inopinata) and other metabolically distinct nitrifiers

Due to global change, increasing nutrient input to ecosystems dramatically affects the nitrogen cycle, especially the nitrification process. Nitrifiers including ammonia-oxidizing archaea (AOAs), ammonia-oxidizing bacteria (AOBs), nitrite-oxidizing bacteria (NOBs), and recently discovered complete ammonia oxidizers (comammoxs) perform nitrification individually or in a community. However, much remains to be learned about their niche differentiation, coexistence, and interactions among those metabolically distinct nitrifiers. Here, we used synthetic microbial ecology approaches to construct synthetic nitrifying communities (SNCs) with different combinations of Nitrospira inopinata as comammox, Nitrososphaera gargensis as AOA, Nitrosomonas communis as AOB, and Nitrospira moscoviensis as NOB. Our results showed that niche differentiation and potential interactions among those metabolically distinct nitrifiers were determined by their kinetic characteristics. The dominant species shifted from N. inopinata to N. communis in the N4 community (with all four types of nitrifiers) as ammonium concentrations increased, which could be well explained by the kinetic difference in ammonia affinity, specific growth rate, and substrate tolerance of nitrifiers in the SNCs. In addition, a conceptual model was developed to infer niche differentiation and possible interactions among the four types of nitrifiers. This study advances our understanding of niche differentiation and provides new strategies to further study their interactions among the four types of nitrifiers.

Accurate identification of taxon-specific molecular markers in plants based on DNA signature sequence

Accurate identification of plants remains a significant challenge for taxonomists and is the basis for plant diversity conservation. Although DNA barcoding methods are commonly used for plant identification, these are limited by the low amplification success and low discriminative power of selected genomic regions. In this study, we developed a k-mer-based approach, the DNA signature sequence (DSS), to accurately identify plant taxon-specific markers, especially at the species level. DSS is a constant-length nucleotide sequence capable of identifying a taxon and distinguishing it from other taxa. In this study, we performed the first large-scale study of DSS markers in plants. DSS candidates of 3,899 angiosperm plant species were calculated based on a chloroplast dataset with 4,356 assemblies. Using Sanger sequencing of PCR amplicons and high-throughput sequencing, DSSs were validated in four and 165 species, respectively. Based on this, the universality of the DSSs was over 79.38%. Several indicators influencing DSS marker identification and detection have also been evaluated, and common criteria for DSS application in plant identification have been proposed.

Bacillus subtilis Produces Amino Acids to Stimulate Protein Synthesis in Ruminal Tissue Explants via the Phosphatidylinositol-4,5-Bisphosphate 3-Kinase Catalytic Subunit Beta–Serine/Threonine Kinase–Mammalian Target of Rapamycin Complex 1 Pathway

Background

Bacillus subtilis is a probiotic strain that is widely used as a feed supplement for ruminants. In this study, one B. subtilis strain isolated from the ruminal fluid of Holstein dairy cows was used for an ex vivo study with ruminal tissue explants. The main goal was to assess the potential endosymbiotic links between B. subtilis and the ruminal epithelium using molecular analyses and amino acid profiling. The explant culture protocol was first optimized to determine the ideal conditions in terms of tissue viability before performing the actual experiments involving active and inactive bacteria with or without protein synthesis inhibitors, such as LY294002 (phosphatidylinositol 3-kinase inhibitor) or rapamycin [mammalian target of rapamycin (mTOR) inhibitor].

Results

The mRNA levels of phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit beta (PIK3CB), serine/threonine kinase (AKT), mTOR, P70S6K1, and eukaryotic translation initiation factor 4E binding protein 1 were the highest (p < 0.01), while those of programmed cell death 4 were the lowest when the tissue was incubated with 10⁷ of B. subtilis. Compared with the inactivated bacteria, the expression levels of PIK3CB and AKT, and overall changes in mTOR and P70S6K1 were greater in rumen explants with living bacteria (p < 0.05). With an increase in B. subtilis concentration, the trends of protein and corresponding gene changes were consistent. There were differences in the concentrations of individual amino acids in the supernatants of living and inactivated bacterial culture groups, with most amino acids enriched in pathways, such as aminoacyl tRNA biosynthesis, cyanoamino acid metabolism, monobactam biosynthesis, or glycine, serine, and threonine metabolism. The addition of psilocybin upregulated the expression levels of PIK3CB and AKT. A significant decrease (p < 0.05) in PIK3CB and mTOR protein expression levels was detected after the addition of LY294002 and rapamycin. In addition, These responses were associated with the downregulation (p < 0.05) of AKT and P70S6K protein expression levels.

Conclusions

We confirmed that the in vivo ruminal tissue culture system is a suitable model for studying probiotic-induced alterations in tissue function. As such, this study provides a means for future mechanistic studies related to microbial regulation and the dietary supply of proteins. In addition, living and inactivated B. subtilis can promote protein synthesis in ruminal tissue explants by altering the expression levels of related factors in the PIK3CB–AKT–mTORC1 pathway, which could further aid in optimizing the feed efficiency and increasing the use of inactivated bacteria as additives in dairy cow farming.

Macroecological distributions of gene variants highlight the functional organization of soil microbial systems

The recent application of macroecological tools and concepts has made it possible to identify consistent patterns in the distribution of microbial biodiversity, which greatly improved our understanding of the microbial world at large scales. However, the distribution of microbial functions remains largely uncharted from the macroecological point of view. Here, we used macroecological models to examine how the genes encoding the functional capabilities of microorganisms are distributed within and across soil systems. Models built using functional gene array data from 818 soil microbial communities showed that the occupancy-frequency distributions of genes were bimodal in every studied site, and that their rank-abundance distributions were best described by a lognormal model. In addition, the relationships between gene occupancy and abundance were positive in all sites. This allowed us to identify genes with high abundance and ubiquitous distribution (core) and genes with low abundance and limited spatial distribution (satellites), and to show that they encode different sets of microbial traits. Common genes encode microbial traits related to the main biogeochemical cycles (C, N, P and S) while rare genes encode traits related to adaptation to environmental stresses, such as nutrient limitation, resistance to heavy metals and degradation of xenobiotics. Overall, this study characterized for the first time the distribution of microbial functional genes within soil systems, and highlight the interest of macroecological models for understanding the functional organization of microbial systems across spatial scales.

Species-specific genomic sequences for classification of bacteria

Modern bacterial classification relies on genomic relatedness. Genetic variation in bacterial populations present a big challenge for taxonomic classification and recently several bacterial species have been reclassified based on the intra-species genome comparison. These were facilitated by next generation sequencing technologies and advances in genome comparison approaches which led to the rearrangement of diverse bacterial species and revolution in the microbial classification system. One of the outcome of these studies is the development of suitable DNA barcodes as reliable and cost-effective method for identifying various bacterial genera. Towards refining this further, we have applied a genome comparison approach in 1104 bacterial genome assemblies (excluding plasmids) to identify unique genomic segments among intra-species genome assemblies. Using extensive bioinformatics analysis, we have identified species-specific genomic regions and designed unique primers for 100 different species (belonging to 62 genera) which includes 62 pathogenic and 13 opportunistic pathogenic bacterial species and built a database (http://slsdb.manipal.edu/bact/). These species-specific genomic regions will have a major impact on in silico and molecular methods aimed at bacterial classification and identification. These may also serve as better DNA barcodes than the markers currently used for delineation of bacteria and may also find application in various translational research programs.

Unique k-mers as Strain-Specific Barcodes for Phylogenetic Analysis and Natural Microbiome Profiling

The need for a comparative analysis of natural metagenomes stimulated the development of new methods for their taxonomic profiling. Alignment-free approaches based on the search for marker k-mers turned out to be capable of identifying not only species, but also strains of microorganisms with known genomes. Here, we evaluated the ability of genus-specific k-mers to distinguish eight phylogroups of Escherichia coli (A, B1, C, E, D, F, G, B2) and assessed the presence of their unique 22-mers in clinical samples from microbiomes of four healthy people and four patients with Crohn's disease. We found that a phylogenetic tree inferred from the pairwise distance matrix for unique 18-mers and 22-mers of 124 genomes was fully consistent with the topology of the tree, obtained with concatenated aligned sequences of orthologous genes. Therefore, we propose strain-specific "barcodes" for rapid phylotyping. Using unique 22-mers for taxonomic analysis, we detected microbes of all groups in human microbiomes; however, their presence in the five samples was significantly different. Pointing to the intraspecies heterogeneity of E. coli in the natural microflora, this also indicates the feasibility of further studies of the role of this heterogeneity in maintaining population homeostasis.

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.

HuMiChip2 for strain level identification and functional profiling of human microbiomes

With the massive data generated by the Human Microbiome Project, how to transform such data into useful information and knowledge remains challenging. Here, with currently available sequencing information (reference genomes and metagenomes), we have developed a comprehensive microarray, HuMiChip2, for strain-level identification and functional characterization of human microbiomes. HuMiChip2 was composed of 29,467 strain-specific probes targeting 2063 microbial strains/species and 133,924 sequence- and group-specific probes targeting 157 key functional gene families involved in various metabolic pathways and host-microbiome interaction processes. Computational evaluation of strain-specific probes suggested that they were not only specific to mock communities of sequenced microorganisms and metagenomes from different human body sites but also to non-sequenced microbial strains. Experimental evaluation of strain-specific probes using single strains/species and mock communities suggested a high specificity of these probes with their corresponding targets. Application of HuMiChip2 to human gut microbiome samples showed the patient microbiomes of alcoholic liver cirrhosis significantly (p < 0.05) shifted their functional structure from the healthy individuals, and the relative abundance of 21 gene families significantly (p < 0.1) differed between the liver cirrhosis patients and healthy individuals. At the strain level, five Bacteroides strains were significantly (p < 0.1) and more frequently detected in liver cirrhosis patients. These results suggest that the developed HuMiChip2 is a useful microbial ecological microarray for both strain-level identification and functional profiling of human microbiomes.

Probe Design Strategies for Oligonucleotide Microarrays

Oligonucleotide microarrays have been widely used for gene detection and/or quantification of gene expression in various samples ranging from a single organism to a complex microbial assemblage. The success of a microarray experiment, however, strongly relies on the quality of designed probes. Consequently, probe design is of critical importance and therefore multiple parameters should be considered for each probe in order to ensure high specificity, sensitivity, and uniformity as well as potentially quantitative power. Moreover, to assess the complete gene repertoire of complex biological samples such as those studied in the field of microbial ecology, exploratory probe design strategies must be also implemented to target not-yet-described sequences. To design such probes, two algorithms, KASpOD and HiSpOD, have been developed and they are available via two user-friendly web services. Here, we describe the use of this software necessary for the design of highly effective probes especially in the context of microbial oligonucleotide microarrays by taking into account all the crucial parameters.

Strategies to develop strain-specific PCR based assays for probiotics

Since health benefits conferred by probiotics are strain-specific, identification to the strain level is mandatory to allow the monitoring of the presence and the abundance of specific probiotic in a product or in a gastrointestinal tract. Compared to standard plate counts, the reduced duration of the assays and higher specificity makes PCR-based methods (standard PCR and quantitative PCR) very appropriate for detection or quantification of probiotics. Development of strain-specific assay consists of 4 main stages: (1) strain-specific marker identification; (2) construction of potential strain-specific primers; (3) validation on DNA from pure cultures of target and related strains; and (4) validation on spiked samples. The most important and also the most challenging step is the identification of strain-specific sequences, which can be subsequently targeted by specific primers or probes. Such regions can be identified on sequences derived from 16S-23S internally transcribed spacers, randomly amplified polymorphic DNA, representational difference analysis and suppression subtractive hybridisation. Already known phenotypic or genotypic characteristics of the target strain can also be used to develop the strain-specific assay. However, the initial stage of strain-specific assay development can be replaced by comparative genomics analysis of target genome with related genomes in public databases. Advances in whole genome sequencing (WGS) have resulted in a cost reduction for bacterial genome sequencing and consequently have made this approach available to most laboratories. In the present paper I reviewed the available literature on PCR and qPCR assays developed for detection of a specific probiotic strain and discussed future WGS and comparative genomics-based approaches.

Nucleic Acid-based approaches for detection of viral hepatitis

CONTEXT

To determining suitable nucleic acid diagnostics for individual viral hepatitis agent, an extensive search using related keywords was done in major medical library and data were collected, categorized, and summarized in different sections.

RESULTS

Various types of molecular biology tools can be used to detect and quantify viral genomic elements and analyze the sequences. These molecular assays are proper technologies for rapidly detecting viral agents with high accuracy, high sensitivity, and high specificity. Nonetheless, the application of each diagnostic method is completely dependent on viral agent.

CONCLUSIONS

Despite rapidity, automation, accuracy, cost-effectiveness, high sensitivity, and high specificity of molecular techniques, each type of molecular technology has its own advantages and disadvantages.

Development of HuMiChip for functional profiling of human microbiomes

Understanding the diversity, composition, structure, function, and dynamics of human microbiomes in individual human hosts is crucial to reveal human-microbial interactions, especially for patients with microbially mediated disorders, but challenging due to the high diversity of the human microbiome. Here we have developed a functional gene-based microarray for profiling human microbiomes (HuMiChip) with 36,802 probes targeting 50,007 protein coding sequences for 139 key functional gene families. Computational evaluation suggested all probes included are highly specific to their target sequences. HuMiChip was used to analyze human oral and gut microbiomes, showing significantly different functional gene profiles between oral and gut microbiome. Obvious shifts of microbial functional structure and composition were observed for both patients with dental caries and periodontitis from moderate to advanced stages, suggesting a progressive change of microbial communities in response to the diseases. Consistent gene family profiles were observed by both HuMiChip and next generation sequencing technologies. Additionally, HuMiChip was able to detect gene families at as low as 0.001% relative abundance. The results indicate that the developed HuMiChip is a useful and effective tool for functional profiling of human microbiomes.