A gene pathogenicity tool "GenePy" identifies missed biallelic diagnoses in the 100,000 Genomes Project.

PURPOSE: The 100,000 Genomes Project diagnosed a quarter of affected participants, but 26% of diagnoses were not on the applied gene panel(s); with many being de novo variants. Assessing biallelic variants without a gene panel is more challenging. METHODS: We sought to identify missed biallelic diagnoses using GenePy, which incorporates allele frequency, zygosity, and a user-defined deleterious metric, generating an aggregate GenePy score per gene, per participant. We calculated GenePy scores for 2862 recessive disease genes in 78,216 100,000 Genomes Project participants. For each gene, we ranked participant GenePy scores and scrutinized affected participants without a diagnosis, whose scores ranked among the top 5 for each gene. In cases which participant phenotypes overlapped with the disease gene of interest, we extracted rare variants and applied phase, ClinVar, and ACMG classification. RESULTS: 3184 affected individuals without a molecular diagnosis had a top-5-ranked GenePy score and 682 of 3184 (21%) had phenotypes overlapping with a top-ranking gene. In 122 of 669 (18%) phenotype-matched cases (excluding 13 withdrawn participants), we identified a putative missed diagnosis (2.2% of all undiagnosed participants). A further 334 of 669 (50%) cases have a possible missed diagnosis but require functional validation. CONCLUSION: Applying GenePy at scale has identified 456 potential diagnoses, demonstrating the value of novel diagnostic strategies.

Targeting de novo loss-of-function variants in constrained disease genes improves diagnostic rates in the 100,000 Genomes Project.

BACKGROUND: Genome sequencing was first offered clinically in the UK through the 100,000 Genomes Project (100KGP). Analysis was restricted to predefined gene panels associated with the patient's phenotype. However, panels rely on clearly characterised phenotypes and risk missing diagnoses outside of the panel(s) applied. We propose a complementary method to rapidly identify pathogenic variants, including those missed by 100KGP methods. METHODS: The Loss-of-function Observed/Expected Upper-bound Fraction (LOEUF) score quantifies gene constraint, with low scores correlated with haploinsufficiency. We applied DeNovoLOEUF, a filtering strategy to sequencing data from 13,949 rare disease trios in the 100KGP, by filtering for rare, de novo, loss-of-function variants in disease genes with a LOEUF score < 0.2. We compared our findings with the corresponding patient's diagnostic reports. RESULTS: 324/332 (98%) of the variants identified using DeNovoLOEUF were diagnostic or partially diagnostic (whereby the variant was responsible for some of the phenotype). We identified 39 diagnoses that were "missed" by 100KGP standard analyses, which are now being returned to patients. CONCLUSION: We have demonstrated a highly specific and rapid method with a 98% positive predictive value that has good concordance with standard analysis, low false-positive rate, and can identify additional diagnoses. Globally, as more patients are being offered genome sequencing, we anticipate that DeNovoLOEUF will rapidly identify new diagnoses and facilitate iterative analyses when new disease genes are discovered.

A novel variant in GATM causes idiopathic renal Fanconi syndrome and predicts progression to end-stage kidney disease.

Renal Fanconi syndrome (RFS) is a generalised disorder of the proximal convoluted tubule. Many genes have been associated with RFS including those that cause systemic disorders such as cystinosis, as well as isolated RFS. We discuss the case of a 10-year-old female who presented with leg pain and raised creatinine on a screening blood test. Her mother has RFS and required a kidney transplant in her thirties. Further investigations confirmed RFS in the daughter. Exome sequencing was performed on the affected mother, child, and unaffected father. We identified a novel variant in GATM; c.965G>C p.(Arg322Pro) segregating dominantly in the mother and daughter. We validated our finding with molecular dynamics simulations and demonstrated a dynamic signature that differentiates our variant and two previously identified pathogenic variants in GATM from wildtype. Genetic testing has uncovered a novel pathogenic variant that predicts progression to end stage kidney failure and has important implications for family planning and cascade testing. We recommend that GATM is screened for in children presenting with RFS, in addition to adults, particularly with kidney failure, who may have had previous negative gene testing.

The Pediatric Crohn Disease Morbidity Index (PCD-MI): Development of a Tool to Assess Long-Term Disease Burden Using a Data-Driven Approach.

