Monophyly of brachiopods and phoronids: reconciliation of molecular evidence with Linnaean classification (the subphylum Phoroniformea nov.)

Molecular phylogenetic analyses of aligned 18S rDNA gene sequences from articulate and inarticulate brachiopods representing all major extant lineages, an enhanced set of phoronids and several unrelated protostome taxa, confirm previous indications that in such data, brachiopod and phoronids form a well-supported clade that (on previous evidence) is unambiguously affiliated with protostomes rather than deuterostomes. Within the brachiopod-phoronid clade, an association between phoronids and inarticulate brachiopods is moderately well supported, whilst a close relationship between phoronids and craniid inarticulates is weakly indicated. Brachiopod-phoronid monophyly is reconciled with the most recent Linnaean classification of brachiopods by abolition of the phylum Phoronida and rediagnosis of the phylum Brachiopoda to include tubiculous, shell-less forms. Recognition that brachiopods and phoronids are close genealogical allies of protostome phyla such as molluscs and annelids, but are much more distantly related to deuterostome phyla such as echinoderms and chordates, implies either (or both) that the morphology and ontogeny of blastopore, mesoderm and coelom formation have been widely misreported or misinterpreted, or that these characters have been subject to extensive homoplasy. This inference, if true, undermines virtually all morphology-based reconstructions of phylogeny made during the past century or more

The eoorthid brachiopod Apheoorthina in the Lower Ordovician of NW Argentina and the dispersal pathways along western Gondwana

The eoorthid brachiopod Apheoorthina is reported for the first time from the Lower Ordovician of NW Argentina. It is represented by a species similar to A. ferrigena from the Tremadocian of the Prague Basin, increasing the faunal affinities between the Central Andean Basin and the South European microcontinents, in particular the Bohemian region (Perunica). Nine out of the fourteen brachiopod genera reported from the Tremadocian of the Central Andean Basin (~64 %) are shared with the Mediterranean region, four of which (~ 28 %) have been recorded in the Prague Basin, and two (Kvania and Apheoorthina) are restricted to the Central Andes and Perunica. Dispersal pathways around Gondwana are analyzed in the light of major factors affecting large-scale distribution of brachiopods (environment, larval capacity for dispersal, oceanic currents).The presence in Apheoorthina aff. ferrigena of a well-preserved larval protegulum measuring 420 μm in width and 210 μm in length strongly suggests that this species had planktotrophic larvae capable of long-distance dispersal. According to recent oceanatmosphere general circulation models for the Ordovician Period, the Central Andean margin was dominated by the cold-water Antarctica Current. Despite the complex nonzonal pattern produced by current deflections around the peri-Gondwanan microcontinents, the general westward circulation sense favoured larval dispersal from the Andean region to North Africa, Avalonia, the Armorican Terrane Assemblage, and Perunica. On the other hand, the eastwards flowing Gondwana Current connected the North Gondwana waters with the South American epicontinental seas, which could explain the reversed migration of some brachiopods.Fil: Muñoz, Diego Fernando. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Córdoba. Centro de Investigaciones en Ciencias de la Tierra. Universidad Nacional de Córdoba. Facultad de Ciencias Exactas Físicas y Naturales. Centro de Investigaciones en Ciencias de la Tierra; ArgentinaFil: Benedetto, Juan Luis Arnaldo. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Córdoba. Centro de Investigaciones en Ciencias de la Tierra. Universidad Nacional de Córdoba. Facultad de Ciencias Exactas Físicas y Naturales. Centro de Investigaciones en Ciencias de la Tierra; Argentin

DNA sequence evidence for speciation, paraphyly and a Mesozoic dispersal of cancellothyridid articulate brachiopods

