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 (14).

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 (59). 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 (2329). 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,3033). 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 (3744).

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.

Table 1.

PODXL primers used for PCR amplification and DNA sequencing

ExonForward primer sequenceReverse primer sequence
15′-CTCCTCTGCCACTGGCTCTG-3′5′-CGTGGATGGTGCAAGGTCAG
2Aa5′-AAGAGCTGGAGAAGCGTGG-3′5′-GGCTGCTTGTGGTGTTAGGT-3′
2Ba5′-AGGCAACCCTACTACCACCA-3′5′-GCTGGGATTACAAGGCATGA-3′
35′-GACCAAAGAACCTGGGAAAA-3′5′-GCGTTGGAGGAAAGAACAGA-3′
45′-TCTGGGTCACCCCTTAGTTG-3′5′-ATGGGAAAGGACCACTTCCT-3′
55′-ACAGCATGCAACCCTTGATT-3′5′-TCCAATGTGGCTTCTGCAT-3′
65′-ATACGGTGCGTCCCTCTCTT-3′5′-TACATGCAGGCACATGACG-3′
75′-TCATGTGCCTGCATGTAGAA-3′5′-TCTTCTGCAACTCGGGAATC-3′
85′-AGAATGACCACCAGGAGCAG-3′5′-ACTTGGGGTGGTTGGTCTG-3′
ExonForward primer sequenceReverse primer sequence
15′-CTCCTCTGCCACTGGCTCTG-3′5′-CGTGGATGGTGCAAGGTCAG
2Aa5′-AAGAGCTGGAGAAGCGTGG-3′5′-GGCTGCTTGTGGTGTTAGGT-3′
2Ba5′-AGGCAACCCTACTACCACCA-3′5′-GCTGGGATTACAAGGCATGA-3′
35′-GACCAAAGAACCTGGGAAAA-3′5′-GCGTTGGAGGAAAGAACAGA-3′
45′-TCTGGGTCACCCCTTAGTTG-3′5′-ATGGGAAAGGACCACTTCCT-3′
55′-ACAGCATGCAACCCTTGATT-3′5′-TCCAATGTGGCTTCTGCAT-3′
65′-ATACGGTGCGTCCCTCTCTT-3′5′-TACATGCAGGCACATGACG-3′
75′-TCATGTGCCTGCATGTAGAA-3′5′-TCTTCTGCAACTCGGGAATC-3′
85′-AGAATGACCACCAGGAGCAG-3′5′-ACTTGGGGTGGTTGGTCTG-3′

aExon 2 was sequenced using two sets of overlapping primers.

Table 1.

PODXL primers used for PCR amplification and DNA sequencing

ExonForward primer sequenceReverse primer sequence
15′-CTCCTCTGCCACTGGCTCTG-3′5′-CGTGGATGGTGCAAGGTCAG
2Aa5′-AAGAGCTGGAGAAGCGTGG-3′5′-GGCTGCTTGTGGTGTTAGGT-3′
2Ba5′-AGGCAACCCTACTACCACCA-3′5′-GCTGGGATTACAAGGCATGA-3′
35′-GACCAAAGAACCTGGGAAAA-3′5′-GCGTTGGAGGAAAGAACAGA-3′
45′-TCTGGGTCACCCCTTAGTTG-3′5′-ATGGGAAAGGACCACTTCCT-3′
55′-ACAGCATGCAACCCTTGATT-3′5′-TCCAATGTGGCTTCTGCAT-3′
65′-ATACGGTGCGTCCCTCTCTT-3′5′-TACATGCAGGCACATGACG-3′
75′-TCATGTGCCTGCATGTAGAA-3′5′-TCTTCTGCAACTCGGGAATC-3′
85′-AGAATGACCACCAGGAGCAG-3′5′-ACTTGGGGTGGTTGGTCTG-3′
ExonForward primer sequenceReverse primer sequence
15′-CTCCTCTGCCACTGGCTCTG-3′5′-CGTGGATGGTGCAAGGTCAG
2Aa5′-AAGAGCTGGAGAAGCGTGG-3′5′-GGCTGCTTGTGGTGTTAGGT-3′
2Ba5′-AGGCAACCCTACTACCACCA-3′5′-GCTGGGATTACAAGGCATGA-3′
35′-GACCAAAGAACCTGGGAAAA-3′5′-GCGTTGGAGGAAAGAACAGA-3′
45′-TCTGGGTCACCCCTTAGTTG-3′5′-ATGGGAAAGGACCACTTCCT-3′
55′-ACAGCATGCAACCCTTGATT-3′5′-TCCAATGTGGCTTCTGCAT-3′
65′-ATACGGTGCGTCCCTCTCTT-3′5′-TACATGCAGGCACATGACG-3′
75′-TCATGTGCCTGCATGTAGAA-3′5′-TCTTCTGCAACTCGGGAATC-3′
85′-AGAATGACCACCAGGAGCAG-3′5′-ACTTGGGGTGGTTGGTCTG-3′

