Accelerating Bayesian hierarchical clustering of time series data with a randomised algorithm

We live in an era of abundant data. This has necessitated the development of new and innovative statistical algorithms to get the most from experimental data. For example, faster algorithms make practical the analysis of larger genomic data sets, allowing us to extend the utility of cutting-edge statistical methods. We present a randomised algorithm that accelerates the clustering of time series data using the Bayesian Hierarchical Clustering (BHC) statistical method. BHC is a general method for clustering any discretely sampled time series data. In this paper we focus on a particular application to microarray gene expression data. We define and analyse the randomised algorithm, before presenting results on both synthetic and real biological data sets. We show that the randomised algorithm leads to substantial gains in speed with minimal loss in clustering quality. The randomised time series BHC algorithm is available as part of the R package BHC, which is available for download from Bioconductor (version 2.10 and above) via http://bioconductor.org/packages/2.10/bioc/html/BHC.html. We have also made available a set of R scripts which can be used to reproduce the analyses carried out in this paper. These are available from the following URL. https://sites.google.com/site/randomisedbhc/

Gene set of nuclear-encoded mitochondrial regulators is enriched for common inherited variation in obesity

There are hints of an altered mitochondrial function in obesity. Nuclear-encoded genes are relevant for mitochondrial function (3 gene sets of known relevant pathways: (1) 16 nuclear regulators of mitochondrial genes, (2) 91 genes for oxidative phosphorylation and (3) 966 nuclear-encoded mitochondrial genes). Gene set enrichment analysis (GSEA) showed no association with type 2 diabetes mellitus in these gene sets. Here we performed a GSEA for the same gene sets for obesity. Genome wide association study (GWAS) data from a case-control approach on 453 extremely obese children and adolescents and 435 lean adult controls were used for GSEA. For independent confirmation, we analyzed 705 obesity GWAS trios (extremely obese child and both biological parents) and a population-based GWAS sample (KORA F4, n = 1,743). A meta-analysis was performed on all three samples. In each sample, the distribution of significance levels between the respective gene set and those of all genes was compared using the leading-edge-fraction-comparison test (cut-offs between the 50(th) and 95(th) percentile of the set of all gene-wise corrected p-values) as implemented in the MAGENTA software. In the case-control sample, significant enrichment of associations with obesity was observed above the 50(th) percentile for the set of the 16 nuclear regulators of mitochondrial genes (p(GSEA,50) = 0.0103). This finding was not confirmed in the trios (p(GSEA,50) = 0.5991), but in KORA (p(GSEA,50) = 0.0398). The meta-analysis again indicated a trend for enrichment (p(MAGENTA,50) = 0.1052, p(MAGENTA,75) = 0.0251). The GSEA revealed that weak association signals for obesity might be enriched in the gene set of 16 nuclear regulators of mitochondrial genes

Emergence of Double- and Triple-Gene Reassortant G1P[8]Rotaviruses Possessing a DS-1-Like Backbone after RotavirusVaccine Introduction in Malawi

