Full-length isoform concatenation sequencing to resolve cancer transcriptome complexity.

BACKGROUND: Cancers exhibit complex transcriptomes with aberrant splicing that induces isoform-level differential expression compared to non-diseased tissues. Transcriptomic profiling using short-read sequencing has utility in providing a cost-effective approach for evaluating isoform expression, although short-read assembly displays limitations in the accurate inference of full-length transcripts. Long-read RNA sequencing (Iso-Seq), using the Pacific Biosciences (PacBio) platform, can overcome such limitations by providing full-length isoform sequence resolution which requires no read assembly and represents native expressed transcripts. A constraint of the Iso-Seq protocol is due to fewer reads output per instrument run, which, as an example, can consequently affect the detection of lowly expressed transcripts. To address these deficiencies, we developed a concatenation workflow, PacBio Full-Length Isoform Concatemer Sequencing (PB_FLIC-Seq), designed to increase the number of unique, sequenced PacBio long-reads thereby improving overall detection of unique isoforms. In addition, we anticipate that the increase in read depth will help improve the detection of moderate to low-level expressed isoforms. RESULTS: In sequencing a commercial reference (Spike-In RNA Variants; SIRV) with known isoform complexity we demonstrated a 3.4-fold increase in read output per run and improved SIRV recall when using the PB_FLIC-Seq method compared to the same samples processed with the Iso-Seq protocol. We applied this protocol to a translational cancer case, also demonstrating the utility of the PB_FLIC-Seq method for identifying differential full-length isoform expression in a pediatric diffuse midline glioma compared to its adjacent non-malignant tissue. Our data analysis revealed increased expression of extracellular matrix (ECM) genes within the tumor sample, including an isoform of the Secreted Protein Acidic and Cysteine Rich (SPARC) gene that was expressed 11,676-fold higher than in the adjacent non-malignant tissue. Finally, by using the PB_FLIC-Seq method, we detected several cancer-specific novel isoforms. CONCLUSION: This work describes a concatenation-based methodology for increasing the number of sequenced full-length isoform reads on the PacBio platform, yielding improved discovery of expressed isoforms. We applied this workflow to profile the transcriptome of a pediatric diffuse midline glioma and adjacent non-malignant tissue. Our findings of cancer-specific novel isoform expression further highlight the importance of long-read sequencing for characterization of complex tumor transcriptomes.

EGFR internal tandem duplications in fusion-negative congenital and neonatal spindle cell tumors.

Next-generation sequencing (NGS) assays can sensitively detect somatic variation, and increasingly can enable the identification of complex structural rearrangements. A subset of infantile spindle cell sarcomas, particularly congenital mesoblastic nephromas with classic or mixed histology, have structural rearrangement in the form of internal tandem duplications (ITD) involving EGFR. We performed prospective analysis to identify EGFR ITD through clinical or research studies, as well as retrospective analysis to quantify the frequency of EGFR ITD in pediatric sarcomas. Within our institution, three tumors with EGFR ITD were prospectively identified, all occurring in patients less than 1 year of age at diagnosis, including two renal tumors and one mediastinal soft tissue tumor. These three cases exhibited both cellular and mixed cellular and classic histology. All patients had no evidence of disease progression off therapy, despite incomplete resection. To extend our analysis and quantify the frequency of EGFR ITD in pediatric sarcomas, we retrospectively analyzed a cohort of tumors (n = 90) that were previously negative for clinical RT-PCR-based fusion testing. We identified EGFR ITD in three analyzed cases, all in patients less than 1 year of age (n = 18; 3/18, 17%). Here we expand the spectrum of tumors with EGFR ITD to congenital soft tissue tumors and report an unusual example of an EGFR ITD in a tumor with cellular congenital mesoblastic nephroma histology. We also highlight the importance of appropriate test selection and bioinformatic analysis for identification of this genomic alteration that is unexpectedly common in congenital and infantile spindle cell tumors.

Pacific Biosciences Fusion and Long Isoform Pipeline for Cancer Transcriptome-Based Resolution of Isoform Complexity.

