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.     

Maternal Transmission Effect of a PDGF-C SNP on Nonsyndromic Cleft Lip with or without Palate from a Chinese Population

Cleft lip with or without palate (CL/P) is a common congenital anomaly with a high birth prevalence in China. Based on a previous linkage signal of nonsyndromic CL/P (NSCL/P) on the chromosomal region 4q31–q32 from the Chinese populations, we screened the 4q31–q32 region for susceptibility genes in 214 trios of Han Chinese. PDGF-C, an important developmental factor, resides in the region and has been implicated in NSCL/P. However, in our family-based association test (transmission disequilibrium test; TDT), we could not conclude an association between PDGF-C and NSCL/P as previously suggested. Instead, we found strong evidence for parent-of-origin effect at a PDGF-C SNP, rs17035464, by a likelihood ratio test (unadjusted p-value = 0.0018; I_m = 2.46). The location of rs17035464 is 13 kb downstream of a previously reported, NSCL/P-associated SNP, rs28999109. Furthermore, a patient from our sample trios was observed with a maternal segmental uniparental isodisomy (UPD) in a region containing rs17035464. Our findings support the involvement of PDGF-C in the development of oral clefts; moreover, the UPD case report contributes to the collective knowledge of rare variants in the human genome.     

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.     

Joint Testing of Genotypic and Gene-Environment Interaction Identified Novel Association for BMP4 with Non-Syndromic CL/P in an Asian Population Using Data from an International Cleft Consortium

Background

Non-syndromic cleft lip with or without cleft palate (NSCL/P) is a common disorder with complex etiology. The Bone Morphogenetic Protein 4 gene (BMP4) has been considered a prime candidate gene with evidence accumulated from animal experimental studies, human linkage studies, as well as candidate gene association studies. The aim of the current study is to test for linkage and association between BMP4 and NSCL/P that could be missed in genome-wide association studies (GWAS) when genotypic (G) main effects alone were considered.

Methodology/Principal Findings

We performed the analysis considering G and interactions with multiple maternal environmental exposures using additive conditional logistic regression models in 895 Asian and 681 European complete NSCL/P trios. Single nucleotide polymorphisms (SNPs) that passed the quality control criteria among 122 genotyped and 25 imputed single nucleotide variants in and around the gene were used in analysis. Selected maternal environmental exposures during 3 months prior to and through the first trimester of pregnancy included any personal tobacco smoking, any environmental tobacco smoke in home, work place or any nearby places, any alcohol consumption and any use of multivitamin supplements. A novel significant association held for rs7156227 among Asian NSCL/P and non-syndromic cleft lip and palate (NSCLP) trios after Bonferroni correction which was not seen when G main effects alone were considered in either allelic or genotypic transmission disequilibrium tests. Odds ratios for carrying one copy of the minor allele without maternal exposure to any of the four environmental exposures were 0.58 (95%CI = 0.44, 0.75) and 0.54 (95%CI = 0.40, 0.73) for Asian NSCL/P and NSCLP trios, respectively. The Bonferroni P values corrected for the total number of 117 tested SNPs were 0.0051 (asymptotic P = 4.39*10⁻⁵) and 0.0065 (asymptotic P = 5.54*10⁻⁵), accordingly. In European trios, no significant association was seen for any SNPs after Bonferroni corrections for the total number of 120 tested SNPs.

Conclusions/Significance

Our findings add evidence from GWAS to support the role of BMP4 in susceptibility to NSCL/P originally identified in linkage and candidate gene association studies.     

Population-Specific Haplotype Association of the Postsynaptic Density Gene DLG4 with Schizophrenia,in Family-Based Association Studies

