Subtypes of Illicit Drug Users: A Latent Class Analysis of Data From an Australian Twin Sample

This article applies methods of latent class analysis (LCA) to data on lifetime illicit drug use in order to determine whether qualitatively distinct classes of illicit drug users can be identified. Self-report data on lifetime illicit drug use (cannabis, stimulants, hallucinogens, sedatives, inhalants, cocaine, opioids and solvents) collected from a sample of 6265 Australian twins (average age 30 years) were analyzed using LCA. Rates of childhood sexual and physical abuse, lifetime alcohol and tobacco dependence, symptoms of illicit drug abuse/dependence and psychiatric comorbidity were compared across classes using multinomial logistic regression. LCA identified a 5-class model: Class 1 (68.5%) had low risks of the use of all drugs except cannabis; Class 2 (17.8%) had moderate risks of the use of all drugs; Class 3 (6.6%) had high rates of cocaine, other stimulant and hallucinogen use but lower risks for the use of sedatives or opioids. Conversely, Class 4 (3.0%) had relatively low risks of cocaine, other stimulant or hallucinogen use but high rates of sedative and opioid use. Finally, Class 5 (4.2%) had uniformly high probabilities for the use of all drugs. Rates of psychiatric comorbidity were highest in the polydrug class although the sedative/opioid class had elevated rates of depression/suicidal behaviors and exposure to childhood abuse. Aggregation of population-level data may obscure important subgroup differences in patterns of illicit drug use and psychiatric comorbidity. Further exploration of a 'self-medicating' subgroup is needed.

Heritability and Stability of Resting Blood Pressure in Australian Twins

In Australian twins participating in three different studies (1979-1996), the contribution of genetic and environmental influences to variation in resting systolic (SBP) and diastolic blood pressure (DBP) was studied. The sample consisted of 368 monozygotic and 335 dizygotic twin pairs with measurements for both individuals. Blood pressure measurements in two studies were available for 115 complete twin pairs, and 49 twin pairs had measurements in three studies. This allowed assessment of blood pressure tracking over an average period of 12 years in the age range of 23 to 45 years. Multivariate analyses showed significant heritability (h(2)) of blood pressure in all studies (SBP h(2) = 19%-56%, DBP h(2) = 37%-52%). In addition, the analyses showed that the blood pressure tracking was explained by the same set of genetic factors. These results replicate an earlier finding in Dutch twins that also showed stability of the contribution of genetic factors to blood pressure tracking.

Migraine With Aura and Migraine Without Aura Are Not Distinct Entities: Further Evidence From a Large Dutch Population Study

It is often debated whether migraine with aura (MA) and migraine without aura (MO) are etiologically distinct disorders. A previous study using latent class analysis (LCA) in Australian twins showed no evidence for separate subtypes of MO and MA. The aim of the present study was to replicate these results in a population of Dutch twins and their parents, siblings and partners (N = 10,144). Latent class analysis of International Headache Society (IHS)-based migraine symptoms resulted in the identification of 4 classes: a class of unaffected subjects (class 0), a mild form of nonmigrainous headache (class 1), a moderately severe type of migraine (class 2), typically without neurological symptoms or aura (8% reporting aura symptoms), and a severe type of migraine (class 3), typically with neurological symptoms, and aura symptoms in approximately half of the cases. Given the overlap of neurological symptoms and nonmutual exclusivity of aura symptoms, these results do not support the MO and MA subtypes as being etiologically distinct. The heritability in female twins of migraine based on LCA classification was estimated at .50 (95% confidence intervals [0CI} .27 -.59), similar to IHS-based migraine diagnosis (h(2) = .49, 95% Cl .19-.57). However, using a dichotomous classification (affected-unaffected) decreased heritability for the IHS-based classification (h(2) = .33, 95% Cl .00-.60), but not the LCA-based classification (h(2) = .51, 95% Cl. 23-.61). Importantly, use of the LCA-based classification increased the number of subjects classified as affected. The heritability of the screening question was similar to more detailed LCA and IHS classifications, suggesting that the screening procedure is an important determining factor in genetic studies of migraine.

