A big data approach to evaluate receipt of optimal care in childhood cerebral palsy

PURPOSE

Through automated electronic health record (EHR) data extraction and analysis, this project systematically quantified actual care delivery for children with cerebral palsy (CP) and evaluated alignment with current evidence-based recommendations.

METHODS

Utilizing EHR data for over 8000 children with CP, we developed an approach to define and quantify receipt of optimal care, and pursued proof-of-concept with two children with unilateral CP, Gross Motor Function Classification System (GMFCS) Level II. Optimal care was codified as a cluster of four components including physical medicine and rehabilitation (PMR) care, spasticity management, physical therapy (PT), and occupational therapy (OT). A Receipt of Care Score (ROCS) quantified the degree of adherence to recommendations and was compared with the Pediatric Outcomes Data Collection Instrument (PODCI) and Pediatric Quality of Life Inventory (PEDS QL).

RESULTS

The two children (12 year old female, 13 year old male) had nearly identical PMR and spasticity component scores while PT and OT scores were more divergent. Functional outcomes were higher for the child who had higher adjusted ROCS.

CONCLUSIONS

ROCSs demonstrate variation in real-world care delivered over time and differentiate between components of care. ROCSs reflect overall function and quality of life. The ROCS methods developed are novel, robust, and scalable and will be tested in a larger sample.IMPLICATIONS FOR REHABILITATIONOptimal practice, with an emphasis on integrated multidisciplinary care, can be defined and quantified utilizing evidence-based recommendations.Receipt of optimal care for childhood cerebral palsy can be scored using existing electronic health record data.Big Data approaches can contribute to the understanding of current care and inform approaches for improved care.

Schizophrenia-associated NRXN1 deletions induce developmental-timing- and cell-type-specific vulnerabilities in human brain organoids

De novo mutations and copy number deletions in NRXN1 (2p16.3) pose a significant risk for schizophrenia (SCZ). It is unclear how NRXN1 deletions impact cortical development in a cell type-specific manner and disease background modulates these phenotypes. Here, we leveraged human pluripotent stem cell-derived forebrain organoid models carrying NRXN1 heterozygous deletions in isogenic and SCZ patient genetic backgrounds and conducted single-cell transcriptomic analysis over the course of brain organoid development from 3 weeks to 3.5 months. Intriguingly, while both deletions similarly impacted molecular pathways associated with ubiquitin-proteasome system, alternative splicing, and synaptic signaling in maturing glutamatergic and GABAergic neurons, SCZ-NRXN1 deletions specifically perturbed developmental trajectories of early neural progenitors and accumulated disease-specific transcriptomic signatures. Using calcium imaging, we found that both deletions led to long-lasting changes in spontaneous and synchronous neuronal networks, implicating synaptic dysfunction. Our study reveals developmental-timing- and cell-type-dependent actions of NRXN1 deletions in unique genetic contexts.

Spatial transcriptome profiling uncovers metabolic regulation of left-right patterning

Left-right patterning disturbance can cause severe birth defects, but it remains least understood of the three body axes. We uncovered an unexpected role for metabolic regulation in left-right patterning. Analysis of the first spatial transcriptome profile of left-right patterning revealed global activation of glycolysis, accompanied by right-sided expression of Bmp7 and genes regulating insulin growth factor signaling. Cardiomyocyte differentiation was left-biased, which may underlie the specification of heart looping orientation. This is consistent with known Bmp7 stimulation of glycolysis and glycolysis suppression of cardiomyocyte differentiation. Liver/lung laterality may be specified via similar metabolic regulation of endoderm differentiation. Myo1d , found to be left-sided, was shown to regulate gut looping in mice, zebrafish, and human. Together these findings indicate metabolic regulation of left-right patterning. This could underlie high incidence of heterotaxy-related birth defects in maternal diabetes, and the association of PFKP, allosteric enzyme regulating glycolysis, with heterotaxy. This transcriptome dataset will be invaluable for interrogating birth defects involving laterality disturbance.

Outpatient hospital utilization after single event multi-level surgery in children with cerebral palsy

PURPOSE

This study aimed to examine outpatient hospital utilization (number of specialties seen and number of visits to each specialty) in the year after single event multi-level surgery (SEMLS) in children with cerebral palsy (CP), and to determine if utilization differs across the medical center in the year after compared to the year before SEMLS.

METHODS

This retrospective cross-sectional study used electronic medical record data of outpatient hospital utilization in children with CP who underwent SEMLS.

RESULTS

Thirty children with CP (Gross Motor Function Classification System Levels I-V, mean age of 9.9 years) were included. In the year after surgery, a significant difference (p = 0.001) was found for the number of specialties seen, with non-ambulatory children seeing more specialties than ambulatory children. No statistically significant difference was found between the number of outpatient visits to each specialty in the year after SEMLS. Compared to the year before SEMLS, fewer therapy visits occurred in the year after SEMLS (p <  0.001) but significantly more visits to orthopaedics (p = 0.001) and radiology (p = 0.001).

CONCLUSION

Children with CP had fewer therapy visits but more orthopaedic and radiology visits the year after SEMLS. Nearly half of the children were non-ambulatory. Examination of care needs in children with CP undergoing SEMLS is justified with consideration of ambulatory status, surgical burden, and post-operative immobilization.

Gross motor function prediction using natural language processing in cerebral palsy

AIM

To predict ambulatory status and Gross Motor Function Classification System (GMFCS) levels in patients with cerebral palsy (CP) by applying natural language processing (NLP) to electronic health record (EHR) clinical notes.

METHOD

Individuals aged 8 to 26 years with a diagnosis of CP in the EHR between January 2009 and November 2020 (~12 years of data) were included in a cross-sectional retrospective cohort of 2483 patients. The cohort was divided into train-test and validation groups. Positive predictive value, sensitivity, specificity, and area under the receiver operating curve (AUC) were calculated for prediction of ambulatory status and GMFCS levels.

RESULTS

The median age was 15 years (interquartile range 10-20 years) for the total cohort, with 56% being male and 75% White. The validation group resulted in 70% sensitivity, 88% specificity, 81% positive predictive value, and 0.89 AUC for predicting ambulatory status. NLP applied to the EHR differentiated between GMFCS levels I-II and III (15% sensitivity, 96% specificity, 46% positive predictive value, and 0.71 AUC); and IV and V (81% sensitivity, 51% specificity, 70% positive predictive value, and 0.75 AUC).

