Eukaryotic initiation factor EIF-3.G augments mRNA translation efficiency to regulate neuronal activity.

The translation initiation complex eIF3 imparts specialized functions to regulate protein expression. However, understanding of eIF3 activities in neurons remains limited despite widespread dysregulation of eIF3 subunits in neurological disorders. Here, we report a selective role of the C. elegans RNA-binding subunit EIF-3.G in shaping the neuronal protein landscape. We identify a missense mutation in the conserved Zinc-Finger (ZF) of EIF-3.G that acts in a gain-of-function manner to dampen neuronal hyperexcitation. Using neuron-type-specific seCLIP, we systematically mapped EIF-3.G-mRNA interactions and identified EIF-3.G occupancy on GC-rich 5'UTRs of a select set of mRNAs enriched in activity-dependent functions. We demonstrate that the ZF mutation in EIF-3.G alters translation in a 5'UTR-dependent manner. Our study reveals an in vivo mechanism for eIF3 in governing neuronal protein levels to control neuronal activity states and offers insights into how eIF3 dysregulation contributes to neurological disorders.

Profiling microRNAs through development of the parasitic nematode Haemonchus identifies nematode-specific miRNAs that suppress larval development.

Parasitic nematodes transition between dramatically different free-living and parasitic stages, with correctly timed development and migration crucial to successful completion of their lifecycle. However little is known of the mechanisms controlling these transitions. microRNAs (miRNAs) negatively regulate gene expression post-transcriptionally and regulate development of diverse organisms. Here we used microarrays to determine the expression profile of miRNAs through development and in gut tissue of the pathogenic nematode Haemonchus contortus. Two miRNAs, mir-228 and mir-235, were enriched in infective L3 larvae, an arrested stage analogous to Caenorhabditis elegans dauer larvae. We hypothesized that these miRNAs may suppress development and maintain arrest. Consistent with this, inhibitors of these miRNAs promoted H. contortus development from L3 to L4 stage, while genetic deletion of C. elegans homologous miRNAs reduced dauer arrest. Epistasis studies with C. elegans daf-2 mutants showed that mir-228 and mir-235 synergise with FOXO transcription factor DAF-16 in the insulin signaling pathway. Target prediction suggests that these miRNAs suppress metabolic and transcription factor activity required for development. Our results provide novel insight into the expression and functions of specific miRNAs in regulating nematode development and identify miRNAs and their target genes as potential therapeutic targets to limit parasite survival within the host.

The C. elegans 3' UTRome v2 resource for studying mRNA cleavage and polyadenylation, 3'-UTR biology, and miRNA targeting.

3' Untranslated regions (3' UTRs) of mRNAs emerged as central regulators of cellular function because they contain important but poorly characterized cis-regulatory elements targeted by a multitude of regulatory factors. The model nematode Caenorhabditis elegans is ideal to study these interactions because it possesses a well-defined 3' UTRome. To improve its annotation, we have used a genome-wide bioinformatics approach to download raw transcriptome data for 1088 transcriptome data sets corresponding to the entire collection of C. elegans trancriptomes from 2015 to 2018 from the Sequence Read Archive at the NCBI. We then extracted and mapped high-quality 3'-UTR data at ultradeep coverage. Here, we describe and release to the community the updated version of the worm 3' UTRome, which we named 3' UTRome v2. This resource contains high-quality 3'-UTR data mapped at single-base ultraresolution for 23,084 3'-UTR isoform variants corresponding to 14,788 protein-coding genes and is updated to the latest release of WormBase. We used this data set to study and probe principles of mRNA cleavage and polyadenylation in C. elegans The worm 3' UTRome v2 represents the most comprehensive and high-resolution 3'-UTR data set available in C. elegans and provides a novel resource to investigate the mRNA cleavage and polyadenylation reaction, 3'-UTR biology, and miRNA targeting in a living organism.

Neuronal Small RNAs Control Behavior Transgenerationally.

