Grill, Diane E.

AbstractMotivationOne objective of human microbiome studies is to identify differentially abundant microbes across biological conditions. Previous statistical methods focus on detecting the shift in the abundance and/or prevalence of the microbes and treat the dispersion (spread of the data) as a nuisance. These methods also assume that the dispersion is the same across conditions, an assumption which may not hold in presence of sample heterogeneity. Moreover, the widespread outliers in the microbiome sequencing data make existing parametric models not overly robust. Therefore, a robust and powerful method that allows covariate-dependent dispersion and addresses outliers is still needed for differential abundance analysis.ResultsWe introduce a novel test for differential distribution analysis of microbiome sequencing data by jointly testing the abundance, prevalence and dispersion. The test is built on a zero-inflated negative binomial regression model and winsorized count data to account for zero-inflation and outliers. Using simulated data and real microbiome sequencing datasets, we show that our test is robust across various biological conditions and overall more powerful than previous methods.Availability and implementationR package is available at https://github.com/jchen1981/MicrobiomeDDA.Supplementary informationSupplementary data are available at Bioinformatics online.

AbstractSubstance P (SP), a neurotransmitter released after injury, has been linked to deregulated tissue repair and fibrosis in musculoskeletal tissues and other organs. Although SP inhibition is an effective treatment for nausea, it has not been previously considered as an anti‐fibrotic therapy. Although there are extensive medical records of individuals who have used SP antagonists, our analysis of human registry data revealed that patients receiving these antagonists and arthroplasty are exceedingly rare, thus precluding a clinical evaluation of their potential effects in the context of arthrofibrosis. Therefore, we pursued in vivo studies to assess the effect of SP inhibition early after injury on pro‐fibrotic gene expression and contractures in an animal model of post‐traumatic joint stiffening. Skeletally mature rabbits (n = 24) underwent surgically induced severe joint contracture, while injected with either fosaprepitant (a selective SP antagonist) or saline (control) early after surgery (3, 6, 12, and 24 h). Biomechanical testing revealed that differences in mean contracture angles between the groups were not statistically significant (P = 0.27), suggesting that the drug neither mitigates nor exacerbates joint contracture. However, microarray gene expression analysis revealed that mRNA levels for proteins related to cell signaling, pro‐angiogenic, pro‐inflammatory, and collagen matrix production were significantly different between control and fosaprepitant treated rabbits (P < 0.05). Hence, our study demonstrates that inhibition of SP alters expression of pro‐fibrotic genes in vivo. This finding will motivate future studies to optimize interventions that target SP to reduce the formation of post‐traumatic joint contractures.

AbstractAimsPlatelets retain cytoplasmic messenger RNA and are capable of protein biosynthesis. Several diseases are known to impact the platelet transcriptome but the effect of non‐valvular atrial fibrillation (NVAF) on platelet RNA transcript is essentially unknown. The aim of this study was to evaluate the impact of NVAF on platelet RNA transcript by measuring platelet genes expression in consecutive NVAF patients before and 3‐4 months after pulmonary vein isolation (PVI) and compared to normal sinus rhythm controls (NSR).Methods and ResultsRNA from isolated platelets were reverse transcribed, assayed against 15 genes using real‐time qPCR, and expressed as mean cycle threshold (ΔCt) using beta‐2‐microglobulin as endogenous control. Expression of all evaluated genes, except cathepsin A gene, was significantly lower (higher ΔCt) in 103 NVAF patients compared to 55 NSR controls. Insulin‐like growth factor binding protein acid labile subunit gene (IGFALS) had expression more than 16 fold‐lower (17.0±2.8 vs 12.5±3.8, P<.001), follow by genes encoding for prostacyclin receptor, and for von Willebrand factor which had fourfold lower expression compared to NSR controls. Gender, type of atrial fibrillation, heart failure, hypertension, prior stroke, diabetes mellitus, and atherosclerosis were associated with different gene expression. Following PVI, expression of four genes significantly increased, particularly IGFALS gene (increased 256‐fold) and ADAMT gene increased 16‐fold); expression of three genes significantly decreased, and expression of eight genes has not changed.ConclusionsPlatelets are capable to respond to the circulatory environment of NVAF by altering transcript and changing prothrombotic status. This shows platelet potential for molecular “reprogramming” possibly induced by flow disturbances of NVAF.