aExon 2 was sequenced using two sets of overlapping primers.

Table 2.

Sequence alterations identified in PODXL

LocationNucleotide changeaAmino acid changebAllele frequencyc (%)
Exon 167–72 del6 (−6 bp)d23–24 del29/34 (27)
67–78 del12 (−12 bp)e23–26 del410/34 (29)
Exon 2237G>ANone5/34 (15)
Exon 2340G>AG114S13/34 (38)
Exon 2587C>TS196L13/34 (38)
Exon 3802C>GP268A2/34 (6)
Exon 4982G>AV328I2/34 (6)
Exon 71320G>ANone6/34 (18)
LocationNucleotide changeaAmino acid changebAllele frequencyc (%)
Exon 167–72 del6 (−6 bp)d23–24 del29/34 (27)
67–78 del12 (−12 bp)e23–26 del410/34 (29)
Exon 2237G>ANone5/34 (15)
Exon 2340G>AG114S13/34 (38)
Exon 2587C>TS196L13/34 (38)
Exon 3802C>GP268A2/34 (6)
Exon 4982G>AV328I2/34 (6)
Exon 71320G>ANone6/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.

Table 2.

Sequence alterations identified in PODXL

LocationNucleotide changeaAmino acid changebAllele frequencyc (%)
Exon 167–72 del6 (−6 bp)d23–24 del29/34 (27)
67–78 del12 (−12 bp)e23–26 del410/34 (29)
Exon 2237G>ANone5/34 (15)
Exon 2340G>AG114S13/34 (38)
Exon 2587C>TS196L13/34 (38)
Exon 3802C>GP268A2/34 (6)
Exon 4982G>AV328I2/34 (6)
Exon 71320G>ANone6/34 (18)
LocationNucleotide changeaAmino acid changebAllele frequencyc (%)
Exon 167–72 del6 (−6 bp)d23–24 del29/34 (27)
67–78 del12 (−12 bp)e23–26 del410/34 (29)
Exon 2237G>ANone5/34 (15)
Exon 2340G>AG114S13/34 (38)
Exon 2587C>TS196L13/34 (38)
Exon 3802C>GP268A2/34 (6)
Exon 4982G>AV328I2/34 (6)
Exon 71320G>ANone6/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.

Table 3.

Allele frequencies for PODXL variants in the family-based case–control studya

AllelesControls (total=479)Cases (total=439)
Low aggressiveHigh aggressive
Exon1
 No del2218287
 del6452187227
 del12272100126
 del181356
Exon2
340G551206249
340A407170201
Exon3
802C919344418
802G693232
AllelesControls (total=479)Cases (total=439)
Low aggressiveHigh aggressive
Exon1
 No del2218287
 del6452187227
 del12272100126
 del181356
Exon2
340G551206249
340A407170201
Exon3
802C919344418
802G693232

aSome cases did not have complete clinical data and therefore were not included in the analyses stratified by disease aggressiveness.

Table 3.

Allele frequencies for PODXL variants in the family-based case–control studya

AllelesControls (total=479)Cases (total=439)
Low aggressiveHigh aggressive
Exon1
 No del2218287
 del6452187227
 del12272100126
 del181356
Exon2
340G551206249
340A407170201
Exon3
802C919344418
802G693232
AllelesControls (total=479)Cases (total=439)
Low aggressiveHigh aggressive
Exon1
 No del2218287
 del6452187227
 del12272100126
 del181356
Exon2
340G551206249
340A407170201
Exon3
802C919344418
802G693232

aSome cases did not have complete clinical data and therefore were not included in the analyses stratified by disease aggressiveness.

Table 4.

Association between the PODXL exon 1 variable in-frame deletion genotype and prostate cancer

