State-of-the art data normalization methods improve NMR-based metabolomic analysis

Extracting biomedical information from large metabolomic datasets by multivariate data analysis is of considerable complexity. Common challenges include among others screening for differentially produced metabolites, estimation of fold changes, and sample classification. Prior to these analysis steps, it is important to minimize contributions from unwanted biases and experimental variance. This is the goal of data preprocessing. In this work, different data normalization methods were compared systematically employing two different datasets generated by means of nuclear magnetic resonance (NMR) spectroscopy. To this end, two different types of normalization methods were used, one aiming to remove unwanted sample-to-sample variation while the other adjusts the variance of the different metabolites by variable scaling and variance stabilization methods. The impact of all methods tested on sample classification was evaluated on urinary NMR fingerprints obtained from healthy volunteers and patients suffering from autosomal polycystic kidney disease (ADPKD). Performance in terms of screening for differentially produced metabolites was investigated on a dataset following a Latin-square design, where varied amounts of 8 different metabolites were spiked into a human urine matrix while keeping the total spike-in amount constant. In addition, specific tests were conducted to systematically investigate the influence of the different preprocessing methods on the structure of the analyzed data. In conclusion, preprocessing methods originally developed for DNA microarray analysis, in particular, Quantile and Cubic-Spline Normalization, performed best in reducing bias, accurately detecting fold changes, and classifying samples.

     

Comparison of the crystal structures of genetically engineered human manganese superoxide dismutase and manganese superoxide dismutase from Thermus thermophilus: differences in dimer-dimer interaction.

The three-dimensional X-ray structure of a recombinant human mitochondrial manganese superoxide dismutase (MnSOD) (chain length 198 residues) was determined by the method of molecular replacement using the related structure of MnSOD from Thermus thermophilus as a search model. This tetrameric human MnSOD crystallizes in space group P2(1)2(1)2 with a dimer in the asymmetric unit (Wagner, U.G., Werber, M.M., Beck, Y., Hartman, J.R., Frolow, F., & Sussman, J.L., 1989, J. Mol. Biol. 206, 787-788). Refinement of the protein structure (3,148 atoms with Mn and no solvents), with restraints maintaining noncrystallographic symmetry, converged at an R-factor of 0.207 using all data from 8.0 to 3.2 A resolution and group thermal parameters. The monomer-monomer interactions typical of bacterial Fe- and Mn-containing SODs are retained in the human enzyme, but the dimer-dimer interactions that form the tetramer are very different from those found in the structure of MnSOD from T. thermophilus. In human MnSOD one of the dimers is rotated by 84 degrees relative to its equivalent in the thermophile enzyme. As a result the monomers are arranged in an approximately tetrahedral array, the dimer-dimer packing is more intimate than observed in the bacterial MnSOD from T. thermophilus, and the dimers interdigitate. The metal-ligand interactions, determined by refinement and verified by computation of omit maps, are identical to those observed in T. thermophilus MnSOD.     

Identification by denaturing high-performance liquid chromatography of numerous polymorphisms in a candidate region for multiple sclerosis susceptibility

Genetic association analysis of candidate regions where evidence of linkage has accumulated is becoming a key issue in the study of complex diseases. A high density of markers, at least one per centimorgan, is required to improve the chances of observing linkage disequilibrium with disease alleles. A recently available single nucleotide polymorphism (SNP) map designed to cover the whole genome provides an average density of one marker per 2 cM. In the present study we show that the number of markers can be approximately doubled in a selected region, thus reaching a density suitable for association studies, by applying a completely automated technique for polymorphism detection, denaturing high-performance liquid chromatography (DHPLC). A systematic search for SNPs was performed in the region 5ptel-q13, where weak but convergent evidence for linkage with multiple sclerosis has accumulated. Screening for polymorphisms was performed on 124 sequence tagged sites (STSs) in the 3'UTR ends of expressed sequence tags totaling about 30,000 bp. Thirty SNPs in 28 STSs were found with less than 10% overlap with the markers already detected in the same region. The data confirm the validity of the approach using DHPLC on expressed gene sequences tagged by a set of standard commercially available primers.     

The 1.15A crystal structure of the Staphylococcus aureus methionyl-aminopeptidase and complexes with triazole based inhibitors

Methionyl aminopeptidases (MetAPs) represent a unique class of protease that are responsible for removing the N-terminal methionine residue from proteins and peptides. There are two major classes of MetAPs (type I and type II) described and each class can be subdivided into two subclasses. Eukaryotes contain both the type I and type II MetAPs, whereas prokaryotes possess only the type I enzyme. Due to the physiological importance of these enzymes there is considerable interest in inhibitors to be used as antiangiogenic and antimicrobial agents. Here, we describe the 1.15A crystal structure of the Staphylococcus aureus MetAP-I as an apo-enzyme and its complexes with various 1,2,4-triazole-based derivatives at high-resolution. The protein has a typical "pita-bread" fold as observed for the other MetAP structures. The inhibitors bind in the active site with the N1 and N2 atoms of the triazole moiety complexing two divalent ions. The 1,2,4-triazols represent a novel class of potent non-peptidic inhibitors for the MetAP-Is.     

