Profiling the human protein-DNA interactome reveals ERK2 as a transcriptional repressor of interferon signaling.

Protein-DNA interactions (PDIs) mediate a broad range of functions essential for cellular differentiation, function, and survival. However, it is still a daunting task to comprehensively identify and profile sequence-specific PDIs in complex genomes. Here, we have used a combined bioinformatics and protein microarray-based strategy to systematically characterize the human protein-DNA interactome. We identified 17,718 PDIs between 460 DNA motifs predicted to regulate transcription and 4,191 human proteins of various functional classes. Among them, we recovered many known PDIs for transcription factors (TFs). We identified a large number of unanticipated PDIs for known TFs, as well as for previously uncharacterized TFs. We also found that over three hundred unconventional DNA-binding proteins (uDBPs)--which include RNA-binding proteins, mitochondrial proteins, and protein kinases--showed sequence-specific PDIs. One such uDBP, ERK2, acts as a transcriptional repressor for interferon gamma-induced genes, suggesting important biological roles for such proteins.

Multiplexed Biomarker Panels Discriminate Zika and Dengue Virus Infection in Humans.

Zika virus (ZIKV) and dengue virus (DENV) are closely related flaviviruses that cause widespread, acute febrile illnesses, notably microcephaly for fetuses of infected pregnant women. Detecting the viral cause of these illnesses is paramount to determine risks to patients, counsel pregnant women, and help fight outbreaks. A combined diagnostic algorithm for ZIKV and DENV requires Reverse transcription polymerase chain reaction (RT-PCR) and IgM antibody detection. RT-PCR differentiates between DENV and ZIKV infections during the acute phases of infection, but differentiation based on IgM antibodies is currently nearly impossible in endemic areas. We have developed a ZIKV/DENV protein array and tested it with serum samples collected from ZIKV- and DENV-infected patients and healthy subjects in Puerto Rico. Our analyses reveal a biomarker panel that are capable of discriminating ZIKV and DENV infections with high accuracy, including Capsid protein from African ZIKV strain MR766, and other 5 pair of proteins (NS1, NS2A, NS3, NS4B and NS5) from ZIKV and DENV respectively. Both sensitivity and specificity of the test for ZIKV from DENV are around 90%. We propose that the ZIKV/DENV protein array will be used in future studies to discriminate patients infected with ZIKV from DENV.

Global Identification of Small Ubiquitin-related Modifier (SUMO) Substrates Reveals Crosstalk between SUMOylation and Phosphorylation Promotes Cell Migration.

Proteomics studies have revealed that SUMOylation is a widely used post-translational modification (PTM) in eukaryotes. However, how SUMO E1/2/3 complexes use different SUMO isoforms and recognize substrates remains largely unknown. Using a human proteome microarray-based activity screen, we identified over 2500 proteins that undergo SUMO E3-dependent SUMOylation. We next constructed a SUMO isoform- and E3 ligase-dependent enzyme-substrate relationship network. Protein kinases were significantly enriched among SUMOylation substrates, suggesting crosstalk between phosphorylation and SUMOylation. Cell-based analyses of tyrosine kinase, PYK2, revealed that SUMOylation at four lysine residues promoted PYK2 autophosphorylation at tyrosine 402, which in turn enhanced its interaction with SRC and full activation of the SRC-PYK2 complex. SUMOylation on WT but not the 4KR mutant of PYK2 further elevated phosphorylation of the downstream components in the focal adhesion pathway, such as paxillin and Erk1/2, leading to significantly enhanced cell migration during wound healing. These studies illustrate how our SUMO E3 ligase-substrate network can be used to explore crosstalk between SUMOylation and other PTMs in many biological processes.

A Human Proteome Array Approach to Identifying Key Host Proteins Targeted by Toxoplasma Kinase ROP18.

Toxoplasma kinase ROP18 is a key molecule responsible for the virulence of Toxoplasma gondii; however, the mechanisms by which ROP18 exerts parasite virulence via interaction with host proteins remain limited to a small number of identified substrates. To identify a broader array of ROP18 substrates, we successfully purified bioactive mature ROP18 and used it to probe a human proteome array. Sixty eight new putative host targets were identified. Functional annotation analysis suggested that these proteins have a variety of functions, including metabolic process, kinase activity and phosphorylation, cell growth, apoptosis and cell death, and immunity, indicating a pleiotropic role of ROP18 kinase. Among these proteins, four candidates, p53, p38, UBE2N, and Smad1, were further validated. We demonstrated that ROP18 targets p53, p38, UBE2N, and Smad1 for degradation. Importantly, we demonstrated that ROP18 phosphorylates Smad1 Ser-187 to trigger its proteasome-dependent degradation. Further functional characterization of the substrates of ROP18 may enhance understanding of the pathogenesis of Toxoplasma infection and provide new therapeutic targets. Similar strategies could be used to identify novel host targets for other microbial kinases functioning at the pathogen-host interface.

