Addressing accuracy and precision issues in iTRAQ quantitation

iTRAQ (isobaric tags for relative or absolute quantitation) is a mass spectrometry technology that allows quantitative comparison of protein abundance by measuring peak intensities of reporter ions released from iTRAQ-tagged peptides by fragmentation during MS/MS. However, current data analysis techniques for iTRAQ struggle to report reliable relative protein abundance estimates and suffer with problems of precision and accuracy. The precision of the data is affected by variance heterogeneity: low signal data have higher relative variability; however, low abundance peptides dominate data sets. Accuracy is compromised as ratios are compressed toward 1, leading to underestimation of the ratio. This study investigated both issues and proposed a methodology that combines the peptide measurements to give a robust protein estimate even when the data for the protein are sparse or at low intensity. Our data indicated that ratio compression arises from contamination during precursor ion selection, which occurs at a consistent proportion within an experiment and thus results in a linear relationship between expected and observed ratios. We proposed that a correction factor can be calculated from spiked proteins at known ratios. Then we demonstrated that variance heterogeneity is present in iTRAQ data sets irrespective of the analytical packages, LC-MS/MS instrumentation, and iTRAQ labeling kit (4-plex or 8-plex) used. We proposed using an additive-multiplicative error model for peak intensities in MS/MS quantitation and demonstrated that a variance-stabilizing normalization is able to address the error structure and stabilize the variance across the entire intensity range. The resulting uniform variance structure simplifies the downstream analysis. Heterogeneity of variance consistent with an additive-multiplicative model has been reported in other MS-based quantitation including fields outside of proteomics; consequently the variance-stabilizing normalization methodology has the potential to increase the capabilities of MS in quantitation across diverse areas of biology and chemistry.

Correlation between extent of osteolytic damage and metastatic burden of human breast cancer metastasis in nude mice : real-time PCR quantitation

Orthotopic or intracardiac injection of human breast cancer cell lines into immunocompromised mice allows study of the molecular basis of breast cancer metastasis. We have established a quantitative real-time PCR approach to analyze metastatic spread of human breast cancer cells inoculated into nude mice via these routes. We employed MDA-MB-231 human breast cancer cells genetically tagged with a bacterial β-galactosidase (Lac-Z) retroviral vector, enabling their detection by TaqMan® real-time PCR. PCR detection was linear, specific, more sensitive than conventional PCR, and could be used to directly quantitate metastatic burden in bone and soft organs. Attesting to the sensitivity and specificity of the PCR detection strategy, as few as several hundred metastatic MDA-MB-231 cells were detectable in 100 μm segments of paraffin-embedded lung tissue, and only in samples adjacent to sections that scored positive by histological detection. Moreover, the measured real-time PCR metastatic burden in the bone environment (mouse hind-limbs, n = 48) displayed a high correlation to the degree of osteolytic damage observed by high resolution X-ray analysis (r2 = 0.972). Such a direct linear relationship to tumor burden and bone damage substantiates the so-called 'vicious cycle' hypothesis in which metastatic tumor cells promote the release of factors from the bone which continue to stimulate the tumor cells. The technique provides a useful tool for molecular and cellular analysis of human breast cancer metastasis to bone and soft organs, can easily be extended to other cell/marker/organ systems, and should also find application in preclinical assessment of anti-metastatic modalities.

Quantitation of gammaH2AX foci in tissue samples

