Filaggrin loss-of-function variants are associated with atopic dermatitis phenotypes in a diverse, early-life prospective cohort.

Loss-of-function (LoF) variants in the filaggrin (FLG) gene are the strongest known genetic risk factor for atopic dermatitis (AD), but the impact of these variants on AD outcomes is poorly understood. We comprehensively identified genetic variants through targeted region sequencing of FLG in children participating in the Mechanisms of Progression of Atopic Dermatitis to Asthma in Children cohort. Twenty FLG LoF variants were identified, including 1 novel variant and 9 variants not previously associated with AD. FLG LoF variants were found in the cohort. Among these children, the presence of 1 or more FLG LoF variants was associated with moderate/severe AD compared with those with mild AD. Children with FLG LoF variants had a higher SCORing for Atopic Dermatitis (SCORAD) and higher likelihood of food allergy within the first 2.5 years of life. LoF variants were associated with higher transepidermal water loss (TEWL) in both lesional and nonlesional skin. Collectively, our study identifies established and potentially novel AD-associated FLG LoF variants and associates FLG LoF variants with higher TEWL in lesional and nonlesional skin.

VExD: a curated resource for human gene expression alterations following viral infection.

Much of the host antiviral response is mediated through changes to host gene expression levels. Likewise, viruses induce changes to host gene expression levels in order to promote the viral life cycle and evade the host immune system. However, there is no resource that specifically collects human gene expression levels pre- and post-virus infection. Further, public gene expression repositories do not contain enough specialized metadata to easily find relevant experiments. Here, we present the Virus Expression Database (VExD), a freely available website and database, that collects human gene expression datasets in response to viral infection. VExD contains ∼8,000 uniformly processed samples obtained from 289 studies examining 51 distinct human viruses. We show that the VExD processing pipeline captures known antiviral responses in the form of interferon-stimulated genes. We further show that the datasets collected in VExD can be used to quickly identify supporting data for experiments performed in human cells or model organisms. VExD is freely available at https://vexd.cchmc.org/.

Single cell transcriptomic analysis of HPV16-infected epithelium identifies a keratinocyte subpopulation implicated in cancer.

Persistent HPV16 infection is a major cause of the global cancer burden. The viral life cycle is dependent on the differentiation program of stratified squamous epithelium, but the landscape of keratinocyte subpopulations which support distinct phases of the viral life cycle has yet to be elucidated. Here, single cell RNA sequencing of HPV16 infected compared to uninfected organoids identifies twelve distinct keratinocyte populations, with a subset mapped to reconstruct their respective 3D geography in stratified squamous epithelium. Instead of conventional terminally differentiated cells, an HPV-reprogrammed keratinocyte subpopulation (HIDDEN cells) forms the surface compartment and requires overexpression of the ELF3/ESE-1 transcription factor. HIDDEN cells are detected throughout stages of human carcinogenesis including primary human cervical intraepithelial neoplasias and HPV positive head and neck cancers, and a possible role in promoting viral carcinogenesis is supported by TCGA analyses. Single cell transcriptome information on HPV-infected versus uninfected epithelium will enable broader studies of the role of individual keratinocyte subpopulations in tumor virus infection and cancer evolution.

Characteristics of children ≤36 months of age with DIPG: A report from the international DIPG registry.