Increased DNA Methylation at the AXIN1 Gene in a Monozygotic Twin from a Pair Discordant for a Caudal Duplication Anomaly

The AXIN1 gene has been implicated in caudal duplication anomalies. Its coding region was sequenced in both members of a monozygotic ( MZ) twin pair discordant for a caudal duplication anomaly, but no mutation was found. Using bisulfite sequencing, we examined methylation at the promoter region of the AXIN1 gene in these twins and in twin and age-matched singleton controls. Methylation of the promoter region in peripheral blood mononucleated cells was variable among individuals, including MZ pairs. In the MZ pair discordant for the caudal duplication, this region of the affected twin was significantly more methylated than that of the unaffected twin (), which was significantly more P < .0001 methylated than those of the controls (). We have confirmed that this CpG island does function as a promoter P = .02 in vitro and that its activity is inversely proportional to the extent of methylation. This finding raises the possibility that hypermethylation of the AXIN1 promoter, by mechanisms as yet undetermined, is associated with the malformation. This case may be paradigmatic for some cases of MZ discordance.

Population differences in finger-length ratios: Ethnicity or latitude?

The relative length of the second and fourth fingers (the 2D:4D ratio) has been taken to be an indicator of prenatal exposure to testosterone, and hence possibly relevant to sexual orientation and other sex-differentiated behaviors. Studies have reported a difference in this ratio between Caucasian males in Britain and in the U.S.: higher average 2D:4D ratios were obtained in Britain. This raises the question of whether differences among different Caucasian gene pools were responsible or whether some environmental variable associated with latitude might be involved (e.g., exposure to sunlight or different day-length patterns). This question was explored by examining 2D:4D ratios for an Australian adolescent sample. The Australians were predominantly of British ancestry, but lived at distances from the equator more like those of the U.S. studies. The Australian 2D:4D ratios resembled those in Britain rather than those in the U.S., tending to exclude hypotheses related to latitude and making differences in gene pools a plausible explanation.

Genome-wide linkage scan for loci influencing plasma triglycerides

Background: Plasma triglyceride concentration is known to be a significant risk factor for cardiovascular disease (CVD). Previous studies have found that the level of triglycerides is strongly influenced by genetic factors. Methods: To identify quantitative trait loci influencing triglycerides, we conducted a genome-wide linkage scan on data from 485 Australian adult dizygotic twin pairs. Prior to linkage analysis, triglyceride values were adjusted for the effects of covariates including age, sex, time since last meal, time of blood collection (CT) and time to plasma separation. Results: The heritability estimate for ln(triglyceride) adjusted for all above fixed effects was 0.49. The highest multipoint LOD score observed was 2.94 (genome-wide p=0.049) on chromosome 7 (at 65cM). This 7p region contains several candidate genes. Two other regions with suggestive multipoint LOD scores were also identified on chromosome 4 (LOD score=2.26 at 62cM) and chromosome X (LOD score=2.01 at 81cM). Conclusions: The linkage peaks found represent newly identified regions for more detailed study, in particular the significant linkage observed on chromosome 7p13. \ (c) 2006 Elsevier B.V. All rights reserved.

A simple method to localise pleiotropic susceptibility loci using univariate linkage analyses of correlated traits

