SEARCH · Search PubMed

Results for “variational inference”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 397 records · Page 22Linked to original sources

PubMed2026

Genetic control of local mutation rates.

Mutations are the source of evolutionary novelty but also the cause of genetic diseases and cancer. Mutation rates are known to be heterogeneous along the genome, however the extent to which local mutation rates vary among individuals in a population and are genetically determined is unknown. To test this, we analyzed the chromosomal distribution of somatic mutations in cell lines from 1,662 individuals, controlling for the confounding effects of DNA replication timing on local mutation rates and of trans-acting modulators on global mutation rates. We describe substantial interindividual variation in mutation rates across the human genome. By comparing mutation-rate variation to individuals' genotypes, we identified 35 instances in which polymorphic alleles in the population associate with somatic mutation rates in their vicinity. We call these mutation quantitative trait loci (mutQTLs). mutQTLs associated with somatic mutations in lymphoblastoid cell lines and in chronic lymphocytic leukemia, and with germline genetic variants. Two of the four mutQTLs inferred to be associated with germline mutation-rate variation were located within large clusters of zinc-finger genes and transposable elements, where they functioned as cis-mutators conferring an increased rate of mutation in their vicinity. mutQTLs provide a portal into the evolution of mutation rate heterogeneity across the genome and across individuals.

Humans↗

PubMed2000

Saddlepoint approximations for small sample logistic regression problems.

Double saddlepoint approximations provide quick and accurate approximations to exact conditional tail probabilities in a variety of situations. This paper describes the use of these approximations in two logistic regression problems. An investigation of regression analysis of the log-odds ratio in a sequence or set of 2x2 tables via simulation studies shows that in practical settings the saddlepoint methods closely approximate exact conditional inference. The double saddlepoint approximation in the test for trend in a sequence of binomial random variates is also shown, via simulation studies, to be an effective approximation to exact conditional inference.

Logistic Models↗

PubMed2004

Coalescent-based association mapping and fine mapping of complex trait loci.

We outline a general coalescent framework for using genotype data in linkage disequilibrium-based mapping studies. Our approach unifies two main goals of gene mapping that have generally been treated separately in the past: detecting association (i.e., significance testing) and estimating the location of the causative variation. To tackle the problem, we separate the inference into two stages. First, we use Markov chain Monte Carlo to sample from the posterior distribution of coalescent genealogies of all the sampled chromosomes without regard to phenotype. Then, averaging across genealogies, we estimate the likelihood of the phenotype data under various models for mutation and penetrance at an unobserved disease locus. The essential signal that these models look for is that in the presence of disease susceptibility variants in a region, there is nonrandom clustering of the chromosomes on the tree according to phenotype. The extent of nonrandom clustering is captured by the likelihood and can be used to construct significance tests or Bayesian posterior distributions for location. A novelty of our framework is that it can naturally accommodate quantitative data. We describe applications of the method to simulated data and to data from a Mendelian locus (CFTR, responsible for cystic fibrosis) and from a proposed complex trait locus (calpain-10, implicated in type 2 diabetes).

Alleles↗

PubMed2005

Approximating the coalescent with recombination.

The coalescent with recombination describes the distribution of genealogical histories and resulting patterns of genetic variation in samples of DNA sequences from natural populations. However, using the model as the basis for inference is currently severely restricted by the computational challenge of estimating the likelihood. We discuss why the coalescent with recombination is so challenging to work with and explore whether simpler models, under which inference is more tractable, may prove useful for genealogy-based inference. We introduce a simplification of the coalescent process in which coalescence between lineages with no overlapping ancestral material is banned. The resulting process has a simple Markovian structure when generating genealogies sequentially along a sequence, yet has very similar properties to the full model, both in terms of describing patterns of genetic variation and as the basis for statistical inference.

Chromosomes↗

PubMed2000

Geographic patterns: how to identify them and why.

