A Krüppel-like factor establishes cellular heterogeneity during schistosome tegumental maintenance.

Schistosomes are blood-dwelling parasitic flatworms that rely on a syncytial surface coat, known as the tegument, for long-term survival and immune evasion in the blood of their human hosts. Previous studies have shown that cells within the tegumental syncytium are perpetually turned over and renewed by somatic stem cells called neoblasts. Yet, little is known about this renewal process on a molecular level. Here, we characterized a Krüppel-like factor 4 (klf4) using a combination of bulk and single cell RNAseq approaches and demonstrate that klf4 is essential for the maintenance of a specific tegumental lineage, resulting in the loss of a subpopulation of molecularly-unique tegument cells. Thus, klf4 is critical for maintaining the balance between different tegumental progenitor pools, thereby fine-tuning the molecular composition of the mature tegument. Understanding these distinct tegumental cell populations is expected to provide insights into parasite defense mechanisms and suggest new avenues for therapeutics.

Horizontal gene transfer of a functional cki homolog in the human pathogen Schistosoma mansoni.

Schistosomes are parasitic flatworms responsible for the neglected tropical disease schistosomiasis, causing devastating morbidity and mortality in the developing world. The parasites are protected by a skin-like tegument, and maintenance of this tegument is controlled by a schistosome ortholog of the tumor suppressor TP53. To understand mechanistically how p53-1 controls tegument production, we identified a cyclin dependent kinase inhibitor homolog (cki) that was co-expressed with p53-1. RNA interference of cki resulted in a hyperproliferation phenotype, that, in combination with p53-1 RNA interference yielded abundant tumor-like growths, indicating that cki and p53-1 are bona fide tumor suppressors in Schistosoma mansoni. Interestingly, cki homologs are widely present throughout parasitic flatworms but evidently absent from their free-living ancestors, suggesting this cki homolog came from an ancient horizontal gene transfer event. This in turn implies that the evolution of parasitism in flatworms may have been aided by a highly unusual means of metazoan genetic inheritance.

Genetic predisposition to altered blood cell homeostasis is associated with glioma risk and survival.

Glioma is a highly fatal brain tumor comprised of molecular subtypes with distinct clinical trajectories. Observational studies have suggested that variability in immune response may play a role in glioma etiology. However, their findings have been inconsistent and susceptible to reverse causation due to treatment effects and the immunosuppressive nature of glioma. We applied genetic variants associated (p<5×10-8) with blood cell traits to a meta-analysis of 3418 glioma cases and 8156 controls. Genetically predicted increase in the platelet to lymphocyte ratio (PLR) was associated with an increased risk of glioma (odds ratio (OR)=1.25, p=0.005), especially in IDH-mutant (IDHmut OR=1.38, p=0.007) and IDHmut 1p/19q non-codeleted (IDHmut-noncodel OR=1.53, p=0.004) tumors. However, reduced glioma risk was observed for higher counts of lymphocytes (IDHmut-noncodel OR=0.70, p=0.004) and neutrophils (IDHmut OR=0.69, p=0.019; IDHmut-noncodel OR=0.60, p=0.009), which may reflect genetic predisposition to enhanced immune-surveillance. In contrast to susceptibility, there was no association with survival in IDHmut-noncodel; however, in IDHmut 1p/19q co-deleted tumors, we observed higher mortality with increasing genetically predicted counts of lymphocytes (hazard ratio (HR)=1.65, 95% CI: 1.24-2.20), neutrophils (HR=1.49, 1.13-1.97), and eosinophils (HR=1.59, 1.18-2.14). Polygenic scores for blood cell traits were also associated with tumor immune microenvironment features, with heterogeneity by IDH status observed for 17 signatures related to interferon signaling, PD-1 expression, and T-cell/Cytotoxic responses. In summary, we identified novel, immune-mediated susceptibility mechanisms for glioma with potential disease management implications.

Convergent evolution of a genotoxic stress response in a parasite-specific p53 homolog.

