Phylogenomic analyses of echinoid diversification prompt a re-evaluation of their fossil record.

Echinoids are key components of modern marine ecosystems. Despite a remarkable fossil record, the emergence of their crown group is documented by few specimens of unclear affinities, rendering their early history uncertain. The origin of sand dollars, one of its most distinctive clades, is also unclear due to an unstable phylogenetic context. We employ 18 novel genomes and transcriptomes to build a phylogenomic dataset with a near-complete sampling of major lineages. With it, we revise the phylogeny and divergence times of echinoids, and place their history within the broader context of echinoderm evolution. We also introduce the concept of a chronospace - a multidimensional representation of node ages - and use it to explore methodological decisions involved in time calibrating phylogenies. We find the choice of clock model to have the strongest impact on divergence times, while the use of site-heterogeneous models and alternative node prior distributions show minimal effects. The choice of loci has an intermediate impact, affecting mostly deep Paleozoic nodes, for which clock-like genes recover dates more congruent with fossil evidence. Our results reveal that crown group echinoids originated in the Permian and diversified rapidly in the Triassic, despite the relative lack of fossil evidence for this early diversification. We also clarify the relationships between sand dollars and their close relatives and confidently date their origins to the Cretaceous, implying ghost ranges spanning approximately 50 million years, a remarkable discrepancy with their rich fossil record.

Cambrian edrioasteroid reveals new mechanism for secondary reduction of the skeleton in echinoderms.

Echinoderms are characterized by a distinctive high-magnesium calcite endoskeleton as adults, but elements of this have been drastically reduced in some groups. Herein, we describe a new pentaradial echinoderm, Yorkicystis haefneri n. gen. n. sp., which provides, to our knowledge, the oldest evidence of secondary non-mineralization of the echinoderm skeleton. This material was collected from the Cambrian Kinzers Formation in York (Pennsylvania, USA) and is dated as ca 510 Ma. Detailed morphological observations demonstrate that the ambulacra (i.e. axial region) are composed of flooring and cover plates, but the rest of the body (i.e. extraxial region) is preserved as a dark film and lacks any evidence of skeletal plating. Moreover, X-ray fluorescence analysis reveals that the axial region is elevated in iron. Based on our morphological and chemical data and on taphonomic comparisons with other fossils from the Kinzers Formation, we infer that the axial region was originally calcified, while the extraxial region was non-mineralized. Phylogenetic analyses recover Yorkicystis as an edrioasteroid, indicating that this partial absence of skeleton resulted from a secondary reduction. We hypothesize that skeletal reduction resulted from lack of expression of the skeletogenic gene regulatory network in the extraxial body wall during development. Secondary reduction of the skeleton in Yorkicystis might have allowed for greater flexibility of the body wall.

The Development and Neuronal Complexity of Bipinnaria Larvae of the Sea Star Asterias rubens.

Free-swimming planktonic larvae are a key stage in the development of many marine phyla, and studies of these organisms have contributed to our understanding of major genetic and evolutionary processes. Although transitory, these larvae often attain a remarkable degree of tissue complexity, with well-defined musculature and nervous systems. Among the best studied are larvae belonging to the phylum Echinodermata, but with work largely focused on the pluteus larvae of sea urchins (class Echinoidea). The greatest diversity of larval strategies among echinoderms is found in the class Asteroidea (sea stars), organisms that are rapidly emerging as experimental systems for genetic and developmental studies. However, the bipinnaria larvae of sea stars have only been studied in detail in a small number of species and although they have been relatively well described neuro-anatomically, they are poorly understood neurochemically. Here, we have analyzed embryonic development and bipinnaria larval anatomy in the common North Atlantic sea star Asterias rubens, using a variety of staining methods in combination with confocal microscopy. Importantly, the chemical complexity of the nervous system of bipinnaria larvae was revealed through use of a diverse set of antibodies, with identification of at least three centers of differing neurochemical signature within the previously described nervous system: the anterior apical organ, oral region, and ciliary bands. Furthermore, the anatomy of the musculature and sites of cell division in bipinnaria larvae was analyzed. Comparisons of developmental progression and molecular anatomy across the Echinodermata provided a basis for hypotheses on the shared evolutionary and developmental processes that have shaped this group of animals. We conclude that bipinnaria larvae appear to be remarkably conserved across ∼200 million years of evolutionary time and may represent a strong evolutionary and/or developmental constraint on species utilizing this larval strategy.

