Search

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.

Survival of patients with acute myeloid leukemia (AML) is inversely associated with age, but the impact of race on outcomes of adolescent and young adult (AYA; range, 18-39 years) patients is unknown. We compared survival of 89 non-Hispanic Black and 566 non-Hispanic White AYA patients with AML treated on frontline Cancer and Leukemia Group B/Alliance for Clinical Trials in Oncology protocols. Samples of 327 patients (50 Black and 277 White) were analyzed via targeted sequencing. Integrated genomic profiling was performed on select longitudinal samples. Black patients had worse outcomes, especially those aged 18 to 29 years, who had a higher early death rate (16% vs 3%; P=.002), lower complete remission rate (66% vs 83%; P=.01), and decreased overall survival (OS; 5-year rates: 22% vs 51%; P

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.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.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.These results further elucidate the genetic basis of structural brain abnormalities leading to focal epilepsy in children and point to new candidate disease genes.

The AURORA US Metastasis Project was established with the goal to identify molecular features associated with metastasis. We assayed 55 females with metastatic breast cancer (51 primary cancers and 102 metastases) by RNA sequencing, tumor/germline DNA exome and low-pass whole-genome sequencing and global DNA methylation microarrays. Expression subtype changes were observed in ~30% of samples and were coincident with DNA clonality shifts, especially involving HER2. Downregulation of estrogen receptor (ER)-mediated cell-cell adhesion genes through DNA methylation mechanisms was observed in metastases. Microenvironment differences varied according to tumor subtype; the ER

Background: It has become increasingly clear that effective treatment of metastatic breast cancer (MBC) requires an in-depth understanding of the molecular differences between primary tumors and metastases. The AURORA US Network was established to collect primary breast cancer-metastasis pairs for multi-platform genomic profiling in order to identify the molecular drivers of metastatic disease. AURORA US has both a retrospective and prospective phase. This is the first report of the retrospective phase. Methods: Archived tissue samples from the primary tumor and at least one distant metastasis were retrospectively collected from 83 MBC patients. Following internal quality assessment, samples from 55 pts, including 105 distinct metastatic lesions, were subject to DNA low pass whole genome and exome sequencing, DNA methylation arrays, and RNA sequencing. Early analyses of these multi-platform data include: DNA methylation, tumor gene expression and microenvironmental signatures, somatic and germline variants, DNA copy number changes, and structural variants between breast primaries and matched metastases. Results: Median age at diagnosis was 49 years (25-76); 32 (58%) were Stage I or II at presentation, 27 (49%) had a family history of breast cancer, and 20 (36%) had a second breast primary. Median disease-free interval before developing MBC was 2 years (range 0-36, 5 patients presented with Stage IV). Median overall survival from initial presentation was 4 years (range 0-37). Median survival after developing MBC was 1 year (range 0-13), with a median of three treatments. Primary tissue samples were banked from 1977-2017 and metastases were banked from 1999-2017. Clinical phenotypes of the primaries included 27 HR+ (49%), 15 triple negative (TNBC, 27%), and 11 HER2+ (20%, 12 missing HER2 status). Intrinsic subtype distribution of the primaries included 17 Basal-like (31%), 17 Luminal A (31%), 7 Normal-like (13%), 5 HER2-enriched (9%), and 1 Luminal B, with 8 pending. All metastases from the Basal-like cases remained Basal-like, while metastases from luminal primaries tended to gain HER2-Enriched subtype features (5/18, p = 0.01). Overall, we identified significant metastasis-enriched alterations in metabolism pathways, an increase in proliferation, and the loss of differentiation signatures and immune infiltrates with progression; the latter being the most pronounced in brain metastases. The most frequent somatic mutations in this cohort were in TP53, NCOR1, and RUNX1. Interestingly, ERBB2, EGFR, and ATM were also mutated in ≥10% of the tumors sequenced. In almost all cases, CpG island hypermethylation was clonally present in the primary tumor and persisted stably in the majority of metastatic lesions. Promoter CpG island hypermethylation was also identified in some metastatic lesions at JAM3, an important cellular adhesion molecule,and this was accompanied by reduced mRNA expression. CONCLUSIONS: Collection of banked primary and metastatic tissue pairs identified a young MBC cohort with a high frequency of breast cancer family history and second breast primaries. Molecular characterization of luminal tumor pairs highlighted acquisition of aggressive traits including increased proliferation and loss of differentiation in the metastases. In contrast, basal-like pairs remained relatively unchanged, except for the loss of immune activation. Ongoing analyses to be presented include clonal heterogeneity and phylogeny, novel metastasis signature discovery, gene fusion, and endogenous retrovirus detection. Citation Format: Tari A King, Minetta C Liu, Marni B McClure, Toshinori Hinoue, Benjamin J Kelly, Chad J Creighton, Jay Bowen, Kristen Leraas, Robyn T Burns, Sara Coppens, Salma Rezk, Amy L Garrett, Justin M Balko, Joel S Parker, Ben H Park, Ian Krop, Carey Anders, Katherine A Hoadley, Julie Gastier-Foster, Mothaffar F Rimawi, Rita Nanda, Nancy U Lin, Claudine Isaacs, P. Kelly Marcom, Anna Maria Storniolo, Fergus J Couch, Elaine R Mardis, Adrian V Lee, Uma Chandran, Peter W Laird, Susan G Hilsenbeck, Larry Norton, Andrea L Richardson, W. Fraser Symmans, Lisa A Carey, Antonio C Wolff, Nancy E Davidson, Charles M Perou, the AURORA US Network. Multiplatform analysis of matched primary and metastatic breast tumors from the AURORA US Network [abstract]. In: Proceedings of the 2019 San Antonio Breast Cancer Symposium; 2019 Dec 10-14; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2020;80(4 Suppl):Abstract nr GS3-08.

Summary Osteosarcoma remains one of the deadliest cancers in pediatrics and young adults. We administered two types of immunotherapies, oncolytic virotherapy and immune checkpoint inhibition, to two murine osteosarcoma models and observed divergent results. Mice bearing F420 showed no response, whereas those with K7M2 showed prolonged survival in response to combination therapy. K7M2 had higher expression of immune-related genes and higher baseline immune cell infiltrates, but there were no significant differences in tumor mutational burden or predicted MHC class I binding of nonsynonymous mutations. Instead, we found several mouse endogenous retrovirus sequences highly expressed in K7M2 compared with F420. T cell tetramer staining for one of them, gp70, was detected in mice with K7M2 but not F420, suggesting that endogenous retrovirus proteins are targets for the anti-tumor immune reaction. Given prior observations of endogenous retrovirus expression in human osteosarcomas, our findings may be translatable to human disease., Graphical abstract, Highlights • K7M2 but not F420 responds to viroimmunotherapy • K7M2 and F420 exhibit similar tumor mutational burdens • K7M2 shows much higher expression of endogenous retrovirus genes including gp70 • Tetramer staining suggest gp70 is the immune target in K7M2, Cancer; Immunology

RLIM, also known as RNF12, is an X-linked E3 ubiquitin ligase acting as a negative regulator of LIM-domain containing transcription factors and participates in X-chromosome inactivation (XCI) in mice. We report the genetic and clinical findings of 84 individuals from nine unrelated families, eight of whom who have pathogenic variants in RLIM (RING finger LIM domain-interacting protein). A total of 40 affected males have X-linked intellectual disability (XLID) and variable behavioral anomalies with or without congenital malformations. In contrast, 44 heterozygous female carriers have normal cognition and behavior, but eight showed mild physical features. All RLIM variants identified are missense changes co-segregating with the phenotype and predicted to affect protein function. Eight of the nine altered amino acids are conserved and lie either within a domain essential for binding interacting proteins or in the C-terminal RING finger catalytic domain. In vitro experiments revealed that these amino acid changes in the RLIM RING finger impaired RLIM ubiquitin ligase activity. In vivo experiments in rlim mutant zebrafish showed that wild type RLIM rescued the zebrafish rlim phenotype, whereas the patient-specific missense RLIM variants failed to rescue the phenotype and thus represent likely severe loss-of-function mutations. In summary, we identified a spectrum of RLIM missense variants causing syndromic XLID and affecting the ubiquitin ligase activity of RLIM, suggesting that enzymatic activity of RLIM is required for normal development, cognition and behavior.

Congenital heart disease (CHD) is a common group of birth defects with a strong genetic contribution to their etiology, but historically the diagnostic yield from exome studies of isolated CHD has been low. Pleiotropy, variable expressivity, and the difficulty of accurately phenotyping newborns contribute to this problem. We hypothesized that performing exome sequencing on selected individuals in families with multiple members affected by left-sided CHD, then filtering variants by population frequency, in silico predictive algorithms, and phenotypic annotations from publicly available databases would increase this yield and generate a list of candidate disease-causing variants that would show a high validation rate. In eight of the nineteen families in our study (42%), we established a well-known gene/phenotype link for a candidate variant or performed confirmation of a candidate variant’s effect on protein function, including variants in genes not previously described or firmly established as disease genes in the body of CHD literature: BMP10, CASZ1, ROCK1 and SMYD1. Two plausible variants in different genes were found to segregate in the same family in two instances suggesting oligogenic inheritance. These results highlight the need for functional validation and demonstrate that in the era of next-generation sequencing, multiplex families with isolated CHD can still bring high yield to the discovery of novel disease genes.

Background: Metastatic breast cancer (MBC) patients often have short survival, and successful treatment represents one of the most challenging aspects of cancer care. This poor prognosis is likely multifactorial, including increased clonal heterogeneity, drug resistance mechanisms, and alterations of the tumor microenvironment. Methods: The primary data source was multi-platform data coming from the AURORA US Project that includes RNAseq, DNAseq, and DNA methylation arrays, which assayed 55 MBC patients representing 51 primary cancers and 102 linked metastatic specimens. In addition, RNA sequencing data from two other datasets of primary tumor-metastasis pairs were also used (i.e. UNC Tumor Donation Program/RAP dataset (24 primary and 74 metastasis specimens), and the GEICAM/2009-03 ConvertHER trial dataset (102 primaries and 102 metastatic pairs)). In total, this combined RNAseq dataset contained 177 primary tumors and 278 metastases, including 28 liver, 18 lung, 12 brain, and 24 lymph node metastases. We used these data to address two pressing questions, namely: 1) do gene expression features vary in primary tumors vs metastases according to PAM50 expression subtype, and 2) do gene expression features vary according to site of metastasis. Results: Using the AURORA multi-platform data, we determined that 17% of metastatic tumors (mainly TNBC/Basal-like) showed reduced expression of HLA-A that was associated with DNA methylation and/or focal DNA deletions near the HLA-A locus; these methylated tumors also showed concomitant lower immune cell infiltrates. Reduced expression of HLA-A gene and immune cell infiltrates were also validated at the RNA level in RAP dataset. Next, RNA expression differences were examined using the combined data set and varied according to tumor subtype. ER+/Luminal metastases had lower fibroblast and endothelial cell content, while triple negative (TNBC)/Basal-like metastases showed a dramatic decrease in T cell and B cell signatures/features. Comparative analyses between primary and site-specific metastasis (i.e., primary vs liver metastasis) or between sites of metastases (i.e., liver vs lung metastasis) revealed that both liver and brain, on average, had low immune cell features regardless of the primary tumor phenotype. Even within the same patient, we detected low immune cell features in brain and liver metastases compared to lung and lymph node metastases. Lastly, liver metastases showed a gain of Luminal B/HER2E gene expression features and MYC targets, and brain TNBC metastasis showed a gain of cell differentiation/Luminal-related gene signatures. Conclusions: These findings could have direct implications for the treatment of MBC patients with immune-based therapies and suggest new therapeutic avenues depending upon the tumor metastasis phenotype, and site of metastasis. Citation Format: Susana Garcia Recio, Toshinori Hinoue, Gregory L. Wheeler, Benjamin J. Kelly, Justin M. Balko, Katherine A. Hoadley, Peter W. Laird, Elaine R. Mardis, Charles M. Perou. Multiplatform analysis of matched primary and metastatic breast tumors from the AURORA US Network identifies microenvironment features as drivers of metastasis [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr LB176.

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.

Background The role of synonymous single-nucleotide variants in human health and disease is poorly understood, yet evidence suggests that this class of “silent” genetic variation plays multiple regulatory roles in both transcription and translation. One mechanism by which synonymous codons direct and modulate the translational process is through alteration of the elaborate structure formed by single-stranded mRNA molecules. While tools to computationally predict the effect of non-synonymous variants on protein structure are plentiful, analogous tools to systematically assess how synonymous variants might disrupt mRNA structure are lacking. Results We developed novel software using a parallel processing framework for large-scale generation of secondary RNA structures and folding statistics for the transcriptome of any species. Focusing our analysis on the human transcriptome, we calculated 5 billion RNA-folding statistics for 469 million single-nucleotide variants in 45,800 transcripts. By considering the impact of all possible synonymous variants globally, we discover that synonymous variants predicted to disrupt mRNA structure have significantly lower rates of incidence in the human population. Conclusions These findings support the hypothesis that synonymous variants may play a role in genetic disorders due to their effects on mRNA structure. To evaluate the potential pathogenic impact of synonymous variants, we provide RNA stability, edge distance, and diversity metrics for every nucleotide in the human transcriptome and introduce a “Structural Predictivity Index” (SPI) to quantify structural constraint operating on any synonymous variant. Because no single RNA-folding metric can capture the diversity of mechanisms by which a variant could alter secondary mRNA structure, we generated a SUmmarized RNA Folding (SURF) metric to provide a single measurement to predict the impact of secondary structure altering variants in human genetic studies.

BackgroundPediatric 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.ResultsOur ensemble fusion detection 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 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 ensemble fusion-calling 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.ConclusionsOur ensemble fusion detection 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.

Retinoblastoma is a childhood cancer of the retina involving germline or somatic alterations of the RB Transcriptional Corepressor 1 gene, RB1. Rare cases of sellar-suprasellar region retinoblastoma without evidence of ocular or pineal tumors have been described. A nine-month-old male presented with a sellar-suprasellar region mass. Histopathology showed an embryonal tumor with focal Flexner-Wintersteiner-like rosettes and loss of retinoblastoma protein (RB1) expression by immunohistochemistry. DNA array-based methylation profiling confidently classified the tumor as pineoblastoma group A/intracranial retinoblastoma. The patient was subsequently enrolled on an institutional translational cancer research protocol and underwent comprehensive molecular profiling, including paired tumor/normal exome and genome sequencing and RNA-sequencing of the tumor. Additionally, Pacific Biosciences (PacBio) Single Molecule Real Time (SMRT) sequencing was performed from comparator normal and disease-involved tissue to resolve complex structural variations. RNA-sequencing revealed multiple fusions clustered within 13q14.1-q21.3, including a novel in-frame fusion of RB1-SIAH3 predicted to prematurely truncate the RB1 protein. SMRT sequencing revealed a complex structural rearrangement spanning 13q14.11-q31.3, including two somatic structural variants within intron 17 of RB1. These events corresponded to the RB1-SIAH3 fusion and a novel RB1 rearrangement expected to correlate with the complete absence of RB1 protein expression. Comprehensive molecular analysis, including DNA array-based methylation profiling and sequencing-based methodologies, were critical for classification and understanding the complex mechanism of RB1 inactivation in this diagnostically challenging tumor.

SummaryRecent characterizations of pioneer transcription factors have led to new insights into their structures and patterns of chromatin recognition that are instructive for understanding their role in cell fate commitment and transformation. Intersecting with these basic science concepts, the 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. In this study, we define the cellular and chromatin localization of the translocation, PAX3-FOXO1, an oncogenic driver of childhood rhabdomyosarcoma (RMS), derived from a genetic fusion of PFs. To quantitatively define its chromatin-targeting functions and capacity to drive epigenetic reprogramming, we developed a new method for ChIP-seq with per-cell normalization (pc-ChIP-seq). Our quantitative localization studies address structural variation in RMS genomes and reveal novel insights into heterochromatin localization of PAX3-FOXO1. From these studies, we report novel pioneer function for the major driver oncogene in RMS, with repressed chromatin binding and nucleosome-motif targeting in human cells.

Meningiomas are a central nervous system tumor primarily afflicting adults, with

Metastatic disease at initial presentation of intracranial ependymoma is an uncommon occurrence with only rare reports of survival and is reportedly more prevalent in the youngest of children. Clinical and molecular characteristics associated with metastatic presentation, their prognostic implications, as well as optimal treatment options for such patients, have not been identified. CASE REPORT: A seven months old child presented with posterior fossa anaplastic ependymoma; following sub-total resection of primary tumor, a spine MRI revealed leptomeningeal dissemination along the cervical spinal cord and nerve roots of the cauda equina. The patient was successfully treated with five cycles of intensive induction chemotherapy (as per Head Start with high-dose methotrexate) followed by three sequential cycles of marrow-ablative chemotherapy and autologous hematopoietic progenitor cell rescue (AuHPCR) without irradiation; he is currently without evidence of the disease now 60 months following initial diagnosis. MOLECULAR/GENOMIC RESULTS: The patient was enrolled on a patient-centric comprehensive molecular profiling protocol, which included paired tumor-normal whole-exome sequencing, RNA sequencing of the disease-involved tissue, and DNA methylation classification. The genomic profile of the tumor was relatively unremarkable, revealing only a terminal gain of chromosome 3p and a terminal deletion of chromosome 22q, suggestive of an unbalanced translocation. Using RNA sequencing, we identified a novel SPECC1L-RAF1 gene fusion. The tumor harbors unique transcriptomic and DNA methylation profiles, failing to discretely classify with well-established ependymoma subgroups. CONCLUSION: Use of genomic profiling techniques provides meaningful information for disease characterization allowing for further expansion of the molecular spectrum associated with malignant disease.

Summary 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., Graphical abstract, Highlights • The fusion oncoprotein PAX3-FOXO1 binds to both active and repressed chromatin • PAX3-FOXO1-binding sites are adjacent to H3K9me3 domains • PAX3-FOXO1 engages partial DNA motifs at early timepoints • PAX3-FOXO1 can bind stably to inaccessible chromatin without inducing accessibility, Molecular Genetics; Systems biology; Transcriptomics

BackgroundIn most organisms the structure of an mRNA molecule is crucial in determining speed of translation, half-life, splicing propensities and final protein configuration. Synonymous variants which distort this wildtype mRNA structure may be pathogenic as a consequence. However, current clinical guidelines classify synonymous or “silent” single nucleotide variants (sSNVs) as largely benign unless a role in RNA splicing can be demonstrated.ResultsWe developed novel software to conduct a global transcriptome study in which RNA folding statistics were computed for 469 million SNVs in 45,800 transcripts using an Apache Spark implementation of ViennaRNA in the cloud. Focusing our analysis on the subset of 17.9 million sSNVs, we discover that variants predicted to disrupt mRNA structure have lower rates of incidence in the human population. Given that the community lacks tools to evaluate the potential pathogenic impact of sSNVs, we introduce a “Structural Predictivity Index” (SPI) to quantify this constraint due to mRNA structure.ConclusionsOur findings support the hypothesis that sSNVs may play a role in genetic disorders due to their effects on mRNA structure. Our RNA-folding scores provide a means of gauging the structural constraint operating on any sSNV in the human genome. Given that the majority of patients with rare or as yet to be diagnosed disease lack a molecular diagnosis, these scores have the potential to enable discovery of novel genetic etiologies. Our RNA Stability Pipeline as well as ViennaRNA structural metrics and SPI scores for all human synonymous variants can be downloaded from GitHub https://github.com/nch-igm/rna-stability.

Summary 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., Graphical Abstract, Highlights • Samovar uses haplotype-specific features from linked reads to call mosaic variants • Samovar quickly evaluates candidates with a random forest over 33 features • Only one sample is needed, with accuracy comparable with paired samples or trios • Samovar finds somatic variants in cancer driving genes in 13 pediatric cancer cases, Biological Sciences; Genomics; Bioinformatics

Klippel-Feil syndrome (KFS) is an exceerlingly rare constitutional disorder in which a paucity of knowledge exists about the disease and its associated morbidity and mortality. We present a 4-year-old male with KFS, who notably was also diagnosed with large-cell anaplastic medulloblastoma. We evaluated the genetic basis of co-occurring KFS and medulloblastoma and the role of MYO18B as related to medulloblastoma. Constitutional and somatic variant and copy number analyses were performed from DNA-based exome studies, along with RNA-sequencing of tumor tissue, to elucidate the genetic etiology of the co-existing disease states. We identified novel constitiitional compound heterozygous frameshift variants (NM_032608.5: p.Leu2257SerfsTerl6 and p.Arg2220SerfsTer74) each encoding a premature stop of translation in MYO18B, consistent with a diagnosis of KFS. We did not identify any somatic variants of known relevance or disease-relevant therapeutic targets in the tumor. The somatic copy number profile was suggestive of Group 3γ medulloblastoma. Relative to pediatric brain tremors, medulloblastoma, particularly, Group 3, had increased gene expression of MYO18B. In summary, coexisting constitutional and somatic diagnoses in this patient enabled the elucidation of the genetic etiology of KFS and provided support for the role of MYO18B in tumor suppression.

We present Samovar, a mosaic single-nucleotide variant (SNV) caller for linked-read whole-genome shotgun sequencing data. Samovar scores candidate sites using a random forest model trained using the input dataset that considers read quality, phasing, and linked-read characteristics. We show Samovar calls mosaic SNVs within a single sample with accuracy comparable to what previously required trios or matched tumor/normal pairs and outperform single-sample mosaic variant callers at MAF 5%-50% with at least 30x coverage. Furthermore, we use Samovar to find somatic variants in whole genome sequencing of both tumor and normal 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.

Proteoglycans have a core polypeptide connected to glycosaminoglycans (GAGs) via a common tetrasaccharide linker region. Defects in enzymes that synthesize the linker result in a group of autosomal recessive conditions called “linkeropathies”. Disease manifests with skeletal and connective tissue features, including short stature, hyperextensible skin, and joint hypermobility. We report a family with three affected pregnancies showing short limbs, cystic hygroma, and perinatal death. Two spontaneously aborted; one survived 1 day after term delivery, and had short limbs, bell-shaped thorax, 11 ribs, absent thumbs, and cleft palate. Exome sequencing of the proband and one affected fetus identified compound heterozygous missense variants, NM_007255.3: c.808C>T (p.(Arg270Cys)) and NM_007255.3: c.398A>G (p.(Gln133Arg)), in B4GALT7, a gene required for GAG linker biosynthesis. Homozygosity for p.(Arg270Cys), associated with partial loss of B4GALT7 function, causes Larsen of Reunion Island syndrome (LRS), however no previous studies have linked p.(Gln133Arg) to disease. The p.(Gln133Arg) and p.(Arg270Cys) variants were transfected into CHO pgsB-618 cells. High protein expression of p.(Gln133Arg) was found, with mislocalization, compared to p.(Arg270Cys) that had a normal Golgi-like pattern. The p.(Gln133Arg) had almost no enzyme activity and little production of heparan sulfate GAGs, while p.(Arg270Cys) only had 17% of wild-type activity. These findings expand the phenotype of B4GALT7-related linkeropathies to include lethal skeletal dysplasia due to more severe loss of function.

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.

Two sisters (ages 16 yr and 15 yr) have been followed by our clinical genetics team for several years. Both girls have severe intellectual disability, hypotonia, seizures, and distinctive craniofacial features. The parents are healthy and have no other children. Oligo array, fragile X testing, and numerous single-gene tests were negative. All four family members underwent research exome sequencing, which revealed a heterozygous nonsense mutation in ASXL3 (p.R1036X) that segregated with disease. Exome data and independent Sanger sequencing confirmed that the variant is de novo, suggesting possible germline mosaicism in one parent. The p.R1036X variant has never been observed in healthy human populations and has been previously reported as a pathogenic mutation. Truncating de novo mutations in ASXL3 cause Bainbridge–Ropers syndrome (BRPS), a developmental disorder with similarities to Bohring–Opitz syndrome. Fewer than 30 BRPS patients have been described in the literature; to our knowledge, this is the first report of the disorder in two related individuals. Our findings lend further support to intellectual disability, absent speech, autistic traits, hypotonia, and distinctive facial appearance as common emerging features of Bainbridge–Ropers syndrome.

Background We report a female patient with endocrine abnormalities, hypogonadotropic hypogonadism and amazia (breasts aplasia/hypoplasia but normal nipples and areolas) in a rare syndrome: Van Maldergem syndrome (VMS). Case presentation Our patient was first evaluated at age 4 for intellectual disability, craniofacial features, and auditory malformations. At age 15, she presented with no breast development and other findings consistent with hypogonadotropic hypogonadism. At age 37, she underwent whole exome sequencing (WES) to identify pathogenic variants. WES revealed compound heterozygous variants in DCHS1 (rs145099391:G > A, p.P197L & rs753548138:G > A, p.T2334 M) [RefSeq NM_003737.3], diagnostic of Van Maldergem syndrome (VMS-1). VMS is a rare autosomal disorder reported in only 13 patients, characterized by intellectual disability, typical craniofacial features, auditory malformations, hearing loss, skeletal and limb malformations, brain abnormalities with periventricular neuronal heterotopia and other variable anomalies. Our patient had similar phenotypic abnormalities. She also had hypogonadotropic hypogonadism and amazia. Based on the clinical findings reported, two previously published patients with VMS may also have been affected by hypogonadotropic hypogonadism, but endocrine abnormalities were not evaluated or mentioned. Conclusion This case highlights an individual with VMS, characterized by compound heterozygous variants in DCHS1. Our observations may provide additional information on the phenotypic spectrum of VMS, including hypogonadotropic hypogonadism and amazia. However, the molecular genetic basis for endocrine anomalies observed in some VMS patients, including ours, remains unexplained. Electronic supplementary material The online version of this article (10.1186/s13633-017-0052-z) contains supplementary material, which is available to authorized users.

Medulloblastoma is a pediatric embryonal tumor that can be classified into four molecular subgroups, each derived from a different progenitor cell. It is estimated that about 30-40% of patients will relapse, typically with recurrence at the primary site and of the same molecular subgroup. We present paired tumor/normal genomic analysis of an 18 year-old male who presented with non-Wnt/non-SHH medulloblastoma at age 12 and relapsed with metastatic disease of the falx cerebri 3 years later. Combination surgery, chemotherapy, and radiation were used in treatment of the primary and recurrent tumor. At a timepoint 6 years from original diagnosis, the patient presented with a cerebellar tumor histologically described as “consistent with recurrent medulloblastoma” with comment recommending genomics to confirm. The diagnosis was made based on near identifical morphology and retention of Neu-N and Synaptophysin in the tumor (confirmed by subsequent genetic analysis). The primary tumor and the tumor occurring 6 years after the primary diagnosis were analyzed by whole exome sequencing (blood and tumor tissues) to assess for germline variants, somatic mutation, and copy number variation. We observed no pathogenic germline variants in cancer predisposition genes. The tumor mutational profiles were distinct, with only 6 (1.8%) shared somatic variants between tumors. Specimen provenance was verified by germline variation and SRY coverage. Two targetable mutations within the RAS-MAPK pathway (PTPN11 p.Glu76Lys and PIK3CA p.Gly1007Arg) were present only in the new CNS tumor. Although the primary tumor harbored isochromosome 17q and a gain of chromosome 4, these somatic chromosomal aberrations were not detected in the new CNS tumor. RNA-seq was performed on both tumors and compared to pediatric CNS tumors from the University of California Santa Cruz Treehouse Initiative (n=434). The primary tumor clustered with the medulloblastoma patients by principal component analysis while the new CNS tumor clustered with a group of gliomas and non-medulloblastoma embryonal tumors. The primary tumor displayed evidence of overexpression of Group 4 medulloblastoma genes (e.g. EOMES, RBM24, SNCAIP, and UNC5D). These genes were not overexpressed in the new CNS tumor. Enrichment of genes commonly found in gliomas (e.g. BCAN, CHI3L2, PDGFRA, and SOX2) were noted in the new CNS tumor only. In summary, tumor genomic profiling of a primary medulloblastoma and the new CNS tumor arising 6 years later revealed two distinct sets of somatic mutations suggestive of second malignancy rather than recurrence in this patient. While second malignancy in the setting of medulloblastoma is a rare event, it has been documented, both in a time period consistent with that described in our patient and in the form of glioma. Thus, tumor profiling refined diagnosis in this patient allowing for a more accurate assessment of treatment and management options. Citation Format: Kathleen M. Schieffer, Katherine E. Miller, Daniel R. Boue, Daniel C. Koboldt, Patrick Brennan, Benjamin J. Kelly, Gregory Wheeler, Vincent Magrini, Amy Wetzel, Elizabeth Varga, Devon Dishman, Kristen Leraas, Vibhuti Agarwal, Mohamed S. AbdelBaki, Jonathan L. Finlay, Jeffrey R. Leonard, Peter White, Julie M. Gastier-Foster, Catherine E. Cottrell, Elaine R. Mardis, Richard K. Wilson. Molecular profiling identifies a second malignancy in a patient with medulloblastoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 484.

Two sisters (ages 16 yr and 15 yr) have been followed by our clinical genetics team for several years. Both girls have severe intellectual disability, hypotonia, seizures, and distinctive craniofacial features. The parents are healthy and have no other children. Oligo array, fragile X testing, and numerous single-gene tests were negative. All four family members underwent research exome sequencing, which revealed a heterozygous nonsense mutation in ASXL3 (p.R1036X) that segregated with disease. Exome data and independent Sanger sequencing confirmed that the variant is de novo, suggesting possible germline mosaicism in one parent. The p.R1036X variant has never been observed in healthy human populations and has been previously reported as a pathogenic mutation. Truncating de novo mutations in ASXL3 cause Bainbridge–Ropers syndrome (BRPS), a developmental disorder with similarities to Bohring–Opitz syndrome. Fewer than 30 BRPS patients have been described in the literature; to our knowledge, this is the first report of the disorder in two related individuals. Our findings lend further support to intellectual disability, absent speech, autistic traits, hypotonia, and distinctive facial appearance as common emerging features of Bainbridge–Ropers syndrome.

While advances in genome sequencing technology make population-scale genomics a possibility, current approaches for analysis of these data rely upon parallelization strategies that have limited scalability, complex implementation and lack reproducibility. Churchill, a balanced regional parallelization strategy, overcomes these challenges, fully automating the multiple steps required to go from raw sequencing reads to variant discovery. Through implementation of novel deterministic parallelization techniques, Churchill allows computationally efficient analysis of a high-depth whole genome sample in less than two hours. The method is highly scalable, enabling full analysis of the 1000 Genomes raw sequence dataset in a week using cloud resources. http://churchill.nchri.org/. Electronic supplementary material The online version of this article (doi:10.1186/s13059-014-0577-x) contains supplementary material, which is available to authorized users.

Erdheim-Chester disease (ECD) is a rare form of non-Langerhans cell histiocytosis. Herein we describe documented skeletal and pericardial involvement by ECD producing cardiac tamponade in a 30-year-old woman. The diagnosis of ECD was established by histopathology, immunocytochemistry, and by radiologic studies demonstrating diffuse, bilateral, symmetrical osteosclerosis of the long bones, sparing the epiphyses and axial skeleton. Scintography using methyl diphosphonate showed increased uptake in involved bone. The patient presented with jaundice and hepatic congestion produced by cardiac tamponade. Pericardial biopsy revealed xanthogranulomatous lesions comprised of foamy and lipid-laden macrophages, multinucleated giant cells, monocytes, and lymphocytes in a mesh of fibrosis. Immunohistochemical staining was positive for CD68 and negative for CD1a, consistent with ECD rather than with the much more common Langerhans cell form of histiocytosis.

Aplasia cutis congenita (ACC) is the absence of localized or widespread areas of skin at birth. A frequently cited classification schema is based on location and the presence of associated anomalies. Histologically it is characterized by dermal fibrosis and absence of adenexal structures. We present a newborn female with extensive truncal ACC associated with fetus papyraceus.

The transcription factor AmrZ regulates genes important for P. aeruginosa virulence, including type IV pili, extracellular polysaccharides, and the flagellum; however, the global effect of AmrZ on gene expression remains unknown, and therefore, AmrZ may directly regulate many additional genes that are crucial for infection. Compared to the wild type strain, a ΔamrZ mutant exhibits a rugose colony phenotype, which is commonly observed in variants that accumulate the intracellular second messenger cyclic diguanylate (c-di-GMP). Cyclic di-GMP is produced by diguanylate cyclases (DGC) and degraded by phosphodiesterases (PDE). We hypothesized that AmrZ limits the intracellular accumulation of c-di-GMP through transcriptional repression of gene(s) encoding a DGC. In support of this, we observed elevated c-di-GMP in the ΔamrZ mutant compared to the wild type strain. Consistent with other strains that accumulate c-di-GMP, when grown as a biofilm, the ΔamrZ mutant formed larger microcolonies than the wild-type strain. This enhanced biofilm formation was abrogated by expression of a PDE. To identify potential target DGCs, a ChIP-Seq was performed and identified regions of the genome that are bound by AmrZ. RNA-Seq experiments revealed the entire AmrZ regulon, and characterized AmrZ as an activator or repressor at each binding site. We identified an AmrZ-repressed DGC-encoding gene (PA4843) from this cohort, which we named AmrZ dependent cyclase A (adcA). PAO1 overexpressing adcA accumulates 29-fold more c-di-GMP than the wild type strain, confirming the cyclase activity of AdcA. In biofilm reactors, a ΔamrZ ΔadcA double mutant formed smaller microcolonies than the single ΔamrZ mutant, indicating adcA is responsible for the hyper biofilm phenotype of the ΔamrZ mutant. This study combined the techniques of ChIP-Seq and RNA-Seq to define the comprehensive regulon of a bifunctional transcriptional regulator. Moreover, we identified a c-di-GMP mediated mechanism for AmrZ regulation of biofilm formation and chronicity., Author Summary Pathogenic bacteria such as Pseudomonas aeruginosa utilize a wide variety of systems to sense and respond to the changing conditions during an infection. When a stress is sensed, signals are transmitted to impact expression of many genes that allow the bacterium to adapt to the changing conditions. AmrZ is a protein that regulates production of several virulence-associated gene products, though we predicted that its role in virulence was more expansive than previously described. Transcription factors such as AmrZ often affect the expression of a gene by binding and promoting or inhibiting expression of the target gene. Two global techniques were utilized to determine where AmrZ binds in the genome, and what effect AmrZ has once bound. This approach revealed that AmrZ represses the production of a signaling molecule called cyclic diguanylate, which is known to induce the formation of difficult to treat communities of bacteria called biofilms. This study also identified many novel targets of AmrZ to promote future studies of this regulator. Collectively, these data can be utilized to develop treatments to inhibit biofilm formation during devastating chronic infections.

Nuclear magnetic resonance (NMR) spectroscopy is a non-invasive method of acquiring a metabolic profile from biofluids. This metabolic information may provide keys to the early detection of exposure to a toxin. A typical NMR toxicology data set has low sample size and high dimensionality. Thus, traditional pattern recognition techniques are not always feasible. In this paper, we evaluate several common alternatives for isolating these biomarkers. The fold test, unpaired t-test, and paired t-test were performed on an NMR-derived toxicological data set and results were compared. The paired t-test method was preferred, due to its ability to attribute statistical significance, to take into consideration consistency of a single subject over a time course, and to mitigate the low sample, high dimensionality problem. We then grouped the resulting statistically salient potential biomarkers based on their significance patterns and compared results to several known metabolites affected by the tested toxin. Based on these results, we present a statistical protocol of sequential t-tests and clustering techniques for identifying putative biomarkers. We then present the results of this protocol applied to a specific real world toxicological data set.

A 75-year-old Caucasian male with papular elastolytic giant cell granuloma was treated successfully with hydroxychloroquine (200 mg p.o. bid for 4 months). Treatment had to be stopped because of early retinopathy. Treatment was started with quinacrine, which accomplished partial control.