Genetic variation in vulnerability to the behavioral effects of neonatal hippocampal damage in rats.

We explored how two independent variables, one genetic (i.e., specific rat strains) and another environmental (i.e., a developmental excitotoxic hippocampal lesion), contribute to phenotypic variation. Sprague-Dawley (SD), Fischer 344 (F344), and Lewis rats underwent two grades of neonatal excitotoxic damage: small and large ventral hippocampal (SVH and LVH) lesions. Locomotion was tested before puberty [postnatal day 35 (P35)] and after puberty (P56) following exposure to a novel environment or administration of amphetamine. The behavioral effects were strain- and lesion-specific. As shown previously, SD rats with LVH lesions displayed enhanced spontaneous and amphetamine-induced locomotion as compared with controls at P56, but not at P35. SVH lesions in SD rats had no effect at any age. In F344 rats with LVH lesions, enhanced spontaneous and amphetamine-induced locomotion appeared early (P35) and was exaggerated at P56. SVH lesions in F344 rats resulted in a pattern of effects analogous to LVH lesions in SD rats--i.e., postpubertal onset of hyperlocomotion (P56). In Lewis rats, LVH lesions had no significant effect on novelty- or amphetamine-induced locomotion at any age. These data show that the degree of genetic predisposition and the extent of early induced hippocampal defect contribute to the particular pattern of behavioral outcome. These results may have implications for modeling interactions of genetic and environmental factors involved in schizophrenia, a disorder characterized by phenotypic heterogeneity, genetic predisposition, a developmental hippocampal abnormality, and vulnerability to environmental stress.

Coreferentiality: A New Method for the Hypothesis-Based Analysis of Phenotypes Characterized by Multivariate Data

Many multifactorial biologic effects, particularly in the context of complex human diseases, are still poorly understood. At the same time, the systematic acquisition of multivariate data has become increasingly easy. The use of such data to analyze and model complex phenotypes, however, remains a challenge. Here, a new analytic approach is described, termed coreferentiality, together with an appropriate statistical test. Coreferentiality is the indirect relation of two variables of functional interest in respect to whether they parallel each other in their respective relatedness to multivariate reference data, which can be informative for a complex effect or phenotype. It is shown that the power of coreferentiality testing is comparable to multiple regression analysis, sufficient even when reference data are informative only to a relatively small extent of 2.5%, and clearly exceeding the power of simple bivariate correlation testing. Thus, coreferentiality testing uses the increased power of multivariate analysis, however, in order to address a more straightforward interpretable bivariate relatedness. Systematic application of this approach could substantially improve the analysis and modeling of complex phenotypes, particularly in the context of human study where addressing functional hypotheses by direct experimentation is often difficult.

Statistical Hurdle Models for Single Cell Gene Expression: Differential Expression and Graphical Modeling

Thesis (Ph.D.)--University of Washington, 2016-06

Multivariate QTL linkage analysis suggests a QTL for platelet count on chromosome 19q

Platelet count is a highly heritable trait with genetic factors responsible for around 80% of the phenotypic variance. We measured platelet count longitudinally in 327 monozygotic and 418 dizygotic twin pairs at 12, 14 and 16 years of age. We also performed a genome-wide linkage scan of these twins and their families in an attempt to localize QTLs that influenced variation in platelet concentrations. Suggestive linkage was observed on chromosome 19q13.13-19q13.31 at 12 (LOD=2.12, P=0.0009), 14 (LOD=2.23, P=0.0007) and 16 (LOD=1.01, P=0.016) years of age and multivariate analysis of counts at all three ages increased the LOD to 2.59 (P=0.0003). A possible candidate in this region is the gene for glycoprotein VI, a receptor involved in platelet aggregation. Smaller linkage peaks were also seen at 2p, 5p, 5q, 10p and 15q. There was little evidence for linkage to the chromosomal regions containing the genes for thrombopoietin (3q27) and the thrombopoietin receptor (1q34), suggesting that polymorphisms in these genes do not contribute substantially to variation in platelet count between healthy individuals.

Using identity-by-de scent information in affected sib pairs to increase the efficiency of genetic association studies

T he aim of this study was to determine whether identity-by-descent (IBD) information for affected sib pairs (ASPs) can be used to select a sample of cases for a genetic case-control study which will provide more power for detecting association with loci in a known linkage region. By modeling the expected frequency of the disease allele in ASPs showing IBD sharing of 0, 1, or 2 alleles, and considering additive, recessive, and dominant disease models, we show that cases selected from IBD 2 families are best for this purpose, followed by those selected from IBD 1 families; least useful are cases selected from IBD 0 families.

