David A. Weese

Phylogenomic Resolution of the Hemichordate and Echinoderm Clade

Ambulacraria, comprising Hemichordata and Echinodermata, is closely related to Chordata, making it integral to understanding chordate origins and polarizing chordate molecular and morphological characters. Unfortunately, relationships within Hemichordata and Echinodermata have remained unresolved, compromising our ability to extrapolate findings from the most closely related molecular and developmental models outside of Chordata (e.g., the acorn worms Saccoglossus kowalevskii and Ptychodera flava and the sea urchin Strongylocentrotus purpuratus). To resolve long-standing phylogenetic issues within Ambulacraria, we sequenced transcriptomes for 14 hemichordates as well as 8 echinoderms and complemented these with existing data for a total of 33 ambulacrarian operational taxonomic units (OTUs). Examination of leaf stability values revealed rhabdopleurid pterobranchs and the enteropneust Stereobalanus canadensis were unstable in placement; therefore, analyses were also run without these taxa. Analyses of 185 genes resulted in reciprocal monophyly of Enteropneusta and Pterobranchia, placed the deep-sea family Torquaratoridae within Ptychoderidae, and confirmed the position of ophiuroid brittle stars as sister to asteroid sea stars (the Asterozoa hypothesis). These results are consistent with earlier perspectives concerning plesiomorphies of Ambulacraria, including pharyngeal gill slits, a single axocoel, and paired hydrocoels and somatocoels. The resolved ambulacrarian phylogeny will help clarify the early evolution of chordate characteristics and has implications for our understanding of major fossil groups, including graptolites and somasteroideans.

Phylogenomics of Lophotrochozoa with consideration of systematic error

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

Rocks and clocks: linking geologic history and rates of genetic differentiation in anchialine organisms

The geologic history of a region can significantly impact the development of its flora and fauna, with past events shaping community patterns and evolutionary trajectories of species. In this context, islands are excellent “natural laboratories” for studying the fundamental processes of evolution due to their discrete geographical nature and dynamic geologic histories. An island system meeting these criteria is the Hawaiian Archipelago, which is ideal for testing how island geologic history influences the processes leading to population genetic variation and differentiation. One Hawaiian endemic whose evolutionary history is closely tied to the geology of the islands is the anchialine atyid shrimp Halocaridina, whose mitochondrial cytochrome oxidase I (COI) gene is hypothesized to be evolving at the rate of 20% per million years. To validate this rapid evolutionary rate, time since divergence estimates between geographically close, yet genetically distinct, populations were calculated for Halocaridina from anchialine habitats on the islands of Hawai’i, Maui, and O’ahu. On the younger (i.e., <1.5 million years) islands of Hawai’i and Maui, where all anchialine habitats occur in basalt, application of the Halocaridina molecular clock identified a strong correlation between levels of genetic divergence and the geologic age of the region inhabited by those populations. In contrast, this relationship weakened when similar analyses were conducted for Halocaridina from limestone anchialine habitats on the older (i.e., >2.75 million years) island of O’ahu. These results suggest geologic age, basin origin and/or composition are important factors that should be taken into consideration when conducting molecular clock analyses on anchialine flora and fauna as well as island populations in general.

THE LONG AND SHORT OF IT: GENETIC VARIATION AND POPULATION STRUCTURE OF THE ANCHIALINE ATYID SHRIMP CARIDINA RUBELLA ON MIYAKO-JIMA, JAPAN

ABSTRACT One of the most threatened ecosystems on many islands may be anchialine habitats, or coastal land-locked water bodies with no surface connection to the sea yet containing brackish water that fluctuates with the tides. To better manage these habitats, it is important to develop a broader understanding of the biodiversity within them since such knowledge plays critical roles in establishing conservation strategies. In this study, the genetic variation and population structure of an anchialine atyid shrimp, Caridina rubella Fujino and Shokita, 1975 was investigated in the Southern Ryukyu Archipelago, Japan. Given a planktotrophic larval stage and its potential amphidromous life cycle, populations of C. rubella on the island of Miyako-jima (to which it is apparently restricted) were hypothesized to have little to no structure across the island. To test this, 61 individuals were collected from four anchialine caves and sequence variation examined at the cytochrome c oxidase subunit I (COI) region of the mitochondrial DNA (mtDNA). Surprisingly, significant genetic structure was exhibited across distances ranging from <20 m to > 10 km. Additionally, deep (∼17% p-distance) genetic divergence correlating with distinct variation in rostrum lengths was found between closely situated, but more-or-less completely isolated, populations. This implies “C. rubella” may actually represent two distinct species on Miyako-jima. Given that this atyid is already listed as a threatened species by the Japanese government, the results presented here are useful in the formulation and implementation of future conservation plans for populations of C. rubella on Miyako-jima.

