Molecular epidemiological studies on genetic predispositions to squamous cell carcinoma of the head and neck

Two molecular epidemiological studies were conducted to examine associations between genetic variation and risk of squamous cell carcinoma of the head and neck (SCCHN). In the first study, we hypothesized that genetic variation in p53 response elements (REs) may play roles in the etiology of SCCHN. We selected and genotyped five polymorphic p53 REs as well as a most frequently studied p53 codon 72 (Arg72Pro, rs1042522) polymorphism in 1,100 non-Hispanic White SCCHN patients and 1,122 age-and sex-matched cancer-free controls recruited at The University of Texas M. D. Anderson Cancer Center. In multivariate logistic regression analysis with adjustment for age, sex, smoking and drinking status, marital status and education level, we observed that the EOMES rs3806624 CC genotype had a significant effect of protection against SCCHN risk (adjusted odds ratio= 0.79, 95% confidence interval =0.64–0.98), compared with the -838TT+CT genotypes. Moreover, a significantly increased risk associated with the combined genotypes of p53 codon 72CC and EOMES -838TT+CT was observed, especially in the subgroup of non-oropharyneal cancer patients. The values of false-positive report probability were also calculated for significant findings. In the second study, we assessed the association between SCCHN risk and four potential regulatory single nucleotide polymorphisms (SNPs) of DEC1 (deleted in esophageal cancer 1) gene, a candidate tumor suppressor gene for esophageal cancer. After adjustment for age, sex, and smoking and drinking status, the variant -606CC (i.e., -249CC) homozygotes had a significantly reduced SCCHN risk (adjusted odds ratio = 0.71, 95% confidence interval = 0.52–0.99), compared with the -606TT homozygotes. Stratification analyses showed that a reduced risk associated with the -606CC genotype was more pronounced in subgroups of non-smokers, non-drinkers, younger subjects (defined as ≤ 57 years), carriers of TP53 Arg/Arg (rs1042522) genotype, patients with oropharyngeal cancer or late-stage SCCHN. Further in silico analysis revealed that the -249 T-to-C change led to a gain of a transcription factor binding site. Additional functional analysis showed that the -249T-to-C change significantly enhanced transcriptional activity of the DEC1 promoter and the DNA-protein binding activity. We conclude that the DEC1 promoter -249 T>C (rs2012775) polymorphism is functional, modulating susceptibility to SCCHN among non-Hispanic Whites. Additional large-scale, preferably population-based studies are needed to validate our findings.^

Genetic variations in the PI3K-AKT-mTOR pathway and bladder cancer susceptibility and clinical outcomes

Bladder cancer is the fourth most common cancer in men in the United States. There is compelling evidence supporting that genetic variations contribute to the risk and outcomes of bladder cancer. The PI3K-AKT-mTOR pathway is a major cellular pathway involved in proliferation, invasion, inflammation, tumorigenesis, and drug response. Somatic aberrations of PI3K-AKT-mTOR pathway are frequent events in several cancers including bladder cancer; however, no studies have investigated the role of germline genetic variations in this pathway in bladder cancer. In this project, we used a large case control study to evaluate the associations of a comprehensive catalogue of SNPs in this pathway with bladder cancer risk and outcomes. Three SNPs in RAPTOR were significantly associated with susceptibility: rs11653499 (OR: 1.79, 95%CI: 1.24–2.60), rs7211818 (OR: 2.13, 95%CI: 1.35–3.36), and rs7212142 (OR: 1.57, 95%CI: 1.19–2.07). Two haplotypes constructed from these 3 SNPs were also associated with bladder cancer risk. In combined analysis, a significant trend was observed for increased risk with an increase in the number of unfavorable genotypes (P for trend<0.001). Classification and regression tree analysis identified potential gene-environment interactions between RPS6KA5 rs11653499 and smoking. In superficial bladder cancer, we found that PTEN rs1234219 and rs11202600, TSC1 rs7040593, RAPTOR rs901065, and PIK3R1 rs251404 were significantly associated with recurrence in patients receiving BCG. In muscle invasive and metastatic bladder cancer, AKT2 rs3730050, PIK3R1 rs10515074, and RAPTOR rs9906827 were associated with survival. Survival tree analysis revealed potential gene-gene interactions: patients carrying the unfavorable genotypes of PTEN rs1234219 and TSC1 rs704059 exhibited a 5.24-fold (95% CI: 2.44–11.24) increased risk of recurrence. In combined analysis, with the increasing number of unfavorable genotypes, there was a significant trend of higher risk of recurrence and death (P for trend<0.001) in Cox proportional hazard regression analysis, and shorter event (recurrence and death) free survival in Kaplan-Meier estimates (P log rank<0.001). This study strongly suggests that genetic variations in PI3K-AKT-mTOR pathway play an important role in bladder cancer development. The identified SNPs, if validated in further studies, may become valuable biomarkers in assessing an individual's cancer risk, predicting prognosis and treatment response, and facilitating physicians to make individualized treatment decisions. ^

