Correction: Sherloc: a comprehensive refinement of the ACMG-AMP variant classification criteria.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Sherloc: a comprehensive refinement of the ACMG-AMP variant classification criteria.

PurposeThe 2015 American College of Medical Genetics and Genomics-Association for Molecular Pathology (ACMG-AMP) guidelines were a major step toward establishing a common framework for variant classification. In practice, however, several aspects of the guidelines lack specificity, are subject to varied interpretations, or fail to capture relevant aspects of clinical molecular genetics. A simple implementation of the guidelines in their current form is insufficient for consistent and comprehensive variant classification.MethodsWe undertook an iterative process of refining the ACMG-AMP guidelines. We used the guidelines to classify more than 40,000 clinically observed variants, assessed the outcome, and refined the classification criteria to capture exceptions and edge cases. During this process, the criteria evolved through eight major and minor revisions.ResultsOur implementation: (i) separated ambiguous ACMG-AMP criteria into a set of discrete but related rules with refined weights; (ii) grouped certain criteria to protect against the overcounting of conceptually related evidence; and (iii) replaced the "clinical criteria" style of the guidelines with additive, semiquantitative criteria.ConclusionSherloc builds on the strong framework of 33 rules established by the ACMG-AMP guidelines and introduces 108 detailed refinements, which support a more consistent and transparent approach to variant classification.

Fine-scale recombination patterns differ between chimpanzees and humans.

Recombination rates seem to vary extensively along the human genome. Pedigree analysis suggests that rates vary by an order of magnitude when measured at the megabase scale, and at a finer scale, sperm typing studies point to the existence of recombination hotspots. These are short regions (1-2 kb) in which recombination rates are 10-1,000 times higher than the background rate. Less is known about how recombination rates change over time. Here we determined to what degree recombination rates are conserved among closely related species by estimating recombination rates from 14 Mb of linkage disequilibrium data in central chimpanzee and human populations. The results suggest that recombination hotspots are not conserved between the two species and that recombination rates in larger (50 kb) genomic regions are only weakly conserved. Therefore, the recombination landscape has changed markedly between the two species.

Possible precision medicine implications from genetic testing using combined detection of sequence and intragenic copy number variants in a large cohort with childhood epilepsy.

OBJECTIVE: Molecular genetic etiologies in epilepsy have become better understood in recent years, creating important opportunities for precision medicine. Building on these advances, detailed studies of the complexities and outcomes of genetic testing for epilepsy can provide useful insights that inform and refine diagnostic approaches and illuminate the potential for precision medicine in epilepsy. METHODS: We used a multi-gene next-generation sequencing (NGS) panel with simultaneous sequence and exonic copy number variant detection to investigate up to 183 epilepsy-related genes in 9769 individuals. Clinical variant interpretation was performed using a semi-quantitative scoring system based on existing professional practice guidelines. RESULTS: Molecular genetic testing provided a diagnosis in 14.9%-24.4% of individuals with epilepsy, depending on the NGS panel used. More than half of these diagnoses were in children younger than 5 years. Notably, the testing had possible precision medicine implications in 33% of individuals who received definitive diagnostic results. Only 30 genes provided 80% of molecular diagnoses. While most clinically significant findings were single-nucleotide variants, ~15% were other types that are often challenging to detect with traditional methods. In addition to clinically significant variants, there were many others that initially had uncertain significance; reclassification of 1612 such variants with parental testing or other evidence contributed to 18.5% of diagnostic results overall and 6.1% of results with precision medicine implications. SIGNIFICANCE: Using an NGS gene panel with key high-yield genes and robust analytic sensitivity as a first-tier test early in the diagnostic process, especially for children younger than 5 years, can possibly enable precision medicine approaches in a significant number of individuals with epilepsy.

DNA hybridization to mismatched templates: a chip study.

