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.     

Genetic and environmental determinants of dental occlusal variation in twins of different nationalities

We have compared 10 occlusal traits in 358 monozygous and dizygous twin pairs in 4 different samples and estimated genetic variances for these features. Variable and frequently nonsignificant genetic variance was noted across samples for incisal overbite and overjet, sagittal molar relationship, posterior crossbite, and rotations and displacements of anterior teeth. Heritability estimates (when appropriately calculated) were low in magnitude (0-40%) and erratic, emphasizing the importance of environmental influences on occlusal variation and the variability of apparent genetic determinants with respect to the environment or population in which they are measured.     

Chorion type as a possible influence on the results and interpretation of twin study data.

The estimation of genetic effects from twin studies usually relies upon the equal environment assumption--that monozygous (MZ) and dizygous (DZ) twin pairs experience equal similarity of their environments from prenatal experiences through adulthood. However, the sharing of a chorion may make a subset of identical twins more similar, or in some cases, more different, than twins that do not share a chorion. Recent studies suggest monochorionic MZ twins resemble one another more than dichorionic MZ twins in cognitive abilities, personality, and risk for psychiatric disorder. To the extent that prenatal environment affects these characteristics, the traditional twin method will yield biased estimates of genetic and environmental influences. We develop models for quantifying this bias and estimating the influence of chorion type on estimates of heritability.     

Bias of heritability estimates from twin data in presence of errors of zygosity determination

Simple heritability estimators of continuous as well as discrete traits from twin data are known to overestimate the degree of genetic determination of the measured traits for several reasons. Errors of zygosity determination will, however, underestimate the true heritability. The bias due to wrong assignment of dizygous twin pairs into monozygous type is evaluated here, and the results indicate that this negative bias has a compensatory effect on the estimate of the degree of genetic determination when other factors of similarity between twin pairs are taken into account. It is shown that when an estimate of zygosity error is available, the bias due to this factor can be evaluated quantitatively, and hence the adjustment for zygosity error can be incorporated in the estimation of the degree of genetic determination of a trait. Although this theory is explicitly developed here for twin studies, the general principle also applies for other types of errors of determining the degree of biological relationships for estimation of heritability, in which case this type of error may be more important than the simple zygosity error.     

Categories of deltaF508 homozygous cystic fibrosis twin and sibling pairs with distinct phenotypic characteristics.

Cystic fibrosis (CF), the most common severe autosomal recessive trait among Caucasians, is caused by molecular lesions in the cystic fibrosis transmembrane conductance regulator gene (CFTR). The course of the multi-organ disease CF is highly variable, suggesting the influence of environmental factors and/or modulating genes other than CFTR on the disease phenotype. To evaluate the cause of CF disease variability, the European CF Twin and Sibling Study collected data on two clinical parameters most sensitive for the course and prognosis of CF, ie weight predicted for height (wfh)% (representative for the nutritional status) and FEVPerc (representative for the pulmonary status) for a cohort of 277 sibling pairs, 12 pairs of dizygous twins and 29 pairs of monozygous twins. Of these 318 CF twin and sib pairs, 114 were reported to be homozygous for the most frequent CF disease-causing lesion, deltaF508. Intra-pair discordance was assessed by the intra-pair differences with wfh% and FEVPerc and by DELTA, a composite parameter defined by linear combination of wfh% and FEVPerc in order to describe discordance with respect to the overall disease severity. Monozygous twins had a significantly lower DELTA than dizygous twins (P = 0.05) indicating that CF disease severity is modulated by an inherited component in addition to the CFTR gene itself. Extreme phenotypes are considered to be more informative for the analysis of any quantitative trait. Thus, we aimed to quantify disease severity and intra-pair discordance in order to select pairs with the extreme phenotypes DIS (discordant patient pairs), CON+ (concordant and mildy affected patient pairs) and CON- (concordant and severely affected patient pairs). The algorithm reliably discriminated between pairs DIS, CON+ and CON- among the cohort of deltaF508 homozygotes. The selected pairs from these categories demonstrated non-overlapping properties for wfh%, FEVPerc and the intra-pair difference of both parameters.     