Genome-wide functional analysis of pathogenicity genes in the rice blast fungus.

Rapid translation of genome sequences into meaningful biological information hinges on the integration of multiple experimental and informatics methods into a cohesive platform. Despite the explosion in the number of genome sequences available, such a platform does not exist for filamentous fungi. Here we present the development and application of a functional genomics and informatics platform for a model plant pathogenic fungus, Magnaporthe oryzae. In total, we produced 21,070 mutants through large-scale insertional mutagenesis using Agrobacterium tumefaciens-mediated transformation. We used a high-throughput phenotype screening pipeline to detect disruption of seven phenotypes encompassing the fungal life cycle and identified the mutated gene and the nature of mutation for each mutant. Comparative analysis of phenotypes and genotypes of the mutants uncovered 202 new pathogenicity loci. Our findings demonstrate the effectiveness of our platform and provide new insights on the molecular basis of fungal pathogenesis. Our approach promises comprehensive functional genomics in filamentous fungi and beyond.

Global analysis of phosphorylation networks in humans.

Phosphorylation-mediated signaling plays a crucial role in nearly every aspect of cellular physiology. A recent study based on protein microarray experiments identified a large number of kinase-substrate relationships (KSRs), and built a comprehensive and reliable phosphorylation network in humans. Analysis of this network, in conjunction with additional resources, revealed several key features. First, comparison of the human and yeast phosphorylation networks uncovered an evolutionarily conserved signaling backbone dominated by kinase-to-kinase relationships. Second, although most of the KSRs themselves are not conserved, the functions enriched in the substrates for a given kinase are often conserved. Third, the prevalence of kinase-transcription factor regulatory modules suggests that phosphorylation and transcriptional regulatory networks are inherently wired together to form integrated regulatory circuits. Overall, the phosphorylation networks described in this work promise to offer new insights into the properties of kinase signaling pathways, at both the global and the protein levels. This article is part of a Special Issue entitled: Computational Proteomics, Systems Biology & Clinical Implications. Guest Editor: Yudong Cai.

PhosphoNetworks: a database for human phosphorylation networks.

SUMMARY: Phosphorylation plays an important role in cellular signal transduction. Current phosphorylation-related databases often focus on the phosphorylation sites, which are mainly determined by mass spectrometry. Here, we present PhosphoNetworks, a phosphorylation database built on a high-resolution map of phosphorylation networks. This high-resolution map of phosphorylation networks provides not only the kinase-substrate relationships (KSRs), but also the specific phosphorylation sites on which the kinases act on the substrates. The database contains the most comprehensive dataset for KSRs, including the relationships from a recent high-throughput project for identification of KSRs using protein microarrays, as well as known KSRs curated from the literature. In addition, the database also includes several analytical tools for dissecting phosphorylation networks. PhosphoNetworks is expected to play a prominent role in proteomics and phosphorylation-related disease research. AVAILABILITY AND IMPLEMENTATION: http://www.phosphonetworks.org

Global expression profiling of transcription factor genes provides new insights into pathogenicity and stress responses in the rice blast fungus.

Because most efforts to understand the molecular mechanisms underpinning fungal pathogenicity have focused on studying the function and role of individual genes, relatively little is known about how transcriptional machineries globally regulate and coordinate the expression of a large group of genes involved in pathogenesis. Using quantitative real-time PCR, we analyzed the expression patterns of 206 transcription factor (TF) genes in the rice blast fungus Magnaporthe oryzae under 32 conditions, including multiple infection-related developmental stages and various abiotic stresses. The resulting data, which are publicly available via an online platform, provided new insights into how these TFs are regulated and potentially work together to control cellular responses to a diverse array of stimuli. High degrees of differential TF expression were observed under the conditions tested. More than 50% of the 206 TF genes were up-regulated during conidiation and/or in conidia. Mutations in ten conidiation-specific TF genes caused defects in conidiation. Expression patterns in planta were similar to those under oxidative stress conditions. Mutants of in planta inducible genes not only exhibited sensitive to oxidative stress but also failed to infect rice. These experimental validations clearly demonstrated the value of TF expression patterns in predicting the function of individual TF genes. The regulatory network of TF genes revealed by this study provides a solid foundation for elucidating how M. oryzae regulates its pathogenesis, development, and stress responses.