BACKGROUND/OBJECTIVE: Heterogeneity and chronicity of Crohn disease (CD) make prediction of outcomes difficult. To date, no longitudinal measure can quantify burden over a patient's disease course, preventing assessment and integration into predictive modeling. Here, we aimed to demonstrate the feasibility of constructing a data driven, longitudinal disease burden score. METHODS: Literature was reviewed for tools used in assessment of CD activity. Themes were identified to construct a pediatric CD morbidity index (PCD-MI). Scores were assigned to variables. Data were extracted automatically from the electronic patient records at Southampton Children's Hospital, diagnosed from 2012 to 2019 (inclusive). PCD-MI scores were calculated, adjusted for duration of follow up and assessed for variation (ANOVA) and distribution (Kolmogorov-Smirnov). RESULTS: Nineteen clinical/biological features across five themes were included in the PCD-MI including blood/fecal/radiological/endoscopic results, medication usage, surgery, growth parameters, and extraintestinal manifestations. Maximal score was 100 after accounting for follow-up duration. PCD-MI was assessed in 66 patients, mean age 12.5 years. Following quality filtering, 9528 blood/fecal test results and 1309 growth measures were included. Mean PCD-MI score was 14.95 (range 2.2-32.5); data were normally distributed ( P = 0.2) with 25% of patients having a PCD-MI < 10. There was no difference in the mean PCD-MI when split by year of diagnosis, F -statistic 1.625, P = 0.147. CONCLUSIONS: PCD-MI is a calculatable measure for a cohort of patients diagnosed over an 8-year period, integrating a wide-range of data with potential to determine high or low disease burden. Future iterations of the PCD-MI require refinement of included features, optimized scores, and validation on external cohorts.

De novo putative loss-of-function variants in TAF4 are associated with a neuro-developmental disorder.

TATA-binding protein associated factor 4 (TAF4) is a subunit of the Transcription Factor IID (TFIID) complex, a central player in transcription initiation. Other members of this multimeric complex have been implicated previously as monogenic disease genes in human developmental disorders. TAF4 has not been described to date as a monogenic disease gene. We here present a cohort of eight individuals, each carrying de novo putative loss-of-function (pLoF) variants in TAF4 and expressing phenotypes consistent with a neuro-developmental disorder (NDD). Common features include intellectual disability, abnormal behavior, and facial dysmorphisms. We propose TAF4 as a novel dominant disease gene for NDD, and coin this novel disorder "TAF4-related NDD" (T4NDD). We place T4NDD in the context of other disorders related to TFIID subunits, revealing shared features of T4NDD with other TAF-opathies.

A gene-to-patient approach uplifts novel disease gene discovery and identifies 18 putative novel disease genes.

PURPOSE: Exome and genome sequencing have drastically accelerated novel disease gene discoveries. However, discovery is still hindered by myriad variants of uncertain significance found in genes of undetermined biological function. This necessitates intensive functional experiments on genes of equal predicted causality, leading to a major bottleneck. METHODS: We apply the loss-of-function observed/expected upper-bound fraction metric of intolerance to gene inactivation to curate a list of predicted haploinsufficient disease genes. Using data from the 100,000 Genomes Project, we adopt a gene-to-patient approach that matches de novo loss-of-function variants in constrained genes to patients with rare disease. Through large-scale aggregation of data, we reduce excess analytical noise currently hindering novel discoveries. RESULTS: Results from 13,949 trios revealed 643 rare, de novo predicted loss-of-function events filtered from 1044 loss-of-function observed/expected upper-bound fraction-constrained genes. A total of 168 variants occurred within 126 genes without a known disease-gene relationship. Of these, 27 genes had >1 kindred affected, and for 18 of these genes, multiple kindreds had overlapping phenotypes. Two years after initial analysis, 11 of 18 (61%) of these genes have been independently published as novel disease gene discoveries. CONCLUSION: Using large cohorts and adopting gene-based approaches can rapidly and objectively accelerate dominantly inherited novel gene discovery by targeting the most appropriate genes for functional validation.

Sequencing era methods for identifying signatures of selection in the genome.

