A Spatial Omnibus Test (SPOT) for Spatial Proteomic Data.

MOTIVATION: Spatial proteomics can reveal the spatial organization of immune cells in the tumor immune microenvironment. Relating measures of spatial clustering, such as Ripley's K or Besag's L, to patient outcomes may offer important clinical insights. However, these measures require pre-specifying a radius in which to quantify clustering, yet no consensus exists on the optimal radius which may be context-specific. RESULTS: We propose a SPatial Omnibus Test (SPOT) which conducts this analysis across a range of candidate radii. At each radius, SPOT evaluates the association between the spatial summary and outcome, adjusting for confounders. SPOT then aggregates results across radii using the Cauchy combination test, yielding an omnibus P-value characterizing the overall degree of association. Using simulations, we verify that the type I error rate is controlled and show SPOT can be more powerful than alternatives. We also apply SPOT to ovarian and lung cancer studies. AVAILABILITY AND IMPLEMENTATION: An R package and tutorial are provided at https://github.com/sarahsamorodnitsky/SPOT.

A multi-bin rarefying method for evaluating alpha diversities in TCR sequencing data.

MOTIVATION: T cell receptors (TCRs) constitute a major component of our adaptive immune system, governing the recognition and response to internal and external antigens. Studying the TCR diversity via sequencing technology is critical for a deeper understanding of immune dynamics. However, library sizes differ substantially across samples, hindering the accurate estimation/comparisons of alpha diversities. To address this, researchers frequently use an overall rarefying approach in which all samples are sub-sampled to an even depth. Despite its pervasive application, its efficacy has never been rigorously assessed. RESULTS: In this paper, we develop an innovative "multi-bin" rarefying approach that partitions samples into multiple bins according to their library sizes, conducts rarefying within each bin for alpha diversity calculations, and performs meta-analysis across bins. Extensive simulations using real-world data highlight the inadequacy of the overall rarefying approach in controlling the confounding effect of library size. Our method proves robust in addressing library size confounding, outperforming competing normalization strategies by achieving better-controlled type-I error rates and enhanced statistical power in association tests. AVAILABILITY AND IMPLEMENTATION: The code is available at https://github.com/mli171/MultibinAlpha. The datasets are freely available at https://doi.org/10.21417/B7001Z and https://doi.org/10.21417/AR2019NC.

Large-scale chemical-genetics yields new M. tuberculosis inhibitor classes.

New antibiotics are needed to combat rising levels of resistance, with new Mycobacterium tuberculosis (Mtb) drugs having the highest priority. However, conventional whole-cell and biochemical antibiotic screens have failed. Here we develop a strategy termed PROSPECT (primary screening of strains to prioritize expanded chemistry and targets), in which we screen compounds against pools of strains depleted of essential bacterial targets. We engineered strains that target 474 essential Mtb genes and screened pools of 100-150 strains against activity-enriched and unbiased compound libraries, probing more than 8.5 million chemical-genetic interactions. Primary screens identified over tenfold more hits than screening wild-type Mtb alone, with chemical-genetic interactions providing immediate, direct target insights. We identified over 40 compounds that target DNA gyrase, the cell wall, tryptophan, folate biosynthesis and RNA polymerase, as well as inhibitors that target EfpA. Chemical optimization yielded EfpA inhibitors with potent wild-type activity, thus demonstrating the ability of PROSPECT to yield inhibitors against targets that would have eluded conventional drug discovery.

Syphilis and cardiovascular risk: a Taiwanese registry.

