A powerful method for pleiotropic analysis under composite null hypothesis identifies novel shared loci between Type 2 Diabetes and Prostate Cancer

There is increasing evidence that pleiotropy, the association of multiple traits with the same genetic variants/loci, is a very common phenomenon. Cross-phenotype association tests are often used to jointly analyze multiple traits from a genome-wide association study (GWAS). The underlying methods, however, are often designed to test the global null hypothesis that there is no association of a genetic variant with any of the traits, the rejection of which does not implicate pleiotropy. In this article, we propose a new statistical approach, PLACO, for specifically detecting pleiotropic loci between two traits by considering an underlying composite null hypothesis that a variant is associated with none or only one of the traits. We propose testing the null hypothesis based on the product of the Z-statistics of the genetic variants across two studies and derive a null distribution of the test statistic in the form of a mixture distribution that allows for fractions of variants to be associated with none or only one of the traits. We borrow approaches from the statistical literature on mediation analysis that allow asymptotic approximation of the null distribution avoiding estimation of nuisance parameters related to mixture proportions and variance components. Simulation studies demonstrate that the proposed method can maintain type I error and can achieve major power gain over alternative simpler methods that are typically used for testing pleiotropy. PLACO allows correlation in summary statistics between studies that may arise due to sharing of controls between disease traits. Application of PLACO to publicly available summary data from two large case-control GWAS of Type 2 Diabetes and of Prostate Cancer implicated a number of novel shared genetic regions: 3q23 (ZBTB38), 6q25.3 (RGS17), 9p22.1 (HAUS6), 9p13.3 (UBAP2), 11p11.2 (RAPSN), 14q12 (AKAP6), 15q15 (KNL1) and 18q23 (ZNF236).     

Genome-wide association meta-analysis identifies pleiotropic risk loci for aerodigestive squamous cell cancers

Squamous cell carcinomas (SqCC) of the aerodigestive tract have similar etiological risk factors. Although genetic risk variants for individual cancers have been identified, an agnostic, genome-wide search for shared genetic susceptibility has not been performed. To identify novel and pleotropic SqCC risk variants, we performed a meta-analysis of GWAS data on lung SqCC (LuSqCC), oro/pharyngeal SqCC (OSqCC), laryngeal SqCC (LaSqCC) and esophageal SqCC (ESqCC) cancers, totaling 13,887 cases and 61,961 controls of European ancestry. We identified one novel genome-wide significant (P_meta<5x10^-8) aerodigestive SqCC susceptibility loci in the 2q33.1 region (rs56321285, TMEM273). Additionally, three previously unknown loci reached suggestive significance (P_meta<5x10^-7): 1q32.1 (rs12133735, near MDM4), 5q31.2 (rs13181561, TMEM173) and 19p13.11 (rs61494113, ABHD8). Multiple previously identified loci for aerodigestive SqCC also showed evidence of pleiotropy in at least another SqCC site, these include: 4q23 (ADH1B), 6p21.33 (STK19), 6p21.32 (HLA-DQB1), 9p21.33 (CDKN2B-AS1) and 13q13.1(BRCA2). Gene-based association and gene set enrichment identified a set of 48 SqCC-related genes rel to DNA damage and epigenetic regulation pathways. Our study highlights the importance of cross-cancer analyses to identify pleiotropic risk loci of histology-related cancers arising at distinct anatomical sites.     

Population-Based Study of Birth Prevalence and Factors Associated with Cleft Lip and/or Palate in Taiwan 2002–2009

Background

Facial cleft deformities, including cleft lip with or without cleft palate (CL/P) and cleft palate (CP), are common congenital birth anomalies, especially in Asia. This study aimed to analyze the prevalence of CL/P and CP and to identify associated factors in Taiwan.

Methods

This population-based epidemiological study retrospectively analyzed birth data obtained from the Department of Health in Taiwan for years 2002–2009. Frequency distribution, percentages and related predictors were investigated, and findings were presented by types of cleft deformities. Logistic regression analysis was performed to identify factors associated with cleft deformities.

