Maximum likelihood estimation of the heterogeneity of substitution rate among nucleotide sites

Genomics and evolutionary analysis of Chlorella variabilis-infecting viruses demarcate criteria for defining species of giant viruses

Chloroviruses exhibit a close relationship with their hosts with the phenotypic aspect of their ability to form lytic plaques having primarily guided the taxonomy. However, with the isolation of viruses that are only able to complete their replication cycle in one strain of Chlorella variabilis, systematic challenges emerged. In this study, we described the genomic features of 53 new chlorovirus isolates and used them to elucidate part of the evolutionary history and taxonomy of this clade. Our analysis revealed new chloroviruses with the largest genomes to date (>400 kbp) and indicated that four genomic features are statistically different in the viruses that only infect the Syngen 2-3 strain of C. variabilis (OSy viruses). We found large regions of dissimilarity in the genomes of viruses PBCV-1 and OSy-NE5 when compared with the other genomes. These regions contained genes related to the interaction with the host cell machinery and viral capsid proteins, which provided insights into the evolution of the replicative and structural modules in these giant viruses. Phylogenetic analysis using hallmark genes of Nucleocytoviricota revealed that OSy-viruses evolved from the NC64A-viruses, possibly emerging as a result of the strict relationship with their hosts. Merging phylogenetics and nucleotide identity analyses, we propose strategies to demarcate viral species, resulting in seven new species of chloroviruses. Collectively, our results show how genomic data can be used as lines of evidence to demarcate viral species. Using the chloroviruses as a case study, we expect that similar initiatives will emerge using the basis exhibited here.IMPORTANCEChloroviruses are a group of giant viruses with long dsDNA genomes that infect different species of Chlorella-like green algae. They are host-specific, and some isolates can only replicate within a single strain of Chlorella variabilis. The genomics of these viruses is still poorly explored, and the characterization of new isolates provides important data on their genetic diversity and evolution. In this work, we describe 53 new chlorovirus genomes, including many isolated from alkaline lakes for the first time. Through comparative genomics and molecular phylogeny, we provide evidence of genomic gigantism in chloroviruses and show that a subset of viruses became highly specific for their hosts at a particular point in evolutionary history. We propose criteria to demarcate species of chloroviruses, paving the way for an update in the taxonomy of other groups of viruses. This study is a new and important piece in the complex puzzle of giant algal viruses.

The Toxoplasma gondii mitochondrial transporter ABCB7L is essential for the biogenesis of cytosolic and nuclear iron-sulfur cluster proteins and cytosolic translation

Iron-sulfur (Fe-S) clusters are ubiquitous inorganic cofactors required for numerous essential cellular pathways. Since they cannot be scavenged from the environment, Fe-S clusters are synthesized de novo in cellular compartments such as the apicoplast, mitochondrion, and cytosol. The cytosolic Fe-S cluster biosynthesis pathway relies on the transport of an intermediate from the mitochondrial pathway. An ATP-binding cassette (ABC) transporter called ABCB7 is responsible for this role in numerous commonly studied organisms, but its role in the medically important apicomplexan parasites has not yet been studied. Here we identify and characterize a Toxoplasma gondii ABCB7 homolog, which we name ABCB7-like (ABCB7L). Genetic depletion shows that it is essential for parasite growth and that its disruption triggers partial stage conversion. Characterization of the knock-down line highlights a defect in the biogenesis of cytosolic and nuclear Fe-S proteins leading to defects in protein translation and other pathways including DNA and RNA replication and metabolism. Our work provides support for a broad conservation of the connection between mitochondrial and cytosolic pathways in Fe-S cluster biosynthesis and reveals its importance for parasite survival.

IMPORTANCE

Iron-sulfur (Fe-S) clusters are inorganic cofactors of proteins that play key roles in numerous essential biological processes, for example, respiration and DNA replication. Cells possess dedicated biosynthetic pathways to assemble Fe-S clusters, including a pathway in the mitochondrion and cytosol. A single transporter, called ABCB7, connects these two pathways, allowing an essential intermediate generated by the mitochondrial pathway to be used in the cytosolic pathway. Cytosolic and nuclear Fe-S proteins are dependent on the mitochondrial pathway, mediated by ABCB7, in numerous organisms studied to date. Here, we study the role of a homolog of ABCB7, which we name ABCB7-like (ABCB7L), in the ubiquitous unicellular apicomplexan parasite Toxoplasma gondii. We generated a depletion mutant of Toxoplasma ABCB7L and showed its importance for parasite fitness. Using comparative quantitative proteomic analysis and experimental validation of the mutants, we show that ABCB7L is required for cytosolic and nuclear, but not mitochondrial, Fe-S protein biogenesis. Our study supports the conservation of a protein homologous to ABCB7 and which has a similar function in apicomplexan parasites and provides insight into an understudied aspect of parasite metabolism.

Identification of key drivers of antimicrobial resistance in Enterococcus using machine learning

With antimicrobial resistance (AMR) rapidly evolving in pathogens, quick and accurate identification of genetic determinants of phenotypic resistance is essential for improving surveillance, stewardship, and clinical mitigation. Machine learning (ML) models show promise for AMR prediction in diagnostics but require a deep understanding of internal processes to use effectively. Our study utilised AMR gene, pangenomic, and predicted plasmid features from 647 Enterococcus faecium and Enterococcus faecalis genomes across the One Health continuum, along with corresponding resistance phenotypes, to develop interpretive ML classifiers. Vancomycin resistance could be predicted with 99% accuracy with AMR gene features, 98% with pangenome features, and 96% with plasmid clusters. Top pangenome features overlapped with the resistance genes of the vanA operon, which are often laterally transmitted via plasmids. Doxycycline resistance prediction achieved approximately 92% accuracy with pangenome features, with the top feature being elements of Tn916 conjugative transposon, a tet(M) carrier. Erythromycin resistance prediction models achieved about 90% accuracy, but top features were negatively correlated with resistance due to the confounding effect of population structure. This work demonstrates the importance of reviewing ML models' features to discern biological relevance even when achieving high-performance metrics. Our workflow offers the potential to propose hypotheses for experimental testing, enhancing the understanding of AMR mechanisms, which are crucial for combating the AMR crisis.

An independent base composition of each rate class for improved likelihood-based phylogeny estimation; the 5rf model

The combination of a time reversible Markov process with a "hidden" mixture of gamma distributed relative site rates plus invariant sites have become the most favoured options for likelihood and other probabilistic models of nucleotide evolution (e.g., tr4gi which approximates a gamma with four rate classes). However, these models assume a homogeneous and stationary distribution of nucleotide (character or base) frequencies. Here, we explore the potential benefits and pitfalls of allowing each rate category (rate class) of a 4gi mixture model to have its own base frequencies. This is achieved by starting each of the five rate classes, at the tree's root, with its own free choice of nucleotide frequencies to create a 4gi5rf model or a 5rf model in shorthand. We assess the practical identifiability of this approach with a BEAST 2 implementation, aiming to determine if it can accurately estimate credibility intervals and expected values for a wide range of plausible parameter values. Practical identifiability, as distinguished from mathematical identifiability, gauges the model's ability to identify parameters in real-world scenarios, as opposed to theoretically with infinite data. One of the most common types of phylogenetic data is mitochondrial DNA (mtDNA) protein coding sequence. It is often assumed current models analyse robustly such data and that higher likelihood/posterior probability models do better. However, this abstract shows that vertebrate mtDNA remains a very difficult type of data to fully model, and that dramatically higher likelihoods do not mean a model is measurably more accurate with respect to recovering key parameters of biological interest (e.g., monophyletic groups, their support and their ages). The 4gi5rf model considerably improves marginal likelihoods and seems to reverse some apparent errors exacerbated by the 4gi model, while introducing others. Problems appear to be linked to non-stationary DNA repair processes that alter the mutation/substitution spectra across lineages and time. We also show such problems are not unique to mtDNA and are encountered in analysing nuclear sequences. Non-stationarity of DNA repair processes mutation/substitution spectra thus pose an active challenge to obtaining reliable inferences of relationships and divergence times near the root of placental mammals, for example. An open source implementation is available under the LGPL 3.0 license in the beastbooster package for BEAST 2, available from https://github.com/rbouckaert/beastbooster.

Mechanistic and evolutionary insights into isoform-specific 'supercharging' in DCLK family kinases

Catalytic signaling outputs of protein kinases are dynamically regulated by an array of structural mechanisms, including allosteric interactions mediated by intrinsically disordered segments flanking the conserved catalytic domain. The doublecortin-like kinases (DCLKs) are a family of microtubule-associated proteins characterized by a flexible C-terminal autoregulatory 'tail' segment that varies in length across the various human DCLK isoforms. However, the mechanism whereby these isoform-specific variations contribute to unique modes of autoregulation is not well understood. Here, we employ a combination of statistical sequence analysis, molecular dynamics simulations, and in vitro mutational analysis to define hallmarks of DCLK family evolutionary divergence, including analysis of splice variants within the DCLK1 sub-family, which arise through alternative codon usage and serve to 'supercharge' the inhibitory potential of the DCLK1 C-tail. We identify co-conserved motifs that readily distinguish DCLKs from all other calcium calmodulin kinases (CAMKs), and a 'Swiss Army' assembly of distinct motifs that tether the C-terminal tail to conserved ATP and substrate-binding regions of the catalytic domain to generate a scaffold for autoregulation through C-tail dynamics. Consistently, deletions and mutations that alter C-terminal tail length or interfere with co-conserved interactions within the catalytic domain alter intrinsic protein stability, nucleotide/inhibitor binding, and catalytic activity, suggesting isoform-specific regulation of activity through alternative splicing. Our studies provide a detailed framework for investigating kinome-wide regulation of catalytic output through cis-regulatory events mediated by intrinsically disordered segments, opening new avenues for the design of mechanistically divergent DCLK1 modulators, stabilizers, or degraders.

