Identifiability of Level-1 Species Networks from Gene Tree Quartets.

When hybridization or other forms of lateral gene transfer have occurred, evolutionary relationships of species are better represented by phylogenetic networks than by trees. While inference of such networks remains challenging, several recently proposed methods are based on quartet concordance factors-the probabilities that a tree relating a gene sampled from the species displays the possible 4-taxon relationships. Building on earlier results, we investigate what level-1 network features are identifiable from concordance factors under the network multispecies coalescent model.

Reaction rebalancing: a novel approach to curating reaction databases.

Purpose: Reaction databases are a key resource for a wide variety of applications in computational chemistry and biochemistry, including Computer-aided Synthesis Planning (CASP) and the large-scale analysis of metabolic networks. The full potential of these resources can only be realized if datasets are accurate and complete. Missing co-reactants and co-products, i.

Phylogenetic and Chemical Probing Information as Soft Constraints in RNA Secondary Structure Prediction.

Extrinsic, experimental information can be incorporated into thermodynamics-based RNA folding algorithms in the form of pseudo-energies. Evolutionary conservation of RNA secondary structure elements is detectable in alignments of phylogenetically related sequences and provides evidence for the presence of certain base pairs that can also be converted into pseudo-energy contributions. We show that the centroid base pairs computed from a consensus folding model such as RNAalifold result in a substantial improvement of the prediction accuracy for single sequences.

Comparative RNA Genomics.

Over the last quarter of a century it has become clear that RNA is much more than just a boring intermediate in protein expression. Ancient RNAs still appear in the core information metabolism and comprise a surprisingly large component in bacterial gene regulation. A common theme with these types of mostly small RNAs is their reliance of conserved secondary structures.

The Theory of Gene Family Histories.

Most genes are part of larger families of evolutionary-related genes. The history of gene families typically involves duplications and losses of genes as well as horizontal transfers into other organisms. The reconstruction of detailed gene family histories, i.

Assessing the Quality of Cotranscriptional Folding Simulations.

Structural changes in RNAs are an important contributor to controlling gene expression not only at the posttranscriptional stage but also during transcription. A subclass of riboswitches and RNA thermometers located in the 5' region of the primary transcript regulates the downstream functional unit - usually an ORF - through premature termination of transcription. Not only such elements occur naturally, but they are also attractive devices in synthetic biology.

Quantifying synergy and redundancy between networks.

Understanding how different networks relate to each other is key for understanding complex systems. We introduce an intuitive yet powerful framework to disentangle different ways in which networks can be (dis)similar and complementary to each other. We decompose the shortest paths between nodes as uniquely contributed by one source network, or redundantly by either, or synergistically by both together.

Cavity approach for the approximation of spectral density of graphs with heterogeneous structures.

Graphs have become widely used to represent and study social, biological, and technological systems. Statistical methods to analyze empirical graphs were proposed based on the graph's spectral density. However, their running time is cubic in the number of vertices, precluding direct application to large instances.

Landau Singularities Revisited: Computational Algebraic Geometry for Feynman Integrals.

We reformulate the analysis of singularities of Feynman integrals in a way that can be practically applied to perturbative computations in the standard model in dimensional regularization. After highlighting issues in the textbook treatment of Landau singularities, we develop an algorithm for classifying and computing them using techniques from computational algebraic geometry. We introduce an algebraic variety called the principal Landau determinant, which captures the singularities even in the presence of massless particles or UV/IR divergences.

Diversity begets stability: Sublinear growth and competitive coexistence across ecosystems.

The worldwide loss of species diversity brings urgency to understanding how diverse ecosystems maintain stability. Whereas early ecological ideas and classic observations suggested that stability increases with diversity, ecological theory makes the opposite prediction, leading to the long-standing "diversity-stability debate." Here, we show that this puzzle can be resolved if growth scales as a sublinear power law with biomass (exponent <1), exhibiting a form of population self-regulation analogous to models of individual ontogeny.

Tighter Constraints of Multi-Qubit Entanglement in Terms of Nonconvex Entanglement Measures LCREN and LCRENoA.

The monogamy property of entanglement is an intriguing feature of multipartite quantum entanglement. Most entanglement measures satisfying the monogamy inequality have turned out to be convex. Whether nonconvex entanglement measures obey the monogamy inequalities remains less known at present.

Relativistic Perturbation Theory for Black-Hole Boson Clouds.

We develop a relativistic perturbation theory for scalar clouds around rotating black holes. We first introduce a relativistic product and corresponding orthogonality relation between modes, extending a recent result for gravitational perturbations. We then derive the analog of time-dependent perturbation theory in quantum mechanics, and apply it to calculate self-gravitational frequency shifts.

Taking identity-by-descent analysis into the wild: Estimating realized relatedness in free-ranging macaques.

Biological relatedness is a key consideration in studies of behavior, population structure, and trait evolution. Except for parent-offspring dyads, pedigrees capture relatedness imperfectly. The number and length of DNA segments that are identical-by-descent (IBD) yield the most precise estimates of relatedness.