Happy families: a twin study of humour.

The objective of this study was to estimate how much of an individual's appreciation of humour is influenced by genetic factors, the shared environment or the individual's unique environment. A population-based classical twin study of 127 pairs of female twins (71 monozygous (MZ) and 56 dizygous (DZ) pairs) aged 20-75 from the St Thomas' UK Adult Twin Registry elicited responses to five 'Far Side' Larson cartoons on a scale of 0-10. Within both MZ and DZ twin pairs, the tetrachoric correlations of responses to all five cartoons were significantly greater than zero. Furthermore, the correlations for MZ and DZ twins were of similar magnitude and in some cases the DZ correlation was greater than that of the MZ twins. This pattern of correlations suggests that shared environment rather then genetic effects contributes to cartoon appreciation. Multivariate model-fitting confirmed that these data were best explained by a model that allowed for the contribution of the shared environment and random environmental factors, but not genetic effects. However, there did not appear to be a general humour factor underlying responses to all five cartoons and no effect of age was seen. The shared environment, rather than genetic factors, explains the familial aggregation of humour appreciation as assessed by the specific 'off the wall' cognitive type of cartoons used in this study.     

Inheritance of quantitative expression of erythrocyte glucose-6-phosphate dehydrogenase activity in the Negro—a twin study

Studies have been conducted on eight sets of monozygous and nine sets of dizygous female Negro twins, both members of whom were heterozygous for G-6-PD deficiency. Twins were studied both by assay of erythrocytic G-6-PD activity and by the methemoglobin elution test (MET). The MET is a procedure which identifies histochemically cells with appreciable G-6-PD activity and permits accurate determination of the percentage of such cells in heterozygotes. Monozygous twins showed significantly less within-pair variation than dizygous twins with both the MET and G-6-PD assay.Concerning the significantly greater agreement in MET results in monozygous twins than dizygous twins, our present working hypothesis is that X-chromosomal inactivation in the Negro female is genetically controlled, rather than random. However, certain alternate hypotheses allowing for random X-inactivation have not been excluded; these include somatic cell selection after random X-inactivation, and cell exchange between identical twins in utero/it. Studies in nontwin related heterozygotes now underway should help differentiate among these various possibilities.In addition to the studies on 17 pairs of female twins heterozygous for G-6-PD deficiency, 26 pairs of nondeficient female Negro twins have been studied by G-6-PD assay. Within-pair variation in monozygous twins was significantly less than within-pair variation in dizygous twins in all cases. The genetic influences detected with the G-6-PD assay in the female twins could theoretically be due to nonrandom X-inactivation, to genetically determined quantitative differences in enzyme activity (e.g., isoalleles), or to both. By appropriate calculations, based on the MET results, we have factored out the effects of X-inactivation on overall enzyme activity in the heterozygous deficient twins. After removal of the effect of X-inactivation, monozygous twins heterozygous for enzyme deficiency continue to show significantly less within-pair variation than dizygous twins. This finding indicates significant genetic influences on quantitative G-6-PD activity other than X-inactivation and other than the deficiency allele. This conclusion has been strengthened by studies on male twins where X-inactivation is not present.Supported by USPHS research grants AM-09381, HE-17544, AM-09919, and HE-03341, by USPHS Career Development Award 1-K3-AM-7959 (Dr. Brewer) and by U.S.A.E.C. Contract (11-1)-1552.     

A rare case of heteropaternal twin calves after natural mating in Brazil

Twin birth is a complex condition observed in most livestock animals, when the female gives birth to two or more offspring, generally out of the same mating. In cattle, it is a rare condition (3 to 5%) and depends on the genetic background and environmental factors. Twin birth is a result of multiple ovulations, being more common in dairy rather than in beef cattle. Calves could be monozygous or dizygous, with the same or of different sexes. When twins are born with different sexes, a sexual condition called Freemartinism occurs in between 90 to 97% of pregnancies, causing infertility in the female calf. Knowing that the twin rate is rare in commercial beef cattle, here we present an even rarer case of twin birth from two different sires after natural mating, also called heteropaternal superfecundation.     

