PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nihpaAbout Author manuscriptsSubmit a manuscriptNIH Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
 
Twin Res Hum Genet. Author manuscript; available in PMC Jun 29, 2010.
Published in final edited form as:
PMCID: PMC2893889
NIHMSID: NIHMS153771
SEARCH FOR GENOMIC ALTERATIONS IN MONOZYGOTIC TWINS DISCORDANT FOR CLEFT LIP AND/OR PALATE
Jane W. Kimani,1 Koh-ichiro Yoshiura,2 Min Shi,3 Astanand Jugessur,4 Danilo Moretti-Ferreira,5 Kaare Christensen,6 and Jeffrey C. Murray1
1Department of Pediatrics, University of Iowa, Iowa City, IA 52242, USA. Present address: Department of Pathology and Lab Medicine, University of North Carolina, 101 Manning Dr., Chapel Hill, NC 27514
2Department of Human Genetics, Nagasaki University Graduate School of Biomedical sciences, Nagasaki, Japan
3Biostatistics Branch, National Institute of Environmental Health Sciences (NIEHS), Research Triangle Park, NC 27709
4Craniofacial Development, Murdoch Children's Research Institute, Royal Children's Hospital, Flemington Rd, Parkville 3052, Australia
5Servico de Aconselhamento, Genetico da Universidade Estadual Paulista, Caixa Postal 529 Botucatu S.P., Brazil 18618-000
6Department of Epidemiology, Institute of Public Health, University of Southern Denmark, J.B. Winslows Vej 9, DK-5000 Odense C, Denmark
Corresponding author: Jeffrey C. Murray, MD, University of Iowa, Department of Pediatrics, S. Grand Avenue, 2182 ML, Iowa City IA 52242 USA, Phone: 1 319 335 6897 Fax: 1 319 335 6970, jeff-murray/at/uiowa.edu
Phenotypically discordant monozygotic twins offer the possibility of gene discovery through delineation of molecular abnormalities in one member of the twin pair. One proposed mechanism of discordance is postzygotically occurring genomic alterations resulting from mitotic recombination and other somatic changes. Detection of altered genomic fragments can reveal candidate gene loci that can be verified through additional analyses. We investigated this hypothesis using array comparative genomic hybridization; the 50K and 250K Affymetrix GeneChip® SNP arrays and an Illumina custom array consisting of 1,536 SNPs, to scan for genomic alterations in a sample of monozygotic twin pairs with discordant cleft lip and/or palate phenotypes. Paired analysis for deletions, amplifications and loss of heterozygosity, along with sequence verification of SNPs with discordant genotype calls did not reveal any genomic discordance between twin pairs in lymphocyte DNA samples. Our results demonstrate that postzygotic genomic alterations are not a common cause of monozygotic twin discordance for isolated cleft lip and/or palate. However, rare or balanced genomic alterations, tissue-specific events and small aberrations beyond the detection level of our experimental approach cannot be ruled out. The stability of genomes we observed in our study samples also suggests that detection of discordant events in other monozygotic twin pairs would be remarkable and of potential disease significance.
Keywords: Monozygotic twins, discordant, cleft lip and palate, genome-wide
Classical twin research studies comparing disease concordance rates between monozygotic (MZ) and dizygotic (DZ) twins have been extensively applied to estimate the contribution of genetic and environmental factors to many complex traits (Boomsma et al., 2002). MZ twin concordance for common diseases and traits rarely reaches 100%, an observation that is often attributed to differential environmental exposures (Wong et al., 2005). However, there are several molecular mechanisms that could underlie phenotypic discordance between MZ twins, such as de novo somatic mutations (Kondo et al., 2002), chromosomal anomalies (Gilbert et al., 2002), skewed X chromosome inactivation (De Gregorio et al., 2005), imprinting defects (Weksberg et al., 2002) and differential gene expression (Mak et al., 2004). The identification of molecular genetic differences between discordant MZ twins suggests that the utility of twin studies could be extended beyond heritability studies to gene discovery.
Postzygotic mitotic recombination and other somatic events such as deletions, nondisjunction, gene conversion, mobile genetic elements, fragile sites and repeat expansions have been proposed as a cause of MZ twin discordance (Cote and Gyftodimou, 1991; Kastern and Kryspin-Sorensen, 1988). Occurrence of somatic genomic changes in the early embryo can have significant implications for the phenotypes of the twins and the twinning process. If a genomic alteration occurs prior to twinning, unequal allocation of cells with the lesion to the two embryos could account for discordant phenotypes, with evidence of mosaicism, while the segregation of two genetically different populations of cells may trigger the twinning process, resulting in co-twins with discordant phenotypes (Machin, 1996). Since detection of chromosomal abnormalities relies on the comparison of a test genome with a reference genome sequence, the identical genomes of MZ twins ensure that any genomic differences can be established with confidence and accuracy. This approach is particularly advantageous for delineating disease causing copy number variants that could arise in one member of the twin pair. Any genomic mismatches between the twins’ genomes could be revealed as amplifications, deletions or loss of heterozygosity. Recent identification of novel cleft lip and/or palate (CLP) candidate genes by genome-wide array-CGH highlights the impact of genomic alterations in the etiology of this complex disorder (Osoegawa et al., 2008).
CLP occurs in approximately 1 out of 700 births worldwide and has a significant clinical and economic impact. Approximately 40% of MZ twins are discordant for the disease phenotype (Christensen and Fogh-Andersen, 1993a). Genetic factors contributing to disease etiology include IRF6 (Zucchero et al., 2004); MSX1 (Jezewski et al., 2003); RYK (Watanabe et al., 2006) and genes in the FGF pathway (Riley et al., 2007). Additionally, non-genetic factors such as maternal cigarette smoking and nutrition also play a role in CLP etiology (Jugessur and Murray, 2005).
The complex nature of CLP complicates traditional mapping approaches such as linkage and association, which rely on population-matched or family-based controls to determine if a sequence variant is disease-causing. Based on identical genetic backgrounds, the unaffected twin in a discordant pair provides a well-matched control for studying a complex disease such as CLP thereby eliminating the need for external controls. We have previously reported our model using MZ twins discordant for CLP in gene discovery (Mansilla et al., 2005). In this report, we describe our search for genomic alterations in MZ twins discordant for CLP using array CGH and high-density SNP genotyping arrays.
Study subjects
A summary of the discordant MZ twin samples and the respective analyses employed for each twin pair are presented in table 1. The variability in the methods used to analyze the different samples resulted from multiple stages of the study being conducted over a period of time. All affected individuals were nonsyndromic cases of CLP. We used DNA previously extracted from peripheral blood lymphocytes. Genotyping a set of DNA markers had previously established that the twins were MZ. Informed consent was obtained for all study participants and the institutional review board approved the study.
Table 1
Table 1
Twin samples analyzed in this study
Array CGH
CGH is a molecular cytogenetic technique for analyzing DNA copy number variations. DNA from a test and reference sample are differentially labeled and hybridized to an array spotted with a genomic representation that allows detection of copy number differences between the two samples at specific genomic locations (Pinkel and Albertson, 2005). For our analysis, Koh-ichiro Yoshiura et al. at Nagasaki University developed an array spotted with 2,173 genomic BAC clones. Samples from six pairs of twins were analyzed: for each pair, the DNA samples were differentially labeled, hybridized, and then scanned. Reverse sample labeling was also performed. The average normalized inter-locus fluorescence ratio (ANILFR) between the affected and unaffected twin samples was calculated. Prior preliminary experiments using five sets of normal/normal control samples had established the thresholds for copy number gain and loss at 0.86 and 1.18 with a standard deviation (SD) of 0.06, so the normal ANILFR range was defined as within ± 2 SD.
SNP genotyping arrays
Large scale SNP genotyping allows detection of allelic imbalances such as loss of heterozygosity (LOH) and copy number changes through hybridization signal intensities. In our analysis, we utilized a custom BeadArray™ platform from Illumina, Inc., (San Diego, CA, USA) and Affymetrix GeneChip® Human Mapping 50K and 250K sets (Santa Clara, CA, USA).
Illumina custom SNP genotyping
The Illumina BeadChip is a platform for performing multiplex gene analyses using oligonucleotides attached to silica beads. Defined SNPs can be chosen for a custom genotyping array based on the interests of the researcher (Steemers and Gunderson, 2007). For our analysis, we selected 1,536 SNPs representing 388 CLP candidate genes for a single genotyping array. The experiments were carried out at the Center for Inherited Disease Research (CIDR) according to the manufacturer’s instructions and included 20 pairs of discordant MZ twins. Data were analyzed using Illumina’s BeadStudio v2 genotyping and LOH plus module software that allow detection of chromosomal aberrations and allelic imbalance in paired samples (www.illumina.com).
Affymetrix SNP genotyping
Affymetrix GeneChips employ a whole genome sampling analysis method to genotype thousands of SNPs on synthetic oligonucleotide arrays by allele-specific hybridization. DNA samples were prepared for analysis and hybridization according to the manufacturer’s instructions (www.affymetrix.com). We used the GeneChip® Human Mapping 50K Xba I (n=2 twin pairs) and the 250K Nsp I (n=10 twin pairs) arrays. SNP allele calls were assigned using a Dynamic Model mapping algorithm, a highly accurate genotype calling method that uses a one-sided Wilcoxon signed rank test to provide a confidence score (CS) for each genotype. The CS measures the reliability of a genotype call thus filtering out SNPs with a high error rate as 'no-calls' (Matsuzaki et al., 2004). Gender status for each sample is inferred based on X chromosome heterozygosity. We performed a paired analysis for LOH and copy number changes using the Affymetrix Chromosome Copy Number Analysis Tool 4.0 (CNAT 4.0).
DNA sequencing
Genotypes generated using the Affymetrix 50K and 250K GeneChips that had discordant allele calls and a confidence score of ≤0.05 were identified as candidates for DNA sequencing to verify the SNP genotypes. After ranking by the significant confidence scores, at least 10 SNPs that did not fall within repeat elements were selected for sequencing in each twin pair. DNA samples from the parents were included in the sequencing analysis to verify Mendelian segregation of alleles. PCR reactions were performed on Applied Biosystems Gene Amp PCR System 9700 with 20ng DNA and 0.25 units of Biolase (Bioline, Randolph, MA). Sequencing was carried out with the Big Dye™ Terminator cycle sequencing and run on ABI Prism 3730 DNA Analalyzer (Applied Biosystems). Sequences were assembled using the Phred-Phrap package and visualized using the Consed program.
Array CGH results
With 2,173 BAC clones spread across the genome, the aCGH provided a resolution of approximately 1-Mb. Table 2 shows 8 BAC clones with fluorescence ratios that deviated significantly (>3 SD) in at least two twin pairs. The regions covered by these BACs do not overlap with any known CLP candidate loci, nor contain genes showing high craniofacial expression according to COGENE - the Craniofacial and Oral Gene Expression Network. Genomic region 19p13 looked especially interesting with 3 BACs showing high signal differences, so we genotyped 2 SNPs of high heterozygosity within each BAC in the twins and the parents, but did not find evidence for allelic imbalance.
Table 2
Table 2
BAC clones with significantly high signal differences from aCGH
Illumina genotyping results
The average SNP call rate on the Illumina genotyping panel was approximately 89%. We did not observe any discordant genotypes between pairs of MZ twins. A paired analysis with the Beadstudio software was performed for LOH and copy number changes at a 1 Mb window size. Although the sparse SNP coverage of ~1,500 SNPs greatly limited the power of this analysis, there was no indication of LOH in the twin pairs. A few genomic regions revealed possible copy number changes for twin pairs 3–6, but the more comprehensive genotype data generated from Affymetrix arrays disproved that observation.
Affymetrix genotyping results
Genotyping data and concordance rates for the samples scanned with the Affymetrix GeneChips are summarized in table 3. The call rate from 58,960 and 262,264 SNPs for the 50K and 250K GeneChips respectively was >90% for all samples. The average proportion of genotypes that were concordant between twin pairs was ~99% for both GeneChips. This was comparable to the 98% degree of genotype concordance observed in comparing two independent 250K array scans of the unaffected individual in twin pair no. 6, which revealed high reproducibility of the genotype calls.
Table 3
Table 3
Summary of Affymetrix GeneChips genotyping results.
Sequencing results
DNA sequencing was carried out for a total of 107 regions surrounding SNPs that had received discordant genotype calls from the Affymetrix GeneChip analyses within twin pairs. Results revealed 181 SNP genotypes that were concordant between twin pairs (Table 4). Additionally, sequencing of DNA samples from the parents showed consistency with Mendelian inheritance.
Table 4
Table 4
Summary of SNPs genotyped by sequencing to confirm concordance between twin pairs.
The advent of genome scanning tools allows comprehensive analysis of chromosomal rearrangements, dependent on the resolution of the experimental approach. The aCGH experiment using genomic BACs provided extensive coverage of the genome at a resolution of ~1 Mb, so any smaller chromosomal aberrations would not be detected. SNP arrays not only provide genotypes for thousands of SNPs, but can also be used to detect copy number changes based on hybridization signal intensities. Genotype data is useful for detecting loss of heterozygosity in chromosomal regions with deletions or uniparental disomy. At an average inter-marker distance of ~60 kb and ~12 kb for the 50K and 250K GeneChips respectively, and additional genotyping through the Illumina SNP panel, our analysis provide sufficient resolution to detect submicroscopic structural variants that are defined in the range of ~10 kb to 3 Mb (Feuk et al., 2006). Several samples were interrogated using more than one method (table 1), thus allowing data comparison and verification.
Our study using aCGH and genotyping arrays did not reveal any genomic alterations within MZ twins discordant for nonsyndromic CLP. Genomic alterations could be confined to specific tissues depending on the timing of the mutational event, resulting in mosaicism. Since acquisition of DNA samples from lip and palatal tissues is difficult, our analysis was limited to analysis of DNA samples collected from peripheral blood lymphocytes, which may not accurately represent the target tissues affected in CLP. Additionally, our experimental approaches were incapable of detecting balanced variants such as those resulting from inversions and translocations, and chromosomal aberrations that involve regions with only homozygous alleles. Routine karyotyping which can generally detect such chromosomal abnormalities could not be performed due to unavailability of living cells from a blood sample.
We are unaware of definitive reports of postzygotic genomic rearrangements underlying MZ twin discordance. A recent report of copy number variants arising between MZ twins during somatic development is consistent with our rationale for using discordant MZ twins in disease gene identification (Bruder et al., 2008). Changes in the somatic genome are well recognized as a source of diversity within the immunoglobulin and T-receptor genes (Kastern and Kryspin-Sorensen, 1988). In disease states, rearrangements can induce a phenotype by directly interrupting a gene sequence, altering gene dosage, or gene expression through position effects (Lupski and Stankiewicz, 2005). Mitotic recombination is especially relevant in tumor development, since it can lead to the expression of recessive tumor suppressor genes and/or amplification of protooncogenes (Gupta et al., 1997). Analysis of MZ twins concordant for cancer can reveal DNA rearrangements that are common to both twins as potential candidates for susceptibility loci (el-Rifai et al., 1999). In addition to oncogenesis, somatic mutations mediated through LOH can potentially contribute to other biological processes such as aging (Grist et al., 1992).
Phenotypic discordance is a common observation in complex genetic diseases, and a postzygotic somatic change is only one of several proposed mechanisms for discordance. Since CLP is a congenital defect, non-genetic intrauterine environmental factors such as unequal cell allocation at twinning and disproportionate placental blood supply may contribute to discordance (Gringras and Chen, 2001). If such developmental influences alone were responsible for MZ twin discordance, then one would expect an excess of CLP in MZ twins compared to singletons. This is not the case as studies in Denmark have demonstrated no significant difference in prevalence of CLP in twins versus singletons (Christensen and Fogh-Andersen, 1993a; Christensen and Fogh-Andersen, 1993b). The absence of genomic differences between MZ twins underscores their use as matched case-controls particularly in studies exploring the environmental component in complex traits. With the availability of array-based techniques to simultaneously scan the whole genome, it is now possible to map genomic alterations at a high resolution and additional experimental tools such as quantitative PCR and multiplex ligation-dependent probe amplification (MLPA) can be applied to independently confirm any observed genomic alterations. This straightforward approach can be applied to other cases of discordant MZ twins where identification of genomic alterations can reveal potential candidate genes, or provide evidence of a gene’s involvement in disease etiology.
ACKNOWLEDGEMENTS
We would like to thank the patients and families who contributed samples and information for this study. We are also grateful to Marla Johnson for technical assistance. This study was supported by the following grants: R37DE008559, P50DE016215, R01DE015197, and by the Intramural Research Program of the NIH, National Institute of Environmental Health Sciences (Z01-ES04007-12).
  • Boomsma D, Busjahn A, Peltonen L. Classical twin studies and beyond. Nature Review Genetics. 2002;3(11):872–882. [PubMed]
  • Bruder CE, Piotrowski A, Gijsbers AA, Andersson R, Erickson S, de Stahl TD, Menzel U, Sandgren J, von Tell D, Poplawski A, Crowley M, Crasto C, Partridge EC, Tiwari H, Allison DB, Komorowski J, van Ommen GJ, Boomsma DI, Pedersen NL, den Dunnen JT, Wirdefeldt K, Dumanski JP. Phenotypically concordant and discordant monozygotic twins display different DNA copy-number-variation profiles. American Journal of Human Genetics. 2008;82(3):763–771. [PubMed]
  • Christensen K, Fogh-Andersen P. Cleft lip (+/− cleft palate) in Danish twins, 1970–1990. American Journal of Medical Genetics. 1993a;47(6):910–916. [PubMed]
  • Christensen K, Fogh-Andersen P. Isolated cleft palate in Danish multiple births, 1970–1990. The Cleft Palate-craniofacial Journal. 1993b;30(5):469-–474. [PubMed]
  • Cote GB, Gyftodimou J. Twinning and mitotic crossing-over: some possibilities and their implications. American Journal of Human Genetics. 1991;49(1):120–130. [PubMed]
  • De Gregorio L, Jinnah HA, Harris JC, Nyhan WL, Schretlen DJ, Trombley LM, O'Neill JP. Lesch-Nyhan disease in a female with a clinically normal monozygotic twin. Molecular Genetics and Metabolism. 2005;85(1):70–77. [PubMed]
  • el-Rifai W, Tarmo L, Hemmer S, Forsti A, Pedersen N, Lichtenstein P, Ahlbom A, Soderberg M, Knuutila S, Hemminki K. DNA copy number losses at 1p32-pter in monozygotic twins concordant for breast cancer. Cancer Genetics and Cytogenetics. 1999;112(2):169–172. [PubMed]
  • Feuk L, Carson AR, Scherer SW. Structural variation in the human genome. Nature Review Genetics. 2006;7(2):85–97. [PubMed]
  • Gilbert B, Yardin C, Briault S, Belin V, Lienhardt A, Aubard Y, Battin J, Servaud M, Philippe HJ, Lacombe D. Prenatal diagnosis of female monozygotic twins discordant for Turner syndrome: implications for prenatal genetic counselling. Prenatal Diagnosis. 2002;22(8):697–702. [PubMed]
  • Gringras P, Chen W. Mechanisms for differences in monozygous twins. Early Human Development. 2001;64(2):105–117. [PubMed]
  • Grist SA, McCarron M, Kutlaca A, Turner DR, Morley AA. In vivo human somatic mutation: frequency and spectrum with age. Mutation Research. 1992;266(2):189–196. [PubMed]
  • Gupta PK, Sahota A, Boyadjiev SA, Bye S, Shao C, O'Neill JP, Hunter TC, Albertini RJ, Stambrook PJ, Tischfield JA. High frequency in vivo loss of heterozygosity is primarily a consequence of mitotic recombination. Cancer Research. 1997;57(6):1188–1193. [PubMed]
  • Jezewski PA, Vieira AR, Nishimura C, Ludwig B, Johnson M, O'Brien SE, Daack-Hirsch S, Schultz RE, Weber A, Nepomucena B, Romitti PA, Christensen K, Orioli IM, Castilla EE, Machida J, Natsume N, Murray JC. Complete sequencing shows a role for MSX1 in non-syndromic cleft lip and palate. Journal of Medical Genetics. 2003;40(6):399–407. [PMC free article] [PubMed]
  • Jugessur A, Murray JC. Orofacial clefting: recent insights into a complex trait. Current Opinion in Genetics & Development. 2005;15(3):270–278. [PMC free article] [PubMed]
  • Kastern W, Kryspin-Sorensen I. Penetrance and low concordance in monozygotic twins in disease: are they the results of alterations in somatic genomes? Molecular Reproduction and Development. 1988;1(1):63–75. [PubMed]
  • Kondo S, Schutte BC, Richardson RJ, Bjork BC, Knight AS, Watanabe Y, Howard E, de Lima RL, Daack-Hirsch S, Sander A, McDonald-McGinn DM, Zackai EH, Lammer EJ, Aylsworth AS, Ardinger HH, Lidral AC, Pober BR, Moreno L, Arcos-Burgos M, Valencia C, Houdayer C, Bahuau M, Moretti-Ferreira D, Richieri-Costa A, Dixon MJ, Murray JC. Mutations in IRF6 cause Van der Woude and popliteal pterygium syndromes. Nature Genetics. 2002;32(2):285–289. [PMC free article] [PubMed]
  • Lupski JR, Stankiewicz P. Genomic disorders: molecular mechanisms for rearrangements and conveyed phenotypes. PLoS Genetics. 2005;1(6):e49. [PMC free article] [PubMed]
  • Machin GA. Some causes of genotypic and phenotypic discordance in monozygotic twin pairs. American Journal of Medical Genetics. 1996;61(3):216–228. [PubMed]
  • Mak YT, Hampson G, Beresford JN, Spector TD. Variations in genome-wide gene expression in identical twins - a study of primary osteoblast-like culture from female twins discordant for osteoporosis. BMC Genetics. 2004;5(14) [PMC free article] [PubMed]
  • Mansilla MA, Kimani J, Mitchell LE, Christensen K, Boomsma DI, Daack-Hirsch S, Nepomucena B, Wyszynski DF, Felix TM, Martin NG, Murray JC. Discordant MZ twins with cleft lip and palate: a model for identifying genes in complex traits. Twin Research and Human Genetics. 2005;8(1):39–46. [PubMed]
  • Matsuzaki H, Dong S, Loi H, Di X, Liu G, Hubbell E, Law J, Berntsen T, Chadha M, Hui H, Yang G, Kennedy GC, Webster TA, Cawley S, Walsh PS, Jones KW, Fodor SP, Mei R. Genotyping over 100,000 SNPs on a pair of oligonucleotide arrays. Nature Methods. 2004;1(2):109–111. [PubMed]
  • Osoegawa K, Vessere GM, Utami KH, Mansilla MA, Johnson MK, Riley BM, L'Heureux J, Pfundt R, Staaf J, van der Vliet WA, Lidral AC, Schoenmakers EF, Borg A, Schutte BC, Lammer EJ, Murray JC, de Jong PJ. Identification of novel candidate genes associated with cleft lip and palate using array comparative genomic hybridisation. Journal of Medical Genetics. 2008;45(2):81–86. [PMC free article] [PubMed]
  • Pinkel D, Albertson DG. Comparative genomic hybridization. Annual Review of Genomics and Human Genetics. 2005;6:331–354. [PubMed]
  • Riley BM, Mansilla MA, Ma J, Daack-Hirsch S, Maher BS, Raffensperger LM, Russo ET, Vieira AR, Dode C, Mohammadi M, Marazita ML, Murray JC. Impaired FGF signaling contributes to cleft lip and palate. Proceedings of the National Academy of Sciences of the United States of America. 2007;104(11):4512–4517. [PubMed]
  • Steemers FJ, Gunderson KL. Whole genome genotyping technologies on the BeadArray platform. Biotechnology Journal. 2007;2(1):41–49. [PubMed]
  • Watanabe A, Akita S, Tin NT, Natsume N, Nakano Y, Niikawa N, Uchiyama T, Yoshiura K. A mutation in RYK is a genetic factor for nonsyndromic cleft lip and palate. The Cleft Palate-craniofacial Journal. 2006;43(3):310–316. [PubMed]
  • Weksberg R, Shuman C, Caluseriu O, Smith AC, Fei YL, Nishikawa J, Stockley TL, Best L, Chitayat D, Olney A, Ives E, Schneider A, Bestor TH, Li M, Sadowski P, Squire J. Discordant KCNQ1OT1 imprinting in sets of monozygotic twins discordant for Beckwith-Wiedemann syndrome. Human Molecular Genetics. 2002;11(11):1317–1325. [PubMed]
  • Wong AH, Gottesman II, Petronis A. Phenotypic differences in genetically identical organisms: the epigenetic perspective. Human Molecular Genetics. 2005;14(Spec No 1):R11–R18. [PubMed]
  • Zucchero TM, Cooper ME, Maher BS, Daack-Hirsch S, Nepomuceno B, Ribeiro L, Caprau D, Christensen K, Suzuki Y, Machida J, Natsume N, Yoshiura K, Vieira AR, Orioli IM, Castilla EE, Moreno L, Arcos-Burgos M, Lidral AC, Field LL, Liu YE, Ray A, Goldstein TH, Schultz RE, Shi M, Johnson MK, Kondo S, Schutte BC, Marazita ML, Murray JC. Interferon regulatory factor 6 (IRF6) gene variants and the risk of isolated cleft lip or palate. The New England Journal of Medicine. 2004;351(8):769–780. [PubMed]