Comparative genomics reveals selective distribution and domain organization of FYVE and PX domain proteins across eukaryotic lineages

PsAF5 functions as an essential adapter for PsPHB2-mediated mitophagy under ROS stress in Phytophthora sojae

Host-derived reactive oxygen species (ROS) are an important defense means to protect against pathogens. Although mitochondria are the main intracellular targets of ROS, how pathogens regulate mitochondrial physiology in response to oxidative stress remains elusive. Prohibitin 2 (PHB2) is an inner mitochondrial membrane (IMM) protein, recognized as a mitophagy receptor in animals and fungi. Here, we find that an ANK and FYVE domain-containing protein PsAF5, is an adapter of PsPHB2, interacting with PsATG8 under ROS stress. Unlike animal PHB2 that can recruit ATG8 directly to mitochondria, PsPHB2 in Phytophthora sojae cannot recruit PsATG8 to stressed mitochondria without PsAF5. PsAF5 deletion impairs mitophagy under ROS stress and increases the pathogen's sensitivity to H2O2, resulting in the attenuation of P. sojae virulence. This discovery of a PsPHB2-PsATG8 adapter (PsAF5) in plant-pathogenic oomycetes reveals that mitophagy induction by IMM proteins is conserved in eukaryotes, but with differences in the details of ATG8 recruitment.

A Novel FYVE Domain-Containing Protein Kinase, PsZFPK1, Plays a Critical Role in Vegetative Growth, Sporangium Formation, Oospore Production, and Virulence in Phytophthora sojae

Proteins containing both FYVE and serine/threonine kinase catalytic (STKc) domains are exclusive to protists. However, the biological function of these proteins in oomycetes has rarely been reported. In the Phytophthora sojae genome database, we identified five proteins containing FYVE and STKc domains, which we named PsZFPK1, PsZFPK2, PsZFPK3, PsZFPK4, and PsZFPK5. In this study, we characterized the biological function of PsZFPK1 using a CRISPR/Cas9-mediated gene replacement system. Compared with the wild-type strain, P6497, the PsZFPK1-knockout mutants exhibited significantly reduced growth on a nutrient-rich V8 medium, while a more pronounced defect was observed on a nutrient-poor Plich medium. The PsZFPK1-knockout mutants also showed a significant increase in sporangium production. Furthermore, PsZFPK1 was found to be essential for oospore production and complete virulence but dispensable for the stress response in P. sojae. The N-terminal region, FYVE and STKc domains, and T602 phosphorylation site were found to be vital for the function of PsZFPK1. Conversely, these domains were not required for the localization of PsZFPK1 protein in the cytoplasm. Our results demonstrate that PsZFPK1 plays a critical role in vegetative growth, sporangium formation, oospore production, and virulence in P. sojae.

The phosphoinositide code is read by a plethora of protein domains

INTRODUCTION

The proteins that decipher nucleic acid- and protein-based information are well known, however, those that read membrane-encoded information remain understudied. Here we report 70 different human, microbial and viral protein folds that recognize phosphoinositides (PIs), comprising the readers of a vast membrane code.

AREAS COVERED

Membrane recognition is best understood for FYVE, PH and PX domains, which exemplify hundreds of PI code readers. Comparable lipid interaction mechanisms may be mediated by kinases, adjacent C1 and C2 domains, trafficking arrestin, GAT and VHS modules, membrane-perturbing annexin, BAR, CHMP, ENTH, HEAT, syntaxin and Tubby helical bundles, multipurpose FERM, EH, MATH, PHD, PDZ, PROPPIN, PTB and SH2 domains, as well as systems that regulate receptors, GTPases and actin filaments, transfer lipids and assembled bacterial and viral particles.

EXPERT OPINION

The elucidation of how membranes are recognized has extended the genetic code to the PI code. Novel discoveries include PIP-stop and MET-stop residues to which phosphates and metabolites are attached to block phosphatidylinositol phosphate (PIP) recognition, memteins as functional membrane protein apparatuses, and lipidons as lipid "codons" recognized by membrane readers. At least 5% of the human proteome senses such membrane signals and allows eukaryotic organelles and pathogens to operate and replicate.

Phosphoinositide Recognition Sites Are Blocked by Metabolite Attachment

