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. Although 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 double-stranded RNA 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., nuclear factor κB and interferon regulatory factor 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.

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.

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.

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.

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.

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.

Mitogenomics reveals phylogenetic relationships of Arcoida (Mollusca, Bivalvia) and multiple independent expansions and contractions in mitochondrial genome size.

Arcoida, comprising about 570 species of blood cockles, is an ecologically and economically important lineage of bivalve molluscs. Current classification of arcoids is largely based on morphology, which shows widespread homoplasy. Despite two recent studies employing multi-locus analyses with broad sampling of Arcoida, evolutionary relationships among major lineages remain controversial. Interestingly, mitochondrial genomes of several ark shell species are 2-3 times larger than those found in most bilaterians, and are among the largest bilaterian mitochondrial genomes reported to date. These results highlight the need of detailed phylogenetic study to explore evolutionary relationships within Arcoida so that the evolution of mitochondrial genome size can be understood. To this end, we sequenced 17 mitochondrial genomes and compared them with publicly available data, including those from other lineages of Arcoida with emphasis on the subclade Arcoidea species. Our phylogenetic analyses indicate that Noetiidae, Cucullaeidae and Glycymerididae are nested within a polyphyletic Arcidae. Moreover, we find multiple independent expansions and potential contractions of mitochondrial genome size, suggesting that the large mitochondrial genome is not a shared ancestral feature in Arcoida. We also examined tandem repeats and inverted repeats in non-coding regions and investigated the presence of such repeats with relation to genome size variation. Our results suggest that tandem repeats might facilitate intraspecific mitochondrial genome size variation, and that inverted repeats, which could be derived from transposons, might be responsible for mitochondrial genome expansions and contractions. We show that mitochondrial genome size in Arcoida is more dynamic than previously understood and provide insights into evolution of mitochondrial genome size variation in metazoans.

The ctenophore genome and the evolutionary origins of neural systems.

The origins of neural systems remain unresolved. In contrast to other basal metazoans, ctenophores (comb jellies) have both complex nervous and mesoderm-derived muscular systems. These holoplanktonic predators also have sophisticated ciliated locomotion, behaviour and distinct development. Here we present the draft genome of Pleurobrachia bachei, Pacific sea gooseberry, together with ten other ctenophore transcriptomes, and show that they are remarkably distinct from other animal genomes in their content of neurogenic, immune and developmental genes. Our integrative analyses place Ctenophora as the earliest lineage within Metazoa. This hypothesis is supported by comparative analysis of multiple gene families, including the apparent absence of HOX genes, canonical microRNA machinery, and reduced immune complement in ctenophores. Although two distinct nervous systems are well recognized in ctenophores, many bilaterian neuron-specific genes and genes of 'classical' neurotransmitter pathways either are absent or, if present, are not expressed in neurons. Our metabolomic and physiological data are consistent with the hypothesis that ctenophore neural systems, and possibly muscle specification, evolved independently from those in other animals.

Toll-like receptor pathway evolution in deuterostomes.

Animals have evolved an array of pattern-recognition receptor families essential for recognizing conserved molecular motifs characteristic of pathogenic microbes. One such family is the Toll-like receptors (TLRs). On pathogen binding, TLRs initiate specialized cytokine signaling catered to the class of invading pathogen. This signaling is pivotal for activating adaptive immunity in vertebrates, suggesting a close evolutionary relationship between innate and adaptive immune systems. Despite significant advances toward understanding TLR-facilitated immunity in vertebrates, knowledge of TLR pathway evolution in other deuterostomes is limited. By analyzing genomes and transcriptomes across 37 deuterostome taxa, we shed light on the evolution and diversity of TLR pathway signaling elements. Here, we show that the deuterostome ancestor possessed a molecular toolkit homologous to that which drives canonical MYD88-dependent TLR signaling in contemporary mammalian lineages. We also provide evidence that TLR3-facilitated antiviral signaling predates the origin of its TCAM1 dependence recognized in the vertebrates. SARM1, a negative regulator of TCAM1-dependent pathways in vertebrates, was also found to be present across all major deuterostome lineages despite the apparent absence of TCAM1 in invertebrate deuterostomes. Whether the presence of SARM1 is the result of its role in immunity regulation, neuron physiology, or a function of both is unclear. Additionally, Bayesian phylogenetic analyses corroborate several lineage-specific TLR gene expansions in urchins and cephalochordates. Importantly, our results underscore the need to sample across taxonomic groups to understand evolutionary patterns of the innate immunity foundation on which complex immunological novelties arose.

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.

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.

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.

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.

Molecular phylogeny of Caudofoveata (Mollusca) challenges traditional views.

