Calculating the statistical significance of rare variants causal for Mendelian and complex disorders

STR mutations on chromosome 15q cause thyrotropin resistance by activating a primate-specific enhancer of MIR7-2/MIR1179

Thyrotropin (TSH) is the master regulator of thyroid gland growth and function. Resistance to TSH (RTSH) describes conditions with reduced sensitivity to TSH. Dominantly inherited RTSH has been linked to a locus on chromosome 15q, but its genetic basis has remained elusive. Here we show that non-coding mutations in a (TTTG)4 short tandem repeat (STR) underlie dominantly inherited RTSH in all 82 affected participants from 12 unrelated families. The STR is contained in a primate-specific Alu retrotransposon with thyroid-specific cis-regulatory chromatin features. Fiber-seq and RNA-seq studies revealed that the mutant STR activates a thyroid-specific enhancer cluster, leading to haplotype-specific upregulation of the bicistronic MIR7-2/MIR1179 locus 35 kb downstream and overexpression of its microRNA products in the participants' thyrocytes. An imbalance in signaling pathways targeted by these micro-RNAs provides a working model for this cause of RTSH. This finding broadens our current knowledge of genetic defects altering pituitary-thyroid feedback regulation.

Idiopathic erythrocytosis: a germline disease?

Polycythemia Vera (PV) is typically caused by V617F or exon 12 JAK2 mutations. Little is known about Polycythemia cases where no JAK2 variants can be detected, and no other causes identified. This condition is defined as idiopathic erythrocytosis (IE). We evaluated clinical-laboratory parameters of a cohort of 56 IE patients and we determined their molecular profile at diagnosis with paired blood/buccal-DNA exome-sequencing coupled with a high-depth targeted OncoPanel to identify a possible underling germline or somatic cause. We demonstrated that most of our cohort (40/56: 71.4%) showed no evidence of clonal hematopoiesis, suggesting that IE is, in large part, a germline disorder. We identified 20 low mutation burden somatic variants (Variant allelic fraction, VAF, < 10%) in only 14 (25%) patients, principally involving DNMT3A and TET2. Only 2 patients presented high mutation burden somatic variants, involving DNMT3A, TET2, ASXL1 and WT1. We identified recurrent germline variants in 42 (75%) patients occurring mainly in JAK/STAT, Hypoxia and Iron metabolism pathways, among them: JAK3-V722I and HIF1A-P582S; a high fraction of patients (48.2%) resulted also mutated in homeostatic iron regulatory gene HFE-H63D or C282Y. By generating cellular models, we showed that JAK3-V722I causes activation of the JAK-STAT5 axis and upregulation of EPAS1/HIF2A, while HIF1A-P582S causes suppression of hepcidin mRNA synthesis, suggesting a major role for these variants in the onset of IE.

A homozygous truncating mutation of FGL2 is associated with immune dysregulation

BACKGROUND

The type II transmembrane protein fibrinogen-like protein 2 (FGL2) plays critical roles in hemostasis and immune regulation. The C-terminal immunoregulatory domain of FGL2 can be secreted and is a mediator of regulatory T (Treg) cell suppression. Fgl2^-/- mice develop autoantibodies and glomerulonephritis and have impaired Treg cell function.

OBJECTIVE

Our aim was to identify the genetic underpinning and immune function in a patient with childhood onset of leukocytoclastic vasculitis, systemic inflammation, and autoantibodies.

METHODS

Whole-exome sequencing was performed on patient genomic DNA. FGL2 protein expression was examined in HEK293 transfected cells by immunoblotting and in PBMCs by flow cytometry. T follicular helper cells and Treg cells were examined by flow cytometry. Treg cell suppression of T-cell proliferation was assessed in vitro.

RESULTS

The patient had a homozygous mutation in FGL2 (c.614_617del:p.V205fs), which led to the expression of a truncated FGL2 protein that preserves the N-terminal domain but lacks the C-terminal immunoregulatory domain. The patient had an increased percentage of circulating T follicular helper and Treg cells. The patient's Treg cells had impaired in vitro suppressive ability that was rescued by the addition of full-length FGL2. Unlike full-length FGL2, the truncated FGL2^V205fs mutant failed to suppress T-cell proliferation.

CONCLUSIONS

We identified a homozygous mutation in FGL2 in a patient with immune dysregulation and impaired Treg cell function. Soluble FGL2 rescued the Treg cell defect, suggesting that it may provide a useful therapy for the patient.

Family-Based Whole-Exome Analysis of Specific Language Impairment (SLI) Identifies Rare Variants in BUD13, a Component of the Retention and Splicing (RES) Complex

