Gu Zhu

Assumption-Free Estimation of Heritability from Genome-Wide Identity-by-Descent Sharing between Full Siblings

The study of continuously varying, quantitative traits is important in evolutionary biology, agriculture, and medicine. Variation in such traits is attributable to many, possibly interacting, genes whose expression may be sensitive to the environment, which makes their dissection into underlying causative factors difficult. An important population parameter for quantitative traits is heritability, the proportion of total variance that is due to genetic factors. Response to artificial and natural selection and the degree of resemblance between relatives are all a function of this parameter. Following the classic paper by R. A. Fisher in 1918, the estimation of additive and dominance genetic variance and heritability in populations is based upon the expected proportion of genes shared between different types of relatives, and explicit, often controversial and untestable models of genetic and non-genetic causes of family resemblance. With genome-wide coverage of genetic markers it is now possible to estimate such parameters solely within families using the actual degree of identity-by-descent sharing between relatives. Using genome scans on 4,401 quasi-independent sib pairs of which 3,375 pairs had phenotypes, we estimated the heritability of height from empirical genome-wide identity-by-descent sharing, which varied from 0.374 to 0.617 (mean 0.498, standard deviation 0.036). The variance in identity-by-descent sharing per chromosome and per genome was consistent with theory. The maximum likelihood estimate of the heritability for height was 0.80 with no evidence for non-genetic causes of sib resemblance, consistent with results from independent twin and family studies but using an entirely separate source of information. Our application shows that it is feasible to estimate genetic variance solely from within-family segregation and provides an independent validation of previously untestable assumptions. Given sufficient data, our new paradigm will allow the estimation of genetic variation for disease susceptibility and quantitative traits that is free from confounding with non-genetic factors and will allow partitioning of genetic variation into additive and non-additive components.

A Major Quantitative-Trait Locus for Mole Density Is Linked to the Familial Melanoma Gene CDKN2A: A Maximum-Likelihood Combined Linkage and Association Analysis in Twins and Their Sibs

Important risk factors for melanoma are densely clustered melanocytic nevi (common moles) and mutations in the p16 (CDKN2A) gene. Nevi may be subclassified as raised or flat. In our sample, raised nevi were 27% of the total, and the two kinds had a correlation of.33. Correlations for total-nevus count (TNC) in 153 MZ and 199 DZ twin pairs were.94 and.60, respectively, which are compatible with a very-high degree of genetic determination. We hypothesized that some of the genetic variance might be due to variation in the p16 gene. Analysis of linkage to a highly polymorphic marker (D9S942), located close to p16, detected quantitative-trait-loci (QTL) effects accounting for 27% of variance in TNC, rising to 33% if flat but not raised moles were considered. Total heritability was higher for raised (.69) than for flat (.42) moles, but QTL linkage was 0 for raised moles, whereas it accounted for 80% of the heritability of flat moles; additionally, family environment accounted for only 15% of variance in raised versus 46% in flat moles. These findings suggest that raised and flat nevi have very different etiologies. Longer alleles at D9S942 were associated with higher flat-mole counts, and a novel modification to a within-sibship association test showed that this association is genuine and not due to population stratification, although it accounts for only 1% of total variance. Since germline mutations in the exons of CDKN2A are rare, it is likely that variants in the noncoding regions of this gene, or in another gene nearby, are responsible for this major determinant of moliness and, hence, of melanoma risk.

Common Variants in the Trichohyalin Gene Are Associated with Straight Hair in Europeans

Hair morphology is highly differentiated between populations and among people of European ancestry. Whereas hair morphology in East Asian populations has been studied extensively, relatively little is known about the genetics of this trait in Europeans. We performed a genome-wide association scan for hair morphology (straight, wavy, curly) in three Australian samples of European descent. All three samples showed evidence of association implicating the Trichohyalin gene (TCHH), which is expressed in the developing inner root sheath of the hair follicle, and explaining approximately 6% of variance (p=1.5x10(-31)). These variants are at their highest frequency in Northern Europeans, paralleling the distribution of the straight-hair EDAR variant in Asian populations.

A Genome-Wide Association Study Identifies Five Loci Influencing Facial Morphology in Europeans

Inter-individual variation in facial shape is one of the most noticeable phenotypes in humans, and it is clearly under genetic regulation; however, almost nothing is known about the genetic basis of normal human facial morphology. We therefore conducted a genome-wide association study for facial shape phenotypes in multiple discovery and replication cohorts, considering almost ten thousand individuals of European descent from several countries. Phenotyping of facial shape features was based on landmark data obtained from three-dimensional head magnetic resonance images (MRIs) and two-dimensional portrait images. We identified five independent genetic loci associated with different facial phenotypes, suggesting the involvement of five candidate genes—PRDM16, PAX3, TP63, C5orf50, and COL17A1—in the determination of the human face. Three of them have been implicated previously in vertebrate craniofacial development and disease, and the remaining two genes potentially represent novel players in the molecular networks governing facial development. Our finding at PAX3 influencing the position of the nasion replicates a recent GWAS of facial features. In addition to the reported GWA findings, we established links between common DNA variants previously associated with NSCL/P at 2p21, 8q24, 13q31, and 17q22 and normal facial-shape variations based on a candidate gene approach. Overall our study implies that DNA variants in genes essential for craniofacial development contribute with relatively small effect size to the spectrum of normal variation in human facial morphology. This observation has important consequences for future studies aiming to identify more genes involved in the human facial morphology, as well as for potential applications of DNA prediction of facial shape such as in future forensic applications.

Effects of HFE C282Y and H63D Polymorphisms and Polygenic Background on Iron Stores in a Large Community Sample of Twins