The post-synaptic density (PSD) of glutamatergic synapses harbors a multitude of proteins critical for maintaining synaptic dynamics. Alteration of protein expression levels in this matrix is a marked phenomenon of neuropsychiatric disorders including schizophrenia, where cognitive functions are impaired. To investigate the genetic relationship of genes expressed in the PSD with schizophrenia, a family-based association analysis of genetic variants in PSD genes such as DLG4, DLG1, PICK1 and MDM2, was performed, using Japanese samples (124 pedigrees, n = 376 subjects). Results showed a significant association of the rs17203281 variant from the DLG4 gene, with preferential transmission of the C allele (p = 0.02), although significance disappeared after correction for multiple testing. Replication analysis of this variant, found no association in a Chinese schizophrenia cohort (293 pedigrees, n = 1163 subjects) or in a Japanese case-control sample (n = 4182 subjects). The DLG4 expression levels between postmortem brain samples from schizophrenia patients showed no significant changes from controls. Interestingly, a five marker haplotype in DLG4, involving rs2242449, rs17203281, rs390200, rs222853 and rs222837, was enriched in a population specific manner, where the sequences A-C-C-C-A and G-C-C-C-A accumulated in Japanese (p = 0.0009) and Chinese (p = 0.0007) schizophrenia pedigree samples, respectively. However, this could not be replicated in case-control samples. None of the variants in other examined candidate genes showed any significant association in these samples. The current study highlights a putative role for DLG4 in schizophrenia pathogenesis, evidenced by haplotype association, and warrants further dense screening for variants within these haplotypes.     

Association of the IL-10 Gene Family Locus on Chromosome 1 with Juvenile Idiopathic Arthritis (JIA)

Background

The cytokine IL-10 and its family members have been implicated in autoimmune diseases and we have previously reported that genetic variants in IL-10 were associated with a rare group of diseases called juvenile idiopathic arthritis (JIA). The aim of this study was to fine map genetic variants within the IL-10 cytokine family cluster on chromosome 1 using linkage disequilibrium (LD)-tagging single nucleotide polymorphisms (tSNPs) approach with imputation and conditional analysis to test for disease associations.

Methodology/Principal Findings

Fifty-three tSNPs were tested for association between Caucasian paediatric cohorts [219 systemic JIA (sJIA), 187 persistent oligoarticular JIA (pOJIA), and 139 extended OJIA (eOJIA) patients], and controls (Wellcome Trust control cohort, WTCCC2). Significant association with sJIA was detected at rs1400986 in the promoter of IL-20 (odds ratio 1.53; 95% CI 1.21–1.93; p = 0.0004), but in no other subtypes. Imputation analysis identified additional associated SNPs for pOJIA at IL-20 and IL-24, including a rare, functional, missense variant at IL-24 with a p = 0.0002. Penalised logistic regression analysis with HyperLasso and conditional analysis identified several further associations with JIA subtypes. In particular, haplotype analysis refined the sJIA association, with a joint effect at rs1400986 and rs4129024 in intron 1 of MAPKAPK2 (p = 3.2E−5). For pOJIA, a 3-SNP haplotype including rs1878672 in intron 3 of IL-10 showed evidence for association (p = 0.0018). In eOJIA, rs10863962 (3′UTR of FCAMR) and rs12409577 (intron of IL-19) haplotype showed some evidence of association (p = 0.0003).

Conclusions

This study supports previous association of IL-20 with sJIA. Haplotype analyses provided stronger association signals than single point analyses, while a penalised logistic regression approach also suggested multiple independent association signals. Replication studies are required to confirm or refute these findings. The results indicate that combined effects with unknown/rare variants remain to be characterised in JIA, and represent a possible example of synthetic association in this region.     

Single Nucleotide Polymorphisms within Interferon Signaling Pathway Genes Are Associated with Colorectal Cancer Susceptibility and Survival

