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.

A genome scan for eye color in 502 twin families: Most variation is due to a QTL on chromosome 15q

We have rated eye color on a 3-point scale (1=blue/grey, 2=hazel/green, 3=brown) in 502 twin families and carried out a 5-10 cM genome scan (400-757 markers). We analyzed eye color as a threshold trait and performed multipoint sib pair linkage analysis using variance components analysis in Mx. A lod of 19.2 was found at the marker D15S1002, less than 1 cM from OCA2, which has been previously implicated in eye color variation. We estimate that 74% of variance in eye color liability is due to this QTL and a further 18% due to polygenic effects. However, a large shoulder on this peak suggests that other loci affecting eye color may be telomeric of OCA2 and inflating the QTL estimate. No other peaks reached genome-wide significance, although lods >2 were seen on 5p and 14q and lods >1 were additionally seen on chromosomes 2, 3, 6, 7, 8, 9, 17 and 18. Most of these secondary peaks were reduced or eliminated when we repeated the scan as a two locus analysis with the 15q linkage included, although this does not necessarily exclude them as false positives. We also estimated the interaction between the 15q QTL and the other marker locus but there was only minor evidence for additive x additive epistasis. Elaborating the analysis to the full two-locus model including non-additive main effects and interactions did not strengthen the evidence for epistasis. We conclude that most variation in eye color in Europeans is due to polymorphism in OCA2 but that there may be modifiers at several other loci.

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.

Interactive effects of MC1R and OCA2 on melanoma risk phenotypes

The relationships between MC1R gene variants and red hair, skin reflectance, degree of freckling and nevus count were investigated in 2331 adolescent twins, their sibs and parents in 645 twin families. Penetrance of each MC1R variant allele was consistent with an allelic model where effects were multiplicative for red hair but additive for skin reflectance. Of nine MC1R variant alleles assayed, four common alleles were strongly associated with red hair and fair skin (Asp84Glu, Arg151Cys, Arg160Trp and Asp294His), with a further three alleles having low penetrance (Val60Leu, Val92Met and Arg163Gln). These variants were separately combined for the purposes of this analysis and designated as strong 'R' (OR=63.3; 95% CI 31.9-139.6) and weak 'r ' (OR=5.1; 95% CI 2.5-11.3) red hair alleles. Red-haired individuals are predominantly seen in the R/R and R/r groups with 67.1 and 10.8%, respectively. To assess the interaction of the brown eye color gene OCA2 on the phenotypic effects of variant MC1R alleles we included eye color as a covariate, and also genotyped two OCA2 SNPs (Arg305Trp and Arg419Gln), which were confirmed as modifying eye color. MC1R genotype effects on constitutive skin color, freckling and mole count were modified by eye color, but not genotype for these two OCA2 SNPs. This is probably due to the association of these OCA2 SNPs with brown/green not blue eye color. Amongst individuals with a R/R genotype (but not R/r), those who also had brown eyes had a mole count twice that of those with blue eyes. This suggests that other OCA2 polymorphisms influence mole count and remain to be described.

A deletion mutation in GDF9 in sisters with spontaneous DZ twins

A loss of function mutation in growth differentiation factor 9 (GDF9) in sheep causes increased ovulation rate and infertility in a dosage-sensitive manner. Spontaneous dizygotic (DZ) twinning in the human is under genetic control and women with a history of DZ twinning have an increased incidence of multiple follicle growth and multiple ovulation. We sequenced the GDF9 coding region in DNA samples from 20 women with DZ twins and identified a four-base pair deletion in GDF9 in two sisters with twins from one family. We screened a further 429 families and did not find the loss of function mutation in any other families. We genotyped eight single nucleotide polymorphisms across the GDF9 locus in 379 families with two sisters who have both given birth to spontaneous DZ twins (1527 individuals) and 226 triad families with mothers of twins and their parents (723 individuals). Using case control analysis and the transmission disequilibrium test we found no evidence for association between common variants in GDF9 and twinning in the families. We conclude that rare mutations in GDF9 may influence twinning, but twinning frequency is not associated with common variation in GDF9.

The Evolving Nature of Small-Scale Forestry in Australia

Two forms of small-scale forestry are developing in Australia, each with different impacts on rural communities. One is based on growing short-rotation Eucalyptus globulus (blue gum) for pulp and the other on production of higher-value products from longer-rotation native hard-woods. Several impediments exist to further development of small-scale forestry, including the lack of a small-scale forestry culture, concerns over harvest rights, lack of market development, the long wait for returns, and satisfaction with current land uses. Nevertheless, the rapid increase in farm woodlot establishment in the past five years has paralleled the strong increase in the private industrial plantation estate. As markets develop and hindrances are overcome, landholders not previously interested in small-scale forestry may consider ita worthwhile land use.