Geographic patterns of genetic diversity allow us to make inferences about population histories and the evolution of inherited disease. The statistical methods describing genetic variation in space, such as estimation of genetic variances, mapping of allele frequencies, and principal components analysis, have opened up the possibility to reconstruct demographic processes whose effects have been tested by a variety of approaches, including spatial autocorrelation, cladistic analyses, and simulations. These studies have significantly contributed to our understanding of human genetic variation; however, the molecular data that have accumulated since the mid-1980s have also created new complications. Reasons include the generally limited sample sizes, but, more generally, it is the nature of molecular variation itself that makes it necessary to develop and apply specific models and methods for the treatment of DNA data. The foreseeable diffusion of laboratory techniques for the rapid typing of many DNA markers will force us to change our approach to the study of human variation anyway, moving from the gene level toward the genome level. Because extensive variation among loci is the rule rather than the exception, an important practical tip is to be skeptical of inferences based on single-locus diversity.

Biological Evolution↗

PubMed1999

Generalized linear mixed models for meta-analysis.

We examine two strategies for meta-analysis of a series of 2 x 2 tables with the odds ratio modelled as a linear combination of study level covariates and random effects representing between-study variation. Penalized quasi-likelihood (PQL), an approximate inference technique for generalized linear mixed models, and a linear model fitted by weighted least squares to the observed log-odds ratios are used to estimate regression coefficients and dispersion parameters. Simulation results demonstrate that both methods perform adequate approximate inference under many conditions, but that neither method works well in the presence of highly sparse data. Under certain conditions with small cell frequencies the PQL method provides better inference.

Humans↗

PubMed1999

Differences in uniquantal amplitude between sites reduce uniquantal variance when few release sites are active.

In many types of central neurons, the coefficient of variation (CV) of stimulus-evoked uniquantal events inferred from quantal analysis is small, frequently less than 20%. In contrast, spontaneous putative uniquantal events (minis) from the same neurons are much more variable in amplitude, having a CV of roughly 50% or more. One explanation for this discrepancy is that, if the variance in mini amplitude were generated by differences between release sites, the small number of sites activated during stimulation would sometimes fortuitously have similar quantal amplitudes. Only in these fortuitous cases where uniquantal variance is small could quantal peaks be resolved, and therefore the uniquantal CV seen in the subset of cells where quantal analysis can be performed would systematically be much smaller than predicted by the mini distribution. We have explored this possibility by Monte Carlo simulation assuming all variance in mini amplitude to be due to intersite differences in uniquantal amplitude. We find that when a small number of release sites are activated under these conditions, there is a reduction in the expected uniquantal variance. However, the expected uniquantal CV is highly variable from one experiment to the next, and low uniquantal CVs are not expected to be seen often enough to account for the high frequency with which quantal peaks with a uniquantal CV < 20% are observed experimentally. We conclude that variance in mini amplitude between release sites cannot account for the small uniquantal CV seen in quantal analysis of many central synapses.

Hippocampus↗

PubMed1996

Leucine side-chain rotamers in a glycophorin A transmembrane peptide as revealed by three-bond carbon-carbon couplings and 13C chemical shifts.

We have used a spin-echo difference NMR pulse sequence to measure three-bond J couplings between lambda- and alpha-carbons of the leucine residues in a micelle-associated helical peptide dimer that corresponds to residues 62-101 of the transmembrane erythrocyte protein glycophorin A. The observed 3J couplings correlate strongly with the 13C chemical shift of the lambda-methyl groups, and within experimental error both the shift distribution of the methyl carbons and the variations in 3J can be accounted for by variations in side-chain rotamer populations. We infer that all leucine side chains in this peptide dimer are in fast exchange among chi 2 rotamers and sample two of the three possible rotameric states, even when the side chain forms part of the dimer interface. The observed correlation of chemical shift with couplings can be traced to a gamma-gauche interaction of methyl and alpha-carbons. This correlation may provide an alternate route to rotamer analysis in some protein systems.

Carbon Isotopes↗

PubMedDate not supplied

Xanthine dehydrogenase (XDH): episodic evolution of a "neutral" protein.

