MiRNA-149 as a Candidate for Facial Clefting and Neural Crest Cell Migration.

Nonsyndromic cleft lip with or without palate (nsCL/P) ranks among the most common human birth defects and has a multifactorial etiology. Human neural crest cells (hNCC) make a substantial contribution to the formation of facial bone and cartilage and are a key cell type in terms of nsCL/P etiology. Based on increasing evidence for the role of noncoding regulatory mechanisms in nsCL/P, we investigated the role of hNCC-expressed microRNAs (miRNA) in cleft development. First, we conducted a systematic analysis of miRNAs expressed in human-induced pluripotent stem cell-derived hNCC using Affymetrix microarrays on cell lines established from 4 unaffected donors. These analyses identified 152 candidate miRNAs. Based on the hypothesis that candidate miRNA loci harbor genetic variation associated with nsCL/P risk, the genomic locations of these candidates were cross-referenced with data from a previous genome-wide association study of nsCL/P. Associated variants were reanalyzed in independent nsCL/P study populations. Jointly, the results suggest that miR-149 is implicated in nsCL/P etiology. Second, functional follow-up included in vitro overexpression and inhibition of miR-149 in hNCC and subsequent analyses at the molecular and phenotypic level. Using 3'RNA-Seq, we identified 604 differentially expressed (DE) genes in hNCC overexpressing miR-149 compared with untreated cells. These included TLR4 and JUNB, which are established targets of miR-149, and NOG, BMP4, and PAX6, which are reported nsCL/P candidate genes. Pathway analyses revealed that DE genes were enriched in pathways including regulation of cartilage development and NCC differentiation. At the cellular level, distinct hNCC migration patterns were observed in response to miR-149 overexpression. Our data suggest that miR-149 is involved in the etiology of nsCL/P via its role in hNCC migration.

Common variants in DLG1 locus are associated with non-syndromic cleft lip with or without cleft palate.

Non-syndromic cleft lip with or without cleft palate (nsCL/P) is a common craniofacial anomaly with a complex and heterogeneous aetiology. Knowledge regarding specific genetic factors underlying this birth defect is still not well understood. Therefore, we conducted an independent replication analysis for the top-associated variants located within the DLG1 locus at chromosome 3q29, which was identified as a novel cleft-susceptibility locus in our genome-wide association study (GWAS). Mega-analysis of the pooled individual data from the GWAS and replication study confirmed that common DLG1 variants are associated with the risk of nsCL/P. Two single nucleotide polymorphisms (SNPs), rs338217 and rs7649443, were statistically significant even at the genome-wide level (Ptrend = 9.70E-10 and Ptrend = 8.96E-09, respectively). Three other SNPs, rs9826379, rs6805920 and rs6583202, reached a suggestive genome-wide significance threshold (Ptrend < 1.00E-05). The location of the strongest individual SNP in the intronic sequence of the gene encoding DLG1 antisense RNA suggests that the true causal variant implicated in the risk of nsCL/P may affect the DLG1 gene expression level rather than structure of the encoded protein. In conclusion, we identified a novel cleft-susceptibility locus at chromosome 3q29 with a DLG1 as a novel candidate gene for this common craniofacial anomaly.

MRPL53, a New Candidate Gene for Orofacial Clefting, Identified Using an eQTL Approach.

A valuable approach to understand how individual and population genetic differences can predispose to disease is to assess the impact of genetic variants on cellular functions (e.g., gene expression) of cell and tissue types related to pathological states. To understand the genetic basis of nonsyndromic cleft lip with or without cleft palate (NSCL/P) susceptibility, a complex and highly prevalent congenital malformation, we searched for genetic variants with a regulatory role in a disease-related tissue, the lip muscle (orbicularis oris muscle [OOM]), of affected individuals. From 46 OOM samples, which are frequently discarded during routine corrective surgeries on patients with orofacial clefts, we derived mesenchymal stem cells and correlated the individual genetic variants with gene expression from these cultured cells. Through this strategy, we detected significant cis-eQTLs (i.e., DNA variants affecting gene expression) and selected a few candidates to conduct an association study in a large Brazilian cohort (624 patients and 668 controls). This resulted in the discovery of a novel susceptibility locus for NSCL/P, rs1063588, the best eQTL for the MRPL53 gene, where evidence for association was mostly driven by the Native American ancestry component of our Brazilian sample. MRPL53 (2p13.1) encodes a 39S protein subunit of mitochondrial ribosomes and interacts with MYC, a transcription factor required for normal facial morphogenesis. Our study illustrates not only the importance of sampling admixed populations but also the relevance of measuring the functional effects of genetic variants over gene expression to dissect the complexity of disease phenotypes.

