Fitting genetic models to twin data with binary and ordered categorical responses: A comparison of structural equation modelling and Bayesian hierarchical models

We compare Bayesian methodology utilizing free-ware BUGS (Bayesian Inference Using Gibbs Sampling) with the traditional structural equation modelling approach based on another free-ware package, Mx. Dichotomous and ordinal (three category) twin data were simulated according to different additive genetic and common environment models for phenotypic variation. Practical issues are discussed in using Gibbs sampling as implemented by BUGS to fit subject-specific Bayesian generalized linear models, where the components of variation may be estimated directly. The simulation study (based on 2000 twin pairs) indicated that there is a consistent advantage in using the Bayesian method to detect a correct model under certain specifications of additive genetics and common environmental effects. For binary data, both methods had difficulty in detecting the correct model when the additive genetic effect was low (between 10 and 20%) or of moderate range (between 20 and 40%). Furthermore, neither method could adequately detect a correct model that included a modest common environmental effect (20%) even when the additive genetic effect was large (50%). Power was significantly improved with ordinal data for most scenarios, except for the case of low heritability under a true ACE model. We illustrate and compare both methods using data from 1239 twin pairs over the age of 50 years, who were registered with the Australian National Health and Medical Research Council Twin Registry (ATR) and presented symptoms associated with osteoarthritis occurring in joints of the hand.

New genetic models and mesenchymal cells as tools to ameliorate anti-tumor immunity

