Denisovan, modern human and mouse TNFAIP3 alleles tune A20 phosphorylation and immunity.

Resisting and tolerating microbes are alternative strategies to survive infection, but little is known about the evolutionary mechanisms controlling this balance. Here genomic analyses of anatomically modern humans, extinct Denisovan hominins and mice revealed a TNFAIP3 allelic series with alterations in the encoded immune response inhibitor A20. Each TNFAIP3 allele encoded substitutions at non-catalytic residues of the ubiquitin protease OTU domain that diminished IκB kinase-dependent phosphorylation and activation of A20. Two TNFAIP3 alleles encoding A20 proteins with partial phosphorylation deficits seemed to be beneficial by increasing immunity without causing spontaneous inflammatory disease: A20 T108A;I207L, originating in Denisovans and introgressed in modern humans throughout Oceania, and A20 I325N, from an N-ethyl-N-nitrosourea (ENU)-mutagenized mouse strain. By contrast, a rare human TNFAIP3 allele encoding an A20 protein with 95% loss of phosphorylation, C243Y, caused spontaneous inflammatory disease in humans and mice. Analysis of the partial-phosphorylation A20 I325N allele in mice revealed diminished tolerance of bacterial lipopolysaccharide and poxvirus inoculation as tradeoffs for enhanced immunity.

The Australian Genomics Mitochondrial Flagship: A National Program Delivering Mitochondrial Diagnoses.

PURPOSE: Families living with mitochondrial diseases (MD) often endure prolonged diagnostic journeys and invasive testing, yet many remain without a molecular diagnosis. The Australian Genomics Mitochondrial flagship, comprising clinicians, diagnostic, and research scientists, conducted a prospective national study to identify the diagnostic utility of singleton genomic sequencing using blood samples. METHODS: 140 children and adults living with suspected MD were recruited using modified Nijmegen criteria (MNC) and randomized to either exome + mtDNA sequencing (ES+mtDNAseq) or genome sequencing (GS). RESULTS: Diagnostic yield was 55% (n=77) with variants in nuclear (n=37) and mtDNA (n=18) MD genes, as well as phenocopy genes (n=22). A nuclear gene etiology was identified in 77% of diagnoses, irrespective of disease onset. Diagnostic rates were higher in pediatric-onset (71%) than adult-onset (31%) cases, and comparable in children with non-European (78%) versus European (67%) ancestry. For children, higher MNC scores correlated with increased diagnostic yield and fewer diagnoses in phenocopy genes. Additionally, three adult patients had a mtDNA deletion discovered in skeletal muscle that was not initially identified in blood. CONCLUSION: Genomic sequencing from blood can simplify the diagnostic pathway for individuals living with suspected MD, especially those with childhood onset diseases and high MNC scores.

Biallelic pathogenic variants in COX11 are associated with an infantile-onset mitochondrial encephalopathy.

Primary mitochondrial diseases are a group of genetically and clinically heterogeneous disorders resulting from oxidative phosphorylation (OXPHOS) defects. COX11 encodes a copper chaperone that participates in the assembly of complex IV and has not been previously linked to human disease. In a previous study, we identified that COX11 knockdown decreased cellular adenosine triphosphate (ATP) derived from respiration, and that ATP levels could be restored with coenzyme Q10 (CoQ10 ) supplementation. This finding is surprising since COX11 has no known role in CoQ10 biosynthesis. Here, we report a novel gene-disease association by identifying biallelic pathogenic variants in COX11 associated with infantile-onset mitochondrial encephalopathies in two unrelated families using trio genome and exome sequencing. Functional studies showed that mutant COX11 fibroblasts had decreased ATP levels which could be rescued by CoQ10 . These results not only suggest that COX11 variants cause defects in energy production but reveal a potential metabolic therapeutic strategy for patients with COX11 variants.

De Novo Variants Disrupting the HX Repeat Motif of ATN1 Cause a Recognizable Non-Progressive Neurocognitive Syndrome.

