A genome-wide association study identified PTPN2 as a population-specific susceptibility gene locus for primary biliary cholangitis.

BACKGROUND AIMS: Previous genome-wide association studies (GWAS) have indicated the involvement of shared (population-non-specific) and non-shared (population-specific) susceptibility genes in the pathogenesis of primary biliary cholangitis (PBC) among European and East-Asian populations. Although a meta-analysis of these distinct populations has recently identified more than 20 novel PBC susceptibility loci, analyses of population-specific genetic architecture are still needed for a more comprehensive search for genetic factors in PBC. APPROACH RESULTS: Protein tyrosine phosphatase non-receptor type 2 (PTPN2) was identified as a novel PBC susceptibility gene locus through a GWAS and subsequent genome-wide meta-analysis involving 2,181 cases and 2,699 controls from the Japanese population (GWAS-lead variant: rs8098858, p=2.6×10-8). In-silico and in-vitro functional analyses indicated that the risk allele of rs2292758, which is a primary functional variant, decreases PTPN2 expression by disrupting Sp1 binding to the PTPN2 promoter in T follicular helper cells (Tfh) and plasmacytoid dendritic cells (pDCs). Infiltration of PTPN2-positive T-cells and pDCs were confirmed in the portal area of the PBC-liver by immunohistochemistry. Furthermore, transcriptomic analysis of PBC-liver samples indicated the presence of a compromised negative feedback loop in-vivo between PTPN2 and IFNG in patients carrying the risk allele of rs2292758. CONCLUSIONS: PTPN2, a novel susceptibility gene for PBC in the Japanese population, may be involved in the pathogenesis of PBC via an insufficient negative feedback loop caused by the PTPN2 risk allele of rs2292758 in IFN signaling. This suggests that PTPN2 could be a potential molecular target for PBC treatment.

rs10924104 in the expression enhancer motif of CD58 confers susceptibility to human autoimmune diseases.

CD58 plays roles in cell adhesion and co-stimulation with antigen presentation from major histocompatibility complex class II on antigen-presenting cells to T-cell antigen receptors on naïve T cells. CD58 reportedly contributes to the development of various human autoimmune diseases. Recently, genome-wide association studies (GWASs) identified CD58 as a susceptibility locus for autoimmune diseases such as systemic lupus erythematosus (SLE), multiple sclerosis (MS), and primary biliary cholangitis (PBC). However, the primary functional variant and molecular mechanisms of susceptibility to autoimmune diseases in the CD58 locus were not clarified. Here, rs10924104, located in the ZNF35-binding motif within the gene expression regulatory motif, was identified as the primary functional variant for SLE, MS, and PBC among genetic variants showing stronger linkage disequilibrium (LD) with GWAS-lead variants in the CD58 locus. Expression-quantitative trait locus (e-QTL) data for each distinct blood cell type and in vitro functional analysis using the CRISPR/Cas9 system corroborated the functional role of rs10924104 in the upregulation of CD58 transcription by the disease-risk allele. Additionally, the strength of disease susceptibility observed in the CD58 locus could be accounted for by the strength of LD between rs10924104 and each GWAS-lead variant. In conclusion, the present study demonstrated for the first time the existence of a shared autoimmune disease-related primary functional variant (i.e., rs10924104) that regulates the expression of CD58. Clarifying the molecular mechanism of disease susceptibility derived from such a shared genetic background is important for understanding human autoimmune diseases and human immunology.

Demographic history of Ryukyu islanders at the southern part of the Japanese Archipelago inferred from whole-genome resequencing data.