Interferon (IFN) signaling has been suggested to play an important role in colorectal carcinogenesis. Our study aimed to examine potentially functional genetic variants in interferon regulatory factor 3 (IRF3), IRF5, IRF7, type I and type II IFN and their receptor genes with respect to colorectal cancer (CRC) risk and clinical outcome. Altogether 74 single nucleotide polymorphisms (SNPs) were covered by the 34 SNPs genotyped in a hospital-based case-control study of 1327 CRC cases and 758 healthy controls from the Czech Republic. We also analyzed these SNPs in relation to overall survival and event-free survival in a subgroup of 483 patients. Seven SNPs in IFNA1, IFNA13, IFNA21, IFNK, IFNAR1 and IFNGR1 were associated with CRC risk. After multiple testing correction, the associations with the SNPs rs2856968 (IFNAR1) and rs2234711 (IFNGR1) remained formally significant (P = 0.0015 and P<0.0001, respectively). Multivariable survival analyses showed that the SNP rs6475526 (IFNA7/IFNA14) was associated with overall survival of the patients (P = 0.041 and event-free survival among patients without distant metastasis at the time of diagnosis, P = 0.034). The hazard ratios (HRs) for rs6475526 remained statistically significant even after adjustment for age, gender, grade and stage (P = 0.029 and P = 0.036, respectively), suggesting that rs6475526 is an independent prognostic marker for CRC. Our data suggest that genetic variation in the IFN signaling pathway genes may play a role in the etiology and survival of CRC and further studies are warranted.     

干扰素调节因子-6基因与非综合征性唇腭裂的关系

先天性唇腭裂常分为综合征性唇腭裂和非综合征性唇腭裂两大类,其中非综合征性唇腭裂（nonsyndromic cleft lip with orwithout cleft palate,NSCL/P）约占先天性唇腭裂的70%-80%。国内外学者在对NSCL/P相关基因进行研究后发现,干扰素调节因子6（Interferon Regulatory Factor 6,IRF6）是迄今发现最有价值的并且与NSCL/P致病有相关性的热点基因之一,但是仍有部分学者通过实验研究后得出了相反的结论,故IRF6基因与NSCL/P之间的相关性说法不一,存在较大的争论,究竟前者是通过何种遗传方式作用于后者、仍然不十分清楚,且需要大样本的研究来证实。本文就IRF6基因与NSCL/P的关系做一综述,为研究两者的关系提供系统性参考。     

Transforming Growth Factor-α (TGFA): Genomic Structure, Boundary Sequences, and Mutation Analysis in Nonsyndromic Cleft Lip/Palate and Cleft Palate Only

Transforming growth factor-α (TGFA) has been proposed as a candidate gene in the etiology of nonsyndromic cleft lip with or without cleft palate (NS-CL/P) and of nonsyndromic cleft palate only (NS-CPO). Biologic support for a role of TGFA arises from its presence at high levels in the epithelial tissue of the medial edge of the palatal shelves at the time of shelf fusion in mice. Genetic support for the role of TGFA in clefting comes from the reported association of TGFA alleles with human NS-CPO and NS-CL/P. In this study we report the sequence and structure of human genomic TGFA and the search for causal TGFA mutations in 250 individuals with NS-CL/P or NS-CPO by conformational analysis of the coding sequence, splice junctions, and a portion of the 3′ untranslated region strongly homologous between human and mouse. We confirm that human TGFA is composed of six exons and here report several new sequence substitutions and their frequencies. Five variants in conserved segments may represent rare causes for clefting in humans and provide support for the role of TGFA in facial morphogenesis.     

Genetic Variants Associated with Serum Thyroid Stimulating Hormone (TSH) Levels in European Americans and African Americans from the eMERGE Network

Thyroid stimulating hormone (TSH) hormone levels are normally tightly regulated within an individual; thus, relatively small variations may indicate thyroid disease. Genome-wide association studies (GWAS) have identified variants in PDE8B and FOXE1 that are associated with TSH levels. However, prior studies lacked racial/ethnic diversity, limiting the generalization of these findings to individuals of non-European ethnicities. The Electronic Medical Records and Genomics (eMERGE) Network is a collaboration across institutions with biobanks linked to electronic medical records (EMRs). The eMERGE Network uses EMR-derived phenotypes to perform GWAS in diverse populations for a variety of phenotypes. In this report, we identified serum TSH levels from 4,501 European American and 351 African American euthyroid individuals in the eMERGE Network with existing GWAS data. Tests of association were performed using linear regression and adjusted for age, sex, body mass index (BMI), and principal components, assuming an additive genetic model. Our results replicate the known association of PDE8B with serum TSH levels in European Americans (rs2046045 p = 1.85×10⁻¹⁷, β = 0.09). FOXE1 variants, associated with hypothyroidism, were not genome-wide significant (rs10759944: p = 1.08×10⁻⁶, β = −0.05). No SNPs reached genome-wide significance in African Americans. However, multiple known associations with TSH levels in European ancestry were nominally significant in African Americans, including PDE8B (rs2046045 p = 0.03, β = −0.09), VEGFA (rs11755845 p = 0.01, β = −0.13), and NFIA (rs334699 p = 1.50×10⁻³, β = −0.17). We found little evidence that SNPs previously associated with other thyroid-related disorders were associated with serum TSH levels in this study. These results support the previously reported association between PDE8B and serum TSH levels in European Americans and emphasize the need for additional genetic studies in more diverse populations.     

