scRNA-seq analysis of cells comprising the amphioxus notochord.

Cephalochordates occupy a key phylogenetic position for deciphering the origin and evolution of chordates, since they diverged earlier than urochordates and vertebrates. The notochord is the most prominent feature of chordates. The amphioxus notochord features coin-shaped cells bearing myofibrils. Notochord-derived hedgehog signaling contributes to patterning of the dorsal nerve cord, as in vertebrates. However, properties of constituent notochord cells remain unknown at the single-cell level. We examined these properties using Iso-seq analysis, single-cell RNA-seq analysis, and in situ hybridization (ISH). Gene expression profiles broadly categorize notochordal cells into myofibrillar cells and non-myofibrillar cells. Myofibrillar cells occupy most of the central portion of the notochord, and some cells extend the notochordal horn to both sides of the ventral nerve cord. Some notochord myofibrillar genes are not expressed in myotomes, suggesting an occurrence of myofibrillar genes that are preferentially expressed in notochord. On the other hand, non-myofibrillar cells contain dorsal, lateral, and ventral Müller cells, and all three express both hedgehog and Brachyury. This was confirmed by ISH, although expression of hedgehog in ventral Müller cells was minimal. In addition, dorsal Müller cells express neural transmission-related genes, suggesting an interaction with nerve cord. Lateral Müller cells express hedgehog and other signaling-related genes, suggesting an interaction with myotomes positioned lateral to the notochord. Ventral Müller cells also expressed genes for FGF- and EGF-related signaling, which may be associated with development of endoderm, ventral to the notochord. Lateral Müller cells were intermediate between dorsal/ventral Müller cells. Since vertebrate notochord contributes to patterning and differentiation of ectoderm (nerve cord), mesoderm (somite), and endoderm, this investigation provides evidence that an ancestral or original form of vertebrate notochord is present in extant cephalochordates.

Transcriptomic evidence for Brachyury expression in the caudal tip region of adult Ptychodera flava (Hemichordata).

Most metazoans have a single copy of the T-box transcription factor gene Brachyury. This gene is expressed in cells of the blastopore of late blastulae and the archenteron invagination region of gastrulae. It appears to be crucial for gastrulation and mesoderm differentiation of embryos. Although this expression pattern is shared by most deuterostomes, Brachyury expression has not been reported in adult stages. Here we show that Brachyury of an indirect developer, the hemichordate acorn worm Ptychodera flava, is expressed not only in embryonic cells, but also in cells of the caudal tip (anus) region of adults. This spatially restricted expression, shown by whole-mount in situ hybridization, was confirmed by Iso-Seq RNA sequencing and single-cell RNA-seq (scRNA-seq) analysis. Iso-Seq analysis showed that gene expression occurs only in the caudal region of adults, but not in anterior regions, including the stomochord. scRNA-seq analysis showed a cluster that contained Brachyury-expressing cells comprising epidermis- and mesoderm-related cells, but which is unlikely to be associated with the nervous system or muscle. Although further investigation is required to examine the roles of Brachyury in adults, this study provides important clues for extending studies on Brachyury expression involved in development of the most posterior region of deuterostomes.

Expression and possible functions of a horizontally transferred glycosyl hydrolase gene, GH6-1, in Ciona embryogenesis.