BACKGROUND AND AIMS: Studies on the impact of syphilis on the cardiovascular system in large populations are limited. This study investigated the effects of syphilis on cardiovascular outcomes. METHODS: Medical records from 2010 to 2015 were retrieved from the Taiwan National Health Insurance Research Database, linked to the Notifiable Infectious Diseases database from the Taiwan Centers for Disease Control. Patients with syphilis were identified, excluding those with missing information, under 20 years of age, or with a history of human immunodeficiency virus infection, acute myocardial infarction, heart failure, aortic regurgitation, replacement of the aortic valve, aneurysm and/or dissection of the aorta, atrial fibrillation, ischaemic stroke, haemorrhagic stroke, and venous thromboembolism. Primary outcomes included new-onset acute myocardial infarction, heart failure, aortic regurgitation, aneurysm and dissection of the aorta, atrial fibrillation, ischaemic stroke, haemorrhagic stroke, venous thromboembolism, cardiovascular death, and all-cause mortality. RESULTS: A total of 28 796 patients with syphilis were identified from 2010 to 2015. After exclusions and frequency matching, 20 601 syphilis patients and 20 601 non-syphilis patients were analysed. The relative rate (RR) was utilized in the analysis, as the competing risk of death was not considered. Compared with patients without syphilis, patients with syphilis had increased risks of acute myocardial infarction (RR 38%, 95% confidence interval [CI] 1.19-1.60, P < .001), heart failure (RR 88%, 95% CI 1.64-2.14, P < .001), aortic regurgitation (RR 81%, 95% CI 1.18-2.75, P = .006), atrial fibrillation (RR 45%, 95% CI 1.20-1.76, P < .001), ischaemic stroke (RR 68%, 95% CI 1.52-1.87, P < .001), haemorrhagic stroke (RR 114%, 95% CI 1.74-2.64, P < .001), venous thromboembolism (RR 67%, 95% CI 1.23-2.26, P = .001), cardiovascular death (RR 155%, 95% CI 2.11-3.08, P < .001), and all-cause death (RR 196%, 95% CI 2.74-3.19, P < .001) but not for aneurysm and dissection of the aorta. CONCLUSIONS: This study demonstrates that patients with syphilis have a higher risk of cardiovascular events and all-cause mortality compared with those without syphilis.

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.

The Study of the Epidemiology of Pediatric Hypertension Registry (SUPERHERO): Rationale and Methods.

Despite increasing prevalence of hypertension in youth and high adult cardiovascular mortality rates, the long-term consequences of youth-onset hypertension remain unknown. This is due to limitations of prior research such as small sample sizes, reliance on manual record review, and limited analytic methods that did not address major biases. The Study of the Epidemiology of Pediatric Hypertension (SUPERHERO) is a multisite retrospective Registry of youth evaluated by subspecialists for hypertension disorders. Sites obtain harmonized electronic health record data using standardized biomedical informatics scripts validated with randomized manual record review. Inclusion criteria are index visit for International Classification of Diseases Diagnostic Codes, 10th Revision (ICD-10 code)-defined hypertension disorder ≥January 1, 2015 and age <19 years. We exclude patients with ICD-10 code-defined pregnancy, kidney failure on dialysis, or kidney transplantation. Data include demographics, anthropomorphics, U.S. Census Bureau tract, histories, blood pressure, ICD-10 codes, medications, laboratory and imaging results, and ambulatory blood pressure. SUPERHERO leverages expertise in epidemiology, statistics, clinical care, and biomedical informatics to create the largest and most diverse registry of youth with newly diagnosed hypertension disorders. SUPERHERO's goals are to (i) reduce CVD burden across the life course and (ii) establish gold-standard biomedical informatics methods for youth with hypertension disorders.

Durability of Original Monovalent mRNA Vaccine Effectiveness Against COVID-19 Omicron-Associated Hospitalization in Children and Adolescents - United States, 2021-2023.

