De novo transcriptome assembly, functional annotation and differential gene expression analysis of juvenile and adult E. fetida, a model oligochaete used in ecotoxicological studies

A Novel Method to Profile Transcripts Encoding SH2 Domains in the Patiria miniata Mature Egg Transcriptome

The critical mechanism to restart zygote metabolism and prevent polyspermy during fertilization is the intracellular Ca2+ increase. All of the signaling molecules leading to the Ca2+ rise are not fully known in any species. In the sea star Patiria miniata, SFK1, SFK3, and PLCγ participate in this fertilization Ca2+ increase. These proteins share common regulatory features, including signaling via tyrosine phosphorylation and their SH2 domains. In this study, we explore two different bioinformatic strategies to identify transcripts in the Patiria miniata mature egg transcriptome (Accession PRJNA398668) that code for proteins possessing an SH2 domain. The first identified the longest open reading frame for each transcript and then utilized similarity searching tools to provide identities for each transcript. The second, novel, method involved a six-frame translation of the entire transcriptome to identify SH2 domain-containing proteins. The identified transcripts were aligned against the NCBI non-redundant database and the SwissProt database. Eighty-two transcripts that encoded SH2 domains were identified. Of these, 33 were only found using the novel method. This work furthers research into egg activation by providing possible target proteins for future experiments and a novel method for identifying specific proteins of interest within a de novo transcriptome.

Identification, characterization, and validation of NBS-encoding genes in grass pea

Grass pea is a promising crop with the potential to provide food and fodder, but its genomics has not been adequately explored. Identifying genes for desirable traits, such as drought tolerance and disease resistance, is critical for improving the plant. Grass pea currently lacks known R-genes, including the nucleotide-binding site-leucine-rich repeat (NBS-LRR) gene family, which plays a key role in protecting the plant from biotic and abiotic stresses. In our study, we used the recently published grass pea genome and available transcriptomic data to identify 274 NBS-LRR genes. The evolutionary relationships between the classified genes on the reported plants and LsNBS revealed that 124 genes have TNL domains, while 150 genes have CNL domains. All genes contained exons, ranging from 1 to 7. Ten conserved motifs with lengths ranging from 16 to 30 amino acids were identified. We found TIR-domain-containing genes in 132 LsNBSs, with 63 TIR-1 and 69 TIR-2, and RX-CCLike in 84 LsNBSs. We also identified several popular motifs, including P-loop, Uup, kinase-GTPase, ABC, ChvD, CDC6, Rnase_H, Smc, CDC48, and SpoVK. According to the gene enrichment analysis, the identified genes undergo several biological processes such as plant defense, innate immunity, hydrolase activity, and DNA binding. In the upstream regions, 103 transcription factors were identified that govern the transcription of nearby genes affecting the plant excretion of salicylic acid, methyl jasmonate, ethylene, and abscisic acid. According to RNA-Seq expression analysis, 85% of the encoded genes have high expression levels. Nine LsNBS genes were selected for qPCR under salt stress conditions. The majority of the genes showed upregulation at 50 and 200 μM NaCl. However, LsNBS-D18, LsNBS-D204, and LsNBS-D180 showed reduced or drastic downregulation compared to their respective expression levels, providing further insights into the potential functions of LsNBSs under salt stress conditions. They provide valuable insights into the potential functions of LsNBSs under salt stress conditions. Our findings also shed light on the evolution and classification of NBS-LRR genes in legumes, highlighting the potential of grass pea. Further research could focus on the functional analysis of these genes, and their potential use in breeding programs to improve the salinity, drought, and disease resistance of this important crop.

Analysis of the toxic mechanisms of fluoxastrobin on the earthworm (Eisenia fetida) using transcriptomics