A study of lactose absorption capacity in twins

Summary Lactose absorption capacity was determined by lactose tolerance tests with breath hydrogen determination in 102 healthy, adult, Hungarian pairs of twins in order to test monogenic Mendelian inheritance of the absorptive lactase phenotypes, lactose absorber and lactose malabsorber. Of the total, 52 pairs were monozygous (MZ) and 50 dizygous (DZ) twins of indentical sex. All MZ twins were concordant with respect to lactase phenotype. Among DZ twins, the distribution of lactase phenotypes was in agreement with Hardy-Weinberg expectations derived from the frequencies of the hypolactasia gene in DZ and MZ twins, and in the general Budapest population.In the second part of the study, three commonly used methods of lactose tolerance testing, the blood glucose, the blood galactose, and the breath hydrogen tests, were compared in 49 pairs of twins concordant for lactase phenotype. Blood galactose concentration showed the greatest and only significant difference between the intrapair correlation coefficients of MZ and DZ, and no overlap between lactose absorbers and lactose malabsorbers. The intrapair correlation coefficients of peak breath hydrogen concentration in MZ and DZ twins did not significantly differ from zero, but the resolution of lactase phenotypes was satisfactory. Differences in glucose absorption and concentration in lactose absorbers and malobsorbers overlapped considerably, and among lactose absorbers correlation coefficients in DZ were higher than in MZ twins. In MZ and DZ twins, the difference in concordance and constancy of lactose intolerance symptoms was not significant.     

Genetic and environmental interactions determine seizure susceptibility in epileptic EL mice

Gene identification has progressed rapidly for monogenic epilepsies, but complex gene-environmental interactions have hindered progress in gene identification for multifactorial epilepsies. We analyzed the role of environmental risk factors in the inheritance of multifactorial idiopathic generalized epilepsy in the EL mouse. Seizure susceptibility was evaluated in the EL (E) and seizure-resistant ABP/LeJ (A) parental mouse strains and in their AEF1 and AEF2 hybrid offspring using a handling-induced seizure test. The seizure test was administered in three environments (environments I, II and III) that differed with respect to the number of seizure tests administered (one test or four tests) and the age of the mice when tested (young or old). The inheritance of seizure susceptibility appeared dominant after repetitive seizure testing in young or old mice, but recessive after a single test in old mice. Heritability was high (0.67-0.77) in each environment. Significant quantitative trait loci (QTL) that were associated with environments I and III (repetitive testing) were found on chromosomes 2 and 9 and colocalized with previously mapped El2 and El4, respectively. The El2 QTL found in environment I associated only with female susceptibility. A novel QTL, El-N, for age-dependent predisposition to seizures was found on proximal chromosome 9 only in environment II. The findings indicate that environmental risk factors determine the genetic architecture of seizure susceptibility in EL mice and suggest that QTL for complex epilepsies should be defined in terms of the environment in which they are expressed.     

Contribution of genetic and environmental factors to variation in body compartments--a twin study in adults

