New computational biology tools for the systematic analysis of the structure and expression of human genes

From the late 1980s, the automation of sequencing techniques and the computer spread gave rise to a flourishing number of new molecular structures and sequences and to proliferation of new databases in which to store them. Here are presented three computational approaches able to analyse the massive amount of publicly avalilable data in order to answer to important biological questions. The first strategy studies the incorrect assignment of the first AUG codon in a messenger RNA (mRNA), due to the incomplete determination of its 5' end sequence. An extension of the mRNA 5' coding region was identified in 477 in human loci, out of all human known mRNAs analysed, using an automated expressed sequence tag (EST)-based approach. Proof-of-concept confirmation was obtained by in vitro cloning and sequencing for GNB2L1, QARS and TDP2 and the consequences for the functional studies are discussed. The second approach analyses the codon bias, the phenomenon in which distinct synonymous codons are used with different frequencies, and, following integration with a gene expression profile, estimates the total number of codons present across all the expressed mRNAs (named here "codonome value") in a given biological condition. Systematic analyses across different pathological and normal human tissues and multiple species shows a surprisingly tight correlation between the codon bias and the codonome bias. The third approach is useful to studies the expression of human autism spectrum disorder (ASD) implicated genes. ASD implicated genes sharing microRNA response elements (MREs) for the same microRNA are co-expressed in brain samples from healthy and ASD affected individuals. The different expression of a recently identified long non coding RNA which have four MREs for the same microRNA could disrupt the equilibrium in this network, but further analyses and experiments are needed.

Translationsregulation des Myelin basischen Proteins in Gliazellen

Im zentralen Nervensystem (ZNS) myelinisieren Oligodendrozyten neuronale Axone, indem sie ihre Zellfortsätze mehrfach um axonale Segmente wickeln. Die Ausbildung dieser multilamellaren Membranstapel ermöglicht eine saltatorische und damit rasche und energie-effiziente Erregungsleitung (Nave, 2010). Eine Schädigung des Myelins beeinträchtigt die Reizweiterleitung und führt zur Degeneration der Axone, wie es zum Beispiel bei der Multiplen Sklerose der Fall ist. Das Myelin basische Protein (MBP) ist ein Hauptbestandteil des Myelin und ist essentiell für die Kompaktierung der Myelinmembran (Wood et al., 1984). Die MBP mRNA wird in hnRNP A2 enthaltenen RNA Granulen in einem translations-inaktiven Zustand zu den distalen Fortsätzen transportiert. Vermittelt durch axonale Signale wird nach axo-glialem Kontakt die Translation von MBP ermöglicht (White et al., 2008). Der genaue Mechanismus der differentiellen Genregulation des MBP Proteins ist bisher nur unzureichend aufgeklärt. In der vorliegenden Arbeit konnte eine kleine regulatorische RNA (sncRNA) identifiziert werden, welche über die seed Region mit der MBP mRNA interagieren und die Translation regulieren kann. In primären Oligodendrozyten führt die Überexpression der sncRNA-715 zu reduzierten MBP Protein Mengen und die Blockierung der endogenen sncRNA-715 führt zu einer gesteigerten MBP Synthese. Interessanterweise korreliert während der Differenzierung der Oligodendrozyten in vitro und in vivo die Synthese des MBP Proteins invers mit der Expression der sncRNA-715. In Oligodendrozyten beeinflusst eine experimentell erhöhte sncRNA-715 Menge die Zellmorphologie und induziert Apoptose. Weiterhin ist sncRNA-715 in zytoplasmatischen granulären Strukturen lokalisiert und assoziiert mit MBP mRNA in hnRNP A2 Transport- Granula. Diese Ergebnisse lassen vermuten, dass sncRNA-715 ein Bestandteil der hnRNP A2 Granula sein könnte und dort spezifisch die Translation der MBP mRNA während des Lokalisationsprozesses inhibiert. In chronischen MS Läsionen sind Olig2+-Zellen zu finden. Obwohl die MBP mRNA in diesen Läsionen nachzuweisen ist, kann kein Protein synthetisiert werden. In dieser Arbeit konnte gezeigt werden, dass in diesen Läsionen die Expression der sncRNA-715 erhöht ist. SncRNA-715 könnte die Translation von MBP verhindern und folglich als Inhibitor der Remyelinisierung während des Krankheitsverlaufs fungieren. Schwann-Zellen sind die myelinisierenden Zellen im peripheren Nervensystem (PNS). Im Zuge der Myelinisierung wird die MBP mRNA in diesen Gliazellen ebenfalls in die distalen Fortsätze transportiert und dort lokal translatiert und in die Myelinmembran eingebaut (Trapp et al., 1987). Im Gegensatz zum ZNS ist im PNS nur wenig über den Transportmechanismus der mRNA bekannt (Masaki, 2012). Es ist es sehr wahrscheinlich, dass in Schwann-Zellen und Oligodendrozyten die Lokalisation und die translationale Hemmung der MBP mRNA ähnlichen Mechanismen unterliegen. In der vorliegenden Arbeit konnte gezeigt werden, dass hnRNP A2 und sncRNA-715 in Schwann-Zellen exprimiert werden und in zytoplasmatischen Granula-ähnlichen Strukturen lokalisiert sind. Während der Differenzierung dieser Gliazellen in vivo und in vitro korreliert die Expression der sncRNA-715 invers mit der Synthese des MBP Proteins. HnRNP A2 und sncRNA-715 scheinen in Schwann-Zellen assoziiert zu sein und könnten wie in Oligodendrozyten den Transport der MBP mRNA vermitteln.