Genotype of affected individualGenotype of controlOR (95%CI)
No deletion (%)One deletiona (%)Two deletionsa (%)
All subjectsNo deletion6 (1.1)12 (2.3)4 (0.8)1.0 (referent)
One deletiona20 (4)101 (20)49 (10)2.14 (1.09–4.20)
Two deletionsa7 (1.4)63 (12)251 (49)2.58 (1.23–5.45)
P-value=0.04
Aggressive diseasebNo deletion4 (1.5)5 (1.9)1 (0.4)1.0 (referent)
One deletiona11 (4)49 (18)24 (9)3.04 (1.01–9.15)
Two deletionsa3 (1.1)37 (14)134 (50)4.42 (1.32–14.85)
P-value=0.02
Genotype of affected individualGenotype of controlOR (95%CI)
No deletion (%)One deletiona (%)Two deletionsa (%)
All subjectsNo deletion6 (1.1)12 (2.3)4 (0.8)1.0 (referent)
One deletiona20 (4)101 (20)49 (10)2.14 (1.09–4.20)
Two deletionsa7 (1.4)63 (12)251 (49)2.58 (1.23–5.45)
P-value=0.04
Aggressive diseasebNo deletion4 (1.5)5 (1.9)1 (0.4)1.0 (referent)
One deletiona11 (4)49 (18)24 (9)3.04 (1.01–9.15)
Two deletionsa3 (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.

Table 4.

Association between the PODXL exon 1 variable in-frame deletion genotype and prostate cancer

Genotype of affected individualGenotype of controlOR (95%CI)
No deletion (%)One deletiona (%)Two deletionsa (%)
All subjectsNo deletion6 (1.1)12 (2.3)4 (0.8)1.0 (referent)
One deletiona20 (4)101 (20)49 (10)2.14 (1.09–4.20)
Two deletionsa7 (1.4)63 (12)251 (49)2.58 (1.23–5.45)
P-value=0.04
Aggressive diseasebNo deletion4 (1.5)5 (1.9)1 (0.4)1.0 (referent)
One deletiona11 (4)49 (18)24 (9)3.04 (1.01–9.15)
Two deletionsa3 (1.1)37 (14)134 (50)4.42 (1.32–14.85)
P-value=0.02
Genotype of affected individualGenotype of controlOR (95%CI)
No deletion (%)One deletiona (%)Two deletionsa (%)
All subjectsNo deletion6 (1.1)12 (2.3)4 (0.8)1.0 (referent)
One deletiona20 (4)101 (20)49 (10)2.14 (1.09–4.20)
Two deletionsa7 (1.4)63 (12)251 (49)2.58 (1.23–5.45)
P-value=0.04
Aggressive diseasebNo deletion4 (1.5)5 (1.9)1 (0.4)1.0 (referent)
One deletiona11 (4)49 (18)24 (9)3.04 (1.01–9.15)
Two deletionsa3 (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.

Table 5.

Association between the PODXL G340A variants and prostate cancera

Genotype of affected individualGenotype of controlOR (95%CI)b
GG (%)GA (%)AA (%)
GG104 (20)57 (11)5 (1)1.0 (referent)
GA66 (13)129 (25)46 (8.9)1.48 (1.02–2.16)
AA14 (2.7)41 (8.9)54 (10)1.47 (0.86–2.51)
AA & GA1.48 (1.02–2.14)
Genotype of affected individualGenotype of controlOR (95%CI)b
GG (%)GA (%)AA (%)
GG104 (20)57 (11)5 (1)1.0 (referent)
GA66 (13)129 (25)46 (8.9)1.48 (1.02–2.16)
AA14 (2.7)41 (8.9)54 (10)1.47 (0.86–2.51)
AA & GA1.48 (1.02–2.14)

aSee Table 4 for description of table and matched analysis.

bAdjusted for age.

Table 5.

Association between the PODXL G340A variants and prostate cancera

Genotype of affected individualGenotype of controlOR (95%CI)b
GG (%)GA (%)AA (%)
GG104 (20)57 (11)5 (1)1.0 (referent)
GA66 (13)129 (25)46 (8.9)1.48 (1.02–2.16)
AA14 (2.7)41 (8.9)54 (10)1.47 (0.86–2.51)
AA & GA1.48 (1.02–2.14)
Genotype of affected individualGenotype of controlOR (95%CI)b
GG (%)GA (%)AA (%)
GG104 (20)57 (11)5 (1)1.0 (referent)
GA66 (13)129 (25)46 (8.9)1.48 (1.02–2.16)
AA14 (2.7)41 (8.9)54 (10)1.47 (0.86–2.51)
AA & GA1.48 (1.02–2.14)

aSee Table 4 for description of table and matched analysis.

bAdjusted for age.

References

1

Klein, E.A., Kupelian, P.A. and Witte, J.S. (

1998
) Does a family history of prostate cancer result in more aggressive disease?
Prostate Cancer Prostatic Dis.
,
1
,
297
–300.

2

Kupelian, P.A., Klein, E.A., Witte, J.S., Kupelian, V.A. and Suh, J.H. (

1997
) Familial prostate cancer: a different disease?
J. Urol.
,
158
,
2197
–2201.

3

Kupelian, P.A., Kupelian, V.A., Witte, J.S., Macklis, R. and Klein, E.A. (

1997
) Family history of prostate cancer in patients with localized prostate cancer: an independent predictor of treatment outcome.
J. Clin. Oncol.
,
15
,
1478
–1480.

4

Rodriguez, C., Calle, E.E., Miracle-McMahill, H.L., Tatham, L.M., Wingo, P.A., Thun, M.J. and Heath, C.W., Jr. (

1997
) Family history and risk of fatal prostate cancer.
Epidemiology
,
8
,
653
–657.

5

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

2000
) Genomewide scan for prostate cancer-aggressiveness loci.
Am. J. Hum. Genet.
,
67
,
92
–99.

6

Neville, P.J., Conti, D.V., Paris, P.L., Levin, H., Catalona, W.J., Suarez, B.K., Witte, J.S. and Casey, G. (

2002
) Prostate cancer aggressiveness locus on chromosome 7q32–q33 identified by linkage and allelic imbalance studies.
Neoplasia
,
4
,
424
–431.

7

Paiss, T., Worner, S., Kurtz, F., Haeussler, J., Hautmann, R.E., Gschwend, J.E., Herkommer, K. and Vogel, W. (

2003
) Linkage of aggressive prostate cancer to chromosome 7q31–33 in German prostate cancer families.
Eur. J. Hum. Genet.
,
11
,
17
–22.

8

Witte, J.S., Suarez, B.K., Theil, B., Lin, J., Yu, A., Banerjee, T.K., Burmester, J.K., Casey, G. and Catalona, W.J. (

2003
) Genome-wide scan of brothers: replication and fine mapping of prostate cancer susceptibility and aggressiveness loci.
Prostate
,
57
,
298
–308.

9

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

2004
) Identification of a prostate cancer susceptibility locus on chromosome 7q11–21 in Jewish families.
Proc. Natl Acad. Sci. USA
,
101
,
1939
–1944.

