Search PubMed⌕ Search

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.

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

At least 235 records · Page 13Linked to original sources

The effects of scale: variation in the APOA1/C3/A4/A5 gene cluster.

While there is considerable appeal to the idea of selecting a few SNPs to represent all, or much, of the DNA sequence variability in a local chromosomal region, it is also important to quantify what detail is lost in adopting such an approach. To address this issue, we compared high- and low-resolution depictions of sequence diversity for the same genomic region, the APOA1/C3/A4/A5 gene cluster on chromosome 11. First, extensive re-sequencing identified all nucleotide and sequence haplotype variation of the linked apolipoprotein genes in 72 individuals from three populations: African-Americans from Jackson, Miss., Europeans from North Karelia, Finland, and European-Americans from Rochester, Minn. We identified 124 SNPs in 17.7 kb and significant differences in variation among genes. APOC3 gene diversity was particularly distinctive at high resolution, showing large allele frequency differences ( F(ST) values >0.250) between Jackson and the other two samples, and divergent population-specific haplotype lineages. Next, we selected haplotype-tagging SNPs (htSNPs) for each gene, at a density of approximately one SNP per kb, using an algorithm suggested by Stram et al. (2003). The 17 htSNPs identified were then used to reconstruct low-resolution haplotypes, from which inferences about the structure of variation were also drawn. This comparison showed that while the htSNPs successfully tagged common haplotype variation, they also left much underlying sequence diversity undetected and failed, in some cases, to co-classify groups of closely related haplotypes. The implications of these findings for other haplotype-based descriptions of human variation are discussed.

Apolipoprotein C-III↗

Relationships of Nematodirus species and Nematodirus battus isolates (Nematoda: Trichostrongyloidea) based on nuclear ribosomal DNA sequences.

Nuclear ribosomal sequence data from the internal transcribed spacers (ITS-1 and ITS-2), 5.8S subunit, and regions of the 18S and 28S genes were used to investigate sequence diversity among geographic samples of Nematodirus battus, and to infer phylogenetic relationships among Nematodirus species. Phylogenetic analysis of these data yielded strong support for relationships among species, depicting Nematodirus helvetianus and Nematodirus spathiger as sister-taxa and a clade of these 2 species and Nematodirus filicollis. This tree is consistent with caprine bovids as ancestral hosts, with a subsequent host shift to Bovinae in N. helvetianus. Eleven of 14 N. battus sequences were unique, with 19 variable sites among sequences representing 5 geographic samples. The lowest number of variable nucleotide sites was observed in samples representing apparently recent introductions to the United States and Canada, which is consistent with a population bottleneck concomitant with translocation. Comparison of directly sequenced polymerase chain reaction products and clones revealed evidence for intraindividual variation at some of the sequence sites, and this pattern of variation and that within geographic samples indicates incomplete rDNA repeat homogenization within species. This pattern of variation is not conducive for inferring phylogenetic relationships among sequences representing N. battus or addressing the putative history of introduction.

Animals↗

Detecting clinical and balanced selection using spatial autocorrelation analysis under kin-structured migration.

Recently spatial autocorrelation has been employed to infer microevolutionary processes from patterns of genetic variation. In theory, different processes should show characteristic signature correlograms; e.g., clinal selection should produce correlograms decreasing from positive to negative autocorrelation, whereas uniform balanced selection should lead to no spatial autocorrelation. The ability of a statistical method such as spatial autocorrelation analysis to distinguish between these selective regimes or even to detect departures from neutrality is dependent on the strength of the evolutionary force and the population structure. Weak selection or migration will not be apparent against the expected background of stochastic noise. Moreover, the population structure may generate sufficient stochastic variation such that even strong evolutionary forces may fail to be detected. This study uses computer simulation to examine the effects of kin-structured migration and three different selective regimes on the shape of spatial correlograms to assess the ability of this technique to detect different microevolutionary processes. Genetic variation among 8 loci is simulated in a linear set of 25 artificial populations. Kin-structured stepping-stone migration among adjacent populations is modeled; directional, balanced, and clinal selection, as well as neutral loci are considered. These experiments show that strong selection produces correlograms of the predicted shape. However, with an anthropologically reasonable population structure, considerable stochastic variation among correlograms for different alleles may still exist. This suggests the need for caution in inferring genetic process from spatial patterns.

