Determining Toxic Potencies of Water-Soluble Contaminants in Wastewater Influents and Effluent Using Gene Expression Profiling in C. elegans as a Bioanalytical Tool

With chemical analysis, it is impossible to qualify and quantify the toxic potency of especially hydrophilic bioactive contaminants. In this study, we applied the nematode C. elegans as a model organism for detecting the toxic potency of whole influent wastewater samples. Gene expression in the nematode was used as bioanalytical tool to reveal the presence, type and potency of molecular pathways induced by 24-h exposure to wastewater from a hospital (H), nursing home (N), community (C), and influent (I) and treated effluent (E) from a local wastewater treatment plant. Exposure to influent water significantly altered expression of 464 genes, while only two genes were differentially expressed in nematodes treated with effluent. This indicates a significant decrease in bioactive pollutant-load after wastewater treatment. Surface water receiving the effluent did not induce any genes in exposed nematodes. A subset of 209 genes was differentially expressed in all untreated wastewaters, including cytochromes P450 and C-type lectins related to the nematode's xenobiotic metabolism and immune response, respectively. Different subsets of genes responded to particular waste streams making them candidates to fingerprint-specific wastewater sources. This study shows that gene expression profiling in C. elegans can be used for mechanism-based identification of hydrophilic bioactive compounds and fingerprinting of specific wastewaters. More comprehensive than with chemical analysis, it can demonstrate the effective overall removal of bioactive compounds through wastewater treatment. This bioanalytical tool can also be applied in the process of identification of the bioactive compounds via a process of toxicity identification evaluation. [GRAPHICS]

Genetic Variation in Complex Traits in Transgenic α-Synuclein Strains of Caenorhabditis elegans.

Different genetic backgrounds can modify the effect of mutated genes. Human α-synuclein (SNCA) gene encodes α-synuclein, and its oligomeric complexes accumulate with age and mediate the disruption of cellular homeostasis, resulting in the neuronal death that is characteristic of Parkinson's Disease. Polymorphic variants modulate this complex pathologic mechanism. Previously, we constructed five transgenic introgression lines of a Caenorhabditis elegans model of α-synuclein using genetic backgrounds that are genetically diverse from the canonical wild-type Bristol N2. A gene expression analysis revealed that the α-synuclein transgene differentially affects genome-wide transcription due to background modifiers. To further investigate how complex traits are affected in these transgenic lines, we measured the α-synuclein transgene expression, the overall accumulation of the fusion protein of α-synuclein and yellow fluorescent protein (YFP), the lysosome-related organelles, and the body size. By using quantitative PCR (qPCR), we demonstrated stable and similar expression levels of the α-synuclein transgene in different genetic backgrounds. Strikingly, we observed that the levels of the a-synuclein:YFP fusion protein vary in different genetic backgrounds by using the COPAS™ biosorter. The quantification of the Nile Red staining assay demonstrates that α-synuclein also affects lysosome-related organelles and body size. Our results show that the same α-synuclein introgression in different C. elegans backgrounds can produces differing effects on complex traits due to background modifiers

Comparative genomics of two inbred lines of the potato cyst nematode Globodera rostochiensis reveals disparate effector family-specific diversification patterns

This work was funded as part of a grant by the Netherlands Organization for Scientific Research (NWO) as part of the Applied and Technical Science domain (TTW) under grant no. 14708. PT: bioinformatics and computational biology analyses were supported by the University of St Andrews Bioinformatics Unit (AMD3BIOINF), funded by Wellcome Trust ISSF award 105621/Z/14/Z.Background: Potato cyst nematodes belong to the most harmful pathogens in potato, and durable management of these parasites largely depends on host-plant resistances. These resistances are pathotype specific. The current Globodera rostochiensis pathotype scheme that defines five pathotypes (Ro1 - Ro5) is both fundamentally and practically of limited value. Hence, resistant potato varieties are used worldwide in a poorly informed manner. Results : We generated two novel reference genomes of G. rostochiensis inbred lines derived from a Ro1 and a Ro5 population. These genome sequences comprise 173 and 189 scaffolds respectively, marking a ≈ 24-fold reduction in fragmentation as compared to the current reference genome. We provide copy number variations for 19 effector families. Four dorsal gland effector families were investigated in more detail. SPRYSECs, known to be implicated in plant defence suppression, constitute by far the most diversified family studied herein with 60 and 99 variants in Ro1 and Ro5 distributed over 18 and 26 scaffolds. In contrast, CLEs, effectors involved in feeding site induction, show strong physical clustering. The 10 and 16 variants cluster on respectively 2 and 1 scaffolds. Given that pathotypes are defined by their effectoromes, we pinpoint the disparate nature of the contributing effector families in terms of sequence diversification and loss and gain of variants. Conclusions : Two novel reference genomes allow for nearly complete inventories of effector diversification and physical organisation within and between pathotypes. Combined with insights we provide on effector family-specific diversification patterns, this constitutes a basis for an effectorome-based virulence scheme for this notorious pathogen.Publisher PDFPeer reviewe

