AMADA: analysis of microarray data

Bioinformatics and Translation Elongation

Codon usage depends on mutation bias, tRNA-mediated selection, and the need for high efficiency and accuracy in translation. One codon in a synonymous codon family is often strongly over-used, especially in highly expressed genes, which often leads to a high dN/dS ratio because dS is very small. Many different codon usage indices have been proposed to measure codon usage and codon adaptation. Sense codon could be misread by release factors and stop codons misread by tRNAs, which also contribute to codon usage in rare cases. This chapter outlines the conceptual framework on codon evolution, illustrates codon-specific and gene-specific codon usage indices, and presents their applications. A new index for codon adaptation that accounts for background mutation bias (Index of Translation Elongation) is presented and contrasted with codon adaptation index (CAI) which does not consider background mutation bias. They are used to re-analyze data from a recent paper claiming that translation elongation efficiency matters little in protein production. The reanalysis disproves the claim.

DAMBE6: New Tools for Microbial Genomics, Phylogenetics, and Molecular Evolution

DAMBE is a comprehensive software workbench for data analysis in molecular biology, phylogenetics, and evolution. Several important new functions have been added since version 5 of DAMBE: 1) comprehensive genomic profiling of translation initiation efficiency of different genes in different prokaryotic species, 2) a new index of translation elongation (I_TE) that takes into account both tRNA-mediated selection and background mutation on codon–anticodon adaptation, 3) a new and accurate phylogenetic approach based on pairwise alignment only, which is useful for highly divergent sequences from which a reliable multiple sequence alignment is difficult to obtain. Many other functions have been updated and improved including PWM for motif characterization, Gibbs sampler for de novo motif discovery, hidden Markov models for protein secondary structure prediction, self-organizing map for nonlinear clustering of transcriptomic data, comprehensive sequence alignment, and phylogenetic functions. DAMBE features a graphic, user-friendly and intuitive interface, and is freely available from http://dambe.bio.uottawa.ca.

Bioinformatics and Drug Discovery

Bioinformatic analysis can not only accelerate drug target identification and drug candidate screening and refinement, but also facilitate characterization of side effects and predict drug resistance. High-throughput data such as genomic, epigenetic, genome architecture, cistromic, transcriptomic, proteomic, and ribosome profiling data have all made significant contribution to mechanismbased drug discovery and drug repurposing. Accumulation of protein and RNA structures, as well as development of homology modeling and protein structure simulation, coupled with large structure databases of small molecules and metabolites, paved the way for more realistic protein-ligand docking experiments and more informative virtual screening. I present the conceptual framework that drives the collection of these high-throughput data, summarize the utility and potential of mining these data in drug discovery, outline a few inherent limitations in data and software mining these data, point out news ways to refine analysis of these diverse types of data, and highlight commonly used software and databases relevant to drug discovery.

DAMBE5: a comprehensive software package for data analysis in molecular biology and evolution

Since its first release in 2001 as mainly a software package for phylogenetic analysis, data analysis for molecular biology and evolution (DAMBE) has gained many new functions that may be classified into six categories: 1) sequence retrieval, editing, manipulation, and conversion among more than 20 standard sequence formats including MEGA, NEXUS, PHYLIP, GenBank, and the new NeXML format for interoperability, 2) motif characterization and discovery functions such as position weight matrix and Gibbs sampler, 3) descriptive genomic analysis tools with improved versions of codon adaptation index, effective number of codons, protein isoelectric point profiling, RNA and protein secondary structure prediction and calculation of minimum folding energy, and genomic skew plots with optimized window size, 4) molecular phylogenetics including sequence alignment, testing substitution saturation, distance-based, maximum parsimony, and maximum-likelihood methods for tree reconstructions, testing the molecular clock hypothesis with either a phylogeny or with relative-rate tests, dating gene duplication and speciation events, choosing the best-fit substitution models, and estimating rate heterogeneity over sites, 5) phylogeny-based comparative methods for continuous and discrete variables, and 6) graphic functions including secondary structure display, optimized skew plot, hydrophobicity plot, and many other plots of amino acid properties along a protein sequence, tree display and drawing by dragging nodes to each other, and visual searching of the maximum parsimony tree. DAMBE features a graphic, user-friendly, and intuitive interface and is freely available from http://dambe.bio.uottawa.ca (last accessed April 16, 2013).