Fluoxastrobin (FLUO), one of the best-selling strobilurin fungicides, could prevent fungal diseases from oilseed crops, fruits, grains, and vegetables. The widespread use of FLUO leads to the continuous accumulation of FLUO in soil. Our previous studies have demonstrated that FLUO exhibited different toxicity in artificial soil and three natural soils (fluvo-aquic soils, black soils, and red clay). The toxicity of FLUO was greater in natural soil than the artificial soil, specifically, showed the highest toxicity in fluvo-aquic soils. To better investigate the mechanism of FLUO toxicity to earthworms (Eisenia fetida), we selected fluvo-aquic soils as representative soil and used transcriptomics to study the gene expression in earthworms after FLUO exposure. The results demonstrated that the differentially expressed genes in earthworms after FLUO exposure mainly presented in pathways involving protein folding, immunity, signal transduction, and cell growth. It may be the reason why FLUO exposure stressed the earthworms and affected their normal growth activities. The present study fills gaps in the literature regarding the soil bio-toxicity of strobilurin fungicides. It also sounds the alarm for the application of such fungicides even at the low concentration (0.1 mg kg^-1).

Species-specific transcriptomic changes upon respiratory syncytial virus infection in cotton rats

The cotton rat (Sigmodon) is the gold standard pre-clinical small animal model for respiratory viral pathogens, especially for respiratory syncytial virus (RSV). However, without a reference genome or a published transcriptome, studies requiring gene expression analysis in cotton rats are severely limited. The aims of this study were to generate a comprehensive transcriptome from multiple tissues of two species of cotton rats that are commonly used as animal models (Sigmodon fulviventer and Sigmodon hispidus), and to compare and contrast gene expression changes and immune responses to RSV infection between the two species. Transcriptomes were assembled from lung, spleen, kidney, heart, and intestines for each species with a contig N50 > 1600. Annotation of contigs generated nearly 120,000 gene annotations for each species. The transcriptomes of S. fulviventer and S. hispidus were then used to assess immune response to RSV infection. We identified 238 unique genes that are significantly differentially expressed, including several genes implicated in RSV infection (e.g., Mx2, I27L2, LY6E, Viperin, Keratin 6A, ISG15, CXCL10, CXCL11, IRF9) as well as novel genes that have not previously described in RSV research (LG3BP, SYWC, ABEC1, IIGP1, CREB1). This study presents two comprehensive transcriptome references as resources for future gene expression analysis studies in the cotton rat model, as well as provides gene sequences for mechanistic characterization of molecular pathways. Overall, our results provide generalizable insights into the effect of host genetics on host-virus interactions, as well as identify new host therapeutic targets for RSV treatment and prevention.

A simple guide to de novo transcriptome assembly and annotation

A transcriptome constructed from short-read RNA sequencing (RNA-seq) is an easily attainable proxy catalog of protein-coding genes when genome assembly is unnecessary, expensive or difficult. In the absence of a sequenced genome to guide the reconstruction process, the transcriptome must be assembled de novo using only the information available in the RNA-seq reads. Subsequently, the sequences must be annotated in order to identify sequence-intrinsic and evolutionary features in them (for example, protein-coding regions). Although straightforward at first glance, de novo transcriptome assembly and annotation can quickly prove to be challenging undertakings. In addition to familiarizing themselves with the conceptual and technical intricacies of the tasks at hand and the numerous pre- and post-processing steps involved, those interested must also grapple with an overwhelmingly large choice of tools. The lack of standardized workflows, fast pace of development of new tools and techniques and paucity of authoritative literature have served to exacerbate the difficulty of the task even further. Here, we present a comprehensive overview of de novo transcriptome assembly and annotation. We discuss the procedures involved, including pre- and post-processing steps, and present a compendium of corresponding tools.

Inferring the genetic responses to acute drought stress across an ecological gradient

BACKGROUND

How do xerophytic species thrive in environments that experience extreme annual drought? Although critical to the survival of many species, the genetic responses to drought stress in many non-model organisms has yet to be explored. We investigated this question in Mentzelia section Bartonia (Loasaceae), which occurs throughout western North America, including arid lands. To better understand the genetic responses to drought stress among species that occur in different habitats, the gene expression levels of three species from Mentzelia were compared across a precipitation gradient. Two de novo reference transcriptomes were generated and annotated. Leaf and root tissues were collected from control and drought shocked plants and compared to one another for differential expression. A target-gene approach was also implemented to better understand how drought-related genes from model and crop species function in non-model systems.

RESULTS

