Functional gene group analysis reveals a role of synaptic heterotrimeric G proteins in cognitive ability.

Although cognitive ability is a highly heritable complex trait, only a few genes have been identified, explaining relatively low proportions of the observed trait variation. This implies that hundreds of genes of small effect may be of importance for cognitive ability. We applied an innovative method in which we tested for the effect of groups of genes defined according to cellular function (functional gene group analysis). Using an initial sample of 627 subjects, this functional gene group analysis detected that synaptic heterotrimeric guanine nucleotide binding proteins (G proteins) play an important role in cognitive ability (P(EMP) = 1.9 x 10(-4)). The association with heterotrimeric G proteins was validated in an independent population sample of 1507 subjects. Heterotrimeric G proteins are central relay factors between the activation of plasma membrane receptors by extracellular ligands and the cellular responses that these induce, and they can be considered a point of convergence, or a "signaling bottleneck." Although alterations in synaptic signaling processes may not be the exclusive explanation for the association of heterotrimeric G proteins with cognitive ability, such alterations may prominently affect the properties of neuronal networks in the brain in such a manner that impaired cognitive ability and lower intelligence are observed. The reported association of synaptic heterotrimeric G proteins with cognitive ability clearly points to a new direction in the study of the genetic basis of cognitive ability.

A variant in LIN28B is associated with 2D:4D finger-length ratio, a putative retrospective biomarker of prenatal testosterone exposure.

The ratio of the lengths of an individual's second to fourth digit (2D:4D) is commonly used as a noninvasive retrospective biomarker for prenatal androgen exposure. In order to identify the genetic determinants of 2D:4D, we applied a genome-wide association approach to 1507 11-year-old children from the Avon Longitudinal Study of Parents and Children (ALSPAC) in whom 2D:4D ratio had been measured, as well as a sample of 1382 12- to 16-year-olds from the Brisbane Adolescent Twin Study. A meta-analysis of the two scans identified a single variant in the LIN28B gene that was strongly associated with 2D:4D (rs314277: p = 4.1 x 10(-8)) and was subsequently independently replicated in an additional 3659 children from the ALSPAC cohort (p = 1.53 x 10(-6)). The minor allele of the rs314277 variant has previously been linked to increased height and delayed age at menarche, but in our study it was associated with increased 2D:4D in the direction opposite to that of previous reports on the correlation between 2D:4D and age at menarche. Our findings call into question the validity of 2D:4D as a simplistic retrospective biomarker for prenatal testosterone exposure.

Genome-wide association study reveals multiple loci associated with primary tooth development during infancy.

Tooth development is a highly heritable process which relates to other growth and developmental processes, and which interacts with the development of the entire craniofacial complex. Abnormalities of tooth development are common, with tooth agenesis being the most common developmental anomaly in humans. We performed a genome-wide association study of time to first tooth eruption and number of teeth at one year in 4,564 individuals from the 1966 Northern Finland Birth Cohort (NFBC1966) and 1,518 individuals from the Avon Longitudinal Study of Parents and Children (ALSPAC). We identified 5 loci at P<5x10(-8), and 5 with suggestive association (P<5x10(-6)). The loci included several genes with links to tooth and other organ development (KCNJ2, EDA, HOXB2, RAD51L1, IGF2BP1, HMGA2, MSRB3). Genes at four of the identified loci are implicated in the development of cancer. A variant within the HOXB gene cluster associated with occlusion defects requiring orthodontic treatment by age 31 years.

A genome-wide association study reveals variants in ARL15 that influence adiponectin levels.

