KMT2A-D pathogenicity, prevalence, and variation according to a population database.

INTRODUCTION: The KMT2 family of genes is essential epigenetic regulators promoting gene expression. The gene family contains three subgroups, each with two paralogues: KMT2A and KMT2B; KMT2C and KMT2D; KMT2F and KMT2G. KMT2A-D are among the most frequent somatically altered genes in several different cancer types. Somatic KMT2A rearrangements are well-characterized in infant leukemia (IL), and growing evidence supports the role of additional family members (KMT2B, KMT2C, and KMT2D) in leukemogenesis. Enrichment of rare heterozygous frameshift variants in KMT2A and C has been reported in acute myeloid leukemia (AML), IL, and solid tumors. Currently, the non-synonymous variation, prevalence, and penetrance of these four genes are unknown. METHODS: This study determined the prevalence of pathogenic/likely pathogenic (P/LP) germline KMT2A-D variants in a cancer-free adult population from the Genome Aggregation Database (gnomAD). Two methods of variant interpretation were utilized: a manual genomic variant interpretation and an automated ACMG pipeline. RESULTS: The ACMG pipeline identified considerably fewer P/LP variants (n = 89) compared to the manual method (n = 660) in all 4 genes. Consequently, the total P/LP prevalence and allele frequency (AF) were higher in the manual method (1:112, AF = 4.46E-03) than in ACMG (1:832, AF = 6.01E-04). Multiple ancestry-exclusive P/LP variants were identified along with an increased frequency in males compared to females. Many of these variants identified in this population database are also associated with severe juvenile conditions. CONCLUSION: These data demonstrate that putatively functional germline variation in these developmentally important genes is more common than previously appreciated and identification in cancer-free adults may indicate incomplete penetrance for many of these variants. Future research should examine a genetic predisposing role in IL and other pediatric cancers.

Transcriptional analysis of muscle tissue and isolated satellite cells in spastic cerebral palsy.

AIM: To analyze transcriptomes from muscle tissue and cells to improve our understanding of differences in gene expression and molecular function in cerebral palsy (CP) muscle. METHOD: In this case-control study, eight participants with CP (five males, three females; mean [SD] age 14y 2mo [1y 8mo]) and 11 comparison individuals (eight males, three females; mean [SD] age 14y 0mo [2y 6mo]) were enrolled after informed consent/assent and skeletal muscle was obtained during surgery. RNA was extracted from tissue and from primary satellite cells grown to form myotubes in vitro. RNA sequencing data were analyzed using validated informatics pipelines. RESULTS: Analysis identified expression of 6308 genes in the tissue samples and 7459 in the cultured cells. Significant differential expression between CP and control was identified in 87 genes in the tissue and 90 genes in isolated satellite cell-derived myotube cultures. INTERPRETATION: Both tissue and cell analyses identified differential expression of genes associated with muscle development and multiple pathways of interest. What this paper adds Expression differences were found in muscle tissue and in isolated muscle cells. There was low variability in expression among cells isolated from different muscles. Expression differences suggest complex functional alterations in spastic cerebral palsy.

Germline Sequencing Identifies Rare Variants in Finnish Subjects with Familial Germ Cell Tumors.

PURPOSE: Pediatric germ cell tumors are rare, representing about 3% of childhood malignancies in children less than 15 years of age, presenting in neonates or adolescents with a greater incidence noted in older adolescents. Aberrations in primordial germ cell proliferation/differentiation can lead to a variety of neoplasms, including teratomas, embryonal carcinoma, choriocarcinoma, and yolk sac tumors. PATIENTS AND METHODS: Three Finnish families with varying familial germ cell tumors were identified, and whole-genome sequencing was performed using an Illumina sequencing platform. In total, 22 unique subjects across the three families were sequenced. Family 1 proband (female) was affected by malignant ovarian teratoma, Family 2 proband (female) was affected by sacrococcygeal teratoma with yolk sac tumor in the setting of Cornelia de Lange syndrome, and Family 3 proband (male) was affected by malignant testicular teratoma. Rare variants were identified using an autosomal recessive or de novo model of inheritance. RESULTS: For family 1 proband (female), an autosomal recessive or de novo model of inheritance identified variants of interest in the following genes: CD109, IKBKB, and CTNNA3, SUPT6H, MUC5AC, and FRG1. Family 2 proband (female) analysis identified gene variants of interest in the following genes: LONRF2, ANO7, HS6ST1, PRB2, and DNM2. Family 3 proband (male) analysis identified the following potential genes: CRIPAK, KRTAP5-7, and CACNA1B. CONCLUSION: Leveraging deep pedigrees and next-generation sequencing, rare germline variants were identified that were enriched in three families from Finland with a history of familial germ cell tumors. The data presented support the importance of germline mutations when analyzing complex cancers with a low somatic mutation landscape.

