Misty Mountain clustering: application to fast unsupervised flow cytometry gating

<p>Abstract</p> <p>Background</p> <p>There are many important clustering questions in computational biology for which no satisfactory method exists. Automated clustering algorithms, when applied to large, multidimensional datasets, such as flow cytometry data, prove unsatisfactory in terms of speed, problems with local minima or cluster shape bias. Model-based approaches are restricted by the assumptions of the fitting functions. Furthermore, model based clustering requires serial clustering for all cluster numbers within a user defined interval. The final cluster number is then selected by various criteria. These supervised serial clustering methods are time consuming and frequently different criteria result in different optimal cluster numbers. Various unsupervised heuristic approaches that have been developed such as affinity propagation are too expensive to be applied to datasets on the order of 10⁶points that are often generated by high throughput experiments.</p> <p>Results</p> <p>To circumvent these limitations, we developed a new, unsupervised density contour clustering algorithm, called Misty Mountain, that is based on percolation theory and that efficiently analyzes large data sets. The approach can be envisioned as a progressive top-down removal of clouds covering a data histogram relief map to identify clusters by the appearance of statistically distinct peaks and ridges. This is a parallel clustering method that finds every cluster after analyzing only once the cross sections of the histogram. The overall run time for the composite steps of the algorithm increases linearly by the number of data points. The clustering of 10⁶data points in 2D data space takes place within about 15 seconds on a standard laptop PC. Comparison of the performance of this algorithm with other state of the art automated flow cytometry gating methods indicate that Misty Mountain provides substantial improvements in both run time and in the accuracy of cluster assignment.</p> <p>Conclusions</p> <p>Misty Mountain is fast, unbiased for cluster shape, identifies stable clusters and is robust to noise. It provides a useful, general solution for multidimensional clustering problems. We demonstrate its suitability for automated gating of flow cytometry data.</p

Search for new phenomena in final states with an energetic jet and large missing transverse momentum in pp collisions at √ s = 8 TeV with the ATLAS detector

Results of a search for new phenomena in final states with an energetic jet and large missing transverse momentum are reported. The search uses 20.3 fb−1 of √ s = 8 TeV data collected in 2012 with the ATLAS detector at the LHC. Events are required to have at least one jet with pT > 120 GeV and no leptons. Nine signal regions are considered with increasing missing transverse momentum requirements between Emiss T > 150 GeV and Emiss T > 700 GeV. Good agreement is observed between the number of events in data and Standard Model expectations. The results are translated into exclusion limits on models with either large extra spatial dimensions, pair production of weakly interacting dark matter candidates, or production of very light gravitinos in a gauge-mediated supersymmetric model. In addition, limits on the production of an invisibly decaying Higgs-like boson leading to similar topologies in the final state are presente

The CAFA challenge reports improved protein function prediction and new functional annotations for hundreds of genes through experimental screens

Background The Critical Assessment of Functional Annotation (CAFA) is an ongoing, global, community-driven effort to evaluate and improve the computational annotation of protein function. Results Here, we report on the results of the third CAFA challenge, CAFA3, that featured an expanded analysis over the previous CAFA rounds, both in terms of volume of data analyzed and the types of analysis performed. In a novel and major new development, computational predictions and assessment goals drove some of the experimental assays, resulting in new functional annotations for more than 1000 genes. Specifically, we performed experimental whole-genome mutation screening in Candida albicans and Pseudomonas aureginosa genomes, which provided us with genome-wide experimental data for genes associated with biofilm formation and motility. We further performed targeted assays on selected genes in Drosophila melanogaster, which we suspected of being involved in long-term memory. Conclusion We conclude that while predictions of the molecular function and biological process annotations have slightly improved over time, those of the cellular component have not. Term-centric prediction of experimental annotations remains equally challenging; although the performance of the top methods is significantly better than the expectations set by baseline methods in C. albicans and D. melanogaster, it leaves considerable room and need for improvement. Finally, we report that the CAFA community now involves a broad range of participants with expertise in bioinformatics, biological experimentation, biocuration, and bio-ontologies, working together to improve functional annotation, computational function prediction, and our ability to manage big data in the era of large experimental screens.Peer reviewe

Sexual selection predicts the rate and direction of colour divergence in a large avian radiation

Sexual selection is proposed to be a powerful driver of phenotypic evolution in animal systems. At macroevolutionary scales, sexual selection can theoretically drive both the rate and direction of phenotypic evolution, but this hypothesis remains contentious. Here, we find that differences in the rate and direction of plumage colour evolution are predicted by a proxy for sexual selection intensity (plumage dichromatism) in a large radiation of suboscine passerine birds (Tyrannida). We show that rates of plumage evolution are correlated between the sexes, but that sexual selection has a strong positive effect on male, but not female, interspecific divergence rates. Furthermore, we demonstrate that rapid male plumage divergence is biased towards carotenoid-based (red/yellow) colours widely assumed to represent honest sexual signals. Our results highlight the central role of sexual selection in driving avian colour divergence, and reveal the existence of convergent evolutionary responses of animal signalling traits under sexual selection

