SNP-based HLA allele tagging, imputation and association with antiepileptic drug-induced cutaneous reactions in Hong Kong Han Chinese.

Human leukocyte antigen (HLA) genes control the regulation of the human immune system and are involved in immune-related diseases. Population surveys on relationships between single nucleotide polymorphisms (SNP) and HLA alleles are essential to conduct genetic association between HLA variants and diseases. Samples were obtained from our in-house database for epilepsy genetics and pharmacogenetics research. Using 184 epilepsy patients with both genome-wide SNP array and HLA-A/B candidate gene sequencing data, we sought tagging SNPs that completely represent sixHLA risk alleles; in addition, a Hong Kong population-specific reference panel was constructed for SNP-based HLA imputation. The performance of our new panel was compared to a recent Han Chinese panel. Finally, genetic associations of HLA variants with mild skin rash were performed on the combined sample of 408 patients. Common SNPs rs2571375 and rs144295468 were found to successfully tag HLA risk alleles A*31:01 and B*13:01, respectively. HLA-B*15:02 can be predicted by rs144012689 with >95% sensitivity and specificity. The imputation reference panel for the Hong Kong population had comparable performance to the Han Chinese panel due to the large sample size for common HLA alleles, though it retained discordance for imputing rare alleles. No significant genetic associations were found between HLA genetic variants and mild skin rash induced by aromatic antiepileptic drugs. This study provides new information on the genetic structure of HLA regions in the Hong Kong population by identifying tagging SNPs and serving as a reference panel. Moreover, our comprehensive genetic analyses revealed no significant association between HLA alleles and mild skin rash in Hong Kong Han Chinese.

Reciprocal causation models of cognitive vs volumetric cerebral intermediate phenotypes for schizophrenia in a pan-European twin cohort.

In aetiologically complex illnesses such as schizophrenia, there is no direct link between genotype and phenotype. Intermediate phenotypes could help clarify the underlying biology and assist in the hunt for genetic vulnerability variants. We have previously shown that cognition shares substantial genetic variance with schizophrenia; however, it is unknown if this reflects pleiotropic effects, direct causality or some shared third factor that links both, for example, brain volume (BV) changes. We quantified the degree of net genetic overlap and tested the direction of causation between schizophrenia liability, brain structure and cognition in a pan-European schizophrenia twin cohort consisting of 1243 members from 626 pairs. Cognitive deficits lie upstream of the liability for schizophrenia with about a quarter of the variance in liability to schizophrenia explained by variation in cognitive function. BV changes lay downstream of schizophrenia liability, with 4% of BV variation explained directly by variation in liability. However, our power to determine the nature of the relationship between BV deviation and schizophrenia liability was more limited. Thus, while there was strong evidence that cognitive impairment is causal to schizophrenia liability, we are not in a position to make a similar statement about the relationship between liability and BV. This is the first study to demonstrate that schizophrenia liability is expressed partially through cognitive deficits. One prediction of the finding that BV changes lie downstream of the disease liability is that the risk loci that influence schizophrenia liability will thereafter influence BV and to a lesser extent. By way of contrast, cognitive function lies upstream of schizophrenia, thus the relevant loci will actually have a larger effect size on cognitive function than on schizophrenia. These are testable predictions.

Genome-wide association study in a Chinese population identifies a susceptibility locus for type 2 diabetes at 7q32 near PAX4.

AIMS/HYPOTHESIS: Most genetic variants identified for type 2 diabetes have been discovered in European populations. We performed genome-wide association studies (GWAS) in a Chinese population with the aim of identifying novel variants for type 2 diabetes in Asians. METHODS: We performed a meta-analysis of three GWAS comprising 684 patients with type 2 diabetes and 955 controls of Southern Han Chinese descent. We followed up the top signals in two independent Southern Han Chinese cohorts (totalling 10,383 cases and 6,974 controls), and performed in silico replication in multiple populations. RESULTS: We identified CDKN2A/B and four novel type 2 diabetes association signals with p < 1 × 10(-5) from the meta-analysis. Thirteen variants within these four loci were followed up in two independent Chinese cohorts, and rs10229583 at 7q32 was found to be associated with type 2 diabetes in a combined analysis of 11,067 cases and 7,929 controls (p meta = 2.6 × 10(-8); OR [95% CI] 1.18 [1.11, 1.25]). In silico replication revealed consistent associations across multiethnic groups, including five East Asian populations (p meta = 2.3 × 10(-10)) and a population of European descent (p = 8.6 × 10(-3)). The rs10229583 risk variant was associated with elevated fasting plasma glucose, impaired beta cell function in controls, and an earlier age at diagnosis for the cases. The novel variant lies within an islet-selective cluster of open regulatory elements. There was significant heterogeneity of effect between Han Chinese and individuals of European descent, Malaysians and Indians. CONCLUSIONS/INTERPRETATION: Our study identifies rs10229583 near PAX4 as a novel locus for type 2 diabetes in Chinese and other populations and provides new insights into the pathogenesis of type 2 diabetes.