Membrane readers take part in trafficking and signaling processes by localizing proteins to organelle surfaces and transducing molecular information. They accomplish this by engaging phosphoinositides (PIs), a class of lipid molecules which are found in different proportions in various cellular membranes. The prototypes are the PX domains, which exhibit a range of specificities for PIs. Our meta-analysis indicates that recognition of membranes by PX domains is specifically controlled by modification of lysine and arginine residues including acetylation, hydroxyisobutyrylation, glycation, malonylation, methylation and succinylation of sidechains that normally bind headgroups of phospholipids including organelle-specific PI signals. Such metabolite-modulated residues in lipid binding elements are named MET-stops here to highlight their roles as erasers of membrane reader functions. These modifications are concentrated in the membrane binding sites of half of all 49 PX domains in the human proteome and correlate with phosphoregulatory sites, as mapped using the Membrane Optimal Docking Area (MODA) algorithm. As these motifs are mutated and modified in various cancers and the responsible enzymes serve as potential drug targets, the discovery of MET-stops as a widespread inhibitory mechanism may aid in the development of diagnostics and therapeutics aimed at the readers, writers and erasers of the PI code.

An FYVE-Domain-Containing Protein, PsFP1, Is Involved in Vegetative Growth, Oxidative Stress Response and Virulence of Phytophthora sojae

Proteins that contain the FYVE zinc-finger domain are recruited to PtdIns3P-containing membranes, participating in numerous biological processes such as membrane trafficking, cytoskeletal regulation, and receptor signaling. However, the genome-wide distribution, evolution, and biological functions of FYVE-containing proteins are rarely reported for oomycetes. By genome mining of Phytophthora sojae, two proteins (PsFP1 and PsFP2) with a combination of the FYVE domain and the PX domain (a major phosphoinositide binding module) were found. To clarify the functions of PsFP1 and PsFP2, the CRISPR/Cas9-mediated gene replacement system was used to knock out the two genes respectively. Only heterozygous deletion mutants of PsFP1 were recovered, and the expression level of PsFP1 in the heterozygous knockout transformants was significantly down-regulated. These PsFP1 mutants showed a decrease in mycelial growth and pathogenicity and were more sensitive to hydrogen peroxide. These phenotypes were recovered to the level of wild-type by overexpression PsFP1 gene in the PsFP1 heterozygous knockout transformant. In contrast, deletion of PsFP2 had no significant effect on vegetative growth, asexual and sexual reproduction, pathogenicity, or oxidative stress sensitivity. PsFP1 was primarily localized in vesicle-like structures and both the FYVE and PX domains are important for its localization. Overall, our results indicate that PsFP1 plays an important role in the vegetative growth and virulence of P. sojae.

Diversity of phosphoinositide binding proteins in Entamoeba histolytica

Phosphatidylinositol phosphates (PIPs, phosphoinositides) are localized to the membranes of all cellular compartments, and play pivotal roles in multiple cellular events. To fulfill their functions, PIPs that are located to specific organelles or membrane domains bind to and recruit various proteins in spatiotemporal specific manner via protein domains that selectively bind to either a single or an array of PIPs. In Entamoeba histolytica, the human intestinal protozoan parasite, PIPs and PIP-binding proteins have been shown to be involved in their virulence-associated mechanisms such as cell motility, vesicular traffic, trogo- and phagocytosis. In silico search of the domains and the signatures implicated in PIP binding in the E. histolytica proteome allows identification of dozens of potential PIP-binding proteins. However, such analysis is often misleading unless the protein domain used as query is cautiously selected and the binding specificity of the proteins are experimentally validated. This is because all the domains initially presumed to bind PIPs in other systems are not always capable of PIP binding, but rather involved in other biological roles. In this review, we carried out in silico survey of proteins which have PIP-binding domains in the E. histolytica genome by utilizing only validated PIP-binding domains that had been experimentally proven to be faithful PIP-binding bioprobes. Our survey has identified that FYVE (Fab1, YOTB1, Vac1, EEA1) and PH (pleckstrin homology) domain containing proteins are the most expanded families in E. histolytica. A few FYVE domain-containing proteins (EhFP4 and 10) and phox homology (PX) domain containing proteins (EhSNX1 and 2) were previously studied in depth in E. histolytica. Furthermore, most of the identified PH domain-containing proteins are annotated as protein kinases and possess protein kinase domains. Overall, PIP-binding domain-containing proteins that can be identified by in silico survey of the genome using the domains from well characterized bioprobes are limited in E. histolytica. However, their domain architectures are often unique, suggesting unique evolution of PIP-binding domain-containing proteins in this organism.

