@article {372, title = {Homozygous and hemizygous CNV detection from exome sequencing data in a Mendelian disease cohort.}, journal = {Nucleic Acids Res}, volume = {45}, year = {2017}, month = {2017 02 28}, pages = {1633-1648}, abstract = {

We developed an algorithm, HMZDelFinder, that uses whole exome sequencing (WES) data to identify rare and intragenic homozygous and hemizygous (HMZ) deletions that may represent complete loss-of-function of the indicated gene. HMZDelFinder was applied to 4866 samples in the Baylor-Hopkins Center for Mendelian Genomics (BHCMG) cohort and detected 773 HMZ deletion calls (567 homozygous or 206 hemizygous) with an estimated sensitivity of 86.5\% (82\% for single-exonic and 88\% for multi-exonic calls) and precision of 78\% (53\% single-exonic and 96\% for multi-exonic calls). Out of 773 HMZDelFinder-detected deletion calls, 82 were subjected to array comparative genomic hybridization (aCGH) and/or breakpoint PCR and 64 were confirmed. These include 18 single-exon deletions out of which 8 were exclusively detected by HMZDelFinder and not by any of seven other CNV detection tools examined. Further investigation of the 64 validated deletion calls revealed at least 15 pathogenic HMZ deletions. Of those, 7 accounted for 17-50\% of pathogenic CNVs in different disease cohorts where 7.1-11\% of the molecular diagnosis solved rate was attributed to CNVs. In summary, we present an algorithm to detect rare, intragenic, single-exon deletion CNVs using WES data; this tool can be useful for disease gene discovery efforts and clinical WES analyses.

}, keywords = {Algorithms, Alternative Splicing, Cohort Studies, Computational Biology, Consanguinity, Datasets as Topic, DNA Copy Number Variations, Exome, Genetic Diseases, Inborn, Hemizygote, High-Throughput Nucleotide Sequencing, Homozygote, Humans, Inheritance Patterns, Models, Genetic, Pedigree, Reproducibility of Results, Sequence Deletion, Workflow}, issn = {1362-4962}, doi = {10.1093/nar/gkw1237}, author = {Gambin, Tomasz and Akdemir, Zeynep C and Yuan, Bo and Gu, Shen and Chiang, Theodore and Carvalho, Claudia M B and Shaw, Chad and Jhangiani, Shalini and Boone, Philip M and Eldomery, Mohammad K and Karaca, Ender and Bayram, Yavuz and Stray-Pedersen, Asbj{\o}rg and Muzny, Donna and Charng, Wu-Lin and Bahrambeigi, Vahid and Belmont, John W and Boerwinkle, Eric and Beaudet, Arthur L and Gibbs, Richard A and Lupski, James R} } @article {535, title = {Localized structural frustration for evaluating the impact of sequence variants.}, journal = {Nucleic Acids Res}, volume = {44}, year = {2016}, month = {2016 12 01}, pages = {10062-10073}, abstract = {

Population-scale sequencing is increasingly uncovering large numbers of rare single-nucleotide variants (SNVs) in coding regions of the genome. The rarity of these variants makes it challenging to evaluate their deleteriousness with conventional phenotype-genotype associations. Protein structures provide a way of addressing this challenge. Previous efforts have focused on globally quantifying the impact of SNVs on protein stability. However, local perturbations may severely impact protein functionality without strongly disrupting global stability (e.g. in relation to catalysis or allostery). Here, we describe a workflow in which localized frustration, quantifying unfavorable local interactions, is employed as a metric to investigate such effects. Using this workflow on the Protein Databank, we find that frustration produces many immediately intuitive results: for instance, disease-related SNVs create stronger changes in localized frustration than non-disease related variants, and rare SNVs tend to disrupt local interactions to a larger extent than common variants. Less obviously, we observe that somatic SNVs associated with oncogenes and tumor suppressor genes (TSGs) induce very different changes in frustration. In particular, those associated with TSGs change the frustration more in the core than the surface (by introducing loss-of-function events), whereas those associated with oncogenes manifest the opposite pattern, creating gain-of-function events.

}, keywords = {Computational Biology, Databases, Protein, Evolution, Molecular, Genes, Tumor Suppressor, Genetic Variation, Humans, Neoplasms, Oncogenes, Polymorphism, Single Nucleotide, Proteins, Workflow}, issn = {1362-4962}, doi = {10.1093/nar/gkw927}, author = {Kumar, Sushant and Clarke, Declan and Gerstein, Mark} } @article {176, title = {PacBio-LITS: a large-insert targeted sequencing method for characterization of human disease-associated chromosomal structural variations.}, journal = {BMC Genomics}, volume = {16}, year = {2015}, month = {2015 Mar 19}, pages = {214}, abstract = {

BACKGROUND: Generation of long (>5 Kb) DNA sequencing reads provides an approach for interrogation of complex regions in the human genome. Currently, large-insert whole genome sequencing (WGS) technologies from Pacific Biosciences (PacBio) enable analysis of chromosomal structural variations (SVs), but the cost to achieve the required sequence coverage across the entire human genome is high.

RESULTS: We developed a method (termed PacBio-LITS) that combines oligonucleotide-based DNA target-capture enrichment technologies with PacBio large-insert library preparation to facilitate SV studies at specific chromosomal regions. PacBio-LITS provides deep sequence coverage at the specified sites at substantially reduced cost compared with PacBio WGS. The efficacy of PacBio-LITS is illustrated by delineating the breakpoint junctions of low copy repeat (LCR)-associated complex structural rearrangements on chr17p11.2 in patients diagnosed with Potocki-Lupski syndrome (PTLS; MIM$\#$610883). We successfully identified previously determined breakpoint junctions in three PTLS cases, and also were able to discover novel junctions in repetitive sequences, including LCR-mediated breakpoints. The new information has enabled us to propose mechanisms for formation of these structural variants.

CONCLUSIONS: The new method leverages the cost efficiency of targeted capture-sequencing as well as the mappability and scaffolding capabilities of long sequencing reads generated by the PacBio platform. It is therefore suitable for studying complex SVs, especially those involving LCRs, inversions, and the generation of chimeric Alu elements at the breakpoints. Other genomic research applications, such as haplotype phasing and small insertion and deletion validation could also benefit from this technology.

}, keywords = {Chromosome Aberrations, Gene Library, Gene Rearrangement, Genetic Association Studies, Genomics, High-Throughput Nucleotide Sequencing, Humans, Workflow}, issn = {1471-2164}, doi = {10.1186/s12864-015-1370-2}, author = {Wang, Min and Beck, Christine R and English, Adam C and Meng, Qingchang and Buhay, Christian and Han, Yi and Doddapaneni, Harsha V and Yu, Fuli and Boerwinkle, Eric and Lupski, James R and Muzny, Donna M and Gibbs, Richard A} }