Candidate Genes for Nonsyndromic Cleft Palate Detected by Exome Sequencing.

Nonsyndromic cleft palate only (nsCPO) is a facial malformation that has a livebirth prevalence of 1 in 2,500. Research suggests that the etiology of nsCPO is multifactorial, with a clear genetic component. To date, genome-wide association studies have identified only 1 conclusive common variant for nsCPO, that is, a missense variant in the gene grainyhead-like-3 ( GRHL3). Thus, the underlying genetic causes of nsCPO remain largely unknown. The present study aimed at identifying rare variants that might contribute to nsCPO risk, via whole-exome sequencing (WES), in multiply affected Central European nsCPO pedigrees. WES was performed in 2 affected first-degree relatives from each family. Variants shared between both individuals were analyzed for their potential deleterious nature and a low frequency in the general population. Genes carrying promising variants were annotated for 1) reported associations with facial development, 2) multiple occurrence of variants, and 3) expression in mouse embryonic palatal shelves. This strategy resulted in the identification of a set of 26 candidate genes that were resequenced in 132 independent nsCPO cases and 623 independent controls of 2 different ethnicities, using molecular inversion probes. No rare loss-of-function mutation was identified in either WES or resequencing step. However, we identified 2 or more missense variants predicted to be deleterious in each of 3 genes ( ACACB, PTPRS, MIB1) in individuals from independent families. In addition, the analyses identified a novel variant in GRHL3 in 1 patient and a variant in CREBBP in 2 siblings. Both genes underlie different syndromic forms of CPO. A plausible hypothesis is that the apparently nonsyndromic clefts in these 3 patients might represent hypomorphic forms of the respective syndromes. In summary, the present study identified rare variants that might contribute to nsCPO risk and suggests candidate genes for further investigation.

The Role of Noncoding Genetic Variation in Isolated Orofacial Clefts.

In the past decade, medical genetic research has generated multiple discoveries, many of which were obtained via genome-wide association studies (GWASs). A major GWAS finding is that the majority of risk variants for complex traits map to noncoding regions. This has resulted in a paradigm shift in terms of the interpretation of human genomic sequence variation, with more attention now being paid to what was previously termed "junk DNA." Translation of genetic findings into biologically meaningful results requires 1) large-scale and cell-specific efforts to annotate non-protein-coding regions and 2) the integration of comprehensive genomic data sets. However, this represents an enormous challenge, particularly in the case of human traits that arise during embryonic development, such as orofacial clefts (OFCs). OFC is a multifactorial trait and ranks among the most common of all human congenital malformations. These 2 attributes apply in particular to its isolated forms (nonsyndromic OFC [nsOFC]). Although genetic studies (including GWASs) have yielded novel insights into the genetic architecture of nsOFC, few data are available concerning causality and affected biological pathways. Reasons for this deficiency include the complex genetic architecture at risk loci and the limited availability of functional data sets from human tissues that represent relevant embryonic sites and time points. The present review summarizes current knowledge of the role of noncoding regions in nsOFC etiology. We describe the identification of genetic risk factors for nsOFC and several of the approaches used to identify causal variants at these loci. These strategies include the use of biological and genetic information from public databases, the assessment of the full spectrum of genetic variability within 1 locus, and comprehensive in vitro and in vivo experiments. This review also highlights the role of the emerging research field "functional genomics" and its increasing contribution to our biological understanding of nsOFC.

Novel IRF6 Mutations Detected in Orofacial Cleft Patients by Targeted Massively Parallel Sequencing.

