Explainable machine learning model for predicting skeletal muscle loss during surgery and adjuvant chemotherapy in ovarian cancer.

BACKGROUND: Skeletal muscle loss during treatment is associated with poor survival outcomes in patients with ovarian cancer. Although changes in muscle mass can be assessed on computed tomography (CT) scans, this labour-intensive process can impair its utility in clinical practice. This study aimed to develop a machine learning (ML) model to predict muscle loss based on clinical data and to interpret the ML model by applying SHapley Additive exPlanations (SHAP) method. METHODS: This study included the data of 617 patients with ovarian cancer who underwent primary debulking surgery and platinum-based chemotherapy at a tertiary centre between 2010 and 2019. The cohort data were split into training and test sets based on the treatment time. External validation was performed using 140 patients from a different tertiary centre. The skeletal muscle index (SMI) was measured from pre- and post-treatment CT scans, and a decrease in SMI ≥ 5% was defined as muscle loss. We evaluated five ML models to predict muscle loss, and their performance was determined using the area under the receiver operating characteristic curve (AUC) and F1 score. The features for analysis included demographic and disease-specific characteristics and relative changes in body mass index (BMI), albumin, neutrophil-to-lymphocyte ratio (NLR), and platelet-to-lymphocyte ratio (PLR). The SHAP method was applied to determine the importance of the features and interpret the ML models. RESULTS: The median (inter-quartile range) age of the cohort was 52 (46-59) years. After treatment, 204 patients (33.1%) experienced muscle loss in the training and test datasets, while 44 (31.4%) patients experienced muscle loss in the external validation dataset. Among the five evaluated ML models, the random forest model achieved the highest AUC (0.856, 95% confidence interval: 0.854-0.859) and F1 score (0.726, 95% confidence interval: 0.722-0.730). In the external validation, the random forest model outperformed all ML models with an AUC of 0.874 and an F1 score of 0.741. The results of the SHAP method showed that the albumin change, BMI change, malignant ascites, NLR change, and PLR change were the most important factors in muscle loss. At the patient level, SHAP force plots demonstrated insightful interpretation of our random forest model to predict muscle loss. CONCLUSIONS: Explainable ML model was developed using clinical data to identify patients experiencing muscle loss after treatment and provide information of feature contribution. Using the SHAP method, clinicians may better understand the contributors to muscle loss and target interventions to counteract muscle loss.

SheddomeDB 2023: A Revision of an Ectodomain Shedding Database Based on a Comprehensive Literature Review and Online Resources.

Ectodomain shedding of membrane proteins is a proteolytic event involved in several biological phenomena, including inflammation, development, diseases, and cancer progression. Though ectodomain shedding is a post-translational modification that plays an important role in cellular regulation, this biological phenomenon is seriously underannotated in public protein databases. Given the importance of the shedding events, we conducted a comprehensive literature review for membrane protein shedding and constructed the database, SheddomeDB in 2017. In response to user feedback, novel shedding findings, more associated biomedical events, and the advance in web technology, we revised SheddomeDB to a new version, SheddomeDB 2023. The revised SheddomeDB 2023 includes 481 protein entries across seven species; all the content was manually verified and curated. The content of SheddomeDB 2023 mainly came from a comprehensive literature survey by our newly developed semiautomated screening tool. We also integrated verified and updated cleavage and secretome information from other databases into the revision. In addition, SheddomeDB 2023 features a graphical presentation of cleavage information and a user-friendly interface for searching and browsing entries in the database. This revised comprehensive database of ectodomain shedding is expected to benefit biomedical researchers across different disciplines.

Association of prognostic nutritional index with muscle loss and survival in patients with ovarian cancer treated with primary debulking surgery and chemotherapy.

