Evolution of myxozoan mitochondrial genomes: insights from myxobolids

BACKGROUND

Myxozoa is a class of cnidarian parasites that encompasses over 2,400 species. Phylogenetic relationships among myxozoans remain highly debated, owing to both a lack of informative morphological characters and a shortage of molecular markers. Mitochondrial (mt) genomes are a common marker in phylogeny and biogeography. However, only five complete myxozoan mt genomes have been sequenced: four belonging to two closely related genera, Enteromyxum and Kudoa, and one from the genus Myxobolus. Interestingly, while cytochrome oxidase genes could be identified in Enteromyxum and Kudoa, no such genes were found in Myxobolus squamalis, and another member of the Myxobolidae (Henneguya salminicola) was found to have lost its entire mt genome. To evaluate the utility of mt genomes to reconstruct myxozoan relationships and to understand if the loss of cytochrome oxidase genes is a characteristic of myxobolids, we sequenced the mt genome of five myxozoans (Myxobolus wulii, M. honghuensis, M. shantungensis, Thelohanellus kitauei and, Sphaeromyxa zaharoni) using Illumina and Oxford Nanopore platforms.

RESULTS

Unlike Enteromyxum, which possesses a partitioned mt genome, the five mt genomes were encoded on single circular chromosomes. An mt plasmid was found in M. wulii, as described previously in Kudoa iwatai. In all new myxozoan genomes, five protein-coding genes (cob, cox1, cox2, nad1, and nad5) and two rRNAs (rnl and rns) were recognized, but no tRNA. We found that Myxobolus and Thelohanellus species shared unidentified reading frames, supporting the view that these mt open reading frames are functional. Our phylogenetic reconstructions based on the five conserved mt genes agree with previously published trees based on the 18S rRNA gene.

CONCLUSIONS

Our results suggest that the loss of cytochrome oxidase genes is not a characteristic of all myxobolids, the ancestral myxozoan mt genome was likely encoded on a single circular chromosome, and mt plasmids exist in a few lineages. Our findings indicate that myxozoan mt sequences are poor markers for reconstructing myxozoan phylogenetic relationships because of their fast-evolutionary rates and the abundance of repeated elements, which complicates assembly.

GenomeFLTR: filtering reads made easy

In the last decade, advances in sequencing technology have led to an exponential increase in genomic data. These new data have dramatically changed our understanding of the evolution and function of genes and genomes. Despite improvements in sequencing technologies, identifying contaminated reads remains a complex task for many research groups. Here, we introduce GenomeFLTR, a new web server to filter contaminated reads. Reads are compared against existing sequence databases from various representative organisms to detect potential contaminants. The main features implemented in GenomeFLTR are: (i) automated updating of the relevant databases; (ii) fast comparison of each read against the database; (iii) the ability to create user-specified databases; (iv) a user-friendly interactive dashboard to investigate the origin and frequency of the contaminations; (v) the generation of a contamination-free file. Availability: https://genomefltr.tau.ac.il/.

Octocorals in the Gulf of Aqaba exhibit high photosymbiont fidelity

Symbiotic associations, widespread in terrestrial and marine ecosystems, are of considerable ecological importance. Many tropical coral species are holobionts, formed by the obligate association between a cnidarian host and endosymbiotic dinoflagellates of the family Symbiodiniaceae. The latter are abundant on coral reefs from very shallow water down to the upper mesophotic zone (30-70 m). The research on scleractinians has revealed that the photosymbiont lineages present in the cnidarian host play an important role in the coral's ability to thrive under different environmental conditions, such as light regime and temperature. However, little is known regarding octocoral photosymbionts, and in particular regarding those found deeper than 30 m. Here, we used ribosomal (ITS2) and chloroplast (23S) markers to uncover, for the first time, the dominant Symbiodiniaceae taxa present in 19 mesophotic octocoral species (30-70 m depth) from the Gulf of Aqaba/Eilat (northern Red Sea). In addition, using high-throughput sequencing of the ITS2 region we characterized both the dominant and the rare Symbiodiniaceae lineages found in several species across depth. The phylogenetic analyses of both markers were in agreement and revealed that most of the studied mesophotic octocorals host the genus Cladocopium. Litophyton spp. and Klyxum utinomii were exceptions, as they harbored Symbiodinium and Durusdinium photosymbionts, respectively. While the dominant algal lineage of each coral species did not vary across depth, the endosymbiont community structure significantly differed between host species, as well as between different depths for some host species. The findings from this study contribute to the growing global-catalogue of Cnidaria-Symbiodiniaceae associations. Unravelling the Symbiodiniaceae composition in octocoral holobionts across environmental gradients, depth in particular, may enable a better understanding of how specialized those associations are, and to what extent coral holobionts are able to modify their photosymbionts.

The Phylogenetic Position of the Enigmatic, Polypodium hydriforme (Cnidaria, Polypodiozoa): Insights from Mitochondrial Genomes

Polypodium hydriforme is an enigmatic parasite that belongs to the phylum Cnidaria. Its taxonomic position has been debated: whereas it was previously suggested to be part of Medusozoa, recent phylogenomic analyses based on nuclear genes support the view that P. hydriforme and Myxozoa form a clade called Endocnidozoa. Medusozoans have linear mitochondrial (mt) chromosomes, whereas myxozoans, as most metazoan species, have circular chromosomes. In this work, we determined the structure of the mt genome of P. hydriforme, using Illumina and Oxford Nanopore Technologies reads, and showed that it is circular. This suggests that P. hydriforme is not nested within Medusozoa, as this would entail linearization followed by recirculation. Instead, our results support the view that P. hydriforme is a sister clade to Myxozoa, and mt linearization in the lineage leading to medusozoans occurred after the divergence of Myxozoa + P. hydriforme. Detailed analyses of the assembled P. hydriforme mt genome show that: (1) it is encoded on a single circular chromosome with an estimated size of ∼93,000 base pairs, making it one of the largest metazoan mt genomes; (2) around 78% of the genome encompasses a noncoding region composed of several repeat types; (3) similar to Myxozoa, no mt tRNAs were identified; (4) the codon TGA is a stop codon and does not encode for tryptophan as in other cnidarians; (5) similar to myxozoan mt genomes, it is extremely fast evolving.

Myxozoan infection in thinlip mullet Chelon ramada (Mugiliformes: Mugilidae) in the Sea of Galilee

Mullets (Mugilidae) are economically important fish in Israel. Two species of mugilids (i.e., the thinlip mullet Chelon ramada and the flathead grey mullet Mugil cephalus) have been stocked in the Sea of Galilee (Lake Kinneret) in order to increase fishermen's income and lake water quality. These catadromous species do not reproduce in the lake, consequently, fingerlings have been introduced every year since 1958. Following a survey of myxozoan infections in the Sea of Galilee, we described Myxobolus pupkoi n. sp. infecting the gill arches, and reported Myxobolus exiguus from visceral peritoneum and gall bladder of C. ramada. The prevalence of infection of both Myxobolus pupkoi n. sp. and M. exiguus were 11.5% (2/23). Our study indicates that the parasites infecting C. ramada belong to a lineage of myxozoans infecting mugilids. This result suggests that the infection took place in the Mediterranean Sea, where the fingerlings were caught, before their introduction into the Sea of Galilee. Since 2018 only farm-raised fingerlings have been introduced. We thus recommend to closely monitor the presence of these parasites in the future to determine if the presence of parasites disappear with the introduction of farm-raised fingerlings.

