Christopher J. Hammond

Human metabolic individuality in biomedical and pharmaceutical research

Genome-wide association studies (GWAS) have identified many risk loci for complex diseases, but effect sizes are typically small and information on the underlying biological processes is often lacking. Associations with metabolic traits as functional intermediates can overcome these problems and potentially inform individualized therapy. Here we report a comprehensive analysis of genotype-dependent metabolic phenotypes using a GWAS with non-targeted metabolomics. We identified 37 genetic loci associated with blood metabolite concentrations, of which 25 show effect sizes that are unusually high for GWAS and account for 10–60% differences in metabolite levels per allele copy. Our associations provide new functional insights for many disease-related associations that have been reported in previous studies, including those for cardiovascular and kidney disorders, type 2 diabetes, cancer, gout, venous thromboembolism and Crohn’s disease. The study advances our knowledge of the genetic basis of metabolic individuality in humans and generates many new hypotheses for biomedical and pharmaceutical research.

The architecture of gene regulatory variation across multiple human tissues: the MuTHER study.

While there have been studies exploring regulatory variation in one or more tissues, the complexity of tissue-specificity in multiple primary tissues is not yet well understood. We explore in depth the role of cis-regulatory variation in three human tissues: lymphoblastoid cell lines (LCL), skin, and fat. The samples (156 LCL, 160 skin, 166 fat) were derived simultaneously from a subset of well-phenotyped healthy female twins of the MuTHER resource. We discover an abundance of cis-eQTLs in each tissue similar to previous estimates (858 or 4.7% of genes). In addition, we apply factor analysis (FA) to remove effects of latent variables, thus more than doubling the number of our discoveries (1,822 eQTL genes). The unique study design (Matched Co-Twin Analysis—MCTA) permits immediate replication of eQTLs using co-twins (93%–98%) and validation of the considerable gain in eQTL discovery after FA correction. We highlight the challenges of comparing eQTLs between tissues. After verifying previous significance threshold-based estimates of tissue-specificity, we show their limitations given their dependency on statistical power. We propose that continuous estimates of the proportion of tissue-shared signals and direct comparison of the magnitude of effect on the fold change in expression are essential properties that jointly provide a biologically realistic view of tissue-specificity. Under this framework we demonstrate that 30% of eQTLs are shared among the three tissues studied, while another 29% appear exclusively tissue-specific. However, even among the shared eQTLs, a substantial proportion (10%–20%) have significant differences in the magnitude of fold change between genotypic classes across tissues. Our results underline the need to account for the complexity of eQTL tissue-specificity in an effort to assess consequences of such variants for complex traits.

Genome-wide association study identifies 12 new susceptibility loci for primary biliary cirrhosis.

In addition to the HLA locus, six genetic risk factors for primary biliary cirrhosis (PBC) have been identified in recent genome-wide association studies (GWAS). To identify additional loci, we carried out a GWAS using 1,840 cases from the UK PBC Consortium and 5,163 UK population controls as part of the Wellcome Trust Case Control Consortium 3 (WTCCC3). We followed up 28 loci in an additional UK cohort of 620 PBC cases and 2,514 population controls. We identified 12 new susceptibility loci (at a genome-wide significance level of P < 5 × 10−8) and replicated all previously associated loci. We identified three further new loci in a meta-analysis of data from our study and previously published GWAS results. New candidate genes include STAT4, DENND1B, CD80, IL7R, CXCR5, TNFRSF1A, CLEC16A and NFKB1. This study has considerably expanded our knowledge of the genetic architecture of PBC.

Common variants near CAV1 and CAV2 are associated with primary open-angle glaucoma

We conducted a genome-wide association study for primary open-angle glaucoma (POAG) in 1,263 affected individuals (cases) and 34,877 controls from Iceland. We identified a common sequence variant at 7q31 (rs4236601[A], odds ratio (OR) = 1.36, P = 5.0 × 10−10). We then replicated the association in sample sets of 2,175 POAG cases and 2,064 controls from Sweden, the UK and Australia (combined OR = 1.18, P = 0.0015) and in 299 POAG cases and 580 unaffected controls from Hong Kong and Shantou, China (combined OR = 5.42, P = 0.0021). The risk variant identified here is located close to CAV1 and CAV2, both of which are expressed in the trabecular meshwork and retinal ganglion cells that are involved in the pathogenesis of POAG.

Genetic influence on early age-related maculopathy: A twin study

