Pedigree tests of transmission disequilibrium

Integrating social science and genetics: news from the political front

There has been growing interest in the use of genetic models to expand the understanding of political preferences, attitudes, and behaviors. Researchers in the social sciences have begun incorporating these models and have revealed that genetic differences account for individual differences in political beliefs, behaviors, and responses to the political environment. The first Integrating Genetics and the Social Sciences Conference, held at Boulder, Colorado in May of 2010, brought together these researchers. As a result, we jointly review the last 5 years of research in this area. In doing so, we explicate the methods, findings, and limitations of behavior genetic approaches, including twin designs, association studies, and genome-wide analyses, in their application toward exploring political preferences.

Genome-wide association study of theta band event-related oscillations identifies serotonin receptor gene HTR7 influencing risk of alcohol dependence

Event-related brain oscillations (EROs) represent highly heritable neuroelectrical correlates of human perception and cognitive performance that exhibit marked deficits in patients with various psychiatric disorders. We report the results of the first genome-wide association study (GWAS) of an ERO endophenotype-frontal theta ERO evoked by visual oddball targets during P300 response in 1,064 unrelated individuals drawn from a study of alcohol dependence. Forty-two SNPs of the Illumina HumanHap 1 M microarray were selected from the theta ERO GWAS for replication in family-based samples (N = 1,095), with four markers revealing nominally significant association. The most significant marker from the two-stage study is rs4907240 located within ARID protein 5A gene (ARID5A) on chromosome 2q11 (unadjusted, Fisher's combined P = 3.68 × 10⁻⁶). However, the most intriguing association to emerge is with rs7916403 in serotonin receptor gene HTR7 on chromosome 10q23 (combined P = 1.53 × 10⁻⁴), implicating the serotonergic system in the neurophysiological underpinnings of theta EROs. Moreover, promising SNPs were tested for association with diagnoses of alcohol dependence (DSM-IV), revealing a significant relationship with the HTR7 polymorphism among GWAS case-controls (P = 0.008). Significant recessive genetic effects were also detected for alcohol dependence in both case-control and family-based samples (P = 0.031 and 0.042, respectively), with the HTR7 risk allele corresponding to theta ERO reductions among homozygotes. These results suggest a role of the serotonergic system in the biological basis of alcohol dependence and underscore the utility of analyzing brain oscillations as a powerful approach to understanding complex genetic psychiatric disorders.

Statistical issues in gene association studies

This chapter reviews statistical issues related to gene association studies. The goal is to review various aspects of study design and analysis for individuals who do not have an extensive statistical background. We will review statistical issues as they relate to both genome-wide and candidate gene studies. Topics reviewed include study design, power and sample size, data checking, statistical methods, population stratification, and multiple testing. We draw examples from the type 2 diabetes genetics literature to illustrate some of the issues discussed.

Pleiotropic effects of DCDC2 and DYX1C1 genes on language and mathematics traits in nuclear families of developmental dyslexia

Converging evidence indicates that developmental problems in oral language and mathematics can predate or co-occur with developmental dyslexia (DD). Substantial genetic correlations have been found between language, mathematics and reading traits, independent of the method of sampling. We tested for association of variants of two DD susceptibility genes, DCDC2 and DYX1C1, in nuclear families ascertained through a proband with DD using concurrent measurements of language and mathematics in both probands and siblings by the Quantitative Transmission Disequilibrium Test. Evidence for significant associations was found between DCDC2 and 'Numerical Facts' (p value = 0.02, with 85 informative families, genetic effect = 0.57) and between 'Mental Calculation' and DYX1C1 markers -3GA (p value = 0.05, with 40 informative families, genetic effect = -0.67) and 1249GT (p value = 0.02, with 49 informative families, genetic effect = -0.65). No statistically significant associations were found between DCDC2 or DYX1C1 and language phenotypes. Both DCDC2 and DYX1C1 DD susceptibility genes appear to have a pleiotropic role on mathematics but not language phenotypes.

Identification of candidate genes for dyslexia susceptibility on chromosome 18

BACKGROUND

Six independent studies have identified linkage to chromosome 18 for developmental dyslexia or general reading ability. Until now, no candidate genes have been identified to explain this linkage. Here, we set out to identify the gene(s) conferring susceptibility by a two stage strategy of linkage and association analysis.

METHODOLOGY/PRINCIPAL FINDINGS

Linkage analysis: 264 UK families and 155 US families each containing at least one child diagnosed with dyslexia were genotyped with a dense set of microsatellite markers on chromosome 18. Association analysis: Using a discovery sample of 187 UK families, nearly 3000 SNPs were genotyped across the chromosome 18 dyslexia susceptibility candidate region. Following association analysis, the top ranking SNPs were then genotyped in the remaining samples. The linkage analysis revealed a broad signal that spans approximately 40 Mb from 18p11.2 to 18q12.2. Following the association analysis and subsequent replication attempts, we observed consistent association with the same SNPs in three genes; melanocortin 5 receptor (MC5R), dymeclin (DYM) and neural precursor cell expressed, developmentally down-regulated 4-like (NEDD4L).

CONCLUSIONS

Along with already published biological evidence, MC5R, DYM and NEDD4L make attractive candidates for dyslexia susceptibility genes. However, further replication and functional studies are still required.

Genetic and molecular functional characterization of variants within TNFSF13B, a positional candidate preeclampsia susceptibility gene on 13q

BACKGROUND

Preeclampsia is a serious pregnancy complication, demonstrating a complex pattern of inheritance. The elucidation of genetic liability to preeclampsia remains a major challenge in obstetric medicine. We have adopted a positional cloning approach to identify maternal genetic components, with linkages previously demonstrated to chromosomes 2q, 5q and 13q in an Australian/New Zealand familial cohort. The current study aimed to identify potential functional and structural variants in the positional candidate gene TNFSF13B under the 13q linkage peak and assess their association status with maternal preeclampsia genetic susceptibility.

METHODOLOGY/PRINCIPAL FINDINGS

The proximal promoter and coding regions of the positional candidate gene TNFSF13B residing within the 13q linkage region was sequenced using 48 proband or founder individuals from Australian/New Zealand families. Ten sequence variants (nine SNPs and one single base insertion) were identified and seven SNPs were successfully genotyped in the total Australian/New Zealand family cohort (74 families/480 individuals). Borderline association to preeclampsia (p = 0.0153) was observed for three rare SNPs (rs16972194, rs16972197 and rs56124946) in strong linkage disequilibrium with each other. Functional evaluation by electrophoretic mobility shift assays showed differential nuclear factor binding to the minor allele of the rs16972194 SNP, residing upstream of the translation start site, making this a putative functional variant. The observed genetic associations were not replicated in a Norwegian case/control cohort (The Nord-Trøndelag Health Study (HUNT2), 851 preeclamptic and 1,440 non-preeclamptic women).

CONCLUSION/SIGNIFICANCE

TNFSF13B has previously been suggested to contribute to the normal immunological adaption crucial for a successful pregnancy. Our observations support TNFSF13B as a potential novel preeclampsia susceptibility gene. We discuss a possible role for TNFSF13B in preeclampsia pathogenesis, and propose the rs16972194 variant as a candidate for further functional evaluation.

Polymorphic cis- and trans-regulation of human gene expression

Using genetic and molecular analyses, we identified over 1,000 polymorphic regulators that regulate expression levels of human genes.

Expression levels of human genes vary extensively among individuals. This variation facilitates analyses of expression levels as quantitative phenotypes in genetic studies where the entire genome can be scanned for regulators without prior knowledge of the regulatory mechanisms, thus enabling the identification of unknown regulatory relationships. Here, we carried out such genetic analyses with a large sample size and identified cis- and trans-acting polymorphic regulators for about 1,000 human genes. We validated the cis-acting regulators by demonstrating differential allelic expression with sequencing of transcriptomes (RNA-Seq) and the trans-regulators by gene knockdown, metabolic assays, and chromosome conformation capture analysis. The majority of the regulators act in trans to the target (regulated) genes. Most of these trans-regulators were not known to play a role in gene expression regulation. The identification of these regulators enabled the characterization of polymorphic regulation of human gene expression at a resolution that was unattainable in the past.

Cellular characteristics and functions are determined largely by gene expression and expression levels differ among individuals, however it is not clear how these levels are regulated. While many cis-acting DNA sequence variants in promoters and enhancers that influence gene expression have been identified, only a few polymorphic trans-regulators of human genes are known. Here, we used human B-cells from individuals belonging to large families and identified polymorphic trans-regulators for about 1,000 human genes. We validated these results by gene knockdown, metabolic perturbation studies and chromosome conformation capture assays. Although these regulatory relationships were identified in cultured B-cells, we show that some of the relationships were also found in primary fibroblasts. The large number of regulators allowed us to better understand gene expression regulation, to uncover new gene functions, and to identify their roles in disease processes. This study shows that genetic variation is a powerful tool not only for gene mapping but also to study gene interaction and regulation.

The R563Q mutation of the epithelial sodium channel beta-subunit is associated with hypertension

BACKGROUND

A high prevalence of the R563Q mutation of the epithelial sodium channel β-subunit has been reported in South African hypertensives compared with unrelated normotensive controls. To delineate the effects of this mutation against a more uniform genetic background, this study investigated the association of the mutation with hypertension within affected kindreds.

METHODS

Forty-five index patients and members of their kindreds were studied. Blood pressure, serum potassium and the presence of the R563Q mutation were determined.

RESULTS

Of the 136 individuals studied, 89 were heterozygous for the R563Q mutation and 47 homozygous RR. The mean arterial pressure was significantly higher in the R563Q heterozygous group (p = 0.005) after adjusting for gender, race, age and kindred membership. Of the R563Q heterozygous subjects, 71 (80%) had hypertension, while 17 (36%) of the R563Q homozygous RR subjects were hypertensive. Six R563Q heterozygous subjects had hypokalaemia and one R563Q homozygous RR subject had hypokalaemia, but the difference was not statistically significant. Two heterozygous patients had Liddle's syndrome, both occurring during pregnancy.

CONCLUSION

The R563Q mutation of β-ENaC is associated with hypertension within affected kindreds, but does not usually cause the full Liddle's syndrome phenotype.

Genetic variants in SLC9A9 are associated with measures of attention-deficit/hyperactivity disorder symptoms in families

OBJECTIVE

A family was previously identified that cosegregates a pericentric inversion, inv(3)(p14 : q21), with an early-onset developmental condition, characterized by impulsive behavior and intellectual deficit. The inversion breakpoints lie within DOCK3 and SLC9A9 at the p-arm and q-arm, respectively. Based on this report, these genes were selected to be evaluated in a family-based attention-deficit/hyperactivity disorder (AD/HD) association study.

METHODS

Conners' Parent (CPRS) and Teacher (CTRS) Rating Scales of AD/HD symptoms and Conners' Continuous Performance Test (CPT) measures were collected and a minimal number of tagging single-nucleotide polymorphisms (SNPs) in each gene were selected for analysis. Analyses were performed on families who met research criteria for AD/HD. Using the program, QTDT, each tagging SNP was tested for association with T-scores from the Diagnostic and Statistical Manual of Mental Disorders, fourth edition (DSM-IV) subscales according to the CTRS and CPRS, and five CPT measures.

RESULTS

After adjusting for multiple testing, a SNP in the 3' UTR of SLC9A9, rs1046706, remained significantly associated (false discovery rate, q value <0.05) with scores on the DSM-IV hyperactive-impulsive and total symptom subscales according to the CTRS and errors of commission on the CPT. In addition, an intronic SLC9A9 SNP, rs2360867, remained significantly associated with errors of commission.

CONCLUSION

Our results suggest that SLC9A9 may be related to hyperactive-impulsive symptoms in AD/HD and the disruption of SLC9A9 may be responsible for the behavioral phenotype observed in the inversion family. The association with SLC9A9 is particularly interesting as it was recently implicated in a genome-wide association study for AD/HD. Further investigation of the role of SLC9A9 in AD/HD and other behavioral disorders is warranted.

ALOX12 polymorphisms are associated with fat mass but not peak bone mineral density in Chinese nuclear families

OBJECTIVE

Arachidonate 12-lipoxygenase (ALOX12) is a member of the lipoxygenase superfamily, which catalyzes the incorporation of molecular oxygen into polyunsaturated fatty acids. The products of ALOX12 reactions serve as endogenous ligands for peroxisome proliferator-activated receptor γ (PPARG). The activation of the PPARG pathway in marrow-derived mesenchymal progenitors stimulates adipogenesis and inhibits osteoblastogenesis. Our objective was to determine whether polymorphisms in the ALOX12 gene were associated with variations in peak bone mineral density (BMD) and obesity phenotypes in young Chinese men.

METHODS

All six tagging single-nucleotide polymorphisms (SNPs) in the ALOX12 gene were genotyped in a total of 1215 subjects from 400 Chinese nuclear families by allele-specific polymerase chain reaction. The BMD at the lumbar spine and hip, total fat mass (TFM) and total lean mass (TLM) were measured using dual-energy X-ray absorptiometry. The pairwise linkage disequilibrium among SNPs was measured, and the haplotype blocks were inferred. Both the individual SNP markers and the haplotypes were tested for an association with the peak BMD, body mass index, TFM, TLM and percentage fat mass (PFM) using the quantitative transmission disequilibrium test (QTDT).

RESULTS

Using the QTDT, significant within-family association was found between the rs2073438 polymorphism in the ALOX12 gene and the TFM and PFM (P=0.007 and 0.012, respectively). Haplotype analyses were combined with our individual SNP results and remained significant even after correction for multiple testing. However, we failed to find significant within-family associations between ALOX12 SNPs and the BMD at any bone site in young Chinese men.

CONCLUSIONS

Our present results suggest that the rs2073438 polymorphism of ALOX12 contributes to the variation of obesity phenotypes in young Chinese men, although we failed to replicate the association with the peak BMD variation in this sample. Further independent studies are needed to confirm our findings.

Polymorphisms in the HOXD4 gene are not associated with peak bone mineral density in Chinese nuclear families

AIM

To determine the associations between HOXD4 gene polymorphisms with peak bone mineral density (BMD) throughing measuring three tagging single nucleotide polymorphisms (tagSNPs), including rs1867863, rs13418078, and rs4972504, in HOXD4.

METHODS

Four hundred Chinese nuclear families with male offspring (1215 subjects) and 401 Chinese nuclear families with female offspring (1260 subjects) were recruited. BMD of the lumbar spine 1-4 (L1-4) and left proximal femur including total hip and femoral neck were measured by dual-energy X-ray absorptiometry. The quantitative transmission disequilibrium test (QTDT) was performed to investigate the association among the tagging SNPs, haplotypes and peak BMD.

RESULTS

Only the CC genotype was identified in rs13418078 in the Chinese population, unlike other populations. We failed to find significant within-family association among these SNPs, haplotypes and peak BMD at any bone site in either male- or female-offspring nuclear families.

CONCLUSION

The results suggest that genetic polymorphisms in HOXD4 may not be a major contributor to the observed variability in peak BMD in the lumbar spine and the hip in Chinese men and women.

Variance component methods for analysis of complex phenotypes

Variance component methods have a long history in human quantitative genetics as well as in agricultural genetics and animal breeding. They are designed for genetic analysis of continuously varying quantitative traits such as body mass index (BMI), cholesterol levels, or intelligence quotient. They can be used to assess the strength of genetic effects on a trait, to localize genes influencing a trait through either linkage or association methods, to assess whether associated variants are likely to be the functional variants behind a given localization signal, to explore whether related traits have shared genetic influences in multivariate analyses, and to characterize the genetic effects on a trait through analyses of gene-gene and gene-environment interactions.

The Ser(326)Cys polymorphism of 8-oxoguanine glycosylase 1 (OGG1) is associated with type 2 diabetes in Mexican Americans

OBJECTIVE

Human 8-oxoguanine glycosylase 1 (OGG1) excises oxidatively damaged promutagenic base 8-oxoguanine, a lesion previously observed in a rat model of type 2 diabetes (T2DM). The objective of the present study is to determine whether genetic variation in OGG1 is associated with type 2 diabetes (T2DM) in a Mexican American cohort.

METHODS

Ten SNPs including two tagging SNPs (rs1052133, rs2072668) across the OGG1 gene region were selected from the Hapmap database and genotyped in the entire cohort (n = 670; 29% diabetes; 39 families) by TaqMan assay. Association analyses between the SNPs and T2DM were performed using the measured genotype approach as implemented in the program SOLAR.

RESULTS

Of the ten SNPs genotyped, only five were polymorphic. The minor allele frequencies of these 5 SNPs ranged from 1-38%. Of the SNPs examined for association, the Ser(326)Cys (rs1052133) exhibited significant association with T2DM (p = 0.016) after accounting for age and sex effects. Another intronic variant (rs2072668), which was in strong linkage disequilibrium (r(2) = 0.96) with Ser(326)Cys also exhibited significant association with T2DM (p = 0.031).

CONCLUSIONS

These results suggest for the first time that the variants in OGG1 could influence diabetes risk in these Mexican American families and support a role for alterations of OGG1 in the pathogenesis of T2DM.

The characterization of Abelson helper integration site-1 in skeletal muscle and its links to the metabolic syndrome

The human Abelson helper integration site-1 (AHI1) gene is associated with both neurologic and hematologic disorders; however, it is also located in a chromosomal region linked to metabolic syndrome phenotypes and was identified as a type 2 diabetes mellitus susceptibility gene from a genomewide association study. To further define a possible role in type 2 diabetes mellitus development, AHI1 messenger RNA expression levels were investigated in a range of tissues and found to be highly expressed in skeletal muscle as well as displaying elevated levels in brain regions and gonad tissues. Further analysis in a rodent polygenic animal model of obesity and type 2 diabetes mellitus identified increased Ahi-1 messenger RNA levels in red gastrocnemius muscle from fasted impaired glucose-tolerant and diabetic rodents compared with healthy animals (P < .002). Moreover, elevated gene expression levels were confirmed in skeletal muscle from fasted obese and type 2 diabetes mellitus human subjects (P < .02). RNAi-mediated suppression of Ahi-1 resulted in increased glucose transport in rat L6 myotubes in both the basal and insulin-stimulated states (P < .01). Finally, single nucleotide polymorphism association studies identified 2 novel AHI1 genetic variants linked with fasting blood glucose levels in Mexican American subjects (P < .037). These findings indicate a novel role for AHI1 in skeletal muscle and identify additional genetic links with metabolic syndrome phenotypes suggesting an involvement of AHI1 in the maintenance of glucose homeostasis and type 2 diabetes mellitus progression.

A polymorphism in the oestrogen receptor gene explains covariance between digit ratio and mating behaviour

In vertebrates, including humans, the relative length of the second to the fourth digit correlates with sex hormone-dependent behavioural, psychological and physiological traits. However, despite a decade of research, the underlying mechanism linking digit ratio to these sex hormone-dependent traits remains unclear. Previous work suggests that during embryo development, circulating levels of plasma androgens or oestrogens may act through their receptors to affect transcription levels of posterior HOX genes in the developing digits, thereby possibly influencing their relative length. The correlation between digit ratio and sex hormone-dependent traits might thus stem from variation in expression or sensitivity of the sex hormone receptors, or from variation in sex hormone levels in the embryo. Here, we show that in a population of 1156 zebra finches Taeniopygia guttata, a polymorphism in the oestrogen receptor α gene (ESR1) explains 11.3 per cent of the variation in digit ratio, and is also associated with male and female-mating behaviour. By contrast, we found no associations between digit ratio or mating behaviours and polymorphisms in the androgen receptor gene. Thus, our results (i) provide an explanation for the observed significant genetic covariance between digit ratio and male and female mating behaviour and (ii) strongly confirm the indicator function of digit ratio through the oestrogen pathway. Finally, we note that the commonly invoked effect of foetal testosterone on human digit ratio seems to be substantially weaker than the effect described here.

Variants within the calpain-10 gene and relationships with type 2 diabetes (T2DM) and T2DM-related traits among Tunisian Arabs

BACKGROUND

Common variations in the calpain 10 (CAPN10) gene variants UCSNP-43, UCSNP-19 and UCSNP-63, and the 112/121 diplotype, are associated with an increased risk of type 2 diabetes (T2DM) and T2DM-related traits.

METHODS

The association of UCSNP-43, -19 and -63 CAPN10 SNPs with T2DM was assessed in 917 Tunisian T2DM patients and 748 ethnically matched non-diabetic controls. CAPN10 genotyping was done by PCR-RFLP.

RESULTS

Significant differences in UCSNP-19 MAF, but not UCSNP-43 or -63, and genotype distribution were seen between patients and controls. Heterogeneity in UCSNP-19, but not UCSNP-43 and -63, genotype distribution was noted according to geographical origin. Obesity was associated with UCSNP-19, while raised fasting glucose was associated with UCSNP-63, and increased HDL was associated with UCSNP-43. Enrichment of homozygous UCSNP-19 2/2 was seen in overweight and obese compared with lean patients; logistic-regression analyses demonstrated a positive association of the 2/2 genotype with overweight [P=0.003; OR (95% CI)=2.07 (1.28-3.33)] and obese [P=0.021; OR (95% CI)=1.83 (1.10-3.07)] patients. Of the six CAPN10 haplotypes identified, significant enrichment of only haplotype 111 was seen in T2DM patients [Pc=0.034; OR (95% CI)=1.22 (1.06-1.41)], while the frequency of all identified CAPN10 diplotypes, including the high-risk 112/121, was comparable between patients and controls.

CONCLUSION

While CAPN10 UCSNP-19 SNP and haplotype 111 contribute to the risk of T2DM in Tunisian subjects, no significant association between CAPN10 diplotypes and T2DM was demonstrated.

Common polymorphisms of calpain-10 and the risk of Type 2 Diabetes in a Tunisian Arab population: a case-control study

Background

Genetic variations in the calpain-10 gene (CAPN10), in particular the at-risk diplotype (112/121), were previously implicated with increased risk of type 2 diabetes (T2D).

Methods

We examined the association of CAPN10 UCSNP-43 (rs3792267), UCSNP-19 (rs3842570), and UCSNP-63 (rs5030952) SNPs with T2D in 917 Tunisian T2D patients and 748 non-diabetic controls. CAPN10 genotyping was done by PCR-RFLP.

Results

Enrichment of UCSNP-19 2R (minor) allele and 2R/2R genotype was found in T2D patients; the allele and genotype distribution of UCSNP-43 and UCSNP-63 alleles and genotypes were not significantly different between patient groups and non-diabetic control subjects. Regression analysis demonstrated progressive increases in T2D risk in 3R/2R [OR (95% CI) = 1.35 (1.08 - 1.68)] and 2R/2R [OR (95% CI) = 1.61 (1.20 - 2.18)] genotypes. Of the six haplotypes detected, enrichment of haplotype 111 (UCSNP-43/UCSNP-19/UCSNP-63) was seen in patients (Pc = 0.034); the distribution of the other haplotypes was comparable between patients and control subjects; neither haplotype 211 nor haplotype 212 was observed. Furthermore, the frequency of all CAPN10 diplotypes identified, including the "high-risk diplotype (112/121) reported for Mexican-Americans and Northern Europeans, were comparable between patients and controls.

Conclusions

CAPN10 UCSNP-19 variant, and the 111 haplotype contribute to the risk of T2D in Tunisian subjects; no significant associations between CAPN10 diplotypes and T2D were demonstrated for Tunisians.

Single-nucleotide polymorphisms in corticotropin releasing hormone receptor 1 gene (CRHR1) are associated with quantitative trait of event-related potential and alcohol dependence

BACKGROUND

Endophenotypes reflect more proximal effects of genes than diagnostic categories, hence providing a more powerful strategy in searching for genes involved in complex psychiatric disorders. There is strong evidence suggesting the P3 amplitude of the event-related potential (ERP) as an endophenotype for the risk of alcoholism and other disinhibitory disorders. Recent studies demonstrated a crucial role of corticotropin releasing hormone receptor 1 (CRHR1) in the environmental stress response and ethanol self-administration in animal models. The aim of the present study was to test the potential associations between single-nucleotide polymorphisms (SNPs) in the CRHR1 gene and the quantitative trait, P3 amplitude during the processing of visual target signals in an oddball paradigm, as well as alcohol dependence diagnosis.

METHODS

We analyzed a sample from the Collaborative Study on the Genetics of Alcoholism (COGA) comprising 1049 Caucasian subjects from 209 families (including 472 alcohol-dependent individuals). Quantitative transmission disequilibrium test (QTDT) and family-based association test (FBAT) were used to test the association, and false discovery rate (FDR) was applied to correct for multiple comparisons.

RESULTS

Significant associations (p < 0.05) were found between the P3 amplitude and alcohol dependence with multiple SNPs in the CRHR1 gene.

CONCLUSIONS

Our results suggest that CRHR1 may be involved in modulating the P3 component of the ERP during information processing and in vulnerability to alcoholism. These findings underscore the utility of electrophysiology and the endophenotype approach in the genetic study of psychiatric disorders.

Catechol-o-methyltransferase genotype and childhood trauma may interact to impact schizotypal personality traits

We attempt to identify gene by childhood abuse interactions which predispose to the development of schizotypal traits in a familial bipolar disorder (BD) sample. Self-report measures of schizotypal personality traits (Schizotypal Personality Scale) and childhood maltreatment (Childhood Trauma Questionnaire) were administered to 222 participants from 44 families with BD. Variants of catechol-o-methyltransferase (COMT) and four other dopamine pathway-related genes: DRD4, DRD2,MAOA, and SLC6A3, were typed. BD type I (BD I) subjects scored significantly higher than their unaffected relatives on the Schizotypal Personality Scale. The val allele of the Val158 Met polymorphism of the COMT gene was associated with increased schizotypal personality trait scores in individuals exposed to higher levels of self-reported childhood trauma (p < 0.05). There was no direct effect of the val158met polymorphism on schizotypal personality traits. Further, no passive correlation between COMT genotype and childhood trauma was found. We raise the possibility that genetically-driven variation in COMT may interact with childhood trauma to contribute to the risk of developing schizotypal personality traits.

Family-based bivariate association tests for quantitative traits

The availability of a large number of dense SNPs, high-throughput genotyping and computation methods promotes the application of family-based association tests. While most of the current family-based analyses focus only on individual traits, joint analyses of correlated traits can extract more information and potentially improve the statistical power. However, current TDT-based methods are low-powered. Here, we develop a method for tests of association for bivariate quantitative traits in families. In particular, we correct for population stratification by the use of an integration of principal component analysis and TDT. A score test statistic in the variance-components model is proposed. Extensive simulation studies indicate that the proposed method not only outperforms approaches limited to individual traits when pleiotropic effect is present, but also surpasses the power of two popular bivariate association tests termed FBAT-GEE and FBAT-PC, respectively, while correcting for population stratification. When applied to the GAW16 datasets, the proposed method successfully identifies at the genome-wide level the two SNPs that present pleiotropic effects to HDL and TG traits.

SCAN: SNP and copy number annotation

MOTIVATION

Genome-wide association studies (GWAS) generate relationships between hundreds of thousands of single nucleotide polymorphisms (SNPs) and complex phenotypes. The contribution of the traditionally overlooked copy number variations (CNVs) to complex traits is also being actively studied. To facilitate the interpretation of the data and the designing of follow-up experimental validations, we have developed a database that enables the sensible prioritization of these variants by combining several approaches, involving not only publicly available physical and functional annotations but also multilocus linkage disequilibrium (LD) annotations as well as annotations of expression quantitative trait loci (eQTLs).

RESULTS

For each SNP, the SCAN database provides: (i) summary information from eQTL mapping of HapMap SNPs to gene expression (evaluated by the Affymetrix exon array) in the full set of HapMap CEU (Caucasians from UT, USA) and YRI (Yoruba people from Ibadan, Nigeria) samples; (ii) LD information, in the case of a HapMap SNP, including what genes have variation in strong LD (pairwise or multilocus LD) with the variant and how well the SNP is covered by different high-throughput platforms; (iii) summary information available from public databases (e.g. physical and functional annotations); and (iv) summary information from other GWAS. For each gene, SCAN provides annotations on: (i) eQTLs for the gene (both local and distant SNPs) and (ii) the coverage of all variants in the HapMap at that gene on each high-throughput platform. For each genomic region, SCAN provides annotations on: (i) physical and functional annotations of all SNPs, genes and known CNVs within the region and (ii) all genes regulated by the eQTLs within the region.

AVAILABILITY

http://www.scandb.org.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Genome-wide linkage scans for type 2 diabetes mellitus in four ethnically diverse populations-significant evidence for linkage on chromosome 4q in African Americans: the Family Investigation of Nephropathy and Diabetes Research Group

BACKGROUND

Previous studies have shown that in addition to environmental influences, type 2 diabetes mellitus (T2DM) has a strong genetic component. The goal of the current study is to identify regions of linkage for T2DM in ethnically diverse populations.

METHODS

Phenotypic and genotypic data were obtained from African American (AA; total number of individuals [N] = 1004), American Indian (AI; N = 883), European American (EA; N = 537), and Mexican American (MA; N = 1634) individuals from the Family Investigation of Nephropathy and Diabetes. Non-parametric linkage analysis, using an average of 4404 SNPs, was performed in relative pairs affected with T2DM in each ethnic group. In addition, family-based tests were performed to detect association with T2DM.

RESULTS

Statistically significant evidence for linkage was observed on chromosome 4q21.1 (LOD = 3.13; genome-wide p = 0.04) in AA. In addition, a total of 11 regions showed suggestive evidence for linkage (estimated at LOD > 1.71), with the highest LOD scores on chromosomes 12q21.31 (LOD = 2.02) and 22q12.3 (LOD = 2.38) in AA, 2p11.1 (LOD = 2.23) in AI, 6p12.3 (LOD = 2.77) in EA, and 13q21.1 (LOD = . 2.24) in MA. While no region overlapped across all ethnic groups, at least five loci showing LOD > 1.71 have been identified in previously published studies.

CONCLUSIONS

The results from this study provide evidence for the presence of genes affecting T2DM on chromosomes 4q, 12q, and 22q in AA; 6p in EA; 2p in AI; and 13q in MA. The strong evidence for linkage on chromosome 4q in AA provides important information given the paucity of diabetes genetic studies in this population.

Evaluation of gremlin 1 (GREM1) as a candidate susceptibility gene for albuminuria-related traits in Mexican Americans with type 2 diabetes mellitus

Several novel genes that are up-regulated in the kidney in diabetes have been identified including GREM1, which encodes gremlin 1. GREM1 maps to human chromosome 15q12, a region previously found to be linked to albumin to creatinine ratio (ACR) in Mexican Americans. The objective of this study is to investigate whether genetic variants in GREM1, a positional candidate gene, contribute to variation in ACR. By sequencing 32 individuals for both exons and 2-kilobase putative promoter region of GREM1, we identified 19 genetic variants including 5 in the promoter region and 13 in the 3' untranslated region. Of 19 polymorphisms identified, 13 polymorphisms were genotyped in the entire cohort (N = 670, 39 large families) either by restriction fragment length polymorphism or by TaqMan (Applied Biosystems, Foster City, CA) assays. Association analyses between the genotypes and ACR, type 2 diabetes mellitus, and related phenotypes were carried out using a measured genotype approach as implemented in the variance component analytical tools (SOLAR). Of the variants examined for association, none exhibited statistically significant association with ACR after accounting for the effects of covariates such as age, sex, diabetes, duration of diabetes, systolic blood pressure, and antihypertensive medications. However, 2 novel variants at the 3' untranslated region showed significant association with estimated glomerular filtration rate (P = .010 and P = .049) and body mass index (P = .013 and P = .019) after accounting for trait-specific covariate influences. Furthermore, a novel variant located in the promoter exhibited a significant association with systolic (P = .038) and diastolic blood pressure (P = .005) after adjusting for the effects of age, sex, diabetes, and antihypertensive medications. In conclusion, the variants examined at GREM1 are not significant contributors to variation in ACR in Mexican Americans, although they appear to minimally influence risk factors related to ACR.

A generalized family-based association test for dichotomous traits

Recent advances in genotyping technology make it possible to utilize large-scale association analysis for disease-gene mapping. Powerful and robust family-based association methods are crucial for successful gene mapping. We propose a family-based association method, the generalized disequilibrium test (GDT), in which the genotype differences of all discordant relative pairs are utilized in assessing association within a family. The improvement of the GDT over existing methods is threefold: (1) information beyond first-degree relatives is incorporated efficiently, yielding substantial gains in power in comparison to existing tests; (2) the GDT statistic is implemented via a robust technique that does not rely on large sample theory, resulting in further power gains, especially at high levels of significance; and (3) covariates and weights based on family size are incorporated. Advantages of the GDT over existing methods are demonstrated by extensive computer simulations and by application to recently published large-scale genome-wide linkage data from the Type 1 Diabetes Genetics Consortium (T1DGC). In our simulations, the GDT consistently outperforms other tests for a common disease and frequently outperforms other tests for a rare disease; the power improvement is > 13% in 6 out of 8 extended pedigree scenarios. All of the six strongest associations identified by the GDT have been reported by other studies, whereas only three or four of these associations can be identified by existing methods. For the T1D association at gene UBASH3A, the GDT resulted in a genome-wide significance (p = 4.3 x 10(-6)), much stronger than the published significance (p = 10(-4)).

Vitamin D receptor gene polymorphisms have negligible effect on human height

Human height is a highly heritable trait, with genetic factors explaining up to 90% of phenotypic variation. Vitamin D levels are known to influence several physiological processes, including skeletal growth. The vitamin D receptor (VDR) gene has been reported as contributing to variation in height. A meta-analysis of 13607 adult individuals found a small but significant association with the rs1544410 (BsmI) polymorphism. In contrast, the meta-analysis found no effect in a sample of 550 children. Two recent studies reported variants with large effect on height elsewhere in VDR (rs10735810 [FokI] and rs7139166 [-1,521] polymorphisms). We genotyped large Caucasian samples from Australia (N = 3,906) and the Netherlands (N = 1,689) for polymorphisms in VDR. The Australian samples were twin families with height measures from 3 time points throughout adolescence. The Dutch samples were adult twins. We use the available family data to perform both within and between family tests of association. We found no significant associations for any of the genotyped variants after multiple testing correction. The (non-significant) effect of rs1544410 in the Australian adolescent cohort was in the same direction and of similar magnitude (additive effect 0.3 cm) to the effect observed in the published adult meta-analysis. An effect of this size explains approximately 0.1% of the phenotypic variance in height - this implies that many, probably hundreds, of such variants are responsible for the observed genetic variation. Our results did not support any role for two other regions (rs10735810, rs7139166) of VDR in explaining variation in height.

Tests of association for quantitative traits in nuclear families using principal components to correct for population stratification

Traditional transmission disequilibrium test (TDT) based methods for genetic association analyses are robust to population stratification at the cost of a substantial loss of power. We here describe a novel method for family-based association studies that corrects for population stratification with the use of an extension of principal component analysis (PCA). Specifically, we adopt PCA on unrelated parents in each family. We then infer principal components for children from those for their parents through a TDT-like strategy. Two test statistics within the variance-components model are proposed for association tests. Simulation results show that the proposed tests have correct type I error rates regardless of population stratification, and have greatly improved power over two popular TDT-based methods: QTDT and FBAT. The application to the Genetic Analysis Workshop 16 (GAW16) data sets attests to the feasibility of the proposed method.

Genome-wide association study of bipolar disorder in European American and African American individuals

To identify bipolar disorder (BD) genetic susceptibility factors, we conducted two genome-wide association (GWA) studies: one involving a sample of individuals of European ancestry (EA; n=1001 cases; n=1033 controls), and one involving a sample of individuals of African ancestry (AA; n=345 cases; n=670 controls). For the EA sample, single-nucleotide polymorphisms (SNPs) with the strongest statistical evidence for association included rs5907577 in an intergenic region at Xq27.1 (P=1.6 x 10(-6)) and rs10193871 in NAP5 at 2q21.2 (P=9.8 x 10(-6)). For the AA sample, SNPs with the strongest statistical evidence for association included rs2111504 in DPY19L3 at 19q13.11 (P=1.5 x 10(-6)) and rs2769605 in NTRK2 at 9q21.33 (P=4.5 x 10(-5)). We also investigated whether we could provide support for three regions previously associated with BD, and we showed that the ANK3 region replicates in our sample, along with some support for C15Orf53; other evidence implicates BD candidate genes such as SLITRK2. We also tested the hypothesis that BD susceptibility variants exhibit genetic background-dependent effects. SNPs with the strongest statistical evidence for genetic background effects included rs11208285 in ROR1 at 1p31.3 (P=1.4 x 10(-6)), rs4657247 in RGS5 at 1q23.3 (P=4.1 x 10(-6)), and rs7078071 in BTBD16 at 10q26.13 (P=4.5 x 10(-6)). This study is the first to conduct GWA of BD in individuals of AA and suggests that genetic variations that contribute to BD may vary as a function of ancestry.

Genetics of human gene expression: mapping DNA variants that influence gene expression

There is extensive natural variation in human gene expression. As quantitative phenotypes, expression levels of genes are heritable. Genetic linkage and association mapping have identified cis- and trans-acting DNA variants that influence expression levels of human genes. New insights into human gene regulation are emerging from genetic analyses of gene expression in cells at rest and following exposure to stimuli. The integration of these genetic mapping results with data from co-expression networks is leading to a better understanding of how expression levels of individual genes are regulated and how genes interact with each other. These findings are important for basic understanding of gene regulation and of diseases that result from disruption of normal gene regulation.

The ordered transmission disequilibrium test: detection of modifier genes

We consider the problem of detection of modifier genes that lead to variations in a disease-related continuous variable (DRCV), such as the age of onset or a measure of disease severity, in a strategy of candidate genes. We propose a novel method, the ordered transmission disequilibrium test (OTDT), to test for a relation between the clinical heterogeneity expressed by a DRCV and marker genotypes of a candidate gene. The OTDT applies to trio families with one patients and his parents, all three genotyped at a bi-allelic marker M. The OTDT aims to find a critical value of the DRCV which separates the sample of families in two subsamples in which the transmission rates are significantly different. We investigate the power of the method by simulations under various genetic models and covariate distributions and compare it with a linear regression analysis.

Genetic aspects of rheumatic diseases

In this review on the genetic aspects of rheumatic diseases, the approach was taken (i) to discuss, in general, important principles in the identification of susceptibility genes and (ii) to focus on five autoimmune rheumatic diseases that have the characteristics of complex diseases and in which important advances have been made in the identification of the genetic component. A decade ago, most reviews on the genetics of rheumatic diseases focused almost exclusively on a discussion of the human leukocyte antigen (HLA) but in this fast-moving field it is now apparent that genes outside the HLA also contribute to susceptibility. Current hypotheses concerning the pathogenesis of autoimmunity have led to the inclusion of hundreds of genes as potential candidates. Almost any molecule involved in immune recognition, cell interaction, intracellular signalling, cytokine pathways or programmed cell death can be reasonably proposed. The identification of the genes involved in complex diseases will contribute to an understanding of disease mechanisms and disease biology. The disease pathways by which the genes exert their effects or functions could lead to the discovery of new therapeutic targets that may be modulated. An increased understanding of the interactions between genes and environment might also be attained.

Association of AKT1 with verbal learning, verbal memory, and regional cortical gray matter density in twins

AKT1, encoding the protein kinase B, has been associated with the genetic etiology of schizophrenia and bipolar disorder. However, minuscule data exist on the role of different alleles of AKT1 in measurable quantitative endophenotypes, such as cognitive abilities and neuroanatomical features, showing deviations in schizophrenia and bipolar disorder. We evaluated the contribution of AKT1 to quantitative cognitive traits and 3D high-resolution neuroanatomical images in a Finnish twin sample consisting of 298 twins: 61 pairs with schizophrenia (8 concordant), 31 pairs with bipolar disorder (5 concordant) and 65 control pairs matched for age, sex and demographics. An AKT1 allele defined by the SNP rs1130214 located in the UTR of the gene revealed association with cognitive traits related to verbal learning and memory (P = 0.0005 for a composite index). This association was further fortified by a higher degree of resemblance of verbal memory capacity in pairs sharing the rs1130214 genotype compared to pairs not sharing the genotype. Furthermore, the same allele was also associated with decreased gray matter density in medial and dorsolateral prefrontal cortex (P < 0.05). Our findings support the role of AKT1 in the genetic background of cognitive and anatomical features, known to be affected by psychotic disorders. The established association of the same allelic variant of AKT1 with both cognitive and neuroanatomical aberrations could suggest that AKT1 exerts its effect on verbal learning and memory via neural networks involving prefrontal cortex.

The ERAP2 gene is associated with preeclampsia in Australian and Norwegian populations

Preeclampsia is a heritable pregnancy disorder that presents new onset hypertension and proteinuria. We have previously reported genetic linkage to preeclampsia on chromosomes 2q, 5q and 13q in an Australian/New Zealand (Aust/NZ) familial cohort. This current study centered on identifying the susceptibility gene(s) at the 5q locus. We first prioritized candidate genes using a bioinformatic tool designed for this purpose. We then selected a panel of known SNPs within ten prioritized genes and genotyped them in an extended set of the Aust/NZ families and in a very large, independent Norwegian case/control cohort (1,139 cases, 2,269 controls). In the Aust/NZ cohort we identified evidence of a genetic association for the endoplasmic reticulum aminopeptidase 1 (ERAP1) gene (rs3734016, P (uncorr) = 0.009) and for the endoplasmic reticulum aminopeptidase 2 (ERAP2) gene (rs2549782, P (uncorr) = 0.004). In the Norwegian cohort we identified evidence of a genetic association for ERAP1 (rs34750, P (uncorr) = 0.011) and for ERAP2 (rs17408150, P (uncorr) = 0.009). The ERAP2 SNPs in both cohorts remained statistically significant (rs2549782, P (corr) = 0.018; rs17408150, P (corr) = 0.039) after corrections at an experiment-wide level. The ERAP1 and ERAP2 genes encode enzymes that are reported to play a role in blood pressure regulation and essential hypertension in addition to innate immune and inflammatory responses. Perturbations within vascular, immunological and inflammatory pathways constitute important physiological mechanisms in preeclampsia pathogenesis. We herein report a novel preeclampsia risk locus, ERAP2, in a region of known genetic linkage to this pregnancy-specific disorder.

Evaluation of polymorphisms in paraoxonase 2 (PON2) gene and their association with cardiovascular-renal disease risk in Mexican Americans

BACKGROUND/AIMS

Genetic polymorphisms in the paraoxonase 2 (PON2) gene are thought to alter its activity and contribute to the development of cardiovascular and renal disease risk. The purpose of this study is to determine whether the Arg148Gly, Cys311Ser and rs12794795 polymorphisms of PON2 examined previously by others, are associated with type 2 diabetes (T2DM), and subclinical measures of cardiovascular and renal disease risk in Mexican Americans.

METHODS

Study participants (n = 848; 21 families) were genotyped for the three polymorphisms by TaqMan assay. Association between the genotypic and phenotypic data was performed by measured genotype approach as implemented in the variance component analytical tools.

RESULTS

The Arg148Gly variant was found to be monomorphic in our dataset. Of the phenotypes examined for association, the A/C variant located in intron-1 (rs12794795) exhibited statistically significant association only with diastolic blood pressure (p = 0.018) after accounting for the trait-specific covariate effects. The Cys311Ser variant failed to show statistically significant association with any of the phenotypes examined.

CONCLUSION

In conclusion, the variants examined at the PON2 locus in Mexican Americans do not appear to be a major contributor to T2DM, cardiovascular or renal disease risk, although they exhibited a small effect on the blood pressure values.

Association between genes of Disrupted in schizophrenia 1 (DISC1) interactors and schizophrenia supports the role of the DISC1 pathway in the etiology of major mental illnesses

BACKGROUND

Disrupted in Schizophrenia 1 (DISC1) is currently one of the most interesting candidate genes for major mental illness, having been demonstrated to associate with schizophrenia, bipolar disorder, major depression, autism, and Asperger's syndrome. We have previously reported a DISC1 haplotype, HEP3, and an NDE1 spanning tag haplotype to associate to schizophrenia in Finnish schizophrenia families. Because both DISC1 and NDE1 display association in our study sample, we hypothesized that other genes interacting with DISC1 might also have a role in the etiology of schizophrenia.

METHODS

We selected 11 additional genes encoding components of the "DISC1 pathway" and studied these in our study sample of 476 families including 1857 genotyped individuals. We performed single nucleotide polymorphism (SNP) and haplotype association analyses in two independent sets of families. For markers and haplotypes found to be consistently associated in both sets, the overall significance was tested with the combined set of families.

RESULTS

We identified three SNPs to be associated with schizophrenia in PDE4D (rs1120303, p = .021), PDE4B (rs7412571, p = .018), and NDEL1 (rs17806986, p = .0038). Greater significance was observed with allelic haplotypes of PDE4D (p = .00084), PDE4B (p = .0022 and p = .029), and NDEL1 (p = .0027) that increased or decreased schizophrenia susceptibility.

CONCLUSIONS

Our findings with other converging lines of evidence support the underlying importance of DISC1-related molecular pathways in the etiology of schizophrenia and other major mental illnesses.

Linkage and association between CA repeat polymorphism of the TNFR2 gene and obesity phenotypes in two independent Caucasian populations

Previously, our group has reported a suggestive linkage evidence of 1p36 with body mass index (BMI) (LOD = 2.09). The tumor necrosis factor receptor 2 (TNFR2) at 1p36 is an excellent positional and functional candidate gene for obesity. In this study, we have investigated the linkage and association between the TNFR2 gene and obesity phenotypes in two large independent samples, using the quantitative transmission disequilibrium tests (QTDT). The first group was made up of 1,836 individuals from 79 multi-generation pedigrees. The second group was a randomly ascertained set of 636 individuals from 157 US Caucasian nuclear families. Obesity phenotypes tested include BMI, fat mass, and percentage fat mass (PFM). A significant result (P = 0.0056) was observed for linkage with BMI in the sample of the multigenerational pedigrees. Our data support the TNFR2 gene as a quantitative trait locus (QTL) underlying BMI variation in the Caucasian populations.

The 27-bp repeat polymorphism in intron 4 (27 bp-VNTR) of endothelial nitric oxide synthase (eNOS) gene is associated with albumin to creatinine ratio in Mexican Americans

The T-786C, Glu298Asp, and 27 bp variable number of tandem repeats (27 bp-VNTR-a/b) polymorphsims of the endothelial nitric oxide synthase (eNOS) gene are thought to alter nitric oxide production and contribute to the development of vascular and renal disease risk. The objective of this study is to investigate whether these three polymorphisms examined previously by others are associated with cardiovascular and renal disease risk in Mexican Americans. Study participants (N = 848; 21 families) were genotyped for T-786C, Glu298Asp, and 27 bp-VNTR-a/b polymorphisms by PCR followed by restriction digestion. Association analyses were performed by a measured genotype approach implemented in the program SOLAR. Of the phenotypes (type 2 diabetes, hypertension, body mass index, waist circumference, total cholesterol, high density lipoprotein cholesterol, triglycerides, systolic and diastolic blood pressure, albumin to creatinine ratio (ACR), and estimated glomerular filtration rate) examined for association, the 27 bp-VNTR-a/b variant exhibited statistically significant association with ACR (P = 0.047) after accounting for the trait specific covariate effects. In addition, the promoter variant (T-786C) showed a significant association with triglycerides (P = 0.034) after accounting for covariate influences. In conclusion, the present study adds evidence to the role of eNOS candidate gene polymorphisms in modulating the risk factors related to cardiovascular-renal disease in Mexican Americans although the magnitude of the genetic effect is small.

Genetic variation in female BMI increases with number of children born but failure to replicate association between GNbeta3 variants and increased BMI in parous females

Objective:As post-pregnancy weight retention in women is a risk factor for obesity later in life, we assessed the changes in magnitude of genetic and environmental variation in BMI due to parity in an Australian female sample, comprising 11,915 female twins and their sisters.Results:Total variance of BMI increased with parity, primarily driven by an increase in nonadditive genetic variance. This finding was of particular interest given Gutersohn et al's (2000) report of recessive association between post-partum weight retention and the 825T allele of theGNβ3gene (rs5443) at 12p13.31. Hence, we attempted to replicate this association and test an additional three SNPs also located near or onGNβ3(rs10744716 (upstream), rs5442 (exon 10), rs5446 (exon 11)) in a sample of 3131 females and 1693 males from 2585 twin families. No association was found between these SNPs among females, even when allowing for a genotype by parity effect. However, a significant association was observed among males for rs10744716 (χ22= 10.22,p= 0.006; effect size = 0.47kg/m2), representing a novel association between a region upstream ofGNβ3and male BMI.

No association between Cholinergic Muscarinic Receptor 2 (CHRM2) genetic variation and cognitive abilities in three independent samples

Cognitive ability has a substantial genetic component and more than 15 candidate genes have been identified over the past 8 years. One gene that has been associated with general cognitive ability is the cholinergic muscarinic 2 receptor (CHRM2). In an attempt to replicate this finding we typed marker rs8191992 (the originally reported CHRM2 SNP) in two population based cohorts-one Scottish aged over 50 years (N = 2,091) and the other English comprising non-demented elderly participants (N = 758)-and a family-based Australian adolescent sample (N = 1,537). An additional 29 SNPs in CHRM2 were typed in the Australian sample and a further seven in the English cohort. No significant association was found between CHRM2 and diverse measures of cognitive ability in any of the samples. In conclusion, this study does not support a role for CHRM2 in cognitive ability.

Family-based mapping of quantitative trait loci in plant breeding populations with resistance to Fusarium head blight in wheat as an illustration

Traditional quantitative trait loci (QTL) mapping approaches are typically based on early or advanced generation analysis of bi-parental populations. A limitation associated with this methodology is the fact that mapping populations rarely give rise to new cultivars. Additionally, markers linked to the QTL of interest are often not immediately available for use in breeding and they may not be useful within diverse genetic backgrounds. Use of breeding populations for simultaneous QTL mapping, marker validation, marker assisted selection (MAS), and cultivar release has recently caught the attention of plant breeders to circumvent the weaknesses of conventional QTL mapping. The first objective of this study was to test the feasibility of using family-pedigree based QTL mapping techniques generally used with humans and animals within plant breeding populations (PBPs). The second objective was to evaluate two methods (linkage and association) to detect marker-QTL associations. The techniques described in this study were applied to map the well characterized QTL, Fhb1 for Fusarium head blight resistance in wheat (Triticum aestivum L.). The experimental populations consisted of 82 families and 793 individuals. The QTL was mapped using both linkage (variance component and pedigree-wide regression) and association (using quantitative transmission disequilibrium test, QTDT) approaches developed for extended family-pedigrees. Each approach successfully identified the known QTL location with a high probability value. Markers linked to the QTL explained 40-50% of the phenotypic variation. These results show the usefulness of a human genetics approach to detect QTL in PBPs and subsequent use in MAS.

Linkage and association analyses of longitudinally measured lipid phenotypes in adolescence

The genetic basis of cardiovascular disease (CVD) is complex and still largely elusive. Plasma lipid concentrations are well-established risk factors for cardiovascular disease (CVD), and have adult heritabilities ranging from 0.48 to 0.87. Estimates for adolescents are slightly higher (range 0.71 to 0.82). To identify loci affecting lipid concentrations across adolescence, we analyzed longitudinal lipid data in a sample of 134 monozygotic and 626 dizygotic twin pairs at ages twelve, fourteen and sixteen, and their siblings, from 760 Australian families. Univariate linkage analysis for each phenotype and time point was supplemented by multivariate analysis across the time points. A genome-wide association scan was also performed on a subset of the subjects (N = 441). The strongest linkage was seen for triglycerides on chromosome 6p24.3 (multivariate -log(10) p = 6.81; equivalent LOD = 6.13; p = 1.55 x 10(-7)). Significant linkage was also found for LDL cholesterol on chromosome 2q35 (multivariate -log(10)p = 5.59; equivalent LOD = 4.53; p = 2.57 x 10(-6)). In the association analysis, rs10503840 on 8p21.1 was significantly associated with total cholesterol levels at age fourteen (p = 8.24 x 10(-7), estimated significance threshold 2.45 x 10(-6)). Association at p < 2.25 x 10(-6) was also found between triglycerides at age 12 and rs10507266, in an intron of THRAP2 (MIM 608771) on 12q24.21; and between HDL-C at age 14 and rs10506325 in an intergenic region of 12q13.13. Suggestive evidence of association at ages twelve and fourteen was found between HDL-C and rs10492859 on 16q23 (p = 2.42 x 10(-5) and 2.77 x 10(-4), respectively). Further longitudinal genetic studies of cardiovascular risk factors, focused on critical periods of development or change, are needed.

Finding genes underlying human disease

With the surge in genome-wide association studies (GWAS), many have asked the question 'Are linkage studies dead?' In this article, we survey the approaches used in mapping human disease genes, reviewing the analysis strategies that preceded and laid the groundwork for GWAS. We note that earlier approaches are still useful and the development of new methodology is warranted.

Genetic association study for RSV bronchiolitis in infancy at the 5q31 cytokine cluster

BACKGROUND

The pathophysiological basis of severe respiratory syncytial virus (RSV) bronchiolitis in infancy is poorly understood and has hindered vaccine development. Studies implicate the cell-mediated immune response in the pathogenesis of the disease. A recent twin study estimated a heritable contribution of 22% to RSV bronchiolitis. Genetic epidemiology provides a new approach to identifying important immune determinants of disease severity.

METHODS

A comprehensive high-density gene-region association study for severe RSV bronchiolitis in infancy at 5q31 across 11 genes including the Th2-cytokine cluster was performed. A haplotype tagging approach was used to analyse genetic variation at 113 single nucleotide polymorphisms (SNPs) in 780 independent cases and 1045 controls. The study had sufficient power to detect small effects, perform extensive haplotype analysis and analyse both a principal phenotype and a refined age-limited phenotype enriched for first-exposure RSV infection.

RESULTS

SNP associations were found at IL4 and a highly significant risk haplotype was identified across IL13 CNS-1 and IL4 (odds ratio 1.69, p<0.0001), present in both case-control and family-based analyses. All associations were strongest for a phenotype limited to <6 months of age, implicating this locus in primary RSV disease. The same risk haplotype has previously been shown to be associated with increased IL13 expression.

CONCLUSIONS

A haplotype at IL13-1L4, which is associated with increased IL13 production, confers an increased risk of severe primary RSV bronchiolitis in early infancy. This study, together with previous studies implicating the same locus in atopic sensitisation, suggests that primary RSV bronchiolitis and atopy share a genetic contribution at the IL13-IL4 locus.

Associations of 25 structural, degradative, and inflammatory candidate genes with lumbar disc desiccation, bulging, and height narrowing

OBJECTIVE

To examine the allelic diversity of structural, inflammatory, and matrix-modifying gene candidates and their association with disc degeneration.

METHODS

Subjects were 588 men ages 35-70 years. We investigated associations of single-nucleotide polymorphisms in AGC1 and in 12 collagen, 8 interleukin, and 4 matrix metalloproteinase genes with quantitative magnetic resonance imaging measurements of disc desiccation and disc bulging and height narrowing scores, after controlling for age and suspected risk factors. Analyses were performed using QTDT software. P values were derived from 1,000 permutations, and empirical P values for global significance also were applied.

RESULTS

Twelve of the 99 variants in 25 selected candidate genes provided evidence of association (P < 0.05) with disc signal intensity in the upper and/or lower lumbar regions. Allelic variants of AGC1 (rs1042631; P = 0.001), COL1A1 (rs2075555; P = 0.005), COL9A1 (rs696990; P = 0.00008), and COL11A2 (rs2076311; P = 0.018) genes provided the most significant evidence of association with disc signal intensity. The same variants of AGC1 (P = 0.010) and COL9A1 (P = 0.014), as well as variants in the COL11A1 gene (rs1463035 [P = 0.004]; rs1337185 [P = 0.015]) were also associated with disc bulging, as was AGC1 with disc height narrowing (rs1516797; P = 0.005). In addition, 4 allelic variants in the immunologic candidate genes (rs2071375 in IL1A [P = 0.027]; rs1420100 in IL18RAP [P = 0.005]) were associated with disc signal intensity.

CONCLUSION

Genetic variants account for interindividual differences in disc matrix synthesis and degradation. The accuracy of the quantitative disc signal intensity measurements we used likely enhanced our ability to observe these associations. Our findings shed light on possible mechanisms of degeneration and support the view that disc degeneration is a polygenetic condition.

Polymorphisms in SLC6A4, PAH, GABRB3, and MAOB and modification of psychotic disorder features

We tested four genes [phenylalanine hydroxylase (PAH), the serotonin transporter (SLC6A4), monoamine oxidase B (MAOB), and the gamma-aminobutyric acid A receptor beta-3 subunit (GABRB3)] for their impact on five schizophrenia symptom factors: delusions, hallucinations, mania, depression, and negative symptoms. In a 90 family subset of the Irish Study of High Density Schizophrenia Families, the PAH 232 bp microsatellite allele demonstrated significant association with the delusions factor using both QTDT (F=8.0, p=.031) and QPDTPHASE (chi-square=12.54, p=.028). Also, a significant association between the GABRB3 191 bp allele and the hallucinations factor was detected using QPDTPHASE (chi-square=15.51, p=.030), but not QTDT (chi-square=2.07, p=.560).

Loss-of-function variants in endothelial lipase are a cause of elevated HDL cholesterol in humans

Elevated plasma concentrations of HDL cholesterol (HDL-C) are associated with protection from atherosclerotic cardiovascular disease. Animal models indicate that decreased expression of endothelial lipase (LIPG) is inversely associated with HDL-C levels, and genome-wide association studies have identified LIPG variants as being associated with HDL-C levels in humans. We hypothesized that loss-of-function mutations in LIPG may result in elevated HDL-C and therefore performed deep resequencing of LIPG exons in cases with elevated HDL-C levels and controls with decreased HDL-C levels. We identified a significant excess of nonsynonymous LIPG variants unique to cases with elevated HDL-C. In vitro lipase activity assays demonstrated that these variants significantly decreased endothelial lipase activity. In addition, a meta-analysis across 5 cohorts demonstrated that the low-frequency Asn396Ser variant is significantly associated with increased HDL-C, while the common Thr111Ile variant is not. Functional analysis confirmed that the Asn396Ser variant has significantly decreased lipase activity both in vitro and in vivo, while the Thr111Ile variant has normal lipase activity. Our results establish that loss-of-function mutations in LIPG lead to increased HDL-C levels and support the idea that inhibition of endothelial lipase may be an effective mechanism to raise HDL-C.

Leptin is associated with blood pressure and hypertension in women from the National Heart, Lung, and Blood Institute Family Heart Study

Leptin is a key neuroendocrine hormone regulating food intake, metabolism, and fat accumulation, and it may also affect blood pressure and contribute to hypertension through sympathetic activation in the vasculature or at the renal level. Although previous studies have shown that the distribution of leptin is significantly different between males and females, as is the risk of hypertension between males and females, results regarding the role of leptin in the gender-specific regulation of blood pressure are controversial. Thus, we performed family-based association analyses in the National Heart, Lung, and Blood Institute Family Heart Study to test the hypothesis that LEPTIN gene variants and the plasma leptin level influence variability in blood pressure and the risk of hypertension differently by gender. We identified significant associations between LEPTIN single nucleotide polymorphisms with blood pressure and hypertension, but in postmenopausal women only. We also identified significant associations between plasma leptin levels and both blood pressure and hypertension in women. The current study supports a role for LEPTIN and plasma leptin levels in blood pressure regulation in women. It also provides insight into the gender differences in hypertension, as well as the differential distribution and activity of leptin in men and women.

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.

Genetic association of the activin A receptor gene (ACVR2A) and pre-eclampsia

Pre-eclampsia is a common serious disorder of human pregnancy, which is associated with significant maternal and perinatal morbidity and mortality. The suspected aetiology of pre-eclampsia is complex, with susceptibility being attributable to multiple environmental factors and a large genetic component. Recently, we reported significant linkage to chromosome 2q22 in 34 Australian/New Zealand (Aust/NZ) pre-eclampsia/eclampsia families, and activin A receptor type IIA (ACVR2A) was identified as a strong positional candidate gene at this locus. In an attempt to identify the putative risk variants, we have now comprehensively re-sequenced the entire coding region of the ACVR2A gene and the conserved non-coding sequences in a subset of 16 individuals from these families. We identified 45 single nucleotide polymorphisms (SNPs), with 9 being novel. These SNPs were genotyped in our total family sample of 480 individuals from 74 Aust/NZ pre-eclampsia families (including the original 34 genome-scanned families). Our best associations between ACVR2A polymorphisms and pre-eclampsia were for rs10497025 (P = 0.025), rs13430086 (P = 0.010) and three novel SNPs: LF004, LF013 and LF020 (all with P = 0.018). After correction for multiple hypothesis testing, none of these associations reached significance (P > 0.05). Based on these data, it remains unclear what role, if any, ACVR2A polymorphisms play in pre-eclampsia risk, at least in these Australian families. However, it would be premature to rule out this gene as significant associations between ACVR2A SNPs and pre-eclampsia have recently been reported in a large Norwegian (HUNT) population sample.

Identification of genomic regions contributing to etoposide-induced cytotoxicity

Etoposide is routinely used in combination-based chemotherapy for testicular cancer and small-cell lung cancer; however, myelosuppression, therapy-related leukemia and neurotoxicity limit its utility. To determine the genetic contribution to cellular sensitivity to etoposide, we evaluated cell growth inhibition in Centre d' Etude du Polymorphisme Humain lymphoblastoid cell lines from 24 multi-generational pedigrees (321 samples) following treatment with 0.02-2.5 microM etoposide for 72 h. Heritability analysis showed that genetic variation contributes significantly to the cytotoxic phenotypes (h (2) = 0.17-0.25, P = 4.9 x 10(-5)-7.3 x 10(-3)). Whole genome linkage scans uncovered 8 regions with peak LOD scores ranging from 1.57 to 2.55, with the most significant signals being found on chromosome 5 (LOD = 2.55) and chromosome 6 (LOD = 2.52). Linkage-directed association was performed on a subset of HapMap samples within the pedigrees to find 22 SNPs significantly associated with etoposide cytotoxicity at one or more treatment concentrations. UVRAG, a DNA repair gene, SEMA5A, SLC7A6 and PRMT7 are implicated from these unbiased studies. Our findings suggest that susceptibility to etoposide-induced cytotoxicity is heritable and using an integrated genomics approach we identified both genomic regions and SNPs associated with the cytotoxic phenotypes.

Genetic variants in COMT and neurocognitive impairment in families of patients with schizophrenia

This study examined the relations of genetic variants in catechol-O-methyltransferase (COMT) gene, including rs737865 in intron 1, rs4680 in exon 4 (Val158Met) and downstream rs165599, to schizophrenia and its related neurocognitive functions in families of patients with schizophrenia. Totally, 680 individuals from 166 simplex (166 affected members and 354 nonpsychotic first-degree relatives) and 46 multiplex families (85 affected members and 75 nonpsychotic first-degree relatives) were interviewed using Diagnostic Interview for Genetic Studies, administered Wisconsin Card Sorting Test (WCST) and Continuous Performance Test (CPT), and drawn for venous blood. Both categorical (dichotomizing families on affected members' neurocognitive performance) and quantitative approaches toward the WCST and CPT performance scores were employed using the family-based association test and the variance components framework, respectively. Both false discovery rate and permutations were used to adjust for multiple testing. The genotypes of rs4680 were associated with both the WCST and CPT performance scores in these families, but not with schizophrenia per se in either whole sample or subgroup analyses. Meanwhile, the other two single nucleotide polymorphisms were differentially associated with the two tasks. For WCST indexes, regardless of subgroup analyses or quantitative approach, only rs737865 exhibited moderate associations. For CPT indexes, rs737865 exhibited association for the subgroup with deficit on CPT reaction time, whereas rs165599 exhibited association for the subgroup with deficit on CPT d' as well as quantitative undegraded d'. Our results indicate that the genetic variants in COMT might be involved in modulation of neurocognitive functions and hence conferring increased risk to schizophrenia.