Insights into Sex Chromosome Evolution and Aging from the Genome of a Short-Lived Fish.

The killifish Nothobranchius furzeri is the shortest-lived vertebrate that can be bred in the laboratory. Its rapid growth, early sexual maturation, fast aging, and arrested embryonic development (diapause) make it an attractive model organism in biomedical research. Here, we report a draft sequence of its genome that allowed us to uncover an intra-species Y chromosome polymorphism representing-in real time-different stages of sex chromosome formation that display features of early mammalian XY evolution "in action." Our data suggest that gdf6Y, encoding a TGF-ß family growth factor, is the master sex-determining gene in N. furzeri. Moreover, we observed genomic clustering of aging-related genes, identified genes under positive selection, and revealed significant similarities of gene expression profiles between diapause and aging, particularly for genes controlling cell cycle and translation. The annotated genome sequence is provided as an online resource (http://www.nothobranchius.info/NFINgb).

Author Correction: Epigenetic stress responses induce muscle stem-cell ageing by Hoxa9 developmental signals.

An Amendment to this paper has been published and can be accessed via a link at the top of the paper.

The HLA ligandome of oropharyngeal squamous cell carcinomas reveals shared tumour-exclusive peptides for semi-personalised vaccination.

BACKGROUND: The immune peptidome of OPSCC has not previously been studied. Cancer-antigen specific vaccination may improve clinical outcome and efficacy of immune checkpoint inhibitors such as PD1/PD-L1 antibodies. METHODS: Mapping of the OPSCC HLA ligandome was performed by mass spectrometry (MS) based analysis of naturally presented HLA ligands isolated from tumour tissue samples (n = 40) using immunoaffinity purification. The cohort included 22 HPV-positive (primarily HPV-16) and 18 HPV-negative samples. A benign reference dataset comprised of the HLA ligandomes of benign haematological and tissue datasets was used to identify tumour-associated antigens. RESULTS: MS analysis led to the identification of naturally HLA-presented peptides in OPSCC tumour tissue. In total, 22,769 peptides from 9485 source proteins were detected on HLA class I. For HLA class II, 15,203 peptides from 4634 source proteins were discovered. By comparative profiling against the benign HLA ligandomic datasets, 29 OPSCC-associated HLA class I ligands covering 11 different HLA allotypes and nine HLA class II ligands were selected to create a peptide warehouse. CONCLUSION: Tumour-associated peptides are HLA-presented on the cell surfaces of OPSCCs. The established warehouse of OPSCC-associated peptides can be used for downstream immunogenicity testing and peptide-based immunotherapy in (semi)personalised strategies.

Analysis, identification and visualization of subgroups in genomics.

MOTIVATION: Cancer is a complex and heterogeneous disease involving multiple somatic mutations that accumulate during its progression. In the past years, the wide availability of genomic data from patients' samples opened new perspectives in the analysis of gene mutations and alterations. Hence, visualizing and further identifying genes mutated in massive sets of patients are nowadays a critical task that sheds light on more personalized intervention approaches. RESULTS: Here, we extensively review existing tools for visualization and analysis of alteration data. We compare different approaches to study mutual exclusivity and sample coverage in large-scale omics data. We complement our review with the standalone software AVAtar ('analysis and visualization of alteration data') that integrates diverse aspects known from different tools into a comprehensive platform. AVAtar supplements customizable alteration plots by a multi-objective evolutionary algorithm for subset identification and provides an innovative and user-friendly interface for the evaluation of concurrent solutions. A use case from personalized medicine demonstrates its unique features showing an application on vaccination target selection. AVAILABILITY: AVAtar is available at: https://github.com/sysbio-bioinf/avatar. CONTACT: hans.kestler@uni-ulm.de, phone: +49 (0) 731 500 24 500, fax: +49 (0) 731 500 24 502.

Epigenetic stress responses induce muscle stem-cell ageing by Hoxa9 developmental signals.

The functionality of stem cells declines during ageing, and this decline contributes to ageing-associated impairments in tissue regeneration and function. Alterations in developmental pathways have been associated with declines in stem-cell function during ageing, but the nature of this process remains poorly understood. Hox genes are key regulators of stem cells and tissue patterning during embryogenesis with an unknown role in ageing. Here we show that the epigenetic stress response in muscle stem cells (also known as satellite cells) differs between aged and young mice. The alteration includes aberrant global and site-specific induction of active chromatin marks in activated satellite cells from aged mice, resulting in the specific induction of Hoxa9 but not other Hox genes. Hoxa9 in turn activates several developmental pathways and represents a decisive factor that separates satellite cell gene expression in aged mice from that in young mice. The activated pathways include most of the currently known inhibitors of satellite cell function in ageing muscle, including Wnt, TGFß, JAK/STAT and senescence signalling. Inhibition of aberrant chromatin activation or deletion of Hoxa9 improves satellite cell function and muscle regeneration in aged mice, whereas overexpression of Hoxa9 mimics ageing-associated defects in satellite cells from young mice, which can be rescued by the inhibition of Hoxa9-targeted developmental pathways. Together, these data delineate an altered epigenetic stress response in activated satellite cells from aged mice, which limits satellite cell function and muscle regeneration by Hoxa9-dependent activation of developmental pathways.

Synergistic targeting and resistance to PARP inhibition in DNA damage repair-deficient pancreatic cancer.

OBJECTIVE: ATM serine/threonine kinase (ATM) is the most frequently mutated DNA damage response gene, involved in homologous recombination (HR), in pancreatic ductal adenocarcinoma (PDAC). DESIGN: Combinational synergy screening was performed to endeavour a genotype-tailored targeted therapy. RESULTS: Synergy was found on inhibition of PARP, ATR and DNA-PKcs (PAD) leading to synthetic lethality in ATM-deficient murine and human PDAC. Mechanistically, PAD-induced PARP trapping, replication fork stalling and mitosis defects leading to P53-mediated apoptosis. Most importantly, chemical inhibition of ATM sensitises human PDAC cells toward PAD with long-term tumour control in vivo. Finally, we anticipated and elucidated PARP inhibitor resistance within the ATM-null background via whole exome sequencing. Arising cells were aneuploid, underwent epithelial-mesenchymal-transition and acquired multidrug resistance (MDR) due to upregulation of drug transporters and a bypass within the DNA repair machinery. These functional observations were mirrored in copy number variations affecting a region on chromosome 5 comprising several of the upregulated MDR genes. Using these findings, we ultimately propose alternative strategies to overcome the resistance. CONCLUSION: Analysis of the molecular susceptibilities triggered by ATM deficiency in PDAC allow elaboration of an efficient mutation-specific combinational therapeutic approach that can be also implemented in a genotype-independent manner by ATM inhibition.

Prospective longitudinal study of immune checkpoint molecule (ICM) expression in immune cell subsets during curative conventional therapy of head and neck squamous cell carcinoma (HNSCC).

Programmed-death-1 (PD1) antibodies are approved for recurrent and metastatic head and neck squamous cell carcinoma. Multiple drugs targeting costimulatory and coinhibitory immune checkpoint molecules (ICM) have been discovered. However, it remains unknown how these ICM are affected by curative conventional therapy on different immune cell subsets during the course of treatment. In the prospective noninterventional clinical study titled "Immune Response Evaluation to Curative conventional Therapy" (NCT03053661), 22 patients were prospectively enrolled. Blood samples were drawn at defined time points throughout curative conventional treatment and follow-up. Immune cells (IC) from the different time points were assessed by multicolor flow cytometry. The following ICM were measured by flow cytometry: PD1, CTLA4, BTLA, CD137, CD27, GITR, OX40, LAG3 and TIM3. Dynamics of ICM expression were assessed using nonparametric paired samples tests. Significant changes were noted for PD1, BTLA and CD27 on multiple IC types during or after radiotherapy. Nonsignificant trends for increased expression of OX40 and GITR from baseline until the end of RT were observed on CD4 T cells and CD4+ CD39+ T cells. In patients with samples at recurrence of disease, a nonsignificant increase of TIM3 and LAG3 positive CD4+ CD39+ T cells was evident, accompanied by an increase of double positive cells for TIM3/LAG3. Potential future targets to be combined with RT in the conventional treatment and anti-PD1/PD-L could be BTLA agonists, or agonistic antibodies to costimulatory ICM like CD137, OX40 or GITR. The combination of cetuximab with CD27 agonistic antibodies enhancing ADCC or the targeting of TIM3/LAG3 may be another promising strategy.

Protein Kinase D1, Reduced in Human Pancreatic Tumors, Increases Secretion of Small Extracellular Vesicles From Cancer Cells That Promote Metastasis to Lung in Mice.

BACKGROUND & AIMS: Pancreatic tumor cells release small extracellular vesicles (sEVs, exosomes) that contain lipids and proteins, RNA, and DNA molecules that might promote formation of metastases. It is not clear what cargo these vesicles contain and how they are released. Protein kinase D1 (PRKD1) inhibits cell motility and is believed to be dysregulated in pancreatic ductal adenocarcinomas. We investigated whether it regulates production of sEVs in pancreatic cancer cells and their ability to form premetastatic niches for pancreatic cancer cells in mice. METHODS: We analyzed data from UALCAN and human pancreatic tissue microarrays to compare levels of PRKD1 between tumor and nontumor tissues. We studied mice with pancreas-specific disruption of Prkd1 (PRKD1KO mice), mice that express oncogenic KRAS (KC mice), and KC mice with disruption of Prkd1 (PRKD1KO-KC mice). Subcutaneous xenograft tumors were grown in NSG mice from Panc1 cells; some mice were then given injections of sEVs. Pancreata and lung tissues from mice were analyzed by histology, immunohistochemistry, and/or quantitative polymerase chain reaction; we performed nanoparticle tracking analysis of plasma sEVs. The Prkd1 gene was disrupted in Panc1 cells using CRISPR-Cas9 or knocked down with small hairpin RNAs, or PRKD1 activity was inhibited with the selective inhibitor CRT0066101. Pancreatic cancer cell lines were analyzed by gene-expression microarray, quantitative polymerase chain reaction, immunoblot, and immunofluorescence analyses. sEVs secreted by Panc1 cell lines were analyzed by flow cytometry, transmission electron microscopy, and mass spectrometry. RESULTS: Levels of PRKD1 were reduced in human pancreatic ductal adenocarcinoma tissues compared with nontumor tissues. PRKD1KO-KC mice developed more pancreatic intraepithelial neoplasia, at a faster rate, than KC mice, and had more lung metastases and significantly shorter average survival time. Serum from PRKD1KO-KC mice had increased levels of sEVs compared with KC mice. Pancreatic cancer cells with loss or inhibition of PRKD1 increased secretion of sEVs; loss of PRKD1 reduced phosphorylation of its substrate, cortactin, resulting in increased F-actin levels at the plasma membrane. sEVs from cells with loss or reduced expression of PRKD1 had altered content, and injection of these sEVs into mice increased metastasis of xenograft tumors to lung, compared with sEVs from pancreatic cells that expressed PRKD1. PRKD1-deficient pancreatic cancer cells showed increased loading of integrin α6ß4 into sEVs-a process that required CD82. CONCLUSIONS: Human pancreatic ductal adenocarcinoma has reduced levels of PRKD1 compared with nontumor pancreatic tissues. Loss of PRKD1 results in reduced phosphorylation of cortactin in pancreatic cancer cell lines, resulting in increased in F-actin at the plasma membrane and increased release of sEVs, with altered content. These sEVs promote metastasis of xenograft and pancreatic tumors to lung in mice.

Metadata Correction: Patient Empowerment During the COVID-19 Pandemic by Ensuring Safe and Fast Communication of Test Results: Implementation and Performance of a Tracking System.

[This corrects the article DOI: 10.2196/27348.].

Patient Empowerment During the COVID-19 Pandemic by Ensuring Safe and Fast Communication of Test Results: Implementation and Performance of a Tracking System.

BACKGROUND: Overcoming the COVID-19 crisis requires new ideas and strategies for online communication of personal medical information and patient empowerment. Rapid testing of a large number of subjects is essential for monitoring and delaying the spread of SARS-CoV-2 in order to mitigate the pandemic's consequences. People who do not know that they are infected may not stay in quarantine and, thus, risk infecting others. Unfortunately, the massive number of COVID-19 tests performed is challenging for both laboratories and the units that conduct throat swabs and communicate the results. OBJECTIVE: The goal of this study was to reduce the communication burden for health care professionals. We developed a secure and easy-to-use tracking system to report COVID-19 test results online that is simple to understand for the tested subjects as soon as these results become available. Instead of personal calls, the system updates the status and the results of the tests automatically. This aims to reduce the delay when informing testees about their results and, consequently, to slow down the virus spread. METHODS: The application in this study draws on an existing tracking tool. With this open-source and browser-based online tracking system, we aim to minimize the time required to inform the tested person and the testing units (eg, hospitals or the public health care system). The system can be integrated into the clinical workflow with very modest effort and avoids excessive load to telephone hotlines. RESULTS: The test statuses and results are published on a secured webpage, enabling regular status checks by patients; status checks are performed without the use of smartphones, which has some importance, as smartphone usage diminishes with age. Stress tests and statistics show the performance of our software. CTest is currently running at two university hospitals in Germany-University Hospital Ulm and University Hospital Tübingen-with thousands of tests being performed each week. Results show a mean number of 10 (SD 2.8) views per testee. CONCLUSIONS: CTest runs independently of existing infrastructures, aims at straightforward integration, and aims for the safe transmission of information. The system is easy to use for testees. QR (Quick Response) code links allow for quick access to the test results. The mean number of views per entry indicates a reduced amount of time for both health care professionals and testees. The system is quite generic and can be extended and adapted to other communication tasks.

Biomarker profile for prediction of response to SMAC mimetic monotherapy in pediatric precursor B-cell acute lymphoblastic leukemia.

Second mitochondria-derived activator of caspase (SMAC) mimetics (SMs) targeting inhibitor of apoptosis proteins (IAPs) activate cell death pathways, and are currently being evaluated in clinical trials. Their successful therapeutic implementation requires upfront identification of patients who could benefit from a SM-based treatment but biomarkers for SM sensitivity have not yet been described. Here, we analyzed the intrinsic activity of two monovalent (AT406 and LCL161) and two bivalent (Birinapant and BV6) SMs on unselected patient-derived pediatric precursor B-cell acute lymphoblastic leukemia (BCP-ALL) identifying a subset of patient samples to be particularly sensitive to SM-induced cell death. This subset was defined by a characteristic gene expression signature with 127 differentially regulated genes, amongst them TNFRSF1A encoding TNFR1, and a critical role of TNFR1 in SM-induced cell death in sensitive BCP-ALL was confirmed on the functional level. Interestingly, samples with intermediate or low sensitivity to SMs were sensitized to SM-induced cell death by inhibition of caspases using zVAD.fmk or Emricasan, a pan-caspase inhibitor in clinical trials. When we compared our expression data to published data sets, we identified an overlap of four genes to be commonly differentially regulated in SM-sensitive BCP-ALL, that is, TSPAN7, DIPK1C, MTX2 and, again, TNFRSF1A. Functional testing revealed that this set of genes identified samples with high sensitivity to SM treatment. In summary, our data suggest using this gene signature as biomarker predicting response to SM treatment and point to the development of new combinatorial treatments consisting of SMs and pan-caspase inhibitors for a successful clinical implementation of SMs in treatment of BCP-ALL.

Patterns of antibody responses to nonviral cancer antigens in head and neck squamous cell carcinoma patients differ by human papillomavirus status.

There have been hints that nonviral cancer antigens are differentially expressed in human papillomavirus (HPV)-positive and HPV-negative head and neck squamous cell carcinoma (HNSCC). Antibody responses (AR) to cancer antigens may be used to indirectly determine cancer antigen expression in the tumor using a noninvasive and tissue-saving liquid biopsy. Here, we set out to characterize AR to a panel of nonviral cancer antigens in HPV-positive and HPV-negative HNSCC patients. A fluorescent microbead multiplex serology to 29 cancer antigens (16 cancer-testis antigens, 5 cancer-retina antigens and 8 oncogenes) and 29 HPV-antigens was performed in 382 HNSCC patients from five independent cohorts (153 HPV-positive and 209 HPV-negative). AR to any of the cancer antigens were found in 272/382 patients (72%). The ten most frequent AR were CT47, cTAGE5a, c-myc, LAGE-1, MAGE-A1, -A3, -A4, NY-ESO-1, SpanX-a1 and p53. AR to MAGE-A3, MAGE-A9 and p53 were found at significantly different prevalences by HPV status. An analysis of AR mean fluorescent intensity values uncovered remarkably different AR clusters by HPV status. To identify optimal antigen selections covering a maximum of patients with ≤10 AR, multiobjective optimization revealed distinct antigen selections by HPV status. We identified that AR to nonviral antigens differ by HPV status indicating differential antigen expression. Multiplex serology may be used to characterize antigen expression using serum or plasma as a tissue-sparing liquid biopsy. Cancer antigen panels should address the distinct antigen repertoire of HPV-positive and HPV-negative HNSCC.

Telomerase abrogates aneuploidy-induced telomere replication stress, senescence and cell depletion.

The causal role of aneuploidy in cancer initiation remains under debate since mutations of euploidy-controlling genes reduce cell fitness but aneuploidy strongly associates with human cancers. Telomerase activation allows immortal growth by stabilizing telomere length, but its role in aneuploidy survival has not been characterized. Here, we analyze the response of primary human cells and murine hematopoietic stem cells (HSCs) to aneuploidy induction and the role of telomeres and the telomerase in this process. The study shows that aneuploidy induces replication stress at telomeres leading to telomeric DNA damage and p53 activation. This results in p53/Rb-dependent, premature senescence of human fibroblast, and in the depletion of hematopoietic cells in telomerase-deficient mice. Endogenous telomerase expression in HSCs and enforced expression of telomerase in human fibroblasts are sufficient to abrogate aneuploidy-induced replication stress at telomeres and the consequent induction of premature senescence and hematopoietic cell depletion. Together, these results identify telomerase as an aneuploidy survival factor in mammalian cells based on its capacity to alleviate telomere replication stress in response to aneuploidy induction.

Central nervous system involvement in acute lymphoblastic leukemia is mediated by vascular endothelial growth factor.

In acute lymphoblastic leukemia (ALL), central nervous system (CNS) involvement is a major clinical concern. Despite nondetectable CNS leukemia in many cases, prophylactic CNS-directed conventional intrathecal chemotherapy is required for relapse-free survival, indicating subclinical CNS manifestation in most patients. However, CNS-directed therapy is associated with long-term sequelae, including neurocognitive deficits and secondary neoplasms. Therefore, molecular mechanisms and pathways mediating leukemia-cell entry into the CNS need to be understood to identify targets for prophylactic and therapeutic interventions and develop alternative CNS-directed treatment strategies. In this study, we analyzed leukemia-cell entry into the CNS using a primograft ALL mouse model. We found that primary ALL cells transplanted onto nonobese diabetic/severe combined immunodeficiency mice faithfully recapitulated clinical and pathological features of meningeal infiltration seen in patients with ALL. ALL cells that had entered the CNS and were infiltrating the meninges were characterized by high expression of vascular endothelial growth factor A (VEGF). Although cellular viability, growth, proliferation, and survival of ALL cells were found to be independent of VEGF, transendothelial migration through CNS microvascular endothelial cells was regulated by VEGF. The importance of VEGF produced by ALL cells in mediating leukemia-cell entry into the CNS and leptomeningeal infiltration was further demonstrated by specific reduction of CNS leukemia on in vivo VEGF capture by the anti-VEGF antibody bevacizumab. Thus, we identified a mechanism of ALL-cell entry into the CNS, which by targeting VEGF signaling may serve as a novel strategy to control CNS leukemia in patients, replacing conventional CNS-toxic treatment.

YAP Activation Drives Liver Regeneration after Cholestatic Damage Induced by Rbpj Deletion.

Liver cholestasis is a chronic liver disease and a major health problem worldwide. Cholestasis is characterised by a decrease in bile flow due to impaired secretion by hepatocytes or by obstruction of bile flow through intra- or extrahepatic bile ducts. Thereby cholestasis can induce ductal proliferation, hepatocyte injury and liver fibrosis. Notch signalling promotes the formation and maturation of bile duct structures. Here we investigated the liver regeneration process in the context of cholestasis induced by disruption of the Notch signalling pathway. Liver-specific deletion of recombination signal binding protein for immunoglobulin kappa j region (Rbpj), which represents a key regulator of Notch signalling, induces severe cholestasis through impaired intra-hepatic bile duct (IHBD) maturation, severe necrosis and increased lethality. Deregulation of the biliary compartment and cholestasis are associated with the change of several signalling pathways including a Kyoto Encyclopedia of Genes and Genomes (KEGG) gene set representing the Hippo pathway, further yes-associated protein (YAP) activation and upregulation of SRY (sex determining region Y)-box 9 (SOX9), which is associated with transdifferentiation of hepatocytes. SOX9 upregulation in cholestatic liver injury in vitro is independent of Notch signalling. We could comprehensively address that in vivo Rbpj depletion is followed by YAP activation, which influences the transdifferentiation of hepatocytes and thereby contributing to liver regeneration.

GiANT: gene set uncertainty in enrichment analysis.

UNLABELLED: Over the past years growing knowledge about biological processes and pathways revealed complex interaction networks involving many genes. In order to understand these networks, analysis of differential expression has continuously moved from single genes towards the study of gene sets. Various approaches for the assessment of gene sets have been developed in the context of gene set analysis (GSA). These approaches are bridging the gap between raw measurements and semantically meaningful terms.We present a novel approach for assessing uncertainty in the definition of gene sets. This is an essential step when new gene sets are constructed from domain knowledge or given gene sets are suspected to be affected by uncertainty. Quantification of uncertainty is implemented in the R-package GiANT. We also included widely used GSA methods, embedded in a generic framework that can readily be extended by custom methods. The package provides an easy to use front end and allows for fast parallelization. AVAILABILITY AND IMPLEMENTATION: The package GiANT is available on CRAN. CONTACTS: hans.kestler@leibniz-fli.de or hans.kestler@uni-ulm.de.

Sputnik: ad hoc distributed computation.

MOTIVATION: In bioinformatic applications, computationally demanding algorithms are often parallelized to speed up computation. Nevertheless, setting up computational environments for distributed computation is often tedious. Aim of this project were the lightweight ad hoc set up and fault-tolerant computation requiring only a Java runtime, no administrator rights, while utilizing all CPU cores most effectively. RESULTS: The Sputnik framework provides ad hoc distributed computation on the Java Virtual Machine which uses all supplied CPU cores fully. It provides a graphical user interface for deployment setup and a web user interface displaying the current status of current computation jobs. Neither a permanent setup nor administrator privileges are required. We demonstrate the utility of our approach on feature selection of microarray data. AVAILABILITY AND IMPLEMENTATION: The Sputnik framework is available on Github http://github.com/sysbio-bioinf/sputnik under the Eclipse Public License. CONTACT: hkestler@fli-leibniz.de or hans.kestler@uni-ulm.de SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Prediction of resistance to bevacizumab plus FOLFOX in metastatic colorectal cancer-Results of the prospective multicenter PERMAD trial.

BACKGROUND: Anti-vascular endothelial growth factor (VEGF) monoclonal antibodies (mAbs) are widely used for tumor treatment, including metastatic colorectal cancer (mCRC). So far, there are no biomarkers that reliably predict resistance to anti-VEGF mAbs like bevacizumab. A biomarker-guided strategy for early and accurate assessment of resistance could avoid the use of non-effective treatment and improve patient outcomes. We hypothesized that repeated analysis of multiple cytokines and angiogenic growth factors (CAFs) before and during treatment using machine learning could provide an accurate and earlier, i.e., 100 days before conventional radiologic staging, prediction of resistance to first-line mCRC treatment with FOLFOX plus bevacizumab. PATIENTS AND METHODS: 15 German and Austrian centers prospectively recruited 50 mCRC patients receiving FOLFOX plus bevacizumab as first-line treatment. Plasma samples were collected every two weeks until radiologic progression (RECIST 1.1) as determined by CT scans performed every 2 months. 102 pre-selected CAFs were centrally analyzed using a cytokine multiplex assay (Luminex, Myriad RBM). RESULTS: Using random forests, we developed a predictive machine learning model that discriminated between the situations of "no progress within 100 days before radiological progress" and "progress within 100 days before radiological progress". We could further identify a combination of ten out of the 102 CAF markers, which fulfilled this task with 78.2% accuracy, 71.8% sensitivity, and 82.5% specificity. CONCLUSIONS: We identified a CAF marker combination that indicates treatment resistance to FOLFOX plus bevacizumab in patients with mCRC within 100 days prior to radiologic progress.