DNA double-strand breaks (DSBs) are particularly lethal and genotoxic lesions, that can arise either by endogenous (physiological or pathological) processes or by exogenous factors, particularly ionizing radiation and radiomimetic compounds. Phosphorylation of the H2A histone variant, H2AX, at the serine-139 residue, in the highly conserved C-terminal SQEY motif, forming γH2AX, is an early response to DNA double-strand breaks1. This phosphorylation event is mediated by the phosphatidyl-inosito 3-kinase (PI3K) family of proteins, ataxia telangiectasia mutated (ATM), DNA-protein kinase catalytic subunit and ATM and RAD3-related (ATR)2. Overall, DSB induction results in the formation of discrete nuclear γH2AX foci which can be easily detected and quantitated by immunofluorescence microscopy2. Given the unique specificity and sensitivity of this marker, analysis of γH2AX foci has led to a wide range of applications in biomedical research, particularly in radiation biology and nuclear medicine. The quantitation of γH2AX foci has been most widely investigated in cell culture systems in the context of ionizing radiation-induced DSBs. Apart from cellular radiosensitivity, immunofluorescence based assays have also been used to evaluate the efficacy of radiation-modifying compounds. In addition, γH2AX has been used as a molecular marker to examine the efficacy of various DSB-inducing compounds and is recently being heralded as important marker of ageing and disease, particularly cancer3. Further, immunofluorescence-based methods have been adapted to suit detection and quantitation of γH2AX foci ex vivo and in vivo4,5. Here, we demonstrate a typical immunofluorescence method for detection and quantitation of γH2AX foci in mouse tissues.

An immunodiagnostic assay for quantitation of specific ige to the major pollen allergen component, pas n 1, of the subtropical bahia grass

Background Pollens of the Panicoideae subfamily of grasses including Bahia (Paspalum notatum) are important allergen sources in subtropical regions of the world. An assay for specific IgE to the major molecular allergenic component, Pas n 1, of Bahia grass pollen (BaGP) would have immunodiagnostic utility for patients with pollen allergy in these regions. Methods Biotinylated Pas n 1 purified from BaGP was coated onto streptavidin ImmunoCAPs. Subjects were assessed by clinical history of allergic rhinitis and skin prick test (SPT) to aeroallergens. Serum total, BaGP-specific and Pas n 1-specific IgE were measured. Results: Pas n 1 IgE concentrations were highly correlated with BaGP SPT (r = 0.795, p < 0.0001) and BaGP IgE (r = 0.915, p < 0.0001). At 0.23 kU/l Pas n 1 IgE, the diagnostic sensitivity (92.4%) and specificity (93.1%) for the detection of BaGP allergy was high (area under receiver operator curve 0.960, p < 0.0001). The median concentrations of Pas n 1 IgE in non-Atopic subjects (0.01 kU/l, n = 67) and those with other allergies (0.02 kU/l, n = 59) showed no inter-group difference, whilst grass pollen-Allergic patients with allergic rhinitis showed elevated Pas n 1 IgE (6.71 kU/l, n = 182, p < 0.0001). The inter-Assay coefficient of variation for the BaGP-Allergic serum pool was 6.92%. Conclusions Pas n 1 IgE appears to account for most of the BaGP-specific IgE. This molecular component immunoassay for Pas n 1 IgE has potential utility to improve the sensitivity and accuracy of diagnosis of BaGP allergy for patients in subtropical regions.

A longitudinal assessment of chronic wound fluid to detect biochemical indicators of healing

