Faculty Bibliography

A wealth of genetic associations for cardiovascular and metabolic phenotypes in humans has been accumulating over the last decade, in particular a large number of loci derived from recent genome wide association studies (GWAS). True complex disease-associated loci often exert modest effects, so their delineation currently requires integration of diverse phenotypic data from large studies to ensure robust meta-analyses. We have designed a gene-centric 50 K single nucleotide polymorphism (SNP) array to assess potentially relevant loci across a range of cardiovascular, metabolic and inflammatory syndromes. The array utilizes a "cosmopolitan" tagging approach to capture the genetic diversity across approximately 2,000 loci in populations represented in the HapMap and SeattleSNPs projects. The array content is informed by GWAS of vascular and inflammatory disease, expression quantitative trait loci implicated in atherosclerosis, pathway based approaches and comprehensive literature searching. The custom flexibility of the array platform facilitated interrogation of loci at differing stringencies, according to a gene prioritization strategy that allows saturation of high priority loci with a greater density of markers than the existing GWAS tools, particularly in African HapMap samples. We also demonstrate that the IBC array can be used to complement GWAS, increasing coverage in high priority CVD-related loci across all major HapMap populations. DNA from over 200,000 extensively phenotyped individuals will be genotyped with this array with a significant portion of the generated data being released into the academic domain facilitating in silico replication attempts, analyses of rare variants and cross-cohort meta-analyses in diverse populations. These datasets will also facilitate more robust secondary analyses, such as explorations with alternative genetic models, epistasis and gene-environment interactions.

Alternative splicing of genes is an efficient means of generating variation in protein function. Several disease states have been associated with rare genetic variants that affect splicing patterns. Conversely, splicing efficiency of some genes is known to vary between individuals without apparent ill effects. What is not clear is whether commonly observed phenotypic variation in splicing patterns, and hence potential variation in protein function, is to a significant extent determined by naturally occurring DNA sequence variation and in particular by single nucleotide polymorphisms (SNPs). In this study, we surveyed the splicing patterns of 250 exons in 22 individuals who had been previously genotyped by the International HapMap Project. We identified 70 simple cassette exon alternative splicing events in our experimental system; for six of these, we detected consistent differences in splicing pattern between individuals, with a highly significant association between splice phenotype and neighbouring SNPs. Remarkably, for five out of six of these events, the strongest correlation was found with the SNP closest to the intron-exon boundary, although the distance between these SNPs and the intron-exon boundary ranged from 2 bp to greater than 1,000 bp. Two of these SNPs were further investigated using a minigene splicing system, and in each case the SNPs were found to exert cis-acting effects on exon splicing efficiency in vitro. The functional consequences of these SNPs could not be predicted using bioinformatic algorithms. Our findings suggest that phenotypic variation in splicing patterns is determined by the presence of SNPs within flanking introns or exons. Effects on splicing may represent an important mechanism by which SNPs influence gene function.

As an anti-inflammatory mediator IL10 is beneficial in certain contexts and deleterious in others. As increased production of IL10 favours protection against inflammatory disease, whereas low production promotes elimination of foreign pathogens by the host, we investigated the possible influence of balancing selection at this locus. We began by resequencing 48 European and 48 African chromosomes across 2.2 kb of the IL10 promoter region, and compared this with four neighbouring gene regions: MK2, IL19, IL20 and IL24. Analysis of nucleotide diversity showed a positive Tajima's D-test for IL10 in Europeans, of borderline statistical significance (1.89, P=0.05). Analysis of F(st) values showed significant population divergence at MK2, IL19, IL20 and IL24 (P<0.01) but not at IL10. Taken together, these findings are consistent with the hypothesis that balancing selection has played a role in the evolution of polymorphisms in the IL10 promoter region.

Dyslexia is one of the most prevalent childhood cognitive disorders, affecting approximately 5% of school-age children. We have recently identified a risk haplotype associated with dyslexia on chromosome 6p22.2 which spans the TTRAP gene and portions of THEM2 and KIAA0319. Here we show that in the presence of the risk haplotype, the expression of the KIAA0319 gene is reduced but the expression of the other two genes remains unaffected. Using in situ hybridization, we detect a very distinct expression pattern of the KIAA0319 gene in the developing cerebral neocortex of mouse and human fetuses. Moreover, interference with rat Kiaa0319 expression in utero leads to impaired neuronal migration in the developing cerebral neocortex. These data suggest a direct link between a specific genetic background and a biological mechanism leading to the development of dyslexia: the risk haplotype on chromosome 6p22.2 down-regulates the KIAA0319 gene which is required for neuronal migration during the formation of the cerebral neocortex.

There is presently much interest in utilizing patterns of linkage disequilibrium (LD) to further genetic association studies. This is particularly pertinent in the class III region of the human major histocompatibility complex (MHC), which has been extensively studied as a disease susceptibility locus in a number of ethnic groups. To date, however, few studies of LD in the MHC have considered non-Caucasian populations. With the advent of large-scale haplotyping of the human genome, the question of utilizing LD patterns across populations has come to the fore. We have previously used LD mapping to direct an MHC class III association study in a UK Caucasian population. As an extension of this, we sought to determine to what extent the pattern of LD observed in that study could be used to conduct a similar study in a West African Gambian population. We found that broad patterns of LD were similar in the two populations, resulting in similar candidate region delineations, but at a higher resolution, marker-specific patterns of LD and population-dependent allele frequencies confounded the choice of regional tagging SNPs. Our results have implications for the applicability of large-scale haplotype maps such as the HapMap to complex regions like the MHC.

Genetic variation at the human LTA locus, encoding lymphotoxin-alpha, is associated with susceptibility to myocardial infarction, asthma and other diseases. By detailed haplotypic analysis of the locus, we identified a single-nucleotide polymorphism (SNP) at LTA+80 as a main predictor of LTA protein production by human B cells. We found that activated B-cell factor-1 (ABF-1) binds to this site in vitro and suppresses reporter gene expression, but only in the presence of the LTA+80A allele. Using haplotype-specific chromatin immunoprecipitation, we confirmed that ABF-1 is preferentially recruited to the low-producer allele in vivo. These findings provide a molecular model of how LTA expression may be genetically regulated by allele-specific recruitment of the transcriptional repressor ABF-1.

In vivo characterization of regulatory polymorphisms is a key requirement for next-generation human genetic analysis. Here we describe haploChIP, a method that uses chromatin immunoprecipitation (ChIP) and mass spectrometry to identify differential protein-DNA binding in vivo associated with allelic variants of a gene. We demonstrate this approach with the imprinted gene SNRPN. HaploChIP showed close correlation between the level of bound phosphorylated RNA polymerase II at the SNRPN locus and allele-specific expression. Application of the approach to the TNF/LTA locus identified functionally important haplotypes that correlate with allele-specific transcription of LTA. The haploChIP method may be useful in high-throughput screening for common DNA polymorphisms that affect gene regulation in vivo.