Association of folate metabolic pathway genes with human spina bifida meningomyelocele

Neural tube defects including spina bifida meningomyelocele (SBMM) are common malformations of the brain and spinal cord, and include all abnormalities resulting from lack of closure of the developing neural tube during embryological development.^ The specific aims of this study were to determine if single nucleotide polymorphic variants (SNPs) in the folate/homocysteine metabolic pathway genes confer a risk for NTD susceptibility within this SBMM population.^ In completion of the first specific aim, two novel SNPs were identified in the FOLR1 gene in Chromosome 11of patients including one in non-coding exon 1 with a C → T transition at nucleotide position 71578317 and another in non-coding exon 3 with a T → G transversion at nucleotide position 71579123. It will be important to determine if these variants are present in the respective parents of these individuals. If they are in fact de novo variants, then these SNPs may be more likely to contribute to the birth defect.^ The second project aim was to analyze genotypes associated with SBMM risk by transmission disequilibrium tests (TDT) and association was detected on several SNPs across the folate metabolic pathway genes in this population. SNPs with significant RC-TDT values were found within the DHFR gene (rs1650723), the MTRR gene (rs327592), the FOLR2 gene (rs13908), four tightly linked variants in the FOLR3 gene (rs7925545, rs7926875, rs7926987, rs7926360) and a variant in the SLC19A1 gene (rs1888530). The product of each of these genes performs a vital function in the folate metabolic pathway. It is conceivable, therefore, that if the individual SNP or SNPs can be proven to perturb the function in some way that they may be involved in the disruption of folate metabolism and in the resulting birth defect. Validating the results of this study in other independent populations will further strengthen the evidence that dysfunction of folate enzymes and receptors may confer SBMM risk in humans. ^

Gene-gene interaction and gene-pathway analysis of genome-wide association for schizophrenia

Schizophrenia (SZ) is a complex disorder with high heritability and variable phenotypes that has limited success in finding causal genes associated with the disease development. Pathway-based analysis is an effective approach in investigating the molecular mechanism of susceptible genes associated with complex diseases. The etiology of complex diseases could be a network of genetic factors and within the genes, interaction may occur. In this work we argue that some genes might be of small effect that by itself are neither sufficient nor necessary to cause the disease however, their effect may induce slight changes to the gene expression or affect the protein function, therefore, analyzing the gene-gene interaction mechanism within the disease pathway would play crucial role in dissecting the genetic architecture of complex diseases, making the pathway-based analysis a complementary approach to GWAS technique. ^ In this study, we implemented three novel linkage disequilibrium based statistics, the linear combination, the quadratic, and the decorrelation test statistics, to investigate the interaction between linked and unlinked genes in two independent case-control GWAS datasets for SZ including participants of European (EA) and African (AA) ancestries. The EA population included 1,173 cases and 1,378 controls with 729,454 genotyped SNPs, while the AA population included 219 cases and 288 controls with 845,814 genotyped SNPs. We identified 17,186 interacting gene-sets at significant level in EA dataset, and 12,691 gene-sets in AA dataset using the gene-gene interaction method. We also identified 18,846 genes in EA dataset and 19,431 genes in AA dataset that were in the disease pathways. However, few genes were reported of significant association to SZ. ^ Our research determined the pathways characteristics for schizophrenia through the gene-gene interaction and gene-pathway based approaches. Our findings suggest insightful inferences of our methods in studying the molecular mechanisms of common complex diseases.^

Genotypic spectrum and genotype-phenotype correlation of trimethylaminuria

