Replication of association of IL1 gene complex members with ankylosing spondylitis in Taiwanese Chinese

Objective: To test the association of interleukin 1 ( IL1) gene family members with ankylosing spondylitis ( AS), previously reported in Europid subjects, in an ethnically remote population. Methods: 200 Taiwanese Chinese AS patients and 200 ethnically matched healthy controls were genotyped for five single nucleotide polymorphisms ( SNPs) and the IL1RN. VNTR, markers previously associated with AS. Allele, genotype, and haplotype frequencies were compared between cases and controls. Results: Association of alleles and genotypes of the markers IL1F10.3, IL1RN. 4, and IL1RN. VNTR was observed with AS ( p < 0.05). Haplotypes of pairs of these markers and of the markers IL1RN.6/1 and IL1RN.6/2 were also significantly associated with AS. The strongest associations observed were with the marker IL1RN. 4, and with the two-marker haplotype IL1RN.4-IL1RN.VNTR ( both p = 0.004). Strong linkage disequilibrium was observed between all marker pairs except those involving IL1B-511 ( D' 0.4 to 0.9, p < 0.01). Conclusions: The IL1 gene cluster is associated with AS in Taiwanese Chinese. This finding provides strong statistical support that the previously observed association of this gene cluster with AS is a true positive finding.

Identification of the semaphorin receptor PLXNA2 as a candidate for susceptibility to schizophrenia

The discovery of genetic factors that contribute to schizophrenia susceptibility is a key challenge in understanding the etiology of this disease. Here, we report the identification of a novel schizophrenia candidate gene on chromosome 1q32, plexin A2 (PLXNA2), in a genome-wide association study using 320 patients with schizophrenia of European descent and 325 matched controls. Over 25 000 single-nucleotide polymorphisms (SNPs) located within approximately 14 000 genes were tested. Out of 62 markers found to be associated with disease status, the most consistent finding was observed for a candidate locus on chromosome 1q32. The marker SNP rs752016 showed suggestive association with schizophrenia (odds ratio (OR) = 1.49, P = 0.006). This result was confirmed in an independent case control sample of European Americans (combined OR = 1.38, P = 0.035) and similar genetic effects were observed in smaller subsets of Latin Americans (OR = 1.26) and Asian Americans (OR = 1.37). Supporting evidence was also obtained from two family-based collections, one of which reached statistical significance (OR = 2.2, P = 0.02). High-density SNP mapping showed that the region of association spans approximately 60 kb of the PLXNA2 gene. Eight out of 14 SNPs genotyped showed statistically significant differences between cases and controls. These results are in accordance with previous genetic findings that identified chromosome 1q32 as a candidate region for schizophrenia. PLXNA2 is a member of the transmembrane semaphorin receptor family that is involved in axonal guidance during development and may modulate neuronal plasticity and regeneration. The PLXNA2 ligand semaphorin 3A has been shown to be upregulated in the cerebellum of individuals with schizophrenia. These observations, together with the genetic results, make PLXNA2 a likely candidate for the 1q32 schizophrenia susceptibility locus.

Teaching principles of genetics through SNP analysis

An understanding of inheritance requires comprehension of genetic processes at all levels, from molecules to populations. Frequently genetics courses are separated into molecular and organismal genetics and students may fail to see the relationships between them. This is particularly true with human genetics, because of the difficulties in designing experimental approaches which are consistent with ethical restrictions, student abilities and background knowledge, and available time and materials. During 2005 we used analysis of single nucleotide polymorphisms (SNPs) in two genetic regions to enhance student learning and provide a practical experience in human genetics. Students scanned databases to discover SNPs in a gene of interest, used software to design PCR primers and a restriction enzyme based assay for the alleles, and carried out an analysis of the SNP on anonymous individual and family DNAs. The project occupied eight to ten hours per week for one semester, with some time spent in the laboratory and some spent in database searching, reading and writing the report. In completing their projects, students acquired a knowledge of Mendel’s first law (through looking at inheritance patterns), Mendel’s second law and the exceptions (the concepts of linkage and linkage disequilibrium), DNA structure (primer design and restriction enzyme analysis) and function (SNPs in coding and non-coding regions), population genetics and the statistical analysis of allele frequencies, genomics, bioinformatics and the ethical issues associated with the use of human samples. They also developed skills in presentation of results by publication and conference participation. Deficiencies in their understanding (for example of inheritance patterns, gene structure, statistical approaches and report writing) were detected and guidance given during the project. SNP analysis was found to be a powerful approach to enhance and integrate student understanding of genetic concepts.

Circadian gene variants and susceptibility to type 2 diabetes:a pilot study

Disruption of endogenous circadian rhythms has been shown to increase the risk of developing type 2 diabetes, suggesting that circadian genes might play a role in determining disease susceptibility. We present the results of a pilot study investigating the association between type 2 diabetes and selected single nucleotide polymorphisms (SNPs) in/near nine circadian genes. The variants were chosen based on their previously reported association with prostate cancer, a disease that has been suggested to have a genetic link with type 2 diabetes through a number of shared inherited risk determinants.

Common variants of the TCF7L2 gene are associated with increased risk of type 2 diabetes mellitus in a UK-resident South Asian population

Background - Recent studies have implicated variants of the transcription factor 7-like 2 (TCF7L2) gene in genetic susceptibility to type 2 diabetes mellitus in several different populations. The aim of this study was to determine whether variants of this gene are also risk factors for type 2 diabetes development in a UK-resident South Asian cohort of Punjabi ancestry. Methods - We genotyped four single nucleotide polymorphisms (SNPs) of TCF7L2 (rs7901695, rs7903146, rs11196205 and rs12255372) in 831 subjects with diabetes and 437 control subjects. Results - The minor allele of each variant was significantly associated with type 2 diabetes; the greatest risk of developing the disease was conferred by rs7903146, with an allelic odds ratio (OR) of 1.31 (95% CI: 1.11 – 1.56, p = 1.96 × 10-3). For each variant, disease risk associated with homozygosity for the minor allele was greater than that for heterozygotes, with the exception of rs12255372. To determine the effect on the observed associations of including young control subjects in our data set, we reanalysed the data using subsets of the control group defined by different minimum age thresholds. Increasing the minimum age of our control subjects resulted in a corresponding increase in OR for all variants of the gene (p ≤ 1.04 × 10-7). Conclusion - Our results support recent findings that TCF7L2 is an important genetic risk factor for the development of type 2 diabetes in multiple ethnic groups.

Entropy Based Approach to Finding Interacting Genes Responsible for Complex Human Disease

2000 Mathematics Subject Classification: 62P10, 92D10, 92D30, 94A17, 62L10.

Overweight/Obesity and HIV Disease Progression in HIV+ Adults in Botswana

Studies indicate that overweight and obesity protect against HIV-disease progression in antiretroviral therapy (ART)-naïve patients. We examined retrospectively the relationship of overweight/obesity with HIV-disease progression in ART-naïve HIV+ adults in Botswana in a case-control study with 18-month follow-up, which included 217 participants, 139 with BMI 18.0-24.9 kg/m2 and 78 with BMI ≥25 kg/m2. Archived plasma samples were used to determine inflammatory markers: leptin and bacterial endotoxin lipopolysaccharide (LPS), and genotype single nucleotide polymorphisms (SNPs) of the Fat Mass and Obesity Associated Gene (FTO). At baseline, BMI was inversely associated with risk for AIDS-defining conditions (HR=0.218; 95%CI=0.068, 0.701, P=0.011), and higher fat mass was associated with reduced risk of the combined outcome of CD4+cell count ≤250/µL and AIDS-defining conditions, whichever occurred earlier (HR=0.918; 95%CI=0.847, 0.994, P=0.036) over 18 months, adjusting for age, gender, marriage, children, and baseline CD4+cell count and HIV-viral load. FTO-SNP rs17817449 was associated with BMI (OR=1.082; 95%CI=1.001, 1.169; P=0.047). Fat mass was associated with the risk alleles of rs1121980 (OR=1.065; 95%CI=1.009, 1.125, P=0.021), rs8050136 (OR=1.078; 95%CI=1.021, 1.140; P=0.007), and rs17817449 (OR=1.086; 95%CI=1.031, 1.145; P=0.002), controlling for age, gender, tribe, total energy intake, and activity. There were no associations of SNPs with markers of disease progression. Leptin levels were positively associated with BMI (β=1.764; 95%CI=0.788, 2.739; P=0.022) and fat mass (β=0.112; 95%CI=0.090, 0.135; P<0.001), but inversely with viral load (β=-0.305; 95%CI=-0.579, -.031; P=0.030). LPS levels were inversely associated with BMI (OR=0.790, 95%CI=0.630, 0.990; P=0.041), and fat mass (OR=0.852, 95%CI=0.757, 0.958; P=0.007) and directly with viral load (OR=2.608, 95%CI=1.111, 6.124; P=0.028), adjusting for age, gender, smoking and %fat mass. In this cohort, overweight/obesity predicted slower HIV-disease progression. Obesity may confer an advantage in maintaining fat stores to support the overactive immune system. FTO-SNPs may contribute to the variation in fat mass; however, they were not associated with HIV-disease progression. Our findings suggest that the obesity paradox may be explained by the association of increased LPS with lower BMI and higher viral load; while viral load decreased with increasing leptin levels. Studies in African populations are needed to clarify whether genetic variation and inflammation mediate the obesity paradox in HIV-disease progression.

Avaliação de polimorfismos funcionais nos genes de reparo XRCC1, APEX1, XPD e XPF em carcinomas de células escamosas orais

Base excision repair (BER) and nucleotide excision repair (NER) pathways play critical role in maintaining genome integrity. Polymorphisms in BER and NER genes which modulate the DNA repair capacity may affect the susceptibility and prognosis of oral cancer. This study was conducted with genomic DNA from 92 patients with oral squamous cell carcinomas (OSCC) and 130 controls. The cases were followed up to explore the associations between BER and NER genes polymorphisms and the risk and prognosis of OSCC. Four single-nucleotide polymorphisms (SNPs) in XRCC1 (rs25487), APEX1 (rs1130409), XPD (rs13181) and XPF (rs1799797) genes were tested by polymerase chain reaction – quantitative real time method. The GraphPad Prism version 6.0.1 statistical software was applied for statistical analysis of association. Odds ratio (OR), hazard ratio (HR), and their 95 % confidence intervals (CIs) were calculated by logistic regression. Kaplan-Meier curve and Cox proportional hazard model were used for prognostic analysis. The presence of polymorphic variants in XRCC1, APEX1, XPD and XPF genes were not associated with an increased risk of OSCC. Gene-environment interactions with smoking were not significant for any polymorphism. The presence of polymorphic variants of the XPD gene in association with alcohol consumption conferred an increased risk of 1.86 (95% CI: 0.86 – 4.01, p=0.03) for OSCC. Only APEX1 was associated with decreased specific survival (HR 3.94, 95% CI: 1.31 – 11.88, p=0.01). These results suggest an interaction between polymorphic variants of the XPF gene and alcohol consumption. Additionally APEX1 may represent a prognostic marker for OSCC.

Valor prognóstico de polimorfismos nos genes de reparo do DNA XRCC3 E RAD51 em pacientes com carcinoma epidermóide oral e de orofaringe

Faults in the genes responsible for repairs to the DNA can influence the onset of cancer or affect the response to treatment. This research evaluated the frequency of three single nucleotide polymorphisms (SNPs) in two repair genes DNA RAD51 172g> T (rs1801321), RAD51 135G> C (rs1801320) and XRCC3 T241M (rs861539) in individuals without cancer (n = 130) and patients with oral squamous cell carcinoma (OSC) and carcinoma oropharyngeal squamous (ORSC) (n = 126) and investigated possible relationships of these findings with clinical and pathological data and clinical outcomes: tumor response to radiotherapy and chemotherapy, disease-free survival, and overall survival. It was found that the allele and genotype frequencies were in equilibrium Hard-Weinberg equilibrium. The presence of at least one polymorphic allele in XRCC3 (rs861539) gene is associated with histological grade (WHO) higher (p = 0.007). We observed a higher recurrence rate trend (p = 0.08) and more advanced stage (p = 0.08) in the group that had at least one polymorphic allele of RAD51 gene (rs1801321). The presence of the analyzed SNPs not proved to be a risk factor for the development of CEO or CEOR; however, when combined with smoking or drinking, increased the risk of developing cancer from three to one hundred and fifty times. The tumor response to radiotherapy and chemotherapy was similar in patients with and without SNPs. No polymorphism showed statistical significance in relation to recurrence-free survival or overall survival. We conclude that the presence of at least one polymorphic allele of the SNPs rs861539 in XRCC3 gene, rs1801320 and rs1801321 in the RAD51 gene increase the risk of development of OSC and ORSC, when associated with the habit of drinking or smoking. Polymorphisms studied in XRCC3 and RAD51 genes are not associated with response to radiation therapy, relapse-free survival or overall survival.

Associations between smoking behavior-related alleles and the risk of melanoma.

Several studies have reported that cigarette smoking is inversely associated with the risk of melanoma. This study further tested whether incorporating genetic factors will provide another level of evaluation of mechanisms underlying the association between smoking and risk of melanoma. We investigated the association between SNPs selected from genome-wide association studies (GWAS) on smoking behaviors and risk of melanoma using 2,298 melanoma cases and 6,654 controls. Among 16 SNPs, three (rs16969968 [A], rs1051730 [A] and rs2036534 [C] in the 15q25.1 region) reached significance for association with melanoma risk in men (0.01 < = P values < = 0.02; 0.85 < = Odds Ratios (ORs) <= 1.20). There was association between the genetic scores based on the number of smoking behavior-risk alleles and melanoma risk with P-trend = 0.005 among HPFS. Further association with smoking behaviors indicating those three SNPs (rs16969968 [A], rs1051730 [A] and rs2036534 [C]) significantly associated with number of cigarettes smoked per day, CPD, with P = 0.009, 0.011 and 0.001 respectively. The SNPs rs215605 in the PDE1C gene and rs6265 in the BDNF gene significantly interacted with smoking status on melanoma risk (interaction P = 0.005 and P = 0.003 respectively). Our study suggests that smoking behavior-related SNPs are likely to play a role in melanoma development and the potential public health importance of polymorphisms in the CHRNA5-A3-B4 gene cluster. Further larger studies are warranted to validate the findings.

Risk Prediction for Epithelial Ovarian Cancer in 11 United States-Based Case-Control Studies: Incorporation of Epidemiologic Risk Factors and 17 Confirmed Genetic Loci.

Previously developed models for predicting absolute risk of invasive epithelial ovarian cancer have included a limited number of risk factors and have had low discriminatory power (area under the receiver operating characteristic curve (AUC) < 0.60). Because of this, we developed and internally validated a relative risk prediction model that incorporates 17 established epidemiologic risk factors and 17 genome-wide significant single nucleotide polymorphisms (SNPs) using data from 11 case-control studies in the United States (5,793 cases; 9,512 controls) from the Ovarian Cancer Association Consortium (data accrued from 1992 to 2010). We developed a hierarchical logistic regression model for predicting case-control status that included imputation of missing data. We randomly divided the data into an 80% training sample and used the remaining 20% for model evaluation. The AUC for the full model was 0.664. A reduced model without SNPs performed similarly (AUC = 0.649). Both models performed better than a baseline model that included age and study site only (AUC = 0.563). The best predictive power was obtained in the full model among women younger than 50 years of age (AUC = 0.714); however, the addition of SNPs increased the AUC the most for women older than 50 years of age (AUC = 0.638 vs. 0.616). Adapting this improved model to estimate absolute risk and evaluating it in prospective data sets is warranted.

Association between the oxytocin receptor (OXTR) gene and mesolimbic responses to rewards.

BACKGROUND: There has been significant progress in identifying genes that confer risk for autism spectrum disorders (ASDs). However, the heterogeneity of symptom presentation in ASDs impedes the detection of ASD risk genes. One approach to understanding genetic influences on ASD symptom expression is to evaluate relations between variants of ASD candidate genes and neural endophenotypes in unaffected samples. Allelic variations in the oxytocin receptor (OXTR) gene confer small but significant risk for ASDs for which the underlying mechanisms may involve associations between variability in oxytocin signaling pathways and neural response to rewards. The purpose of this preliminary study was to investigate the influence of allelic variability in the OXTR gene on neural responses to monetary rewards in healthy adults using functional magnetic resonance imaging (fMRI). METHODS: The moderating effects of three single nucleotide polymorphisms (SNPs) (rs1042778, rs2268493 and rs237887) of the OXTR gene on mesolimbic responses to rewards were evaluated using a monetary incentive delay fMRI task. RESULTS: T homozygotes of the rs2268493 SNP demonstrated relatively decreased activation in mesolimbic reward circuitry (including the nucleus accumbens, amygdala, insula, thalamus and prefrontal cortical regions) during the anticipation of rewards but not during the outcome phase of the task. Allelic variation of the rs1042778 and rs237887 SNPs did not moderate mesolimbic activation during either reward anticipation or outcomes. CONCLUSIONS: This preliminary study suggests that the OXTR SNP rs2268493, which has been previously identified as an ASD risk gene, moderates mesolimbic responses during reward anticipation. Given previous findings of decreased mesolimbic activation during reward anticipation in ASD, the present results suggest that OXTR may confer ASD risk via influences on the neural systems that support reward anticipation.

The differential impact of oxytocin receptor (OXTR) genotypes on the risk of autism spectrum disorder (ASD) and resulting social communication deficits

Background: Autism spectrum disorder (ASD) is multifactorial and is likely the result of complex interactions between multiple environmental and genetic factors. Recently, it has been suggested that each symptom cluster of the disorder, such as poor social communication, may be mediated by different genetic influences. Genes in the oxytocin pathway, which mediates social behaviours in humans, have been studied with single nucleotide polymorphisms (SNPs) in the oxytocin receptor gene (OXTR) being implicated in ASD. This thesis examines the presence of different oxytocin receptor genotypes, and their associations with ASD and resulting social communication deficits. Methods: The relationship between four OXTR variants and ASD was evaluated in 607 ASD simplex (SPX) families. Cases were compared to their unaffected siblings using a conditional logistic approach. Odds ratios and associated 95 percent confidence intervals were obtained. A second sample of 235 individuals with a diagnosis of ASD was examined to evaluate whether these four OXTR variants were associated with social communication scores on the Autism Diagnostic Interview – Revised (ADI-R). Parameter estimates and associated 95 percent confidence intervals were generated using a linear regression approach. Multiple testing issues were addressed using false discovery adjustments. Results: The rs53576 AG genotype was significantly associated with a lower risk of ASD (OR = 0.707, 95% CI: 0.512-0.975). A single genotype (AG) provided by the rs2254298 marker was found to be significantly associated with higher social communication scores (Parameter estimate = 1.833, SE = 0.762, p = 0.0171). This association was also seen in a Caucasian only and mothers as the respondent samples. No association was significant following false discovery rate adjustments. Conclusion: The findings from these studies provide limited support for the role of OXTR SNPs in ASD, especially in social communication skills. The clinical significance of these associations remains unknown, however, it is likely that these associations do not play a role in the severity of symptoms associated with ASD. Rather, they may be important in the appearance of social deficits due to the rs2254298 markers association with enlarged amygdalas.

Evidence for three genetic loci involved in both anorexia nervosa risk and variation of body mass index

The maintenance of normal body weight is disrupted in patients with anorexia nervosa (AN) for prolonged periods of time. Prior to the onset of AN, premorbid body mass index (BMI) spans the entire range from underweight to obese. After recovery, patients have reduced rates of overweight and obesity. As such, loci involved in body weight regulation may also be relevant for AN and vice versa. Our primary analysis comprised a cross-trait analysis of the 1000 single-nucleotide polymorphisms (SNPs) with the lowest P-values in a genome-wide association meta-analysis (GWAMA) of AN (GCAN) for evidence of association in the largest published GWAMA for BMI (GIANT). Subsequently we performed sex-stratified analyses for these 1000 SNPs. Functional ex vivo studies on four genes ensued. Lastly, a look-up of GWAMA-derived BMI-related loci was performed in the AN GWAMA. We detected significant associations (P-values <5 × 10(-5), Bonferroni-corrected P<0.05) for nine SNP alleles at three independent loci. Interestingly, all AN susceptibility alleles were consistently associated with increased BMI. None of the genes (chr. 10: CTBP2, chr. 19: CCNE1, chr. 2: CARF and NBEAL1; the latter is a region with high linkage disequilibrium) nearest to these SNPs has previously been associated with AN or obesity. Sex-stratified analyses revealed that the strongest BMI signal originated predominantly from females (chr. 10 rs1561589; Poverall: 2.47 × 10(-06)/Pfemales: 3.45 × 10(-07)/Pmales: 0.043). Functional ex vivo studies in mice revealed reduced hypothalamic expression of Ctbp2 and Nbeal1 after fasting. Hypothalamic expression of Ctbp2 was increased in diet-induced obese (DIO) mice as compared with age-matched lean controls. We observed no evidence for associations for the look-up of BMI-related loci in the AN GWAMA. A cross-trait analysis of AN and BMI loci revealed variants at three chromosomal loci with potential joint impact. The chromosome 10 locus is particularly promising given that the association with obesity was primarily driven by females. In addition, the detected altered hypothalamic expression patterns of Ctbp2 and Nbeal1 as a result of fasting and DIO implicate these genes in weight regulation.

Increasing the power of genome wide association studies in natural populations using repeated measures : evaluation and implementation

1. Genomewide association studies (GWAS) enable detailed dissections of the genetic basis for organisms' ability to adapt to a changing environment. In long-term studies of natural populations, individuals are often marked at one point in their life and then repeatedly recaptured. It is therefore essential that a method for GWAS includes the process of repeated sampling. In a GWAS, the effects of thousands of single-nucleotide polymorphisms (SNPs) need to be fitted and any model development is constrained by the computational requirements. A method is therefore required that can fit a highly hierarchical model and at the same time is computationally fast enough to be useful. 2. Our method fits fixed SNP effects in a linear mixed model that can include both random polygenic effects and permanent environmental effects. In this way, the model can correct for population structure and model repeated measures. The covariance structure of the linear mixed model is first estimated and subsequently used in a generalized least squares setting to fit the SNP effects. The method was evaluated in a simulation study based on observed genotypes from a long-term study of collared flycatchers in Sweden. 3. The method we present here was successful in estimating permanent environmental effects from simulated repeated measures data. Additionally, we found that especially for variable phenotypes having large variation between years, the repeated measurements model has a substantial increase in power compared to a model using average phenotypes as a response. 4. The method is available in the R package RepeatABEL. It increases the power in GWAS having repeated measures, especially for long-term studies of natural populations, and the R implementation is expected to facilitate modelling of longitudinal data for studies of both animal and human populations.