To combat the high burden of rotavirus gastroenteritis, multiple African countries have introduced rotavirus vaccines into their childhood immunisation programmes. Malawi incorporated a G1P[8] rotavirus vaccine (Rotarix™) into its immunisation schedule in 2012. Utilising a surveillance platform of hospitalised rotavirus gastroenteritis cases, we examined the phylodynamics of G1P[8] rotavirus strains that circulated in Malawi before (1998 - 2012) and after (2013 - 2014) vaccine introduction. Analysis of whole genomes obtained through next generation sequencing revealed that all randomly-selected pre-vaccine G1P[8] strains sequenced (n=32) possessed a Wa-like genetic constellation, whereas post-vaccine G1P[8] strains (n=18) had a DS-1-like constellation. Phylodynamic analyses indicated that post-vaccine G1P[8] strains emerged through reassortment events between human Wa- and DS-1-like rotaviruses that circulated in Malawi from the 1990's, hence classified as atypical DS-1-like reassortants. The time to the most recent common ancestor for G1P[8] strains was from 1981-1994; their evolutionary rates ranged from 9.7 x 10(-4)-4.1 x 10(-3) nucleotide/substitutions/site/year. Three distinct G1P[8] lineages chronologically replaced each other between 1998 and 2014. Genetic drift was the likely driver for lineage turnover in 2005, whereas replacement in 2013 was due to reassortment. Amino acid substitution within the outer glycoprotein VP7 of G1P[8] strains had no impact on the structural conformation of the antigenic regions, suggesting that it is unlikely that they would affect recognition by vaccine-induced neutralizing antibodies. While the emergence of DS-1-like G1P[8] rotavirus reassortants in Malawi was therefore likely due to natural genotype variation, vaccine effectiveness against such strains needs careful evaluation.ImportanceThe error-prone RNA-dependent RNA polymerase and the segmented RNA genome predispose rotaviruses to genetic mutation and genome reassortment, respectively. These evolutionary mechanisms generate novel strains and have the potential to lead to the emergence of vaccine-escape mutants. While multiple African countries have introduced rotavirus vaccine, there are few data describing the evolution of rotaviruses that circulated before and after vaccine introduction. We report the emergence of atypical DS-1-like G1P[8] strains during the post-vaccine era in Malawi. Three distinct G1P[8] lineages circulated chronologically from 1998-2014; mutation and reassortment drove lineage turnover in 2005 and 2013, respectively. Amino acid substitutions within the outer capsid VP7 glycoprotein did not affect the structural conformation of mapped antigenic sites, suggesting limited effect in recognition of G1 specific vaccine-derived antibodies. The genes that constitute the remaining genetic backbone may play important roles in immune evasion, and vaccine effectiveness against such atypical strains needs careful evaluation

The Universal One-Loop Effective Action

We present the universal one-loop effective action for all operators of dimension up to six obtained by integrating out massive, non-degenerate multiplets. Our general expression may be applied to loops of heavy fermions or bosons, and has been checked against partial results available in the literature. The broad applicability of this approach simplifies one-loop matching from an ultraviolet model to a lower-energy effective field theory (EFT), a procedure which is now reduced to the evaluation of a combination of matrices in our universal expression, without any loop integrals to evaluate. We illustrate the relationship of our results to the Standard Model (SM) EFT, using as an example the supersymmetric stop and sbottom squark Lagrangian and extracting from our universal expression the Wilson coefficients of dimension-six operators composed of SM fields.Comment: 30 pages, v2 contains additional comments and corrects typos, version accepted for publication in JHE

Quantum Gravity in 2+1 Dimensions: The Case of a Closed Universe

In three spacetime dimensions, general relativity drastically simplifies, becoming a ``topological'' theory with no propagating local degrees of freedom. Nevertheless, many of the difficult conceptual problems of quantizing gravity are still present. In this review, I summarize the rather large body of work that has gone towards quantizing (2+1)-dimensional vacuum gravity in the setting of a spatially closed universe.Comment: 61 pages, draft of review for Living Reviews; comments, criticisms, additions, missing references welcome; v2: minor changes, added reference

Emergence of double- and triple-gene reassortant G1P[8] rotaviruses possessing a DS-1-like backbone after rotavirus vaccine introduction in Malawi

