Variance components models for gene-environment interaction in twin analysis.

Gene-environment interaction is likely to be a common and important source of variation for complex behavioral traits. Often conceptualized as the genetic control of sensitivity to the environment, it can be incorporated in variance components twin analyses by partitioning genetic effects into a mean part, which is independent of the environment, and a part that is a linear function of the environment. The model allows for one or more environmental moderator variables (that possibly interact with each other) that may i). be continuous or binary ii). differ between twins within a pair iii). interact with residual environmental as well as genetic effects iv) have nonlinear moderating properties v). show scalar (different magnitudes) or qualitative (different genes) interactions vi). be correlated with genetic effects acting upon the trait, to allow for a test of gene-environment interaction in the presence of gene-environment correlation. Aspects and applications of a class of models are explored by simulation, in the context of both individual differences twin analysis and, in a companion paper (Purcell & Sham, 2002) sibpair quantitative trait locus linkage analysis. As well as elucidating environmental pathways, consideration of gene-environment interaction in quantitative and molecular studies will potentially direct and enhance gene-mapping efforts.     

Bias, precision and heritability of self-reported and clinically measured height in Australian twins

Many studies of quantitative and disease traits in human genetics rely upon self-reported measures. Such measures are based on questionnaires or interviews and are often cheaper and more readily available than alternatives. However, the precision and potential bias cannot usually be assessed. Here we report a detailed quantitative genetic analysis of stature. We characterise the degree of measurement error by utilising a large sample of Australian twin pairs (857 MZ, 815 DZ) with both clinical and self-reported measures of height. Self-report height measurements are shown to be more variable than clinical measures. This has led to lowered estimates of heritability in many previous studies of stature. In our twin sample the heritability estimate for clinical height exceeded 90%. Repeated measures analysis shows that 2–3 times as many self-report measures are required to recover heritability estimates similar to those obtained from clinical measures. Bivariate genetic repeated measures analysis of self-report and clinical height measures showed an additive genetic correlation >0.98. We show that the accuracy of self-report height is upwardly biased in older individuals and in individuals of short stature. By comparing clinical and self-report measures we also showed that there was a genetic component to females systematically reporting their height incorrectly; this phenomenon appeared to not be present in males. The results from the measurement error analysis were subsequently used to assess the effects of error on the power to detect linkage in a genome scan. Moderate reduction in error (through the use of accurate clinical or multiple self-report measures) increased the effective sample size by 22%; elimination of measurement error led to increases in effective sample size of 41%.     

A twin study of the etiology of prolonged fatigue and immune activation.

Risk factors to prolonged fatigue syndromes (PFS) are controversial. Pre-morbid and/or current psychiatric disturbance, and/or disturbed cell-mediated immunity (CMI), have been proposed as etiologic factors. Self-report measures of fatigue and psychologic distress and three in vitro measures of CMI were collected from 124 twin pairs. Crosstwin-crosstrait correlations were estimated for the complete monozygotic (MZ; 79 pairs) and dizygotic (DZ; 45 pairs) twin groups. Multivariate genetic and environmental models were fitted to explore the patterns of covariation between etiologic factors. For fatigue, the MZ correlation was more than double the DZ correlation (0.49 versus 0.16) indicating strong genetic control of familial aggregation. By contrast, for in vitro immune activation measures MZ and DZ correlations were similar (0.49-0.69 versus 0.42-0.53) indicating the etiologic role of shared environments. As small univariate associations were noted between prolonged fatigue and the in vitro immune measures (r = -0.07 to -0.12), multivariate models were fitted. Relevant etiologic factors included: a common genetic factor accounting for 48% of the variance in fatigue which also accounted for 4%, 6% and 8% reductions in immune activation; specific genetic factors for each of the in vitro immune measures; a shared environment factor influencing the three immune activation measures; and, most interestingly, unique environmental influences which increased fatigue but also increased markers of immune activation. PFS that are associated with in vitro measures of immune activation are most likely to be the consequence of current environmental rather than genetic factors. Such environmental factors could include physical agents such as infection and/or psychologic stress.     

Genetic and environmental correlations among serum lipids and apolipoproteins in elderly twins reared together and apart.

