Pedigree tests of transmission disequilibrium

Identification of common genetic variants that account for transcript isoform variation between human populations

In addition to the differences between populations in transcriptional and translational regulation of genes, alternative pre-mRNA splicing (AS) is also likely to play an important role in regulating gene expression and generating variation in mRNA and protein isoforms. Recently, the genetic contribution to transcript isoform variation has been reported in individuals of recent European descent. We report here results of an investigation of the differences in AS patterns between human populations. AS patterns in 176 HapMap lymphoblastoid cell lines derived from individuals of European and African ancestry were evaluated using the Affymetrix GeneChip Human Exon 1.0 ST Array. A variety of biological processes such as response to stimulus and transcription were found to be enriched among the differentially spliced genes. The differentially spliced genes also include some involved in human diseases that have different prevalence or susceptibility between populations. The genetic contribution to the population differences in transcript isoform variation was then evaluated by a genome-wide association using the HapMap genotypic data on single nucleotide polymorphisms (SNPs). The results suggest that local and distant genetic variants account for a substantial fraction of the observed transcript isoform variation between human populations. Our findings provide new insights into the complexity of the human genome as well as the health disparities between the two populations.

Mouse to human comparative genetics reveals a novel immunoglobulin E-controlling locus on Hsa8q12

Atopy is a predisposition to hyperproduction of immunoglobulin E (IgE) against common environmental allergens. It is often associated with development of allergic diseases such as asthma, rhinitis, and dermatitis. Production of IgE is influenced by genetic and environmental factors. In spite of progress in the study of heredity of atopy, the genetic mechanisms of IgE regulation have not yet been completely elucidated. The analysis of complex traits can benefit considerably from integration of human and mouse genetics. Previously, we mapped a mouse IgE-controlling locus Lmr9 on chromosome 4 to a segment of <9 Mb. In this study, we tested levels of total IgE and 25 specific IgEs against inhalant and food allergens in 67 Czech atopic families. In the position homologous to Lmr9 on chromosome 8q12 marked by D8S285, we demonstrated a novel human IgE-controlling locus exhibiting suggestive linkage to composite inhalant allergic sensitization (limit of detection, LOD = 2.11, P = 0.0009) and to nine specific IgEs, with maximum LOD (LOD = 2.42, P = 0.0004) to plantain. We also tested 16 markers at previously reported chromosomal regions of atopy. Linkage to plant allergens exceeding the LOD > 2.0 was detected at 5q33 (D5S1507, LOD = 2.11, P = 0.0009) and 13q14 (D13S165, LOD = 2.74, P = 0.0002). The significant association with plant allergens (quantitative and discrete traits) was found at 7p14 (D7S2250, corrected P = 0.026) and 12q13 (D12S1298, corrected P = 0.043). Thus, the finding of linkage on chromosome 8q12 shows precision and predictive power of mouse models in the investigation of complex traits in humans. Our results also confirm the role of loci at 5q33, 7p14, 12q14, and 13q13 in control of IgE.

Associations of ADH and ALDH2 gene variation with self report alcohol reactions, consumption and dependence: an integrated analysis

Alcohol dependence (AD) is a complex disorder with environmental and genetic origins. The role of two genetic variants in ALDH2 and ADH1B in AD risk has been extensively investigated. This study tested for associations between nine polymorphisms in ALDH2 and 41 in the seven ADH genes, and alcohol-related flushing, alcohol use and dependence symptom scores in 4597 Australian twins. The vast majority (4296) had consumed alcohol in the previous year, with 547 meeting DSM-IIIR criteria for AD. There were study-wide significant associations (P<2.3 x 10(-4)) between ADH1B-Arg48His (rs1229984) and flushing and consumption, but only nominally significant associations (P<0.01) with dependence. Individuals carrying the rs1229984 G-allele (48Arg) reported a lower prevalence of flushing after alcohol (P=8.2 x 10(-7)), consumed alcohol on more occasions (P=2.7 x 10(-6)), had a higher maximum number of alcoholic drinks in a single day (P=2.7 x 10(-6)) and a higher overall alcohol consumption (P=8.9 x 10(-8)) in the previous year than those with the less common A-allele (48His). After controlling for rs1229984, an independent association was observed between rs1042026 (ADH1B) and alcohol intake (P=4.7 x 10(-5)) and suggestive associations (P<0.001) between alcohol consumption phenotypes and rs1693482 (ADH1C), rs1230165 (ADH5) and rs3762894 (ADH4). ALDH2 variation was not associated with flushing or alcohol consumption, but was weakly associated with AD measures. These results bridge the gap between DNA sequence variation and alcohol-related behavior, confirming that the ADH1B-Arg48His polymorphism affects both alcohol-related flushing in Europeans and alcohol intake. The absence of study-wide significant effects on AD results from the low P-value required when testing multiple single nucleotide polymorphisms and phenotypes.

Genetic variants in the renin-angiotensin system genes are associated with cardiovascular-renal-related risk factors in Mexican Americans

The aim of this study is to examine whether the ACE-I/D, AGT-M235T, and AT1R-A1166C polymorphisms of the renin-angiotensin system (RAS) genes are associated with cardiovascular and renal-related risk factors in Mexican Americans. Study participants (N = 848) were genotyped by Taqman assays. Association analyses were performed by measured genotype approach. Of the phenotypes examined, the ACE-I/D, AGT-M235T, and AT1R-A1166C polymorphisms exhibited significant association with systolic blood pressure, glomerular filtration rate and body mass index, respectively. The data suggest that the polymorphisms examined in the RAS may modulate the risk factors associated with cardiovascular-renal disease.

Calculation of IBD probabilities with dense SNP or sequence data

The probabilities that two individuals share 0, 1, or 2 alleles identical by descent (IBD) at a given genotyped marker locus are quantities of fundamental importance for disease gene and quantitative trait mapping and in family-based tests of association. Until recently, genotyped markers were sufficiently sparse that founder haplotypes could be modelled as having been drawn from a population in linkage equilibrium for the purpose of estimating IBD probabilities. However, with the advent of high-throughput single nucleotide polymorphism genotyping assays, this is no longer a reasonable assumption. Indeed, the imminent arrival of individual sequencing will enable high-density single nucleotide polymorphism genotyping on a scale for which current algorithms are not equipped. In this paper, we present a simple new model in which founder haplotypes are modelled as a Markov chain. Another important innovation is that genotyping errors are explicitly incorporated into the model. We compare results obtained using the new model to those obtained using the popular genetic linkage analysis package Merlin, with and without using the cluster model of linkage disequilibrium that is incorporated into that program. We find that the new model results in accuracy approaching that of Merlin with haplotype blocks, but achieves this with orders of magnitude faster run times. Moreover, the new algorithm scales linearly with number of markers, irrespective of density, whereas Merlin scales supralinearly. We also confirm a previous finding that ignoring linkage disequilibrium in founder haplotypes can cause errors in the calculation of IBD probabilities.

Personality endophenotypes for bipolar affective disorder: a family-based genetic association analysis