Chronic venous leg ulcers are a detrimental health issue plaguing our society, resulting in long term pain, immobility and decreased quality of life for a large proportion of sufferers. The frequency of these chronic wounds has led current research to focus on the wound environment to provide important information regarding the prolonged, fluctuated or static healing patterns of these wounds. Disruption to the normal wound healing process results in release of multiple factors in the wound environment that could correlate to wound chronicity. These biochemical factors can often be detected through non-invasively sampling chronic wound fluid (CWF) from the site of injury. Of note, whilst there are numerous studies comparing acute and chronic wound fluids, there have not been any reports in the literature employing a longitudinal study in order to track biochemical changes in wound fluid as patients transition from a non-healing to healed state. Initially the objective of this study was to identify biochemical changes in CWF associated with wound healing using a proteomic approach. The proteomic approach incorporated a multi-dimensional liquid chromatography fractionation technique coupled with mass spectrometry (MS) to enable identification of proteins present in lower concentrations in CWF. Not surprisingly, many of the proteins identified in wound fluid were acute phase proteins normally expressed during the inflammatory phase of healing. However, the number of proteins positively identified by MS was quite low. This was attributed to the diverse range in concentration of protein species in CWF making it challenging to detect the diagnostically relevant low molecular weight proteins. In view of this, SELDI-TOF MS was also explored as a means to target low molecular weight proteins in sequential patient CWF samples during the course of healing. Unfortunately, the results generated did not yield any peaks of interest that were altered as wounds transitioned to a healed state. During the course of proteomic assessment of CWF, it became evident that a fraction of non-proteinaceous compounds strongly absorbed at 280 nm. Subsequent analyses confirmed that most of these compounds were in fact part of the purine catabolic pathway, possessing distinctive aromatic rings and which results in high absorbance at 254 nm. The accumulation of these purinogenic compounds in CWF suggests that the wound bed is poorly oxygenated resulting in a switch to anaerobic metabolism and consequently ATP breakdown. In addition, the presence of the terminal purine catabolite, uric acid (UA), indicates that the enzyme xanthine oxidoreductase (XOR) catalyses the reaction of hypoxanthine to xanthine and finally to UA. More importantly, the studies provide evidence for the first time of the exogenous presence of XOR in CWF. XOR is the only enzyme in humans capable of catalysing the production of UA in conjunction with a burst of the highly reactive superoxide radical and other oxidants like H2O2. Excessive release of these free radicals in the wound environment can cause cellular damage disrupting the normal wound healing process. In view of this, a sensitive and specific assay was established for monitoring low concentrations of these catabolites in CWF. This procedure involved combining high performance liquid chromatography (HPLC) with tandem mass spectrometry and multiple reaction monitoring (MRM). This application was selective, using specific MRM transitions and HPLC separations for each analyte, making it ideal for the detection and quantitation of purine catabolites in CWF. The results demonstrated that elevated levels of UA were detected in wound fluid obtained from patients with clinically worse ulcers. This suggests that XOR is active in the wound site generating significant amounts of reactive oxygen species (ROS). In addition, analysis of the amount of purine precursors in wound fluid revealed elevated levels of purine precursors in wound fluid from patients with less severe ulcers. Taken together, the results generated in this thesis suggest that monitoring changes of purine catabolites in CWF is likely to provide valuable information regarding the healing patterns of chronic venous leg ulcers. XOR catalysis of purine precursors not only provides a method for monitoring the onset, prognosis and progress of chronic venous leg ulcers, but also provides a potential therapeutic target by inhibiting XOR, thus blocking UA and ROS production. Targeting a combination of these purinogenic compounds and XOR could lead to the development of novel point of care diagnostic tests. Therefore, further investigation of these processes during wound healing will be worthwhile and may assist in elucidating the pathogenesis of this disease state, which in turn may lead to the development of new diagnostics and therapies that target these processes.

Prostate cancer biomarkers and the emerging proteomic techniques used in biomarker research

Prostate cancer (CaP) is the second leading cause of cancer-related deaths in North American males and the most common newly diagnosed cancer in men world wide. Biomarkers are widely used for both early detection and prognostic tests for cancer. The current, commonly used biomarker for CaP is serum prostate specific antigen (PSA). However, the specificity of this biomarker is low as its serum level is not only increased in CaP but also in various other diseases, with age and even body mass index. Human body fluids provide an excellent resource for the discovery of biomarkers, with the advantage over tissue/biopsy samples of their ease of access, due to the less invasive nature of collection. However, their analysis presents challenges in terms of variability and validation. Blood and urine are two human body fluids commonly used for CaP research, but their proteomic analyses are limited both by the large dynamic range of protein abundance making detection of low abundance proteins difficult and in the case of urine, by the high salt concentration. To overcome these challenges, different techniques for removal of high abundance proteins and enrichment of low abundance proteins are used. Their applications and limitations are discussed in this review. A number of innovative proteomic techniques have improved detection of biomarkers. They include two dimensional differential gel electrophoresis (2D-DIGE), quantitative mass spectrometry (MS) and functional proteomic studies, i.e., investigating the association of post translational modifications (PTMs) such as phosphorylation, glycosylation and protein degradation. The recent development of quantitative MS techniques such as stable isotope labeling with amino acids in cell culture (SILAC), isobaric tags for relative and absolute quantitation (iTRAQ) and multiple reaction monitoring (MRM) have allowed proteomic researchers to quantitatively compare data from different samples. 2D-DIGE has greatly improved the statistical power of classical 2D gel analysis by introducing an internal control. This chapter aims to review novel CaP biomarkers as well as to discuss current trends in biomarker research from two angles: the source of biomarkers (particularly human body fluids such as blood and urine), and emerging proteomic approaches for biomarker research.