Genomic profiling using short-read sequencing has utility in detecting disease-associated variation in both DNA and RNA. However, given the frequent occurrence of structural variation in cancer, molecular profiling using long-read sequencing improves the resolution of such events. For example, the Pacific Biosciences long-read RNA-sequencing (Iso-Seq) transcriptome protocol provides full-length isoform characterization, discernment of allelic phasing, and isoform discovery, and identifies expressed fusion partners. The Pacific Biosciences Fusion and Long Isoform Pipeline (PB_FLIP) incorporates a suite of RNA-sequencing software analysis tools and scripts to identify expressed fusion partners and isoforms. In addition, sequencing of a commercial reference (Spike-In RNA Variants) with known isoform complexity was performed and demonstrated high recall of the Iso-Seq and PB_FLIP workflow to benchmark our protocol and analysis performance. This study describes the utility of Iso-Seq and PB_FLIP analysis in improving deconvolution of complex structural variants and isoform detection within an institutional pediatric and adolescent/young adult translational cancer research cohort. The exemplar case studies demonstrate that Iso-Seq and PB_FLIP discover novel expressed fusion partners, resolve complex intragenic alterations, and discriminate between allele-specific expression profiles.

Detection of brain somatic variation in epilepsy-associated developmental lesions.

OBJECTIVE: Epilepsy-associated developmental lesions, including malformations of cortical development and low-grade developmental tumors, represent a major cause of drug-resistant seizures requiring surgical intervention in children. Brain-restricted somatic mosaicism has been implicated in the genetic etiology of these lesions; however, many contributory genes remain unidentified. METHODS: We enrolled 50 children who were undergoing epilepsy surgery into a translational research study. Resected tissue was divided for clinical neuropathologic evaluation and genomic analysis. We performed exome and RNA sequencing to identify somatic variation and we confirmed our findings using high-depth targeted DNA sequencing. RESULTS: We uncovered candidate disease-causing somatic variation affecting 28 patients (56%), as well as candidate germline variants affecting 4 patients (8%). In agreement with previous studies, we identified somatic variation affecting solute carrier family 35 member A2 (SLC35A2) and mechanistic target of rapamycin kinase (MTOR) pathway genes in patients with focal cortical dysplasia. Somatic gains of chromosome 1q were detected in 30% (3 of 10) of patients with Type I focal cortical dysplasia (FCD)s. Somatic variation in mitogen-activated protein kinase (MAPK) pathway genes (i.e., fibroblast growth factor receptor 1 [FGFR1], FGFR2, B-raf proto-oncogene, serine/threonine kinase [BRAF], and KRAS proto-oncogene, GTPase [KRAS]) was associated with low-grade epilepsy-associated developmental tumors. RNA sequencing enabled the detection of somatic structural variation that would have otherwise been missed, and which accounted for more than one-half of epilepsy-associated tumor diagnoses. Sampling across multiple anatomic regions revealed that somatic variant allele fractions vary widely within epileptogenic tissue. Finally, we identified putative disease-causing variants in genes not yet associated with focal cortical dysplasia. SIGNIFICANCE: These results further elucidate the genetic basis of structural brain abnormalities leading to focal epilepsy in children and point to new candidate disease genes.

Murine cytomegalovirus promotes renal allograft inflammation via Th1/17 cells and IL-17A.

Human cytomegalovirus (HCMV) infection is associated with renal allograft failure. Allograft damage in animal models is accelerated by CMV-induced T helper 17 (Th17) cell infiltrates. However, the mechanisms whereby CMV promotes Th17 cell-mediated pathological organ inflammation are uncharacterized. Here we demonstrate that murine CMV (MCMV)-induced intragraft Th17 cells have a Th1/17 phenotype co-expressing IFN-γ and/or TNF-α, but only a minority of these cells are MCMV specific. Instead, MCMV promotes intragraft expression of CCL20 and CXCL10, which are associated with recruitment of CCR6+ CXCR3+ Th17 cells. MCMV also enhances Th17 cell infiltrates after ischemia-reperfusion injury, independent of allogeneic responses. Pharmacologic inhibition of the Th17 cell signature cytokine, IL-17A, ameliorates MCMV-associated allograft damage without increasing intragraft viral loads or reducing MCMV-specific Th1 cell infiltrates. Clinically, HCMV DNAemia is associated with higher serum IL-17A among renal transplant patients with acute rejection, linking HCMV reactivation with Th17 cell cytokine expression. In summary, CMV promotes allograft damage via cytokine-mediated Th1/17 cell recruitment, which may be pharmacologically targeted to mitigate graft injury while preserving antiviral T cell immunity.