BACKGROUND: The Tunicata or Urochordata is the only animal group with the ability to synthesize cellulose directly and cellulose is a component of the tunic that covers the entire tunicate body. The genome of Ciona intestinalis type A contains a cellulose synthase gene, CesA, that it acquired via an ancient, horizontal gene transfer. CesA is expressed in embryonic epidermal cells and functions in cellulose production. Ciona CesA is composed of both a glycosyltransferase domain, GT2, and a glycosyl hydrolase domain, GH6, which shows a mutation at a key position and seems functionless. Interestingly, the Ciona genome contains a glycosyl hydrolase gene, GH6-1, in which the GH6 domain seems intact. This suggests expression and possible functions of GH6-1 during Ciona embryogenesis. Is GH6-1 expressed during embryogenesis? If so, in what tissues is the gene expressed? Does GH6-1 serve a function? If so, what is it? Answers to these questions may advance our understanding of evolution of this unique animal group. RESULTS: Quantitative reverse transcription PCR and in situ hybridization revealed that GH6-1 is expressed in epidermis of tailbud embryos and in early swimming larvae, a pattern similar to that of CesA. Expression is downregulated at later stages and becomes undetectable in metamorphosed juveniles. The GH6-1 expression level is higher in the anterior-trunk region and caudal-tip regions of late embryos. Single-cell RNA sequencing analysis of the late tailbud stage showed that cells of three clusters with epidermal identity express GH6-1, and that some of them co-express CesA. TALEN-mediated genome editing was used to generate GH6-1 knockout Ciona larvae. Around half of TALEN-electroporated larvae showed abnormal development of adhesive papillae and altered distribution of surface cellulose. In addition, three-fourths of TALEN-electroporated animals failed to complete larval metamorphosis. CONCLUSIONS: This study showed that tunicate GH6-1, a gene that originated by horizontal gene transfer of a prokaryote gene, is recruited into the ascidian genome, and that it is expressed and functions in epidermal cells of ascidian embryos. Although further research is required, this observation demonstrates that both CesA and GH6-1 are involved in tunicate cellulose metabolism, impacting tunicate morphology and ecology.

An environmental DNA metabarcoding survey reveals generic-level occurrence of scleractinian corals at reef slopes of Okinawa Island.

Coral reefs have the highest biodiversity of all marine ecosystems in tropical and subtropical oceans. However, scleractinian corals, keystone organisms of reef productivity, are facing a crisis due to climate change and anthropogenic activities. A broad survey of reef-building corals is essential for worldwide reef preservation. To this end, direct observations made by coral-specialist divers might be supported by another robust method. We improved a recently devised environmental DNA (eDNA) metabarcoding method to identify more than 43 scleractinian genera by sampling 2 l of surface seawater above reefs. Together with direct observations by divers, we assessed the utility of eDNA at 63 locations spanning approximately 250 km near Okinawa Island. Slopes of these islands are populated by diverse coral genera, whereas shallow 'moats' sustain fewer and less varied coral taxa. Major genera recorded by divers included Acropora, Pocillopora, Porites and Montipora, the presence of which was confirmed by eDNA analyses. In addition, eDNA identified more genera than direct observations and documented the presence of previously unrecorded species. This scleractinian coral-specific eDNA method promises to be a powerful tool to survey coral reefs broadly, deeply and robustly.

Active Expression of Genes for Protein Modification Enzymes in Habu Venom Glands.

Genes encoding snake venom toxins have been studied extensively. However, genes involved in the modification and functioning of venom proteins are little known. Protobothrops is a genus of pit vipers, which are venomous and inhabit the Nansei (Southwest) islands of Japan, Taiwan China, Vietnam, Thailand, Myanmar, Nepal, Bhutan, and India. Our previous study decoded the genome of Protobothrops flavoviridis, a species endemic to the Nansei Islands, Japan, and revealed unique evolutionary processes of some venom genes. In this study, we analyzed genes that are highly expressed in venom glands to survey genes for candidate enzymes or chaperone proteins involved in toxin folding and modification. We found that, in addition to genes that encode venom proteins and ribosomal proteins, genes that encode protein disulfide isomerase (PDI) family members (orthologs of human P4HB and PDIA3), Selenoprotein M (SELENOM), and Calreticulin (CALR) are highly expressed in venom glands. Since these enzymes or chaperones are involved in protein modification and potentially possess protein folding functions, we propose that P4HB, SELENOM, CALR, and PDIA3 encode candidate enzymes or chaperones to confer toxic functions upon the venom transcriptome.

