Whole-genome sequencing analysis of suicide deaths integrating brain-regulatory eQTLs data to identify risk loci and genes

Recent large-scale genome-wide association studies (GWAS) have started to identify potential genetic risk loci associated with risk of suicide; however, a large portion of suicide-associated genetic factors affecting gene expression remain elusive. Dysregulated gene expression, not assessed by GWAS, may play a significant role in increasing the risk of suicide death. We performed the first comprehensive genomic association analysis prioritizing brain expression quantitative trait loci (eQTLs) within regulatory regions in suicide deaths from the Utah Suicide Genetic Risk Study (USGRS). 440,324 brain-regulatory eQTLs were obtained by integrating brain eQTLs, histone modification ChIP-seq, ATAC-seq, DNase-seq, and Hi-C results from publicly available data. Subsequent genomic analyses were conducted in whole-genome sequencing (WGS) data from 986 suicide deaths of non-Finnish European (NFE) ancestry and 415 ancestrally matched controls. Additional independent USGRS suicide deaths with genotyping array data (n = 4657) and controls from the Genome Aggregation Database were explored for WGS result replication. One significant eQTL locus, rs926308 (p = 3.24e-06), was identified. The rs926308-T is associated with lower expression of RFPL3S, a gene important for neocortex development and implicated in arousal. Gene-based analyses performed using Sherlock Bayesian statistical integrative analysis also detected 20 genes with expression changes that may contribute to suicide risk. From analyzing publicly available transcriptomic data, ten of these genes have previous evidence of differential expression in suicide death or in psychiatric disorders that may be associated with suicide, including schizophrenia and autism (ZNF501, ZNF502, CNN3, IGF1R, KLHL36, NBL1, PDCD6IP, SNX19, BCAP29, and ARSA). Electronic health records (EHR) data was further merged to evaluate if there were clinically relevant subsets of suicide deaths associated with genetic variants. In summary, our study identified one risk locus and ten genes associated with suicide risk via gene expression, providing new insight into possible genetic and molecular mechanisms leading to suicide.

Extended familial risk of suicide death is associated with younger age at death and elevated polygenic risk of suicide

Suicide accounts for >800,000 deaths annually worldwide; prevention is an urgent public health issue. Identification of risk factors remains challenging due to complexity and heterogeneity. The study of suicide deaths with increased extended familial risk provides an avenue to reduce etiological heterogeneity and explore traits associated with increased genetic liability. Using extensive genealogical records, we identified high-risk families where distant relatedness of suicides implicates genetic risk. We compared phenotypic and polygenic risk score (PRS) data between suicides in high-risk extended families (high familial risk (HFR), n = 1,634), suicides linked to genealogical data not in any high-risk families (low familial risk (LFR), n = 147), and suicides not linked to genealogical data with unknown familial risk (UFR, n = 1,865). HFR suicides were associated with lower age at death (mean = 39.34 years), more suicide attempts, and more PTSD and trauma diagnoses. For PRS tests, we included only suicides with >90% European ancestry and adjusted for residual ancestry effects. HFR suicides showed markedly higher PRS of suicide death (calculated using cross-validation), supporting specific elevation of genetic risk of suicide in this subgroup, and also showed increased PRS of PTSD, suicide attempt, and risk taking. LFR suicides were substantially older at death (mean = 49.10 years), had fewer psychiatric diagnoses of depression and pain, and significantly lower PRS of depression. Results suggest extended familiality and trauma/PTSD may provide specificity in identifying individuals at genetic risk for suicide death, especially among younger ages, and that LFR of suicide warrants further study regarding the contribution of demographic and medical risks.

A support vector machine for identification of single-nucleotide polymorphisms from next-generation sequencing data

MOTIVATION

Accurate determination of single-nucleotide polymorphisms (SNPs) from next-generation sequencing data is a significant challenge facing bioinformatics researchers. Most current methods use mechanistic models that assume nucleotides aligning to a given reference position are sampled from a binomial distribution. While such methods are sensitive, they are often unable to discriminate errors resulting from misaligned reads, sequencing errors or platform artifacts from true variants.

RESULTS

To enable more accurate SNP calling, we developed an algorithm that uses a trained support vector machine (SVM) to determine variants from .BAM or .SAM formatted alignments of sequence reads. Our SVM-based implementation determines SNPs with significantly greater sensitivity and specificity than alternative platforms, including the UnifiedGenotyper included with the Genome Analysis Toolkit, samtools and FreeBayes. In addition, the quality scores produced by our implementation more accurately reflect the likelihood that a variant is real when compared with those produced by the Genome Analysis Toolkit. While results depend on the model used, the implementation includes tools to easily build new models and refine existing models with additional training data.

AVAILABILITY

Source code and executables are available from github.com/brendanofallon/SNPSVM/

Protein and steroid profiles in follicular fluid after ovarian hyperstimulation as potential biomarkers of IVF outcome

Controlled ovarian hyperstimulation is performed to assist with generation of multiple mature oocytes for use in in vitro fertilization (IVF). The goal of our study was to evaluate differences in protein and steroid profiles in ovarian follicular fluid (hFF) samples obtained during oocyte retrieval from women undergoing IVF treatment and to identify physiological pathways associated with the proteins. The hFF samples were depleted of abundant proteins, fractionated by ultrafiltration, digested, and analyzed by nano-LC-QTOF. Concentrations of 15 endogenous steroids were determined in the samples using LC-MS/MS methods. The total number of proteins identified in the samples was 75, of which 4, 7, and 2 were unique to the samples from women with viable pregnancy, miscarriage, and no pregnancy, respectively. Identified proteins were associated with the acute response signaling, coagulation system, intrinsic and extrinsic prothrombin activation, complement system, neuroprotective role of THOP1, FXR/RXR activation, role of tissue factor, and growth hormone pathways. A greater number of proteins associated with biosynthesis was found in hFF samples corresponding to the oocytes resulting in pregnancy. The abundance of seven proteins was found to be associated with steroidogenesis. The obtained data will contribute to better understanding of the pathogenesis and development of noninvasive markers for assessment of oocytes viability.

