Whole exome sequencing of distant relatives in multiplex families implicates rare variants in candidate genes for oral clefts

A dozen genes/regions have been confirmed as genetic risk factors for oral clefts in human association and linkage studies, and animal models argue even more genes may be involved. Genomic sequencing studies should identify specific causal variants and may reveal additional genes as influencing risk to oral clefts, which have a complex and heterogeneous etiology. We conducted a whole exome sequencing (WES) study to search for potentially causal variants using affected relatives drawn from multiplex cleft families. Two or three affected second, third, and higher degree relatives from 55 multiplex families were sequenced. We examined rare single nucleotide variants (SNVs) shared by affected relatives in 348 recognized candidate genes. Exact probabilities that affected relatives would share these rare variants were calculated, given pedigree structures, and corrected for the number of variants tested. Five novel and potentially damaging SNVs shared by affected distant relatives were found and confirmed by Sanger sequencing. One damaging SNV in CDH1, shared by three affected second cousins from a single family, attained statistical significance (P = 0.02 after correcting for multiple tests). Family-based designs such as the one used in this WES study offer important advantages for identifying genes likely to be causing complex and heterogeneous disorders.

Mutations in Alström protein impair terminal differentiation of cardiomyocytes

Cardiomyocyte cell division and replication in mammals proceed through embryonic development and abruptly decline soon after birth. The process governing cardiomyocyte cell cycle arrest is poorly understood. Here we carry out whole-exome sequencing in an infant with evidence of persistent postnatal cardiomyocyte replication to determine the genetic risk factors. We identify compound heterozygous ALMS1 mutations in the proband, and confirm their presence in her affected sibling, one copy inherited from each heterozygous parent. Next, we recognize homozygous or compound heterozygous truncating mutations in ALMS1 in four other children with high levels of postnatal cardiomyocyte proliferation. Alms1 mRNA knockdown increases multiple markers of proliferation in cardiomyocytes, the percentage of cardiomyocytes in G2/M phases, and the number of cardiomyocytes by 10% in cultured cells. Homozygous Alms1-mutant mice have increased cardiomyocyte proliferation at 2 weeks postnatal compared with wild-type littermates. We conclude that deficiency of Alström protein impairs postnatal cardiomyocyte cell cycle arrest.

A genome-wide study of de novo deletions identifies a candidate locus for non-syndromic isolated cleft lip/palate risk

Background

Copy number variants (CNVs) may play an important part in the development of common birth defects such as oral clefts, and individual patients with multiple birth defects (including clefts) have been shown to carry small and large chromosomal deletions. In this paper we investigate de novo deletions defined as DNA segments missing in an oral cleft proband but present in both unaffected parents. We compare de novo deletion frequencies in children of European ancestry with an isolated, non-syndromic oral cleft to frequencies in children of European ancestry from randomly sampled trios.

Results

We identified a genome-wide significant 62 kilo base (kb) non-coding region on chromosome 7p14.1 where de novo deletions occur more frequently among oral cleft cases than controls. We also observed wider de novo deletions among cleft lip and palate (CLP) cases than seen among cleft palate (CP) and cleft lip (CL) cases.

Conclusions

This study presents a region where de novo deletions appear to be involved in the etiology of oral clefts, although the underlying biological mechanisms are still unknown. Larger de novo deletions are more likely to interfere with normal craniofacial development and may result in more severe clefts. Study protocol and sample DNA source can severely affect estimates of de novo deletion frequencies. Follow-up studies are needed to further validate these findings and to potentially identify additional structural variants underlying oral clefts.

Evidence of gene-environment interaction for two genes on chromosome 4 and environmental tobacco smoke in controlling the risk of nonsyndromic cleft palate