A single-cell RNA-seq analysis of Brachyury-expressing cell clusters suggests a morphogenesis-associated signal center of oral ectoderm in sea urchin embryos.

Brachyury is a T-box family transcription factor and plays pivotal roles in morphogenesis. In sea urchin embryos, Brachyury is expressed in the invaginating endoderm, and in the oral ectoderm of the invaginating mouth opening. The oral ectoderm is hypothesized to serve as a signaling center for oral (ventral)-aboral (dorsal) axis formation and to function as a ventral organizer. Our previous results of a single-cell RNA-seq (scRNA-seq) atlas of early Strongylocentrotus purpuratus embryos categorized the constituent cells into 22 clusters, in which the endoderm consists of three clusters and the oral ectoderm four clusters (Foster et al., 2020). Here we examined which clusters of cells expressed Brachyury in relation to the morphogenesis and the identity of the ventral organizer. Our results showed that cells of all three endoderm clusters expressed Brachyury in blastulae. Based on expression profiles of genes involved in the gene regulatory networks (GRNs) of sea urchin embryos, the three clusters are distinguishable, two likely derived from the Veg2 tier and one from the Veg1 tier. On the other hand, of the four oral-ectoderm clusters, cells of two clusters expressed Brachyury at the gastrula stage and genes that are responsible for the ventral organizer at the late blastula stage, but the other two clusters did not. At a single-cell level, most cells of the two oral-ectoderm clusters expressed organizer-related genes, nearly a half of which coincidently expressed Brachyury. This suggests that the ventral organizer contains Brachyury-positive cells which invaginate to form the stomodeum. This scRNA-seq study therefore highlights significant roles of Brachyury-expressing cells in body-plan formation of early sea urchin embryos, though cellular and molecular mechanisms for how Brachyury functions in these processes remain to be elucidated in future studies.

Genome and transcriptome assemblies of the kuruma shrimp, Marsupenaeus japonicus.

The kuruma shrimp Marsupenaeus japonicus (order Decapoda, family Penaeidae) is an economically important crustacean that occurs in shallow, warm seas across the Indo-Pacific. Here, using a combination of Illumina and Oxford Nanopore Technologies platforms, we produced a draft genome assembly of M. japonicus (1.70 Gbp; 18,210 scaffolds; scaffold N50 = 234.9 kbp; 34.38% GC, 93.4% BUSCO completeness) and a complete mitochondrial genome sequence (15,969 bp). As with other penaeid shrimp genomes, the M. japonicus genome is extremely rich in simple repeats, which occupies 27.4% of the assembly. A total of 26,381 protein-coding gene models (94.7% BUSCO completeness) were predicted, of which 18,005 genes (68.2%) were assigned functional description by at least one method. We also produced an Illumina-based transcriptome shotgun assembly (40,991 entries; 93.0% BUSCO completeness) and a PacBio Iso-Seq transcriptome assembly (25,415 entries; 67.5% BUSCO completeness). We envision that the M. japonicus genome and transcriptome assemblies will serve as useful resources for the basic research, fisheries management, and breeding programs of M. japonicus.

A Preliminary Single-Cell RNA-Seq Analysis of Embryonic Cells That Express Brachyury in the Amphioxus, Branchiostoma japonicum.

