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. ^

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. ^

Single nucleotide polymorphism (SNP) selection for genotype-phenotype association studies

With hundreds of single nucleotide polymorphisms (SNPs) in a candidate gene and millions of SNPs across the genome, selecting an informative subset of SNPs to maximize the ability to detect genotype-phenotype association is of great interest and importance. In addition, with a large number of SNPs, analytic methods are needed that allow investigators to control the false positive rate resulting from large numbers of SNP genotype-phenotype analyses. This dissertation uses simulated data to explore methods for selecting SNPs for genotype-phenotype association studies. I examined the pattern of linkage disequilibrium (LD) across a candidate gene region and used this pattern to aid in localizing a disease-influencing mutation. The results indicate that the r2 measure of linkage disequilibrium is preferred over the common D′ measure for use in genotype-phenotype association studies. Using step-wise linear regression, the best predictor of the quantitative trait was not usually the single functional mutation. Rather it was a SNP that was in high linkage disequilibrium with the functional mutation. Next, I compared three strategies for selecting SNPs for application to phenotype association studies: based on measures of linkage disequilibrium, based on a measure of haplotype diversity, and random selection. The results demonstrate that SNPs selected based on maximum haplotype diversity are more informative and yield higher power than randomly selected SNPs or SNPs selected based on low pair-wise LD. The data also indicate that for genes with small contribution to the phenotype, it is more prudent for investigators to increase their sample size than to continuously increase the number of SNPs in order to improve statistical power. When typing large numbers of SNPs, researchers are faced with the challenge of utilizing an appropriate statistical method that controls the type I error rate while maintaining adequate power. We show that an empirical genotype based multi-locus global test that uses permutation testing to investigate the null distribution of the maximum test statistic maintains a desired overall type I error rate while not overly sacrificing statistical power. The results also show that when the penetrance model is simple the multi-locus global test does as well or better than the haplotype analysis. However, for more complex models, haplotype analyses offer advantages. The results of this dissertation will be of utility to human geneticists designing large-scale multi-locus genotype-phenotype association studies. ^

Influence of host interleukin-10 polymorphisms on development of traveler's diarrhea due to heat-labile enterotoxin-producing Escherichia coli in travelers from the United States who are visiting Mexico.

Up to 60% of U.S. visitors to Mexico develop traveler's diarrhea (TD), mostly due to enterotoxigenic Escherichia coli (ETEC) strains that produce heat-labile (LT) and/or heat-stable (ST) enterotoxins. Distinct single-nucleotide polymorphisms (SNPs) within the interleukin-10 (IL-10) promoter have been associated with high, intermediate, or low production of IL-10. We conducted a prospective study to investigate the association of SNPs in the IL-10 promoter and the occurrence of TD in ETEC LT-exposed travelers. Sera from U.S. travelers to Mexico collected on arrival and departure were studied for ETEC LT seroconversion by using cholera toxin as the antigen. Pyrosequencing was performed to genotype IL-10 SNPs. Stools from subjects who developed diarrhea were also studied for other enteropathogens. One hundred twenty-one of 569 (21.3%) travelers seroconverted to ETEC LT, and among them 75 (62%) developed diarrhea. Symptomatic seroconversion was more commonly seen in subjects who carried a genotype producing high levels of IL-10; it was seen in 83% of subjects with the GG genotype versus 54% of subjects with the AA genotype at IL-10 gene position -1082 (P, 0.02), in 71% of those with the CC genotype versus 33% of those with the TT genotype at position -819 (P, 0.005), and in 71% of those with the CC genotype versus 38% of those with the AA genotype at position -592 (P, 0.02). Travelers with the GCC haplotype were more likely to have symptomatic seroconversion than those with the ATA haplotype (71% versus 38%; P, 0.002). Travelers genetically predisposed to produce high levels of IL-10 were more likely to experience symptomatic ETEC TD.

Prognostic significance of xCT polymorphisms and expression in patients with advanced pancreatic cancer treated with chemotherapy