The adipocyte-derived protein adiponectin is highly heritable and inversely associated with risk of type 2 diabetes mellitus (T2D) and coronary heart disease (CHD). We meta-analyzed 3 genome-wide association studies for circulating adiponectin levels (n = 8,531) and sought validation of the lead single nucleotide polymorphisms (SNPs) in 5 additional cohorts (n = 6,202). Five SNPs were genome-wide significant in their relationship with adiponectin (P< or =5x10(-8)). We then tested whether these 5 SNPs were associated with risk of T2D and CHD using a Bonferroni-corrected threshold of P< or =0.011 to declare statistical significance for these disease associations. SNPs at the adiponectin-encoding ADIPOQ locus demonstrated the strongest associations with adiponectin levels (P-combined = 9.2x10(-19) for lead SNP, rs266717, n = 14,733). A novel variant in the ARL15 (ADP-ribosylation factor-like 15) gene was associated with lower circulating levels of adiponectin (rs4311394-G, P-combined = 2.9x10(-8), n = 14,733). This same risk allele at ARL15 was also associated with a higher risk of CHD (odds ratio [OR] = 1.12, P = 8.5x10(-6), n = 22,421) more nominally, an increased risk of T2D (OR = 1.11, P = 3.2x10(-3), n = 10,128), and several metabolic traits. Expression studies in humans indicated that ARL15 is well-expressed in skeletal muscle. These findings identify a novel protein, ARL15, which influences circulating adiponectin levels and may impact upon CHD risk.

Common variants in the region around Osterix are associated with bone mineral density and growth in childhood.

Peak bone mass achieved in adolescence is a determinant of bone mass in later life. In order to identify genetic variants affecting bone mineral density (BMD), we performed a genome-wide association study of BMD and related traits in 1518 children from the Avon Longitudinal Study of Parents and Children (ALSPAC). We compared results with a scan of 134 adults with high or low hip BMD. We identified associations with BMD in an area of chromosome 12 containing the Osterix (SP7) locus, a transcription factor responsible for regulating osteoblast differentiation (ALSPAC: P = 5.8 x 10(-4); Australia: P = 3.7 x 10(-4)). This region has previously shown evidence of association with adult hip and lumbar spine BMD in an Icelandic population, as well as nominal association in a UK population. A meta-analysis of these existing studies revealed strong association between SNPs in the Osterix region and adult lumbar spine BMD (P = 9.9 x 10(-11)). In light of these findings, we genotyped a further 3692 individuals from ALSPAC who had whole body BMD and confirmed the association in children as well (P = 5.4 x 10(-5)). Moreover, all SNPs were related to height in ALSPAC children, but not weight or body mass index, and when height was included as a covariate in the regression equation, the association with total body BMD was attenuated. We conclude that genetic variants in the region of Osterix are associated with BMD in children and adults probably through primary effects on growth.

Strong evidence that GNB1L is associated with schizophrenia.

Evidence that a gene or genes on chromosome 22 is involved in susceptibility to schizophrenia comes from two sources: the increased incidence of schizophrenia in individuals with 22q11 deletion syndrome (22q11DS) and genetic linkage studies. In mice, hemizygous deletion of either Tbx1 or Gnb1l can cause deficits in pre-pulse inhibition, a sensory motor gating defect which is associated with schizophrenia. We tested the hypothesis that variation at this locus confers risk of schizophrenia and related disorders in a series of case-control association studies. First, we found evidence for a male-specific genotypic association (P = 0.00017) TBX1/GNB1L in 662 schizophrenia cases and 1416 controls from the UK. Moreover, we replicated this finding in two independent case-control samples (additional 746 cases and 1330 controls) (meta analysis P = 1.8 x 10(-5)) and also observed significant evidence for genotypic association in an independent sample of 480 schizophrenia parent-proband trios from Bulgaria with markers at this locus, which was again strongest in the male probands (P = 0.004). Genotyping the most significant SNPs in a sample of 83 subjects with 22q11DS with and without psychosis again revealed a significant allelic association with psychosis in males with 22q11DS (P = 0.01). Finally, using allele specific expression analysis, we have shown that the markers associated with psychosis are also correlated with alterations in GNB1L expression, raising the hypothesis that the risk to develop psychosis at this locus could be mediated in a dose sensitive manner via gene expression. However, other explanations are possible, and further analyses will be required to clarify the correct functional mechanism.

Analysis of copy number variation using quantitative interspecies competitive PCR.

