Represent Putative Virulence Factors That Are Shared by Animal Pathogenic Oomycetes, but Are Absent in Phytopathogens
Total Page:16
File Type:pdf, Size:1020Kb
Nova Southeastern University NSUWorks Biology Faculty Articles Department of Biological Sciences 10-6-2016 Glycoside Hydrolases Family 20 (GH20) Represent Putative Virulence Factors That Are Shared by Animal Pathogenic Oomycetes, but Are Absent in Phytopathogens Isabel Olivera Katrina Fins Sara Rodriguez Sumayyah K. Abiff Jaime Tartar See next page for additional authors Follow this and additional works at: https://nsuworks.nova.edu/cnso_bio_facarticles Part of the Biology Commons Authors Isabel Olivera, Katrina Fins, Sara Rodriguez, Sumayyah K. Abiff, Jaime Tartar, and Aurelien Tartar Olivera et al. BMC Microbiology (2016) 16:232 DOI 10.1186/s12866-016-0856-7 RESEARCH ARTICLE Open Access Glycoside hydrolases family 20 (GH20) represent putative virulence factors that are shared by animal pathogenic oomycetes, but are absent in phytopathogens Isabel E. Olivera1, Katrina C. Fins1, Sara A. Rodriguez1, Sumayyah K. Abiff1, Jaime L. Tartar2 and Aurélien Tartar1* Abstract Background: Although interest in animal pathogenic oomycetes is increasing, the molecular basis mediating oomycete-animal relationships remains virtually unknown. Crinkler (CRN) genes, which have been traditionally associated with the cytotoxic activity displayed by plant pathogenic oomycetes, were recently detected in transcriptome sequences from the entomopathogenic oomycete Lagenidium giganteum, suggesting that these genes may represent virulence factors conserved in both animal and plant pathogenic oomycetes. In order to further characterize the L. giganteum pathogenome, an on-going genomic survey was mined to reveal novel putative virulence factors, including canonical oomycete effectors Crinkler 13 (CRN13) orthologs. These novel sequences provided a basis to initiate gene expression analyses and determine if the proposed L. giganteum virulence factors are differentially expressed in the presence of mosquito larvae (Aedes aegypti). Results: Sequence analyses revealed that L. giganteum express CRN13 transcripts. The predicted proteins, like other L. giganteum CRNs, contained a conserved LYLA motif at the N terminal, but did not display signal peptides. In contrast, other potential virulence factors, such as Glycoside Hydrolases family 20 (hexosaminidase) and 37 (trehalase) proteins (GH20 and GH37), contained identifiable signal peptides. Genome mining demonstrated that GH20 genes are absent from phytopathogenic oomycete genomes, and that the L. giganteum GH20 sequence is the only reported peronosporalean GH20 gene. All other oomycete GH20 homologs were retrieved from animal pathogenic, saprolegnialean genomes. Furthermore, phylogenetic analyses demonstrated that saprolegnialean and peronosporalean GH20 protein sequences clustered in unrelated clades. The saprolegnialean GH20 sequences appeared as a strongly supported, monophyletic group nested within an arthropod-specific clade, suggesting that this gene was acquired via a lateral gene transfer event from an insect or crustacean genome. In contrast, the L. giganteum GH20 protein sequence appeared as a sister taxon to a plant-specific clade that included exochitinases with demonstrated insecticidal activities. Finally, gene expression analyses demonstrated that the L. giganteum GH20 gene expression level is significantly modulated in the presence of mosquito larvae. In agreement with the protein secretion predictions, CRN transcripts did not show any differential expression. Conclusions: These results identified GH20 enzymes, and not CRNs, as potential pathogenicity factors shared by all animal pathogenic oomycetes. Keywords: Crinkler, CRN, Effector, Lagenidium giganteum, Entomopathogen, Hexosaminidase * Correspondence: [email protected] 1Department of Biological Sciences, Nova Southeastern University, Fort Lauderdale, FL, USA Full list of author information is available at the end of the article © 2016 The Author(s). Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated. Olivera et al. BMC Microbiology (2016) 16:232 Page 2 of 11 Background bioactive compounds with potential as bioinsecticides Oomycetes are fungal-like heterokonts that are princi- against mosquitoes [6]. pally known as plant pathogens [1]. The most exten- An alternative hypothesis that reconciles the presence sively studied oomycete genus, Phytophthora, includes of CRN genes in the L. giganteum transcriptome and its the Irish potato famine pathogen, P. infestans, and re- pathogenicity to animals, proposes that CRN proteins mains responsible for serious economic losses to crops may play a role during mosquito infection [6]. This hy- worldwide [1]. Although over 60 % of all described pothesis is mainly supported by the fact that CRN genes oomycetes are recognized as plant pathogens, recent evi- have been detected in the genome of the chytrid fungus dence has suggested that oomycetes evolved from a mar- Batrachochytrium dendrobatidis (Bd), which is predom- ine animal pathogen [2], and that phytopathogenicity inantly known as a devastating frog pathogen that was acquired independently as a derived, apomorphic threatens natural amphibian populations worldwide [9]. trait in multiple oomycete clades [3]. The transition to Differential gene expression analyses indicated that plant pathogenicity has been associated with a dramatic BdCRN genes were up-regulated in the presence of frog expansion of the Crinkler (CRN) gene family in the skin [9], suggesting that CRN proteins may represent genomes of phylogenetically distant phytopathogenic pathogenicity factors that are also active on animal cells. oomycetes, such as the Peronosporalean P. infestans and Parallel investigations of CRN proteins originating from the Saprolegnialean Aphanomyces euteiches [4]. In ac- Bd and the plant pathogen oomycete Aphanomyces cordance with this hypothesis, CRN genes were shown euteiches have focused on one specific CRN sequence, to be absent from genome sequences generated from known as CRN13 [10], and have demonstrated that basal (non-plant pathogen) oomycetes [5]. Canonical homologous proteins of different origins shared similar CRN effector proteins are characterized by the con- functions at the molecular level (DNA binding), concen- served LxLYLA or LxLFLA motifs at the N terminal, trate to similar location in the respective host cells which have been implicated in the transport of these (nucleus), and lead to similar outcomes (host cell death). proteins in the host cells during plant infection. Follow- The L. giganteum CRN proteins have yet to be included ing translocation in the host cells, CRN proteins accu- in such comparative analyses, and it remains unclear if mulate in the nucleus, where they induce cell death [4]. they include a CRN13 homolog, and if they have any The entomopathogen Lagenidium giganteum repre- role in the pathogenicity process. Overall, the molecular sents a unique, extant animal pathogenic oomycete that processes mediating mosquito infections by L. giganteum has been shown to express canonical CRN oomycete ef- remain uncharacterized. A recent study reported that L. fector genes [6]. This observation, combined with phylo- giganteum secretes GH5_27 enzymes that appear to be genetic analyses that placed L. giganteum nested within mostly specific to cuticle degrading entomopathogens, a Peronosporalean clade of plant pathogens, suggested including not only entomopathogenic oomycetes but that L. giganteum evolved from a plant pathogenic also Fungi [6]. Earlier reports have proposed trehalases ancestor, and may have reverted back towards a as potential pathogenicity factors, based on the observa- plesiomorphic-like state. The presence of CRN sequences tions that trehalose is the most abundant sugar source in the L. giganteum transcriptome contrasted with historic in insects’ hemolymph, and that trehalose depletion may observations describing this oomycete as a mosquito contribute to the L. giganteum infection process and pathogen with a narrow range of invertebrate hosts, and ultimate death of the insect host [11]. little impact on plant tissues [6]. Interestingly, the true na- In this study, an on-going survey of the L. giganteum ture of the L. giganteum host range has also been recently genome and transcriptome was used to characterize the challenged by reports of L. giganteum infections in mam- L. giganteum gene sequences for CRN13 and trehalase mals [7]. These reports have diminished the potential of L. (Glycoside Hydrolase family 37, or GH37) orthologs. giganteum as a biocontrol agent against mosquitoes. How- The CRN13 and GH37 nucleotide sequences provided a ever, they also have contributed to reinforce the original basis to initiate differential gene expression studies, and assertions that L. giganteum is defined as a mosquito estimate if these genes are up-regulated in the presence pathogen, since this characteristic phenotypical feature of Aedes aegypti mosquito larvae, which are demon- was used to complement molecular-based phylogenetic strated natural hosts for L. giganteum [12]. The differen- analyses, and validate