Genetic analyses of complex conditions such as bipolar disorder (BD) may be facilitated by the use of intermediate phenotypes. Various personality traits are overrepresented in people with BD and their unaffected relatives, and may constitute genetically transmitted risk factors or endophenotypes of the illness. In this study, we administered a battery of seven different personality questionnaires comprising 19 subscales to 31 Caucasian BD families (n = 241). Ten of these personality traits showed significant evidence of heritability and were therefore selected as candidate endophenotypes. In addition, a principal components analysis produced two heritable components (negative affect and appetitive drive), which accounted for a considerable proportion of the variance (29% + 12%) and were also used in the analysis. A family-based quantitative association study was carried out using the orthogonal model from the quantitative transmission disequilibrium tests (QTDT) program. Monte Carlo permutations (M = 500), which allow for non-normal data and provide a global P value, corrected for multiple testing, were used to calculate empirical P values for the within-family component of association. The 3' untranslated region repeat polymorphism of the dopamine transporter gene (SLC6A3) was associated with self-directedness (P < 0.0001) and negative affect (P = 0.010). The short allele of the serotonin transporter gene (SLC6A4) promoter polymorphism showed a trend toward association with higher harm avoidance (P = 0.016) and negative affect (P = 0.028). The catechol-o-methyltransferase val158met polymorphism was weakly associated with the personality traits, 'Spirituality' (P = 0.040) and irritable temperament (P = 0.022). Furthermore, the met allele of the brain-derived neurotrophic factor val66met polymorphism was associated with higher hyperthymic temperament scores. We raise the possibility that the 10R allele of the SLC6A3 repeat polymorphism and the short allele of the SLC6A4 promoter variant constitute risk factors for irritable-aggressive and anxious-dysthymic subtypes of BD, respectively.

Effects of GABRA2 variation on physiological, psychomotor and subjective responses in the alcohol challenge twin study

Multiple reports have identified variation in the GABRA2 gene as contributing to the genetic susceptibility to alcohol dependence. However, both the mechanism behind this association, and the range of alcohol-related phenotypes affected by variation in this gene, are currently undefined. Other data suggest that the risk of alcohol dependence is increased by relative insensitivity to alcohol's intoxicating effects. We have therefore tested whether GABRA2 variation is associated with variation in the subjective and objective effects of a standard dose of alcohol in humans. Data on responses to alcohol from the Alcohol Challenge Twin Study (Martin et al., 1985) have been tested against allelic and haplotype information obtained by typing 41 single-nucleotide polymorphisms in or close to the GABRA2 gene. Nominally significant allelic associations (p<.05, without correction for multiple testing) were found for body sway, motor coordination, pursuit rotor and arithmetical computation tasks, and for the personality dimension of Neuroticism. Because of the large number of phenotypes tested, these possibly significant findings will need to be confirmed in further studies.

Ovarian morphology is a marker of heritable biochemical traits in sisters with polycystic ovaries

CONTEXT

Polycystic ovary syndrome (PCOS) is a common endocrinopathy of uncertain etiology but with strong evidence for a genetic contribution.

OBJECTIVE

The objective of the study was to test the hypothesis that the typical polycystic ovarian morphology is a marker of inherited biochemical features in families of women with PCOS.

DESIGN

A study of probands with PCOS and their sisters.

PATIENTS

Patients included 125 probands and 214 sisters. All probands had PCOS, defined by symptoms of anovulation and/or hyperandrogenism with polycystic ovaries on ultrasound. Affected sisters were defined by polycystic ovaries, regardless of symptoms, and unaffected sisters defined by normal ovarian morphology.

SETTING

This was a clinic-based study.

MAIN OUTCOME MEASURES

Clinical, endocrine, and metabolic features in the various groups were compared, and estimates of broad-sense heritability were obtained using the quantitative transmission disequilibrium test.

RESULTS

Although affected sisters had fewer symptoms than probands (30% had no symptoms of PCOS), serum testosterone, androstenedione, LH, and fasting insulin and insulin sensitivity were similar in the two groups with polycystic ovaries but significantly different from those in unaffected sisters or controls. We observed moderate to high heritabilities for all traits studied in affected sister pairs, whereas heritabilities calculated from discordant siblings were substantially lower.

CONCLUSIONS

These data provide further evidence for a genetic basis of PCOS. The high heritability of biochemical features in probands and affected sisters, despite wide variation in symptoms, shows that not only are these biochemical traits strongly influenced by genetic factors but also, importantly, that polycystic ovarian morphology is an index of inherited traits in families with PCOS.

The role of GABRA2 in alcohol dependence, smoking, and illicit drug use in an Australian population sample

BACKGROUND

Multiple studies have shown that genetic variation in the alpha-2 subunit of the gamma-aminobutyric acid type A (GABA(A)) receptor (GABRA2) is associated with risk for alcohol dependence. Recent reports have suggested that GABRA2 may exert its influence on dependence through factors such as sensitivity to alcohol's intoxicating effects and that GABRA2 may also contribute to a common underlying genetic vulnerability to both alcohol and drug dependence. The present study tested for association between GABRA2 and alcohol dependence, smoking, and illicit drug use within the Australian population.

METHODS

We genotyped 11 single nucleotide polymorphisms (SNPs) within or flanking GABRA2 in 4597 subjects (34.6% males) from 2618 families comprising 814 monozygotic pairs, 1177 dizygotic pairs, and 627 twins whose co-twin did not participate. Family-based association tests were conducted for binary and quantitative measures of alcohol dependence, smoking, and cannabis and other illicit drug use.

RESULTS

We observed evidence of association (p < 0.05) between multiple GABRA2 SNPs and quantitative measures of alcohol dependence, including symptom scores and principal component factor scores from the 9 criteria for DSM-IV alcohol dependence, in the opposite direction to that previously reported. In contrast, GABRA2 was not associated overall with dichotomous measure of alcohol dependence nor with smoking, cannabis, or illicit drug use.

CONCLUSIONS

The GABRA2 allelic associations found in clinical case-control studies have detectable but minor effects on DSM-defined alcohol dependence in the general community. Systematic comparisons of allelic effects on alcohol dependence in clinical cases and in the general community are required.

Polymorphisms of the tumor necrosis factor-alpha receptor 2 gene are associated with obesity phenotypes among 405 Caucasian nuclear families

The plasma level of the tumor necrosis factor-alpha receptor 2 (TNFR2) is associated with obesity phenotypes. However, the genetic polymorphisms for such an association have rarely been explored and are generally unknown. In this study, by employing a large sample of 1,873 subjects from 405 Caucasian nuclear families, we explored the association of 12 SNPs of the TNFR2 gene and obesity-related phenotypes, including body mass index (BMI), fat mass, and percentage fat mass (PFM). The within-family quantitative transmission disequilibrium test, which is robust to sample stratification, was implemented to evaluate the association of TNFR2 gene with obesity phenotypes. Evidence of association was obtained at SNP9 (rs5746059) with fat mass (P = 0.0002), BMI (P = 0.002), and PFM (P = 0.0006). The contribution of this polymorphism to the variation of fat mass and PFM was 6.24 and 7.82%, respectively. Individuals carrying allele A at the SNP9 site had a 4.6% higher fat mass and a 2.5% increased PFM compared to noncarriers. The results remained significant even after correction for multiple testing. Evidence of association between the TNFR2 gene and obesity phenotypes are also found in 700 independent Chinese Han and 1,000 random Caucasians samples. The results suggest that the TNFR2 gene polymorphisms contribute to the variation of obesity phenotypes.

Association of genetic variation in ENPP1 with obesity-related phenotypes