Pediatric COVID-19 vaccination is effective in preventing COVID-19-related hospitalization, but duration of protection of the original monovalent vaccine during SARS-CoV-2 Omicron predominance merits evaluation, particularly given low coverage with updated COVID-19 vaccines. During December 19, 2021-October 29, 2023, the Overcoming COVID-19 Network evaluated vaccine effectiveness (VE) of ≥2 original monovalent COVID-19 mRNA vaccine doses against COVID-19-related hospitalization and critical illness among U.S. children and adolescents aged 5-18 years, using a case-control design. Too few children and adolescents received bivalent or updated monovalent vaccines to separately evaluate their effectiveness. Most case-patients (persons with a positive SARS-CoV-2 test result) were unvaccinated, despite the high frequency of reported underlying conditions associated with severe COVID-19. VE of the original monovalent vaccine against COVID-19-related hospitalizations was 52% (95% CI = 33%-66%) when the most recent dose was administered <120 days before hospitalization and 19% (95% CI = 2%-32%) if the interval was 120-364 days. VE of the original monovalent vaccine against COVID-19-related hospitalization was 31% (95% CI = 18%-43%) if the last dose was received any time within the previous year. VE against critical COVID-19-related illness, defined as receipt of noninvasive or invasive mechanical ventilation, vasoactive infusions, extracorporeal membrane oxygenation, and illness resulting in death, was 57% (95% CI = 21%-76%) when the most recent dose was received <120 days before hospitalization, 25% (95% CI = -9% to 49%) if it was received 120-364 days before hospitalization, and 38% (95% CI = 15%-55%) if the last dose was received any time within the previous year. VE was similar after excluding children and adolescents with documented immunocompromising conditions. Because of the low frequency of children who received updated COVID-19 vaccines and waning effectiveness of original monovalent doses, these data support CDC recommendations that all children and adolescents receive updated COVID-19 vaccines to protect against severe COVID-19.

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.

Testing microbiome association using integrated quantile regression models.

MOTIVATION: Most existing microbiome association analyses focus on the association between microbiome and conditional mean of health or disease-related outcomes, and within this vein, vast computational tools and methods have been devised for standard binary or continuous outcomes. However, these methods tend to be limited either when the underlying microbiome-outcome association occurs somewhere other than the mean level, or when distribution of the outcome variable is irregular (e.g. zero-inflated or mixtures) such that conditional outcome mean is less meaningful. We address this gap by investigating association analysis between microbiome compositions and conditional outcome quantiles. RESULTS: We introduce a new association analysis tool named MiRKAT-IQ within the Microbiome Regression-based Kernel Association Test framework using Integrated Quantile regression models to examine the association between microbiome and the distribution of outcome. For an individual quantile, we utilize the existing kernel machine regression framework to examine the association between that conditional outcome quantile and a group of microbial features (e.g. microbiome community compositions). Then, the goal of examining microbiome association with the whole outcome distribution is achieved by integrating all outcome conditional quantiles over a process, and thus our new MiRKAT-IQ test is robust to both the location of association signals (e.g. mean, variance, median) and the heterogeneous distribution of the outcome. Extensive numerical simulation studies have been conducted to show the validity of the new MiRKAT-IQ test. We demonstrate the potential usefulness of MiRKAT-IQ with applications to actual biological data collected from a previous microbiome study. AVAILABILITY AND IMPLEMENTATION: R codes to implement the proposed methodology is provided in the MiRKAT package, which is available on CRAN. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Effectiveness of Maternal mRNA COVID-19 Vaccination During Pregnancy Against COVID-19-Associated Hospitalizations in Infants Aged <6 Months During SARS-CoV-2 Omicron Predominance - 20 States, March 9, 2022-May 31, 2023.