Polyglutamine expansions in the transcriptional co-repressor Atrophin-1, encoded by ATN1, cause the neurodegenerative condition dentatorubral-pallidoluysian atrophy (DRPLA) via a proposed novel toxic gain of function. We present detailed phenotypic information on eight unrelated individuals who have de novo missense and insertion variants within a conserved 16-amino-acid "HX repeat" motif of ATN1. Each of the affected individuals has severe cognitive impairment and hypotonia, a recognizable facial gestalt, and variable congenital anomalies. However, they lack the progressive symptoms typical of DRPLA neurodegeneration. To distinguish this subset of affected individuals from the DRPLA diagnosis, we suggest using the term CHEDDA (congenital hypotonia, epilepsy, developmental delay, digit abnormalities) to classify the condition. CHEDDA-related variants alter the particular structural features of the HX repeat motif, suggesting that CHEDDA results from perturbation of the structural and functional integrity of the HX repeat. We found several non-homologous human genes containing similar motifs of eight to 10 HX repeat sequences, including RERE, where disruptive variants in this motif have also been linked to a separate condition that causes neurocognitive and congenital anomalies. These findings suggest that perturbation of the HX motif might explain other Mendelian human conditions.

Seave: a comprehensive web platform for storing and interrogating human genomic variation.

Motivation: Genome sequencing has had a remarkable impact on our ability to study the effects of human genetic variation, however, variant interpretation remains the major bottleneck. Understanding the potential impact of variants, including structural variants, requires extensive annotation from disparate sources of knowledge, and in silico prediction algorithms. Results: We introduce Seave, an intuitive web platform that enables all types of variants to be securely stored, annotated and filtered. Variants are annotated with allele frequencies and pathogenicity assessments from many popular databases and in silico pathogenicity prediction scores. Seave enables filtering of variants with specific inheritance patterns, including somatic variants, by quality, allele frequencies and gene lists which can be curated and saved. Seave was made for whole genome data and is capable of storing and querying copy number and structural variants. Availability and implementation: To demo Seave with public data, see https://www.seave.bio. Source code is available at http://code.seave.bio and extensive documentation is available at http://documentation.seave.bio. Seave can be locally installed on an Apache server with PHP and MySQL, or we provide an Amazon Machine Image for quick deployment. For commercial and clinical diagnostic licensing, contact the corresponding author. Supplementary information: Supplementary data are available at Bioinformatics online.

The diagnostic utility of genome sequencing in a pediatric cohort with suspected mitochondrial disease.

PURPOSE: The utility of genome sequencing (GS) in the diagnosis of suspected pediatric mitochondrial disease (MD) was investigated. METHODS: An Australian cohort of 40 pediatric patients with clinical features suggestive of MD were classified using the modified Nijmegen mitochondrial disease severity scoring into definite (17), probable (17), and possible (6) MD groups. Trio GS was performed using DNA extracted from patient and parent blood. Data were analyzed for single-nucleotide variants, indels, mitochondrial DNA variants, and structural variants. RESULTS: A definitive MD gene molecular diagnosis was made in 15 cases and a likely MD molecular diagnosis in a further five cases. Causative mitochondrial DNA (mtDNA) variants were identified in four of these cases. Three potential novel MD genes were identified. In seven cases, causative variants were identified in known disease genes with no previous evidence of causing a primary MD. Diagnostic rates were higher in patients classified as having definite MD. CONCLUSION: GS efficiently identifies variants in MD genes of both nuclear and mitochondrial origin. A likely molecular diagnosis was identified in 67% of cases and a definitive molecular diagnosis achieved in 55% of cases. This study highlights the value of GS for a phenotypically and genetically heterogeneous disorder like MD.

Reanalysis and optimisation of bioinformatic pipelines is critical for mutation detection.

