Mutations in LNPK, Encoding the Endoplasmic Reticulum Junction Stabilizer Lunapark, Cause a Recessive Neurodevelopmental Syndrome

Rare variation in non-coding regions with evolutionary signatures contributes to autism spectrum disorder risk

Little is known about the role of non-coding regions in the etiology of autism spectrum disorder (ASD). We examined three classes of non-coding regions: human accelerated regions (HARs), which show signatures of positive selection in humans; experimentally validated neural VISTA enhancers (VEs); and conserved regions predicted to act as neural enhancers (CNEs). Targeted and whole-genome analysis of >16,600 samples and >4,900 ASD probands revealed that likely recessive, rare, inherited variants in HARs, VEs, and CNEs substantially contribute to ASD risk in probands whose parents share ancestry, which enriches for recessive contributions, but modestly contribute, if at all, in simplex family structures. We identified multiple patient variants in HARs near IL1RAPL1 and in VEs near OTX1 and SIM1 and showed that they change enhancer activity. Our results implicate both human-evolved and evolutionarily conserved non-coding regions in ASD risk and suggest potential mechanisms of how regulatory changes can modulate social behavior.

Fam134c and Fam134b shape axonal endoplasmic reticulum architecture in vivo

Endoplasmic reticulum (ER) remodeling is vital for cellular organization. ER-phagy, a selective autophagy targeting ER, plays an important role in maintaining ER morphology and function. The FAM134 protein family, including FAM134A, FAM134B, and FAM134C, mediates ER-phagy. While FAM134B mutations are linked to hereditary sensory and autonomic neuropathy in humans, the physiological role of the other FAM134 proteins remains unknown. To address this, we investigate the roles of FAM134 proteins using single and combined knockouts (KOs) in mice. Single KOs in young mice show no major phenotypes; however, combined Fam134b and Fam134c deletion (Fam134b/cdKO), but not the combination including Fam134a deletion, leads to rapid neuromuscular and somatosensory degeneration, resulting in premature death. Fam134b/cdKO mice show rapid loss of motor and sensory axons in the peripheral nervous system. Long axons from Fam134b/cdKO mice exhibit expanded tubular ER with a transverse ladder-like appearance, whereas no obvious abnormalities are present in cortical ER. Our study unveils the critical roles of FAM134C and FAM134B in the formation of tubular ER network in axons of both motor and sensory neurons.

RREB1 regulates neuronal proteostasis and the microtubule network

Transcription factors play vital roles in neuron development; however, little is known about the role of these proteins in maintaining neuronal homeostasis. Here, we show that the transcription factor RREB1 (Ras-responsive element-binding protein 1) is essential for neuron survival in the mammalian brain. A spontaneous mouse mutation causing loss of a nervous system-enriched Rreb1 transcript is associated with progressive loss of cerebellar Purkinje cells and ataxia. Analysis of chromatin immunoprecipitation and sequencing, along with RNA sequencing data revealed dysregulation of RREB1 targets associated with the microtubule cytoskeleton. In agreement with the known role of microtubules in dendritic development, dendritic complexity was disrupted in Rreb1-deficient neurons. Analysis of sequencing data also suggested that RREB1 plays a role in the endomembrane system. Mutant Purkinje cells had fewer numbers of autophagosomes and lysosomes and contained P62- and ubiquitin-positive inclusions. Together, these studies demonstrate that RREB1 functions to maintain the microtubule network and proteostasis in mammalian neurons.

Toward a CRISPR understanding of gene function in human brain development

Genome-wide association studies pinpoint genetic risk factors for neurodevelopmental disorders (NDDs), but the next challenge is to understand the mechanisms through which these genes affect brain development. Two recent CRISPR screens in human brain organoids1,2 interrogate the function of risk genes for autism spectrum disorder and other NDDs.

Assembloid CRISPR screens reveal impact of disease genes in human neurodevelopment