Genetic and environmental correlations among five serum-lipid measures were examined in the Swedish Adoption/Twin Study of Aging. The sample included 302 twin pairs; 146 of these twin pairs were separated at an early age and were reared apart. The lipid measures examined include total cholesterol, HDL-cholesterol, triglycerides, and apolipoproteins A-I and B. Genetic and environmental correlations were evaluated for two different age groups, formed by dividing the sample at the median. The younger group included individuals 41.8-65.4 years of age at the midpoint of testing, although only 24 individuals were < 50 years of age. The older group included all those > 65.4 years of age, up to age 87 years of age. Substantial genetic correlations were found within each age group, although there is no evidence for a single genetic factor common to all five lipids. The comparison of twins reared together with twins reared apart allowed estimation of the effects of shared rearing environment; however, shared rearing environment only appears to be a significant mediator of the phenotypic correlation between apolipoprotein B and cholesterol in the older group. Examination of the genetic and environmental covariances suggests that the relative contributions of genetic factors are lower in the older group. Nonshared environmental factors are relatively more important mediators of phenotypic correlations among the serum lipids in individuals > 65.4 years of age than they are for the younger group. Sex differences in the mediation of these serum lipids were not as clear.     

Power of the classical twin design revisited.

Statistical power of the classical twin design was revisited. The approximate sampling variances of a least-squares estimate of the heritability in a univariate analysis and estimate of the genetic correlation coefficient in a bivariate analysis were derived analytically for the ACE model. Statistical power to detect additive genetic variation under the ACE model was derived analytically for least-squares, goodness-of-fit and maximum likelihood-based test statistics. The noncentrality parameter for the likelihood ratio test statistic is shown to be a simple function of the MZ and DZ intraclass correlation coefficients and the proportion of MZ and DZ twin pairs in the sample. All theoretical results were validated using simulation. The derived expressions can be used to calculate power of the classical twin design in a simple and rapid manner.     

Are there differences in brain morphometry between twins and unrelated singletons? A pediatric MRI study

Twins provide a unique capacity to explore relative genetic and environmental contributions to brain development, but results are applicable to non‐twin populations only to the extent that twin and singleton brains are alike. A reason to suspect differences is that as a group twins are more likely than singletons to experience adverse prenatal and perinatal events that may affect brain development. We sought to assess whether this increased risk leads to differences in child or adolescent brain anatomy in twins who do not experience behavioral or neurological sequelae during the perinatal period. Brain MRI scans of 185 healthy pediatric twins (mean age = 11.0, SD = 3.6) were compared to scans of 167 age‐ and sex‐matched unrelated singletons on brain structures measured, which included gray and white matter lobar volumes, ventricular volume, and area of the corpus callosum. There were no significant differences between groups for any structure, despite sufficient power for low type II (i.e. false negative) error. The implications of these results are twofold: (1) within this age range and for these measures, it is appropriate to include healthy twins in studies of typical brain development, and (2) findings regarding heritability of brain structures obtained from twin studies can be generalized to non‐twin populations.     

Glucocorticoid receptor binding in twin pairs is affected by shared environment but not by shared genes

We set out to determine whether glucocorticoid receptor activity is affected mainly by genetic or environmental factors. The affinity and capacity of the glucocorticoid receptor was measured using dexamethasone binding in whole leukocytes from 53 monozygotic and 48 dizygotic twin pairs. Receptor binding characteristics assayed from twin pairs on the same day were highly correlated within twin pairs irrespective of zygosity. Apparent Kd was negatively correlated with environmental temperature (R²=0.13, P<0.0001) but this did not confound the intra-pair correlation, suggesting a strong familial component independent of zygosity. Receptor binding parameters were not more closely correlated in monozygotic twins than dizygotic twin pairs indicating that there is no major genetic contribution to receptor binding and that environmental influences predominate. The close similarity in binding between twin pairs in adulthood raises the possibility that familial, non-genetic, factors such as shared early life environment may programme the glucocorticoid receptor.     

Using bivariate models to understand between- and within-cluster regression coefficients, with application to twin data