Mechanistic and evolutionary insights into isoform-specific 'supercharging' in DCLK family kinases

Catalytic signaling outputs of protein kinases are dynamically regulated by an array of structural mechanisms, including allosteric interactions mediated by intrinsically disordered segments flanking the conserved catalytic domain. The Doublecortin Like Kinases (DCLKs) are a family of microtubule-associated proteins characterized by a flexible C-terminal autoregulatory 'tail' segment that varies in length across the various human DCLK isoforms. However, the mechanism whereby these isoform-specific variations contribute to unique modes of autoregulation is not well understood. Here, we employ a combination of statistical sequence analysis, molecular dynamics simulations and in vitro mutational analysis to define hallmarks of DCLK family evolutionary divergence, including analysis of splice variants within the DCLK1 sub-family, which arise through alternative codon usage and serve to 'supercharge' the inhibitory potential of the DCLK1 C-tail. We identify co-conserved motifs that readily distinguish DCLKs from all other Calcium Calmodulin Kinases (CAMKs), and a 'Swiss-army' assembly of distinct motifs that tether the C-terminal tail to conserved ATP and substrate-binding regions of the catalytic domain to generate a scaffold for auto-regulation through C-tail dynamics. Consistently, deletions and mutations that alter C-terminal tail length or interfere with co-conserved interactions within the catalytic domain alter intrinsic protein stability, nucleotide/inhibitor-binding, and catalytic activity, suggesting isoform-specific regulation of activity through alternative splicing. Our studies provide a detailed framework for investigating kinome-wide regulation of catalytic output through cis-regulatory events mediated by intrinsically disordered segments, opening new avenues for the design of mechanistically-divergent DCLK1 modulators, stabilizers or degraders.

First report of the beetle Henosepilachna nana (Kapur, 1950) (Coleoptera: Coccinellidae) from Maharashtra with special reference to molecular phylogeny and host plants

A ladybird beetle, collected from different localities of Kolhapur and Satara districts (Maharashtra) was identified as Henosepilachna nana (Kapur, 1950). The presence of this species in Maharashtra considerably extends its range to the north-west by about 700 km. Since this species is found on vegetables, Pumpkin (Cucurbita spp.), Cucumber (Cucumis spp.), and Karit fruit plant (Cucumis spp.). It is being recorded for the first time that this species is a pest of these vegetables. Also, molecular phylogeny has been studied for the first time in this species in which, this species is the sister taxon of Henosepilachna boisduvali. This species has been described briefly with colour photographs of male genitalia, female coxites and the 6th sternal plate of female, tarsi, & pronotum. Henosepilachna nana is being reported for the first time from Maharashtra and now the molecular data of this species is available.

dN/dS-H, a New Test to Distinguish Different Selection Modes in Protein Evolution and Cancer Evolution

One of the most popular measures in the analysis of protein sequence evolution is the ratio of nonsynonymous distance (d_N) to synonymous distance (d_S). Under the assumption that synonymous substitutions in the coding region are selectively neutral, the d_N/d_S ratio can be used to statistically detect the adaptive evolution (or purifying selection) if d_N/d_S > 1 (or d_N/d_S < 1) significantly. However, due to strong structural constraints and/or variable functional constraints imposed on amino acid sites, most encoding genes in most species have demonstrated d_N/d_S < 1. Consequently, the statistical power for testing d_N/d_S = 1 may be insufficient to distinguish between different selection modes. In this paper, we propose a more powerful test, called d_N/d_S-H, in which a new parameter H, a relative measure of rate variation among sites, was introduced. Given the condition of strong purifying selections at some sites, the d_N/d_S-H model predicts d_N/d_S = 1-H for neutral evolution, d_N/d_S < 1-H for nearly neutral selection, and d_N/d_S > 1-H for adaptive evolution. The potential of this new method for resolving the neutral-adaptive debates is illustrated by the protein sequence evolution in vertebrates, Drosophila and yeasts, as well as somatic cancer evolution (specialized as the C_N/C_S-H test).

Determining the International Spread of B.1.1.523 SARS-CoV-2 Lineage with a Set of Mutations Highly Associated with Reduced Immune Neutralization

Here, we report the emergence of the variant lineage B.1.1.523 that contains a set of mutations including 156_158del, E484K and S494P in the spike protein. E484K and S494P are known to significantly reduce SARS-CoV-2 neutralization by convalescent and vaccinated sera and are considered as mutations of concern. Lineage B.1.1.523 presumably originated in the Russian Federation and spread across European countries with the peak of transmission in April-May 2021. The B.1.1.523 lineage has now been reported from 31 countries. In this article, we analyze the possible origin of this mutation subset and its immune response using in silico methods.

Multifunctional Enzyme with Endoglucanase and Alginase/Glucuronan Lyase Activities from Bacterium Cellulophaga lytica

Cellulophaga lytica is a Gram-negative aerobic bacterium in the genome of which there are many genes encoding polysaccharide degrading enzymes. One of the enzymes named ClGP contains a glycoside hydrolase domain from the GH5 family and a polysaccharide lyase domain from the PL31 family. The enzyme also contains the TAT signaling peptide and the TIGR04183 domain that indicates extracellular nature of the enzyme. Phylogenetic analysis has shown that the enzymes most closely related to ClGP and containing all four domains (TAT, GH5, PL31, TIGR04183) are widespread among bacterial species belonging to the Flavobacteriaceae family. ClGP produced by the recombinant strain of E. coli was purified and characterized. ClGP exhibited activity of endoglucanase (EC 3.2.1.4) and catalyzed hydrolysis of β-D-glucan, carboxymethyl cellulose sodium salt (CMC-Na), and amorphous cellulose, but failed to hydrolyze microcrystalline cellulose and xylan. Products of CMC hydrolysis were cellobiose and cellotriose, whereas β-D-glucan was hydrolyzed to glucose, cellobiose, cellotetraose, and cellopentaose. ClGP was more active against the poly-β-D-mannuronate blocks than against the poly-α-L-glucuronate blocks of alginic acid. This indicates that the enzyme is a polyM lyase (EC 4.2.2.3). ClGP was active against polyglucuronic acid, so it displayed a glucuronan lyase (EC 4.2.2.14) activity. The enzyme had a neutral pH-optimum, was stable in the pH range 6.0-8.0, and displayed moderate thermal stability. ClGP effectively saccharified two species of brown algae, Saccharina latissima and Laminaria digitata, that suggests its potential for use in the production of biofuel from macroalgae.

A New Species of Paradiplozoon (Monogenea: Diplozoidae), A Gill Parasite of the Schizothorax Fish (Cyprinidae: Schizothoracinae) from the Yarkand River, Xinjiang, China

INTRODUCTION

Monogeneans of the genus Paradiplozoon were found on the gills of specimens of five species of schizothoracid caught using fyke nets in the upper stream of the Yarkand River, Xinjiang, China in May-August 2019.

METHODS

The preserved parasite were stained with boric acid magenta and hematoxylin, respectively. Morphological observations, line drawings, photomicrographs and measurements were made in Nikon ECLIPSE E200 imaging optical microscope and digitally edited. The molecular analysis included the study of the sequence of the second internal transcribed spacer (ITS 2) of the ribosomal DNA region, calculation and analysis of genetic distance, with phylogenetic reconstructions based on the Bayesian inference and Maximum Likelihood analysis.

RESULTS

The natural infection rate of host fish was 10-88%. Morphological analysis indicated that the average length of the new species was 2.125 mm while the width was 0.69 mm. The anterior part was 1.387 mm in length and the average length of the posterior part was 0.545 mm. The vitellaria was well-developed and located in the front of the body. A single ovary (oval shaped) was located at the back end of the reproductive binding area. A testis (irregular mass) was located behind or parallelled to the ovary. The new species can be distinguished from all the recorded Paradiplozoon species in terms of morphological characteristics such as haptor, clamp and central hook morphology, intestine shape and body size. In addition, the second internal transcribed spacer (ITS 2) of the ribosomal DNA region of the diplozoid was compared with that of known diplozoids previously published. It indicated that there were significant differences between the new species and the published diplozoids.

CONCLUSION

Both morphological and molecular analysis support that the diplozoid is a new species. Based on the sampling location, the new species was named Paradiplozoon yarkandense n. sp.

Evolutionary history and divergence times of Odonata (dragonflies and damselflies) revealed through transcriptomics