Combination of expression levels of miR-21 and miR-126 is associated with cancer-specific survival in clear-cell renal cell carcinoma

BACKGROUND Renal cell carcinoma (RCC) is marked by high mortality rate. To date, no robust risk stratification by clinical or molecular prognosticators of cancer-specific survival (CSS) has been established for early stages. Transcriptional profiling of small non-coding RNA gene products (miRNAs) seems promising for prognostic stratification. The expression of miR-21 and miR-126 was analysed in a large cohort of RCC patients; a combined risk score (CRS)-model was constructed based on expression levels of both miRNAs. METHODS Expression of miR-21 and miR-126 was evaluated by qRT-PCR in tumour and adjacent non-neoplastic tissue in n = 139 clear cell RCC patients. Relation of miR-21 and miR-126 expression with various clinical parameters was assessed. Parameters were analysed by uni- and multivariate COX regression. A factor derived from the z-score resulting from the COX model was determined for both miRs separately and a combined risk score (CRS) was calculated multiplying the relative expression of miR-21 and miR-126 by this factor. The best fitting COX model was selected by relative goodness-of-fit with the Akaike information criterion (AIC). RESULTS RCC with and without miR-21 up- and miR-126 downregulation differed significantly in synchronous metastatic status and CSS. Upregulation of miR-21 and downregulation of miR-126 were independently prognostic. A combined risk score (CRS) based on the expression of both miRs showed high sensitivity and specificity in predicting CSS and prediction was independent from any other clinico-pathological parameter. Association of CRS with CSS was successfully validated in a testing cohort containing patients with high and low risk for progressive disease. CONCLUSIONS A combined expression level of miR-21 and miR-126 accurately predicted CSS in two independent RCC cohorts and seems feasible for clinical application in assessing prognosis.

tRNA-derived fragments target the small ribosomal subunit to fine-tune translation

Post-transcriptional cleavage of RNA molecules to generate smaller fragments is a widespread mechanism that enlarges the structural and functional complexity of cellular RNomes. In particular, fragments deriving from both precursor and mature tRNAs represent one of the rapidly growing classes of post-transcriptional RNA pieces. Importantly, these tRNA-derived fragments (tRFs) possess distinct expression patterns, abundance, cellular localizations, or biological roles compared with their parental tRNA molecules [1]. Here we present evidence that tRFs from the halophilic archaeon Haloferax volcanii directly bind to ribosomes. In a previous genomic screen for ribosome-associated small RNAs we have identified a 26 residue long fragment originating from the 5’ part of valine tRNA (Val-tRF) to be by far the most abundant tRF in H. volcanii [2]. The Val-tRF is processed in a stress-dependent manner and was found to primarily target the small ribosomal subunit in vitro and in vivo. Translational activity was markedly reduced in the presence of Val-tRF, while control RNA fragments of similar length did not show inhibition of protein biosynthesis. Crosslinking experiments and subsequent primer extension analyses revealed the Val-tRF interaction site to surround the mRNA path in the 30S subunit. In support of this, binding experiments demonstrated that Val-tRF does compete with mRNAs for ribosome binding. Therefore this tRF represents a ribosome-associated non-coding RNA (rancRNA) capable of regulating gene expression in H. volcanii under environmental stress conditions probably by fine-tuning the rate of protein production [3].

