(C) PLOS One [1]. This unaltered content originally appeared in journals.plosone.org. Licensed under Creative Commons Attribution (CC BY) license. url:https://journals.plos.org/plosone/s/licenses-and-copyright ------------ The Schistosoma mansoni nuclear receptor FTZ-F1 maintains esophageal gland function via transcriptional regulation of meg-8.3 ['Aracely A. Romero', 'Department Of Pharmacology', 'Ut Southwestern Medical Center', 'Dallas', 'Texas', 'United States Of America', 'Sarah A. Cobb', 'Julie N. R. Collins', 'Steven A. Kliewer', 'Department Of Molecular Biology'] Date: 2022-02 Schistosomes infect over 200 million of the world’s poorest people, but unfortunately treatment relies on a single drug. Nuclear hormone receptors are ligand-activated transcription factors that regulate diverse processes in metazoans, yet few have been functionally characterized in schistosomes. During a systematic analysis of nuclear receptor function, we found that an FTZ-F1-like receptor was essential for parasite survival. Using a combination of transcriptional profiling and chromatin immunoprecipitation (ChIP), we discovered that the micro-exon gene meg-8.3 is a transcriptional target of SmFTZ-F1. We found that both Smftz-f1 and meg-8.3 are required for esophageal gland maintenance as well as integrity of the worm’s head. Together, these studies define a new role for micro-exon gene function in the parasite and suggest that factors associated with the esophageal gland could represent viable therapeutic targets. Schistosomes kill an estimated 250,000 people every year and cause severe morbidity in millions more. Therefore, it is critical we continue to understand fundamental processes in the parasite so new therapies can be developed. Here, we identify a transcriptional regulator protein SmFTZ-F1 that is critical for parasite survival. We find that this protein controls the levels of a schistosome-specific gene (meg-8.3) and that meg-8.3 is critical for the maintenance of the worm’s esophageal gland. Loss of \ either meg-8.3 or Smftz-f1 leads to degeneration of the parasites’ head tissues. These studies suggest that molecules associated with esophageal gland function could represent therapeutic targets. Funding: This work was supported by a National Institutes of Health grants 5T32GM007062 to A.R., R01-AI150776 and R01-AI121037 to J.J.C.; Robert A. Welch Foundation grants I-1558 to S.A.K., I-1275 to D.J.M. and I-1948 to J.J.C.; and the Howard Hughes Medical Institute to D.J.M. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. NRs possess two distinct domains: a zinc finger DNA-binding domain and a ligand-binding domain. These domains, especially the DNA binding domain, are highly conserved across species and can be used to identify genes within the superfamily. In comparing known NR sequences from other metazoan species to schistosome transcriptomic data, we generated a list of 21 putative schistosome NRs. While multiple studies have been conducted on these receptors in S. mansoni [ 8 – 26 ], their function largely remains uncharacterized in schistosomes. We cloned full-length sequences for 17 out of 21 of these receptors and functionally characterized these receptors using RNA interference and in situ hybridizations. We found that mRNA depletion in S. mansoni of a previously described homolog of NR5A2 (SmFTZ-F1 (Smp_328000)) [ 12 , 22 ] leads to curling of the parasite’s head and eventually death in vitro. First identified in Drosophila, the FTZ-F1 receptor governs a number of early morphological processes, highlighting its role as a crucial developmental transcription factor [ 27 ]. RNA-seq of RNAi-treated parasites revealed transcriptional changes following loss of SmFTZ-F1 and led to the discovery of meg-8.3 as direct transcriptional target of SmFTZ-F1. meg-8.3 is expressed exclusively in the worm’s esophageal gland, an enigmatic tissue that has recently been shown to play a critical role in defending the worm from host-attack [ 28 ]. We find that loss of either Smftz-f1 or meg-8.3 results in defects in the maintenance of the esophageal gland and ultimately leads to degeneration of the worm’s head. Together these data suggest that disrupting esophageal gland function is catastrophic for the worm, highlighting the potential of targeting this tissue for the development of therapeutics. Currently, praziquantel is the only drug used to treat schistosomiasis and there is a growing need for novel therapeutics [ 2 ]. Unfortunately, few processes within these parasites have been shown to be critical for parasite survival, hindering target-based approaches for therapeutic development. Nuclear hormone receptors (NRs) are ligand-activated transcription factors that regulate diverse processes in metazoans [ 3 ], including homeostasis, reproduction, and cellular differentiation [ 4 – 7 ]. Given their diverse and conserved roles in metazoan biology, together with a paucity of molecular studies examining their function in schistosomes, we reasoned that examining this class of proteins could reveal biological processes critical to parasite survival. Schistosomiasis is a disease caused by parasitic flatworms of the genus Schistosoma that affects over 200 million people. The main cause of disease pathology is the egg production of the parasites, which can lay several hundred eggs per day while living in the host. The eggs that do not exit into the environment get swept into the circulatory system and lodge in the host’s liver or bladder, depending on the schistosome species, which then leads to inflammation that can disrupt organ function [ 1 ]. ChIP assays were performed using SimpleChIP Plus Enzymatic Chromatin IP Kit (Cell Signaling) according to the manufacturer’s instructions. Before harvesting, male worm heads were amputated and processed for chromatin immunoprecipitation. Briefly, worms were fixed with formaldehyde (0.75%) to cross-link DNA and proteins, DNA was digested using micrococcal nuclease, chromatin was sheared using a Biorupter (8 cycles of sonication: 30s each, 0.5 sec on, 0.5 sec off; amplitude = 30%,), chromatin was incubated with antibodies (SmFTZ-F1 at 5 μl, preimmune serum negative control at 5 μl), and immunoprecipitates were bound to protein G magnetic beads (30 μl). The protein-DNA cross-linking was reversed, DNA was purified, and enrichment of DNA sequences was detected using qPCR and primers listed in S1 Table . Data were normalized and analyzed using percent input analysis [ 37 ]. The parasites were fixed in methacarn (60 mL of methanol, 30 mL chloroform, and 10 mL glacial acetic acid) for 1 hr. Parasites were dehydrated in MeOH and stored at −20°C until processing. Worms were then processed similarly to colorimetric and fluorescence in situ hybridization following rehydration steps which were performed as previously described [ 23 , 29 , 31 ] with the following modifications. After post-fixation with 4% formaldehyde in PBSTx, the worms were then washed twice in PBSTx then were blocked with FISH block (0.1 M Tris pH 7.5, 0.15 M NaCl and 0.1% Tween-20 with 5% Horse Serum and 0.5% Roche Western Blocking Reagent [ 36 ]) for 1 hr and incubated overnight with 2.5 μg/mL of anti-SmFTZ-F1 antibody. After 6X PBSTx washes for 20 minutes each, the worms were then incubated in secondary antibody (1:1000 Goat anti-Rabbit IgG (H+L) Secondary Antibody, Alexa Fluor 488 conjugate) in FISH block for 4 hours at 4°C. After several washes in PBSTx, parasites were counterstained with DAPI for 30 minutes then placed in 80% glycerol and were then mounted on slides in Vectashield (Vector Laboratories). To generate protein lysates, RNAi-treated male parasites were separated with 0.5% tricaine, homogenized in 100 μl of sample buffer (2.5 mL of Solution D (0.47 M Tris, pH 6.7 (with 1 M phosphoric acid)) 1 mL of 100% glycerol, 1 mL 20% SDS, 0.5 mL H 2 O, and 50 μL of 1 M DTT) and protease inhibitors (Roche cOmplete, Mini, EDTA-free)). Homogenized samples were sonicated using a Biorupter system at 30% amplitude for 5 minutes, 0.5 sec on, 0.5 sec off. Protein concentrations were determined by Pierce Detergent Compatible Bradford Assay Kit. 40 mg of lysate was separated by SDS-PAGE, proteins were transferred to nitrocellulose membranes, membranes were blocked in Li-Cor Odyssey Blocking Buffer, incubated in rabbit anti-SmFTZ-F1 (0.5 mg/mL) and mouse anti-Actin (0.05mg/ml, Monoclonal JLA20, Developmental Studies Hybridoma Bank), washed in TBST, and incubated in secondary antibodies (1:10,000 goat anti-rabbit IR Dye 680 RD, 1:15,000 goat anti-mouse IgM IR Dye 800CW, Li-Cor). Blots were imaged using a Li-Cor Odyssey Infrared Imager. Biotin end-labeled probes were prepared using Pierce Biotin 3’ End DNA Labeling Kit following manufacturer’s instructions with the following modifications. Complementary oligos were end-labeled separately and then annealed before use. Oligos were annealed by mixing together equal amounts of labeled complementary oligos, denaturing the mixture at 90°C for 1 minute, then slowly cooling and incubating at the melting temperature for 30 minutes. The oligos were then frozen and thawed on ice for use. SmFTZ-F1 protein was prepared by in vitro translation with PURExpress (NEB). EMSA was performed using the LightShift Chemiluminescent EMSA Kit (ThermoFisher Scientific, USA) according to the manufacturer’s instructions. DNA binding reactions were performed in a 20 μL volume containing biotin-labeled oligonucleotides and in vitro translated protein from the PURExpress kit. Reaction products were then separated by electrophoresis. Thereafter, the protein-DNA complexes were transferred onto a positively charged nylon membrane (Millipore, USA) and detected by chemiluminescence using a GE Healthcare ImageQuant LAS 4000. For competitive binding experiments, 20- or 200-fold excesses of unlabeled DNA probes were also included in the binding reaction. To generate an anti-SmFTZ-F1 antibody, an N-terminal fragment of SmFTZ-F1 corresponding to AA 1–150 was amplified and subcloned into the pET28 vector with an N-terminal SUMO-His6 tag for expression in Escherichia coli. This fragment was purified from transformed BL21 cells grown in LB medium and induced with 1mM Isopropyl-β-D-thiogalactoside for 16 hrs at 18°C. Cells were pelleted and resuspended into lysis buffer containing 50 mM NaH 2 PO 4 , 500 mM NaCl, 20 mM Imidazole, and protease inhibitors (Roche cOmplete, Mini, EDTA-free). The suspension was sonicated, lysate was centrifuged for 25 min at 14,000 x g, and rotated with resin for 1 hr. The resin was washed, and protein was eluted in lysis buffer containing 250mM imidazole. The sample was then dialyzed to remove the NaCl and imidazole. The SUMO-His6tag was then cleaved using purified ULP-1. The protein was then concentrated using an Amicon concentrator, 10-kDa cut-off (Millipore). Rabbit polyclonal antibodies were generated by YenZym (California). To examine gene expression changes following loss of Smp_328000 (Smftz-f1), worms were cultured as previously described [ 33 ] and only male worms were harvested after 5 or 9 days in culture. We predicted that transcriptional changes that occur at earlier timepoints following RNAi are most likely to represent direct targets of SmFTZ-F1, whereas changes at later time points would represent both the direct and indirect effects of SmFTZ-F1 loss. As controls, worms cultured in parallel were treated with a non-specific dsRNA [ 31 ]. For RNA extraction, male worms were collected, excess media removed, and 100 μL of TRIzol was added to the worms. Total cellular RNA was isolated from cells using TRIzol reagent (Invitrogen) following the manufacturer’s protocol and processed for Illumina sequencing. Libraries for RNAseq analysis were mapped to the S. mansoni genome and analyzed as previously described [ 33 ]. RNAseq datasets for the Smftz-f1(RNAi) experiments are available at NCBI through the accession number (GSE188736). dsRNA was generated as previously described [ 30 , 31 ]. For both the initial NR screen and the RNA-seq “screen” RNAi experiments, RNAi experiments were performed as previously described using paired worms [ 32 ]. Subsequent experiments were performed using only male parasites. The parasites were fixed as previously described [ 31 ]. Throughout the 20-day experiment, worms were monitored for gross physical phenotypes (i.e., detachment from the plate, loss of movement, head curling, etc). Results/discussion An SmFTZ-F1 homolog is essential for parasite vitality To investigate potential roles that nuclear receptors play in schistosome biology, we set out to examine their expression and function using both whole mount in situ hybridization (WISH) and RNA interference (RNAi). We cloned either full or partial length cDNAs for 17 of 21 receptors from adult S. mansoni cDNA. We were unable to clone the remaining four due to low expression in the adult stage [8]. Expression analyses revealed that a majority of the NRs were expressed in the reproductive systems of both the males and females (S1 Fig). VF1, a parasite-specific NR [40], was expressed in cells of the female vitellaria (a tissue responsible for making “yolk” cells) as well as what appeared to be primordial vitellaria of male worms (S1 Fig). Although many NRs were broadly expressed throughout the parasite tissues (e.g., Smftz-f1), the expression of some NRs was restricted to smaller sets of cells (e.g., SmNR4A4 and SmCOUP-TF1). To explore the function of the schistosome NR complement, we performed RNAi for each NR and examined whether any of these treatments compromised parasite vitality in vitro (Fig 1A). We found that RNAi-mediated depletion of the mRNA encoding the schistosome ftz-f1 homolog [12,22] resulted in a fully penetrant degenerative phenotype (Fig 1B and 1C). By day 9 (D9), male and female Smftz-f1(RNAi) parasites became unpaired and lost the ability to attach to the surface of the tissue culture dish, and by D20, movement in Smftz-f1(RNAi) parasites was virtually nonexistent in male worms (S1 Video). The mRNA of Smftz-f1 was found to be expressed throughout both sexes of the worm (Figs 1D and S2A). Similar phenotypes were observed with dsRNAs targeting two distinct regions of the Smftz-f1 gene, indicating these effects are specific to the reduction of Smftz-f1 level and not due to off-target effects (S2B–S2D Fig). PPT PowerPoint slide PNG larger image TIFF original image Download: Fig 1. An SmFTZ-F1 homolog is essential for parasite vitality. (A) A schematic diagram depicting RNAi screen strategy B and C, Images of control and Smftz-f1(RNAi) worms during in vitro culture. At day 20 of in vitro culture, control(RNAi) animals remain as male and female worm pairs, are firmly attached to the bottom of the plate, and move freely. By D9 of in vitro culture, Smftz-f1(RNAi) animals become unpaired and fail to firmly attach to the bottom of the place. Over time, movement of Smftz-f1(RNAi) worms diminishes and they stop moving altogether. At D9, the heads are curled (arrow). At D20 the worm’s full body is curled. n = 10 worm pairs, > 3 biological replicates (D) Whole mount in situ hybridization showing mRNA expression of Smftz-f1 in both male and female parasites. The expression is found throughout the bodies of both sexes and does not seem to be specific for a tissue. Scale Bars: B and C, 500 μm. D,100 μm. https://doi.org/10.1371/journal.ppat.1010140.g001 Smftz-f1 and meg-8.3 share similar RNAi phenotypes We could envision two non-mutually exclusive scenarios for the phenotype observed following Smftz-f1 knockdown. In the first scenario, loss of Smftz-f1 disrupts the expression of multiple target genes and this broad transcriptional misregulation results in parasite degeneration and death. Alternatively, loss of Smftz-f1 leads to misregulation of individual target genes resulting in some or all of the observed phenotype. To distinguish between these models, we performed RNAi on each of the seven of the putative Smftz-f1 target genes. Consistent with the second scenario, we observed that RNAi of Smp_331380 led to detachment of parasites from the tissue culture dish beginning at D12 and a head curling phenotype identical to that observed following Smftz-f1(RNAi) (Fig 3B and 3C). Based on BLAST homology, Smp_331380 appeared to be a member of the schistosome-specific MEG-8 protein family, sharing complete sequence identity with the previously described MEG-8.3 protein [44,45]. Thus, we will refer to Smp_331380 as meg-8.3 from this point forward. MEG proteins are small proteins of largely unknown function that are encoded by small exons (~ 6 to 36 bp) that were uncovered during the sequencing of the S. mansoni genome [46,47]. The meg-8 protein family includes four paralogs (meg-8.1, 8.2, 8.3, and 8.4), whose homolog in Schistosoma japonicum MEG-8.2 appears to interact with host leukocytes as they pass by the worm’s esophageal gland [48]. The meg-8.3 gene is predicted to encode an 11 kD protein with a predicted signal peptide (S4A Fig) generated from 10 exons [49], some of which are as small as 15 bp in length (S4B Fig). In principle, the observed loss of meg-8.3 expression in the esophageal gland following Smftz-f1 loss could be due to either misregulation of meg-8.3 expression or loss of the esophageal gland tissue entirely. To explore this, we queried scRNAseq data [23] and examined the expression of a panel of genes expressed in the esophageal gland following Smftz-f1(RNAi) treatment. In contrast to meg-8.3 expression that was rapidly depleted following Smftz-f1 knockdown (Fig 3A), the expression of six esophageal gland marker genes were not obviously different in Smftz-f1(RNAi) parasites (Fig 3D). This suggests that loss of Smftz-f1 appears to rapidly and specifically lead to the reduction of meg-8.3 expression in the esophageal gland. We next evaluated the general morphology of the esophageal gland using the lectin PNA [50]. Under control conditions, the worms maintained a normal esophageal gland, however loss of either Smftz-f1 or meg-8.3 resulted in a loss of PNA staining in the esophageal gland with little labeling observed by D18 of RNAi (Fig 3E). Simultaneously, we noted a significant reduction in the expression of the esophageal gland marker, meg-4 between D15 and D18 of meg-8.3 RNAi treatment (S4C and S4D Fig). Taken together, these data are consistent with the model that SmFTZ-F1 regulates the expression of meg-8.3 in the esophageal gland, and that loss of meg-8.3 expression in the tissue results in loss of esophageal gland function that is at least partially responsible for the detachment and head curling observed in Smftz-f1(RNAi) parasites. Given that the S. mansoni genome contains three other MEG-8 proteins with sequence identity to MEG-8.3, we wanted to rule out the possibility that our phenotype was due to off-target effects. Therefore, we measured the levels of meg-8.3 and its closest relative meg-8.2 following meg-8.3 RNAi. As anticipated, RNAi of meg-8.3 resulted in large decreases in meg-8.3 mRNA levels (S4E Fig). Quantitative analysis of genes expressed in the esophageal gland following meg-8.3(RNAi) is complicated by the fact that loss of meg-8.3 appears to lead to loss of esophageal gland tissue over time (Figs 3E, S4C and S4D). Therefore, to examine the potential of RNAi off-target effects on meg-8.2 levels following meg-8.3 RNAi, we examined the expression of meg-8.2 and an unrelated esophageal gland marker gene cystatin (Smp_328980) [23]. When compared to meg-8.3 levels, we observed modest decreases in the levels of both meg-8.2 and cystatin following meg-8.3(RNAi) (S4E Fig), suggesting the reduction in transcript levels is due to tissue dysfunction rather than off target RNAi effects. This observation, along with the fact meg-8.2 and meg-8.3 share limited sequence identity on the nucleotide level [51], indicates that meg-8.3 RNAi-treatment specifically targets meg-8.3 levels. SmFTZ-F1 controls meg-8.3 expression To determine if meg-8.3 is indeed a direct target of SmFTZ-F1 in the parasite, we generated an antibody against SmFTZ-F1 and performed chromatin immune precipitation experiments. Our SmFTZ-F1 antiserum recognizes an ~110 kD protein in whole worm lysates (S5A Fig) and the levels of this protein are depleted upon Smftz-f1 dsRNA treatment (S5B Fig), confirming the specificity of the antibody. As expected for a DNA-binding transcription factor, we observed by immunofluorescence that this antibody labeled nuclei throughout the worm body (S5C Fig). This signal was also lost following Smftz-f1(RNAi) treatment (S5D Fig). Given the specificity of the anti-SmFTZ-F1 antibody, we set out to determine if meg-8.3 is a direct target of SmFTZ-F1 in vivo by ChIP. For these studies, we amputated the heads of male worms to enrich for cells expressing meg-8.3, performed ChIP using either anti-SmFTZ-F1 or pre-immune serum, and then used quantitative PCR to evaluate relative SmFTZ-F1 binding to various sites in the genome (Fig 3F and 3G). As an initial step, we performed various controls evaluating SmFTZ-F1 binding at distal regions 5 and 10 kb upstream of meg-8.3 and at a region surrounding a SmFTZ-F1 consensus response element found upstream of foxA, a gene encoding a transcription factor that appears to be a master regulator of esophageal gland identity [28]; we did not note significant SmFTZ-F1 binding to any of these sites (S6A–S6D Fig). We next evaluated the consensus SmFTZ-F1 response element upstream of meg-8.3 that was bound by SmFTZ-F1 in vitro (Fig 2D). Evaluation of the sequence surrounding this consensus site revealed that it was repeated at several sites within the S. mansoni genome (S1 File), ruling out possibility of designing primers that could specifically amplify this region of the genome. However, since a large number of transcription factor binding sites are known to exist in repetitive elements [52–54], we proceeded with ChIP studies. Using qPCR primers that amplify the region surrounding the putative SmFTZ-F1 binding site upstream of meg-8.3, we noted a statistically-significant 10-fold increase when we performed ChIP using the anti-SmFTZ-F1 (Fig 3H). ChIP-qPCR using Smftz-f1 RNAi-treated worms showed that the enrichment of SmFTZ-F1 at the response element in the meg-8.3 promoter was abolished to levels similar to those observed using preimmune serum (S6E Fig). These data, along with our RNAseq studies (Fig 2B) and EMSA data (Fig 2D), strongly support the model that SmFTZ-F1 binds to a response element found in the meg-8.3 promoter. [END] [1] Url: https://journals.plos.org/plospathogens/article?id=10.1371/journal.ppat.1010140 (C) Plos One. "Accelerating the publication of peer-reviewed science." Licensed under Creative Commons Attribution (CC BY 4.0) URL: https://creativecommons.org/licenses/by/4.0/ via Magical.Fish Gopher News Feeds: gopher://magical.fish/1/feeds/news/plosone/