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.

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.

Mutations in interferon regulatory factor 6 (IRF6) account for ∼70% of cases of Van der Woude syndrome (VWS), the most common syndromic form of cleft lip and palate. In 8 of 45 VWS-affected families lacking a mutation in IRF6, we found coding mutations in grainyhead-like 3 (GRHL3). According to a zebrafish-based assay, the disease-associated GRHL3 mutations abrogated periderm development and were consistent with a dominant-negative effect, in contrast to haploinsufficiency seen in most VWS cases caused by IRF6 mutations. In mouse, all embryos lacking Grhl3 exhibited abnormal oral periderm and 17% developed a cleft palate. Analysis of the oral phenotype of double heterozygote (Irf6(+/-);Grhl3(+/-)) murine embryos failed to detect epistasis between the two genes, suggesting that they function in separate but convergent pathways during palatogenesis. Taken together, our data demonstrated that mutations in two genes, IRF6 and GRHL3, can lead to nearly identical phenotypes of orofacial cleft. They supported the hypotheses that both genes are essential for the presence of a functional oral periderm and that failure of this process contributes to VWS.

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.

The FaceBase Consortium consists of ten interlinked research and technology projects whose goal is to generate craniofacial research data and technology for use by the research community through a central data management and integrated bioinformatics hub. Funded by the National Institute of Dental and Craniofacial Research (NIDCR) and currently focused on studying the development of the middle region of the face, the Consortium will produce comprehensive datasets of global gene expression patterns, regulatory elements and sequencing; will generate anatomical and molecular atlases; will provide human normative facial data and other phenotypes; conduct follow up studies of a completed genome-wide association study; generate independent data on the genetics of craniofacial development, build repositories of animal models and of human samples and data for community access and analysis; and will develop software tools and animal models for analyzing and functionally testing and integrating these data. The FaceBase website (http://www.facebase.org) will serve as a web home for these efforts, providing interactive tools for exploring these datasets, together with discussion forums and other services to support and foster collaboration within the craniofacial research community.

Cleft lip and/or palate (CL/P) increase mortality and morbidity risks for affected infants especially in less developed countries. This study aimed at assessing the effects of systematic pediatric care on neonatal mortality and hospitalizations of infants with cleft lip and/or palate (CL/P) in South America.

The objective of this study was to investigate whether skewed X chromosome inactivation (XCI) is associated with IVF. A retrospective cohort study was performed comprising 30 female infants conceived by IVF and 44 naturally conceived control infants matched for gestational age and sex. Cord blood DNA samples were obtained and XCI patterns were analysed using a methylation-sensitive assay. Eight IVF samples and 13 control samples were excluded from the study because they were either homozygous or alleles were too similar for the assay to determine skewing. Mildly skewed XCI (80-90% inactivation of one allele) was present in two of 22 (9.1%) IVF samples and two of 31 (6.5%) control samples. Extremely skewed XCI (>90% inactivation of one allele) was found in two of 22 (9.1%) IVF samples and none of 31 control samples. Neither difference was statistically significant. However, the mean degree of skewed XCI in the IVF group was 72.0% and in the control group was 62.4% (P=0.002). Larger studies are needed to clarify the relationship between IVF and skewed XCI.

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.

Orofacial clefts (OFCs) are congenital malformations of the face and palate, with an incidence of 1 per 700 live births. Clubfoot or congenital talipes equinovarus (CTEV) is a three-dimensional abnormality of the leg, ankle, and feet that leads to the anomalous positioning of foot and ankle joints and has an incidence of 1 per 1000 live births. OFCs and CTEV may occur together or separately in certain genetic syndromes in addition to other congenital abnormalities. Here, we sought to decipher the genetic etiology of OFC and CTEV that occurred together in six probands.

Analyses of genome-wide association study (GWAS) data have revealed that detectable genetic mosaicism involving large (>2 Mb) structural autosomal alterations occurs in a fraction of individuals. We present results for a set of 24,849 genotyped individuals (total GWAS set II [TGSII]) in whom 341 large autosomal abnormalities were observed in 168 (0.68%) individuals. Merging data from the new TGSII set with data from two prior reports (the Gene-Environment Association Studies and the total GWAS set I) generated a large dataset of 127,179 individuals; we then conducted a meta-analysis to investigate the patterns of detectable autosomal mosaicism (n = 1,315 events in 925 [0.73%] individuals). Restricting to events >2 Mb in size, we observed an increase in event frequency as event size decreased. The combined results underscore that the rate of detectable mosaicism increases with age (p value = 5.5 × 10(-31)) and is higher in men (p value = 0.002) but lower in participants of African ancestry (p value = 0.003). In a subset of 47 individuals from whom serial samples were collected up to 6 years apart, complex changes were noted over time and showed an overall increase in the proportion of mosaic cells as age increased. Our large combined sample allowed for a unique ability to characterize detectable genetic mosaicism involving large structural events and strengthens the emerging evidence of non-random erosion of the genome in the aging population.

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.

Mutations in the transcription factor IRF6 cause allelic autosomal dominant clefting syndromes, Van der Woude syndrome, and popliteal pterygium syndrome. We compared the distribution of IRF6 coding and splice-site mutations from 549 families with Van der Woude syndrome or popliteal pterygium syndrome with that of variants from the 1000 Genomes and National Heart, Lung, and Blood Institute Exome Sequencing Projects.

The aim of this study was to identify genetic variants contributing to preterm birth (PTB) using a linkage candidate gene approach.

Preterm delivery (PTD) is an adverse birth outcome associated with increased infant mortality and negative lifelong health consequences. PTD may be the result of interactions between genetics and maternal/fetal environmental factors including smoking exposure (SMK). A common deletion in the GSTT1 gene was previously reported to affect birth outcomes in smokers. In this study, we dissect the associations among SMK, birth outcomes, and copy number variations (CNVs) in the GSTT1/GSTT2 region. A preterm birth case-control dataset of 1937 mothers was part of the GENEVA preterm birth study, which included genome-wide genotyping used to identify CNVs. We examined the association of SMK with birth outcomes, detected CNVs within the GSTT1/GSTT2 region using PennCNV, and examined associations of the identified CNVs with preterm birth and with birth weight (BW) in full term birth controls, including interactions with SMK. Finally, we tested the association of CNVs in GSTT1/GSTT2 with SMK. We confirmed the association of smoking with low BW and PTD. We identified 2 CNVs in GSTT2 (GSTT2 (a) and GSTT2 (b) ), 1 CNV in GSTTP1 and 2 CNVs in GSTT1 (GSTT1 (a) and GSTT1 (b) ). The GSTT2 (a) deletion was associated with reduced BW (-284 g, p = 2.50E-7) in smokers, and was more common in smokers [odds ratio(OR) = 1.30, p = 0.04]. We found that the size of the reported common deletion CNV in GSTT1 was larger than previously shown. The GSTTP1 and GSTT1 (b) null genotypes were in high linkage disequilibrium (LD) (D' = 0.89) and less common in smokers (OR = 0.68, p = 0.019 and OR = 0.73, p = 0.055, respectively). These two deletions were in partial LD with GSTT2 (a) and GSTT2 (b) duplications. All 5 CNVs seem to be associated with increased risk of preterm birth before 35 completed weeks. CNVs in the GSTTT1/GSTT2 region appear associated with low BW and PTD outcomes, but LD complicated these CNVs in GSTT1/GSTT2. In genetic association studies of BW, multiple CNVs in this region need to be investigated instead of a single polymorphism.

Preterm birth (PTB) is a complex disorder associated with significant neonatal mortality and morbidity and long-term adverse health consequences. Multiple lines of evidence suggest that genetic factors play an important role in its etiology. This study was designed to identify genetic variation associated with PTB in oxytocin pathway genes whose role in parturition is well known.

Epidermolytic ichthyosis (EI) is a rare skin disorder characterized by generalized erythroderma and cutaneous blistering at birth, which is substituted by hyperkeratosis later in life. It is caused by autosomal dominant mutations in highly conserved regions of KRT1 and KRT10. To date, only 4 mutations with autosomal recessive inheritance of EI have been described in consanguineous families. All of them affect the 2B domain of KRT10. In the present study we describe four patients with EI (including one lethal case) born from unaffected parents in a consanguineous family of a native Venezuelan community. The objective of this study was to characterize the clinical, genetic and morphological aspects of the disease in this family, as well as understand its functional implications. Genomic DNA was sequenced for KRT10 and KRT1. Immunofluoresence for keratin expression was performed on cutaneous biopsies. After examination of cutaneous biopsies histology, our results showed hyperkeratosis and acantholysis with an expanded granular layer. Sequencing of KRT10 demonstrated a non-sense mutation (p.Tyr282Ter.) corresponding to the 1B domain of the protein in patients and a heterozygous pattern in other family members, resulting in complete absence of K10. The loss of K10 was compensated by upregulation of K14 and K17. In conclusion, this novel mutation in KRT10 is the first recessive genetic variation that is not located in the so called "hot spot" for recessive EI, suggesting that other areas of the gene are also susceptible for such mutations.

There is a large literature showing the detrimental effects of prenatal smoking on birth and childhood health outcomes. It is somewhat unclear though, whether these effects are causal or reflect other characteristics and choices by mothers who choose to smoke that may also affect child health outcomes or biased reporting of smoking. In this paper we use genetic markers that predict smoking behaviors as instruments to address the endogeneity of smoking choices in the production of birth and childhood health outcomes. Our results indicate that prenatal smoking produces more dramatic declines in birth weight than estimates that ignore the endogeneity of prenatal smoking, which is consistent with previous studies with non-genetic instruments. We use data from two distinct samples from Norway and the United States with different measured instruments and find nearly identical results. The study provides a novel application that can be extended to study several behavioral impacts on health and social and economic outcomes.

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.

Axenfeld-Rieger syndrome (ARS) is associated with mutations in the PITX2 gene that encodes a homeobox transcription factor. Several intronic PITX2 mutations have been reported in Axenfeld-Rieger patients but their effects on gene expression have not been tested.

Maternal smoking is a recognized risk factor for orofacial clefts. Maternal or fetal pharmacogenetic variants are plausible modulators of this risk. In this work, we studied 5,427 DNA samples, including 1,244 from subjects in Denmark and Iowa with facial clefting and 4,183 from parents, siblings, or unrelated population controls. We examined 25 single-nucleotide polymorphisms in 16 genes in pathways for detoxification of components of cigarette smoke, to look for evidence of gene-environment interactions. For genes identified as related to oral clefting, we studied gene-expression profiles in fetal development in the relevant tissues and time intervals. Maternal smoking was a significant risk factor for clefting and showed dosage effects, in both the Danish and Iowan data. Suggestive effects of variants in the fetal NAT2 and CYP1A1 genes were observed in both the Iowan and the Danish participants. In an expanded case set, NAT2 continued to show significant overtransmission of an allele to the fetus, with a final P value of .00003. There was an interaction between maternal smoking and fetal inheritance of a GSTT1-null deletion, seen in both the Danish (P=.03) and Iowan (P=.002) studies, with a Fisher's combined P value of <.001, which remained significant after correction for multiple comparisons. Gene-expression analysis demonstrated expression of GSTT1 in human embryonic craniofacial tissues during the relevant developmental interval. This study benefited from two large samples, involving independent populations, that provided substantial power and a framework for future studies that could identify a susceptible population for preventive health care.