Systems Biology Approaches to Probe Gene Regulation in Bacteria

Mechanisms that allow pathogens to colonize the host are not the product of isolated genes, but instead emerge from the concerted operation of regulatory networks. Therefore, identifying components and the systemic behavior of networks is necessary to a better understanding of gene regulation and pathogenesis. To this end, I have developed systems biology approaches to study transcriptional and post-transcriptional gene regulation in bacteria, with an emphasis in the human pathogen Mycobacterium tuberculosis (Mtb). First, I developed a network response method to identify parts of the Mtb global transcriptional regulatory network utilized by the pathogen to counteract phagosomal stresses and survive within resting macrophages. As a result, the method unveiled transcriptional regulators and associated regulons utilized by Mtb to establish a successful infection of macrophages throughout the first 14 days of infection. Additionally, this network-based analysis identified the production of Fe-S proteins coupled to lipid metabolism through the alkane hydroxylase complex as a possible strategy employed by Mtb to survive in the host. Second, I developed a network inference method to infer the small non-coding RNA (sRNA) regulatory network in Mtb. The method identifies sRNA-mRNA interactions by integrating a priori knowledge of possible binding sites with structure-driven identification of binding sites. The reconstructed network was useful to predict functional roles for the multitude of sRNAs recently discovered in the pathogen, being that several sRNAs were postulated to be involved in virulence-related processes. Finally, I applied a combined experimental and computational approach to study post-transcriptional repression mediated by small non-coding RNAs in bacteria. Specifically, a probabilistic ranking methodology termed rank-conciliation was developed to infer sRNA-mRNA interactions based on multiple types of data. The method was shown to improve target prediction in Escherichia coli, and therefore is useful to prioritize candidate targets for experimental validation.

Análise transcriptômica de miRNAs em pacientes sob dupla terapia de antiagregação

A terapia antiagregante é comumente indicada na prevenção e tratamento de doenças cardiovasculares. A dupla antiagregação com clopidrogrel e ácido acetilsalicílico (AAS) tem sido frequentemente adotada em pacientes com Doença Arterial Coronariana (DAC), mas apresenta ineficácia em uma parcela significativa da população com genótipo de respondedores. Essa falha terapêutica nos leva a questionar se outros mecanismos moleculares podem estar influenciando na resposta a esses fármacos. Recentes estudos sugerem que pequenas sequências de RNA não codificantes denominadas microRNAs (miRNAs) podem estar fortemente relacionadas com resposta ao tratamento fármaco-terapêutico, controlando as proteínas envolvidas na farmacocinética e farmacodinâmica. Entretanto, os principais miRNAs que atuam na dinâmica da resposta medicamentosa ainda não foram bem definidos. O objetivo deste estudo foi avaliar o perfil de miRNAs no sangue total periférico, procurando melhor esclarecer os mecanismos envolvidos na resposta aos antiagregantes plaquetários AAS e clopidogrel. Para isso, selecionou-se pacientes com DAC, os quais apresentavam diferentes respostas à dupla terapia de antiagregação determinadas pelo teste de agregação plaquetária. Baseados nos fenótipos, os perfis de expressão de miRNAs foram comparados entre os valores da taxa de agregação categorizados em tercis (T) de resposta. O grupo T1 foi constituído de pacientes respondedores, o T2 de respondedores intermediários e o T3 de não respondedores. Os perfis de miRNAs foram obtidos após sequenciamento de última geração e os dados obtidos foram analisados pelo pacote Deseq2. Os resultados mostraram 18 miRNAs diferentemente expressos entre os dois tercis extremos. Dentre esses miRNAs, 10 deles apresentaram importantes alvos relacionados com vias de ativação e agregação plaquetária quando analisados pelo software Ingenuity®. Dos 10 miRNAs, 4 deles, os quais apresentaram-se menos expressos no sequenciamento, demonstraram os mesmos perfis de expressão quando analisados pela reação em cadeia pela polimerase quantitativa (qPCR): hsa-miR-423-3p, hsa-miR-744-5p, hsa-miR- 30a-5p e hsa-let-7g-5p. A partir das análises de predição de alvos, pôde-se observar que os quatro miRNAs, quando menos expressos simultaneamente, predizem ativação da agregação plaquetária. Além disso, os miRNAs hsa-miR- 423-5p, hsa-miR-744-5p e hsa-let-7g-5p mostraram correlação com o perfil lipídico dos pacientes que, por sua vez, apresentou influência nos valores de agregação compreendidos no T3 de resposta a ambos os medicamentos. Sendo assim, conclui-se que maiores taxas de agregação plaquetária podem estar indiretamente relacionadas com os padrões de expressão de hsa-miR- 423-3p, hsa-miR-744-5p e hsa-let-7g-5p. Sugere-se que a avaliação do perfil de expressão destes 3 miRNAs no sangue periférico de pacientes com DAC possa predizer resposta terapêutica inadequada ao AAS e ao clopidogrel