Two novel myxosporean parasite species of Ceratomyxa Thélohan, 1892 from the banded cusk-eel Raneya brasiliensis (Kaup) (Ophidiiformes: Ophidiidae) off Patagonia, Argentina

We described two novel myxozoan parasite species Ceratomyxa argentina n. sp. and Ceratomyxa raneyae n. sp. from the gall bladder of Raneya brasiliensis (Kaup) from the Patagonian coast of Argentina. Both species can be distinguished from other ceratomyxids by myxospore and polar capsule (nematocyst) morphology and morphometry, fish host and geographic locality. Phylogenetic reconstruction using ssrDNA gene sequences showed that the two new species are placed in a long-branching ceratomyxid clade which also include Ceratomyxa appendiculata Thélohan, 1892, Ceratomyxa anko Freeman, Yokoyama and Ogawa, 2008, Ceratomyxa pantherini Gunter, Burger and Adlard, 2010 and Pseudoalataspora kovalevae Kalavati, MacKenzie, Collins, Hemmingsen and Brickle, 2013. This study documents additional biodiversity of marine myxozoans in the South Atlantic, a region still largely unexplored for this group of parasitic cnidarians.

Myxosporea (Myxozoa, Cnidaria) Lack DNA Cytosine Methylation

DNA cytosine methylation is central to many biological processes, including regulation of gene expression, cellular differentiation, and development. This DNA modification is conserved across animals, having been found in representatives of sponges, ctenophores, cnidarians, and bilaterians, and with very few known instances of secondary loss in animals. Myxozoans are a group of microscopic, obligate endoparasitic cnidarians that have lost many genes over the course of their evolution from free-living ancestors. Here, we investigated the evolution of the key enzymes involved in DNA cytosine methylation in 29 cnidarians and found that these enzymes were lost in an ancestor of Myxosporea (the most speciose class of Myxozoa). Additionally, using whole-genome bisulfite sequencing, we confirmed that the genomes of two distant species of myxosporeans, Ceratonova shasta and Henneguya salminicola, completely lack DNA cytosine methylation. Our results add a notable and novel taxonomic group, the Myxosporea, to the very short list of animal taxa lacking DNA cytosine methylation, further illuminating the complex evolutionary history of this epigenetic regulatory mechanism.

And Then There Were Three…: Extreme Regeneration Ability of the Solitary Chordate Polycarpa mytiligera

Extensive regenerative ability is a common trait of animals capable of asexual development. The current study reveals the extraordinary regeneration abilities of the solitary ascidian Polycarpa mytiligera. Dissection of a single individual into separate fragments along two body axes resulted in the complete regeneration of each fragment into an independent, functional individual. The ability of a solitary ascidian, incapable of asexual development, to achieve bidirectional regeneration and fully regenerate all body structures and organs is described here for the first time. Amputation initiated cell proliferation in proximity to the amputation line. Phylogenetic analysis demonstrated the close affinity of P. mytiligera to colonial species. This evolutionary proximity suggests the ability for regeneration as an exaptation feature for colonial lifestyle. P. mytiligera’s exceptional regenerative abilities and phylogenetic position highlight its potential to serve as a new comparative system for studies seeking to uncover the evolution of regeneration and coloniality among the chordates.

Molecular relationships of the Israeli shrews (Eulipotyphla: Soricidae) based on cytochrome b sequences

The number of shrew species in Israel has been and still is the subject of debate. In this work we used for the first time a molecular marker, the cytochrome b gene, to investigate the number and identity of shrew species in Israel. Our molecular results confirmed the presence of four species: Crocidura leucodon, Crocidura suaveolens gueldenstaedtii, Crocidura ramona, and Suncus etruscus. The C. ramona sequences were found to differ from all other Crocidura species sequenced to date, supporting its status as a distinct species. Whether it is conspecific with Crocidura portali (described in 1920 from Israel and usually synonymized with C. suaveolens), will require additional study. The sequences of Israeli C. suaveolens were found to be very similar to those of Iran, Turkey, and Georgia (i.e., C. suaveolens gueldenstaedtii), in agreement with previous studies. The Israeli C. leucodon sequences, however, formed a distinct clade among C. leucodon. Finally, the S. etruscus sequences clustered with sequences from France, Italy, and Iran.

Illuminating Genetic Mysteries of the Dead Sea Scrolls

The discovery of the 2,000-year-old Dead Sea Scrolls had an incomparable impact on the historical understanding of Judaism and Christianity. "Piecing together" scroll fragments is like solving jigsaw puzzles with an unknown number of missing parts. We used the fact that most scrolls are made from animal skins to "fingerprint" pieces based on DNA sequences. Genetic sorting of the scrolls illuminates their textual relationship and historical significance. Disambiguating the contested relationship between Jeremiah fragments supplies evidence that some scrolls were brought to the Qumran caves from elsewhere; significantly, they demonstrate that divergent versions of Jeremiah circulated in parallel throughout Israel (ancient Judea). Similarly, patterns discovered in non-biblical scrolls, particularly the Songs of the Sabbath Sacrifice, suggest that the Qumran scrolls represent the broader cultural milieu of the period. Finally, genetic analysis divorces debated fragments from the Qumran scrolls. Our study demonstrates that interdisciplinary approaches enrich the scholar's toolkit.

A new mitochondrial gene order in the banded cusk-eel Raneya brasiliensis (Actinopterygii, Ophidiiformes)

The complete mitochondrial genome of the banded cusk-eel, Raneya brasilensis (Kaup, 1856), was obtained using next-generation sequencing approaches. The genome sequence was 16,881 bp and exhibited a novel gene order for a vertebrate. Specifically, the WANCY and the nd6 – D-loop regions were re-ordered, supporting the hypothesis that these two regions are hotspots for gene rearrangements in Actinopterygii. Phylogenetic reconstructions confirmed that R. brasiliensis is nested within Ophidiiformes. Mitochondrial genomes are required from additional ophidiins to determine whether the gene rearrangements that we observed are specific to the genus Raneya or to the subfamily Ophidiinae.

A genome wide survey reveals multiple nematocyst-specific genes in Myxozoa

Background

Myxozoa represents a diverse group of microscopic endoparasites whose life cycle involves two hosts: a vertebrate (usually a fish) and an invertebrate (usually an annelid worm). Despite lacking nearly all distinguishing animal characteristics, given that each life cycle stage consists of no more than a few cells, molecular phylogenetic studies have revealed that myxozoans belong to the phylum Cnidaria, which includes corals, sea anemones, and jellyfish. Myxozoa, however, do possess a polar capsule; an organelle that is homologous to the stinging structure unique to Cnidaria: the nematocyst. Previous studies have identified in Myxozoa a number of protein-coding genes that are specific to nematocytes (the cells producing nematocysts) and thus restricted to Cnidaria. Determining which other genes are also homologous with the myxozoan polar capsule genes could provide insight into both the conservation and changes that occurred during nematocyst evolution in the transition to endoparasitism.

Results

