shinyExprPortal: a configurable 'shiny' portal for sharing analysis of molecular expression data.

MOTIVATION: The scale of omics research presents many obstacles to full sharing and access to analysis results. Current publication models impose limits on the number of pages and figures, requiring careful preparation and selection of content. At the same time, depositing data in open repositories significantly shifts the burden of access and reproduction to readers, who may include people who are not programmers or analysts. RESULTS: We introduce shinyExprPortal, an R package that implements omics web portals with minimal coding effort. The portals allow exploration of transcriptomic or proteomic expression data and phenotypes, showcasing results of various types of analysis including differential expression, co-expression and pathways analysis. The integration with bioinformatics workflows enables researchers to focus on their results and share findings using interactive and publication-quality plots. AVAILABILITY AND IMPLEMENTATION: The shinyExprPortal package is available to download and install from CRAN and https://github.com/C4TB/shinyExprPortal.

Expression quantitative trait loci analysis in rheumatoid arthritis identifies tissue specific variants associated with severity and outcome.

OBJECTIVE: Genome-wide association studies have successfully identified more than 100 loci associated with susceptibility to rheumatoid arthritis (RA). However, our understanding of the functional effects of genetic variants in causing RA and their effects on disease severity and response to treatment remains limited. METHODS: In this study, we conducted expression quantitative trait locus (eQTL) analysis to dissect the link between genetic variants and gene expression comparing the disease tissue against blood using RNA-Sequencing of synovial biopsies (n=85) and blood samples (n=51) from treatment-naïve patients with RA from the Pathobiology of Early Arthritis Cohort. RESULTS: This identified 898 eQTL genes in synovium and genes loci in blood, with 232 genes in common to both synovium and blood, although notably many eQTL were tissue specific. Examining the HLA region, we uncovered a specific eQTL at HLA-DPB2 with the critical triad of single-nucleotide polymorphisms (SNPs) rs3128921 driving synovial HLA-DPB2 expression, and both rs3128921 and HLA-DPB2 gene expression correlating with clinical severity and increasing probability of the lympho-myeloid pathotype. CONCLUSIONS: This analysis highlights the need to explore functional consequences of genetic associations in disease tissue. HLA-DPB2 SNP rs3128921 could potentially be used to stratify patients to more aggressive treatment immediately at diagnosis.

The benefits and pitfalls of machine learning for biomarker discovery.

Prospects for the discovery of robust and reproducible biomarkers have improved considerably with the development of sensitive omics platforms that can enable measurement of biological molecules at an unprecedented scale. With technical barriers to success lowering, the challenge is now moving into the analytical domain. Genome-wide discovery presents a problem of scale and multiple testing as standard statistical methods struggle to distinguish signal from noise in increasingly complex biological systems. Machine learning and AI methods are good at finding answers in large datasets, but they have a tendency to overfit solutions. It may be possible to find a local answer or mechanism in a specific patient sample or small group of samples, but this may not generalise to wider patient populations due to the high likelihood of false discovery. The rise of explainable AI offers to improve the opportunity for true discovery by providing explanations for predictions that can be explored mechanistically before proceeding to costly and time-consuming validation studies. This review aims to introduce some of the basic concepts of machine learning and AI for biomarker discovery with a focus on post hoc explanation of predictions. To illustrate this, we consider how explainable AI has already been used successfully, and we explore a case study that applies AI to biomarker discovery in rheumatoid arthritis, demonstrating the accessibility of tools for AI and machine learning. We use this to illustrate and discuss some of the potential challenges and solutions that may enable AI to critically interrogate disease and response mechanisms.

The successes and challenges of harmonising juvenile idiopathic arthritis (JIA) datasets to create a large-scale JIA data resource.