When comparing the drought-shock treatment plants to their respective control plants, we identified 165 differentially expressed clusters across all three species. Differentially expressed genes including those associated with water movement, photosynthesis, and delayed senescence. The transcriptome profiling approach was coupled with a target genes approach that measured expression of 90 genes associated with drought tolerance in model organisms. Comparing differentially expressed genes with a ≥ 2 log-fold value between species and tissue types showed significant differences in drought response. In pairwise comparisons, species that occurred in drier environments differentially expressed greater genes in leaves when drought shocked than those from wetter environments, but expression in the roots mostly produced opposite results.

CONCLUSIONS

Arid-adapted species mount greater genetic responses compared to the mesophytic species, which has likely evolved in response to consistent annual drought exposure across generations. Drought responses also depended on organ type. Xerophytes, for example, mounted a larger response in leaves to downregulate photosynthesis and senescence, while mobilizing carbon and regulating water in the roots. The complexity of drought responses in Mentzelia suggest that whole organism responses need to be considered when studying drought and, in particular, the physiological mechanisms in which plants regulate water, carbon, cell death, metabolism, and secondary metabolites.

The transcriptome of anterior regeneration in earthworm Eudrilus eugeniae

The oligochaete earthworm, Eudrilus eugeniae is capable of regenerating both anterior and posterior segments. The present study focuses on the transcriptome analysis of earthworm E. eugeniae to identify and functionally annotate the key genes supporting the anterior blastema formation and regulating the anterior regeneration of the worm. The Illumina sequencing generated a total of 91,593,182 raw reads which were assembled into 105,193 contigs using CLC genomics workbench. In total, 40,946 contigs were annotated against the NCBI nr and SwissProt database and among them, 15,702 contigs were assigned to 14,575 GO terms. Besides a total of 9389 contigs were mapped to 416 KEGG biological pathways. The RNA-Seq comparison study identified 10,868 differentially expressed genes (DEGs) and of them, 3986 genes were significantly upregulated in the anterior regenerated blastema tissue samples of the worm. The GO enrichment analysis showed angiogenesis and unfolded protein binding as the top enriched functions and the pathway enrichment analysis denoted TCA cycle as the most significantly enriched pathway associated with the upregulated gene dataset of the worm. The identified DEGs and their function and pathway information can be effectively utilized further to interpret the key cellular, genetic and molecular events associated with the regeneration of the worm.

Tri-n-butyl phosphate induced earthworm intestinal damage by influencing nutrient absorption and energy homeostasis of intestinal epithelial cells

Tri-n-butyl phosphate (TnBP) is a typical alkyl organophosphate ester that has been used for decades in various products. However, toxicity on terrestrial organisms induced by TnBP has been rarely reported though soil is a predominant sink for hydrophobic organic compounds. The objective of this study was to investigate the TnBP-induced intestinal toxicity mechanism on earthworm Eisenia fetida as well as the potential role of gut bacteria on host's health. TnBP was found to have high bioconcentrations in earthworm intestinal tract. Digestive tract degradation and digestive enzyme activities disruption associated with nutrients absorption were noticed. Using multi-omics approaches, detailed intestinal toxic mechanism of earthworms under TnBP exposure was provided. Tight junctions between small intestinal epithelial cells and osmotic equilibrium were destroyed under 10 mg/kg TnBP, leading to nutrient absorption disturbance. To satisfy the excessive energy requirements induced by TnBP, amino acids gluconeogenesis and protein degradation were detected. Moreover, TnBP significantly decreased the diversity of gut microbiota and changed their structure and function involved in hosts' health and nutrients supply. Overall, this study provides insight into the molecular mechanism of intestinal toxicity by which earthworms respond to TnBP exposure and offer important information for risk assessment of organophosphate esters on soil ecosystems.

De novo assembly and characterization of transcriptome towards understanding molecular mechanism associated with MYMIV-resistance in Vigna mungo - A computational study.. [DOI: 10.1101/844639] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