Previous studies have examined the phylogeny of two cnidarian-restricted gene families: minicollagens and nematogalectins. Here we identify and characterize seven additional cnidarian-restricted genes in myxozoan genomes using a phylogenetic approach. Four of the seven had never previously been identified as cnidarian-specific and none have been studied in a phylogenetic context. A majority of the proteins appear to be involved in the structure of the nematocyst capsule and tubule. No venom proteins were identified among the cnidarian-restricted genes shared by myxozoans.

Conclusions

Given the highly divergent forms that comprise Cnidaria, obtaining insight into the processes underlying their ancient diversification remains challenging. In their evolutionary transition to microscopic endoparasites, myxozoans lost nearly all traces of their cnidarian ancestry, with the one prominent exception being their nematocysts (or polar capsules). Thus nematocysts, and the genes that code for their structure, serve as rich sources of information to support the cnidarian origin of Myxozoa.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1253-7) contains supplementary material, which is available to authorized users.

Extensive mitochondrial gene rearrangements in Ctenophora: insights from benthic Platyctenida

BACKGROUND

Complete mitochondrial (mt) genomes have been sequenced for thousands of animals and represent a molecule of choice for many evolutionary studies. Nevertheless, some animal groups have remained under-sampled. Ctenophora (comb jellies) is one such example, with only two complete mt sequences determined hitherto for this phylum, which encompasses ca. 150-200 described species. This lack of data derives from the extremely fast mt evolutionary rate in this lineage, complicating primer design and DNA amplification. Indeed, in the two ctenophore mt genomes sequenced to date, i.e. those of Mnemiopsis leidyi (order Lobata) and Pleurobrachia bachei (order Cydippida), both rRNA and protein coding genes exhibit an extraordinary size reduction and have highly derived sequences. Additionally, all tRNAs, and the atp6 and atp8 genes are absent. In order to determine whether these characteristics are shared by other ctenophores, we obtained the complete mt genomes of three benthic ctenophores belonging to the so far unsampled order of Platyctenida: Coeloplana loyai, Coeloplana yulianicorum and Vallicula multiformis.

RESULTS

The mt genomes of benthic ctenophores reveal the same peculiarities found in Mnemiopsis and Pleurobrachia, demonstrating that the fast evolutionary rate is a general trait of the ctenophore mt genomes. Our results also indicate that this high evolutionary rate not only affects the nucleotide substitution but also gene rearrangements. Indeed, gene order was highly rearranged among representatives of the different taxonomic orders in which it was close to random, but also quite variable within Platyctenida, in which the genera Coeloplana and Vallicula share only four conserved synteny blocks. However, the two congeneric Coeloplana species display exactly the same gene order. Because of the extreme evolutionary rate, our phylogenetic analyses were unable to resolve the phylogenetic position of ctenophores within metazoans or the relationships among the different Ctenophora orders. Comparative sequence-analyses allowed us to correct the annotation of the Pleurobrachia mt genome, confirming the absence of tRNAs, the presence of both rRNA genes, and the existence of a reassignment of codon TGA from tryptophan to serine for this species.

CONCLUSIONS

Since Platyctenida is an early diverging lineage among Ctenophora, our findings suggest that the mt traits described above are ancestral characteristics of this phylum.

The First Mitochondrial Genomics and Evolution SMBE-Satellite Meeting: A New Scientific Symbiosis