PURPOSE: Sarcopenia is prevalent in ovarian cancer and contributes to poor survival. This study is aimed at investigating the association of prognostic nutritional index (PNI) with muscle loss and survival outcomes in patients with ovarian cancer. METHODS: This retrospective study analyzed 650 patients with ovarian cancer treated with primary debulking surgery and adjuvant platinum-based chemotherapy at a tertiary center from 2010 to 2019. PNI-low was defined as a pretreatment PNI of < 47.2. Skeletal muscle index (SMI) was measured on pre- and posttreatment computed tomography (CT) at L3. The cut-off for the SMI loss associated with all-cause mortality was calculated using maximally selected rank statistics. RESULTS: The median follow-up was 4.2 years, and 226 deaths (34.8%) were observed. With a median duration of 176 days (interquartile range: 166-187) between CT scans, patients experienced an average decrease in SMI of 1.7% (P < 0.001). The cut-off for SMI loss as a predictor of mortality was - 4.2%. PNI-low was independently associated with SMI loss (odds ratio: 1.97, P = 0.001). On multivariable analysis of all-cause mortality, PNI-low and SMI loss were independently associated with all-cause mortality (hazard ratio: 1.43, P = 0.017; hazard ratio: 2.27, P < 0.001, respectively). Patients with both SMI loss and PNI-low (vs. neither) had triple the risk of all-cause mortality (hazard ratio: 3.10, P < 0.001). CONCLUSIONS: PNI is a predictor of muscle loss during treatment for ovarian cancer. PNI and muscle loss are additively associated with poor survival. PNI can help clinicians guide multimodal interventions to preserve muscle and optimize survival outcomes.

Putative markers for the detection of early-stage bladder cancer selected by urine metabolomics.

BACKGROUND: Early detection of bladder cancer remains challenging because patients with early-stage bladder cancer usually have no incentive to take cytology or cystoscopy tests if they are asymptomatic. Our goal is to find non-invasive marker candidates that may help us gain insight into the metabolism of early-stage bladder cancer and be examined in routine health checks. RESULTS: We acquired urine samples from 124 patients diagnosed with early-stage bladder cancer or hernia (63 cancer patients and 61 controls). In which 100 samples were included in our marker discovery cohort, and the remaining 24 samples were included in our independent test cohort. We obtained metabolic profiles of 922 compounds of the samples by gas chromatography-mass spectrometry. Based on the metabolic profiles of the marker discovery cohort, we selected marker candidates using Wilcoxon rank-sum test with Bonferroni correction and leave-one-out cross-validation; we further excluded compounds detected in less than 60% of the bladder cancer samples. We finally selected eight putative markers. The abundance of all the eight markers in bladder cancer samples was high but extremely low in hernia samples. Moreover, the up-regulation of these markers might be in association with sugars and polyols metabolism. CONCLUSIONS: In the present study, comparative urine metabolomics selected putative metabolite markers for the detection of early-stage bladder cancer. The suggested relations between early-stage bladder cancer and sugars and polyols metabolism may create opportunities for improving the detection of bladder cancer.

Morpho-physiological properties and connectivity of vasoactive intestinal polypeptide-expressing interneurons in the mouse hippocampal dentate gyrus.

The hippocampus is a key brain structure for cognitive and emotional functions. Among the hippocampal subregions, the dentate gyrus (DG) is the first station that receives multimodal sensory information from the cortex. Local-circuit inhibitory GABAergic interneurons (INs) regulate the excitation-inhibition balance in the DG principal neurons (PNs) and therefore are critical for information processing. Similar to PNs, GABAergic INs also receive distinct inhibitory inputs. Among various classes of INs, vasoactive intestinal polypeptide-expressing (VIP+ ) INs preferentially target other INs in several brain regions and thereby directly modulate the GABAergic system. However, the morpho-physiological characteristics and postsynaptic targets of VIP+ INs in the DG are poorly understood. Here, we report that VIP+ INs in the mouse DG are highly heterogeneous based on their morpho-physiological characteristics. In approximately two-thirds of morphologically reconstructed cells, their axons ramify in the hilus. The remaining cells project their axons exclusively to the molecular layer (15%), to both the molecular layer and hilus (10%), or throughout the entire DG layers (8%). Generally, VIP+ INs display variable intrinsic properties and discharge patterns without clear correlation with their morphologies. Finally, VIP+ INs are recruited with a long latency in response to theta-band cortical inputs and preferentially innervate GABAergic INs over glutamatergic PNs. In summary, VIP+ INs in the DG are composed of highly diverse subpopulations and control the DG output via disinhibition.