Association of CD247 with systemic lupus erythematosus in Asian populations.

OBJECTIVE: Systemic lupus erythematosus (SLE) is a prototypic autoimmune disease with complex genetic inheritance. CD247 (CD3Z, TCRZ) plays a vital role in antigen recognition and signal transduction in antigen-specific immune responses, and is known to be involved in SLE pathogenesis. Weak disease association was reported for genetic variants in this gene in Caucasian studies for SLE, Crohn's disease and systemic sclerosis, but its role as a genetic risk factor was never firmly established. METHODS: In this study, using a collection of 612 SLE patients and 2193 controls of Chinese ethnicity living in Hong Kong in a genome-wide study, single nucleotide polymorphisms (SNPs) in and around CD247 were identified as being associated with SLE. The two most significant SNPs in this locus were selected for further replication using TaqMan genotyping assay in 3339 Asian patients from Hong Kong, Mainland China, and Thailand, as well as 4737 ethnically and geographically matched controls. RESULTS: The association of CD247 with SLE in Asian populations was confirmed (rs704853: odds ratio [OR] = 0. 81, p = 2.47 × 10(-7); rs858543: OR = 1.10, p = 0.0048). Patient-only analysis suggested that rs704853 is also linked to oral ulcers, hematologic disorders and anti-double-stranded DNA (dsDNA) antibody production. CONCLUSION: A significant association between variants in CD247 and SLE was demonstrated in Asian populations. Understanding the involvement of CD247 in SLE may shed new light on disease mechanisms and development of new treatment paradigms.

C-reactive protein as a predictor of hypertension in the Hong Kong Cardiovascular Risk Factor Prevalence Study (CRISPS) cohort.

Inflammation contributes to the development of hypertension. Whether C-reactive protein (CRP) has a causal role in hypertension remains unknown. We studied the relationship between circulating CRP levels and hypertension. The role of single-nucleotide polymorphisms (SNPs) in the CRP gene as determinants of its plasma levels and the propensity to develop hypertension was investigated. Plasma CRP and genotypes of nine SNPs were determined in 1925 unrelated subjects from the Hong Kong Cardiovascular Risk Factor Prevalence Study-2 (CRISPS-2) in 2000-2004. Among 1378 subjects normotensive in CRISPS-2, 1115 subjects had been followed up in CRISPS-3 after a median interval of 5.3 years, 236 of whom had developed hypertension. Plasma CRP was independently associated with the development of hypertension in CRISPS-3 (odds ratio per quartile=1.26, P=0.010). Six SNPs were associated with plasma CRP (all P<0.001). However, none of the SNPs was significantly associated with blood pressure, prevalent or incident hypertension, or change in blood pressure. In conclusion, plasma CRP predicts the development of hypertension. Genetic variants in the CRP gene are significantly associated with plasma CRP but not with hypertension. The future risk of hypertension is therefore more related to plasma CRP than SNPs in the CRP gene in this population.

Two missense variants in UHRF1BP1 are independently associated with systemic lupus erythematosus in Hong Kong Chinese.

UHRF1BP1 encodes a highly conserved protein with unknown function. Previously, a coding variant in this gene was found to be associated with systemic lupus erythematosus (SLE) in populations of European ancestry (rs11755393, R454Q, P=2.22 x 10⁻8, odds ratio=1.17). In this study, by a combination of genome-wide study and replication involving a total of 1230 patients and 3144 controls, we confirmed the association of this coding variant to SLE in Hong Kong Chinese. We also identified another coding variant in this gene that independently contributes to SLE susceptibility (rs13205210, M1098T, P=4.44 x 10⁻9, odds ratio=1.49). Cross-population confirmation establishes the involvement of this locus in SLE and indicates that distinct alleles are contributing to disease susceptibility.

Replication study of SNP associations for colorectal cancer in Hong Kong Chinese.