Genetic and environmental sources of covariance between reading tests used in neuropsychological assessment and IQ subtests

In this study, we examined genetic and environmental influences on covariation among two reading tests used in neuropsychological assessment (Cambridge Contextual Reading Test [CCRT], [Beardsall, L., and Huppert, F. A. ( 1994). J. Clin. Exp. Neuropsychol. 16: 232 - 242], Schonell Graded Word Reading Test [SGWRT], [ Schonell, F. J., and Schonell, P. E. ( 1960). Diagnostic and attainment testing. Edinburgh: Oliver and Boyd.]) and among a selection of IQ subtests from the Multidimensional Aptitude Battery (MAB), [Jackson, D. N. (1984). Multidimensional aptitude battery, Ontario: Research Psychologists Press.] and the Wechsler Adult Intelligence Scale-Revised (WAIS-R) [Wechsler, D. (1981). Manual for the Wechsler Adult Intelligence Scale-Revised (WAIS-R). San Antonio: The Psychological Corporation]. Participants were 225 monozygotic and 275 dizygotic twin pairs aged from 15 years to 18 years ( mean, 16 years). For Verbal IQ subtests, phenotypic correlations with the reading tests ranged from 0.44 to 0.65. For Performance IQ subtests, phenotypic correlations with the reading tests ranged from 0.23 to 0.34. Results of Structural Equation Modeling (SEM) supported a model with one genetic General factor and three genetic group factors ( Verbal, Performance, Reading). Reading performance was influenced by the genetic General factor ( accounting for 13% and 20% of the variance for the CCRT and SGWRT, respectively), the genetic Verbal factor ( explaining 17% and 19% of variance for the CCRT and SGWRT), and the genetic Reading factor ( explaining 21% of the variance for both the CCRT and SGWRT). A common environment factor accounted for 25% and 14% of the CCRT and SGWRT variance, respectively. Genetic influences accounted for more than half of the phenotypic covariance between the reading tests and each of the IQ subtests. The heritabilities of the CCRT and SGWRT were 0.54 and 0.65, respectively. Observable covariance between reading assessments used by neuropsychologists to estimate IQ and IQ subtests appears to be largely due to genetic effects.

Latent class and genetic analysis does not support migraine with aura and migraine without aura as separate entities

Latent class and genetic analyses were used to identify subgroups of migraine sufferers in a community sample of 6,265 Australian twins (55% female) aged 25-36 who had completed an interview based on International Headache Society UHS) criteria. Consistent with prevalence rates from other population-based studies, 703 (20%) female and 250 (9%) male twins satisfied the IHS criteria for migraine without aura (MO), and of these, 432 (13%) female and 166 (6%) male twins satisfied the criteria for migraine with aura (MA) as indicated by visual symptoms. Latent class analysis (LCA) of IHS symptoms identified three major symptomatic classes, representing 1) a mild form of recurrent nonmigrainous headache, 2) a moderately severe form of migraine, typically without visual aura symptoms (although 40% of individuals in this class were positive for aura), and 3) a severe form of migraine typically with visual aura symptoms (although 24% of individuals were negative for aura). Using the LCA classification, many more individuals were considered affected to some degree than when using IHS criteria (35% vs. 13%). Furthermore, genetic model fitting indicated a greater genetic contribution to migraine using the LCA classification (heritability, h(2) =0.40; 95% CI, 0.29-0.46) compared with the IHS classification (h(2)=0.36; 95% CI, 0.22-0.42). Exploratory latent class modeling, fitting up to 10 classes, did not identify classes corresponding to either the IHS MO or MA classification. Our data indicate the existence of a continuum of severity, with MA more severe but not etiologically distinct from MO. In searching for predisposing genes, we should therefore expect to find some genes that may underlie all major recurrent headache subtypes, with modifying genetic or environmental factors that may lead to differential expression of the liability for migraine. (C) 2004 Wiley-Liss, Inc.

Genetic diversity and population structure of Tasmanian devils, the largest marsupial carnivore

