Speciation of pelagic zooplankton: Invisible boundaries can drive isolation of oceanic ctenophores.

The study of evolution and speciation in non-model systems provides us with an opportunity to expand our understanding of biodiversity in nature. Connectivity studies generally focus on species with obvious boundaries to gene flow, but in open-ocean environments, such boundaries are difficult to identify. Due to the lack of obvious boundaries, speciation and population subdivision in the pelagic environment remain largely unexplained. Comb jellies (Phylum Ctenophora) are mostly planktonic gelatinous invertebrates, many of which are considered to have freely interbreeding distributions worldwide. It is thought that the lobate ctenophore Bolinopsis infundibulum is distributed throughout cooler northern latitudes and B. vitrea warmer. Here, we examined the global population structure for species of Bolinopsis with genetic and morphological data. We found distinct evolutionary patterns within the genus, where B. infundibulum had a broad distribution from northern Pacific to Atlantic waters despite many physical barriers, while other species were geographically segregated despite few barriers. Divergent patterns of speciation within the genus suggest that oceanic currents, sea-level, and geological changes over time can act as either barriers or aids to dispersal in the pelagic environment. Further, we used population genomic data to examine evolution in the open ocean of a distinct lineage of Bolinopsis ctenophores from the North Eastern Pacific. Genetic information and morphological observations validated this as a separate species, Bolinopsis microptera, which was previously described but has recently been called B. infundibulum. We found that populations of B. microptera from California were in cytonuclear discordance, which indicates a secondary contact zone for previously isolated populations. Discordance at this scale is rare, especially in a continuous setting.

Atolla reynoldsi sp. nov. (Cnidaria, Scyphozoa, Coronatae, Atollidae): A New Species of Coronate Scyphozoan Found in the Eastern North Pacific Ocean.

We have observed and collected unusual specimens of what we recognize as undescribed types of the genus Atolla over the past 15 years. Of these, there appear to be three potentially different types. One of these has now been genetically sequenced and compared both morphologically and molecularly with five other Atolla species that have been found in the eastern Pacific. This new variant is so morphologically distinct from other previously described Atolla species that we believe it can be described as a new species, Atolla reynoldsi sp. nov. This species along with two additional types may comprise a new genus. It is also clear that a more accurate and diagnostic morphological key for the genus Atolla needs to be developed. This paper will also provide some potential starting points for a new key to the genus.

Hidden diversity of Ctenophora revealed by new mitochondrial COI primers and sequences.

The mitochondrial gene cytochrome-c-oxidase subunit 1 (COI) is useful in many taxa for phylogenetics, population genetics, metabarcoding, and rapid species identifications. However, the phylum Ctenophora (comb jellies) has historically been difficult to study due to divergent mitochondrial sequences and the corresponding inability to amplify COI with degenerate and standard COI "barcoding" primers. As a result, there are very few COI sequences available for ctenophores, despite over 200 described species in the phylum. Here, we designed new primers and amplified the COI fragment from members of all major groups of ctenophores, including many undescribed species. Phylogenetic analyses of the resulting COI sequences revealed high diversity within many groups that was not evident from more conserved 18S rDNA sequences, in particular among the Lobata (Ctenophora; Tentaculata; Lobata). The COI phylogenetic results also revealed unexpected community structure within the genus Bolinopsis, suggested new species within the genus Bathocyroe, and supported the ecological and morphological differences of some species such as Lampocteis cruentiventer and similar undescribed lobates (Lampocteis sp. "V" stratified by depth, and "A" differentiated by colour). The newly designed primers reported herein provide important tools to enable researchers to illuminate the diversity of ctenophores worldwide via quick molecular identifications, improve the ability to analyse environmental DNA by improving reference libraries and amplifications, and enable a new breadth of population genetic studies.

A chromosome-scale genome assembly and karyotype of the ctenophore Hormiphora californensis.

Here, we present a karyotype, a chromosome-scale genome assembly, and a genome annotation from the ctenophore Hormiphora californensis (Ctenophora: Cydippida: Pleurobrachiidae). The assembly spans 110 Mb in 44 scaffolds and 99.47% of the bases are contained in 13 scaffolds. Chromosome micrographs and Hi-C heatmaps support a karyotype of 13 diploid chromosomes. Hi-C data reveal three large heterozygous inversions on chromosome 1, and one heterozygous inversion shares the same gene order found in the genome of the ctenophore Pleurobrachia bachei. We find evidence that H. californensis and P. bachei share thirteen homologous chromosomes, and the same karyotype of 1n = 13. The manually curated PacBio Iso-Seq-based genome annotation reveals complex gene structures, including nested genes and trans-spliced leader sequences. This chromosome-scale assembly is a useful resource for ctenophore biology and will aid future studies of metazoan evolution and phylogenetics.

Conserved novel ORFs in the mitochondrial genome of the ctenophore Beroe forskalii.