Résumé Le but final de ce projet est d'utiliser des cellules T ou des cellules souches mésenchymateuses modifiées génétiquement afin de surexprimer localement les deux chémokines CXCL13 et CCL2 ensemble ou chacune séparément à l'intérieur d'une tumeur solide. CXCL13 est supposé induire des structures lymphoïdes ectopiques. Un niveau élevé de CCL2 est présumé initier une inflammation aiguë. La combinaison des deux effets amène à un nouveau modèle d'étude des mécanismes régulateur de la tolérance périphérique et de l'immunité tumorale. Les connaissances acquises grâce à ce modèle pourraient permettre le développement ou l'amélioration des thérapies immunes du cancer. Le but premier de ce travail a été l'établissement d'un modèle génétique de la souris permettant d'exprimer spécifiquement dans la tumeur les deux chémokines d'intérêt à des niveaux élevés. Pour accomplir cette tâche, qui est en fait une thérapie génétique de tumeurs solides, deux types de cellules porteuses potentielles ont été évaluées. Des cellules CD8+ T et des cellules mésenchymateuses de la moelle osseuse transférées dans des receveurs portant une tumeur. Si on pouvait répondre aux besoins de la thérapie génétique, indépendamment de la thérapie immune envisagée, on posséderait là un outil précieux pour bien d'autres approches thérapeutiques. Plusieurs lignées de souris transgéniques ont été générées comme source de cellules CD8+ T modifiées afin d'exprimer les chémokines d'intérêt. Dans une approche doublement transgénique les propriétés de deux promoteurs spécifiques de cellules T ont été combinées en utilisant la technologie Cre-loxP. Le promoteur de granzyme B confère une dépendance d'activation et le promoteur distal de lck assure une forte expression constitutive dès que les cellules CD8+ T ont été activées. Les transgènes construits ont montré une bonne performance in vivo et des souris qui expriment CCL2 dans des cellules CD8+ T activées ont été obtenues. Ces cellules peuvent maintenant être utilisées avec différents protocoles pour transférer des cellules T cytotoxiques (CTL) dans des receveurs porteur d'une tumeur, permettant ainsi d'évaluer leur capacité en tant que porteuse de chémokine d'infiltrer la tumeur. L'établissement de souris transgéniques, qui expriment pareillement CXCL13 est prévu dans un avenir proche. L'évaluation de cellules mésenchymateuses de la moelle osseuse a démontré que ces cellules se greffent efficacement dans le stroma tumoral suite à la co-injection avec des cellules tumorales. Cela représente un outil précieux pour la recherche, vu qu'il permet d'introduire des cellules manipulées dans un modèle tumoral. Les résultats confirment partiellement d'autres résultats rapportés dans un modèle amélioré. Cependant, l'efficacité et la spécificité suggérées de la migration systémique de cellules mésenchymateuses de la moelle osseuse dans une tumeur n'ont pas été observées dans notre modèle, ce qui indique, que ces cellules ne se prêtent pas à une utilisation thérapeutique. Un autre résultat majeur de ce travail est l'établissement de cultures de cellules mésenchymateuses de la moelle osseuse in vitro conditionnées par des tumeurs, ce qui a permis à ces cellules de s'étendre plus rapidement en gardant leur capacité de migration et de greffe. Cela offre un autre outil précieux, vu que la culture in vitro est un pas nécessaire pour une manipulation thérapeutique. Abstract The ultimate aim of the presented project is to use genetically modified T cells or mesenchymal stem cells to locally overexpress the two chemokines CXCL13 and CCL2 together or each one alone inside a solid tumor. CXCL13 is supposed to induce ectopic lymphoid structures and a high level of CCL2 is intended to trigger acute inflamation. The combination of these two effects represents a new model for studying mechanisms that regulate peripheral tolerance and tumor immunity. Gained insights may help developing or improving immunotherapy of cancer. The primary goal of the executed work was the establishment of a genetic mouse model that allows tumor-specific expression of high levels of the two chemokines of interest. For accomplishing this task, which represents gene therapy of solid tumors, two types of potentially useful carrier cells were evaluated. CD8+ T cells and mesenchymal bone marrow cells to be used in adoptive cell transfers into tumor-bearing mice. Irrespectively of the envisaged immunotherapy, satisfaction of so far unmet needs of gene therapy would be a highly valuable tool that may be employed by many other therapeutic approaches, too. Several transgenic mouse lines were generated as a source of CD8+ T cells modified to express the chemokines of interest. In a double transgenic approach the properties of two T cell-specific promoters were combined using Cre-loxP technology. The granzyme B promoter confers activation-dependency and the lck distal promoter assures strong constitutive expression once the CD8+ T cell has been activated. The constructed transgenes showed a good performance in vivo and mice expressing CCL2 in activated CD8+ T cells were obtained. These cells can now be used with different protocols for adoptively transferring cytotoxic T cells (CTL) into tumor-bearing recipients, thus allowing to study their capacity as tumor-infiltrating chemokine carrier. The establishment of transgenic mice likewisely expressing CXCL13 is expected in the near future. In addition, T cells from generated single transgenic mice that have high expression of an EGFP reporter in both CD4+ and CD8+ cells can be easily traced in vivo when setting up adoptive transfer conditions. The evaluation of mesenchymal bone marrow cells demonstrated that these cells can efficiently engraft into tumor stroma upon local coinjection with tumor cells. This represents a valuable tool for research purposes as it allows to introduce manipulated stromal cells into a tumor model. Therefore, the established engraftment model is suited for studying the envisaged immunotherapy. These results confirm to some extend previously reported results in an improved model, however, the suggested systemic tumor homing efficiency and specificity of mesenchymal bone marrow cells was not observed in our model indicating that these cells may not be suited for therapeutic use. Another major result of the presented work is the establishment oftumor-conditioned in vitro culture of mesenchymal bone marrow cells, which allowed to more rapidly expand these cells while maintaining their tumor homing and engrafting capacities. This offers another valuable tool as in vitro culture is a necessary step for therapeutic manipulations.

Effect of race, genetic population structure, and genetic models in two-locus association studies: clustering of functional renin-angiotensin system gene variants in hypertension association studies