It is unknown whether the activity of the nervous system can be inherited. In Caenorhabditis elegans nematodes, parental responses can transmit heritable small RNAs that regulate gene expression transgenerationally. In this study, we show that a neuronal process can impact the next generations. Neurons-specific synthesis of RDE-4-dependent small RNAs regulates germline amplified endogenous small interfering RNAs (siRNAs) and germline gene expression for multiple generations. Further, the production of small RNAs in neurons controls the chemotaxis behavior of the progeny for at least three generations via the germline Argonaute HRDE-1. Among the targets of these small RNAs, we identified the conserved gene saeg-2, which is transgenerationally downregulated in the germline. Silencing of saeg-2 following neuronal small RNA biogenesis is required for chemotaxis under stress. Thus, we propose a small-RNA-based mechanism for communication of neuronal processes transgenerationally.

Detection of microRNA-Target Interactions by Chimera PCR (ChimP).

MicroRNAs (miRNAs) regulate gene expression by directing Argonaute proteins to target RNAs, which usually results in destabilization and translational inhibition of the target RNA. The prediction of animal miRNA target sites has remained a challenge due to the ability of miRNAs to bind target RNAs through imperfect base pairing. Recently, several labs have established methods to produce biochemical evidence of miRNA-target interactions by generating chimeric reads where the miRNA is ligated to its target RNA. Despite the insights that can be gained from chimera producing methods, the current approaches are inefficient, labor intensive and require computational expertise. Here we describe a method, called Chimera PCR (ChimP), for the validation or testing of specific miRNA-target interactions. This method allows for focused experiments to analyze miRNA targeting in a variety of conditions.

GTSF-1 is required for formation of a functional RNA-dependent RNA Polymerase complex in Caenorhabditis elegans.

Argonaute proteins and their associated small RNAs (sRNAs) are evolutionarily conserved regulators of gene expression. Gametocyte-specific factor 1 (Gtsf1) proteins, characterized by two tandem CHHC zinc fingers and an unstructured C-terminal tail, are conserved in animals and have been shown to interact with Piwi clade Argonautes, thereby assisting their activity. We identified the Caenorhabditis elegans Gtsf1 homolog, named it gtsf-1 and characterized it in the context of the sRNA pathways of C. elegans We report that GTSF-1 is not required for Piwi-mediated gene silencing. Instead, gtsf-1 mutants show a striking depletion of 26G-RNAs, a class of endogenous sRNAs, fully phenocopying rrf-3 mutants. We show, both in vivo and in vitro, that GTSF-1 interacts with RRF-3 via its CHHC zinc fingers. Furthermore, we demonstrate that GTSF-1 is required for the assembly of a larger RRF-3 and DCR-1-containing complex (ERIC), thereby allowing for 26G-RNA generation. We propose that GTSF-1 homologs may act to drive the assembly of larger complexes that act in sRNA production and/or in imposing sRNA-mediated silencing activities.

Organ specific gene expression in the regenerating tail of Macrostomum lignano.

Temporal and spatial characterization of gene expression is a prerequisite for the understanding of cell-, tissue-, and organ-differentiation. In a multifaceted approach to investigate gene expression in the tail plate of the free-living marine flatworm Macrostomum lignano, we performed a posterior-region-specific in situ hybridization screen, RNA sequencing (RNA-seq) of regenerating animals, and functional analyses of selected tail-specific genes. The in situ screen revealed transcripts expressed in the antrum, cement glands, adhesive organs, prostate glands, rhabdite glands, and other tissues. Next we used RNA-seq to characterize temporal expression in the regenerating tail plate revealing a time restricted onset of both adhesive organs and copulatory apparatus regeneration. In addition, we identified three novel previously unannotated genes solely expressed in the regenerating stylet. RNA interference showed that these genes are required for the formation of not only the stylet but the whole male copulatory apparatus. RNAi treated animals lacked the stylet, vesicula granulorum, seminal vesicle, false seminal vesicle, and prostate glands, while the other tissues of the tail plate, such as adhesive organs regenerated normally. In summary, our findings provide a large resource of expression data during homeostasis and regeneration of the morphologically complex tail regeneration and pave the way for a better understanding of organogenesis in M. lignano.

Comparative genome analysis of programmed DNA elimination in nematodes.