Over recent years small submicroscopic DNA copy-number variants (CNVs) have been highlighted as an important source of variation in the human genome, human phenotypic diversity and disease susceptibility. Consequently, there is a pressing need for the development of methods that allow the efficient, accurate and cheap measurement of genomic copy number polymorphisms in clinical cohorts. We have developed a simple competitive PCR based method to determine DNA copy number which uses the entire genome of a single chimpanzee as a competitor thus eliminating the requirement for competitive sequences to be synthesized for each assay. This results in the requirement for only a single reference sample for all assays and dramatically increases the potential for large numbers of loci to be analysed in multiplex. In this study we establish proof of concept by accurately detecting previously characterized mutations at the PARK2 locus and then demonstrating the potential of quantitative interspecies competitive PCR (qicPCR) to accurately genotype CNVs in association studies by analysing chromosome 22q11 deletions in a sample of previously characterized patients and normal controls.

Comparison of methods for combining case-control and family-based association studies.

OBJECTIVES: Combining the analysis of family-based samples with unrelated individuals can enhance the power of genetic association studies. Various combined analysis techniques have been recently developed; as yet, there have been no comparisons of their power, or robustness to confounding factors. We investigated empirically the power of up to six combined methods using simulated samples of trios and unrelated cases/controls (TDTCC), trios and unrelated controls (TDTC), and affected sibpairs with parents and unrelated cases/controls (ASPFCC). METHODS: We simulated multiplicative, dominant and recessive models with varying risk parameters in single samples. Additionally, we studied false-positive rates and investigated, if possible, the coverage of the true genetic effect (TDTCC). RESULTS/CONCLUSIONS: Under the TDTCC design, we identified four approaches with equivalent power and false-positive rates. Combined statistics were more powerful than single-sample statistics or a pooled chi(2)-statistic when risk parameters were similar in single samples. Adding parental information to the CC part of the joint likelihood increased the power of generalised logistic regression under the TDTC but not the TDTCC scenario. Formal testing of differences between risk parameters in subsamples was the most sensitive approach to avoid confounding in combined analysis. Non-parametric analysis based on Monte-Carlo testing showed the highest power for ASPFCC samples.

Analysis of ProDH, COMT and ZDHHC8 risk variants does not support individual or interactive effects on schizophrenia susceptibility.

Synergistic interaction between genes on chromosome 22q11 recently has been proposed as a possible mechanism which could confer increased risk for schizophrenia. Based on this hypothesis, our study aimed to explore main, cis- and trans-interacting effects of three candidate genes on 22q11, ProDH, COMT and ZDHHC8. We selected four putative risk variants, residing within these genes, ProDH 1945, ProDH 2026, COMT ValMet and ZDHHC8 rs175174, and studied these in a large family-based schizophrenia association sample of European origin (488 Bulgarian parent-offspring trios). The presence of interaction between the variants was tested by conditional logistic regression analysis based on a case-pseudocontrol design. Our study did not find statistical evidence for allelic (investigation of ProDH markers only), genotypic, haplotypic, or interactive effects between ProDH, COMT and ZDHHC8. Our data do not support the hypothesis that an interaction between these genes influences susceptibility to schizophrenia.

No association between the putative functional ZDHHC8 single nucleotide polymorphism rs175174 and schizophrenia in large European samples.

BACKGROUND: It has been recently reported that a functional variant in the ZDHHC8 gene encoding a putative palmitoyltransferase directly confers susceptibility to schizophrenia in females (). METHODS: We investigated the putative risk allele (rs175174) in four schizophrenia association samples including a Bulgarian proband and parent sample (474 trios) and three case-control panels of European origin (1028 patients/1253 control subjects) in an attempt to replicate these findings. RESULTS: Our results do not support the hypothesis that genetic variation at rs175174 is associated with increased risk for schizophrenia nor do they suggest the presence of gender-specific differences. CONCLUSIONS: Our data suggest that the reported genetic association by either represents type I error resulting from sampling variance or that rs175174 is in linkage disequilibrium (LD) with the functional variant for schizophrenia and different LD patterns obscure the detection of association.

No association between schizophrenia and polymorphisms in COMT in two large samples.