Identification of novel susceptibility loci for non‐syndromic cleft lip with or without cleft palate

Although several genome‐wide association studies (GWAS) of non‐syndromic cleft lip with or without cleft palate (NSCL/P) have been reported, more novel association signals are remained to be exploited. Here, we performed an in‐depth analysis of our previously published Chinese GWAS cohort study with replication in an extra dbGaP case‐parent trios and another in‐house Nanjing cohort, and finally identified five novel significant association signals (rs11119445: 3’ of SERTAD4, P = 6.44 × 10⁻¹⁴; rs227227 and rs12561877: intron of SYT14, P = 5.02 × 10⁻¹³ and 2.80 × 10⁻¹¹, respectively; rs643118: intron of TRAF3IP3, P = 4.45 × 10⁻⁶; rs2095293: intron of NR6A1, P = 2.98 × 10⁻⁵). The mean (standard deviation) of the weighted genetic risk score (wGRS) from these SNPs was 1.83 (0.65) for NSCL/P cases and 1.58 (0.68) for controls, respectively (P = 2.67 × 10⁻¹⁶). Rs643118 was identified as a shared susceptible factor of NSCL/P among Asians and Europeans, while rs227227 may contribute to the risk of NSCL/P as well as NSCPO. In addition, sertad4 knockdown zebrafish models resulted in down‐regulation of sox2 and caused oedema around the heart and mandibular deficiency, compared with control embryos. Taken together, this study has improved our understanding of the genetic susceptibility to NSCL/P and provided further clues to its aetiology in the Chinese population.     

Genome-Wide Association Study Identifies GPC5 as a Novel Genetic Locus Protective against Sudden Cardiac Arrest

Background

Existing studies indicate a significant genetic component for sudden cardiac arrest (SCA) and genome-wide association studies (GWAS) provide an unbiased approach for identification of novel genes. We performed a GWAS to identify genetic determinants of SCA.

Methodology/Principal Findings

We used a case-control design within the ongoing Oregon Sudden Unexpected Death Study (Oregon-SUDS). Cases (n = 424) were SCAs with coronary artery disease (CAD) among residents of Portland, OR (2002–07, population ∼1,000,000) and controls (n = 226) were residents with CAD, but no history of SCA. All subjects were of White-European ancestry and GWAS was performed using Affymetrix 500K/5.0 and 6.0 arrays. High signal markers were genotyped in SCA cases (n = 521) identified from the Atherosclerosis Risk in Communities Study (ARIC) and the Cardiovascular Health Study (CHS) (combined n = 19,611). No SNPs reached genome-wide significance (p<5×10⁻⁸). SNPs at 6 loci were prioritized for follow-up primarily based on significance of p<10⁻⁴ and proximity to a known gene (CSMD2, GPR37L1, LIN9, B4GALNT3, GPC5, and ZNF592). The minor allele of GPC5 (GLYPICAN 5, rs3864180) was associated with a lower risk of SCA in Oregon-SUDS, an effect that was also observed in ARIC/CHS whites (p<0.05) and blacks (p<0.04). In a combined Cox proportional hazards model analysis that adjusted for race, the minor allele exhibited a hazard ratio of 0.85 (95% CI 0.74 to 0.98; p<0.01).

Conclusions/Significance

A novel genetic locus for SCA, GPC5, was identified from Oregon-SUDS and successfully validated in the ARIC and CHS cohorts. Three other members of the Glypican family have been previously implicated in human disease, including cardiac conditions. The mechanism of this specific association requires further study.     