Insights into genetic loci which are under selection and their functional roles contribute to increased understanding of the patterns of phenotypic variation we observe today. The availability of whole-genome sequence data, for humans and other species, provides opportunities to investigate adaptation and evolution at unprecedented resolution. Many analytical methods have been developed to interrogate these large data sets and characterize signatures of selection in the genome. We review here recently developed methods and consider the impact of increased computing power and data availability on the detection of selection signatures. Consideration of demography, recombination and other confounding factors is important, and use of a range of methods in combination is a powerful route to resolving different forms of selection in genome sequence data. Overall, a substantial improvement in methods for application to whole-genome sequencing is evident, although further work is required to develop robust and computationally efficient approaches which may increase reproducibility across studies.

Growth failure is rare in a contemporary cohort of paediatric inflammatory bowel disease patients.

AIM: We assessed growth in a paediatric inflammatory bowel disease (PIBD) cohort. METHODS: Paediatric inflammatory bowel disease patients were eligible if they were diagnosed at Southampton Children's Hospital from 2011 to 2018. Weight and height standard deviation scores (SDS) were retrieved. Mean SDS values, SDS change and anti-TNF status were analysed at diagnosis and during follow-up. RESULTS: Four hundred and ninety patients were included, 313 with Crohn's disease (CD). CD patients presented with mean height SDS -0.13, -0.1 at 1-year, -0.11 at 2-years and -0.03 at 5 years, reflecting preserved linear growth. There was no significant height-SDS change from diagnosis to 5-year follow-up, +0.12, 95%-CI: 0.48 to -0.24. Mean weight-SDS at diagnosis was -0.39, driven by CD patients (-0.65). Mean weight-SDS approached 0 after 1 year and remained at the 50th centile throughout follow-up. Growth in ulcerative colitis was maintained. In multivariable regression males had worse height growth from diagnosis to transition (P = .036). Anti-TNF treatment (P = .013) and surgical resection (P = .005) were also associated with poorer linear growth. Patients treated with anti-TNF therapy had lower height-SDS compared to those never treated with anti-TNF at 1 year (-0.2 vs -0.01, P = .22), 2-years (-0.27 vs -0.01, P = .07) and 5 years (-0.21 vs 0.25, P = .051). CONCLUSION: Height was generally maintained in Crohn's disease, and impaired linear growth was rare in this cohort.

Immunological Profiling of Paediatric Inflammatory Bowel Disease Using Unsupervised Machine Learning.

OBJECTIVES: The current classification of inflammatory bowel disease (IBD) is based on clinical phenotypes, which is blind to the molecular basis of the disease. The aim of this study was to stratify a treatment-naïve paediatric IBD cohort through specific innate immunity pathway profiling and application of unsupervised machine learning (UML). METHODS: In order to test the molecular integrity of biological pathways implicated in IBD, innate immune responses were assessed at diagnosis in 22 paediatric patients and 10 age-matched controls. Peripheral blood mononuclear cells (PBMCs) were selectively stimulated for assessing the functionality of upstream activation receptors including NOD2, toll-like receptor (TLR) 1-2 and TLR4, and the downstream cytokine responses (IL-10, IL-1ß, IL-6, and TNF-α) using multiplex assays. Cytokine data generated were subjected to hierarchical clustering to assess for patient stratification. RESULTS: Combined immune responses in patients across 12 effector responses were significantly reduced compared with controls (P = 0.003) and driven primarily by "hypofunctional" TLR responses (P values 0.045, 0.010, and 0.018 for TLR4-mediated IL-10, IL-1ß, and TNF-α, respectively; 0.018 and 0.015 for TLR1-2 -mediated IL-10 and IL-1ß). Hierarchical clustering generated 3 distinct clusters of patients and a fourth group of "unclustered" individuals. No relationship was observed between the observed immune clusters and the clinical disease phenotype. CONCLUSIONS: Although a clinically useful outcome was not observed through hierarchical clustering, our study provides a rationale for using an UML approach to stratify patients. The study also highlights the predominance of hypo-inflammatory innate immune responses as a key mechanism in the pathogenesis of IBD.

Inactivation of AMMECR1 is associated with growth, bone, and heart alterations.

