Genome-wide Enrichment of De Novo Coding Mutations in Orofacial Cleft Trios.

Although de novo mutations (DNMs) are known to increase an individual's risk of congenital defects, DNMs have not been fully explored regarding orofacial clefts (OFCs), one of the most common human birth defects. Therefore, whole-genome sequencing of 756 child-parent trios of European, Colombian, and Taiwanese ancestry was performed to determine the contributions of coding DNMs to an individual's OFC risk. Overall, we identified a significant excess of loss-of-function DNMs in genes highly expressed in craniofacial tissues, as well as genes associated with known autosomal dominant OFC syndromes. This analysis also revealed roles for zinc-finger homeobox domain and SOX2-interacting genes in OFC etiology.

EIF2AK2 Missense Variants Associated with Early Onset Generalized Dystonia.

OBJECTIVE: The study was undertaken to identify a monogenic cause of early onset, generalized dystonia. METHODS: Methods consisted of genome-wide linkage analysis, exome and Sanger sequencing, clinical neurological examination, brain magnetic resonance imaging, and protein expression studies in skin fibroblasts from patients. RESULTS: We identified a heterozygous variant, c.388G>A, p.Gly130Arg, in the eukaryotic translation initiation factor 2 alpha kinase 2 (EIF2AK2) gene, segregating with early onset isolated generalized dystonia in 5 patients of a Taiwanese family. EIF2AK2 sequencing in 191 unrelated patients with unexplained dystonia yielded 2 unrelated Caucasian patients with an identical heterozygous c.388G>A, p.Gly130Arg variant, occurring de novo in one case, another patient carrying a different heterozygous variant, c.413G>C, p.Gly138Ala, and one last patient, born from consanguineous parents, carrying a third, homozygous variant c.95A>C, p.Asn32Thr. These 3 missense variants are absent from gnomAD, and are located in functional domains of the encoded protein. In 3 patients, additional neurological manifestations were present, including intellectual disability and spasticity. EIF2AK2 encodes a kinase (protein kinase R [PKR]) that phosphorylates eukaryotic translation initiation factor 2 alpha (eIF2α), which orchestrates the cellular stress response. Our expression studies showed abnormally enhanced activation of the cellular stress response, monitored by PKR-mediated phosphorylation of eIF2α, in fibroblasts from patients with EIF2AK2 variants. Intriguingly, PKR can also be regulated by PRKRA (protein interferon-inducible double-stranded RNA-dependent protein kinase activator A), the product of another gene causing monogenic dystonia. INTERPRETATION: We identified EIF2AK2 variants implicated in early onset generalized dystonia, which can be dominantly or recessively inherited, or occur de novo. Our findings provide direct evidence for a key role of a dysfunctional eIF2α pathway in the pathogenesis of dystonia. ANN NEUROL 2021;89:485-497.

The SNP rs560426 Within ABCA4-ARHGAP29 Locus and the Risk of Nonsyndromic Oral Clefts.

OBJECTIVE: Nonsyndromic oral clefts are common birth defect with complex etiology. In the present study, we attempt to further validate the possible role for ABCA4 and ARHGAP29 in the susceptibility to nonsyndromic oral clefts. DESIGN: We performed allelic transmission disequilibrium test analysis, on 10 eligible single nucleotide polymorphisms (SNPs) and SNP haplotypes using the Family-Based Association Test. PARTICIPANTS: The study sample consisted of 334 case-parent trios of nonsyndromic oral clefts from Taiwanese population, separated into nonsyndromic cleft lip with or without cleft palate (NSCL/P) and nonsyndromic cleft palate only (NSCPO) groups. RESULTS: We found only the SNP rs560426 within the ABCA4 gene showed strong association with NSCPO (P = .03498; Permuted P = .05382). No association between other 9 selected SNPs in ABCA4-ARHGAP29 region and the risk of nonsyndromic oral clefts was found. For the haplotype analyses, we found only haplotype T-C (rs570926 and rs3789431) in ABCA4 block 2 showed significant association with nonsyndromic NSCL/P in these Taiwanese trios. CONCLUSIONS: We used a family-based analysis in 334 Taiwanese case-parent trios to validate the possible role for ABCA4 and ARHGAP29 in the susceptibility to nonsyndromic oral clefts. This study provides a new evidence for an association between the intron variant rs560426 within ABCA4 and nonsyndromic cleft palate which may contribute their regulatory role in craniofacial development.

