Genome-wide association analyses of physical activity and sedentary behavior provide insights into underlying mechanisms and roles in disease prevention

Although physical activity and sedentary behavior are moderately heritable, little is known about the mechanisms that influence these traits. Combining data for up to 703,901 individuals from 51 studies in a multi-ancestry meta-analysis of genome-wide association studies yields 99 loci that associate with self-reported moderate-to-vigorous intensity physical activity during leisure time (MVPA), leisure screen time (LST) and/or sedentary behavior at work. Loci associated with LST are enriched for genes whose expression in skeletal muscle is altered by resistance training. A missense variant in ACTN3 makes the alpha-actinin-3 filaments more flexible, resulting in lower maximal force in isolated type IIA muscle fibers, and possibly protection from exercise-induced muscle damage. Finally, Mendelian randomization analyses show that beneficial effects of lower LST and higher MVPA on several risk factors and diseases are mediated or confounded by body mass index (BMI). Our results provide insights into physical activity mechanisms and its role in disease prevention.

Genome-wide meta-analysis of 241,258 adults accounting for smoking behaviour identifies novel loci for obesity traits

Few genome-wide association studies (GWAS) account for environmental exposures, like smoking, potentially impacting the overall trait variance when investigating the genetic contribution to obesity-related traits. Here, we use GWAS data from 51,080 current smokers and 190,178 nonsmokers (87% European descent) to identify loci influencing BMI and central adiposity, measured as waist circumference and waist-to-hip ratio both adjusted for BMI. We identify 23 novel genetic loci, and 9 loci with convincing evidence of gene-smoking interaction (GxSMK) on obesity-related traits. We show consistent direction of effect for all identified loci and significance for 18 novel and for 5 interaction loci in an independent study sample. These loci highlight novel biological functions, including response to oxidative stress, addictive behaviour, and regulatory functions emphasizing the importance of accounting for environment in genetic analyses. Our results suggest that tobacco smoking may alter the genetic susceptibility to overall adiposity and body fat distribution.

Missing heritability: is the gap closing? An analysis of 32 complex traits in the Lifelines Cohort Study

Despite the recent explosive rise in number of genetic markers for complex disease traits identified in genome-wide association studies, there is still a large gap between the known heritability of these traits and the part explained by these markers. To gauge whether this 'heritability gap' is closing, we first identified genome-wide significant SNPs from the literature and performed replication analyses for 32 highly relevant traits from five broad disease areas in 13 436 subjects of the Lifelines Cohort. Next, we calculated the variance explained by multi-SNP genetic risk scores (GRSs) for each trait, and compared it to their broad- and narrow-sense heritabilities captured by all common SNPs. The majority of all previously-associated SNPs (median=75%) were significantly associated with their respective traits. All GRSs were significant, with unweighted GRSs generally explaining less phenotypic variance than weighted GRSs, for which the explained variance was highest for height (15.5%) and varied between 0.02 and 6.7% for the other traits. Broad-sense common-SNP heritability estimates were significant for all traits, with the additive effect of common SNPs explaining 48.9% of the variance for height and between 5.6 and 39.2% for the other traits. Dominance effects were uniformly small (0-1.5%) and not significant. On average, the variance explained by the weighted GRSs accounted for only 10.7% of the common-SNP heritability of the 32 traits. These results indicate that GRSs may not yet be ready for accurate personalized prediction of complex disease traits limiting widespread adoption in clinical practice.

A principal component meta-analysis on multiple anthropometric traits identifies novel loci for body shape

Large consortia have revealed hundreds of genetic loci associated with anthropometric traits, one trait at a time. We examined whether genetic variants affect body shape as a composite phenotype that is represented by a combination of anthropometric traits. We developed an approach that calculates averaged PCs (AvPCs) representing body shape derived from six anthropometric traits (body mass index, height, weight, waist and hip circumference, waist-to-hip ratio). The first four AvPCs explain >99% of the variability, are heritable, and associate with cardiometabolic outcomes. We performed genome-wide association analyses for each body shape composite phenotype across 65 studies and meta-analysed summary statistics. We identify six novel loci: LEMD2 and CD47 for AvPC1, RPS6KA5/C14orf159 and GANAB for AvPC3, and ARL15 and ANP32 for AvPC4. Our findings highlight the value of using multiple traits to define complex phenotypes for discovery, which are not captured by single-trait analyses, and may shed light onto new pathways.

Correction: The Influence of Age and Sex on Genetic Associations with Adult Body Size and Shape: A Large-Scale Genome-Wide Interaction Study

[This corrects the article DOI: 10.1371/journal.pgen.1005378.].

Genome-wide association study of lifetime cannabis use based on a large meta-analytic sample of 32 330 subjects from the International Cannabis Consortium

Cannabis is the most widely produced and consumed illicit psychoactive substance worldwide. Occasional cannabis use can progress to frequent use, abuse and dependence with all known adverse physical, psychological and social consequences. Individual differences in cannabis initiation are heritable (40-48%). The International Cannabis Consortium was established with the aim to identify genetic risk variants of cannabis use. We conducted a meta-analysis of genome-wide association data of 13 cohorts (N=32 330) and four replication samples (N=5627). In addition, we performed a gene-based test of association, estimated single-nucleotide polymorphism (SNP)-based heritability and explored the genetic correlation between lifetime cannabis use and cigarette use using LD score regression. No individual SNPs reached genome-wide significance. Nonetheless, gene-based tests identified four genes significantly associated with lifetime cannabis use: NCAM1, CADM2, SCOC and KCNT2. Previous studies reported associations of NCAM1 with cigarette smoking and other substance use, and those of CADM2 with body mass index, processing speed and autism disorders, which are phenotypes previously reported to be associated with cannabis use. Furthermore, we showed that, combined across the genome, all common SNPs explained 13-20% (P<0.001) of the liability of lifetime cannabis use. Finally, there was a strong genetic correlation (rg=0.83; P=1.85 × 10(-8)) between lifetime cannabis use and lifetime cigarette smoking implying that the SNP effect sizes of the two traits are highly correlated. This is the largest meta-analysis of cannabis GWA studies to date, revealing important new insights into the genetic pathways of lifetime cannabis use. Future functional studies should explore the impact of the identified genes on the biological mechanisms of cannabis use.

The Influence of Age and Sex on Genetic Associations with Adult Body Size and Shape: A Large-Scale Genome-Wide Interaction Study

Genome-wide association studies (GWAS) have identified more than 100 genetic variants contributing to BMI, a measure of body size, or waist-to-hip ratio (adjusted for BMI, WHR_adjBMI), a measure of body shape. Body size and shape change as people grow older and these changes differ substantially between men and women. To systematically screen for age- and/or sex-specific effects of genetic variants on BMI and WHR_adjBMI, we performed meta-analyses of 114 studies (up to 320,485 individuals of European descent) with genome-wide chip and/or Metabochip data by the Genetic Investigation of Anthropometric Traits (GIANT) Consortium. Each study tested the association of up to ~2.8M SNPs with BMI and WHR_adjBMI in four strata (men ≤50y, men >50y, women ≤50y, women >50y) and summary statistics were combined in stratum-specific meta-analyses. We then screened for variants that showed age-specific effects (G x AGE), sex-specific effects (G x SEX) or age-specific effects that differed between men and women (G x AGE x SEX). For BMI, we identified 15 loci (11 previously established for main effects, four novel) that showed significant (FDR<5%) age-specific effects, of which 11 had larger effects in younger (<50y) than in older adults (≥50y). No sex-dependent effects were identified for BMI. For WHR_adjBMI, we identified 44 loci (27 previously established for main effects, 17 novel) with sex-specific effects, of which 28 showed larger effects in women than in men, five showed larger effects in men than in women, and 11 showed opposite effects between sexes. No age-dependent effects were identified for WHR_adjBMI. This is the first genome-wide interaction meta-analysis to report convincing evidence of age-dependent genetic effects on BMI. In addition, we confirm the sex-specificity of genetic effects on WHR_adjBMI. These results may provide further insights into the biology that underlies weight change with age or the sexually dimorphism of body shape.