Results

Overall prevalence of cleft deformities among 1,705,192 births was 0.1% for CL/P and 0.04% for CP over the 8-year study period. Higher prevalence of CL/P or CP was observed with multiple pregnancies, being male for CL/P, being female for CP, gestational age ≤37 weeks and lower birth weight (<1.5 kg). Both CL/P and CP were significantly associated with gestational age <37 weeks and birth weight<1.5 kg (all P <0.0001). CL/P was significantly associated with multiple parities (P = 0.0004–0.002). Male newborns and female newborns were significantly associated with CL/P and CP, respectively (both P<0.0001).

Conclusions

Overall prevalence for congenital cleft deformities in study subjects was 0.1%, in keeping with high rates in Asia. Results suggest the need for awareness and early identification of those at high risk for cleft deformities, including newborns with gestational age <37 weeks, weighing <1.5 kg at birth and women with multiple parities, as a potential strategy to counter long-term adverse effects on speech and language in this population.     

Evaluating Genome-Wide Association Study-Identified Breast Cancer Risk Variants in African-American Women

Genome-wide association studies (GWAS), conducted mostly in European or Asian descendants, have identified approximately 67 genetic susceptibility loci for breast cancer. Given the large differences in genetic architecture between the African-ancestry genome and genomes of Asians and Europeans, it is important to investigate these loci in African-ancestry populations. We evaluated index SNPs in all 67 breast cancer susceptibility loci identified to date in our study including up to 3,300 African-American women (1,231 cases and 2,069 controls), recruited in the Southern Community Cohort Study (SCCS) and the Nashville Breast Health Study (NBHS). Seven SNPs were statistically significant (P≤0.05) with the risk of overall breast cancer in the same direction as previously reported: rs10069690 (5p15/TERT), rs999737 (14q24/RAD51L1), rs13387042 (2q35/TNP1), rs1219648 (10q26/FGFR2), rs8170 (19p13/BABAM1), rs17817449 (16q12/FTO), and rs13329835 (16q23/DYL2). A marginally significant association (P<0.10) was found for three additional SNPs: rs1045485 (2q33/CASP8), rs4849887 (2q14/INHBB), and rs4808801 (19p13/ELL). Three additional SNPs, including rs1011970 (9p21/CDKN2A/2B), rs941764 (14q32/CCDC88C), and rs17529111 (6q14/FAM46A), showed a significant association in analyses conducted by breast cancer subtype. The risk of breast cancer was elevated with an increasing number of risk variants, as measured by quintile of the genetic risk score, from 1.00 (reference), to 1.75 (1.30–2.37), 1.56 (1.15–2.11), 2.02 (1.50–2.74) and 2.63 (1.96–3.52), respectively, (P = 7.8×10^–10). Results from this study highlight the need for large genetic studies in AAs to identify risk variants impacting this population.     

Medical sequencing of candidate genes for nonsyndromic cleft lip and palate

Nonsyndromic or isolated cleft lip with or without cleft palate (CL/P) occurs in wide geographic distribution with an average birth prevalence of 1/700. We used direct sequencing as an approach to study candidate genes for CL/P. We report here the results of sequencing on 20 candidate genes for clefts in 184 cases with CL/P selected with an emphasis on severity and positive family history. Genes were selected based on expression patterns, animal models, and/or role in known human clefting syndromes. For seven genes with identified coding mutations that are potentially etiologic, we performed linkage disequilibrium studies as well in 501 family triads (affected child/mother/father). The recently reported MSX1 P147Q mutation was also studied in an additional 1,098 cleft cases. Selected missense mutations were screened in 1,064 controls from unrelated individuals on the Centre d'Étude du Polymorphisme Humain (CEPH) diversity cell line panel. Our aggregate data suggest that point mutations in these candidate genes are likely to contribute to 6% of isolated clefts, particularly those with more severe phenotypes (bilateral cleft of the lip with cleft palate). Additional cases, possibly due to microdeletions or isodisomy, were also detected and may contribute to clefts as well. Sequence analysis alone suggests that point mutations in FOXE1, GLI2, JAG2, LHX8, MSX1, MSX2, SATB2, SKI, SPRY2, and TBX10 may be rare causes of isolated cleft lip with or without cleft palate, and the linkage disequilibrium data support a larger, as yet unspecified, role for variants in or near MSX2, JAG2, and SKI. This study also illustrates the need to test large numbers of controls to distinguish rare polymorphic variants and prioritize functional studies for rare point mutations.     