Programmed DNA elimination is a developmentally regulated process leading to the reproducible loss of specific genomic sequences. DNA elimination occurs in unicellular ciliates and a variety of metazoans, including invertebrates and vertebrates. In metazoa, DNA elimination typically occurs in somatic cells during early development, leaving the germline genome intact. Reference genomes for metazoa that undergo DNA elimination are not available. Here, we generated germline and somatic reference genome sequences of the DNA eliminating pig parasitic nematode Ascaris suum and the horse parasite Parascaris univalens. In addition, we carried out in-depth analyses of DNA elimination in the parasitic nematode of humans, Ascaris lumbricoides, and the parasitic nematode of dogs, Toxocara canis. Our analysis of nematode DNA elimination reveals that in all species, repetitive sequences (that differ among the genera) and germline-expressed genes (approximately 1000-2000 or 5%-10% of the genes) are eliminated. Thirty-five percent of these eliminated genes are conserved among these nematodes, defining a core set of eliminated genes that are preferentially expressed during spermatogenesis. Our analysis supports the view that DNA elimination in nematodes silences germline-expressed genes. Over half of the chromosome break sites are conserved between Ascaris and Parascaris, whereas only 10% are conserved in the more divergent T. canis. Analysis of the chromosomal breakage regions suggests a sequence-independent mechanism for DNA breakage followed by telomere healing, with the formation of more accessible chromatin in the break regions prior to DNA elimination. Our genome assemblies and annotations also provide comprehensive resources for analysis of DNA elimination, parasitology research, and comparative nematode genome and epigenome studies.

Comparative sequence analysis reveals regulation of genes in developing schistosomula of Schistosoma mansoni exposed to host portal serum.

Once inside a vertebrate host after infection, individual schistosomula of the parasite Schistosoma mansoni find a new and complex environment, which requires quick adjustments for survival, such as those that allow it to avoid the innate immune response of the host. Thus, it is very important for the parasite to remain within the skin after entering the host for a period of about 3 days, at which time it can then reach the venous system, migrate to the lungs and, by the end of eighth day post-infection, it reach the portal venous system, while undergoing minimal changes in morphology. However, after just a few days in the portal blood system, the parasite experiences an extraordinary increase in biomass and significant morphological alterations. Therefore, determining the constituents of the portal venous system that may trigger these changes that causes the parasite to consolidate its development inside the vertebrate host, thus causing the disease schistosomiasis, is essential. The present work simulated the conditions found in the portal venous system of the vertebrate host by exposing schistosomula of S. mansoni to in vitro culture in the presence of portal serum of the hamster, Mesocricetus auratus. Two different incubation periods were evaluated, one of 3 hours and one of 12 hours. These time periods were used to mimic the early contact of the parasite with portal serum during the course of natural infection. As a control, parasites were incubated in presence of hamster peripheral serum, in order to compare gene expression signatures between the two conditions. The mRNA obtained from parasites cultured under both conditions were submitted to a whole transcriptome library preparation and sequenced with a next generation platform. On average, nearly 15 million reads were produced per sample and, for the purpose of gene expression quantification, only reads mapped to one location of the transcriptome were considered. After statistical analysis, we found 103 genes differentially expressed by schistosomula cultured for 3 hours and 12 hours in the presence of hamster portal serum. After the subtraction of a second list of genes, also differentially expressed between schistosomula cultured for 3 hours and 12 hours in presence of peripheral serum, a set of 58 genes was finally established. This pattern was further validated for a subset of 17 genes, by measuring gene expression through quantitative real time polymerase chain reaction (qPCR). Processes that were activated by the portal serum stimulus include response to stress, membrane transport, protein synthesis and folding/degradation, signaling, cytoskeleton arrangement, cell adhesion and nucleotide synthesis. Additionally, a smaller number of genes down-regulated under the same condition act on cholinergic signaling, inorganic cation and organic anion membrane transport, cell adhesion and cytoskeleton arrangement. Considering the role of these genes in triggering processes that allow the parasite to quickly adapt, escape the immune response of the host and start maturation into an adult worm after contact with the portal serum, this work may point to unexplored molecular targets for drug discovery and vaccine development against schistosomiasis.

Meis/UNC-62 isoform dependent regulation of CoupTF-II/UNC-55 and GABAergic motor neuron subtype differentiation.