Adult body size and body shape differ substantially between men and women and change over time. More than 100 genetic variants that influence body mass index (measure of body size) or waist-to-hip ratio (measure of body shape) have been identified. While there is evidence that some genetic loci affect body shape differently in men than in women, little is known about whether genetic effects differ in older compared to younger adults, and whether such changes differ between men and women. Therefore, we conducted a systematic genome-wide search, including 114 studies (>320,000 individuals), to specifically identify genetic loci with age- and or sex-dependent effects on body size and shape. We identified 15 loci of which the effect on BMI was different in older compared to younger adults, whereas we found no evidence for loci with different effects in men compared to women. The opposite was seen for body shape as we identified 44 loci of which the effect on waist-to-hip ratio differed between men and women, but no difference between younger and older adults were observed. Our observations may provide new insights into the biology that underlies weight change with age or the sexual dimorphism of body shape.

New genetic loci link adipose and insulin biology to body fat distribution

Body fat distribution is a heritable trait and a well-established predictor of adverse metabolic outcomes, independent of overall adiposity. To increase our understanding of the genetic basis of body fat distribution and its molecular links to cardiometabolic traits, we conducted genome-wide association meta-analyses of waist and hip circumference-related traits in up to 224,459 individuals. We identified 49 loci (33 new) associated with waist-to-hip ratio adjusted for body mass index (WHRadjBMI) and an additional 19 loci newly associated with related waist and hip circumference measures (P<5×10⁻⁸). Twenty of the 49 WHRadjBMI loci showed significant sexual dimorphism, 19 of which displayed a stronger effect in women. The identified loci were enriched for genes expressed in adipose tissue and for putative regulatory elements in adipocytes. Pathway analyses implicated adipogenesis, angiogenesis, transcriptional regulation, and insulin resistance as processes affecting fat distribution, providing insight into potential pathophysiological mechanisms.

Defining the role of common variation in the genomic and biological architecture of adult human height

Using genome-wide data from 253,288 individuals, we identified 697 variants at genome-wide significance that together explained one-fifth of the heritability for adult height. By testing different numbers of variants in independent studies, we show that the most strongly associated ∼2,000, ∼3,700 and ∼9,500 SNPs explained ∼21%, ∼24% and ∼29% of phenotypic variance. Furthermore, all common variants together captured 60% of heritability. The 697 variants clustered in 423 loci were enriched for genes, pathways and tissue types known to be involved in growth and together implicated genes and pathways not highlighted in earlier efforts, such as signaling by fibroblast growth factors, WNT/β-catenin and chondroitin sulfate-related genes. We identified several genes and pathways not previously connected with human skeletal growth, including mTOR, osteoglycin and binding of hyaluronic acid. Our results indicate a genetic architecture for human height that is characterized by a very large but finite number (thousands) of causal variants.

Comparing the hub-and-spoke model practices of the LifeLines Study in the Netherlands and the H3Africa Initiative

Biobank sample storage is critical in population health and epidemiology studies. Biobanks bridge two very different worlds: they connect to the participants and patients at an individual level, but they also aggregate information and represent the cutting edge of scientific discovery. In this brief report, we describe how the LifeLines study in the Netherlands manages its resources for communication and services, and how it can serve as a model for the Human Heredity and Health in Africa Initiative (H3Africa Initiative).

Common genes underlying asthma and COPD? Genome-wide analysis on the Dutch hypothesis

Asthma and chronic obstructive pulmonary disease (COPD) are thought to share a genetic background ("Dutch hypothesis"). We investigated whether asthma and COPD have common underlying genetic factors, performing genome-wide association studies for both asthma and COPD and combining the results in meta-analyses. Three loci showed potential involvement in both diseases: chr2p24.3, chr5q23.1 and chr13q14.2, containing DDX1, COMMD10 (both participating in the nuclear factor (NF) κβ pathway) and GNG5P5, respectively. Single nucleotide polymorphisms (SNPs) rs9534578 in GNG5P5 reached genome-wide significance after first replication phase (p=9.96×10(-9)). The second replication phase, in seven independent cohorts, provided no significant replication. Expression quantitative trait loci (eQTL) analysis in blood cells and lung tissue on the top 20 associated SNPs identified two SNPs in COMMD10 that influenced gene expression. Inflammatory processes differ in asthma and COPD and are mediated by NF-κβ, which could be driven by the same underlying genes, COMMD10 and DDX1. None of the SNPs reached genome-wide significance. Our eQTL studies support a functional role for two COMMD10 SNPs, since they influence gene expression in both blood cells and lung tissue. Our findings suggest that there is either no common genetic component in asthma and COPD or, alternatively, different environmental factors, e.g. lifestyle and occupation in different countries and continents, which may have obscured the genetic common contribution.

Pleiotropic effects of obesity-susceptibility loci on metabolic traits: a meta-analysis of up to 37,874 individuals

AIMS/HYPOTHESIS

Genetic pleiotropy may contribute to the clustering of obesity and metabolic conditions. We assessed whether genetic variants that are robustly associated with BMI and waist-to-hip ratio (WHR) also influence metabolic and cardiovascular traits, independently of obesity-related traits, in meta-analyses of up to 37,874 individuals from six European population-based studies.

METHODS

We examined associations of 32 BMI and 14 WHR loci, individually and combined in two genetic predisposition scores (GPSs), with glycaemic traits, blood lipids and BP, with and without adjusting for BMI and/or WHR.

RESULTS

We observed significant associations of BMI-increasing alleles at five BMI loci with lower levels of 2 h glucose (RBJ [also known as DNAJC27], QPTCL: effect sizes -0.068 and -0.107 SD, respectively), HDL-cholesterol (SLC39A8: -0.065 SD, MTCH2: -0.039 SD), and diastolic BP (SLC39A8: -0.069 SD), and higher and lower levels of LDL- and total cholesterol (QPTCL: 0.041 and 0.042 SDs, respectively, FLJ35779 [also known as POC5]: -0.042 and -0.041 SDs, respectively) (all p < 2.4 × 10(-4)), independent of BMI. The WHR-increasing alleles at two WHR loci were significantly associated with higher proinsulin (GRB14: 0.069 SD) and lower fasting glucose levels (CPEB4: -0.049 SD), independent of BMI and WHR. A higher GPS-BMI was associated with lower systolic BP (-0.005 SD), diastolic BP (-0.006 SD) and 2 h glucose (-0.013 SD), while a higher GPS-WHR was associated with lower HDL-cholesterol (-0.015 SD) and higher triacylglycerol levels (0.014 SD) (all p < 2.9 × 10(-3)), independent of BMI and/or WHR.