Because the classification of extant and fossil articulate brachiopods is based largely upon shell characters observable in fossils, it identifies morphotaxa whose biological status can, in practice, best be inferred from estimates of genetic divergence. Allozyme polymorphism and restriction fragment length polymorphism of mitochondrial DNA (mtDNA RFLP) have been used to show that nuclear and mitochondrial genetic divergence between samples of the cancellothyridid brachiopods Terebratulina septentrionalis from Canada and T. retusa from Europe is compatible with biological speciation, but the genetic distances obtained were biased by methodological limitations. Here, we report estimates of divergence in 12S rDNA mitochondrial sequences within and between samples of these brachiopods. The sequence-based genetic distance between these samples (5.98-0.07% SE) is at least 10 times greater than within them and, since they also differ in a complex life-history trait, their species status is considered to be securely established. Divergence levels between 12S rDNA genes of three other cancellothyridids, T. unguicula from Alaska, T. crossei from near Japan, and Cancellothyris hedleyi from near Australia are higher than between the two North Atlantic species, and the mean nucleotide distance between all these cancellothyrids is similar to the mean distance between species of Littorina (Mollusca: Gastropoda). Sequences of both 12S and 16S genes from cancellothyridids and other short-looped brachiopod species show neither saturation nor lineage-specific rate differences and, when analysed with different outgroups, either separately or together, yield one unexpected, but well-supported, tree with Alaskan T. unguicula basal and C. hedleyi nested within Terebratulina, i.e. these genera are paraphyletic. A geologically dated divergence between Antarctic and New Zealand species of the short-looped brachiopod Liothyrella is used to calibrate the rate of 12S divergence at ca. 0.1% per million years (MY), and this rate is used to infer that T. septentrionalis and T. retusa have been diverging for ca. 60 MY and that they and T. unguicula have been diverging from their last common ancestor for ca. 100 MY. This indicates a Mesozoic origin for the present-day distribution of cancellothyridids and the basal position of T. unguicula suggests a possible North Pacific centre of origin, with separate Atlantic and Pacific radiations. The inclusion of Cancellothyris within Terebratulina also shows that adult shell characters such as umbo, foramen and symphytium shape, whilst probably indispensible for the practical classification of fossils, are not reliable guides to genealogy

Comparison of articulate brachiopod nuclear and mitochondrial gene trees leads to a clade-based redefinition of protostomes (Protostomozoa) and deuterostomes (Deuterostomozoa)

Nuclear and mtDNA sequences from selected short-looped terebratuloid (terebratulacean) articulate brachiopods yield congruent and genetically independent phylogenetic reconstructions by parsimony, neighbor-joining and maximum likelihood methods, suggesting that both sources of data are reliable guides to brachiopod species phylogeny. The present-day genealogical relationships and geographical distributions of the tested terebratuloid brachiopods are consistent with a tethyan dispersal and subsequent radiation. Concordance of nuclear and mitochondrial gene phylogenies reinforces previous indications that articulate brachiopods, inarticulate brachiopods, phoronids and ectoprocts cluster with other organisms generally regarded as protostomes. Since ontogeny and morphology in brachiopods, ectoprocts and phoronids depart in important respects from those features supposedly diagnostic of protostomes, this demonstrates that the operational definition of protostomy by the usual ontological characters must be misleading or unreliable. New, molecular, operational definitions are proposed to replace the traditional criteria for the recognition of protostomes and deuterostomes, and the clade-based terms 'Protostomozoa' and 'Deuterostomozoa' are proposed to replace the existing terms 'Protostomia' and 'Deuterostomia'

Molecular evidence that phoronids are a subtaxon of brachiopods (Brachiopoda: Phoronata) and that genetic divergence of metazoan phyla began long before the early Cambrian

Concatenated SSU (18S) and partial LSU (28S) sequences (~2 kb) from 12 ingroup taxa, comprising 2 phoronids, 2 members of each of the craniid, discinid, and lingulid inarticulate brachiopod lineages, and 4 rhynchonellate, articulate brachiopods (2 rhynchonellides, 1 terebratulide and 1 terebratellide) were aligned with homologous sequences from 6 protostome, deuterostome and sponge outgroups (3964 sites). Regions of potentially ambiguous alignment were removed, and the resulting data (3275 sites, of which 377 were parsimony-informative and 635 variable) were analysed by parsimony, and by maximum and Bayesian likelihood using objectively selected models. There was no base composition heterogeneity. Relative rate tests led to the exclusion (from most analyses) of the more distant outgroups, with retention of the closer pectinid and polyplacophoran (chiton). Parsimony and likelihood bootstrap and Bayesian clade support values were generally high, but only likelihood analyses recovered all brachiopod indicator clades designated a priori. All analyses confirmed the monophyly of (brachiopods+phoronids) and identified phoronids as the sister-group of the three inarticulate brachiopod lineages. Consequently, a revised Linnean classification is proposed in which the subphylum Linguliformea comprises three classes: Lingulata, ‘Phoronata’ (the phoronids), and ‘Craniata’ (the current subphylum Craniiformea). Divergence times of all nodes were estimated by regression from node depths in non-parametrically rate-smoothed and other chronograms, calibrated against palaeontological data, with probable errors not less than 50 My. Only three predicted brachiopod divergence times disagree with palaeontological ages by more than the probable error, and a reasonable explanation exists for at least two. Pruning long-branched ingroups made scant difference to predicted divergence time estimates. The palaeontological age calibration and the existence of Lower Cambrian fossils of both main brachiopod clades together indicate that initial genetic divergence between brachiopod and molluscan (chiton) lineages occurred well before the Lower Cambrian, suggesting that much divergence between metazoan phyla took place in the Proterozoic