The central role of the mitochondrion for cellular and organismal metabolism is well known, yet its functional role in evolution has rarely been featured in leading international conferences. Moreover, the contribution of mitochondrial genetics to complex disease phenotypes is particularly important, and although major advances have been made in the field of genomics, mitochondrial genomic data have in many cases been overlooked. Accumulating data and new knowledge support a major contribution of this maternally inherited genome, and its interactions with the nucleus, to both major evolutionary processes and diverse disease phenotypes. These advances encouraged us to assemble the first Mitochondrial Genomics and Evolution (MGE) meeting-an SMBE satellite and Israeli Science foundation international conference (Israel, September 2017). Here, we report the content and outcome of the MGE meeting (https://www.mge2017.com/; last accessed November 5, 2017).

The Multipartite Mitochondrial Genome of Enteromyxum leei (Myxozoa): Eight Fast-Evolving Megacircles

Myxozoans are a large group of poorly characterized cnidarian parasites. To gain further insight into their evolution, we sequenced the mitochondrial (mt) genome of Enteromyxum leei and reevaluate the mt genome structure of Kudoa iwatai. Although the typical animal mt genome is a compact, 13-25 kb, circular chromosome, the mt genome of E. leei was found to be fragmented into eight circular chromosomes of ∼23 kb, making it the largest described animal mt genome. Each chromosome was found to harbor a large noncoding region (∼15 kb), nearly identical between chromosomes. The protein coding genes show an unusually high rate of sequence evolution and possess little similarity to their cnidarian homologs. Only five protein coding genes could be identified and no tRNA genes. Surprisingly, the mt genome of K. iwatai was also found to be composed of two chromosomes. These observations confirm the remarkable plasticity of myxozoan mt genomes.

Eastern Mediterranean Mobility in the Bronze and Early Iron Ages: Inferences from Ancient DNA of Pigs and Cattle

The Late Bronze of the Eastern Mediterranean (1550–1150 BCE) was a period of strong commercial relations and great prosperity, which ended in collapse and migration of groups to the Levant. Here we aim at studying the translocation of cattle and pigs during this period. We sequenced the first ancient mitochondrial and Y chromosome DNA of cattle from Greece and Israel and compared the results with morphometric analysis of the metacarpal in cattle. We also increased previous ancient pig DNA datasets from Israel and extracted the first mitochondrial DNA for samples from Greece. We found that pigs underwent a complex translocation history, with links between Anatolia with southeastern Europe in the Bronze Age, and movement from southeastern Europe to the Levant in the Iron I (ca. 1150–950 BCE). Our genetic data did not indicate movement of cattle between the Aegean region and the southern Levant. We detected the earliest evidence for crossbreeding between taurine and zebu cattle in the Iron IIA (ca. 900 BCE). In light of archaeological and historical evidence on Egyptian imperial domination in the region in the Late Bronze Age, we suggest that Egypt attempted to expand dry farming in the region in a period of severe droughts.

Mitochondrial group I and group II introns in the sponge orders Agelasida and Axinellida

Background

Self-splicing introns are present in the mitochondria of members of most eukaryotic lineages. They are divided into Group I and Group II introns, according to their secondary structure and splicing mechanism. Being rare in animals, self-splicing introns were only described in a few sponges, cnidarians, placozoans and one annelid species. In sponges, three types of mitochondrial Group I introns were previously described in two demosponge families (Tetillidae, and Aplysinellidae) and in the homoscleromorph family Plakinidae. These three introns differ in their insertion site, secondary structure and in the sequence of the LAGLIDADG gene they encode. Notably, no group II introns have been previously described in sponges.

Results

We report here the presence of mitochondrial introns in the cytochrome oxidase subunit 1 (COI) gene of three additional sponge species from three different families: Agelas oroides (Agelasidae, Agelasida), Cymbaxinella^pverrucosa (Hymerhabdiidae, Agelasida) and Axinella polypoides (Axinellidae, Axinellida). We show, for the first time, that sponges can also harbour Group II introns in their COI gene, whose presence in animals’ mitochondria has so far been described in only two phyla, Placozoa and Annelida. Surprisingly, two different Group II introns were discovered in the COI gene of C. verrucosa. Phylogenetic analysis indicates that the Group II introns present in C. verrucosa are related to red algae (Rhodophyta) introns.

Conclusions

The differences found among intron secondary structures and the phylogenetic inferences support the hypothesis that the introns originated from independent horizontal gene transfer events. Our results thus suggest that self-splicing introns are more diverse in the mitochondrial genome of sponges than previously anticipated.

Electronic supplementary material

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

Indel reliability in indel-based phylogenetic inference

It is often assumed that it is unlikely that the same insertion or deletion (indel) event occurred at the same position in two independent evolutionary lineages, and thus, indel-based inference of phylogeny should be less subject to homoplasy compared with standard inference which is based on substitution events. Indeed, indels were successfully used to solve debated evolutionary relationships among various taxonomical groups. However, indels are never directly observed but rather inferred from the alignment and thus indel-based inference may be sensitive to alignment errors. It is hypothesized that phylogenetic reconstruction would be more accurate if it relied only on a subset of reliable indels instead of the entire indel data. Here, we developed a method to quantify the reliability of indel characters by measuring how often they appear in a set of alternative multiple sequence alignments. Our approach is based on the assumption that indels that are consistently present in most alternative alignments are more reliable compared with indels that appear only in a small subset of these alignments. Using simulated and empirical data, we studied the impact of filtering and weighting indels by their reliability scores on the accuracy of indel-based phylogenetic reconstruction. The new method is available as a web-server at http://guidance.tau.ac.il/RELINDEL/.

Diversity and evolution of myxozoan minicollagens and nematogalectins

Background

Myxozoa are a diverse group of metazoan parasites with a very simple organization, which has for decades eluded their evolutionary origin. Their most prominent and characteristic feature is the polar capsule: a complex intracellular structure of the myxozoan spore, which plays a role in host infection. Striking morphological similarities have been found between myxozoan polar capsules and nematocysts, the stinging structures of cnidarians (corals, sea anemones and jellyfish) leading to the suggestion that Myxozoa and Cnidaria share a more recent common ancestry. This hypothesis has recently been supported by phylogenomic evidence and by the identification of a nematocyst specific minicollagen gene in the myxozoan Tetracapsuloides bryosalmonae. Here we searched genomes and transcriptomes of several myxozoan taxa for the presence of additional cnidarian specific genes and characterized these genes within a phylogenetic context.

Results

Illumina assemblies of transcriptome or genome data of three myxozoan species (Enteromyxum leei, Kudoa iwatai, and Sphaeromyxa zaharoni) and of the enigmatic cnidarian parasite Polypodium hydriforme (Polypodiozoa) were mined using tBlastn searches with nematocyst-specific proteins as queries. Several orthologs of nematogalectins and minicollagens were identified. Our phylogenetic analyses indicate that myxozoans possess three distinct minicollagens. We found that the cnidarian repertoire of nematogalectins is more complex than previously thought and we identified additional members of the nematogalectin family. Cnidarians were found to possess four nematogalectin/ nematogalectin-related genes, while in myxozoans only three genes could be identified.

Conclusions

Our results demonstrate that myxozoans possess a diverse array of genes that are taxonomically restricted to Cnidaria. Characterization of these genes provide compelling evidence that polar capsules and nematocysts are homologous structures and that myxozoans are highly degenerate cnidarians. The diversity of minicollagens was higher than previously thought, with the presence of three minicollagen genes in myxozoans. Our phylogenetic results suggest that the different myxozoan sequences are the results of ancient divergences within Cnidaria and not of recent specializations of the polar capsule. For both minicollagen and nematogalectin, our results show that myxozoans possess less gene copies than their cnidarian counter parts, suggesting that the polar capsule gene repertoire was simplified with their reduced body plan.

Electronic supplementary material

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

The complete mitochondrial genome of the devil firefish Pterois miles (Bennett, 1828) (Scorpaenidae)

The complete mitochondrial genome of the devil firefish Pterois miles (Bennett, 1828) was obtained using next generation sequencing approaches. The genome sequence was comprised of 16,497 bp exhibiting the standard vertebrate mitochondrial gene arrangement. Regions of gene overlap, tRNA lengths, as well as start and stop codons were similar to those observed in closely related families (i.e. Sebastidae, Peristediidae). Phylogenetic reconstructions support the polyphyly of Scorpaeniformes, and confirm the close relationship of Scorpaenidae and Sebastidae.

The complete mitochondrial genome of the gilthead seabream Sparus aurata L. (Sparidae)

The complete mitochondrial genome of the gilthead seabream Sparus aurata Linnaeus 1758, one of the world's most important mariculture species, was sequenced using next generation sequencing technology. The genome sequence is comprised of 16,652 bp exhibiting the canonical vertebrate mitochondria gene order. Regions of gene overlap, tRNA length, as well as start and stop codon were similar to those observed in other Sparidae. Phylogenetic reconstructions based on mitochondrial protein coding genes corroborate the view that Sparidae is paraphyletic and includes Centracanthidae.

Biology of a new xenoma-forming gonadotropic microsporidium in the invasive blotchfin dragonet Callionymus filamentosus

A gonadotropic microsporidian parasite, Obruspora papernae gen. et sp. nov. (Microsporidia: Enterocytozoonidae), is described from Callionymus filamentosus (Teleostei: Callionymidae) in the Mediterranean Sea. The host, a Red Sea invasive species which entered the Mediterranean through the Suez Canal, was first collected in the Levant Basin in 1953, whereas its parasite went unobserved until 2008. Analysis of partial small subunit ribosomal gene sequences (SSU rDNA) placed the new species within the Nucleospora, Desmozoon, and Paranucleospora clade, and as it differs from each of them, it is assigned to a new genus. The development of the parasite is described, and the biological mechanisms underlying this parasite-host system are analyzed. Prevalence of infection approached 80% in female samples throughout most of the year. Males showed no signs of infection, but parasite rDNA was detected in male internal organs. The parasite-induced xenomas progressively occupied and eventually replaced much of the ovary, in some cases producing effective castration. Despite high levels of parasite infection, current trawl fishery statistics indicate that the abundance of Mediterranean populations of the host remains high. The parasite impact on the host population dynamics is unclear. Possible effects of the new microsporidian parasite on the reproductive effort of C. filamentosus and the potential role of another parasite, the ectoparasitic copepod Lernanthropus callionymicola, as an additional host in the life cycle of O. papernae, require further investigation.

Deep sequencing of mixed total DNA without barcodes allows efficient assembly of highly plastic ascidian mitochondrial genomes

Ascidians or sea squirts form a diverse group within chordates, which includes a few thousand members of marine sessile filter-feeding animals. Their mitochondrial genomes are characterized by particularly high evolutionary rates and rampant gene rearrangements. This extreme variability complicates standard polymerase chain reaction (PCR) based techniques for molecular characterization studies, and consequently only a few complete Ascidian mitochondrial genome sequences are available. Using the standard PCR and Sanger sequencing approach, we produced the mitochondrial genome of Ascidiella aspersa only after a great effort. In contrast, we produced five additional mitogenomes (Botrylloides aff. leachii, Halocynthia spinosa, Polycarpa mytiligera, Pyura gangelion, and Rhodosoma turcicum) with a novel strategy, consisting in sequencing the pooled total DNA samples of these five species using one Illumina HiSeq 2000 flow cell lane. Each mitogenome was efficiently assembled in a single contig using de novo transcriptome assembly, as de novo genome assembly generally performed poorly for this task. Each of the new six mitogenomes presents a different and novel gene order, showing that no syntenic block has been conserved at the ordinal level (in Stolidobranchia and in Phlebobranchia). Phylogenetic analyses support the paraphyly of both Ascidiacea and Phlebobranchia, with Thaliacea nested inside Phlebobranchia, although the deepest nodes of the Phlebobranchia-Thaliacea clade are not well resolved. The strategy described here thus provides a cost-effective approach to obtain complete mitogenomes characterized by a highly plastic gene order and a fast nucleotide/amino acid substitution rate.

Ancient DNA and population turnover in southern levantine pigs--signature of the sea peoples migration?

Near Eastern wild boars possess a characteristic DNA signature. Unexpectedly, wild boars from Israel have the DNA sequences of European wild boars and domestic pigs. To understand how this anomaly evolved, we sequenced DNA from ancient and modern pigs from Israel. Pigs from Late Bronze Age (until ca. 1150 BCE) in Israel shared haplotypes of modern and ancient Near Eastern pigs. European haplotypes became dominant only during the Iron Age (ca. 900 BCE). This raises the possibility that European pigs were brought to the region by the Sea Peoples who migrated to the Levant at that time. Then, a complete genetic turnover took place, most likely because of repeated admixture between local and introduced European domestic pigs that went feral. Severe population bottlenecks likely accelerated this process. Introductions by humans have strongly affected the phylogeography of wild animals, and interpretations of phylogeography based on modern DNA alone should be taken with caution.

Phylogeny of Tetillidae (Porifera, Demospongiae, Spirophorida) based on three molecular markers

Tetillidae are spherical to elliptical cosmopolitan demosponges. The family comprises eight genera: namely, Acanthotetilla Burton, 1959, Amphitethya Lendenfeld, 1907, CinachyraSollas, 1886, CinachyrellaWilson, 1925, Craniella Schmidt, 1870, Fangophilina Schmidt, 1880, Paratetilla Dendy, 1905, and Tetilla Schmidt, 1868. These genera are characterized by few conflicting morphological characters, resulting in an ambiguity of phylogenetic relationships. The phylogeny of tetillid genera was investigated using the cox1, 18S rRNA and 28S rRNA (C1-D2 domains) genes in 88 specimens (8 genera, 28 species). Five clades were identified: (i) Cinachyrella, Paratetilla and Amphitethya species, (ii) Cinachyrella levantinensis, (iii) Tetilla, (iv) Craniella, Cinachyra and Fangophilina and (v) Acanthotetilla. Consequently, the phylogenetic analysis supports the monophyly of Tetilla, a genus lacking any known morphological synapomorphy. Acanthotetilla is also recovered. In contrast, within the first clade, species of the genera Paratetilla and Amphitethya were nested within Cinachyrella. Similarly, within the fourth clade, species of the genera Cinachyra and Fangophilina were nested within Craniella. As previously postulated by taxonomists, the loss of ectodermal specialization (i.e., a cortex) has occurred several times independently. Nevertheless, the presence or absence of a cortex and its features carry a phylogenetic signal. Surprisingly, the common view that assumes close relationships among sponges with porocalices (i.e., surface depressions) is refuted.

Mitochondrial and morphological variation of Tilapia zillii in Israel

BACKGROUND

Tilapia zillii is widespread in the East Levant inland aquatic systems as well as in artificial water reservoirs. In this study we explore the genetic and morphological variation of this widespread species, using mitochondrial control region sequences and meristic characters. We examine the hypothesis that T. zillii's population structure corresponds to the four Israeli aquatic systems.

RESULTS

Out of seven natural water bodies, only two were found to possess genetically divergent populations of T. zillii. In addition to its presence in fish farms, the species was found in two artificial recreational ponds which were supposed to have been stocked only with other fish species. In these two artificial habitats, the haplotype frequencies diverged significantly from those of natural populations. Finally, fish from the Dead Sea springs of Ne'ot HaKikar appear to differ both genetically and morphologically from fish of the same aquatic system but not from fish of other water systems.

CONCLUSIONS

Our results show that the population structure of T. zillii does not match the geography of the Israeli water-basins, with the exception of the Dead Sea and Kishon River, when considering natural populations only. The absence of a significant divergence between basins is discussed. Our results and observations suggest that the Ne'ot HaKikar Dead Sea population and those of artificial ponds could have originated from the "hitchhiking" of T. zillii, at the expense of some other cultivated tilapiine species.

GLOOME: gain loss mapping engine

The evolutionary analysis of presence and absence profiles (phyletic patterns) is widely used in biology. It is assumed that the observed phyletic pattern is the result of gain and loss dynamics along a phylogenetic tree. Examples of characters that are represented by phyletic patterns include restriction sites, gene families, introns and indels, to name a few. Here, we present a user-friendly web server that accurately infers branch-specific and site-specific gain and loss events. The novel inference methodology is based on a stochastic mapping approach utilizing models that reliably capture the underlying evolutionary processes. A variety of features are available including the ability to analyze the data with various evolutionary models, to infer gain and loss events using either stochastic mapping or maximum parsimony, and to estimate gain and loss rates for each character analyzed.

AVAILABILITY

Freely available for use at http://gloome.tau.ac.il/.

Diversity of sponge mitochondrial introns revealed by cox 1 sequences of Tetillidae

Background

Animal mitochondrial introns are rare. In sponges and cnidarians they have been found in the cox 1 gene of some spirophorid and homosclerophorid sponges, as well as in the cox 1 and nad 5 genes of some Hexacorallia. Their sporadic distribution has raised a debate as to whether these mobile elements have been vertically or horizontally transmitted among their hosts. The first sponge found to possess a mitochondrial intron was a spirophorid sponge from the Tetillidae family. To better understand the mode of transmission of mitochondrial introns in sponges, we studied cox 1 intron distribution among representatives of this family.

Results

Seventeen tetillid cox 1 sequences were examined. Among these sequences only six were found to possess group I introns. Remarkably, three different forms of introns were found, named introns 714, 723 and 870 based on their different positions in the cox 1 alignment. These introns had distinct secondary structures and encoded LAGLIDADG ORFs belonging to three different lineages. Interestingly, sponges harboring the same intron form did not always form monophyletic groups, suggesting that their introns might have been transferred horizontally. To evaluate whether the introns were vertically or horizontally transmitted in sponges and cnidarians we used a host parasite approach. We tested for co-speciation between introns 723 (the introns with the highest number of sponge representatives) and their nesting cox 1 sequences. Reciprocal AU tests indicated that the intron and cox 1 tree are significantly different, while a likelihood ratio test was not significant. A global test of co-phylogeny had significant results; however, when cnidarian sequences were analyzed separately the results were not significant.

Conclusions

The co-speciation analyses thus suggest that a vertical transmission of introns in the ancestor of sponges and cnidarians, followed by numerous independent losses, cannot solely explain the current distribution of metazoan group I introns. An alternative scenario that includes horizontal gene transfer events appears to be more suitable to explain the incongruence between the intron 723 and the cox 1 topologies. In addition, our results suggest that three different intron forms independently colonized the cox 1 gene of tetillids. Among sponges, the Tetillidae family seems to be experiencing an unusual number of intron insertions.

Dasyatispora levantinae gen. et sp. nov., a new microsporidian parasite from the common stingray Dasyatis pastinaca in the eastern Mediterranean

A new microsporidian infecting the Mediterranean common stingray Dasyatis pastinaca (Linnaeus, 1758) is described from Iskenderun Bay, Turkey. The parasite invades the disc muscles, producing slender, spindle-shaped subcutaneous swellings that develop into massive, elongated, tumor-like protuberances measuring up to 11 x 4 cm. Severity of the infection may vary from light (1 or 2 small lesions) to intense, with large parts of the dorsal surface covered with lumps and protrusions. These masses contained a yellowish-white caseous substance consisting of degraded host tissue and microsporidian sporophorous vesicles, which in turn contained developing sporonts, sporoblasts and spores. The ripe spore contained a uni-nucleate sporoplasm and large posterior vacuole, and measured 3.8-4.3 x 2.6-2.8 microm. Infection prevalence was 21% in a sample of 143 host individuals examined. All the infected stingray individuals were within the weight class of 300 to 800 g (200 to 305 mm disc width). Phylogenetic analyses of rDNA sequences indicate that this microsporidian belongs to the Pleistophoridae and clusters with species of the genera Ovipleistophora Pekkarinen, Lom & Nilsen, 2002 and Heterosporis Schubert, 1969. However, the morphology, development and host differ distinctly from all reported species, including those belonging to these 2 genera, and it is thus assigned to a newly erected genus and named Dasyatispora levantinae gen. et sp. nov. This is the first record of a microsporidian infection in a batoid. It is also the first microsporidian species to be formally described from an elasmobranch.

Tunicate mitogenomics and phylogenetics: peculiarities of the Herdmania momus mitochondrial genome and support for the new chordate phylogeny

Background

Tunicates represent a key metazoan group as the sister-group of vertebrates within chordates. The six complete mitochondrial genomes available so far for tunicates have revealed distinctive features. Extensive gene rearrangements and particularly high evolutionary rates have been evidenced with regard to other chordates. This peculiar evolutionary dynamics has hampered the reconstruction of tunicate phylogenetic relationships within chordates based on mitogenomic data.

Results

In order to further understand the atypical evolutionary dynamics of the mitochondrial genome of tunicates, we determined the complete sequence of the solitary ascidian Herdmania momus. This genome from a stolidobranch ascidian presents the typical tunicate gene content with 13 protein-coding genes, 2 rRNAs and 24 tRNAs which are all encoded on the same strand. However, it also presents a novel gene arrangement, highlighting the extreme plasticity of gene order observed in tunicate mitochondrial genomes. Probabilistic phylogenetic inferences were conducted on the concatenation of the 13 mitochondrial protein-coding genes from representatives of major metazoan phyla. We show that whereas standard homogeneous amino acid models support an artefactual sister position of tunicates relative to all other bilaterians, the CAT and CAT+BP site- and time-heterogeneous mixture models place tunicates as the sister-group of vertebrates within monophyletic chordates. Moreover, the reference phylogeny indicates that tunicate mitochondrial genomes have experienced a drastic acceleration in their evolutionary rate that equally affects protein-coding and ribosomal-RNA genes.

Conclusion

This is the first mitogenomic study supporting the new chordate phylogeny revealed by recent phylogenomic analyses. It illustrates the beneficial effects of an increased taxon sampling coupled with the use of more realistic amino acid substitution models for the reconstruction of animal phylogeny.

Large-scale parsimony analysis of metazoan indels in protein-coding genes

Insertions and deletions (indels) are considered to be rare evolutionary events, the analysis of which may resolve controversial phylogenetic relationships. Indeed, indel characters are often assumed to be less homoplastic than amino acid and nucleotide substitutions and, consequently, more reliable markers for phylogenetic reconstruction. In this study, we analyzed indels from over 1,000 metazoan orthologous genes. We studied the impact of different species sampling, ortholog data sets, lengths of included indels, and indel-coding methods on the resulting metazoan tree. Our results show that, similar to sequence substitutions, indels are homoplastic characters, and their analysis is sensitive to the long-branch attraction artifact. Furthermore, improving the taxon sampling and choosing a closely related outgroup greatly impact the phylogenetic inference. Our indel-based inferences support the Ecdysozoa hypothesis over the Coelomata hypothesis and suggest that sponges are a sister clade to other animals.

An updated 18S rRNA phylogeny of tunicates based on mixture and secondary structure models

BACKGROUND

Tunicates have been recently revealed to be the closest living relatives of vertebrates. Yet, with more than 2500 described species, details of their evolutionary history are still obscure. From a molecular point of view, tunicate phylogenetic relationships have been mostly studied based on analyses of 18S rRNA sequences, which indicate several major clades at odds with the traditional class-level arrangements. Nonetheless, substantial uncertainty remains about the phylogenetic relationships and taxonomic status of key groups such as the Aplousobranchia, Appendicularia, and Thaliacea.

RESULTS

Thirty new complete 18S rRNA sequences were acquired from previously unsampled tunicate species, with special focus on groups presenting high evolutionary rate. The updated 18S rRNA dataset has been aligned with respect to the constraint on homology imposed by the rRNA secondary structure. A probabilistic framework of phylogenetic reconstruction was adopted to accommodate the particular evolutionary dynamics of this ribosomal marker. Detailed Bayesian analyses were conducted under the non-parametric CAT mixture model accounting for site-specific heterogeneity of the evolutionary process, and under RNA-specific doublet models accommodating the occurrence of compensatory substitutions in stem regions. Our results support the division of tunicates into three major clades: 1) Phlebobranchia + Thaliacea + Aplousobranchia, 2) Appendicularia, and 3) Stolidobranchia, but the position of Appendicularia could not be firmly resolved. Our study additionally reveals that most Aplousobranchia evolve at extremely high rates involving changes in secondary structure of their 18S rRNA, with the exception of the family Clavelinidae, which appears to be slowly evolving. This extreme rate heterogeneity precluded resolving with certainty the exact phylogenetic placement of Aplousobranchia. Finally, the best fitting secondary-structure and CAT-mixture models suggest a sister-group relationship between Salpida and Pyrosomatida within Thaliacea.