P53 is a widely studied tumor suppressor that plays important roles in cell-cycle regulation, cell death, and DNA damage repair. P53 is found throughout metazoans, even in invertebrates that do not develop malignancies. The prevailing theory for why these invertebrates possess a tumor suppressor is that P53 originally evolved to protect the germline of early metazoans from genotoxic stress such as ultraviolet radiation. This theory is largely based upon functional data from only three invertebrates, omitting important groups of animals including flatworms. Previous studies in the freshwater planarian flatworm Schmidtea mediterranea suggested that flatworm P53 plays an important role in stem cell maintenance and skin production, but these studies did not directly test for any tumor suppressor functions. To better understand the function of P53 homologs across diverse flatworms, we examined the function of two different P53 homologs in the parasitic flatworm Schistosoma mansoni. The first P53 homolog (p53-1) is orthologous to S. mediterranea P53(Smed-p53) and human TP53 and regulates flatworm stem cell maintenance and skin production. The second P53 homolog (p53-2) is a parasite-specific paralog that is conserved across parasitic flatworms and is required for the normal response to genotoxic stress in S. mansoni. We then found that Smed-p53 does not seem to play any role in the planarian response to genotoxic stress. The existence of this parasite-specific paralog that bears a tumor suppressor-like function in parasitic flatworms implies that the ability to respond to genotoxic stress in parasitic flatworms may have arisen from convergent evolution.

Schistosome Sulfotransferases: Mode of Action, Expression and Localization.

Oxamniquine (OXA) is a prodrug activated by a sulfotransferase (SULT) that was only active against Schistosoma mansoni. We have reengineered OXA to be effective against S. haematobium and S. japonicum. Three derivatives stand out, CIDD-0066790, CIDD-0072229, and CIDD-0149830 as they kill all three major human schistosome species. However, questions remain. Is the OXA mode of action conserved in derivatives? RNA-interference experiments demonstrate that knockdown of the SmSULT, ShSULT, and SjSULT results in resistance to CIDD-0066790. Confirming that the OXA-derivative mode of action is conserved. Next is the level of expression of the schistosome SULTs in each species, as well as changes in SULT expression throughout development in S. mansoni. Using multiple tools, our data show that SmSULT has higher expression compared to ShSULT and SjSULT. Third, is the localization of SULT in the adult, multicellular eucaryotic schistosome species. We utilized fluorescence in situ hybridization and uptake of radiolabeled OXA to determine that multiple cell types throughout the adult schistosome worm express SULT. Thus, we hypothesize the ability of many cells to express the sulfotransferase accounts for the ability of the OXA derivatives to kill adult worms. Our studies demonstrate that the OXA derivatives are able to kill all three human schistosome species and thus will be a useful complement to PZQ.

The immunogenetics of viral antigen response is associated with subtype-specific glioma risk and survival.

Glioma is a highly fatal cancer with prognostically significant molecular subtypes and few known risk factors. Multiple studies have implicated infections in glioma susceptibility, but evidence remains inconsistent. Genetic variants in the human leukocyte antigen (HLA) region modulate host response to infection and have been linked to glioma risk. In this study, we leveraged genetic predictors of antibody response to 12 viral antigens to investigate the relationship with glioma risk and survival. Genetic reactivity scores (GRSs) for each antigen were derived from genome-wide-significant (p < 5 × 10-8) variants associated with immunoglobulin G antibody response in the UK Biobank cohort. We conducted parallel analyses of glioma risk and survival for each GRS and HLA alleles imputed at two-field resolution by using data from 3,418 glioma-affected individuals subtyped by somatic mutations and 8,156 controls. Genetic reactivity scores to Epstein-Barr virus (EBV) ZEBRA and EBNA antigens and Merkel cell polyomavirus (MCV) VP1 antigen were associated with glioma risk and survival (Bonferroni-corrected p < 0.01). GRSZEBRA and GRSMCV were associated in opposite directions with risk of IDH wild-type gliomas (ORZEBRA = 0.91, p = 0.0099/ORMCV = 1.11, p = 0.0054). GRSEBNA was associated with both increased risk for IDH mutated gliomas (OR = 1.09, p = 0.040) and improved survival (HR = 0.86, p = 0.010). HLA-DQA1∗03:01 was significantly associated with decreased risk of glioma overall (OR = 0.85, p = 3.96 × 10-4) after multiple testing adjustment. This systematic investigation of the role of genetic determinants of viral antigen reactivity in glioma risk and survival provides insight into complex immunogenomic mechanisms of glioma pathogenesis. These results may inform applications of antiviral-based therapies in glioma treatment.

Analysis of Schistosoma mansoni Extracellular Vesicles Surface Glycans Reveals Potential Immune Evasion Mechanism and New Insights on Their Origins of Biogenesis.