Univariate linkage analysis is used routinely to localise genes for human complex traits. Often, many traits are analysed but the significance of linkage for each trait is not corrected for multiple trait testing, which increases the experiment-wise type-I error rate. In addition, univariate analyses do not realise the full power provided by multivariate data sets. Multivariate linkage is the ideal solution but it is computationally intensive, so genome-wide analysis and evaluation of empirical significance are often prohibitive. We describe two simple methods that efficiently alleviate these caveats by combining P-values from multiple univariate linkage analyses. The first method estimates empirical pointwise and genome-wide significance between one trait and one marker when multiple traits have been tested. It is as robust as an appropriate Bonferroni adjustment, with the advantage that no assumptions are required about the number of independent tests performed. The second method estimates the significance of linkage between multiple traits and one marker and, therefore, it can be used to localise regions that harbour pleiotropic quantitative trait loci (QTL). We show that this method has greater power than individual univariate analyses to detect a pleiotropic QTL across different situations. In addition, when traits are moderately correlated and the QTL influences all traits, it can outperform formal multivariate VC analysis. This approach is computationally feasible for any number of traits and was not affected by the residual correlation between traits. We illustrate the utility of our approach with a genome scan of three asthma traits measured in families with a twin proband.

Rapid screening of 4000 individuals for germ-line variations in the BRAF gene

Background: The BRAF gene is frequently somatically altered in malignant melanoma. A majority of variations are at the valine 600 residue leading to a V600E substitution that constitutively activates the kinase. We screened 4000 patient and control DNAs for germ-line variations at the valine 600 residue. Methods: We developed a novel assay by adapting single-base variation assays and software for MALDI-TOF (matrix-assisted laser desorption/ionization time-of-flight) mass spectrometry to screen for all 5 reported variants at codon 600 of the BRAF gene. We screened a case-control collection comprising samples from 1082 melanoma patients and 154 of their unaffected relatives from 1278 families and from 2744 individuals from 659 unselected twin families with no history of melanoma. A panel of 66 melanoma cell lines was used for variation-positive controls. Results: All melanoma cell lines that we had found previously to carry a codon 600 variation were verified in this study. Three of the 4 possible variants (V600E n = 47, V600K n = 2, V600R n = 1) were detected, but no case of V600D was available. No germ-line variants were found in the samples from the 3980 melanoma patients or from the control individuals. Conclusions: This new assay is a high-throughput, automated alternative to standard sequencing and can be used as a rapid initial screen for somatic variants associated with melanoma. Germ-line variants at valine 600 are unlikely to exist and do not contribute to the reported role of the BRAF gene in melanoma predisposition. (c) 2006 American Association for Clinical Chemistry.

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%.

Novel variants in growth differentiation factor 9 in mothers of dizygotic twins

Context: Genes from the ovarian bone morphogenetic signaling pathway (GDF9 and BMP15) are critical for normal human fertility. We previously identified a deletion mutation in GDF9 in sisters with spontaneous dizygotic (DZ) twins, but the prevalence of rare GDF9 variants in twinning families is unknown. Objective: The objective was to evaluate the frequency of rare variants in GDF9 in families with a history of DZ twinning. Design and Subjects: We recruited 3450 individuals from 915 DZ twinning families (1693 mothers of twins) and 1512 controls of Caucasian origin. One mother of DZ twins was selected from 279 of the 915 families, and a DNA sample was screened for rare variants in GDF9 using denaturant HPLC. Variants were confirmed by DNA sequencing and genotyped in the entire sample by matrix-assisted laser desorption ionization time of flight (MALDI-TOF) mass spectrometry. Results: We found two novel insertion/deletions (c.392-393insT, c.1268-1269delAA) and four missense alterations in the GDF9 sequence in mothers of twins. Two of the missense variants (c.307C > T, p.Pro103Ser and c.362C > T, p.Thr121Leu) were located in the proregion of GDF9 and two (c.1121C > T, p.Pro374Leu and c.1360C > T, p.Arg454Cys) in the mature protein region. For each variant, the frequencies were higher in cases compared with controls. The proportion of mothers of DZ twins carrying any variant (4.12%) was significantly higher (P < 0.0001) than the proportion of carriers in controls (2.29%). Conclusion: We describe new variants in the GDF9 gene that are significantly more common in mothers of DZ twins than controls, suggesting that rare GDF9 variants contribute to the likelihood of DZ twinning.

HLA and Genomewide Allele Sharing in Dizygotic Twins

