PubMed Health⌕ Search

Biomedical subjects

Elja Arjas

Publications and source records attributed to Elja Arjas.

5 recordsLinked to original sources

Bayesian hierarchical model for correcting signal saturation in microarrays using pixel intensities.

Pixel saturation occurs when the pixel intensity exceeds the scanner upper threshold of detection and the recorded pixel intensity is then truncated at the threshold. Truncation of the pixel intensity causes the estimates of gene expression (i.e., intensity) to be biased. Microarray experiments are commonly affected by saturated pixels; as a result all higher level analyses are made on these biased gene expression estimates. In this paper, we propose a method for improving the quality of the signal for cDNA microarrays by making use of several scans at varying scanner sensitivities. For each spot, pixel level intensity readings are given as input to a Bayesian hierarchical model. The model uses the pixel intensities of the spot to provide a posterior distribution of the true expression level of the corresponding genes. The parameters of the hierarchical model are estimated jointly with these expression levels, thus performing an integrated analysis of the measurement data. The method improves in all ranges the accuracy with which intensities can be estimated and extends the dynamic range of measured gene expression at the high end. The method is generic and can be applied to data from any organism and for imaging with any scanner. Results from a real data set illustrate an improved precision in the estimation of the expression of genes compared to what can be achieved by applying standard methods and using only a single scan.

Bayes Theorem↗

Backward simulation of ancestors of sampled individuals.

If the population is large and the sampling mechanism is random, the coalescent is commonly used to model the haplotypes in the sample. Ordered genotypes can then be formed by random matching of the derived haplotypes. However, this approach is not realistic when (1) there is departure from random mating (e.g., dominant individuals in breeding populations or monogamy in humans), or (2) the population is small and/or the individuals in the sample are ascertained by applying some particular non-random sampling scheme, as is usually the case when considering the statistical modeling and analysis of pedigree data. For such situations, we present here a data generation method where an ancestral graph with non-overlapping generations is first generated backwards in time, using ideas from coalescent theory. Alleles are randomly assigned to the founders, and subsequently the gene flow over the entire genome is simulated forwards in time by dropping alleles down the graph according to recombination model without interference. The parameters controlling the mating behavior of generated individuals in the graph (degree of monogamy) can be tuned in order to match a particular demographic situation, without restriction to simple random mating. The performance of the approach is illustrated with a simulation example. The software (written in C-language) is freely available for research purposes at http://www.rni.helsinki.fi/~dag/.

Female↗

Bayesian integrated functional analysis of microarray data.

MOTIVATION: The statistical analysis of microarray data usually proceeds in a sequential manner, with the output of the previous step always serving as the input of the next one. However, the methods currently used in such analyses do not properly account for the fact that the intermediate results may not always be correct, then leading to cumulating error in the inferences drawn based on such steps. RESULTS: Here we show that, by an application of hierarchical Bayesian methodology, this sequential procedure can be replaced by a single joint analysis, while systematically accounting for the uncertainties in this process. Moreover, we can also integrate relevant functional information available from databases into such an analysis, thereby increasing the reliability of the biological conclusions that are drawn. We illustrate these points by analysing real data and by showing that the genes can be divided into categories of interest, with the defining characteristic depending on the biological question that is considered. We contend that the proposed method has advantages at two levels. First, there are gains in the statistical and biological results from the analysis of this particular dataset. Second, it opens up new possibilities in analysing microarray data in general.

Algorithms↗

Increasing incidence of Type 1 diabetes--role for genes?

BACKGROUND: The incidence of Type 1 diabetes (T1DM) is increasing fast in many populations. The reasons for this are not known, although an increase in the penetrance of the diabetes-associated alleles, through changes in the environment, might be the most plausible mechanism. After the introduction of insulin treatment in 1930s, an increase in the pool of genetically susceptible individuals has been suggested to contribute to the increase in the incidence of Type 1 diabetes. RESULTS: To explore this hypothesis, the authors formulate a simple population genetic model for the incidence change driven by non-Mendelian transmission of a single susceptibility factor, either allele(s) or haplotype(s). A Poisson mixture model is used to model the observed number of cases. Model parameters were estimated by maximizing the log-likelihood function. Based on the Finnish incidence data 1965-1996 the point estimate of the transmission probability was 0.998. Given our current knowledge of the penetrance of the most diabetic gene variants in the HLA region and their transmission probabilities, this value is exceedingly unrealistic. CONCLUSIONS: As a consequence, non-Mendelian transmission of diabetic allele(s)/haplotype(s) if present, could explain only a small part of the increase in incidence in Finland. Hence, the importance of other, probably environmental factors modifying the disease incidence is emphasized.

Diabetes Mellitus, Type 1↗