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.

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

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.

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.

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.

Molecular, phenotypic, and sample-associated data to describe pluripotent stem cell lines and derivatives.

The use of induced pluripotent stem cells (iPSC) derived from independent patients and sources holds considerable promise to improve the understanding of development and disease. However, optimized use of iPSC depends on our ability to develop methods to efficiently qualify cell lines and protocols, monitor genetic stability, and evaluate self-renewal and differentiation potential. To accomplish these goals, 57 stem cell lines from 10 laboratories were differentiated to 7 different states, resulting in 248 analyzed samples. Cell lines were differentiated and characterized at a central laboratory using standardized cell culture methodologies, protocols, and metadata descriptors. Stem cell and derived differentiated lines were characterized using RNA-seq, miRNA-seq, copy number arrays, DNA methylation arrays, flow cytometry, and molecular histology. All materials, including raw data, metadata, analysis and processing code, and methodological and provenance documentation are publicly available for re-use and interactive exploration at https://www.synapse.org/pcbc. The goal is to provide data that can improve our ability to robustly and reproducibly use human pluripotent stem cells to understand development and disease.

Integrated Genomic Analysis of Diverse Induced Pluripotent Stem Cells from the Progenitor Cell Biology Consortium.

The rigorous characterization of distinct induced pluripotent stem cells (iPSC) derived from multiple reprogramming technologies, somatic sources, and donors is required to understand potential sources of variability and downstream potential. To achieve this goal, the Progenitor Cell Biology Consortium performed comprehensive experimental and genomic analyses of 58 iPSC from ten laboratories generated using a variety of reprogramming genes, vectors, and cells. Associated global molecular characterization studies identified functionally informative correlations in gene expression, DNA methylation, and/or copy-number variation among key developmental and oncogenic regulators as a result of donor, sex, line stability, reprogramming technology, and cell of origin. Furthermore, X-chromosome inactivation in PSC produced highly correlated differences in teratoma-lineage staining and regulator expression upon differentiation. All experimental results, and raw, processed, and metadata from these analyses, including powerful tools, are interactively accessible from a new online portal at https://www.synapse.org to serve as a reusable resource for the stem cell community.

IRAK1 is a novel DEK transcriptional target and is essential for head and neck cancer cell survival.

The chromatin-binding DEK protein was recently reported to promote the growth of HPV+ and HPV- head and neck squamous cell carcinomas (HNSCCs). Relevant cellular and molecular mechanism(s) controlled by DEK in HNSCC remain poorly understood. While DEK is known to regulate specific transcriptional targets, global DEK-dependent gene networks in HNSCC are unknown. To identify DEK transcriptional signatures we performed RNA-Sequencing (RNA-Seq) in HNSCC cell lines that were either proficient or deficient for DEK. Bioinformatic analyses and subsequent validation revealed that IRAK1, a regulator of inflammatory signaling, and IRAK1-dependent regulatory networks were significantly repressed upon DEK knockdown in HNSCC. According to TCGA data, 14% of HNSCC specimens overexpressed IRAK1, thus supporting possible oncogenic functions. Furthermore, genetic or pharmacologic inhibition of IRAK1 in HNSCC cell lines was sufficient to attenuate downstream signaling such as ERK1/2 and to induce HNSCC cell death by apoptosis. Finally, targeting DEK and IRAK1 simultaneously enhanced cell death as compared to targeting either alone. Our findings reveal that IRAK1 promotes cell survival and is an attractive therapeutic target in HNSCC cells. Thus, we propose a model wherein IRAK1 stimulates tumor signaling and phenotypes both independently and in conjunction with DEK.

Intrinsic Age-Dependent Changes and Cell-Cell Contacts Regulate Nephron Progenitor Lifespan.

During fetal development, nephrons of the metanephric kidney form from a mesenchymal progenitor population that differentiates en masse before or shortly after birth. We explored intrinsic and extrinsic mechanisms controlling progenitor lifespan in a transplantation assay that allowed us to compare engraftment of old and young progenitors into the same young niche. The progenitors displayed an age-dependent decrease in proliferation and concomitant increase in niche exit rates. Single-cell transcriptome profiling revealed progressive age-dependent changes, with heterogeneity increasing in older populations. Age-dependent elevation in mTor and reduction in Fgf20 could contribute to increased exit rates. Importantly, 30% of old progenitors remained in the niche for up to 1 week post engraftment, a net gain of 50% to their lifespan, but only if surrounded by young neighbors. We provide evidence in support of a model in which intrinsic age-dependent changes affect inter-progenitor interactions that drive cessation of nephrogenesis.