Post-metamorphic skeletal growth in the sea urchin Paracentrotus lividus and implications for body plan evolution.

BACKGROUND: Understanding the molecular and cellular processes that underpin animal development are crucial for understanding the diversity of body plans found on the planet today. Because of their abundance in the fossil record, and tractability as a model system in the lab, skeletons provide an ideal experimental model to understand the origins of animal diversity. We herein use molecular and cellular markers to understand the growth and development of the juvenile sea urchin (echinoid) skeleton. RESULTS: We developed a detailed staging scheme based off of the first ~ 4 weeks of post-metamorphic life of the regular echinoid Paracentrotus lividus. We paired this scheme with immunohistochemical staining for neuronal, muscular, and skeletal tissues, and fluorescent assays of skeletal growth and cell proliferation to understand the molecular and cellular mechanisms underlying skeletal growth and development of the sea urchin body plan. CONCLUSIONS: Our experiments highlight the role of skeletogenic proteins in accretionary skeletal growth and cell proliferation in the addition of new metameric tissues. Furthermore, this work provides a framework for understanding the developmental evolution of sea urchin body plans on macroevolutionary timescales.

A Total-Evidence Dated Phylogeny of Echinoidea Combining Phylogenomic and Paleontological Data.

Phylogenomic and paleontological data constitute complementary resources for unraveling the phylogenetic relationships and divergence times of lineages, yet few studies have attempted to fully integrate them. Several unique properties of echinoids (sea urchins) make them especially useful for such synthesizing approaches, including a remarkable fossil record that can be incorporated into explicit phylogenetic hypotheses. We revisit the phylogeny of crown group Echinoidea using a total-evidence dating approach that combines the largest phylogenomic data set for the clade, a large-scale morphological matrix with a dense fossil sampling, and a novel compendium of tip and node age constraints. To this end, we develop a novel method for subsampling phylogenomic data sets that selects loci with high phylogenetic signal, low systematic biases, and enhanced clock-like behavior. Our results demonstrate that combining different data sources increases topological accuracy and helps resolve conflicts between molecular and morphological data. Notably, we present a new hypothesis for the origin of sand dollars, and restructure the relationships between stem and crown echinoids in a way that implies a long stretch of undiscovered evolutionary history of the crown group in the late Paleozoic. Our efforts help bridge the gap between phylogenomics and phylogenetic paleontology, providing a model example of the benefits of combining the two. [Echinoidea; fossils; paleontology; phylogenomics; time calibration; total evidence.].

Evolution and Development at the Origin of a Phylum.

Quantifying morphological evolution is key to determining the patterns and processes underlying the origin of phyla. We constructed a hierarchical morphological character matrix to characterize the radiation and establishment of echinoderm body plans during the early Paleozoic. This showed that subphylum-level clades diverged gradually through the Cambrian, and the distinctiveness of the resulting body plans was amplified by the extinction of transitional forms and obscured by convergent evolution during the Ordovician. Higher-order characters that define these body plans were not fixed at the origin of the phylum, countering hypotheses regarding developmental processes governing the early evolution of animals. Instead, these burdened characters were flexible, enabling continued evolutionary innovation throughout the clades' history.

Echinoids from the Tesero Member (Werfen Formation) of the Dolomites (Italy): implications for extinction and survival of echinoids in the aftermath of the end-Permian mass extinction.

The end-Permian mass extinction (∼252 Ma) was responsible for high rates of extinction and evolutionary bottlenecks in a number of animal groups. Echinoids, or sea urchins, were no exception, and the Permian to Triassic represents one of the most significant intervals of time in their macroevolutionary history. The extinction event was responsible for significant turnover, with the Permian-Triassic representing the transition from stem group echinoid-dominated faunas in the Palaeozoic to Mesozoic faunas dominated by crown group echinoids. This turnover is well-known, however, the environmental and taxonomic distribution of echinoids during the latest Permian and Early Triassic is not. Here we report on an echinoid fauna from the Tesero Member, Werfen Formation (latest Permian to Early Triassic) of the Dolomites (northern Italy). The fauna is largely known from disarticulated ossicles, but consists of both stem group taxa, and a new species of crown group echinoid, Eotiaris teseroensis n. sp. That these stem group echinoids were present in the Tesero Member indicates that stem group echinoids did not go extinct in the Dolomites coincident with the onset of extinction, further supporting other recent work indicating that stem group echinoids survived the end-Permian extinction. Furthermore, the presence of Eotiaris across a number of differing palaeoenvironments in the Early Triassic may have had implications for the survival of cidaroid echinoids during the extinction event.