BACKGROUND: CLUSTER is a UK consortium focussed on precision medicine research in JIA/JIA-Uveitis. As part of this programme, a large-scale JIA data resource was created by harmonizing and pooling existing real-world studies. Here we present challenges and progress towards creation of this unique large JIA dataset. METHODS: Four real-world studies contributed data; two clinical datasets of JIA patients starting first-line methotrexate (MTX) or tumour necrosis factor inhibitors (TNFi) were created. Variables were selected based on a previously developed core dataset, and encrypted NHS numbers were used to identify children contributing similar data across multiple studies. RESULTS: Of 7013 records (from 5435 individuals), 2882 (1304 individuals) represented the same child across studies. The final datasets contain 2899 (MTX) and 2401 (TNFi) unique patients; 1018 are in both datasets. Missingness ranged from 10 to 60% and was not improved through harmonisation. CONCLUSIONS: Combining data across studies has achieved dataset sizes rarely seen in JIA, invaluable to progressing research. Losing variable specificity and missingness, and their impact on future analyses requires further consideration.

The impact of expression vector position on transgene transcription allows for rational expression vector design in a targeted integration system.

Site-specific integration (SSI) technology has emerged as an effective approach by the pharmaceutical industry for the development of recombinant Chinese hamster ovary (CHO) cell lines. While SSI systems have been demonstrated to be effective for the development of CHO cell lines, they can be limiting in terms of both transgene expression and in the case of multi-specifics, the ability to generate the correct product of interest. To maximize the performance of Pfizer's dual SSI expression system for expressing monoclonal and multi-specific antibodies, we used a novel approach to investigate the positional effect of transgenes within expression vectors by engineering nucleotide polymorphisms (NP)s to use as biomarkers to track the level of transcript output from each expression vector position. We observed differences in transcript level for two different transgenes across all four expression vector positions interrogated. We then applied these learnings to rationally design expression vectors for six different mAbs and a multi-specific antibody. We showed enhanced productivity and optimal product quality when compared to a conventional expression vector topology. The learnings gained here can potentially aid in the determination of optimal vector topologies for several IgG-like multi-specific formats.

nestedcv: an R package for fast implementation of nested cross-validation with embedded feature selection designed for transcriptomics and high-dimensional data.

Motivation: Although machine learning models are commonly used in medical research, many analyses implement a simple partition into training data and hold-out test data, with cross-validation (CV) for tuning of model hyperparameters. Nested CV with embedded feature selection is especially suited to biomedical data where the sample size is frequently limited, but the number of predictors may be significantly larger (P ≫ n). Results: The nestedcv R package implements fully nested k × l-fold CV for lasso and elastic-net regularized linear models via the glmnet package and supports a large array of other machine learning models via the caret framework. Inner CV is used to tune models and outer CV is used to determine model performance without bias. Fast filter functions for feature selection are provided and the package ensures that filters are nested within the outer CV loop to avoid information leakage from performance test sets. Measurement of performance by outer CV is also used to implement Bayesian linear and logistic regression models using the horseshoe prior over parameters to encourage a sparse model and determine unbiased model accuracy. Availability and implementation: The R package nestedcv is available from CRAN: https://CRAN.R-project.org/package=nestedcv.

Local delivery of nitric oxide prevents endothelial dysfunction in periodontitis.

AIMS: Increased cardiovascular disease risk underlies elevated rates of mortality in individuals with periodontitis. A key characteristic of those with increased cardiovascular risk is endothelial dysfunction, a phenomenon synonymous with deficiencies of bioavailable nitric oxide (NO), and prominently expressed in patients with periodontitis. Also, inorganic nitrate can be reduced to NO in vivo to restore NO levels, leading us to hypothesise that its use may be beneficial in reducing periodontitis-associated endothelial dysfunction. Herein we sought to determine whether inorganic nitrate improves endothelial function in the setting of periodontitis and if so to determine the mechanisms underpinning any responses seen. METHODS AND RESULTS: Periodontitis was induced in mice by placement of a ligature for 14 days around the second molar. Treatment in vivo with potassium nitrate, either prior to or following establishment of experimental periodontitis, attenuated endothelial dysfunction, as determined by assessment of acetylcholine-induced relaxation of aortic rings, compared to control (potassium chloride treatment). These beneficial effects were associated with a suppression of vascular wall inflammatory pathways (assessed by quantitative-PCR), increases in the anti-inflammatory cytokine interleukin (IL)-10 and reduced tissue oxidative stress due to attenuation of xanthine oxidoreductase-dependent superoxide generation. In patients with periodontitis, plasma nitrite levels were not associated with endothelial function indicating dysfunction. CONCLUSION: Our results suggest that inorganic nitrate protects against, and can partially reverse pre-existing, periodontitis-induced endothelial dysfunction through restoration of nitrite and thus NO levels. This research highlights the potential of dietary nitrate as adjunct therapy to target the associated negative cardiovascular outcomes in patients with periodontitis.

