Colocalisation of predicted exonic splicing enhancers in BRCA2 with reported sequence variants.

Disruption of the breast cancer susceptibility gene BRCA2 is associated with increased risk of developing breast and ovarian cancer. Over 1800 sequence changes in BRCA2 have been reported, although for many the pathogenicity is unclear. Classifying these changes remains a challenge, as they may disrupt regulatory sequences as well as the primary protein coding sequence. Sequence changes located in the splice site consensus sequences often disrupt splicing, however sequence changes located within exons are also able to alter splicing patterns. Unfortunately, the presence of these exonic splicing enhancers (ESEs) and the functional effect of variants within ESEs it is currently difficult to predict. We have previously developed a method of predicting which sequence changes within exons are likely to affect splicing, using BRCA1 as an example. In this paper, we have predicted ESEs in BRCA2 using the web-based tool ESEfinder and incorporated the same series of filters (increased threshold, 125 nt limit and evolutionary conservation of the motif) in order to identify predicted ESEs that are more likely to be functional. Initially 1114 ESEs were predicted for BRCA2, however after all the filters were included, this figure was reduced to 31, 3% of the original number of predicted ESEs. Reported unclassified sequence variants in BRCA2 were found to colocalise to 55% (17/31) of these conserved ESEs, while polymorphisms colocalised to 0 of the conserved ESEs. In summary, we have identified a subset of unclassified sequence variants in BRCA2 that may adversely affect splicing and thereby contribute to BRCA2 disruption.

Genetic and histopathologic evaluation of BRCA1 and BRCA2 DNA sequence variants of unknown clinical significance.

Classification of rare missense variants as neutral or disease causing is a challenge and has important implications for genetic counseling. A multifactorial likelihood model for classification of unclassified variants in BRCA1 and BRCA2 has previously been developed, which uses data on co-occurrence of the unclassified variant with pathogenic mutations in the same gene, cosegregation of the unclassified variant with affected status, and Grantham analysis of the fit between the missense substitution and the evolutionary range of variation observed at its position in the protein. We have further developed this model to take into account relevant features of BRCA1- and BRCA2-associated tumors, such as the characteristic histopathology and immunochemical profiles associated with pathogenic mutations in BRCA1, and the fact that approximately 80% of tumors from BRCA1 and BRCA2 carriers undergo inactivation of the wild-type allele by loss of heterozygosity. We examined 10 BRCA1 and 15 BRCA2 unclassified variants identified in Australian, multiple-case breast cancer families. By a combination of genetic, in silico, and histopathologic analyses, we were able to classify one BRCA1 variant as pathogenic and six BRCA1 and seven BRCA2 variants as neutral. Five of these neutral variants were also found in at least 1 of 180 healthy controls, suggesting that screening a large number of appropriate controls might be a useful adjunct to other methods for evaluation of unclassified variants.

Identification of BRCA1 missense substitutions that confer partial functional activity: potential moderate risk variants?

INTRODUCTION: Many of the DNA sequence variants identified in the breast cancer susceptibility gene BRCA1 remain unclassified in terms of their potential pathogenicity. Both multifactorial likelihood analysis and functional approaches have been proposed as a means to elucidate likely clinical significance of such variants, but analysis of the comparative value of these methods for classifying all sequence variants has been limited. METHODS: We have compared the results from multifactorial likelihood analysis with those from several functional analyses for the four BRCA1 sequence variants A1708E, G1738R, R1699Q, and A1708V. RESULTS: Our results show that multifactorial likelihood analysis, which incorporates sequence conservation, co-inheritance, segregation, and tumour immunohistochemical analysis, may improve classification of variants. For A1708E, previously shown to be functionally compromised, analysis of oestrogen receptor, cytokeratin 5/6, and cytokeratin 14 tumour expression data significantly strengthened the prediction of pathogenicity, giving a posterior probability of pathogenicity of 99%. For G1738R, shown to be functionally defective in this study, immunohistochemistry analysis confirmed previous findings of inconsistent 'BRCA1-like' phenotypes for the two tumours studied, and the posterior probability for this variant was 96%. The posterior probabilities of R1699Q and A1708V were 54% and 69%, respectively, only moderately suggestive of increased risk. Interestingly, results from functional analyses suggest that both of these variants have only partial functional activity. R1699Q was defective in foci formation in response to DNA damage and displayed intermediate transcriptional transactivation activity but showed no evidence for centrosome amplification. In contrast, A1708V displayed an intermediate transcriptional transactivation activity and a normal foci formation response in response to DNA damage but induced centrosome amplification. CONCLUSION: These data highlight the need for a range of functional studies to be performed in order to identify variants with partially compromised function. The results also raise the possibility that A1708V and R1699Q may be associated with a low or moderate risk of cancer. While data pooling strategies may provide more information for multifactorial analysis to improve the interpretation of the clinical significance of these variants, it is likely that the development of current multifactorial likelihood approaches and the consideration of alternative statistical approaches will be needed to determine whether these individually rare variants do confer a low or moderate risk of breast cancer.