BACKGROUND: Recent genome-wide association studies of colorectal cancer (CRC) have identified common single-nucleotide polymorphisms (SNPs) mapping to 10 independent loci that confer modest increased risk. These studies have been conducted in European populations and it is unclear whether these observations generalise to populations with different ethnicities and rates of CRC. METHODS: An association study was performed on 892 CRC cases and 890 controls recruited from the Hong Kong Chinese population, genotyping 32 SNPs, which were either associated with CRC in previous studies or are in close proximity to previously reported risk SNPs. RESULTS: Twelve of the SNPs showed evidence of an association. The strongest associations were provided by rs10795668 on 10p14, rs4779584 on 15q14 and rs12953717 on 18q21.2. There was significant linear association between CRC risk and the number of independent risk variants possessed by an individual (P=2.29 × 10(-5)). CONCLUSION: These results indicate that some previously reported SNP associations also impact on CRC risk in the Chinese population. Possible reasons for failure of replication for some loci include inadequate study power, differences in allele frequency, linkage disequilibrium structure or effect size between populations. Our results suggest that many associations for CRC are likely to generalise across populations.

Fine mapping of the 9q31 Hirschsprung's disease locus.

Hirschsprung's disease (HSCR) is a congenital disorder characterised by the absence of ganglia along variable lengths of the intestine. The RET gene is the major HSCR gene. Reduced penetrance of RET mutations and phenotypic variability suggest the involvement of additional modifying genes in the disease. A RET-dependent modifier locus was mapped to 9q31 in families bearing no coding sequence (CDS) RET mutations. Yet, the 9q31 causative locus is to be identified. To fine-map the 9q31 region, we genotyped 301 tag-SNPs spanning 7 Mb on 137 HSCR Dutch trios. This revealed two HSCR-associated regions that were further investigated in 173 Chinese HSCR patients and 436 controls using the genotype data obtained from a genome-wide association study recently conducted. Within one of the two identified regions SVEP1 SNPs were found associated with Dutch HSCR patients in the absence of RET mutations. This ratifies the reported linkage to the 9q31 region in HSCR families with no RET CDS mutations. However, this finding could not be replicated. In Chinese, HSCR was found associated with IKBKAP. In contrast, this association was stronger in patients carrying RET CDS mutations with p = 5.10 x 10(-6) [OR = 3.32 (1.99, 5.59)] after replication. The HSCR-association found for IKBKAP in Chinese suggests population specificity and implies that RET mutation carriers may have an additional risk. Our finding is supported by the role of IKBKAP in the development of the nervous system.

Application of genome-wide SNP data for uncovering pairwise relationships and quantitative trait loci.

The genetic analysis of quantitative traits in humans is changing as a result of the availability of whole-genome SNP data. Heritability analysis can make use of actual genetic sharing between pairs of individuals estimated from the genotype data, rather than the expected genetic sharing implied by their family relationship. This could provide more accurate heritability estimates and help to overcome the equal environment assumption. Quantitative trait locus (QTL) linkage mapping can make use of local genetic sharing inferred from very dense local genotype data from pedigree members or individuals not previously known to be related. This approach may be particularly suited for detecting loci that contain rare variants with major effect on the phenotype. Finally, whole-genome SNP data can be used to measure the genetic similarity between individuals to provide matched sets for association studies, in order to avoid spurious association from population stratification.

Repeat instability in the 27-39 CAG range of the HD gene in the Venezuelan kindreds: Counseling implications.

The instability of the CAG repeat size of the HD gene when transmitted intergenerationally has critical implications for genetic counseling practices. In particular, CAG repeats between 27 and 35 have been the subject of debate based on small samples. To address this issue, we analyzed allelic instability in the Venezuelan HD kindreds, the largest and most informative families ascertained for HD. We identified 647 transmissions. Our results indicate that repeats in the 27-35 CAG range are highly stable. Out of 69 transmitted alleles in this range, none expand into any penetrant ranges. Contrastingly, 14% of alleles transmitted from the incompletely penetrant range (36-39 CAGs) expand into the completely penetrant range, characterized by alleles with 40 or more CAG repeats. At least 12 of the 534 transmissions from the completely penetrant range contract into the incompletely penetrant range of 36-39 CAG repeats. In these kindreds, none of the individuals with 27-39 CAGs were symptomatic, even though they ranged in age from 11 to 82 years. We expect these findings to be helpful in updating genetic counseling practices.

Genetic influences on the difference in variability of height, weight and body mass index between Caucasian and East Asian adolescent twins.