This study aimed at analyzing the contribution of genetic and environmental factors on phenotypic variation of various traits of body composition. Subjects were 30 same-sexed pairs of twins including 20 monozygous (MZ) and 10 dizygous (DZ) pairs, aged 19-62 years. Zygosity was determined by DNA typing and morphological diagnosis. Body composition parameters (fat mass FM, lean body mass LBM, body cell mass BCM, extracellular mass ECM, total body water TBW, extracellular water ECW, and intracellular water ICW) were estimated by tetrapolar bioelectrical impedance analysis. Potential environmental factors influencing body composition (number of children, sporting activity and smoking behaviour) were determined by questionnaires. Heritabilities for traits of body composition were calculated by use of the twin method. Intraclass correlation is > 0.80 for the variation of LBM, BCM, ECM, TBW, ECW, and ICW in both MZ and DZ twins. Estimated heritability (h2) for FM, LBM, BCM, ECW, TBW, ECW, and ICW is 65%, 77%, 79%, 83%, 76%, 68%, and 82%, respectively. The h2 values for FM and LBM are consistent with those reported in other twin studies. For BCM, ECM, ECW and ICW, no comparative h2 estimates exist. Within-pair differences in body compartments do not change with increasing age in MZ and DZ twin pairs (p > 0.05). Stepwise multiple regression analyses indicate that zygosity, age, sex, number of children, sporting level and smoking behaviour do not significantly predict within-pair differences for weight, BMI, FM, LBM, TBW, ECW and ICW (each, p > 0.05). In contrast, sex and the number of children explain together 27% of observed within-pair differences for BCM. Zygosity is the only significant predictor of within-pair differences for ECM and height, explaining 20% (p = 0.008) and 36% of variance, respectively (p < 0.0001). Results indicate that genetic factors exert stronger influences on body composition than the considered environmental traits.     

Weight Change in Monozygotic Twins Treated with Valproate

Objective: To evaluate whether genetic factors contribute to weight gain associated with valproate (VPA) therapy. Research Methods and Procedures: We retrospectively and prospectively evaluated five pairs of monozygotic twins concordant for epilepsy and treated with VPA regarding weight course. Results: In all twin pairs, both twins showed similar weight courses under therapy with VPA. Discussion: These results suggest that genetic factors may have an influence on the weight change induced by VPA. Further studies are necessary to obtain heritability estimates regarding this effect of VPA therapy and to identify the relevant genes.     

Use of chromosome substitution strains to identify seizure susceptibility loci in mice

Seizure susceptibility varies among inbred mouse strains. Chromosome substitution strains (CSS), in which a single chromosome from one inbred strain (donor) has been transferred onto a second strain (host) by repeated backcrossing, may be used to identify quantitative trait loci (QTLs) that contribute to seizure susceptibility. QTLs for susceptibility to pilocarpine-induced seizures, a model of temporal lobe epilepsy, have not been reported, and CSS have not previously been used to localize seizure susceptibility genes. We report QTLs identified using a B6 (host) × A/J (donor) CSS panel to localize genes involved in susceptibility to pilocarpine-induced seizures. Three hundred fifty-five adult male CSS mice, 58 B6, and 39 A/J were tested for susceptibility to pilocarpine-induced seizures. Highest stage reached and latency to each stage were recorded for all mice. B6 mice were resistant to seizures and slower to reach stages compared to A/J mice. The CSS for Chromosomes 10 and 18 progressed to the most severe stages, diverging dramatically from the B6 phenotype. Latencies to stages were also significantly shorter for CSS10 and CSS18 mice. CSS mapping suggests seizure susceptibility loci on mouse Chromosomes 10 and 18. This approach provides a framework for identifying potentially novel homologous candidate genes for human temporal lobe epilepsy.     

Genetic and environmental influence on the asymmetry of dermatoglyphic traits

Fluctuating asymmetry (FA) is defined as random deviations from bilateral symmetry of the body. Thus, its magnitude is often used to evaluate developmental homeostasis. In this study we evaluate the following hypotheses: 1) FA of dermatoglyphic traits has a significant genetic component; 2) prenatal maternal environment (PME) has a significant effect on the FA of dermatoglyphic traits in developmentally healthy individuals; and 3) genetic or environmental factors affect FA on organismal or systemic levels. Therefore, their effect is better seen in composite scores of FA rather than in FA indices for single traits. We analyzed 15 dermatoglyphic traits from 140 pairs of monozygous twins, 120 pairs of dizygous twins, and 106 pairs of mothers and daughters. All individuals were developmentally healthy. The influence of genetic and environmental factors on FA was evaluated by analysis of variance and regression analysis. For a majority of the traits in our study, FA showed significant but weak heritabilities, with values falling within the 0.20-0.35 range. None of the traits taken separately demonstrated the effect of PME on FA to be significantly greater than zero. The composite score of FA tended to have greater heritability values than individual traits. One of them, obtained in principal components analysis, showed a significant PME effect, supporting the hypothesis that FA is a systemic property.     

