Evidence of gene–environment interaction for the IRF6 gene and maternal multivitamin supplementation in controlling the risk of cleft lip with/without cleft palate

Although multiple genes have been identified as genetic risk factors for isolated, non-syndromic cleft lip with/without cleft palate (CL/P), a complex and heterogeneous birth defect, interferon regulatory factor 6 gene (IRF6) is one of the best documented genetic risk factors. In this study, we tested for association between markers in IRF6 and CL/P in 326 Chinese case–parent trios, considering gene–environment interaction for two common maternal exposures, and parent-of-origin effects. CL/P case–parent trios from three sites in mainland China and Taiwan were genotyped for 22 single nucleotide polymorphisms (SNPs) in IRF6. The transmission disequilibrium test was used to test for marginal effects of individual SNPs. We used PBAT to screen the SNPs and haplotypes for gene–environment (G × E) interaction and conditional logistic regression models to quantify effect sizes for SNP–environment interaction. After Bonferroni correction, 14 SNPs showed statistically significant association with CL/P. Evidence of G × E interaction was found for both maternal exposures, multivitamin supplementation and environmental tobacco smoke (ETS). Two SNPs showed evidence of interaction with multivitamin supplementation in conditional logistic regression models (rs2076153 nominal P = 0.019, rs17015218 nominal P = 0.012). In addition, rs1044516 yielded evidence for interaction with maternal ETS (nominal P = 0.041). Haplotype analysis using PBAT also suggested interaction between SNPs in IRF6 and both multivitamin supplementation and ETS. However, no evidence for maternal genotypic effects or significant parent-of-origin effects was seen in these data. These results suggest IRF6 gene may influence risk of CL/P through interaction with multivitamin supplementation and ETS in the Chinese population.     

Analysis of candidate genes on chromosome 2 in oral cleft case-parent trios from three populations

Isolated oral clefts, including cleft lip with/without cleft palate (CL/P) and cleft palate (CP), have a complex and heterogeneous etiology. Case-parent trios from three populations were used to study genes spanning chromosome 2, where single nucleotide polymorphic (SNP) markers were analyzed individually and as haplotypes. Case-parent trios from three populations (74 from Maryland, 64 from Singapore and 95 from Taiwan) were genotyped for 962 SNPs in 104 genes on chromosome 2, including two well-recognized candidate genes: TGFA and SATB2. Individual SNPs and haplotypes (in sliding windows of 2–5 SNPs) were used to test for linkage and disequilibrium separately in CL/P and CP trios. A novel candidate gene (ZNF533) showed consistent evidence of linkage and disequilibrium in all three populations for both CL/P and CP. SNPs in key regions of ZNF533 showed considerable variability in estimated genotypic odds ratios and their significance, suggesting allelic heterogeneity. Haplotype frequencies for regions of ZNF533 were estimated and used to partition genetic variance into among-and within-population components. Wright’s fixation index, a measure of genetic diversity, showed little difference between Singapore and Taiwan compared with Maryland. The tensin-1 gene (TNS1) also showed evidence of linkage and disequilibrium among both CL/P and CP trios in all three populations, albeit at a lower level of significance. Additional genes (VAX2, GLI2, ZHFX1B on 2p; WNT6–WNT10A and COL4A3–COL4A4 on 2q) showed consistent evidence of linkage and disequilibrium only among CL/P trios in all three populations, and TGFA showed significant evidence in two of three populations.Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.     

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.     

Interaction between IRF6 and TGFA Genes Contribute to the Risk of Nonsyndromic Cleft Lip/Palate

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⁻⁶). IRF6 was also associated with cleft palate (CP) with impaction of permanent teeth (p<10⁻⁶). 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⁻⁶ 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.     

Maternal Transmission Effect of a PDGF-C SNP on Nonsyndromic Cleft Lip with or without Palate from a Chinese Population

