Characterization of <i>BRCA1/2</i> mutations in patients with family history of breast cancer in Armenia

Sofi Atshemyan; Andranik Chavushyan; Nerses Berberian; Arthur Sahakyan; Roksana Zakharyan; Arsen Arakelyan

doi:10.12688/f1000research.10434.1

Home Browse Characterization of BRCA1/2 mutations in patients with family history...

ALL Metrics

-

Views

-

Downloads

Get PDF

Get XML

Export

▬

✚

Research Article

Characterization of BRCA1/2 mutations in patients with family history of breast cancer in Armenia

[version 1; peer review: 2 approved]

Sofi Atshemyan¹, Andranik Chavushyan¹, Nerses Berberian², Arthur Sahakyan², Roksana Zakharyan¹, Arsen Arakelyan ¹

Sofi Atshemyan¹, Andranik Chavushyan¹, [...] Nerses Berberian², Arthur Sahakyan², Roksana Zakharyan¹, Arsen Arakelyan ¹

PUBLISHED 10 Jan 2017

Author details Author details

¹ Institute of Molecular Biology, Armenian National Academy of Sciences, Yerevan, Armenia
² ARTMED Medical Rehabilitation Center (CJSC), Yerevan, Armenia

OPEN PEER REVIEW

REVIEWER STATUS

Abstract

Background. Breast cancer is one of the most common cancers in women worldwide. The germline mutations of the BRCA1 and BRCA2 genes are the most significant and well characterized genetic risk factors for hereditary breast cancer. Intensive research in the last decades has demonstrated that the incidence of mutations varies widely among different populations. In this study we attempted to perform a pilot study for identification and characterization of mutations in BRCA1 and BRCA2 genes among Armenian patients with family history of breast cancer and their healthy relatives. Methods. We performed targeted exome sequencing for BRCA1 and BRCA2 genes in 6 patients and their healthy relatives. After alignment of short reads to the reference genome, germline single nucleotide variation and indel discovery was performed using GATK software. Functional implications of identified variants were assessed using ENSEMBL Variant Effect Predictor tool. Results. In total, 39 single nucleotide variations and 4 indels were identified, from which 15 SNPs and 3 indels were novel. No known pathogenic mutations were identified, but 2 SNPs causing missense amino acid mutations had significantly increased frequencies in the study group compared to the 1000 Genome populations. Conclusions. Our results demonstrate the importance of screening of BRCA1 and BRCA2 gene variants in the Armenian population in order to identity specifics of mutation spectrum and frequencies and enable accurate risk assessment of hereditary breast cancers.

Keywords

breast cancer, BRCA1, BRCA2, mutation screening, targeted exome sequencing

Corresponding author: Arsen Arakelyan

Competing interests: No competing interests were disclosed.