We investigated the evolution of xanthine dehydrogenase (Xdh) in 34 species from the three multicellular kingdoms, including one plant, two fungi, and three animal phyla, two classes of vertebrates, four orders of mammals, and two orders of insects. We adopted a model-based maximum-likelihood framework of inference. After accounting for among-site rate variation and heterogeneous nucleotide composition of the sequences using the discrete gamma distribution, and using nonhomogeneous nonstationary representations of the substitution process, the rate of amino acid replacement is 30.4 x 10(-10)/site/year when Drosophila species are compared but only approximately 18 x 10(-10)/site/year when comparisons are made between mammal orders, between insect orders, or between different animal phyla and approximately 11 x 10(-10)/site/year when comparisons are made between birds and mammals, between fungi, or between the three multicellular kingdoms. To account for these observations, the rate of amino acid replacement must have been eight or more times higher in some lineages and at some times than in others. Spastic evolution of Xdh appears to be related to the particularities of the genomes in which the locus is embedded.

Amino Acid Substitution↗

PubMed2004

Germs, genomes and genealogies.

Genetic diversity in pathogen species contains information about evolutionary and epidemiological processes, including the origins and history of disease, the nature of the selective forces acting on pathogen genes and the role of recombination in generating genetic novelty. Here, we review recent developments in these fields and compare the use of population genetic, or population-model based, approaches to phylogenetic, or population-model free, methodologies. We show how simple epidemiological models can be related to the ancestral, or coalescent, process underlying samples from pathogen species, enabling detailed inference about pathogen biology from patterns of molecular variation.

Journal Article↗

PubMed1998

Systematic relationships of some members of the genera Oesophagostomum and Chabertia (Nematoda: Chabertiidae) based on ribosomal DNA sequence data.

The present study characterised seven species of the Chabertiidae (Nematoda: Strongyloidea) belonging to either the subfamily Oesophagostominae (Oesophagostomum radiatum, Oesophagostomum venulosum, Oesophagostomum dentatum, Oesophagostomum quadrispinulatum, Oesophagostomum columbianum, Oesophagostomum bifurcum) or to the subfamily Chabertiinae (Chabertia ovina) by their second internal transcribed spacer rDNA sequence, assessed the extent of intraspecific variation and interspecific differences in the sequence, and inferred the phylogenetic relationship of C. ovina with respect to members of the Oesophagostominae. In both the phenetic and cladistic analyses of the sequence data, Chabertia was nested within Oesophagostomum, suggesting either that the species examined represent members of the same genus, or alternatively, that Oesophagostomum may represent more than one genus.

Animals↗

PubMed2002

Household economic strategies and nutritional anthropometry of women in American Samoa and highland Bolivia.

This study compares findings from research projects involving different genetic, environmental, and cultural contexts: a study of lifestyle and health from American Samoa (ASLS) and the Bolivian project. Reproduction and Ecology in Provincia Aroma (REPA). This paper presents analyses of varying economic strategies and their association with nutritional status indicators in each population. The ASLS sample includes 66 Samoan women and the REPA sample includes 210 Aymara women. Principle components analysis of household economic resources within each sample extracted two significant factors: one represents modernizing influences including education and occupational status, and the other represents ethnographically salient traditional economic behavior. The traditional pattern includes adding household members in Samoa and selling agricultural products in Bolivia. This analysis places each woman along two continua, traditional and modern, based on her household mobilization of economic resources, permitting an understanding of the patterns underlying household economic behavior that is not possible in univariate analyses of socioeconomic variables. For the Bolivian women the strategy involving more education and higher occupational status was associated with higher measures of several nutritional status indicators, including body mass index, arm muscle area, and peripheral skinfolds. But among the Samoan women, where substantial obesity was the norm, there were no significant differences in anthropometric measurements based on economic strategies. These data argue for the importance of directly measuring the potential consequences of variation in household economic strategies rather than merely inferring such, and of assessing ethnographically relevant aspects of household economic production rather than limiting analyses to non-context-specific economic indicators such as income. This focus on household strategy is likely to be fruitful especially where economic and nutritional conditions are marginal. The findings from Bolivia also support efforts in developing countries to improve girls' education, and thereby occupational prospects, as a means to improve their health status as women.

