Reversible Sterilization of Channel Catfish via Overexpression of Glutamic Acid Decarboxylase Gene.

The confinement of transgenic fish is essential to prevent their escape and reproduction in natural ecosystems. Reversible transgenic sterilization is a promising approach to control the reproduction of transgenic fish. Therefore, the present study was conducted to develop a reversibly sterile channel catfish (Ictalurus punctatus) via the transgenic overexpression of the goldfish (Carassius auratus) glutamic acid decarboxylase (GAD) gene driven by the common carp (Cyprinus carpio) ß-actin promoter to disrupt normal gamma-aminobutyric acid (GABA) regulation. Three generations of GAD-transgenic fish were produced. All studied generations showed repressed reproductive performance; however, this was not always statistically significant. In F1, 5.4% of the transgenic fish showed a sexual maturity score ≥ 4 (maximum = 5) at five years of age, which was lower (p = 0.07) than that of the control group (16.8%). In the spawning experiments conducted on F1 transgenic fish at six and nine years of age, 45.5% and 20.0% of fish spawned naturally, representing lower values (p = 0.09 and 0.12, respectively) than the percentages in the sibling control fish of the same age (83.3% and 66.7%, respectively). Four of six pairs of the putative infertile six-year-old fish spawned successfully after luteinizing hormone-releasing hormone analog (LHRHa) therapy. Similar outcomes were noted in the three-year-old F2 fish, with a lower spawning percentage in transgenic fish (20.0%) than in the control (66.7%). In one-year-old F2-generation transgenic fish, the observed mean serum gonadotropin-releasing hormone (GnRH) levels were 9.23 ± 2.49 and 8.14 ± 2.21 ng/mL for the females and males, respectively. In the control fish, the mean levels of GnRH were 11.04 ± 4.06 and 9.03 ± 2.36 ng/mL for the females and males, respectively, which did not differ significantly from the control (p = 0.15 and 0.27 for females and males, respectively). There was no significant difference in the estradiol levels of the female transgenic and non-transgenic fish in the one- and four-year-old F2-generation fish. The four-year-old F2-generation male transgenic fish exhibited significantly (p < 0.05) lower levels of GnRH and testosterone than the control fish. In conclusion, while overexpressing GAD repressed the reproductive abilities of channel catfish, it did not completely sterilize transgenic fish. The sterilization rate might be improved through selection in future generations.

Identification and Expressional Analysis of Putative PRDI-BF1 and RIZ Homology Domain-Containing Transcription Factors in Mulinia lateralis.

Mollusca represents one of the ancient bilaterian groups with high morphological diversity, while the formation mechanisms of the precursors of all germ cells, primordial germ cells (PGCs), have not yet been clarified in mollusks. PRDI-BF1 and RIZ homology domain-containing proteins (PRDMs) are a group of transcriptional repressors, and PRDM1 (also known as BLIMP1) and PRDM14 have been reported to be essential for the formation of PGCs. In the present study, we performed a genome-wide retrieval in Mulinia lateralis and identified 11 putative PRDMs, all of which possessed an N-terminal PR domain. Expressional profiles revealed that all these prdm genes showed specifically high expression levels in the given stages, implying that all PRDMs played important roles during early development stages. Specifically, Ml-prdm1 was highly expressed at the gastrula stage, the key period when PGCs arise, and was specifically localized in the cytoplasm of two or three cells of blastula, gastrula, or trochophore larvae, matching the typical characteristics of PGCs. These results suggested that Ml-prdm1-positive cells may be PGCs and that Ml-prdm1 could be a candidate marker for tracing the formation of PGCs in M. lateralis. In addition, the expression profiles of Ml-prdm14 hinted that it may not be associated with PGCs of M. lateralis. The present study provides insights into the evolution of the PRDM family in mollusks and offers a better understanding of the formation of PGCs in mollusks.

CRISPR/Cas9-mediated knock-in of masu salmon (Oncorhyncus masou) elongase gene in the melanocortin-4 (mc4r) coding region of channel catfish (Ictalurus punctatus) genome.

