Markers of genetic susceptibility in human environmental hygiene and toxicology : the role of selected CYP, NAT and GST genes

Inherited genetic traits co-determine the susceptibility of an individual to a toxic chemical. Special emphasis has been put on individual responses to environmental and industrial carcinogens, but other chronic diseases are of increasing interest. Polymorphisms of relevant xenobiotic metabolising enzymes may be used as toxicological susceptibility markers. A growing number of genes encoding enzymes involved in biotransformation of toxicants and in cellular defence against toxicant-induced damage to the cells has been identified and cloned, leading to increased knowledge of allelic variants of genes and genetic defects that may result in a differential susceptibility toward environmental toxicants. "Low penetrating" polymorphisms in metabolism genes tend to be much more common in the population than allelic variants of "high penetrating" cancer genes, and are therefore of considerable importance from a public health point of view. Positive associations between cancer and CYP1A1 alleles, in particular the *2C I462V allele, were found for tissues following the aerodigestive tract. Again, in most cases, the effect of the variant CYP1A1 allele becomes apparent or clearer in connection with the GSTM1 null allele. The CYP1B1 codon 432 polymorphism (CYP1B1*3) has been identified as a susceptibility factor in smoking-related head-and-neck squameous cell cancer. The impact of this polymorphic variant of CYP1B1 on cancer risk was also reflected by an association with the frequency of somatic mutations of the p53 gene. Combined genotype analysis of CYP1B1 and the glutathione transferases GSTM1 or GSTT1 has also pointed to interactive effects. Of particular interest for the industrial and environmental field is the isozyme CYP2E1. Several genotypes of this isozyme have been characterised which seem to be associated with different levels of expression of enzyme activity. The acetylator status for NAT2 can be determined by genotyping or by phenotyping. In the pathogenesis of human bladder cancer due to occupational exposure to "classical" aromatic amines (benzidine, 4-aminodiphenyl, 1-naphthylamine) acetylation by NAT2 is regarded as a detoxication step. Interestingly, the underlying European findings of a higher susceptibility of slow acetylators towards aromatic amines are in contrast to findings in Chinese workers occupationally exposed to aromatic amines which points to different mechanisms of susceptibility between European and Chinese populations. Regarding human bladder cancer, the hypothesis has been put forward that genetic polymorphism of GSTM1 might be linked with the occurrence of this tumour type. This supports the hypothesis that exposure to PAH might causally be involved in urothelial cancers. The human polymorphic GST catalysing conjugation of halomethanes, dihalomethanes, ethylene oxide and a number of other industrial compounds could be characterised as a class theta enzyme (GSTT1) by means of molecular biology. "Conjugator" and "non-conjugator" phenotypes are coincident with the presence and absence of the GSTT1 gene. There are wide variations in the frequencies of GSTT1 deletion (GSTT1 *0/0) among different ethnicities. Human phenotyping is facilitated by the GST activity towards methyl bromide or ethylene oxide in erythrocytes which is representative of the metabolic GSTT1 competence of the entire organism. Inter-individual variations in xenobiotic metabolism capacities may be due to polymorphisms of the genes coding for the enzymes themselves or of the genes coding for the receptors or transcription factors which regulate the expression of the enzymes. Also, polymorphisms in several regions of genes may cause altered ligand affinity, transactivation activity or expression levels of the receptor subsequently influencing the expression of the downstream target genes. Studies of individual susceptibility to toxicants and gene-environment interaction are now emerging as an important component of molecular epidemiology.

Fine-mapping the HOXB region detects common variants tagging a rare coding allele : evidence for synthetic association in prostate cancer

The HOXB13 gene has been implicated in prostate cancer (PrCa) susceptibility. We performed a high resolution fine-mapping analysis to comprehensively evaluate the association between common genetic variation across the HOXB genetic locus at 17q21 and PrCa risk. This involved genotyping 700 SNPs using a custom Illumina iSelect array (iCOGS) followed by imputation of 3195 SNPs in 20,440 PrCa cases and 21,469 controls in The PRACTICAL consortium. We identified a cluster of highly correlated common variants situated within or closely upstream of HOXB13 that were significantly associated with PrCa risk, described by rs117576373 (OR 1.30, P = 2.62×10(-14)). Additional genotyping, conditional regression and haplotype analyses indicated that the newly identified common variants tag a rare, partially correlated coding variant in the HOXB13 gene (G84E, rs138213197), which has been identified recently as a moderate penetrance PrCa susceptibility allele. The potential for GWAS associations detected through common SNPs to be driven by rare causal variants with higher relative risks has long been proposed; however, to our knowledge this is the first experimental evidence for this phenomenon of synthetic association contributing to cancer susceptibility.