Genome-wide Comparative Analysis of Atopic Dermatitis and Psoriasis Gives Insight into Opposing Genetic Mechanisms

Atopic dermatitis and psoriasis are the two most common immune-mediated inflammatory disorders affecting the skin. Genome-wide studies demonstrate a high degree of genetic overlap, but these diseases have mutually exclusive clinical phenotypes and opposing immune mechanisms. Despite their prevalence, atopic dermatitis and psoriasis very rarely co-occur within one individual. By utilizing genome-wide association study and ImmunoChip data from >19,000 individuals and methodologies developed from meta-analysis, we have identified opposing risk alleles at shared loci as well as independent disease-specific loci within the epidermal differentiation complex (chromosome 1q21.3), the Th2 locus control region (chromosome 5q31.1), and the major histocompatibility complex (chromosome 6p21–22). We further identified previously unreported pleiotropic alleles with opposing effects on atopic dermatitis and psoriasis risk in PRKRA and ANXA6/TNIP1. In contrast, there was no evidence for shared loci with effects operating in the same direction on both diseases. Our results show that atopic dermatitis and psoriasis have distinct genetic mechanisms with opposing effects in shared pathways influencing epidermal differentiation and immune response. The statistical analysis methods developed in the conduct of this study have produced additional insight from previously published data sets. The approach is likely to be applicable to the investigation of the genetic basis of other complex traits with overlapping and distinct clinical features.     

Linkage analysis in a large Brazilian family with van der Woude syndrome suggests the existence of a susceptibility locus for cleft palate at 17p11.2-11.1.

van der Woude syndrome (VWS), which has been mapped to 1q32-41, is characterized by pits and/or sinuses of the lower lip, cleft lip/palate (CL/P), cleft palate (CP), bifid uvula, and hypodontia (H). The expression of VWS, which has incomplete penetrance, is highly variable. Both the occurrence of CL/P and CP within the same genealogy and a recurrence risk <40% for CP among descendants with VWS have suggested that the development of clefts in this syndrome is influenced by modifying genes at other loci. To test this hypothesis, we have conducted linkage analysis in a large Brazilian kindred with VWS, considering as affected the individuals with CP, regardless of whether it is associated with other clinical signs of VWS. Our results suggest that a gene at 17p11.2-11.1, together with the VWS gene at 1p32-41, enhances the probability of CP in an individual carrying the two at-risk genes. If this hypothesis is confirmed in other VWS pedigrees, it will represent one of the first examples of a gene, mapped through linkage analysis, which modifies the expression of a major gene. It will also have important implications for genetic counseling, particularly for more accurately predicting recurrence risks of clefts among the offspring of patients with VWS.     

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.     

Etiology and genetic factors in clefts of lip and/or palate reported at children’s hospital,Lahore, Pakistan