Cleft lip with or without palate (CL/P) is a common congenital anomaly with a high birth prevalence in China. Based on a previous linkage signal of nonsyndromic CL/P (NSCL/P) on the chromosomal region 4q31–q32 from the Chinese populations, we screened the 4q31–q32 region for susceptibility genes in 214 trios of Han Chinese. PDGF-C, an important developmental factor, resides in the region and has been implicated in NSCL/P. However, in our family-based association test (transmission disequilibrium test; TDT), we could not conclude an association between PDGF-C and NSCL/P as previously suggested. Instead, we found strong evidence for parent-of-origin effect at a PDGF-C SNP, rs17035464, by a likelihood ratio test (unadjusted p-value = 0.0018; I_m = 2.46). The location of rs17035464 is 13 kb downstream of a previously reported, NSCL/P-associated SNP, rs28999109. Furthermore, a patient from our sample trios was observed with a maternal segmental uniparental isodisomy (UPD) in a region containing rs17035464. Our findings support the involvement of PDGF-C in the development of oral clefts; moreover, the UPD case report contributes to the collective knowledge of rare variants in the human genome.     

NOS2A, TLR4, and IFNGR1 interactions influence pulmonary tuberculosis susceptibility in African-Americans

Tuberculosis (TB) has substantial mortality worldwide with 5–10% of those exposed progressing to active TB disease. Studies in mice and humans indicate that the inducible nitric oxide synthase (iNOS) molecule plays an important role in immune response to TB. A mixed case–control association study of individuals with TB, relatives, or close contact controls was performed in 726 individuals (279 case and 166 control African-Americans; 198 case and 123 control Caucasians). Thirty-nine single nucleotide polymorphisms (SNPs) were selected from the NOS2A gene for single SNP, haplotype, and multilocus interaction analyses with other typed candidate genes using generalized estimating equations. In African-Americans, ten NOS2A SNPs were associated with TB. The strongest associations were observed at rs2274894 (odds ratio (OR) = 1.84, 95% confidence interval (CI) [1.23–2.77], p = 0.003) and rs7215373 (OR = 1.67, 95% CI [1.17–2.37], p = 0.004), both of which passed a false discovery rate correction for multiple comparisons (q* = 0.20). The strongest gene–gene interactions were observed between NOS2A rs2248814 and IFNGR1 rs1327474 (p = 0.0004) and NOS2A rs944722 and IFNGR1 rs1327474 (p = 0.0006). Three other SNPs in NOS2A interacted with TLR4 rs5030729 and five other NOS2A SNPs interacted with IFNGR1 rs1327474. No significant associations were observed in Caucasians. These results suggest that NOS2A variants may contribute to TB susceptibility, particularly in individuals of African descent, and may act synergistically with SNPs in TLR4 and IFNGR1.     

Novel Approach Identifies SNPs in SLC2A10 and KCNK9 with Evidence for Parent-of-Origin Effect on Body Mass Index

The phenotypic effect of some single nucleotide polymorphisms (SNPs) depends on their parental origin. We present a novel approach to detect parent-of-origin effects (POEs) in genome-wide genotype data of unrelated individuals. The method exploits increased phenotypic variance in the heterozygous genotype group relative to the homozygous groups. We applied the method to >56,000 unrelated individuals to search for POEs influencing body mass index (BMI). Six lead SNPs were carried forward for replication in five family-based studies (of ∼4,000 trios). Two SNPs replicated: the paternal rs2471083-C allele (located near the imprinted KCNK9 gene) and the paternal rs3091869-T allele (located near the SLC2A10 gene) increased BMI equally (beta = 0.11 (SD), P<0.0027) compared to the respective maternal alleles. Real-time PCR experiments of lymphoblastoid cell lines from the CEPH families showed that expression of both genes was dependent on parental origin of the SNPs alleles (P<0.01). Our scheme opens new opportunities to exploit GWAS data of unrelated individuals to identify POEs and demonstrates that they play an important role in adult obesity.     