Dragonflies and damselflies are among the earliest flying insects with extant representatives. However, unraveling details of their long evolutionary history, such as egg laying (oviposition) strategies, is impeded by unresolved phylogenetic relationships, particularly in damselflies. Here we present a transcriptome-based phylogenetic reconstruction of Odonata, analyzing 2,980 protein-coding genes in 105 species representing nearly all the order's families. All damselfly and most dragonfly families are recovered as monophyletic. Our data suggest a sister relationship between dragonfly families of Gomphidae and Petaluridae. According to our divergence time estimates, both crown-Zygoptera and -Anisoptera arose during the late Triassic. Egg-laying with a reduced ovipositor apparently evolved in dragonflies during the late Jurassic/early Cretaceous. Lastly, we also test the impact of fossil choice and placement, particularly, of the extinct fossil species, †Triassolestodes asiaticus, and †Proterogomphus renateae on divergence time estimates. We find placement of †Proterogomphus renateae to be much more impactful than †Triassolestodes asiaticus.

Genomic and Phenotypic Analysis of COVID-19-Associated Pulmonary Aspergillosis Isolates of Aspergillus fumigatus

The ongoing global pandemic caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for coronavirus disease 2019 (COVID-19), first described in Wuhan, China. A subset of COVID-19 patients has been reported to have acquired secondary infections by microbial pathogens, such as opportunistic fungal pathogens from the genus Aspergillus. To gain insight into COVID-19-associated pulmonary aspergillosis (CAPA), we analyzed the genomes and characterized the phenotypic profiles of four CAPA isolates of Aspergillus fumigatus obtained from patients treated in the area of North Rhine-Westphalia, Germany. By examining the mutational spectrum of single nucleotide polymorphisms, insertion-deletion polymorphisms, and copy number variants among 206 genes known to modulate A. fumigatus virulence, we found that CAPA isolate genomes do not exhibit significant differences from the genome of the Af293 reference strain. By examining a number of factors, including virulence in an invertebrate moth model, growth in the presence of osmotic, cell wall, and oxidative stressors, secondary metabolite biosynthesis, and the MIC of antifungal drugs, we found that CAPA isolates were generally, but not always, similar to A. fumigatus reference strains Af293 and CEA17. Notably, CAPA isolate D had more putative loss-of-function mutations in genes known to increase virulence when deleted. Moreover, CAPA isolate D was significantly more virulent than the other three CAPA isolates and the A. fumigatus reference strains Af293 and CEA17, but similarly virulent to two other clinical strains of A. fumigatus. These findings expand our understanding of the genomic and phenotypic characteristics of isolates that cause CAPA. IMPORTANCE The global pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the etiological agent of coronavirus disease 2019 (COVID-19), has already killed millions of people. COVID-19 patient outcome can be further complicated by secondary infections, such as COVID-19-associated pulmonary aspergillosis (CAPA). CAPA is caused by Aspergillus fungal pathogens, but there is little information about the genomic and phenotypic characteristics of CAPA isolates. We conducted genome sequencing and extensive phenotyping of four CAPA isolates of Aspergillus fumigatus from Germany. We found that CAPA isolates were often, but not always, similar to other reference strains of A. fumigatus across 206 genetic determinants of infection-relevant phenotypes, including virulence. For example, CAPA isolate D was more virulent than other CAPA isolates and reference strains in an invertebrate model of fungal disease, but similarly virulent to two other clinical strains. These results expand our understanding of COVID-19-associated pulmonary aspergillosis.

Exploring Codon Adjustment Strategies towards Escherichia coli-Based Production of Viral Proteins Encoded by HTH1, a Novel Prophage of the Marine Bacterium Hypnocyclicus thermotrophus

Marine viral sequence space is immense and presents a promising resource for the discovery of new enzymes interesting for research and biotechnology. However, bottlenecks in the functional annotation of viral genes and soluble heterologous production of proteins hinder access to downstream characterization, subsequently impeding the discovery process. While commonly utilized for the heterologous expression of prokaryotic genes, codon adjustment approaches have not been fully explored for viral genes. Herein, the sequence-based identification of a putative prophage is reported from within the genome of Hypnocyclicus thermotrophus, a Gram-negative, moderately thermophilic bacterium isolated from the Seven Sisters hydrothermal vent field. A prophage-associated gene cluster, consisting of 46 protein coding genes, was identified and given the proposed name Hypnocyclicus thermotrophus phage H1 (HTH1). HTH1 was taxonomically assigned to the viral family Siphoviridae, by lowest common ancestor analysis of its genome and phylogeny analyses based on proteins predicted as holin and DNA polymerase. The gene neighbourhood around the HTH1 lytic cassette was found most similar to viruses infecting Gram-positive bacteria. In the HTH1 lytic cassette, an N-acetylmuramoyl-L-alanine amidase (Amidase_2) with a peptidoglycan binding motif (LysM) was identified. A total of nine genes coding for enzymes putatively related to lysis, nucleic acid modification and of unknown function were subjected to heterologous expression in Escherichia coli. Codon optimization and codon harmonization approaches were applied in parallel to compare their effects on produced proteins. Comparison of protein yields and thermostability demonstrated that codon optimization yielded higher levels of soluble protein, but codon harmonization led to proteins with higher thermostability, implying a higher folding quality. Altogether, our study suggests that both codon optimization and codon harmonization are valuable approaches for successful heterologous expression of viral genes in E. coli, but codon harmonization may be preferable in obtaining recombinant viral proteins of higher folding quality.

An evolutionary genomic approach reveals both conserved and species-specific genetic elements related to human disease in closely related Aspergillus fungi

Aspergillosis is an important opportunistic human disease caused by filamentous fungi in the genus Aspergillus. Roughly 70% of infections are caused by Aspergillus fumigatus, with the rest stemming from approximately a dozen other Aspergillus species. Several of these pathogens are closely related to A. fumigatus and belong in the same taxonomic section, section Fumigati. Pathogenic species are frequently most closely related to nonpathogenic ones, suggesting Aspergillus pathogenicity evolved multiple times independently. To understand the repeated evolution of Aspergillus pathogenicity, we performed comparative genomic analyses on 18 strains from 13 species, including 8 species in section Fumigati, which aimed to identify genes, both ones previously connected to virulence as well as ones never before implicated, whose evolution differs between pathogens and nonpathogens. We found that most genes were present in all species, including approximately half of those previously connected to virulence, but a few genes were section- or species-specific. Evolutionary rate analyses identified over 1700 genes whose evolutionary rate differed between pathogens and nonpathogens and dozens of genes whose rates differed between specific pathogens and the rest of the taxa. Functional testing of deletion mutants of 17 transcription factor-encoding genes whose evolution differed between pathogens and nonpathogens identified eight genes that affect either fungal survival in a model of phagocytic killing, host survival in an animal model of fungal disease, or both. These results suggest that the evolution of pathogenicity in Aspergillus involved both conserved and species-specific genetic elements, illustrating how an evolutionary genomic approach informs the study of fungal disease.

Felsenstein Phylogenetic Likelihood

In 1981, the Journal of Molecular Evolution (JME) published an article entitled "Evolutionary trees from DNA sequences: A maximum likelihood approach" by Joseph (Joe) Felsenstein (J Mol Evol 17:368-376, 1981). This groundbreaking work laid the foundation for the emerging field of statistical phylogenetics, providing a tractable way of finding maximum likelihood (ML) estimates of evolutionary trees from DNA sequence data. This paper is the second most cited (more than 9000 citations) in JME after Kimura's (J Mol Evol 16:111-120, 1980) seminal paper on a model of nucleotide substitution (with nearly 20,000 citations). On the occasion of the 50th anniversary of JME, we elaborate on the significance of Felsenstein's ML approach to estimating phylogenetic trees.

IQ-TREE 2: New Models and Efficient Methods for Phylogenetic Inference in the Genomic Era