Ectonucleotide pyrophosphatase phosphodiesterase (ENPP1) is a positional candidate gene at chromosome 6q23 where we previously detected strong linkage with fasting-specific plasma insulin and obesity in Mexican Americans from the San Antonio Family Diabetes Study (SAFDS). We genotyped 106 single-nucleotide polymorphisms (SNPs) within ENPP1 in all 439 subjects from the linkage study, and measured association with obesity and metabolic syndrome (MS)-related traits. Of 72 polymorphic SNPs, 24 were associated, using an additive model, with at least one of eight key metabolic traits. Three traits were associated with at least four SNPs. They were high-density lipoprotein cholesterol (HDL-C), leptin, and fasting plasma glucose (FPG). HDL-C was associated with seven SNPs, of which the two most significant P values were 0.0068 and 0.0096. All SNPs and SNP combinations were analyzed for functional contribution to the traits using the Bayesian quantitative-trait nucleotide (BQTN) approach. With this SNP-prioritization analysis, HDL-C was the most strongly associated trait in a four-SNP model (P=0.00008). After accounting for multiple testing, we conclude that ENPP1 is not a major contributor to our previous linkage peak with MS-related traits in Mexican Americans. However, these results indicate that ENPP1 is a genetic determinant of these traits in this population, and are consistent with multiple positive association findings in independent studies in diverse human populations.

Endothelial nitric oxide synthase (eNOS) gene polymorphisms and their association with type 2 diabetes-related traits in Mexican Americans

Genetic variants of the endothelial nitric oxide synthase (eNOS) gene such as T-786C, Glu298Asp and 27bp-VNTR have been examined for their association with type 2 diabetes (T2DM)-related traits in various populations but not in Mexican Americans. However, the results from such studies have been controversial. This study investigated whether these three polymorphisms are associated with T2DM and its related traits in Mexican Americans, a population at high risk for T2DM and its complications. The study participants (n=670; 39 families) were genotyped for the three polymorphisms using polymerase chain reaction followed by restriction fragment length polymorphism assay. Association analyses between these polymorphisms and T2DM and its related phenotypes were carried out using a measured genotype approach as implemented in the computer program SOLAR. Of the variants examined, only the 27bp-VNTR variant exhibited significant association with high-density lipoprotein cholesterol (HDL-C) (p=0.04) and diastolic blood pressure (DBP) levels (p=0.02) after accounting for trait-specific covariates. The carriers of the rare allele (27bp-VNTR-4a) were found to have decreased HDL-C and increased DBP levels. In conclusion, of the genetic polymorphisms examined at the eNOS locus, only 27bp-VNTR appears to be a minor contributor to the variation in T2DM-related traits in Mexican Americans.

A regulatory SNP of the BICD1 gene contributes to telomere length variation in humans

Telomeres are repetitive sequences of variable length at the ends of chromosomes involved in maintaining their integrity. Telomere dysfunction is associated with increased risk of cancer and other age-related diseases. Telomere length is an important determinant of telomere function and has a strong genetic basis. We previously carried out a genome-wide linkage analysis of mean leukocyte telomere length, and identified a 12 cm quantitative-trait locus affecting telomere length on human chromosome 12. In the present study we confirmed linkage to this locus in an extended sample (380 families, 520 sib-pairs, maximum LOD score 4.3). Fine-mapping identified a 51 kb region of association within intron 1 of the Bicaudal-D homolog 1 (BICD1, MIM 602204) gene. The strongest association (P = 1.9 x 10(-5)) was with SNP rs2630578 where the minor allele C (frequency 0.21) was associated with telomeres that were shorter by 604 (+/-204) base pairs, equivalent to approximately 15-20 years of age-related attrition in telomere length. Subjects carrying the C allele for rs2630778 had 44% lower BICD1 mRNA levels in their leukocytes compared with GG homozygotes (P = 0.004). BICD1 is involved in Golgi-to-endoplasmic reticulum vacuolar transport. Previous studies have implicated vacuolar genes in telomere length homeostasis in yeast. Our study indicates that BICD1 plays a similar role in humans.

Association genetics in Pinus taeda L. II. Carbon isotope discrimination

Dissection of complex traits that influence fitness is not only a central topic in evolutionary research but can also assist breeding practices for economically important plant species, such as loblolly pine (Pinus taeda L). In this study, 46 single nucleotide polymorphisms (SNPs) from 41 disease and abiotic stress-inducible genes were tested for their genetic association with carbon isotope discrimination (CID), a time-integrated trait measure of stomatal conductance. A family-based approach to detect genotype/phenotype genetic association was developed for the first time in plants by applying the quantitative transmission disequilibrium test on an association population of 961 clones from 61 families (adopted from previous breeding programs) evaluated for phenotypic expression of CID at two sites. Two particularly promising candidates for their genetic effects on CID are: dhn-1, involved in stabilization of cell structures, and lp5-like, a glycine rich protein putatively related to cell wall reinforcement proteins, both of which were shown in previous studies to be water-deficit inducible. Moreover, association in lp5-like involves a nonsynonymous mutation in linkage disequilibrium with two other nonsynonymous polymorphisms that could, by acting together, enhance overall phenotypic effects. This study highlights the complexity of dissecting CID traits and provides insights for designing second-generation association studies based on candidate gene approaches in forest trees.

A new transmission test for affected sib-pair families

Family-based association approaches such as the transmission-disequilibrium test (TDT) are used extensively in the study of genetic traits because they are generally robust to the presence of population structure. However, these approaches necessarily involve recruitment of families, which is more costly and time-consuming than sampling unrelated individuals in the population-based approaches. Therefore, a family-based approach, which has high power, would be appealing because of the gain in time and cost due to the reduced sample size that is required to attain adequate power. Here we introduce a new family-based transmission test using the joint transmission status from affected sib pairs. We show that by including the transmission status of both siblings, our method gives higher power than the TDT design, while maintaining the correct type I error rate. We use the simulated data from affected sib-pair families with rheumatoid arthritis provided by Genetic Analysis Workshop 15 to illustrate our approach.

Human-mouse quantitative trait locus concordance and the dissection of a human neuroticism locus

BACKGROUND

Exploiting synteny between mouse and human disease loci has been proposed as a cost-effective method for the identification of human susceptibility genes. Here we explore its utility in an analysis of a human personality trait, neuroticism, which can be modeled in mice by tests of emotionality. We investigated a mouse emotionality locus on chromosome 1 that contains no annotated genes but abuts four regulators of G protein signaling, one of which (rgs2) has been previously identified as a quantitative trait gene for emotionality. This locus is syntenic with a human region that has been consistently implicated in the genetic aetiology of neuroticism.

METHODS

The functional candidacy of 29 murine sequence variants was tested by a combination of gel shift and transient transfection assays. Murine sequences that contained functional variants and exhibited significant cross-species conservation were prioritized for investigation in humans. Genetic association with neuroticism was tested in 1869 high and 2032 low unrelated individuals scored for neuroticism, selected from the extremes of 88,141 people from southwest England.

RESULTS

Fifteen sequence variants contributed to variation in the expression of rgs18, the gene lying at the edge of the quantitative trait loci (QTL) interval. There was no evidence of association between neuroticism and single nucleotide polymorphisms (SNPs) lying in the human regions homologous to those of mouse functional variants. One SNP, rs6428058, in a region of sequence conservation 644 kb upstream of RGS18, showed significant association (p = .000631).

CONCLUSIONS

It is unlikely that a single variant is responsible for the mouse emotionality locus on chromosome 1. This level of underlying genetic complexity means that although cross-species QTL concordance may be invaluable for the identification of human disease loci, it is unlikely to be as informative in the identification of human disease-causing variants.

Genetic architecture of transcript-level variation in humans

