Independent Innexin Radiation Shaped Signaling in Ctenophores.

Innexins facilitate cell-cell communication by forming gap junctions or nonjunctional hemichannels, which play important roles in metabolic, chemical, ionic, and electrical coupling. The lack of knowledge regarding the evolution and role of these channels in ctenophores (comb jellies), the likely sister group to the rest of animals, represents a substantial gap in our understanding of the evolution of intercellular communication in animals. Here, we identify and phylogenetically characterize the complete set of innexins of four ctenophores: Mnemiopsis leidyi, Hormiphora californensis, Pleurobrachia bachei, and Beroe ovata. Our phylogenetic analyses suggest that ctenophore innexins diversified independently from those of other animals and were established early in the emergence of ctenophores. We identified a four-innexin genomic cluster, which was present in the last common ancestor of these four species and has been largely maintained in these lineages. Evidence from correlated spatial and temporal gene expression of the M. leidyi innexin cluster suggests that this cluster has been maintained due to constraints related to gene regulation. We describe the basic electrophysiological properties of putative ctenophore hemichannels from muscle cells using intracellular recording techniques, showing substantial overlap with the properties of bilaterian innexin channels. Together, our results suggest that the last common ancestor of animals had gap junctional channels also capable of forming functional innexin hemichannels, and that innexin genes have independently evolved in major lineages throughout Metazoa.

A reference genome for ecological restoration of the sunflower sea star, Pycnopodia helianthoides.

Wildlife diseases, such as the sea star wasting (SSW) epizootic that outbroke in the mid-2010s, appear to be associated with acute and/or chronic abiotic environmental change; dissociating the effects of different drivers can be difficult. The sunflower sea star, Pycnopodia helianthoides, was the species most severely impacted during the SSW outbreak, which overlapped with periods of anomalous atmospheric and oceanographic conditions, and there is not yet a consensus on the cause(s). Genomic data may reveal underlying molecular signatures that implicate a subset of factors and, thus, clarify past events while also setting the scene for effective restoration efforts. To advance this goal, we used Pacific Biosciences HiFi long sequencing reads and Dovetail Omni-C proximity reads to generate a highly contiguous genome assembly that was then annotated using RNA-seq-informed gene prediction. The genome assembly is 484 Mb long, with contig N50 of 1.9 Mb, scaffold N50 of 21.8 Mb, BUSCO completeness score of 96.1%, and 22 major scaffolds consistent with prior evidence that sea star genomes comprise 22 autosomes. These statistics generally fall between those of other recently assembled chromosome-scale assemblies for two species in the distantly related asteroid genus Pisaster. These novel genomic resources for P. helianthoides will underwrite population genomic, comparative genomic, and phylogenomic analyses-as well as their integration across scales-of SSW and environmental stressors.

Comparative transcriptomics reveals altered species interaction between the bioeroding sponge Cliona varians and the coral Porites furcata under ocean acidification.

Bioeroding sponges interact and compete with corals on tropical reefs. Experimental studies have shown global change alters this biotic interaction, often in favour of the sponge. Ocean acidification in particular increases sponge bioerosion and reduces coral calcification, yet little is known about the molecular basis of these changes. We used RNA-Seq data to understand how acidification impacts the interaction between the bioeroding sponge, Cliona varians, and the coral, Porites furcata, at the transcriptomic level. Replicate sponge and coral genets were exposed to ambient (8.1 pH) and acidified (7.6 pH) conditions in isolation and in treatments where they were joined for 48 h. The coral had a small gene expression response (tens of transcripts) to the sponge, suggesting it does little at the transcriptomic level to deter sponge overgrowth. By contrast, the sponge differentially expressed 7320 transcripts in response to the coral under ambient conditions and 3707 transcripts in response to acidification. Overlap in the responses to acidification and the coral, 2500 transcripts expressed under both treatments, suggests a similar physiological response to both cues. The sponge expressed 50× fewer transcripts in response to the coral under acidification, suggesting energetic costs of bioerosion, and other cellular processes, are lower for sponges under acidification. Our results suggest how acidification drives ecosystem-level changes in the accretion/bioerosion balance on coral reefs. This shift is not only the result of changes to the thermodynamic balance of these chemical reactions but also the result of active physiological responses of organisms to each other and their abiotic environment.

