Abstract Hirschsprung disease (HSCR) is a congenital neurodevelopmental disorder characterized by segmental aganglionosis due to impaired developmental processes of enteric neural crest cells (NCCs). Despite being the leading genetic cause of functional intestinal obstruction in early childhood, HSCR represents a paradigmatic challenge in precision medicine: its multifactorial etiology, complex gene-environment interactions and limited resolution of single-modality analyses have long hindered mechanistic understanding and therapeutic translation. Here, we applied an integrative multi⍰omics approach combining genetic, phenotypic, epigenomic and transcriptomic analyses of matched ganglionic and aganglionic formalin⍰fixed paraffin⍰embedded (FFPE) patient tissues, complemented by patient-specific in vitro models. Beyond established genetic contributors, our integrative approach reveals novel regulatory pathways predominantly affecting enteric NCC differentiation, with convergent evidence pointing to epigenetic dysregulation as a primary disease mechanism. Notably, we identified over 1,300 differentially methylated positions between ganglionic and aganglionic FFPE samples, with HAND2 emerging as a key candidate due to multiple hypermethylated sites and consistently reduced expression levels in aganglionic tissues and in vitro models, suggesting a potential role in HSCR pathophysiology. We propose that our multi-omics approach offers a powerful and comprehensive framework for dissecting disease mechanisms. Beyond advancing biological understanding, this strategy holds promise for paving the way for molecularly informed patient stratification and supporting the development of personalized treatment and postoperative management strategies.

Short stature homeobox (SHOX) and its close paralogue SHOX2 are important transcriptional regulators of developmental processes. To understand how they function both individually and as an ensemble, we examined their mutual interdependence using zebrafish as a model. As both genes maintain important roles in various tissues, we assessed their expression patterns in the whole embryo as well as separately in fins, heart and brain. Fibroblast growth factor receptor 3 (FGFR3) is a direct target of SHOX, and FGFR3 mutations cause achondroplasia, a condition treatable with C-type natriuretic peptide (CNP) analogues. Building on this knowledge, we investigated how altered shox and shox2 expression relates to fgfr3 expression and the natriuretic peptide family (nppa, nppb, nppc) in zebrafish. Here, we demonstrate that shox and shox2 mutually regulate each other’s expression in a tissue-specific manner. In brain and whole embryos their effects are additive and synergistic, while in fins and heart their combined actions become partially antagonistic. Distinct tissue-specific effects of shox and shox2 were also noted on natriuretic peptide gene activity. In the fin, shox knockdown reduces nppc transcript levels, whereas in heart tissue, deficiency of shox and shox2, particularly in the double knockdown, results in an increased expression of nppb, a known marker for cardiac stress. We also show that shox2 deficiency increases fgfr3 levels, which are reduced by CNP treatment. These results in zebrafish may have implications for other vertebrates, including humans with SHOX deficiency, who should be monitored for cardiac disease in later life, and those who have an inadequate response to growth hormone may profit from CNP therapy.

KCNQ1 encodes a voltage-gated potassium channel implicated in various peripheral and neurological disorders, yet its role during human neuronal development remains unclear. To investigate this, we generated KCNQ1 knockouts (KO) in human induced pluripotent stem cell lines and differentiated them into neural stem cells (NSCs) and cortical neurons. KCNQ1 -deficient NSCs showed impaired neurite outgrowth, linked to reduced cell adhesion and disrupted neural cell adhesion molecule (NCAM) signalling. This phenotype was reproduced in wild-type and heterozygous lines by pharmacological KCNQ1 inhibition. Whole transcriptome, proteome, and follow-up analyses revealed mitochondrial dysfunction in KO NSCs, including reduced mitochondrial copy number and ATP synthase expression. Additionally, evidence was obtained for an impairment of insulin signalling in NSCs and neurons, with diminished insulin receptor gene expression and perturbation of key downstream signalling pathways (RAS-MAPK, PI3K-AKT). In neurons, KCNQ1 loss resulted in decreased synaptic activity and a more immature gene expression profile. Overall, our work reveals a novel role for KCNQ1 in human neurodevelopment by regulating cell adhesion, mitochondrial function, and insulin signalling. This work increases our understanding of KCNQ1 function in neurons and its contribution to neurological phenotypes observed in patients with KCNQ1 -related diseases.