CONCLUSIONS/INTERPRETATION

These pleiotropic effects of obesity-susceptibility loci provide novel insights into mechanisms that link obesity with metabolic abnormalities.

Meta-analysis of gene-level associations for rare variants based on single-variant statistics

Meta-analysis of genome-wide association studies (GWASs) has led to the discoveries of many common variants associated with complex human diseases. There is a growing recognition that identifying "causal" rare variants also requires large-scale meta-analysis. The fact that association tests with rare variants are performed at the gene level rather than at the variant level poses unprecedented challenges in the meta-analysis. First, different studies may adopt different gene-level tests, so the results are not compatible. Second, gene-level tests require multivariate statistics (i.e., components of the test statistic and their covariance matrix), which are difficult to obtain. To overcome these challenges, we propose to perform gene-level tests for rare variants by combining the results of single-variant analysis (i.e., p values of association tests and effect estimates) from participating studies. This simple strategy is possible because of an insight that multivariate statistics can be recovered from single-variant statistics, together with the correlation matrix of the single-variant test statistics, which can be estimated from one of the participating studies or from a publicly available database. We show both theoretically and numerically that the proposed meta-analysis approach provides accurate control of the type I error and is as powerful as joint analysis of individual participant data. This approach accommodates any disease phenotype and any study design and produces all commonly used gene-level tests. An application to the GWAS summary results of the Genetic Investigation of ANthropometric Traits (GIANT) consortium reveals rare and low-frequency variants associated with human height. The relevant software is freely available.

Genome-wide meta-analysis identifies 11 new loci for anthropometric traits and provides insights into genetic architecture

Approaches exploiting trait distribution extremes may be used to identify loci associated with common traits, but it is unknown whether these loci are generalizable to the broader population. In a genome-wide search for loci associated with the upper versus the lower 5th percentiles of body mass index, height and waist-to-hip ratio, as well as clinical classes of obesity, including up to 263,407 individuals of European ancestry, we identified 4 new loci (IGFBP4, H6PD, RSRC1 and PPP2R2A) influencing height detected in the distribution tails and 7 new loci (HNF4G, RPTOR, GNAT2, MRPS33P4, ADCY9, HS6ST3 and ZZZ3) for clinical classes of obesity. Further, we find a large overlap in genetic structure and the distribution of variants between traits based on extremes and the general population and little etiological heterogeneity between obesity subgroups.

Genome-wide association analyses identify 18 new loci associated with serum urate concentrations

Elevated serum urate concentrations can cause gout, a prevalent and painful inflammatory arthritis. By combining data from >140,000 individuals of European ancestry within the Global Urate Genetics Consortium (GUGC), we identified and replicated 28 genome-wide significant loci in association with serum urate concentrations (18 new regions in or near TRIM46, INHBB, SFMBT1, TMEM171, VEGFA, BAZ1B, PRKAG2, STC1, HNF4G, A1CF, ATXN2, UBE2Q2, IGF1R, NFAT5, MAF, HLF, ACVR1B-ACVRL1 and B3GNT4). Associations for many of the loci were of similar magnitude in individuals of non-European ancestry. We further characterized these loci for associations with gout, transcript expression and the fractional excretion of urate. Network analyses implicate the inhibins-activins signaling pathways and glucose metabolism in systemic urate control. New candidate genes for serum urate concentration highlight the importance of metabolic control of urate production and excretion, which may have implications for the treatment and prevention of gout.

Large-scale association analyses identify new loci influencing glycemic traits and provide insight into the underlying biological pathways

Through genome-wide association meta-analyses of up to 133,010 individuals of European ancestry without diabetes, including individuals newly genotyped using the Metabochip, we have increased the number of confirmed loci influencing glycemic traits to 53, of which 33 also increase type 2 diabetes risk (q < 0.05). Loci influencing fasting insulin concentration showed association with lipid levels and fat distribution, suggesting impact on insulin resistance. Gene-based analyses identified further biologically plausible loci, suggesting that additional loci beyond those reaching genome-wide significance are likely to represent real associations. This conclusion is supported by an excess of directionally consistent and nominally significant signals between discovery and follow-up studies. Functional analysis of these newly discovered loci will further improve our understanding of glycemic control.

Large-scale gene-centric meta-analysis across 39 studies identifies type 2 diabetes loci

To identify genetic factors contributing to type 2 diabetes (T2D), we performed large-scale meta-analyses by using a custom ∼50,000 SNP genotyping array (the ITMAT-Broad-CARe array) with ∼2000 candidate genes in 39 multiethnic population-based studies, case-control studies, and clinical trials totaling 17,418 cases and 70,298 controls. First, meta-analysis of 25 studies comprising 14,073 cases and 57,489 controls of European descent confirmed eight established T2D loci at genome-wide significance. In silico follow-up analysis of putative association signals found in independent genome-wide association studies (including 8,130 cases and 38,987 controls) performed by the DIAGRAM consortium identified a T2D locus at genome-wide significance (GATAD2A/CILP2/PBX4; p = 5.7 × 10(-9)) and two loci exceeding study-wide significance (SREBF1, and TH/INS; p < 2.4 × 10(-6)). Second, meta-analyses of 1,986 cases and 7,695 controls from eight African-American studies identified study-wide-significant (p = 2.4 × 10(-7)) variants in HMGA2 and replicated variants in TCF7L2 (p = 5.1 × 10(-15)). Third, conditional analysis revealed multiple known and novel independent signals within five T2D-associated genes in samples of European ancestry and within HMGA2 in African-American samples. Fourth, a multiethnic meta-analysis of all 39 studies identified T2D-associated variants in BCL2 (p = 2.1 × 10(-8)). Finally, a composite genetic score of SNPs from new and established T2D signals was significantly associated with increased risk of diabetes in African-American, Hispanic, and Asian populations. In summary, large-scale meta-analysis involving a dense gene-centric approach has uncovered additional loci and variants that contribute to T2D risk and suggests substantial overlap of T2D association signals across multiple ethnic groups.

Comparative analysis of the human hepatic and adipose tissue transcriptomes during LPS-induced inflammation leads to the identification of differential biological pathways and candidate biomarkers

BACKGROUND

Insulin resistance (IR) is accompanied by chronic low grade systemic inflammation, obesity, and deregulation of total body energy homeostasis. We induced inflammation in adipose and liver tissues in vitro in order to mimic inflammation in vivo with the aim to identify tissue-specific processes implicated in IR and to find biomarkers indicative for tissue-specific IR.

METHODS

Human adipose and liver tissues were cultured in the absence or presence of LPS and DNA Microarray Technology was applied for their transcriptome analysis. Gene Ontology (GO), gene functional analysis, and prediction of genes encoding for secretome were performed using publicly available bioinformatics tools (DAVID, STRING, SecretomeP). The transcriptome data were validated by proteomics analysis of the inflamed adipose tissue secretome.

RESULTS