To combat the high burden of rotavirus gastroenteritis, multiple African countries have introduced rotavirus vaccines into their childhood immunization programs. Malawi incorporated a G1P[8] rotavirus vaccine (Rotarix) into its immunization schedule in 2012. Utilizing a surveillance platform of hospitalized rotavirus gastroenteritis cases, we examined the phylodynamics of G1P[8] rotavirus strains that circulated in Malawi before (1998 to 2012) and after (2013 to 2014) vaccine introduction. Analysis of whole genomes obtained through next-generation sequencing revealed that all randomly selected prevaccine G1P[8] strains sequenced (n = 32) possessed a Wa-like genetic constellation, whereas postvaccine G1P[8] strains (n = 18) had a DS-1-like constellation. Phylodynamic analyses indicated that postvaccine G1P[8] strains emerged through reassortment events between human Wa- and DS-1-like rotaviruses that circulated in Malawi from the 1990s and hence were classified as atypical DS-1-like reassortants. The time to the most recent common ancestor for G1P[8] strains was from 1981 to 1994; their evolutionary rates ranged from 9.7 × 10−4 to 4.1 × 10−3 nucleotide substitutions/site/year. Three distinct G1P[8] lineages chronologically replaced each other between 1998 and 2014. Genetic drift was the likely driver for lineage turnover in 2005, whereas replacement in 2013 was due to reassortment. Amino acid substitution within the outer glycoprotein VP7 of G1P[8] strains had no impact on the structural conformation of the antigenic regions, suggesting that it is unlikely that they would affect recognition by vaccine-induced neutralizing antibodies. While the emergence of DS-1-like G1P[8] rotavirus reassortants in Malawi was therefore likely due to natural genotype variation, vaccine effectiveness against such strains needs careful evaluation. IMPORTANCE: The error-prone RNA-dependent RNA polymerase and the segmented RNA genome predispose rotaviruses to genetic mutation and genome reassortment, respectively. These evolutionary mechanisms generate novel strains and have the potential to lead to the emergence of vaccine escape mutants. While multiple African countries have introduced a rotavirus vaccine, there are few data describing the evolution of rotaviruses that circulated before and after vaccine introduction. We report the emergence of atypical DS-1-like G1P[8] strains during the postvaccine era in Malawi. Three distinct G1P[8] lineages circulated chronologically from 1998 to 2014; mutation and reassortment drove lineage turnover in 2005 and 2013, respectively. Amino acid substitutions within the outer capsid VP7 glycoprotein did not affect the structural conformation of mapped antigenic sites, suggesting a limited effect on the recognition of G1-specific vaccine-derived antibodies. The genes that constitute the remaining genetic backbone may play important roles in immune evasion, and vaccine effectiveness against such atypical strains needs careful evaluation

“I Think I Became a Swimmer Rather than Just Someone with a Disability Swimming Up and Down”: Paralympic Athletes Perceptions of Self and Identity Development

This is an Accepted Manuscript of an article published by Taylor & Francis Group in Disability and Rehabilitation on 27 September 2016, available online at:DOI: https://doi.org/10.1080/09638288.2016.1217074.Purpose: The purpose of this study was to explore the role of swimming on Paralympic athletes’ perceptions of self and identity development. Method: A hermeneutic phenomenological approach was taken. During semi-structured interviews five Paralympic swimmers (aged 20-24 years) were asked questions about their swimming career, perceptions of self, integration, and impairment. Interviews were audio-recorded and transcribed verbatim. Results: An Interpretative Phenomenological Analysis1 yielded three superordinate themes: a) ‘One of the crowd’; none of the participants viewed themselves as disabled, nor as supercrips; these perceptions stemmed from family-, school-, and swimming- related experiences, b) ‘Becoming me’; participation in swimming facilitated self- and social-acceptance, and identity development, and c) ‘A badge of honour’; swimming presented opportunity to present and reinforce a positive identity. Conclusions: Swimming experiences enabled the participants to enhance personal and social identities, integrate through pro-social mechanisms, and to develop a career path following retirement from competition.through pro-social mechanisms, and to develop a career path following retirement from competition.Peer reviewe

Self-reported medication side effects in an older cohort living independently in the community - the Melbourne Longitudinal Study on Health Ageing (MELSHA) : cross-sectional analysis of prevalence and risk factors