We report five individuals with loss-of-function of the X-linked AMMECR1: a girl with a balanced X-autosome translocation and inactivation of the normal X-chromosome; two boys with maternally inherited and de novo nonsense variants; and two half-brothers with maternally inherited microdeletion variants. They present with short stature, cardiac and skeletal abnormalities, and hearing loss. Variants of unknown significance in AMMECR1 in four male patients from two families with partially overlapping phenotypes were previously reported. AMMECR1 is coexpressed with genes implicated in cell cycle regulation, five of which were previously associated with growth and bone alterations. Our knockdown of the zebrafish orthologous gene resulted in phenotypes reminiscent of patients' features. The increased transcript and encoded protein levels of AMMECR1L, an AMMECR1 paralog, in the t(X;9) patient's cells indicate a possible partial compensatory mechanism. AMMECR1 and AMMECR1L proteins dimerize and localize to the nucleus as suggested by their nucleic acid-binding RAGNYA folds. Our results suggest that AMMECR1 is potentially involved in cell cycle control and linked to a new syndrome with growth, bone, heart, and kidney alterations with or without elliptocytosis.

Endoscopic and Histological Assessment of Paediatric Inflammatory Bowel Disease Over a 3-Year Follow-up Period.

OBJECTIVES: Discrepancies between inflammatory bowel disease (IBD) endoscopic/histological extent are documented at diagnosis. It is unclear whether these differences persist through disease course, with potential impact on categorization and management. We aimed to analyze the progression of disease over a 3-year period. METHODS: Patients younger than 17 years, diagnosed between 2010 and 2013 at Southampton Children's Hospital and followed-up for 3 years were eligible. Primary outcome was disease extent at diagnosis and follow-up. Data are presented as percentage of patients undergoing endoscopy. Paris classification (PC) and PC using histological, rather than endoscopic disease, were determined. RESULTS: One hundred and twenty-five patients were included, 66 boys; Crohn's disease (CD) 74, ulcerative colitis (UC) 40, IBD unclassified (IBDU) 11. All had endoscopy at diagnosis. One hundred and two patients underwent ≥1 repeat endoscopies.Disease extent reduced from diagnosis to first follow-up endoscopy for both endoscopic and histological disease extent (CD/UC/IBDU, all P < 0.00006). Histological extent remained greater than endoscopic in CD with significant differences in stomach, ileum, and large bowel at all follow-up points (P =  < 0.045). Endoscopic matched histological extent in UC/IBDU. Applying a modified PC resulted in significant changes for CD (L3 27.4%-53.2%, P = 0.006, L3 + L4A 21%-50%, P = 0.001, and upper gastrointestinal disease 50%-80.6%, P = 0.0006) but not UC. CD height (-0.37 to -0.25) and weight (-1.09 to -0.19) standard deviation scores increased from diagnosis to follow-up. CONCLUSIONS: Histological disease is greater than endoscopic extent at diagnosis and during follow-up in CD, although not in UC/IBDU. Classification of disease extent in CD should be based on both endoscopic and histological criteria.

AMMECR1: a single point mutation causes developmental delay, midface hypoplasia and elliptocytosis.

BACKGROUND: Deletions in the Xq22.3-Xq23 region, inclusive of COL4A5, have been associated with a contiguous gene deletion syndrome characterised by Alport syndrome with intellectual disability (Mental retardation), Midface hypoplasia and Elliptocytosis (AMME). The extrarenal biological and clinical significance of neighbouring genes to the Alport locus has been largely speculative. We sought to discover a genetic cause for two half-brothers presenting with nephrocalcinosis, early speech and language delay and midface hypoplasia with submucous cleft palate and bifid uvula. METHODS: Whole exome sequencing was undertaken on maternal half-siblings. In-house genomic analysis included extraction of all shared variants on the X chromosome in keeping with X-linked inheritance. Patient-specific mutants were transfected into three cell lines and microscopically visualised to assess the nuclear expression pattern of the mutant protein. RESULTS: In the affected half-brothers, we identified a hemizygous novel non-synonymous variant of unknown significance in AMMECR1 (c.G530A; p.G177D), a gene residing in the AMME disease locus. Transfected cell lines with the p.G177D mutation showed aberrant nuclear localisation patterns when compared with the wild type. Blood films revealed the presence of elliptocytes in the older brother. CONCLUSIONS: Our study shows that a single missense mutation in AMMECR1 causes a phenotype of midface hypoplasia, mild intellectual disability and the presence of elliptocytes, previously reported as part of a contiguous gene deletion syndrome. Functional analysis confirms mutant-specific protein dysfunction. We conclude that AMMECR1 is a critical gene in the pathogenesis of AMME, causing midface hypoplasia and elliptocytosis and contributing to early speech and language delay, infantile hypotonia and hearing loss, and may play a role in dysmorphism, nephrocalcinosis and submucous cleft palate.