10

Doyonnas, R., Kershaw, D.B., Duhme, C., Merkens, H., Chelliah, S., Graf, T. and McNagny, K.M. (

2001
) Anuria, omphalocele, and perinatal lethality in mice lacking the CD34-related protein podocalyxin.
J. Exp. Med.
,
194
,
13
–27.

11

Takeda, T., Go, W.Y., Orlando, R.A. and Farquhar, M.G. (

2000
) Expression of podocalyxin inhibits cell–cell adhesion and modifies junctional properties in Madin–Darby canine kidney cells.
Mol. Biol. Cell
,
11
,
3219
–3232.

12

Sassetti, C., Tangemann, K., Singer, M.S., Kershaw, D.B. and Rosen, S.D. (

1998
) Identification of podocalyxin-like protein as a high endothelial venule ligand for L-selectin: parallels to CD34.
J. Exp. Med.
,
187
,
1965
–1975.

13

Palmer, R.E., Kotsianti, A., Cadman, B., Boyd, T., Gerald, W. and Haber, D.A. (

2001
) WT1 regulates the expression of the major glomerular podocyte membrane protein podocalyxin.
Curr. Biol.
,
11
,
1805
–1809.

14

Guo, J.K., Menke, A.L., Gubler, M.C., Clarke, A.R., Harrison, D., Hammes, A., Hastie, N.D. and Schedl, A. (

2002
) WT1 is a key regulator of podocyte function: reduced expression levels cause crescentic glomerulonephritis and mesangial sclerosis.
Hum. Mol. Genet.
,
11
,
651
–659.

15

