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/

UNCLES: Method for the identification of genes differentially consistently co-expressed in a specific subset of datasets

Background: Collective analysis of the increasingly emerging gene expression datasets are required. The recently proposed binarisation of consensus partition matrices (Bi-CoPaM) method can combine clustering results from multiple datasets to identify the subsets of genes which are consistently co-expressed in all of the provided datasets in a tuneable manner. However, results validation and parameter setting are issues that complicate the design of such methods. Moreover, although it is a common practice to test methods by application to synthetic datasets, the mathematical models used to synthesise such datasets are usually based on approximations which may not always be sufficiently representative of real datasets. Results: Here, we propose an unsupervised method for the unification of clustering results from multiple datasets using external specifications (UNCLES). This method has the ability to identify the subsets of genes consistently co-expressed in a subset of datasets while being poorly co-expressed in another subset of datasets, and to identify the subsets of genes consistently co-expressed in all given datasets. We also propose the M-N scatter plots validation technique and adopt it to set the parameters of UNCLES, such as the number of clusters, automatically. Additionally, we propose an approach for the synthesis of gene expression datasets using real data profiles in a way which combines the ground-truth-knowledge of synthetic data and the realistic expression values of real data, and therefore overcomes the problem of faithfulness of synthetic expression data modelling. By application to those datasets, we validate UNCLES while comparing it with other conventional clustering methods, and of particular relevance, biclustering methods. We further validate UNCLES by application to a set of 14 real genome-wide yeast datasets as it produces focused clusters that conform well to known biological facts. Furthermore, in-silico-based hypotheses regarding the function of a few previously unknown genes in those focused clusters are drawn. Conclusions: The UNCLES method, the M-N scatter plots technique, and the expression data synthesis approach will have wide application for the comprehensive analysis of genomic and other sources of multiple complex biological datasets. Moreover, the derived in-silico-based biological hypotheses represent subjects for future functional studies.The National Institute for Health Research (NIHR) under its Programme Grants for Applied Research Programme (Grant Reference Number RP-PG-0310-1004)

Magnetic Control of Valley Pseudospin in Monolayer WSe2

Local energy extrema of the bands in momentum space, or valleys, can endow electrons in solids with pseudo-spin in addition to real spin. In transition metal dichalcogenides this valley pseudo-spin, like real spin, is associated with a magnetic moment which underlies the valley-dependent circular dichroism that allows optical generation of valley polarization, intervalley quantum coherence, and the valley Hall effect. However, magnetic manipulation of valley pseudospin via this magnetic moment, analogous to what is possible with real spin, has not been shown before. Here we report observation of the valley Zeeman splitting and magnetic tuning of polarization and coherence of the excitonic valley pseudospin, by performing polarization-resolved magneto-photoluminescence on monolayer WSe2. Our measurements reveal both the atomic orbital and lattice contributions to the valley orbital magnetic moment; demonstrate the deviation of the band edges in the valleys from an exact massive Dirac fermion model; and reveal a striking difference between the magnetic responses of neutral and charged valley excitons which is explained by renormalization of the excitonic spectrum due to strong exchange interactions

Pharmacomechanical Catheter-Directed Thrombolysis for Deep-Vein Thrombosis.

