Posttranslational Arginylation Enzyme Arginyltransferase1 Shows Genetic Interactions With Specific Cellular Pathways in vivo

Arginyltransferase1 (ATE1) is a conserved enzyme in eukaryotes mediating posttranslational arginylation, the addition of an extra arginine to an existing protein. In mammals, the dysregulations of the ATE1 gene (ate1) is shown to be involved in cardiovascular abnormalities, cancer, and aging-related diseases. Although biochemical evidence suggested that arginylation may be involved in stress response and/or protein degradation, the physiological role of ATE1 in vivo has never been systematically determined. This gap of knowledge leads to difficulties for interpreting the involvements of ATE1 in diseases pathogenesis. Since ate1 is highly conserved between human and the unicellular organism Schizosaccharomyces pombe (S. pombe), we take advantage of the gene-knockout library of S. pombe, to investigate the genetic interactions between ate1 and other genes in a systematic and unbiased manner. By this approach, we found that ate1 has a surprisingly small and focused impact size. Among the 3659 tested genes, which covers nearly 75% of the genome of S. pombe, less than 5% of them displayed significant genetic interactions with ate1. Furthermore, these ate1-interacting partners can be grouped into a few discrete clustered categories based on their functions or their physical interactions. These categories include translation/transcription regulation, biosynthesis/metabolism of biomolecules (including histidine), cell morphology and cellular dynamics, response to oxidative or metabolic stress, ribosomal structure and function, and mitochondrial function. Unexpectedly, inconsistent to popular belief, very few genes in the global ubiquitination or degradation pathways showed interactions with ate1. Our results suggested that ATE1 specifically regulates a handful of cellular processes in vivo, which will provide critical mechanistic leads for studying the involvements of ATE1 in normal physiologies as well as in diseased conditions.