Protein glycosylation is an essential post-translational modification in all domains of life. Its impairment in humans can result in severe diseases named congenital disorders of glycosylation (CDGs). Most of the glycosyltransferases (GTs) responsible for proper glycosylation are polytopic membrane proteins that represent challenging targets in proteomics. We established a multiple reaction monitoring (MRM) assay to comprehensively quantify GTs involved in the processes of N-glycosylation and O- and C-mannosylation in the endoplasmic reticulum. High robustness was achieved by using an enriched membrane protein fraction of isotopically labeled HEK 293T cells as an internal protein standard. The analysis of primary skin fibroblasts from eight CDG type I patients with impaired ALG1, ALG2, and ALG11 genes, respectively, revealed a substantial reduction in the corresponding protein levels. The abundance of the other GTs, however, remained unchanged at the transcript and protein levels, indicating that there is no fail-safe mechanism for the early steps of glycosylation in the endoplasmic reticulum. The established MRM assay was shared with the scientific community via the commonly used open source Skyline software environment, including Skyline Batch for automated data analysis. We demonstrate that another research group could easily reproduce all analysis steps, even while using different LC-MS hardware.

Very tall people attract much attention and represent a clinically and genetically heterogenous group of individuals. Identifying the genetic etiology can provide important insights into the molecular mechanisms regulating linear growth. We studied a three-generation pedigree with five isolated (non-syndromic) tall members and one individual with normal stature by whole exome sequencing; the tallest man had a height of 211 cm. Six heterozygous gene variants predicted as damaging were shared among the four genetically related tall individuals and not present in a family member with normal height. To gain insight into the putative role of these candidate genes in bone growth, we assessed the transcriptome of murine growth plate by microarray and RNA Seq. Two (Ift140, Nav2) of the six genes were well-expressed in the growth plate. Nav2 (p-value 1.91E-62) as well as Ift140 (p-value of 2.98E-06) showed significant downregulation of gene expression between the proliferative and hypertrophic zone, suggesting that these genes may be involved in the regulation of chondrocyte proliferation and/or hypertrophic differentiation. IFT140, NAV2 and SCAF11 have also significantly associated with height in GWAS studies. Pathway and network analysis indicated functional connections between IFT140, NAV2 and SCAF11 and previously associated (tall) stature genes. Knockout of the all-trans retinoic acid responsive gene, neuron navigator 2 NAV2, in Xenopus supports its functional role as a growth promotor. Collectively, our data expand the spectrum of genes with a putative role in tall stature phenotypes and, among other genes, highlight NAV2 as an interesting gene to this phenotype.

Background Epithelial to mesenchymal transition (EMT) is a key process in carcinogenesis of head and neck squamous cell carcinoma (HNSCC), contributing to tumor invasiveness, distant metastasis, and recurrence. Exosomes are known mediators and regulators of EMT. Here, we analyze the impact of exosomes that were primed by conventional therapy on EMT modulation. Methods Plasmas of n = 22 HNSCC patients were collected before and after standard of care surgery and adjuvant or primary (chemo)radiotherapy. Exosomes were isolated by size exclusion chromatography. Upon co-incubation of exosomes with HNSCC cells, the cellular EMT profile was analyzed by flow cytometry and RT-qPCR. Wound healing assays were performed to evaluate migratory potential of exosome-treated cells. Results Reduction of total exosome protein after therapy and in vitro exosome induced EMT profiles were dependent on the type of treatment. Exosomal TFG-β and miRNA cargo were partly responsible for observed exosome induced EMT changes. Exosomes from recurrent patients induced higher tumor cell migration after therapy than exosomes from disease-free patients. Conclusions HNSCC patients’ exosomes from timepoints before and after therapy were able to confer therapy induced EMT modulation in vitro and have the potential to monitor the EMT process. Exosome induced changes in migratory potential emerged as discriminants of therapy outcome.