LPS treatment significantly affected 667 and 483 genes in adipose and liver tissues respectively. The GO analysis revealed that during inflammation adipose tissue, compared to liver tissue, had more significantly upregulated genes, GO terms, and functional clusters related to inflammation and angiogenesis. The secretome prediction led to identification of 399 and 236 genes in adipose and liver tissue respectively. The secretomes of both tissues shared 66 genes and the remaining genes were the differential candidate biomarkers indicative for inflamed adipose or liver tissue. The transcriptome data of the inflamed adipose tissue secretome showed excellent correlation with the proteomics data.

CONCLUSIONS

The higher number of altered proinflammatory genes, GO processes, and genes encoding for secretome during inflammation in adipose tissue compared to liver tissue, suggests that adipose tissue is the major organ contributing to the development of systemic inflammation observed in IR. The identified tissue-specific functional clusters and biomarkers might be used in a strategy for the development of tissue-targeted treatment of insulin resistance in patients.

Gene expression analysis predicts insect venom anaphylaxis in indolent systemic mastocytosis

BACKGROUND

Anaphylaxis to insect venom (Hymenoptera) is most severe in patients with mastocytosis and may even lead to death. However, not all patients with mastocytosis suffer from anaphylaxis. The aim of the study was to analyze differences in gene expression between patients with indolent systemic mastocytosis (ISM) and a history of insect venom anaphylaxis (IVA) compared to those patients without a history of anaphylaxis, and to determine the predictive use of gene expression profiling.

METHODS

Whole-genome gene expression analysis was performed in peripheral blood cells.

RESULTS

Twenty-two adults with ISM were included: 12 with a history of IVA and 10 without a history of anaphylaxis of any kind. Significant differences in single gene expression corrected for multiple testing were found for 104 transcripts (P < 0.05). Gene ontology analysis revealed that the differentially expressed genes were involved in pathways responsible for the development of cancer and focal and cell adhesion suggesting that the expression of genes related to the differentiation state of cells is higher in patients with a history of anaphylaxis. Based on the gene expression profiles, a naïve Bayes prediction model was built identifying patients with IVA.

CONCLUSIONS

In ISM, gene expression profiles are different between patients with a history of IVA and those without. These findings might reflect a more pronounced mast cells dysfunction in patients without a history of anaphylaxis. Gene expression profiling might be a useful tool to predict the risk of anaphylaxis on insect venom in patients with ISM. Prospective studies are needed to substantiate any conclusions.

SERPINE1 -675 4G/5G polymorphism is associated with asthma severity and inhaled corticosteroid response

Asthma is characterised by chronic airway inflammation and remodelling, which can be (partially) suppressed by inhaled corticosteroids (ICSs). Plasminogen activator inhibitor-1, encoded by the SERPINE1 gene, is the key inhibitor of the plasminogen activator system, which affects tissue repair and remodelling. We studied associations between a functional SERPINE1 -675 4G/5G promoter polymorphism and asthma development, severity and response to ICSs. Longitudinal cohorts of 281 asthmatics and their nonasthmatic spouses, and the general population (n=1,390) were studied. No significant associations were found with asthma development and immunoglobulin (Ig)E levels, or with forced expiratory volume in 1 s (FEV₁) in nonasthmatic controls. Asthmatic subjects carrying the SERPINE1 5G allele had higher IgE and lower lung function levels at follow-up, lower maximally attained lung function levels, and faster lung function decline compared with individuals with the 4G/4G genotype. ICS treatment showed an immediate improvement in FEV₁ in asthmatics carrying the 5G allele. However, these asthmatics still had the fastest rate of FEV₁ decline after initiating ICS treatment. Finally, the 5G allele was associated with a lower prevalence of complete asthma remission at follow-up. These findings suggest that SERPINE1 is not an asthma susceptibility gene, but rather affects the severity, progression and long-term ICS response in asthma.

Human primary adipocytes exhibit immune cell function: adipocytes prime inflammation independent of macrophages

Background

Obesity promotes inflammation in adipose tissue (AT) and this is implicated in pathophysiological complications such as insulin resistance, type 2 diabetes and cardiovascular disease. Although based on the classical hypothesis, necrotic AT adipocytes (ATA) in obese state activate AT macrophages (ATM) that then lead to a sustained chronic inflammation in AT, the link between human adipocytes and the source of inflammation in AT has not been in-depth and systematically studied. So we decided as a new hypothesis to investigate human primary adipocytes alone to see whether they are able to prime inflammation in AT.

Methods and Results