Cell type phylogenetics informs the evolutionary origin of echinoderm larval skeletogenic cell identity.

The multiplicity of cell types comprising multicellular organisms begs the question as to how cell type identities evolve over time. Cell type phylogenetics informs this question by comparing gene expression of homologous cell types in distantly related taxa. We employ this approach to inform the identity of larval skeletogenic cells of echinoderms, a clade for which there are phylogenetically diverse datasets of spatial gene expression patterns. We determined ancestral spatial expression patterns of alx1, ets1, tbr, erg, and vegfr, key components of the skeletogenic gene regulatory network driving identity of the larval skeletogenic cell. Here we show ancestral state reconstructions of spatial gene expression of extant eleutherozoan echinoderms support homology and common ancestry of echinoderm larval skeletogenic cells. We propose larval skeletogenic cells arose in the stem lineage of eleutherozoans during a cell type duplication event that heterochronically activated adult skeletogenic cells in a topographically distinct tissue in early development.

A new ophiocistioid with soft-tissue preservation from the Silurian Herefordshire Lagerstätte, and the evolution of the holothurian body plan.

Reconstructing the evolutionary assembly of animal body plans is challenging when there are large morphological gaps between extant sister taxa, as in the case of echinozoans (echinoids and holothurians). However, the inclusion of extinct taxa can help bridge these gaps. Here we describe a new species of echinozoan, Sollasina cthulhu, from the Silurian Herefordshire Lagerstätte, UK. Sollasina cthulhu belongs to the ophiocistioids, an extinct group that shares characters with both echinoids and holothurians. Using physical-optical tomography and computer reconstruction, we visualize the internal anatomy of S. cthulhu in three dimensions, revealing inner soft tissues that we interpret as the ring canal, a key part of the water vascular system that was previously unknown in fossil echinozoans. Phylogenetic analyses strongly suggest that Sollasina and other ophiocistioids represent a paraphyletic group of stem holothurians, as previously hypothesized. This allows us to reconstruct the stepwise reduction of the skeleton during the assembly of the holothurian body plan, which may have been controlled by changes in the expression of biomineralization genes.

A new stem group echinoid from the Triassic of China leads to a revised macroevolutionary history of echinoids during the end-Permian mass extinction.

The Permian-Triassic bottleneck has long been thought to have drastically altered the course of echinoid evolution, with the extinction of the entire echinoid stem group having taken place during the end-Permian mass extinction. The Early Triassic fossil record of echinoids is, however, sparse, and new fossils are paving the way for a revised interpretation of the evolutionary history of echinoids during the Permian-Triassic crisis and Early Mesozoic. A new species of echinoid, Yunnanechinus luopingensis n. sp. recovered from the Middle Triassic (Anisian) Luoping Biota fossil Lagerstätte of South China, displays morphologies that are not characteristic of the echinoid crown group. We have used phylogenetic analyses to further demonstrate that Yunnanechinus is not a member of the echinoid crown group. Thus a clade of stem group echinoids survived into the Middle Triassic, enduring the global crisis that characterized the end-Permian and Early Triassic. Therefore, stem group echinoids did not go extinct during the Palaeozoic, as previously thought, and appear to have coexisted with the echinoid crown group for at least 23 million years. Stem group echinoids thus exhibited the Lazarus effect during the latest Permian and Early Triassic, while crown group echinoids did not.

Paleogenomics of echinoids reveals an ancient origin for the double-negative specification of micromeres in sea urchins.