Nonsyndromic cleft palate (CP) is one of the most common human birth defects and both genetic and environmental risk factors contribute to its etiology. We conducted a genome-wide association study (GWAS) using 550 CP case-parent trios ascertained in an international consortium. Stratified analysis among trios with different ancestries was performed to test for GxE interactions with common maternal exposures using conditional logistic regression models. While no single nucleotide polymorphism (SNP) achieved genome-wide significance when considered alone, markers in SLC2A9 and the neighboring WDR1 on chromosome 4p16.1 gave suggestive evidence of gene-environment interaction with environmental tobacco smoke (ETS) among 259 Asian trios when the models included a term for GxE interaction. Multiple SNPs in these two genes were associated with increased risk of nonsyndromic CP if the mother was exposed to ETS during the peri-conceptual period (3 months prior to conception through the first trimester). When maternal ETS was considered, fifteen of 135 SNPs mapping to SLC2A9 and 9 of 59 SNPs in WDR1 gave P values approaching genome-wide significance (10(-6)<P<10(-4)) in a test for GxETS interaction. SNPs rs3733585 and rs12508991 in SLC2A9 yielded P = 2.26×10(-7) in a test for GxETS interaction. SNPs rs6820756 and rs7699512 in WDR1 also yielded P = 1.79×10(-7) and P = 1.98×10(-7) in a 1 df test for GxE interaction. Although further replication studies are critical to confirming these findings, these results illustrate how genetic associations for nonsyndromic CP can be missed if potential GxE interaction is not taken into account, and this study suggest SLC2A9 and WDR1 should be considered as candidate genes for CP.

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.

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.

Genome-wide study identifies two loci associated with lung function decline in mild to moderate COPD

Accelerated lung function decline is a key COPD phenotype; however, its genetic control remains largely unknown. We performed a genome-wide association study using the Illumina Human660W-Quad v.1_A BeadChip. Generalized estimation equations were used to assess genetic contributions to lung function decline over a 5-year period in 4,048 European American Lung Health Study participants with largely mild COPD. Genotype imputation was performed using reference HapMap II data. To validate regions meeting genome-wide significance, replication of top SNPs was attempted in independent cohorts. Three genes (TMEM26, ANK3 and FOXA1) within the regions of interest were selected for tissue expression studies using immunohistochemistry. Two intergenic SNPs (rs10761570, rs7911302) on chromosome 10 and one SNP on chromosome 14 (rs177852) met genome-wide significance after Bonferroni. Further support for the chromosome 10 region was obtained by imputation, the most significantly associated imputed SNPs (rs10761571, rs7896712) being flanked by observed markers rs10761570 and rs7911302. Results were not replicated in four general population cohorts or a smaller cohort of subjects with moderate to severe COPD; however, we show novel expression of genes near regions of significantly associated SNPS, including TMEM26 and FOXA1 in airway epithelium and lung parenchyma, and ANK3 in alveolar macrophages. Levels of expression were associated with lung function and COPD status. We identified two novel regions associated with lung function decline in mild COPD. Genes within these regions were expressed in relevant lung cells and their expression related to airflow limitation suggesting they may represent novel candidate genes for COPD susceptibility.

Fast detection of de novo copy number variants from SNP arrays for case-parent trios

Background

In studies of case-parent trios, we define copy number variants (CNVs) in the offspring that differ from the parental copy numbers as de novo and of interest for their potential functional role in disease. Among the leading array-based methods for discovery of de novo CNVs in case-parent trios is the joint hidden Markov model (HMM) implemented in the PennCNV software. However, the computational demands of the joint HMM are substantial and the extent to which false positive identifications occur in case-parent trios has not been well described. We evaluate these issues in a study of oral cleft case-parent trios.

Results

Our analysis of the oral cleft trios reveals that genomic waves represent a substantial source of false positive identifications in the joint HMM, despite a wave-correction implementation in PennCNV. In addition, the noise of low-level summaries of relative copy number (log R ratios) is strongly associated with batch and correlated with the frequency of de novo CNV calls. Exploiting the trio design, we propose a univariate statistic for relative copy number referred to as the minimum distance that can reduce technical variation from probe effects and genomic waves. We use circular binary segmentation to segment the minimum distance and maximum a posteriori estimation to infer de novo CNVs from the segmented genome. Compared to PennCNV on simulated data, MinimumDistance identifies fewer false positives on average and is comparable to PennCNV with respect to false negatives. Genomic waves contribute to discordance of PennCNV and MinimumDistance for high coverage de novo calls, while highly concordant calls on chromosome 22 were validated by quantitative PCR. Computationally, MinimumDistance provides a nearly 8-fold increase in speed relative to the joint HMM in a study of oral cleft trios.

Conclusions