Clinically aggressive pediatric spinal ependymoma with novel MYC amplification demonstrates molecular and histopathologic similarity to newly described MYCN-amplified spinal ependymomas.

Primary spinal cord tumors contribute to ≤ 10% of central nervous system tumors in individuals of pediatric or adolescent age. Among intramedullary tumors, spinal ependymomas make up ~ 30% of this rare tumor population. A twelve-year-old male presented with an intradural, extramedullary mass occupying the dorsal spinal canal from C6 through T2. Gross total resection and histopathology revealed a World Health Organization (WHO) grade 2 ependymoma. He recurred eleven months later with extension from C2 through T1-T2. Subtotal resection was achieved followed by focal proton beam irradiation and chemotherapy. Histopathology was consistent with WHO grade 3 ependymoma. Molecular profiling of the primary and recurrent tumors revealed a novel amplification of the MYC (8q24) gene, which was confirmed by fluorescence in situ hybridization studies. Although MYC amplification in spinal ependymoma is exceedingly rare, a newly described classification of spinal ependymoma harboring MYCN (2p24) amplification (SP-MYCN) has been defined by DNA methylation-array based profiling. These individuals typically present with a malignant progression and dismal outcomes, contrary to the universally excellent survival outcomes seen in other spinal ependymomas. DNA methylation array-based classification confidently classified this tumor as SP-MYCN ependymoma. Notably, among the cohort of 52 tumors comprising the SP-MYCN methylation class, none harbor MYC amplification, highlighting the rarity of this genomic amplification in spinal ependymoma. A literature review comparing our individual to reported SP-MYCN tumors (n = 26) revealed similarities in clinical, histopathologic, and molecular features. Thus, we provide evidence from a single case to support the inclusion of MYC amplified spinal ependymoma within the molecular subgroup of SP-MYCN.

Discovery of clinically relevant fusions in pediatric cancer.

BACKGROUND: Pediatric cancers typically have a distinct genomic landscape when compared to adult cancers and frequently carry somatic gene fusion events that alter gene expression and drive tumorigenesis. Sensitive and specific detection of gene fusions through the analysis of next-generation-based RNA sequencing (RNA-Seq) data is computationally challenging and may be confounded by low tumor cellularity or underlying genomic complexity. Furthermore, numerous computational tools are available to identify fusions from supporting RNA-Seq reads, yet each algorithm demonstrates unique variability in sensitivity and precision, and no clearly superior approach currently exists. To overcome these challenges, we have developed an ensemble fusion calling approach to increase the accuracy of identifying fusions. RESULTS: Our Ensemble Fusion (EnFusion) approach utilizes seven fusion calling algorithms: Arriba, CICERO, FusionMap, FusionCatcher, JAFFA, MapSplice, and STAR-Fusion, which are packaged as a fully automated pipeline using Docker and Amazon Web Services (AWS) serverless technology. This method uses paired end RNA-Seq sequence reads as input, and the output from each algorithm is examined to identify fusions detected by a consensus of at least three algorithms. These consensus fusion results are filtered by comparison to an internal database to remove likely artifactual fusions occurring at high frequencies in our internal cohort, while a "known fusion list" prevents failure to report known pathogenic events. We have employed the EnFusion pipeline on RNA-Seq data from 229 patients with pediatric cancer or blood disorders studied under an IRB-approved protocol. The samples consist of 138 central nervous system tumors, 73 solid tumors, and 18 hematologic malignancies or disorders. The combination of an ensemble fusion-calling pipeline and a knowledge-based filtering strategy identified 67 clinically relevant fusions among our cohort (diagnostic yield of 29.3%), including RBPMS-MET, BCAN-NTRK1, and TRIM22-BRAF fusions. Following clinical confirmation and reporting in the patient's medical record, both known and novel fusions provided medically meaningful information. CONCLUSIONS: The EnFusion pipeline offers a streamlined approach to discover fusions in cancer, at higher levels of sensitivity and accuracy than single algorithm methods. Furthermore, this method accurately identifies driver fusions in pediatric cancer, providing clinical impact by contributing evidence to diagnosis and, when appropriate, indicating targeted therapies.

Transcriptomic Profiles in Children With Septic Shock With or Without Immunoparalysis.