We report here the results of testing the pairwise association of 12,747 transcriptional gene-expression values with more than two million single-nucleotide polymorphisms (SNPs) in samples of European (CEPH from Utah; CEU) and African (Yoruba from Ibadan; YRI) ancestry. We found 4,677 and 5,125 significant associations between expression quantitative nucleotides (eQTNs) and transcript clusters in the CEU and the YRI samples, respectively. The physical distance between an eQTN and its associated transcript cluster was referred to as the intrapair distance. An association with 4 Mb or less intrapair distance was defined as local; otherwise, it was defined as distant. The enrichment analysis of functional categories shows that genes harboring the local eQTNs are enriched in the categories related to nucleosome and chromatin assembly; the genes harboring the distant eQTNs are enriched in the categories related to transmembrane signal transduction, suggesting that these biological pathways are likely to play a significant role in regulation of gene expression. We highlight in the EPHX1 gene a deleterious nonsynonymous SNP that is distantly associated with gene expression of ORMDL3, a susceptibility gene for asthma.

A general test for gene-environment interaction in sib pair-based association analysis of quantitative traits

Several association studies support the hypothesis that genetic variants can modify the influence of environmental factors on behavioral outcomes, i.e., G × E interaction. The case-control design used in these studies is powerful, but population stratification with respect to allele frequencies can give rise to false positive or false negative associations. Stratification with respect to the environmental factors can lead to false positives or false negatives with respect to environmental main effects and G × E interaction effects as well. Here we present a model based on Fulker et al. (1999) and Purcell (2002) for the study of G × E interaction in family-based association designs, in which the effects of stratification can be controlled. Simulations illustrate the power to detect genetic and environmental main effects, and G × E interaction effects for the sib pair design. The power to detect interaction was studied in eight different situations, both with and without the presence of population stratification, and for categorical and continuous environmental factors. Results show that the power to detect genetic and environmental main effects, and G × E interaction effects, depends on the allele frequencies and the distribution of the environmental moderator. Admixture effects of realistic effect size lead only to very small stratification effects in the G × E component, so impractically large numbers of sib pairs are required to detect such stratification.

Review and evaluation of methods correcting for population stratification with a focus on underlying statistical principles

When two or more populations have been separated by geographic or cultural boundaries for many generations, drift, spontaneous mutations, differential selection pressures and other factors may lead to allele frequency differences among populations. If these 'parental' populations subsequently come together and begin inter-mating, disequilibrium among linked markers may span a greater genetic distance than it typically does among populations under panmixia [see glossary]. This extended disequilibrium can make association studies highly effective and more economical than disequilibrium mapping in panmictic populations since less marker loci are needed to detect regions of the genome that harbor phenotype-influencing loci. However, under some circumstances, this process of intermating (as well as other processes) can produce disequilibrium between pairs of unlinked loci and thus create the possibility of confounding or spurious associations due to this population stratification. Accordingly, researchers are advised to employ valid statistical tests for linkage disequilibrium mapping allowing conduct of genetic association studies that control for such confounding. Many recent papers have addressed this need. We provide a comprehensive review of advances made in recent years in correcting for population stratification and then evaluate and synthesize these methods based on statistical principles such as (1) randomization, (2) conditioning on sufficient statistics, and (3) identifying whether the method is based on testing the genotype-phenotype covariance (conditional upon familial information) and/or testing departures of the marginal distribution from the expected genotypic frequencies.

The application of the entropy-based statistic for genomic association study of QTL

An entropy-based statistic T(PE) has been proposed for genomic association study for disease-susceptibility locus. The statistic T(PE) may be directly adopted and/or extended to quantitative-trait locus (QTL) mapping for quantitative traits. In this article, the statistic T(PE) was extended and applied to quantitative trait for association analysis of QTL by means of selective genotyping. The statistical properties (the type I error rate and the power) were examined under a range of parameters and population-sampling strategies (e.g., various genetic models, various heritabilities, and various sample-selection threshold values) by simulation studies. The results indicated that the statistic T(PE) is robust and powerful for genomic association study of QTL. A simulation study based on the haplotype frequencies of 10 single nucleotide polymorphisms (SNPs) of angiotensin-I converting enzyme genes was conducted to evaluate the performance of the statistic T(PE) for genetic association study.

Evaluation of genetic variation contributing to differences in gene expression between populations

Gene expression is a complex quantitative trait partially regulated by genetic variation in DNA sequence. Population differences in gene expression could contribute to some of the observed differences in susceptibility to common diseases and response to drug treatments. We characterized gene expression in the full set of HapMap lymphoblastoid cell lines derived from individuals of European and African ancestry for 9156 transcript clusters (gene-level) evaluated with the Affymetrix GeneChip Human Exon 1.0 ST Array. Gene expression was found to differ significantly between these samples for 383 transcript clusters. Biological processes including ribosome biogenesis and antimicrobial humoral response were found to be enriched in these differential genes, suggesting their possible roles in contributing to the population differences at a higher level than that of mRNA expression and in response to environmental information. Genome-wide association studies for local or distant genetic variants that correlate with the differentially expressed genes enabled identification of significant associations with one or more single-nucleotide polymorphisms (SNPs), consistent with the hypothesis that genetic factors and not simply population identity or other characteristics (age of cell lines, length of culture, etc.) contribute to differences in gene expression in these samples. Our results provide a comprehensive view of the genes differentially expressed between populations and the enriched biological processes involved in these genes. We also provide an evaluation of the contributions of genetic variation and nongenetic factors to the population differences in gene expression.

Power calculations using exact data simulation: a useful tool for genetic study designs

Statistical power calculations constitute an essential first step in the planning of scientific studies. If sufficient summary statistics are available, power calculations are in principle straightforward and computationally light. In designs, which comprise distinct groups (e.g., MZ & DZ twins), sufficient statistics can be calculated within each group, and analyzed in a multi-group model. However, when the number of possible groups is prohibitively large (say, in the hundreds), power calculations on the basis of the summary statistics become impractical. In that case, researchers may resort to Monte Carlo based power studies, which involve the simulation of hundreds or thousands of replicate samples for each specified set of population parameters. Here we present exact data simulation as a third method of power calculation. Exact data simulation involves a transformation of raw data so that the data fit the hypothesized model exactly. As in power calculation with summary statistics, exact data simulation is computationally light, while the number of groups in the analysis has little bearing on the practicality of the method. The method is applied to three genetic designs for illustrative purposes.

Combined linkage and association mapping of quantitative trait Loci with missing completely at random genotype data

In genetics study, the genotypes or phenotypes can be missing due to various reasons. In this paper, the impact of missing genotypes is investigated for high resolution combined linkage and association mapping of quantitative trait loci (QTL). We assume that the genotype data are missing completely at random (MCAR). Two regression models, "genotype effect model" and "additive effect model", are proposed to model the association between the markers and the trait locus. If the marker genotype is not missing, the model is exactly the same as those of our previous study, i.e., the number of genotype or allele is used as weight to model the effect of the genotype or allele in single marker case. If the marker genotype is missing, the expected number of genotype or allele is used as weight to model the effect of the genotype or allele. By analytical formulae, we show that the "genotype effect model" can be used to model the additive and dominance effects simultaneously, and the "additive effect model" can only be used to model the additive effect. Based on the two models, F-test statistics are proposed to test association between the QTL and markers. The non-centrality parameter approximations of F-test statistics are derived to calculate power and to compare power, which show that the power of the F-tests is reduced due to the missingness. By simulation study, we show that the two models have reasonable type I error rates for a dataset of moderate sample size. However, the type I error rates can be very slightly inflated if all individuals with missing genotypes are removed from analysis. Hence, the proposed method can help to get correct type I error rates although it does not improve power. As a practical example, the method is applied to analyze the angiotensin-1 converting enzyme (ACE) data.