Trimethylaminuria (TMAU) or Fish odor syndrome is an autosomal recessive disease that is characterized by pungent body odor with subsequent psychosocial complications. There are limited studies of the sequence variants causing TMAU in the literature with most studies describing only one or two patients and lacking genotype-phenotype correlations. Also to date, there is no laboratory in the US or Europe that offers TMA genetic testing on a clinical basis. We have recently validated genetic testing in the University of Colorado DNA Diagnostic Laboratory. We have a database of a few dozen patients with a biochemical diagnosis of TMA at the University of Colorado at Denver Health Sciences Center (UCDHSC) which includes a few patients with the classical form of the disease. We have used the newly established clinical test in our institution to attempt to characterize the genotype (sequence variants including mutations and polymorphisms) of classical TMAU patients and to establish a genotype-phenotype (biochemical and clinical) association. The questionnaire results confirmed most of the previously reported epidemiological findings of TMAU and also indicated that TMAU patients use multiple intervention measures in attempt to control their symptoms with dietary control being most effective. Despite the complexity of intervention, most patients did not have any medical follow up and there was underutilization of specialist care. In a set of our patients, two deleterious mutations were identified in 4/12 patients including a novel T237P sequence variant, while the majority of our patients (8/12) did not reveal any mutations. Some of the latter were double heterozygous for the E158K and E308G polymorphisms which could explain a mild phenotype while others had only the E158K variant which raised the question of undetected mutations. These results indicate that further experiments are needed to further delineate the full mutational spectrum of the FMO3 gene. ^

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

Identification of diabetic retinopathy genes through a genome-wide association study among Mexican-Americans from Starr County, Texas

To identify genetic susceptibility loci for severe diabetic retinopathy, 286 Mexican-Americans with type 2 diabetes from Starr County, Texas completed detailed physical and ophthalmologic examinations including fundus photography for diabetic retinopathy grading. 103 individuals with moderate-to-severe non-proliferative diabetic retinopathy or proliferative diabetic retinopathy were defined as cases for this study. DNA samples extracted from study subjects were genotyped using the Affymetrix GeneChip® Human Mapping 100K Set, which includes 116,204 single nucleotide polymorphisms (SNPs) across the whole genome. Single-marker allelic tests and 2- to 8-SNP sliding-window Haplotype Trend Regression implemented in HelixTreeTM were first performed with these direct genotypes to identify genes/regions contributing to the risk of severe diabetic retinopathy. An additional 1,885,781 HapMap Phase II SNPs were imputed from the direct genotypes to expand the genomic coverage for a more detailed exploration of genetic susceptibility to diabetic retinopathy. The average estimated allelic dosage and imputed genotypes with the highest posterior probabilities were subsequently analyzed for associations using logistic regression and Fisher's Exact allelic tests, respectively. To move beyond these SNP-based approaches, 104,572 directly genotyped and 333,375 well-imputed SNPs were used to construct genetic distance matrices based on 262 retinopathy candidate genes and their 112 related biological pathways. Multivariate distance matrix regression was then used to test hypotheses with genes and pathways as the units of inference in the context of susceptibility to diabetic retinopathy. This study provides a framework for genome-wide association analyses, and implicated several genes involved in the regulation of oxidative stress, inflammatory processes, histidine metabolism, and pancreatic cancer pathways associated with severe diabetic retinopathy. Many of these loci have not previously been implicated in either diabetic retinopathy or diabetes. In summary, CDC73, IL12RB2, and SULF1 had the best evidence as candidates to influence diabetic retinopathy, possibly through novel biological mechanisms related to VEGF-mediated signaling pathway or inflammatory processes. While this study uncovered some genes for diabetic retinopathy, a comprehensive picture of the genetic architecture of diabetic retinopathy has not yet been achieved. Once fully understood, the genetics and biology of diabetic retinopathy will contribute to better strategies for diagnosis, treatment and prevention of this disease.^

Systematic review of genetic risk score in coronary heart disease and other diseases

