Ariella Sasson

Filling the gaps – The generation of full genomic sequences for 15 common and well-documented HLA class I alleles using next-generation sequencing technology

Many common and well-documented (CWD) HLA alleles have only been partially characterized. The DNA sequence of these incomplete alleles, as published in the IMGT/HLA database, is most often limited to exons that code for the extracellular domains of the mature protein. Here we describe the application of next-generation sequencing technology to obtain full length genomic sequence from a single long-range PCR amplicon for 15 common and well-documented HLA Class I alleles. This technology is well suited to fill in the gaps of the current HLA allele sequence database which is largely incomplete. A more comprehensive catalog of HLA allele sequences would be beneficial in the evaluation of mismatches in transplantation, studies of population genetics, the evolution of HLAs, regulatory mechanisms and HLA expression, and issues related to the genomic organization of the MHC.

Efficient digest of high-throughput sequencing data in a reproducible report

BackgroundHigh-throughput sequencing (HTS) technologies are spearheading the accelerated development of biomedical research. Processing and summarizing the large amount of data generated by HTS presents a non-trivial challenge to bioinformatics. A commonly adopted standard is to store sequencing reads aligned to a reference genome in SAM (Sequence Alignment/Map) or BAM (Binary Alignment/Map) files. Quality control of SAM/BAM files is a critical checkpoint before downstream analysis. The goal of the current project is to facilitate and standardize this process.ResultsWe developed bamchop, a robust program to efficiently summarize key statistical metrics of HTS data stored in BAM files, and to visually present the results in a formatted report. The report documents information about various aspects of HTS data, such as sequencing quality, mapping to a reference genome, sequencing coverage, and base frequency. Bamchop uses the R language and Bioconductor packages to calculate statistical matrices and the Sweave utility and associated LaTeX markup for documentation. Bamchops efficiency and robustness were tested on BAM files generated by local sequencing facilities and the 1000 Genomes Project. Source code, instruction and example reports of bamchop are freely available from https://github.com/CBMi-BiG/bamchop.ConclusionsBamchop enables biomedical researchers to quickly and rigorously evaluate HTS data by providing a convenient synopsis and user-friendly reports.

Characterizing reduced coverage regions through comparison of exome and genome sequencing data across 10 centers

PurposeAs massively parallel sequencing is increasingly being used for clinical decision making, it has become critical to understand parameters that affect sequencing quality and to establish methods for measuring and reporting clinical sequencing standards. In this report, we propose a definition for reduced coverage regions and describe a set of standards for variant calling in clinical sequencing applications.MethodsTo enable sequencing centers to assess the regions of poor sequencing quality in their own data, we optimized and used a tool (ExCID) to identify reduced coverage loci within genes or regions of particular interest. We used this framework to examine sequencing data from 500 patients generated in 10 projects at sequencing centers in the National Human Genome Research Institute/National Cancer Institute Clinical Sequencing Exploratory Research Consortium.ResultsThis approach identified reduced coverage regions in clinically relevant genes, including known clinically relevant loci that were uniquely missed at individual centers, in multiple centers, and in all centers.ConclusionThis report provides a process road map for clinical sequencing centers looking to perform similar analyses on their data.

Utility and limitations of exome sequencing as a genetic diagnostic tool for children with hearing loss

PurposeHearing loss (HL) is the most common sensory disorder in children. Prompt molecular diagnosis may guide screening and management, especially in syndromic cases when HL is the single presenting feature. Exome sequencing (ES) is an appealing diagnostic tool for HL as the genetic causes are highly heterogeneous.MethodsES was performed on a prospective cohort of 43 probands with HL. Sequence data were analyzed for primary and secondary findings. Capture and coverage analysis was performed for genes and variants associated with HL.ResultsThe diagnostic rate using ES was 37.2%, compared with 15.8% for the clinical HL panel. Secondary findings were discovered in three patients. For 247 genes associated with HL, 94.7% of the exons were targeted for capture and 81.7% of these exons were covered at 20× or greater. Further analysis of 454 randomly selected HL-associated variants showed that 89% were targeted for capture and 75% were covered at a read depth of at least 20×.ConclusionES has an improved yield compared with clinical testing and may capture diagnoses not initially considered due to subtle clinical phenotypes. Technical challenges were identified, including inadequate capture and coverage of HL genes. Additional considerations of ES include secondary findings, cost, and turnaround time.