Specific language impairment (SLI) is a common neurodevelopmental disorder (NDD) that displays high heritability estimates. Genetic studies have identified several loci, but the molecular basis of SLI remains unclear. With the aim to better understand the genetic architecture of SLI, we performed whole-exome sequencing (WES) in a single family (ID: 489; n = 11). We identified co-segregating rare variants in three new genes: BUD13, APLP2, and NDRG2. To determine the significance of these genes in SLI, we Sanger sequenced all coding regions of each gene in unrelated individuals with SLI (n = 175). We observed 13 additional rare variants in 18 unrelated individuals. Variants in BUD13 reached genome-wide significance (p-value < 0.01) upon comparison with similar variants in the 1000 Genomes Project, providing gene level evidence that BUD13 is involved in SLI. Additionally, five BUD13 variants showed cohesive variant level evidence of likely pathogenicity. Bud13 is a component of the retention and splicing (RES) complex. Additional supportive evidence from studies of an animal model (loss-of-function mutations in BUD13 caused a profound neural phenotype) and individuals with an NDD phenotype (carrying a CNV spanning BUD13), indicates BUD13 could be a target for investigation of the neural basis of language.

Genome sequencing data analysis for rare disease gene discovery

Rare diseases occur in a smaller proportion of the general population, which is variedly defined as less than 200 000 individuals (US) or in less than 1 in 2000 individuals (Europe). Although rare, they collectively make up to approximately 7000 different disorders, with majority having a genetic origin, and affect roughly 300 million people globally. Most of the patients and their families undergo a long and frustrating diagnostic odyssey. However, advances in the field of genomics have started to facilitate the process of diagnosis, though it is hindered by the difficulty in genome data analysis and interpretation. A major impediment in diagnosis is in the understanding of the diverse approaches, tools and datasets available for variant prioritization, the most important step in the analysis of millions of variants to select a few potential variants. Here we present a review of the latest methodological developments and spectrum of tools available for rare disease genetic variant discovery and recommend appropriate data interpretation methods for variant prioritization. We have categorized the resources based on various steps of the variant interpretation workflow, starting from data processing, variant calling, annotation, filtration and finally prioritization, with a special emphasis on the last two steps. The methods discussed here pertain to elucidating the genetic basis of disease in individual patient cases via trio- or family-based analysis of the genome data. We advocate the use of a combination of tools and datasets and to follow multiple iterative approaches to elucidate the potential causative variant.

Utility of Gene Panels for the Diagnosis of Inborn Errors of Metabolism in a Metabolic Reference Center

Next-generation sequencing (NGS) technologies have been proposed as a first-line test for the diagnosis of inborn errors of metabolism (IEM), a group of genetically heterogeneous disorders with overlapping or nonspecific phenotypes. Over a 3-year period, we prospectively analyzed 311 pediatric patients with a suspected IEM using four targeted gene panels. The rate of positive diagnosis was 61.86% for intermediary metabolism defects, 32.84% for complex molecular defects, 19% for hypoglycemic/hyperglycemic events, and 17% for mitochondrial diseases, and a conclusive molecular diagnosis was established in 2-4 weeks. Forty-one patients for whom negative results were obtained with the mitochondrial diseases panel underwent subsequent analyses using the NeuroSeq panel, which groups all genes from the individual panels together with genes associated with neurological disorders (1870 genes in total). This achieved a diagnostic rate of 32%. We next evaluated the utility of a tool, Phenomizer, for differential diagnosis, and established a correlation between phenotype and molecular findings in 39.3% of patients. Finally, we evaluated the mutational architecture of the genes analyzed by determining z-scores, loss-of-function observed/expected upper bound fraction (LOEUF), and haploinsufficiency (HI) scores. In summary, targeted gene panels for specific groups of IEMs enabled rapid and effective diagnosis, which is critical for the therapeutic management of IEM patients.

Nonsense mutations in alpha-II spectrin in three families with juvenile onset hereditary motor neuropathy

Distal hereditary motor neuropathies are a rare subgroup of inherited peripheral neuropathies hallmarked by a length-dependent axonal degeneration of lower motor neurons without significant involvement of sensory neurons. We identified patients with heterozygous nonsense mutations in the αII-spectrin gene, SPTAN1, in three separate dominant hereditary motor neuropathy families via next-generation sequencing. Variable penetrance was noted for these mutations in two of three families, and phenotype severity differs greatly between patients. The mutant mRNA containing nonsense mutations is broken down by nonsense-mediated decay and leads to reduced protein levels in patient cells. Previously, dominant-negative αII-spectrin gene mutations were described as causal in a spectrum of epilepsy phenotypes.

Targeted resequencing identifies genes with recurrent variation in cerebral palsy