BackgroundExosomes contribute to immunosuppression in head and neck squamous cell carcinoma (HNSCC), a tumor entity which lacks specific tumor biomarkers. Plasma-derived exosomes from HNSCC patients correlate with clinical parameters and have potential as liquid biopsy. Here, we investigate the cargo and functional profile of saliva-derived exosomes from HNSCC patients and their potential as non-invasive biomarkers for disease detection and immunomodulation.MethodsExosomes were isolated from saliva of HNSCC patients (n = 21) and healthy donors (HD, n = 12) by differential ultracentrifugation. Surface values of immune checkpoints and tumor associated antigens on saliva-derived exosomes were analyzed by bead-based flow cytometry using CD63 capture. Upon co-incubation with saliva-derived exosomes, activity and proliferation of T cells were assessed by flow cytometry (CD69 expression, CFSE assay). Adenosine levels were measured by mass spectrometry after incubation of saliva-derived exosomes with exogenous ATP. miRNA profiling of saliva-derived exosomes was performed using the nCounter® SPRINT system.ResultsSaliva-derived, CD63-captured exosomes from HNSCC patients carried high amounts of CD44v3, PDL1 and CD39. Compared to plasma, saliva was rich in tumor-derived, CD44v3+ exosomes and poor in hematopoietic cell-derived, CD45+ exosomes. CD8+ T cell activity was attenuated by saliva-derived exosomes from HNSCC patients, while proliferation of CD4+ T cells was not affected. Further, saliva-derived exosomes produced high levels of immunosuppressive adenosine. 62 HD- and 31 HNSCC-exclusive miRNAs were identified. Samples were grouped in “Healthy” and “Cancer” based on their saliva-derived exosomal miRNA profile, which was further found to be involved in RAS/MAPK, NF-κB complex, Smad2/3, and IFN-α signaling.ConclusionsSaliva-derived exosomes from HNSCC patients were enriched in tumor-derived exosomes whose cargo and functional profile reflected an immunosuppressive TME. Surface values of CD44v3, PDL1 and CD39 on CD63-captured exosomes, adenosine production and the miRNA cargo of saliva-derived exosomes emerged as discriminators of disease and emphasized their potential as liquid biomarkers specific for HNSCC.

Clinical studies in glioblastoma and pancreatic carcinoma patients strongly support the further development of H-1 protoparvovirus (H-1PV)-based anticancer therapies. The identification of cellular factors involved in the H-1PV life cycle may provide the knowledge to improve H-1PV anticancer potential. Recently, we showed that sialylated laminins mediate H-1PV attachment at the cell membrane. In this study, we revealed that H-1PV also interacts at the cell surface with galectin-1 and uses this glycoprotein to enter cancer cells. Indeed, knockdown/out of LGALS1, the gene encoding galectin-1, strongly decreases the ability of H-1PV to infect and kill cancer cells. This ability is rescued by the re-introduction of LGALS1 into cancer cells. Pre-treatment with lactose, which is able to bind to galectins and modulate their cellular functions, decreased H-1PV infectivity in a dose dependent manner. In silico analysis reveals that LGALS1 is overexpressed in various tumours including glioblastoma and pancreatic carcinoma. We show by immunohistochemistry analysis of 122 glioblastoma biopsies that galectin-1 protein levels vary between tumours, with levels in recurrent glioblastoma higher than those in primary tumours or normal tissues. We also find a direct correlation between LGALS1 transcript levels and H-1PV oncolytic activity in 53 cancer cell lines from different tumour origins. Strikingly, the addition of purified galectin-1 sensitises poorly susceptible GBM cell lines to H-1PV killing activity by rescuing cell entry. Together, these findings demonstrate that galectin-1 is a crucial determinant of the H-1PV life cycle.

