Cancer treatment decisions are increasingly guided by which specific genes are mutated within each patient’s tumor. For example, agents inhibiting the epidermal growth factor receptor (EGFR) benefit many colorectal cancer (CRC) patients, with the general exception of those whose tumor includes a KRAS mutation. However, among the various KRAS mutations, that which encodes the G13D mutant protein (KRASG13D) behaves differently; for unknown reasons, KRASG13D CRC patients benefit from the EGFR-blocking antibody cetuximab. Controversy surrounds this observation, because it contradicts the well-established mechanisms of EGFR signaling with regard to RAS mutations. Here, we identified a systems-level, mechanistic explanation for why KRASG13D cancers respond to EGFR inhibition. A computational model of RAS signaling revealed that the biophysical differences between the three most common KRAS mutants were sufficient to generate different sensitivities to EGFR inhibition. Integrated computation with experimentation then revealed a nonintuitive, mutant-specific dependency of wild-type RAS activation by EGFR that is determined by the interaction strength between KRAS and the tumor suppressor neurofibromin (NF1). KRAS mutants that strongly interacted with and competitively inhibited NF1 drove wild-type RAS activation in an EGFR-independent manner, whereas KRASG13D weakly interacted with and could not competitively inhibit NF1 and, thus, KRASG13D cells remained dependent on EGFR for wild-type RAS activity. Overall, our work demonstrates how systems approaches enable mechanism-based inference in genomic medicine and can help identify patients for selective therapeutic strategies.

Background Microbial source tracking methods are used to determine the origin of contaminating bacteria and other microorganisms, particularly in contaminated water systems. The Bayesian SourceTracker approach uses deep-sequencing marker gene libraries (16S ribosomal RNA) to determine the proportional contributions of bacteria from many potential source environments to a given sink environment simultaneously. Since its development, SourceTracker has been applied to an extensive diversity of studies, from beach contamination to human behavior. Methods Here, we demonstrate a novel application of SourceTracker to work with metagenomic datasets and tested this approach using sink samples from a study of coastal marine environments. Source environment metagenomes were obtained from metagenomics studies of gut, freshwater, marine, sand and soil environments. As part of this effort, we implemented features for determining the stability of source proportion estimates, including precision visualizations for performance optimization, and performed domain-specific source-tracking analyses (i.e., Bacteria, Archaea, Eukaryota and viruses). We also applied SourceTracker to metagenomic libraries generated from samples collected from the International Space Station (ISS). Results SourceTracker proved highly effective at predicting the composition of known sources using shotgun metagenomic libraries. In addition, we showed that different taxonomic domains sometimes presented highly divergent pictures of environmental source origins for both the coastal marine and ISS samples. These findings indicated that applying SourceTracker to separate domains may provide a deeper understanding of the microbial origins of complex, mixed-source environments, and further suggested that certain domains may be preferable for tracking specific sources of contamination.

Compositional data analysis (CoDA) methods have increased in popularity as a new framework for analyzing next-generation sequencing (NGS) data. CoDA methods, such as the centered log-ratio (clr) transformation, adjust for the compositional nature of NGS counts, which is not addressed by traditional normalization methods. CoDA has only been sparsely applied to NGS data generated from microbial communities or to multiple ‘omics’ datasets. In this study, we applied CoDA methods to analyze NGS and untargeted metabolomic datasets obtained from bacterial and fungal communities. Specifically, we used clr transformation to reanalyze NGS amplicon and metabolomics data from a study investigating the effects of building material type, moisture and time on microbial and metabolomic diversity. Compared to analysis of untransformed data, analysis of clr-transformed data revealed novel relationships and stronger associations between sample conditions and microbial and metabolic community profiles.

Introduction. Tobacco smoke contains numerous toxic chemicals that accumulate in indoor environments creating thirdhand smoke (THS). We investigated if THS-polluted homes differed in children’s human and built environment microbiomes as compared to THS-free homes. Methods. Participants were N=19 THS exposed children and N=10 unexposed children (≤ 5 years) and their parents. Environmental and biological samples were analyzed for THS pollutants and exposure. Swab samples were collected from the built environment (floor, table, armrest, bedframe) and child (finger, nose, mouth, and ear canal) and 16S ribosomal RNA genes were analyzed for bacterial taxa using high-throughput DNA sequencing. Results. Phylogenetic α-diversity was significantly higher for the built environmental microbiomes in THS-polluted homes compared to THS-free homes (p<0.014). Log2 fold comparison found differences between THS-polluted and THS-free homes for specific genera at built environment (e.g., Acinetobacter, Bradyrhizobium, Corynebacterium, Gemella, Neisseria, Staphylococcus, Streptococcus, Veillonella) and in samples from children s (esp. Corynebacterium, Gemella, Lautropia, Neisseria, Rothia, Staphylococcus, and Veillonella). Conclusion. When exposed to THS, indoor and children microbiomes are altered in an environment-specific manner. Changes are similar to those reported in previous studies for smokers and secondhand smoke-exposed persons. THS-induced changes in child and environmental microbiome may play a role in clinical outcomes in children.

Periodontal disease (PD) is a chronic, progressive polymicrobial disease that induces a strong host immune response. Culture-independent methods, such as next-generation sequencing (NGS) of bacteria 16S amplicon and shotgun metagenomic libraries, have greatly expanded our understanding of PD biodiversity, identified novel PD microbial associations, and shown that PD biodiversity increases with pocket depth. NGS studies have also found PD communities to be highly host-specific in terms of both biodiversity and the response of microbial communities to periodontal treatment. As with most microbiome work, the majority of PD microbiome studies use standard data normalization procedures that do not account for the compositional nature of NGS microbiome data. Here, we apply recently developed compositional data analysis (CoDA) approaches and software tools to reanalyze multiomics (16S, metagenomics, and metabolomics) data generated from previously published periodontal disease studies. CoDA methods, such as centered log-ratio (clr) transformation, compensate for the compositional nature of these data, which can not only remove spurious correlations but also allows for the identification of novel associations between microbial features and disease conditions. We validated many of the studies’ original findings, but also identified new features associated with periodontal disease, including the genera Schwartzia and Aerococcus and the cytokine C-reactive protein (CRP). Furthermore, our network analysis revealed a lower connectivity among taxa in deeper periodontal pockets, potentially indicative of a more “random” microbiome. Our findings illustrate the utility of CoDA techniques in multiomics compositional data analysis of the oral microbiome.