BACKGROUND: The post-thrombotic syndrome frequently develops in patients with proximal deep-vein thrombosis despite treatment with anticoagulant therapy. Pharmacomechanical catheter-directed thrombolysis (hereafter pharmacomechanical thrombolysis ) rapidly removes thrombus and is hypothesized to reduce the risk of the post-thrombotic syndrome. METHODS: We randomly assigned 692 patients with acute proximal deep-vein thrombosis to receive either anticoagulation alone (control group) or anticoagulation plus pharmacomechanical thrombolysis (catheter-mediated or device-mediated intrathrombus delivery of recombinant tissue plasminogen activator and thrombus aspiration or maceration, with or without stenting). The primary outcome was development of the post-thrombotic syndrome between 6 and 24 months of follow-up. RESULTS: Between 6 and 24 months, there was no significant between-group difference in the percentage of patients with the post-thrombotic syndrome (47% in the pharmacomechanical-thrombolysis group and 48% in the control group; risk ratio, 0.96; 95% confidence interval [CI], 0.82 to 1.11; P=0.56). Pharmacomechanical thrombolysis led to more major bleeding events within 10 days (1.7% vs. 0.3% of patients, P=0.049), but no significant difference in recurrent venous thromboembolism was seen over the 24-month follow-up period (12% in the pharmacomechanical-thrombolysis group and 8% in the control group, P=0.09). Moderate-to-severe post-thrombotic syndrome occurred in 18% of patients in the pharmacomechanical-thrombolysis group versus 24% of those in the control group (risk ratio, 0.73; 95% CI, 0.54 to 0.98; P=0.04). Severity scores for the post-thrombotic syndrome were lower in the pharmacomechanical-thrombolysis group than in the control group at 6, 12, 18, and 24 months of follow-up (P CONCLUSIONS: Among patients with acute proximal deep-vein thrombosis, the addition of pharmacomechanical catheter-directed thrombolysis to anticoagulation did not result in a lower risk of the post-thrombotic syndrome but did result in a higher risk of major bleeding. (Funded by the National Heart, Lung, and Blood Institute and others; ATTRACT ClinicalTrials.gov number, NCT00790335 .)

Predicting Cellular Growth from Gene Expression Signatures

Maintaining balanced growth in a changing environment is a fundamental systems-level challenge for cellular physiology, particularly in microorganisms. While the complete set of regulatory and functional pathways supporting growth and cellular proliferation are not yet known, portions of them are well understood. In particular, cellular proliferation is governed by mechanisms that are highly conserved from unicellular to multicellular organisms, and the disruption of these processes in metazoans is a major factor in the development of cancer. In this paper, we develop statistical methodology to identify quantitative aspects of the regulatory mechanisms underlying cellular proliferation in Saccharomyces cerevisiae. We find that the expression levels of a small set of genes can be exploited to predict the instantaneous growth rate of any cellular culture with high accuracy. The predictions obtained in this fashion are robust to changing biological conditions, experimental methods, and technological platforms. The proposed model is also effective in predicting growth rates for the related yeast Saccharomyces bayanus and the highly diverged yeast Schizosaccharomyces pombe, suggesting that the underlying regulatory signature is conserved across a wide range of unicellular evolution. We investigate the biological significance of the gene expression signature that the predictions are based upon from multiple perspectives: by perturbing the regulatory network through the Ras/PKA pathway, observing strong upregulation of growth rate even in the absence of appropriate nutrients, and discovering putative transcription factor binding sites, observing enrichment in growth-correlated genes. More broadly, the proposed methodology enables biological insights about growth at an instantaneous time scale, inaccessible by direct experimental methods. Data and tools enabling others to apply our methods are available at http://function.princeton.edu/growthrate

Gravitational Waves From Known Pulsars: Results From The Initial Detector Era

We present the results of searches for gravitational waves from a large selection of pulsars using data from the most recent science runs (S6, VSR2 and VSR4) of the initial generation of interferometric gravitational wave detectors LIGO (Laser Interferometric Gravitational-wave Observatory) and Virgo. We do not see evidence for gravitational wave emission from any of the targeted sources but produce upper limits on the emission amplitude. We highlight the results from seven young pulsars with large spin-down luminosities. We reach within a factor of five of the canonical spin-down limit for all seven of these, whilst for the Crab and Vela pulsars we further surpass their spin-down limits. We present new or updated limits for 172 other pulsars (including both young and millisecond pulsars). Now that the detectors are undergoing major upgrades, and, for completeness, we bring together all of the most up-to-date results from all pulsars searched for during the operations of the first-generation LIGO, Virgo and GEO600 detectors. This gives a total of 195 pulsars including the most recent results described in this paper.United States National Science FoundationScience and Technology Facilities Council of the United KingdomMax-Planck-SocietyState of Niedersachsen/GermanyAustralian Research CouncilInternational Science Linkages program of the Commonwealth of AustraliaCouncil of Scientific and Industrial Research of IndiaIstituto Nazionale di Fisica Nucleare of ItalySpanish Ministerio de Economia y CompetitividadConselleria d'Economia Hisenda i Innovacio of the Govern de les Illes BalearsNetherlands Organisation for Scientific ResearchPolish Ministry of Science and Higher EducationFOCUS Programme of Foundation for Polish ScienceRoyal SocietyScottish Funding CouncilScottish Universities Physics AllianceNational Aeronautics and Space AdministrationOTKA of HungaryLyon Institute of Origins (LIO)National Research Foundation of KoreaIndustry CanadaProvince of Ontario through the Ministry of Economic Development and InnovationNational Science and Engineering Research Council CanadaCarnegie TrustLeverhulme TrustDavid and Lucile Packard FoundationResearch CorporationAlfred P. Sloan FoundationAstronom

Conserved molecular interactions in centriole-to-centrosome conversion.

Centrioles are required to assemble centrosomes for cell division and cilia for motility and signalling. New centrioles assemble perpendicularly to pre-existing ones in G1-S and elongate throughout S and G2. Fully elongated daughter centrioles are converted into centrosomes during mitosis to be able to duplicate and organize pericentriolar material in the next cell cycle. Here we show that centriole-to-centrosome conversion requires sequential loading of Cep135, Ana1 (Cep295) and Asterless (Cep152) onto daughter centrioles during mitotic progression in both Drosophila melanogaster and human. This generates a molecular network spanning from the inner- to outermost parts of the centriole. Ana1 forms a molecular strut within the network, and its essential role can be substituted by an engineered fragment providing an alternative linkage between Asterless and Cep135. This conserved architectural framework is essential for loading Asterless or Cep152, the partner of the master regulator of centriole duplication, Plk4. Our study thus uncovers the molecular basis for centriole-to-centrosome conversion that renders daughter centrioles competent for motherhood.J.F., Z.L., S.S. and N.S.D. are supported from Programme Grant to D.M.G. from Cancer Research UK. H.R. is supported from MRC Programme Grant to D.M.G. J.F. thank the British Academy and the Royal Society for Newton International Fellowship and Z.L. thanks the Federation of European Biochemical Societies for the Long-Term postdoctoral Fellowship. The authors thank Nicola Lawrence and Alex Sossick for assistance with 3D-SIM.This is the author accepted manuscript. The final version is available from NPG via http://dx.doi.org/10.1038/ncb327

Measurement of the cross-section and charge asymmetry of WW bosons produced in proton-proton collisions at s=8\sqrt{s}=8 TeV with the ATLAS detector

This paper presents measurements of the $W^+ \rightarrow \mu^+\nu$ and $W^- \rightarrow \mu^-\nu$ cross-sections and the associated charge asymmetry as a function of the absolute pseudorapidity of the decay muon. The data were collected in proton--proton collisions at a centre-of-mass energy of 8 TeV with the ATLAS experiment at the LHC and correspond to a total integrated luminosity of 20.2~\mbox{fb^{-1}}. The precision of the cross-section measurements varies between 0.8% to 1.5% as a function of the pseudorapidity, excluding the 1.9% uncertainty on the integrated luminosity. The charge asymmetry is measured with an uncertainty between 0.002 and 0.003. The results are compared with predictions based on next-to-next-to-leading-order calculations with various parton distribution functions and have the sensitivity to discriminate between them.Comment: 38 pages in total, author list starting page 22, 5 figures, 4 tables, submitted to EPJC. All figures including auxiliary figures are available at https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/STDM-2017-13

Search for chargino-neutralino production with mass splittings near the electroweak scale in three-lepton final states in √s=13 TeV pp collisions with the ATLAS detector

A search for supersymmetry through the pair production of electroweakinos with mass splittings near the electroweak scale and decaying via on-shell W and Z bosons is presented for a three-lepton final state. The analyzed proton-proton collision data taken at a center-of-mass energy of √s=13  TeV were collected between 2015 and 2018 by the ATLAS experiment at the Large Hadron Collider, corresponding to an integrated luminosity of 139  fb−1. A search, emulating the recursive jigsaw reconstruction technique with easily reproducible laboratory-frame variables, is performed. The two excesses observed in the 2015–2016 data recursive jigsaw analysis in the low-mass three-lepton phase space are reproduced. Results with the full data set are in agreement with the Standard Model expectations. They are interpreted to set exclusion limits at the 95% confidence level on simplified models of chargino-neutralino pair production for masses up to 345 GeV