Consensus: a framework for evaluation of uncertain gene variants in laboratory test reporting

Accurate interpretation of gene testing is a key component in customizing patient therapy. Where confirming evidence for a gene variant is lacking, computational prediction may be employed. A standardized framework, however, does not yet exist for quantitative evaluation of disease association for uncertain or novel gene variants in an objective manner. Here, complementary predictors for missense gene variants were incorporated into a weighted Consensus framework that includes calculated reference intervals from known disease outcomes. Data visualization for clinical reporting is also discussed.

Utility of gene-specific algorithms for predicting pathogenicity of uncertain gene variants

The rapid advance of gene sequencing technologies has produced an unprecedented rate of discovery of genome variation in humans. A growing number of authoritative clinical repositories archive gene variants and disease phenotypes, yet there are currently many more gene variants that lack clear annotation or disease association. To date, there has been very limited coverage of gene-specific predictors in the literature. Here the evaluation is presented of "gene-specific" predictor models based on a naïve Bayesian classifier for 20 gene-disease datasets, containing 3986 variants with clinically characterized patient conditions. The utility of gene-specific prediction is then compared with "all-gene" generalized prediction and also with existing popular predictors. Gene-specific computational prediction models derived from clinically curated gene variant disease datasets often outperform established generalized algorithms for novel and uncertain gene variants.

Utility of characteristic QTOF MS/MS fragmentation for MHC class I peptides

Systematic investigation of cellular process by mass spectrometric detection of peptides obtained from proteins digestion or directly from immuno-purification can be a powerful tool when used appropriately. The true sequence of these peptides is defined by the interpretation of spectral data using a variety of available algorithms. However peptide match algorithm scoring is typically based on some, but not all, of the mechanisms of peptide fragmentation. Although algorithm rules for soft ionization techniques generally fit very well to tryptic peptides, manual validation of spectra is often required for endogenous peptides such as MHC class I molecules where traditional trypsin digest techniques are not used. This study summarizes data mining and manual validation of hundreds of peptide sequences from MHC class I molecules in publically available data files. We herein describe several important features to improve and quantify manual validation for these endogenous peptides--post automated algorithm searching. Important fragmentation patterns are discussed for the studied MHC Class I peptides. These findings lead to practical rules that are helpful when performing manual validation. Furthermore, these observations may be useful to improve current peptide search algorithms or development of novel software tools.

Predicting phenotypic severity of uncertain gene variants in the RET proto-oncogene

Although reported gene variants in the RET oncogene have been directly associated with multiple endocrine neoplasia type 2 and hereditary medullary thyroid carcinoma, other mutations are classified as variants of uncertain significance (VUS) until the associated clinical phenotype is made clear. Currently, some 46 non-synonymous VUS entries exist in curated archives. In the absence of a gold standard method for predicting phenotype outcomes, this follow up study applies feature selected amino acid physical and chemical properties feeding a Bayes classifier to predict disease association of uncertain gene variants into categories of benign and pathogenic. Algorithm performance and VUS predictions were compared to established phylogenetic based mutation prediction algorithms. Curated outcomes and unpublished RET gene variants with known disease association were used to benchmark predictor performance. Reliable classification of RET uncertain gene variants will augment current clinical information of RET mutations and assist in improving prediction algorithms as knowledge increases.

The Alport syndrome COL4A5 variant database

Alport Syndrome is a progressive renal disease with cochlear and ocular involvement. The most common form ( approximately 80%) is inherited in an X-linked pattern. X-linked Alport Syndrome (XLAS) is caused by mutations in the type IV collagen alpha chain 5 (COL4A5). We have developed a curated disease-specific database containing reported sequence variants in COL4A5. Currently the database archives a total of 520 sequence variants, verified for their position within the COL4A5 gene and named following standard nomenclature. Sequence variants are reported with accompanying information on protein effect, classification of mutation vs. polymorphism, mutation type based on the first description in the literature, and links to pertinent publications. In addition, features of this database include disease information, relevant links for Alport syndrome literature, reference sequence information, and ability to query by various criteria. On-line submission for novel gene variants or updating information on existing database entries is also possible. This free online scientific resource was developed with the clinical laboratory in mind to serve as a reference and repository for COL4A5 variants.

Endoplasmic reticulum aminopeptidase associated with antigen processing defines the composition and structure of MHC class I peptide repertoire in normal and virus-infected cells

The MHC class I (MHC-I) molecules ferry a cargo of peptides to the cell surface as potential ligands for CD8(+) cytotoxic T cells. For nearly 20 years, the cargo has been described as a collection of short 8-9 mer peptides, whose length and sequences were believed to be primarily determined by the peptide-binding groove of MHC-I molecules. Yet the mechanisms for producing peptides of such optimal length and composition have remained unclear. In this study, using mass spectrometry, we determined the amino acid sequences of a large number of naturally processed peptides in mice lacking the endoplasmic reticulum aminopeptidase associated with Ag processing (ERAAP). We find that ERAAP-deficiency changed the oeuvre and caused a marked increase in the length of peptides normally presented by MHC-I. Furthermore, we observed similar changes in the length of viral peptides recognized by CD8(+) T cells in mouse CMV-infected ERAAP-deficient mice. In these mice, a distinct CD8(+) T cell population was elicited with specificity for an N-terminally extended epitope. Thus, the characteristic length, as well as the composition of MHC-I peptide cargo, is determined not only by the MHC-I peptide-binding groove but also by ERAAP proteolysis in the endoplasmic reticulum.

Multiple endocrine neoplasia type 2 RET protooncogene database: repository of MEN2-associated RET sequence variation and reference for genotype/phenotype correlations