CONCLUSION

An updated phylogenetic framework for tunicates is provided based on phylogenetic analyses using the most realistic evolutionary models currently available for ribosomal molecules and an unprecedented taxonomic sampling. Detailed analyses of the 18S rRNA gene allowed a clear definition of the major tunicate groups and revealed contrasting evolutionary dynamics among major lineages. The resolving power of this gene nevertheless appears limited within the clades composed of Phlebobranchia + Thaliacea + Aplousobranchia and Pyuridae + Styelidae, which were delineated as spots of low resolution. These limitations underline the need to develop new nuclear markers in order to further resolve the phylogeny of this keystone group in chordate evolution.

Rodent phylogeny revised: analysis of six nuclear genes from all major rodent clades

BACKGROUND

Rodentia is the most diverse order of placental mammals, with extant rodent species representing about half of all placental diversity. In spite of many morphological and molecular studies, the family-level relationships among rodents and the location of the rodent root are still debated. Although various datasets have already been analyzed to solve rodent phylogeny at the family level, these are difficult to combine because they involve different taxa and genes.

RESULTS

We present here the largest protein-coding dataset used to study rodent relationships. It comprises six nuclear genes, 41 rodent species, and eight outgroups. Our phylogenetic reconstructions strongly support the division of Rodentia into three clades: (1) a "squirrel-related clade", (2) a "mouse-related clade", and (3) Ctenohystrica. Almost all evolutionary relationships within these clades are also highly supported. The primary remaining uncertainty is the position of the root. The application of various models and techniques aimed to remove non-phylogenetic signal was unable to solve the basal rodent trifurcation.