Evidence of statistical epistasis between DISC1, CIT and NDEL1 impacting risk for schizophrenia: biological validation with functional neuroimaging

The etiology of schizophrenia likely involves genetic interactions. DISC1, a promising candidate susceptibility gene, encodes a protein which interacts with many other proteins, including CIT, NDEL1, NDE1, FEZ1 and PAFAH1B1, some of which also have been associated with psychosis. We tested for epistasis between these genes in a schizophrenia case–control study using machine learning algorithms (MLAs: random forest, generalized boosted regression and Monte Carlo logic regression). Convergence of MLAs revealed a subset of seven SNPs that were subjected to 2-SNP interaction modeling using likelihood ratio tests for nested unconditional logistic regression models. Of the ⁷C₂ = 21 interactions, four were significant at the α = 0.05 level: DISC1 rs1411771–CIT rs10744743 OR = 3.07 (1.37, 6.98) p = 0.007; CIT rs3847960–CIT rs203332 OR = 2.90 (1.45, 5.79) p = 0.003; CIT rs3847960–CIT rs440299 OR = 2.16 (1.04, 4.46) p = 0.038; one survived Bonferroni correction (NDEL1 rs4791707–CIT rs10744743 OR = 4.44 (2.22, 8.88) p = 0.00013). Three of four interactions were validated via functional magnetic resonance imaging (fMRI) in an independent sample of healthy controls; risk associated alleles at both SNPs predicted prefrontal cortical inefficiency during the N-back task, a schizophrenia-linked intermediate biological phenotype: rs3847960–rs440299; rs1411771–rs10744743, rs4791707–rs10744743 (SPM5 p < 0.05, corrected), although we were unable to statistically replicate the interactions in other clinical samples. Interestingly, the CIT SNPs are proximal to exons that encode the DISC1 interaction domain. In addition, the 3′ UTR DISC1 rs1411771 is predicted to be an exonic splicing enhancer and the NDEL1 SNP is ~3,000 bp from the exon encoding the region of NDEL1 that interacts with the DISC1 protein, giving a plausible biological basis for epistasis signals validated by fMRI.     

Joint Testing of Genotypic and Gene-Environment Interaction Identified Novel Association for BMP4 with Non-Syndromic CL/P in an Asian Population Using Data from an International Cleft Consortium

Background

Non-syndromic cleft lip with or without cleft palate (NSCL/P) is a common disorder with complex etiology. The Bone Morphogenetic Protein 4 gene (BMP4) has been considered a prime candidate gene with evidence accumulated from animal experimental studies, human linkage studies, as well as candidate gene association studies. The aim of the current study is to test for linkage and association between BMP4 and NSCL/P that could be missed in genome-wide association studies (GWAS) when genotypic (G) main effects alone were considered.

Methodology/Principal Findings

We performed the analysis considering G and interactions with multiple maternal environmental exposures using additive conditional logistic regression models in 895 Asian and 681 European complete NSCL/P trios. Single nucleotide polymorphisms (SNPs) that passed the quality control criteria among 122 genotyped and 25 imputed single nucleotide variants in and around the gene were used in analysis. Selected maternal environmental exposures during 3 months prior to and through the first trimester of pregnancy included any personal tobacco smoking, any environmental tobacco smoke in home, work place or any nearby places, any alcohol consumption and any use of multivitamin supplements. A novel significant association held for rs7156227 among Asian NSCL/P and non-syndromic cleft lip and palate (NSCLP) trios after Bonferroni correction which was not seen when G main effects alone were considered in either allelic or genotypic transmission disequilibrium tests. Odds ratios for carrying one copy of the minor allele without maternal exposure to any of the four environmental exposures were 0.58 (95%CI = 0.44, 0.75) and 0.54 (95%CI = 0.40, 0.73) for Asian NSCL/P and NSCLP trios, respectively. The Bonferroni P values corrected for the total number of 117 tested SNPs were 0.0051 (asymptotic P = 4.39*10⁻⁵) and 0.0065 (asymptotic P = 5.54*10⁻⁵), accordingly. In European trios, no significant association was seen for any SNPs after Bonferroni corrections for the total number of 120 tested SNPs.