Parasitic helminths are master manipulators of host immunity. Their strategy is complex and involves the release of excreted/secreted products, including extracellular vesicles (EVs). The protein and miRNA contents of EVs have been characterised for many parasitic helminths but, despite reports suggesting the importance of EV surface carbohydrate structures (glycans) in the interactions with target cells and thus subsequent effector functions, little is known about parasite EV glycomics. Using lectin microarrays, we identified several lectins that exhibit strong adhesion to Schistosoma mansoni EVs, suggesting the presence of multiple glycan structures on these vesicles. Interestingly, SNA-I, a lectin that recognises structures with terminal sialic acid, displayed strong affinity for S. mansoni EVs, which was completely abolished by neuraminidase treatment, suggesting sialylation in the EV sample. This finding is of interest, as sialic acids play important roles in the context of infection by aiding immune evasion, affecting target recognition, cell entry, etc., but are not thought to be synthesised by helminths. These data were validated by quantitative analysis of free sialic acid released from EVs following treatment with neuraminidase. Lectin histochemistry and fluorescence in situ hybridisation analyses on whole adult worms suggest the involvement of sub-tegumental cell bodies, as well as the digestive and excretory systems, in the release of EVs. These results support previous reports of EV biogenesis diversity in trematodes and potentially highlight new means of immune modulation and evasion employed by schistosomes.

SchistoCyte Atlas: A Single-Cell Transcriptome Resource for Adult Schistosomes.

Schistosomes cause untold disease and disability in the developing world. Here, we introduce SchistoCyte Atlas, a web-based platform for exploring gene expression at single-cell resolution in adult Schistosoma mansoni. Similar resources accessible to non-specialists across the globe will expedite our ability to understand the biology of these devastating parasites.

The landscape of host genetic factors involved in immune response to common viral infections.

BACKGROUND: Humans and viruses have co-evolved for millennia resulting in a complex host genetic architecture. Understanding the genetic mechanisms of immune response to viral infection provides insight into disease etiology and therapeutic opportunities. METHODS: We conducted a comprehensive study including genome-wide and transcriptome-wide association analyses to identify genetic loci associated with immunoglobulin G antibody response to 28 antigens for 16 viruses using serological data from 7924 European ancestry participants in the UK Biobank cohort. RESULTS: Signals in human leukocyte antigen (HLA) class II region dominated the landscape of viral antibody response, with 40 independent loci and 14 independent classical alleles, 7 of which exhibited pleiotropic effects across viral families. We identified specific amino acid (AA) residues that are associated with seroreactivity, the strongest associations presented in a range of AA positions within DRß1 at positions 11, 13, 71, and 74 for Epstein-Barr virus (EBV), Varicella zoster virus (VZV), human herpesvirus 7, (HHV7), and Merkel cell polyomavirus (MCV). Genome-wide association analyses discovered 7 novel genetic loci outside the HLA associated with viral antibody response (P < 5.0 × 10-8), including FUT2 (19q13.33) for human polyomavirus BK (BKV), STING1 (5q31.2) for MCV, and CXCR5 (11q23.3) and TBKBP1 (17q21.32) for HHV7. Transcriptome-wide association analyses identified 114 genes associated with response to viral infection, 12 outside of the HLA region, including ECSCR: P = 5.0 × 10-15 (MCV), NTN5: P = 1.1 × 10-9 (BKV), and P2RY13: P = 1.1 × 10-8 EBV nuclear antigen. We also demonstrated pleiotropy between viral response genes and complex diseases, from autoimmune disorders to cancer to neurodegenerative and psychiatric conditions. CONCLUSIONS: Our study confirms the importance of the HLA region in host response to viral infection and elucidates novel genetic determinants beyond the HLA that contribute to host-virus interaction.

A single-cell RNA-seq atlas of Schistosoma mansoni identifies a key regulator of blood feeding.

Schistosomiasis is a neglected tropical disease that infects 240 million people. With no vaccines and only one drug available, new therapeutic targets are needed. The causative agents, schistosomes, are intravascular flatworm parasites that feed on blood and lay eggs, resulting in pathology. The function of the parasite's various tissues in successful parasitism are poorly understood, hindering identification of therapeutic targets. Using single-cell RNA sequencing (RNA-seq), we characterize 43,642 cells from the adult schistosome and identify 68 distinct cell populations, including specialized stem cells that maintain the parasite's blood-digesting gut. These stem cells express the gene hnf4, which is required for gut maintenance, blood feeding, and pathology in vivo. Together, these data provide molecular insights into the organ systems of this important pathogen and identify potential therapeutic targets.

The landscape of host genetic factors involved in immune response to common viral infections.