Association between real-world home blood pressure measurement patterns and blood pressure variability among older individuals with hypertension: A community-based blood pressure variability study.

Home blood pressure (BP) monitoring is a useful tool for hypertension management. BP variability (BPV) has been associated with an increased risk of cardiovascular events. However, little is known about the correlation between BPV and different measurement patterns of long-term home BP monitoring. This longitudinal cohort study aimed to assess the associations between dynamic BP measurement patterns and BPV. A total of 1128 participants (mean age, 77.4 ± 9.3 years; male, 51%) with 23 269 behavior measuring units were included. We used sliding window sampling to classify the home BP data with a regular 6-month interval into units in a sliding manner until the data are not continuous. Three measurement patterns (stable frequent [SF], stable infrequent [SI], and unstable [US]) were assessed based on the home BP data obtained within the first 3 months of the study, and the data in the subsequent 3 months were used to assess the BPV of that unit. We used linear mixed-effects model to assess the association between BP measurement patterns and BPV with adjustment for possible confounding factors including average BP. Average real variability and coefficient variability were used as measures of the BPV. No significant differences were observed in average BP between the SF, SI, and US patterns. However, BPV in the SF group was significantly lower than that in the US and SI groups (all p-values < .05). The BPV in SI and US groups was not significantly different. A stable and frequent BP measuring pattern was independently associated with a lower BPV.

Stout camphor tree genome fills gaps in understanding of flowering plant genome evolution.

We present reference-quality genome assembly and annotation for the stout camphor tree (Cinnamomum kanehirae (Laurales, Lauraceae)), the first sequenced member of the Magnoliidae comprising four orders (Laurales, Magnoliales, Canellales and Piperales) and over 9,000 species. Phylogenomic analysis of 13 representative seed plant genomes indicates that magnoliid and eudicot lineages share more recent common ancestry than monocots. Two whole-genome duplication events were inferred within the magnoliid lineage: one before divergence of Laurales and Magnoliales and the other within the Lauraceae. Small-scale segmental duplications and tandem duplications also contributed to innovation in the evolutionary history of Cinnamomum. For example, expansion of the terpenoid synthase gene subfamilies within the Laurales spawned the diversity of Cinnamomum monoterpenes and sesquiterpenes.

Combining membrane proteomics and computational three-way pathway analysis revealed signalling pathways preferentially regulated in human iPSCs and human ESCs.

Owing to the clinical potential of human induced pluripotent stem cells (hiPSCs) in regenerative medicine, a thorough examination of the similarities and differences between hiPSCs and human embryonic stem cells (hESCs) has become indispensable. Moreover, as the important roles of membrane proteins in biological signalling, functional analyses of membrane proteome are therefore promising. In this study, a pathway analysis by the bioinformatics tool GSEA was first performed to identify significant pathways associated with the three comparative membrane proteomics experiments: hiPSCs versus precursor human foreskin fibroblasts (HFF), hESCs versus precursor HFF, and hiPSCs versus hESCs. A following three-way pathway comparison was conducted to identify the differentially regulated pathways that may contribute to the differences between hiPSCs and hESCs. Our results revealed that pathways related to oxidative phosphorylation and focal adhesion may undergo incomplete regulations during the reprogramming process. This hypothesis was supported by another public proteomics dataset to a certain degree. The identified pathways and their core enriched proteins could serve as the starting point to explore the possible ways to make hiPSCs closer to hESCs.

A knowledge-based T2-statistic to perform pathway analysis for quantitative proteomic data.