Transcriptional analysis of the response of \u3ci\u3eC. elegans\u3c/i\u3e to ethanol exposure

Ethanol-induced transcriptional changes underlie important physiological responses to ethanol that are likely to contribute to the addictive properties of the drug. We examined the transcriptional responses of Caenorhabditis elegans across a timecourse of ethanol exposure, between 30 min and 8 h, to determine what genes and genetic pathways are regulated in response to ethanol in this model. We found that short exposures to ethanol (up to 2 h) induced expression of metabolic enzymes involved in metabolizing ethanol and retinol, while longer exposure (8 h) had much more profound effects on the transcriptome. Several genes that are known to be involved in the physiological response to ethanol, including direct ethanol targets, were regulated at 8 h of exposure. This longer exposure to ethanol also resulted in the regulation of genes involved in cilia function, which is consistent with an important role for the effects of ethanol on cilia in the deleterious effects of chronic ethanol consumption in humans. Finally, we found that food deprivation for an 8-h period induced gene expression changes that were somewhat ameliorated by the presence of ethanol, supporting previous observations that worms can use ethanol as a calorie source

Chemical differentiation of Bolivian Cedrela species as a tool to trace illegal timber trade

Combating illegal timber trade requires the ability to identify species and verify geographic origin of timber. Forensic techniques that independently verify the declared species and geographic origin are needed, as current legality procedures are based on certificates and documents that can be falsified. Timber from the genus Cedrela is among the most economically valued tropical timbers worldwide. Three Cedrela species are included in the Appendix III of CITES: C. fissilis, C. odorata and C. angustifolia (listed as C. lilloi). Cedrela timber is currently traded with false origin declarations and under a different species name, but tools to verify this are lacking. We used Direct Analysis in Real Time Time-of-Flight Mass Spectrometry (DART-TOFMS) to chemically identify Cedrela species and sites of origin. Heartwood samples from six Cedrela species (the three CITES-listed species plus C. balansae, C. montana and C. saltensis) were collected at 11 sites throughout Bolivia. Mass spectra detected by DART-TOFMS comprised 1062 compounds; their relative intensities were analysed using Principal Component Analyses, Kernel Discriminant Analysis (KDA) and Random Forest analyses to check discrimination potential among species and sites. Species were identified with a mean discrimination error of 15–19 per cent, with substantial variation in discrimination accuracy among species. The lowest error was observed in C. fissilis (mean = 4.4 per cent). Site discrimination error was considerably higher: 43–54 per cent for C. fissilis and 42–48 per cent for C. odorata. These results provide good prospects to differentiate C. fissilis from other species, but at present there is no scope to do so for other tested species. Thus, discrimination is highly species specific. Our findings for tests of geographic origin suggest no potential to discriminate at the studied scale and for the studied species. Cross-checking results from different methods (KDA and Random Forest) reduced discrimination errors. In all, the DART-TOFMS technique allows independent verification of claimed identity of certain Cedrela species in timber trade

Comparative genomics among cyst nematodes reveals distinct evolutionary histories among effector families and an irregular distribution of effector-associated promoter motifs

