Cystic fibrosis testing in a referral laboratory: results and lessons from a six-year period
© Ridge et al.; licensee BioMed Central Ltd. 2013
Received: 31 July 2012
Accepted: 8 January 2013
Published: 23 January 2013
The recent introduction of high throughput sequencing technologies into clinical genetics has made it practical to simultaneously sequence many genes. In contrast, previous technologies limited sequencing based tests to only a handful of genes. While the ability to more accurately diagnose inherited diseases is a great benefit it introduces specific challenges. Interpretation of missense mutations continues to be challenging and the number of variants of uncertain significance continues to grow.
We leveraged the data available at ARUP Laboratories, a major reference laboratory, for the CFTR gene to explore specific challenges related to variant interpretation, including a focus on understanding ethnic-specific variants and an evaluation of existing databases for clinical interpretation of variants. In this study we analyzed 555 patients representing eight different ethnic groups. We observed 184 different variants, most of which were ethnic group specific. Eighty-five percent of these variants were present in the Cystic Fibrosis Mutation Database, whereas the Human Mutation Database and dbSNP/1000 Genomes had far fewer of the observed variants. Finally, 21 of the variants were novel and we report these variants and their clinical classifications.
Based on our analyses of data from six years of CFTR testing at ARUP Laboratories a more comprehensive, clinical grade database is needed for the accurate interpretation of observed variants. Furthermore, there is a particular need for more and better information regarding variants from individuals of non-Caucasian ethnicity.
KeywordsCystic fibrosis CFTR Novel variants Next-generation sequencing Interpretation of variants
Over the past few years, large scale sequencing efforts have provided a greater understanding of the variability of the human genome. Notably, whole genome sequencing studies have shown that each individual harbors 2.7–4.2 million single nucleotide variants (SNVs) that differ from the human reference genome , whereas exome sequencing typically identifies 18–24,000 coding region based SNVs per individual [2–7]. With regard to SNVs in coding regions, the findings generated by whole genome or exome sequencing parallel those observed when sequencing individual genes. First, a majority of identified variants are present in dbSNP [8, 9] and therefore represent more common variation. For example, exome and genome sequencing reports have shown that 88–99% of observed SNVs reside in dbSNP [1–3, 5, 10]. Second, the number of SNVs is dependent on individual genetic variability, ethnicity, and the reference sequence to which results are aligned and compared. At present, most studies utilize the human genome reference sequence for alignment (hg18/GRCh36 or hg19/GRCh37), which shares greatest similarity to Caucasian individuals of Northern European ancestry. As a consequence, the number of SNVs observed can vary considerably depending on the ethnic background of samples. Third, at the individual gene level, variants may have already been described and classified in a gene specific database. Not infrequently, however, novel SNVs are identified, even in genes that have been extensively studied through clinical research or diagnostic testing. While guidelines exist to assist in SNV annotation and functional prediction [11–14], many novel variants continue to be classified as variants of uncertain significance (VUS). As a greater number of exomes and genomes are sequenced, having a more comprehensive catalogue of human genetic variation will facilitate individual gene variant classification.
In the context of the above observations, we proposed that analysis of a dataset of variants identified in a single gene would yield insights into what will be revealed by large scale sequencing studies going forward. In our referral laboratory setting, we chose to study the cystic fibrosis transmembrane conductance regulator (CFTR) gene, representing a high volume full gene sequencing diagnostic assay. It should be noted, however, that at ARUP Laboratories, most cases have previously undergone testing with a 32-mutation panel identifying the most common disease-causing alleles, before sequencing. Thus, sequencing results are enriched for rare CFTR mutations. CFTR (NM_000492) is located at 7q31.2 and consists of 27 exons coding for a 1480 amino acid protein, which is a member of the ATP-binding cassette (ABC) transporter superfamily. Mutations in CFTR are known to result in multiple conditions, ranging from classic cystic fibrosis (CF) to monosymptomatic diseases such as congenital absence of the vas deferens, pancreatitis, or chronic bronchiectasis.
Classic CF, a recessively inherited genetic disorder, has an incidence of one in 2500–3200 in Caucasians making it one of the most common lethal genetic disorders . CF occurs with different frequencies in different ethnic groups with estimated carrier rates of one in 28, 29, 46, 65, and 90 in Caucasians, Ashkenazi Jews, Hispanics, African Americans, and Asians, respectively [16, 17]. The American College of Medical Genetics recommends carrier screening for CF in expectant individuals or those planning a pregnancy by testing for 23 known disease-causing mutations ; between 48% and 84% of clinically diagnosed CF patients have at least one of these mutations . The most common CFTR gene mutation is a three base pair deletion, p.Phe508del (prevalence of 24%-88% depending on ethnic background [17, 19–21]), which is associated with a more severe phenotype when present in a homozygous state . Similarly, other variants have variable frequency in different populations . In all ethnic groups, the majority of CFTR variants are of unknown clinical significance [22, 23]. Several databases have reported variants in CFTR including dbSNP [8, 9], the Cystic Fibrosis Mutation Database (CFMDB) , and the Human Gene Mutation Database (HGMD) . Variants in HGMD are assumed to be disease causing, but there are exceptions. Variants in dbSNP, on the other hand, are often assumed to be benign; however, that is not always the case. The CFMDB contains both disease causing and benign variants.
Herein we present results from a six-year period of CFTR diagnostic testing, including 21 novel variants, during which samples from 1407 individuals were referred to ARUP Laboratories for full gene CFTR sequencing. We focus on the need to develop a more complete understanding of variants in non-Caucasian ethnic groups, evaluate the usefulness and completeness of databases for clinical testing, and report novel variants observed at ARUP with ethnicity and clinical classifications.
Description of dataset
The dataset for the current study was comprised of variants identified through CFTR gene sequencing from 555 patients referred to the ARUP Laboratories from 2004 to March of 2011. Indications for testing included carrier testing in healthy individuals, confirmation or diagnostic testing in patient with classically affected CF, and diagnostic testing for patients with potentially CF-related symptoms, but without a diagnosis of CF. Patients self-reported ethnicity from the following categories: African American, Ashkenazi Jewish, Asian/Oriental, Caucasian, Hispanic, Mediterranean, Mid-Eastern, Native American, Other, or any combination of the above categories. All patient information (including demographic information, clinical symptoms, and laboratory reports including sequencing results) was stored in a Progeny database .
Samples were sequenced bi-directionally by a Sanger method using dye-terminator chemistry (BigDye® Direct Cycle Sequencing kit; Life Technologies, Carlsbad, CA) with M13-tailed primers. Sequencing products were separated and detected by capillary electrophoresis (ABI 3730; Life Technologies). All 27 exons of the CFTR gene and intron/exon boundaries were interrogated. The sequence was analyzed with Mutation Surveyor® (SoftGenetics, State College, PA) by two independent reviewers and a third and final review by a board certified (ABMG) clinical molecular geneticist.
Correlation with existing databases
Sorting and correlation of sequence changes between ethnic groups and with databases, and functional classifications of variants were completed using a series of custom-built Java programs. Databases queried included dbSNP build 132 [8, 9], which includes variants from the 1000 Genomes Project , CFMDB , downloaded May 2011, and the professional version of HGMD  (as of May 2011).
Within our institution, clinical molecular geneticists used the following steps to classify variants. First, relevant databases were interrogated to see if the sequence change had been previously observed, and if present, how it was reported. For this study, a local Progeny database, dbSNP, the Human Gene Mutation Database (HGMD), and the Cystic Fibrosis Mutation Database (CFMDB) were searched. Next, literature searches were performed to see if functional consequences had been previously reported for the variant. If the variant was not previously reported in the literature its frequency was determined by cross-referencing dbSNP/1000 Genomes data. In addition, in silico prediction algorithms (PolyPhen , SIFT , PSAAP , Human Splicing Finder , MaxEntScan , etc.) were used. While in silico prediction was employed in the overall evaluation of certain variants, no classification was based solely on in silico predictors. Lastly, unless evidence existed to the contrary, common, synonymous, and deep intronic (more than 20 base pairs into the intron) SNVs were typically reported as benign or suspected benign (depending on the frequency).
Total and ethnic specific variants
To determine variant distribution within the CFTR gene and in separate ethnic groups, we analyzed observed variants from 555 patient referral samples for which ethnicity data was available and at least one single nucleotide variant (SNV) was identified. The nine self-reported ethnicities selectable by patients were Mediterranean, Ashkenazi Jewish, Native American, Mid-Eastern, Asian/Oriental, Hispanic, Other/Mixed, African American, and Caucasian (in order of increasing numbers of patients).
A total of 184 different SNVs were observed in our dataset, 107 of which were observed only a single time (if a patient were homozygous for the variant, this was designated as two observations) and 16 were observed 10 or more times. The four most common variants: c.1408A > G (p.M470V), c.2562 T > G (synonymous, p.Thr854Thr), c.4389 G > A (synonymous, p.Gln1463Gln), and c.869 + 11C > T (intronic) were observed 955, 728, 427, and 236 times, respectively, and each is classified as a common polymorphism.
Shown are numbers of patients with variants, total SNVs observed in a particular ethnic group, ethnic specific SNVs (where ethnic specific SNVs are SNVs seen in only one ethnic group in this study), and SNVs that were both novel (not present in any of the three databases in Table 2 ) and ethnic specific
Novel & Ethnic Specific
Novel and database variants
The number of SNVs found in each of three databases, as well as totals is reported
SNVs in database (%)
SNVs not in database (%)
Cystic Fibrosis Mutation Database
Human Gene Mutation Database
Each of the 21 novel variants was only observed a single time in our dataset. The novel variants were classified by type with 11 variants residing in exons (8 missense, 2 nonsense, and 1 synonymous), 9 residing in introns, and a single variant upstream of the translational start site. Because these variants were confirmed bi-directionally by Sanger sequencing, which has a specificity over 99% , these variants are considered true positives. Three variants, inclusive of those mentioned, were assigned to two classes (intronic and splice site, or missense and splice site). Novel variants were classified as described in Methods with 8 variants classified as benign, 3 as suspected benign, 4 as disease causing, and 6 of unknown significance.
Functional annotation and variant significance
Number of variants
These data were derived from six years of CFTR gene sequence analysis at a national reference laboratory. While Sanger-based sequencing assays of a single gene or a few genes are contemporary practice, a shift is ongoing towards the development and implementation of larger gene panels performed by next generation sequencing. Furthermore, it is anticipated that exome and whole genome sequencing will also transition from a research to a diagnostic tool, with early examples already reported . The current study results are relevant in the context of current single-gene Sanger-based assays and larger scale gene sequencing.
Based on published whole genome and exome studies showing high percentages of observed SNVs in dbSNP, we expected the majority of SNVs in the CFTR gene to be present in dbSNP. Surprisingly, only 40% of observed SNVs in our dataset were in dbSNP. This relatively low number of observed SNVs in dbSNP can be easily explained. Common variants are more likely to be present in dbSNP than rare variants, and our dataset is enriched for rare variation. This is because, in addition to full gene sequencing, ARUP also offers a common mutation panel for CFTR, and sequenced patients are typically those with symptoms of CF, whose common mutation panel did not identify two pathogenic mutations. dbSNP is an effective catalogue for common variation, but has limited information about private variants (estimated to be thousands per individual ) and is not considered a clinical grade database.
In contrast to dbSNP, the majority of SNVs were present in CFMDB, the locus specific database. This demonstrates two key points. First, in order to interpret SNVs it is important to have sequence information for a large number of individuals. In our dataset, even after six years and more than 1000 individuals, the majority of SNVs (107 of 184) were only observed a single time. Second, locus specific mutation databases are important resources for variant interpretation as they often feature a better representation of the rare variation present at a particular locus and they generally include phenotype information about the variants.
A challenge common to sequencing studies that range from analyzing a single locus to analyzing an entire genome is variant interpretation. The first step in any study is to determine the location of variants relative to a chosen reference sequence, and the choice of reference sequence has specific implications. For example, the reference sequence may contain minor alleles in certain positions, rather than the true wild type allele. Thus variants can be identified which are, in reality, representative of the more common alleles . Next, it has been reported that most of the DNA used for sequencing of the NCBI reference genome came from a single, anonymous male donor (RPCI-11) from Buffalo, NY . In our dataset, the majority of observed SNVs were specific to certain ethnicities. For interpretation, it is necessary to determine whether these represent normal polymorphic variants within the ethnic groups or rare pathogenic alleles. This is an important distinction to make since the functional significance of a SNV may depend on the genetic background of the individual, meaning that SNVs can have different consequences in different ethnic groups [37, 38]. The importance of considering ethnic background was demonstrated in a study utilizing different sets of variant panels for CF testing . In this study, panels of 70 and 86 CF causing mutations yielded a detection rate of 85% in Caucasians and 95% in Ashkenazi Jews; however, only 58% and 62% in Hispanics and African Americans, respectively. When Hispanic and African American specific mutations were included in the panel, detection rates increased to ~95% in both groups . Additionally, by utilizing a unique reference sequence based on ethnic-specific allele frequencies, the error rate in identifying disease-associated variants can be substantially reduced .
As a result of our analyses, we have identified a number of challenges involving variant interpretation, such as identification of novel variants, choice of reference sequence, and ethnic background of the individual. These challenges are likely to extend, in greater magnitude, to gene panels, exomes, and genomes. As the scale of genomic information increases, the need for highly curated, clinical grade gene databases, such as the new ClinVar database, http://www.ncbi.nlm.nih.gov/clinvar/, will be increasingly pressing to facilitate interpretation.
Single nucleotide variants
Variants of uncertain significance
Cystic fibrosis transmembrane conductance regulator
Cystic Fibrosis Mutation Database
Human Gene Mutation Database.
The authors gratefully acknowledge the ARUP Institute for Clinical and Experimental Pathology for financial support and the medical technologists at ARUP Laboratories who performed the sequencing tests of CFTR. The ARUP Institute for Clinical and Experimental Pathology had no role in the design, interpretation, or reporting of these results.
- Moore B, Hu H, Singleton M, Reese MG, De La Vega FM, Yandell M: Global analysis of disease-related DNA sequence variation in 10 healthy individuals: implications for whole genome-based clinical diagnostics. Genet Med. 2011, 13: 210-217. 10.1097/GIM.0b013e31820ed321.PubMed CentralView ArticlePubMed
- Ng SB, Turner EH, Robertson PD, Flygare SD, Bigham AW, Lee C, Shaffer T, Wong M, Bhattacharjee A, Eichler EE, et al: Targeted capture and massively parallel sequencing of 12 human exomes. Nature. 2009, 461: 272-276. 10.1038/nature08250.PubMed CentralView ArticlePubMed
- Musunuru K, Pirruccello JP, Do R, Peloso GM, Guiducci C, Sougnez C, Garimella KV, Fisher S, Abreu J, Barry AJ, et al: Exome sequencing, ANGPTL3 mutations, and familial combined hypolipidemia. N Engl J Med. 2010, 363: 2220-2227. 10.1056/NEJMoa1002926.PubMed CentralView ArticlePubMed
- Bolze A, Byun M, McDonald D, Morgan NV, Abhyankar A, Premkumar L, Puel A, Bacon CM, Rieux-Laucat F, Pang K, et al: Whole-Exome-Sequencing-Based Discovery of Human FADD Deficiency. Am J Hum Genet. 2010, 87: 873-881. 10.1016/j.ajhg.2010.10.028.PubMed CentralView ArticlePubMed
- Byun M, Abhyankar A, Lelarge V, Plancoulaine S, Palanduz A, Telhan L, Boisson B, Picard C, Dewell S, Zhao C, et al: Whole-exome sequencing-based discovery of STIM1 deficiency in a child with fatal classic Kaposi sarcoma. J Exp Med. 2010, 207: 2307-2312. 10.1084/jem.20101597.PubMed CentralView ArticlePubMed
- Vissers LE, de Ligt J, Gilissen C, Janssen I, Steehouwer M, de Vries P, van Lier B, Arts P, Wieskamp N, Del Rosario M, et al: A de novo paradigm for mental retardation. Nat Genet. 2010, 42: 1109-1112. 10.1038/ng.712.View ArticlePubMed
- Gilissen C, Arts HH, Hoischen A, Spruijt L, Mans DA, Arts P, van Lier B, Steehouwer M, van Reeuwijk J, Kant SG, et al: Exome sequencing identifies WDR35 variants involved in Sensenbrenner syndrome. Am J Hum Genet. 2010, 87: 418-423. 10.1016/j.ajhg.2010.08.004.PubMed CentralView ArticlePubMed
- Sherry ST, Ward M, Sirotkin K: dbSNP-database for single nucleotide polymorphisms and other classes of minor genetic variation. Genome Res. 1999, 9: 677-679.PubMed
- Sherry ST, Ward MH, Kholodov M, Baker J, Phan L, Smigielski EM, Sirotkin K: dbSNP: the NCBI database of genetic variation. Nucleic Acids Res. 2001, 29: 308-311. 10.1093/nar/29.1.308.PubMed CentralView ArticlePubMed
- Bonnefond A, Durand E, Sand O, De Graeve F, Gallina S, Busiah K, Lobbens S, Simon A, Bellanne-Chantelot C, Letourneau L, et al: Molecular diagnosis of neonatal diabetes mellitus using next-generation sequencing of the whole exome. PLoS One. 2010, 5: e13630-10.1371/journal.pone.0013630.PubMed CentralView ArticlePubMed
- Richards CS, Bale S, Bellissimo DB, Das S, Grody WW, Hegde MR, Lyon E, Ward BE: ACMG recommendations for standards for interpretation and reporting of sequence variations: Revisions 2007. Genet Med. 2008, 10: 294-300. 10.1097/GIM.0b013e31816b5cae.View ArticlePubMed
- Bell J, Bodmer D, Sistermans E, Ramsden SC: Practice guidelines for the interpretation and reporting of unclassified variants (UVs) in Clinical Molecular Genetics. Clin Mol Genet Soc. 2007
- Vink GR, van Asperen CJ, Devilee P, Breuning MH, Bakker E: Unclassified variants in disease-causing genes: nonuniformity of genetic testing and counselling, a proposal for guidelines. Eur J Hum Genet. 2005, 13: 525-527. 10.1038/sj.ejhg.5201379.View ArticlePubMed
- Plon SE, Eccles DM, Easton D, Foulkes WD, Genuardi M, Greenblatt MS, Hogervorst FB, Hoogerbrugge N, Spurdle AB, Tavtigian SV: Sequence variant classification and reporting: recommendations for improving the interpretation of cancer susceptibility genetic test results. Hum Mutat. 2008, 29: 1282-1291. 10.1002/humu.20880.PubMed CentralView ArticlePubMed
- Wilson JA: Population-based carrier screening for Cystic Fibrosis - clinical, molecular, and practical considerations for the obstetrician. US Respir Dis. 2007
- Grody WW, Cutting GR, Klinger KW, Richards CS, Watson MS, Desnick RJ: Laboratory standards and guidelines for population-based cystic fibrosis carrier screening. Genet Med. 2001, 3: 149-154. 10.1097/00125817-200103000-00010.View ArticlePubMed
- Moskowitz SM, Chmiel JF, Sternen DL, Cheng E, Gibson RL, Marshall SG, Cutting GR: Clinical practice and genetic counseling for cystic fibrosis and CFTR-related disorders. Genet Med. 2008, 10: 851-868. 10.1097/GIM.0b013e31818e55a2.PubMed CentralView ArticlePubMed
- Amos J, Feldman GL, Grody WW, Monaghan K, Palomaki GE, Prior TW, Richards CS, Watson MS: Technical Standards and Guidelines for CFTR Mutation Testing. Am Coll Med Genet. 2005
- Watson MS, Cutting GR, Desnick RJ, Driscoll DA, Klinger K, Mennuti M, Palomaki GE, Popovich BW, Pratt VM, Rohlfs EM, et al: Cystic fibrosis population carrier screening: 2004 revision of American College of Medical Genetics mutation panel. Genet Med. 2004, 6: 387-391. 10.1097/01.GIM.0000139506.11694.7C.PubMed CentralView ArticlePubMed
- Alper OM, Wong LJ, Young S, Pearl M, Graham S, Sherwin J, Nussbaum E, Nielson D, Platzker A, Davies Z, et al: Identification of novel and rare mutations in California Hispanic and African American cystic fibrosis patients. Hum Mutat. 2004, 24: 353-View ArticlePubMed
- McKone EF, Emerson SS, Edwards KL, Aitken ML: Effect of genotype on phenotype and mortality in cystic fibrosis: a retrospective cohort study. Lancet. 2003, 361: 1671-1676. 10.1016/S0140-6736(03)13368-5.View ArticlePubMed
- Farrell PM, Rosenstein BJ, White TB, Accurso FJ, Castellani C, Cutting GR, Durie PR, Legrys VA, Massie J, Parad RB, et al: Guidelines for diagnosis of cystic fibrosis in newborns through older adults: Cystic Fibrosis Foundation consensus report. J Pediatr. 2008, 153: S4-S14. 10.1016/j.jpeds.2008.05.005.PubMed CentralView ArticlePubMed
- Grody WW, Cutting GR, Watson MS: The Cystic Fibrosis mutation "arms race": when less is more. Genet Med. 2007, 9: 739-744. 10.1097/GIM.0b013e318159a331.View ArticlePubMed
- Cystic Fibrosis Mutation Database. 2011, Available at http://www.genet.sickkids.on.ca/Home.html
- Stenson PD, Ball EV, Mort M, Phillips AD, Shiel JA, Thomas NS, Abeysinghe S, Krawczak M, Cooper DN: Human Gene Mutation Database (HGMD): 2003 update. Hum Mutat. 2003, 21: 577-581. 10.1002/humu.10212.View ArticlePubMed
- Genetic Pedigree Software. [http://www.progenysoftware.com]
- Consortium GP: A map of human genome variation from population-scale sequencing. Nature. 2010, 467: 1061-1073. 10.1038/nature09534.View Article
- Adzhubei IA, Schmidt S, Peshkin L, Ramensky VE, Gerasimova A, Bork P, Kondrashov AS, Sunyaev SR: A method and server for predicting damaging missense mutations. Nat Methods. 2010, 7: 248-249. 10.1038/nmeth0410-248.PubMed CentralView ArticlePubMed
- Kumar P, Henikoff S, Ng PC: Predicting the effects of coding non-synonymous variants on protein function using the SIFT algorithm. Nat Protoc. 2009, 4: 1073-1081. 10.1038/nprot.2009.86.View ArticlePubMed
- Crockett DK, Piccolo SR, Ridge PG, Margraf RL, Lyon E, Williams MS, Mitchell JA: Predicting phenotypic severity of uncertain gene variants in the RET proto-oncogene. PLoS One. 2011, 6: e18380-10.1371/journal.pone.0018380.PubMed CentralView ArticlePubMed
- Desmet FO, Hamroun D, Lalande M, Collod-Beroud G, Claustres M, Beroud C: Human Splicing Finder: an online bioinformatics tool to predict splicing signals. Nucleic Acids Res. 2009, 37: e67-10.1093/nar/gkp215.PubMed CentralView ArticlePubMed
- Yeo G, Burge CB: Maximum entropy modeling of short sequence motifs with applications to RNA splicing signals. J Comput Biol. 2004, 11: 377-394. 10.1089/1066527041410418.View ArticlePubMed
- Ewing B, Green P: Base-calling of automated sequencer traces using phred. II. Error probabilities. Genome Res. 1998, 8: 186-194.View ArticlePubMed
- Worthey EA, Mayer AN, Syverson GD, Helbling D, Bonacci BB, Decker B, Serpe JM, Dasu T, Tschannen MR, Veith RL, et al: Making a definitive diagnosis: successful clinical application of whole exome sequencing in a child with intractable inflammatory bowel disease. Genet Med. 2011, 13: 255-262. 10.1097/GIM.0b013e3182088158.View ArticlePubMed
- Cooper GM, Shendure J: Needles in stacks of needles: finding disease-causal variants in a wealth of genomic data. Nat Rev Genet. 2011, 12: 628-640. 10.1038/nrg3046.View ArticlePubMed
- Osoegawa K, Mammoser AG, Wu C, Frengen E, Zeng C, Catanese JJ, de Jong PJ: A bacterial artificial chromosome library for sequencing the complete human genome. Genome Res. 2001, 11: 483-496. 10.1101/gr.169601.PubMed CentralView ArticlePubMed
- Block R, Corsetti J, Goldenberg I, Vorobiof G, McNitt S, Ryan D, Wojciech Z, Moss AJ: The common apolipoprotein A-1 polymorphism -75A > G is associated with ethnic differences in recurrent coronary events after recovery from an acute myocardial infarction. Heart Int. 2009, 4: e8-PubMed CentralView ArticlePubMed
- Linder CC: The influence of genetic background on spontaneous and genetically engineered mouse models of complex diseases. Lab Anim (NY). 2001, 30: 34-39.
- Dewey FE, Chen R, Cordero SP, Ormond KE, Caleshu C, Karczewski KJ, Whirl-Carrillo M, Wheeler MT, Dudley JT, Byrnes JK, et al: Phased whole-genome genetic risk in a family quartet using a major allele reference sequence. PLoS Genet. 2011, 7: e1002280-10.1371/journal.pgen.1002280.PubMed CentralView ArticlePubMed
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.