HLA-DR typing in identical twins with insulin-dependent diabetes: difference between concordant and discordant pairs.

A total of 106 pairs of identical twins, of whom 56 were concordant and 50 discordant for insulin-dependent diabetes, were typed for HLA-DR. In both the concordant and discordant groups there was a high prevalence of the antigens DR3 and DR4, a low prevalence of DR5 and DR7, and a virtual absence of DR2. The heterozygous phenotype DR3,DR4 was more prevalent in concordant than discordant pairs. This was therefore the first demonstration of a genetic difference between concordant and discordant identical twin pairs. These findings suggest that possession of both DR3 and DR4 antigens confers a greater genetic predisposition to insulin-dependent diabetes than does the possession of either antigen alone.     

Fourier analysis of facial profiles of young twins

Twins studies provide a powerful approach to determining the relative contribution of genetics and environment to observed variation. Such studies assume trait differences in monozygous (MZ) twins are due to environmental factors and those in dizygous (DZ) twins are due to both genetic and environmental factors. This study quantitated facial profiles of twins using Fourier equations, determining their value in profile analysis and the assessment of the genetic contribution to facial shape. Standardized profile slide photographs of 79 pairs of 4-6 year-old twins (37 MZ pairs, 42 DZ pairs) were scanned and x and y coordinates were extracted from each profile using sellion and Camper's plane as references. The coordinates were subjected to Fourier analysis and the normalised vertex projection coefficients were studied. The means of the differences between coefficients for MZ co-twins did not differ significantly from that of DZ co-twins, although the DZ group showed higher mean differences in the higher harmonics. Subjective examination of superimposed reconstructions showed wider variation between DZ co-twins than MZ co-twins. Correct classification of twins by discriminant function analysis using Fourier coefficients was similar for both groups (MZ: 70.3%; DZ: 73.8%). Fourier analysis could quantitate facial profiles of young children and differentiate some details, but was unable to discriminate between genetic and environmental influences, and any possible interactions between these influences, on their overall facial profiles at this developmental stage.     

Newly established low seizure susceptible and seizure-prone inbred strains of Mongolian gerbil.

Two inbred strains of the Mongolian gerbil with different phenotypes in seizure behaviour and coat color were newly established. LSAG/Nu has low seizure susceptibility and albino phenotypes, whereas SPBG/Nu has seizure-prone and black coat color phenotypes. LSAG was compared with SPBG as to seizure incidence and grade. Mean ages at seizure onset of LSAG and SPBG were 6 and 3 months, respectively. Seizure incidences in over 9 months old LSAG and SPBG gerbils were 37.3% (66/177) and 95.2% (118/124), respectively. LSAG has a significantly lower incidence (p < 0.001) and grade (p < 0.001) of seizures than SPBG. Only a few seizing LSAG gerbils exhibited myoclonus to tonic-clonic seizure progression. These results suggest that LSAG has some mechanisms which delay the onset of seizures and prevent them from becoming serious. Both strains of gerbils can be expected to be useful animal models for the study of human idiopathic generalized epilepsy.     

Dynamic Imaging of Coherent Sources Reveals Different Network Connectivity Underlying the Generation and Perpetuation of Epileptic Seizures