Conclusions/Significance

Our findings add evidence from GWAS to support the role of BMP4 in susceptibility to NSCL/P originally identified in linkage and candidate gene association studies.     

Expression of “brown-in-white” adipocyte biomarkers shows gender differences and the influence of early dietary exposure

Fish oil supplementation provides an inconsistent degree of protection from cardiovascular disease (CVD), which may be attributed to genetic variation. Single nucleotide polymorphisms (SNPs) in the elongation-of-very-long-chain-fatty-acids-2 (ELOVL2) gene have been strongly associated with plasma proportions of n-3 long-chain polyunsaturated fatty acids (LC-PUFA). We investigated the effect of genotype interaction with fish oil dosage on plasma n-3 LC-PUFA proportions in a parallel double-blind controlled trial, involving 367 subjects randomised to treatment with 0.45, 0.9 and 1.8 g/day eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA) (1.51:1) or olive oil placebo for 6 months. We genotyped 310 subjects for ELOVL2 gene SNPs rs3734398, rs2236212 and rs953413. At baseline, carriers of all minor alleles had lower proportions of plasma DHA than non-carriers (P = 0.021–0.030). Interaction between genotype and treatment was a significant determinant of plasma EPA (P < 0.0001) and DHA (P = 0.004–0.032). After the 1.8 g/day dose, carriers of ELOVL2 SNP minor alleles had approximately 30 % higher proportions of EPA (P = 0.002–0.004) and 9 % higher DHA (P = 0.013–0.017) than non-carriers. Minor allele carriers could therefore particularly benefit from a high intake of EPA and DHA in maintaining high levels of plasma n-3 PUFA conducive to protection from CVD.     

Comparison of multimarker logistic regression models, with application to a genomewide scan of schizophrenia

Background

Imprinted genes show expression from one parental allele only and are important for development and behaviour. This extreme mode of allelic imbalance has been described for approximately 56 human genes. Imprinting status is often disrupted in cancer and dysmorphic syndromes. More subtle variation of gene expression, that is not parent-of-origin specific, termed 'allele-specific gene expression' (ASE) is more common and may give rise to milder phenotypic differences. Using two allele-specific high-throughput technologies alongside bioinformatics predictions, normal term human placenta was screened to find new imprinted genes and to ascertain the extent of ASE in this tissue.

Results

Twenty-three family trios of placental cDNA, placental genomic DNA (gDNA) and gDNA from both parents were tested for 130 candidate genes with the Sequenom MassArray system. Six genes were found differentially expressed but none imprinted. The Illumina ASE BeadArray platform was then used to test 1536 SNPs in 932 genes. The array was enriched for the human orthologues of 124 mouse candidate genes from bioinformatics predictions and 10 human candidate imprinted genes from EST database mining. After quality control pruning, a total of 261 informative SNPs (214 genes) remained for analysis. Imprinting with maternal expression was demonstrated for the lymphocyte imprinted gene ZNF331 in human placenta. Two potential differentially methylated regions (DMRs) were found in the vicinity of ZNF331. None of the bioinformatically predicted candidates tested showed imprinting except for a skewed allelic expression in a parent-specific manner observed for PHACTR2, a neighbour of the imprinted PLAGL1 gene. ASE was detected for two or more individuals in 39 candidate genes (18%).

Conclusions

Both Sequenom and Illumina assays were sensitive enough to study imprinting and strong allelic bias. Previous bioinformatics approaches were not predictive of new imprinted genes in the human term placenta. ZNF331 is imprinted in human term placenta and might be a new ubiquitously imprinted gene, part of a primate-specific locus. Demonstration of partial imprinting of PHACTR2 calls for re-evaluation of the allelic pattern of expression for the PHACTR2-PLAGL1 locus. ASE was common in human term placenta.     

