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.