The Ryukyu Islands are located in the southernmost part of the Japanese Archipelago and consist of several island groups. Each island group has its own history and culture, which differ from those of mainland Japan. People of the Ryukyu Islands are genetically subdivided; however, their detailed demographic history remains unclear. We report the results of a whole-genome sequencing analysis of a total of 50 Ryukyu islanders, focusing on genetic differentiation between Miyako and Okinawa islanders. We confirmed that Miyako and Okinawa islanders cluster differently in principal component analysis and ADMIXTURE analysis and that there is a population structure among Miyako islanders. The present study supports the hypothesis that population differentiation is primarily caused by genetic drift rather than by differences in the rate of migration from surrounding regions, such as the Japanese main islands or Taiwan. In addition, the genetic cline observed among Miyako and Okinawa islanders can be explained by recurrent migration beyond the bounds of these islands. Our analysis also suggested that the presence of multiple subpopulations during the Neolithic Ryukyu Jomon period is not crucial to explain the modern Ryukyu populations. However, the assumption of multiple subpopulations during the time of admixture with mainland Japanese is necessary to explain the modern Ryukyu populations. Our findings add insights that could help clarify the complex history of populations in the Ryukyu Islands.

rs2013278 in the multiple immunological-trait susceptibility locus CD28 regulates the production of non-functional splicing isoforms.

BACKGROUND: Ligation of CD28 with ligands such as CD80 or CD86 provides a critical second signal alongside antigen presentation by class II major histocompatibility complex expressed on antigen-presenting cells through the T cell antigen receptor for naïve T cell activation. A number of studies suggested that CD28 plays an important role in the pathogenesis of various human diseases. Recent genome-wide association studies (GWASs) identified CD28 as a susceptibility locus for lymphocyte and eosinophil counts, multiple sclerosis, ulcerative colitis, celiac disease, rheumatoid arthritis, asthma, and primary biliary cholangitis. However, the primary functional variant and molecular mechanisms of disease susceptibility in this locus remain to be elucidated. This study aimed to identify the primary functional variant from thousands of genetic variants in the CD28 locus and elucidate its functional effect on the CD28 molecule. RESULTS: Among the genetic variants exhibiting stronger linkage disequilibrium (LD) with all GWAS-lead variants in the CD28 locus, rs2013278, located in the Rbfox binding motif related to splicing regulation, was identified as a primary functional variant related to multiple immunological traits. Relative endogenous expression levels of CD28 splicing isoforms (CD28i and CD28Δex2) compared with full-length CD28 in allele knock-in cell lines generated using CRISPR/Cas9 were directly regulated by rs2013278 (P < 0.05). Although full-length CD28 protein expressed on Jurkat T cells showed higher binding affinity for CD80/CD86, both CD28i and CD28Δex2 encoded loss-of-function isoforms. CONCLUSION: The present study demonstrated for the first time that CD28 has a shared disease-related primary functional variant (i.e., rs2013278) that regulates the CD28 alternative splicing that generates loss-of-function isoforms. They reduce disease risk by inducing anergy of effector T cells that over-react to autoantigens and allergens.

X Chromosome Contribution to the Genetic Architecture of Primary Biliary Cholangitis.

BACKGROUND & AIMS: Genome-wide association studies in primary biliary cholangitis (PBC) have failed to find X chromosome (chrX) variants associated with the disease. Here, we specifically explore the chrX contribution to PBC, a sexually dimorphic complex autoimmune disease. METHODS: We performed a chrX-wide association study, including genotype data from 5 genome-wide association studies (from Italy, United Kingdom, Canada, China, and Japan; 5244 case patients and 11,875 control individuals). RESULTS: Single-marker association analyses found approximately 100 loci displaying P < 5 × 10-4, with the most significant being a signal within the OTUD5 gene (rs3027490; P = 4.80 × 10-6; odds ratio [OR], 1.39; 95% confidence interval [CI], 1.028-1.88; Japanese cohort). Although the transethnic meta-analysis evidenced only a suggestive signal (rs2239452, mapping within the PIM2 gene; OR, 1.17; 95% CI, 1.09-1.26; P = 9.93 × 10-8), the population-specific meta-analysis showed a genome-wide significant locus in East Asian individuals pointing to the same region (rs7059064, mapping within the GRIPAP1 gene; P = 6.2 × 10-9; OR, 1.33; 95% CI, 1.21-1.46). Indeed, rs7059064 tags a unique linkage disequilibrium block including 7 genes: TIMM17B, PQBP1, PIM2, SLC35A2, OTUD5, KCND1, and GRIPAP1, as well as a superenhancer (GH0XJ048933 within OTUD5) targeting all these genes. GH0XJ048933 is also predicted to target FOXP3, the main T-regulatory cell lineage specification factor. Consistently, OTUD5 and FOXP3 RNA levels were up-regulated in PBC case patients (1.75- and 1.64-fold, respectively). CONCLUSIONS: This work represents the first comprehensive study, to our knowledge, of the chrX contribution to the genetics of an autoimmune liver disease and shows a novel PBC-related genome-wide significant locus.