High-density oligonucleotide arrays are among the most rapidly expanding technologies in biology today. In the GeneChip system, the reconstruction of the sample mRNA concentrations depends upon the differential signal generated by hybridizing the RNA to two nearly identical templates: a perfect match probe (PM) containing the exact biological sequence; and a single mismatch (MM) differing from the PM by a single base substitution. It has been observed that a large fraction of MMs repetitively bind targets better than the PMs, against the obvious expectation of sequence specificity. We examine this problem via statistical analysis of a large set of microarray experiments. We classify the probes according to their signal to noise (S/N) ratio, defined as the eccentricity of a (PM,MM) pair's "trajectory" across many experiments. Of those probes having large S/N (>3) only a fraction behave consistently with the commonly assumed hybridization model. Our results imply that the physics of DNA hybridization in microarrays is more complex than expected, and suggest estimators for the target RNA concentration.

In vivo transcriptional profile analysis reveals RNA splicing and chromatin remodeling as prominent processes for adult neurogenesis.

Neural stem cells and neurogenesis persist in the adult mammalian brain subventricular zone (SVZ). Cells born in the rodent SVZ migrate to the olfactory bulb (Ob) where they differentiate into interneurons. To determine the gene expression and functional profile of SVZ neurogenesis, we performed three complementary sets of transcriptional analysis experiments using Affymetrix GeneChips: (1) comparison of adult mouse SVZ and Ob gene expression profiles with those of the striatum, cerebral cortex, and hippocampus; (2) profiling of SVZ stem cells and ependyma isolated by fluorescent-activated cell sorting (FACS); and (3) analysis of gene expression changes during in vivo SVZ regeneration after anti-mitotic treatment. Gene Ontology (GO) analysis of data from these three separate approaches showed that in adult SVZ neurogenesis, RNA splicing and chromatin remodeling are biological processes as statistically significant as cell proliferation, transcription, and neurogenesis. In non-neurogenic brain regions, RNA splicing and chromatin remodeling were not prominent processes. Fourteen mRNA splicing factors including Sf3b1, Sfrs2, Lsm4, and Khdrbs1/Sam68 were detected along with 9 chromatin remodeling genes including Mll, Bmi1, Smarcad1, Baf53a, and Hat1. We validated the transcriptional profile data with Northern blot analysis and in situ hybridization. The data greatly expand the catalogue of cell cycle components, transcription factors, and migration genes for adult SVZ neurogenesis and reveal RNA splicing and chromatin remodeling as prominent biological processes for these germinal cells.

Empirical characterization of the expression ratio noise structure in high-density oligonucleotide arrays.

BACKGROUND: High-density oligonucleotide arrays (HDONAs) are a powerful tool for assessing differential mRNA expression levels. To establish the statistical significance of an observed change in expression, one must take into account the noise introduced by the enzymatic and hybridization steps, called type I noise. We undertake an empirical characterization of the experimental repeatability of results by carrying out statistical analysis of a large number of duplicate HDONA experiments. RESULTS: We assign scoring functions for expression ratios and associated quality measures. Both the perfect-match (PM) probes and the differentials between PM and single-mismatch (MM) probes are considered as raw intensities. We then calculate the log-ratio of the noise structure using robust estimates of their intensity-dependent variance. The noise structure in the log-ratios follows a local log-normal distribution in both the PM and PM-MM cases. Significance relative to the type I noise can therefore be quantified reliably using the local standard deviation (SD). We discuss the intensity dependence of the SD and show that ratio scores greater than 1.25 are significant in the mid- to high-intensity range. CONCLUSIONS: The noise inherent in HDONAs is characteristically dependent on intensity and can be well described in terms of local normalization of log-ratio distributions. Therefore, robust estimates of the local SD of these distributions provide a simple and powerful way to assess significance (relative to type I noise) in differential gene expression, and will be helpful in practice for improving the reliability of predictions from hybridization experiments.

Characterization of the expression ratio noise structure in high-density oligonucleotide arrays.