Establishing a timeline for the evolution of novelties is a common, unifying goal at the intersection of evolutionary and developmental biology. Analyses of gene regulatory networks (GRNs) provide the ability to understand the underlying genetic and developmental mechanisms responsible for the origin of morphological structures both in the development of an individual and across entire evolutionary lineages. Accurately dating GRN novelties, thereby establishing a timeline for GRN evolution, is necessary to answer questions about the rate at which GRNs and their subcircuits evolve, and to tie their evolution to paleoenvironmental and paleoecological changes. Paleogenomics unites the fossil record and all aspects of deep time, with modern genomics and developmental biology to understand the evolution of genomes in evolutionary time. Recent work on the regulatory genomic basis of development in cidaroid echinoids, sand dollars, heart urchins, and other nonmodel echinoderms provides an ideal dataset with which to explore GRN evolution in a comparative framework. Using divergence time estimation and ancestral state reconstructions, we have determined the age of the double-negative gate (DNG), the subcircuit which specifies micromeres and skeletogenic cells in Strongylocentrotus purpuratus We have determined that the DNG has likely been used for euechinoid echinoid micromere specification since at least the Late Triassic. The innovation of the DNG thus predates the burst of post-Paleozoic echinoid morphological diversification that began in the Early Jurassic. Paleogenomics has wide applicability for the integration of deep time and molecular developmental data, and has wide utility in rigorously establishing timelines for GRN evolution.

Neither Hematocrit Normalization nor Exercise Training Restores Oxygen Consumption to Normal Levels in Hemodialysis Patients.

Patients treated with hemodialysis develop severely reduced functional capacity, which can be partially ameliorated by correcting anemia and through exercise training. In this study, we determined perturbations of an erythroid-stimulating agent and exercise training to examine if and where limitation to oxygen transport exists in patients on hemodialysis. Twenty-seven patients on hemodialysis completed a crossover study consisting of two exercise training phases at two hematocrit (Hct) values: 30% (anemic) and 42% (physiologic; normalized by treatment with erythroid-stimulating agent). To determine primary outcome measures of peak power and oxygen consumption (VO2) and secondary measures related to components of oxygen transport and utilization, all patients underwent numerous tests at five time points: baseline, untrained at Hct of 30%, after training at Hct of 30%, untrained at Hct of 42%, and after training at Hct of 42%. Hct normalization, exercise training, or the combination thereof significantly improved peak power and VO2 relative to values in the untrained anemic phase. Hct normalization increased peak arterial oxygen and arteriovenous oxygen difference, whereas exercise training improved cardiac output, citrate synthase activity, and peak tissue diffusing capacity. However, although the increase in arterial oxygen observed in the combination phase reached a value similar to that in healthy sedentary controls, the increase in peak arteriovenous oxygen difference did not. Muscle biopsy specimens showed markedly thickened endothelium and electron-dense interstitial deposits. In conclusion, exercise and Hct normalization had positive effects but failed to normalize exercise capacity in patients on hemodialysis. This effect may be caused by abnormalities identified within skeletal muscle.

Ancestral state reconstruction by comparative analysis of a GRN kernel operating in echinoderms.

Diverse sampling of organisms across the five major classes in the phylum Echinodermata is beginning to reveal much about the structure and function of gene regulatory networks (GRNs) in development and evolution. Sea urchins are the most studied clade within this phylum, and recent work suggests there has been dramatic rewiring at the top of the skeletogenic GRN along the lineage leading to extant members of the euechinoid sea urchins. Such rewiring likely accounts for some of the observed developmental differences between the two major subclasses of sea urchins-cidaroids and euechinoids. To address effects of topmost rewiring on downstream GRN events, we cloned four downstream regulatory genes within the skeletogenic GRN and surveyed their spatiotemporal expression patterns in the cidaroid Eucidaris tribuloides. We performed phylogenetic analyses with homologs from other non-vertebrate deuterostomes and characterized their spatiotemporal expression by quantitative polymerase chain reaction (qPCR) and whole-mount in situ hybridization (WMISH). Our data suggest the erg-hex-tgif subcircuit, a putative GRN kernel, exhibits a mesoderm-specific expression pattern early in Eucidaris development that is directly downstream of the initial mesodermal GRN circuitry. Comparative analysis of the expression of this subcircuit in four echinoderm taxa allowed robust ancestral state reconstruction, supporting hypotheses that its ancestral function was to stabilize the mesodermal regulatory state and that it has been co-opted and deployed as a unit in mesodermal subdomains in distantly diverged echinoderms. Importantly, our study supports the notion that GRN kernels exhibit structural and functional modularity, locking down and stabilizing clade-specific, embryonic regulatory states.