Pharmacokinetic studies of a Chinese triple herbal drug formula

Yin Chen Hao Tang preparation (YCHTP) is a classic traditional Chinese medicine formula, which is commonly used for clinical treatment of hepatological diseases. In this study, a rapid and validated high-performance liquid chromatography (HPLC) method was developed to simultaneously identify 6,7-dimethylesculetin and geniposide in rat plasma. This assay was performed on a Dikma Diamonsil RP(18) column (200 mmx4.6 mm, 5 mum) with acetonitrile-methanol-water (0.1% formic acid) as the mobile phase, showing acceptable linearity, intra- and inter-day precision and accuracy (R.S.D.=5%), and absolute recovery for two analytes (74%); the limits of quantitation were 0.4 and 1.12 mug/ml, and the limits of detection were 0.06 and 0.09 mug/ml for two analytes. The developed method was successfully applied to study the effect of formula compatibility on the pharmacokinetics of 6,7-dimethylesculetin and geniposide in YCHTP when orally administrating an effective human daily dose of YCHTP to rats. We surmise that formula compatibility can significantly influence the pharmacokinetics of YCHTP, and we have elucidated and validated the compatible administration of YCHTP.

Simultaneous determination of 6,7-dimethylesculetin and geniposide in rat plasma and its application to pharmacokinetic studies of Yin Chen Hao Tang preparation

A method for the rapid and simultaneous determination of 6,7-dimethylesculetin (CAS 120-08-1) and geniposide (CAS 24512-63-8) in rat plasma has been developed, using validated high performance liquid chromatography (HPLC) with solid phase extraction (SPE). The HPLC analysis was performed on a commercially available column (200 mm x 4.6 mm, 5 microm) with acetonitrile-methanol-0.1% aqueous formic acid as mobile phase and the UV detection at 343 nm and 238 nm for 6,7-dimethylesculetin and geniposide, respectively. The calibration curves for 6,7-dimethylesculetin and geniposide were linear over the range 0.4-25.6 microg/mL and 1.12-71.68 microg/mL, respectively. The lower limits of quantitation were 0.40 microg/ mL and 1.12 microg/mL, and the lower limits of detection were 0.06 microg/mL and 0.09 microg/ mL, respectively. The intra-day and inter-day precision for 6,7-dimethylesculetin and geniposide were < 5%, whereas the absolute recovery percentages were > 74%. A successful application of the developed HPLC analysis was demonstrated for the pharmacokinetic study of a Traditional Chinese Medicine formula of Yin Chen Hao Tang preparation.

Detection of mRNA levels for the estrogen alpha, estrogen beta and androgen nuclear receptor genes in archival breast cancer tissue

Previous studies in our laboratory have shown association of nuclear receptor expression and histological breast cancer grade. To further investigate these findings, it was the objective of this study to determine if expression levels of the estrogen alpha, estrogen beta and androgen nuclear receptor genes varied in different breast cancer grades. RNA extracted from paraffin embedded archival breast tumour tissue was converted into cDNA and cDNA underwent PCR to enable quantitation of mRNA expression. Expression data was normalised against the 18S ribosomal gene multiplex and analysed using ANOVA. Analysis indicated a significant alteration of expression for the androgen receptor in different cancer grades (P=0.014), as well as in tissues that no longer possess estrogen receptor alpha proteins (P=0.025). However, expression of estrogen receptors alpha and beta did not vary significantly with cancer grade (P=0.057 and 0.622, respectively). Also, the expression of estrogen receptor alpha or beta did not change, regardless of the presence of estrogen receptor alpha protein in the tissue (P=0.794 and 0.716, respectively). Post-hoc tests indicate that the expression of the androgen receptor is increased in estrogen receptor negative tissue as well as in grade 2 and grade 3 tumours, compared to control tissue. This increased expression in late stage breast tumours may have implications to the treatment of breast tumours, particularly those lacking expression of other nuclear receptor genes.