Gene regulatory networks orchestrate the assembly of functionally related cells within a cellular network. Subtle differences often exist among functionally related cells within such networks. How differences are created among cells with similar functions has been difficult to determine due to the complexity of both the gene and the cellular networks. In Caenorhabditis elegans, the DD and VD motor neurons compose a cross-inhibitory, GABAergic network that coordinates dorsal and ventral muscle contractions during locomotion. The Pitx2 homologue, UNC-30, acts as a terminal selector gene to create similarities and the Coup-TFII homologue, UNC-55, is necessary for creating differences between the two motor neuron classes. What is the organizing gene regulatory network responsible for initiating the expression of UNC-55 and thus creating differences between the DD and VD motor neurons? We show that the unc-55 promoter has modules that contain Meis/UNC-62 binding sites. These sites can be subdivided into regions that are capable of activating or repressing UNC-55 expression in different motor neurons. Interestingly, different isoforms of UNC-62 are responsible for the activation and the stabilization of unc-55 transcription. Furthermore, specific isoforms of UNC-62 are required for proper synaptic patterning of the VD motor neurons. Isoform specific regulation of differentiating neurons is a relatively unexplored area of research and presents a mechanism for creating differences among functionally related cells within a network.

A mir-231-Regulated Protection Mechanism against the Toxicity of Graphene Oxide in Nematode Caenorhabditis elegans.

Recently, several dysregulated microRNAs (miRNAs) have been identified in organisms exposed to graphene oxide (GO). However, their biological functions and mechanisms of the action are still largely unknown. Here, we investigated the molecular mechanism of mir-231 in the regulation of GO toxicity using in vivo assay system of Caenorhabditis elegans. We found that GO exposure inhibited the expression of mir-231::GFP in multiple tissues, in particular in the intestine. mir-231 acted in intestine to regulate the GO toxicity, and overexpression of mir-231 in intestine caused a susceptible property of nematodes to GO toxicity. smk-1 encoding a homologue to mammalian SMEK functioned as a targeted gene for mir-231, and was also involved in the intestinal regulation of GO toxicity. Mutation of smk-1 gene induced a susceptible property to GO toxicity, whereas the intestinal overexpression of smk-1 resulted in a resistant property to GO toxicity. Moreover, mutation of smk-1 gene suppressed the resistant property of mir-231 mutant to GO toxicity. In nematodes, SMK-1 further acted upstream of the transcriptional factor DAF-16/FOXO in insulin signaling pathway to regulate GO toxicity. Therefore, mir-231 may encode a GO-responsive protection mechanism against the GO toxicity by suppressing the function of the SMK-1 - DAF-16 signaling cascade in nematodes.

DsRNA-mediated silencing of Nudix hydrolase in Trichinella spiralis inhibits the larval invasion and survival in mice.

The aim of this study was to investigate the functions of Trichinella spiralis Nudix hydrolase (TsNd) during the larval invasion of intestinal epithelial cells (IECs), development and survival in host by RNAi. The TsNd-specific double-stranded RNA (dsRNA) was designed to silence the expression of TsNd in T. spiralis larvae. DsRNA were delivered to the larvae by soaking incubation or electroporation. Silencing effect of TsNd transcription and expression was determined by real-time PCR and Western blotting, respectively. The infectivity of larvae treated with dsRNA was investigated by the in vitro larval invasion of IECs and experimental infection in mice. After being soaked with 40 ng/µl of dsRNA-TsNd, the transcription and expression level of TsNd gene was inhibited 65.8% and 56.4%, respectively. After being electroporated with 40 ng/µl of dsRNA-TsNd, the transcription and expression level of TsNd gene was inhibited 74.2% and 58.2%, respectively. Silencing TsNd expression by both soaking and electroporation inhibited significantly the larval invasion of IECs in a dose-dependent manner (r1 = -0.96798, r2 = -0.98707). Compared with the mice inoculated with untreated larvae, mice inoculated with larvae soaked with TsNd dsRNA displayed a 49.9% reduction in adult worms and 39.9% reduction in muscle larvae, while mice inoculated with larvae electroporated with TsNd dsRNA displayed a 83.4% reduction in adult worms and 69.5% reduction in muscle larvae, indicating that electroporation has a higher efficiency than soaking in inhibiting the larval development and survival in mice. Our results showed that silencing TsNd expression in T. spiralis inhibited significantly the larval invasion and survival in host.