OBJECTIVE: A valine/methionine polymorphism in the catechol O-methyltransferase (COMT) gene has been proposed to influence susceptibility to schizophrenia, as has a COMT haplotype in Ashkenazi Jewish and Irish subjects. The authors examined these hypotheses. METHOD: They reviewed data from more than 2,800 individuals, including almost 1,200 with schizophrenia, from case-control and family-based European association samples. RESULTS: The authors found no support for the hypothesis that a valine/methionine polymorphism in the COMT gene influences susceptibility to schizophrenia or the hypothesis that a COMT haplotype influences susceptibility to schizophrenia in Ashkenazi Jewish and Irish subjects. CONCLUSIONS: The data suggest that the valine allele of COMT does not increase susceptibility to schizophrenia in Europeans and that the Ashkenazi or Irish haplotype does not increase susceptibility. Ethnic variation in the linkage disequilibrium structure at COMT means that the haplotype data may not generalize across populations. However, the authors' examination of the hypothesis that the valine allele confers susceptibility, with a particularly strong effect in Europeans, reveals that no such caveat applies.

Localization of bipolar susceptibility locus by molecular genetic analysis of the chromosome 12q23-q24 region in two pedigrees with bipolar disorder and Darier's disease.

OBJECTIVE: The authors previously reported two families (pedigrees 324 and 5501) in which Darier's disease-a rare, autosomal dominant skin disease-and bipolar disorder cosegregate. In each of these families there is complete cosegregation of mood disorder with a segment of chromosome 12q23-q24, consistent with the existence of a highly penetrant dominant variant. Here molecular genetic analyses aimed at localizing and identifying the susceptibility gene in this region are reported. METHOD: In the two families, the authors undertook 1) linkage and haplotype studies using 45 highly polymorphic molecular genetic markers in order to delineate the region of interest and 2) direct analysis of genes within this region. RESULTS: Linkage and haplotype information from the most severely affected individuals defined a region of interest that spanned two neighboring regions of 19 megabases (Mb) (D12S362-D12S1646) and 7 Mb (D12S1718-D12S837). Information from all individuals refined the region of interest to 6.5 Mb (D12S127-D12S1646). Systematic study of the coding and flanking intronic regions of 25 known genes within this latter region failed to identify any highly penetrant autosomal dominant disease-conferring mutations in these pedigrees. CONCLUSIONS: This linkage and haplotype analysis, together with data from several other linkage studies, provides compelling evidence for the existence in the 12q23-q24 region of one or more genes involved in the pathogenesis of bipolar disorder. Further molecular genetic analysis of this region is required to identify the gene(s).

The moderating role of close friends in the relationship between conduct problems and adolescent substance use.

PURPOSE: Conduct problems and peer effects are among the strongest risk factors for adolescent substance use and problem use. However, it is unclear to what extent the effects of conduct problems and peer behavior interact, and whether adolescents' capacity to refuse the offer of substances may moderate such links. This study was conducted to examine relationships between conduct problems, close friends' substance use, and refusal assertiveness with adolescents' alcohol use problems, tobacco, and marijuana use. METHODS: We studied a population-based sample of 1,237 individuals from the Cardiff Study of All Wales and North West of England Twins aged 11-18 years. Adolescent and mother-reported information was obtained. Statistical analyses included cross-sectional and prospective logistic regression models and family-based permutations. RESULTS: Conduct problems and close friends' substance use were associated with increased adolescents' substance use, whereas refusal assertiveness was associated with lower use of cigarettes, alcohol, and marijuana. Peer substance use moderated the relationship between conduct problems and alcohol use problems, such that conduct problems were only related to increased risk for alcohol use problems in the presence of substance-using friends. This effect was found in both cross-sectional and prospective analyses and confirmed using the permutation approach. CONCLUSIONS: Reduced opportunities for interaction with alcohol-using peers may lower the risk of alcohol use problems in adolescents with conduct problems.

Clear detection of ADIPOQ locus as the major gene for plasma adiponectin: results of genome-wide association analyses including 4659 European individuals.

