Integrating transcriptomic and proteomic data for accurate assembly and annotation of genomes.

Complementing genome sequence with deep transcriptome and proteome data could enable more accurate assembly and annotation of newly sequenced genomes. Here, we provide a proof-of-concept of an integrated approach for analysis of the genome and proteome of Anopheles stephensi, which is one of the most important vectors of the malaria parasite. To achieve broad coverage of genes, we carried out transcriptome sequencing and deep proteome profiling of multiple anatomically distinct sites. Based on transcriptomic data alone, we identified and corrected 535 events of incomplete genome assembly involving 1196 scaffolds and 868 protein-coding gene models. This proteogenomic approach enabled us to add 365 genes that were missed during genome annotation and identify 917 gene correction events through discovery of 151 novel exons, 297 protein extensions, 231 exon extensions, 192 novel protein start sites, 19 novel translational frames, 28 events of joining of exons, and 76 events of joining of adjacent genes as a single gene. Incorporation of proteomic evidence allowed us to change the designation of more than 87 predicted "noncoding RNAs" to conventional mRNAs coded by protein-coding genes. Importantly, extension of the newly corrected genome assemblies and gene models to 15 other newly assembled Anopheline genomes led to the discovery of a large number of apparent discrepancies in assembly and annotation of these genomes. Our data provide a framework for how future genome sequencing efforts should incorporate transcriptomic and proteomic analysis in combination with simultaneous manual curation to achieve near complete assembly and accurate annotation of genomes.

A draft map of the human proteome.