CONCLUSION

Sequencing and analyzing a large sequence dataset enabled us to resolve most of the evolutionary relationships among Rodentia. Our findings suggest that the uncertainty regarding the position of the rodent root reflects the rapid rodent radiation that occurred in the Paleocene rather than the presence of conflicting phylogenetic and non-phylogenetic signals in the dataset.

Bird mitochondrial gene order: insight from 3 warbler mitochondrial genomes

Two main gene orders exist in birds: the ancestral gene order and the remnant control region (CR) 2 gene order. These gene orders differ by the presence of 1 or 2 copies of the CR, respectively. Among songbirds, Oscines were thought to follow the ancestral gene order, with the exception of the lyrebird and Phylloscopus warblers. Here, we determined the complete mitochondrial genome sequence of 3 non-Phylloscopus warblers species and found that the blackcap (Sylvia atricapilla) and the reed warbler (Acrocephalus scirpaceus) have 2 almost identical copies of the CR, whereas the eastern orphean warbler (Sylvia crassirostris) follows the remnant CR 2 gene order. Our results contradict previous studies suggesting that Acrocephalus and most sylvioid warblers exhibit the ancestral gene order. We were able to trace this contradiction to a misidentification of gene order from polymerase chain reaction length determination. We thus suggest that passerine gene order evolution needs to be revised.