Reorganization of sea urchin gene regulatory networks at least 268 million years ago as revealed by oldest fossil cidaroid echinoid.

Echinoids, or sea urchins, are rare in the Palaeozoic fossil record, and thus the details regarding the early diversification of crown group echinoids are unclear. Here we report on the earliest probable crown group echinoid from the fossil record, recovered from Permian (Roadian-Capitanian) rocks of west Texas, which has important implications for the timing of the divergence of crown group echinoids. The presence of apophyses and rigidly sutured interambulacral areas with two columns of plates indicates this species is a cidaroid echinoid. The species, Eotiaris guadalupensis, n. sp. is therefore the earliest stem group cidaroid. The occurrence of this species in Roadian strata pushes back the divergence of cidaroids and euechinoids, the clades that comprise all living echinoids, to at least 268.8 Ma, ten million years older than the previously oldest known cidaroid. Furthermore, the genomic regulation of development in echinoids is amongst the best known, and this new species informs the timing of large-scale reorganization in echinoid gene regulatory networks that occurred at the cidaroid-euechinoid divergence, indicating that these changes took place by the Roadian stage of the Permian.

Juvenile skeletogenesis in anciently diverged sea urchin clades.

Mechanistic understanding of evolutionary divergence in animal body plans devolves from analysis of those developmental processes that, in forms descendant from a common ancestor, are responsible for their morphological differences. The last common ancestor of the two extant subclasses of sea urchins, i.e., euechinoids and cidaroids, existed well before the Permian/Triassic extinction (252 mya). Subsequent evolutionary divergence of these clades offers in principle a rare opportunity to solve the developmental regulatory events underlying a defined evolutionary divergence process. Thus (i) there is an excellent and fairly dense (if yet incompletely analyzed) fossil record; (ii) cladistically confined features of the skeletal structures of modern euechinoid and cidaroid sea urchins are preserved in fossils of ancestral forms; (iii) euechinoids and cidaroids are among current laboratory model systems in molecular developmental biology (here Strongylocentrotus purpuratus [Sp] and Eucidaris tribuloides [Et]); (iv) skeletogenic specification in sea urchins is uncommonly well understood at the causal level of interactions of regulatory genes with one another, and with known skeletogenic effector genes, providing a ready arsenal of available molecular tools. Here we focus on differences in test and perignathic girdle skeletal morphology that distinguish all modern euechinoid from all modern cidaroid sea urchins. We demonstrate distinct canonical test and girdle morphologies in juveniles of both species by use of SEM and X-ray microtomography. Among the sharply distinct morphological features of these clades are the internal skeletal structures of the perignathic girdle to which attach homologous muscles utilized for retraction and protraction of Aristotles׳ lantern and its teeth. We demonstrate that these structures develop de novo between one and four weeks after metamorphosis. In order to study the underlying developmental processes, a method of section whole mount in situ hybridization was adapted. This method displays current gene expression in the developing test and perignathic girdle skeletal elements of both Sp and Et juveniles. Active, specific expression of the sm37 biomineralization gene in these muscle attachment structures accompanies morphogenetic development of these clade-specific features in juveniles of both species. Skeletogenesis at these clade-specific muscle attachment structures displays molecular earmarks of the well understood embryonic skeletogenic GRN: thus the upstream regulatory gene alx1 and the gene encoding the vegfR signaling receptor are both expressed at the sites where they are formed. This work opens the way to analysis of the alternative spatial specification processes that were installed at the evolutionary divergence of the two extant subclasses of sea urchins.

Evaluation of methods to estimate understory fruit biomass.