Previous genetic association studies have overlooked the potential for biased results when analyzing different population structures in ethnically diverse populations. The purpose of the present study was to quantify this bias in two-locus association studies conducted on an admixtured urban population. We studied the genetic structure distribution of angiotensin-converting enzyme insertion/deletion (ACE I/D) and angiotensinogen methionine/threonine (M/T) polymorphisms in 382 subjects from three subgroups in a highly admixtured urban population. Group I included 150 white subjects; group II, 142 mulatto subjects, and group III, 90 black subjects. We conducted sample size simulation studies using these data in different genetic models of gene action and interaction and used genetic distance calculation algorithms to help determine the population structure for the studied loci. Our results showed a statistically different population structure distribution of both ACE I/D (P = 0.02, OR = 1.56, 95% CI = 1.05-2.33 for the D allele, white versus black subgroup) and angiotensinogen M/T polymorphism (P = 0.007, OR = 1.71, 95% CI = 1.14-2.58 for the T allele, white versus black subgroup). Different sample sizes are predicted to be determinant of the power to detect a given genotypic association with a particular phenotype when conducting two-locus association studies in admixtured populations. In addition, the postulated genetic model is also a major determinant of the power to detect any association in a given sample size. The present simulation study helped to demonstrate the complex interrelation among ethnicity, power of the association, and the postulated genetic model of action of a particular allele in the context of clustering studies. This information is essential for the correct planning and interpretation of future association studies conducted on this population.

Comparison of voluntary ethanol intake by two pairs of rat lines used as genetic models of anxiety

The relationship between anxiety-related behaviors and voluntary ethanol intake was examined in two pairs of rat lines by the oral ethanol self-administration procedure. Floripa high (H) and low (L) rats selectively bred for contrasting anxiety responses in the open-field test, and two inbred strains, spontaneously hypertensive rats (SHR) and Lewis rats which are known to differ significantly when submitted to several behavioral tests of anxiety/emotionality, were used (9-10 animals/line/sex). No differences in the choice of ethanol solutions (2%, days 1-4, and 4%, days 5-8, respectively) in a 2-bottle paradigm were detected between Floripa H and L rats (1.94 ± 0.37 vs 1.61 ± 0.37 g/kg for ethanol intake on day 8 by the Floripa H and L rat lines, respectively). Contrary to expectations, the less anxious SHR rats consumed significantly more ethanol than Lewis rats (respective intake of 2.30 ± 0.45 and 0.72 ± 0.33 g/kg on day 8) which are known to be both addiction-prone and highly anxious. Regardless of strain, female rats consumed more ethanol than males (approximately 46%). The results showed no relationship between high anxiety and voluntary intake of ethanol for Floripa H and L rats. A negative association between these two variables, however, was found for SHR and Lewis rat strains. Data from the literature regarding the association between anxiety and alcohol intake in animal models are not conclusive, but the present results indicate that factors other than increased inborn anxiety probably lead to the individual differences in ethanol drinking behavior.

Population genetic models can be used to study the evolution of the interacting behaviours of parents and their progeny

An example of the evolution of the interacting behaviours of parents and progeny is studied using iterative equations linking the frequencies of the gametes produced by the progeny to the frequencies of the gametes in the parental generation. This population genetics approach shows that a model in which both behaviours are determined by a single locus can lead to a stable equilibrium in which the two behaviours continue to segregate. A model in which the behaviours are determined by genes at two separate loci leads eventually to fixation of the alleles at both loci but this can take many generations of selection. Models of the type described in this paper will be needed to understand the evolution of complex behaviour when genomic or experimental information is available about the genetic determinants of behaviour and the selective values of different genomes. (c) 2007 Elsevier Inc. All rights reserved.

Exact moment calculations for genetic models with migration, mutation, and drift

Using properties of moment stationarity we develop exact expressions for the mean and covariance of allele frequencies at a single locus for a set of populations subject to drift, mutation, and migration. Some general results can be obtained even for arbitrary mutation and migration matrices, for example: (1) Under quite general conditions, the mean vector depends only on mutation rates, not on migration rates or the number of populations. (2) Allele frequencies covary among all pairs of populations connected by migration. As a result, the drift, mutation, migration process is not ergodic when any finite number of populations is exchanging genes. in addition, we provide closed form expressions for the mean and covariance of allele frequencies in Wright's finite-island model of migration under several simple models of mutation, and we show that the correlation in allele frequencies among populations can be very large for realistic rates of mutation unless an enormous number of populations are exchanging genes. As a result, the traditional diffusion approximation provides a poor approximation of the stationary distribution of allele frequencies among populations. Finally, we discuss some implications of our results for measures of population structure based on Wright's F-statistics.

