Algorithms in nature: the convergence of systems biology and computational thinking

Biologists rely on computational methods to analyze and integrate large data sets, while several computational methods were inspired by the high-level design principles of biological systems. This Perspectives discusses the recent convergence of these two ways of thinking

Canalization of Gene Expression and Domain Shifts in the Drosophila Blastoderm by Dynamical Attractors

The variation in the expression patterns of the gap genes in the blastoderm of the fruit fly Drosophila melanogaster reduces over time as a result of cross regulation between these genes, a fact that we have demonstrated in an accompanying article in PLoS Biology (see Manu et al., doi:10.1371/journal.pbio.1000049). This biologically essential process is an example of the phenomenon known as canalization. It has been suggested that the developmental trajectory of a wild-type organism is inherently stable, and that canalization is a manifestation of this property. Although the role of gap genes in the canalization process was established by correctly predicting the response of the system to particular perturbations, the stability of the developmental trajectory remains to be investigated. For many years, it has been speculated that stability against perturbations during development can be described by dynamical systems having attracting sets that drive reductions of volume in phase space. In this paper, we show that both the reduction in variability of gap gene expression as well as shifts in the position of posterior gap gene domains are the result of the actions of attractors in the gap gene dynamical system. Two biologically distinct dynamical regions exist in the early embryo, separated by a bifurcation at 53% egg length. In the anterior region, reduction in variation occurs because of stability induced by point attractors, while in the posterior, the stability of the developmental trajectory arises from a one-dimensional attracting manifold. This manifold also controls a previously characterized anterior shift of posterior region gap domains. Our analysis shows that the complex phenomena of canalization and pattern formation in the Drosophila blastoderm can be understood in terms of the qualitative features of the dynamical system. The result confirms the idea that attractors are important for developmental stability and shows a richer variety of dynamical attractors in developmental systems than has been previously recognized

Control of Canalization and Evolvability by Hsp90

Partial reduction of Hsp90 increases expression of morphological novelty in qualitative traits of Drosophila and Arabidopsis, but the extent to which the Hsp90 chaperone also controls smaller and more likely adaptive changes in natural quantitative traits has been unclear. To determine the effect of Hsp90 on quantitative trait variability we deconstructed genetic, stochastic and environmental components of variation in Drosophila wing and bristle traits of genetically matched flies, differing only by Hsp90 loss-of-function or wild-type alleles. Unexpectedly, Hsp90 buffering was remarkably specific to certain normally invariant and highly discrete quantitative traits. Like the qualitative trait phenotypes controlled by Hsp90, highly discrete quantitative traits such as scutellor and thoracic bristle number are threshold traits. When tested across genotypes sampled from a wild population or in laboratory strains, the sensitivity of these traits to many types of variation was coordinately controlled, while continuously variable bristle types and wing size, and critically invariant left-right wing asymmetry, remained relatively unaffected. Although increased environmental variation and developmental noise would impede many types of selection response, in replicate populations in which Hsp90 was specifically impaired, heritability and ‘extrinsic evolvability’, the expected response to selection, were also markedly increased. However, despite the overall buffering effect of Hsp90 on variation in populations, for any particular individual or genotype in which Hsp90 was impaired, the size and direction of its effects were unpredictable. The trait and genetic-background dependence of Hsp90 effects and its remarkable bias toward invariant or canalized traits support the idea that traits evolve independent and trait-specific mechanisms of canalization and evolvability through their evolution of non-linearity and thresholds. Highly non-linear responses would buffer variation in Hsp90-dependent signaling over a wide range, while over a narrow range of signaling near trait thresholds become more variable with increasing probability of triggering all-or-none developmental responses

Studies of new Higgs boson interactions through nonresonant HH production in the b¯bγγ fnal state in pp collisions at √s = 13 TeV with the ATLAS detector

A search for nonresonant Higgs boson pair production in the b ¯bγγ fnal state is performed using 140 fb−1 of proton-proton collisions at a centre-of-mass energy of 13 TeV recorded by the ATLAS detector at the CERN Large Hadron Collider. This analysis supersedes and expands upon the previous nonresonant ATLAS results in this fnal state based on the same data sample. The analysis strategy is optimised to probe anomalous values not only of the Higgs (H) boson self-coupling modifer κλ but also of the quartic HHV V (V = W, Z) coupling modifer κ2V . No signifcant excess above the expected background from Standard Model processes is observed. An observed upper limit µHH < 4.0 is set at 95% confdence level on the Higgs boson pair production cross-section normalised to its Standard Model prediction. The 95% confdence intervals for the coupling modifers are −1.4 < κλ < 6.9 and −0.5 < κ2V < 2.7, assuming all other Higgs boson couplings except the one under study are fxed to the Standard Model predictions. The results are interpreted in the Standard Model efective feld theory and Higgs efective feld theory frameworks in terms of constraints on the couplings of anomalous Higgs boson (self-)interactions