Segmenting eukaryotic genomes with the generalized Gibbs sampler

Eukaryotic genomes display segmental patterns of variation in various properties, including GC content and degree of evolutionary conservation. DNA segmentation algorithms are aimed at identifying statistically significant boundaries between such segments. Such algorithms may provide a means of discovering new classes of functional elements in eukaryotic genomes. This paper presents a model and an algorithm for Bayesian DNA segmentation and considers the feasibility of using it to segment whole eukaryotic genomes. The algorithm is tested on a range of simulated and real DNA sequences, and the following conclusions are drawn. Firstly, the algorithm correctly identifies non-segmented sequence, and can thus be used to reject the null hypothesis of uniformity in the property of interest. Secondly, estimates of the number and locations of change-points produced by the algorithm are robust to variations in algorithm parameters and initial starting conditions and correspond to real features in the data. Thirdly, the algorithm is successfully used to segment human chromosome 1 according to GC content, thus demonstrating the feasibility of Bayesian segmentation of eukaryotic genomes. The software described in this paper is available from the author's website (www.uq.edu.au/similar to uqjkeith/) or upon request to the author.

Novel expression of microRNAs in serum samples of Iraqi breast cancer women

Although a lot of hard work against cancer to reduces its spread but it still continues to kill with abandon. The need for a biomarker for cancer early detection becomes the most mind concentrated scientists. MicroRNAs the tiny non coding RNA molecules opened new path for the scientists to determine the cancer in its early stages. Expression of microRNAs profiles has been investigated to be involved in cancer development. Here we determined the expression of microRNAs in serum of Iraqi healthy volunteers and other women diagnosed with breast cancer. MicroRNAs expression has been determined by using real time qPCR and delta method has been used. Four of thirteen microRNAs were shown to be expressed in serum of Iraqi breast cancer women. Let-7a and miR-21 were shown to be significantly over expressed in serum of breast cancer compared with healthy serum volunteers (P= 0.022 and 0.026) respectively. While miR-26b and miR-429 found to be significantly down expressed in serum of breast cancer women (P= 0.0034 and 0.031) respectively. The result concluded that these expressed microRNAs in serum of breast cancer women could be used as a first indicator of breast cancer occurrence.

MicroRNA expression and function in neonatal hypoxic ischaemic encephalopathy

Hypoxic ischaemic encephalopathy (HIE) is a devastating neonatal condition which affects 2-3 per 1000 infants annually. The current gold standard of treatment - induced hypothermia, has the ability to reduce neonatal mortality and improve neonatal morbidity. However, to be effective it needs to be initiated within the therapeutic window which exists following initial insult until approximately 6 hours after birth. Current methods of assessment which are relied upon to identify infants with HIE are subjective and unreliable. To overcome this issue, an early and reliable biomarker of HIE severity must be identified. MicroRNA (miRNA) are a class of small non-coding RNA molecules which have potential as biomarkers of disease state and potential therapeutic targets. These tiny molecules can modulate gene expression by inhibiting translation of messenger RNA (mRNA) and as a result, can regulate protein synthesis. These miRNA are understood to be released into the circulation during cellular stress, where they are highly stable and relatively easy to quantify. Therefore, these miRNAs may be ideal candidates for biomarkers of HIE severity and may aid in directing the clinical management of these infants. By using both transcriptomic and proteomic approaches to analyse the expression of miRNAs and their potential targets in the umbilical cord blood, I have confirmed that infants with perinatal asphyxia and HIE have a significantly different UCB miRNA signature compared to UCB samples from healthy controls. Finally, I have identified and investigated 2 individual miRNAs; both of which show some potential as classifiers of HIE severity and predictors of long term outcome, particularly when coupled with their downstream targets. While this work will need to be validated and expanded in a new and larger cohort of infants, it suggests the potential of miRNA as biomarkers of neonatal pathological conditions such as HIE.

MiRNA 34a: a therapeutic target for castration-resistant prostate cancer

