An R package for RAIDA can be accessed from http://cals.arizona.edu/%7Eanling/sbg/software.htm.

SUMMARY: Distance-based ordination methods, such as principal coordinates analysis (PCoA), are widely used in the analysis of microbiome data. However, these methods are prone to pose a potential risk of misinterpretation about the compositional difference in samples across different populations if there is a difference in dispersion effects. Accounting for high sparsity and overdispersion of microbiome data, we propose a GLM-based Ordination Method for Microbiome Samples (GOMMS) in this paper. This method uses a zero-inflated quasi-Poisson (ZIQP) latent factor model. An EM algorithm based on the quasi-likelihood is developed to estimate parameters. It performs comparatively to the distance-based approach when dispersion effects are negligible and consistently better when dispersion effects are strong, where the distance-based approach sometimes yields undesirable results. The estimated latent factors from GOMMS can be used to associate the microbiome community with covariates or outcomes using the standard multivariate tests, which can be investigated in future confirmatory experiments. We illustrate the method in simulations and an analysis of microbiome samples from nasopharynx and oropharynx.

Motivated by recent advances in causal mediation analysis and problems in the analysis of microbiome data, we consider the setting where the effect of a treatment on an outcome is transmitted through perturbing the microbial communities or compositional mediators. Compositional and high-dimensional nature of such mediators makes the standard mediation analysis not directly applicable to our setting. We propose a sparse compositional mediation model that can be used to estimate the causal direct and indirect (or mediation) effects utilizing the algebra for compositional data in the simplex space. We also propose tests of total and component-wise mediation effects using bootstrap. We conduct extensive simulation studies to assess the performance of the proposed method and apply the method to a real metagenomic dataset to investigate the effect of fat intake on body mass index mediated through the gut microbiome composition.

Motivated by recent advances in causal mediation analysis and problems in the analysis of microbiome data, we consider the setting where the effect of a treatment on an outcome is transmitted through perturbing the microbial communities or compositional mediators. The compositional and highdimensional nature of such mediators makes the standard mediation analysis not directly applicable to our setting. We propose a sparse compositional mediation model that can be used to estimate the causal direct and indirect (or mediation) effects utilizing the algebra for compositional data in the simplex space. We also propose tests of total and component-wise mediation effects. We conduct extensive simulation studies to assess the performance of the proposed method and apply the method to a real microbiome dataset to investigate an effect of fat intake on body mass index mediated through the gut microbiome.

Staphylococcus aureus is a major cause of ocular infections, often resulting in devastating vision loss. Despite the significant morbidity associated with these infections, little is yet known regarding the specific strain types that may have a predilection for ocular tissues nor the set of virulence factors that drive its pathogenicity in this specific biological niche. Whole genome sequencing (WGS) can provide valuable insight in this regard by providing a prospective, comprehensive assessment of the strain types and virulence factors driving disease among specific subsets of clinical isolates. As such, a set of 163-member S. aureus ocular clinical strains were sequenced and assessed for both common strain types (multilocus sequence type (MLST), spa, agr) associated with ocular infections as well as the presence/absence of 235 known virulence factors in a high throughput manner. This ocular strain set was then directly compared to a fully sequenced 116-member non-ocular S. aureus strain set curated from NCBI in order to identify key differences between ocular and non-ocular S. aureus isolates. The most common sequence types found among ocular S. aureus isolates were ST5, ST8 and ST30, generally reflecting circulating non-ocular pathogenic S. aureus strains. However, importantly, ocular isolates were found to be significantly enriched for a set of enterotoxins, suggesting a potential role for this class of virulence factors in promoting ocular disease. Further genomic analysis revealed that these enterotoxins are located on mobile pathogenicity islands, thus horizontal gene transfer may promote the acquisition of enterotoxins, potentially amplifying S. aureus virulence in ocular tissues.