rs9459874 and rs1012656 in CCR6/FGFR1OP confer susceptibility to primary biliary cholangitis.

Primary biliary cholangitis (PBC) is a chronic cholestatic autoimmune liver disease that appears to be strongly influenced by genetic factors. Recently, an international meta-analysis of genome-wide association studies (GWAS) identified CC-Motif Chemokine Receptor-6 (CCR6) and FGFR1 Oncogene-Partner (FGFR1OP) as PBC-susceptibility genes. However, the lead single nucleotide polymorphisms (SNPs) of CCR6/FGFR1OP showed low linkage disequilibrium with each other in East Asian and European populations. Additionally, the primary functional variants and the molecular mechanisms responsible for PBC-susceptibility remain unclear. Here, among the PBC-susceptibility SNPs identified by high-density association mapping in our previous meta-GWAS (Patients: n = 10,516; healthy controls: n = 20,772) within the CCR6/FGFR1OP locus, rs9459874 and rs1012656 were identified as primary functional variants. These functional variants accounted for the effects of GWAS-identified lead SNPs in CCR6/FGFR1OP. Additionally, the roles of rs9459874 and rs1012656 in regulating FGFR1OP transcription and CCR6 translation, respectively, were supported by expression quantitative trait loci (eQTL) analysis and gene editing technology using the CRISPR/Cas9 system. Immunohistochemistry showed higher expression of CCR6 protein in the livers of patients with PBC than in those of a non-diseased control. In conclusion, we identified primary functional variants in CCR6/FGFR1OP and revealed the molecular mechanisms by which these variants confer PBC-susceptibility in an eQTL-dependent or -independent manner. The approach in this study is applicable for the elucidation of the pathogenesis of other autoimmune disorders in which CCR6/FGFR1OP is known as a susceptibility locus, as well as PBC.

An international genome-wide meta-analysis of primary biliary cholangitis: Novel risk loci and candidate drugs.

BACKGROUNDS & AIMS: Primary biliary cholangitis (PBC) is a chronic liver disease in which autoimmune destruction of the small intrahepatic bile ducts eventually leads to cirrhosis. Many patients have inadequate response to licensed medications, motivating the search for novel therapies. Previous genome-wide association studies (GWAS) and meta-analyses (GWMA) of PBC have identified numerous risk loci for this condition, providing insight into its aetiology. We undertook the largest GWMA of PBC to date, aiming to identify additional risk loci and prioritise candidate genes for in silico drug efficacy screening. METHODS: We combined new and existing genotype data for 10,516 cases and 20,772 controls from 5 European and 2 East Asian cohorts. RESULTS: We identified 56 genome-wide significant loci (20 novel) including 46 in European, 13 in Asian, and 41 in combined cohorts; and a 57th genome-wide significant locus (also novel) in conditional analysis of the European cohorts. Candidate genes at newly identified loci include FCRL3, INAVA, PRDM1, IRF7, CCR6, CD226, and IL12RB1, which each play key roles in immunity. Pathway analysis reiterated the likely importance of pattern recognition receptor and TNF signalling, JAK-STAT signalling, and differentiation of T helper (TH)1 and TH17 cells in the pathogenesis of this disease. Drug efficacy screening identified several medications predicted to be therapeutic in PBC, some of which are well-established in the treatment of other autoimmune disorders. CONCLUSIONS: This study has identified additional risk loci for PBC, provided a hierarchy of agents that could be trialled in this condition, and emphasised the value of genetic and genomic approaches to drug discovery in complex disorders. LAY SUMMARY: Primary biliary cholangitis (PBC) is a chronic liver disease that eventually leads to cirrhosis. In this study, we analysed genetic information from 10,516 people with PBC and 20,772 healthy individuals recruited in Canada, China, Italy, Japan, the UK, or the USA. We identified several genetic regions associated with PBC. Each of these regions contains several genes. For each region, we used diverse sources of evidence to help us choose the gene most likely to be involved in causing PBC. We used these 'candidate genes' to help us identify medications that are currently used for treatment of other conditions, which might also be useful for treatment of PBC.

