Gholson J. Lyon

Low concordance of multiple variant-calling pipelines: practical implications for exome and genome sequencing

BackgroundTo facilitate the clinical implementation of genomic medicine by next-generation sequencing, it will be critically important to obtain accurate and consistent variant calls on personal genomes. Multiple software tools for variant calling are available, but it is unclear how comparable these tools are or what their relative merits in real-world scenarios might be.MethodsWe sequenced 15 exomes from four families using commercial kits (Illumina HiSeq 2000 platform and Agilent SureSelect version 2 capture kit), with approximately 120X mean coverage. We analyzed the raw data using near-default parameters with five different alignment and variant-calling pipelines (SOAP, BWA-GATK, BWA-SNVer, GNUMAP, and BWA-SAMtools). We additionally sequenced a single whole genome using the sequencing and analysis pipeline from Complete Genomics (CG), with 95% of the exome region being covered by 20 or more reads per base. Finally, we validated 919 single-nucleotide variations (SNVs) and 841 insertions and deletions (indels), including similar fractions of GATK-only, SOAP-only, and shared calls, on the MiSeq platform by amplicon sequencing with approximately 5000X mean coverage.ResultsSNV concordance between five Illumina pipelines across all 15 exomes was 57.4%, while 0.5 to 5.1% of variants were called as unique to each pipeline. Indel concordance was only 26.8% between three indel-calling pipelines, even after left-normalizing and intervalizing genomic coordinates by 20 base pairs. There were 11% of CG variants falling within targeted regions in exome sequencing that were not called by any of the Illumina-based exome analysis pipelines. Based on targeted amplicon sequencing on the MiSeq platform, 97.1%, 60.2%, and 99.1% of the GATK-only, SOAP-only and shared SNVs could be validated, but only 54.0%, 44.6%, and 78.1% of the GATK-only, SOAP-only and shared indels could be validated. Additionally, our analysis of two families (one with four individuals and the other with seven), demonstrated additional accuracy gained in variant discovery by having access to genetic data from a multi-generational family.ConclusionsOur results suggest that more caution should be exercised in genomic medicine settings when analyzing individual genomes, including interpreting positive and negative findings with scrutiny, especially for indels. We advocate for renewed collection and sequencing of multi-generational families to increase the overall accuracy of whole genomes.

Exfoliatin-Producing Strains Define a Fourth agr Specificity Group in Staphylococcus aureus

The staphylococcal virulon is activated by the density-sensing agr system, which is autoinduced by a short peptide (autoinducing peptide [AIP]) processed from a propeptide encoded by agrD. A central segment of the agr locus, consisting of the C-terminal two-thirds of AgrB (the putative processing enzyme), AgrD, and the N-terminal half of AgrC (the receptor), shows striking interstrain variation. This finding has led to the division of Staphylococcus aureus isolates into three different agr specificity groups and to the division of non-aureus staphylococci into a number of others. The AIPs cross-inhibit the agr responses between groups. We have previously shown that most menstrual toxic shock strains belong to agr specificity group III but that no strong clinical identity has been associated with strains of the other two groups. In the present report, we demonstrate a fourth agr specificity group among S. aureus strains and show that most exfoliatin-producing strains belong to this group. A striking common feature of group IV strains is activation of the agr response early in exponential phase, at least 2 h earlier than in strains of the other groups. This finding raises the question of the biological significance of the agr autoinduction threshold.

Mice deficient for the secreted glycoprotein SPARC/osteonectin/BM40 develop normally but show severe age-onset cataract formation and disruption of the lens.

SPARC (secreted protein acidic and rich in cysteine, also known as osteonectin/BM40) is a secreted Ca2+‐binding glycoprotein that interacts with a range of extracellular matrix molecules, including collagen IV. It is widely expressed during embryogenesis, and in vitro studies have suggested roles in the regulation of cell adhesion and proliferation, and in the modulation of cytokine activity. In order to analyse the function of this protein in vivo, the endogenous Sparc locus was disrupted by homologous recombination in murine embryonic stem cells. SPARC‐deficient mice (Sparctm1Cam) appear normal and fertile until around 6 months of age, when they develop severe eye pathology characterized by cataract formation and rupture of the lens capsule. The first sign of lens pathology occurs in the equatorial bow region where vacuoles gradually form within differentiating epithelial cells and fibre cells. The lens capsule, however, shows no qualitative changes in the major basal lamina proteins laminin, collagen IV, perlecan or entactin. These mice are an excellent resource for further studies on how SPARC affects cell behaviour in vivo.

Genome-wide copy number variation study associates metabotropic glutamate receptor gene networks with attention deficit hyperactivity disorder