A chromosome-level reference genome for the giant pink sea star, Pisaster brevispinus, a species severely impacted by wasting.

Efforts to protect the ecologically and economically significant California Current Ecosystem from global change will greatly benefit from data about patterns of local adaptation and population connectivity. To facilitate that work, we present a reference-quality genome for the giant pink sea star, Pisaster brevispinus, a species of ecological importance along the Pacific west coast of North America that has been heavily impacted by environmental change and disease. We used Pacific Biosciences HiFi long sequencing reads and Dovetail Omni-C proximity reads to generate a highly contiguous genome assembly of 550 Mb in length. The assembly contains 127 scaffolds with a contig N50 of 4.6 Mb and a scaffold N50 of 21.4 Mb; the BUSCO completeness score is 98.70%. The P. brevispinus genome assembly is comparable to the genome of the congener species P. ochraceus in size and completeness. Both Pisaster assemblies are consistent with previously published karyotyping results showing sea star genomes are organized into 22 autosomes. The reference genome for P. brevispinus is an important first step toward the goal of producing a comprehensive, population genomics view of ecological and evolutionary processes along the California coast. This resource will help scientists, managers, and policy makers in their task of understanding and protecting critical coastal regions from the impacts of global change.

Reference genome for the California ribbed mussel, Mytilus californianus, an ecosystem engineer.

The California ribbed mussel, Mytilus californianus, is an ecosystem engineer crucial for the survival of many marine species inhabiting the intertidal zone of California. Here, we describe the first reference genome for M. californianus and compare it to previously published genomes from three other Mytilus species: M. edulis, M. coruscus, and M. galloprovincialis. The M. californianus reference genome is 1.65 Gb in length, with N50 sequence length of 118 Mb, and an estimated 86.0% complete single copy genes. Compared with the other three Mytilus species, the M. californianus genome assembly is the longest, has the highest N50 value, and the highest percentage complete single copy genes. This high-quality genome assembly provides a foundation for population genetic analyses that will give insight into future conservation work along the coast of California.

ANISEED 2019: 4D exploration of genetic data for an extended range of tunicates.

