Identification of candidate genes for prostate cancer-risk SNPs utilizing a normal prostate tissue eQTL data set.

Multiple studies have identified loci associated with the risk of developing prostate cancer but the associated genes are not well studied. Here we create a normal prostate tissue-specific eQTL data set and apply this data set to previously identified prostate cancer (PrCa)-risk SNPs in an effort to identify candidate target genes. The eQTL data set is constructed by the genotyping and RNA sequencing of 471 samples. We focus on 146 PrCa-risk SNPs, including all SNPs in linkage disequilibrium with each risk SNP, resulting in 100 unique risk intervals. We analyse cis-acting associations where the transcript is located within 2 Mb (±1 Mb) of the risk SNP interval. Of all SNP-gene combinations tested, 41.7% of SNPs demonstrate a significant eQTL signal after adjustment for sample histology and 14 expression principal component covariates. Of the 100 PrCa-risk intervals, 51 have a significant eQTL signal and these are associated with 88 genes. This study provides a rich resource to study biological mechanisms underlying genetic risk to PrCa.

Pharmacogenomics of selective serotonin reuptake inhibitor treatment for major depressive disorder: genome-wide associations and functional genomics.

A genome-wide association (GWA) study of treatment outcomes (response and remission) of selective serotonin reuptake inhibitors (SSRIs) was conducted using 529 subjects with major depressive disorder. While no SNP associations reached the genome-wide level of significance, 14 SNPs of interest were identified for functional analysis. The rs11144870 SNP in the riboflavin kinase (RFK) gene on chromosome 9 was associated with 8-week treatment response (odds ratio (OR)=0.42, P=1.04 × 10⁻6). The rs915120 SNP in the G protein-coupled receptor kinase 5 (GRK5) gene on chromosome 10 was associated with 8-week remission (OR=0.50, P=1.15 × 10⁻5). Both SNPs were shown to influence transcription by a reporter gene assay and to alter nuclear protein binding using an electrophoretic mobility shift assay. This report represents an example of joining functional genomics with traditional GWA study results derived from a GWA analysis of SSRI treatment outcomes. The goal of this analytical strategy is to provide insights into the potential relevance of biologically plausible observed associations.

Glycine and a glycine dehydrogenase (GLDC) SNP as citalopram/escitalopram response biomarkers in depression: pharmacometabolomics-informed pharmacogenomics.

Major depressive disorder (MDD) is a common psychiatric disease. Selective serotonin reuptake inhibitors (SSRIs) are an important class of drugs used in the treatment of MDD. However, many patients do not respond adequately to SSRI therapy. We used a pharmacometabolomics-informed pharmacogenomic research strategy to identify citalopram/escitalopram treatment outcome biomarkers. Metabolomic assay of plasma samples from 20 escitalopram remitters and 20 nonremitters showed that glycine was negatively associated with treatment outcome (P = 0.0054). This observation was pursued by genotyping tag single-nucleotide polymorphisms (SNPs) for genes encoding glycine synthesis and degradation enzymes, using 529 DNA samples from SSRI-treated MDD patients. The rs10975641 SNP in the glycine dehydrogenase (GLDC) gene was associated with treatment outcome phenotypes. Genotyping for rs10975641 was carried out in 1,245 MDD patients in the Sequenced Treatment Alternatives to Relieve Depression (STAR*D) study, and its presence was significant (P = 0.02) in DNA taken from these patients. These results highlight a possible role for glycine in SSRI response and illustrate the use of pharmacometabolomics to "inform" pharmacogenomics.

SLC6A4 variation and citalopram response.

