Methods for detecting gene x gene interaction in multiplex extended pedigrees.

Complex diseases are multifactorial in nature and can involve multiple loci with gene x gene and gene x environment interactions. Research on methods to uncover the interactions between those genes that confer susceptibility to disease has been extensive, but many of these methods have only been developed for sibling pairs or sibships. In this report, we assess the performance of two methods for finding gene x gene interactions that are applicable to arbitrarily sized pedigrees, one based on correlation in per-family nonparametric linkage scores and another that incorporates candidate loci genotypes as covariates into an affected relative pair linkage analysis. The power and type I error rate of both of these methods was addressed using the simulated Genetic Analysis Workshop 14 data. In general, we found detection of the interacting loci to be a difficult problem, and though we experienced some modest success there is a clear need to continue developing new methods and approaches to the problem.

Identifying genomic regions for fine-mapping using genome scan meta-analysis (GSMA) to identify the minimum regions of maximum significance (MRMS) across populations.

In order to detect linkage of the simulated complex disease Kofendrerd Personality Disorder across studies from multiple populations, we performed a genome scan meta-analysis (GSMA). Using the 7-cM microsatellite map, nonparametric multipoint linkage analyses were performed separately on each of the four simulated populations independently to determine p-values. The genome of each population was divided into 20-cM bin regions, and each bin was rank-ordered based on the most significant linkage p-value for that population in that region. The bin ranks were then averaged across all four studies to determine the most significant 20-cM regions over all studies. Statistical significance of the averaged bin ranks was determined from a normal distribution of randomly assigned rank averages. To narrow the region of interest for fine-mapping, the meta-analysis was repeated two additional times, with each of the 20-cM bins offset by 7 cM and 13 cM, respectively, creating regions of overlap with the original method. The 6-7 cM shared regions, where the highest averaged 20-cM bins from each of the three offsets overlap, designated the minimum region of maximum significance (MRMS). Application of the GSMA-MRMS method revealed genome wide significance (p-values refer to the average rank assigned to the bin) at regions including or adjacent to all of the simulated disease loci: chromosome 1 (p < 0.0001 for 160-167 cM, including D1), chromosome 3 (p-value < 0.0000001 for 287-294 cM, including D2), chromosome 5 (p-value < 0.001 for 0-7 cM, including D3), and chromosome 9 (p-value < 0.05 for 7-14 cM, the region adjacent to D4). This GSMA analysis approach demonstrates the power of linkage meta-analysis to detect multiple genes simultaneously for a complex disorder. The MRMS method enhances this powerful tool to focus on more localized regions of linkage.

Genome scan for loci involved in nonsyndromic cleft lip with or without cleft palate in families from West Bengal, India.

In order to identify genes or regions involved in nonsyndromic cleft lip with or without cleft palate (CL/P) in families from India, we analyzed 38 multiplex families (DNA from 272 individuals, 82 affected with CL/P, 190 unaffected) for 285 genome-wide markers (average spacing 12.6 cM), including markers in six candidate loci or regions on chromosomes 2, 4, 6, 14, 17, and 19 that have been implicated in other studies of CL/P. LOD scores (two-point and multipoint), and model-free association (TDT) and linkage (NPL) statistics, were calculated between each of the markers and a hypothetical CL/P susceptibility locus. The most statistically significant two-point linkage results were with markers on chromosome 7 (LOD = 1.89 with D7S435, 7p15, 47 cM), chromosome 5 (LOD = 1.76 with D5S407, 5q11, 65 cM), chromosome 15 (LOD = 1.55 with D15S652, 15q26, 90 cM), and chromosome 20 (LOD = 1.46 with STS155130, 20q13, 54 cM). The most significant multipoint linkage result was on chromosome 5q, again near D5S407 (HLOD = 1.40). Regions on chromosomes 1p, 1q, 7q, 12q, 16q, 18q, and Xp also had a LOD or HLOD > or = 1.0. Of seven candidate markers and regions with previous positive reports in the literature (TGFA, MSX1, D4S175, F13A1, TGFB3, D17S250, and APOC2), none had a significant linkage result, but one (the APOC2 region) had a significant association result and three others (TGFA, MSX1, F13A1) had suggestive results. The results are consistent with the involvement of multiple loci in CL/P expression in this West Bengal population, which concurs with results found in other CL/P study populations.