Approaches to identify significant pathways from high-throughput quantitative data have been developed in recent years. Still, the analysis of proteomic data stays difficult because of limited sample size. This limitation also leads to the practice of using a competitive null as common approach; which fundamentally implies genes or proteins as independent units. The independent assumption ignores the associations among biomolecules with similar functions or cellular localization, as well as the interactions among them manifested as changes in expression ratios. Consequently, these methods often underestimate the associations among biomolecules and cause false positives in practice. Some studies incorporate the sample covariance matrix into the calculation to address this issue. However, sample covariance may not be a precise estimation if the sample size is very limited, which is usually the case for the data produced by mass spectrometry. In this study, we introduce a multivariate test under a self-contained null to perform pathway analysis for quantitative proteomic data. The covariance matrix used in the test statistic is constructed by the confidence scores retrieved from the STRING database or the HitPredict database. We also design an integrating procedure to retain pathways of sufficient evidence as a pathway group. The performance of the proposed T2-statistic is demonstrated using five published experimental datasets: the T-cell activation, the cAMP/PKA signaling, the myoblast differentiation, and the effect of dasatinib on the BCR-ABL pathway are proteomic datasets produced by mass spectrometry; and the protective effect of myocilin via the MAPK signaling pathway is a gene expression dataset of limited sample size. Compared with other popular statistics, the proposed T2-statistic yields more accurate descriptions in agreement with the discussion of the original publication. We implemented the T2-statistic into an R package T2GA, which is available at https://github.com/roqe/T2GA.

Exploring the expression bar code of SAA variants for gastric cancer detection.

We reported an integrated platform to explore serum protein variant pattern in cancer and its utility as a new class of biomarker panel for diagnosis. On the model study of serum amyloid A (SAA), we employed nanoprobe-based affinity mass spectrometry for enrichment, identification and quantitation of SAA variants from serum of 105 gastric cancer patients in comparison with 54 gastritis patients, 54 controls, and 120 patients from other cancer. The result revealed surprisingly heterogeneous and most comprehensive SAA bar code to date, which comprises 24 SAA variants including SAA1- and SAA2-encoded products, polymorphic isoforms, N-terminal-truncated forms, and three novel SAA oxidized isotypes, in which the variant-specific peptide sequence were also confirmed by LC-MS/MS. A diagnostic model was developed for dimension reduction and computational classification of the 24 SAA-variant bar code, providing good discrimination (AUC = 0.85 ± 3.2E-3) for differentiating gastric cancer group from gastritis and normal groups (sensitivity, 0.76; specificity, 0.81) and was validated with external validation cohort (sensitivity, 0.71; specificity, 0.74). Our platform not only shed light on the occurrence and modification extent of under-represented serum protein variants in cancer, but also suggested a new concept of diagnostic platform by serum protein variant profile.

SheddomeDB: the ectodomain shedding database for membrane-bound shed markers.

BACKGROUND: A number of membrane-anchored proteins are known to be released from cell surface via ectodomain shedding. The cleavage and release of membrane proteins has been shown to modulate various cellular processes and disease pathologies. Numerous studies revealed that cell membrane molecules of diverse functional groups are subjected to proteolytic cleavage, and the released soluble form of proteins may modulate various signaling processes. Therefore, in addition to the secreted protein markers that undergo secretion through the secretory pathway, the shed membrane proteins may comprise an additional resource of noninvasive and accessible biomarkers. In this context, identifying the membrane-bound proteins that will be shed has become important in the discovery of clinically noninvasive biomarkers. Nevertheless, a data repository for biological and clinical researchers to review the shedding information, which is experimentally validated, for membrane-bound protein shed markers is still lacking. RESULTS: In this study, the database SheddomeDB was developed to integrate publicly available data of the shed membrane proteins. A comprehensive literature survey was performed to collect the membrane proteins that were verified to be cleaved or released in the supernatant by immunological-based validation experiments. From 436 studies on shedding, 401 validated shed membrane proteins were included, among which 199 shed membrane proteins have not been annotated or validated yet by existing cleavage databases. SheddomeDB attempted to provide a comprehensive shedding report, including the regulation of shedding machinery and the related function or diseases involved in the shedding events. In addition, our published tool ShedP was embedded into SheddomeDB to support researchers for predicting the shedding event on unknown or unrecorded membrane proteins. CONCLUSIONS: To the best of our knowledge, SheddomeDB is the first database for the identification of experimentally validated shed membrane proteins and currently may provide the most number of membrane proteins for reviewing the shedding information. The database included membrane-bound shed markers associated with numerous cellular processes and diseases, and some of these markers are potential novel markers because they are not annotated or validated yet in other databases. SheddomeDB may provide a useful resource for discovering membrane-bound shed markers. The interactive web of SheddomeDB is publicly available at http://bal.ym.edu.tw/SheddomeDB/ .