Channel catfish, Ictalurus punctatus, have limited ability to synthesize Ω-3 fatty acids. The ccßA-msElovl2 transgene containing masu salmon, Oncorhynchus masou, elongase gene driven by the common carp, Cyprinus carpio, ß-actin promoter was inserted into the channel catfish melanocortin-4 receptor (mc4r) gene site using the two-hit two-oligo with plasmid (2H2OP) method. The best performing sgRNA resulted in a knockout mutation rate of 92%, a knock-in rate of 54% and a simultaneous knockout/knock-in rate of 49%. Fish containing both the ccßA-msElovl2 transgene knock-in and mc4r knockout (Elovl2) were 41.8% larger than controls at 6 months post-hatch (p = 0.005). Mean eicosapentaenoic acid (EPA, C20:5n-3) levels in Elov2 mutants and mc4r knockout mutants (MC4R) were 121.6% and 94.1% higher than in controls, respectively (p = 0.045; p = 0.025). Observed mean docosahexaenoic acid (DHA, C22:6n-3) and total EPA + DHA content was 32.8% and 45.1% higher, respectively, in Elovl2 transgenic channel catfish than controls (p = 0.368; p = 0.025). To our knowledge this is the first example of genome engineering to simultaneously target transgenesis and knock-out a gene in a commercially important aquaculture species for multiple improved performance traits. With a high transgene integration rate, improved growth, and higher omega-3 fatty acid content, the use of Elovl2 transgenic channel catfish appears beneficial for application on commercial farms.

Gene Editing of the Follicle-Stimulating Hormone Gene to Sterilize Channel Catfish, Ictalurus punctatus, Using a Modified Transcription Activator-like Effector Nuclease Technology with Electroporation.

Follicle-stimulating hormone (fsh) plays an important role in sexual maturation in catfish. Knocking out the fsh gene in the fish zygote should suppress the reproduction of channel catfish (Ictalurus punctatus). In this study, transcription activator-like effector nuclease (TALEN) plasmids targeting the fsh gene were electroporated into fertilized eggs with the standard double electroporation technique. Targeted fsh cleavage efficiency was 63.2% in P1fsh-knockout catfish. Ten of fifteen (66.7%) control pairs spawned, and their eggs had 32.3-74.3% average hatch rates in 2016 and 2017. Without hormone therapy, the spawning rates of P1 mutants ranged from 33.3 to 40.0%, with an average egg hatching rate of 0.75%. After confirmation of the low fertility of P1 mutants in 2016, human chorionic gonadotropin (HCG) hormone therapy improved the spawning rates by 80% for female mutants and 88.9% for male mutants, and the mean hatch rate was 35.0% for F1 embryos, similar to that of the controls (p > 0.05). Polymerase chain reaction (PCR) identification showed no potential TALEN plasmid integration into the P1 channel catfish genome. Neither the P1 nor the F1 mutant fish showed any noticeable changes in in body weight, survival rate, and hatching rate when the reproductive gene was knocked out. F1 families had a mean inheritance rate of 50.3%. The results brought us one step closer to allowing implementation of certain genetic techniques to aquaculture and fisheries management, while essentially eliminating the potential environment risk posed by transgenic, hybrid, and exotic fish as well as domestic fish.

A New Strategy for Increasing Knock-in Efficiency: Multiple Elongase and Desaturase Transgenes Knock-in by Targeting Long Repeated Sequences.

CRISPR/Cas9-mediated knock-in (KI) has a wide application in gene therapy, gene function study, and transgenic breeding programs. Unlike gene therapy, which requires accurate KI to correct gene mutation, transgenic breeding programs can accept robust KI as long as integration does not interrupt normal gene functions and result in any negative pleiotropic effects. High KI efficiency is required to reduce the breeding cost and shorten the breeding period, especially in transferring multiple foreign genes to a single individual. To elevate the KI efficacy and achieve multiple gene KIs simultaneously, we introduced a new strategy that enables transgene integration into numerous sites of the genome by targeting long repeated sequences (LRSs). Using this simple strategy, for the first time we successfully generated transgenic fish carrying the masu salmon (Oncorhynchus masou) elovl2 gene and rabbitfish (Siganus canaliculatus) Δ4 fad and Δ6 fad genes, and achieved robust target KI of elovl2 and Δ6 fad genes at multiple sites of LRS1 and LRS3, respectively, in the initial generation. This demonstrated that donor plasmid homology arms, which were nearly identical but not completely the same as the genome sequence, still led to on-target KI. Although the target KI efficiencies at LRS1, LRS2, and LRS3 sites were still relatively low in the current study, it is very promising that 100% KI efficiency in the future could be realized and perfected by selection of better LRSs and optimization of sgRNAs.