Rapid advances in genomic technologies have facilitated the identification pathogenic variants causing human disease. We report siblings with developmental and epileptic encephalopathy due to a novel, shared heterozygous pathogenic 13 bp duplication in SYNGAP1 (c.435_447dup, p.(L150Vfs*6)) that was identified by whole genome sequencing (WGS). The pathogenic variant had escaped earlier detection via two methodologies: whole exome sequencing and high-depth targeted sequencing. Both technologies had produced reads carrying the variant, however, they were either not aligned due to the size of the insertion or aligned to multiple major histocompatibility complex (MHC) regions in the hg19 reference genome, making the critical reads unavailable for variant calling. The WGS pipeline followed different protocols, including alignment of reads to the GRCh37 reference genome, which lacks the additional MHC contigs. Our findings highlight the benefit of using orthogonal clinical bioinformatic pipelines and all relevant inheritance patterns to re-analyze genomic data in undiagnosed patients.

MECR Mutations Cause Childhood-Onset Dystonia and Optic Atrophy, a Mitochondrial Fatty Acid Synthesis Disorder.

Mitochondrial fatty acid synthesis (mtFAS) is an evolutionarily conserved pathway essential for the function of the respiratory chain and several mitochondrial enzyme complexes. We report here a unique neurometabolic human disorder caused by defective mtFAS. Seven individuals from five unrelated families presented with childhood-onset dystonia, optic atrophy, and basal ganglia signal abnormalities on MRI. All affected individuals were found to harbor recessive mutations in MECR encoding the mitochondrial trans-2-enoyl-coenzyme A-reductase involved in human mtFAS. All six mutations are extremely rare in the general population, segregate with the disease in the families, and are predicted to be deleterious. The nonsense c.855T>G (p.Tyr285∗), c.247_250del (p.Asn83Hisfs∗4), and splice site c.830+2_830+3insT mutations lead to C-terminal truncation variants of MECR. The missense c.695G>A (p.Gly232Glu), c.854A>G (p.Tyr285Cys), and c.772C>T (p.Arg258Trp) mutations involve conserved amino acid residues, are located within the cofactor binding domain, and are predicted by structural analysis to have a destabilizing effect. Yeast modeling and complementation studies validated the pathogenicity of the MECR mutations. Fibroblast cell lines from affected individuals displayed reduced levels of both MECR and lipoylated proteins as well as defective respiration. These results suggest that mutations in MECR cause a distinct human disorder of the mtFAS pathway. The observation of decreased lipoylation raises the possibility of a potential therapeutic strategy.

Beyond the panel: preconception screening in consanguineous couples using the TruSight One "clinical exome".

PURPOSE: To provide proof of concept by broadening preconception screening beyond targeted testing to inform reproductive risk in consanguineous couples. METHODS: Consanguineous couples were screened for autosomal recessive and X-linked disorders using the TruSight One panel of 4,813 genes associated with human disease. RESULTS: We recruited 22 couples, of whom 15 elected to have sequencing. We found four couples to be at risk of autosomal recessive disorders, including one with a child affected by Poretti-Boltshauser syndrome (a diagnosis not made prior to the study) and another previously known to carry a ß-globin variant. Two couples were found to carry variants unrelated to known family history. These variants were in the genes C5orf42 (associated with Joubert syndrome and orofaciodigital syndrome) and GYS2 (associated with glycogen synthase deficiency). One known variant was not detected-a single exon deletion in FAM20C. We would not expect to identify this variant with the methodology employed. Of the four variants identified, only the ß-globin variant would have been found using available commercial preconception screening panels. CONCLUSION: Preconception screening of consanguineous couples for recessive and X-linked disorders using genomic sequencing is practicable, and is likely to detect many more at-risk couples than any targeted panel could achieve. A couples-based approach greatly reduces the associated analysis and counselling burden.

Genome sequencing as a first-line genetic test in familial dilated cardiomyopathy.