The influence of genetic variations in SLC6A4 (serotonin transporter gene) on citalopram treatment of depression using the Sequenced Treatment to Relieve Depression (STAR*D) sample was assessed. Of primary interest were three previously studied polymorphisms: 1) the VNTR variation of the second intron, 2) the indel promoter polymorphism (5HTTLPR or SERT), and 3) a single nucleotide polymorphism (SNP) rs25531. Additionally, SLC6A4 was resequenced to identify new SNPs for exploratory analyses. DNA from 1914 subjects in the STAR*D study were genotyped for the intron 2 VNTR region, the indel promoter polymorphism, and rs25531. Associations of these variants with remission of depressive symptoms were evaluated following citalopram treatment. In white non-Hispanic subjects, variations in the intron 2 VNTR (point-wise P = 0.041) and the indel promoter polymorphism (point-wise P = 0.039) were associated with remission following treatment with citalopram. The haplotype composed of the three candidate loci was also associated with remission, with a global p-value of 0.040 and a maximum statistic simulation p-value of 0.0031 for the S-a-12 haplotype, under a dominant model. One SNP identified through re-sequencing the SLC6A4 gene, Intron7-83-TC, showed point-wise evidence of association, which did not remain significant after correction for the number of SNPs evaluated in this exploratory analysis. No associations between these SLC6A4 variations and remission were found in the white Hispanic or black subjects. These findings suggest that multiple variations in the SLC6A4 gene are associated with remission in white non-Hispanic depressed adults treated with citalopram. The mechanism of action of these variants remains to be determined.

Genetic diversity and function in the human cytosolic sulfotransferases.

Amino-acid substitutions, which result from common nonsynonymous (NS) polymorphisms, may dramatically alter the function of the encoded protein. Gaining insight into how these substitutions alter function is a step toward acquiring predictability. In this study, we incorporated gene resequencing, functional genomics, amino-acid characterization and crystal structure analysis for the cytosolic sulfotransferases (SULTs) to attempt to gain predictability regarding the function of variant allozymes. Previously, four SULT genes were resequenced in 118 DNA samples. With additional resequencing of the remaining eight SULT family members in the same DNA samples, a total of 217 polymorphisms were revealed. Of 64 polymorphisms identified within 8785 bp of coding regions from SULT genes examined, 25 were synonymous and 39 were NS. Overall, the proportion of synonymous changes was greater than expected from a random distribution of mutations, suggesting the presence of a selective pressure against amino-acid substitutions. Functional data for common variants of five SULT genes have been previously published. These data, together with the SULT1A1 variant allozyme data presented in this paper, showed that the major mechanism by which amino acid changes altered function in a transient expression system was through decreases in immunoreactive protein rather than changes in enzyme kinetics. Additional insight with regard to mechanisms by which NS single nucleotide polymorphisms alter function was sought by analysis of evolutionary conservation, physicochemical properties of the amino-acid substitutions and crystal structure analysis. Neither individual amino-acid characteristics nor structural models were able to accurately and reliably predict the function of variant allozymes. These results suggest that common amino-acid substitutions may not dramatically alter the protein structure, but affect interactions with the cellular environment that are currently not well understood.

Estimation and tests of haplotype-environment interaction when linkage phase is ambiguous.

In the study of complex traits, the utility of linkage analysis and single marker association tests can be limited for researchers attempting to elucidate the complex interplay between a gene and environmental covariates. For these purposes, tests of gene-environment interactions are needed. In addition, recent studies have indicated that haplotypes, which are specific combinations of nucleotides on the same chromosome, may be more suitable as the unit of analysis for statistical tests than single genetic markers. The difficulty with this approach is that, in standard laboratory genotyping, haplotypes are often not directly observable. Instead, unphased marker phenotypes are collected. In this article, we present a method for estimating and testing haplotype-environment interactions when linkage phase is potentially ambiguous. The method builds on the work of Schaid et al. [2002] and is applicable to any trait that can be placed in the generalized linear model framework. Simulations were run to illustrate the salient features of the method. In addition, the method was used to test for haplotype-smoking exposure interaction with data from the Childhood Asthma Management Program.

Validity of family history data on PD: evidence for a family information bias.

