Search

Based on epidemiologic and embryologic patterns, nonsyndromic orofacial clefts- the most common craniofacial birth defects in humans- are commonly categorized into cleft lip with or without cleft palate (CL/P) and cleft palate alone (CP), which are traditionally considered to be etiologically distinct. However, some evidence of shared genetic risk in IRF6, GRHL3 and ARHGAP29 regions exists; only FOXE1 has been recognized as significantly associated with both CL/P and CP in genome-wide association studies (GWAS). We used a new statistical approach, PLACO (pleiotropic analysis under composite null), on a combined multi-ethnic GWAS of 2,771 CL/P and 611 CP case-parent trios. At the genome-wide significance threshold of 5 × 10-8, PLACO identified 1 locus in 1q32.2 (IRF6) that appears to increase risk for one OFC subgroup but decrease risk for the other. At a suggestive significance threshold of 10-6, we found 5 more loci with compelling candidate genes having opposite effects on CL/P and CP: 1p36.13 (PAX7), 3q29 (DLG1), 4p13 (LIMCH1), 4q21.1 (SHROOM3) and 17q22 (NOG). Additionally, we replicated the recognized shared locus 9q22.33 (FOXE1), and identified 2 loci in 19p13.12 (RAB8A) and 20q12 (MAFB) that appear to influence risk of both CL/P and CP in the same direction. We found locus-specific effects may vary by racial/ethnic group at these regions of genetic overlap, and failed to find evidence of sex-specific differences. We confirmed shared etiology of the two OFC subtypes comprising CL/P, and additionally found suggestive evidence of differences in their pathogenesis at 2 loci of genetic overlap. Our novel findings include 6 new loci of genetic overlap between CL/P and CP; 3 new loci between pairwise OFC subtypes; and 4 loci not previously implicated in OFCs. Our in-silico validation showed PLACO is robust to subtype-specific effects, and can achieve massive power gains over existing approaches for identifying genetic overlap between disease subtypes. In summary, we found suggestive evidence for new genetic regions and confirmed some recognized OFC genes either exerting shared risk or with opposite effects on risk to OFC subtypes., Competing Interests: The authors have declared that no competing interests exist.

Background: Many maternal traits are associated with a neonate's gestational duration, birth weight, and birth length. These birth outcomes are subsequently associated with late-onset health conditions. The causal mechanisms and the relative contributions of maternal and fetal genetic effects behind these observed associations are unresolved., Methods and Findings: Based on 10,734 mother-infant duos of European ancestry from the UK, Northern Europe, Australia, and North America, we constructed haplotype genetic scores using single-nucleotide polymorphisms (SNPs) known to be associated with adult height, body mass index (BMI), blood pressure (BP), fasting plasma glucose (FPG), and type 2 diabetes (T2D). Using these scores as genetic instruments, we estimated the maternal and fetal genetic effects underlying the observed associations between maternal phenotypes and pregnancy outcomes. We also used infant-specific birth weight genetic scores as instrument and examined the effects of fetal growth on pregnancy outcomes, maternal BP, and glucose levels during pregnancy. The maternal nontransmitted haplotype score for height was significantly associated with gestational duration (p = 2.2 × 10-4). Both maternal and paternal transmitted height haplotype scores were highly significantly associated with birth weight and length (p < 1 × 10-17). The maternal transmitted BMI scores were associated with birth weight with a significant maternal effect (p = 1.6 × 10-4). Both maternal and paternal transmitted BP scores were negatively associated with birth weight with a significant fetal effect (p = 9.4 × 10-3), whereas BP alleles were significantly associated with gestational duration and preterm birth through maternal effects (p = 3.3 × 10-2 and p = 4.5 × 10-3, respectively). The nontransmitted haplotype score for FPG was strongly associated with birth weight (p = 4.7 × 10-6); however, the glucose-increasing alleles in the fetus were associated with reduced birth weight through a fetal effect (p = 2.2 × 10-3). The haplotype scores for T2D were associated with birth weight in a similar way but with a weaker maternal effect (p = 6.4 × 10-3) and a stronger fetal effect (p = 1.3 × 10-5). The paternal transmitted birth weight score was significantly associated with reduced gestational duration (p = 1.8 × 10-4) and increased maternal systolic BP during pregnancy (p = 2.2 × 10-2). The major limitations of the study include missing and heterogenous phenotype data in some data sets and different instrumental strength of genetic scores for different phenotypic traits., Conclusions: We found that both maternal height and fetal growth are important factors in shaping the duration of gestation: genetically elevated maternal height is associated with longer gestational duration, whereas alleles that increase fetal growth are associated with shorter gestational duration. Fetal growth is influenced by both maternal and fetal effects and can reciprocally influence maternal phenotypes: taller maternal stature, higher maternal BMI, and higher maternal blood glucose are associated with larger birth size through maternal effects; in the fetus, the height- and metabolic-risk-increasing alleles are associated with increased and decreased birth size, respectively; alleles raising birth weight in the fetus are associated with shorter gestational duration and higher maternal BP. These maternal and fetal genetic effects may explain the observed associations between the studied maternal phenotypes and birth outcomes, as well as the life-course associations between these birth outcomes and adult phenotypes., Competing Interests: I have read the journal's policy and the authors of this manuscript have the following competing interests: DAL has received support, via her University, from Roche Diagnostics and Medtronic Ltd for biomarker research unrelated to this study. LJM consults for Mirvie, Inc., a biotech company involved in preterm birth diagnostics. The content of this paper is unrelated to the areas of their diagnostics development. GDS is a member of the Editorial Board of PLOS Medicine. The other authors report no conflicts.