OBJECTIVE: Twin studies are useful for investigating the causes of trait variation between as well as within a population. The goals of the present study were two-fold: First, we aimed to compare the total phenotypic, genetic and environmental variances of height, weight and BMI between Caucasians and East Asians using twins. Secondly, we intended to estimate the extent to which genetic and environmental factors contribute to differences in variability of height, weight and BMI between Caucasians and East Asians. DESIGN: Height and weight data from 3735 Caucasian and 1584 East Asian twin pairs (age: 13-15 years) from Australia, China, Finland, Japan, the Netherlands, South Korea, Taiwan and the United States were used for analyses. Maximum likelihood twin correlations and variance components model-fitting analyses were conducted to fulfill the goals of the present study. RESULTS: The absolute genetic variances for height, weight and BMI were consistently greater in Caucasians than in East Asians with corresponding differences in total variances for all three body measures. In all 80 to 100% of the differences in total variances of height, weight and BMI between the two population groups were associated with genetic differences. CONCLUSION: Height, weight and BMI were more variable in Caucasian than in East Asian adolescents. Genetic variances for these three body measures were also larger in Caucasians than in East Asians. Variance components model-fitting analyses indicated that genetic factors contributed to the difference in variability of height, weight and BMI between the two population groups. Association studies for these body measures should take account of our findings of differences in genetic variances between the two population groups.

OpenADAM: an open source genome-wide association data management system for Affymetrix SNP arrays.

BACKGROUND: Large scale genome-wide association studies have become popular since the introduction of high throughput genotyping platforms. Efficient management of the vast array of data generated poses many challenges. DESCRIPTION: We have developed an open source web-based data management system for the large amount of genotype data generated from the Affymetrix GeneChip Mapping Array and Affymetrix Genome-Wide Human SNP Array platforms. The database supports genotype calling using DM, BRLMM, BRLMM-P or Birdseed algorithms provided by the Affymetrix Power Tools. The genotype and corresponding pedigree data are stored in a relational database for efficient downstream data manipulation and analysis, such as calculation of allele and genotype frequencies, sample identity checking, and export of genotype data in various file formats for analysis using commonly-available software. A novel method for genotyping error estimation is implemented using linkage disequilibrium information from the HapMap project. All functionalities are accessible via a web-based user interface. CONCLUSION: OpenADAM provides an open source database system for management of Affymetrix genome-wide association SNP data.

Association of F11 receptor gene polymorphisms with central obesity and blood pressure.

OBJECTIVES: F11 receptor, also known as junctional adhesion molecule-1, in the autonomic nervous system is implicated in the development of hypertension in spontaneous hypertensive rats. We investigated the association of single nucleotide polymorphisms (SNPs) in the F11 receptor gene (F11R) with hypertension and central obesity in Hong Kong Chinese. METHODS: Seven tagging SNPs were identified in the HapMap database. Genotyping was performed using Sequenom MassArray in 263 hypertensive subjects and 393 normotensive controls, of whom 263 matched the cases in age and sex. RESULTS: When subjects on anti-hypertensive medication were excluded, rs790056 and rs2774276 were associated with lower systolic blood pressure (TT:124.8 +/- 18.3 mmHg vs. TC + CC: 120.2 +/- 15.5 mmHg, P = 0.004 and CC: 124.7 +/- 18.5 mmHg vs. CG+GG: 120.5 +/- 15.1 mmHg, P = 0.007 respectively). Comparing 213 subjects with central obesity with 213 controls matched for sex and age, rs2481084 and rs3737787 were associated with lower odds of central obesity (odds ratio = 0.516, P = 0.002 and odds ratio = 0.540, P = 0.005 respectively). All these associations remained significant after correction for multiple testing. Analysis of statistically similar SNPs suggested that the causative variants for systolic blood pressure were located in F11R, whilst those for central obesity could be due to causative variants in the transcription factor 1 gene immediately upstream. CONCLUSIONS: F11 receptor plays a role in blood pressure regulation, not only in rats but also in man. The link between F11 receptor and central obesity merits further investigation.

Replication of twelve association studies for Huntington's disease residual age of onset in large Venezuelan kindreds.