Among chordate taxa, the cephalochordates diverged earlier than urochordates and vertebrates; thus, they retain unique, primitive developmental features. In particular, the amphioxus notochord has muscle-like properties, a feature not seen in urochordates or vertebrates. Amphioxus contains two Brachyury genes, Bra1 and Bra2. Bra2 is reportedly expressed in the blastopore, notochord, somites, and tail bud, in contrast to a low level of Bra1 expression only in notochord. To distinguish the expression profiles of the two Brachyury genes at the single-cell level, we carried out single-cell RNA-seq (scRNA-seq) analysis using the amphioxus, Branchiostoma japonicum. This scRNA-seq analysis classified B. japonicum embryonic cells into 15 clusters at developmental stages from midgastrula to early swimming larva. Brachyury was expressed in cells of clusters 4, 5, 8, and 9. We first confirmed that cluster 8 comprises cells that form somites since this cluster specifically expresses four myogenic factor genes. Cluster 9 contains a larger number of cells with high levels of Bra2 expression and a smaller number of cells with Bra1 expression. Simultaneous expression in cluster 9 of tool-kit genes, including FoxA, Goosecoid, and hedgehog, showed that this cluster comprises cells that form the notochord. Expression of Bra2, but not Bra1, in cells of clusters 4 and 5 at the gastrula stage together with expression of Wnt1 and Caudal indicates that clusters 4 and 5 comprise cells of the blastopore, which contiguously form the tail bud. In addition, Hox1, Hox3, and Hox4 were highly expressed in Bra2-expressing clusters 4, 5, 8, and 9 in a temporally coordinated manner, suggesting roles of anterior Hox genes in specification of mesodermal organs, including somites, notochord, and tail bud. This scRNA-seq analysis therefore highlights differences between the two Brachyury genes in relation to embryonic regions in which they are expressed and their levels of expression. Bra2 is the ancestral Brachyury in amphioxus, since expression in the blastopore is shared with other deuterostomes. On the other hand, Bra1 is a duplicate copy and likely evolved a supplementary function in notochord and somite formation in the Branchiostoma lineage.

Expansion and Diversification of Fluorescent Protein Genes in Fifteen Acropora Species during the Evolution of Acroporid Corals.

In addition to a purple, non-fluorescent chromoprotein (ChrP), fluorescent proteins (FPs) account for the vivid colors of corals, which occur in green (GFP), cyan (CFP), and red (RFP) FPs. To understand the evolution of the coral FP gene family, we examined the genomes of 15 Acropora species and three confamilial taxa. This genome-wide survey identified 219 FP genes. Molecular phylogeny revealed that the 15 Acropora species each have 9-18 FP genes, whereas the other acroporids examined have only two, suggesting a pronounced expansion of the FP genes in the genus Acropora. The data estimates of FP gene duplication suggest that the last common ancestor of the Acropora species that survived in the period of high sea surface temperature (Paleogene period) has already gained 16 FP genes. Different evolutionary histories of lineage-specific duplication and loss were discovered among GFP/CFPs, RFPs, and ChrPs. Synteny analysis revealed core GFP/CFP, RFP, and ChrP gene clusters, in which a tandem duplication of the FP genes was evident. The expansion and diversification of Acropora FPs may have contributed to the present-day richness of this genus.

Chromosomal Inversion Polymorphisms in Two Sympatric Ascidian Lineages.

Chromosomal rearrangements can reduce fitness of heterozygotes and can thereby prevent gene flow. Therefore, such rearrangements can play a role in local adaptation and speciation. In particular, inversions are considered to be a major potential cause for chromosomal speciation. There are two closely related, partially sympatric lineages of ascidians in the genus Ciona, which we call type-A and type-B animals in the present study. Although these invertebrate chordates are largely isolated reproductively, hybrids can be found in wild populations, suggesting incomplete prezygotic barriers. Although the genome of type-A animals has been decoded and widely used, the genome for type-B animals has not been decoded at the chromosomal level. In the present study, we sequenced the genomes of two type-B individuals from different sides of the English Channel (in the zone of sympatry with type-A individuals) and compared them at the chromosomal level with the type-A genome. Although the overall structures were well conserved between type A and type B, chromosomal alignments revealed many inversions differentiating these two types of Ciona; it is probable that the frequent inversions have contributed to separation between these two lineages. In addition, comparisons of the genomes between the two type-B individuals revealed that type B had high rates of inversion polymorphisms and nucleotide polymorphisms, and thus type B might be in the process of differentiation into multiple new types or species. Our results suggest an important role of inversions in chromosomal speciation of these broadcasting spawners.