A Guillain-Barré syndrome-associated SIGLEC10 rare variant impairs its recognition of gangliosides.

Guillain-Barré syndrome (GBS), including its variant Miller Fisher syndrome (MFS), is an acute peripheral neuropathy that involves autoimmune mechanisms leading to the production of autoantibodies to gangliosides; sialic acid-containing glycosphingolipids. Although association with various genetic polymorphisms in the major histocompatibility complex (MHC) is shown in other autoimmune diseases, GBS is an exception, showing no such link. No significant association was found by genome wide association studies, suggesting that GBS is not associated with common variants. To address the involvement of rare variants in GBS, we analyzed Siglec-10, a sialic acid-recognizing inhibitory receptor expressed on B cells. Here we demonstrate that two rare variants encoding R47Q and A108V substitutions in the ligand-binding domain are significantly accumulated in patients with GBS. Because of strong linkage disequilibrium, there was no patient carrying only one of them. Recombinant Siglec-10 protein containing R47Q but not A108V shows impaired binding to gangliosides. Homology modeling revealed that the R47Q substitution causes marked alteration in the ligand-binding site. Thus, GBS is associated with a rare variant of the SIGLEC10 gene that impairs ligand binding of Siglec-10. Because Siglec-10 regulates antibody production to sialylated antigens, our finding suggests that Siglec-10 regulates development of GBS by suppressing antibody production to gangliosides, with defects in its function predisposing to disease.

Common risk variants in NPHS1 and TNFSF15 are associated with childhood steroid-sensitive nephrotic syndrome.

To understand the genetics of steroid-sensitive nephrotic syndrome (SSNS), we conducted a genome-wide association study in 987 childhood SSNS patients and 3,206 healthy controls with Japanese ancestry. Beyond known associations in the HLA-DR/DQ region, common variants in NPHS1-KIRREL2 (rs56117924, P=4.94E-20, odds ratio (OR) =1.90) and TNFSF15 (rs6478109, P=2.54E-8, OR=0.72) regions achieved genome-wide significance and were replicated in Korean, South Asian and African populations. Trans-ethnic meta-analyses including Japanese, Korean, South Asian, African, European, Hispanic and Maghrebian populations confirmed the significant associations of variants in NPHS1-KIRREL2 (Pmeta=6.71E-28, OR=1.88) and TNFSF15 (Pmeta=5.40E-11, OR=1.33) loci. Analysis of the NPHS1 risk alleles with glomerular NPHS1 mRNA expression from the same person revealed allele specific expression with significantly lower expression of the transcript derived from the risk haplotype (Wilcox test p=9.3E-4). Because rare pathogenic variants in NPHS1 cause congenital nephrotic syndrome of the Finnish type (CNSF), the present study provides further evidence that variation along the allele frequency spectrum in the same gene can cause or contribute to both a rare monogenic disease (CNSF) and a more complex, polygenic disease (SSNS).

Inhibitory effects of biseokeaniamide A against lipopolysaccharide-induced signal transduction.

