A region of suppressed recombination misleads neoavian phylogenomics.

Genomes are typically mosaics of regions with different evolutionary histories. When speciation events are closely spaced in time, recombination makes the regions sharing the same history small, and the evolutionary history changes rapidly as we move along the genome. When examining rapid radiations such as the early diversification of Neoaves 66 Mya, typically no consistent history is observed across segments exceeding kilobases of the genome. Here, we report an exception. We found that a 21-Mb region in avian genomes, mapped to chicken chromosome 4, shows an extremely strong and discordance-free signal for a history different from that of the inferred species tree. Such a strong discordance-free signal, indicative of suppressed recombination across many millions of base pairs, is not observed elsewhere in the genome for any deep avian relationships. Although long regions with suppressed recombination have been documented in recently diverged species, our results pertain to relationships dating circa 65 Mya. We provide evidence that this strong signal may be due to an ancient rearrangement that blocked recombination and remained polymorphic for several million years prior to fixation. We show that the presence of this region has misled previous phylogenomic efforts with lower taxon sampling, showing the interplay between taxon and locus sampling. We predict that similar ancient rearrangements may confound phylogenetic analyses in other clades, pointing to a need for new analytical models that incorporate the possibility of such events.

Complexity of avian evolution revealed by family-level genomes.

Despite tremendous efforts in the past decades, relationships among main avian lineages remain heavily debated without a clear resolution. Discrepancies have been attributed to diversity of species sampled, phylogenetic method and the choice of genomic regions1-3. Here we address these issues by analysing the genomes of 363 bird species4 (218 taxonomic families, 92% of total). Using intergenic regions and coalescent methods, we present a well-supported tree but also a marked degree of discordance. The tree confirms that Neoaves experienced rapid radiation at or near the Cretaceous-Palaeogene boundary. Sufficient loci rather than extensive taxon sampling were more effective in resolving difficult nodes. Remaining recalcitrant nodes involve species that are a challenge to model due to either extreme DNA composition, variable substitution rates, incomplete lineage sorting or complex evolutionary events such as ancient hybridization. Assessment of the effects of different genomic partitions showed high heterogeneity across the genome. We discovered sharp increases in effective population size, substitution rates and relative brain size following the Cretaceous-Palaeogene extinction event, supporting the hypothesis that emerging ecological opportunities catalysed the diversification of modern birds. The resulting phylogenetic estimate offers fresh insights into the rapid radiation of modern birds and provides a taxon-rich backbone tree for future comparative studies.

ClockstaRX: Testing Molecular Clock Hypotheses With Genomic Data.

Phylogenomic data provide valuable opportunities for studying evolutionary rates and timescales. These analyses require theoretical and statistical tools based on molecular clocks. We present ClockstaRX, a flexible platform for exploring and testing evolutionary rate signals in phylogenomic data. Here, information about evolutionary rates in branches across gene trees is placed in Euclidean space, allowing data transformation, visualization, and hypothesis testing. ClockstaRX implements formal tests for identifying groups of loci and branches that make a large contribution to patterns of rate variation. This information can then be used to test for drivers of genomic evolutionary rates or to inform models for molecular dating. Drawing on the results of a simulation study, we recommend forms of data exploration and filtering that might be useful prior to molecular-clock analyses.

Confusion will be my epitaph: genome-scale discordance stifles phylogenetic resolution of Holothuroidea.

Sea cucumbers (Holothuroidea) are a diverse clade of echinoderms found from intertidal waters to the bottom of the deepest oceanic trenches. Their reduced skeletons and limited number of phylogenetically informative traits have long obfuscated morphological classifications. Sanger-sequenced molecular datasets have also failed to constrain the position of major lineages. Noteworthy, topological uncertainty has hindered a resolution for Neoholothuriida, a highly diverse clade of Permo-Triassic age. We perform the first phylogenomic analysis of Holothuroidea, combining existing datasets with 13 novel transcriptomes. Using a highly curated dataset of 1100 orthologues, our efforts recapitulate previous results, struggling to resolve interrelationships among neoholothuriid clades. Three approaches to phylogenetic reconstruction (concatenation under both site-homogeneous and site-heterogeneous models, and coalescent-aware inference) result in alternative resolutions, all of which are recovered with strong support and across a range of datasets filtered for phylogenetic usefulness. We explore this intriguing result using gene-wise log-likelihood scores and attempt to correlate these with a large set of gene properties. While presenting novel ways of exploring and visualizing support for alternative trees, we are unable to discover significant predictors of topological preference, and our efforts fail to favour one topology. Neoholothuriid genomes seem to retain an amalgam of signals derived from multiple phylogenetic histories.