Gametic selection during fertilization or the effects of specific genotypes on the viability of embryos may cause a skewed transmission of chromosomes to surviving offspring. A recent analysis of transmission distortion in humans reported significant excess sharing among full siblings. Dizygotic (DZ) twin pairs are a special case of the simultaneous survival of two genotypes, and there have been reports of DZ pairs with excess allele sharing around the HLA locus, a candidate locus for embryo survival. We performed an allele-sharing study of 1,592 DZ twin pairs from two independent Australian cohorts, of which 1,561 pairs were informative for linkage on chromosome 6. We also analyzed allele sharing in 336 DZ twin pairs from The Netherlands. We found no evidence of excess allele sharing, either at the HLA locus or in the rest of the genome. In contrast, we found evidence of a small but significant (P = .003 for the Australian sample) genomewide deficit in the proportion of two alleles shared identical by descent among DZ twin pairs. We reconciled conflicting evidence in the literature for excess genomewide allele sharing by performing a simulation study that shows how undetected genotyping errors can lead to an apparent deficit or excess of allele sharing among sibling pairs, dependent on whether parental genotypes are known. Our results imply that gene-mapping studies based on affected sibling pairs that include DZ pairs will not suffer from false-positive results due to loci involved in embryo survival.

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

KRAS variation and risk of endometriosis

Endometriosis is a common gynaecological disease with symptoms of pelvic pain and infertility which affects 7-10% of women in their reproductive years. Activation of an oncogenic allele of Kirsten rat sarcoma viral oncogene homologue (KRAS) in the reproductive tract of mice resulted in the development of endometriosis. We hypothesized that variation in KRAS may influence risk of endometriosis in humans. Thirty tagSNPs spanning a region of 60.7 kb across the KRAS locus were genotyped using iPLEX chemistry on a MALDI-TOF MassARRAY platform in 959 endometriosis cases and 959 unrelated controls, and data were analysed for association with endometriosis. Genotypes were obtained for most individuals with a mean completion rate of 99.1%. We identified six haplotype blocks across the KRAS locus in our sample. There were no significant differences between cases and controls in the frequencies of individual single-nucleotide polymorphisms (SNPs) or haplotypes. We also developed a rapid method to screen for 11 common KRAS and BRAF mutations on the Sequenom MassARRAY system. The assay detected all mutations previously identified by direct sequencing in a panel of positive controls. No germline variants for KRAS or BRAF were detected. Our results demonstrate that any risk of endometriosis in women because of common variation in KRAS must be very small.

Genetic and Cultural Transmission of Smoking Initiation: An Extended Twin Kinship Model

Background Considerable evidence from twin and adoption studies indicates that genetic and shared environmental factors play a significant role in the initiation of smoking behavior. Although twin and adoption designs are powerful to detect genetic and environmental influences, they do not provide information on the processes of assortative mating and parent–offspring transmission and their contribution to the variability explained by genetic and/or environmental factors. Methods We examined the role of genetic and environmental factors for smoking initiation using an extended kinship design. This design allows the simultaneous testing of additive and non-additive genetic, shared and individual-specific environmental factors, as well as sex differences in the expression of genes and environment in the presence of assortative mating and combined genetic and cultural transmission. A dichotomous lifetime smoking measure was obtained from twins and relatives in the Virginia 30,000 sample. Results Results demonstrate that both genetic and environmental factors play a significant role in the liability to smoking initiation. Major influences on individual differences appeared to be additive genetic and unique environmental effects, with smaller contributions from assortative mating, shared sibling environment, twin environment, cultural transmission and resulting genotype–environment covariance. The finding of negative cultural transmission without dominance led us to investigate more closely two possible mechanisms for the lower parent–offspring correlations compared to the sibling and DZ twin correlations in subsets of the data: (i) age × gene interaction, and (ii) social homogamy. Neither mechanism provided a significantly better explanation of the data, although age regression was significant. Conclusions This study showed significant heritability, partly due to assortment, and significant effects of primarily non-parental shared environment on smoking initiation.