Linking great apes genome evolution across time scales using polymorphism-aware phylogenetic models

The genomes of related species contain valuable information on the history of the considered taxa. Great apes in particular exhibit variation of evolutionary patterns along their genomes. However, the great ape data also bring new challenges, such as the presence of incomplete lineage sorting and ancestral shared polymorphisms. Previous methods for genome-scale analysis are restricted to very few individuals or cannot disentangle the contribution of mutation rates and fixation biases. This represents a limitation both for the understanding of these forces as well as for the detection of regions affected by selection. Here, we present a new model designed to estimate mutation rates and fixation biases from genetic variation within and between species. We relax the assumption of instantaneous substitutions, modeling substitutions as mutational events followed by a gradual fixation. Hence, we straightforwardly account for shared ancestral polymorphisms and incomplete lineage sorting. We analyze genome-wide synonymous site alignments of human, chimpanzee, and two orangutan species. From each taxon, we include data from several individuals. We estimate mutation rates and GC-biased gene conversion intensity. We find that both mutation rates and biased gene conversion vary with GC content. We also find lineage-specific differences, with weaker fixation biases in orangutan species, suggesting a reduced historical effective population size. Finally, our results are consistent with directional selection acting on coding sequences in relation to exonic splicing enhancers.Publisher PDFPeer reviewe

PoMo : an allele frequency-based approach for species tree estimation

This work was supported by a grant from the Austrian Science Fund (FWF, P24551-B25 to C.K.). N.D.M. and D.S. were members of the Vienna Graduate School of Population Genetics which is supported by a grant of the Austrian Science Fund (FWF, W1225-B20). N.D.M. was partially supported by the Institute for Emerging Infections, funded by the Oxford Martin School.Incomplete lineage sorting can cause incongruencies of the overall species-level phylogenetic tree with the phylogenetic trees for individual genes or genomic segments. If these incongruencies are not accounted for, it is possible to incur several biases in species tree estimation. Here, we present a simple maximum likelihood approach that accounts for ancestral variation and incomplete lineage sorting. We use a POlymorphisms-aware phylogenetic MOdel (PoMo) that we have recently shown to efficiently estimate mutation rates and fixation biases from within and between-species variation data. We extend this model to perform efficient estimation of species trees. We test the performance of PoMo in several different scenarios of incomplete lineage sorting using simulations and compare it with existing methods both in accuracy and computational speed. In contrast to other approaches, our model does not use coalescent theory but is allele frequency based. We show that PoMo is well suited for genome-wide species tree estimation and that on such data it is more accurate than previous approaches.Publisher PDFPeer reviewe

Disease and Disorders of Freshwater Unionid Mussels: A Brief Overview of Recent Studies

The use of aquatic invertebrates in biomedical research and as environmental sentinels has dramatically grown in recent decades, with an increased need in understanding of comparative pathology. The Unionids freshwater mussels are a group of worldwide distributed bivalves residing small ditches and ponds, lakes, canals and rivers, often used as animal test in eco-toxicological studies. Once one of the most abundant bivalve molluscs in ancient rivers around the world, now many of them are declining in many countries and consequently are nearly extinct in many areas. The causes of this decline are not fully understood but alteration and degradation of the freshwater habitat seemed to play a central role. To date, link causality to the observed losses during episode of mussel die-offs has been more difficult to establish, and disease and pathogen presence have been scarcely considered. In this article we provide a brief overview of unionids freshwater mussel conservation status, also describing reported diseases and pathogens and illustrating a few relatively well-documented studies

Estimating empirical codon hidden Markov models

Empirical codon models (ECMs) estimated from a large number of globular protein families outperformed mechanistic codon models in their description of the general process of protein evolution. Among other factors, ECMs implicitly model the influence of amino acid properties and multiple nucleotide substitutions (MNS). However, the estimation of ECMs requires large quantities of data, and until recently, only few suitable data sets were available. Here, we take advantage of several new Drosophila species genomes to estimate codon models from genome-wide data. The availability of large numbers of genomes over varying phylogenetic depths in the Drosophila genus allows us to explore various divergence levels. In consequence, we can use these data to determine the appropriate level of divergence for the estimation of ECMs, avoiding overestimation of MNS rates caused by saturation. To account for variation in evolutionary rates along the genome, we develop new empirical codon hidden Markov models (ecHMMs). These models significantly outperform previous ones with respect to maximum likelihood values, suggesting that they provide a better fit to the evolutionary process. Using ECMs and ecHMMs derived from genome-wide data sets, we devise new likelihood ratio tests (LRTs) of positive selection. We found classical LRTs very sensitive to the presence of MNSs, showing high false-positive rates, especially with small phylogenies. The new LRTs are more conservative than the classical ones, having acceptable false-positive rates and reduced power.Publisher PDFPeer reviewe