RMS2 Encoding a GDSL Lipase Mediates Lipid Homeostasis in Anthers to Determine Rice Male Fertility

Plant male gametogenesis is a coordinated effort involving both reproductive tissues and sporophytic tissues, in which lipid metabolism plays an essential role. Although GDSL esterases/lipases have been well known as key enzymes for many plant developmental processes and stress responses, their functions in reproductive development remain unclear. Here, we report the identification of a rice male sterile2 (rms2) mutant in rice (Oryza sativa), which is completely male sterile due to the defects in tapetum degradation, cuticle formation in sporophytic tissues, and impaired exine and central vacuole development in pollen grains. RMS2 was map-based cloned as an endoplasmic reticulum-localized GDSL lipase gene, which is predominantly transcribed during early anther development. In rms2, a three-nucleotide deletion and one base substitution (TTGT to A) occurred within the GDSL domain, which reduced the lipid hydrolase activity of the resulting protein and led to significant changes in the content of 16 lipid components and numerous other metabolites, as revealed by a comparative metabolic analysis. Furthermore, RMS2 is directly targeted by the male fertility regulators Undeveloped Tapetum1 and Persistent Tapetal Cell1 both in vitro and in vivo, suggesting that RMS2 may serve as a key node in the rice male fertility regulatory network. These findings shed light on the function of GDSLs in reproductive development and provide a promising gene resource for hybrid rice breeding.

Evolutionary insights into FYVE and PHOX effector proteins from the moss Physcomitrella patens

Genome-wide identification, together with gene expression patterns and promoter region analysis of FYVE and PHOX proteins in Physcomitrella patens, emphasized their importance in regulating mainly developmental processes in P. patens. Phosphatidylinositol 3-phosphate (PtdIns3P) is a signaling phospholipid, which regulates several aspects of plant growth and development, as well as responses to biotic and abiotic stresses. The mechanistic insights underlying PtdIns3P mode of action, specifically through effector proteins have been partially explored in plants, with main focus on Arabidopsis thaliana. In this study, we searched for genes coding for PtdIns3P-binding proteins such as FYVE and PHOX domain-containing sequences from different photosynthetic organisms to gather evolutionary insights on these phosphoinositide binding domains, followed by an in silico characterization of the FYVE and PHOX gene families in the moss Physcomitrella patens. Phylogenetic analysis showed that PpFYVE proteins can be grouped in 7 subclasses, with an additional subclass whose FYVE domain was lost during evolution to higher plants. On the other hand, PpPHOX proteins are classified into 5 subclasses. Expression analyses based on RNAseq data together with the analysis of cis-acting regulatory elements and transcription factor (TF) binding sites in promoter regions suggest the importance of these proteins in regulating stress responses but mainly developmental processes in P. patens. The results provide valuable information and robust candidate genes for future functional analysis aiming to further explore the role of this signaling pathway mainly during growth and development of tip growing cells and during the transition from 2 to 3D growth. These studies would identify ancestral regulatory players undertaken during plant evolution.

Analysis of the protein domain and domain architecture content in fungi and its application in the search of new antifungal targets