Background Medication side effects are an important cause of morbidity, mortality and costs in older people. The aim of our study was to examine prevalence and risk factors for self-reported medication side effects in an older cohort living independently in the community.Methods The Melbourne Longitudinal Study on Healthy Ageing (MELSHA), collected information on those aged 65 years or older living independently in the community and commenced in 1994. Data on medication side effects was collected from the baseline cohort (n = 1000) in face-to-face baseline interviews in 1994 and analysed as cross-sectional data. Risk factors examined were: socio-demographics, health status and medical conditions; medication use and health service factors. Analysis included univariate logistic regression to estimate unadjusted risk and multivariate logistic regression analysis to assess confounding and estimate adjusted risk.Results Self-reported medication side effects were reported by approximately 6.7% (67/1000) of the entire baseline MELSHA cohort, and by 8.5% (65/761) of those on medication. Identified risk factors were increased education level, co-morbidities and health service factors including recency of visiting the pharmacist, attending younger doctors, and their doctor\u27s awareness of their medications. The greatest increase in risk for medication side effects was associated with liver problems and their doctor\u27s awareness of their medications. Aging and gender were not risk factors.Conclusion Prevalence of self-reported medication side effects was comparable with that reported in adults attending General Practices in a primary care setting in Australia. The prevalence and identified risk factors provide further insight and opportunity to develop strategies to address the problem of medication side effects in older people living independently in the community setting. <br /

ILC3 function as a double-edged sword in inflammatory bowel diseases

Inflammatory bowel diseases (IBD), composed mainly of Crohn’s disease (CD) and ulcerative colitis (UC), are strongly implicated in the development of intestinal inflammation lesions. Its exact etiology and pathogenesis are still undetermined. Recently accumulating evidence supports that group 3 innate lymphoid cells (ILC3) are responsible for gastrointestinal mucosal homeostasis through moderate generation of IL-22, IL-17, and GM-CSF in the physiological state. ILC3 contribute to the progression and aggravation of IBD while both IL-22 and IL-17, along with IFN-γ, are overexpressed by the dysregulation of NCR− ILC3 or NCR+ ILC3 function and the bias of NCR+ ILC3 towards ILC1 as well as regulatory ILC dysfunction in the pathological state. Herein, we feature the group 3 innate lymphoid cells’ development, biological function, maintenance of gut homeostasis, mediation of IBD occurrence, and potential application to IBD therapy

Nonsecretor Histo-blood Group Antigen Phenotype Is Associated With Reduced Risk of Clinical Rotavirus Vaccine Failure in Malawian Infants

Background Histo–blood group antigen (HBGA) Lewis/secretor phenotypes predict genotype-specific susceptibility to rotavirus gastroenteritis (RVGE). We tested the hypothesis that nonsecretor/Lewis-negative phenotype leads to reduced vaccine take and lower clinical protection following vaccination with G1P[8] rotavirus vaccine (RV1) in Malawian infants Methods A cohort study recruited infants receiving RV1 at age 6 and 10 weeks. HBGA phenotype was determined by salivary enzyme-linked immunosorbent assay (ELISA). RV1 vaccine virus shedding was detected by quantitative real-time polymerase chain reaction (qRT-PCR) in stool collected on alternate days for 10 days post-immunization. Plasma rotavirus–specific immunoglobulin A was determined by ELISA pre- and post-immunization. In a case-control study, HBGA phenotype distribution was compared between RV1-vaccinated infants with RVGE and 1:1 age-matched community controls. Rotavirus genotype was determined by RT-PCR. Results In 202 cohort participants, neither overall vaccine virus fecal shedding nor seroconversion differed by HBGA phenotype. In 238 case-control infants, nonsecretor phenotype was less common in infants with clinical vaccine failure (odds ratio [OR], 0.39; 95% confidence interval [CI], 0.20–0.75). Nonsecretor phenotype was less common in infants with P[8] RVGE (OR, 0.12; 95% CI, 0.03–0.50) and P[4] RVGE (OR, 0.17; 95% CI, 0.04–0.75). Lewis-negative phenotype was more common in infants with P[6] RVGE (OR, 3.2; 95% CI, 1.4–7.2). Conclusions Nonsecretor phenotype was associated with reduced risk of rotavirus vaccine failure. There was no significant association between HBGA phenotype and vaccine take. These data refute the hypothesis that high prevalence of nonsecretor/Lewis-negative phenotypes contributes to lower rotavirus vaccine effectiveness in Malawi