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

Discovery and characterization of the α-amylases cDNAs from Enchytraeus albidus shed light on the evolution of "Enchytraeus-Eisenia type" Amy homologs in Annelida

Although enchytraeids have gained popularity in scientific research, fundamental questions regarding their feeding ecology and biology remain largely unexplored. This study investigates α-amylases, major digestive enzymes responsible for hydrolyzing starch and similar polysaccharides into sugars, in Enchytraeus albidus. Genetic data related to α-amylases is currently lacking for the family Enchytraeidae but also for the entire Annelida. To detect and identify coding sequences of the expressed α-amylase genes in COI-monohaplotype culture (PL-A strain) of E. albidus, we used classical "gene fishing" and transcriptomic approaches. We also compared coding sequence variants of α-amylase retrieved from transcriptomic data related to freeze-tolerant strains. Our results reveal that E. albidus possesses two distinct α-amylase genes (Amy I and Amy II) that are homologs to earthworm Eisenia fetida Ef-Amy genes. Different strains of E. albidus possess distinctive alleles of α-amylases with unique SNP patterns specific to a particular strain. Unlike Amy II, Amy I seems to be a highly polymorphic and multicopy gene. The domain architecture of the putative Amy proteins was found the same as for classical animal α-amylases with ABC-domains. A characteristic feature of Amy II is the lack of GHGA motif in the flexible loop region, similarly to many insect amylases. We identified "Enchytraeus-Eisenia type" α-amylase homologs in other clitellates and polychaetes, indicating the ancestral origin of Amy I/II proteins in Annelida. This study provides the first insight into the endogenous non-proteolytic digestive enzyme genes in potworms, discusses the evolution of Amy α-amylases in Annelida, and explores phylogenetic implications.

Genome-wide profiling of DNA methylome and transcriptome reveals epigenetic regulation of Urechis unicinctus response to sulfide stress

Sulfide is a well-known environmental pollutant that can have detrimental effects on most organisms. However, few metazoans living in sulfide-rich environments have developed mechanisms to tolerate and adapt to sulfide stress. Epigenetic mechanisms, including DNA methylation, have been shown to play a vital role in environmental stress adaptation. Nevertheless, the precise function of DNA methylation in biological sulfide adaptation remains unclear. Urechis unicinctus, a benthic organism inhabiting sulfide-rich intertidal environments, is an ideal model organism for studying adaptation to sulfide environments. In this study, we conducted a comprehensive analysis of the DNA methylome and transcriptome of U. unicinctus after exposure to 50 μM sulfide. The results revealed dynamic changes in the DNA methylation (5-methylcytosine) landscape in response to sulfide stress, with U. unicinctus exhibiting elevated DNA methylation levels following stress exposure. Integrating differentially expressed genes (DEGs) and differentially methylated regions (DMRs), we identified a crucial role of gene body methylation in predicting gene expression. Furthermore, using a DNA methyltransferase inhibitor, we validated the involvement of DNA methylation in the sulfide stress response and the gene regulatory network influenced by DNA methylation. The results indicated that by modulating DNA methylation levels during sulfide stress, the expression of glutathione S-transferase, glutamyl aminopeptidase, and cytochrome c oxidase could be up-regulated, thereby facilitating the metabolism and detoxification of exogenous sulfides. Moreover, DNA methylation was found to regulate and enhance the oxidative phosphorylation pathway, including NADH dehydrogenase, isocitrate dehydrogenase, and ATP synthase. Additionally, DNA methylation influenced the regulation of Cytochrome P450 and macrophage migration inhibitory factor, both of which are closely associated with oxidative stress and stress resistance. Our findings not only emphasize the role of DNA methylation in sulfide adaptation but also provide novel insights into the potential mechanisms through which marine organisms adapt to environmental changes.

Genomic loss of the HSP70cA gene in the vertebrate lineage

Metazoan 70 kDa heat shock protein (HSP70) genes have been classified into four lineages: cytosolic A (HSP70cA), cytosolic B (HSP70cB), endoplasmic reticulum (HSP70er), and mitochondria (HSP70m). Because previous studies have identified no HSP70cA genes in vertebrates, we hypothesized that this gene was lost on the evolutionary path to vertebrates. To test this hypothesis, the present study conducted a comprehensive database search followed by phylogenetic and synteny analyses. HSP70cA genes were found in invertebrates and in two of the three subphyla of Chordata, Cephalochordata (lancelets) and Tunicata (tunicates). However, no HSP70cA gene was found in the genomes of Craniata (another subphylum of Chordata; lamprey, hagfish, elephant shark, and coelacanth), suggesting the loss of the HSP70cA gene in the early period of vertebrate evolution. Synteny analysis using available genomic resources indicated that the synteny around the HSP70 genes was generally conserved between tunicates but was largely different between tunicates and lamprey. These results suggest the presence of dynamic chromosomal rearrangement in early vertebrates that possibly caused the loss of the HSP70cA gene in the vertebrate lineage.

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.

Identification and validation of the reference genes in the echiuran worm Urechis unicinctus based on transcriptome data

BACKGROUND

Real-time quantitative PCR (RT-qPCR) is a crucial and widely used method for gene expression analysis. Selecting suitable reference genes is extremely important for the accuracy of RT-qPCR results. Commonly used reference genes are not always stable in various organisms or under different environmental conditions. With the increasing application of high-throughput sequencing, transcriptome analysis has become an effective method for identifying novel stable reference genes.

RESULTS

In this study, we identified candidate reference genes based on transcriptome data covering embryos and larvae of early development, normal adult tissues, and the hindgut under sulfide stress using the coefficient of variation (CV) method in the echiuran Urechis unicinctus, resulting in 6834 (15.82%), 7110 (16.85%) and 13880 (35.87%) candidate reference genes, respectively. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses revealed that the candidate reference genes were significantly enriched in cellular metabolic process, protein metabolic process and ribosome in early development and normal adult tissues as well as in cellular localization and endocytosis in the hindgut under sulfide stress. Subsequently, ten genes including five new candidate reference genes and five commonly used reference genes, were validated by RT-qPCR. The expression stability of the ten genes was analyzed using four methods (geNorm, NormFinder, BestKeeper, and ∆Ct). The comprehensive results indicated that the new candidate reference genes were more stable than most commonly used reference genes. The commonly used ACTB was the most unstable gene. The candidate reference genes STX12, EHMT1, and LYAG were the most stable genes in early development, normal adult tissues, and hindgut under sulfide stress, respectively. The log₂(TPM) of the transcriptome data was significantly negatively correlated with the Ct values of RT-qPCR (Ct =  - 0.5405 log₂(TPM) + 34.51), which made it possible to estimate the Ct value before RT-qPCR using transcriptome data.

CONCLUSION

Our study is the first to select reference genes for RT-qPCR from transcriptome data in Echiura and provides important information for future gene expression studies in U. unicinctus.

Phylogenetic annotation of Caenorhabditis elegans heat shock protein 70 genes

Annotation of the 70 kDa heat shock proteins (Hsp70s) has been chaotic especially in invertebrates. In this study, we validated an emerging nomenclature of Hsp70s, which can be potentially applied to all metazoan Hsp70s, by conducting a genome-wide annotation of Caenorhabditis elegans Hsp70s. Using the phylogenetic annotation, the seven canonical C. elegans Hsp70s were successfully classified into four known lineages, cytosolic A, cytosolic B, endoplasmic reticulum, and mitochondria. Motifs specific to each lineage were all conserved in the C. elegans Hsp70s. From these results, we propose new aliases of C. elegans Hsp70s that should help future annotation of this important molecular chaperone.