Abstract
Small ubiquitin-like modifier 1 (SUMO1) and environmental factors have been shown to be associated with nonsyndromic cleft lip with or without cleft palate (NSCL/P) in several populations. This study aimed at confirming the contribution of SUMO1 gene and environmental factors to nonsyndromic orofacial clefts risk in western Han Chinese. Four single-nucleotide polymorphisms were investigated in 212 case trios in western China using conditional logistic regression models and the transmission disequilibrium test under a case–parent trio design. Strong evidence of linkage and linkage disequilibrium was found between these markers and the disease in both single-nucleotide polymorphism analysis (T allele at rs6761234 [p = 0.0005, odds ratio [OR] = 1.82, 95% confidence interval [CI]:1.30–2.57) and C allele at rs12470401 (p < 0.0001, OR = 2.82, 95% CI: 1.90–4.19)] and sliding window haplotype analysis (T-T-T for rs6761234-rs12470401-rs7599810 [p = 0.018], C-C-G for rs12470401-rs7599810-rs6435133 [0.0033], C-T-T-T for rs6761234-rs12470401-rs7599810-rs6435133 [p = 0.018] among others). Interactions between mothers' passive smoking during the first trimester and C/C genotype of rs12470401 showed statistical significance (OR0 = 2.53 and OR1 = 8.83). Risk factors identified in this study may provide a better understanding of the etiological role of SUMO1 gene in NSCL/P incidence.
Introduction
Small ubiquitin-like modifier 1 (SUMO1) gene, which is located on 2q33, encodes a protein that is a member of the SUMO protein family. It functions in a manner similar to ubiquitin in that it is bound to target proteins as part of a post-translational modification system (NCBI).
We considered SUMO1 a potential candidate (gene), because it was abundantly expressed on mouse embryonic day 13.5 in the upper lip, primary palate, and medial edge epithelia of the secondary palate. When the SUMO1 gene is disrupted by balanced translocation, it may lead to cleft lip and palate (Alkuraya et al., 2006). Moreover, case–control and family-based association studies have identified linkage disequililibrum (LD) between polymorphic alleles at the SUMO1 locus and NSCL/P (Song et al., 2008; Carter et al., 2010).
All these reports suggest that SUMO1 may predispose to CL/P. The risk of CL/P associated with particular variants in SUMO1 may differ among different ethnic groups, and it is rarely reported to be associated with NSCL/P in the western Han Chinese.
It is known that NSCL/P is a complex and heterogeneous congenital malformation; its occurrence is highly associated with several environmental and genetic risk factors. Many epidemiological investigations and animal experiments have confirmed the exposure to environmental risk factors during early pregnancy, such as smoking, drinking, lack of vitamins, viral infection, and mothers' illness and medication during the first trimester can increase the risk of NSCL/P (Little et al., 2004; Shaw et al., 2005; Nandor et al., 2006; Grewal et al., 2008; Jia et al., 2009). In this study, transmission disequilibrium test (TDT) and conditional logistic regression models were used in case–parent trio design to find the relationship among environmental exposure factors during the first trimester and four single-nucleotide polymorphisms (SNPs) in SUMO1 gene in 212 case–parent trios from western China.
Materials and Methods
Subjects
Our sample consists of 212 complete case trios, all of whom were recruited between 2008 and 2009 from the Cleft Surgery Department of West China College of Stomatology, Sichuan University. Informed consent was obtained from each participant before enrollment in the study, as approved by the local ethics committee. To assess nonsyndromic status of cases, all probands were screened for the presence of associated anomalies or syndromes by a physician, and only those determined to have an isolated CL/P were included in this study. Patients were asked about family history of oral clefts among first and second relatives. All subjects were self-identified as western Han Chinese. Table 1 shows gender and family history among NSCL/P cases.
NSCL/P, nonsyndromic cleft lip with or without cleft palate.
Data collection
A brief interview of mothers and fathers was conducted at their initial visit, including family history, fathers' drinking and smoking history, mothers' abortion history, and environmental exposures during the first trimester: medication history (if they took any medication against viral infection), illness history, vitamins (including folic acid), supplementation history, and smoking history. Mothers were asked whether they or their partner smoked at the beginning of pregnancy, whether they currently smoked, and about daily and/or weekly cigarette consumption. A separate question was asked about the total number of cigarettes smoked indoors. Smoking was coded as a dichotomous variable: If they smoked more than one cigarette within a week, they were coded as smokers; if the answer was “no,” they were coded as nonsmokers; and if the answers to all the questions on smoking were missing, they were coded as missing.
Genotyping
Venous blood samples were drawn from participants at their initial visit. DNA was extracted by phenol chloroform extraction protocol and genotyped for SNP markers using the restriction fragment length polymorphism–polymerase chain reaction (Supplementary Table S1; Supplementary Data are available online at
Statistical analysis
Hardy–Weinberg equilibrium at each SNP was computed among parents. The TDT were performed by the family-based association test (FBAT) program (
Results
Hardy–Weinberg equilibrium and LD
Hardy–Weinberg equilibrium was assessed for each of the four SNPs among the parents of cases. Significant deviation from Hardy–Weinberg expectations could reflect genotyping errors or true heterogeneity in the population and could bias our statistical tests. All χ 2 tests using the genotype frequency of parents showed that the four SNPs included in the study were in Hardy–Weinberg equilibrium. LD between markers was calculated by D′ and r 2 statistics from haploview analysis. Among the four SNPs, one group of markers (rs6761234-rs7599810) showed tight LD (D′ = 0.74); whereas the other groups of markers showed relatively weak LD (D′ < 0.8 and r 2 < 0.8) (Supplementary Table S2).
TDT analyses of individual markers and haplotypes
TDT results of the 2 × 2 χ 2 test carried out on heterozygous case–parent trios showed that the T allele at rs6761234 (p = 0.0005, OR = 1.82, 95% confidence interval (CI): 1.30–2.57) and C allele at rs12470401 (p < 0.0001, OR = 2.82, 95% CI: 1.90–4.19) were overtransmitted from parents to affected progeny. This suggests that these alleles may significantly increase the risk of being NSCL/P cases (Table 2).
95% CI: 95% confidence interval; OR, odds ratio.
<0.01.
<0.0001.
SNP, single-nucleotide polymorphism; T/NT, transmitted or nontransmitted counts from heterozygous parents.
Bold numbers are statistically significant.
Results of the conditional logistic regress analysis showed that the T/T genotype at rs6435133 may decrease the risk of being NSCL/P cases (p = 0.043, OR = 0.055, 95% CI: 0.022–0.68) (Table 3).
Pseudo control.
<0.05.
Bold numbers are statistically significant.
FBAT analysis showed that the C allele (p = 0.0038), C/C genotype (p = 0.00039) at rs12470401, and G/G genotype at rs6435133 (p = 0.023) were over-represented in the affected progeny. FBAT analysis also showed that the T/C genotype at rs12470401 was under-represented among affected progeny, which suggests that this genotype may decrease the risk of being affected (Table 4).
0.005.
<0.01.
S, score; E(S), expected score; Var (S), variance (score).
Bold numbers are statistically significant.
In Table 5, haplotypes consisting of two SNPs (rs6761234-rs12470401 and rs7599810-rs6435133) yielded statistical significance among NSCL/P cases (p = 0.03 and 0.031); haplotypes consisting of three SNPs (rs6761234-rs12470401-rs7599810 and rs12470401-rs7599810-rs6435133) also yielded statistical significance among NSCL/P cases (p = 0.018 and 0.0033); and haplotypes consisting of four SNPs yielded statistical significance among NSCL/P cases (p = 0.018). Interestingly then, two-SNP haplotypes (e.g., T-C for rs6761234-rs12470401 or C-G composed of rs7599810 and rs6435133) were more informative than three- or four-SNP haplotypes. Haplotype analysis showed that some haplotypes, for example, T-C (for rs6761234-rs12470401), C-G (for rs7599810-rs6435133) were far more common, and these haplotypes had a stronger impact on risk to NSCL/P.
Number of informative case–parent trios.
The p-value for a statistical model with whole haplotypes consisted of two, three, or four SNPs.
Interactions between SUMO1 polymorphism and environmental factors
Jia et al. (2010) used single-factor χ 2 analysis to test the association between nonsyndromic orofacial clefts and environmental factors. Their results showed that mothers'smoking, passive smoking, drinking, illness, medication, abortion history, and vitamins (including folic acid) supplementation during the first trimester were significantly associated with nonsyndromic orofacial clefts in western Han Chinese.
In this study, the previous “significant” environmental factors (mothers'smoking, passive smoking, drinking, illness, medication, abortion history, vitamins [including folic acid] supplementation during the first trimester), and genotypic risk factors (C/C genotype at rs12470401) were used to test for interactions between markers in SUMO1 and environmental factors from case–parental trio designs. The interaction between infant's C/C genotype at rs12470401 and mothers' drinking during the first trimester was significant (OR0 = 2.53 and OR1 = 8.83) (Table 6).
Case–parental control.
0, T/T genotype of rs12470401; 1, C/C and T/C genotype of rs12470401.
Discussion
NSCL/P is a complex and heterogenerous trait with no obvious mode of inheritance. Numerous studies have failed to identify genes with any major influence on risk to NSCL/P (Jugessur and Murray et al., 2005). The SUMO1 gene was analyzed in the present study, because it was expressed on mouse embryonic day 13.5 in the upper lip, primary palate, and medial edge epithelia of the secondary palate; additional studies on animal models confirmed the role of SUMO1 in palate formation (Alkuraya et al., 2006).
The four SNPs examined here are located in the intron of the SUMO1 gene. Although intronic SNPs do not typically alter protein structure, associations with intronic variants have been reported for a number of complex diseases. Among these four SNPs in our entire sample, C allele and C/C genotype at rs12470401, T allele at rs6761234, and G/G genotype at rs6435133 were overtransmitted from parents to affected progeny with a significant increase being NSCL/P cases.
In the studied sample of case–parent trios, haplotypes carrying the G allele at rs6435133 and C allele at rs12470401 were found to be overtransmitted from parents to affected progeny, whereas those haplotypes carrying the T allele at rs12470401 were negatively associated. Here, the alleles were reported on the forward strand of the chromosome (NCBI build 36.3), although the gene is transcribed from the reverse strand. Allele designations need to be reversed when compared with published reports that used the transcription strand as the reference.
Many studies identified a relationship between environmental factors and NSCL/P (Little et al., 2004; Johnson and Little, 2008; Jia et al., 2009, 2010). In this study, statistically significant evidence was found from case–parental trio design for the gene- environment interaction between C/C genotype at rs12470401 and mothers' drinking during the first trimester (OR0 = 2.53 and OR1 = 8.83) (Table 6), which provided new estimate of the possible impact of environmental factors and SUMO1 polymorphism on oral clefts. This argues that genes such as SUMO1 may become more critical in the etiology of NSCL/P in the presence of certain maternal exposures.
Footnotes
Acknowledgments
This research was supported by the Key Program of National Natural Science Foundation of China (Grant No. 30530730), Scientific Research Funds for young teachers of Sichuan University (No. 2008067), and Sichuan University (No. 2009SCU11159).
Disclosure Statement
No competing financial interests exist.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