Association Studies Between Regulatory Regions of IRF6/TP63 Genes and Nonsyndromic Oral Clefts.

OBJECTIVE: To evaluate genetic variants within the regulatory regions of interferon regulatory factor 6 (IRF6) and TP63 for the etiology of nonsyndromic oral clefts risk factors. DESIGN: We performed allelic transmission disequilibrium test analysis on 5 eligible single-nucleotide polymorphisms (SNPs) and SNP haplotypes using the Family-Based Association Test. PARTICIPANTS: The study sample consisted of 334 case-parent trios of nonsyndromic oral clefts from Taiwanese population, separated into nonsyndromic cleft lip/palate (NSCL/P) and nonsyndromic cleft palate only (NSCPO) groups. RESULTS: We found all 3 selected SNPs of the IRF6 gene show significant association with nonsyndromic oral clefts (rs2235371, P = 5.10E-07; rs642961, P = .00194; and rs77542756, P = 9.08E-07). Haplotype analyses identified 3 possible SNP combination haplotypes in the IRF6 gene and found that C-G-G showed significant undertransmission (P = .058), whereas 2 other haplotypes, T-G-A and C-A-G (P = 2.71E-06 and P = 5.00E-04, respectively), were significantly overtransmitted to the NSCL/P children but not to the NSCPO children. For the TP63 gene, we failed to detect evidence of nonsyndromic oral cleft association in the 2 SNPs within the TP63 large intron 1 region. CONCLUSIONS: We used a family-based analysis in 334 Taiwanese case-parent trios to evaluate selected SNPs of IRF6 genes and TP63 genes for a risk of orofacial clefting. This study provides additional evidence for an association between IRF6 and NSCL/P, including the genetic variants within the 5'-noncoding region of the gene. We also confirmed that NSCL/P and NSCPO individuals belong to different groups. For the TP63, our data did not favor the direct involvement of TAp63 isoforms during orofacial development.

Targeted Next Generation Sequencing for Genetic Mutations of Dilated Cardiomyopathy.

BACKGROUND: Approximately one-third of cases of dilated cardiomyopathy (DCM) are caused by genetic mutations. With new sequencing technologies, numerous variants have been associated with this inherited cardiomyopathy, however the prevalence and genotype-phenotype correlations in different ethnic cohorts remain unclear. This study aimed to investigate the variants in Chinese DCM patients and correlate them with clinical presentations and prognosis. METHODS AND RESULTS: From September 2013 to December 2016, 70 index patients underwent DNA sequencing for 12 common disease-causing genes with next generation sequencing. Using a bioinformatics filtering process, 12 rare truncating variants (7 nonsense variants, 4 frameshift variants, and 1 splice site variant) and 29 rare missense variants were identified. Of these, 3 patients were double heterozygotes and 10 patients were compound heterozygotes. Overall, 47.1% (33/70) of the index patients had the seputatively pathogenic variants. The majority (33/41, 80.4%) of these variants were located in titin (TTN). More than 80% of the TTN variants (27/33, 81.8%) were distributed in the A band region of the sarcomere. Patients carrying these variants did not have a different phenotype in disease severity, clinical outcome and reversibility of ventricular function compared with non-carriers. CONCLUSIONS: Several new rare variants were identified in a Chinese population in this study, indicating that there are ethnic differences in genetic mutations in DCM patients. TTN remains the major disease-causing gene. Our results could be a reference for future genetic tests in Chinese populations. No specific genotype-phenotype correlations were found, however a prospective large cohort study may be needed to confirm our findings.

Genetic diagnosis of neurofibromatosis type 1: targeted next- generation sequencing with Multiple Ligation-Dependent Probe Amplification analysis.