Metabolite marker discovery for the detection of bladder cancer by comparative metabolomics.

Bladder cancer is one of the most common urinary tract carcinomas in the world. Urine metabolomics is a promising approach for bladder cancer detection and marker discovery since urine is in direct contact with bladder epithelia cells; metabolites released from bladder cancer cells may be enriched in urine samples. In this study, we applied ultra-performance liquid chromatography time-of-flight mass spectrometry to profile metabolite profiles of 87 samples from bladder cancer patients and 65 samples from hernia patients. An OPLS-DA classification revealed that bladder cancer samples can be discriminated from hernia samples based on the profiles. A marker discovery pipeline selected six putative markers from the metabolomic profiles. An LLE clustering demonstrated the discriminative power of the chosen marker candidates. Two of the six markers were identified as imidazoleacetic acid whose relation to bladder cancer has certain degree of supporting evidence. A machine learning model, decision trees, was built based on the metabolomic profiles and the six marker candidates. The decision tree obtained an accuracy of 76.60%, a sensitivity of 71.88%, and a specificity of 86.67% from an independent test.

Wiring Specificity and Synaptic Diversity in the Mouse Lateral Central Amygdala.

The central amygdala (CeA) nucleus, a subcortical structure composed of mostly GABA-releasing (GABAergic) neurons, controls fear expression via projections to downstream targets in the hypothalamus and brainstem. The CeA consists of the lateral (CeL) and medial (CeM) subdivisions. The CeL strongly gates information transfer to the CeM, the main output station of the amygdala, but little is known about the functional organization of local circuits in this region. Using cluster analysis, we identified two major electrophysiologically distinct CeL neuron classes in mouse amygdala slices, the early-spiking (ES) and late-spiking (LS) neurons. These two classes displayed distinct autaptic transmission. Compared with LS neurons, ES neurons had strong and depressing autapses, which enhanced spike-timing precision. With multiple patch-clamp recordings, we found that CeL neurons made chemical, but not electrical, synapses. Analysis of individual connections revealed cannabinoid type 1 receptor-mediated suppression of the ES, but not of the LS cell output synapse. More interestingly, the efficacy of the ES→LS or LS→ES synapse was ~2-fold greater than that of the LS→LS or ES→ES synapse. When tested at 20 Hz, synapses between different neurons, but not within the same class, were markedly depressing and were more powerful to sculpt activity of postsynaptic neurons. Moreover, neurons of different classes also form synapses with higher degree of connectivity. We demonstrate that ES and LS neurons represent two functionally distinct cell classes in the CeL and interactions between presynaptic and postsynaptic neurons dictate synaptic properties between neurons. SIGNIFICANCE STATEMENT: The central lateral amygdala (CeL) is a key node in fear circuits, but the functional organization of local circuits in this region is largely unknown. The CeL consists of mostly GABAergic inhibitory neurons with different functional and molecular features. Here, we report that the presynaptic cell class determines functional properties of autapses and cannabinoid-mediated modulation of synaptic transmission between neurons, whereas presynaptic versus postsynaptic cell classes dictate the connectivity, efficacy, and dynamics of GABAergic synapses between any two neurons. The wiring specificity and synaptic diversity have a great impact on neuronal output in amygdala inhibitory networks. Such synaptic organizing principles advance our understanding of the significance of physiologically defined neuronal phenotypes in amygdala inhibitory networks.

Prioritization of cancer marker candidates based on the immunohistochemistry staining images deposited in the human protein atlas.