The plasma membrane xc- cystine/glutamate transporter mediates cellular uptake of cystine in exchange for intracellular glutamate and is highly expressed by pancreatic cancer cells. The xCT gene, encoding the cystine-specific xCT protein subunit of xc-, is important in regulating intracellular glutathione (GSH) levels, critical for cancer cell protection against oxidative stress, tumor growth and resistance to chemotherapeutic agents including platinum. We examined 4 single nucleotide polymorphisms (SNPs) of the xCT gene in 269 advanced pancreatic cancer patients who received first line gemcitabine with or without cisplatin or oxaliplatin. Genotyping was performed using Taqman real-time PCR assays. A statistically significant correlation was noted between the 3' untranslated region (UTR) xCT SNP rs7674870 and overall survival (OS): Median survival time (MST) was 10.9 and 13.6 months, respectively, for the TT and TC/CC genotypes (p = 0.027). Stratified analysis showed the genotype effect was significant in patients receiving gemcitabine in combination with platinum therapy (n = 145): MST was 10.5 versus 14.1 months for the TT and TC/CC genotypes, respectively (p = 0.013). The 3' UTR xCT SNP rs7674870 may correlate with OS in pancreatic cancer patients receiving gemcitabine and platinum combination therapy. Paraffin-embedded core and surgical biopsy tumor specimens from 98 patients with metastatic pancreatic adenocarcinoma were analyzed by immunohistochemistry using an xCT specific antibody. xCT protein IHC expression scores were analyzed in relation to overall survival in 86 patients and genotype in 12 patients and no statistically significant association was found between the level of xCT IHC expression score and overall survival (p = 0.514). When xCT expression was analyzed in terms of treatment response, no statistically significant associations could be determined (p = 0.908). These data suggest that polymorphic variants of xCT may have predictive value, and that the xc- transporter may represent an important target for therapy in pancreatic cancer.

Association of p21/p27 polymorphisms with risk of HPV16-associated oral squamous cell carcinoma

The increasing incidence of oral squamous cell carcinoma (OSCC) among young adults has been associated with sexually transmitted infection of human papillomavirus (HPV), particularly HPV16. Given the roles of p21 (WAF1/Cip1/CDKN1A) and p27 (Kip1/CDKNIB) in cell-cycle regulation and of HPV16 E6 and E7 oncoproteins in p53 degradation and pRb inactivation, the effect of HPV16 L1 seropositivity and three putatively functional single-nucleotide polymorphisms (SNPs) of p21 (p21 C70T and p21 C98A) and p27 (p27 T109G), individually and in combination, on the risk of OSCC was evaluated in a hospital-based case-control study of 327 cases and 401 cancer-free controls who were frequency-matched on age, gender and smoking status. Individuals with HPV16 L1 seropositivity had an overall 3-fold increased risk of having OSCC than those with HPV16 seronegativity. The increased risk of HPV16-associated OSCC was particularly found among younger people (aged ≤ 50 years), males, never smokers, never drinkers and oropharynx cancer patients. None of three p21 and p27 polymorphisms alone was significantly associated with risk of OSCC. Individuals with variant genotypes for both p21 polymorphisms were more likely to have OSCC than individuals with wild-type genotypes or variant genotypes for either one of the p21 polymorphisms (adjusted OR, 1.4; 95% CI, 0.9-2.1). There was a borderline significant or significant interaction between the p21 C70T, combined p21 and combined p21/p27 genotypes and HPV16 L1 seropositivity on risk of OSCC. The three studied p21 and p27 polymorphisms, individually or in combination, did not appear to have an effect on HPV16-related clinical outcomes (overall and disease-free survival and tumor recurrence). Despite the fact that the exact biological mechanism remains to be explored, these findings suggest possible involvement of p21variants, particularly the p21 C70T variant genotypes (CT/TT), in the etiology of HPV16-associated OPSCC. Further large and functional studies are required to validate the findings.^

Role of obesity, weight gain, physical activity, and polymorphisms in the mTOR pathway and the risk of renal cell carcinoma

