MetaPIGA v2.0: maximum likelihood large phylogeny estimation using the metapopulation genetic algorithm and other stochastic heuristics

Cloning, characterization, and evolutionary patterns of KCNQ4 genes in anurans

Acoustic communication plays important roles in the survival and reproduction of anurans. The perception and discrimination of conspecific sound signals of anurans were always affected by masking background noise. Previous studies suggested that some frogs evolved the high-frequency hearing to minimize the low-frequency noise. However, the molecular mechanisms underlying the high-frequency hearing in anurans have not been well explored. Here, we cloned and obtained the coding regions of a high-frequency hearing-related gene (KCNQ4) from 11 representative anuran species and compared them with orthologous sequences from other four anurans. The sequence characteristics and evolutionary analyses suggested the highly conservation of the KCNQ4 gene in anurans, which supported their functional importance. Branch-specific analysis showed that KCNQ4 genes were under different evolutionary forces in anurans and most anuran lineages showed a generally strong purifying selection. Intriguingly, one significantly positively selected site was identified in the anuran KCNQ4 gene based on FEL model. Positive selection was also found along the common ancestor of Ranidae and Rhacophoridae as well as the ancestral O. tianmuii based on the branch-site analysis, and the positively selected sites identified were involved in or near the N-terminal ion transport and the potassium ion channel functional domain of the KCNQ4 genes. The present study revealed valuable information regarding the KCNQ4 genes in anurans and provided some new insights for the underpinnings of the high-frequency hearing in frogs.

A Guide to Phylogenomic Inference

Phylogenomics aims at reconstructing the evolutionary histories of organisms taking into account whole genomes or large fractions of genomes. Phylogenomics has significant applications in fields such as evolutionary biology, systematics, comparative genomics, and conservation genetics, providing valuable insights into the origins and relationships of species and contributing to our understanding of biological diversity and evolution. This chapter surveys phylogenetic concepts and methods aimed at both gene tree and species tree reconstruction while also addressing common pitfalls, providing references to relevant computer programs. A practical phylogenomic analysis example including bacterial genomes is presented at the end of the chapter.

An evolution strategy approach for the balanced minimum evolution problem

MOTIVATION

The Balanced Minimum Evolution (BME) is a powerful distance based phylogenetic estimation model introduced by Desper and Gascuel and nowadays implemented in popular tools for phylogenetic analyses. It was proven to be computationally less demanding than more sophisticated estimation methods, e.g. maximum likelihood or Bayesian inference while preserving the statistical consistency and the ability to run with almost any kind of data for which a dissimilarity measure is available. BME can be stated in terms of a nonlinear non-convex combinatorial optimization problem, usually referred to as the Balanced Minimum Evolution Problem (BMEP). Currently, the state-of-the-art among approximate methods for the BMEP is represented by FastME (version 2.0), a software which implements several deterministic phylogenetic construction heuristics combined with a local search on specific neighbourhoods derived by classical topological tree rearrangements. These combinations, however, may not guarantee convergence to close-to-optimal solutions to the problem due to the lack of solution space exploration, a phenomenon which is exacerbated when tackling molecular datasets characterized by a large number of taxa.

RESULTS

To overcome such convergence issues, in this article, we propose a novel metaheuristic, named PhyloES, which exploits the combination of an exploration phase based on Evolution Strategies, a special type of evolutionary algorithm, with a refinement phase based on two local search algorithms. Extensive computational experiments show that PhyloES consistently outperforms FastME, especially when tackling larger datasets, providing solutions characterized by a shorter tree length but also significantly different from the topological perspective.

AVAILABILITY AND IMPLEMENTATION

The software and the data are available at https://github.com/andygaspar/PHYLOES.

Mitochondrial Lineage Diversity and Phylogeography of Daphnia (Daphnia) (Crustacea: Cladocera) in North-East Russia

The variability of the 12S gene fragment of the mtDNA for taxa belonging to subgenus Daphnia (Daphnia) O.F. Müller, 1776 (Crustacea: Cladocera) in NE Russia is studied, and their phylogenetic analysis performed. We identified (based both on morphological and molecular data) nine species belonging to four species complexes, namely: (A) D. longispina s.l.: (1) D. longispina O.F. Müller, 1776; (2) D. dentifera Forbes, 1893; (3) D. galeata Sars, 1864; (4) D. umbra Taylor, Hebert et Colbourne, 1996; (B) D. cristata s.l.: (5) D. cristata Sars, 1862; (6) D. longiremis Sars, 1862; (C) D. curvirostris s.l.: (7) D. curvirostris Eylmann, 1887; (D) D. pulex s.l.: (8) D. pulex Leydig, 1860; (9) D. middendorffiana Fischer, 1851. Rare arcto-mountainous taxon D. umbra was found in the mountains of the Sakha (Yakutia) Republic for the first time. Species diversity in NE Asia is relatively low, and the most revealed taxa are trans-Beringian. We also performed a phylogeographic analysis of D. dentifera and D. pulex s.l., the two most common species in NE Russia. Our new data allow us to assume that the daphniids of NE Asia have undergone various evolutionary scenarios during the Pleistocene period: survival is within some local refugia, and re-colonization from these areas and from North America through the Beringian land bridge, etc. We agree with previous authors who revealed that the patterns in the studied species groups are relatively recent (of Late Pleistocene or even Holocene age), although the main phylogenetic daphniid lineages (mainly congruent with the biological species) are very old. Our results provide convincing evidence for the hypothesis that NE Russia is a very important source of modern haplotypic diversity for the cladocerans.

OUP accepted manuscript

Motivation

In recent years, full-genome sequences have become increasingly available and as a result many modern phylogenetic analyses are based on very long sequences, often with over 100 000 sites. Phylogenetic reconstructions of large-scale alignments are challenging for likelihood-based phylogenetic inference programs and usually require using a powerful computer cluster. Current tools for alignment trimming prior to phylogenetic analysis do not promise a significant reduction in the alignment size and are claimed to have a negative effect on the accuracy of the obtained tree.

Results

Here, we propose an artificial-intelligence-based approach, which provides means to select the optimal subset of sites and a formula by which one can compute the log-likelihood of the entire data based on this subset. Our approach is based on training a regularized Lasso-regression model that optimizes the log-likelihood prediction accuracy while putting a constraint on the number of sites used for the approximation. We show that computing the likelihood based on 5% of the sites already provides accurate approximation of the tree likelihood based on the entire data. Furthermore, we show that using this Lasso-based approximation during a tree search decreased running-time substantially while retaining the same tree-search performance.

