A retrospective comparative exploratory study on two methylentetrahydrofolate reductase (MTHFR) polymorphisms in esophagogastric cancer: the A1298C MTHFR polymorphism is an independent prognostic factor only in neoadjuvantly treated gastric cancer patients.

BACKGROUND Methylentetrahydrofolate reductase (MTHFR) plays a major role in folate metabolism and consequently could be an important factor for the efficacy of a treatment with 5-fluorouracil. Our aim was to evaluate the prognostic and predictive value of two well characterized constitutional MTHFR gene polymorphisms for primarily resected and neoadjuvantly treated esophagogastric adenocarcinomas. METHODS 569 patients from two centers were analyzed (gastric cancer: 218, carcinoma of the esophagogastric junction (AEG II, III): 208 and esophagus (AEG I): 143). 369 patients received neoadjuvant chemotherapy followed by surgery, 200 patients were resected without preoperative treatment. The MTHFR C677T and A1298C polymorphisms were determined in DNA from peripheral blood lymphozytes. Associations with prognosis, response and clinicopathological factors were analyzed retrospectively within a prospective database (chi-square, log-rank, cox regression). RESULTS Only the MTHFR A1298C polymorphisms had prognostic relevance in neoadjuvantly treated patients but it was not a predictor for response to neoadjuvant chemotherapy. The AC genotype of the MTHFR A1298C polymorphisms was significantly associated with worse outcome (p = 0.02, HR 1.47 (1.06-2.04). If neoadjuvantly treated patients were analyzed based on their tumor localization, the AC genotype of the MTHFR A1298C polymorphisms was a significant negative prognostic factor in patients with gastric cancer according to UICC 6th edition (gastric cancer including AEG type II, III: HR 2.0, 95% CI 1.3-2.0, p = 0.001) and 7th edition (gastric cancer without AEG II, III: HR 2.8, 95% CI 1.5-5.7, p = 0.003), not for AEG I. For both definitions of gastric cancer the AC genotype was confirmed as an independent negative prognostic factor in cox regression analysis. In primarily resected patients neither the MTHFR A1298C nor the MTHFR C677T polymorphisms had prognostic impact. CONCLUSIONS The MTHFR A1298C polymorphisms was an independent prognostic factor in patients with neoadjuvantly treated gastric adenocarcinomas (according to both UICC 6th or 7th definitions for gastric cancer) but not in AEG I nor in primarily resected patients, which confirms the impact of this enzyme on chemotherapy associated outcome.

PTX3 Polymorphisms and Invasive Mold Infections after Solid Organ Transplantation

Donor PTX3 polymorphisms were shown to influence the risk of invasive aspergillosis among hematopoietic stem cell transplant recipients. Here, we show that PTX3 polymorphisms are independent risk factors for invasive mold infection among 1101 solid organ transplant recipients, thereby strengthening their role in mold infection pathogenesis and patient's risk stratification.

Polymorphisms in MIR27A Associated with Early-Onset Toxicity in Fluoropyrimidine-Based Chemotherapy

PURPOSE The microRNA miR-27a was recently shown to directly regulate dihydropyrimidine dehydrogenase (DPD), the key enzyme in fluoropyrimidine catabolism. A common polymorphism (rs895819A>G) in the miR-27a genomic region (MIR27A) was associated with reduced DPD activity in healthy volunteers, but the clinical relevance of this effect is still unknown. Here, we assessed the association of MIR27A germline variants with early-onset fluoropyrimidine toxicity. EXPERIMENTAL DESIGN MIR27A was sequenced in 514 patients with cancer receiving fluoropyrimidine-based chemotherapy. Associations of MIR27A polymorphisms with early-onset (cycles 1-2) fluoropyrimidine toxicity were assessed in the context of known risk variants in the DPD gene (DPYD) and additional covariates associated with toxicity. RESULTS The association of rs895819A>G with early-onset fluoropyrimidine toxicity was strongly dependent on DPYD risk variant carrier status (Pinteraction = 0.0025). In patients carrying DPYD risk variants, rs895819G was associated with a strongly increased toxicity risk [OR, 7.6; 95% confidence interval (CI), 1.7-34.7; P = 0.0085]. Overall, 71% (12/17) of patients who carried both rs895819G and a DPYD risk variant experienced severe toxicity. In patients without DPYD risk variants, rs895819G was associated with a modest decrease in toxicity risk (OR, 0.62; 95% CI, 0.43-0.9; P = 0.012). CONCLUSIONS These results indicate that miR-27a and rs895819A>G may be clinically relevant for further toxicity risk stratification in carriers of DPYD risk variants. Our data suggest that direct suppression of DPD by miR-27a is primarily relevant in the context of fluoropyrimidine toxicity in patients with reduced DPD activity. However, miR-27a regulation of additional targets may outweigh its effect on DPD in patients without DPYD risk variants.

Bovine prion protein gene (PRNP) promoter polymorphisms modulate PRNP expression and may be responsible for differences in bovine spongiform encephalopathy susceptibility

The susceptibility of humans to the variant Creutzfeldt-Jakob disease is greatly influenced by polymorphisms within the human prion protein gene (PRNP). Similar genetic differences exist in sheep, in which PRNP polymorphisms modify the susceptibility to scrapie. However, the known coding polymorphisms within the bovine PRNP gene have little or no effect on bovine spongiform encephalopathy (BSE) susceptibility in cattle. We have recently found a tentative association between PRNP promoter polymorphisms and BSE susceptibility in German cattle (Sander, P., Hamann, H., Pfeiffer, I., Wemheuer, W., Brenig, B., Groschup, M., Ziegler, U., Distl, O., and Leeb, T. (2004) Neurogenetics 5, 19-25). A plausible hypothesis explaining this observation could be that the bovine PRNP promoter polymorphisms cause changes in PRNP expression that might be responsible for differences in BSE incubation time and/or BSE susceptibility. To test this hypothesis, we performed a functional promoter analysis of the different bovine PRNP promoter alleles by reporter gene assays in vitro and by measuring PRNP mRNA levels in calves with different PRNP genotypes in vivo. Two variable sites, a 23-bp insertion/deletion (indel) polymorphism containing a RP58-binding site and a 12-bp indel polymorphism containing an SP1-binding site, were investigated. Band shift assays indicated differences in transcription factor binding to the different alleles at the two polymorphisms. Reporter gene assays demonstrated an interaction between the two postulated transcription factors and lower expression levels of the ins/ins allele compared with the del/del allele. The in vivo data revealed substantial individual variation of PRNP expression in different tissues. In intestinal lymph nodes, expression levels differed between the different PRNP genotypes.

Polymorphisms within the canine MLPH gene are associated with dilute coat color in dogs

BACKGROUND Pinschers and other dogs with coat color dilution show a characteristic pigmentation phenotype. The fur colors are a lighter shade, e.g. silvery grey (blue) instead of black and a sandy color (Isabella fawn) instead of red or brown. In some dogs the coat color dilution is sometimes accompanied by hair loss and recurrent skin inflammation, the so called color dilution alopecia (CDA) or black hair follicular dysplasia (BHFD). In humans and mice a comparable pigmentation phenotype without any documented hair loss is caused by mutations within the melanophilin gene (MLPH). RESULTS We sequenced the canine MLPH gene and performed a mutation analysis of the MLPH exons in 6 Doberman Pinschers and 5 German Pinschers. A total of 48 sequence variations was identified within and between the breeds. Three families of dogs showed co-segregation for at least one polymorphism in an MLPH exon and the dilute phenotype. No single polymorphism was identified in the coding sequences or at splice sites that is likely to be causative for the dilute phenotype of all dogs examined. In 18 German Pinschers a mutation in exon 7 (R199H) was consistently associated with the dilute phenotype. However, as this mutation was present in homozygous state in four dogs of other breeds with wildtype pigmentation, it seems unlikely that this mutation is truly causative for coat color dilution. In Doberman Pinschers as well as in Large Munsterlanders with BHFD, a set of single nucleotide polymorphisms (SNPs) around exon 2 was identified that show a highly significant association to the dilute phenotype. CONCLUSION This study provides evidence that coat color dilution is caused by one or more mutations within or near the MLPH gene in several dog breeds. The data on polymorphisms that are strongly associated with the dilute phenotype will allow the genetic testing of Pinschers to facilitate the breeding of dogs with defined coat colors and to select against Large Munsterlanders carrying BHFD.

Can ancient colour polymorphisms explain why some cichlid lineages speciate rapidly under disruptive sexual selection?

It is not sufficiently understood why some lineages of cichlid fishes have proliferated in the Great Lakes of East Africa much more than anywhere else in the world, and much faster than other cichlid lineages or any other group of freshwater fish. Recent field and experimental work on Lake Victoria haplochromines suggests that mate choice-mediated disruptive sexual selection on coloration, that can cause speciation even in the absence of geographical isolation, may explain it. We summarize the evidence and propose a hypothesis for the genetics of coloration that may help understand the phenomenon. By detl ning colour patterns by hue and arrangement of hues on the body, we could assign almost all observed phenotypes of Lake Victoria cichlids to one of three female («plain», «orange blotched», «black and white») and three male («blue», «red-ventrum», «reddorsum») colour patterns. These patterns diagnose species but frequently eo-occur also as morphs within the same population, where they are associated with variation in mate preferences, and appear to be transient stages in speciation. Particularly the male patterns occur in almost every genus of the species flock. We propose that the patterns and their association into polymorphisms express an ancestral trait that is retained across speciation. Our model for male colour pattern assumes two structural loci. When both are switched off, the body is blue. When switched on by a cascade of polymorphic regulatory genes, one expresses a yellow to red ventrum, the other one a yellow to red dorsum. The expression of colour variation initiates speciation. The blue daughter species will inherit the variation at the regulatory genes that can, without new mutational events, purely by recombination, again expose the colour polymorphism, starting the process anew. Very similar colour patterns also dominate among the Mbuna of Lake Malawi. In contrast, similar colour polymorphisms do not exist in the lineages that have not proliferated in the Great Lakes. The colour pattern polymorphism may be an ancient trait in the lineage (or lineages) that gave rise to the two large haplochromine radiations. We propose two tests of our hypothesis.

Genetic polymorphisms associated with the inflammatory response in bacterial meningitis.

BACKGROUND Bacterial meningitis (BM) is an infectious disease that results in high mortality and morbidity. Despite efficacious antibiotic therapy, neurological sequelae are often observed in patients after disease. Currently, the main challenge in BM treatment is to develop adjuvant therapies that reduce the occurrence of sequelae. In recent papers published by our group, we described the associations between the single nucleotide polymorphisms (SNPs) AADAT +401C > T, APEX1 Asn148Glu, OGG1 Ser326Cys and PARP1 Val762Ala and BM. In this study, we analyzed the associations between the SNPs TNF -308G > A, TNF -857C > T, IL-8 -251A > T and BM and investigated gene-gene interactions, including the SNPs that we published previously. METHODS The study was conducted with 54 BM patients and 110 healthy volunteers (as the control group). The genotypes were investigated via primer-introduced restriction analysis-polymerase chain reaction (PIRA-PCR) or polymerase chain reaction-based restriction fragment length polymorphism (PCR-RFLP) analysis. Allelic and genotypic frequencies were also associated with cytokine and chemokine levels, as measured with the x-MAP method, and cell counts. We analyzed gene-gene interactions among SNPs using the generalized multifactor dimensionality reduction (GMDR) method. RESULTS We did not find significant association between the SNPs TNF -857C > T and IL-8 -251A > T and the disease. However, a higher frequency of the variant allele TNF -308A was observed in the control group, associated with changes in cytokine levels compared to individuals with wild type genotypes, suggesting a possible protective role. In addition, combined inter-gene interaction analysis indicated a significant association between certain genotypes and BM, mainly involving the alleles APEX1 148Glu, IL8 -251 T and AADAT +401 T. These genotypic combinations were shown to affect cyto/chemokine levels and cell counts in CSF samples from BM patients. CONCLUSIONS In conclusion, this study revealed a significant association between genetic variability and altered inflammatory responses, involving important pathways that are activated during BM. This knowledge may be useful for a better understanding of BM pathogenesis and the development of new therapeutic approaches.

CCDC26, CDKN2BAS, RTEL1 and TERT Polymorphisms in pediatric brain tumor susceptibility.

The role of genetic polymorphisms in pediatric brain tumor (PBT) etiology is poorly understood. We hypothesized that single nucleotide polymorphisms (SNPs) identified in genome-wide association studies (GWAS) on adult glioma would also be associated with PBT risk. The study is based on the Cefalo study, a population-based multicenter case-control study. Saliva DNA from 245 cases and 489 controls, aged 7-19 years at diagnosis/reference date, was extracted and genotyped for 29 SNPs reported by GWAS to be significantly associated with risk of adult glioma. Data were analyzed using unconditional logistic regression. Stratified analyses were performed for two histological subtypes: astrocytoma alone and the other tumor types combined. The results indicated that four SNPs, CDKN2BAS rs4977756 (p = 0.036), rs1412829 (p = 0.037), rs2157719 (p = 0.018) and rs1063192 (p = 0.021), were associated with an increased susceptibility to PBTs, whereas the TERT rs2736100 was associated with a decreased risk (p = 0.018). Moreover, the stratified analyses showed a decreased risk of astrocytoma associated with RTEL1 rs6089953, rs6010620 and rs2297440 (p trend = 0.022, p trend = 0.042, p trend = 0.029, respectively) as well as an increased risk of this subtype associated with RTEL1 rs4809324 (p trend = 0.033). In addition, SNPs rs10464870 and rs891835 in CCDC26 were associated with an increased risk of non-astrocytoma tumor subtypes (p trend = 0.009, p trend = 0.007, respectively). Our findings indicate that SNPs in CDKN2BAS, TERT, RTEL1 and CCDC26 may be associated with the risk of PBTs. Therefore, we suggest that pediatric and adult brain tumors might share common genetic risk factors and similar etiological pathways.

The impact of ABCC11 polymorphisms on the risk of early-onset fluoropyrimidine toxicity

A missense variant (c.1637C>T, T546M) in ABCC11 encoding the MRP8 (multidrug resistance protein 8), a transporter of 5-fluorodeoxyuridine monophosphate, has been associated with an increased risk of 5-fluorouracil-related severe leukopenia. To validate this association, we investigated the impact of the ABCC11 variants c.1637C>T, c.538G>A and c.395+1087C>T on the risk of early-onset fluoropyrimidine-related toxicity in 514 cancer patients. The ABCC11 variant c.1637C>T was strongly associated with severe leukopenia in patients carrying risk variants in DPYD, encoding the key fluoropyrimidine-metabolizing enzyme dihydropyrimidine dehydrogenase (odds ratio (OR): 71.0; 95% confidence interval (CI): 2.5-2004.8; Pc.1637C>T*DPYD=0.013). In contrast, in patients without DPYD risk variants, no association with leukopenia (OR: 0.95; 95% CI: 0.34-2.6) or overall fluoropyrimidine-related toxicity (OR: 1.02; 95% CI: 0.5-2.1) was observed. Our study thus suggests that c.1637C>T affects fluoropyrimidine toxicity to leukocytes particularly in patients with high drug exposure, for example, because of reduced fluoropyrimidine catabolism.

Breast cancer risk and heterocyclic amines and effect modification by NAT gene polymorphisms

Breast cancer is the most common cancer in women in the United States and is a leading cause of cancer-related deaths (1). Recently, dietary heterocyclic amines (HCAs) have been proposed to be a risk factor for breast cancer (2). This study uses the data collected for a case-control study conducted at the M.D. Anderson Cancer Center to assess the association between breast cancer risk and HCAs {2-amino-1-methyl-6-phenylimidazole [4,5-b] pyridine (PhIP), 2-amino-3,8-dimethylimidazo [4,5-f] quinoxaline (MeIQx), 2-amino-3,4,8-trimethylimidazo [4,5-f] quinoxaline (DiMeIQx) and mutagenicity of HCAs} and to examine if this association is modified by genetic polymorphisms of N-acetyl transferases (NAT1/NAT2). The NAT1/2 genotype was determined using Taqman technology. HCAs were estimated by using a meat preparation questionnaire on meat type, cooking method, and doneness, combined with a quantitative HCA database. Three hundred and fifty patients with breast cancer attending the Diagnostic Radiology Clinic at M. D. Anderson Cancer Center and fulfilling the eligibility criteria were compared to three hundred and fifty patients attending the same clinic for benign breast lesions to answer these questions. Logistic regression models were used to control for known risk factors and showed no statistically significant association between breast cancer versus benign breast cancer lesions and dietary intake of heterocyclic amines. There was no clear difference in their effect after subgroup analyses in different acetylator strata of NAT1/2 and no statistical interactions were found between NAT1/2 genotypes and HCAs, suggesting no effect modification by NAT1/2 acetylator status. These results suggest the need for further research to analyze if these null associations were because of the benign breast lesions sharing the risk factors with breast cancer or any other factors which haven't been explored yet.^

DNA nucleotide excision repair gene single nucleotide polymorphisms and hereditary nonpolyposis colorectal cancer

Hereditary nonpolyposis colorectal cancer (HNPCC) is an autosomal dominant disease caused by germline mutations in DNA mismatch repair(MMR) genes. The nucleotide excision repair(NER) pathway plays a very important role in cancer development. We systematically studied interactions between NER and MMR genes to identify NER gene single nucleotide polymorphism (SNP) risk factors that modify the effect of MMR mutations on risk for cancer in HNPCC. We analyzed data from polymorphisms in 10 NER genes that had been genotyped in HNPCC patients that carry MSH2 and MLH1 gene mutations. The influence of the NER gene SNPs on time to onset of colorectal cancer (CRC) was assessed using survival analysis and a semiparametric proportional hazard model. We found the median age of onset for CRC among MMR mutation carriers with the ERCC1 mutation was 3.9 years earlier than patients with wildtype ERCC1(median 47.7 vs 51.6, log-rank test p=0.035). The influence of Rad23B A249V SNP on age of onset of HNPCC is age dependent (likelihood ratio test p=0.0056). Interestingly, using the likelihood ratio test, we also found evidence of genetic interactions between the MMR gene mutations and SNPs in ERCC1 gene(C8092A) and XPG/ERCC5 gene(D1104H) with p-values of 0.004 and 0.042, respectively. An assessment using tree structured survival analysis (TSSA) showed distinct gene interactions in MLH1 mutation carriers and MSH2 mutation carriers. ERCC1 SNP genotypes greatly modified the age onset of HNPCC in MSH2 mutation carriers, while no effect was detected in MLH1 mutation carriers. Given the NER genes in this study play different roles in NER pathway, they may have distinct influences on the development of HNPCC. The findings of this study are very important for elucidation of the molecular mechanism of colon cancer development and for understanding why some mutation carriers of the MSH2 and MLH1 gene develop CRC early and others never develop CRC. Overall, the findings also have important implications for the development of early detection strategies and prevention as well as understanding the mechanism of colorectal carcinogenesis in HNPCC. ^

Characterization of glucocorticoid receptor (NR3C1) gene polymorphisms and a family-based association study with hypertension

Hypertension is usually defined as having values of systolic blood pressure ≥140 mmHg, diastolic blood pressure ≥90 mmHg. Hypertension is one of the main adverse effects of glucocorticoid on the cardiovascular system. Glucocorticoids are essential hormones, secreted from adrenal glands in circadian fashion. Glucocorticoid's effect on blood pressure is conveyed by the glucocorticoid receptor (NR3C1), an omnipresent nuclear transcription factor. Although polymorphisms in this gene have long been implicated to be a causal factor for cardiovascular diseases such as hypertension, no study has yet thoroughly interrogated the gene's polymorphisms for their effect on blood pressure levels. Therefore, I have first resequenced ∼30 kb of the gene, encompassing all exons, promoter regions, 5'/3' UTRs as well as at least 1.5 kb of the gene's flanking regions from 114 chromosome 5 monosomic cell lines, comprised of three major American ethnic groups—European American, African American and Mexican American. I observed 115 polymorphisms and 14 common molecularly phased haplotypes. A subset of markers was chosen for genotyping study populations of GENOA (Genetic Epidemiology Network of Atherosclerosis; 1022 non-Hispanic whites, 1228 African Americans and 954 Mexican Americans). Since these study populations include sibships, the family-based association test was performed on 4 blood pressure-related quantitative variables—pulse, systolic blood pressure, diastolic blood pressure and mean arterial pressure. Using these analyses, multiple correlated SNPs are significantly protective against high systolic blood pressure in non-Hispanic whites, which includes rsb198, a SNP formerly associated with beneficial body compositions. Haplotype association analysis also supports this finding and all p-values remained significant after permutation tests. I therefore conclude that multiple correlated SNPs on the gene may confer protection against high blood pressure in non-Hispanic whites. ^

Folate pathway polymorphisms and neuropsychological impairment after leukemia therapy

Neuropsychological impairment occurs in 20%-40% of childhood acute lymphoblastic leukemia (ALL) survivors, possibly mediated by folate depletion following methotrexate chemotherapy. We evaluated the relationship between two folate pathway polymorphisms and neuropsychological impairment after childhood ALL chemotherapy. Eighty-six childhood ALL survivors were recruited between 2004-2007 at Texas Children's Hospital after exclusion for central nervous system leukemia, cranial irradiation, and age<1 year at diagnosis. Neuropsychological evaluation at a median of 5.3 years off therapy included a parental questionnaire and the following child performance measures: Trail Making Tests A and B, Grooved Pegboard Test Dominant-Hand and Nondominant-Hand, and Digit Span subtest. We performed genotyping for polymorphisms in two folate pathway genes: reduced folate carrier (RFC1 80G>A, rs1051266) and dihydrofolate reductase (DHFR Intron-1 19bp deletion). Fisher exact test, logistic regression, Student's t-test, and ANOVA were used to compare neuropsychological test scores by genotype, using a dominant model to group genotypes. In univariate analysis, survivors with cumulative methotrexate exposure ≥9000 mg/m2 had an increased risk of attention disorder (OR=6.2, 95% CI 1.2 – 31.3), compared to survivors with methotrexate exposure <9000 mg/m2. On average, female survivors scored 8.5 points higher than males on the Digit Span subtest, a test of working memory (p=0.02). The RFC1 80G>A and DHFR Intron-1 deletion polymorphisms were not related to attention disorder or impairment on tests of attention, processing speed, fine motor speed, or memory. These data imply a strong relationship between methotrexate dose intensity and impairment in attention after childhood ALL therapy. We did not find an association between the RFC1 80G>A or DHFR Intron-1 deletion polymorphisms and long-term neuropsychological impairment in childhood ALL survivors.^

Neurocognitive impairment in relation to glutathione S-transferase enzyme polymorphisms among medulloblastoma patients

Background. Medulloblastoma is a type of brain cancer that accounts for approximately 7-8% of all intracranial tumors and 20-30% of pediatric brain tumors. It is the most common type of malignant brain tumor in childhood. It was reported that majority of survivors with medulloblastoma have social problems, endocrine deficits, and neurological complications. Furthermore, all had significant deficits in neurocognitive functioning. Glutathione S-transferases belong to a family of isoenzymes that catalyze the glutathione conjugation of a variety of electrophilic compounds. ^ Objective. We aimed to determine whether the development of neurocognitive impairment is associated with GST polymorphisms among children and adolescents diagnosed with medulloblastoma (MB) after radiation therapy. ^ Methods. A pilot study composing of 16 children and adolescents diagnosed with MB at Texas Children's Cancer Center was conducted. The t-test was used to determine if the GST polymorphisms were related to neurocognitive impairment and logistic regression was performed to explore association between GST polymorphisms and gender, age at diagnosis, race/ethnicity, and risk group. ^ Results. An association was observed between GSTT1 polymorphism and cognitive impairment one year after radiation and GSTM1 polymorphism two years after radiation. It was observed that patients with GSTT1 null genotype have lower performance IQ (p=0.03) and full scale IQ (p=0.02) one year after radiation and patients with GSTM1 null genotype have lower verbal IQ (p=0.02) two years after radiation. Patients under age 8 have a statistically non-significant higher risk of having not null genotypes compared to those older than age 8 (OR= 7.5, 95%CI: 0.62-90.65 and OR= 2.63, 95%CI: 0.30-23.00 for GSTT1 and GSTM1 respectively). ^ Conclusion. There was a significant association between GSTT1 polymorphism and cognitive impairment one year after radiation and between GSTM1 polymorphism and cognitive impairment two years after radiation. Further large scale studies may be needed to confirm this finding and to examine the underlying mechanism of neurocognitive impairments after treatment of medulloblastoma patients.^

Single nucleotide polymorphisms (SNPs) associated with TGF-beta pathway and their significance in systemic sclerosis - A multilevel analysis

Systemic sclerosis (SSc) or Scleroderma is a complex disease and its etiopathogenesis remains unelucidated. Fibrosis in multiple organs is a key feature of SSc and studies have shown that transforming growth factor-β (TGF-β) pathway has a crucial role in fibrotic responses. For a complex disease such as SSc, expression quantitative trait loci (eQTL) analysis is a powerful tool for identifying genetic variations that affect expression of genes involved in this disease. In this study, a multilevel model is described to perform a multivariate eQTL for identifying genetic variation (SNPs) specifically associated with the expression of three members of TGF-β pathway, CTGF, SPARC and COL3A1. The uniqueness of this model is that all three genes were included in one model, rather than one gene being examined at a time. A protein might contribute to multiple pathways and this approach allows the identification of important genetic variations linked to multiple genes belonging to the same pathway. In this study, 29 SNPs were identified and 16 of them located in known genes. Exploring the roles of these genes in TGF-β regulation will help elucidate the etiology of SSc, which will in turn help to better manage this complex disease. ^