In order to better take advantage of the abundant results from large-scale genomic association studies, investigators are turning to a genetic risk score (GRS) method in order to combine the information from common modest-effect risk alleles into an efficient risk assessment statistic. The statistical properties of these GRSs are poorly understood. As a first step toward a better understanding of GRSs, a systematic analysis of recent investigations using a GRS was undertaken. GRS studies were searched in the areas of coronary heart disease (CHD), cancer, and other common diseases using bibliographic databases and by hand-searching reference lists and journals. Twenty-one independent case-control studies, cohort studies, and simulation studies (12 in CHD, 9 in other diseases) were identified. The underlying statistical assumptions of the GRS using the experience of the Framingham risk score were investigated. Improvements in the construction of a GRS guided by the concept of composite indicators are discussed. The GRS will be a promising risk assessment tool to improve prediction and diagnosis of common diseases.^

Genetic variants in the mammalian target of rapamycin (mTOR) signaling pathway as predictors of survival and clinical response in women with ovarian cancer

Background. The mTOR pathway is commonly altered in human tumors and promotes cell survival and proliferation. Preliminary evidence suggests this pathway's involvement in chemoresistance to platinum and taxanes, first line therapy for epithelial ovarian cancer. A pathway-based approach was used to identify individual germline single nucleotide polymorphisms (SNPs) and cumulative effects of multiple genetic variants in mTOR pathway genes and their association with clinical outcome in women with ovarian cancer. ^ Methods. The case-series was restricted to 319 non-Hispanic white women with high grade ovarian cancer treated with surgery and platinum-based chemotherapy. 135 SNPs in 20 representative genes in the mTOR pathway were genotyped. Hazard ratios (HRs) for death and Odds ratios (ORs) for failure to respond to primary therapy were estimated for each SNP using the multivariate Cox proportional hazards model and multivariate logistic regression model, respectively, while adjusting for age, stage, histology and treatment sequence. A survival tree analysis of SNPs with a statistically significant association (p<0.05) was performed to identify higher order gene-gene interactions and their association with overall survival. ^ Results. There was no statistically significant difference in survival by tumor histology or treatment regimen. The median survival for the cohort was 48.3 months. Seven SNPs were significantly associated with decreased survival. Compared to those with no unfavorable genotypes, the HR for death increased significantly with the increasing number of unfavorable genotypes and women in the highest risk category had HR of 4.06 (95% CI 2.29–7.21). The survival tree analysis also identified patients with different survival patterns based on their genetic profiles. 13 SNPs on five different genes were found to be significantly associated with a treatment response, defined as no evidence of disease after completion of primary therapy. Rare homozygous genotype of SNP rs6973428 showed a 5.5-fold increased risk compared to the wild type carrying genotypes. In the cumulative effect analysis, the highest risk group (individuals with ≥8 unfavorable genotypes) was significantly less likely to respond to chemotherapy (OR=8.40, 95% CI 3.10–22.75) compared to the low risk group (≤4 unfavorable genotypes). ^ Conclusions. A pathway-based approach can demonstrate cumulative effects of multiple genetic variants on clinical response to chemotherapy and survival. Therapy targeting the mTOR pathway may modify outcome in select patients.^

Longitudinal and logistic analysis of endocrine hormone profiles, physical maturation, and candidate genes related to pubertal events in children