Phylogeny of coral-inhabiting barnacles (Cirripedia; Thoracica; Pyrgomatidae) based on 12S, 16S and 18S rDNA analysis

The traditional phylogeny of the coral-inhabiting barnacles, the Pyrgomatidae, is based on morphological characteristics, mainly of the hard parts. It has been difficult to establish the phylogenetic relationships among Pyrgomatidae because of the apparent convergence of morphological characteristics, and due to the use of non-cladistic systematics, which emphasize ancestor-descendant relationships rather than sister-clade relationships. We used partial sequences of two mithochondrial genes, 12S rDNA and 16S rDNA, and a nuclear gene, 18S rDNA, to infer the molecular phylogeny of the pyrgomatids. Our phylogenetic results allowed us to reject previous classifications of Pyrgomatidae based on morphological characteristics. Our results also suggested the possibility of paraphyly of the Pyrgomatidae. The hydrocoral barnacle Wanella is not found on the same clade as the other pyrgomatids, but rather, with the free-living balanids. The basal position of Megatrema and Ceratoconcha is supported. The archeaobalanid Armatobalanus is grouped with Cantellius at the base of the Indo-Pacific pyrgomatines. Fusion of the shell plate and modification of the opercular valves are homoplasious features that occurred more than three times on different clades. The monophyly of the "Savignium" group, comprising four nominal genera, is also not supported, and the different taxa are placed on different clades.

Multiple molecular evidences for a living mammalian fossil

Laonastes aenigmamus is an enigmatic rodent first described in 2005. Molecular and morphological data suggested that it is the sole representative of a new mammalian family, the Laonastidae, and a member of the Hystricognathi. However, the validity of this family is controversial because fossil-based phylogenetic analyses suggest that Laonastes is a surviving member of the Diatomyidae, a family considered to have been extinct for 11 million years. According to these data, Laonastes and Diatomyidae are the sister clade of extant Ctenodactylidae (i.e., gundies) and do not belong to the Hystricognathi. To solve the phylogenetic position of Laonastes, we conducted a large-scale molecular phylogeny of rodents. The analysis includes representatives of all major rodent taxonomic groups and was based on 5.5 kb of sequence data from four nuclear and two mitochondrial genes. To further validate the obtained results, a short interspersed element insertion analysis including 11 informative loci was also performed. Our molecular data based on sequence and short interspersed element analyses unambiguously placed Laonastes as a sister clade of gundies. All alternative hypotheses were significantly rejected based on Shimodaira-Hasegawa tests, supporting the idea that Laonastes does not belong to the Hystricognathi. Molecular dating analysis also supports an ancient divergence, approximately 44 Mya ago, between Ctenodactylidae and Laonastes. These combined analyses support the hypothesis that Laonastes is indeed a living fossil. Protection of this surviving species would conserve an ancient mammalian family.

Automated scanning for phylogenetically informative transposed elements in rodents

Transposed elements constitute an attractive, useful source of phylogenetic markers to elucidate the evolutionary history of their hosts. Frequent and successive amplifications over evolutionary time are important requirements for utilizing their presence or absence as landmarks of evolution. Although transposed elements are well distributed in rodent taxa, the generally high degree of genomic sequence divergence among species complicates our access to presence/absence data. With this in mind we developed a novel, high-throughput computational strategy, called CPAL (Conserved Presence/Absence Locus-finder), to identify genome-wide distributed, phylogenetically informative transposed elements flanked by highly conserved regions. From a total of 232 extracted chromosomal mouse loci we randomly selected 14 of these plus 2 others from previous test screens and attempted to amplify them via PCR in representative rodent species. All loci were amplifiable and ultimately contributed 31 phylogenetically informative markers distributed throughout the major groups of Rodentia.

Putative cross-kingdom horizontal gene transfer in sponge (Porifera) mitochondria

BACKGROUND

The mitochondrial genome of Metazoa is usually a compact molecule without introns. Exceptions to this rule have been reported only in corals and sea anemones (Cnidaria), in which group I introns have been discovered in the cox1 and nad5 genes. Here we show several lines of evidence demonstrating that introns can also be found in the mitochondria of sponges (Porifera).

RESULTS

A 2,349 bp fragment of the mitochondrial cox1 gene was sequenced from the sponge Tetilla sp. (Spirophorida). This fragment suggests the presence of a 1143 bp intron. Similar to all the cnidarian mitochondrial introns, the putative intron has group I intron characteristics. The intron is present in the cox1 gene and encodes a putative homing endonuclease. In order to establish the distribution of this intron in sponges, the cox1 gene was sequenced from several representatives of the demosponge diversity. The intron was found only in the sponge order Spirophorida. A phylogenetic analysis of the COI protein sequence and of the intron open reading frame suggests that the intron may have been transmitted horizontally from a fungus donor.

CONCLUSION

Little is known about sponge-associated fungi, although in the last few years the latter have been frequently isolated from sponges. We suggest that the horizontal gene transfer of a mitochondrial intron was facilitated by a symbiotic relationship between fungus and sponge. Ecological relationships are known to have implications at the genomic level. Here, an ecological relationship between sponge and fungus is suggested based on the genomic analysis.

