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

RACK1 Associates with Muscarinic Receptors and Regulates M2 Receptor Trafficking

Receptor internalization from the cell surface occurs through several mechanisms. Some of these mechanisms, such as clathrin coated pits, are well understood. The M2 muscarinic acetylcholine receptor undergoes internalization via a poorly-defined clathrin-independent mechanism. We used isotope coded affinity tagging and mass spectrometry to identify the scaffolding protein, receptor for activated C kinase (RACK1) as a protein enriched in M2-immunoprecipitates from M2-expressing cells over those of non-M2 expressing cells. Treatment of cells with the agonist carbachol disrupted the interaction of RACK1 with M2. We further found that RACK1 overexpression inhibits the internalization and subsequent down regulation of the M2 receptor in a receptor subtype-specific manner. Decreased RACK1 expression increases the rate of agonist internalization of the M2 receptor, but decreases the extent of subsequent down-regulation. These results suggest that RACK1 may both interfere with agonist-induced sequestration and be required for subsequent targeting of internalized M2 receptors to the degradative pathway

A toolkit for genetics providers in follow‐up of patients with non‐diagnostic exome sequencing

There are approximately 7,000 rare diseases affecting 25-30 million Americans, with 80% estimated to have a genetic basis. This presents a challenge for genetics practitioners to determine appropriate testing, make accurate diagnoses, and conduct up-to-date patient management. Exome sequencing (ES) is a comprehensive diagnostic approach, but only 25%-41% of the patients receive a molecular diagnosis. The remaining three-fifths to three-quarters of patients undergoing ES remain undiagnosed. The Stanford Center for Undiagnosed Diseases (CUD), a clinical site of the Undiagnosed Diseases Network, evaluates patients with undiagnosed and rare diseases using a combination of methods including ES. Frequently these patients have non-diagnostic ES results, but strategic follow-up techniques identify diagnoses in a subset. We present techniques used at the CUD that can be adopted by genetics providers in clinical follow-up of cases where ES is non-diagnostic. Solved case examples illustrate different types of non-diagnostic results and the additional techniques that led to a diagnosis. Frequent approaches include segregation analysis, data reanalysis, genome sequencing, additional variant identification, careful phenotype-disease correlation, confirmatory testing, and case matching. We also discuss prioritization of cases for additional analyses

Extracutaneous manifestations in phacomatosis cesioflammea and cesiomarmorata: Case series and literature review

Phacomatosis pigmentovascularis (PPV) comprises a family of rare conditions that feature vascular abnormalities and melanocytic lesions that can be solely cutaneous or multisystem in nature. Recently published work has demonstrated that both vascular and melanocytic abnormalities in PPV of the cesioflammea and cesiomarmorata subtypes can result from identical somatic mosaic activating mutations in the genes GNAQ and GNA11. Here, we present three new cases of PPV with features of the cesioflammea and/or cesiomarmorata subtypes and mosaic mutations in GNAQ or GNA11. To better understand the risk of potentially occult complications faced by such patients we additionally reviewed 176 cases published in the literature. We report the frequency of clinical findings, their patterns of co-occurrence as well as published recommendations for surveillance after diagnosis. Features assessed include: capillary malformation; dermal and ocular melanocytosis; glaucoma; limb asymmetry; venous malformations; and central nervous system (CNS) anomalies, such as ventriculomegaly and calcifications. We found that ocular findings are common in patients with phacomatosis cesioflammea and cesiomarmorata. Facial vascular involvement correlates with a higher risk of seizures (p = .0066). Our genetic results confirm the role of mosaic somatic mutations in GNAQ and GNA11 in phacomatosis cesioflammea and cesiomarmorata. Their clinical and molecular findings place these conditions on a clinical spectrum encompassing other GNAQ and GNA11 related disorders and inform recommendations for their management

Identification of rare and common regulatory variants in pluripotent cells using population-scale transcriptomics.

Induced pluripotent stem cells (iPSCs) are an established cellular system to study the impact of genetic variants in derived cell types and developmental contexts. However, in their pluripotent state, the disease impact of genetic variants is less well known. Here, we integrate data from 1,367 human iPSC lines to comprehensively map common and rare regulatory variants in human pluripotent cells. Using this population-scale resource, we report hundreds of new colocalization events for human traits specific to iPSCs, and find increased power to identify rare regulatory variants compared with somatic tissues. Finally, we demonstrate how iPSCs enable the identification of causal genes for rare diseases

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

Spectrum of neurodevelopmental disease associated with the GNAO1 guanosine triphosphate-binding region

Objective To characterize the phenotypic spectrum associated with GNAO1 variants and establish genotype-protein structure-phenotype relationships. Methods We evaluated the phenotypes of 14 patients with GNAO1 variants, analyzed their variants for potential pathogenicity, and mapped them, along with those in the literature, on a three-dimensional structural protein model. Results The 14 patients in our cohort, including one sibling pair, had 13 distinct, heterozygous GNAO1 variants classified as pathogenic or likely pathogenic. We attributed the same variant in two siblings to parental mosaicism. Patients initially presented with seizures beginning in the first 3 months of life (8/14), developmental delay (4/14), hypotonia (1/14), or movement disorder (1/14). All patients had hypotonia and developmental delay ranging from mild to severe. Nine had epilepsy, and nine had movement disorders, including dystonia, ataxia, chorea, and dyskinesia. The 13 GNAO1 variants in our patients are predicted to result in amino acid substitutions or deletions in the GNAO1 guanosine triphosphate (GTP)-binding region, analogous to those in previous publications. Patients with variants affecting amino acids 207-221 had only movement disorder and hypotonia. Patients with variants affecting the C-terminal region had the mildest phenotypes.