OBJECTIVE: To study the validity of information provided by case and control subjects (or their proxies) about PD among their first-degree relatives. METHODS: Secondary cases of PD were assessed both through a single informant (family history method) and through the study of each relative (family study method). The family study method was considered as the standard for comparison, and the sensitivity and specificity of the family history method were studied. RESULTS: A total of 133 population-based case subjects and their 655 relatives were recruited, and 119 population-based control subjects and their 511 relatives. Sensitivity was 68% (95% CI = 47 to 85) for cases and 45% (95% CI = 17 to 77) for controls. Specificity was 99% (95% CI = 98 to 99) for cases and 100% (95% CI = 99 to 100) for controls. The odds ratio (OR) for family history of PD was 4.34 (95% CI = 1.63 to 11.58, p = 0.003) using the family history method and 1.86 (95% CI = 0.78 to 4.44, p = 0.16) using the family study method. The former significant OR more than doubled the latter not significant OR (relative bias = 133%). Bias was more pronounced for proxy interviews and for women informants, and when the relatives were siblings, were living, and were examined or had medical record documentation. CONCLUSIONS: Case subjects with PD (or their proxies) are more aware of PD among their first-degree relatives than control subjects (or their proxies); however, they overreport PD in relatives who are not affected. This causes a substantial family information bias.

Progression of familial and non-familial dilated cardiomyopathy: long term follow up.

BACKGROUND: It is unknown whether progression of familial idiopathic dilated cardiomyopathy differs from progression in the non-familial form. It has been suggested that familial disease indicates a worse prognosis, and that this should be considered when planning the timing of heart transplantation. OBJECTIVE: To compare five year survival or time to heart transplantation in an unselected series of patients with dilated cardiomyopathy who had been evaluated for familial v non-familial disease through the echocardiographic investigation of first degree relatives. DESIGN: Medical records were reviewed and questionnaires were mailed to all patients who had previously participated in a family based study of dilated cardiomyopathy. Information was gathered about survival, heart transplantation, and left ventricular ejection fraction (LVEF) measurements. Survival data were censored at the time of cardiac transplantation. RESULTS: Follow up data were obtained for 99 of 101 patients (69 with non-familial and 30 with familial disease). Five year survival was 55% for non-familial and 51% for familial patients (NS). The main predictor of mortality was an LVEF of < 30%. Familial status did not predict mortality. There was no significant difference in follow up LVEF values between the groups. CONCLUSIONS: Five year survival is not significantly different in the familial and non-familial forms of dilated cardiomyopathy.

Candidate-gene association studies with pedigree data: controlling for environmental covariates.

Case-control studies provide an important epidemiological tool to evaluate candidate genes. There are many different study designs available. We focus on a more recently proposed design, which we call a multiplex case-control (MCC) design. This design compares allele frequencies between related cases, each of whom are sampled from multiplex families, and unrelated controls. Since within-family genotype correlations will exist, statistical methods will need to take this into account. Moreover, there is a need to develop methods to simultaneously control for potential confounders in the analysis. Generalized estimating equations (GEE) are one approach to analyze this type of data; however, this approach can have singularity problems when estimating the correlation matrix. To allow for modeling of other covariates, we extend our previously developed method to a more general model-based approach. Our proposed methods use the score statistic, derived from a composite likelihood. We propose three different approaches to estimate the variance of this statistic. Under random ascertainment of pedigrees, score tests have correct type I error rates; however, pedigrees are not randomly ascertained. Thus, through simulations, we test the validity and power of the score tests under different ascertainment schemes, and an illustration of our methods, applied to data from a prostate cancer study, is presented. We find that our robust score statistic has estimated type I error rates within the expected range for all situations we considered whereas the other two statistics have inflated type I error rates under nonrandom ascertainment schemes. We also find GEE to fail at least 5% of the time for each simulation configuration; at times, the failure rate reaches above 80%. In summary, our robust method may be the only current regression analysis method available for MCC data.

The gene for HMSN2C maps to 12q23-24: a region of neuromuscular disorders.