INTRODUCTION
Protein posttranslational modifications (PTM) change protein properties without requiring de novo synthesis. Thus, PTMs are frequently relied upon to respond to acute stress or cellular signaling, and are often used to activate lateral response factors such as transcription or epigenetic modulations.
For these reasons, dysregulation of PTMs are often indicated in cardiovascular diseases where stress response or cellular signaling play critical roles.
In eukaryotes, many proteins appear to be subjected to N-terminal arginylation (referred to as arginylation in this proposal), a ribosome-independent addition of one extra arginine on the N-terminus of a protein. Arginylation is catalyzed by the family of arginyltransferase (ATE). While plants contain ATE1 and ATE2, fungi and metazoans only contain ATE1, which is highly conserved across different species. Multiple lines of genetic studies have shown an important role of the ATE1 gene (ate1) in the cardiovascular and/or metabolic diseases in animals. For example, a genomic deletion of ATE1 was found to result in complete embryonic lethality in mice during the mid-gestation stage (E9-14) with severe defects in angiogenetic remodeling and cardiac development (Kwon et al., 2002). Moreover, organ/tissue-specific deletions of ATE1 in animals were found to cause a variety of abnormalities in morphology and function. Examples include a progressive dilated cardiomyopathy in mice when ate1 knockout is driven by the cardiac myosin heavy chain promoter (Kurosaka et al., 2010;Kaji and Kaji, 2012;Wang et al., 2017b). Also, inducible systematic deletion of ate1 appears to cause rapid weight loss, damaged spermatogenesis, neurological perturbations, and early lethality in adult mice (Brower and Varshavsky, 2009). In addition to these involvement in cardiovascular and metabolic abnormalities, a dysregulation of ATE1 is indicated in cancer as well. For example, reports from us and other groups showed that ATE1 is often downregulated in high grade cancer cases and is associated with poorer outcomes (Rai et al., 2015;Birnbaum et al., 2019), and that an inhibition of ATE1-mediated arginylation confers cancer cell resistance to apoptosis-induced by radiation (Masdehors et al., 2000). Furthermore, mounting evidence is also starting to indicate the involvement of ATE1 in aging-related conditions (Brower et al., 2013;Wang et al., 2017a). Unfortunately, the physiological role of ATE1 (and its arginylation activity) remains poorly understood, which adds to the difficulty of interpreting its involvements in normal conditions or diseases.
The studies about ATE1 and arginylation are still relatively scarce. Our understandings for these topics are being continuously reshaped with emerging new evidence. Arginylation has been found to take place on nearly a hundred eukaryotic proteins and the list is still expanding on a daily basis (Decca et al., 2006;Wong et al., 2007;Piatkov et al., 2012). Considering that a wide range of proteins are substrates of arginylation, it is reasonable to speculate that ATE1 may act as a root regulator of multiple cellular processes.
A popular theory about arginylation is its involvement in global protein degradation. Arginylation was shown to promote hyper-ubiquitination of the substrate proteins, which is then shown to be degraded by proteasome or autophagy (Saha and Kashina, 2011;Varshavsky, 2011). Based on mostly artificial substrates and in vitro data, arginylation was proposed to take place on proteins bearing certain amino acids on the 2nd residue on the N-terminus. These include the amino acids Asp, Glu, Asn, Gln (in fungi and mammals), and Cys (in mammal only) (Varshavsky, 2011). By this rule, in any given eukaryotic organism, at least 20-25% of its proteome would be estimated to be degraded by arginylation. Based on this assumption, arginylation was proposed as a central component in the generic protein degradation machinery (Varshavsky, 2011). However, the exact impact of arginylation on global homeostasis of proteins in vivo remain undetermined. In a more recent report, based on comparison of the sizes of protein spots on 2D-gels, the knockout of ate1 in mouse cells appears to affect ∼20% of those spots on the steady-state levels. However, much of these effects appear to be derived from transcriptional changes. Also, proteasome inhibitors can only reverse the ate1-dependent reduction of less than 3% of the observed protein spots (Wong et al., 2007). Such a small impact size is inconsistent with the proposed role of arginylation as a global degradation machinery.
In addition to protein degradation, arginylation was also suggested to be involved in many other processes. These may include cellular response to various types of stressors such as those that are closely related to cardiovascular stresses. For example, the activity of arginylation in cells or animal tissues is altered during injury, high temperature, or exposures to high concentrations of oxidants or salt (Zanakis et al., 1984;Chakraborty et al., 1986;Shyne-Athwal et al., 1986, 1988Luo et al., 1990;Jack et al., 1992;Chakraborty and Ingoglia, 1993;Xu et al., 1993;Wang and Ingoglia, 1997;Bongiovanni et al., 1999;Kumar et al., 2016). Furthermore, other lines of studies also suggested that ATE1 may regulate the dynamics of cytoskeleton (Karakozova et al., 2006;Saha et al., 2010Kaji and Kaji, 2012).
While multiple lines of researches are starting to establish ATE1 as a root regulator for multiple processes in the cell, contradicting conclusions are very often seen among reports from different groups. While many of these discrepancies may arise from differences in test conditions, variations may also be generated because many past studies were focused on the role of individual arginylation substrates without considering the effects of arginylation on the other known or potential substrates. To further understand the function of ATE1 in normal or diseased conditions, it is desirable to investigate the physiological role of ate1 in a systematic manner. Among the currently available tools for studying functional genomics, approaches based on a yeast gene-deletion library remain as one of the most robust and straight-forward methods (Boone et al., 2007;Giaever and Nislow, 2014). However, the role of ATE1 gene (ate1) has never been studied with such an unbiased manner, either as a specific query subject or as part of a comprehensive screening.
In this study, we took advantage of the Schizosaccharomyces pombe (S. pombe) single-gene knockout library. S. pombe is a commonly used test model for eukaryotic genes due to its similarity to metazoan organisms, while preserving its ease of genetic operation as a microbe. The meiosis and mating process of this organism also greatly facilitate the combination of different gene deletions to examine their synthetic effects, which can be quantitated by monitoring the growth rates of the resulting cells to deduce the genetic interactions between these two genes (Spirek et al., 2010;Wiley et al., 2014). The advantage of S. pombe as a model system is further demonstrated by the fact that nearly 70% of its gene have orthologs in the human genome, which is higher than S. cerevisiae, the other commonly used yeast as test model (Yanagida, 2002;Hoffman et al., 2015). The conserved genes including ate1, which is nearly 70% homologous in amino acid sequence in the core domain (∼200 residues) compared to its counterpart in human. However, unlike mammalian cells, S. pombe allows > 75% of its genome to be individually knocked out for functional tests (Spirek et al., 2010). These facts made S. pombe a desirable test model to determine the physiological role of ate1 in vivo.
In this study, we examined the effects of combining ate1deletion with the deletions of other genes in S. pombe. The library we employed covered 3659 genes in S. pombe, which account for nearly 75% of the predicted open reading frames (up to 4940) in this organism (Decottignies et al., 2003). Surprisingly, only 173 of these genes, which is less than 5% of the effective library size, showed significant genetic interactions with ate1. Furthermore, many of these genes can be clustered into a few discrete groups in relation to cellular pathways. These include ribosomal component and translation regulation, gene transcription, oxidative stress response, cytoskeletal/structural components, mitochondrial function, and synthesis/metabolism of organic molecules and amino acids (such as histidine). Also, unexpectedly no significant interactions were observed between ate1 and genes involved in the ubiquitin/proteome system (UPS). Our results indicate that ATE1 may specifically regulate several physiological pathways in vivo. Many of these interactions can be used to provide satisfactory explanations for many observed involvements of ATE1-mediated arginylation in cardiovascular/metabolic abnormalities and other diseased conditions. These novel findings will also provide clues for designing approaches intervening ATE1-related phenotypes in various diseases including cancer and metabolic dysregulations.