Multiple endocrine neoplasia type 2 (MEN2) is an inherited, autosomal-dominant disorder caused by deleterious mutations within the RET protooncogene. MEN2 RET mutations are mainly heterozygous, missense sequence changes found in RET exons 10, 11, and 13-16. Our group has developed the publicly available, searchable MEN2 RET database to aid in genotype/phenotype correlations, using Human Genome Variation Society recommendations for sequence variation nomenclature and database content. The MEN2 RET database catalogs all RET sequence variation relevant to the MEN2 syndromes, with associated clinical information. Each database entry lists a RET sequence variation's location within the RET gene, genotype, pathogenicity classification, MEN2 phenotype, first literature reference, and comments (which may contain information on other clinical features, complex genotypes, and additional literature references). The MEN2 phenotype definitions were derived from the International RET Mutation Consortium guidelines for classification of MEN2 disease phenotypes. Although nearly all of the 132 RET sequence variation entries initially cataloged in the database were from literature reports, novel sequence variation and updated phenotypic information for any existing database entry can be submitted electronically on the database website. The database website also contains links to selected MEN2 literature reviews, gene and protein information, and RET reference sequences. The MEN2 RET database (www.arup.utah.edu/database/MEN2/MEN2_welcome.php) will serve as a repository for MEN2-associated RET sequence variation and reference for RET genotype/MEN2 phenotype correlations.

A depletion strategy for improved detection of human proteins from urine

With rapidly growing interest in the urine proteome, methods for reducing sample complexity are becoming increasingly important. Depletion strategies for removal of high-abundance proteins from human urine have not been reported. A commercial kit designed for depletion of abundant proteins from plasma was evaluated for removing top proteins from urine of patients with proteinuria. The number of low-abundance proteins identified in urine after depletion increased nearly 2.5-fold.

Analysis of gene expression profile of TPM3-ALK positive anaplastic large cell lymphoma reveals overlapping and unique patterns with that of NPM-ALK positive anaplastic large cell lymphoma

Anaplastic large cell lymphoma (ALCL) comprises a group of non-Hodgkin lymphomas characterized by the expression of the CD30/Ki-1 antigen. A subset of ALCL is characterized by chromosomal translocations involving the anaplastic lymphoma kinase (ALK) gene on chromosome 2. While the most common translocation is the t(2;5)(p23;q35) involving the nucleophosmin (NPM) gene on chromosome 5, up to 12 other translocations partners of the ALK gene have been identified. One of these is the t(1;2)(q25;p23) which results in the formation of the chimeric protein TPM3-ALK. While several of the signaling pathways induced by NPM-ALK have been elucidated, those involved in ALCLs harboring TPM3-ALK are largely unknown. In order to investigate the expression profiles of ALCLs carrying the NPM-ALK and TPM3-ALK fusions, we carried out cDNA microarray analysis of two ALCL tissue samples, one expressing the NPM-ALK fusion protein and the other the TPM3-ALK fusion protein. RNA was extracted from snap-frozen tissues, labeled with fluorescent dyes and analyzed using cDNAs microarray containing approximately 9,200 genes and expressed sequence tags (ESTs). Quantitative fluorescence RT-PCR was performed to validate the cDNA microarray data on nine selected gene targets. Our results show a significant overlap of genes deregulated in the NPM-ALK and TPM-ALK positive lymphomas. These deregulated genes are involved in diverse cellular functions, such as cell cycle regulation, apoptosis, proliferation, and adhesion. Interestingly, a subset of the genes was distinct in their expression pattern in the two types of lymphomas. More importantly, many genes that were not previously associated with ALK positive lymphomas were identified. Our results demonstrate the overlapping and unique transcriptional patterns associated with the NPM-ALK and TPM3-ALK fusions in ALCL.

Comprehensive identification of proteins in Hodgkin lymphoma-derived Reed-Sternberg cells by LC-MS/MS

Mass spectrometry-based proteomics in conjunction with liquid chromatography and bioinformatics analysis provides a highly sensitive and high-throughput approach for the identification of proteins. Hodgkin lymphoma is a form of malignant lymphoma characterized by the proliferation of Reed-Sternberg cells and background reactive lymphocytes. Comprehensive analysis of proteins expressed and released by Reed-Sternberg cells would assist in the discovery of potential biomarkers and improve our understanding of its pathogenesis. The subcellular proteome of the three cellular compartments from L428 and KMH2 Hodgkin lymphoma-derived cell lines were fractionated, and analyzed by reverse-phase liquid chromatography coupled with electrospray ionization tandem mass spectrometry. Additionally, proteins released by Hodgkin lymphoma-derived L428 cells were extracted from serum-free culture media and analyzed. Peptide spectra were analyzed using TurboSEQUEST against the UniProt protein database (5.26.05; 188 712 entries). A subset of the identified proteins was validated by Western blot analysis, immunofluorescence microscopy and immunohistochemistry. A total of 1945 proteins were identified with 785 from the cytosolic fraction, 305 from the membrane fraction, 441 from the nuclear fraction and 414 released proteins using a minimum of two peptide identifications per protein and an error rate of <5.0%. Identification of proteins from diverse functional groups reflected the functional complexity of the Reed-Sternberg proteome. Proteins with previously reported oncogenic function in other cancers and from signaling pathways implicated in Hodgkin lymphoma were identified. Selected proteins without previously demonstrated expression in Hodgkin lymphoma were validated by Western blot analysis (B-RAF, Erb-B3), immunofluorescence microscopy (Axin1, Tenascin-X, Mucin-2) and immunohistochemistry using a tissue microarray (BRAF, PIM1). This study represents the first comprehensive inventory of proteins expressed by Reed-Sternberg cells of Hodgkin lymphoma and demonstrates the utility of combining cellular subfractionation, protein precipitation, tandem mass spectrometry and bioinformatics analysis for comprehensive identification of proteins that may represent potential biomarkers of the disease.

Mutation database for the galactose-1-phosphate uridyltransferase (GALT) gene