BACKGROUND: Hereditary motor and sensory neuropathy type 2C (HMSN2C, Charcot-Marie-Tooth 2C [CMT2C]) is an autosomal dominant motor and sensory neuropathy involving limb, diaphragm, vocal cord, and intercostal muscles. OBJECTIVE: To identify the chromosome localization for this disorder in one large American family of English and Scottish ethnicity. METHODS: Variable clinical severity led the authors to combine several approaches to accurately identify affected patients. Genome-wide two-point linkage analysis, high-definition mapping, and multipoint and recombinant haplotype analyses were performed. Mutation analysis of the triplet repeat region of ataxin-2 was also carried out. RESULTS: The initial genome-wide scan identified a region at 12q24, and fine mapping provided a maximal lod score of 4.73 (D12S1645 and D12S1583 at theta = 0.01 and 0, respectively). With multipoint analysis, a higher lod score of 5.17 was obtained and localized to the same region at 119.0 cM. Haplotype analysis narrowed the region to approximately 5.0 cM between D12S1646,D12S1330 and D12S105,D12S1339 (12q23.3-24.21). Ataxin-2, the gene responsible for spinocerebellar ataxia type 2 (SCA2), localizes to this region, but no triplet repeat expansion or point mutations within the repeat were found. CONCLUSIONS: The gene for HMSN2C maps to 12q23-24. This region is associated with SCA2, scapuloperoneal spinal muscular atrophy, and congenital distal spinal muscular atrophy. Further studies are needed to demonstrate the specific gene alteration and its relationship with nearby genes.

Confirmation of linkage of prostate cancer aggressiveness with chromosome 19q.

Regions on chromosomes 7 and 19 were recently reported to contain susceptibility loci that regulate tumor aggressiveness of prostate cancer. To confirm these findings, we analyzed genome scan data from 161 pedigrees affected with prostate cancer. Using the Gleason score as a quantitative measure of tumor aggressiveness, we regressed the squared trait difference, as well as the mean-corrected cross product, on the estimated proportion of alleles shared identical-by-descent at each marker position. Our results confirm the previous linkage results for chromosome 19q (D19S902, P<.00001). In addition, we report suggestive evidence for linkage on chromosome 4 (D4S403, P=.00012). The results of previous findings, together with our results, provide strong evidence that chromosome 19 harbors a gene for tumor aggressiveness.

Hysterectomy, menopause, and estrogen use preceding Parkinson's disease: an exploratory case-control study.

We studied the association of Parkinson's disease (PD) with type of menopause (natural or surgical), age at menopause, and postmenopausal estrogen replacement therapy using a case-control design. We used the medical records-linkage system of the Rochester Epidemiology Project to identify 72 women who developed PD in Olmsted County, MN, during the twenty years 1976-1995. Each incident case was matched by age (+/- 1 year) to a general population control subject. We collected exposure data through review of the complete medical records of cases and control subjects in the system. PD cases had undergone hysterectomy (with or without unilateral oophorectomy) significantly more often than control subjects (odds ratio [OR] = 3.36; 95% confidence interval [CI] = 1.05-10.77). In addition, PD cases had experienced early menopause (< or = 46 years) more commonly than control subjects (OR = 2.18; 95% CI = 0.88-5.39). Finally, PD cases had used estrogens orally or parenterally for at least 6 months after menopause less frequently (8%) than control subjects (14%; OR = 0.47; 95% CI = 0.12-1.85). However, the findings for early menopause and estrogen replacement therapy were not statistically significant. Despite the limited sample size of this exploratory study, we hypothesize that there is an increased risk of PD in conditions causing an early reduction in endogenous estrogen. This hypothesis needs to be confirmed in a larger study.

Incidence of reading disability in a population-based birth cohort, 1976-1982, Rochester, Minn.

OBJECTIVE: To report the incidence of reading disability among school-aged children. SUBJECTS AND METHODS: In this population-based, retrospective birth cohort study, subjects included all 5718 children born between 1976 and 1982 who remained in Rochester, Minn, after the age of 5 years. Based on records from all public and nonpublic schools, medical facilities, and private tutorial services and on results of all individually administered IQ and achievement tests, extensive medical, educational, and socioeconomic information were abstracted. Reading disability was established with use of research criteria based on 4 formulas (2 regression-based discrepancy, 1 non-regression-based discrepancy, and 1 low achievement). RESULTS: Cumulative incidence rates of reading disability varied from 5.3% to 11.8% depending on the formula used. Boys were 2 to 3 times more likely to be affected than girls, regardless of the identification methods applied. CONCLUSIONS: In this population-based birth cohort, reading disability was common among school-aged children and significantly more frequent among boys than girls, regardless of definition.