The effect of BCMO1 gene variants on macular pigment optical density in young healthy Caucasians

Background Serum lutein (L) and zeaxanthin (Z) positively correlate with macular pigment optical density (MPOD), hence the latter is a valuable indirect tool for measuring L and Z content in the macula. L and Z have been attributed antioxidant capacity and protection from certain retinal diseases but their uptake within the eye is thought to depend on genetic, age and environmental factors. In particular gene variants within beta-carotene monooxygenase (BCMO1) are thought to modulate MPOD in the macula. Objectives: To determine the effect of BCMO1 single nucleotide polymorphisms (SNPs) rs11645428, rs6420424 and rs6464851 on macular pigment optical density (MPOD) in a cohort of young healthy participants of Caucasian origin with normal ocular health. Design In this cohort study, MPOD was assessed in 46 healthy participants (22 male and 24 female) with a mean age of 24 ± 4.0 years (range 19-33). The three SNPs, rs11645428, rs6420424, rs6564851 that have established associations with MPOD were determined using MassEXTEND (hME) Sequenom assay. One-way analysis of variance (ANOVA) was performed on groups segregated into homozygous and heterozygous BCMO1 genotypes. Correlations between body mass index (BMI), iris colour, gender, central retinal thickness (CRT), diet and MPOD were investigated. Results MPOD did not significantly vary with BCMO1 rs11645428 (F2,41 = 0.700, p = 0.503), rs6420424 (F2,41 = 0.210, p = 0.801) nor rs6464851 homozygous or heterozygous genotypes (F2,41 = 0,13, p = 0.88), in this young healthy cohort. The combination of these three SNPs into triple genotypes based on plasma conversion efficiency did not affect MPOD (F2,41 = 0.07, p = 0.9). There was a significant negative correlation with MPOD and central retinal thickness (r = - 0.39, p = 0.01) but no significant correlation between BMI, iris colour, gender and MPOD. Conclusion Our results indicate that macular pigment deposition within the central retina is not dependent on BCMO1 gene variants in young healthy people. We propose that MPOD is saturated in younger persons and/or other gene variant combinations determine its deposition.

Fine-mapping of the HNF1B multicancer locus identifies candidate variants that mediate endometrial cancer risk

Common variants in the hepatocyte nuclear factor 1 homeobox B (HNF1B) gene are associated with the risk of Type II diabetes and multiple cancers. Evidence to date indicates that cancer risk may be mediated via genetic or epigenetic effects on HNF1B gene expression. We previously found single-nucleotide polymorphisms (SNPs) at the HNF1B locus to be associated with endometrial cancer, and now report extensive fine-mapping and in silico and laboratory analyses of this locus. Analysis of 1184 genotyped and imputed SNPs in 6608 Caucasian cases and 37 925 controls, and 895 Asian cases and 1968 controls, revealed the best signal of association for SNP rs11263763 (P = 8.4 × 10−14, odds ratio = 0.86, 95% confidence interval = 0.82–0.89), located within HNF1B intron 1. Haplotype analysis and conditional analyses provide no evidence of further independent endometrial cancer risk variants at this locus. SNP rs11263763 genotype was associated with HNF1B mRNA expression but not with HNF1B methylation in endometrial tumor samples from The Cancer Genome Atlas. Genetic analyses prioritized rs11263763 and four other SNPs in high-to-moderate linkage disequilibrium as the most likely causal SNPs. Three of these SNPs map to the extended HNF1B promoter based on chromatin marks extending from the minimal promoter region. Reporter assays demonstrated that this extended region reduces activity in combination with the minimal HNF1B promoter, and that the minor alleles of rs11263763 or rs8064454 are associated with decreased HNF1B promoter activity. Our findings provide evidence for a single signal associated with endometrial cancer risk at the HNF1B locus, and that risk is likely mediated via altered HNF1B gene expression.