Infants aged <6 months are not eligible for COVID-19 vaccination. Vaccination during pregnancy has been associated with protection against infant COVID-19-related hospitalization. The Overcoming COVID-19 Network conducted a case-control study during March 9, 2022-May 31, 2023, to evaluate the effectiveness of maternal receipt of a COVID-19 vaccine dose (vaccine effectiveness [VE]) during pregnancy against COVID-19-related hospitalization in infants aged <6 months and a subset of infants aged <3 months. VE was calculated as (1 - adjusted odds ratio) x 100% among all infants aged <6 months and <3 months. Case-patients (infants hospitalized for COVID-19 outside of birth hospitalization and who had a positive SARS-CoV-2 test result) and control patients (infants hospitalized for COVID-19-like illness with a negative SARS-CoV-2 test result) were compared. Odds ratios were determined using multivariable logistic regression, comparing the odds of receipt of a maternal COVID-19 vaccine dose (completion of a 2-dose vaccination series or a third or higher dose) during pregnancy with maternal nonvaccination between case- and control patients. VE of maternal vaccination during pregnancy against COVID-19-related hospitalization was 35% (95% CI = 15%-51%) among infants aged <6 months and 54% (95% CI = 32%-68%) among infants aged <3 months. Intensive care unit admissions occurred in 23% of all case-patients, and invasive mechanical ventilation was more common among infants of unvaccinated (9%) compared with vaccinated mothers (1%) (p = 0.02). Maternal vaccination during pregnancy provides some protection against COVID-19-related hospitalizations among infants, particularly those aged <3 months. Expectant mothers should remain current with COVID-19 vaccination to protect themselves and their infants from hospitalization and severe outcomes associated with COVID-19.

Para-conduit diaphragmatic hernia following esophagectomy-the new price of minimally invasive surgery?

Esophageal Cancer is the seventh commonest cancer worldwide with poor overall survival. Significant morbidity related to open esophagectomy has driven practice toward hybrid, totally minimally invasive and robotic procedures. With the increase in minimally invasive approaches, it has been suggested that there might be an increased incidence of subsequent para-conduit diaphragmatic hernia. To assess the incidence, modifiable risk factors and association with operative approach of this emerging complication, we evaluated outcomes following esophagectomy from two Australian Centers. Prospectively collected databases were examined to identify patients who developed versus did not develop a para-conduit hernia. Patient characteristics, disease factors, treatment factors, operative and post-operative factors were compared for these two groups. A total of 24 of 297 patients who underwent esophagectomy were diagnosed with a symptomatic para-conduit diaphragmatic hernia (8.1%). The significant risk factor for hernia was a minimally invasive abdominal approach (70.8% vs. 35.5%; P = 0.004, odds ratio = 12.876, 95% CI 2.214-74.89). Minimally invasive thoracic approaches were not associated with increased risk. Minimally invasive abdominal approaches to esophagectomy doubled the risk of developing a para-conduit diaphragmatic hernia. Effective operative solutions to address this complication are required.

Thinner biological tissues induce leaflet flutter in aortic heart valve replacements.

Valvular heart disease has recently become an increasing public health concern due to the high prevalence of valve degeneration in aging populations. For patients with severely impacted aortic valves that require replacement, catheter-based bioprosthetic valve deployment offers a minimally invasive treatment option that eliminates many of the risks associated with surgical valve replacement. Although recent percutaneous device advancements have incorporated thinner, more flexible biological tissues to streamline safer deployment through catheters, the impact of such tissues in the complex, mechanically demanding, and highly dynamic valvular system remains poorly understood. The present work utilized a validated computational fluid-structure interaction approach to isolate the behavior of thinner, more compliant aortic valve tissues in a physiologically realistic system. This computational study identified and quantified significant leaflet flutter induced by the use of thinner tissues that initiated blood flow disturbances and oscillatory leaflet strains. The aortic flow and valvular dynamics associated with these thinner valvular tissues have not been previously identified and provide essential information that can significantly advance fundamental knowledge about the cardiac system and support future medical device innovation. Considering the risks associated with such observed flutter phenomena, including blood damage and accelerated leaflet deterioration, this study demonstrates the potentially serious impact of introducing thinner, more flexible tissues into the cardiac system.

Vaginal and Extra-Vaginal Bacterial Colonization and Risk for Incident Bacterial Vaginosis in a Population of Women Who Have Sex With Men.