Background: Severe innate immune suppression, termed immunoparalysis, is associated with increased risks of nosocomial infection and mortality in children with septic shock. Currently, immunoparalysis cannot be clinically diagnosed in children, and mechanisms remain unclear. Transcriptomic studies identify subsets of septic children with downregulation of genes within adaptive immune pathways, but assays of immune function have not been performed as part of these studies, and little is known about transcriptomic profiles of children with immunoparalysis. Methods: We performed a nested case-control study to identify differences in RNA expression patterns between children with septic shock with immunoparalysis (defined as lipopolysaccharide (LPS)-induced tumor necrosis factor (TNF)α response < 200 pg/ml) vs those with normal LPS-induced TNFα response. Children were enrolled within 48 hours of the onset of septic shock and divided into two groups based on LPS-induced TNFα response. RNA was extracted from whole blood for RNAseq, differential expression analyses using DESeq2 software, and pathway analyses using Ingenuity Pathway Analysis. Results: 32 children were included in analyses. Comparing those with immunoparalysis (n =19) to those with normal TNFα response (n = 13), 2,303 transcripts were differentially expressed with absolute value fold change ≥ 1.5 and false discovery rate ≤ 0.05. The majority of downregulated pathways in children with immunoparalysis were pathways that involved interactions between innate and adaptive immune cells necessary for cell-mediated immunity, crosstalk between dendritic cells and natural killer cells, and natural killer cell signaling pathways. Upregulated pathways included those involved in humoral immunity (T helper cell type 2), corticotropin signaling, platelet activation (GP6 signaling), and leukocyte migration and extravasation. Conclusions: Our study suggests that gene expression data might be useful to identify children with immunoparalysis and identifies several key differentially regulated pathways involved in both innate and adaptive immunity. Our ongoing work in this area aims to dissect interactions between innate and adaptive immunity in septic children and to more fully elucidate patient-specific immunologic pathophysiology to guide individualized immunotherapeutic targets.

Defining the AHR-regulated transcriptome in NK cells reveals gene expression programs relevant to development and function.

The aryl hydrocarbon receptor (AHR) is a ligand-activated transcription factor that regulates cellular processes in cancer and immunity, including innate immune cell development and effector function. However, the transcriptional repertoire through which AHR mediates these effects remains largely unexplored. To elucidate the transcriptional elements directly regulated by AHR in natural killer (NK) cells, we performed RNA and chromatin immunoprecipitation sequencing on NK cells exposed to AHR agonist or antagonist. We show that mature peripheral blood NK cells lack AHR, but its expression is induced by Stat3 during interleukin-21-driven activation and proliferation, coincident with increased NCAM1 (CD56) expression resulting in a CD56bright phenotype. Compared with control conditions, NK cells expanded in the presence of the AHR antagonist, StemRegenin-1, were unaffected in proliferation or cytotoxicity, had no increase in NCAM1 transcription, and maintained the CD56dim phenotype. However, it showed altered expression of 1004 genes including those strongly associated with signaling pathways. In contrast, NK cells expanded in the presence of the AHR agonist, kynurenine, showed decreased cytotoxicity and altered expression of 97 genes including those strongly associated with oxidative stress and cellular metabolism. By overlaying these differentially expressed genes with AHR chromatin binding, we identified 160 genes directly regulated by AHR, including hallmark AHR targets AHRR and CYP1B1 and known regulators of phenotype, development, metabolism, and function such as NCAM1, KIT, NQO1, and TXN. In summary, we define the AHR transcriptome in NK cells, propose a model of AHR and Stat3 coregulation, and identify potential pathways that may be targeted to overcome AHR-mediated immune suppression.

Evidence of pioneer factor activity of an oncogenic fusion transcription factor.

Recent characterizations of pioneer transcription factors provide insights into their structures and patterns of chromatin recognition associated with their roles in cell fate commitment and transformation. Intersecting with these basic science concepts, identification of pioneer factors (PFs) fused together as driver translocations in childhood cancers raises questions of whether these fusions retain the fundamental ability to invade repressed chromatin, consistent with their monomeric PF constituents. This study defines the cellular and chromatin localization of PAX3-FOXO1, an oncogenic driver of childhood rhabdomyosarcoma (RMS), derived from a fusion of PFs. To quantitatively define its chromatin-targeting functions and capacity to drive epigenetic reprogramming, we developed a ChIP-seq workflow with per-cell normalization (pc-ChIP-seq). Our quantitative localization studies address structural variation in RMS genomes and reveal insights into inactive chromatin localization of PAX3-FOXO1. Taken together, our studies are consistent with pioneer function for a driver oncoprotein in RMS, with repressed chromatin binding and nucleosome-motif targeting.