Echiuran Hox genes provide new insights into the correspondence between Hox subcluster organization and collinearity pattern.

In many bilaterians, Hox genes are generally clustered along the chromosomes and expressed in spatial and temporal order. In vertebrates, the expression of Hox genes follows a whole-cluster spatio-temporal collinearity (WSTC) pattern, whereas in some invertebrates the expression of Hox genes exhibits a subcluster-level spatio-temporal collinearity pattern. In bilaterians, the diversity of collinearity patterns and the cause of collinearity differences in Hox gene expression remain poorly understood. Here, we investigate genomic organization and expression pattern of Hox genes in the echiuran worm Urechis unicinctus (Annelida, Echiura). Urechis unicinctus has a split cluster with four subclusters divided by non-Hox genes: first subcluster (Hox1 and Hox2), second subcluster (Hox3), third subcluster (Hox4, Hox5, Lox5, Antp and Lox4), fourth subcluster (Lox2 and Post2). The expression of U. unicinctus Hox genes shows a subcluster-based whole-cluster spatio-temporal collinearity (S-WSTC) pattern: the anterior-most genes in each subcluster are activated in a spatially and temporally colinear manner (reminiscent of WSTC), with the subsequent genes in each subcluster then being very similar to their respective anterior-most subcluster gene. Combining genomic organization and expression profiles of Hox genes in different invertebrate lineages, we propose that the spatio-temporal collinearity of invertebrate Hox is subcluster-based.

Gene Editing of the Catfish Gonadotropin-Releasing Hormone Gene and Hormone Therapy to Control the Reproduction in Channel Catfish, Ictalurus punctatus.

Transcription activator-like effector nuclease (TALEN) plasmids targeting the channel catfish gonadotropin-releasing hormone (cfGnRH) gene were delivered into fertilized eggs with double electroporation to sterilize channel catfish (Ictalurus punctatus). Targeted cfGnRH fish were sequenced and base deletion, substitution, and insertion were detected. The gene mutagenesis was achieved in 52.9% of P1 fish. P1 mutants (individuals with human-induced sequence changes at the cfGnRH locus) had lower spawning rates (20.0−50.0%) when there was no hormone therapy compared to the control pairs (66.7%) as well as having lower average egg hatch rates (2.0% versus 32.3−74.3%) except for one cfGnRH mutated female that had a 66.0% hatch rate. After low fertility was observed in 2016, application of luteinizing hormone-releasing hormone analog (LHRHa) hormone therapy resulted in good spawning and hatch rates for mutants in 2017, which were not significantly different from the controls (p > 0.05). No exogenous DNA fragments were detected in the genome of mutant P1 fish, indicating no integration of the plasmids. No obvious effects on other economically important traits were observed after the knockout of the reproductive gene in the P1 fish. Growth rates, survival, and appearance between mutant and control individuals were not different. While complete knock-out of reproductive output was not achieved, as these were mosaic P1 brood stock, gene editing of channel catfish for the reproductive confinement of gene-engineered, domestic, and invasive fish to prevent gene flow into the natural environment appears promising.

CRISPR/Cas9-Mediated Transgenesis of the Masu Salmon (Oncorhynchus masou) elovl2 Gene Improves n-3 Fatty Acid Content in Channel Catfish (Ictalurus punctatus).