Classical galactosemia is an autosomal recessive disorder caused by mutations in the galactose-1-phosphate uridyltransferase (GALT) gene. Our group developed a disease-specific database containing all of the reported sequence variants in GALT (Available at: http://arup.utah.edu/database/galactosemia/GALT_welcome.php; Last accessed: 13 April 2007). Currently the database contains a total of 229 sequence variants, of which 196 are mutations (including nine novel mutations identified in our laboratory), 31 polymorphisms in both introns and exons, and two variants of unknown or uncertain significance. All sequence variants have been verified for their position within the GALT gene and named following standard nomenclature. Sequence variants are reported with accompanying information on protein effect, classification of mutation vs. polymorphism, mutation type (when applicable) based on how each was first described in the literature, and accompanying link to pertinent publication. Unpublished variants are described with relevant clinical information that supports their classification as causative of the disease vs. polymorphisms. Other features of this database include disease information, relevant links for galactosemia and literature, reference sequences, ability to query by various criteria, and submit of novel variations to the database. This free online scientific resource was developed with the clinical laboratory in mind to serve as a reference and repository for novel findings that are periodically collected, verified, and updated into the database.

Proteome-wide changes induced by the Hsp90 inhibitor, geldanamycin in anaplastic large cell lymphoma cells

The molecular chaperone heat shock protein 90 (Hsp90) affects the function of many oncogenic signaling proteins including nucleophosmin-anaplastic lymphoma kinase (NPM-ALK) expressed in anaplastic large cell lymphoma (ALCL). While ALK-positive ALCL cells are sensitive to the Hsp90 inhibitor and the geldanamycin (GA) analog, 17-allylamino-17-demethoxygeldanamycin (17-AAG), the proteomic effects of these drugs on ALK-positive ALCL cells are unpublished. In this study, we investigated the cellular, biologic, and proteomic changes occurring in ALK-positive ALCL cells in response to GA treatment. GA induced G2/M cell cycle arrest and caspase-3-mediated apoptosis. Furthermore, quantitative proteomic changes analyzed by cleavable isotope-coded affinity tag-LC-MS/MS (cICAT-LC-MS/MS) identified 176 differentially expressed proteins. Out of these, 49 were upregulated 1.5-fold or greater and 70 were downregulated 1.5-fold or greater in GA-treated cells. Analysis of biological functions of differentially expressed proteins revealed diverse changes, including induction of proteins involved in the 26S proteasome as well as downregulation of proteins involved in signal transduction and protein and nucleic acid metabolism. Pathway analysis revealed changes in MAPK, WNT, NF-kappaB, TGFbeta, PPAR, and integrin signaling components. Our studies reveal some of the molecular and proteomic consequences of Hsp90 inhibition in ALK-positive ALCL cells and provide novel insights into the mechanisms of its diverse cellular effects.

Global proteome profiling of NPM/ALK-positive anaplastic large cell lymphoma

OBJECTIVE

Constitutive overexpression of nucleophosmin/anaplastic lymphoma kinase (NPM/ALK) is a key oncogenic event in anaplastic large cell lymphomas (ALCL) that carry the t(2;5)(p23;q35) translocation. Global proteomic analysis of NPM/ALK-positive ALCL would improve understanding of the disease pathogenesis and yield new candidate targets for novel treatment and diagnostic strategies.

MATERIALS AND METHODS

To comprehensively determine the inventory of proteins from NPM/ALK-positive ALCL SUDHL-1 cells, the membrane, cytoplasm, and nuclear subcellular fractions were resolved by one-dimensional sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) and analyzed using liquid chromatography-tandem mass spectrometry (LC-MS/MS). The MS spectra were interpreted using SEQUEST to search the electronic UniProt protein database, and analyzed by ProteinProphet and INTERACT.

RESULTS

A total of 623 proteins consisting of 210 membrane, 229 cytoplasm, and 184 nuclear proteins were identified with a <or=5% error rate. Extensive annotation and systematic examination of the literature for information on 209 representative proteins indicated that 19.9% were reported to be expressed in T cells and 44.7% were reported to have important function in cancers, while only 4.3% were reported to be involved in ALCL pathogenesis. Categorization of proteins into functional groups was performed using GOMiner. A subset of the identified proteins was confirmed by Western blots and immunohistochemistry of tissue samples.

CONCLUSION

We present an extensive catalog of proteins expressed by NPM/ALK-positive ALCL. This study illustrates the potential for novel pathogenetic discovery in NPM/ALK-positive ALCL and the utility of combining cellular subfractionation, 1D SDS-PAGE, and LC-MS/MS for the comprehensive protein analysis of lymphoma.

Evaluation of enrichment techniques for mass spectrometry: identification of tyrosine phosphoproteins in cancer cells

Phosphorylation of tyrosine residues by protein tyrosine kinases mediates numerous cellular processes. Deregulated tyrosine phosphorylation underlies constitutive activation of signaling pathways leading to oncogenesis. Analytical techniques for evaluation of the global phosphoproteome level are challenging and can be improved on to enhance yields. Here, we evaluated several approaches to enrich for tyrosine phosphoproteins in cancer cells for subsequent liquid chromatography-tandem mass spectrometry analysis using lysates from SU-DHL-1 cells, which express the nucleophosmin-anaplastic lymphoma kinase tyrosine kinase as a model system. Cells were grown in the presence or absence of the phosphatase inhibitor sodium orthovanadate, and tyrosine phosphoproteins were subsequently enriched by immunoprecipitation or immunoaffinity chromatography and protein identification performed by liquid chromatography-tandem mass spectrometry. Our results show that sodium orthovanadate improves enrichment and thus detection of tyrosine phosphoproteins. Immunoprecipitation of tyrosine phosphoproteins using two different antiphosphotyrosine antibodies increased the number of protein identifications. Finally, peptides from proteins enriched by immunoprecipitation were more abundant (n=338) than those enriched by immunoaffinity chromatography (n=138), and relatively few proteins were found in common (n=43). Our data demonstrate the utility of an enrichment strategy for the mass spectrometry-based identification of tyrosine phosphoproteins and show the advantage of complementary techniques for greater protein identification.

Analytical characteristics of cleavable isotope-coded affinity tag-LC-tandem mass spectrometry for quantitative proteomic studies