PURPOSE: We evaluated genome sequencing (GS) as an alternative to multigene panel sequencing (PS) for genetic testing in dilated cardiomyopathy (DCM). METHODS: Forty-two patients with familial DCM underwent PS and GS, and detection rates of rare single-nucleotide variants and small insertions/deletions in panel genes were compared. Loss-of-function variants in 406 cardiac-enriched genes were evaluated, and an assessment of structural variation was performed. RESULTS: GS provided broader and more uniform coverage than PS, with high concordance for rare variant detection in panel genes. GS identified all PS-identified pathogenic or likely pathogenic variants as well as two additional likely pathogenic variants: one was missed by PS due to low coverage, the other was a known disease-causing variant in a gene not included on the panel. No loss-of-function variants in the extended gene set met clinical criteria for pathogenicity. One BAG3 structural variant was classified as pathogenic. CONCLUSION: Our data support the use of GS for genetic testing in DCM, with high variant detection accuracy and a capacity to identify structural variants. GS provides an opportunity to go beyond suites of established disease genes, but the incremental yield of clinically actionable variants is limited by a paucity of genetic and functional evidence for DCM association.

Cryptic intronic NBAS variant reveals the genetic basis of recurrent liver failure in a child.

BACKGROUND: In almost half of patients with acute liver failure the cause is unknown, making targeted treatment and decisions about liver transplantation a challenge. Monogenic disorders may contribute to a significant proportion of these undiagnosed patients, and so the incorporation of technologies such as next generation sequencing (NGS) in the clinic could aid in providing a definitive diagnosis. However, this technology may present a major challenge in interpretation of sequence variants, particularly those in non-coding regions. RESULTS: In this report we describe a case of Infantile liver failure syndrome 2 (ILFS2; MIM 616483) due to novel bi-allelic variants in the NBAS gene. A missense variant NM_015909.3(NBAS):c.2617C > T, NP_056993.2(NBAS):p.(Arg873Trp) was identified by whole genome sequencing (WGS). By combining WGS and reverse transcription-polymerase chain reaction (RT-PCR) we were able to identify a novel deep intronic variant, NM_015909.3(NBAS):c.2423 + 404G > C, leading to the inclusion of a pseudo-exon. This mechanism has not been described previously in this syndrome. CONCLUSIONS: This study highlights the utility of analyzing NGS data in conjunction with investigating complementary DNA (cDNA) using techniques such as RT-PCR for detection of variants that otherwise would be likely to be missed in common NGS bioinformatic analysis pipelines. Combining these approaches, particularly when the phenotype match is strong, could lead to an increase in the diagnostic yield in acute liver failure and thus aid in targeted treatment, accurate genetic counseling and restoration of reproductive confidence.

Increased Diagnostic Yield of Spastic Paraplegia with or Without Cerebellar Ataxia Through Whole-Genome Sequencing.

Inherited disorders of spasticity or ataxia exist on a spectrum with overlapping causative genes and phenotypes. We investigated the use of whole-genome sequencing (WGS) to detect a genetic cause when considering this spectrum of disorders as a single group. We recruited 18 Korean individuals with spastic paraplegia with or without cerebellar ataxia in whom common causes of hereditary cerebellar ataxia and hereditary spastic paraplegia had been excluded. We performed WGS with analysis for single nucleotide variants, small insertions and deletions, copy number variants (CNVs), structural variants (SVs) and intronic variants. Disease-relevant variants were identified in ABCD1 (n = 3), CAPN1 (n = 2), NIPA1 (n = 1) and PLA2G6 (n = 1) for 7/18 patients (38.9%). A 'reverse phenotyping' approach was used to clarify the diagnosis in individuals with PLA2G6 and ABCD1 variants. One of the ABCD1 disease-relevant variants was detected on analysis for intronic variants. No CNV or SV causes were found. The two males with ABCD1 variants were initiated on monitoring for adrenal dysfunction. This is one of only a few studies to analyse spastic-ataxias as a continuous spectrum using a single approach. The outcome was improved diagnosis of unresolved cases for which common genetic causes had been excluded. This includes the detection of ABCD1 variants which had management implications. Therefore, WGS may be particularly relevant to diagnosing spastic ataxias given the large number of genes associated with this condition and the relatively high diagnostic yield.

Brief Report: Potent clinical and radiological response to larotrectinib in TRK fusion-driven high-grade glioma.

