Glycine decarboxylase activity drives non-small cell lung cancer tumor-initiating cells and tumorigenesis.

Identification of the factors critical to the tumor-initiating cell (TIC) state may open new avenues in cancer therapy. Here we show that the metabolic enzyme glycine decarboxylase (GLDC) is critical for TICs in non-small cell lung cancer (NSCLC). TICs from primary NSCLC tumors express high levels of the oncogenic stem cell factor LIN28B and GLDC, which are both required for TIC growth and tumorigenesis. Overexpression of GLDC and other glycine/serine enzymes, but not catalytically inactive GLDC, promotes cellular transformation and tumorigenesis. We found that GLDC induces dramatic changes in glycolysis and glycine/serine metabolism, leading to changes in pyrimidine metabolism to regulate cancer cell proliferation. In the clinic, aberrant activation of GLDC correlates with poorer survival in lung cancer patients, and aberrant GLDC expression is observed in multiple cancer types. This link between glycine metabolism and tumorigenesis may provide novel targets for advancing anticancer therapy.

Alterations to DNA methylation patterns induced by chemotherapy treatment are associated with negative impacts on the olfactory pathway.

BACKGROUND: Exposure to cytotoxic chemotherapy treatment may alter DNA methylation (DNAm) in breast cancer patients. METHODS: We performed DNAm analysis in 125 breast cancer patients with blood drawn before and after chemotherapy, using the Illumina MethylationEPIC array. DNAm changes of 588,798 individual CpGs (including 41,207 promoter regions) were evaluated using linear regression models adjusted for monocyte proportion. Gene set enrichment analyses (GSEA) were conducted to identify key Gene Ontology (GO) biological processes or Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways associated with chemotherapy. Results were validated in a separate cohort of breast cancer patients who were treated (n = 1273) and not treated (n = 872) by chemotherapy (1808 blood, 337 saliva). RESULTS: A total of 141 differentially methylated CpGs and 11 promoters were significantly associated with chemotherapy after multiple testing corrections in both the paired sample and single time point analyses. GSEA of promoter regions (pre-ranked by test statistics) identified six suppressed biological processes (p < 4.67e-8) related to sensory perception and detection of chemical stimuli, including smell perception (GO:0007606, GO:0007608, GO:0009593, GO:0050906, GO:0050907, and GO:0050911). The same six biological processes were significantly suppressed in the validation dataset (p < 9.02e-14). The KEGG pathway olfactory transduction (hsa04740) was also found to be significantly suppressed (ppaired-samples = 1.72e-9, psingle-timepoint-blood = 2.03e-15 and psingle-timepoint-saliva = 7.52e-56). CONCLUSION: The enrichment of imprinted genes within biological processes and pathways suggests a biological mechanism by which chemotherapy could affect the perception of smell.

oFlowSeq: a quantitative approach to identify protein coding mutations affecting cell type enrichment using mosaic CRISPR-Cas9 edited cerebral organoids.

Cerebral organoids are comprised of diverse cell types found in the developing human brain, and can be leveraged in the identification of critical cell types perturbed by genetic risk variants in common, neuropsychiatric disorders. There is great interest in developing high-throughput technologies to associate genetic variants with cell types. Here, we describe a high-throughput, quantitative approach (oFlowSeq) by utilizing CRISPR-Cas9, FACS sorting, and next-generation sequencing. Using oFlowSeq, we found that deleterious mutations in autism-associated gene KCTD13 resulted in increased proportions of Nestin+ cells and decreased proportions of TRA-1-60+ cells within mosaic cerebral organoids. We further identified that a locus-wide CRISPR-Cas9 survey of another 18 genes in the 16p11.2 locus resulted in most genes with > 2% maximum editing efficiencies for short and long indels, suggesting a high feasibility for an unbiased, locus-wide experiment using oFlowSeq. Our approach presents a novel method to identify genotype-to-cell type imbalances in an unbiased, high-throughput, quantitative manner.

Breast cancer risk stratification using genetic and non-genetic risk assessment tools for 246,142 women in the UK Biobank.