INTRODUCTION: Humans and viruses have co-evolved for millennia resulting in a complex host genetic architecture. Understanding the genetic mechanisms of immune response to viral infection provides insight into disease etiology and therapeutic opportunities. METHODS: We conducted a comprehensive study including genome-wide and transcriptome-wide association analyses to identify genetic loci associated with immunoglobulin G antibody response to 28 antigens for 16 viruses using serological data from 7924 European ancestry participants in the UK Biobank cohort. RESULTS: Signals in human leukocyte antigen (HLA) class II region dominated the landscape of viral antibody response, with 40 independent loci and 14 independent classical alleles, 7 of which exhibited pleiotropic effects across viral families. We identified specific amino acid (AA) residues that are associated with seroreactivity, the strongest associations presented in a range of AA positions within DRß1 at positions 11, 13, 71, and 74 for Epstein-Barr Virus (EBV), Varicella Zoster Virus (VZV), Human Herpes virus 7, (HHV7) and Merkel cell polyomavirus (MCV). Genome-wide association analyses discovered 7 novel genetic loci outside the HLA associated with viral antibody response (P<5.×10-8), including FUT2 (19q13.33) for human polyomavirus BK (BKV), STING1 (5q31.2) for MCV, as well as CXCR5 (11q23.3) and TBKBP1 (17q21.32) for HHV7. Transcriptome-wide association analyses identified 114 genes associated with response to viral infection, 12 outside of the HLA region, including ECSCR: P=5.0×10-15 (MCV), NTN5: P=1.1×10-9 (BKV), and P2RY13: P=1.1×10-8 EBV nuclear antigen. We also demonstrated pleiotropy between viral response genes and complex diseases; from autoimmune disorders to cancer to neurodegenerative and psychiatric conditions. CONCLUSIONS: Our study confirms the importance of the HLA region in host response to viral infection and elucidates novel genetic determinants beyond the HLA that contribute to host-virus interaction.

MGMS2: Membrane glycolipid mass spectrum simulator for polymicrobial samples.

RATIONALE: Polymicrobial samples present unique challenges for mass spectrometric identification. A recently developed glycolipid technology has the potential to accurately identify individual bacterial species from polymicrobial samples. In order to develop and validate bacterial identification algorithms (e.g. machine learning) using this glycolipid technology, generating a large number of various polymicrobial samples can be beneficial, but it is costly and labor-intensive. Here, we propose an alternative cost-effective approach that generates realistic in silico polymicrobial glycolipid mass spectra. METHODS: We introduce MGMS2 (membrane glycolipid mass spectrum simulator) as a simulation software package that generates in silico polymicrobial membrane glycolipid matrix-assisted laser desorption/ionization time-of-flight mass spectra. Unlike currently available simulation algorithms for polymicrobial mass spectra, the proposed algorithm considers errors in m/z values and variances of intensity values, occasions of missing signature ions, and noise peaks. To our knowledge, this is the first stand-alone bacterial membrane glycolipid mass spectral simulator. MGMS2 software and its manual are freely available as an R package. An interactive MGSM2 app that helps users explore various simulation parameter options is also available. RESULTS: We demonstrated the performance of MGSM2 using six microbes. The software generated in silico glycolipid mass spectra that are similar to real polymicrobial glycolipid mass spectra. The maximum correlation between in silico mass spectra generated by MGMS2 and the real polymicrobial mass spectrum was about 87%. CONCLUSIONS: We anticipate that MGMS2, which considers spectrum-to-spectrum variation, will advance the bacterial algorithm development for polymicrobial samples.

Labeling of the Schistosoma mansoni Tegument.

Schistosomes are deadly pathogens responsible for the neglected tropical disease schistosomiasis. The parasite's virulence is aided by a skin-like tissue called the tegument. The study of the tegument is hampered by a lack of tools suitable for visualizing the tissue. Here we describe a novel methodology employing fluorophore-conjugated dextrans that allows specific fluorescent labeling of the tegument and that is compatible with downstream fluorescence-labeling techniques including phalloidin labeling, RNA FISH, and immunofluorescence.

Model-Based Spectral Library Approach for Bacterial Identification via Membrane Glycolipids.

By circumventing the need for a pure colony, MALDI-TOF mass spectrometry of bacterial membrane glycolipids (lipid A) has the potential to identify microbes more rapidly than protein-based methods. However, currently available bioinformatics algorithms (e.g., dot products) do not work well with glycolipid mass spectra such as those produced by lipid A, the membrane anchor of lipopolysaccharide. To address this issue, we propose a spectral library approach coupled with a machine learning technique to more accurately identify microbes. Here, we demonstrate the performance of the model-based spectral library approach for microbial identification using approximately a thousand mass spectra collected from multi-drug-resistant bacteria. At false discovery rates < 1%, our approach identified many more bacterial species than the existing approaches such as the Bruker Biotyper and characterized over 97% of their phenotypes accurately. As the diversity in our glycolipid mass spectral library increases, we anticipate that it will provide valuable information to more rapidly treat infected patients.

