Ocular manifestations in a cohort of 43 patients with KBG syndrome

Ophthalmological conditions are underreported in patients with KBG syndrome, which is classically described as presenting with dental, developmental, intellectual, skeletal, and craniofacial abnormalities. This study analyzed the prevalence of four ophthalmological conditions (strabismus, astigmatism, myopia, hyperopia) in 43 patients with KBG syndrome carrying variants in ANKRD11 or deletions in 16q24.3 and compared it to the literature. Forty-three patients were recruited via self-referral or a private Facebook group hosted by the KBG Foundation, with 40 of them having pathogenic or likely pathogenic variants. Virtual interviews were conducted to collect a comprehensive medical history verified by medical records. From these records, data analysis was performed to calculate the prevalence of ophthalmological conditions. Out of the 40 participants with pathogenic or likely pathogenic variants, strabismus was reported in 9 (22.5%) participants, while astigmatism, myopia, and hyperopia were reported in 11 (27.5%), 6 (15.0%), and 8 (20.0%) participants, respectively. Other reported conditions include anisometropia, amblyopia, and nystagmus. When compared to the literature, the prevalence of strabismus and refractive errors is higher than other studies. However, more research is needed to determine if variants in ANKRD11 play a role in abnormal development of the visual system. In patients with established KBG syndrome, screening for misalignment or refractive errors should be done, as interventions in patients with these conditions can improve functioning and quality of life.

GestaltMatcher Database - A global reference for the facial phenotypic variability of rare human diseases