Predicting post-contrast information from contrast agent free cardiac MRI using machine learning: Challenges and methods.

Objectives: Currently, administering contrast agents is necessary for accurately visualizing and quantifying presence, location, and extent of myocardial infarction (MI) with cardiac magnetic resonance (CMR). In this study, our objective is to investigate and analyze pre- and post-contrast CMR images with the goal of predicting post-contrast information using pre-contrast information only. We propose methods and identify challenges. Methods: The study population consists of 272 retrospectively selected CMR studies with diagnoses of MI (n = 108) and healthy controls (n = 164). We describe a pipeline for pre-processing this dataset for analysis. After data feature engineering, 722 cine short-axis (SAX) images and segmentation mask pairs were used for experimentation. This constitutes 506, 108, and 108 pairs for the training, validation, and testing sets, respectively. We use deep learning (DL) segmentation (UNet) and classification (ResNet50) models to discover the extent and location of the scar and classify between the ischemic cases and healthy cases (i.e., cases with no regional myocardial scar) from the pre-contrast cine SAX image frames, respectively. We then capture complex data patterns that represent subtle signal and functional changes in the cine SAX images due to MI using optical flow, rate of change of myocardial area, and radiomics data. We apply this dataset to explore two supervised learning methods, namely, the support vector machines (SVM) and the decision tree (DT) methods, to develop predictive models for classifying pre-contrast cine SAX images as being a case of MI or healthy. Results: Overall, for the UNet segmentation model, the performance based on the mean Dice score for the test set (n = 108) is 0.75 (±0.20) for the endocardium, 0.51 (±0.21) for the epicardium and 0.20 (±0.17) for the scar. For the classification task, the accuracy, F1 and precision scores of 0.68, 0.69, and 0.64, respectively, were achieved with the SVM model, and of 0.62, 0.63, and 0.72, respectively, with the DT model. Conclusion: We have presented some promising approaches involving DL, SVM, and DT methods in an attempt to accurately predict contrast information from non-contrast images. While our initial results are modest for this challenging task, this area of research still poses several open problems.

Interferon-α-mediated therapeutic resistance in early rheumatoid arthritis implicates epigenetic reprogramming.

OBJECTIVES: An interferon (IFN) gene signature (IGS) is present in approximately 50% of early, treatment naive rheumatoid arthritis (eRA) patients where it has been shown to negatively impact initial response to treatment. We wished to validate this effect and explore potential mechanisms of action. METHODS: In a multicentre inception cohort of eRA patients (n=191), we examined the whole blood IGS (MxA, IFI44L, OAS1, IFI6, ISG15) with reference to circulating IFN proteins, clinical outcomes and epigenetic influences on circulating CD19+ B and CD4+ T lymphocytes. RESULTS: We reproduced our previous findings demonstrating a raised baseline IGS. We additionally showed, for the first time, that the IGS in eRA reflects circulating IFN-α protein. Paired longitudinal analysis demonstrated a significant reduction between baseline and 6-month IGS and IFN-α levels (p<0.0001 for both). Despite this fall, a raised baseline IGS predicted worse 6-month clinical outcomes such as increased disease activity score (DAS-28, p=0.025) and lower likelihood of a good EULAR clinical response (p=0.034), which was independent of other conventional predictors of disease activity and clinical response. Molecular analysis of CD4+ T cells and CD19+ B cells demonstrated differentially methylated CPG sites and dysregulated expression of disease relevant genes, including PARP9, STAT1, and EPSTI1, associated with baseline IGS/IFNα levels. Differentially methylated CPG sites implicated altered transcription factor binding in B cells (GATA3, ETSI, NFATC2, EZH2) and T cells (p300, HIF1α). CONCLUSIONS: Our data suggest that, in eRA, IFN-α can cause a sustained, epigenetically mediated, pathogenic increase in lymphocyte activation and proliferation, and that the IGS is, therefore, a robust prognostic biomarker. Its persistent harmful effects provide a rationale for the initial therapeutic targeting of IFN-α in selected patients with eRA.