The shell-less, worm-shaped Caudofoveata (=Chaetodermomorpha) is one of the least known groups of molluscs. The taxon consists of 141 recognized species found from intertidal environments to the deep-sea where they live burrowing in sediment. Evolutionary relationships of the group have been debated, but few studies based on morphological or molecular data have investigated the phylogeny of the group. Here we use molecular phylogenetics to resolve relationships among and within families of Caudofoveata. Phylogenetic analyses were performed using selected mitochondrial and nuclear genes from species from all recognized families of Caudofoveata. In resulting trees and contrary to traditional views, Prochaetodermatidae forms the sister clade to a clade containing the other two currently recognized families, Chaetodermatidae and Limifossoridae. The monophyly of Prochaetodermatidae is highly supported, but Limifossoridae and Chaetodermatidae are not recovered as monophyletic. Most of the caudofoveate genera are also not recovered as monophyletic in our analyses. Thus results from our molecular data suggest that the current classification of Caudofoveata is in need of revision, and indicate evolutionary scenarios that differ from previously proposed hypotheses based on morphology.

Na+ /K+ -ATPase gene duplications in clitellate annelids are associated with freshwater colonization.

Major habitat transitions, such as those from marine to freshwater habitats or from aquatic to terrestrial habitats, have occurred infrequently in animal evolution and may represent a barrier to diversification. Identifying genomic events associated with these transitions can help us better understand mechanisms that allow animals to cross these barriers and diversify in new habitats. Study of the Capitella telata and Helobdella robusta genomes allows examination of one such habitat transition (marine to freshwater) in Annelida. Initial examination of these genomes indicated that the freshwater leech H. robusta contains many more copies (12) of the sodium-potassium pump alpha-subunit (Na+ /K+ -ATPase) gene than does the marine polychaete C. telata (2). The sodium-potassium pump plays a key role in maintenance of cellular ionic balance and osmoregulation, and Na+ /K+ -ATPase duplications may have helped annelids invade and diversify in freshwater habitats. To assess whether the timing of Na+ /K+ -ATPase duplications coincided with the marine-to-freshwater transition in Clitellata, we used transcriptomic data from 18 annelid taxa, along with the two genomes, to infer a species phylogeny and identified Na+ /K+ -ATPase gene transcripts in order to infer the timing of gene duplication events using tree-based methods. The inferred timing of Na+ /K+ -ATPase duplication events is consistent with the timing of the initial marine-to-freshwater transition early in the history of clitellate annelids, supporting the hypothesis that gene duplications may have played a role in the annelid diversification into freshwater habitats.

Spirorchis spp. (Digenea: Schistosomatoidea) infecting map turtles (Cryptodira: Emydidae: Graptemys spp.) in southeastern North America: A new species, molecular phylogenies, and key to species.

Black-knobbed map turtles (Graptemys nigrinoda Cagle) and Alabama map turtles (Graptemys pulchra Baur) were infected with several blood flukes in Alabama (southeastern North America). Spirorchis paraminutus Roberts & Bullard n. sp. differs from its congeners by having a body that is 12-24× longer than wide, a testicular column of 10 testes that is 1/5-1/4 of the body length and located far posterior to the caecal bifurcation (the anterior-most testis is located in the posterior body half), and a common genital pore that is ventral to the ovary and 1/4-1/3 of the body length from the posterior extremity. These turtles and an Escambia map turtle (Graptemys ernsti Lovich & McCoy) were infected with Spirorchis elegans Stunkard, 1923, Spirorchis scripta Stunkard, 1923 and two innominate species of Spirorchis MacCallum, 1918. Phylogenetic analyses of the internal transcribed spacer region (ITS2) and large subunit rDNA (28S) recovered a monophyletic Spirorchis and the new species sister to Spirorchis collinsi Roberts & Bullard, 2016.

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.

Phylogenomics offers resolution of major tunicate relationships.

Tunicata, a diverse clade of approximately 3000 described species of marine, filter-feeding chordates, is of great interest to researchers because tunicates are the closest living relatives of vertebrates and they facilitate comparative studies of our own biology. The group also includes numerous invasive species that cause considerable economic damage and some species of tunicates are edible. Despite their diversity and importance, relationships among major lineages of Tunicata are not completely resolved. Here, we supplemented public data with transcriptomes from seven species spanning the diversity of Tunicata and conducted phylogenomic analyses on data sets of up to 798 genes. Sensitivity analyses were employed to examine the influences of reducing compositional heterogeneity and branch-length heterogeneity. All analyses maximally supported a monophyletic Tunicata within Olfactores (Vertebrata + Tunicata). Within Tunicata, all analyses recovered Appendicularia sister to the rest of Tunicata and confirmed (with maximal support) that Thaliacea is nested within Ascidiacea. Stolidobranchia is the sister taxon to all other tunicates except Appendicularia. In most analyses, phlebobranch tunicates were recovered paraphyletic with respect to Aplousobranchia. Support for this topology varied but was strong in some cases. However, when only the 50 best genes based on compositional heterogeneity were analysed, we recovered Phlebobranchia and Aplousobranchia reciprocally monophyletic with strong support, consistent with most traditional morphology-based hypotheses. Examination of internode certainty also cast doubt on results of phlebobranch paraphyly, which may be due to limited taxon sampling. Taken together, these results provide a higher-level phylogenetic framework for our closest living invertebrate relatives.

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