Gastroblastoma with a novel EWSR1-CTBP1 fusion presenting in adolescence.

Gastroblastomas are rare tumors with a biphasic epithelioid/spindle cell morphology that typically present in early adulthood and have recurrent MALAT1-GLI1 fusions. We describe an adolescent patient with Wiskott-Aldrich syndrome who presented with a large submucosal gastric tumor with biphasic morphology. Despite histologic features consistent with gastroblastoma, a MALAT1-GLI1 fusion was not found in this patient's tumor; instead, comprehensive molecular profiling identified a novel EWSR1-CTBP1 fusion and no other significant genetic alterations. The tumor also overexpressed NOTCH and FGFR by RNA profiling. The novel fusion and expression profile suggest a role for epithelial-mesenchymal transition in this tumor, with potential implications for the pathogenesis of biphasic gastric tumors such as gastroblastoma.

PTEN somatic mutations contribute to spectrum of cerebral overgrowth.

Phosphatase and tensin homologue (PTEN) regulates cell growth and survival through inhibition of the mammalian target of rapamycin (MTOR) signalling pathway. Germline genetic variation of PTEN is associated with autism, macrocephaly and PTEN hamartoma tumour syndromes. The effect of developmental PTEN somatic mutations on nervous system phenotypes is not well understood, although brain somatic mosaicism of MTOR pathway genes is an emerging cause of cortical dysplasia and epilepsy in the paediatric population. Here we report two somatic variants of PTEN affecting a single patient presenting with intractable epilepsy and hemimegalencephaly that varied in clinical severity throughout the left cerebral hemisphere. High-throughput sequencing analysis of affected brain tissue identified two somatic variants in PTEN. The first variant was present in multiple cell lineages throughout the entire hemisphere and associated with mild cerebral overgrowth. The second variant was restricted to posterior brain regions and affected the opposite PTEN allele, resulting in a segmental region of more severe malformation, and the only neurons in which it was found by single-nuclei RNA-sequencing had a unique disease-related expression profile. This study reveals brain mosaicism of PTEN as a disease mechanism of hemimegalencephaly and furthermore demonstrates the varying effects of single- or bi-allelic disruption of PTEN on cortical phenotypes.

Novel morphologic findings in PLAG1-rearranged soft tissue tumors.

Oncogenesis in PLAG1-rearranged tumors often results from PLAG1 transcription factor overexpression driven by promoter-swapping between constitutively expressed fusion partners. PLAG1-rearranged tumors demonstrate diverse morphologies. This study adds to this morphologic heterogeneity by introducing two tumors with PLAG1 rearrangements that display distinct histologic features. The first arose in the inguinal region of a 3-year-old, appeared well-circumscribed with a multinodular pattern, and harbored two fusions: ZFHX4-PLAG1 and CHCHD7-PLAG1. The second arose in the pelvic cavity of a 15-year-old girl, was extensively infiltrative and vascularized with an adipocytic component, and demonstrated a COL3A1-PLAG1 fusion. Both showed low-grade cytomorphology, scarce mitoses, no necrosis, and expression of CD34 and desmin. The ZFHX4-/CHCHD7-PLAG1-rearranged tumor showed no evidence of recurrence after 5 months. By contrast, the COL3A1-PLAG1-rearranged tumor quickly recurred following primary excision with positive margins; subsequent re-excision with adjuvant chemotherapy resulted in no evidence of recurrence after 2 years. While both tumors show overlap with benign and malignant fibroblastic and fibrovascular neoplasms, they also display divergent features. These cases highlight the importance of appropriate characterization in soft tissue tumors with unusual clinical and histologic characteristics.

The Sensor Kinase QseC Regulates the Unlinked PmrA Response Regulator and Downstream Gene Expression in Francisella.