Cancer marker discovery is an emerging topic in high-throughput quantitative proteomics. However, the omics technology usually generates a long list of marker candidates that requires a labor-intensive filtering process in order to screen for potentially useful markers. Specifically, various parameters, such as the level of overexpression of the marker in the cancer type of interest, which is related to sensitivity, and the specificity of the marker among cancer groups, are the most critical considerations. Protein expression profiling on the basis of immunohistochemistry (IHC) staining images is a technique commonly used during such filtering procedures. To systematically investigate the protein expression in different cancer versus normal tissues and cell types, the Human Protein Atlas is a most comprehensive resource because it includes millions of high-resolution IHC images with expert-curated annotations. To facilitate the filtering of potential biomarker candidates from large-scale omics datasets, in this study we have proposed a scoring approach for quantifying IHC annotation of paired cancerous/normal tissues and cancerous/normal cell types. We have comprehensively calculated the scores of all the 17219 tested antibodies deposited in the Human Protein Atlas based on their accumulated IHC images and obtained 457110 scores covering 20 different types of cancers. Statistical tests demonstrate the ability of the proposed scoring approach to prioritize cancer-specific proteins. Top 100 potential marker candidates were prioritized for the 20 cancer types with statistical significance. In addition, a model study was carried out of 1482 membrane proteins identified from a quantitative comparison of paired cancerous and adjacent normal tissues from patients with colorectal cancer (CRC). The proposed scoring approach demonstrated successful prioritization and identified four CRC markers, including two of the most widely used, namely CEACAM5 and CEACAM6. These results demonstrate the potential of this scoring approach in terms of cancer marker discovery and development. All the calculated scores are available at http://bal.ym.edu.tw/hpa/.

SecretePipe: a screening pipeline for secreted proteins with competence to identify potential membrane-bound shed markers.

The identification of secreted protein markers has been receiving great attention as part of the trend toward noninvasive biomarker discovery. In addition, certain cell membrane proteins are known to be released into the extracellular milieu via ectodomain shedding. As membrane proteins play an essential role in signaling pathways and because most of the cancer biomarkers approved by the FDA today are membrane shed proteins, a tool that can correctly predict these class shed proteins is valuable. In this study, an in-house predictor, ShedP, was developed to predict the ectodomain shedding events of membrane proteins. ShedP is the first computational method to our knowledge to allow shed membrane protein prediction. By integrating ShedP with other state-of-the-art predictors, a screening pipeline, SecretePipe, has been created that is able to identify secreted nonmembrane proteins on the basis of signal peptides and to identify released membrane proteins on the basis of ectodomain shedding. The predictive results using secretome data sets revealed that SecretePipe outperformed other state-of-the-art secreted protein predictors. When evaluated against released membrane proteins, SecretePipe performed better than other predictors in identifying membrane-bound released proteins due to the presence of ShedP. SecretePipe showed a great potential in assisting the identification of membrane-bound shed markers in biomarker discovery.

A linear-time algorithm for reconstructing zero-recombinant haplotype configuration on a pedigree.

BACKGROUND: When studying genetic diseases in which genetic variations are passed on to offspring, the ability to distinguish between paternal and maternal alleles is essential. Determining haplotypes from genotype data is called haplotype inference. Most existing computational algorithms for haplotype inference have been designed to use genotype data collected from individuals in the form of a pedigree. A haplotype is regarded as a hereditary unit and therefore input pedigrees are preferred that are free of mutational events and have a minimum number of genetic recombinational events. These ideas motivated the zero-recombinant haplotype configuration (ZRHC) problem, which strictly follows the Mendelian law of inheritance, namely that one haplotype of each child is inherited from the father and the other haplotype is inherited from the mother, both without any mutation. So far no linear-time algorithm for ZRHC has been proposed for general pedigrees, even though the number of mating loops in a human pedigree is usually very small and can be regarded as constant. RESULTS: Given a pedigree with n individuals, m marker loci, and k mating loops, we proposed an algorithm that can provide a general solution to the zero-recombinant haplotype configuration problem in O(kmn + k2m) time. In addition, this algorithm can be modified to detect inconsistencies within the genotype data without loss of efficiency. The proposed algorithm was subject to 12000 experiments to verify its performance using different (n, m) combinations. The value of k was uniformly distributed between zero and six throughout all experiments. The experimental results show a great linearity in terms of execution time in relation to input size when both n and m are larger than 100. For those experiments where n or m are less than 100, the proposed algorithm runs very fast, in thousandth to hundredth of a second, on a personal desktop computer. CONCLUSIONS: We have developed the first deterministic linear-time algorithm for the zero-recombinant haplotype configuration problem. Our experimental results demonstrated the linearity of its execution time in relation to the input size. The proposed algorithm can be modified to detect inconsistency within the genotype data without loss of efficiency and is expected to be able to handle recombinant and missing data with further extension.