Variation in the RAD51 gene and familial breast cancer.

INTRODUCTION: Human RAD51 is a homologue of the Escherichia coli RecA protein and is known to function in recombinational repair of double-stranded DNA breaks. Mutations in the lower eukaryotic homologues of RAD51 result in a deficiency in the repair of double-stranded DNA breaks. Loss of RAD51 function would therefore be expected to result in an elevated mutation rate, leading to accumulation of DNA damage and, hence, to increased cancer risk. RAD51 interacts directly or indirectly with a number of proteins implicated in breast cancer, such as BRCA1 and BRCA2. Similar to BRCA1 mice, RAD51-/- mice are embryonic lethal. The RAD51 gene region has been shown to exhibit loss of heterozygosity in breast tumours, and deregulated RAD51 expression in breast cancer patients has also been reported. Few studies have investigated the role of coding region variation in the RAD51 gene in familial breast cancer, with only one coding region variant--exon 6 c.449G>A (p.R150Q)--reported to date. METHODS: All nine coding exons of the RAD51 gene were analysed for variation in 46 well-characterised, BRCA1/2-negative breast cancer families using denaturing high-performance liquid chromatography. Genotyping of the exon 6 p.R150Q variant was performed in an additional 66 families. Additionally, lymphoblastoid cell lines from breast cancer patients were subjected to single nucleotide primer extension analysis to assess RAD51 expression. RESULTS: No coding region variation was found, and all intronic variation detected was either found in unaffected controls or was unlikely to have functional consequences. Single nucleotide primer extension analysis did not reveal any allele-specific changes in RAD51 expression in all lymphoblastoid cell lines tested. CONCLUSION: Our study indicates that RAD51 is not a major familial breast cancer predisposition gene.

Evolutionary conservation analysis increases the colocalization of predicted exonic splicing enhancers in the BRCA1 gene with missense sequence changes and in-frame deletions, but not polymorphisms.

INTRODUCTION: Aberrant pre-mRNA splicing can be more detrimental to the function of a gene than changes in the length or nature of the encoded amino acid sequence. Although predicting the effects of changes in consensus 5' and 3' splice sites near intron:exon boundaries is relatively straightforward, predicting the possible effects of changes in exonic splicing enhancers (ESEs) remains a challenge. METHODS: As an initial step toward determining which ESEs predicted by the web-based tool ESEfinder in the breast cancer susceptibility gene BRCA1 are likely to be functional, we have determined their evolutionary conservation and compared their location with known BRCA1 sequence variants. RESULTS: Using the default settings of ESEfinder, we initially detected 669 potential ESEs in the coding region of the BRCA1 gene. Increasing the threshold score reduced the total number to 464, while taking into consideration the proximity to splice donor and acceptor sites reduced the number to 211. Approximately 11% of these ESEs (23/211) either are identical at the nucleotide level in human, primates, mouse, cow, dog and opossum Brca1 (conserved) or are detectable by ESEfinder in the same position in the Brca1 sequence (shared). The frequency of conserved and shared predicted ESEs between human and mouse is higher in BRCA1 exons (2.8 per 100 nucleotides) than in introns (0.6 per 100 nucleotides). Of conserved or shared putative ESEs, 61% (14/23) were predicted to be affected by sequence variants reported in the Breast Cancer Information Core database. Applying the filters described above increased the colocalization of predicted ESEs with missense changes, in-frame deletions and unclassified variants predicted to be deleterious to protein function, whereas they decreased the colocalization with known polymorphisms or unclassified variants predicted to be neutral. CONCLUSION: In this report we show that evolutionary conservation analysis may be used to improve the specificity of an ESE prediction tool. This is the first report on the prediction of the frequency and distribution of ESEs in the BRCA1 gene, and it is the first reported attempt to predict which ESEs are most likely to be functional and therefore which sequence variants in ESEs are most likely to be pathogenic.

Identification of SQSTM1 mutations in familial Paget's disease in Australian pedigrees.

We have conducted a genome-wide scan on a pedigree containing 372 adult members, of whom 49 have PDB. In the present study, we report linkage of a large pedigree to the PDB3 region on chromosome 5q35-qter with a peak multipoint LOD score of 6.77. Sequestosome 1 (SQSTM/p62) has been identified as the causative PDB gene in this region. Six sequestosome 1 mutations have been described to date. Four mutations have been identified in exon 8, 1210delT and 1215delC both resulting in premature stop codon at amino acid 394, 1215C to T (P392L), 1224insT (E396X), one mutation in exon 7, 1200C to T (P387L) and a G to A splice junction mutation at IVS7+1. These mutations cluster in the C terminus of the protein and are predicted to disrupt the ubiquitin binding properties of sequestosome 1. Sequence analysis of the gene encoding sequestosome 1 revealed a single base pair deletion (1215delC) segregating with the majority of affected members in the pedigree. This deletion introduces a stop codon at position 394, resulting in premature termination of the protein (L394X) and loss of the ubiquitin-associated binding domain. Screening of affected members from 10 further PDB families identified the previously reported P392L mutation in one family. No SQSTM1/p62 coding mutations were found in the remaining 9 families or in 113 age-matched controls.

Prediction of BRCA1 and BRCA2 mutation status using post-irradiation assays of lymphoblastoid cell lines is compromised by inter-cell-line phenotypic variability.

BACKGROUND AND PURPOSE: Assays to determine the pathogenicity of unclassified sequence variants in disease-associated genes include the analysis of lymphoblastoid cell lines (LCLs). We assessed the ability of several assays of LCLs to distinguish carriers of germline BRCA1 and BRCA2 gene mutations from mutation-negative controls to determine their utility for use in a diagnostic setting. MATERIALS AND METHODS: Post-ionising radiation cell viability and micronucleus formation, and telomere length were assayed in LCLs carrying BRCA1 or BRCA2 mutations, and in unaffected mutation-negative controls. RESULTS: Post-irradiation cell viability and micronucleus induction assays of LCLs from individuals carrying pathogenic BRCA1 mutations, unclassified BRCA1 sequence variants or wildtype BRCA1 sequence showed significant phenotypic heterogeneity within each group. Responses were not consistent with predicted functional consequences of known pathogenic or normal sequences. Telomere length was also highly heterogeneous within groups of LCLs carrying pathogenic BRCA1 or BRCA2 mutations, and normal BRCA1 sequences, and was not predictive of mutation status. CONCLUSION: Given the significant degree of phenotypic heterogeneity of LCLs after gamma-irradiation, and the lack of association with BRCA1 or BRCA2 mutation status, we conclude that the assays evaluated in this study should not be used as a means of differentiating pathogenic and non-pathogenic sequence variants for clinical application. We suggest that a range of normal controls must be included in any functional assays of LCLs to ensure that any observed differences between samples reflect the genotype under investigation rather than generic inter-individual variation.

A genomewide search for type 2 diabetes-susceptibility genes in indigenous Australians.

The prevalence of type 2 diabetes among Australian residents is 7.5%; however, prevalence rates up to six times higher have been reported for indigenous Australian communities. Epidemiological evidence implicates genetic factors in the susceptibility of indigenous Australians to type 2 diabetes and supports the hypothesis of the "thrifty genotype," but, to date, the nature of the genetic predisposition is unknown. We have ascertained clinical details from a community of indigenous Australian descent in North Stradbroke Island, Queensland. In this population, the phenotype is characterized by severe insulin resistance. We have conducted a genomewide scan, at an average resolution of 10 cM, for type 2 diabetes-susceptibility genes in a large multigeneration pedigree from this community. Parametric linkage analysis undertaken using FASTLINK version 4.1p yielded a maximum two-point LOD score of +2.97 at marker D2S2345. Multipoint analysis yielded a peak LOD score of +3.9 <1 cM from marker D2S2345, with an 18-cM 3-LOD support interval. Secondary peak LOD scores were noted on chromosome 3 (+1.8 at recombination fraction [theta] 0.05, at marker D3S1311) and chromosome 8 (+1.77 at theta=0.0, at marker D8S549). These chromosomal regions are likely to harbor novel susceptibility genes for type 2 diabetes in the indigenous Australian population.