Our results indicate that batch effects and genomic waves are important considerations for case-parent studies of de novo CNV, and that the minimum distance is an effective statistic for reducing technical variation contributing to false de novo discoveries. Coupled with segmentation and maximum a posteriori estimation, our algorithm compares favorably to the joint HMM with MinimumDistance being much faster.

Mutations in the TGF-β repressor SKI cause Shprintzen-Goldberg syndrome with aortic aneurysm

Increased transforming growth factor beta (TGF-β) signaling has been implicated in the pathogenesis of syndromic presentations of aortic aneurysm, including Marfan syndrome (MFS) and Loeys-Dietz syndrome (LDS)^1-4. However, the location and character of many of the causal mutations in LDS would intuitively infer diminished TGF-β signaling⁵. Taken together, these data have engendered controversy regarding the specific role of TGF-β in disease pathogenesis. Shprintzen-Goldberg syndrome (SGS) has considerable phenotypic overlap with MFS and LDS, including aortic aneurysm^6-8. We identified causative variation in 10 patients with SGS in the proto-oncogene SKI, a known repressor of TGF-β activity^9,10. Cultured patient dermal fibroblasts showed enhanced activation of TGF-β signaling cascades and increased expression of TGF-β responsive genes. Morpholino-induced silencing of SKI paralogs in zebrafish recapitulated abnormalities seen in SGS patients. These data support the conclusion that increased TGF-β signaling is the mechanism underlying SGS and contributes to multiple syndromic presentations of aortic aneurysm.

Genome-wide meta-analyses of nonsyndromic cleft lip with or without cleft palate identify six new risk loci

We have conducted the first meta-analyses for nonsyndromic cleft lip with or without cleft palate (NSCL/P) using data from the two largest genome-wide association studies published to date. We confirmed associations with all previously identified loci and identified six additional susceptibility regions (1p36, 2p21, 3p11.1, 8q21.3, 13q31.1 and 15q22). Analysis of phenotypic variability identified the first specific genetic risk factor for NSCLP (nonsyndromic cleft lip plus palate) (rs8001641; P(NSCLP) = 6.51 × 10(-11); homozygote relative risk = 2.41, 95% confidence interval (CI) 1.84-3.16).

Detectable clonal mosaicism from birth to old age and its relationship to cancer

We detected clonal mosaicism for large chromosomal anomalies (duplications, deletions and uniparental disomy) using SNP microarray data from over 50,000 subjects recruited for genome-wide association studies. This detection method requires a relatively high frequency of cells with the same abnormal karyotype (>5-10%; presumably of clonal origin) in the presence of normal cells. The frequency of detectable clonal mosaicism in peripheral blood is low (<0.5%) from birth until 50 years of age, after which it rapidly rises to 2-3% in the elderly. Many of the mosaic anomalies are characteristic of those found in hematological cancers and identify common deleted regions with genes previously associated with these cancers. Although only 3% of subjects with detectable clonal mosaicism had any record of hematological cancer before DNA sampling, those without a previous diagnosis have an estimated tenfold higher risk of a subsequent hematological cancer (95% confidence interval = 6-18).

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.

BMP4 was associated with NSCL/P in an Asian population

BACKGROUND

The Bone Morphogenetic Protein 4 gene (BMP4) is located in chromosome 14q22-q23 which has shown evidence of linkage for isolated nonsyndromic cleft lip with or without cleft palate (NSCL/P) in a genome wide linkage analysis of human multiplex families. BMP4 has been shown to play crucial roles in lip and palatal development in animal models. Several candidate gene association analyses also supported its potential risk for NSCL/P, however, results across these association studies have been inconsistent. The aim of the current study was to test for possible association between markers in and around the BMP4 gene and NSCL/P in Asian and Maryland trios.

METHODOLOGY/PRINCIPAL FINDINGS

Family Based Association Test was used to test for deviation from Mendelian assortment for 12 SNPs in and around BMP4. Nominal significant evidence of linkage and association was seen for three SNPs (rs10130587, rs2738265 and rs2761887) in 221 Asian trios and for one SNP (rs762642) in 76 Maryland trios. Statistical significance still held for rs10130587 after Bonferroni correction (corrected p = 0.019) among the Asian group. Estimated odds ratio for carrying the apparent high risk allele at this SNP was 1.61 (95%CI = 1.20, 2.18).