Genes encoding TRK are oncogenic drivers in multiple tumour types including infantile fibrosarcoma, papillary thyroid cancer and high-grade gliomas (HGG). TRK fusions have a critical role in tumourigenesis in 40% of infant HGG. Here we report the first case of a TRK fusion-driven HGG treated with larotrectinib-the first selective pan-TRK inhibitor in clinical development. This 3-year-old girl had failed multiple therapies including chemotherapy and radiotherapy. Tumour profiling confirmed an ETV6-NTRK3 fusion. Treatment with larotrectinib led to rapid clinical improvement with near total resolution of primary and metastatic lesions on MRI imaging. This is the first report of a TRK fusion glioma successfully treated with a TRK inhibitor.

Whole-exome sequencing reanalysis at 12 months boosts diagnosis and is cost-effective when applied early in Mendelian disorders.

PURPOSE: Whole-exome sequencing (WES) has revolutionized Mendelian diagnostics, however, there is no consensus on the timing of data review in undiagnosed individuals and only preliminary data on the cost-effectiveness of this technology. We aimed to assess the utility of WES data reanalysis for diagnosis in Mendelian disorders and to analyze the cost-effectiveness of this technology compared with a traditional diagnostic pathway. METHODS: WES was applied to a cohort of 54 patients from 37 families with a variety of Mendelian disorders to identify the genetic etiology. Reanalysis was performed after 12 months with an improved WES diagnostic pipeline. A comparison was made between costs of a modeled WES pathway and a traditional diagnostic pathway in a cohort with intellectual disability (ID). RESULTS: Reanalysis of WES data at 12 months improved diagnostic success from 30 to 41% due to interim publication of disease genes, expanded phenotype data from referrer, and an improved bioinformatics pipeline. Cost analysis on the ID cohort showed average cost savings of US$586 (AU$782) for each additional diagnosis. CONCLUSION: Early application of WES in Mendelian disorders is cost-effective and reanalysis of an undiagnosed individual at a 12-month time point increases total diagnoses by 11%.

A SLC39A8 variant causes manganese deficiency, and glycosylation and mitochondrial disorders.

SLC39A8 variants have recently been reported to cause a type II congenital disorder of glycosylation (CDG) in patients with intellectual disability and cerebellar atrophy. Here we report a novel SLC39A8 variant in siblings with features of Leigh-like mitochondrial disease. Two sisters born to consanguineous Lebanese parents had profound developmental delay, dystonia, seizures and failure to thrive. Brain MRI of both siblings identified bilateral basal ganglia hyperintensities on T2-weighted imaging and cerebral atrophy. CSF lactate was elevated in patient 1 and normal in patient 2. Respiratory chain enzymology was only performed on patient 1 and revealed complex IV and II + III activity was low in liver, with elevated complex I activity. Complex IV activity was borderline low in patient 1 muscle and pyruvate dehydrogenase activity was reduced. Whole genome sequencing identified a homozygous Chr4(GRCh37):g.103236869C>G; c.338G>C; p.(Cys113Ser) variant in SLC39A8, located in one of eight regions identified by homozygosity mapping. SLC39A8 encodes a manganese and zinc transporter which localises to the cell and mitochondrial membranes. Patient 2 blood and urine manganese levels were undetectably low. Transferrin electrophoresis of patient 2 serum revealed a type II CDG defect. Oral supplementation with galactose and uridine led to improvement of the transferrin isoform pattern within 14 days of treatment initiation. Oral manganese has only recently been added to the treatment. These results suggest SLC39A8 deficiency can cause both a type II CDG and Leigh-like syndrome, possibly via reduced activity of the manganese-dependent enzymes ß-galactosyltransferase and mitochondrial manganese superoxide dismutase.

Global Phosphoproteomic Mapping of Early Mitotic Exit in Human Cells Identifies Novel Substrate Dephosphorylation Motifs.