Exposing extinction risk analysis to pathogens: Is disease just another form of density dependence?

In the United States and several other countries., the development of population viability analyses (PVA) is a legal requirement of any species survival plan developed for threatened and endangered species. Despite the importance of pathogens in natural populations, little attention has been given to host-pathogen dynamics in PVA. To study the effect of infectious pathogens on extinction risk estimates generated from PVA, we review and synthesize the relevance of host-pathogen dynamics in analyses of extinction risk. We then develop a stochastic, density-dependent host-parasite model to investigate the effects of disease on the persistence of endangered populations. We show that this model converges on a Ricker model of density dependence under a suite of limiting assumptions, including. a high probability that epidemics will arrive and occur. Using this modeling framework, we then quantify: (1) dynamic differences between time series generated by disease and Ricker processes with the same parameters; (2) observed probabilities of quasi-extinction for populations exposed to disease or self-limitation; and (3) bias in probabilities of quasi-extinction estimated by density-independent PVAs when populations experience either form of density dependence. Our results suggest two generalities about the relationships among disease, PVA, and the management of endangered species. First, disease more strongly increases variability in host abundance and, thus, the probability of quasi-extinction, than does self-limitation. This result stems from the fact that the effects and the probability of occurrence of disease are both density dependent. Second, estimates of quasi-extinction are more often overly optimistic for populations experiencing disease than for those subject to self-limitation. Thus, although the results of density-independent PVAs may be relatively robust to some particular assumptions about density dependence, they are less robust when endangered populations are known to be susceptible to disease. If potential management actions involve manipulating pathogens, then it may be useful to. model disease explicitly.

Genetic time-series analysis identifies a major QTL for in vivo alcohol metabolism not predicted by in vitro studies of structural protein polymorphism at the ADH1B or ADH1C loci

After ingestion of a standardized dose of ethanol, alcohol concentrations were assessed, over 3.5 hours from blood (six readings) and breath (10 readings) in a sample of 412 MZ and DZ twins who took part in an Alcohol Challenge Twin Study (ACTS). Nearly all participants were subsequently genotyped on two polymorphic SNPs in the ADH1B and ADH1C loci known to affect in vitro ADH activity. In the DZ pairs, 14 microsatellite markers covering a 20.5 cM region on chromosome 4 that includes the ADH gene family were assessed, Variation in the timed series of autocorrelated blood and breath alcohol readings was studied using a bivariate simplex design. The contribution of a quantitative trait locus (QTL) or QTL's linked to the ADH region was estimated via a mixture of likelihoods weighted by identity-by-descent probabilities. The effects of allelic substitution at the ADH1B and ADH1C loci were estimated in the means part of the model simultaneously with the effects sex and age. There was a major contribution to variance in alcohol metabolism due to a QTL which accounted for about 64% of the additive genetic covariation common to both blood and breath alcohol readings at the first time point. No effects of the ADH1B*47His or ADH1C*349Ile alleles on in vivo metabolism were observed, although these have been shown to have major effects in vitro. This implies that there is a major determinant of variation for in vivo alcohol metabolism in the ADH region that is not accounted for by these polymorphisms. Earlier analyses of these data suggested that alcohol metabolism is related to drinking behavior and imply that this QTL may be protective against alcohol dependence.

Double-strand break repair gene polymorphisms and risk of breast or ovarian cancer