[This corrects the article DOI: 10.1371/journal.pgen.1007394.].

Preterm birth is a leading cause of morbidity and mortality in infants. Genetic and environmental factors play a role in the susceptibility to preterm birth, but despite many investigations, the genetic basis for preterm birth remain largely unknown. Our objective was to identify rare, possibly damaging, nucleotide variants in mothers from families with recurrent spontaneous preterm births (SPTB). DNA samples from 17 Finnish mothers who delivered at least one infant preterm were subjected to whole exome sequencing. All mothers were of northern Finnish origin and were from seven multiplex families. Additional replication samples of European origin consisted of 93 Danish sister pairs (and two sister triads), all with a history of a preterm delivery. Rare exonic variants (frequency <1%) were analyzed to identify genes and pathways likely to affect SPTB susceptibility. We identified rare, possibly damaging, variants in genes that were common to multiple affected individuals. The glucocorticoid receptor signaling pathway was the most significant (p<1.7e-8) with genes containing these variants in a subgroup of ten Finnish mothers, each having had 2-4 SPTBs. This pathway was replicated among the Danish sister pairs. A gene in this pathway, heat shock protein family A (Hsp70) member 1 like (HSPA1L), contains two likely damaging missense alleles that were found in four different Finnish families. One of the variants (rs34620296) had a higher frequency in cases compared to controls (0.0025 vs. 0.0010, p = 0.002) in a large preterm birth genome-wide association study (GWAS) consisting of mothers of general European ancestry. Sister pairs in replication samples also shared rare, likely damaging HSPA1L variants. Furthermore, in silico analysis predicted an additional phosphorylation site generated by rs34620296 that could potentially affect chaperone activity or HSPA1L protein stability. Finally, in vitro functional experiment showed a link between HSPA1L activity and decidualization. In conclusion, rare, likely damaging, variants in HSPA1L were observed in multiple families with recurrent SPTB., Competing Interests: The 23andMe authors are employees of and hold stock or stock options in 23andMe, Inc.

Cleft palate only is a common birth defect with high heritability. Only a small fraction of this heritability is explained by the genetic variants identified so far, underscoring the need to investigate other disease mechanisms, such as gene-environment (GxE) interactions and parent-of-origin (PoO) effects. Furthermore, PoO effects may vary across exposure levels (PoOxE effects). Such variation is the focus of this study. We upgraded the R-package Haplin to enable direct tests of PoOxE effects at the genome-wide level. From a previous GWAS, we had genotypes for 550 case-parent trios, of mainly European and Asian ancestry, and data on three maternal exposures (smoking, alcohol, and vitamins). Data were analyzed for Europeans and Asians separately, and also for all ethnicities combined. To account for multiple testing, a false discovery rate method was used, where q-values were generated from the p-values. In the Europeans-only analyses, interactions with maternal smoking yielded the lowest q-values. Two SNPs in the 'Interactor of little elongation complex ELL subunit 1' (ICE1) gene had a q-value of 0.14, and five of the 20 most significant SNPs were in the 'N-acetylated alpha-linked acidic dipeptidase-like 2' (NAALADL2) gene. No evidence of PoOxE effects was found in the other analyses. The connections to ICE1 and NAALADL2 are novel and warrant further investigation. More generally, the new methodology presented here is easily applicable to other traits and exposures in which a family-based study design has been implemented.