IQ-TREE (http://www.iqtree.org, last accessed February 6, 2020) is a user-friendly and widely used software package for phylogenetic inference using maximum likelihood. Since the release of version 1 in 2014, we have continuously expanded IQ-TREE to integrate a plethora of new models of sequence evolution and efficient computational approaches of phylogenetic inference to deal with genomic data. Here, we describe notable features of IQ-TREE version 2 and highlight the key advantages over other software.

QMaker: Fast and accurate method to estimate empirical models of protein evolution

Amino acid substitution models play a crucial role in phylogenetic analyses. Maximum likelihood (ML) methods have been proposed to estimate amino acid substitution models; however, they are typically complicated and slow. In this article, we propose QMaker, a new ML method to estimate a general time-reversible \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$Q$\end{document} matrix from a large protein data set consisting of multiple sequence alignments. QMaker combines an efficient ML tree search algorithm, a model selection for handling the model heterogeneity among alignments, and the consideration of rate mixture models among sites. We provide QMaker as a user-friendly function in the IQ-TREE software package (http://www.iqtree.org) supporting the use of multiple CPU cores so that biologists can easily estimate amino acid substitution models from their own protein alignments. We used QMaker to estimate new empirical general amino acid substitution models from the current Pfam database as well as five clade-specific models for mammals, birds, insects, yeasts, and plants. Our results show that the new models considerably improve the fit between model and data and in some cases influence the inference of phylogenetic tree topologies.[Amino acid replacement matrices; amino acid substitution models; maximum likelihood estimation; phylogenetic inferences.]

Morphological Characters Can Strongly Influence Early Animal Relationships Inferred from Phylogenomic Data Sets

There are considerable phylogenetic incongruencies between morphological and phylogenomic data for the deep evolution of animals. This has contributed to a heated debate over the earliest-branching lineage of the animal kingdom: the sister to all other Metazoa (SOM). Here, we use published phylogenomic data sets ($\sim $45,000-400,000 characters in size with $\sim $15-100 taxa) that focus on early metazoan phylogeny to evaluate the impact of incorporating morphological data sets ($\sim $15-275 characters). We additionally use small exemplar data sets to quantify how increased taxon sampling can help stabilize phylogenetic inferences. We apply a plethora of common methods, that is, likelihood models and their "equivalent" under parsimony: character weighting schemes. Our results are at odds with the typical view of phylogenomics, that is, that genomic-scale data sets will swamp out inferences from morphological data. Instead, weighting morphological data 2-10$\times $ in both likelihood and parsimony can in some cases "flip" which phylum is inferred to be the SOM. This typically results in the molecular hypothesis of Ctenophora as the SOM flipping to Porifera (or occasionally Placozoa). However, greater taxon sampling improves phylogenetic stability, with some of the larger molecular data sets ($>$200,000 characters and up to $\sim $100 taxa) showing node stability even with $\geqq100\times $ upweighting of morphological data. Accordingly, our analyses have three strong messages. 1) The assumption that genomic data will automatically "swamp out" morphological data is not always true for the SOM question. Morphological data have a strong influence in our analyses of combined data sets, even when outnumbered thousands of times by molecular data. Morphology therefore should not be counted out a priori. 2) We here quantify for the first time how the stability of the SOM node improves for several genomic data sets when the taxon sampling is increased. 3) The patterns of "flipping points" (i.e., the weighting of morphological data it takes to change the inferred SOM) carry information about the phylogenetic stability of matrices. The weighting space is an innovative way to assess comparability of data sets that could be developed into a new sensitivity analysis tool. [Metazoa; Morphology; Phylogenomics; Weighting.].

Prevalence and genetic diversity of avian haemosporidian parasites in wild bird species of the order Columbiformes

Diseases can play a role in species decline. Among them, haemosporidian parasites, vector-transmitted protozoan parasites, are known to constitute a risk for different avian species. However, the magnitude of haemosporidian infection in wild columbiform birds, including strongly decreasing European turtle doves, is largely unknown. We examined the prevalence and diversity of haemosporidian parasites Plasmodium, Leucocytozoon and subgenera Haemoproteus and Parahaemoproteus in six species of the order Columbiformes during breeding season and migration by applying nested PCR, one-step multiplex PCR assay and microscopy. We detected infections in 109 of the 259 screened individuals (42%), including 15 distinct haemosporidian mitochondrial cytochrome b lineages, representing five H. (Haemoproteus), two H. (Parahaemoproteus), five Leucocytozoon and three Plasmodium lineages. Five of these lineages have never been described before. We discriminated between single and mixed infections and determined host species-specific prevalence for each parasite genus. Observed differences among sampled host species are discussed with reference to behavioural characteristics, including nesting and migration strategy. Our results support previous suggestions that migratory birds have a higher prevalence and diversity of blood parasites than resident or short-distance migratory species. A phylogenetic reconstruction provided evidence for H. (Haemoproteus) as well as H. (Parahaemoproteus) infections in columbiform birds. Based on microscopic examination, we quantified parasitemia, indicating the probability of negative effects on the host. This study provides a large-scale baseline description of haemosporidian infections of wild birds belonging to the order Columbiformes sampled in the northern hemisphere. The results enable the monitoring of future changes in parasite transmission areas, distribution and diversity associated with global change, posing a potential risk for declining avian species as the European turtle dove.

Comparative transcriptomics of ice-crawlers demonstrates cold specialization constrains niche evolution in a relict lineage

Key changes in ecological niche space are often critical to understanding how lineages diversify during adaptive radiations. However, the converse, or understanding why some lineages are depauperate and relictual, is more challenging, as many factors may constrain niche evolution. In the case of the insect order Grylloblattodea, highly conserved thermal breadth is assumed to be closely tied to their relictual status, but has not been formerly tested. Here, we investigate whether evolutionary constraints in the physiological tolerance of temperature can help explain relictualism in this lineage. Using a comparative transcriptomics approach, we investigate gene expression following acute heat and cold stress across members of Grylloblattodea and their sister group, Mantophasmatodea. We additionally examine patterns of protein evolution, to identify candidate genes of positive selection. We demonstrate that cold specialization in Grylloblattodea has been accompanied by the loss of the inducible heat shock response under both acute heat and cold stress. Additionally, there is widespread evidence of selection on protein-coding genes consistent with evolutionary constraints due to cold specialization. This includes positive selection on genes involved in trehalose transport, metabolic function, mitochondrial function, oxygen reduction, oxidative stress, and protein synthesis. These patterns of molecular adaptation suggest that Grylloblattodea have undergone evolutionary trade-offs to survive in cold habitats and should be considered highly vulnerable to climate change. Finally, our transcriptomic data provide a robust backbone phylogeny for generic relationships within Grylloblattodea and Mantophasmatodea. Major phylogenetic splits in each group relate to arid conditions driving biogeographical patterns, with support for a sister-group relationship between North American Grylloblatta and Altai-Sayan Grylloblattella, and a range disjunction in Namibia splitting major clades within Mantophasmatodea.

Classification and phylogeny for the annotation of novel eukaryotic GNAT acetyltransferases

The enzymes of the GCN5-related N-acetyltransferase (GNAT) superfamily count more than 870 000 members through all kingdoms of life and share the same structural fold. GNAT enzymes transfer an acyl moiety from acyl coenzyme A to a wide range of substrates including aminoglycosides, serotonin, glucosamine-6-phosphate, protein N-termini and lysine residues of histones and other proteins. The GNAT subtype of protein N-terminal acetyltransferases (NATs) alone targets a majority of all eukaryotic proteins stressing the omnipresence of the GNAT enzymes. Despite the highly conserved GNAT fold, sequence similarity is quite low between members of this superfamily even when substrates are similar. Furthermore, this superfamily is phylogenetically not well characterized. Thus functional annotation based on sequence similarity is unreliable and strongly hampered for thousands of GNAT members that remain biochemically uncharacterized. Here we used sequence similarity networks to map the sequence space and propose a new classification for eukaryotic GNAT acetyltransferases. Using the new classification, we built a phylogenetic tree, representing the entire GNAT acetyltransferase superfamily. Our results show that protein NATs have evolved more than once on the GNAT acetylation scaffold. We use our classification to predict the function of uncharacterized sequences and verify by in vitro protein assays that two fungal genes encode NAT enzymes targeting specific protein N-terminal sequences, showing that even slight changes on the GNAT fold can lead to change in substrate specificity. In addition to providing a new map of the relationship between eukaryotic acetyltransferases the classification proposed constitutes a tool to improve functional annotation of GNAT acetyltransferases.

Enzymes of the GCN5-related N-acetyltransferase (GNAT) superfamily transfer an acetyl group from one molecule to another. This reaction is called acetylation and is one of the most common reactions inside the cell. The GNAT superfamily counts more than 870 000 members through all kingdoms of life. Despite sharing the same fold the GNAT superfamily is very diverse in terms of amino acid sequence and substrates. The eight N-terminal acetyltransferases (NatA, NatB, etc.. to NatH) are a GNAT subtype which acetylates the free amine group of polypeptide chains. This modification is called N-terminal acetylation and is one of the most abundant protein modifications in eukaryotic cells. This subtype is also characterized by a high sequence diversity even though they share the same substrate. In addition, the phylogeny of the superfamily is not characterized. This hampers functional annotation based on sequence similarity, and discovery of novel NATs. In this work we set out to solve the problem of the classification of eukaryotic GCN5-related acetyltransferases and report the first classification framework of the superfamily. This framework can be used as a tool for annotation of all GCN5-related acetyltransferases. As an example of what can be achieved we report in this paper the computational prediction and in vitro verification of the function of two previously uncharacterized N-terminal acetyltransferases. We also report the first acetyltransferase phylogenetic tree of the GCN5 superfamily. It indicates that N-terminal acetyltransferases do not constitute one homogeneous protein family, but that the ability to bind and acetylate protein N-termini had evolved more than once on the same acetylation scaffold. We also show that even small changes in key positions can lead to altered enzyme specificity.

Quest for the Best Evolutionary Model