Testing for biases in selection on avian reproductive traits and partitioning direct and indirect selection using quantitative genetic models

Key life history traits such as breeding time and clutch size are frequently both heritable and under directional selection, yet many studies fail to document micro-evolutionary responses. One general explanation is that selection estimates are biased by the omission of correlated traits that have causal effects on fitness, but few valid tests of this exist. Here we show, using a quantitative genetic framework and six decades of life-history data on two free-living populations of great tits Parus major, that selection estimates for egg-laying date and clutch size are relatively unbiased. Predicted responses to selection based on the Robertson-Price Identity were similar to those based on the multivariate breeder’s equation, indicating that unmeasured covarying traits were not missing from the analysis. Changing patterns of phenotypic selection on these traits (for laying date, linked to climate change) therefore reflect changing selection on breeding values, and genetic constraints appear not to limit their independent evolution. Quantitative genetic analysis of correlational data from pedigreed populations can be a valuable complement to experimental approaches to help identify whether apparent associations between traits and fitness are biased by missing traits, and to parse the roles of direct versus indirect selection across a range of environments.

Estimating Two-Stage Models for Genetic Influences on Alcohol, Tobacco or Drug Use Initiation and Dependence Vulnerability in Twin and Family Data

Genetic research on risk of alcohol, tobacco or drug dependence must make allowance for the partial overlap of risk-factors for initiation of use, and risk-factors for dependence or other outcomes in users. Except in the extreme cases where genetic and environmental risk-factors for initiation and dependence overlap completely or are uncorrelated, there is no consensus about how best to estimate the magnitude of genetic or environmental correlations between Initiation and Dependence in twin and family data. We explore by computer simulation the biases to estimates of genetic and environmental parameters caused by model misspecification when Initiation can only be defined as a binary variable. For plausible simulated parameter values, the two-stage genetic models that we consider yield estimates of genetic and environmental variances for Dependence that, although biased, are not very discrepant from the true values. However, estimates of genetic (or environmental) correlations between Initiation and Dependence may be seriously biased, and may differ markedly under different two-stage models. Such estimates may have little credibility unless external data favor selection of one particular model. These problems can be avoided if Initiation can be assessed as a multiple-category variable (e.g. never versus early-onset versus later onset user), with at least two categories measurable in users at risk for dependence. Under these conditions, under certain distributional assumptions., recovery of simulated genetic and environmental correlations becomes possible, Illustrative application of the model to Australian twin data on smoking confirmed substantial heritability of smoking persistence (42%) with minimal overlap with genetic influences on initiation.

Genetic effects on preweaning weight gain of Nelore-Hereford calves according to different models and estimation methods

Additive and nonadditive genetic effects on preweaning weight gain (PWG) of a commercial crossbred population were estimated using different genetic models and estimation methods. The data set consisted of 103,445 records on purebred and crossbred Nelore-Hereford calves raised under pasture conditions on farms located in south, southeast, and middle west Brazilian regions. In addition to breed additive and dominance effects, the models including different epistasis covariables were tested. Models considering joint additive and environment (latitude) by genetic effects interactions were also applied. In a first step, analyses were carried out under animal models. In a second step, preadjusted records were analyzed using ordinary least squares (OLS) and ridge regression (RR). The results reinforced evidence that breed additive and dominance effects are not sufficient to explain the observed variability in preweaning traits of Bos taurus x Bos indicus calves, and that genotype x environment interaction plays an important role in the evaluation of crossbred calves. Data were ill-conditioned to estimate the effects of genotype x environment interactions. Models including these effects presented multicolinearity problems. In this case, RR seemed to be a powerful tool for obtaining more plausible and stable estimates. Estimated prediction error variances and variance inflation factors were drastically reduced, and many effects that were not significant under ordinary least squares became significant under RR. Predictions of PWG based on RR estimates were more acceptable from a biological perspective. In temperate and subtropical regions, calves with intermediate genetic compositions (close to 1/2 Nelore) exhibited greater predicted PWG. In the tropics, predicted PWG increased linearly as genotype got closer to Nelore. ©2006 American Society of Animal Science. All rights reserved.

