Next Article in Journal
Heat Shock Factor 1 Depletion Sensitizes A172 Glioblastoma Cells to Temozolomide via Suppression of Cancer Stem Cell-Like Properties
Next Article in Special Issue
Hotspot Selective Preference of the Chimeric Sequences Formed in Multiple Displacement Amplification
Previous Article in Journal
The Temporal Pattern, Flux, and Function of Autophagy in Spinal Cord Injury
Previous Article in Special Issue
Next-Generation Sequencing in Oncology: Genetic Diagnosis, Risk Prediction and Cancer Classification
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 18
Issue 2
10.3390/ijms18020467
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Exome Sequencing in a Family with Luminal-Type Breast Cancer Underpinned by Variation in the Methylation Pathway

by
Nicole Van der Merwe
1,
Armand V. Peeters
1,
Fredrieka M. Pienaar
2,
Juanita Bezuidenhout
1,
Susan J. Van Rensburg
3 and
Maritha J. Kotze
3,4,*
1
Division of Anatomical Pathology, Department of Pathology, Faculty of Medicine and Health Sciences, Stellenbosch University, Tygerberg 7500, South Africa
2
GVI Oncology, Mediclinic Panorama, Cape Town 8000, South Africa
3
Division of Chemical Pathology, Department of Pathology, Faculty of Medicine and Health Sciences, Stellenbosch University, Tygerberg 7500, South Africa
4
National Health Laboratory Service, Tygerberg Hospital, Tygerberg 7500, South Africa
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2017, 18(2), 467; https://doi.org/10.3390/ijms18020467
Submission received: 2 January 2017 / Revised: 31 January 2017 / Accepted: 10 February 2017 / Published: 22 February 2017
(This article belongs to the Special Issue Next-Generation Sequencing for Clinical Application)
Browse Figures
Versions Notes

Abstract

:
Panel-based next generation sequencing (NGS) is currently preferred over whole exome sequencing (WES) for diagnosis of familial breast cancer, due to interpretation challenges caused by variants of uncertain clinical significance (VUS). There is also no consensus on the selection criteria for WES. In this study, a pathology-supported genetic testing (PSGT) approach was used to select two BRCA1/2 mutation-negative breast cancer patients from the same family for WES. Homozygosity for the MTHFR 677 C>T mutation detected during this PSGT pre-screen step was considered insufficient to cause bilateral breast cancer in the index case and her daughter diagnosed with early-onset breast cancer (<30 years). Extended genetic testing using WES identified the RAD50 R385C missense mutation in both cases. This rare variant with a minor allele frequency (MAF) of <0.001 was classified as a VUS after exclusion in an affected cousin and extended genotyping in 164 unrelated breast cancer patients and 160 controls. Detection of functional polymorphisms (MAF > 5%) in the folate pathway in all three affected family members is consistent with inheritance of the luminal-type breast cancer in the family. PSGT assisted with the decision to pursue extended genetic testing and facilitated clinical interpretation of WES aimed at reduction of recurrence risk.

1. Introduction

Oncogenomics is a rapidly evolving integrative discipline that enables a novel approach to the molecular characterization of breast cancer pathology. Understanding the relationship between germline and tumour genetics is important in patients where early-onset familial breast cancer or differences in drug response and survival outcome cannot be explained by mutations in the major tumour suppressor genes, BRCA1 and BRCA2 [1,2,3]. In these circumstances, genome-scale microarray analysis and/or next generation sequencing (NGS) provide viable options for comprehensive testing in a well-defined target population. Application of these advanced molecular technologies in the diverse South African population evolved into a pathology-supported genetic testing (PSGT) service used to facilitate the move from single- to multi-gene testing and NGS [4]. Utilization of PSGT may extrapolate the capacity of NGS from genotype to phenotype.
NGS includes targeted sequencing of preselected gene panels as well as whole genome/exome sequencing, RNA sequencing, and bisulfite sequencing for assessment of DNA methylation patterns [5,6]. Although targeted sequencing panels are currently the method of choice in breast cancer diagnostics, we favour the use of whole exome sequencing (WES) preceded by PSGT to facilitate patient selection and clinical interpretation [7]. WES is not limited by the number of genes that can be evaluated for risk of familial breast cancer and tailored therapeutic intervention in a single test. However, the large amount of information that is generated throughout the coding region of the genome may outpace our ability to adequately interpret variant pathogenicity for clinical application. Comprehensive studies on the pathology of hereditary breast cancer is needed to characterize variants of uncertain clinical significance (VUS) increasingly uncovered by WES [8,9], which may use lower thresholds for breast cancer mutation screening compared to targeted BRCA1/2 mutation screening. There is currently no consensus on eligibility criteria for WES in BRCA1/2 mutation-negative familial breast cancer patients.
Grant et al. [10] demonstrated the value of a pre-screen algorithm using PSGT to reduce chemotherapy overtreatment in South African breast cancer patients eligible for microarray analysis using the 70-gene MammaPrint test. A similar approach was validated in the Microarray in Node-negative Disease may Avoid ChemoTherapy (MINDACT) trial, which provided level 1A evidence for the clinical utility of MammaPrint. Early-stage breast cancer patients classified as clinical high-risk based on tumour size, grade, nodal status, estrogen receptor (ER) and human epidermal growth factor receptor 2 (HER2) status, are eligible for microarray analysis [11]. The fact that approximately 50% of early-stage breast cancer patients at high risk of recurrence could safely forego chemotherapy based on a low-risk MammaPrint profile, confirmed the important role of genomics to stratify breast cancer patients into different treatment groups. The finding that susceptibility to a particular molecular subtype of breast cancer is associated with familial aggregation [8,12] furthermore emphasised the importance of MammaPrint used in conjunction with the 80-gene Blueprint microarray to subdivide tumours into the Luminal A, Luminal B, HER2-enriched and basal-types [13,14].
Germline genotyping interpreted in the context of tumour subtype may help to determine whether breast cancer is caused by genetic or environmental factors, or a combination of both, in order to direct treatment accordingly. Naushad et al. [15] found that the Luminal A subtype is associated with lack of family history of breast cancer, late age of onset and high body mass index (BMI), in contrast to the basal-type known to be strongly associated with high-penetrance BRCA1 familial breast cancer [8]. In line with tumour pathology, genes involved in breast cancer can broadly be categorized into those with high, moderate, and low penetrance of clinical manifestation. High-penetrance alleles typically confer lifetime risks of breast cancer above 50%, while moderate penetrance alleles confer risks greater than 20% and low-penetrance alleles increases of 10%–20% lifetime risk [16]. Risk above 50% is sufficient to recommend prophylactic surgery or magnetic resonance imaging screening in a familial context. Risk reduction strategies based on low-moderate penetrance genes are focused more on individual breast cancer survivors, where identification of biomarkers representing critical biological processes may be targeted at the gene–environment level to reduce the incidence of breast cancer recurrence.
It is estimated that known breast cancer susceptibility loci including BRCA1 and BRCA2 account for approximately 20% of familial risk, with the remaining variation caused by mutations with low-moderate penetrance that may be inherited in a subtype specific manner [17,18]. Failure to show significant survival differences between patients with familial and sporadic breast cancer [19] may partly be explained by polymorphic variation in drug metabolizing enzymes such as cytochrome P450 (CYP) [20], or modifier genes dependent on nutrient co-factors. Folate and other B-vitamins are important determinants of genome stability, yet might also contribute to the mechanism of action of anti-folates and fluoropyrimidines used in cancer treatment [21,22]. Xu et al. [22] identified one-carbon metabolism as a key pathway to be targeted for improved breast cancer survival.
Germline variation in folate pathway genes is associated with aberrant patterns of tumour pathology in cancer patients [23], with dysfunction of one-carbon metabolism found to be one of the three most important disease mechanisms in the pathogenesis of cancer [24]. Growing evidence of the complex interaction between methyl-related nutrients, alcohol and genetic variation in the methylene tetrahydrofolate reductase (MTHFR) gene support individualised dietary recommendations to reduce cancer risk in the general population [22,25]. These include substitution of low folate foods with folate-dense fruits and vegetables in countries where there is no mandatory folic acid fortification. Folate provides protection against carcinogenesis, especially in MTHFR mutation-positive cases, yet high doses of the synthetic form (folic acid) are suspected to promote tumour growth when administered after diagnosis of cancer. It is therefore reassuring that recent evidence obtained in postmenopausal patients with increased folate intake through diet and/or supplementation showed significantly higher circulating folate concentrations associated with improved survival after breast cancer diagnosis [26]. Since the negative effects of low intakes of the methyl-related nutrients with high intakes of alcohol are additive, Bailey [25] suggested that alcohol intake should be restricted to less than 15 g per day, if consumed.
The low-penetrance MTHFR 677 C>T mutation found to be associated with hormone receptor positive breast cancer [15], has previously been linked with increased BMI in South African patients diagnosed with major depressive disorder [27]. A low folate score, as reflected by the nutrition questionnaire used by Delport et al. [27], correlated with increased BMI and risk of depression in South African patients with the MTHFR 677 C>T mutation. Depression is a major comorbidity in breast cancer survivors and impacts on both the quality of life and treatment decision making. In this context, Van der Merwe et al. [28] recommended CYP2D6 genotyping in breast cancer patients who (1) receive tamoxifen treatment; and (2) are at increased recurrence risk due to a family history of breast cancer/BRCA-mutation positivity; or (3) are required to take potential competing antidepressants that may inhibit CYP2D6 enzyme function. CYP2D6 genotyping combined with the unique spectrum of BRCA1/2 founder mutations identified in the South African population [29,30], forms part of the pathology-supported pharmacogenetics test developed by van der Merwe [31] for improved clinical management of patients with breast cancer and associated comorbidities. Combining clinico-pathological features with pharmacogenetics in this chronic disease screening test is supported by the findings of Phipps et al. [32,33], who demonstrated subtype association with a family history of cancer and high BMI.
In this study, two BRCA1/2 mutation-negative breast cancer patients were selected from the same family for WES, followed by extended mutation screening in 164 unrelated cases stratified according to ER status. Lifestyle risk factors such as BMI were assessed in study participants due to the link with cardiovascular disease (CVD), predicted to equal or exceed recurrence risk of breast cancer in ER-positive postmenopausal women [34]. This approach can be extended into a PSGT pre-screen algorithm for WES, to help distinguish between familial breast cancer eligible for NGS and patients with lifestyle-related breast cancer who may reduce their recurrence risk by maintaining a healthy body weight, smoking cessation and alcohol restriction [35]. These lifestyle factors may act in combination with low-penetrance mutations underlying the folate-vitamin B12-methylation pathway evaluated as part of a CVD multi-gene assay applied in this study before commencing WES [4]. Testing for MTHFR 677 C>T as an invariant marker of increased risk for many chronic diseases associated with dysfunctional methylation may add value beyond biochemical analysis of homocysteine and folate levels, which are influenced by various environmental factors. For this reason, MTHFR is not only a feasible anticancer drug target, but may also be considered an important factor in prevention of cumulative disease risk across clinical entities.