BACKGROUND: Bacterial vaginosis (BV) is a common cause of vaginal discharge and associated with vaginal acquisition of BV-associated bacteria (BVAB). METHODS: We used quantitative polymerase chain reaction assays to determine whether presence or concentrations of BVAB in the mouth, anus, vagina, or labia before BV predict risk of incident BV in 72 women who have sex with men. RESULTS: Baseline vaginal and extra-vaginal colonization with Gardnerella spp, Megasphaera spp, Sneathia spp, BVAB-2, Dialister sp type 2, and other BVAB was more common among subjects with incident BV. CONCLUSIONS: Prior colonization with BVAB is a consistent risk for BV.

TCR-L: an analysis tool for evaluating the association between the T-cell receptor repertoire and clinical phenotypes.

BACKGROUND: T cell receptors (TCRs) play critical roles in adaptive immune responses, and recent advances in genome technology have made it possible to examine the T cell receptor (TCR) repertoire at the individual sequence level. The analysis of the TCR repertoire with respect to clinical phenotypes can yield novel insights into the etiology and progression of immune-mediated diseases. However, methods for association analysis of the TCR repertoire have not been well developed. METHODS: We introduce an analysis tool, TCR-L, for evaluating the association between the TCR repertoire and disease outcomes. Our approach is developed under a mixed effect modeling, where the fixed effect represents features that can be explicitly extracted from TCR sequences while the random effect represents features that are hidden in TCR sequences and are difficult to be extracted. Statistical tests are developed to examine the two types of effects independently, and then the p values are combined. RESULTS: Simulation studies demonstrate that (1) the proposed approach can control the type I error well; and (2) the power of the proposed approach is greater than approaches that consider fixed effect only or random effect only. The analysis of real data from a skin cutaneous melanoma study identifies an association between the TCR repertoire and the short/long-term survival of patients. CONCLUSION: The TCR-L can accommodate features that can be extracted as well as features that are hidden in TCR sequences. TCR-L provides a powerful approach for identifying association between TCR repertoire and disease outcomes.

Multisystem Inflammatory Syndrome in Children During Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Delta and Omicron Variant Circulation-United States, July 2021-January 2022.

We describe 2116 multisystem inflammatory syndrome in children (MIS-C) cases reported to the Centers for Disease Control and Prevention during Delta and Omicron circulation from July 2021 through January 2022. Half of MIS-C patients were aged 5-11 years, 52% received intensive care unit-level care, and 1.1% died. Only 3.0% of eligible patients were fully vaccinated prior to MIS-C onset.

Trends in Treatments for Multisystem Inflammatory Syndrome in Children (MIS-C), United States, February 2020 - July 2021.

BACKGROUND: Multisystem inflammatory syndrome in children (MIS-C) is a novel severe postinfectious condition associated with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The purpose of this report is to describe nationwide trends in the evolving clinical management of MIS-C. METHODS: Patients with MIS-C were reported from state and local jurisdictions to the Centers for Disease Control and Prevention's (CDC's) MIS-C national surveillance system. Patients' case reports were reviewed to ensure that they met the CDC MIS-C case definition and had sufficient data for analysis. The prevalence of use of treatments for MIS-C, temporal trends in use of these treatments, and frequency of administration of different treatment combinations were analyzed. RESULTS: There were 4470 patients meeting the MIS-C case definition with onset dates from 19 February 2020 to 31 July 2021. The proportion of patients admitted to an intensive care unit (ICU) has declined over time, from 78.7% in April 2020 to 57.5% in June 2021 (P = .001). The most common treatments were intravenous immunoglobulin (IVIG), given to 85.6% of patients; steroids (77.7%), and antiplatelet medications (73.7%); use of each of these treatments has increased over time, particularly in patients not requiring admission to an ICU (all P < .001). Older patients and non-Hispanic Black patients were more likely to receive additional modes of therapy including vasoactive medication, noninvasive respiratory support, anticoagulation medication, and intubation/mechanical ventilation. CONCLUSIONS: IVIG, steroids, and antiplatelet medication have become increasingly utilized as standard treatment for MIS-C patients, while the use of other treatments may be contingent on the type and severity of clinical findings.