Fleshy fruit is consumed by many wildlife species and is a critical component of forest ecosystems. Because fruit production may change quickly during forest succession, frequent monitoring of fruit biomass may be needed to better understand shifts in wildlife habitat quality. Yet, designing a fruit sampling protocol that is executable on a frequent basis may be difficult, and knowledge of accuracy within monitoring protocols is lacking. We evaluated the accuracy and efficiency of 3 methods to estimate understory fruit biomass (Fruit Count, Stem Density, and Plant Coverage). The Fruit Count method requires visual counts of fruit to estimate fruit biomass. The Stem Density method uses counts of all stems of fruit producing species to estimate fruit biomass. The Plant Coverage method uses land coverage of fruit producing species to estimate fruit biomass. Using linear regression models under a censored-normal distribution, we determined the Fruit Count and Stem Density methods could accurately estimate fruit biomass; however, when comparing AIC values between models, the Fruit Count method was the superior method for estimating fruit biomass. After determining that Fruit Count was the superior method to accurately estimate fruit biomass, we conducted additional analyses to determine the sampling intensity (i.e., percentage of area) necessary to accurately estimate fruit biomass. The Fruit Count method accurately estimated fruit biomass at a 0.8% sampling intensity. In some cases, sampling 0.8% of an area may not be feasible. In these cases, we suggest sampling understory fruit production with the Fruit Count method at the greatest feasible sampling intensity, which could be valuable to assess annual fluctuations in fruit production.

Hyperkalemia in free-ranging white-tailed deer (Odocoileus virginianus).

Sixty adult and yearling female white-tailed deer (Odocoileus virginianus) were collected in July 2008 (n=30) and March 2009 (n=30) from eastern North Carolina as part of a population health assessment. During July 2008, standard serum analyses revealed hyperkalemia in all deer sampled. In March, the effect of processing time as a possible source of the hyperkalemia was investigated. For a subset of deer (n=10), blood tubes were centrifuged and processed at four time points (0, 30, 60, and 120 min) postcollection. Delayed centrifugation and plasma separation did not affect potassium (K(+)) concentration over time, indicating that a shift in intracellular K(+) did not occur and the hyperkalemia was not due to improper sample handling. Potassium levels were negatively correlated with age and varied across collection periods. Also, K(+) levels were positively correlated with glucose and not correlated with creatine kinase (CK). No single variable indicated a strong enough relationship to explain the hyperkalemia in the study.

A population-based study of the effects of birth weight on early developmental delay or disability in children.

Improving medical treatment of extremely low-birth-weight infants over the last 20 to 30 years resulted in increased survival rates. The developmental sequela of salvaged infants is of great interest to perinatologists. The primary purposes of the current study were to assess the effect of birth weight (BW) on developmental delay or disability (DDD) in the first three years of life and determine whether there is a BW threshold below which all infants should be evaluated to determine if intervention services for children with DDD should be received. Three statewide databases were merged: 1998 Birth Vital Statistics; 1997-1998 Medicaid eligibility files; and 1998-2001 Children's Medical Services' Early Intervention Program (CMS-EIP) data. Infants who died within the first year of life and plural births were excluded. The final dataset consisted of 170,874 records. A child was determined to have a DDD if a developmental delay, or an established condition, such as sensory impairment, genetic, metabolic, neurological, or severe attachment disorders, was diagnosed through a multidisciplinary evaluation. Logistic regression models were used to relate BW to DDD, controlling for sociodemographic, behavioral, and perinatal variables. Adjusted odds ratios (OR) were calculated to describe the effects of BW on DDD. There was a significant effect of BW on DDD (Adjusted OR &equals 97.50, 40.01, 15.84, 3.29, 1.39, 1.00, 1.52 for BW categories 450-749, 750-999, 1000-1499, 1500- 2499, 2500-2999, 3000-4749, 4750-6050 g, respectively). In these categories, 70%, 56%, 36%, 11%, 4%, 3%, and 6% of surviving singleton infants, respectively, suffered a DDD in their first 3 years of life. Four medical, five sociodemographic, and two behavioral factors were significant in addition to BW. An equation for predicting the probability of DDD given these factors was obtained, and its use exemplified. BW is strongly associated with DDD. Over 60% of infants weighing < 1000 g and nearly half (46%) of those weighing < 1500 g at birth are diagnosed with a DDD before 3 years of age. The probability of DDD for a specific infant also varies by sociodemographic, other perinatal, and behavioral factors. The results of this paper suggest that all surviving infants of BW < 1000 g, and perhaps < 1500 g, should be automatically referred for evaluation.