2. Results

2.1. Pathology-Supported Genetic Testing (PSGT)

Figure 1 illustrates the PSGT approach used to select the index case for WES. Despite bilateral and early-onset breast cancer in the family, this patient tested negative for high-penetrance mutations in the BRCA1 and BRCA2 genes. She was initially referred for the MammaPrint test and could safely forego chemotherapy based on a low genomic risk for breast cancer recurrence. Microarray analysis was preceded by immunohistochemistry (IHC) testing indicating ER-positive (ER+), progesterone receptor-positive (PR+) and HER2-negative (HER2−) status. The index case remained disease-free since 2008 despite termination of tamoxifen after one year due to severe side effects, including. Use of the CVD multi-gene test excluded the most prevalent CYP2D6 null-allele in Caucasians (CYP2D6*4), while homozygosity of MTHFR 677 C>T was detected in both the index case and her daughter with the luminal-type breast cancer. It was uncertain whether this functional single nucleotide polymorphism (SNP) expressed in a high-risk environment, as reflected by high BMI (>27 kg/m2) and inadequate folate intake (score <14) in both cases, was sufficient for development of breast cancer in this family.

2.2. Comprehensive Cancer Panel Screen Using Whole Exome Sequencing

WES was used to search for potentially causative gene variants in 409 known cancer-related genes included in the Ion AmpliSeq™ Comprehensive Cancer Panel (ThermoFisher, Waltham, MA, USA). Rare variants in the double strand break repair protein (RAD50 R385C, Refseq mRNA accession NM_005732.3) and mucin 1 (MUC1 Q85E, Refseq mRNA accession NM_002456.5) genes were identified as potential causative mutations. These missense mutations were confirmed using Sanger sequencing (Figure 2 and Figure 3).
The RAD50 R385C mutation was detected in both the index case and her affected daughter, while the MUC1 Q85E mutation was only identified in the daughter. Both mutations were absent in a cousin of the index case diagnosed with breast cancer at the age of 42 years. Unlike RAD50, MUC1 has not been classified as a driver gene in breast cancer and was therefore not included in further analysis. Bioinformatic analysis of RAD50 R385C with a minor allele frequency (MAF) of <0.001 is shown in Table 1. Three of the seven software tools predicted a disease-causing damaging effect and one a moderate impact.

2.3. Extended Mutation Analysis

Extended mutation screening for RAD50 R385C using allele-specific real-time polymerase chain reaction (PCR) methodology excluded this rare variant in an additional 163 unrelated breast cancer patients stratified according to ER status (49 ER-negative, 114 ER-positive). This missense mutation was however detected in one (aged 36 years without a family history of cancer) of the 160 control individuals screened, but absent in the cousin of the index case diagnosed with lobular invasive carcinoma at the age of 42 years.
Failure to identify high impact pathogenic mutations in all affected family members prompted extended analysis for functional polymorphisms (minor allele frequency > 1%) in the folate pathway genes, methionine synthase (MTR) and methionine synthase reductase (MTRR). These genes are included in the Ion AmpliSeq™ Comprehensive Cancer Panel used to identify the RAD50 and MUC1 rare variants. Two functional SNPs were detected in the MTR (2756 A>G, rs1805087) and MTRR (66 A>G, rs1801394) genes, in addition to the low-penetrance MTHFR 677 C>T (rs1801133) mutation previously identified using the testing algorithm illustrated in Figure 1. MTHFR 1298 A>C (rs1801131) known to be in linkage disequilibrium with MTHFR 677 C>T was not detected. Extended analysis of the CYP2D6 gene in the index case due to premature termination of tamoxifen, confirmed a fully functional enzyme based on the gene region covered by WES. Figure 4 illustrates the clinical, histopathological and genetic heterogeneity in the index patient, her daughter and cousin.

3. Discussion