INTERPRETATION

NLP applied to the EHR demonstrated excellent differentiation between ambulatory and non-ambulatory status, and good differentiation between GMFCS levels I-II and III, and IV and V. Clinical use of NLP may help to individualize functional characterization and management.

WHAT THIS PAPER ADDS

Natural language processing (NLP) applied to the electronic health record (EHR) can predict ambulatory status in children with cerebral palsy (CP). NLP provides good prediction of Gross Motor Function Classification System level in children with CP using the EHR. NLP methods described could be integrated in an EHR system to provide real-time information.

CASK loss of function differentially regulates neuronal maturation and synaptic function in human induced cortical excitatory neurons

Loss-of-function (LOF) mutations in CASK cause severe developmental phenotypes, including microcephaly with pontine and cerebellar hypoplasia, X-linked intellectual disability, and autism. Unraveling the pathological mechanisms of CASK-related disorders has been challenging owing to limited human cellular models to study the dynamic roles of this molecule during neuronal maturation and synapse development. Here, we investigate cell-autonomous functions of CASK in cortical excitatory induced neurons (iNs) generated from CASK knockout (KO) isogenic human embryonic stem cells (hESCs) using gene expression, morphometrics, and electrophysiology. While immature CASK KO iNs show robust neuronal outgrowth, mature CASK KO iNs display severe defects in synaptic transmission and synchronized network activity without compromising neuronal morphology and synapse numbers. In the developing human cortical excitatory neurons, CASK functions to promote both structural integrity and establishment of cortical excitatory neuronal networks. These results lay the foundation for future studies identifying suppressors of such phenotypes relevant to human patients.

Uncompensated mitochondrial oxidative stress underlies heart failure in an iPSC-derived model of congenital heart disease

Hypoplastic left heart syndrome (HLHS) is a severe congenital heart disease with 30% mortality from heart failure (HF) in the first year of life, but the cause of early HF remains unknown. Induced pluripotent stem-cell-derived cardiomyocytes (iPSC-CM) from patients with HLHS showed that early HF is associated with increased apoptosis, mitochondrial respiration defects, and redox stress from abnormal mitochondrial permeability transition pore (mPTP) opening and failed antioxidant response. In contrast, iPSC-CM from patients without early HF showed normal respiration with elevated antioxidant response. Single-cell transcriptomics confirmed that early HF is associated with mitochondrial dysfunction accompanied with endoplasmic reticulum (ER) stress. These findings indicate that uncompensated oxidative stress underlies early HF in HLHS. Importantly, mitochondrial respiration defects, oxidative stress, and apoptosis were rescued by treatment with sildenafil to inhibit mPTP opening or TUDCA to suppress ER stress. Together these findings point to the potential use of patient iPSC-CM for modeling clinical heart failure and the development of therapeutics.

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.

Response to comment on 'SARS-CoV-2 suppresses anticoagulant and fibrinolytic gene expression in the lung'

Early in the SARS-CoV-2 pandemic, we compared transcriptome data from hospitalized COVID-19 patients and control patients without COVID-19. We found changes in procoagulant and fibrinolytic gene expression in the lungs of COVID-19 patients (Mast et al., 2021). These findings have been challenged based on issues with the samples (Fitzgerald and Jamieson, 2022). We have revisited our previous analyses in the light of this challenge and find that these new analyses support our original conclusions.

Betacoronavirus-specific alternate splicing

Viruses can subvert a number of cellular processes including splicing in order to block innate antiviral responses, and many viruses interact with cellular splicing machinery. SARS-CoV-2 infection was shown to suppress global mRNA splicing, and at least 10 SARS-CoV-2 proteins bind specifically to one or more human RNAs. Here, we investigate 17 published experimental and clinical datasets related to SARS-CoV-2 infection, datasets from the betacoronaviruses SARS-CoV and MERS, as well as Streptococcus pneumonia, HCV, Zika virus, Dengue virus, influenza H3N2, and RSV. We show that genes showing differential alternative splicing in SARS-CoV-2 have a similar functional profile to those of SARS-CoV and MERS and affect a diverse set of genes and biological functions, including many closely related to virus biology. Additionally, the differentially spliced transcripts of cells infected by coronaviruses were more likely to undergo intron-retention, contain a pseudouridine modification, and have a smaller number of exons as compared with differentially spliced transcripts in the control groups. Viral load in clinical COVID-19 samples was correlated with isoform distribution of differentially spliced genes. A significantly higher number of ribosomal genes are affected by differential alternative splicing and gene expression in betacoronavirus samples, and the betacoronavirus differentially spliced genes are depleted for binding sites of RNA-binding proteins. Our results demonstrate characteristic patterns of differential splicing in cells infected by SARS-CoV-2, SARS-CoV, and MERS. The alternative splicing changes observed in betacoronaviruses infection potentially modify a broad range of cellular functions, via changes in the functions of the products of a diverse set of genes involved in different biological processes.

Identification of distinct tumor cell populations and key genetic mechanisms through single cell sequencing in hepatoblastoma

Hepatoblastoma (HB) is the most common primary liver malignancy of childhood, and molecular investigations are limited and effective treatment options for chemoresistant disease are lacking. There is a knowledge gap in the investigation of key driver cells of HB in tumor. Here we show single cell ribonucleic acid sequencing (scRNAseq) analysis of human tumor, background liver, and patient derived xenograft (PDX) to demonstrate gene expression patterns within tumor and to identify intratumor cell subtype heterogeneity to define differing roles in pathogenesis based on intracellular signaling in pediatric HB. We have identified a driver tumor cell cluster in HB by genetic expression which can be examined to define disease mechanism and treatments. Identification of both critical mechanistic pathways combined with unique cell populations provide the basis for discovery and investigation of novel treatment strategies in vitro and in vivo.

Betacoronavirus-specific alternate splicing

