Divergence time estimates for the hypoxia-inducible factor-1 alpha (HIF1α) reveal an ancient emergence of animals in low-oxygen environments

Unveiling the tempo and mode of animal evolution is necessary to understand the links between environmental changes and biological innovation. Although the earliest unambiguous metazoan fossils date to the late Ediacaran period, molecular clock estimates agree that the last common ancestor (LCA) of all extant animals emerged ~850 Ma, in the Tonian period, before the oldest evidence for widespread ocean oxygenation at ~635-560 Ma in the Ediacaran period. Metazoans are aerobic organisms, that is, they are dependent on oxygen to survive. In low-oxygen conditions, most animals have an evolutionarily conserved pathway for maintaining oxygen homeostasis that triggers physiological changes in gene expression via the hypoxia-inducible factor (HIFa). However, here we confirm the absence of the characteristic HIFa protein domain responsible for the oxygen sensing of HIFa in sponges and ctenophores, indicating the LCA of metazoans lacked the functional protein domain as well, and so could have maintained their transcription levels unaltered under the very low-oxygen concentrations of their environments. Using Bayesian relaxed molecular clock dating, we inferred that the ancestral gene lineage responsible for HIFa arose in the Mesoproterozoic Era, ~1273 Ma (Credibility Interval 957-1621 Ma), consistent with the idea that important genetic machinery associated with animals evolved much earlier than the LCA of animals. Our data suggest at least two duplication events in the evolutionary history of HIFa, which generated three vertebrate paralogs, products of the two successive whole-genome duplications that occurred in the vertebrate LCA. Overall, our results support the hypothesis of a pre-Tonian emergence of metazoans under low-oxygen conditions, and an increase in oxygen response elements during animal evolution.

Multi-omics for studying and understanding polar life

Polar ecosystems are experiencing amongst the most rapid rates of regional warming on Earth. Here, we discuss 'omics' approaches to investigate polar biodiversity, including the current state of the art, future perspectives and recommendations. We propose a community road map to generate and more fully exploit multi-omics data from polar organisms. These data are needed for the comprehensive evaluation of polar biodiversity and to reveal how life evolved and adapted to permanently cold environments with extreme seasonality. We argue that concerted action is required to mitigate the impact of warming on polar ecosystems via conservation efforts, to sustainably manage these unique habitats and their ecosystem services, and for the sustainable bioprospecting of novel genes and compounds for societal gain.

Genomic characterization of a novel, widely distributed Mycoplasma species "Candidatus Mycoplasma mahonii" associated with the brittlestar Gorgonocephalus chilensis

Symbiotic relationships are ubiquitous throughout the world's oceans, yet for many marine organisms, including those in the high latitudes, little is understood about symbiotic associations and functional relationships. From a recently determined genome sequence of a filter-feeding basket star from Argentina, Gorgonocephalus chilensis, we discovered a novel Mycoplasma species with a 796Kb genome (CheckM completeness of 97.9%, G+C content = 30.1%). Similar to other Mycoplasma spp. within Mycoplasmatota, genomic analysis of the novel organism revealed reduced metabolic pathways including incomplete biosynthetic pathways, suggesting an obligate association with their basket star host. Results of 16S rRNA and multi-locus phylogenetic analyses revealed that this organism belonged to a recently characterized non-free-living lineage of Mycoplasma spp. specifically associated with marine invertebrate animals. Thus, the name "Candidatus Mycoplasma mahonii" is proposed for this novel species. Based on 16S rRNA PCR-screening, we found that Ca. M. mahonii also occurs in Gorgonocephalus eucnemis from the Northwest Pacific and other Gorgonocephalus chilensis from Argentinian waters. The level of sequence conservation within Ca. M. mahonii is considerable between widely disparate high-latitude Gorgonocephalus species, suggesting that oceanic dispersal of this microbe may be greater than excepted.

Unsupervised Deep Learning Can Identify Protein Functional Groups from Unaligned Sequences

Interpreting protein function from sequence data is a fundamental goal of bioinformatics. However, our current understanding of protein diversity is bottlenecked by the fact that most proteins have only been functionally validated in model organisms, limiting our understanding of how function varies with gene sequence diversity. Thus, accuracy of inferences in clades without model representatives is questionable. Unsupervised learning may help to ameliorate this bias by identifying highly complex patterns and structure from large datasets without external labels. Here we present DeepSeqProt, an unsupervised deep learning program for exploring large protein sequence datasets. DeepSeqProt is a clustering tool capable of distinguishing between broad classes of proteins while learning local and global structure of functional space. DeepSeqProt is capable of learning salient biological features from unaligned, unannotated sequences. DeepSeqProt is more likely to capture complete protein families and statistically significant shared ontologies within proteomes than other clustering methods. We hope this framework will prove of use to researchers and provide a preliminary step in further developing unsupervised deep learning in molecular biology.

Unexpected Distribution of Chitin and Chitin Synthase across Soft-Bodied Cnidarians

Cnidarians are commonly recognized as sea jellies, corals, or complex colonies such as the Portuguese man-of-war. While some cnidarians possess rigid internal calcareous skeletons (e.g., corals), many are soft-bodied. Intriguingly, genes coding for the chitin-biosynthetic enzyme, chitin synthase (CHS), were recently identified in the model anemone Nematostella vectensis, a species lacking hard structures. Here we report the prevalence and diversity of CHS across Cnidaria and show that cnidarian chitin synthase genes display diverse protein domain organizations. We found that CHS is expressed in cnidarian species and/or developmental stages with no reported chitinous or rigid morphological structures. Chitin affinity histochemistry indicates that chitin is present in soft tissues of some scyphozoan and hydrozoan medusae. To further elucidate the biology of chitin in cnidarian soft tissues, we focused on CHS expression in N. vectensis. Spatial expression data show that three CHS orthologs are differentially expressed in Nematostella embryos and larvae during development, suggesting that chitin has an integral role in the biology of this species. Understanding how a non-bilaterian lineage such as Cnidaria employs chitin may provide new insight into hitherto unknown functions of polysaccharides in animals, as well as their role in the evolution of biological novelty.

Induced immune reaction in the acorn worm, Saccoglossus kowalevskii, informs the evolution of antiviral immunity

Evolutionary perspectives on the deployment of immune factors following infection have been shaped by studies on a limited number of biomedical model systems with a heavy emphasis on vertebrate species. Though their contributions to contemporary immunology cannot be understated, a broader phylogenetic perspective is needed to understand the evolution of immune systems across Metazoa. In our study, we leverage differential gene expression analyses to identify genes implicated in the antiviral immune response of the acorn worm hemichordate, Saccoglossus kowalevskii, and place them in the context of immunity evolution within deuterostomes - the animal clade composed of chordates, hemichordates, and echinoderms. Following acute exposure to the synthetic viral dsRNA analog, poly(I:C), we show that S. kowalevskii responds by regulating the transcription of genes associated with canonical innate immunity signaling pathways (e.g., NF-κB and IRF signaling) and metabolic processes (e.g., lipid metabolism), as well as many genes without clear evidence of orthology with those of model species. Aggregated across all experimental time point contrasts, we identify 423 genes that are differentially expressed in response to poly(I:C). We also identify 147 genes with altered temporal patterns of expression in response to immune challenge. By characterizing the molecular toolkit involved in hemichordate antiviral immunity, our findings provide vital evolutionary context for understanding the origins of immune systems within Deuterostomia.