Lipopolysaccharides (LPS) are associated with various inflammatory diseases; therefore, the inhibition of LPS-induced nitric oxide (NO) production may have extensive therapeutic applications. We searched for inhibitors of NO production in the LPS-stimulated murine macrophage-like cell line RAW264.7 from MeOH extracts of marine organisms. The MeOH extract of the marine cyanobacterium Okeania sp., collected in Okinawa, Japan, showed inhibitory activity. Biseokeaniamide A was isolated from the MeOH extract by chromatographic separation. Biseokeaniamide A inhibited NO production without cytotoxicity. It reduced inducible nitric oxide synthase levels and suppressed the expression of IL-1ß in LPS-stimulated RAW264.7 cells. Biseokeaniamide A did not inhibit IκBα degradation but inhibited IκBα expression. Thus, biseokeaniamide A, a naturally occurring lipopeptide, was identified as a selective inhibitor of LPS signal transduction.

Genome-wide association studies identify PRKCB as a novel genetic susceptibility locus for primary biliary cholangitis in the Japanese population.

A previous genome-wide association study (GWAS) performed in 963 Japanese individuals (487 primary biliary cholangitis [PBC] cases and 476 healthy controls) identified TNFSF15 (rs4979462) and POU2AF1 (rs4938534) as strong susceptibility loci for PBC. In this study, we performed GWAS in additional 1,923 Japanese individuals (894 PBC cases and 1,029 healthy controls), and combined the results with the previous data. This GWAS, together with a subsequent replication study in an independent set of 7,024 Japanese individuals (512 PBC cases and 6,512 healthy controls), identified PRKCB (rs7404928) as a novel susceptibility locus for PBC (odds ratio [OR] = 1.26, P = 4.13 × 10-9). Furthermore, a primary functional variant of PRKCB (rs35015313) was identified by genotype imputation using a phased panel of 1,070 Japanese individuals from a prospective, general population cohort study and subsequent in vitro functional analyses. These results may lead to improved understanding of the disease pathways involved in PBC, forming a basis for prevention of PBC and development of novel therapeutics.

De novo mutations in epileptic encephalopathies.

Epileptic encephalopathies are a devastating group of severe childhood epilepsy disorders for which the cause is often unknown. Here we report a screen for de novo mutations in patients with two classical epileptic encephalopathies: infantile spasms (n = 149) and Lennox-Gastaut syndrome (n = 115). We sequenced the exomes of 264 probands, and their parents, and confirmed 329 de novo mutations. A likelihood analysis showed a significant excess of de novo mutations in the ∼4,000 genes that are the most intolerant to functional genetic variation in the human population (P = 2.9 × 10(-3)). Among these are GABRB3, with de novo mutations in four patients, and ALG13, with the same de novo mutation in two patients; both genes show clear statistical evidence of association with epileptic encephalopathy. Given the relevant site-specific mutation rates, the probabilities of these outcomes occurring by chance are P = 4.1 × 10(-10) and P = 7.8 × 10(-12), respectively. Other genes with de novo mutations in this cohort include CACNA1A, CHD2, FLNA, GABRA1, GRIN1, GRIN2B, HNRNPU, IQSEC2, MTOR and NEDD4L. Finally, we show that the de novo mutations observed are enriched in specific gene sets including genes regulated by the fragile X protein (P < 10(-8)), as has been reported previously for autism spectrum disorders.

Strong Association of the HLA-DR/DQ Locus with Childhood Steroid-Sensitive Nephrotic Syndrome in the Japanese Population.