Efficacy of bilateral prophylactic mastectomy in BRCA1 and BRCA2 gene mutation carriers.

BACKGROUND: In women with a family history of breast cancer, bilateral prophylactic mastectomy is associated with a decreased risk of subsequent breast cancer of approximately 90%. We examined the association between bilateral prophylactic mastectomy and breast cancer risk in women at high risk for breast cancer who also had mutations in BRCA1 and BRCA2 genes. METHODS: We obtained blood samples from 176 of the 214 high-risk women who participated in our previous retrospective cohort study of bilateral prophylactic mastectomy. We used conformation-sensitive gel electrophoresis and direct sequence analysis of the blood specimens to identify women with mutations in BRCA1 and BRCA2. The carriers' probabilities of developing breast cancer were estimated from two different penetrance models. RESULTS: We identified 26 women with an alteration in BRCA1 or BRCA2. Eighteen of the mutations were considered to be deleterious and eight to be of uncertain clinical significance. None of the 26 women has developed breast cancer after a median of 13.4 years of follow-up (range, 5.8-28.5 years). Three of the 214 women are known to have developed a breast cancer after prophylactic mastectomy. For two of these women, BRCA1 and BRCA2 screening was negative, and no blood specimen was available for the third. Estimations of the effectiveness of prophylactic mastectomy were performed, considering this woman as both a mutation carrier and a noncarrier. These calculations predicted that six to nine breast cancers should have developed among the mutation carriers, which translates into a risk reduction, after bilateral prophylactic mastectomy, of 89.5%-100% (95% confidence interval = 41.4% to 100%). CONCLUSIONS: Prophylactic mastectomy is associated with a substantial reduction in the incidence of subsequent breast cancer not only in women identified as being at high risk on the basis of a family history of breast cancer but also in known BRCA1 or BRCA2 mutation carriers.

Genetic variation in the transforming growth factor beta1 gene in multiple sclerosis.

Transforming growth factor beta1 (TGFbeta1) is a Th2 cytokine encoded on chromosome 19q13, a region possibly linked to multiple sclerosis (MS). TGFbeta1 exerts favorable effects on experimental allergic encephalomyelitis. We performed a comprehensive search for genetic variants in this gene in 122 population-based sporadic cases of MS. We detected six variants, including three missense variants. We tested for association of the variants with susceptibility and course of MS and for linkage and transmission disequilibrium in a family series consisting of 395 samples in 59 pedigrees. Genetic variation in TGFB1 does not appear to contribute in a major way to susceptibility to MS.

Efficacy of contralateral prophylactic mastectomy in women with a personal and family history of breast cancer.

PURPOSE: To estimate the efficacy of contralateral prophylactic mastectomy in women with a personal and family history of breast cancer. PATIENTS AND METHODS: We followed the course of 745 women with a first breast cancer and a family history of breast and/or ovarian cancer who underwent contralateral prophylactic mastectomy at the Mayo Clinic between 1960 and 1993. Family history information and cancer follow-up information were obtained from the medical record, a study-specific questionnaire, and telephone follow-up. Life-tables for contralateral breast cancers, which consider age at first breast cancer, current age, and type of family history, were used to calculate the number of breast cancers expected in our cohort had they not had a prophylactic mastectomy. RESULTS: Of the 745 women in our cohort, 388 were premenopausal (age < 50 years) and 357 were post- menopausal. Eight women developed a contralateral breast cancer. Six events were observed among the premenopausal women, compared with 106.2 predicted, resulting in a risk reduction of 94.4% (95% confidence interval [CI], 87.7% to 97.9%). For the 357 postmenopausal women, 50.3 contralateral breast cancers were predicted, whereas only two were observed, representing a 96.0% risk reduction (95% CI, 85.6% to 99.5%). CONCLUSION: The incidence of contralateral breast cancer seems to be reduced significantly after contralateral prophylactic mastectomy in women with a personal and family history of breast cancer.

