Coalescence Model Research Articles

Museum Skins Enable Identification of Introgression Associated with Cytonuclear Discordance.

Increased sampling of genomes and populations across closely related species has revealed that levels of genetic exchange during and after speciation are higher than previously thought. One obvious manifestation of such exchange is strong cytonuclear discordance, where the divergence in mitochondrial DNA (mtDNA) differs from that for nuclear genes more (or less) than expected from differences between mtDNA and nuclear DNA (nDNA) in population size and mutation rate. Given genome-scale data sets and coalescent modeling, we can now confidently identify cases of strong discordance and test specifically for historical or recent introgression as the cause. Using population sampling, combining exon capture data from historical museum specimens and recently collected tissues we showcase how genomic tools can resolve complex evolutionary histories in the brachyotis group of rock-wallabies (Petrogale). In particular, applying population and phylogenomic approaches we can assess the role of demographic processes in driving complex evolutionary patterns and assess a role of ancient introgression and hybridization. We find that described species are well supported as monophyletic taxa for nDNA genes, but not for mtDNA, with cytonuclear discordance involving at least 4 operational taxonomic units across 4 species which diverged 183-278 kya. ABC modeling of nDNA gene trees supports introgression during or after speciation for some taxon pairs with cytonuclear discordance. Given substantial differences in body size between the species involved, this evidence for gene flow is surprising. Heterogenous patterns of introgression were identified but do not appear to be associated with chromosome differences between species. These and previous results suggest that dynamic past climates across the monsoonal tropics could have promoted reticulation among related species.

Violins we see, species we don’t… Species delimitation of the spider genus Loxosceles Heineken & Lowe (Araneae: Sicariidae) from North America using morphological and molecular evidence

In modern systematics, different sources of evidence are commonly used for the discovery, identification, and delimitation of species, especially when morphology fails to delineate between species or in underestimated species complexes or cryptic species. In this study, morphological data and two DNA barcoding markers—cytochrome c oxidase subunit I (COI) and internal transcribed spacer 2 (ITS2)—were used to delimit species in the spider genus Loxosceles from North America. The molecular species delimitation analyses were carried out using three different methods under the corrected p-distance Neighbor-Joining (NJ) criteria: 1) Assemble Species by Automatic Partitioning (ASAP), 2) General Mixed Yule Coalescent model (GMYC), and 3) Bayesian Poisson Tree Processes (bPTP). The analyses incorporated 192 terminals corresponding to 43 putative species of Loxosceles, of which 15 are newly recognized herein, as putative new species, based on morphology and congruence between molecular methods with COI. The average intraspecific genetic distance (p-distance) was <2%, whereas the average interspecific genetic distance was 15.6%. The GMYC and bPTP molecular methods recovered 65-79 and 69 species respectively, overestimating the diversity in comparison with morphology, whereas the ASAP method delimited 60 species. The morphology of primary sexual structures (males palps and female seminal receptacles) was congruent with most of the molecular methods mainly with COI, showing that they are robust characters for identification at the species level. For species delimitation COI was more informative than ITS2. The diversity of Loxosceles species is still underestimated for North America, particularly in Mexico which holds the highest diversity of this genus worldwide.

Mass splitting and spin alignment for ϕ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\phi $$\\end{document} mesons in a magnetic field in NJL model

Based on the Nambu–Jona–Lasinio (NJL) model, we develop a framework for calculating the spin alignment of vector mesons and applied it to study ϕ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\phi $$\\end{document} mesons in a magnetic field. We calculate mass spectra for ϕ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\phi $$\\end{document} mesons and observe mass splitting between the longitudinally polarized state and transversely polarized states. The ϕ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\phi $$\\end{document} meson in a thermal equilibrium system is preferred to occupy the state with spin λ=0\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\lambda =0$$\\end{document} than those with spin λ=±1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\lambda =\\pm 1$$\\end{document}, because the former state has a smaller energy. As a consequence, we conclude that the spin alignment will be larger than 1/3 if one measures along the direction of the magnetic field, which is qualitatively consistent with the recent STAR data. Around the critical temperature TC=150MeV\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$T_{C}=150~\\hbox {MeV}$$\\end{document}, the positive deviation from 1/3 is proportional to the square of the magnetic field strength, which agrees with the result from the non-relativistic coalescence model. Including the anomalous magnetic moments for quarks will modify the dynamical masses of quarks and thus affect the mass spectra and spin alignment of ϕ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\phi $$\\end{document} mesons. The discussion of spin alignment in the NJL model may help us better understand the formation of hadron’s spin structure during the chiral phase transition.