The interplay between obesity, physical activity, weight gain and genetic variants in mTOR pathway have not been studied in renal cell carcinoma (RCC). We examined the associations between obesity, weight gain, physical activity and RCC risk. We also analyzed whether genetic variants in the mTOR pathway could modify the association. Incident renal cell carcinoma cases and healthy controls were recruited from the University of Texas MD Anderson Cancer Center in Houston, Texas. Cases and controls were frequency-matched by age (±5 years), ethnicity, sex, and county of residence. Epidemiologic data were collected via in-person interview. A total of 577 cases and 593 healthy controls (all white) were included. One hundred ninety-two (192) SNPs from 22 genes were available and their genotyping data were extracted from previous genome-wide association studies. Logistic regression and regression spline were performed to obtain odds ratios. Obesity at age 20, 40, and 3 years prior to diagnosis/recruitment, and moderate and large weight gain from age 20 to 40 were each significantly associated with increased RCC risk. Low physical activity was associated with a 4.08-fold (95% CI: 2.92-5.70) increased risk. Five single nucleotide polymorphisms (SNPs) were significantly associated with RCC risk and their cumulative effect increased the risk by up to 72% (95% CI: 1.20-2.46). Strata specific effects for weight change and genotyping cumulative groups were observed. However, no interaction was suggested by our study. In conclusion, energy balance related risk factors and genetic variants in the mTOR pathway may jointly influence susceptibility to RCC. ^

Association of retinoic acid receptor genes with meningomyelocele.

BACKGROUND: Neural tube defects (NTDs) occur in as many as 0.5-2 per 1000 live births in the United States. One of the most common and severe neural tube defects is meningomyelocele (MM) resulting from failed closure of the caudal end of the neural tube. MM has been induced by retinoic acid teratogenicity in rodent models. We hypothesized that genetic variants influencing retinoic acid (RA) induction via retinoic acid receptors (RARs) may be associated with risk for MM. METHODS: We analyzed 47 single nucleotide polymorphisms (SNPs) that span across the three retinoic acid receptor genes using the SNPlex genotyping platform. Our cohort consisted of 610 MM families. RESULTS: One variant in the RARA gene (rs12051734), three variants in the RARB gene (rs6799734, rs12630816, rs17016462), and a single variant in the RARG gene (rs3741434) were found to be statistically significant at p < 0.05. CONCLUSION: RAR genes were associated with risk for MM. For all associated SNPs, the rare allele conferred a protective effect for MM susceptibility.

HLA-DPB1 and DPB2 are genetic loci for systemic sclerosis: a genome-wide association study in Koreans with replication in North Americans.

OBJECTIVE: To identify systemic sclerosis (SSc) susceptibility loci via a genome-wide association study. METHODS: A genome-wide association study was performed in 137 patients with SSc and 564 controls from Korea using the Affymetrix Human SNP Array 5.0. After fine-mapping studies, the results were replicated in 1,107 SSc patients and 2,747 controls from a US Caucasian population. RESULTS: The single-nucleotide polymorphisms (SNPs) (rs3128930, rs7763822, rs7764491, rs3117230, and rs3128965) of HLA-DPB1 and DPB2 on chromosome 6 formed a distinctive peak with log P values for association with SSc susceptibility (P=8.16x10(-13)). Subtyping analysis of HLA-DPB1 showed that DPB1*1301 (P=7.61x10(-8)) and DPB1*0901 (P=2.55x10(-5)) were the subtypes most susceptible to SSc in Korean subjects. In US Caucasians, 2 pairs of SNPs, rs7763822/rs7764491 and rs3117230/rs3128965, showed strong association with SSc patients who had either circulating anti-DNA topoisomerase I (P=7.58x10(-17)/4.84x10(-16)) or anticentromere autoantibodies (P=1.12x10(-3)/3.2x10(-5)), respectively. CONCLUSION: The results of our genome-wide association study in Korean subjects indicate that the region of HLA-DPB1 and DPB2 contains the loci most susceptible to SSc in a Korean population. The confirmatory studies in US Caucasians indicate that specific SNPs of HLA-DPB1 and/or DPB2 are strongly associated with US Caucasian patients with SSc who are positive for anti-DNA topoisomerase I or anticentromere autoantibodies.

Genes in glucose metabolism and association with spina bifida.

The authors test single nucleotide polymorphisms (SNPs) in coding sequences of 12 candidate genes involved in glucose metabolism and obesity for associations with spina bifida. Genotyping was performed on 507 children with spina bifida and their parents plus anonymous control DNAs from Hispanic and Caucasian individuals. The transmission disequilibrium test was performed to test for genetic associations between transmission of alleles and spina bifida in the offspring (P < .05). A statistically significant association between Lys481 of HK1 (G allele), Arg109Lys of LEPR (G allele), and Pro196 of GLUT1 (A allele) was found ( P = .019, .039, and .040, respectively). Three SNPs on 3 genes involved with glucose metabolism and obesity may be associated with increased susceptibility to spina bifida.