Children who experience early pubertal development have an increased risk of developing cancer (breast, ovarian, and testicular), osteoporosis, insulin resistance, and obesity as adults. Early pubertal development has been associated with depression, aggressiveness, and increased sexual prowess. Possible explanations for the decline in age of pubertal onset include genetics, exposure to environmental toxins, better nutrition, and a reduction in childhood infections. In this study we (1) evaluated the association between 415 single nucleotide polymorphisms (SNPs) from hormonal pathways and early puberty, defined as menarche prior to age 12 in females and Tanner Stage 2 development prior to age 11 in males, and (2) measured endocrine hormone trajectories (estradiol, testosterone, and DHEAS) in relation to age, race, and Tanner Stage in a cohort of children from Project HeartBeat! At the end of the 4-year study, 193 females had onset of menarche and 121 males had pubertal staging at age 11. African American females had a younger mean age at menarche than Non-Hispanic White females. African American females and males had a lower mean age at each pubertal stage (1-5) than Non-Hispanic White females and males. African American females had higher mean BMI measures at each pubertal stage than Non-Hispanic White females. Of the 415 SNPs evaluated in females, 22 SNPs were associated with early menarche, when adjusted for race ( p<0.05), but none remained significant after adjusting for multiple testing by False Discovery Rate (p<0.00017). In males, 17 SNPs were associated with early pubertal development when adjusted for race (p<0.05), but none remained significant when adjusted for multiple testing (p<0.00017). ^ There were 4955 hormone measurements taken during the 4-year study period from 632 African American and Non-Hispanic White males and females. On average, African American females started and ended the pubertal process at a younger age than Non-Hispanic White females. The mean age of Tanner Stage 2 breast development in African American and Non-Hispanic White females was 9.7 (S.D.=0.8) and 10.2 (S.D.=1.1) years, respectively. There was a significant difference by race in mean age for each pubertal stage, except Tanner Stage 1 for pubic hair development. Both Estradiol and DHEAS levels in females varied significantly with age, but not by race. Estradiol and DHEAS levels increased from Tanner Stage 1 to Tanner Stage 5.^ African American males had a lower mean age at each Tanner Stage of development than Non-Hispanic White males. The mean age of Tanner Stage 2 genital development in African American and Non-Hispanic White males was 10.5 (S.D.=1.1) and 10.8 (S.D.=1.1) years, respectively, but this difference was not significant (p=0.11). Testosterone levels varied significantly with age and race. Non-Hispanic White males had higher levels of testosterone than African American males from Tanner Stage 1-4. Testosterone levels increased for both races from Tanner Stage 1 to Tanner Stage 5. Testosterone levels had the steepest increase from ages 11-15 for both races. DHEAS levels in males varied significantly with age, but not by race. DHEAS levels had the steepest increase from ages 14-17. ^ In conclusion, African American males and females experience pubertal onset at a younger age than Non-Hispanic White males and females, but in this study, we could not find a specific gene that explained the observed variation in age of pubertal onset. Future studies with larger study populations may provide a better understanding of the contribution of genes in early pubertal onset.^

Pathway analysis of genome-wide association study data in primary biliary cirrhosisa

Genome-wide association studies (GWAS) have successfully identified several genetic loci associated with inherited predisposition to primary biliary cirrhosis (PBC), the most common autoimmune disease of the liver. Pathway-based tests constitute a novel paradigm for GWAS analysis. By evaluating genetic variation across a biological pathway (gene set), these tests have the potential to determine the collective impact of variants with subtle effects that are individually too weak to be detected in traditional single variant GWAS analysis. To identify biological pathways associated with the risk of development of PBC, GWAS of PBC from Italy (449 cases and 940 controls) and Canada (530 cases and 398 controls) were independently analyzed. The linear combination test (LCT), a recently developed pathway-level statistical method was used for this analysis. For additional validation, pathways that were replicated at the P <0.05 level of significance in both GWAS on LCT analysis were also tested for association with PBC in each dataset using two complementary GWAS pathway approaches. The complementary approaches included a modification of the gene set enrichment analysis algorithm (i-GSEA4GWAS) and Fisher's exact test for pathway enrichment ratios. Twenty-five pathways were associated with PBC risk on LCT analysis in the Italian dataset at P<0.05, of which eight had an FDR<0.25. The top pathway in the Italian dataset was the TNF/stress related signaling pathway (p=7.38×10 -4, FDR=0.18). Twenty-six pathways were associated with PBC at the P<0.05 level using the LCT in the Canadian dataset with the regulation and function of ChREBP in liver pathway (p=5.68×10-4, FDR=0.285) emerging as the most significant pathway. Two pathways, phosphatidylinositol signaling system (Italian: p=0.016, FDR=0.436; Canadian: p=0.034, FDR=0.693) and hedgehog signaling (Italian: p=0.044, FDR=0.636; Canadian: p=0.041, FDR=0.693), were replicated at LCT P<0.05 in both datasets. Statistically significant association of both pathways with PBC genetic susceptibility was confirmed in the Italian dataset on i-GSEA4GWAS. Results for the phosphatidylinositol signaling system were also significant in both datasets on applying Fisher's exact test for pathway enrichment ratios. This study identified a combination of known and novel pathway-level associations with PBC risk. If functionally validated, the findings may yield fresh insights into the etiology of this complex autoimmune disease with possible preventive and therapeutic application.^

Identification of genetic risk factors for cerebellar mutism in pediatric brain tumor patients