Power of transmission/disequilibrium tests in admixed populations

The power of transmission/disequilibrium tests (TDTs) for detecting disease susceptibility loci is expected to be influenced by population admixture through its impact on the degree of linkage disequilibrium (LD) between the genetic marker and the DSL. However, few studies have been done to systematically examine this behavior of the TDTs in admixed populations. In the present study, extensive computer simulations were conducted to explore how population admixture affects the power of TDTs. It was found that (1) in newly admixed populations, the LD due to admixture makes no contribution to the power of TDTs, and it is the averaged background LD in the parental populations that determines the power of TDTs; but (2) after random mating between the admixed populations, the LD due to admixture becomes effective in increasing or decreasing the power of the tests, and (3) incomplete random mating can prolong the time for the LD due to admixture to become effective. This study clarifies the potential influence of population admixture on the performance of TDTs.

Haplotype-association analysis

Association methods based on linkage disequilibrium (LD) offer a promising approach for detecting genetic variations that are responsible for complex human diseases. Although methods based on individual single nucleotide polymorphisms (SNPs) may lead to significant findings, methods based on haplotypes comprising multiple SNPs on the same inherited chromosome may provide additional power for mapping disease genes and also provide insight on factors influencing the dependency among genetic markers. Such insights may provide information essential for understanding human evolution and also for identifying cis-interactions between two or more causal variants. Because obtaining haplotype information directly from experiments can be cost prohibitive in most studies, especially in large scale studies, haplotype analysis presents many unique challenges. In this chapter, we focus on two main issues: haplotype inference and haplotype-association analysis. We first provide a detailed review of methods for haplotype inference using unrelated individuals as well as related individuals from pedigrees. We then cover a number of statistical methods that employ haplotype information in association analysis. In addition, we discuss the advantages and limitations of different methods.

LRRTM1 on chromosome 2p12 is a maternally suppressed gene that is associated paternally with handedness and schizophrenia

Left-right asymmetrical brain function underlies much of human cognition, behavior and emotion. Abnormalities of cerebral asymmetry are associated with schizophrenia and other neuropsychiatric disorders. The molecular, developmental and evolutionary origins of human brain asymmetry are unknown. We found significant association of a haplotype upstream of the gene LRRTM1 (Leucine-rich repeat transmembrane neuronal 1) with a quantitative measure of human handedness in a set of dyslexic siblings, when the haplotype was inherited paternally (P=0.00002). While we were unable to find this effect in an epidemiological set of twin-based sibships, we did find that the same haplotype is overtransmitted paternally to individuals with schizophrenia/schizoaffective disorder in a study of 1002 affected families (P=0.0014). We then found direct confirmatory evidence that LRRTM1 is an imprinted gene in humans that shows a variable pattern of maternal downregulation. We also showed that LRRTM1 is expressed during the development of specific forebrain structures, and thus could influence neuronal differentiation and connectivity. This is the first potential genetic influence on human handedness to be identified, and the first putative genetic effect on variability in human brain asymmetry. LRRTM1 is a candidate gene for involvement in several common neurodevelopmental disorders, and may have played a role in human cognitive and behavioral evolution.

An entropy-based measure for QTL mapping using extreme samples of population

Quantitative trait locus (QTL) mapping can be accomplished through the method of selective genotyping, which is based on the differences of frequencies between an upper sample and a lower sample in population. However, amplifying the differences in marker allele frequencies in extreme samples may increase the probability for QTL mapping. Shannon entropy, which is a nonlinear function of allele frequencies, can be used to amplify the differences in marker allele frequencies. In this paper, we present a novel measure for linkage disequilibrium (LD) between a marker and single QTL, that is based on the comparison of the entropy and conditional entropy in a marker in extreme samples of population. This measure of LD between the marker and the trait locus can be used when the marker allele frequencies are known in the extreme samples of a population. We investigate the mapping performance in both analytic and simulation scenarios of a single QTL linked to a single marker. Our results show that the measure has very reasonable performance. In addition, a simulation study is performed on the basis of the haplotype frequencies of 10 SNPs of angiotensin-I converting enzyme (ACE) genes.

A chromosome 11q quantitative-trait locus influences change of blood-pressure measurements over time in Mexican Americans of the San Antonio Family Heart Study

Although previous genome scans have searched for quantitative-trait loci (QTLs) influencing variation in blood pressure (BP), few have investigated the rate of change in BP over time as a phenotype. Here, we compare results from genomewide scans to localize QTLs for systolic, diastolic, and mean arterial BPs (SBP, DBP, and MBP, respectively) and for rates of change in systolic, diastolic, and mean arterial BPs (rSBP, rDBP, and rMBP, respectively), with use of the longitudinal data collected about Mexican Americans of the San Antonio Family Heart Study (SAFHS). Significant evidence of linkage was found for rSBP (LOD 4.15) and for rMBP (LOD 3.94) near marker D11S4464 located on chromosome 11q24.1. This same chromosome 11q region also shows suggestive linkage to SBP (LOD 2.23) and MBP (LOD 2.37) measurements collected during the second clinic visit. Suggestive evidence of linkage to chromosome 5 was also found for rMBP, to chromosome 16 for rSBP, and to chromosomes 1, 5, 6, 7, and 21 for the single-time-point BP traits collected at the first two SAFHS clinic visits. We also present results from fine mapping the chromosome 11 QTL with use of SNP-association analysis within candidate genes identified from a bioinformatic search of the region and from whole-genome transcriptional expression data collected from 1,240 SAFHS participants. Our results show that the use of longitudinal BP data to calculate the rate of change in BP over time provides more information than do the single-time measurements, since they reveal physiological trends in the subjects that a single-time measurement could never capture. Further investigation of this region is necessary for the identification of the genetic variation responsible for QTLs influencing the rate of change in BP.

Identification of promoter variants in baboon endothelial lipase that regulate high-density lipoprotein cholesterol levels

BACKGROUND

High-density lipoprotein cholesterol (HDL) levels are a major risk factor for cardiovascular disease. Previously we identified a quantitative trait locus on baboon chromosome 18 that regulates HDL. From positional cloning studies and expression studies, we identified the endothelial lipase gene (LIPG) as the primary candidate gene for the quantitative trait locus. The mechanism by which LIPG variation influences HDL levels has not been determined.

METHODS AND RESULTS

We identified 164 LIPG polymorphisms in a panel of sibling baboons discordant for HDL1 and genotyped putative regulatory polymorphisms in a population of 951 pedigreed baboons. With the use of quantitative trait nucleotide analysis we identified 3 polymorphisms in the LIPG promoter associated with variation in serum HDL1 levels. In addition, we demonstrated that these 3 polymorphisms affect LIPG promoter activity in vitro. In silico analysis was used to identify putative transcription factors that differentially bind the functional promoter polymorphisms.

CONCLUSIONS

These results reveal LIPG variants that specifically contribute to HDL1 levels and demonstrate mechanisms by which these polymorphisms may regulate LIPG promoter activity. Results from the present study provide a mechanism, namely variation in LIPG promoter activity possibly caused by altered transcription factor binding, by which LIPG variation affects HDL levels.

A genome-wide scan for genes involved in primary vesicoureteric reflux

BACKGROUND

Vesicoureteric reflux (VUR) is the retrograde flow of urine from the bladder into the ureters. It is the most common urological anomaly in children, and a major cause of end-stage renal failure and hypertension in both children and adults. VUR is seen in approximately 1-2% of Caucasian newborns and is frequently familial.