Bias↗

Evolution of microsatellite alleles in four species of mice (genus Apodemus).

Microsatellite length variation was investigated at a highly variable microsatellite locus in four species of Apodemus. Information obtained from microsatellite allele sequences was contrasted with allele sizes, which included 18 electromorphs. Additional analysis of a 400-bp unique sequence in the flanking region identified 26 different haplotype sequences or "true" alleles in the sample. Three molecular mechanisms, namely, (1) addition/deletion of repeats, (2) substitutions and indels in the flanking region, and (3) mutations interrupting the repeat, contributed to the generation of allelic variation. Size homoplasy can be inferred for alleles within populations, from different populations of the same species, and from different species. We propose that microsatellite flanking sequences may be informative markers for investigating mutation processes in microsatellite repeats as well as phylogenetic relationships among alleles, populations, and species.

Alleles↗

Estimation of mutation rates in cultured mammalian cells.

The factors that affect reliable estimations of mutation rates (mu) in cultured mammalian somatic cell populations by fluctuation analysis are studied experimentally and statistically. We analyze the differential effect of the final cell population size in each culture (Nt) and the number of parallel cultures (C) on the variation in the rate estimates (mu) inferred from the P0 method. The analysis can be made after the derivation of the variance of mu, which is a measure of variation of mu for a given combination of Nt and C in a number of repeat experiments. The variance of mu is inversely proportional to C and to the square of Nt . Nt determines the probability of occurrence of mutations in a cell culture. By influencing the size of P0, Nt also determines whether a rate estimate is obtainable from the experiment. Since Po is estimated from the fraction of cultures containing no mutation in a set of C cultures, C becomes a determining factor for the accuracy of mu. The rate estimated from P0 is biased, but the bias is in general 2 orders of magnitude smaller than mu. By the selection of an appropriate combination of Nt and C for the experiment, this bias can be reduced even further. Based on the notion of comparing two proportions, we propose a test statistic and have applied it to experimental results for a test of equality of mutation rates in different cell lines. This development places the comparison of mutation rates on a statistical basis.

Animals↗

Dynamic Trees for unsupervised segmentation and matching of image regions.

We present a probabilistic framework--namely, multiscale generative models known as Dynamic Trees (DT)--for unsupervised image segmentation and subsequent matching of segmented regions in a given set of images. Beyond these novel applications of DTs, we propose important additions for this modeling paradigm. First, we introduce a novel DT architecture, where multilayered observable data are incorporated at all scales of the model. Second, we derive a novel probabilistic inference algorithm for DTs--Structured Variational Approximation (SVA)--which explicitly accounts for the statistical dependence of node positions and model structure in the approximate posterior distribution, thereby relaxing poorly justified independence assumptions in previous work. Finally, we propose a similarity measure for matching dynamic-tree models, representing segmented image regions, across images. Our results for several data sets show that DTs are capable of capturing important component-subcomponent relationships among objects and their parts, and that DTs perform well in segmenting images into plausible pixel clusters. We demonstrate the significantly improved properties of the SVA algorithm--both in terms of substantially faster convergence rates and larger approximate posteriors for the inferred models--when compared with competing inference algorithms. Furthermore, results on unsupervised object recognition demonstrate the viability of the proposed similarity measure for matching dynamic-structure statistical models.

Algorithms↗

Evolutionary relationships among the eukaryotic crown taxa taking into account site-to-site rate variation in 18S rRNA.