Rituximab versus tocilizumab in rheumatoid arthritis: synovial biopsy-based biomarker analysis of the phase 4 R4RA randomized trial.

Patients with rheumatoid arthritis (RA) receive highly targeted biologic therapies without previous knowledge of target expression levels in the diseased tissue. Approximately 40% of patients do not respond to individual biologic therapies and 5-20% are refractory to all. In a biopsy-based, precision-medicine, randomized clinical trial in RA (R4RA; n = 164), patients with low/absent synovial B cell molecular signature had a lower response to rituximab (anti-CD20 monoclonal antibody) compared with that to tocilizumab (anti-IL6R monoclonal antibody) although the exact mechanisms of response/nonresponse remain to be established. Here, in-depth histological/molecular analyses of R4RA synovial biopsies identify humoral immune response gene signatures associated with response to rituximab and tocilizumab, and a stromal/fibroblast signature in patients refractory to all medications. Post-treatment changes in synovial gene expression and cell infiltration highlighted divergent effects of rituximab and tocilizumab relating to differing response/nonresponse mechanisms. Using ten-by-tenfold nested cross-validation, we developed machine learning algorithms predictive of response to rituximab (area under the curve (AUC) = 0.74), tocilizumab (AUC = 0.68) and, notably, multidrug resistance (AUC = 0.69). This study supports the notion that disease endotypes, driven by diverse molecular pathology pathways in the diseased tissue, determine diverse clinical and treatment-response phenotypes. It also highlights the importance of integration of molecular pathology signatures into clinical algorithms to optimize the future use of existing medications and inform the development of new drugs for refractory patients.

GEOexplorer: a webserver for gene expression analysis and visualisation.

Gene Expression Omnibus (GEO) is a database repository hosting a substantial proportion of publicly available high throughput gene expression data. Gene expression analysis is a powerful tool to gain insight into the mechanisms and processes underlying the biological and phenotypic differences between sample groups. Despite the wide availability of gene expression datasets, their access, analysis, and integration are not trivial and require specific expertise and programming proficiency. We developed the GEOexplorer webserver to allow scientists to access, integrate and analyse gene expression datasets without requiring programming proficiency. Via its user-friendly graphic interface, users can easily apply GEOexplorer to perform interactive and reproducible gene expression analysis of microarray and RNA-seq datasets, while producing a wealth of interactive visualisations to facilitate data exploration and interpretation, and generating a range of publication ready figures. The webserver allows users to search and retrieve datasets from GEO as well as to upload user-generated data and combine and harmonise two datasets to perform joint analyses. GEOexplorer, available at https://geoexplorer.rosalind.kcl.ac.uk, provides a solution for performing interactive and reproducible analyses of microarray and RNA-seq gene expression data, empowering life scientists to perform exploratory data analysis and differential gene expression analysis on-the-fly without informatics proficiency.

Bayesian network analysis of panomic biological big data identifies the importance of triglyceride-rich LDL in atherosclerosis development.