Using mRNA expression, human preadipocytes and adipocytes express the cytokines/chemokines and their receptors, MHC II molecule genes and 14 acute phase reactants including C-reactive protein. Using multiplex ELISA revealed the expression of 50 cytokine/chemokine proteins by human adipocytes. Upon lipopolysaccharide stimulation, most of these adipocyte-associated cytokines/chemokines and immune cell modulating receptors were up-regulated and a few down-regulated such as (ICAM-1, VCAM-1, MCP-1, IP-10, IL-6, IL-8, TNF-α and TNF-β highly up-regulated and IL-2, IL-7, IL-10, IL-13 and VEGF down-regulated. In migration assay, human adipocyte-derived chemokines attracted significantly more CD4+ T cells than controls and the number of migrated CD4+ cells was doubled after treating the adipocytes with LPS. Neutralizing MCP-1 effect produced by adipocytes reduced CD4+ migration by approximately 30%.

Conclusion

Human adipocytes express many cytokines/chemokines that are biologically functional. They are able to induce inflammation and activate CD4+ cells independent of macrophages. This suggests that the primary event in the sequence leading to chronic inflammation in AT is metabolic dysfunction in adipocytes, followed by production of immunological mediators by these adipocytes, which is then exacerbated by activated ATM, activation and recruitment of immune cells. This study provides novel knowledge about the prime of inflammation in human obese adipose tissue, opening a new avenue of investigations towards obesity-associated type 2 diabetes.

Changes in gene expression caused by insect venom immunotherapy responsible for the long-term protection of insect venom-allergic patients

BACKGROUND

Insect venom immunotherapy (VIT) is the only causative treatment of insect venom allergy (IVA). The immunological mechanism(s) responsible for long-term protection achieved by VIT are largely unknown. A better understanding is relevant for improving the diagnosis, prediction of anaphylaxis, and monitoring and simplifying treatment of IVA.

OBJECTIVE

To find genes that are differentially expressed during the maintenance phase of VIT and after stopping, to get clues about the pathways involved in the long-term protective effect of immunotherapy.

METHODS

Whole genome gene expression analysis was performed on RNA samples from 50 patients treated with VIT and 43 healthy controls. Patients were divided into three groups: (1) before the start of VIT; (2) on maintenance phase of VIT for at least 3 years still receiving injections; and (3) after VIT.

RESULTS

Of all 48,804 probes present in the array, 48,773 transcripts had sufficient data for further analysis. The list of genes that were differentially expressed (at least log2 FC > 2; P < .05 corrected for multiple testing) during the maintenance phase of VIT as well as after successful VIT contains 89 entities. The function of these genes affects cell signaling, cell differentiation, and ion transport.

CONCLUSION

This study shows that a group of genes is differentially expressed both during and after VIT in comparison with gene expression in patients before VIT. Although the results of this study should be confirmed prospectively, the relevance of these findings is supported by the fact that they are related to putative mechanisms of immunotherapy.

Gene expression profile, pathways, and transcriptional system regulation in indolent systemic mastocytosis

BACKGROUND

  Mastocytosis is an uncommon disease resulting from proliferation of abnormal mast cells infiltrating skin, bone marrow, liver, and other tissues. The aim of this study was to find differences in gene expression in peripheral blood cells of patients with indolent systemic mastocytosis compared to healthy controls. The second aim was to define a specific gene expression profile in patients with mastocytosis.

METHODS

  Twenty-two patients with indolent systemic mastocytosis and 43 healthy controls were studied. Whole genome gene expression analysis was performed on RNA samples isolated from the peripheral blood. For amplification and labelling of the RNA, the Illumina TotalPrep 96 RNA Amplification Kit was used. Human HT-12_V3_expression arrays were processed. Data analysis was performed using GeneSpring, Genecodis, and Transcriptional System Regulators.

RESULTS

  Comparison of gene expression between patients and controls revealed a significant difference (P < 0.05 corrected for multiple testing) and the fold change difference >2 in gene expression in 2303 of the 48.794 analysed transcripts. Functional annotation indicated that the main pathways in which the differently expressed genes were involved are ubiquitin-mediated proteolysis, MAPK signalling pathway, pathways in cancer, and Jak-STAT signalling. The expression distributions for both groups did not overlap at all, indicating that many genes are highly differentially expressed in both groups.

CONCLUSION

  We were able to find abnormalities in gene expression in peripheral blood cells of patients with indolent systemic mastocytosis and to construct a gene expression profile which may be useful in clinical practice to predict the presence of mastocytosis and in further research of novel drugs.

Gene expression studies in cells from primary ciliary dyskinesia patients identify 208 potential ciliary genes

Cilia are small cellular projections that either act as sensors (primary cilia) or propel fluid over the epithelia of various organs (motile cilia). The organellum has gained much attention lately because of its involvement in a group of human diseases called ciliopathies. Primary ciliary dyskinesia (PCD) is an autosomal recessive ciliopathy caused by mutations in cilia motility genes. The disease is characterized by recurrent respiratory tract infections due to the lack of an efficient mucociliary clearance. We performed whole-genome gene expression profiling in bronchial biopsies from PCD patients. We used the quality threshold clustering algorithm to identify groups of genes that revealed highly correlated RNA expression patterns in the biopsies. The largest cluster contained 372 genes and was significantly enriched for genes related to cilia. The database and literature search showed that 164 genes in this cluster were known cilia genes, strongly indicating that the remaining 208 genes were likely to be new cilia genes. The tissue expression pattern of the 208 new cilia genes and the 164 known genes was consistent with the presence of motile cilia in a given tissue. The analysis of the upstream promotor sequences revealed evidence for RFX transcription factors binding site motif in both subgroups. Based on the correlated expression patterns in PCD-affected tissues, we identified 208 genes that we predict to be involved in cilia biology. Our predictions are based directly on the human material and not on model organisms. This list of genes provides candidate genes for PCD and other ciliopathies.

Co-expressed immune and metabolic genes in visceral and subcutaneous adipose tissue from severely obese individuals are associated with plasma HDL and glucose levels: a microarray study

Background

Excessive accumulation of body fat, in particular in the visceral fat depot, is a major risk factor to develop a variety of diseases such as type 2 diabetes. The mechanisms underlying the increased risk of obese individuals to develop co-morbid diseases are largely unclear.

We aimed to identify genes expressed in subcutaneous adipose tissue (SAT) and visceral adipose tissue (VAT) that are related to blood parameters involved in obesity co-morbidity, such as plasma lipid and glucose levels, and to compare gene expression between the fat depots.

Methods

Whole-transcriptome SAT and VAT gene expression levels were determined in 75 individuals with a BMI >35 kg/m². Modules of co-expressed genes likely to be functionally related were identified and correlated with BMI, plasma levels of glucose, insulin, HbA_1c, triglycerides, non-esterified fatty acids, ALAT, ASAT, C-reactive protein, and LDL- and HDL cholesterol.

Results

Of the approximately 70 modules identified in SAT and VAT, three SAT modules were inversely associated with plasma HDL-cholesterol levels, and a fourth module was inversely associated with both plasma glucose and plasma triglyceride levels (p < 5.33 × 10^-5). These modules were markedly enriched in immune and metabolic genes. In VAT, one module was associated with both BMI and insulin, and another with plasma glucose (p < 4.64 × 10^-5). This module was also enriched in inflammatory genes and showed a marked overlap in gene content with the SAT modules related to HDL. Several genes differentially expressed in SAT and VAT were identified.

Conclusions

In obese subjects, groups of co-expressed genes were identified that correlated with lipid and glucose metabolism parameters; they were enriched with immune genes. A number of genes were identified of which the expression in SAT correlated with plasma HDL cholesterol, while their expression in VAT correlated with plasma glucose. This underlines both the singular importance of these genes for lipid and glucose metabolism and the specific roles of these two fat depots in this respect.

Identification of PCDH1 as a novel susceptibility gene for bronchial hyperresponsiveness

RATIONALE

Asthma is a chronic inflammatory airway disease that affects more than 300 million individuals worldwide. Asthma is caused by interaction of genetic and environmental factors. Bronchial hyperresponsiveness (BHR) is a hallmark of asthma and results from increased sensitivity of the airways to physical or chemical stimulants. BHR and asthma are linked to chromosome 5q31-q33.

OBJECTIVES

To identify a gene for BHR on chromosome 5q31-q33.

METHODS

In 200 Dutch families with asthma, linkage analysis and fine mapping were performed, and the Protocadherin 1 gene (PCDH1) was identified. PCDH1 was resequenced in 96 subjects from ethnically diverse populations to identify novel sequence variants. Subsequent replication studies were undertaken in seven populations from The Netherlands, the United Kingdom, and the United States, including two general population samples, two family samples, and three case-control samples. PCDH1 mRNA and protein expression was investigated using polymerase chain reaction, Western blotting, and immunohistochemistry.

MEASUREMENTS AND MAIN RESULTS

In seven out of eight populations (n = 6,168) from The Netherlands, United Kingdom, and United States, PCHD1 gene variants were significantly associated with BHR (P values, 0.005-0.05) This association was present in both families with asthma and general populations. PCDH1 mRNA and protein were expressed in airway epithelial cells and in macrophages.

CONCLUSIONS

PCDH1 is a novel gene for BHR in adults and children. The identification of PCDH1 as a BHR susceptibility gene may suggest that a structural defect in the integrity of the airway epithelium, the first line of defense against inhaled substances, contributes to the development of BHR.

Coeliac disease-associated risk variants in TNFAIP3 and REL implicate altered NF-kappaB signalling

OBJECTIVE

Our previous coeliac disease genome-wide association study (GWAS) implicated risk variants in the human leucocyte antigen (HLA) region and eight novel risk regions. To identify more coeliac disease loci, we selected 458 single nucleotide polymorphisms (SNPs) that showed more modest association in the GWAS for genotyping and analysis in four independent cohorts.

DESIGN

458 SNPs were assayed in 1682 cases and 3258 controls from three populations (UK, Irish and Dutch). We combined the results with the original GWAS cohort (767 UK cases and 1422 controls); six SNPs showed association with p<1 x 10(-04) and were then genotyped in an independent Italian coeliac cohort (538 cases and 593 controls).

RESULTS

We identified two novel coeliac disease risk regions: 6q23.3 (OLIG3-TNFAIP3) and 2p16.1 (REL), both of which reached genome-wide significance in the combined analysis of all 2987 cases and 5273 controls (rs2327832 p = 1.3 x 10(-08), and rs842647 p = 5.2 x 10(-07)). We investigated the expression of these genes in the RNA isolated from biopsies and from whole blood RNA. We did not observe any changes in gene expression, nor in the correlation of genotype with gene expression.

CONCLUSIONS

Both TNFAIP3 (A20, at the protein level) and REL are key mediators in the nuclear factor kappa B (NF-kappaB) inflammatory signalling pathway. For the first time, a role for primary heritable variation in this important biological pathway predisposing to coeliac disease has been identified. Currently, the HLA risk factors and the 10 established non-HLA risk factors explain approximately 40% of the heritability of coeliac disease.

Complex nature of SNP genotype effects on gene expression in primary human leucocytes

Background

Genome wide association studies have been hugely successful in identifying disease risk variants, yet most variants do not lead to coding changes and how variants influence biological function is usually unknown.

Methods

We correlated gene expression and genetic variation in untouched primary leucocytes (n = 110) from individuals with celiac disease – a common condition with multiple risk variants identified. We compared our observations with an EBV-transformed HapMap B cell line dataset (n = 90), and performed a meta-analysis to increase power to detect non-tissue specific effects.

Results

In celiac peripheral blood, 2,315 SNP variants influenced gene expression at 765 different transcripts (< 250 kb from SNP, at FDR = 0.05, cis expression quantitative trait loci, eQTLs). 135 of the detected SNP-probe effects (reflecting 51 unique probes) were also detected in a HapMap B cell line published dataset, all with effects in the same allelic direction. Overall gene expression differences within the two datasets predominantly explain the limited overlap in observed cis-eQTLs. Celiac associated risk variants from two regions, containing genes IL18RAP and CCR3, showed significant cis genotype-expression correlations in the peripheral blood but not in the B cell line datasets. We identified 14 genes where a SNP affected the expression of different probes within the same gene, but in opposite allelic directions. By incorporating genetic variation in co-expression analyses, functional relationships between genes can be more significantly detected.

Conclusion

In conclusion, the complex nature of genotypic effects in human populations makes the use of a relevant tissue, large datasets, and analysis of different exons essential to enable the identification of the function for many genetic risk variants in common diseases.

Restenosis after percutaneous coronary intervention is associated with the angiotensin-II type-1 receptor 1166A/C polymorphism but not with polymorphisms of angiotensin-converting enzyme, angiotensin-II receptor, angiotensinogen or heme oxygenase-1

OBJECTIVES

The renin-angiotensin system (RAS) is thought to play a major role in the pathophysiology of de-novo restenotic lesions and in-stent restenosis after percutaneous coronary intervention (PCI). Heme oxygenase-1 (HO-1), is thought to beneficially influence these processes. We examined the effect of pharmacologic as well as genetic RAS interactions on restenosis in a large population of consecutive patients undergoing PCI, and evaluated possible gene-gene interactions in both systems.

METHODS

The GENDER project is a multicenter prospective follow-up study, including 3146 patients after successful PCI. Genotyping in these patients was performed for the ACE gene insertion/deletion, the angiotensinogen 235Met/Thr, T174M and A(-6)G, the angiotensin-II type 1 receptor (AT1R) 1166A/C and T810A, the angiotensin-II type 2 receptor (AT2R) 1675G/A and 3123A polymorphisms and the length polymorphism in the HO-1 promoter region.

RESULTS

A total of 3104 patients were followed for 10 months. In 2975 patients at least one of the nine genotypes could be determined. The AT1R 1166 CC genotype showed a significant association with TVR; the other polymorphisms did not. RAS-inhibitory drugs were not associated with the incidence of TVR, nor did they interact with any of the investigated polymorphisms. Patients with the ACE I/I polymorphism showed a trend towards a better outcome if they had a short number of repeats in the HO-1 promoter. This relationship was inversely present in carriers of the ACE D/D polymorphism.

CONCLUSION

We could only establish a role for the AT1R 1166A/C polymorphism in restenosis after PCI. However, significant gene-gene interaction was suggested for the ACE gene and the HO-1 promotor. The RAS and HO-1 relation in restenosis merits further investigation.

Decorin and TGF-beta1 polymorphisms and development of COPD in a general population

Background

Decorin, an extracellular matrix (ECM) proteoglycan, and TGF-β₁ are both involved in lung ECM turnover. Decorin and TGF-β₁ expression are decreased respectively increased in COPD lung tissue. Interestingly, they act as each other's feedback regulator. We investigated whether single nucleotide polymorphisms (SNPs) in decorin and TGF-β₁ underlie accelerated decline in FEV₁ and development of COPD in the general population.

Methods

We genotyped 1390 subjects from the Vlagtwedde/Vlaardingen cohort. Lung function was measured every 3 years for a period of 25 years. We tested whether five SNPs in decorin (3'UTR and four intron SNPs) and three SNPs in TGF-β₁ (3'UTR rs6957, C-509T rs1800469 and Leu10Pro rs1982073), and their haplotypes, were associated with COPD (last survey GOLD stage = II). Linear mixed effects models were used to analyze genotype associations with FEV₁ decline.

Results

We found a significantly higher prevalence of carriers of the minor allele of the TGF-β₁ rs6957 SNP (p = 0.001) in subjects with COPD. Additionally, we found a significantly lower prevalence of the haplotype with the major allele of rs6957 and minor alleles for rs1800469 and rs1982073 SNPs in TGF-β₁ in subjects with COPD (p = 0.030), indicating that this association is due to the rs6957 SNP. TGF-β₁ SNPs were not associated with FEV₁ decline. SNPs in decorin, and haplotypes constructed of both TGF-β₁ and decorin SNPs were not associated with development of COPD or with FEV₁ decline.

Conclusion

Our study shows for the first time that SNPs in decorin on its own or in interaction with SNPs in TGF-β₁ do not underlie the disturbed balance in expression between these genes in COPD. TGF-β₁ SNPs are associated with COPD, yet not with accelerated FEV₁ decline in the general population.

No increased susceptibility to breast cancer from combined CHEK2 1100delC genotype and the HLA class III region risk factors

CHEK2 is low-penetrance breast cancer susceptibility gene. The 1100delC mutation may interact with variants/mutations in other breast cancer susceptibility loci. We identified a risk haplotype in the HLA class III region in breast cancer patients [de Jong MM, Nolte IM, de Vries EGE, et al. The HLA class III subregion is responsible for an increased breast cancer risk. Hum Mol Genet 2003, 12, 2311-2319] and tested whether it interacted with 1100delC mutation. The CHEK2 1100delC mutation was analysed in the same series of patients and controls as in the HLA breast cancer study. In 962 unselected breast cancer patients, the 1100delC mutation was observed in 2.9% and in 367 controls in 1.4% (NS). The highest 1100delC frequency occurred in high-risk (4.4%), followed by moderate-risk (3.8%), and lowest in low genetic risk patients (2.4%, P(trend) 0.029). In HLA risk haplotype carriers no increased breast cancer risk was observed in the presence of 1100delC mutation. Patients more often had one than both genetic risk factors. The 1100delC mutation and the HLA risk haplotype confer increased breast cancer risks, but an interactive effect on breast cancer between both factors is unlikely. In contrast, the effect of 1100delC mutation on breast cancer risk was limited to individuals without HLA risk haplotype, suggesting a mutual excluding effect between these risk factors.

A disintegrin and metalloprotease 33 polymorphisms and lung function decline in the general population

RATIONALE

A disintegrin and metalloprotease 33 (ADAM33) has been identified as a susceptibility gene for asthma and single nucleotide polymorphisms (SNPs) in this gene have been associated with excessive decline of lung function in individuals with asthma.

OBJECTIVES

To assess whether SNPs in ADAM33 are associated with accelerated lung function loss in the general population and with chronic obstructive pulmonary disease (COPD).

METHODS

DNA was collected from subjects of the Vlagtwedde-Vlaardingen cohort participating in the last survey in 1989-1990 after a follow-up of 25 years. Information was collected every 3 years, including lung function measurements. We defined COPD as GOLD stage 2 or higher at the last survey. A total of 1,390 subjects from the cohort was genotyped for the following SNPs in ADAM33: F+1, Q-1, S_1, S_2, T_1, T_2, V_4, and ST+5. Differences in prevalence of SNPs were analyzed with chi(2) tests. Linear mixed effects models were used to analyze FEV(1) decline according to genotype.

MEASUREMENTS AND MAIN RESULTS

In the whole population, mean adjusted decline was 18.7 and 12.7 ml/year in females and males, respectively. Individuals homozygous for minor alleles of SNPs S_2 and Q-1 and heterozygous for SNP S_1 had a significantly accelerated decline in FEV(1) of, respectively, 4.9, 9.6, and 3.6 ml/year compared with wild type. We found a significantly higher prevalence of SNPs F+1, S_1, S_2, and T_2 in subjects with COPD.

CONCLUSIONS

We demonstrated that SNPs in ADAM33 are associated with accelerated lung function decline in the general population. These SNPs are also risk factors for COPD.

Polymorphisms in SPINK5 are not associated with asthma in a Dutch population

BACKGROUND

Asthma and allergic phenotypes are complex genetic diseases with known linkage to chromosome 5q. This region has many candidate genes, including serine protease inhibitor Kazal type 5 (SPINK5), which has been associated with asthma and atopic dermatitis in family-based studies of children with atopic dermatitis.

OBJECTIVE

We sought to investigate whether single nucleotide polymorphisms in SPINK5 are associated with asthma, atopic phenotypes, and atopic dermatitis.

METHODS

We investigated whether single nucleotide polymorphisms in SPINK5 (ie, -785 A/G, Asn368Ser, and Lys420Glu) are associated with asthma, atopic phenotypes, and atopic dermatitis in 200 families ascertained by a proband with asthma (nonaffected spouses served as a matched control population) and an independent set of 252 trios with asthma.

RESULTS

We found no association with asthma, atopic phenotypes, and atopic dermatitis after correction for multiple testing.

CONCLUSION

The negative results in this study suggest that SPINK5 is not associated with asthma or atopic phenotypes in individuals ascertained by a proband with asthma. This is consistent with the finding that SPINK5 is not expressed in the lung. Because our patients were ascertained for asthma, a role of SPINK5 in atopic dermatitis cannot be excluded.

Haplotype sharing analysis identifies a retroviral dUTPase as candidate susceptibility gene for psoriasis

The psoriasis susceptibility locus 1 (PSORS1) mutation is assumed to reside within a region around human leukocyte antigen-C spanning 250 kb, termed risk haplotype (RH) 1/2. By re-analyzing a published data set with a previously developed method, the haplotype sharing statistic, we confirm localization of PSORS1 to the RH1 region and refine its location to marker M6S168. We replicate this result in an independent patient sample. The target region harbors fragments of a human endogenous retrovirus K (HERV-K) endogenous retrovirus. Two single-nucleotide polymorphisms with alleles differing between high- and low-risk haplotypes are located within the HERV-K dUTPase. One of these encodes a predicted non-conserved Glu-Arg exchange. The HERV-K dUTPase is expressed in peripheral blood and in normal as well as lesional psoriatic skin. Our results indicate that an endogenous retroviral dUTPase constitutes a candidate gene for the PSORS1 mutation.

The human leukocyte antigen region and colorectal cancer risk

PURPOSE

Recently, we found a certain haplotype in the human leukocyte antigen Class III subregion to be associated with breast cancer. Epidemiologic studies have shown that breast cancer and colorectal cancer have several risk factors in common. In view of these studies and because polymorphisms located in the human leukocyte antigen III region have been found to be associated with colorectal cancer, we wondered whether the same region also is involved in colorectal cancer susceptibility.

METHODS

The human leukocyte antigen region was genotyped with 14 microsatellite markers in germline DNA from 643 colorectal cancer patients and 841 family-based controls. Association analyses and the Haplotype Sharing Statistic were used to search for differences between patients and controls. Subgroup analyses were performed for gender, age at diagnosis, and localization of the tumor.

RESULTS

The Haplotype Sharing Statistic analysis revealed neither a difference in mean haplotype sharing between all patients and controls, nor in any of the subgroups. The single allele, genotype, and two-locus association analyses for all patients and for the different subgroups did not show an association with colorectal cancer for the 14 microsatellite markers. Also, no association was observed between the tumor necrosis factor-beta polymorphism and colorectal cancer.

CONCLUSIONS

No association was observed between commonly occurring haplotypes and alleles in the human leukocyte antigen region and colorectal cancer risk.

Efficiency Control in Large-Scale Genotyping Using Analysis of Variance

The efficiency of the genotyping process is determined by many simultaneous factors. In actual genotyping, a production run is often preceded by small-scale experiments to find optimal conditions. We propose to use statistical analysis of production run data as well, to gain insight into factors important for the outcome of genotyping. As an example, we show that analysis of variance (ANOVA) applied to the first-pass results of a genetic study reveals important determinants of genotyping success. The largest factor limiting genotyping appeared to be interindividual variation among DNA samples, explaining 20% of the variance, and a smaller reaction volume, sizing failure, and differences among markers all explained approximately 10%. Other potentially important factors, such as sample position within the plate and reusing electrophoresis matrix, appeared to be of minor influence. About 55% of the total variance could be explained by systematic factors. These results show that ANOVA can provide valuable feedback to improve genotyping efficiency. We propose to adjust genotype production runs using principles of experimental design in order to maximize genotyping efficiency at little additional cost.

Colorectal cancer and theCHEK2 1100delC mutation

The CHEK2 1100delC mutation was recently identified as a low-penetrance breast cancer susceptibility allele. The mutation occurred more frequently in families with clustering of breast and colorectal cancers (CRCs) than in families with clustering of breast cancer only. Hence, the 1100delC mutation could also be a low-penetrance CRC susceptibility allele. To test this hypothesis, we examined the mutation in 629 unselected CRC cases, 230 controls, and 105 selected CRCs diagnosed in patients before age 50. The mutation was observed in 1.6% of unselected patients and in 0.3% of controls (Not significant (NS)). After stratifying unselected patients according to defined genetic risk (on the basis of age at diagnosis and family history of colorectal and endometrial cancer), the highest frequency was observed in high-risk patients (12.5%), followed by moderate-risk patients (3.3%), and was lowest in low-risk patients (1.0%, P(trend) 0.014). In selected patients, 1.6% carried the mutation (NS). Subgroup analyses for tumor localization, gender, and age at diagnosis did not reveal an association with the 1100delC genotype. In addition, a pooled analysis, combining data of one published study in unselected CRC cases and our study, also did not reveal an association. In conclusion, the frequency of the 1100delC genotype was neither significantly increased in unselected CRC patients nor in selected CRC patients diagnosed before age 50. However, after stratifying unselected CRC patients according to defined genetic risk, a significant trend of increasing frequency was observed. Together, the results are consistent with a low-penetrance effect (OR 1.5-2.0) of the CHEK2 1100delC on CRC risk. Large case-control studies are required to clarify the exact role of the CHEK2 1100delC mutation in CRC.

No association between the Arg201Gly polymorphism of the DCC gene and colorectal cancer

BACKGROUND AND AIMS

In one small study, the DCC Arg201Gly polymorphism has been observed more frequently in colorectal cancer cases compared with controls. We wondered whether these results could be replicated in a much larger study.

METHODOLOGY

The DCC Arg201 Gly polymorphism was genotyped in 625 unselected Caucasian colorectal cancer patients and 220 controls. Association analysis was used to search for a difference between patients and controls. Subgroup analyses were performed for site of tumour, gender, age at diagnosis, family history of colorectal cancer and modified Dukes classification.

RESULTS

The association analyses revealed no difference in Arg201Gly genotype frequency between patients and controls, neither overall nor for different subgroups according to site of tumour, gender, age at diagnosis, family history of colorectal cancer and modified Dukes classification.

CONCLUSION

No association was observed between the Arg201Gly polymorphism of DCC and colorectal cancer risk.

Evaluation of the IRF-2 Gene as a Candidate for PSORS3

Type 1 interferon can trigger flares of psoriasis. Hypersensitivity to type 1 interferon signaling causes a psoriasis-like skin disease in mice deficient for the transcription factor interferon regulatory factor 2 (IRF2). The human IRF2 gene is located at a previously identified candidate psoriasis susceptibility locus on chromosome 4q (PSORS3 at D4S1535). Therefore, we tested association of psoriasis with IRF2. We generated a sample consisting of 157 families with a total of 521 individuals. Five novel microsatellite markers were developed and typed, and complemented with three known markers to yield a set of eight markers spaced within 600 kb around the IRF2 gene, three of which are located in the gene. We detected association of IRF2 with type 1 psoriasis at two markers in the IRF2 gene. Haplotype sharing analysis confirmed association of IRF2 with type 1 psoriasis (p=0.0017; pcorr=0.03). The 921G/A SNP in exon 9 was found to obliterate a predicted exon splice enhancer in an allele-specific manner. There was a suggestive increase of homozygosity for the splicing-deficient allele in type 1 psoriasis patients. Our data identify IRF2 as a potential susceptibility gene for psoriasis.

The HLA class III subregion is responsible for an increased breast cancer risk

BRCA1 and BRCA2 germline mutations account for <5% of breast cancer cases. Less penetrant breast cancer susceptibility genes are likely to exist. Earlier studies have suggested involvement of the HLA region. The HLA region was genotyped with 24 microsatellite markers and markers for two single nucleotide polymorphisms (SNPs) in TNFalpha and TNFbeta, in germline DNA from 956 breast cancer patients and 1271 family-based controls. Association analyses and the haplotype sharing statistic (HSS) were used to search for differences in haplotype sharing between patients and controls. Based on criteria known to influence genetic breast cancer risk, patients were divided into groups of high, moderate and low risk. The HSS revealed a significant difference in mean haplotype sharing between patients and controls for four consecutive markers (D6S2671, TNFa, D6S2672 and MICA), the highest being at D6S2671 (P=0.017). Subgroup analyses showed that moderate-risk patients were responsible for this difference, with the strongest association for D6S2672 (P=0.0009). A single haplotype was more frequent and longer in moderate-risk patients than in controls. The results were confirmed with association analyses. Individuals homozygous for haplotype 110-184 (D6S2672-MICA) were observed in 9.0% of moderate-risk patients and 1.5% of controls [odds ratio (OR)=7.14], while heterozygotes were at a lower risk (OR=1.41), suggesting a recessive effect. No association was observed between the two SNPs in TNFalpha (-308) and TNFbeta (intron 1) and breast cancer risk. The results reveal a potential role of the HLA class III subregion in susceptibility to breast cancer in patients at moderate familial risk.

Testicular carcinoma and HLA Class II genes

BACKGROUND

The association with histocompatibility antigens (HLA), in particular Class II genes (DQB1, DRB1), has recently been suggested to be one of the genetic factors involved in testicular germ cell tumor (TGCT) development. The current study, which uses genotyping of microsatellite markers, was designed to replicate previous associations.

METHODS

In 151 patients, along with controls comprising parents or spouses, the HLA region (particularly Class II) on chromosome 6p21 was genotyped for a set of 15 closely linked microsatellite markers.

RESULTS

In both patients and controls, strong linkage disequilibrium was observed in the genotyped region, indicating that similar haplotypes are likely to be identical by descent. However, association analysis and the transmission disequilibrium test did not show significant results. Haplotype sharing statistics, a haplotype method that derives extra information from phase and single marker tests, did not show differences in haplotype sharing between patients and controls.

CONCLUSION

The current genotyping study did not confirm the previously reported association between HLA Class II genes and TGCT. As the HLA alleles for which associations were reported are also prevalent in the Dutch populations, these associations are likely to be nonexistent or much weaker than previously reported.

Mapping of a susceptibility gene for multiple sclerosis to the 51 kb interval between G511525 and D6S1666 using a new method of haplotype sharing analysis

Multiple sclerosis (MS) is a complex disease that is partly genetic in origin. Although an association of MS with specific human leukocyte antigen (HLA) types has been known for almost 30 years, the nature of this relationship has remained unclear. Furthermore, genetic resolution sufficient to implicate a specific gene in the HLA region has not been achieved. Many loci in the HLA region have been found to be significantly associated with MS, which is largely explained by the extended haplotype sharing and varying marker informativity of the region. We have determined 248 haplotypes of MS patients from the population of the northern Netherlands and 226 haplotypes of their relatives as controls using a set of 22 microsatellite markers covering the HLA region. The data were analyzed using standard association methods and a new statistical method, haplotype sharing statistics (HSS). Haplotype sharing statistics determines the extent of haplotype sharing for all pairs of haplotypes of patients and of controls and calculates the difference in mean haplotype sharing between patients and controls. Haplotype sharing was found to be significantly greater among patients than among controls in a region of 1.1 Mb between markers G511525 and TNFalpha. The involvement of this region is also supported by association analysis and the transmission/disequilibrium test (TDT). Within this region, HSS, which is largely independent of association and TDT, indicated the interval of 51 kb between G511525 and D6S1666 as that most likely to contain a susceptibility gene for MS. As DQB1 is the sole gene known in this interval at present, the results of our analysis suggest that this gene plays a role in the pathogenesis of MS.