BACKGROUND: Children ≤36 months with diffuse intrinsic pontine glioma (DIPG) have increased long-term survival (LTS, overall survival (OS) ≥24 months). Understanding distinguishing characteristics in this population is critical to improving outcomes. METHODS: Patients ≤36 months at diagnosis enrolled on the International DIPG Registry (IDIPGR) with central imaging confirmation were included. Presentation, clinical course, imaging, pathology and molecular findings were analyzed. RESULTS: Among 1183 patients in IDIPGR, 40 were eligible (median age: 29 months). Median OS was 15 months. Twelve patients (30%) were LTS, 3 (7.5%) very long-term survivors ≥5 years. Among 8 untreated patients, median OS was 2 months. Patients enrolled in the registry but excluded from our study by central radiology review or tissue diagnosis had median OS of 7 months. All but 1 LTS received radiation. Among 32 treated patients, 1-, 2-, 3-, and 5-year OS rates were 68.8%, 31.2%, 15.6% and 12.5%, respectively. LTS had longer duration of presenting symptoms (P = .018). No imaging features were predictive of outcome. Tissue and genomic data were available in 18 (45%) and 10 patients, respectively. Among 9 with known H3K27M status, 6 had a mutation. CONCLUSIONS: Children ≤36 months demonstrated significantly more LTS, with an improved median OS of 15 months; 92% of LTS received radiation. Median OS in untreated children was 2 months, compared to 17 months for treated children. LTS had longer duration of symptoms. Excluded patients demonstrated a lower OS, contradicting the hypothesis that children ≤36 months with DIPG show improved outcomes due to misdiagnosis.

Single-cell RNA-Seq of human esophageal epithelium in homeostasis and allergic inflammation.

Inflammation of the esophageal epithelium is a hallmark of eosinophilic esophagitis (EoE), an emerging chronic allergic disease. Herein, we probed human esophageal epithelial cells at single-cell resolution during homeostasis and EoE. During allergic inflammation, the epithelial differentiation program was blocked, leading to loss of KRT6hi differentiated populations and expansion of TOP2hi proliferating, DSPhi transitioning, and SERPINB3hi transitioning populations; however, there was stability of the stem cell-enriched PDPNhi basal epithelial compartment. This differentiation program blockade was associated with dysregulation of transcription factors, including nuclear receptor signalers, in the most differentiated epithelial cells and altered NOTCH-related cell-to-cell communication. Each epithelial population expressed genes with allergic disease risk variants, supporting their functional interplay. The esophageal epithelium differed notably between EoE in histologic remission and controls, indicating that remission is a transitory state poised to relapse. Collectively, our data uncover the dynamic nature of the inflamed human esophageal epithelium and provide a framework to better understand esophageal health and disease.

The genetic architecture of pediatric cardiomyopathy.

To understand the genetic contribution to primary pediatric cardiomyopathy, we performed exome sequencing in a large cohort of 528 children with cardiomyopathy. Using clinical interpretation guidelines and targeting genes implicated in cardiomyopathy, we identified a genetic cause in 32% of affected individuals. Cardiomyopathy sub-phenotypes differed by ancestry, age at diagnosis, and family history. Infants < 1 year were less likely to have a molecular diagnosis (p < 0.001). Using a discovery set of 1,703 candidate genes and informatic tools, we identified rare and damaging variants in 56% of affected individuals. We see an excess burden of damaging variants in affected individuals as compared to two independent control sets, 1000 Genomes Project (p < 0.001) and SPARK parental controls (p < 1 × 10-16). Cardiomyopathy variant burden remained enriched when stratified by ancestry, variant type, and sub-phenotype, emphasizing the importance of understanding the contribution of these factors to genetic architecture. Enrichment in this discovery candidate gene set suggests multigenic mechanisms underlie sub-phenotype-specific causes and presentations of cardiomyopathy. These results identify important information about the genetic architecture of pediatric cardiomyopathy and support recommendations for clinical genetic testing in children while illustrating differences in genetic architecture by age, ancestry, and sub-phenotype and providing rationale for larger studies to investigate multigenic contributions.

Characteristics of patients ≥10 years of age with diffuse intrinsic pontine glioma: a report from the International DIPG/DMG Registry.