The etiology of cleft lip (CL) and/or cleft palate (CP) has been extensively studied in industrialized countries and is suggested to be heterogeneous with increasing evidence that both genetic and environmental factors are operating. To evaluate this assertion in a developing country like Pakistan, a case finding cross-sectional study was completed from 1^st July 2010 to 31^st May 2011 for 100 cases of CL and/or CP referred to the Genetic Clinic of the Children’s Hospital, Lahore, Pakistan. A clinical examination followed by necessary diagnostic work-up was completed for each case. The cause of CL and/or CP was clear in 18% of the children (n = 18). Environmental causes were found in 6 children (four mothers developed hyperthermia during the 2^nd month of gestation, one mother was diabetic, and one mother was a known case of epilepsy and took sodium valproate throughout her pregnancy). Six children were suffering from known genetic malformation syndromes (each with Jarcho-Levin syndrome, Oral-Facial-Digital syndrome type XI, Oral-Duplication syndrome, Kabuki syndrome, Fronto-nasal dysplasia and Nager syndrome). Novel chromosomal aberrations were identified in 2 children. In 82% of the children (n = 82) the cause of oro-facial clefts remained unknown. Impact of gender and consanguinity on the development of CL and/or CP was also studied. Prevalence of CP was significantly more among female children as compared to that in males (P < 0.05). Associated anomalies were present in 18% of the cases, anomalies of the craniofacial region being the most common. These findings were compared with regional and international studies.     

Genome-Wide Association Studies in Dogs and Humans Identify ADAMTS20 as a Risk Variant for Cleft Lip and Palate

Cleft lip with or without cleft palate (CL/P) is the most commonly occurring craniofacial birth defect. We provide insight into the genetic etiology of this birth defect by performing genome-wide association studies in two species: dogs and humans. In the dog, a genome-wide association study of 7 CL/P cases and 112 controls from the Nova Scotia Duck Tolling Retriever (NSDTR) breed identified a significantly associated region on canine chromosome 27 (unadjusted p=1.1 x 10^-13; adjusted p= 2.2 x 10^-3). Further analysis in NSDTR families and additional full sibling cases identified a 1.44 Mb homozygous haplotype (chromosome 27: 9.29 – 10.73 Mb) segregating with a more complex phenotype of cleft lip, cleft palate, and syndactyly (CLPS) in 13 cases. Whole-genome sequencing of 3 CLPS cases and 4 controls at 15X coverage led to the discovery of a frameshift mutation within ADAMTS20 (c.1360_1361delAA (p.Lys453Ilefs*3)), which segregated concordant with the phenotype. In a parallel study in humans, a family-based association analysis (DFAM) of 125 CL/P cases, 420 unaffected relatives, and 392 controls from a Guatemalan cohort, identified a suggestive association (rs10785430; p =2.67 x 10^-6) with the same gene, ADAMTS20. Sequencing of cases from the Guatemalan cohort was unable to identify a causative mutation within the coding region of ADAMTS20, but four coding variants were found in additional cases of CL/P. In summary, this study provides genetic evidence for a role of ADAMTS20 in CL/P development in dogs and as a candidate gene for CL/P development in humans.     

Genetic risk factors for nonsyndromic cleft lip with or without cleft palate in a Mesoamerican population: Evidence for IRF6 and variants at 8q24 and 10q25

INTRODUCTION: Nonsyndromic cleft lip with or without cleft palate (NSCL/P) is one of the most common of all birth defects. NSCL/P has a multifactorial etiology that includes both genetic and environmental factors. The IRF6 gene and three further susceptibility loci at 8q24, 10q25, and 17q22, which were identified by a recent genome‐wide association scan (GWAS), are confirmed genetic risk factors for NSCL/P in patients of European descent. METHODS: A case‐control association study was performed to investigate whether these four risk loci contribute to NSCL/P in a Mesoamerican population using four single nucleotide polymorphisms to represent IRF6 and the three novel susceptibility loci. A total of 149 NSCL/P patients and 303 controls of Mayan origin were included. RESULTS: Single marker analysis revealed a significant association between NSCL/P and risk variants in IRF6 and the 8q24 and 10q25 loci. In contrast to previous findings, the association at the 8q24 locus was driven solely by homozygote carriers of the risk allele. This suggests that this locus might act in a recessive manner in the Mayan population. No evidence for association was found at the 17q22 locus. This may have been attributable to the limited power of the sample. CONCLUSION: These results suggest that IRF6 and the 10q25 and 8q24 loci confer a risk for the development of NSCL/P in persons of Mayan origin. Birth Defects Research (Part A), 2010. © 2010 Wiley‐Liss, Inc.     