Combined measurement of differential and total cross sections in the H → γγ and the H → ZZ* → 4ℓ decay channels at s=13 TeV with the ATLAS detector

A combined measurement of differential and inclusive total cross sections of Higgs boson production is performed using 36.1 fb−1 of 13 TeV proton–proton collision data produced by the LHC and recorded by the ATLAS detector in 2015 and 2016. Cross sections are obtained from measured H→γγ and H→ZZ*(→4ℓ event yields, which are combined taking into account detector efficiencies, resolution, acceptances and branching fractions. The total Higgs boson production cross section is measured to be 57.0−5.9 +6.0 (stat.) −3.3 +4.0 (syst.) pb, in agreement with the Standard Model prediction. Differential cross-section measurements are presented for the Higgs boson transverse momentum distribution, Higgs boson rapidity, number of jets produced together with the Higgs boson, and the transverse momentum of the leading jet. The results from the two decay channels are found to be compatible, and their combination agrees with the Standard Model predictions

Operation and performance of the ATLAS Tile Calorimeter in Run 1

The Tile Calorimeter is the hadron calorimeter covering the central region of the ATLAS experiment at the Large Hadron Collider. Approximately 10,000 photomultipliers collect light from scintillating tiles acting as the active material sandwiched between slabs of steel absorber. This paper gives an overview of the calorimeter’s performance during the years 2008–2012 using cosmic-ray muon events and proton–proton collision data at centre-of-mass energies of 7 and 8TeV with a total integrated luminosity of nearly 30 fb−1. The signal reconstruction methods, calibration systems as well as the detector operation status are presented. The energy and time calibration methods performed excellently, resulting in good stability of the calorimeter response under varying conditions during the LHC Run 1. Finally, the Tile Calorimeter response to isolated muons and hadrons as well as to jets from proton–proton collisions is presented. The results demonstrate excellent performance in accord with specifications mentioned in the Technical Design Report

Performance of missing transverse momentum reconstruction with the ATLAS detector using proton–proton collisions at √s = 13 TeV

The performance of the missing transverse momentum (EmissT) reconstruction with the ATLAS detector is evaluated using data collected in proton–proton collisions at the LHC at a centre-of-mass energy of 13 TeV in 2015. To reconstruct EmissT, fully calibrated electrons, muons, photons, hadronically decaying τ -leptons, and jets reconstructed from calorimeter energy deposits and charged-particle tracks are used. These are combined with the soft hadronic activity measured by reconstructed charged-particle tracks not associated with the hard objects. Possible double counting of contributions from reconstructed charged-particle tracks from the inner detector, energy deposits in the calorimeter, and reconstructed muons from the muon spectrometer is avoided by applying a signal ambiguity resolution procedure which rejects already used signals when combining the various EmissT contributions. The individual terms as well as the overall reconstructed EmissT are evaluated with various performance metrics for scale (linearity), resolution, and sensitivity to the data-taking conditions. The method developed to determine the systematic uncertainties of the EmissT scale and resolution is discussed. Results are shown based on the full 2015 data sample corresponding to an integrated luminosity of 3.2 fb−1

Measurement of VH, H → b b ¯ production as a function of the vector-boson transverse momentum in 13 TeV pp collisions with the ATLAS detector

Cross-sections of associated production of a Higgs boson decaying into bottom-quark pairs and an electroweak gauge boson, W or Z, decaying into leptons are measured as a function of the gauge boson transverse momentum. The measurements are performed in kinematic fiducial volumes defined in the `simplified template cross-section' framework. The results are obtained using 79.8 fb-1 of proton-proton collisions recorded by the ATLAS detector at the Large Hadron Collider at a centre-of-mass energy of 13 TeV. All measurements are found to be in agreement with the Standard Model predictions, and limits are set on the parameters of an effective Lagrangian sensitive to modifications of the Higgs boson couplings to the electroweak gauge bosons

Measurement of the t¯tZ and t¯tW cross sections in proton-proton collisions at √s=13 TeV with the ATLAS detector

A measurement of the associated production of a top-quark pair (t¯t) with a vector boson (W, Z) in proton-proton collisions at a center-of-mass energy of 13 TeV is presented, using 36.1  fb−1 of integrated luminosity collected by the ATLAS detector at the Large Hadron Collider. Events are selected in channels with two same- or opposite-sign leptons (electrons or muons), three leptons or four leptons, and each channel is further divided into multiple regions to maximize the sensitivity of the measurement. The t¯tZ and t¯tW production cross sections are simultaneously measured using a combined fit to all regions. The best-fit values of the production cross sections are σt¯tZ=0.95±0.08stat±0.10syst pb and σt¯tW=0.87±0.13stat±0.14syst pb in agreement with the Standard Model predictions. The measurement of the t¯tZ cross section is used to set constraints on effective field theory operators which modify the t¯tZ vertex