Estimating between-country migration in pneumococcal populations.

Streptococcus pneumoniae (the pneumococcus) is a globally distributed, human obligate opportunistic bacterial pathogen which, although often carried commensally, is also a significant cause of invasive disease. Apart from multi-drug resistant and virulent clones, the rate and direction of pneumococcal dissemination between different countries remains largely unknown. The ability for the pneumococcus to take a foothold in a country depends on existing population configuration, the extent of vaccine implementation, as well as human mobility since it is a human obligate bacterium. To shed light on its international movement, we used extensive genome data from the Global Pneumococcal Sequencing project and estimated migration parameters between multiple countries in Africa. Data on allele frequencies of polymorphisms at housekeeping-like loci for multiple different lineages circulating in the populations of South Africa, Malawi, Kenya, and The Gambia were used to calculate the fixation index (Fst) between countries. We then further used these summaries to fit migration coalescent models with the likelihood-free inference algorithms available in the ELFI software package. Synthetic datawere additionally used to validate the inference approach. Our results demonstrate country-pair specific migration patterns and heterogeneity in the extent of migration between different lineages. Our approach demonstrates that coalescent models can be effectively used for inferring migration rates for bacterial species and lineages provided sufficiently granular population genomics surveillance data. Further, it can demonstrate the connectivity of respiratory disease agents between countries to inform intervention policy in the longer term.

Integrative methods reveal multiple drivers of diversification in rice paddy snakes

Divergence dating analyses in systematics provide a framework to develop and test biogeographic hypotheses regarding speciation. However, as molecular datasets grow from multilocus to genomic, sample sizes decrease due to computational burdens, and the testing of fine-scale biogeographic hypotheses becomes difficult. In this study, we use coalescent demographic models to investigate the diversification of poorly known rice paddy snakes from Southeast Asia (Homalopsidae: Hypsiscopus), which have conflicting dates of origin based on previous studies. We use coalescent modeling to test the hypothesis that Hypsiscopus diversified 2.5 mya during the Khorat Plateau uplift in Thailand. Additionally, we use ecological niche analyses to identify potential differences in the niche space of the two most widely distributed species in the past and present. Our results suggest Hypsiscopus diversified ~ 2.4 mya, supporting that the Khorat Plateau may have initiated the diversification of rice paddy snakes. We also find significant niche differentiation and shifts between species of Hypsiscopus, indicating that environmental differences may have sustained differentiation of this genus after the Khorat Plateau uplift. Our study expands on the diversification history of snakes in Southeast Asia, and highlights how results from smaller multilocus datasets can be useful in developing and testing biogeographic hypotheses alongside genomic datasets.

Machine learning approaches delimit cryptic taxa in a previously intractable species complex