BACKGROUND: Neurofibromatosis type 1 (NF1) is a dominantly inherited tumor predisposition syndrome that targets the peripheral nervous system. It is caused by mutations of the NF1 gene which serve as a negative regulator of the cellular Ras/MAPK (mitogen-activated protein kinases) signaling pathway. Owing to the complexity in some parts of clinical diagnoses and the need for better understanding of its molecular relationships, a genetic characterization of this disorder will be helpful in the clinical setting. METHODS: In this study, we present a customized targeted gene panel of NF1/KRAS/BRAF/p53 and SPRED1 genes combined with Multiple Ligation-Dependent Probe Amplification analysis for the NF1 mutation screening in a cohort of patients clinically suspected as NF1. RESULTS: In this study, we identified 73 NF1 mutations and two BRAF novel variants from 100 NF1 patients who were suspected as having NF1. These genetic alterations are heterogeneous and distribute in a complicated way without clustering in either cysteine-serine-rich domain or within the GAP-related domain. We also detected fifteen multi-exon deletions within the NF1 gene by MLPA Analysis. CONCLUSIONS: Our results suggested that a genetic screening using a NGS panel with high coverage of Ras-signaling components combined with Multiple Ligation-Dependent Probe Amplification analysis will enable differential diagnosis of patients with overlapping clinical features.

SLC22A5 Mutations in a Patient With Systemic Primary Carnitine Deficiency and Cleft Palate-Successful Perioperative Management.

BACKGROUND: Primary systemic carnitine deficiency (SCD) is an autosomal-recessive disorder caused by SLC22A5 gene mutation resulting in defective cellular carnitine transporter organic cation transporter 2. Defective carnitine transporter causes renal carnitine wasting and low serum carnitine. Carnitine is an essential cofactor for the transportation of long-chain fatty acids into the mitochondria. Lacking of carnitine may cause metabolic decompensation and sudden death when the patient is exposed to prolonged fasting before an operation. METHODS: An asymptomatic 9-month-old boy with SCD diagnosed by local hospital was referred to the authors' hospital for incomplete cleft palate plastic surgery. Due to potential metabolic decompensation from prolonged fasting before the surgery, the patient underwent proper perioperative management. RESULTS: The operation was successful and subsequent clinical course was fine. The patient was discharged on postoperative day 3. CONCLUSION: With proper perioperative management, patients with SCD and cleft palate can survive from prolonged fasting time before and during operation without metabolic decompensation manifestations. Early recognition of SCD and perioperative management can be lifesaving in preoperative infants with SCD.

Gene-gene interaction of single nucleotide polymorphisms in 16p13.3 may contribute to the risk of non-syndromic cleft lip with or without cleft palate in Chinese case-parent trios.

Non-syndromic cleft lip with or without cleft palate (NSCL/P) is a common birth defect with a complex and heterogeneous etiology. A recent genome-wide association study (GWAS) among Chinese populations has identified a new region at 16p13.3 as being associated with NSCL/P, which requires further replication. Here, we attempted to replicate and further clarify the genetic association between this region and NSCL/P, as well as testing for potential gene-gene (G × G) and gene-environment (G × E) interactions. We conducted transmission disequilibrium tests on 69 single nucleotide polymorphisms (SNPs) mapping to 16p13.3 among 806 Chinese case-parent trios ascertained through an international consortium where a GWAS of oral clefts was conducted. G × G, as well as G × E interactions involving maternal environmental tobacco smoke (ETS) and multivitamin supplementation, were explored using conditional logistic regression model. We applied Cordell's method as implemented in the R package TRIO to test for possible interactions. While no SNPs showed evidence of linkage and association with NSCL/P after Bonferroni correction, we found signals of G × G interactions between SNPs in 16p13.3. Nine pairs of SNP-SNP interactions attained significance after Bonferroni correction, among which the most significant interaction was found between rs2072346 (ADCY9) and rs11646137 (intergenic region, P = 7.2 × 10-5 ). Linkage disequilibrium (LD) analysis revealed only low level of LD between these SNPs. This study failed to confirm the significant association between SNPs within 16p13.3 and the risk of NSCL/P, but underlined the importance of taking into account potential G × G interactions for the genetic association analysis of NSCL/P.

Nonsynonymous variants in MYH9 and ABCA4 are the most frequent risk loci associated with nonsyndromic orofacial cleft in Taiwanese population.

