Onchidoris muricata is a widespread shell-less species of nudibranch molluscs, which has unique for Gastropoda skeletal elements - subepidermal calcite spicules. The general and fine morphology of the spicules, as well as their maturation process in ontogenesis, have been studied in detail by authors. The uniqueness of spicules lies in their intracellular formation and location under the ectodermal epithelium, which is more typical for deuterostomes. We present O. muricata as a potentially new model species for studying calcification of intracellular protein structure. A total of 96 individuals were collected in the Kandalaksha Bay of the White Sea, both manually and by scuba diving. All individuals were divided into three groups based on morphological characteristics such as specimens\' size, spicule condition etc. This division suggests the existence of three stages in postembryonic ontogenesis of O. muricata reflecting the maturation of the spicule complex. Total RNA samples were isolated from three size groups of molluscs in three biological replicates. Libraries were prepared from the polyadenylated RNA fraction and sequenced at NovaSeq6000 (Illumina), yielding a total of 112.8 Gb of 150 bp paired-end reads, corresponding to almost 1,000-fold coverage of the transcriptome. Representative transcriptome assembled de novo with Trinity. In addition to obtaining the transcriptome sequences of O. muricata, differential expression analysis was also performed for these three size groups. This allows us to trace the dynamics of molecular and biological processes during the life of a mollusc. The obtained data can then be used as a reference transcriptome for closely related species, to study specific expressed genes, to identify various unique sequences, including protein-coding ones, to understand biological processes, including biomineralization and much more.

Cnidarians thriving in biofouling communities on aquaculture net pens represent a significant health risk for farmed finfish due to their stinging cells. The toxins coming into contact with the fish, during net cleaning, can adversely affect their behavior, welfare, and survival, with a particularly serious health risk for the gills, causing direct tissue damage such as formation of thrombi and increasing risks of secondary infections. The hydroid Ectopleura larynx is one of the most common fouling organisms in Northern Europe. However, despite its significant economic, environmental, and operational impact on finfish aquaculture, biological information on this species is scarce and its venom composition has never been investigated. In this study, we generated a whole transcriptome of E. larynx, and identified its putative expressed venom toxin proteins (predicted toxin proteins, not functionally characterized) based on in silico transcriptome annotation mining and protein sequence analysis. The results uncovered a broad and diverse repertoire of putative toxin proteins for this hydroid species. Its toxic arsenal appears to include a wide and complex selection of toxin proteins, covering a large panel of potential biological functions that play important roles in envenomation. The putative toxins identified in this species, such as neurotoxins, GTPase toxins, metalloprotease toxins, ion channel impairing toxins, hemorrhagic toxins, serine protease toxins, phospholipase toxins, pore-forming toxins, and multifunction toxins may cause various major deleterious effects in prey, predators, and competitors. These results provide valuable new insights into the venom composition of cnidarians, and venomous marine organisms in general, and offer new opportunities for further research into novel and valuable bioactive molecules for medicine, agronomics and biotechnology.

Spider venoms have evolved over thousands of years, optimizing feeding and defense mechanisms. Venom components show pharmacological and biotechnological potential, rising interest in their study. However, the isolation of spider toxins for experimental evaluation poses significant challenges. To address this, transcriptomic analysis combined with computational tools has emerged as an appealing approach to characterizing spider venoms. However, many sequences remain unidentified after automatic annotation. In this study, we manually curated a subset of previously unannotated sequences from the Phoneutria nigriventer transcriptome and identified new putative venom components. Our manual analysis revealed 29 % of the analyzed sequences were potential venom components, 29 % hypothetical/uncharacterized proteins, and 17 % cellular function proteins. Only 25 % of the originally unannotated dataset remained without any identification. Most reclassified components were cysteine-rich peptides, including 23 novel putative toxins. We also found glycine-rich peptides (GRP), corroborating the previous description of GRPs in Phoneutria pertyi venom glands. Furthermore, to emphasize the recurrence of the lack of annotation in spider venom glands transcripts, we provide a survey of the percentage of unidentified sequences in several published spider venom transcriptomics studies. In conclusion, our study highlights the importance of manual curation in uncovering novel venom components and underscores the need for improved annotation strategies to fully exploit the medical and biotechnological potential of spider venoms.