A back migration from Asia to sub-Saharan Africa is supported by high-resolution analysis of human Y-chromosome haplotypes

The variation of 77 biallelic sites located in the nonrecombining portion of the Y chromosome was examined in 608 male subjects from 22 African populations. This survey revealed a total of 37 binary haplotypes, which were combined with microsatellite polymorphism data to evaluate internal diversities and to estimate coalescence ages of the binary haplotypes. The majority of binary haplotypes showed a nonuniform distribution across the continent. Analysis of molecular variance detected a high level of interpopulation diversity (PhiST=0.342), which appears to be partially related to the geography (PhiCT=0.230). In sub-Saharan Africa, the recent spread of a set of haplotypes partially erased pre-existing diversity, but a high level of population (PhiST=0.332) and geographic (PhiCT=0.179) structuring persists. Correspondence analysis shows that three main clusters of populations can be identified: northern, eastern, and sub-Saharan Africans. Among the latter, the Khoisan, the Pygmies, and the northern Cameroonians are clearly distinct from a tight cluster formed by the Niger-Congo-speaking populations from western, central western, and southern Africa. Phylogeographic analyses suggest that a large component of the present Khoisan gene pool is eastern African in origin and that Asia was the source of a back migration to sub-Saharan Africa. Haplogroup IX Y chromosomes appear to have been involved in such a migration, the traces of which can now be observed mostly in northern Cameroon.     

Frequentist estimation of coalescence times from nucleotide sequence data using a tree-based partition

This article proposes a method of estimating the time to the most recent common ancestor (TMRCA) of a sample of DNA sequences. The method is based on the molecular clock hypothesis, but avoids assumptions about population structure. Simulations show that in a wide range of situations, the point estimate has small bias and the confidence interval has at least the nominal coverage probability. We discuss conditions that can lead to biased estimates. Performance of this estimator is compared with existing methods based on the coalescence theory. The method is applied to sequences of Y chromosomes and mtDNAs to estimate the coalescent times of human male and female populations.     

A decade of high-resolution liquid chromatography of nucleic acids on styrene-divinylbenzene copolymers

The introduction of alkylated, nonporous poly-(styrene-divinylbenzene) microparticles in 1992 enabled the subsequent development of denaturing HPLC that has emerged as the most sensitive screening method for mutations to date. Denaturing HPLC has provided unprecedented insight into human origins and prehistoric migrations, accelerated the cloning of genes involved in mono- and polygenic traits, and facilitated the mutational analysis of more than a hundred candidate genes of human disease. A significant step toward increased sample-throughput and information content was accomplished by the recent introduction of monolithic poly(styrene-divinylbenzene) capillary columns. They have enabled the construction of capillary arrays amenable to multiplex analysis of fluorescent dye-labeled nucleic acids by laser-induced fluorescence detection. Hyphenation of denaturing HPLC with electrospray ionization mass spectrometry, on the other hand, has allowed the direct elucidation of the chemical nature of DNA variation and determination of phase of multiple alleles on a chromosome.     

Sequence variation and the biological function of genes: methodological and biological considerations

Single nucleotide polymorphisms (SNPs) are expected to facilitate the chromosomal mapping and eventual cloning of genetic determinants of complex quantitative phenotypes. To date, more than 2.5 million non-redundant human SNPs have been reported in the public domain, of which approximately 100000 have been validated by either independent investigators or by independent methods. Equally impressive is the myriad of methods developed for allelic discrimination. Nevertheless, reports of successful applications of SNPs to genome-wide linkage analysis of both mono- and polygenic traits are rare and limited to a few model organisms, that provide affordable platforms to test both novel methodological and biological concepts at a whole-genome scale under conditions that can be reasonably controlled. Progress in the analysis of SNPs needs to be complemented by methods that allow the systematic elucidation of both primary and secondary phenotypes of genes. Importantly, observations made in one species may very well be of immediate applicability to other species including human. This is particularly true for conserved biological processes such as mitochondrial respiration and DNA repair.     

A simple procedure for the analysis of single nucleotide polymorphisms facilitates map-based cloning in Arabidopsis

We developed a modified allele-specific PCR procedure for assaying single nucleotide polymorphisms (SNPs) and used the procedure (called SNAP for single-nucleotide amplified polymorphisms) to generate 62 Arabidopsis mapping markers. SNAP primers contain a single base pair mismatch within three nucleotides from the 3' end of one allele (the specific allele) and in addition have a 3' mismatch with the nonspecific allele. A computer program called SNAPER was used to facilitate the design of primers that generate at least a 1,000-fold difference in the quantity of the amplification products from the specific and nonspecific SNP alleles. Because SNAP markers can be readily assayed by electrophoresis on standard agarose gels and because a public database of over 25,000 SNPs is available between the Arabidopsis Columbia and Landsberg erecta ecotypes, the SNAP method greatly facilitates the map-based cloning of Arabidopsis genes defined by a mutant phenotype.