The facultative intracellular bacterial pathogen Francisella tularensis is the causative agent of tularemia in humans and animals. Gram-negative bacteria utilize two-component regulatory systems (TCS) to sense and respond to their changing environment. No classical, tandemly arranged sensor kinase and response regulator TCS genes exist in the human virulent Francisella tularensis subsp. tularensis, but orphaned members are present. PmrA is an orphan response regulator responsible for intramacrophage growth and virulence; however, the regulation of PmrA activity is not understood. We and others have shown that PmrA represses the expression of priM, described to encode an antivirulence determinant. By screening a mutant library for increased priM promoter activity, we identified the sensor kinase homolog QseC as an upstream regulator of priM expression, and this regulation is in part dependent upon the aspartate phosphorylation site of PmrA (D51). Several examined environmental signals, including epinephrine, which is reported to activate QseC in other bacteria, do not affect priM expression in a manner dependent on PmrA. Intramacrophage survival assays also question the finding that PriM is an antivirulence factor. Thus, these data suggest that the PmrA-regulated gene priM is modulated by the QseC-PmrA (QseB) TCS in FrancisellaIMPORTANCE The disease tularemia is caused by the highly infectious Gram-negative pathogen Francisella tularensis This bacterium encodes few regulatory factors (e.g., two-component systems [TCS]). PmrA, required for intramacrophage survival and virulence in the mouse model, is encoded by an orphan TCS response regulator gene. It is unclear how PmrA is responsive to environmental signals to regulate loci, including the PmrA-repressed gene priM We identify an orphan sensor kinase (QseC) that is required for priM repression and further explore both environmental signals that might regulate the QseC-PmrA TCS and the function of PriM.

Whole-blood transcriptomic responses to lumacaftor/ivacaftor therapy in cystic fibrosis.

BACKGROUND: Cystic fibrosis (CF) remains without a definitive cure. Novel therapeutics targeting the causative defect in the cystic fibrosis transmembrane conductance regulator (CFTR) gene are in clinical use. Lumacaftor/ivacaftor is a CFTR modulator approved for patients homozygous for the CFTR variant p.Phe508del, but there are wide variations in treatment responses preventing prediction of patient responses. We aimed to determine changes in gene expression related to treatment initiation and response. METHODS: Whole-blood transcriptomics was performed using RNA-Seq in 20 patients with CF pre- and 6 months post-lumacaftor/ivacaftor (drug) initiation and 20 non-CF healthy controls. Correlation of gene expression with clinical variables was performed by stratification via clinical responses. RESULTS: We identified 491 genes that were differentially expressed in CF patients (pre-drug) compared with non-CF controls and 36 genes when comparing pre-drug to post-drug profiles. Both pre- and post-drug CF profiles were associated with marked overexpression of inflammation-related genes and apoptosis genes, and significant under-expression of T cell and NK cell-related genes compared to non-CF. CF patients post-drug demonstrated normalized protein synthesis expression, and decreased expression of cell-death genes compared to pre-drug profiles, irrespective of clinical response. However, CF clinical responders demonstrated changes in eIF2 signaling, oxidative phosphorylation, IL-17 signaling, and mitochondrial function compared to non-responders. Top overexpressed genes (MMP9 and SOCS3) that decreased post-drug were validated by qRT-PCR. Functional assays demonstrated that CF monocytes normalized calcium (increases MMP9 expression) concentrations post-drug. CONCLUSIONS: Transcriptomics revealed differentially regulated pathways in CF patients at baseline compared to non-CF, and in clinical responders to lumacaftor/ivacaftor.

Samovar: Single-Sample Mosaic Single-Nucleotide Variant Calling with Linked Reads.

Linked-read sequencing enables greatly improves haplotype assembly over standard paired-end analysis. The detection of mosaic single-nucleotide variants benefits from haplotype assembly when the model is informed by the mapping between constituent reads and linked reads. Samovar evaluates haplotype-discordant reads identified through linked-read sequencing, thus enabling phasing and mosaic variant detection across the entire genome. Samovar trains a random forest model to score candidate sites using a dataset that considers read quality, phasing, and linked-read characteristics. Samovar calls mosaic single-nucleotide variants (SNVs) within a single sample with accuracy comparable with what previously required trios or matched tumor/normal pairs and outperforms single-sample mosaic variant callers at minor allele frequency 5%-50% with at least 30X coverage. Samovar finds somatic variants in both tumor and normal whole-genome sequencing from 13 pediatric cancer cases that can be corroborated with high recall with whole exome sequencing. Samovar is available open-source at https://github.com/cdarby/samovar under the MIT license.