Methylmercury exposure increases lipocalin related (lpr) and decreases activated in blocked unfolded protein response (abu) genes and specific miRNAs in Caenorhabditis elegans.

Methylmercury (MeHg) is a persistent environmental and dietary contaminant that causes serious adverse developmental and physiologic effects at multiple cellular levels. In order to understand more fully the consequences of MeHg exposure at the molecular level, we profiled gene and miRNA transcripts from the model organism Caenorhabditis elegans. Animals were exposed to MeHg (10 µM) from embryo to larval 4 (L4) stage and RNAs were isolated. RNA-seq analysis on the Illumina platform revealed 541 genes up- and 261 genes down-regulated at a cutoff of 2-fold change and false discovery rate-corrected significance q < 0.05. Among the up-regulated genes were those previously shown to increase under oxidative stress conditions including hsp-16.11 (2.5-fold), gst-35 (10.1-fold), and fmo-2 (58.5-fold). In addition, we observed up-regulation of 6 out of 7 lipocalin related (lpr) family genes and down regulation of 7 out of 15 activated in blocked unfolded protein response (abu) genes. Gene Ontology enrichment analysis highlighted the effect of genes related to development and organism growth. miRNA-seq analysis revealed 6-8 fold down regulation of mir-37-3p, mir-41-5p, mir-70-3p, and mir-75-3p. Our results demonstrate the effects of MeHg on specific transcripts encoding proteins in oxidative stress responses and in ER stress pathways. Pending confirmation of these transcript changes at protein levels, their association and dissociation characteristics with interaction partners, and integration of these signals, these findings indicate broad and dynamic mechanisms by which MeHg exerts its harmful effects.

CEL-Seq: single-cell RNA-Seq by multiplexed linear amplification.

High-throughput sequencing has allowed for unprecedented detail in gene expression analyses, yet its efficient application to single cells is challenged by the small starting amounts of RNA. We have developed CEL-Seq, a method for overcoming this limitation by barcoding and pooling samples before linearly amplifying mRNA with the use of one round of in vitro transcription. We show that CEL-Seq gives more reproducible, linear, and sensitive results than a PCR-based amplification method. We demonstrate the power of this method by studying early C. elegans embryonic development at single-cell resolution. Differential distribution of transcripts between sister cells is seen as early as the two-cell stage embryo, and zygotic expression in the somatic cell lineages is enriched for transcription factors. The robust transcriptome quantifications enabled by CEL-Seq will be useful for transcriptomic analyses of complex tissues containing populations of diverse cell types.

Temporal quantification of Ts43 gene expression of Trichinella spiralis using real-time RT-qPCR.

Determination of the transcription level of Trichinella spiralis (T. spiralis) Ts43 gene is essential for understanding its poorly explained role in different developmental stage. Herein we report on measurement of T. spiralis Ts43 mRNA by adapting a relative quantitative real-time reverse transcription and polymerase chain reaction (RT-PCR). Total RNA was isolated from the samples of adult worms, newborn larvae and muscle larvae of T. spiralis and Ts43 mRNA was quantified by real-time RT-PCR using an internally calibrated standard curve constructed with the recombinant Ts43 and 18S rRNA plasmid. The results indicated that the Ts43 gene was expressed at any stage and the expression levels in adults and newborn larvae were lower than those in muscle larvae in general. The expression of the Ts43 gene increased from 18 days post-infection (p.i.), reached a peak at 22-30 days p.i., and then decreased to a low level by 58 days p.i. but it was still higher than that in adults and newborn larvae. The study provided the first quantitative of Ts43 mRNA expression in T. spiralis for further studies of Ts43 function.

Small regulatory RNAs inhibit RNA polymerase II during the elongation phase of transcription.