Adult↗

PubMed2002

Terminology-driven literature mining and knowledge acquisition in biomedicine.

In this paper we describe Tagged Information Management System (TIMS), an integrated knowledge management system for the domain of molecular biology and biomedicine, in which terminology-driven literature mining, knowledge acquisition (KA), knowledge integration (KI), and XML-based knowledge retrieval are combined using tag information management and ontology inference. The system integrates automatic terminology acquisition, term variation management, hierarchical term clustering, tag-based information extraction (IE), and ontology-based query expansion. TIMS supports introducing and combining different types of tags (linguistic and domain-specific, manual and automatic). Tag-based interval operations and a query language are introduced in order to facilitate KA and retrieval from XML documents. Through KA examples, we illustrate the way in which literature mining techniques can be utilised for knowledge discovery from documents.

Artificial Intelligence↗

PubMed1992

Heparin-induced conformational change and activation of mucus proteinase inhibitor.

Low molecular mass heparin (5.1 kDa) forms a tight complex with mucus proteinase inhibitor, the physiologic neutrophil elastase inhibitor of the upper respiratory tract. This binding strongly enhances the intrinsic fluorescence of the inhibitor and the rate of neutrophil elastase inhibitor association. One mole of this heparin fragment binds 1 mol of inhibitor with a Kd of 50 nM. From the variation of Kd with ionic strength, it is inferred that (i) 85% of the heparin--inhibitor binding energy i due to electrostatic interactions, (ii) about seven ionic interactions are involved in heparin--inhibitor binding. strength, it is inferred that (i) 85% of the heparin--inhibitor binding energy is due to electrostatic interactions, (ii) about seven ionic interactions are involved in heparin--inhibitor binding. and (iii), about one-third of low quantum yield of Trp30, the single tryptophan residue of the inhibitor, blue-shifts its maximum emission wavelength by 6 nm, decreases the acrylamide quenching rate constant by a factor of 4, and increases the mean intensity weighted lifetime by a factor of 2.5. These important spectroscopic changes evidence a heparin--induced conformational change of the inhibitor which buries Trp30 in a very hydrophobic environment. Heparin accelerates the inhibition of elastase in a concentration-dependent manner. When both enzyme and inhibitor are saturated by the polymer, the second-order association rate constant is 7.7 x 10(7) M-1 s-1, a value that is 27-fold higher than that measured with the free partners. This finding may have important physiologic and therapeutic bearing.

Amino Acid Sequence↗

PubMed2003

Determination of critical micelle concentration by hyper-rayleigh scattering.

The critical micelle concentration (CMC) of several surfactants that contain an NLO chromophore, either at the hydrocarbon tail, or at the hydrophilic headgroup, or even as a counterion, was determined by hyper-Rayleigh scattering (HRS). In all cases, the HRS signal exhibited a similar variation with surfactant concentration, wherein the CMC is inferred from a rather unprecedented drop in the signal intensity. This drop is attributed to the formation of small pre-micellar aggregates, whose concentrations become negligible above CMC. In addition, a probe molecule, which upon protonation yielded a species with significantly enhanced HRS intensity, was developed and its utility for the determination of the CMC of simple fatty acids was demonstrated.

Journal Article↗

PubMed2004

Expression and nucleotide sequence of an INS (3) P1 synthase gene associated with low-phytate kernels in maize (Zea mays L.).