RESULTS
The Arginyltransferase Gene ate1 Showed Interactions With Only a Small Subset of Genes While ate1 is an essential gene for mammals, it can be knocked out in S. pombe without causing lethality. Although S. pombe usually exists in a haploid form and reproduces by symmetric division (fission), they can also be induced to perform mating, during which chromosome recombination proceeds in a relatively high rate. These unique properties allow the knockout of ate1 to be easily combined with the knockout of other individual genes in the library, provided that different selection markers are being used to trace the knockout library and the query knockout (of ate1). Based on the growth of the products, compared to the parental strains, at least two types of interactions can be measured: (1) phenotype-enhancement, in which the crossing of the query results in a lower growth rate (also referred to as synthetic lethal); (2) phenotype-suppression, with a faster growth rate (also known as synthetic suppression).
By using the above-mentioned synthetic knockout approach, we examined the effects of combining ate1-knockout with 3721 individual knockouts in a S. pombe library, in which most of these genes were functionally annotated by either experimental evidence or prediction based on known orthologs. After excluding knockouts that lead to sterilization of the yeast, totally 3659 effective crossings were examined. Out of these crossings, we only found 173 of them resulted in significant phenotype enhancement (see Table 1A) or suppression (see Table 1B). This number is less than 5% of the library size.