Common variants in the CYP2C19 gene are associated with susceptibility to endometriosis

Objective: To follow-up previous studies highlighting a possible role for cytochrome P450, family 2, subfamily C, 19 (CYP2C19) in susceptibility to endometriosis by searching for additional variants in the CYP2C19 gene that may be associated with the disease. Design Case-control study. Setting Academic research. Subject(s) The cases comprised 2,271 women with surgically confirmed endometriosis; the controls comprised 939 women with self-report of no endometriosis and 1,770 unscreened population samples. Intervention(s) Sequencing of the CYP2C19 region and follow-up of 80 single nucleotide polymorphisms (SNPs) in two case-control samples. Main Outcome Measure(s) Allele frequency differences between cases and controls. Result(s) Sequencing of the CYP2C19 gene region resulted in the detection of a large number of known and novel SNPs. Genotyping of 80 polymorphic SNPs in 901 endometriosis cases and 939 controls resulted in study-wide significant association signals for SNPs in moderate or complete linkage disequilibrium with rs4244285, a functional SNP in exon 5 that abrogates CYP2C19 function through the creation of an alternative splice site. Evidence of association was also detected for another functional SNP in the CYP2C19 promoter, rs12248560, which was highlighted in our previous study. Conclusion(s) Functional variants in CYP2C19 may contribute to endometriosis susceptibility in both familial and sporadic cases. © 2014 by American Society for Reproductive Medicine.

The future for genetic studies in reproduction

Genetic factors contribute to risk of many common diseases affecting reproduction and fertility. In recent years, methods for genome-wide association studies(GWAS) have revolutionized gene discovery forcommontraits and diseases. Results of GWAS are documented in the Catalog of Published Genome-Wide Association Studies at the National Human Genome Research Institute and report over 70 publications for 32 traits and diseases associated with reproduction. These include endometriosis, uterine fibroids, age at menarche and age at menopause. Results that pass appropriate stringent levels of significance are generally well replicated in independent studies. Examples of genetic variation affecting twinning rate, infertility, endometriosis and age at menarche demonstrate that the spectrum of disease-related variants for reproductive traits is similar to most other common diseases.GWAS 'hits' provide novel insights into biological pathways and the translational value of these studies lies in discovery of novel gene targets for biomarkers, drug development and greater understanding of environmental factors contributing to disease risk. Results also show that genetic data can help define sub-types of disease and co-morbidity with other traits and diseases. To date, many studies on reproductive traits have used relatively small samples. Future genetic marker studies in large samples with detailed phenotypic and clinical information will yield new insights into disease risk, disease classification and co-morbidity for many diseases associated with reproduction and infertility.

Genetic clustering on the hippocampal surface for genome-wide association studies

Imaging genetics aims to discover how variants in the human genome influence brain measures derived from images. Genome-wide association scans (GWAS) can screen the genome for common differences in our DNA that relate to brain measures. In small samples, GWAS has low power as individual gene effects are weak and one must also correct for multiple comparisons across the genome and the image. Here we extend recent work on genetic clustering of images, to analyze surface-based models of anatomy using GWAS. We performed spherical harmonic analysis of hippocampal surfaces, automatically extracted from brain MRI scans of 1254 subjects. We clustered hippocampal surface regions with common genetic influences by examining genetic correlations (r(g)) between the normalized deformation values at all pairs of surface points. Using genetic correlations to cluster surface measures, we were able to boost effect sizes for genetic associations, compared to clustering with traditional phenotypic correlations using Pearson's r.

The ENIGMA Consortium: Large-scale collaborative analyses of neuroimaging and genetic data