MetaMSD: meta analysis for mass spectrometry data.

Mass spectrometry-based proteomics facilitate disease understanding by providing protein abundance information about disease progression. For the same type of disease studies, multiple mass spectrometry datasets may be generated. Integrating multiple mass spectrometry datasets can provide valuable information that a single dataset analysis cannot provide. In this article, we introduce a meta-analysis software, MetaMSD (Meta Analysis for Mass Spectrometry Data) that is specifically designed for mass spectrometry data. Using Stouffer's or Pearson's test, MetaMSD detects significantly more differential proteins than the analysis based on the single best experiment. We demonstrate the performance of MetaMSD using simulated data, urinary proteomic data of kidney transplant patients, and breast cancer proteomic data. Noting the common practice of performing a pilot study prior to a main study, this software will help proteomics researchers fully utilize the benefit of multiple studies (or datasets), thus optimizing biomarker discovery. MetaMSD is a command line tool that automatically outputs various graphs and differential proteins with confidence scores. It is implemented in R and is freely available for public use at https://github.com/soyoungryu/MetaMSD. The user manual and data are available at the site. The user manual is written in such a way that scientists who are not familiar with R software can use MetaMSD.

To ERV Is Human: A Phenotype-Wide Scan Linking Polymorphic Human Endogenous Retrovirus-K Insertions to Complex Phenotypes.

Approximately 8% of the human genome is comprised of endogenous retroviral insertions (ERVs) originating from historic retroviral integration into germ cells. The function of ERVs as regulators of gene expression is well established. Less well studied are insertional polymorphisms of ERVs and their contribution to the heritability of complex phenotypes. The most recent integration of ERV, HERV-K, is expressed in a range of complex human conditions from cancer to neurologic diseases. Using an in-house computational pipeline and whole-genome sequencing data from the diverse 1,000 Genomes Phase 3 population (n = 2,504), we identified 46 polymorphic HERV-K insertions that are tagged by adjacent single nucleotide polymorphisms (SNPs). To test the potential role of polymorphic HERV-K in the heritability of complex diseases, existing databases were queried for enrichment of established relationships between the HERV-K insertion-associated SNPs (hiSNPs), and tissue specific gene expression and disease phenotypes. Overall, hiSNPs for the 46 polymorphic HERV-K sites were statistically enriched (p < 1.0E-16) for eQTLs across 44 human tissues. Fifteen of the 46 HERV-K insertions had hiSNPs annotated in the EMBL-EBI GWAS Catalog and cumulatively associated with >100 phenotypes. Experimental factor ontology enrichment analysis suggests that polymorphic HERV-K specifically contribute to neurologic and immunologic disease phenotypes, including traits related to intra cranial volume (FDR 2.00E-09), Parkinson's disease (FDR 1.80E-09), and autoimmune diseases (FDR 1.80E-09). These results provide strong candidates for context-specific study of polymorphic HERV-K insertions in disease-related traits, serving as a roadmap for future studies of the heritability of complex disease.

Flatworm-specific transcriptional regulators promote the specification of tegumental progenitors in Schistosoma mansoni.

Schistosomes infect more than 200 million people. These parasitic flatworms rely on a syncytial outer coat called the tegument to survive within the vasculature of their host. Although the tegument is pivotal for their survival, little is known about maintenance of this tissue during the decades schistosomes survive in the bloodstream. Here, we demonstrate that the tegument relies on stem cells (neoblasts) to specify fusogenic progenitors that replace tegumental cells lost to turnover. Molecular characterization of neoblasts and tegumental progenitors led to the discovery of two flatworm-specific zinc finger proteins that are essential for tegumental cell specification. These proteins are homologous to a protein essential for neoblast-driven epidermal maintenance in free-living flatworms. Therefore, we speculate that related parasites (i.e., tapeworms and flukes) employ similar strategies to control tegumental maintenance. Since parasitic flatworms infect every vertebrate species, understanding neoblast-driven tegumental maintenance could identify broad-spectrum therapeutics to fight diseases caused by these parasites.

Transplantation of schistosome sporocysts between host snails: A video guide.

Schistosomiasis is an important parasitic disease, touching roughly 200 million people worldwide. The causative agents are different Schistosoma species. Schistosomes have a complex life cycle, with a freshwater snail as intermediate host. After infection, sporocysts develop inside the snail host and give rise to human dwelling larvae. We present here a detailed step-by-step video instruction in English, French, Spanish and Portuguese that shows how these sporocysts can be manipulated and transferred from one snail to another. This procedure provides a technical basis for different types of ex vivo modifications, such as those used in functional genomics studies.