In this study we constructed a bootstrapped distance tree of 500 small subunit ribosomal RNA sequences from organisms belonging to the so-called crown of eukaryote evolution. Taking into account the substitution rate of the individual nucleotides of the rRNA sequence alignment, our results suggest that (1) animals, true fungi, and choanoflagellates share a common origin: The branch joining these taxa is highly supported by bootstrap analysis (bootstrap support [BS] > 90%), (2) stramenopiles and alveolates are sister groups (BS = 75%), (3) within the alveolates, dinoflagellates and apicomplexans share a common ancestor BS > 95%), while in turn they both share a common origin with the ciliates (BS > 80%), and (4) within the stramenopiles, heterokont algae, hyphochytriomycetes, and oomycetes form a monophyletic grouping well supported by bootstrap analysis (BS > 85%), preceded by the well-supported successive divergence of labyrinthulomycetes and bicosoecids. On the other hand, many evolutionary relationships between crown taxa are still obscure on the basis of 18S rRNA. The branching order between the animal-fungal-choanoflagellates clade and the chlorobionts, the alveolates and stramenopiles, red algae, and several smaller groups of organisms remains largely unresolved.When among-site rate variation is not considered, the inferred tree topologies are inferior to those where the substitution rate spectrum for the 18S rRNA is taken into account. This is primarily indicated by the erroneous branching of fast-evolving sequences. Moreover, when different substitution rates among sites are not considered, the animals no longer appear as a monophyletic grouping in most distance trees.

Animals↗

Genetic variation of Ginkgo biloba L. (Ginkgoaceae) based on cpDNA PCR-RFLPs: inference of glacial refugia.

Ginkgo biloba, a famous living fossil, is the sole survivor of the genus Ginkgo. To make inferences about the glacial refugia that harbored G. biloba, we examined the genetic structure of eight potential refugial populations and plantations using chloroplast DNA (cpDNA) with eight size variants in the trnK1-trnK2 fragment. The data consist of haplotypes from 158 trees collected from eight localities. The majority of the cpDNA haplotypes are restricted to minor portions of the geographical range. Our results suggest that refugia of G. biloba were located in southwestern China. This area is a current biodiversity hotspot of global importance, and may have been protected from the extremes of climatic fluctuations during the Pleistocene. The Ginkgos on West Tianmu Mountain, which were previously considered to be wild by many researchers, may, instead, have been introduced by Buddhist monks.

China↗

Geographic structure of mitochondrial and nuclear gene polymorphisms in Australian green turtle populations and male-biased gene flow.

The genetic structure of green turtle (Chelonia mydas) rookeries located around the Australian coast was assessed by (1) comparing the structure found within and among geographic regions, (2) comparing microsatellite loci vs. restriction fragment length polymorphism analyses of anonymous single copy nuclear DNA (ascnDNA) loci, and (3) comparing the structure found at nuclear DNA markers to that of previously analyzed mitochondrial (mtDNA) control region sequences. Significant genetic structure was observed over all regions at both sets of nuclear markers, though the microsatellite data provided greater resolution in identifying significant genetic differences in pairwise tests between regions. Inferences about population structure and migration rates from the microsatellite data varied depending on whether statistics were based on the stepwise mutation or infinite allele model, with the latter being more congruent with geography. Estimated rates of gene flow were generally higher than expected for nuclear DNA (nDNA) in comparison to mtDNA, and this difference was most pronounced in comparisons between the northern and southern Great Barrier Reef (GBR). The genetic data combined with results from physical tagging studies indicate that the lack of nuclear gene divergence through the GBR is likely due to the migration of sGBR turtles through the courtship area of the nGBR population, rather than male-biased dispersal. This example highlights the value of combining comparative studies of molecular variation with ecological data to infer population processes.

Alleles↗

Genetic diversity among a complex of cereal cyst nematodes inferred from RFLP analysis of the ribosomal internal transcribed spacer region.

This study examined the restriction polymorphism (RFLP) of the nuclear ribosomal DNA in Heterodera avenae, H. filipjevi, H. mani, H. latipons, and the taxonomically unclear Gotland strain in order to establish a molecular characterization and phylogenetic relationships in the complex of cereal cyst nematodes (CCN). The internal transcribed spacer (ITS) and 5.8S rDNA were amplified by PCR from a single female or a cyst of 27 different geographic isolates of the CCN complex and one population of H. schachtii, used as outgroup. The amplified product was 1.2 kb long and 14 of 15 enzymes produced restriction fragments for each isolate. Relationships between populations were determined from UPGMA analysis based on distance values calculated from RFLP data. Digestions with TaqI clearly differentiated H. avenae, H. latipons, and a group composed of H. filipjevi and the Gotland strain. Six endonucleases (HaeIII, HinfI, ItaI, PstI, TaqI, and Tru9I) produced the same restriction pattern with H. filipjevi and the Gotland strain, and both were clearly separated from H. avenae with PstI. Restriction sites have revealed a mixture of the species H. latipons and H. avenae, and possible infraspecific variation in H. avenae. The inferred phylogenetic relationships of species in the CCN complex are in agreement with their morphological characterization.