Background Nephrotic syndrome is the most common cause of chronic glomerular disease in children. Most of these patients develop steroid-sensitive nephrotic syndrome (SSNS), but the loci conferring susceptibility to childhood SSNS are mainly unknown.Methods We conducted a genome-wide association study (GWAS) in the Japanese population; 224 patients with childhood SSNS and 419 adult healthy controls were genotyped using the Affymetrix Japonica Array in the discovery stage. Imputation for six HLA genes (HLA-A, -C, -B, -DRB1, -DQB1, and -DPB1) was conducted on the basis of Japanese-specific references. We performed genotyping for HLA-DRB1/-DQB1 using a sequence-specific oligonucleotide-probing method on a Luminex platform. Whole-genome imputation was conducted using a phased reference panel of 2049 healthy Japanese individuals. Replication was performed in an independent Japanese sample set including 216 patients and 719 healthy controls. We genotyped candidate single-nucleotide polymorphisms using the DigiTag2 assay.Results The most significant association was detected in the HLA-DR/DQ region and replicated (rs4642516 [minor allele G], combined Pallelic=7.84×10-23; odds ratio [OR], 0.33; 95% confidence interval [95% CI], 0.26 to 0.41; rs3134996 [minor allele A], combined Pallelic=1.72×10-25; OR, 0.29; 95% CI, 0.23 to 0.37). HLA-DRB1*08:02 (Pc=1.82×10-9; OR, 2.62; 95% CI, 1.94 to 3.54) and HLA-DQB1*06:04 (Pc=2.09×10-12; OR, 0.10; 95% CI, 0.05 to 0.21) were considered primary HLA alleles associated with childhood SSNS. HLA-DRB1*08:02-DQB1*03:02 (Pc=7.01×10-11; OR, 3.60; 95% CI, 2.46 to 5.29) was identified as the most significant genetic susceptibility factor.Conclusions The most significant association with childhood SSNS was detected in the HLA-DR/DQ region. Further HLA allele/haplotype analyses should enhance our understanding of molecular mechanisms underlying SSNS.

Significance of functional disease-causal/susceptible variants identified by whole-genome analyses for the understanding of human diseases.

Human genome variation may cause differences in traits and disease risks. Disease-causal/susceptible genes and variants for both common and rare diseases can be detected by comprehensive whole-genome analyses, such as whole-genome sequencing (WGS), using next-generation sequencing (NGS) technology and genome-wide association studies (GWAS). Here, in addition to the application of an NGS as a whole-genome analysis method, we summarize approaches for the identification of functional disease-causal/susceptible variants from abundant genetic variants in the human genome and methods for evaluating their functional effects in human diseases, using an NGS and in silico and in vitro functional analyses. We also discuss the clinical applications of the functional disease causal/susceptible variants to personalized medicine.

IKZF1, a new susceptibility gene for cold medicine-related Stevens-Johnson syndrome/toxic epidermal necrolysis with severe mucosal involvement.

BACKGROUND: Stevens-Johnson syndrome (SJS) and its severe form, toxic epidermal necrolysis (TEN), are acute inflammatory vesiculobullous reactions of the skin and mucous membranes, including the ocular surface, oral cavity, and genitals. These reactions are very rare but are often associated with inciting drugs, infectious agents, or both. OBJECTIVE: We sought to identify susceptibility loci for cold medicine-related SJS/TEN (CM-SJS/TEN) with severe mucosal involvement (SMI). METHODS: A genome-wide association study was performed in 808 Japanese subjects (117 patients with CM-SJS/TEN with SMI and 691 healthy control subjects), and subsequent replication studies were performed in 204 other Japanese subjects (16 cases and 188 control subjects), 117 Korean subjects (27 cases and 90 control subjects), 76 Indian subjects (20 cases and 56 control subjects), and 174 Brazilian subjects (39 cases and 135 control subjects). RESULTS: In addition to the most significant susceptibility region, HLA-A, we identified IKZF1, which encodes Ikaros, as a novel susceptibility gene (meta-analysis, rs4917014 [G vs. T]; odds ratio, 0.5; P = 8.5 × 10(-11)). Furthermore, quantitative ratios of the IKZF1 alternative splicing isoforms Ik1 and Ik2 were significantly associated with rs4917014 genotypes. CONCLUSION: We identified IKZF1 as a susceptibility gene for CM-SJS/TEN with SMI not only in Japanese subjects but also in Korean and Indian subjects and showed that the Ik2/Ik1 ratio might be influenced by IKZF1 single nucleotide polymorphisms, which were significantly associated with susceptibility to CM-SJS/TEN with SMI.

Human primary biliary cirrhosis-susceptible allele of rs4979462 enhances TNFSF15 expression by binding NF-1.