THE EXPRESSION PROFILE OF A GWAS DETECTED GENE (NINJ2) IN THE BRAIN OF A MOUSE MODEL OF ISCHEMIC STROKE

Stroke is the third leading cause of death and a major debilitating disease in the United States. Multiple factors, including genetic factors, contribute to the development of the disease. Genome-wide association studies (GWAS) have contributed to the identification of genetic loci influencing risk for complex diseases, such as stroke. In 2010, a GWAS of incident stroke was performed in four large prospective cohorts from the USA and Europe and identified an association of two Single Nucleotide Polymorphisms (SNPs) on chromosome 12p13 with a greater risk of ischemic stroke in individuals of European and African-American ancestry. These SNPs are located 11 Kb upstream of the nerve injury-induced gene 2, Ninjurin2 (NINJ2), suggesting that this gene may be involved in stroke pathogenesis. NINJ2 is a cell adhesion molecule induced in the distal glial cells from a sciatic-nerve injury at 7-days after injury. In an effort to ascribe a possible role of NINJ2 in stroke, we have assessed changes in the level of gene and protein expression of NINJ2 following a time-course from a transiently induced middle cerebral artery ischemic stroke in mice brains. We report an increase in the gene expression of NINJ2 in the ischemic and peri-infarct (ipsilateral) cortical tissues at 7 and 14-days after stroke. We also report an increase in the protein expression of NINJ2 in the cortex of both the ipsilateral and contralateral cortical tissues at the same time-points. We conclude that the expression of NINJ2 is regulated by an ischemic stroke in the cortex and is consistent with NINJ2 being involved in the recovery time-points of stroke.

Theoretical and experimental studies of linkage disequilibrium in human populations

Linkage disequilibrium (LD) is defined as the nonrandom association of alleles at two or more loci in a population and may be a useful tool in a diverse array of applications including disease gene mapping, elucidating the demographic history of populations, and testing hypotheses of human evolution. However, the successful application of LD-based approaches to pertinent genetic questions is hampered by a lack of understanding about the forces that mediate the genome-wide distribution of LD within and between human populations. Delineating the genomic patterns of LD is a complex task that will require interdisciplinary research that transcends traditional scientific boundaries. The research presented in this dissertation is predicated upon the need for interdisciplinary studies and both theoretical and experimental projects were pursued. In the theoretical studies, I have investigated the effect of genotyping errors and SNP identification strategies on estimates of LD. The primary importance of these two chapters is that they provide important insights and guidance for the design of future empirical LD studies. Furthermore, I analyzed the allele frequency distribution of 26,530 single nucleotide polymorphisms (SNPs) in three populations and generated the first-generation natural selection map of the human genome, which will be an important resource for explaining and understanding genomic patterns of LD. Finally, in the experimental study, I describe a novel and simple, low-cost, and high-throughput SNP genotyping method. The theoretical analyses and experimental tools developed in this dissertation will facilitate a more complete understanding of patterns of LD in human populations. ^

Nonlinear tests for genetic studies of complex disease

Linkage and association studies are major analytical tools to search for susceptibility genes for complex diseases. With the availability of large collection of single nucleotide polymorphisms (SNPs) and the rapid progresses for high throughput genotyping technologies, together with the ambitious goals of the International HapMap Project, genetic markers covering the whole genome will be available for genome-wide linkage and association studies. In order not to inflate the type I error rate in performing genome-wide linkage and association studies, multiple adjustment for the significant level for each independent linkage and/or association test is required, and this has led to the suggestion of genome-wide significant cut-off as low as 5 × 10 −7. Almost no linkage and/or association study can meet such a stringent threshold by the standard statistical methods. Developing new statistics with high power is urgently needed to tackle this problem. This dissertation proposes and explores a class of novel test statistics that can be used in both population-based and family-based genetic data by employing a completely new strategy, which uses nonlinear transformation of the sample means to construct test statistics for linkage and association studies. Extensive simulation studies are used to illustrate the properties of the nonlinear test statistics. Power calculations are performed using both analytical and empirical methods. Finally, real data sets are analyzed with the nonlinear test statistics. Results show that the nonlinear test statistics have correct type I error rates, and most of the studied nonlinear test statistics have higher power than the standard chi-square test. This dissertation introduces a new idea to design novel test statistics with high power and might open new ways to mapping susceptibility genes for complex diseases. ^