Numerous lines of evidence point to a genetic basis for facial morphology in humans, yet little is known about how specific genetic variants relate to the phenotypic expression of many common facial features. We conducted genome-wide association meta-analyses of 20 quantitative facial measurements derived from the 3D surface images of 3118 healthy individuals of European ancestry belonging to two US cohorts. Analyses were performed on just under one million genotyped SNPs (Illumina OmniExpress+Exome v1.2 array) imputed to the 1000 Genomes reference panel (Phase 3). We observed genome-wide significant associations (p < 5 x 10-8) for cranial base width at 14q21.1 and 20q12, intercanthal width at 1p13.3 and Xq13.2, nasal width at 20p11.22, nasal ala length at 14q11.2, and upper facial depth at 11q22.1. Several genes in the associated regions are known to play roles in craniofacial development or in syndromes affecting the face: MAFB, PAX9, MIPOL1, ALX3, HDAC8, and PAX1. We also tested genotype-phenotype associations reported in two previous genome-wide studies and found evidence of replication for nasal ala length and SNPs in CACNA2D3 and PRDM16. These results provide further evidence that common variants in regions harboring genes of known craniofacial function contribute to normal variation in human facial features. Improved understanding of the genes associated with facial morphology in healthy individuals can provide insights into the pathways and mechanisms controlling normal and abnormal facial morphogenesis., Competing Interests: The authors have declared that no competing interests exist.

[This corrects the article DOI: 10.1371/journal.pgen.1005378.].

Background: The aim of this study was to investigate the risk of psychiatric diagnoses in individuals with non-syndromic oral clefts (OC) compared with individuals without OC, including ages from 1 to 76 years., Methods: Linking four Danish nationwide registers, we investigated the risk of psychiatric diagnoses at Danish psychiatric hospitals during the period 1969-2012 for individuals born with non-syndromic OC in Denmark 1936-2009 compared with a cohort of 10 individuals without OC per individual with OC, matched by sex and birth year. The sample included 8,568 individuals with OC, observed for 247,821 person-years, and 85,653 individuals without OC followed for 2,501,129 person-years., Results: A total of 953 (11.1%) of the individuals with OC (9.6% for cleft lip (CL), 10.8% for cleft lip and palate (CLP) and 13.1% for cleft palate (CP)) and 8,117 (9.5%) in the comparison group had at least one psychiatric diagnosis. Cox proportional hazard regression model revealed that individuals with OC had significantly higher risk of a psychiatric diagnosis (hazard ratio (HR) = 1.19, 95% CI: 1.12-1.28). When examining cleft type, no difference was found for CL (HR = 1.03, 95% CI: 0.90-1.17), but CLP was associated with a small increased risk (HR = 1.13, 95% CI: 1.01-1.26), whereas individuals with CP had the largest increased risk (HR = 1.45, 95% CI: 1.30-1.62). The largest differences were found in schizophrenia-like disorders, mental retardation and pervasive developmental disorders, but we found no increased risk of mood disorders and anxiety-related disorders., Conclusion: Individuals with non-syndromic OC had significantly higher risk of psychiatric diagnoses compared with individuals without OC. However, the elevated risk was observed for individuals with CLP and CP but not for individuals with CL and the absolute risk increase was modest.

