- Split View
-
Views
-
Cite
Cite
Graham Casey, Phillippa J. Neville, Xin Liu, Sarah J. Plummer, Mine S. Cicek, Lisa M. Krumroy, Anthony P. Curran, Michelle R. McGreevy, William J. Catalona, Eric A. Klein, John S. Witte, Podocalyxin variants and risk of prostate cancer and tumor aggressiveness, Human Molecular Genetics, Volume 15, Issue 5, 1 March 2006, Pages 735–741, https://doi.org/10.1093/hmg/ddi487
- Share Icon Share
Abstract
We previously reported linkage of a prostate cancer tumor aggressiveness locus to chromosome 7q32–q33, a region also associated with a high frequency of allelic imbalance in prostate tumors. The smallest region of allelic imbalance contains the podocalyxin-like (PODXL) gene, which we evaluate here as a candidate prostate cancer aggressiveness gene mapping to 7q32–q33. DNA from probands of linked families was examined for germ-line mutations in PODXL. A variable in-frame deletion, four missense variants and two nonsense variants were identified in linked men. Variants that affected amino acid sequence were further evaluated for association with risk of prostate cancer and tumor aggressiveness in a family-based case–control population (439 cases and 479 sibling controls). The presence of any single in-frame deletion was positively associated with prostate cancer [odds ratio (OR)=2.14, 95% confidence interval (95%CI)=1.09–4.20, P=0.03] and the presence of two copies of any deletion further increased risk (OR=2.58, 95%CI=1.23–5.45, P=0.01). This finding was strengthened when stratifying among men with more aggressive disease (high grade or stage): OR=3.04 for one deletion (95%CI=1.01–9.15) and OR=4.42 for two deletions (95%CI=1.32–14.85, P=0.02). A weak positive association was also observed between prostate cancer risk and PODXL variant 340A (in linkage disequilibrium with another variant, 587T) (OR=1.48, 95%CI=1.02–2.14, P=0.04). These results implicate PODXL as a candidate prostate cancer tumor aggressiveness gene mapping to chromosome 7q32–q33.
INTRODUCTION
The aggressiveness of prostate cancer varies widely among patients, with some cancers rapidly progressing to invasive, life-threatening disease, whereas others remain indolent for many years. The mechanisms underlying this variability in disease presentation and course are not known, but there is evidence that genetics may play a role, as prostate cancer tumor aggressiveness is influenced by a family history of this disease (1–4).
To identify genes that affect the aggressiveness of prostate cancer, we previously performed a genome-wide scan of 594 affected siblings, using Gleason score as a surrogate for tumor aggressiveness, and reported strong linkage to chromosome 7q32–q33 (5,6). Linkage to this region has been replicated in two independent studies (7,8) and a third report has linked a neighboring region on chromosome 7q to prostate cancer susceptibility among 36 Jewish families (9). Chromosome 7q32–q33 has also been shown to exhibit a high frequency of allelic imbalance in prostate tumors, with the smallest region of allelic imbalance spanning approximately 1.1 Mb and corresponding to the peak linkage marker from our affected sibling pair study (5,6,8).
The 1.1 Mb region on chromosome 7q32–q33 contains the podocalyxin-like (PODXL) gene, which encodes a 528 amino acid protein. PODXL, also known as podocalyxin, podocalyxin-like protein 1 (PCLP-1), Myb-Ets-transformed progenitor 21 (MEP21) and thrombomucin, is a transmembrane sialomucin structurally related to the vascular endothelium marker CD34 (10). Studies reveal two diverse roles for the podocalyxin protein. In podocytes of the renal glomerulae, podocalyxin functions as an anti-adhesion molecule that maintains an open filtration pathway between neighboring podocyte foot processes by charge repulsion (11). In the vasculature, podocalyxin acts as an adhesion molecule on high endothelial venules, supporting the tethering and rolling of lymphocytes through binding to L-selectin (12).
Podocalyxin has been implicated in the development of cancer and, in particular, more aggressive cancer. Podocalyxin is a downstream target gene of the WT1 tumor suppressor (13,14) and has been implicated in beta catenin signaling (15). Studies have also shown that podocalyxin expression is negatively regulated by TP53 in kidney cells (16). Recent data have shown a correlation between high levels of podocalyxin protein expression and a more aggressive breast cancer (17). Furthermore, gene expression studies have implicated podocalyxin in angiogenesis, where a positive correlation was seen between expression of PODXL and ETS1, a regulator of angiogenesis, in human vascular endothelial cells (18). Additionally, podocalyxin was recently reported commonly expressed by blasts in both acute myeloid and acute lymphoblastic leukemias (19).
The linkage and allelic imbalance studies noted previously, along with the known functions of podocalyxin, make it an attractive candidate gene for the involvement in tumor aggressiveness. To assess the potential role of podocalyxin in prostate cancer risk, we first performed PODXL germ-line mutation analyses of probands from families linked to chromosome 7q32–q33 (5). We then evaluated the relationship between these PODXL variants and prostate cancer and tumor aggressiveness in a family-based association study. A variable in-frame deletion and a missense variant of PODXL were associated with risk of prostate cancer and tumor aggressiveness.
RESULTS
PODXL variants in families linked to chromosome 7q32–q33
Mutation analyses of PODXL were performed by DNA sequencing, using genomic DNA from white blood cells of probands of the 17 families that showed linkage to chromosome 7q32–q33. Primers were designed to amplify all the eight exons of PODXL including the intron splice donor and acceptor regions (Table 1). A number of coding sequence variants in PODXL were identified (Table 2). Specifically, we identified a common variable in-frame deletion, four common missense variants, two silent variants and two intronic substitutions (g.IVS4-51 T>C and g.IVS6+6 C>T). The variable in-frame deletion in exon 1 was 6 (TCGCCG) or 12 bp (TCGCCGTCGCCG) nucleotides in length (herein termed del6 and del12) and resulted in the loss of one or two serine–proline repeats (two or four amino acid residues, respectively). The in-frame deletion variant and missense variants 340G>A, 587C>T and 802C>G were genotyped in the case–control population. A reliable assay for the 982G>A variant could not be developed and was not investigated further.
Allele frequencies for PODXL variants in a family-based case–control population
Allele frequencies of the PODXL variants in the case–control population are shown in Table 3. Three cases were not successfully genotyped in the case–control population and were not included in the results. Note that a third form of the in-frame deletion, a deletion of 18 bp (TCGCCGTCGCCGTCGCCG) (herein termed del18), was identified in our case–control population but was not observed in our original 17 linked families. The del18 deletion occurs at the same locus as the del6 and del12 deletions and results in the loss of an additional serine–proline repeats (for a total of six amino acid residues). This deletion was almost exclusively seen in African American men, with 20/76 (26%) African American men carrying at least one copy of the del18 deletion. In contrast, 1/831 (<0.1%) Caucasian men in the study carried this deletion.
Note that as the 340G and 587C alleles and the 340A and 587T alleles were in complete linkage disequilibrium, we present findings only for the former. Alleles for G340A/587T and C802G were in Hardy–Weinberg equilibrium among controls within ethnic groups (P>0.3).
Associations between the PODXL variants and prostate cancer
No specific model was assumed in our initial regression analysis, and odds ratios (ORs) comparing each observed variant genotype to the wild-type homozygote were calculated. For individual in-frame deletions, positive associations with prostate cancer were observed for almost all of the comparisons (including deletion heterozygotes, homozygotes and compound heterozygotes), with ORs ranging from 2.19 to 3.71 (corresponding P-values ≤0.05). The exceptions were no deletion/del12 heterozygotes where only a weak positive association (OR=1.67) was seen and del18 carriers where no association was found (not shown).
Combining the in-frame deletion genotypes in our family-matched analysis, the presence of any single deletion gave an OR=2.14 [95% confidence interval (CI)=1.09–4.20, P=0.03], whereas the presence of two PODXL deletions further increased this risk, giving an OR=2.58 (95%CI=1.23–5.45, P=0.01) (Table 4, upper panels).
These findings for the in-frame deletion were strengthened among men with more advanced disease at diagnosis. Specifically, for men with more aggressive prostate cancer, positive associations were found for individual in-frame deletions for almost all of the comparisons, with ORs ranging from 4.00 to 7.39 (corresponding to P≤0.03). Once again, the exceptions were a weak positive association for the no deletion/del12 heterozygotes (OR=2.15) and no association for carriers of the del18 (not shown).
Combining the deletion genotypes, the presence of any single deletion gave an OR=3.04 (95%CI=1.01–9.15, P=0.05); carrying two copies of any deletion gave an OR=4.42 (95%CI=1.32–14.85, P=0.02) (Table 4, lower panels). No association was found between any deletion and prostate cancer among men with less aggressive disease (not shown).
Analysis of the G340A variant suggested that carrying one or two of the G340A variants increased prostate cancer risk by approximately 50% (Table 5). Comparing individuals carrying either one or two G340A variants versus none gave an OR=1.48 (95%CI=1.02–2.14, P=0.04). Restricting this analysis to men with more aggressive disease did not strengthen these results. Restricting our analyses to Caucasians only did not materially change any of the results (not shown). No associations were found for other variants.
Two additional variants within the PODXL gene region were genotyped to determine the specificity of the observed associations. The variants are located on either side of the PODXL transcript. One variant (rs15956) lies approximately 8 kb 5′ to the PODXL deletion variant and the other (rs1733877) approximately 54 kb 3′ to the G340A variant. Neither variant exhibited any association with the risk of prostate cancer (OR=1.31, 95%CI=0.85–2.02, P=0.23 for rs15956; OR=1.16, 95%CI=0.74–1.81, P=0.52 for rs1733877) or advanced tumor aggressiveness (OR=0.99, 95%CI=0.54–1.83, P=0.98 for rs15956; OR=1.21, 95%CI=0.67–2.17, P=0.53 for rs1733877) in this population.
DISCUSSION
Both linkage and allelic imbalance studies support the mapping of a gene influencing prostate cancer and tumor aggressiveness on (or near) chromosome 7q32–q33 (5–9). The sialomucin podocalyxin maps within this region, and here we provide evidence that two novel PODXL variants (a variable in-frame deletion and G340A) are associated with the risk of both prostate cancer and tumor aggressiveness in a moderately large-sized family-based prostate cancer case–control population.
We observed a positive association between men carrying one or two deletion copies and prostate cancer, which was strengthened among men with more advanced disease. The corresponding ORs suggest that carrying two copies of the variable in-frame deletion doubles a man's risk of developing more aggressive prostate cancer. The impact of each individual deletion varied slightly, with the strongest association seen with the 6 bp deletion. This may simply reflect the imprecision in the estimates [i.e. because the CIs for the different deletions overlapped] or suggest that different deletion variants have different biological consequences on the podocalyxin protein. In addition, we observed positive associations between men carrying at least one copy of the G340A variant (which was in linkage disequilibrium with the PODXL variant C587T) and prostate cancer risk. Two adjacent variants outside the PODXL transcript region did not show any association with prostate cancer or more aggressive disease, implying that the variable in-frame deletion and G340A variants or other variants on each of their haplotypes (i.e. in linkage disequilibrium) are risk factors for both prostate cancer and aggressive disease.
We also identified an 18 bp deletion (del18) in PODXL, which was found predominantly in African American men and appears to be a founder mutation in this population. African American men are of particular interest because they are at higher risk of developing more aggressive forms of prostate cancer when compared with other ethnic groups (20). Unfortunately, our African American population was too small to determine any associations with this variant, and additional studies focused on clarifying the relationship between the del18 variant and prostate cancer risk in African American men are warranted.
The mechanism by which these variants may affect podocalyxin function and, subsequently, the risk of prostate cancer progression is not known. However, insight into a possible mechanism may come from studies of podocalyxin in the kidney where it plays a critical role in kidney architecture through the maintenance of foot processes and filtration slits through its highly negatively charged extracellular domain. This essential role of podocalyxin in the kidney was demonstrated in Podxl−/− null mice where pups died within 24 h of birth with severe podocyte dysfunction and renal failure (10). This membrane stabilization of podocytes is maintained through a membrane-cytoskeleton-associated complex involving PODXL, Na+/H+-exchanger regulatory factor 1/2 (NHERF1/2), the metastasis-associated protein EZRIN and F-actin (21,22).
Ezrin plays a critical role in a number of signaling pathways implicated in metastasis (23–29). Furthermore, increased ezrin expression correlated with the metastatic potential of both rhabodsarcomas and osteosarcomas (28,30) and with a metastatic phenotype or poor outcome in a number of other cancers (17,26,30–33). Ezrin overexpression has been reported in prostatic intraepithelial neoplasia and prostate tumors (34). However, it should be noted that a role in prostate cancer for podocalyxin and other members of this complex has not yet been reported.
PODXL gene variants may influence prostate cancer and tumor aggressiveness through neutralization of the negative charge of the extracellular domain of podocalyxin. Charge neutralization has previously been shown to disrupt the association between PODXL/NHERF1/2/EZRIN and F-actin (21,22). Alternatively, PODXL variants may perturb downstream targets such as ezrin or NHERF1/2 and contribute to tumor development and metastasis as a result of cell depolarization and increased motility and invasiveness (35). However, it remains to be proven whether the variable in-frame deletion we report in PODXL has any effect on the function of the podocalyxin protein.
There is a precedent for the association of a common in-frame deletion in a protein and the risk of developing cancer: transforming growth factor-beta receptor, Type I (TGFBR1) (36). A common variant (TGFBR1*6A) resulting in the loss of three alanine amino acids within a nine alanine repeat at the end of exon 1 of TGFBR1 results in a reduction in the anti-proliferative effects of TGF beta (37,38). This variant (TGFBR1*6A) has been associated with increased risk of a number of malignancies, including breast, colon and other cancers (37–44).
In summary, we report an association between two PODXL variants and increased risk of developing prostate cancer and increased tumor aggressiveness. Our data suggest that individually or in concert, PODXL variants or variants on their haplotypes may alter the expression or function of the podocalyxin protein possibly through its interactions with the metastasis-related gene EZRIN, leading to an increase in prostate cancer risk. Further studies will be needed to clarify the relationship between podocalyxin and ezrin in the prostate and the role of PODXL variants in prostate cancer aggressiveness.
MATERIALS AND METHODS
Study subjects used in mutation studies
Seventeen probands from brother-pairs that showed strong linkage with D7S1804, the peak linkage marker in our original sibling-pair study (5), were chosen for mutation analysis. Marker D7S1804 was chosen because it also mapped to the distal boundary of the smallest common region of allelic imbalance found in prostate tumors (6). The 17 families showed mean allele sharing between brothers of greater than 70% (with 11 of the 17 families sharing at >95%) at marker D7S1804 with respect to Gleason score. Mutation analyses were performed on DNA isolated from blood leukocytes. Details about recruitment of the subjects that were included in the original linkage analysis are reported in Suarez et al. (45).
Study subjects used in association studies
The family-based case–control study population was recruited through the major medical institutions in the Greater Cleveland, Ohio area and from the Henry Ford Health System in Detroit, Michigan between 1998 and 2001. Men with histologically confirmed prostate cancer were invited to join the study if they fulfilled specific criteria, that is, they were diagnosed at 73 years or younger and had a living, unaffected brother who was either older or at most eight years younger than their age at diagnosis. These criteria were selected in an attempt to increase the likelihood of detecting genetic factors influencing prostate cancer and to minimize potential for misclassification of prostate cancer status among controls. We chose to use this sibling-based study design to ensure that our control group was selected from the same genetic source population as the cases, and thus not to be subject to population stratification from factors such as ethnicity.
The clinical characteristics (e.g. Gleason score, tumor stage) of each case were obtained from medical records. Prostate specific antigen (PSA) levels were measured for all potential controls to screen for possible asymptomatic disease. If a control had PSA levels above 4 ng/ml, they were informed and advised by our collaborating urologists that they should schedule an appointment with their primary care physicians to further evaluate this finding. Such individuals were retained in the study as controls unless a subsequent diagnosis of prostate cancer was made. Their inclusion in the study is important because automatically excluding men with elevated PSA levels regardless of their ultimate prostate cancer status can lead to biased estimates of association (46,47).
A total of 918 men (439 cases and 479 controls) were included in the present study. Each proband had at least one unaffected brother, with some families having multiple affected and unaffected brothers. Ninety-one percent of the subjects were Caucasian, 8% were African American and 1% were Asian or Latino. The mean age at diagnosis of the cases was 61 years and the mean age of controls at enrollment was 62.7 years. This study was approved by the institutional review boards of all participating institutions and all study participants provided written informed consent.
Mutation analysis of PODXL
Genomic DNA samples from the probands of 17 families that showed linkage to chromosome 7q32–q33 were subjected to mutation analyses of the PODXL gene by DNA sequencing. Primers used to amplify all eight exons including the intron splice donor and acceptor regions are shown (Table 1). PCR amplification was performed separately for each exon. Each 15 µl PCR reaction contained 20 ng of DNA, 1.25 mm of each dNTP, 67 mm Tris–HCL (pH 8.8), 6.7 mm magnesium chloride, 16.6 mm ammonium sulfate, 10 mm β-mercaptoethanol, 10% DMSO, 0.5 units Platinum Taq Polymerase (Invitrogen, Carlsbad, CA, USA) and 0.5 µM of forward and reverse primers. The PCR amplification reaction was performed using a 5 min denaturing step at 94°C, 35 cycles of 94°C for 45 s, optimal annealing temperature for 1 min, 72°C for 1 min, followed by a final 5 min extension at 72°C. PCR products were purified using the QIAquick PCR purification kit from QIAGEN (Valencia, CA, USA) according to the manufacturer's instructions. Sequencing was performed using an ABI 377 (Applied Biosystems, Foster City, CA, USA) automated sequencer or a MegaBACE 1000 (Amersham Biosciences, Piscataway, NJ, USA) automated DNA analysis system. Sequences were aligned with and compared with reference PODXL sequence (NM_005397) and analyzed using Sequencher 4.1 (Gene Codes Corp., Ann Arbor, MI, USA).
PODXL genotyping
Several coding sequence variants were identified in PODXL and these are shown in Table 2. Genotyping was performed on those variants that affected amino acid sequence. The PODXL variable in-frame deletion in exon 1 was genotyped following PCR amplification, using the primers 5′-CTCCTCTGCCACTGGCTCTG-3′ and 5′-CGTGGATGGTG CAAGGTCAG. The forward primer included a 5′ fluorophore for analysis using a MegaBACE 1000 DNA analysis system. Allele sizes were assessed using Genetic Profiler v1.5 software (Amersham Biosciences).
The PODXL 340G>A variant was detected following PCR amplification, using the primers 5′-AAGAGCTGGAG AAGCGTGG-3′ and 5′-GGCTGCTTGTGGTGTTAGGT-3′. The 414 bp PCR product was digested with AciI (New England Biolabs, Beverly, MA, USA) for 1 h at 37°C. Digested products were resolved on a 6% polyacrylamide gel. PCR products containing the 340A variant are not cleaved by restriction enzyme AciI, whereas 340G fragments are cleaved into 301 and 113 bp fragments. The PODXL 587C>T variant was detected following PCR amplification, using the primers 5′-AGGCAACCCTACTACCACCA-3′ and 5′-GCTGGGATTACAAGGCATGA-3′. The 459 bp product was digested with Taqα1 (New England Biolabs) at 65°C for 1 h and fragments were separated on a 6% polyacrylamide gel. The 587T PCR product is cleaved into two fragments of 438 and 21 bp, whereas the 587C PCR product is cleaved into three fragments of 221, 217 and 21 bp. The PODXL 802C>G variant was detected by PCR amplification, using the primers 5′-GACCAAAGAACCTGGGAAAA-3′ and 5′-GCGTTGGAGGAAAGAACAGA-3′. The 412 bp product was digested with HaeIII (New England Biolabs) for 1 h at 37°C, resulting in two fragments of 248 and 164 bp for the 802C genotype and three fragments of 221, 164 and 27 bp for the 802G genotype. A reliable assay for the 982G>A variant could not be developed and no further analyses were performed on this polymorphism.
Two additional variants (rs15956 and rs1733877) within the PODXL gene region were genotyped using pre-made TaqMan assays (Applied Biosystems). Reactions were performed according to the manufacturer's instructions, using the RealMasterMix Probe (+ROX) kit (Eppendorf), and genotypes determined using a Prism 7900HT instrument (Applied Biosystems).
Genotyping calls were made independently by two investigators. A randomly selected 10% of samples were re-genotyped and an additional 1% of samples were sequenced using an ABI 377 automated sequencer (Applied Biosystems). No discrepancies in the genotype calls were identified.
Statistical analysis
Genotype frequencies were calculated for each of the alterations being studied, and tested for Hardy–Weinberg equilibrium. Conditional logistic regression (using family as the matching variable, and a robust variance estimator that incorporates familial correlations) was used to estimate ORs and corresponding 95%CIs for the relationship between genotypes and prostate cancer. Sibling sets are treated appropriately by conditional logistic regression in which each family is analyzed as a matched set, so all men are counted only once in the analysis. The regression model was adjusted for potential confounding by age. To investigate the potential effect of genotype on disease aggressiveness, we stratified the analyses by the case's clinical characteristics at diagnosis. Aggressiveness was defined as ‘low’ if a case's Gleason score was <7 and their tumor category was <T2c, and ‘high’ if their Gleason score was ≥7 or their tumor category was ≥T2c. The tumor category reflects the TNM System. All P-values are from two-sided tests, and all analyses were undertaken with SAS software (version 8.2, SAS Institute, Cary, NC, USA).
ACKNOWLEDGEMENTS
This work was supported in part by National Institutes of Health grants RO-1 CA989683 (G.C.), RO-1 CA94211 (J.S.W.) and RO-1 CA88164 (J.S.W.) and by the Urologic Research Foundation (G.C. and J.S.W.).
Conflict of Interest statement. The authors declare no conflicts of interest.
Exon . | Forward primer sequence . | Reverse primer sequence . |
---|---|---|
1 | 5′-CTCCTCTGCCACTGGCTCTG-3′ | 5′-CGTGGATGGTGCAAGGTCAG |
2Aa | 5′-AAGAGCTGGAGAAGCGTGG-3′ | 5′-GGCTGCTTGTGGTGTTAGGT-3′ |
2Ba | 5′-AGGCAACCCTACTACCACCA-3′ | 5′-GCTGGGATTACAAGGCATGA-3′ |
3 | 5′-GACCAAAGAACCTGGGAAAA-3′ | 5′-GCGTTGGAGGAAAGAACAGA-3′ |
4 | 5′-TCTGGGTCACCCCTTAGTTG-3′ | 5′-ATGGGAAAGGACCACTTCCT-3′ |
5 | 5′-ACAGCATGCAACCCTTGATT-3′ | 5′-TCCAATGTGGCTTCTGCAT-3′ |
6 | 5′-ATACGGTGCGTCCCTCTCTT-3′ | 5′-TACATGCAGGCACATGACG-3′ |
7 | 5′-TCATGTGCCTGCATGTAGAA-3′ | 5′-TCTTCTGCAACTCGGGAATC-3′ |
8 | 5′-AGAATGACCACCAGGAGCAG-3′ | 5′-ACTTGGGGTGGTTGGTCTG-3′ |
Exon . | Forward primer sequence . | Reverse primer sequence . |
---|---|---|
1 | 5′-CTCCTCTGCCACTGGCTCTG-3′ | 5′-CGTGGATGGTGCAAGGTCAG |
2Aa | 5′-AAGAGCTGGAGAAGCGTGG-3′ | 5′-GGCTGCTTGTGGTGTTAGGT-3′ |
2Ba | 5′-AGGCAACCCTACTACCACCA-3′ | 5′-GCTGGGATTACAAGGCATGA-3′ |
3 | 5′-GACCAAAGAACCTGGGAAAA-3′ | 5′-GCGTTGGAGGAAAGAACAGA-3′ |
4 | 5′-TCTGGGTCACCCCTTAGTTG-3′ | 5′-ATGGGAAAGGACCACTTCCT-3′ |
5 | 5′-ACAGCATGCAACCCTTGATT-3′ | 5′-TCCAATGTGGCTTCTGCAT-3′ |
6 | 5′-ATACGGTGCGTCCCTCTCTT-3′ | 5′-TACATGCAGGCACATGACG-3′ |
7 | 5′-TCATGTGCCTGCATGTAGAA-3′ | 5′-TCTTCTGCAACTCGGGAATC-3′ |
8 | 5′-AGAATGACCACCAGGAGCAG-3′ | 5′-ACTTGGGGTGGTTGGTCTG-3′ |
aExon 2 was sequenced using two sets of overlapping primers.
Exon . | Forward primer sequence . | Reverse primer sequence . |
---|---|---|
1 | 5′-CTCCTCTGCCACTGGCTCTG-3′ | 5′-CGTGGATGGTGCAAGGTCAG |
2Aa | 5′-AAGAGCTGGAGAAGCGTGG-3′ | 5′-GGCTGCTTGTGGTGTTAGGT-3′ |
2Ba | 5′-AGGCAACCCTACTACCACCA-3′ | 5′-GCTGGGATTACAAGGCATGA-3′ |
3 | 5′-GACCAAAGAACCTGGGAAAA-3′ | 5′-GCGTTGGAGGAAAGAACAGA-3′ |
4 | 5′-TCTGGGTCACCCCTTAGTTG-3′ | 5′-ATGGGAAAGGACCACTTCCT-3′ |
5 | 5′-ACAGCATGCAACCCTTGATT-3′ | 5′-TCCAATGTGGCTTCTGCAT-3′ |
6 | 5′-ATACGGTGCGTCCCTCTCTT-3′ | 5′-TACATGCAGGCACATGACG-3′ |
7 | 5′-TCATGTGCCTGCATGTAGAA-3′ | 5′-TCTTCTGCAACTCGGGAATC-3′ |
8 | 5′-AGAATGACCACCAGGAGCAG-3′ | 5′-ACTTGGGGTGGTTGGTCTG-3′ |
Exon . | Forward primer sequence . | Reverse primer sequence . |
---|---|---|
1 | 5′-CTCCTCTGCCACTGGCTCTG-3′ | 5′-CGTGGATGGTGCAAGGTCAG |
2Aa | 5′-AAGAGCTGGAGAAGCGTGG-3′ | 5′-GGCTGCTTGTGGTGTTAGGT-3′ |
2Ba | 5′-AGGCAACCCTACTACCACCA-3′ | 5′-GCTGGGATTACAAGGCATGA-3′ |
3 | 5′-GACCAAAGAACCTGGGAAAA-3′ | 5′-GCGTTGGAGGAAAGAACAGA-3′ |
4 | 5′-TCTGGGTCACCCCTTAGTTG-3′ | 5′-ATGGGAAAGGACCACTTCCT-3′ |
5 | 5′-ACAGCATGCAACCCTTGATT-3′ | 5′-TCCAATGTGGCTTCTGCAT-3′ |
6 | 5′-ATACGGTGCGTCCCTCTCTT-3′ | 5′-TACATGCAGGCACATGACG-3′ |
7 | 5′-TCATGTGCCTGCATGTAGAA-3′ | 5′-TCTTCTGCAACTCGGGAATC-3′ |
8 | 5′-AGAATGACCACCAGGAGCAG-3′ | 5′-ACTTGGGGTGGTTGGTCTG-3′ |
aExon 2 was sequenced using two sets of overlapping primers.
Location . | Nucleotide changea . | Amino acid changeb . | Allele frequencyc (%) . |
---|---|---|---|
Exon 1 | 67–72 del6 (−6 bp)d | 23–24 del2 | 9/34 (27) |
67–78 del12 (−12 bp)e | 23–26 del4 | 10/34 (29) | |
Exon 2 | 237G>A | None | 5/34 (15) |
Exon 2 | 340G>A | G114S | 13/34 (38) |
Exon 2 | 587C>T | S196L | 13/34 (38) |
Exon 3 | 802C>G | P268A | 2/34 (6) |
Exon 4 | 982G>A | V328I | 2/34 (6) |
Exon 7 | 1320G>A | None | 6/34 (18) |
Location . | Nucleotide changea . | Amino acid changeb . | Allele frequencyc (%) . |
---|---|---|---|
Exon 1 | 67–72 del6 (−6 bp)d | 23–24 del2 | 9/34 (27) |
67–78 del12 (−12 bp)e | 23–26 del4 | 10/34 (29) | |
Exon 2 | 237G>A | None | 5/34 (15) |
Exon 2 | 340G>A | G114S | 13/34 (38) |
Exon 2 | 587C>T | S196L | 13/34 (38) |
Exon 3 | 802C>G | P268A | 2/34 (6) |
Exon 4 | 982G>A | V328I | 2/34 (6) |
Exon 7 | 1320G>A | None | 6/34 (18) |
aNucleotides are numbered from the A(+1) of the ATG-translation initiation codon. Note that the numbering of variants was based on the ‘no deletion’ (full length) variant, whereas the UCSC browser numbers according to the −6 bp deletion variant and differs by two amino acids.
bAmino acids are numbered from the translation initiation methionine (+1). Also see footnote ‘a’.
cAllele frequency reflects the frequency of the variant allele in the total number of chromosomes tested (percentage).
dDeleted amino acids were SerPro.
eDeleted amino acids were SerProSerPro.
Location . | Nucleotide changea . | Amino acid changeb . | Allele frequencyc (%) . |
---|---|---|---|
Exon 1 | 67–72 del6 (−6 bp)d | 23–24 del2 | 9/34 (27) |
67–78 del12 (−12 bp)e | 23–26 del4 | 10/34 (29) | |
Exon 2 | 237G>A | None | 5/34 (15) |
Exon 2 | 340G>A | G114S | 13/34 (38) |
Exon 2 | 587C>T | S196L | 13/34 (38) |
Exon 3 | 802C>G | P268A | 2/34 (6) |
Exon 4 | 982G>A | V328I | 2/34 (6) |
Exon 7 | 1320G>A | None | 6/34 (18) |
Location . | Nucleotide changea . | Amino acid changeb . | Allele frequencyc (%) . |
---|---|---|---|
Exon 1 | 67–72 del6 (−6 bp)d | 23–24 del2 | 9/34 (27) |
67–78 del12 (−12 bp)e | 23–26 del4 | 10/34 (29) | |
Exon 2 | 237G>A | None | 5/34 (15) |
Exon 2 | 340G>A | G114S | 13/34 (38) |
Exon 2 | 587C>T | S196L | 13/34 (38) |
Exon 3 | 802C>G | P268A | 2/34 (6) |
Exon 4 | 982G>A | V328I | 2/34 (6) |
Exon 7 | 1320G>A | None | 6/34 (18) |
aNucleotides are numbered from the A(+1) of the ATG-translation initiation codon. Note that the numbering of variants was based on the ‘no deletion’ (full length) variant, whereas the UCSC browser numbers according to the −6 bp deletion variant and differs by two amino acids.
bAmino acids are numbered from the translation initiation methionine (+1). Also see footnote ‘a’.
cAllele frequency reflects the frequency of the variant allele in the total number of chromosomes tested (percentage).
dDeleted amino acids were SerPro.
eDeleted amino acids were SerProSerPro.
Alleles . | Controls (total=479) . | Cases (total=439) . | |
---|---|---|---|
. | . | Low aggressive . | High aggressive . |
Exon1 | |||
No del | 221 | 82 | 87 |
del6 | 452 | 187 | 227 |
del12 | 272 | 100 | 126 |
del18 | 13 | 5 | 6 |
Exon2 | |||
340G | 551 | 206 | 249 |
340A | 407 | 170 | 201 |
Exon3 | |||
802C | 919 | 344 | 418 |
802G | 69 | 32 | 32 |
Alleles . | Controls (total=479) . | Cases (total=439) . | |
---|---|---|---|
. | . | Low aggressive . | High aggressive . |
Exon1 | |||
No del | 221 | 82 | 87 |
del6 | 452 | 187 | 227 |
del12 | 272 | 100 | 126 |
del18 | 13 | 5 | 6 |
Exon2 | |||
340G | 551 | 206 | 249 |
340A | 407 | 170 | 201 |
Exon3 | |||
802C | 919 | 344 | 418 |
802G | 69 | 32 | 32 |
aSome cases did not have complete clinical data and therefore were not included in the analyses stratified by disease aggressiveness.
Alleles . | Controls (total=479) . | Cases (total=439) . | |
---|---|---|---|
. | . | Low aggressive . | High aggressive . |
Exon1 | |||
No del | 221 | 82 | 87 |
del6 | 452 | 187 | 227 |
del12 | 272 | 100 | 126 |
del18 | 13 | 5 | 6 |
Exon2 | |||
340G | 551 | 206 | 249 |
340A | 407 | 170 | 201 |
Exon3 | |||
802C | 919 | 344 | 418 |
802G | 69 | 32 | 32 |
Alleles . | Controls (total=479) . | Cases (total=439) . | |
---|---|---|---|
. | . | Low aggressive . | High aggressive . |
Exon1 | |||
No del | 221 | 82 | 87 |
del6 | 452 | 187 | 227 |
del12 | 272 | 100 | 126 |
del18 | 13 | 5 | 6 |
Exon2 | |||
340G | 551 | 206 | 249 |
340A | 407 | 170 | 201 |
Exon3 | |||
802C | 919 | 344 | 418 |
802G | 69 | 32 | 32 |
aSome cases did not have complete clinical data and therefore were not included in the analyses stratified by disease aggressiveness.
. | Genotype of affected individual . | Genotype of control . | OR (95%CI) . | ||
---|---|---|---|---|---|
. | . | No deletion (%) . | One deletiona (%) . | Two deletionsa (%) . | . |
All subjects | No deletion | 6 (1.1) | 12 (2.3) | 4 (0.8) | 1.0 (referent) |
One deletiona | 20 (4) | 101 (20) | 49 (10) | 2.14 (1.09–4.20) | |
Two deletionsa | 7 (1.4) | 63 (12) | 251 (49) | 2.58 (1.23–5.45) | |
P-value=0.04 | |||||
Aggressive diseaseb | No deletion | 4 (1.5) | 5 (1.9) | 1 (0.4) | 1.0 (referent) |
One deletiona | 11 (4) | 49 (18) | 24 (9) | 3.04 (1.01–9.15) | |
Two deletionsa | 3 (1.1) | 37 (14) | 134 (50) | 4.42 (1.32–14.85) | |
P-value=0.02 |
. | Genotype of affected individual . | Genotype of control . | OR (95%CI) . | ||
---|---|---|---|---|---|
. | . | No deletion (%) . | One deletiona (%) . | Two deletionsa (%) . | . |
All subjects | No deletion | 6 (1.1) | 12 (2.3) | 4 (0.8) | 1.0 (referent) |
One deletiona | 20 (4) | 101 (20) | 49 (10) | 2.14 (1.09–4.20) | |
Two deletionsa | 7 (1.4) | 63 (12) | 251 (49) | 2.58 (1.23–5.45) | |
P-value=0.04 | |||||
Aggressive diseaseb | No deletion | 4 (1.5) | 5 (1.9) | 1 (0.4) | 1.0 (referent) |
One deletiona | 11 (4) | 49 (18) | 24 (9) | 3.04 (1.01–9.15) | |
Two deletionsa | 3 (1.1) | 37 (14) | 134 (50) | 4.42 (1.32–14.85) | |
P-value=0.02 |
For simplification, the genotype counts shown reflect the number of sibling pairs with specific case–control combinations of genotypes. If there is more than one case and one control in a family, all possible combinations of the sibling-pairs are included. Bold counts indicate pairs where a case carries more copies of the variant than the control and bold and italicized counts the opposite. The difference between these drives the association with prostate cancer. For the ORs shown here, sibling sets are treated appropriately by conditional logistic regression in which each family is analyzed as a matched set, so all men are counted only once in the analysis. This regression provides estimate of ORs and 95% CIs comparing men with one or two PODXL exon 1 in frame deletions to those with none (adjusted for age).
aDeletions of any length.
bIncludes cases with Gleason score ≥7 or tumor stage ≥T2C and their brothers.
. | Genotype of affected individual . | Genotype of control . | OR (95%CI) . | ||
---|---|---|---|---|---|
. | . | No deletion (%) . | One deletiona (%) . | Two deletionsa (%) . | . |
All subjects | No deletion | 6 (1.1) | 12 (2.3) | 4 (0.8) | 1.0 (referent) |
One deletiona | 20 (4) | 101 (20) | 49 (10) | 2.14 (1.09–4.20) | |
Two deletionsa | 7 (1.4) | 63 (12) | 251 (49) | 2.58 (1.23–5.45) | |
P-value=0.04 | |||||
Aggressive diseaseb | No deletion | 4 (1.5) | 5 (1.9) | 1 (0.4) | 1.0 (referent) |
One deletiona | 11 (4) | 49 (18) | 24 (9) | 3.04 (1.01–9.15) | |
Two deletionsa | 3 (1.1) | 37 (14) | 134 (50) | 4.42 (1.32–14.85) | |
P-value=0.02 |
. | Genotype of affected individual . | Genotype of control . | OR (95%CI) . | ||
---|---|---|---|---|---|
. | . | No deletion (%) . | One deletiona (%) . | Two deletionsa (%) . | . |
All subjects | No deletion | 6 (1.1) | 12 (2.3) | 4 (0.8) | 1.0 (referent) |
One deletiona | 20 (4) | 101 (20) | 49 (10) | 2.14 (1.09–4.20) | |
Two deletionsa | 7 (1.4) | 63 (12) | 251 (49) | 2.58 (1.23–5.45) | |
P-value=0.04 | |||||
Aggressive diseaseb | No deletion | 4 (1.5) | 5 (1.9) | 1 (0.4) | 1.0 (referent) |
One deletiona | 11 (4) | 49 (18) | 24 (9) | 3.04 (1.01–9.15) | |
Two deletionsa | 3 (1.1) | 37 (14) | 134 (50) | 4.42 (1.32–14.85) | |
P-value=0.02 |
For simplification, the genotype counts shown reflect the number of sibling pairs with specific case–control combinations of genotypes. If there is more than one case and one control in a family, all possible combinations of the sibling-pairs are included. Bold counts indicate pairs where a case carries more copies of the variant than the control and bold and italicized counts the opposite. The difference between these drives the association with prostate cancer. For the ORs shown here, sibling sets are treated appropriately by conditional logistic regression in which each family is analyzed as a matched set, so all men are counted only once in the analysis. This regression provides estimate of ORs and 95% CIs comparing men with one or two PODXL exon 1 in frame deletions to those with none (adjusted for age).
aDeletions of any length.
bIncludes cases with Gleason score ≥7 or tumor stage ≥T2C and their brothers.
Genotype of affected individual . | Genotype of control . | OR (95%CI)b . | ||
---|---|---|---|---|
. | GG (%) . | GA (%) . | AA (%) . | . |
GG | 104 (20) | 57 (11) | 5 (1) | 1.0 (referent) |
GA | 66 (13) | 129 (25) | 46 (8.9) | 1.48 (1.02–2.16) |
AA | 14 (2.7) | 41 (8.9) | 54 (10) | 1.47 (0.86–2.51) |
AA & GA | 1.48 (1.02–2.14) |
Genotype of affected individual . | Genotype of control . | OR (95%CI)b . | ||
---|---|---|---|---|
. | GG (%) . | GA (%) . | AA (%) . | . |
GG | 104 (20) | 57 (11) | 5 (1) | 1.0 (referent) |
GA | 66 (13) | 129 (25) | 46 (8.9) | 1.48 (1.02–2.16) |
AA | 14 (2.7) | 41 (8.9) | 54 (10) | 1.47 (0.86–2.51) |
AA & GA | 1.48 (1.02–2.14) |
aSee Table 4 for description of table and matched analysis.
bAdjusted for age.
Genotype of affected individual . | Genotype of control . | OR (95%CI)b . | ||
---|---|---|---|---|
. | GG (%) . | GA (%) . | AA (%) . | . |
GG | 104 (20) | 57 (11) | 5 (1) | 1.0 (referent) |
GA | 66 (13) | 129 (25) | 46 (8.9) | 1.48 (1.02–2.16) |
AA | 14 (2.7) | 41 (8.9) | 54 (10) | 1.47 (0.86–2.51) |
AA & GA | 1.48 (1.02–2.14) |
Genotype of affected individual . | Genotype of control . | OR (95%CI)b . | ||
---|---|---|---|---|
. | GG (%) . | GA (%) . | AA (%) . | . |
GG | 104 (20) | 57 (11) | 5 (1) | 1.0 (referent) |
GA | 66 (13) | 129 (25) | 46 (8.9) | 1.48 (1.02–2.16) |
AA | 14 (2.7) | 41 (8.9) | 54 (10) | 1.47 (0.86–2.51) |
AA & GA | 1.48 (1.02–2.14) |
aSee Table 4 for description of table and matched analysis.
bAdjusted for age.
References
Klein, E.A., Kupelian, P.A. and Witte, J.S. (
Kupelian, P.A., Klein, E.A., Witte, J.S., Kupelian, V.A. and Suh, J.H. (
Kupelian, P.A., Kupelian, V.A., Witte, J.S., Macklis, R. and Klein, E.A. (
Rodriguez, C., Calle, E.E., Miracle-McMahill, H.L., Tatham, L.M., Wingo, P.A., Thun, M.J. and Heath, C.W., Jr. (
Witte, J.S., Goddard, K.A., Conti, D.V., Elston, R.C., Lin, J., Suarez, B.K., Broman, K.W., Burmester, J.K., Weber, J.L. and Catalona, W.J. (
Neville, P.J., Conti, D.V., Paris, P.L., Levin, H., Catalona, W.J., Suarez, B.K., Witte, J.S. and Casey, G. (
Paiss, T., Worner, S., Kurtz, F., Haeussler, J., Hautmann, R.E., Gschwend, J.E., Herkommer, K. and Vogel, W. (
Witte, J.S., Suarez, B.K., Theil, B., Lin, J., Yu, A., Banerjee, T.K., Burmester, J.K., Casey, G. and Catalona, W.J. (
Friedrichsen, D.M., Stanford, J.L., Isaacs, S.D., Janer, M., Chang, B.L., Deutsch, K., Gillanders, E., Kolb, S., Wiley, K.E., Badzioch, M.D. et al. (
Doyonnas, R., Kershaw, D.B., Duhme, C., Merkens, H., Chelliah, S., Graf, T. and McNagny, K.M. (
Takeda, T., Go, W.Y., Orlando, R.A. and Farquhar, M.G. (
Sassetti, C., Tangemann, K., Singer, M.S., Kershaw, D.B. and Rosen, S.D. (
Palmer, R.E., Kotsianti, A., Cadman, B., Boyd, T., Gerald, W. and Haber, D.A. (
Guo, J.K., Menke, A.L., Gubler, M.C., Clarke, A.R., Harrison, D., Hammes, A., Hastie, N.D. and Schedl, A. (
Naishiro, Y., Yamada, T., Idogawa, M., Honda, K., Takada, M., Kondo, T., Imai, K. and Hirohashi, S. (
Stanhope-Baker, P., Kessler, P.M., Li, W., Agarwal, M.L. and Williams, B.R. (
Somasiri, A., Nielsen, J.S., Makretsov, N., McCoy, M.L., Prentice, L., Gilks, C.B., Chia, S.K., Gelmon, K.A., Kershaw, D.B., Huntsman, D.G. et al. (
Teruyama, K., Abe, M., Nakano, T., Takahashi, S., Yamada, S. and Sato, Y. (
Kelley, T.W., Huntsman, D., McNagny, K.M., Roskelley, C.D. and Hsi, E.D. (
Takeda, T., McQuistan, T., Orlando, R.A. and Farquhar, M.G. (
Takeda, T. (
Crepaldi, T., Gautreau, A., Comoglio, P.M., Louvard, D. and Arpin, M. (
Gautreau, A., Poullet, P., Louvard, D. and Arpin, M. (
Hiscox, S. and Jiang, W.G. (
Martin, T.A., Harrison, G., Mansel, R.E. and Jiang, W.G. (
Pujuguet, P., Del Maestro, L., Gautreau, A., Louvard, D. and Arpin, M. (
Yu, Y., Khan, J., Khanna, C., Helman, L., Meltzer, P.S. and Merlino, G. (
Khanna, C., Wan, X., Bose, S., Cassaday, R., Olomu, O., Mendoza, A., Yeung, C., Gorlick, R., Hewitt, S.M. and Helman, L.J. (
Makitie, T., Carpen, O., Vaheri, A. and Kivela, T. (
Elliott, B.E., Meens, J.A., Sen Gupta, S.K., Louvard, D. and Arpin, M. (
Ilmonen, S., Vaheri, A., Asko-Seljavaara, S. and Carpen, O. (
Pang, S.T., Fang, X., Valdman, A., Norstedt, G., Pousette, A., Egevad, L. and Ekman, P. (
McClatchey, A.I. (
Pasche, B., Luo, Y., Rao, P.H., Nimer, S.D., Dmitrovsky, E., Caron, P., Luzzatto, L., Offit, K., Cordon-Cardo, C., Renault, B. et al. (
Chen, T., de Vries, E.G., Hollema, H., Yegen, H.A., Vellucci, V.F., Strickler, H.D., Hildesheim, A. and Reiss, M. (
Pasche, B., Kolachana, P., Nafa, K., Satagopan, J., Chen, Y.G., Lo, R.S., Brener, D., Yang, D., Kirstein, L., Oddoux, C. et al. (
Baxter, S.W., Choong, D.Y., Eccles, D.M. and Campbell, I.G. (
Kaklamani, V.G., Baddi, L., Liu, J., Rosman, D., Phukan, S., Bradley, C., Hegarty, C., McDaniel, B., Rademaker, A., Oddoux, C. et al. (
Kaklamani, V.G., Hou, N., Bian, Y., Reich, J., Offit, K., Michel, L.S., Rubinstein, W.S., Rademaker, A. and Pasche, B. (
Pasche, B., Kaklamani, V., Hou, N., Young, T., Rademaker, A., Peterlongo, P., Ellis, N., Offit, K., Caldes, T., Reiss, M. et al. (
Bian, Y., Caldes, T., Wijnen, J., Franken, P., Vasen, H., Kaklamani, V., Nafa, K., Peterlongo, P., Ellis, N., Baron, J.A. et al. (
Pasche, B., Knobloch, T.J., Bian, Y., Liu, J., Phukan, S., Rosman, D., Kaklamani, V., Baddi, L., Siddiqui, F.S., Frankel, W. et al. (
Suarez, B.K., Lin, J., Burmester, J.K., Broman, K.W., Weber, J.L., Banerjee, T.K., Goddard, K.A., Witte, J.S., Elston, R.C. and Catalona, W.J. (
Lubin, J.H. and Hartge, P. (