After excluding high-penetrance mutations in the coding regions of the BRCA1 and BRCA2 genes, as well as screening for relatively common low-penetrance SNPs in key disease pathways underlying genomic instability (MTHFR 677 C>T and 1298 A>C) and drug metabolism (CYP2D6*4), WES was performed in the index patient and her daughter. An important aim was to determine whether homozygosity for the MTHFR 677 C>T mutation may be sufficient to cause ER-positive breast cancer in a high-risk environment as reflected by high BMI in both patients. WES identified a rare missense mutation R385C in the DNA mismatch repair gene RAD50, in both the mother and daughter. This rare variant was excluded in a cousin diagnosed with invasive lobular carcinoma at the age of 42 years, as well as a further 163 unrelated breast cancer patients selected from the genomics database. The index case was diagnosed with bilateral breast cancer, invasive ductal and lobular carcinoma, while the daughter was diagnosed with invasive ductal carcinoma. It is therefore possible that RAD50 R385C may only be associated with invasive ductal carcinoma. However, detection of RAD50 R385C in one of 160 controls screened and conflicting results obtained with use of several bioinformatics tools to predict functionality, classifies this missense mutation as a VUS.
Aloraifi et al. [37] identified RAD50 mutations in approximately 3% of genetically uncharacterised breast cancer in a cohort of non-BRCA cases analysed. These authors reported potentially pathogenic variants that may explain 35% of the missing heritability. Bodmer and Tomlinson [38] suggested that whether or not rare variants end up filling the heritability gap, it is imperative to search for rare variants alongside common variants. RAD50 has been documented as an intermediate-risk gene, attributable to its critical role in the repair of DNA double strand breaks [39]. The RAD50 variant (rs139372231) is a C to T base change at nucleotide position 1153 that results in an amino acid change from arginine to cysteine at codon 385. The RAD50 protein forms a complex with the meiotic recombination 11 (MRE11) and nibrin (NBN) double-strand break repair proteins, which bind to strands of damaged DNA and hold broken ends together during the repair process. This MRE11-RAD50-NBN complex is important to prevent accumulation of DNA damage that may trigger uncontrolled cellular division, in order to promote maintenance of genome stability [40]. Pathogenic mutations in the RAD50 gene may lead to the production of an abnormally small, nonfunctional or missing RAD50 protein that may fail to respond effectively to DNA damage [41].
This study has demonstrated the feasibility of WES to achieve a high level of comprehensive genomic testing. The rare missense mutation Q85E, identified in the MUC1 gene in the daughter of the index case with early-onset breast cancer, was absent in the mother and her affected cousin. MUC1 is a transmembrane heterodimeric oncoprotein that is highly expressed in various cancers and correlates with potential for malignancy [42]. Yamada et al. [41] reported that MUC1 gene expression is regulated by DNA methylation and histone modification of the MUC1 promoter. They postulated that an understanding of the epigenetic changes of MUC1 may be of importance for prediction of risk and outcome for patients diagnosed with breast cancer. Deleterious germline mutations in the MUC1 gene is a well-documented cause of medullary cystic kidney disease, also known as mucin-1 kidney disease [43]. In breast cancer, however, mutations in the MUC1 gene may only present as somatic changes within tumour DNA, and are unlikely to be clinically useful in prediction of familial risk. Thus, detection of this germline MUC1 mutation in the daughter may be considered an incidental finding in relation to kidney disease, if proven to be pathogenic in future studies.
A shift in our focus to low-penetrance gene variants highlighted the potential significance of folate pathway SNPs as a risk determinant for hormone receptor-positive breast cancer. Previous studies have demonstrated that the genotype distribution and allele frequencies of SNPs in the MTHFR gene differ significantly between population groups in South Africa [44]. The relatively high frequency of MTHFR 677 C>T is not only of concern for cancer risk, but also in development of other chronic diseases such as CVD. The SNP modifies the effect of hormone replacement therapy (HRT) on lipid and metabolic parameters in postmenopausal women [45]. Homozygosity (TT genotype) for MTHFR 677 C>T is associated with a statistically significant increase in both total and low-density lipoprotein (LDL) cholesterol levels, while presence of the C-allele is associated with a reduction in these lipid levels after one year of HRT. The favourable effect of HRT on lipid levels was reversed in MTHFR TT homozygotes.
The correlation of a given genomic risk profile with clinical and therapeutic outcomes as well as risk for drug toxicity is complicated by clinico-pathological heterogeneity as observed in the family studied here. The low-penetrance MTHFR 677 C>T mutation has previously been linked to the luminal subtype of breast cancer [15], in a similar way that BRCA1 mutations predominate in patients with hormone receptor-negative breast cancer of the basal-type. MTHFR 677 C>T has also been implicated as a causative genetic determinant in altered expression of microRNAs (miRNAs) [46]. These short ≈22 nucleotide RNA regulators are involved in a wide spectrum of cellular processes involving epigenetic changes [47]. The overall importance of polymorphic variation remains high due to multiplicative effects, such that individuals possessing several functional polymorphisms may have a significantly increased risk of breast cancer [48,49,50]. The presence of RAD50 R385C, in combination with MTHFR 677 C>T previously identified as a BRCA-modifier gene [51], is consistent with the modest family history reported by the index case.
Although functional SNPs do not confer significant individual risks for breast cancer, cumulative data supported by WES suggest that co-inheritance of multiple low-penetrance mutations underlying cancer pathways and rare variants in driver genes may explain the majority of non-BRCA familial disease [52,53]. Variation in folate-pathway genes such as MTHFR, MTR and MTRR may disrupt homeostasis and confer significant risk of breast cancer, attributable to impaired DNA repair [54,55]. Sangrajrang et al. [56] highlighted the role of genetic variation in the MTR and MTRR genes in relation to menopausal status of breast cancer patients. Co-inheritance of MTHFR 677C>T with high-penetrance causative mutations in the BRCA1/2 genes may increase the risk of early-onset breast/ovarian cancer, while mutation 1298 A>C in the same gene is associated with sporadic breast cancer [51]. Elucidation of the role of MTHFR as a modifier gene in familial breast cancer is of particular relevance in South Africa, given the high frequency of founder BRCA1/2 mutations in Afrikaner, Coloured and Xhosa breast cancer patients [29,30].
Current literature does not suggest a particular strategy with regards to surveillance and prophylactic management in patients with low-moderate risk variants [16]. However, combining pathway-based SNP genotyping with BRCA founder mutation screening where appropriate, may provide a viable risk reduction strategy for non-communicable diseases characterized by shared disease pathways [4,7]. Due to constantly evolving genetic knowledge, the ideal breast cancer susceptibility test remains elusive as it requires continuous development and refinement as new information on VUS classification and clinical relevance of functional polymorphisms become available in the scientific literature. The identification of actionable functional polymorphisms in important disease pathways covered to a certain extent in the CVD multi-gene assay, which was recently extended to include CYP2D6 pharmacogenetics [28,31], caused a paradigm shift in our understanding of the pieces of the susceptibility puzzle. Implications for genetic counseling remains uncertain and can only be comprehended by incorporating fixed genetic and modifiable environmental risk factors into an existing body of knowledge. The PSGT strategy, which incorporates lifestyle factors with pathology and genetic test results, strives to facilitate patient selection and clinical interpretation [4]. Improved understanding of the genetic overlap between risk factors for ER-positive breast cancer and other chronic diseases of lifestyle warrants consideration of a polygenic approach to disease prevention.
Limitations of the study include the fact that WES evaluates only the subset of DNA regions that encode proteins (known as exons, or collectively the exome), which constitutes about 1% of the human genome. There are approximately 20,000 genes encoded by an estimated 180,000 exons, but due to the technological limitations of NGS previously reviewed [7], it is not feasible to completely cover every single gene with 100% certainty. Clinically relevant mutations could therefore be missed in cases where the gene region containing a disease-causing variant is not represented in the WES reads. Complex genetic insertions or deletions are also not reliably detected by WES. We have previously illustrated these problems in relation to the CYP2D6 gene, which could be addressed by using a combination of different mutation detection methods to obtain reliable results [4]. A similar approach was followed in this study. Furthermore, it has been shown that every healthy person carries several deleterious mutations, without clinical manifestation of the associated medical condition due to gene–gene and gene–environment interaction [57]. For this reason, information on VUS confirmed by Sanger sequencing is stored in a database to gain knowledge from a larger number of samples in future.
The strength of this study lies in capturing genetic, pathology and lifestyle information before and after WES to formulate pathway-based algorithms for breast cancer treatment and prevention. The important role of gene–gene and gene–environment interaction has been emphasised previously to explain the missing heritability in breast cancer. However, to our knowledge, this is the first study to demonstrate the importance of functional SNPs in the MTHFR gene that could be further explored by WES to confirm or support clinical relevance.

4. Materials and Methods

4.1. Ethics Approval

In response to a needs assessment survey amongst healthcare professionals in South Africa, the PSGT platform was established to combine service delivery with health outcome studies [58]. Approval for this study and use of the genomics database was obtained from the Human Research Ethics Committee of Stellenbosch University (reference number N09/08/224).

4.2. Study Population

