100,000 Genomes Pilot on Rare-Disease Diagnosis in Health Care - Preliminary Report

BACKGROUND

The U.K. 100,000 Genomes Project is in the process of investigating the role of genome sequencing in patients with undiagnosed rare diseases after usual care and the alignment of this research with health care implementation in the U.K. National Health Service. Other parts of this project focus on patients with cancer and infection.

METHODS

We conducted a pilot study involving 4660 participants from 2183 families, among whom 161 disorders covering a broad spectrum of rare diseases were present. We collected data on clinical features with the use of Human Phenotype Ontology terms, undertook genome sequencing, applied automated variant prioritization on the basis of applied virtual gene panels and phenotypes, and identified novel pathogenic variants through research analysis.

RESULTS

Diagnostic yields varied among family structures and were highest in family trios (both parents and a proband) and families with larger pedigrees. Diagnostic yields were much higher for disorders likely to have a monogenic cause (35%) than for disorders likely to have a complex cause (11%). Diagnostic yields for intellectual disability, hearing disorders, and vision disorders ranged from 40 to 55%. We made genetic diagnoses in 25% of the probands. A total of 14% of the diagnoses were made by means of the combination of research and automated approaches, which was critical for cases in which we found etiologic noncoding, structural, and mitochondrial genome variants and coding variants poorly covered by exome sequencing. Cohortwide burden testing across 57,000 genomes enabled the discovery of three new disease genes and 19 new associations. Of the genetic diagnoses that we made, 25% had immediate ramifications for clinical decision making for the patients or their relatives.

CONCLUSIONS

Our pilot study of genome sequencing in a national health care system showed an increase in diagnostic yield across a range of rare diseases. (Funded by the National Institute for Health Research and others.).

Heterozygous ANKRD17 loss-of-function variants cause a syndrome with intellectual disability, speech delay, and dysmorphism

ANKRD17 is an ankyrin repeat-containing protein thought to play a role in cell cycle progression, whose ortholog in Drosophila functions in the Hippo pathway as a co-factor of Yorkie. Here, we delineate a neurodevelopmental disorder caused by de novo heterozygous ANKRD17 variants. The mutational spectrum of this cohort of 34 individuals from 32 families is highly suggestive of haploinsufficiency as the underlying mechanism of disease, with 21 truncating or essential splice site variants, 9 missense variants, 1 in-frame insertion-deletion, and 1 microdeletion (1.16 Mb). Consequently, our data indicate that loss of ANKRD17 is likely the main cause of phenotypes previously associated with large multi-gene chromosomal aberrations of the 4q13.3 region. Protein modeling suggests that most of the missense variants disrupt the stability of the ankyrin repeats through alteration of core structural residues. The major phenotypic characteristic of our cohort is a variable degree of developmental delay/intellectual disability, particularly affecting speech, while additional features include growth failure, feeding difficulties, non-specific MRI abnormalities, epilepsy and/or abnormal EEG, predisposition to recurrent infections (mostly bacterial), ophthalmological abnormalities, gait/balance disturbance, and joint hypermobility. Moreover, many individuals shared similar dysmorphic facial features. Analysis of single-cell RNA-seq data from the developing human telencephalon indicated ANKRD17 expression at multiple stages of neurogenesis, adding further evidence to the assertion that damaging ANKRD17 variants cause a neurodevelopmental disorder.

Leukoencephalopathy with calcifications and cysts: Genetic and phenotypic spectrum

Biallelic mutations in SNORD118, encoding the small nucleolar RNA U8, cause leukoencephalopathy with calcifications and cysts (LCC). Given the difficulty in interpreting the functional consequences of variants in nonprotein encoding genes, and the high allelic polymorphism across SNORD118 in controls, we set out to provide a description of the molecular pathology and clinical spectrum observed in a cohort of patients with LCC. We identified 64 affected individuals from 56 families. Age at presentation varied from 3 weeks to 67 years, with disease onset after age 40 years in eight patients. Ten patients had died. We recorded 44 distinct, likely pathogenic, variants in SNORD118. Fifty two of 56 probands were compound heterozygotes, with parental consanguinity reported in only three families. Forty nine of 56 probands were either heterozygous (46) or homozygous (three) for a mutation involving one of seven nucleotides that facilitate a novel intramolecular interaction between the 5' end and 3' extension of precursor-U8. There was no obvious genotype-phenotype correlation to explain the marked variability in age at onset. Complementing recently published functional analyses in a zebrafish model, these data suggest that LCC most often occurs due to combinatorial severe and milder mutations, with the latter mostly affecting 3' end processing of precursor-U8.

