Viral genomic variation and the severity of genital HSV-2 infection as quantified by shedding rate: a viral genome-wide association study.

BACKGROUND: The clinical severity of genital HSV-2 infection varies widely among infected persons with some experiencing frequent genital lesions while others are asymptomatic. The viral genital shedding rate is closely associated with and has been established as a surrogate marker of clinical severity. METHODS: To assess the relationship between viral genetics and shedding, we assembled a set of 145 persons who had the severity of their genital herpes quantified through determination of their HSV genital shedding rate. An HSV-2 sample from each person was sequenced and biallelic variants among these genomes were identified. RESULTS: We found no association between metrics of genome-wide variation in HSV-2 and shedding rate. A viral genome-wide association study (vGWAS) identified the minor alleles of three individual unlinked variants as significantly associated with higher shedding rate (p<8.4x10-5): C44973T (A512T), a non-synonymous variant in UL22 (glycoprotein H); A74534G, a synonymous variant in UL36 (large tegument protein); and T119283C, an intergenic variant. We also found an association between the total number of minor alleles for the significant variants and shedding rate (p=6.6x10-7). CONCLUSIONS: These results add to a growing body of literature for HSV suggesting a connection between viral genetic variation and clinically important phenotypes of infection.

Limitation of permutation-based differential correlation analysis.

The comparison of biological systems, through the analysis of molecular changes under different conditions, has played a crucial role in the progress of modern biological science. Specifically, differential correlation analysis (DCA) has been employed to determine whether relationships between genomic features differ across conditions or outcomes. Because ascertaining the null distribution of test statistics to capture variations in correlation is challenging, several DCA methods utilize permutation which can loosen parametric (e.g., normality) assumptions. However, permutation is often problematic for DCA due to violating the assumption that samples are exchangeable under the null. Here, we examine the limitations of permutation-based DCA and investigate instances where the permutation-based DCA exhibits poor performance. Experimental results show that the permutation-based DCA often fails to control the type I error under the null hypothesis of equal correlation structures.

Accommodating multiple potential normalizations in microbiome associations studies.

BACKGROUND: Microbial communities are known to be closely related to many diseases, such as obesity and HIV, and it is of interest to identify differentially abundant microbial species between two or more environments. Since the abundances or counts of microbial species usually have different scales and suffer from zero-inflation or over-dispersion, normalization is a critical step before conducting differential abundance analysis. Several normalization approaches have been proposed, but it is difficult to optimize the characterization of the true relationship between taxa and interesting outcomes.  RESULTS: To avoid the challenge of picking an optimal normalization and accommodate the advantages of several normalization strategies, we propose an omnibus approach. Our approach is based on a Cauchy combination test, which is flexible and powerful by aggregating individual p values. We also consider a truncated test statistic to prevent substantial power loss. We experiment with a basic linear regression model as well as recently proposed powerful association tests for microbiome data and compare the performance of the omnibus approach with individual normalization approaches. Experimental results show that, regardless of simulation settings, the new approach exhibits power that is close to the best normalization strategy, while controling the type I error well.  CONCLUSIONS: The proposed omnibus test releases researchers from choosing among various normalization methods and it is an aggregated method that provides the powerful result to the underlying optimal normalization, which requires tedious trial and error. While the power may not exceed the best normalization, it is always much better than using a poor choice of normalization.

Association between vaginal microbiota and vaginal inflammatory immune markers in postmenopausal women.