175 Background: The efficient delivery of adjuvant and salvage therapy after radical prostatectomy in patients with prostate cancer is hampered by a lack of biomarkers to assess the risk of clinically significant recurrence and progression. Methods: Mayo Clinic Radical Prostatectomy Registry (RP) patient specimens were selected from a case-control cohort with 14 years median follow-up for training and initial validation of an expression biomarker genomic classifier (GC). An independent, blinded case-cohort study of high-risk RP subjects was used to validate GC, comparing the performance of GC to a multivariate logistic regression clinical model (CM) and GC combined with clinical variables (genomic-clinical classifier, GCC) for predicting clinical recurrence (defined as positive bone or CT scan within 5 years after biochemical recurrence). The concordance index (c-index) and Cox model were used to evaluate discrimination and estimate the risk of clinical recurrence. Results: In the training subset (n=359), both GC and GCC had a c-index of 0.90 whereas CM had a c-index of 0.76. In the internal validation set (n=186), GC and GCC had a c-index of 0.76 and 0.75, while CM had a c-index of 0.69. In an independent high-risk study (n=219), GC and GCC had a c-index of 0.77 and 0.76, while CM had a c-index of 0.68. In subset analysis of Gleason score 7 patients within the high-risk group, GC and GCC showed improved discrimination with c-index of 0.78 and 0.76, respectively compared to 0.70 for CM. In the high-risk group, the risk of recurrence by GC model score quartiles at 5 years after RP was estimated at 1%, 5%, 5% and 18%. Conclusions: The GC model shows improved performance over CM in the prediction of clinical recurrence in a high-risk cohort and in subset analysis of Gleason score 7 patients. The addition of clinical variables to the GC model did not significantly contribute to classifier performance in patients with high-risk features. We are further testing the performance of the GC and GCC models and their usefulness in guiding decision-making (e.g., for the adjuvant therapy setting) in additional studies of prostate cancer clinical risk groups.

SummaryInfluenza causes significant morbidity and mortality annually. Although vaccination offers a considerable amount of protection, it is far from perfect, especially in aging populations. This is due to age‐related defects in immune function, a process called immunosenescence. To date, there are no assays or methods to predict or explain variations in an individual's level of response to influenza vaccination. In this study, we measured levels of several immune cell subsets at baseline (Day 0) and at Days 3 and 28 post‐vaccination using flow cytometry. Statistical modelling was performed to assess correlations between levels of cell subsets and Day 28 immune responses – haemagglutination inhibition (HAI) assay, virus neutralizing antibody (VNA) assay, and memory B cell ELISPOT. Changes in several groups of cell types from Day 0 to Day 28 and Day 3 to Day 28 were found to be significantly associated with immune response. Baseline levels of several immune cell subsets, including B cells and regulatory T cells, were able to partially explain variation in memory B‐cell ELISPOT results. Increased expression of HLA‐DR on plasmacytoid dendritic cells after vaccination was correlated with increased HAI and VNA responses. Our data suggest that the expression of activation markers (HLA‐DR and CD86) on various immune cell subsets, as well as the relative distribution of cell subsets, both have value in predicting immune responses to influenza vaccination in older individuals.

Introduction: Comprehensive evaluation of measles-specific humoral immunity after vaccination is important for determining new and/or additional correlates of vaccine immunogenicity and efficacy. Methods: We used a novel proteome microarray technology and statistical modeling to identify factors and models associated with measles-specific functional protective immunity in 150 measles vaccine recipients representing the extremes of neutralizing antibody response after two vaccine doses. Results: Our findings demonstrate a high seroprevalence of antibodies directed to the measles virus (MV) phosphoprotein (P), nucleoprotein (N), as well as antibodies to the large polymerase (L) protein (fragment 1234 to 1900 AA). Antibodies to these proteins, in addition to anti-F antibodies (and, to a lesser extent, anti-H antibodies), were correlated with neutralizing antibody titer and/or were associated with and predictive of neutralizing antibody response. Conclusion: Our results identify antibodies to specific measles virus proteins and statistical models for monitoring and assessment of measles-specific functional protective immunity in vaccinated individuals.

