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 machine-learning-based alternative to phylogenetic bootstrap.

MOTIVATION: Currently used methods for estimating branch support in phylogenetic analyses often rely on the classic Felsenstein's bootstrap, parametric tests, or their approximations. As these branch support scores are widely used in phylogenetic analyses, having accurate, fast, and interpretable scores is of high importance. RESULTS: Here, we employed a data-driven approach to estimate branch support values with a probabilistic interpretation. To this end, we simulated thousands of realistic phylogenetic trees and the corresponding multiple sequence alignments. Each of the obtained alignments was used to infer the phylogeny using state-of-the-art phylogenetic inference software, which was then compared to the true tree. Using these extensive data, we trained machine-learning algorithms to estimate branch support values for each bipartition within the maximum-likelihood trees obtained by each software. Our results demonstrate that our model provides fast and more accurate probability-based branch support values than commonly used procedures. We demonstrate the applicability of our approach on empirical datasets. AVAILABILITY AND IMPLEMENTATION: The data supporting this work are available in the Figshare repository at https://doi.org/10.6084/m9.figshare.25050554.v1, and the underlying code is accessible via GitHub at https://github.com/noaeker/bootstrap_repo.

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/.

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.

A cnidarian parasite of salmon (Myxozoa: Henneguya) lacks a mitochondrial genome.

Although aerobic respiration is a hallmark of eukaryotes, a few unicellular lineages, growing in hypoxic environments, have secondarily lost this ability. In the absence of oxygen, the mitochondria of these organisms have lost all or parts of their genomes and evolved into mitochondria-related organelles (MROs). There has been debate regarding the presence of MROs in animals. Using deep sequencing approaches, we discovered that a member of the Cnidaria, the myxozoan Henneguya salminicola, has no mitochondrial genome, and thus has lost the ability to perform aerobic cellular respiration. This indicates that these core eukaryotic features are not ubiquitous among animals. Our analyses suggest that H. salminicola lost not only its mitochondrial genome but also nearly all nuclear genes involved in transcription and replication of the mitochondrial genome. In contrast, we identified many genes that encode proteins involved in other mitochondrial pathways and determined that genes involved in aerobic respiration or mitochondrial DNA replication were either absent or present only as pseudogenes. As a control, we used the same sequencing and annotation methods to show that a closely related myxozoan, Myxobolus squamalis, has a mitochondrial genome. The molecular results are supported by fluorescence micrographs, which show the presence of mitochondrial DNA in M. squamalis, but not in H. salminicola. Our discovery confirms that adaptation to an anaerobic environment is not unique to single-celled eukaryotes, but has also evolved in a multicellular, parasitic animal. Hence, H. salminicola provides an opportunity for understanding the evolutionary transition from an aerobic to an exclusive anaerobic metabolism.

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.

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.

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 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.

Genomic insights into the evolutionary origin of Myxozoa within Cnidaria.

The Myxozoa comprise over 2,000 species of microscopic obligate parasites that use both invertebrate and vertebrate hosts as part of their life cycle. Although the evolutionary origin of myxozoans has been elusive, a close relationship with cnidarians, a group that includes corals, sea anemones, jellyfish, and hydroids, is supported by some phylogenetic studies and the observation that the distinctive myxozoan structure, the polar capsule, is remarkably similar to the stinging structures (nematocysts) in cnidarians. To gain insight into the extreme evolutionary transition from a free-living cnidarian to a microscopic endoparasite, we analyzed genomic and transcriptomic assemblies from two distantly related myxozoan species, Kudoa iwatai and Myxobolus cerebralis, and compared these to the transcriptome and genome of the less reduced cnidarian parasite, Polypodium hydriforme. A phylogenomic analysis, using for the first time to our knowledge, a taxonomic sampling that represents the breadth of myxozoan diversity, including four newly generated myxozoan assemblies, confirms that myxozoans are cnidarians and are a sister taxon to P. hydriforme. Estimations of genome size reveal that myxozoans have one of the smallest reported animal genomes. Gene enrichment analyses show depletion of expressed genes in categories related to development, cell differentiation, and cell-cell communication. In addition, a search for candidate genes indicates that myxozoans lack key elements of signaling pathways and transcriptional factors important for multicellular development. Our results suggest that the degeneration of the myxozoan body plan from a free-living cnidarian to a microscopic parasitic cnidarian was accompanied by extreme reduction in genome size and gene content.