Delineation of dominant and recessive forms of LZTR1-associated Noonan syndrome

Noonan syndrome (NS) is characterised by distinctive facial features, heart defects, variable degrees of intellectual disability and other phenotypic manifestations. Although the mode of inheritance is typically dominant, recent studies indicate LZTR1 may be associated with both dominant and recessive forms. Seeking to describe the phenotypic characteristics of LZTR1-associated NS, we searched for likely pathogenic variants using two approaches. First, scrutiny of exomes from 9624 patients recruited by the Deciphering Developmental Disorders (DDDs) study uncovered six dominantly-acting mutations (p.R97L; p.Y136C; p.Y136H, p.N145I, p.S244C; p.G248R) of which five arose de novo, and three patients with compound-heterozygous variants (p.R210*/p.V579M; p.R210*/p.D531N; c.1149+1G>T/p.R688C). One patient also had biallelic loss-of-function mutations in NEB, consistent with a composite phenotype. After removing this complex case, analysis of human phenotype ontology terms indicated significant phenotypic similarities (P = 0.0005), supporting a causal role for LZTR1. Second, targeted sequencing of eight unsolved NS-like cases identified biallelic LZTR1 variants in three further subjects (p.W469*/p.Y749C, p.W437*/c.-38T>A and p.A461D/p.I462T). Our study strengthens the association of LZTR1 with NS, with de novo mutations clustering around the KT1-4 domains. Although LZTR1 variants explain ~0.1% of cases across the DDD cohort, the gene is a relatively common cause of unsolved NS cases where recessive inheritance is suspected.

De Novo and Inherited Loss-of-Function Variants in TLK2: Clinical and Genotype-Phenotype Evaluation of a Distinct Neurodevelopmental Disorder

Next-generation sequencing is a powerful tool for the discovery of genes related to neurodevelopmental disorders (NDDs). Here, we report the identification of a distinct syndrome due to de novo or inherited heterozygous mutations in Tousled-like kinase 2 (TLK2) in 38 unrelated individuals and two affected mothers, using whole-exome and whole-genome sequencing technologies, matchmaker databases, and international collaborations. Affected individuals had a consistent phenotype, characterized by mild-borderline neurodevelopmental delay (86%), behavioral disorders (68%), severe gastro-intestinal problems (63%), and facial dysmorphism including blepharophimosis (82%), telecanthus (74%), prominent nasal bridge (68%), broad nasal tip (66%), thin vermilion of the upper lip (62%), and upslanting palpebral fissures (55%). Analysis of cell lines from three affected individuals showed that mutations act through a loss-of-function mechanism in at least two case subjects. Genotype-phenotype analysis and comparison of computationally modeled faces showed that phenotypes of these and other individuals with loss-of-function variants significantly overlapped with phenotypes of individuals with other variant types (missense and C-terminal truncating). This suggests that haploinsufficiency of TLK2 is the most likely underlying disease mechanism, leading to a consistent neurodevelopmental phenotype. This work illustrates the power of international data sharing, by the identification of 40 individuals from 26 different centers in 7 different countries, allowing the identification, clinical delineation, and genotype-phenotype evaluation of a distinct NDD caused by mutations in TLK2.

Short interspersed element (SINE) depletion and long interspersed element (LINE) abundance are not features universally required for imprinting

Genomic imprinting is a form of gene dosage regulation in which a gene is expressed from only one of the alleles, in a manner dependent on the parent of origin. The mechanisms governing imprinted gene expression have been investigated in detail and have greatly contributed to our understanding of genome regulation in general. Both DNA sequence features, such as CpG islands, and epigenetic features, such as DNA methylation and non-coding RNAs, play important roles in achieving imprinted expression. However, the relative importance of these factors varies depending on the locus in question. Defining the minimal features that are absolutely required for imprinting would help us to understand how imprinting has evolved mechanistically. Imprinted retrogenes are a subset of imprinted loci that are relatively simple in their genomic organisation, being distinct from large imprinting clusters, and have the potential to be used as tools to address this question. Here, we compare the repeat element content of imprinted retrogene loci with non-imprinted controls that have a similar locus organisation. We observe no significant differences that are conserved between mouse and human, suggesting that the paucity of SINEs and relative abundance of LINEs at imprinted loci reported by others is not a sequence feature universally required for imprinting.