Attention deficit hyperactivity disorder (ADHD) is a common, heritable neuropsychiatric disorder of unknown etiology. We performed a whole-genome copy number variation (CNV) study on 1,013 cases with ADHD and 4,105 healthy children of European ancestry using 550,000 SNPs. We evaluated statistically significant findings in multiple independent cohorts, with a total of 2,493 cases with ADHD and 9,222 controls of European ancestry, using matched platforms. CNVs affecting metabotropic glutamate receptor genes were enriched across all cohorts (P = 2.1 × 10−9). We saw GRM5 (encoding glutamate receptor, metabotropic 5) deletions in ten cases and one control (P = 1.36 × 10−6). We saw GRM7 deletions in six cases, and we saw GRM8 deletions in eight cases and no controls. GRM1 was duplicated in eight cases. We experimentally validated the observed variants using quantitative RT-PCR. A gene network analysis showed that genes interacting with the genes in the GRM family are enriched for CNVs in ∼10% of the cases (P = 4.38 × 10−10) after correction for occurrence in the controls. We identified rare recurrent CNVs affecting glutamatergic neurotransmission genes that were overrepresented in multiple ADHD cohorts.

The Human Phenotype Ontology in 2017

Deep phenotyping has been defined as the precise and comprehensive analysis of phenotypic abnormalities in which the individual components of the phenotype are observed and described. The three components of the Human Phenotype Ontology (HPO; www.human-phenotype-ontology.org) project are the phenotype vocabulary, disease-phenotype annotations and the algorithms that operate on these. These components are being used for computational deep phenotyping and precision medicine as well as integration of clinical data into translational research. The HPO is being increasingly adopted as a standard for phenotypic abnormalities by diverse groups such as international rare disease organizations, registries, clinical labs, biomedical resources, and clinical software tools and will thereby contribute toward nascent efforts at global data exchange for identifying disease etiologies. This update article reviews the progress of the HPO project since the debut Nucleic Acids Research database article in 2014, including specific areas of expansion such as common (complex) disease, new algorithms for phenotype driven genomic discovery and diagnostics, integration of cross-species mapping efforts with the Mammalian Phenotype Ontology, an improved quality control pipeline, and the addition of patient-friendly terminology.

Using VAAST to identify an X-linked disorder resulting in lethality in male infants due to N-terminal acetyltransferase deficiency.

We have identified two families with a previously undescribed lethal X-linked disorder of infancy; the disorder comprises a distinct combination of an aged appearance, craniofacial anomalies, hypotonia, global developmental delays, cryptorchidism, and cardiac arrhythmias. Using X chromosome exon sequencing and a recently developed probabilistic algorithm aimed at discovering disease-causing variants, we identified in one family a c.109T>C (p.Ser37Pro) variant in NAA10, a gene encoding the catalytic subunit of the major human N-terminal acetyltransferase (NAT). A parallel effort on a second unrelated family converged on the same variant. The absence of this variant in controls, the amino acid conservation of this region of the protein, the predicted disruptive change, and the co-occurrence in two unrelated families with the same rare disorder suggest that this is the pathogenic mutation. We confirmed this by demonstrating a significantly impaired biochemical activity of the mutant hNaa10p, and from this we conclude that a reduction in acetylation by hNaa10p causes this disease. Here we provide evidence of a human genetic disorder resulting from direct impairment of N-terminal acetylation, one of the most common protein modifications in humans.

SNVer: a statistical tool for variant calling in analysis of pooled or individual next-generation sequencing data

We develop a statistical tool SNVer for calling common and rare variants in analysis of pooled or individual next-generation sequencing (NGS) data. We formulate variant calling as a hypothesis testing problem and employ a binomial–binomial model to test the significance of observed allele frequency against sequencing error. SNVer reports one single overall P-value for evaluating the significance of a candidate locus being a variant based on which multiplicity control can be obtained. This is particularly desirable because tens of thousands loci are simultaneously examined in typical NGS experiments. Each user can choose the false-positive error rate threshold he or she considers appropriate, instead of just the dichotomous decisions of whether to ‘accept or reject the candidates’ provided by most existing methods. We use both simulated data and real data to demonstrate the superior performance of our program in comparison with existing methods. SNVer runs very fast and can complete testing 300 K loci within an hour. This excellent scalability makes it feasible for analysis of whole-exome sequencing data, or even whole-genome sequencing data using high performance computing cluster. SNVer is freely available at http://snver.sourceforge.net/.

Lifetime Prevalence, Age of Risk, and Genetic Relationships of Comorbid Psychiatric Disorders in Tourette Syndrome

