Genome-wide meta-analyses of nonsyndromic cleft lip with or without cleft palate identify six new risk loci

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; P(NSCLP) = 6.51 × 10(-11); homozygote relative risk = 2.41, 95% confidence interval (CI) 1.84-3.16).     

A common haplotype of interferon regulatory factor 5 (IRF5) regulates splicing and expression and is associated with increased risk of systemic lupus erythematosus

Systemic lupus erythematosus (SLE) is a complex autoimmune disease characterized by activation of the type I interferon (IFN) pathway. Here we convincingly replicate association of the IFN regulatory factor 5 (IRF5) rs2004640 T allele with SLE in four independent case-control cohorts (P = 4.4 x 10(-16)) and by family-based transmission disequilibrium test analysis (P = 0.0006). The rs2004640 T allele creates a 5' donor splice site in an alternate exon 1 of IRF5, allowing expression of several unique IRF5 isoforms. We also identify an independent cis-acting variant associated with elevated expression of IRF5 and linked to the exon 1B splice site. Haplotypes carrying the variant associated with elevated expression and lacking the exon 1B donor site do not confer risk of SLE. Thus, a common IRF5 haplotype driving elevated expression of multiple unique isoforms of IRF5 is an important genetic risk factor for SLE, establishing a causal role for type I IFN pathway genes in human autoimmunity.     

Mutation of PVRL1 is associated with sporadic, non-syndromic cleft lip/palate in northern Venezuela

Non-syndromic cleft lip with or without cleft palate (CL/P, MIM 119530) is among the most common of major birth defects. Homozygosity for a nonsense mutation of PVRL1, W185X, results in an autosomal recessive CL/P syndrome on Margarita Island, CLPED1 (ref. 1). Here we demonstrate highly significant association between heterozygosity for this mutation and sporadic, non-syndromic CL/P in northern Venezuela.     

Mutations of PVRL1, encoding a cell-cell adhesion molecule/herpesvirus receptor, in cleft lip/palate-ectodermal dysplasia

Cleft lip, with or without cleft palate (CL/P), is one of the most common birth defects, occurring in 0.4 to 2.0 per 1,000 infants born alive. Approximately 70% of CL/P cases are non-syndromic (MIM 119530), but CL/P also occurs in many single-gene syndromes, each affecting a protein critical for orofacial development. Here we describe positional cloning of the gene responsible for an autosomal recessive CL/P-ectodermal dysplasia (ED) syndrome (CLPED1; previously ED4; ref. 2), which we identify as PVRL1, encoding nectin-1, an immunoglobulin (Ig)-related transmembrane cell-cell adhesion molecule that is part of the NAP cell adhesion system. Nectin-1 is also the principal cell surface receptor for alpha-herpesviruses (HveC; ref. 7), and the high frequency of CLPED1 on Margarita Island in the Caribbean Sea might result from resistance of heterozygotes to infection by these viruses.     

Ryk-deficient mice exhibit craniofacial defects associated with perturbed Eph receptor crosstalk

Secondary palate formation is a complex process that is frequently disturbed in mammals, resulting in the birth defect cleft palate. Gene targeting has identified components of cytokine/growth factor signalling systems such as Tgf-alpha/Egfr, Eph receptors B2 and B3 (Ephb2 and Ephb3, respectively), Tgf-beta2, Tgf-beta3 and activin-betaA (ref. 3) as regulators of secondary palate development. Here we demonstrate that the mouse orphan receptor 'related to tyrosine kinases' (Ryk) is essential for normal development and morphogenesis of craniofacial structures including the secondary palate. Ryk belongs to a subclass of catalytically inactive, but otherwise distantly related, receptor protein tyrosine kinases (RTKs). Mice homozygous for a null allele of Ryk have a distinctive craniofacial appearance, shortened limbs and postnatal mortality due to feeding and respiratory complications associated with a complete cleft of the secondary palate. Consistent with cleft palate phenocopy in Ephb2/Ephb3-deficient mice and the role of a Drosophila melanogaster Ryk orthologue, Derailed, in the transduction of repulsive axon pathfinding cues, our biochemical data implicate Ryk in signalling mediated by Eph receptors and the cell-junction-associated Af-6 (also known as Afadin). Our findings highlight the importance of signal crosstalk between members of different RTK subfamilies.     

Epidermal growth factor receptor function is necessary for normal craniofacial development and palate closure.