The aim of this study was to assess and to compare the role of HFE polymorphisms and other genetic factors in variation in iron stores. Blood samples were obtained from 3,375 adult male and female twins (age range 29-82 years) recruited from the Australian Twin Registry. There were 1,233 complete pairs (562 monozygotic and 571 dizygotic twins). Serum iron, transferrin, transferrin saturation with iron, and ferritin were measured, and the HFE C282Y and H63D genotypes were determined. The frequency of the C282Y allele was.072, and that of the H63D allele was.141. Significant sources of variation in the indices of iron status included age, sex, age-sex interaction, body-mass index, and both the C282Y and H63D genotypes. The iron, transferrin, and saturation values of CC and CY subjects differed significantly, but the ferritin values did not. After correction for age and body-mass index, 23% and 31% of the variance in iron, 66% and 49% of the variance in transferrin, 33% and 47% of the variance in transferrin saturation, and 47% and 47% of the variance in ferritin could be explained by additive genetic factors, for men and women, respectively. HFE C282Y and H63D variation accounted for <5% of the corrected phenotypic variance, except for saturation (12% in women and 5% in men). We conclude that HFE CY and HD heterozygotes differ in iron status from the CC and HH homozygotes and that serum transferrin saturation is more affected than is serum ferritin. There are highly significant effects of other as-yet-unidentified genes on iron stores, in addition to HFE genotype.

A Study of the Genetic and Environmental Etiology of Stuttering in a Selected Twin Sample

Stuttering is a developmental disorder of speech production that usually emerges in childhood. In this study, a large population-based twin sample from the Australian Twin Registry (1567 pairs and 634 singles aged 17–29 years) was screened to identify twin pairs in which one or both members reported themselves to be affected by stuttering. Telephone interview-based diagnoses were obtained for 457 of these individuals (self-reported affected cases, cotwins, and controls) to determine whether the self-report was correct. To correct for ascertainment bias we carried out a bivariate analysis of the final diagnosis in the selected sample with the screening item in the full sample, using the categorical raw data option of Mx 1.47c. After correcting for ascertainment bias, approximately 70% (95% confidence interval: 39–86%) of the variance in liability to stuttering was found to be attributable to additive genetic effects, with the remainder due to nonshared environmental effects.

Digital Quantification of Human Eye Color Highlights Genetic Association of Three New Loci

Previous studies have successfully identified genetic variants in several genes associated with human iris (eye) color; however, they all used simplified categorical trait information. Here, we quantified continuous eye color variation into hue and saturation values using high-resolution digital full-eye photographs and conducted a genome-wide association study on 5,951 Dutch Europeans from the Rotterdam Study. Three new regions, 1q42.3, 17q25.3, and 21q22.13, were highlighted meeting the criterion for genome-wide statistically significant association. The latter two loci were replicated in 2,261 individuals from the UK and in 1,282 from Australia. The LYST gene at 1q42.3 and the DSCR9 gene at 21q22.13 serve as promising functional candidates. A model for predicting quantitative eye colors explained over 50% of trait variance in the Rotterdam Study. Over all our data exemplify that fine phenotyping is a useful strategy for finding genes involved in human complex traits.

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.

Genome Partitioning of Genetic Variation for Height from 11,214 Sibling Pairs

Height has been used for more than a century as a model by which to understand quantitative genetic variation in humans. We report that the entire genome appears to contribute to its additive genetic variance. We used genotypes and phenotypes of 11,214 sibling pairs from three countries to partition additive genetic variance across the genome. Using genome scans to estimate the proportion of the genomes of each chromosome from siblings that were identical by descent, we estimated the heritability of height contributed by each of the 22 autosomes and the X chromosome. We show that additive genetic variance is spread across multiple chromosomes and that at least six chromosomes (i.e., 3, 4, 8, 15, 17, and 18) are responsible for the observed variation. Indeed, the data are not inconsistent with a uniform spread of trait loci throughout the genome. Our estimate of the variance explained by a chromosome is correlated with the number of times suggestive or significant linkage with height has been reported for that chromosome. Variance due to dominance was not significant but was difficult to assess because of the high sampling correlation between additive and dominance components. Results were consistent with the absence of any large between-chromosome epistatic effects. Notwithstanding the proposed architecture of complex traits that involves widespread gene-gene and gene-environment interactions, our results suggest that variation in height in humans can be explained by many loci distributed over all autosomes, with an additive mode of gene action.

The genetics of coronary heart disease: The contribution of twin studies

Despite the decline in coronary heart disease in many European countries, the disease remains an enormous public health problem. Although we know a great deal about environmental risk factors for coronary heart disease, a heritable component was recognized a long time ago. The earliest and best known examples of how our genetic constitution may determine cardiovascular risk relate to lipoprotein(a), familial hypercholesterolaemia and apolipoprotein E. In the past 20 years a fair number of polymorphisms assessed singly have shown strong associations with the disease but most are subject to poor repeatability. Twins constitute a compelling natural experiment to establish the genetic contribution to coronary heart disease and its risk factors. GenomEUtwin, a recently funded Framework 5 Programme of the European Community, affords the opportunity of comparing the heritability of risk factors in different European Twin Registries. As an illustration we present the heritabilities of systolic and diastolic blood pressure, based on data from over 4000 twin pairs from six different European countries and Australia. Heritabilities for systolic blood pressure are between 52 and 66% and for diastolic blood pressure between 44 and 66%. There is no evidence of sex differences in heritability estimates and very little to no evidence for a significant contribution of shared family environment. A non-twin based prospective case/cohort study of coronary heart disease and stroke (MORGAM) will allow hypotheses relating to cardiovascular disease, generated in the twin cohorts, to be tested prospectively in adult populations. Twin studies have also contributed to our understanding of the life course hypothesis, and GenomEUtwin has the potential to add to this.