PADI4 and HLA-DRB1 Are Genetic Risks for Radiographic Progression in RA Patients,Independent of ACPA Status: Results from the IORRA Cohort Study

Introduction

Rheumatoid arthritis (RA) is a systemic, chronic inflammatory disease influenced by both genetic and environmental factors, leading to joint destruction and functional impairment. Recently, a large-scaled GWAS meta-analysis using more than 37,000 Japanese samples were conducted and 13 RA susceptibility loci were identified. However, it is not clear whether these loci have significant impact on joint destruction or not. This is the first study focused on the 13 loci to investigate independent genetic risk factors for radiographic progression in the first five years from onset of RA.

Methods

Sharp/van der Heijde score of hands at 5-year disease duration, which represents joint damage, were measured retrospectively and used as an outcome variable in 865 Japanese RA patients. Genetic factors regarded as putative risk factors were RA-susceptible polymorphisms identified by the Japanese GWAS meta-analysis, including HLA-DRB1 (shared epitope, SE), rs2240340 (PADI4), rs2230926 (TNFAIP3), rs3093024 (CCR6), rs11900673 (B3GNT2), rs2867461 (ANXA3), rs657075 (CSF2), rs12529514 (CD83), rs2233434 (NFKBIE), rs10821944 (ARID5B), rs3781913 (PDE2A-ARAP1), rs2841277 (PLD4) and rs2847297 (PTPN2). These putative genetic risk factors were assessed by a stepwise multiple regression analysis adjusted for possible non-genetic risk factors: autoantibody positivity (anti-citrullinated peptide antibody [ACPA] and rheumatoid factor), history of smoking, gender and age at disease onset.

Results

The number of SE alleles (P = 0.002) and risk alleles of peptidyl arginine deiminase type IV gene (PADI4, P = 0.04) had significant impact on progressive joint destruction, as well as following non-genetic factors: ACPA positive (P = 0.0006), female sex (P = 0.006) and younger age of onset (P = 0.02).

Conclusions

In the present study, we found that PADI4 risk allele and HLA-DRB1 shared epitope are independent genetic risks for radiographic progression in Japanese rheumatoid arthritis patients. The results of this study give important knowledge of the risks on progressive joint damage in RA patients.     

Analysis of candidate genes on chromosome 2 in oral cleft case-parent trios from three populations

Isolated oral clefts, including cleft lip with/without cleft palate (CL/P) and cleft palate (CP), have a complex and heterogeneous etiology. Case-parent trios from three populations were used to study genes spanning chromosome 2, where single nucleotide polymorphic (SNP) markers were analyzed individually and as haplotypes. Case-parent trios from three populations (74 from Maryland, 64 from Singapore and 95 from Taiwan) were genotyped for 962 SNPs in 104 genes on chromosome 2, including two well-recognized candidate genes: TGFA and SATB2. Individual SNPs and haplotypes (in sliding windows of 2–5 SNPs) were used to test for linkage and disequilibrium separately in CL/P and CP trios. A novel candidate gene (ZNF533) showed consistent evidence of linkage and disequilibrium in all three populations for both CL/P and CP. SNPs in key regions of ZNF533 showed considerable variability in estimated genotypic odds ratios and their significance, suggesting allelic heterogeneity. Haplotype frequencies for regions of ZNF533 were estimated and used to partition genetic variance into among-and within-population components. Wright’s fixation index, a measure of genetic diversity, showed little difference between Singapore and Taiwan compared with Maryland. The tensin-1 gene (TNS1) also showed evidence of linkage and disequilibrium among both CL/P and CP trios in all three populations, albeit at a lower level of significance. Additional genes (VAX2, GLI2, ZHFX1B on 2p; WNT6–WNT10A and COL4A3–COL4A4 on 2q) showed consistent evidence of linkage and disequilibrium only among CL/P trios in all three populations, and TGFA showed significant evidence in two of three populations.Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.     