A coalescent analysis for modeling the mutation process in colorectal cancer

Colorectal cancer is the forth most common diagnosed cancer in the United States. Every year about a hundred forty-seven thousand people will be diagnosed with colorectal cancer and fifty-six thousand people lose their lives due to this disease. Most of the hereditary nonpolyposis colorectal cancer (HNPCC) and 12% of the sporadic colorectal cancer show microsatellite instability. Colorectal cancer is a multistep progressive disease. It starts from a mutation in a normal colorectal cell and grows into a clone of cells that further accumulates mutations and finally develops into a malignant tumor. In terms of molecular evolution, the process of colorectal tumor progression represents the acquisition of sequential mutations. ^ Clinical studies use biomarkers such as microsatellite or single nucleotide polymorphisms (SNPs) to study mutation frequencies in colorectal cancer. Microsatellite data obtained from single genome equivalent PCR or small pool PCR can be used to infer tumor progression. Since tumor progression is similar to population evolution, we used an approach known as coalescent, which is well established in population genetics, to analyze this type of data. Coalescent theory has been known to infer the sample's evolutionary path through the analysis of microsatellite data. ^ The simulation results indicate that the constant population size pattern and the rapid tumor growth pattern have different genetic polymorphic patterns. The simulation results were compared with experimental data collected from HNPCC patients. The preliminary result shows the mutation rate in 6 HNPCC patients range from 0.001 to 0.01. The patients' polymorphic patterns are similar to the constant population size pattern which implies the tumor progression is through multilineage persistence instead of clonal sequential evolution. The results should be further verified using a larger dataset. ^

Diabetes genes and risk of prostate cancer in the Atherosclerosis Risk in Communities study

Prostate cancer (PrCa) is a leading cause of morbidity and mortality, yet the etiology remains uncertain. Meta-analyses show that PrCa risk is reduced by 16% in men with type 2 diabetes (T2D), but the mechanism is unknown. Recent genome-wide association studies and meta-analyses have found single nucleotide polymorphisms (SNPs) that consistently predict T2D risk. We evaluated associations of incident PrCa with 14 T2D SNPs in the Atherosclerosis Risk in Communities (ARIC) study. From 1987-2000, there were 397 incident PrCa cases ascertained from state or local cancer registries among 6,642 men (1,560 blacks and 5,082 whites) aged 45-64 years at baseline. Genotypes were determined by TaqMan assay. Cox proportional hazards models were used to assess the association between PrCa and increasing number of T2D risk-raising alleles for individual SNPs and for genetic risk scores (GRS) comprised of the number of T2D risk-raising alleles across SNPs. Two-way gene-gene interactions were evaluated with likelihood ratio tests. Using additive genetic models, the T2D risk-raising allele was associated with significantly reduced risk of PrCa for IGF2BP2 rs4402960 (hazard ratio [HR]=0.79; P=0.07 among blacks only), SLC2A2 rs5400 (race-adjusted HR=0.85; P=0.05) and UCP2 rs660339 (race-adjusted HR=0.84; P=0.02), but significantly increased risk of PrCa for CAPN10 rs3792267 (race-adjusted HR=1.20; P=0.05). No other SNPs were associated with PrCa using an additive genetic model. However, at least one copy of the T2D risk-raising allele for TCF7L2 rs7903146 was associated with reduced PrCa risk using a dominant genetic model (race-adjusted HR=0.79; P=0.03). These results imply that the T2D-PrCa association may be partly due to shared genetic variation, but these results should be verified since multiple tests were performed. When the combined, additive effects of these SNPs were tested using a GRS, there was nearly a 10% reduction in risk of PrCa per T2D risk-raising allele (race-adjusted HR=0.92; P=0.02). SNPs in IGF2BP2, KCNJ11 and SLC2A2 were also involved in multiple synergistic gene-gene interactions on a multiplicative scale. In conclusion, it appears that the T2D-PrCa association may be due, in part, to common genetic variation. Further knowledge of T2D gene-PrCa mechanisms may improve understanding of PrCa etiology and may inform PrCa prevention and treatment.^