Demographic response of cutlassfish ( Trichiurus japonicus and T. nanhaiensis ) to fluctuating palaeo-climate and regional oceanographic conditions in the China seas

Glacial cycles of the Quaternary have heavily influenced the demographic history of various species. To test the evolutionary impact of palaeo-geologic and climatic events on the demographic history of marine taxa from the coastal Western Pacific, we investigated the population structure and demographic history of two economically important fish (Trichiurus japonicus and T. nanhaiensis) that inhabit the continental shelves of the East China and northern South China Seas using the mitochondrial cytochrome b sequences and Bayesian Skyline Plot analyses. A molecular rate of 2.03% per million years, calibrated to the earliest flooding of the East China Sea shelf (70–140 kya), revealed a strong correlation between population sizes and primary production. Furthermore, comparison of the demographic history of T. japonicus populations from the East China and South China Seas provided evidence of the postglacial development of the Changjiang (Yangtze River) Delta. In the South China Sea, interspecific comparisons between T. japonicus and T. nanhaiensis indicated possible evolutionary responses to changes in palaeo-productivity that were influenced by East Asian winter monsoons. This study not only provides insight into the demographic history of cutlassfish but also reveals potential clues regarding the historic productivity and regional oceanographic conditions of the Western Pacific marginal seas.

Multiple colonizations lead to cryptic biodiversity in an island ecosystem: comparative phylogeography of anchialine shrimp species in the Ryukyu Archipelago, Japan.

Archipelagos of the Indo-West Pacific are considered to be among the richest in the world in biodiversity, and phylogeographic studies generally support either the center of origin or the center of accumulation hypothesis to explain this pattern. To differentiate between these competing hypotheses for organisms from the Indo-West Pacific anchialine ecosystem, defined as coastal bodies of mixohaline water fluctuating with the tides but having no direct oceanic connections, we investigated the genetic variation, population structure, and evolutionary history of three caridean shrimp species (Antecaridina lauensis, Halocaridinides trigonophthalma, and Metabetaeus minutus) in the Ryukyu Archipelago, Japan. We used two mitochondrial genes—cytochrome c oxidase subunit I (COI) and large ribosomal subunit (16S-rDNA)—complemented with genetic examination of available specimens from the same or closely related species from the Indian and Pacific Oceans. In the Ryukyus, each species encompassed 2–3 divergent (9.52%–19.2% COI p-distance) lineages, each having significant population structure and varying geographic distributions. Phylogenetically, the A. lauensis and M. minutus lineages in the Ryukyus were more closely related to ones from outside the archipelago than to one another. These results, when interpreted in the context of Pacific oceanographic currents and geologic history of the Ryukyus, imply multiple colonizations of the archipelago by the three species, consistent with the center of accumulation hypothesis. While this study contributes toward understanding the biodiversity, ecology, and evolution of organisms in the Ryukyus and the Indo-West Pacific, it also has potential utility in establishing conservation strategies for anchialine fauna of the Pacific Basin in general.

Reproduction and Development in Halocaridina rubra Holthuis, 1963 (Crustacea: Atyidae) Clarifies Larval Ecology in the Hawaiian Anchialine Ecosystem

Larvae in aquatic habitats often develop in environments different from those they inhabit as adults. Shrimp in the Atyidae exemplify this trend, as larvae of many species require salt or brackish water for development, while adults are freshwater-adapted. An exception within the Atyidae family is the “anchialine clade,” which are euryhaline as adults and endemic to habitats with subterranean fresh and marine water influences. Although the Hawaiian anchialine atyid Halocaridina rubra is a strong osmoregulator, its larvae have never been observed in nature. Moreover, larval development in anchialine species is poorly studied. Here, reproductive trends in laboratory colonies over a 5-y period are presented from seven genetic lineages and one mixed population of H. rubra; larval survivorship under varying salinities is also discussed. The presence and number of larvae differed significantly among lineages, with the mixed population being the most prolific. Statistical differences in reproduction attributable to seasonality also were identified. Larval survivorship was lowest (12% settlement rate) at a salinity approaching fresh water and significantly higher in brackish and seawater (88% and 72%, respectively). Correlated with this finding, identifiable gills capable of ion transport did not develop until metamorphosis into juveniles. Thus, early life stages of H. rubra are apparently excluded from surface waters, which are characterized by lower and fluctuating salinities. Instead, these stages are restricted to the subterranean (where there is higher and more stable salinity) portion of Hawaii’s anchialine habitats due to their inability to tolerate low salinities. Taken together, these data contribute to the understudied area of larval ecology in the anchialine ecosystem.