A molecular genetic approach for forensic animal species identification

This study investigated potential markers within chromosomal, mitochondrial DNA (mtDNA) and ribosomal RNA (rRNA) with the aim of developing a DNA based method to allow differentiation between animal species. Such discrimination tests may have important applications in the forensic science, agriculture, quarantine and customs fields. DNA samples from five different animal individuals within the same species for 10 species of animal (including human) were analysed. DNA extraction and quantitation followed by PCR amplification and GeneScan visualisation formed the basis of the experimental analysis. Five gene markers from three different types of genes were investigated. These included genomic markers for the β-actin and TP53 tumor suppressor gene. Mitochondrial DNA markers, designed by Bataille et al. [Forensic Sci. Int. 99 (1999) 165], examined the Cytochrome b gene and Hypervariable Displacement Loop (D-Loop) region. Finally, a ribosomal RNA marker for the 28S rRNA gene optimised by Naito et al. [J. Forensic Sci. 37 (1992) 396] was used as a possible marker for speciation. Results showed a difference of only several base pairs between all species for the β-actin and 28S markers, with the exception of Sus scrofa (pig) β-actin fragment length, which produced a significantly smaller fragment. Multiplexing of Cytochrome b and D-Loop markers gave limited species information, although positive discrimination of human DNA was evident. The most specific and discriminatory results were shown using the TP53 gene since this marker produced greatest fragment size differences between animal species studied. Sample differentiation for all species was possible following TP53 amplification, suggesting that this gene could be used as a potential animal species identifier.

Expression of glucocorticoid and progesterone nuclear receptor genes in archival breast cancer tissue.

BACKGROUND: Previous studies in our laboratory have shown associations of specific nuclear receptor gene variants with sporadic breast cancer. In order to investigate these findings further, we conducted the present study to determine whether expression levels of the progesterone and glucocorticoid nuclear receptor genes vary in different breast cancer grades. METHODS: RNA was extracted from paraffin-embedded archival breast tumour tissue and converted into cDNA. Sample cDNA underwent PCR using labelled primers to enable quantitation of mRNA expression. Expression data were normalized against the 18S ribosomal gene multiplex and analyzed using analysis of variance. RESULTS: Analysis of variance indicated a variable level of expression of both genes with regard to breast cancer grade (P = 0.00033 for glucocorticoid receptor and P = 0.023 for progesterone receptor). CONCLUSION: Statistical analysis indicated that expression of the progesterone nuclear receptor is elevated in late grade breast cancer tissue.

Quantification of gammaH2AX foci in response to ionising radiation

DNA double-strand breaks (DSBs), which are induced by either endogenous metabolic processes or by exogenous sources, are one of the most critical DNA lesions with respect to survival and preservation of genomic integrity. An early response to the induction of DSBs is phosphorylation of the H2A histone variant, H2AX, at the serine-139 residue, in the highly conserved C-terminal SQEY motif, forming gammaH2AX(1). Following induction of DSBs, H2AX is rapidly phosphorylated by the phosphatidyl-inosito 3-kinase (PIKK) family of proteins, ataxia telangiectasia mutated (ATM), DNA-protein kinase catalytic subunit and ATM and RAD3-related (ATR)(2). Typically, only a few base-pairs (bp) are implicated in a DSB, however, there is significant signal amplification, given the importance of chromatin modifications in DNA damage signalling and repair. Phosphorylation of H2AX mediated predominantly by ATM spreads to adjacent areas of chromatin, affecting approximately 0.03% of total cellular H2AX per DSB(2,3). This corresponds to phosphorylation of approximately 2000 H2AX molecules spanning approximately 2 Mbp regions of chromatin surrounding the site of the DSB and results in the formation of discrete gammaH2AX foci which can be easily visualized and quantitated by immunofluorescence microscopy(2). The loss of gammaH2AX at DSB reflects repair, however, there is some controversy as to what defines complete repair of DSBs; it has been proposed that rejoining of both strands of DNA is adequate however, it has also been suggested that re-instatement of the original chromatin state of compaction is necessary(4-8). The disappearence of gammaH2AX involves at least in part, dephosphorylation by phosphatases, phosphatase 2A and phosphatase 4C(5,6). Further, removal of gammaH2AX by redistribution involving histone exchange with H2A.Z has been implicated(7,8). Importantly, the quantitative analysis of gammaH2AX foci has led to a wide range of applications in medical and nuclear research. Here, we demonstrate the most commonly used immunofluorescence method for evaluation of initial DNA damage by detection and quantitation of gammaH2AX foci in gamma-irradiated adherent human keratinocytes(9)