Range-wide population genomics of common seadragons shows secondary contact over a former barrier and insights on illegal capture.

BACKGROUND: Common seadragons (Phyllopteryx taeniolatus, Syngnathidae) are an emblem of the diverse endemic fauna of Australia's southern rocky reefs, the newly recognized "Great Southern Reef." A lack of assessments spanning this global biodiversity hotspot in its entirety is currently hampering an understanding of the factors that have contributed to its diversity. The common seadragon has a wide range across Australia's entire temperate south and includes a geogenetic break over a former land bridge, which has called its status as a single species into question. As a popular aquarium display that sells for high prices, common seadragons are also vulnerable to illegal capture. RESULTS: Here, we provide range-wide nuclear sequences (986 variable Ultraconserved Elements) for 198 individuals and mitochondrial genomes for 140 individuals to assess species status, identify genetic units and their diversity, and trace the source of two poached individuals. Using published data of the other two seadragon species, we found that lineages of common seadragons have diverged relatively recently (< 0.63 Ma). Within common seadragons, we found pronounced genetic structure, falling into three major groups in the western, central, and eastern parts of the range. While populations across the Bassian Isthmus were divergent, there is also evidence for secondary contact since the passage opened. We found a strong cline of genetic diversity from the range center tapering symmetrically towards the range peripheries. Based on their genetic similarities, the poached individuals were inferred to have originated from around Albany in southwestern Australia. CONCLUSIONS: We conclude that common seadragons constitute a single species with strong geographic structure but coherence through gene flow. The low genetic diversity on the east and west coasts is concerning given that these areas are projected to face fast climate change. Our results suggest that in addition to their life history, geological events and demographic expansions have all played a role in shaping populations in the temperate south. These insights are an important step towards understanding the historical determinants of the diversity of species endemic to the Great Southern Reef.

Evolution of the germline mutation rate across vertebrates.

The germline mutation rate determines the pace of genome evolution and is an evolving parameter itself1. However, little is known about what determines its evolution, as most studies of mutation rates have focused on single species with different methodologies2. Here we quantify germline mutation rates across vertebrates by sequencing and comparing the high-coverage genomes of 151 parent-offspring trios from 68 species of mammals, fishes, birds and reptiles. We show that the per-generation mutation rate varies among species by a factor of 40, with mutation rates being higher for males than for females in mammals and birds, but not in reptiles and fishes. The generation time, age at maturity and species-level fecundity are the key life-history traits affecting this variation among species. Furthermore, species with higher long-term effective population sizes tend to have lower mutation rates per generation, providing support for the drift barrier hypothesis3. The exceptionally high yearly mutation rates of domesticated animals, which have been continually selected on fecundity traits including shorter generation times, further support the importance of generation time in the evolution of mutation rates. Overall, our comparative analysis of pedigree-based mutation rates provides ecological insights on the mutation rate evolution in vertebrates.

Molecular exploration of fossil eggshell uncovers hidden lineage of giant extinct bird.

The systematics of Madagascar's extinct elephant birds remains controversial due to large gaps in the fossil record and poor biomolecular preservation of skeletal specimens. Here, a molecular analysis of 1000-year-old fossil eggshells provides the first description of elephant bird phylogeography and offers insight into the ecology and evolution of these flightless giants. Mitochondrial genomes from across Madagascar reveal genetic variation that is correlated with eggshell morphology, stable isotope composition, and geographic distribution. The elephant bird crown is dated to ca. 30 Mya, when Madagascar is estimated to have become less arid as it moved northward. High levels of between-clade genetic variation support reclassifying Mullerornis into a separate family. Low levels of within-clade genetic variation suggest there were only two elephant bird genera existing in southern Madagascar during the Holocene. However, we find an eggshell collection from Madagascar's far north that represents a unique lineage of Aepyornis. Furthermore, divergence within Aepyornis coincides with the aridification of Madagascar during the early Pleistocene ca. 1.5 Ma, and is consistent with the fragmentation of populations in the highlands driving diversification and the evolution of extreme gigantism over shorts timescales. We advocate for a revision of their taxonomy that integrates palaeogenomic and palaeoecological perspectives.