Availability and implementation

The code was implemented in Python version 3.8 and is available through GitHub (https://github.com/noaeker/lasso_positions_sampling). The datasets used in this paper were retrieved from Zhou et al. (2018) as described in section 3.

Supplementary information

Supplementary data are available at Bioinformatics online.

How to choose the right real-time RT-PCR primer sets for the SARS-CoV-2 genome detection?

OBJECTIVES

The SARS-CoV-2 pandemic has created an unprecedented need for rapid large-scale diagnostic testing to prompt clinical and public health interventions. Currently, several quantitative reverse-transcription polymerase chain reaction (RT-qPCR) assays recommended by the World Health Organization are being used by clinical and public health laboratories and typically target regions of the RNA-dependent RNA polymerase (RdRp), envelope (E) and nucleocapsid (N) coding region. However, it is currently unclear if results from different tests are comparable. This study aimed to clarify the clinical performances of the primer/probe sets designed by US CDC and Charité/Berlin to help clinical laboratories in assay selection for SARS-CoV-2 routine detection.

METHODS

We compared the clinical performances of the recommended primer/probe sets using one hundred nasopharyngeal swab specimens from patients who were clinically diagnosed with COVID-19. An additional 30 "pre-intervention screening" samples from patients who were not suspected of COVID-19 were also included in this study. We also performed sequence alignment between 31064 European SARS-CoV-2 and variants of concern genomes and the recommended primer/probe sets.

RESULTS

The present study demonstrates substantial differences in SARS-CoV-2 RNA detection sensitivity among the primer/probe sets recommended by the World Health Organization especially for low-level viral loads. The alignment of thousands of SARS-CoV-2 sequences reveals that the genetic diversity remains relatively low at the primer/probe binding sites. However, multiple nucleotide mismatches might contribute to false negatives.

CONCLUSION

An understanding of the limitations depending on the targeted genes and primer/probe sets may influence the selection of molecular detection assays by clinical laboratories.

Harnessing machine learning to guide phylogenetic-tree search algorithms

Inferring a phylogenetic tree is a fundamental challenge in evolutionary studies. Current paradigms for phylogenetic tree reconstruction rely on performing costly likelihood optimizations. With the aim of making tree inference feasible for problems involving more than a handful of sequences, inference under the maximum-likelihood paradigm integrates heuristic approaches to evaluate only a subset of all potential trees. Consequently, existing methods suffer from the known tradeoff between accuracy and running time. In this proof-of-concept study, we train a machine-learning algorithm over an extensive cohort of empirical data to predict the neighboring trees that increase the likelihood, without actually computing their likelihood. This provides means to safely discard a large set of the search space, thus potentially accelerating heuristic tree searches without losing accuracy. Our analyses suggest that machine learning can guide tree-search methodologies towards the most promising candidate trees.

Bioluminescence in lanternsharks: Insight from hormone receptor localization

As part of the study of their bioluminescence, the deep-sea lanternshark Etmopterus spinax and Etmopterus molleri (Chondrichthyes, Etmopteridae) received growing interest over the past ten years. These mesopelagic sharks produce light thanks to a finely tuned hormonal control involving melatonin, adrenocorticotropic hormone and α-melanocyte-stimulating hormone. Receptors of these hormones, respectively the melatonin receptors and the melanocortin receptors, are all members of the G-protein coupled receptor family i.e. coupled with specific G proteins involved in the preliminary steps of their transduction pathways. The present study highlights the specific localization of the hormonal receptors, as well as of their associated G-proteins within the light organs, the so-called photophores, in E. spinax and E. molleri through immunohistofluorescence technic. Our results allow gaining insight into the molecular actors and mechanisms involved in the control of the light emission in Etmopterid sharks.

Distance-based phylogenetic inference from typing data: a unifying view

Typing methods are widely used in the surveillance of infectious diseases, outbreaks investigation and studies of the natural history of an infection. Moreover, their use is becoming standard, in particular with the introduction of high-throughput sequencing. On the other hand, the data being generated are massive and many algorithms have been proposed for a phylogenetic analysis of typing data, addressing both correctness and scalability issues. Most of the distance-based algorithms for inferring phylogenetic trees follow the closest pair joining scheme. This is one of the approaches used in hierarchical clustering. Moreover, although phylogenetic inference algorithms may seem rather different, the main difference among them resides on how one defines cluster proximity and on which optimization criterion is used. Both cluster proximity and optimization criteria rely often on a model of evolution. In this work, we review, and we provide a unified view of these algorithms. This is an important step not only to better understand such algorithms but also to identify possible computational bottlenecks and improvements, important to deal with large data sets.

Analysing ambiguities in trypanosomatids taxonomy by barcoding

BACKGROUND

Biodiversity screens and phylogenetic studies are dependent on reliable DNA sequences in public databases. Biological collections possess vouchered specimens with a traceable history. Therefore, DNA sequencing of samples available at institutional collections can greatly contribute to taxonomy, and studies on evolution and biodiversity.

METHODS

We sequenced part of the glycosomal glyceraldehyde phosphate dehydrogenase (gGAPDH) and the SSU rRNA (V7/V8) genes from 102 trypanosomatid cultures, which are available on request at www.colprot.fiocruz.br.

OBJECTIVE

The main objective of this work was to use phylogenetic inferences, using the obtained DNA sequences and those from representatives of all Trypanosomatidae genera, to generate phylogenetic trees that can simplify new isolates screenings.

FINDINGS

A DNA sequence is provided for the first time for several isolates, the phylogenetic analysis allowed the classification or reclassification of several specimens, identification of candidates for new genera and species, as well as the taxonomic validation of several deposits.

MAIN CONCLUSIONS

This survey aimed at presenting a list of validated species and their associated DNA sequences combined with a short historical overview of each isolate, which can support taxonomic and biodiversity research and promote culture collections.

Optimizing taxon addition order and branch lengths in the construction of phylogenetic trees using maximum likelihood

Taxon addition order and branch lengths are optimized by genetic algorithms (GAS) within the fastDNAml algorithm for constructing phylogenetic trees of high likelihood. Results suggest that optimizing the order in which taxa are added improves the likelihood of the resulting trees.

Morelli K, L. C. Santos H, Brouillard P, Grellier P, Gala JL, Vikkula M. First Draft Genome of the Trypanosomatid Herpetomonas muscarum ingenoplastis through MinION Oxford Nanopore Technology and Illumina Sequencing

Here, we present first draft genome sequence of the trypanosomatid Herpetomonas muscarum ingenoplastis. This parasite was isolated repeatedly in the black blowfly, Phormia regina, and it forms a phylogenetically distinct clade in the Trypanosomatidae family.

A comprehensive approach uncovers hidden diversity in freshwater mussels (Bivalvia: Unionidae) with the description of a novel species

Major geological processes have shaped biogeographical patterns of riverine biota. The Edwards Plateau of central Texas, USA, exhibits unique aquatic communities and endemism, including several species of freshwater mussels. Lampsilis bracteata (Gould, 1855) is endemic to the Edwards Plateau region; however, its phylogenetic relationship with other species in the Gulf coastal rivers and Mississippi River basin is unknown. We evaluated phylogenetic relationships, shell morphologies and soft anatomy characters of L. bracteata and a closely related congener, Lampsilis hydiana (Lea, 1838) throughout their ranges. Our results showed the presence of an undescribed species: Lampsilis bergmanni sp.n. Lampsilis bracteata and L. bergmanni sp.n. share similar shell morphologies and soft anatomy characters; however, they are genetically distinct. Geological processes, such as faulting and sea-level changes during the Miocene to Pliocene, are likely to have facilitated diversification of Lampsilis species, resulting in isolation of L. bracteata on the Edwards Plateau and diversification between L. bergmanni sp.n. and L. hydiana. We conclude that L. bracteata range is restricted to the Colorado River basin, whereas L. bergmanni sp.n. occurs only in upstream reaches of the Guadalupe River basin. Conservation actions are warranted for both species due to their restricted distributions and potential anthropogenic threats.

Reversals in complex traits uncovered as reticulation events: Lessons from the evolution of parity-mode, chromosome morphology, and maternal resource transfer

Complex traits include, among many others, the evolution of eyes, wings, body forms, reproductive modes, human intelligence, social behavior, diseases, and chromosome morphology. Dollo's law states that the evolution of complex traits is irreversible. However, potential exceptions have been proposed. Here, we investigated whether reticulation, a simple and elegant means by which complex characters may be reacquired, could account for suggested reversals in the evolution of complex characters using two datasets with sufficient genetic coverage and a total of five potential reversals. Our analyses uncovered a potential reversal in the evolution of parity mode and a potential reversal in the evolution of placentotrophy of fish (Cyprinodontiformes) as reticulation events. Moreover, in a reptile that exhibits a potential reversal from viviparity to oviparity (Zootoca vivipara), reticulation provided the most parsimonious explanation for sex chromosome evolution. Therefore, three of the five studied potential reversals were unraveled as reticulation events. This constitutes the first evidence that accounting for reticulation can fundamentally influence the interpretation of the evolution of complex traits, that testing for reticulation is crucial for obtaining robust phylogenies, and that complex ancestral characters may be reacquired through hybridization with a lineage that still exhibits the trait. Hybridization, rather than reappearance of ancestral traits by means of small evolutionary steps, may thus account for suggested exceptions to Dollo's law. Consequently, ruling out reticulation is required to claim the evolutionary reversal of complex characters and potential exceptions to Dollo's rule.

Interspecies comparison of sea star adhesive proteins

Sea stars use adhesive secretions to attach their numerous tube feet strongly and temporarily to diverse surfaces. After detachment of the tube feet, the adhesive material stays bound to the substrate as so-called 'footprints'. In the common sea star species Asterias rubens, the adhesive material has been studied extensively and the first sea star footprint protein (Sfp1) has been characterized. We identified Sfp1-like sequences in 17 additional sea star species, representing different taxa and tube foot morphologies, and analysed the evolutionary conservation of this protein. In A. rubens, we confirmed the expression of 34 footprint proteins in the tube foot adhesive epidermis, with 22 being exclusively expressed in secretory cells of the adhesive epidermis and 12 showing an additional expression in the stem epidermis. The sequences were used for BLAST searches in seven asteroid transcriptomes providing a first insight in the conservation of footprint proteins among sea stars. Our results highlighted a high conservation of the large proteins making up the structural core of the footprints, whereas smaller, potential surface-binding proteins might be more variable among sea star species. This article is part of the theme issue 'Transdisciplinary approaches to the study of adhesion and adhesives in biological systems'.

Unexpected endemism in the Daphnia longispina complex (Crustacea: Cladocera) in Southern Siberia

The biological significance of regional cladoceran morphotypes in the montane regions of the central Palearctic remains poorly understood. In the Holarctic Daphnia longispina complex (Cladocera: Daphniidae), several variants, lineages and species have been proposed as endemic for Southern Siberia. Daphnia turbinata Sars, for example, named after its unusual head shape, is known only from Southern Siberia. Here we sequence DNA of Daphnia from three mitochondrial genes (12S rRNA, 16S rRNA, and NADH dehydrogenase subunit 2, ND2) from 57 localities in Russia and Mongolia (the majority being from Southern Siberia) and place them in evolutionary context with existing data. Our aim was to examine regional endemism of the Daphnia longispina complex in Southern Siberian; to improve the phylogenetic understanding with improved taxonomic and regional sampling, and to better understand the influence of Pleistocene glaciation on the biogeography of these lineages. At least three lineages showed genetic evidence for endemism in Southern Siberia. There was strong support for D. turbinata as a sister lineage to to D. longispina/D. dentifera. Another endemic, Siberian D. cf. longispina, is a sister group to the longispina group in general. Within D. longispina s. str. there was an endemic Siberian clade with a western range boundary near the Yenisei River Basin. Gene flow estimates among populations (based on F_ST values) were very low for clades of D. longispina on a regional (the original 12S dataset), and on a pan-Eurasian (the extended 12S dataset) scale. Negative values of Fu’s F_S and Tajima’s D tests prevailed for the species examined with significant values found for two D. longispina clades, D. dentifera, D. galeata and D. cristata. Our results support the notion that Southern Siberia is an important biogeographic region for cladocerans as it contained unexpected diversity of endemics (such as D. turbinata, D. cf. longispina and lineages of D. umbra and D. longsipina s.str.) and from being the geographic meeting place of expanding postglacial lineages from eastern and western refugia.

The AB loop of oncostatin M (OSM) determines species-specific signaling in humans and mice

The pleiotropic interleukin-6 (IL-6)-type cytokine oncostatin M (OSM) signals in multiple cell types, affecting processes such as cell differentiation, hematopoiesis, and inflammation. In humans, OSM exerts its effects through activation of either of two different heterodimeric receptor complexes, formed by glycoprotein 130 (gp130) and either OSM receptor (OSMR) or leukemia inhibitory factor receptor (LIFR). In contrast, the mouse OSM orthologue acts mainly through dimers containing OSMR and gp130 and shows limited activity through mouse LIFR. Despite their structural similarity, neither human nor mouse OSM signal through the other species' OSMR. The molecular basis for such species-specific signaling, however, remains poorly understood. To identify key molecular features of OSM that determine receptor activation in humans and mice, we generated chimeric mouse-human cytokines. Replacing regions within binding site III of murine OSM with the human equivalents showed that the cytokine's AB loop was critical for receptor selection. Substitutions of individual amino acids within this region demonstrated that residues Asn-37, Thr-40, and Asp-42 of the murine cytokine were responsible for limited LIFR activation and absence of human OSMR/LIFR signaling. In human OSM, Lys-44 appeared to be the main residue preventing mouse OSMR activation. Our data reveal that individual amino acids within the AB loop of OSM determine species-specific activities. These mutations might reflect a key step in the evolutionary process of this cytokine, in which receptor promiscuity gives way to ligand-receptor specialization.

Selective human inhibitors of ATR and ATM render Leishmania major promastigotes sensitive to oxidative damage

All cellular processes, including those involved in normal cell metabolism to those responsible for cell proliferation or death, are finely controlled by cell signaling pathways, whose core proteins constitute the family of phosphatidylinositol 3-kinase-related kinases (PIKKs). Ataxia Telangiectasia Mutated (ATM) and Ataxia Telangiectasia and Rad3 related (ATR) are two important PIKK proteins that act in response to DNA damage, phosphorylating a large number of proteins to exert control over genomic integrity. The genus Leishmania belongs to a group of early divergent eukaryotes in evolution and has a highly plastic genome, probably owing to the existence of signaling pathways designed to maintain genomic integrity. The objective of this study was to evaluate the use of specific human inhibitors of ATR and ATM in Leishmania major. Bioinformatic analyses revealed the existence of the putative PIKK genes ATR and ATM, in addition to mTOR and DNA-PKcs in Leishmania spp. Moreover, it was possible to suggest that the inhibitors VE-821 and KU-55933 have binding affinity for the catalytic sites of putative L. major ATR and ATM, respectively. Promastigotes of L. major exposed to these inhibitors show slight growth impairment and minor changes in cell cycle and morphology. It is noteworthy that treatment of promastigotes with inhibitors VE-821 and KU-55933 enhanced the oxidative damage caused by hydrogen peroxide. These inhibitors could significantly reduce the number of surviving L. major cells following H₂O₂ exposure whilst also decreasing their evaluated IC₅₀ to H₂O₂ to less than half of that observed for non-treated cells. These results suggest that the use of specific inhibitors of ATR and ATM in Leishmania interferes in the signaling pathways of this parasite, which can impair its tolerance to DNA damage and affect its genome integrity. ATR and ATM could constitute novel targets for drug development and/or repositioning for treatment of leishmaniases.

Genes for Membrane Transport Proteins: Not So Rare in Viruses

Some viruses have genes encoding proteins with membrane transport functions. It is unknown if these types of proteins are rare or are common in viruses. In particular, the evolutionary origin of some of the viral genes is obscure, where other viral proteins have homologs in prokaryotic and eukaryotic organisms. We searched virus genomes in databases looking for transmembrane proteins with possible transport function. This effort led to the detection of 18 different types of putative membrane transport proteins indicating that they are not a rarity in viral genomes. The most abundant proteins are K⁺ channels. Their predicted structures vary between different viruses. With a few exceptions, the viral proteins differed significantly from homologs in their current hosts. In some cases the data provide evidence for a recent gene transfer between host and virus, but in other cases the evidence indicates a more complex evolutionary history.

Isolation drives increased diversification rates in freshwater amphipods

Vicariance and dispersal events affect current biodiversity patterns in desert springs. Whether major diversification events are due to environmental changes leading to radiation or due to isolation resulting in relict species is largely unknown. We seek to understand whether the Gammarus pecos species complex underwent major diversification events due to environmental changes in the area leading either to radiation into novel habitats, or formation of relicts due to isolation. Specifically, we tested the hypothesis that Gammarus in the northern Chihuahuan Desert of New Mexico and Texas, USA are descendants of an ancient marine lineage now containing multiple undescribed species. We sequenced a nuclear (28S) and two mitochondrial (16S, COI) genes from gammarid amphipods representing 16 desert springs in the northern Chihuahuan Desert. We estimated phylogenetic relationships, divergence times, and diversification rates of the Gammarus pecos complex. Our results revealed that the region contained two evolutionarily independent lineages: a younger Freshwater Lineage that shared a most-recent-common-ancestor with an older Saline Lineage ∼66.3  MYA (95.6-42.4  MYA). Each spring system generally formed a monophyletic clade based on the concatenated dataset. Freshwater Lineage diversification rates were 2.0-9.8 times higher than rates of the Saline Lineage. A series of post-Cretaceous colonizations by ancestral Gammarus taxa was likely followed by isolation. Paleo-geological, hydrological, and climatic events in the Neogene-to-Quaternary periods (23.03  MYA - present) in western North America promoted allopatric speciation of both lineages. We suggest that Saline Lineage populations include two undescribed Gammarus species, while the Freshwater Lineage shows repetition of fine-scale genetic structure in all major clades suggesting incipient speciation. Such ongoing speciation suggests that this region will continue to be a biodiversity hotspot for amphipods and other freshwater taxa.

The cereal pathogen Fusarium pseudograminearum produces a new class of active cytokinins during infection

The fungal pathogen Fusarium pseudograminearum causes important diseases of wheat and barley. During a survey of secondary metabolites produced by this fungus, a novel class of cytokinins, herein termed Fusarium cytokinins, was discovered. Cytokinins are known for their growth-promoting and anti-senescence activities, and the production of a cytokinin mimic by what was once considered as a necrotrophic pathogen that promotes cell death and senescence challenges the simple view that this pathogen invades its hosts by employing a barrage of lytic enzymes and toxins. Through genome mining, a gene cluster in the F. pseudograminearum genome for the production of Fusarium cytokinins was identified and the biosynthetic pathway was established using gene knockouts. The Fusarium cytokinins could activate plant cytokinin signalling, demonstrating their genuine hormone mimicry. In planta analysis of the transcriptional response to one Fusarium cytokinin suggests extensive reprogramming of the host environment by these molecules, possibly through crosstalk with defence hormone signalling pathways.

New Mycobacterium tuberculosis Complex Sublineage, Brazzaville, Congo

Tuberculosis is a leading cause of illness and death in Congo. No data are available about the population structure and transmission dynamics of the Mycobacterium tuberculosis complex strains prevalent in this central Africa country. On the basis of single-nucleotide polymorphisms detected by whole-genome sequencing, we phylogenetically characterized 74 MTBC isolates from Brazzaville, the capital of Congo. The diversity of the study population was high; most strains belonged to the Euro-American lineage, which split into Latin American Mediterranean, Uganda I, Uganda II, Haarlem, X type, and a new dominant sublineage named Congo type (n = 26). Thirty strains were grouped in 5 clusters (each within 12 single-nucleotide polymorphisms), from which 23 belonged to the Congo type. High cluster rates and low genomic diversity indicate recent emergence and transmission of the Congo type, a new Euro-American sublineage of MTBC.

The Pleurobemini (Bivalvia : Unionida) revisited: molecular species delineation using a mitochondrial DNA gene reveals multiple conspecifics and undescribed species

The Pleurobemini (Bivalvia: Unionida) represent approximately one-third of freshwater mussel diversity in North America. Species identification within this group is challenging due to morphological convergence and phenotypic plasticity. Accurate species identification, including characterisation of currently unrecognised taxa, is required to develop effective conservation strategies because many species in the group are imperiled. We examined 575 cox1 sequences from 110 currently recognised species (including 13 Fusconaia and 21 Pleurobema species) to understand phylogenetic relationships among pleurobemine species (mainly Fusconaia and Pleurobema) and to delineate species boundaries. The results of phylogenetic analyses showed no geographic structure within widespread species and illustrated a close relationship between Elliptio lanceolata and Parvaspina collina. Constraint tests supported monophyly of the genera Fusconaia and Pleurobema, including the subgenus P. (Sintoxia). Furthermore, results revealed multiple conspecifics, including P. hanleyianum and P. troschelianum, P. chattanoogaense and P. decisum, P. clava and P. oviforme, P. rubrum and P. sintoxia, F. askewi and F. lananensis, and F. cerina and F. flava. Species delimitation analyses identified three currently unrecognised taxa (two in Fusconaia and one in Pleurobema). Further investigation using additional genetic markers and other lines of evidence (e.g. morphology, life history, ecology) are necessary before any taxonomic changes are formalised.

Using MOEA with Redistribution and Consensus Branches to Infer Phylogenies

In recent years, to infer phylogenies, which are NP-hard problems, more and more research has focused on using metaheuristics. Maximum Parsimony and Maximum Likelihood are two effective ways to conduct inference. Based on these methods, which can also be considered as the optimal criteria for phylogenies, various kinds of multi-objective metaheuristics have been used to reconstruct phylogenies. However, combining these two time-consuming methods results in those multi-objective metaheuristics being slower than a single objective. Therefore, we propose a novel, multi-objective optimization algorithm, MOEA-RC, to accelerate the processes of rebuilding phylogenies using structural information of elites in current populations. We compare MOEA-RC with two representative multi-objective algorithms, MOEA/D and NAGA-II, and a non-consensus version of MOEA-RC on three real-world datasets. The result is, within a given number of iterations, MOEA-RC achieves better solutions than the other algorithms.

Culture and Next-generation sequencing-based drug susceptibility testing unveil high levels of drug-resistant-TB in Djibouti: results from the first national survey

Djibouti is a small country in the Horn of Africa with a high TB incidence (378/100,000 in 2015). Multidrug-resistant TB (MDR-TB) and resistance to second-line agents have been previously identified in the country but the extent of the problem has yet to be quantified. A national survey was conducted to estimate the proportion of MDR-TB among a representative sample of TB patients. Sputum was tested using XpertMTB/RIF and samples positive for MTB and resistant to rifampicin underwent first line phenotypic susceptibility testing. The TB supranational reference laboratory in Milan, Italy, undertook external quality assurance, genotypic testing based on whole genome and targeted-deep sequencing and phylogenetic studies. 301 new and 66 previously treated TB cases were enrolled. MDR-TB was detected in 34 patients: 4.7% of new and 31% of previously treated cases. Resistance to pyrazinamide, aminoglycosides and capreomycin was detected in 68%, 18% and 29% of MDR-TB strains respectively, while resistance to fluoroquinolones was not detected. Cluster analysis identified transmission of MDR-TB as a critical factor fostering drug resistance in the country. Levels of MDR-TB in Djibouti are among the highest on the African continent. High prevalence of resistance to pyrazinamide and second-line injectable agents have important implications for treatment regimens.

Convergent origination of a Drosophila-like dosage compensation mechanism in a reptile lineage

Sex chromosomes differentiated from different ancestral autosomes in various vertebrate lineages. Here, we trace the functional evolution of the XY Chromosomes of the green anole lizard (Anolis carolinensis), on the basis of extensive high-throughput genome, transcriptome and histone modification sequencing data and revisit dosage compensation evolution in representative mammals and birds with substantial new expression data. Our analyses show that Anolis sex chromosomes represent an ancient XY system that originated at least ≈160 million years ago in the ancestor of Iguania lizards, shortly after the separation from the snake lineage. The age of this system approximately coincides with the ages of the avian and two mammalian sex chromosomes systems. To compensate for the almost complete Y Chromosome degeneration, X-linked genes have become twofold up-regulated, restoring ancestral expression levels. The highly efficient dosage compensation mechanism of Anolis represents the only vertebrate case identified so far to fully support Ohno's original dosage compensation hypothesis. Further analyses reveal that X up-regulation occurs only in males and is mediated by a male-specific chromatin machinery that leads to global hyperacetylation of histone H4 at lysine 16 specifically on the X Chromosome. The green anole dosage compensation mechanism is highly reminiscent of that of the fruit fly, Drosophila melanogaster. Altogether, our work unveils the convergent emergence of a Drosophila-like dosage compensation mechanism in an ancient reptilian sex chromosome system and highlights that the evolutionary pressures imposed by sex chromosome dosage reductions in different amniotes were resolved in fundamentally different ways.

Real-time imaging of the growth-inhibitory effect of JS399-19 on Fusarium

Real-time imaging was used to study the effects of a novel Fusarium-specific cyanoacrylate fungicide (JS399-19) on growth and morphology of four Fusarium sp. This fungicide targets the motor domain of type I myosin. Fusarium graminearum PH-1, Fusarium solani f. sp. pisi 77-13-4, Fusarium avenaceum IBT8464, and Fusarium avenaceum 05001, which has a K216Q amino-acid substitution at the resistance-implicated site in its myosin type I motor domain, were analyzed. Real-time imaging shows that JS399-19 inhibits fungal growth but not to the extent previously reported. The fungicide causes the hypha to become entangled and unable to extend vertically. This implies that type I myosin in Fusarium is essential for hyphal and mycelia propagation. The K216Q substitution correlates with reduced susceptibility in F. avenaceum.

On simulated annealing phase transitions in phylogeny reconstruction

Phylogeny reconstruction with global criteria is NP-complete or NP-hard, hence in general requires a heuristic search. We investigate the powerful, physically inspired, general-purpose heuristic simulated annealing, applied to phylogeny reconstruction. Simulated annealing mimics the physical process of annealing, where a liquid is gently cooled to form a crystal. During the search, periods of elevated specific heat occur, analogous to physical phase transitions. These simulated annealing phase transitions play a crucial role in the outcome of the search. Nevertheless, they have received comparably little attention, for phylogeny or other optimisation problems. We analyse simulated annealing phase transitions during searches for the optimal phylogenetic tree for 34 real-world multiple alignments. In the same way in which melting temperatures differ between materials, we observe distinct specific heat profiles for each input file. We propose this reflects differences in the search landscape and can serve as a measure for problem difficulty and for suitability of the algorithm's parameters. We discuss application in algorithmic optimisation and as a diagnostic to assess parameterisation before computationally costly, large phylogeny reconstructions are launched. Whilst the focus here lies on phylogeny reconstruction under maximum parsimony, it is plausible that our results are more widely applicable to optimisation procedures in science and industry.

Plant VAP27 proteins: domain characterization, intracellular localization and role in plant development

The endoplasmic reticulum (ER) is connected to the plasma membrane (PM) through the plant-specific NETWORKED protein, NET3C, and phylogenetically conserved vesicle-associated membrane protein-associated proteins (VAPs). Ten VAP homologues (VAP27-1 to 27-10) can be identified in the Arabidopsis genome and can be divided into three clades. Representative members from each clade were tagged with fluorescent protein and expressed in Nicotiana benthamiana. Proteins from clades I and III localized to the ER as well as to ER/PM contact sites (EPCSs), whereas proteins from clade II were found only at the PM. Some of the VAP27-labelled EPCSs localized to plasmodesmata, and we show that the mobility of VAP27 at EPCSs is influenced by the cell wall. EPCSs closely associate with the cytoskeleton, but their structure is unaffected when the cytoskeleton is removed. VAP27-labelled EPCSs are found in most cell types in Arabidopsis, with the exception of cells in early trichome development. Arabidopsis plants expressing VAP27-GFP fusions exhibit pleiotropic phenotypes, including defects in root hair morphogenesis. A similar effect is also observed in plants expressing VAP27 RNAi. Taken together, these data indicate that VAP27 proteins used at EPCSs are essential for normal ER-cytoskeleton interaction and for plant development.

Terrace Aware Data Structure for Phylogenomic Inference from Supermatrices

In phylogenomics the analysis of concatenated gene alignments, the so-called supermatrix, is commonly accompanied by the assumption of partition models. Under such models each gene, or more generally partition, is allowed to evolve under its own evolutionary model. Although partition models provide a more comprehensive analysis of supermatrices, missing data may hamper the tree search algorithms due to the existence of phylogenetic (partial) terraces. Here, we introduce the phylogenetic terrace aware (PTA) data structure for the efficient analysis under partition models. In the presence of missing data PTA exploits (partial) terraces and induced partition trees to save computation time. We show that an implementation of PTA in IQ-TREE leads to a substantial speedup of up to 4.5 and 8 times compared with the standard IQ-TREE and RAxML implementations, respectively. PTA is generally applicable to all types of partition models and common topological rearrangements thus can be employed by all phylogenomic inference software.

From commensalism to parasitism in Carapidae (Ophidiiformes): heterochronic modes of development?

Phenotypic variations allow a lineage to move into new regions of the adaptive landscape. The purpose of this study is to analyse the life history of the pearlfishes (Carapinae) in a phylogenetic framework and particularly to highlight the evolution of parasite and commensal ways of life. Furthermore, we investigate the skull anatomy of parasites and commensals and discuss the developmental process that would explain the passage from one form to the other. The genus Carapus forms a paraphyletic grouping in contrast to the genus Encheliophis, which forms a monophyletic cluster. The combination of phylogenetic, morphologic and ontogenetic data clearly indicates that parasitic species derive from commensal species and do not constitute an iterative evolution from free-living forms. Although the head morphology of Carapus species differs completely from Encheliophis, C. homei is the sister group of the parasites. Interestingly, morphological characteristics allowing the establishment of the relation between Carapus homei and Encheliophis spp. concern the sound-producing mechanism, which can explain the diversification of the taxon but not the acquisition of the parasite morphotype. Carapus homei already has the sound-producing mechanism typically found in the parasite form but still has a commensal way of life and the corresponding head structure. Moreover, comparisons between the larval and adult Carapini highlight that the adult morphotype "Encheliophis" is obtained by going beyond the adult stage reached by Carapus. The entrance into the new adaptive landscape could have been realised by at least two processes: paedomorphosis and allometric repatterning.

Identification of the non-ribosomal peptide synthetase responsible for biosynthesis of the potential anti-cancer drug sansalvamide in Fusarium solani

Sansalvamide is a cyclic pentadepsipeptide produced by Fusarium solani and has shown promising results as potential anti-cancer drug. The biosynthetic pathway has until now remained unidentified, but here we used an Agrobacterium tumefaciens-mediated transformation (ATMT) approach to generate knockout mutants of two candidate non-ribosomal peptide synthetases (NRPS29 and NRPS30). Comparative studies of secondary metabolites in the two deletion mutants and wild type confirmed the absence of sansalvamide in the NRPS30 deletion mutant, implicating this synthetase in the biosynthetic pathway for sansalvamide. Sansalvamide is structurally related to the cyclic hexadepsipeptide destruxin, which both contain an α-hydroxyisocaproic acid (HICA) unit. A gene cluster responsible for destruxin production has previously been identified in Metarhizium robertsii together with a hypothetical biosynthetic pathway. Using comparative bioinformatic analyses of the catalytic domains in the destruxin and sansalvamide NRPSs, we were able to propose a model for sansalvamide biosynthesis. Orthologues of the gene clusters were also identified in species from several other genera including Acremonium chrysogenum and Trichoderma virens, which suggests that the ability to produce compounds related to destruxin and sansalvamide is widespread.

The role of defensive ecological interactions in the evolution of conotoxins

Venoms comprise of complex mixtures of peptides evolved for predation and defensive purposes. Remarkably, some carnivorous cone snails can inject two distinct venoms in response to predatory or defensive stimuli, providing a unique opportunity to study separately how different ecological pressures contribute to toxin diversification. Here, we report the extraordinary defensive strategy of the Rhizoconus subgenus of cone snails. The defensive venom from this worm-hunting subgenus is unusually simple, almost exclusively composed of αD-conotoxins instead of the ubiquitous αA-conotoxins found in the more complex defensive venom of mollusc- and fish-hunting cone snails. A similarly compartmentalized venom gland as those observed in the other dietary groups facilitates the deployment of this defensive venom. Transcriptomic analysis of a Conus vexillum venom gland revealed the αD-conotoxins as the major transcripts, with lower amounts of 15 known and four new conotoxin superfamilies also detected with likely roles in prey capture. Our phylogenetic and molecular evolution analysis of the αD-conotoxins from five subgenera of cone snails suggests they evolved episodically as part of a defensive strategy in the Rhizoconus subgenus. Thus, our results demonstrate an important role for defence in the evolution of conotoxins.

Torsins: not your typical AAA+ ATPases

Torsin ATPases (Torsins) belong to the widespread AAA+ (ATPases associated with a variety of cellular activities) family of ATPases, which share structural similarity but have diverse cellular functions. Torsins are outliers in this family because they lack many characteristics of typical AAA+ proteins, and they are the only members of the AAA+ family located in the endoplasmic reticulum and contiguous perinuclear space. While it is clear that Torsins have essential roles in many, if not all metazoans, their precise cellular functions remain elusive. Studying Torsins has significant medical relevance since mutations in Torsins or Torsin-associated proteins result in a variety of congenital human disorders, the most frequent of which is early-onset torsion (DYT1) dystonia, a severe movement disorder. A better understanding of the Torsin system is needed to define the molecular etiology of these diseases, potentially enabling corrective therapy. Here, we provide a comprehensive overview of the Torsin system in metazoans, discuss functional clues obtained from various model systems and organisms and provide a phylogenetic and structural analysis of Torsins and their regulatory cofactors in relation to disease-causative mutations. Moreover, we review recent data that have led to a dramatically improved understanding of these machines at a molecular level, providing a foundation for investigating the molecular defects underlying the associated movement disorders. Lastly, we discuss our ideas on how recent progress may be utilized to inform future studies aimed at determining the cellular role(s) of these atypical molecular machines and their implications for dystonia treatment options.

New data on two remarkable Antarctic species Amblydorylaimusisokaryon (Loof, 1975) Andrássy, 1998 and Pararhyssocolpusparadoxus (Loof, 1975), gen. n., comb. n. (Nematoda, Dorylaimida)

The taxonomic position of two antarctic dorylaimid species Amblydorylaimusisokaryon (Loof, 1975) Andrássy, 1998 and Pararhyssocolpusparadoxus (Loof, 1975), gen. n., comb. n. are discussed on the basis of morphological, including SEM study, morphometric, postembryonic and sequence data of 18S rDNA and the D2-D3 expansion fragments of large subunit rDNA. The evolutionary trees inferred from 18S sequences show insufficient resolution to determine the assignment of the two species to particular families, moreover Pararhyssocolpusparadoxus gen. n., comb. n. (=Rhyssocolpusparadoxus) previously regarded as a member of Nordiidae or Qudsianematidae, showed distant relationship both to Rhyssocolpusvinciguerrae and Eudorylaimus spp. The phylogram inferred from 28S sequences revealed that Amblydorylaimusisokaryon is a member of a well-supported group comprised of several Aporcelaimellus spp., while, no close relationships could be revealed for the Pararhyssocolpusparadoxus gen. n., comb. n. to any nematode genus. On the basis of molecular data and morphological characteristics, some taxonomic changes are proposed. Amblydorylaimusisokaryon is transferred from family Qudsianematidae to family Aporcelaimidae, and a new monotypic genus Pararhyssocolpus gen. n. is proposed, attributed to Pararhyssocolpidaefam. n. The diagnosis of the new family is provided together with emended diagnosis of the genera Amblydorylaimus and Pararhyssocolpus gen. n. Data concerning distribution of these endemic genera in the Antarctic region are also given.

Against All Odds: Trehalose-6-Phosphate Synthase and Trehalase Genes in the Bdelloid Rotifer Adineta vaga Were Acquired by Horizontal Gene Transfer and Are Upregulated during Desiccation

The disaccharide sugar trehalose is essential for desiccation resistance in most metazoans that survive dryness; however, neither trehalose nor the enzymes involved in its metabolism have ever been detected in bdelloid rotifers despite their extreme resistance to desiccation. Here we screened the genome of the bdelloid rotifer Adineta vaga for genes involved in trehalose metabolism. We discovered a total of four putative trehalose-6-phosphate synthase (TPS) and seven putative trehalase (TRE) gene copies in the genome of this ameiotic organism; however, no trehalose-6-phosphate phosphatase (TPP) gene or domain was detected. The four TPS copies of A. vaga appear more closely related to plant and fungi proteins, as well as to some protists, whereas the seven TRE copies fall in bacterial clades. Therefore, A. vaga likely acquired its trehalose biosynthesis and hydrolysis genes by horizontal gene transfers. Nearly all residues important for substrate binding in the predicted TPS domains are highly conserved, supporting the hypothesis that several copies of the genes might be functional. Besides, RNAseq library screening showed that trehalase genes were highly expressed compared to TPS genes, explaining probably why trehalose had not been detected in previous studies of bdelloids. A strong overexpression of their TPS genes was observed when bdelloids enter desiccation, suggesting a possible signaling role of trehalose-6-phosphate or trehalose in this process.

Detection and characterisation of novel bocavirus (genus Bocaparvovirus) and gastroenteritis viruses from asymptomatic pigs in Ireland

Background

Livestock animals have been the assumed source of several human epidemics in recent years, for example, influenza H1N1, rotavirus G8/G9, and MERS-CoV. Surveillance of novel viruses in animals is essential to evaluate the risk to human and animal health and to determine any economic impact, for example, failure to thrive. There is a paucity of data regarding detection and characterisation of gastroenteritis viruses, particularly novel viruses, in porcines in Ireland. Recently, a number of small novel porcine DNA viruses have emerged globally, for example, torque teno sus virus, porcine bocavirus, and parvoviruses 2 & 4, and little is known about the biology and potential pathogenicity of these viruses. Bocaparvovirus is a genetically distinct group of viruses which has been recently detected in humans and animals.

Methods

In this study, the presence of gastroenteritis viruses (rotavirus A, porcine circovirus, adenovirus, and porcine bocavirus) was investigated in a selection of archived faecal samples from asymptomatic piglets from a commercial farm in Ireland. A total of 104 specimens were pooled and screened using conventional molecular techniques (PCR and RT-PCR), a subset of specimens (n=44) were then examined individually. Viral diversity was then investigated using statistical and phylogenetic techniques.

Results

Initial screening showed a high prevalence of PBoV in this farm, with the formation of three distinct groups in phylogenetic analysis. Other viruses were also investigated in this study with the first report of PCV, PAdV and lineage I G5 RVA in Ireland. Some specimens contained >1 virus, with statistical analysis indicating a strong correlation for mixed infections of PBoV and PAdV on this farm.

Conclusion

Investigating the diversity of circulating enteric viruses on Irish porcine farms is important to improve the prophylactic tools available and to facilitate the early detection of changes in circulating viruses.

Characterization and distribution of mating-type genes of the turfgrass pathogen Sclerotinia homoeocarpa on a global scale

Sclerotinia homoeocarpa F.T. Bennett is a filamentous member of Ascomycota that causes dollar spot, the most economically important disease of turfgrass worldwide. We sequenced and characterized the mating-type (MAT) locus of four recently-collected contemporary strains causing dollar spot, four historical type strains used to describe the fungus, and three species of Rutstroemiaceae. Moreover, we developed a multiplex PCR assay to screen 1019 contemporary isolates for mating-type. The organization of the MAT loci of all strains examined could be classified into one of four categories: (1) putatively heterothallic, as exemplified by all contemporary strains and three of four historical type strains; (2) putatively heterothallic with a deleted putative gene in the MAT1-2 idiomorph, as detected in strains from two recently-collected populations in the United Kingdom that show more similarity to historical strains; (3) putatively homothallic with close physical linkage between MAT1-1-1 and MAT1-2-1, as found in one historical type strain of S. homoeocarpa and two strains of Rutstroemia cuniculi; and (4) an unresolved but apparently homothallic organization in which strains contained both MAT1-1-1 and MAT1-2-1 but linkage between these genes and between the two flanking genes could not be confirmed, as identified in R. paludosa and Poculum henningsianum. In contemporary S. homoeocarpa populations there was no significant difference in the frequency of the two mating types in clone-corrected samples when analyzed on regional and local scales, suggesting sex may be possible in this pathogen. However, two isolates from Italy and twenty from California were heterokaryotic for both complete heterothallic MAT idiomorphs. Results from this study contribute to knowledge about mating systems in filamentous fungi and enhance our understanding of the evolution and biology of an important plant pathogen.

The Complete Sequence of the Acacia ligulata Chloroplast Genome Reveals a Highly Divergent clpP1 Gene

Legumes are a highly diverse angiosperm family that include many agriculturally important species. To date, 21 complete chloroplast genomes have been sequenced from legume crops confined to the Papilionoideae subfamily. Here we report the first chloroplast genome from the Mimosoideae, Acacia ligulata, and compare it to the previously sequenced legume genomes. The A. ligulata chloroplast genome is 158,724 bp in size, comprising inverted repeats of 25,925 bp and single-copy regions of 88,576 bp and 18,298 bp. Acacia ligulata lacks the inversion present in many of the Papilionoideae, but is not otherwise significantly different in terms of gene and repeat content. The key feature is its highly divergent clpP1 gene, normally considered essential in chloroplast genomes. In A. ligulata, although transcribed and spliced, it probably encodes a catalytically inactive protein. This study provides a significant resource for further genetic research into Acacia and the Mimosoideae. The divergent clpP1 gene suggests that Acacia will provide an interesting source of information on the evolution and functional diversity of the chloroplast Clp protease complex.

The systematics of Echinorhynchus Zoega in Müller, 1776 (Acanthocephala, Echinorhynchidae) elucidated by nuclear and mitochondrial sequence data from eight European taxa

The acanthocephalan genus Echinorhynchus Zoega in Müller, 1776 (sensuYamaguti 1963) is a large and widespread group of parasites of teleost fish and malacostracan crustaceans, distributed from the Arctic to the Antarctic in habitats ranging from freshwaters to the deep-sea. A total of 52 species are currently recognised based on the conventional morphological species concept; however, the true diversity in the genus is masked by cryptic speciation. The considerable diversity within Echinorhynchus is an argument for subdividing the genus if monophyletic groups with supporting morphological characters can be identified. With this objective in mind, partial sequences of two genes with different rates of evolution and patterns of inheritance (nuclear 28S rRNA and mitochondrial cytochrome c oxidase subunit I) were used to infer the phylogenetic relationships among eight taxa of Echinorhynchus. These included representatives of each of three genus group taxa proposed in a controversial revision of the genus based on cement gland pattern, namely Echinorhynchus (sensu stricto), Metechinorhynchus Petrochenko, 1956 and Pseudoechinorhynchus Petrochenko, 1956. These groupings have previously been rejected by some authorities, because the diagnostic character is poorly defined; this study shows that Echinorhynchus (sensu stricto) and Metechinorhynchus are not natural, monophyletic groups. A revision of Echinorhynchus will require tandem molecular phylogenetic and morphological analyses of a larger sample of taxa, but this study has identified two morhological characters that might potentially be used to define new genera. The estimated phylogeny also provides insight into the zoogeographical history of Echinorhynchus spp. We postulate that the ancestral Echinorhynchus had a freshwater origin and the genus subsequently invaded the sea, probably several times. The freshwater taxa of the Echinorhynchusbothniensis Zdzitowiecki & Valtonen, 1987 clade may represent a reinvasion of freshwater by one or more ancestral marine species.