Genetic variants in the JAK1 gene confer higher risk of Behcet’s disease with ocular involvement in Han Chinese

Recent surveys have identified SLC22A4, SLC22A5, RUNX1, JAK1 as susceptibility genes for various immune-related diseases. An association study was performed in 738 Behcet’s patients with ocular involvement and 1,873 controls using the iPLEX system method. The first-stage study for 30 SNPs showed that SNPs rs2780815, rs310241, rs3790532 in JAK1 were associated with Behcet’s disease in Han Chinese (Pc_{(Bonferroni correction)} = 0.022–7.7 × 10^?3). The G allele and AA genotype of SNP rs2834643 in RUNX1 (Pc = 0.041–1.75 × 10^?3), but none of the other SNPs, were associated with Behcet’s disease. Haplotype analysis for the SLC22A4, SLC22A5 genes showed an increased tendency for AGTCTGCCGC frequency in patients compared with controls; however, the significance was lost after Bonferroni correction (P = 0.004, Pc > 0.05). Subsequently, we further replicated the significantly associated SNPs using another independent cohort. Replication and combining studies showed that three SNPs rs2780815, rs310241, rs3790532 in JAK1, but not SNP rs2834643 in RUNX1, were consistently associated with Behcet’s disease (replication: Pc = 0.012–9.60 × 10^?4; combining: Pc = 0.030–1.90 × 10^?4). SNPs rs2780815, rs310241, rs3790532 were estimated to confer a population attributable risk of 35.0, 28.0, 27.0 %, respectively. We found a strong association between HLA-B51 with Behcet’s disease in Chinese Han population (P = 1.35 × 10^?73; OR = 5.15; 95 % CI 4.28–6.19). GMDR analysis showed that no gene–gene interaction was detectable between JAK1 and HLA-B51. Logistic analysis indicated that the JAK1 gene was an independent risk factor for Behcet’s disease (P > 0.05). Real-time PCR analysis showed that no difference on the expression of JAK1 in PBMCs or LPS-stimulated PBMCs between individuals with the different rs1762780815 genotypes studied (P > 0.05). In conclusion, this study suggests that JAK1, but not SLC22A4, SLC22A5 and RUNX1, contributes to the genetic susceptibility to Behcet’s disease with ocular involvement.     

The novel human MRC1 gene polymorphisms are associated with susceptibility to pulmonary tuberculosis in Chinese Uygur and Kazak populations

The MRC1 gene, encoding the human mannose receptor (MR), is a member of the C-type lectin receptors family. MR can recognize and bind to Mycobacterium tuberculosis by the extracellular structure, and play a role in antigen-presenting and maintaining a stable internal environment. This study aimed to investigate potential associations of SNPs in exon 7 of the MRC1 gene with pulmonary tuberculosis (TB). G1186A, G1195A, T1212C, C1221G, C1303T and C1323T were genotyped using PCR and DNA sequencing in 595 Chinese Uygur and 513 Kazak subjects. In the Uygur, the frequency of allele G (P = 0.031, OR = 1.29, 95 % CI = 1.02–1.62) and AA genotype (P = 0.033, OR = 1.64, 95 % CI = 1.04–2.60) for G1186A was lower in the pulmonary TB than healthy control and were significantly correlated with pulmonary TB. After adjustment for age and gender, G1186A was found to be additive models in association with pulmonary TB (P = 0.04, OR = 1.27, 95 % CI = 1.01–1.60). By calculating linkage disequilibrium, the frequency of haplotype GGTCCT (P = 0.032, OR = 0.75, 95 % CI = 0.57–0.97) and GGTCCC (P = 0.044, OR = 0.57, 95 % CI = 0.33–0.99) was significantly associated with pulmonary TB. No association was found between other SNPs and pulmonary TB. In the Kazak, all SNPs were not associated with pulmonary TB. Our results suggest that genetic factors play an important role in susceptibility to pulmonary TB at the individual level, and provide an experimental basis to clarify the pathogenesis of pulmonary TB.     