In the early 1980s, DNA sequencing became a routine and the increasing computing power opened the door to reconstruct molecular phylogenies using probabilistic approaches. DNA sequence alignments provided a large number of positions containing phylogenetic information, which could be extracted using explicit statistical models that described the mutation process using appropriate parameters. Consequently, an active quest started for building increasingly improved (more realistic) statistical models of nucleotide substitution. The simplest model assumed that nucleotide frequencies were in equilibrium and one single category of substitutions. Subsequent models allowed either unequal nucleotide frequencies or separate rates for transitions and transversions. The HKY85 model (Hasegawa et al. in J Mol Evol 22:160, 1985) combined elegantly both options into a single model, which became one of the most useful ones and has been the choice in many molecular phylogenetic studies ever since. The use of improved substitution models such as HKY85 allows reconstructing more accurate and reliable phylogenies, which in turn provide robust frameworks for understanding how biological diversity evolved and for performing a wealth of comparative studies in different disciplines such as ecology, biogeography, developmental biology, biochemistry, genomics, epidemiology, and biomedicine.

Genomic and phenotypic analysis of COVID-19-associated pulmonary aspergillosis isolates of Aspergillus fumigatus

The ongoing global pandemic caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for the coronavirus disease 2019 (COVID-19) first described from Wuhan, China. A subset of COVID-19 patients has been reported to have acquired secondary infections by microbial pathogens, such as fungal opportunistic pathogens from the genus Aspergillus . To gain insight into COVID-19 associated pulmonary aspergillosis (CAPA), we analyzed the genomes and characterized the phenotypic profiles of four CAPA isolates of Aspergillus fumigatus obtained from patients treated in the area of North Rhine-Westphalia, Germany. By examining the mutational spectrum of single nucleotide polymorphisms, insertion-deletion polymorphisms, and copy number variants among 206 genes known to modulate A. fumigatus virulence, we found that CAPA isolate genomes do not exhibit major differences from the genome of the Af293 reference strain. By examining virulence in an invertebrate moth model, growth in the presence of osmotic, cell wall, and oxidative stressors, and the minimum inhibitory concentration of antifungal drugs, we found that CAPA isolates were generally, but not always, similar to A. fumigatus reference strains Af293 and CEA17. Notably, CAPA isolate D had more putative loss of function mutations in genes known to increase virulence when deleted (e.g., in the FLEA gene, which encodes a lectin recognized by macrophages). Moreover, CAPA isolate D was significantly more virulent than the other three CAPA isolates and the A. fumigatus reference strains tested. These findings expand our understanding of the genomic and phenotypic characteristics of isolates that cause CAPA.

Blood Parasites in Endangered Wildlife-Trypanosomes Discovered During a Survey of Haemoprotozoa from the Tasmanian Devil

The impact of emerging infectious diseases is increasingly recognised as a major threat to wildlife. Wild populations of the endangered Tasmanian devil, Sarcophilus harrisii, are experiencing devastating losses from a novel transmissible cancer, devil facial tumour disease (DFTD); however, despite the rapid decline of this species, there is currently no information on the presence of haemoprotozoan parasites. In the present study, 95 Tasmanian devil blood samples were collected from four populations in Tasmania, Australia, which underwent molecular screening to detect four major groups of haemoprotozoa: (i) trypanosomes, (ii) piroplasms, (iii) Hepatozoon, and (iv) haemosporidia. Sequence results revealed Trypanosoma infections in 32/95 individuals. Trypanosoma copemani was identified in 10 Tasmanian devils from three sites and a second Trypanosoma sp. was identified in 22 individuals that were grouped within the poorly described T. cyclops clade. A single blood sample was positive for Babesia sp., which most closely matched Babesia lohae. No other blood protozoan parasite DNA was detected. This study provides the first insight into haemoprotozoa from the Tasmanian devil and the first identification of Trypanosoma and Babesia in this carnivorous marsupial.

Adaptive evolution shapes the present-day distribution of the thermal sensitivity of population growth rate

Developing a thorough understanding of how ectotherm physiology adapts to different thermal environments is of crucial importance, especially in the face of global climate change. A key aspect of an organism's thermal performance curve (TPC)-the relationship between fitness-related trait performance and temperature-is its thermal sensitivity, i.e., the rate at which trait values increase with temperature within its typically experienced thermal range. For a given trait, the distribution of thermal sensitivities across species, often quantified as "activation energy" values, is typically right-skewed. Currently, the mechanisms that generate this distribution are unclear, with considerable debate about the role of thermodynamic constraints versus adaptive evolution. Here, using a phylogenetic comparative approach, we study the evolution of the thermal sensitivity of population growth rate across phytoplankton (Cyanobacteria and eukaryotic microalgae) and prokaryotes (bacteria and archaea), 2 microbial groups that play a major role in the global carbon cycle. We find that thermal sensitivity across these groups is moderately phylogenetically heritable, and that its distribution is shaped by repeated evolutionary convergence throughout its parameter space. More precisely, we detect bursts of adaptive evolution in thermal sensitivity, increasing the amount of overlap among its distributions in different clades. We obtain qualitatively similar results from evolutionary analyses of the thermal sensitivities of 2 physiological rates underlying growth rate: net photosynthesis and respiration of plants. Furthermore, we find that these episodes of evolutionary convergence are consistent with 2 opposing forces: decrease in thermal sensitivity due to environmental fluctuations and increase due to adaptation to stable environments. Overall, our results indicate that adaptation can lead to large and relatively rapid shifts in thermal sensitivity, especially in microbes for which rapid evolution can occur at short timescales. Thus, more attention needs to be paid to elucidating the implications of rapid evolution in organismal thermal sensitivity for ecosystem functioning.

Marked mitochondrial genetic variation in individuals and populations of the carcinogenic liver fluke Clonorchis sinensis

Clonorchiasis is a neglected tropical disease caused by the Chinese liver fluke, Clonorchis sinensis, and is often associated with a malignant form of bile duct cancer (cholangiocarcinoma). Although some aspects of the epidemiology of clonorchiasis are understood, little is known about the genetics of C. sinensis populations. Here, we conducted a comprehensive genetic exploration of C. sinensis from endemic geographic regions using complete mitochondrial protein gene sets. Genomic DNA samples from C. sinensis individuals (n = 183) collected from cats and dogs in China (provinces of Guangdong, Guangxi, Hunan, Heilongjiang and Jilin) as well as from rats infected with metacercariae from cyprinid fish from the Russian Far East (Primorsky Krai region) were deep sequenced using the BGISEQ-500 platform. Informatic analyses of mitochondrial protein gene data sets revealed marked genetic variation within C. sinensis; significant variation was identified within and among individual worms from distinct geographical locations. No clear affiliation with a particular location or host species was evident, suggesting a high rate of dispersal of the parasite across endemic regions. The present work provides a foundation for future biological, epidemiological and ecological studies using mitochondrial protein gene data sets, which could aid in elucidating associations between particular C. sinensis genotypes/haplotypes and the pathogenesis or severity of clonorchiasis and its complications (including cholangiocarcinoma) in humans.

Clonorchiasis is an important neglected tropical disease caused by the Chinese liver fluke, Clonorchis sinensis, which can induce malignant bile duct cancer (cholangiocarcinoma). Little precise information is available on the biology, epidemiology and population genetics of C. sinensis. For this reason, we explored here the genetic composition of C. sinensis populations in distinct endemic areas in China and Russia. Using a deep sequencing-informatic approach, we revealed marked mitochondrial genetic variation within and between individuals and populations of C. sinensis, with no particular affiliation with geographic or host origin. These molecular findings and the methodology established should underpin future genetic studies of C. sinensis causing human clonorchiasis and associated complications (cancer) as well as transmission patterns in endemic regions.

The Aspergillus fumigatus transcription factor RglT is important for gliotoxin biosynthesis and self-protection, and virulence

Aspergillus fumigatus is an opportunistic fungal pathogen that secretes an array of immune-modulatory molecules, including secondary metabolites (SMs), which contribute to enhancing fungal fitness and growth within the mammalian host. Gliotoxin (GT) is a SM that interferes with the function and recruitment of innate immune cells, which are essential for eliminating A. fumigatus during invasive infections. We identified a C6 Zn cluster-type transcription factor (TF), subsequently named RglT, important for A. fumigatus oxidative stress resistance, GT biosynthesis and self-protection. RglT regulates the expression of several gli genes of the GT biosynthetic gene cluster, including the oxidoreductase-encoding gene gliT, by directly binding to their respective promoter regions. Subsequently, RglT was shown to be important for virulence in a chemotherapeutic murine model of invasive pulmonary aspergillosis (IPA). Homologues of RglT and GliT are present in eurotiomycete and sordariomycete fungi, including the non-GT-producing fungus A. nidulans, where a conservation of function was described. Phylogenetically informed model testing led to an evolutionary scenario in which the GliT-based resistance mechanism is ancestral and RglT-mediated regulation of GliT occurred subsequently. In conclusion, this work describes the function of a previously uncharacterised TF in oxidative stress resistance, GT biosynthesis and self-protection in both GT-producing and non-producing Aspergillus species.

Animal Disease Surveillance in the 21st Century: Applications and Robustness of Phylodynamic Methods in Recent U.S. Human-Like H3 Swine Influenza Outbreaks