Interferon regulatory factor 6 (IRF6) gene variants and the risk of isolated cleft lip or palate.

BACKGROUND: Cleft lip or palate (or the two in combination) is a common birth defect that results from a mixture of genetic and environmental factors. We searched for a specific genetic factor contributing to this complex trait by examining large numbers of affected patients and families and evaluating a specific candidate gene. METHODS: We identified the gene that encodes interferon regulatory factor 6 (IRF6) as a candidate gene on the basis of its involvement in an autosomal dominant form of cleft lip and palate, Van der Woude's syndrome. A single-nucleotide polymorphism in this gene results in either a valine or an isoleucine at amino acid position 274 (V274I). We carried out transmission-disequilibrium testing for V274I in 8003 individual subjects in 1968 families derived from 10 populations with ancestry in Asia, Europe, and South America, haplotype and linkage analyses, and case-control analyses, and determined the risk of cleft lip or palate that is associated with genetic variation in IRF6. RESULTS: Strong evidence of overtransmission of the valine (V) allele was found in the entire population data set (P<10(-9)); moreover, the results for some individual populations from South America and Asia were highly significant. Variation at IRF6 was responsible for 12 percent of the genetic contribution to cleft lip or palate and tripled the risk of recurrence in families that had already had one affected child. CONCLUSIONS: DNA-sequence variants associated with IRF6 are major contributors to cleft lip, with or without cleft palate. The contribution of variants in single genes to cleft lip or palate is an important consideration in genetic counseling.

Meta-analysis of 13 genome scans reveals multiple cleft lip/palate genes with novel loci on 9q21 and 2q32-35.

Isolated or nonsyndromic cleft lip with or without cleft palate (CL/P) is a common birth defect with a complex etiology. A 10-cM genome scan of 388 extended multiplex families with CL/P from seven diverse populations (2,551 genotyped individuals) revealed CL/P genes in six chromosomal regions, including a novel region at 9q21 (heterogeneity LOD score [HLOD]=6.6). In addition, meta-analyses with the addition of results from 186 more families (six populations; 1,033 genotyped individuals) showed genomewide significance for 10 more regions, including another novel region at 2q32-35 (P=.0004). These are the first genomewide significant linkage results ever reported for CL/P, and they represent an unprecedented demonstration of the power of linkage analysis to detect multiple genes simultaneously for a complex disorder.

Genome-scan for loci involved in cleft lip with or without cleft palate in consanguineous families from Turkey.

Cleft lip with or without cleft palate (CL/P) is a common congenital anomaly, with birth prevalence ranging from 1/500 to 1/1,000. A number of genetic loci have shown positive linkage or association results in European Caucasian populations. The purpose of the current study was to assess whether any of those loci have positive results in Turkish Caucasian CL/P families, and to perform a 10 cM genome scan to identify other regions potentially containing cleft susceptibility loci. Eighteen affected individuals with consanguineous parents were identified as part of our on-going studies of orofacial clefts in Ankara, Turkey. Genotyped were 383 genome-scan markers, and 70 additional markers, including markers in six candidate loci or regions on chromosomes 2, 4, 6, 14, 17, and 19 (TGFA, D4S175, F13A1, TGFB3, D17S250, and APOC2) that have been implicated in other studies of families with orofacial clefting. LOD scores (two point and multiple point) and family-based association statistics (TDT) were calculated between each of the markers and CL/P. For the LOD score calculations, an autosomal recessive model was assumed for the inheritance of CL/P. Of the six candidate markers, significant TDT results were obtained with TGFA (P = 0.05). The most statistically significant multipoint results from the linkage genome scan were between putative genes controlling risk of CL/P and regions on chromosomes 4, 10, 12, and 15 (maximum multipoint HLOD's of 1.25, 1.30, 2.73, and 1.28 respectively). These results demonstrate the power of small numbers of families with inbred probands to detect linkage and association.