Pitfalls of the most commonly used models of context dependent substitution.
Ontology highlight
ABSTRACT:
SUBMITTER: Lindsay H
PROVIDER: S-EPMC2662811 | biostudies-other | 2009
REPOSITORIES: biostudies-other
Similar Datasets
Project description:BackgroundNeighboring nucleotides exert a striking influence on mutation, with the hypermutability of CpG dinucleotides in many genomes being an exemplar. Among the approaches employed to measure the relative importance of sequence neighbors on molecular evolution have been continuous-time Markov process models for substitutions that treat sequences as a series of independent tuples. The most widely used examples are the codon substitution models. We evaluated the suitability of derivatives of the nucleotide frequency weighted (hereafter NF) and tuple frequency weighted (hereafter TF) models for measuring sequence context dependent substitution. Critical properties we address are their relationships to an independent nucleotide process and the robustness of parameter estimation to changes in sequence composition. We then consider the impact on inference concerning dinucleotide substitution processes from application of these two forms to intron sequence alignments from primates.ResultsWe prove that the NF form always nests the independent nucleotide process and that this is not true for the TF form. As a consequence, using TF to study context effects can be misleading, which is shown by both theoretical calculations and simulations. We describe a simple example where a context parameter estimated under TF is confounded with composition terms unless all sequence states are equi-frequent. We illustrate this for the dinucleotide case by simulation under a nucleotide model, showing that the TF form identifies a CpG effect when none exists. Our analysis of primate introns revealed that the effect of nucleotide neighbors is over-estimated under TF compared with NF. Parameter estimates for a number of contexts are also strikingly discordant between the two model forms.ConclusionOur results establish that the NF form should be used for analysis of independent-tuple context dependent processes. Although neighboring effects in general are still important, prominent influences such as the elevated CpG transversion rate previously identified using the TF form are an artifact. Our results further suggest as few as 5 parameters may account for approximately 85% of neighboring nucleotide influence.
| S-EPMC2628887 | biostudies-literature
Project description:MotivationThe modeling of conservation patterns in genomic DNA has become increasingly popular for a number of bioinformatic applications. While several systems developed to date incorporate context-dependence in their substitution models, the impact on computational complexity and generalization ability of the resulting higher order models invites the question of whether simpler approaches to context modeling might permit appreciable reductions in model complexity and computational cost, without sacrificing prediction accuracy.ResultsWe formulate several alternative methods for context modeling based on windowed Bayesian networks, and compare their effects on both accuracy and computational complexity for the task of discriminating functionally distinct segments in vertebrate DNA. Our results show that substantial reductions in the complexity of both the model and the associated inference algorithm can be achieved without reducing predictive accuracy.
| S-EPMC2732293 | biostudies-literature
Project description:The mouse is a useful preclinical species for evaluating disease etiology due to the availability of a wide variety of genetically modified strains and the ability to perform disease-modifying manipulations. In order to establish an atrial filtration (AF) model in our laboratory, we profiled several commonly used murine AF models. We initially evaluated a pharmacological model of acute carbachol (CCh) treatment plus atrial burst pacing in C57BL/6 mice. In an effort to observe micro-reentrant circuits indicative of authentic AF, we employed optical mapping imaging in isolated mouse hearts. While CCh reduced atrial refractoriness and increased atrial tachyarrhythmia vulnerability, the left atrial (LA) excitation patterns were rather regular without reentrant circuits or wavelets. Therefore, the atrial tachyarrhythmia resembled high frequency atrial flutter, not typical AF per se. We next examined both a chronic angiotensin II (Ang II) infusion model and the surgical model of transverse aortic constriction (TAC), which have both been reported to induce atrial and ventricular structural changes that serve as a substrates for micro-reentrant AF. Although we observed some extent of atrial remodeling such as fibrosis or enlarged LA diameter, burst pacing-induced atrial tachyarrhythmia vulnerability did not differ from control mice in either model. This again suggested that an AF-like pathophysiology is difficult to demonstrate in the mouse. To continue searching for a valid murine AF model, we studied mice with a cardiac-specific deficiency (KO) in liver kinase B1 (Cardiac-LKB1), which has been reported to exhibit spontaneous AF. Indeed, the electrocardiograms (ECG) of conscious Cardiac-LKB1 KO mice exhibited no P waves and had irregular RR intervals, which are characteristics of AF. Histological evaluation of Cardiac-LKB1 KO mice revealed dilated and fibrotic atria, again consistent with AF. However, atrial electrograms and optical mapping revealed that electrical activity was limited to the sino-atrial node area with no electrical conduction into the atrial myocardium beyond. Thus, Cardiac-LKB1 KO mice have severe atrial myopathy or atrial standstill, but not AF. In summary, the atrial tachyarrhythmias we observed in the four murine models were distinct from typical human AF, which often exhibits micro- or macro-reentrant atrial circuits. Our results suggest that the four murine AF models we examined may not reflect human AF well, and raise a cautionary note for use of those murine models to study AF.
| S-EPMC8741011 | biostudies-literature
Project description:BackgroundModels of sequence evolution typically assume that different nucleotide positions evolve independently. This assumption is widely appreciated to be an over-simplification. The best known violations involve biases due to adjacent nucleotides. There have also been suggestions that biases exist at larger scales, however this possibility has not been systematically explored.ResultsTo address this we have developed a method which identifies over- and under-represented substitution patterns and assesses their overall impact on the evolution of genome composition. Our method is designed to account for biases at smaller pattern sizes, removing their effects. We used this method to investigate context bias in the human lineage after the divergence from chimpanzee. We examined bias effects in substitution patterns between 2 and 5 bp long and found significant effects at all sizes. This included some individual three and four base pair patterns with relatively large biases. We also found that bias effects vary across the genome, differing between transposons and non-transposons, between different classes of transposons, and also near and far from genes.ConclusionsWe found that nucleotides beyond the immediately adjacent one are responsible for substantial context effects, and that these biases vary across the genome.
| S-EPMC2945941 | biostudies-literature
Project description:This SuperSeries is composed of the SubSeries listed below.
2022-07-31 | GSE197330 | GEO
Project description:IntroductionUltrasonic aspiration (UA) devices are commonly used for resecting intracranial tumors, as they allow for internal debulking of large tumors, hereby avoiding damage to adjacent brain tissue during the dissection. Little is known about their comparative safety profiles.Methods and materialsWe analyzed data from a prospective patient registry. Procedures using one of the following UA models were included: Integra® CUSA, Söring®, and Stryker® Sonopet. The primary endpoint was morbidity at discharge, defined as significant worsening on the Karnofsky Performance Scale. Secondary endpoints included morbidity and mortality until 3 months postoperative (M3), occurrence, type, and etiology of complications.ResultsOf n = 1028 procedures, the CUSA was used in n = 354 (34.4 %), the Söring in n = 461 (44.8 %), and the Sonopet in n = 213 (20.7 %). There was some heterogeneity of study groups. In multivariable analysis, patients in the Söring (adjusted odds ratio (aOR) 1.29; 95 % confidence interval (CI), 0.80-2.08; p = 0.299), and Sonopet group (aOR, 0.86; 95 % CI, 0.46-1.61; p = 0.645) were as likely as patients in the CUSA group to experience discharge morbidity. At M3, patients in the Söring (aOR, 1.20; 95 % CI, 0.78-1.86; p = 0.415) and Sonopet group (aOR, 0.53; 95 % CI, 0.26-1.08; p = 0.080) were as likely as patients in the CUSA group to experience morbidity. There were also no differences for M3 morbidity in subgroup analyses for gliomas, meningiomas, and metastases. The grade (p = 0.608) and etiology (p = 0.849) of postoperative complications were similar.ConclusionsNeurosurgeons select UA types with regard to certain case-specific characteristics. The safety profiles of three commonly used UA types appear mostly similar.
| S-EPMC6739453 | biostudies-literature
Project description:Background: Despite advances in therapeutics, outcomes for hepatocellular carcinoma (HCC) remain poor and there is an urgent need for efficacious systemic therapy. Unfortunately, drugs that are successful in preclinical studies often fail in the clinical setting, and we hypothesize that this is due to functional differences between primary tumors and commonly used preclinical models. In this study, we attempt to answer this question by comparing tumor morphology and gene expression profiles between primary tumors, xenografts and HCC cell lines. Methods: Hep G2 cell lines and tumor cells from patient tumor explants were subcutaneously (ectopically) injected into the flank and orthotopically into liver parenchyma of Mus Musculus SCID mice. The mice were euthanized after two weeks. RNA was extracted from the tumors, and gene expression profiling was performed using the Gene Chip Human Genome U133 Plus 2.0. Principal component analyses (PCA) and construction of dendrograms were conducted using Partek genomics suite. Results: PCA showed that the commonly used HepG2 cell line model and its xenograft counterparts were vastly different from all fresh, primary tumors. Expression profiles of primary tumors were also significantly divergent from their counterpart patient-derived xenograft (PDX) models, regardless of the site of implantation. Xenografts from the same primary tumors were more likely to cluster together regardless of site of implantation, although heat maps showed distinct differences in gene expression profiles between orthotopic and ectopic models. Conclusions: The data presented here challenges the utility of routinely used preclinical models. Models using HepG2 were vastly different from primary tumors and PDXs, suggesting that this is not clinically representative. Surprisingly, site of implantation (orthotopic versus ectopic) resulted in limited impact on gene expression profiles, and in both scenarios xenografts differed significantly from the original primary tumors, challenging the long-held notion that orthotopic PDX model is the gold standard preclinical model for HCC.
2015-11-17 | GSE72981 | GEO
Project description:MYC is one of the most dysregulated oncogenes and is thought to be fundamental to tumor formation and/or maintenance in many cancer types. This dominant pro-tumor activity makes MYC an attractive target for cancer therapy. However, MYC is a transcription factor lacking enzymatic activity, and the structure of one of its two domains is unknown e.g., its transactivation domain. Consequently, few direct MYC-targeting therapies have been developed, and none have been successful in the clinic. Nevertheless, significant effort has been devoted to understanding the mechanisms of oncogenic MYC activity with the objective of uncovering novel vulnerabilities of MYC-dependent cancers. These extensive investigations have revealed in detail how MYC translocation, amplification, and other upstream perturbations contribute to MYC activity in cancer. However, missense mutations of the MYC gene have remained relatively understudied for their potential role in MYC-mediated oncogenesis. While the function of several low-frequency mutations in MYC have been described, our understanding of other equally or more frequent mutations is incomplete. Herein, we define the function of a recurrent missense mutation in MYC resulting in the substitution S146L. This mutation enhances the interaction between MYC and its cofactor TRRAP and may enhance oncogenic MYC activity in certain cellular contexts.
| S-EPMC9417018 | biostudies-literature
Project description:BackgroundEcological communities tend to be spatially structured due to environmental gradients and/or spatially contagious processes such as growth, dispersion and species interactions. Data transformation followed by usage of algorithms such as Redundancy Analysis (RDA) is a fairly common approach in studies searching for spatial structure in ecological communities, despite recent suggestions advocating the use of Generalized Linear Models (GLMs). Here, we compared the performance of GLMs and RDA in describing spatial structure in ecological community composition data. We simulated realistic presence/absence data typical of many β-diversity studies. For model selection we used standard methods commonly used in most studies involving RDA and GLMs.MethodsWe simulated communities with known spatial structure, based on three real spatial community presence/absence datasets (one terrestrial, one marine and one freshwater). We used spatial eigenvectors as explanatory variables. We varied the number of non-zero coefficients of the spatial variables, and the spatial scales with which these coefficients were associated and then compared the performance of GLMs and RDA frameworks to correctly retrieve the spatial patterns contained in the simulated communities. We used two different methods for model selection, Forward Selection (FW) for RDA and the Akaike Information Criterion (AIC) for GLMs. The performance of each method was assessed by scoring overall accuracy as the proportion of variables whose inclusion/exclusion status was correct, and by distinguishing which kind of error was observed for each method. We also assessed whether errors in variable selection could affect the interpretation of spatial structure.ResultsOverall GLM with AIC-based model selection (GLM/AIC) performed better than RDA/FW in selecting spatial explanatory variables, although under some simulations the methods performed similarly. In general, RDA/FW performed unpredictably, often retaining too many explanatory variables and selecting variables associated with incorrect spatial scales. The spatial scale of the pattern had a negligible effect on GLM/AIC performance but consistently affected RDA's error rates under almost all scenarios.ConclusionWe encourage the use of GLM/AIC for studies searching for spatial drivers of species presence/absence patterns, since this framework outperformed RDA/FW in situations most likely to be found in natural communities. It is likely that such recommendations might extend to other types of explanatory variables.
| S-EPMC7474884 | biostudies-literature
Project description:BackgroundDespite advances in therapeutics, outcomes for hepatocellular carcinoma (HCC) remain poor and there is an urgent need for efficacious systemic therapy. Unfortunately, drugs that are successful in preclinical studies often fail in the clinical setting, and we hypothesize that this is due to functional differences between primary tumors and commonly used preclinical models. In this study, we attempt to answer this question by comparing tumor morphology and gene expression profiles between primary tumors, xenografts and HCC cell lines.MethodsHep G2 cell lines and tumor cells from patient tumor explants were subcutaneously (ectopically) injected into the flank and orthotopically into liver parenchyma of Mus Musculus SCID mice. The mice were euthanized after two weeks. RNA was extracted from the tumors, and gene expression profiling was performed using the Gene Chip Human Genome U133 Plus 2.0. Principal component analyses (PCA) and construction of dendrograms were conducted using Partek genomics suite.ResultsPCA showed that the commonly used HepG2 cell line model and its xenograft counterparts were vastly different from all fresh primary tumors. Expression profiles of primary tumors were also significantly divergent from their counterpart patient-derived xenograft (PDX) models, regardless of the site of implantation. Xenografts from the same primary tumors were more likely to cluster together regardless of site of implantation, although heat maps showed distinct differences in gene expression profiles between orthotopic and ectopic models.ConclusionsThe data presented here challenges the utility of routinely used preclinical models. Models using HepG2 were vastly different from primary tumors and PDXs, suggesting that this is not clinically representative. Surprisingly, site of implantation (orthotopic versus ectopic) resulted in limited impact on gene expression profiles, and in both scenarios xenografts differed significantly from the original primary tumors, challenging the long-held notion that orthotopic PDX model is the gold standard preclinical model for HCC.
| S-EPMC4628260 | biostudies-literature