A comparative analysis of mitochondrial ORFs provides new insights on expansion of mitochondrial genome size in Arcidae

BACKGROUND

Arcidae, comprising about 260 species of ark shells, is an ecologically and economically important lineage of bivalve mollusks. Interestingly, mitochondrial genomes of several Arcidae species are 2-3 times larger than those of most bilaterians, and are among the largest bilaterian mitochondrial genomes reported to date. The large mitochondrial genome size is mainly due to expansion of unassigned regions (regions that are functionally unassigned). Previous work on unassigned regions of Arcidae mtDNA genomes has focused on nucleotide-level analyses to observe sequence characteristics, however the origin of expansion remains unclear.

RESULTS

We assembled six new mitogenomes and sequenced six transcriptomes of Scapharca broughtonii to identify conserved functional ORFs that are transcribed in unassigned regions. Sixteen lineage-specific ORFs with different copy numbers were identified from seven Arcidae species, and 11 of 16 ORFs were expressed and likely biologically active. Unassigned regions of 32 Arcidae mitogenomes were compared to verify the presence of these novel mitochondrial ORFs and their distribution. Strikingly, multiple structural analyses and functional prediction suggested that these additional mtDNA-encoded proteins have potential functional significance. In addition, our results also revealed that the ORFs have a strong connection to the expansion of Arcidae mitochondrial genomes and their large-scale duplication play an important role in multiple expansion events. We discussed the possible origin of ORFs and hypothesized that these ORFs may originate from duplication of mitochondrial genes.

CONCLUSIONS

The presence of lineage-specific mitochondrial ORFs with transcriptional activity and potential functional significance supports novel features for Arcidae mitochondrial genomes. Given our observation and analyses, these ORFs may be products of mitochondrial gene duplication. These findings shed light on the origin and function of novel mitochondrial genes in bivalves and provide new insights into evolution of mitochondrial genome size in metazoans.

Structure-Function Relationships of Oxygen Transport Proteins in Marine Invertebrates Enduring Higher Temperatures and Deoxygenation

AbstractPredictions for climate change-to lesser and greater extents-reveal a common scenario in which marine waters are characterized by a deadly trio of stressors: higher temperatures, lower oxygen levels, and acidification. Ectothermic taxa that inhabit coastal waters, such as shellfish, are vulnerable to rapid and prolonged environmental disturbances, such as heatwaves, pollution-induced eutrophication, and dysoxia. Oxygen transport capacity of the hemolymph (blood equivalent) is considered the proximal driver of thermotolerance and respiration in many invertebrates. Moreover, maintaining homeostasis under environmental duress is inextricably linked to the activities of the hemolymph-based oxygen transport or binding proteins. Several protein groups fulfill this role in marine invertebrates: copper-based extracellular hemocyanins, iron-based intracellular hemoglobins and hemerythrins, and giant extracellular hemoglobins. In this brief text, we revisit the distribution and multifunctional properties of oxygen transport proteins, notably hemocyanins, in the context of climate change, and the consequent physiological reprogramming of marine invertebrates.

Different phylogenomic methods support monophyly of enigmatic 'Mesozoa' (Dicyemida + Orthonectida, Lophotrochozoa)

Dicyemids and orthonectids were traditionally classified in a group called Mesozoa, but their placement in a single clade has been contested and their position(s) within Metazoa is uncertain. Here, we assembled a comprehensive matrix of Lophotrochozoa (Metazoa) and investigated the position of Dicyemida (= Rhombozoa) and Orthonectida, employing multiple phylogenomic approaches. We sequenced seven new transcriptomes and one draft genome from dicyemids (Dicyema, Dicyemennea) and two transcriptomes from orthonectids (Rhopalura). Using these and published data, we assembled and analysed contamination-filtered datasets with up to 987 genes. Our results recover Mesozoa monophyletic and as a close relative of Platyhelminthes or Gnathifera. Because of the tendency of the long-branch mesozoans to group with other long-branch taxa in our analyses, we explored the impact of approaches purported to help alleviate long-branch attraction (e.g. taxon removal, coalescent inference, gene targeting). None of these were able to break the association of Orthonectida with Dicyemida in the maximum-likelihood trees. Contrastingly, the Bayesian analysis and site-specific frequency model in maximum-likelihood did not recover a monophyletic Mesozoa (but only when using a specific 50 gene matrix). The classic hypothesis on monophyletic Mesozoa is possibly reborn and should be further tested.

Contrasting modes of mitochondrial genome evolution in sister taxa of wood-eating marine bivalves (Teredinidae and Xylophagaidae)

The bivalve families Teredinidae and Xylophagaidae include voracious consumers of wood in shallow and deep-water marine environments, respectively. The taxa are sister clades whose members consume wood as food with the aid of intracellular cellulolytic endosymbionts housed in their gills. This combination of adaptations is found in no other group of animals and was likely present in the common ancestor of both families. Despite these commonalities, the two families have followed dramatically different evolutionary paths with respect to anatomy, life history and distribution. Here we present 42 new mitochondrial genome sequences from Teredinidae and Xylophagaidae and show that distinct trajectories have also occurred in the evolution and organization of their mitochondrial genomes. Teredinidae display significantly greater rates of amino acid substitution but absolute conservation of protein-coding gene order, whereas Xylophagaidae display significantly less amino acid change but have undergone numerous and diverse changes in genome organization since their divergence from a common ancestor. As with many bivalves, these mitochondrial genomes encode two ribosomal RNAs, 12 protein coding genes, and 22 tRNAs; atp8 was not detected. We further show that their phylogeny, as inferred from amino acid sequences of 12 concatenated mitochondrial protein-coding genes, is largely congruent with those inferred from their nuclear genomes based on 18S and 28S ribosomal RNA sequences. Our results provide a robust phylogenetic framework to explore the tempo and mode of mitochondrial genome evolution and offer directions for future phylogenetic and taxonomic studies of wood-boring bivalves.

Molecular dating of the blood pigment hemocyanin provides new insight into the origin of animals

The Neoproterozoic included changes in oceanic redox conditions, the configuration of continents and climate, extreme ice ages (Sturtian and Marinoan), and the rise of complex life forms. A much-debated topic in geobiology concerns the influence of atmospheric oxygenation on Earth and the origin and diversification of animal lineages, with the most widely popularized hypotheses relying on causal links between oxygen levels and the rise of animals. The vast majority of extant animals use aerobic metabolism for growth and homeostasis; hence, the binding and transportation of oxygen represent a vital physiological task. Considering the blood pigment hemocyanin (Hc) is present in sponges and ctenophores, and likely to be present in the common ancestor of animals, we investigated the evolution and date of Hc emergence using bioinformatics approaches on both transcriptomic and genomic data. Bayesian molecular dating suggested that the ancestral animal Hc gene arose approximately 881 Ma during the Tonian Period (1000-720 Ma), prior to the extreme glaciation events of the Cryogenian Period (720-635 Ma). This result is corroborated by a recently discovered fossil of a putative sponge ~890 Ma and modern molecular dating for the origin of metazoans of ~1,000-650 Ma (but does contradict previous inferences regarding the origin of Hc ~700-600 Ma). Our data reveal that crown-group animals already possessed hemocyanin-like blood pigments, which may have enhanced the oxygen-carrying capacity of these animals in hypoxic environments at that time or acted in the transport of hormones, detoxification of heavy metals, and immunity pathways.