In the regression analysis of clustered data it is important to allow for the possibility of distinct between- and within-cluster exposure effects on the outcome measure, represented, respectively, by regression coefficients for the cluster mean and the deviation of the individual-level exposure value from this mean. In twin data, the within-pair regression effect represents association conditional on exposures shared within pairs, including any common genetic or environmental influences on the outcome measure. It has therefore been proposed that a comparison of the within-pair regression effects between monozygous (MZ) and dizygous (DZ) twins can be used to examine whether the association between exposure and outcome has a genetic origin. We address this issue by proposing a bivariate model for exposure and outcome measurements in twin-pair data. The between- and within-pair regression coefficients are shown to be weighted averages of ratios of the exposure and outcome variances and covariances, from which it is straightforward to determine the conditions under which the within-pair regression effect in MZ pairs will be different from that in DZ pairs. In particular, we show that a correlation structure in twin pairs for exposure and outcome that appears to be due to genetic factors will not necessarily be reflected in distinct MZ and DZ values for the within-pair regression coefficients. We illustrate these results in a study of female twin pairs from Australia and North America relating mammographic breast density to weight and body mass index.     

Parent ratings of temperament in twins: explaining the 'too low' DZ correlations.

Twin studies of child temperament using objective measures consistently suggest moderate heritability for most dimensions. However, parent rating measures produce unusual patterns of results. Intraclass correlations for identical (MZ) twins are typically high, whereas fraternal (DZ) twin intraclass correlations are much lower than would be predicted from an additive genetic model. The 'too low' DZ correlations can be explained by parent-rating biases that either exaggerate the differences between DZ twins (contrast effects) or that inflate the similarity of MZ twins (assimilation effects), or by the presence of non-additive genetic variance. To evaluate the three possible explanations, we used model-fitting procedures applied to parent-rating data averaged across 14, 20, 24, and 36 months of age in a sample of 196 twin pairs participating in the MacArthur Longitudinal Twin Study. The data were best described by a model that included contrast effects. Implications for non-twin research are discussed.     

Effects of birth weight on later intelligence

Abstract

Identical twins are used to show the long‐term effects of birth weight when genetic and many environmental factors are controlled. Female twins between the ages of 6 and 10 are evaluated for performance IQ and other measures of activity. No relationship was found between any of these measures and birth weight except for IQ. The heavier cotwin had the higher IQ score even when the birth weight of each twin was above 2,500 grams. Birth order was not associated with intelligence within the twin pairs.     

Analysis of functional abilities for elderly Danish twins using GEE models.

In this work we present a new method for genetic analysis of twin data which is based on generalized estimating equations and allows for analysis of various response types (e.g., continuous, binary, counts) combined with estimation of residual correlations. The new approach allows for control of covariates of any kind (e.g., continuous, counts) by modeling the dependence of mean and variance on background variables. The proposed method was applied to identify the covariates that have a significant influence on elderly people's functional abilities, and find the estimates for the correlation coefficients of residuals for MZ and DZ twins in a sample of 2401 Danish twin 75 years of age or older. The bootstrap method was used to obtain standard errors for correlation coefficients. It was shown, that the chosen covariates have similar effects on MZ and DZ twins, and that the residual correlation in MZ twins is significantly higher than in DZ twins, which indicates that genetic factors play an etiological role in the determination of physical status of elderly people, controlled for 10 background variables.     

Genetic and environmental influences of surfactant protein D serum levels

The collectin surfactant protein D (SP-D) is an important component of the pulmonary innate immune system, but SP-D is also present on extrapulmonary epithelial surfaces and in serum, where it has been used as a biomarker for pulmonary disease states. In this study, we investigate the mechanisms defining the constitutional serum level of SP-D and determine the magnitude of the genetic contribution to serum SP-D in the adult population. Recent studies have demonstrated that serum SP-D concentrations in children are genetically determined and that a single nucleotide polymorphism (SNP) located in the NH(2)-terminal region (Met11Thr) of the mature protein is significantly associated with the serum SP-D levels. A classic twin study was performed on a twin population including 1,476 self-reported healthy adults. The serum SP-D levels increased with male sex, age, and smoking status. The intraclass correlation was significantly higher for monozygotic (MZ) twin pairs than for dizygotic (DZ) twin pairs. Serum SP-D variance was influenced by nonshared environmental effects and additive genetic effects. Multivariate analysis of MZ and DZ covariance matrixes showed significant genetic correlation among serum SP-D and metabolic variables. The Met11Thr variant explained a significant part of the heritability indicating that serum SP-D variance could be decomposed into non-shared environmental effects (e(2) = 0.19), additive genetic effects (h(2) = 0.42), and the effect of the Met11Thr variations (q(2) = 0.39).     

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.     