Genome-Wide Association Study of Celiac Disease in North America Confirms FRMD4B as New Celiac Locus

We performed a genome-wide association study (GWAS) of 1550 North American celiac disease cases and 3084 controls. Twelve SNPs, distributed across four regions (3p21.31, 4q27, 6q15, 6q25), were significantly associated with disease (p-value <1.0×10⁻⁷), and a further seven SNPs, across four additional regions (1q24.3, 10p15.1, 6q22.31, 17q21.32) had suggestive evidence (1.0×10⁻⁷ < p-value < 1.0×10⁻⁶). This study replicated a previous suggestive association within FRMD4B (3p14.1), confirming it as a celiac disease locus. All four regions with significant associations and two regions with suggestive results (1q24.3, 10p15.1) were known disease loci. The 6q22.31 and 10p11.23 regions were not replicated. A total of 410 SNPs distributed across the eight significant and suggestive regions were tested for association with dermatitis herpetiformis and microscopic colitis. Preliminary, suggestive statistical evidence for association with the two traits was found at chromosomes 3p21.31, 6q15, 6q25, 1q24.3 and 10p11.23, with future studies being required to validate the reported associations.     

Interaction between IRF6 and TGFA Genes Contribute to the Risk of Nonsyndromic Cleft Lip/Palate

Previous evidence from tooth agenesis studies suggested IRF6 and TGFA interact. Since tooth agenesis is commonly found in individuals with cleft lip/palate (CL/P), we used four large cohorts to evaluate if IRF6 and TGFA interaction contributes to CL/P. Markers within and flanking IRF6 and TGFA genes were tested using Taqman or SYBR green chemistries for case-control analyses in 1,000 Brazilian individuals. We looked for evidence of gene-gene interaction between IRF6 and TGFA by testing if markers associated with CL/P were overtransmitted together in the case-control Brazilian dataset and in the additional family datasets. Genotypes for an additional 142 case-parent trios from South America drawn from the Latin American Collaborative Study of Congenital Malformations (ECLAMC), 154 cases from Latvia, and 8,717 individuals from several cohorts were available for replication of tests for interaction. Tgfa and Irf6 expression at critical stages during palatogenesis was analyzed in wild type and Irf6 knockout mice. Markers in and near IRF6 and TGFA were associated with CL/P in the Brazilian cohort (p<10⁻⁶). IRF6 was also associated with cleft palate (CP) with impaction of permanent teeth (p<10⁻⁶). Statistical evidence of interaction between IRF6 and TGFA was found in all data sets (p = 0.013 for Brazilians; p = 0.046 for ECLAMC; p = 10⁻⁶ for Latvians, and p = 0.003 for the 8,717 individuals). Tgfa was not expressed in the palatal tissues of Irf6 knockout mice. IRF6 and TGFA contribute to subsets of CL/P with specific dental anomalies. Moreover, this potential IRF6-TGFA interaction may account for as much as 1% to 10% of CL/P cases. The Irf6-knockout model further supports the evidence of IRF6-TGFA interaction found in humans.     

Evaluating eight newly identified susceptibility loci for nonsyndromic cleft lip with or without cleft palate in a Mesoamerican population