Quantitative proteomic studies using cleavable isotope-coded affinity tags (cICAT) in concert with tandem mass spectrometry (MS/MS) permit unbiased comparisons between biologically distinct samples. We sought to determine the analytical characteristics of cICAT-based studies by examining the cumulative results of multiple, separate cICAT-based experiments involving human lymphoma-derived cells. We found that the number of identified proteins increased with larger numbers of fractions analyzed. The majority of proteins were identified by single peptides. Only 24 to 41% of the peptides contained cysteine residues, but 85% of the cysteine-containing peptides yielded quantification data. Approximately 28% of all identified proteins yielded quantification data, with 57% of these being differentially expressed by at least 1.5-fold. The quantification ratios of peptides for proteins with multiple quantified peptides were concordant in trend in 87% of instances. cICAT-labeled peptides identified proteins in all subcellular compartments without significant bias. Analysis of the flow-through fraction did not increase the number of peptides identified per protein. Our studies indicate that cICAT-LC-MS/MS yields quantifications primarily based on single peptides, and analysis of flow-through peptides does not contribute signifi-cantly to the results. Nevertheless, identifications based on single cICAT-labeled peptides with tryptic ends provide sufficiently reliable protein identifications and quantification information in cICAT-LC-MS/MS-based proteomic studies.

Identification of proteins released by follicular lymphoma-derived cells using a mass spectrometry-based approach

The recent advent of mass spectrometry-based methodologies for the analysis of complex protein mixtures opens new opportunities for the discovery of biomarkers that may aid in the diagnostic work-up of cancer. Follicular lymphoma (FL) is the most common form of low-grade non-Hodgkin lymphoma in the Western Hemisphere. Identification of tumor markers that facilitate early disease detection would be a significant advance in the management of FL. We have employed a strategy that entailed propagation of a follicular-derived cell line in serum-free media, protein extraction, and reverse-phase liquid chromatography, with subsequent electrospray ionization and tandem mass spectrometry analysis for the identification of proteins that are released by FL. Using a two-peptide minimum per protein and standard criteria, 209 proteins (5.6% maximum predicted error rate) released from the FL cells were identified. The released proteins included several growth factors, cytokines, acute phase reactants and cellular components previously known to be present in FL cells. Importantly, a greater proportion of proteins previously unassociated with FL were identified with high statistical confidence. Our studies provide a list of proteins, which may be candidates for early screening, diagnosis and therapeutic monitoring of patients with a suspected or biopsy-confirmed diagnosis of FL.

Proteomic identification of oncogenic chromosomal translocation partners encoding chimeric anaplastic lymphoma kinase fusion proteins

The anaplastic lymphoma kinase (ALK) on 2p23 is a tyrosine kinase that forms chimeric fusions with numerous translocation partners. We describe a mass spectrometry-based approach for the identification of ALK fusion partners. This approach accurately identified the nucleophosmin (NPM)-ALK fusion protein in an anaplastic large cell lymphoma (ALCL)-derived cell line carrying the t(2;5)(p23;q35), and the TPM3-ALK in a clinical biopsy of inflammatory myofibroblastic tumor (IMT) carrying the t(1;2)(q21;p23). This study shows the ability of mass spectrometry to identify oncogenic chimeric proteins resulting from chromosomal rearrangements. This strategy can be adapted for the identification of known and unknown translocation partners of chimeric ALK fusion proteins involved in oncogenesis.

Sequence alignment by cross-correlation

Many recent advances in biology and medicine have resulted from DNA sequence alignment algorithms and technology. Traditional approaches for the matching of DNA sequences are based either on global alignment schemes or heuristic schemes that seek to approximate global alignment algorithms while providing higher computational efficiency. This report describes an approach using the mathematical operation of cross-correlation to compare sequences. It can be implemented using the fast fourier transform for computational efficiency. The algorithm is summarized and sample applications are given. These include gene sequence alignment in long stretches of genomic DNA, finding sequence similarity in distantly related organisms, demonstrating sequence similarity in the presence of massive (approximately 90%) random point mutations, comparing sequences related by internal rearrangements (tandem repeats) within a gene, and investigating fusion proteins. Application to RNA and protein sequence alignment is also discussed. The method is efficient, sensitive, and robust, being able to find sequence similarities where other alignment algorithms may perform poorly.

Annotated proteome of a human T-cell lymphoma

As the reliable identification of proteins by tandem mass spectrometry becomes increasingly common, the full characterization of large data sets of proteins remains a difficult challenge. Our goal was to survey the proteome of a human T-cell lymphoma-derived cell line in a single set of experiments and present an automated method for the annotation of lists of proteins. A downstream application of these data includes the identification of novel pathogenetic and candidate diagnostic markers of T-cell lymphoma. Total protein isolated from cytoplasmic, membrane, and nuclear fractions of the SUDHL-1 T-cell lymphoma cell line was resolved by SDS-PAGE, and the entire gel lanes digested and analyzed by tandem mass spectrometry. Acquired data files were searched against the UniProt protein database using the SEQUEST algorithm. Search results for each subcellular fraction were analyzed using INTERACT and ProteinProphet. All protein identifications with an error rate of less than 10% were directly exported into excel and analyzed using GOMiner (NIH/NCI). The Gene ontology molecular function and cell location data were summarized for the identified proteins and results exported as user-interactive directed acyclic graphs. A total of 1105 unique proteins were identified and fully annotated, including numerous proteins that had not been previously characterized in lymphoma, in functional categories such as cell adhesion, migration, signaling, and stress response. This study demonstrates the utility of currently available bioinformatics tools for the robust identification and annotation of large numbers of proteins in a batchwise fashion.

Analysis of BCL6-interacting Proteins by Tandem Mass Spectrometry