Conceptual analysis of methods applied to assessment of diversity within and distance between populations with asexual or mixed mode of reproduction

Measures of diversity within populations, and distance between populations, are compared for organisms with an asexual or mixed mode of reproduction. Examples are drawn from studies of plant pathogenic fungi based on binary traits including presence/absence of DNA bands or virulence/avirulence to differential hosts. Commonly used measures of population diversity or genetic distance consider either genotype frequencies or allele frequencies. Kosman's diversity and distance measures are the most suitable for populations with an asexual or mixed mode of reproduction, because by considering genetic patterns of all individuals they take into account not just the genotype frequencies but also the genetic similarities between genotypes in the populations. The Kosman distance and diversity measures for populations can be calculated using different measures of dissimilarity between individuals (the simple mismatch, Jaccard and Dice coefficients of dissimilarity). Kosman's distances based on the simple mismatch and Jaccard dissimilarities are metrics. Comparisons of diversity indices for hypothetical examples as well as for actual data sets are presented to demonstrate that inferences from diversity analysis of populations can be driven by techniques of diversity and distance assessments and not only data driven.     

Effects of English admixture and geographic distance on anthropometric variation and genetic structure in 19th-century Ireland

The analysis of anthropometric data often allows investigation of patterns of genetic structure in historical populations. This paper focuses on interpopulational anthropometric variation in seven populations in Ireland using data collected in the 1890s. The seven populations were located within a 120-km range along the west coast of Ireland and include islands and mainland isolates. Two of the populations (the Aran Islands and Inishbofin) have a known history of English admixture in earlier centuries. Ten anthropometric measures (head length, breadth, and height; nose length and breadth; bizygomatic and bigonial breadth; stature; hand length; and forearm length) on 259 adult Irish males were analyzed following age adjustment. Discriminant and canonical variates analysis were used to determine the degree and pattern of among-group variation. Mahalanobis' distance measure, D2, was computed between each pair of populations and compared to distance measures based on geographic distance and English admixture (a binary measure indicating whether either of a pair of populations had historical indications of admixture). In addition, surname frequencies were used to construct distance measures based on random isonymy. Correlations were computed between distance measures, and their probabilities were derived using the Mantel matrix permutation method. English admixture has the greatest effect on anthropometric variation among these populations, followed by geographic distance. The correlation between anthropometric distance and geographic distance is not significant (r = -0.081, P = .590), but the correlation of admixture and anthropometric distance is significant (r = 0.829, P = .047). When the two admixed populations are removed from the analysis the correlation between geographic and anthropometric distance becomes significant (r = 0.718, P = .025). Isonymy distance shows a significant correlation with geographic distance (r = 0.425, P = .046) but not with admixture distance (r = -0.052, P = .524). The fact that anthropometrics show past patterns of gene flow and surnames do not reflects the greater impact of stochastic processes on surnames, along with the continued extinction of surnames. This study shows that 1) anthropometrics can be extremely useful in assessing population structure and history, 2) differential gene flow into populations can have a major impact on local genetic structure, and 3) microevolutionary processes can have different effects on biological characters and surnames.     

A twin study of genetic influences on epilepsy outcome.

The identification of genetic factors that confer susceptibility to the epilepsies has to date been the focus of genetic efforts in this field. Few studies have assessed the genetic contribution to disease course in epilepsy, yet an understanding of the genetic influences on epilepsy outcome is key to developing new therapeutic strategies. The aim of this study was to assess the genetic contributions to epilepsy outcome in twin pairs concordant for epilepsy. We studied 37 epilepsy concordant twin pairs (27 monozygotic, 10 dizygotic) in whom there were no recognized environmental contributions (e.g., acquired brain injury) to epilepsy, and in whom the most likely cause for epilepsy was a shared genetic susceptibility. Clinical outcome was determined using the binary measure of Seizure Status (seizure remission or recurrence) and on a six-category ordinal Outcome Scale. Epilepsy outcome was independent of age of seizure onset, age at assessment and major epilepsy syndrome diagnosis. The proportion of twin pairs concordant for Seizure Status was 0.81 (22/27) for monozygous and 1.0 (10/10) for dizygous pairs, p = 0.3. Within-pair correlation in outcome (Outcome Scale) was 0.60 (95% CI: 0.32, 0.78) in monozygous and 0.78 (0.48, 0.92) in dizygous pairs. These data provide no evidence for genetic influences on epilepsy outcome independent of those that contribute to disease susceptibility. The observed high correlations for outcome suggest that, for epilepsy, susceptibility genes also have a major influence on outcome.     