Viruses can subvert a number of cellular processes in order to block innate antiviral responses, and many viruses interact with cellular splicing machinery. SARS-CoV-2 infection was shown to suppress global mRNA splicing, and at least 10 SARS-CoV-2 proteins bind specifically to one or more human RNAs. Here, we investigate 17 published experimental and clinical datasets related to SARS-CoV-2 infection as well as datasets from the betacoronaviruses SARS-CoV and MERS as well as Streptococcus pneumonia, HCV, Zika virus, Dengue virus, influenza H3N2, and RSV. We show that genes showing differential alternative splicing in SARS-CoV-2 have a similar functional profile to those of SARS-CoV and MERS and affect a diverse set of genes and biological functions, including many closely related to virus biology. Additionally, the differentially spliced transcripts of cells infected by coronaviruses were more likely to undergo intron-retention, contain a pseudouridine modification and a smaller number of exons than differentially spliced transcripts in the control groups. Viral load in clinical COVID-19 samples was correlated with isoform distribution of differentially spliced genes. A significantly higher number of ribosomal genes are affected by DAS and DGE in betacoronavirus samples, and the betacoronavirus differentially spliced genes are depleted for binding sites of RNA-binding proteins. Our results demonstrate characteristic patterns of differential splicing in cells infected by SARS-CoV-2, SARS-CoV, and MERS, potentially modifying a broad range of cellular functions and affecting a diverse set of genes and biological functions.

Abstract 908: Preclinical in vivo characterization of a first-in-class humanized antibody targeting alternatively spliced tissue factor

Tissue Factor, the initiator of the extrinsic pathway of coagulation, and its isoform alternatively spliced Tissue Factor (asTF) are often overexpressed in cancer cells. Previously, we have shown that RabMab1, a rabbit monoclonal antibody specific for human asTF, disrupts the binding of asTF to beta-integrins and thereby inhibits the growth of breast cancer and pancreatic ductal adenocarcinoma (PDAC) cell lines in vitro and in vivo. Here, we report on humanization of RabMab1, assessment of its binding characteristics, and determination of its in vivo properties. The variable regions of the heavy and light chains of hybridoma derived RabMab1 were cloned and used to generate a rabbit/human chimera (cRabMab1). Substitutions of species-specific residues were then carried out to produce humanized RabMab1 (hRabMab1; IgG1 isotype). Antigen binding was assessed via anti-human IgG Fc capture biosensor assays. cRabMab1 was found to have a KD of 4.24 nM whereas hRabMab1 was found to have a KD in the low picomolar range, which could not be precisely measured due to an extremely slow off-rate. The in vivo half-lives of each were determined to be 280 and 908 hours, respectively. Orthotopic co-implantation of cRabMab1 with high-grade, human PDAC cell line Pt45.P1 in athymic nude mice resulted in tumors that were 82% and 92% smaller than tumors in the vehicle and IgG control group groups, respectively, with no statistical difference between vehicle and IgG isotype control groups. Intravenous administration of hRabMab1 slowed the growth of pre-formed orthotopic Pt45.P1 tumors in athymic nude mice by 64% and 61% compared to vehicle and isotype control groups, respectively; again, there was no difference between vehicle and IgG isotype control groups. Immunohistochemical analysis of tumor tissue revealed a statistically significant 68% reduction in neovascularization (CD31), a 58% reduction in M2-polarized macrophages (CD206), and a 24% reduction in proliferating cells (Ki67) in the hRabMab1 cohort. RNAseq analysis of tumor tissue revealed that components of the focal adhesion system were the most affected by hRabMab1 treatment. In addition, a significant suppression of pathways that promote mitosis and the cell cycle was seen in hRabMab1-treated tumors. By qRT-PCR, we validated a greater than 82% reduction in the expression of HBEGF, a gene encoding a cell surface mitogen, and a greater than 68% reduction in the expression of STN1, a gene encoding a component of the CST-complex, a telomere replication and maintenance structure. This is the first proof-of-concept study whereby a novel biologic that inhibits asTF has been used as a systemically administered single agent, with encouraging results. Because hRabMab1 has a favorable pharmacokinetic profile and is able to suppress PDAC tumor cell growth in vivo, it is an attractive candidate for further clinical development.