Evaluation of the spatial distribution of γH2AX following ionizing radiation

An early molecular response to DNA double-strand breaks (DSBs) is phosphorylation of the Ser-139 residue within the terminal SQEY motif of the histone H2AX1,2. This phosphorylation of H2AX is mediated by the phosphatidyl-inosito 3-kinase (PI3K) family of proteins, ataxia telangiectasia mutated (ATM), DNA-protein kinase catalytic subunit and ATM and RAD3-related (ATR)3. The phosphorylated form of H2AX, referred to as γH2AX, spreads to adjacent regions of chromatin from the site of the DSB, forming discrete foci, which are easily visualized by immunofluorecence microscopy3. Analysis and quantitation of γH2AX foci has been widely used to evaluate DSB formation and repair, particularly in response to ionizing radiation and for evaluating the efficacy of various radiation modifying compounds and cytotoxic compounds Given the exquisite specificity and sensitivity of this de novo marker of DSBs, it has provided new insights into the processes of DNA damage and repair in the context of chromatin. For example, in radiation biology the central paradigm is that the nuclear DNA is the critical target with respect to radiation sensitivity. Indeed, the general consensus in the field has largely been to view chromatin as a homogeneous template for DNA damage and repair. However, with the use of γH2AX as molecular marker of DSBs, a disparity in γ-irradiation-induced γH2AX foci formation in euchromatin and heterochromatin has been observed5-7. Recently, we used a panel of antibodies to either mono-, di- or tri- methylated histone H3 at lysine 9 (H3K9me1, H3K9me2, H3K9me3) which are epigenetic imprints of constitutive heterochromatin and transcriptional silencing and lysine 4 (H3K4me1, H3K4me2, H3K4me3), which are tightly correlated actively transcribing euchromatic regions, to investigate the spatial distribution of γH2AX following ionizing radiation8. In accordance with the prevailing ideas regarding chromatin biology, our findings indicated a close correlation between γH2AX formation and active transcription9. Here we demonstrate our immunofluorescence method for detection and quantitation of γH2AX foci in non-adherent cells, with a particular focus on co-localization with other epigenetic markers, image analysis and 3Dmodeling.

Droplet digital PCR as a useful tool for the quantitative detection of Enterovirus 71

Hand, foot and mouth disease (HFMD) is a contagious viral disease that frequently affects infants and children and present with blisters and flu-like symptoms. This disease is caused by a group of enteroviruses such as enterovirus 71 (EV71) and coxsackievirus A16 (CA16). However, unlike other HFMD causing enteroviruses, EV71 have also been shown to be associated with more severe clinical manifestation such as aseptic meningitis, brainstem and cerebellar encephalitis which may lead to cardiopulmonary failure and death. Clinically, HFMD caused by EV71 is indistinguishable from other HFMD causing enteroviruses such as CA16. Molecular diagnosis methods such as the use of real-time PCR has been used commonly for the identification of EV71. In this study, two platforms namely the real-time PCR and the droplet digital PCR were compared for the detection quantitation of known EV71 viral copy number. The results reveal accurate and consistent results between the two platforms. In summary, the droplet digital PCR was demonstrated to be a promising technology for the identification and quantitation of EV71 viral copy number.