Phosphorylation of the chromatin binding domain of KSHV LANA.

The Kaposi sarcoma associated herpesvirus (KSHV) latency associated nuclear antigen (LANA) is expressed in all KSHV associated malignancies and is essential for maintenance of KSHV genomes in infected cells. To identify kinases that are potentially capable of modifying LANA, in vitro phosphorylation assays were performed using an Epstein Barr virus plus LANA protein microarray and 268 human kinases purified in active form from yeast. Interestingly, of the Epstein-Barr virus proteins on the array, the EBNA1 protein had the most similar kinase profile to LANA. We focused on nuclear kinases and on the N-terminus of LANA (amino acids 1-329) that contains the LANA chromatin binding domain. Sixty-three nuclear kinases phosphorylated the LANA N-terminus. Twenty-four nuclear kinases phosphorylated a peptide covering the LANA chromatin binding domain (amino acids 3-21). Alanine mutations of serine 10 and threonine 14 abolish or severely diminish chromatin and histone binding by LANA. However, conversion of these residues to the phosphomimetic glutamic acid restored histone binding suggesting that phosphorylation of serine 10 and threonine 14 may modulate LANA function. Serine 10 and threonine 14 were validated as substrates of casein kinase 1, PIM1, GSK-3 and RSK3 kinases. Short-term treatment of transfected cells with inhibitors of these kinases found that only RSK inhibition reduced LANA interaction with endogenous histone H2B. Extended treatment of PEL cell cultures with RSK inhibitor caused a decrease in LANA protein levels associated with p21 induction and a loss of PEL cell viability. The data indicate that RSK phosphorylation affects both LANA accumulation and function.

Construction of human activity-based phosphorylation networks.