Most of the phosphorus (P) in maize (Zea mays L.) kernels is in the form of phytic acid. A low phytic acid (lpa) maize mutant, lpa1-1, displays levels reduced by 66%. A goal of genetic breeding is development of low phytic acid feedstocks to improve P nutrition of nonruminant animals and reduce the adverse environmental impacts of excess P in manure. The genetic basis of the lpa1-1 mutation is not known, but previous genetic mapping has shown both the mutant phenotype and the Ins (3) P(1) synthase (MIPS) gene, which encodes the first enzyme, myo-inositol phosphate synthase, in the phytic acid biosynthetic pathway, map to the same chromosomal region in maize. Research was conducted to determine the expression of the MIPS gene in lpa1-1 and wild-type kernels with similar genetic backgrounds and to ascertain if variation in the MIPS coding sequence could be inferred to be the basis of the mutation. MIPS enzyme activity determined by TLC was reduced 2-3-fold in mutant kernels. RT-PCR-based experiments using primers specific for the 1S-MIPS sequence indicated gene expression was reduced 50-60% in the mutant. Sequence analysis of the MIPS genomic sequence revealed 10 exons and 9 introns that are identical in both mutant and wild-type developing kernels. These findings support an association between reduced MIPS gene activity and low phytic acid content, but additional research should examine the promoter, the 5'UTR, or transcriptional controlling elements of the MIPS gene to ascertain the possible presence of a genetic lesion in those regions.

Base Sequence↗

PubMed2001

Detection of low levels of the amorphous phase in crystalline pharmaceutical materials by thermally stimulated current spectrometry.

PURPOSE: To demonstrate the applicability of thermally stimulated current (TSC) spectrometry for the detection of low levels of the amorphous phase in crystalline pharmaceutical materials. METHODS: A crystalline drug substance was melt quenched to produce an amorphous material. Blends of the crystalline and amorphous phases in different ratios (from 75:25 to 99:01) were prepared by serial dilution. TSC studies were performed by applying an electric field at a temperature above the glass transition temperature (Tg) to orient the dipoles, rapidly cooling to 0 degrees C, short circuiting for 1 min, and scanning at 7 degrees C/min to measure the depolarization current. The temperature of the peak in the spectrum corresponds to the Tg of the amorphous phase. Modulated differential scanning calorimtery (DSC) studies were performed using three different test protocols (varying linear heating rate, modulation amplitude, and time period). Powder X-ray diffraction (XRD) studies were performed using a Siemens D500 diffractometer. RESULTS: The ability to detect the amorphous phase by powder XRD is beset with problems due to indirect inference, orientation effects, and instrument-related intensity variations. Even using a consistent sampling procedure and an internal standard, the XRD could quantify the amorphous phase at a level of 5%. In the conventional or modulated DSC, the amorphous phase manifests itself as a shift in the baseline. Using modulated DSC it was possible to detect the amorphous phase at a level of 5% when tested at a heating rate of 2 degrees C/min and an amplitude of +/-1.0 degrees C with a period of 30 s. The moisture sorption method appears to have a similar detection capability. In TSC scans, the glass transition event due to molecular/segmental mobility in the amorphous phase was manifested as a peak/shoulder on the low-temperature side of the depolarization peak of the crystalline phase. The amorphous phase was unambiguously detected at 2% with a lower detection limit of approximately 1%. CONCLUSIONS: On the basis of the results of this preliminary investigation, TSC appears to be capable of detecting the amorphous phase at as low as approximately 1% in crystalline pharmaceuticals, thus offering a much needed capability in discerning factors.

Antihypertensive Agents↗

PubMed2006

Mechanistic aspects of proton chain transfer in the green fluorescent protein. Part II. A comparison of minimal quantum chemical models.

In this paper we report the results of extensive quantum chemical reaction pathway calculations for the electronic ground state of several different cluster models that mimic the proton chain transfer path within the green fluorescent protein (GFP). Our principal objective is to establish the robustness with respect to variations in the model of our recent mechanistic inferences for the ground state proton chain transfer [S. Wang and S. C. Smith, J. Phys. Chem. B, 2006, 110, 5084]. Additionally, comparison of our ground state results with the excited state proton transfer (ESPT) study by Vendrell et al. [O. Vendrell, R. Gelabert, M. Moreno and J. M. Lluch, J. Am. Chem. Soc., 2006, 128, 3564] leads to the conclusion that the mechanism of proton chain transfer may be expected to be analogous in ground and excited states, principally because in both cases the loss of the chromophore's phenolic proton contributes strongly to the reaction coordinate only late in the reaction path.

Computer Simulation↗