Novel caries loci in children and adults implicated by genome-wide analysis of families.

BACKGROUND: Dental caries is a common chronic disease among children and adults alike, posing a substantial health burden. Caries is affected by multiple genetic and environmental factors, and prior studies have found that a substantial proportion of caries susceptibility is genetically inherited. METHODS: To identify such genetic factors, we conducted a genome-wide linkage scan in 464 extended families with 2616 individuals from Iowa, Pennsylvania and West Virginia for three dental caries phenotypes: (1) PRIM: dichotomized as zero versus one or more affected primary teeth, (2) QTOT1: age-adjusted quantitative caries measure for both primary and permanent dentitions including pre-cavitated lesions, and (3) QTOT2: age-adjusted quantitative caries excluding pre-cavitated lesions. Genotyping was conducted for approximately 600,000 SNPs on an Illumina platform, pruned to 127,511 uncorrelated SNPs for the analyses reported here. RESULTS: Multipoint non-parametric linkage analyses generated peak LOD scores exceeding 2.0 for eight genomic regions, but no LOD scores above 3.0 were observed. The maximum LOD score for each of the three traits was 2.90 at 1q25.3 for PRIM, 2.38 at 6q25.3 for QTOT1, and 2.76 at 5q23.3 for QTOT2. Some overlap in linkage regions was observed among the phenotypes. Genes with a potential role in dental caries in the eight chromosomal regions include CACNA1E, LAMC2, ALMS1, STAMBP, GXYLT2, SLC12A2, MEGF10, TMEM181, ARID1B, and, as well as genes in several immune gene families. Our results are also concordant with previous findings from association analyses on chromosomes 11 and 19. CONCLUSIONS: These multipoint linkage results provide evidence in favor of novel chromosomal regions, while also supporting earlier association findings for these data. Understanding the genetic etiology of dental caries will allow designing personalized treatment plans based on an individual's genetic risk of disease.

CRISPLD2 variants including a C471T silent mutation may contribute to nonsyndromic cleft lip with or without cleft palate.