Genetic diversity and population structure were investigated across the core range of Tasmanian devils (Sarcophilus laniarius; Dasyuridae), a wide-ranging marsupial carnivore restricted to the island of Tasmania. Heterozygosity (0.386-0.467) and allelic diversity (2.7-3.3) were low in all subpopulations and allelic size ranges were small and almost continuous, consistent with a founder effect. Island effects and repeated periods of low population density may also have contributed to the low variation. Within continuous habitat, gene flow appears extensive up to 50 km (high assignment rates to source or close neighbour populations; nonsignificant values of pairwise F-ST), in agreement with movement data. At larger scales (150-250 km), gene flow is reduced (significant pairwise F-ST) but there is no evidence for isolation by distance. The most substantial genetic structuring was observed for comparisons spanning unsuitable habitat, implying limited dispersal of devils between the well-connected, eastern populations and a smaller northwestern population. The genetic distinctiveness of the northwestern population was reflected in all analyses: unique alleles; multivariate analyses of gene frequency (multidimensional scaling, minimum spanning tree, nearest neighbour); high self-assignment (95%); two distinct populations for Tasmania were detected in isolation by distance and in Bayesian model-based clustering analyses. Marsupial carnivores appear to have stronger population subdivisions than their placental counterparts.

Experimental evidence for multivariate stabilizing sexual selection

Stabilizing selection is a fundamental concept in evolutionary biology. In the presence of a single intermediate optimum phenotype (fitness peak) on the fitness surface, stabilizing selection should cause the population to evolve toward such a peak. This prediction has seldom been tested, particularly for suites of correlated traits. The lack of tests for an evolutionary match between population means and adaptive peaks may be due, at least in part, to problems associated with empirically detecting multivariate stabilizing selection and with testing whether population means are at the peak of multivariate fitness surfaces. Here we show how canonical analysis of the fitness surface, combined with the estimation of confidence regions for stationary points on quadratic response surfaces, may be used to define multivariate stabilizing selection on a suite of traits and to establish whether natural populations reside on the multivariate peak. We manufactured artificial advertisement calls of the male cricket Teleogryllus commodus and played them back to females in laboratory phonotaxis trials to estimate the linear and nonlinear sexual selection that female phonotactic choice imposes on male call structure. Significant nonlinear selection on the major axes of the fitness surface was convex in nature and displayed an intermediate optimum, indicating multivariate stabilizing selection. The mean phenotypes of four independent samples of males, from the same population as the females used in phonotaxis trials, were within the 95% confidence region for the fitness peak. These experiments indicate that stabilizing sexual selection may play an important role in the evolution of male call properties in natural populations of T. commodus.

A reassessment of genetic limits to evolutionary change

An absence of genetic variance in traits under selection is perhaps the oldest explanation for a limit to evolutionary change, but has also been the most easily dismissed. We review a range of theoretical and empirical results covering single traits to more complex multivariate systems, and show that an absence of genetic variance may be more common than is currently appreciated. From a single-trait perspective, we highlight that it is becoming clear that some trait types do not display significant levels of genetic variation, and we raise the possibility that species with restricted ranges may differ qualitatively from more widespread species in levels of genetic variance in ecologically important traits. A common misconception in many life-history studies is that a lack of genetic variance in single traits, and genetic constraints as a consequence of bivariate genetic correlations, are different causes of selection limits. We detail how interpretations of bivariate patterns are unlikely to demonstrate genetic limits to selection in many cases. We advocate a multivariate definition of genetic constraints that emphasizes the presence (or otherwise) of genetic variance in the multivariate direction of selection. For multitrait systems, recent results using longer term studies of organisms, in which more is understood concerning what traits may be under selection, have indicated that selection may exhaust genetic variance, resulting in a limit to the selection response.

Simulation of Brugada syndrome using cellular and three-dimensional whole-heart modeling approaches