Gene × environment interactions in schizophrenia: evidence from genetic mouse models

The study of gene × environment, as well as epistatic interactions in schizophrenia, has provided important insight into the complex etiopathologic basis of schizophrenia. It has also increased our understanding of the role of susceptibility genes in the disorder and is an important consideration as we seek to translate genetic advances into novel antipsychotic treatment targets. This review summarises data arising from research involving the modelling of gene × environment interactions in schizophrenia using preclinical genetic models. Evidence for synergistic effects on the expression of schizophrenia-relevant endophenotypes will be discussed. It is proposed that valid and multifactorial preclinical models are important tools for identifying critical areas, as well as underlying mechanisms, of convergence of genetic and environmental risk factors, and their interaction in schizophrenia.

Genetic Variants of Diabetes Risk and Incident Cardiovascular Events in Chronic Coronary Artery Disease

Objective: To determine whether information from genetic risk variants for diabetes is associated with cardiovascular events incidence. Methods: From the about 30 known genes associated with diabetes, we genotyped single-nucleotide polymorphisms at the 10 loci most associated with type-2 diabetes in 425 subjects from the MASS-II Study, a randomized study in patients with multi-vessel coronary artery disease. The combined genetic information was evaluated by number of risk alleles for diabetes. Performance of genetic models relative to major cardiovascular events incidence was analyzed through Kaplan-Meier curve comparison and Cox Hazard Models and the discriminatory ability of models was assessed for cardiovascular events by calculating the area under the ROC curve. Results: Genetic information was able to predict 5-year incidence of major cardiovascular events and overall-mortality in non-diabetic individuals, even after adjustment for potential confounders including fasting glycemia. Non-diabetic individuals with high genetic risk had a similar incidence of events then diabetic individuals (cumulative hazard of 33.0 versus 35.1% of diabetic subjects). The addition of combined genetic information to clinical predictors significantly improved the AUC for cardiovascular events incidence (AUC = 0.641 versus 0.610). Conclusions: Combined information of genetic variants for diabetes risk is associated to major cardiovascular events incidence, including overall mortality, in non-diabetic individuals with coronary artery disease.

Strategies for genetic model specification in the screening of genome-wide meta-analysis signals for further replication

Background Meta-analysis is increasingly being employed as a screening procedure in large-scale association studies to select promising variants for follow-up studies. However, standard methods for meta-analysis require the assumption of an underlying genetic model, which is typically unknown a priori. This drawback can introduce model misspecifications, causing power to be suboptimal, or the evaluation of multiple genetic models, which augments the number of false-positive associations, ultimately leading to waste of resources with fruitless replication studies. We used simulated meta-analyses of large genetic association studies to investigate naive strategies of genetic model specification to optimize screenings of genome-wide meta-analysis signals for further replication. Methods Different methods, meta-analytical models and strategies were compared in terms of power and type-I error. Simulations were carried out for a binary trait in a wide range of true genetic models, genome-wide thresholds, minor allele frequencies (MAFs), odds ratios and between-study heterogeneity (tau(2)). Results Among the investigated strategies, a simple Bonferroni-corrected approach that fits both multiplicative and recessive models was found to be optimal in most examined scenarios, reducing the likelihood of false discoveries and enhancing power in scenarios with small MAFs either in the presence or in absence of heterogeneity. Nonetheless, this strategy is sensitive to tau(2) whenever the susceptibility allele is common (MAF epsilon 30%), resulting in an increased number of false-positive associations compared with an analysis that considers only the multiplicative model. Conclusion Invoking a simple Bonferroni adjustment and testing for both multiplicative and recessive models is fast and an optimal strategy in large meta-analysis-based screenings. However, care must be taken when examined variants are common, where specification of a multiplicative model alone may be preferable.