CONCLUSIONS

Our results provided further evidence of association between BMP4 and NSCL/P.

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.

ROR2 gene is associated with risk of non-syndromic cleft palate in an Asian population

BACKGROUND

The receptor tyrosine kinase-like orphan receptor 2 (ROR2) gene has been recently shown to play important roles in palatal development in animal models and resides in the chromosomal region linked to non syndromic cleft lip with or without cleft palate in humans. The aim of this study was to investigate the possible association between ROR2 gene and non-syndromic oral clefts.

METHODS

Here we tested 38 eligible single-nucleotide polymorphisms (SNPs) in ROR2 gene in 297 non-syndromic cleft lip with or without cleft palate and in 82 non-syndromic cleft palate case parent trios recruited from Asia and Maryland. Family Based Association Test was used to test for deviation from Mendelian inheritance. Plink software was used to test potential parent of origin effect. Possible maternally mediated in utero effects were assessed using the TRIad Multi-Marker approach under an assumption of mating symmetry in the population.

RESULTS

Significant evidence of linkage and association was shown for 3 SNPs (rs7858435, rs10820914 and rs3905385) among 57 Asian non-syndromic cleft palate trios in Family Based Association Tests. P values for these 3 SNPs equaled to 0.000068, 0.000115 and 0.000464 respectively which were all less than the significance level (0.05/38 = 0.0013) adjusted by strict Bonferroni correction. Relevant odds ratios for the risk allele were 3.42 (1.80 - 6.50), 3.45 (1.75 - 6.67) and 2.94 (1.56 - 5.56), respectively. Statistical evidence of linkage and association was not shown for study groups other than non-syndromic cleft palate. Neither evidence for parent-of-origin nor maternal genotypic effect was shown for any of the ROR2 markers in our analysis for all study groups.

CONCLUSION

Our results provided evidence of linkage and association between the ROR2 gene and a gene controlling risk to non-syndromic cleft palate.

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.

The FGF and FGFR Gene Family and Risk of Cleft Lip With or Without Cleft Palate

Background : Isolated, nonsyndromic cleft lip with or without cleft palate is a common human congenital malformation with a complex and heterogeneous etiology. Genes coding for fibroblast growth factors and their receptors (FGF/FGFR genes) are excellent candidate genes. Methods : We tested single-nucleotide polymorphic markers in 10 FGF/FGFR genes (including FGFBP1, FGF2, FGF10, FGF18, FGFR1, FGFR2, FGF19, FGF4, FGF3, and FGF9) for genotypic effects, interactions with one another, and with common maternal environmental exposures in 221 Asian and 76 Maryland case-parent trios ascertained through a child with isolated, nonsyndromic cleft lip with or without cleft palate. Results : Both FGFR1 and FGF19 yielded evidence of linkage and association in the transmission disequilibrium test, confirming previous evidence. Haplotypes of three single-nucleotide polymorphisms in FGFR1 were nominally significant among Asian trios. Estimated odds ratios for individual single-nucleotide polymorphic markers and haplotypes of multiple markers in FGF19 ranged from 1.31 to 1.87. We also found suggestive evidence of maternal genotypic effects for markers in FGF2 and FGF10 among Asian trios. Tests for gene-environment (G × E) interaction between markers in FGFR2 and maternal smoking or multivitamin supplementation yielded significant evidence of G × E interaction separately. Tests of gene-gene (G × G) interaction using Cordell's method yielded significant evidence between single-nucleotide polymorphisms in FGF9 and FGF18, which was confirmed in an independent sample of trios from an international consortium. Conclusion : Our results suggest several genes in the FGF/FGFR family may influence risk for isolated, nonsyndromic cleft lip with or without cleft palate through distinct biological mechanisms.

Performance assessment of copy number microarray platforms using a spike-in experiment

MOTIVATION

Changes in the copy number of chromosomal DNA segments [copy number variants (CNVs)] have been implicated in human variation, heritable diseases and cancers. Microarray-based platforms are the current established technology of choice for studies reporting these discoveries and constitute the benchmark against which emergent sequence-based approaches will be evaluated. Research that depends on CNV analysis is rapidly increasing, and systematic platform assessments that distinguish strengths and weaknesses are needed to guide informed choice.