Over the past several years fungal infections have shown an increasing incidence in the susceptible population, and caused high mortality rates. In parallel, multi-resistant fungi are emerging in human infections. Therefore, the identification of new potential antifungal targets is a priority. The first task of this study was to analyse the protein domain and domain architecture content of the 137 fungal proteomes (corresponding to 111 species) available in UniProtKB (UniProt KnowledgeBase) by January 2013. The resulting list of core and exclusive domain and domain architectures is provided in this paper. It delineates the different levels of fungal taxonomic classification: phylum, subphylum, order, genus and species. The analysis highlighted Aspergillus as the most diverse genus in terms of exclusive domain content. In addition, we also investigated which domains could be considered promiscuous in the different organisms. As an application of this analysis, we explored three different ways to detect potential targets for antifungal drugs. First, we compared the domain and domain architecture content of the human and fungal proteomes, and identified those domains and domain architectures only present in fungi. Secondly, we looked for information regarding fungal pathways in public repositories, where proteins containing promiscuous domains could be involved. Three pathways were identified as a result: lovastatin biosynthesis, xylan degradation and biosynthesis of siroheme. Finally, we classified a subset of the studied fungi in five groups depending on their occurrence in clinical samples. We then looked for exclusive domains in the groups that were more relevant clinically and determined which of them had the potential to bind small molecules. Overall, this study provides a comprehensive analysis of the available fungal proteomes and shows three approaches that can be used as a first step in the detection of new antifungal targets.

Some fungi have become pathogenic to plants and in a lesser extent to animals. Under certain conditions their presence in the human body can prove a threat for human health, especially for immunocompromised patients. Yet, some fungi can also infect healthy individuals. The low sensitivity of the antifungal drugs available together with the clinically observed resistance of some fungi raises the demand for new alternative treatments. Proteins are biological molecules which perform essential functions within the living organisms. Many of those functions are attributed to the varying folded structure of each protein. These configurations are composed of functional units -also called domains- each one independently responsible for a fraction of the overall biological function. Understanding how the different block combinations are distributed across members of the same or similar families of organisms is important. For instance, exclusive domain combinations can hold particular acquired functions. Blocks displaying a high mobility can play major roles for the organism's survival. The biological goal of this study was to analyse the functional implications of protein domains and domain combinations in the available fungal proteomes. This information can be used to highlight proteins and pathways that could be potentially used as drug targets.

Involvement of Phosphatidylinositol 3-kinase in the regulation of proline catabolism in Arabidopsis thaliana

Plant adaptation to abiotic stresses such as drought and salinity involves complex regulatory processes. Deciphering the signaling components that are involved in stress signal transduction and cellular responses is of importance to understand how plants cope with salt stress. Accumulation of osmolytes such as proline is considered to participate in the osmotic adjustment of plant cells to salinity. Proline accumulation results from a tight regulation between its biosynthesis and catabolism. Lipid signal components such as phospholipases C and D have previously been shown to be involved in the regulation of proline metabolism in Arabidopsis thaliana. In this study, we demonstrate that proline metabolism is also regulated by class-III Phosphatidylinositol 3-kinase (PI3K), VPS34, which catalyses the formation of phosphatidylinositol 3-phosphate (PI3P) from phosphatidylinositol. Using pharmacological and biochemical approaches, we show that the PI3K inhibitor, LY294002, affects PI3P levels in vivo and that it triggers a decrease in proline accumulation in response to salt treatment of A. thaliana seedlings. The lower proline accumulation is correlated with a lower transcript level of Pyrroline-5-carboxylate synthetase 1 (P5CS1) biosynthetic enzyme and higher transcript and protein levels of Proline dehydrogenase 1 (ProDH1), a key-enzyme in proline catabolism. We also found that the ProDH1 expression is induced in a pi3k-hemizygous mutant, further demonstrating that PI3K is involved in the regulation of proline catabolism through transcriptional regulation of ProDH1. A broader metabolomic analysis indicates that LY294002 also reduced other metabolites, such as hydrophobic and aromatic amino acids and sugars like raffinose.

Evolutionary reconstruction of the retromer complex and its function in Trypanosoma brucei

Intracellular trafficking and protein sorting are mediated by various protein complexes, with the retromer complex being primarily involved in retrograde traffic from the endosome or lysosome to the Golgi complex. Here, comparative genomics, cell biology and phylogenetics were used to probe the early evolution of retromer and its function. Retromer subunits Vps26, Vps29 and Vps35 are near universal, and, by inference, the complex was an ancient feature of eukaryotic cells. Surprisingly, we found DSCR3, a Vps26 paralogue in humans associated with Down's syndrome, in at least four eukaryotic supergroups, implying a more ancient origin than previously suspected. By contrast, retromer cargo proteins showed considerable interlineage variability, with lineage-specific and broadly conserved examples found. Vps10 trafficking probably represents an ancestral role for the complex. Vps5, the BAR-domain-containing membrane-deformation subunit, was found in diverse eukaryotes, including in the divergent eukaryote Trypanosoma brucei, where it is the first example of a BAR-domain protein. To determine functional conservation, an initial characterisation of retromer was performed in T. brucei; the endosomal localisation and its role in endosomal targeting are conserved. Therefore retromer is identified as a further feature of the sophisticated intracellular trafficking machinery of the last eukaryotic common ancestor, with BAR domains representing a possible third independent mechanism of membrane-deformation arising in early eukaryotes.