Animals↗

Spatial ecology of mating success in a sexually polymorphic plant.

The spatial context of reproduction is of crucial importance to plants because of their sessile habit. Since pollen and seed dispersal is often restricted, mating success is likely to depend on the quantity and quality of mates in local neighbourhoods. Here we use neighbourhood models to investigate the spatial ecology of pollination and mating in Narcissus assoanus, a sexually polymorphic plant with two mating morphs that differ in style length. By mapping individuals in eight populations from southwestern France, we investigated the influence of the density and morph identity of plants at different spatial scales on variation in female fertility. By using inferences on the expected patterns of pollen transfer based on floral morphology, we were able to predict the quantitative relations between local morph ratios and variation in fertility. Our analyses revealed differences in the spatial clustering of morphs and in their response to plant density and morph identity within local neighbourhoods. Mating success in N. assoanus was characterized by both density- and frequency-dependent processes, a condition that may be a general feature of the spatial ecology of plant mating.

Ecosystem↗

Phylogeography of Mabuya maculilabris (Reptilia) from São Tomé Island (Gulf of Guinea) inferred from mtDNA sequences.

The pattern of genetic variation of the lizard Mabuya maculilabris from São Tomé Island (Gulf of Guinea) was investigated using a combination of three mitochondrial DNA gene fragments. Forty-eight haplotypes were recovered among 66 individuals covering the whole island. The genealogy inferred from the most parsimonious network of haplotypes allows us to detect two main and long branches departing from the putative group of oldest haplotypes. The tips of these branches exhibit star-like phylogenies, which may indicate of recently expanded populations, most probably from a small number of founders. A nested clade analysis suggests a complex pattern of past events that gave rise to the extant geographical pattern found in the haplotype distribution: past and allopatric fragmentation, range expansion, restricted gene flow and long-distance dispersal. These results are consistent with the complex geological history of the island where important volcanic activity with extensive lava flows has occurred during several periods. Mismatch-distribution analysis and AMOVA also support these conclusions. Substantial genetic structuring among these lizards was detected as well as high levels of differentiation between the southern edge populations (particularly those from the Rolas Islet) and the remaining ones. However, variation is low relative to the geological age of the island. Our results indicate that patterns of variation observed in reptiles in other oceanic islands are not indicative of those observed in the islands of the Gulf of Guinea.

Animals↗

Investigating deep phylogenetic relationships among cyanobacteria and plastids by small subunit rRNA sequence analysis.

Small subunit rRNA sequence data were generated for 27 strains of cyanobacteria and incorporated into a phylogenetic analysis of 1,377 aligned sequence positions from a diverse sampling of 53 cyanobacteria and 10 photosynthetic plastids. Tree inference was carried out using a maximum likelihood method with correction for site-to-site variation in evolutionary rate. Confidence in the inferred phylogenetic relationships was determined by construction of a majority-rule consensus tree based on alternative topologies not considered to be statistically significantly different from the optimal tree. The results are in agreement with earlier studies in the assignment of individual taxa to specific sequence groups. Several relationships not previously noted among sequence groups are indicated, whereas other relationships previously supported are contradicted. All plastids cluster as a strongly supported monophyletic group arising near the root of the cyanobacterial line of descent.

Cyanobacteria↗

Mammalian male mutation bias: impacts of generation time and regional variation in substitution rates.