Much of our current knowledge regarding the association of FOXP2 with speech and language development comes from singleton and small family studies where a small number of rare variants have been identified. However, neither genome-wide nor gene-specific studies have provided evidence that common polymorphisms in the gene contribute to individual differences in language development in the general population. One explanation for this inconsistency is that previous studies have been limited to relatively small samples of individuals with low language abilities, using low density gene coverage. The current study examined the association between common variants in FOXP2 and a quantitative measure of language ability in a population-based cohort of European decent (n = 812). No significant associations were found for a panel of 13 SNPs that covered the coding region of FOXP2 and extended into the promoter region. Power analyses indicated we should have been able to detect a QTL variance of 0.02 for an associated allele with MAF of 0.2 or greater with 80% power. This suggests that, if a common variant associated with language ability in this gene does exist, it is likely of small effect. Our findings lead us to conclude that while genetic variants in FOXP2 may be significant for rare forms of language impairment, they do not contribute appreciably to individual variation in the normal range as found in the general population.

Genome-wide association studies (GWAS) have identified more than 100 genetic variants contributing to BMI, a measure of body size, or waist-to-hip ratio (adjusted for BMI, WHRadjBMI), a measure of body shape. Body size and shape change as people grow older and these changes differ substantially between men and women. To systematically screen for age- and/or sex-specific effects of genetic variants on BMI and WHRadjBMI, we performed meta-analyses of 114 studies (up to 320,485 individuals of European descent) with genome-wide chip and/or Metabochip data by the Genetic Investigation of Anthropometric Traits (GIANT) Consortium. Each study tested the association of up to ~2.8M SNPs with BMI and WHRadjBMI in four strata (men ≤50y, men >50y, women ≤50y, women >50y) and summary statistics were combined in stratum-specific meta-analyses. We then screened for variants that showed age-specific effects (G x AGE), sex-specific effects (G x SEX) or age-specific effects that differed between men and women (G x AGE x SEX). For BMI, we identified 15 loci (11 previously established for main effects, four novel) that showed significant (FDR<5%) age-specific effects, of which 11 had larger effects in younger (<50y) than in older adults (≥50y). No sex-dependent effects were identified for BMI. For WHRadjBMI, we identified 44 loci (27 previously established for main effects, 17 novel) with sex-specific effects, of which 28 showed larger effects in women than in men, five showed larger effects in men than in women, and 11 showed opposite effects between sexes. No age-dependent effects were identified for WHRadjBMI. This is the first genome-wide interaction meta-analysis to report convincing evidence of age-dependent genetic effects on BMI. In addition, we confirm the sex-specificity of genetic effects on WHRadjBMI. These results may provide further insights into the biology that underlies weight change with age or the sexually dimorphism of body shape.

