Search

Supplementary Methods, Tables 1-3, Figure 1 from A Risk Variant in an miR-125b Binding Site in BMPR1B Is Associated with Breast Cancer Pathogenesis

Supplementary Figure and Table Legends from A Risk Variant in an miR-125b Binding Site in BMPR1B Is Associated with Breast Cancer Pathogenesis

MicroRNAs regulate diverse cellular processes and play an integral role in cancer pathogenesis. Genomic variation within miRNA target sites may therefore be important sources for genetic differences in cancer risk. To investigate this possibility, we mapped HapMap single nucleotide polymorphisms (SNP) to putative miRNA recognition sites within genes dysregulated in estrogen receptor–stratified breast tumors and used local linkage disequilibirum patterns to identify high-ranking SNPs in the Cancer Genetic Markers of Susceptibility (CGEMS) breast cancer genome-wide association study for further testing. Two SNPs, rs1970801 and rs11097457, scoring in the top 100 from the CGEMS study, were in strong linkage disequilibrium with rs1434536, an SNP that resides within a miR-125b target site in the 3′ untranslated region of the bone morphogenic receptor type 1B (BMPR1B) gene encoding a transmembrane serine/threonine kinase. We validated the CGEMS association findings for rs1970801 in an independent cohort of admixture-corrected cases identified from families with multiple case histories. Subsequent association testing of rs1434536 for these cases and CGEMS controls with imputed genotypes supported the association. Furthermore, luciferase reporter assays and overexpression of miR-125b–mimics combined with quantitative reverse transcription-PCR showed that BMPR1B transcript is a direct target of miR-125b and that miR-125b differentially regulates the C and T alleles of rs1434536. These results suggest that allele-specific regulation of BMPR1B by miR-125b explains the observed disease risk. Our approach is general and can help identify and explain the mechanisms behind disease association for alleles that affect miRNA regulation. [Cancer Res 2009;69(18):7459–65]

Supplementary Figure from Genetic Linkage of Prostate Cancer Risk to the Chromosome 3 Region Bearing FHIT