BACKGROUND:There is increasing evidence that Parkinson's disease (PD) might start in the gut, thus involving and compromising also the enteric nervous system (ENS). At the clinical onset of the disease the majority of dopaminergic neurons in the midbrain is already destroyed, so that the lack of early biomarkers for the disease represents a major challenge for developing timely treatment interventions. Here, we use a transgenic A30P-α-synuclein-overexpressing PD mouse model to identify appropriate candidate markers in the gut before hallmark symptoms begin to manifest.METHODS:Based on a gait analysis and striatal dopamine levels, we defined 2-month-old A30P mice as pre-symptomatic (psA30P), since they are not showing any motoric impairments of the skeletal neuromuscular system and no reduced dopamine levels, but an intestinal α-synuclein pathology. Mice at this particular age were further used to analyze functional and molecular alterations in both, the gastrointestinal tract and the ENS, to identify early pathological changes. We examined the gastrointestinal motility, the molecular composition of the ENS, as well as the expression of regulating miRNAs. Moreover, we applied A30P-α-synuclein challenges in vitro to simulate PD in the ENS.RESULTS:A retarded gut motility and early molecular dysregulations were found in the myenteric plexus of psA30P mice. We found that i.e. neurofilament light chain, vesicle-associated membrane protein 2 and calbindin 2, together with the miRNAs that regulate them, are significantly altered in the psA30P, thus representing potential biomarkers for early PD. Many of the dysregulated miRNAs found in the psA30P mice are reported to be changed in PD patients as well, either in blood, cerebrospinal fluid or brain tissue. Interestingly, the in vitro approaches delivered similar changes in the ENS cultures as seen in the transgenic animals, thus confirming the data from the mouse model.CONCLUSIONS:These findings provide an interesting and novel approach for the identification of appropriate biomarkers in men.

Congenital primary hypothyroidism (CH; OMIM 218700) is characterized by an impaired thyroid development, or dyshormonogenesis, and can lead to intellectual disability and growth retardation if untreated. Most of the children with congenital hypothyroidism present thyroid dysgenesis, a developmental anomaly of the thyroid. Various genes have been associated with thyroid dysgenesis, but all known genes together can only explain a small number of cases. To identify novel genetic causes for congenital hypothyroidism, we performed trio whole-exome sequencing in an affected newborn and his unaffected parents. A predicted damaging de novo missense mutation was identified in the ZBTB26 gene (Zinc Finger A and BTB Domain containing 26). An additional cohort screening of 156 individuals with congenital thyroid dysgenesis identified two additional ZBTB26 gene variants of unknown significance. To study the underlying disease mechanism, morpholino knock-down of zbtb26 in Xenopus laevis was carried out, which demonstrated significantly smaller thyroid anlagen in knock-down animals at tadpole stage. Marker genes expressed in thyroid tissue precursors also indicated a specific reduction in the Xenopus ortholog of human Paired-Box-Protein PAX8, a transcription factor required for thyroid development, which could be rescued by adding zbtb26. Pathway and network analysis indicated network links of ZBTB26 to PAX8 and other genes involved in thyroid genesis and function. GWAS associations of ZBTB26 were found with height. Together, our study added a novel genetic risk factor to the list of genes underlying congenital primary hypothyroidism and provides additional support that de novo mutations, together with inherited variants, might contribute to the genetic susceptibility to CH.

Human growth is a complex trait. A considerable number of gene defects have been shown to cause short stature, but there are only few examples of genetic causes of non-syndromic tall stature. Besides rare variants with large effects and common risk alleles with small effect size, oligogenic effects may contribute to this phenotype. Exome sequencing was carried out in a tall male (height 3.5 SDS) and his parents. Filtered damaging variants with high CADD scores were validated by Sanger sequencing in the trio and three other affected and one unaffected family members. Network analysis was carried out to assess links between the candidate genes, and the transcriptome of murine growth plate was analyzed by microarray as well as RNA Seq. Heterozygous gene variants in CEP104, CROCC, NEK1, TOM1L2, and TSTD2 predicted as damaging were found to be shared between the four tall family members. Three of the five genes (CEP104, CROCC, and NEK1) belong to the ciliary gene family. All genes are expressed in mouse growth plate. Pathway and network analyses indicated close functional connections. Together, these data expand the spectrum of genes with a role in linear growth and tall stature phenotypes.