Metagenomics of Antarctic Marine Sediment Reveals Potential for Diverse Chemolithoautotrophy

The microbial biogeochemical processes occurring in marine sediment in Antarctica remain underexplored due to limited access. Further, these polar habitats are unique, as they are being exposed to significant changes in their climate. To explore how microbes drive biogeochemistry in these sediments, we performed a shotgun metagenomic survey of marine surficial sediment (0 to 3 cm of the seafloor) collected from 13 locations in western Antarctica and assembled 16 high-quality metagenome assembled genomes for focused interrogation of the lifestyles of some abundant lineages. We observe an abundance of genes from pathways for the utilization of reduced carbon, sulfur, and nitrogen sources. Although organotrophy is pervasive, nitrification and sulfide oxidation are the dominant lithotrophic pathways and likely fuel carbon fixation via the reverse tricarboxylic acid and Calvin cycles. Oxygen-dependent terminal oxidases are common, and genes for reduction of oxidized nitrogen are sporadically present in our samples. Our results suggest that the underlying benthic communities are well primed for the utilization of settling organic matter, which is consistent with findings from highly productive surface water. Despite the genetic potential for nitrate reduction, the net catabolic pathway in our samples remains aerobic respiration, likely coupled to the oxidation of sulfur and nitrogen imported from the highly productive Antarctic water column above. IMPORTANCE The impacts of climate change in polar regions, like Antarctica, have the potential to alter numerous ecosystems and biogeochemical cycles. Increasing temperature and freshwater runoff from melting ice can have profound impacts on the cycling of organic and inorganic nutrients between the pelagic and benthic ecosystems. Within the benthos, sediment microbial communities play a critical role in carbon mineralization and the cycles of essential nutrients like nitrogen and sulfur. Metagenomic data collected from sediment samples from the continental shelf of western Antarctica help to examine this unique system and document the metagenomic potential for lithotrophic metabolisms and the cycles of both nitrogen and sulfur, which support not only benthic microbes but also life in the pelagic zone.

Assessing genomic diversity, connectivity, and riverscape genetics hypotheses in the endangered Rough Hornsnail, Pleurocera foremani, following habitat disruption

The southeastern United States is home to some of the richest biodiversity in the world. Over the last 200 years, however, rapid industrialization and urbanization have threatened many natural areas, including freshwater habitats. River impoundments have also rapidly altered freshwater habitats, often resulting in species extirpation or extinction. The Coosa River in Alabama experienced one of the largest faunal declines in modern history after impoundment, making it an ideal system for studying how invertebrate species are affected by reservoir creation. One such species, the Rough Hornsnail, Pleurocera foremani, is an endangered freshwater snail in the family Pleuroceridae. We sampled all known localities of P. foremani and used 2bRAD-seq to measure genetic diversity. We assessed riverscape genomic patterns across the current range of P. foremani and measured gene flow within and between impoundments. We also investigated the degree to which P. foremani displays an isolation by distance pattern and conforms to broad hypotheses that have been put forth for population genetics of riverine species like the Mighty Headwater Hypothesis that predicts greater genetic diversity in headwater reaches compared with mainstem populations. Like most other freshwater species, a pattern of isolation by distance was observed in P. foremani. We also found that Coosa River dams are a barrier to gene flow, and genetic fragmentation of P. foremani is likely to increase. However, gene flow appeared common within reservoirs and tributaries. Additionally, we found that spatial genetic structure of P. foremani deviates from what is expected under the Mighty Headwaters Hypothesis, adding to a growing body of research suggesting that the majority of genetic diversity in low-dispersing gastropods is found in mainstem populations.

TIAMMAt: Leveraging biodiversity to revise protein domain models, evidence from innate immunity

Sequence annotation is fundamental for studying the evolution of protein families, particularly when working with nonmodel species. Given the rapid, ever-increasing number of species receiving high-quality genome sequencing, accurate domain modeling that is representative of species diversity is crucial for understanding protein family sequence evolution and their inferred function(s). Here, we describe a bioinformatic tool called Taxon-Informed Adjustment of Markov Model Attributes (TIAMMAt) which revises domain profile hidden Markov models (HMMs) by incorporating homologous domain sequences from underrepresented and nonmodel species. Using innate immunity pathways as a case study, we show that revising profile HMM parameters to directly account for variation in homologs among underrepresented species provides valuable insight into the evolution of protein families. Following adjustment by TIAMMAt, domain profile HMMs exhibit changes in their per-site amino acid state emission probabilities and insertion/deletion probabilities while maintaining the overall structure of the consensus sequence. Our results show that domain revision can heavily impact evolutionary interpretations for some families (i.e., NLR’s NACHT domain), whereas impact on other domains (e.g., rel homology domain and interferon regulatory factor domains) is minimal due to high levels of sequence conservation across the sampled phylogenetic depth (i.e., Metazoa). Importantly, TIAMMAt revises target domain models to reflect homologous sequence variation using the taxonomic distribution under consideration by the user. TIAMMAt’s flexibility to revise any subset of the Pfam database using a user-defined taxonomic pool will make it a valuable tool for future protein evolution studies, particularly when incorporating (or focusing) on nonmodel species.

Genome-wide characterization of LTR retrotransposons in the non-model deep-sea annelid Lamellibrachia luymesi

Background

Long Terminal Repeat retrotransposons (LTR retrotransposons) are mobile genetic elements composed of a few genes between terminal repeats and, in some cases, can comprise over half of a genome’s content. Available data on LTR retrotransposons have facilitated comparative studies and provided insight on genome evolution. However, data are biased to model systems and marine organisms, including annelids, have been underrepresented in transposable elements studies. Here, we focus on genome of Lamellibrachia luymesi, a vestimentiferan tubeworm from deep-sea hydrocarbon seeps, to gain knowledge of LTR retrotransposons in a deep-sea annelid.

Results

We characterized LTR retrotransposons present in the genome of L. luymesi using bioinformatic approaches and found that intact LTR retrotransposons makes up about 0.1% of L. luymesi genome. Previous characterization of the genome has shown that this tubeworm hosts several known LTR-retrotransposons. Here we describe and classify LTR retrotransposons in L. luymesi as within the Gypsy, Copia and Bel-pao superfamilies. Although, many elements fell within already recognized families (e.g., Mag, CSRN1), others formed clades distinct from previously recognized families within these superfamilies. However, approximately 19% (41) of recovered elements could not be classified. Gypsy elements were the most abundant while only 2 Copia and 2 Bel-pao elements were present. In addition, analysis of insertion times indicated that several LTR-retrotransposons were recently transposed into the genome of L. luymesi, these elements had identical LTR’s raising possibility of recent or ongoing retrotransposon activity.