Deficiencies in DNA repair have been hypothesized to increase cancer risk and excess cancer incidence is a feature of inherited diseases caused by defects in DNA damage recognition and repair. We investigated, using a case-control design, whether the double-strand break repair gene polymorphisms RAD51 5' untranslated region -135 G > C, XRCC2 R188H G > A, and XRCC3 T241M C > T were associated with risk of breast or ovarian cancer in Australian women. Sample sets included 1,456 breast cancer cases and 793 age-matched controls ages under 60 years of age, 549 incident ovarian cancer cases, and 335 controls of similar age distribution. For the total sample and the subsample of Caucasian women, there were no significant differences in genotype distribution between breast cancer cases and controls or between ovarian cancer cases and combined control groups. The crude odds ratios (OR) and 95% confidence intervals (95% CI) associated with the RAD51 GC/CC genotype frequency was OR, 1.10; 95% CI, 0.80-1.41 for breast cancer and OR, 1.22; 95% CI, 0.92-1.62 for ovarian cancer. Similarly, there were no increased risks associated with the XRCC2 GA/AA genotype (OR, 0.98; 95% CI, 0.76-1.26 for breast cancer and OR, 0.93; 95% CI, 0.69-1.25 for ovarian cancer) or the XRCC3 CT/TT genotype (OR, 0.92; 95% Cl, 0.77-1.10 for breast cancer and OR, 0.87; 95% CI, 0.71-1.08 for ovarian cancer). Results were little changed after adjustment for age and other measured risk factors. Although there was little statistical power to detect modest increases in risk for the homozygote variant genotypes, particularly for the rare RAD51 and XRCC2 variants, the data suggest that none of these variants play a major role in the etiology of breast or ovarian cancer.

Association between polymorphisms in the progesterone receptor gene and endometriosis

The progesterone receptor (PR) is a candidate gene for the development of endometriosis, a complex disease with strong hormonal features, common in women of reproductive age. We typed the 306 base pair Alu insertion (AluIns) polymorphism in intron G of PR in 101 individuals, estimated linkage disequilibrium (LD) between five single-nucleotide polymorphisms (SNPs) across the PR locus in 980 Australian triads (endometriosis case and two parents) and used transmission disequilibrium testing (TDT) for association with endometriosis. The five SNPs showed strong pairwise LD, and the AluIns was highly correlated with proximal SNPs rs1042839 ({Delta}2 = 0.877, D9 = 1.00, P < 0.0001) and rs500760 ({Delta}2 = 0.438, D9 = 0.942, P < 0.0001). TDT showed weak evidence of allelic association between endometriosis and rs500760 (P = 0.027) but not in the expected direction. We identified a common susceptibility haplotype GGGCA across the five SNPs (P = 0.0167) in the whole sample, but likelihood ratio testing of haplotype transmission and non-transmission of the AluIns and flanking SNPs showed no significant pattern. Further, analysis of our results pooled with those from two previous studies suggested that neither the T2 allele of the AluIns nor the T1/T2 genotype was associated with endometriosis.

New concepts for distinguishing the hidden patterns of linkage disequilibrium which underlie association between genotypes and complex phenotypes

A disappointing feature of conventional methods for detecting association between DNA variation and a phenotype of interest is that they tell us little about the hidden pattern of linkage disequilibrium (LD) with the functional variant that is actually responsible for the association. This limitation applies to case-control studies and also to the transmission/disequilibrium test (TDT) and other family-based association methods. Here we present a fresh perspective on genetic association based on two novel concepts called 'LD squares' and 'equi-risk alleles'. These describe and characterize the different patterns of gametic LD which underlie genetic association. These concepts lead to a general principle - the Equi-Risk Allele Segregation Principle - which captures the way in which underlying LD patterns affect the transmission patterns of genetic variants associated with a phenotype. This provides a basis for distinguishing the hidden LD patterns and might help to locate the functional variants responsible for the association.

Genomewide linkage study in 1,176 affected sister pair families identifies a significant susceptibility locus for endometriosis on chromosome 10q26

Endometriosis is a common gynecological disease that affects up to 10% of women in their reproductive years. It causes pelvic pain, severe dysmenorrhea, and subfertility. The disease is defined as the presence of tissue resembling endometrium in sites outside the uterus. Its cause remains uncertain despite 150 years of hypothesis-driven research, and thus the therapeutic options are limited. Disease predisposition is inherited as a complex genetic trait, which provides an alternative route to understanding the disease. We seek to identify susceptibility loci, using a positional-cloning approach that starts with linkage analysis to identify genomic regions likely to harbor these genes. We conducted a linkage study of 1,176 families ( 931 from an Australian group and 245 from a U. K. group), each with at least two members-mainly affected sister pairs-with surgically diagnosed disease. We have identified a region of significant linkage on chromosome 10q26 ( maximum LOD score [MLS] of 3.09; genomewide P = .047) and another region of suggestive linkage on chromosome 20p13 MLS p 2.09). Minor peaks with MLS > 1.0) were found on chromosomes 2, 6, 7, 8, 12, 14, 15, and 17. This is the first report of linkage to a major locus for endometriosis. The findings will facilitate discovery of novel positional genetic variants that influence the risk of developing this debilitating disease. Greater understanding of the aberrant cellular and molecular mechanisms involved in the etiology and pathophysiology of endometriosis should lead to better diagnostic methods and targeted treatments.