Identification of an association between HLA class II alleles and low antibody levels after measles immunization.

This is the first large cohort study to report a genetic association between humoral antibody level after measles vaccine and the HLA class II genes. The WHO goal to eradicate measles world-wide magnifies the importance of data relating to the influence of immunogenetics on measles vaccine-induced antibody responses. We present here the analysis of 242 individuals who received one dose of measles-mumps-rubella-II (MMR-II) vaccine at the age of 15 months and were genotyped for HLA class II alleles. These subjects fit into one of three categories; 72 were classified as seronegative, 93 were seropositive and 77 were serohyperpositive. HLA-DRB1*03 (odds ratio (OR), 2.22) and HLA-DPA1*0201 (OR, 1.71) were significantly associated with measles vaccine seronegativity, while additional alleles provided suggestive evidence of association with seronegativity: DQA1*0201, DQB1*0201, and DQA1*0501. The alleles DRB1*03 and DQA1*0201 remained statistically significant after accounting for the effects of other alleles. These findings are crucial in designing both measles eradication by the use of vaccine, and future vaccines to be used in genetically heterozygous populations.

Case-control studies of genetic markers: power and sample size approximations for Armitage's test for trend.

The association of a candidate gene with disease can be efficiently evaluated by a case-control study in which allele frequencies are compared for diseased cases and unaffected controls. However, when the distribution of genotypes in the population deviates from Hardy-Weinberg proportions, the frequency of genotypes--rather than alleles--should be compared by the Armitage test for trend. We present formulas for power and sample size for studies that use Armitage's trend test. The formulas make no assumptions about Hardy-Weinberg equilibrium, but do assume random ascertainment of cases and controls, all of whom are independent of one another. We demonstrate the accuracy of the formulas by simulations.

The frequency of hereditary defective mismatch repair in a prospective series of unselected colorectal carcinomas.

A comprehensive analysis of somatic and germline mutations related to DNA mismatch-repair (MMR) genes can clarify the prevalence and mechanism of inactivation in colorectal carcinoma (CRC). In the present study, 257 unselected patients referred for CRC resection were examined for evidence of defective DNA MMR. In particular, we sought to determine the frequency of hereditary defects in DNA MMR in this cohort of patients. MMR status was assessed by testing of tumors for the presence or absence of hMLH1, hMSH2, and hMSH6 protein expression and for microsatellite instability (MSI). Of the 257 patients, 51 (20%) had evidence of defective MMR, demonstrating high levels of MSI (MSI-H) and an absence of either hMLH1 (n=48) or hMSH2 (n=3). All three patients lacking hMSH2, as well as one patient lacking hMLH1, also demonstrated an absence of hMSH6. DNA sequence analysis of the 51 patients with defective MMR revealed seven germline mutations-four in hMLH1 (two truncating and two missense) and three in hMSH2 (all truncating). A detailed family history was available for 225 of the 257 patients. Of the seven patients with germline mutations, only three had family histories consistent with hereditary nonpolyposis colorectal cancer. Of the remaining patients who had tumors with defective MMR, eight had somatic mutations in hMLH1. In addition, hypermethylation of the hMLH1 gene promoter was present in 37 (88%) of the 42 hMLH1-negative cases available for study and in all MSI-H tumors that showed loss of hMLH1 expression but no detectable hMLH1 mutations. Our results suggest that, although defective DNA MMR occurs in approximately 20% of unselected patients presenting for CRC resection, hereditary CRC due to mutations in the MMR pathway account for only a small proportion of patients. Of the 257 patients, only 5 (1.9%) appear to have unequivocal evidence of hereditary defects in MMR. The epigenetic (nonhereditary) mechanism of hMLH1 promoter hypermethylation appears to be responsible for the majority of the remaining patients whose tumors are characterized by defective DNA MMR.