Genes That Genetically Interact With ate1 Can Be Clustered by Biological Function
Among the genes that showed significant interactions with ate1, we found specific enrichment of several types of functional relevancies. For example, simply by applying functional annotations such as gene ontology (GO), gene expression category, or Fission Yeast Phenotype Ontology (FYPO), we found that several terms are enriched or present in higher frequencies among genes interacting with ate1 compared to the whole library (Figure 1 and Supplementary Table S1). As a further validation of the functional interactions observed between ate1 and these genes, many of them also have known physical interactions or associations of their products among each other (Figure 2).
The enriched or increased functional annotations can be grouped into several clustered categories. Based on the frequency of presence in the hit list, the most abundant genes (nearly 50% of the hits) are those related to biosynthesis of biological molecules or proteins (Supplementary Table S1A). These include many regulators of translation or transcription (17-39%). Interestingly, among these categories, those related to histidine/amidazole synthesis/metabolism appear to be particularly impacted by ate1, as most of the genes associated with these pathways (5 out of 7∼9) showed interaction with ate1 (Figure 1 and Supplementary Table S1A).
The second most abundant genre (up to 40%) is those related to cell morphology as well as cellular dynamics (fusion, conjugation, or division). These include many known regulators of cytoskeleton (Figure 1 and Supplementary Table S1B).
The third type (10∼15%) are those related to oxidative stress response, which include many highly expressed genes coding for chaperones or redox regulators (Figure 1 and Supplementary Table S1C).
The fifth are those related to nutrient or metabolite stresses, such as genes related to nitrogen starvation, and genes involved in responses to metabolism interferers (caffeine, rapamycin) or metabolite analogs of lipids or nucleotides (Figure 1 and Supplementary Tables S1B,C).
In addition to the above analysis by functional annotations, the genes that have genetic interactions with ate1 can also be grouped into a handful of clustered functional pathways or categories when being analyzed with PANTHER Classification The genetic interaction was determined by the effects of combing two gene knockouts on the growth rate of the resulting yeast colony compared to the individual knockouts. Sexual mating of S. pombe was used to generate the crossings between the ate1-knockout and the other gene-knockouts in a library. While the original library contains 3721 different knockout strains, some of these deletion strains appear to be sterile and cannot be mated. Thus, the effective library size was 3659 in our tests. Four replicates were used to determine the effects on growth of the crossed, double-knockout colony compared to the parental strain with a single knockout (without ate1 deletion). Based on the comparison, a result of slower growth is considered as phenotype-enhancement (in Table 1A). Vice versa, those that grow faster than the parental strain is considered phenotype-suppression (Table 1B). A confidence value (p-values) of 0.05 is the minimum for significance.
Frontiers in Physiology | www.frontiersin.org System 1 (Mi et al., 2013(Mi et al., , 2019, as shown in Supplementary Table S2. As an example to validate the observations reached in the above non-biased study, we examined the sensitivity of S. pombe to exogenous histidine, which is not an essential amino acid for this organism. While a lower amount of exogenous histidine usually promotes the growth of yeasts, at high concentrations it is known to generate cytotoxicity and requires the actions of the histidine/amidazole synthesis/metabolism pathway for mitigation (Winkler and Ramos-Montanez, 2009;Watanabe et al., 2014;Duncan et al., 2018). When we compared the growth rates of the yeasts hosting the ate1-deletion to the control cells, we found that, although there is no apparent difference in their proliferation in the absence of exogenous histidine 1 www.pantherdb.org ( Figure 3A), ate1-deleted yeast grow much slower than the control in high concentrations of histidine (Figures 3B-E). To ensure these observations are not specific to the colony of yeast used in the test, we repeated these experiments with two additional clones of the ate1-deleted S. pombe and reached similar conclusions ( Figure 3F). Therefore, these data suggested that ate1-deletion indeed possesses an interaction with the histidine/amidazole synthesis/metabolism pathway.

Many Genes Related to Mitochondria or Energy Production Genetically Interact With ate1
Compared to the relative abundance of literatures showing the effects of arginylation in protein degradation or cytoskeletal regulations, the potential involvements of FIGURE 2 | Known protein-protein interactions (PPI) between products of genes that have genetic interactions with ate1 in S. pombe. The PPI database being utilized is Version 11 of STRING (https://string-db.org). The interaction map was generated with Cytoscape using the STRING App for creating the PPI information and enrichment mapping. The thickness of the connecting line represents the confidence of the experimental protein-protein interaction. A thicker line represents a higher confidence (a minimum of 0.4% confidence was used). Each rounded shape represent the product of a gene (with gene name labeled on the side). The assigned color of the rounded shape represents gene category (by GO terms), while the color of the ring (red or green) represent the direction of the genetic interaction (phenotype enhancing or suppressing). arginylation or ATE1 in other cellular processes are less known. Interestingly, many of the observed phenotypes in animals resulted from ate1 dysregulation may be at least partly attributed to a disruption of mitochondrial function, which supplies the majority of energy molecules (ATP, NADH, and NADPH, etc.) required for many biosynthesis processes, and also constitute a major source of oxidative stressors in the cell. For example, postnatal systemic knockout of ate1 appears to lead to drastic loss of fat and infertility, which are common consequences expected from a compromise of mitochondrial function. Consistent to this possibility, we observed a significant FIGURE 3 | The deletion of ate1 leads to lower tolerance of exogenous histidine in the media. The yeast strains were grown in PMG liquid media supplemented with adenine, leucine, and uracil plus histidine as indicated (no H = no histidine; 1/2H = 112 µg/ml histidine; H = 225 µg/ml histidine), while the culture density was monitor by absorbance (OD 595 nm). In the absence of histidine (referred as to "no H"), both the ate1 and the control (ate1 +) strains grow in a similar rate (see A). The addition of moderate amount of histidine (112 µg/ml; "1/2H") increases the growth of control strain (ate1 +), while a higher concentration (225 µg/ml; "H") exhibits cytotoxic effect as anticipated (see B). However, in ate1 strain, the addition of either 112 or 225 mg/ml of exogenous histidine both lead to slower growth in a dose-dependent manner (see C). Direct comparisons of the growth of ate1 and the control strain (ate1 +) at different concentrations of histidine are presented at panel (D) and (E). To exclude the possibility of clone-specific, two different clones of ate1 strain (ate1 −C2, and −C3) isolated from the same knockout process were subjected to the challenge of 112 µg/ml exogenous histidine, compared to when no histidine was added (see F). The curves were generated with non-linear regression. Error Bars represent standard deviation from 6 replicates.
number (at least 19) of genes related to mitochondria or energy production showing interactions with ate1 ( Table 2). These genes represent more than 10% of the total hits, although they were not present as "enriched" because the frequency of them in the hit-list is not higher than that in the library. The numbers of those genes in the effective library (non-sterile) and the hit list are being shown. The names (or systematic ID) of the genes in the hit list are also shown.

Few Genes in Global Protein Degradation Pathways Showed Genetic Interactions With ate1
Another surprising finding in the interacting partners of ate1 is the rarity of genes involved in global protein ubiquitination or degradation machineries. In principle, the uncertainties about the role of arginylation in proteome homeostasis can at least be partly addressed by the genetic interaction screening. If ate1 is involved in global protein ubiquitination or degradation, it would be expected to have genetic interactions with the other components in these pathways. However, while many genes with known roles in global ubiquitination or degradation are present in the S. pombe knockout library we employed (Supplementary Table S3), only two (ubi5 and atg14) showed significant interactions (either enhancing or suppressing) with ate1 (Table 3). Among these two genes, ubi5 is a fusion gene of ubiquitin and ribosomal component. As such, its genetic interaction with ate1 may derive from the ribosomal component and not necessarily the ubiquitin, since no interaction was observed between ate1 and other genes coding for ubiquitin. The rarity of interactions between ate1 and degradation-related genes is actually consistent with the selectivity of the interacting partners of ate1 as described above (Figure 2 and Tables 1-3). If arginylation is a generic degradation pathway for up to 20% proteins as predicted, then the expected interactions between ate1 and other genes should be much larger than the observed number (< 5%), because many components in the pathways regulated by these arginylation-target proteins would be expected to show genetic interaction with ate1. Interestingly, the small impact size (< 5%) of ate1 in genetic interactions is well consistent with the impact size (∼3%) of arginylation on the degradation of endogenous proteins as estimated in 2D-gels (Wong et al., 2007). These evidences suggest that ATE1-mediated arginylation may not be a major degradation pathway in vivo in the experimental condition we employed.

DISCUSSION
Despite that ate1 gene has been identified for 30 years, it's in vivo role has remained poorly defined. For the first time, by using a systematic approach, our results showed that ate1 possesses significant genetic interactions with a small and focused subset of genes concerning multiple critical cellular processes. The results The numbers of those genes in the effective library (non-sterile) and the hit list are being shown. The names of the genes in the hit list are also shown.
from this study can also provide important leads for mechanistic investigations about the role of ate1 or arginylation in normal or diseased conditions. The power of the systematic approach employed in our study is demonstrated by the fact that many results from this unbiased investigation are highly consistent or complementing to data from past reports. For example, the observed genetic interactions between ate1 and regulators of cell morphology are well consistent with previous reports showing the impact of arginylation on many cytoskeletal proteins (Wong et al., 2007) and cytoskeletal dynamics (Karakozova et al., 2006;Rai et al., 2008;Saha et al., 2010;Zhang et al., 2010Zhang et al., , 2012Kurosaka et al., 2012). The interactions between ate1 and genes involved in oxidative stress response also support the proposed role of arginylation in stress response process (Zanakis et al., 1984;Chakraborty et al., 1986;Shyne-Athwal et al., 1986, 1988Luo et al., 1990;Jack et al., 1992;Chakraborty and Ingoglia, 1993;Xu et al., 1993;Wang and Ingoglia, 1997;Bongiovanni et al., 1999;Kumar et al., 2016). Overall, it appears that our non-biased screening study was able to recapitulate many past observations performed in cell or animal.
Many results of this study will provide new clues for investigate the role of arginylation in physiological processes. Arginylation was shown to be involved in stress response, but the exact mechanism still awaits clarification and may benefit from the genetic interactions revealed in our study. For example, while several redox regulators, including manganese superoxide dismutase (SOD) appear to have phenotypeenhancing relationship with ate1, many genes related to mitochondria, the main source of reactive oxygen species (ROS) in the cell, display synthetic suppression relationship. These differences imply that ATE1 may act as a scavenger of oxidatively damaged proteins in stress response. Also, the broad interactions between ate1 and many transcription regulators and histone modulators may help to elucidate the observed but unexplained effects of ATE1 on global transcriptional landscape (Lee et al., 2012;Eisenach et al., 2014;Deka et al., 2016). Particularly, many genes of histone modulators appear to have a phenotypesuppression relationship with ate1. This is highly intriguing considering that previous studies have found that histone proteins are subjected to arginylation modification, which may also affect the other PTMs on histone (Wong et al., 2007;. The interactions between ate1 and genes related to biomolecule synthesis/metabolism are also of high interest because emerging evidence indicated a potential involvement of ATE1 in metabolism. As an example, the enrichment of the genes in histidine synthesis/metabolism pathways in principle is consistent with one of our previous findings about arginylation of phosphoribosyl pyrophosphate synthetase (PRPPS). This is because PRPPS is responsible for production of Phosphoribosyl Diphosphate, which serves as a precursor substrate for biosynthesis of histidine . Due to the conserved nature of ATE1-mediated arginylation in eukaryotes (McGary et al., 2010), many of these findings may provide mechanistic insights for the role of ATE1/arginylation in cardiovascular diseases, metabolic dysregulations, and cancer in human. These exciting possibilities will become important directions for future studies.
The lack of interactions between ate1 and genes related to global ubiquitination or protein degradation, as revealed in our study, is quite surprising. Based on the loose consensus on amino acid sequence of known substrates, arginylation was long hypothesized as a signal for the ubiquitination and degradation of at least 20-25% members of the proteome in yeast or metazoan. Such a broad range of substrates would predict a very large impact size of genetic interactions of ate1. However, less than 5% of the tested genes, which cover nearly 75% of the genome in S. pombe yeast, showed interactions with ate1. Consistently, we found that ate1 has very few interactions with components of the ubiquitination and degradation pathways. While these findings appear to be at odds with popular theories, it is worthy pointing out that the exact nature of arginylation in the degradation of endogenous proteins has not been decisively determined and the current theories about arginylation were mainly built on studies with artificial substrates (Bachmair et al., 1986;Varshavsky, 2011). It is possible that the substrate preference of ATE1 is more complexed than originally expected. For example, recent evidence suggested that the efficiency of arginylation may be affected by at least 11 residues on the N-terminus of a peptide, and the substrate preference of ATE1 may also be influenced by additional in vivo factors (Wang et al., 2011(Wang et al., , 2018. Also, since the majority of arginylation takes place on the N-terminus of a protein, it must compete with N-terminally acetylation in vivo, which is a dominant PTM in most eukaryotes (Wang et al., 2011). Furthermore, while arginylation in many cases indeed promotes ubiquitination and degradation, exceptions are also abundant on endogenous arginylated proteins (Karakozova et al., 2006;Zhang et al., 2010Zhang et al., , 2012Zhang et al., , 2015. Maybe not coincidentally, past experimental attempts examining the impact of arginylation on the whole proteome also revealed a relatively small size. In one of such reports, the usage of 2D gels showed that less than 3% of individual proteins appear to be affected by arginylation on proteasome-dependent degradation (Karakozova et al., 2006;Wong et al., 2007). As such, alternative interpretation for the role of arginylation in protein ubiquitination/degradation may be needed and our study may provide clues for that. It is possible that arginylation is not activated under resting state but is specifically utilized for protein degradation during certain conditions such as stress response or nutrient deprivation. This possibility is at least indirectly supported by existing evidence showing the activation of arginylation during stress response and that arginylation preferentially takes place on oxidatively damaged proteins (Zhang et al., 1998;Kumar et al., 2016). It is also supported by the extensive genetic interactions between ate1 and genes related to oxidative and metabolic stress as revealed in our study. However, this current study was conducted in a non-stressed condition and therefore cannot directly test this possibility, which will require further investigations for validation.
The other unexpected finding is the extensive interaction between ate1 and ribosome-related genes. While the biological meaning of this interaction is still unclear, existing literatures may provide a few clues. Particularly, the ATE1 protein was found to at least partially co-localize with ribosome during cellular fractionation (Wang et al., 2011). Also, many ribosomeassociated proteins were found to be arginylated (Wong et al., 2007;Wang et al., 2011). Based on this evidence, it is likely ATE1 may regulate the function of ribosome. In addition, current evidence suggest that arginylation may take place during co-translational stage and the outcome in arginylation-mediated degradation may be dependent on the dynamic of co-translational folding . As such, it is also likely that arginylation is part of the mechanism for quality check of nascent peptide synthesis. This unexpected connection between arginylation and ribosome, as well as many other newly information about arginylation uncovered in this study, constitute intriguing directions for future research endeavors.

CONCLUSION
By using a systematic approach, we found that the gene for arginyltransferase1 has specific interactions with a small subset of genes in the eukaryotic genome, which fall into a few clustered functional categories. Our data suggest that ATE1 may specifically regulate a few cellular pathways in vivo. These results will provide novel mechanistic clues to understand the role of protein arginylation in the development of cardiovascular system and the pathogenesis of related diseases.

Query Yeast Strain Creation
Growth conditions and genetic manipulations of S. pombe were performed as previously described (Moreno et al., 1991). The query strain containing ate1-deletion (ate1 ) was prepared by targeted deletion with a linear DNA containing a nourseothricinresistance gene flanked with sequence derived from the ate1 gene (systemic ID: SPAC3C7.07c) loci. The 5 -region (394 bp) flanking the resistance gene was cloned with these primers, with underline to indicate the regions matching the genome sequence while the rest primes to the resistance cassette: 5 -FWD: TAGAACTTGGTGGATGGTATCGTGG 5 -REV: GGGGATCCGTCGACCTGCAGCGTACGAACTAT TGTTTGAAAATTTCCCTGTTTAAT The 3 -region (385 bp) flanking the resistance gene was cloned with these primers:

Yeast Synthetic Genetic Array
The ate1 strain was grown to exponential phase in liquid PMG media (Sunrise Scientific Cat. #2060) supplemented with adenine, leucine and uracil. It was then pinned to agar plates in a 384matrix format using a robotic platform RoToR HDA (Singer Instruments). The query strain was then crossed with individual strains in the S. pombe haploid deletion library containing individual gene deletions marked by a kanamycin resistance gene (Bioneer, version 4.0 equivalent) utilizing a modified SGA procedure (Dixon et al., 2008). This procedure was described in detailed in our previous publication (Wiley et al., 2014). In brief, the crossing was induced on a SPAS mating media (1% glucose, 7.3mM KH 2 PO 4 , with 45mg/L supplements of adenine, histidine, leucine, uracil and lysine-HCl, and vitamin supplement for pantothenic acid, nicotinic acid, inositol, and biotin; see this website for detailed recipe: https://dornsife.usc. edu/pombenet/media/). For germination, four replicates of each crossing product were pinned to a 1536 format on selective PAU + G418 media, which is the PMG media containing adenine (225 mg/L, Sigma Cat. #A8751), leucine (225 mg/L, Sigma Cat. #L8912), uracil (225 mg/L, Sigma Cat. #U0750), and the antibiotics G418 and nourseothricin. Colony growth was monitored for 3 days utilizing a flatbed scanner, and then quantified and compared using ScreenMill according to published protocols (Dittmar et al., 2010). Double mutants with a significant growth rate difference compared to the corresponding library gene deletion alone were scored as either slowed (referred to as "phenotype-enhancement") or accelerated (referred to as "phenotype-suppression") growth.

Hit Analysis and Bio-Informatics
The gene feature enrichments (over-representation) in the hit list versus the library were examined by AnGeLi 2 (Bitton et al., 2015). The hit list was compared directly to the genes screened in our assay using the false discovery rate setting for multiple testing with a p-value setting of 0.01. In addition, for this enrichment we performed pairwise interaction enrichment with 1000 permutations and allowed for the p-value to adjust. As a negative control, we also looked for underrepresented terms and we cannot find significant underrepresented entries based on inputs, which suggesting that our analysis did not create artifacts.
The image of protein-protein interaction was created using Cytoscape with the STRING App (Szklarczyk et al., 2019). The PPI database being utilized is Version 11 of STRING 4 .

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation, to any qualified researcher.

AUTHOR CONTRIBUTIONS
DW performed most of the experiments and the analysis of data. GD'U was involved in supervising the performance of the experiments. FZ conceptualized and designed the project and also supervised the analysis of the data.

FUNDING
This project was partly supported by a pilot funding from the Florida Breast Cancer Foundation, and an American Cancer Society Institutional Research Grant 98-277-13, to FZ.