Autosomal dominant tubulointerstitial kidney disease-UMOD is the most frequent non polycystic genetic kidney disease.

BACKGROUND: Autosomal dominant tubulointerstitial kidney disease (ADTKD) caused by mutations in the UMOD gene (ADTKD-UMOD) is considered rare and often remains unrecognised. We aimed to establish the prevalence of genetic kidney diseases, ADTKD and ADTKD-UMOD in adult chronic kidney disease (CKD) patients, and to investigate characteristic features. METHODS: We sent questionnaires on family history to all patients with CKD stages 3-5 in our tertiary renal centre to identify patients with inherited renal disease. Details on clinical and family history were obtained from patient interviews and clinical records. Sanger sequencing of the UMOD gene was performed from blood or saliva samples. RESULTS: 2027 of 3770 sent questionnaires were returned. 459 patients reported a family history, which was consistent with inherited kidney disease in 217 patients. 182 non-responders with inherited kidney diseases were identified through a database search. Of these 399 individuals, 252 had autosomal dominant polycystic kidney disease (ADPKD), 28 had ADTKD, 25 had Alports, and 44 were unknown, resulting in 11% of CKD 3-5 patients and 19% of end-stage renal disease patients with genetic kidney diseases. Of the unknown, 40 were genotyped, of whom 31 had findings consistent with ADTKD. 30% of unknowns and 39% of unknowns with ADTKD had UMOD mutations. Altogether, 35 individuals from 18 families were found to have ten distinct UMOD mutations (three novel), making up 1% of patients with CKD 3-5, 2% of patients with end-stage renal disease, 9% of inherited kidney diseases and 56% with ADTKD. ADTKD-UMOD was the most common genetic kidney disease after ADPKD with a population prevalence of 9 per million. Less proteinuria and haematuria, but not hyperuricaemia or gout were predictive of ADTKD-UMOD. The main limitations of the study are the single-centre design and a predominantly Caucasian population. CONCLUSIONS: The prevalence of genetic kidney diseases and ADTKD-UMOD is significantly higher than previously described. Clinical features poorly predicted ADTKD-UMOD, highlighting the need for genetic testing guided by family history alone.

ADCY5-related dyskinesia presenting as familial myoclonus-dystonia.

We describe a family with an autosomal dominant familial dyskinesia resembling myoclonus-dystonia associated with a novel missense mutation in ADCY5, found through whole-exome sequencing. A tiered analytical approach was used to analyse whole-exome sequencing data from an affected grandmother-granddaughter pair. Whole-exome sequencing identified 18,000 shared variants, of which 46 were non-synonymous changes not present in a local cohort of control exomes (n = 422). Further filtering based on predicted splicing effect, minor allele frequency in the 1000 Genomes Project and on phylogenetic conservation yielded 13 candidate variants, of which the heterozygous missense mutation c.3086T>G, p. M1029R in ADCY5 most closely matched the observed phenotype. This report illustrates the utility of whole-exome sequencing in cases of undiagnosed movement disorders with clear autosomal dominant inheritance. Moreover, ADCY5 mutations should be considered in cases with apparent myoclonus-dystonia, particularly where SCGE mutations have been excluded. ADCY5-related dyskinesia may manifest variable expressivity within a single family, and affected individuals may be initially diagnosed with differing neurological phenotypes.

Exome sequence read depth methods for identifying copy number changes.