An extension of the transmission disequilibrium test incorporating imprinting

The recombination rates in meioses of females and males are often different. Some genes that affect development and behavior in mammals are known to be imprinted, and >1% of all mammalian genes are believed to be imprinted. When the gene is imprinted and the recombination fractions are sex specific, the conventional transmission disequilibrium test (TDT) is shown to be still valid for testing for linkage. The power function of the TDT is derived, and the effect of the degree of imprinting on the power of the TDT is investigated. It is learned that imprinting has little effect on the power when the female and male recombination rates are equal. On the basis of case-parents trios, the transmissions from the heterozygous fathers/mothers to their affected children are separated as paternal and maternal, and two TDT-like statistics, TDT(p) and TDT(m), are consequently constructed. It is found that the TDT(p) possesses a higher power than the TDT for maternal imprinting genes, and the TDT(m) is more powerful than the TDT for paternal imprinting genes. On the basis of the parent-of-origin effects test statistic (POET), a novel statistic, TDT incorporating imprinting (TDTI) is proposed to test for linkage in the presence of linkage disequilibrium, which is shown to be more powerful than the TDT when parent-of-origin effects are significant but slightly less powerful than the TDT when parent-of-origin effects are negligible. The validity of the TDT and TDTI is assessed by simulation. The power approximation formulas for the TDT and TDTI are derived and the simulation results show that they are accurate. The simulation study on power comparison shows that the TDTI outperforms the TDT for imprinted genes. The improvement can be substantial in the case of complete paternal/maternal imprinting.     

PCMT1 gene polymorphisms, maternal folate metabolism, and neural tube defects: a case–control study in a population with relatively low folate intake

The PCMT1 gene encodes the protein repair enzyme protein-l-isoaspartate (d-aspartate) O-methyltransferase, which is known to protect certain neural cells against Bax-induced apoptosis. Previous studies have produced inconsistent results regarding the effects of PCMT1 (rs4816 and rs4552) polymorphisms on neural tube defects (NTDs). Reduced maternal plasma folate levels and/or elevated homocysteine (Hcy) levels are considered to be risk factors for NTDs. In order to clarify the key factors contributing to the apparent discrepancy and investigate gene–environment interaction, we conducted a case–control study including 121 cases and 146 matched controls to investigate the association between the two PCMT1 polymorphisms in fetuses and the risk of NTDs in the Chinese population of Lvliang, which has low folate intake. Maternal plasma folate and Hcy levels were also measured, and the interaction between fetal PCMT1 gene status and maternal folate metabolites was assessed. Maternal plasma folate concentrations in the NTD group were lower than in controls (10.23 vs. 13.08 nmol/L, adjusted P = 0.059), and Hcy concentrations were significantly higher (14.46 vs. 11.65 μmol/L, adjusted P = 0.026). Fetuses carrying the rs4816 AG + GG genotype, combined with higher maternal plasma Hcy, had a 6.46-fold (95 % CI 1.15–36.46) increased risk of anencephaly. The results of this study imply that the fetal PCMT1 rs4816 polymorphism may play only a weak role in NTD formation and that gene–environment interactions might be more significant.     

TGFA and IRF6 Contribute to the Risk of Nonsyndromic Cleft Lip with or without Cleft Palate in Northeast China