Short-range template switching in great ape genomes explored using pair hidden Markov models.

Many complex genomic rearrangements arise through template switch errors, which occur in DNA replication when there is a transient polymerase switch to an alternate template nearby in three-dimensional space. While typically investigated at kilobase-to-megabase scales, the genomic and evolutionary consequences of this mutational process are not well characterised at smaller scales, where they are often interpreted as clusters of independent substitutions, insertions and deletions. Here we present an improved statistical approach using pair hidden Markov models, and use it to detect and describe short-range template switches underlying clusters of mutations in the multi-way alignment of hominid genomes. Using robust statistics derived from evolutionary genomic simulations, we show that template switch events have been widespread in the evolution of the great apes' genomes and provide a parsimonious explanation for the presence of many complex mutation clusters in their phylogenetic context. Larger-scale mechanisms of genome rearrangement are typically associated with structural features around breakpoints, and accordingly we show that atypical patterns of secondary structure formation and DNA bending are present at the initial template switch loci. Our methods improve on previous non-probabilistic approaches for computational detection of template switch mutations, allowing the statistical significance of events to be assessed. By specifying realistic evolutionary parameters based on the genomes and taxa involved, our methods can be readily adapted to other intra- or inter-species comparisons

CMAPLE: Efficient Phylogenetic Inference in the Pandemic Era

We have recently introduced MAPLE (MAximum Parsimonious Likelihood Estimation), a new pandemic-scale phylogenetic inference method exclusively designed for genomic epidemiology. In response to the need for enhancing MAPLE's performance and scalability, here we present two key components: (i) CMAPLE software, a highly optimized C++ reimplementation of MAPLE with many new features and advancements, and (ii) CMAPLE library, a suite of application programming interfaces to facilitate the integration of the CMAPLE algorithm into existing phylogenetic inference packages. Notably, we have successfully integrated CMAPLE into the widely used IQ-TREE 2 software, enabling its rapid adoption in the scientific community. These advancements serve as a vital step toward better preparedness for future pandemics, offering researchers powerful tools for large-scale pathogen genomic analysis

The surgical approach of late-onset tracheoesophageal fistula in a tracheostomized COVID-19 patient

In the COVID-19 era the tracheal complications due to prolonged mechanical ventilation have significantly increased. Acquired tracheoesophageal fistula is one of those in ventilated COVID-19 patients. Thus, the knowledge of their management in such fragile patient is crucial. We report a case of tracheoesophageal fistula in a 56-year-old female under prolonged mechanical ventilation for COVID-19 bilateral pneumonia and discuss its management. A surgical approach was proposed. By a collar-shaped transverse cervicotomic access, we transected the trachea at level of fistula en-bloc with the tracheostoma. The esophageal lesion was longitudinally repaired in two-layers. Protective left strap muscle was sandwiched between esophagus and trachea. The tracheal end-to-end anastomosis was completed without a re-tracheostoma. Even if surgical approach of tracheoesophageal fistula in COVID-19 patients has not been tested before, surgery remains the treatment of choice according to the multidisciplinary board

Reversible polymorphism-aware phylogenetic models and their application to tree inference

This work is supported by the Austrian Science Fund (FWF-P24551 and I-2805-B29) and partially by the Vienna Graduate School of Population Genetics (FWF-W1225).We present a reversible Polymorphism-Aware Phylogenetic Model (revPoMo) for species tree estimation from genome-wide data. revPoMo enables the reconstruction of large scale species trees for many within-species samples. It expands the alphabet of DNA substitution models to include polymorphic states, thereby, naturally accounting for incomplete lineage sorting. We implemented revPoMo in the maximum likelihood software IQ-TREE. A simulation study and an application to great apes data show that the runtimes of our approach and standard substitution models are comparable but that revPoMo has much better accuracy in estimating trees, divergence times and mutation rates. The advantage of revPoMo is that an increase of sample size per species improves estimations but does not increase runtime. Therefore, revPoMo is a valuable tool with several applications, from speciation dating to species tree reconstruction.Publisher PDFPeer reviewe

Haplotype characterization of a stranded common minke whale calf (Balaenoptera acutorostrata lacépède, 1804): Is the mediterranean sea a potential calving or nursery ground for the species?

The stranding of a suckling calf of Common Minke Whale (Balaenoptera acutorostrata) on the coast near Salerno (Campania, Southern Italy) is reported. The molecular analysis of a partial sequence of the mitochondrial DNA control region shows that the animal bore a haplotype identical to haplotype Ba169 considered as typical of individuals from North Atlantic population. Historical data and our results suggest the possibility that the Mediterranean Sea might be a potential calving or nursery ground for this species