Background: Nonsyndromic cleft lip with or without cleft palate (NSCL/P) is among the most frequently occurring congenital malformations worldwide. The number of genetic loci identified as being involved in NSCL/P etiology was recently increased by a large genome‐wide meta‐analysis of European and Asian samples. This meta‐analysis confirmed all six previously recognized genetic susceptibility loci and identified six novel ones. Methods: To investigate which of these 12 loci contribute to NSCL/P risk in an independent sample of distinct ethnicity, we performed a case–control association analysis in a sample of the Mesoamerican population. A total of 153 individuals with NSCL/P (cases) and 337 unaffected controls were included. Top single‐nucleotide polymorphisms (SNPs) at 8 of the 12 loci (1p22.1, 1p36, 2p21, 3p11.1, 8q21.3, 13q31.1, 15q22, and 20q12) were analyzed using mass spectroscopy and restriction‐length‐fragment polymorphism analyses. In a previous study, we had analyzed the remaining four NSCL/P susceptibility regions (IRF6, 8q24, 10q25, and 17q22) in the same sample. Results: Single‐marker association analyses applying allelic, dominant, and recessive models revealed nominal significant associations for four of the eight loci, with two additional loci showing at least a trend of association in the hypothesized direction. Conclusion: In combination with results from our previous study using the same sample, our data suggest that the majority of the known NSCL/P susceptibility regions identified to date also confer risk for this malformation in the Mesoamerican population. Birth Defects Research (Part A) 100:43–47, 2014. © 2013 Wiley Periodicals, Inc.     

BCL3 gene role in facial morphology

BACKGROUND: Cleft lip (CL) with or without palate (CLP) and isolated cleft palate (CP) are etiologically complex diseases with interactions among various environmental and genetic factors. The aim of the current study was to identify association with genetic markers and phenotypic craniofacial data in patients with CL/CLP/CP parents. METHODS: Posteroanterior and lateral digital radiographs of the cranium were obtained from 74 parents of patients with CL/CLP/CP. One hundred seventy‐three patients with CL/CLP/CP and 190 controls were enrolled in the study for the association test. Five genetic markers of the IRF6 gene and 14 markers of the 19q13 locus were genotyped. Linear regression analysis was performed for the relationship of cephalometric measurements with genotype data adjusted for age, gender, and cleft type. Chi‐square and transmission disequilibrium tests were performed to evaluate differences in alleles of the BCL3 gene. Positive findings were replicated in an independent sample (n = 95) of patients with CL/CLP/CP parents. RESULTS: Genetic markers of the BCL3 gene at 19q13, rs7257231, and rs1979377 in the familial association test and rs10401176 in the case‐control association test, were associated with craniofacial phenotype. Carriers of BCL3 allele rs7257231T had longer posterior cranial bases than noncarriers (p_adjusted = 0.0028), and in the familial‐based association test showed the statistically strongest relationship (p_adjusted = 0.05) to phenotype. Relation of rs7257231 to facial formation was confirmed in the replication group (p = 0.0024). CONCLUSIONS: The results indicate that BCL3, which has functions related to cell adhesion and whose downregulation can cause disruption of ectodermal development, is likely to be important in facial formation. Birth Defects Research (Part A), 2012. © 2012 Wiley Periodicals, Inc.     

Genome-Wide Association Study of Breast Cancer in the Japanese Population

Breast cancer is the most common malignancy among women in worldwide including Japan. Several studies have identified common genetic variants to be associated with the risk of breast cancer. Due to the complex linkage disequilibrium structure and various environmental exposures in different populations, it is essential to identify variants associated with breast cancer in each population, which subsequently facilitate the better understanding of mammary carcinogenesis. In this study, we conducted a genome-wide association study (GWAS) as well as whole-genome imputation with 2,642 cases and 2,099 unaffected female controls. We further examined 13 suggestive loci (P<1.0×10⁻⁵) using an independent sample set of 2,885 cases and 3,395 controls and successfully validated two previously-reported loci, rs2981578 (combined P-value of 1.31×10⁻¹², OR = 1.23; 95% CI = 1.16–.30) on chromosome 10q26 (FGFR2), rs3803662 (combined P-value of 2.79×10⁻¹¹, OR = 1.21; 95% CI = 1.15–.28) and rs12922061 (combined P-value of 3.97×10⁻¹⁰, OR = 1.23; 95% CI = 1.15–.31) on chromosome 16q12 (TOX3-LOC643714). Weighted genetic risk score on the basis of three significantly associated variants and two previously reported breast cancer associated loci in East Asian population revealed that individuals who carry the most risk alleles in category 5 have 2.2 times higher risk of developing breast cancer in the Japanese population than those who carry the least risk alleles in reference category 1. Although we could not identify additional loci associated with breast cancer, our study utilized one of the largest sample sizes reported to date, and provided genetic status that represent the Japanese population. Further local and international collaborative study is essential to identify additional genetic variants that could lead to a better, accurate prediction for breast cancer.     