Cryptic species are not diagnosable via morphological criteria, but can be detected through analysis of DNA sequences. A number of methods have been developed for identifying species based on genetic data; however, these methods are prone to over-splitting taxa with extreme population structure, such as dispersal-limited organisms. Machine learning methodologies have the potential to overcome this challenge. Here, we apply such approaches, using a large dataset generated through hybrid target enrichment of ultraconserved elements (UCEs). Our study taxon is the Aoraki denticulata species complex, a lineage of extremely low-dispersal arachnids endemic to the South Island of Aotearoa New Zealand. This group of mite harvesters has been the subject of previous species delimitation studies using smaller datasets generated through Sanger sequencing and analytical approaches that rely on multispecies coalescent models and barcoding gap discovery. Those analyses yielded a number of putative cryptic species that seems unrealistic and extreme, based on what we know about species’ geographic ranges and genetic diversity in non-cryptic mite harvesters. We find that machine learning approaches, on the other hand, identify cryptic species with geographic ranges that are similar to those seen in other morphologically diagnosable mite harvesters in Aotearoa New Zealand’s South Island. We performed both unsupervised and supervised machine learning analyses, the latter with training data drawn either from animals broadly (vagile and non-vagile) or from a custom training dataset from dispersal-limited harvesters. We conclude that applying machine learning approaches to the analysis of UCE-derived genetic data is an effective method for delimiting species in complexes of low-vagility cryptic species, and that the incorporation of training data from biologically relevant analogues can be critically informative.

Traditional multilocus phylogeny fails to fully resolve Palearctic ground squirrels (Spermophilus) relationships but reveals a new species endemic to West Siberia

Previous efforts to reconstruct evolutionary history of Palearctic ground squirrels within the genus Spermophilus have primarily relied on a single mitochondrial marker for phylogenetic data. In this study, we present the first phylogeny with comprehensive taxon sampling of Spermophilus via a conventional multilocus approach utilizing five mitochondrial and five nuclear markers. Through application of the multispecies coalescent model, we constructed a species tree revealing four distinct clades that diverged during the Late Miocene. These clades are 1) S. alaschanicus and S. dauricus from East Asia; 2) S. musicus and S. pygmaeus from East Europe and northwestern Central Asia; 3) the subgenus Colobotis found across Central Asia and its adjacent regions and encompassing S. brevicauda, S. erythrogenys, S. fulvus, S. major, S. pallidicauda, S. ralli, S. relictus, S. selevini, and S. vorontsovi sp. nov.; and 4) a Central/Eastern Europe and Asia Minor clade comprising S. citellus, S. taurensis, S. xanthoprymnus, S. suslicus, and S. odessanus. The latter clade lacked strong support owing to uncertainty of taxonomic placement of S. odessanus and S. suslicus. Resolving relationships within the subgenus Colobotis, which radiated rapidly, remains challenging likely because of incomplete lineage sorting and introgressive hybridization. Most of modern Spermophilus species diversified during the Early–Middle Pleistocene (2.2–1.0 million years ago). We propose a revised taxonomic classification for the genus Spermophilus by recognizing 18 species including a newly identified one (S. vorontsovi sp. nov.), which is found only in a limited area in the southeast of West Siberia. Employing genome-wide single-nucleotide polymorphism genotyping, we substantiated the role of the Ob River as a major barrier ensuring robust isolation of this taxon from S. erythrogenys. Despite its inherent limitations, the traditional multilocus approach remains a valuable tool for resolving relationships and can provide important insights into otherwise poorly understood groups. It is imperative to recognize that additional efforts are needed to definitively determine phylogenetic relationships between certain species of Palearctic ground squirrels.

Multi-charmed and singled charmed hadrons from coalescence: yields and ratios in different collision systems at LHC

We study the production of charmed and multi-charmed hadrons in ultra-relativistic Heavy Ion Collisions coupling the transport approach for charm dynamics in the medium to an hybrid hadronization model of coalescence plus fragmentation. In this paper, we mainly discuss the particle yields for single charmed and multi-charmed baryons focusing mainly on the production of Ξcc\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\Xi _{cc}$$\\end{document} and Ωccc\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\Omega _{ccc}$$\\end{document}. We provide first predictions for PbPb collision in 0-10%\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$0 \\! - \\! 10\\%\\ $$\\end{document}centrality class and then we explore the system size dependence through KrKr , to ArAr and OO collisions, planned within the ALICE3 experiment. In these cases, a monotonic behavior for the yields emerges which can be tested in future experimental data. We found about three order of magnitude increase in the production of Ωccc\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\Omega _{ccc}$$\\end{document} in PbPb collisions compared with the yield in small collision systems like OO collisions. Furthermore, we investigate the effects on the Ωccc\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\Omega _{ccc}$$\\end{document} particle yield and spectra coming from the modification of the charm quark distribution due to the different size of the collision systems also comparing it to the case of thermalized charm distributions. These results suggest that observation on the Ωccc\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\Omega _{ccc}$$\\end{document} spectra and their evolution across system size can give novel information about the partial thermalization of the charm quark distribution as well as to its wave function width. Furthermore, we find that the Ωccc\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\Omega _{ccc}$$\\end{document} /D0\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$D^0$$\\end{document} ratio is an observable more sensitive with respect to Λc\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\Lambda _c$$\\end{document} /D0\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$D^0$$\\end{document} , this ratio is predicted to span over two order of magnitude from large to small systems.