Conclusions

Our analysis contributes to knowledge on diversity of LTR-retrotransposons in marine settings and also serves as an important step to assist our understanding of the potential role of retroelements in marine organisms. We find that many LTR retrotransposons, which have been inserted in the last few million years, are similar to those found in terrestrial model species. However, several new groups of LTR retrotransposons were discovered suggesting that the representation of LTR retrotransposons may be different in marine settings. Further study would improve understanding of the diversity of retrotransposons across animal groups and environments.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07749-1.

The impact of aquaculture on the genetics and distribution of the onuphid annelid Diopatra biscayensis

AIM

Evolutionary history of natural populations can be confounded by human intervention such as the case of decorator worm species Diopatra (Onuphidae), which have a history of being transported through anthropogenic activities. Because they build tubes and act as ecosystem engineers, they can have a large impact on the overall ecosystem in which they occur. One conspicuous member, Diopatra biscayensis, which was only described in 2012, has a fragmented distribution that includes the Bay of Biscay and the Normanno-Breton Gulf in the English Channel. This study explores the origin of these worms in the Normanno-Breton region, which has been debated to either be the result of a historic range contraction from a relic continuous population or a more recent introduction.

LOCATION

Northeastern Atlantic, the Bay of Biscay, and the Normanno-Breton Gulf.

METHODS

We utilized a RAD-tag-based SNP approach to create a reduced genomic data set to recover fine-scale population structure and infer which hypothesis best describes the D. biscayensis biogeographic distribution. The reduced genomic data set was used to calculate standard genetic diversities and genetic differentiation statistics, and utilized various clustering analyses, including PCAs, DAPC, and admixture.

RESULTS

Clustering analyses were consistent with D. biscayensis as a single population spanning the Bay of Biscay to the Normanno-Breton Gulf in the English Channel, although unexpected genetic substructure was recovered from Arcachon Bay, in the middle of its geographic range. Consistent with a hypothesized introduction, the isolated Sainte-Anne locality in the Normanno-Breton Gulf was recovered to be a subset of the diversity found in the rest of the Bay of Biscay.

MAIN CONCLUSIONS

These results are congruent with previous simulations that did not support connectivity from the Bay of Biscay to the Normanno-Breton Gulf by natural dispersal. These genomic findings, with support from previous climatic studies, further support the hypothesis that D. biscayensis phylogeographic connectivity is the result of introductions, likely through the regions' rich shellfish aquaculture, and not of a historically held range contraction.

The mitochondrial genome of the bone-eating worm Osedaxrubiplumus(Annelida, Siboglinidae)

Osedaxrubiplumus(Annelida, Siboglinidae)uses heterotrophic bacteria to feed onvertebrate carcasses and is currently found in the Pacific, Antarctic and Indian Ocean.Here, we report its nearly complete mitochondrial genomes assembled for 2 individuals, one from the East Pacific and the other from the Southwest Indian Ocean. Recoveredmitogenomes were 15591 and 15972 bp in length, with both consisting of 37 typical metazoan mitochondrial genes. All genes were transcribed from the same strand, and arranged in the same order as the other siboglinids, revealing conserved gene arrangement withinSiboglinidae. Phylogeneticanalysis of 13 protein coding genes confirms the placement of Osedaxsister to the Vestimentifera+Sclerolinum clade.

Evolutionary History of the Globin Gene Family in Annelids

Animals depend on the sequential oxidation of organic molecules to survive; thus, oxygen-carrying/transporting proteins play a fundamental role in aerobic metabolism. Globins are the most common and widespread group of respiratory proteins. They can be divided into three types: circulating intracellular, noncirculating intracellular, and extracellular, all of which have been reported in annelids. The diversity of oxygen transport proteins has been underestimated across metazoans. We probed 250 annelid transcriptomes in search of globin diversity in order to elucidate the evolutionary history of this gene family within this phylum. We report two new globin types in annelids, namely androglobins and cytoglobins. Although cytoglobins and myoglobins from vertebrates and from invertebrates are referred to by the same name, our data show they are not genuine orthologs. Our phylogenetic analyses show that extracellular globins from annelids are more closely related to extracellular globins from other metazoans than to the intracellular globins of annelids. Broadly, our findings indicate that multiple gene duplication and neo-functionalization events shaped the evolutionary history of the globin family.

Phylogenomics of Aplacophora (Mollusca, Aculifera) and a solenogaster without a foot

Recent molecular phylogenetic investigations strongly supported the placement of the shell-less, worm-shaped aplacophoran molluscs (Solenogastres and Caudofoveata) and chitons (Polyplacophora) in a clade called Aculifera, which is the sister taxon of all other molluscs. Thus, understanding the evolutionary history of aculiferan molluscs is important for understanding early molluscan evolution. In particular, fundamental questions about evolutionary relationships within Aplacophora have long been unanswered. Here, we supplemented the paucity of available data with transcriptomes from 25 aculiferans and conducted phylogenomic analyses on datasets with up to 525 genes and 75 914 amino acid positions. Our results indicate that aplacophoran taxonomy requires revision as several traditionally recognized groups are non-monophyletic. Most notably, Cavibelonia, the solenogaster taxon defined by hollow sclerites, is polyphyletic, suggesting parallel evolution of hollow sclerites in multiple lineages. Moreover, we describe Apodomenia enigmatica sp. nov., a bizarre new species that appears to be a morphological intermediate between Solenogastres and Caudofoveata. This animal is not a missing link, however; molecular and morphological studies show that it is a derived solenogaster that lacks a foot, mantle cavity and radula. Taken together, these results shed light on the evolutionary history of Aplacophora and reveal a surprising degree of morphological plasticity within the group.

Patterns of gene evolution following duplications and speciations in vertebrates

BACKGROUND

Eukaryotic genes typically form independent evolutionary lineages through either speciation or gene duplication events. Generally, gene copies resulting from speciation events (orthologs) are expected to maintain similarity over time with regard to sequence, structure and function. After a duplication event, however, resulting gene copies (paralogs) may experience a broader set of possible fates, including partial (subfunctionalization) or complete loss of function, as well as gain of new function (neofunctionalization). This assumption, known as the Ortholog Conjecture, is prevalent throughout molecular biology and notably plays an important role in many functional annotation methods. Unfortunately, studies that explicitly compare evolutionary processes between speciation and duplication events are rare and conflicting.

METHODS

To provide an empirical assessment of ortholog/paralog evolution, we estimated ratios of nonsynonymous to synonymous substitutions (ω = dN/dS) for 251,044 lineages in 6,244 gene trees across 77 vertebrate taxa.

RESULTS

Overall, we found ω to be more similar between lineages descended from speciation events (p < 0.001) than lineages descended from duplication events, providing strong support for the Ortholog Conjecture. The asymmetry in ω following duplication events appears to be largely driven by an increase along one of the paralogous lineages, while the other remains similar to the parent. This trend is commonly associated with neofunctionalization, suggesting that gene duplication is a significant mechanism for generating novel gene functions.