Reconstruction of Cyclooxygenase Evolution in Animals Suggests Variable, Lineage-Specific Duplications, and Homologs with Low Sequence Identity

Cyclooxygenase (COX) enzymatically converts arachidonic acid into prostaglandin G/H in animals and has importance during pregnancy, digestion, and other physiological functions in mammals. COX genes have mainly been described from vertebrates, where gene duplications are common, but few studies have examined COX in invertebrates. Given the increasing ease in generating genomic data, as well as recent, although incomplete descriptions of potential COX sequences in Mollusca, Crustacea, and Insecta, assessing COX evolution across Metazoa is now possible. Here, we recover 40 putative COX orthologs by searching publicly available genomic resources as well as ~250 novel invertebrate transcriptomic datasets. Results suggest the common ancestor of Cnidaria and Bilateria possessed a COX homolog similar to those of vertebrates, although such homologs were not found in poriferan and ctenophore genomes. COX was found in most crustaceans and the majority of molluscs examined, but only specific taxa/lineages within Cnidaria and Annelida. For example, all octocorallians appear to have COX, while no COX homologs were found in hexacorallian datasets. Most species examined had a single homolog, although species-specific COX duplications were found in members of Annelida, Mollusca, and Cnidaria. Additionally, COX genes were not found in Hemichordata, Echinodermata, or Platyhelminthes, and the few previously described COX genes in Insecta lacked appreciable sequence homology (although structural analyses suggest these may still be functional COX enzymes). This analysis provides a benchmark for identifying COX homologs in future genomic and transcriptomic datasets, and identifies lineages for future studies of COX.

Looking for Needles in a Haystack: Molecular Identification of Anchialine Crustacean Larvae (Decapoda: Caridea) From the Shiokawa Spring, Okinawa Island, Ryukyu Islands, Japan

Little is known about the larval ecology of organisms from the anchialine ecosystem, tidally influenced land-locked bodies of mixohaline water possessing simultaneous subterranean connections to the ocean and groundwater aquifer. Larval stages have not been observed in the epigeal (surface) portion of such environments, implying reproduction and development is restricted to the hypogeal (subterranean) component. To test this hypothesis, larvae were collected from the hypogeal-derived waters of Shiokawa Spring on Okinawa Island, Ryukyu Islands, Japan and species identity determined through analyses of sequence variation in the mitochondrial cytochrome oxidase subunit I (COI) gene. The majority of larvae were identified as belonging to one of two species of anchialine atyid shrimps ( Caridina rubella Fujino and Shokita, 1975 and Halocaridinides trigonophthalma Fujino and Shokita, 1975), as well as a palaemonid shrimp ( Macrobrachium miyakoense Fujino and Shokita, 1975) previously reported as being endemic to a single island (Miyako) in the Ryukyus. This study represents the first report of anchialine caridean larvae being collected from the epigeal portion of the anchialine habitat and supports the long-standing hypothesis that reproduction and development, as well as the larval habitat, of many anchialine taxa occurs in, and is, the hypogeal component of the ecosystem. The high species and genetic diversity among the sampled larvae was unexpected given the strong population structure previously reported for anchialine caridean taxa in the Ryukyus. These findings suggest a reproductively active and underappreciated anchialine community inhabiting the hypogeal waters of the island and highlight the significance of Okinawa in the context of the phylogeography of anchialine taxa in the Ryukyus.

Diversity and the environmental drivers of spatial variation in Bacteria and micro-Eukarya communities from the Hawaiian anchialine ecosystem

Little is currently known regarding microbial community structure, and the environmental factors influencing it, within the anchialine ecosystem, defined as near-shore, land-locked water bodies with subsurface connections to the ocean and groundwater aquifer. The Hawaiian Archipelago is home to numerous anchialine habitats, with some on the islands of Maui and Hawaii harboring unique, laminated orange cyanobacterial–bacterial crusts that independently assembled in relatively young basalt fields. Here, benthic and water column bacterial and micro-eukaryotic communities from nine anchialine habitats on Oahu, Maui, and Hawaii were surveyed using high-throughput amplicon sequencing of the V6 (Bacteria-specific) and V9 (Eukarya-biased) hypervariable regions of the 16S- and 18S-rDNA genes, respectively. While benthic communities from habitats with cyanobacterial–bacterial crusts were more similar to each other than to ones lacking it on the same island, each habitat had distinct benthic and water column microbial communities. Analyses of the survey data in the context of environmental factors identified salinity, site, aquifer, and watershed as having the highest explanatory power for the observed variation in microbial diversity and community structure, with lesser drivers being annual rainfall, longitude, ammonium, and dissolved organic carbon. Our results epitomize the abiotic and biotic uniqueness characteristic of individual habitats comprising the Hawaiian anchialine ecosystem.