Multisystem Inflammatory Syndrome in Children-United States, February 2020-July 2021.

BACKGROUND: Multisystem inflammatory syndrome in children (MIS-C) is a severe hyperinflammatory condition in persons aged <21 years associated with antecedent severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. Our objective was to describe MIS-C cases reported to Centers for Disease Control and Prevention's (CDC's) national surveillance since the coronavirus disease 2019 (COVID-19) pandemic began. METHODS: We included patients meeting the MIS-C case definition with onset date from 19 February 2020 through 31 July 2021, using CDC's MIS-C case report form, which collects information on demographics, clinical presentation, and laboratory results. Trends over time across 3 MIS-C pandemic waves were assessed using Cochran-Armitage test for categorical and Jonckheere-Terpstra test for continuous variables. RESULTS: Of 4901 reported cases, 4470 met inclusion criteria. Median patient age increased over time (P < .001), with a median of 9 years (interquartile range, 5-13 years) during the most recent (third) wave. Male predominance also increased (62% in third wave, P < .001). A significant (P < .001) increase in severe hematologic and gastrointestinal involvement was observed across the study period. Frequency of several cardiovascular complications (ie, cardiac dysfunction, myocarditis, and shock/vasopressor receipt) and renal failure declined (P < .001). Provision of critical care including mechanical ventilation (P < .001) and extracorporeal membrane oxygenation (ECMO; P = .046) decreased, as did duration of hospitalization and mortality (each P < .001). CONCLUSIONS: Over the first 3 pandemic waves of MIS-C in the United States, cardiovascular complications and clinical outcomes including length of hospitalization, receipt of ECMO, and death decreased over time. These data serve as a baseline for monitoring future trends associated with SARS-CoV-2 B.1.617.2 (Delta) or other variants and increased COVID-19 vaccination among children.

A method for subtype analysis with somatic mutations.

MOTIVATION: Cancer is a highly heterogeneous disease, and virtually all types of cancer have subtypes. Understanding the association between cancer subtypes and genetic variations is fundamental to the development of targeted therapies for patients. Somatic mutation plays important roles in tumor development and has emerged as a new type of genetic variations for studying the association with cancer subtypes. However, the low prevalence of individual mutations poses a tremendous challenge to the related statistical analysis. RESULTS: In this article, we propose an approach, subtype analysis with somatic mutations (SASOM), for the association analysis of cancer subtypes with somatic mutations. Our approach tests the association between a set of somatic mutations (from a genetic pathway) and subtypes, while incorporating functional information of the mutations into the analysis. We further propose a robust p-value combination procedure, DAPC, to synthesize statistical significance from different sources. Simulation studies show that the proposed approach has correct type I error and tends to be more powerful than possible alternative methods. In a real data application, we examine the somatic mutations from a cutaneous melanoma dataset, and identify a genetic pathway that is associated with immune-related subtypes. AVAILABILITY AND IMPLEMENTATION: The SASOM R package is available at https://github.com/rksyouyou/SASOM-pkg. R scripts and data are available at https://github.com/rksyouyou/SASOM-analysis. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

MiRKAT: kernel machine regression-based global association tests for the microbiome.

SUMMARY: Distance-based tests of microbiome beta diversity are an integral part of many microbiome analyses. MiRKAT enables distance-based association testing with a wide variety of outcome types, including continuous, binary, censored time-to-event, multivariate, correlated and high-dimensional outcomes. Omnibus tests allow simultaneous consideration of multiple distance and dissimilarity measures, providing higher power across a range of simulation scenarios. Two measures of effect size, a modified R-squared coefficient and a kernel RV coefficient, are incorporated to allow comparison of effect sizes across multiple kernels. AVAILABILITY AND IMPLEMENTATION: MiRKAT is available on CRAN as an R package. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.