PURPOSE: The benefit of using individual risk prediction tools to identify high-risk individuals for breast cancer (BC) screening is uncertain, despite the personalized approach of risk-based screening. METHODS: We studied the overlap of predicted high-risk individuals among 246,142 women enrolled in the UK Biobank. Risk predictors assessed include the Gail model (Gail), BC family history (FH, binary), BC polygenic risk score (PRS), and presence of loss-of-function (LoF) variants in BC predisposition genes. Youden J-index was used to select optimal thresholds for defining high-risk. RESULTS: In total, 147,399 were considered at high risk for developing BC within the next 2 years by at least 1 of the 4 risk prediction tools examined (Gail2-year > 0.5%: 47%, PRS2-yea r > 0.7%: 30%, FH: 6%, and LoF: 1%); 92,851 (38%) were flagged by only 1 risk predictor. The overlap between individuals flagged as high-risk because of genetic (PRS) and Gail model risk factors was 30%. The best-performing combinatorial model comprises a union of high-risk women identified by PRS, FH, and, LoF (AUC2-year [95% CI]: 62.2 [60.8 to 63.6]). Assigning individual weights to each risk prediction tool increased discriminatory ability. CONCLUSION: Risk-based BC screening may require a multipronged approach that includes PRS, predisposition genes, FH, and other recognized risk factors.

Multi-center evaluation of artificial intelligent imaging and clinical models for predicting neoadjuvant chemotherapy response in breast cancer.

BACKGROUND: Neoadjuvant chemotherapy (NAC) plays an important role in the management of locally advanced breast cancer. It allows for downstaging of tumors, potentially allowing for breast conservation. NAC also allows for in-vivo testing of the tumors' response to chemotherapy and provides important prognostic information. There are currently no clearly defined clinical models that incorporate imaging with clinical data to predict response to NAC. Thus, the aim of this work is to develop a predictive AI model based on routine CT imaging and clinical parameters to predict response to NAC. METHODS: The CT scans of 324 patients with NAC from multiple centers in Singapore were used in this study. Four different radiomics models were built for predicting pathological complete response (pCR): first two were based on textural features extracted from peri-tumoral and tumoral regions, the third model based on novel space-resolved radiomics which extract feature maps using voxel-based radiomics and the fourth model based on deep learning (DL). Clinical parameters were included to build a final prognostic model. RESULTS: The best performing models were based on space-resolved and DL approaches. Space-resolved radiomics improves the clinical AUCs of pCR prediction from 0.743 (0.650 to 0.831) to 0.775 (0.685 to 0.860) and our DL model improved it from 0.743 (0.650 to 0.831) to 0.772 (0.685 to 0.853). The tumoral radiomics model performs the worst with no improvement of the AUC from the clinical model. The peri-tumoral combined model gives moderate performance with an AUC of 0.765 (0.671 to 0.855). CONCLUSIONS: Radiomics features extracted from diagnostic CT augment the predictive ability of pCR when combined with clinical features. The novel space-resolved radiomics and DL radiomics approaches outperformed conventional radiomics techniques.

Overlap of high-risk individuals predicted by family history, and genetic and non-genetic breast cancer risk prediction models: implications for risk stratification.

BACKGROUND: Family history, and genetic and non-genetic risk factors can stratify women according to their individual risk of developing breast cancer. The extent of overlap between these risk predictors is not clear. METHODS: In this case-only analysis involving 7600 Asian breast cancer patients diagnosed between age 30 and 75 years, we examined identification of high-risk patients based on positive family history, the Gail model 5-year absolute risk [5yAR] above 1.3%, breast cancer predisposition genes (protein-truncating variants [PTV] in ATM, BRCA1, BRCA2, CHEK2, PALB2, BARD1, RAD51C, RAD51D, or TP53), and polygenic risk score (PRS) 5yAR above 1.3%. RESULTS: Correlation between 5yAR (at age of diagnosis) predicted by PRS and the Gail model was low (r=0.27). Fifty-three percent of breast cancer patients (n=4041) were considered high risk by one or more classification criteria. Positive family history, PTV carriership, PRS, or the Gail model identified 1247 (16%), 385 (5%), 2774 (36%), and 1592 (21%) patients who were considered at high risk, respectively. In a subset of 3227 women aged below 50 years, the four models studied identified 470 (15%), 213 (7%), 769 (24%), and 325 (10%) unique patients who were considered at high risk, respectively. For younger women, PRS and PTVs together identified 745 (59% of 1276) high-risk individuals who were not identified by the Gail model or family history. CONCLUSIONS: Family history and genetic and non-genetic risk stratification tools have the potential to complement one another to identify women at high risk.

The Genetic Landscape of Diamond-Blackfan Anemia.