Three members of the pedigree of the index case were included in this study. The index patient diagnosed with breast cancer at the age of 57 years was referred for the MammaPrint test performed in 2008. She forms part of the cases included in the study of Grant et al. [10], who made special mention of this patient due to a low-risk MammaPrint result for both tumours (invasive ductal and lobular carcinoma). Her daughter was diagnosed with invasive ductal carcinoma approximately one year later, at the age of 29 years. A cousin was diagnosed with invasive lobular carcinoma at the age of 42 years. The index case experienced severe drug side-effects with use of tamoxifen, including vaginal atrophy and hemorrhagic cystitis. She also reported folate supplementation to alleviate side-effects during methotrexate treatment for arthritis.
A total of 164 unrelated female breast cancer patients (Table 2) and 160 controls without a family history of breast cancer were included for extended mutation screening of variants detected by WES. Except for a high BMI, the index case did not report any of the lifestyle risk factors listed in Table 2. Periodic high alcohol consumption was reported in the affected daughter with early-onset breast cancer treated with both tamoxifen and chemotherapy.
Patients with benign neoplasms of the breast or other malignancies were excluded from the study. Only unrelated breast cancer patients with known ER status who signed informed consent for genetic studies were selected for high-throughput genotyping, comprising 60 (37%) Coloured patients of mixed ancestry and 104 (63%) Caucasians, including the index case. Patients with histopathology confirmed breast cancer were recruited at the Tygerberg Hospital Breast Cancer Clinic or referred from private practicing clinicians, as previously described by van der Merwe et al. [28]. Patients referred by clinicians in private practice participated in a chronic disease screen using the CVD multi-gene assay [59] as the core genetic component, recently extended to include CYP2D6 pharmacogenetics [31]. An online questionnaire was used to obtain information on personal and family medical conditions, pathology, medication use/side effects and lifestyle factors to enable clinical interpretation of the genetic findings. Both the index patient and her affected 29-year old daughter tested positive for MTHFR 677 C>T, investigated in this study as a potential contributing factor for breast cancer development against a possible genetic background of rare variants that may be uncovered by WES. In this family with inherited breast cancer, both the daughter (BMI 27.5 kg/m2) and her mother (BMI 27.8 kg/m2) were overweight with suboptimal intake of folate-rich food in the diet (folate score < 14) according to the nutrition and lifestyle assessment performed as previously described by Delport et al. [27]. The presence of BRCA1 and BRCA2 mutations was largely excluded as the cause of familial breast cancer using NGS and multiplex ligation-dependent probe amplification (MLPA) at a commercial laboratory (NewGene, Newcastle, UK).

4.3. DNA Extraction

DNA was extracted from whole blood using the spin protocol of QIAGEN QIAamp® DNA Blood Mini Kit (Hilden, Germany) or Oragene reagents (Ottawa, ON, Canada) for saliva specimens, according to the instructions provided with these commercially available kits. The NanoDrop® ND-1000 Spectrophotometer (Nanodrop Technologies, Wilmington, NC, USA) and its v3.5.2 software package was used to measure the DNA quality and quantity. All genomic DNA samples were diluted to a final concentration of 10 ng/µL using nuclease-free water.

4.4. Whole Exome Sequencing

WES was performed at the Central Analytical Facility (Stellenbosch University, Stellenbosch, South Africa) using the Ion Proton™ System (Thermo Fisher Scientific, Waltham, MA, USA) for NGS. The Ion AmpliSeq™ Exome RDY Library Preparation protocol was used (Publication Number MAN0010084) and template amplification was performed using the Ion PI™ Template OT2 200 Kit v3 (Thermo Fisher Scientific, Waltham, MA, USA) (Publication Number MAN0009133). Semi-conductor sequencing on the Ion Proton system was performed using the Ion PI™ Sequencing 200 Kit v3 (Publication Number MAN0009136) with the Ion PI™ Chip Kit v2. This method is designed to target all human exons. The sequencing run time on the Ion Proton™ System was approximately 4 h which enabled the entire workflow from library construction to primary data analysis and result generation to be done in under 24 h. The Torrent Suite™ (Thermo Fisher Scientific, Waltham, MA, USA) proprietary software is the standard program through which this NGS apparatus operates and offered simple run setup using predesigned workflows and run monitoring in real-time. WES analysis comprised several intermediate steps involved in the read-processing pipeline, as illustrated in Figure 5.
First, sequence reads were aligned to the human reference genome (hg19) using the torrent mapping alignment (TMAP) program [60]. This was followed by read analysis and processing using the preinstalled Torrent Suite™ Software and variant calling with the Torrent Variant Caller (TVC) plug-in. Mapped reads were processed in binary alignment/map (BAM) and indexed sequence alignment/map (SAM) file format. TVC was used for functional annotation and filtering of variant output in variant call files (VCFs). The finding that hg19 contains minor alleles at >1.5 million loci led to repeat mapping of the reads to a synthetic Caucasian major allele reference sequence (CEU-MARS) [61], following the same process illustrated in Figure 5. Resulting VCFs were then processed using GeneTalk [62], a web-based tool for filtering and annotation of uploadable VCFs.
A coverage depth of 15× was used for detection of potentially causative gene variants, prior to confirmation by Sanger sequencing. Coverage provides counts of read depth (DP) at two different levels; the sample level where the DP value is the count of reads that passed the caller’s internal quality control metrics and at the site level where this value is the unfiltered depth over all samples. In addition to the aforementioned quality filters applied, rare variants were filtered on a population frequency of <0.01% with variant function set to exclude all synonymous variants. A list containing 409 genes from the Ion AmpliSeq™ Comprehensive Cancer Panel (Thermo Fisher Scientific, Waltham, MA, USA; catalog number 4477685) was used to search for high impact mutations as well as polymorphic variation in folate-related genes included in the panel. Bioinformatics tools freely available on the internet (Table 1) were used to predict the functional effects for variants of interest.

4.5. Sanger Sequencing and Real-Time Polymerase Chain Reaction (PCR)

Following the identification of potentially causative variants or functional SNPs involved in the folate pathway with use of WES, confirmation was performed by Sanger sequencing. Selected variants were screened for in controls (n = 160) and the total breast cancer sample (n = 164) using real-time PCR. Two microliters of template DNA (10 ng/µL) was used in a total reaction volume of 10 µL. Applied Biosystems® (ABI™) (Foster City, CA, USA) TaqMan® SNP Genotyping Assays were used for the real-time PCR runs according to standard TaqMan® SNP Genotyping protocols. Pre-designed 40× assay mix consisted of unlabelled primers and TaqMan® Minor Groove Binder (MGB) probes (VIC® and FAM™ dye-labelled). These assays were used for end-point genotyping by allelic discrimination analysis for SNPs on the Corbett Rotor-Gene™ 6000 (Mortlake, New South Wales, Australia) or Roche LightCycler® 480 II (Basel, Switzerland). The oligonucleotide primer sequences used for Sanger sequencing and the TaqMan assays employed in this study are presented in Table 3.

5. Conclusions

This study describes a framework for WES performed alongside clinical and pathology assessments to provide a clear strategy for the way forward. We moved from single-gene analysis of BRCA1/2 founder mutations in the high-risk South African population to multi-gene panel testing, meant to assist clinicians in the decision to pursue WES. The selection of the index patient was primarily based on the need to establish the role of MTHFR 677 C>T in familial risk, due to the potential for cancer prevention at the population level. Our results support previous findings indicating that the majority of genetically uncharacterised breast cancer cases may be caused by a combination of low–moderate penetrance mutations exerting their effect against a high-risk clinical background. To our knowledge, this is the first study using WES to investigate the significance of folate pathway SNPs as risk reduction targets beyond BRCA1/2 in familial risk. WES preceded by PSGT facilitated the identification and clinical interpretation of genetic risk factors of relevance to both cancer development and tailored therapeutic intervention in a single test.

Acknowledgments

Research reported in this publication was supported by the Strategic Health Innovation Partnerships (SHIP) Unit of the South African Medical Research Council (MRC) with funds received from the South African Department of Science and Technology (Research grant number S003665). This work is also based on the research supported by the Cancer Association of South Africa and the National Research Foundation (NRF) of South Africa (UID 83962). The grant holders acknowledge that opinions, findings and conclusions or recommendations expressed in this publication are those of the authors. Carel van Heerden of the Central Analytical Facility, supported by NRF grant UID 88063, is thanked for exome sequencing and critical reading of the manuscript.

Author Contributions

Nicole van der Merwe wrote the first draft of the manuscript from data included in her doctoral thesis, participated in the sample collection and development of the genomics database, performed high-throughput genotyping and bioinformatics. Armand V. Peeters verified the bioinformatic analysis and performed variant confirmation by real-time PCR and Sanger sequencing. Fredrieka M. Pienaar referred the index patient for genetic testing and provided clinical information of the extended patient sample from medical records. Juanita Bezuidenhout provided oversight of histopathology assessments. Susan J. van Rensburg participated in the study design and successful completion of the study. Maritha J. Kotze conceived the study, enabled development of the database used for clinical interpretation and edited the manuscript into a logical flow of information. All authors read and approved the final manuscript.

Conflicts of Interest

Maritha J. Kotze is a director and shareholder of Gknowmix (Pty) Ltd. that has developed a database tool for research translation under the auspices of the South African Medical Research Council. Susan J. van Rensburg is a scientific advisor of Gknowmix, with no potential conflict of interest reported by the other authors. Maritha J. Kotze is the inventor of a South African Medical Research Council patent 2001/5419 (2 July 2001) entitled “A method of diagnosing patients with cardiovascular disease or a genetic predisposition for cardiovascular disease”. The funding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results.