The landscape of human phosphorylation networks has not been systematically explored, representing vast, unchartered territories within cellular signaling networks. Although a large number of in vivo phosphorylated residues have been identified by mass spectrometry (MS)-based approaches, assigning the upstream kinases to these residues requires biochemical analysis of kinase-substrate relationships (KSRs). Here, we developed a new strategy, called CEASAR, based on functional protein microarrays and bioinformatics to experimentally identify substrates for 289 unique kinases, resulting in 3656 high-quality KSRs. We then generated consensus phosphorylation motifs for each of the kinases and integrated this information, along with information about in vivo phosphorylation sites determined by MS, to construct a high-resolution map of phosphorylation networks that connects 230 kinases to 2591 in vivo phosphorylation sites in 652 substrates. The value of this data set is demonstrated through the discovery of a new role for PKA downstream of Btk (Bruton's tyrosine kinase) during B-cell receptor signaling. Overall, these studies provide global insights into kinase-mediated signaling pathways and promise to advance our understanding of cellular signaling processes in humans.

Regulation of CK2 by phosphorylation and O-GlcNAcylation revealed by semisynthesis.

Protein serine-threonine kinase casein kinase II (CK2) is involved in a myriad of cellular processes including cell growth and proliferation through its phosphorylation of hundreds of substrates, yet how CK2 function is regulated is poorly understood. Here we report that the CK2 catalytic subunit CK2α is modified by O-linked ß-N-acetyl-glucosamine (O-GlcNAc) on Ser347, proximal to a cyclin-dependent kinase phosphorylation site (Thr344). We use protein semisynthesis to show that phosphorylation of Thr344 increases the cellular stability of CK2α by strengthening its interaction with Pin1, whereas glycosylation of Ser347 seems to be antagonistic to Thr344 phosphorylation and permissive to proteasomal degradation. By performing kinase assays with site-specifically phospho- and glyco-modified CK2α in combination with CK2ß and Pin1 binding partners on human protein microarrays, we show that the kinase substrate selectivity of CK2 is modulated by these specific post-translational modifications. This study suggests how a promiscuous protein kinase can be regulated at multiple levels to achieve particular biological outputs.

Rapid identification of monospecific monoclonal antibodies using a human proteome microarray.

To broaden the range of tools available for proteomic research, we generated a library of 16,368 unique full-length human ORFs that are expressible as N-terminal GST-His(6) fusion proteins. Following expression in yeast, these proteins were then individually purified and used to construct a human proteome microarray. To demonstrate the usefulness of this reagent, we developed a streamlined strategy for the production of monospecific monoclonal antibodies that used immunization with live human cells and microarray-based analysis of antibody specificity as its central components. We showed that microarray-based analysis of antibody specificity can be performed efficiently using a two-dimensional pooling strategy. We also demonstrated that our immunization and selection strategies result in a large fraction of monospecific monoclonal antibodies that are both immunoblot and immunoprecipitation grade. Our data indicate that the pipeline provides a robust platform for the generation of monoclonal antibodies of exceptional specificity.

Homeobox transcription factors are required for conidiation and appressorium development in the rice blast fungus Magnaporthe oryzae.

The appropriate development of conidia and appressoria is critical in the disease cycle of many fungal pathogens, including Magnaporthe oryzae. A total of eight genes (MoHOX1 to MoHOX8) encoding putative homeobox transcription factors (TFs) were identified from the M. oryzae genome. Knockout mutants for each MoHOX gene were obtained via homology-dependent gene replacement. Two mutants, DeltaMohox3 and DeltaMohox5, exhibited no difference to wild-type in growth, conidiation, conidium size, conidial germination, appressorium formation, and pathogenicity. However, the DeltaMohox1 showed a dramatic reduction in hyphal growth and increase in melanin pigmentation, compared to those in wild-type. DeltaMohox4 and DeltaMohox6 showed significantly reduced conidium size and hyphal growth, respectively. DeltaMohox8 formed normal appressoria, but failed in pathogenicity, probably due to defects in the development of penetration peg and invasive growth. It is most notable that asexual reproduction was completely abolished in DeltaMohox2, in which no conidia formed. DeltaMohox2 was still pathogenic through hypha-driven appressoria in a manner similar to that of the wild-type. However, DeltaMohox7 was unable to form appressoria either on conidial germ tubes, or at hyphal tips, being non-pathogenic. These factors indicate that M. oryzae is able to cause foliar disease via hyphal appressorium-mediated penetration, and MoHOX7 is mutually required to drive appressorium formation from hyphae and germ tubes. Transcriptional analyses suggest that the functioning of M. oryzae homeobox TFs is mediated through the regulation of gene expression and is affected by cAMP and Ca(2+) signaling and/or MAPK pathways. The divergent roles of this gene set may help reveal how the genome and regulatory pathways evolved within the rice blast pathogen and close relatives.

COVID-BIOCHIP: A Web Tool to Analyse COVID-19 Antigen Microarrays.

During the COVID-19 pandemic, there was a growing need to characterise the disease. A very important aspect is the ability to measure the immunisation extent, which can be achieved using antigen microarrays that quantitively measure the presence of COVID-related antibodies. A significant limitation for these tests was the complexity of manually analysing the results, and the limited availability of software for its analysis. In this paper, we describe the development of COVID-BIOCHIP, an ad-hoc web-based solution for the automatic analysis and visualisation of COVID-19 antigen microarray data results.

Differential roles of the phospholipase C genes in fungal development and pathogenicity of Magnaporthe oryzae.

Calcium plays a critical role in a variety of cellular processes in cells. However, relatively little is known about the biological effects of Ca²+ signaling on morphogenesis and pathogenesis in the rice blast fungus Magnaporthe oryzae compared to other signaling pathways. We have previously demonstrated that MoPLC1-mediated calcium regulation is important for infection-related development and pathogenicity in M. oryzae. In the present study, four genes encoding phospholipase C (PLC) isozymes (MoPLC2 to MoPLC5), which differ from MoPLC1 in their domain organization, were additionally identified. The C2 domain involved in Ca²+-dependent membrane binding is found only in MoPLC2 and MoPLC3. Detailed functional analysis using deletion mutants for MoPLC2 and MoPLC3 indicated that MoPLC2 and MoPLC3 play essential roles in development. The two deletion mutants for MoPLC2 and MoPLC3 showed reduced conidiation and a defect in appressorium-mediated penetration. Reintroduction of the genes restored defects of ΔMoplc2 and ΔMoplc3. Notably, ΔMoplc2 and ΔMoplc3 mutants developed multiple appressoria on separate germ tubes of a conidium, indicating that MoPLC2- and MoPLC3-regulated signaling suppresses a feedback loop of a pathway for appressorial development. The similarity in phenotypic defects between the two mutants indicates that both MoPLC2 and MoPLC3 are important for regulation of appropriate levels of signaling molecules in a similar manner. Comparative analysis indicated that the two MoPLCs-mediated signaling pathways have interrelated, but distinct, roles in the development of M. oryzae.

An Integrated Systems Biology Approach Identifies the Proteasome as A Critical Host Machinery for ZIKV and DENV Replication.

The Zika virus (ZIKV) and dengue virus (DENV) flaviviruses exhibit similar replicative processes but have distinct clinical outcomes. A systematic understanding of virus-host protein-protein interaction networks can reveal cellular pathways critical to viral replication and disease pathogenesis. Here we employed three independent systems biology approaches toward this goal. First, protein array analysis of direct interactions between individual ZIKV/DENV viral proteins and 20,240 human proteins revealed multiple conserved cellular pathways and protein complexes, including proteasome complexes. Second, an RNAi screen of 10,415 druggable genes identified the host proteins required for ZIKV infection and uncovered that proteasome proteins were crucial in this process. Third, high-throughput screening of 6016 bioactive compounds for ZIKV inhibition yielded 134 effective compounds, including six proteasome inhibitors that suppress both ZIKV and DENV replication. Integrative analyses of these orthogonal datasets pinpoint proteasomes as critical host machinery for ZIKV/DENV replication. Our study provides multi-omics datasets for further studies of flavivirus-host interactions, disease pathogenesis, and new drug targets.

A putative MAP kinase kinase kinase, MCK1, is required for cell wall integrity and pathogenicity of the rice blast fungus, Magnaporthe oryzae.

Insertional mutagenesis of Magnaporthe oryzae led to the identification of MCK1, a pathogenicity gene predicted to encode mitogen-activated protein kinase kinase kinase (MAPKKK) homologous to BCK1 in Saccharomyces cerevisiae. Targeted disruption of MCK1 resulted in the fungus undergoing autolysis and showing hypersensitivity to cell-wall-degrading enzyme. The mck1 produced significantly reduced numbers of conidia and developed appressoria in a slightly retarded manner compared with the wild type. Appressorium of the mck1 mutant was unable to penetrate into plant tissues, thereby rendering the mutant nonpathogenic. Cytorrhysis assay and monitoring of lipid mobilization suggested that the appressorial wall was altered, presumably affecting the level of turgor pressure within appressorium. Furthermore, the mck1 mutant failed to grow inside plant tissue. Complementation of the mutated gene restored its ability to cause disease symptoms, demonstrating that MCK1 is required for fungal pathogenicity. Taken together, our results suggest that MCK1 is an MAPKKK involved in maintaining cell wall integrity of M. oryzae, and that remodeling of the cell wall in response to host environments is essential for fungal pathogenesis.

Agrobacterium tumefaciens-mediated transformation in Colletotrichum falcatum and C. acutatum.

Agrobacterium tumefaciens-mediated transformation (ATMT) is becoming an effective system as an insertional mutagenesis tool in filamentous fungi. We developed and optimized ATMT for two Colletotrichum species, C. falcatum and C. acutatum, which are the causal agents of sugarcane red rot and pepper anthracnose, respectively. A. tumefaciens strain SK1044, carrying a hygromycin phosphotransferase gene (hph) and a green fluorescent protein (GFP) gene, was used to transform the conidia of these two Colletotrichum species. Transformation efficiency was correlated with cocultivation time and bacterial cell concentration and was higher in C. falcatum than in C. acutatum. Southern blot analysis indicated that about 65% of the transformants had a single copy of the T-DNA in both C. falcatum and C. acutatum and that T-DNA integrated randomly in both fungal genomes. T-DNA insertions were identified in transformants through thermal asymmetrical interlaced PCR (TAIL-PCR) followed by sequencing. Our results suggested that ATMT can be used as a molecular tool to identify and characterize pathogenicity-related genes in these two economically important Colletotrichum species.

Zika-Virus-Encoded NS2A Disrupts Mammalian Cortical Neurogenesis by Degrading Adherens Junction Proteins.

Zika virus (ZIKV) directly infects neural progenitors and impairs their proliferation. How ZIKV interacts with the host molecular machinery to impact neurogenesis in vivo is not well understood. Here, by systematically introducing individual proteins encoded by ZIKV into the embryonic mouse cortex, we show that expression of ZIKV-NS2A, but not Dengue virus (DENV)-NS2A, leads to reduced proliferation and premature differentiation of radial glial cells and aberrant positioning of newborn neurons. Mechanistically, in vitro mapping of protein-interactomes and biochemical analysis suggest interactions between ZIKA-NS2A and multiple adherens junction complex (AJ) components. Functionally, ZIKV-NS2A, but not DENV-NS2A, destabilizes the AJ complex, resulting in impaired AJ formation and aberrant radial glial fiber scaffolding in the embryonic mouse cortex. Similarly, ZIKA-NS2A, but not DENV-NS2A, reduces radial glial cell proliferation and causes AJ deficits in human forebrain organoids. Together, our results reveal pathogenic mechanisms underlying ZIKV infection in the developing mammalian brain.