OBJECTIVE: To assess the association between nonsyndromic (NS) cleft lip with or without cleft palate (CL(P)) and single-nucleotide polymorphisms (SNPs) within the CRISPLD2 gene (cysteine-rich secretory protein LCCL domain containing 2). DESIGN: Four SNPs within the CRISPLD2 gene domain (rs1546124, rs8061351, rs2326398, rs4783099) were genotyped to test for association via family-based association methods. PARTICIPANTS: A total of 5826 individuals from 1331 families in which one or more family member is affected with CL(P). RESULTS: Evidence of association was seen for SNP rs1546124 in U.S. (p  =  .02) and Brazilian (p  =  .04) Caucasian cohorts. We also found association of SNP rs1546124 with cleft palate alone (CP) in South Americans (Guatemala and ECLAMC) and combined Hispanics (Guatemala, ECLAMC, and Texas Hispanics; p  =  .03 for both comparisons) and with both cleft lip with cleft palate (CLP; p  =  .04) and CL(P) (p  =  .02) in North Americans. Strong evidence of association was found for SNP rs2326398 with CP in Asian populations (p  =  .003) and with CL(P) in Hispanics (p  =  .03) and also with bilateral CL(P) in Brazilians (p  =  .004). In Brazilians, SNP rs8061351 showed association with cleft subgroups incomplete CL(P) (p  =  .004) and unilateral incomplete CL(P) (p  =  .003). Prediction of SNP functionality revealed that the C allele in the C471T silent mutation (overrepresented in cases with CL(P) presents two putative exonic splicing enhancer motifs and creates a binding site AP-2 alpha, a transcription factor involved in craniofacial development. CONCLUSIONS: Our results support the hypothesis that variants in the CRISPLD2 gene may be involved in the etiology of NS CL(P).

Follow-up association studies of chromosome region 9q and nonsyndromic cleft lip/palate.

Cleft lip/palate comprises a large fraction of all human birth defects, and is notable for its significant lifelong morbidity and complex etiology. Several studies have shown that genetic factors appear to play a significant role in the etiology of cleft lip/palate. Human chromosomal region 9q21 has been suggested in previous reports to contain putative cleft loci. Moreover, a specific region (9q22.3-34.1) was suggested to present a approximately 45% probability of harboring a cleft susceptibility gene. Fine mapping of 50 SNPs across the 9q22.3-34.11 region was performed to test for association with cleft lip/palate in families from United States, Spain, Turkey, Guatemala, and China. We performed family-based analyses and found evidence of association of cleft lip/palate with STOM (rs306796) in Guatemalan families (P = 0.004) and in all multiplex families pooled together (P = 0.002). This same SNP also showed borderline association in the US families (P = 0.04). Under a nominal value of 0.05, other SNPs also showed association with cleft lip/palate and cleft subgroups. SNPs in STOM and PTCH genes and nearby FOXE1 were further associated with cleft phenotypes in Guatemalan and Chinese families. Gene prioritization analysis revealed PTCH and STOM ranking among the top fourteen candidates for cleft lip/palate among 339 genes present in the region. Our results support the hypothesis that the 9q22.32-34.1 region harbors cleft susceptibility genes. Additional studies with other populations should focus on these loci to further investigate the participation of these genes in human clefting.

Genome scan, fine-mapping, and candidate gene analysis of non-syndromic cleft lip with or without cleft palate reveals phenotype-specific differences in linkage and association results.

OBJECTIVES: Non-syndromic orofacial clefts, i.e. cleft lip (CL) and cleft palate (CP), are among the most common birth defects. The goal of this study was to identify genomic regions and genes for CL with or without CP (CL/P). METHODS: We performed linkage analyses of a 10 cM genome scan in 820 multiplex CL/P families (6,565 individuals). Significant linkage results were followed by association analyses of 1,476 SNPs in candidate genes and regions, utilizing a weighted false discovery rate (wFDR) approach to control for multiple testing and incorporate the genome scan results. RESULTS: Significant (multipoint HLOD >or=3.2) or genome-wide-significant (HLOD >or=4.02) linkage results were found for regions 1q32, 2p13, 3q27-28, 9q21, 12p11, 14q21-24 and 16q24. SNPs in IRF6 (1q32) and in or near FOXE1 (9q21) reached formal genome-wide wFDR-adjusted significance. Further, results were phenotype dependent in that the IRF6 region results were most significant for families in which affected individuals have CL alone, and the FOXE1 region results were most significant in families in which some or all of the affected individuals have CL with CP. CONCLUSIONS: These results highlight the importance of careful phenotypic delineation in large samples of families for genetic analyses of complex, heterogeneous traits such as CL/P.

Genetic Predictors of Knee Pain in Persons With Mild to Moderate Osteoarthritis.

The purpose of this study was to examine genetic variability and knee pain in persons with osteoarthritis (OA). Seventy-five participants with medial compartment knee OA were recruited from a large Midwestern tertiary care center. Participants exhibited a mean age of 56.3 years; females comprised 61% of the sample. Measures of pain included subjective pain intensity at rest and with movement, cutaneous mechanical sensation and pain testing, heat pain threshold, and pressure pain threshold. Seventy-four participants were genotyped for 25 genetic variants across 15 candidate genes for central or peripheral pain pathways. Analysis suggests a role for four genes (EDNRA, COMT, BDRKB1, and IL1B) in several components of pain in persons with knee OA. The results from this study will help guide the development and evaluation of tailored strategies to decrease pain, improve function, and prevent the development of new chronic pain syndromes in older adults experiencing OA. [Research in Gerontological Nursing, xx(x), xx-xx.].

Effects of genotype on TENS effectiveness in controlling knee pain in persons with mild to moderate osteoarthritis.

BACKGROUND: This study examined the extent to which genetic variability modifies Transcutaneous Electrical Nerve Stimulation (TENS) effectiveness in osteoarthritic knee pain. METHODS: Seventy-five participants with knee osteoarthritis were randomly assigned to either: (a) High-frequency TENS, (b) Low-frequency TENS or (c) Transient Placebo TENS. Pain measures were collected pre- and post-treatment. Participants were genotyped on genes implicated in central or peripheral pain pathways: NGFB, NTRK1, EDNRA, EDNRB, EDN1, OPRM1, TAC1, TACR1, BDNF, BDKRB1, 5HTT, COMT, ESR2, IL6 and IL1B. Genetic association using linear regression modelling was performed separately for the transient placebo TENS subjects, and within the High-frequency TENS + Low-frequency TENS participants, including TENS level as a covariate. RESULTS: In the placebo group, SNPs rs165599 (COMT) was significantly associated with an increased heat pain threshold (ß = -1.87; p = .003) and rs6827096 (EDNRA) with an increased resting pain (ß = 2.68; p = .001). Within the treatment groups, TENS effectiveness was reduced by the SNP rs6537485 (EDNRA) minor allele in relationship to mechanical sensation (ß = 184.13; p = 5.5E-9). Individuals with the COMT rs4680 minor allele reported lowered pain at rest after TENS (ß = -42.30; p = .001), with a higher magnitude of pain reduction (28 unit difference) in the low-frequency TENS group compared to the high-frequency TENS group (ß = 28.37; p = .0004). CONCLUSIONS: EDNRA and COMT are implicated in osteoarthritic knee pain and provide a basis for tailoring TENS interventions according to individual characteristics. SIGNIFICANCE: Findings from this study demonstrate that genetic variation within the COMT and EDNRA genes influences the effectiveness of TENS, a non-pharmacologic pain-reduction intervention, in the context of osteoarthritic knee pain. Evidence such as this may contribute to risk models that provide a clinically useful tool for personalizing TENS interventions according to individual characteristics in order to best control pain and maximize functional status.

Aggressive periodontitis is likely influenced by a few small effect genes.

AIM: To evaluate the inheritance mode of aggressive periodontitis in a collection of families with a similar geographic origin. MATERIALS AND METHODS: Segregation analysis was performed in pedigree data from 74 families by the use of the SEGREG program of SAGE v.5.4.2. Homogeneous no transmission, homogeneous Mendelian transmission, homogeneous general transmission, semi-general transmission and heterogeneous general transmission models were tested assuming the prevalence of aggressive periodontitis as 1% and no deviations from Hardy-Weinberg equilibrium. The parameters of the model were estimated by the method of maximum likelihood, which provides the overall ln (likelihood), -2ln and the AIC (Akaike's score) for each model. The likelihood ratio test (LRT) was used to compare each model against a fully general model (p>0.05). RESULTS: The most parsimonious mode of inheritance was the semi-general transmission model that allows the heterozygote transmission probability to vary. CONCLUSION: This result provides strong support for the hypothesis that genetic factors play a role in aggressive periodontitis and that a few loci, each with relatively small effects, contribute to aggressive periodontitis, with or without interaction with environmental factors.

MLIP: using multiple processors to compute the posterior probability of linkage.

BACKGROUND: Localization of complex traits by genetic linkage analysis may involve exploration of a vast multidimensional parameter space. The posterior probability of linkage (PPL), a class of statistics for complex trait genetic mapping in humans, is designed to model the trait model complexity represented by the multidimensional parameter space in a mathematically rigorous fashion. However, the method requires the evaluation of integrals with no functional form, making it difficult to compute, and thus further test, develop and apply. This paper describes MLIP, a multiprocessor two-point genetic linkage analysis system that supports statistical calculations, such as the PPL, based on the full parameter space implicit in the linkage likelihood. RESULTS: The fundamental question we address here is whether the use of additional processors effectively reduces total computation time for a PPL calculation. We use a variety of data - both simulated and real - to explore the question "how close can we get?" to linear speedup. Empirical results of our study show that MLIP does significantly speed up two-point log-likelihood ratio calculations over a grid space of model parameters. CONCLUSION: Observed performance of the program is dependent on characteristics of the data including granularity of the parameter grid space being explored and pedigree size and structure. While work continues to further optimize performance, the current version of the program can already be used to efficiently compute the PPL. Thanks to MLIP, full multidimensional genome scans are now routinely being completed at our centers with runtimes on the order of days, not months or years.

Identifying genetic risk loci for diabetic complications and showing evidence for heterogeneity of type 1 diabetes based on complications risk.

There is a growing body of evidence suggesting that type 1 diabetes (T1D) is a genetically heterogeneous disease. However, the extent of this heterogeneity, and what observations may distinguish different forms, is unclear. One indicator may be T1D-related microvascular complications (MVCs), which are familial, but occur in some families, and not others. We tested the hypothesis that T1D plus MVC is genetically distinct from T1D without MCV. We studied 415 families (2,462 individuals, 896 with T1D) using genome-wide linkage analysis, comparing families with and without MVC. We also tested for interaction between identified loci and alleles at the HLA-DRB1 locus. We found significant linkage scores at 1p36.12, 1q32.1, 8q21.3, 12p11.21 and 22q11.21. In all regions except 1p36.12, linkage scores differed between MVC-based phenotype groups, suggesting that families with MVCs express different genetic influences than those without. Our linkage results also suggested gene-gene interaction between the above putative loci and the HLA region; HLA-based strata produced significantly increased linkage scores in some strata, but not others within a phenotype group. We conclude that families with type 1 diabetes plus MVCs are genetically different from those with diabetes alone.

Vitamin D metabolic loci and preeclampsia risk in multi-ethnic pregnant women.

Allelic variants in vitamin D metabolism genes may increase the risk of preeclampsia, but few studies have systematically tested this hypothesis. Our objective was to evaluate the relationship between maternal allelic variants in three vitamin D metabolism genes and risk of preeclampsia. Samples were from two case-control studies of pregnant women who delivered in Pittsburgh, PA from 1999 to 2010 and twelve recruiting sites across the United States from 1959 to 1965. Single-nucleotide polymorphisms (SNPs) were genotyped 50 kilobases up- and down-stream in three genes (VDR, GC, and CYP27B1) in the samples from both studies, for a total of 744 preeclampsia cases and 2411 controls. Using multivariable logistic regression, we estimated the associations between allelic variation in each locus and preeclampsia risk by maternal race and study. Meta-analysis was used to estimate the association across race-study groups for each SNP. Minor allele of a noncoding region of the VDR gene was significantly associated with preeclampsia risk, which was verified in the meta-analysis [odds ratio (OR), 95% confidence intervals (CI)] after adjusting for multiple comparisons [rs12831006:1.5 (1.2, 2.0), P < 0.0001]. The meta-analysis identified associations for one intron GC variant [rs843010:1.4 (1.1, 1.9) P < 0.05] and two variants of the flanking region of GC [rs842991:1.5 (1.1, 2.0) P < 0.05; rs16846876:0.75 (0.58, 0.98) P < 0.05]. There were no statistically significant associations for CYP27B1 SNPs. Our results provide additional support for a biological role of vitamin D in preeclampsia.

Vitamin D metabolic loci and vitamin D status in Black and White pregnant women.

BACKGROUND: Several candidate genes and genome wide association studies have reported significant associations between vitamin D metabolism genes and 25-hydroxyvitamin D. Few studies have examined these relationships in pregnancy. OBJECTIVE: We evaluated the relationship between maternal allelic variants in three vitamin D metabolism genes and 25-hydroxyvitamin D (25(OH)D) concentration in pregnancy. STUDY DESIGN: In two case-control studies, samples were drawn from women who delivered at Magee Womens Hospital in Pittsburgh, PA from 1999 to 2010 and twelve recruiting sites across the United States from 1959 to 65. For 882 Black and 1796 White pregnant women from these studies, 25(OH)D concentration was measured and single nucleotide polymorphisms (SNPs) were genotyped 50 kilobases up- and down-stream in three genes (VDR, GC, and CYP27B1). Using multivariable linear regression, we estimated the associations between allelic variation of each locus and log-transformed 25(OH)D concentration separately by race and study group. Meta-analysis was used to estimate the association across the four groups for each SNP. RESULTS: Minor alleles of several variants in VDR, GC, and CYP27B1 were associated with differences in log-transformed 25(OH)D concentration compared to the corresponding major alleles [beta, 95% confidence intervals (CI)]. The meta-analysis confirmed the associations for differences in log-transformed 25(OH)D by allelic loci for one intron VDR variant [rs2853559 0.08 (0.02, 0.13), p<0.01] and a variant in the GC flanking region [rs13150174: 0.04 (0.02, 0.07), p<0.01], and a GC missense mutation [rs7041 0.05 (0.01, 0.09), p<0.01]. The meta-analysis also revealed possible associations for SNPs in linkage disequilibrium with variants in the VDR 3-prime untranslated region, another GC missense variant (rs4588), and a variant of the 3-prime untranslated region of CYP27B1. CONCLUSION: We observed associations between VDR, GC, and CYP27B1 variants and maternal 25-hydroxyvitamin D concentration. Our results provide additional support for a possible role of genetic variation in vitamin D metabolism genes on vitamin D status during pregnancy.

Caries Experience Differs between Females and Males across Age Groups in Northern Appalachia.

Sex disparities in dental caries have been observed across many populations, with females typically exhibiting higher prevalence and more affected teeth. In this study we assessed the sex disparities in two Northern Appalachian populations from West Virginia (WV, N = 1997) and Pennsylvania (PA, N = 1080) by comparing caries indices between males and females across four phases of dental development: primary dentition in children aged 1-5 years, mixed dentition in children aged 6-11 years, permanent dentition in adolescents aged 12-17 years, and permanent dentition in adults aged 18-59 years. No significant sex differences were observed for children aged 1-5 years. Contrary to national and international trends, WV girls aged 6-11 years had 1.5 fewer affected teeth than boys (p < 0.001). However, by ages 12-17, caries indices in the WV girls matched those in boys. In both WV and PA adults, women and men had similar total counts of affected teeth (i.e., DMFT), although women had more dental restorations (p < 0.001) and men had more current decay (p < 0.001). These results suggest that in some Appalachian populations, young girls benefit from protection against caries that is lost during adolescence and that adult women utilize dental health care to a greater degree than men.

Oral Health in a Sample of Pregnant Women from Northern Appalachia (2011-2015).

Background. Chronic poor oral health has a high prevalence in Appalachia, a large region in the eastern USA. The Center for Oral Health Research in Appalachia (COHRA) has been enrolling pregnant women and their babies since 2011 in the COHRA2 study of genetic, microbial, and environmental factors involved in oral health in Northern Appalachia. Methods. The COHRA2 protocol is presented in detail, including inclusion criteria (healthy, adult, pregnant, US Caucasian, English speaking, and nonimmunocompromised women), recruiting (two sites: Pittsburgh, Pennsylvania, and West Virginia, USA), assessments (demographic, medical, dental, psychosocial/behavioral, and oral microbial samples and DNA), timelines (longitudinal from pregnancy to young childhood), quality control, and retention rates. Results. Preliminary oral health and demographic data are presented in 727 pregnant women, half from the greater Pittsburgh region and half from West Virginia. Despite similar tooth brushing and flossing habits, COHRA2 women in West Virginia have significantly worse oral health than the Pittsburgh sample. Women from Pittsburgh are older and more educated and have less unemployment than the West Virginia sample. Conclusions. We observed different prevalence of oral health and demographic variables between pregnant women from West Virginia (primarily rural) and Pittsburgh (primarily urban). These observations suggest site-specific differences within Northern Appalachia that warrant future studies.

The PDGF-C regulatory region SNP rs28999109 decreases promoter transcriptional activity and is associated with CL/P.

Human linkage and association studies suggest a gene(s) for nonsyndromic cleft lip with or without cleft palate (CL/P) on chromosome 4q31-q32 at or near the platelet-derived growth factor-C (PDGF-C) locus. The mouse Pdgfc(-/-) knockout shows that PDGF-C is essential for palatogenesis. To evaluate the role of PDGF-C in human clefting, we performed sequence analysis and SNP genotyping using 1048 multiplex CL/P families and 1000 case-control samples from multiple geographic origins. No coding region mutations were identified, but a novel -986 C>T SNP (rs28999109) was significantly associated with CL/P (P=0.01) in cases from Chinese families yielding evidence of linkage to 4q31-q32. Significant or near-significant association was also seen for this and several other PDGF-C SNPs in families from the United States, Spain, India, Turkey, China, and Colombia, whereas no association was seen in families from the Philippines, and Guatemala, and case-controls from Brazil. The -986T allele abolished six overlapping potential transcription regulatory motifs. Transfection assays of PDGF-C promoter reporter constructs show that the -986T allele is associated with a significant decrease (up to 80%) of PDGF-C gene promoter activity. This functional polymorphism acting on a susceptible genetic background may represent a component of human CL/P etiology.

Disruption of an AP-2alpha binding site in an IRF6 enhancer is associated with cleft lip.

Previously we have shown that nonsyndromic cleft lip with or without cleft palate (NSCL/P) is strongly associated with SNPs in IRF6 (interferon regulatory factor 6). Here, we use multispecies sequence comparisons to identify a common SNP (rs642961, G>A) in a newly identified IRF6 enhancer. The A allele is significantly overtransmitted (P = 1 x 10(-11)) in families with NSCL/P, in particular those with cleft lip but not cleft palate. Further, there is a dosage effect of the A allele, with a relative risk for cleft lip of 1.68 for the AG genotype and 2.40 for the AA genotype. EMSA and ChIP assays demonstrate that the risk allele disrupts the binding site of transcription factor AP-2alpha and expression analysis in the mouse localizes the enhancer activity to craniofacial and limb structures. Our findings place IRF6 and AP-2alpha in the same developmental pathway and identify a high-frequency variant in a regulatory element contributing substantially to a common, complex disorder.

Model selection and Bayesian methods in statistical genetics: summary of group 11 contributions to Genetic Analysis Workshop 15.

The research presented in group 11 of the Genetic Analysis Workshop 15 (GAW15) falls into two major themes: Model selection approaches for gene mapping (both Bayesian and Frequentist); and other Bayesian methods. These methods either allow relaxation of some of the common assumptions, such as mode of inheritance, for studying complicated genetic systems, or allow incorporation of additional information into the model. Over half of the groups applied model selection methods on all three data sets, using models in which genetic markers were used as predictors for linkage, phenotype expression, or transmission to an affected offspring. Most groups employed variations of Stochastic Search Variable Selection as the model selection method of choice. A brief review of this class of methods is given in this summary paper, followed by highlights of other methods and overall summaries of each contribution to the GAW15 presentation group 11. These group contributions exhibit the value of framing genetic problems in terms of model selection, and highlight the impact of variable selection for gene mapping.