Abbreviations

ABIApplied Biosystems
BAMBinary alignment/map
BMIBody mass index
CEU-MARSCaucasian major allele reference sequence
CVDCardiovascular disease
CYPCytochrome P450
dnSNPDatabase of Single Nucleotide Polymorphisms
EREstrogen receptor
hgHuman reference genome
HRTHormone replacement therapy
IGVIntegrative Genome Viewer
MINDACTMicroarray In Node negative Disease may Avoid Chemo Therapy
MRCMedical Research Council
MRE11Meiotic recombination 11
MTHFRMethylene tetrahydrofolate reductase
MTRMethionine synthase
MTRRMethionine synthase reductase
MUC1Mucin 1
NBNNibrin
NGSNext generation sequencing
NRFNational Research Foundation
PCRPolymerase chain reaction
PRProgesterone receptor
PSGTPathology-supported genetic testing
RAD50Double strand break repair protein
SAMSequence alignment/map
SHIPStrategic Health Innovation Partnerships
SNPsSingle nucleotide polymorphisms
TMAPTorrent mapping alignment
TVCTorrent Variant Caller
VCFVariant call files
VUSVariant of uncertain clinical significance
WESWhole exome sequencing

References

  1. Nordgard, S.H.; Alnaes, G.I.; Hihn, B.; Lingjaerde, O.C.; Liestøl, K.; Tsalenko, A.; Sørlie, T.; Lønning, P.E.; Børresen-Dale, A.L.; Kristensen, V.N. Pathway based analysis of SNPs with relevance to 5-FU therapy: Relation to intratumoral mRNA expression and survival. Int. J. Cancer 2008, 123, 577–585. [Google Scholar] [CrossRef] [PubMed]
  2. Babyshkina, N.; Malinovskaya, E.; Nazarenko, M.; Koval, M.; Gervas, P.; Potapova, O.; Slonimskaya, E.; Cherdyntseva, N. The effect of folate-related SNPs on clinicopathological features, response to neoadjuvant treatment and survival in pre- and postmenopausal breast cancer patients. Gene 2013, 518, 397–404. [Google Scholar] [CrossRef] [PubMed]
  3. Seymour, H.J.; Wainstein, T.; Macaulay, S.; Haw, T.; Krause, A. Breast cancer in high-risk Afrikaner families: Is BRCA founder mutation testing sufficient? S. Afr. Med. J. 2016, 106, 264–267. [Google Scholar] [CrossRef] [PubMed]
  4. Kotze, M.J. Application of advanced molecular technology in the diagnosis and management of genetic disorders in South Africa. S. Afr. Med. J. 2016, 106, S114–S118. [Google Scholar] [CrossRef] [PubMed]
  5. Bahassi, E.M.; Stambrook, P.J. Next-generation sequencing technologies: Breaking the sound barrier of human genetics. Mutagenesis 2014, 29, 303–310. [Google Scholar] [CrossRef] [PubMed]
  6. Masser, D.R.; Stanford, D.R.; Freeman, W.M. Targeted DNA methylation analysis by next-generation sequencing. J. Vis. Exp. 2015, 96, 52488. [Google Scholar] [CrossRef] [PubMed]
  7. Kotze, M.J.; Lückhoff, H.K.; Peeters, A.V.; Baatjes, K.; Schoeman, M.; van der Merwe, L.; Grant, K.A.; Fisher, L.R.; van der Merwe, N.; Pretorius, J.; et al. Genomic medicine and risk prediction across the disease spectrum. Crit. Rev. Clin. Lab. Sci. 2015, 52, 120–137. [Google Scholar] [CrossRef] [PubMed]
  8. Van der Groep, P.; van der Wall, E.; van Diest, P.J. Pathology of hereditary breast cancer. Cell. Oncol. 2011, 34, 71–88. [Google Scholar] [CrossRef] [PubMed]
  9. Cheon, J.Y.; Mozersky, J.; Cook-Deegan, R. Variants of uncertain significance in BRCA: A harbinger of ethical and policy issues to come? Genome Med. 2014, 6, 121–131. [Google Scholar] [CrossRef] [PubMed]
  10. Grant, K.A.; Apffelstaed, J.P.; Wright, C.; Myburgh, E.; Pienaar, R.; De Klerk, M.; Kotze, M.J. MammaPrint pre-screen algorithm (MPA) reduces chemotherapy in patients with early-stage breast cancer. S. Afr. Med. J. 2013, 103, 522–526. [Google Scholar] [CrossRef] [PubMed]
  11. Cardoso, F.; van’t Veer, L.J.; Bogaerts, J.; Slaets, L.; Viale, G.; Delaloge, S.; Pierga, J.Y.; Brain, E.; Causeret, S.; DeLorenzi, M.; et al. 70-gene signature as an aid to treatment decisions in early-stage breast cancer. N. Engl. J. Med. 2016, 375, 717–729. [Google Scholar] [CrossRef] [PubMed]
  12. Larsen, M.J.; Thomassen, M.; Tan, Q.; Lænkholm, A.V.; Bak, M.; Sørensen, K.P.; Andersen, M.K.; Kruse, T.A.; Gerdes, A.M. RNA profiling reveals familial aggregation of molecular subtypes in non-BRCA1/2 breast cancer families. BMC Med. Genom. 2014, 7, 9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Yao, K.; Goldschmidt, R.; Turk, M.; Wesseling, J.; Stork-Sloots, L.; de Snoo, F.; Cristofanilli, M. Molecular subtyping improves diagnostic stratification of patients with primary breast cancer into prognostically defined risk groups. Breast Cancer Res. Treat. 2015, 154, 81–88. [Google Scholar] [CrossRef] [PubMed]
  14. Myburgh, E.J.; Langenhoven, L.; Grant, K.A.; van der Merwe, L.; Kotze, M.J. Clinical overestimation of HER2 positivity in early estrogen and progesterone receptor–positive breast cancer and the value of molecular subtyping using BluePrint. J. Glob. Oncol. 2016. [Google Scholar] [CrossRef]
  15. Naushad, S.M.; Pavani, A.; Rupasree, Y.; Divyya, S.; Deepti, S.; Digumarti, R.R.; Gottumukkala, S.R.; Prayaga, A.; Kutala, V.K. Association of aberrations in one-carbon metabolism with molecular phenotype and grade of breast cancer. Mol. Carcinog. 2012, 51, E32–E41. [Google Scholar] [CrossRef] [PubMed]
  16. Ghoussaini, M.; Pharoah, P.D.P.; Easton, D.F. Inherited genetic susceptibility to breast cancer—The beginning of the end or the end of the beginning? Am. J. Pathol. 2013, 183, 1038–1051. [Google Scholar] [CrossRef] [PubMed]
  17. Easton, D.F.; Pooley, K.A.; Dunning, A.M.; Pharoah, P.D.; Thompson, D.; Ballinger, D.G.; Struewing, J.P.; Morrison, J.; Field, H.; Luben, R.; et al. Genome-wide association study identifies novel breast cancer susceptibility loci. Nature 2007, 447, 1087–1093. [Google Scholar] [CrossRef] [PubMed]
  18. Easton, D.F.; Pharoah, P.D.; Antoniou, A.C.; Tischkowitz, M.; Tavtigian, S.V.; Nathanson, K.L.; Devilee, P.; Meindl, A.; Couch, F.J.; Southey, M.; et al. Gene-panel sequencing and the prediction of breast-cancer risk. N. Engl. J. Med. 2015, 372, 2243–2257. [Google Scholar] [CrossRef] [PubMed]
  19. Bordeleau, L.; Panchal, S.; Goodwin, P. Prognosis of BRCA-associated breast cancer: A summary of evidence. Breast Cancer Res. Treat. 2010, 119, 13–24. [Google Scholar] [CrossRef] [PubMed]
  20. Bijl, M.J.; van Schaik, R.H.; Lammers, L.A.; Hofman, A.; Vulto, A.G.; van Gelder, T.; Stricker, B.H.; Visser, L.E. The CYP2D6*4 polymorphism affects breast cancer survival in tamoxifen users. Breast Cancer Res. Treat. 2009, 18, 125–130. [Google Scholar] [CrossRef] [PubMed]
  21. Kimura, M.; Umegaki, K.; Higuchi, M.; Thomas, P.; Fenech, M. Methylenetetrahydrofolate reductase C677T polymorphism, folic acid and riboflavin are important determinants of genome stability in cultured human lymphocytes. J. Nutr. 2004, 134, 48–56. [Google Scholar] [PubMed]
  22. Xu, X.; Gammon, M.D.; Wetmur, J.G.; Bradshaw, P.T.; Teitelbaum, S.L.; Neugut, A.I.; Santella, R.M.; Chen, J. B-vitamin intake, one-carbon metabolism, and survival in a population-based study of women with breast cancer. Cancer Epidemiol. Biomark. Prev. 2008, 17, 2109–2116. [Google Scholar] [CrossRef] [PubMed]
  23. Paz, M.F.; Avila, S.; Fraga, M.F.; Pollan, M.; Capella, G.; Peinado, M.A.; Sanchez-Cespedes, M.; Herman, J.G.; Esteller, M. Germ-line variants in methyl-group metabolism genes and susceptibility to DNA methylation in normal tissues and human primary tumours. Cancer Res. 2002, 62, 4519–4524. [Google Scholar] [PubMed]
  24. Wishart, D.S. Is cancer a genetic disease or a metabolic disease? EBioMedicine 2015, 2, 478–479. [Google Scholar] [CrossRef] [PubMed]
  25. Bailey, L.B. Folate, methyl-related nutrients, alcohol, and the MTHFR 677 C>T polymorphism affect cancer risk: Intake recommendations. J. Nutr. 2003, 133, 3748S–3753S. [Google Scholar] [PubMed]
  26. McEligor, A.J.; Ziogas, A.; Pfeiffer, C.M.; Fazil, Z.; Anton-Culver, H. The association between circulating total folate and folate vitamers with overall survival after postmenopausal breast cancer diagnosis. Nutr. Cancer 2015, 67, 442–448. [Google Scholar] [CrossRef] [PubMed]
  27. Delport, D.; Schoeman, R.; van der Merwe, N.; van der Merwe, L.; Fisher, L.R.; Geiger, D.; Kotze, M.J. Significance of dietary folate intake, homocysteine levels and MTHFR 677 C>T genotyping in South African patients diagnosed with depression: Test development for clinical application. Metab. Brain Dis. 2014, 29, 377–384. [Google Scholar] [CrossRef] [PubMed]
  28. Van der Merwe, N.; Bouwens, C.S.; Pienaar, R.; van der Merwe, L.; Yako, Y.Y.; Geiger, D.H.; Kotze, M.J. CYP2D6 genotyping and use of antidepressants in breast cancer patients: Test development for clinical application. Metab. Brain Dis. 2012, 27, 319–326. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Reeves, M.D.; Yawitch, T.M.; van der Merwe, N.C.; van den Berg, H.J.; Dreyer, G.; van Rensburg, E.J. BRCA1 mutations in South African breast and/or ovarian cancer families: Evidence of a novel founder mutation in Afrikaner families. Int. J. Cancer 2004, 110, 677–682. [Google Scholar] [CrossRef] [PubMed]
  30. Van der Merwe, N.C.; Hamel, N.; Schneider, S.R.; Apffelstaedt, J.P.; Wijnen, J.T.; Foulkes, W.D. A founder BRCA2 mutation in non-Afrikaner breast cancer patients of the Western Cape of South Africa. Clin. Genet. 2012, 81, 179–184. [Google Scholar] [CrossRef] [PubMed]
  31. Van der Merwe, N. Development and Application of a Pathology Supported Pharmacogenetic Test for Improved Clinical Management of South African Patients with Breast Cancer and Associated Comorbidities. Ph.D. Thesis, Stellenbosch University, Stellenbosch, South Africa, 2016. [Google Scholar]
  32. Phipps, A.I.; Buist, D.S.; Malone, K.E.; Barlow, W.E.; Porter, P.L.; Kerlikowske, K.; Li, C.I. Family history of breast cancer in first-degree relatives and triple-negative breast cancer risk. Breast Cancer Res. Treat. 2011, 126, 671–678. [Google Scholar] [CrossRef] [PubMed]
  33. Phipps, A.I.; Buist, D.S.; Malone, K.E.; Barlow, W.E.; Porter, P.L.; Kerlikowske, K.; O’Meara, E.S.; Li, C.I. Breast density, body mass index, and risk of tumour marker-defined subtypes of breast cancer. Ann. Epidemiol. 2012, 22, 340–348. [Google Scholar] [CrossRef] [PubMed]
  34. Bardia, A.; Arieas, E.T.; Zhang, Z.; Defilippis, A.; Tarpinian, K.; Jeter, S.; Nguyen, A.; Henry, N.L.; Flockhart, D.A.; Hayes, D.F.; et al. Comparison of breast cancer recurrence risk and cardiovascular disease incidence risk among postmenopausal women with breast cancer. Breast Cancer Res. Treat. 2012, 131, 907–914. [Google Scholar] [CrossRef] [PubMed]
  35. Li, C.I.; Daling, J.R.; Porter, P.L.; Tang, M.T.; Malone, K.E. Relationship between potentially modifiable lifestyle factors and risk of second primary contralateral breast cancer among women diagnosed with estrogen receptor-positive invasive breast cancer. J. Clin. Oncol. 2009, 27, 5312–5318. [Google Scholar] [CrossRef] [PubMed]
  36. Integrative Genomics Viewer. Available online: http://software.broadinstitute.org/software/igv/ (accessed on 31 July 2015).
  37. Aloraifi, F.; Boland, M.R.; Green, A.J.; Geraghty, J.G. Gene analysis techniques and susceptibility gene discovery in non-BRCA1/BRCA2 familial breast cancer. Surg. Oncol. 2015, 24, 100–109. [Google Scholar] [CrossRef] [PubMed]
  38. Bodmer, W.; Tomlinson, I. Rare genetic variants and the risk of cancer. Curr. Opin. Genet. Dev. 2010, 20, 262–267. [Google Scholar] [CrossRef] [PubMed]
  39. Damiola, F.; Pertesi, M.; Oliver, J.; Le Calvez-Kelm, F.; Voegele, C.; Young, E.L.; Robinot, N.; Forey, N.; Durand, G.; Vallée, M.P.; et al. Rare key functional domain missense substitutions in MRE11A, RAD50, and NBN contribute to breast cancer susceptibility: Results from a Breast Cancer Family Registry case-control mutation-screening study. Breast Cancer Res. 2014, 16, R58. [Google Scholar] [CrossRef] [PubMed]
  40. Heikkinen, K.; Rapakko, K.; Karppinen, S.M.; Erkko, H.; Knuutila, S.; Lundán, T.; Mannermaa, A.; Børresen-Dale, A.L.; Borg, A.; Barkardottir, R.B.; et al. RAD50 and NBS1 are breast cancer susceptibility genes associated with genomic instability. Carcinogenesis 2006, 27, 1593–1599. [Google Scholar] [CrossRef] [PubMed]
  41. Bartkova, J.; Tommiska, J.; Oplustilova, L.; Aaltonen, K.; Tamminen, A.; Heikkinen, T.; Mistrik, M.; Aittomäki, K.; Blomqvist, C.; Heikkilä, P. Aberrations of the MRE11–RAD50–NBS1 DNA damage sensor complex in human breast cancer: MRE11 as a candidate familial cancer predisposing gene. Mol. Oncol. 2008, 2, 296–316. [Google Scholar] [CrossRef] [PubMed]
  42. Yamada, N.; Nishida, Y.; Tsutsumida, H.; Hamada, T.; Goto, M.; Higashi, M.; Nomoto, M.; Yonezawa, S. MUC1 expression is regulated by DNA methylation and histone H3 lysine 9 modification in cancer cells. Cancer Res. 2008, 68, 2708–2716. [Google Scholar] [CrossRef] [PubMed]
  43. Kirby, A.; Gnirke, A.; Jaffe, D.B.; Barešová, V.; Pochet, N.; Blumenstiel, B.; Ye, C.; Aird, D.; Stevens, C.; Robinson, J.T. Mutations causing medullary cystic kidney disease type 1 lie in a large VNTR in MUC1 missed by massively parallel sequencing. Nat. Genet. 2013, 45, 299–303. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Scholtz, C.L.; Odendaal, H.J.; Thiart, R.; Loubser, L.; Hillermann, R.; Delport, R.; Vermaak, W.J.; Kotze, M.J. Analysis of two mutations in the MTHFR gene associated with mild hyperhomocysteinaemia: Heterogeneous distribution in the South African population. S. Afr. Med. J. 2002, 92, 464–467. [Google Scholar] [PubMed]
  45. Lambrinoudaki, I.; Papadimitriou, D.; Kaparos, G.; Rizos, D.; Panoulis, C.; Deligeoroglou, E.; Alexandrou, A.; Auguolea, A.; Apostolakis, M.; Creatsa, M.; et al. MTHFR C677 T polymorphism modifies the effect of HRT on metabolic parameters in postmenopausal women. Climacteric 2013, 16, 568–575. [Google Scholar] [CrossRef] [PubMed]
  46. Stone, N.; Pangilinan, F.; Molloy, A.M.; Shane, B.; Scott, J.M.; Ueland, P.M.; Mills, J.L.; Kirke, P.N.; Sethupathy, P.; Brody, L.C. Bioinformatic and genetic association analysis of microRNA target sites in one-carbon metabolism genes. PLoS ONE 2011, 6, e21851. [Google Scholar] [CrossRef] [PubMed]
  47. Heijmans, B.T.; Kremer, D.; Tobi, E.W.; Boomsma, D.I.; Slagboom, P.E. Heritable rather than age-related environmental and stochastic factors dominate variation in DNA methylation of the human IGF2/H19 locus. Hum. Mol. Genet. 2007, 16, 547–554. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Antoniou, A.C.; Spurdle, A.B.; Sinilnikova, O.M.; Healey, S.; Pooley, K.A.; Schmutzler, R.K.; Versmold, B.; Engel, C.; Meindl, A.; Arnold, N.; et al. Common breast cancer-predisposition alleles are associated with breast cancer risk in BRCA1 and BRCA2 mutation carriers. Am. J. Hum. Genet. 2008, 82, 937–948. [Google Scholar] [CrossRef] [PubMed]
  49. Pharoah, P.D.; Antoniou, A.C.; Easton, D.F.; Ponder, B.A. Polygenes, risk prediction, and targeted prevention of breast cancer. N. Engl. J. Med. 2008, 358, 2796–2803. [Google Scholar] [CrossRef] [PubMed]
  50. Ricks-Santi, L.J.; Nie, J.; Marian, C.; Ochs-Balcom, H.M.; Trevisan, M.; Edge, S.B.; Kanaan, Y.; Freudenheim, J.L.; Shields, P.G. BRCA1 Polymorphisms and breast cancer epidemiology in the Western New York Exposures and Breast Cancer [WEB] Study. Genet. Epidemiol. 2013, 37, 504–511. [Google Scholar] [CrossRef] [PubMed]
  51. Pepe, C.; Guidugli, L.; Sensi, E.; Aretini, P.; D’Andrea, E.; Montagna, M.; Manoukian, S.; Ottini, L.; Radice, P.; Viel, A.; et al. Methyl group metabolism gene polymorphisms as modifiers of breast cancer risk in Italian BRCA1/2 carriers. Breast Cancer Res. Treat. 2007, 103, 29–36. [Google Scholar] [CrossRef] [PubMed]
  52. Antoniou, A.C.; Easton, D.F. Models of genetic susceptibility to breast cancer. Oncogene 2006, 25, 5898–5905. [Google Scholar] [CrossRef] [PubMed]
  53. Gracia-Aznarez, F.J.; Fernandez, V.; Pita, G.; Peterlongo, P.; Dominguez, O.; de la Hoya, M.; Duran, M.; Osorio, A.; Moreno, L.; Gonzalez-Neira, A.; et al. Whole exome sequencing suggests much of non-BRCA1/BRCA2 familial breast cancer is due to moderate and low penetrance susceptibility alleles. PLoS ONE 2003, 8, e55681. [Google Scholar] [CrossRef] [PubMed]
  54. Macis, D.; Maisonneuve, P.; Johansson, H.; Bonanni, B.; Botteri, E.; Iodice, S.; Santillo, B.; Penco, S.; Gucciardo, G.; D’Aiuto, G.; et al. Methylenetetrahydrofolate reductase (MTHFR) and breast cancer risk: A nested-case-control study and a pooled meta-analysis. Breast Cancer Res. Treat. 2007, 106, 263–271. [Google Scholar] [CrossRef] [PubMed]
  55. Weiwei, Z.; Liping, C.; Dequan, L. Association between dietary intake of folate, vitamin B6, B12 & MTHFR, MTR Genotype and breast cancer risk. Pak. J. Med. Sci. 2014, 30, 106–110. [Google Scholar] [PubMed]
  56. Sangrajrang, S.; Sato, Y.; Sakamoto, H.; Ohnami, S.; Khuhaprema, T.; Yoshida, T. Genetic polymorphisms in folate and alcohol metabolism and breast cancer risk: A case-control study in Thai women. Breast Cancer Res. Treat. 2010, 123, 885–893. [Google Scholar] [CrossRef] [PubMed]
  57. Xue, Y.; Chen, Y.; Ayub, Q.; Huang, N.; Ball, E.V.; Mort, M.; Phillips, A.D.; Shaw, K.; Stenson, P.D.; Cooper, D.N.; et al. Deleterious- and disease-allele prevalence in healthy individuals: Insights from current predictions, mutation databases, and population-scale resequencing. Am. J. Hum. Genet. 2012, 91, 1022–1032. [Google Scholar] [CrossRef] [PubMed]
  58. Kotze, M.J.; van Velden, D.P.; Botha, K.; Badenhorst, C.H.; Avenant, H.J.; van Rensburg, S.; Cronje, F. Pathology-supported genetic testing directed at shared disease pathways for optimized health in later life. Personal. Med. 2013, 10, 497–507. [Google Scholar] [CrossRef]
  59. Kotze, M.J.; Thiart, R. Genetics of dyslipidaemia. CME J. 2003, 21, 399–402. [Google Scholar]
  60. Github. Available online: https://github.com/iontorrent/TS/tree/master/Analysis/TMAP (accessed on 7 April 2015).
  61. Dewey, F.E.; Chen, R.; Cordero, S.P.; Ormond, K.E.; Caleshu, C.; Karczewski, K.J.; Whirl-Carrillo, M.; Wheeler, M.T.; Dudley, J.T.; Byrnes, J.K.; et al. Phased whole-genome genetic risk in a family quartet using a major allele reference sequence. PLoS Genet. 2011, 7, e1002280. [Google Scholar] [CrossRef] [PubMed]
  62. GeneTalk. Available online: https://www.gene-talk.de/ (accessed on 12 April 2015).