Genomic adaptations to chemosymbiosis in the deep-sea seep-dwelling tubeworm Lamellibrachia luymesi

BACKGROUND

Symbiotic relationships between microbes and their hosts are widespread and diverse, often providing protection or nutrients, and may be either obligate or facultative. However, the genetic mechanisms allowing organisms to maintain host-symbiont associations at the molecular level are still mostly unknown, and in the case of bacterial-animal associations, most genetic studies have focused on adaptations and mechanisms of the bacterial partner. The gutless tubeworms (Siboglinidae, Annelida) are obligate hosts of chemoautotrophic endosymbionts (except for Osedax which houses heterotrophic Oceanospirillales), which rely on the sulfide-oxidizing symbionts for nutrition and growth. Whereas several siboglinid endosymbiont genomes have been characterized, genomes of hosts and their adaptations to this symbiosis remain unexplored.

RESULTS

Here, we present and characterize adaptations of the cold seep-dwelling tubeworm Lamellibrachia luymesi, one of the longest-lived solitary invertebrates. We sequenced the worm's ~ 688-Mb haploid genome with an overall completeness of ~ 95% and discovered that L. luymesi lacks many genes essential in amino acid biosynthesis, obligating them to products provided by symbionts. Interestingly, the host is known to carry hydrogen sulfide to thiotrophic endosymbionts using hemoglobin. We also found an expansion of hemoglobin B1 genes, many of which possess a free cysteine residue which is hypothesized to function in sulfide binding. Contrary to previous analyses, the sulfide binding mediated by zinc ions is not conserved across tubeworms. Thus, the sulfide-binding mechanisms in sibgolinids need to be further explored, and B1 globins might play a more important role than previously thought. Our comparative analyses also suggest the Toll-like receptor pathway may be essential for tolerance/sensitivity to symbionts and pathogens. Several genes related to the worm's unique life history which are known to play important roles in apoptosis, cell proliferation, and aging were also identified. Last, molecular clock analyses based on phylogenomic data suggest modern siboglinid diversity originated in 267 mya (± 70 my) support previous hypotheses indicating a Late Mesozoic or Cenozoic origins of approximately 50-126 mya for vestimentiferans.

CONCLUSIONS

Here, we elucidate several specific adaptations along various molecular pathways that link phenome to genome to improve understanding of holobiont evolution. Our findings of adaptation in genomic mechanisms to reducing environments likely extend to other chemosynthetic symbiotic systems.

Neotropical Turtle Blood Flukes: Two New Genera and Species from the Amazon River Basin with a Key to Genera and Comments on a Marine-Derived Parasite Lineage in South America

Two new genera and species of freshwater turtle blood flukes (TBFs) are described herein based on specimens infecting the nephritic and mesenteric blood vessels of "matamatas" (a side-necked turtle, Chelus fimbriata [Schneider, 1783] [Pleurodira: Chelidae]) from the Amazon River Basin, Peru. These taxa comprise the first-named species and the first-proposed genera of freshwater TBFs from the continent of South America. A new comparison of all TBF genera produced 6 morphologically diagnosed groups that are discussed in light of previous TBF classification schemes and a novel phylogenetic hypothesis based on the nuclear large subunit ribosomal DNA (28S). Considering external and internal anatomical features, species of the new genera (Atamatam Bullard and Roberts n. gen., Paratamatam Bullard and Roberts n. gen.) are most similar to each other and are together most similar to those of several marine TBF genera. The 28S phylogenetic analysis supported the monophyly of all 6 morphologically diagnosed groups of genera. Most notably, the freshwater TBFs of South America comprise a derived group nested within the clade that includes the paraphyletic marine TBFs. Not surprisingly in light of morphology, another marine TBF lineage (Neospirorchis Price, 1934) clustered with the freshwater TBFs of Baracktrema Roberts, Platt, and Bullard, 2016 and Unicaecum Stunkard, 1925. Our results, including an ancestral state reconstruction, indicated that (1) freshwater TBFs have colonized marine turtles twice independently and that (2) the South American freshwater TBFs comprise a marine-derived lineage. This is the first evidence that TBFs have twice independently transitioned from a marine to freshwater definitive host. Marine incursion is considered as a possible mechanism affecting the natural history of marine-derived freshwater TBFs in South America. A dichotomous key to accepted TBF genera is provided.

Sequencing Disparity in the Genomic Era

Advances in sequencing technology have resulted in the expectation that genomic studies will become more representative of organismal diversity. To test this expectation, we explored species representation of nonhuman eukaryotes in the Sequence Read Archive. Though species richness has been increasing steadily, species evenness is decreasing over time. Moreover, the top 1% most studied organisms increasingly represent a larger proportion of total experiments, demonstrating growing bias in favor of a small minority of species. To better understand molecular processes and patterns, genomic studies should reverse current trends by adopting more comparative approaches.

Riverscape genetic variation, migration patterns, and morphological variation of the threatened Round Rocksnail, Leptoxis ampla

Within riverine systems, headwater populations are hypothesized to harbour higher amounts of genetic distinctiveness than populations in the main stem of a river and display increased genetic diversity in large, downstream habitats. However, these hypotheses were mostly developed with insects and fish, and they have not been tested on many invertebrate lineages. Pleuroceridae gastropods are of particular ecological importance to rivers of eastern North America, sometimes comprising over 90% of macroinvertebrate biomass. Yet, virtually nothing is known of pleurocerid landscape genetics, including whether genetic diversity follows predictions made by hypotheses developed on more mobile species. Moreover, the commonly repeated hypothesis that intraspecific morphological variation in gastropods results from ecophenotypic plasticity has not been well tested on pleurocerids. Using 2bRAD-seq to discover single nucleotide polymorphisms, we show that the threatened, Cahaba River endemic pleurocerid, Leptoxis ampla, has limited gene flow among populations and that migration is downstream-biased, conflicting with previous hypotheses. Both tributary and main stem populations harbour unique genomic profiles, and genetic diversity was highest in downstream populations. Furthermore, L. ampla shell morphology was more correlated with genetic differences among individuals and populations than habitat characteristics. We anticipate similar genetic and demographic patterns to be seen in other pleurocerids, and hypotheses about gene flow and population demographics that were based on more mobile taxa often, but not always, apply to freshwater gastropods. From a conservation standpoint, genetic structure of L. ampla populations suggests distinctive genetic diversity is lost with localized extirpation, a phenomenon common across the range of Pleuroceridae.

Newly Discovered Occurrences and Gene Tree of the Extracellular Globins and Linker Chains from the Giant Hexagonal Bilayer Hemoglobin in Metazoans