The Enhancing NeuroImaging Genetics through Meta-Analysis (ENIGMA) Consortium is a collaborative network of researchers working together on a range of large-scale studies that integrate data from 70 institutions worldwide. Organized into Working Groups that tackle questions in neuroscience, genetics, and medicine, ENIGMA studies have analyzed neuroimaging data from over 12,826 subjects. In addition, data from 12,171 individuals were provided by the CHARGE consortium for replication of findings, in a total of 24,997 subjects. By meta-analyzing results from many sites, ENIGMA has detected factors that affect the brain that no individual site could detect on its own, and that require larger numbers of subjects than any individual neuroimaging study has currently collected. ENIGMA's first project was a genome-wide association study identifying common variants in the genome associated with hippocampal volume or intracranial volume. Continuing work is exploring genetic associations with subcortical volumes (ENIGMA2) and white matter microstructure (ENIGMA-DTI). Working groups also focus on understanding how schizophrenia, bipolar illness, major depression and attention deficit/hyperactivity disorder (ADHD) affect the brain. We review the current progress of the ENIGMA Consortium, along with challenges and unexpected discoveries made on the way.

Case-control study of ADARB1 and ADARB2 gene variants in migraine

Background: Migraine causes crippling attacks of severe head pain along with associated nausea, vomiting, photophobia and/or phonophobia. The aim of this study was to investigate single nucleotide polymorphisms (SNPs) in the adenosine deaminase, RNA-specific, B1 (ADARB1)and adenosine deaminase, RNA specific, B2 (ADARB2) genes in an Australian case-control Caucasian population for association with migraine. Both candidate genes are highly expressed in the central nervous system (CNS) and fit criteria for migraine neuropathology. SNPs in the ADARB2 gene were previously found to be positively associated with migraine in a pedigree-based GWAS using the genetic isolate of Norfolk Island, Australia. The ADARB1 gene was also chosen for investigation due to its important function in editing neurotransmitter receptor transcripts. Methods: Four SNPs in ADARB1 and nine in ADARB2 were selected by inspecting blocks of LD in Haploview for genotyping using either TaqMan or Sequenom assays. These SNPs were genotyped in two-hundred and ninety one patients who satisfied the International Classification of Headache Disorders, ICHD-II 2004 diagnostic criteria for migraine and three-hundred and fourteen controls and PLINK was used for association testing. Results: Chi-square (χ2) analysis found no significant association between any of the SNPs tested in the ADARB1 and ADARB2 genes in this study and the occurrence of migraine. Conclusions: In contrast to findings that SNPs in the ADARB2 gene were positively associated with migraine in the Norfolk Island population, we find no evidence to support the involvement of RNA editing genes in migraine susceptibility in an Australian Caucasian population.

Identity-by-Descent Mapping to Detect Rare Variants Conferring Susceptibility to Multiple Sclerosis

Genome-wide association studies (GWAS) have identified around 60 common variants associated with multiple sclerosis (MS), but these loci only explain a fraction of the heritability of MS. Some missing heritability may be caused by rare variants that have been suggested to play an important role in the aetiology of complex diseases such as MS. However current genetic and statistical methods for detecting rare variants are expensive and time consuming. 'Population-based linkage analysis' (PBLA) or so called identity-by-descent (IBD) mapping is a novel way to detect rare variants in extant GWAS datasets. We employed BEAGLE fastIBD to search for rare MS variants utilising IBD mapping in a large GWAS dataset of 3,543 cases and 5,898 controls. We identified a genome-wide significant linkage signal on chromosome 19 (LOD = 4.65; p = 1.9×10-6). Network analysis of cases and controls sharing haplotypes on chromosome 19 further strengthened the association as there are more large networks of cases sharing haplotypes than controls. This linkage region includes a cluster of zinc finger genes of unknown function. Analysis of genome wide transcriptome data suggests that genes in this zinc finger cluster may be involved in very early developmental regulation of the CNS. Our study also indicates that BEAGLE fastIBD allowed identification of rare variants in large unrelated population with moderate computational intensity. Even with the development of whole-genome sequencing, IBD mapping still may be a promising way to narrow down the region of interest for sequencing priority. © 2013 Lin et al.

Association of variants in MMEL1 and CTLA4 with rheumatoid arthritis in the Han Chinese population