Omega-3 polyunsaturated fatty acids (n-3 PUFAs), particularly eicosapentaenoic acid (EPA, 20:5n-3) and docosahexaenoic acid (DHA, 22:6n-3), play a very important role in human health. Channel catfish (Ictalurus punctatus) is one of the leading freshwater aquaculture species in the USA, but has low levels of EPA and DHA compared to some fish such as salmon. To improve EPA and DHA content, a modification of the n-3 PUFA biosynthetic pathway was achieved through the insertion of an elovl2 transgene isolated from masu salmon (Oncorhynchus masou) driven by a carp ß-actin promoter using a two-hit by gRNA and two oligos with a targeting plasmid (2H2OP) CRISPR/Cas9 approach. Integration rate of the transgene was high (37.5%) and detected in twelve different tissues of P1 transgenic fish with tissue-specific gene expression. Liver and muscle had relative high gene expression (13.4- and 9.2-fold change, respectively). Fatty acid analysis showed DHA content in the muscle from transgenic fish was 1.62-fold higher than in non-transgenic fish (P < 0.05). Additionally, total n-3 PUFAs and omega-6 polyunsaturated fatty acids (n-6 PUFAs) increased to 1.41-fold and 1.50-fold, respectively, suggesting the ß-actin-elovl2 transgene improved biosynthesis of PUFAs in channel catfish as a whole. The n-9 fatty acid level decreased in the transgenic fish compared to the control. Morphometric analysis showed that there were significant differences between injected fish with sgRNAs (including positive and negative fish) and sham-injected controls (P < 0.001). Potential off-target effects are likely the major factor responsible for morphological deformities. Optimization of sgRNA design to maximize activity and reduce off-target effects of CRISPR/Cas9 should be examined in future transgenic research, but this research shows a promising first step in the improvement of n-3 PUFAs in channel catfish.

Genome-Wide Analyses of Heat Shock Protein Superfamily Provide New Insights on Adaptation to Sulfide-Rich Environments in Urechis unicinctus (Annelida, Echiura).

The intertidal zone is a transitional area of the land-sea continuum, in which physical and chemical properties vary during the tidal cycle and highly toxic sulfides are rich in sediments due to the dynamic regimes. As a typical species thriving in this habitat, Urechis unicinctus presents strong sulfide tolerance and is expected to be a model species for sulfide stress research. Heat shock proteins (HSPs) consist of a large group of highly conserved molecular chaperones, which play important roles in stress responses. In this study, we systematically analyzed the composition and expression of HSPs in U. unicinctus. A total of eighty-six HSP genes from seven families were identified, in which two families, including sHSP and HSP70, showed moderate expansion, and this variation may be related to the benthic habitat of the intertidal zone. Furthermore, expression analysis revealed that almost all the HSP genes in U. unicinctus were significantly induced under sulfide stress, suggesting that they may be involved in sulfide stress response. Weighted gene co-expression network analysis (WGCNA) showed that 12 HSPs, including 5 sHSP and 4 HSP70 family genes, were highly correlated with the sulfide stress response which was distributed in steelblue and green modules. Our data indicate that HSPs, especially sHSP and HSP70 families, may play significant roles in response to sulfide stress in U. unicinctus. This systematic analysis provides valuable information for further understanding of the function of the HSP gene family for sulfide adaptation in U. unicinctus and contributes a better understanding of the species adaptation strategies of marine benthos in the intertidal zone.

Identification and Characterization of MicroRNAs Involving in Initial Sex Differentiation of Chlamys farreri Gonads.

Research on expressional regulation of genes at the initial sex differentiation of gonads will help to elucidate the mechanisms of sex determination and differentiation in animals. However, information on initial sex differentiation of gonads is limited in bivalves. MicroRNAs (miRNAs) are a class of endogenous small noncoding RNAs that can regulate the target gene expression at the posttranscription level by degrading the mRNA or repressing the mRNA translation. In the present study, we investigated the small RNAs transcriptome using the testes and ovaries of Zhikong scallop Chlamys farreri juveniles with a shell height of 5.0 mm, a critical stage of initial sex differentiation of gonads. A total of 75 known mature miRNAs and 103 novel miRNAs were identified. By comparing the expression of miRNAs between the ovary and testis, 11 miRNAs were determined to be differentially expressed. GO annotations and KEGG analyses indicated that many putative target genes that matched to these differentially expressed miRNAs participated in the regulation of sex differentiation. Furthermore, two selected miRNAs, cfa-novel_miR65 and cfa-miR-87a-3p_1, were confirmed to downregulate expressions of Foxl2 (a female-critical gene) and Klf4 (a male-critical gene), respectively, using a dual-luciferase reporter analysis. Our findings provided new insights into the initial sex differentiation of gonads regulated by miRNAs in bivalves.

Screening and Identification of Transcription Factors Potentially Regulating Foxl2 Expression in Chlamys farreri Ovary.