Emerging and endemic animal viral diseases continue to impose substantial impacts on animal and human health. Most current and past molecular surveillance studies of animal diseases investigated spatio-temporal and evolutionary dynamics of the viruses in a disjointed analytical framework, ignoring many uncertainties and made joint conclusions from both analytical approaches. Phylodynamic methods offer a uniquely integrated platform capable of inferring complex epidemiological and evolutionary processes from the phylogeny of viruses in populations using a single Bayesian statistical framework. In this study, we reviewed and outlined basic concepts and aspects of phylodynamic methods and attempted to summarize essential components of the methodology in one analytical pipeline to facilitate the proper use of the methods by animal health researchers. Also, we challenged the robustness of the posterior evolutionary parameters, inferred by the commonly used phylodynamic models, using hemagglutinin (HA) and polymerase basic 2 (PB2) segments of the currently circulating human-like H3 swine influenza (SI) viruses isolated in the United States and multiple priors. Subsequently, we compared similarities and differences between the posterior parameters inferred from sequence data using multiple phylodynamic models. Our suggested phylodynamic approach attempts to reduce the impact of its inherent limitations to offer less biased and biologically plausible inferences about the pathogen evolutionary characteristics to properly guide intervention activities. We also pinpointed requirements and challenges for integrating phylodynamic methods in routine animal disease surveillance activities.

Fulvivirga aurantia sp. nov. and Xanthovirga aplysinae gen. nov., sp. nov., marine bacteria isolated from the sponge Aplysina fistularis, and emended description of the genus Fulvivirga

Two Gram-stain-negative, strictly aerobic, marine bacteria, designated as strains RKSG066^T and RKSG123^T, were isolated from a sponge Aplysina fistularis collected at a depth of 15 m off the west coast of San Salvador, The Bahamas. Investigation of nearly full-length 16S rRNA gene and whole genome-based phylogenies revealed that both strains belong to the order Cytophagales within the class Cytophagia and phylum Bacteroidetes. Strain RKSG066^T formed a monophyletic clade with described members of the genus Fulvivirga, while strain RKSG123^T formed a well-supported paraphyletic branch apart from this and other related genera within the family Flammeovirgaceae. For both RKSG066^T and RKSG123^T, optimal growth parameters were 30-37 °C, pH 7-8 and 2-3 % (w/v) NaCl; cells were catalase- and oxidase-positive, and flexirubin-type pigments were absent. The predominant fatty acids were iso-C_15 : 0, C_16 : 0, C_18 : 0, iso-C_17 : 0 3-OH, C_16 : 1 ω5c, iso-C_15 : 0 3-OH, C_18 : 1 ω9c and iso-C_15 : 1 G for RKSG066^T, and iso-C_17 : 0 3-OH, C_16 : 1 ω5c, iso-C_15 : 0, C_16 : 0 3-OH and summed feature 4 (iso-C_17 : 1 I and/or anteiso-C_17 : 1 B) for RKSG123^T. Menaquinone-7 was the major respiratory quinone for both strains. The DNA G+C contents of RKSG066^T and RKSG123^T were 39.5 and 36.7 mol%, respectively. On the basis of phylogenetic distinctiveness and polyphasic analysis, the type strain RKSG066^T (=TSD-73^T=LMG 29870^T) is proposed to represent a novel species of the genus Fulvivirga, for which the name Fulvivirga aurantia sp. nov. is proposed. The type strain RKSG123^T (=TSD-75^T=LMG 30075^T) is proposed to represent the type species of a novel genus and species with the proposed name Xanthovirga aplysinae gen. nov., sp. nov. Additionally, the genus Fulvivirga is emended to include strains of orange-pigmented colonies that contain the predominant cellular fatty acids C_16 : 0, C_18 : 0, C_{16  :  1} ω5c and C_{18  :  1} ω9c.

Phytoplankton thermal responses adapt in the absence of hard thermodynamic constraints

To better predict how populations and communities respond to climatic temperature variation, it is necessary to understand how the shape of the response of fitness-related rates to temperature evolves (the thermal performance curve). Currently, there is disagreement about the extent to which the evolution of thermal performance curves is constrained. One school of thought has argued for the prevalence of thermodynamic constraints through enzyme kinetics, whereas another argues that adaptation can-at least partly-overcome such constraints. To shed further light on this debate, we perform a phylogenetic meta-analysis of the thermal performance curves of growth rate of phytoplankton-a globally important functional group-controlling for environmental effects (habitat type and thermal regime). We find that thermodynamic constraints have a minor influence on the shape of the curve. In particular, we detect a very weak increase of maximum performance with the temperature at which the curve peaks, suggesting a weak "hotter-is-better" constraint. Also, instead of a constant thermal sensitivity of growth across species, as might be expected from strong constraints, we find that all aspects of the thermal performance curve evolve along the phylogeny. Our results suggest that phytoplankton thermal performance curves adapt to thermal environments largely in the absence of hard thermodynamic constraints.

Unique k-mers as Strain-Specific Barcodes for Phylogenetic Analysis and Natural Microbiome Profiling

The need for a comparative analysis of natural metagenomes stimulated the development of new methods for their taxonomic profiling. Alignment-free approaches based on the search for marker k-mers turned out to be capable of identifying not only species, but also strains of microorganisms with known genomes. Here, we evaluated the ability of genus-specific k-mers to distinguish eight phylogroups of Escherichia coli (A, B1, C, E, D, F, G, B2) and assessed the presence of their unique 22-mers in clinical samples from microbiomes of four healthy people and four patients with Crohn's disease. We found that a phylogenetic tree inferred from the pairwise distance matrix for unique 18-mers and 22-mers of 124 genomes was fully consistent with the topology of the tree, obtained with concatenated aligned sequences of orthologous genes. Therefore, we propose strain-specific "barcodes" for rapid phylotyping. Using unique 22-mers for taxonomic analysis, we detected microbes of all groups in human microbiomes; however, their presence in the five samples was significantly different. Pointing to the intraspecies heterogeneity of E. coli in the natural microflora, this also indicates the feasibility of further studies of the role of this heterogeneity in maintaining population homeostasis.

TMRS: an algorithm for computing the time to the most recent substitution event from a multiple alignment column

Background 

As the number of sequenced genomes grows, researchers have access to an increasingly rich source for discovering detailed evolutionary information. However, the computational technologies for inferring biologically important evolutionary events are not sufficiently developed.

Results 

We present algorithms to estimate the evolutionary time (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$t_{\text {MRS}}$$\end{document}tMRS) to the most recent substitution event from a multiple alignment column by using a probabilistic model of sequence evolution. As the confidence in estimated \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$t_{\text {MRS}}$$\end{document}tMRS values varies depending on gap fractions and nucleotide patterns of alignment columns, we also compute the standard deviation \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma$$\end{document}σ of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$t_{\text {MRS}}$$\end{document}tMRS by using a dynamic programming algorithm. We identified a number of human genomic sites at which the last substitutions occurred between two speciation events in the human lineage with confidence. A large fraction of such sites have substitutions that occurred between the concestor nodes of Hominoidea and Euarchontoglires. We investigated the correlation between tissue-specific transcribed enhancers and the distribution of the sites with specific substitution time intervals, and found that brain-specific transcribed enhancers are threefold enriched in the density of substitutions in the human lineage relative to expectations.

Conclusions 

We have presented algorithms to estimate the evolutionary time (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$t_{\text {MRS}}$$\end{document}tMRS) to the most recent substitution event from a multiple alignment column by using a probabilistic model of sequence evolution. Our algorithms will be useful for Evo-Devo studies, as they facilitate screening potential genomic sites that have played an important role in the acquisition of unique biological features by target species.

Electronic supplementary material

The online version of this article (10.1186/s13015-019-0158-3) contains supplementary material, which is available to authorized users.

Mobile genetic elements explain size variation in the mitochondrial genomes of four closely-related Armillaria species

Background

Species in the genus Armillaria (fungi, basidiomycota) are well-known as saprophytes and pathogens on plants. Many of them cause white-rot root disease in diverse woody plants worldwide. Mitochondrial genomes (mitogenomes) are widely used in evolutionary and population studies, but despite the importance and wide distribution of Armillaria, the complete mitogenomes have not previously been reported for this genus. Meanwhile, the well-supported phylogeny of Armillaria species provides an excellent framework in which to study variation in mitogenomes and how they have evolved over time.

Results

Here we completely sequenced, assembled, and annotated the circular mitogenomes of four species: A. borealis, A. gallica, A. sinapina, and A. solidipes (116,443, 98,896, 103,563, and 122,167 bp, respectively). The variation in mitogenome size can be explained by variable numbers of mobile genetic elements, introns, and plasmid-related sequences. Most Armillaria introns contained open reading frames (ORFs) that are related to homing endonucleases of the LAGLIDADG and GIY-YIG families. Insertions of mobile elements were also evident as fragments of plasmid-related sequences in Armillaria mitogenomes. We also found several truncated gene duplications in all four mitogenomes.

Conclusions

Our study showed that fungal mitogenomes have a high degree of variation in size, gene content, and genomic organization even among closely related species of Armillara. We suggest that mobile genetic elements invading introns and intergenic sequences in the Armillaria mitogenomes have played a significant role in shaping their genome structure. The mitogenome changes we describe here are consistent with widely accepted phylogenetic relationships among the four species.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5732-z) contains supplementary material, which is available to authorized users.