Color morphs of the coral, Acropora tenuis, show different responses to environmental stress and different expression profiles of fluorescent-protein genes.

Corals of the family Acroporidae are key structural components of reefs that support the most diverse marine ecosystems. Due to increasing anthropogenic stresses, coral reefs are in decline. Along the coast of Okinawa, Japan, three different color morphs of Acropora tenuis have been recognized for decades. These include brown (N morph), yellow green (G), and purple (P) forms. The tips of axial polyps of each morph exhibit specific fluorescence spectra. This attribute is inherited asexually, and color morphs do not change seasonally. In Okinawa Prefecture, during the summer of 2017, N and P morphs experienced bleaching, in which many N morphs died. Dinoflagellates (Symbiodiniaceae) are essential partners of scleractinian corals, and photosynthetic activity of symbionts was reduced in N and P morphs. In contrast, G morphs successfully withstood the stress. Examination of the clade and type of Symbiodiniaceae indicated that the three color-morphs host similar sets of Clade-C symbionts, suggesting that beaching of N and P morphs is unlikely attributable to differences in the clade of Symbiodiniaceae the color morphs hosted. Fluorescent proteins play pivotal roles in physiological regulation of corals. Since the A. tenuis genome has been decoded, we identified five genes for green fluorescent proteins (GFPs), two for cyan fluorescent proteins (CFPs), three for red fluorescent proteins (RFPs), and seven genes for chromoprotein (ChrP). A summer survey of gene expression profiles under outdoor aquarium conditions demonstrated that (a) expression of CFP and REP was quite low during the summer in all three morphs, (b) P morphs expressed higher levels of ChrP than N and G morphs, (c) both N and G morphs expressed GFP more highly than P morphs, and (d) GFP expression in N morphs was reduced during summer whereas G morphs maintained high levels of GFP expression throughout the summer. Although further studies are required to understand the biological significance of these color morphs of A. tenuis, our results suggest that thermal stress resistance is modified by genetic mechanisms that coincidentally lead to diversification of color morphs of this coral.

A New Dinoflagellate Genome Illuminates a Conserved Gene Cluster Involved in Sunscreen Biosynthesis.

Photosynthetic dinoflagellates of the Family Symbiodiniaceae live symbiotically with many organisms that inhabit coral reefs and are currently classified into fifteen groups, including seven genera. Draft genomes from four genera, Symbiodinium, Breviolum, Fugacium, and Cladocopium, which have been isolated from corals, have been reported. However, no genome is available from the genus Durusdinium, which occupies an intermediate phylogenetic position in the Family Symbiodiniaceae and is well known for thermal tolerance (resistance to bleaching). We sequenced, assembled, and annotated the genome of Durusdinium trenchii, isolated from the coral, Favia speciosa, in Okinawa, Japan. Assembled short reads amounted to 670 Mb with ∼47% GC content. This GC content was intermediate among taxa belonging to the Symbiodiniaceae. Approximately 30,000 protein-coding genes were predicted in the D. trenchii genome, fewer than in other genomes from the Symbiodiniaceae. However, annotations revealed that the D. trenchii genome encodes a cluster of genes for synthesis of mycosporine-like amino acids, which absorb UV radiation. Interestingly, a neighboring gene in the cluster encodes a glucose-methanol-choline oxidoreductase with a flavin adenine dinucleotide domain that is also found in Symbiodinium tridacnidorum. This conservation seems to partially clarify an ancestral genomic structure in the Symbiodiniaceae and its loss in late-branching lineages, including Breviolum and Cladocopium, after splitting from the Durusdinium lineage. Our analysis suggests that approximately half of the taxa in the Symbiodiniaceae may maintain the ability to synthesize mycosporine-like amino acids. Thus, this work provides a significant genomic resource for understanding the genomic diversity of Symbiodiniaceae in corals.