Epistatic Effects Between Amino Acid Insertions and Substitutions Mediate Toxin resistance of Vertebrate Na+,K+-ATPases.

The recurrent evolution of resistance to cardiotonic steroids (CTS) across diverse animals most frequently involves convergent amino acid substitutions in the H1-H2 extracellular loop of Na+,K+-ATPase (NKA). Previous work revealed that hystricognath rodents (e.g., chinchilla) and pterocliform birds (sandgrouse) have convergently evolved amino acid insertions in the H1-H2 loop, but their functional significance was not known. Using protein engineering, we show that these insertions have distinct effects on CTS resistance in homologs of each of the two species that strongly depend on intramolecular interactions with other residues. Removing the insertion in the chinchilla NKA unexpectedly increases CTS resistance and decreases NKA activity. In the sandgrouse NKA, the amino acid insertion and substitution Q111R both contribute to an augmented CTS resistance without compromising ATPase activity levels. Molecular docking simulations provide additional insight into the biophysical mechanisms responsible for the context-specific mutational effects on CTS insensitivity of the enzyme. Our results highlight the diversity of genetic substrates that underlie CTS insensitivity in vertebrate NKA and reveal how amino acid insertions can alter the phenotypic effects of point mutations at key sites in the same protein domain.

Ancient proteins resolve controversy over the identity of Genyornis eggshell.

The realization that ancient biomolecules are preserved in "fossil" samples has revolutionized archaeological science. Protein sequences survive longer than DNA, but their phylogenetic resolution is inferior; therefore, careful assessment of the research questions is required. Here, we show the potential of ancient proteins preserved in Pleistocene eggshell in addressing a longstanding controversy in human and animal evolution: the identity of the extinct bird that laid large eggs which were exploited by Australia's indigenous people. The eggs had been originally attributed to the iconic extinct flightless bird Genyornis newtoni (†Dromornithidae, Galloanseres) and were subsequently dated to before 50 ± 5 ka by Miller et al. [Nat. Commun. 7, 10496 (2016)]. This was taken to represent the likely extinction date for this endemic megafaunal species and thus implied a role of humans in its demise. A contrasting hypothesis, according to which the eggs were laid by a large mound-builder megapode (Megapodiidae, Galliformes), would therefore acquit humans of their responsibility in the extinction of Genyornis. Ancient protein sequences were reconstructed and used to assess the evolutionary proximity of the undetermined eggshell to extant birds, rejecting the megapode hypothesis. Authentic ancient DNA could not be confirmed from these highly degraded samples, but morphometric data also support the attribution of the eggshell to Genyornis. When used in triangulation to address well-defined hypotheses, paleoproteomics is a powerful tool for reconstructing the evolutionary history in ancient samples. In addition to the clarification of phylogenetic placement, these data provide a more nuanced understanding of the modes of interactions between humans and their environment.

Phylogenomic analysis of Syngnathidae reveals novel relationships, origins of endemic diversity and variable diversification rates.

BACKGROUND: Seahorses, seadragons, pygmy pipehorses, and pipefishes (Syngnathidae, Syngnathiformes) are among the most recognizable groups of fishes because of their derived morphology, unusual life history, and worldwide distribution. Despite previous phylogenetic studies and recent new species descriptions of syngnathids, the evolutionary relationships among several major groups within this family remain unresolved. RESULTS: Here, we provide a reconstruction of syngnathid phylogeny based on genome-wide sampling of 1314 ultraconserved elements (UCEs) and expanded taxon sampling to assess the current taxonomy and as a basis for macroevolutionary insights. We sequenced a total of 244 new specimens across 117 species and combined with published UCE data for a total of 183 species of Syngnathidae, about 62% of the described species diversity, to compile the most data-rich phylogeny to date. We estimated divergence times using 14 syngnathiform fossils, including nine fossils with newly proposed phylogenetic affinities, to better characterize current and historical biogeographical patterns, and to reconstruct diversification through time. We present a phylogenetic hypothesis that is well-supported and provides several notable insights into syngnathid evolution. We found nine non-monophyletic genera, evidence for seven cryptic species, five potentially invalid synonyms, and identified a novel sister group to the seahorses, the Indo-Pacific pipefishes Halicampus macrorhynchus and H. punctatus. In addition, the morphologically distinct southwest Pacific seahorse Hippocampus jugumus was recovered as the sister to all other non-pygmy seahorses. As found in many other groups, a high proportion of syngnathid lineages appear to have originated in the Central Indo-Pacific and subsequently dispersed to adjoining regions. Conversely, we also found an unusually high subsequent return of lineages from southern Australasia to the Central Indo-Pacific. Diversification rates rose abruptly during the Middle Miocene Climate Transition and peaked after the closure of the Tethys Sea. CONCLUSIONS: Our results reveal a previously underappreciated diversity of syngnathid lineages. The observed biogeographic patterns suggest a significant role of the southern Australasian region as a source and sink of lineages. Shifts in diversification rates imply possible links to declining global temperatures, the separation of the Atlantic and Pacific faunas, and the environmental changes associated with these events.