Dysmorphologists sometimes encounter challenges in recognizing disorders due to phenotypic variability influenced by factors such as age and ethnicity. Moreover, the performance of Next Generation Phenotyping Tools such as GestaltMatcher is dependent on the diversity of the training set. Therefore, we developed GestaltMatcher Database (GMDB) - a global reference for the phenotypic variability of rare diseases that complies with the FAIR-principles. We curated dysmorphic patient images and metadata from 2,224 publications, transforming GMDB into an online dynamic case report journal. To encourage clinicians worldwide to contribute, each case can receive a Digital Object Identifier (DOI), making it a citable micro-publication. This resulted in a collection of 2,312 unpublished images, partly with longitudinal data. We have compiled a collection of 10,189 frontal images from 7,695 patients representing 683 disorders. The web interface enables gene- and phenotype-centered queries for registered users (https://db.gestaltmatcher.org/). Despite the predominant European ancestry of most patients (59%), our global collaborations have facilitated the inclusion of data from frequently underrepresented ethnicities, with 17% Asian, 4% African, and 6% with other ethnic backgrounds. The analysis has revealed a significant enhancement in GestaltMatcher performance across all ethnic groups, incorporating non-European ethnicities, showcasing a remarkable increase in Top-1-Accuracy by 31.56% and Top-5-Accuracy by 12.64%. Importantly, this improvement was achieved without altering the performance metrics for European patients. GMDB addresses dysmorphology challenges by representing phenotypic variability and including underrepresented groups, enhancing global diagnostic rates and serving as a vital clinician reference database.

The Challenge of Examining Social Determinants of Health in People Living With Tourette Syndrome

BACKGROUND

To examine the association between race, ethnicity, and parental educational attainment on tic-related outcomes among Tourette Syndrome (TS) participants in the Tourette Association of America International Consortium for Genetics (TAAICG) database.

METHODS

723 participants in the TAAICG dataset aged ≤21 years were included. The relationships between tic-related outcomes and race and ethnicity were examined using linear and logistic regressions. Parametric and nonparametric tests were performed to examine the association between parental educational attainment and tic-related outcomes.

RESULTS

Race and ethnicity were collapsed as non-Hispanic white (N=566, 88.0%) versus Other (N=77, 12.0%). Tic symptom onset was earlier by 1.1 years (P < 0.0001) and TS diagnosis age was earlier by 0.9 years (P = 0.0045) in the Other group (versus non-Hispanic white). Sex and parental education as covariates did not contribute to the differences observed in TS diagnosis age. There were no significant group differences observed across the tic-related outcomes in parental education variable.

CONCLUSIONS

Our study was limited by the low number of nonwhite or Hispanic individuals in the cohort. Racial and ethnic minoritized groups experienced an earlier age of TS diagnosis than non-Hispanic white individuals. Tic severity did not differ between the two groups, and parental educational attainment did not affect tic-related outcomes. There remain significant disparities and gaps in knowledge regarding TS and associated comorbid conditions. Our study suggests the need for more proactive steps to engage individuals with tic disorders from all racial and ethnic minoritized groups to participate in research studies.

Longitudinal Adaptive Behavioral Outcomes in Ogden Syndrome by Seizure Status and Therapeutic Intervention

Ogden syndrome, also known as NAA10-related neurodevelopmental syndrome, is a rare genetic condition associated with pathogenic variants in the NAA10 N-terminal acetylation family of proteins. The condition was initially described in 2011, and is characterized by a range of neurologic symptoms, including intellectual disability and seizures, as well as developmental delays, psychiatric symptoms, congenital heart abnormalities, hypotonia and others. Previously published articles have described the etiology and phenotype of Ogden syndrome, mostly with retrospective analyses; herein, we report prospective data concerning its progress over time. Additionally, we describe the nature of seizures in this condition in greater detail, as well as investigate how already-available non-pharmaceutical therapies impact individuals with NAA10-related neurodevelopmental syndrome. Using Vineland-3 scores, we show decline in cognitive function over time in individuals with Ogden syndrome. Sub-domain analysis found the decline to be present across all modalities. Additional investigation between seizure and non-seizure groups showed no significant difference in adaptive behavior outcomes. Therapy investigation showed speech therapy to be the most commonly used therapy by individuals with NAA10-related neurodevelopmental syndrome, followed by occupational and physical therapy. with more severely affected individuals receiving more types of therapy than their less-severe counterparts. Early intervention analysis was only significantly effective for speech therapy, with analyses of all other therapies being non-significant. Our study portrays the decline in cognitive function over time of individuals within our cohort, independent of seizure status and therapies being received, and highlights the urgent need for the development of effective treatments for Ogden syndrome.

Ophthalmic Manifestations of NAA10-Related and NAA15-Related Neurodevelopmental Syndrome: Analysis of Cortical Visual Impairment and Refractive Errors

NAA10-related and NAA15-related neurodevelopmental syndrome, otherwise known as Ogden Syndrome, is known to present with varying degrees of intellectual disability, hypotonia, congenital cardiac abnormalities, seizures, and delayed speech and motor development. However, the ophthalmic manifestations of NAA10 and NAA15 mutations are not yet fully characterized or understood. This study analyzed the prevalence of six ophthalmic conditions (cortical visual impairment, myopia, hyperopia, strabismus, nystagmus, and astigmatism) in 67 patients with pathogenic mutations in the NAA10 cohort (54 inherited, 10 de novo; 65 missense, 2 frameshift) and 19 patients with pathogenic mutations in the NAA15 cohort (18 de novo; 8 frameshift, 4 missense, 4 nonsense, and 1 splice site). Patients were interviewed virtually or in-person to collect a comprehensive medical history verified by medical records. These records were then analyzed to calculate the prevalence of these ophthalmic manifestations in each cohort. Analysis revealed a higher prevalence of ophthalmic conditions in our NAA10 cohort compared to existing literature (myopia 25.4% vs. 4.7%; astigmatism 37.3% vs. 13.2%; strabismus 28.4% vs. 3.8%; CVI 22.4% vs. 8.5%, respectively). No statistically significant differences were identified between the NAA10 and NAA15 mutations. Our study includes novel neuroimaging of 13 NAA10 and 5 NAA15 probands, which provides no clear correlation between globe size and severity of comorbid ophthalmic disease. Finally, anecdotal evidence was compiled to underscore the importance of early ophthalmologic evaluations and therapeutic interventions.

Evaluating possible maternal effect lethality and genetic background effects in Naa10 knockout mice

Amino-terminal (Nt-) acetylation (NTA) is a common protein modification, affecting approximately 80% of all human proteins. The human essential X-linked gene, NAA10, encodes for the enzyme NAA10, which is the catalytic subunit in the N-terminal acetyltransferase A (NatA) complex. There is extensive genetic variation in humans with missense, splice-site, and C-terminal frameshift variants in NAA10. In mice, Naa10 is not an essential gene, as there exists a paralogous gene, Naa12, that substantially rescues Naa10 knockout mice from embryonic lethality, whereas double knockouts (Naa10-/Y Naa12-/-) are embryonic lethal. However, the phenotypic variability in the mice is nonetheless quite extensive, including piebaldism, skeletal defects, small size, hydrocephaly, hydronephrosis, and neonatal lethality. Here we replicate these phenotypes with new genetic alleles in mice, but we demonstrate their modulation by genetic background and environmental effects. We cannot replicate a prior report of "maternal effect lethality" for heterozygous Naa10-/X female mice, but we do observe a small amount of embryonic lethality in the Naa10-/Y male mice on the inbred genetic background in this different animal facility.

PhenoScore quantifies phenotypic variation for rare genetic diseases by combining facial analysis with other clinical features using a machine-learning framework

Several molecular and phenotypic algorithms exist that establish genotype-phenotype correlations, including facial recognition tools. However, no unified framework that investigates both facial data and other phenotypic data directly from individuals exists. We developed PhenoScore: an open-source, artificial intelligence-based phenomics framework, combining facial recognition technology with Human Phenotype Ontology data analysis to quantify phenotypic similarity. Here we show PhenoScore's ability to recognize distinct phenotypic entities by establishing recognizable phenotypes for 37 of 40 investigated syndromes against clinical features observed in individuals with other neurodevelopmental disorders and show it is an improvement on existing approaches. PhenoScore provides predictions for individuals with variants of unknown significance and enables sophisticated genotype-phenotype studies by testing hypotheses on possible phenotypic (sub)groups. PhenoScore confirmed previously known phenotypic subgroups caused by variants in the same gene for SATB1, SETBP1 and DEAF1 and provides objective clinical evidence for two distinct ADNP-related phenotypes, already established functionally.

Documentation and prevalence of prenatal and neonatal outcomes in a cohort of individuals with KBG syndrome

Ankyrin Repeat Domain 11 (ANKRD11) gene mutations are associated with KBG syndrome, a developmental disability that affects multiple organ systems. The function of ANKRD11 in human growth and development is not clear, but gene knockout or mutation are lethal in mice embryos and/or pups. In addition, it plays a vital role in chromatin regulation and transcription. Individuals with KBG syndrome are often misdiagnosed or remain undiagnosed until later in life. This is largely due to KBG syndrome's varying and nonspecific phenotypes as well as a lack of accessible genetic testing and prenatal screening. This study documents perinatal outcomes for individuals with KBG syndrome. We obtained data from 42 individuals through videoconferences, medical records, and emails. 45.2% of our cohort was born by C-section, 33.3% had a congenital heart defect, 23.8% were born prematurely, 23.8% were admitted to the NICU, 14.3% were small for gestational age, and 14.3% of the families had a history of miscarriage. These rates were higher in our cohort compared to the overall population, including non-Hispanic and Hispanic populations. Other reports included feeding difficulties (21.4%), neonatal jaundice (14.3%), decreased fetal movement (7.1%), and pleural effusions in utero (4.7%). Comprehensive perinatal studies about KBG syndrome and updated documentation of its phenotypes are important in ensuring prompt diagnosis and can facilitate correct management.

Expanding the phenotypic spectrum of NAA10-related neurodevelopmental syndrome and NAA15-related neurodevelopmental syndrome

Amino-terminal (Nt-) acetylation (NTA) is a common protein modification, affecting 80% of cytosolic proteins in humans. The human essential gene, NAA10, encodes for the enzyme NAA10, which is the catalytic subunit in the N-terminal acetyltransferase A (NatA) complex, also including the accessory protein, NAA15. The full spectrum of human genetic variation in this pathway is currently unknown. Here we reveal the genetic landscape of variation in NAA10 and NAA15 in humans. Through a genotype-first approach, one clinician interviewed the parents of 56 individuals with NAA10 variants and 19 individuals with NAA15 variants, which were added to all known cases (N = 106 for NAA10 and N = 66 for NAA15). Although there is clinical overlap between the two syndromes, functional assessment demonstrates that the overall level of functioning for the probands with NAA10 variants is significantly lower than the probands with NAA15 variants. The phenotypic spectrum includes variable levels of intellectual disability, delayed milestones, autism spectrum disorder, craniofacial dysmorphology, cardiac anomalies, seizures, and visual abnormalities (including cortical visual impairment and microphthalmia). One female with the p.Arg83Cys variant and one female with an NAA15 frameshift variant both have microphthalmia. The frameshift variants located toward the C-terminal end of NAA10 have much less impact on overall functioning, whereas the females with the p.Arg83Cys missense in NAA10 have substantial impairment. The overall data are consistent with a phenotypic spectrum for these alleles, involving multiple organ systems, thus revealing the widespread effect of alterations of the NTA pathway in humans.

Phenotypic variability and gastrointestinal manifestations/interventions for growth in NAA10-related neurodevelopmental syndrome

Our study of 61 children with NAA10-related neurodevelopmental syndrome, an X-linked disorder due to NAA10 gene variants, demonstrated a high prevalence of growth failure, with weight and height percentiles often in the failure-to-thrive diagnostic range; however, dramatic weight fluctuations and phenotypic variability is evidenced in the growth parameters of this population. Although never previously explored in depth, the gastrointestinal pathology associated with NAA10-related neurodevelopmental syndrome includes feeding difficulties in infancy, dysphagia, GERD/silent reflux, vomiting, constipation, diarrhea, bowel incontinence, and presence of eosinophils on esophageal endoscopy, in order from most to least prevalent. Additionally, the gastrointestinal symptom profile for children with this syndrome has been expanded to include eosinophilic esophagitis, cyclic vomiting syndrome, Mallory Weiss tears, abdominal migraine, esophageal dilation, and subglottic stenosis. Although the exact cause of poor growth in NAA10-related neurodevelopmental syndrome probands is unclear and the degree of contribution to this problem by GI symptomatology remains uncertain, an analysis including nine G-tube or GJ-tube fed probands demonstrates that G/GJ-tubes are overall efficacious with respect to improvements in weight gain and caregiving. The choice to insert a gastrostomy or gastrojejunal tube to aid with weight gain is often a challenging decision to make for parents, who may alternatively choose to rely on oral feeding, caloric supplementation, calorie tracking, and feeding therapy. In this case, if NAA10-related neurodevelopmental syndrome children are not tracking above the failure to thrive (FTT) range past 1 year of age despite such efforts, the treating physicians should be consulted regarding possibly undergoing G-tube placement to avoid prolonged growth failure. If G-tubes are not immediately inducing weight gain after insertion, recommendations could include altering formula, increasing caloric input, or exchanging a G-tube for a GJ-tube by means of a minimally invasive procedure.

KBG syndrome: videoconferencing and use of artificial intelligence driven facial phenotyping in 25 new patients

Genetic variants in Ankyrin Repeat Domain 11 (ANKRD11) and deletions in 16q24.3 are known to cause KBG syndrome, a rare syndrome associated with craniofacial, intellectual, and neurobehavioral anomalies. We report 25 unpublished individuals from 22 families with molecularly confirmed diagnoses. Twelve individuals have de novo variants, three have inherited variants, and one is inherited from a parent with low-level mosaicism. The mode of inheritance was unknown for nine individuals. Twenty are truncating variants, and the remaining five are missense (three of which are found in one family). We present a protocol emphasizing the use of videoconference and artificial intelligence (AI) in collecting and analyzing data for this rare syndrome. A single clinician interviewed 25 individuals throughout eight countries. Participants' medical records were reviewed, and data was uploaded to the Human Disease Gene website using Human Phenotype Ontology (HPO) terms. Photos of the participants were analyzed by the GestaltMatcher and DeepGestalt, Face2Gene platform (FDNA Inc, USA) algorithms. Within our cohort, common traits included short stature, macrodontia, anteverted nares, wide nasal bridge, wide nasal base, thick eyebrows, synophrys and hypertelorism. Behavioral issues and global developmental delays were widely present. Neurologic abnormalities including seizures and/or EEG abnormalities were common (44%), suggesting that early detection and seizure prophylaxis could be an important point of intervention. Almost a quarter (24%) were diagnosed with attention deficit hyperactivity disorder and 28% were diagnosed with autism spectrum disorder. Based on the data, we provide a set of recommendations regarding diagnostic and treatment approaches for KBG syndrome.

SvAnna: efficient and accurate pathogenicity prediction of coding and regulatory structural variants in long-read genome sequencing

Structural variants (SVs) are implicated in the etiology of Mendelian diseases but have been systematically underascertained owing to sequencing technology limitations. Long-read sequencing enables comprehensive detection of SVs, but approaches for prioritization of candidate SVs are needed. Structural variant Annotation and analysis (SvAnna) assesses all classes of SVs and their intersection with transcripts and regulatory sequences, relating predicted effects on gene function with clinical phenotype data. SvAnna places 87% of deleterious SVs in the top ten ranks. The interpretable prioritizations offered by SvAnna will facilitate the widespread adoption of long-read sequencing in diagnostic genomics. SvAnna is available at https://github.com/TheJacksonLaboratory/SvAnn a .

Autosomal recessive SLC30A9 variants in a Proband with a Cerebro-Renal Syndrome and No Parental Consanguinity

An SLC30A9-associated cerebrorenal syndrome was first reported in consanguineous Bedouin kindred by Perez et al. in 2017. Although the function of the gene has not yet been fully elucidated, it may be implicated in Wnt signaling and nuclear regulation, as well as in cell and mitochondrial zinc regulation. In this research report, we present a female proband with two distinct, inherited autosomal recessive loss-of-function SLC30A9 variants from unrelated parents. To our knowledge, this is the first reported case of a possible SLC30A9-associated cerebrorenal syndrome in a nonconsanguineous family. Furthermore, a limited statistical analysis was conducted to identify possible allele frequency differences between populations. Our findings provide further support for an SLC30A9-associated cerebrorenal syndrome and may help clarify the gene's function through its possible disease association.

Naa12 compensates for Naa10 in mice in the amino-terminal acetylation pathway

Amino-terminal acetylation is catalyzed by a set of N-terminal acetyltransferases (NATs). The NatA complex (including X-linked Naa10 and Naa15) is the major acetyltransferase, with 40-50% of all mammalian proteins being potential substrates. However, the overall role of amino-terminal acetylation on a whole-organism level is poorly understood, particularly in mammals. Male mice lacking Naa10 show no globally apparent in vivo amino-terminal acetylation impairment and do not exhibit complete embryonic lethality. Rather Naa10 nulls display increased neonatal lethality, and the majority of surviving undersized mutants exhibit a combination of hydrocephaly, cardiac defects, homeotic anterior transformation, piebaldism, and urogenital anomalies. Naa12 is a previously unannotated Naa10-like paralog with NAT activity that genetically compensates for Naa10. Mice deficient for Naa12 have no apparent phenotype, whereas mice deficient for Naa10 and Naa12 display embryonic lethality. The discovery of Naa12 adds to the currently known machinery involved in amino-terminal acetylation in mice.

Is Persistent Motor or Vocal Tic Disorder a Milder Form of Tourette Syndrome?

Background

Persistent motor or vocal tic disorder (PMVT) has been hypothesized to be a forme fruste of Tourette syndrome (TS). Although the primary diagnostic criterion for PMVT (presence of motor or vocal tics, but not both) is clear, less is known about its clinical presentation.

Objective

The goals of this study were to compare the prevalence and number of comorbid psychiatric disorders, tic severity, age at tic onset, and family history for TS and PMVT.

Methods

We analyzed data from two independent cohorts using generalized linear equations and confirmed our findings using meta‐analyses, incorporating data from previously published literature.

Results

Rates of obsessive–compulsive disorder (OCD) and attention deficit hyperactivity disorder (ADHD) were lower in PMVT than in TS in all analyses. Other psychiatric comorbidities occurred with similar frequencies in PMVT and TS in both cohorts, although meta‐analyses suggested lower rates of most psychiatric disorders in PMVT compared with TS. ADHD and OCD increased the odds of comorbid mood, anxiety, substance use, and disruptive behaviors, and accounted for observed differences between PMVT and TS. Age of tic onset was approximately 2 years later, and tic severity was lower in PMVT than in TS. First‐degree relatives had elevated rates of TS, PMVT, OCD, and ADHD compared with population prevalences, with rates of TS equal to or greater than PMVT rates.

Conclusions

Our findings support the hypothesis that PMVT and TS occur along a clinical spectrum in which TS is a more severe and PMVT a less severe manifestation of a continuous neurodevelopmental tic spectrum disorder. © 2021 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society

Synaptic processes and immune-related pathways implicated in Tourette syndrome

Tourette syndrome (TS) is a neuropsychiatric disorder of complex genetic architecture involving multiple interacting genes. Here, we sought to elucidate the pathways that underlie the neurobiology of the disorder through genome-wide analysis. We analyzed genome-wide genotypic data of 3581 individuals with TS and 7682 ancestry-matched controls and investigated associations of TS with sets of genes that are expressed in particular cell types and operate in specific neuronal and glial functions. We employed a self-contained, set-based association method (SBA) as well as a competitive gene set method (MAGMA) using individual-level genotype data to perform a comprehensive investigation of the biological background of TS. Our SBA analysis identified three significant gene sets after Bonferroni correction, implicating ligand-gated ion channel signaling, lymphocytic, and cell adhesion and transsynaptic signaling processes. MAGMA analysis further supported the involvement of the cell adhesion and trans-synaptic signaling gene set. The lymphocytic gene set was driven by variants in FLT3, raising an intriguing hypothesis for the involvement of a neuroinflammatory element in TS pathogenesis. The indications of involvement of ligand-gated ion channel signaling reinforce the role of GABA in TS, while the association of cell adhesion and trans-synaptic signaling gene set provides additional support for the role of adhesion molecules in neuropsychiatric disorders. This study reinforces previous findings but also provides new insights into the neurobiology of TS.

Early-onset cerebellar ataxia in a patient with CMT2A2

A 9-yr 8-mo-old right-handed female presented with a history of gait difficulties, which first became apparent at age 9 mo of age, along with slurred speech and hand tremors while holding a tray. Her past medical history was significant for global developmental delay, and she was attending fourth grade special education classes. On examination, she had an ataxic gait, dysarthria, absent deep tendon reflexes, and flexor plantar responses. There were no signs of optic atrophy or hearing loss. Nerve conduction studies were consistent with an axonal neuropathy. A fascicular sural nerve biopsy showed a marked decrease of myelinated fibers larger than 6 µm in diameter as compared with an age-matched control. By electron microscopy, clusters of degenerating axonal mitochondria in both myelinated and unmyelinated fibers were frequently found. Whole-exome sequencing revealed a heterozygous c.314C > T (p.Thr105Met) missense variant in MFN2 in the patient but not in her mother. The father was unavailable for testing. The phenotypes with MFN2 variants can be quite variable, including intellectual disability, optic atrophy, auditory impairment, spinal atrophy with or without hydromyelia, and hydrocephalus. We report here that early onset ataxia with intellectual disability can also be associated with MFN2-related Charcot–Marie–Tooth, Type 2A2A diagnosis, the most common type of autosomal dominant axonal neuropathy.

Phen2Gene: rapid phenotype-driven gene prioritization for rare diseases

Human Phenotype Ontology (HPO) terms are increasingly used in diagnostic settings to aid in the characterization of patient phenotypes. The HPO annotation database is updated frequently and can provide detailed phenotype knowledge on various human diseases, and many HPO terms are now mapped to candidate causal genes with binary relationships. To further improve the genetic diagnosis of rare diseases, we incorporated these HPO annotations, gene-disease databases and gene-gene databases in a probabilistic model to build a novel HPO-driven gene prioritization tool, Phen2Gene. Phen2Gene accesses a database built upon this information called the HPO2Gene Knowledgebase (H2GKB), which provides weighted and ranked gene lists for every HPO term. Phen2Gene is then able to access the H2GKB for patient-specific lists of HPO terms or PhenoPacket descriptions supported by GA4GH (http://phenopackets.org/), calculate a prioritized gene list based on a probabilistic model and output gene-disease relationships with great accuracy. Phen2Gene outperforms existing gene prioritization tools in speed and acts as a real-time phenotype-driven gene prioritization tool to aid the clinical diagnosis of rare undiagnosed diseases. In addition to a command line tool released under the MIT license (https://github.com/WGLab/Phen2Gene), we also developed a web server and web service (https://phen2gene.wglab.org/) for running the tool via web interface or RESTful API queries. Finally, we have curated a large amount of benchmarking data for phenotype-to-gene tools involving 197 patients across 76 scientific articles and 85 patients' de-identified HPO term data from the Children's Hospital of Philadelphia.

VAC14 syndrome in two siblings with retinitis pigmentosa and neurodegeneration with brain iron accumulation

Whole-exome sequencing was used to identify the genetic etiology of a rapidly progressing neurological disease present in two of six siblings with early childhood onset of severe progressive spastic paraparesis and learning disabilities. A homozygous mutation (c.2005G>T, p, V669L) was found in VAC14, and the clinical phenotype is consistent with the recently described VAC14-related striatonigral degeneration, childhood-onset syndrome (SNDC) (MIM#617054). However, the phenotype includes a distinct clinical presentation of retinitis pigmentosa (RP), which has not previously been reported in association with VAC14 mutations. Brain magnetic resonance imaging (MRI) revealed abnormal magnetic susceptibility in the globus pallidus, which can be seen in neurodegeneration with brain iron accumulation (NBIA). RP is a group of inherited retinal diseases with phenotypic/genetic heterogeneity, and the pathophysiologic basis of RP is not completely understood but is thought to be due to a primary retinal photoreceptor cell degenerative process. Most cases of RP are seen in isolation (nonsyndromic); this is a report of RP in two siblings with VAC14-associated syndrome, and it is suggested that a connection between RP and VAC14-associated syndrome should be explored in future studies.

From Molecular Understanding to Organismal Biology of N-Terminal Acetyltransferases

In this issue of Structure, Deng et al. (2019) determine the structure of the yeast N-terminal acetyltransferases Naa10 and Naa50 in complex with Naa15 and demonstrate that Naa50 has negligible catalytic activity on its own but modulates Naa10/Naa15. This study provides insights into mechanisms involving amino-terminal acetylation of proteins.

Missense variants in TAF1 and developmental phenotypes: challenges of determining pathogenicity

We recently described a new neurodevelopmental syndrome (TAF1/MRXS33 intellectual disability syndrome) (MIM# 300966) caused by pathogenic variants involving the X-linked gene TAF1, which participates in RNA polymerase II transcription. The initial study reported eleven families, and the syndrome was defined as presenting early in life with hypotonia, facial dysmorphia, and developmental delay that evolved into intellectual disability (ID) and/or autism spectrum disorder (ASD). We have now identified an additional 27 families through a genotype-first approach. Familial segregation analysis, clinical phenotyping, and bioinformatics were capitalized on to assess potential variant pathogenicity, and molecular modelling was performed for those variants falling within structurally characterized domains of TAF1. A novel phenotypic clustering approach was also applied, in which the phenotypes of affected individuals were classified using 51 standardized Human Phenotype Ontology (HPO) terms. Phenotypes associated with TAF1 variants show considerable pleiotropy and clinical variability, but prominent among previously unreported effects were brain morphological abnormalities, seizures, hearing loss, and heart malformations. Our allelic series broadens the phenotypic spectrum of TAF1/MRXS33 intellectual disability syndrome and the range of TAF1 molecular defects in humans. It also illustrates the challenges for determining the pathogenicity of inherited missense variants, particularly for genes mapping to chromosome X. This article is protected by copyright. All rights reserved.

Phenotypic and biochemical analysis of an international cohort of individuals with variants in NAA10 and NAA15

N-alpha-acetylation is one of the most common co-translational protein modifications in humans and is essential for normal cell function. NAA10 encodes for the enzyme NAA10, which is the catalytic subunit in the N-terminal acetyltransferase A (NatA) complex. The auxiliary and regulatory subunits of the NatA complex are NAA15 and Huntington-interacting protein (HYPK), respectively. Through a genotype-first approach with exome sequencing, we identified and phenotypically characterized 30 individuals from 30 unrelated families with 17 different de novo or inherited, dominantly acting missense variants in NAA10 or NAA15. Clinical features of affected individuals include variable levels of intellectual disability, delayed speech and motor milestones and autism spectrum disorder. Additionally, some subjects present with mild craniofacial dysmorphology, congenital cardiac anomalies and seizures. One of the individuals is an 11-year-old boy with a frameshift variant in exon 7 of NAA10, who presents most notably with microphthalmia, which confirms a prior finding with a single family with Lenz microphthalmia syndrome. Biochemical analyses of variants as part of the human NatA complex, as well as enzymatic analyses with and without the HYPK regulatory subunit, help to explain some of the phenotypic differences seen among the different variants.

Paralog Studies Augment Gene Discovery: DDX and DHX Genes

Members of a paralogous gene family in which variation in one gene is known to cause disease are eight times more likely to also be associated with human disease. Recent studies have elucidated DHX30 and DDX3X as genes for which pathogenic variant alleles are involved in neurodevelopmental disorders. We hypothesized that variants in paralogous genes encoding members of the DExD/H-box RNA helicase superfamily might also underlie developmental delay and/or intellectual disability (DD and/or ID) disease phenotypes. Here we describe 15 unrelated individuals who have DD and/or ID, central nervous system (CNS) dysfunction, vertebral anomalies, and dysmorphic features and were found to have probably damaging variants in DExD/H-box RNA helicase genes. In addition, these individuals exhibit a variety of other tissue and organ system involvement including ocular, outer ear, hearing, cardiac, and kidney tissues. Five individuals with homozygous (one), compound-heterozygous (two), or de novo (two) missense variants in DHX37 were identified by exome sequencing. We identified ten total individuals with missense variants in three other DDX/DHX paralogs: DHX16 (four individuals), DDX54 (three individuals), and DHX34 (three individuals). Most identified variants are rare, predicted to be damaging, and occur at conserved amino acid residues. Taken together, these 15 individuals implicate the DExD/H-box helicases in both dominantly and recessively inherited neurodevelopmental phenotypes and highlight the potential for more than one disease mechanism underlying these disorders.

Enabling Global Clinical Collaborations on Identifiable Patient Data: The Minerva Initiative

The clinical utility of computational phenotyping for both genetic and rare diseases is increasingly appreciated; however, its true potential is yet to be fully realized. Alongside the growing clinical and research availability of sequencing technologies, precise deep and scalable phenotyping is required to serve unmet need in genetic and rare diseases. To improve the lives of individuals affected with rare diseases through deep phenotyping, global big data interrogation is necessary to aid our understanding of disease biology, assist diagnosis, and develop targeted treatment strategies. This includes the application of cutting-edge machine learning methods to image data. As with most digital tools employed in health care, there are ethical and data governance challenges associated with using identifiable personal image data. There are also risks with failing to deliver on the patient benefits of these new technologies, the biggest of which is posed by data siloing. The Minerva Initiative has been designed to enable the public good of deep phenotyping while mitigating these ethical risks. Its open structure, enabling collaboration and data sharing between individuals, clinicians, researchers and private enterprise, is key for delivering precision public health.

Interrogating the Genetic Determinants of Tourette's Syndrome and Other Tic Disorders Through Genome-Wide Association Studies

OBJECTIVE

Tourette's syndrome is polygenic and highly heritable. Genome-wide association study (GWAS) approaches are useful for interrogating the genetic architecture and determinants of Tourette's syndrome and other tic disorders. The authors conducted a GWAS meta-analysis and probed aggregated Tourette's syndrome polygenic risk to test whether Tourette's and related tic disorders have an underlying shared genetic etiology and whether Tourette's polygenic risk scores correlate with worst-ever tic severity and may represent a potential predictor of disease severity.

METHODS

GWAS meta-analysis, gene-based association, and genetic enrichment analyses were conducted in 4,819 Tourette's syndrome case subjects and 9,488 control subjects. Replication of top loci was conducted in an independent population-based sample (706 case subjects, 6,068 control subjects). Relationships between Tourette's polygenic risk scores (PRSs), other tic disorders, ascertainment, and tic severity were examined.

RESULTS

GWAS and gene-based analyses identified one genome-wide significant locus within FLT3 on chromosome 13, rs2504235, although this association was not replicated in the population-based sample. Genetic variants spanning evolutionarily conserved regions significantly explained 92.4% of Tourette's syndrome heritability. Tourette's-associated genes were significantly preferentially expressed in dorsolateral prefrontal cortex. Tourette's PRS significantly predicted both Tourette's syndrome and tic spectrum disorders status in the population-based sample. Tourette's PRS also significantly correlated with worst-ever tic severity and was higher in case subjects with a family history of tics than in simplex case subjects.

CONCLUSIONS

Modulation of gene expression through noncoding variants, particularly within cortico-striatal circuits, is implicated as a fundamental mechanism in Tourette's syndrome pathogenesis. At a genetic level, tic disorders represent a continuous spectrum of disease, supporting the unification of Tourette's syndrome and other tic disorders in future diagnostic schemata. Tourette's PRSs derived from sufficiently large samples may be useful in the future for predicting conversion of transient tics to chronic tic disorders, as well as tic persistence and lifetime tic severity.

Analysis of shared heritability in common disorders of the brain

Disorders of the brain can exhibit considerable epidemiological comorbidity and often share symptoms, provoking debate about their etiologic overlap. We quantified the genetic sharing of 25 brain disorders from genome-wide association studies of 265,218 patients and 784,643 control participants and assessed their relationship to 17 phenotypes from 1,191,588 individuals. Psychiatric disorders share common variant risk, whereas neurological disorders appear more distinct from one another and from the psychiatric disorders. We also identified significant sharing between disorders and a number of brain phenotypes, including cognitive measures. Further, we conducted simulations to explore how statistical power, diagnostic misclassification, and phenotypic heterogeneity affect genetic correlations. These results highlight the importance of common genetic variation as a risk factor for brain disorders and the value of heritability-based methods in understanding their etiology.

Truncating Variants in NAA15 Are Associated with Variable Levels of Intellectual Disability, Autism Spectrum Disorder, and Congenital Anomalies

N-alpha-acetylation is a common co-translational protein modification that is essential for normal cell function in humans. We previously identified the genetic basis of an X-linked infantile lethal Mendelian disorder involving a c.109T>C (p.Ser37Pro) missense variant in NAA10, which encodes the catalytic subunit of the N-terminal acetyltransferase A (NatA) complex. The auxiliary subunit of the NatA complex, NAA15, is the dimeric binding partner for NAA10. Through a genotype-first approach with whole-exome or genome sequencing (WES/WGS) and targeted sequencing analysis, we identified and phenotypically characterized 38 individuals from 33 unrelated families with 25 different de novo or inherited, dominantly acting likely gene disrupting (LGD) variants in NAA15. Clinical features of affected individuals with LGD variants in NAA15 include variable levels of intellectual disability, delayed speech and motor milestones, and autism spectrum disorder. Additionally, mild craniofacial dysmorphology, congenital cardiac anomalies, and seizures are present in some subjects. RNA analysis in cell lines from two individuals showed degradation of the transcripts with LGD variants, probably as a result of nonsense-mediated decay. Functional assays in yeast confirmed a deleterious effect for two of the LGD variants in NAA15. Further supporting a mechanism of haploinsufficiency, individuals with copy-number variant (CNV) deletions involving NAA15 and surrounding genes can present with mild intellectual disability, mild dysmorphic features, motor delays, and decreased growth. We propose that defects in NatA-mediated N-terminal acetylation (NTA) lead to variable levels of neurodevelopmental disorders in humans, supporting the importance of the NatA complex in normal human development.

Scikit-ribo Enables Accurate Estimation and Robust Modeling of Translation Dynamics at Codon Resolution

Ribosome profiling (Ribo-seq) is a powerful technique for measuring protein translation; however, sampling errors and biological biases are prevalent and poorly understood. Addressing these issues, we present Scikit-ribo (https://github.com/schatzlab/scikit-ribo), an open-source analysis package for accurate genome-wide A-site prediction and translation efficiency (TE) estimation from Ribo-seq and RNA sequencing data. Scikit-ribo accurately identifies A-site locations and reproduces codon elongation rates using several digestion protocols (r = 0.99). Next, we show that the commonly used reads per kilobase of transcript per million mapped reads-derived TE estimation is prone to biases, especially for low-abundance genes. Scikit-ribo introduces a codon-level generalized linear model with ridge penalty that correctly estimates TE, while accommodating variable codon elongation rates and mRNA secondary structure. This corrects the TE errors for over 2,000 genes in S. cerevisiae, which we validate using mass spectrometry of protein abundances (r = 0.81), and allows us to determine the Kozak-like sequence directly from Ribo-seq. We conclude with an analysis of coverage requirements needed for robust codon-level analysis and quantify the artifacts that can occur from cycloheximide treatment.

Genetic and phenotypic overlap of specific obsessive-compulsive and attention-deficit/hyperactive subtypes with Tourette syndrome

BACKGROUND

The unique phenotypic and genetic aspects of obsessive-compulsive (OCD) and attention-deficit/hyperactivity disorder (ADHD) among individuals with Tourette syndrome (TS) are not well characterized. Here, we examine symptom patterns and heritability of OCD and ADHD in TS families.

METHOD

OCD and ADHD symptom patterns were examined in TS patients and their family members (N = 3494) using exploratory factor analyses (EFA) for OCD and ADHD symptoms separately, followed by latent class analyses (LCA) of the resulting OCD and ADHD factor sum scores jointly; heritability and clinical relevance of the resulting factors and classes were assessed.

RESULTS

EFA yielded a 2-factor model for ADHD and an 8-factor model for OCD. Both ADHD factors (inattentive and hyperactive/impulsive symptoms) were genetically related to TS, ADHD, and OCD. The doubts, contamination, need for sameness, and superstitions factors were genetically related to OCD, but not ADHD or TS; symmetry/exactness and fear-of-harm were associated with TS and OCD while hoarding was associated with ADHD and OCD. In contrast, aggressive urges were genetically associated with TS, OCD, and ADHD. LCA revealed a three-class solution: few OCD/ADHD symptoms (LC1), OCD & ADHD symptoms (LC2), and symmetry/exactness, hoarding, and ADHD symptoms (LC3). LC2 had the highest psychiatric comorbidity rates (⩾50% for all disorders).

CONCLUSIONS

Symmetry/exactness, aggressive urges, fear-of-harm, and hoarding show complex genetic relationships with TS, OCD, and ADHD, and, rather than being specific subtypes of OCD, transcend traditional diagnostic boundaries, perhaps representing an underlying vulnerability (e.g. failure of top-down cognitive control) common to all three disorders.

KBG syndrome involving a single-nucleotide duplication in ANKRD11

KBG syndrome is a rare autosomal dominant genetic condition characterized by neurological involvement and distinct facial, hand, and skeletal features. More than 70 cases have been reported; however, it is likely that KBG syndrome is underdiagnosed because of lack of comprehensive characterization of the heterogeneous phenotypic features. We describe the clinical manifestations in a male currently 13 years of age, who exhibited symptoms including epilepsy, severe developmental delay, distinct facial features, and hand anomalies, without a positive genetic diagnosis. Subsequent exome sequencing identified a novel de novo heterozygous single base pair duplication (c.6015dupA) in ANKRD11, which was validated by Sanger sequencing. This single-nucleotide duplication is predicted to lead to a premature stop codon and loss of function in ANKRD11, thereby implicating it as contributing to the proband's symptoms and yielding a molecular diagnosis of KBG syndrome. Before molecular diagnosis, this syndrome was not recognized in the proband, as several key features of the disorder were mild and were not recognized by clinicians, further supporting the concept of variable expressivity in many disorders. Although a diagnosis of cerebral folate deficiency has also been given, its significance for the proband's condition remains uncertain.

SCN8A mutation in a child presenting with seizures and developmental delays

The SCN8A gene encodes the sodium voltage-gated channel alpha subunit 8. Mutations in this gene have been associated with early infantile epileptic encephalopathy type 13. With the use of whole-exome sequencing, a de novo missense mutation in SCN8A was identified in a 4-yr-old female who initially exhibited symptoms of epilepsy at the age of 5 mo that progressed to a severe condition with very little movement, including being unable to sit or walk on her own.

Autism Spectrum Symptoms in a Tourette's Disorder Sample

OBJECTIVE

Tourette's disorder (TD) and autism spectrum disorder (ASD) share clinical features and possibly an overlapping etiology. The aims of this study were to examine ASD symptom rates in participants with TD, and to characterize the relationships between ASD symptom patterns and TD, obsessive-compulsive disorder (OCD), and attention-deficit/hyperactivity disorder (ADHD).

METHOD

Participants with TD (n = 535) and their family members (n =234) recruited for genetic studies reported TD, OCD, and ADHD symptoms and completed the Social Responsiveness Scale Second Edition (SRS), which was used to characterize ASD symptoms.

RESULTS

SRS scores in participants with TD were similar to those observed in other clinical samples but lower than in ASD samples (mean SRS total raw score = 51; SD = 32.4). More children with TD met cut-off criteria for ASD (22.8%) than adults with TD (8.7%). The elevated rate in children was primarily due to high scores on the SRS Repetitive and Restricted Behaviors (RRB) subscale. Total SRS scores were correlated with TD (r = 0.27), OCD (r = 0.37), and ADHD (r = 0.44) and were higher among individuals with OCD symptom-based phenotypes than for those with tics alone.

CONCLUSION

Higher observed rates of ASD among children affected by TD may in part be due to difficulty in discriminating complex tics and OCD symptoms from ASD symptoms. Careful examination of ASD-specific symptom patterns (social communication vs. repetitive behaviors) is essential. Independent of ASD, the SRS may be a useful tool for identifying patients with TD with impairments in social communication that potentially place them at risk for bullying and other negative sequelae.

Rare Copy Number Variants in NRXN1 and CNTN6 Increase Risk for Tourette Syndrome

Tourette syndrome (TS) is a model neuropsychiatric disorder thought to arise from abnormal development and/or maintenance of cortico-striato-thalamo-cortical circuits. TS is highly heritable, but its underlying genetic causes are still elusive, and no genome-wide significant loci have been discovered to date. We analyzed a European ancestry sample of 2,434 TS cases and 4,093 ancestry-matched controls for rare (< 1% frequency) copy-number variants (CNVs) using SNP microarray data. We observed an enrichment of global CNV burden that was prominent for large (> 1 Mb), singleton events (OR = 2.28, 95% CI [1.39-3.79], p = 1.2 × 10-3) and known, pathogenic CNVs (OR = 3.03 [1.85-5.07], p = 1.5 × 10-5). We also identified two individual, genome-wide significant loci, each conferring a substantial increase in TS risk (NRXN1 deletions, OR = 20.3, 95% CI [2.6-156.2]; CNTN6 duplications, OR = 10.1, 95% CI [2.3-45.4]). Approximately 1% of TS cases carry one of these CNVs, indicating that rare structural variation contributes significantly to the genetic architecture of TS.

Identification of Two Heritable Cross-Disorder Endophenotypes for Tourette Syndrome

OBJECTIVE

Phenotypic heterogeneity in Tourette syndrome is partly due to complex genetic relationships among Tourette syndrome, obsessive-compulsive disorder (OCD), and attention deficit hyperactivity disorder (ADHD). Identifying symptom-based endophenotypes across diagnoses may aid gene-finding efforts.

METHOD

Assessments for Tourette syndrome, OCD, and ADHD symptoms were conducted in a discovery sample of 3,494 individuals recruited for genetic studies. Symptom-level factor and latent class analyses were conducted in Tourette syndrome families and replicated in an independent sample of 882 individuals. Classes were characterized by comorbidity rates and proportion of parents included. Heritability and polygenic load associated with Tourette syndrome, OCD, and ADHD were estimated.

RESULTS

The authors identified two cross-disorder symptom-based phenotypes across analyses: symmetry (symmetry, evening up, checking obsessions; ordering, arranging, counting, writing-rewriting compulsions, repetitive writing tics) and disinhibition (uttering syllables/words, echolalia/palilalia, coprolalia/copropraxia, and obsessive urges to offend/mutilate/be destructive). Heritability estimates for both endophenotypes were high and statistically significant (disinhibition factor=0.35, SE=0.03; symmetry factor=0.39, SE=0.03; symmetry class=0.38, SE=0.10). Mothers of Tourette syndrome probands had high rates of symmetry (49%) but not disinhibition (5%). Polygenic risk scores derived from a Tourette syndrome genome-wide association study (GWAS) were significantly associated with symmetry, while risk scores derived from an OCD GWAS were not. OCD polygenic risk scores were significantly associated with disinhibition, while Tourette syndrome and ADHD risk scores were not.

CONCLUSIONS

The analyses identified two heritable endophenotypes related to Tourette syndrome that cross traditional diagnostic boundaries. The symmetry phenotype correlated with Tourette syndrome polygenic load and was present in otherwise Tourette-unaffected mothers, suggesting that this phenotype may reflect additional Tourette syndrome (rather than OCD) genetic liability that is not captured by traditional DSM-based diagnoses.

Whole genome sequencing of one complex pedigree illustrates challenges with genomic medicine

Background

Human Phenotype Ontology (HPO) has risen as a useful tool for precision medicine by providing a standardized vocabulary of phenotypic abnormalities to describe presentations of human pathologies; however, there have been relatively few reports combining whole genome sequencing (WGS) and HPO, especially in the context of structural variants.

Methods

We illustrate an integrative analysis of WGS and HPO using an extended pedigree, which involves Prader–Willi Syndrome (PWS), hereditary hemochromatosis (HH), and dysautonomia-like symptoms. A comprehensive WGS pipeline was used to ensure reliable detection of genomic variants. Beyond variant filtering, we pursued phenotypic prioritization of candidate genes using Phenolyzer.

Results

Regarding PWS, WGS confirmed a 5.5 Mb de novo deletion of the parental allele at 15q11.2 to 15q13.1. Phenolyzer successfully returned the diagnosis of PWS, and pinpointed clinically relevant genes in the deletion. Further, Phenolyzer revealed how each of the genes is linked with the phenotypes represented by HPO terms. For HH, WGS identified a known disease variant (p.C282Y) in HFE of an affected female. Analysis of HPO terms alone fails to provide a correct diagnosis, but Phenolyzer successfully revealed the phenotype-genotype relationship using a disease-centric approach. Finally, Phenolyzer also revealed the complexity behind dysautonomia-like symptoms, and seven variants that might be associated with the phenotypes were identified by manual filtering based on a dominant inheritance model.

Conclusions

The integration of WGS and HPO can inform comprehensive molecular diagnosis for patients, eliminate false positives and reveal novel insights into undiagnosed diseases. Due to extreme heterogeneity and insufficient knowledge of human diseases, it is also important that phenotypic and genomic data are standardized and shared simultaneously.

Electronic supplementary material

The online version of this article (doi:10.1186/s12920-017-0246-5) contains supplementary material, which is available to authorized users.

Proteomic and genomic characterization of a yeast model for Ogden syndrome

Naa10 is an N^α -terminal acetyltransferase that, in a complex with its auxiliary subunit Naa15, co-translationally acetylates the α-amino group of newly synthetized proteins as they emerge from the ribosome. Roughly 40-50% of the human proteome is acetylated by Naa10, rendering this an enzyme one of the most broad substrate ranges known. Recently, we reported an X-linked disorder of infancy, Ogden syndrome, in two families harbouring a c.109 T > C (p.Ser37Pro) variant in NAA10. In the present study we performed in-depth characterization of a yeast model of Ogden syndrome. Stress tests and proteomic analyses suggest that the S37P mutation disrupts Naa10 function and reduces cellular fitness during heat shock, possibly owing to dysregulation of chaperone expression and accumulation. Microarray and RNA-seq revealed a pseudo-diploid gene expression profile in ΔNaa10 cells, probably responsible for a mating defect. In conclusion, the data presented here further support the disruptive nature of the S37P/Ogden mutation and identify affected cellular processes potentially contributing to the severe phenotype seen in Ogden syndrome. Data are available via GEO under identifier GSE86482 or with ProteomeXchange under identifier PXD004923. © 2016 The Authors. Yeast published by John Wiley & Sons, Ltd.

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.

Social disinhibition is a heritable subphenotype of tics in Tourette syndrome

OBJECTIVE

To identify heritable symptom-based subtypes of Tourette syndrome (TS).

METHODS

Forty-nine motor and phonic tics were examined in 3,494 individuals (1,191 TS probands and 2,303 first-degree relatives). Item-level exploratory factor and latent class analyses (LCA) were used to identify tic-based subtypes. Heritabilities of the subtypes were estimated, and associations with clinical characteristics were examined.

RESULTS

A 6-factor exploratory factor analysis model provided the best fit, which paralleled the somatotopic representation of the basal ganglia, distinguished simple from complex tics, and separated out socially disinhibited and compulsive tics. The 5-class LCA model best distinguished among the following groups: unaffected, simple tics, intermediate tics without social disinhibition, intermediate with social disinhibition, and high rates of all tic types. Across models, a phenotype characterized by high rates of social disinhibition emerged. This phenotype was associated with increased odds of comorbid psychiatric disorders, in particular, obsessive-compulsive disorder and attention-deficit/hyperactivity disorder, earlier age at TS onset, and increased tic severity. The heritability estimate for this phenotype based on the LCA was 0.53 (SE 0.08, p 1.7 × 10(-18)).

CONCLUSIONS

Expanding on previous modeling approaches, a series of TS-related phenotypes, including one characterized by high rates of social disinhibition, were identified. These phenotypes were highly heritable and may reflect underlying biological networks more accurately than traditional diagnoses, thus potentially aiding future genetic, imaging, and treatment studies.

Long-read sequencing and de novo assembly of a Chinese genome

Short-read sequencing has enabled the de novo assembly of several individual human genomes, but with inherent limitations in characterizing repeat elements. Here we sequence a Chinese individual HX1 by single-molecule real-time (SMRT) long-read sequencing, construct a physical map by NanoChannel arrays and generate a de novo assembly of 2.93 Gb (contig N50: 8.3 Mb, scaffold N50: 22.0 Mb, including 39.3 Mb N-bases), together with 206 Mb of alternative haplotypes. The assembly fully or partially fills 274 (28.4%) N-gaps in the reference genome GRCh38. Comparison to GRCh38 reveals 12.8 Mb of HX1-specific sequences, including 4.1 Mb that are not present in previously reported Asian genomes. Furthermore, long-read sequencing of the transcriptome reveals novel spliced genes that are not annotated in GENCODE and are missed by short-read RNA-Seq. Our results imply that improved characterization of genome functional variation may require the use of a range of genomic technologies on diverse human populations.

Short-read sequencing has inherent limitations in the characterisation of long repeat elements. Shi and Guo et al. combine single-molecule real-time sequencing and IrysChip to construct a Chinese reference genome that fills many gaps in the reference genome, and identify novel spliced genes.

Genome-wide variant analysis of simplex autism families with an integrative clinical-bioinformatics pipeline

Autism spectrum disorders (ASDs) are a group of developmental disabilities that affect social interaction and communication and are characterized by repetitive behaviors. There is now a large body of evidence that suggests a complex role of genetics in ASDs, in which many different loci are involved. Although many current population-scale genomic studies have been demonstrably fruitful, these studies generally focus on analyzing a limited part of the genome or use a limited set of bioinformatics tools. These limitations preclude the analysis of genome-wide perturbations that may contribute to the development and severity of ASD-related phenotypes. To overcome these limitations, we have developed and utilized an integrative clinical and bioinformatics pipeline for generating a more complete and reliable set of genomic variants for downstream analyses. Our study focuses on the analysis of three simplex autism families consisting of one affected child, unaffected parents, and one unaffected sibling. All members were clinically evaluated and widely phenotyped. Genotyping arrays and whole-genome sequencing were performed on each member, and the resulting sequencing data were analyzed using a variety of available bioinformatics tools. We searched for rare variants of putative functional impact that were found to be segregating according to de novo, autosomal recessive, X-linked, mitochondrial, and compound heterozygote transmission models. The resulting candidate variants included three small heterozygous copy-number variations (CNVs), a rare heterozygous de novo nonsense mutation in MYBBP1A located within exon 1, and a novel de novo missense variant in LAMB3. Our work demonstrates how more comprehensive analyses that include rich clinical data and whole-genome sequencing data can generate reliable results for use in downstream investigations.

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.

SeqHBase: a big data toolset for family based sequencing data analysis

Background

Whole-genome sequencing (WGS) and whole-exome sequencing (WES) technologies are increasingly used to identify disease-contributing mutations in human genomic studies. It can be a significant challenge to process such data, especially when a large family or cohort is sequenced. Our objective was to develop a big data toolset to efficiently manipulate genome-wide variants, functional annotations and coverage, together with conducting family based sequencing data analysis.

Methods

Hadoop is a framework for reliable, scalable, distributed processing of large data sets using MapReduce programming models. Based on Hadoop and HBase, we developed SeqHBase, a big data-based toolset for analysing family based sequencing data to detect de novo, inherited homozygous, or compound heterozygous mutations that may contribute to disease manifestations. SeqHBase takes as input BAM files (for coverage at every site), variant call format (VCF) files (for variant calls) and functional annotations (for variant prioritisation).

Results

We applied SeqHBase to a 5-member nuclear family and a 10-member 3-generation family with WGS data, as well as a 4-member nuclear family with WES data. Analysis times were almost linearly scalable with number of data nodes. With 20 data nodes, SeqHBase took about 5 secs to analyse WES familial data and approximately 1 min to analyse WGS familial data.

Conclusions

These results demonstrate SeqHBase's high efficiency and scalability, which is necessary as WGS and WES are rapidly becoming standard methods to study the genetics of familial disorders.

Accounting for uncertainty in DNA sequencing data

Science is defined in part by an honest exposition of the uncertainties that arise in measurements and propagate through calculations and inferences, so that the reliabilities of its conclusions are made apparent. The recent rapid development of high-throughput DNA sequencing technologies has dramatically increased the number of measurements made at the biochemical and molecular level. These data come from many different DNA-sequencing technologies, each with their own platform-specific errors and biases, which vary widely. Several statistical studies have tried to measure error rates for basic determinations, but there are no general schemes to project these uncertainties so as to assess the surety of the conclusions drawn about genetic, epigenetic, and more general biological questions. We review here the state of uncertainty quantification in DNA sequencing applications, describe sources of error, and propose methods that can be used for accounting and propagating these errors and their uncertainties through subsequent calculations.

Biochemical and cellular analysis of Ogden syndrome reveals downstream Nt-acetylation defects

The X-linked lethal Ogden syndrome was the first reported human genetic disorder associated with a mutation in an N-terminal acetyltransferase (NAT) gene. The affected males harbor an Ser37Pro (S37P) mutation in the gene encoding Naa10, the catalytic subunit of NatA, the major human NAT involved in the co-translational acetylation of proteins. Structural models and molecular dynamics simulations of the human NatA and its S37P mutant highlight differences in regions involved in catalysis and at the interface between Naa10 and the auxiliary subunit hNaa15. Biochemical data further demonstrate a reduced catalytic capacity and an impaired interaction between hNaa10 S37P and Naa15 as well as Naa50 (NatE), another interactor of the NatA complex. N-Terminal acetylome analyses revealed a decreased acetylation of a subset of NatA and NatE substrates in Ogden syndrome cells, supporting the genetic findings and our hypothesis regarding reduced Nt-acetylation of a subset of NatA/NatE-type substrates as one etiology for Ogden syndrome. Furthermore, Ogden syndrome fibroblasts display abnormal cell migration and proliferation capacity, possibly linked to a perturbed retinoblastoma pathway. N-Terminal acetylation clearly plays a role in Ogden syndrome, thus revealing the in vivo importance of N-terminal acetylation in human physiology and disease.

"Genotype-first" approaches on a curious case of idiopathic progressive cognitive decline

BACKGROUND

In developing countries, many cases with rare neurological diseases remain undiagnosed due to limited diagnostic experience. We encountered a case in China where two siblings both began to develop idiopathic progressive cognitive decline starting from age six, and were suspected to have an undiagnosed neurological disease.

METHODS

Initial clinical assessments included review of medical history, comprehensive physical examination, genetic testing for metabolic diseases, blood tests and brain imaging. We performed exome sequencing with Agilent SureSelect exon capture and Illumina HiSeq2000 platform, followed by variant annotation and selection of rare, shared mutations that fit a recessive model of inheritance. To assess functional impacts of candidate variants, we performed extensive biochemical tests in blood and urine, and examined their possible roles by protein structure modeling.

RESULTS

Exome sequencing identified NAGLU as the most likely candidate gene with compound heterozygous mutations (chr17:40695717C > T and chr17:40693129A > G in hg19 coordinate), which were documented to be pathogenic. Sanger sequencing confirmed the recessive patterns of inheritance, leading to a genetic diagnosis of Sanfilippo syndrome (mucopolysaccharidosis IIIB). Biochemical tests confirmed the complete loss of activity of alpha-N-acetylglucosaminidase (encoded by NAGLU) in blood, as well as significantly elevated dermatan sulfate and heparan sulfate in urine. Structure modeling revealed the mechanism on how the two variants affect protein structural stability.

CONCLUSIONS

Successful diagnosis of a rare genetic disorder with an atypical phenotypic presentation confirmed that such "genotype-first" approaches can particularly succeed in areas of the world with insufficient medical genetics expertise and with cost-prohibitive in-depth phenotyping.

Reducing INDEL calling errors in whole genome and exome sequencing data

Background

INDELs, especially those disrupting protein-coding regions of the genome, have been strongly associated with human diseases. However, there are still many errors with INDEL variant calling, driven by library preparation, sequencing biases, and algorithm artifacts.

Methods

We characterized whole genome sequencing (WGS), whole exome sequencing (WES), and PCR-free sequencing data from the same samples to investigate the sources of INDEL errors. We also developed a classification scheme based on the coverage and composition to rank high and low quality INDEL calls. We performed a large-scale validation experiment on 600 loci, and find high-quality INDELs to have a substantially lower error rate than low-quality INDELs (7% vs. 51%).

Results

Simulation and experimental data show that assembly based callers are significantly more sensitive and robust for detecting large INDELs (>5 bp) than alignment based callers, consistent with published data. The concordance of INDEL detection between WGS and WES is low (53%), and WGS data uniquely identifies 10.8-fold more high-quality INDELs. The validation rate for WGS-specific INDELs is also much higher than that for WES-specific INDELs (84% vs. 57%), and WES misses many large INDELs. In addition, the concordance for INDEL detection between standard WGS and PCR-free sequencing is 71%, and standard WGS data uniquely identifies 6.3-fold more low-quality INDELs. Furthermore, accurate detection with Scalpel of heterozygous INDELs requires 1.2-fold higher coverage than that for homozygous INDELs. Lastly, homopolymer A/T INDELs are a major source of low-quality INDEL calls, and they are highly enriched in the WES data.

Conclusions

Overall, we show that accuracy of INDEL detection with WGS is much greater than WES even in the targeted region. We calculated that 60X WGS depth of coverage from the HiSeq platform is needed to recover 95% of INDELs detected by Scalpel. While this is higher than current sequencing practice, the deeper coverage may save total project costs because of the greater accuracy and sensitivity. Finally, we investigate sources of INDEL errors (for example, capture deficiency, PCR amplification, homopolymers) with various data that will serve as a guideline to effectively reduce INDEL errors in genome sequencing.

Electronic supplementary material

The online version of this article (doi:10.1186/s13073-014-0089-z) contains supplementary material, which is available to authorized users.