Comparative transcriptome analysis of Lactococcus lactis subsp. cremoris strains under conditions simulating Cheddar cheese manufacture

Gene expression in response to technological variations can influence fermentation and flavor generation in Cheddar cheese, and can vary from one lactococcal strain to another, perceived as differences in starter performance. The aim of this study was to determine the influence of cheese cooking temperature at 38 °C and salting on the transcriptional profiles of four closely related strains of L. lactis subsp. cremoris under simulated conditions of Cheddar cheese manufacture. Two responses could be distinguished, a core gene expression, corresponding to the common response of all strains and strain-specific response during the Cheddar simulating process. For the core gene expression after heating of inoculated milk at 38 °C, two groups of differentially expressed genes were identified: i) stress response and ii) carbohydrate and amino acid metabolism. The response to combined stresses of heat, acid and salt resulted in: i) general decrease of functions linked to cell division and metabolism, ii) specific responses related to stress such as the induction of genes coding for chaperones and proteases and iii) expression of prophage lytic systems for certain strains. Strain-specific responses were mainly observed in three of the four tested strains. These responses were the induction of genes related to osmotic stress or the release of CodY repression leading to the activation of oligopeptide transporters as well as the bcaT gene, related to amino acid degradation for the production of flavor. Comparing transcriptomes provides a core expression profile that contributes to understanding gene expression responses to environmental variations. The strain-specific responses identify predictive markers for the transcriptional state of starter strains before they enter the cheese ripening phase.

Daily genetic profiling indicates JAK/STAT signaling promotes early hepatic stellate cell transdifferentiation

AIM: To identify signaling pathways and genes that initiate and commit hepatic stellate cells (HSCs) to transdifferentiation.

METHODS: Primary HSCs were isolated from male Sprague-Dawley rats and cultured on plastic for 0-10 d. Gene expression was assessed daily (quiescent to day 10 culture-activation) by real time polymerase chain reaction and data clustered using AMADA software. The significance of JAK/STAT signaling to HSC transdifferentiation was determined by treating cells with a JAK2 inhibitor.

RESULTS: Genetic cluster analyses, based on expression of these 21 genes, showed similar expression profiles on days 1-3, days 5 and 6, and days 7-10, while freshly isolated cells (day Q) and day 4 cells were genotypically distinct from any of the other days. Additionally, gene expression clustering revealed strong upregulation of interleukin-6, JAK2 and STAT3 mRNA in the early stages of activation. Inhibition of the JAK/STAT signaling pathway impeded the morphological transdifferentiation of HSCs which correlated with decreased mRNA expression of several profibrotic genes including collagens, α-SMA, PDGFR and TGFβR.

CONCLUSION: These data demonstrate unique clustered genetic profiles during the daily progression of HSC transdifferentiation and that JAK/STAT signaling may be critical in the early stages of transdifferentiation.

A preliminary matrix-assisted laser desorption/ionization time-of-flight approach for the characterization of Italian lentil varieties

Matrix-assisted laser desorption/ionization mass spectrometry (MALDI-MS) was used in this study to obtain protein fingerprints of seven different lentil varieties, to characterize their differences and similarities. Two different matrices have been tested in order to obtain reproducible and significant mass spectra. Extraction with water containing 0.1% of trifluoroacetic acid has been used as preparative step to obtain hydrophilic protein samples of lentil seeds. The obtained MALDI protein profiles identified clear differences between the seven studied lentil varieties. Moreover, considering the high complexity of the obtained MALDI spectra, multivariate techniques of data analysis were employed to find further classification details. These multivariate analyses confirmed the possibility of a clear classification of the seven lentil varieties, indicating that the proposed procedure can be a valid taxonomic tool, and a method to certify the origin of lentils, useful for high added value lentils (Italian lentils).