BACKGROUND: The major determinant of age of onset in Huntington's disease is the length of the causative triplet CAG repeat. Significant variance remains, however, in residual age of onset even after repeat length is factored out. Many genetic polymorphisms have previously shown evidence of association with age of onset of Huntington's disease in several different populations. OBJECTIVE: To replicate these genetic association tests in 443 affected people from a large set of kindreds from Venezuela. METHODS: Previously tested polymorphisms were analysed in the HD gene itself (HD), the GluR6 kainate glutamate receptor (GRIK2), apolipoprotein E (APOE), the transcriptional coactivator CA150 (TCERG1), the ubiquitin carboxy-terminal hydrolase L1 (UCHL1), p53 (TP53), caspase-activated DNase (DFFB), and the NR2A and NR2B glutamate receptor subunits (GRIN2A, GRIN2B). RESULTS: The GRIN2A single-nucleotide polymorphism explains a small but considerable amount of additional variance in residual age of onset in our sample. The TCERG1 microsatellite shows a trend towards association but does not reach statistical significance, perhaps because of the uninformative nature of the polymorphism caused by extreme allele frequencies. We did not replicate the genetic association of any of the other genes. CONCLUSIONS: GRIN2A and TCERG1 may show true association with residual age of onset for Huntington's disease. The most surprising negative result is for the GRIK2 (TAA)(n) polymorphism, which has previously shown association with age of onset in four independent populations with Huntington's disease. The lack of association in the Venezuelan kindreds may be due to the extremely low frequency of the key (TAA)(16) allele in this population.

Introduction to the special issue: human linkage studies for behavioral traits.

In the post Genome era, the aim of behavior genetics has shifted from estimating the relative contributions of genes and environmental factors to (co-)variation in human complex traits, to localization of genes and identification of functional genetic variants. This special issue reflects this transition and presents fifteen papers that report on genome-wide linkage scans for complex traits in humans and on methodological tools and innovations. Six papers focus on cognition and report overlapping linkage peaks on chromosomes 6p and 14p. Papers on addictive behavior, i.e. smoking and alcohol dependence and its endophenotypes, find moderate LOD scores on chromosomes 6p, 5q, 4p and 7q, respectively. Three papers concentrate on emotionality, depression and loneliness and examine chromosomes 2q and 12q. The papers in this issue represent a summary of the first large scale linkage enterprises of human behavioral traits.

Genetic Power Calculator: design of linkage and association genetic mapping studies of complex traits.

SUMMARY: A website for performing power calculations for the design of linkage and association genetic mapping studies of complex traits. AVAILABILITY: The package is made available athttp://statgen.iop.kcl.ac.uk/gpc/.

GRR: graphical representation of relationship errors.

SUMMARY: A graphical tool for verifying assumed relationships between individuals in genetic studies is described. GRR can detect many common errors using genotypes from many markers. AVAILABILITY: GRR is available at http://bioinformatics.well.ox.ac.uk/GRR.

Optimal sibship selection for genotyping in quantitative trait locus linkage analysis.

In this paper we present a novel method for selecting optimally informative sibships of any size for quantitative trait locus (QTL) linkage analysis. The method allocates a quantitative index of potential informativeness to each sibship on the basis of observed trait scores and an assumed true QTL model. Any sample of phenotypically screened sibships can therefore be easily rank-ordered for selective genotyping. The quantitative index is the sibship's expected contribution to the non-centrality parameter. This expectation represents the weighted sum of chi(2) test statistics that would be obtained given the observed trait values over all possible sibship genotypic configurations; each configuration is weighted by the likelihood of it occurring given the assumed true genetic model. The properties of this procedure are explored in relation to the accuracy of the assumed true genetic model and sibship size. In comparison to previous methods of selecting phenotypically extreme sibships for genotyping, the proposed method is considerably more efficient and is robust with regard to the specification of the genetic model.

A powerful and rapid approach to human genome scanning using small quantities of genomic DNA.

Dense maps of short-tandem-repeat polymorphisms (STRPs) have allowed genome-wide searches for genes involved in a great variety of diseases with genetic influences, including common complex diseases. Generally for this purpose, marker sets with a 10 cM spacing are genotyped in hundreds of individuals. We have performed power simulations to estimate the maximum possible intermarker distance that still allows for sufficient power. In this paper we further report on modifications of previously published protocols, resulting in a powerful screening set containing 229 STRPs with an average spacing of 18.3 cM. A complete genome scan using our protocol requires only 80 multiplex PCR reactions which are all carried out using one set of conditions and which do not contain overlapping marker allele sizes. The multiplex PCR reactions are grouped by sets of chromosomes, which enables on-line statistical analysis of a set of chromosomes, as sets of chromosomes are being genotyped. A genome scan following this modified protocol can be performed using a maximum amount of 2.5 micrograms of genomic DNA per individual, isolated from either blood or from mouth swabs.