To date, five ctenophore species' mitochondrial genomes have been sequenced, and each contains open reading frames (ORFs) that if translated have no identifiable orthologs. ORFs with no identifiable orthologs are called unidentified reading frames (URFs). If truly protein-coding, ctenophore mitochondrial URFs represent a little understood path in early-diverging metazoan mitochondrial evolution and metabolism. We sequenced and annotated the mitochondrial genomes of three individuals of the beroid ctenophore Beroe forskalii and found that in addition to sharing the same canonical mitochondrial genes as other ctenophores, the B. forskalii mitochondrial genome contains two URFs. These URFs are conserved among the three individuals but not found in other sequenced species. We developed computational tools called pauvre and cuttlery to determine the likelihood that URFs are protein coding. There is evidence that the two URFs are under negative selection, and a novel Bayesian hypothesis test of trinucleotide frequency shows that the URFs are more similar to known coding genes than noncoding intergenic sequence. Protein structure and function prediction of all ctenophore URFs suggests that they all code for transmembrane transport proteins. These findings, along with the presence of URFs in other sequenced ctenophore mitochondrial genomes, suggest that ctenophores may have uncharacterized transmembrane proteins present in their mitochondria.

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.

Symplectin evolved from multiple duplications in bioluminescent squid.

The squid Sthenoteuthis oualaniensis, formerly Symplectoteuthis oualaniensis, generates light using the luciferin coelenterazine and a unique enzyme, symplectin. Genetic information is limited for bioluminescent cephalopod species, so many proteins, including symplectin, occur in public databases only as sequence isolates with few identifiable homologs. As the distribution of the symplectin/pantetheinase protein family in Metazoa remains mostly unexplored, we have sequenced the transcriptomes of four additional luminous squid, and make use of publicly available but unanalyzed data of other cephalopods, to examine the occurrence and evolution of this protein family. While the majority of spiralians have one or two copies of this protein family, four well-supported groups of proteins are found in cephalopods, one of which corresponds to symplectin. A cysteine that is critical for symplectin functioning is conserved across essentially all members of the protein family, even those unlikely to be used for bioluminescence. Conversely, active site residues involved in pantetheinase catalysis are also conserved across essentially all of these proteins, suggesting that symplectin may have multiple functions including hydrolase activity, and that the evolution of the luminous phenotype required other changes in the protein outside of the main binding pocket.

Non-excitable fluorescent protein orthologs found in ctenophores.

BACKGROUND: Fluorescent proteins are optically active proteins found across many clades in metazoans. A fluorescent protein was recently identified in a ctenophore, but this has been suggested to derive from a cnidarian, raising again the question of origins of this group of proteins. RESULTS: Through analysis of transcriptome data from 30 ctenophores, we identified a member of an orthologous group of proteins similar to fluorescent proteins in each of them, as well as in the genome of Mnemiopsis leidyi. These orthologs lack canonical residues involved in chromophore formation, suggesting another function. CONCLUSIONS: The phylogenetic position of the ctenophore protein family among fluorescent proteins suggests that this gene was present in the common ancestor of all ctenophores and that the fluorescent protein previously found in a ctenophore actually derives from a siphonophore.

Occurrence of Isopenicillin-N-Synthase Homologs in Bioluminescent Ctenophores and Implications for Coelenterazine Biosynthesis.

The biosynthesis of the luciferin coelenterazine has remained a mystery for decades. While not all organisms that use coelenterazine appear to make it themselves, it is thought that ctenophores are a likely producer. Here we analyze the transcriptome data of 24 species of ctenophores, two of which have published genomes. The natural precursors of coelenterazine have been shown to be the amino acids L-tyrosine and L-phenylalanine, with the most likely biosynthetic pathway involving cyclization and further modification of the tripeptide Phe-Tyr-Tyr ("FYY"). Therefore, we searched the ctenophore transcriptome data for genes with the short peptide "FYY" as part of their coding sequence. We recovered a group of candidate genes for coelenterazine biosynthesis in the luminous species which encode a set of highly conserved non-heme iron oxidases similar to isopenicillin-N-synthase. These genes were absent in the transcriptomes and genome of the two non-luminous species. Pairwise identities and substitution rates reveal an unusually high degree of identity even between the most unrelated species. Additionally, two related groups of non-heme iron oxidases were found across all ctenophores, including those which are non-luminous, arguing against the involvement of these two gene groups in luminescence. Important residues for iron-binding are conserved across all proteins in the three groups, suggesting this function is still present. Given the known functions of other members of this protein superfamily are involved in heterocycle formation, we consider these genes to be top candidates for laboratory characterization or gene knockouts in the investigation of coelenterazine biosynthesis.

A comparison across non-model animals suggests an optimal sequencing depth for de novo transcriptome assembly.