The assembly of cortical circuits involves the generation and migration of interneurons from the ventral to the dorsal forebrain1-3, which has been challenging to study at inaccessible stages of late gestation and early postnatal human development4. Autism spectrum disorder and other neurodevelopmental disorders (NDDs) have been associated with abnormal cortical interneuron development5, but which of these NDD genes affect interneuron generation and migration, and how they mediate these effects remains unknown. We previously developed a platform to study interneuron development and migration in subpallial organoids and forebrain assembloids6. Here we integrate assembloids with CRISPR screening to investigate the involvement of 425 NDD genes in human interneuron development. The first screen aimed at interneuron generation revealed 13 candidate genes, including CSDE1 and SMAD4. We subsequently conducted an interneuron migration screen in more than 1,000 forebrain assembloids that identified 33 candidate genes, including cytoskeleton-related genes and the endoplasmic reticulum-related gene LNPK. We discovered that, during interneuron migration, the endoplasmic reticulum is displaced along the leading neuronal branch before nuclear translocation. LNPK deletion interfered with this endoplasmic reticulum displacement and resulted in abnormal migration. These results highlight the power of this CRISPR-assembloid platform to systematically map NDD genes onto human development and reveal disease mechanisms.

Lunapark deficiency leads to an autosomal recessive neurodevelopmental phenotype with a degenerative course, epilepsy and distinct brain anomalies