Nonsyndromic cleft lip with or without cleft palate (NSCL/P) are common birth defects with a complex etiology. Multiple interacting loci and possible environmental factors influence the risk of NSCL/P. 12 single nucleotide polymorphisms (SNPs) in 7 candidate genes were tested using an allele-specific primer extension for case-control and case-parent analyses in northeast China in 236 unrelated patients, 185 mothers and 154 fathers, including 128 complete trios, and 400 control individuals. TGFA and IRF6 genes showed a significant association with NSCL/P. In IRF6, statistical evidence of an association between rs2235371 (p = 0.003), rs2013162 (p<0.0001) and NSCL/P was observed in case-control analyses. Family based association tests (FBATs) showed over-transmission of the C allele at the rs2235371 polymorphism (p = 0.007). In TGFA, associations between rs3771494, rs3771523 (G3822A), rs11466285 (T3851C) and NSCL/P were observed in case-control and FBAT analyses. Associations between other genes (BCL3, TGFB3, MTHFR, PVRL1 and SUMO1) and NSCL/P were not detected.     

Limited Evidence for Parent-of-Origin Effects in Inflammatory Bowel Disease Associated Loci

Background

Genome-wide association studies of two main forms of inflammatory bowel diseases (IBD), Crohn’s disease (CD) and ulcerative colitis (UC), have identified 99 susceptibility loci, but these explain only ∼23% of the genetic risk. Part of the ‘hidden heritability’ could be in transmissible genetic effects in which mRNA expression in the offspring depends on the parental origin of the allele (genomic imprinting), since children whose mothers have CD are more often affected than children with affected fathers. We analyzed parent-of-origin (POO) effects in Dutch and Indian cohorts of IBD patients.

Methods

We selected 28 genetic loci associated with both CD and UC, and tested them for POO effects in 181 Dutch IBD case-parent trios. Three susceptibility variants in NOD2 were tested in 111 CD trios and a significant finding was re-evaluated in 598 German trios. The UC-associated gene, BTNL2, reportedly imprinted, was tested in 70 Dutch UC trios. Finally, we used 62 independent Indian UC trios to test POO effects of five established Indian UC risk loci.

Results

We identified POO effects for NOD2 (L1007fs; OR = 21.0, P-value = 0.013) for CD; these results could not be replicated in an independent cohort (OR = 0.97, P-value = 0.95). A POO effect in IBD was observed for IL12B (OR = 3.2, P-value = 0.019) and PRDM1 (OR = 5.6, P-value = 0.04). In the Indian trios the IL10 locus showed a POO effect (OR = 0.2, P-value = 0.03).

Conclusions

Little is known about the effect of genomic imprinting in complex diseases such as IBD. We present limited evidence for POO effects for the tested IBD loci. POO effects explain part of the hidden heritability for complex genetic diseases but need to be investigated further.     

Replication and extension of association of choline acetyltransferase with nicotine dependence in European and African American smokers

Choline acetyltransferase is critical in the synthesis of acetylcholine and regulation of cholinergic neuron functions. We recently reported association of the encoding gene ChAT with both smoking cessation and nicotine dependence (ND) in two independent European American (EA) samples; however, in the replication sample, only limited SNPs partially covering the gene were examined. In this study, we examined the association of 14 SNPs, which cover the entire gene, with ND, assessed by smoking quantity (SQ), heaviness of smoking index (HSI), and Fagerström Test for ND (FTND), in 2,037 subjects from 602 families of African American (AA) or EA origin. Individual SNP-based association analysis revealed that five SNPs showed nominal association with at least one ND measure in one of the samples (P = 0.022–0.042); none remained significant after correction for multiple testing. Haplotype-based association analysis revealed that haplotypes G–G–A–C, formed by rs1880676–rs3810950–rs10082479–rs8178990 (P = 0.005–0.0178), and G–G–T–C–G–C, formed by rs1880676–rs3810950–rs10082479–rs8178990–rs3793790–rs12266458 (P = 0.00247–0.00468), displayed significant association with all three ND measures in the AA sample, as did haplotype T–C–G–A–T, formed by rs12266458–rs11101191–rs8178991–rs4838544–rs4838547 (P = 0.00741–0.0103), in the EA sample. All these detected haplotype-based associations remained significant after correction for all major haplotypes for a given SNP combination. Together, our findings, in conjunction with the previous report of the association, warrant further investigation of ChAT in ND.