article: The functions of the majority of trypanosomatid-specific proteins are unknown, hindering our understanding of the biology and pathogenesis of Trypanosomatida. While protein-protein interactions are highly informative about protein function, a global map of protein interactions and complexes is still lacking for these important human parasites. Here, benefiting from in-depth biochemical fractionation, we systematically interrogated the co-complex interactions of more than 3354 protein groups in procyclic life stage of Trypanosoma brucei, the protozoan parasite responsible for human African trypanosomiasis. Using a rigorous methodology, our analysis led to identification of 128 high-confidence complexes encompassing 716 protein groups, including 635 protein groups that lacked experimental annotation. These complexes correlate well with known pathways as well as for proteins co-expressed across the T. brucei life cycle, and provide potential functions for a large number of previously uncharacterized proteins. We validated the functions of several novel proteins associated with the RNA-editing machinery, identifying a candidate potentially involved in the mitochondrial post-transcriptional regulation of T. brucei. Our data provide an unprecedented view of the protein complex map of T. brucei, and serve as a reliable resource for further characterization of trypanosomatid proteins. The presented results in this study are available at: www. TrypsNetDB. org.
  The unicellular flagellated parasite Trypanosoma brucei is the causative agent of the fatal human African trypanosomiasis (HAT), also known as sleeping sickness, and the economically devastating disease Nagana in cattle. The life cycle of T. brucei comprise an insect stage (procyclic stage) and a mammalian stage. Trypanosomatid parasites, including T. brucei, are highly diverged from well-studied eukaryotes such as yeast and mammals [1–3], resulting in a lack of sequence similarity with model organisms for the majority of their proteins. This has led to inapplicability of homology-based approaches for identification of potential functions of more than half of all trypanosomatid genes [4]. Protein interaction maps offer an invaluable resource for functional annotation of proteins [5]. Current methodological/instrumental advances have led to the development of several ex vivo, in vivo, and in silico approaches to systematically chart protein interactions and complexes [6]. Optimized yeast two-hybrid (Y2H) approaches have been employed to infer pairwise interactions among proteins [7,8]. Immunoprecipitation [9], biochemical fractionation [10–13], and affinity purification (AP) -based approaches [14–16] are widely used for the identification of protein complexes in a specific cell context. Additionally, functional association of proteins can be predicted computationally using data types such as transcriptomics data [17], synthetic lethality [18], and chemical sensitivity [19]. However, each of these approaches has limitations and is inherently associated, in varying degrees, with false positive and negative results. In AP-based approaches, for example, tagging the protein may affect the binding partners of the tagged protein by inactivation, capping the binding site, or changing the localization of the protein. Highly expressed proteins are also often co-purified with the tagged protein as a false positive contaminant. Moreover, transient interactions are likely to be lost if stringent conditions are used for the purification of the tagged protein. In biochemical fractionation strategies, fortuitous interactions can arise because confounding protein complexes can still be present in the same fraction regardless of in-depth fractionation [10]. In addition to non-negligible false positive rates, the Y2H system is relatively weak at detection of co-complex associations, although it works well at capturing binary, particularly transient, interactions [20]. Therefore, the integration of data from different approaches has been shown to improve the precision of protein maps [21]. To explore the protein complexes underlying the survival and pathogenesis of T. brucei, we performed four high resolution fractionation experiments benefiting from two orthogonal, complementary biochemical approaches. High-resolution mass spectrometry analysis of the fractions led to the construction of a global co-fractionation network for T. brucei in procyclic life stage. Evaluation of the constructed network demonstrated that it has topological and biological characteristics that are similar to those observed in the sampled networks of model organisms from previous large-scale studies. Importantly, our results demonstrated significantly higher precision for those interactions that were supported by the two orthogonal fractionation approaches compared to those that co-fractionated only in one approach. To extract a high-confidence core network, we combined the fractionation-derived network with other orthogonal resources of protein-protein interaction data. This high-confidence network predicts the network context of 866 protein groups, including many hypothetical and experimentally unannotated proteins. Clustering of this high-confidence network led to the assignment of 716 protein groups to 128 complexes. To our knowledge, this study presents the largest proteomics-based interaction map of trypanosomatid parasites to date, providing a useful resource for formulating new biological hypothesises and further experimental leads. To showcase the utility of this protein complex map, we used it to reveal novel factors involved in the mitochondrial post-transcriptional regulation of T. brucei, and validated them by several independent experiments.
  Late log phase (~2 x 107 cells/ml) T. brucei procyclic form IsTaR 1. 7 A cells (a derivative of EATRO 164) were grown in 225 cm2 flasks to obtain 4×109 cells. The cells were harvested by centrifugation at 6000 x g for 10 min at 4°C, and washed once with ice-cold glucose-supplemented PBS (6mM glucose). The cells were then resuspended in 500 μl lysis buffer [10 mM Tris-HCl pH 7. 2,10 mM MgCl2,100 mM KCl, 1 μg/ml pepstatin, 1 mM DTT, 1% triton X-100 and 1 x EDTA-free protease inhibitor cocktail (Roche) ] and incubated on a tube rotator for 15 min at 4°C. The lysate was treated with 40 units of RNase-free DNase I (Roche) for 1 h on ice and cleared twice by centrifugation at 16,000x g for 15 min at 4°C. The methods for preparing the extracts were adapted from conventional purification techniques [22–24]. Cell pellets were washed with ice-cold glucose-supplemented PBS as above, resuspended in 30 ml DTE buffer [1 mM Tris-HCl pH 8. 0 and 1 mM EDTA] containing a protease inhibitor cocktail (Roche), and lysed using a 40 ml sterile tight-fitting Dounce-homogenizer on ice. The lysate was immediately made isotonic by adding a 2 M sucrose stock solution to a final concentration of 250 mM, and the mitochondria were centrifuged at 15,800 x g for 10 min at 4°C. The mitochondrial pellet was treated with 0. 3 mM CaCl2 and 40 Us of RNase-free DNase I (Roche) in 4. 6 ml STM buffer [20 mM Tris-HCl pH 8. 0,250 mM sucrose and 2 mM MgCl2] for 1 hr on ice and precipitated again. The mitochondrial lysate was prepared in 500 μl lysis buffer as described above for whole cells. The supernatant obtained upon collecting mitochondria at 15,800 x g, represents a crude preparation of the cytosol. This was centrifuged at 100,000 x g for 1 h at 4°C to eliminate small organelle contamination. Whole cell and mitochondrial extracts were resolved on 10–30% linear glycerol gradients [40 mM HEPES pH 7. 9,20 mM Mg (OAc) 2,100 mM KCl, and 2 mM EDTA] at 178,000x g for 6 hr at 4°C, and fractionated into 46–48 fractions (250 μl each), as described elsewhere [25]. Protein separation on these gradients were standardized with known amounts of BSA (Bovine serum albumin), catalase and IgM (Immunoglobulin M), with apparent masses of 66 KDa, 230 KDa and 970 KDa, respectively. The glycerol gradient fractions were analyzed using an Agilent 1100 series capillary LC system and a ThermoScientific LTQ linear ion trap mass spectrometer. Fifty μl portions of each fraction were concentrated to approximately 10 μl volume by vacuum centrifugation, and 30 μl of 8M urea, 1. 0 M Tris (pH 8. 5), 8 mM CaCl2, and 0. 2 M methylamine added. Cysteines were then reduced by 14 mM dithiothreitol (DTT) at 50°C for 15 min, and alkylated by 33 mM iodoacetamide at room temperature for 30 min. The urea concentration in the samples was then diluted to 2M concentration and trypsin (Sigma, T6567) added at an approximate ratio of 1: 25 enzyme: substrate during an overnight digestion. Digests were then acidified by addition of formic acid at a final 5% concentration, vacuum concentrated to approximately 25 μl and injected for LC-MS/MS analysis. Electrospray ionization was performed with no sheath gas using an ion max source fitted with a 34 gauge metal needle. Samples were applied at 20 μl/min to a trap cartridge (Michrom BioResources, Auburn, California, USA) and then switched onto a 0. 5×250 mm Zorbax SB-C18 column (Agilent) using a mobile phase containing 0. 1% formic acid and a 7–30% acetonitrile (ACN) gradient over 95 min, with a 10 μl/min flow rate. Peptide sequences were determined by comparing the observed MS/MS spectra to theoretical MS/MS spectra of peptides generated from a T. brucei protein database containing 9,826 entries (www. tritrypdb. org, release 3. 2) comprised of both forward and reversed sequences using the program Sequest (Version 27, rev. 12, ThermoScientific). SEQUEST was configured as follows: parent ion average mass tolerance of 2. 5 Da, fragment ion monoisotopic mass tolerance of 1. 0 Da, tryptic enzymatic cleavage, a maximum of 2 missed cleavages, static C+57 modifications, and variable M+16 modifications. Peptides were filtered at 2% false discovery rate (FDR) for whole cell cytosolic samples and 5% FDR for mitochondrial samples using the PAW pipeline [26] and the target/decoy method. It should be noted that the reported peptide-to-spectrum match (PSM) FDRs are for the entire set of fractions in each experiment; and the FDR per fraction is much less than 1% even for the mitochondrial glycerol gradient fractionation data. During the protein identification, we created protein groups by grouping the paralogous proteins (proteins that were generated by duplication events within the genome) with nearly identical protein sequences together. The total number of protein groups identified in this study was 1526, using a minimum of 2 unique peptides per protein, at a false discovery rate of about 1% for the combined fraction analysis and about 3% for the individual fraction analysis based on comparison of hits identified for forward and reverse sequences. Protein abundances were based on MS2 intensity-weighted spectral counting, which have been shown to provide accurate protein abundance estimates from LTQ linear ion trap mass spectrometer [27]. Cytosolic and mitochondrial extracts were resolved by liquid chromatography using tandem cationic (S) –anionic (Q) exchange columns (UNOsphere, Biorad), adapted from [10,28]. The mobile phase consisted of buffer A [10 mM Tris-HCl pH 7. 8,10 mM MgCl2,50 mM KCl and 1 mM DTT] and buffer B [buffer A + 950 mM KCl]. Chromatography was performed using Beckman-Coulter Gold high performance liquid chromatography system. Samples were passed through a 0. 22 μm membrane, loaded on the columns and then washed for 15 mins with buffer A, followed by 0–50% buffer B (1%/min), a 10 min wash and then 50–100% buffer B (5%/min). The flow-rate was maintained at 400 μl/min, and approximately 40 fractions (800 μl each) were collected. The ion exchange separated proteins were analyzed using a Dionex NCP-3200RS UltiMate RSLCnano UPLC system and Thermo Scientific Orbitrap Fusion mass spectrometer. Pairs of adjacent fractions were pooled to create approximately 20 final fractions from each cytosol and mitochondria preparation and up to 200 μl of each were dried by vacuum concentration. Samples were then digested by trypsin (Sigma, T6567) as described above, except 50 μl of urea digestion buffer and an enzyme: substrate ratio of 1: 10 was used. Following vacuum centrifugation to dryness, samples were dissolved in up to 100 μl of 5% formic acid, centrifuged at 16,000 x g for 10 min to remove particulates, and approximately 2 μg of each digest injected onto an Acclaim PepMap 100 μm x 2 cm NanoViper C18,5 μm peptide trap on a switching valve. After 10 min of loading at 5μl/min, the trap column was switched on-line to a PepMap RSLC C18,2 μm, 75 μm x 25 cm EasySpray column (Thermo Scientific). Peptides were then separated using a 7. 5–30% ACN gradient over 60 min in mobile phase containing 0. 1% formic acid at a 300 nl/min flow rate and ionized using an EasySpray NanoSource (Thermo Scientific). Survey scans were performed by the Orbitrap mass analyzer at a resolution of 120,000, and data-dependent MS2 scans acquired by the linear ion trap in rapid mode. Data-dependant scanning used a 30 sec exclusion time, repeat count of 1, exclusion of +1 charge ions, top speed mode, and a 3 sec dwell time between survey scans. MaxQuant (version 1. 5. 2. 8) was used for peptide and protein identification from Orbitrap Fusion data. MaxQuant was configured as follows: FASTA database was a T. brucei protein database containing 11,567 entries (www. tritrypdb. org, release 24. 0), default decoy sequence generation, tryptic enzymatic cleavage, default parent ion and fragment ion mass tolerances, a maximum of 2 missed cleavages, static C+57. 02 modifications, variable M+15. 995 modifications, and target peptide FDR of 1% (cytosolic) or 2% (mitochondrial). A high protein FDR target was chosen because MaxQuant does not compute protein FDR per sample in multi-sample experiment designs; per sample protein FDRs are much less than overall experiment-wide protein FDR when there are many samples. Protein quantification used the protein total intensity values in the MaxQuant protein summary files. We also invoked a two peptide per protein per sample (one LC run) criteria to reduce protein identification noise. The RNA interference (RNAi) vector was prepared by amplifying 433 bp fragment from the open reading frame of Tb927. 10. 7910 gene. The amplified fragment was introduced into p2T7-177 [29]. Linearized plasmids were transfected into procyclic form 29–13 cells by electroporation, and transfected cells were selected by adding 2. 5 μg/ml pheleomycin. Growth effect in RNAi cell lines was monitored for up to 10 days in the presence and absence of 2. 5 μg/ml tetracyclin. RNA was collected from uninduced and Tet-induced RNAi Tb927. 10. 7910 cell line at day-3 post induction, selected based on the observed growth defect. RNA was isolated from ~108 cells using TRIzol (Invitrogen). Ten micrograms of total RNA were treated with DNA-free DNase Kit (Ambion) to remove any residual DNA. Purified RNA was reverse transcribed in 25 μl RT-PCR reactions including TagMan reverse transcriptase and random hexamer primers using TagMan Reverse Transcription Kit (Applied Biosciences). Different RT-PCR reactions were performed with primers specific to pre-edited, edited, and never-edited mitochondrial transcripts, as described previously [30] and primers that flank the junctions of the adjacent genes, 9S/ND8, Cyb/A6, and RPS12/ND5 [31] using a Corbett rotor gene 3000. Data were normalized against the 18S rRNA as an internal reference, and were confirmed for several genes considering the β-tubulin as the internal reference. Each mRNA target was analyzed in two biological replicates and three technical replicates. Relative changes in RNA abundance were calculated using ΔΔCt, as described before [32]. The sequence information for primers used in this study is provided in S1 Table. To generate C-myc-tagged cell lines, the open reading frames of Tb927. 1. 1730, Tb927. 10. 1730, and Tb927. 10. 7910 were amplified and cloned into BamHI and HpaI sites (for Tb927. 1. 1730) or BamHI and HindIII sites (for Tb927. 10. 1730 and Tb927. 10. 7910) of pHD-1700 [33]. All constructs were digested with NotI before transfection into 13–13 procyclic form cells. C-myc tagged cell lines were induced with 0. 5 μg/ml tetracycline for 48 hr and then used to prepare slides for immunofluorescence analysis. Briefly, mid-log phase cells were fixed with 4% paraformaldehyde in PBS and placed on poly-L-lysine-treated round coverslips (Fisher). Cells were permeabilized with 0. 2% Triton-X-100 in PBS and blocked with 3% BSA in PBS. Anti-myc antibody (1: 500) was used to visualize myc-tagged cells. Mitotracker (Invitrogen) was used to stain the organelles and DAPI was used for DNA staining. Slides were analyzed using a Nikon up-right microscope. Hypotonic purification of mitochondria from ~5 x 107 uninduced and tet-induced PF cmyc-Tb927. 1. 1730, cmyc-Tb927. 10. 1730, and cmyc-Tb927. 10. 7910 cells was carried out as described above. The lysis buffer was prepared either with 40 U of RNaseOUT (Invitrogen) or 200μg/ml of RNase A from bovine pancreas (Sigma). Anti c-myc agarose conjugated beads (Sigma) were washed five times with 1ml ice-cold PBS at 4°C and subsequently washed once with 1ml of ice-cold immunoprecipitation wash buffer (Tris-HCL, pH 8. 0,10mM, NaCl 100mM, NP-40 0. 1%, 1X complete EDTA-free protease inhibitor (Roch) and 1% PBS). After the last wash, 50 μl of beads were re-suspended for each reaction in 1 ml of ice-cold wash buffer and incubated for 1 hour at 4°C on a tube rotator. After adding the mitochondrial lysate to the beads, mixture was rotated for 2 hours at 4°C followed by centrifugation at 500rpm for 1 min at 4°C. After removing the supernatant (unbound proteins), the beads were washed four times with 1ml of immunoprecipitation wash buffer and then resuspended in SDS-PAGE loading dye. Aliquots of lysate, bound and unbound fractions were loaded on 10% SDS-PAGE gel. Proteins were transferred onto nitrocellulose membrane and probed with polyclonal antibodies against MRB 8170 and TbRGG2 (a generous gift from Laurie Read, state University of New York at Buffalo, USA). Although the stoichiometric ratio among proteins that are involved in one complex is expected to be linear, this relationship can be more complicated for those participating in several complexes (As an example, see the fractionation patterns for proteins related to the RNA editing process presented in the validation section). For stringent analysis of data, we used a previously developed mutual information-based approach, termed context likelihood of relatedness (CLR) [34], to infer pairwise interactions among the proteins based on the observed fractionation patterns. In contrast to association measures such as Pearson correlation coefficient, mutual information is a measure of association that does not have strong prior assumption about the association type [35]. After calculating fractionation pattern similarity scores for each possible interaction using mutual information, CLR determines significant pairwise interactions by comparing the similarity score of each protein pair to a background joint distribution, which was obtained by the calculation of similarities for all possible interactions that each of the interacting proteins can be involved in. To construct the co-fractionation networks, we set the false discovery rate cut-off threshold at 0. 05. False discovery rates were calculated using the provided functions in the CLR package [34]. However, since mutual information does not discriminate the type of association (positive or negative), we considered only those interactions as valid that their fractionation patterns were non-negatively correlated with each other as judged by Pearson correlation coefficient. Comparison of results for the mutual information-based version of CLR with the correlation-based version of this algorithm as well as conventional Pearson correlation coefficient-based method indicated that the former is more consistent with the previous knowledge on T. brucei protein complexes. To assess whether the observed co-fractionation patterns for low abundant proteins is significant or unreliable, we generated 100 noisy datasets from each fractionation experiments. Expecting a uniform distribution of noises in the dataset, added noise to each cell in a dataset were modeled by Poisson distribution with the lambda equal to the value of the cell plus the noise term of [the lowest identified ion intensity in the dataset]/[No. of fractions in the dataset]. Those protein pairs that were significantly co-fractionated (FDR ≤0. 05) in the original dataset, but lost their significant co-fractionation (p-value >0. 05) in at least 10% of noisy datasets were discarded from the corresponding co-fractionation network. Modulation score examines the density of connections among a pre-specified group of nodes (e. g. , proteins) in a given-network with the interaction density that is expected to occur by chance [36]. The density of interactions for a group of nodes is defined as the number of within-group interactions divided by the total number of interactions of the members of the group. This ratio is then compared to the density distribution generated over random groups with the same number of nodes in the network. In the original implementation of the modulation score, the density distribution was estimated by normal distribution assumption [36]; however, the distribution may not necessarily be normal. Therefore, we developed a modified version of the algorithm that uses a kernel smoothing function to estimate the distribution. To find high-confidence interactions, we considered following available orthogonal resources that contain protein-protein interaction data: To estimate the precision (percentage of interactions that are true) and recall (the fraction of overall true interactions that are present in the networks) of the constructed networks, benefiting from our extensive literature search, we created a gold standard set of >200 proteins in 20 distinct protein complexes (S3 Table), allowing us to evaluate our network against current literature. It should be noted that some of these complexes are partially identified and there might be some subunits that are still uncharacterized. We also ignored putative subunits that were suggested solely based on a single experiment like pull down without further experimental characterizations and/or verifications. The interactions between proteins in the same complex were defined as the positive set. Negative set were defined as interactions between proteins of different complexes. For more stringent analysis, we added to the negative set the other interactions that the subunits of the 20 complexes had with the rest of proteins in the network, unless they had external experimental evidences (e. g. , pull down). Clearly, the defined negative set can contain some of genuine interactions that were not identified previously, leading to an underestimation of the calculated precision for the networks. Precision was defined as the number of interactions in a network that are in the positive set, divided by the total number of interactions in the network that belong to either positive or negative sets. For each network, two types of recalls were defined: 1) Recall for proteins: defined as the number of distinct subunits of the twenty complexes that are present in the network divided by the total number of known subunits identified in our mass spectrometry experiments; 2) Recall for the interactions: defined as the total number of interactions that the identified subunits of the same complex have in the network divided by the maximum number of interactions that they can have theoretically. Since in a network, the subunits of the same complex are often much sparser than the theoretically assumed fully-connected model, the calculated recalls for the interactions are usually small [10,13]. All networks were visualized using Cytoscape, a network visualization tool for Genome Space workflows [39]. NetworkAnalyzer, a plugin of Cytoscape, was used for the topological analysis of networks, including: node-degree distribution, shortest path-length distribution, and topological coefficient distribution [40]. To examine the gene ontology-biological process similarity of interacting genes in the constructed networks, we used Resnik’s approach to quantify semantic similarity between gene ontology terms [41]. Semantic similarities between GO categories of biological processes and cellular compartments were calculated for each interacting protein pair present in a network, if both proteins were annotated in the uniprot database [42]. When calculating semantic similarities, we ignored those terms with evidence codes of NR (Not recorded), ND (No biological data available), and IEA (Inferred from Electronic Annotation). The GO-BP sematic similarities were calculated by GOssTO tool using default parameters [43]. MATLAB R2014b software (The MathWorks Inc. , Natick, MA) was used for the CLR score false discovery rate estimation and kernel-based p-value estimation. C# programing language was used for the calculation of the KEGG pathways modulation scores as well as construction of random networks. Other statistical analysis was performed using the R programing language. The mass spectrometry data have been deposited to the ProteomeXchange Consortium (http: //proteomecentral. proteomexchange. org) via the PRIDE partner repository [44] with the dataset identifier PXD002640. Full instrument settings, search parameters, pipeline processing details, and dataset statistics can be found in the PRIDE submission.
  Biochemical fractionation techniques have been widely applied to trypanosomatid organisms to test the possibility of physical associations among a set of pre-specified proteins. Fractionation approaches allow dissection of protein complexes based on different biochemical properties. In the glycerol gradient (GG) fractionation approach, protein complexes become separated according to their shape/density. Alternatively, complexes were fractionated on the basis of their overall charge in ion exchange high performance liquid chromatography (IEX) experiments. As summarized in Fig 1A, by coupling GG and IEX deep fractionation techniques with semi-quantitative, ultra-sensitive, mass spectrometry, we generalized the approach to chart a proteome-scale T. brucei interaction network. We were able to observe the fractionation pattern of 3354 protein groups (paralogous proteins with nearly identical sequences were grouped together) across total of 133 separate fractions from whole cell-GG (48 fractions), mitochondrial-GG (46 fractions), cytosolic-IEX (19 fractions), and mitochondrial-IEX (20 fractions) experiments on T. brucei procyclic form cells (Fig 2A). Due to complementary design of experiments, we were able to observe the fractionation patterns of 1398 (42% of total) proteins in both GG and IEX fractionation approaches, providing a global picture of protein complexes present in T. brucei procyclic cells (Fig 2B). Comparison of our results from mitochondrial experiments with a previously reported repository of mitochondrial proteins [25] revealed that 79% of the detected proteins in our experiments were independently supported to be present in the mitochondria. Comparing the mitochondrial signal probability of the proteins identified in our mitochondrial samples with the above mentioned list and also a list of transcripts that are expressed in procyclic life stage of T. brucei [45] suggested that the identified proteins were significantly enriched (comparing to proteins identified in [25]; p-value < 1. 4E-08, Wilcoxon-Mann-Whitney rank sum test) for the mitochondrial proteins (Fig 2C). The same analysis on the proteins identified in the cytosolic experiment demonstrated a significant depletion (comparing to total procyclic transcripts; p-value <1. 4E-07, Wilcoxon-Mann-Whitney rank sum test) of mitochondrial proteins in the sample. These results demonstrate the accuracy of the employed compartment enrichment procedures. In total, proteins identified in our fractionation experiments cover 43% of all expressed protein coding genes in T. brucei procyclic stage [45] and 77–127% of number of proteins reported in previous proteome-wide SILAC studies [47–49]. To test for bias in the data, we examined the distribution of detected proteins in terms of cellular component and number of putative transmembrane domains. As expected, soluble proteins were preferentially detected by the employed fractionation approaches, while the membrane associated proteins were under-represented (Fig 2D and S1 Fig). The detected proteins did not show bias in terms of protein length (S1 Fig). Moreover, observing the transcriptome responses of the identified proteins in different life stages [50] and over the differentiation process [51,52] indicated that they show over-expression trend in the procyclic life stage, but biased towards more abundant transcripts in the cell (S2 Fig). Because paralogous proteins with nearly identical sequences are expected to play similar functions in the cell, we randomly selected one representative protein from each protein group to construct a protein interaction map based on the observed fractionation patterns. The relative abundance of proteins in each fraction were measured based on label-free MS2 ion intensity-based approaches (Detailed in the method section). To decipher physical associations from the data, we reasoned physically interacting proteins would show similar patterns across the biochemical fractions. The employed bioinformatics approach for the inference of co-fractionation networks based on the observed patterns is summarized in Fig 1B. As described below, for stringent analysis of data, physically interacting protein pairs were predicted based on five criteria: 1) showing significant similarity in their fractionation patterns; 2) being robust to the addition of noise in their observed co-fractionation similarity; 3) showing shared-peak in at least one fraction; 4) being stably interacting with protein complexes as judge by GG-derived patterns; 5) reproducibility of co-fractionation when the same fractionation technique was used. Initially, we constructed four separate co-fractionation networks by application of the context likelihood relatedness (CLR) algorithm [34] on each of four datasets, separately. As described in the methods section, the CLR algorithm is an information theoretic-based approach that can predict the association of two random variables (proteins) based on their observed values (fractionation patterns). It predicts two proteins as interacting only if the similarity of their patterns is significantly higher than what expected based on the background (estimated Benjamini-corrected FDR of 5%). Low abundance proteins might not have reliable patterns in the dataset and therefore show random similarities to other proteins in their fractionation patterns. To confirm that observed similarity is not due to presence of noise in the dataset, co-fractionation similarities were recalculated after the generation of Poisson noise models from each dataset. In each network, we discarded those interactions that lost their significant similarity due to the addition of noise (Detailed in the method section). To further reduce the possibility of chance co-fractionation, we only kept those interactions for which both interacting proteins had a peak in at least one shared fraction. A protein was deemed to have a peak in a fraction if its ion intensity for that fraction was at least 80% of its second-highest intensity. We observed that filtering the interactions by this more stringent criterion led to an increased reproducibility rate of results and, most likely, the elimination of putative false positive interactions from the networks (Fig 3A). In GG fractionation experiments, early fractions are highly enriched for monomeric proteins in the cell (e. g. , many enzymes) or those that are not stably involved in the complexes. Since interacting partners for these groups of proteins cannot be reliably identified by our approach, proteins that had a peak only in the top two fractions of GG experiments were discarded from the networks. Cytosolic-IEX and mitochondrial-IEX networks were significantly depleted for proteins of one another (p-value < 6E-35), most likely because of the enrichment procedures used for these two experiments. To check the technical reproducibility of the results, we focused on interactions occurring among proteins that were detected in both GG experiments. As illustrated in Fig 3B, we found that more than half of co-sedimented protein pairs (940 reproducible interactions) in one GG experiment were also co-sedimented in the other (FDR ≤0. 05 in one experiment and p-value ≤0. 05 in the other). Moreover, comparing the number of common interactions between whole cell-GG and mitochondrial-GG networks with random networks with the same structural characteristics revealed that these two networks are significantly enriched for reproducible interactions (S3 Fig). These results indicated that significantly co-sedimented protein pairs in the whole cell-GG experiment are typically co-sedimented in the mitochondrial-GG experiment and vice versa. To increase the accuracy, we removed the interactions among those proteins that were not consistently co-fractionated with each other in both experiments (i. e. , despite the detection of proteins in both experiments, they co-fractionated in only one experiment). Overall, merging GG networks led to a network composed of 12,196 interactions among 1,417 proteins. The IEX experiments led to the identification of 1,708 interactions among 1,261 proteins. The smaller size of IEX network compared to those of GG is because former experiments had smaller number of fractions and, therefore, were less informative on protein complexes compared to GG experiments. We next assessed the agreements between IEX-derived and GG-derived networks. As shown in Fig 3C, comparison of cytosolic-IEX network with whole cell-GG network indicated that interactions have peak in the reproducible region (median correlation of 0. 47 and 0. 83 for cytosolic-IEX and whole cell-GG, respectively), but skewed toward the region with negative IEX correlations. Importantly, none of significantly co-eluted protein pairs in cytosolic-IEX experiment had correlations less than 0. 47 in the whole cell-GG dataset. Comparison of mitochondrial-IEX network with mitochondrial-GG network led to a similar result (Fig 3D); i. e. , interactions had peak in the reproducible region (median correlation of 0. 53 and 0. 73 in mitochondrial-IEX and mitochondrial-GG datasets, respectively). Moreover, comparison with random graphs demonstrated that the IEX-derived and GG-derived networks were significantly enriched for the reproducible interactions (S4 Fig). Due to separation of protein complexes in IEX and GG experiments based on different biochemical properties, the confounding complexes in one approach have lower chance to co-fractionate in the other approach as well. However, the observed discrepancies in co-fractionation patterns are not only because of confounding elements, but rather reflect the differences in the nature of the two fractionation experiments as well. As an illustration, we observed that members of translation initiation complex co-sedimented in the whole cell-GG network with the median correlation of 0. 85. However, their median correlation in the cytosolic-IEX experiments reduced to 0. 53. Likewise, comparison of the co-fractionation patterns among a set of proteins enriched for RNA-dependent interactions demonstrated their strong co-sedimentation in the GG experiments, but separation in the IEX experiment (S5 Fig). This discrepancy can be due to the use of salt gradient in the latter experiment and disruption of the less stable ionic interaction complexes. However, these effects were minimal on complexes that are known to form more stable complexes such as ribosome, proteasome, F0F1 ATPase, and core editosome (discussed more on the validation part). These results suggested that our GG fractionation experiments better preserved the less stable interactions (i. e. ionic-based interactions), but our IEX experiments favored more stable interactions. For the follow-up analysis, we merged the GG and IEX networks together. This global physical map, named TbCF net (T. brucei co-fractionation network), was constructed by distinguishing two types of interactions (S2 Table): 1) Those protein pairs that were reproducibly co-fractionated in both GG and IEX networks (FDR ≤0. 05 in one experiment and p-value ≤0. 05 in the other); and 2) Those protein pairs that were co-fractionated only in one experiment. The TbCF net connects 2,151 proteins with 13,865 interactions. The orthogonal reproducible part of TbCF network, termed TbCFOR net, was composed of 2,601 (19% of total) interactions among 828 (38% of total) proteins. To test the validity of TbCF net, we first examined the topological properties of the constructed network (Fig 4A). In agreement with protein interaction networks of model organisms, TbCF net has scales-free architecture [53]; i. e. , while most proteins interact with a small number of proteins, some of them (known as hubs) are highly linked to the other proteins (Fig 4B). Moreover, as shown in Fig 4C, each protein typically can be reached from every other protein by a small traverse in the network (i. e. , small path length), as expected from a network with a small-world property [54]. Additionally, we observed that expect for ribosomal proteins which make a large, densely intra-connected module in the TbCF net, the topological coefficient also decreases with the number of neighbors (Fig 4D), reflecting that the number of common neighbors for hub proteins compared with the other proteins in the network is relatively low. This feature indicates that highly interacting proteins are not sporadically connected to each other. We next assessed different biological features that can be expected from a protein interaction network. To this end, we first examined whether interacting proteins in the network tend to be involved in the same biological process. As shown in Fig 5A and S6 Fig, we observed that higher similarity scores of fractionation patterns consistently led to an increased probability that two proteins are involved in the same biological process, as judged by gene ontology-biological process (GO-BP) analysis. Importantly, this analysis suggested the higher precision of TbCFOR net compared to the TbCF net. KEGG pathway [55] analysis also led to similar results, i. e. , the chance of two proteins participating in the same KEGG pathway increases by having a higher similarity score (S7 Fig). Although proteins in the same KEGG pathway do not necessarily interact with each other as they can have functional rather than physical associations, protein complexes are involved in some KEGG pathways such as ribosome and proteasome. Importantly, we found that KEGG pathway interactions that were also supported by the TbCF net had significantly higher (p-value <3. 5E-33, Wilcoxon-Mann-Whitney rank sum test) co-localization scores compared with the other KEGG interactions (Fig 5B), suggesting the physical nature of associations for the captured interactions. In a more stringent analysis, we next investigated whether or not proteins involved in a same KEGG pathway tend to form a significant module in the TbCF net. To this end, we used a score, termed modulation score [36], that assess the density of connections among a set of pre-specified proteins; e. g. , proteins with a shared KEGG attribute (See methods for details). This analysis revealed that members of 15 KEGG pathways form densely connected modules (p-value<0. 05) in the TbCF net (Fig 5C). Interestingly, we observed a significant modulation score for many pathways that are well-known to involve protein complexes (ribosomes, proteasome, RNA transport, oxidative phosphorylation, etc.). The same analysis on the TbCFOR net indicated that this network is not biased towards a specific process and most KEGG pathways remain or even become significantly connected in the TbCFOR net. However, as discussed earlier, TbCFOR net were depleted for those pathways that involved less stable and salt sensitive interactions (Fig 5C). We also analyzed the TbCF net in terms of gene expression responses. In model organisms, it has been shown that transcriptional co-regulation of proteins plays a major role in the efficient control of the cell homeostasis in different environments [17]. It has also been shown in T. brucei that functionally-related proteins tend to be co-expressed with each other [56]. Indeed, gene expression analysis indicated that interacting proteins in the TbCF net show significantly high (p-value < 2. 2E-16, Wilcoxon-Mann-Whitney rank sum test) co-expression trends in different life stages and also during the differentiation process (Fig 5D and S8 Fig). Consistent with GO-BP semantic similarity analysis, co-expression analysis also indicated a significantly higher (p-value < 2E-07, Wilcoxon-Mann-Whitney rank sum test) co-expression trend for interactions of TbCFOR net compared to those of TbCF net, suggesting a higher precision of the former network. To estimate precision of TbCF net, we focused on twenty complexes that have been identified experimentally in T. brucei, composed of more than 200 proteins (S3 Table). This analysis demonstrated a precision of 34% in TbCF net with a recall of 80% on proteins and 17% on interactions, which is comparable to previous high-throughput studies of protein interactions based on biochemical fractionation [10,13] (see methods for details). However, the estimated precision can be an underestimate as in many cases, subunits of the complexes had been partially identified because of application of stringent conditions. Importantly, consistent with the expected higher precision of TbCFOR net as judged by GO-BP and co-expression analysis, the same analysis on TbCFOR net with literature interactions suggested a precision of 59% (~2-fold increase in precision). As mentioned in the introduction, the integration of different data sources is highly recommended to reduce false positive results from an interactome. However, due to the lack of a global, unbiased protein interaction map for T. brucei, we were not able to utilize the state-of-art machine learning approaches to systematically integrate additional data sources. To experimentally examine the extent to which the precision of TbCF network can be further improved by other orthogonal information, we performed TAP-TAG TEV-elution of RNA editing ligase 1 (REL1) protein. We selected this protein because it is a subunit of a well-studied T. brucei complex, the core editosome. TEV-elution is a low stringency condition in which transient interactors along with many contaminants can co-purify with the tagged protein. TEV-elution of REL1 protein led to the co-purification of 83 proteins as putative interacting partners (S4 Table). We extracted the sub-network from TbCF net that was restricted to the co-purified proteins. Consistent with that, REL1 was directly connected to its known interacting partners, while well separated from other contaminant proteins in the network (Fig 6A). It should be noted that since only a small subset (i. e. , five proteins) of previously known interacting partners of REL1 protein were co-purified in the pull down experiment, and the constructed sub-network was restricted to the identified proteins in that experiment, the other previously known interacting proteins that were also connected to REL1 in TbCF net are not represented in the sub-network, indicating the importance of queried proteins on obtaining a comprehensive sub-network. To further assess the TbCF net, we considered six separate pull down experiments of aminoacyl-tRNA synthetase (aaRS) proteins in different T. brucei life stages [57]. In these pull down experiments, 262 proteins were co-purified as candidate proteins involved in T. brucei tRNA-synthesis with a minimum overall mass spectral count of two. Following the same procedure as that used for REL1, the sub-network of these proteins suggested the existence of several distinct complexes among the proteins co-purified with aaRS proteins (Fig 6B). Interestingly, one of these complexes was highly enriched for proteins involved in tRNA-synthesis including the recently identified members of MARS complex [57] and two additional hypothetical proteins. Therefore, the TbCF net was able to successfully distinguish contaminant proteins from direct interactors in both cases. Integration of various TAP-TAG experiments with the TbCF net also suggested that the chance of retrieving a high false positive rate with a targeted search, i. e. , a list of putatively interacting proteins, is very low. In our next attempt, we considered the functional protein interaction network of T. brucei deposited in the STRING database [58]. In the STRING database, most of available functional association data for T. brucei is inferred based on the indirect approaches (such as text mining and co-expression) rather than direct experimental evidence (S9 Fig). We extracted the STRING network using a medium confidence level threshold (the default threshold set by STRING) for the proteins present in the TbCF net. The retained STRING network was densely connected with 19,119 interactions among 1,402 proteins (S10 Fig). Next, we derived a secondary network composed of the interactions that were supported by both STRING and TbCF net, termed TbCFSTRING net. The TbCFSTRING net was composed of 2,413 interactions among 449 proteins (that included 13% of interactions and 32% of proteins in the primary STRING network, S10 Fig). To assess the validity of the network, the TbCFSTRING net was clustered using the clusterOne algorithm [59]. As shown in Fig 6C (and presented in more details in S1 File), many previously known protein complexes were recovered by this approach, indicating the high accuracy of the TbCFSTRING net. This result indicated that integration of our co-fractionation network with other independent sources (like AP and STRING in this case) leads to the elimination of false positive interactions from both sources. Based on the findings in the previous step, we partitioned the interactions in the TbCF net into two parts; those with high confidence, TbCFHC, and those with no external evidence, TbCFNE. As schematically shown in Fig 1A, the high confidence group was composed of interactions that were supported by at least one of the orthogonal resources including KEGG pathways, the STRING database, interlog-mapping, extensive literature searches, and orthogonal reproducibility (see methods for details). As discussed earlier, interaction data in each of these orthogonal resources (e. g. KEGG, STRING, AP, etc.) suffers from false positives and also, for some sources, does not imply physical interaction among protein pairs. However, we expected to observe elimination of the false positive interactions by integration of these data with our fractionation-derived network. TbCFHC network was composed of 4,726 interactions among 866 proteins (S2 Table), encompassing 34% of the total interactions in TbCF net (Table 1). Analysis of TbCFHC net indicated its improvement over TbCF net in terms of GO-BP (S12 Fig; p-value <2E-08, Wilcoxon-Mann-Whitney rank sum test), GO-CC (S12 Fig; p-value <6E-07, Wilcoxon-Mann-Whitney rank sum test), and co-expression (S12 Fig; p-value <5E-117, Wilcoxon-Mann-Whitney rank sum test). As judged by the analysis of the same twenty protein complexes discussed above, the estimated precision for TbCFHC net was 80% and 68% without and with the exclusion of literature-derived data from the network, respectively. The estimated precisions for both TbCFOR and TbCFHC networks suggest that integration of TbCF net with other orthogonal resources leads to the overall false discovery rate of less than 40%. Clustering TbCFHC net using ClusterOne algorithm [59] led to the prediction of 128 protein complexes among 716 proteins (S5 Table). The predicted complexes varied in size between three and 70 with the median size of four proteins per complex that is similar to reports from other organisms [10]. Many of the predicted complexes were significantly enriched for subunits of previously known T. brucei complexes. For example, as illustrated in S13 Fig, we were able to successfully identify complexes related to cytoplasmic and mitochondrial ribosomes, proteasome, T-Complex, intraflagellar transport, translation initiation, and mitochondrial RNA editing. To further assess the quality of predicted complexes, we associated them with the available large-scale RNAi screening data for T. brucei [60]. Because each protein complex act as a functional unit in the cell, it is expected that the essential proteins to be over- or under-represented in a protein complex, depending on the function of the complex [61]. Indeed, enrichment analysis of predicted complexes indicated over representation of complexes with enriched or depleted fraction of essential proteins, while complexes with random distribution of essential proteins were under represented (Fig 7A). The predicted essential complexes recapitulated the previous findings for T. brucei. As an illustration, the results indicated that the cytoplasmic ribosomal complex (Complex 1), proteasome (Complex 3), and T-complex (Complexes 10 and 11) are essential in all life stages of the parasite (Fig 7B). Conversely, the complex 7 which is highly enriched for intraflagellar proteins (Fig 7B), was not essential in the procyclic life stage of T. brucei [62]. Additionally, complexes related to the mitochondrial ribosome (Complex 2 and 8) were depleted for the essential proteins in the bloodstream life stages (Fig 7B). Clustering of TbCFHC net also led to the prediction of complex membership for 188 protein groups currently annotated as hypothetical, 350 protein groups with the annotated name as putative, and 635 protein groups lacked experimental GO-BP annotation (evidence codes of EXP, IDA, IPI, IMP, IGI, and IEP according to TriTrypDB v5). To experimentally verify the quality of predictions, we focused on the overlapping complexes that were highly enriched for proteins involved in the Kinetoplastid RNA editing process. Mitochondrial gene regulation is a highly interesting, yet not fully understood, process in T. brucei. In this process, mitochondrial genes are transcribed as polycistronic units. After cleavage, mitochondrial transcripts (mtRNAs) become stabilized by the addition of short tails at their 3' ends [63]. Intriguingly, most of the produced mtRNAs originally do not possess correct open reading frames (ORFs) and require to be edited before the translation. In the editing process, small RNAs, known as guide RNAs (gRNAs), dictate the insertion, and less frequently deletion, of a defined number of uridine nucleotides at pre-specified positions in mtRNAs [64]. The edited mtRNAs become marked for the translation by the addition of long A/U tails [65]. This highly complicated and intertwined process provides the parasite multiple post-transcriptional regulatory layers over mtRNAs. Mitochondrial post-transcriptional gene regulation is essential for the survival of the parasite in both bloodstream and procyclic life stages, although the precise role of this process is not well understood in the bloodstream form [66,67]. Experimental evidence confirms differential regulation of some mtRNAs in at least editing step during the parasite’s life cycle [68,69]. For example, while the Cytochrome b (Cyb; a subunit of complex III) and cytochrome oxidase subunit II (COII) transcripts are preferentially edited in procyclic form, some subunits of complex I such as NADH-ubiquinone oxidoreductase subunit 8 (ND8) and NADH-ubiquinone oxidoreductase subunit 7 (ND7) are mostly edited in the bloodstream stage. This developmental regulation in RNA editing is coordinated with the activities of the trypanosome mitochondria during its life cycle to allow the adaptation and survival of the organism in changing environmental conditions [70]. Consistent with the developmental regulation of the RNA editing process, comparative sedimentation analysis of the RNA editing machinery demonstrated that the complexes associated with the RNA editing machinery of bloodstream and procyclic forms are not identical [71]. The TbCFHC net suggested the physical association of 49 protein groups (50 proteins) with the RNA editing machinery of T. brucei in the procyclic life stage (S6 Table). Many of these proteins are well-known to be involved in the RNA editing process. For example, among the 50 predicted proteins, 17 are known to be involved in the core editosome complex [72], while 21 function in the MRB1 complex [73]. For detailed analysis of their interactions, we distinguished predicted interactions originated from IEX experiments to those of GG experiment. This analysis indicated that members of core editosome and some of the accessory elements were reproducibly co-fractionated in both approaches (S15 Fig). However, about 70% of significantly co-sedimenting protein pairs in the GG experiments were not significantly co-eluted in the IEX experiments (S15 Fig). Visual inspection of the fractionation patterns of these 50 proteins confirmed that they mainly co-sedimented together in mitochondrial-GG experiment, but dissociated into different protein clusters in the mitochondrial-IEX experiment (Fig 7C). It should be pointed out that many proteins that are functionally associated with the RNA editing machinery, mediate their functions through binding to the RNA and GG and IEX mitochondrial fractionation experiments were not RNase treated. However, technical differences between the IEX and GG experiments provided a high resolution picture of the distinct complexes involved in the mitochondrial post-transcriptional regulation of T. brucei. For example, while members of core editosome were reproducibly co-fractionated in both approaches, comparison of IEX and GG data suggested existence of at least three distinct groups of proteins among the accessory elements. These three groups were co-sedimented with each other and with core editosome in the mitochondrial-GG experiment. However, they failed to co-fractionate in the mitochondrial-IEX experiment possibly due to the increased salt concentration, suggesting RNA-dependent or less stable nature of interactions between these groups. Consistent with previous reports [73], observed fractionation patterns supported a direct interaction between Tb927. 8. 8170 and Tb927. 11. 16860, while suggested RNA-dependent interaction of Tb927. 2. 6070 with MRB core proteins (i. e. , Tb927. 5. 3010, Tb927. 10. 11870, and Tb927. 10. 10130). However, comparison with previous findings on the interactome of proteins related to the RNA editing machinery suggested that the lack of co-elution in the IEX data could be also because of less stable direct protein-protein interactions. For example, although MRB8170 (Tb927. 8. 8170) and TbRGG2 (Tb927. 10. 10830) did not co-eluted in the IEX experiment, they were reported to directly interact in Y2H assays [73,74], or AP-based studies [75]. Importantly, a previous study has demonstrated the TEV co-elution of the two proteins in a RNA-enhanced manner [76]. Hence, the lack of co-elution in the IEX fractionation experiment for MRB8170 and TbRGG2 proteins is likely due to less stable interaction of the two proteins. Moreover, fractionation data indicated that members of KPAP1 polyadenylation complex reproducibly fractionated differently from those of core editosome, and suggested their association with mitochondrial ribosome, which is consistent with their functional role that couple the mitochondrial editing with the translation process [65]. Also, the fractionation patterns successfully captured interactions of KPAP1 protein (Tb927. 11. 7960) with both editing and ribosomal complexes (Fig 7C), recapitulating the results obtained by mass spectrometric and immunochemical experiments [65]. The TbCFHC net also suggested that six new proteins (Tb927. 1. 3010, Tb927. 10. 7910, Tb927. 1. 1730, Tb927. 10. 5830, Tb927. 6. 1200, and Tb927. 10. 1730) play a role in the mitochondrial post-transcriptional gene regulation. Fractionation patterns of GG experiment for these candidates indicated their co-sedimentation with core editosome and accessory elements, but separation from subunits of KPAP1 polyadenylation complex. Importantly, IEX fractionation patterns suggested that two of these candidates (Tb927. 6. 1200 and Tb927. 10. 1730) strikingly co-fractionated with members of the core-editosome and some accessory elements, but the co-elution of the other four proteins were, in varying degrees, sensitive to the presence of salt with Tb927. 10. 7910 being the most salt sensitive one. Examining the localization pattern for three of our candidate proteins (Tb927. 10. 7910, Tb927. 10. 1730, and Tb927. 1. 1730) clearly confirmed their exclusive presence in the mitochondria of the parasite (S16 Fig). We also performed immunoprecipitation experiments on tagged versions of these three proteins. These experiments verified their interactions with TbRGG2 and MRB8170 proteins. Consistent with the observed IEX fractionation patterns, we found that the interaction of all three candidates with MRB8170 was abolished following the RNase treatment (Fig 8A). Likewise, we found that only Tb927. 10. 1730 remains still bound to TbRGG2 after RNase treatment (Fig 8A). Consistent with this work, another study reported four out of our six candidates (including Tb927. 1. 3010, Tb927. 1. 1730, Tb927. 6. 1200, Tb927. 10. 1730) form a novel complex involved in the post-transcriptional regulation of mtRNAs, termed polyadenylation mediator complex, confirming the predictions on these four proteins [77]. To further assess the role of the candidate proteins in the RNA editing process, we performed tetracycline (Tet) -inducible RNA interference (RNAi) knockdown of Tb927. 10. 7910 in the procyclic form of T. brucei. RNAi induction for this gene led to growth-defect phenotype, reflecting its essential role in normal growth of the parasite in the procyclic life stage (Fig 8B). Follow up quantitative RT-PCR verified the knockdown of the candidate transcript compared with the control, uninduced cells (Fig 8C). We next quantified the relative changes in mitochondrial-encoded pre-edited, edited, and never-edited transcripts in the RNAi-knock down background. We also considered three precursor RNAs to examine whether this protein play a role in the precursor RNA processing. This experiment indicated that Tb927. 10. 7910 affect the RNA editing process as judged by the accumulation or reduction of pre-edited or edited transcripts for different target RNAs (Fig 8C). Interestingly, our results suggested that Tb927. 10. 7910 affect the editing process of the Cyb transcript (i. e. upregulation of pre-edited mRNA and down-regulation of the edited mRNA). The knock down of Tb927. 10. 7910 also led to the down regulation of RPS12 edited transcript, and also accumulation of MURF2 pre-edited as well as COIII and A6 edited mtRNAs, suggesting the multiple functionality of Tb927. 10. 7910 in the mitochondrial post transcriptional regulatory network of T. brucei. Consistently, a previous study has suggested the essentiality of Tb927. 10. 7910 protein in the bloodstream life stage of the parasite [60]. We have presented a systematic study of protein complexes in T. brucei using two complementary biochemical fractionation approaches. Our results led to the assignment of many previously uncharacterized proteins to complexes. The quality of predictions was verified by independent follow up experiments on newly characterized proteins associated with RNA editing machinery. Interestingly, we found that at least five out of six predictions of TbCFHC net are truly associated with the RNA editing machinery, and also that one of them preferentially affect the editing process of Cyb transcript, a developmentally regulated mitochondrial mRNA. Further experiments are required to clarify the roles of these proteins in the mitochondrial gene regulatory network of T. brucei. Despite the employment of stringent filtration criteria, our preliminary network (TbCF net) contains false positive interactions (estimated precision of 34%) that arose from limitations of the employed fractionation approaches. However, integration of results from two different fractionation approaches led to significant boosting of the precision (~2-fold), reflecting the importance of data integration for accurate predictions. Due to the lack of an unbiased, high-confidence, and large-enough protein complex map for T. brucei, we were not able to apply sophisticated machine-learning approaches to increase the precision of TbCF net by incorporating other data sources such as transcriptome. The KEGG interactions imply functional rather than physical associations. The interactions deposited in the STRING database were mostly inferred based on the indirect evidences that support functional associations. The literature-derived network was highly biased towards specific protein complexes and also contaminated with non-specific interactors. According to the primary literature-derived network, for example, there were ~360 co-purified protein groups with known members of RNA editing machinery. Of these, 242 protein groups were detected in our experiments. However, the TbCF net predicted strikingly similar co-fractionation pattern for only 49 of these protein groups, reflecting the high-contamination rate of literature-derived network. By stringent filtration of data, we focused on the high confidence sub-network of TbCF net (TbCFHC net), composed of interactions whose existences were supported by at least one other independent resource. The TbCFHC net refines and extracts new information from previous data and computational predictions on the interactome of T. brucei. However, this additional step can lead to the loss of information on proteins for which these types of evidences were not available. Therefore, we can expect that the number of high-confidence interactions will increase with the availability of more experimental data on T. brucei protein-protein interactions. Consistent with this, we have found that integration of results from other experiments (e. g. AP or immunoprecipitation) with TbCF net (which is a mixture of high- and low-confidence interactions) leads to the elimination of false positive results from both sources. Importantly, our analysis suggested that integration of TbCF net with other orthogonal resources leads to the overall false discovery rate of less than 40%. This finding has a great impact since most protein interactions in T. brucei are inferred by applying stringent conditions at the expense of an increased false negative rate, i. e. , losing transient interactions. Our results suggest that this criterion can be relaxed by considering TbCF net as an orthogonal validation resource. Mass spectrometry-based experiments are known to have limitations in the detection of low abundance proteins (including many regulator proteins), as they become masked in the sample by more abundant proteins. This issue was partially addressed in our approach by fractionating cell-extracts before mass spectrometry. Although comparable with previous SILAC experiments, our experiments cover 42% of proteins associated with expressed mRNAs in the procyclic stage [45]. Despite the fact that some of these transcripts can be translationally silent, the fraction of identified proteins is still potentially low. Our fractionation experiments demonstrated that enrichment for specific cellular compartments offers a viable solution to this issue. Hence, a greater depth could be obtained by performing similar experiments for other subcellular compartments of the cell such as the nucleus or mRNA enriched extracts. Moreover, our analysis clearly indicated an under-representation of membrane proteins in the mass spectrometry data that stems from the employed experimental procedure. Recovery of this class of proteins could likely be increased by use of detergents that can solubilize various membrane associated complexes and protein sub-domains. The results of this work are deposited in a prototype version of a database available at: www. trypsNetDB. org.