OBJECTIVE: Plasma adiponectin is strongly associated with various components of metabolic syndrome, type 2 diabetes and cardiovascular outcomes. Concentrations are highly heritable and differ between men and women. We therefore aimed to investigate the genetics of plasma adiponectin in men and women. METHODS: We combined genome-wide association scans of three population-based studies including 4659 persons. For the replication stage in 13795 subjects, we selected the 20 top signals of the combined analysis, as well as the 10 top signals with p-values less than 1.0 x 10(-4) for each the men- and the women-specific analyses. We further selected 73 SNPs that were consistently associated with metabolic syndrome parameters in previous genome-wide association studies to check for their association with plasma adiponectin. RESULTS: The ADIPOQ locus showed genome-wide significant p-values in the combined (p=4.3 x 10(-24)) as well as in both women- and men-specific analyses (p=8.7 x 10(-17) and p=2.5 x 10(-11), respectively). None of the other 39 top signal SNPs showed evidence for association in the replication analysis. None of 73 SNPs from metabolic syndrome loci exhibited association with plasma adiponectin (p>0.01). CONCLUSIONS: We demonstrated the ADIPOQ gene as the only major gene for plasma adiponectin, which explains 6.7% of the phenotypic variance. We further found that neither this gene nor any of the metabolic syndrome loci explained the sex differences observed for plasma adiponectin. Larger studies are needed to identify more moderate genetic determinants of plasma adiponectin.

Variants in ADCY5 and near CCNL1 are associated with fetal growth and birth weight.

To identify genetic variants associated with birth weight, we meta-analyzed six genome-wide association (GWA) studies (n = 10,623 Europeans from pregnancy/birth cohorts) and followed up two lead signals in 13 replication studies (n = 27,591). rs900400 near LEKR1 and CCNL1 (P = 2 x 10(-35)) and rs9883204 in ADCY5 (P = 7 x 10(-15)) were robustly associated with birth weight. Correlated SNPs in ADCY5 were recently implicated in regulation of glucose levels and susceptibility to type 2 diabetes, providing evidence that the well-described association between lower birth weight and subsequent type 2 diabetes has a genetic component, distinct from the proposed role of programming by maternal nutrition. Using data from both SNPs, we found that the 9% of Europeans carrying four birth weight-lowering alleles were, on average, 113 g (95% CI 89-137 g) lighter at birth than the 24% with zero or one alleles (P(trend) = 7 x 10(-30)). The impact on birth weight is similar to that of a mother smoking 4-5 cigarettes per day in the third trimester of pregnancy.

Visuospatial working memory in children and adolescents with 22q11.2 deletion syndrome; an fMRI study.

22q11.2 deletion syndrome (22q11DS) is a genetic disorder associated with a microdeletion of chromosome 22q11. In addition to high rates of neuropsychiatric disorders such as schizophrenia and attention deficit hyperactivity disorder, children with 22q11DS have a specific neuropsychological profile with particular deficits in visuospatial and working memory. However, the neurobiological substrate underlying these deficits is poorly understood. We investigated brain function during a visuospatial working memory (SWM) task in eight children with 22q11DS and 13 healthy controls, using fMRI. Both groups showed task-related activation in dorsolateral prefrontal cortex (DLPFC) and bilateral parietal association cortices. Controls activated parietal and occipital regions significantly more than those with 22q11DS but there was no significant between-group difference in DLPFC. In addition, while controls had a significant age-related increase in the activation of posterior brain regions and an age-related decrease in anterior regions, the 22q11DS children showed the opposite pattern. Genetically determined differences in the development of specific brain systems may underpin the cognitive deficits in 22q11DS, and may contribute to the later development of neuropsychiatric disorders.

Genetic variation in LIN28B is associated with the timing of puberty.

The timing of puberty is highly variable. We carried out a genome-wide association study for age at menarche in 4,714 women and report an association in LIN28B on chromosome 6 (rs314276, minor allele frequency (MAF) = 0.33, P = 1.5 × 10(-8)). In independent replication studies in 16,373 women, each major allele was associated with 0.12 years earlier menarche (95% CI = 0.08-0.16; P = 2.8 × 10(-10); combined P = 3.6 × 10(-16)). This allele was also associated with earlier breast development in girls (P = 0.001; N = 4,271); earlier voice breaking (P = 0.006, N = 1,026) and more advanced pubic hair development in boys (P = 0.01; N = 4,588); a faster tempo of height growth in girls (P = 0.00008; N = 4,271) and boys (P = 0.03; N = 4,588); and shorter adult height in women (P = 3.6 × 10(-7); N = 17,274) and men (P = 0.006; N = 9,840) in keeping with earlier growth cessation. These studies identify variation in LIN28B, a potent and specific regulator of microRNA processing, as the first genetic determinant regulating the timing of human pubertal growth and development.