Twin removal in genetic algorithms for protein structure prediction using low-resolution model

This paper presents the impact of twins and the measures for their removal from the population of genetic algorithm (GA) when applied to effective conformational searching. It is conclusively shown that a twin removal strategy for a GA provides considerably enhanced performance when investigating solutions to complex ab initio protein structure prediction (PSP) problems in low-resolution model. Without twin removal, GA crossover and mutation operations can become ineffectual as generations lose their ability to produce significant differences, which can lead to the solution stalling. The paper relaxes the definition of chromosomal twins in the removal strategy to not only encompass identical, but also highly correlated chromosomes within the GA population, with empirical results consistently exhibiting significant improvements solving PSP problems.     

Longevity studies in GenomEUtwin.

Previous twin studies have indicated that approximately 25% of the variation in life span can be attributed to genetic factors and recent studies have also suggested a moderate clustering of extreme longevity within families. Here we discuss various definitions of extreme longevity and some analytical approaches with special attention to the challenges due to censored data. Lexis diagrams are provided for the Danish, Dutch, Finnish, Italian, Norwegian, and Swedish Twin registries hereby outlining possibilities for longevity studies within GenomEUtwin. We extend previous analyses of lifespan for the Danish 1870-1900 twin cohorts to include the new 1901-1910 cohorts, which are consistent with the previous findings. The size of the twin cohorts in GenomEUtwin and the existence of population-based, nationwide health and death registers make epidemiological studies of longevity very powerful. The combined GenomEUtwin sample will also allow detailed age-specific heritability analyses of lifespan. Finally, it will provide a resource for identifying unusual sibships (i.e., dizygotic twin pairs) where both survived to extreme ages, as a basis for discovering genetic variants of importance for extreme survival.     

Increased genetic variance of BMI with a higher prevalence of obesity

Background and objectives

There is no doubt that the dramatic worldwide increase in obesity prevalence is due to changes in environmental factors. However, twin studies suggest that genetic differences are responsible for the major part of the variation in body mass index (BMI) and other measures of body fatness within populations. Several recent studies suggest that the genetic effects on adiposity may be stronger when combined with presumed risk factors for obesity. We tested the hypothesis that a higher prevalence of obesity and overweight and a higher BMI mean is associated with a larger genetic variation in BMI.

Methods

The data consisted of self-reported height and weight from two Danish twin surveys in 1994 and 2002. A total of 15,017 monozygotic and dizygotic twin pairs were divided into subgroups by year of birth (from 1931 through 1982) and sex. The genetic and environmental variance components of BMI were calculated for each subgroup using the classical twin design. Likewise, the prevalence of obesity, prevalence of overweight and the mean of the BMI distribution was calculated for each subgroup and tested as explanatory variables in a random effects meta-regression model with the square root of the additive genetic variance (equal to the standard deviation) as the dependent variable.

Results

The size of additive genetic variation was positively and significantly associated with obesity prevalence (p = 0.001) and the mean of the BMI distribution (p = 0.015). The association with prevalence of overweight was positive but not statistically significant (p = 0.177).

Conclusion

The results suggest that the genetic variation in BMI increases as the prevalence of obesity, prevalence of overweight and the BMI mean increases. The findings suggest that the genes related to body fatness are expressed more aggressively under the influence of an obesity-promoting environment.     

Familial resemblance for anthropometric traits in dizygotic twins and siblings

Familial resemblance for fifty anthropometric traits was studied on a sample of 45 MZ, 101 DZ twin pairs and their 125 singleton siblings. Intraclass correlation coefficients were significant for all the traits. However, resemblance within DZ twin pairs was significantly greater than within sibs for 22 variables, showing that the former had a more correlated environment than the latter. The study also showed that head and facial traits were relatively more stable to the environmental factors than the body traits and hence more suitable for cross-cultural comparisons. The study listed measures of girth and skinfold, thickness as the most labile traits.