IMPORTANCE Tourette syndrome (TS) is characterized by high rates of psychiatric comorbidity; however, few studies have fully characterized these comorbidities. Furthermore, most studies have included relatively few participants (<200), and none has examined the ages of highest risk for each TS-associated comorbidity or their etiologic relationship to TS. OBJECTIVE To characterize the lifetime prevalence, clinical associations, ages of highest risk, and etiology of psychiatric comorbidity among individuals with TS. DESIGN, SETTING, AND PARTICIPANTS Cross-sectional structured diagnostic interviews conducted between April 1, 1992, and December 31, 2008, of participants with TS (n = 1374) and TS-unaffected family members (n = 1142). MAIN OUTCOMES AND MEASURES Lifetime prevalence of comorbid DSM-IV-TR disorders, their heritabilities, ages of maximal risk, and associations with symptom severity, age at onset, and parental psychiatric history. RESULTS The lifetime prevalence of any psychiatric comorbidity among individuals with TS was 85.7%; 57.7% of the population had 2 or more psychiatric disorders. The mean (SD) number of lifetime comorbid diagnoses was 2.1 (1.6); the mean number was 0.9 (1.3) when obsessive-compulsive disorder (OCD) and attention-deficit/hyperactivity disorder (ADHD) were excluded, and 72.1% of the individuals met the criteria for OCD or ADHD. Other disorders, including mood, anxiety, and disruptive behavior, each occurred in approximately 30% of the participants. The age of greatest risk for the onset of most comorbid psychiatric disorders was between 4 and 10 years, with the exception of eating and substance use disorders, which began in adolescence (interquartile range, 15-19 years for both). Tourette syndrome was associated with increased risk of anxiety (odds ratio [OR], 1.4; 95% CI, 1.0-1.9; P = .04) and decreased risk of substance use disorders (OR, 0.6; 95% CI, 0.3-0.9; P = .02) independent from comorbid OCD and ADHD; however, high rates of mood disorders among participants with TS (29.8%) may be accounted for by comorbid OCD (OR, 3.7; 95% CI, 2.9-4.8; P < .001). Parental history of ADHD was associated with a higher burden of non-OCD, non-ADHD comorbid psychiatric disorders (OR, 1.86; 95% CI, 1.32-2.61; P < .001). Genetic correlations between TS and mood (RhoG, 0.47), anxiety (RhoG, 0.35), and disruptive behavior disorders (RhoG, 0.48), may be accounted for by ADHD and, for mood disorders, by OCD. CONCLUSIONS AND RELEVANCE This study is, to our knowledge, the most comprehensive of its kind. It confirms the belief that psychiatric comorbidities are common among individuals with TS, demonstrates that most comorbidities begin early in life, and indicates that certain comorbidities may be mediated by the presence of comorbid OCD or ADHD. In addition, genetic analyses suggest that some comorbidities may be more biologically related to OCD and/or ADHD rather than to TS.

Genome-wide association study of Tourette's syndrome

Tourettes syndrome (TS) is a developmental disorder that has one of the highest familial recurrence rates among neuropsychiatric diseases with complex inheritance. However, the identification of definitive TS susceptibility genes remains elusive. Here, we report the first genome-wide association study (GWAS) of TS in 1285 cases and 4964 ancestry-matched controls of European ancestry, including two European-derived population isolates, Ashkenazi Jews from North America and Israel and French Canadians from Quebec, Canada. In a primary meta-analysis of GWAS data from these European ancestry samples, no markers achieved a genome-wide threshold of significance (P<5 × 10−8); the top signal was found in rs7868992 on chromosome 9q32 within COL27A1 (P=1.85 × 10−6). A secondary analysis including an additional 211 cases and 285 controls from two closely related Latin American population isolates from the Central Valley of Costa Rica and Antioquia, Colombia also identified rs7868992 as the top signal (P=3.6 × 10−7 for the combined sample of 1496 cases and 5249 controls following imputation with 1000 Genomes data). This study lays the groundwork for the eventual identification of common TS susceptibility variants in larger cohorts and helps to provide a more complete understanding of the full genetic architecture of this disorder.

Detection of secreted peptides by using hypothesis-driven multistage mass spectrometry

A method is presented for the rapid detection and characterization of trace amounts of peptides secreted from microorganisms, including pheromones, virulence factors, and quorum-sensing peptides. The procedure, based on targeted multistage MS, uses a novel matrix-assisted laser desorption/ionization-ion trap mass spectrometer to overcome limitations of current MS methods (limited dynamic range, signal suppression effects, and chemical noise) that impair observation of low abundance peptides from complex biological matrixes. Here, secreted peptides that are hypothesized to be present in the supernatant, but that may not be sufficiently abundant to be observed in single-stage mass spectra, are subjected to multistage MS. Highly specific fragmentation signatures enable unambiguous identification of the peptides of interest and differentiation of the signals from the background. As examples, we demonstrate the rapid (<1 min) determination of the mating type of cells in colonies of Saccharomyces cerevisiae and the elucidation of autoinducing peptides (AIPs) from supernatants of pathogenic Staphylococci. We confirm the primary structures of the agrD encoded cyclic AIPs of Staphylococcus aureus for groups I, II, and IV and provide direct evidence that the native group-III AIP is a heptapeptide (INCDFLL). We also show that the homologous peptide from Staphylococcus intermedius is a nonapeptide (RIPTSTGFF) with a lactone ring formed through condensation of the serine side chain with the C terminus of the peptide. This is the first demonstration of cyclization in a staphylococcal AIP that occurs via lactone formation. These examples demonstrate the analytical power of the present procedure for characterizing secreted peptides and its potential utility for identifying microorganisms.