BACKGROUND:Congenital primary hypothyroidism (CH) is the most common endocrine disorder in neonates.METHODS:To identify novel genes, we performed whole exome sequencing (WES) in 6 patients with CH due to thyroid dysgenesis (TD). The potential effects of the most relevant variants were analyzed using in silico prediction tools. The most promising candidate gene, transient receptor potential channel 4-associated protein (TRPC4AP), was sequenced in 179 further patients with TD. Expression of TRPC4AP in human thyroid was investigated using RT-PCR. Trpc4ap- functional analysis was performed in Xenopus laevis using Morpholino (MO) antisense oligomers.RESULTS:WES identified a likely damaging mutation in TRPC4AP leading to a de novo stop codon p.Q552*. Targeted sequencing of TRPC4AP demonstrated gene variants with predicted damaging potential in 5 patients resulting each in an amino acid exchange (p.P706S, p.F729L, p.S777C, and p.N229S). We demonstrated that TRPC4AP is expressed in human thyroid gland tissue. Using Xenopus laevis, we showed that the volume of the tadpole thyroid anlage was reduced by 20% in Trpc4ap MO knockdowns compared to controls and by 41% in "Clustered Regularly Interspaced Short Palindromic Repeats"/Cas9-mediated gene knockout experiments.DISCUSSION:A recognized interaction of TRPC4AP and the NF-kappa-B-essential-modulator encoded by IKBKG gene was identified by IPA analysis. IKBKG plays a role in activation of the NF-κB-signaling pathway and regulates genes involved in proliferation and survival of thyrocytes and expression of key enzymes of thyroid hormone synthesis.CONCLUSION:TRPC4AP was identified as a novel candidate gene in TD, but further studies are needed to validate its role in thyroid function.

BACKGROUND:The enteric nervous system (ENS), a complex network of neurons and glial cells, coordinates major gastrointestinal functions. Impaired development or secondary aberrations cause severe enteric neuropathies. Neural crest-derived stem cells as well as enteric neuronal progenitor cells, which form enteric neurospheres, represent a promising tool to unravel molecular pathomechanisms and to develop novel therapy options. However, so far little is known about the detailed cellular composition and the proportional distribution of enteric neurospheres. Comprehensive knowledge will not only be essential for basic research but also for prospective cell replacement therapies to restore or to improve enteric neuronal dysfunction.METHODS:Human enteric neurospheres were generated from three individuals with varying age. For detailed molecular characterization, nCounter target gene expression analyses focusing on stem, progenitor, neuronal, glial, muscular, and epithelial cell markers were performed. Corresponding archived paraffin-embedded individuals' specimens were analyzed accordingly.KEY RESULTS:Our data revealed a remarkable molecular complexity of enteric neurospheres and archived specimens. Amongst the expression of multipotent stem cell, progenitor cell, neuronal, glial, muscle and epithelial cell markers, moderate levels for the pluripotency marker POU5F1 were observed. Furthermore, besides the interindividual variability, we identified highly distinct intraindividual expression profiles.CONCLUSIONS & INFERENCES:Our results emphasize the assessment of molecular signatures to be essential for standardized use, optimization of experimental approaches, and elimination of potential risk factors, as the formation of tumors. Our study pipeline may serve as a blueprint implemented into the characterization procedure of enteric neurospheres for various future applications.

