Jürgen Claesen

Identification of novel causative genes determining the complex trait of high ethanol tolerance in yeast using pooled-segregant whole-genome sequence analysis

High ethanol tolerance is an exquisite characteristic of the yeast Saccharomyces cerevisiae, which enables this microorganism to dominate in natural and industrial fermentations. Up to now, ethanol tolerance has only been analyzed in laboratory yeast strains with moderate ethanol tolerance. The genetic basis of the much higher ethanol tolerance in natural and industrial yeast strains is unknown. We have applied pooled-segregant whole-genome sequence analysis to map all quantitative trait loci (QTL) determining high ethanol tolerance. We crossed a highly ethanol-tolerant segregant of a Brazilian bioethanol production strain with a laboratory strain with moderate ethanol tolerance. Out of 5974 segregants, we pooled 136 segregants tolerant to at least 16% ethanol and 31 segregants tolerant to at least 17%. Scoring of SNPs using whole-genome sequence analysis of DNA from the two pools and parents revealed three major loci and additional minor loci. The latter were more pronounced or only present in the 17% pool compared to the 16% pool. In the locus with the strongest linkage, we identified three closely located genes affecting ethanol tolerance: MKT1, SWS2, and APJ1, with SWS2 being a negative allele located in between two positive alleles. SWS2 and APJ1 probably contained significant polymorphisms only outside the ORF, and lower expression of APJ1 may be linked to higher ethanol tolerance. This work has identified the first causative genes involved in high ethanol tolerance of yeast. It also reveals the strong potential of pooled-segregant sequence analysis using relatively small numbers of selected segregants for identifying QTL on a genome-wide scale.

BRAIN: A Universal Tool for High-Throughput Calculations of the Isotopic Distribution for Mass Spectrometry