RESULTS

We evaluated the sensitivity and specificity of six platforms, provided by four leading vendors, using a spike-in experiment. NimbleGen and Agilent platforms outperformed Illumina and Affymetrix in accuracy and precision of copy number dosage estimates. However, Illumina and Affymetrix algorithms that leverage single nucleotide polymorphism (SNP) information make up for this disadvantage and perform well at variant detection. Overall, the NimbleGen 2.1M platform outperformed others, but only with the use of an alternative data analysis pipeline to the one offered by the manufacturer.

AVAILABILITY

The data is available from http://rafalab.jhsph.edu/cnvcomp/.

CONTACT

pevsner@jhmi.edu; fspencer@jhmi.edu; rafa@jhu.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

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.

Association between genes on chromosome 4p16 and non-syndromic oral clefts in four populations

Isolated cleft lip with or without cleft palate and cleft palate are among the most common human birth defects. Several candidate gene studies on MSX1 have shown significant association between markers in MSX1 and risk of oral clefts, and re-sequencing studies have identified multiple mutations in MSX1 in a small minority of cases, which may account for 1-2% of all isolated oral clefts cases. We explored the 2-Mb region around MSX1, using a marker map of 393 single nucleotide polymorphisms (SNPs) in 297 cleft lip, with or without cleft palate, case-parent trios and 84 cleft palate trios from Maryland, Taiwan, Singapore, and Korea. Both individual markers and haplotypes of two to five SNPs showed several regions yielding statistical evidence for linkage and disequilibrium. Two genes (STK32B and EVC) yielded consistent evidence from cleft lip, with or without cleft palate, trios in all four populations. These two genes plus EVC2 also yielded suggestive evidence for linkage and disequilibrium among cleft palate trios. This analysis suggests that several genes, not just MSX1, in this region may influence risk of oral clefts.

Variants in the gene encoding C3 are associated with asthma and related phenotypes among African Caribbean families