Molecular phylogeny of brachiopods and phoronids based on nuclear-encoded small subunit ribosomal RNA gene sequences

Brachiopod and phoronid phylogeny is inferred from SSU rDNA sequences of 28 articulate and nine inarticulate brachiopods, three phoronids, two ectoprocts and various outgroups, using gene trees reconstructed by weighted parsimony, distance and maximum likelihood methods. Of these sequences, 33 from brachiopods, two from phoronids and one each from an ectoproct and a priapulan are newly determined. The brachiopod sequences belong to 31 different genera and thus survey about 10% of extant genus-level diversity. Sequences determined in different laboratories and those from closely related taxa agree well, but evidence is presented suggesting that one published phoronid sequence (GenBank accession UO12648) is a brachiopod-phoronid chimaera, and this sequence is excluded from the analyses. The chiton, Acanthopleura, is identified as the phenetically proximal outgroup; other selected outgroups were chosen to allow comparison with recent, non-molecular analyses of brachiopod phylogeny. The different outgroups and methods of phylogenetic reconstruction lead to similar results, with differences mainly in the resolution of weakly supported ancient and recent nodes, including the divergence of inarticulate brachiopod sub-phyla, the position of the rhynchonellids in relation to long- and short-looped articulate brachiopod clades and the relationships of some articulate brachiopod genera and species. Attention is drawn to the problem presented by nodes that are strongly supported by non-molecular evidence but receive only low bootstrap resampling support. Overall, the gene trees agree with morphology-based brachiopod taxonomy, but novel relationships are tentatively suggested for thecideidine and megathyrid brachiopods. Articulate brachiopods are found to be monophyletic in all reconstructions, but monophyly of inarticulate brachiopods and the possible inclusion of phoronids in the inarticulate brachiopod clade are less strongly established. Phoronids are clearly excluded from a sister-group relationship with articulate brachiopods, this proposed relationship being due to the rejected, chimaeric sequence (GenBank UO12648). Lineage relative rate tests show no heterogeneity of evolutionary rate among articulate brachiopod sequences, but indicate that inarticulate brachiopod plus phoronid sequences evolve somewhat more slowly. Both brachiopods and phoronids evolve slowly by comparison with other invertebrates. A number of palaeontologically dated times of earliest appearance are used to make upper and lower estimates of the global rate of brachiopod SSU rDNA evolution, and these estimates are used to infer the likely divergence times of other nodes in the gene tree. There is reasonable agreement between most inferred molecular and palaeontological ages. The estimated rates of SSU rDNA sequence evolution suggest that the last common ancestor of brachiopods, chitons and other protostome invertebrates (Lophotrochozoa and Ecdysozoa) lived deep in Precambrian time. Results of this first DNA-based, taxonomically representative analysis of brachiopod phylogeny are in broad agreement with current morphology-based classification and systematics and are largely consistent with the hypothesis that brachiopod shell ontogeny and morphology are a good guide to phylogeny

The Mg/Ca–temperature relationship in brachiopod shells: calibrating a potential palaeoseasonality proxy

Brachiopods are long-lived, long-ranging, extant organisms, of which some groups precipitate a relatively diagenetically stable low magnesium calcite shell. Previous work has suggested that the incorporation of Mg into brachiopod calcite may be controlled by temperature (Brand et al., 2013). Here we build upon this work by using laser ablation sampling to define the intra-shell variations in two modern brachiopod species,Terebratulina retusa (Linnaeus, 1758) and Liothyrella neozelanica (Thomson, 1918). We studied three T. retusa shells collected live from the Firth of Lorne, Scotland, which witnessed annual temperature variations on the order of 7 °C, in addition to four L. neozelanica shells, which were dredged from a water depth transect (168–1488 m) off the north coast of New Zealand. The comparison of intra-shell Mg/Ca profiles with shell δ<sup>18</sup>O confirms a temperature control on brachiopod Mg/Ca and supports the use of brachiopod Mg/Ca as a palaeoseasonality indicator. Our preliminary temperature calibrations are Mg/Ca = 1.76 ± 0.27 e<sup>(0.16 ± 0.03)T</sup>, R<sup>2</sup> = 0.75, for T. retusa and Mg/Ca = 0.49 ± 1.27 e<sup>(0.2 ± 0.11)T</sup>, R<sup>2</sup> = 0.32, for L. neozelanica (errors are 95% confidence intervals)