Genetic Association Study of TNFAIP3, IFIH1, IRF5 Polymorphisms with Polymyositis/Dermatomyositis in Chinese Han Population

Background

Single-nucleotide polymorphisms (SNPs) in the TNFAIP3, IFIH1, and IRF5 genes have been associated with several auto-inflammation diseases, while the susceptibility between these genes and idiopathic inflammatory myopathies (IIMs) were not reported. This study aimed to investigate whether TNFAIP3, IFIH1, and IRF5 gene polymorphisms confer susceptibility for the IIMs in Chinese Han population.

Methods

A large case–control study of Chinese subjects with polymyositis (PM) (n = 298) and dermatomyositis (DM) (n = 530) was accomplished. 968 healthy and ethnically matched controls were available for comparison. Six SNPs in the TNFAIP3 region (rs2230926 and rs5029939), the IFIH1 gene (rs1990760 and rs3747517) and the IRF5 region (rs4728142 and rs729302) were assessed and genotyped using the Sequenom MassArray iPLEX platform.

Results

Our study indicated a strong allele association was observed in PM/DM and PM patients for rs2230926 (OR: 1.61, 95%CI: 1.20–2.16, P_c = 7.5×10⁻³; OR: 1.88, 95%CI: 1.30–2.74, P_c = 4.0×10⁻³, respectively) and rs5029939 (OR: 1.64, 95%CI: 1.21–2.21, P_c = 6.0×10⁻³; OR: 1.88, 95%CI: 1.28–2.76, P_c = 5.5×10⁻³,respectively). And rs2230926 and rs5029939 were significantly associated with interstitial lung disease (ILD) in PM/DM and PM patients (P_c = 0.04 and P_c = 0.016; P_c = 0.02 and P_c = 0.03, respectively). In addition, rs4728142 allele and genotype had significant association with PM/DM patients (P_c = 0.026 and P_c = 0.048, respectively). Further analysis with three logistic regression genetic models revealed statistically significant difference in the genotypic distribution in the PM/DM, PM or DM patients when the additive and dominant models were used.

Conclusions

This was the first study to reveal TNFAIP3 and IRF5 polymorphisms were associated with PM/DM patients or these patients with ILD, indicating that TNFAIP3 and IRF5 might be the susceptibility gene for PM/DM patients in Chinese Han population.     

Genome Wide Meta-analysis Highlights the Role of Genetic Variation in RARRES2 in the Regulation of Circulating Serum Chemerin

Chemerin is an adipokine proposed to link obesity and chronic inflammation of adipose tissue. Genetic factors determining chemerin release from adipose tissue are yet unknown. We conducted a meta-analysis of genome-wide association studies (GWAS) for serum chemerin in three independent cohorts from Europe: Sorbs and KORA from Germany and PPP-Botnia from Finland (total N = 2,791). In addition, we measured mRNA expression of genes within the associated loci in peripheral mononuclear cells by micro-arrays, and within adipose tissue by quantitative RT-PCR and performed mRNA expression quantitative trait and expression-chemerin association studies to functionally substantiate our loci. Heritability estimate of circulating chemerin levels was 16.2% in the Sorbs cohort. Thirty single nucleotide polymorphisms (SNPs) at chromosome 7 within the retinoic acid receptor responder 2 (RARRES2)/Leucine Rich Repeat Containing (LRRC61) locus reached genome-wide significance (p<5.0×10⁻⁸) in the meta-analysis (the strongest evidence for association at rs7806429 with p = 7.8×10⁻¹⁴, beta = −0.067, explained variance 2.0%). All other SNPs within the cluster were in linkage disequilibrium with rs7806429 (minimum r² = 0.43 in the Sorbs cohort). The results of the subgroup analyses of males and females were consistent with the results found in the total cohort. No significant SNP-sex interaction was observed. rs7806429 was associated with mRNA expression of RARRES2 in visceral adipose tissue in women (p<0.05 after adjusting for age and body mass index). In conclusion, the present meta-GWAS combined with mRNA expression studies highlights the role of genetic variation in the RARRES2 locus in the regulation of circulating chemerin concentrations.     