OBJECTIVE Age-related macular degeneration (AMD) is the most common cause of blindness in industrialized countries. There has been considerable interest in the genetics of early age-related maculopathy (ARM) and AMD, because they have phenotypes similar to inherited diseases where mutations have been identified, but the heritability of ARM and AMD is unknown. DESIGN A classical twin study was performed to compare the concordance in monozygotic (MZ) and dizygotic (DZ) twins in an unselected sample of female volunteer twins. PARTICIPANTS Five hundred six twin pairs, 226 MZ and 280 DZ, with a mean age of 62 years, were examined. METHODS ARM was graded from stereoscopic macular photographs of 501 of the twin pairs (99%) according to the International ARM Epidemiologic Study Group grading system. The casewise concordance was calculated for twin pairs from 2 x 2 contingency tables of affected/unaffected twins, and these tables were used in maximum likelihood genetic modeling to estimate the heritabilities of phenotypes graded. MAIN OUTCOME MEASURES Prevalence of ARM; concordance in MZ and DZ twins of the phenotypes of ARM, soft drusen >63 microm and > or =125 microm diameter, pigmentary changes and hard drusen (<20 and > or =20 in number); heritability of ARM and subphenotypes. RESULTS The overall prevalence of ARM was 14.6% (95% confidence interval [CI], 12.4%-16.8%). The concordance for ARM in MZ twins was 0.37 compared with 0.19 in DZ twins, suggesting a role for genes. Modeling confirmed a genetic effect for phenotypes of ARM, soft drusen, pigmentary changes, and > or =20 hard drusen, although there was little genetic effect for scattered (<20) hard drusen. The heritability of ARM was estimated as 45% (95% CI, 35%-53%). The most heritable phenotypes were soft drusen > or =125 microm (57%) and > or =20 hard drusen (81%), with the latter being dominantly inherited. CONCLUSIONS This study confirms a significant genetic influence in ARM and suggests that future genetic studies should examine phenotypes of large (> or =125 microm) soft drusen and > or =20 hard drusen, because these seem to be the most heritable components.

Genetic and environmental factors in age-related nuclear cataracts in monozygotic and dizygotic twins

BACKGROUND Age-related cataracts are a major public health problem. The relative importance of genes and environment in the causation of nuclear cataracts, the most common form of age-related cataracts, is not known. METHODS We studied 506 pairs of female twins (226 monozygotic and 280 dizygotic) who were 50 to 79 years old (mean, 62). The amount of nuclear cataract in the right and left eyes was determined objectively by analysis of Scheimpflug lens photographs (yielding three measures) and subjectively with use of the Oxford Clinical Cataract Classification and Grading System (yielding one measure). All eight measures (four in each eye) were subsequently combined in one summary measure of nuclear cataract for each woman. A univariate maximum-likelihood model was used to estimate the variance of the genetic and environmental contributions to each of the measures. RESULTS The different measures of cataract formation were highly correlated (correlation coefficients, 0.71 to 0.94). The mean scores were similar for the right and left eyes and for monozygotic and dizygotic twins. Quantitative genetic modeling of each of the nuclear-cataract scores invariably resulted in a best-fitting model that involved additive genetic effects, unique environmental effects, and age. The common environmental and dominant genetic effects could be removed from the models without significant loss of fit. The overall heritability in the combined nuclear-cataract score (the proportion of the variance explained by genetic factors) was 48 percent (95 percent confidence interval, 42 to 54 percent); age accounted for 38 percent of the variance (95 percent confidence interval, 31 to 44 percent) and unique environmental effects for 14 percent (95 percent confidence interval, 12 to 18 percent). CONCLUSIONS Genetic effects are important even in such a clearly age-related disease as nuclear cataract, explaining almost 50 percent of the variation in the severity of this disease.

A Susceptibility Locus for Myopia in the Normal Population Is Linked to the PAX6 Gene Region on Chromosome 11: A Genomewide Scan of Dizygotic Twins