BACKGROUND: Nonsyndromic orofacial cleft is a common birth defect with a complex etiology, including multiple genetic and environmental risk factors. Recent whole genome analyses suggested associations between nonsyndromic orofacial cleft and up to 18 genetic risk loci (ABCA4, BMP4, CRISPLD2, GSTT1, FGF8, FGFR2, FOXE1, IRF6, MAFB, MSX1, MTHFR, MYH9, PDGFC, PVRL1, SUMO1, TGFA, TGFB3, and VAX1), each of which confers a different relative risk in different populations. We evaluate the nonsynonymous variants in these 18 genetic risk loci in nonsyndromic orofacial clefts and normal controls to clarify the specific variants in Taiwanese population. METHODS: We evaluated these 18 genetic risk loci in 103 cases of nonsyndromic orofacial clefts and 100 normal controls using a next-generation sequencing (NGS) customized panel and manipulated a whole-exon targeted-sequencing study based on the NGS system of an Ion Torrent Personal Genome Machine (IT-PGM). IT-PGM data processing, including alignment with the human genome build 19 reference genome (hg19), base calling, trimming of barcoded adapter sequences, and filtering of poor signal reads, was performed using the IT platform-specific pipeline software Torrent Suite, version 4.2, with the plug-in "variant caller" program. Further advanced annotation was facilitated by uploading the exported VCF file from Variant Caller to the commercial software package Ion Reporter; the free online annotation software Vanno and Mutation Taster. Benign or tolerated amino acid changes were excluded after analysis using sorting intolerant from tolerant and polymorphism phenotyping. Sanger sequencing was used to validate the significant variants identified by NGS. Furthermore, each variant was confirmed in asymptomatic controls using the Sequenom MassARRAY (San Diego, CA, USA). RESULTS: We identified totally 22 types of nonsynonymous variants specific in nonsyndromic orofacial clefts, including 19 single nucleotide variants, 2 deletions, and 1 duplication in 10 studied genes(ABCA4, MYH9, MTHFR, CRISPLD2, FGF8, PVRL1, FOXE1, VAX1, FGFR2, and IRF6). Nonsynonymous variants in MYH9 and ABCA4, which were detected in 6 and 5 individuals, respectively, were identified to be the most frequent risk loci in nonsyndromic orofacial clefts in the Taiwanese population. CONCLUSIONS: Nonsynonymous variants in MYH9 and ABCA4 were identified to be the most frequent risk loci in nonsyndromic orofacial clefts in the Taiwanese population. These findings in our study have provided additional information regarding specific variants associated with nonsyndromic orofacial clefts in different population and demonstrate the power of our customized NGS panel, which is clinically useful for the simultaneous detection of multiple genes associated with nonsyndromic orofacial clefts.

Paroxysmal exercise-induced dystonia within the phenotypic spectrum of ECHS1 deficiency.

BACKGROUND: ECHS1 encodes a mitochondrial enzyme involved in the degradation of essential amino acids and fatty acids. Recently, ECHS1 mutations were shown to cause a new severe metabolic disorder presenting as Leigh or Leigh-like syndromes. The objective of this study was to describe a family with 2 siblings affected by different dystonic disorders as a resulting phenotype of ECHS1 mutations. METHODS: Clinical evaluation, MRI imaging, genome-wide linkage, exome sequencing, urine metabolite profiling, and protein expression studies were performed. RESULTS: The first sibling is 17 years old and presents with generalized dystonia and severe bilateral pallidal MRI lesions after 1 episode of infantile subacute metabolic encephalopathy (Leigh-like syndrome). In contrast, the younger sibling (15 years old) only suffers from paroxysmal exercise-induced dystonia and has very mild pallidal MRI abnormalities. Both patients carry compound heterozygous ECHS1 mutations: c.232G>T (predicted protein effect: p.Glu78Ter) and c.518C>T (p.Ala173Val). Linkage analysis, exome sequencing, cosegregation, expression studies, and metabolite profiling support the pathogenicity of these mutations. Expression studies in patients' fibroblasts showed mitochondrial localization and severely reduced levels of ECHS1 protein. Increased urinary S-(2-carboxypropyl)cysteine and N-acetyl-S-(2-carboxypropyl)cysteine levels, proposed metabolic markers of this disorder, were documented in both siblings. Sequencing ECHS1 in 30 unrelated patients with paroxysmal dyskinesias revealed no further mutations. CONCLUSIONS: The phenotype associated with ECHS1 mutations might be milder than reported earlier, compatible with prolonged survival, and also includes isolated paroxysmal exercise-induced dystonia. ECHS1 screening should be considered in patients with otherwise unexplained paroxysmal exercise-induced dystonia, in addition to those with Leigh and Leigh-like syndromes. Diet regimens and detoxifying agents represent potential therapeutic strategies. © 2016 International Parkinson and Movement Disorder Society.