Study on condensation of oily particles under the influence of water vapor phase transition

To address the problem that oily fine particles are too small to be easily captured by traditional dust removal equipment, which hinders improving the purification efficiency, this paper proposes using the technical principle of water vapor phase transition to conduct experiments and simulations. A set of water vapor experimental devices for fine particles of oil mist is designed, and the data are analyzed by a Grimm spectrometer. The experimental results show the phase change of water vapor has a remarkable effect on the agglomeration of oily fine particles. Phase change condensation makes the particle size of fine particles increase by 5–10 times in a truly brief time, and the polymerization and removal efficiency are greatly improved. To analyze the changes of the coalescence characteristics in a larger range of parameter changes, the Euler multiphase flow model is adopted, the particle condensation growth rate function is introduced by a user-defined function to establish a coalescence model and the heterogeneous condensation growth characteristics of oily particles are analyzed in a supersaturated water vapor environment. The simulation results show that the condensing chamber temperature of the condensing chamber is 283 K, and the lower the air flow rate of the condensing chamber, the more favorable the heterogeneous condensation growth of oil particles. The change rule of agglomeration characteristics of oily particles obtained by experiment and simulation proves the feasibility of improving the removal efficiency of oily fine particles.

Bursts of coalescence within population pedigrees whenever big families occur.

We consider a simple diploid population-genetic model with potentially high variability of offspring numbers among individuals. Specifically, against a backdrop of Wright-Fisher reproduction and no selection, there is an additional probability that a big family occurs, meaning that a pair of individuals has a number of offspring on the order of the population size. We study how the pedigree of the population generated under this model affects the ancestral genetic process of a sample of size two at a single autosomal locus without recombination. Our population model is of the type for which multiple-merger coalescent processes have been described. We prove that the conditional distribution of the pairwise coalescence time given the random pedigree converges to a limit law as the population size tends to infinity. This limit law may or may not be the usual exponential distribution of the Kingman coalescent, depending on the frequency of big families. But because it includes the number and times of big families, it differs from the usual multiple-merger coalescent models. The usual multiple-merger coalescent models are seen as describing the ancestral process marginal to, or averaging over, the pedigree. In the limiting ancestral process conditional on the pedigree, the intervals between big families can be modeled using the Kingman coalescent but each big family causes a discrete jump in the probability of coalescence. Analogous results should hold for larger samples and other population models. We illustrate these results with simulations and additional analysis, highlighting their implications for inference and understanding of multilocus data.

Momentum correlation of light nuclei in Au + Au collisions at [formula omitted] GeV

Within the Ultra-relativistic Quantum Molecular Dynamics model (UrQMD) coupled with nucleon coalescence model and Mini-Spanning-Tree model, the yields of light nuclei have been stimulated in Au + Au collisions over an energy range of sNN=2.0∼7.7GeV and the momentum correlation functions of light nuclei pairs have been calculated by both the Lednický-Lyuboshitz and the Correlation After Burner methods. We compare the yields of light nuclei and their momentum correlation functions at midrapidity in this energy region with experimental data. It is found that there are differences between the results of the two models, and the coalescence method seems less valid at low collision energy. Furthermore, both the peak values of proton-proton correlation functions and the transition point of elliptic flows from out-of-plane to in-plane emission show a turning point around 3-4 GeV, which suggests that there is a relation between momentum correlation function and collective flow of particles.