Arrival and diversification of caviomorph rodents and platyrrhine primates in South America

Platyrrhine primates and caviomorph rodents are clades of mammals that colonized South America during its period of isolation from the other continents, between 100 and 3 million years ago (Mya). Until now, no molecular study investigated the timing of the South American colonization by these two lineages with the same molecular data set. Using sequences from three nuclear genes (ADRA2B, vWF, and IRBP, both separate and combined) from 60 species, and eight fossil calibration constraints, we estimated the times of origin and diversification of platyrrhines and caviomorphs via a Bayesian relaxed molecular clock approach. To account for the possible effect of an accelerated rate of evolution of the IRBP gene along the branch leading to the anthropoids, we performed the datings with and without IRBP (3768 sites and 2469 sites, respectively). The time window for the colonization of South America by primates and by rodents is demarcated by the dates of origin (upper bound) and radiation (lower bound) of platyrrhines and caviomorphs. According to this approach, platyrrhine primates colonized South America between 37.0 +/- 3.0 Mya (or 38.9 +/- 4.0 Mya without IRBP) and 16.8 +/- 2.3 (or 20.1 +/- 3.3) Mya, and caviomorph rodents between 45.4 +/- 4.1 (or 43.7 +/- 4.8) Mya and 36.7 +/- 3.7 (or 35.8 +/- 4.3) Mya. Considering both the fossil record and these molecular datings, the favored scenarios are a trans-Atlantic migration of primates from Africa at the end of the Eocene or beginning of the Oligocene, and a colonization of South America by rodents during the Middle or Late Eocene. Based on our nuclear DNA data, we cannot rule out the possibility of a concomitant arrival of primates and rodents in South America. The caviomorphs radiated soon after their arrival, before the Oligocene glaciations, and these early caviomorph lineages persisted until the present. By contrast, few platyrrhine fossils are known in the Oligocene, and the present-day taxa are the result of a quite recent, Early Miocene diversification.

16S rRNA phylogeny of sponge-associated cyanobacteria

Phylogenetic analyses of 16S rRNA sequences of sponge-associated cyanobacteria showed them to be polyphyletic, implying that they derived from multiple independent symbiotic events. Most of the symbiont sequences were affiliated to a group of Synechococcus and Prochlorococcus species. However, other symbionts were related to different groups, such as the Oscillatoriales.

Rabbits, if anything, are likely Glires

Rodentia (e.g., mice, rats, dormice, squirrels, and guinea pigs) and Lagomorpha (e.g., rabbits, hares, and pikas) are usually grouped into the Glires. Status of this controversial superorder has been evaluated using morphology, paleontology, and mitochondrial plus nuclear DNA sequences. This growing corpus of data has been favoring the monophyly of Glires. Recently, Misawa and Janke [Mol. Phylogenet. Evol. 28 (2003) 320] analyzed the 6441 amino acids of 20 nuclear proteins for six placental mammals (rat, mouse, rabbit, human, cattle, and dog) and two outgroups (chicken and xenopus), and observed a basal position of the two murine rodents among the former. They concluded that "the Glires hypothesis was rejected." We here reanalyzed [loc. cit.] data set under maximum likelihood and Bayesian tree-building approaches, using phylogenetic models that take into account among-site variation in evolutionary rates and branch-length variation among proteins. Our observations support both the association of rodents and lagomorphs and the monophyly of Euarchontoglires (=Supraprimates) as the most likely explanation of the protein alignments. We conducted simulation studies to evaluate the appropriateness of lissamphibian and avian outgroups to root the placental tree. When the outgroup-to-ingroup evolutionary distance increases, maximum parsimony roots the topology along the long Mus-Rattus branch. Maximum likelihood, in contrast, roots the topology along different branches as a function of their length. Maximum likelihood appears less sensitive to the "long-branch attraction artifact" than is parsimony. Our phylogenetic conclusions were confirmed by the analysis of a different protein data set using a similar sample of species but different outgroups. We also tested the effect of the addition of afrotherian and xenarthran taxa. Using the linearized tree method, [loc. cit.] estimated that mice and rats diverged about 35 million years ago. Molecular dating based on the Bayesian relaxed molecular clock method suggests that the 95% credibility interval for the split between mice and rats is 7-17 Mya. We here emphasize the need for appropriate models of sequence evolution (matrices of amino acid replacement, taking into account among-site rate variation, and independent parameters across independent protein partitions) and for a taxonomically broad sample, and conclude on the likelihood that rodents and lagomorphs together constitute a monophyletic group (Glires).

Local molecular clocks in three nuclear genes: divergence times for rodents and other mammals and incompatibility among fossil calibrations

Reconstructing the chronology of mammalian evolution is a debated issue between molecule- and fossil-based inferences. A methodological limitation of molecules is the evolutionary rate variation among lineages, precluding the application of the global molecular clock. We considered 2422 first and second codon positions of the combined ADRA2B, IRBP, and vWF nuclear genes for a well-documented set of placentals including an extensive sampling of rodents. Using seven independent calibration points and a maximum-likelihood framework, we evaluated whether molecular and paleontological estimates of mammalian divergence dates may be reconciled by the local molecular clocks approach, allowing local constancy of substitution rates with variations at larger phylogenetic scales. To handle the difficulty of choosing among all possible rate assignments for various lineages, local molecular clocks were based on the results of branch-length and two-cluster tests. Extensive lineage-specific variation of evolutionary rates was detected, even among rodents. Cross-calibrations indicated some incompatibilities between divergence dates based on different paleontological references. To decrease the impact of a single calibration point, estimates derived from independent calibrations displaying only slight reciprocal incompatibility were averaged. The divergence dates inferred for the split between mice and rats (approximately 13-19 Myr) was younger than previously published molecular estimates. The most recent common ancestors of rodents, primates and rodents, boreoeutherians, and placentals were estimated to be, respectively, approximately 60, 70, 75, and 78 Myr old. Global clocks, local clocks, and quartet dating analyses suggested a Late Cretaceous origin of the crown placental clades followed by a Tertiary radiation of some placental orders like rodents.

Combining multiple data sets in a likelihood analysis: which models are the best?

Until recently, phylogenetic analyses have been routinely based on homologous sequences of a single gene. Given the vast number of gene sequences now available, phylogenetic studies are now based on the analysis of multiple genes. Thus, it has become necessary to devise statistical methods to combine multiple molecular data sets. Here, we compare several models for combining different genes for the purpose of evaluating the likelihood of tree topologies. Three methods of branch length estimation were studied: assuming all genes have the same branch lengths (concatenate model), assuming that branch lengths are proportional among genes (proportional model), or assuming that each gene has a separate set of branch lengths (separate model). We also compared three models of among-site rate variation: the homogenous model, a model that assumes one gamma parameter for all genes, and a model that assumes one gamma parameter for each gene. On the basis of two nuclear and one mitochondrial amino acid data sets, our results suggest that, depending on the data set chosen, either the separate model or the proportional model represents the most appropriate method for branch length analysis. For all the data sets examined, one gamma parameter for each gene represents the best model for among-site rate variation. Using these models we analyzed alternative mammalian tree topologies, and we describe the effect of the assumed model on the maximum likelihood tree. We show that the choice of the model has an impact on the best phylogeny obtained.