Integrated omics unveil the secondary metabolic landscape of a basal dinoflagellate.

BACKGROUND: Some dinoflagellates cause harmful algal blooms, releasing toxic secondary metabolites, to the detriment of marine ecosystems and human health. Our understanding of dinoflagellate toxin biosynthesis has been hampered by their unusually large genomes. To overcome this challenge, for the first time, we sequenced the genome, microRNAs, and mRNA isoforms of a basal dinoflagellate, Amphidinium gibbosum, and employed an integrated omics approach to understand its secondary metabolite biosynthesis. RESULTS: We assembled the ~ 6.4-Gb A. gibbosum genome, and by probing decoded dinoflagellate genomes and transcriptomes, we identified the non-ribosomal peptide synthetase adenylation domain as essential for generation of specialized metabolites. Upon starving the cells of phosphate and nitrogen, we observed pronounced shifts in metabolite biosynthesis, suggestive of post-transcriptional regulation by microRNAs. Using Iso-Seq and RNA-seq data, we found that alternative splicing and polycistronic expression generate different transcripts for secondary metabolism. CONCLUSIONS: Our genomic findings suggest intricate integration of various metabolic enzymes that function iteratively to synthesize metabolites, providing mechanistic insights into how dinoflagellates synthesize secondary metabolites, depending upon nutrient availability. This study provides insights into toxin production associated with dinoflagellate blooms. The genome of this basal dinoflagellate provides important clues about dinoflagellate evolution and overcomes the large genome size, which has been a challenge previously.

Comparative genomics of four strains of the edible brown alga, Cladosiphon okamuranus.

BACKGROUND: The brown alga, Cladosiphon okamuranus (Okinawa mozuku), is one of the most important edible seaweeds, and it is cultivated for market primarily in Okinawa, Japan. Four strains, denominated S, K, O, and C, with distinctively different morphologies, have been cultivated commercially since the early 2000s. We previously reported a draft genome of the S-strain. To facilitate studies of seaweed biology for future aquaculture, we here decoded and analyzed genomes of the other three strains (K, O, and C). RESULTS: Here we improved the genome of the S-strain (ver. 2, 130 Mbp, 12,999 genes), and decoded the K-strain (135 Mbp, 12,511 genes), the O-strain (140 Mbp, 12,548 genes), and the C-strain (143 Mbp, 12,182 genes). Molecular phylogenies, using mitochondrial and nuclear genes, showed that the S-strain diverged first, followed by the K-strain, and most recently the C- and O-strains. Comparisons of genome architecture among the four strains document the frequent occurrence of inversions. In addition to gene acquisitions and losses, the S-, K-, O-, and C-strains possess 457, 344, 367, and 262 gene families unique to each strain, respectively. Comprehensive Blast searches showed that most genes have no sequence similarity to any entries in the non-redundant protein sequence database, although GO annotation suggested that they likely function in relation to molecular and biological processes and cellular components. CONCLUSIONS: Our study compares the genomes of four strains of C. okamuranus and examines their phylogenetic relationships. Due to global environmental changes, including temperature increases, acidification, and pollution, brown algal aquaculture is facing critical challenges. Genomic and phylogenetic information reported by the present research provides useful tools for isolation of novel strains.

An Investigation into the Genetic History of Japanese Populations of Three Starfish, Acanthaster planci, Linckia laevigata, and Asterias amurensis, Based on Complete Mitochondrial DNA Sequences.