A digenic cause of cleft lip in A-strain mice and definition of candidate genes for the two loci

BACKGROUND: Nonsyndromic cleft lip with or without cleft palate, CL(P), is a common human birth defect with a complex unknown genetic cause. The mouse model is the "A/-" strains. Our previous studies mapped two loci: clf1 on Chr11 and clf2 on Chr13--with a strong genetic maternal effect on the level of risk. Here we test the hypothesis that CL(P) is digenic and identify candidate genes for clf1 and clf2. METHODS: We observed E14 CL(P) frequencies in backcross (BC1) embryos from a new cross of A/WySn to AXB-4/Pgn and from test crosses of three new "congenic RI" lines. Using new polymorphic markers from genes and our mapping panels of segregants and RI strains, we identified the candidate genes for clf1 and clf2. We sequenced the coding region of Ptch in A/WySn cDNA. RESULTS: Seventy new BC1 CL(P) segregants (4%) were obtained, as predicted. All three new congenic RI lines homozygous for both clf1 and clf2 had A/WySn-level CL(P) frequencies (10-30%) in test crosses. The clf1 region contains 10 known genes (Arf2, Cdc27, Crhr1, Gosr2, Itgb3, Mapt, Myl4, Nsf, Wnt3, and Wnt9b). The clf2 region contains 17 known genes with human orthologs. Both regions contain additional potential genes. No causal mutation in Ptch coding sequence was found. CONCLUSIONS: In A-strain mice, nonsyndromic CL(P) is digenic, suggesting that nonsyndromic human CL(P) may also be digenic. The orthologous human genes are on 17q (clf1) and 9q, 8q and 5p (clf2), and good candidate genes are WNT3 or WNT9B (17q), and PTCH (9q) or MTRR (5p).     

Combined analysis of genome-wide association studies for Crohn disease and psoriasis identifies seven shared susceptibility loci

Psoriasis (PS) and Crohn disease (CD) have been shown to be epidemiologically, pathologically, and therapeutically connected, but little is known about their shared genetic causes. We performed meta-analyses of five published genome-wide association studies on PS (2,529 cases and 4,955 controls) and CD (2,142 cases and 5,505 controls), followed up 20 loci that showed strongest evidence for shared disease association and, furthermore, tested cross-disease associations for previously reported PS and CD risk alleles in additional 6,115 PS cases, 4,073 CD cases, and 10,100 controls. We identified seven susceptibility loci outside the human leukocyte antigen region (9p24 near JAK2, 10q22 at ZMIZ1, 11q13 near PRDX5, 16p13 near SOCS1, 17q21 at STAT3, 19p13 near FUT2, and 22q11 at YDJC) shared between PS and CD with genome-wide significance (p < 5 × 10⁻⁸) and confirmed four already established PS and CD risk loci (IL23R, IL12B, REL, and TYK2). Three of the shared loci are also genome-wide significantly associated with PS alone (10q22 at ZMIZ1, p_rs1250544 = 3.53 × 10⁻⁸, 11q13 near PRDX5, p_rs694739 = 3.71 × 10⁻⁰⁹, 22q11 at YDJC, p_rs181359 = 8.02 × 10⁻¹⁰). In addition, we identified one susceptibility locus for CD (16p13 near SOCS1, p_rs4780355 = 4.99 × 10⁻⁸). Refinement of association signals identified shared genome-wide significant associations for exonic SNPs at 10q22 (ZMIZ1) and in silico expression quantitative trait locus analyses revealed that the associations at ZMIZ1 and near SOCS1 have a potential functional effect on gene expression. Our results show the usefulness of joint analyses of clinically distinct immune-mediated diseases and enlarge the map of shared genetic risk loci.