Phylogenomics resolves ambiguous relationships within Aciculata (Errantia, Annelida).

Aciculata (Eunicida + Phyllodocida) is among the largest clades of annelids, comprising almost half of the known diversity of all marine annelids. Despite the group's large size and biological importance, most phylogenomic studies on Annelida to date have had a limited sampling of this clade. The phylogenetic placement of many clades within Phyllodocida in particular has remained poorly understood. To resolve the relationships within Aciculata we conducted a large-scale phylogenomic analysis based on 24 transcriptomes (13 new), chosen to represent many family-ranked taxa that have never been included in a broad phylogenomic study. Our sampling also includes several enigmatic taxa with challenging phylogenetic placement, such as Histriobdella, Struwela, Lacydonia, Pilargis and the holopelagic worms Lopadorrhynchus, Travisiopsis and Tomopteris. Our robust phylogeny allows us to name and place some of these problematic clades and has significant implications on the systematics of the group. Within Eunicida we reinstate the names Eunicoidea and Oenonoidea. Within Phyllodocida we delineate Phyllodociformia, Glyceriformia, Nereidiformia, Nephtyiformia and Aphroditiformia. Phyllodociformia now includes: Lacydonia, Typhloscolecidae, Lopadorrhynchidae and Phyllodocidae. Nephtyiformia includes Nephtyidae and Pilargidae. We also broaden the delineation of Glyceriformia to include Sphaerodoridae, Tomopteridae and Glyceroidea (Glyceridae + Goniadidae). Furthermore, our study demonstrates and explores how conflicting, yet highly supported topologies can result from confounding signals in gene trees.

Phylogenomic analyses of mud dragons (Kinorhyncha).

Mud dragons (Kinorhyncha) are microscopic invertebrates, inhabiting marine sediments across the globe from intertidal to hadal depths. They are segmented, moulting animals like arthropods, but grouping with the unsegmented priapulans and loriciferans within Ecdysozoa. There are more than 300 species of kinorhynchs described within 31 genera and 11 families, however, their evolutionary relationships have so far only been investigated using morphology and a few molecular markers. Here we aim to resolve the relationships and classification of major clades within Kinorhyncha using transcriptomic data. In addition, we wish to revisit the position of three indistinctly segmented, aberrant genera in order to reconstruct the evolution of distinct segmentation within the group. We conducted a phylogenomic analysis of Kinorhyncha including 21 kinorhynch transcriptomes (of which 18 are new) representing 15 genera, and seven outgroups including priapulan, loriciferan, nematode and nematomorph transcriptomes. Results show a congruent and robust tree that supports the division of Kinorhyncha into two major clades: Cyclorhagida and Allomalorhagida. Cyclorhagida is composed of three subclades: Xenosomata, Kentrorhagata comb. nov. (including the aberrant Zelinkaderes) and Echinorhagata. Allomalorhagida is composed of two subclades: Pycnophyidae and Anomoirhaga nom. nov. Anomoirhaga nom. nov. accommodates the aberrant genera Cateria (previously nested within Cyclorhagida) and Franciscideres together with five additional genera. The distant and derived positions of the aberrant Zelinkaderes, Cateria and Franciscideres species suggest that their less distinct trunk segmentation evolved convergently, and that segmentation evolved among kinorhynch stem groups.

Multiple origins of a frameshift insertion in a mitochondrial gene in birds and turtles.