Citation Format: Clayton S. Lewis, Aniruddha Karve, Kateryna Matiash, Timothy Stone, Jingxing Li, Jordan Wang, Henri Versteeg, Bruce Aronow, Syed Ahmad, Pankaj Desai, Vladimir Bogdanov. Preclinical in vivo characterization of a first-in-class humanized antibody targeting alternatively spliced tissue factor [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 908.

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.

Describing the Delivery of Evidence-Based Physical Therapy Intervention to Individuals With Cerebral Palsy

PURPOSE

To characterize by evidence grades and examine variation in type of physical therapy intervention delivered in routine clinical care in individuals with cerebral palsy (CP).

METHODS

Retrospective data collection from the electronic record over 1 year at a tertiary care pediatric outpatient therapy division.

RESULTS

Four hundred sixty-five individuals with CP received 28 344 interventions during 4335 treatment visits. Sixty-six percent of interventions were evidence-based interventions (EBIs). Significant variation was demonstrated across Gross Motor Function Classification System levels, with children classified as level V receiving the least and level III the most. The most frequent EBIs delivered were caregiver education, motor control, functional strengthening, ankle-foot orthoses, treadmill training, and fit of adaptive equipment.

CONCLUSIONS

Further work is needed to determine whether amount of EBI is related to better outcomes. Combining this information with other aspects of dose (intensity, time, and frequency) may elucidate the contribution of each with outcomes.

SARS-CoV-2 suppresses anticoagulant and fibrinolytic gene expression in the lung

Extensive fibrin deposition in the lungs and altered levels of circulating blood coagulation proteins in COVID-19 patients imply local derangement of pathways that limit fibrin formation and/or promote its clearance. We examined transcriptional profiles of bronchoalveolar lavage fluid (BALF) samples to identify molecular mechanisms underlying these coagulopathies. mRNA levels for regulators of the kallikrein-kinin (C1-inhibitor), coagulation (thrombomodulin, endothelial protein C receptor), and fibrinolytic (urokinase and urokinase receptor) pathways were significantly reduced in COVID-19 patients. While transcripts for several coagulation proteins were increased, those encoding tissue factor, the protein that initiates coagulation and whose expression is frequently increased in inflammatory disorders, were not increased in BALF from COVID-19 patients. Our analysis implicates enhanced propagation of coagulation and decreased fibrinolysis as drivers of the coagulopathy in the lungs of COVID-19 patients.

Ontology-guided segmentation and object identification for developmental mouse lung immunofluorescent images

BACKGROUND

Immunofluorescent confocal microscopy uses labeled antibodies as probes against specific macromolecules to discriminate between multiple cell types. For images of the developmental mouse lung, these cells are themselves organized into densely packed higher-level anatomical structures. These types of images can be challenging to segment automatically for several reasons, including the relevance of biomedical context, dependence on the specific set of probes used, prohibitive cost of generating labeled training data, as well as the complexity and dense packing of anatomical structures in the image. The use of an application ontology helps surmount these challenges by combining image data with its metadata to provide a meaningful biological context, modeled after how a human expert would make use of contextual information to identify histological structures, that constrains and simplifies the process of segmentation and object identification.

RESULTS

We propose an innovative approach for the semi-supervised analysis of complex and densely packed anatomical structures from immunofluorescent images that utilizes an application ontology to provide a simplified context for image segmentation and object identification. We describe how the logical organization of biological facts in the form of an ontology can provide useful constraints that facilitate automatic processing of complex images. We demonstrate the results of ontology-guided segmentation and object identification in mouse developmental lung images from the Bioinformatics REsource ATlas for the Healthy lung database of the Molecular Atlas of Lung Development (LungMAP1) program CONCLUSION: We describe a novel ontology-guided approach to segmentation and classification of complex immunofluorescence images of the developing mouse lung. The ontology is used to automatically generate constraints for each image based on its biomedical context, which facilitates image segmentation and classification.

The elephant in the lung: Integrating lineage-tracing, molecular markers, and single cell sequencing data to identify distinct fibroblast populations during lung development and regeneration

During lung development, the mesenchyme and epithelium are dependent on each other for instructive morphogenic cues that direct proliferation, cellular differentiation and organogenesis. Specification of epithelial and mesenchymal cell lineages occurs in parallel, forming cellular subtypes that guide the formation of both transitional developmental structures and the permanent architecture of the adult lung. While epithelial cell types and lineages have been relatively well-defined in recent years, the definition of mesenchymal cell types and lineage relationships has been more challenging. Transgenic mouse lines with permanent and inducible lineage tracers have been instrumental in identifying lineage relationships among epithelial progenitor cells and their differentiation into distinct airway and alveolar epithelial cells. Lineage tracing experiments with reporter mice used to identify fibroblast progenitors and their lineage trajectories have been limited by the number of cell specific genes and the developmental timepoint when the lineage trace was activated. In this review, we discuss major developmental mesenchymal lineages, focusing on time of origin, major cell type, and other lineage derivatives, as well as the transgenic tools used to find and define them. We describe lung fibroblasts using function, location, and molecular markers in order to compare and contrast cells with similar functions. The temporal and cell-type specific expression of fourteen "fibroblast lineage" genes were identified in single-cell RNA-sequencing data from LungMAP in the LGEA database. Using these lineage signature genes as guides, we clustered murine lung fibroblast populations from embryonic day 16.5 to postnatal day 28 (E16.5-PN28) and generated heatmaps to illustrate expression of transcription factors, signaling receptors and ligands in a temporal and population specific manner.

A mechanistic model and therapeutic interventions for COVID-19 involving a RAS-mediated bradykinin storm

Neither the disease mechanism nor treatments for COVID-19 are currently known. Here, we present a novel molecular mechanism for COVID-19 that provides therapeutic intervention points that can be addressed with existing FDA-approved pharmaceuticals. The entry point for the virus is ACE2, which is a component of the counteracting hypotensive axis of RAS. Bradykinin is a potent part of the vasopressor system that induces hypotension and vasodilation and is degraded by ACE and enhanced by the angiotensin1-9 produced by ACE2. Here, we perform a new analysis on gene expression data from cells in bronchoalveolar lavage fluid (BALF) from COVID-19 patients that were used to sequence the virus. Comparison with BALF from controls identifies a critical imbalance in RAS represented by decreased expression of ACE in combination with increases in ACE2, renin, angiotensin, key RAS receptors, kinogen and many kallikrein enzymes that activate it, and both bradykinin receptors. This very atypical pattern of the RAS is predicted to elevate bradykinin levels in multiple tissues and systems that will likely cause increases in vascular dilation, vascular permeability and hypotension. These bradykinin-driven outcomes explain many of the symptoms being observed in COVID-19.

TCGA Kidney Cell Atlas: A differential gene expression database for the dissection of tumor-specific gene modules

e17078

Background: The Cancer Genome Atlas project has become a leading source for data that has allowed the identification of a broad range of human cancer tumor types and subtypes and has revealed deep complexity with respect to the differentiation, or lack thereof, among human cancers. In particular, differential gene expression analyses have revealed a wealth of active oncogenic pathways, underlying gene mutation drivers, discriminative markers, and candidate therapeutic targets. Despite its rich composition, several factors have led to it not attaining the utility it would seem to offer. Methods: To study this, we dissected molecular subtypes in the TCGA and used the Pan-Kidney (n = 1022 samples) Portion within it to determine where obstacles seem to limit its utility. We re-clustered the renal carcinomas to create more appropriate histology annotations for these samples. The molecular subtypes were then found through K-means using differentially expressed known developmental regulators per histological annotation. After deriving these new annotations, the histology and molecular subtypes were compared to one another via T-test to generate gene modules that characterize these classes/subclasses. Results: We identified a number of factors that include inconsistent metadata attributes, apparent misclassification of histological subtypes, and molecular subtypes that do not match with that obtained by focused approaches to rederive principle subclasses. Our gene modules showed a molecular subtype of clear cell renal carcinoma that was enriched for vascular development and nephron development. In general, the clear cell renal carcinoma and papillary renal cell carcinoma cohorts both showed significant co-expression with atlases that were enriched for genes involved in kidney development. Conclusions: Our atlas highlights the limitations of the current TCGA atlas and provides another tool to capture the rich insights from the TCGA repository through the efforts explained above, highlighted by its use in kidney carcinoma.

Proteasomal adaptations underlying carfilzomib-resistance in human bone marrow plasma cells

Donor-specific antibodies (DSAs) have a deleterious effect on allografts and remain a major immunologic barrier in transplantation. Current therapies to eliminate DSAs are ineffective in highly HLA-sensitized patients. Proteasome inhibitors have been employed as a strategy to target bone marrow plasma cells (BMPCs), the source of long-term antibody production; however, their efficacy has been limited by poorly defined drug-resistance mechanisms. Here, we performed transcriptomic profiling of CD138⁺ BMPCs that survived in vivo desensitization therapy with the proteasome inhibitor carfilzomib to identify mechanisms of drug resistance. The results revealed a genomic signature that included increased expression of the immunoproteasome, a highly specialized proteasomal variant. Western blotting and functional studies demonstrated that catalytically active immunoproteasomes and the immunoproteasome activator PA28 were upregulated in carfilzomib-resistant BMPCs. Carfilzomib-resistant BMPCs displayed reduced sensitivity to the proteasome inhibitors carfilzomib, bortezomib, and ixazomib, but enhanced sensitivity to an immunoproteasome-specific inhibitor ONX-0914. Finally, in vitro carfilzomib treatment of BMPCs from HLA-sensitized patients increased levels of the immunoproteasome β5i (PSMB8) catalytic subunit suggesting that carfilzomib therapy directly induces an adaptive immunoproteasome response. Taken together, our results indicate that carfilzomib induces structural changes in proteasomes and immunoproteasome formation.

The Molecular Signature of Megakaryocyte-Erythroid Progenitors Reveals a Role for the Cell Cycle in Fate Specification

Megakaryocytic-erythroid progenitors (MEPs) give rise to the cells that produce red blood cells and platelets. Although the mechanisms underlying megakaryocytic (MK) and erythroid (E) maturation have been described, those controlling their specification from MEPs are unknown. Single-cell RNA sequencing of primary human MEPs, common myeloid progenitors (CMPs), megakaryocyte progenitors, and E progenitors revealed a distinct transitional MEP signature. Inferred regulatory transcription factors (TFs) were associated with differential expression of cell cycle regulators. Genetic manipulation of selected TFs validated their role in lineage specification and demonstrated coincident modulation of the cell cycle. Genetic and pharmacologic modulation demonstrated that cell cycle activation is sufficient to promote E versus MK specification. These findings, obtained from healthy human cells, lay a foundation to study the mechanisms underlying benign and malignant disease states of the megakaryocytic and E lineages.

Bipotent megakaryocytic-erythroid progenitors (MEPs) produce megakaryocytic and erythroid cells. Using single-cell RNA sequencing of primary human MEPs and their upstream and downstream progenitors, Lu et al. show that MEPs are a unique transitional population. Functional and molecular studies show that MEP lineage fate is toggled by cell cycle speed.

Discovery of SERPINA3 as a candidate urinary biomarker of lupus nephritis activity

Objectives

We used an unbiased proteomics approach to identify candidate urine biomarkers (CUBMs) predictive of LN chronicity and pursued their validation in a larger cohort.

Methods

In this cross-sectional pilot study, we selected urine collected at kidney biopsy from 20 children with varying levels of LN damage (discovery cohort) and performed proteomic analysis using isobaric tags for relative and absolute quantification (iTRAQ). We identified differentially excreted proteins based on degree of LN chronicity and sought to distinguish markers exhibiting different relative expression patterns using hierarchically clustered log10-normalized relative abundance data with linked and distinct functions by biological network analyses. For each CUBM, we performed specific ELISAs on urine from a validation cohort (n = 41) and analysis of variance to detect differences between LN chronicity, with LN activity adjustment. We evaluated for CUBM expression in LN biopsies with immunohistochemistry.

Results

iTRAQ detected 112 proteins in urine from the discovery cohort, 51 quantifiable in all replicates. Simple analysis of variance revealed four differentially expressed, chronicity-correlated proteins (P-values < 0.05). Further correlation and network analyses led to selection of seven CUBMs for LN chronicity. In the validation cohort, none of the CUBMs distinguished LN chronicity degree; however, urine SERPINA3 demonstrated a moderate positive correlation with LN histological activity. Immunohistochemistry further demonstrated SERPINA3 staining in proximal tubular epithelial and endothelial cells.

Conclusion

We identified SERPINA3, a known inhibitor of neutrophil cathepsin G and angiotensin II production, as a potential urine biomarker to help quantify LN activity.

Synergistic effects of common schizophrenia risk variants

The mechanisms by which common risk variants of small effect interact to contribute to complex genetic disorders are unclear. Here, we apply a genetic approach, using isogenic human induced pluripotent stem cells, to evaluate the effects of schizophrenia (SZ)-associated common variants predicted to function as SZ expression quantitative trait loci (eQTLs). By integrating CRISPR-mediated gene editing, activation and repression technologies to study one putative SZ eQTL (FURIN rs4702) and four top-ranked SZ eQTL genes (FURIN, SNAP91, TSNARE1 and CLCN3), our platform resolves pre- and postsynaptic neuronal deficits, recapitulates genotype-dependent gene expression differences and identifies convergence downstream of SZ eQTL gene perturbations. Our observations highlight the cell-type-specific effects of common variants and demonstrate a synergistic effect between SZ eQTL genes that converges on synaptic function. We propose that the links between rare and common variants implicated in psychiatric disease risk constitute a potentially generalizable phenomenon occurring more widely in complex genetic disorders.

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.

Inhibition of STEP61 ameliorates deficits in mouse and hiPSC-based schizophrenia models

The brain-specific tyrosine phosphatase, STEP (STriatal-Enriched protein tyrosine Phosphatase) is an important regulator of synaptic function. STEP normally opposes synaptic strengthening by increasing N-methyl D-aspartate glutamate receptor (NMDAR) internalization through dephosphorylation of GluN2B and inactivation of the kinases extracellular signal-regulated kinase 1/2 and Fyn. Here we show that STEP₆₁ is elevated in the cortex in the Nrg1^+/- knockout mouse model of schizophrenia (SZ). Genetic reduction or pharmacological inhibition of STEP prevents the loss of NMDARs from synaptic membranes and reverses behavioral deficits in Nrg1^+/- mice. STEP₆₁ protein is also increased in cortical lysates from the central nervous system-specific ErbB2/4 mouse model of SZ, as well as in human induced pluripotent stem cell (hiPSC)-derived forebrain neurons and Ngn2-induced excitatory neurons, from two independent SZ patient cohorts. In these selected SZ models, increased STEP₆₁ protein levels likely reflect reduced ubiquitination and degradation. These convergent findings from mouse and hiPSC SZ models provide evidence for STEP₆₁ dysfunction in SZ.

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.

Obesity alters the long-term fitness of the hematopoietic stem cell compartment through modulation of Gfi1 expression

Lee et al. show that established obesity alters the composition and long-term fitness of the hematopoietic stem cell (HSC) compartment, in part through a Gfi1-dependent HSC regulatory program that is activated by the chronic oxidative stress associated with this condition.

Obesity is a chronic organismal stress that disrupts multiple systemic and tissue-specific functions. In this study, we describe the impact of obesity on the activity of the hematopoietic stem cell (HSC) compartment. We show that obesity alters the composition of the HSC compartment and its activity in response to hematopoietic stress. The impact of obesity on HSC function is progressively acquired but persists after weight loss or transplantation into a normal environment. Mechanistically, we establish that the oxidative stress induced by obesity dysregulates the expression of the transcription factor Gfi1 and that increased Gfi1 expression is required for the abnormal HSC function induced by obesity. These results demonstrate that obesity produces durable changes in HSC function and phenotype and that elevation of Gfi1 expression in response to the oxidative environment is a key driver of the altered HSC properties observed in obesity. Altogether, these data provide phenotypic and mechanistic insight into durable hematopoietic dysregulations resulting from obesity.

A Novel Pkhd1 Mutation Interacts with the Nonobese Diabetic Genetic Background To Cause Autoimmune Cholangitis

We previously reported that NOD.c3c4 mice develop spontaneous autoimmune biliary disease (ABD) with anti-mitochondrial Abs, histopathological lesions, and autoimmune T lymphocytes similar to human primary biliary cholangitis. In this article, we demonstrate that ABD in NOD.c3c4 and related NOD ABD strains is caused by a chromosome 1 region that includes a novel mutation in polycystic kidney and hepatic disease 1 (Pkhd1). We show that a long terminal repeat element inserted into intron 35 exposes an alternative polyadenylation site, resulting in a truncated Pkhd1 transcript. A novel NOD congenic mouse expressing aberrant Pkhd1, but lacking the c3 and c4 chromosomal regions (NOD.Abd3), reproduces the immunopathological features of NOD ABD. RNA sequencing of NOD.Abd3 common bile duct early in disease demonstrates upregulation of genes involved in cholangiocyte injury/morphology and downregulation of immunoregulatory genes. Consistent with this, bone marrow chimera studies show that aberrant Pkhd1 must be expressed in the target tissue (cholangiocytes) and the immune system (bone marrow). Mutations of Pkhd1 produce biliary abnormalities in mice but have not been previously associated with autoimmunity. In this study, we eliminate clinical biliary disease by backcrossing this Pkhd1 mutation onto the C57BL/6 genetic background; thus, the NOD genetic background (which promotes autoimmunity) is essential for disease. We propose that loss of functional Pkhd1 on the NOD background produces early bile duct abnormalities, initiating a break in tolerance that leads to autoimmune cholangitis in NOD.Abd3 congenic mice. This model is important for understanding loss of tolerance to cholangiocytes and is relevant to the pathogenesis of several human cholangiopathies.

A core program of gene expression characterizes cancer metastases

While aberrant expression or splicing of metastasis genes conveys to cancers the ability to break through tissue barriers and disseminate, the genetic basis for organ preference in metastasis formation has remained incompletely understood. Utilizing the gene expression profiles from 653 GEO datasets, we investigate whether the signatures by diverse cancers in various metastatic sites display common features. We corroborate the meta-analysis in a murine model. Metastases are generally characterized by a core program of gene expression that induces the oxidative metabolism, activates vascularization/tissue remodeling, silences extracellular matrix interactions, and alters ion homeostasis. This program distinguishes metastases from their originating primary tumors as well as from their target host tissues. Site-selectivity is accomplished through specific components that adjust to the target micro-environment. The same functional groups of gene expression programs are activated in the metastases of B16-F10 cells to various target organs. It remains to be investigated whether these genetic signatures precede implantation and thus determine organ preference or are shaped by the target site and are thus a consequence of implantation. Conceivably, chemotherapy of disseminated cancer might be more efficacious if selected to match the genetic makeup of the metastases rather than the organ of origin by the primary tumor.

LncPRESS1 Is a p53-Regulated LncRNA that Safeguards Pluripotency by Disrupting SIRT6-Mediated De-acetylation of Histone H3K56

Recent evidence suggests that lncRNAs play an integral regulatory role in numerous functions, including determination of cellular identity. We determined global expression (RNA-seq) and genome-wide profiles (ChIP-seq) of histone post-translational modifications and p53 binding in human embryonic stem cells (hESCs) undergoing differentiation to define a high-confidence set of 40 lncRNAs, which are p53 transcriptional targets. We focused on lncRNAs highly expressed in pluripotent hESCs and repressed by p53 during differentiation to identify lncPRESS1 as a p53-regulated transcript that maintains hESC pluripotency in concert with core pluripotency factors. RNA-seq of hESCs depleted of lncPRESS1 revealed that lncPRESS1 controls a gene network that promotes pluripotency. Further, we found that lncPRESS1 physically interacts with SIRT6 and prevents SIRT6 chromatin localization, which maintains high levels of histone H3K56 and H3K9 acetylation at promoters of pluripotency genes. In summary, we describe a p53-regulated, pluripotency-specific lncRNA that safeguards the hESC state by disrupting SIRT6 activity.

Large scale adverse event data mining for targeted therapies development

2538

Background: Targeted anti-cancer small molecule drugs & immune therapies have had a dramatic impact in improving outcomes & the approach to clinical trials. Increasingly, regulatory approvals are expedited with small studies designed to identify strong efficacy signals. However, this may limit the extent of safety profiling. The use of large scale/big data meta-analyses can identify novel safety & efficacy signals in "real-world" medical settings. Methods: We used AERSMine, an open-source data mining platform to identify drug toxicity signatures in the FDA’s Adverse Event Reporting System of 8.6 million patients. We identified patients (n = 732,198) who received either traditional and targeted cancer therapy & identified therapy-specific toxicity patterns. Patients were classified based on exposures: anthracyclines (n = 83,179), platinum (117,993), antimetabolites (93,062), alkylators (81,507), antimicrotubule agents (97,726), HER2 inhibitors (40,040), VEGFis (79,144), VEGF-TKis (90,734), multi TKis (34,457), anaplastic lymphoma Kis (7,635), PI3K-AKT-mTOR inhibitors (33,864), Bruton TKis (9,247), MEKis (4,018), immunomodulatory agents (174,810), proteasome inhibitors (44,681), immune checkpoint inhibitors (20,287). Pharmacovigilance metrics [Relative Risks & safety signals] were used to establish statistical correlation & toxicity signatures were differentiated using the Kolmogorov–Smirnov test. Results: To validate the use of the AERSMine to detect AEs, we focused on cardiotoxicity. It identified classic drug associated AEs (e.g. ventricular dysfunction with anthracyclines, HER2is & VEGFis; VEGFi hypertension & vascular toxicity; multi TKIs vascular events). AERSMine also identified recently reported uncommon toxicities of myositis/myocarditis with immune checkpoint inhibitors. It indicated a higher frequency of myositis/myocarditis with combination immune checkpoint therapy, paralleling industry corporate safety databases. These toxicities were reported at higher frequencies in patients > 65 yrs. Conclusions: AERSMine “big data” analyses provide a sensitive tool to detect potential new patterns of AEs simultaneously across multiple clinical trials & in the real-world setting.

Lung Gene Expression Analysis (LGEA): an integrative web portal for comprehensive gene expression data analysis in lung development

'LungGENS', our previously developed web tool for mapping single-cell gene expression in the developing lung, has been well received by the pulmonary research community. With continued support from the 'LungMAP' consortium, we extended the scope of the LungGENS database to accommodate transcriptomics data from pulmonary tissues and cells from human and mouse at different stages of lung development. Lung Gene Expression Analysis (LGEA) web portal is an extended version of LungGENS useful for the analysis, display and interpretation of gene expression patterns obtained from single cells, sorted cell populations and whole lung tissues. The LGEA web portal is freely available at http://research.cchmc.org/pbge/lunggens/mainportal.html.

Abstract A06: Characterization of a novel p53-regulated embryonic stem cell-specific lncRNA that safeguards pluripotency

Recent work shows that long non-coding RNAs (lncRNAs), transcribed from 80% of our genomes, play important roles in gene regulation and are associated with disease pathogenesis. However, there is a broad gap in our knowledge regarding any contributions of lncRNAs in mediating p53 response or regulating stem cell state. Here, we report our discovery of lncRNAs that fine-tune p53's transcriptional outcomes in differentiating human embryonic stem cells (hESCs). By integrating RNA-Seq with p53 ChIP-Seq analyses along with active or repressive histone marks of hESCs undergoing differentiation, we defined a high-confidence set of 40 lncRNAs that are p53 transcriptional targets. We found that p53 directly regulates nuclear lncRNAs, which are either activated (HOTAIRM1) or repressed (lncPRESS-1) during hESC differentiation. Single-cell gene expression analysis during lineage-specific hESC differentiation revealed that lncPRESS-1 is a novel pluripotency-specific, p53-regulated transcript with expression tightly correlated with OCT4, NANOG and PRDM14. Loss-of-function of lncPRESS-1 results in diminished pluripotency and spontaneous differentiation of hESCs. In addition, mass cytometry-based multi-parametric analysis at a single cell level revealed that depletion of lncPRESS-1 results in overall loss of hESC pluripotency with marked induction of an ectoderm specific gene-signature. Further, we show that lncPRESS-1 physically interacts with SIRT6 and restrains its chromatin localization to maintain high levels of histone H3 Lys-56 (H3K56) and Lys-9 (H3K9) acetylation at promoters of pluripotency genes to sustain a chromatin architecture that favors pluripotency. In summary, we describe a novel pluripotency-specific lncRNA that safeguards the hESC state by disrupting SIRT6 activity. Since lncPRESS-1 is repressed by p53, we believe that cancers that harbor mutant TP53 will have robust expression of this lncRNA. Further investigations are underway to determine the disease relevance of lncPRESS-1 and correlation with p53-mutations.

Citation Format: Abhinav K. Jain, Kadir C. Akdemir, Yuanxin Xi, Ryan L. McCarthy, Kendra Allton, Bruce Aronow, Chunru Lin, Wei Li, Liuqing Yang, Michelle C. Barton. Characterization of a novel p53-regulated embryonic stem cell-specific lncRNA that safeguards pluripotency. [abstract]. In: Proceedings of the AACR Special Conference on Noncoding RNAs and Cancer: Mechanisms to Medicines ; 2015 Dec 4-7; Boston, MA. Philadelphia (PA): AACR; Cancer Res 2016;76(6 Suppl):Abstract nr A06.

Spatial genomic heterogeneity in diffuse intrinsic pontine and midline high-grade glioma: implications for diagnostic biopsy and targeted therapeutics

Introduction

Diffuse intrinsic pontine glioma (DIPG) and midline high-grade glioma (mHGG) are lethal childhood brain tumors. Spatial genomic heterogeneity has been well-described in adult HGG but has not been comprehensively characterized in pediatric HGG. We performed whole exome sequencing on 38-matched primary, contiguous, and metastatic tumor sites from eight children with DIPG (n = 7) or mHGG (n = 1) collected using a unique MRI-guided autopsy protocol. Validation was performed using Sanger sequencing, Droplet Digital polymerase-chain reaction, immunohistochemistry, and fluorescent in-situ hybridization.

Results

Median age at diagnosis was 6.1 years (range: 2.9–23.3 years). Median overall survival was 13.2 months (range: 11.2–32.2 months). Contiguous tumor infiltration and distant metastases were observed in seven and six patients, respectively, including leptomeningeal dissemination in three DIPGs. Histopathological heterogeneity was evident in seven patients, including intra-pontine heterogeneity in two DIPGs, ranging from World Health Organization grade II to IV astrocytoma. We found conservation of heterozygous K27M mutations in H3F3A (n = 4) or HIST1H3B (n = 3) across all primary, contiguous, and metastatic tumor sites in all DIPGs. ACVR1 (n = 2), PIK3CA (n = 2), FGFR1 (n = 2), and MET (n = 1) were also intra-tumorally conserved. ACVR1 was co-mutated with HIST1H3B (n = 2). In contrast, PDGFRA amplification and mutation were spatially heterogeneous, as were mutations in BCOR (n = 1), ATRX (n = 2), and MYC (n = 1). TP53 aberrations (n = 3 patients) varied by type and location between primary and metastatic tumors sites but were intra-tumorally conserved.

Conclusion

Spatial conservation of prognostically-relevant and therapeutically-targetable somatic mutations in DIPG and mHGG contrasts the significant heterogeneity of driver mutations seen in adult HGG and supports uniform implementation of diagnostic biopsy in DIPG and mHGG to classify molecular risk groups and guide therapeutic strategy.

Electronic supplementary material

The online version of this article (doi:10.1186/s40478-015-0269-0) contains supplementary material, which is available to authorized users.

Dissecting Allele Architecture of Early Onset IBD Using High-Density Genotyping

Background

The inflammatory bowel diseases (IBD) are common, complex disorders in which genetic and environmental factors are believed to interact leading to chronic inflammatory responses against the gut microbiota. Earlier genetic studies performed in mostly adult population of European descent identified 163 loci affecting IBD risk, but most have relatively modest effect sizes, and altogether explain only ~20% of the genetic susceptibility. Pediatric onset represents about 25% of overall incident cases in IBD, characterized by distinct disease physiology, course and risks. The goal of this study is to compare the allelic architecture of early onset IBD with adult onset in population of European descent.

Methods

We performed a fine mapping association study of early onset IBD using high-density Immunochip genotyping on 1008 pediatric-onset IBD cases (801 Crohn’s disease; 121 ulcerative colitis and 86 IBD undetermined) and 1633 healthy controls. Of the 158 SNP genotypes obtained (out of the 163 identified in adult onset), this study replicated 4% (5 SNPs out of 136) of the SNPs identified in the Crohn’s disease (CD) cases and 0.8% (1 SNP out of 128) in the ulcerative colitis (UC) cases. Replicated SNPs implicated the well known NOD2 and IL23R. The point estimate for the odds ratio (ORs) for NOD2 was above and outside the confidence intervals reported in adult onset. A polygenic liability score weakly predicted the age of onset for a larger collection of CD cases (p< 0.03, R²= 0.007), but not for the smaller number of UC cases.

Conclusions

The allelic architecture of common susceptibility variants for early onset IBD is similar to that of adult onset. This immunochip genotyping study failed to identify additional common variants that may explain the distinct phenotype that characterize early onset IBD. A comprehensive dissection of genetic loci is necessary to further characterize the genetic architecture of early onset IBD.

Divergence of transcriptional landscape occurs early in B cell activation

Background

Signaling via B cell receptor (BCR) and Toll-like receptors (TLRs) results in activation of B cells with distinct physiological outcomes, but transcriptional regulatory mechanisms that drive activation and distinguish these pathways remain unknown.

Results

Two hours after ligand exposure RNA-seq, ChIP-seq and computational methods reveal that BCR- or TLR-mediated activation of primary resting B cells proceeds via a large set of shared and a smaller subset of distinct signal-selective transcriptional responses. BCR stimulation resulted in increased global recruitment of RNA Pol II to promoters that appear to transit slowly to downstream regions. Conversely, lipopolysaccharide (LPS) stimulation involved an enhanced RNA Pol II transition from initiating to elongating mode accompanied by greater H3K4me3 activation markings compared to BCR stimulation. These rapidly diverging transcriptomic landscapes also show distinct repressing (H3K27me3) histone signatures, mutually exclusive transcription factor binding in promoters, and unique miRNA profiles.

Conclusions

Upon examination of genome-wide transcription and regulatory elements, we conclude that the B cell commitment to different activation states occurs much earlier than previously thought and involves a multi-faceted receptor-specific transcriptional landscape.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-015-0012-x) contains supplementary material, which is available to authorized users.