Introduction: We sought to explore biomarkers of coronary atherosclerosis in an unbiased fashion. Methods: We analyzed 665 patients (mean ± SD age, 56 ± 11 years; 47% male) from the GLOBAL clinical study (NCT01738828). Cases were defined by the presence of any discernable atherosclerotic plaque based on comprehensive cardiac computed tomography (CT). De novo Bayesian networks built out of 37,000 molecular measurements and 99 conventional biomarkers per patient examined the potential causality of specific biomarkers. Results: Most highly ranked biomarkers by gradient boosting were interleukin-6, symmetric dimethylarginine, LDL-triglycerides [LDL-TG], apolipoprotein B48, palmitoleic acid, small dense LDL, alkaline phosphatase, and asymmetric dimethylarginine. In Bayesian analysis, LDL-TG was directly linked to atherosclerosis in over 95% of the ensembles. Genetic variants in the genomic region encoding hepatic lipase (LIPC) were associated with LIPC gene expression, LDL-TG levels and with atherosclerosis. Discussion: Triglyceride-rich LDL particles, which can now be routinely measured with a direct homogenous assay, may play an important role in atherosclerosis development. Clinical trial registration: GLOBAL clinical study (Genetic Loci and the Burden of Atherosclerotic Lesions); [https://clinicaltrials.gov/ct2/show/NCT01738828?term=NCT01738828&rank=1], identifier [NCT01738828].

Whole exome sequencing identifies deleterious rare variants in CCDC141 in familial self-limited delayed puberty.

Developmental abnormalities of the gonadotropin-releasing hormone (GnRH) neuronal network result in a range of conditions from idiopathic hypogonadotropic hypogonadism to self-limited delayed puberty. We aimed to discover important underlying regulators of self-limited delayed puberty through interrogation of GnRH pathways. Whole exome sequencing (WES) data consisting of 193 individuals, from 100 families with self-limited delayed puberty, was analysed using a virtual panel of genes related to GnRH development and function (n = 12). Five rare predicted deleterious variants in Coiled-Coil Domain Containing 141 (CCDC141) were identified in 21 individuals from 6 families (6% of the tested cohort). Homology modeling predicted all five variants to be deleterious. CCDC141 mutant proteins showed atypical subcellular localization associated with abnormal distribution of acetylated tubulin, and expression of mutants resulted in a significantly delayed cell migration, demonstrated in transfected HEK293 cells. These data identify mutations in CCDC141 as a frequent finding in patients with self-limited delayed puberty. The mis-localization of acetylated tubulin and reduced cell migration seen with mutant CCDC141 suggests a role of the CCDC141-microtubule axis in GnRH neuronal migration, with heterozygous defects potentially impacting the timing of puberty.

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.).

Enhanced NF-κB signaling in type-2 dendritic cells at baseline predicts non-response to adalimumab in psoriasis.

Biologic therapies have transformed the management of psoriasis, but clinical outcome is variable leaving an unmet clinical need for predictive biomarkers of response. Here we perform in-depth immunomonitoring of blood immune cells of 67 patients with psoriasis, before and during therapy with the anti-TNF drug adalimumab, to identify immune mediators of clinical response and evaluate their predictive value. Enhanced NF-κBp65 phosphorylation, induced by TNF and LPS in type-2 dendritic cells (DC) before therapy, significantly correlates with lack of clinical response after 12 weeks of treatment. The heightened NF-κB activation is linked to increased DC maturation in vitro and frequency of IL-17+ T cells in the blood of non-responders before therapy. Moreover, lesional skin of non-responders contains higher numbers of dermal DC expressing the maturation marker CD83 and producing IL-23, and increased numbers of IL-17+ T cells. Finally, we identify and clinically validate LPS-induced NF-κBp65 phosphorylation before therapy as a predictive biomarker of non-response to adalimumab, with 100% sensitivity and 90.1% specificity in an independent cohort. Our study uncovers important molecular and cellular mediators underpinning adalimumab mechanisms of action in psoriasis and we propose a blood biomarker for predicting clinical outcome.

Meeting Report: Psoriasis Stratification to Optimize Relevant Therapy Showcase.

A stratified medicine approach for the treatment of psoriasis promises greater certainty of clinical decision making through prediction of response on the basis of clinical, pharmacological, and -omics data from an individual patient. As yet, there is no predictive model for treatment response in routine clinical use for psoriasis. The Psoriasis Stratification to Optimise Relevant Therapy (PSORT) Consortium is a United Kingdom Medical Research Council‒funded, academic‒industrial stratified medicine consortium established with the objective of discovering the predictors and stratifiers of response of psoriasis to biologic therapies. A showcase meeting was convened and attended by 80 stakeholders at the Royal College of Physicians, London, United Kingdom on 18 November 2019. The purpose was to disseminate the research findings from the PSORT consortium discovered thus far. This report summarizes the presentations made on the day and the significant advances made by PSORT toward a stratified medicine approach to the management of psoriasis.