In mammals, males undergo a greater number of germline cell divisions compared with females. Thus, the male germline accumulates more DNA replication errors, which result in male mutation bias-a higher mutation rate for males than for females. The phenomenon of male mutation bias has been investigated mostly for rodents and primates, however, it has not been studied in detail for other mammalian orders. Here we sequenced and analyzed five introns of three genes (DBX/DBY, UTX/UTY, and ZFX/ZFY) homologous between X and Y chromosomes in several species of perissodactyls (horses and rhinos) and of primates. Male mutation bias was evident: substitution rate was higher for a Y chromosome intron than for its X chromosome homologue for all five intron pairs studied. Substitution rates varied regionally among introns sequenced on the same chromosome and this variation influenced male mutation bias inferred from each intron pair. Interestingly, we observed a positive correlation in substitution rates between homologous X and homologous Y introns as well as between orthologous primate and perissodactyl introns. The male-to-female mutation rate ratio estimated from concatenated sequences of five perissodactyl introns was 3.88 (95% CI = 2.90-6.07). Using the data generated here and estimates available in the literature, we compared male mutation bias among several mammalian orders. We conclude that male mutation bias is significantly higher for organisms with long generation times (primates, perissodactyls, and felids) than for organisms with short generation times (e.g., rodents) since the former undergo a greater number of male germline cell divisions.

Aging↗

Individual variation in measures of voice.

Measures of inferred subglottal air pressure, glottal airflow waveform characteristics, sound pressure level (SPL) and the acoustic spectral slope were studied for individual speakers with normal voices. Combinations of different levels of subglottal air pressure and varying glottal configurations could result in the same SPL. Relatively high air pressure levels were associated with a steep spectral slope, reflecting a more sinusoidal glottal waveform and a relatively abducted membranous glottis, which would result in damping of F1. Data suggested that the interarytenoid glottal opening could vary without systematically affecting SPL or voice quality. The results indicate that the principles of production-related economy of effort and physiological, acoustic and perceptual constraints may apply to voice production.

Air Pressure↗

Inference of demographic history from genealogical trees using reversible jump Markov chain Monte Carlo.

BACKGROUND: Coalescent theory is a general framework to model genetic variation in a population. Specifically, it allows inference about population parameters from sampled DNA sequences. However, most currently employed variants of coalescent theory only consider very simple demographic scenarios of population size changes, such as exponential growth. RESULTS: Here we develop a coalescent approach that allows Bayesian non-parametric estimation of the demographic history using genealogies reconstructed from sampled DNA sequences. In this framework inference and model selection is done using reversible jump Markov chain Monte Carlo (MCMC). This method is computationally efficient and overcomes the limitations of related non-parametric approaches such as the skyline plot. We validate the approach using simulated data. Subsequently, we reanalyze HIV-1 sequence data from Central Africa and Hepatitis C virus (HCV) data from Egypt. CONCLUSIONS: The new method provides a Bayesian procedure for non-parametric estimation of the demographic history. By construction it additionally provides confidence limits and may be used jointly with other MCMC-based coalescent approaches.

Algorithms↗

Selective neutrality and enzyme kinetics.

This article appeals to a recent theory of enzyme evolution to show that the properties, neutral or adaptive, which characterize the observed allelic variation in natural populations can be inferred from the functional parameters, substrate specificity, and reaction rate. This study delineates the following relations between activity variables, and the forces--adaptive or neutral--determining allelic variation: (1) Enzymes with broad substrate specificity: The observed polymorphism is adaptive; mutations in this class of enzymes can result in increased fitness of the organism and hence be relevant for positive selection. (2) Enzymes with absolute substrate specificity and diffusion-controlled rates: Observed allelic variation will be absolutely neutral; mutations in this class of enzymes will be either deleterious or have no effect on fitness. (3) Enzymes with absolute or group specificity and nondiffusion-controlled rates: Observed variation will be partially neutral; mutants which are selectively neutral may become advantageous under an appropriate environmental condition or different genetic background. We illustrate each of the relations between kinetic properties and evolutionary states with examples drawn from enzymes whose evolutionary dynamics have been intensively studied.

Adaptation, Physiological↗

Genetic variation of FUT2 in a Ghanaian population: identification of four novel mutations and inference of balancing selection.

The FUT2 is rich in polymorphisms that show population-specific patterns in its coding sequence. In this study, we found four novel substitutions in 121 Ghanaian samples. In addition, statistical tests considering population expansion scenarios suggested that balancing selection might be responsible for the genetic diversity of FUT2 in this population.

Base Sequence↗