Deducing the stage of origin of Wilms' tumours from a developmental series of Wt1-mutant mice

Wilms' tumours, paediatric kidney cancers, are the archetypal example of tumours caused through the disruption of normal development. The genetically best-defined subgroup of Wilms' tumours is the group caused by biallelic loss of the WT1 tumour suppressor gene. Here, we describe a developmental series of mouse models with conditional loss of Wt1 in different stages of nephron development before and after the mesenchymal-to-epithelial transition (MET). We demonstrate that Wt1 is essential for normal development at all kidney developmental stages under study. Comparison of genome-wide expression data from the mutant mouse models with human tumour material of mutant or wild-type WT1 datasets identified the stage of origin of human WT1-mutant tumours, and emphasizes fundamental differences between the two human tumour groups due to different developmental stages of origin.

Summary: The comparison of different nephron-specific Wt1-knockout mouse models identifies the stage of origin of human WT1-mutant Wilms' tumours.

Construction and characterization of stable, constitutively expressed, chromosomal green and red fluorescent transcriptional fusions in the select agents, Bacillus anthracis, Yersinia pestis, Burkholderia mallei, and Burkholderia pseudomallei

Here, we constructed stable, chromosomal, constitutively expressed, green and red fluorescent protein (GFP and RFP) as reporters in the select agents, Bacillus anthracis, Yersinia pestis, Burkholderia mallei, and Burkholderia pseudomallei. Using bioinformatic approaches and other experimental analyses, we identified P0253 and P1 as potent promoters that drive the optimal expression of fluorescent reporters in single copy in B. anthracis and Burkholderia spp. as well as their surrogate strains, respectively. In comparison, Y. pestis and its surrogate strain need two chromosomal copies of cysZK promoter (P2cysZK) for optimal fluorescence. The P0253-, P2cysZK-, and P1-driven GFP and RFP fusions were first cloned into the vectors pRP1028, pUC18R6KT-mini-Tn7T-Km, pmini-Tn7-gat, or their derivatives. The resultant constructs were delivered into the respective surrogates and subsequently into the select agent strains. The chromosomal GFP- and RFP-tagged strains exhibited bright fluorescence at an exposure time of less than 200 msec and displayed the same virulence traits as their wild-type parental strains. The utility of the tagged strains was proven by the macrophage infection assays and lactate dehydrogenase release analysis. Such strains will be extremely useful in high-throughput screens for novel compounds that could either kill these organisms, or interfere with critical virulence processes in these important bioweapon agents and during infection of alveolar macrophages.