ANISEED (https://www.aniseed.cnrs.fr) is the main model organism database for the worldwide community of scientists working on tunicates, the vertebrate sister-group. Information provided for each species includes functionally-annotated gene and transcript models with orthology relationships within tunicates, and with echinoderms, cephalochordates and vertebrates. Beyond genes the system describes other genetic elements, including repeated elements and cis-regulatory modules. Gene expression profiles for several thousand genes are formalized in both wild-type and experimentally-manipulated conditions, using formal anatomical ontologies. These data can be explored through three complementary types of browsers, each offering a different view-point. A developmental browser summarizes the information in a gene- or territory-centric manner. Advanced genomic browsers integrate the genetic features surrounding genes or gene sets within a species. A Genomicus synteny browser explores the conservation of local gene order across deuterostome. This new release covers an extended taxonomic range of 14 species, including for the first time a non-ascidian species, the appendicularian Oikopleura dioica. Functional annotations, provided for each species, were enhanced through a combination of manual curation of gene models and the development of an improved orthology detection pipeline. Finally, gene expression profiles and anatomical territories can be explored in 4D online through the newly developed Morphonet morphogenetic browser.

Phylotocol: Promoting Transparency and Overcoming Bias in Phylogenetics.

The integrity of science requires that the process be based on sound experimental design and objective methodology. Strategies that increase reproducibility and transparency in science protect this integrity by reducing conscious and unconscious biases. Given the large number of analysis options and the constant development of new methodologies in phylogenetics, this field is one that would particularly benefit from more transparent research design. Herein, we introduce phylotocol (fi lo 'ta kôl), an a priori protocol-driven approach in which all analyses are planned and documented at the start of a project. The phylotocol template is simple and the implementation options are flexible to reduce administrative burdens and allow researchers to adapt it to their needs without restricting scientific creativity. While the primary goal of phylotocol is to increase transparency and accountability, it has a number of auxiliary benefits including improving study design and reproducibility, enhancing collaboration and education, and increasing the likelihood of project completion. Our goal with this Point of View article is to encourage a dialog about transparency in phylogenetics and the best strategies to bring transparent research practices to our field.

Integrating Embryonic Development and Evolutionary History to Characterize Tentacle-Specific Cell Types in a Ctenophore.

The origin of novel traits can promote expansion into new niches and drive speciation. Ctenophores (comb jellies) are unified by their possession of a novel cell type: the colloblast, an adhesive cell found only in the tentacles. Although colloblast-laden tentacles are fundamental for prey capture among ctenophores, some species have tentacles lacking colloblasts and others have lost their tentacles completely. We used transcriptomes from 36 ctenophore species to identify gene losses that occurred specifically in lineages lacking colloblasts and tentacles. We cross-referenced these colloblast- and tentacle-specific candidate genes with temporal RNA-Seq during embryogenesis in Mnemiopsis leidyi and found that both sets of candidates are preferentially expressed during tentacle morphogenesis. We also demonstrate significant upregulation of candidates from both data sets in the tentacle bulb of adults. Both sets of candidates were enriched for an N-terminal signal peptide and protein domains associated with secretion; among tentacle candidates we also identified orthologs of cnidarian toxin proteins, presenting tantalizing evidence that ctenophore tentacles may secrete toxins along with their adhesive. Finally, using cell lineage tracing, we demonstrate that colloblasts and neurons share a common progenitor, suggesting the evolution of colloblasts involved co-option of a neurosecretory gene regulatory network. Together these data offer an initial glimpse into the genetic architecture underlying ctenophore cell-type diversity.

Phenotypic and transcriptomic responses to salinity stress across genetically and geographically divergent Tigriopus californicus populations.

Species inhabiting the North American west coast intertidal must tolerate an extremely variable environment, with large fluctuations in both temperature and salinity. Uncovering the mechanisms for this tolerance is key to understanding species' persistence. We tested for differences in salinity tolerance between populations of Tigriopus californicus copepods from locations in northern (Bodega Reserve) and southern (San Diego) California known to differ in temperature, precipitation and humidity. We also tested for differences between populations in their transcriptomic responses to salinity. Although these two populations have ~20% mtDNA sequence divergence and differ strongly in other phenotypic traits, we observed similarities in their phenotypic and transcriptomic responses to low and high salinity stress. Salinity significantly affected respiration rate (increased under low salinity and reduced under high salinity), but we found no significant effect of population on respiration or a population by salinity interaction. Under high salinity, there was no population difference in knock-down response, but northern copepods had a smaller knock-down under low salinity stress, corroborating previous results for T. californicus. Northern and southern populations had a similar transcriptomic response to salinity based on a principle components analysis, although differential gene expression under high salinity stress was three times lower in the northern population compared to the southern population. Transcripts differentially regulated under salinity stress were enriched for "amino acid transport" and "ion transport" annotation categories, supporting previous work demonstrating that the accumulation of free amino acids is important for osmotic regulation in T. californicus.

Plastic and Evolved Responses to Global Change: What Can We Learn from Comparative Transcriptomics?

Physiological plasticity and adaptive evolution may facilitate persistence in a changing environment. As a result, there is an interest in understanding species' capacities for plastic and evolved responses, and the mechanisms by which these responses occur. Transcriptome sequencing has become a powerful tool for addressing these questions, providing insight into otherwise unobserved effects of changing conditions on organismal physiology and variation in these effects among individuals and populations. Here, we review recent studies using comparative transcriptomics to understand plastic and evolutionary responses to changing environments. We focus on 2 areas where transcriptomics has played an important role: first, in understanding the genetic basis for local adaptation to current gradients as a proxy for future adaptation, and second, in understanding organismal responses to multiple stressors. We find most studies examining multiple stressors have tested the effects of each stressor individually; the few studies testing multiple stressors simultaneously have found synergistic effects on gene expression that would not have been predicted from single stressor studies. We discuss the importance of robust experimental design to allow for a more sophisticated characterization of transcriptomic responses and conclude by offering recommendations for future research, including integrating genomics with transcriptomics, testing gene regulatory networks, and comparing the equivalence of transcription to translation and the effects of environmental stress on the proteome.

Mitochondrial DNA Variation in Southeastern Pre-Columbian Canids.

The taxonomic status of the red wolf (Canis rufus) is heavily debated, but could be clarified by examining historic specimens from the southeastern United States. We analyzed mitochondrial DNA (mtDNA) from 3 ancient (350-1900 year olds) putative wolf samples excavated from middens and sinkholes within the historic red wolf range. We detected 3 unique mtDNA haplotypes, which grouped with the coyote mtDNA clade, suggesting that the canids inhabiting southeastern North America prior to human colonization from Europe were either coyotes, which would vastly expand historic coyote distributions, an ancient coyote-wolf hybrid, or a North American evolved red wolf lineage related to coyotes. Should the red wolf prove to be a distinct species, our results support the idea of either an ancient hybrid origin for red wolves or a shared common ancestor between coyotes and red wolves.

Evaluating summary statistics used to test for incomplete lineage sorting: mito-nuclear discordance in the reef sponge Callyspongia vaginalis.

Conflicting patterns of population differentiation between the mitochondrial and nuclear genomes (mito-nuclear discordance) have become increasingly evident as multilocus data sets have become easier to generate. Incomplete lineage sorting (ILS) of nucDNA is often implicated as the cause of such discordance, stemming from the large effective population size of nucDNA relative to mtDNA. However, selection, sex-biased dispersal and historical demography can also lead to mito-nuclear discordance. Here, we compare patterns of genetic diversity and subdivision for six nuclear protein-coding gene regions to those for mtDNA in a common Caribbean coral reef sponge, Callyspongia vaginalis, along the Florida reef tract. We also evaluated a suite of summary statistics to determine which are effective metrics for comparing empirical and simulated data when testing drivers of mito-nuclear discordance in a statistical framework. While earlier work revealed three divergent and geographically subdivided mtDNACOI haplotypes separated by 2.4% sequence divergence, nuclear alleles were admixed with respect to mitochondrial clade and geography. Bayesian analysis showed that substitution rates for the nuclear loci were up to 7 times faster than for mitochondrial COI. Coalescent simulations and neutrality tests suggested that mito-nuclear discordance in C. vaginalis is not the result of ILS in the nucDNA or selection on the mtDNA but is more likely caused by changes in population size. Sperm-mediated gene flow may also influence patterns of population subdivision in the nucDNA.

Morphological and dietary changes encoded in the genome of Beroe ovata, a ctenophore-eating ctenophore.

As the sister group to all other animals, ctenophores (comb jellies) are important for understanding the emergence and diversification of numerous animal traits. Efforts to explore the evolutionary processes that promoted diversification within Ctenophora are hindered by undersampling genomic diversity within this clade. To address this gap, we present the sequence, assembly and initial annotation of the genome of Beroe ovata. Beroe possess unique morphology, behavior, ecology and development. Unlike their generalist carnivorous kin, beroid ctenophores feed exclusively on other ctenophores. Accordingly, our analyses revealed a loss of chitinase, an enzyme critical for the digestion of most non-ctenophore prey, but superfluous for ctenophorivores. Broadly, our genomic analysis revealed that extensive gene loss and changes in gene regulation have shaped the unique biology of B. ovata. Despite the gene losses in B. ovata, our phylogenetic analyses on photosensitive opsins and several early developmental regulatory genes show that these genes are conserved in B. ovata. This additional sampling contributes to a more complete reconstruction of the ctenophore ancestor and points to the need for extensive comparisons within this ancient and diverse clade of animals. To promote further exploration of these data, we present BovaDB (http://ryanlab.whitney.ufl.edu/bovadb/), a portal for the B. ovata genome.

Inferring Tunicate Relationships and the Evolution of the Tunicate Hox Cluster with the Genome of Corella inflata.

Tunicates, the closest living relatives of vertebrates, have served as a foundational model of early embryonic development for decades. Comparative studies of tunicate phylogeny and genome evolution provide a critical framework for analyzing chordate diversification and the emergence of vertebrates. Toward this goal, we sequenced the genome of Corella inflata (Ascidiacea, Phlebobranchia), so named for the capacity to brood self-fertilized embryos in a modified, "inflated" atrial chamber. Combining the new genome sequence for Co. inflata with publicly available tunicate data, we estimated a tunicate species phylogeny, reconstructed the ancestral Hox gene cluster at important nodes in the tunicate tree, and compared patterns of gene loss between Co. inflata and Ciona robusta, the prevailing tunicate model species. Our maximum-likelihood and Bayesian trees estimated from a concatenated 210-gene matrix were largely concordant and showed that Aplousobranchia was nested within a paraphyletic Phlebobranchia. We demonstrated that this relationship is not an artifact due to compositional heterogeneity, as had been suggested by previous studies. In addition, within Thaliacea, we recovered Doliolida as sister to the clade containing Salpida and Pyrosomatida. The Co. inflata genome provides increased resolution of the ancestral Hox clusters of key tunicate nodes, therefore expanding our understanding of the evolution of this cluster and its potential impact on tunicate morphological diversity. Our analyses of other gene families revealed that several cardiovascular associated genes (e.g., BMP10, SCL2A12, and PDE2a) absent from Ci. robusta, are present in Co. inflata. Taken together, our results help clarify tunicate relationships and the genomic content of key ancestral nodes within this phylogeny, providing critical insights into tunicate evolution.

Population Genomic Analyses of the Sea Urchin Echinometra sp. EZ across an Extreme Environmental Gradient.

Extreme environmental gradients represent excellent study systems to better understand the variables that mediate patterns of genomic variation between populations. They also allow for more accurate predictions of how future environmental change might affect marine species. The Persian/Arabian Gulf is extreme in both temperature and salinity, whereas the adjacent Gulf of Oman has conditions more typical of tropical oceans. The sea urchin Echinometra sp. EZ inhabits both of these seas and plays a critical role in coral reef health as a grazer and bioeroder, but, to date, there have been no population genomic studies on this or any urchin species in this unique region. E sp. EZ's life history traits (e.g., large population sizes, large reproductive clutches, and long life spans), in theory, should homogenize populations unless nonneutral processes are occurring. Here, we generated a draft genome and a restriction site-associated DNA sequencing data set from seven populations along an environmental gradient across the Persian/Arabian Gulf and the Gulf of Oman. The estimated genome size of E. sp. EZ was 609 Mb and the heterozygosity was among the highest recorded for an echinoderm at 4.5%. We recovered 918 high-quality SNPs from 85 individuals which we then used in downstream analyses. Population structure analyses revealed a high degree of admixture between all sites, although there was population differentiation and significant pairwise FST values between the two seas. Preliminary results suggest migration is bidirectional between the seas and nine candidate loci were identified as being under putative natural selection, including one collagen gene. This study is the first to investigate the population genomics of a sea urchin from this extreme environmental gradient and is an important contribution to our understanding of the complex spatial patterns that drive genomic divergence.

Sponges Lack ParaHox Genes.

Addressing the origin of axial-patterning machinery is essential for understanding the evolution of animal form. Historically, sponges, a lineage that branched off early in animal evolution, were thought to lack Hox and ParaHox genes, suggesting that these critical axial-patterning genes arose after sponges diverged. However, a recent study has challenged this long-held doctrine by claiming to identify ParaHox genes (Cdx family) in two calcareous sponge species, Sycon ciliatum and Leucosolenia complicata. We reanalyzed the main data sets in this paper and analyzed an additional data set that expanded the number of bilaterians represented and removed outgroup homeodomains. As in the previous study, our Neighbor-Joining analyses of the original data sets recovered a clade that included sponge and Cdx genes, whereas Bayesian analyses placed these sponge genes within the NKL subclass of homeodomains. Unlike the original study, only one of our two maximum-likelihood analyses was congruent with Cdx genes in sponges. Our analyses of our additional data set led to the sponge genes consistently being placed within the NKL subclass of homeodomains regardless of method or model. Our results show more support for these sponge genes belonging to the NKL subclass, and therefore imply that Hox and ParaHox genes arose after Porifera diverged from the rest of animals.

Variable levels of drift in tunicate cardiopharyngeal gene regulatory elements.

BACKGROUND: Mutations in gene regulatory networks often lead to genetic divergence without impacting gene expression or developmental patterning. The rules governing this process of developmental systems drift, including the variable impact of selective constraints on different nodes in a gene regulatory network, remain poorly delineated. RESULTS: Here we examine developmental systems drift within the cardiopharyngeal gene regulatory networks of two tunicate species, Corella inflata and Ciona robusta. Cross-species analysis of regulatory elements suggests that trans-regulatory architecture is largely conserved between these highly divergent species. In contrast, cis-regulatory elements within this network exhibit distinct levels of conservation. In particular, while most of the regulatory elements we analyzed showed extensive rearrangements of functional binding sites, the enhancer for the cardiopharyngeal transcription factor FoxF is remarkably well-conserved. Even minor alterations in spacing between binding sites lead to loss of FoxF enhancer function, suggesting that bound trans-factors form position-dependent complexes. CONCLUSIONS: Our findings reveal heterogeneous levels of divergence across cardiopharyngeal cis-regulatory elements. These distinct levels of divergence presumably reflect constraints that are not clearly associated with gene function or position within the regulatory network. Thus, levels of cis-regulatory divergence or drift appear to be governed by distinct structural constraints that will be difficult to predict based on network architecture.

The complete mitochondrial genome of the sea urchin, Echinometra sp. EZ.

The complete mitogenome of Echinometra sp. EZ has been described and fully annotated in this study. Phylogenetic analysis of cytochrome c oxidase subunit I (COI) from six Echinometra species confirms that our sample is E. sp. EZ. The mitogenome is 15,698 bp in length and contains 13 protein-coding genes, 22 tRNAs, 2 rRNAs, and a non-coding region with an identical organization to other Echinoidea. The E. sp. EZ mitogenome shared ∼99.1% identity to the published Echinometra mathaei mitogenome, differing by 147 SNPs. The E. sp. EZ mitogenome will serve as a resource that can be applied to disentangling the Echinometra species complex and to future population genetic studies of this ecologically important sea urchin species.

Adaptation to climate change: trade-offs among responses to multiple stressors in an intertidal crustacean.

Trade-offs may influence both physiological and evolutionary responses to co-occurring stressors, but their effects on both plastic and adaptive responses to climate change are poorly understood. To test for genetic and physiological trade-offs incurred in tolerating multiple stressors, we hybridized two populations of the intertidal copepod Tigriopus californicus that were divergent for both heat and salinity tolerance. Starting in the F2 generation, we selected for increased tolerance of heat, low salinity, and high salinity in replicate lines. After five generations of selection, heat-selected lines had greater heat tolerance but lower fecundity, indicating an energetic cost to tolerance. Lines selected for increased salinity tolerance did not show evidence of adaptation to their respective environments; however, hypo-osmotic selection lines showed substantial loss of tolerance to hyperosmotic stress. Neither of the salinity selection regimes resulted in diminished heat tolerance at ambient salinity; however, simultaneous exposure to heat and hypo-osmotic stress led to decreased heat tolerance, implying a physiological trade-off in tolerance to the two stressors. When we quantified the transcriptomic response to heat and salinity stress via RNA sequencing, we observed little overlap in the stress responses, suggesting the observed synergistic effects of heat and salinity stress were driven by competing energetic demands, rather than shared stress response pathways.

Discordance between morphological and molecular species boundaries among Caribbean species of the reef sponge Callyspongia.

Sponges are among the most species-rich and ecologically important taxa on coral reefs, yet documenting their diversity is difficult due to the simplicity and plasticity of their morphological characters. Genetic attempts to identify species are hampered by the slow rate of mitochondrial sequence evolution characteristic of sponges and some other basal metazoans. Here we determine species boundaries of the Caribbean coral reef sponge genus Callyspongia using a multilocus, model-based approach. Based on sequence data from one mitochondrial (COI), one ribosomal (28S), and two single-copy nuclear protein-coding genes, we found evolutionarily distinct lineages were not concordant with current species designations in Callyspongia. While C. fallax,C. tenerrima, and C. plicifera were reciprocally monophyletic, four taxa with different morphologies (C. armigera,C. longissima,C. eschrichtii, and C. vaginalis) formed a monophyletic group and genetic distances among these taxa overlapped distances within them. A model-based method of species delimitation supported collapsing these four into a single evolutionary lineage. Variation in spicule size among these four taxa was partitioned geographically, not by current species designations, indicating that in Callyspongia, these key taxonomic characters are poor indicators of genetic differentiation. Taken together, our results suggest a complex relationship between morphology and species boundaries in sponges.