[Association study between candidate genes on transforming growth factor-ß signaling pathway and the risk of non-syndromic cleft lip with or without cleft palate in Chinese populations].

OBJECTIVE: To explore the association between 10 candidate genes on transforming growth factor-ß (TGFB) signaling pathway and non-syndromic cleft lip with or without cleft palate (NSCL/P) among Chinese populations, and to study the gene-environment interaction. METHODS: A total of 806 Chinese Han NSCL/P trios were ascertained from an international consortium, which conducted a genome-wide association study using a case-parent trio design to investigate the genes affecting risk to NSCL/P. The transmission disequilibrium test (TDT) was used to test for effects of 343 single nucleotide polymorphisms (SNPs) in 10 genes on TGFB signaling pathway including DCN, TGFB1, TGFB2, TGFB3, TGFBR1, TGFBR2, BAMBI, SMAD2, SMAD3 and SMAD4. The conditional regression models were used to test for gene-environment interaction. RESULTS: For TDT, although 19 SNPs showed nominal significant association with NSCL/P, no significant evidence of association was seen for all SNPs in 806 NSCL/P trios after Bonferroni correction. The interactions between genes and maternal smoking, environmental tobacco smoke, alcohol consumption and multi-vitamin supplementation during pregnancy did not attain statistical significance after correction for multiple comparisons. CONCLUSION: No evidence for SNP effect of genes on TGFB signaling pathway and significant gene-environment interaction was seen in our data.

Lower lip pits: van der woude or kabuki syndrome?

Kabuki syndrome (KS) is a multiple congenital anomaly/mental retardation syndrome with characteristic facial features. Despite more than 350 documented cases and recent correlation of MLL2 mutations as a genetic cause, its full clinical spectrum is still being defined. This report describes two patients who were initially diagnosed with Van der Woude syndrome (VWS) based on the presence of lower lip pits. However, this finding can occur with KS, albeit infrequently. For patients with lower lip pits, a thorough evaluation should be made to distinguish between VWS and KS, as there are differences in long-term prognosis.

Examining markers in 8q24 to explain differences in evidence for association with cleft lip with/without cleft palate between Asians and Europeans.

In a recent genome-wide association study (GWAS) from an international consortium, evidence of linkage and association in chr8q24 was much stronger among nonsyndromic cleft lip/palate (CL/P) case-parent trios of European ancestry than among trios of Asian ancestry. We examined marker information content and haplotype diversity across 13 recruitment sites (from Europe, United States, and Asia) separately, and conducted principal components analysis (PCA) on parents. As expected, PCA revealed large genetic distances between Europeans and Asians, and a north-south cline from Korea to Singapore in Asia, with Filipino parents forming a somewhat distinct Southeast Asian cluster. Hierarchical clustering of SNP heterozygosity revealed two major clades consistent with PCA results. All genotyped SNPs giving P < 10(-6) in the allelic transmission disequilibrium test (TDT) showed higher heterozygosity in Europeans than Asians. On average, European ancestry parents had higher haplotype diversity than Asians. Imputing additional variants across chr8q24 increased the strength of statistical evidence among Europeans and also revealed a significant signal among Asians (although it did not reach genome-wide significance). Tests for SNP-population interaction were negative, indicating the lack of strong signal for 8q24 in families of Asian ancestry was not due to any distinct genetic effect, but could simply reflect low power due to lower allele frequencies in Asians.

A combined targeted mutation analysis of IRF6 gene would be useful in the first screening of oral facial clefts.