JvS, MH and SvdE were supported by a grant from the Applied and Technical Science domain (TTW) of the Netherlands Organization for Scientific Research (NWO) under grant no. 14708. PT received support from the University of St Andrews Bioinformatics Unit (AMD3BIOINF), funded by Wellcome Trust ISSF award 105621/Z/14/Z. MS benefitted from funding by a VENI grant (17282) from the NWO domain Applied and Engineering Sciences.Potato cyst nematodes (PCNs), an umbrella term used for two species, Globodera pallida and G. rostochiensis, belong worldwide to the most harmful pathogens of potato. Pathotype-specific host plant resistances are an essential handle for PCN control. However, the poor delineation of G. pallida pathotypes hampers the efficient use of available host plant resistances. Long-read sequencing technology allowed us to generate a new reference genome of G. pallida population D383 and, as compared to the current reference, the new genome assembly is 42 times less fragmented. For comparison of diversification patterns of six effector families between G. pallida and G. rostochiensis, an additional reference genome was generated for an outgroup, the beet cyst nematode Heterodera schachtii (IRS population). Large evolutionary contrasts in effector family topologies were observed. While VAPs diversified before the split between the three cyst nematode species, the families GLAND5 and GLAND13 only expanded in PCN after their separation from the genus Heterodera. Although DNA motifs in the promoter regions thought to be involved in the orchestration of effector expression ('DOG boxes') were present in all three cyst nematode species, their presence is not a necessity for dorsal gland-produced effectors. Notably, DOG box dosage was only loosely correlated with expression level of individual effector variants. Comparison of the G. pallida genome with those of two other cyst nematodes underlined the fundamental differences in evolutionary history between effector families. Re-sequencing of PCN populations with deviant virulence characteristics will allow for the linking of these characteristics with the composition of the effector repertoire as well as for the mapping of PCN diversification patterns resulting from extreme anthropogenic range expansion.Publisher PDFPeer reviewe

The TIR-NB-LRR pair DSC1 and WRKY19 contributes to basal immunity of Arabidopsis to the root-knot nematode Meloidogyne incognita

BackgroundRoot-knot nematodes transform vascular host cells into permanent feeding structures to withdraw nutrients from the host plant. Ecotypes of Arabidopsis thaliana can display large quantitative variation in susceptibility to the root-knot nematode Meloidogyne incognita, which is thought to be independent of dominant major resistance genes. However, in an earlier genome-wide association study of the interaction between Arabidopsis and M. incognita we identified a quantitative trait locus harboring homologs of dominant resistance genes but with minor effect on susceptibility to the M. incognita population tested.ResultsHere, we report on the characterization of two of these genes encoding the TIR-NB-LRR immune receptor DSC1 (DOMINANT SUPPRESSOR OF Camta 3 NUMBER 1) and the TIR-NB-LRR-WRKY-MAPx protein WRKY19 in nematode-infected Arabidopsis roots. Nematode infection studies and whole transcriptome analyses using the Arabidopsis mutants showed that DSC1 and WRKY19 co-regulate susceptibility of Arabidopsis to M. incognita.ConclusionGiven the head-to-head orientation of DSC1 and WRKY19 in the Arabidopsis genome our data suggests that both genes may function as a TIR-NB-LRR immune receptor pair. Unlike other TIR-NB-LRR pairs involved in dominant disease resistance in plants, DSC1 and WRKY19 most likely regulate basal levels of immunity to root-knot nematodes

Cryptic genetic variation of expression quantitative trait locus architecture revealed by genetic perturbation in Caenorhabditis elegans

Genetic perturbation in different genetic backgrounds can cause a range of phenotypes within a species. These phenotypic differences can be the result of the interaction between the genetic background and the perturbation. Previously, we reported that perturbation of gld-1, an important player in the developmental control of Caenorhabditis elegans, released cryptic genetic variation (CGV) affecting fitness in different genetic backgrounds. Here, we investigated the change in transcriptional architecture. We found 414 genes with a cis-expression quantitative trait locus (eQTL) and 991 genes with a trans-eQTL that were specifically found in the gld-1 RNAi treatment. In total, we detected 16 eQTL hotspots, of which 7 were only found in the gld-1 RNAi treatment. Enrichment analysis of those 7 hotspots showed that the regulated genes were associated with neurons and the pharynx. Furthermore, we found evidence of accelerated transcriptional aging in the gld-1 RNAi-treated nematodes. Overall, our results illustrate that studying CGV leads to the discovery of hidden polymorphic regulators

Contribution of trans regulatory eQTL to cryptic genetic variation in C-elegans