Craniofacial malformations are among the most frequent congenital birth defects in humans; cleft palate, that is inadequate fusion of the palatal shelves, occurs with an annual incidence of 1 in 700 to 1 in 1,000 live births among individuals of European descent. The secondary palate arises as bilateral outgrowths from the maxillary processes, and its formation depends on the coordinated development of craniofacial structures including the Meckel's cartilage and the mandible. Cleft lip and palate syndromes in humans are associated with polymorphisms in the gene (TGFA) encoding transforming growth factor-alpha (TGF-alpha), an epidermal growth factor receptor (EGFR) ligand made by most epithelia. Here we have characterized craniofacial development in Egfr-deficient (Egfr-/-) mice. Newborn Egfr-/- mice have facial mediolateral defects including narrow, elongated snouts, underdeveloped lower jaw and a high incidence of cleft palate. Palatal shelf explants from Egfr-/- mice fused, but frequently had residual epithelium in the midline. In addition, morphogenesis of Meckel's cartilage was deficient in cultured mandibular processes from Egfr-/- embryos. The secretion of matrix metalloproteinases (MMPs) was diminished in Egfr-/- explants, consistent with the ability of EGF to increase MMP secretion and with the decreased MMP expression caused by inhibition of Egfr signalling in wild-type explants. Accordingly, inactivation of MMPs in wild-type explants phenocopied the defective morphology of Meckel's cartilage seen in Egfr-/- explants. Our results indicate that EGFR signalling is necessary for normal craniofacial development and that its role is mediated in part by its downstream targets, the MMPs, and may explain the genetic correlation of human cleft palate with polymorphisms in TGFA.     

MicroRNA Mirn140 modulates Pdgf signaling during palatogenesis

Disruption of signaling pathways such as those mediated by sonic hedgehog (Shh) or platelet-derived growth factor (Pdgf) causes craniofacial abnormalities, including cleft palate. The role that microRNAs play in modulating palatogenesis, however, is completely unknown. We show that, in zebrafish, the microRNA Mirn140 negatively regulates Pdgf signaling during palatal development, and we provide a mechanism for how disruption of Pdgf signaling causes palatal clefting. The pdgf receptor alpha (pdgfra) 3' UTR contained a Mirn140 binding site functioning in the negative regulation of Pdgfra protein levels in vivo. pdgfra mutants and Mirn140-injected embryos shared a range of facial defects, including clefting of the crest-derived cartilages that develop in the roof of the larval mouth. Concomitantly, the oral ectoderm beneath where these cartilages develop lost pitx2 and shha expression. Mirn140 modulated Pdgf-mediated attraction of cranial neural crest cells to the oral ectoderm, where crest-derived signals were necessary for oral ectodermal gene expression. Mirn140 loss of function elevated Pdgfra protein levels, altered palatal shape and caused neural crest cells to accumulate around the optic stalk, a source of the ligand Pdgfaa. These results suggest that the conserved regulatory interactions of mirn140 and pdgfra define an ancient mechanism of palatogenesis, and they provide candidate genes for cleft palate.     

Common variants near MC4R are associated with fat mass, weight and risk of obesity

To identify common variants influencing body mass index (BMI), we analyzed genome-wide association data from 16,876 individuals of European descent. After previously reported variants in FTO, the strongest association signal (rs17782313, P = 2.9 x 10(-6)) mapped 188 kb downstream of MC4R (melanocortin-4 receptor), mutations of which are the leading cause of monogenic severe childhood-onset obesity. We confirmed the BMI association in 60,352 adults (per-allele effect = 0.05 Z-score units; P = 2.8 x 10(-15)) and 5,988 children aged 7-11 (0.13 Z-score units; P = 1.5 x 10(-8)). In case-control analyses (n = 10,583), the odds for severe childhood obesity reached 1.30 (P = 8.0 x 10(-11)). Furthermore, we observed overtransmission of the risk allele to obese offspring in 660 families (P (pedigree disequilibrium test average; PDT-avg) = 2.4 x 10(-4)). The SNP location and patterns of phenotypic associations are consistent with effects mediated through altered MC4R function. Our findings establish that common variants near MC4R influence fat mass, weight and obesity risk at the population level and reinforce the need for large-scale data integration to identify variants influencing continuous biomedical traits.     

Interleukin 7 receptor alpha chain (IL7R) shows allelic and functional association with multiple sclerosis

Multiple sclerosis is a demyelinating neurodegenerative disease with a strong genetic component. Previous genetic risk studies have failed to identify consistently linked regions or genes outside of the major histocompatibility complex on chromosome 6p. We describe allelic association of a polymorphism in the gene encoding the interleukin 7 receptor alpha chain (IL7R) as a significant risk factor for multiple sclerosis in four independent family-based or case-control data sets (overall P = 2.9 x 10(-7)). Further, the likely causal SNP, rs6897932, located within the alternatively spliced exon 6 of IL7R, has a functional effect on gene expression. The SNP influences the amount of soluble and membrane-bound isoforms of the protein by putatively disrupting an exonic splicing silencer.     

Genetic variation near IRS1 associates with reduced adiposity and an impaired metabolic profile