Common variants in interferon regulatory factor 6 ( IRF6) have been associated with nonsyndromic cleft lip with or without cleft palate (NSCL/P) as well as with tooth agenesis (TA). These variants contribute a small risk towards the 2 congenital conditions and explain only a small percentage of heritability. On the other hand, many IRF6 mutations are known to be a monogenic cause of disease for syndromic orofacial clefting (OFC). We hypothesize that IRF6 mutations in some rare instances could also cause nonsyndromic OFC. To find novel rare variants in IRF6 responsible for nonsyndromic OFC and TA, we performed targeted multiplex sequencing using molecular inversion probes (MIPs) in 1,072 OFC patients, 67 TA patients, and 706 controls. We identified 3 potentially pathogenic de novo mutations in OFC patients. In addition, 3 rare missense variants were identified, for which pathogenicity could not unequivocally be shown, as all variants were either inherited from an unaffected parent or the parental DNA was not available. Retrospective investigation of the patients with these variants revealed the presence of lip pits in one of the patients with a de novo mutation suggesting a Van der Woude syndrome (VWS) phenotype, whereas, in other patients, no lip pits were identified.

Strong association of variants around FOXE1 and orofacial clefting.

Nonsyndromic orofacial clefting (nsOFC) is a common, complex congenital disorder. The most frequent forms are nonsyndromic cleft lip with or without cleft palate (nsCL/P) and nonsyndromic cleft palate only (nsCPO). Although they are generally considered distinct entities, a recent study has implicated a region around the FOXE1 gene in both nsCL/P and nsCPO. To investigate this hypothesis, we analyzed the 2 most strongly associated markers (rs3758249 and rs4460498) in 2 independent samples of differing ethnicities: Central European (949 nsCL/P cases, 155 nsCPO cases, 1163 controls) and Mayan Mesoamerican (156 nsCL/P cases, 10 nsCPO cases, 338 controls). While highly significant associations for both single-nucleotide polymorphisms were obtained in nsCL/P (rs4460498: p Europe = 6.50 × 10(-06), p Mayan = .0151; rs3758249: p Europe = 2.41 × 10(-05), p Mayan = .0299), no association was found in nsCPO (p > .05). Genotyping of rs4460498 in 472 independent European trios revealed significant associations for nsCL/P (p = .016) and nsCPO (p = .043). A meta-analysis of all data revealed a genomewide significant result for nsCL/P (p = 1.31 × 10(-08)), which became more significant when nsCPO cases were added (p nsOFC = 1.56 × 10(-09)). These results strongly support the FOXE1 locus as a risk factor for nsOFC. With the data of the initial study, there is now considerable evidence that this locus is the first conclusive risk factor shared between nsCL/P and nsCPO.

Evidence for a polygenic contribution to androgenetic alopecia.

BACKGROUND: Male pattern baldness (androgenetic alopecia, AGA) is a highly heritable trait and the most common form of hair loss in humans. Eight genome-wide significant risk loci for AGA have been identified. OBJECTIVES: To determine whether a polygenic component contributes to the genetic risk for AGA. METHODS: This study used a German case-control sample for AGA, which comprised 581 severely affected patients and 617 controls, to determine the contribution of polygenic variance to AGA risk. The sample was divided evenly into discovery and test samples. An additive polygenic risk score was calculated from risk alleles with increasingly liberal P-values in the discovery dataset, which was then used to test for the enrichment of AGA risk score alleles in the independent test samples. RESULTS: The polygenic score analysis provided significant evidence for a polygenic contribution to AGA where the amount of variance explained was 1·4-4·5%. CONCLUSION: This study provides evidence for the specific contribution of a polygenic component to the overall heritable risk for AGA. To some degree, the polygenic architecture of AGA might reflect the complexity of the biological pathways involved. Further analyses and strategies that complement conventional genome-wide association studies are needed to identify these factors. These may include pathway-based analyses, the analysis of functional candidate genes and tests for epistatic effects with known loci.

Confirming genes influencing risk to cleft lip with/without cleft palate in a case-parent trio study.

A collection of 1,108 case-parent trios ascertained through an isolated, nonsyndromic cleft lip with or without cleft palate (CL/P) was used to replicate the findings from a genome-wide association study (GWAS) conducted by Beaty et al. (Nat Genet 42:525-529, 2010), where four different genes/regions were identified as influencing risk to CL/P. Tagging SNPs for 33 different genes were genotyped (1,269 SNPs). All four of the genes originally identified as showing genome-wide significance (IRF6, ABCA4 and MAF, plus the 8q24 region) were confirmed in this independent sample of trios (who were primarily of European and Southeast Asian ancestry). In addition, eight genes classified as 'second tier' hits in the original study (PAX7, THADA, COL8A1/FILIP1L, DCAF4L2, GADD45G, NTN1, RBFOX3 and FOXE1) showed evidence of linkage and association in this replication sample. Meta-analysis between the original GWAS trios and these replication trios showed PAX7, COL8A1/FILIP1L and NTN1 achieved genome-wide significance. Tests for gene-environment interaction between these 33 genes and maternal smoking found evidence for interaction with two additional genes: GRID2 and ELAVL2 among European mothers (who had a higher rate of smoking than Asian mothers). Formal tests for gene-gene interaction (epistasis) failed to show evidence of statistical interaction in any simple fashion. This study confirms that many different genes influence risk to CL/P.