Crown-of-thorns starfish, Acanthaster planci (COTS), are common in coral reefs of Indo-Pacific Ocean. Since they are highly fecund predators of corals, periodic outbreaks of COTS cause substantial loss of healthy coral reefs. Using complete mitochondrial DNA sequences, we here examined how COTS outbreaks in the Ryukyu Archipelago, Japan are reflected by the profile of their population genetics. Population genetics of the blue starfish, Linckia laevigata, which lives in the Ryukyu Archipelago, but not break out and the northern Pacific sea star, Asterias amurensis, which lives in colder seawater around the main Islands of Japan, were also examined as controls. Our results showed that As. amurensis has at least two local populations that diverged approximately 4.7 million years ago (MYA), and no genetic exchanges have occurred between the populations since then. Linckia laevigata shows two major populations in the Ryukyu Archipelago that likely diverged ∼6.8 MYA. The two populations, each comprised of individuals collected from coast of the Okinawa Island and those from the Ishigaki Island, suggest the presence of two cryptic species in the Ryukyu Archipelago. On the other hand, population genetics of COTS showed a profile quite different from those of Asterias and Linckia At least five lineages of COTS have arisen since their divergence ∼0.7 MYA, and each of the lineages is present at the Okinawa Island, Miyako Island, and Ishigaki Island. These results suggest that COTS have experienced repeated genetic bottlenecks that may be associated with or caused by repeated outbreaks.

Alternative mRNA Splicing in Three Venom Families Underlying a Possible Production of Divergent Venom Proteins of the Habu Snake, Protobothrops flavoviridis.

Snake venoms are complex mixtures of toxic proteins encoded by various gene families that function synergistically to incapacitate prey. A huge repertoire of snake venom genes and proteins have been reported, and alternative splicing is suggested to be involved in the production of divergent gene transcripts. However, a genome-wide survey of the transcript repertoire and the extent of alternative splicing still remains to be determined. In this study, the comprehensive analysis of transcriptomes in the venom gland was achieved by using PacBio sequencing. Extensive alternative splicing was observed in three venom protein gene families, metalloproteinase (MP), serine protease (SP), and vascular endothelial growth factors (VEGF). Eleven MP and SP genes and a VEGF gene are expressed as a total of 81, 61, and 8 transcript variants, respectively. In the MP gene family, individual genes are transcribed into different classes of MPs by alternative splicing. We also observed trans-splicing among the clustered SP genes. No other venom genes as well as non-venom counterpart genes exhibited alternative splicing. Our results thus indicate a potential contribution of mRNA alternative and trans-splicing in the production of highly variable transcripts of venom genes in the habu snake.

A Nearly Complete Genome of Ciona intestinalis Type A (C. robusta) Reveals the Contribution of Inversion to Chromosomal Evolution in the Genus Ciona.

Since its initial publication in 2002, the genome of Ciona intestinalis type A (Ciona robusta), the first genome sequence of an invertebrate chordate, has provided a valuable resource for a wide range of biological studies, including developmental biology, evolutionary biology, and neuroscience. The genome assembly was updated in 2008, and it included 68% of the sequence information in 14 pairs of chromosomes. However, a more contiguous genome is required for analyses of higher order genomic structure and of chromosomal evolution. Here, we provide a new genome assembly for an inbred line of this animal, constructed with short and long sequencing reads and Hi-C data. In this latest assembly, over 95% of the 123 Mb of sequence data was included in the chromosomes. Short sequencing reads predicted a genome size of 114-120 Mb; therefore, it is likely that the current assembly contains almost the entire genome, although this estimate of genome size was smaller than previous estimates. Remapping of the Hi-C data onto the new assembly revealed a large inversion in the genome of the inbred line. Moreover, a comparison of this genome assembly with that of Ciona savignyi, a different species in the same genus, revealed many chromosomal inversions between these two Ciona species, suggesting that such inversions have occurred frequently and have contributed to chromosomal evolution of Ciona species. Thus, the present assembly greatly improves an essential resource for genome-wide studies of ascidians.

A draft nuclear-genome assembly of the acoel flatworm Praesagittifera naikaiensis.