Low molecular weight proteins in urines from healthy subjects as well as diabetic, nephropathic and diabetic-nephropathic patients: a MALDI study

Urine samples from healthy subjects as well as diabetic, nephropathic and diabetic-nephropathic patients were analyzed by matrix assisted laser desorption/ionization (MALDI) mass spectrometry in order to establish evidence of some possible differences in the peptide profile related to the pathological states. Multivariate analysis suggested the possibility of a distinction among the considered groups of patients. Some differences have been found, in particular, in the relative abundances of three ions at m/z 1912, 1219 and 2049. For these reasons, further investigation was carried out by MALDI/TOF/TOF to determine the sequence of these peptides and, consequently, to individuate their possible origin. By this approach, the peptide at m/z 1912 was found to originate from uromodulin, and its lower expression in the case of nephropathy can be well related to the pathological condition. Ions at m/z 2049 and 1219 originate from the collagen alpha-1(I) chain precursor and from the collagen alpha-5 (IV) chain precursor, respectively, and, also in this case, their different expressions can be related to the pathologies under investigation. The obtained data seem to indicate that urine is an interesting biological fluid to investigate on the peptide profile and to obtain, consequently, information on the dismetabolism activated by specific pathologies.

Utilizing a D-Amino Acid as a Drug Carrier for Antineoplastic Nitrogen Mustard Groups

The evaluation of an alkylating nitrogen mustard agent that utilizes D-alanine as a drug carrier for three chloroethyl substituents (CICH2CH2-) is shown. Various important pharmacological properties were determined including polar surface area, partition coefficient, molar volume, polarizability, numbers of -OH and -NH2 groups, and aqueous solubility. The synthetic approach utilizes 1,2-dichloroethane reaction with the primary amine of D-alanine resulting in chloroethyl (CICH2CH2-) substituents. A nitrogen mustard group results, and this agent showed alkylation activity in aqueous solution directed toward a nucleophilic primary amine group. Kinetics of alkylation activity is determined by placing p-chloroaniline and the nitrogen mustard agent in sodium bicarbonate buffered aqueous solution at physiological pH 7.4 and 37 degrees C. Samples taken from the test solution are injected with fluorescamine that reacts specifically with primary amine functional groups. Absorbance measurements obtained at 400 nm by UV-Vis spectrometer indicates the relative amounts of nonalkylated p-chloroaniline in the test solution. The D-alanine nitrogen mustard agent effectively alkylated the nucleophilic primary amine of p-chloroaniline with zero-order kinetics. The rate constant was determined to be 7.445E-04 mol/l/min. Formula weight, polar surface area, Log P, molar volume, violations of Rule of 5 for D-alanine mustard agent are 276.59, 29.543 angstroms2, 1.605, 222.3 cm3, and zero violations, respectively. Cluster analysis of molecular properties showed D-alanine mustard agent to be quite similar to cyclophosphamide. This agent showed good druglikeness and zero violations of the Rule of 5, indicating good bioavailability. Molecular properties calculated for the mustard agent are numerically comparable to some clinical anticancer drugs.

Global identification of bursicon-regulated genes in Drosophila melanogaster

Background

Bursicon is a heterodimer neuropeptide responsible for regulating cuticle sclerotization and wing expansion in several insect species. Recent studies indicate that the action of bursicon is mediated by a specific G protein-coupled receptor DLGR2 and the cAMP/PKA signaling pathway. However, little is known regarding the genes that are regulated by bursicon. The identification of bursicon-regulated genes is the focus of this investigation.

Results

We used DNA microarray analysis to identify bursicon-regulated genes in neck-ligated flies (Drosophila melanogaster) that received recombinant bursicon (r-bursicon). Fifty four genes were found to be regulated by bursicon 1 h post r-bursicon injection, 52 being up-regulated and 2 down-regulated while 33 genes were influenced by r-bursicon 3 h post-injection (24 up-regulated and 9 down-regulated genes). Analysis of these genes by inference from the fly database revealed that these genes encode proteins with diverse functions, including cell signaling, gene transcription, DNA/RNA binding, ion trafficking, proteolysis-peptidolysis, metabolism, cytoskeleton formation, immune response and cell-adhesion. Twenty eight genes randomly selected from the microarray-identified list were verified by real time PCR (qPCR) which supported the microarray data. Temporal response studies of 13 identified and verified genes by qPCR revealed that the temporal expression patterns of these genes are consistent with the microarray data.