A common variant in myosin-18B contributes to mathematical abilities in children with dyslexia and intraparietal sulcus variability in adults.

The ability to perform mathematical tasks is required in everyday life. Although heritability estimates suggest a genetic contribution, no previous study has conclusively identified a genetic risk variant for mathematical performance. Research has shown that the prevalence of mathematical disabilities is increased in children with dyslexia. We therefore correlated genome-wide data of 200 German children with spelling disability, with available quantitative data on mathematic ability. Replication of the top findings in additional dyslexia samples revealed that rs133885 was a genome-wide significant marker for mathematical abilities (P(comb) = 7.71 × 10(-10), n = 699), with an effect size of 4.87%. This association was also found in a sample from the general population (P = 0.048, n = 1080), albeit with a lower effect size. The identified variant encodes an amino-acid substitution in MYO18B, a protein with as yet unknown functions in the brain. As areas of the parietal cortex, in particular the intraparietal sulcus (IPS), are involved in numerical processing in humans, we investigated whether rs133885 was associated with IPS morphology using structural magnetic resonance imaging data from 79 neuropsychiatrically healthy adults. Carriers of the MYO18B risk-genotype displayed a significantly lower depth of the right IPS. This validates the identified association between rs133885 and mathematical disability at the level of a specific intermediate phenotype.

Evidence for the involvement of ZNF804A in cognitive processes of relevance to reading and spelling.

Previous studies have shown that individuals with schizophrenia and dyslexia display common neurocognitive abnormalities. The aim of the present study was to determine whether known schizophrenia-risk genes contribute to dyslexia risk or to disease-relevant cognitive functions. For this purpose, we genotyped the schizophrenia-associated risk variants within zinc-finger protein 804A (ZNF804A), transcription-factor 4 and neurogranin in a large dyslexia case-control sample. We tested all variants for association with dyslexia (927 cases, 1096 controls), and with eight language-relevant cognitive processes (1552 individuals). We observed six significant associations between language-relevant traits and the ZNF804A-variant rs1344706. Interestingly, the ZNF804A schizophrenia risk variant was associated with a better cognitive performance in our data set. This finding might be consistent with a previously reported ZNF804A association in schizophrenia, in which patients carrying the schizophrenia-risk allele at rs1344706 showed a better performance in two memory tests. In conclusion, the present study provides evidence that ZNF804A might have a role in cognitive traits of relevance to reading and spelling, and underlines the phenotypic complexity that might be associated with ZNF804A.

First genome-wide association scan on neurophysiological endophenotypes points to trans-regulation effects on SLC2A3 in dyslexic children.

Dyslexia is one of the most common learning disorders affecting about 5% of all school-aged children. It has been shown that event-related potential measurements reveal differences between dyslexic children and age-matched controls. This holds particularly true for mismatch negativity (MMN), which reflects automatic speech deviance processing and is altered in dyslexic children. We performed a whole-genome association analysis in 200 dyslexic children, focusing on MMN measurements. We identified rs4234898, a marker located on chromosome 4q32.1, to be significantly associated with the late MMN component. This association could be replicated in an independent second sample of 186 dyslexic children, reaching genome-wide significance in the combined sample (P = 5.14e-08). We also found an association between the late MMN component and a two-marker haplotype of rs4234898 and rs11100040, one of its neighboring single nucleotide polymorphisms (SNPs). In the combined sample, this marker combination withstands correction for multiple testing (P = 6.71e-08). Both SNPs lie in a region devoid of any protein-coding genes; however, they both show significant association with mRNA-expression levels of SLC2A3 on chromosome 12, the predominant facilitative glucose transporter in neurons. Our results suggest a possible trans-regulation effect on SLC2A3, which might lead to glucose deficits in dyslexic children and could explain their attenuated MMN in passive listening tasks.