INTRODUCTION: Development of a therapy for bone metastases is of paramount importance for castration-resistant prostate cancer (CRPC). The osteomimetic properties of CRPC confer a propensity to metastasize to osseous sites. Micro-ribonucleic acid (miRNA) is non-coding RNA that acts as a post-transcriptional regulator of multiple proteins and associated pathways. Therefore identification of miRNAs could reveal a valid third generation therapy for CRPC. Areas covered: miR34a has been found to play an integral role in the progression of prostate cancer, particularly in the regulation of metastatic genes involved in migration, intravasation, extravasation, bone attachment and bone homeostasis. The correlation between miR34a down-regulation and metastatic progression has generated substantial interest in this field. Expert opinion: Examination of the evidence reveals that miR34a is an ideal target for gene therapy for metastatic CRPC. We also conclude that future studies should focus on the effects of miR34a upregulation in CRPC with respect to migration, translocation to bone micro-environment and osteomimetic phenotype development. The success of miR34a as a therapeutic is reliant on the development of appropriate delivery systems and targeting to the bone micro-environment. In tandem with any therapeutic studies, biomarker serum levels should also be ascertained as an indicator of successful miR34a delivery.

miRNAs as therapeutic agents in neurodegeneration, a pilot study

L’échec des différents essais cliniques souligne la nécessité de développer des nouvelles thérapies pour la maladie d’Alzheimer (MA), la cause la plus commune de démence. Les microARNs (miARNs) sont les ARNs non-codants les plus étudiés et ils jouent un rôle important dans la modulation de l’expression des gènes et de multiples voies de signalisation. Des études antérieures, dont celles de mon laboratoire d’accueil, ont permis de développer l’hypothèse que certains membres de la famille miR-15/107 (c.-à-d. miR-15ab, miR-16, miR-195, miR-424, and miR-497) pourraient être utilisés comme agents thérapeutiques dans MA. En effet, cette famille avait le potentiel de réguler de multiples gènes associés à MA, tels que la protéine précurseur de l’amyloïde (APP), la β-secrétase (BACE1), et la protéine Tau. Tel que démontré dans ce projet de thèse, j’ai choisi miR-16 comme cible thérapeutique potentielle pour MA parmi tous les membres de la famille. L’essai luciférase dans ce projet confirme que miR-16 peut réguler simultanément APP et BACE1, directement par une interaction avec la région non-codante en 3’ de l’ARNm). Notamment, nous observons aussi une réduction de la production des peptides amyloïdes et de la phosphorylation de Tau après une augmentation de miR-16 en cellule. J’ai ensuite validé mes résultats in vivo dans la souris en utilisant une méthode de livraison de miR-16 via une pompe osmotique implanté dans le cerveau. Dans ce cas, l’expression des protéines d’intérêts (APP, BACE1, Tau) a été mesurée par immunobuvardage et PCR à temps réel. Après validation, ces résultats ont été complémentés par une étude protéomique (iTRAQ) du tronc cérébral et de l’hippocampe, deux régions associées à la maladie. Ces données m’ont permis d’identifier d’autres protéines régulées par miR-16 in vivo, incluant α-Synucléine, Transferrine receptor1, et SRm300. Une autre observation intéressante : les voies régulées par miR-16 in vivo sont directement en lien avec le stress oxydatif et la neurodégénération. En résumé, ce travail démontre l’efficacité et la faisabilité d’utiliser un miARN comme outil thérapeutique pour la maladie d’Alzheimer. Ces résultats rentrent dans un cadre plus vaste de découvrir de nouvelles cibles pour MA, et en particulier la forme sporadique de la maladie qui représente plus de 95% de tous les cas. Évidemment, la découverte d’une molécule pouvant cibler simultanément les deux pathologies de la maladie (plaques amyloïdes et hyper phosphorylation de tau) est nouvelle et intéressante, et ce domaine de recherche ouvre la porte aux autres petits ARNs non-codants dans MA et les maladies neurodégénératives connexes.

Genetic control of development in the parasitic nematode Haemonchus contortus by microRNAs

The parasitic nematode Haemonchus contortus has a major impact on the welfare and economic sustainability of small ruminant farming throughout the world. Increasing drug resistance requires the development of novel therapeutic agents. To further this process, we examined the fundamental biology of development in H. contortus, specifically, the potential role of microRNAs (miRNAs). miRNAs are short, non-coding RNA molecules that negatively regulate gene expression. In the free-living nematode Caenorhabditis elegans, miRNAs regulate a variety of genes including those involved in development. This thesis describes the expression patterns, potential targets and possible functions of miRNAs in H. contortus throughout development.