BACKGROUND: Interferon Regulatory Factor 6 (IRF6) is a member of the IRF family of transcription factors. It has been suggested to be an important contributor to orofacial development since mutations of the IRF6 gene has been found in Van der Woude (VWS) and popliteal pterygium syndromes (PPS), two disorders that can present with isolated cleft lip and palate. The association between IRF6 gene and cleft lip and palate has also been independently replicated in many populations. METHODS: We screened a total of 155 Taiwanese patients with cleft lip with or without cleft palate (CL/P); 31 syndromic (including 19 VWS families), 44 non-syndromic families with at least two affected members, and 80 non-syndromic patients through a combined targeted, polymerase chain reaction (PCR)-based mutation analysis for the entire coding regions of IRF6 gene. RESULTS: We found 11 mutations in 57.89% (11/19) of the VWS patients and no IRF6 mutation in 44 of the non-syndromic multiplex families and 80 non-syndromic oral cleft patients. In this IRF6 gene screening, five of these mutations (c.290 A>G, p.Tyr97Cys; c.360-375 16 bp deletion, p.Gln120HisfsX24; c.411_412 insA, p.Glu136fsX3; c.871 A>C, p.Thr291Pro; c.969 G>A, and p.Trp323X) have not been reported in the literature previously. Exon deletion was not detected in this series of IRF6 gene screening. CONCLUSIONS: Our results confirm the crucial role of IRF6 in the VWS patients and further work is needed to explore for its function in the non-syndromic oral cleft with vary clinical features.

Novel PKD1 and PKD2 mutations in Taiwanese patients with autosomal dominant polycystic kidney disease.

Autosomal dominant polycystic kidney disease (ADPKD) is a heterogeneous disease caused by mutations in PKD1 and PKD2. The genotype-phenotype correlations are not completely understood. We performed direct PCR-sequencing plus multiplex ligation-dependent probe amplification for PKD1 and PKD2 in 46 unrelated patients. Disease-causing mutations were identified in 30 (65%) patients: 23 (77%) patients have mutations in PKD1 and 7 (23%) have mutations in PKD2. Nonsense, splicing or frame-shifting mutations were found in 18 patients, exon duplication in 1 and missense mutations in 11 patients. Two likely PKD1 hypomorphic alleles (p.Arg2477His and p.Arg3439Trp) segregated with mild disease in a family. A total of 34 mutations were identified and 17 (50%) of which are novel. The median age at onset of dialysis was significantly earlier in patients with PKD1 mutations (52 years) than in patients with PKD2 mutations (65.5 years) and those with an undetermined genotype (67 years) by survival analysis (log-rank test, P=0.014). Patients carrying PKD1-truncating mutations have a trend toward earlier initiation of dialysis compared with carriers of non-truncating mutations (52 years vs 57 years, P=0.061). A family history of dialysis before 55 years was more common in PKD1 patients than in others (P<0.05). In conclusion, this study identified novel mutations in PKD1 and PKD2 and demonstrated the presence of PKD1 hypomorphic alleles in Taiwanese patients. Patients carrying PKD1 mutations, especially those with truncating mutations, could have a more rapidly progressive disease than others. These results might have implications for diagnosis and risk stratification in patients with ADPKD.

X-linked markers in the Duchenne muscular dystrophy gene associated with oral clefts.

As part of an international consortium, case-parent trios were collected for a genome-wide association study of isolated, non-syndromic oral clefts, including cleft lip (CL), cleft palate (CP), and cleft lip and palate (CLP). Non-syndromic oral clefts have a complex and heterogeneous etiology. Risk is influenced by genes and environmental factors, and differs markedly by gender. Family-based association tests (FBAT) were used on 14,486 single nucleotide polymorphisms (SNPs) spanning the X chromosome, stratified by type of cleft and racial group. Significant results, even after multiple-comparisons correction, were obtained for the Duchenne muscular dystrophy (DMD) gene, the largest single gene in the human genome, among CL/P (i.e., both CL and CLP combined) trios. When stratified into groups of European and Asian ancestry, stronger signals were obtained for Asian subjects. Although conventional sliding-window haplotype analysis showed no increase in significance, selected combinations of the 25 most significant SNPs in the DMD gene identified four SNPs together that attained genome-wide significance among Asian CL/P trios, raising the possibility of interaction between distant SNPs within the DMD gene.

[Association study for gene polymorphism of folic acid/homocysteine metabolic pathway and nonsyndromic cleft lip with or without cleft palate in Chinese populations].