Background: The genome-wide association study era has made great progress in identifying susceptibility genes and genetic loci for rheumatoid arthritis (RA) in populations of White European ancestry. However, few studies have tried to dissect disease aetiopathogenesis in other ethnic populations. Objective: To investigate these associations in the Han Chinese population. Methods: Haplotypes from the HapMap database Chinese population were used to select tag-single-nucleotide polymorphisms (SNPs) (r2 =0.8) across 19 distinct RA genomic regions. A two phase case-control association study was performed, with 169 SNPs genotyped in phase I (n=571 cases, n=880 controls), and 64 SNPs achieving p<0.2 in the first phase being genotyped in phase II (n=464 cases, n=822 controls). Association statistics were calculated using permutation tests both unadjusted and adjusted for the number of markers studied. Results: Robust association was detected for MMEL1 and CTLA4 , and modest association was identified for another six loci: PADI4 , STAT4 , PRDM1 , CDK6 , TRAF1-C5 and KIF5A-PIP4K2C. All three markers genotyped in MMEL1 demonstrated association, with peak signal for rs3890745 (p=2.6×10 -5unadjusted, p=0.003 adjusted, OR=0.79). For CTLA4 , significance was detected for three of five variants showing association, with peak association for marker rs12992492 (p=4.3×10-5 unadjusted, p=0.0021 adjusted, OR=0.77). Lack of association of common variants in PTPN22 with RA in Han Chinese was confirmed. Conclusion: This study identifies MMEL1 and CTLA4 as RA susceptibility genes, provides suggestive evidence of association for a further six loci in the Han Chinese population and confirms lack of PTPN22 association in Asian populations. It also confirms the value of multiethnic population studies to help dissect disease aetiopathogenesis.

Role of NOD2 variants in spondylarthritis

Objective. To investigate the role of the gene NOD2 in susceptibility to, and clinical manifestations of, ankylosing spondylitis (AS). Methods. A case-control study of NOD2 polymorphisms known to be associated with Crohn's disease (CD) (Pro268 Ser, Arg702 Trp, GlY908 Arg, and Len1007fsinsC) was performed in 229 cases of primary AS with no diagnosed inflammatory bowel disease (IBD), 197 cases of AS associated with IBD (referred to as colitic spondylarthritis; comprising 78 with CD and 119 with ulcerative colitis [UC]), and 229 ethnically matched, healthy controls. Associations between NOD2 polymorphisms and several clinical features of AS, including disease severity assessed by questionnaire and age at spondylarthritis onset, were also investigated. Exclusion linkage mapping of chromosome 16 was performed in a separate group of 185 multicase families with AS. Results. An association was identified between Gly908 Arg and UC spondylarthritis (P = 0.016, odds ratio [OR] 4.6, 95% confidence interval [95% CI] 1.316), and a nonsignificant trend with a similar magnitude was observed in association with CD spondylarthritis (P = 0.08, OR 3.9, 95% CI 0.8-18). The Pro268Ser variant was inversely associated with UC spondylarthritis (P = 0.003, OR 0.55, 95% CI 0.37-0.82), but not with CD spondylarthritis. No association was demonstrated between NOD2 variants and primary AS, or between other variants of NOD2 and either UC or CD spondylarthritis. Carriage of the Pro268 Ser polymorphism was associated with greater disease activity as measured by the Bath Ankylosing Spondylitis Disease Activity Index (P = 0.002). Although patients with CD had a younger age at spondylarthritis onset than did those with UC (22.4 years versus 26.4 years; P = 0.01), no association was noted between the NOD2 variants linked with CD and age at spondylarthritis onset. In primary AS, the presence of a gene with a magnitude of association >2.0 was excluded (exclusion logarithm of odds score less than -2.0), and no association was observed with the microsatellite D16S3136. Conclusion. NOD2 variants do not significantly affect the risk of developing primary AS, but may influence susceptibility to, and clinical manifestations of, colitic spondylarthritis.

Is there a higher genetic load of susceptibility loci in familial ankylosing spondylitis?