Diamond-Blackfan anemia (DBA) is a rare bone marrow failure disorder that affects 7 out of 1,000,000 live births and has been associated with mutations in components of the ribosome. In order to characterize the genetic landscape of this heterogeneous disorder, we recruited a cohort of 472 individuals with a clinical diagnosis of DBA and performed whole-exome sequencing (WES). We identified relevant rare and predicted damaging mutations for 78% of individuals. The majority of mutations were singletons, absent from population databases, predicted to cause loss of function, and located in 1 of 19 previously reported ribosomal protein (RP)-encoding genes. Using exon coverage estimates, we identified and validated 31 deletions in RP genes. We also observed an enrichment for extended splice site mutations and validated their diverse effects using RNA sequencing in cell lines obtained from individuals with DBA. Leveraging the size of our cohort, we observed robust genotype-phenotype associations with congenital abnormalities and treatment outcomes. We further identified rare mutations in seven previously unreported RP genes that may cause DBA, as well as several distinct disorders that appear to phenocopy DBA, including nine individuals with biallelic CECR1 mutations that result in deficiency of ADA2. However, no new genes were identified at exome-wide significance, suggesting that there are no unidentified genes containing mutations readily identified by WES that explain >5% of DBA-affected case subjects. Overall, this report should inform not only clinical practice for DBA-affected individuals, but also the design and analysis of rare variant studies for heterogeneous Mendelian disorders.

An enhanced CRISPR repressor for targeted mammalian gene regulation.

The RNA-guided endonuclease Cas9 can be converted into a programmable transcriptional repressor, but inefficiencies in target-gene silencing have limited its utility. Here we describe an improved Cas9 repressor based on the C-terminal fusion of a rationally designed bipartite repressor domain, KRAB-MeCP2, to nuclease-dead Cas9. We demonstrate the system's superiority in silencing coding and noncoding genes, simultaneously repressing a series of target genes, improving the results of single and dual guide RNA library screens, and enabling new architectures of synthetic genetic circuits.

Electrocardiographic predictors of successful resynchronization of left bundle branch block by His bundle pacing.

BACKGROUND: His bundle pacing (HBP) is an alternative to biventricular pacing (BVP) for delivering cardiac resynchronization therapy (CRT) in patients with heart failure and left bundle branch block (LBBB). It is not known whether ventricular activation times and patterns achieved by HBP are equivalent to intact conduction systems and not all patients with LBBB are resynchronized by HBP. OBJECTIVE: To compare activation times and patterns of His-CRT with BVP-CRT, LBBB and intact conduction systems. METHODS: In patients with LBBB, noninvasive epicardial mapping (ECG imaging) was performed during BVP and temporary HBP. Intrinsic activation was mapped in all subjects. Left ventricular activation times (LVAT) were measured and epicardial propagation mapping (EPM) was performed, to visualize epicardial wavefronts. Normal activation pattern and a normal LVAT range were determined from normal subjects. RESULTS: Forty-five patients were included, 24 with LBBB and LV impairment, and 21 with normal 12-lead ECG and LV function. In 87.5% of patients with LBBB, His-CRT successfully shortened LVAT by ≥10 ms. In 33.3%, His-CRT resulted in complete ventricular resynchronization, with activation times and patterns indistinguishable from normal subjects. EPM identified propagation discontinuity artifacts in 83% of patients with LBBB. This was the best predictor of whether successful resynchronization was achieved by HBP (logarithmic odds ratio, 2.19; 95% confidence interval, 0.07-4.31; p = .04). CONCLUSION: Noninvasive electrocardiographic mapping appears to identify patients whose LBBB can be resynchronized by HBP. In contrast to BVP, His-CRT may deliver the maximum potential ventricular resynchronization, returning activation times, and patterns to those seen in normal hearts.

Evaluation of phospho-histone H3 in Asian triple-negative breast cancer using multiplex immunofluorescence.

PURPOSE: We used multiplex immunofluorescence (mIF) to determine whether mitotic rate represents an independent prognostic marker in triple-negative breast cancer (TNBC). Secondary aims were to confirm the prognostic significance of immune cells in TNBC, and to investigate the relationship between immune cells and proliferating tumour cells. METHODS: A retrospective Asian cohort of 298 patients with TNBC diagnosed from 2003 to 2015 at the Singapore General Hospital was used in the present study. Formalin-fixed, paraffin-embedded breast cancer samples were analysed on tissue microarrays using mIF, which combined phospho-histone H3 (pHH3) expression with cytokeratin (CK) and leukocyte common antigen (CD45) expression to identify tumour and immune cells, respectively. RESULTS: Multivariate analysis showed that a high pHH3 index was associated with significantly improved overall survival (OS; p = 0.004), but this was not significantly associated with disease-free survival (DFS; p = 0.22). Similarly, multivariate analysis also revealed that a pHH3 positive count of > 1 cell per high-power field in the malignant epithelial compartment was an independent favourable prognostic marker for OS (p = 0.033) but not for DFS (p = 0.250). Furthermore, a high CD45 index was an independent favourable prognostic marker for DFS (p = 0.018), and there was a significant positive correlation between CD45 and pHH3 index (Spearman rank correlation coefficient, 0.250; p < 0.001). CONCLUSIONS: Mitotic rates as determined by pHH3 expression in epithelial cells are significantly associated with improved survival in TNBC. mIF analysis of pHH3 in combination with CK and CD45 could help clinicians in prognosticating patients with TNBC.