B-cell lymphoma 6 (BCL6) is a 95-kDa nuclear phosphoprotein and member of the Pox virus zinc finger/bric-a-brac, tramtrack, broad complex (POZ/BTB) family of transcription factors. BCL6 is a transcriptional repressor required for germinal center formation, and the gene encoding it is frequently altered in diffuse large B-cell and follicular lymphomas. The dysregulation of BCL6 has therefore been implicated in lymphomagenesis. A limited number of proteins is known to interact with BCL6 and modulate its activity or participate in its role in transcriptional regulation. Identification of additional BCL6-binding proteins could reveal potential signaling targets and previously undescribed functional roles for BCL6. We used a functional proteomic approach to determine the identity of proteins that interact with BCL6. Proteins were isolated by co-immunoprecipitation with an anti-BCL6 antibody and identified using MS/MS. We identified 61 proteins in the BCL6 immunocomplex from the following Gene Ontology categories: transcription regulator activity (n = 18), binding activity (n = 11), signal transducer activity (n = 10), catalytic activity (n = 8), structural molecule activity (n = 3), enzyme regulator activity (n = 3), transporter activity (n = 2), motor activity (n = 2), chaperone activity (n = 1), and unknown function (n = 3). Importantly we identified BCL6 and several previously reported BCL6-interacting proteins in the BCL6 immunocomplex. The remaining proteins have not been shown previously to be associated with BCL6. MS/MS results were validated on four proteins using immunoprecipitation and Western blotting. Two of these protein interactions were further confirmed by reciprocal immunoprecipitation. This study demonstrates the utility of antibody immunoprecipitation and subsequent peptide identification by MS/MS for the elucidation of BCL6-binding proteins. Many of the novel proteins identified in this study suggest additional functional roles for BCL6 beyond transcriptional repression.

Identification of proteins from formalin-fixed paraffin-embedded cells by LC-MS/MS

There exists a need for robust approaches for tandem mass spectrometry (MS/MS)-based identification of proteins in formalin-fixed paraffin-embedded (FFPE) material. We demonstrate herein the identification of proteins in FFPE material using enzymatic cleavage for extraction of peptides from the FFPE specimen and liquid chromatography (LC) followed by MS/MS. We identified 324 proteins from a 3-year-old FFPE cell-block of a human lymphoma cell line. The identified proteins were assigned to the membrane, cytosol and nucleus, with diverse cellular functions. The results were comparable to those obtained with lysates from a fresh specimen of the lymphoma cell line. Western blotting analysis and immunofluorescence microscopy confirmed the expression of selected proteins. The functional significance of one protein (PKC eta) was validated using a PKC inhibitory peptide which resulted in lymphoma cell death in vitro. The ability to identify proteins from FFPE specimens has significant implications for MS/MS-based proteomics of vast repositories of archival primary tissue samples for disease-related discovery research.

Analysis of phosphatase and tensin homolog tumor suppressor interacting proteins by in vitro and in silico proteomics

The phosphatase and tensin homolog (PTEN) tumor suppressor is a multifunctional protein deregulated in many types of cancer. To date, a comprehensive documentation of PTEN interacting proteins has not been performed. The goal of our study was to characterize the PTEN interactome using affinity pull-down and tandem mass spectrometry (MS/MS). Wild-type PTEN cDNA was inserted into pTRC-His2 vector to create a 6-His tagged protein, which was expressed in Escherichia coli. Lysate from a human lymphoma cell line was used in pull-down assays, utilizing affinity for nickel-agarose beads. Bound proteins were eluted with imidazole, digested and analyzed on an LCQ DecaXP ion trap mass spectrometer. The nickel affinity pull-down efficiency was evaluated by sodium dodecyl sulfate-polyacrylamide gel electrophoresis and Western blot analysis. Acquired data were searched against the NCBI nr.fasta nonredundant protein database using the SEQUEST algorithm and screened using INTERACT and ProteinProphet. All experiments were performed in duplicate with 6-His-lacZ serving as control. A total of 79 proteins were identified in the wild-type 6-His-PTEN pull-down by MS/MS. We further validated a subset of the proteins present in the PTEN interactome by performing immunoprecipitation using an anti-PTEN antibody and establishing the presence of the proteins in the immunocomplex by Western blot analysis. A search of published PTEN interactions was also performed using Online Mendelian Inheritance in Man, Human Protein Reference Database, the IntAct Project database, and PubMed. This in silico analysis confirmed 42 out of 79 (53%) of the proteins identified by MS/MS. The remaining 37 proteins represent probable PTEN interactions not previously documented in public databases or reported in the literature. These results highlight the value of combining both in vitro biochemical approaches with in silico analyses for a comprehensive study of protein-protein interactions.

Quantitative proteomic and transcriptional analysis of the response to the p38 mitogen-activated protein kinase inhibitor SB203580 in transformed follicular lymphoma cells

The p38 mitogen-activated protein kinase (MAPK) is a key mediator of stress, extracellular-, growth factor-, and cytokine-induced signaling, and has been implicated in the development of cancer. Our previous work showed evidence for p38 MAPK activation in a subset of transformed follicular lymphomas (Elenitoba-Johnson et al. (2003) Proc. Natl. Acad. Sci. U.S.A. 100, 7259). We demonstrated that inhibition of p38 MAPK by SB203580 resulted in dose- and time-dependent caspase-3-mediated apoptosis. In order to further elucidate the basis of the cellular effects of SB203580, we have employed a systems biologic approach involving cDNA microarray and quantitative proteomic analysis of transformed follicular lymphoma derived-cells (OCI Ly-1) treated with SB203580. Gene expression profiling revealed differential expression (>/=1.5-fold) of 374 genes/ESTs in cells treated for 3 h and 515 genes/ESTs in cells treated for 21 h. The majority (52% at 3 h and 91% at 21 h) were down-regulated, including genes encoding growth cytokines, transcriptional regulators and cytoskeletal proteins. Quantitative proteomic analysis using ICAT-LC-MS/MS identified 277 differentially expressed proteins at 3 h and 350 proteins at 21 h of treatment with SB203580, the majority of which were also down-regulated. Analysis of functional groups of the differentially expressed proteins implicated components of diverse overlapping pathways including the IL-6/phosphatidylinositol 3-kinase, insulin-like growth factor 2/Ras/Raf, WNT8d/Frizzled, MAPK-activated protein kinase 2, and nuclear factor kappaB. The differential phosphorylation status of selected kinase-active proteins was validated by Western blotting analysis. Our complementary genomic and proteomic approach reveal the global cellular consequences of SB203580 treatment and provide insights into its growth inhibitory effect on transformed follicular lymphoma cells.