BACKGROUND: During evolutionary history, molecular mechanisms have emerged to cope with deleterious mutations. Frameshift insertions in protein-coding sequences are extremely rare because they disrupt the reading frame. There are a few known examples of their correction through translational frameshifting, a process that enables ribosomes to skip nucleotides during translation to regain proper reading frame. Corrective frameshifting has been proposed to act on the single base pair insertion at position 174 of the mitochondrial NADH dehydrogenase subunit 3 gene (ND3) that has been observed in several turtles and birds. However, the relatively sparse taxonomic representation has hampered our understanding of the evolution of this insertion in vertebrates. RESULTS: Here, we analyzed 87,707 ND3 sequences from 10,309 vertebrate taxa to reveal the evolutionary history of this insertion and its common genomic characteristics. We confirmed that the insertion only appears in turtles and birds and reconstructed that it evolved independently in both groups with complex patterns of gains and losses. The insertion was observed in almost all bird orders but was absent in all members of the diverse Passeriformes. We found strong conservation in the nucleotides surrounding the insertion in both turtles and birds, which implies that the insertion enforces structural constraints that could be involved in its correction. CONCLUSIONS: Our study demonstrates that frameshifts can be widespread and can be retained for millions of years if they are embedded in a conserved sequence theme.

Using ultraconserved elements to track the influence of sea-level change on leafy seadragon populations.

During the Last Glacial Maximum (LGM), global sea levels were 120-130 m lower than today, resulting in the emergence of most continental shelves and extirpation of subtidal organisms from these areas. During the interglacial periods, rapid inundation of shelf regions created a dynamic environment for coastal organisms, such as the charismatic leafy seadragon (Phycodurus eques, Syngnathidae), a brooder with low dispersal ability inhabiting kelp beds in temperate Australia. Reconstructions of the palaeoshoreline revealed that the increase of shallow areas since the LGM was not uniform across the species' range and we investigated the effects of these asymmetries on genetic diversity and structuring. Using targeted capture of 857 variable ultraconserved elements (UCEs, 2,845 single nucleotide polymorphisms) in 68 individuals, we found that the regionally different shelf topographies were paralleled by contrasting population genetic patterns. In the west, populations may not have persisted through sea-level lows because shallow seabed was very limited. Shallow genetic structure, weak expansion signals and a westward cline in genetic diversity indicate a postglacial recolonization of the western part of the range from a more eastern location following sea-level rise. In the east, shallow seabed persisted during the LGM and increased considerably after the flooding of large bays, which resulted in strong demographic expansions, deeper genetic structure and higher genetic diversity. This study suggests that postglacial flooding with rising sea levels produced locally variable signatures in colonizing populations.

Progressive Cactus is a multiple-genome aligner for the thousand-genome era.

New genome assemblies have been arriving at a rapidly increasing pace, thanks to decreases in sequencing costs and improvements in third-generation sequencing technologies1-3. For example, the number of vertebrate genome assemblies currently in the NCBI (National Center for Biotechnology Information) database4 increased by more than 50% to 1,485 assemblies in the year from July 2018 to July 2019. In addition to this influx of assemblies from different species, new human de novo assemblies5 are being produced, which enable the analysis of not only small polymorphisms, but also complex, large-scale structural differences between human individuals and haplotypes. This coming era and its unprecedented amount of data offer the opportunity to uncover many insights into genome evolution but also present challenges in how to adapt current analysis methods to meet the increased scale. Cactus6, a reference-free multiple genome alignment program, has been shown to be highly accurate, but the existing implementation scales poorly with increasing numbers of genomes, and struggles in regions of highly duplicated sequences. Here we describe progressive extensions to Cactus to create Progressive Cactus, which enables the reference-free alignment of tens to thousands of large vertebrate genomes while maintaining high alignment quality. We describe results from an alignment of more than 600 amniote genomes, which is to our knowledge the largest multiple vertebrate genome alignment created so far.

Dense sampling of bird diversity increases power of comparative genomics.

Whole-genome sequencing projects are increasingly populating the tree of life and characterizing biodiversity1-4. Sparse taxon sampling has previously been proposed to confound phylogenetic inference5, and captures only a fraction of the genomic diversity. Here we report a substantial step towards the dense representation of avian phylogenetic and molecular diversity, by analysing 363 genomes from 92.4% of bird families-including 267 newly sequenced genomes produced for phase II of the Bird 10,000 Genomes (B10K) Project. We use this comparative genome dataset in combination with a pipeline that leverages a reference-free whole-genome alignment to identify orthologous regions in greater numbers than has previously been possible and to recognize genomic novelties in particular bird lineages. The densely sampled alignment provides a single-base-pair map of selection, has more than doubled the fraction of bases that are confidently predicted to be under conservation and reveals extensive patterns of weak selection in predominantly non-coding DNA. Our results demonstrate that increasing the diversity of genomes used in comparative studies can reveal more shared and lineage-specific variation, and improve the investigation of genomic characteristics. We anticipate that this genomic resource will offer new perspectives on evolutionary processes in cross-species comparative analyses and assist in efforts to conserve species.