The concept of focal epilepsies includes a seizure origin in brain regions with hyper synchronous activity (epileptogenic zone and seizure onset zone) and a complex epileptic network of different brain areas involved in the generation, propagation, and modulation of seizures. The purpose of this work was to study functional and effective connectivity between regions involved in networks of epileptic seizures. The beginning and middle part of focal seizures from ictal surface EEG data were analyzed using dynamic imaging of coherent sources (DICS), an inverse solution in the frequency domain which describes neuronal networks and coherences of oscillatory brain activities. The information flow (effective connectivity) between coherent sources was investigated using the renormalized partial directed coherence (RPDC) method. In 8/11 patients, the first and second source of epileptic activity as found by DICS were concordant with the operative resection site; these patients became seizure free after epilepsy surgery. In the remaining 3 patients, the results of DICS / RPDC calculations and the resection site were discordant; these patients had a poorer post-operative outcome. The first sources as found by DICS were located predominantly in cortical structures; subsequent sources included some subcortical structures: thalamus, Nucl. Subthalamicus and cerebellum. DICS seems to be a powerful tool to define the seizure onset zone and the epileptic networks involved. Seizure generation seems to be related to the propagation of epileptic activity from the primary source in the seizure onset zone, and maintenance of seizures is attributed to the perpetuation of epileptic activity between nodes in the epileptic network. Despite of these promising results, this proof of principle study needs further confirmation prior to the use of the described methods in the clinical praxis.     

Contributions of astrocytes to epileptogenesis following status epilepticus: Opportunities for preventive therapy?

Status epilepticus (SE) is a life threatening condition that often precedes the development of epilepsy. Traditional treatments for epilepsy have been focused on targeting neuronal mechanisms contributing to hyperexcitability, however, approximately 30% of patients with epilepsy do not respond to existing neurocentric pharmacotherapies. A growing body of evidence has demonstrated that profound changes in the morphology and function of astrocytes accompany SE and persist in epilepsy. Astrocytes are increasingly recognized for their diverse roles in modulating neuronal activity, and understanding the changes in astrocytes following SE could provide important clues about the mechanisms underlying seizure generation and termination. By understanding the contributions of astrocytes to the network changes underlying epileptogenesis and the development of epilepsy, we will gain a greater appreciation of the contributions of astrocytes to dynamic circuit changes, which will enable us to develop more successful therapies to prevent and treat epilepsy. This review summarizes changes in astrocytes following SE in animal models and human temporal lobe epilepsy and addresses the functional consequences of those changes that may provide clues to the process of epileptogenesis.     

Connexin 26 35delG does not represent a mutational hotspot

Recent evidence suggests that the susceptibility to respiratory distress syndrome (RDS) is partly explained by genetic variation in the surfactant proteins (SP) SP-A and SP-B. The present study was designed to evaluate the concordance difference method and candidate gene analysis, in parallel, for the investigation of genetic susceptibility to RDS. We studied 100 same-sex twin pairs with established RDS in at least one twin. The difference in RDS concordance rates between the monozygotic (MZ) and dizygotic (DZ) twin pairs as evidence of a genetic influence was evaluated, and the SP-A and SP-B genes were investigated for potential associations with the susceptibility to RDS. The concordance rates of RDS were 54 and 44% in the MZ and DZ pairs, respectively. The concordance difference of 10% was not significant [95% confidence interval (CI) -0.1 to +0.3, P=0.32], suggesting a low hereditary impact. However, the SP-B Ile131Thr polymorphism was associated with RDS. The threonine allele was associated with an increased risk of RDS [odds ratio (OR) 2.2, 95% CI 1.4-3.5, P=0.0014]. This was particularly apparent in first-born male infants (OR 6.2, 95% CI 2.4-16.3, P<0.001). The degree of prematurity (<32 weeks OR 2.0, 95% CI 1.1-3.7, P=0.021) and birth order (second-born OR 3.1, 95% CI 1.3-7.4, P=0.009) were the clinical variables affecting the risk of RDS. An association between the SP-B Ile131Thr polymorphism and RDS was found. The threonine allele was associated with the risk of RDS, particularly in the first-born twin infants. The concordance difference between MZ and DZ twin pairs underestimates the genetic impact on the risk of RDS. The traditional twin concordance study is insufficient to evaluate genetic predisposition to RDS or other diseases that are confounded by the birth order or multiple pregnancy in itself.