This Letter presents the R-package implementation of the recently introduced polynomial method for calculating the aggregated isotopic distribution called BRAIN (Baffling Recursive Algorithm for Isotopic distributioN calculations). The algorithm is simple, easy to understand, highly accurate, fast, and memory-efficient. The method is based on the application of the Newton-Girard theorem and Viètes formulae to the polynomial coding of different aggregated isotopic variants. As a result, an elegant recursive equation is obtained for computing the occurrence probabilities of consecutive aggregated isotopic peaks. Additionally, the algorithm also allows calculating the center-masses of the aggregated isotopic variants. We propose an implementation which is suitable for high-throughput processing and easily customizable for application in different areas of mass spectral data analyses. A case study demonstrates how the R-package can be applied in the context of protein research, but the software can be also used for calculating the isotopic distribution in the context of lipidomics, metabolomics, glycoscience, or even space exploration. More materials, i.e., reference manual, vignette, and the package itself are available at Bioconductor online (http://www.bioconductor.org/packages/release/bioc/html/BRAIN.html).

Computational methods and challenges in hydrogen/deuterium exchange mass spectrometry

Hydrogen/Deuterium exchange (HDX) has been applied, since the 1930s, as an analytical tool to study the structure and dynamics of (small) biomolecules. The popularity of using HDX to study proteins increased drastically in the last two decades due to the successful combination with mass spectrometry (MS). Together with this growth in popularity, several technological advances have been made, such as improved quenching and fragmentation. As a consequence of these experimental improvements and the increased use of protein-HDXMS, large amounts of complex data are generated, which require appropriate analysis. Computational analysis of HDXMS requires several steps. A typical workflow for proteins consists of identification of (non-)deuterated peptides or fragments of the protein under study (local analysis), or identification of the deuterated protein as a whole (global analysis); determination of the deuteration level; estimation of the protection extent or exchange rates of the labile backbone amide hydrogen atoms; and a statistically sound interpretation of the estimated protection extent or exchange rates. Several algorithms, specifically designed for HDX analysis, have been proposed. They range from procedures that focus on one specific step in the analysis of HDX data to complete HDX workflow analysis tools. In this review, we provide an overview of the computational methods and discuss outstanding challenges.

Simultaneous mapping of multiple gene loci with pooled segregants.

The analysis of polygenic, phenotypic characteristics such as quantitative traits or inheritable diseases remains an important challenge. It requires reliable scoring of many genetic markers covering the entire genome. The advent of high-throughput sequencing technologies provides a new way to evaluate large numbers of single nucleotide polymorphisms (SNPs) as genetic markers. Combining the technologies with pooling of segregants, as performed in bulked segregant analysis (BSA), should, in principle, allow the simultaneous mapping of multiple genetic loci present throughout the genome. The gene mapping process, applied here, consists of three steps: First, a controlled crossing of parents with and without a trait. Second, selection based on phenotypic screening of the offspring, followed by the mapping of short offspring sequences against the parental reference. The final step aims at detecting genetic markers such as SNPs, insertions and deletions with next generation sequencing (NGS). Markers in close proximity of genomic loci that are associated to the trait have a higher probability to be inherited together. Hence, these markers are very useful for discovering the loci and the genetic mechanism underlying the characteristic of interest. Within this context, NGS produces binomial counts along the genome, i.e., the number of sequenced reads that matches with the SNP of the parental reference strain, which is a proxy for the number of individuals in the offspring that share the SNP with the parent. Genomic loci associated with the trait can thus be discovered by analyzing trends in the counts along the genome. We exploit the link between smoothing splines and generalized mixed models for estimating the underlying structure present in the SNP scatterplots.

A hidden Markov-model for gene mapping based on whole-genome next generation sequencing data

Abstract The analysis of polygenic, phenotypic characteristics such as quantitative traits or inheritable diseases requires reliable scoring of many genetic markers covering the entire genome. The advent of high-throughput sequencing technologies provides a new way to evaluate large numbers of single nucleotide polymorphisms as genetic markers. Combining the technologies with pooling of segregants, as performed in bulk segregant analysis, should, in principle, allow the simultaneous mapping of multiple genetic loci present throughout the genome. We propose a hidden Markov-model to analyze the marker data obtained by the bulk segregant next generation sequencing. The model includes several states, each associated with a different probability of observing the same/different nucleotide in an offspring as compared to the parent. The transitions between the molecular markers imply transitions between the states of the model. After estimating the transition probabilities and state-related probabilities of nucleotide (dis)similarity, the most probable state for each SNP is selected. The most probable states can then be used to indicate which genomic regions may be likely to contain trait-related genes. The application of the model is illustrated on the data from a study of ethanol tolerance in yeast. Software is written in R. R-functions, R-scripts and documentation are available on www.ibiostat.be/software/bioinformatics.

A nonhomogeneous hidden markov model for gene mapping based on next-generation sequencing data.

The analysis of polygenetic characteristics for mapping quantitative trait loci (QTL) remains an important challenge. QTL analysis requires two or more strains of organisms that differ substantially in the (poly-)genetic trait of interest, resulting in a heterozygous offspring. The offspring with the trait of interest is selected and subsequently screened for molecular markers such as single-nucleotide polymorphisms (SNPs) with next-generation sequencing. Gene mapping relies on the co-segregation between genes and/or markers. Genes and/or markers that are linked to a QTL influencing the trait will segregate more frequently with this locus. For each identified marker, observed mismatch frequencies between the reads of the offspring and the parental reference strains can be modeled by a multinomial distribution with the probabilities depending on the state of an underlying, unobserved Markov process. The states indicate whether the SNP is located in a (vicinity of a) QTL or not. Consequently, genomic loci associated with the QTL can be discovered by analyzing hidden states along the genome. The aforementioned hidden Markov model assumes that the identified SNPs are equally distributed along the chromosome and does not take the distance between neighboring SNPs into account. The distance between the neighboring SNPs could influence the chance of co-segregation between genes and markers. To address this issue, we propose a nonhomogeneous hidden Markov model with a transition matrix that depends on a set of distance-varying observed covariates. The application of the model is illustrated on the data from a study of ethanol tolerance in yeast.

Comparison of the Mahalanobis Distance and Pearson's χ 2 Statistic as Measures of Similarity of Isotope Patterns

AbstractTo extract a genuine peptide signal from a mass spectrum, an observed series of peaks at a particular mass can be compared with the isotope distribution expected for a peptide of that mass. To decide whether the observed series of peaks is similar to the isotope distribution, a similarity measure is needed. In this short communication, we investigate whether the Mahalanobis distance could be an alternative measure for the commonly employed Pearson’s χ2 statistic. We evaluate the performance of the two measures by using a controlled MALDI-TOF experiment. The results indicate that Pearson’s χ2 statistic has better discriminatory performance than the Mahalanobis distance and is a more robust measure. Figureᅟ

The Analysis of Peptide-Centric Mass-Spectrometry Data Utilizing Information About the Expected Isotope Distribution

In shotgun proteomics, much attention and instrument time is dedicated to the generation of tandem mass spectra. These spectra contain information about the fragments of, ideally, one peptide and are used to infer the amino acid sequence of the scrutinized peptide. This type of spectrum acquisition is called a product ion scan, tandem MS, or MS2 spectrum. Another type of spectrum is the, often overlooked, precursor ion scan or MS1 spectrum that catalogs all ionized analytes present in a mass spectrometer. While MS2 spectra are important to identify the peptides and proteins in the sample, MS1 spectra provide valuable information about the quantity of the analyte. In this chapter, we describe some properties of MS1 spectra, such as the isotope distribution, and how these properties can be employed for low-level signal processing to reduce data complexity and as a tool for quality assurance. Furthermore, we describe some cases in which advanced modeling of the isotope distribution can be used in quantitative proteomics analysis.

POPPeT: a New Method to Predict the Protection Factor of Backbone Amide Hydrogens

AbstractHydrogen exchange (HX) has become an important tool to monitor protein structure and dynamics. The interpretation of HX data with respect to protein structure requires understanding of the factors that influence exchange. Simulated protein structures can be validated by comparing experimental deuteration profiles with the profiles derived from the modeled protein structure. To do this, we propose here a new method, POPPeT, for protection factor prediction based on protein motions that enable HX. By comparing POPPeT with two existing methods, the phenomenological approximation and COREX, we show enhanced predictability measured at both protection factor and deuteration level. This method can be subsequently used by modeling strategies for protein structure prediction. Graphical Abstractᅟ

Gene Expression Signature Differentiates Histology But Not Progression Status of Early-Stage NSCLC

Advances in molecular analyses based on high-throughput technologies can contribute to a more accurate classification of non–small cell lung cancer (NSCLC), as well as a better prediction of both the disease course and the efficacy of targeted therapies. Here we set out to analyze whether global gene expression profiling performed in a group of early-stage NSCLC patients can contribute to classifying tumor subtypes and predicting the disease prognosis. Gene expression profiling was performed with the use of the microarray technology in a training set of 108 NSCLC samples. Subsequently, the recorded findings were validated further in an independent cohort of 44 samples. We demonstrated that the specific gene patterns differed significantly between lung adenocarcinoma (AC) and squamous cell lung carcinoma (SCC) samples. Furthermore, we developed and validated a novel 53-gene signature distinguishing SCC from AC with 93% accuracy. Evaluation of the classifier performance in the validation set showed that our predictor classified the AC patients with 100% sensitivity and 88% specificity. We revealed that gene expression patterns observed in the early stages of NSCLC may help elucidate the histological distinctions of tumors through identification of different gene-mediated biological processes involved in the pathogenesis of histologically distinct tumors. However, we showed here that the gene expression profiles did not provide additional value in predicting the progression status of the early-stage NSCLC. Nevertheless, the gene expression signature analysis enabled us to perform a reliable subclassification of NSCLC tumors, and it can therefore become a useful diagnostic tool for a more accurate selection of patients for targeted therapies.