Following posterior fossa surgery for resection of childhood medulloblastoma and primitive neuroectodermal tumor (M/PNET), cerebellar mutism (CM) may develop. This is a condition of absent or diminished speech in a conscious patient with no evidence of oral apraxia, which can be accompanied by other symptoms of the posterior fossa syndrome complex, which includes ataxia and hypotonia. Little is known about the etiology. Therefore, we conducted a SNP, gene, and pathway-level analysis to assess the role of host genetic variation on the risk of CM in M/PNET subjects following treatment. Cases (n= 20) and controls (n= 53) were recruited from the Childhood Cancer Epidemiology and Prevention Center, in Houston, TX. DNA samples were genotyped using the Illumina Human 1M Quad SNP chip. Ten pathways were identified from logistic regression used to identify the marginal effect of each SNP on CM risk. The minP test was conducted to identify associations between SNPs categorized to genes and CM risk. Pathways were assessed to determine if there was a significant enrichment of genes in the pathway compared to all other pathways. There were 78 genes that reached the threshold of min P ≤0.05 in 948 genes. The Neurotoxicity pathway was the most significant pathway after adjusting for multiple comparisons (q=0.040 and q=0.005, using Fisher's exact test and a test of proportions, respectively). Most genes within the Neurotoxicity pathway that reached a threshold of minP ≤0.05 were known to have an apoptosis function, possibly inducing neuronal apoptosis in the dentatothalamocortical pathway, and may be important in CM etiology in this population. This is the first study to assess the potential role of genetic risk factors on CM. As an exploratory study, these results should be replicated in a larger sample. ^

A molecular study of the rhodopsin locus in patients with retinitis pigmentosa

DNA for this study was collected from a sample of 133 retinitis pigmentosa (RP) patients and the rhodopsin locus molecularly analyzed by linkage and for disease specific mutations. The cohort of patients consisted of 85 individuals diagnosed with autosomal dominant RP (adRP), and 48 patients representing other forms of retinitis pigmentosa or retinal dystrophy related disease. In three large families with adRP rhodopsin was excluded from linkage to the disease locus. A search for subtle mutations in the rhodopsin coding region using single strand conformational polymorphisms (SSCP) and sequencing detected a total of 14 unique sequence variants in 24 unrelated patients. These variants included one splicing variant, 5168 -1G-A, one deletion variant of 17 base pairs causing a frame shift at codon 332, and 12 misense variants: Pro23His, Leu46Arg, Gly106Trp, Arg135Pro, Pro171Glu, Pro180Ala, Glu181Lys, Asp190Asn, His211Arg, Ser270Arg, Leu328Pro and Pro347Thr. All but three of the missense variants change amino acids that are evolutionarily conserved. The Pro23His mutation was found in 10 unrelated individuals with family histories of adRP and not in any normal controls (over 80 chromosomes tested). The Pro180Ala mutation was present in a patient with simplex RP and probably represents a new mutation. Three normal polymorphic nucleotide substitutions, A-269-G, T-3982-C, and G-5145-A, were also identified. We conclude, based on this study, that 25% of adRP cases are attributable to rhodopsin mutations.^ Clinical data, including ERG results and visual field testing, was available for patients with eleven different mutations. The eleven patients were all diagnosed with RP, however the severity of the disease varied with five patients mildly affected and diagnosed with type II adRP and 5 patients severely affected and diagnosed with type I adRP. The patient with simplex RP was mildly affected. The location of the mutations within the rhodopsin protein was randomly associated with the severity of the disease in those patients evaluated. However, four mutations, Pro23His, Leu46Arg, Pro347Thr, and 5168 -1G-A, are particularly interesting. The Pro23His mutation appears to have radiated from a recent common ancestor of the affected patients as all of them share a common haplotype at the rhodopsin locus. The Leu46Arg mutation causes an unusually severe form of RP. Hydropathy analysis of the mutated sequence revealed a marked change in the hydrophobicity of this first transmembrane spanning region. Codon 347 has been the target of multiple mutations with at least six documented changes at the position, significantly more than expected by a random distribution of mutations. Finally the splice-site variant is extremely variable in its expression in the family studied. Similar mutations have been reported in other cases of adRP and postulated to be involved in autosomal recessive RP (arRP). Mechanisms to account for the variable expression of rhodopsin mutations in relation to RP heterogeneity are discussed. (Abstract shortened by UMI.) ^

Structure-function analyses of {\it PAX6\/} and the molecular bases of aniridia

PAX6 is a transcription activator that regulates eye development in animals ranging from Drosophila to human. The C-terminal region of PAX6 is proline/serine/threonine-rich (PST) and functions as a potent transactivation domain when attached to a heterologous DNA-binding domain of the yeast transcription factor, GAL4. The PST region comprises 152 amino acids encoded by four exons. The transactivation function of the PST region has not been defined and characterized in detail by in vitro mutagenesis. I dissected the PST domain in two independent systems, a heterologous system using a GAL4 DNA-binding site and the native system of PAX6. In both systems, the results show consistently that all four constituent exons of the PST domain are responsible for the transactivation function. The four exon fragments act cooperatively to stimulate transcription, although none of them can function individually as an independent transactivation domain. Combinations of two or more exon fragments can reconstitute substantial transactivation activity when fused to the DNA-binding domain of GAL4, but they surprisingly do not produce much activity in the context of native PAX6 even though the mutant PAX6 proteins are stable and their DNA-binding function remains unaffected. I conclude that the PAX6 protein contains an unusually large transactivation domain that is evolutionarily conserved to a high degree, and that its full transactivation activity relies on the cooperative action of the four exon fragments.^ Most PAX6 mutations detected in patients with aniridia result in truncations of the protein. Some of the truncation mutations occur in the PST region of PAX6, resulting in mutant proteins that retain their DNA-binding ability but have no significant transactivation activity. It is not clear whether such mutants are true loss-of-function or dominant-negative mutants. I show that these mutants are dominant-negative if they are coexpressed with wild-type PAX6 in cultured cells and that the dominant-negative effects result from enhanced DNA-binding ability of these mutants due to removal of the PST domain. These mutants are able to repress the wild-type PAX6 activity not only at target genes with paired domain binding sites but also at target genes with homeodomain binding sites.^ Mutations in the human PAX6 gene produce various phenotypes, including aniridia, Peters' anomaly, autosomal dominant keratitis, and familial foveal dysplasia. The various phenotypes may arise from different mutations in the same gene. To test this theory, I performed a functional analysis of two missense mutations in the paired domain: the R26G mutation reported in a case of Peters' anomaly, and the I87R mutation identified in a patient with aniridia. While both the R26 and the I87 positions are conserved in the paired boxes of all known PAX genes, X-ray crystallography has shown that only R26 makes contact with DNA. I found that the R26G mutant failed to bind a subset of paired domain binding sites but, surprisingly, bound other sites and successfully transactivated promoters containing those sites. In contrast, the I87R mutant had lost the ability to bind DNA at all tested sites and failed to transactivate promoters. My data support the haploinsufficiency hypothesis of aniridia, and the hypothesis that R26G is a hypomorphic allele. ^

Population and molecular evolution studies on the Melanocortin 1 receptor locus implicate its role in human pigmentation variation

Human pigmentation is a complex trait with the observed variation caused by the varied production of eumelanin (brown/black melanins) and phaeomelanin (red/yellow melanins) by the melanocytes. The melanocortin 1 receptor (MC1R), a G protein-coupled receptor expressed in the melanocytes, is a regulator eu- and phaeomelanin synthesis, and MC1R mutations causing skin and coat color changes are known in many mammals. To understand the role of MC1R in human pigmentation variation, I have sequenced the MC1R gene in 121 individuals sampled from world populations. In addition, I have sequenced the MC1R gene in common and pygmy chimpanzees, gorilla, orangutan, and baboon to study the evolution of MC1R and to infer the ancestral human MC1R sequence. The ancestral MC1R sequence is observed in all 25 African individuals studied, but at lower frequencies in the other populations examined, especially in East and Southeast Asians. The Arg163Gln variant is absent in the Africans studied, almost absent in Europeans, and at a low frequency in Indians, but is at an exceptionally high frequency (70%) in East and Southeast Asians. To further evaluate the role of MC1R variants in human pigmentation variation, I have combined these molecular evolution and population studies with functional assays on MC1R variants and primate MC1Rs. ^