Machine Learning Models for Genetic Risk Assessment of Infants with Non-syndromic Orofacial Cleft

The isolated type of orofacial cleft, termed non-syndromic cleft lip with or without cleft palate (NSCL/P), is the second most common birth defect in China, with Asians having the highest incidence in the world. NSCL/P involves multiple genes and complex interactions between genetic and environmental factors, imposing difficulty for the genetic assessment of the unborn fetus carrying multiple NSCL/P-susceptible variants. Although genome-wide association studies (GWAS) have uncovered dozens of single nucleotide polymorphism (SNP) loci in different ethnic populations, the genetic diagnostic effectiveness of these SNPs requires further experimental validation in Chinese populations before a diagnostic panel or a predictive model covering multiple SNPs can be built. In this study, we collected blood samples from control and NSCL/P infants in Han and Uyghur Chinese populations to validate the diagnostic effectiveness of 43 candidate SNPs previously detected using GWAS. We then built predictive models with the validated SNPs using different machine learning algorithms and evaluated their prediction performance. Our results showed that logistic regression had the best performance for risk assessment according to the area under curve. Notably, defective variants in MTHFR and RBP4, two genes involved in folic acid and vitamin A biosynthesis, were found to have high contributions to NSCL/P incidence based on feature importance evaluation with logistic regression. This is consistent with the notion that folic acid and vitamin A are both essential nutritional supplements for pregnant women to reduce the risk of conceiving an NSCL/P baby. Moreover, we observed a lower predictive power in Uyghur than in Han cases, likely due to differences in genetic background between these two ethnic populations. Thus, our study highlights the urgency to generate the HapMap for Uyghur population and perform resequencing-based screening of Uyghur-specific NSCL/P markers.     

Host Genetic Factors and Vaccine-Induced Immunity to HBV Infection: Haplotype Analysis

Hepatitis B virus (HBV) infection remains a significant health burden world-wide, although vaccines help decrease this problem. We previously identified associations of single nucleotide polymorphisms in several candidate genes with vaccine-induced peak antibody level (anti-HBs), which is predictive of long-term vaccine efficacy and protection against infection and persistent carriage; here we report on a haplotype-based analysis. A total of 688 SNPs from 117 genes were examined for a two, three and four sliding window haplotype analysis in a Gambian cohort. Analysis was performed on 197 unrelated individuals, 454 individuals from 174 families, and the combined sample (N = 651). Global and individual haplotype association tests were carried out (adjusted for covariates), employing peak anti-HBs level as outcome. Five genes (CD44, CD58, CDC42, IL19 and IL1R1) had at least one significant haplotype in the unrelated or family analysis as well as the combined analysis. Previous single locus results were confirmed for CD44 (combined global p = 9.1×10⁻⁵ for rs353644-rs353630-rs7937602) and CD58 (combined global p = 0.008 for rs1414275-rs11588376-rs1016140). Haplotypes in CDC42, IL19 and IL1R1 also associated with peak anti-HBs level. We have identified strong haplotype effects on HBV vaccine-induced antibody level in five genes, three of which, CDC42, IL19 and IL1R1, did not show evidence of association in a single SNP analyses and corroborated the majority of these effects in two datasets. The haplotype analysis identified associations with HBV vaccine-induced immunity in several new genes.     

A Meta-Analysis of P2X7 Gene-762T/C Polymorphism and Pulmonary Tuberculosis Susceptibility