The fast climate change affects yield in Vigna mungo via enhancing both biotic and abiotic stresses. Out of all factors, the yellow mosaic disease has the most damaging effect. However, due to lack of reference genome of Vigna mungo, the complete mechanism associated with MYMIV (Mungbean Yellow Mosaic Indian Virus) resistance in Vigna mungo remain elusive to date. Considering this, the authors made an attempt to release new transcriptome and its annotation by employing computational approaches. Quality assessment of the generated transcriptomes reveals that it successfully aligned with 99.03% of the raw reads and hence can be employed for future research. Functional annotation of the transcriptome reveals that 31% and ∼14% of the total transcripts encode lncRNAs and protein-coding sequences, respectively. Further, analysis reveals that, out of total transcripts, only 4536 and 78808 are significantly down and up-regulated during MYMIV infection in Vigna mungo, respectively. These significant transcripts are mainly associated with ribosome, spliceosome, glycolysis /gluconeogenesis, RNA transport, oxidative phosphorylation, protein processing in the endoplasmic reticulum, MAPK signaling pathway - plant, methionine and cysteine metabolism, purine metabolism and RNA degradation. Unlike the previous study, this is for the first time, the present study identified these pathways may play key role in MYMIV resistance in Vigna mungo. Thus, information and transcriptomes data available in the present study make a significant contribution to understanding the genomic structure of Vigna mungo, enabling future analyses as well as downstream applications of gene expression, sequence evolution, and genome annotation.

De novo assembly and characterization of transcriptome towards understanding molecular mechanism associated with MYMIV-resistance in Vigna mungo - A computational study.. [DOI: 10.1101/844639]

The fast climate change affects yield in Vigna mungo via enhancing both biotic and abiotic stresses. Out of all factors, the yellow mosaic disease has the most damaging effect. However, due to lack of reference genome of Vigna mungo, the complete mechanism associated with MYMIV (Mungbean Yellow Mosaic Indian Virus) resistance in Vigna mungo remain elusive to date. Considering this, the authors made an attempt to release new transcriptome and its annotation by employing computational approaches. Quality assessment of the generated transcriptomes reveals that it successfully aligned with 99.03% of the raw reads and hence can be employed for future research. Functional annotation of the transcriptome reveals that 31% and ∼14% of the total transcripts encode lncRNAs and protein-coding sequences, respectively. Further, analysis reveals that, out of total transcripts, only 4536 and 78808 are significantly down and up-regulated during MYMIV infection in Vigna mungo, respectively. These significant transcripts are mainly associated with ribosome, spliceosome, glycolysis /gluconeogenesis, RNA transport, oxidative phosphorylation, protein processing in the endoplasmic reticulum, MAPK signaling pathway - plant, methionine and cysteine metabolism, purine metabolism and RNA degradation. Unlike the previous study, this is for the first time, the present study identified these pathways may play key role in MYMIV resistance in Vigna mungo. Thus, information and transcriptomes data available in the present study make a significant contribution to understanding the genomic structure of Vigna mungo, enabling future analyses as well as downstream applications of gene expression, sequence evolution, and genome annotation.

De novo assembly and characterization of transcriptome towards understanding molecular mechanism associated with MYMIV-resistance in Vigna mungo - A computational study.. [DOI: 10.1101/844639] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

The fast climate change affects yield in Vigna mungo via enhancing both biotic and abiotic stresses. Out of all factors, the yellow mosaic disease has the most damaging effect. However, due to lack of reference genome of Vigna mungo, the complete mechanism associated with MYMIV (Mungbean Yellow Mosaic Indian Virus) resistance in Vigna mungo remain elusive to date. Considering this, the authors made an attempt to release new transcriptome and its annotation by employing computational approaches. Quality assessment of the generated transcriptomes reveals that it successfully aligned with 99.03% of the raw reads and hence can be employed for future research. Functional annotation of the transcriptome reveals that 31% and ∼14% of the total transcripts encode lncRNAs and protein-coding sequences, respectively. Further, analysis reveals that, out of total transcripts, only 4536 and 78808 are significantly down and up-regulated during MYMIV infection in Vigna mungo, respectively. These significant transcripts are mainly associated with ribosome, spliceosome, glycolysis /gluconeogenesis, RNA transport, oxidative phosphorylation, protein processing in the endoplasmic reticulum, MAPK signaling pathway - plant, methionine and cysteine metabolism, purine metabolism and RNA degradation. Unlike the previous study, this is for the first time, the present study identified these pathways may play key role in MYMIV resistance in Vigna mungo. Thus, information and transcriptomes data available in the present study make a significant contribution to understanding the genomic structure of Vigna mungo, enabling future analyses as well as downstream applications of gene expression, sequence evolution, and genome annotation.

