Reconciliation Revisited: Handling Multiple Optima when Reconciling with Duplication, Transfer, and Loss

Phylogenetic tree reconciliation is a powerful approach for inferring evolutionary events like gene duplication, horizontal gene transfer, and gene loss, which are fundamental to our understanding of molecular evolution. While duplication–loss (DL) reconciliation leads to a unique maximum-parsimony solution, duplication-transfer-loss (DTL) reconciliation yields a multitude of optimal solutions, making it difficult to infer the true evolutionary history of the gene family. This problem is further exacerbated by the fact that different event cost assignments yield different sets of optimal reconciliations. Here, we present an effective, efficient, and scalable method for dealing with these fundamental problems in DTL reconciliation. Our approach works by sampling the space of optimal reconciliations uniformly at random and aggregating the results. We show that even gene trees with only a few dozen genes often have millions of optimal reconciliations and present an algorithm to efficiently sample the space of optimal reconciliations uniformly at random in O(mn[superscript 2]) time per sample, where m and n denote the number of genes and species, respectively. We use these samples to understand how different optimal reconciliations vary in their node mappings and event assignments and to investigate the impact of varying event costs. We apply our method to a biological dataset of approximately 4700 gene trees from 100 taxa and observe that 93% of event assignments and 73% of mappings remain consistent across different multiple optima. Our analysis represents the first systematic investigation of the space of optimal DTL reconciliations and has many important implications for the study of gene family evolution.National Science Foundation (U.S.) (CAREER Award 0644282)National Institutes of Health (U.S.) (Grant RC2 HG005639)National Science Foundation (U.S.). Assembling the Tree of Life (Program) (Grant 0936234

J-PLUS: The javalambre photometric local universe survey

ABSTRACT: TheJavalambrePhotometric Local UniverseSurvey (J-PLUS )isanongoing 12-band photometricopticalsurvey, observingthousands of squaredegrees of theNorthernHemispherefromthededicated JAST/T80 telescope at the Observatorio Astrofísico de Javalambre (OAJ). The T80Cam is a camera with a field of view of 2 deg2 mountedon a telescopewith a diameter of 83 cm, and isequippedwith a uniquesystem of filtersspanningtheentireopticalrange (3500–10 000 Å). Thisfiltersystemis a combination of broad-, medium-, and narrow-band filters, optimallydesigned to extracttherest-framespectralfeatures (the 3700–4000 Å Balmer break region, Hδ, Ca H+K, the G band, and the Mg b and Ca triplets) that are key to characterizingstellartypes and delivering a low-resolutionphotospectrumforeach pixel of theobservedsky. With a typicaldepth of AB ∼21.25 mag per band, thisfilter set thusallowsforanunbiased and accuratecharacterization of thestellarpopulation in our Galaxy, itprovidesanunprecedented 2D photospectralinformationforall resolved galaxies in the local Universe, as well as accuratephoto-z estimates (at the δ z/(1 + z)∼0.005–0.03 precisionlevel) formoderatelybright (up to r ∼ 20 mag) extragalacticsources. Whilesomenarrow-band filters are designedforthestudy of particular emissionfeatures ([O II]/λ3727, Hα/λ6563) up to z < 0.017, theyalsoprovidewell-definedwindowsfortheanalysis of otheremissionlines at higherredshifts. As a result, J-PLUS has thepotential to contribute to a widerange of fields in Astrophysics, both in thenearbyUniverse (MilkyWaystructure, globular clusters, 2D IFU-likestudies, stellarpopulations of nearby and moderate-redshiftgalaxies, clusters of galaxies) and at highredshifts (emission-line galaxies at z ≈ 0.77, 2.2, and 4.4, quasi-stellarobjects, etc.). Withthispaper, wereleasethefirst∼1000 deg2 of J-PLUS data, containingabout 4.3 millionstars and 3.0 milliongalaxies at r <  21mag. With a goal of 8500 deg2 forthe total J-PLUS footprint, thesenumbers are expected to rise to about 35 millionstars and 24 milliongalaxiesbytheend of thesurvey.Funding for the J-PLUS Project has been provided by the Governments of Spain and Aragón through the Fondo de Inversiones de Teruel, the Spanish Ministry of Economy and Competitiveness (MINECO; under grants AYA2017-86274-P, AYA2016-77846-P, AYA2016-77237-C3-1-P, AYA2015-66211-C2-1-P, AYA2015-66211-C2-2, AYA2012-30789, AGAUR grant SGR-661/2017, and ICTS-2009-14), and European FEDER funding (FCDD10-4E-867, FCDD13-4E-2685

Study of the lineshape of the chi(c1) (3872) state

A study of the lineshape of the chi(c1) (3872) state is made using a data sample corresponding to an integrated luminosity of 3 fb(-1) collected in pp collisions at center-of-mass energies of 7 and 8 TeV with the LHCb detector. Candidate chi(c1)(3872) and psi(2S) mesons from b-hadron decays are selected in the J/psi pi(+)pi(-) decay mode. Describing the lineshape with a Breit-Wigner function, the mass splitting between the chi(c1 )(3872) and psi(2S) states, Delta m, and the width of the chi(c1 )(3872) state, Gamma(Bw), are determined to be (Delta m=185.598 +/- 0.067 +/- 0.068 Mev,)(Gamma BW=1.39 +/- 0.24 +/- 0.10 Mev,) where the first uncertainty is statistical and the second systematic. Using a Flatte-inspired model, the mode and full width at half maximum of the lineshape are determined to be (mode=3871.69+0.00+0.05 MeV.)(FWHM=0.22-0.04+0.13+0.07+0.11-0.06-0.13 MeV, ) An investigation of the analytic structure of the Flatte amplitude reveals a pole structure, which is compatible with a quasibound D-0(D) over bar*(0) state but a quasivirtual state is still allowed at the level of 2 standard deviations

Measurement of the CKM angle γγ in B±→DK±B^\pm\to D K^\pm and B±→Dπ±B^\pm \to D π^\pm decays with D→KS0h+h−D \to K_\mathrm S^0 h^+ h^-

A measurement of $CP$-violating observables is performed using the decays $B^\pm\to D K^\pm$ and $B^\pm\to D \pi^\pm$, where the $D$ meson is reconstructed in one of the self-conjugate three-body final states $K_{\mathrm S}\pi^+\pi^-$ and $K_{\mathrm S}K^+K^-$ (commonly denoted $K_{\mathrm S} h^+h^-$). The decays are analysed in bins of the $D$-decay phase space, leading to a measurement that is independent of the modelling of the $D$-decay amplitude. The observables are interpreted in terms of the CKM angle $\gamma$. Using a data sample corresponding to an integrated luminosity of $9\,\text{fb}^{-1}$ collected in proton-proton collisions at centre-of-mass energies of $7$, $8$, and $13\,\text{TeV}$ with the LHCb experiment, $\gamma$ is measured to be $\left(68.7^{+5.2}_{-5.1}\right)^\circ$. The hadronic parameters $r_B^{DK}$, $r_B^{D\pi}$, $\delta_B^{DK}$, and $\delta_B^{D\pi}$, which are the ratios and strong-phase differences of the suppressed and favoured $B^\pm$ decays, are also reported

Measurement of forward charged hadron flow harmonics in peripheral PbPb collisions at √sNN = 5.02 TeV with the LHCb detector

Flow harmonic coefficients, v n , which are the key to studying the hydrodynamics of the quark-gluon plasma (QGP) created in heavy-ion collisions, have been measured in various collision systems and kinematic regions and using various particle species. The study of flow harmonics in a wide pseudorapidity range is particularly valuable to understand the temperature dependence of the shear viscosity to entropy density ratio of the QGP. This paper presents the first LHCb results of the second- and the third-order flow harmonic coefficients of charged hadrons as a function of transverse momentum in the forward region, corresponding to pseudorapidities between 2.0 and 4.9, using the data collected from PbPb collisions in 2018 at a center-of-mass energy of 5.02 TeV . The coefficients measured using the two-particle angular correlation analysis method are smaller than the central-pseudorapidity measurements at ALICE and ATLAS from the same collision system but share similar features

Helium identification with LHCb

The identification of helium nuclei at LHCb is achieved using a method based on measurements of ionisation losses in the silicon sensors and timing measurements in the Outer Tracker drift tubes. The background from photon conversions is reduced using the RICH detectors and an isolation requirement. The method is developed using pp collision data at √(s) = 13 TeV recorded by the LHCb experiment in the years 2016 to 2018, corresponding to an integrated luminosity of 5.5 fb-1. A total of around 105 helium and antihelium candidates are identified with negligible background contamination. The helium identification efficiency is estimated to be approximately 50% with a corresponding background rejection rate of up to O(10^12). These results demonstrate the feasibility of a rich programme of measurements of QCD and astrophysics interest involving light nuclei

Curvature-bias corrections using a pseudomass method

Momentum measurements for very high momentum charged particles, such as muons from electroweak vector boson decays, are particularly susceptible to charge-dependent curvature biases that arise from misalignments of tracking detectors. Low momentum charged particles used in alignment procedures have limited sensitivity to coherent displacements of such detectors, and therefore are unable to fully constrain these misalignments to the precision necessary for studies of electroweak physics. Additional approaches are therefore required to understand and correct for these effects. In this paper the curvature biases present at the LHCb detector are studied using the pseudomass method in proton-proton collision data recorded at centre of mass energy √(s)=13 TeV during 2016, 2017 and 2018. The biases are determined using Z→μ + μ - decays in intervals defined by the data-taking period, magnet polarity and muon direction. Correcting for these biases, which are typically at the 10-4 GeV-1 level, improves the Z→μ + μ - mass resolution by roughly 18% and eliminates several pathological trends in the kinematic-dependence of the mean dimuon invariant mass

Study of CP violation in B0 → DK⋆(892)0 decays with D → Kπ(ππ), ππ(ππ), and KK final states

A measurement of CP-violating observables associated with the interference of B0 → D0K⋆ (892)0 and B0 → D¯ 0K⋆ (892)0 decay amplitudes is performed in the D0 → K∓π ±(π +π −), D0 → π +π −(π +π −), and D0 → K+K− fnal states using data collected by the LHCb experiment corresponding to an integrated luminosity of 9 fb−1 . CP-violating observables related to the interference of B0 s → D0K¯ ⋆ (892)0 and B0 s → D¯ 0K¯ ⋆ (892)0 are also measured, but no evidence for interference is found. The B0 observables are used to constrain the parameter space of the CKM angle γ and the hadronic parameters r DK⋆ B0 and δ DK⋆ B0 with inputs from other measurements. In a combined analysis, these measurements allow for four solutions in the parameter space, only one of which is consistent with the world average

Study of the doubly charmed tetraquark T+cc

Quantum chromodynamics, the theory of the strong force, describes interactions of coloured quarks and gluons and the formation of hadronic matter. Conventional hadronic matter consists of baryons and mesons made of three quarks and quark-antiquark pairs, respectively. Particles with an alternative quark content are known as exotic states. Here a study is reported of an exotic narrow state in the D0D0π+ mass spectrum just below the D*+D0 mass threshold produced in proton-proton collisions collected with the LHCb detector at the Large Hadron Collider. The state is consistent with the ground isoscalar T+cc tetraquark with a quark content of ccu⎯⎯⎯d⎯⎯⎯ and spin-parity quantum numbers JP = 1+. Study of the DD mass spectra disfavours interpretation of the resonance as the isovector state. The decay structure via intermediate off-shell D*+ mesons is consistent with the observed D0π+ mass distribution. To analyse the mass of the resonance and its coupling to the D*D system, a dedicated model is developed under the assumption of an isoscalar axial-vector T+cc state decaying to the D*D channel. Using this model, resonance parameters including the pole position, scattering length, effective range and compositeness are determined to reveal important information about the nature of the T+cc state. In addition, an unexpected dependence of the production rate on track multiplicity is observed

The 2006 NESCent phyloinformatics hackathon: A field report

10.1177/117693430700300016Evolutionary Bioinformatics3287-29