Entry into mitosis is driven by the coordinated phosphorylation of thousands of proteins. For the cell to complete mitosis and divide into two identical daughter cells it must regulate dephosphorylation of these proteins in a highly ordered, temporal manner. There is currently a lack of a complete understanding of the phosphorylation changes that occur during the initial stages of mitotic exit in human cells. Therefore, we performed a large unbiased, global analysis to map the very first dephosphorylation events that occur as cells exit mitosis. We identified and quantified the modification of >16,000 phosphosites on >3300 unique proteins during early mitotic exit, providing up to eightfold greater resolution than previous studies. The data have been deposited to the ProteomeXchange with identifier PXD001559. Only a small fraction (∼ 10%) of phosphorylation sites were dephosphorylated during early mitotic exit and these occurred on proteins involved in critical early exit events, including organization of the mitotic spindle, the spindle assembly checkpoint, and reformation of the nuclear envelope. Surprisingly this enrichment was observed across all kinase consensus motifs, indicating that it is independent of the upstream phosphorylating kinase. Therefore, dephosphorylation of these sites is likely determined by the specificity of phosphatase/s rather than the activity of kinase/s. Dephosphorylation was significantly affected by the amino acids at and surrounding the phosphorylation site, with several unique evolutionarily conserved amino acids correlating strongly with phosphorylation status. These data provide a potential mechanism for the specificity of phosphatases, and how they co-ordinate the ordered events of mitotic exit. In summary, our results provide a global overview of the phosphorylation changes that occur during the very first stages of mitotic exit, providing novel mechanistic insight into how phosphatase/s specifically regulate this critical transition.

Defining the genetic basis of early onset hereditary spastic paraplegia using whole genome sequencing.

We performed whole genome sequencing (WGS) in nine families from India with early-onset hereditary spastic paraplegia (HSP). We obtained a genetic diagnosis in 4/9 (44 %) families within known HSP genes (DDHD2 and CYP2U1), as well as perixosomal biogenesis disorders (PEX16) and GM1 gangliosidosis (GLB1). In the remaining patients, no candidate structural variants, copy number variants or predicted splice variants affecting an extended candidate gene list were identified. Our findings demonstrate the efficacy of using WGS for diagnosing early-onset HSP, particularly in consanguineous families (4/6 diagnosed), highlighting that two of the diagnoses would not have been made using a targeted approach.

Saccharomyces eubayanus and Saccharomyces arboricola reside in North Island native New Zealand forests.

Saccharomyces is one of the best-studied microbial genera, but our understanding of the global distributions and evolutionary histories of its members is relatively poor. Recent studies have altered our view of Saccharomyces' origin, but a lack of sampling from the vast majority of the world precludes a holistic perspective. We evaluate alternate Gondwanan and Far East Asian hypotheses concerning the origin of these yeasts. Being part of Gondwana, and only colonized by humans in the last ∼1000 years, New Zealand represents a unique environment for testing these ideas. Genotyping and ribosomal sequencing of samples from North Island native forest parks identified a widespread population of Saccharomyces. Whole genome sequencing identified the presence of S. arboricola and S. eubayanus in New Zealand, which is the first report of S. arboricola outside Far East Asia, and also expands S. eubayanus' known distribution to include the Oceanic region. Phylogenomic approaches place the S. arboricola population as significantly diverged from the only other sequenced Chinese isolate but indicate that S. eubayanus might be a recent migrant from South America. These data tend to support the Far East Asian origin of the Saccharomyces, but the history of this group is still far from clear.

European derived Saccharomyces cerevisiae colonisation of New Zealand vineyards aided by humans.

Humans have acted as vectors for species and expanded their ranges since at least the dawn of agriculture. While relatively well characterised for macrofauna and macroflora, the extent and dynamics of human-aided microbial dispersal is poorly described. We studied the role which humans have played in manipulating the distribution of Saccharomyces cerevisiae, one of the world's most important microbes, using whole genome sequencing. We include 52 strains representative of the diversity in New Zealand to the global set of genomes for this species. Phylogenomic approaches show an exclusively European origin of the New Zealand population, with a minimum of 10 founder events mostly taking place over the last 1000 years. Our results show that humans have expanded the range of S. cerevisiae and transported it to New Zealand where it was not previously present, where it has now become established in vineyards, but radiation to native forests appears limited.