Foxl2 is an evolutionarily conserved female sex gene, which is specifically expressed in the ovary and mainly involved in oogenesis and ovarian function maintenance. However, little is known about the mechanism that regulates Foxl2 specific expression during the ovary development. In the present study, we constructed the gonadal yeast one-hybrid (Y1H) library of Chlamysfarreri with ovaries and testes at different developmental stages using the Gateway technology. The library capacity was more than 1.36 × 107 CFU, and the length of the inserted fragment was 0.75 Kb~2 Kb, which fully met the demand of yeast library screening. The highly transcriptional activity promoter sequence of C. farreri Foxl2 (Cf-Foxl2) was determined at -1000~-616 bp by dual-luciferase reporter (DLR) assay and was used as bait to screen possible transcription factors from the Y1H library. Eleven candidate factors, including five unannotated factors, were selected based on Y1H as well as their expressional differences between ovaries and testes and were verified for the first time to be involved in the transcriptional regulation of Cf-Foxl2 by RT-qPCR and DLR. Our findings provided valuable data for further studying the specific regulation mechanism of Foxl2 in the ovary.

Growth Differences of Growth Hormone Transgenic Female and Male Channel Catfish, Ictalurus punctatus, Grown in Earthen Ponds to Sexual Maturation.

This study compared growth performance between female and male transgenic channel catfish, Ictalurus punctatus, containing channel catfish growth hormone full-length cDNA driven by the ocean pout antifreeze protein promoter, opAFP-ccGH, the rainbow trout metallothionein promoter, rtMT-ccGH, or both constructs, and their non-transgenic siblings in earthen ponds at 16 and 48 months of age. Body weight between the transgenic and their non-transgenic siblings differed (P < 0.001) at all ages. Transgenic F2 opAFP-ccGH grew 1.51- to 2.58-, F2 rtMT-ccGH grew 1.44- to 2.99- and F1fish transgenic for both constructs grew 1.36- to 2.92- fold larger than their non-transgenic sibling controls, depending upon age and sex. Body weight of the transgenic GH males was significantly higher than those of the transgenic GH females at 16 months of age (P < 0.001). However, body weight of the transgenic GH females was significantly higher (P < 0.001) compared with those of the transgenic GH males at 48 months of age, but not for the double transgenics (P > 0.05). In the case of non-transgenic GH siblings, males were larger than females at both 16 and 48 months of age (P < 0.001). Sexually dimorphic responses to GH transgenes were the opposite after sexual maturation. When critically low dissolved oxygen levels were encountered, survival of transgenic male and female opAFP-ccGH channel catfish was lower than that of controls (P = 0.004), as well as rtMT-ccGH females (P = 0.11), which is not surprising since the largest fish are most likely to succumb during an oxygen depletion.

Direct and pleiotropic effects of the Masou Salmon Delta-5 Desaturase transgene in F1 channel catfish (Ictalurus punctatus).

Channel catfish (Ictalurus punctatus) is the primary culture species in the US along with its hybrid made with male blue catfish, I. furcatus. In an effort to improve the nutritional value of channel catfish, the masou salmon Δ5-desaturase like gene (D5D) driven by the common carp beta-actin promoter (ßactin) was inserted into channel catfish. The objectives of this study were to determine the effectiveness of ßactin-D5D for improving n-3 fatty acid production in F1 transgenic channel catfish, as well as examine pleiotropic effects on growth, proximate analysis, disease resistance, and other performance traits. Transgenic F1 channel catfish showed a 33% increase in the relative proportion of n-3 fatty acids coupled with a 15% decrease in n-6 fatty acids and a 17% decrease in n-9 fatty acids when compared to non-transgenic full-siblings (P < 0.01, P < 0.01, P < 0.01). However, while the relative proportion of n-3 fatty acids was achieved, the total amount of fatty acids in the transgenic fish decreased resulting in a reduction of all fatty acids. Insertion of the ßactin-D5D transgene into channel catfish also had large effects on the body composition, and growth of channel catfish. Transgenic channel catfish grew faster, were more disease resistant, had higher protein and moisture percentage, but lower fat percentage than full-sib controls. There were sex effects as performance changes were more dramatic and significant in males. The ßactin-D5D transgenic channel catfish were also more uniform in their fatty acid composition, growth and other traits.

A novel transcription factor MRPS27 up-regulates the expression of sqr, a key gene of mitochondrial sulfide metabolism in echiuran worm Urechis unicinctus.