Strong concordance between RNA structural and single nucleotide variants identified via next generation sequencing techniques in primary pediatric leukemia and patient-derived xenograft samples.

Acute leukemia represents the most common pediatric malignancy comprising diverse subtypes with varying prognosis and treatment outcomes. New and targeted treatment options are warranted for this disease. Patient-derived xenograft (PDX) models are increasingly being used for preclinical testing of novel treatment modalities. A novel approach involving targeted error-corrected RNA sequencing using ArcherDX HemeV2 kit was employed to compare 25 primary pediatric acute leukemia samples and their corresponding PDX samples. A comparison of the primary samples and PDX samples revealed a high concordance between single nucleotide variants and gene fusions whereas other complex structural variants were not as consistent. The presence of gene fusions representing the major driver mutations at similar allelic frequencies in PDX samples compared to primary samples and over multiple passages confirms the utility of PDX models for preclinical drug testing. Characterization and tracking of these novel cryptic fusions and exonal variants in PDX models is critical in assessing response to potential new therapies.

Error-corrected sequencing strategies enable comprehensive detection of leukemic mutations relevant for diagnosis and minimal residual disease monitoring.

BACKGROUND: Pediatric leukemias have a diverse genomic landscape associated with complex structural variants, including gene fusions, insertions and deletions, and single nucleotide variants. Routine karyotype and fluorescence in situ hybridization (FISH) techniques lack sensitivity for smaller genomic alternations. Next-generation sequencing (NGS) assays are being increasingly utilized for assessment of these various lesions. However, standard NGS lacks quantitative sensitivity for minimal residual disease (MRD) surveillance due to an inherently high error rate. METHODS: Primary bone marrow samples from pediatric leukemia (n = 32) and adult leukemia subjects (n = 5), cell line MV4-11, and an umbilical cord sample were utilized for this study. Samples were sequenced using molecular barcoding with targeted DNA and RNA library enrichment techniques based on anchored multiplexed PCR (AMP®) technology, amplicon based error-corrected sequencing (ECS) or a human cancer transcriptome assay. Computational analyses were performed to quantitatively assess limit of detection (LOD) for various DNA and RNA lesions, which could be systematically used for MRD assays. RESULTS: Matched leukemia patient samples were analyzed at three time points; diagnosis, end of induction (EOI), and relapse. Similar to flow cytometry for ALL MRD, the LOD for point mutations by these sequencing strategies was ≥0.001. For DNA structural variants, FLT3 internal tandem duplication (ITD) positive cell line and patient samples showed a LOD of ≥0.001 in addition to previously unknown copy number losses in leukemia genes. ECS in RNA identified multiple novel gene fusions, including a SPANT-ABL gene fusion in an ALL patient, which could have been used to alter therapy. Collectively, ECS for RNA demonstrated a quantitative and complex landscape of RNA molecules with 12% of the molecules representing gene fusions, 12% exon duplications, 8% exon deletions, and 68% with retained introns. Droplet digital PCR validation of ECS-RNA confirmed results to single mRNA molecule quantities. CONCLUSIONS: Collectively, these assays enable a highly sensitive, comprehensive, and simultaneous analysis of various clonal leukemic mutations, which can be tracked across disease states (diagnosis, EOI, and relapse) with a high degree of sensitivity. The approaches and results presented here highlight the ability to use NGS for MRD tracking.

Mapping Citrullinated Sites in Multiple Organs of Mice Using Hypercitrullinated Library.