A genome-wide association study (GWAS) identified tumor necrosis factor superfamily member 15 (TNFSF15) as the strongest associated gene with susceptibility to primary biliary cirrhosis (PBC) outside the HLA loci in the Japanese population. However, causal functional variants of the TNFSF15 locus and the molecular mechanism underlying disease susceptibility have not been clarified. Here, to identify the functional causal variants of the TNFSF15 locus, integrated analysis comprising in silico analysis, a case-control association study and in vitro functional analysis was performed. Initially, 32 functional candidate single-nucleotide polymorphisms (SNPs) in the expression regulatory motifs, the coding region, or the untranslated regions (UTRs) of the TNFSF15 locus were selected by in silico analysis. By the case-control association studies using PBC patients (n = 1279) and healthy controls (n = 1091) in the Japanese population, rs4979462 [P = 1.85 × 10(-14) (our previous study)], rs56211063 (P = 2.21 × 10(-14)), and rs55768522 (r(2) = 1 with rs4979462) were likely candidates for causal variants. Among these SNPs, rs4979462 was identified as the causal variant by in vitro functional analysis using luciferase assay and electrophoretic mobility shift assay (EMSA). Super-shift assay clarified that PBC-susceptible allele of rs4979462 generated a novel NF-1 binding site. Moreover, higher endogenous TNFSF15 protein and mRNA expression levels were observed in individuals with the PBC-susceptible allele of rs4979462. This study identified the causal variant for PBC susceptibility in the TNFSF15 locus and clarified its underlying molecular mechanism. TNFSF15 and NF-1 are considered to be potential targets for the treatment of PBC.

Mutation in TECPR2 reveals a role for autophagy in hereditary spastic paraparesis.

We studied five individuals from three Jewish Bukharian families affected by an apparently autosomal-recessive form of hereditary spastic paraparesis accompanied by severe intellectual disability, fluctuating central hypoventilation, gastresophageal reflux disease, wake apnea, areflexia, and unique dysmorphic features. Exome sequencing identified one homozygous variant shared among all affected individuals and absent in controls: a 1 bp frameshift TECPR2 deletion leading to a premature stop codon and predicting significant degradation of the protein. TECPR2 has been reported as a positive regulator of autophagy. We thus examined the autophagy-related fate of two key autophagic proteins, SQSTM1 (p62) and MAP1LC3B (LC3), in skin fibroblasts of an affected individual, as compared to a healthy control, and found that both protein levels were decreased and that there was a more pronounced decrease in the lipidated form of LC3 (LC3II). siRNA knockdown of TECPR2 showed similar changes, consistent with aberrant autophagy. Our results are strengthened by the fact that autophagy dysfunction has been implicated in a number of other neurodegenerative diseases. The discovered TECPR2 mutation implicates autophagy, a central intracellular mechanism, in spastic paraparesis.

Whole-exome sequencing in undiagnosed genetic diseases: interpreting 119 trios.

PURPOSE: Despite the recognized clinical value of exome-based diagnostics, methods for comprehensive genomic interpretation remain immature. Diagnoses are based on known or presumed pathogenic variants in genes already associated with a similar phenotype. Here, we extend this paradigm by evaluating novel bioinformatics approaches to aid identification of new gene-disease associations. METHODS: We analyzed 119 trios to identify both diagnostic genotypes in known genes and candidate genotypes in novel genes. We considered qualifying genotypes based on their population frequency and in silico predicted effects we also characterized the patterns of genotypes enriched among this collection of patients. RESULTS: We obtained a genetic diagnosis for 29 (24%) of our patients. We showed that patients carried an excess of damaging de novo mutations in intolerant genes, particularly those shown to be essential in mice (P = 3.4 × 10(-8)). This enrichment is only partially explained by mutations found in known disease-causing genes. CONCLUSION: This work indicates that the application of appropriate bioinformatics analyses to clinical sequence data can also help implicate novel disease genes and suggest expanded phenotypes for known disease genes. These analyses further suggest that some cases resolved by whole-exome sequencing will have direct therapeutic implications.