BACKGROUND: The lack of genomic resources can present challenges for studies of non-model organisms. Transcriptome sequencing offers an attractive method to gather information about genes and gene expression without the need for a reference genome. However, it is unclear what sequencing depth is adequate to assemble the transcriptome de novo for these purposes. RESULTS: We assembled transcriptomes of animals from six different phyla (Annelids, Arthropods, Chordates, Cnidarians, Ctenophores, and Molluscs) at regular increments of reads using Velvet/Oases and Trinity to determine how read count affects the assembly. This included an assembly of mouse heart reads because we could compare those against the reference genome that is available. We found qualitative differences in the assemblies of whole-animals versus tissues. With increasing reads, whole-animal assemblies show rapid increase of transcripts and discovery of conserved genes, while single-tissue assemblies show a slower discovery of conserved genes though the assembled transcripts were often longer. A deeper examination of the mouse assemblies shows that with more reads, assembly errors become more frequent but such errors can be mitigated with more stringent assembly parameters. CONCLUSIONS: These assembly trends suggest that representative assemblies are generated with as few as 20 million reads for tissue samples and 30 million reads for whole-animals for RNA-level coverage. These depths provide a good balance between coverage and noise. Beyond 60 million reads, the discovery of new genes is low and sequencing errors of highly-expressed genes are likely to accumulate. Finally, siphonophores (polymorphic Cnidarians) are an exception and possibly require alternate assembly strategies.

A photoactivatable green-fluorescent protein from the phylum Ctenophora.

Genes for the family of green-fluorescent proteins (GFPs) have been found in more than 100 species of animals, with some species containing six or more copies producing a variety of colours. Thus far, however, these species have all been within three phyla: Cnidaria, Arthropoda and Chordata. We have discovered GFP-type fluorescent proteins in the phylum Ctenophora, the comb jellies. The ctenophore proteins share the xYG chromophore motif of all other characterized GFP-type proteins. These proteins exhibit the uncommon property of reversible photoactivation, in which fluorescent emission becomes brighter upon exposure to light, then gradually decays to a non-fluorescent state. In addition to providing potentially useful optical probes with novel properties, finding a fluorescent protein in one of the earliest diverging metazoans adds further support to the possibility that these genes are likely to occur throughout animals.

Spectral diversity of fluorescent proteins from the anthozoan Corynactis californica.

Color morphs of the temperate, nonsymbiotic corallimorpharian Corynactis californica show variation in pigment pattern and coloring. We collected seven distinct color morphs of C. californica from subtidal locations in Monterey Bay, California, and found that tissue- and color-morph-specific expression of at least six different genes is responsible for this variation. Each morph contains at least three to four distinct genetic loci that code for these colors, and one morph contains at least five loci. These genes encode a subfamily of new GFP-like proteins, which fluoresce across the visible spectrum from green to red, while sharing between 75% to 89% pairwise amino-acid identity. Biophysical characterization reveals interesting spectral properties, including a bright yellow protein, an orange protein, and a red protein exhibiting a "fluorescent timer" phenotype. Phylogenetic analysis indicates that the FP genes from this species evolved together but that diversification of anthozoan fluorescent proteins has taken place outside of phylogenetic constraints, especially within the Corallimorpharia. The discovery of more examples of fluorescent proteins in a non-bioluminescent, nonsymbiotic anthozoan highlights possibilities of adaptive ecological significance unrelated to light regulation for algal symbionts. The patterns and colors of fluorescent proteins in C. californica and similar species may hold meaning for organisms that possess the visual pigments to distinguish them.

Proteorhodopsin genes are distributed among divergent marine bacterial taxa.

Proteorhodopsin (PR) is a retinal-binding bacterial integral membrane protein that functions as a light-driven proton pump. The gene encoding this photoprotein was originally discovered on a large genome fragment derived from an uncultured marine gamma-proteobacterium of the SAR86 group. Subsequently, many variants of the PR gene have been detected in marine plankton, via PCR-based gene surveys. It has not been clear, however, whether these different PR genes are widely distributed among different bacterial groups, or whether they have a restricted taxonomic distribution. We report here comparative analyses of PR-bearing genomic fragments recovered directly from planktonic bacteria inhabiting the California coast, the central Pacific Ocean, and waters offshore the Antarctica Peninsula. Sequence analysis of an Antarctic genome fragment harboring PR (ANT32C12) revealed moderate conservation in gene order and identity, compared with a previously reported PR-containing genome fragment from a Monterey Bay gamma-proteobacterium (EBAC31A08). Outside the limited region of synteny shared between these clones, however, no significant DNA or protein identity was evident. Analysis of a third PR-containing genome fragment (HOT2C01) from the North Pacific subtropical gyre showed even more divergence from the gamma-proteobacterial PR-flanking region. Subsequent phylogenetic and comparative genomic analyses revealed that the Central North Pacific PR-containing genome fragment (HOT2C01) originated from a planktonic alpha-proteobacterium. These data indicate that PR genes are distributed among a variety of divergent marine bacterial taxa, including both alpha- and gamma-proteobacteria. Our analyses also demonstrate the utility of cultivation-independent comparative genomic approaches for assessing gene content and distribution in naturally occurring microbes.