Ijms 18 00467 g001 550
Figure 1. Outline of the pathology-supported genetic testing (PSGT) strategy used to select the index case with bilateral BRCA1/2 mutation-negative breast cancer for whole exome sequencing (WES). Estrogen receptor (ER), progesterone receptor (PR) and human epidermal growth factor receptor-2 (HER2) status was determined using immunohistochemistry (IHC). These features, including tumour size, grade and nodal status were evaluated to determine eligibility for MammaPrint. Based on the low-risk recurrence score obtained with this microarray test, tamoxifen treatment was administered without addition of chemotherapy. Termination of tamoxifen after one year due to adverse drug response led to pharmacogenetic CYP2D6 genotyping, which was performed as part of a chronic disease screen incorporating the cardiovascular disease (CVD) multi-gene assay. Detection of the methylenetetrahydrofolate reductase (MTHFR) 677 C>T mutation identified dysfunction of the folate pathway as a potential target for recurrence risk reduction in this overweight patient. Evaluation of family history, clinical characteristics, and laboratory test results at the protein (IHC), RNA (MammaPrint) and DNA (BRCA1/2, CVD multi-gene assay) levels supported the decision to perform extended genetic testing using WES. BMI: body mass index; CYP: cytochrome P450.
Figure 1. Outline of the pathology-supported genetic testing (PSGT) strategy used to select the index case with bilateral BRCA1/2 mutation-negative breast cancer for whole exome sequencing (WES). Estrogen receptor (ER), progesterone receptor (PR) and human epidermal growth factor receptor-2 (HER2) status was determined using immunohistochemistry (IHC). These features, including tumour size, grade and nodal status were evaluated to determine eligibility for MammaPrint. Based on the low-risk recurrence score obtained with this microarray test, tamoxifen treatment was administered without addition of chemotherapy. Termination of tamoxifen after one year due to adverse drug response led to pharmacogenetic CYP2D6 genotyping, which was performed as part of a chronic disease screen incorporating the cardiovascular disease (CVD) multi-gene assay. Detection of the methylenetetrahydrofolate reductase (MTHFR) 677 C>T mutation identified dysfunction of the folate pathway as a potential target for recurrence risk reduction in this overweight patient. Evaluation of family history, clinical characteristics, and laboratory test results at the protein (IHC), RNA (MammaPrint) and DNA (BRCA1/2, CVD multi-gene assay) levels supported the decision to perform extended genetic testing using WES. BMI: body mass index; CYP: cytochrome P450.
Ijms 18 00467 g001
Ijms 18 00467 g002 550
Figure 2. Detection of the RAD50 missense mutation R385C (rs139372231, Refseq mRNA accession NM_005732.3) using whole exome sequencing (A); and visualized using the Integrative Genome Viewer software tool [36]. The C to T base change at nucleotide position 1153 was confirmed by Sanger sequencing (B) as indicated by the arrow.
Figure 2. Detection of the RAD50 missense mutation R385C (rs139372231, Refseq mRNA accession NM_005732.3) using whole exome sequencing (A); and visualized using the Integrative Genome Viewer software tool [36]. The C to T base change at nucleotide position 1153 was confirmed by Sanger sequencing (B) as indicated by the arrow.
Ijms 18 00467 g002
Ijms 18 00467 g003 550
Figure 3. Detection of the MUC1 missense mutation Q85E (rs773704188, Refseq mRNA accession NM_002456.5) using WES (A); and visualized using the Integrative Genome Viewer software tool [36]. The C to G base change at nucleotide position 325 was confirmed by Sanger sequencing (B) as indicated by the arrow.
Figure 3. Detection of the MUC1 missense mutation Q85E (rs773704188, Refseq mRNA accession NM_002456.5) using WES (A); and visualized using the Integrative Genome Viewer software tool [36]. The C to G base change at nucleotide position 325 was confirmed by Sanger sequencing (B) as indicated by the arrow.
Ijms 18 00467 g003
Ijms 18 00467 g004 550
Figure 4. Pedigree of the index patient (arrow) subjected to whole exome sequencing, depicting the clinical (age of onset), pathology (tumour characteristics) and genetic heterogeneity in the family. Homozygosity for a minor allele is indicated by (+/+) and heterozygosity by (+/−). ER+: estrogen receptor-positive.
Figure 4. Pedigree of the index patient (arrow) subjected to whole exome sequencing, depicting the clinical (age of onset), pathology (tumour characteristics) and genetic heterogeneity in the family. Homozygosity for a minor allele is indicated by (+/+) and heterozygosity by (+/−). ER+: estrogen receptor-positive.
Ijms 18 00467 g004
Ijms 18 00467 g005 550
Figure 5. A simplified representation of the data analysis pipeline used by Torrent Suite for the identification of variants in WES data. Base calling, read filtering and trimming were performed by the Torrent Suite software. The torrent mapping alignment (TMAP) [60] program was used to map the resulting reads to the human reference genome. Variant calling was performed by the Torrent Variant Caller (TVC) in three steps: generate candidates, evaluate candidates and filter candidates. Final filtering and annotation of variant call files (VCF) were performed using GeneTalk [62] web-based bioinformatics tool.
Figure 5. A simplified representation of the data analysis pipeline used by Torrent Suite for the identification of variants in WES data. Base calling, read filtering and trimming were performed by the Torrent Suite software. The torrent mapping alignment (TMAP) [60] program was used to map the resulting reads to the human reference genome. Variant calling was performed by the Torrent Variant Caller (TVC) in three steps: generate candidates, evaluate candidates and filter candidates. Final filtering and annotation of variant call files (VCF) were performed using GeneTalk [62] web-based bioinformatics tool.
Ijms 18 00467 g005
Table
Table 1. Bioinformatics analysis of missense mutation RAD50 R385C as a variant of uncertain clinical significance.
Table 1. Bioinformatics analysis of missense mutation RAD50 R385C as a variant of uncertain clinical significance.
Bioinformatics Tool AppliedPredicted Biological Effect
Sorting Intolerant From Tolerant (SIFT)Damaging (rank score = 0.7209)
PolyphenBenign (score = 0.339)
Variant Effect PredictorModerate impact
MutationTasterDisease causing (rank score = 0.8103)
MetalRTolerated (rank score = 0.2311)
ProveanNeutral (rank score = 0.4681)
Likelihood Ratio TestDeleterious (rank score = 0.5373)
Table
Table 2. Characteristics of 164 female breast cancer patients subdivided according to tumour pathology into estrogen receptor (ER) positive (n = 49) and negative (n = 115) cases.
Table 2. Characteristics of 164 female breast cancer patients subdivided according to tumour pathology into estrogen receptor (ER) positive (n = 49) and negative (n = 115) cases.
VariablesER-NegativeER-Positive
Median (Range)Number49115 *
Age (median, range)16449 (30–77)54 (31–83)
Body mass index14625 (17–47)26 (17–41)
Count (Percentage)
Ethnicity †16422 (45)38 (33)
Family history of cancer16330 (61)59 (52)
Hormone replacement therapy1639 (19)15 (13)
Oral contraceptives16320 (42)32 (28)
Current smoking16017 (35)34 (30)
High alcohol consumption15625 (52)62 (57)
* Includes the index case; † Percentage of non-Caucasian patients are indicated in brackets.
Table
Table 3. Oligonucleotide primers used for conventional polymerase chain reaction (PCR) and Sanger sequencing, and the identification numbers of commercially available TaqMan assays used for real-time PCR.
Table 3. Oligonucleotide primers used for conventional polymerase chain reaction (PCR) and Sanger sequencing, and the identification numbers of commercially available TaqMan assays used for real-time PCR.
GenedbSNP ID#PrimerOligonucleotide Primers for Sanger SequencingSize (bp)TaqMan Assay ID Numbers
MTHFRrs1801133ForwardATCCCTCGCCTTGAACA256C_1202889_20
ReverseTCACCTGGATGGGAAAGAT
MTRrs1805087ForwardGAACATCCCAAGCCCAC595C_12005959_10
ReverseCACCTGTTTCCCTGCTG
MTRRrs1801394ForwardGTTTCATTCGTACACTCTCC616C_3068176_10
ReverseCAGCATATGCTACTTCTGTC
RAD50rs139372231ForwardATCCACATGCTCAGGGGTAC528C_171053490_10
ReverseGCCAAAATGGAGTCCAACC
MUC1rs773704188ForwardATTCCCAGCCACCACTCTGA493Not available
ReverseCCCAACCTTAAGTGCACCAGT
bp: base pair, dnSNP: database of single nucleotide polymorphisms.