Copy number variants (CNVs) play important roles in a number of human diseases and in pharmacogenetics. Powerful methods exist for CNV detection in whole genome sequencing (WGS) data, but such data are costly to obtain. Many disease causal CNVs span or are found in genome coding regions (exons), which makes CNV detection using whole exome sequencing (WES) data attractive. If reliably validated against WGS-based CNVs, exome-derived CNVs have potential applications in a clinical setting. Several algorithms have been developed to exploit exome data for CNV detection and comparisons made to find the most suitable methods for particular data samples. The results are not consistent across studies. Here, we review some of the exome CNV detection methods based on depth of coverage profiles and examine their performance to identify problems contributing to discrepancies in published results. We also present a streamlined strategy that uses a single metric, the likelihood ratio, to compare exome methods, and we demonstrated its utility using the VarScan 2 and eXome Hidden Markov Model (XHMM) programs using paired normal and tumour exome data from chronic lymphocytic leukaemia patients. We use array-based somatic CNV (SCNV) calls as a reference standard to compute prevalence-independent statistics, such as sensitivity, specificity and likelihood ratio, for validation of the exome-derived SCNVs. We also account for factors known to influence the performance of exome read depth methods, such as CNV size and frequency, while comparing our findings with published results.

Mutations specific to the Rac-GEF domain of TRIO cause intellectual disability and microcephaly.

BACKGROUND: Neurodevelopmental disorders have challenged clinical genetics for decades, with over 700 genes implicated and many whose function remains unknown. The application of whole-exome sequencing is proving pivotal in closing the genotype/phenotype gap through the discovery of new genes and variants that help to unravel the pathogenic mechanisms driving neuropathogenesis. One such discovery includes TRIO, a gene recently implicated in neurodevelopmental delay. Trio is a Dbl family guanine nucleotide exchange factor (GEF) and a major regulator of neuronal development, controlling actin cytoskeleton dynamics by activating the GTPase Rac1. METHODS: Whole-exome sequencing was undertaken on a family presenting with global developmental delay, microcephaly and mild dysmorphism. Father/daughter exome analysis was performed, followed by confirmatory Sanger sequencing and segregation analysis on four individuals. Three further patients were recruited through the deciphering developmental disorders (DDD) study. Functional studies were undertaken using patient-specific Trio protein mutations. RESULTS: We identified a frameshift deletion in TRIO that segregated autosomal dominantly. By scrutinising data from DDD, we further identified three unrelated children with a similar phenotype who harboured de novo missense mutations in TRIO. Biochemical studies demonstrated that in three out of four families, the Trio mutations led to a markedly reduced Rac1 activation. CONCLUSIONS: We describe an inherited global developmental delay phenotype associated with a frameshift deletion in TRIO. Additionally, we identify pathogenic de novo missense mutations in TRIO associated with the same consistent phenotype, intellectual disability, microcephaly and dysmorphism with striking digital features. We further functionally validate the importance of the GEF domain in Trio protein function. Our study demonstrates how genomic technologies are yet again proving prolific in diagnosing and advancing the understanding of neurodevelopmental disorders.

Collagen (COL4A) mutations are the most frequent mutations underlying adult focal segmental glomerulosclerosis.

BACKGROUND: Multiple genes underlying focal segmental glomerulosclerosis (FSGS) and/or steroid-resistant nephrotic syndrome (SRNS) have been identified, with the recent inclusion of collagen IV mutations responsible for Alport disease (AD) or thin basement membrane nephropathy (TBMN). We aimed to investigate the distribution of gene mutations in adult patients with primary FSGS/SRNS by targeted next generation sequencing (NGS). METHODS: Eighty-one adults from 76 families were recruited; 24 families had a history of renal disease. A targeted NGS panel was designed and applied, covering 39 genes implicated in FSGS/SRNS including COL4A3-5. RESULTS: Confirmed pathogenic mutations were found in 10 patients (6 with family history) from 9 families (diagnostic rate 12%). Probably pathogenic mutations were identified in an additional six patients (combined diagnostic rate 20%). Definitely pathogenic mutations were identified in 22% of patients with family history and 10% without. Mutations in COL4A3-5 were present in eight patients from six families, representing 56% of definitely pathogenic mutations, and establishing a diagnosis of AD in six patients and TBMN in two patients. Collagen mutations were identified in 38% of families with familial FSGS, and 3% with sporadic FSGS, with over half the mutations occurring in COL4A5. Patients with collagen mutations were younger at presentation and more likely to have family history, haematuria and glomerular basement membrane abnormalities. CONCLUSIONS: We show that collagen IV mutations, including COL4A5, frequently underlie FSGS and should be considered, particularly with a positive family history. Targeted NGS improves diagnostic efficiency by investigating many candidate genes in parallel.