Eukaryotic cells express a wide variety of endogenous small regulatory RNAs that regulate heterochromatin formation, developmental timing, defence against parasitic nucleic acids and genome rearrangement. Many small regulatory RNAs are thought to function in nuclei. For instance, in plants and fungi, short interfering RNA (siRNAs) associate with nascent transcripts and direct chromatin and/or DNA modifications. To understand further the biological roles of small regulatory RNAs, we conducted a genetic screen to identify factors required for RNA interference (RNAi) in Caenorhabditis elegans nuclei. Here we show that the gene nuclear RNAi defective-2 (nrde-2) encodes an evolutionarily conserved protein that is required for siRNA-mediated silencing in nuclei. NRDE-2 associates with the Argonaute protein NRDE-3 within nuclei and is recruited by NRDE-3/siRNA complexes to nascent transcripts that have been targeted by RNAi. We find that nuclear-localized siRNAs direct an NRDE-2-dependent silencing of pre-messenger RNAs (pre-mRNAs) 3' to sites of RNAi, an NRDE-2-dependent accumulation of RNA polymerase (RNAP) II at genomic loci targeted by RNAi, and NRDE-2-dependent decreases in RNAP II occupancy and RNAP II transcriptional activity 3' to sites of RNAi. These results define NRDE-2 as a component of the nuclear RNAi machinery and demonstrate that metazoan siRNAs can silence nuclear-localized RNAs co-transcriptionally. In addition, these results establish a novel mode of RNAP II regulation: siRNA-directed recruitment of NRDE factors that inhibit RNAP II during the elongation phase of transcription.

Additive effects of plant expressed double-stranded RNAs on root-knot nematode development.

Ectopically expressed double-stranded RNAs (dsRNAs) have recently been shown to suppress parasitic success of Meloidogyne spp. in plants. We have targeted two genes from the root-knot nematode Meloidogyne incognita; a dual oxidase gene implicated in the tyrosine cross-linking of the developing cuticle and a subunit of signal peptidase, a protein complex required for the processing of secreted proteins. While these genes are involved in different aspects of nematode development, the phenotypic consequences of RNA interference (RNAi) were similar with >or=50% reduction in nematode numbers in the roots and retardation of development to the egg-producing saccate females. Expression of processed dsRNA was observed, but no evidence of detectable levels of small interfering RNAs (siRNAs) was found in the transgenic plants. We show, to our knowledge for the first time, that combining expression of these dsRNAs by crossing appropriate Arabidopsis thaliana lines resulted in an additive effect that further reduced nematode numbers and developmental capacity. Combining RNAi target genes has the potential to enhance the efficacy of RNAi and may allow control of different nematode species or genera in the crop of interest.

Cloning, expression and evaluation of the efficacy of a recombinant Baylisascaris schroederi Bs-Ag3 antigen in mice.

The gene Bs-Ag3 enconding a antigen of 37kDa from Baylisascaris schroederi (giant panda isolates), as well as the recombinant Bs-Ag3, obtained by cloning and expression of the Bs-Ag3 gene in heterologous host Escherichia coli BL-21 (DE3), were used to evaluate their ability to induce immune protective responses in BALB/c mice against L3-challenge infection in a mouse-B. schroederi model. There was a significant reduction (62.91%) of mice vaccinated with rBs-Ag3 coupled with Freund's complete adjuvant (FCA) in recovery of challenged B. schroederi L3 compared with either controls or mice vaccinated with FCA alone. Our data indicate recombinant Bs-Ag3 may be a potential target as a vaccine antigen for giant panda ascariasis.

RETRACTED: Secondary siRNAs result from unprimed RNA synthesis and form a distinct class.

In Caenorhabditis elegans, an effective RNA interference (RNAi) response requires the production of secondary short interfering RNAs (siRNAs) by RNA-directed RNA polymerases (RdRPs). We cloned secondary siRNAs from transgenic C. elegans lines expressing a single 22-nucleotide primary siRNA. Several secondary siRNAs start a few nucleotides downstream of the primary siRNA, indicating that non-RISC (RNA-induced silencing complex)-cleaved mRNAs are substrates for secondary siRNA production. In lines expressing primary siRNAs with single-nucleotide mismatches, secondary siRNAs do not carry the mismatch but contain the nucleotide complementary to the mRNA. We infer that RdRPs perform unprimed RNA synthesis. Secondary siRNAs are only of antisense polarity, carry 5' di- or triphosphates, and are only in the minority associated with RDE-1, the RNAi-specific Argonaute protein. Therefore, secondary siRNAs represent a distinct class of small RNAs. Their biogenesis depends on RdRPs, and we propose that each secondary siRNA is an individual RdRP product.