BACKGROUND: Acoels are primitive bilaterians with very simple soft bodies, in which many organs, including the gut, are not developed. They provide platforms for studying molecular and developmental mechanisms involved in the formation of the basic bilaterian body plan, whole-body regeneration, and symbiosis with photosynthetic microalgae. Because genomic information is essential for future research on acoel biology, we sequenced and assembled the nuclear genome of an acoel, Praesagittifera naikaiensis. FINDINGS: To avoid sequence contamination derived from symbiotic microalgae, DNA was extracted from embryos that were free of algae. More than 290x sequencing coverage was achieved using a combination of Illumina (paired-end and mate-pair libraries) and PacBio sequencing. RNA sequencing and Iso-Seq data from embryos, larvae, and adults were also obtained. First, a preliminary ∼17-kilobase pair (kb) mitochondrial genome was assembled, which was deleted from the nuclear sequence assembly. As a result, a draft nuclear genome assembly was ∼656 Mb in length, with a scaffold N50 of 117 kb and a contig N50 of 57 kb. Although ∼70% of the assembled sequences were likely composed of repetitive sequences that include DNA transposons and retrotransposons, the draft genome was estimated to contain 22,143 protein-coding genes, ∼99% of which were substantiated by corresponding transcripts. We could not find horizontally transferred microalgal genes in the acoel genome. Benchmarking Universal Single-Copy Orthologs analyses indicated that 77% of the conserved single-copy genes were complete. Pfam domain analyses provided a basic set of gene families for transcription factors and signaling molecules. CONCLUSIONS: Our present sequencing and assembly of the P. naikaiensis nuclear genome are comparable to those of other metazoan genomes, providing basic information for future studies of genic and genomic attributes of this animal group. Such studies may shed light on the origins and evolution of simple bilaterians.

A siphonous macroalgal genome suggests convergent functions of homeobox genes in algae and land plants.

Genome evolution and development of unicellular, multinucleate macroalgae (siphonous algae) are poorly known, although various multicellular organisms have been studied extensively. To understand macroalgal developmental evolution, we assembled the ∼26 Mb genome of a siphonous green alga, Caulerpa lentillifera, with high contiguity, containing 9,311 protein-coding genes. Molecular phylogeny using 107 nuclear genes indicates that the diversification of the class Ulvophyceae, including C. lentillifera, occurred before the split of the Chlorophyceae and Trebouxiophyceae. Compared with other green algae, the TALE superclass of homeobox genes, which expanded in land plants, shows a series of lineage-specific duplications in this siphonous macroalga. Plant hormone signalling components were also expanded in a lineage-specific manner. Expanded transport regulators, which show spatially different expression, suggest that the structural patterning strategy of a multinucleate cell depends on diversification of nuclear pore proteins. These results not only imply functional convergence of duplicated genes among green plants, but also provide insight into evolutionary roots of green plants. Based on the present results, we propose cellular and molecular mechanisms involved in the structural differentiation in the siphonous alga.

Diversified secondary metabolite biosynthesis gene repertoire revealed in symbiotic dinoflagellates.

Symbiodiniaceae dinoflagellates possess smaller nuclear genomes than other dinoflagellates and produce structurally specialized, biologically active, secondary metabolites. Till date, little is known about the evolution of secondary metabolism in dinoflagellates as comparative genomic approaches have been hampered by their large genome sizes. Here, we overcome this challenge by combining genomic and metabolomics approaches to investigate how chemical diversity arises in three decoded Symbiodiniaceae genomes (clades A3, B1 and C). Our analyses identify extensive diversification of polyketide synthase and non-ribosomal peptide synthetase genes from two newly decoded genomes of Symbiodinium tridacnidorum (A3) and Cladocopium sp. (C). Phylogenetic analyses indicate that almost all the gene families are derived from lineage-specific gene duplications in all three clades, suggesting divergence for environmental adaptation. Few metabolic pathways are conserved among the three clades and we detect metabolic similarity only in the recently diverged clades, B1 and C. We establish that secondary metabolism protein architecture guides substrate specificity and that gene duplication and domain shuffling have resulted in diversification of secondary metabolism genes.