Large-scale linkage analysis of 1302 affected relative pairs with rheumatoid arthritis.

Rheumatoid arthritis is the most common systematic autoimmune disease and its etiology is believed to have both strong genetic and environmental components. We demonstrate the utility of including genetic and clinical phenotypes as covariates within a linkage analysis framework to search for rheumatoid arthritis susceptibility loci. The raw genotypes of 1302 affected relative pairs were combined from four large family-based samples (North American Rheumatoid Arthritis Consortium, United Kingdom, European Consortium on Rheumatoid Arthritis Families, and Canada). The familiality of the clinical phenotypes was assessed. The affected relative pairs were subjected to autosomal multipoint affected relative-pair linkage analysis. Covariates were included in the linkage analysis to take account of heterogeneity within the sample. Evidence of familiality was observed with age at onset (p << 0.001) and rheumatoid factor (RF) IgM (p << 0.001), but not definite erosions (p = 0.21). Genome-wide significant evidence for linkage was observed on chromosome 6. Genome-wide suggestive evidence for linkage was observed on chromosomes 13 and 20 when conditioning on age at onset, chromosome 15 conditional on gender, and chromosome 19 conditional on RF IgM after allowing for multiple testing of covariates.

Combining linkage data sets for meta-analysis and mega-analysis: the GAW15 rheumatoid arthritis data set.

We have used the genome-wide marker genotypes from Genetic Analysis Workshop 15 Problem 2 to explore joint evidence for genetic linkage to rheumatoid arthritis across several samples. The data consisted of four high-density genome scans on samples selected for rheumatoid arthritis. We cleaned the data, removed intermarker linkage disequilibrium, and assembled the samples onto a common genetic map using genome sequence positions as a reference for map interpolation. The individual studies were combined first at the genotype level (mega-analysis) prior to a multipoint linkage analysis on the combined sample, and second using the genome scan meta-analysis method after linkage analysis of each sample. The two approaches were compared, and give strong support to the HLA locus on chromosome 6 as a susceptibility locus. Other regions of interest include loci on chromosomes 11, 2, and 12.

Analyses of single marker and pairwise effects of candidate loci for rheumatoid arthritis using logistic regression and random forests.

Using parametric and nonparametric techniques, our study investigated the presence of single locus and pairwise effects between 20 markers of the Genetic Analysis Workshop 15 (GAW15) North American Rheumatoid Arthritis Consortium (NARAC) candidate gene data set (Problem 2), analyzing 463 independent patients and 855 controls. Specifically, our work examined the correspondence between logistic regression (LR) analysis of single-locus and pairwise interaction effects, and random forest (RF) single and joint importance measures. For this comparison, we selected small but stable RFs (500 trees), which showed strong correlations (r~0.98) between their importance measures and those by RFs grown on 5000 trees. Both RF importance measures captured most of the LR single-locus and pairwise interaction effects, while joint importance measures also corresponded to full LR models containing main and interaction effects. We furthermore showed that RF measures were particularly sensitive to data imputation. The most consistent pairwise effect on rheumatoid arthritis was found between two markers within MAP3K7IP2/SUMO4 on 6q25.1, although LR and RFs assigned different significance levels.Within a hypothetical two-stage design, pairwise LR analysis of all markers with significant RF single importance would have reduced the number of possible combinations in our small data set by 61%, whereas joint importance measures would have been less efficient for marker pair reduction. This suggests that RF single importance measures, which are able to detect a wide range of interaction effects and are computationally very efficient, might be exploited as pre-screening tool for larger association studies. Follow-up analysis, such as by LR, is required since RFs do not indicate high-risk genotype combinations.