Variance components analyses of multiple asthma traits in a large sample of Australian families ascertained through a twin proband

Background: Intermediate phenotypes are often measured as a proxy for asthma. It is largely unclear to what extent the same set of environmental or genetic factors regulate these traits. Objective: Estimate the environmental and genetic correlations between self-reported and clinical asthma traits. Methods: A total of 3073 subjects from 802 families were ascertained through a twin proband. Traits measured included self-reported asthma, airway histamine responsiveness (AHR), skin prick response to common allergens including house dust mite (Dermatophagoides pteronyssinus [D. pter]), baseline lung function, total serum immunoglobulin E (IgE) and eosinophilia. Bivariate and multivariate analyses of eight traits were performed with adjustment for ascertainment and significant covariates. Results: Overall 2716 participants completed an asthma questionnaire and 2087 were clinically tested, including 1289 self-reported asthmatics (92% previously diagnosed by a doctor). Asthma, AHR, markers of allergic sensitization and eosinophilia had significant environmental correlations with each other (range: 0.23-0.89). Baseline forced expiratory volume in 1 s (FEV1) showed low environmental correlations with most traits. Fewer genetic correlations were significantly different from zero. Phenotypes with greatest genetic similarity were asthma and atopy (0.46), IgE and eosinophilia (0.44), AHR and D. pter (0.43) and AHR and airway obstruction (-0.43). Traits with greatest genetic dissimilarity were FEV1 and atopy (0.05), airway obstruction and IgE (0.07) and FEV1 and D. pter (0.11). Conclusion: These results suggest that the same set of environmental factors regulates the variation of many asthma traits. In addition, although most traits are regulated to great extent by specific genetic factors, there is still some degree of genetic overlap that could be exploited by multivariate linkage approaches.

A genome scan for epidermal skin pattern in adolescent twins reveals suggestive linkage on 12p13.31

Deterioration in stratum corneum reticular patterning (skin pattern or skin wrinkling) has been associated with increased rates of solar keratoses and skin cancer. A previous analysis of data from the twin sample used in this investigation has shown that 86% of the variation in skin pattern is genetic at age 12 and 62% in an adult sample (mean age 47.5). Variation due to genetic influences is likely to be influenced by more than one locus. Here, we present results of a genome-wide linkage scan of skin pattern in adolescent twins and siblings from 428 nuclear twin families. Sib-pair linkage analysis was performed on skin pattern data collected from twins at age 12 (378 informative families) and 14 (316 families). Suggestive linkage was found at marker D12S397 (12p13.31, logarithm of the odds (lod) 1.94), when the effect of the trait locus was modelled to influence the skin pattern equally at both ages 12 and 14. In the same analysis, a peak was seen at 4q23 with a lod score of 1.55. A possible candidate for the peak at 12p13.31 is the protease inhibitor, alpha-2-macroglobulin.

Replicated linkage for eye color on 15q using comparative ratings of sibling pairs

The aim of the study was to perform a genetic linkage analysis for eye color, for comparative data. Similarity in eye color of mono- and dizygotic twins was rated by the twins' mother, their father and/or the twins themselves. For 4748 twin pairs the similarity in eye color was available on a three point scale (not at all alike-somewhat alike-completely alike), absolute eye color on individuals was not assessed. The probability that twins were alike for eye color was calculated as a weighted average of the different responses of all respondents on several different time points. The mean probability of being alike for eye color was 0.98 for MZ twins (2167 pairs), whereas the mean probability for DZ twins was 0.46 (2537 pairs), suggesting very high heritability for eye color. For 294 DZ twin pairs genome-wide marker data were available. The probability of being alike for eye color was regressed on the average amount of IBD sharing. We found a peak LOD-score of 2.9 at chromosome 15q, overlapping with the region recently implicated for absolute ratings of eye color in Australian twins [Zhu, G., Evans, D. M., Duffy, D. L., Montgomery, G. W., Medland, S. E., Gillespie, N. A., Ewen, K. R., Jewell, M., Liew, Y. W., Hayward, N. K., Sturm, R. A., Trent, J. M., and Martin, N. G. (2004). Twin Res. 7:197-210] and containing the OCA2 gene, which is the major candidate gene for eye color [Sturm, R. A. Teasdale, R. D, and Box, N. F. (2001). Gene 277:49-62]. Our results demonstrate that comparative measures on relatives can be used in genetic linkage analysis.