Morphological Observation and Transcriptome Analysis of Ciliogenesis in Urechis unicinctus (Annelida, Echiura).

During the early development of marine invertebrates, planktic larvae usually occur, and their body surfaces often form specific types of cilia that are involved in locomotion and feeding. The echiuran worm Urechis unicinctus sequentially undergoes the formation and disappearance of different types of body surface cilia during embryonic and larval development. The morphological characteristics and molecular mechanisms involved in the process remain unclear. In this study, we found that body surface cilia in U. unicinctus embryos and larvae can be distinguished into four types: body surface short cilia, apical tufts, circumoral cilia and telotrochs. Further, distribution and genesis of the body surface cilia were characterized using light microscope and electron microscope. To better understand the molecular mechanism during ciliogenesis, we revealed the embryonic and larval transcriptome profile of the key stages of ciliogenesis in U. unicinctus using RNA-Seq technology. A total of 29,158 differentially expressed genes (DEGs) were obtained from 24 cDNA libraries by RNA-Seq. KEGG pathway enrichment results showed that Notch, Wnt and Ca2+ signaling pathways were significantly enriched during the occurrence of apical tufts and circumoral cilia. Furthermore, all DEGs were classified according to their expression pattern, and DEGs with similar expression pattern were grouped into a module. All DEG co-expression modules were correlated with traits (body surface short cilia, apical tufts, circumoral cilia and telotrochs) by WGCNA, the results showed DEGs were divided into 13 modules by gene expression patterns and that the genes in No. 7, No. 8 and No. 10 modules were to be highly correlated with the occurrence of apical tufts, circumoral cilia and telotrochs. The top 10 hub genes in the above three modules were identified to be highly correlated with ciliogenesis, including the reported cilium-related gene Cnbd2 and unreported cilium-related candidate genes FAM181B, Capsl, Chst3, TMIE and Innexin. Notably, Innexin was included in the top10 hub genes of the two modules (No. 7 and No. 8), suggesting that Innexin may play an important role in U. unicinctus apical tufts, circumoral cilia and telotrochs genesis. This study revealed the characteristics of ciliogenesis on the body surface of U. unicinctus embryos and larvae, providing basic data for exploring the molecular mechanism of ciliogenesis on the body surface.

Genome-wide identification and expression analyses of Toll-like receptors provide new insights on adaptation to intertidal benthic environments in Urechis unicinctus (Annelida, Echiura).

Toll-like receptors (TLR) are an important class of pattern recognition receptors involved in innate immunity that recognize pathogen-associated and damage-associated molecular patterns. Although the role of TLRs in immunity has been extensively studied, a systematic investigation of their function in environmental adaptation is still in its infancy. In this study, a genome-wide search was conducted to systematically investigate TLR family members of Urechis unicinctus, a typical benthic organism in intertidal mudflats. A total of 28 TLR genes were identified in the U. unicinctus genome, and their fundamental physiological and biochemical properties were characterized. Gene copy number analysis among species in different habitats indicated that TLR family gene expansion may be probably related with benthic environmental adaptation. To further investigate the expression patterns of TLR members under environmental stress, transcriptome data was analyzed from different developmental stages and the hindgut under sulfide stress. Transcriptome analysis of different developmental stages showed that most TLR genes were highly expressed during a key period of benthic environment adaptation (worm-shaped larva). Transcriptome analysis of the hindgut under sulfide stress showed that the expression of 12 TLR members was significantly induced under sulfide stress. These results indicate that the regulation of TLR gene expression may be probably involved in the adaptation of U. unicinctus to the benthic intertidal zone environment. Taken together, this study may lay the foundation for future functional analysis of the specific role of TLRs in host immune responses against sulfide exposure and benthic environmental stress in annelid.

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.

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.

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.

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.

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.

SOX2 participates in spermatogenesis of Zhikong scallop Chlamys farreri.

As an important transcription factor, SOX2 involves in embryogenesis, maintenance of stem cells and proliferation of primordial germ cell (PGC). However, little was known about its function in mature gonads. Herein, we investigated the SOX2 gene profiles in testis of scallop, Chlamys farreri. The level of C. farreri SOX2 (Cf-SOX2) mRNA increased gradually along with gonadal development and reached the peak at mature stage, and was located in all germ cells, including spermatogonia, spermatocytes, spermatids and spermatozoa. Knockdown of Cf-SOX2 using RNAi leaded to a mass of germ cells lost, and only a few spermatogonia retained in the nearly empty testicular acini after 21 days. TUNEL assay showed that apoptosis occurred in spermatocytes. Furthermore, transcriptome profiles of the testes were compared between Cf-SOX2 knockdown and normal scallops, 131,340 unigenes were obtained and 2,067 differential expression genes (DEGs) were identified. GO and KEGG analysis showed that most DEGs were related to cell apoptosis (casp2, casp3, casp8), cell proliferation (samd9, crebzf, iqsec1) and spermatogenesis (htt, tusc3, zmynd10, nipbl, mfge8), and enriched in p53, TNF and apoptosis pathways. Our study revealed Cf-SOX2 is essential in spermatogenesis and testis development of C. farreri and provided important clues for better understanding of SOX2 regulatory mechanisms in bivalve testis.

A novel transcription factor Rwdd1 and its SUMOylation inhibit the expression of sqr, a key gene of mitochondrial sulfide metabolism in Urechis unicinctus.

Sulfide-quinone oxidoreductase (SQR) is a key enzyme of sulfide metabolism in metazoans, and responsible for oxidizing sulfide into thiosulfate and transmitting the generated electrons to the ubiquinone. It has been revealed that the sqr mRNA level increases significantly in echiuran worm Urechis unicinctus exposed to sulfide, and HSF1, NF1 and Sp1 have been verified to participate in its transcriptional regulation. In this study, we obtained 23 potential transcription factors interacting possibly with the proximal region (-391 to +50) of sqr promoter, and focused on the RWD domain-containing 1 (Rwdd1), a protein with the maximum number of clones in yeast one-hybrid (Y1H) screening, to investigate its transcriptional regulation to U. unincitus sqr. The ChIP and EMSA assays identified that the Rwdd1 can bind directly to the promoter (+18/+36) of U. unicinctus sqr. The point mutation and transient transfection experiments discovered that TACG was the key sequence of the DNA element bound by the Rwdd1. Furthermore, the U. unicinctus Rwdd1 (UuRwdd1) was identified to be a transcription repressor inhibiting the sqr promoter activity, and the SUMOylation of UuRwdd1 at the lysine of 90th enhanced its inhibitory effect on sqr transcription further. Western blotting found Rwdd1 responded to sulfide in hindguts from U. unincitus, and the protein content showed a remarkable drop in hindgut nuclei in the early sulfide exposure, and then increased significantly both in the total protein and the nuclear protein extract. We suggested that the Rwdd1 is a novel transcription factor, and these data improve our understanding of the sqr transcriptional regulation and the mitochondrial sulfide metabolism.