Objective Several genetic risk variants for ankylosing spondylitis (AS) have been identified in genome-wide association studies. Our objective was to examine whether familial AS cases have a higher genetic load of these susceptibility variants. Methods Overall, 502 AS patients were examined, consisting of 312 patients who had first-degree relatives (FDRs) with AS (familial) and 190 patients who had no FDRs with AS or spondylarthritis (sporadic). All patients and affected FDRs fulfilled the modified New York criteria for AS. The patients were recruited from 2 US cohorts (the North American Spondylitis Consortium and the Prospective Study of Outcomes in Ankylosing Spondylitis) and from the UK-Oxford cohort. The frequencies of AS susceptibility loci in IL-23R, IL1R2, ANTXR2, ERAP-1, 2 intergenic regions on chromosomes 2p15 and 21q22, and HLA-B27 status as determined by the tag single-nucleotide polymorphism (SNP) rs4349859 were compared between familial and sporadic cases of AS. Association between SNPs and multiplex status was assessed by logistic regression controlling for sibship size. Results HLA-B27 was significantly more prevalent in familial than sporadic cases of AS (odds ratio 4.44 [95% confidence interval 2.06, 9.55], P = 0.0001). Furthermore, the AS risk allele at chromosome 21q22 intergenic region showed a trend toward higher frequency in the multiplex cases (P = 0.08). The frequency of the other AS risk variants did not differ significantly between familial and sporadic cases, either individually or combined. Conclusion HLA-B27 is more prevalent in familial than sporadic cases of AS, demonstrating higher familial aggregation of AS in patients with HLA-B27 positivity. The frequency of the recently described non-major histocompatibility complex susceptibility loci is not markedly different between the sporadic and familial cases of AS.

Evidence of genetic association between TNFRSF1A encoding the p55 tumour necrosis factor receptor, and ankylosing spondylitis in UK Caucasians

Objectives: To replicate the possible genetic association between ankylosing spondylitis (AS) and TNFRSF1A. Methods: TNFRSF1A was re-sequenced in 48 individuals with AS to identify novel polymorphisms. Nine single nucleotide polymorphisms (SNPs) in TNFRSF1A and 5 SNPs in the neighbouring gene SCNN1A were genotyped in 1604 UK Caucasian individuals with AS and 1019 matched controls. An extended study was implemented using additional genotype data on 8 of these SNPs from 1400 historical controls from the 1958 British Birth Cohort. A meta-analysis of previously published results was also undertaken. Results: One novel variant in intron 6 was identified but no new coding variants. No definite associations were seen in the initial study but in the extended study there were weak associations with rs4149576 (p=0.04) and rs4149577 (p=0.007). In the metaanalysis consistent, somewhat stronger associations were seen with rs4149577 (p=0.002) and rs4149578 (p=0.006). Conclusions: These studies confirm the weak genetic associations between AS and TNFRSF1A. In view of the previously reported associations of TNFRSF1A with AS, in Caucasians and Chinese, and the biological plausibility of this candidate gene, replication of this finding in well powered studies is clearly indicated.

Genome-wide association analysis identifies 11 risk variants associated with the asthma with hay fever phenotype

Background To date, no genome-wide association study (GWAS) has considered the combined phenotype of asthma with hay fever. Previous analyses of family data from the Tasmanian Longitudinal Health Study provide evidence that this phenotype has a stronger genetic cause than asthma without hay fever. Objective We sought to perform a GWAS of asthma with hay fever to identify variants associated with having both diseases. Methods We performed a meta-analysis of GWASs comparing persons with both physician-diagnosed asthma and hay fever (n = 6,685) with persons with neither disease (n = 14,091). Results At genome-wide significance, we identified 11 independent variants associated with the risk of having asthma with hay fever, including 2 associations reaching this level of significance with allergic disease for the first time: ZBTB10 (rs7009110; odds ratio [OR], 1.14; P = 4 × 10−9) and CLEC16A (rs62026376; OR, 1.17; P = 1 × 10−8). The rs62026376:C allele associated with increased asthma with hay fever risk has been found to be associated also with decreased expression of the nearby DEXI gene in monocytes. The 11 variants were associated with the risk of asthma and hay fever separately, but the estimated associations with the individual phenotypes were weaker than with the combined asthma with hay fever phenotype. A variant near LRRC32 was a stronger risk factor for hay fever than for asthma, whereas the reverse was observed for variants in/near GSDMA and TSLP. Single nucleotide polymorphisms with suggestive evidence for association with asthma with hay fever risk included rs41295115 near IL2RA (OR, 1.28; P = 5 × 10−7) and rs76043829 in TNS1 (OR, 1.23; P = 2 × 10−6). Conclusion By focusing on the combined phenotype of asthma with hay fever, variants associated with the risk of allergic disease can be identified with greater efficiency.