Brugada syndrome (BS) is a genetic disease identified by an abnormal electrocardiogram ( ECG) ( mainly abnormal ECGs associated with right bundle branch block and ST-elevation in right precordial leads). BS can lead to increased risk of sudden cardiac death. Experimental studies on human ventricular myocardium with BS have been limited due to difficulties in obtaining data. Thus, the use of computer simulation is an important alternative. Most previous BS simulations were based on animal heart cell models. However, due to species differences, the use of human heart cell models, especially a model with three-dimensional whole-heart anatomical structure, is needed. In this study, we developed a model of the human ventricular action potential (AP) based on refining the ten Tusscher et al (2004 Am. J. Physiol. Heart Circ. Physiol. 286 H1573 - 89) model to incorporate newly available experimental data of some major ionic currents of human ventricular myocytes. These modified channels include the L-type calcium current (ICaL), fast sodium current (I-Na), transient outward potassium current (I-to), rapidly and slowly delayed rectifier potassium currents (I-Kr and I-Ks) and inward rectifier potassium current (I-Ki). Transmural heterogeneity of APs for epicardial, endocardial and mid-myocardial (M) cells was simulated by varying the maximum conductance of IKs and Ito. The modified AP models were then used to simulate the effects of BS on cellular AP and body surface potentials using a three-dimensional dynamic heart - torso model. Our main findings are as follows. (1) BS has little effect on the AP of endocardial or mid-myocardial cells, but has a large impact on the AP of epicardial cells. (2) A likely region of BS with abnormal cell AP is near the right ventricular outflow track, and the resulting ST-segment elevation is located in the median precordium area. These simulation results are consistent with experimental findings reported in the literature. The model can reproduce a variety of electrophysiological behaviors and provides a good basis for understanding the genesis of abnormal ECG under the condition of BS disease.

Stochastic models for regulatory networks of the genetic toggle switch

Bistability arises within a wide range of biological systems from the A phage switch in bacteria to cellular signal transduction pathways in mammalian cells. Changes in regulatory mechanisms may result in genetic switching in a bistable system. Recently, more and more experimental evidence in the form of bimodal population distributions indicates that noise plays a very important role in the switching of bistable systems. Although deterministic models have been used for studying the existence of bistability properties under various system conditions, these models cannot realize cell-to-cell fluctuations in genetic switching. However, there is a lag in the development of stochastic models for studying the impact of noise in bistable systems because of the lack of detailed knowledge of biochemical reactions, kinetic rates, and molecular numbers. in this work, we develop a previously undescribed general technique for developing quantitative stochastic models for large-scale genetic regulatory networks by introducing Poisson random variables into deterministic models described by ordinary differential equations. Two stochastic models have been proposed for the genetic toggle switch interfaced with either the SOS signaling pathway or a quorum-sensing signaling pathway, and we have successfully realized experimental results showing bimodal population distributions. Because the introduced stochastic models are based on widely used ordinary differential equation models, the success of this work suggests that this approach is a very promising one for studying noise in large-scale genetic regulatory networks.

Studying phenotypic evolution using multivariate quantitative genetics

Quantitative genetics provides a powerful framework for studying phenotypic evolution and the evolution of adaptive genetic variation. Central to the approach is G, the matrix of additive genetic variances and covariances. G summarizes the genetic basis of the traits and can be used to predict the phenotypic response to multivariate selection or to drift. Recent analytical and computational advances have improved both the power and the accessibility of the necessary multivariate statistics. It is now possible to study the relationships between G and other evolutionary parameters, such as those describing the mutational input, the shape and orientation of the adaptive landscape, and the phenotypic divergence among populations. At the same time, we are moving towards a greater understanding of how the genetic variation summarized by G evolves. Computer simulations of the evolution of G, innovations in matrix comparison methods, and rapid development of powerful molecular genetic tools have all opened the way for dissecting the interaction between allelic variation and evolutionary process. Here I discuss some current uses of G, problems with the application of these approaches, and identify avenues for future research.

Genetic Influences on Exercise Participation in 37.051 Twin Pairs from Seven Countries

Background A sedentary lifestyle remains a major threat to health in contemporary societies. To get more insight in the relative contribution of genetic and environmental influences on individual differences in exercise participation, twin samples from seven countries participating in the GenomEUtwin project were used. Methodology Self-reported data on leisure time exercise behavior from Australia, Denmark, Finland, Norway, the Netherlands, Sweden and United Kingdom were used to create a comparable index of exercise participation in each country (60 minutes weekly at a minimum intensity of four metabolic equivalents). Principal Findings Modest geographical variation in exercise participation was revealed in 85,198 subjects, aged 19–40 years. Modeling of monozygotic and dizygotic twin resemblance showed that genetic effects play an important role in explaining individual differences in exercise participation in each country. Shared environmental effects played no role except for Norwegian males. Heritability of exercise participation in males and females was similar and ranged from 48% to 71% (excluding Norwegian males). Conclusions Genetic variation is important in individual exercise behavior and may involve genes influencing the acute mood effects of exercise, high exercise ability, high weight loss ability, and personality. This collaborative study suggests that attempts to find genes influencing exercise participation can pool exercise data across multiple countries and different instruments