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.

High-resolution mapping of quantitative trait loci in outbred mice.

Screening the whole genome of a cross between two inbred animal strains has proved to be a powerful method for detecting genetic loci underlying quantitative behavioural traits, but the level of resolution offered by quantitative trait loci (QTL) mapping is still too coarse to permit molecular cloning of the genetic determinants. To achieve high-resolution mapping, we used an outbred stock of mice for which the entire genealogy is known. The heterogeneous stock (HS) was established 30 years ago from an eight-way cross of C57BL/6, BALB/c, RIII, AKR, DBA/2, I, A/J and C3H inbred mouse strains. At the time of the experiment reported here, the HS mice were at generation 58, theoretically offering at least a 30-fold increase in resolution for QTL mapping compared with a backcross or an F2 intercross. Using the HS mice we have mapped a QTL influencing a psychological trait in mice to a 0.8-cM interval on chromosome 1. This method allows simultaneous fine mapping of multiple QTLs, as shown by our report of a second QTL on chromosome 12. The high resolution possible with this approach makes QTLs accessible to positional cloning.

Linkage between sexual orientation and chromosome Xq28 in males but not in females.

We have extended our analysis of the role of the long arm of the X chromosome (Xq28) in sexual orientation by DNA linkage analyses of two newly ascertained series of families that contained either two gay brothers or two lesbian sisters as well as heterosexual siblings. Linkage between the Xq28 markers and sexual orientation was detected for the gay male families but not for the lesbian families or for families that failed to meet defined inclusion criteria for the study of sex-linked sexual orientation. Our results corroborate the previously reported linkage between Xq28 and male homosexuality in selected kinships and suggest that this region contains a locus that influences individual variations in sexual orientation in men but not in women.

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.

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.

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.

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.

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.

The impact of genotyping error on family-based analysis of quantitative traits.

Errors in genotyping can substantially influence the power to detect linkage using affected sib-pairs, but it is not clear what effect such errors have on quantitative trait analyses. Here we use Monte Carlo simulation to examine the influence of genotyping error on multipoint vs two-point analysis, variable map density, locus effect size and allele frequency in quantitative trait linkage and association studies of sib-pairs. The analyses are conducted using variance components methods. We contrast the effects of error on quantitative trait analyses with those on the affected sib-pair design. The results indicate that genotyping error influences linkage studies of affected sib pairs more severely than studies of quantitative traits in unselected sibs. In situations of modest effect size, 5% genotyping error eliminates all supporting evidence for linkage to a true susceptibility locus in affected pairs, but may only result in a loss of 15% of linkage information in random pairs. Multipoint analysis does not suffer substantially more than two-point analysis; for moderate error rates (< 5%), multipoint analysis with error is more powerful than two-point with no error. Map density does not appear to be an important factor for linkage analysis. QTL association analyses of common alleles are reasonably robust to genotyping error but power can be affected dramatically with rare alleles.

The familial aggregation of depressive symptoms, antisocial behavior, and alcohol abuse.

This study describes results from an ongoing family study of adolescent boys and their families designed to investigate potential risk factors for substance abuse. The adolescent treatment probands have severe drug and alcohol related problems and were recruited through a residential rehabilitation program. To date, the sample includes 251 individuals: 39 male probands and their families and 34 control families matched for age and geographic location (zip code). Probands and participating family members are given a structured interview which assesses alcohol and drug problems, and various psychiatric symptoms. The purpose of the present study was to examine the coaggregation of depressive symptoms, antisocial behavior, and alcohol misuse. Multivariate pedigree analyses were performed using a model that allowed for the estimation of vertical familial transmission, residual sibling resemblance, and assortative mating. Spouse correlations were estimated at .57, .21, and .31 for antisocial behavior, depressive symptoms, and alcohol abuse, respectively. Residual sibling environment (i.e., sibling resemblance unaccounted for by parent-offspring transmission) was not found for alcohol problem symptoms, but did contribute to resemblance for antisocial behavior and depressive symptoms. The proportion of variance accounted for by vertical familial transmission was estimated at approximately 30 to 40%. More important, correlations among the transmissible family factors for these psychiatric syndromes ranged from .58 to .73, suggesting substantial overlap among the underlying familial antecedents for these disorders.

Exploring the genetic etiology of low general cognitive ability from 14 to 36 months.

The genetic and environmental etiology of low general cognitive ability (g) during infancy and early childhood has not previously been investigated. The current study examined the genetic etiology of low cognitive ability at 14, 20, 24, and 36 months with twins from the MacArthur Longitudinal Twin Study. Low g groups were formed from the lowest 10th percentile at each age. Univariate probandwise concordance rates and DeFries-Fulker (J. C. DeFries & D. W. Fulker, 1985, 1988) multiple regression techniques suggest genetic etiology in low general cognitive ability groups. The stability of low general cognitive ability over time also appears to be primarily due to genetic factors. Although replication is necessary, these results suggest that the genetic etiology of low g during infancy and early childhood is at least as great as the heritability of g in the unselected population.

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.