Identification of NPM-ALK interacting proteins by tandem mass spectrometry

Constitutive overexpression of nucleophosmin-anaplastic lymphoma kinase (NPM-ALK) is a key oncogenic event in anaplastic large-cell lymphomas with the characteristic chromosomal aberration t(2;5)(p23;q35). Proteins that interact with ALK tyrosine kinase play important roles in mediating downstream cellular signals, and are potential targets for novel therapies. Using a functional proteomic approach, we determined the identity of proteins that interact with the ALK tyrosine kinase by co-immunoprecipitation with anti-ALK antibody, followed by electrospray ionization and tandem mass spectrometry (MS/MS). A total of 46 proteins were identified as unique to the ALK immunocomplex using monoclonal and polyclonal antibodies, while 11 proteins were identified in the NPM immunocomplex. Previously reported proteins in the ALK signal pathway were identified including PI3-K, Jak2, Jak3, Stat3, Grb2, IRS, and PLCgamma1. More importantly, many proteins previously not recognized to be associated with NPM-ALK, but with potential NPM-ALK interacting protein domains, were identified. These include adaptor molecules (SOCS, Rho-GTPase activating protein, RAB35), kinases (MEK kinase 1 and 4, PKC, MLCK, cyclin G-associated kinase, EphA1, JNK kinase, MAP kinase 1), phosphatases (meprin, PTPK, protein phosphatase 2 subunit), and heat shock proteins (Hsp60 precursor). Proteins identified by MS were confirmed by Western blotting and reciprocal immunoprecipitation. This study demonstrates the utility of antibody immunoprecipitation and subsequent peptide identification by tandem mass spectrometry for the elucidation of ALK-binding proteins, and its potential signal transduction pathways.

High-throughput analysis of protein/peptide complexes by immunoprecipitation and automated LC-MS/MS

BACKGROUND

The identification of interacting proteins within protein complexes is key to understanding the transduction and regulation of cell signaling pathways, and is also a useful tool for identifying novel disease markers. Immunoprecipitation of protein complexes followed by liquid chromatography-tandem mass spectrometry (LC-MS/MS) has been used to identify targets that bind to a protein of interest. Here, we report a high-throughput nanoflow LC-MS/MS method for analysis of protein/peptide complexes.

APPROACH

To overcome the large dwell volume from connecting lines and the injector when using an autosampler in microcapillary LC-MS/MS, we employed a valve-controlled variable flow method with a peptide trap. This method enables fast trapping of peptides from samples injected by an autosampler within minutes followed by a split-controlled nanoflow of acetonitrile gradient through a microcapillary C18 column to separate and elute peptides into the ion-trap MS/MS. Over 40 protein/peptide samples at femtomole levels could be analyzed continuously using the same in-house packed microcapillary C18 column. The p38 mitogen-activated protein kinase (p38 MAP kinase) is a signaling protein that is involved in transduction of extracellular signals including oxidative stress, growth factors, and cytokines, with diverse cellular consequences such as cell proliferation, differentiation, or apoptosis. Immunocomplex of monoclonal anti-p38 beta antibodies from 2 mg total lysates from the OCI LY-1 lymphoma cell line was resolved in 1D-PAGE gel followed by silver staining of the gel. Protein bands were excised and digested with trypsin. Peptides were extracted and analyzed using the automated LC-MS/MS method.

RESULTS AND CONCLUSIONS

From 37 excised protein bands in the 1D-PAGE gel, we identified more than 50 proteins, including cytoskeletal proteins, ribosomal proteins, transcription factors, and KIAA potential signaling proteins. These proteins are the potential targets that may interact directly or indirectly with the p38 MAP kinase proteins. Our studies demonstrate the utility of automated nanoflow LC-MS/MS for sensitive and high-throughput analysis of protein/peptide complexes. Utilization of methods based on this principle would greatly facilitate peptide interaction mapping and other proteomic studies requiring high-throughput pre-analytical sample preparation.

Rapid analysis of metanephrine and normetanephrine in urine by gas chromatography-mass spectrometry

BACKGROUND

Widely used HPLC methods for quantification of metanephrine and normetanephrine in urine often have long analysis times and are frequently plagued by drug interferences. We describe a gas chromatography-mass spectrometry method designed to overcome these limitations.

METHODS

Metanephrine and normetanephrine conjugates were converted to unconjugated metanephrine and normetanephrine by acid hydrolysis. To avoid the rapid decomposition of the deuterated internal standards (metanephrine-d(3) and normetanephrine-d(3)) under hydrolysis conditions, the internal standards were added after hydrolysis. Solid-phase extraction was used to isolate the hydrolyzed metanephrines from urine. Samples were concentrated by evaporation, then derivatized simultaneously with N-methyl-N-(trimethylsilyl)trifluoroacetamide and N-methyl-bis-heptafluoro-butryamide at room temperature.

RESULTS

The assay was linear from 25 to 7000 microg/L. The intraassay CVs were < 5 % and the interassay CVs < 12%. Comparison with a routine HPLC method (n = 192) by Deming regression yielded a slope of 1.00 +/- 0.02 microg/L, an intercept of -5.8 +/- 7.8 micro/L, and S(y/x) = 50.6 microg/L for metanephrine and a slope of 0.94 +/- 0.03, intercept of 19 +/- 11 microg/L, and S(y/x) = 60 microg/L for normetanephrine. The correlation coefficients (r) were calculated after log transformation of the data and gave r = 0.97 for metanephrine and r = 0.97 for normetanephrine. Interference from common medications or drug metabolites was seen in <1% of samples. The time between sequential injections was < 7 min.

CONCLUSIONS

This new gas chromatography-mass spectrometry assay for total fractionated metanephrines is rapid, compares well with a standard HPLC assay, and avoids most drug interferences that commonly affect HPLC assays for urine metanephrines.

Rapid Analysis of Metanephrine and Normetanephrine in Urine by Gas Chromatography-Mass Spectrometry