OBJECTIVE AND METHODS

In order to search for genetic loci involved in VUR, we performed a genome-wide linkage scan using 4710 single-nucleotide polymorphisms (SNPs) in 609 individuals from 129 Irish families with >1 affected member.

RESULTS

Nonparametric linkage (NPL) analysis of the dataset yielded moderately suggestive linkage at chromosome 2q37 (NPL(max) = 2.67, p<0.001). Analysis of a subset without any additional features, such as duplex kidneys, yielded a maximum NPL score of 4.1 (p = 0.001), reaching levels of genome-wide statistical significance. Suggestive linkage was also seen at 10q26 and 6q27, and there were several smaller peaks.

CONCLUSION

Our results confirm the previous conclusion that VUR is genetically heterogeneous, and support the identification of several disease-associated regions indicated by smaller studies, as well as indicating new regions of interest for investigation.

Accounting for relatedness in family based genetic association studies

OBJECTIVE

Assess the differences in point estimates, power and type 1 error rates when accounting for and ignoring family structure in genetic tests of association.

METHODS

We compare by simulation the performance of analytic models using variance components to account for family structure and regression models that ignore relatedness for a range of possible family based study designs (i.e., sib pairs vs. large sibships vs. nuclear families vs. extended families).

RESULTS

Our analyses indicate that effect size estimates and power are not significantly affected by ignoring family structure. Type 1 error rates increase when family structure is ignored, as density of family structures increases, and as trait heritability increases. For discrete traits with moderate levels of heritability and across many common sampling designs, type 1 error rates rise from a nominal 0.05 to 0.11.

CONCLUSION

Ignoring family structure may be useful in screening although it comes at a cost of a increased type 1 error rate, the magnitude of which depends on trait heritability and pedigree configuration.

A linkage and family-based association analysis of a potential neurocognitive endophenotype of bipolar disorder

The identification of the genetic variants underpinning bipolar disorder (BPD) has been impeded by a complex pattern of inheritance characterized by genetic and phenotypic heterogeneity, genetic epistasis, and gene-environment interactions. In this paper two strategies were used to ameliorate these confounding factors. A unique South African sample including 190 individuals of the relatively, reproductively isolated Afrikaner population was assessed with a battery of neuropsychological tests in an attempt to identify a BPD-associated quantitative trait or endophenotype. BPD individuals performed significantly worse than their unaffected relatives on visual and verbal memory tasks, a finding congruent with the literature. Afocused linkage and family-based association study was carried out using this memory-related endophenotype. In the largest 77-strong Afrikaner pedigree significant evidence for linkage was detected on chromosome 22q11, a region previously implicated in BPD. The quantitative transmission disequilibrium tests-based association analysis suggested that functional variants of the DRD4 and MAO-A genes modulate memory-related cognition. We speculate that polymorphisms at these loci may predispose to a subtype of BPD characterized by memory-related deficits.

Effect of the BDNF V166M polymorphism on working memory in healthy adolescents

Brain-derived neurotrophic factor (BDNF) may play a role in modulating memory function and there is growing evidence that the BDNF V166M polymorphism may influence episodic memory in humans. However, previous association studies examining this polymorphism and working memory are inconsistent. The current study examined this association in a large sample of adolescent twin-pairs and siblings (785 individuals from 439 families). A range of measures (event-related potential, general performance and reaction time) was obtained from a delayed-response working-memory task and total association was examined using the quantitative transmission disequilibrium tests (QTDT) program. Analyses had approximately 93-97% power (alpha= 0.05) to detect an association accounting for as little as 2% of the variance in the phenotypes examined. Results indicated that the BDNF V166M polymorphism is not associated with variation in working memory in healthy adolescents.

Association of neuropeptide Y receptor Y5 polymorphisms with dyslipidemia in Mexican Americans

We examined the genetic association of neuropeptide Y receptor Y5 (NPY5R) single nucleotide polymorphisms (SNPs) with measures of the insulin resistance (metabolic) syndrome. We genotyped 10 NPY5R SNPs in 439 Mexican American individuals (age=43.3+/-17.3 years and BMI=30.0+/-6.7 kg/m2) distributed across 27 pedigrees from the San Antonio Family Diabetes Study and performed association analyses using the measured genotype approach as implemented in Sequential Oligogenic Linkage Analysis Routines (SOLAR). Minor alleles for five (rs11100493, rs12501691, P1, rs11100494, rs12512687) of the NPY5R SNPs were found to be significantly (p<0.05) associated with fasting plasma triglyceride concentrations and decreased high-density lipoprotein concentrations. In addition, the minor allele for SNP P2 was significantly associated (p=0.031) with a decreased homeostasis model assessment of beta-cell function (HOMA-%beta). Linkage disequilibrium between SNP pairs indicated one haplotype block of five SNPs (rs11100493, rs12501691, P1, rs11100494, rs12512687) that were highly correlated (r2>0.98). These preliminary results provide evidence for association of SNPs in the NPY5R gene with dyslipidemia (elevated triglyceride concentrations and reduced high-density lipoprotein levels) in our Mexican American population.

An investigation of genome-wide associations of hypertension with microsatellite markers in the family blood pressure program (FBPP)

The Family Blood Pressure Program (FBPP) has data on 387 microsatellite markers in 13,524 subjects from four major ethnic groups. We investigated genetic association with hypertension of the linkage markers. Family-based methods were used to test association of the 387 loci with resting blood pressures (BPs) [systolic blood pressure (SBP) and diastolic blood pressure (DBP)] and the hypertension status (HT). We applied a vote-counting approach to pool results across the three correlated traits, network samples, and ethnic groups to refine the selection of susceptibility loci. The association analyses captured signals missed by previous linkage scans. We found 71 loci associated with at least one of the three traits in at least one of the four ethnic groups at the significance level of 0.01. After validation across multiple samples and related traits, we identified by vote-counting 21 candidate loci for hypertension. Two loci, D3S2459 and D10S1412 confirmed findings in Network-specific linkage scans (GENOA and SAPPHIRe). Many of the candidate loci were reported by others in linkage to BPs, body weight, heart disease, and diabetes. We also observed frequent presence of quantitative trait loci (QTLs) involved in autoimmune and neurological disorders (e.g., NOD2). The vote-counting method of pooling results recognizes the potential that a gene may be involved in varying ways among different samples, which we believe is responsible for identifying genes in the less explored inflammatory pathways to hypertension.

Interpretation of simultaneous linkage and family-based association tests in genome screens

Linkage and association analyses have played important roles in identifying susceptibility genes for complex diseases. Linkage tests and family-based tests of association are often applied in the same data to help fine-map disease loci or validate results. This paradigm increases efficiency by making maximal use of family data sets. However, it is not intuitively clear under what conditions association and linkage tests performed in the same data set may be correlated. Understanding this relationship is important for interpreting the combined results of both tests. We used computer simulations and theoretical statements to estimate the correlation between linkage statistics (affected sib pair maximum LOD scores) and family-based association statistics (pedigree disequilibrium test (PDT) and association in the pressure of linkage (APL)) under various hypotheses. Different types of pedigrees were studied: nuclear families with affected sib pairs, extended pedigrees and incomplete pedigrees. Both simulation and theoretical results showed that when there is no linkage or no association, the linkage and association tests are not correlated. When there is linkage and association in the data, the two tests have a positive correlation. We concluded that when linkage and association tests are applied in the same data, the type I error rate of neither test will be affected and that power can be increased by applying tests conditionally.