Anomalous networks under the multispecies coalescent: theory and prevalence.

Reticulations in a phylogenetic network represent processes such as gene flow, admixture, recombination and hybrid speciation. Extending definitions from the tree setting, an anomalous network is one in which some unrooted tree topology displayed in the network appears in gene trees with a lower frequency than a tree not displayed in the network. We investigate anomalous networks under the Network Multispecies Coalescent Model with possible correlated inheritance at reticulations. Focusing on subsets of 4 taxa, we describe a new algorithm to calculate quartet concordance factors on networks of any level, faster than previous algorithms because of its focus on 4 taxa. We then study topological properties required for a 4-taxon network to be anomalous, uncovering the key role of [Formula: see text]-cycles: cycles of 3 edges parent to a sister group of 2 taxa. Under the model of common inheritance, that is, when each gene tree coalesces within a species tree displayed in the network, we prove that 4-taxon networks are never anomalous. Under independent and various levels of correlated inheritance, we use simulations under realistic parameters to quantify the prevalence of anomalous 4-taxon networks, finding that truly anomalous networks are rare. At the same time, however, we find a significant fraction of networks close enough to the anomaly zone to appear anomalous, when considering the quartet concordance factors observed from a few hundred genes. These apparent anomalies may challenge network inference methods.

Experimental Investigation of Void Coalescence in XTral-728 Plate Containing Three-Void Cluster

Void clustering effect is investigated by performing uniaxial tensile test on commercial 95% aluminum alloy XTral 728 by strategically placing an additional hole. Incorporating additional void induces strong stress and strain localization and reduces incipient coalescence strain approximately 17%. Also, two different cluster orientations with respect to applied loading and hole spacing ratios have been considered. The experimental revealed that material ductility is significantly decreased with increasing hole spacing ratios and for certain cluster orientation. To evaluate existing void coalescence models, numerical simulations are also performed and found that existing models overestimate the incipient coalescence strains for considered three-void cluster.

Phylogenomic analysis reveals five independently evolved African forage grass clades in the genus Urochloa.

The grass genus Urochloa (Brachiaria) sensu lato includes forage crops that are important for beef and dairy industries in tropical and sub-tropical Africa, South America and Oceania/Australia. Economically important species include U. brizantha, U. decumbens, U. humidicola, U. mutica, U. arrecta, U. trichopus, U. mosambicensis and Megathyrsus maximus, all native to the African continent. Perennial growth habits, large, fast growing palatable leaves, intra- and interspecific morphological variability, apomictic reproductive systems and frequent polyploidy are widely shared within the genus. The combination of these traits probably favoured the selection for forage domestication and weediness, but trait emergence across Urochloa cannot be modelled, as a robust phylogenetic assessment of the genus has not been conducted. We aim to produce a phylogeny for Urochloa that includes all important forage species, and identify their closest wild relatives (crop wild relatives). Finally, we will use our phylogeny and available trait data to infer the ancestral states of important forage traits across Urochloa s.l. and model the evolution of forage syndromes across the genus. Using a target enrichment sequencing approach (Angiosperm 353), we inferred a species-level phylogeny for Urochloa s.l., encompassing 54 species (~40 % of the genus) and outgroups. Phylogenies were inferred using a multispecies coalescent model and maximum likelihood method. We determined the phylogenetic placement of agriculturally important species and identified their closest wild relatives, or crop wild relatives, based on well-supported monophyly. Further, we mapped key traits associated with Urochloa forage crops to the species tree and estimated ancestral states for forage traits along branch lengths for continuous traits and at ancestral nodes in discrete traits. Agricultural species belong to five independent clades, including U. brizantha and U. decumbens lying in a previously defined species complex. Crop wild relatives were identified for these clades supporting previous sub-generic groupings in Urochloa based on morphology. Using ancestral trait estimation models, we find that five morphological traits that correlate with forage potential (perennial growth habits, culm height, leaf size, a winged rachis and large seeds) independently evolved in forage clades. Urochloa s.l. is a highly diverse genus that contains numerous species with agricultural potential, including crop wild relatives that are currently underexploited. All forage species and their crop wild relatives naturally occur on the African continent and their conservation across their native distributions is essential. Genomic and phenotypic diversity in forage clade species and their wild relatives need to be better assessed both to develop conservation strategies and to exploit the diversity in the genus for improved sustainability in Urochloa cultivar production.