Transcriptomic responses of regenerating earthworms (Eisenia foetida) to retinoic acid reveals the role of pluripotency genes

Exogenous retinoic acid (RA) delays and disturbs the regeneration of Eisenia foetida and inhibits the expression of pluripotent gene Sox2. However, studies of E. foetida conducted at the molecular level have been unable to elucidate its regeneration and mechanisms of RA effects on its regeneration. We merged existing transcriptomic data for E. foetida to generate a high-confidence set of transcriptomes. The de novo assembly of transcriptomes was performed by using the Trinity method, and functional annotations were analysed. We performed RNA-seq on four samples of regenerating tail fragments, three across a time-course (0, 3 and 7 days post amputation) and the fourth sample exposed to RA (7 days post amputation). E. foetida regeneration genes underwent significant upregulation and downregulation over the examined time periods, which may have been caused by a shared regulatory programme controlled by multiple gene families. The inhibition of RA against earthworm regeneration is likely related to the expression of these genes. Using annotation data and clustering, we also identified specific transcripts of 6 gene superfamilies enriched among genes exhibiting differential expression during regeneration periods and exhibiting the same expression patterns as those of the Sox2 gene. The regeneration transcriptome of tail fragment regeneration serves as a strong resource for investigating global expression changes that occur during regeneration and the toxicity of RA. This study offers insight for better understanding the regeneration of lower animals and molecular mechanisms of RA toxicity in invertebrates.

The up-regulation of two identified wound healing specific proteins-HSP70 and lysozyme in regenerated Eisenia fetida through transcriptome analysis

ETHNOPHARMACOLOGICAL RELEVANCE

Disposed earthworm has been used to treat various common ailments including burns, arthritis, itching, and inflammation for thousands of years in China. As their remarkable ability to fully regenerate in a scar-free manner, regenerated tissue homogenate of amputated Eisenia fetida (E. fetida) have been considered as an excellent wound repair therapy in our previous study. We have demonstrated that regenerated earthworm (G-90') can perform higher wound repair ability to non-regeneration tissue (G-90) through significant promotion of cutaneous wound repair in mice after their administration into wound beds.

OBJECTIVE

In the present study, we aimed to reveal the mechanism of G-90' and to explore a potential wound healing accelerated strategy.

METHODS AND RESULTS

Two functional proteins- HSP70 and lysozyme in G-90' were confirmed by cross-identification of LC-MS/MS and transcriptome analyses. Followed with semi-quantitative PCR and western blot, their expression were validated to up-regulate in 3-day regenerated tissues (G-90').

CONCLUSION

This study implies the therapeutic potency of G-90' for wound recovery and provides a new strategy to assess other natural materials targeting wound healing with the tail-amputated E .fetida as a model organism.

Annotation of nerve cord transcriptome in earthworm Eisenia fetida

In annelid worms, the nerve cord serves as a crucial organ to control the sensory and behavioral physiology. The inadequate genome resource of earthworms has prioritized the comprehensive analysis of their transcriptome dataset to monitor the genes express in the nerve cord and predict their role in the neurotransmission and sensory perception of the species. The present study focuses on identifying the potential transcripts and predicting their functional features by annotating the transcriptome dataset of nerve cord tissues prepared by Gong et al., 2010 from the earthworm Eisenia fetida. Totally 9762 transcripts were successfully annotated against the NCBI nr database using the BLASTX algorithm and among them 7680 transcripts were assigned to a total of 44,354 GO terms. The conserve domain analysis indicated the over representation of P-loop NTPase domain and calcium binding EF-hand domain. The COG functional annotation classified 5860 transcript sequences into 25 functional categories. Further, 4502 contig sequences were found to map with 124 KEGG pathways. The annotated contig dataset exhibited 22 crucial neuropeptides having considerable matches to the marine annelid Platynereis dumerilii, suggesting their possible role in neurotransmission and neuromodulation. In addition, 108 human stem cell marker homologs were identified including the crucial epigenetic regulators, transcriptional repressors and cell cycle regulators, which may contribute to the neuronal and segmental regeneration. The complete functional annotation of this nerve cord transcriptome can be further utilized to interpret genetic and molecular mechanisms associated with neuronal development, nervous system regeneration and nerve cord function.