A growing body of evidence points to a considerable and heterogeneous genetic aetiology of cerebral palsy (CP). To identify recurrently variant CP genes, we designed a custom gene panel of 112 candidate genes. We tested 366 clinically unselected singleton cases with CP, including 271 cases not previously examined using next-generation sequencing technologies. Overall, 5.2% of the naïve cases (14/271) harboured a genetic variant of clinical significance in a known disease gene, with a further 4.8% of individuals (13/271) having a variant in a candidate gene classified as intolerant to variation. In the aggregate cohort of individuals from this study and our previous genomic investigations, six recurrently hit genes contributed at least 4% of disease burden to CP: COL4A1, TUBA1A, AGAP1, L1CAM, MAOB and KIF1A. Significance of Rare VAriants (SORVA) burden analysis identified four genes with a genome-wide significant burden of variants, AGAP1, ERLIN1, ZDHHC9 and PROC, of which we functionally assessed AGAP1 using a zebrafish model. Our investigations reinforce that CP is a heterogeneous neurodevelopmental disorder with known as well as novel genetic determinants.

Confirmation that variants in TTI2 are responsible for autosomal recessive intellectual disability

TTI2 (MIM 614126) has been described as responsible for autosomal recessive intellectual disability (ID; MRT39, MIM:615541) in only two inbred families. Here, we give an account of two individuals from two unrelated outbred families harbouring compound heterozygous TTI2 pathogenic variants. Together with severe ID, progressive microcephaly, scoliosis and sleeping disorder are the most striking features in the two individuals concerned. TTI2, together with TTI1 and TELO2, encode proteins that constitute the triple T heterotrimeric complex. This TTT complex interacts with the HSP90 and R2TP to form a super-complex that has a chaperone function stabilising and maturing a number of kinases, such as ataxia-telangiectasia mutated and mechanistic target of rapamycin, which are key regulators of cell proliferation and genome maintenance. Pathogenic variants in TTI2 logically result in a phenotype close to that caused by TELO2 variants.

Identification of rare-disease genes using blood transcriptome sequencing and large control cohorts

It is estimated that 350 million individuals worldwide suffer from rare diseases, which are predominantly caused by mutation in a single gene1. The current molecular diagnostic rate is estimated at 50%, with whole-exome sequencing (WES) among the most successful approaches2-5. For patients in whom WES is uninformative, RNA sequencing (RNA-seq) has shown diagnostic utility in specific tissues and diseases6-8. This includes muscle biopsies from patients with undiagnosed rare muscle disorders6,9, and cultured fibroblasts from patients with mitochondrial disorders7. However, for many individuals, biopsies are not performed for clinical care, and tissues are difficult to access. We sought to assess the utility of RNA-seq from blood as a diagnostic tool for rare diseases of different pathophysiologies. We generated whole-blood RNA-seq from 94 individuals with undiagnosed rare diseases spanning 16 diverse disease categories. We developed a robust approach to compare data from these individuals with large sets of RNA-seq data for controls (n = 1,594 unrelated controls and n = 49 family members) and demonstrated the impacts of expression, splicing, gene and variant filtering strategies on disease gene identification. Across our cohort, we observed that RNA-seq yields a 7.5% diagnostic rate, and an additional 16.7% with improved candidate gene resolution.

Comprehensive Analysis of Rare Variants of 101 Autism-Linked Genes in a Hungarian Cohort of Autism Spectrum Disorder Patients

Background

Autism spectrum disorder (ASD) is genetically and phenotypically heterogeneous. Former genetic studies suggested that both common and rare genetic variants play a role in the etiology. In this study, we aimed to analyze rare variants detected by next generation sequencing (NGS) in an autism cohort from Hungary.

Methods

We investigated the yield of NGS panel sequencing of an unselected ASD cohort (N = 174 ) for the detection of ASD associated syndromes. Besides, we analyzed rare variants in a common disease-rare variant framework and performed rare variant burden analysis and gene enrichment analysis in phenotype based clusters.

Results

We have diagnosed 13 molecularly proven syndromic autism cases. Strongest indicators of syndromic autism were intellectual disability, epilepsy or other neurological plus symptoms. Rare variant analysis on a cohort level confirmed the association of five genes with autism (AUTS2, NHS, NSD1, SLC9A9, and VPS13). We found no correlation between rare variant burden and number of minor malformation or autism severity. We identified four phenotypic clusters, but no specific gene was enriched in a given cluster.

Conclusion

Our study indicates that NGS panel gene sequencing can be useful, where the clinical picture suggests a clinically defined syndromic autism. In this group, targeted panel sequencing may provide reasonable diagnostic yield. Unselected NGS panel screening in the clinic remains controversial, because of uncertain utility, and difficulties of the variant interpretation. However, the detected rare variants may still significantly influence autism risk and subphenotypes in a polygenic model, but to detect the effects of these variants larger cohorts are needed.

Diagnosing rare diseases after the exome

High-throughput sequencing has ushered in a diversity of approaches for identifying genetic variants and understanding genome structure and function. When applied to individuals with rare genetic diseases, these approaches have greatly accelerated gene discovery and patient diagnosis. Over the past decade, exome sequencing has emerged as a comprehensive and cost-effective approach to identify pathogenic variants in the protein-coding regions of the genome. However, for individuals in whom exome-sequencing fails to identify a pathogenic variant, we discuss recent advances that are helping to reduce the diagnostic gap.