Comparative phylogeography shows congruent co‐divergence in Neotropical seasonally dry forest and savanna tree species

AbstractAimThe response of South American seasonally dry tropical forest (SDTF) and savanna species to the Quaternary climate changes is still poorly understood. Here we use multi‐model inference approach to compare the phylogeography and demographic history of five trumpet tree species (also known as pau‐d'arco) and draw general biogeographical patterns: Handroanthus impetiginosus, H. ochraceus, H. serratifolius, Tabebuia aurea and T. roseoalba.LocationAmazonia, Atlantic Forest, Cerrado and Caatinga ecoregions, South America and Brazil.TaxaBignoniaceae, Tabebuia Alliance, trees.MethodsWe compiled genetic data, occurrence records and ecological niche modelling (ENM) from published articles. We used coalescent modelling to test patterns of co‐divergence, and analysed niche breadth and overlap and spatial patterns in genetic diversity using quantile regressions.ResultsSpecies showed concordant lineage divergence times in the Pleistocene. Lineage diversification and phylogeographical patterns are more related to species niche breadth than to biome (savanna or SDTF). Differences in genetic diversity among species matches the expected for the demographic scenarios obtained with coalescent simulations and ENMs: species with wider range expansion during the Last Glacial Maximum had higher diversity than species with range retraction or multiple refugia.Main ConclusionsDiversification burst in the Pliocene, and major diversification of lineages during the Pleistocene seems to be a general pattern in SDTF and savannas in the South America. The high number of haplotypes at the Central Brazil, corresponding to the Cerrado ecoregion, highlights the importance of Protected Areas for SDTFs and savannas, highly threatened biomes in Brazil.

Robustness of steady state and stochastic cyclicity in generalized coalescence-fragmentation models

Processes of coalescence and fragmentation are used to understand the time-evolution of the mass distribution of various systems and may result in a steady state or in stable deterministic or stochastic cycles. Motivated by applications in insurgency warfare we investigate coalescence-fragmentation systems. We begin with a simple model of size-biased coalescence accompanied by shattering into monomers. Depending on the parameters this model has an approximately power-law-distributed steady state or stochastic cycles of alternating gelation and shattering. We conduct stochastic simulations of this model and its generalizations to include different kernel types, accretion and erosion, and various distributions of non-shattering fragmentation. Our central aim is to explore the robustness of the steady state and gel-shatter stochastic cycles to these variations. We show that an approximate power-law steady state persists with the addition of accretion and erosion, and with partial rather than total shattering. However, broader distributions of fragment sizes typically vitiate both the power law steady state and gel-shatter cyclicity. This work clarifies features shown in coalescence/fragmentation model simulations and elucidates the relationship between the microscopic dynamics and observed phenomena in this widely applicable interdisciplinary model type.Graphical abstract

χc1(4274) multiplicity in heavy-ion collisions

In a previous work we computed the thermally averaged cross sections for the production and absorption reactions of the χc1(4274) state in the hot hadron gas formed in heavy-ion collisions. In the present work we estimate the final yield of the χc1(4274) state in these collisions. We use the coalescence model to fix the initial multiplicity of the χc1(4274), which is treated as a P-wave bound state of DsD¯s0 and also as a compact tetraquark. The Bjorken picture is used to model the hydrodynamic expansion and cooling. We also consider the Hubble cooling, which is faster and mimics the effect of the transverse expansion. Then, the kinetic equation is solved to evaluate the time evolution of the χc1(4274) yield during the hot hadron gas phase. Since the χc1(4274) decay width is large it might decay inside the hadron gas. Therefore, we also include the χc1(4274) decay and regeneration terms by means of an effective coupling, estimated from the available data. The combined effects of hadronic interactions and the χc1(4274) decay have a strong impact on the final yield. Also, predictions of the χc1(4274) multiplicity as a function of centrality and of the charged hadron multiplicity measured at midrapidity [dNch/dη(η<0.5)] are presented. Finally, we calculate the yield of a proposed P-wave molecular state of DsD¯s0, Y′(4274), characterized by a smaller width and smaller coupling constant obtained from the Weinberg compositeness condition. Published by the American Physical Society 2024

Past volcanic activity predisposes an endemic threatened seabird to negative anthropogenic impacts

Humans are regularly cited as the main driver of current biodiversity extinction, but the impact of historic volcanic activity is often overlooked. Pre-human evidence of wildlife abundance and diversity are essential for disentangling anthropogenic impacts from natural events. Réunion Island, with its intense and well-documented volcanic activity, endemic biodiversity, long history of isolation and recent human colonization, provides an opportunity to disentangle these processes. We track past demographic changes of a critically endangered seabird, the Mascarene petrel Pseudobulweria aterrima, using genome-wide SNPs. Coalescent modeling suggested that a large ancestral population underwent a substantial population decline in two distinct phases, ca. 125,000 and 37,000 years ago, coinciding with periods of major eruptions of Piton des Neiges. Subsequently, the ancestral population was fragmented into the two known colonies, ca. 1500 years ago, following eruptions of Piton de la Fournaise. In the last century, both colonies declined significantly due to anthropogenic activities, and although the species was initially considered extinct, it was rediscovered in the 1970s. Our findings suggest that the current conservation status of wildlife on volcanic islands should be firstly assessed as a legacy of historic volcanic activity, and thereafter by the increasing anthropogenic impacts, which may ultimately drive species towards extinction.

Ultrafast learning of four-node hybridization cycles in phylogenetic networks using algebraic invariants.

The abundance of gene flow in the Tree of Life challenges the notion that evolution can be represented with a fully bifurcating process which cannot capture important biological realities like hybridization, introgression, or horizontal gene transfer. Coalescent-based network methods are increasingly popular, yet not scalable for big data, because they need to perform a heuristic search in the space of networks as well as numerical optimization that can be NP-hard. Here, we introduce a novel method to reconstruct phylogenetic networks based on algebraic invariants. While there is a long tradition of using algebraic invariants in phylogenetics, our work is the first to define phylogenetic invariants on concordance factors (frequencies of four-taxon splits in the input gene trees) to identify level-1 phylogenetic networks under the multispecies coalescent model. Our novel hybrid detection methodology is optimization-free as it only requires the evaluation of polynomial equations, and as such, it bypasses the traversal of network space, yielding a computational speed at least 10 times faster than the fastest-to-date network methods. We illustrate our method's performance on simulated and real data from the genus Canis. We present an open-source publicly available Julia package PhyloDiamond.jl available at https://github.com/solislemuslab/PhyloDiamond.jl with broad applicability within the evolutionary community.

Computational understanding of the coalescence of metallic nanoparticles: a mini review.

Metallic nanoparticles exhibit extraordinary properties that differ from those of bulk materials due to their large surface area to volume ratios. Coalescence of metallic nanoparticles has a huge impact on their properties. Remarkable progress has been made by using computational methods for understanding nanoparticle coalescence. This work aims to provide a mini review on the state-of-the-art modelling and simulation of nanoparticle coalescence. First, we will discuss the outstanding performances and coalescence behaviors of metallic nanoparticles, and list some challenges in the coalescence of metallic nanoparticles. Next, we will introduce the applications of molecular dynamics and the Monte Carlo method in nanoparticle coalescence. Furthermore, we will discuss the coalescence kinetics and mechanisms of metal nanoparticles with the same element and different elements, alloy nanoparticles and metal oxide nanoparticles. Finally, we will present our perspective and conclusion.