ATP-binding cassette (ABC) transporters are the members of the efflux pumps that are responsible for the removal of cytotoxic substances by active transport. ABCB11, the bile salt efflux pump of hepatocytes, coordinates cellular excretion of numerous conjugated bile salts into the bile canaliculi, whereas ABCB4 acts as an ATP-dependent floppase translocating phosphatidylcholine from the inner to the outer leaflet of the bile canalicular membrane. Loss of functional ABCB11 and ABCB4 proteins causes early-onset refractory cholestasis or cholangiopathy. In this study, we investigated the expression and localization pattern of ABCB11 and ABCB4 using immunohistochemistry and RNA profiling in liver samples from patients with different types and stages of chronic cholestatic liver disease, with emphasis on primary sclerosing cholangitis (PSC), compared to a variety of cholestatic and noncholestatic hepatopathies. Therefore, ABCB11 and ABCB4 expressions were investigated on formalin-fixed and paraffin-embedded (FFPE) material in a patient cohort of total 43 patients with or without cholestatic liver diseases, on protein level using immunohistochemistry and on RNA level using nanoString technology. Intriguingly, our results demonstrated increased expression of ABCB11 and ABCB4 on protein as well as RNA level in PSC, and the expression pattern correlated with disease progression. We concluded from our study that patients with PSC demonstrate altered expression levels and pattern of ABCB11 and ABCB4 which correlated with disease progression; thereby, ABCB11 and ABCB4 analysis may be a useful tool for assessment of disease stages in PSC.

Height is a complex quantitative trait with a high heritability. Short stature is diagnosed when height is significantly below the average of the general population for that person's age and sex. We have recently found that the retinoic acid degrading enzyme CYP26C1 modifies SHOX deficiency phenotypes toward more severe clinical manifestations. Here, we asked whether damaging variants in CYP26C1 alone could lead to short stature. We performed exome and Sanger sequencing to analyze 856 individuals with short stature where SHOX deficiency was previously excluded. Three different damaging missense variants and one splicing variant were identified in six independent individuals; the functional significance of the identified variants was tested in vitro or in vivo using zebrafish as a model. The genetic and functional data reported here indicate that CYP26C1 represents a novel gene underlying growth disorders and that damaging variants in the absence of SHOX variants can lead to short stature.

Objective Micro-RNAs (miRNAs) play a crucial role in controlling intestinal epithelial barrier function partly by modulating the expression of tight junction (TJ) proteins. We have previously shown differential messenger RNA (mRNA) expression correlated with ultrastructural abnormalities of the epithelial barrier in patients with diarrhoea-predominant IBS (IBS-D). However, the participation of miRNAs in these differential mRNA-associated findings remains to be established. Our aims were (1) to identify miRNAs differentially expressed in the small bowel mucosa of patients with IBS-D and (2) to explore putative target genes specifically involved in epithelial barrier function that are controlled by specific dysregulated IBS-D miRNAs. Design Healthy controls and patients meeting Rome III IBS-D criteria were studied. Intestinal tissue samples were analysed to identify potential candidates by: (a) miRNA-mRNA profiling; (b) miRNA-mRNA pairing analysis to assess the co-expression profile of miRNA-mRNA pairs; (c) pathway analysis and upstream regulator identification; (d) miRNA and target mRNA validation. Candidate miRNA-mRNA pairs were functionally assessed in intestinal epithelial cells. Results IBS-D samples showed distinct miRNA and mRNA profiles compared with healthy controls. TJ signalling was associated with the IBS-D transcriptional profile. Further validation of selected genes showed consistent upregulation in 75% of genes involved in epithelial barrier function. Bioinformatic analysis of putative miRNA binding sites identified hsa-miR-125b-5p and hsa-miR-16 as regulating expression of the TJ genes CGN (cingulin) and CLDN2 (claudin-2), respectively. Consistently, protein expression of CGN and CLDN2 was upregulated in IBS-D, while the respective targeting miRNAs were downregulated. In addition, bowel dysfunction, perceived stress and depression and number of mast cells correlated with the expression of hsa-miR-125b-5p and hsa-miR-16 and their respective target proteins. Conclusions Modulation of the intestinal epithelial barrier function in IBS-D involves both transcriptional and post-transcriptional mechanisms. These molecular mechanisms include miRNAs as master regulators in controlling the expression of TJ proteins and are associated with major clinical symptoms.