Conclusion

Using r-bursicon, we identified 87 genes that are regulated by bursicon, 30 of which have no previously known function. Most importantly, all genes randomly selected from the microarray-identified list were verified by real time PCR. Temporal analysis of 13 verified genes revealed that the expression of these genes was indeed induced by bursicon and correlated well with the cuticle sclerotization process. The composite data suggest that these genes play important roles in regulating the cuticle sclerotization and wing expansion processes. The data obtained here will form the basis for future studies aimed at elucidating the exact mechanisms upstream from the secretion of bursicon and its binding to target cells.

Can crude oils be distinguished by different component distribution obtained by laser desorption ionization mass spectrometry and evaluated by chemometrics?

Fourteen different crude oil samples were analyzed as collected by mass spectrometry. For two of them the results obtained by means of different techniques, i.e. electrospray ionization, atmospheric pressure chemical ionization, laser desorption/ionization, were compared. The latter technique leads to the best results: even if unable to give specific information on heteroatom-containing components, it allows a general view to be obtained of the panorama of the oil composition in terms of molecular weight distribution. The statistical evaluation of the mass spectrometry data by multivariate techniques, such as cluster analysis (Average Linkage) and correspondence analysis, allows evidence for the differences and similarities among the crude oils under study.

Mass spectrometry in the characterization of ambers. I. Studies of amber samples of different origin and ages by laser desorption ionization, atmospheric pressure chemical ionization and atmospheric pressure photoionization mass spectrometry

Amber is a fossil resin constituted of organic polymers derived through complex maturation processes of the original plant resin. A classification of eight samples of amber of different geological age (Miocene to Triassic) and geographical origin is here proposed using direct mass spectrometric techniques, i.e. laser desorption ionization (LDI), atmospheric pressure chemical ionization (APCI) and atmospheric pressure photoionization (APPI), in order to obtain a fingerprint related to the amber origin. Differences and similarities were detected among the spectra with the four methods, showing quite complex spectra, full of ionic species in the mass range investigated (up to m/z 2000). The evaluation required statistical analysis involving multivariate techniques. Cluster analysis or principal component analysis (PCA) generally did not show a clear clustering with respect to the age of samples, except for the APPI method, which allowed a satisfying clustering. Using the total ion current (TIC) obtained by the different analytical approaches on equal quantities of the different amber samples and plotted against the age, the only significant correlation appeared to be that involving APPI. To validate the method, four amber samples from Cretaceous of Spain were analyzed. Also in this case a significant correlation with age was found only with APPI data. PCA obtained with TIC values from all the MS methods showed a fair grouping of samples, according to their age. Three main clusters could be detected, belonging to younger, intermediate and older fossil resins, respectively. This MS analysis on crude amber, either solid or extract, followed by appropriate multivariate statistical evaluation, can provide useful information on amber age. The best results are those obtained by APPI, indicating that the quantity of amber soluble components that can be photoionized decreases with increasing age, in agreement with the formation of highly stable, insoluble polymers.

Overview of microarray analysis of gene expression and its applications to cervical cancer investigation