AimWe performed a comprehensive meta-analysis to determine the association between P2X7 -762T/C polymorphism and pulmonary tuberculosis susceptibility.MethodologyBased on comprehensive searches of the PubMed, SCI, Elsevier, China National Knowledge Infrastructure (CNKI) and Wanfang Database, we identified eligible studies about the association between P2X7 -762T/C polymorphism and pulmonary tuberculosis risk. Pooled odds ratio (ORs) and 95% confidence intervals (95%CIs) were calculated in random-effects model.ResultsA total of 2207 tuberculosis cases and 2220 controls in 8 case-control studies were included in this meta-analysis. Allele model (C vs. T: p = 0.15; OR = 0.83, 95% CI = 0.65–1.07), homozygous model (CC vs. TT: p = 0.23; OR = 0.73, 95% CI = 0.44 to 1.22), and heterozygous model (CT vs. TT: p = 0.57; OR = 0.92, 95% CI = 0.68 to 1.24) did not show increased risk of developing pulmonary tuberculosis. Similarly, dominant model (CC+CT vs. TT: p = 0.32; OR = 0.84, 95% CI = 0.59 to 1.19) and recessive model (CC vs. CT+TT: p = 0.08; OR = 0.77, 95% CI = 0.57 to 1.04) failed to show increased risk of developing pulmonary tuberculosis. Subgroup analysis by ethnicity did not detect any significant association between P2X7–762T/C polymorphism and pulmonary tuberculosis susceptibility.ConclusionsP2X7 -762T/C gene polymorphism is not associated with pulmonary tuberculosis susceptibility.     

A Novel MMP12 Locus Is Associated with Large Artery Atherosclerotic Stroke Using a Genome-Wide Age-at-Onset Informed Approach

Genome-wide association studies (GWAS) have begun to identify the common genetic component to ischaemic stroke (IS). However, IS has considerable phenotypic heterogeneity. Where clinical covariates explain a large fraction of disease risk, covariate informed designs can increase power to detect associations. As prevalence rates in IS are markedly affected by age, and younger onset cases may have higher genetic predisposition, we investigated whether an age-at-onset informed approach could detect novel associations with IS and its subtypes; cardioembolic (CE), large artery atherosclerosis (LAA) and small vessel disease (SVD) in 6,778 cases of European ancestry and 12,095 ancestry-matched controls. Regression analysis to identify SNP associations was performed on posterior liabilities after conditioning on age-at-onset and affection status. We sought further evidence of an association with LAA in 1,881 cases and 50,817 controls, and examined mRNA expression levels of the nearby genes in atherosclerotic carotid artery plaques. Secondly, we performed permutation analyses to evaluate the extent to which age-at-onset informed analysis improves significance for novel loci. We identified a novel association with an MMP12 locus in LAA (rs660599; p = 2.5×10⁻⁷), with independent replication in a second population (p = 0.0048, OR(95% CI) = 1.18(1.05–1.32); meta-analysis p = 2.6×10⁻⁸). The nearby gene, MMP12, was significantly overexpressed in carotid plaques compared to atherosclerosis-free control arteries (p = 1.2×10⁻¹⁵; fold change = 335.6). Permutation analyses demonstrated improved significance for associations when accounting for age-at-onset in all four stroke phenotypes (p<0.001). Our results show that a covariate-informed design, by adjusting for age-at-onset of stroke, can detect variants not identified by conventional GWAS.     

Genetic Variation in Circadian Rhythm Genes CLOCK and ARNTL as Risk Factor for Male Infertility

Background

The circadian system has a major role in maintaining homeostasis and proper body functions including reproductive capacity. The aim of this study was to examine whether there is an association between genetic variability in the primary clock genes CLOCK and ARNTL and male infertility in humans.

Methodology/Principal Findings

We performed a case-control study, where we searched for an association between polymorphisms of CLOCK and ARNTL genes and male infertility in 961 Slovenian and Serbian Caucasian men. The study group consisted of 517 patients with idiopathic infertility and a control group of 444 fertile men. A statistically significant difference was found in genotype distribution between the two groups in the CLOCK gene: rs11932595 (p = 6·10⁻⁵, q = 4·10⁻⁴, OR equaled 1.9 with 95% CI 1.4–2.7), rs6811520 (p = 2·10⁻³, q = 8·10⁻³, OR = 1.7 with 95% CI 1.2–2.2) and rs6850524 (p = 0.01, q = 0.02, OR = 1.4 with 95% CI 1.1–1.9). Further analyses of haplotypes were consistent with genotyping results.

Conclusions/Significance

We provide evidence that genetic variability in the CLOCK gene might be associated with male infertility warranting further confirmation and mechanistic investigations.