Macrophages and CD4 T-cells in rheumatoid arthritis and their modulation by JAK inhibitors

Background: Rheumatoid arthritis (RA) is a chronic inflammatory arthritis that causes significant morbidity and mortality and has no cure. Although early treatment strategies and biologic therapies such as TNFα blocking antibodies have revolutionised treatment, there still remains considerable unmet need. JAK kinase inhibitors, which target multiple inflammatory cytokines, have shown efficacy in treating RA although their exact mechanism of action remains to be determined. Stratified medicine promises to deliver the right drug to the right patient at the right time by using predictive ‘omic biomarkers discovered using bioinformatic and “Big Data” techniques. Therefore, knowledge across the realms of clinical rheumatology, applied immunology, bioinformatics and data science is required to realise this goal. Aim: To use bioinformatic tools to analyse the transcriptome of CD14 macrophages derived from patients with inflammatory arthritis and define a JAK/STAT signature. Thereafter to investigate the role of JAK inhibition on inflammatory cytokine production in a macrophage cell contact activation assay. Finally, to investigate JAK inhibition, following RA synovial fluid stimulation of monocytes. Methods and Results: Using bioinformatic software such as limma from the Bioconductor repository, I determined that there was a JAK/STAT signature in synovial CD14 macrophages from patients with RA and this differed from psoriatic arthritis samples. JAK inhibition using a JAK1/3 inhibitor tofacitinib reduced TNFα production when macrophages were cell contact activated by cytokine stimulated CD4 T-cells. Other pro-inflammatory cytokines such as IL-6 and chemokines such as IP-10 were also reduced. RA synovial fluid failed to stimulate monocytes to phosphorylate STAT1, 3 or 6 but CD4 T-cells activated STAT3 with this stimulus. RNA sequencing of synovial fluid stimulated CD4 T-cells showed an upregulation of SOCS3, BCL6 and SBNO2, a gene associated with RA but with unknown function and tofacitinib reversed this. Conclusion: These studies demonstrate that tofacitinib is effective at reducing inflammatory mediator production in a macrophage cell contact assay and also affects soluble factor mediated stimulation of CD4 T-cells. This suggests that the effectiveness of JAK inhibition is due to inhibition of multiple cytokine pathways such as IL-6, IL-15 and interferon. RNA sequencing is a useful tool to identify non-coding RNA transcripts that are associated with synovial fluid stimulation and JAK inhibition but these require further validation. SBNO2, a gene that is associated with RA, may be biomarker of tofacitinib treatment but requires further investigation and validation in wider disease cohorts.

La perturbation du locus Nr2f1-K12 entraine une différenciation gliale précoce dans un nouveau modèle murin de mégacôlon aganglionnaire

La maladie de Hirschsprung est une affection congénitale de la motilité intestinale caractérisée par un segment aganglionnaire dans le côlon terminal. Un criblage génétique par mutation insertionnelle aléatoire chez la souris nous a permis d’identifier la lignée transgénique Spot dont les homozygotes souffrent de mégacôlon aganglionnaire. L’analyse d’intestins d’embryons mutants a révélé une baisse de prolifération et un délai de migration des cellules de la crête neurale entériques (CCNe) progénitrices dus à leur différenciation gliale précoce, entrainant un défaut de colonisation de l’intestin et une aganglionose du côlon. Le séquençage du génome Spot indique que le transgène s’est inséré à l’intérieur du locus K12-Nr2f1 sur le chromosome 13, une région dépourvue de gènes préalablement associés à la maladie, perturbant également une séquence non-codante très conservée dans l’évolution. K12 est un gène d’ARN long non codant (ARNlnc) et antisens du gène Nr2f1, lui-même impliqué dans la gliogénèse du système nerveux central. Le séquençage du transcriptome des CCN a montré une surexpression de Nr2f1 et des formes courtes de K12 chez Spot et des essais luciférase ont révélé l’activité répressive de l’élément conservé. Nous avons observé l’expression de K12 dans les CCNe et sa localisation subcellulaire dans des zones transcriptionnellement actives du noyau. Avec l’émergence des ARNlnc régulateurs, ces données nous permettent de pointer deux nouveaux gènes candidats associés à une différenciation gliale prématurée du SNE menant au mégacôlon aganglionnaire, en supposant que la régulation de Nr2f1 se fait par son antisens, K12.