Cleft lip with or without cleft palate (CL/P) is the most commonly occurring craniofacial birth defect. We provide insight into the genetic etiology of this birth defect by performing genome-wide association studies in two species: dogs and humans. In the dog, a genome-wide association study of 7 CL/P cases and 112 controls from the Nova Scotia Duck Tolling Retriever (NSDTR) breed identified a significantly associated region on canine chromosome 27 (unadjusted p=1.1 x 10(-13); adjusted p= 2.2 x 10(-3)). Further analysis in NSDTR families and additional full sibling cases identified a 1.44 Mb homozygous haplotype (chromosome 27: 9.29 - 10.73 Mb) segregating with a more complex phenotype of cleft lip, cleft palate, and syndactyly (CLPS) in 13 cases. Whole-genome sequencing of 3 CLPS cases and 4 controls at 15X coverage led to the discovery of a frameshift mutation within ADAMTS20 (c.1360&#95;1361delAA (p.Lys453Ilefs*3)), which segregated concordant with the phenotype. In a parallel study in humans, a family-based association analysis (DFAM) of 125 CL/P cases, 420 unaffected relatives, and 392 controls from a Guatemalan cohort, identified a suggestive association (rs10785430; p =2.67 x 10-6) with the same gene, ADAMTS20. Sequencing of cases from the Guatemalan cohort was unable to identify a causative mutation within the coding region of ADAMTS20, but four coding variants were found in additional cases of CL/P. In summary, this study provides genetic evidence for a role of ADAMTS20 in CL/P development in dogs and as a candidate gene for CL/P development in humans.

Cleft palate (CP) is one of the most commonly occurring craniofacial birth defects in humans. In order to study cleft palate in a naturally occurring model system, we utilized the Nova Scotia Duck Tolling Retriever (NSDTR) dog breed. Micro-computed tomography analysis of CP NSDTR craniofacial structures revealed that these dogs exhibit defects similar to those observed in a recognizable subgroup of humans with CP: Pierre Robin Sequence (PRS). We refer to this phenotype in NSDTRs as CP1. Individuals with PRS have a triad of birth defects: shortened mandible, posteriorly placed tongue, and cleft palate. A genome-wide association study in 14 CP NSDTRs and 72 unaffected NSDTRs identified a significantly associated region on canine chromosome 14 (24.2 Mb-29.3 Mb; p(raw )= 4.64 × 10(-15)). Sequencing of two regional candidate homeobox genes in NSDTRs, distal-less homeobox 5 (DLX5) and distal-less homeobox 6 (DLX6), identified a 2.1 kb LINE-1 insertion within DLX6 in CP1 NSDTRs. The LINE-1 insertion is predicted to insert a premature stop codon within the homeodomain of DLX6. This prompted the sequencing of DLX5 and DLX6 in a human cohort with CP, where a missense mutation within the highly conserved DLX5 homeobox of a patient with PRS was identified. This suggests the involvement of DLX5 in the development of PRS. These results demonstrate the power of the canine animal model as a genetically tractable approach to understanding naturally occurring craniofacial birth defects in humans.

Nonsyndromic cleft palate (CP) is one of the most common human birth defects and both genetic and environmental risk factors contribute to its etiology. We conducted a genome-wide association study (GWAS) using 550 CP case-parent trios ascertained in an international consortium. Stratified analysis among trios with different ancestries was performed to test for GxE interactions with common maternal exposures using conditional logistic regression models. While no single nucleotide polymorphism (SNP) achieved genome-wide significance when considered alone, markers in SLC2A9 and the neighboring WDR1 on chromosome 4p16.1 gave suggestive evidence of gene-environment interaction with environmental tobacco smoke (ETS) among 259 Asian trios when the models included a term for GxE interaction. Multiple SNPs in these two genes were associated with increased risk of nonsyndromic CP if the mother was exposed to ETS during the peri-conceptual period (3 months prior to conception through the first trimester). When maternal ETS was considered, fifteen of 135 SNPs mapping to SLC2A9 and 9 of 59 SNPs in WDR1 gave P values approaching genome-wide significance (10(-6)

We investigated the housekeeping cytochrome P450 CYP51A1 encoding lanosterol 14α-demethylase from cholesterol synthesis that was so far not directly linked to human disorders. By direct sequencing of CYP51A1 in 188 women with spontaneous preterm delivery and 188 unrelated preterm infants (gestational age <37 weeks) we identified 22 variants where 10 are novel and rare. In infants there were two novel CYP51A1 variants where damaging effects of p.Tyr145Asp from the substrate recognition region, but not p.Asn193Asp, were predicted by PolyPhen2 and SIFT. This was confirmed by molecular modeling showing that Tyr145Asp substitution results in changed electrostatic potential of the CYP51 protein surface and lengthened distance to the heme which prevents hydrogen bonding. The CYP51 Tyr145Asp mutation is rare and thus very interesting for further structure/function relationship studies. From the 12 identified known variants rs6465348 was chosen for family based association studies due to its high minor allele frequency. Interestingly, this CYP51A1 common variant associates with small for gestational age weight in newborns (p = 0.028) and lower blood total cholesterol and low density lipoprotein cholesterol levels in mothers in 2nd trimester of pregnancy (p = 0.042 and p = 0.046 respectively). Our results indicate a new link between a cholesterol synthesis gene CYP51A1 and pregnancy pathologies.

Orofacial clefts are among the most common birth defects and result in an improper formation of the mouth or the roof of the mouth. Monosomy of the distal aspect of human chromosome 6p has been recognized as causative in congenital malformations affecting the brain and cranial skeleton including orofacial clefts. Among the genes located in this region is PAK1IP1, which encodes a nucleolar factor involved in ribosomal stress response. Here, we report the identification of a novel mouse line that carries a point mutation in the Pak1ip1 gene. Homozygous mutants show severe developmental defects of the brain and craniofacial skeleton, including a median orofacial cleft. We recovered this line of mice in a forward genetic screen and named the allele manta-ray (mray). Our findings prompted us to examine human cases of orofacial clefting for mutations in the PAK1IP1 gene or association with the locus. No deleterious variants in the PAK1IP1 gene coding region were recognized, however, we identified a borderline association effect for SNP rs494723 suggesting a possible role for the PAK1IP1 gene in human orofacial clefting.

Background: Recent epidemiological studies suggest that the maternal genome is an important contributor to spontaneous preterm delivery (PTD). There is also a significant excess of males among preterm born infants, which may imply an X-linked mode of inheritance for a subset of cases. To explore this, we examined the effect of maternal and fetal X-chromosomal single nucleotide polymorphisms (SNPs) on the risk of PTD in two independent genome-wide association studies and one replication study., Methods: Participants were recruited from the Danish National Birth Cohort and the Norwegian Mother and Child cohort studies. Data from these two populations were first analyzed independently, and then combined in a meta-analysis. Overall, we evaluated 12,211 SNPs in 1,535 case-mother dyads and 1,487 control-mother dyads. Analyses were done using a hybrid design that combines case-mother dyads and control-mother dyads, as implemented in the Haplin statistical software package. A sex-stratified analysis was performed for the fetal SNPs. In the replication study, 10 maternal and 16 fetal SNPs were analyzed using case-parent triads from independent studies of PTD in the United States, Argentina and Denmark., Results: In the meta-analysis, the G allele at the maternal SNP rs2747022 in the FERM domain containing 7 gene (FRMD7) increased the risk of spontaneous PTD by 1.2 (95% confidence interval (CI): 1.1, 1.4). Although an association with this SNP was confirmed in the replication study, it was no longer statistically significant after a Bonferroni correction for multiple testing., Conclusion: We did not find strong evidence in our data to implicate X-chromosomal SNPs in the etiology of spontaneous PTD. Although non-significant after correction for multiple testing, the mother's G allele at rs2747022 in FRMD7 increased the risk of spontaneous PTD across all populations in this study, thus warranting further investigation in other populations.

Previous evidence from tooth agenesis studies suggested IRF6 and TGFA interact. Since tooth agenesis is commonly found in individuals with cleft lip/palate (CL/P), we used four large cohorts to evaluate if IRF6 and TGFA interaction contributes to CL/P. Markers within and flanking IRF6 and TGFA genes were tested using Taqman or SYBR green chemistries for case-control analyses in 1,000 Brazilian individuals. We looked for evidence of gene-gene interaction between IRF6 and TGFA by testing if markers associated with CL/P were overtransmitted together in the case-control Brazilian dataset and in the additional family datasets. Genotypes for an additional 142 case-parent trios from South America drawn from the Latin American Collaborative Study of Congenital Malformations (ECLAMC), 154 cases from Latvia, and 8,717 individuals from several cohorts were available for replication of tests for interaction. Tgfa and Irf6 expression at critical stages during palatogenesis was analyzed in wild type and Irf6 knockout mice. Markers in and near IRF6 and TGFA were associated with CL/P in the Brazilian cohort (p<10(-6)). IRF6 was also associated with cleft palate (CP) with impaction of permanent teeth (p<10(-6)). Statistical evidence of interaction between IRF6 and TGFA was found in all data sets (p = 0.013 for Brazilians; p = 0.046 for ECLAMC; p = 10(-6) for Latvians, and p = 0.003 for the 8,717 individuals). Tgfa was not expressed in the palatal tissues of Irf6 knockout mice. IRF6 and TGFA contribute to subsets of CL/P with specific dental anomalies. Moreover, this potential IRF6-TGFA interaction may account for as much as 1% to 10% of CL/P cases. The Irf6-knockout model further supports the evidence of IRF6-TGFA interaction found in humans.

Background: Orofacial clefts are common birth defects of complex etiology, with an excess of males among babies with cleft lip and palate, and an excess of females among those with cleft palate only. Although genes on the X chromosome have been implicated in clefting, there has been no association analysis of X-linked markers., Methodology/principal Findings: We added new functionalities in the HAPLIN statistical software to enable association analysis of X-linked markers and an exploration of various causal scenarios relevant to orofacial clefts. Genotypes for 48 SNPs in 18 candidate genes on the X chromosome were analyzed in two population-based samples from Scandinavia (562 Norwegian and 235 Danish case-parent triads). For haplotype analysis, we used a sliding-window approach and assessed isolated cleft lip with or without cleft palate (iCL/P) separately from isolated cleft palate only (iCPO). We tested three statistical models in HAPLIN, allowing for: i) the same relative risk in males and females, ii) sex-specific relative risks, and iii) X-inactivation in females. We found weak but consistent associations with the oral-facial-digital syndrome 1 (OFD1) gene (formerly known as CXORF5) in the Danish iCL/P samples across all models, but not in the Norwegian iCL/P samples. In sex-specific analyses, the association with OFD1 was in male cases only. No analyses showed associations with iCPO in either the Norwegian or the Danish sample., Conclusions: The association of OFD1 with iCL/P is plausible given the biological relevance of this gene. However, the lack of replication in the Norwegian samples highlights the need to verify these preliminary findings in other large datasets. More generally, the novel analytic methods presented here are widely applicable to investigations of the role of X-linked genes in complex traits.

