ABSTRACT
The pathophysiology of Gulf War Illness (GWI) remains elusive even after three decades. The persistence of multiple complex symptoms along with metabolic disorders such as obesity worsens the health of present Gulf War (GW) Veterans often by the interactions of the host gut microbiome and inflammatory mediators. In this study, we hypothesized that the administration of a Western diet might alter the host metabolomic profile, which is likely associated with the altered bacterial species. Using a five-month symptom persistence GWI model in mice and whole-genome sequencing, we characterized the species-level dysbiosis and global metabolomics, along with heterogenous co-occurrence network analysis, to study the bacteriome-metabolomic association. Microbial analysis at the species level showed a significant alteration of beneficial bacterial species. The beta diversity of the global metabolomic profile showed distinct clustering due to the Western diet, along with the alteration of metabolites associated with lipid, amino acid, nucleotide, vitamin, and xenobiotic metabolism pathways. Network analysis showed novel associations of gut bacterial species with metabolites and biochemical pathways that could be used as biomarkers or therapeutic targets to ameliorate symptom persistence in GW Veterans.
Subject(s)
Dysbiosis , Gastrointestinal Microbiome , Mice , Animals , Gulf War , Diet, Western , Gastrointestinal Microbiome/physiology , Bacteria , ObesityABSTRACT
This study aims to investigate the presence of SARS-CoV-2 in public spaces and assess the utility of inexpensive air purifiers equipped with high-efficiency particulate air (HEPA) filters for viral detection. Samples were collected from six community-based organizations in underserved minority neighborhoods in Northwest Miami, Florida, from February to May 2022. Reverse transcription-quantitative polymerase chain reaction (RT-qPCR) was used to detect SARS-CoV-2 in air purifier filters and surface swabs. Among 32 filters tested, three yielded positive results, while no positive surface swabs were found. Notably, positive samples were obtained exclusively from child daycare centers. These findings highlight the potential for airborne transmission of SARS-CoV-2 in indoor air, particularly in child daycare centers. Moreover, the study demonstrates the effectiveness of readily available HEPA filters in detecting the virus. Improving indoor ventilation and implementing air filtration systems are crucial in reducing COVID-19 transmission where people gather. Air filtration systems incorporating HEPA filters offer a valuable approach to virus detection and reducing transmission risks. Future research should explore the applicability of this technology for early identification and mitigation of viral outbreaks.
Subject(s)
Air Filters , Air Pollution, Indoor , COVID-19 , Child , Humans , SARS-CoV-2 , Air Pollution, Indoor/analysis , Environmental Monitoring , Dust/analysisABSTRACT
Following the initial surges of the Alpha (B.1.1.7) and the Beta (B.1.351) variants, a more infectious Delta variant (B.1.617.2) is now surging, further deepening the health crises caused by the pandemic. The sharp rise in cases attributed to the Delta variant has made it especially disturbing and is a variant of concern. Fortunately, current vaccines offer protection against known variants of concern, including the Delta variant. However, the Delta variant has exhibited some ability to dodge the immune system as it is found that neutralizing antibodies from prior infections or vaccines are less receptive to binding with the Delta spike protein. Here, we investigated the structural changes caused by the mutations in the Delta variant's receptor-binding interface and explored the effects on binding with the ACE2 receptor as well as with neutralizing antibodies. We find that the receptor-binding ß-loop-ß motif adopts an altered but stable conformation causing separation in some of the antibody binding epitopes. Our study shows reduced binding of neutralizing antibodies and provides a possible mechanism for the immune evasion exhibited by the Delta variant.
Subject(s)
Angiotensin-Converting Enzyme 2/immunology , COVID-19/immunology , Immune Evasion/immunology , Mutation/immunology , SARS-CoV-2/immunology , Spike Glycoprotein, Coronavirus/immunology , Amino Acids/genetics , Amino Acids/immunology , Amino Acids/metabolism , Angiotensin-Converting Enzyme 2/genetics , Angiotensin-Converting Enzyme 2/metabolism , Antibodies, Viral/immunology , Binding Sites/genetics , Binding Sites/immunology , COVID-19/metabolism , COVID-19/virology , Humans , Immune Evasion/genetics , Molecular Dynamics Simulation , Mutation/genetics , Neutralization Tests , Protein Binding , Protein Domains , SARS-CoV-2/genetics , SARS-CoV-2/metabolism , Spike Glycoprotein, Coronavirus/chemistry , Spike Glycoprotein, Coronavirus/geneticsABSTRACT
BACKGROUND: Long non-coding RNA plays a vital role in changing the expression profiles of various target genes that lead to cancer development. Thus, identifying prognostic lncRNAs related to different cancers might help in developing cancer therapy. METHOD: To discover the critical lncRNAs that can identify the origin of different cancers, we propose the use of the state-of-the-art deep learning algorithm concrete autoencoder (CAE) in an unsupervised setting, which efficiently identifies a subset of the most informative features. However, CAE does not identify reproducible features in different runs due to its stochastic nature. We thus propose a multi-run CAE (mrCAE) to identify a stable set of features to address this issue. The assumption is that a feature appearing in multiple runs carries more meaningful information about the data under consideration. The genome-wide lncRNA expression profiles of 12 different types of cancers, with a total of 4768 samples available in The Cancer Genome Atlas (TCGA), were analyzed to discover the key lncRNAs. The lncRNAs identified by multiple runs of CAE were added to a final list of key lncRNAs that are capable of identifying 12 different cancers. RESULTS: Our results showed that mrCAE performs better in feature selection than single-run CAE, standard autoencoder (AE), and other state-of-the-art feature selection techniques. This study revealed a set of top-ranking 128 lncRNAs that could identify the origin of 12 different cancers with an accuracy of 95%. Survival analysis showed that 76 of 128 lncRNAs have the prognostic capability to differentiate high- and low-risk groups of patients with different cancers. CONCLUSION: The proposed mrCAE, which selects actual features, outperformed the AE even though it selects the latent or pseudo-features. By selecting actual features instead of pseudo-features, mrCAE can be valuable for precision medicine. The identified prognostic lncRNAs can be further studied to develop therapies for different cancers.
Subject(s)
Algorithms , Biomarkers, Tumor/genetics , Deep Learning , Gene Expression Regulation, Neoplastic , Neoplasms/pathology , Neural Networks, Computer , RNA, Long Noncoding/genetics , Humans , Neoplasms/classification , Neoplasms/genetics , Precision Medicine , Prognosis , Survival RateABSTRACT
BACKGROUND: Partial Least-Squares Discriminant Analysis (PLS-DA) is a popular machine learning tool that is gaining increasing attention as a useful feature selector and classifier. In an effort to understand its strengths and weaknesses, we performed a series of experiments with synthetic data and compared its performance to its close relative from which it was initially invented, namely Principal Component Analysis (PCA). RESULTS: We demonstrate that even though PCA ignores the information regarding the class labels of the samples, this unsupervised tool can be remarkably effective as a feature selector. In some cases, it outperforms PLS-DA, which is made aware of the class labels in its input. Our experiments range from looking at the signal-to-noise ratio in the feature selection task, to considering many practical distributions and models encountered when analyzing bioinformatics and clinical data. Other methods were also evaluated. Finally, we analyzed an interesting data set from 396 vaginal microbiome samples where the ground truth for the feature selection was available. All the 3D figures shown in this paper as well as the supplementary ones can be viewed interactively at http://biorg.cs.fiu.edu/plsda CONCLUSIONS: Our results highlighted the strengths and weaknesses of PLS-DA in comparison with PCA for different underlying data models.
Subject(s)
Computational Biology , Discriminant Analysis , Least-Squares Analysis , Machine Learning , Principal Component AnalysisABSTRACT
BACKGROUND: Microbe-microbe and host-microbe interactions in a microbiome play a vital role in both health and disease. However, the structure of the microbial community and the colonization patterns are highly complex to infer even under controlled wet laboratory conditions. In this study, we investigate what information, if any, can be provided by a Bayesian Network (BN) about a microbial community. Unlike the previously proposed Co-occurrence Networks (CoNs), BNs are based on conditional dependencies and can help in revealing complex associations. RESULTS: In this paper, we propose a way of combining a BN and a CoN to construct a signed Bayesian Network (sBN). We report a surprising association between directed edges in signed BNs and known colonization orders. CONCLUSIONS: BNs are powerful tools for community analysis and extracting influences and colonization patterns, even though the analysis only uses an abundance matrix with no temporal information. We conclude that directed edges in sBNs when combined with negative correlations are consistent with and strongly suggestive of colonization order.
Subject(s)
Microbiota , Bayes TheoremABSTRACT
MOTIVATION: Bacterial metagenomics profiling for metagenomic whole sequencing (mWGS) usually starts by aligning sequencing reads to a collection of reference genomes. Current profiling tools are designed to work against a small representative collection of genomes, and do not scale very well to larger reference genome collections. However, large reference genome collections are capable of providing a more complete and accurate profile of the bacterial population in a metagenomics dataset. In this paper, we discuss a scalable, efficient and affordable approach to this problem, bringing big data solutions within the reach of laboratories with modest resources. RESULTS: We developed Flint, a metagenomics profiling pipeline that is built on top of the Apache Spark framework, and is designed for fast real-time profiling of metagenomic samples against a large collection of reference genomes. Flint takes advantage of Spark's built-in parallelism and streaming engine architecture to quickly map reads against a large (170 GB) reference collection of 43 552 bacterial genomes from Ensembl. Flint runs on Amazon's Elastic MapReduce service, and is able to profile 1 million Illumina paired-end reads against over 40 K genomes on 64 machines in 67 s-an order of magnitude faster than the state of the art, while using a much larger reference collection. Streaming the sequencing reads allows this approach to sustain mapping rates of 55 million reads per hour, at an hourly cluster cost of $8.00 USD, while avoiding the necessity of storing large quantities of intermediate alignments. AVAILABILITY AND IMPLEMENTATION: Flint is open source software, available under the MIT License (MIT). Source code is available at https://github.com/camilo-v/flint. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.
Subject(s)
Cloud Computing , High-Throughput Nucleotide Sequencing , Microbiota , Algorithms , Metagenomics , Sequence Analysis, DNA , SoftwareABSTRACT
BACKGROUND: Computing centrality is a foundational concept in social networking that involves finding the most "central" or important nodes. In some biological networks defining importance is difficult, which then creates challenges in finding an appropriate centrality algorithm. RESULTS: We instead generalize the results of any k centrality algorithms through our iterative algorithm MATRIA, producing a single ranked and unified set of central nodes. Through tests on three biological networks, we demonstrate evident and balanced correlations with the results of these k algorithms. We also improve its speed through GPU parallelism. CONCLUSIONS: Our results show iteration to be a powerful technique that can eliminate spatial bias among central nodes, increasing the level of agreement between algorithms with various importance definitions. GPU parallelism improves speed and makes iteration a tractable problem for larger networks.
Subject(s)
Algorithms , Animals , Bacteria/genetics , Computer Graphics , Gene Regulatory Networks , Ostreidae/genetics , Time FactorsABSTRACT
Motivation: Software pipelines have become almost standardized tools for microbiome analysis. Currently many pipelines are available, often sharing some of the same algorithms as stages. This is largely because each pipeline has its own source language and file formats, making it typically more economical to reinvent the wheel than to learn and interface to an existing package. We present Plugin-Based Microbiome Analysis (PluMA), which addresses this problem by providing a lightweight back end that can be infinitely extended using dynamically loaded plugin extensions. These can be written in one of many compiled or scripting languages. With PluMA and its online plugin pool, algorithm designers can easily plug-and-play existing pipeline stages with no knowledge of their underlying implementation, allowing them to efficiently test a new algorithm alongside these stages or combine them in a new and creative way. Results: We demonstrate the usefulness of PluMA through an example pipeline (P-M16S) that expands an obesity study involving gut microbiome samples from the mouse, by integrating multiple plugins using a variety of source languages and file formats, and producing new results. Availability and implementation: Links to github repositories for the PluMA source code and P-M16S, in addition to the plugin pool are available from the Bioinformatics Research Group (BioRG) at: http://biorg.cis.fiu.edu/pluma.
Subject(s)
Microbiota , Algorithms , Animals , Gastrointestinal Microbiome , Mice , SoftwareABSTRACT
BACKGROUND: The notion of centrality is used to identify "important" nodes in social networks. Importance of nodes is not well-defined, and many different notions exist in the literature. The challenge of defining centrality in meaningful ways when network edges can be positively or negatively weighted has not been adequately addressed in the literature. Existing centrality algorithms also have a second shortcoming, i.e., the list of the most central nodes are often clustered in a specific region of the network and are not well represented across the network. METHODS: We address both by proposing Ablatio Triadum (ATria), an iterative centrality algorithm that uses the concept of "payoffs" from economic theory. RESULTS: We compare our algorithm with other known centrality algorithms and demonstrate how ATria overcomes several of their shortcomings. We demonstrate the applicability of our algorithm to synthetic networks as well as biological networks including bacterial co-occurrence networks, sometimes referred to as microbial social networks. CONCLUSIONS: We show evidence that ATria identifies three different kinds of "important" nodes in microbial social networks with different potential roles in the community.
Subject(s)
Algorithms , Computational Biology/methods , Models, Biological , Bacterial Physiological Phenomena , SoftwareABSTRACT
Using microarray and bioinformatics, we examined the gene expression profiles in transgenic mouse hearts expressing mutations in the myosin regulatory light chain shown to cause hypertrophic cardiomyopathy (HCM). We focused on two malignant RLC-mutations, Arginine 58âGlutamine (R58Q) and Aspartic Acid 166 â Valine (D166V), and one benign, Lysine 104 â Glutamic Acid (K104E)-mutation. Datasets of differentially expressed genes for each of three mutants were compared to those observed in wild-type (WT) hearts. The changes in the mutant vs. WT samples were shown as fold-change (FC), with stringency FC ≥ 2. Based on the gene profiles, we have identified the major signaling pathways that underlie the R58Q-, D166V- and K104E-HCM phenotypes. The correlations between different genotypes were also studied using network-based algorithms. Genes with strong correlations were clustered into one group and the central gene networks were identified for each HCM mutant. The overall gene expression patterns in all mutants were distinct from the WT profiles. Both malignant mutations shared certain classes of genes that were up or downregulated, but most similarities were noted between D166V and K104E mice, with R58Q hearts showing a distinct gene expression pattern. Our data suggest that all three HCM mice lead to cardiomyopathy in a mutation-specific manner and thus develop HCM through diverse mechanisms.
Subject(s)
Cardiomyopathy, Hypertrophic/genetics , Cardiomyopathy, Hypertrophic/metabolism , Gene Expression Regulation , Mutation , Myosin Light Chains/metabolism , Algorithms , Animals , Arginine/chemistry , Computational Biology , Gene Expression Profiling , Glutamic Acid/chemistry , Glutamine/chemistry , Lysine/chemistry , Mice , Mice, Transgenic , Multigene Family , Myocardium/metabolism , Myosin Light Chains/genetics , Oligonucleotide Array Sequence Analysis , Phenotype , Principal Component Analysis , Valine/chemistryABSTRACT
Pathogenicity of Pseudomonas aeruginosa, a major cause of many acute and chronic human infections, is determined by tightly regulated expression of multiple virulence factors. Quorum sensing (QS) controls expression of many of these pathogenic determinants. Previous microarray studies have shown that the AmpC ß-lactamase regulator AmpR, a member of the LysR family of transcription factors, also controls non-ß-lactam resistance and multiple virulence mechanisms. Using RNA-Seq and complementary assays, this study further expands the AmpR regulon to include diverse processes such as oxidative stress, heat shock and iron uptake. Importantly, AmpR affects many of these phenotypes, in part, by regulating expression of non-coding RNAs such as rgP32, asRgsA, asPrrF1 and rgRsmZ. AmpR positively regulates expression of the major QS regulators LasR, RhlR and MvfR, and genes of the Pseudomonas quinolone system. Chromatin immunoprecipitation (ChIP)-Seq and ChIP-quantitative real-time polymerase chain reaction studies show that AmpR binds to the ampC promoter both in the absence and presence of ß-lactams. In addition, AmpR directly binds the lasR promoter, encoding the QS master regulator. Comparison of the AmpR-binding sequences from the transcriptome and ChIP-Seq analyses identified an AT-rich consensus-binding motif. This study further attests to the role of AmpR in regulating virulence and physiological processes in P. aeruginosa.
Subject(s)
Bacterial Proteins/metabolism , Gene Expression Regulation, Bacterial , Pseudomonas aeruginosa/genetics , RNA, Small Untranslated/metabolism , Regulon , Transcription Factors/metabolism , Bacterial Proteins/genetics , Gene Expression Profiling , Heat-Shock Response/genetics , High-Throughput Nucleotide Sequencing , Iron/metabolism , Oligonucleotide Array Sequence Analysis , Operon , Oxidative Stress/genetics , Phenazines/metabolism , Pseudomonas aeruginosa/metabolism , Pseudomonas aeruginosa/pathogenicity , Quorum Sensing , Sequence Analysis, RNA , Trans-Activators/geneticsABSTRACT
BACKGROUND: It is well understood that distinct communities of bacteria are present at different sites of the body, and that changes in the structure of these communities have strong implications for human health. Yet, challenges remain in understanding the complex interconnections between the bacterial taxa within these microbial communities and how they change during the progression of diseases. Many recent studies attempt to analyze the human microbiome using traditional ecological measures and cataloging differences in bacterial community membership. In this paper, we show how to push metagenomic analyses beyond mundane questions related to the bacterial taxonomic profiles that differentiate one sample from another. METHODS: We develop tools and techniques that help us to investigate the nature of social interactions in microbial communities, and demonstrate ways of compactly capturing extensive information about these networks and visually conveying them in an effective manner. We define the concept of bacterial "social clubs", which are groups of taxa that tend to appear together in many samples. More importantly, we define the concept of "rival clubs", entire groups that tend to avoid occurring together in many samples. We show how to efficiently compute social clubs and rival clubs and demonstrate their utility with the help of examples including a smokers' dataset and a dataset from the Human Microbiome Project (HMP). RESULTS: The tools developed provide a framework for analyzing relationships between bacterial taxa modeled as bacterial co-occurrence networks. The computational techniques also provide a framework for identifying clubs and rival clubs and for studying differences in the microbiomes (and their interactions) of two or more collections of samples. CONCLUSIONS: Microbial relationships are similar to those found in social networks. In this work, we assume that strong (positive or negative) tendencies to co-occur or co-infect is likely to have biological, physiological, or ecological significance, possibly as a result of cooperation or competition. As a consequence of the analysis, a variety of biological interpretations are conjectured. In the human microbiome context, the pattern of strength of interactions between bacterial taxa is unique to body site.
Subject(s)
Bacteria/genetics , Bacterial Physiological Phenomena , Metagenomics/methods , Bacteria/classification , Female , Humans , Male , Microbiota , Middle Aged , SmokingABSTRACT
Pseudomonas aeruginosa is a dreaded pathogen in many clinical settings. Its inherent and acquired antibiotic resistance thwarts therapy. In particular, derepression of the AmpC ß-lactamase is a common mechanism of ß-lactam resistance among clinical isolates. The inducible expression of ampC is controlled by the global LysR-type transcriptional regulator (LTTR) AmpR. In the present study, we investigated the genetic and structural elements that are important for ampC induction. Specifically, the ampC (PampC) and ampR (PampR) promoters and the AmpR protein were characterized. The transcription start sites (TSSs) of the divergent transcripts were mapped using 5' rapid amplification of cDNA ends-PCR (RACE-PCR), and strong σ(54) and σ(70) consensus sequences were identified at PampR and PampC, respectively. Sigma factor RpoN was found to negatively regulate ampR expression, possibly through promoter blocking. Deletion mapping revealed that the minimal PampC extends 98 bp upstream of the TSS. Gel shifts using membrane fractions showed that AmpR binds to PampC in vitro whereas in vivo binding was demonstrated using chromatin immunoprecipitation-quantitative PCR (ChIP-qPCR). Additionally, site-directed mutagenesis of the AmpR helix-turn-helix (HTH) motif identified residues critical for binding and function (Ser38 and Lys42) and critical for function but not binding (His39). Amino acids Gly102 and Asp135, previously implicated in the repression state of AmpR in the enterobacteria, were also shown to play a structural role in P. aeruginosa AmpR. Alkaline phosphatase fusion and shaving experiments suggest that AmpR is likely to be membrane associated. Lastly, an in vivo cross-linking study shows that AmpR dimerizes. In conclusion, a potential membrane-associated AmpR dimer regulates ampC expression by direct binding.
Subject(s)
Bacterial Proteins/metabolism , Gene Expression Regulation, Bacterial/physiology , Pseudomonas aeruginosa/metabolism , Amino Acid Motifs , Anti-Bacterial Agents/pharmacology , Bacterial Proteins/chemistry , Bacterial Proteins/genetics , Chromosome Mapping , Chromosomes, Bacterial , Consensus Sequence , Drug Resistance, Bacterial , Promoter Regions, Genetic , Protein Binding , Protein Conformation , Pseudomonas aeruginosa/drug effects , Pseudomonas aeruginosa/genetics , beta-Lactam Resistance , beta-Lactamases/genetics , beta-Lactamases/metabolism , beta-Lactams/pharmacologyABSTRACT
A key challenge in the analysis of microbiome data is the integration of multi-omic datasets and the discovery of interactions between microbial taxa, their expressed genes, and the metabolites they consume and/or produce. In an effort to improve the state of the art in inferring biologically meaningful multi-omic interactions, we sought to address some of the most fundamental issues in causal inference from longitudinal multi-omics microbiome data sets. We developed METALICA, a suite of tools and techniques that can infer interactions between microbiome entities. METALICA introduces novel unrolling and de-confounding techniques used to uncover multi-omic entities that are believed to act as confounders for some of the relationships that may be inferred using standard causal inferencing tools. The results lend support to predictions about biological models and processes by which microbial taxa interact with each other in a microbiome. The unrolling process helps identify putative intermediaries (genes and/or metabolites) to explain the interactions between microbes; the de-confounding process identifies putative common causes that may lead to spurious relationships to be inferred. METALICA was applied to the networks inferred by existing causal discovery, and network inference algorithms were applied to a multi-omics data set resulting from a longitudinal study of IBD microbiomes. The most significant unrollings and de-confoundings were manually validated using the existing literature and databases. IMPORTANCE: We have developed a suite of tools and techniques capable of inferring interactions between microbiome entities. METALICA introduces novel techniques called unrolling and de-confounding that are employed to uncover multi-omic entities considered to be confounders for some of the relationships that may be inferred using standard causal inferencing tools. To evaluate our method, we conducted tests on the inflammatory bowel disease (IBD) dataset from the iHMP longitudinal study, which we pre-processed in accordance with our previous work. From this dataset, we generated various subsets, encompassing different combinations of metagenomics, metabolomics, and metatranscriptomics datasets. Using these multi-omics datasets, we demonstrate how the unrolling process aids in the identification of putative intermediaries (genes and/or metabolites) to explain the interactions between microbes. Additionally, the de-confounding process identifies potential common causes that may give rise to spurious relationships to be inferred. The most significant unrollings and de-confoundings were manually validated using the existing literature and databases.
ABSTRACT
Unlike chronological age, biological age is a strong indicator of health of an individual. However, the molecular fingerprint associated with biological age is ill-defined. To define a high-resolution signature of biological age, we analyzed metabolome, circulating senescence-associated secretome (SASP)/inflammation markers and the interaction between them, from a cohort of healthy and rapid agers. The balance between two fatty acid oxidation mechanisms, ß-oxidation and ω-oxidation, associated with the extent of functional aging. Furthermore, a panel of 25 metabolites, Healthy Aging Metabolic (HAM) index, predicted healthy agers regardless of gender and race. HAM index was also validated in an independent cohort. Causal inference with machine learning implied three metabolites, ß-cryptoxanthin, prolylhydroxyproline, and eicosenoylcarnitine as putative drivers of biological aging. Multiple SASP markers were also elevated in rapid agers. Together, our findings reveal that a network of metabolic pathways underlie biological aging, and the HAM index could serve as a predictor of phenotypic aging in humans.
Subject(s)
Cellular Senescence , Secretome , Humans , Aging/genetics , Aging/metabolism , Metabolome , Biomarkers/metabolismABSTRACT
A key challenge in the analysis of microbiome data is the integration of multi-omic datasets and the discovery of interactions between microbial taxa, their expressed genes, and the metabolites they consume and/or produce. In an effort to improve the state-of-the-art in inferring biologically meaningful multi-omic interactions, we sought to address some of the most fundamental issues in causal inference from longitudinal multi-omics microbiome data sets. We developed METALICA, a suite of tools and techniques that can infer interactions between microbiome entities. METALICA introduces novel unrolling and de-confounding techniques used to uncover multi-omic entities that are believed to act as confounders for some of the relationships that may be inferred using standard causal inferencing tools. The results lend support to predictions about biological models and processes by which microbial taxa interact with each other in a microbiome. The unrolling process helps to identify putative intermediaries (genes and/or metabolites) to explain the interactions between microbes; the de-confounding process identifies putative common causes that may lead to spurious relationships to be inferred. METALICA was applied to the networks inferred by existing causal discovery and network inference algorithms applied to a multi-omics data set resulting from a longitudinal study of IBD microbiomes. The most significant unrollings and de-confoundings were manually validated using the existing literature and databases.
ABSTRACT
Abundance profiles from metagenomic sequencing data synthesize information from billions of sequenced reads coming from thousands of microbial genomes. Analyzing and understanding these profiles can be a challenge since the data they represent are complex. Particularly challenging is their visualization, as existing techniques are inadequate when the taxa number is in the thousands. We present a technique, and accompanying software, for the visualization of metagenomic abundance profiles using a space-filling curve that transforms a profile into an interactive 2D image. We created Jasper, an easy to use tool for the visualization and exploration of metagenomic profiles from DNA sequencing data. It orders taxa using a space-filling Hilbert curve, and creates a "Microbiome Map", where each position in the image represents the abundance of a single taxon from a reference collection. Jasper can order taxa in multiple ways, and the resulting microbiome maps can highlight "hot spots" of microbes that are dominant in taxonomic clades or biological conditions. We use Jasper to visualize samples from a variety of microbiome studies, and discuss ways in which microbiome maps can be an invaluable tool to visualize spatial, temporal, disease, and differential profiles. Our approach can create detailed microbiome maps involving hundreds of thousands of microbial reference genomes with the potential to unravel latent relationships (taxonomic, spatio-temporal, functional, and other) that could remain hidden using traditional visualization techniques. The maps can also be converted into animated movies that bring to life the dynamicity of microbiomes.
ABSTRACT
Attention Deficit Hyperactivity Disorder (ADHD) is an increasingly prevalent neuropsychiatric disorder characterized by hyperactivity, inattention, and impulsivity. Symptoms emerge from underlying deficiencies in neurocircuitry, and recent research has suggested a role played by the gut microbiome. The gut microbiome is an ecosystem of interdependent taxa involved in an exponentially complex web of interactions, plus host gene and reaction pathways, some of which involve neurotransmitters with roles in ADHD neurocircuitry. Studies have analyzed the ADHD gut microbiome using macroscale metrics such as diversity and differential abundance, and have proposed several taxa as elevated or reduced in ADHD compared to Control. Few studies have delved into the complex underlying dynamics ultimately responsible for the emergence of such metrics, leaving a largely incomplete, sometimes contradictory, and ultimately inconclusive picture. We aim to help complete this picture by venturing beyond taxa abundances and into taxa relationships (i.e. cooperation and competition), using a publicly available gut microbiome dataset (targeted 16S, v3-4 region, qPCR) from an observational, case-control study of 30 Control (15 female, 15 male) and 28 ADHD (15 female, 13 male) undergraduate students. We first perform the same macroscale analyses prevalent in ADHD gut microbiome literature (diversity, differential abundance, and composition) to observe the degree of correspondence, or any new trends. We then estimate two-way ecological relationships by producing Control and ADHD Microbial Co-occurrence Networks (MCNs), using SparCC correlations (p ≤ 0.01). We perform community detection to find clusters of taxa estimated to mutually cooperate along with their centroids, and centrality calculations to estimate taxa most vital to overall gut ecology. We finally summarize our results, providing conjectures on how they can guide future experiments, some methods for improving our experiments, and general implications for the field.
Subject(s)
Attention Deficit Disorder with Hyperactivity , Gastrointestinal Microbiome , Humans , Female , Male , Gastrointestinal Microbiome/genetics , Case-Control Studies , Ecosystem , BenchmarkingABSTRACT
Introduction. The role of the microbiome in health and disease continues to be increasingly recognized. However, there is significant variability in the bioinformatic protocols for analysing genomic data. This, in part, has impeded the potential incorporation of microbiomics into the clinical setting and has challenged interstudy reproducibility. In microbial compositional analysis, there is a growing recognition for the need to move away from a one-size-fits-all approach to data processing.Gap Statement. Few evidence-based recommendations exist for setting parameters of programs that infer microbiota community profiles despite these parameters significantly impacting the accuracy of taxonomic inference.Aim. To compare three commonly used programs (DADA2, QIIME2, and mothur) and optimize them into four user-adapted pipelines for processing paired-end amplicon reads. We aim to increase the accuracy of compositional inference and help standardize microbiomic protocol.Methods. Two key parameters were isolated across four pipelines: filtering sequence reads based on a whole-number error threshold (maxEE) and truncating read ends based on a quality score threshold (QTrim). Closeness of sample inference was then evaluated using a mock community of known composition.Results. We observed that raw genomic data lost were proportionate to how stringently parameters were set. Exactly how much data were lost varied by pipeline. Accuracy of sample inference correlated with increased sequence read retention. Falsely detected taxa and unaccounted for microbial constituents were unique to pipeline and parameter. Implementation of optimized parameter values led to better approximation of the known mock community.Conclusions. Microbial compositions generated based on the 16S rRNA marker gene should be interpreted with caution. To improve microbial community profiling, bioinformatic protocols must be user-adapted. Analysis should be performed with consideration for the select target amplicon, pipelines and parameters used, and taxa of interest.