Long noncoding RNAs (lncRNAs) are increasingly being recognized as modulators in various biological processes. However, due to their low expression, their systematic characterization is difficult to determine. Here, we performed transcript annotation by a newly developed computational pipeline, termed RNA-seq and small RNA-seq combined strategy (RSCS), in a wide variety of cellular contexts. Thousands of high-confidence potential novel transcripts were identified by the RSCS, and the reliability of the transcriptome was verified by analysis of transcript structure, base composition, and sequence complexity. Evidenced by the length comparison, the frequency of the core promoter and the polyadenylation signal motifs, and the locations of transcription start and end sites, the transcripts appear to be full length. Furthermore, taking advantage of our strategy, we identified a large number of endogenous retrovirus-associated lncRNAs, and a novel endogenous retrovirus-lncRNA that was functionally involved in control of Yap1 expression and essential for early embryogenesis was identified. In summary, the RSCS can generate a more complete and precise transcriptome, and our findings greatly expanded the transcriptome annotation for the mammalian community.

BACKGROUND: RNA-seq followed by de novo transcriptome assembly has been a transformative technique in biological research of non-model organisms, but the computational processing of RNA-seq data entails many different software tools. The complexity of these de novo transcriptomics workflows therefore presents a major barrier for researchers to adopt best-practice methods and up-to-date versions of software.
RESULTS: Here we present a streamlined and universal de novo transcriptome assembly and annotation pipeline, transXpress, implemented in Snakemake. transXpress supports two popular assembly programs, Trinity and rnaSPAdes, and allows parallel execution on heterogeneous cluster computing hardware.
CONCLUSIONS: transXpress simplifies the use of best-practice methods and up-to-date software for de novo transcriptome assembly, and produces standardized output files that can be mined using SequenceServer to facilitate rapid discovery of new genes and proteins in non-model organisms.

The advancement of hybrid sequencing technologies is increasingly expanding genome assemblies that are often annotated using hybrid sequencing transcriptomics, leading to improved genome characterization and the identification of novel genes and isoforms in a wide variety of organisms.
We developed an easy-to-use genome-guided transcriptome annotation pipeline that uses assembled transcripts from hybrid sequencing data as input and distinguishes between coding and long non-coding RNAs by integration of several bioinformatic approaches, including gene reconciliation with previous annotations in GTF format. We demonstrated the efficiency of this approach by correctly assembling and annotating all exons from the chicken SCO-spondin gene (containing more than 105 exons), including the identification of missing genes in the chicken reference annotations by homology assignments.
Our method helps to improve the current transcriptome annotation of the chicken brain. Our pipeline, implemented on Anaconda/Nextflow and Docker is an easy-to-use package that can be applied to a broad range of species, tissues, and research areas helping to improve and reconcile current annotations. The code and datasets are publicly available at https://github.com/cfarkas/annotate_my_genomes.

The annotation of animal genomes plays an important role in elucidating molecular mechanisms behind the genetic control of economically important traits. Here, we employed long-read sequencing technology, Oxford Nanopore Technology, to annotate the pig transcriptome across 17 tissues from two Yorkshire littermate pigs. More than 9.8 million reads were obtained from a single flow cell, and 69 781 unique transcripts at 50 108 loci were identified. Of these transcripts, 16 255 were found to be novel isoforms, and 22 344 were found at loci that were novel and unannotated in the Ensembl (release 102) and NCBI (release 106) annotations. Novel transcripts were mostly expressed in cerebellum, followed by lung, liver, spleen, and hypothalamus. By comparing the unannotated transcripts to existing databases, there were 21 285 (95.3%) transcripts matched to the NT database (v5) and 13 676 (61.2%) matched to the NR database (v5). Moreover, there were 4324 (19.4%) transcripts matched to the SwissProt database (v5), corresponding to 11 356 proteins. Tissue-specific gene expression analyses showed that 9749 transcripts were highly tissue-specific, and cerebellum contained the most tissue-specific transcripts. As the same samples were used for the annotation of cis-regulatory elements in the pig genome, the transcriptome annotation generated by this study provides an additional and complementary annotation resource for the Functional Annotation of Animal Genomes effort to comprehensively annotate the pig genome.