Capsulimonas corticalis gen. nov., sp. nov., an aerobic capsulated bacterium, of a novel bacterial order, Capsulimonadales ord. nov., of the class Armatimonadia of the phylum Armatimonadetes

An aerobic bacterial strain designated AX-7^T was isolated from the trunk surface of a Japanese beech (Fagus crenata). Cells of strain AX-7^T were Gram-stain-negative, non-spore-forming, non-motile rods (1.0-1.2 µm in width and 1.2-3.0 µm in length) with peritrichous fimbriae. Cells were capsulated, and a number of them were surrounded by a thick slime layer. During growth, large aggregates formed, and the culture medium became viscous probably owing to exopolysaccharide release from the slime layer. The temperature range for growth was 10-37 °C, with an optimum at 30 °C. The pH range for growth was 5.0-7.0, with an optimum at pH 6.0. Strain AX-7^T used various sugars, including polysaccharides, and yeast extract as growth substrates. Strain AX-7^T contained menaquinones MK-9 and MK-10 as the respiratory quinones, and C16 : 1ω5c, C16 : 1ω11c, C16 : 0 and C14 : 0 as the major cellular fatty acids. Four unidentified phospholipids and 11 unidentified polar lipids constituted the polar lipids. The DNA G+C content was 61.0 mol%. The cell-wall peptidoglycan contained ll-diaminopimelic acid. Phylogenetic analysis based on 16S rRNA gene sequences revealed that strain AX-7^T belonged to the class Armatimonadia, its closest relative being Armatimonas rosea YO-36^T, with sequence similarity of 88.1%. Based on data from this polyphasic study, we propose that strain AX-7^T represents a new genus of a novel species within the novel order Capsulimonadales ord. nov. of the class Armatimonadia, for which the name Capsulimonas corticalis gen. nov., sp. nov. is proposed. The type strain of C. corticalis is AX-7^T (=DSM 105890^T=NBRC 113044^T).

Using insects to detect, monitor and predict the distribution of Xylella fastidiosa: a case study in Corsica

We sampled ca 2500 specimens of Philaenus spumarius (Hemiptera: Aphrophoridae) throughout Corsica without a priori knowledge on the presence of symptoms on plants. We screened 448 specimens for the presence of Xylella fastidiosa (Xf) using qPCR and a custom nested PCR. qPCR appeared versatile and under-estimated the prevalence of Xf. Nested PCR showed that Xf was present in all populations. Molecular results were validated by prediction on the distribution of Xf made from tests conducted on plants, which shows the pertinence of using vectors in risk assessment studies. Xf was detected in tenerals and adults. Thus, P. spumarius could acquire Xf from its host plant, mostly Cistus monspeliensis in Corsica, which may act as reservoir for the next season. This contrasts with other observations and suggests that management strategies may have to be adapted on a case-by-case basis. At least two genetic entities and several variants of Xf not yet identified on plants were present in the insects, which suggests ancient introductions of Xf and a probable underestimation of the current diversity of the strains present in Corsica. Interestingly 6% of the specimens carried two subspecies of Xf. Studies are required to better characterize the strains present in Corsica and to determine how the disease was introduced, spread and why no sign of a potential epidemic was detected earlier. This study shows that, when sensitive enough methods are implemented, spittlebugs (and more specifically P. spumarius for which species distribution modelling shows it could be a good sentinel for Europe) can be used to predict and better assess the exact distribution of Xf. Furthermore, Xf multiply only in their foregut and does not become circulative, which facilitates its detection.

Complex Models of Sequence Evolution Require Accurate Estimators as Exemplified with the Invariable Site Plus Gamma Model

The invariable site plus $\Gamma$ model (I$+\Gamma)$ is widely used to model rate heterogeneity among alignment sites in maximum likelihood and Bayesian phylogenetic analyses. The proof that the I$+$ continuous $\Gamma$ model is identifiable (model parameters can be inferred correctly given enough data) has increased the creditability of its application to phylogeny reconstruction. However, most phylogenetic software implement the I$+$ discrete $\Gamma$ model, whose identifiability is likely but unproven. How well the parameters of the I$+$ discrete $\Gamma$ model are estimated is still disputed. Especially the correlation between the fraction of invariable sites and the fractions of sites with a slow evolutionary rate is discussed as being problematic. We show that optimization heuristics as implemented in frequently used phylogenetic software (PhyML, RAxML, IQ-TREE, and MrBayes) cannot always reliably estimate the shape parameter, the proportion of invariable sites, and the tree length. Here, we propose an improved optimization heuristic that accurately estimates the three parameters. While research efforts mainly focus on tree search methods, our results signify the equal importance of verifying and developing effective estimation methods for complex models of sequence evolution.

A biologist's guide to Bayesian phylogenetic analysis

Bayesian methods have become very popular in molecular phylogenetics due to the availability of user-friendly software implementing sophisticated models of evolution. However, Bayesian phylogenetic models are complex, and analyses are often carried out using default settings, which may not be appropriate. Here, we summarize the major features of Bayesian phylogenetic inference and discuss Bayesian computation using Markov chain Monte Carlo (MCMC), the diagnosis of an MCMC run, and ways of summarising the MCMC sample. We discuss the specification of the prior, the choice of the substitution model, and partitioning of the data. Finally, we provide a list of common Bayesian phylogenetic software and provide recommendations as to their use.

Streptomyces aridus sp. nov., isolated from a high altitude Atacama Desert soil and emended description of Streptomyces noboritoensis Isono et al. 1957

A polyphasic study was undertaken to determine the taxonomic status of a Streptomyces strain which had been isolated from a high altitude Atacama Desert soil and shown to have bioactive properties. The strain, isolate H9^T, was found to have chemotaxonomic, cultural and morphological properties that place it in the genus Streptomyces. 16S rRNA gene sequence analyses showed that the isolate forms a distinct branch at the periphery of a well-delineated subclade in the Streptomyces 16S rRNA gene tree together with the type strains of Streptomyces crystallinus, Streptomyces melanogenes and Streptomyces noboritoensis. Multi-locus sequence analysis (MLSA) based on five house-keeping gene alleles showed that isolate H9^T is closely related to the latter two type strains and to Streptomyces polyantibioticus NRRL B-24448^T. The isolate was distinguished readily from the type strains of S. melanogenes, S. noboritoensis and S. polyantibioticus using a combination of phenotypic properties. Consequently, the isolate is considered to represent a new species of Streptomyces for which the name Streptomyces aridus sp. nov. is proposed; the type strain is H9^T (=NCIMB 14965^T=NRRL B65268^T). In addition, the MLSA and phenotypic data show that the S. melanogenes and S. noboritoensis type strains belong to a single species, it is proposed that S. melanogenes be recognised as a heterotypic synonym of S. noboritoensis for which an emended description is given.

High-throughput sequencing of multiple amplicons for barcoding and integrative taxonomy

Until now, the potential of NGS for the construction of barcode libraries or integrative taxonomy has been seldom realised. Here, we amplified (two-step PCR) and simultaneously sequenced (MiSeq) multiple markers from hundreds of fig wasp specimens. We also developed a workflow for quality control of the data. Illumina and Sanger sequences accumulated in the past years were compared. Interestingly, primers and PCR conditions used for the Sanger approach did not require optimisation to construct the MiSeq library. After quality controls, 87% of the species (76% of the specimens) had a valid MiSeq sequence for each marker. Importantly, major clusters did not always correspond to the targeted loci. Nine specimens exhibited two divergent sequences (up to 10%). In 95% of the species, MiSeq and Sanger sequences obtained from the same sampling were similar. For the remaining 5%, species were paraphyletic or the sequences clustered into divergent groups on the Sanger + MiSeq trees (>7%). These problematic cases may represent coding NUMTS or heteroplasms. Our results illustrate that Illumina approaches are not artefact-free and confirm that Sanger databases can contain non-target genes. This highlights the importance of quality controls, working with taxonomists and using multiple markers for DNA-taxonomy or species diversity assessment.

bModelTest: Bayesian phylogenetic site model averaging and model comparison

Background

Reconstructing phylogenies through Bayesian methods has many benefits, which include providing a mathematically sound framework, providing realistic estimates of uncertainty and being able to incorporate different sources of information based on formal principles. Bayesian phylogenetic analyses are popular for interpreting nucleotide sequence data, however for such studies one needs to specify a site model and associated substitution model. Often, the parameters of the site model is of no interest and an ad-hoc or additional likelihood based analysis is used to select a single site model.

Results

bModelTest allows for a Bayesian approach to inferring and marginalizing site models in a phylogenetic analysis. It is based on trans-dimensional Markov chain Monte Carlo (MCMC) proposals that allow switching between substitution models as well as estimating the posterior probability for gamma-distributed rate heterogeneity, a proportion of invariable sites and unequal base frequencies. The model can be used with the full set of time-reversible models on nucleotides, but we also introduce and demonstrate the use of two subsets of time-reversible substitution models.

Conclusion

