We conducted a genome-wide association study for nonsyndromic cleft lip with or without cleft palate (NSCL/P) in 401 affected individuals and 1,323 controls, with replication in an independent sample of 793 NSCL/P triads. We report two new loci associated with NSCL/P at 17q22 (rs227731, combined P = 1.07 × 10 −8 , relative risk in homozygotes = 1.84, 95% CI 1.34-2.53) and 10q25.3 (rs7078160, combined P = 1.92 × 10 −8 , relative risk in homozygotes = 2.17, 95% CI 1.32-3.56).NSCL/P is one of the most common human birth defects. In European populations, NSCL/P has a prevalence ranging from 1 in 700 to 1 in 1,000. We recently reported a susceptibility locus for NSCL/P at chromo some 8q24.21 from a genome wide association study in 224 individuals with NSCL/P (cases) and 383 population based controls 1 . This locus is the second susceptibility locus to have been unequivocally identified for NSCL/P to date, the first being the IRF6 locus 2 .To identify additional cleft susceptibility loci, we enlarged our sample by genotyping an additional set of 177 NSCL/P cases and adding the genotypes of 940 population based controls of central European origin. Genotyping was performed using Illumina BeadChips (Human610 Quad and HumanHap 550k).Following quality control (Supplementary Methods and Supplementary Fig. 1), association analysis of 521,288 SNPs having a minor allele frequency (MAF) of ≥1% in controls was performed in 399 cases and 1,318 controls.After excluding markers from the previously described 8q24.21 locus, 20 SNPs with P < 10 −5 remained. Five chromosomal loci (8q12.3, 10q25.3, 13q31.1, 15q13.3 and 17q22) were located within these 20 top SNPs, and the associations at these loci were further supported by at least three more SNPs with P < 10 −4 ( Supplementary Fig. 2 and Supplementary Table 1). Two additional regions were considered to be promising NSCL/P susceptibility loci (6p22.1, 11q14.2), as they contained at least four markers with P < 10 −4 .To replicate the genome wide association study (GWAS) findings, we selected the 20 top SNPs (P < 10 −5 ) as well as additional backup markers for each of the seven previously mentioned loci, resulting in two replication assays. We included additional SNPs with P < 10 −4 in the two replication assays, giving highest priority to SNPs with the lowest P values. Thus, a total of 56 markers were genotyped in a replication sample of 793 NSCL/P triads of European origin. Genotyping using matrix assisted laser desorption/ionization time of flight (MALDI TOF) mass spectrometry (Sequenom Inc.) was successful for 45 markers (representing 32 different loci), which were then analyzed by the transmission disequilibrium test in 665 triads (128 triads were excluded after quality control, Supplementary Methods).Of the 45 SNPs successfully genotyped, 11 (representing six differ ent loci) showed P < 0.05 in the replication sample (Supplementary Table 2). Two of these SNPs remained significant after correction for multiple testing by a conservative Bonferroni procedure (17q22: rs227731, P corr ...

We conducted a genome-wide association study involving 224 cases and 383 controls of Central European origin to identify susceptibility loci for nonsyndromic cleft lip with or without cleft palate (NSCL/P). A 640-kb region at chromosome 8q24.21 was found to contain multiple markers with highly significant evidence for association with the cleft phenotype, including three markers that reached genome-wide significance. The 640-kb cleft-associated region was saturated with 146 SNP markers and then analyzed in our entire NSCL/P sample of 462 unrelated cases and 954 controls. In the entire sample, the most significant SNP (rs987525) had a P value of 3.34 x 10(-24). The odds ratio was 2.57 (95% CI = 2.02-3.26) for the heterozygous genotype and 6.05 (95% CI = 3.88-9.43) for the homozygous genotype. The calculated population attributable risk for this marker is 0.41, suggesting that this study has identified a major susceptibility locus for NSCL/P.

We have conducted the first meta-analyses for nonsyndromic cleft lip with or without cleft palate (NSCL/P) using data from the two largest genome-wide association studies published to date. We confirmed associations with all previously identified loci and identified six additional susceptibility regions (1p36, 2p21, 3p11.1, 8q21.3, 13q31.1 and 15q22). Analysis of phenotypic variability identified the first specific genetic risk factor for NSCLP (nonsyndromic cleft lip plus palate) (rs8001641; PNSCLP = 6.51 × 10−11; homozygote relative risk = 2.41, 95% confidence interval (CI) 1.84–3.16).

12The cerebral cortex underlies our complex cognitive capabilities, yet we know little about the specific genetic loci influencing human cortical structure. To identify genetic variants, including structural variants, impacting cortical structure, we conducted a genome-wide association meta-analysis of brain MRI data from 51,662 individuals. We analysed the surface area and average thickness of the whole cortex and 34 regions with known functional specialisations. We identified 255 nominally significant loci (P ≤ 5 x 10 -8 ); 199 survived multiple testing correction (P ≤ 8.3 x 10 -10 ; 187 surface area; 12 thickness). We found significant enrichment for loci influencing total surface area within regulatory elements active during prenatal cortical development, supporting the radial unit hypothesis. Loci impacting regional surface area cluster near genes in Wnt signalling pathways, known to influence progenitor expansion and areal identity. Variation in cortical structure is genetically correlated with cognitive function, Parkinson's disease, insomnia, depression and ADHD.One Sentence Summary: Common genetic variation is associated with inter-individual variation in the structure of the human cortex, both globally and within specific regions, and is shared with genetic risk factors for some neuropsychiatric disorders.The human cerebral cortex is the outer grey matter layer of the brain, which is implicated in multiple aspects of higher cognitive function. Its distinct folding pattern is characterised by convex (gyral) and concave (sulcal) regions. Computational brain mapping approaches use the consistent folding patterns across individual cortices to label brain regions(1). During fetal development excitatory neurons, the predominant neuronal cell-type in the cortex, are generated from neural progenitor cells in the developing germinal zone(2). The radial unit hypothesis(3) posits that the expansion of cortical surface area (SA) is driven by the proliferation of these neural progenitor cells, whereas thickness (TH) is determined by the number of neurogenic divisions. Variation in global and regional measures of cortical SA and TH are associated with neuropsychiatric disorders and psychological traits(4) ( Table S1). Twin and family-based brain imaging studies show that SA and TH measurements are highly heritable and are largely influenced by independent genetic factors(5). Despite extensive studies of genes impacting cortical structure in model organisms (6), our current understanding of genetic variation impacting human cortical size and patterning is limited to rare, highly penetrant variants (7,8). These variants often disrupt cortical development, leading to altered post-natal structure. However, little is known about how common genetic variants impact human cortical SA and TH.To address this, we conducted genome-wide association meta-analyses of cortical SA and TH measures in 51,662 individuals from 60 cohorts from around the world (Tables S2-S4). Cortical measures were extracted from structural brain MRI scan...

Genetic variants in a gene on 6p22.3, dysbindin, have been shown recently to be associated with schizophrenia (Straub et al. 2002a). There is no doubt that replication in other independent samples would enhance the significance of this finding considerably. Since the gene is located in the center of the linkage peak on chromosome 6p that we reported earlier, we decided to test six of the most positive DNA polymorphisms in a sib-pair sample and in an independently ascertained sample of triads comprising 203 families, including the families for which we detected linkage on chromosome 6p. Evidence for association was observed in the two samples separately as well as in the combined sample (P=.00068 for SNP rs760761). Multilocus haplotype analysis increased the significance further to .00002 for a two-locus haplotype and to .00001 for a three-locus haplotype. Estimation of frequencies for six-locus haplotypes revealed one common haplotype with a frequency of 73.4% in transmitted, and only 57.6% in nontransmitted, parental haplotypes. All other six-locus haplotypes occurring at a frequency of >1% were less often transmitted than nontransmitted. Our results represent a first successful replication of linkage disequilibrium in psychiatric genetics detected in a region with previous evidence of linkage and will encourage the search for causes of schizophrenia by the genetic approach.

The cerebral cortex underlies our complex cognitive capabilities, yet little is known about the specific genetic loci that influence human cortical structure. To identify genetic variants that affect cortical structure, we conducted a genome-wide association meta-analysis of brain magnetic resonance imaging data from 51,665 individuals. We analyzed the surface area and average thickness of the whole cortex and 34 regions with known functional specializations. We identified 199 significant loci and found significant enrichment for loci influencing total surface area within regulatory elements that are active during prenatal cortical development, supporting the radial unit hypothesis. Loci that affect regional surface area cluster near genes in Wnt signaling pathways, which influence progenitor expansion and areal identity. Variation in cortical structure is genetically correlated with cognitive function, Parkinson’s disease, insomnia, depression, neuroticism, and attention deficit hyperactivity disorder.

The transmission/disequilibrium test (TDT) is a popular method for detection of the genetic basis of a disease. Investigators planning such studies require computation of sample size and power, allowing for a general genetic model. Here, a rigorous method is presented for obtaining the power approximations of the TDT for samples consisting of families with either a single affected child or affected sib pairs. Power calculations based on simulation show that these approximations are quite precise. By this method, it is also shown that a previously published power approximation of the TDT is erroneous.

The importance of haplotype analysis in the context of association fine mapping of disease genes has grown steadily over the last years. Since experimental methods to determine haplotypes on a large scale are not available, phase has to be inferred statistically. For individual genotype data, several reconstruction techniques and many implementations of the expectation-maximization (EM) algorithm for haplotype frequency estimation exist. Recent research work has shown that incorporating available genotype information of related individuals largely increases the precision of haplotype frequency estimates. We, therefore, implemented a highly flexible program written in C, called FAMHAP, which calculates maximum likelihood estimates (MLEs) of haplotype frequencies from general nuclear families with an arbitrary number of children via the EM-algorithm for up to 20 SNPs. For more loci, we have implemented a locus-iterative mode of the EM-algorithm, which gives reliable approximations of the MLEs for up to 63 SNP loci, or less when multi-allelic markers are incorporated into the analysis. Missing genotypes can be handled as well. The program is able to distinguish cases (haplotypes transmitted to the first affected child of a family) from pseudo-controls (non-transmitted haplotypes with respect to the child). We tested the performance of FAMHAP and the accuracy of the obtained haplotype frequencies on a variety of simulated data sets. The implementation proved to work well when many markers were considered and no significant differences between the estimates obtained with the usual EM-algorithm and those obtained in its locus-iterative mode were observed. We conclude from the simulations that the accuracy of haplotype frequency estimation and reconstruction in nuclear families is very reliable in general and robust against missing genotypes.