Genome-wide Analysis of Body Proportion Classifies Height-Associated Variants by Mechanism of Action and Implicates Genes Important for Skeletal Development.

Human height is a composite measurement, reflecting the sum of leg, spine, and head lengths. Many common variants influence total height, but the effects of these or other variants on the components of height (body proportion) remain largely unknown. We studied sitting height ratio (SHR), the ratio of sitting height to total height, to identify such effects in 3,545 African Americans and 21,590 individuals of European ancestry. We found that SHR is heritable: 26% and 39% of the total variance of SHR can be explained by common variants in European and African Americans, respectively, and global European admixture is negatively correlated with SHR in African Americans (r(2) ≈ 0.03). Six regions reached genome-wide significance (p < 5 × 10(-8)) for association with SHR and overlapped biological candidate genes, including TBX2 and IGFBP3. We found that 130 of 670 height-associated variants are nominally associated (p < 0.05) with SHR, more than expected by chance (p = 5 × 10(-40)). At these 130 loci, the height-increasing alleles are associated with either a decrease (71 loci) or increase (59 loci) in SHR, suggesting that different height loci disproportionally affect either leg length or spine/head length. Pathway analyses via DEPICT revealed that height loci affecting SHR, and especially those affecting leg length, show enrichment of different biological pathways (e.g., bone/cartilage/growth plate pathways) than do loci with no effect on SHR (e.g., embryonic development). These results highlight the value of using a pair of related but orthogonal phenotypes, in this case SHR with height, as a prism to dissect the biology underlying genetic associations in polygenic traits and diseases.

Association of mitochondrial DNA content in peripheral blood with cancer-related fatigue and chemotherapy-related cognitive impairment in early-stage breast cancer patients: a prospective cohort study.

PURPOSE: Cancer-related fatigue (CRF) and chemotherapy-related cognitive impairment (CRCI) are reported to be associated with mitochondrial dysfunction. Hence, mitochondrial DNA (mtDNA) content, a biomarker of mitochondrial dysfunction, is hypothesized to correlate with the onset of CRF and CRCI. This study aims to evaluate the association between peripheral blood mtDNA content reduction and severity of CRF and CRCI in patients receiving chemotherapy. METHODS: This was a prospective cohort study. Early-stage breast cancer patients receiving anthracycline- or taxane-based chemotherapy were recruited. CRF was assessed using MFSI-SF, and CRCI was assessed using FACT-Cog and CANTAB at two timepoints: baseline (T1; prior to treatment) and 6 weeks after initiation of treatment (T2). mtDNA content was measured at both timepoints using real-time quantitative polymerase chain reaction. Multiple logistic regression was utilized to evaluate the association between mtDNA reduction and worsening of CRF and CRCI, adjusting for age, anxiety, insomnia, plasma cytokines concentrations, and other clinically important covariates. RESULTS: A total of 108 patients (age 52.0 ± 9.2 years; 82.4% Chinese; 64.8% receiving anthracycline-based chemotherapy) were recruited. Proportions of patients with worsening of CRF increased from the lower to the upper quartiles of mtDNA reduction (22.2, 33.3, 55.6, and 63.0% in quartiles 1, 2, 3, and 4, respectively, p = 0.001 for trend). Reduction of mtDNA content was significantly greater among those with worsening of CRF and CRCI compared to those without CRF [mean reduction (± SD): 36.5 (46.1) vs. 9.4 (34.5), p < 0.001]. After adjusting for covariates, every 1-unit reduction of the mtDNA content was associated with a 4% increased risk for worsening of CRF (95% CI, 1-6%; p = 0.009). CONCLUSIONS: This is the first study to show that the reduction of mtDNA content in peripheral blood is associated with the onset of CRF in patients receiving chemotherapy. Further validation studies are required to confirm the findings.