Hydrogen sulfide is a natural, widely distributed, poisonous substance and sulfide: quinone oxidoreductase (SQR) is responsible for oxidizing hydrogen sulfide to less toxic sulfur compounds. The increase of SQR mRNA level is an important mechanism for organisms to adapt to hydrogen sulfide-rich environments. However, its transcriptional regulation mechanism is not very clear. In this study, a mitochondrial 28S ribosomal protein S27 (MRPS27), which has never been reported as a transcription factor, was screened by yeast one-hybrid experiment from the echiuran worm Urechis unicinctus, a benthic organism living in marine sediments. Western blotting indicated that UuMRPS27 contents increased significantly in the nuclear extract of hindgut under exposed to 150 µM sulfide. ChIP and EMSA assays demonstrated that UuMRPS27 did bind to the sqr proximal promoter, the key binding sequence was CTAGAG (+12 to +17 of the promoter) detected by DNase I footprinting assay as well as transient transfection experiments. Furthermore, UuMRPS27, as a transcription activator, exhibited the highest transcription activity compared with other reported sqr transcription factors. Our data revealed for the first time the role of MRPS27 acting as a transcription factor which expanded the understanding of sqr transcriptional regulation in sulfide metabolism mechanism.

CRISPR/Cas9-mediated knock-in of alligator cathelicidin gene in a non-coding region of channel catfish genome.

CRISPR/Cas9-based gene knockout in animal cells, particularly in teleosts, has proven to be very efficient with regards to mutation rates, but the precise insertion of exogenous DNA or gene knock-in via the homology-directed repair (HDR) pathway has seldom been achieved outside of the model organisms. Here, we succeeded in integrating with high efficiency an exogenous alligator cathelicidin gene into a targeted non-coding region of channel catfish (Ictalurus punctatus) chromosome 1 using two different donor templates (synthesized linear dsDNA and cloned plasmid DNA constructs). We also tested two different promoters for driving the gene, zebrafish ubiquitin promoter and common carp ß-actin promoter, harboring a 250-bp homologous region flanking both sides of the genomic target locus. Integration rates were found higher in dead fry than in live fingerlings, indicating either off-target effects or pleiotropic effects. Furthermore, low levels of mosaicism were detected in the tissues of P1 individuals harboring the transgene, and high transgene expression was observed in the blood of some P1 fish. This can be an indication of the localization of cathelicidin in neutrophils and macrophage granules as also observed in most antimicrobial peptides. This study marks the first use of CRISPR/Cas9 HDR for gene integration in channel catfish and may contribute to the generation of a more efficient system for precise gene integration in catfish and other aquaculture species, and the development of gene-edited, disease-resistant fish.

Identification of the neuropeptide precursor genes potentially involved in the larval settlement in the Echiuran worm Urechis unicinctus.

BACKGROUND: In marine invertebrate life cycles, which often consist of planktonic larval and benthonic adult stages, settlement of the free-swimming larva to the sea floor in response to environmental cues is a key life cycle transition. Settlement is regulated by a specialized sensory-neurosecretory system, the larval apical organ. The neuroendocrine mechanisms through which the apical organ transduces environmental cues into behavioral responses during settlement are not fully understood yet. RESULTS: In this study, a total of 54 neuropeptide precursors (pNPs) were identified in the Urechis unicinctus larva and adult transcriptome databases using local BLAST and NpSearch prediction, of which 10 pNPs belonging to the ancient eumetazoa, 24 pNPs belonging to the ancient bilaterian, 3 pNPs belonging to the ancient protostome, 9 pNPs exclusive in lophotrochozoa, 3 pNPs exclusive in annelid, and 5 pNPs only found in U. unicinctus. Furthermore, four pNPs (MIP, FRWamide, FxFamide and FILamide) which may be associated with the settlement and metamorphosis of U. unicinctus larvae were analysed by qRT-PCR. Whole-mount in situ hybridization results showed that all the four pNPs were expressed in the region of the apical organ of the larva, and the positive signals were also detected in the ciliary band and abdomen chaetae. We speculated that these pNPs may regulate the movement of larval cilia and chaeta by sensing external attachment signals. CONCLUSIONS: This study represents the first comprehensive identification of neuropeptides in Echiura, and would contribute to a complete understanding on the roles of various neuropeptides in larval settlement of most marine benthonic invertebrates.