Protein citrullination (or deimination), an irreversible post-translational modification, has been implicated in several physiological and pathological processes, including gene expression regulation, apoptosis, rheumatoid arthritis, and Alzheimer's disease. Several research studies have been carried out on citrullination under many conditions. However, until now, challenges in sample preparation and data analysis have made it difficult to confidently identify a citrullinated protein and assign the citrullinated site. To overcome these limitations, we generated a mouse hyper-citrullinated spectral library and set up coordinates to confidently identify and validate citrullinated sites. Using this workflow, we detect a four-fold increase in citrullinated proteome coverage across six mouse organs compared with the current state-of-the art techniques. Our data reveal that the subcellular distribution of citrullinated proteins is tissue-type-dependent and that citrullinated targets are involved in fundamental physiological processes, including the metabolic process. These data represent the first report of a hyper-citrullinated library for the mouse and serve as a central resource for exploring the role of citrullination in this organism.

Subphenotype meta-analysis of testicular cancer genome-wide association study data suggests a role for RBFOX family genes in cryptorchidism susceptibility.

STUDY QUESTION: Can subphenotype analysis of genome-wide association study (GWAS) data from subjects with testicular germ cell tumor (TGCT) provide insight into cryptorchidism (undescended testis, UDT) susceptibility? SUMMARY ANSWER: Suggestive intragenic GWAS signals common to UDT, TGCT case-case and TGCT case-control analyses occur in genes encoding RBFOX RNA-binding proteins (RBPs) and their neurodevelopmental targets. WHAT IS KNOWN ALREADY: UDT is a strong risk factor for TGCT, but while genetic risk factors for TGCT are well-known, genetic susceptibility to UDT is poorly understood and appears to be more complex. STUDY DESIGN, SIZE, DURATION: We performed a secondary subphenotype analysis of existing GWAS data from the Testicular Cancer Consortium (TECAC) and compared these results with our previously published UDT GWAS data, and with data previously acquired from studies of the fetal rat gubernaculum. PARTICIPANTS/MATERIALS, SETTING, METHODS: Studies from the National Cancer Institute (NCI), United Kingdom (UK) and University of Pennsylvania (Penn) that enrolled white subjects were the source of the TGCT GWAS data. We completed UDT subphenotype case-case (TGCT/UDT vs TGCT/non-UDT) and case-control (TGCT/UDT vs control), collectively referred to as 'TECAC' analyses, followed by a meta-analysis comprising 129 TGCT/UDT cases, 1771 TGCT/non-UDT cases, and 3967 unaffected controls. We reanalyzed our UDT GWAS results comprising 844 cases and 2718 controls by mapping suggestive UDT and TECAC signals (defined as P < 0.001) to genes using Ingenuity Pathway Analysis (IPA®). We compared associated pathways and enriched gene categories common to all analyses after Benjamini-Hochberg multiple testing correction, and analyzed transcript levels and protein expression using qRT-PCR and rat fetal gubernaculum confocal imaging, respectively. MAIN RESULTS AND THE ROLE OF CHANCE: We found suggestive signals within 19 genes common to all three analyses, including RBFOX1 and RBFOX3, neurodevelopmental paralogs that encode RBPs targeting (U)GCATG-containing transcripts. Ten of the 19 genes participate in neurodevelopment and/or contribute to risk of neurodevelopmental disorders. Experimentally predicted RBFOX gene targets were strongly overrepresented among suggestive intragenic signals for the UDT (117 of 628 (19%), P = 3.5 × 10-24), TECAC case-case (129 of 711 (18%), P = 2.5 × 10-27) and TECAC case-control (117 of 679 (17%), P = 2 × 10-21) analyses, and a majority of the genes common to all three analyses (12 of 19 (63%), P = 3 × 10-9) are predicted RBFOX targets. Rbfox1, Rbfox2 and their encoded proteins are expressed in the rat fetal gubernaculum. Predicted RBFOX targets are also enriched among transcripts differentially regulated in the fetal gubernaculum during normal development (P = 3 × 10-31), in response to in vitro hormonal stimulation (P = 5 × 10-45) and in the cryptorchid LE/orl rat (P = 2 × 10-42). LARGE SCALE DATA: GWAS data included in this study are available in the database of Genotypes and Phenotypes (dbGaP accession numbers phs000986.v1.p1 and phs001349.v1p1). LIMITATIONS, REASONS FOR CAUTION: These GWAS data did not reach genome-wide significance for any individual analysis. UDT appears to have a complex etiology that also includes environmental factors, and such complexity may require much larger sample sizes than are currently available. The current methodology may also introduce bias that favors false discovery of larger genes. WIDER IMPLICATIONS OF THE FINDINGS: Common suggestive intragenic GWAS signals suggest that RBFOX paralogs and other neurodevelopmental genes are potential UDT risk candidates, and potential TGCT susceptibility modifiers. Enrichment of predicted RBFOX targets among differentially expressed transcripts in the fetal gubernaculum additionally suggests a role for this RBP family in regulation of testicular descent. As RBFOX proteins regulate alternative splicing of Calca to generate calcitonin gene-related peptide, a protein linked to development and function of the gubernaculum, additional studies that address the role of these proteins in UDT are warranted. STUDY FUNDING/COMPETING INTEREST(S): The Eunice Kennedy Shriver National Institute for Child Health and Human Development (R01HD060769); National Center for Research Resources (P20RR20173), National Institute of General Medical Sciences (P20GM103464), Nemours Biomedical Research, the Testicular Cancer Consortium (U01CA164947), the Intramural Research Program of the NCI, a support services contract HHSN26120130003C with IMS, Inc., the Abramson Cancer Center at Penn, National Cancer Institute (CA114478), the Institute of Cancer Research, UK and the Wellcome Trust Case-Control Consortium (WTCCC) 2. None of the authors reports a conflict of interest.