A novel approach to mapping the atrial ganglionated plexus network by generating a distribution probability atlas.

INTRODUCTION: The ganglionated plexuses (GPs) of the intrinsic cardiac autonomic system are implicated in arrhythmogenesis. GP localization by stimulation of the epicardial fat pads to produce atrioventricular dissociating (AVD) effects is well described. We determined the anatomical distribution of the left atrial GPs that influence atrioventricular (AV) dissociation. METHODS AND RESULTS: High frequency stimulation was delivered through a Smart-Touch catheter in the left atrium of patients undergoing atrial fibrillation (AF) ablation. Three dimensional locations of points tested throughout the entire chamber were recorded on the CARTO™ system. Impact on the AV conduction was categorized as ventricular asystole, bradycardia, or no effect. CARTO maps were exported, registered, and transformed onto a reference left atrial geometry using a custom software, enabling data from multiple patients to be overlaid. In 28 patients, 2108 locations were tested and 283 sites (13%) demonstrated (AVD-GP) effects. There were 10 AVD-GPs (interquartile range, 11.5) per patient. Eighty percent (226) produced asystole and 20% (57) showed bradycardia. The distribution of the two groups was very similar. Highest probability of AVD-GPs (>20%) was identified in: inferoseptal portion (41%) and right inferior pulmonary vein base (30%) of the posterior wall, right superior pulmonary vein antrum (31%). CONCLUSION: It is feasible to map the entire left atrium for AVD-GPs before AF ablation. Aggregated data from multiple patients, producing a distribution probability atlas of AVD-GPs, identified three regions with a higher likelihood for finding AVD-GPs and these matched the histological descriptions. This approach could be used to better characterize the autonomic network.

Patterns and rates of exonic de novo mutations in autism spectrum disorders.

Autism spectrum disorders (ASD) are believed to have genetic and environmental origins, yet in only a modest fraction of individuals can specific causes be identified. To identify further genetic risk factors, here we assess the role of de novo mutations in ASD by sequencing the exomes of ASD cases and their parents (n = 175 trios). Fewer than half of the cases (46.3%) carry a missense or nonsense de novo variant, and the overall rate of mutation is only modestly higher than the expected rate. In contrast, the proteins encoded by genes that harboured de novo missense or nonsense mutations showed a higher degree of connectivity among themselves and to previous ASD genes as indexed by protein-protein interaction screens. The small increase in the rate of de novo events, when taken together with the protein interaction results, are consistent with an important but limited role for de novo point mutations in ASD, similar to that documented for de novo copy number variants. Genetic models incorporating these data indicate that most of the observed de novo events are unconnected to ASD; those that do confer risk are distributed across many genes and are incompletely penetrant (that is, not necessarily sufficient for disease). Our results support polygenic models in which spontaneous coding mutations in any of a large number of genes increases risk by 5- to 20-fold. Despite the challenge posed by such models, results from de novo events and a large parallel case-control study provide strong evidence in favour of CHD8 and KATNAL2 as genuine autism risk factors.

An excess of risk-increasing low-frequency variants can be a signal of polygenic inheritance in complex diseases.

In most complex diseases, much of the heritability remains unaccounted for by common variants. It has been postulated that lower-frequency variants contribute to the remaining heritability. Here, we describe a method to test for polygenic inheritance from lower-frequency variants by using GWAS summary association statistics. We explored scenarios with many causal low-frequency variants and showed that there is more power to detect risk variants than to detect protective variants, resulting in an increase in the ratio of detected risk to protective variants (R/P ratio). Such an excess can also occur if risk variants are present and kept at lower frequencies because of negative selection. The R/P ratio can be falsely elevated because of reasons unrelated to polygenic inheritance, such as uneven sample sizes or asymmetric population stratification, so precautions to correct for these confounders are essential. We tested our method on published GWAS results and observed a strong signal in some diseases (schizophrenia and type 2 diabetes) but not others. We also explored the shared genetic component in overlapping phenotypes related to inflammatory bowel disease (Crohn disease [CD] and ulcerative colitis [UC]) and diabetic nephropathy (macroalbuminuria and end-stage renal disease [ESRD]). Although the signal was still present when both CD and UC were jointly analyzed, the signal was lost when macroalbuminuria and ESRD were jointly analyzed, suggesting that these phenotypes should best be studied separately. Thus, our method may also help guide the design of future genetic studies of various traits and diseases.