Publisher Correction: Evolutionary transcriptomics of metazoan biphasic life cycle supports a single intercalation origin of metazoan larvae.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Evolutionary transcriptomics of metazoan biphasic life cycle supports a single intercalation origin of metazoan larvae.

The transient larva-bearing biphasic life cycle is the hallmark of many metazoan phyla, but how metazoan larvae originated remains a major enigma in animal evolution. There are two hypotheses for larval origin. The 'larva-first' hypothesis suggests that the first metazoans were similar to extant larvae, with later evolution of the adult-added biphasic life cycle; the 'adult-first' hypothesis suggests that the first metazoans were adult forms, with the biphasic life cycle arising later via larval intercalation. Here, we investigate the evolutionary origin of primary larvae by conducting ontogenetic transcriptome profiling for Mollusca-the largest marine phylum characterized by a trochophore larval stage and highly variable adult forms. We reveal that trochophore larvae exhibit rapid transcriptome evolution with extraordinary incorporation of novel genes (potentially contributing to adult shell evolution), and that cell signalling/communication genes (for example, caveolin and innexin) are probably crucial for larval evolution. Transcriptome age analysis of eight metazoan species reveals the wide presence of young larval transcriptomes in both trochozoans and other major metazoan lineages, therefore arguing against the prevailing larva-first hypothesis. Our findings support an adult-first evolutionary scenario with a single metazoan larval intercalation, and suggest that the first appearance of proto-larva probably occurred after the divergence of direct-developing Ctenophora from a metazoan ancestor.

Transcriptome Analysis of Larval Segment Formation and Secondary Loss in the Echiuran Worm Urechis unicinctus.

The larval segment formation and secondary loss in echiurans is a special phenomenon, which is considered to be one of the important characteristics in the evolutionary relationship between the Echiura and Annelida. To better understand the molecular mechanism of this phenomenon, we revealed the larval transcriptome profile of the echiuran worm Urechis unicinctus using RNA-Seq technology. Twelve cDNA libraries of U. unicinctus larvae, late-trochophore (LT), early-segmentation larva (ES), segmentation larva (SL), and worm-shaped larva (WL) were constructed. Totally 243,381 unigenes were assembled with an average length of 1125 bp and N50 of 1836 bp, and 149,488 unigenes (61.42%) were annotated. We obtained 70,517 differentially expressed genes (DEGs) by pairwise comparison of the larval transcriptome data at different developmental stages and clustered them into 20 gene expression profiles using STEM software. Based on the typical profiles during the larval segment formation and secondary loss, eight signaling pathways were enriched, and five of which, mTOR, PI3K-AKT, TGF-ß, MAPK, and Dorso-ventral axis formation signaling pathway, were proposed for the first time to be involved in the segment formation. Furthermore, we identified 119 unigenes related to the segment formation of annelids, arthropods, and chordates, in which 101 genes were identified in Drosophila and annelids. The function of most segment polarity gene homologs (hedgehog, wingless, engrailed, etc.) was conserved in echiurans, annelids, and arthropods based on their expression profiles, while the gap and pair-rule gene homologs were not. Finally, we verified that strong positive signals of Hedgehog were indeed located on the boundary of larval segments using immunofluorescence. Data in this study provide molecular evidence for the understanding of larval segment development in echiurans and may serve as a blueprint for segmented ancestors in future research.

Analysis of chromosome karyotype and genome size in echiuran Urechisunicinctus Drasche, 1880 (Polychaeta, Urechidae).

Karyotype and genome size are two primary cytogenetic characteristics of species, which are of great significance to the study of cytogenetics, taxonomy, phylogenesis, evolution as well as molecular biology. However, this basic cytogenetic information in echiurans is lacking. Therefore, we analyzed characteristics of karyotype and genome size in the echiuran worm Urechisunicinctus Drasche, 1880. In this study, coelomic cells of U.unicinctus were used for analyzing the genome size by a flow cytometry with chicken erythrocytes as DNA standard, and the 2C DNA content was determined to be 1.85 pg, which was corresponded to the genome size of 904.58 Mbp approximately. Furthermore, trochophores of U.unicinctus were dissociated and cells were utilized for preparing the chromosomes stained with DAPI, and the karyotype was determined as 2n = 30 (10m + 6sm + 6st + 8t), FN=52. Our data provided the basic cytogenetic information of U.unicinctus, which could be utilized in taxonomic study and whole-genome sequencing in future.