An Indo-Pacific Humpback Dolphin Genome Reveals Insights into Chromosome Evolution and the Demography of a Vulnerable Species.

The Indo-Pacific humpback dolphin (Sousa chinensis) is a small inshore species of odontocete cetacean listed as Vulnerable on the IUCN Red List. Here, we report on the evolution of S. chinensis chromosomes from its cetruminant ancestor and elucidate the evolutionary history and population genetics of two neighboring S. chinensis populations. We found that breakpoints in ancestral chromosomes leading to S. chinensis could have affected the function of genes related to kidney filtration, body development, and immunity. Resequencing of individuals from two neighboring populations in the northwestern South China Sea, Leizhou Bay and Sanniang Bay, revealed genetic differentiation, low diversity, and small contemporary effective population sizes. Demographic analyses showed a marked decrease in the population size of the two investigated populations over the last ~4,000 years, possibly related to climatic oscillations. This study implies a high risk of extinction and strong conservation requirement for the Indo-Pacific humpback dolphin.

More is needed-Thousands of loci are required to elucidate the relationships of the 'flowers of the sea' (Sabellida, Annelida).

Sabellida is a well-known clade containing tube-dwelling annelid worms with a radiolar crown. Iterative phylogenetic analyses over three decades have resulted in three main clades being recognized; Fabriciidae, Serpulidae and Sabellidae, with Fabriciidae proposed as the sister group to Serpulidae. However, relationships within Sabellidae have remained poorly understood, with a proliferation of genera. In order to obtain a robust phylogeny with optimal support, we conducted a large-scale phylogenomic analysis with 19 new sabellid transcriptomes for a total of 21 species. In contrast to earlier findings based on limited DNA data, our results support the position of Fabriciidae as sister taxon to a Sabellidae + Serpulidae clade. Our large sampling within Sabellidae also allows us to establish a stable phylogeny within this clade. We restrict Sabellinae to a subclade of Sabellidae and broaden the previously monotypic Myxicolinae to include Amphicorina and Chone. We tested the robustness of species tree reconstruction by subsampling increasing numbers of genes to uncover hidden support of alternative topologies. Our results show that inclusion of more genes leads to a more stable topology with higher support, and also that including higher divergence genes leads to stronger resolution.

Spaghetti to a Tree: A Robust Phylogeny for Terebelliformia (Annelida) Based on Transcriptomes, Molecular and Morphological Data.

Terebelliformia-"spaghetti worms" and their allies-are speciose and ubiquitous marine annelids but our understanding of how their morphological and ecological diversity evolved is hampered by an uncertain delineation of lineages and their phylogenetic relationships. Here, we analyzed transcriptomes of 20 terebelliforms and an outgroup to build a robust phylogeny of the main lineages grounded on 12,674 orthologous genes. We then supplemented this backbone phylogeny with a denser sampling of 121 species using five genes and 90 morphological characters to elucidate fine-scale relationships. The monophyly of six major taxa was supported: Pectinariidae, Ampharetinae, Alvinellidae, Trichobranchidae, Terebellidae and Melinninae. The latter, traditionally a subfamily of Ampharetidae, was unexpectedly the sister to Terebellidae, and hence becomes Melinnidae, and Ampharetinae becomes Ampharetidae. We found no support for the recently proposed separation of Telothelepodidae, Polycirridae and Thelepodidae from Terebellidae. Telothelepodidae was nested within Thelepodinae and is accordingly made its junior synonym. Terebellidae contained the subfamily-ranked taxa Terebellinae and Thelepodinae. The placement of the simplified Polycirridae within Terebellinae differed from previous hypotheses, warranting the division of Terebellinae into Lanicini, Procleini, Terebellini and Polycirrini. Ampharetidae (excluding Melinnidae) were well-supported as the sister group to Alvinellidae and we recognize three clades: Ampharetinae, Amaginae and Amphicteinae. Our analysis found several paraphyletic genera and undescribed species. Morphological transformations on the phylogeny supported the hypothesis of an ancestor that possessed both branchiae and chaetae, which is at odds with proposals of a "naked" ancestor. Our study demonstrates how a robust backbone phylogeny can be combined with dense taxon coverage and morphological traits to give insights into the evolutionary history and transformation of traits.