Lateralization of hand skill in bipolar affective disorder

Diverse strands of evidence suggest that schizophrenia is associated with an excess of left and mixed handedness, reflecting anomalous cerebral lateralization. Genetic studies have indicated a degree of overlap between bipolar disorder (BPD) and schizophrenia. Nevertheless, pattern of handedness and degree of lateralization have not been explicitly tested in BPD. We measured handedness, footedness and relative manual dexterity in a sample of 47 families comprising BPD probands and their bipolar-spectrum and unaffected relatives (N = 240). The BPD I sample (N = 55) was significantly more lateralized on handedness, footedness and relative manual dexterity than their unaffected relatives (N = 66). They were also more lateralized than their relatives with other psychiatric diagnoses. No evidence of excess mixed handedness or footedness was observed in the BPD I sample. We raise the possibility that schizophrenia and BPD I differ in that disproportionate left-hemisphere dominance in BPD I is associated with right-hemisphere dysfunction leading to deficits in emotional regulation. Given our results, we hypothesized that degree of lateralization may be a phenotypic marker or endophenotype for BPD I. We therefore conducted a family-based genetic association analysis with this quantitative trait. Relative hand skill was significantly associated with a functional variant in the catechol-O-methyltransferase gene. We speculate that this polymorphism may influence brain lateralization.

P2 promoter variants of the hepatocyte nuclear factor 4alpha gene are associated with type 2 diabetes in Mexican Americans

Common and rare variants of the hepatocyte nuclear factor 4alpha (HNF4A) gene have been associated with type 2 diabetes and related traits in several populations suggesting the involvement of this transcription factor in diabetes pathogenesis. Single nucleotide polymorphisms (SNPs) within a large haplotype block surrounding the alternate P2 promoter, located approximately 45 kb upstream from the coding region, have been investigated in several populations of varying ethnicity with inconsistent results. Additionally, SNPs located within the P1 promoter and coding region have also been inconsistently associated with type 2 diabetes. Characterization of variation across this gene region in Mexican-American populations has not been reported. We therefore examined polymorphisms across the HNF4A gene in a cohort of Mexican-American pedigrees and assessed their association with type 2 diabetes. We observed evidence for association of SNPs in the P2 promoter region with type 2 diabetes (P = 0.003) and its age at diagnosis (P = 0.003). The risk allele frequency (53%) was intermediate to that reported in Caucasian populations (20-27%) and Pima Indians (83%). No other SNPs were associated with either trait. These results support the possibility that a variant in the P2 promoter region of HNF4A, or variants in linkage disequilibrium within this region, contributes to susceptibility to type 2 diabetes in many ethnic populations including Mexican Americans.

Haplotypes of transcription factor 7-like 2 (TCF7L2) gene and its upstream region are associated with type 2 diabetes and age of onset in Mexican Americans

TCF7L2 acts as both a repressor and transactivator of genes, as directed by the Wnt signaling pathway. Recently, several highly correlated sequence variants located within a haplotype block of the TCF7L2 gene were observed to associate with type 2 diabetes in three Caucasian cohorts. We previously reported linkage of type 2 diabetes in the San Antonio Family Diabetes Study (SAFADS) cohort consisting of extended pedigrees of Mexican Americans to the region of chromosome 10q harboring TCF7L2. We therefore genotyped 11 single nucleotide polymorphisms (SNPs) from nine haplotype blocks across the gene in 545 SAFADS subjects (178 diabetic) to investigate their role in diabetes pathogenesis. We observed nominal association between four SNPs (rs10885390, rs7903146, rs12255372, and rs3814573) in three haplotype blocks and type 2 diabetes, age at diagnosis, and 2-h glucose levels (P = 0.001-0.055). Furthermore, we identified a common protective haplotype defined by these four SNPs that was significantly associated with type 2 diabetes and age at diagnosis (P = 4.2 x 10(-5), relative risk [RR] 0.69; P = 6.7 x 10(-6), respectively) and a haplotype that confers diabetes risk that contains the rare alleles at SNPs rs10885390 and rs12255372 (P = 0.02, RR 1.64). These data provide evidence that variation in the TCF7L2 genomic region may affect risk for type 2 diabetes in Mexican Americans, but the attributable risk may be lower than in Caucasian populations.

Family Based Association Analyses between the Serotonin Transporter Gene Polymorphism (5-HTTLPR) and Neuroticism, Anxiety and Depression

We studied the association between the short/long promotor-based length polymorphism of the serotonin transporter gene (5-HTTLPR) and neuroticism, anxiety and depression. Subjects included twins, their siblings and parents from the Netherlands Twin Register (559 parents and 1,245 offspring). Subjects had participated between one and five times in a survey study measuring neuroticism, anxiety and depression. Offspring of these families were also approached to participate in a psychiatric interview diagnosing DSM-IV major depression. Within-family and total association between 5-HTTLPR and these traits were tested. Only three of the 36 tests showed a significant effect of 5-HTTLPR (P<0.05). These effects were in opposite directions, i.e. both negative and positive regression coefficients were found for the s allele. No additive effect of the s allele was found for DSM-IV depression. Our results strongly suggest that there is no straightforward association between 5-HTTLPR and neuroticism, anxiety and depression.

X-APL: an improved family-based test of association in the presence of linkage for the X chromosome

Family-based association methods have been developed primarily for autosomal markers. The X-linked sibling transmission/disequilibrium test (XS-TDT) and the reconstruction-combined TDT for X-chromosome markers (XRC-TDT) are the first association-based methods for testing markers on the X chromosome in family data sets. These are valid tests of association in family triads or discordant sib pairs but are not theoretically valid in multiplex families when linkage is present. Recently, XPDT and XMCPDT, modified versions of the pedigree disequilibrium test (PDT), were proposed. Like the PDT, XPDT compares genotype transmissions from parents to affected offspring or genotypes of discordant siblings; however, the XPDT can have low power if there are many missing parental genotypes. XMCPDT uses a Monte Carlo sampling approach to infer missing parental genotypes on the basis of true or estimated population allele frequencies. Although the XMCPDT was shown to be more powerful than the XPDT, variability in the statistic due to the use of an estimate of allele frequency is not properly accounted for. Here, we present a novel family-based test of association, X-APL, a modification of the test for association in the presence of linkage (APL) test. Like the APL, X-APL can use singleton or multiplex families and properly infers missing parental genotypes in linkage regions by considering identity-by-descent parameters for affected siblings. Sampling variability of parameter estimates is accounted for through a bootstrap procedure. X-APL can test individual marker loci or X-chromosome haplotypes. To allow for different penetrances in males and females, separate sex-specific tests are provided. Using simulated data, we demonstrated validity and showed that the X-APL is more powerful than alternative tests. To show its utility and to discuss interpretation in real-data analysis, we also applied the X-APL to candidate-gene data in a sample of families with Parkinson disease.

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.

Type I error rates in association versus joint linkage/association tests in related individuals

Positional gene discovery on pedigree data typically involves initial gross localization by linkage analysis with subsequent finer localization by association analysis in areas that show evidence of linkage. We examine the effect of including linkage information when testing for association in the context of variance-components-based pedigree analysis. We present simulation experiments showing that, at least in the extreme case of a rare private allele, failing to include the linkage variance component in the association model results in excessive Type I error that increases with allele copy number and/or quantitative trait locus (QTL) effect size. Joint estimation of the linkage variance component in the association model reduces Type I error to nominal expectations. This holds whether allele-sharing probabilities are estimated from a polymorphic marker or from the very single-nucleotide polymorphism (SNP) being tested for association, although the latter provides much less power. These results support the idea that an appropriate association analysis must test both the random effect of shared marker alleles (linkage) and the mean effects of the marker genotypes (association).