The sequence and timing of permanent tooth eruption is thought to be highly heritable and can have important implications for the risk of malocclusion, crowding, and periodontal disease. We conducted a genome-wide association study of number of permanent teeth erupted between age 6 and 14 years, analyzed as age-adjusted standard deviation score averaged over multiple time points, based on childhood records for 5,104 women from the Danish National Birth Cohort. Four loci showed association at P<5×10(-8) and were replicated in four independent study groups from the United States and Denmark with a total of 3,762 individuals; all combined P-values were below 10(-11). Two loci agreed with previous findings in primary tooth eruption and were also known to influence height and breast cancer, respectively. The two other loci pointed to genomic regions without any previous significant genome-wide association study results. The intronic SNP rs7924176 in ADK could be linked to gene expression in monocytes. The combined effect of the four genetic variants was most pronounced between age 10 and 12 years, where children with 6 to 8 delayed tooth eruption alleles had on average 3.5 (95% confidence interval: 2.9-4.1) fewer permanent teeth than children with 0 or 1 of these alleles., Competing Interests: TW declares having received consultancy and lecture fees from Lundbeck A/S.

Preterm birth is the major cause of neonatal death and serious morbidity. Most preterm births are due to spontaneous onset of labor without a known cause or effective prevention. Both maternal and fetal genomes influence the predisposition to spontaneous preterm birth (SPTB), but the susceptibility loci remain to be defined. We utilized a combination of unique population structures, family-based linkage analysis, and subsequent case-control association to identify a susceptibility haplotype for SPTB. Clinically well-characterized SPTB families from northern Finland, a subisolate founded by a relatively small founder population that has subsequently experienced a number of bottlenecks, were selected for the initial discovery sample. Genome-wide linkage analysis using a high-density single-nucleotide polymorphism (SNP) array in seven large northern Finnish non-consanginous families identified a locus on 15q26.3 (HLOD 4.68). This region contains the IGF1R gene, which encodes the type 1 insulin-like growth factor receptor IGF-1R. Haplotype segregation analysis revealed that a 55 kb 12-SNP core segment within the IGF1R gene was shared identical-by-state (IBS) in five families. A follow-up case-control study in an independent sample representing the more general Finnish population showed an association of a 6-SNP IGF1R haplotype with SPTB in the fetuses, providing further evidence for IGF1R as a SPTB predisposition gene (frequency in cases versus controls 0.11 versus 0.05, P = 0.001, odds ratio 2.3). This study demonstrates the identification of a predisposing, low-frequency haplotype in a multifactorial trait using a well-characterized population and a combination of family and case-control designs. Our findings support the identification of the novel susceptibility gene IGF1R for predisposition by the fetal genome to being born preterm., Competing Interests: The authors have declared that no competing interests exist.