Back to solitude: Solving the phylogenetic position of the Diazonidae using molecular and developmental characters.

The order Aplousobranchia (Chordata, Ascidiacea) contains approximately 1500 species distributed worldwide. Their phylogeny, however, remains unclear, with unresolved family relationships. While most Aplousobranchia are colonial, debates exist concerning the phylogenetic position of families such as the Diazonidae and Cionidae, which exhibit a solitary lifestyle and share morphological characteristics with both Aplousobranchia and Phlebobranchia orders. To clarify the phylogenetic position of the Diazonidae and Cionidae, we determined the complete mitochondrial sequence of the solitary diazonid Rhopalaea idoneta. The phylogenetic reconstruction based on the 13 mitochondrial protein coding genes strongly supports a positioning of Diazonidae well-nested within the Aplousobranchia rather than a positioning as a sister clade of the Aplousobranchia. In addition, we examined the regenerative ability of R. idoneta. Similar to colonial Aplousobranchia, R. idoneta was found to be able to completely regenerate its thorax. Ciona, also known to possess high regenerative abilities, is the Aplousobranchia sister clade rather than a member of the Phlebobranchia. Our results thus indicate that the colonial lifestyle was acquired in the Aplousobranchia, starting from a Ciona-like solitary ancestor and secondarily lost in Diazonidae representatives such as Rhopalaea. The solitary lifestyle of Rhopalaea is thus a derived characteristic rather than an ancestral trait.

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 (p) verrucosa (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.

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.

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.

Evolutionary Insights from the Mitochondrial Genome of Oikopleura dioica: Sequencing Challenges, RNA Editing, Gene Transfers to the Nucleus, and tRNA Loss.

Sequencing the mitochondrial genome of the tunicate Oikopleura dioica is a challenging task due to the presence of long poly-A/T homopolymer stretches, which impair sequencing and assembly. Here, we report on the sequencing and annotation of the majority of the mitochondrial genome of O. dioica by means of combining several DNA and amplicon reads obtained by Illumina and MinIon Oxford Nanopore Technologies (ONT) with public RNA sequences. We document extensive RNA editing, since all homopolymer stretches present in the mitochondrial DNA correspond to 6U-regions in the mitochondrial RNA. Out of the 13 canonical protein-coding genes, we were able to detect eight, plus an unassigned ORF that lacked sequence similarity to canonical mitochondrial protein-coding genes. We show that the nad3 gene has been transferred to the nucleus and acquired a mitochondria-targeting signal. In addition to two very short rRNAs, we could only identify a single tRNA (tRNA-Met), suggesting multiple losses of tRNA genes, supported by a corresponding loss of mitochondrial aminoacyl-tRNA synthetases in the nuclear genome. Based on the eight canonical protein-coding genes identified, we reconstructed maximum likelihood and Bayesian phylogenetic trees and inferred an extreme evolutionary rate of this mitochondrial genome. The phylogenetic position of appendicularians among tunicates, however, could not be accurately determined.

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.

ALG11--a new variable DNA marker for sponge phylogeny: comparison of phylogenetic performances with the 18S rDNA and the COI gene.

Phylogenetic relationships within sponge classes are highly debated. The low phylogenetic signal observed with some current molecular data can be attributed to the use of few markers, usually slowly-evolving, such as the nuclear rDNA genes and the mitochondrial COI gene. In this study, we conducted a bioinformatics search for a new molecular marker. We sought a marker that (1) is likely to have no paralogs; (2) evolves under a fast evolutionary rate; (3) is part of a continuous exonic region; and (4) is flanked by conserved regions. Our search suggested the nuclear ALG11 as a potential suitable marker. We next demonstrated that this marker can indeed be used for solving phylogenetic relationships within sponges. Specifically, we successfully amplified the ALG11 gene from DNA samples of representatives from all four sponge classes as well as from several cnidarian classes. We also amplified the 18S rDNA and the COI gene for these species. Finally, we analyzed the phylogenetic performance of ALG11 to solve sponge relationships compared to and in combination with the nuclear 18S rDNA and the COI mtDNA genes. Interestingly, the ALG11 marker seems to be superior to the widely-used COI marker. Our work thus indicates that the ALG11 marker is a relevant marker which can complement and corroborate the phylogenetic inferences observed with nuclear ribosomal genes. This marker is also expected to contribute to resolving evolutionary relationships of other apparently slow-evolving animal phyla, such as cnidarians.

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.

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.

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.