IDENTIFIABILITY OF LEVEL-1 SPECIES NETWORKS FROM GENE TREE QUARTETS.

When hybridization or other forms of lateral gene transfer have occurred, evolutionary relationships of species are better represented by phylogenetic networks than by trees. While inference of such networks remains challenging, several recently proposed methods are based on quartet concordance factors - the probabilities that a tree relating a gene sampled from the species displays the possible 4-taxon relationships. Building on earlier results, we investigate what level-1 network features are identifiable from concordance factors under the network multispecies coalescent model.

Geometry and convergence of natural policy gradient methods.

We study the convergence of several natural policy gradient (NPG) methods in infinite-horizon discounted Markov decision processes with regular policy parametrizations. For a variety of NPGs and reward functions we show that the trajectories in state-action space are solutions of gradient flows with respect to Hessian geometries, based on which we obtain global convergence guarantees and convergence rates. In particular, we show linear convergence for unregularized and regularized NPG flows with the metrics proposed by Kakade and Morimura and co-authors by observing that these arise from the Hessian geometries of conditional entropy and entropy respectively.

Proteinortho6: pseudo-reciprocal best alignment heuristic for graph-based detection of (co-)orthologs.

Proteinortho is a widely used tool to predict (co)-orthologous groups of genes for any set of species. It finds application in comparative and functional genomics, phylogenomics, and evolutionary reconstructions. With a rapidly increasing number of available genomes, the demand for large-scale predictions is also growing.

Conditioning in Tropical Probability Theory.

We define a natural operation of conditioning of tropical diagrams of probability spaces and show that it is Lipschitz continuous with respect to the asymptotic entropy distance.

Arrow Contraction and Expansion in Tropical Diagrams.

Arrow contraction applied to a tropical diagram of probability spaces is a modification of the diagram, replacing one of the morphisms with an isomorphism while preserving other parts of the diagram. It is related to the rate regions introduced by Ahlswede and Körner. In a companion article, we use arrow contraction to derive information about the shape of the entropic cone.

Master regulators of biological systems in higher dimensions.

A longstanding goal of biology is to identify the key genes and species that critically impact evolution, ecology, and health. Network analysis has revealed keystone species that regulate ecosystems and master regulators that regulate cellular genetic networks. Yet these studies have focused on pairwise biological interactions, which can be affected by the context of genetic background and other species present, generating higher-order interactions.

MuDoGeR: Multi-Domain Genome recovery from metagenomes made easy.

Several computational frameworks and workflows that recover genomes from prokaryotes, eukaryotes and viruses from metagenomes exist. Yet, it is difficult for scientists with little bioinformatics experience to evaluate quality, annotate genes, dereplicate, assign taxonomy and calculate relative abundance and coverage of genomes belonging to different domains. MuDoGeR is a user-friendly tool tailored for those familiar with Unix command-line environment that makes it easy to recover genomes of prokaryotes, eukaryotes and viruses from metagenomes, either alone or in combination.

Generalized Bayes approach to inverse problems with model misspecification.

We propose a general framework for obtaining probabilistic solutions to PDE-based inverse problems. Bayesian methods are attractive for uncertainty quantification but assume knowledge of the likelihood model or data generation process. This assumption is difficult to justify in many inverse problems, where the specification of the data generation process is not obvious.

Relative timing information and orthology in evolutionary scenarios.

Background: Evolutionary scenarios describing the evolution of a family of genes within a collection of species comprise the mapping of the vertices of a gene tree T to vertices and edges of a species tree S. The relative timing of the last common ancestors of two extant genes (leaves of T) and the last common ancestors of the two species (leaves of S) in which they reside is indicative of horizontal gene transfers (HGT) and ancient duplications. Orthologous gene pairs, on the other hand, require that their last common ancestors coincides with a corresponding speciation event.

The AnimalAssociatedMetagenomeDB reveals a bias towards livestock and developed countries and blind spots in functional-potential studies of animal-associated microbiomes.

Background: Metagenomic data can shed light on animal-microbiome relationships and the functional potential of these communities. Over the past years, the generation of metagenomics data has increased exponentially, and so has the availability and reusability of data present in public repositories. However, identifying which datasets and associated metadata are available is not straightforward.

Evolution as a result of resource flow in ecosystems: Ecological dynamics can drive evolution.

To see how the flow of energy across ecosystems can derive evolution, I introduce a framework in which individuals interact with their peers and environment to accumulate resources, and use the resources to pay for their metabolic costs, grow and reproduce. I show that two conservation principles determine the system's equilibrium state: conservation of resources- a physical principle stating that in the equilibrium, resource production and consumption should balance, and payoff equality- an economic principle, stating that the payoffs of different types in equilibrium should equal. Besides the equilibrium state, the system shows non-equilibrium fluctuations derived by the exponential growth of the individuals in which the payoff equality principle does not hold.

The human cell count and size distribution.

Cell size and cell count are adaptively regulated and intimately linked to growth and function. Yet, despite their widespread relevance, the relation between cell size and count has never been formally examined over the whole human body. Here, we compile a comprehensive dataset of cell size and count over all major cell types, with data drawn from >1,500 published sources.