Electron and photon energy calibration with the ATLAS detector using LHC Run 2 data

This paper presents the electron and photon energy calibration obtained with the ATLAS detector using 140 fb−1 of LHC proton-proton collision data recorded at root(s) = 13 TeV between 2015 and 2018. Methods for the measurement of electron and photon energies are outlined, along with the current knowledge of the passive material in front of the ATLAS electromagnetic calorimeter. The energy calibration steps are discussed in detail, with emphasis on the improvements introduced in this paper. The absolute energy scale is set using a large sample of Z-boson decays into electron-positron pairs, and its residual dependence on the electron energy is used for the first time to further constrain systematic uncertainties. The achieved calibration uncertainties are typically 0.05% for electrons from resonant Z-boson decays, 0.4% at ET tilde 10 GeV, and 0.3% at ET tilde 1 TeV; for photons at ET tilde 60 GeV, they are 0.2% on average. This is more than twice as precise as the previous calibration. The new energy calibration is validated using J/psi -> ee and radiative Z-boson decays

Measurement of exclusive pion pair production in proton–proton collisions at √s=7 TeV with the ATLAS detector

The exclusive production of pion pairs in the process pp→ ppπ+π- has been measured at s=7TeV with the ATLAS detector at the LHC, using 80μb-1 of low-luminosity data. The pion pairs were detected in the ATLAS central detector while outgoing protons were measured in the forward ATLAS ALFA detector system. This represents the first use of proton tagging to measure an exclusive hadronic final state at the LHC. A cross-section measurement is performed in two kinematic regions defined by the proton momenta, the pion rapidities and transverse momenta, and the pion–pion invariant mass. Cross-section values of 4.8±1.0(stat)-0.2+0.3(syst)μb and 9±6(stat)-2+2(syst)μb are obtained in the two regions; they are compared with theoretical models and provide a demonstration of the feasibility of measurements of this type

A search for the dimuon decay of the Standard Model Higgs boson with the ATLAS detector

A search for the dimuon decay of the Standard Model (SM) Higgs boson is performed using data corresponding to an integrated luminosity of 139 fb(-1) collected with the ATLAS detector in Run 2 pp collisions at root s = 13 TeV at the Large Hadron Collider. The observed (expected) significance over the background-only hypothesis for a Higgs boson with a mass of 125.09 GeV is 2.0 sigma (1.7 sigma). The observed upper limit on the cross section times branching ratio for pp -> H -> mu mu is 2.2 times the SM prediction at 95% confidence level, while the expected limit on a H -> mu mu signal assuming the absence (presence) of a SM signal is 1.1(2.0). The best-fit value of the signal strength parameter, defined as the ratio of the observed signal yield to the one expected in the SM, is mu = 1.2 +/- 0.6. (C) 2020 The Author(s). Published by Elsevier B.V

Measurement of the H → γ γ and H → ZZ∗ → 4 cross-sections in pp collisions at √s = 13.6 TeV with the ATLAS detector

The inclusive Higgs boson production cross section is measured in the di-photon and the Z Z∗ → 4 decay channels using 31.4 and 29.0 fb−1 of pp collision data respectively, collected with the ATLAS detector at a centre of-mass energy of √s = 13.6 TeV. To reduce the model dependence, the measurement in each channel is restricted to a particle-level phase space that closely matches the chan nel’s detector-level kinematic selection, and it is corrected for detector effects. These measured fiducial cross-sections are σfid,γ γ = 76+14 −13 fb, and σfid,4 = 2.80 ± 0.74 fb, in agreement with the corresponding Standard Model predic tions of 67.6±3.7 fb and 3.67±0.19 fb. Assuming Standard Model acceptances and branching fractions for the two chan nels, the fiducial measurements are extrapolated to the full phase space yielding total cross-sections of σ (pp → H) = 67+12 −11 pb and 46±12 pb at 13.6 TeV from the di-photon and Z Z∗ → 4 measurements respectively. The two measure ments are combined into a total cross-section measurement of σ (pp → H) = 58.2±8.7 pb, to be compared with the Stan dard Model prediction of σ (pp → H)SM = 59.9 ± 2.6 p