Background: Fetal conditions can in principle be affected by the mother's genotype working through the prenatal environment., Methodology/principal Findings: Genotypes for 1536 SNPs in 357 cleft candidate genes were available from a previous analysis in which we focused on fetal gene effects. After data-cleaning, genotypes for 1315 SNPs in 334 autosomal genes were available for the current analysis of maternal gene effects. Two complementary statistical methods, TRIMM and HAPLIN, were used to detect multi-marker effects in population-based samples from Norway (562 case-parent and 592 control-parent triads) and Denmark (235 case-parent triads). We analyzed isolated cleft lip with or without cleft palate (iCL/P) and isolated cleft palate only (iCP) separately and assessed replication by looking for genes detected in both populations by both methods. In iCL/P, neither TRIMM nor HAPLIN detected more genes than expected by chance alone; furthermore, the selected genes were not replicated across the two methods. In iCP, however, FLNB was identified by both methods in both populations. Although HIC1 and ZNF189 did not fully satisfy our stringency criterion for replication, they were strongly associated with iCP in TRIMM analyses of the Norwegian triads., Conclusion/significance: Except for FLNB, HIC1 and ZNF189, maternal genes did not appear to influence the risk of clefting in our data. This is consistent with recent epidemiological findings showing no apparent difference between mother-to-offspring and father-to-offspring recurrence of clefts in these two populations. It is likely that fetal genes make the major genetic contribution to clefting risk in these populations, but we cannot rule out the possibility that maternal genes can affect risk through interactions with specific teratogens or fetal genes.