Whole Transcriptome Analysis of Renal Intercalated Cells Predicts Lipopolysaccharide Mediated Inhibition of Retinoid X Receptor alpha Function.

The renal collecting duct consists of intercalated cells (ICs) and principal cells (PCs). We have previously demonstrated that collecting ducts have a role in the innate immune defense of the kidney. Transcriptomics is an important tool used to enhance systems-level understanding of cell biology. However, transcriptomics performed on whole kidneys provides limited insight of collecting duct cell gene expression, because these cells comprise a small fraction of total kidney cells. Recently we generated reporter mouse models to enrich collecting duct specific PC and ICs and reported targeted gene expression of anti-microbial peptide genes. Here we report transcriptomics on enriched ICs and PCs and performed a pilot study sequencing four single ICs. We identified 3,645 genes with increased relative expression in ICs compared to non-ICs. In comparison to non-PCs, 2,088 genes had higher relative expression in PCs. IC associated genes included the innate interleukin 1 receptor, type 1 and the antimicrobial peptide(AMP) adrenomedullin. The top predicted canonical pathway for enriched ICs was lipopolysaccharide/Interleukin 1 mediated inhibition of Retinoid X Receptor alpha function and decreased Retinoid X Receptor expression was confirmed to occur 1-hour post experimental murine UTI in ICs but not in non-ICs.

Genome sequencing identifies somatic BRAF duplication c.1794_1796dupTAC;p.Thr599dup in pediatric patient with low-grade ganglioglioma.

Gangliogliomas (WHO grade I) are rare tumors affecting the central nervous system and are most frequently observed in children. Next-generation sequencing of tumors is being utilized at an increasing rate in both research and clinical settings to characterize the genetic factors that drive tumorigenesis. Here, we report a rare BRAF somatic mutation (NM_004333.4:c.1794_1796dupTAC; p.Thr599dup) in the tumor genome from a pediatric patient in her late teens, who was initially diagnosed with low-grade ganglioglioma at age 13. This duplication of 3 nt introduces a second threonine residue at amino acid 599 of the BRAF protein. Based on previous studies, this variant is likely to increase kinase activity, similar to the well-characterized BRAF p.Val600Glu (V600E) pathogenic variant. In addition, although the p.T599dup somatic mutation has been documented rarely in human cancers, the variant has not been previously reported in ganglioglioma. The identification of this variant presents an opportunity to consider targeted therapy (e.g., BRAF inhibitor) for this patient.

Safety Profile of Good Manufacturing Practice Manufactured Interferon γ-Primed Mesenchymal Stem/Stromal Cells for Clinical Trials.

Mesenchymal stem/stromal cells (MSCs) are widely studied by both academia and industry for a broad array of clinical indications. The collective body of data provides compelling evidence of the clinical safety of MSC therapy. However, generally accepted proof of therapeutic efficacy has not yet been reported. In an effort to generate a more effective therapeutic cell product, investigators are focused on modifying MSC processing protocols to enhance the intrinsic biologic activity. Here, we report a Good Manufacturing Practice-compliant two-step MSC manufacturing protocol to generate MSCs or interferon γ (IFNγ) primed MSCs which allows freshly expanded cells to be infused in patients on a predetermined schedule. This protocol eliminates the need to infuse cryopreserved, just thawed cells which may reduce the immune modulatory activity. Moreover, using (IFNγ) as a prototypic cytokine, we demonstrate the feasibility of priming the cells with any biologic agent. We then characterized MSCs and IFNγ primed MSCs prepared with our protocol, by karyotype, in vitro potential for malignant transformation, biodistribution, effect on engraftment of transplanted hematopoietic cells, and in vivo toxicity in immune deficient mice including a complete post-mortem examination. We found no evidence of toxicity attributable to the MSC or IFNγ primed MSCs. Our data suggest that the clinical risk of infusing MSCs or IFNγ primed MSCs produced by our two-step protocol is not greater than MSCs currently in practice. While actual proof of safety requires phase I clinical trials, our data support the use of either cell product in new clinical studies. Stem Cells Translational Medicine 2017;6:1868-1879.