LNPK encodes a conserved membrane protein that stabilizes the junctions of the tubular endoplasmic reticulum network playing crucial roles in diverse biological functions. Recently, homozygous variants in LNPK were shown to cause a neurodevelopmental disorder (OMIM#618090) in four patients displaying developmental delay, epilepsy and nonspecific brain malformations including corpus callosum hypoplasia and variable impairment of cerebellum. We sought to delineate the molecular and phenotypic spectrum of LNPK-related disorder. Exome or genome sequencing was carried out in 11 families. Thorough clinical and neuroradiological evaluation was performed for all the affected individuals, including review of previously reported patients. We identified 12 distinct homozygous loss-of-function variants in 16 individuals presenting with moderate to profound developmental delay, cognitive impairment, regression, refractory epilepsy and a recognizable neuroimaging pattern consisting of corpus callosum hypoplasia and signal alterations of the forceps minor ('ear-of-the-lynx' sign), variably associated with substantia nigra signal alterations, mild brain atrophy, short midbrain and cerebellar hypoplasia/atrophy. In summary, we define the core phenotype of LNPK-related disorder and expand the list of neurological disorders presenting with the 'ear-of-the-lynx' sign suggesting a possible common underlying mechanism related to endoplasmic reticulum-phagy dysfunction.

Lysosomal release of amino acids at ER three-way junctions regulates transmembrane and secretory protein mRNA translation

One-third of the mammalian proteome is comprised of transmembrane and secretory proteins that are synthesized on endoplasmic reticulum (ER). Here, we investigate the spatial distribution and regulation of mRNAs encoding these membrane and secretory proteins (termed "secretome" mRNAs) through live cell, single molecule tracking to directly monitor the position and translation states of secretome mRNAs on ER and their relationship to other organelles. Notably, translation of secretome mRNAs occurred preferentially near lysosomes on ER marked by the ER junction-associated protein, Lunapark. Knockdown of Lunapark reduced the extent of secretome mRNA translation without affecting translation of other mRNAs. Less secretome mRNA translation also occurred when lysosome function was perturbed by raising lysosomal pH or inhibiting lysosomal proteases. Secretome mRNA translation near lysosomes was enhanced during amino acid deprivation. Addition of the integrated stress response inhibitor, ISRIB, reversed the translation inhibition seen in Lunapark knockdown cells, implying an eIF2 dependency. Altogether, these findings uncover a novel coordination between ER and lysosomes, in which local release of amino acids and other factors from ER-associated lysosomes patterns and regulates translation of mRNAs encoding secretory and membrane proteins.

Heteromeric clusters of ubiquitinated ER-shaping proteins drive ER-phagy

Membrane-shaping proteins characterized by reticulon homology domains play an important part in the dynamic remodelling of the endoplasmic reticulum (ER). An example of such a protein is FAM134B, which can bind LC3 proteins and mediate the degradation of ER sheets through selective autophagy (ER-phagy)¹. Mutations in FAM134B result in a neurodegenerative disorder in humans that mainly affects sensory and autonomic neurons². Here we report that ARL6IP1, another ER-shaping protein that contains a reticulon homology domain and is associated with sensory loss³, interacts with FAM134B and participates in the formation of heteromeric multi-protein clusters required for ER-phagy. Moreover, ubiquitination of ARL6IP1 promotes this process. Accordingly, disruption of Arl6ip1 in mice causes an expansion of ER sheets in sensory neurons that degenerate over time. Primary cells obtained from Arl6ip1-deficient mice or from patients display incomplete budding of ER membranes and severe impairment of ER-phagy flux. Therefore, we propose that the clustering of ubiquitinated ER-shaping proteins facilitates the dynamic remodelling of the ER during ER-phagy and is important for neuronal maintenance.

Understanding neural development and diseases using CRISPR screens in human pluripotent stem cell-derived cultures

The brain is arguably the most complex part of the human body in form and function. Much remains unclear about the molecular mechanisms that regulate its normal and pathological physiology. This lack of knowledge largely stems from the inaccessible nature of the human brain, and the limitation of animal models. As a result, brain disorders are difficult to understand and even more difficult to treat. Recent advances in generating human pluripotent stem cells (hPSCs)-derived 2-dimensional (2D) and 3-dimensional (3D) neural cultures have provided an accessible system to model the human brain. Breakthroughs in gene editing technologies such as CRISPR/Cas9 further elevate the hPSCs into a genetically tractable experimental system. Powerful genetic screens, previously reserved for model organisms and transformed cell lines, can now be performed in human neural cells. Combined with the rapidly expanding single-cell genomics toolkit, these technological advances culminate to create an unprecedented opportunity to study the human brain using functional genomics. This review will summarize the current progress of applying CRISPR-based genetic screens in hPSCs-derived 2D neural cultures and 3D brain organoids. We will also evaluate the key technologies involved and discuss their related experimental considerations and future applications.

Novel LNPK variant causes progressive cerebral atrophy: Expanding the clinical phenotype

The biallelic variations of the LNPK gene are associated with the "neurodevelopmental disorder with epilepsy and hypoplasia of the corpus callosum" phenotype [MIM:618090] in the Online Mendelian Inheritance In Men database, and so far, two families have been identified in the literature. A third family with novel clinical features, who bears a novel variant in LNPK (NM_030650.3: c.770delA, p.D257fs*31) is described in the present study. The coexistence of psychomotor regression and neurodegeneration in brain magnetic resonance imaging was found for the first time in the present study, thanks to the long-term follow-up data on the case, which contributed to the phenotypic and mutation spectrum by means of the novel variation.

Replication of gene polymorphisms associated with periodontitis-related traits in an elderly cohort: the Washington Heights/Inwood Community Aging Project Ancillary Study of Oral Health

AIM

We sought to replicate findings from published genome-wide association studies (GWAS), linking specific candidate gene loci with periodontitis-related clinical/microbial traits.

MATERIALS AND METHODS

In the published GWAS, a total of 2196 single nucleotide polymorphisms associated with periodontitis-related traits at a p ≤ 5 × 10^-6 and mapped to 136 gene loci. The replication cohort included 1124 individuals, 65-98 years old (67% female, 45% Hispanic, 30% Black, 23% White) with available genome-wide genotypes and full-mouth periodontal status. Microbial profiles using checkerboard DNA-DNA hybridization and 16SrRNA sequencing were available from 912 and 739 participants, respectively.

RESULTS

Using gene-specific p-values after linkage disequilibrium pruning, the following gene/phenotype associations replicated successfully: CLEC19A with edentulism and %teeth with pocket depth (PD) ≥4 mm; IL37, HPVC1, TRPS1, ABHD12B, LDLRAD4 (C180rF1), TGM3, and GRK5 with %teeth with PD ≥4 mm; DAB2IP with presence of PD ≥6 mm; KIAA1715(LNPK), ROBO2, RAB28, LINC01017, NELL1, LDLRAD4(C18orF1), and CRYBB2P1 with %teeth with clinical attachment level (CAL) ≥3 mm; RUNX2 and LAMA2 with %teeth with CAL ≥5 mm; and KIAA1715(LNPK) with high colonization by Aggregatibacter actinomycetemcomitans. In addition, CLEC19A, IQSEC1, and EMR1 associated with microbial abundance based on checkerboard data, LBP and NCR2 with abundance based on sequencing data, and NCR2 with microbial diversity based on sequencing data.

CONCLUSIONS

Several gene loci identified in published GWAS as associated with periodontitis-related phenotypes replicated successfully in an elderly cohort.

Inferring mammalian tissue-specific regulatory conservation by predicting tissue-specific differences in open chromatin

BACKGROUND

Evolutionary conservation is an invaluable tool for inferring functional significance in the genome, including regions that are crucial across many species and those that have undergone convergent evolution. Computational methods to test for sequence conservation are dominated by algorithms that examine the ability of one or more nucleotides to align across large evolutionary distances. While these nucleotide alignment-based approaches have proven powerful for protein-coding genes and some non-coding elements, they fail to capture conservation of many enhancers, distal regulatory elements that control spatial and temporal patterns of gene expression. The function of enhancers is governed by a complex, often tissue- and cell type-specific code that links combinations of transcription factor binding sites and other regulation-related sequence patterns to regulatory activity. Thus, function of orthologous enhancer regions can be conserved across large evolutionary distances, even when nucleotide turnover is high.

RESULTS

We present a new machine learning-based approach for evaluating enhancer conservation that leverages the combinatorial sequence code of enhancer activity rather than relying on the alignment of individual nucleotides. We first train a convolutional neural network model that can predict tissue-specific open chromatin, a proxy for enhancer activity, across mammals. Next, we apply that model to distinguish instances where the genome sequence would predict conserved function versus a loss of regulatory activity in that tissue. We present criteria for systematically evaluating model performance for this task and use them to demonstrate that our models accurately predict tissue-specific conservation and divergence in open chromatin between primate and rodent species, vastly out-performing leading nucleotide alignment-based approaches. We then apply our models to predict open chromatin at orthologs of brain and liver open chromatin regions across hundreds of mammals and find that brain enhancers associated with neuron activity have a stronger tendency than the general population to have predicted lineage-specific open chromatin.

CONCLUSION

The framework presented here provides a mechanism to annotate tissue-specific regulatory function across hundreds of genomes and to study enhancer evolution using predicted regulatory differences rather than nucleotide-level conservation measurements.

ER-phagy: mechanisms, regulation and diseases connected to the lysosomal clearance of the endoplasmic reticulum

ER-phagy (reticulo-phagy) defines the degradation of portions of the endoplasmic reticulum (ER) within lysosomes or vacuoles. It is part of the self-digestion (i.e., auto-phagic) programs recycling cytoplasmic material and organelles, which rapidly mobilize metabolites in cells confronted with nutrient shortage. Moreover, selective clearance of ER subdomains participates to the control of ER size and activity during ER stress, the re-establishment of ER homeostasis after ER stress resolution and the removal of ER parts, in which aberrant and potentially cytotoxic material has been segregated. ER-phagy relies on the individual and/or concerted activation of the ER-phagy receptors, ER peripheral or integral membrane proteins that share the presence of LC3/Atg8-binding motifs in their cytosolic domains. ER-phagy involves the physical separation of portions of the ER from the bulk ER network, and their delivery to the endolysosomal/vacuolar catabolic district. This last step is accomplished by a variety of mechanisms including macro-ER-phagy (in which ER fragments are sequestered by double-membrane autophagosomes that eventually fuse with lysosomes/vacuoles), micro-ER-phagy (in which ER fragments are directly engulfed by endosomes/lysosomes/vacuoles), or direct fusion of ER-derived vesicles with lysosomes/vacuoles. ER-phagy is dysfunctional in specific human diseases and its regulators are subverted by pathogens, highlighting its crucial role for cell and organism life.

Autism-Related Transcription Factors Underlying the Sex-Specific Effects of Prenatal Bisphenol A Exposure on Transcriptome-Interactome Profiles in the Offspring Prefrontal Cortex

Bisphenol A (BPA) is an environmental risk factor for autism spectrum disorder (ASD). BPA exposure dysregulates ASD-related genes in the hippocampus and neurological functions of offspring. However, whether prenatal BPA exposure has an impact on genes in the prefrontal cortex, another brain region highly implicated in ASD, and through what mechanisms have not been investigated. Here, we demonstrated that prenatal BPA exposure disrupts the transcriptome-interactome profiles of the prefrontal cortex of neonatal rats. Interestingly, the list of BPA-responsive genes was significantly enriched with known ASD candidate genes, as well as genes that were dysregulated in the postmortem brain tissues of ASD cases from multiple independent studies. Moreover, several differentially expressed genes in the offspring's prefrontal cortex were the targets of ASD-related transcription factors, including AR, ESR1, and RORA. The hypergeometric distribution analysis revealed that BPA may regulate the expression of such genes through these transcription factors in a sex-dependent manner. The molecular docking analysis of BPA and ASD-related transcription factors revealed novel potential targets of BPA, including RORA, SOX5, TCF4, and YY1. Our findings indicated that prenatal BPA exposure disrupts ASD-related genes in the offspring's prefrontal cortex and may increase the risk of ASD through sex-dependent molecular mechanisms, which should be investigated further.

Advances in ER-Phagy and Its Diseases Relevance

As an important form of selective autophagy in cells, ER-phagy (endoplasmic reticulum-selective autophagy), the autophagic degradation of endoplasmic reticulum (ER), degrades ER membranes and proteins to maintain cellular homeostasis. The relationship between ER-phagy and human diseases, including neurodegenerative disorders, cancer, and other metabolic diseases has been unveiled by extensive research in recent years. Starting with the catabolic process of ER-phagy and key mediators in this pathway, this paper reviews the advances in the mechanism of ER-phagy and its diseases relevance. We hope to provide some enlightenment for further study on ER-phagy and the development of novel therapeutic strategies for related diseases.

Machine learning compensates fold-change method and highlights oxidative phosphorylation in the brain transcriptome of Alzheimer's disease

Alzheimer's disease (AD) is a neurodegenerative disorder causing 70% of dementia cases. However, the mechanism of disease development is still elusive. Despite the availability of a wide range of biological data, a comprehensive understanding of AD's mechanism from machine learning (ML) is so far unrealized, majorly due to the lack of needed data density. To harness the AD mechanism's knowledge from the expression profiles of postmortem prefrontal cortex samples of 310 AD and 157 controls, we used seven predictive operators or combinations of RapidMiner Studio operators to establish predictive models from the input matrix and to assign a weight to each attribute. Besides, conventional fold-change methods were also applied as controls. The identified genes were further submitted to enrichment analysis for KEGG pathways. The average accuracy of ML models ranges from 86.30% to 91.22%. The overlap ratio of the identified genes between ML and conventional methods ranges from 19.7% to 21.3%. ML exclusively identified oxidative phosphorylation genes in the AD pathway. Our results highlighted the deficiency of oxidative phosphorylation in AD and suggest that ML should be considered as complementary to the conventional fold-change methods in transcriptome studies.

ER-Phagy, ER Homeostasis, and ER Quality Control: Implications for Disease

Lysosomal degradation of endoplasmic reticulum (ER) fragments by autophagy, termed ER-phagy or reticulophagy, occurs under normal as well as stress conditions. The recent discovery of multiple ER-phagy receptors has stimulated studies on the roles of ER-phagy. We discuss how the ER-phagy receptors and the cellular components that work with these receptors mediate two important functions: ER homeostasis and ER quality control. We highlight that ER-phagy plays an important role in alleviating ER expansion induced by ER stress, and acts as an alternative disposal pathway for misfolded proteins. We suggest that the latter function explains the emerging connection between ER-phagy and disease. Additional ER-phagy-associated functions and important unanswered questions are also discussed.

A Plasmodium homolog of ER tubule-forming proteins is required for parasite virulence

Reticulon and REEP family of proteins stabilize the high curvature of endoplasmic reticulum (ER) tubules. Plasmodium berghei Yop1 (PbYop1) is a REEP5 homolog in Plasmodium. Here, we characterize its function using a gene-knockout (Pbyop1∆). Pbyop1∆ asexual stage parasites display abnormal ER architecture and an enlarged digestive vacuole. The erythrocytic cycle of Pbyop1∆ parasites is severely attenuated and the incidence of experimental cerebral malaria is significantly decreased in Pbyop1∆-infected mice. Pbyop1∆ sporozoites have reduced speed, are slower to invade host cells but give rise to equal numbers of infected HepG2 cells, as WT sporozoites. We propose that PbYOP1's disruption may lead to defects in trafficking and secretion of a subset of proteins required for parasite development and invasion of erythrocytes. Furthermore, the maintenance of ER morphology in different parasite stages is likely to depend on different proteins.

A COPII subunit acts with an autophagy receptor to target endoplasmic reticulum for degradation

The COPII-cargo adaptor complex Lst1-Sec23 selectively sorts proteins into vesicles that bud from the endoplasmic reticulum (ER) and traffic to the Golgi. Improperly folded proteins are prevented from exiting the ER and are degraded. ER-phagy is an autophagic degradation pathway that uses ER-resident receptors. Working in yeast, we found an unexpected role for Lst1-Sec23 in ER-phagy that was independent from its function in secretion. Up-regulation of the stress-inducible ER-phagy receptor Atg40 induced the association of Lst1-Sec23 with Atg40 at distinct ER domains to package ER into autophagosomes. Lst1-mediated ER-phagy played a vital role in maintaining cellular homeostasis by preventing the accumulation of an aggregation-prone protein in the ER. Lst1 function appears to be conserved because its mammalian homolog, SEC24C, was also required for ER-phagy.

Axonal Endoplasmic Reticulum Dynamics and Its Roles in Neurodegeneration

The physical continuity of axons over long cellular distances poses challenges for their maintenance. One organelle that faces this challenge is endoplasmic reticulum (ER); unlike other intracellular organelles, this forms a physically continuous network throughout the cell, with a single membrane and a single lumen. In axons, ER is mainly smooth, forming a tubular network with occasional sheets or cisternae and low amounts of rough ER. It has many potential roles: lipid biosynthesis, glucose homeostasis, a Ca2+ store, protein export, and contacting and regulating other organelles. This tubular network structure is determined by ER-shaping proteins, mutations in some of which are causative for neurodegenerative disorders such as hereditary spastic paraplegia (HSP). While axonal ER shares many features with the tubular ER network in other contexts, these features must be adapted to the long and narrow dimensions of axons. ER appears to be physically continuous throughout axons, over distances that are enormous on a subcellular scale. It is therefore a potential channel for long-distance or regional communication within neurons, independent of action potentials or physical transport of cargos, but involving its physiological roles such as Ca2+ or organelle homeostasis. Despite its apparent stability, axonal ER is highly dynamic, showing features like anterograde and retrograde transport, potentially reflecting continuous fusion and breakage of the network. Here we discuss the transport processes that must contribute to this dynamic behavior of ER. We also discuss the model that these processes underpin a homeostatic process that ensures both enough ER to maintain continuity of the network and repair breaks in it, but not too much ER that might disrupt local cellular physiology. Finally, we discuss how failure of ER organization in axons could lead to axon degenerative diseases, and how a requirement for ER continuity could make distal axons most susceptible to degeneration in conditions that disrupt ER continuity.

Recurrent homozygous damaging mutation in TMX2, encoding a protein disulfide isomerase, in four families with microlissencephaly

BACKGROUND

Protein disulfide isomerase (PDI) proteins are part of the thioredoxin protein superfamily. PDIs are involved in the formation and rearrangement of disulfide bonds between cysteine residues during protein folding in the endoplasmic reticulum and are implicated in stress response pathways.

METHODS

Eight children from four consanguineous families residing in distinct geographies within the Middle East and Central Asia were recruited for study. All probands showed structurally similar microcephaly with lissencephaly (microlissencephaly) brain malformations. DNA samples from each family underwent whole exome sequencing, assessment for repeat expansions and confirmatory segregation analysis.

RESULTS

An identical homozygous variant in TMX2 (c.500G>A), encoding thioredoxin-related transmembrane protein 2, segregated with disease in all four families. This variant changed the last coding base of exon 6, and impacted mRNA stability. All patients presented with microlissencephaly, global developmental delay, intellectual disability and epilepsy. While TMX2 is an activator of cellular C9ORF72 repeat expansion toxicity, patients showed no evidence of C9ORF72 repeat expansions.

CONCLUSION

The TMX2 c.500G>A allele associates with recessive microlissencephaly, and patients show no evidence of C9ORF72 expansions. TMX2 is the first PDI implicated in a recessive disease, suggesting a protein isomerisation defect in microlissencephaly.

Expression-based decision tree model reveals distinct microRNA expression pattern in pediatric neuronal and mixed neuronal-glial tumors

Background

The understanding of the molecular biology of pediatric neuronal and mixed neuronal-glial brain tumors is still insufficient due to low frequency and heterogeneity of those lesions which comprise several subtypes presenting neuronal and/or neuronal-glial differentiation. Important is that the most frequent ganglioglioma (GG) and dysembryoplastic neuroepithelial tumor (DNET) showed limited number of detectable molecular alterations. In such cases analyses of additional genomic mechanisms seem to be the most promising. The aim of the study was to evaluate microRNA (miRNA) profiles in GGs, DNETs and pilocytic asytrocytomas (PA) and test the hypothesis of plausible miRNA connection with histopathological subtypes of particular pediatric glial and mixed glioneronal tumors.

Methods

The study was designed as the two-stage analysis. Microarray testing was performed with the use of the miRCURY LNA microRNA Array technology in 51 cases. Validation set comprised 107 samples used during confirmation of the profiling results by qPCR bioinformatic analysis.

Results

Microarray data was compared between the groups using an analysis of variance with the Benjamini-Hochberg procedure used to estimate false discovery rates. After filtration 782 miRNAs were eligible for further analysis. Based on the results of 10 × 10-fold cross-validation J48 algorithm was identified as the most resilient to overfitting. Pairwise comparison showed the DNETs to be the most divergent with the largest number of miRNAs differing from either of the two comparative groups. Validation of array analysis was performed for miRNAs used in the classification model: miR-155-5p, miR-4754, miR-4530, miR-628-3p, let-7b-3p, miR-4758-3p, miRPlus-A1086 and miR-891a-5p. Model developed on their expression measured by qPCR showed weighted AUC of 0.97 (95% CI for all classes ranging from 0.91 to 1.00). A computational analysis was used to identify mRNA targets for final set of selected miRNAs using miRWalk database. Among genomic targets of selected molecules ZBTB20, LCOR, PFKFB2, SYNJ2BP and TPD52 genes were noted.

Conclusions

Our data showed the existence of miRNAs which expression is specific for different histological types of tumors. miRNA expression analysis may be useful in in-depth molecular diagnostic process of the tumors and could elucidate their origins and molecular background.

Electronic supplementary material

The online version of this article (10.1186/s12885-019-5739-5) contains supplementary material, which is available to authorized users.