With the new method the site model can be inferred (and marginalized) during the MCMC analysis and does not need to be pre-determined, as is now often the case in practice, by likelihood-based methods. The method is implemented in the bModelTest package of the popular BEAST 2 software, which is open source, licensed under the GNU Lesser General Public License and allows joint site model and tree inference under a wide range of models.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0890-6) contains supplementary material, which is available to authorized users.

Sulfated Polysaccharides in the Freshwater Green Macroalga Cladophora surera Not Linked to Salinity Adaptation

The presence of sulfated polysaccharides in cell walls of seaweeds is considered to be a consequence of the physiological adaptation to the high salinity of the marine environment. Recently, it was found that sulfated polysaccharides were present in certain freshwater Cladophora species and some vascular plants. Cladophora (Ulvophyceae, Chlorophyta) is one of the largest genera of green algae that are able to grow in both, seas and freshwater courses. Previous studies carried out on the water-soluble polysaccharides of the marine species C. falklandica established the presence of sulfated xylogalactoarabinans constituted by a backbone of 4-linked β-L-arabinopyranose units partially sulfated mainly on C3 and also on C2 with partial glycosylation, mostly on C2, with terminal β-D-xylopyranose or β-D-galactofuranose units. Besides, minor amounts of 3-, 6- and/or 3,6-linked β-D-galactan structures, with galactose in the pyranosic form were detected. In this work, the main water soluble cell wall polysaccharides from the freshwater alga Cladophora surera were characterized. It was found that this green alga biosynthesizes sulfated polysaccharides, with a structure similar to those found in marine species of this genus. Calibration of molecular clock with fossil data suggests that colonization of freshwater environments occurred during the Miocene by its ancestor. Therefore, the presence of sulfated polysaccharides in the freshwater green macroalga C. surera could be, in this case, an adaptation to transient desiccation and changes in ionic strength. Retention of sulfated polysaccharides at the cell walls may represent a snapshot of an evolutionary event, and, thus constitutes an excellent model for further studies on the mechanisms of sulfation on cell wall polysaccharides and environmental stress co-evolution.

Evolution of a Planktonic Foraminifer during Environmental Changes in the Tropical Oceans

Ecological adaptation to environmental changes is a strong driver of evolution, enabling speciation of pelagic plankton in the open ocean without the presence of effective physical barriers to gene flow. The tropical ocean environment, which plays an important role in shaping marine biodiversity, has drastically and frequently changed since the Pliocene. Nevertheless, the evolutionary history of tropical pelagic plankton has been poorly understood, as phylogeographic investigations are still in the developing state and paleontological approaches are insufficient to obtain a sequential record from the deep-sea sediments. The planktonic foraminifer Pulleniatina obliquiloculata is widely distributed in the tropical area throughout the world’s oceans, and its phylogeography is well established. It is thus one of the best candidates to examine how past environmental changes may have shifted the spatial distribution and affected the diversification of tropical pelagic plankton. Such an examination requires the divergence history of the planktonic foraminifers, yet the gene marker (partial small subunit (SSU) rDNA) previously used for phylogeographic studies was not powerful enough to achieve a high accuracy in estimating the divergence times. The present study focuses on improving the precision of divergence time estimates for the splits between sibling species (genetic types) of planktonic foraminifers by increasing the number of genes as well as the number of nucleotide bases used for molecular clock estimates. We have amplified the entire coding regions of two ribosomal RNA genes (SSU rDNA and large subunit (LSU) rDNA) of three genetic types of P. obliquiloculata and two closely related species for the first time and applied them to the Bayesian relaxed clock method. The comparison of the credible intervals of the four datasets consisting either of sequences of the partial SSU rDNA, the complete SSU rDNA, LSU rDNA, or a combination of both genes (SSU+LSU) clearly demonstrated that the two-gene dataset improved the accuracy of divergence time estimates. The P. obliquiloculata lineage diverged twice, first at the end of the Pliocene (3.1 Ma) and again in the middle Pleistocene (1.4 Ma). Both timings coincided with the environmental changes, which indirectly involved geographic separation of populations. The habitat of P. obliquiloculata was expanded toward the higher latitudinal zones during the stable warm periods and subsequently placed on the steep environmental gradients following the global cooling. Different environmental conditions in the stable warm tropics and unstable higher latitudes may have triggered ecological divergence among the populations, leading to adaptive differentiation and eventually speciation. A comprehensive analysis of divergence time estimates combined with phylogeography enabled us to reveal the evolutionary history of the pelagic plankton and to find the potential paleoenvironmental events, which could have changed their biogeography and ecology.

Modeling DNA methylation dynamics with approaches from phylogenetics

Motivation: Methylation of CpG dinucleotides is a prevalent epigenetic modification that is required for proper development in vertebrates. Genome-wide DNA methylation assays have become increasingly common, and this has enabled characterization of DNA methylation in distinct stages across differentiating cellular lineages. Changes in CpG methylation are essential to cellular differentiation; however, current methods for modeling methylation dynamics do not account for the dependency structure between precursor and dependent cell types.

Results: We developed a continuous-time Markov chain approach, based on the observation that changes in methylation state over tissue differentiation can be modeled similarly to DNA nucleotide changes over evolutionary time. This model explicitly takes precursor to descendant relationships into account and enables inference of CpG methylation dynamics. To illustrate our method, we analyzed a high-resolution methylation map of the differentiation of mouse stem cells into several blood cell types. Our model can successfully infer unobserved CpG methylation states from observations at the same sites in related cell types (90% correct), and this approach more accurately reconstructs missing data than imputation based on neighboring CpGs (84% correct). Additionally, the single CpG resolution of our methylation dynamics estimates enabled us to show that DNA sequence context of CpG sites is informative about methylation dynamics across tissue differentiation. Finally, we identified genomic regions with clusters of highly dynamic CpGs and present a likely functional example. Our work establishes a framework for inference and modeling that is well suited to DNA methylation data, and our success suggests that other methods for analyzing DNA nucleotide substitutions will also translate to the modeling of epigenetic phenomena.

Availability and implementation: Source code is available at www.kostkalab.net/software.

Contact: tony.capra@vanderbilt.edu or kostka@pitt.edu

Phylogeography and systematics of the westernmost Italian Dolichopoda species (Orthoptera, Rhaphidophoridae)

The genus Dolichopoda (Orthoptera; Rhaphidopohoridae) is present in Italy with 9 species distributed from northwestern Italy (Piedmont and Liguria) to the southernmost Apennines (Calabria), occurring also in the Tyrrhenian coastal areas and in Sardinia. Three morphologically very close taxa have been described in Piedmont and Liguria, i.e., D. ligustica ligustica, D. ligustica septentrionalis and D. azami azami. To investigate the delimitation of the northwestern species of Dolichopoda, we performed both morphological and molecular analyses. Morphological analysis was carried out by considering diagnostic characters generally used to distinguish different taxa, as the shape of epiphallus in males and the subgenital fig in females. Molecular analysis was performed by sequencing three mitochondrial genes, 12S rRNA, 16S rRNA, partially sequenced and the entire gene of COI. Results from both morphological and molecular analyses highlighted a very homogeneous group of populations, although genetically structured. Three haplogroups geographically distributed could be distinguished and based on these results we suggest a new taxonomic arrangement. All populations, due to the priority of description, should be assigned to D. azami azami Saulcy, 1893 and to preserve the names ligustica and septentrionalis, corresponding to different genetic haplogroups, we assign them to D. azami ligustica stat. n. Baccetti & Capra, 1959 and to D. azami septentrionalis stat. n. Baccetti & Capra, 1959.

πBUSS: a parallel BEAST/BEAGLE utility for sequence simulation under complex evolutionary scenarios

Background

Simulated nucleotide or amino acid sequences are frequently used to assess the performance of phylogenetic reconstruction methods. BEAST, a Bayesian statistical framework that focuses on reconstructing time-calibrated molecular evolutionary processes, supports a wide array of evolutionary models, but lacked matching machinery for simulation of character evolution along phylogenies.

Results

We present a flexible Monte Carlo simulation tool, called πBUSS, that employs the BEAGLE high performance library for phylogenetic computations to rapidly generate large sequence alignments under complex evolutionary models. πBUSS sports a user-friendly graphical user interface (GUI) that allows combining a rich array of models across an arbitrary number of partitions. A command-line interface mirrors the options available through the GUI and facilitates scripting in large-scale simulation studies. πBUSS may serve as an easy-to-use, standard sequence simulation tool, but the available models and data types are particularly useful to assess the performance of complex BEAST inferences. The connection with BEAST is further strengthened through the use of a common extensible markup language (XML), allowing to specify also more advanced evolutionary models. To support simulation under the latter, as well as to support simulation and analysis in a single run, we also add the πBUSS core simulation routine to the list of BEAST XML parsers.

Conclusions

πBUSS offers a unique combination of flexibility and ease-of-use for sequence simulation under realistic evolutionary scenarios. Through different interfaces, πBUSS supports simulation studies ranging from modest endeavors for illustrative purposes to complex and large-scale assessments of evolutionary inference procedures. Applications are not restricted to the BEAST framework, or even time-measured evolutionary histories, and πBUSS can be connected to various other programs using standard input and output format.