Enzymatic degradation of RNA causes widespread protein aggregation in cell and tissue lysates.

Most proteins in cell and tissue lysates are soluble. We show here that in lysate from human neurons, more than 1,300 proteins are maintained in a soluble and functional state by association with endogenous RNA, as degradation of RNA invariably leads to protein aggregation. The majority of these proteins lack conventional RNA-binding domains. Using synthetic oligonucleotides, we identify the importance of nucleic acid structure, with single-stranded pyrimidine-rich bulges or loops surrounded by double-stranded regions being particularly efficient in the maintenance of protein solubility. These experiments also identify an apparent one-to-one protein-nucleic acid stoichiometry. Furthermore, we show that protein aggregates isolated from brain tissue from Amyotrophic Lateral Sclerosis patients can be rendered soluble after refolding by both RNA and synthetic oligonucleotides. Together, these findings open new avenues for understanding the mechanism behind protein aggregation and shed light on how certain proteins remain soluble.

Investigation of genetically regulated gene expression and response to treatment in rheumatoid arthritis highlights an association between IL18RAP expression and treatment response.

OBJECTIVES: In this study, we sought to investigate whether there was any association between genetically regulated gene expression (as predicted using various reference panels) and anti-tumour necrosis factor (anti-TNF) treatment response (change in erythrocyte sedimentation rate (ESR)) using 3158 European ancestry patients with rheumatoid arthritis. METHODS: The genetically regulated portion of gene expression was estimated in the full cohort of 3158 subjects (as well as within a subcohort consisting of 1575 UK patients) using the PrediXcan software package with three different reference panels. Estimated expression was tested for association with anti-TNF treatment response. As a replication/validation experiment, we also investigated the correlation between change in ESR with measured gene expression at the Interleukin 18 Receptor Accessory Protein (IL18RAP) gene in whole blood and synovial tissue, using an independent replication data set of patients receiving conventional synthetic disease modifying anti-rheumatic drugs, with directly measured (via RNA sequencing) gene expression. RESULTS: We found that predicted expression of IL18RAP showed a consistent signal of association with treatment response across the reference panels. In our independent replication data set, IL18RAP expression in whole blood showed correlation with the change in ESR between baseline and follow-up (r=-0.35, p=0.0091). Change in ESR was also correlated with the expression of IL18RAP in synovial tissue (r=-0.28, p=0.02). CONCLUSION: Our results suggest that IL18RAP expression is worthy of further investigation as a potential predictor of treatment response in rheumatoid arthritis that is not specific to a particular drug type.

LGR4 deficiency results in delayed puberty through impaired Wnt/ß-catenin signaling.

The initiation of puberty is driven by an upsurge in hypothalamic gonadotropin-releasing hormone (GnRH) secretion. In turn, GnRH secretion upsurge depends on the development of a complex GnRH neuroendocrine network during embryonic life. Although delayed puberty (DP) affects up to 2% of the population, is highly heritable, and is associated with adverse health outcomes, the genes underlying DP remain largely unknown. We aimed to discover regulators by whole-exome sequencing of 160 individuals of 67 multigenerational families in our large, accurately phenotyped DP cohort. LGR4 was the only gene remaining after analysis that was significantly enriched for potentially pathogenic, rare variants in 6 probands. Expression analysis identified specific Lgr4 expression at the site of GnRH neuron development. LGR4 mutant proteins showed impaired Wnt/ß-catenin signaling, owing to defective protein expression, trafficking, and degradation. Mice deficient in Lgr4 had significantly delayed onset of puberty and fewer GnRH neurons compared with WT, whereas lgr4 knockdown in zebrafish embryos prevented formation and migration of GnRH neurons. Further, genetic lineage tracing showed strong Lgr4-mediated Wnt/ß-catenin signaling pathway activation during GnRH neuron development. In conclusion, our results show that LGR4 deficiency impairs Wnt/ß-catenin signaling with observed defects in GnRH neuron development, resulting in a DP phenotype.