BackgroundCryptic genetic variation (CGV) is the hidden genetic variation that can be unlocked by perturbing normal conditions. CGV can drive the emergence of novel complex phenotypes through changes in gene expression. Although our theoretical understanding of CGV has thoroughly increased over the past decade, insight into polymorphic gene expression regulation underlying CGV is scarce. Here we investigated the transcriptional architecture of CGV in response to rapid temperature changes in the nematode Caenorhabditis elegans. We analyzed regulatory variation in gene expression (and mapped eQTL) across the course of a heat stress and recovery response in a recombinant inbred population.ResultsWe measured gene expression over three temperature treatments: i) control, ii) heat stress, and iii) recovery from heat stress. Compared to control, exposure to heat stress affected the transcription of 3305 genes, whereas 942 were affected in recovering animals. These affected genes were mainly involved in metabolism and reproduction. The gene expression pattern in recovering animals resembled both the control and the heat-stress treatment. We mapped eQTL using the genetic variation of the recombinant inbred population and detected 2626 genes with an eQTL in the heat-stress treatment, 1797 in the control, and 1880 in the recovery. The cis-eQTL were highly shared across treatments. A considerable fraction of the trans-eQTL (40–57%) mapped to 19 treatment specific trans-bands. In contrast to cis-eQTL, trans-eQTL were highly environment specific and thus cryptic. Approximately 67% of the trans-eQTL were only induced in a single treatment, with heat-stress showing the most unique trans-eQTL.ConclusionsThese results illustrate the highly dynamic pattern of CGV across three different environmental conditions that can be evoked by a stress response over a relatively short time-span (2 h) and that CGV is mainly determined by response related trans regulatory eQTL.<br/

Transcriptomic and proteomic analysis of arbovirus-infected tick cells

Ticks are important vectors of a wide variety of pathogens including protozoa, bacteria and viruses. Many of the viruses transmitted by ticks are of medical or veterinary importance including tick-borne encephalitis virus (TBEV) and Crimean- Congo hemorrhagic fever virus causing disease in humans, and African swine fever virus and Nairobi sheep disease virus affecting livestock. Although several studies have elucidated tick antimicrobial mechanisms including cellular immune responses such as nodulation, encapsulation and phagocytosis and humoral immune responses such as the JAK/STAT pathway, complement-like proteins, antimicrobial peptides, lectin like pattern-recognition molecules and lysozymes, very little is known about the innate immune response of ticks towards viral infection. This study therefore aimed to identify molecules that might be involved in the response of ticks to viral infection. The hypothesis was that TBEV infection leads to changes in the expression of immunity-related transcripts and proteins in Ixodes spp. tick cells and that at least some of these might be antiviral. Ixodes scapularis-derived cell lines IDE8 and ISE6 were chosen since I. scapularis is currently the only tick species with a sequenced genome and an Ixodes ricinus-derived cell line, IRE/CTVM19, was used because I. ricinus is the natural vector of TBEV. Basic parameters required to study the responses of tick cells to infection were determined, including levels of virus infection, kinetics of virus replication and production, formation of replication complexes and uptake of dsRNA or siRNA. The cell lines IDE8, ISE6 and IRE/CTVM19 were infected with either of two tick-borne flaviviruses, TBEV and Langat virus (LGTV), or with the mosquito-borne alphavirus Semliki Forest virus (SFV). Infection was characterised using techniques including plaque assay, luciferase assay, immunostaining and conventional, confocal and electron microscopy. Two time points for transcriptomics and proteomics analysis of TBEVinfected IDE8 and IRE/CTVM19 cells were selected: day 2 post-infection (p.i.) when virus production was increasing and day 6 p.i. when virus production was decreasing. RNA and protein were isolated from TBEV-infected and mock-infected tick cells at days 2 and 6 p.i. and RNA-Seq and mass spectrometric technologies were used to identify changes in, respectively, transcript and protein abundance. Differential expression of transcripts was determined using the data analysis package DESeq resulting in a total of 43 statistically significantly differentially expressed transcripts in IDE8 cells and 83 in IRE/CTVM19 cells, while differential protein representation using Χ2 test statistics with Bonferroni correction in IDEG6 software resulted in 76 differentially represented proteins in IDE8 cells and 129 in IRE/CTVM19 cells. These included transcripts and proteins which could affect stages of the virus infection, including virus entry, replication, maturation and protein trafficking, and also innate immune responses such as phagocytosis, RNA interference (RNAi), the complement system, the ubiquitin-proteasome pathway, cell stress and the endoplasmic reticulum (ER) stress response. After verification of sequencing data by qRT-PCR, the ability of several of the identified transcripts or proteins to affect virus infection was determined by knockdown experiments in IDE8 and IRE/CTVM19 cells using wild type LGTV, LGTV replicons or TBEV replicons. Knockdown of genes encoding proteins including the ER chaperone gp96 and the heat-shock protein HSP90 resulted in increased virus production in both cell lines, hinting at an antiviral role. In contrast, knockdown of calreticulin, another ER chaperone, resulted in a decrease in virus production in IRE/CTVM19 cells but not in IDE8 cells, implying a requirement for virus production. This functional genomics approach has identified possible novel genes/proteins involved in the interaction between flaviviruses and tick cells and also revealed that there might be antiviral innate immune pathways present in ticks additional to the exogenous RNAi pathway