A survey of proteomic biomarkers for heterotopic ossification in blood serum.

BACKGROUND: Heterotopic ossification (HO) is a significant problem for wounded warriors surviving high-energy blast injuries; however, currently, there is no biomarker panel capable of globally characterizing, diagnosing, and monitoring HO progression. The aim of this study was to identify biomarkers for HO using proteomic techniques and blood serum. METHODS: Isobaric tags for relative and absolute quantitation (iTRAQ) was used to generate a semi-quantitative global proteomics survey of serum from patients with and without heterotopic ossification. Leveraging the iTRAQ data, a targeted selection reaction monitoring mass spectrometry (SRM-MS) assay was developed for 10 protein candidates: alkaline phosphatase, osteocalcin, alpha-2 type I collagen, collagen alpha-1(V) chain isoform 2 preprotein, bone sialoprotein 2, phosphatidate phosphatase LPIN2, osteomodulin, protein phosphatase 1J, and RRP12-like protein. RESULTS: The proteomic survey of serum from both healthy and disease patients includes 1220 proteins and was enriched for proteins involved in the response to elevated platelet Ca+2, wound healing, and extracellular matrix organization. Proteolytic peptides from three of the ten SRM-MS proteins, osteocalcin preprotein, osteomodulin precursor, and collagen alpha-1(v) chain isoform 2 preprotein from serum, are potential clinical biomarkers for HO. CONCLUSIONS: This study is the first reported SRM-MS analysis of serum from individuals with and without heterotopic ossification, and differences in the serum proteomic profile between healthy and diseased subjects were identified. Furthermore, our results indicate that normal wound healing signals can impact the ability to identify biomarkers, and a multi-protein panel assay, including osteocalcin preproprotein, osteomodulin precursor, and collagen alpha-1(v) chain isoform 2 preprotein, may provide a solution for HO detection and monitoring.

Mapping Biological Networks from Quantitative Data-Independent Acquisition Mass Spectrometry: Data to Knowledge Pipelines.

Data-independent acquisition mass spectrometry (DIA-MS) strategies and applications provide unique advantages for qualitative and quantitative proteome probing of a biological sample allowing constant sensitivity and reproducibility across large sample sets. These advantages in LC-MS/MS are being realized in fundamental research laboratories and for clinical research applications. However, the ability to translate high-throughput raw LC-MS/MS proteomic data into biological knowledge is a complex and difficult task requiring the use of many algorithms and tools for which there is no widely accepted standard and best practices are slowly being implemented. Today a single tool or approach inherently fails to capture the full interpretation that proteomics uniquely supplies, including the dynamics of quickly reversible chemically modified states of proteins, irreversible amino acid modifications, signaling truncation events, and, finally, determining the presence of protein from allele-specific transcripts. This chapter highlights key steps and publicly available algorithms required to translate DIA-MS data into knowledge.

Dual Labeling Biotin Switch Assay to Reduce Bias Derived From Different Cysteine Subpopulations: A Method to Maximize S-Nitrosylation Detection.