OBJECTIVE: In premenopausal individuals, vaginal microbiota diversity and lack of Lactobacillus dominance are associated with greater mucosal inflammation, which is linked to a higher risk of cervical dysplasia and infections. It is not known if the association between the vaginal microbiota and inflammation is present after menopause, when the vaginal microbiota is generally higher-diversity and fewer people have Lactobacillus dominance. METHODS: This is a post hoc analysis of a subset of postmenopausal individuals enrolled in a randomized trial for treatment of moderate-severe vulvovaginal discomfort that compared vaginal moisturizer, estradiol, or placebo. Vaginal fluid samples from 0, 4, and 12 weeks were characterized using 16S rRNA gene sequencing (microbiota) and MesoScale Discovery (vaginal fluid immune markers: IL-1b, IL-1a, IL-2, IL-6, IL-18, IL-10, IL-9, IL-13, IL-8, IP10, MIP1a, MIP1b, MIP3a). Global associations between cytokines and microbiota (assessed by relative abundance of individual taxa and Shannon index for alpha, or community, diversity) were explored, adjusting for treatment arm, using linear mixed models, principal component analysis, and Generalized Linear Mixed Model + Microbiome Regression-based Kernel Association Test (GLMM-MiRKAT). RESULTS: A total of 119 individuals with mean age of 61 years were included. At baseline, 29.5% of participants had a Lactobacillus -dominant vaginal microbiota. Across all timepoints, alpha diversity (Shannon index, P = 0.003) was highly associated with immune markers. Individual markers that were associated with Lactobacillus dominance were similar to those observed in premenopausal people: IL-10, IL-1b, IL-6, IL-8 (false discovery rate [FDR] < 0.01), IL-13 (FDR = 0.02), and IL-2 (FDR = 0.09). Over 12 weeks, change in alpha diversity was associated with change in cytokine concentration (Shannon, P = 0.018), with decreased proinflammatory cytokine concentrations observed with decreasing alpha diversity. CONCLUSIONS: In this cohort of postmenopausal individuals, Lactobacillus dominance and lower alpha diversity were associated with lower concentrations of inflammatory immune markers, as has been reported in premenopausal people. This suggests that after menopause lactobacilli continue to have beneficial effects on vaginal immune homeostasis, despite lower prevalence.

A fast kernel independence test for cluster-correlated data.

Cluster-correlated data receives a lot of attention in biomedical and longitudinal studies and it is of interest to assess the generalized dependence between two multivariate variables under the cluster-correlated structure. The Hilbert-Schmidt independence criterion (HSIC) is a powerful kernel-based test statistic that captures various dependence between two random vectors and can be applied to an arbitrary non-Euclidean domain. However, the existing HSIC is not directly applicable to cluster-correlated data. Therefore, we propose a HSIC-based test of independence for cluster-correlated data. The new test statistic combines kernel information so that the dependence structure in each cluster is fully considered and exhibits good performance under high dimensions. Moreover, a rapid p value approximation makes the new test fast applicable to large datasets. Numerical studies show that the new approach performs well in both synthetic and real world data.

Batch effects removal for microbiome data via conditional quantile regression.

Batch effects in microbiome data arise from differential processing of specimens and can lead to spurious findings and obscure true signals. Strategies designed for genomic data to mitigate batch effects usually fail to address the zero-inflated and over-dispersed microbiome data. Most strategies tailored for microbiome data are restricted to association testing or specialized study designs, failing to allow other analytic goals or general designs. Here, we develop the Conditional Quantile Regression (ConQuR) approach to remove microbiome batch effects using a two-part quantile regression model. ConQuR is a comprehensive method that accommodates the complex distributions of microbial read counts by non-parametric modeling, and it generates batch-removed zero-inflated read counts that can be used in and benefit usual subsequent analyses. We apply ConQuR to simulated and real microbiome datasets and demonstrate its advantages in removing batch effects while preserving the signals of interest.

Preventing Failures by Dataset Shift Detection in Safety-Critical Graph Applications.

Dataset shift refers to the problem where the input data distribution may change over time (e.g., between training and test stages). Since this can be a critical bottleneck in several safety-critical applications such as healthcare, drug-discovery, etc., dataset shift detection has become an important research issue in machine learning. Though several existing efforts have focused on image/video data, applications with graph-structured data have not received sufficient attention. Therefore, in this paper, we investigate the problem of detecting shifts in graph structured data through the lens of statistical hypothesis testing. Specifically, we propose a practical two-sample test based approach for shift detection in large-scale graph structured data. Our approach is very flexible in that it is suitable for both undirected and directed graphs, and eliminates the need for equal sample sizes. Using empirical studies, we demonstrate the effectiveness of the proposed test in detecting dataset shifts. We also corroborate these findings using real-world datasets, characterized by directed graphs and a large number of nodes.