A novel test for recessive contributions to complex diseases implicates Bardet-Biedl syndrome gene BBS10 in idiopathic type 2 diabetes and obesity.

Rare-variant association studies in common, complex diseases are customarily conducted under an additive risk model in both single-variant and burden testing. Here, we describe a method to improve detection of rare recessive variants in complex diseases termed RAFT (recessive-allele-frequency-based test). We found that RAFT outperforms existing approaches when the variant influences disease risk in a recessive manner on simulated data. We then applied our method to 1,791 Finnish individuals with type 2 diabetes (T2D) and 2,657 matched control subjects. In BBS10, we discovered a rare variant (c.1189A>G [p.Ile397Val]; rs202042386) that confers risk of T2D in a recessive state (p = 1.38 × 10(-6)) and would be missed by conventional methods. Testing of this variant in an established in vivo zebrafish model confirmed the variant to be pathogenic. Taken together, these data suggest that RAFT can effectively reveal rare recessive contributions to complex diseases overlooked by conventional association tests.

Distribution and medical impact of loss-of-function variants in the Finnish founder population.

Exome sequencing studies in complex diseases are challenged by the allelic heterogeneity, large number and modest effect sizes of associated variants on disease risk and the presence of large numbers of neutral variants, even in phenotypically relevant genes. Isolated populations with recent bottlenecks offer advantages for studying rare variants in complex diseases as they have deleterious variants that are present at higher frequencies as well as a substantial reduction in rare neutral variation. To explore the potential of the Finnish founder population for studying low-frequency (0.5-5%) variants in complex diseases, we compared exome sequence data on 3,000 Finns to the same number of non-Finnish Europeans and discovered that, despite having fewer variable sites overall, the average Finn has more low-frequency loss-of-function variants and complete gene knockouts. We then used several well-characterized Finnish population cohorts to study the phenotypic effects of 83 enriched loss-of-function variants across 60 phenotypes in 36,262 Finns. Using a deep set of quantitative traits collected on these cohorts, we show 5 associations (p<5×10⁻8) including splice variants in LPA that lowered plasma lipoprotein(a) levels (P = 1.5×10⁻¹¹7). Through accessing the national medical records of these participants, we evaluate the LPA finding via Mendelian randomization and confirm that these splice variants confer protection from cardiovascular disease (OR = 0.84, P = 3×10⁻4), demonstrating for the first time the correlation between very low levels of LPA in humans with potential therapeutic implications for cardiovascular diseases. More generally, this study articulates substantial advantages for studying the role of rare variation in complex phenotypes in founder populations like the Finns and by combining a unique population genetic history with data from large population cohorts and centralized research access to National Health Registers.

Ataxia, dementia, and hypogonadotropism caused by disordered ubiquitination.

BACKGROUND: The combination of ataxia and hypogonadism was first described more than a century ago, but its genetic basis has remained elusive. METHODS: We performed whole-exome sequencing in a patient with ataxia and hypogonadotropic hypogonadism, followed by targeted sequencing of candidate genes in similarly affected patients. Neurologic and reproductive endocrine phenotypes were characterized in detail. The effects of sequence variants and the presence of an epistatic interaction were tested in a zebrafish model. RESULTS: Digenic homozygous mutations in RNF216 and OTUD4, which encode a ubiquitin E3 ligase and a deubiquitinase, respectively, were found in three affected siblings in a consanguineous family. Additional screening identified compound heterozygous truncating mutations in RNF216 in an unrelated patient and single heterozygous deleterious mutations in four other patients. Knockdown of rnf216 or otud4 in zebrafish embryos induced defects in the eye, optic tectum, and cerebellum; combinatorial suppression of both genes exacerbated these phenotypes, which were rescued by nonmutant, but not mutant, human RNF216 or OTUD4 messenger RNA. All patients had progressive ataxia and dementia. Neuronal loss was observed in cerebellar pathways and the hippocampus; surviving hippocampal neurons contained ubiquitin-immunoreactive intranuclear inclusions. Defects were detected at the hypothalamic and pituitary levels of the reproductive endocrine axis. CONCLUSIONS: The syndrome of hypogonadotropic hypogonadism, ataxia, and dementia can be caused by inactivating mutations in RNF216 or by the combination of mutations in RNF216 and OTUD4. These findings link disordered ubiquitination to neurodegeneration and reproductive dysfunction and highlight the power of whole-exome sequencing in combination with functional studies to unveil genetic interactions that cause disease. (Funded by the National Institutes of Health and others.).