No contribution of angiotensin-converting enzyme (ACE) gene variants to severe obesity: a model for comprehensive case/control and quantitative cladistic analysis of ACE in human diseases

Candidate gene analyses are often inconclusive owing to genetic or phenotypic heterogeneity, low statistical power, selection of nonfunctional SNPs, and inadequate statistical analysis of the genetic architecture. Angiotensin-converting enzyme (ACE) is involved in adipocyte growth and function and the ACE-processed angiotensin II inhibits adipocyte differentiation. Associations between body mass index (BMI) and ACE polymorphisms have been reported in general populations, but the contribution to severe obesity of this gene, which is located under an obesity genome-scan linkage peak on 17q23, is unknown. ACE is one of the most studied genes and markers responsible for variation in circulating ACE enzyme levels have been extensively characterised. Eight of these variants were genotyped in 1054 severely obese cases and 918 nonobese controls, as well as 116 nuclear families from the genome scan (n=447), enabling the known clades to be inferred. Qualitative analysis of individual single-nucleotide polymorphisms (SNPs), haplotypes, clades, and diploclades demonstrated no significant associations (P<0.05) after minimal correction for multiple testing. Quantitative analysis of clades and diploclades for BMI, waist-to-hip ratio, or ZBMI in children were also not significant. This rigorous, large-scale study of common, well-defined, severe polygenic obesity provides strong evidence that functionally relevant sequence variation in ACE, whether it is defined at the level of SNPs, haplotypes, or clades, is not associated with severe obesity in French Caucasians. Such a study design exemplifies the strategy needed to clearly define the contribution of the ACE gene to the plethora of complex genetic diseases where weak associations have been previously reported.

Autosomal linkage analysis for the level of response to alcohol

BACKGROUND

The level of response (LR) to alcohol is a genetically-influenced phenotype related to the alcoholism risk. Usually measured by evaluating psychological and physiological changes that follow the administration of alcohol, the heritability of LR is estimated to be between 0.4 and 0.6, and efforts are being made to find genes related to this phenotype. This paper presents data from a family-based genome with linkage analysis focusing on alcohol challenge determinants of LR.

METHODS

The subjects were 18-to-29-year-old sibling pairs with at least one parent who was alcohol-dependent and who had experience with alcohol but were not yet alcohol-dependent themselves. Both members of the sibling pairs were given oral alcohol challenges (0.75-0.90 ml/kg of ethanol for females and males, respectively), with LR established using the Subjective High Assessment Scale (SHAS) and changes in body sway (BS) repeatedly over a 3.5-hr. period. Blood samples from siblings and at least one parent were genotyped using 811 microsatellite markers, with results evaluated using several related variance component approaches as implemented in SOLAR for continuous traits. In addition, association was tested using single nucleotide polymorphisms (SNPs) within the KCNMA1, HTR7 and SLC18A2 genes that may relate to a finding on chromosome 10.

RESULTS

Data were generated from 238 sib-pairs representing 365 individuals (41.6% were males) from 165 families. The most consistent results across methods and samples were observed for SHAS on chromosome 10 between 120 and 140 cM (with a maximum LOD score of 2.6 at 122 cM), and a second region of possible interest at 173 cM (LOD = 1.2). Statistical analysis with the KCNMA1, HTR7 and SLC18A2 genes, which lie in the support region of interest revealed no evidence for association after correction for multiple comparisons.

CONCLUSIONS

These evaluations from the largest known alcohol challenge-based genetic study to date highlight the potential importance of genes on chromosome 10 as possible contributors to the low LR to alcohol as a risk factor for alcoholism.

Association between variants in the genes for adiponectin and its receptors with insulin resistance syndrome (IRS)-related phenotypes in Mexican Americans

AIMS/HYPOTHESIS

The aim of this study was to examine whether genetic variation in ADIPOQ, ADIPOR1 and ADIPOR2 may contribute to increased susceptibility to components of the insulin resistance syndrome (IRS).

MATERIALS AND METHODS

We genotyped single-nucleotide polymorphisms (SNPs) in ADIPOQ, ADIPOR1 and ADIPOR2 in Mexican American subjects (N=439) and performed an association analysis of IRS-related traits.

RESULTS

Of the eight SNPs examined in the ADIPOQ gene, rs4632532 and rs182052 exhibited significant associations with BMI (p=0.029 and p=0.032), fasting specific insulin (p=0.023 and p=0.026), sum of skin folds (SS) (p=0.0089 and p=0.0084) and homeostasis model assessment of insulin sensitivity (HOMA-%S) (p=0.015 and p=0.016). Two other SNPs, rs266729 and rs2241767, were significantly associated with SS (p=0.036 and p=0.013). SNP rs7539542 of ADIPOR1 was significantly associated with BMI, SS and waist circumference (p=0.025, p=0.047 and p=0.0062). Fourteen of the ADIPOR2 SNPs were found to be significantly (p<0.05) associated with fasting plasma triglyceride concentrations. Four of these SNPs (rs10848569, rs929434, rs3809266 and rs12342) were in high pairwise linkage disequilibrium (r (2)=0.99) and were strongly associated with fasting triglyceride levels (p=0.00029, p=0.00016, p=0.00027 and p=0.00021). Adjusting for the effects of BMI and HOMA-%S on triglyceride concentrations increased significance to p=0.000060 for SNP rs929434. Bayesian quantitative trait nucleotide analysis was used to examine all possible models of gene action. Again, SNP rs929434 provided the strongest statistical evidence of an effect on triglyceride concentrations.

CONCLUSIONS/INTERPRETATION

These results provide evidence for association of SNPs in ADIPOQ and its receptors with multiple IRS-related phenotypes. Specifically, several genetic variants in ADIPOR2 were strongly associated with decreased triglyceride levels.

A single nucleotide polymorphism fine mapping study of chromosome 1q42.1 reveals the vulnerability genes for schizophrenia, GNPAT and DISC1: Association with impairment of sustained attention

BACKGROUND

The marker D1S251 of chromosome 1q42.1 showed significant association with schizophrenia in a Taiwanese sample. We used single nucleotide polymorphism (SNP) fine mapping to search for the vulnerability genes of schizophrenia.

METHODS

We selected 120 SNPs covering 1 Mb around D1S251 from the public database. These selected SNPs were initially validated if allele frequency was >10%. Forty-seven validated SNPs were genotyped in 102 families with at least 2 siblings affected with schizophrenia.

RESULTS

Two SNP blocks showed significant association with schizophrenia. Block 1 (five-SNP), located between intron 2 and intron 13 of the glyceronephosphate O-acyltransferase (GNPAT) gene, showed the most significant associations using single-locus TDT (z = -2.07, p = .038, df = 1) and haplotype association analyses (z = -1.99, p = .046, df = 1). Block 2 (two-SNP), located between intron 4 and intron 5 of the disrupted-in-schizophrenia 1 (DISC1) gene, also showed the most significant results in both the single-locus (z = -3.22, p = .0013, df = 1) and haplotype association analyses (z = 3.35, p = .0008, df = 1). The association of the DISC1 gene with schizophrenia was mainly in the patient group with sustained attention deficits as assessed by the Continuous Performance Test.

CONCLUSIONS

Chromosome 1q42.1 harbors GNPAT and DISC1 as candidate genes for schizophrenia, and DISC1 is associated with sustained attention deficits.