Myopia is a common, complex trait with considerable economic and social impact and, in highly affected individuals, ocular morbidity. We performed a classic twin study of 506 unselected twin pairs and inferred the heritability of refractive error to be 0.89 (95% confidence interval 0.86-0.91). A genomewide scan of 221 dizygotic twin pairs, analyzed by use of optimal Haseman-Elston regression methods implemented by use of generalized linear modeling, showed significant linkage (LOD >3.2) to refractive error at four loci, with a maximum LOD score of 6.1 at 40 cM on chromome 11p13. Evidence of linkage at this locus, as well as at the other linkage peaks at chromosomes 3q26 (LOD 3.7), 8p23 (LOD 4.1), and 4q12 (LOD 3.3), remained the same or became stronger after model fit was checked and outliers were downweighted. Examination of potential candidate genes showed the PAX6 gene directly below the highest peak at the 11p13 locus. PAX6 is fundamental to identity and growth of the eye, but reported mutations usually result in catastrophic congenital phenotypes such as aniridia. Haplotype tagging of 17 single-nucleotide polymorphisms (SNPs), which covered the PAX6 gene and had common minor allele frequencies, identified 5 SNPs that explained 0.999 of the haplotype diversity. Linkage and association analysis of the tagging SNPs showed strong evidence of linkage for all markers with a minimum chi 21 of 7.5 (P=.006) but no association. This suggests that PAX6 may play a role in myopia development, possibly because of genetic variation in an upstream promoter or regulator, although no definite association between PAX6 common variants and myopia was demonstrated in this study.

Genome-wide association analyses identify multiple loci associated with central corneal thickness and keratoconus

Central corneal thickness (CCT) is associated with eye conditions including keratoconus and glaucoma. We performed a meta-analysis on >20,000 individuals in European and Asian populations that identified 16 new loci associated with CCT at genome-wide significance (P < 5 × 10−8). We further showed that 2 CCT-associated loci, FOXO1 and FNDC3B, conferred relatively large risks for keratoconus in 2 cohorts with 874 cases and 6,085 controls (rs2721051 near FOXO1 had odds ratio (OR) = 1.62, 95% confidence interval (CI) = 1.4–1.88, P = 2.7 × 10−10, and rs4894535 in FNDC3B had OR = 1.47, 95% CI = 1.29–1.68, P = 4.9 × 10−9). FNDC3B was also associated with primary open-angle glaucoma (P = 5.6 × 10−4; tested in 3 cohorts with 2,979 cases and 7,399 controls). Further analyses implicate the collagen and extracellular matrix pathways in the regulation of CCT.

A Genome-Wide Association Study of Optic Disc Parameters

The optic nerve head is involved in many ophthalmic disorders, including common diseases such as myopia and open-angle glaucoma. Two of the most important parameters are the size of the optic disc area and the vertical cup-disc ratio (VCDR). Both are highly heritable but genetically largely undetermined. We performed a meta-analysis of genome-wide association (GWA) data to identify genetic variants associated with optic disc area and VCDR. The gene discovery included 7,360 unrelated individuals from the population-based Rotterdam Study I and Rotterdam Study II cohorts. These cohorts revealed two genome-wide significant loci for optic disc area, rs1192415 on chromosome 1p22 (p = 6.72×10−19) within 117 kb of the CDC7 gene and rs1900004 on chromosome 10q21.3-q22.1 (p = 2.67×10−33) within 10 kb of the ATOH7 gene. They revealed two genome-wide significant loci for VCDR, rs1063192 on chromosome 9p21 (p = 6.15×10−11) in the CDKN2B gene and rs10483727 on chromosome 14q22.3-q23 (p = 2.93×10−10) within 40 kbp of the SIX1 gene. Findings were replicated in two independent Dutch cohorts (Rotterdam Study III and Erasmus Rucphen Family study; N = 3,612), and the TwinsUK cohort (N = 843). Meta-analysis with the replication cohorts confirmed the four loci and revealed a third locus at 16q12.1 associated with optic disc area, and four other loci at 11q13, 13q13, 17q23 (borderline significant), and 22q12.1 for VCDR. ATOH7 was also associated with VCDR independent of optic disc area. Three of the loci were marginally associated with open-angle glaucoma. The protein pathways in which the loci of optic disc area are involved overlap with those identified for VCDR, suggesting a common genetic origin.

Cohort Profile: TwinsUK and Healthy Ageing Twin Study

The UKs largest registry of adult twins, or TwinsUK Registry, started in 1992 and encompasses about 12000 volunteer twins from all over the United Kingdom. More than 70% of the registered twins have filled at least one detailed health questionnaire and about half of them undergone a baseline comprehensive assessment and two follow-up clinical evaluations. The most recent follow-up visit, known as Healthy Ageing Twin Study (HATS), involved 3125 female twins aged >40 years with at least one previous clinical assessment to enable inspection of longitudinal changes in ageing traits and their genetic and environmental components. The study benefits from several state-of-the-art OMICs studies including genome-wide association, next-generation genome and transcriptome sequencing, and epigenetic and metabolomic profiles. This makes our cohort as one of the most deeply phenotyped and genotyped in the world. Several collaborative projects in the field of epidemiology of complex disorders are ongoing in our cohort and interested researchers are encouraged to get in contact for future collaborations.