Genome-wide association studies have identified 32 loci influencing body mass index, but this measure does not distinguish lean from fat mass. To identify adiposity loci, we meta-analyzed associations between ～2.5 million SNPs and body fat percentage from 36,626 individuals and followed up the 14 most significant (P < 10(-6)) independent loci in 39,576 individuals. We confirmed a previously established adiposity locus in FTO (P = 3 × 10(-26)) and identified two new loci associated with body fat percentage, one near IRS1 (P = 4 × 10(-11)) and one near SPRY2 (P = 3 × 10(-8)). Both loci contain genes with potential links to adipocyte physiology. Notably, the body-fat-decreasing allele near IRS1 is associated with decreased IRS1 expression and with an impaired metabolic profile, including an increased visceral to subcutaneous fat ratio, insulin resistance, dyslipidemia, risk of diabetes and coronary artery disease and decreased adiponectin levels. Our findings provide new insights into adiposity and insulin resistance.     

Myosin IXB variant increases the risk of celiac disease and points toward a primary intestinal barrier defect

Celiac disease is probably the best-understood immune-related disorder. The disease presents in the small intestine and results from the interplay between multiple genes and gluten, the triggering environmental factor. Although HLA class II genes explain 40% of the heritable risk, non-HLA genes accounting for most of the familial clustering have not yet been identified. Here we report significant and replicable association (P = 2.1 x 10(-6)) to a common variant located in intron 28 of the gene myosin IXB (MYO9B), which encodes an unconventional myosin molecule that has a role in actin remodeling of epithelial enterocytes. Individuals homozygous with respect to the at-risk allele have a 2.3-times higher risk of celiac disease (P = 1.55 x 10(-5)). This result is suggestive of a primary impairment of the intestinal barrier in the etiology of celiac disease, which may explain why immunogenic gluten peptides are able to pass through the epithelial barrier.     

Genome-wide association identifies a susceptibility locus for coronary artery disease in the Chinese Han population

Coronary artery disease (CAD) causes more than 700,000 deaths each year in China. Previous genome-wide association studies (GWAS) in populations of European ancestry identified several genetic loci for CAD, but no such study has yet been reported in the Chinese population. Here we report a three-stage GWAS in the Chinese Han population. We identified a new association between rs6903956 in a putative gene denoted as C6orf105 on chromosome 6p24.1 and CAD (P = 5.00 × 10?3, stage 2 validation; P = 3.00 × 10?3, P = 1.19 × 10?? and P = 4.00 × 10?3 in three independent stage 3 replication populations; P = 4.87 × 10?12, odds ratio = 1.51 in the combined population). The minor risk allele A of rs6903956 is associated with decreased C6orf105 mRNA expression. We report the first GWAS for CAD in the Chinese Han population and identify a SNP, rs6903956, in C6orf105 associated with susceptibility to CAD in this population.     

A genome-wide association study identifies six susceptibility loci for chronic lymphocytic leukemia

We conducted a genome-wide association study of 299,983 tagging SNPs for chronic lymphocytic leukemia (CLL) and performed validation in two additional series totaling 1,529 cases and 3,115 controls. We identified six previously unreported CLL risk loci at 2q13 (rs17483466; P = 2.36 x 10(-10)), 2q37.1 (rs13397985, SP140; P = 5.40 x 10(-10)), 6p25.3 (rs872071, IRF4; P = 1.91 x 10(-20)), 11q24.1 (rs735665; P = 3.78 x 10(-12)), 15q23 (rs7176508; P = 4.54 x 10(-12)) and 19q13.32 (rs11083846, PRKD2; P = 3.96 x 10(-9)). These data provide the first evidence for the existence of common, low-penetrance susceptibility to a hematological malignancy and new insights into disease causation in CLL.     

Genetic variation in PNPLA3 confers susceptibility to nonalcoholic fatty liver disease

Nonalcoholic fatty liver disease (NAFLD) is a burgeoning health problem of unknown etiology that varies in prevalence among ancestry groups. To identify genetic variants contributing to differences in hepatic fat content, we carried out a genome-wide association scan of nonsynonymous sequence variations (n = 9,229) in a population comprising Hispanic, African American and European American individuals. An allele in PNPLA3 (rs738409[G], encoding I148M) was strongly associated with increased hepatic fat levels (P = 5.9 x 10(-10)) and with hepatic inflammation (P = 3.7 x 10(-4)). The allele was most common in Hispanics, the group most susceptible to NAFLD; hepatic fat content was more than twofold higher in PNPLA3 rs738409[G] homozygotes than in noncarriers. Resequencing revealed another allele of PNPLA3 (rs6006460[T], encoding S453I) that was associated with lower hepatic fat content in African Americans, the group at lowest risk of NAFLD. Thus, variation in PNPLA3 contributes to ancestry-related and inter-individual differences in hepatic fat content and susceptibility to NAFLD.