Multicellular organisms depend on oxygen-carrying proteins to transport oxygen throughout the body; therefore, proteins such as hemoglobins (Hbs), hemocyanins, and hemerythrins are essential for maintenance of tissues and cellular respiration. Vertebrate Hbs are among the most extensively studied proteins; however, much less is known about invertebrate Hbs. Recent studies of hemocyanins and hemerythrins have demonstrated that they have much wider distributions than previously thought, suggesting that oxygen-binding protein diversity is underestimated across metazoans. Hexagonal bilayer hemoglobin (HBL-Hb), a blood pigment found exclusively in annelids, is a polymer comprised up to 144 extracellular globins and 36 linker chains. To further understand the evolutionary history of this protein complex, we explored the diversity of linkers and extracellular globins from HBL-Hbs using in silico approaches on 319 metazoan and one choanoflagellate transcriptomes. We found 559 extracellular globin and 414 linker genes transcribed in 171 species from ten animal phyla with new records in Echinodermata, Hemichordata, Brachiopoda, Mollusca, Nemertea, Bryozoa, Phoronida, Platyhelminthes, and Priapulida. Contrary to previous suggestions that linkers and extracellular globins emerged in the annelid ancestor, our findings indicate that they have putatively emerged before the protostome-deuterostome split. For the first time, we unveiled the comprehensive evolutionary history of metazoan HBL-Hb components, which consists of multiple episodes of gene gains and losses. Moreover, because our study design surveyed linkers and extracellular globins independently, we were able to cross-validate our results, significantly reducing the rate of false positives. We confirmed that the distribution of HBL-Hb components has until now been underestimated among animals.

Mitogenomics Reveals a Novel Genetic Code in Hemichordata

The diverse array of codon reassignments demonstrate that the genetic code is not universal in nature. Exploring mechanisms underlying codon reassignment is critical for understanding the evolution of the genetic code during translation. Hemichordata, comprising worm-like Enteropneusta and colonial filter-feeding Pterobranchia, is the sister taxon of echinoderms and is more distantly related to chordates. However, only a few hemichordate mitochondrial genomes have been sequenced, hindering our understanding of mitochondrial genome evolution within Deuterostomia. In this study, we sequenced four mitochondrial genomes and two transcriptomes, including representatives of both major hemichordate lineages and analyzed together with public available data. Contrary to the current understanding of the mitochondrial genetic code in hemichordates, our comparative analyses suggest that UAA encodes Tyr instead of a "Stop" codon in the pterobranch lineage Cephalodiscidae. We also predict that AAA encodes Lys in pterobranch and enteropneust mitochondrial genomes, contradicting the previous assumption that hemichordates share the same genetic code with echinoderms for which AAA encodes Asn. Thus, we propose a new mitochondrial genetic code for Cephalodiscus and a revised code for enteropneusts. Moreover, our phylogenetic analyses are largely consistent with previous phylogenomic studies. The only exception is the phylogenetic position of the enteropneust Stereobalanus, whose placement as sister to all other described enteropneusts. With broader taxonomic sampling, we provide evidence that evolution of mitochondrial gene order and genetic codes in Hemichordata are more dynamic than previously thought and these findings provide insights into mitochondrial genome evolution within this clade.

Discovery of Novel Hemocyanin-Like Genes in Metazoans

Among animals, two major groups of oxygen-binding proteins are found: proteins that use iron to bind oxygen (hemoglobins and hemerythrins) and two non-homologous hemocyanins that use copper. Although arthropod and mollusc hemocyanins bind oxygen in the same manner, they are distinct in their molecular structures. In order to better understand the range of natural variation in hemocyanins, we searched for them in a diverse array of metazoan transcriptomes by using bioinformatics tools to examine hemocyanin evolutionary history and to consequently revive the discussion about whether all metazoan hemocyanins shared a common origin with frequent losses or whether they originated separately after the divergence of Lophotrochozoa and Ecdysozoa. We confirm that the distribution of hemocyanin-like genes is more widespread than previously reported, including five putative novel mollusc hemocyanin genes in two annelid species from Chaetopteridae. For arthropod hemocyanins, 16 putative novel genes were retained, and the presence of arthropod hemocyanins in 11 annelid species represents a novel observation. Interestingly, Annelida is the lineage that presents the greatest repertoire of oxygen transport proteins reported to date, possessing all the main superfamily proteins, which could be explained partially by the immense variability of lifestyles and habitats. Work presented here contradicts the canonical view that hemocyanins are restricted to molluscs and arthropods, suggesting that the occurrence of copper-based blood pigments in metazoans has been underestimated. Our results also support the idea of the presence of oxygen carrier hemocyanins being widespread across metazoans with an evolutionary history characterized by frequent losses.

Challenging Dogma Concerning Biogeographic Patterns of Antarctica and the Southern Ocean

Antarctica is enormous, cold, remote, and particularly sensitive to climate change. Most biological research below 60°S has focused on the isolated nature of the biota and how organisms have adapted to the cold and ice. However, biogeographic patterns in Antarctica and the Southern Ocean, and the processes explaining how those patterns came about, still await adequate explanation. Both terrestrial and marine organisms have been influenced by climatic change (e.g., glaciation), physical phenomena (e.g., oceanic currents), and/or potential barriers to gene flow (e.g., steep thermal gradients). Whereas the Antarctic region contains diverse and complex marine communities, terrestrial systems tend to be comparatively simple with limited diversity. Here, we challenge the current dogma used to explain the diversity and biogeographic patterns present in the Antarctic. We assert that relatively modern processes within the last few million years, rather than geo-logical events that occurred in the Eocene and Miocene, account for present patterns of biodiversity in the region. Additionally, reproductive life history stages appear to have little influence in structuring genetic patterns in the Antarctic, as currents and glacial patterns are noted to be more important drivers of organismal patterns of distribution. Finally, we highlight the need for additional sampling, high-throughput genomic approaches, and broad, multinational cooperation for addressing outstanding questions of Antarctic biogeography and biodiversity.

River Flow Impacts Bacterial and Archaeal Community Structure in Surface Sediments in the Northern Gulf of Mexico

Meiobenthic community structure in the northern Gulf of Mexico has been shown to be driven by geographical differences due to inshore-offshore gradients and location relative to river discharge. Samples collected along three transects spanning Mobile Bay, Alabama, showed significant differences in meiobenthic communities east of the bay compared to those sampled from the west. In contrast, analysis of bacterial and archaeal communities from the same sediment samples shows that the inshore-offshore gradient has minimal impact on their community structure. Significant differences in community structure were observed for Bacteria and Archaea between the east and west samples, but there was no difference in richness or diversity. Grouped by sediment type, higher richness was observed in silty samples compared to sandy samples. Significant differences were also observed among sediment types for community structure with bacteria communities in silty samples having more anaerobic sulfate reducers compared to aerobic heterotrophs, which had higher abundances in sandy sediments. This is likely due to increased organic matter in the silty sediments from the overlying river leading to low oxygen habitats. Most archaeal sequences represented poorly characterized high-level taxa, limiting interpretation of their distributions. Overlap between groups based on transect and sediment characteristics made determining which factor is more important in structuring bacterial and archaeal communities difficult. However, both factors are driven by discharge from the Mobile River. Although inshore-offshore gradients do not affect Bacteria or Archaea to the same extent as the meiobenthic communities, all three groups are strongly affected by sediment characteristics.