Naishiro, Y., Yamada, T., Idogawa, M., Honda, K., Takada, M., Kondo, T., Imai, K. and Hirohashi, S. (

2005
) Morphological and transcriptional responses of untransformed intestinal epithelial cells to an oncogenic beta-catenin protein.
Oncogene
,
24
,
3141
–3153.

16

Stanhope-Baker, P., Kessler, P.M., Li, W., Agarwal, M.L. and Williams, B.R. (

2004
) The Wilms tumor suppressor-1 target gene podocalyxin is transcriptionally repressed by p53.
J. Biol. Chem.
,
279
,
33575
–33585.

17

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

2004
) Overexpression of the anti-adhesin podocalyxin is an independent predictor of breast cancer progression.
Cancer Res.
,
64
,
5068
–5073.

18

Teruyama, K., Abe, M., Nakano, T., Takahashi, S., Yamada, S. and Sato, Y. (

2001
) Neurophilin-1 is a downstream target of transcription factor Ets-1 in human umbilical vein endothelial cells.
FEBS Lett.
,
504
,
1
–4.

19

Kelley, T.W., Huntsman, D., McNagny, K.M., Roskelley, C.D. and Hsi, E.D. (

2005
) Podocalyxin: a marker of blasts in acute leukemia.
Am. J. Clin. Pathol.
,
124
,
134
–142.

20

American Cancer Society. (

2005
)
Cancer Facts and Figures 2005
. Atlanta: American Cancer Society.

21

Takeda, T., McQuistan, T., Orlando, R.A. and Farquhar, M.G. (

2001
) Loss of glomerular foot processes is associated with uncoupling of podocalyxin from the actin cytoskeleton.
J. Clin. Invest.
,
108
,
289
–301.

22

Takeda, T. (

2003
) Podocyte cytoskeleton is connected to the integral membrane protein podocalyxin through Na+/H+-exchanger regulatory factor 2 and ezrin.
Clin. Exp. Nephrol.
,
7
,
260
–269.

23

Crepaldi, T., Gautreau, A., Comoglio, P.M., Louvard, D. and Arpin, M. (

1997
) Ezrin is an effector of hepatocyte growth factor-mediated migration and morphogenesis in epithelial cells.
J. Cell Biol.
,
138
,
423
–434.

24

Gautreau, A., Poullet, P., Louvard, D. and Arpin, M. (

1999
) Ezrin, a plasma membrane-microfilament linker, signals cell survival through the phosphatidylinositol 3-kinase/Akt pathway.
Proc. Natl Acad. Sci. USA
,
96
,
7300
–7305.

25

Hiscox, S. and Jiang, W.G. (

1999
) Ezrin regulates cell–cell and cell-matrix adhesion, a possible role with E-cadherin/beta-catenin.
J. Cell Sci.
,
112
,
3081
–3090.

26

Martin, T.A., Harrison, G., Mansel, R.E. and Jiang, W.G. (

2003
) The role of the CD44/ezrin complex in cancer metastasis.
Crit. Rev. Oncol. Hematol.
,
46
,
165
–186.

27

Pujuguet, P., Del Maestro, L., Gautreau, A., Louvard, D. and Arpin, M. (

2003
) Ezrin regulates E-cadherin-dependent adherens junction assembly through Rac1 activation.
Mol. Biol. Cell
,
14
,
2181
–2191.

28

Yu, Y., Khan, J., Khanna, C., Helman, L., Meltzer, P.S. and Merlino, G. (

2004
) Expression profiling identifies the cytoskeletal organizer ezrin and the developmental homeoprotein Six-1 as key metastatic regulators.
Nat. Med.
,
10
,
175
–181.

29

Hunter, K.W. (

2004
) Ezrin, a key component in tumor metastasis.
Trends Mol. Med.
,
10
,
201
–204.

30

Khanna, C., Wan, X., Bose, S., Cassaday, R., Olomu, O., Mendoza, A., Yeung, C., Gorlick, R., Hewitt, S.M. and Helman, L.J. (

2004
) The membrane-cytoskeleton linker ezrin is necessary for osteosarcoma metastasis.
Nat. Med.
,
10
,
182
–186.

31

Makitie, T., Carpen, O., Vaheri, A. and Kivela, T. (

2001
) Ezrin as a prognostic indicator and its relationship to tumor characteristics in uveal malignant melanoma.
Invest. Ophthalmol. Vis. Sci.
,
42
,
2442
–2449.

32