Egg-laying rate is mainly determined by ovarian function and regulated by the hypothalamic-pituitary-gonadal axis; however, the mechanism by which the ovary regulates the egg-laying rate is still poorly understood. The purpose of this study was to compare the differences in the transcriptomes of the ovary of Lingyun black-bone chickens with relatively high and low egg-laying rates and screen candidate genes related to the egg-laying rate. RNA-sequencing (RNA-Seq) was conducted to explore the chicken transcriptome from the ovarian tissue of six Lingyun black-bone chickens with high (group G, n = 3) and low (group D, n = 3) egg-laying rates. The results showed that 235 differentially expressed genes (DEGs) were identified between the chickens with high and low egg-laying rates; among them, 209 DEGs were up-regulated and 26 DEGs were down-regulated. Gene Ontology analysis showed that the up-regulated 209 DEGs were enriched in 50 GO terms and the down-regulated 26 DEGs were enriched in 40 GO terms. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that up-regulated DEGs were significantly enriched in 25 pathways and down-regulated DEGs were significantly enriched in three pathways. Among the pathways, we found the longevity regulating pathway-multiple species pathway, Estrogen signalling pathway and PPAR signalling pathway may have an essential function in regulating the egg-laying rate. The results highlighted DEGs in the ovarian tissues of relatively high and low laying Lingyun black-bone chicken and identified essential candidate genes related to the egg-laying rate, thereby providing a theoretical basis for improving the egg-laying rate of Lingyun black-bone chicken.

Asian citrus psyllid Diaphorina citri Kuwayama is an important economic pest of citrus, as it transmits Candidatus Liberibacter asiaticus, the causative agent of huanglongbing. In this study, we used RNA-seq to identify novel genes and provide the first high-resolution view of the of D. citri transcriptome throughout development. The transcriptomes of D. citri during eight developmental stages, including the egg, five instars, and male and female adults were sequenced. In total, 115 million clean reads were obtained and assembled into 354,726 unigenes with an average length of 925.65 bp and an N50 length of 1733 bp. Clusters of Orthologous Groups, Gene Ontology, and Kyoto Encyclopedia of Genes and Genomes analyses were conducted to functionally annotate the genes. Differential expression analysis highlighted developmental stage-specific expression patterns. Furthermore, two trehalase genes were characterized with lower expression in adults compared to that in the other stages. The RNA interference (RNAi)-mediated suppression of the two trehalase genes resulted in significantly high D. citri mortality. This study enriched the genomic information regarding D. citri. Importantly, these data represent the most comprehensive transcriptomic resource currently available for D. citri and will facilitate functional genomics studies of this notorious pest.

Henosepilachna vigintioctopunctata (Coleoptera: Coccinellidae) is a major pest affecting Solanaceae plants in Asian countries. In this study, we sequenced the ovary and testis transcriptomes of H. vigintioctopunctata to identify gonad-related genes. Comparison of the unigene sequences in ovary and testis libraries identified 1,421 and 5,315 ovary- and testis-specific genes, respectively. Among the ovary-specific genes, we selected the RC2-like and PSHS-like genes to investigate the effects of gene silencing on the mortality, percentage infertility, pre-oviposition period, fecundity, daily number of eggs laid, and hatching rate in female adults. Although the percentage mortality and infertility of females did not differ significantly among dsRNA treatments, fecundity was significantly reduced in the dsRC2-like and dsPSHS-like treatment groups. Moreover, the pre-oviposition period was markedly prolonged in response to dsPSHS-like treatment. This is the first reported RNA sequencing of H. vigintioctopunctata. The transcriptome sequences and gene expression profiles of the ovary and testis libraries will provide useful information for the identification of gonad-related genes in H. vigintioctopunctata and facilitate further research on the reproductive biology of this species. Moreover, the gonad-specific genes identified may represent candidate target genes for inhibiting the population growth of H. vigintioctopunctata.