BACKGROUND: High-density oligonucleotide microarrays provide a powerful tool for assessing differential mRNA expression levels. Characterizing the noise resulting from the enzymatic and hybridization steps, called type I noise, is essential for attributing significance measures to the differential expression scores. We introduce scoring functions for expression ratios, and associated quality measures. Both the PM (Perfect Match) probes and PM-MM differentials (MM is the single MisMatch) are considered as raw intensities. We then characterize the log-ratio noise structure using robust estimates of their intensity dependent variance. RESULTS: We show the relationships between the obtained ratios and their quality measures. The complementarity of PM and PM-MM methods is emphasized by the probe sets signal to noise measures. Using a large set of replicate experiments, we demonstrate that the noise structure in the log-ratios very closely follows a local log-normal distribution for both the PM and PM-MM cases. Therefore, significance relative to the type I noise can be quantified reliably using the local STD. We discuss the intensity dependence of the STD and show that ratio scores >1.25 are significant in the mid- to high-intensity range. CONCLUSIONS: The ratio noise structure inherent to high-density oligonucleotide arrays can be well described in terms of local log-normal ratio distributions with characteristic intensity dependence. Therefore, robust estimates of the local STD of these distributions provide a simple and powerful way for assessing significance (relative to type I noise) in differential gene expression. This approach will be helpful for improving the reliability of predictions from hybridization experiments in general.

Genomic DNA insertions and deletions occur frequently between humans and nonhuman primates.

Comparative DNA sequence studies between humans and nonhuman primates will be important for understanding the genetic basis of the phenotypic differences between these species. Here we compare approximately 27 Mb of human chromosome 21 with chimpanzee DNA sequences identifying 57 genomic rearrangements (deletions and insertions ranging in size from 0.2 to 8.0 kb) between the two species. These rearrangements are distributed along the entire length of chromosome 21, with approximately 35% found in genomic intervals encoding genes (genic intervals), and have occurred in the genomes of both humans and chimpanzees. Comparison of approximately 9 Mb of human chromosome 21 with orangutan, rhesus macaque, and woolly monkey DNA sequences identified a combined total of 114 genomic rearrangements between humans and nonhuman primates. Analysis of these rearrangements revealed that they are randomly distributed with respect to genic and nongenic intervals and identified one deletion that has likely resulted in the inactivation of a gene (beta1,3-galactosyltransferase) in the woolly monkey. Our data show that genomic rearrangements have occurred frequently during primate genome evolution and significantly contribute to the DNA differences between these species. These DNA rearrangements are commonly found in genic intervals, and thus provide natural starting points for focused investigations of qualitative and quantitative gene expression differences between humans and other primates.

Matching strategies for genetic association studies in structured populations.

Association studies in populations that are genetically heterogeneous can yield large numbers of spurious associations if population subgroups are unequally represented among cases and controls. This problem is particularly acute for studies involving pooled genotyping of very large numbers of single-nucleotide-polymorphism (SNP) markers, because most methods for analysis of association in structured populations require individual genotyping data. In this study, we present several strategies for matching case and control pools to have similar genetic compositions, based on ancestry information inferred from genotype data for approximately 300 SNPs tiled on an oligonucleotide-based genotyping array. We also discuss methods for measuring the impact of population stratification on an association study. Results for an admixed population and a phenotype strongly confounded with ancestry show that these simple matching strategies can effectively mitigate the impact of population stratification.

Segmental phylogenetic relationships of inbred mouse strains revealed by fine-scale analysis of sequence variation across 4.6 mb of mouse genome.

High-density SNP screening of panels of inbred mouse strains has been proposed as a method to accelerate the identification of genes associated with complex biomedical phenotypes. To evaluate the potential of these studies, a more detailed understanding of the fine structure of sequence variation across inbred mouse strains is needed. Here, we use high-density oligonucleotide arrays to discover an extremely dense set of SNPs in 13 classical and two wild-derived inbred strains in five genomic intervals totaling 4.6 Mb of DNA sequence, and then analyze the segmental haplotype structure defined by these high-density SNPs. This analysis reveals segments ranging from 12 to 608 kb in length within which the inbred strains have a simple and distinct phylogenetic relationship with typically two or three clades accounting for the 13 classical strains examined. The phylogenetic relationships among strains change abruptly and unpredictably from segment to segment, and are distinct in each of the five genomic regions examined. The data suggest that at least 12 strains would need to be resequenced for exhaustive SNP discovery in every region of the mouse genome, that approximately 97% of the variation among inbred strains is ancestral (between clades) and approximately 3% private (within clades), and provides critical insights into the proposed use of panels of inbred strains to identify genes underlying quantitative trait loci.