Finding Quantitative Trait Loci Genes with Collaborative Targeted Maximum Likelihood Learning.
Ontology highlight
ABSTRACT:
SUBMITTER: Wang H
PROVIDER: S-EPMC3090625 | biostudies-literature | 2011 Jul
REPOSITORIES: biostudies-literature
Publications
Finding Quantitative Trait Loci Genes with Collaborative Targeted Maximum Likelihood Learning.
Wang Hui H Rose Sherri S van der Laan Mark J MJ
Statistics & probability letters 20110701 7
Quantitative trait loci mapping is focused on identifying the positions and effect of genes underlying an an observed trait. We present a collaborative targeted maximum likelihood estimator in a semi-parametric model using a newly proposed 2-part super learning algorithm to find quantitative trait loci genes in listeria data. Results are compared to the parametric composite interval mapping approach. ...[more]
Similar Datasets
Project description:We present a novel semiparametric method for quantitative trait loci (QTL) mapping in experimental crosses. Conventional genetic mapping methods typically assume parametric models with Gaussian errors and obtain parameter estimates through maximum-likelihood estimation. In contrast with univariate regression and interval-mapping methods, our model requires fewer assumptions and also accommodates various machine-learning algorithms. Estimation is performed with targeted maximum-likelihood learning methods. We demonstrate our semiparametric targeted learning approach in a simulation study and a well-studied barley data set.
| S-EPMC4256757 | biostudies-literature
Project description:Evolutionary processes during plant polyploidization and speciation have led to extensive presence-absence variation (PAV) in crop genomes, and there is increasing evidence that PAV associates with important traits. Today, high-resolution genetic analysis in major crops frequently implements simple, cost-effective, high-throughput genotyping from single nucleotide polymorphism (SNP) hybridization arrays; however, these are normally not designed to distinguish PAV from failed SNP calls caused by hybridization artefacts. Here, we describe a strategy to recover valuable information from single nucleotide absence polymorphisms (SNaPs) by population-based quality filtering of SNP hybridization data to distinguish patterns associated with genuine deletions from those caused by technical failures. We reveal that including SNaPs in genetic analyses elucidate segregation of small to large-scale structural variants in nested association mapping populations of oilseed rape (Brassica napus), a recent polyploid crop with widespread structural variation. Including SNaP markers in genomewide association studies identified numerous quantitative trait loci, invisible using SNP markers alone, for resistance to two major fungal diseases of oilseed rape, Sclerotinia stem rot and blackleg disease. Our results indicate that PAV has a strong influence on quantitative disease resistance in B. napus and that SNaP analysis using cost-effective SNP array data can provide extensive added value from 'missing data'. This strategy might also be applicable for improving the precision of genetic mapping in many important crop species.
| S-EPMC6230954 | biostudies-literature
Project description:Sparse modeling, a feature selection method widely used in the machine-learning community, has been recently applied to identify associations in genetic studies including expression quantitative trait locus (eQTL) mapping. These genetic studies usually involve high dimensional data where the number of features is much larger than the number of samples. The high dimensionality of genetic data introduces a problem that there exist multiple solutions for optimizing a sparse model. In such situations, a single optimization result provides only an incomplete view of the data and lacks power to find alternative features associated with the same trait. In this article, we propose a novel method aimed to detecting alternative eQTLs where two genetic variants have alternative relationships regarding their associations with the expression of a particular gene. Our method accomplishes this goal by exploring multiple solutions sampled from the solution space. We proved our method theoretically and demonstrated its usage on simulated data. We then applied our method to a real eQTL data and identified a set of alternative eQTLs with potential biological insights. Additionally, these alternative eQTLs implicate a network view of understanding gene regulation.
| S-EPMC4010169 | biostudies-literature
Project description:Trauma- and stress-related disorders are clinically heterogeneous and associated with substantial genetic risk. Understanding the biological origins of heterogeneity of key intermediate phenotypes such as cognition and emotion can provide novel mechanistic insights into disorder pathogenesis. Performing quantitative genetics in animal models is a tractable strategy for examining both the genetic basis of intermediate phenotypes and functional testing of candidate quantitative traits genes (QTGs). Here, existing and newly collected data were used for collaborative genome-wide mapping of cued fear acquisition and expression in 65 mouse strains from the BXD genetic reference panel. For fear acquisition, we identified a significant locus on chromosome (Chr) 10 and eight suggestive loci on Chr 2, 4, 5, 11, 13, and 15. For fear expression, we identified one significant and another highly suggestive locus on Chr 13, as well as four suggestive loci on Chr 10, 11, and X. Across these loci, 60 putative QTGs were identified. The quantitative trait locus on distal Chr 13 contained a single, highly promising gene at the location of the peak likelihood ratio statistic score. The gene, hyperpolarization-activated cyclic nucleotide-gated channel 1 (Hcn1), regulates neuronal excitability. Validation experiments using behavioral pharmacology revealed that functional Hcn channels in the basolateral amygdala are necessary for conditioned fear acquisition and expression. Hcn1, together with the other candidate QTGs, thus provide new targets for neurobiological and treatment studies of fear learning and trauma- and stress-related disorders.There is a knowledge gap in understanding the genetic contributions to behavioral heterogeneity in typical and atypical populations. Mouse genetic reference panels (GRPs) provide one approach for identifying genetic sources of variation. Here, we identified three loci for conditioned fear acquisition and expression in a mouse GRP. Each locus contained candidate quantitative trait genes (QTGs). One locus had a single QTG, Hcn1 (hyperpolarization-activated cyclic nucleotide-gated channel 1), which has been implicated in neuronal excitability and learning. This discovery was validated using behavioral pharmacology, revealing that Hcn channels in the basolateral amygdala are required for fear acquisition and expression. The study thus identifies novel candidate QTGs that may contribute to variation in emotional learning and highlight the utility of mouse GRPs for the identification of genes underlying complex traits.
| S-EPMC4899527 | biostudies-other
Project description:Estimators of the conditional expectation, i.e., prediction, function involve a global bias-variance trade off. In some cases, an estimator that yields unbiased estimates of the conditional expectation for a particular partitioning of the data may be desirable. Such estimators are calibrated with respect to the partitioning. We identify the conditional expectation given a particular partitioning as a smooth parameter of the distribution of the data, where the partitioning may be defined on the covariate space or on the prediction space of the estimator. We propose a targeted maximum likelihood estimation (TMLE) procedure that updates an initial prediction estimator such that the updated estimator yields an unbiased and efficient estimator of this smooth parameter in the nonparametric statistical model. When the partitioning is defined on the prediction space of the estimator, our TMLE involves enforcing an implicit constraint on the estimator itself. We show that our resulting estimator of the smooth parameter is equal to the empirical estimator, which is also known to be unbiased and efficient in the nonparametric statistical model. We derive the TMLE for single time-point prediction and also time-dependent prediction in a counting process framework.
| S-EPMC6053284 | biostudies-literature
Project description:ObjectivesTo compare the performance of a targeted maximum likelihood estimator (TMLE) and a collaborative TMLE (CTMLE) to other estimators in a drug safety analysis, including a regression-based estimator, propensity score (PS)-based estimators, and an alternate doubly robust (DR) estimator in a real example and simulations.Study design and settingThe real data set is a subset of observational data from Kaiser Permanente Northern California formatted for use in active drug safety surveillance. Both the real and simulated data sets include potential confounders, a treatment variable indicating use of one of two antidiabetic treatments and an outcome variable indicating occurrence of an acute myocardial infarction (AMI).ResultsIn the real data example, there is no difference in AMI rates between treatments. In simulations, the double robustness property is demonstrated: DR estimators are consistent if either the initial outcome regression or PS estimator is consistent, whereas other estimators are inconsistent if the initial estimator is not consistent. In simulations with near-positivity violations, CTMLE performs well relative to other estimators by adaptively estimating the PS.ConclusionEach of the DR estimators was consistent, and TMLE and CTMLE had the smallest mean squared error in simulations.
| S-EPMC3818128 | biostudies-literature
Project description:A genome-wide association study (GWAS) was conducted to examine expression quantitative trait loci (eQTLs) for histone genes. We examined common eQTLs for multiple histone genes in 373 European lymphoblastoid cell lines (LCLs). A linear regression model was employed to identify single-nucleotide polymorphisms (SNPs) associated with expression of the histone genes, and the number of eQTLs was determined by linkage disequilibrium analysis. Additional associations of the identified eQTLs with other genes were also examined. We identified 31 eQTLs for 29 histone genes through genome-wide analysis using 29 histone genes (P < 2.97 × 10-10). Among them, 12 eQTLs were associated with the expression of multiple histone genes. Transcriptome-wide association analysis using the identified eQTLs showed their associations with additional 80 genes (P < 4.75 × 10-6). In particular, expression of RPPH1, SCARNA2, and SCARNA7 genes was associated with 26, 25, and 23 eQTLs, respectively. This study suggests that histone genes shared 12 common eQTLs that might regulate cell cycle-dependent transcription of histone and other genes. Further investigations are needed to elucidate the transcriptional mechanisms of these genes.
| S-EPMC5591967 | biostudies-literature
Project description:DNA methylation is a reaction that results in the formation of 5-methylcytosine when a methyl group is added to the cytosine’s C5 position. As organisms age, DNA methylation patterns change in a reproducible fashion. This phenomenon has established DNA methylation as a valuable biomarker in aging studies. Epigenetic clocks based on weighted combinations of methylation sites have been developed to accurately predict the age of an individual from their methylome. However, many epigenetic clocks, particularly those that utilize penalized regression, model the changes in methylation linearly with age. Moreover, these models, which use methylation levels as features, are not robust to missing data and do not account for the count-based nature of bisulfite sequence data. Additionally, the models are generally not interpretable. To overcome these challenges, we present BayesAge, an extension of the previously developed scAge approach that was developed for the analysis of single cell DNA methylation datasets. BayesAge utilizes maximum likelihood estimation (MLE) to infer ages, models count data using binomial distributions, and uses LOWESS smoothing to capture the non-linear dynamics between methylation and age. Our approach is designed for use with bulk bisulfite sequencing datasets. BayesAge outperforms scAge in several respects. Specifically, BayesAge’s age residuals are not age associated, thus providing a less biased representation of epigenetic age variation across populations. Moreover, BayesAge enables the estimation of error bounds on age inference and, when run on downsampled data, its coefficient of determination between predicted and actual ages surpasses both scAge and penalized regression.
2024-03-20 | GSE261769 | GEO
Project description:Standard statistical practice used for determining the relative importance of competing causes of disease typically relies on ad hoc methods, often byproducts of machine learning procedures (stepwise regression, random forest, etc.). Causal inference framework and data-adaptive methods may help to tailor parameters to match the clinical question and free one from arbitrary modeling assumptions. Our focus is on implementations of such semiparametric methods for a variable importance measure (VIM). We propose a fully automated procedure for VIM based on collaborative targeted maximum likelihood estimation (cTMLE), a method that optimizes the estimate of an association in the presence of potentially numerous competing causes. We applied the approach to data collected from traumatic brain injury patients, specifically a prospective, observational study including three US Level-1 trauma centers. The primary outcome was a disability score (Glasgow Outcome Scale - Extended (GOSE)) collected three months post-injury. We identified clinically important predictors among a set of risk factors using a variable importance analysis based on targeted maximum likelihood estimators (TMLE) and on cTMLE. Via a parametric bootstrap, we demonstrate that the latter procedure has the potential for robust automated estimation of variable importance measures based upon machine-learning algorithms. The cTMLE estimator was associated with substantially less positivity bias as compared to TMLE and larger coverage of the 95% CI. This study confirms the power of an automated cTMLE procedure that can target model selection via machine learning to estimate VIMs in complicated, high-dimensional data.
| S-EPMC5589499 | biostudies-literature
Project description:Genetic dissection of the S rat genome has provided strong evidence for the presence of two interacting blood pressure (BP) quantitative trait loci (QTLs), termed QTL1 and QTL2, on rat chromosome 5. However, the identities of the underlying interacting genetic factors remain unknown. Further experiments targeted to identify the interacting genetic factors by the substitution mapping approach alone are difficult because of the interdependency of natural recombinations to occur at the two QTLs. We hypothesized that the interacting genetic factors underlying these two QTLs may interact at the level of gene transcription and thereby represent expression QTLs (eQTLs). To detect these interacting eQTLs, a custom QTL chip containing the annotated genes within QTL1 and QTL2 was developed and used to conduct a transcriptional profiling study of S and two congenic strains that retain either one or both the QTLs. The results uncovered an interaction between two transcription factors, DMRTA2 and NFIA. Further, the âbiological signatureâ elicited by these two transcription factors was differential between the congenic strain that retained LEW alleles at both QTL1 and 2 compared to the congenic strain that retained LEW alleles at QTL1 alone. A network of transcription factors potentially affecting BP could be traced, lending support to our hypothesis. Pairs of Cy5 and Cy3 labeled targets were co-hybridized onto either a custom long oligonucleotide microarray for the interrogation of 231 genes encompassed by QTL1 and QTL2, or a TIGR rat cDNA array consisting of 26,401 probe elements representing 20,465 unique non-QTL genes. A âflip-dyeâ or âbalanced blockâ design was used as the experimental method of choice to account for potential dye-bias labeling effects. Six âbalanced blockâ normalized files are submitted for the long oligonucleotide array interrogating the hearts from S versus S.LEW(5)x6x9 animals, fourteen âflip-dyeâ hybridizations are submitted for the cDNA array interrogating the hearts from S versus S.LEW(5)x6x9 animals, and twelve âflip-dyeâ hybridizations are submitted for the cDNA array interrogating the hearts from S versus S.LEW(5)x6x11 animals. GPR and MEV files cannot be located.
2010-06-25 | E-GEOD-7628 | biostudies-arrayexpress