BACKGROUND: Diffuse intrinsic pontine gliomas (DIPG) generally occur in young school-age children, although can occur in adolescents and young adults. The purpose of this study was to describe clinical, radiological, pathologic, and molecular characteristics in patients ≥10 years of age with DIPG enrolled in the International DIPG Registry (IDIPGR). METHODS: Patients ≥10 years of age at diagnosis enrolled in the IDIPGR with imaging confirmed DIPG diagnosis were included. The primary outcome was overall survival (OS) categorized as long-term survivors (LTS) (≥24 months) or short-term survivors (STS) (<24 months). RESULTS: Among 1010 patients, 208 (21%) were ≥10 years of age at diagnosis; 152 were eligible with a median age of 12 years (range 10-26.8). Median OS was 13 (2-82) months. The 1-, 3-, and 5-year OS was 59.2%, 5.3%, and 3.3%, respectively. The 18/152 (11.8%) LTS were more likely to be older (P < .01) and present with longer symptom duration (P < .01). Biopsy and/or autopsy were performed in 50 (33%) patients; 77%, 61%, 33%, and 6% of patients tested had H3K27M (H3F3A or HIST1H3B), TP53, ATRX, and ACVR1 mutations/genome alterations, respectively. Two of 18 patients with IDH1 testing were IDH1-mutant and 1 was a LTS. The presence or absence of H3 alterations did not affect survival. CONCLUSION: Patients ≥10 years old with DIPG have a median survival of 13 months. LTS present with longer symptom duration and are likely to be older at presentation compared to STS. ATRX mutation rates were higher in this population than the general DIPG population.

Genetic Causes of Cardiomyopathy in Children: First Results From the Pediatric Cardiomyopathy Genes Study.

Background Pediatric cardiomyopathy is a genetically heterogeneous disease with substantial morbidity and mortality. Current guidelines recommend genetic testing in children with hypertrophic, dilated, or restrictive cardiomyopathy, but practice variations exist. Robust data on clinical testing practices and diagnostic yield in children are lacking. This study aimed to identify the genetic causes of cardiomyopathy in children and to investigate clinical genetic testing practices. Methods and Results Children with familial or idiopathic cardiomyopathy were enrolled from 14 institutions in North America. Probands underwent exome sequencing. Rare sequence variants in 37 known cardiomyopathy genes were assessed for pathogenicity using consensus clinical interpretation guidelines. Of the 152 enrolled probands, 41% had a family history of cardiomyopathy. Of 81 (53%) who had undergone clinical genetic testing for cardiomyopathy before enrollment, 39 (48%) had a positive result. Genetic testing rates varied from 0% to 97% between sites. A positive family history and hypertrophic cardiomyopathy subtype were associated with increased likelihood of genetic testing (P=0.005 and P=0.03, respectively). A molecular cause was identified in an additional 21% of the 63 children who did not undergo clinical testing, with positive results identified in both familial and idiopathic cases and across all phenotypic subtypes. Conclusions A definitive molecular genetic diagnosis can be made in a substantial proportion of children for whom the cause and heritable nature of their cardiomyopathy was previously unknown. Practice variations in genetic testing are great and should be reduced. Improvements can be made in comprehensive cardiac screening and predictive genetic testing in first-degree relatives. Overall, our results support use of routine genetic testing in cases of both familial and idiopathic cardiomyopathy. Registration URL: https://www.clinicaltrials.gov; Unique identifier: NCT01873963.

Mucosal Genomics Implicate Lymphocyte Activation and Lipid Metabolism in Refractory Environmental Enteric Dysfunction.