We conducted linkage analysis of 80 candidate genes in 201 brother pairs affected with prostatic adenocarcinoma. Markers representing two adjacent candidate genes on chromosome 3p, CDC25A and FHIT, showed suggestive evidence for linkage with single-point identity-by-descent allele-sharing statistics. Fine-structure multipoint linkage analysis yielded a maximum LOD score of 3.17 (P = 0.00007) at D3S1234 within FHIT intron 5. For a subgroup of 38 families in which three or more affected brothers were reported, the LOD score was 3.83 (P = 0.00001). Further analysis reported herein suggested a recessive mode of inheritance. Association testing of 16 single nucleotide polymorphisms (SNP) spanning a 381-kb interval surrounding D3S1234 in 202 cases of European descent with 143 matched, unrelated controls revealed significant evidence for association between case status and the A allele of single nucleotide polymorphism rs760317, located within intron 5 of FHIT (Pearson's χ2 = 8.54, df = 1, P = 0.0035). Our results strongly suggest involvement of germline variations of FHIT in prostate cancer risk.

Posttranslational modification by the ubiquitin homologue, small ubiquitin-like modifier 1 (SUMO-1), has been established as an important regulatory mechanism. However, in most cases it is not clear how sumoylation regulates various cellular functions. Emerging evidence suggests that sumoylation may play a general role in regulating protein-protein interactions, as shown in RanBP2/Nup358 and RanGAP1 interaction. In this study, we have defined an amino acid sequence motif that binds SUMO. This motif, V/I-X-V/I-V/I, was identified by NMR spectroscopic characterization of interactions among SUMO-1 and peptides derived from proteins that are known to bind SUMO or sumoylated proteins. This motif binds all SUMO paralogues (SUMO-1-3). Using site-directed mutagenesis, we also show that this SUMO-binding motif in RanBP2/Nup358 is responsible for the interaction between RanBP2/Nup358 and sumoylated RanGAP1. The SUMO-binding motif exists in nearly all proteins known to be involved in SUMO-dependent processes, suggesting its general role in sumoylation-dependent cellular functions.

Linkage and association analyses were performed to identify loci affecting disease susceptibility by scoring previously characterized sequence variations such as microsatellites and single nucleotide polymorphisms. Lack of markers in regions of interest, as well as difficulty in adapting various methods to high-throughput settings, often limits the effectiveness of the analyses. We have adapted the Escherichia coli mismatch detection system, employing the factors MutS, MutL and MutH, for use in PCR-based, automated, high-throughput genotyping and mutation detection of genomic DNA. Optimal sensitivity and signal-to-noise ratios were obtained in a straightforward fashion because the detection reaction proved to be principally dependent upon monovalent cation concentration and MutL concentration. Quantitative relationships of the optimal values of these parameters with length of the DNA test fragment were demonstrated, in support of the translocation model for the mechanism of action of these enzymes, rather than the molecular switch model. Thus, rapid, sequence-independent optimization was possible for each new genomic target region. Other factors potentially limiting the flexibility of mismatch scanning, such as positioning of dam recognition sites within the target fragment, have also been investigated. We developed several strategies, which can be easily adapted to automation, for limiting the analysis to intersample heteroduplexes. Thus, the principal barriers to the use of this methodology, which we have designated PCR candidate region mismatch scanning, in cost-effective, high-throughput settings have been removed.

The t(14,18) chromosomal translocation that occurs in human follicular lymphoma constitutively activates the BCL2 gene and disrupts control of apoptosis. Interestingly, 70% of the t(14,18) translocations are confined to three 15-bp clusters positioned within a 150-bp region (major breakpoint region or [MBR]) in the untranslated portion of terminal exon 3. We analyzed DNA-protein interactions in the MBR, as these may play some role in targeting the translocation to this region. An 87-bp segment (87MBR) immediately 3′ to breakpoint cluster 3 was essential for DNA-protein interaction monitored with mobility shift assays. We further delineated a core binding region within 87MBR: a 33-bp, very AT-rich sequence highly conserved between the human and mouse BCL2 gene (37MBR). We have purified and identified one of the core factors as the matrix attachment region (MAR) binding protein, SATB1, which is known to bind to AT-rich sequences with a high propensity to unwind. Additional factors in nuclear extracts, which we have not yet characterized further, increased SATB1 affinity for the 37MBR target four- to fivefold. Specific binding activity within 37MBR displayed cell cycle regulation in Jurkat T cells, while levels of SATB1 remained constant throughout the cell cycle. Finally, we demonstrated in vivo binding of SATB1 to the MBR, strongly suggesting the BCL2 major breakpoint region is a MAR. We discuss the potential consequences of our observations for both MBR fragility and regulatory function.

DNA sequence analysis of 130 alleles of the HRAS1 minisatellite has demonstrated that breast cancer-associated variants arise as a consequence of both replication errors and gene conversions. Unlike mutations at other variable number of tandem repeats (VNTRs), high-risk variants of the HRAS1 minisatellite do not demonstrate positional polarity. Instead, most mutations occur at three hotspots, with replication errors confined to one hotspot, gene conversions to a second and a mixed pattern of mutation at the third. DNA sequence analysis of 66 low-risk a1 alleles revealed no evidence for hypermutation. Therefore, while the HRAS1 minisatellite may serve as a reporter for a broad-based group of mutational mechanisms, these results are consistent with a direct pathogenetic contribution by high-risk alleles as the biological basis underlying cancer association of this VNTR.

We have tested a simple procedure, disease association by locus stratification, for identifying breast cancer patients with pathogenetic allelic variants at several candidate loci. The strategy was based on the assumption of epistatic interactions of the candidates. We analyzed 66 independent cases from sib pairs affected with breast cancer that had previously been collected during an investigation of pathogenetic-allele-sharing at the HRAS1 mini-satellite locus. An exon 24 polymorphism of ATM, substituting arginine for proline was associated with breast cancer in these cases with an overall odds ratio of 4.5 (95% confidence interval, 1.2–20.5, nominal p = 0.02, 2–tail Fisher exact test). In the presence of a rare HRAS1 allele, the odds ratio increased to 6.9 (95% CI, 1.2–38.3, p = 0.03). Thus, our procedure identified at least one allelic variant of ATM associated with breast cancer, and indicated that the ATM locus may interact with HRAS1.

The majority of translocations involving BCL2 are very narrowly targeted to three breakpoint clusters evenly spaced over a 100-bp region of the gene's terminal exon. We have recently shown that the immediate upstream boundary of this major breakpoint region (mbr) is a specific recognition site for single-strand DNA (ssDNA) binding proteins on the sense and antisense strands. The downstream flank of the mbr is a helicase binding site. In this report we demonstrate that the helicase and ssDNA binding proteins show reciprocal changes in binding activity over the cell cycle. The helicase is maximally active in G1 and early S phases; the ssDNA binding proteins are maximally active in late S and G2/M phases. An inhibitor of helicase binding appears in late S and G2/M. Finally, at least one component of the helicase binding complex is the Ku antigen. Thus, a protein with helicase activity implicated in repair of double-strand breaks, variable (diversity) joining recombination, and, potentially, cell-cycle regulation is targeted to the BCL2 mbr.

We describe three statistical results that we have found to be useful in case-control genetic association testing. All three involve combining the discovery of novel genetic variants, usually by sequencing, with genotyping methods that recognize previously discovered variants. We first consider expanding the list of known variants by concentrating variant-discovery in cases. Although the naive inclusion of cases-only sequencing data would create a bias, we show that some sequencing data may be retained, even if controls are not sequenced. Furthermore, for alleles of intermediate frequency, cases-only sequencing with bias-correction entails little if any loss of power, compared to dividing the same sequencing effort among cases and controls. Secondly, we investigate more strongly focused variant discovery to obtain a greater enrichment for disease-related variants. We show how case status, family history, and marker sharing enrich the discovery set by increments that are multiplicative with penetrance, enabling the preferential discovery of high-penetrance variants. A third result applies when sequencing is the primary means of counting alleles in both cases and controls, but a supplementary pooled genotyping sample is used to identify the variants that are very rare. We show that this raises no validity issues, and we evaluate a less expensive and more adaptive approach to judging rarity, based on group-specific variants. We demonstrate the important and unusual caveat that this method requires equal sample sizes for validity. These three results can be used to more efficiently detect the association of rare genetic variants with disease.

We have examined an unselected series of 72 lymphomas of diverse histologies with a panel of mouse monoclonal antibodies specific for human B-lymphoma-derived Ig idiotypes (anti-ids) to determine the nature and extent of id/anti-id crossreactivity. The anti-id antibodies were prepared from Ig isolated from seven follicular center cell lymphomas by heterohybridoma technique. Thirty-six of 75 individual anti-ids obtained in this manner were further selected based on their reactivity with highly restricted or private idiotypic determinants. Twelve of the 72 (17%) lymphoma biopsies reacted with one or more of the 36 anti-ids that detect private determinants. A pool consisting of five individual antibodies would have detected 11 of the 12 crossreacting tumors. Staining of tumor cell populations was homogeneous in the positive cases, suggesting uniform idiotype expression. If there was a segregated staining pattern, it was generally related to the presence of CD3+ T cells in the section. These follicular center cell-derived anti-ids cross-reacted with follicular center cell tumors of all histologic grades with frequencies ranging from 13% to 50%. The structural basis for the crossreactivity of our lymphoma-derived private anti-ids is as yet not known. However, the reactivity of certain anti-ids with both kappa- and lambda-expressing tumors suggests a biased use of V gene segments in these crossreactive clones that is probably related to the VH gene. These data suggest the possibility that lymphoma may develop in a highly restricted pool of normal differentiated B cells or in B-cell subsets that express a limited repertoire of unmutated VH gene segments.

Examination of 64 translocations involving the major breakpoint region (mbr) of the BCL2 oncogene and the immunoglobulin heavy chain locus identified three short (14, 16, and 18 bp) segments within the mbr at which translocations occurred with very high frequency. Each of these clusters was associated with a 15-bp region of sequence homology, the principal one containing an octamer related to chi, the procaryotic activator of recombination. The presence of short deletions and N nucleotide additions at the breakpoints, as well as involvement of JH and DH coding regions, suggested that these sequences served as signals capable of interacting with the VDJ recombinase complex, even though no homology with the traditional heptamer/spacer/nonamer (IgRSS) existed. Furthermore, the BCL2 signal sequences were employed in a bidirectional fashion and could mediate recombination of one mbr region with another. Segments homologous to the BCL2 signal sequences flanked individual members of the XP family of diversity gene segments, which were themselves highly overrepresented in the reciprocal products (18q-) of BCL2 translocation. We propose that the chi-like signal sequences of BCL2 represent a distinct class of recognition sites for the recombinase complex, responsible for initiating interactions between regions of DNA separated by great distances, and that BCL2 translocation begins by a recombination event between mbr and DXP chi signals. Since recombinant joints containing chi, not IgRSS, occur in brain cells expressing RAG-1 (Matsuoka, M., F. Nagawa, K. Okazaki, L. Kingsbury, K. Yoshida, U. Muller, D. T. Larue, J. A. Winer, and H. Sakano. 1991. Science [Wash. DC]. 254:81; reference 1), we further suggest that the product of this gene could mediate both BCL2 translocation and the first step of normal DJ assembly through the creation of chi joints, rather than signal or coding joints.

Proto-oncogenes (H-ras-1 and L-myc) and tumor-suppressor gene (p53) loci have been implicated in lung carcinogenesis. DNA restriction fragment length polymorphisms at these gene loci are being evaluated in a case-control study as markers predictive of risk for cancer or of prognosis when cancer is present. The cases and controls had a cigarette-smoking history of 40 or more pack years or other abnormalities in pulmonary function tests, their ages were closely matched (64 years for cases and 61 years for controls) and the ratio of Caucasians to African Americans was close to unity (cases, 0.95:1.00, controls, 1.00:0.88). The H-ras-1 gene contains an insertion deletion polymorphism. Inheritance of rare H-ras-1 alleles, defined by MspI digestion, confers a relative risk for lung cancer of 2.0 (95% confidence interval, 0.5-7.3) for Caucasians and 3.2 (0.9-11.6) for African Americans (74 cases, 67 controls). The L-myc gene sequence has a restriction site (EcoR1) polymorphism between the second and third exons. Inheritance of restriction site-present alleles was reported to confer poor prognosis (presence of lymph node metastases) in Japanese lung cancer patients. This hypothesis was tested in both case-control study subjects (56 cases, 55 controls) and additional surgical cases (40), but no evidence was found to support the hypothesis in the U.S. population. The p53 gene is a tumor-suppressor gene that can encode either a proline or an arginine in the 72nd residue. No associations was found between the minor allele (proline) and diagnosis of lung cancer (76 cases, 68 controls).(ABSTRACT TRUNCATED AT 250 WORDS) Images FIGURE 1. A FIGURE 1. B FIGURE 1. C FIGURE 2. ABC FIGURE 2. DEF

BACKGROUND:Fibromyalgia syndrome (FMS), a common, chronic, widespread musculoskeletal pain disorder found in 2% of the general population and with a preponderance of 85% in females, has both genetic and environmental contributions. Patients and their parents have high plasma levels of the chemokines MCP-1 and eotaxin, providing evidence for both a genetic and an immunological/inflammatory origin for the syndrome (Zhang et al., 2008, Exp. Biol. Med. 233: 1171-1180). METHODS AND FINDINGS:In a search for a candidate gene affecting inflammatory pathways, among five screened in our patient samples (100 probands with FMS and their parents), we found 10 rare and one common alleles for MEFV, a gene in which various compound heterozygous mutations lead to Familial Mediterranean Fever (FMF). A total of 2.63 megabases of genomic sequence of the MEFV gene were scanned by direct sequencing. The collection of rare missense mutations (all heterozygotes and tested in the aggregate) had a significant elevated frequency of transmission to affecteds (p = 0.0085, one-sided, exact binomial test). Our data provide evidence that rare missense variants of the MEFV gene are, collectively, associated with risk of FMS and are present in a subset of 15% of FMS patients. This subset had, on average, high levels of plasma IL-1beta (p = 0.019) compared to FMS patients without rare variants, unaffected family members with or without rare variants, and unrelated controls of unknown genotype. IL-1beta is a cytokine associated with the function of the MEFV gene and thought to be responsible for its symptoms of fever and muscle aches. CONCLUSIONS:Since misregulation of IL-1beta expression has been predicted for patients with mutations in the MEFV gene, we conclude that patients heterozygous for rare missense variants of this gene may be predisposed to FMS, possibly triggered by environmental factors.

MicroRNAs regulate diverse cellular processes and play an integral role in cancer pathogenesis. Genomic variation within miRNA target sites may therefore be important sources for genetic differences in cancer risk. To investigate this possibility, we mapped HapMap single nucleotide polymorphisms (SNP) to putative miRNA recognition sites within genes dysregulated in estrogen receptor–stratified breast tumors and used local linkage disequilibirum patterns to identify high-ranking SNPs in the Cancer Genetic Markers of Susceptibility (CGEMS) breast cancer genome-wide association study for further testing. Two SNPs, rs1970801 and rs11097457, scoring in the top 100 from the CGEMS study, were in strong linkage disequilibrium with rs1434536, an SNP that resides within a miR-125b target site in the 3′ untranslated region of the bone morphogenic receptor type 1B (BMPR1B) gene encoding a transmembrane serine/threonine kinase. We validated the CGEMS association findings for rs1970801 in an independent cohort of admixture-corrected cases identified from families with multiple case histories. Subsequent association testing of rs1434536 for these cases and CGEMS controls with imputed genotypes supported the association. Furthermore, luciferase reporter assays and overexpression of miR-125b–mimics combined with quantitative reverse transcription-PCR showed that BMPR1B transcript is a direct target of miR-125b and that miR-125b differentially regulates the C and T alleles of rs1434536. These results suggest that allele-specific regulation of BMPR1B by miR-125b explains the observed disease risk. Our approach is general and can help identify and explain the mechanisms behind disease association for alleles that affect miRNA regulation. [Cancer Res 2009;69(18):7459–65]

SATB1 (special AT-rich sequence-binding protein-1) provides a key link between DNA loop organization, chromatin modification/remodeling, and association of transcription factors at matrix attachment regions (MARs). To investigate the role of SATB1 in cellular events, we performed a yeast two-hybrid screen that identified SUMO-1, Ubc9, and protein inhibitor of activated STAT (PIAS) family members as SATB1 interaction partners. These proteins, working in concert, enhanced SUMO conjugation to lysine-744 of SATB1. Overexpression of SUMO or PIAS in Jurkat cells, which express high levels of endogenous SATB1, exhibited enhanced caspase cleavage of this MAR-associating protein. Sumoylation-deficient SATB1 (SATB1(K744R)) failed to display the characteristic caspase cleavage pattern; however, fusion of SUMO in-frame to SATB1(K744R) restored cleavage. A SUMO-independent interaction of inactive caspase-6 and SATB1 was noted. A subset of total cellular SATB1 localized into promyelocytic leukemia nuclear bodies where enhanced SATB1 cleavage was detected subsequent to caspase activation. These results reveal a novel sumoylation-directed caspase cleavage of this key regulatory molecule. The role of regulated proteolysis of SATB1 may be to control transcription in immune cells during normal cell functions or to assist in efficient and rapid clearance of nonfunctional or potentially damaging immune cells.

Previously, a candidate gene linkage approach on brother pairs affected with prostate cancer identified a locus of prostate cancer susceptibility at D3S1234 within the fragile histidine triad gene (FHIT), a tumor suppressor that induces apoptosis. Subsequent association tests on 16 SNPs spanning approximately 381 kb surrounding D3S1234 in Americans of European descent revealed significant evidence of association for a single SNP within intron 5 of FHIT. In the current study, re-sequencing and genotyping within a 28.5 kb region surrounding this SNP further delineated the association with prostate cancer risk to a 15 kb region. Multiple SNPs in sequences under evolutionary constraint within intron 5 of FHIT defined several related haplotypes with an increased risk of prostate cancer in European-Americans. Strong associations were detected for a risk haplotype defined by SNPs 138543, 142413, and 152494 in all cases (Pearson's chi(2) = 12.34, df 1, P = 0.00045) and for the homozygous risk haplotype defined by SNPs 144716, 142413, and 148444 in cases that shared 2 alleles identical by descent with their affected brothers (Pearson's chi(2) = 11.50, df 1, P = 0.00070). In addition to highly conserved sequences encompassing SNPs 148444 and 152413, population studies revealed strong signatures of natural selection for a 1 kb window covering the SNP 144716 in two human populations, the European American (pi = 0.0072, Tajima's D = 3.31, 14 SNPs) and the Japanese (pi = 0.0049, Fay & Wu's H = 8.05, 14 SNPs), as well as in chimpanzees (Fay & Wu's H = 8.62, 12 SNPs). These results strongly support the involvement of the FHIT intronic region in an increased risk of prostate cancer.

A reexamination of human minisatellite (hypervariable) regions following the cloning and sequencing of the new minisatellite, VTR1.1, revealed that many of these structures possessed a strongly conserved copy of the chi-like octamer, GC[A/T]GG[A/T]GG. In oncogene translocations apparently created by aberrant VDJ recombinase activity, this VTR octamer was often found within a few bases of the breakpoint (p less than 10(-10)). Three bcl2 rearrangements which occurred within 2 bp of one another were located precisely adjacent to this consensus; it defined the 5' border of that oncogene's major breakpoint cluster. Several c-myc translocations also occurred within 2 bp of this sequence. While the appearance of a chi-like element in polymorphic minisatellite sequences is consistent with a role promoting either recombination or replication slippage, the existence of such elements at sites of somatic translocations suggests chi function in site-specific recombination, perhaps as a subsidiary recognition signal in immunoglobulin gene rearrangement. We discuss the implications of these observations for mechanisms by which oncogene translocations and minisatellite sequences are generated.

We have employed a large prostate cancer affected sibling pair cohort for candidate gene based linkage analyses. In this work we sought to enlarge a pre-existing cohort of CaP (Prostate Cancer) ASP with continued institutional recruitment of brothers affected with disease. We performed candidate gene based fine structure linkage analysis on approximately 2 dozen genes previously implicated in CaP risk. We also tested gene x gene interactions with a new paradigm based upon allele sharing enrichment. Our major finding was the localization of a susceptibility locus to intron 5 of the FHIT gene. By utilizing a combination of extensive mutation/single nucleotide polymorphism (SNP) discovery efforts in select disease cases in conjunction with linkage disequilibrium (LD) mapping and association testing we identified a SNP, rs760317, showing strong association with disease in affected brothers sharing 2 alleles identify by descent (IBD). The findings were published in 2005 and have recently been replicated by independent researchers in both a family-based Caucasian patient cohort and an African American patient cohort. Our efforts represent a significant accomplishment in the identification of a new gene associated with CaP risk as quite often promising initial linkage or association results fail to be replicated in independent studies. We continue our efforts today with the hope of finding the causative allele(s) in FHIT and it/their possible function using population genetic tools. This represents extreme challenges as the mechanistic basis for how disease alleles residing deep within the introns contribute to disease risk.

The study of genetic variation will promote our understanding of the differential predisposition to common diseases and variation in drug responses of individuals and ethnic populations. Such genetic variation is intrinsically structured into blocks of haplotypes in populations. Therefore, a comprehensive haplotype map based on the most abundant form of genetic variation, single nucleotide polymorphisms, will be useful. At the present time, however, our knowledge of the similarities and differences of haplotype structure among different ancestral populations is still inadequate.To determine whether common underlying haplotype patterns existed across ethnic populations, we analyzed data derived from African and European Americans for twenty-two genes spanning a total of 516 kb and the HapMap ENCODE data across 500 kb on chromosome 2p16.3 from three major world populations.We observed that strong pairwise linkage disequilibrium (LD) between SNPs selected from populations having African ancestry was highly conserved across other non-African populations. Common haplotypes described by these LD-selected SNPs demonstrated a simple evolutionary structure with up to three major frameworks, which were likely ancestral backgrounds upon which more recent mutations have been superimposed. Also, haplotype block boundaries defined in populations having African ancestry revealed completely concordant recombinant haplotypes across all populations, providing a consistent definition of block structure. Finally, a large fraction of regulatory polymorphisms described in the literature appeared to tag these conserved haplotype frameworks, strongly suggesting their significance for disease association and pharmacogenetic studies.

There is some evidence that women with a higher number of CAG repeat lengths on the androgen receptor (AR) gene have increased breast cancer risk. We evaluated the association between AR-CAG repeat length and mammographic density, a strong breast cancer risk factor, in 404 African-American and Caucasian breast cancer patients. In postmenopausal estrogen progestin therapy users, carriers of the less active AR-CAG had statistically significantly higher mean percentage of density (41.4%) than carriers of the more active AR-CAG (25.7%; P = 0.04). Our results raise the question of whether the number of AR-CAG repeats predicts breast cancer risk in estrogen progestin therapy users.

Rare alleles at the HRAS1 variable number of tandem repeats (VNTRs) locus have been implicated in breast cancer risk. We assessed the association of rare HRAS1 alleles and breast cancer in a case-control study nested within the Nurses' Health Study cohort. Using PCR-based methods, 717 incident breast cancer cases and 798 controls were genotyped for the HRAS1 VNTRs. The prevalence of the rare alleles in breast cancer cases was not different compared with controls (10.7 versus 12.0%, respectively; P = 0.45, two-sided Cochran-Mantel-Haenzel chi(2) test). There was no evidence that women heterozygous (multivariate odds ratio, 0.97; 95% confidence interval, 0.73-1.27) or homozygous (multivariate odds ratio, 0.83; 95% confidence interval, 0.32-2.14) for rare alleles were at an increased risk of breast cancer or that a positive gene-dose effect existed. The results did not vary by menopausal status. Although as a group the rare alleles were not associated with breast cancer, one class of rare alleles between the common alleles of a3 and a4 was associated with a significantly increased risk. These results suggest that there is no overall association between rare alleles of the HRAS1 VNTR and breast cancer.

PML bodies play an important role in multiple pathways of growth control, such as transformation suppression, apoptosis, and Ras-induced senescence. However, the molecular and biological bases for these physiological phenomena are not well understood. The findings that viruses transcribe their genomes adjacent to PML bodies and that nascent RNA accumulates at their periphery have suggested that PML bodies are transcription-permissible domains. To investigate the role of PML bodies in regulation of cell transformation and apoptosis-related gene transcription, we employed the immuno-FISH method to examine the relationship between PML bodies and the TP53 and BCL2 gene loci. PML bodies were found to localize specifically with the TP53 locus in about 50% of Jurkat interphase nuclei, but never in proximity with the BCL2 locus. We speculate that PML bodies may interact directly with the TP53 DNA sequence to regulate TP53 gene expression.

A yeast two-hybrid screen of a Jurkat (T cell) derived cDNA library, using SATB1 (a matrix attachment region binding protein) as the bait, yielded four independent isolates of a novel variant of the DNA directed RNA polymerase II, subunit 11 (RPB11). Absence of lysine-17 from the amino terminus of this variant cannot be explained by alternative mRNA splicing. Instead, allele-specific PCR, combined with GenBank database searches, suggests that a recent gene duplication event has resulted in distinct loci encoding three variant forms of RPB11. Differential splicing of mRNA transcripts accounts for unique carboxy termini among the RPB11 proteins. The exclusive association of SATB1 with one form of RPB11 is influenced primarily by the N-terminal amino acid disparity, as deletion of the entire C terminus does not alter interaction affinity. Association of RPB11 with SATB1 maps between amino acids 58 and 222 of SATB1, a region that includes a PDZ-like dimerization motif.

Translocation of the BCL2 gene from chromosome 18 to chromosome 14 results in constitutive expression of the gene. We have recently demonstrated that the major breakpoint region (mbr) of BCL2, which is implicated in 70% of t(14;18) translocations present in human follicular lymphoma, is a matrix attachment region. Since these regions are implicated in control of both transcription and replication, we wished to determine whether BCL2 translocation was also accompanied by changes in replication timing of the translocated allele. Using both fluorescence in situ hybridization and allele-specific PCR, we have demonstrated that the translocated allele replicates at the G1/S boundary, while the wildtype allele continues to replicate as usual in mid-S phase. These differences are accompanied by allele-specific changes in BCL2 expression. Since the net structural effect of t(14;18) translocations within the mbr is to disrupt the BCL2 MAR and replace it with the IGH MARs located just downstream of each breakpoint, we conclude that MAR exchange is a significant, selectable outcome of these translocations. We propose that subsequent changes of replication and transcriptional patterns for the translocated BCL2 allele result from this exchange and represent important early steps in lymphomagenesis.

Inherited mutations in the cystatin B gene ( CSTB ) are responsible for progressive myoclonus epilepsy type 1 (EPM1; MIM 254800). This autosomal recessive disease is characterized by variable progression to mental retardation, dementia and ataxia. The majority of EPM1 alleles identified to date contain expansions of a dodecamer repeat located upstream of the transcription start site of the CSTB gene. Normal alleles contain two or three copies of the repeat, whereas pathogenic alleles contain >40 repeats. We examined the meiotic stability of pathogenic, expanded EPM1 alleles from 17 EPM1 families by employing a fluorescence-based PCR-based genotyping assay capable of detecting single dodecamer repeat unit differences on an automated DNA sequencer. We followed 74 expanded allele transmissions to 30 affected individuals and 22 carriers. Thirty-five of 74 expanded allele transmissions demonstrated either contraction or expansion of the minisatellite, typically by a single repeat unit. Thus expanded alleles of the EPM1 minisatellite demonstrate a mutation rate of 47%, the highest yet observed for pathogenetic alleles of a human minisatellite.

Fibromyalgia syndrome (FMS) is a chronic musculoskeletal pain disorder affecting 2% to 5% of the general population. Both genetic and environmental factors may be involved. To ascertain in an unbiased manner which genes play a role in the disorder, we performed complete exome sequencing on a subset of FMS patients. Out of 150 nuclear families (trios) DNA from 19 probands was subjected to complete exome sequencing. Since >80,000 SNPs were found per proband, the data were further filtered, including analysis of those with stop codons, a rare frequency (

A substantial proportion of familial colorectal cancer (CRC) is not a consequence of known susceptibility loci, such as mismatch repair (MMR) genes, supporting the existence of additional loci. To identify novel CRC loci, we conducted a genome-wide linkage scan in 356 white families with no evidence of defective MMR (i.e., no loss of tumor expression of MMR proteins, no microsatellite instability (MSI)-high tumors, or no evidence of linkage to MMR genes). Families were ascertained via the Colon Cancer Family Registry multi-site NCI-supported consortium (Colon CFR), the City of Hope Comprehensive Cancer Center, and Memorial University of Newfoundland. A total of 1,612 individuals (average 5.0 per family including 2.2 affected) were genotyped using genome-wide single nucleotide polymorphism linkage arrays; parametric and non-parametric linkage analysis used MERLIN in a priori-defined family groups. Five lod scores greater than 3.0 were observed assuming heterogeneity. The greatest were among families with mean age of diagnosis less than 50 years at 4q21.1 (dominant HLOD = 4.51, α = 0.84, 145.40 cM, rs10518142) and among all families at 12q24.32 (dominant HLOD = 3.60, α = 0.48, 285.15 cM, rs952093). Among families with four or more affected individuals and among clinic-based families, a common peak was observed at 15q22.31 (101.40 cM, rs1477798; dominant HLOD = 3.07, α = 0.29; dominant HLOD = 3.03, α = 0.32, respectively). Analysis of families with only two affected individuals yielded a peak at 8q13.2 (recessive HLOD = 3.02, α = 0.51, 132.52 cM, rs1319036). These previously unreported linkage peaks demonstrate the continued utility of family-based data in complex traits and suggest that new CRC risk alleles remain to be elucidated.

BACKGROUND The role of mutations in protooncogenes and their regulatory sequences in the pathogenesis of cancer is under close scrutiny. Minisatellites are unstable repetitive sequences of DNA that are present throughout the human genome. The highly polymorphic HRAS1 minisatellite locus just downstream from the protooncogene H-ras-1 consists of four common progenitor alleles and several dozen rare alleles, which apparently derive from mutations of the progenitors. We previously observed an association of the rare mutant alleles with many forms of cancer, and we undertook the present study to pursue this observation further. METHODS We conducted a case-control study, typing 736 HRAS1 alleles from patients with cancer and 652 from controls by Southern blotting of leukocyte DNA. We also carried out a meta-analysis of this study and 22 other published studies, estimating the relative risk of cancer (such as bladder, breast, or colorectal cancer) when one of the rare HRAS1 alleles was present. RESULTS Both the present case-control study (odds ratio, 1.83; 95 percent confidence interval, 1.28 to 2.67; P = 0.002) and the present study combined with our previous study (odds ratio, 2.07; 95 percent confidence interval, 1.47 to 2.92; P < 0.001), as well as the meta-analysis of all 23 studies (odds ratio, 1.93; 95 percent confidence interval, 1.63 to 2.30; chi-square = 57.58; P < 0.001), replicated our original finding and demonstrated a significant association of rare HRAS1 alleles with cancer. We found significant associations for four types of cancer: carcinomas of the breast, colorectum, and urinary bladder and acute leukemia. We also identified suggestive but not statistically significant associations for cancers of the lung and prostate and for non-Hodgkin's lymphoma. CONCLUSIONS Mutant alleles of the HRAS1 minisatellite locus represent a major risk factor for common types of cancer. Although the relative risk associated with the presence of one rare allele is moderate, the aggregate prevalence of one rare allele is moderate, the aggregate prevalence of this class of mutant alleles implies an extremely important attributable risk: 1 in 11 cancers of the breast, colorectum, and bladder.

We have recently shown that constitutively expressed members of the rel /NF-κB family of transcription factors bind the HRAS1 minisatellite, VTR HRAS1 . We now report that, like other NF-κB binding sites, VTR HRAS1 displays pleiotropic transcriptional regulatory activity that is promoter- and cell-type-specific. Both enhancement and suppression are restricted to the human bladder carcinoma cell line EJ, in which we have previously defined a unique form of NF-κB p50. We also observe allelic variation in functional activity: the rare a2.1 allele, one member of a class of HRAS1 alleles overrepresented in the genomes of cancer patients, possesses twofold greater enhancer activity than the low-risk alleles, a0.1, a1, and a2. Finally, VTR HRAS1 enhancer activity is upregulated by the adenovirus E1A 13S gene product, demonstrating the potential of the minisatellite for influencing gene expression through several distinct interactions with the transcriptional apparatus.

The 50bp repeat unit of the minisatellite within the DH-JH interval of the human immunoglobulin heavy chain locus binds a nuclear factor present in a wide variety of cell types. The binding site contains the myc/HLH motif, CACGTG, and represents a 15 of 17 base match for the USF/MLTF binding site adjacent to the adenovirus major late promoter (MLP). Unlike the USF/MLTF site, the IGH minisatellite possesses no enhancer activity. However, it can significantly suppress, in cis and in trans, USF-site-mediated transcriptional activation of the MLP. In murine myeloma cells, the IGH minisatellite can suppress, in trans, MLP activation by the murine heavy chain gene enhancer, E mu. These activities potentially represent a DNA-based form of squelching.

We have cloned and characterized a 55-kb region of DNA surrounding HRAS1. It contains a cluster of two, and possibly three, genes associated with CpG islands within the 32 kb immediately upstream of HRAS1. We have sequenced cDNAs representing one of these genes, provisionally designated HRC1. The locus, which is located 29 kb upstream of HRAS1, is divergently transcribed. HRC1 cDNA probe recognizes fragments on Southern blots of DNA from other vertebrate species. In human DNA, multiple homologous fragments are detected in addition to the predicted ones containing HRC1. Therefore, this locus may represent a member of an evolutionarily conserved gene family. HRC1 expression is upregulated with HRAS1 in the EJ bladder carcinoma cell line, suggesting the possibility of coordinate regulation. The deduced translational product of the longest open reading frame (1119 nucleotides, 373 amino acids) predicts a protein with regions rich in glutamine and proline and a region similar to the helix-loop-helix motif adjacent to a carboxy-terminal leucine zipper dimerization motif with four heptad repeats. Alternate splicing of terminal exons occurs, resulting in the truncation of one proline-rich domain and preservation of the leucine zipper. Thus, a biologically important region of chromosome 11p consists of a gene cluster. At least one of these genes, in addition to HRAS1, may be involved in regulation of cell growth or differentiation.

The 28 base pair repeat unit of a minisatellite 1000 bp downstream from the human HRAS1 gene (VTRHRAS1) bound four proteins (p45, p50, p72 and p85) in nuclear extracts from a variety of human cell lines which were indistinguishable from several members of the rel/NF-kappa B family of transcriptional regulatory factors. VTRHRAS1 bound the constitutively expressed, but not the inducible, forms of these proteins. Analysis of partially purified binding factors from different cell lines demonstrated qualitative differences in the p50 subunit; phosphocellulose fractionation also revealed considerable heterogeneity in the p72 and p85 subunits. These results suggest the possibility that the HRAS1 minisatellite, in serving as a tandem array of rel/NF-kappa B binding sites, may function in the transcriptional regulation of HRAS1 and nearby genes.

Restriction fragment length polymorphism analyses of the Ha-ras-1 proto-oncogene were undertaken in white and black populations residing in the Baltimore-Washington metropolitan area to address whether specific rare alleles of the Ha-ras-1 proto-oncogene locus vary in their distribution among different racial groups. High-molecular-weight genomic DNA samples from the lungs of 80 lung cancer patients and 92 accident victims were digested with appropriate restriction enzymes and subjected to Southern analysis using the 6.6-kb BamHI human Ha-ras-1 recombinant fragment from the plasmid pEC. Thirty allelomorphs of different sizes were detected among the 172 study subjects. An association was observed between race and specific alleles. Rare alleles were more frequent in black cancer patients and trauma victims than in whites. Within each racial category, lung cancer patients had an excess of rare alleles. These data indicate the importance of controlling for racial variation when designing studies to determine human cancer risk factors.

Variation in gene expression may give rise to a significant fraction of inter-individual phenotypic variation. Studies searching for the underlying genetic controls for such variation have been conducted in model organisms and humans in recent years. In our previous effort of assessing conserved underlying haplotype patterns across ethnic populations, we constructed common haplotypes using SNPs having conserved linkage disequilibrium (LD) across ethnic populations. These common haplotypes cluster into a simple evolutionary structure based on their frequencies, defining only up to three conserved clusters termed 'haplotype frameworks'. One intriguing preliminary finding was that a significant portion of reported variants strongly associated with cis-regulation tags these globally conserved haplotype frameworks. Here we expand the investigation by collecting genes showing stringently determined cis-association between genotypes and expression phenotypes from major studies. We conducted phylogenetic analysis of current major haplotypes along with the corresponding haplotypes derived from chimpanzee reference sequences. Our analysis reveals that, for the vast majority of such cis-regulatory genes, the tagging SNPs showing the strongest association also tag the haplotype lineages directly separated from ancestry, inferred from either chimpanzee reference sequences or the allele frequency-derived haplotype frameworks, suggesting that the differentially expressed phenotypes were evolved relatively early in human history. Such evolutionary signatures provide keys for a more effective identification of globally-conserved candidate regulatory haplotypes across human genes in future epidemiologic and pharmacogenetic studies.

New methods are needed for the identification of pathogenic alleles of candidate genes that may increase cancer susceptibility. Such risk alleles are expected to be of low penetrance and may act alone or modify the effects of other genes. We have developed a method that enriches for pathogenetic disease variants contingent on gene−gene interactions. Candidate gene pairs are chosen based on previous evidence demonstrating genetic or biochemical interaction. Utilizing a cohort of sibling pairs affected with breast cancer, we tested this paradigm by identifying risk variants of cell cycle control gene CDKN1A by means of interactions with TP53 and BRCA1. This approach, which we call disease association by locus stratification, first stratified affected pairs based on sharing of both alleles at a microsatellite marker linked to CDKN1A. The second stratification was based on microsatellite marker sharing at TP53 or BRCA1. We identified subsets of affected pairs sharing both alleles at both loci as screening targets. Utilizing this approach, we were able to enrich for two noncoding disease haplotypes of CDKN1A by virtue of BRCA1 interactions. We defined each haplotype by single-nucleotide polymorphisms at two positions potentially important in CDKN1A transcriptional activation by both p53 and BRCA1p. Our results indicated that an approach based on allele sharing and gene−gene interactions will be valuable not only in identifying risk alleles but also in elucidating their mechanism of action.

A rapid and exciting accumulation of data about cellular oncogenes in human tumors has resulted from convergent research on DNA-mediated gene transfer, retroviruses, and tumor cytogenetics. Such work promises to increase our understanding of the genetic events that predispose to, and result in, malignant disease. This knowledge may quickly find clinical application in tumor classification and prediction of risk. Ultimately, therapeutic benefits may be achieved as we begin to explore the mechanisms by which transforming gene products act to defeat the normal regulatory processes of cells.

A highly polymorphic locus associated with the variable tandem repetition of a 35 bp consensus sequence was mapped to chromosome 10, band q26. Examination of leukocyte DNA from a cancer patient revealed the twenty-fold amplification of one allelic fragment of this locus, while the other allelic fragment demonstrated a normal copy number. In another patient, Southern blotting of leukocyte DNA detected the deletion of the 3'-flanking region from one tandem repeat allele. These results indicate that variable tandem repeats may mark highly unstable regions of DNA in the human genome which can be altered by changes more extensive than simple tandem repeat variation.

White blood cell DNA from cancer patients and DNA from tumours and tumour-derived cell lines frequently demonstrated allelic restriction fragments of the Harvey ras oncogene locus not found in the unaffected population. The presence of such unusual alleles may be linked to susceptibility to cancer.

A variable tandem repeat (VTR) is responsible for the hyperallelism one kilobase 3' to the human c-Ha-ras-1 (Ha-ras) gene. Thirty-two distinct restriction fragments, comprising 3 allelic classes by frequency of occurrence, have thus far been detected in a sample size of approximately 800 caucasians. Rare Ha-ras alleles, 21 in all, are almost exclusively confined to the genomes of cancer patients (p less than 0.001). From our data we have computed the relative cancer risk associated with possession of a rare Ha-ras allele to be 27. To understand the molecular basis for this phenomenon, we have begun to clone Ha-ras fragments from nontumor DNA of cancer patients. We report here the weak activation, as detected by transfection and transformation of NIH 3T3 mouse cells, of two Ha-ras genes which were obtained from lymphocyte DNA of a melanoma patient. We have mapped the regions that confer this transforming activity to the fragment containing the VTR in one Ha-ras clone and the fragment containing gene coding sequences in the other. Images FIGURE 1.