Multi-Q: a fully automated tool for multiplexed protein quantitation.

The iTRAQ labeling method combined with shotgun proteomic techniques represents a new dimension in multiplexed quantitation for relative protein expression measurement in different cell states. To expedite the analysis of vast amounts of spectral data, we present a fully automated software package, called Multi-Q, for multiplexed iTRAQ-based quantitation in protein profiling. Multi-Q is designed as a generic platform that can accommodate various input data formats from search engines and mass spectrometer manufacturers. To calculate peptide ratios, the software automatically processes iTRAQ's signature peaks, including peak detection, background subtraction, isotope correction, and normalization to remove systematic errors. Furthermore, Multi-Q allows users to define their own data-filtering thresholds based on semiempirical values or statistical models so that the computed results of fold changes in peptide ratios are statistically significant. This feature facilitates the use of Multi-Q with various instrument types with different dynamic ranges, which is an important aspect of iTRAQ analysis. The performance of Multi-Q is evaluated with a mixture of 10 standard proteins and human Jurkat T cells. The results are consistent with expected protein ratios and thus demonstrate the high accuracy, full automation, and high-throughput capability of Multi-Q as a large-scale quantitation proteomics tool. These features allow rapid interpretation of output from large proteomic datasets without the need for manual validation. Executable Multi-Q files are available on Windows platform at http://ms.iis.sinica.edu.tw/Multi-Q/.

RIBRA--an error-tolerant algorithm for the NMR backbone assignment problem.

We develop an iterative relaxation algorithm called RIBRA for NMR protein backbone assignment. RIBRA applies nearest neighbor and weighted maximum independent set algorithms to solve the problem. To deal with noisy NMR spectral data, RIBRA is executed in an iterative fashion based on the quality of spectral peaks. We first produce spin system pairs using the spectral data without missing peaks, then the data group with one missing peak, and finally, the data group with two missing peaks. We test RIBRA on two real NMR datasets, hbSBD and hbLBD, and perfect BMRB data (with 902 proteins) and four synthetic BMRB data which simulate four kinds of errors. The accuracy of RIBRA on hbSBD and hbLBD are 91.4% and 83.6%, respectively. The average accuracy of RIBRA on perfect BMRB datasets is 98.28%, and 98.28%, 95.61%, 98.16%, and 96.28% on four kinds of synthetic datasets, respectively.

GANA--a genetic algorithm for NMR backbone resonance assignment.

NMR data from different experiments often contain errors; thus, automated backbone resonance assignment is a very challenging issue. In this paper, we present a method called GANA that uses a genetic algorithm to automatically perform backbone resonance assignment with a high degree of precision and recall. Precision is the number of correctly assigned residues divided by the number of assigned residues, and recall is the number of correctly assigned residues divided by the number of residues with known human curated answers. GANA takes spin systems as input data and uses two data structures, candidate lists and adjacency lists, to assign the spin systems to each amino acid of a target protein. Using GANA, almost all spin systems can be mapped correctly onto a target protein, even if the data are noisy. We use the BioMagResBank (BMRB) dataset (901 proteins) to test the performance of GANA. To evaluate the robustness of GANA, we generate four additional datasets from the BMRB dataset to simulate data errors of false positives, false negatives and linking errors. We also use a combination of these three error types to examine the fault tolerance of our method. The average precision rates of GANA on BMRB and the four simulated test cases are 99.61, 99.55, 99.34, 99.35 and 98.60%, respectively. The average recall rates of GANA on BMRB and the four simulated test cases are 99.26, 99.19, 98.85, 98.87 and 97.78%, respectively. We also test GANA on two real wet-lab datasets, hbSBD and hbLBD. The precision and recall rates of GANA on hbSBD are 95.12 and 92.86%, respectively, and those of hbLBD are 100 and 97.40%, respectively.