Background: Facial clefts are common birth defects with a strong genetic component. To identify fetal genetic risk factors for clefting, 1536 SNPs in 357 candidate genes were genotyped in two population-based samples from Scandinavia (Norway: 562 case-parent and 592 control-parent triads; Denmark: 235 case-parent triads)., Methodology/principal Findings: We used two complementary statistical methods, TRIMM and HAPLIN, to look for associations across these two national samples. TRIMM tests for association in each gene by using multi-SNP genotypes from case-parent triads directly without the need to infer haplotypes. HAPLIN on the other hand estimates the full haplotype distribution over a set of SNPs and estimates relative risks associated with each haplotype. For isolated cleft lip with or without cleft palate (I-CL/P), TRIMM and HAPLIN both identified significant associations with IRF6 and ADH1C in both populations, but only HAPLIN found an association with FGF12. For isolated cleft palate (I-CP), TRIMM found associations with ALX3, MKX, and PDGFC in both populations, but only the association with PDGFC was identified by HAPLIN. In addition, HAPLIN identified an association with ETV5 that was not detected by TRIMM., Conclusion/significance: Strong associations with seven genes were replicated in the Scandinavian samples and our approach effectively replicated the strongest previously known association in clefting--with IRF6. Based on two national cleft cohorts of similar ancestry, two robust statistical methods and a large panel of SNPs in the most promising cleft candidate genes to date, this study identified a previously unknown association with clefting for ADH1C and provides additional candidates and analytic approaches to advance the field.

Nonsyndromic or isolated cleft lip with or without cleft palate (CL/P) occurs in wide geographic distribution with an average birth prevalence of 1/700. We used direct sequencing as an approach to study candidate genes for CL/P. We report here the results of sequencing on 20 candidate genes for clefts in 184 cases with CL/P selected with an emphasis on severity and positive family history. Genes were selected based on expression patterns, animal models, and/or role in known human clefting syndromes. For seven genes with identified coding mutations that are potentially etiologic, we performed linkage disequilibrium studies as well in 501 family triads (affected child/mother/father). The recently reported MSX1 P147Q mutation was also studied in an additional 1,098 cleft cases. Selected missense mutations were screened in 1,064 controls from unrelated individuals on the Centre d'Etude du Polymorphisme Humain (CEPH) diversity cell line panel. Our aggregate data suggest that point mutations in these candidate genes are likely to contribute to 6% of isolated clefts, particularly those with more severe phenotypes (bilateral cleft of the lip with cleft palate). Additional cases, possibly due to microdeletions or isodisomy, were also detected and may contribute to clefts as well. Sequence analysis alone suggests that point mutations in FOXE1, GLI2, JAG2, LHX8, MSX1, MSX2, SATB2, SKI, SPRY2, and TBX10 may be rare causes of isolated cleft lip with or without cleft palate, and the linkage disequilibrium data support a larger, as yet unspecified, role for variants in or near MSX2, JAG2, and SKI. This study also illustrates the need to test large numbers of controls to distinguish rare polymorphic variants and prioritize functional studies for rare point mutations., Competing Interests: Competing interests. The authors have declared that no competing interests exist.