summary: Due to high evolutionary divergence of trypanosomatid pathogens from other eukaryotes, accurate prediction of functional roles for most of their proteins is not feasible based on homology-based approaches. Although protein co-complex maps provide a compelling tool for the functional annotation of proteins, as subunits of a complex are expected to be involved in similar biological processes, the current knowledge about these maps is still rudimentary. Here, we systematically examined the protein co-complex membership of more than one third of T. brucei proteome using two orthogonal fractionation approaches. A high-confidence network of co-complex relationships predicts the network context of 866 proteins, including many hypothetical and experimentally unannotated proteins. To our knowledge, this study presents the largest proteomics-based interaction map of trypanosomatid parasites to date, providing a useful resource for formulating new biological hypothesises and further experimental leads.
keywords: protein interactions
  protein interaction networks
  parasitic protozoans
  trypanosoma brucei gambiense
  protozoans
  network analysis
  mitochondria
  bioenergetics
  cellular structures and organelles
  extraction techniques
  research and analysis methods
  separation processes
  computer and information sciences
  protein extraction
  proteins
  fractionation
  proteomics
  biochemistry
  rna
  trypanosoma
  protein complexes
  cell biology
  nucleic acids
  rna editing
  biology and life sciences
  energy-producing organelles
  organisms
title: A Protein Complex Map of Trypanosoma brucei