Cervical cancer is one of the leading female cancers in Taiwan and ranks as the fifth cause of cancer death in the female population. Human papillomavirus has been established as the causative agent for cervical neoplasia and cervical cancer. However, the tumor biology involved in the prognoses of different cell types in early cancers and tumor responses to radiation in advanced cancers remain largely unknown. The introduction of microarray technologies in the 1990s has provided genome-wide strategies for searching tens of thousands of genes simultaneously. In this review, we first summarize the two types of microarrays: oligonucleotides microarray and cDNA microarray. Then, we review the studies of functional genomics in cervical cancer. Gene expression studies that involved cervical cancer cell lines, cervical cells of cancer versus normal ectocervix, cancer tissues of different histology, radioresistant versus radiosensitive patients, and the combinatorial gene expression associated with chromosomal amplifications are discussed. In particular, CEACAM5 , TACSTD1 , S100P , and MSLN have shown to be upregulated in adenocarcinoma, and increased expression levels of CEACAM5 and TACSTD1 were significantly correlated with poorer patient outcomes. On the other hand, 35 genes, including apoptotic genes (e.g. BIK , TEGT , SSI-3 ), hypoxia-inducible genes (e.g. HIF1A , CA12 ), and tumor cell invasion and metastasis genes (e.g. CTSL , CTSB , PLAU , CD44 ), have been noted to echo the hypothesis that increased tumor hypoxia leads to radiation resistance in cervical cancer during radiation.

Multivariate analysis of matrix-assisted laser desorption/ionization mass spectrometric data related to glycoxidation products of human globins in nephropathic patients

To clarify the possible pathogenetic role of oxidation products originated from the glycation of proteins, human globins from nephropathic patients have been studied by matrix-assisted laser desorption/ionization mass spectrometry (MALDI), revealing not only unglycated and monoglycated globins, but also a series of different species. For the last ones, structural assignments were tentatively done on the basis of observed masses and expectations for the Maillard reaction pattern. Consequently, they must be considered only propositive, and the discussion which will follow must be considered in this view. In our opinion this approach does not seem to compromise the intended diagnostic use of the data because distinctions are valid even if the assignments are uncertain. We studied nine healthy subjects and 19 nephropathic patients and processed the data obtained from the MALDI spectra using a multivariate analysis. Our results showed that multivariate analytical techniques enable differential aspects of the profile of molecular species to be identified in the blood of end stage nephropathic patients. A correct grouping can be achieved by principal component analysis (PCA) and the results suggest that several products involved in carbonyl stress exist in nephropathic patients. These compounds may have a relevant role as specific markers of the pathological state.

Advanced glycation end products and antioxidant status in type 2 diabetic patients with and without peripheral artery disease

OBJECTIVE

Advanced glycation end products (AGEs), pentosidine and malondialdehyde (MDA), are elevated in type 2 diabetic subjects with coronary and carotid angiopathy. We investigated the relationship of AGEs, MDA, total reactive antioxidant potentials (TRAPs), and vitamin E in type 2 diabetic patients with and without peripheral artery disease (PAD).

RESEARCH DESIGN AND METHODS

AGEs, pentosidine, MDA, TRAP, vitamin E, and ankle-brachial index (ABI) were measured in 99 consecutive type 2 diabetic subjects and 20 control subjects.

RESULTS

AGEs, pentosidine, and MDA were higher and vitamin E and TRAP were lower in patients with PAD (ABI <0.9) than in patients without PAD (ABI >0.9) (P < 0.001). After multiple regression analysis, a correlation between AGEs and pentosidine, as independent variables, and ABI, as the dependent variable, was found in both patients with and without PAD (r = 0.9198, P < 0.001 and r = 0.5764, P < 0.001, respectively) but not in control subjects. When individual regression coefficients were evaluated, only that due to pentosidine was confirmed as significant. For patients with PAD, considering TRAP, vitamin E, and MDA as independent variables and ABI as the dependent variable produced an overall significant regression (r = 0.6913, P < 0.001). The regression coefficients for TRAP and vitamin E were not significant, indicating that the model is best explained by a single linear regression between MDA and ABI. These findings were also confirmed by principal component analysis.

CONCLUSIONS

Results show that pentosidine and MDA are strongly associated with PAD in type 2 diabetic patients.

Multivariate analysis approach to the plasma protein profile of patients with advanced colorectal cancer