BACKGROUND & AIMS: Environmental enteric dysfunction (EED) limits the Sustainable Development Goals of improved childhood growth and survival. We applied mucosal genomics to advance our understanding of EED. METHODS: The Study of Environmental Enteropathy and Malnutrition (SEEM) followed 416 children from birth to 24 months in a rural district in Pakistan. Biomarkers were measured at 9 months and tested for association with growth at 24 months. The duodenal methylome and transcriptome were determined in 52 undernourished SEEM participants and 42 North American controls and patients with celiac disease. RESULTS: After accounting for growth at study entry, circulating insulin-like growth factor-1 (IGF-1) and ferritin predicted linear growth, whereas leptin correlated with future weight gain. The EED transcriptome exhibited suppression of antioxidant, detoxification, and lipid metabolism genes, and induction of anti-microbial response, interferon, and lymphocyte activation genes. Relative to celiac disease, suppression of antioxidant and detoxification genes and induction of antimicrobial response genes were EED-specific. At the epigenetic level, EED showed hyper-methylation of epithelial metabolism and barrier function genes, and hypo-methylation of immune response and cell proliferation genes. Duodenal coexpression modules showed association between lymphocyte proliferation and epithelial metabolic genes and histologic severity, fecal energy loss, and wasting (weight-for-length/height Z < -2.0). Leptin was associated with expression of epithelial carbohydrate metabolism and stem cell renewal genes. Immune response genes were attenuated by giardia colonization. CONCLUSIONS: Children with reduced circulating IGF-1 are more likely to experience stunting. Leptin and a gene signature for lymphocyte activation and dysregulated lipid metabolism are implicated in wasting, suggesting new approaches for EED refractory to nutritional intervention. ClinicalTrials.gov, Number: NCT03588013. (https://clinicaltrials.gov/ct2/show/NCT03588013).

A pilot radiogenomic study of DIPG reveals distinct subgroups with unique clinical trajectories and therapeutic targets.

An adequate understanding of the relationships between radiographic and genomic features in diffuse intrinsic pontine glioma (DIPG) is essential, especially in the absence of universal biopsy, to further characterize the molecular heterogeneity of this disease and determine which patients are most likely to respond to biologically-driven therapies. Here, a radiogenomics analytic approach was applied to a cohort of 28 patients with DIPG. Tumor size and imaging characteristics from all available serial MRIs were evaluated by a neuro-radiologist, and patients were divided into three radiographic response groups (partial response [PR], stable disease [SD], progressive disease [PD]) based on MRI within 2 months of radiotherapy (RT) completion. Whole genome and RNA sequencing were performed on autopsy tumor specimens. We report several key, therapeutically-relevant findings: (1) Certain radiologic features on first and subsequent post-RT MRIs are associated with worse overall survival, including PD following irradiation as well as present, new, and/or increasing peripheral ring enhancement, necrosis, and diffusion restriction. (2) Upregulation of EMT-related genes and distant tumor spread at autopsy are observed in a subset of DIPG patients who exhibit poorer radiographic response to irradiation and/or higher likelihood of harboring H3F3A mutations, suggesting possible benefit of upfront craniospinal irradiation. (3) Additional genetic aberrations were identified, including DYNC1LI1 mutations in a subgroup of patients with PR on post-RT MRI; further investigation into potential roles in DIPG tumorigenesis and/or treatment sensitivity is necessary. (4) Whereas most DIPG tumors have an immunologically "cold" microenvironment, there appears to be a subset which harbor a more inflammatory genomic profile and/or higher mutational burden, with a trend toward improved overall survival and more favorable radiographic response to irradiation, in whom immunotherapy should be considered. This study has begun elucidating relationships between post-RT radiographic response with DIPG molecular profiles, revealing radiogenomically distinct subgroups with unique clinical trajectories and therapeutic targets.

Mucosal Inflammatory and Wound Healing Gene Programs Reveal Targets for Stricturing Behavior in Pediatric Crohn's Disease.

BACKGROUND AND AIMS: Ileal strictures are the major indication for resective surgery in Crohn's disease (CD). We aimed to define ileal gene programs present at diagnosis linked with future stricturing behavior during five year follow-up, and to identify potential small molecules to reverse these gene signatures. METHODS: Antimicrobial serologies and pre-treatment ileal gene expression were assessed in a representative subset of 249 CD patients within the RISK multicenter pediatric CD inception cohort study, including 113 that are unique to this report. These data were used to define genes associated with stricturing behavior and for model testing to predict stricturing behavior. A bioinformatics approach to define small molecules which may reverse the stricturing gene signature was applied. RESULTS: 19 of the 249 patients developed isolated B2 stricturing behavior during follow-up, while 218 remained B1 inflammatory. Using deeper RNA sequencing than in our prior report, we have now defined an inflammatory gene signature including an oncostatin M co-expression signature, tightly associated with extra-cellular matrix (ECM) gene expression in those who developed stricturing complications. We further computationally prioritize small molecules targeting macrophage and fibroblast activation and angiogenesis which may reverse the stricturing gene signature. A model containing ASCA and CBir1 serologies and a refined eight ECM gene set was significantly associated with stricturing development by year five after diagnosis (AUC (95th CI) = 0.82 (0.7-0.94)). CONCLUSION: An ileal gene program for macrophage and fibroblast activation is linked to stricturing complications in treatment naïve pediatric CD, and may inform novel small molecule therapeutic approaches.

AB569, a nontoxic chemical tandem that kills major human pathogenic bacteria.

Antibiotic-resistant superbug bacteria represent a global health problem with no imminent solutions. Here we demonstrate that the combination (termed AB569) of acidified nitrite (A-NO2-) and Na2-EDTA (disodium ethylenediaminetetraacetic acid) inhibited all Gram-negative and Gram-positive bacteria tested. AB569 was also efficacious at killing the model organism Pseudomonas aeruginosa in biofilms and in a murine chronic lung infection model. AB569 was not toxic to human cell lines at bactericidal concentrations using a basic viability assay. RNA-Seq analyses upon treatment of P. aeruginosa with AB569 revealed a catastrophic loss of the ability to support core pathways encompassing DNA, RNA, protein, ATP biosynthesis, and iron metabolism. Electrochemical analyses elucidated that AB569 produced more stable SNO proteins, potentially explaining one mechanism of bacterial killing. Our data implicate that AB569 is a safe and effective means to kill pathogenic bacteria, suggesting that simple strategies could be applied with highly advantageous therapeutic/toxicity index ratios to pathogens associated with a myriad of periepithelial infections and related disease scenarios.

NF1 patient missense variants predict a role for ATM in modifying neurofibroma initiation.

In Neurofibromatosis type 1, NF1 gene mutations in Schwann cells (SC) drive benign plexiform neurofibroma (PNF), and no additional SC changes explain patient-to-patient variability in tumor number. Evidence from twin studies suggests that variable expressivity might be caused by unidentified modifier genes. Whole exome sequencing of SC and fibroblast DNA from the same resected PNFs confirmed biallelic SC NF1 mutations; non-NF1 somatic SC variants were variable and present at low read number. We identified frequent germline variants as possible neurofibroma modifier genes. Genes harboring variants were validated in two additional cohorts of NF1 patients and by variant burden test. Genes including CUBN, CELSR2, COL14A1, ATR and ATM also showed decreased gene expression in some neurofibromas. ATM-relevant DNA repair defects were also present in a subset of neurofibromas with ATM variants, and in some neurofibroma SC. Heterozygous ATM G2023R or homozygous S707P variants reduced ATM protein expression in heterologous cells. In mice, genetic Atm heterozygosity promoted Schwann cell precursor self-renewal and increased tumor formation in vivo, suggesting that ATM variants contribute to neurofibroma initiation. We identify germline variants, rare in the general population, overrepresented in NF1 patients with neurofibromas. ATM and other identified genes are candidate modifiers of PNF pathogenesis.

cellHarmony: cell-level matching and holistic comparison of single-cell transcriptomes.

To understand the molecular pathogenesis of human disease, precision analyses to define alterations within and between disease-associated cell populations are desperately needed. Single-cell genomics represents an ideal platform to enable the identification and comparison of normal and diseased transcriptional cell populations. We created cellHarmony, an integrated solution for the unsupervised analysis, classification, and comparison of cell types from diverse single-cell RNA-Seq datasets. cellHarmony efficiently and accurately matches single-cell transcriptomes using a community-clustering and alignment strategy to compute differences in cell-type specific gene expression over potentially dozens of cell populations. Such transcriptional differences are used to automatically identify distinct and shared gene programs among cell-types and identify impacted pathways and transcriptional regulatory networks to understand the impact of perturbations at a systems level. cellHarmony is implemented as a python package and as an integrated workflow within the software AltAnalyze. We demonstrate that cellHarmony has improved or equivalent performance to alternative label projection methods, is able to identify the likely cellular origins of malignant states, stratify patients into clinical disease subtypes from identified gene programs, resolve discrete disease networks impacting specific cell-types, and illuminate therapeutic mechanisms. Thus, this approach holds tremendous promise in revealing the molecular and cellular origins of complex disease.

Single cell RNA sequencing reveals regional heterogeneity of hepatobiliary innate lymphoid cells in a tissue-enriched fashion.

IL-33 promotes type 2 immunity, epithelial repair, and tissue fibrosis by activating group 2 innate lymphoid cells (ILC2). ILC2 lack all known surface markers of mature T, B, NK, and myeloid cell lineages (Linneg), express the IL-33 receptor ST2, and release type 2 cytokines which contribute to cholangiocyte proliferation and activation of hepatic stellate cells. This pathway results in massive proliferation of the extrahepatic bile duct (EHBD) but also exacerbates liver fibrosis, suggesting that there may be tissue-specific subpopulations of IL-33-induced ILC. To determine the tissue-specific subsets of ILC in the hepatobiliary system, we analyzed CD45+Linneg mononuclear cells from IL-33 treated adult Balb/c mouse liver or EHBD by single cell RNA sequencing. Principal component analysis identified 6 major CD45+Linneg cell classes, two of which were restricted to the EHBD. One of these classes, biliary immature myeloid (BIM) cells, was predicted to interact with ILC2 by a network of shared receptor-ligand pairs. BIM highly expressed Gp49 and ST2 receptors on the cell surface while lacking surface expression of markers for mature myeloid cells. In conclusion, single cell RNA sequencing identified IL-33 responsive cell groups regionally confined to the liver or extrahepatic bile duct, including a novel population of CD45+Linneg Gp49-expressing mononuclear cells.

Ulcerative colitis mucosal transcriptomes reveal mitochondriopathy and personalized mechanisms underlying disease severity and treatment response.

Molecular mechanisms driving disease course and response to therapy in ulcerative colitis (UC) are not well understood. Here, we use RNAseq to define pre-treatment rectal gene expression, and fecal microbiota profiles, in 206 pediatric UC patients receiving standardised therapy. We validate our key findings in adult and paediatric UC cohorts of 408 participants. We observe a marked suppression of mitochondrial genes and function across cohorts in active UC, and that increasing disease severity is notable for enrichment of adenoma/adenocarcinoma and innate immune genes. A subset of severity genes improves prediction of corticosteroid-induced remission in the discovery cohort; this gene signature is also associated with response to anti-TNFα and anti-α4ß7 integrin in adults. The severity and therapeutic response gene signatures were in turn associated with shifts in microbes previously implicated in mucosal homeostasis. Our data provide insights into UC pathogenesis, and may prioritise future therapies for nonresponders to current approaches.

Age-of-diagnosis dependent ileal immune intensification and reduced alpha-defensin in older versus younger pediatric Crohn Disease patients despite already established dysbiosis.

Age-of-diagnosis associated variation in disease location and antimicrobial sero-reactivity has suggested fundamental differences in pediatric Crohn Disease (CD) pathogenesis. This variation may be related to pubertal peak incidence of ileal involvement and Peyer's patches maturation, represented by IFNγ-expressing Th1 cells. However, direct mucosal evidence is lacking. We characterize the global pattern of ileal gene expression and microbial communities in 525 treatment-naive pediatric CD patients and controls (Ctl), stratifying samples by their age-of-diagnosis. We show a robust ileal gene signature notable for higher expression of specific immune genes including GM-CSF and INFγ, and reduced expression of antimicrobial Paneth cell α-defensins, in older compared to younger patients. Reduced α-defensin expression in older patients was associated with higher IFNγ expression. By comparison, the CD-associated ileal dysbiosis, characterized by expansion of Enterobacteriaceae and contraction of Lachnospiraceae and Ruminococcaceae, was already established within the younger group and did not vary systematically with increasing age-of-diagnosis. Multivariate analysis considering individual taxa, however did demonstrate negative associations between Lachnospiraceae and IFNγ, and positive associations between Bacteroides and α-defensin expression. These data provide evidence for maturation of mucosal Th1 immune responses and loss of epithelial antimicrobial α-defensins which are associated with specific taxa with increasing age-of-diagnosis in pediatric CD.

Whole-exome sequencing uncovers oxidoreductases DHTKD1 and OGDHL as linkers between mitochondrial dysfunction and eosinophilic esophagitis.

Eosinophilic esophagitis (EoE) is an allergic inflammatory esophageal disorder with a complex underlying genetic etiology often associated with other comorbidities. Using whole-exome sequencing (WES) of 63 patients with EoE and 60 unaffected family members and family-based trio analysis, we sought to uncover rare coding variants. WES analysis identified 5 rare, damaging variants in dehydrogenase E1 and transketolase domain-containing 1 (DHTKD1). Rare variant burden analysis revealed an overabundance of putative, potentially damaging DHTKD1 mutations in EoE (P = 0.01). Interestingly, we also identified 7 variants in the DHTKD1 homolog oxoglutarate dehydrogenase-like (OGDHL). Using shRNA-transduced esophageal epithelial cells and/or patient fibroblasts, we further showed that disruption of normal DHTKD1 or OGDHL expression blunts mitochondrial function. Finally, we demonstrated that the loss of DHTKD1 expression increased ROS production and induced the expression of viperin, a gene previously shown to be involved in production of Th2 cytokines in T cells. Viperin had increased expression in esophageal biopsies of EoE patients compared with control individuals and was upregulated by IL-13 in esophageal epithelial cells. These data identify a series of rare genetic variants implicating DHTKD1 and OGDHL in the genetic etiology of EoE and underscore a potential pathogenic role for mitochondrial dysfunction in EoE.

Long ncRNA Landscape in the Ileum of Treatment-Naive Early-Onset Crohn Disease.

Background: Long noncoding RNAs (lncRNA) are key regulators of gene transcription and many show tissue-specific expression. We previously defined a novel inflammatory and metabolic ileal gene signature in treatment-naive pediatric Crohn disease (CD). We now extend our analyses to include potential regulatory lncRNA. Methods: Using RNAseq, we systematically profiled lncRNAs and protein-coding gene expression in 177 ileal biopsies. Co-expression analysis was used to identify functions and tissue-specific expression. RNA in situ hybridization was used to validate expression. Real-time polymerase chain reaction was used to test lncRNA regulation by IL-1ß in Caco-2 enterocytes. Results: We characterize widespread dysregulation of 459 lncRNAs in the ileum of CD patients. Using only the lncRNA in discovery and independent validation cohorts showed patient classification as accurate as the protein-coding genes, linking lncRNA to CD pathogenesis. Co-expression and functional annotation enrichment analyses across several tissues and cell types 1showed that the upregulated LINC01272 is associated with a myeloid pro-inflammatory signature, whereas the downregulated HNF4A-AS1 exhibits association with an epithelial metabolic signature. We confirmed tissue-specific expression in biopsies using in situ hybridization, and validated regulation of prioritized lncRNA upon IL-1ß exposure in differentiated Caco-2 cells. Finally, we identified significant correlations between LINC01272 and HNF4A-AS1 expression and more severe mucosal injury. Conclusions: We systematically define differentially expressed lncRNA in the ileum of newly diagnosed pediatric CD. We show lncRNA utility to correctly classify disease or healthy states and demonstrate their regulation in response to an inflammatory signal. These lncRNAs, after mechanistic exploration, may serve as potential new tissue-specific targets for RNA-based interventions.