RATIONALE: S-nitrosylation (SNO), an oxidative post-translational modification of cysteine residues, responds to changes in the cardiac redox-environment. Classic biotin-switch assay and its derivatives are the most common methods used for detecting SNO. In this approach, the labile SNO group is selectively replaced with a single stable tag. To date, a variety of thiol-reactive tags have been introduced. However, these methods have not produced a consistent data set, which suggests an incomplete capture by a single tag and potentially the presence of different cysteine subpopulations. OBJECTIVE: To investigate potential labeling bias in the existing methods with a single tag to detect SNO, explore if there are distinct cysteine subpopulations, and then, develop a strategy to maximize the coverage of SNO proteome. METHODS AND RESULTS: We obtained SNO-modified cysteine data sets for wild-type and S-nitrosoglutathione reductase knockout mouse hearts (S-nitrosoglutathione reductase is a negative regulator of S-nitrosoglutathione production) and nitric oxide-induced human embryonic kidney cell using 2 labeling reagents: the cysteine-reactive pyridyldithiol and iodoacetyl based tandem mass tags. Comparison revealed that <30% of the SNO-modified residues were detected by both tags, whereas the remaining SNO sites were only labeled by 1 reagent. Characterization of the 2 distinct subpopulations of SNO residues indicated that pyridyldithiol reagent preferentially labels cysteine residues that are more basic and hydrophobic. On the basis of this observation, we proposed a parallel dual-labeling strategy followed by an optimized proteomics workflow. This enabled the profiling of 493 SNO sites in S-nitrosoglutathione reductase knockout hearts. CONCLUSIONS: Using a protocol comprising 2 tags for dual-labeling maximizes overall detection of SNO by reducing the previously unrecognized labeling bias derived from different cysteine subpopulations.

Development of Bioinformatics Pipeline for Analyzing Clinical Pediatric NGS Data.

Using an Illumina exome sequencing dataset generated from pediatric Acute Myeloid Leukemia patients (AML; type FLT3/ITD+) a comprehensive bioinformatics pipeline was developed to aid in a better clinical understanding of the genetic data associated with the clinical phenotype. The pipeline starts with raw next generation sequencing reads and using both publicly available resources and custom scripts, analyzes the genomic data for variants associated with pediatric AML. By incorporating functional information such as Gene Ontology annotation and protein-protein interactions, the methodology prioritizes genomic variants and returns disease specific results and knowledge maps. Furthermore, it compares the somatic mutations at diagnosis with the somatic mutations at relapse and outputs variants and functional annotations that are specific for the relapse state.

Citrullination of myofilament proteins in heart failure.

AIMS: Citrullination, the post-translational conversion of arginine to citrulline by the enzyme family of peptidylarginine deiminases (PADs), is associated with several diseases, and specific citrullinated proteins have been shown to alter function while others act as auto-antigens. In this study, we identified citrullinated proteins in human myocardial samples, from healthy and heart failure patients, and determined several potential functional consequences. Further we investigated PAD isoform cell-specific expression in the heart. METHODS AND RESULTS: A citrullination-targeted proteomic strategy using data-independent (SWATH) acquisition method was used to identify the modified cardiac proteins. Citrullinated-induced sarcomeric proteins were validated using two-dimensional gel electrophoresis and investigated using biochemical and functional assays. Myocardial PAD isoforms were confirmed by RT-PCR with PAD2 being the major isoform in myocytes. In total, 304 citrullinated sites were identified that map to 145 proteins among the three study groups: normal, ischaemia, and dilated cardiomyopathy. Citrullination of myosin (using HMM fragment) decreased its intrinsic ATPase activity and inhibited the acto-HMM-ATPase activity. Citrullinated TM resulted in stronger F-actin binding and inhibited the acto-HMM-ATPase activity. Citrullinated TnI did not alter the binding to F-actin or acto-HMM-ATPase activity. Overall, citrullination of sarcomeric proteins caused a decrease in Ca(2+) sensitivity in skinned cardiomyocytes, with no change in maximal calcium-activated force or hill coefficient. CONCLUSION: Citrullination unique to the cardiac proteome was identified. Our data indicate important structural and functional alterations to the cardiac sarcomere and the contribution of protein citrullination to this process.