Elliott, B.E., Meens, J.A., Sen Gupta, S.K., Louvard, D. and Arpin, M. (

2005
) The membrane cytoskeletal crosslinker ezrin is required for metastasis of breast carcinoma cells.
Breast Cancer Res.
,
7
,
365
–373.

33

Ilmonen, S., Vaheri, A., Asko-Seljavaara, S. and Carpen, O. (

2005
) Ezrin in primary cutaneous melanoma.
Mod. Pathol.
,
18
,
503
–510.

34

Pang, S.T., Fang, X., Valdman, A., Norstedt, G., Pousette, A., Egevad, L. and Ekman, P. (

2004
) Expression of ezrin in prostatic intraepithelial neoplasia.
Urology
,
63
,
609
–612.

35

McClatchey, A.I. (

2003
) Merlin and ERM proteins: unappreciated roles in cancer development?
Nat. Rev. Cancer
,
3
,
877
–883.

36

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

1998
) Type I transforming growth factor beta receptor maps to 9q22 and exhibits a polymorphism and a rare variant within a polyalanine tract.
Cancer Res.
,
58
,
2727
–2732.

37

Chen, T., de Vries, E.G., Hollema, H., Yegen, H.A., Vellucci, V.F., Strickler, H.D., Hildesheim, A. and Reiss, M. (

1999
) Structural alterations of transforming growth factor-beta receptor genes in human cervical carcinoma.
Int. J. Cancer
,
82
,
43
–51.

38

Pasche, B., Kolachana, P., Nafa, K., Satagopan, J., Chen, Y.G., Lo, R.S., Brener, D., Yang, D., Kirstein, L., Oddoux, C. et al. (

1999
) TbetaR-I(6A) is a candidate tumor susceptibility allele.
Cancer Res.
,
59
,
5678
–5682.

39

Baxter, S.W., Choong, D.Y., Eccles, D.M. and Campbell, I.G. (

2002
) Transforming growth factor beta receptor 1 polyalanine polymorphism and exon 5 mutation analysis in breast and ovarian cancer.
Cancer Epidemiol. Biomarkers Prev.
,
11
,
211
–214.

40

Kaklamani, V.G., Baddi, L., Liu, J., Rosman, D., Phukan, S., Bradley, C., Hegarty, C., McDaniel, B., Rademaker, A., Oddoux, C. et al. (

2005
) Combined genetic assessment of transforming growth factor-beta signaling pathway variants may predict breast cancer risk.
Cancer Res.
,
65
,
3454
–3461.

41

Kaklamani, V.G., Hou, N., Bian, Y., Reich, J., Offit, K., Michel, L.S., Rubinstein, W.S., Rademaker, A. and Pasche, B. (

2003
) TGFBR1*6A and cancer risk: a meta-analysis of seven case–control studies.
J. Clin. Oncol.
,
21
,
3236
–3243.

42

Pasche, B., Kaklamani, V., Hou, N., Young, T., Rademaker, A., Peterlongo, P., Ellis, N., Offit, K., Caldes, T., Reiss, M. et al. (

2004
) TGFBR1*6A and cancer: a meta-analysis of 12 case–control studies.
J. Clin. Oncol.
,
22
,
756
–758.

43

Bian, Y., Caldes, T., Wijnen, J., Franken, P., Vasen, H., Kaklamani, V., Nafa, K., Peterlongo, P., Ellis, N., Baron, J.A. et al. (

2005
) TGFBR1*6A may contribute to hereditary colorectal cancer.
J. Clin. Oncol.
,
23
,
3074
–3078.

44

Pasche, B., Knobloch, T.J., Bian, Y., Liu, J., Phukan, S., Rosman, D., Kaklamani, V., Baddi, L., Siddiqui, F.S., Frankel, W. et al. (

2005
) Somatic acquisition and signaling of TGFBR1*6A in cancer.
JAMA
,
294
,
1634
–1646.

45

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

2000
) A genome screen of multiplex sibships with prostate cancer.
Am. J. Hum. Genet.
,
66
,
933
–944.

46

Lubin, J.H. and Hartge, P. (

1984
) Excluding controls: misapplications in case–control studies.
Am. J. Epidemiol.
,
120
,
791
–793.

47

Poole, C. (

1999
) Controls who experienced hypothetical causal intermediates should not be excluded from case–control studies.
Am. J. Epidemiol.
,
150
,
547
–551.