The aim of the present study was to identify the pattern of plasma protein species of interest as markers of colorectal cancer (CRC). Using matrix-assisted laser desorption/ionization-mass spectrometry (MALDI-MS), the plasma protein profile was determined in nine stage IV CRC patients (study group) and nine clean-colon healthy subjects (control group). Multivariate analysis methods were employed to identify distinctive disease patterns at protein spectrum. In the study and control groups, cluster analysis (CA) on the complete MALDI-MS spectra plasma protein profile showed a distinction between CRC patients and healthy subjects, thus allowing the identification of the most discriminating ionic species. Principal component analysis (PCA) and linear discriminant analysis (LDA) yielded similar grouping results. LDA with leave-one-out cross validation achieved a correct classification rate of 89% in both the patients and the healthy subjects.

Metabolomic fingerprinting of plant extracts

The standardization and quality control of plant extracts is an important topic, in particular, when such extracts are used for medicinal purposes. Consequently, the development of fast and effective analytical methods for metabolomic fingerprinting of plant extracts is of high interest. In this investigation, electrospray mass spectrometry (ESI-MS) and (1)H NMR techniques were employed with further statistical analyses of the acquired data. The results showed that negative ion mode ESI-MS is particularly effective for characterization of plant extracts. Different samples of the same species appear well-clustered and separated from the other species. To verify the effectiveness of the method, two other batches of extracts from a species, in which the principal components were already identified (Cynara scolymus), were analyzed, and the components that were verified by the principal component analysis (PCA) were found to be within the region identified as characteristic of Cynara Scolymus extracts. The data from extracts of the other species were well separated from those pertaining to the species previously characterized. Only the case of a species that was strictly correlated from a botanical point of view, with extracts that were previously analyzed, showed overlapping.

Applying pattern recognition methods and structure property correlations to determine drug carrier potential of nicotinic acid and analogize to dihydropyridine

Multivariate methods are utilized to compare nicotinic acid and dihydropyridine as a drug carrier. Nicotinic acid and dihydropyridine form ester groups on 10 beta-lactam antibiotics with an oxymethyl group forming a linkage between the antibiotic and the drug carrier (nicotinic acid or dihydropyridine). Calculated molecular properties are analyzed by self-organizing tree algorithm (SOTA), bivariate regression, cluster analysis, factor analysis, discriminant analysis, hierarchical classification, and principal coordinates analysis. Ten important pharmacological properties for each of the nicotinic acid and dihydropyridine antibiotic derivatives are numerically similar and highly correlated. Calculated molecular properties include molar refractivity, molar volume, parachor, index of refraction, partition coefficient (log P), polarizability, and polar surface area. Dermal permeability coefficients (Kp) for nicotinic acid derivatives are similar to values for dihydropyridine derivatives. Dermal permeabliity coefficients analyzed by hierarchical classification and SOTA analysis were shown to be closely interrelated and highly correlated. Ten properties of the nicotinic acid and dihydropyridine were compared by Passing-Bablok regression analysis and shown to be highly correlated (r=0.9879). Box plot analysis of 10 properties, inclusive of both groups of derivatives, showed narrow ranges in values. Cluster analysis of derivative properties showed the nicotinic acid derivatives to be highly similar to the dihydropyridine derivatives of the same antibiotics. Cluster analysis was performed by single linkage, complete linkage, and centroid linkage. Factor analysis showed the nicotinic acid derivatives to be interrelated and similar to the dihydropyridine derivatives. Discriminant analysis performed on all derivatives formed a single highly cohesive and non-differentiated cluster, demonstrating strong similarity among nicotinic acid and dihydropyridine derivatives. Principal coordinates analysis (determines similarity) of Kp values showed high similarity between the nicotinic acid and dihydropyridine antibiotic derivatives.

Thematic review series: The pathogenesis of atherosclerosis. Toward a biological network for atherosclerosis

The goal of systems biology is to define all of the elements present in a given system and to create an interaction network between these components so that the behavior of the system, as a whole and in parts, can be explained under specified conditions. The elements constituting the network that influences the development of atherosclerosis could be genes, pathways, transcript levels, proteins, or physiologic traits. In this review, we discuss how the integration of genetics and technologies such as transcriptomics and proteomics, combined with mathematical modeling, may lead to an understanding of such networks.