Background: Widely used HPLC methods for quantification of metanephrine and normetanephrine in urine often have long analysis times and are frequently plagued by drug interferences. We describe a gas chromatography-mass spectrometry method designed to overcome these limitations.

Methods: Metanephrine and normetanephrine conjugates were converted to unconjugated metanephrine and normetanephrine by acid hydrolysis. To avoid the rapid decomposition of the deuterated internal standards (metanephrine-d3 and normetanephrine-d3) under hydrolysis conditions, the internal standards were added after hydrolysis. Solid-phase extraction was used to isolate the hydrolyzed metanephrines from urine. Samples were concentrated by evaporation, then derivatized simultaneously with N-methyl-N-(trimethylsilyl)trifluoroacetamide and N-methyl-bis-heptafluoro-butryamide at room temperature.

Results: The assay was linear from 25 to 7000 μg/L. The intraassay CVs were &lt;5% and the interassay CVs &lt;12%. Comparison with a routine HPLC method (n = 192) by Deming regression yielded a slope of 1.00 ± 0.02 μg/L, an intercept of −5.8 ± 7.8 μg/L, and Sy|x = 50.6 μg/L for metanephrine and a slope of 0.94 ± 0.03, intercept of 19 ± 11 μg/L, and Sy|x = 60 μg/L for normetanephrine. The correlation coefficients (r) were calculated after log transformation of the data and gave r = 0.97 for metanephrine and r = 0.97 for normetanephrine. Interference from common medications or drug metabolites was seen in &lt;1% of samples. The time between sequential injections was &lt;7 min.

Conclusions: This new gas chromatography-mass spectrometry assay for total fractionated metanephrines is rapid, compares well with a standard HPLC assay, and avoids most drug interferences that commonly affect HPLC assays for urine metanephrines.

Solid-phase extraction of 11-nor-delta9-tetrahydrocannabinol-9-carboxylic acid from urine drug-testing specimens with the cerex polycrom-THC column

Confirmation of drugs of abuse by gas chromatography-mass spectrometry (GC-MS) is the most time-consuming process used by drug-testing laboratories. Cost effectiveness and competitive turnaround times for testing results demand fast, efficient, and reliable extraction methods. We applied the Cerex Polycrom-THC solid-phase extraction (SPE) column to the extraction of 11-nor-delta9-tetrahydrocannabinol carboxylic acid (9-THCA) from urine. This column uses an anion exchange divinyl-benzene copolymer, which requires no pH adjustment after hydrolysis of the THC-glucuronide with base, as is necessary with many other SPE columns. With urine, no preconditioning of the SPE column was necessary. After extraction, trimethylsilation derivatization was performed with N-methyl-N-(trimethylsilyl) trifluoroacetamide. This further improved efficiency because no heated incubation was required. A method correlation to an existing liquid-liquid extraction was performed. Analysis was on a Finnigan Voyager GC-MS with a 2.5-min run time per injection. Using 9-THCA-d9 deuterated internal standard, the assay was linear from 2 to 2000 ng/mL. Total precision at the 15-ng/mL cutoff concentration was 4.4%. The Cerex column was also evaluated for interference using two common drug-test adulterants, Klear and Urine Luck. Considerably less interference was observed when compared to an existing liquid-liquid method.

Comparison of four derivatizing reagents for 6-acetylmorphine GC-MS analysis

The propionyl, trimethylsilyl, trifluroacetyl, and heptafluoroacyl derivatives of 6-acetylmorphine (6-AM) were evaluated with respect to optimal method performance, derivative stability, and methods characterization for use in gas chromatographic-mass spectrometric (GC-MS) analysis with electron ionization mode and selected ion monitoring. The most common potential interferences and compatibility with other derivatives when used on the same GC-MS were determined for the derivatizing reagents. The propionyl, trimethylsilyl, and trifluroacetyl derivatives produced adequate stability, accuracy, and precision for the method. The 6-AM derivatization with commercially available propionic anhydride generated a relatively small amount of 6-AM-propionyl derivative from the free morphine present in a specimen. The trimethylsilyl derivative obtained by the reaction with MSTFA did not require incubation, was the easiest to prepare, and had the highest potential for use on an automated sample-preparation device. An important advantage of derivatization with MSTFA is elimination of the possibility of heroin decomposition to 6-AM that is due to incubation at elevated temperature.

Nitrite adulteration of workplace urine drug-testing specimens. I. Sources and associated concentrations of nitrite in urine and distinction between natural sources and adulteration

The active ingredient in the commercial workplace urine drug-testing adulterant, Klear, was previously determined to be nitrite ion. Nitrite adulteration compromises the confirmation of some drugs, notably the marijuana metabolite. A previously reported bisulfite step overcomes some nitrite adulteration, but it cannot do so in every case, which leaves the laboratory to report the specimen as not suitable for testing. Unlike many other adulterants, nitrite is found in normal urine at low concentrations. In order to defend a report of nitrite adulteration, it is necessary to provide evidence that the amount of nitrite in a workplace urine specimen could not arise by normal means. The objectives of this study were to identify all sources of nitrite in urine and the range of concentrations associated with these sources and to determine if nitrite adulteration can be supported based upon a quantitative result. The scientific literature was reviewed for internal and external sources of nitrite and their concentration ranges and are reported. The following specimens were obtained and nitrite concentrations measured by a spectrophotometric method: clinical specimens nitrite positive by test strip (< 15 micrograms/mL); specimens culture positive for nitrate-reducing microorganisms (< 36 micrograms/mL); specimens from patients on medications that may metabolize to nitrite (< 6 micrograms/mL); and drug-test specimens, both negative (< 130 micrograms/mL) and others that appeared to be adulterated with nitrite (range 1910-12,200 micrograms/mL, mean 5910). The literature and the nitrite measurements of this study indicate a substantial difference between concentrations from natural sources compared with adulteration. A quantitative measurement of nitrite by a well-structured assay can provide scientifically valid and forensically defensible proof of adulteration with a nitrite-containing substance.