Phylogenomics of phosphoinositide lipid kinases: perspectives on the evolution of second messenger signaling and drug discovery

Background

Phosphoinositide lipid kinases (PIKs) generate specific phosphorylated variants of phosatidylinositols (PtdIns) that are critical for second messenger signaling and cellular membrane remodeling. Mammals have 19 PIK isoforms spread across three major families: the PtIns 3-kinases (PI3Ks), PtdIns 4-kinases (PI4Ks), and PtdIns-P (PIP) kinases (PIPKs). Other eukaryotes have fewer yet varying PIK complements. PIKs are also an important, emerging class of drug targets for many therapeutic areas including cancer, inflammatory and metabolic diseases and host-pathogen interactions. Here, we report the genomic occurrences and evolutionary relationships or phylogenomics of all three PIK families across major eukaryotic groups and suggest potential ramifications for drug discovery.

Results

Our analyses reveal four core eukaryotic PIKs which are type III PIK4A and PIK4B, and at least one homolog each from PI3K (possibly PIK3C3 as the ancestor) and PIP5K families. We also applied evolutionary analyses to PIK disease ontology and drug discovery. Mutated PIK3CA are known to be oncogenic and several inhibitors are in anti-cancer clinical trials. We found conservation of activating mutations of PIK3CA in paralogous isoforms suggesting specific functional constraints on these residues. By mapping published compound inhibition data (IC50s) onto a phylogeny of PI3Ks, type II PI4Ks and distantly related, MTOR, ATM, ATR and PRKDC kinases, we also show that compound polypharmacology corresponds to kinase evolutionary relationships. Finally, we extended the rationale for drugs targeting PIKs of malarial Plasmodium falciparum, and the parasites, Leishmania sp. and Trypanosoma sp. by identifying those PIKs highly divergent from human homologs.

Conclusion

Our phylogenomic analysis of PIKs provides new insights into the evolution of second messenger signaling. We postulate two waves of PIK diversification, the first in metazoans with a subsequent expansion in cold-blooded vertebrates that was post-emergence of Deutrostomia\Chordata but prior to the appearance of mammals. Reconstruction of the evolutionary relationships among these lipid kinases also adds to our understanding of their roles in various diseases and assists in their development as potential drug targets.

Defining the role of a FYVE domain in the localization and activity of a cAMP phosphodiesterase implicated in osmoregulation in Trypanosoma cruzi

Intracellular levels of cyclic nucleotide second messengers are regulated predominantly by a large superfamily of phosphodiesterases (PDEs). Trypanosoma cruzi, the causative agent of Chagas disease, encodes four different PDE families. One of these PDEs, T. cruzi PDE C2 (TcrPDEC2) has been characterized as a FYVE domain containing protein. Here, we report a novel role for TcrPDEC2 in osmoregulation in T. cruzi and reveal the relevance of its FYVE domain. Our data show that treatment of epimastigotes with TcrPDEC2 inhibitors improves their regulatory volume decrease, whereas cells overexpressing this enzyme are unaffected by the same inhibitors. Consistent with these results, TcrPDEC2 localizes to the contractile vacuole complex, showing strong labelling in the region corresponding to the spongiome. Furthermore, transgenic parasites overexpressing a truncated version of TcrPDEC2 without the FYVE domain show a failure in its targeting to the contractile vacuole complex and a marked decrease in PDE activity, supporting the importance of this domain to the localization and activity of TcrPDEC2. Taking together, the results here presented are consistent with the importance of the cyclic AMP signalling pathway in regulatory volume decrease and implicate TcrPDEC2 as a specifically localized PDE involved in osmoregulation in T. cruzi.