Conservation of mitochondrial genome arrangements in brittle stars (Echinodermata, Ophiuroidea)

Brittle stars are conspicuous members of benthic ecosystems, fill many ecological niches and are the most speciose of all classes of echinoderms. With high levels of biodiversity, elucidating the evolutionary history of this group is important. Understanding of higher-level relationships within Ophiuroidea has been aided by multilocus nuclear data and DNA barcoding. However, the degree of consistency between mitochondrial and nuclear data within ophiuroids remains unclear and deserves further assessment. In this study, 17 mitochondrial genomes spanning the taxonomic breadth of Ophiuroidea were utilized to explore evolutionary relationships through maximum likelihood analyses, Bayesian inference and comparative assessment of gene order. Our phylogenetic analyses, based on both nucleotide and amino acid residues, support recent findings based on multilocus nuclear data and morphology, in that the brittle star clades Ophintegrida and Euryophiurida were recovered as monophyletic with the latter comprising Euyalida, Ophiuridae and Ophiopyrgidae. Only three different arrangements of the 13 protein coding and 2 ribosomal RNA genes were observed. As expected, tRNA genes were more likely to have undergone rearrangement but the order of all 37 genes was found to be conserved in all sampled Euryalida and Ophiuridae. Both Euryalida and the clade comprised of Ophiuridae and Ophiopyrgidae, each had their own conserved rearrangement of protein coding genes and ribosomal genes, after divergence from their last common ancestor. Euryalida has a rearrangement of the two ribosomal RNA genes, rrnS and rrnL, in contrast to Ophiuridae and Ophiopyrgidae, which had an inversion of the genes nad1, nad2, and cob relative to Ophintegrida. Further, our data support the gene order found in all sampled Euryalida as the most likely ancestral order for all Ophiuroidea.

Population Genomics of Nymphon australe Hodgson, 1902 (Pycnogonida, Nymphonidae) in the Western Antarctic

Within the Southern Ocean, the Antarctic Circumpolar Current is hypothesized to facilitate a circumpolar distribution for many taxa, even though some, such as pycnogonids, are assumed to have limited ability to disperse, based on brooding life histories and adult ambulatory capabilities. With a number of contradictions to circumpolarity reported in the literature for other pycnogonids, alternative hypotheses have been explored, particularly for Nymphon australe, the most common species of Pycnogonida (sea spider) in the Southern Ocean. Glacial events have been hypothesized to impact the capacity of organisms to colonize suitable areas without ice coverage as refuge and without the eurybathic capacity to colonize deeper areas. In this study, we examine populations of one presumed circumpolar species, the pycnogonid N. australe, from throughout the Western Antarctic, using a 2b-RAD approach to detect genetic variation with single-nucleotide polymorphisms. Using this approach, we found that N. australe included two distinct groups from within >5000-km sampling region. By using a discriminant analysis of principle components, sparse nonnegative matrix factorization, and admixture coefficient analysis, two distinctive populations were revealed in the Western Antarctic: one covered distances greater than 5000 km (Weddell, Western Antarctic Peninsula, and Ross Sea), and the other shared limited connectivity entrained within the Amundsen Sea. Under further scrutiny of the 3086 single-nucleotide polymorphisms in the data set, only 78 loci had alignment stacks between the two populations. We propose that the populations analyzed are divergent enough to constitute two different species from within this common Antarctic genus known for its phenotypic plasticity.

Mitogenomics reveals phylogenetic relationships of caudofoveate aplacophoran molluscs

The worm-shaped, shell-less aplacophoran molluscs Caudofoveata and Solenogastres have recently received attention as part of the clade Aculifera, but relationships within these two lineages are still largely unknown. Here, we use complete mitochondrial genomes to shed light on higher-level relationships within Caudofoveata. Mitochondrial genomes have been sequenced for many diverse molluscs, but only two mitochondrial genomes from aplacophoran molluscs (the caudofoveates Scutopus ventrolineatus and Chaetoderma nitidulum) are available to date. We sequenced and assembled complete or near complete mitochondrial genomes of five additional species of Caudofoveata (Falcidens acutargatus, Falcidens halanychi, Scutopus robustus, Psilodens balduri and Spathoderma clenchi) and one species of Solenogastres (Neomenia carinata) for comparison to available mitochondrial genomes of aculiferans. Comparison of mitochondrial gene order among different lineages revealed a highly conserved order of protein coding genes corresponding to the hypothesized ancestral gene order for Mollusca. Unique arrangements of tRNAs were found among lineages of aculiferan molluscs as well as among caudofoveate taxa. Phylogenetic analyses of amino acid sequences for the 13 protein-coding genes recovered a monophyletic Aplacophora. Within Caudofoveata, Chaetodermatidae, but not Limifossoridae, was recovered monophyletic. Taken together, our results suggest that mitochondrial genomes can serve as useful molecular markers for aculiferan phylogenetics, especially for more recent phylogenetic events.

Broad Phylogenetic Occurrence of the Oxygen-Binding Hemerythrins in Bilaterians

Animal tissues need to be properly oxygenated for carrying out catabolic respiration and, as such, natural selection has presumably favored special molecules that can reversibly bind and transport oxygen. Hemoglobins, hemocyanins, and hemerythrins (Hrs) fulfill this role, with Hrs being the least studied. Knowledge of oxygen-binding proteins is crucial for understanding animal physiology. Hr genes are present in the three domains of life, Archaea, Bacteria, and Eukaryota; however, within Animalia, Hrs has been reported only in marine species in six phyla (Annelida, Brachiopoda, Priapulida, Bryozoa, Cnidaria, and Arthropoda). Given this observed Hr distribution, whether all metazoan Hrs share a common origin is circumspect. We investigated Hr diversity and evolution in metazoans, by employing in silico approaches to survey for Hrs from of 120 metazoan transcriptomes and genomes. We found 58 candidate Hr genes actively transcribed in 36 species distributed in 11 animal phyla, with new records in Echinodermata, Hemichordata, Mollusca, Nemertea, Phoronida, and Platyhelminthes. Moreover, we found that “Hrs” reported from Cnidaria and Arthropoda were not consistent with that of other metazoan Hrs. Contrary to previous suggestions that Hr genes were absent in deuterostomes, we find Hr genes present in deuterostomes and were likely present in early bilaterians, but not in nonbilaterian animal lineages. As expected, the Hr gene tree did not mirror metazoan phylogeny, suggesting that Hrs evolutionary history was complex and besides the oxygen carrying capacity, the drivers of Hr evolution may also consist of secondary functional specializations of the proteins, like immunological functions.

Who Let the CAT Out of the Bag? Accurately Dealing with Substitutional Heterogeneity in Phylogenomic Analyses