The availability of human genome sequence has transformed biomedical research over the past decade. However, an equivalent map for the human proteome with direct measurements of proteins and peptides does not exist yet. Here we present a draft map of the human proteome using high-resolution Fourier-transform mass spectrometry. In-depth proteomic profiling of 30 histologically normal human samples, including 17 adult tissues, 7 fetal tissues and 6 purified primary haematopoietic cells, resulted in identification of proteins encoded by 17,294 genes accounting for approximately 84% of the total annotated protein-coding genes in humans. A unique and comprehensive strategy for proteogenomic analysis enabled us to discover a number of novel protein-coding regions, which includes translated pseudogenes, non-coding RNAs and upstream open reading frames. This large human proteome catalogue (available as an interactive web-based resource at http://www.humanproteomemap.org) will complement available human genome and transcriptome data to accelerate biomedical research in health and disease.

Widespread somatic L1 retrotransposition occurs early during gastrointestinal cancer evolution.

Somatic L1 retrotransposition events have been shown to occur in epithelial cancers. Here, we attempted to determine how early somatic L1 insertions occurred during the development of gastrointestinal (GI) cancers. Using L1-targeted resequencing (L1-seq), we studied different stages of four colorectal cancers arising from colonic polyps, seven pancreatic carcinomas, as well as seven gastric cancers. Surprisingly, we found somatic L1 insertions not only in all cancer types and metastases but also in colonic adenomas, well-known cancer precursors. Some insertions were also present in low quantities in normal GI tissues, occasionally caught in the act of being clonally fixed in the adjacent tumors. Insertions in adenomas and cancers numbered in the hundreds, and many were present in multiple tumor sections, implying clonal distribution. Our results demonstrate that extensive somatic insertional mutagenesis occurs very early during the development of GI tumors, probably before dysplastic growth.

Tissue matrix arrays for high-throughput screening and systems analysis of cell function.

Cell and protein arrays have demonstrated remarkable utility in the high-throughput evaluation of biological responses; however, they lack the complexity of native tissue and organs. Here we spotted tissue extracellular matrix (ECM) particles as two-dimensional (2D) arrays or incorporated them with cells to generate three-dimensional (3D) cell-matrix microtissue arrays. We then investigated the responses of human stem, cancer and immune cells to tissue ECM arrays originating from 11 different tissues. We validated the 2D and 3D arrays as representative of the in vivo microenvironment by means of quantitative analysis of tissue-specific cellular responses, including matrix production, adhesion and proliferation, and morphological changes after culture. The biological outputs correlated with tissue proteomics, and network analysis identified several proteins linked to cell function. Our methodology enables broad screening of ECMs to connect tissue-specific composition with biological activity, providing a new resource for biomaterials research and further understanding of regeneration and disease mechanisms.

Phosphoproteomic Profiling Reveals Epstein-Barr Virus Protein Kinase Integration of DNA Damage Response and Mitotic Signaling.

Epstein-Barr virus (EBV) is etiologically linked to infectious mononucleosis and several human cancers. EBV encodes a conserved protein kinase BGLF4 that plays a key role in the viral life cycle. To provide new insight into the host proteins regulated by BGLF4, we utilized stable isotope labeling by amino acids in cell culture (SILAC)-based quantitative proteomics to compare site-specific phosphorylation in BGLF4-expressing Akata B cells. Our analysis revealed BGLF4-mediated hyperphosphorylation of 3,046 unique sites corresponding to 1,328 proteins. Frequency analysis of these phosphosites revealed a proline-rich motif signature downstream of BGLF4, indicating a broader substrate recognition for BGLF4 than its cellular ortholog cyclin-dependent kinase 1 (CDK1). Further, motif analysis of the hyperphosphorylated sites revealed enrichment in ATM, ATR and Aurora kinase substrates while functional analyses revealed significant enrichment of pathways related to the DNA damage response (DDR), mitosis and cell cycle. Phosphorylation of proteins associated with the mitotic spindle assembly checkpoint (SAC) indicated checkpoint activation, an event that inactivates the anaphase promoting complex/cyclosome, APC/C. Furthermore, we demonstrated that BGLF4 binds to and directly phosphorylates the key cellular proteins PP1, MPS1 and CDC20 that lie upstream of SAC activation and APC/C inhibition. Consistent with APC/C inactivation, we found that BGLF4 stabilizes the expression of many known APC/C substrates. We also noted hyperphosphorylation of 22 proteins associated the nuclear pore complex, which may contribute to nuclear pore disassembly and SAC activation. A drug that inhibits mitotic checkpoint activation also suppressed the accumulation of extracellular EBV virus. Taken together, our data reveal that, in addition to the DDR, manipulation of mitotic kinase signaling and SAC activation are mechanisms associated with lytic EBV replication. All MS data have been deposited in the ProteomeXchange with identifier PXD002411 (http://proteomecentral.proteomexchange.org/dataset/PXD002411).

Phosphoproteomic Analysis Identifies Focal Adhesion Kinase 2 (FAK2) as a Potential Therapeutic Target for Tamoxifen Resistance in Breast Cancer.

Tamoxifen, an estrogen receptor-α (ER) antagonist, is an important agent for the treatment of breast cancer. However, this therapy is complicated by the fact that a substantial number of patients exhibit either de novo or acquired resistance. To characterize the signaling mechanisms underlying this resistance, we treated the MCF7 breast cancer cell line with tamoxifen for over six months and showed that this cell line acquired resistance to tamoxifen in vitro and in vivo. We performed SILAC-based quantitative phosphoproteomic profiling on the tamoxifen resistant and vehicle-treated sensitive cell lines to quantify the phosphorylation alterations associated with tamoxifen resistance. From >5600 unique phosphopeptides identified, 1529 peptides exhibited hyperphosphorylation and 409 peptides showed hypophosphorylation in the tamoxifen resistant cells. Gene set enrichment analysis revealed that focal adhesion pathway was one of the most enriched signaling pathways activated in tamoxifen resistant cells. Significantly, we showed that the focal adhesion kinase FAK2 was not only hyperphosphorylated but also transcriptionally up-regulated in tamoxifen resistant cells. FAK2 suppression by specific siRNA knockdown or a small molecule inhibitor repressed cellular proliferation in vitro and tumor formation in vivo. More importantly, our survival analysis revealed that high expression of FAK2 is significantly associated with shorter metastasis-free survival in estrogen receptor-positive breast cancer patients treated with tamoxifen. Our studies suggest that FAK2 is a potential therapeutic target for the management of hormone-refractory breast cancers.

Annotation of the zebrafish genome through an integrated transcriptomic and proteomic analysis.

Accurate annotation of protein-coding genes is one of the primary tasks upon the completion of whole genome sequencing of any organism. In this study, we used an integrated transcriptomic and proteomic strategy to validate and improve the existing zebrafish genome annotation. We undertook high-resolution mass-spectrometry-based proteomic profiling of 10 adult organs, whole adult fish body, and two developmental stages of zebrafish (SAT line), in addition to transcriptomic profiling of six organs. More than 7,000 proteins were identified from proteomic analyses, and ∼ 69,000 high-confidence transcripts were assembled from the RNA sequencing data. Approximately 15% of the transcripts mapped to intergenic regions, the majority of which are likely long non-coding RNAs. These high-quality transcriptomic and proteomic data were used to manually reannotate the zebrafish genome. We report the identification of 157 novel protein-coding genes. In addition, our data led to modification of existing gene structures including novel exons, changes in exon coordinates, changes in frame of translation, translation in annotated UTRs, and joining of genes. Finally, we discovered four instances of genome assembly errors that were supported by both proteomic and transcriptomic data. Our study shows how an integrative analysis of the transcriptome and the proteome can extend our understanding of even well-annotated genomes.

Quantitative phosphoproteomic analysis of IL-33-mediated signaling.

Interleukin-33 (IL-33) is a novel member of the IL-1 family of cytokines that plays diverse roles in the regulation of immune responses. IL-33 exerts its effects through a heterodimeric receptor complex resulting in the production and release of proinflammatory cytokines. A detailed understanding of the signaling pathways activated by IL-33 is still unclear. To gain insights into the IL-33-mediated signaling mechanisms, we carried out a SILAC-based global quantitative phosphoproteomic analysis that resulted in the identification of 7191 phosphorylation sites derived from 2746 proteins. We observed alterations in the level of phosphorylation in 1050 sites corresponding to 672 proteins upon IL-33 stimulation. We report, for the first time, phosphorylation of multiple protein kinases, including mitogen-activated protein kinase activated protein kinase 2 (Mapkapk2), receptor (TNFRSF) interacting serine-threonine kinase 1 (Ripk1), and NAD kinase (Nadk) that are induced by IL-33. In addition, we observed IL-33-induced phosphorylation of several protein phosphatases including protein tyrosine phosphatase, nonreceptor-type 12 (Ptpn12), and inositol polyphosphate-5-phosphatase D (Inpp5d), which have not been reported previously. Network analysis revealed an enrichment of actin binding and cytoskeleton reorganization that could be important in macrophage activation induced by IL-33. Our study is the first quantitative analysis of IL-33-regulated phosphoproteome. Our findings significantly expand the understanding of IL-33-mediated signaling events and have the potential to provide novel therapeutic targets pertaining to immune-related diseases such as asthma where dysregulation of IL-33 is observed. All MS data have been deposited in the ProteomeXchange with identifier PXD000984 (http://proteomecentral.proteomexchange.org/dataset/PXD000984).

Phosphoproteomic analysis reveals compensatory effects in the piriform cortex of VX nerve agent exposed rats.

To gain insights into the toxicity induced by the nerve agent VX, an MS-based phosphoproteomic analysis was carried out on the piriform cortex region of brains from VX-treated rats. Using isobaric tag based TMT labeling followed by titanium dioxide enrichment strategy, we identified 9975 unique phosphosites derived from 3287 phosphoproteins. Temporal changes in the phosphorylation status of peptides were observed over a time period of 24 h in rats exposed to a 1× LD50, intravenous (i.v.) dose with the most notable changes occurring at the 1 h postexposure time point. Five major functional classes of proteins exhibited changes in their phosphorylation status: (i) ion channels/transporters, including ATPases, (ii) kinases/phosphatases, (iii) GTPases, (iv) structural proteins, and (v) transcriptional regulatory proteins. This study is the first quantitative phosphoproteomic analysis of VX toxicity in the brain. Understanding the toxicity and compensatory signaling mechanisms will improve the understanding of the complex toxicity of VX in the brain and aid in the elucidation of novel molecular targets that would be important for development of improved countermeasures. All MS data have been deposited in the ProteomeXchange with identifier PXD001184 (http://proteomecentral.proteomexchange.org/dataset/PXD001184).

Functional annotation of proteome encoded by human chromosome 22.

As part of the chromosome-centric human proteome project (C-HPP) initiative, we report our progress on the annotation of chromosome 22. Chromosome 22, spanning 51 million base pairs, was the first chromosome to be sequenced. Gene dosage alterations on this chromosome have been shown to be associated with a number of congenital anomalies. In addition, several rare but aggressive tumors have been associated with this chromosome. A number of important gene families including immunoglobulin lambda locus, Crystallin beta family, and APOBEC gene family are located on this chromosome. On the basis of proteomic profiling of 30 histologically normal tissues and cells using high-resolution mass spectrometry, we show protein evidence of 367 genes on chromosome 22. Importantly, this includes 47 proteins, which are currently annotated as "missing" proteins. We also confirmed the translation start sites of 120 chromosome 22-encoded proteins. Employing a comprehensive proteogenomics analysis pipeline, we provide evidence of novel coding regions on this chromosome which include upstream ORFs and novel exons in addition to correcting existing gene structures. We describe tissue-wise expression of the proteins and the distribution of gene families on this chromosome. These data have been deposited to ProteomeXchange with the identifier PXD000561.

Proteomic analysis of human osteoarthritis synovial fluid.

BACKGROUND: Osteoarthritis is a chronic musculoskeletal disorder characterized mainly by progressive degradation of the hyaline cartilage. Patients with osteoarthritis often postpone seeking medical help, which results in the diagnosis being made at an advanced stage of cartilage destruction. Sustained efforts are needed to identify specific markers that might help in early diagnosis, monitoring disease progression and in improving therapeutic outcomes. We employed a multipronged proteomic approach, which included multiple fractionation strategies followed by high resolution mass spectrometry analysis to explore the proteome of synovial fluid obtained from osteoarthritis patients. In addition to the total proteome, we also enriched glycoproteins from synovial fluid using lectin affinity chromatography. RESULTS: We identified 677 proteins from synovial fluid of patients with osteoarthritis of which 545 proteins have not been previously reported. These novel proteins included ADAM-like decysin 1 (ADAMDEC1), alanyl (membrane) aminopeptidase (ANPEP), CD84, fibulin 1 (FBLN1), matrix remodelling associated 5 (MXRA5), secreted phosphoprotein 2 (SPP2) and spondin 2 (SPON2). We identified 300 proteins using lectin affinity chromatography, including the glycoproteins afamin (AFM), attractin (ATRN), fibrillin 1 (FBN1), transferrin (TF), tissue inhibitor of metalloproteinase 1 (TIMP1) and vasorin (VSN). Gene ontology analysis confirmed that a majority of the identified proteins were extracellular and are mostly involved in cell communication and signaling. We also confirmed the expression of ANPEP, dickkopf WNT signaling pathway inhibitor 3 (DKK3) and osteoglycin (OGN) by multiple reaction monitoring (MRM) analysis of osteoarthritis synovial fluid samples. CONCLUSIONS: We present an in-depth analysis of the synovial fluid proteome from patients with osteoarthritis. We believe that the catalog of proteins generated in this study will further enhance our knowledge regarding the pathophysiology of osteoarthritis and should assist in identifying better biomarkers for early diagnosis.

Whole exome sequencing uncovers HRAS mutations as potential mediators of resistance to metronomic chemotherapy.

OBJECTIVES: The aim of this pilot study is to identify the genetic factors that contribute to the response of metronomic chemotherapy in head and neck squamous cell carcinoma (HNSCC) patients using whole-exome sequencing (WES). This study would facilitate the identification of predictive biomarkers, which would enable personalized treatment strategies and improve treatment outcomes for patients with HNSCC. MATERIALS AND METHODS: We have selected patients with recurrent head and neck cancer who underwent metronomic chemotherapy. Sequential tumor biopsies were collected from the patients at different stages of treatment to capture the genomic alterations and tumor evolution during metronomic chemotherapy and sequenced using WES. RESULTS: We identified several known HNSCC hallmark genes reported in COSMIC, including KMT2B, NOTCH1, FAT1, TP53, HRAS, CASP8, and CDKN2A. Copy number alteration analysis revealed amplifications and deletions in several oncogenic and tumor suppressor genes. COSMIC Mutational Signature 15 associated with defective DNA mismatch repair was enriched in 73% of HNSCC samples. Further, the comparison of genomic alterations between responders and non-responders identified HRAS gene uniquely mutated in non-responders that could potentially contribute to resistance against metronomic chemotherapy. DISCUSSION: Our findings corroborate the molecular heterogeneity of recurrent HNSCC tumors and establish an association between HRAS mutations and resistance to metronomic chemotherapy, suggesting HRAS as a potential therapeutic target. Combining HRAS inhibitors with metronomic regimens could improve treatment sensitivity in HRAS-mutated HNSCC patients. Further studies are needed to fully elucidate the genomic mechanisms underlying the response to metronomic chemotherapy.

Pan-cancer proteomic map of 949 human cell lines.

The proteome provides unique insights into disease biology beyond the genome and transcriptome. A lack of large proteomic datasets has restricted the identification of new cancer biomarkers. Here, proteomes of 949 cancer cell lines across 28 tissue types are analyzed by mass spectrometry. Deploying a workflow to quantify 8,498 proteins, these data capture evidence of cell-type and post-transcriptional modifications. Integrating multi-omics, drug response, and CRISPR-Cas9 gene essentiality screens with a deep learning-based pipeline reveals thousands of protein biomarkers of cancer vulnerabilities that are not significant at the transcript level. The power of the proteome to predict drug response is very similar to that of the transcriptome. Further, random downsampling to only 1,500 proteins has limited impact on predictive power, consistent with protein networks being highly connected and co-regulated. This pan-cancer proteomic map (ProCan-DepMapSanger) is a comprehensive resource available at https://cellmodelpassports.sanger.ac.uk.

PIONEER: Pipeline for Generating High-Quality Spectral Libraries for DIA-MS Data.

Data-independent-acquisition mass spectrometry (DIA-MS) is a state-of-the-art proteomic technique for high-throughput identification and quantification of peptides and proteins. Interpretation of DIA-MS data relies on the use of a spectral library, which is optimally created from data acquired from the same samples in data-dependent acquisition (DDA) mode. As DIA-MS quantification relies on the spectral libraries, having a high-quality, non-redundant, and comprehensive spectral library is essential. This article describes the major steps for creating a high-quality spectral library using a combination of multiple complementary search engines. We discuss appropriate strategies to control the false discovery rate for the final spectral library as a result of merging multiple searches. © 2021 The Authors Current Protocols © 2021 Wiley Periodicals LLC. Basic Protocol 1: Searching DDA-MS files with multiple search engines Basic Protocol 2: Merging results from multiple search engines Basic Protocol 3: Creating spectral libraries from merged results Alternate Protocol: Using CLI for automating tasks Support Protocol: Creating concatenated FASTA files.

Comprehensive analysis of the secreted proteome of adult Necator americanus hookworms.

The human hookworm Necator americanus infects more than 400 million people worldwide, contributing substantially to the poverty in these regions. Adult stage N. americanus live in the small intestine of the human host where they inject excretory/secretory (ES) products into the mucosa. ES products have been characterized at the proteome level for a number of animal hookworm species, but until now, the difficulty in obtaining sufficient live N. americanus has been an obstacle in characterizing the secretome of this important human pathogen. Herein we describe the ES proteome of N. americanus and utilize this information along with RNA Seq data to conduct the first proteogenomic analysis of a parasitic helminth, significantly improving the available genome and thereby generating a robust description of the parasite secretome. The genome annotation resulted in a revised prediction of 3,425 fewer genes than initially reported, accompanied by a significant increase in the number of exons and introns, total gene length and the percentage of the genome covered by genes. Almost 200 ES proteins were identified by LC-MS/MS with SCP/TAPS proteins, 'hypothetical' proteins and proteases among the most abundant families. These proteins were compared to commonly used model species of human parasitic infections, including Ancylostoma caninum, Nippostrongylus brasiliensis and Heligmosomoides polygyrus. SCP/TAPS proteins are immunogenic in nematode infections, so we expressed four of those identified in this study in recombinant form and showed that they are all recognized to varying degrees by serum antibodies from hookworm-infected subjects from a disease-endemic area of Brazil. Our findings provide valuable information on important families of proteins with both known and unknown functions that could be instrumental in host-parasite interactions, including protein families that might be key for parasite survival in the onslaught of robust immune responses, as well as vaccine and diagnostic targets.

Strategies to enable large-scale proteomics for reproducible research.

Reproducible research is the bedrock of experimental science. To enable the deployment of large-scale proteomics, we assess the reproducibility of mass spectrometry (MS) over time and across instruments and develop computational methods for improving quantitative accuracy. We perform 1560 data independent acquisition (DIA)-MS runs of eight samples containing known proportions of ovarian and prostate cancer tissue and yeast, or control HEK293T cells. Replicates are run on six mass spectrometers operating continuously with varying maintenance schedules over four months, interspersed with ~5000 other runs. We utilise negative controls and replicates to remove unwanted variation and enhance biological signal, outperforming existing methods. We also design a method for reducing missing values. Integrating these computational modules into a pipeline (ProNorM), we mitigate variation among instruments over time and accurately predict tissue proportions. We demonstrate how to improve the quantitative analysis of large-scale DIA-MS data, providing a pathway toward clinical proteomics.

Proteomics of human aqueous humor.

The aqueous humor is a colorless, transparent fluid that fills the anterior chamber of the eye. It plays an important role in maintaining the intraocular pressure and providing nourishment to the lens and cornea. The constitution of the aqueous humor is controlled by the blood-aqueous barrier. Though this ocular fluid has been extensively studied, its role in ocular physiology is still not completely understood. In this study, aqueous humor samples were collected from 250 patients undergoing cataract surgery, subjected to multiple fractionation strategies and analyzed on a Fourier transform LTQ-Orbitrap Velos mass spectrometer. In all, we identified 763 proteins, of which 386 have been identified for the first time in this study. Sorbitol dehydrogenase (SORD), filensin (BFSP1), and phakinin (BFSP2) are some of the proteins that have not been previously reported in the aqueous humor. Gene Ontology analysis revealed 35% of the identified proteins to be extracellular, with a majority of them involved in cell communication and signal transduction. This study comprehensively reports 386 novel proteins that have important potential as biomarker candidates for future research into personalized medicine and diagnostics aimed towards improving visual health.

Identification of differentially expressed serum proteins in gastric adenocarcinoma.

Gastric adenocarcinoma is an aggressive cancer with poor prognosis. Blood based biomarkers of gastric cancer have the potential to improve diagnosis and monitoring of these tumors. Proteins that show altered levels in the circulation of gastric cancer patients could prove useful as putative biomarkers. We used an iTRAQ-based quantitative proteomic approach to identify proteins that show altered levels in the sera of patients with gastric cancer. Our study resulted in identification of 643 proteins, of which 48 proteins showed increased levels and 11 proteins showed decreased levels in serum from gastric cancer patients compared to age and sex matched healthy controls. Proteins that showed increased expression in gastric cancer included inter-alpha-trypsin inhibitor heavy chain H4 (ITIH4), Mannose-binding protein C (MBL2), sex hormone-binding globulin (SHBG), insulin-like growth factor-binding protein 2 (IGFBP2), serum amyloid A protein (SAA1), Orosomucoid 1 (ORM1) and extracellular superoxide dismutase [Cu-Zn] (SOD3). We used multiple reaction monitoring assays and validated elevated levels of ITIH4 and SAA1 proteins in serum from gastric cancer patients. BIOLOGICAL SIGNIFICANCE: Gastric cancer is a highly aggressive cancer associated with high mortality. Serum-based biomarkers are of considerable interest in diagnosis and monitoring of various diseases including cancers. Gastric cancer is often diagnosed at advanced stages resulting in poor prognosis and high mortality. Pathological diagnosis using biopsy specimens remains the gold standard for diagnosis of gastric cancer. Serum-based biomarkers are of considerable importance as they are minimally invasive. In this study, we carried out quantitative proteomic profiling of serum from gastric cancer patients to identify proteins that show altered levels in gastric cancer patients. We identified more than 50 proteins that showed altered levels in gastric cancer patient sera. Validation in a large cohort of well classified patient samples would prove useful in identifying novel blood based biomarkers for gastric cancers. This article is part of a Special Issue entitled: Proteomics in India.

A multi-omic analysis of human naïve CD4+ T cells.

BACKGROUND: Cellular function and diversity are orchestrated by complex interactions of fundamental biomolecules including DNA, RNA and proteins. Technological advances in genomics, epigenomics, transcriptomics and proteomics have enabled massively parallel and unbiased measurements. Such high-throughput technologies have been extensively used to carry out broad, unbiased studies, particularly in the context of human diseases. Nevertheless, a unified analysis of the genome, epigenome, transcriptome and proteome of a single human cell type to obtain a coherent view of the complex interplay between various biomolecules has not yet been undertaken. Here, we report the first multi-omic analysis of human primary naïve CD4+ T cells isolated from a single individual. RESULTS: Integrating multi-omics datasets allowed us to investigate genome-wide methylation and its effect on mRNA/protein expression patterns, extent of RNA editing under normal physiological conditions and allele specific expression in naïve CD4+ T cells. In addition, we carried out a multi-omic comparative analysis of naïve with primary resting memory CD4+ T cells to identify molecular changes underlying T cell differentiation. This analysis provided mechanistic insights into how several molecules involved in T cell receptor signaling are regulated at the DNA, RNA and protein levels. Phosphoproteomics revealed downstream signaling events that regulate these two cellular states. Availability of multi-omics data from an identical genetic background also allowed us to employ novel proteogenomics approaches to identify individual-specific variants and putative novel protein coding regions in the human genome. CONCLUSIONS: We utilized multiple high-throughput technologies to derive a comprehensive profile of two primary human cell types, naïve CD4+ T cells and memory CD4+ T cells, from a single donor. Through vertical as well as horizontal integration of whole genome sequencing, methylation arrays, RNA-Seq, miRNA-Seq, proteomics, and phosphoproteomics, we derived an integrated and comparative map of these two closely related immune cells and identified potential molecular effectors of immune cell differentiation following antigen encounter.