OBJECTIVE: To explore the association between 18 candidate genes encoding enzymes on the folate/homocysteine metabolism pathway and non-syndromic cleft lip with or without cleft palate (NSCL/P) in Chinese populations. METHODS: A total of 806 NSCL/P trios were drawn by an international consortium, which conducted a genome-wide association study using a case-parent trio design to investigate genes affecting risks to NSCL/P. The transmission disequilibrium test (TDT) was used for deviation from Mendelian expectations for 257 SNPs in 18 folate/homocysteine metabolism-related genes. The interactions between markers in these gene and environmental risk factors were also tested using conditional Logistic regressions. RESULTS: Although four SNPs (rs6428977, rs12060264, rs7730643 and rs4920037) showed nominal significant association with NSCL/P in the TDT on 806 NSCL/P trios (P<0.05), no significant evidence of linkage and association remained in all the SNPs after Bonferroni correction. Similar tests for interactions between genes and maternal smoking, environmental tobacco smoke, alcohol consumption and multi-vitamin supplementation during pregnancy did not attain statistical significance after correction for multiple comparisons. CONCLUSION: Folate/homocysteine metabolism-related genes could not influence the risk of NSCL/P.

Evidence for gene-environment interaction in a genome wide study of nonsyndromic cleft palate.

Nonsyndromic cleft palate (CP) is a common birth defect with a complex and heterogeneous etiology involving both genetic and environmental risk factors. We conducted a genome-wide association study (GWAS) using 550 case-parent trios, ascertained through a CP case collected in an international consortium. Family-based association tests of single nucleotide polymorphisms (SNP) and three common maternal exposures (maternal smoking, alcohol consumption, and multivitamin supplementation) were used in a combined 2 df test for gene (G) and gene-environment (G × E) interaction simultaneously, plus a separate 1 df test for G × E interaction alone. Conditional logistic regression models were used to estimate effects on risk to exposed and unexposed children. While no SNP achieved genome-wide significance when considered alone, markers in several genes attained or approached genome-wide significance when G × E interaction was included. Among these, MLLT3 and SMC2 on chromosome 9 showed multiple SNPs resulting in an increased risk if the mother consumed alcohol during the peri-conceptual period (3 months prior to conception through the first trimester). TBK1 on chr. 12 and ZNF236 on chr. 18 showed multiple SNPs associated with higher risk of CP in the presence of maternal smoking. Additional evidence of reduced risk due to G × E interaction in the presence of multivitamin supplementation was observed for SNPs in BAALC on chr. 8. These results emphasize the need to consider G × E interaction when searching for genes influencing risk to complex and heterogeneous disorders, such as nonsyndromic CP.

Genome wide study of maternal and parent-of-origin effects on the etiology of orofacial clefts.

We performed a genome wide association analysis of maternally-mediated genetic effects and parent-of-origin (POO) effects on risk of orofacial clefting (OC) using over 2,000 case-parent triads collected through an international cleft consortium. We used log-linear regression models to test individual SNPs. For SNPs with a P-value <10(-5) for maternal genotypic effects, we also applied a haplotype-based method, TRIMM, to extract potential information from clusters of correlated SNPs. None of the SNPs were significant at the genome wide level. Our results suggest neither maternal genome nor POO effects play major roles in the etiology of OC in our sample. This finding is consistent with previous genetic studies and recent population-based cohort studies in Norway and Denmark, which showed no apparent difference between mother-to-offspring and father-to-offspring recurrence of clefting. We, however, cannot completely rule out maternal genome or POO effects as risk factors because very small effects might not be detectable with our sample size, they may influence risk through interactions with environmental exposures or may act through a more complex network of interacting genes. Thus, the most promising SNPs identified by this study may still be worth further investigation.

Evidence of gene-environment interaction for the RUNX2 gene and environmental tobacco smoke in controlling the risk of cleft lip with/without cleft palate.

This study examined the association between 49 markers in the Runt-related transcription factor 2 (RUNX2) gene and nonsyndromic cleft lip with/without cleft palate (CL/P) among 326 Chinese case-parent trios, while considering gene-environment (GxE) interaction and parent-of-origin effects. Five single-nucleotide polymorphisms (SNPs) showed significant evidence of linkage and association with CL/P and these results were replicated in an independent European sample of 825 case-parent trios. We also report compelling evidence for interaction between markers in RUNX2 and environmental tobacco smoke (ETS). Although most marginal SNP effects (i.e., ignoring maternal exposures) were not statistically significant, eight SNPs were significant when considering possible interaction with ETS when testing for gene (G) and GxE interaction simultaneously or when considering GxE alone. Independent samples from European populations showed consistent evidence of significant GxETS interaction at two SNPs (rs6904353 and rs7748231). Our results suggest genetic variation in RUNX2 may influence susceptibility to CL/P through interacting with ETS.