Grant information: This research received the grant from Armenian National Science & Education Fund (ANSEF) [#molbio-4334] to AC and SA.
The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Copyright: © 2017 Atshemyan S et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite: Atshemyan S, Chavushyan A, Berberian N et al. Characterization of BRCA1/2 mutations in patients with family history of breast cancer in Armenia [version 1; peer review: 2 approved]. F1000Research 2017, 6:29 (https://doi.org/10.12688/f1000research.10434.1) First published: 10 Jan 2017, 6:29 (https://doi.org/10.12688/f1000research.10434.1) Latest published: 10 Jan 2017, 6:29 (https://doi.org/10.12688/f1000research.10434.1)

Introduction

Breast cancer (BC) is one of the most common cancers in females worldwide¹ and particularly in Armenia². Despite the high prevalence of this disease in developed countries, it has become highly prevalent in developing countries (50% of all cancer cases) and is characterized by high mortality rate (58% of all breast cancer related deaths)³.

The germline mutations of the BRCA1⁴ and BRCA2⁵ genes are the most significant and well characterized genetic risk factors for hereditary breast cancer, which constitutes about 5–10% of all cases⁶. Inherited mutations in BRCA1 and BRCA2 genes account for 30–50% of all known mutations associated with this disease^7,8. Women who carry BRCA1 mutations are particularly susceptible to the development of breast cancer before the age of 35–40 with a probability rate of 45%–60%, whereas women who inherit a BRCA2 mutation have a 25%–40% risk of developing breast cancer^7,8. The association of BRCA1/BRCA2 gene mutations with breast cancer was first well described in Ashkenazi Jews^8–11. Intensive research in the last decades has demonstrated that the incidence of mutations in high-risk families varies widely among different populations⁶. For example, the mutations in BRCA1 and BRCA2 were each estimated to account for 45–50% of families with multiple cases of breast and ovarian cancer in UK and USA^3,12, whereas mutation prevalence among African–Americans with family breast and ovarian cancer history was 16.3% for BRCA1 and 11.3–14.4% for BRCA2^13,14, which is significantly lower compared to Caucasian populations. Identification of the BRCA1/BRCA2 mutations in different populations and ethnic groups is an important endeavor, which enables geneticists and oncologists to make more specific choices in genetic testing of members of high-risk families^15–17.

Here we have attempted to perform a pilot study for identification and characterization of mutations in BRCA1 and BRCA2 genes among Armenian patients with family history of breast cancer and their healthy relatives.

Materials and methods

Samples

Six patients with confirmed family history of breast cancer (at least two cases in a family) and their first-degree healthy relatives were recruited in this study (except for the BC10 patient, see Table 1). Patients were admitted to the National Center of Oncology MH RA and ARTMED Medical Rehabilitation CJSC. Written informed consent forms were obtained from all the study participants. This study was approved by the Institutional Review Board (IRB00004079) of the Institute of Molecular Biology NAS RA.

Table 1. Family structure of the studied subjects.

Family	Sample	Agec (Age at Diagnosis)	Disease	Family history	Healthy relative in study
FAM1	BC01	67	BC	Yes (daughter)	BC02
FAM1	BC02		HC
FAM2	BC03	41	BC	Yes (sister [BC04], grandmother)	BC05
FAM2	BC04	40	BC	Yes (sister [BC03], grandmother)	BC05
FAM2	BC05		HC
FAM3	BC06	42	BC	Yes (mother, aunt)	BC07
FAM3	BC07		HC
FAM4	BC08	36	BC	Yes (aunt)	BC09
FAM4	BC09		HC
FAM5	BC10	38	BC	Yes (mother, aunt)	NA

Blood samples were collected in EDTA-containing tubes and genomic DNA was extracted according to the protocol described elsewhere¹⁸. A260/A280 ratio measured for evaluation of quality and quantity of extracted DNA was in the range of 1.8–2.

Exome sequencing

BRCA1 and BRCA2 exome sequencing was performed by an external service provider (Admera Health LLC, South Plainfield, NJ, USA) using the proprietary breast cancer panel iBRCA^TM, which detects genetic variations in all exons of BRCA1 and BRCA2. According to the service provider’s description, this panel utilizes the targeted amplicon (166 amplicons) sequencing method, based on Seq-Ready™ TE Panels protocol (WaferGen Biosystems Inc, Freemont, CA, USA). Reagent cocktails and samples were aliquoted into a 384-well sample source plate. The source plate and BRCA1/2 SmartChip™ were pre-dispensed with Seq-Ready™ TE BRCA1/2 Primers and were placed into the SmartChip™ Multisample Nanodispenser. The SmartChip™ was then amplified with Bio-Rad T100 SmartChip™ TE Cycler. PCR product was then purified with Agencourt AMPure XP (Beckman Coulter, Inc.), according to manufacturer’s instructions. Samples were then quantified with Qubit® 2.0 Fluorometer (Thermo Fisher Scientific, Inc.) and quality analyzed with Tapestation (Agilent Technologies). Sequencing was performed with Illumina MySeq platform on a single lane. Raw reads for each sequenced sample were stored in separate fastq files. DNA samples were shipped on ice to avoid degradation and were passed internal quality check before processing.

Short-read alignment

For each sample, raw sequences were aligned to the human reference genome sequence (hg19, see Public genome data section) using Burrows-Wheeler Aligner (BWA) version 0.7.10 with default parameters. The resulting bam files were used in downstream variant discovery analysis.

Variant discovery

Variant discovery was performed using Genome Analysis Tool Kit (GATK) version 3.6 according to recommended workflows for germline single nucleotide variations (SNVs) and indel discovery in whole genome and exome sequencing data¹⁹. Base quality score recalibration, indel realignment and mate pair fixing were performed in bam files. Variant calling was performed without duplicate read removal. SNV and indel discovery and genotyping were performed simultaneously across all samples using standard hard filtering parameters¹⁹.

Public genome data

For the alignment, we have used the human reference genome sequence (NCBI build 36.1/hg19) from the UCSC (University of California, Santa Cruz) database (http://genome.ucsc.edu). Known SNPs (single nucleotide polymorphisms) were annotated using the UCSC database (single nucleotide polymorphism database, dbSNP version 135). 1000 Genomes phase 1 genotype data was used for human genetic variations filtration (ftp://ftp.1000genomes.ebi.ac.uk/vol1/ftp/). Allelic frequencies of detected variants were compared against 1000 Genomes phase 3 genotypes, as well as with the genome-wide association study (GWAS) data from 54 healthy Armenian females that were genotyped in the framework of population genetics study by Harber et al.²⁰ (ftp://ngs.sanger.ac.uk/scratch/project/team19/Armenian). The Data on clinically significant BRCA1 and BRCA2 variants were obtained from Breast Cancer Core DataBase maintained by National Human Genome Research Institute (https://research.nhgri.nih.gov/bic/).

Statistical analysis and functional annotation

Comparison of allele frequency distributions in the study group with 1000 Genomes and healthy Armenians was performed using Fisher’s exact test available in R 3.3.2 base package. Variant functional annotation was performed using ENSEMBL Variant Effect Predictor tool²¹.

Results

In this study we have performed exome sequencing of BRCA1 and BRCA2 genes in patients with a positive family history of breast cancer and their healthy relatives of Armenian origin. Patients’ clinical data and family structure of the studied subjects are presented in the Table 1. The aligned sequencing data is available in the NCBI Sequence Read Archive (SRA, https://www.ncbi.nlm.nih.gov/sra/) under accession SRP095082. For each sample, a total of 166 different primer pairs were used to amplify all the coding regions of BRCA1 and BRCA2 (as described in the Methods section). The average sequencing depth per base per sample was 6696±606. Detailed NGS statistics are presented in Table 2 and Supplementary file S1.

Table 2. Overall NGS statistics.

Number of samples	10
Total aligned reads (percent aligned reads)	1106492 (90%)
Target aligned reads (percent aligned reads)	1041136 (94%)
Mean coverage depth	6696.9
Target coverage at 1×	100%
Target coverage at 10×	99.99%
Target coverage at 50×	99.95%

In total, variant calling resulted in detection of 232 sequence variations (200 SNVs and 32 indels, Supplementary datasets S2 and S3). Thirty-nine SNVs and 4 indels passed the thresholds after applying hard filters (Table 3).

Table 3. Polymorphic variants in BRCA1 and BRCA2 genes in patients and their healthy relatives.

This table provides functional annotation of mutations in BRCA1 and BRCA2 genes that passed filters during variant calling with GATK.

HGVSg	Consequence	Impact	Gene	HGVSp	SIFT	PolyPhen	Clinical significance
13:g.32889775 G>T	5' UTR	MD	BRCA2	-	-	-	-
13:g.32889792 A>G	5' UTR	MD	BRCA2	-	-	-	-
13:g.32890572 G>A	5' UTR	MD	BRCA2	-	-	-	-
13:g.32899159 C>T	intronic	MD	BRCA2	-	-	-	uncertain significance, not provided
13:g.32900933 T>A	intronic	MD	BRCA2	-	-	-	-
13:g.32906729 A>C	missense	MO	BRCA2	p.Asn372His	tolerated	benign	not provided, benign
13:g.32910561 A>G	missense	MO	BRCA2	p.Glu690Gly	deleterious	possibly damaging	-
13:g.32910594 T>G	missense	MO	BRCA2	p.Phe701Cys	tolerated	benign	-
13:g.32911888 A>G	synonymous	L	BRCA2	p.Lys1132Lys	-	-	benign
13:g.32913055 A>G	synonymous	L	BRCA2	p.Leu1521Leu	-	-	benign, likely benign
13:g.32913081 A>G	missense	MO	BRCA2	p.Lys1530Arg	deleterious	probably damaging	not provided
13:g.32913609 A>G	missense	MO	BRCA2	p.Asn1706Ser	tolerated	benign	uncertain significance
13:g.32914236 C>T	missense	MO	BRCA2	p.Thr1915Met	tolerated	benign	benign
13:g.32918825 T>C	intronic	MD	BRCA2	-	-	-	-
13:g.32920905 T>C	intronic	MD	BRCA2	-	-	-	-
13:g.32929232 A>G	synonymous	L	BRCA2	p.Ser2414Ser	-	-	benign
13:g.32929387 T>C	missense	MO	BRCA2	p.Val2466Ala	tolerated	benign	uncertain significance, benign
13:g.32929451 A>G	intronic	MD	BRCA2	-	-	-	-
13:g.32936646 T>C	intronic	MD	BRCA2	-	-	-	benign
13:g.32953388 T>C	intronic	MD	BRCA2	-	-	-	not provided
13:g.32973012 A>C	3' UTR	MD	BRCA2	-	-	-	uncertain_significance
17:g.41216021 G>A	synonymous	MD	BRCA1	-	-	-	uncertain significance
17:g.41223094 T>C	missense	MO	BRCA1	p.Ser1634Gly	tolerated	benign	not provided, benign, likely benign
17:g.41226601 G>C	intronic	MD	BRCA1	-	-	-	-
17:g.41231516 C>T	intronic	MD	BRCA1	-	-	-	-
17:g.41234470 A>G	synonymous	L	BRCA1	p.Ser1436Ser	-	-	-
17:g.41244000 T>C	missense	MO	BRCA1	p.Lys1183Arg	tolerated	benign	benign
17:g.41244429 C>T	missense	MO	BRCA1	p.Ser1040Asn	tolerated	benign	benign
17:g.41244434 T>C	synonymous	L	BRCA1	p.Glu1038Glu	-	-	-
17:g.41244644 T>A	synonymous	L	BRCA1	p.Pro968Pro	-	-	-
17:g.41244936 G>A	missense	MO	BRCA1	p.Pro871Leu	tolerated	benign	not provided
17:g.41245237 A>G	synonymous	L	BRCA1	p.Leu771Leu	-	-	benign
17:g.41245466 G>A	synonymous	L	BRCA1	p.Ser694Ser	-	-	benign
17:g.41245471 C>T	missense	MO	BRCA1	p.Asp693Asn	deleterious	benign	benign
17:g.41249324 A>G	intronic	MD	BRCA1	-	-	-	uncertain significance, likely benign
17:g.41251906 T>A	intronic	MD	BRCA1	-	-	-	-
17:g.41251931 G>A	intronic	MD	BRCA1	-	-	-	not provided
17:g.41256032 G>C	intronic	MD	BRCA1	-	-	-	-
17:g.41256037 C>T	intronic	MD	BRCA1	-	-	-	-
13:g.32913172 delC	frameshift	H	BRCA2	-	-	-	-
17:g.41249364 delA	intronic	MD	BRCA1	-	-	-	-
17:g.41256076 delA	intronic	MD	BRCA1	-	-	-	-
17:g.41256087_41256101 delGAAAAAAAAAAGAAA	intronic	MD	BRCA1	-	-	-	-

HGVSg – genomic position of mutation notation by Human Genome Variation Society; Consequence – consequence of mutation; Impact – functional impact of mutation (MD – modifier, MO – moderate, L – low, H – high); HGVSp - protein sequence name notation by Human Genome Variation Society; SIFT - prediction of protein function change depending on amino acid substitution using SIFT software (http://sift.jcvi.org/); PolyPhen - prediction of protein function change depending on amino acid substitution using PolyPhen software (genetics.bwh.harvard.edu/pph2/).

From these variants, 18 were novel (15 SNV and 3 indels), and the rest have already been described in 1000 Genomes populations (Table 4). The novel variants were detected only in one or two subjects (8 in healthy relatives and 7 in patients). We identified 12 missense variants (5 in BRCA1 and 7 in BRCA2), 8 synonymous variants (5 in BRCA1 and 3 in BRCA2), 15 intronic variants (8 in BRCA1 and 7 in BRCA2) and 4 in untranslated regions of BRCA2. The frequency distributions of known BRCA1/2 variants were similar to those in 1000 Genomes populations and/or GWAS of healthy Armenians, except for the g.32914236 C>T (p_Fisher=8.35E-24 vs Armenians, p_Fisher=0.013 vs 1000 Genomes) and g.41245471 C>T (p_Fisher=0.013 vs Armenians, p_Fisher=4.7-E05). No known clinically significant variants were detected in breast cancer patients and their healthy relatives.

Table 4. Distribution of identified variants in healthy Armenians and in 1000 genomes populations.

The frequency distributions of identified mutations in the study group were compared with data from 1000 Genomes population, as well as the genome-wide association study from 54 healthy Armenian females²⁰.

HGVSg	RAF	MAF	RAF 1000 Genomes	MAF 1000 Genomes	RAF Armenians	MAF Armenians
13:g.32889775 G>T	0,90	0,10	-	-	-	-
13:g.32889792 A>G	0,75	0,25	0,85	0,15	-	-
13:g.32890572 G>A	0,75	0,25	0,79	0,21	0,73	0,27
13:g.32899159 C>T	0,95	0,05	-	-	-	-
13:g.32900933 T>A	0,70	0,30	0,65	0,35	-	-
13:g.32906729 A>C	0,70	0,30	0,75	0,25	0,72	0,28
13:g.32910561 A>G	0,95	0,05	-	-	-	-
13:g.32910594 T>G	0,95	0,05	-	-	-	-
13:g.32911888 A>G	0,70	0,30	0,73	0,27	0,71	0,29
13:g.32913055 A>G	0	1	0,03	0,97	-	-
13:g.32913081 A>G	0,95	0,05	-	-	-	-
13:g.32913609 A>G	0,95	0,05	-	-	-	-
13:g.32914236 C>T	0,90	0,10	0,99	0,01	0,98	0,02
13:g.32918825 T>C	0,95	0,05	-	-	-	-
13:g.32920905 T>C	0,95	0,05	-	-	-	-
13:g.32929232 A>G	0,70	0,30	0,77	0,23	-	-
13:g.32929387 T>C	0	1	0,02	0,98	0	1
13:g.32929451 A>G	0,95	0,05	-	-	-	-
13:g.32936646 T>C	0,60	0,40	0,47	0,53	0,64	0,36
13:g.32953388 T>C	0,60	0,40	0,49	0,51	0,63	0,37
13:g.32973012 A>C	0,60	0,40	0,84	0,16	0,70	0,30
17:g.41216021 G>A	0,95	0,05	0,99	0,01	-	-
17:g.41223094 T>C	0,50	0,50	0,64	0,36	-	-
17:g.41226601 G>C	0,50	0,50	0,65	0,35	-	-
17:g.41231516 C>T	0,50	0,50	0,65	0,35	-	-
17:g.41234470 A>G	0,50	0,50	0,66	0,34	0,55	0,45
17:g.41244000 T>C	0,50	0,50	0,65	0,35	0,55	0,45
17:g.41244429 C>T	0,95	0,05	0,99	0,01	0,98	0,02
17:g.41244434 T>C	0,90	0,10	-	-	-	-
17:g.41244644 T>A	0,95	0,05	-	-	-	-
17:g.41244936 G>A	0,50	0,50	0,46	0,54	0,53	0,47
17:g.41245237 A>G	0,50	0,50	0,66	0,34	0,55	0,45
17:g.41245466 G>A	0,50	0,50	0,66	0,34	0,61	0,39
17:g.41245471 C>T	0,70	0,30	0,97	0,03	0,92	0,08
17:g.41249324 A>G	0,95	0,05	-	-	-	-
17:g.41251906 T>A	0,95	0,05	-	-	-	-
17:g.41251931 G>A	0,80	0,20	0,90	0,10	0,78	0,22
17:g.41256032 G>C	0,95	0,05	-	-	-	-
17:g.41256037 C>T	0,95	0,05	-	-	-	-
13:g.32913172 delC	0,95	0,05	-	-	-	-
17:g.41249364 delA	0,50	0,50	0,67	0,33	-	-
17:g.41256076 delA	0,25	0,50	-	-	-	-
17:g.41256087_41256101 delGAAAAAAAAAAGAAA	0,50	0,18	-	-	-	-

MAF – minor allele frequency; RAF – reference allele frequency.

Discussion

This study provides preliminary characterization of variations in BRCA1 and BRCA2 genes in Armenian patients with family history of breast cancer. Our data suggest that no known clinically significant variants²² contribute to the disease development in these patients. Meanwhile, two other frequent mutations were identified that cause missense substitutions in coding regions of BRCA1 and BRCA2 and were predicted as having pathogenic consequence. The results of this study are in agreement with a a previous report, which also failed to identify known high risk mutations of BRCA1 and BRCA2 genes in Armenian patients using high-resolution melting PCR approach^23,24.

Mutations in BRCA1 and BRCA2 genes are known markers for hereditary breast/ovarian cancer²⁵. Currently more than 100 clinically important mutations and polymorphisms have been described. Genetic testing of these mutations was among the first included in the guidelines for cancer prognostics^3,4. Nowadays, in many countries genetic testing is routinely prescribed to patients in high-risk groups for hereditary breast and ovarian cancer^26–28. However, it has also become apparent that the distribution and appearance of particular risk alleles in BRCA1 and BRCA2 genes is population dependent, and in many cases population specific mutations are being identified^8–11. This is especially relevant to populations that have for a long time remained culturally and genetically isolated^8–11, as in the case of Armenians. Recent research has demonstrated that the genetic structure of Armenians “stabilized” about 4000 years ago and has remained almost unchanged since that time²⁰. Furthermore, our own data indicate that the frequencies of genetic variations associated with various complex human diseases share similarities both with European and Asian populations^29–31. From the other side, Armenian genomes are highly underrepresented in the current human genome sequencing initiatives and little is known about genetic predisposition to complex diseases in this particular population.

In conclusion, despite the small sample size limitation, our results demonstrate the importance of screening of BRCA1 and BRCA2 gene variants in the Armenian population in order to identity specifics of mutation spectra and frequencies and enable accurate assessment of the risk of hereditary breast cancers.

Data availability

The aligned sequencing data is available in the NCBI Sequence Read Archive (SRA) under accession number SRP095082 (https://www.ncbi.nlm.nih.gov/sra/?term=SRP095082). Scripts and vcf files with called and filtered genotypes are available: DOI, 10.5281/zenodo.215615³².

Author contributions

AA conceived the study, performed data analysis and drafted the manuscript. SA, AC and RZ performed experiments, data analysis and participated in drafting. NB and SA were responsible for patient selection, data analysis and contributed to manuscript writing.

Competing interests

No competing interests were disclosed.

Grant information

This research received the grant from Armenian National Science & Education Fund (ANSEF) [#molbio-4334] to AC and SA.

The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Supplementary materials

Supplementary File 1. Sequencing statistics: Coverage and target enrichment statistics.

This file contains details on sequencing coverage and enrichment, which were extracted from the QC report compiled by Admera Health LLC, South Plainfield, NJ, USA.

Click here to download the file.

F1000 recommended

References

1. Porter P: “Westernizing” women’s risks? Breast cancer in lower-income countries. N Engl J Med. 2008; 358(3): 213–6. PubMed Abstract | Publisher Full Text
2. Wright HZ, Simonsen K, Cheng Y: High breast cancer-related mortality in Armenia: Examining the breast cancer knowledge gap. Ann Glob Health. 2014; 80(3): 230. Publisher Full Text
3. Ferlay J, Soerjomataram I, Ervik M, et al.: GLOBOCAN 2012 v1.0. Cancer Incidence and Mortality Worldwide: IARC CancerBase No. 11 [Internet]. International Agency for Research on Cancer. Lyon, France: 2013. Reference Source
4. Miki Y, Swensen J, Shattuck-Eidens D, et al.: A strong candidate for the breast and ovarian cancer susceptibility gene BRCA1. Science. 1994; 266(5182): 66–71. PubMed Abstract | Publisher Full Text
5. Wooster R, Neuhausen SL, Mangion J, et al.: Localization of a breast cancer susceptibility gene, BRCA2, to chromosome 13q12-13. Science. 1994; 265(5181): 2088–90. PubMed Abstract | Publisher Full Text
6a. Szabo CI, King MC: Population genetics of BRCA1 and BRCA2. Am J Hum Genet. 1997; 60(5): 1013–20. PubMed Abstract | Free Full Text
6b. Ferla R, Calò V, Cascio S, et al.: Founder mutations in BRCA1 and BRCA2 genes. Ann Oncol. 2007; 18(Suppl 6): vi93–vi98. PubMed Abstract | Publisher Full Text
7. Antoniou A, Pharoah PD, Narod S, et al.: Average risks of breast and ovarian cancer associated with BRCA1 or BRCA2 mutations detected in case series unselected for family history: A combined analysis of 22 studies. Am J Hum Genet. 2003; 72(5): 1117–30. PubMed Abstract | Publisher Full Text | Free Full Text
8. Neuhausen S, Gilewski T, Norton L, et al.: Recurrent BRCA2 6174delT mutations in Ashkenazi Jewish women affected by breast cancer. Nat Genet. 1996; 13(1): 126–8. PubMed Abstract | Publisher Full Text
9. Friedman LS, Szabo CI, Ostermeyer EA, et al.: Novel inherited mutations and variable expressivity of BRCA1 alleles, including the founder mutation 185delAG in Ashkenazi Jewish families. Am J Hum Genet. 1995; 57(6): 1284–97. PubMed Abstract | Free Full Text
10. Tonin P, Weber B, Offit K, et al.: Frequency of recurrent BRCA1 and BRCA2 mutations in Ashkenazi Jewish breast cancer families. Nat Med. 1996; 2(11): 1179–83. PubMed Abstract | Publisher Full Text
11. Neuhausen S, Gilewski T, Norton L, et al.: Recurrent BRCA2 6174delT mutations in Ashkenazi Jewish women affected by breast cancer. Nat Genet. 1996; 13(1): 126–8. PubMed Abstract | Publisher Full Text
12. Easton DF, Bishop DT, Ford D, et al.: Genetic linkage analysis in familial breast and ovarian cancer: results from 214 families. The Breast Cancer Linkage Consortium. Am J Hum Genet. 1993; 52(4): 678–701. PubMed Abstract | Free Full Text
13. Nanda R, Schumm LP, Cummings S, et al.: Genetic testing in an ethnically diverse cohort of high-risk women: a comparative analysis of BRCA1 and BRCA2 mutations in American families of European and African ancestry. JAMA. 2005; 294(15): 1925–33. PubMed Abstract | Publisher Full Text
14. Huo D, Senie RT, Daly M, et al.: Prediction of BRCA Mutations Using the BRCAPRO Model in Clinic-Based African American, Hispanic, and Other Minority Families in the United States. J Clin Oncol. 2009; 27(8): 1184–90. PubMed Abstract | Publisher Full Text | Free Full Text
15. Kauff ND, Domchek SM, Friebel TM, et al.: Risk-reducing salpingo-oophorectomy for the prevention of BRCA1- and BRCA2-associated breast and gynecologic cancer: a multicenter, prospective study. J Clin Oncol. 2008; 26(8): 1331–37. PubMed Abstract | Publisher Full Text | Free Full Text
16. Rebbeck TR, Kauff ND, Domchek SM: Meta-analysis of risk reduction estimates associated with risk-reducing salpingo-oophorectomy in BRCA1 or BRCA2 mutation carriers. J Natl Cancer Inst. 2009; 101(2): 80–7. PubMed Abstract | Publisher Full Text | Free Full Text
17. Evans DG, Baildam AD, Anderson E, et al.: Risk reducing mastectomy: outcomes in 10 European centres. J Med Genet. 2009; 46(4): 254–8. PubMed Abstract | Publisher Full Text
18. Sambrook J, Russell DW: Molecular Cloning: A Laboratory Manual. 3rd ed. New York: Cold Spring Harbor Laboratory Press; 2001. Reference Source
19. Van der Auwera GA, Carneiro M, Hartl C, et al.: From FastQ data to high confidence variant calls: the Genome Analysis Toolkit best practices pipeline. Curr Protoc Bioinformatics. 2013; 43(1110): 11.10.1–33. PubMed Abstract | Publisher Full Text | Free Full Text
20. Haber M, Mezzavilla M, Xue Y, et al.: Genetic evidence for an origin of the Armenians from Bronze Age mixing of multiple populations. Eur J Hum Genet. 2016; 24(6): 931–6. PubMed Abstract | Publisher Full Text | Free Full Text
21. McLaren W, Gil L, Hunt SE, et al.: The Ensembl Variant Effect Predictor. Genome Biol. 2016; 17(1): 122. PubMed Abstract | Publisher Full Text | Free Full Text
22. Peto J, Collins N, Barfoot R, et al.: Prevalence of BRCA1 and BRCA2 gene mutations in patients with early-onset breast cancer. J Natl Cancer Inst. 1999; 91(11): 943–9. PubMed Abstract | Publisher Full Text
23. Babikyan DT, Sarkisian TF: Preliminary genetic investigation of high-risk breast cacner patients in Armenia. Eur J Hum Genet. 2009; 17(2): 191.
24. Mkrtchyan AG: Genetic analisys of hereditary breast cancer (review). Proceedings of Yerevan State Medical University post-graduate students research. 2004; 2: 73–80.
25. Janavičius R: Founder BRCA1/2 mutations in the Europe: implications for hereditary breast-ovarian cancer prevention and control. EPMA J. 2010; 1(3): 397–412. PubMed Abstract | Publisher Full Text | Free Full Text
26. George A, Riddell D, Seal S, et al.: Implementing rapid, robust, cost-effective, patient-centred, routine genetic testing in ovarian cancer patients. Sci Rep. 2016; 6: 29506. PubMed Abstract | Publisher Full Text | Free Full Text
27. Jung J, Kang E, Gwak JM, et al.: Association between basal-like phenotype and BRCA1/2 germline mutations in Korean breast cancer patients. Curr Oncol. 2016; 23(5): 298–303. PubMed Abstract | Publisher Full Text | Free Full Text
28. Nakamura S, Kwong A, Kim SW, et al.: Current Status of the Management of Hereditary Breast and Ovarian Cancer in Asia: First Report by the Asian BRCA Consortium. Public Health Genomics. 2016; 19(1): 53–60. PubMed Abstract | Publisher Full Text
29. de Bruin MA, Kwong A, Goldstein BA, et al.: Breast cancer risk factors differ between Asian and white women with BRCA1/2 mutations. Fam Cancer. 2012; 11(3): 429–39. PubMed Abstract | Publisher Full Text
30. Nair AK, Baier LJ: Complex Genetics of Type 2 Diabetes and Effect Size: What have We Learned from Isolated Populations? Rev Diabet Stud. 2015; 12(3–4): 299–319. PubMed Abstract | Publisher Full Text
31. Zou WB, Boulling A, Masamune A, et al.: No Association Between CEL-HYB Hybrid Allele and Chronic Pancreatitis in Asian Populations. Gastroenterology. 2016; 150(7): 1558–60.e5. PubMed Abstract | Publisher Full Text
32. Arakelyan A: Raw BRCA1/2 variants in breast cancer patients and healthy relatives produced with GATK. [Data set]. Zenodo. 2016. Data Source

Comments on this article Comments (0)

Version 1

VERSION 1 PUBLISHED 10 Jan 2017

Author details Author details

¹ Institute of Molecular Biology, Armenian National Academy of Sciences, Yerevan, Armenia
² ARTMED Medical Rehabilitation Center (CJSC), Yerevan, Armenia

Competing interests

No competing interests were disclosed.

Grant information

This research received the grant from Armenian National Science & Education Fund (ANSEF) [#molbio-4334] to AC and SA.
The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Article Versions (1)

version 1

Published: 10 Jan 2017, 6:29

https://doi.org/10.12688/f1000research.10434.1

Copyright

© 2017 Atshemyan S et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Download

Export To

metrics

	Views	Downloads
F1000Research	-	-
PubMed Central Data from PMC are received and updated monthly.	-	-

Citations

0

SEE MORE DETAILS

CITE

how to cite this article

Atshemyan S, Chavushyan A, Berberian N et al. Characterization of BRCA1/2 mutations in patients with family history of breast cancer in Armenia [version 1; peer review: 2 approved]. F1000Research 2017, 6:29 (https://doi.org/10.12688/f1000research.10434.1)

NOTE: If applicable, it is important to ensure the information in square brackets after the title is included in all citations of this article.

track

receive updates on this article

Track an article to receive email alerts on any updates to this article.

Open Peer Review

Version 1

VERSION 1

PUBLISHED 10 Jan 2017

Views

7

Reviewer Report 17 Mar 2017

Lusine Nazaryan-Petersen, Department of Cellular and Molecular Medicine (ICMM), Faculty of Health Science, University of Copenhagen, Copenhagen, Denmark

Approved

https://doi.org/10.5256/f1000research.11244.r21044

This is a very important pilot study characterizing variations in BRCA1 and BRCA2 genes in Armenian patients with family history of breast cancer. It provides a good background for further large-scale study in Armenia.

I have few notes to consider: ... Continue reading

This is a very important pilot study characterizing variations in BRCA1 and BRCA2 genes in Armenian patients with family history of breast cancer. It provides a good background for further large-scale study in Armenia.

I have few notes to consider:

In the section Methods, the authors used 1000 Genomes phase 1 genotype data for variations filtration. Is there any reason why they prefer phase1 data but not phase 3, which they used for assessing allelic frequencies?
I noticed that the authors did not verify the NGS detected variants by other methods, e.g. by Sanger sequencing. It is especially important to confirm the detected novel mutations to exclude that they could be false positive.
In the Table 3, the authors report a frameshift variant 13:g.32913172delC, which has a high functional impact on BRCA2. Is it detected in a patient or in a healthy relative? Could it be a novel mutation specific for the Armenian population? It is known that PolyPhen and SIFT may fail to predict the impact for some variants. The authors might consider to verify this mutation by other methods, e.g. Sanger sequencing, and report it to the appropriate databases. I would suggest to mention about this variant under the section Discussion.

Competing Interests: No competing interests were disclosed.

I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard.

CITE

Report a concern

Respond or Comment

Views

9

Reviewer Report 01 Mar 2017

David A. Goukassian, PhD Center for Translational Medicine, Temple University School of Medicine, Philadelphia, PA, USA

Approved

https://doi.org/10.5256/f1000research.11244.r20365

The breast cancer is an important health problem in Armenia and identifying specific genetic factors that may predispose to breast cancer development, especially in the families of patient that were already diagnosed with this condition may improve significantly the dire ... Continue reading

The breast cancer is an important health problem in Armenia and identifying specific genetic factors that may predispose to breast cancer development, especially in the families of patient that were already diagnosed with this condition may improve significantly the dire situation with breast cancer prevention in Armenia. Although in the small number of patients and family members, the manuscript presents a good step forward and sets an example how genetic studies in the larger cohort of breast cancer patients and members of their families could identify clinically relevant variants in BRCA1/2 gene mutations that known elsewhere outside of Armenian population as well identify variants that could be specific for Armenian populations only.

The title of manuscript is appropriate and the abstract summarizes well the reported findings. Study design is appropriate, albeit with small number of patients. Materials and methods and data analyses are suitable for the design and conclusion are justified. Methodology provides sufficient information and references for replication of the experiments as well as to build up the data base with the larger cohort of patients and their family members.

A few suggestions to make the discussion of the results better:

In this study there were no known clinically relevant variants identified. Could this be because of small number of patients in addition to the conceived notion that genetic structure of Armenians is “stabilized” 4000 years ago? Could there be other predisposing factors, as well? Needs a bit more discussion.
What is the value of the novel variants identified in this study? Could these novel variants be specific for Armenian population? Are there any other “close ethnic groups” that have shown novel variants that are not clinically relevant for the "mainstream population" but became relevant for the specific ethnic group. Brief discussion will suffice.

Competing Interests: No competing interests were disclosed.

I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard.

CITE

Report a concern

Respond or Comment

Comments on this article Comments (0)

Version 1

VERSION 1 PUBLISHED 10 Jan 2017

Open Peer Review

Reviewer Status

Reviewer Reports

	Invited Reviewers
	1	2
Version 1 10 Jan 17	read	read

David A. Goukassian, Temple University School of Medicine, Philadelphia, USA
Lusine Nazaryan-Petersen, University of Copenhagen, Copenhagen, Denmark

Comments on this article

All Comments(0)

Add a comment

Sign up for content alerts

Browse by related subjects

Back to all reports

Reviewer Report

7 Views

17 Mar 2017 | for Version 1

Lusine Nazaryan-Petersen, Department of Cellular and Molecular Medicine (ICMM), Faculty of Health Science, University of Copenhagen, Copenhagen, Denmark

7 Views Cite this report Responses(0)

Approved

This is a very important pilot study characterizing variations in BRCA1 and BRCA2 genes in Armenian patients with family history of breast cancer. It provides a good background for further large-scale study in Armenia.

I have few notes to consider:

In the section Methods, the authors used 1000 Genomes phase 1 genotype data for variations filtration. Is there any reason why they prefer phase1 data but not phase 3, which they used for assessing allelic frequencies?
I noticed that the authors did not verify the NGS detected variants by other methods, e.g. by Sanger sequencing. It is especially important to confirm the detected novel mutations to exclude that they could be false positive.
In the Table 3, the authors report a frameshift variant 13:g.32913172delC, which has a high functional impact on BRCA2. Is it detected in a patient or in a healthy relative? Could it be a novel mutation specific for the Armenian population? It is known that PolyPhen and SIFT may fail to predict the impact for some variants. The authors might consider to verify this mutation by other methods, e.g. Sanger sequencing, and report it to the appropriate databases. I would suggest to mention about this variant under the section Discussion.

Competing Interests

No competing interests were disclosed.

I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard.

Respond to this report

Responses (0)

Back to all reports

Reviewer Report

9 Views

01 Mar 2017 | for Version 1

David A. Goukassian, PhD Center for Translational Medicine, Temple University School of Medicine, Philadelphia, PA, USA

9 Views Cite this report Responses(0)

Approved

The breast cancer is an important health problem in Armenia and identifying specific genetic factors that may predispose to breast cancer development, especially in the families of patient that were already diagnosed with this condition may improve significantly the dire situation with breast cancer prevention in Armenia. Although in the small number of patients and family members, the manuscript presents a good step forward and sets an example how genetic studies in the larger cohort of breast cancer patients and members of their families could identify clinically relevant variants in BRCA1/2 gene mutations that known elsewhere outside of Armenian population as well identify variants that could be specific for Armenian populations only.

The title of manuscript is appropriate and the abstract summarizes well the reported findings. Study design is appropriate, albeit with small number of patients. Materials and methods and data analyses are suitable for the design and conclusion are justified. Methodology provides sufficient information and references for replication of the experiments as well as to build up the data base with the larger cohort of patients and their family members.

A few suggestions to make the discussion of the results better:

In this study there were no known clinically relevant variants identified. Could this be because of small number of patients in addition to the conceived notion that genetic structure of Armenians is “stabilized” 4000 years ago? Could there be other predisposing factors, as well? Needs a bit more discussion.
What is the value of the novel variants identified in this study? Could these novel variants be specific for Armenian population? Are there any other “close ethnic groups” that have shown novel variants that are not clinically relevant for the "mainstream population" but became relevant for the specific ethnic group. Brief discussion will suffice.

Competing Interests

No competing interests were disclosed.

I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard.

Respond to this report

Responses (0)

[1] 1. Porter P: “Westernizing” women’s risks? Breast cancer in lower-income countries. N Engl J Med. 2008; 358(3): 213–6. PubMed Abstract | Publisher Full Text

[2] 2. Wright HZ, Simonsen K, Cheng Y: High breast cancer-related mortality in Armenia: Examining the breast cancer knowledge gap. Ann Glob Health. 2014; 80(3): 230. Publisher Full Text

[3] 3. Ferlay J, Soerjomataram I, Ervik M, et al.: GLOBOCAN 2012 v1.0. Cancer Incidence and Mortality Worldwide: IARC CancerBase No. 11 [Internet]. International Agency for Research on Cancer. Lyon, France: 2013. Reference Source

[4] 4. Miki Y, Swensen J, Shattuck-Eidens D, et al.: A strong candidate for the breast and ovarian cancer susceptibility gene BRCA1. Science. 1994; 266(5182): 66–71. PubMed Abstract | Publisher Full Text

[5] 5. Wooster R, Neuhausen SL, Mangion J, et al.: Localization of a breast cancer susceptibility gene, BRCA2, to chromosome 13q12-13. Science. 1994; 265(5181): 2088–90. PubMed Abstract | Publisher Full Text

[6] 6a. Szabo CI, King MC: Population genetics of BRCA1 and BRCA2. Am J Hum Genet. 1997; 60(5): 1013–20. PubMed Abstract | Free Full Text

[7] 6b. Ferla R, Calò V, Cascio S, et al.: Founder mutations in BRCA1 and BRCA2 genes. Ann Oncol. 2007; 18(Suppl 6): vi93–vi98. PubMed Abstract | Publisher Full Text

[8] 7. Antoniou A, Pharoah PD, Narod S, et al.: Average risks of breast and ovarian cancer associated with BRCA1 or BRCA2 mutations detected in case series unselected for family history: A combined analysis of 22 studies. Am J Hum Genet. 2003; 72(5): 1117–30. PubMed Abstract | Publisher Full Text | Free Full Text

[9] 8. Neuhausen S, Gilewski T, Norton L, et al.: Recurrent BRCA2 6174delT mutations in Ashkenazi Jewish women affected by breast cancer. Nat Genet. 1996; 13(1): 126–8. PubMed Abstract | Publisher Full Text

[10] 9. Friedman LS, Szabo CI, Ostermeyer EA, et al.: Novel inherited mutations and variable expressivity of BRCA1 alleles, including the founder mutation 185delAG in Ashkenazi Jewish families. Am J Hum Genet. 1995; 57(6): 1284–97. PubMed Abstract | Free Full Text

[11] 10. Tonin P, Weber B, Offit K, et al.: Frequency of recurrent BRCA1 and BRCA2 mutations in Ashkenazi Jewish breast cancer families. Nat Med. 1996; 2(11): 1179–83. PubMed Abstract | Publisher Full Text

[12] 11. Neuhausen S, Gilewski T, Norton L, et al.: Recurrent BRCA2 6174delT mutations in Ashkenazi Jewish women affected by breast cancer. Nat Genet. 1996; 13(1): 126–8. PubMed Abstract | Publisher Full Text

[13] 12. Easton DF, Bishop DT, Ford D, et al.: Genetic linkage analysis in familial breast and ovarian cancer: results from 214 families. The Breast Cancer Linkage Consortium. Am J Hum Genet. 1993; 52(4): 678–701. PubMed Abstract | Free Full Text

[14] 13. Nanda R, Schumm LP, Cummings S, et al.: Genetic testing in an ethnically diverse cohort of high-risk women: a comparative analysis of BRCA1 and BRCA2 mutations in American families of European and African ancestry. JAMA. 2005; 294(15): 1925–33. PubMed Abstract | Publisher Full Text

[15] 14. Huo D, Senie RT, Daly M, et al.: Prediction of BRCA Mutations Using the BRCAPRO Model in Clinic-Based African American, Hispanic, and Other Minority Families in the United States. J Clin Oncol. 2009; 27(8): 1184–90. PubMed Abstract | Publisher Full Text | Free Full Text

[16] 15. Kauff ND, Domchek SM, Friebel TM, et al.: Risk-reducing salpingo-oophorectomy for the prevention of BRCA1- and BRCA2-associated breast and gynecologic cancer: a multicenter, prospective study. J Clin Oncol. 2008; 26(8): 1331–37. PubMed Abstract | Publisher Full Text | Free Full Text

[17] 16. Rebbeck TR, Kauff ND, Domchek SM: Meta-analysis of risk reduction estimates associated with risk-reducing salpingo-oophorectomy in BRCA1 or BRCA2 mutation carriers. J Natl Cancer Inst. 2009; 101(2): 80–7. PubMed Abstract | Publisher Full Text | Free Full Text

[18] 17. Evans DG, Baildam AD, Anderson E, et al.: Risk reducing mastectomy: outcomes in 10 European centres. J Med Genet. 2009; 46(4): 254–8. PubMed Abstract | Publisher Full Text

[19] 18. Sambrook J, Russell DW: Molecular Cloning: A Laboratory Manual. 3rd ed. New York: Cold Spring Harbor Laboratory Press; 2001. Reference Source

[20] 19. Van der Auwera GA, Carneiro M, Hartl C, et al.: From FastQ data to high confidence variant calls: the Genome Analysis Toolkit best practices pipeline. Curr Protoc Bioinformatics. 2013; 43(1110): 11.10.1–33. PubMed Abstract | Publisher Full Text | Free Full Text

[21] 20. Haber M, Mezzavilla M, Xue Y, et al.: Genetic evidence for an origin of the Armenians from Bronze Age mixing of multiple populations. Eur J Hum Genet. 2016; 24(6): 931–6. PubMed Abstract | Publisher Full Text | Free Full Text

[22] 21. McLaren W, Gil L, Hunt SE, et al.: The Ensembl Variant Effect Predictor. Genome Biol. 2016; 17(1): 122. PubMed Abstract | Publisher Full Text | Free Full Text

[23] 22. Peto J, Collins N, Barfoot R, et al.: Prevalence of BRCA1 and BRCA2 gene mutations in patients with early-onset breast cancer. J Natl Cancer Inst. 1999; 91(11): 943–9. PubMed Abstract | Publisher Full Text

[24] 23. Babikyan DT, Sarkisian TF: Preliminary genetic investigation of high-risk breast cacner patients in Armenia. Eur J Hum Genet. 2009; 17(2): 191.

[25] 24. Mkrtchyan AG: Genetic analisys of hereditary breast cancer (review). Proceedings of Yerevan State Medical University post-graduate students research. 2004; 2: 73–80.

[26] 25. Janavičius R: Founder BRCA1/2 mutations in the Europe: implications for hereditary breast-ovarian cancer prevention and control. EPMA J. 2010; 1(3): 397–412. PubMed Abstract | Publisher Full Text | Free Full Text

[27] 26. George A, Riddell D, Seal S, et al.: Implementing rapid, robust, cost-effective, patient-centred, routine genetic testing in ovarian cancer patients. Sci Rep. 2016; 6: 29506. PubMed Abstract | Publisher Full Text | Free Full Text

[28] 27. Jung J, Kang E, Gwak JM, et al.: Association between basal-like phenotype and BRCA1/2 germline mutations in Korean breast cancer patients. Curr Oncol. 2016; 23(5): 298–303. PubMed Abstract | Publisher Full Text | Free Full Text

[29] 28. Nakamura S, Kwong A, Kim SW, et al.: Current Status of the Management of Hereditary Breast and Ovarian Cancer in Asia: First Report by the Asian BRCA Consortium. Public Health Genomics. 2016; 19(1): 53–60. PubMed Abstract | Publisher Full Text

[30] 29. de Bruin MA, Kwong A, Goldstein BA, et al.: Breast cancer risk factors differ between Asian and white women with BRCA1/2 mutations. Fam Cancer. 2012; 11(3): 429–39. PubMed Abstract | Publisher Full Text

[31] 30. Nair AK, Baier LJ: Complex Genetics of Type 2 Diabetes and Effect Size: What have We Learned from Isolated Populations? Rev Diabet Stud. 2015; 12(3–4): 299–319. PubMed Abstract | Publisher Full Text

[32] 31. Zou WB, Boulling A, Masamune A, et al.: No Association Between CEL-HYB Hybrid Allele and Chronic Pancreatitis in Asian Populations. Gastroenterology. 2016; 150(7): 1558–60.e5. PubMed Abstract | Publisher Full Text

[33] 32. Arakelyan A: Raw BRCA1/2 variants in breast cancer patients and healthy relatives produced with GATK. [Data set]. Zenodo. 2016. Data Source

Characterization of BRCA1/2 mutations in patients with family history of breast cancer in Armenia

Abstract

Keywords

Introduction

Materials and methods

Samples

Table 1. Family structure of the studied subjects.

Exome sequencing

Short-read alignment

Variant discovery

Public genome data

Statistical analysis and functional annotation

Results

Table 2. Overall NGS statistics.

Table 3. Polymorphic variants in BRCA1 and BRCA2 genes in patients and their healthy relatives.

Table 4. Distribution of identified variants in healthy Armenians and in 1000 genomes populations.

Discussion

Data availability

Author contributions

Competing interests

Grant information

Supplementary materials

References

Comments on this article Comments (0)

Open Peer Review

Comments on this article Comments (0)

Open Peer Review

Reviewer Status

Reviewer Reports

Comments on this article

Browse by related subjects

Competing Interests Policy

Stay Updated