Share and Cite

MDPI and ACS Style

Van der Merwe, N.; Peeters, A.V.; Pienaar, F.M.; Bezuidenhout, J.; Van Rensburg, S.J.; Kotze, M.J. Exome Sequencing in a Family with Luminal-Type Breast Cancer Underpinned by Variation in the Methylation Pathway. Int. J. Mol. Sci. 2017, 18, 467. https://doi.org/10.3390/ijms18020467

AMA Style

Van der Merwe N, Peeters AV, Pienaar FM, Bezuidenhout J, Van Rensburg SJ, Kotze MJ. Exome Sequencing in a Family with Luminal-Type Breast Cancer Underpinned by Variation in the Methylation Pathway. International Journal of Molecular Sciences. 2017; 18(2):467. https://doi.org/10.3390/ijms18020467

Chicago/Turabian Style

Van der Merwe, Nicole, Armand V. Peeters, Fredrieka M. Pienaar, Juanita Bezuidenhout, Susan J. Van Rensburg, and Maritha J. Kotze. 2017. "Exome Sequencing in a Family with Luminal-Type Breast Cancer Underpinned by Variation in the Methylation Pathway" International Journal of Molecular Sciences 18, no. 2: 467. https://doi.org/10.3390/ijms18020467

APA Style

Van der Merwe, N., Peeters, A. V., Pienaar, F. M., Bezuidenhout, J., Van Rensburg, S. J., & Kotze, M. J. (2017). Exome Sequencing in a Family with Luminal-Type Breast Cancer Underpinned by Variation in the Methylation Pathway. International Journal of Molecular Sciences, 18(2), 467. https://doi.org/10.3390/ijms18020467

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop