- Open Access
An international effort towards developing standards for best practices in analysis, interpretation and reporting of clinical genome sequencing results in the CLARITY Challenge
Genome Biologyvolume 15, Article number: R53 (2014)
There is tremendous potential for genome sequencing to improve clinical diagnosis and care once it becomes routinely accessible, but this will require formalizing research methods into clinical best practices in the areas of sequence data generation, analysis, interpretation and reporting. The CLARITY Challenge was designed to spur convergence in methods for diagnosing genetic disease starting from clinical case history and genome sequencing data. DNA samples were obtained from three families with heritable genetic disorders and genomic sequence data were donated by sequencing platform vendors. The challenge was to analyze and interpret these data with the goals of identifying disease-causing variants and reporting the findings in a clinically useful format. Participating contestant groups were solicited broadly, and an independent panel of judges evaluated their performance.
A total of 30 international groups were engaged. The entries reveal a general convergence of practices on most elements of the analysis and interpretation process. However, even given this commonality of approach, only two groups identified the consensus candidate variants in all disease cases, demonstrating a need for consistent fine-tuning of the generally accepted methods. There was greater diversity of the final clinical report content and in the patient consenting process, demonstrating that these areas require additional exploration and standardization.
The CLARITY Challenge provides a comprehensive assessment of current practices for using genome sequencing to diagnose and report genetic diseases. There is remarkable convergence in bioinformatic techniques, but medical interpretation and reporting are areas that require further development by many groups.
The transition of genomics from research into clinical practice has begun, predicated on rapidly improving technology, data analysis methods, and more recently and importantly, standardization [1, 2]. Methods and tools for genomic diagnostics have quickly evolved to encompass all of the processes from consenting, through data generation and analysis, to interpretation, prioritization, and revisable reporting . Nonetheless, there is not currently a widely accepted set of published standards to enable the consistent and widespread use of genomics in the practice of medicine.
There have been a growing number of publicized successes in the application of genomic sequencing and interpretations for children with rare diseases of unknown etiology and patients with refractory cancers [4–11]. This has led to a growing expectation that clinical whole exome sequencing (WES) or whole genome sequencing (WGS) services will soon be standard practice for a much larger population of patients. Unlike other data-intensive diagnostic modalities, such as magnetic resonance imaging (MRI), there are no standards for the use of computational tools to analyze the outputs of different next-generation sequencing (NGS) technologies for patient care . There is a large methodological armamentarium for assembling genomic reads into a sequence, detecting variation, interpreting the clinical significance of specific sequence variants, and compiling a clinically usable report. Yet just how these methods are used in context, and in what combination, all critically impact the quality of genomically informed diagnoses. For example, many studies have utilized WES datasets essentially as large gene panels, interrogating data for only a small set of candidate genes determined based on clinical presentations , while others have utilized the entire datasets to identify and qualify mutations anywhere in the genome .
The present study was initially conceived at the 2010 Clinical Bioinformatics Summit hosted in Boston by Harvard University, the Children’s Hospital Informatics Program, and Harvard Medical School Center for Biomedical Informatics. The conference was attended by a wide range of stakeholders who discussed what it would take to attain a consistent and safe standard for clinical-grade genome-wide data interpretation. One of the consensus outcomes of this conference was the catalytic effect that a full clinical-grade genomic diagnostic challenge contest would have upon the emergence of both de facto and formal standards for genome-scale diagnostics.
This contest – dubbed the CLARITY Challenge (Children’s Leadership Award for the Reliable Interpretation and Appropriate Transmission of Your Genomic Information) – was hosted by the Manton Center for Orphan Disease Research at Boston Children’s Hospital and the Center for Biomedical Informatics at Harvard Medical School . Prizes totaling USD 25,000 were made available to the team or teams that could best analyze, interpret and report, in a clinically meaningful format, the results of parallel WES and WGS. The inspiration for CLARITY arose from the marked success of contests as a technique to focus a community on a particularly interesting and high-impact problem (e.g., various X Prizes). Successful competitions have accelerated progress in protein folding, including the MATLAB Protein Folding Contest  and the International Protein Folding Competition (CASP) , gene identification, such as EGASP , and in silico tools for predicting variant pathogenicity such as the CAGI experiment . Contests have been used to evoke ‘co-opetition’ – a collaboration centered on competition – in the hopes of crystallizing best practices and, thereby, accelerating the field. Comparative analysis is not new to this field either, as projects such as the 1000 Genomes Project  have provided the opportunity to compare technological and analytic methods across platforms and pipelines; its Exon Pilot project compared technologies from 454 Life Sciences, a Roche company (Branford, CT, USA), Applied Biosystems (Carlsbad, CA, USA), and Illumina Inc (San Diego, CA, USA), comparing capture biases, coverage fluctuations, indel alignment issues, population biases, and sequencing errors . More recently, a prominent paper compared the accuracy and sensitivity of results obtained using an Illumina Hiseq 2000 instrument and Complete Genomics’ WGS service . But there has not been a competition that has focused on the entire front-to-back process of applying NGS to patient care in a manner suitable for large-scale clinical adoption.
Admittedly, there are limitations to this method. To keep the scope of the competition manageable, it was focused largely on assessing the processes of variant annotation and subsequent medical interpretation and reporting, and no attempt was made to represent a range of clinical conditions and genetic models, or deal with the challenges of assessing clinical similarities amongst different presentations. Thus, the contest did not fully assess the real world challenges of finding causal mutations, but instead focused on comparative methods by which variants are called and assessed bioinformatically. Also outside the scope of the CLARITY Challenge are issues related to the importance of direct experimental evaluation of the functional consequence of mutation, which is a key part of the interpretation of novel variants and where improvement is also needed.
We present here a survey of the various methods used in the Challenge and summarize the opinions and attitudes of the contestants after the fact regarding the practice of clinical-grade genome-scale diagnostics for clinical practice.
Results and discussion
Three families were identified by the Manton Center for Orphan Disease Research to serve as test cases for the CLARITY challenge on the basis of having a child with clinical manifestations and/or pedigree structure suggestive of a likely genetic disease (Table 1). The clinical study reported here was performed under the auspices of the Boston Children’s Hospital Institutional Review Board (IRB) under Protocol IRB-P00000167. The organizing team worked closely with the IRB to define a protocol that protected the families’ interests, as well as the patients’ rights and prerogatives, yet allowed them to share their de-identified medical histories and DNA sequences with teams of qualified competitors around the world.
DNA samples and medical records from 12 individuals in total were collected under informed consent. Probands and their parents (i.e., trios) were enrolled from Families 1 and 3, and two affected first cousins and their parents were enrolled for Family 2. WES for all 12 participants was performed and donated by Life Technologies (Carlsbad, CA, USA), using standard protocols for the LIFE Library Builder, and sequenced with Exact Call Chemistry on SOLiD 5500xl machines. Both raw reads (XSQ format) and aligned reads (BAM format, generated with LifeScope ) were provided.
WGS for ten individuals (excluding an affected male cousin of the Family 2 proband and the cousin’s unaffected mother, for whom sufficient DNA was not available) were donated by Complete Genomics Incorporated (Mountain View, CA, USA) utilizing their standard proprietary protocols and generated using their Standard Pipeline v. 2.0. Variant call files along with aligned reads in Complete’s proprietary format, ‘masterVarBeta’, were provided.
Comprehensive clinical summaries providing clinical and diagnostic data for the presenting complaints and significant secondary findings were prepared by Manton Center staff from the primary medical records and made available on a secure server to the contestants, together with the genomic data described above.
Contestants were solicited from around the world via professional contacts, word of mouth, and an external website . Forty teams applied to participate in the Challenge, 32 of the most experienced multidisciplinary groups were invited to compete, and 30 accepted the offer. Participants – working either independently or as teams – were tasked with working toward an analysis, interpretation, and report suitable for use in a clinical setting.
At the conclusion of the Challenge, 23 teams successfully submitted entries that included descriptive reports of their bioinformatic analytical strategies with rationale, examples of data output and tables of variants, and clinical diagnostic reports for each family. Some groups also provided examples of their patient education materials, informed consent forms, preference setting documents, plans for revisable reporting, and protocols for dealing with incidental findings. Reasons given by four of the seven non-completing teams for dropping out were: technical and management issues, personnel changes within the team, inability to finish on time, or difficulty re-aligning the WES datasets (N = 1 each). The other three teams gave no reason.
The 23 completed entries represented a diverse group of approaches and treatments, with some groups focusing almost entirely on bioinformatic issues, others on clinical and ethical considerations. The most compelling entries including a detailed description of the bioinformatic pipelines coupled with clear, concise, and understandable clinical reports. Among the 23 entries, multiple genes were listed as possibly causative for all families (25 for Family 1, 42 for Family 2 and 29 for Family 3). Nevertheless, a consensus was achieved regarding probable pathogenic variants in two of the families. In Family 1, mutations of the titin gene, TTN [Online Mendelian Inheritance in Man (OMIM) 188840/603689], recently reported to cause a form of centronuclear myopathy , were identified as possibly or likely pathogenic by 8/23 groups, and 6/23 groups reported GJB2 (OMIM 121011/220290) variants as the likely cause of the hearing loss in the proband. Similarly, 13/23 groups identified and reported a variant in TRPM4 (OMIM 606936/604559)  as likely responsible for the cardiac conduction defects in Family 2. Although no convincing pathogenic variants were identified for Family 3, there were two plausible candidates requiring further study, OBSCN and TTN, mentioned by six groups each (Table 2).
Following the independent review and discussion by the panel of judges, one ‘winner’, the multi-institution team led by Brigham and Woman’s Hospital, Division of Genetics, et al. (Boston) was selected, largely on the basis of having a solid pipeline that correctly identified most of the genes judged to be likely pathogenic, as well as for having clear and concise clinical reports that were judged to be best at conveying the complex genetic information in a clinically meaningful and understandable format. Two runners-up were also cited. The first was a combined team from Genomatix (Munich, Germany), CeGaT (Tübingen, Germany) and the University Hospital of Bonn (Bonn, Germany), which had a robust pipeline that correctly identified every relevant gene in clear clinical reports. The second was a team from the Iowa Institute of Human Genetics at the University of Iowa, which had an outstanding array of patient education materials, procedures for patient preference setting and dealing with incidental findings, and policies for transfer of results of uncertain significance to an appropriate research setting if so desired by the patients. The content of the three winning entries is available as Additional files 1, 2 and 3. Five additional teams were cited for ‘honorable mention’ for having pipelines that identified one or more of the likely ‘correct’ genes and for providing clear clinical reporting (Table 3). These eight teams recognized by the judges are defined as ‘finalists’ in the text and for purposes of statistical analysis.
Criterion 1 (pipeline): what methods did each team use to analyze and interpret the genome sequences?
The particulars of the bioinformatic pipelines, variant annotation and report generation approaches employed by the contestants are summarized in Table 4.
The majority of contestants chose to use the supplied alignments of the data. This is not surprising since the read data from Complete Genomics and SOLiD require special handling due to the nature of sequencing, split reads in the former, and potential for color-space reads in the latter. However, three teams were unable to read the data formats provided and did not submit complete entries.
Alignments were recomputed for the Complete Genomics data by 5 out of 21 teams, with only one team reporting use of the aligner DNAnexus (Palo Alto, CA, USA), while 8 out of 21 teams recomputed alignments for the SOLiD data. For the SOLiD data, five teams recomputed alignments with software aware of color-space, and two teams indicated that they compared their color-space results against a base-space aligner. Reported aligners used for SOLiD data included the LifeScope aligner, BFAST , BWA [26–28], Novocraft’s novoalignCS (Selangor, Malaysia) and the Genomatix aligner (Munich, Germany), with some teams utilizing multiple tools for comparison. One team performed error correction prior to alignment for the SOLiD data using LifeScope’s SAET (SOLiD Accuracy Enhancement Tool, Carlsbad, CA, USA).
Prior to variant calling, many teams removed read duplicates using Picard  or SAMtools , while some teams omitted this step due to the danger of removing non-duplicate reads from single-end data. Using WGS and WES data together gave an additional way to account for PCR duplication. Limited quality control (QC) was performed prior to variant calling, with a single team using BEDTools  to analyze coverage QC metrics, and one other team reporting custom mapping QC filters.
O’Rawe et al. suggested that the choice of pipeline might be a significant source of variability in the outcome of NGS analyses . Of the teams, 40% used both the Gene Analysis Toolkit (GATK) [33, 34] and SAMtools  for variant calling, with the majority using at least one or the other. This indicates that while there is not complete consensus, using GATK, SAMtools or both resulted in acceptable results for the challenge. While GATK and SAMtools are the most popular variant callers used today and reported in this survey, their relative performance has been shown to vary with the sequencing depth [35, 36], and direct comparison of variant calls resulting from a parallel analysis of the same raw data by different variant-calling pipelines has revealed remarkably low concordance , leading to words of caution in interpreting individual genomes for genomic medicine.
SAMtools was used by some teams to jointly call SNPs and indels while recalibrating quality scores, while other teams used GATK to call SNPs and indels separately. Teams using GATK typically followed the Broad Institute’s best practice guidelines, performing indel realignment prior to indel calling, base quality score recalibration prior to SNP calling, and variant-calling score recalibration after variant calling. Some teams ignored GATK’s base quality score recalibration, mentioning that at the time GATK did not support SOLiD error profiles. LifeScope software containing DiBayes was also used on SOLiD data to call SNPs, and with local realignment to call small indels. In some cases, multiple variant-calling methods were used and compared, with all but one using GATK, SAMtools or some combination thereof. Other tools used with one mention each include: the DNAnexus variant caller, FreeBayes  and Avadis NGS (v1.3.1). A number of teams utilized the WGS results from Complete Genomics to look for potentially pathogenic de novo copy number variants, but none were found.
A significant source of variation among the different entries was the number of de novo mutations reported. Less than five de novo mutations per exome, and only about 75 de novo mutations per genome, are expected for each trio [38, 39], yet some groups reported much higher numbers, recognizing that many of these changes fell within areas with low or poor coverage. Groups that used a family-aware zygosity calling approach, such as the GATK module ‘Phase by Transmission’, developed much more refined lists of only a few potential de novo variants per proband, demonstrating the importance of this approach. However, several teams reported problems using the SOLiD data for this analysis as the BAM format provided by SOLiD was different from that expected by GATK, limiting the analysis to Complete Genomics data in those cases.
Variant filtering or recalibration after initial variant calls was performed by 16 out of 20 teams. Six teams used GATK variant quality score recalibration, with other teams reporting use of custom tools. Some teams used BEDTools for coverage QC metrics, but there was no consensus on tools to report sequencing and analysis QC metrics for post-alignment and variant calling.
Teams were asked if they employed any reference datasets in calling variants or comparing datasets to known variants (e.g., batched variant calls, known variant lists, etc.). The most common reference data reported included variants from the 1000 Genomes Project, dbSNP , HapMap Project , NHLBI Grand Opportunity Exome Sequencing Project (Bethesda, MD, USA), and the GATK Resource Bundle (distributed with GATK). Other reference datasets mentioned were the Mills Indel Gold Standard , NCBI ClinVar (Bethesda, MD, USA) as well as public sequencing data produced from the technologies used in this challenge.
One limitation of exome and genome sequencing is that the low/no coverage regions can lead to false positive or false negative results (sometimes 7% to 10% of the exons of the genes of interest have insufficient sequence reads to make a variant call ). Only 42% of teams quantified and reported on regions with insufficient coverage or data quality, though 50% of the finalists and two of the top three teams did.
Many clinical diagnostic protocols still require independent confirmation of NGS results, often by Sanger-based resequencing studies, to validate clinically relevant findings. Although this was not possible in the context of a competition where the contestants did not have access to DNA from the participants, 11 groups took advantage of the independently derived WES and WGS datasets to cross-check and validate their findings. In every instance except two, the teams reported concordance between the variant calls for the TTN, GJB2, and TRPM4 mutations that were considered likely pathogenic. The exceptions were both related to calls that were considered false positives in the SOLiD data due to poor quality or coverage at the GJB2 and TRPM4 loci, respectively. The GJB2 findings had previously been clinically confirmed and the contest organizers subsequently arranged for independent research and clinical testing, which confirmed the TTN and TRPM4 variants as well.
Medical interpretation of variant lists
The most frequent methods used to annotate variants reported were Annovar  (52%), in-house developed software (17%), and Ingenuity (Redwood City, CA, USA) (12%). Other tools reported were Variant Tools , KggSeq , SG-ADVISER (Scripps Genome Annotation and Distributed Variant Interpretation Server, La Jolla, CA, USA), Genome Trax (Wolfenbüttel, Germany), VAAST (Variant Annotation and Search Tool) , Omicia Opal , MapSNPs , in-house pipelines, and combinations thereof. There were a large variety of annotation sources (see Table 4), including but not limited to: OMIM , Uniprot , SeattleSeq , SNPedia , NCBI ClinVar, PharmGKB , Human Gene Mutation Database , dbNSFP , and in-house annotations. More importantly, most teams (14/20, 70%) performed their own curation of annotations, for example, by performing a medical literature review or by checking for errors in externally accessed databases. Thus, a manual review of annotations was deemed necessary by most contestants. Many teams considered the family pedigree structure as an important input for evaluating variants, as this allowed identification of potential de novo mutations, filtering for dominant inheritance in Family 2, ensuring Mendelian segregation and carrier status in parents for recessive mutations, etc. The function was largely performed manually, but use of automated tools such as the GATK module ‘Phase by Transmission’ was considered by some groups although the underlying structure of the SOLiD data led to problems with the analysis.
Reasons given for why teams did not report each of the likely pathogenic variants in Families 1 and 2 varied by gene and by team, but in many instances, were due to decisions made during the medical interpretation phase of analysis. Of the 15 teams that did not report the TTN variants for whom survey data were available, the variant calls generated by three failed to identify them. Twelve groups reported that their variant callers identified the two variants, but in six of these, automatic filters eliminated the gene from further consideration because the frequency of potentially pathogenic variants in this enormous gene was considered too high to be credible as a likely disease gene. Of the six instances where the automated pipelines reported the variants as potentially pathogenic, five were subsequently manually eliminated from further consideration because medical consultants lacked the clinical expertise or did not believe the published association with cardio- or skeletal-myopathy because of the high frequency of missense changes in the normal population. Notably, in none of the exclusions based on the high degree of heterogeneity of the gene was a distinction made between predicted truncating mutations, which are much rarer, versus more common missense changes. In one instance, a simple programming error prevented TTN from rising to the top of the candidate gene list in an automated expert system, and subsequent correction of this mistake resulted in a correct call of likely pathogenicity for the TTN variants in Family 1.
Seventeen teams reported not flagging the GJB2 mutations as likely causative for hearing loss in the proband of Family 1. Remarkably, the variant callers employed by ten teams failed to identify these changes despite the fact that seven of these teams used either GATK and/or SAMtools. Among the remaining seven teams, two ignored the findings because they were considered irrelevant to the ‘primary phenotype’ of skeletal myopathy and two reported a lack of clinical expertise necessary to recognize that hearing loss was a distinct phenotype. The remaining three teams reported that one of the previously published known pathogenic variants was automatically filtered out due to its high minor allele frequency in normal populations.
The TRPM4 variant in Family 2 was clinically reported by 13 of the 23 teams. Only two teams cited failure of their variant callers to identify this mutation, but five more reported that the variant was discarded due to poor quality data (low depth and noisy location with multiple non-reference alleles at that location in the SOLiD data) in one of more of the individuals, which led to inconsistent calls among the different affected family members. Two groups failed to recognize the likely pathogenicity of this variant; one reported it as a variant of unknown significance while the last one’s computational genetic predictive scoring simply failed to weight this gene highly enough to pass the cutoff given their entered phenotypic parameters. The remaining group identified the TRPM4 variant, but strongly favored another variant in the NOS3 gene as a better explanation for the structural heart defects.
Pathogenicity prediction of missense variants
The most common tools to tackle the problem of determining the effect of amino acid substitutions on protein function for missense mutations were SIFT  and Polyphen . While 80% of teams used both SIFT and Polyphen to predict pathogenicity, there was no significant difference in the success of the teams using both SIFT and Polyphen and those who used one or the other or some other tool entirely. Other tools listed by teams were PhyloP , likelihood ratio test scores (LRT) , MutationTaster , GERP , and in-house developed tools. Also of note: 45% of teams attempted to assess the statistical confidence of assignment of pathogenicity (63% of finalists). Methods named included custom in-house methods (N = 3), considering gene size (N = 2), utilizing known predictions of pathogenicity (N = 3) and allele frequencies (N = 2), assessing commonly mutated segments (N = 2), and using true positive and neutral datasets within a Bayesian framework (N = 1).
Use of splice prediction tools is particularly important, as approximately 14% to 15% of all hereditary disease alleles are annotated as splicing mutations . Groups that utilized a suite of splice prediction tools, such as the maximum entropy model MAXENT , ExonScan  or positional distribution analysis [64, 65], were more likely to have identified potentially pathogenic mutations, particularly in the TTN gene in Family 1.
It was well recognized by all groups that allele frequency is an important consideration in assessing pathogenicity (though specific cutoffs were not mentioned). All groups also agreed that conservation of amino acid sequence across species is useful for interpretation of missense variants. Half of the teams (63% of finalists) took advantage of the whole genomic sequences to analyze non-coding variants, but none of the teams reported potential pathogenic changes in deep intronic or intergenic regions, even for Family 3, likely largely due to the undefined and uncertain status of such variants. Of teams that reported methods for predicting pathogenicity of non-coding variants, the most frequently used methods were splicing prediction algorithms (85%) and transcription factor binding site prediction (46%), with 23% also considering changes in known promoter/enhancer elements, and one team each assessing evolutionary conservation, DNase hypersensitivity sites and microRNA-binding sites.
Medical interpretation and correlation of pathogenic variants with the clinical presentations
Almost all entrants performed a clinical correlation at the level of a single general diagnosis such as ‘myopathy’, ‘centronuclear myopathy’ or ‘nemaline myopathy’ with a list of predetermined candidate genes. From a clinical perspective, this reduces clinical diagnostic decision support to a list or panel and counts on that subset being complete for maximum sensitivity. However, in the case of Family 1, for example, the likely pathogenic gene was not generally recognized as causative for centronuclear myopathy at the time of the contest. In contrast, one entrant used clinically driven diagnostic decision support  in which the clinical analysis was carried out based on a description of the patient’s various pertinent positive and pertinent negative findings, including their age of onset. This was then paired to the genome analysis in a way that used a novel pertinence calculation to find the one or more genes among those with described phenotypes that best explains the set of pertinent positive and negative findings . As they become refined and validated, such automated approaches will become a critical aid in the future for reducing the analysis times to a manageable level necessary to support the higher throughputs required in a clinical diagnostic setting. Indeed, the reported range of person-hours per case required for medical interpretation of each case was 1 to 50 hours, with the automated approach requiring less than 4 hours on average to complete.
Attitudes and remarks
Three teams were unable to read the data formats provided and did not submit complete applications. This likely reflects the unique nature and format of SOLiD and Complete Genomics data and suggests that greater adoption of standard formats (FASTQ, SAM/BAM and VCF) for bioinformatics tools is required.
We observed that finalists were significantly more likely to express a preference for generating their own sequencing data instead of having it generated by an external sequencing provider (75% versus 27%, P = 0.041). The main reason expressed for in-house data generation was control over the sequencing process to ensure production and assessment of high quality data. Other reasons expressed included cost, turnaround time, and ability for reanalysis. This preference may also reflect a tendency for the most experienced groups to have a legacy capacity to generate sequence data, and thus a bias towards using their own capacity. However, it also raises the reasonable possibility that integrated control of the process from sequence generation through variant calling is important for producing the highest quality variant calls.
Overall, the teams when asked for reasons for their preference in their preferred sequencing technology mentioned accuracy and standardized software tools, highlighting the need for standard methods and tools for primary bioinformatics analysis. Furthermore, the majority of teams (13/18) felt that NGS should be combined with classical techniques (e.g. Sanger sequencing and PCR methods) for confirmatory testing in clinical situations. However, a few recognized that with increasing depth of coverage and accuracy of alignment, NGS, particularly of less complex libraries such as gene panels and possibly exomes, had potential to be utilized as a stand-alone test once QC studies demonstrate sufficient concordance with traditional methods.
Interestingly, all four of the finalists that did not report low-coverage or uncallable regions reported that they were going to begin doing so, whereas one of the non-finalists mentioned that they were going to add coverage quality to their reports. Regions in which sequencing technology or reference-genome-specific difficulties exist are important considerations for accurate variant detection. Moreover, it is critical to provide locations in which variant calling is not possible due to lapses in coverage.
Teams had different opinions on the level of coverage they felt was necessary for accurate variant calling from NGS of whole genomes. The finalists reported that they felt a higher level of coverage was necessary (59× average) than the rest of the teams (38× average). Similarly, the finalists differed on the coverage required for whole exomes (74× versus 49×) or gene panels (121× versus 69×).
A large majority of the teams used SIFT and Polyphen to predict the pathogenicity of a variant, which is a sound strategy given the programs do not always agree in protein predictions, and in both, specificity is reported to be high but sensitivity low .
When asked about their process used to validate pathogenicity predictions, 58% of teams reported that they did not use any validation method, or did not have any datasets to compare estimates against. The finalists were more likely to have had in-house datasets to work against, which may be due to differences in analytical resources that could be devoted to this problem. Overall, this process was reported as manual for the majority of the teams.
The diversity of approaches to preparing the contest entries made direct comparisons of methods difficult, so the post-contest survey was designed to elicit a more homogeneous dataset. Nevertheless, several contestants neglected to respond to some of the questions, and the responses to others was variable, indicating some confusion on the part of respondents regarding the intent of the query.
Criterion 2: were the methods used efficient, scalable and replicable?
There are still some manual elements to many pipelines that inhibit scalability. For an average case, teams reported that the interpretation process ranges from 1 to 50 hours (mean 15 ± 16 hours). For the CLARITY challenge, the time spent was much greater: each case took from 1 to 200 hours (mean 63 ± 59 hours). The average CPU time required for the analyses was difficult to estimate as contestants utilized different approaches, and not every entry was normalized for the number of parallel processors, but contestants reported utilizing 306 ± 965 CPU hours per case (range 6 to 8,700 hours). Reported costs to run the pipeline also varied considerably ranging from USD 100 to USD 16,000 (average USD 3,754 ± 4,589), but some contestants were unable to calculate salary costs leading to some lower estimates. Although costs have fallen dramatically, and computational resources are becoming increasingly available, the requirement for manual curation and interpretation of variant lists remains a considerable barrier to scalability, which could inhibit widespread use of NGS exome and genome diagnostics in the clinic if well-validated and substantially automated annotation tools do not emerge.
Criterion 3: was the interpretive report produced from genomic sequencing understandable and clinically useful?
Consent and return of results
When asked about their approach to consenting and return of results in the survey, teams’ responses varied considerably. The question was irrelevant for a number of contestants (9/21) whose activities were restricted to research or contract sequencing without direct patient contact. Finalists were more likely to ask patients undergoing WES/WGS to sign a specific consent form or provide specific explanatory materials for the methodology (P = 0.057). Finalists were much more likely to detail how they were going to handle incidental (i.e., unanticipated) results (P = 0.002). However, only 35% of teams reported that their consent materials include an option for patients to express their preferences around the return of incidental results. Most teams (76%) reported that they did not provide examples of consent and/or explanatory materials for patients with their CLARITY submissions, and since patient interaction was not allowed for the challenge, a number of contestants simply considered the issue moot. However, upon reflection, many teams agreed that including consent and explanatory materials would have strengthened their entries.
Overall, it is notable that most teams’ submissions did not include specific consent and explanatory materials, did not detail a predetermined approach for handling incidental results, and did not describe any options for patient preferences. In some cases, survey responses indicated that such materials and plans are used in practice but were not included in the CLARITY Challenge submission because it was not clear that such content was in the scope of the challenge. In other cases, teams reported that they have not developed these materials and plans or they do not routinely focus on this aspect of the process. These findings highlight the fact that these components, though they are essential for the patient-facing implementation of clinical sequencing, are not consistently prioritized or highlighted by many groups involved in the clinical use of NGS.
Reporting methods were not uniform amongst teams. Reporting the accession number for cDNA reference sequences was significantly more frequent in finalists than in non-finalists (87% versus 22%, P = 0.009). However, teams did converge on some items: reporting zygosity was standard, with 88% of responding teams doing so. Reporting the genome build was also specified by 72%. That said, the genome build reporting was problematic even among the winning teams; two of the finalists submitted elegant reports, clearly stating the variants found, summarizing the location, the classification and the parental inheritance, with a short interpretation (Figure 1). However, the accession numbers reported were different: a different build was used in each report and not specified, so it would take considerable effort to discern whether the two reports were truly referring to the same variants.
Finalists were more likely to present a clinical summary report with their entry, with the trend approaching significance (100% versus 69%, P = 0.089). Perhaps in response to recently published guidelines , there was striking concordance in interpretation and reporting philosophy, with all finalist and most non-finalist teams gearing their reports towards a clinical geneticist, genetic counselor or non-geneticist clinician. Almost all teams agreed that a non-geneticist clinician should be the target audience of clinical summary reports (75% of finalists and 89% of non-finalists). Finalists were more likely to feel that their clinical summary report could be used in clinical care (100% versus 67%, P = 0.08), though there was overall agreement that it was important that NGS studies produce a clinical summary report that can be implemented in the clinic (95% ranked this as ‘important’ or ‘extremely important’). Most of the teams (80%) filtered their variant list by relevance to phenotype, with more successful teams more likely to do so (P = 0.074). All teams but one finalist (95%) agreed that filtering the variant list by relevance to phenotype is an appropriate method for communicating information to clinicians.
It is still not commonplace to consult with an expert physician during report preparation, but doing so clearly correlated with success. Only 61% of teams routinely consult with a medical doctor in a relevant disease area. Finalists were significantly more likely to involve clinicians on a regular basis (100% versus 36%, P = 0.001). Perhaps related, in their reports prior to the survey, all but one of the finalists considered the hearing loss to be a separate phenotype from the myopathy in Family 1, while only 36% of the less successful teams did (P = 0.059). Of those who considered the separate phenotype, 75% of finalists and 63% of non-finalists considered its genetic basis.
Overall convergence and agreement across the finalists
Overall concordance among the teams in the development of variant lists was remarkable given the dozens of available measurement and analytical components of NGS pipelines and the hundreds of thousands of variants harbored by the genomes of the families. Despite the many paths that could be taken, the finalists utilized much the same philosophy and tools in processing the data and generating variant calls, and there were often minimal differences between finalist and non-finalist teams in the large lists of potentially pathogenic variants. A caveat of our study design was the choice of sequencing technologies, as Illumina platforms now account for a greater proportion of clinical studies than either SOLiD or Complete Genomics-based studies. Eight groups analyzed only the SOLiD WES data and four restricted their analysis to the Complete Genomics WGS data, often because of real or perceived difficulties with converting the extensible sequence format from the SOLiD runs into generic FASTQ files that would run on BWA, or unfamiliarity with the proprietary Complete Genomics data formats. However, as many aspects of the analytical pipelines, including variant calling and annotation, pathogenicity prediction, medical interpretation and reporting methods, are platform independent, most results discussed here should be generally applicable even as sequencing technology continues to evolve.
A number of teams preferred to recompute alignments, even though vendor alignment data was supplied, showing a preference for control over the analysis process and methods, and to ensure high quality results. Furthermore, a subset of teams for the same reasons expressed a preference for generating sequencing data in-house with higher coverage.
The selection of bioinformatic tools used by the teams did not appear to differ greatly. Tools for variant calling centered on GATK and/or SAMtools. Of the teams, 80% performed variant filtering or recalibration after initial calls were made. It is difficult to evaluate the need for recomputing alignment, performing indel realignment, variant filtering, or recalibration, given the small number of samples in this exercise. Fewer teams reported regions with insufficient coverage or data quality, only 42% overall. Without this information, it is impossible to evaluate the sensitivity of any NGS-based testing, making this an area requiring further development throughout the field.
Use of reference datasets (1000 Genomes, dbSNP, HapMap, NHLBI Go ESP and OMNI), and annotation databases (OMIM, Uniprot, SeattleSeq, SNPedia, ClinVar, PharmGKB, Human Gene Mutation Database, dbNSFP and in-house annotations) revealed considerable consensus and uniformity across entries. This shows the preference for a wide variety of rich data sources to maximize power to understand how to prioritize and contextualize variants in the presence of known information. Annovar was the most common annotation tool, with Ingenuity also used frequently. SIFT and Polyphen were overwhelmingly used to predict pathogenicity of missense changes.
Supplementary analyses that were more likely to be employed by successful teams included consideration of allele frequency, conservation of amino acid sequence across species (for coding variants), use of splicing prediction algorithms, and assessment of transcription factor binding sites (non-coding variants). Finalists were more likely to have in-house datasets to validate pathogenicity estimates. The use of in-house datasets to serve as validation sets for estimates of pathogenicity shows the need for a large, publicly available database for this purpose.
Methods and results diverged more widely in the medical interpretation of the variant lists and correlation of variants with the clinical presentations and the medical literature. Nearly half of the teams rated their process to determine pathogenicity as ‘manual’, while the mean time per case was over 10 hours, underscoring the need for standardized automated processes. Some teams have made progress towards automating this process – e.g., Genomatix’s automated literature search tool; LitInspector  was noted by judges and other teams alike as being best in class. Some teams mentioned a desire to utilize such methods in their own pipelines. SimulConsult was able to determine most variants with minimal manual effort and less hours per case than average, providing a tremendous potential advantage in high throughput clinical environments. The ability to automate the genome–phenome correlations is a key capability that can make the difference between an analysis that can become part of clinical care and an analysis that is only practical in a research setting of gene discovery.
Questions of patient preference and the responsibilities of laboratories to return incidental findings are a controversial and rapidly evolving area . The team from Iowa highlighted the importance of patient preferences in defining the style of their reports. This represents an open challenge to the medical community to decide whether future reports should take into account patient preferences or defer to a more paternalistic model of clinically indicated disclosure. In terms of clinical reports and return of results, finalists were more likely to have consent or explanatory materials, and have a plan for incidental result return. Regardless, upon being surveyed, there was general agreement amongst all teams that clinical reports should be geared towards a clinical geneticist, genetic counselor, or non-geneticist clinician.
Variability of detection power
The fact that only two teams identified all the likely causative mutations, despite using generally similar approaches, demonstrates the need for consistency and rigor in approaches to variant interpretation. There is room for tuning the tradeoffs in sensitivity, specificity and number of etiologic hypotheses being tested that would benefit many teams performing NGS interpretation. Currently, there is little consensus on the thresholds used by various teams to determine pathogenicity of potential disease-causing variants. In some cases contestants explicitly excluded variants as potentially causative due to the belief that they were likely sequencing or variant-calling false positives or benign variants that, although occurring naturally, are not disease causing or not solely disease causing. Several groups, for example, noted that in Family 3 the proband carries multiple variants in the OBSCN gene, and that any diagnosis based upon variants in this gene must therefore be viewed cautiously.
The titin gene, TTN, presented a similar dilemma as multiple potentially pathogenic variants were detected in both Families 1 and 3. Nevertheless, successful teams recognized the probable causative nature of the TTN variants in Family 1 based on the fact that one was a published pathogenic change previously reported to cause dilated cardiomyopathy  and the second mutation was predicted to alter splicing. The winning team also cited a conference abstract, then available on the web  and now published , describing a parallel study of a cohort of patients with centronuclear myopathy with validated mutations in the TTN gene. Thus, the ability to correlate genomic results with emerging literature, almost in real time, provided the determining factor between making the correct call or not, and highlights the potential power of retrospectively revising reports as new research results become available: i.e., the concept of ‘revisibility’.
The two GJB2 gene variants identified as causative for sensorineural hearing loss for the proband in Family 1 had been clinically confirmed prior to the contest, but were not disclosed to the participants, and therefore served as a validated disease-causing variant set. Six groups identified and reported these mutations as likely responsible for the sensorineural hearing loss. The way teams dealt with the reported hearing loss in Family 1 is illustrative of variation in their understanding of the clinical phenotypes, as well as their views on reporting incidental findings. Two groups considered that the defect was likely part of the myopathic phenotype, while seven others considered the GJB2 mutations to be incidental, and hence did not look for or report them, because, even though the audiometry results were detailed in the clinical records, the hearing deficit was not listed as part of the primary diagnosis.
Pre-test differential diagnosis is needed
Fourteen of 19 teams reported having a medical geneticist on board and another included a physician partner, but four teams among the non-finalists did not have a medical expert. The fact that many teams did not appreciate the significance of GJB2 mutations for Patient 1 suggests that additional detailed input from medical experts reviewing the clinical data would have been beneficial, highlighting the need to have a clinician with genetics expertise involved in preparing a carefully considered pre-test differential diagnosis.
Emergence of standard of care
Implied by the convergent methods across the leading contestants is that there is a de facto consensus of experts for interpretation of NGS. This represents a signal opportunity to codify and make this consensus explicit to ensure the greater safety and accelerated commoditization of NGS. Aspects that still need attention and further development before becoming part of the standard of care include robust family-aware zygosity calling, coverage estimation and reporting, splice site prediction and analysis, and automation of genome–phenome interpretation.
While there has been rapid progress in the development and characterization of each of the individual components of the analysis, interpretation, and reporting pipeline, there is not yet a set of best practices that can be applied to the entire ‘end-to-end’ process of genomic measurement and interpretation. Genomic medicine will require such consensus and standardization to achieve widespread, routine, and reliable clinical use. While, eventually, organizations such as the American College of Medical Genetics and the College of American Pathologists will promulgate standards to be used in the management and accreditation of laboratories, it was the intention of the CLARITY challenge to help identify the emerging forerunners of such standards, and accelerate their development. The general feedback among contestants has been very positive and the stimulus for these groups and the entire industry to generate more and better tools and reports for molecular diagnosis has truly been achieved, also clearly documented by the number of participants.
In summary, the contest highlighted: a) the relative uniformity of methods employed for alignment, variant calling, and pathogenicity prediction; b) the need to continue developing publicly available reference genome databases; c) the need for more attention to coverage analysis and estimation of false negative rates for candidate genes; d) the need for greater attention to the development of clear, concise clinical reports, with common elements such as use of reference accession numbers and genome builds, consistent criteria for definition of pathogenicity (or degree of uncertainty); e) the value of input from medical experts who could correlate the reported phenotypic elements with the expanding literature on genes and gene function; and f) the importance of clinical genetics expertise in identifying candidate families for testing. Given the labor-intensive nature of variant analysis and clinical report generation, attention to automated genome–phenome analysis based on methods for literature mining and curation, as well as variant assessment, is a pressing need that will improve reproducibility and scalability of genomic-level analyses in the future.
Materials and methods
Subject recruitment and informed consent
Probands with rare medical conditions of apparent, but unknown, genetic etiology were identified through the Manton Center for Orphan Disease Research and their families were approached about participation in the contest. Every subject who provided clinical information and DNA specimens for analysis first provided informed consent through Protocol IRB-P00000167 under the supervision of the Boston Children’s Hospital IRB. Under the terms of this protocol, the distribution of the complete genome and exome sequences was restricted to contest organizers and qualified contestants, who all signed legal agreements to protect the privacy of the participants and pledges to return or destroy the sequences at the conclusion of the contest. Because of the risk of detection of incidental findings not related to the specific medical conditions identified in the clinical descriptions, and the fact that some participants might be publicly identified through publicity related to the Challenge, the IRB precluded any possibility of public dissemination of the raw genomic sequences. All clinical and molecular datasets were de-identified prior to distribution to the contestants, and any identifiers included in the contest entries and additional files are pseudonyms or codes with no relationships to the participants’ actual protected health information as defined by the HIPAA Privacy Rule of the US Department of Health and Human Services .
Contest entries were evaluated by an independent group of six judges not affiliated with the contest organizers (ISK, AHB and DMM). Judges represented a diverse array of disciplines, including computer science and bioinformatics (PN, DM Jr and PS), medical/human genetics (J Majzoub and HFW), and clinical diagnostics (EL). Judges were asked to evaluate all aspects of the entries, but to pay particular attention to their areas of expertise. Final selection of winners was achieved by consensus among the six independent judges and was largely based on evaluation of three main criteria:
What methods did each team use to analyze and interpret the genome sequences?
Were the methods used efficient, scalable and replicable?
Was the interpretive report produced from genomic sequencing understandable and clinically useful?
Although identification of the ‘correct’ likely causative mutations for each family was considered, this was not an overriding factor, especially in light of the fact that the mutations for each family were not previously known and in some cases the results remain uncertain and fall into the realm of ongoing research. As it was, multiple genes were listed as possibly causative for all families (25 for Family 1, 42 for Family 2 and 29 for Family 3).
Post-contest data collection and analysis
After the finalists and winners were declared, all teams were sent a packet including a structured survey of contestants’ methods and practices and copies of the winning three teams’ entries. The purpose of the survey was to provide uniformity in data for summarization and allow for self-assessment of each team’s entries relative to the winning entries. Of 23 groups that submitted contest entries, 21 (91%) returned the survey. A follow-up survey in response to reviewers’ suggestions resulted in a 100% response rate for the 23 contestants. The complete set of survey questions and aggregate responses are provided as Additional file 4. Statistical analyses were performed using the computing environment R  and all reported P values are from unpaired t-tests.
Children’s Leadership Award for the Reliable Interpretation and Appropriate Transmission of Your Genomic Information
Gene Analysis Toolkit
Institutional Review Board
next generation sequencing
Online Mendelian Inheritance in Man
polymerase chain reaction
single nucleotide polymorphism
whole exome sequencing
whole genome sequencing.
Gonzaga-Jauregui C, Lupski JR, Gibbs RA: Human genome sequencing in health and disease. Annu Rev Med. 2012, 63: 35-61.
Green ED, Guyer MS: Charting a course for genomic medicine from base pairs to bedside. Nature. 2011, 470: 204-213.
Biesecker LG: Opportunities and challenges for the integration of massively parallel genomic sequencing into clinical practice: lessons from the ClinSeq project. Genet Med. 2012, 14: 393-398.
Worthey EA, Mayer AN, Syverson GD, Helbling D, Bonacci BB, Decker B, Serpe JM, Dasu T, Tschannen MR, Veith RL, Basehore MJ, Broeckel U, Tomita-Mitchell A, Arca MJ, Casper JT, Margolis DA, Bick DP, Hessner MJ, Routes JM, Verbsky JW, Jacob HJ, Dimmock DP: 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.
Lupski JR, Reid JG, Gonzaga-Jauregui C, Rio Deiros D, Chen DC, Nazareth L, Bainbridge M, Dinh H, Jing C, Wheeler DA, McGuire AL, Zhang F, Stankiewicz P, Halperin JJ, Yang C, Gehman C, Guo D, Irikat RK, Tom W, Fantin NJ, Muzny DM, Gibbs RA: Whole-genome sequencing in a patient with Charcot-Marie-Tooth neuropathy. N Engl J Med. 2010, 362: 1181-1191.
Choi M, Scholl UI, Ji W, Liu T, Tikhonova IR, Zumbo P, Nayir A, Bakkaloglu A, Ozen S, Sanjad S, Nelson-Williams C, Farhi A, Mane S, Lifton RP: Genetic diagnosis by whole exome capture and massively parallel DNA sequencing. Proc Natl Acad Sci U S A. 2009, 106: 19096-19101.
Dixon-Salazar TJ, Silhavy JL, Udpa N, Schroth J, Bielas S, Schaffer AE, Olvera J, Bafna V, Zaki MS, Abdel-Salam GH, Mansour LA, Selim L, Abdel-Hadi S, Marzouki N, Ben-Omran T, Al-Saana NA, Sonmez FM, Celep F, Azam M, Hill KJ, Collazo A, Fenstermaker AG, Novarino G, Akizu N, Garimella KV, Sougnez C, Russ C, Gabriel SB, Gleeson JG: Exome sequencing can improve diagnosis and alter patient management. Sci Transl Med. 2012, 4: 138ra178-
Choi BO, Koo SK, Park MH, Rhee H, Yang SJ, Choi KG, Jung SC, Kim HS, Hyun YS, Nakhro K, Lee HJ, Woo HM, Chung KW: Exome sequencing is an efficient tool for genetic screening of Charcot-Marie-Tooth Disease. Hum Mutat. 2012, 33: 1610-1615.
Need AC, Shashi V, Hitomi Y, Schoch K, Shianna KV, McDonald MT, Meisler MH, Goldstein DB: Clinical application of exome sequencing in undiagnosed genetic conditions. J Med Genet. 2012, 49: 353-361.
Ellis MJ, Ding L, Shen D, Luo J, Suman VJ, Wallis JW, Van Tine BA, Hoog J, Goiffon RJ, Goldstein TC, Ng S, Lin L, Crowder R, Snider J, Ballman K, Weber J, Chen K, Koboldt DC, Kandoth C, Schierding WS, McMichael JF, Miller CA, Lu C, Harris CC, McLellan MD, Wendl MC, DeSchryver K, Allred DC, Esserman L, Unzeitig G, et al: Whole-genome analysis informs breast cancer response to aromatase inhibition. Nature. 2012, 486: 353-360.
O’Daniel JM, Lee K: Whole-genome and whole-exome sequencing in hereditary cancer: impact on genetic testing and counseling. Cancer J. 2012, 18: 287-292.
Oetting WS: Exome and genome analysis as a tool for disease identification and treatment: the, Human Genome Variation Society scientific meeting. Hum Mutat. 2011, 2012: 586-590.
Yang Y, Muzny DM, Reid JG, Bainbridge MN, Willis A, Ward PA, Braxton A, Beuten J, Xia F, Niu Z, Hardison M, Person R, Bekheirnia MR, Leduc MS, Kirby A, Pham P, Scull J, Wang M, Ding Y, Plon SE, Lupski JR, Beaudet AL, Gibbs RA, Eng CM: Clinical whole-exome sequencing for the diagnosis of Mendelian disorders. N Engl J Med. 2013, 369: 1502-1511.
The Boston Children’s Hospital CLARITY Challenge. [http://genes.childrenshospital.org/]
Gulley N: In praise of tweaking: a wiki-like programming contest. Interactions – Personalized Shared Devices. 2004, 11: 18-23.
Moult J, Fidelis K, Kryshtafovych A, Tramontano A: Critical assessment of methods of protein structure prediction (CASP) – round IX. Proteins. 2011, 79: 1-5.
Guigo R, Flicek P, Abril JF, Reymond A, Lagarde J, Denoeud F, Antonarakis S, Ashburner M, Bajic VB, Birney E, Castelo R, Eyras E, Ucla C, Gingeras TR, Harrow J, Hubbard T, Lewis SE, Reese MG: EGASP: the human ENCODE Genome Annotation Assessment Project. Genome Biol. 2006, 7: 1-31.
The Critical Assessment of Genome Interpretation (CAGI) experiment. [https://genomeinterpretation.org]
Abecasis GR, Auton A, Brooks LD, DePristo MA, Durbin RM, Handsaker RE, Kang HM, Marth GT, McVean GA, Genomes Project C: An integrated map of genetic variation from 1,092 human genomes. Nature. 2012, 491: 56-65.
Marth GT, Yu F, Indap AR, Garimella K, Gravel S, Leong WF, Tyler-Smith C, Bainbridge M, Blackwell T, Zheng-Bradley X, Chen Y, Challis D, Clarke L, Ball EV, Cibulskis K, Cooper DN, Fulton B, Hartl C, Koboldt D, Muzny D, Smith R, Sougnez C, Stewart C, Ward A, Yu J, Xue Y, Altshuler D, Bustamante CD, Clark AG, Daly M, et al: The functional spectrum of low-frequency coding variation. Genome Biol. 2011, 12: R84-
Lam HY, Clark MJ, Chen R, Natsoulis G, O’Huallachain M, Dewey FE, Habegger L, Ashley EA, Gerstein MB, Butte AJ, Ji HP, Snyder M: Performance comparison of whole-genome sequencing platforms. Nat Biotechnol. 2012, 30: 78-82.
Sakarya O, Breu H, Radovich M, Chen Y, Wang YN, Barbacioru C, Utiramerur S, Whitley PP, Brockman JP, Vatta P, Zhang Z, Popescu L, Muller MW, Kudlingar V, Garg N, Li CY, Kong BS, Bodeau JP, Nutter RC, Gu J, Bramlett KS, Ichikawa JK, Hyland FC, Siddiqui AS: RNA-seq mapping and detection of gene fusions with a suffix array algorithm. PLoS Comput Biol. 2012, 8: e1002464-
Ceyhan-Birsoy O, Agrawal PB, Hidalgo C, Schmitz-Abe K, DeChene ET, Swanson LC, Soemedi R, Vasli N, Iannaccone ST, Shieh PB, Shur N, Dennison JM, Lawlor MW, Laporte J, Markianos K, Fairbrother WG, Granzier H, Beggs AH: Recessive truncating titin gene, TTN, mutations presenting as centronuclear myopathy. Neurology. 2013, 81: 1205-1214.
Stallmeyer B, Zumhagen S, Denjoy I, Duthoit G, Hebert JL, Ferrer X, Maugenre S, Schmitz W, Kirchhefer U, Schulze-Bahr E, Guicheney P: Mutational spectrum in the Ca2+–activated cation channel gene TRPM4 in patients with cardiac conductance disturbances. Hum Mutat. 2012, 33: 109-117.
Homer N, Merriman B, Nelson SF: BFAST: an alignment tool for large scale genome resequencing. PLoS One. 2009, 4: e7767-
Li H, Durbin R: Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009, 25: 1754-1760.
Li H, Durbin R: Fast and accurate long-read alignment with Burrows-Wheeler transform. Bioinformatics. 2010, 26: 589-595.
Li H: Exploring single-sample SNP and INDEL calling with whole-genome de novo assembly. Bioinformatics. 2012, 28: 1838-1844.
Picard software. [http://picard.sourceforge.net]
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R, Genome Project Data Processing S: The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009, 25: 2078-2079.
Quinlan AR, Hall IM: BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010, 26: 841-842.
O’Rawe J, Jiang T, Sun G, Wu Y, Wang W, Hu J, Bodily P, Tian L, Hakonarson H, Johnson WE, Wei Z, Wang K, Lyon GJ: Low concordance of multiple variant-calling pipelines: practical implications for exome and genome sequencing. Genome Med. 2013, 5: 28-
McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, Garimella K, Altshuler D, Gabriel S, Daly M, DePristo MA: The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010, 20: 1297-1303.
DePristo MA, Banks E, Poplin R, Garimella KV, Maguire JR, Hartl C, Philippakis AA, del Angel G, Rivas MA, Hanna M, McKenna A, Fennell TJ, Kernytsky AM, Sivachenko AY, Cibulskis K, Gabriel SB, Altshuler D, Daly MJ: A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat Genet. 2011, 43: 491-498.
Liu Q, Guo Y, Li J, Long J, Zhang B, Shyr Y: Steps to ensure accuracy in genotype and SNP calling from Illumina sequencing data. BMC Genomics. 2012, 13: S8-
Nielsen R, Paul JS, Albrechtsen A, Song YS: Genotype and SNP calling from next-generation sequencing data. Nat Rev Genet. 2011, 12: 443-451.
Garrison E, Marth G: Haplotype-based variant detection from short-read sequencing. Cornell University Library. 2012, eprint arXiv:1207.3907
Neale BM, Kou Y, Liu L, Ma’ayan A, Samocha KE, Sabo A, Lin CF, Stevens C, Wang LS, Makarov V, Polak P, Yoon S, Maguire J, Crawford EL, Campbell NG, Geller ET, Valladares O, Schafer C, Liu H, Zhao T, Cai G, Lihm J, Dannenfelser R, Jabado O, Peralta Z, Nagaswamy U, Muzny D, Reid JG, Newsham I, Wu Y, et al: Patterns and rates of exonic de novo mutations in autism spectrum disorders. Nature. 2012, 485: 242-245.
Kong A, Frigge ML, Masson G, Besenbacher S, Sulem P, Magnusson G, Gudjonsson SA, Sigurdsson A, Jonasdottir A, Jonasdottir A, Wong WSW, Sigurdsson G, Walters GB, Steinberg S, Helgason H, Thorleifsson G, Gudbjartsson DF, Helgason A, Magnusson OT, Thorsteinsdottir U, Stefansson K: Rate of de novo mutations and the importance of father’s age to disease risk. Nature. 2012, 488: 471-475.
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.
International HapMap C: The International HapMap Project. Nature. 2003, 426: 789-796.
Mills RE, Pittard WS, Mullaney JM, Farooq U, Creasy TH, Mahurkar AA, Kemeza DM, Strassler DS, Ponting CP, Webber C, Devine SE: Natural genetic variation caused by small insertions and deletions in the human genome. Genome Res. 2011, 21: 830-839.
Bai R, Higgs J, Suchy S, Gibellini F, Knight M, Buchholz S, Benhamed S, Arjona D, Chinault C, Brandon R, Smaoui N, Richard G, Bale S: PCR-Based Enrichment and Next-Generation Sequencing of 101 Nuclear Genes for the Diagnosis of Mitochondrial Disorders. 2012, Washington DC, USA: United Mitochondrial Disease Foundation
Wang K, Li M, Hakonarson H: ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010, 38: e164-
San Lucas FA, Wang G, Scheet P, Peng B: Integrated annotation and analysis of genetic variants from next-generation sequencing studies with variant tools. Bioinformatics. 2012, 28: 421-422.
Li MX, Gui HS, Kwan JS, Bao SY, Sham PC: A comprehensive framework for prioritizing variants in exome sequencing studies of Mendelian diseases. Nucleic Acids Res. 2012, 40: e53-
Yandell M, Huff C, Hu H, Singleton M, Moore B, Xing J, Jorde LB, Reese MG: A probabilistic disease-gene finder for personal genomes. Genome Res. 2011, 21: 1529-1542.
Coonrod EM, Margraf RL, Russell A, Voelkerding KV, Reese MG: Clinical analysis of whole genome NGS data using the Omicia platform. Expert Rev Mol Diagn. 2013, 13: 529-540.
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.
Amberger J, Bocchini C, Hamosh A: A new face and new challenges for Online Mendelian Inheritance in Man (OMIM(R)). Hum Mutat. 2011, 32: 564-567.
UniProt C: Reorganizing the protein space at the Universal Protein Resource (UniProt). Nucleic Acids Res. 2012, 40: D71-D75.
Ng SB, Turner EH, Robertson PD, Flygare SD, Bigham AW, Lee C, Shaffer T, Wong M, Bhattacharjee A, Eichler EE, Bamshad M, Nickerson DA, Shendure J: Targeted capture and massively parallel sequencing of 12 human exomes. Nature. 2009, 461: 272-276.
Cariaso M, Lennon G: SNPedia: a wiki supporting personal genome annotation, interpretation and analysis. Nucleic Acids Res. 2012, 40: D1308-D1312.
Hewett M, Oliver DE, Rubin DL, Easton KL, Stuart JM, Altman RB, Klein TE: PharmGKB: the Pharmacogenetics Knowledge Base. Nucleic Acids Res. 2002, 30: 163-165.
Stenson PD, Mort M, Ball EV, Howells K, Phillips AD, Thomas NS, Cooper DN: The Human Gene Mutation Database: 2008 update. Genome Med. 2009, 1: 13-
Liu X, Jian X, Boerwinkle E: dbNSFP: a lightweight database of human nonsynonymous SNPs and their functional predictions. Hum Mutat. 2011, 32: 894-899.
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.
Pollard KS, Hubisz MJ, Rosenbloom KR, Siepel A: Detection of nonneutral substitution rates on mammalian phylogenies. Genome Res. 2010, 20: 110-121.
Chun S, Fay JC: Identification of deleterious mutations within three human genomes. Genome Res. 2009, 19: 1553-1561.
Schwarz JM, Rodelsperger C, Schuelke M, Seelow D: MutationTaster evaluates disease-causing potential of sequence alterations. Nat Methods. 2010, 7: 575-576.
Cooper GM, Stone EA, Asimenos G, Program NCS, Green ED, Batzoglou S, Sidow A: Distribution and intensity of constraint in mammalian genomic sequence. Genome Res. 2005, 15: 901-913.
Yeo G, Burge CB: Maximum entropy modeling of short sequence motifs with applications to RNA splicing signals. J Comput Biol. 2004, 11: 377-394.
Wang Z, Rolish ME, Yeo G, Tung V, Mawson M, Burge CB: Systematic identification and analysis of exonic splicing silencers. Cell. 2004, 119: 831-845.
Lim KH, Fairbrother WG: Spliceman – a computational web server that predicts sequence variations in pre-mRNA splicing. Bioinformatics. 2012, 28: 1031-1032.
Lim KH, Ferraris L, Filloux ME, Raphael BJ, Fairbrother WG: Using positional distribution to identify splicing elements and predict pre-mRNA processing defects in human genes. Proc Natl Acad Sci U S A. 2011, 108: 11093-11098.
Segal MM, Williams MS, Gropman AL, Torres AR, Forsyth R, Connolly AM, El-Hattab AW, Perlman SJ, Samanta D, Parikh S, Pavlakis SG, Feldman LK, Betensky RA, Gospe SM: Evidence-based decision support for neurological diagnosis reduces errors and unnecessary workup. J Child Neurol. 2013, in press (doi:10.1177/0883073813483365)
Flanagan SE, Patch AM, Ellard S: Using SIFT and PolyPhen to predict loss-of-function and gain-of-function mutations. Genet Test Mol Biomarkers. 2010, 14: 533-537.
Scheuner MT, Hilborne L, Brown J, Lubin IM: A report template for molecular genetic tests designed to improve communication between the clinician and laboratory. Genet Test Mol Biomarkers. 2012, 16: 761-769.
LitInspector software module in the Genomatix Software Suite. [http://www.litinspector.org]
Green RC, Berg JS, Grody WW, Kalia SS, Korf BR, Martin CL, McGuire AL, Nussbaum RL, O’Daniel JM, Ormond KE, Rehm HL, Watson MS, Williams MS, Biesecker LG: ACMG recommendations for reporting of incidental findings in clinical exome and genome sequencing. Genet Med. 2013, 15: 565-574.
Herman DS, Lam L, Taylor MRG, Wang L, Teekakirikul P, Christodoulou D, Conner L, DePalma SR, McDonough B, Sparks E, Teodorescu DL, Cirino AL, Banner NR, Pennell DJ, Graw S, Merlo M, Di Lenarda A, Sinagra G, Bos JM, Ackerman MJ, Mitchell RN, Murry CE, Lakdawala NK, Ho CY, Barton PJR, Cook SA, Mestroni L, Seidman JG, Seidman CE: Truncations of titin causing dilated cardiomyopathy. N Engl J Med. 2012, 366: 619-628.
Agrawal PB, Schmitz K, DeChene ET, Ceyhan O, Mercier M, Viola M, Markianos K, Beggs AH: Complete genetic analysis by whole exome sequencing of a cohort with centronuclear myopathy identifies titin gene mutations. Neuromuscul Disord. 2012, 22: 840-
HIPAA privacy rules. [http://www.hhs.gov/ocr/privacy/index.html]
Team RDC: R: A Language and Environment for Statistical Computing, Reference Index Version 2.x.x. 2005, Vienna: R Foundation for Statistical Computing
This work was supported by funds provided through the Gene Partnership and the Manton Center for Orphan Disease Research at Boston Children’s Hospital and the Center for Biomedical Informatics at Harvard Medical School and by generous donations in-kind of genomic sequencing services by Life Technologies (Carlsbad, CA, USA) and Complete Genomics (Mountain View, CA, USA). All the authors would especially like to thank the families for their participation and resulting critical contributions to this study. The study organizers would also like to thank Elizabeth Andrews and Keri Steadman for seminal contributions to the contest planning, Peter Kang for assistance tabulating patient phenotypes, Lindsay Swanson for assistance with return of results, and Elizabeth Torosian for help with manuscript preparation. This manuscript is dedicated to the memory of David Newsom.
BM is an employee of, and holds an equity stake in, Life Technologies, Inc. EL is an employee of the University of Utah, with an assignment as medical director at ARUP Laboratories and receives consulting fees from Complete Genomics. DM Jr is an employee of, and holds an equity stake in, Cerner Corp. HD, BHF, MSL, IL, HMM and HLR are employees of the Laboratory for Molecular Medicine of the Partners Healthcare Center for Personalized Genetic Medicine. HLR sits on advisory boards for BioBase, Clinical Future, Complete Genomics, GenomeQuest, Ingenuity, Knome and Omicia. BHF holds an equity stake in InVitae Corporation. C Gugenmus, A Hahn and BK are or were employees of Genomatix GmbH. MS, J Supper and M Menzel are employees of CeGaT GmbH. SB is an employee of, and holds an equity stake in, CeGaT GmbH. PF is an employee of the Children’s Hospital Reutlingen and is on the scientific advisory board of CeGaT GmbH. SP is a member of the scientific advisory board of CeGaT GmbH, Tübingen, Germany. MMS is an employee of and holds an equity stake in SimulConsult. AT holds an equity stake in Cypher Genomics, Inc. WZ was an employee of Sanofi SA. LZ, K Blair and DK are employees of, and hold equity stakes in, Seven Bridges Genomics, Inc. M Cariaso and GGL are principals of River Road Bio. PCN holds an equity stake in Illumina, Inc. DRR is an employee of Ingenuity Systems. YB and YH are employees of, and hold equity stakes in, Regeneron Pharmaceuticals, Inc. MGR and EK are employees of, and hold an equity stake in, Omicia. M Cargil, RKH and JMS are employees of, and hold equity stakes in, InVitae Corporation. ZIA is an employee of, and holds an equity stake in, Novocraft Technologies Sdn Bhd. The remaining authors have not declared any competing interests.
ISK, AHB and DMM conceived of, organized and ran the Challenge. CAB conducted and analyzed the post-challenge survey. NH, BM, TWY and LJL analyzed bioinformatic approaches. KCF coordinated the Challenge and collected data. ETD, MCT and AHB ascertained and recruited the subject families. CAB, MCT, SKS and ENP analyzed contest entries. IAH analyzed aspects of the entries regarding reporting of consent and return of results. EL, J Majzoub, PN, DM Jr, PS and HFW judged the Challenge and reported on the entries. NJM, RT, RSF, SWY and LM referred the families and provided and analyzed clinical data. SRS, IA, CACassa, PIWdB, HD, WF, LF, BHF, MAG, REH, KL, MSL, ML, IL, DGM, HMM, MFM, PPP, SR, HLR, R Soemedi, NOS, SV, J Supper, C Gugenmus, BK, AHahn, MS, M Menzel, SB, PF, MD, MB, SP, RJHS, JLA, JH, KR, VCS, EMS, TB, EAB-Z, TAB, BD, APD, DLK, TES, AES, R Sompallae, KW, AGB, EE, KM, SAM, OAS, PT, AB, CACampbell, JWH, AK, TM, KP, TW, DVD, HA, K Booth, NM, MMS, MSW, GT, PW, DC, SFB, GH, DLT, KLM, D Newsom, CRP, ATR, AM, LL, AP, BP, AT, AW, MH, AA, JML, M Magnusson, D Nilsson, HS, FT, C Gilissen, AHoischen, BVB, HY, MN, WZ, J Sager, LZ, K Blair, DK, MCariaso, GGL, AJ, SA, PCN, KSS, SK, VV, OI, EH, EF, NS, GG, JCR, JC, HCC, DM Jr, SAA, LR, DRR, FASL, MLGG, CTC, YB, YH, FF, YZ, ZW, JB, JMGL, DGL, JL, MCR, IV, MGR, FMDLV, EK, MCargill, RKH, JMS, GJL, DAS, BEB, BMM, KE, MY, HZ, LH, XC, XY, MChen, CL, CY, MG, PL, YK, ACA, ZIA, KMB, DEB, PMKG, AMI, BMK, J Majewski, CRM, JSP, SLS, MES and J Schwartzentruber were contestants who participated in the Challenge, reported on their standard practices, and contributed comments and material for the manuscript. CAB and AHB drafted the manuscript. NH, SKS, IAH, EL, J Majzoub, PN, DM Jr, PS, HFW, NJM, ISK and DMM contributed to drafting the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 4: Individual and aggregated results from questions in the structured surveys of contestants’ practices. Responses are broken down into separate sheets according to category as follows: PART A: Consenting and explanatory materials for whole exome/genome sequencing technology. PART B: About your summary clinical report. PART C: Interpretive reports. PART D: Revisible reporting. PART E: Variant identification. PART F: Data analysis. PART G: Validation of analytical tools. PART H: Methods predicting variant pathogenicity. PART I: From variants to phenotype. PART J: Overall impressions and team composition. PART K: Follow-up questions, costs and sensitivity. (XLSX 135 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.