Proinflammatory and immunoregulatory products from C3 play a major role in phagocytosis, respiratory burst, and airways inflammation. C3 is critical in adaptive immunity; studies in mice deficient in C3 demonstrate that features of asthma are significantly attenuated in the absence of C3. To test the hypothesis that the C3 gene on chromosome 19p13.3-p13.2 contains variants associated with asthma and related phenotypes, we genotyped 25 single nucleotide polymorphism (SNP) markers distributed at intervals of approximately 1.9 kb within the C3 gene in 852 African Caribbean subjects from 125 nuclear and extended pedigrees. We used the multiallelic test in the family-based association test program to examine sliding windows comprised of 2-6 SNPs. A five-SNP window between markers rs10402876 and rs366510 provided strongest evidence for linkage in the presence of linkage disequilibrium for asthma, high log[total IgE], and high log[IL-13]/[log[IFN-gamma] in terms of global P-values (P = 0.00027, 0.00013, and 0.003, respectively). A three-SNP haplotype GGC for the first three of these markers showed best overall significance for the three phenotypes (P = 0.003, 0.007, 0.018, respectively) considering haplotype-specific tests. Taken together, these results implicate the C3 gene as a priority candidate controlling risk for asthma and allergic disease in this population of African descent.

Maternal transmission effects of the PAX genes among cleft case-parent trios from four populations

Isolated cleft lip with or without cleft palate (CL/P) is among the most common human birth defects, with a prevalence of 1 in 700 live births. The paired box (PAX) genes have been suggested as candidate genes for CL/P based largely on mouse models; however, few human studies have focused on this gene family. This study tests for association between markers in four PAX genes and CL/P using a case-parent trio design considering parent-of-origin effects. Trios from four populations (76 from Maryland, 146 from Taiwan, 35 from Singapore, and 40 from Korea) were genotyped for 34 single nucleotide polymorphisms (SNPs) in the PAX3, PAX6, PAX7, and PAX9 genes. We performed the transmission disequilibrium test (TDT) on individual SNPs. Parent-of-origin effects were assessed using the transmission asymmetry test (TAT) and the parent-of-origin likelihood ratio test (PO-LRT). TDT analysis showed one SNP (rs766325) in PAX7 yielding evidence of linkage and association when parent-of-origin was not considered, with an OR(transmission)=1.62 (P=0.003), and five SNPs in PAX6 (including two pairs in near perfect linkage disequilibrium). TAT analysis of all trios revealed two SNPs in PAX7 and four SNPs in PAX3 showing significant excess maternal transmission. For these six SNPs, the maternal OR(transmission) ranged between 1.74 and 2.40, and PO-LRT was also significant (P-values=0.035-0.012). When this analysis was limited to trios with male cases, SNPs in PAX7 showed higher maternal OR(transmission) and greater significance. PAX genes may influence the risk of CL/P through maternal effects, possibly imprinting, which seems to be stronger among male cases.

Differential parental transmission of markers in RUNX2 among cleft case-parent trios from four populations

Isolated cleft lip with or without cleft palate (CL/P) is among the most common human birth defects, with a prevalence around 1 in 700 live births. The Runt-related transcription factor 2 (RUNX2) gene has been suggested as a candidate gene for CL/P based largely on mouse models; however, no human studies have focused on RUNX2 as a risk factor for CL/P. This study examines the association between markers in RUNX2 and isolated, nonsyndromic CL/P using a case-parent trio design, while considering parent-of-origin effects. Case-parent trios from four populations (77 from Maryland, 146 from Taiwan, 35 from Singapore, and 40 from Korea) were genotyped for 24 single nucleotide polymorphisms (SNPs) in the RUNX2 gene. We performed the transmission disequilibrium test on individual SNPs. Parent-of-origin effects were assessed using the transmission asymmetry test and the parent-of-origin likelihood ratio test (PO-LRT). When all trios were combined, the transmission asymmetry test revealed a block of 11 SNPs showing excess maternal transmission significant at the P<0.01 level, plus one SNP (rs1934328) showing excess paternal transmission (P=0.002). For the 11 SNPs showing excess maternal transmission, odds ratios of being transmitted to the case from the mother ranged between 3.00 and 4.00. The parent-of-origin likelihood ratio tests for equality of maternal and paternal transmission were significant for three individual SNPs (rs910586, rs2819861, and rs1934328). Thus, RUNX2 appears to influence risk of CL/P through a parent-of-origin effect with excess maternal transmission.

McKusick's Online Mendelian Inheritance in Man (OMIM)

McKusick's Online Mendelian Inheritance in Man (OMIM®; http://www.ncbi.nlm.nih.gov/omim), a knowledgebase of human genes and phenotypes, was originally published as a book, Mendelian Inheritance in Man, in 1966. The content of OMIM is derived exclusively from the published biomedical literature and is updated daily. It currently contains 18 961 full-text entries describing phenotypes and genes. To date, 2239 genes have mutations causing disease, and 3770 diseases have a molecular basis. Approximately 70 new entries are added and 700 entries are updated per month. OMIM® is expanding content and organization in response to shifting biological paradigms and advancing biotechnology.

Excess maternal transmission of markers in TCOF1 among cleft palate case-parent trios from three populations

Isolated cleft palate is among the most common human birth defects. The TCOF1 gene has been suggested as a candidate gene for cleft palate based on animal models. This study tests for association between markers in TCOF1 and isolated, nonsyndromic cleft palate using a case-parent trio design considering parent-of-origin effects. Case-parent trios from three populations (comprising a total of 81 case-parent trios) were genotyped for single nucleotide polymorphisms (SNPs) in the TCOF1 gene. We used the transmission disequilibrium test and the transmission asymmetry test on individual SNPs. When all trios were combined, the odds ratio for transmission of the minor allele, OR(transmission), was significant for SNP rs15251 (OR = 2.88, P = 0.007), as well as rs2255796 and rs2569062 (OR = 2.08, P = 0.03; OR = 2.43, P = 0.041; respectively) when parent of origin was not considered. The transmission asymmetry test also revealed one SNP (rs15251) showing excess maternal transmission significant at the P = 0.005 level (OR = 6.50). Parent-of-origin effects were assessed using the parent-of-origin likelihood ratio test on both SNPs and haplotypes. While the parent-of-origin likelihood ratio test was only marginally significant for this SNP (P = 0.136), analysis of haplotypes of rs2255796 and rs15251 suggested excess maternal transmission. Therefore, these data suggest TCOF1 may influence risk of cleft palate through a parent-of-origin effect.