As phylogenetic datasets have increased in size, site-heterogeneous substitution models such as CAT-F81 and CAT-GTR have been advocated in favor of other models because they purportedly suppress long-branch attraction (LBA). These models are two of the most commonly used models in phylogenomics, and they have been applied to a variety of taxa, ranging from Drosophila to land plants. However, many arguments in favor of CAT models have been based on tenuous assumptions about the true phylogeny, rather than rigorous testing with known trees via simulation. Moreover, CAT models have not been compared to other approaches for handling substitutional heterogeneity such as data partitioning with site-homogeneous substitution models. We simulated amino acid sequence datasets with substitutional heterogeneity on a variety of tree shapes including those susceptible to LBA. Data were analyzed with both CAT models and partitioning to explore model performance; in total over 670,000 CPU hours were used, of which over 97% was spent running analyses with CAT models. In many cases, all models recovered branching patterns that were identical to the known tree. However, CAT-F81 consistently performed worse than other models in inferring the correct branching patterns, and both CAT models often overestimated substitutional heterogeneity. Additionally, reanalysis of two empirical metazoan datasets supports the notion that CAT-F81 tends to recover less accurate trees than data partitioning and CAT-GTR. Given these results, we conclude that partitioning and CAT-GTR perform similarly in recovering accurate branching patterns. However, computation time can be orders of magnitude less for data partitioning, with commonly used implementations of CAT-GTR often failing to reach completion in a reasonable time frame (i.e., for Bayesian analyses to converge). Practices such as removing constant sites and parsimony uninformative characters, or using CAT-F81 when CAT-GTR is deemed too computationally expensive, cannot be logically justified. Given clear problems with CAT-F81, phylogenies previously inferred with this model should be reassessed. [Data partitioning; phylogenomics, simulation, site-heterogeneity, substitution models.].

Phylogenomics of Lophotrochozoa with Consideration of Systematic Error

Phylogenomic studies have improved understanding of deep metazoan phylogeny and show promise for resolving incongruences among analyses based on limited numbers of loci. One region of the animal tree that has been especially difficult to resolve, even with phylogenomic approaches, is relationships within Lophotrochozoa (the animal clade that includes molluscs, annelids, and flatworms among others). Lack of resolution in phylogenomic analyses could be due to insufficient phylogenetic signal, limitations in taxon and/or gene sampling, or systematic error. Here, we investigated why lophotrochozoan phylogeny has been such a difficult question to answer by identifying and reducing sources of systematic error. We supplemented existing data with 32 new transcriptomes spanning the diversity of Lophotrochozoa and constructed a new set of Lophotrochozoa-specific core orthologs. Of these, 638 orthologous groups (OGs) passed strict screening for paralogy using a tree-based approach. In order to reduce possible sources of systematic error, we calculated branch-length heterogeneity, evolutionary rate, percent missing data, compositional bias, and saturation for each OG and analyzed increasingly stricter subsets of only the most stringent (best) OGs for these five variables. Principal component analysis of the values for each factor examined for each OG revealed that compositional heterogeneity and average patristic distance contributed most to the variance observed along the first principal component while branch-length heterogeneity and, to a lesser extent, saturation contributed most to the variance observed along the second. Missing data did not strongly contribute to either. Additional sensitivity analyses examined effects of removing taxa with heterogeneous branch lengths, large amounts of missing data, and compositional heterogeneity. Although our analyses do not unambiguously resolve lophotrochozoan phylogeny, we advance the field by reducing the list of viable hypotheses. Moreover, our systematic approach for dissection of phylogenomic data can be applied to explore sources of incongruence and poor support in any phylogenomic data set. [Annelida; Brachiopoda; Bryozoa; Entoprocta; Mollusca; Nemertea; Phoronida; Platyzoa; Polyzoa; Spiralia; Trochozoa.].

Mitochondrial genome of Dinophilus gyrociliatus (Annelida: Dinophilidae)

Here we report the 14,678 bp mitochondrial genome of the annelid Dinophilus gyrociliatus, the first mitochondrial genome from Dinophilidae. We recovered 13 protein-coding genes, two rRNA, and 21 tRNA, the order of which is different from other annelid species. Interestingly, trnS1 was not recovered. The GC% across the genome was 34.20%.

Ctenophore relationships and their placement as the sister group to all other animals

Ctenophora, compromising approximately 200 described species, is an important lineage for understanding metazoan evolution and is of great ecological and economic importance. Ctenophore diversity includes species with unique colloblasts used for prey capture, smooth and striated muscles, benthic and pelagic lifestyles, and locomotion with ciliated paddles or muscular propulsion. However, ancestral states of traits are debated and relationships among many lineages are unresolved. Here, using 27 newly sequenced ctenophore transcriptomes, publicly available data, and methods to control systematic error we establish the placement of Ctenophora as the sister group to all other animals and refine phylogenetic relationships within ctenophores. Molecular clock analyses suggest modern ctenophore diversity originated approximately 350MYA ± 88 MY, conflicting with previous hypotheses of approximately 65 MYA. We recover Euplokamis dunlapae, a species with striated muscles, as the sister lineage to other sampled ctenophores. Ancestral state reconstruction shows the most recent common ancestor of extant ctenophores was pelagic, possessed tentacles, was bioluminescent, and did not have separate sexes. Our results imply at least two transitions from a pelagic to a benthic lifestyle within Ctenophora, suggesting such transitions were more common in animal diversification than appreciated.

Emendation and new species of Hapalorhynchus Stunkard, 1922 (Digenea: Schistosomatoidea) from musk turtles (Kinosternidae: Sternotherus) in Alabama and Florida rivers

Hapalorhynchus Stunkard, 1922 is emended based on morphological study of existing museum specimens (type and voucher specimens) and newly-collected specimens infecting musk turtles (Testudines: Kinosternidae: Sternotherus spp.) from rivers in Alabama and Florida (USA). Hapalorhynchus conecuhensis n. sp. is described from an innominate musk turtle, Sternotherus cf. minor, (type host) from Blue Spring (31°5'27.64″N, 86°30'53.21″W; Pensacola Bay Basin, Alabama) and the loggerhead musk turtle, Sternotherus minor (Agassiz, 1857) from the Wacissa River (30°20'24.73″N, 83°59'27.56″W; Apalachee Bay Basin, Florida). It differs from congeners by lacking a body constriction at level of the ventral sucker, paired anterior caeca, and a transverse ovary as well as by having a small ventral sucker, proportionally short posterior caeca, nearly equally-sized anterior and posterior testes, a small cirrus sac, and a uterus extending dorsal to the ovary and the anterior testis. Specimens of Hapalorhynchus reelfooti Byrd, 1939 infected loggerhead musk turtles, stripe-necked musk turtles (Sternotherus peltifer Smith and Glass, 1947), Eastern musk turtles (Sternotherus odoratus [Latreille in Sonnini and Latreille, 1801]), and S. cf. minor. Those of Hapalorhynchus cf. stunkardi infected S. minor and S. odoratus. Sternothorus minor, S. peltifer, and S. cf. minor plus S. minor and S. odoratus are new host records for H. reelfooti and H. cf. stunkardi, respectively. This is the first report of an infected musk turtle from the Coosa and Tallapoosa Rivers (Mobile-Tensaw River Basin), Pensacola Bay Basin, or Apalachee Bay Basin. Sequence analysis of the large subunit rDNA (28S) showed a strongly-supported clade for Hapalorhynchus.