Objectives The goal of this study was to determine whether intra-articular administration of the potentially anti-fibrotic agent decorin influences the expression of genes involved in the fibrotic cascade, and ultimately leads to less contracture, in an animal model. Methods A total of 18 rabbits underwent an operation on their right knees to form contractures. Six limbs in group 1 received four intra-articular injections of decorin; six limbs in group 2 received four intra-articular injections of bovine serum albumin (BSA) over eight days; six limbs in group 3 received no injections. The contracted limbs of rabbits in group 1 were biomechanically and genetically compared with the contracted limbs of rabbits in groups 2 and 3, with the use of a calibrated joint measuring device and custom microarray, respectively. Results There was no statistical difference in the flexion contracture angles between those limbs that received intra-articular decorin versus those that received intra-articular BSA (66° vs 69°; p = 0.41). Likewise, there was no statistical difference between those limbs that received intra-articular decorin versus those who had no injection (66° vs 72°; p = 0.27). When compared with BSA, decorin led to a statistically significant increase in the mRNA expression of 12 genes (p < 0.01). In addition, there was a statistical change in the mRNA expression of three genes, when compared with those without injection. Conclusions In this model, when administered intra-articularly at eight weeks, 2 mg of decorin had no significant effect on joint contractures. However, our genetic analysis revealed a significant alteration in several fibrotic genes. Cite this article: Bone Joint Res 2014;3:82–8.

Abstract Background: First in human studies of Z-endoxifen hydrochloride (E), the active metabolite of tamoxifen (T), are underway in metastatic breast cancer (BC). Previous data have demonstrated the superiority of aromatase inhibitors (AI's) over T in estrogen receptor (ER) + BC. Using an in vivo aromatase expressing model (MCF7/AC1), we compared the antitumor activity of E with T and Letrozole (L), as well as the antitumor activity and global gene expression changes of E with T in an L-resistant model. Methods: MCF7/AC1 tumors were stimulated with androstenedione. Once tumor size reached 300 mm3, mice (30/group) were randomly assigned to one of five treatment groups: control (daily, po), T (500 μg/day, sc), endoxifen 25 mg/kg/day p.o.(LDE) endoxifen 75 mg/kg/day p.o. (HDE) or letrozole, 10 μg/day s.c for 4 weeks. Tumors were harvested from control, T, and E groups while the L group continued treatment until the development of resistance defined as an increase in tumor volume of at least 300% from day 1. Mice with L-resistant tumors were randomly assigned to T (n=4) or E (n=5) for 4 weeks and then sacrificed. Gene expression in L-resistant tumors was quantified using Affymetrix U133+2 and changes in gene expression profiles [comparing T and E with L-resistant (n=3)] were analyzed. Genes identified as significantly different were confirmed by real-time RT-PCR assays. Results: At the 4 week time point, both doses of E and L resulted in greater anti-tumor activity than control (Wilcoxon rank sum test: all p &lt; 0.0001); however, tumor burden did not differ between T and control (p=0.095). HDE resulted in significantly less tumor burden than T (p=0.002) but was similar to L. In mice that continued on L, resistance developed at 24 weeks in 9/25 mice. These mice were randomly assigned to either T (n=4) or E (n=5) for 4 weeks. Tumor volume (expressed as a% of its size prior to randomization) was significantly different comparing E (73.3%; range: 69.3 to 80.75%) versus T (148.39%; range: 114.07 to 165.99%) (Wilcoxon rank sum test p=0.016). Compared to control, microarray studies identified 1518 unique probe sets regulated by E (p&lt;0.001) compared to 441 for T including estrogen-regulated genes such as progesterone receptor (PGR) and amphiregulin (AREG) that were significantly down-regulated in the E group [PGR (−6.2 fold, p=0.000008) and AREG (−3.2 fold, p=.0006) but unchanged or up-regulated in the T group (PGR unchanged and AREG +9.2 fold p=0.00002). These findings were confirmed by RT-PCR. Conclusions: Using the MCF7/AC1 model previously used to show the superiority of AI's over T, HDE demonstrated similar antitumor activity to L and was superior to T. In cells resistant to L, E was superior to T and gene expression changes demonstrate that E down-regulates while T activates estrogen regulated genes. These findings support the ongoing development of E for the treatment of ER+ BC. Citation Information: Cancer Res 2011;71(24 Suppl):Abstract nr PD01-06.