Search Results - EGA European Genome-Phenome Archive

UK10K RARE COLOBOMA

In the UK10K project we propose a series of complementary genetic approaches to find new low frequency/rare variants contributing to disease phenotypes. These will be based on obtaining the genome wide sequence of 4000 samples from the TwinsUK and ALSPAC cohorts (at 6x sequence coverage), and the exome sequence (protein coding regions and related conserved sequence) of 6000 samples selected for extreme phenotypes. Our studies will focus primarily on cardiovascular-related quantitative traits, obesity and related metabolic traits, neurodevelopmental disorders and a limited number of extreme clinical phenotypes that will provide proof-of-concept for future familial trait sequencing. We will analyse directly quantitative traits in the cohorts and the selected traits in the extreme samples, and also use imputation down to 0.1% allele frequency to extend the analyses to further sample sets with genome wide genotype data. In each case we will investigate indels and larger structural variants as well as SNPs, and use statistical methods that combine rare variants in a locus or pathway as well as single-variant approaches. Ocular coloboma is the most common significant developmental eye defect with an incidence of ~1 in 5,000 live births. It results from failure of optic fissure closure during embryogenesis. The position and extent of the fusion failure dictates the clinical appearance and functional effect. ~30% of coloboma cases are associated with other systemic malformations. These UK10K samples will mostly comprise isolated coloboma cases without systemic involvement (aka non-syndromal coloboma). There is strong evidence from family studies that coloboma has a major genetic component with autosomal dominance being the most common pattern of inheritance. However, many cases are isolated or show complex patterns of familial clustering. The genes responsible for isolated coloboma are largely unknown, but in a small number of families mutations in SHH, CHX10, and PAX6 have been identified indicating marked genetic heterogeneity. Thus in addition to the clinical benefits of achieving a molecular diagnosis there are also major scientific advantages to identifying coloboma genes, as these are likely to provide insights into the complex process of optic fissure closure, that is critical to normal eye development. In the longer term, understanding the molecular basis of the disease may provide clues to therapeutic strategies.For further information with regard to this cohort please contact David Fitzpatrick (david.fitzpatrick@ed.ac.uk).

Study EGAS00001000127

RNAseq dataset for MALT1 in MCL

Manuscript Title: Co-targeting of BTK and MALT1 overcomes resistance to BTK inhibitors in mantle cell lymphoma Journal: Journal of Clinical Investigation Authors Vivian Changying Jiang1, Yang Liu1, Junwei Lian1, Shengjian Huang1, Alexa Jordan1, Qingsong Cai1, Fangfang Yan3, Joseph Mitchell McIntosh1, Yijing Li1, Yuxuan Che1, Zhihong Chen1, Jovanny Vargas1, Maria Badillo1, JohnNelson Bigcal1, Heng-Huan Lee1, Wei Wang1, Yixin Yao1, Lei Nie1, Christopher Flowers1, and Michael Wang1, 2* Abstract Bruton’s tyrosine kinase (BTK) is a proven target in mantle cell lymphoma (MCL), an aggressive subtype of non-Hodgkin lymphoma. However, resistance to BTK inhibitors is a major clinical challenge. We here report that MALT1 is one of the top overexpressed genes in ibrutinib-resistant MCL cells, while expression of CARD11, which is upstream of MALT1, is decreased. MALT1 genetic knockout or inhibition produced dramatic defects in MCL cell growth regardless of ibrutinib sensitivity. Conversely, CARD11 knockout cells showed anti-tumor effects only in ibrutinib-sensitive cells, suggesting that MALT1 overexpression could drive ibrutinib resistance via bypassing BTK-CARD11 signaling. Additionally, BTK knockdown and MALT1 knockout markedly impaired MCL tumor migration and dissemination, and MALT1 pharmacological inhibition decreased MCL cell viability, adhesion, and migration by suppressing NF-κB, PI3K-ATK-mTOR, and integrin signaling. Importantly, co-targeting MALT1 with safimaltib and BTK with pirtobrutinib induced potent anti-MCL activity in ibrutinib-resistant MCL cell lines and patient-derived xenografts. Therefore, we conclude that MALT1 overexpression associates with resistance to BTK inhibitors in MCL, targeting abnormal MALT1 activity could be a promising therapeutic strategy to overcome BTK inhibitor resistance, and co-targeting of MALT1 and BTK should improve MCL treatment efficacy and durability as well as patient outcomes. Dataset description: The bulk RNA-seq dataset was generated for the cell lines below and used for two major purposes: 1. DEG analysis and GSEA analysis comparing IBN-R and IBN-S cells 2. DEG analysis and GSEA analysis comparing MCL cells with/without MI-2 treatment. sample Cell MI-2 Ibrutinib (IBN) Venetoclax (VEN) Used for IBN-R vs IBN-S comparison Used for MI-2 vs untreated (DMSO) H9 Granta519 - R S yes H21 Granta519 - R S yes H33 Granta519 - R S yes H10 Granta519-VEN-R - R R yes H22 Granta519-VEN-R - R R yes H34 Granta519-VEN-R - R R yes H3 JeKo BTK KD_1 - R R yes yes H15 JeKo BTK KD_1 - R R yes yes H27 JeKo BTK KD_1 - R R yes yes H5 JeKo BTK KD_2 - R R yes yes H17 JeKo BTK KD_2 - R R yes yes H29 JeKo BTK KD_2 - R R yes yes H1 JeKo-1 - S R yes yes H13 JeKo-1 - S R yes yes H25 JeKo-1 - S R yes yes H7 Mino - S S yes H19 Mino - S S yes H31 Mino - S S yes H8 Mino-VEN-R - S R yes H20 Mino-VEN-R - S R yes H32 Mino-VEN-R - S R yes H11 Rec-1 - S S yes H23 Rec-1 - S S yes H12 Rec-VEN-R - S S yes H24 Rec-VEN-R - S R yes H36 Rec-VEN-R - S R yes H35 Rec-1 -- S R yes H4 JeKo BTK KD_1 + MI-2 + yes H16 JeKo BTK KD_1 + MI-2 + yes H28 JeKo BTK KD_1 + MI-2 + yes H6 JeKo BTK KD_2 + MI-2 + yes H18 JeKo BTK KD_2 + MI-2 + yes H30 JeKo BTK KD_2 + MI-2 + yes H2 JeKo-1 + MI-2 + yes H14 JeKo-1 + MI-2 + yes H26 JeKo-1 + MI-2 + yes

Dataset EGAD00001009771

The patterns and dynamics of genomic instability in metastatic pancreatic cancer

Pancreatic cancer is an aggressive malignancy with a five-year mortality of 97â€“98%, usually due to widespread metastatic disease. Previous studies indicate that this disease has a complex genomic landscape, with frequent copy number changes and point mutations, but genomic rearrangements have not been characterized in detail. Despite the clinical importance of metastasis, there remain fundamental questions about the clonal structures of metastatic tumours, including phylogenetic relationships among metastases, the scale of ongoing parallel evolution in metastatic and primary sites, and how the tumour disseminates. Here we harness advances in DNA sequencing to annotate genomic rearrangements in 13 patients with pancreatic cancer and explore clonal relationships among metastases. We find that pancreatic cancer acquires rearrangements indicative of telomere dysfunction and abnormal cell-cycle control, namely dysregulated G1-to-S-phase transition with intact G2â€“M checkpoint. These initiate amplification of cancer genes and occur predominantly in early cancer development rather than the later stages of the disease. Genomic instability frequently persists after cancer dissemination, resulting in ongoing, parallel and even convergent evolution among different metastases. We find evidence that there is genetic heterogeneity among metastasis-initiating cells, that seeding metastasis may require driver mutations beyond those required for primary tumours, and that phylogenetic trees across metastases show organ-specific branches. These data attest to the richness of genetic variation in cancer, brought about by the tandem forces of genomic instability and evolutionary selection.

Study EGAS00000000064

Phase II Therapeutic Trial of Mexiletine in Non-Dystrophic Myotonia (IND #77,021)

Non-dystrophic myotonias are a group of rare muscle disorders caused by abnormalities in different muscle cell membrane proteins. Patients experience delayed muscle relaxation that causes impaired physical activity, stiffness, and pain. Mexiletine has been shown to be beneficial for symptoms of myotonia in case reports of patients with non-dystrophic myotonia, however, current management of these patients primarily depends on individual practitioners' preference. This study is a placebo-controlled, double-blind, cross-over treatment trial using mexiletine in non-dystrophic myotonias. The results, if positive, will permit the identification of mexiletine as an effective therapy for myotonia, making it a first-line therapy for patients with non-dystrophic myotonia. The specific aim of this proposal is to perform a randomized, double-blind, placebo-controlled cross-over study to assess whether mexiletine improves both quantitative and qualitative measures of myotonia in patients with non-dystrophic myotonia. Given the rarity of non-dystrophic myotonia, large randomized controlled treatment trials in these disorders have not been feasible in the past. Primary Outcome Measure: Patient-assessed symptom: stiffness as measured by an Interactive Voice Response Diary (IVR). Secondary Outcome Measures: a) patient-assessed pain, weakness, and fatigue as measured by IVR; b) clinical myotonia assessment; c) quality of life as measured by INQoL, SF36; d) a quantitative measure of grip myotonia; e) measurement of CMAP after short and long exercise; f) grading of myotonia on needle EMG.

Study phs001311

Genetics of Disorders Affecting Tooth Structure, Number, Morphology and Eruption

The long-term objective of this research project is to discover the complete spectrum of genes that contribute to the etiologies of genetic disorders of dentition. Inherited tooth defects can be isolated (the only symptom) or syndromic (having additional symptoms), and can be caused by mutations in many different genes. Clinically, it can be difficult to distinguish among the diverse forms. Making an accurate diagnosis would be greatly improved if reliable tests could identify the specific genetic defect that causes the disease in each case. Current technology can identify the numerous sequence variations throughout a patient's exome, but the disease-causing variation among them cannot be identified unless it is already on the list of genes that potentially cause the disease. A diagnosis based upon the genetic defect that causes the disorder distinguishes between isolated and syndromic tooth defects. A precise diagnosis permits formulation of a more accurate prognosis, which informs treatment decisions. The benefits of accomplishing our aims are to 1) enable genetic testing to achieve a precise and accurate diagnosis, and thereby optimize treatment, and 2) improve our understanding of normal and pathological tooth formation, leading to future improvements in treatment and eventually cures. Our short-term objectives are to 1) recruit kindreds with inherited disorders affecting their dentition, characterize their dental phenotype (and non-dental phenotypes when applicable) and the distribution of the disorder throughout the family, and 2) determine the genetic etiology of their disorders.

Study phs001491

National Institute on Alcohol Abuse and Alcoholism (NIAAA) Postmortem Prefrontal Cortex eQTL and mQTL Study

Mapping expression and methylation quantitative trait loci (eQTLs and mQTLs) using genome-wide genotype, transcriptome, and DNA methylome data can provide an unbiased approach to identify novel functional variants. One critical issue of the eQTL/mQTL analysis is that the role of genetic variants in regulating gene expression (or DNA methylation) can be tissue- or cell-specific. To unravel causal variants from disease-associated loci, it would be necessary to map eQTLs and mQTLs in both patients and controls. Those eQTLs and mQTLs in patients are more likely to be causal variants for the disease. Identification of disease- and tissue-specific eQTLs and mQTLs will help disentangle the role of the specific tissue in disease development as well. In the present study, we mapped eQTLs and mQTLs in postmortem prefrontal cortex of patients with alcoholism or alcohol use disorder (AUD) and matched control subjects by integrative analyses of genome-wide SNP genotype, mRNA transcriptome, and DNA methylome data. As AUD is comorbid with schizophrenia (SCZ), bipolar disorder (BIP), major depressive disorder (MDD), and attention deficit/hyperactivity disorder (ADHD), we further examined the co-localization of PFC eQTLs/mQTLs and GWAS-identified SNPs (mainly noncoding SNPs) that were associated with AUD as well as comorbid SCZ, BIP, MDD, and ADHD. Our study demonstrated the enrichment of PFC eQTLs/mQTLs in the above psychiatric disorder-associated genetic variants identified by genome-wide association studies (GWAS).

Study phs001981

A Phase I/II Study of Revlimid (lenalidomide) in Combination with Vidaza (azacitidine) in Patients with Advanced Myelodysplastic Syndrome (MDS)

This single-arm study design has been selected to investigate the safety of combination therapy with Vidaza® (azacitidine) and Revlimid® (lenalidomide) in a population of 18 participants with advanced MDS for the Phase I portion, and 18 participants for the Phase II portion. The multi-center nature of the study should enhance the general applicability of the results. Participants with advanced MDS were chosen because it is expected that the combination of the two drugs will by myelosuppressive, which can be considered a necessary toxicity to effect improvements in PR and CR rates in this advanced population, akin to participants with AML. Finally, the combination was chosen to be studied because currently available therapies in this population, and their effectiveness, are limited, and single-agent Vidaza® (azacitidine), while effective at delaying transformation to AML or death, is not curative. The doses of Revlimid® (lenalidomide) in this study are similar to those that were demonstrated to have an effect and to be well-tolerated in participants with MDS. Doses for Vidaza® (azacitidine) are similar to those used in the pivotal trial, though for a 5-day treatment period instead of 7 to minimize toxicity for the first 3 blocks of participants, and then at a lower dose over a protracted period of time to explore the safety of a 10-day treatment course.

Study phs001318

University of Texas PDX Development and Trial Center Grant

The goal for the University Texas PDX Development and Trial Center (UTPDTC) is to optimize personalized biomarker-based cancer therapy and identify effective targeted drugs based on the molecular characteristics of each tumor. Our short-term goals are to establish a biobank of clinically, and molecularly-annotated Patient-Derived Xenografts (PDXs) and to use PDXs as a platform for preclinical drug development and biomarker discovery. The primary goal for UTPDTC investigators will be to develop PDX trial strategies for preclinical testing of single agents and drug combinations. These models will allow the determination of the optimal treatments (single drugs or combinations) that should be tested in clinical trials in increasingly individualized, molecularly defined subsets of tumors. The goal of the Patient-Derived Xenograft Core is to provide high-quality clinically relevant and molecularly annotated PDX models for the research projects proposed in the University of Texas PDX Development and Trial Center (UTPDTC) grant application and to the research activity of the NCI PDX Development and Trial Centers Research Network (PDXNet) by leveraging PDX resources at our institutions and developing new PDX models from human cancer specimens using rigorous quality standards so that the models can be used to guide clinical trial development. The PDX models developed and/or characterized by the PDX Core will be available to other cancer researchers through the PDXNet and NCI Patient-Derived Models Repository (PDMR).

Study phs001980

Whole Genome Sequencing Analysis in a Family of Discordant Twins With Non-Syndromic Microtia and Hemifacial Microsomia: Identification of Novel Candidate Genes and Variants

This research is a non-NIH effort to discover the pathogenesis of non-syndromic microtia and to detect potential candidate genes. Microtia (OMIM 600674, 612290, 251800) is a spectrum of congenital anomalies and facial abnormalities, 90% of the cases followed by hearing loss. Two types of the disease are described in terms of connection to the affected syndromes: anomalies (syndromic) and only clinical manifestation (non-syndromic) associated with other symptoms. Though the exact pathogenesis of microtia remains unknown, pathogenic genes of majority of the microtia syndromic types have been identified. However, the genetic mutation of the non-syndromic microtia has not been confirmed. In this study, whole genome sequencing (WGS) was utilized to screen mutations in monozygotic discordant twins of congenital microtia-atresia and their family members. One of the twins, along with her grandfather, was discordant for microtia-atresia. By bioinformatic analysis of the WGS data from monozygotic twins with microtia and their non-microtia relatives, we proposed that among 28 identified potential genes, 4 of them (CNOT1, ODAD4, TBX15, and GIPC3) could be more responsible for non-syndromic microtia development and phenotypic features. Interested investigators may request to access the DNA genomic sequencing data of the study participants via dbGaP authorized access. In addition to raw genomic sequencing data, datasets of germline mutations and candidate gene mutations may also be available for investigators to request.

Study phs003216

Small Intestine Neuroendocrine Tumors (Carcinoid Tumors)

The diagnosed incidence of small bowel neuroendocrine tumors (NETs) is increasing. While patients with localized disease can be treated surgically, those with metastatic disease currently have few treatment options. The success of biologically targeted therapies in other malignancies has led to interest in the molecular alterations underlying the pathogenesis of these NETs. To identify genetic aberrations in small intestine NETs, we generated copy number profiles from 31 primary and metastatic tumors and performed whole-exome sequencing on a subset of 29 primary small intestine NETs and 24 metastatic NETs in parallel with normal blood DNA. Whole-genome sequencing data was generated on 15 tumor/normal pairs and 5 primary/metastasis/normal trios. The global genetic landscape of small bowel NETs is relatively quiet. Consistent with previous studies, the overwhelming majority of tumors were characterized by loss of chromosome 18 and, to a lesser extent, other chromosome arm gains and losses. In stark contrast to arm-level alterations, recurrent high-level focal amplifications and deletions were much less prevalent in these tumors. High-throughput mutation screening and exome sequencing revealed similarly low rates of somatic mutation in NETs (median of 0.77 non-silent mutations per megabase (Mb) of coding DNA) compared to other recent cancer exome sequencing efforts. Our analysis of this cohort identified only a single, statistically significant recurrent somatic mutation targeting the cyclin-dependent kinase inhibitor gene, CDKN1B, encoding p27.

Study phs000579

The Genetics of Type 2 Diabetes Consortium (GoT2D): Low-Pass Sequencing and High-Density SNP Genotyping for Type 2 Diabetes

The Genetics of Type 2 Diabetes Consortium (GoT2D) is a collaboration between the University of Michigan, the Broad Institute and the Wellcome Trust Centre for Human Genetics. The overall aim is to extend upon recent efforts, such as genome-wide association studies (GWAS) and large scale meta-analyses. While they have proved successful at mapping genomic loci that influence human diseases, like type 2 diabetes, much of the heritability remains unexplained. In this study, we use next generation sequencing and genotyping technologies to query for lower frequency variants in the human genome. Thereby, allowing a deeper characterization of the spectrum of alleles associated with type 2 diabetes risk, and a better assessment of the genes that play a role in the etiology of type 2 diabetes development. We studied 1,326 T2D cases and 1,331 normoglycemic controls from Northern and Central Europe (Sweden, Finland, UK, and Germany). To efficiently characterize the entire genome sequence of each individual, we performed low-coverage (~5x) whole-genome sequencing, augmented by deep coverage (~100x) sequencing of the exome (Fuchsberger et al, 2016), and dense (2.5M) single nucleotide polymorphism (SNP) genotyping using the HumanOmni2.5 array. The data deposited in dbGaP will include all the Swedish, Finnish, and UK samples, but the German data will be deposited in the European Genome-phenome Archive (EGA), by virtue of the project specific funding requirements.

Study phs000840

Use of Whole Genome Sequencing for Diagnosis and Discovery in the Cancer Genetics Clinic

Despite the potential of whole-genome sequencing (WGS) to improve patient diagnosis and care, the empirical value of WGS in the cancer genetics clinic is unknown. We performed WGS on members of two cohorts of cancer genetics patients: those with BRCA1/2 mutations (n = 176) and those without (n = 82). Initial analysis of potentially pathogenic variants (PPVs, defined as nonsynonymous variants with allele frequency < 1% in ESP6500) in 163 clinically-relevant genes suggested that WGS will provide useful clinical results. This is despite the fact that a majority of PPVs were novel missense variants likely to be classified as variants of unknown significance (VUS). Furthermore, previously reported pathogenic missense variants did not always associate with their predicted diseases in our patients. This suggests that the clinical use of WGS will require large-scale efforts to consolidate WGS and patient data to improve accuracy of interpretation of rare variants. While loss-of-function (LoF) variants represented only a small fraction of PPVs, WGS identified additional cancer risk LoF PPVs in patients with known BRCA1/2 mutations and led to cancer risk diagnoses in 21% of non-BRCA cancer genetics patients after expanding our analysis to 3209 ClinVar genes. These data illustrate how WGS can be used to improve our ability to discover patients' cancer genetic risks. "Reprinted from doi:10.1016/j.ebiom.2014.12.003, with permission from EBioMedicine."

Study phs000942

Uncovering Inversion Formation in the Human Genome and its Impact to Disease

The relevance of inversions for disease causation, speciation and adaptation is broadly recognized, although the prevalence of inversions is unknown. In humans, de novo inversions are associated with congenital anomalies in approximately 9.6% of patients. Yet, despite the biological relevance of inversions, their molecular features, formation mechanism, impact to the genomic structure of carriers, as well as their contribution to clinical phenotypes, have not been further explored. Inversions are typically classified as balanced reciprocal events generated by ectopic recombination, although recent studies reveal a distinct picture whereby inversions originate from mechanisms that concomitantly generate copy number variants (CNVs). Surprisingly, those complex inversions underlie as much as 30% of neurodevelopmental-associated CNVs. The hypothesis of this project are: inversions are often generated de novo by mechanisms other than ectopic recombination a relevant fraction of inversions are associated with complex genomic rearrangements (CGRs), often overlooked in sporadic diseases, and inversions are a "hidden" type of structural variation for which contribution to a clinical phenotypes has been under assessed due to the lack of appropriate detection tools. A combined strategy using multiple genomic tools will be applied to characterize inversions and associated genomic alterations in individuals with neurodevelopmental disorders and controls: whole genome sequencing (WGS) with short-and long-reads; genome mapping; classical cytogenetics; and array CGH. This project will establish common ground to bridge studies of rare and common diseases, human evolution, and population genetics.

Study phs002999

Genomic and Phenotypic Profile of Sickle Cell Disease in Human Population in Cameroon

Sickle cell disease (SCD) is caused by a biallelic single nucleotide substitution in the beta-globin gene (HBB), resulting in an amino acid substitution, Glu7Val, formerly known as Glu6Val. An estimated 300,000 babies with SCD are born worldwide each year, with nearly 75% of these births being in sub-Saharan Africa. In Africa, at least 30-50% of children with untreated SCD die before the age of 5 years. Therefore, accelerating the path for novel therapies for SCD through genomics research on fetal hemoglobin (HbF) is critical. Variants in the currently known HbF‐modulating genes/loci, i.e., BCL11A, HBS1L-MYB, and XmnI-HBG2, explain only 10-20% of the variation of HbF levels in African individuals with SCD, compared with nearly 50% of the variation in HbF levels among Europeans. Expanding genomic research in populations of African ancestry could uncover the missing heritability of HbF-promoting loci. This study used the Human Heredity and Health (H3Africa) consortium SNP genotyping array developed from whole genome enriched for common variants in sub-Saharan Africans to investigate genomic variations associated with HbF levels in individuals with sickle cell anemia of the HbSS genotype from Cameroon. Two of the previously described loci, BCL11A and HBS1L-MYB were replicated, while a novel locus on chromosome 13 mapping to FLT1 was uncovered. The study therefore opens up a novel avenue for exploring SCD gene therapy.

Study phs003748

Biopsy-Derived Organoids in Personalised Early Breast Cancer Care: Challenges of Tumour Purity and Normal Cell Overgrowth Cap Their Practical Utility

The ability to establish organoids composed exclusively of tumour rather than healthy cells is essential for their implementation into clinical practice. Organoids have in recent years emerged as a powerful tool to expand patient material in culture and generate modifiable 3D models derived from humans or animal models. For translational research, they enable the creation of model systems for an ever-increasing number of cell types and diseases. And in personalised medicine, they potentially allow for functional drug testing with high predictive power in certain settings. We found that using biopsy material from untreated, early-stage primary breast cancer patients poses significant challenges for consistently culturing tumour cells as organoids. Specifically, we observed frequent outgrowth of genetically normal, non-cancerous epithelial cells. We analysed >100 biopsy samples from early-stage breast cancer and present our large collection of >70 organoid lines. We also show methods of assessing successful tumour cell culture in a time and cost-efficient manner, proving a high rate (>85%) of normal cell overgrowth in early-stage breast cancer organoids. Finally, we show a number of successful attempts to culture cancer organoids from mastectomy-derived tissue of advanced, metastatic breast cancer. We conclude that the usefulness of organoids from early breast cancer for translational research and personalised medicine, especially guidance of adjuvant or post-surgical maintenance therapy, is strongly limited by the low success rate of culturing cancerous cells under organoid conditions.

Study EGAS50000000605

Multi-omics bulk and single-cell profiling of epithelioid sarcoma

Epithelioid sarcoma (EpS) is an aggressive sarcoma, characterized by the loss of SMARCB1 expression. EpS is traditionally classified as distal or proximal according to clinicopathological features, but its molecular characteristics remain largely unknown. To establish an EpS molecular classification and uncover determinants of inter- and intra-patient heterogeneity in EpS, we used multi-omics profiling and integrated the genomic, transcriptional and methylome landscapes with single-cell RNA sequencing on fresh samples as well as spatial transcriptomics. We identified two molecular subtypes of EpS: “distal-like” and “proximal-like”, which were distinct from the histological subtypes. Distal-like tumors include all morphological distal subtypes and a subset of proximal and hybrid EpS. They express a specific molecular single-cell derived epithelial-to-mesenchymal transition signature, which associates with improved patient survival and includes Desmoglein 2 (DSG2), which we identify as a potential routine immunohistochemical biomarker to improve diagnosis. Distal-like EpS also display increased peri-tumoral CD8+ T cell infiltrates and specific tumor-immune cell interactions. Conversely, proximal-like EpS harbor a higher inter-tumoral molecular heterogeneity with some cases resembling other SMARCB1-deficient tumors by DNA methylation profiling, and increased intra-tumoral pro-tumoral macrophage infiltrates. Our study introduces a novel molecular classification of EpS with prognostic value, associated diagnostic biomarkers and distinct tumor immune microenvironments, which extends beyond traditional clinicopathological classification. This paves the way for precision medicine-based therapeutic strategies for patients with EpS.

Study EGAS50000000973

Collagen XVII Promotes Pancreatic Cancer Through Regulation of PIK3R5

Pancreatic ductal adenocarcinoma (PDAC) contains a desmoplastic stroma that promotes the aggressive biology of this disease. While the varied mechanisms by which fibrillar collagens contribute to PDAC have been defined, little is known about the contribution of transmembrane collagens. Among this class of collagens, Collagen XVII uniquely contains an extensive intracellular domain that connects intracellular cytoplasmic components with the extracellular matrix. Collagen XVII is best characterized as the defining component of the type I hemidesmosome, an anchoring complex that attaches epithelial cells to the underlying basement membrane. In this study, we characterize the expression of collagen XVII in PDAC and the mechanisms by which it contributes to PDAC biology. We demonstrate that crosstalk between PDAC cells and cancer associated fibroblasts drive collagen XVII expression in vitro and xenograft tumor models in vivo. Genetic depletion of collagen XVII in low-passage patient derived cell lines diminishes their metastatic potential and ability to form tumors in mice. Whole transcriptome sequencing of tumors formed by cells with diminished collagen XVII levels revealed dramatically lower levels of PIK3R5, a regulatory subunit of class IB phosphatidylinositol-3-kinase (PI3K) that interacts with G protein-coupled receptors. Consistent with this, collagen XVII-deficient tumors exhibited lower levels of phosphorylated AKT, a downstream effector of PI3K. The studies highlight the importance of collagen XVII signaling in PDAC and define a novel mechanism by which collagen XVII mediates activation of PI3K.

Study phs003641

Field studies of Cryptosporidiosis and Enteropathogens in Bangladesh

Cryptosporidiosis causes severe diarrhea in infants in the developing world. There is no vaccine to prevent it, and little in the way of treatment. This study on Bangladeshi urban slum children aims to support the design of a vaccine, both by determining how the immune system protects from infection and by identifying the genotypes of the parasite that should be included in a vaccine, as well as aid in development of therapies by identifying human genes that control infection. The primary objective of the study was to determine the incidence and contribution to disease of the different species and genotypes of cryptosporidia. Secondary objectives were designed to determine acquired immune response to cryptosporidiosis and identify human genes that influence susceptibility to cryptosporidiosis. This is an observational study. Children were recruited from Mirpur Dhaka slum (Cohort 1) and rural Mirzapur (Cohort 2) and followed for cryptosporidium infection longitudinally from birth through age 4 years in Cohort 1 and birth through age 2 years in Cohort 2. Biweekly household visits for diarrheal surveillance were made in addition to anthropometric measurements of mother and child, blood samples collected two times each year from the child, a work-up of diarrheal stools and non-diarrheal surveillance stools for cryptosporidium and other enteropathogens, and blood and breast milk samples from the mother. A detailed description of the study design and procedures can be obtained from publication: Kevin L Steiner, et al., 2018, PMID: 29897482.

Study phs001665

RNA Transcriptomic and DNA Methylation Landscape in FTLD-TDP and Controls

Frontotemporal lobar degeneration (FTLD) is one of the leading causes of dementia in individuals younger than 65 years, with the aggregation of TDP-43 as one of the most common forms of FTLD. FTLD-TDP is clinically, genetically and pathologically heterogeneous, with GRN and C9orf72 as the most common genetic forms (although more than 50% of the cases are genetically unexplained), and A, B and C as the most common pathological subtypes. To investigate the molecular differences between the different forms of FTLD-TDP, we performed bulk RNA sequencing (RNAseq) and reduced representation bisulfite sequencing (RRBS) from different brain regions of FTLD-TDP patients and controls. The RNAseq dataset includes samples from the frontal cortex of 149 individuals, comprising 127 FTLD-TDP cases (genetically unexplained types A (27), B (20) and C (22), 24 GRN mutation carriers, and 34 C9orf72 expansion carriers) and 22 controls. The RRBS dataset includes samples from the frontal cortex of 158 individuals —119 FTLD-TDP cases (genetically unexplained types A (25), B (24), and C (25), 20 GRN mutation carriers, and 25 C9orf72 expansion carriers) and 39 controls; and from the cerebellum of 151 individuals —113 FTLD-TDP cases (genetically unexplained types A (24), B (21), and C (23), 20 GRN mutation carriers, and 25 C9orf72 expansion carriers) and 38 controls. In total, the study consists of data from 183 unique individuals, with some overlap across the different datasets.

Study phs004075

Gut metagenome associations with extensive digital health data in a volunteer-based EstMB cohort second timepoint data

As part of the Estonian Biobank, we established the Estonian Microbiome Cohort which includes stool, oral and plasma samples from 2,509 participants and is supplemented with multi-omic measurements, questionnaires, and regular linkages to national electronic health records. Out of the 2,509 Estonian Microbiome Project participants, a sub-cohort of over 300 individuals provided an additional stool sample after a median follow-up period of 4.4 years. All participants of the EstMB cohort gave informed consent for the data and samples to be used for scientific purposes. In the second time point, the participants were instructed to send the samples by post immediately after taking the sample and time their sample collection to avoid shipping on the weekends. In the study, by utilising retrospective medication usage data from the electronic health records and the Estonian microbiome cohort shotgun metagenomics data set (n = 2,509), we systematically evaluate the long-term effects of antibiotics and human-targeted medications on the gut microbiome. We show that past usage of medications is associated with the gut microbiome. For example, the effects of antibiotics, psycholeptics, antidepressants, proton pump inhibitors, and beta-blockers are detectable several years after use. Furthermore, by analysing a subcohort (n >300 individuals) from the second time point, we show that similar changes in the gut microbiome occur after treatment initiation or discontinuation, possibly indicating causal effects.

Study EGAS50000001181

Genetic alterations that deregulate RB and PDGFRA signaling pathways drive tumor progression in IDH2-mutant astrocytoma

In IDH-mutant astrocytomas, IDH2 mutation is quite rare and biological mechanisms underlying tumor progression in IDH2-mutant astrocytoma remain elusive. Here, we report a unique case of IDH2 mutant astrocytoma, CNS WHO grade 3 that developed tumor progression. We performed a comprehensive genomic and epigenomic analysis for primary and recurrent tumors and found that both tumors harbored recurrent IDH2R172K and TP53R248W mutation with CDKN2A/B hemizygous deletion. We also found amplifications of CDK4 and MDM2 with PDGFRA gain in the recurrent tumor and upregulated protein expressions of these genes. We further developed, for the first time, a xenograft mouse model of IDH2R172K mutant astrocytoma from the recurrent tumor, but not from the primary tumor. Consistent with parent recurrent tumor cells, amplifications of CDK4 and MDM2 with PDGFRA gain was found, while CDKN2A/B was identified as homozygous deletion in the xenografts, qualifying for integrated diagnosis of astrocytoma, IDH2-mutant, CNS WHO grade 4. Cell viability assay found that CDK4/6 inhibitor and PDGFR inhibitor potently decreased cell viability in recurrent tumor cells, as compared to primary tumor cells. These findings suggest that gene alterations that activate retinoblastoma (RB) signaling pathways and PDGFR may drive tumor progression and xenograft formation in IDH2-mutant astrocytoma, which is equivalent to progressive IDH1-mutant astrocytoma. Also, our findings suggest that these genomic alterations may represent therapeutic targets in IDH2-mutant astrocytoma.

Study JGAS000646

Comprehensive molecular and clinicopathological profiling of lung adenocarcinoma in Japanese never or light smokers

Studies have yet to characterize the differences in molecular profiles of lung adenocarcinoma (LUAD) among divergent ethnic groups. Herein, we conducted comprehensive molecular profiling of LUAD in never or light smokers from Asia in order to discover novel targetable mutations and prognostic biomarkers of this distinct disease entity. Methods: We analyzed 996 cases of Asian LUAD, and performed whole-exome sequencing and/or RNA-seq in 125 cases of Asian LUAD negative for the driver oncogenes defined by conventional laboratory testing. We also investigated the clinical and pathological characteristics among the 996 cases. Results: Driver oncogenes were identified in 88 cases (70.4%) with specific hotspot mutations differing from those in The Cancer Genome Atlas (TCGA) study. Clustering based on mRNA expression profiles, but not genetic mutational ones, could predict patient prognosis. The risk score generated by the expression of a three-gene set was a strong prognostic marker for overall survival and progression-free survival in our cohort, and was further validated using TCGA cohort. Among the 996 cases, each driver alteration is distributed across all histological subtypes. Adenocarcinoma in situ was identified to harbor driver mutations, suggesting that these alterations are early events in the pathogenesis of LUAD. This study indicates the value of applying gene expression profiling for predicting the prognosis after a surgical operation, and that the identification of actionable mutations is important for optimizing targeted drugs in Japanese LUAD.

Study JGAS000215

DNA methylation at HBV integrants and flanking host genomes

Integration of DNA viruses into the human genome plays an important role in various types of tumors, including hepatitis B virus (HBV)-related hepatocellular carcinoma. However, the molecular details and clinical impact of HBV integration on either the human or HBV epigenomes are unknown. Here, we show that methylation of the integrated HBV DNA is related to the methylation status of the flanking human genome. We developed a next-generation sequencing-based method for structural methylation analysis of integrated viral genomes (denoted G-NaVI). This method is a novel approach that enables enrichment of viral fragments for sequencing using unique baits based on the sequence of the HBV genome. We detected integrated HBV sequences in the genome of the PLC/PRF/5 cell line and found variable levels of methylation within the integrated HBV genomes. Allele-specific methylation analysis revealed that the HBV genome often became significantly methylated when integrated into highly methylated host sites. After integration into unmethylated human genome regions such as promoters, however, the HBV DNA remains unmethylated and may eventually play an important role in tumorigenesis. The observed dynamic changes in DNA methylation of the host and viral genomes may functionally affect the biological behavior of HBV. These findings may impact public health given that millions of people worldwide are carriers of HBV. We also believe our assay will be a powerful tool to increase our understanding of the various types of DNA virus-associated tumorigenesis.

Study JGAS000015

Relapse CHL study

Classical Hodgkin Lymphoma (CHL) is characterized by a complex tumor microenvironment (TME) that supports disease progression. While immune cell recruitment by Hodgkin and Reed–Sternberg (HRS) cells is well-documented, the role of non-malignant B cells in relapse remains unclear. Using single-cell RNA sequencing (scRNA-seq) on paired diagnostic and relapsed CHL samples, we identified distinct shifts in B cell populations, particularly an enrichment of naïve B cells and a reduction of memory B cells in early-relapse compared to late-relapse and newly diagnosed CHL. The enrichment of naïve B cells in early relapse biopsies was confirmed in independent validation cohorts using scRNA-seq and immunohistochemistry. Notably, naïve B cells in early-relapse samples exhibited high expression of LGALS9, an immunosuppressive molecule known to interact with T cell immunoglobulin and mucin domain 3 (TIM-3) on regulatory T cells (Tregs). Cell-cell interaction analysis revealed the importance of interactions between LGALS9+ naïve B cells and TIM-3+ Tregs in the early-relapse setting. Spatial analysis by imaging mass cytometry confirmed close proximity of LGALS9+ naïve B cells with TIM-3⁺ CD4⁺ T cells and HRS cells, pointing to their role in shaping an immunosuppressive niche. Our findings highlight a previously unrecognized population of LGALS9⁺ naïve B cells with immunoregulatory potential in early-relapse CHL and provide new insights into the spatial and transcriptional architecture of the relapsed TME in CHL.

Study EGAS00001008222

Noninvasive prenatal molecular karyotyping from maternal plasma

Fetal DNA is present in the plasma of pregnant women. Massively parallel sequencing of maternal plasma DNA has been used to detect fetal trisomies 21, 18, 13 and selected sex chromosomal aneuploidies noninvasively. Case reports describing the detection of fetal microdeletions from maternal plasma using massively parallel sequencing have been reported. However, these previous reports were either polymorphism-dependent or used statistical analyses which were confined to one or a small number of selected parts of the genome. In this report, we reported a procedure for performing noninvasive prenatal karyotyping at 3 Mb resolution across the whole genome through the massively parallel sequencing of maternal plasma DNA. This method has been used to analyze the plasma obtained from 6 cases. In 5 cases, fetal microduplications or microdeletions have been detected successfully from maternal plasma. The two cases with fetal microduplications represented the first noninvasive prenatal detection of such changes from maternal plasma. In the remaining case, the plasma DNA sequencing result was consistent with the pregnant mother being a carrier of a microduplication. Simulation analyses were performed for determining the number of plasma DNA molecules that would need to be sequenced and aligned for enhancing the diagnostic resolution of noninvasive prenatal karyotyping to 2 Mb and 1 Mb. In conclusion, noninvasive prenatal molecular karyotyping from maternal plasma by massively parallel sequencing is feasible and would enhance the diagnostic spectrum of noninvasive prenatal testing.

Study EGAS00001000439

Targeting the p53 pathway to treat Malignant Rhabdoid and High-Risk Atypical Teratoid Rhabdoid Tumors

Rhabdoid tumors (RTs) predominantly affect young children and are among the deadliest pediatric solid tumors. Despite multimodal therapy consisting of surgery, radiation, and chemotherapy, children with these tumors have median survival of less than one year. RTs can arise throughout the body, including the central nervous system (CNS) where they are called atypical teratoid rhabdoid tumors (AT/RTs), and in extra-CNS locations such as the kidneys and other soft tissues where they are designated malignant RTs (MRTs). We previously identified MDM2 and MDM4 as therapeutic vulnerabilities in RTs and showed that treatment with the MDM2 inhibitor idasanutlin increased survival in mice bearing MRT xenografts. However, the therapeutic potential of idasanutlin in CNS RT tumors is unknown. Single agent therapies are prone to resistance and show limited clinical benefit on their own. Therefore, we sought to identify combination strategies incorporating idasanutlin that would be effective in both CNS and extra-CNS RTs. Here, we show that the XPO1 inhibitor selinexor increased nuclear retention of p53 and potentiated idasanutlin-induced p53 pathway activation and cytotoxicity in AT/RT and MRT cell lines in vitro. Importantly, combination therapy limited acquired resistance through TP53 mutation. In vivo, combination therapy was well-tolerated, reduced tumor burden, and increased survival in orthotopic models of both AT/RT and MRT. Our results demonstrate that combining idasanutlin with selinexor is a promising therapeutic strategy for children with rhabdoid tumors.

Study EGAS00001007680

Genome-wide analysis of genetic risk factors for rheumatic heart disease in Aboriginal Australians provides support for pathogenic molecular mimicry

Background. Rheumatic heart disease (RHD) following Group A Streptococcus (GAS) infections is heritable and prevalent in Indigenous populations. Molecular mimicry between human and GAS proteins triggers pro-inflammatory cardiac valve-reactive T-cells. Methods. Genome-wide genetic analysis was undertaken in 1263 Aboriginal Australians (398 RHD cases; 865 controls). Single nucleotide polymorphisms (SNPs) were genotyped using Illumina HumanCoreExome BeadChips. Direct typing and imputation was used to fine-map the human leukocyte antigen (HLA) region. Epitope binding affinities were mapped for human cross-reactive GAS proteins, including M5 and M6. Results. The strongest genetic association was intronic to HLA-DQA1 (rs9272622; P=1.86x10-7). Conditional analyses showed rs9272622 and/or DQA1*AA16 account for the HLA signal. HLA-DQA1*0101_DQB1*0503 (OR 1.44, 95%CI 1.09-1.90, P=9.56x10-3) and HLA-DQA1*0103_DQB1*0601 (OR 1.27, 95%CI 1.07-1.52, P=7.15x10-3) were risk haplotypes; HLA_DQA1*0301-DQB1*0402 (OR 0.30, 95%CI 0.14-0.65, P=2.36x10-3) was protective. Human myosin cross-reactive N-terminal and B repeat epitopes of GAS M5/M6 bind with higher affinity to DQA1/DQB1 alpha/beta dimers for the two risk haplotypes than the protective haplotype. Conclusions. Variation at HLA_DQA1-DQB1 is the major genetic risk factor for RHD in Aboriginal Australians studied here. Cross-reactive epitopes bind with higher affinity to alpha/beta dimers formed by risk haplotypes, supporting molecular mimicry as the key mechanism of RHD pathogenesis.

Study EGAS00001002678

Bioinformatic Methods and Bridging of Assay Results for Reliable Tumor Mutational Burden Assessment in Non-Small Cell Lung Cancer

Tumor mutational burden (TMB) has emerged as a predictive biomarker of response to immune checkpoint inhibitors. Standardization of TMB measurement is essential for implementing diagnostic tools to guide treatment. Here we evaluate bioinformatic TMB analysis by whole exome sequencing (WES) in formalin-fixed, paraffin-embedded samples. In CheckMate 026, TMB was retrospectively assessed in 312 patients with non-small cell lung cancer (58% of the intent-to-treat population) who received first-line nivolumab treatment or chemotherapy. We examined the sensitivity of TMB assessment to bioinformatic filtering methods and assessed concordance between TMB data derived by WES and the FoundationOne CDx™ assay. TMB scores comprising synonymous, indel, frameshift, and nonsense mutations (all mutations) were 3.1-fold higher than data including missense mutations only, but values were highly correlated (Spearman’s r = 0.99). Scores from CheckMate 026 samples including missense mutations only were similar to those generated from data in The Cancer Genome Atlas, but those including all mutations were generally higher. Using databases for germline subtraction (instead of matched controls) showed a trend for race-dependent increases in TMB scores. Parameter variation can therefore impact TMB calculations, highlighting the need for standardization. Encouragingly, WES and FoundationOne CDx outputs were highly correlated (Spearman’s r = 0.90) and differences could be accounted for by empirical calibration, suggesting that reliable TMB assessment across assays, platforms and centers is achievable.

Study EGAS00001003661

Transcriptomic profiles of chronic lymphocytic leukemia before and after frontline therapy: 5-year results from the randomized CLL14 study

CLL14 (NCT02242942) is a multinational, open-label Phase III study designed to compare venetoclax-obinutuzumab (Ven-Obi) versus chlorambucil-obinutuzumab (Clb-Obi) in previously untreated patients with CLL and co-existing conditions. Between Aug 7, 2015, and Aug 4, 2016, 432 patients were enrolled and randomly assigned to receive either venetoclax plus obinutuzumab (Ven-Obi, n=216) or chlorambucil plus obinutuzumab (Clb-Obi, n=216). In both study arms, patients received fixed-duration treatment of overall 12 cycles, each cycle 28 days long, and entered subsequent post-treatment surveillance. Peripheral blood samples were collected before the start of therapy, as well as every 3-6 months during and after therapy, including at the time of disease progression. The primary endpoint of the study was progression-free survival, secondary endpoints included overall response rate, minimal residual disease response, and overall survival. The primary read-out was conducted after a median follow-up of 28.1 months, when the primary endpoint was met with 30 primary end-point events (disease progression or death) in the Ven-Obi arm and 77 in the Clb-Obi arm (hazard ratio, 0.35; 95% confidence interval [CI], 0.23 to 0.53;P<0.001). Follow-up is ongoing until 9 years after the last patient enrollment. In an exploratory analysis, bulk RNAseq was conducted with CD19-enriched blood from patients before start of therapy (418 patients) and at first relapse (44 patients).

Study EGAS00001006596

RB1 Gene Inactivation by Chromothripsis in Human Retinoblastoma

Retinoblastoma is a rare childhood cancer of the developing retina. Most retinoblastomas initiate with biallelic inactivation of the RB1 gene through diverse mechanisms including point mutations, nucleotide insertions, deletions, loss of heterozygosity and promoter hypermethylation. Recently, a novel mechanism of retinoblastoma initiation was proposed. Gallie and colleagues discovered that a small proportion of retinoblastomas lack RB1 mutations and had MYCN amplification [1]. In this study, we identified recurrent chromosomal, regional and focal genomic lesions in 94 primary retinoblastomas with their matched normal DNA using SNP 6.0 chips. We also analyzed the RB1 gene mutations and compared the mechanism of RB1 inactivation to the recurrent copy number variations in the retinoblastoma genome. In addition to the previously described focal amplification of MYCN and deletions in RB1 and BCOR, we also identified recurrent focal amplification of OTX2, a transcription factor required for retinal photoreceptor development. We identified 10 retinoblastomas in our cohort that lacked RB1 point mutations or indels. We performed whole genome sequencing on those 10 tumors and their corresponding germline DNA. In one of the tumors, the RB1 gene was unaltered, the MYCN gene was amplified and RB1 protein was expressed in the nuclei of the tumor cells. In addition, several tumors had complex patterns of structural variations and we identified 3 tumors with chromothripsis at the RB1 locus. This is the first report of chromothripsis as a mechanism for RB1 gene inactivation in cancer.

Study EGAS00001000598

Identification_of_synthetic_lethal_genes_by_CRISPR_Cas9_library

Sequencing is performed on PCR products from the amplification of gRNA segments introduced into cell lines (some are commercial, some are derived by collaborators, as described in HMDMC 13/076). Cellular DNA or RNA are not sequenced. This data is defined as "Other Biological/Biochemical Assay Data" under the WTSI datasharing guidelines and consequently will be released at the time of publication and/or through faculty websites. SEQUENCING DETAILS: in October 2016 3 DNA samples will be submitted for PCR amplification and purification through CGP pipeline using Thompson's custom primer set and PCR settings and sequencing with Hiseq using a custom protocol (with 3 forward reads to read sgRNA1, sgRNA2 and the index, optimized with Illumina Bespoke Sanger team); read length 50bp. The following set of samples will be submitted in later 2016 and ealry 2017 as cell pellet to be processed through CGP Cas9 pipeline (Thompson's oligos and protocols) and sequences on Hiseq (same custom sequencing protocol described above) We are utilizing a custom CRISPR/Cas9 library to identify synthetic lethal gene pairs. Cells will be transduced with LV expressing Cas9, then infected with LV expressign the library of 42'000 sgRNA pairs. Cell pellet will be collected from the cell population at different timepoints. The samples will be processed by CGP pipeline customized for the specific protocol and then sequenced to identify the 19bp of the sgRNA and quantify their representation.

Study EGAS00001002117

H3.3G34R/V-transformed interneuron progenitors co-opt PDGFRA for gliomagenesis

Histone H3.3 glycine 34 to arginine/valine (H3.3G34R/V) mutations occur in deadly hemispheric high-grade gliomas. These tumors show exquisite regional and temporal specificity, suggesting a developmental context permissive to the effects of G34R/V mutations. Here we present the molecular landscape of G34R/V gliomas (n=85) and show that 50% bear activating mutations in PDGFRA, with strong selection pressure for PDGFRAMUT clones at recurrence. We show that G34R/V tumors arise in interneuron progenitors of the foetal ventral forebrain expressing GSX2 and the DLX family of homeobox transcription factors, where terminal neuronal differentiation is impaired through aberrant G34R/V-mediated H3K27me3. Frequent co-occurrence of G34R/V & PDGFRAMUT is facilitated in this interneuron lineage-of-origin as PDGFRA forms an aberrant chromatin loop with the adjacent GSX2, hijacking its active chromatin conformation. At the single-cell level, G34R/V tumours entirely lack oligodendroglial transcriptional programs prominent in other glioma entities, and instead harbour dual neuronal and astroglial compartments. CRISPR-removal of H3.3G34R/V does not impact tumorigenicity suggesting this mutation becomes dispensable, while PDGFRAMUT are potently oncogenic regardless of G34 mutation. Collectively, our results suggest that G34R/V gliomas arise in foetal interneuron progenitors unable to undergo terminal differentiation, enabling co-option of PDGFRA through inappropriate expression and activating mutations to promote gliogenesis and oncogenicity. Reliance on PDGFRA for oncogenesis may be of therapeutic opportunity in G34R/V glioma.

Study EGAS00001004301

The spatio-temporal evolution of lymph node spread in early breast cancer

Purpose: The most significant prognostic factor in early breast cancer is lymph node involvement. This stage between localised and systemic disease is key to understanding breast cancer progression, however our knowledge of the evolution of lymph node malignant invasion remains limited, as most currently available data derive from primary tumours. Experimental design: In 11 treatment-naïve node positive early breast cancer patients without clinical evidence of distant metastasis, we investigated lymph node evolution using spatial multi-region sequencing (n=78 samples) of primary and lymph node deposits and genomic profiling of matched longitudinal circulating tumour DNA (ctDNA). Results: Linear evolution from primary to lymph node was rare (1/11) whereas the majority of cases displayed either early divergence between primary and nodes (4/11), or no detectable divergence (6/11) where both primary and nodal cells belonged to a single recent expansion of a metastatic clone. Divergence of metastatic subclones was driven in part by APOBEC. Longitudinal ctDNA samples from 2 subjects taken peri-operatively reflected the two major evolutionary patterns and demonstrate that private mutations can be detected even from early metastatic nodal deposits. Moreover, node removal resulted in disappearance of private lymph node mutations in ctDNA.Conclusions: This study sheds new light on a crucial evolutionary step in the natural history of breast cancer, demonstrating early establishment of axillary lymph node metastasis in a substantial proportion of patients.

Study EGAS00001002947

ICU_transcriptomics__Assessing_the_role_of_the_host_immune_response_in_patients_with_ventilator_associated_pneumonia

The aim of the proposed project is to provide transcriptomic sequencing data analysis of cells isolated from bronchioalveolar lavage and blood samples from critically ill patients with pneumonia on an intensive care unit in order to investigate the host contribution to disease. The lung microbiome of these patients will be analysed in parallel by deep sequencing and 16S sequencing. These data will provide important information to the development of strategies aimed at improving patient outcomes, reducing the development of antimicrobial resistance and facilitating infection control measures in hospital. Specifically, the proposed work will be a proof of principle study of 100 patients. We will use transcriptomic data to define the host contribution to the symptoms of these critically ill patients diagnosed with pneumonia. The data will be correlated with cytokine/chemokine concentrations, routine microbiology analysis and pathogen genomic data from the same samples to determine the ecology of the lungs of these patients. We will look for associations between clinical features, host response and pathogen genomics in order to identify new diagnostic markers. This study should provide important information to improve the speed and accuracy of diagnosis - resulting in more targeted and successful treatment of critically ill patients in the ICU. This data is part of a pre-publication release. For information on the proper use of pre-publication data shared by the Wellcome Trust Sanger Institute (including details of any publication moratoria), please see http://www.sanger.ac.uk/datasharing/

Study EGAS00001003074

Exploiting evolutionary steering in cancer therapy

Drug resistance, mediated by intra-tumour heterogeneity and clonal evolution, is arguably the biggest problem in cancer therapy today. However, evolving resistance to one drug may come at the cost of increased sensitivity to another due to so-called evolutionary trade-offs. This weakness can be exploited in the clinic using an approach called ‘evolutionary herding’ that aims at controlling the tumour cell population to delay or prevent resistance. However, model systems able to recapitulate cancer evolution experimentally are lacking and current in vitro techniques based on small populations, re-plating and escalating dose, are unsuitable to develop evolutionary herding strategies. We present a novel methodology for evolutionary herding in vitro and ex vivo using patient-derived organoids. Our approach is based on a combination of single-cell barcoding, very large populations of 108-109 cells grown without re-plating, realistic high drug doses, time-course monitoring of cancer clones, and mathematical modelling of tumour evolutionary dynamics. We demonstrate evolutionary herding in non-small cell lung cancer in vitro and in patient-derived colorectal cancer organoids (PDO). We show that herding causes controlled evolutionary bottlenecks that lead to collateral sensitivity. Through genomic analysis, we were also able to determine the mechanisms that drive such sensitivity. Our approach allows modelling evolutionary trade-offs experimentally to test patient-specific evolutionary herding strategies that can be translated into the clinic to control treatment resistance.

Study EGAS00001003200

Efficacy of JAK/STAT pathway inhibition in murine xenograft models of early T-cell precursor (ETP) acute lymphoblastic leukemia

Early T-cell precursor (ETP) acute lymphoblastic leukemia (ALL) is a recently described subtype of T-ALL characterized by a unique immunophenotype and genomic profile as well as a high rate of induction failure. Frequent mutations in cytokine receptor and JAK/STAT signaling pathways led us to hypothesize that ETP-ALL is dependent on JAK/STAT signaling. Here we demonstrate aberrant activation of the JAK/STAT pathway in ETP-ALL blasts relative to non-ETP T-ALL. Moreover, ETP-ALL showed hyperactivation of STAT5 in response to IL7, an effect that was abrogated by the JAK1/2 inhibitor ruxolitinib. In vivo, ruxolitinib displayed activity in 6/6 patient-derived murine xenograft models of ETP-ALL, with profound single-agent efficacy in 5 models. Ruxolitinib treatment decreased peripheral blast counts relative to pre-treatment levels and compared to control (P<0.01) in 5/6 ETP-ALL xenografts, with marked reduction in mean splenic blast counts (P<0.01) in 6/6 samples. Surprisingly, both JAK/STAT pathway activation and ruxolitinib efficacy were independent of the presence of JAK/STAT pathway mutations, raising the possibility that the therapeutic potential of ruxolitinib in ETP-ALL extends beyond those cases with JAK mutations. These findings establish the preclinical in vivo efficacy of ruxolitinib in ETP-ALL, a biologically distinct subtype for which novel therapies are needed.

Study EGAS00001001146

RRBS sequencing of 7 tumour regions and a normal sample from a single TRACERx patient.

Background: Epigenetic heterogeneity within a tumour can play an important role in tumour evolution and the emergence of resistance to treatment. It is increasingly recognised that the study of DNA methylation (DNAm) patterns along the genome - so-called `epialleles' - offers greater insight into epigenetic dynamics than conventional analyses which examine DNAm marks individually. Results: We have developed a Bayesian model to infer which epialleles are present in multiple regions of the same tumour. We apply our method to reduced representation bisulfite sequencing (RRBS) data from multiple regions of one lung cancer tumour and a matched normal sample. The model borrows information from all tumour regions to leverage greater statistical power. The total number of epialleles, the epiallele DNAm patterns, and a noise hyperparameter are all automatically inferred from the data. Uncertainty as to which epiallele an observed sequencing read originated from is explicitly incorporated by marginalising over the appropriate posterior densities. The degree to which tumour samples are contaminated with normal tissue can be estimated and corrected for. By tracing the distribution of epialleles throughout the tumour we can infer the phylogenetic history of the tumour, identify epialleles that differ between normal and cancer tissue, and define a measure of global epigenetic disorder. Conclusions: Detection and comparison of epialleles within multiple tumour regions enables phylogenetic analyses, identification of differentially expressed epialleles, and provides a measure of epigenetic heterogeneity.

Study EGAS00001002484

Combining transcription factor binding affinities with open chromatin data for accurate gene expression prediction

The binding and contribution of transcription factors (TF) to cell specific gene expression is often deduced from open-chromatin measurements to avoid cost and labour intensive TF ChIP-seq assays.It is important to develop reliable and fast computational methods for accurate TF binding prediction in open-chromatin regions (OCRs). Here, we report a novel segmentation-based method, TEPIC, to predict TF binding by combining sets of OCRs with position weight matrices.TEPIC can be applied to various open-chromatin data, e.g. DNaseI-seq and NOMe-seq, using either peaks or footprints as input.In addition to open-chromatin data, also Histone-Marks (HMs) can be used in TEPIC to identify candidate TF binding sites.TEPIC computes TF affinities and uses open-chromatin/HM signal intensity as quantitative measures of TF binding strength.Using machine learning techniques, we show that incorporating low affinity binding sites improves our ability to explain gene expression variability compared to the standard presence/absence classification of binding sites.Further, we show that both footprints and peaks capture essential TF binding events and lead to a good prediction performance.In our application, gene-based scores computed by TEPIC with one open-chromatin assay nearly reach the quality of several TF ChIP-seq datasets.Finally, we show that these scores correctly predict known transcriptional regulators as illustrated by the application to novel DNaseI-seq and NOMe-seq data for primary human hepatocytes and CD4+ T-cells, respectively.

Study EGAS00001002073

Atypical 3q26/MECOM rearrangements genocopy inv(3)/t(3;3) in acute myeloid leukemia

Acute myeloid leukemia (AML) with inv(3)/t(3;3)(q21q26) is a distinct WHO recognized entity, characterized by its aggressive course and poor prognosis. In this subtype of AML, the translocation of a GATA2 enhancer (3q21) to MECOM (3q26) results in overexpression of the MECOM isoform EVI1 and monoallelic expression of GATA2 from the unaffected allele. The full-length MECOM transcript, MDS1-EVI1, is not expressed as the result of the 3q26 rearrangement. Besides the classical inv(3)/t(3;3), a number of other 3q26/MECOM rearrangements with poor treatment response have been reported in AML. Here we demonstrate, in a group of 33 AML patients with atypical 3q26 rearrangements, MECOM involvement with EVI1 overexpression, but no or low MDS1-EVI1 levels. Moreover, the 3q26 translocations in these AML patients often involve super-enhancers of genes active in myeloid development (e.g. CD164, PROM1, CDK6 or MYC). In more than 50% of these cases allele specific GATA2 expression was observed, either by copy number loss or by an unexplained allelic imbalance. Altogether, atypical 3q26 recapitulate the main leukemic mechanism of inv(3)/t(3;3) AML, namely EVI1 overexpression driven by enhancer hijacking, absent MDS1-EVI1 expression and potential GATA2 involvement. Therefore, we conclude that both atypical 3q26/MECOM and inv(3)/t(3;3) can be classified as a single entity of 3q26-rearranged AMLs. Routine analyses determining MECOM rearrangements, EVI1 and MDS1-EVI1 expression are required to recognize 3q-rearranged AML cases.

Study EGAS00001004325

Partner-independent fusion gene detection by multiplexed CRISPR/Cas9 enrichment and long-read Nanopore sequencing

Fusion genes are hallmarks of various cancer types and important determinants for diagnosis, prognosis and treatment possibilities. The promiscuity of fusion genes with respect to partner choice and exact breakpoint-positions restricts their detection in the diagnostic setting, even for known and recurrent fusion gene configurations. To accurately identify these gene fusions in an unbiased manner, we developed FUDGE: a FUsion gene Detection assay from Gene Enrichment. FUDGE couples target-selected and strand-specific CRISPR/Cas9 activity for enrichment and detection of fusion gene drivers (e.g. BRAF, EWSR1, KMT2A/MLL) - without prior knowledge of fusion partner or breakpoint-location - to long-read Nanopore sequencing. FUDGE encompasses a dedicated bioinformatics approach (NanoFG) to detect fusion genes from Nanopore sequencing data. Our strategy is flexible with respect to target choice and enables multiplexed enrichment for simultaneous analysis of several genes in multiple samples in a single sequencing run. We observe on average a 508-fold on-target enrichment and identify fusion breakpoints at nucleotide resolution - all within two days. We demonstrate that FUDGE effectively identifies fusion genes in cancer cell lines, tumor samples and on whole genome amplified DNA irrespective of partner gene or breakpoint-position in 100% of cases. Furthermore, we show that FUDGE is superior to routine diagnostic methods for fusion gene detection. In summary, we have developed a rapid and versatile fusion gene detection assay, providing an unparalleled opportunity for pan-cancer detection of fusion genes in routine diagnostics.

Study EGAS00001003964

Identification of the mutational consequences of precancerous liver disease (including alcohol abuse) on the genomes of human adult stem cells.

Stem cells are believed to be the cells-of-origin cancer. Nevertheless, their genome remains remarkably stable during life. To increase insight into how these cells can drive tumorigenesis, we set out to determine the genomic integrity of human adult stem cells from patients with liver disease associated with increased cancer risk. To this end, we characterized genome-wide mutation profiles of single ASCs from livers of patients with alcoholic liver disease, non-alcoholic liver disease and primary sclerosing cholangitis. The patients included in our study received a liver transplant due to their liver disease. Strikingly, we found that ASCs in these livers maintain a stable genome and have comparable mutation rates as healthy human liver ASCs. We did observe more nonsense and nonsynonymous mutations in COSMIC cancer genes that are hit by dominant mutations in tumors in comparison to healthy ASCs (n.s.), including two ASCs in alcoholic liver disease patients with a nonsense mutation in PTPRK. Next, we sequenced 5 biopsies of an alcoholic tumour to identify the mutational patterns in the tumour-initiating ASC. We found that the mutational patterns are, again, similar to healthy liver, but the mutational load is slightly higher. This suggests that the increased cancer risk in the liver disease patients is not caused by an increased mutational load, but rather that the precancerous liver environment promotes outgrowth of cells that would normally be selected against.

Study EGAS00001002983

Small variants in mtDNA Canary Islands - WES Illumina (ITER)

The mitochondrial genome (mtDNA) is of interest for a range of fields including evolutionary, forensic, and medical genetics. Human mitogenomes can be classified into evolutionary related haplogroups that provide ancestral information and pedigree relationships. Because of this and the advent of high-throughput sequencing (HTS) technology, there is a diversity of bioinformatic tools for haplogroup classification. We present a benchmarking of the 11 most salient tools for human mtDNA classification using empirical whole-genome (WGS) and whole-exome (WES) short-read sequencing data from 36 unrelated donors. Besides, because of its relevance, we also assess the best performing tool in third-generation long noisy read WGS data obtained with nanopore technology for a subset of the donors. We found that, for short-read WGS, most of the tools exhibit high accuracy for haplogroup classification irrespective of the input file used for the analysis. However, for short- read WES, Haplocheck and MixEmt were the most accurate tools. Based on the performance shown for WGS and WES, and the accompanying qualitative assessment, Haplocheck stands out as the most complete tool. For third-generation HTS data, we also showed that Haplocheck was able to accurately retrieve mtDNA haplogroups for all samples assessed, although only after following assembly-based approaches (either based on a referenced-based assembly or a hybrid de novo assembly). Taken together, our results provide guidance for researchers to select the most suitable tool to conduct the mtDNA analyses from HTS data.

Study EGAS00001005678

Small variants in mtDNA Canary Islands - WGS Illumina (ITER)

The mitochondrial genome (mtDNA) is of interest for a range of fields including evolutionary, forensic, and medical genetics. Human mitogenomes can be classified into evolutionary related haplogroups that provide ancestral information and pedigree relationships. Because of this and the advent of high-throughput sequencing (HTS) technology, there is a diversity of bioinformatic tools for haplogroup classification. We present a benchmarking of the 11 most salient tools for human mtDNA classification using empirical whole-genome (WGS) and whole-exome (WES) short-read sequencing data from 36 unrelated donors. Besides, because of its relevance, we also assess the best performing tool in third-generation long noisy read WGS data obtained with nanopore technology for a subset of the donors. We found that, for short-read WGS, most of the tools exhibit high accuracy for haplogroup classification irrespective of the input file used for the analysis. However, for short- read WES, Haplocheck and MixEmt were the most accurate tools. Based on the performance shown for WGS and WES, and the accompanying qualitative assessment, Haplocheck stands out as the most complete tool. For third-generation HTS data, we also showed that Haplocheck was able to accurately retrieve mtDNA haplogroups for all samples assessed, although only after following assembly-based approaches (either based on a referenced-based assembly or a hybrid de novo assembly). Taken together, our results provide guidance for researchers to select the most suitable tool to conduct the mtDNA analyses from HTS data.

Study EGAS00001005679

DNA methylation analysis on PBL obtained from male patients with UC-PSC, UC and HC

Primary sclerosing cholangitis (PSC) is a chronic inflammatory liver disease affecting the intra- and extrahepatic bile ducts, and is strongly associated with ulcerative colitis (UC). In this study, we explored the peripheral blood DNA methylome and its immune cell composition in patients with PSC-UC, UC, and healthy controls (HC) with the aim to develop a predictive assay in distinguishing patients with PSC-UC from those with UC alone. The peripheral blood DNA methylome of male patients with PSC and concomitant UC, UC and HCs was profiled using the Illumina HumanMethylation Infinium EPIC BeadChip (850K) array. Differentially methylated CpG position (DMP) and region (DMR) analyses were performed alongside gradient boosting classification analyses to discern PSC-UC from UC patients. As observed differences in the DNA methylome could be the result of differences in cellular populations, we additionally employed mass cytometry to characterize the immune cell compositions. Genome wide methylation analysis revealed no differentially methylated positions between PSC-UC and UC patients. Nonetheless, using gradient boosting we were capable of discerning PSC-UC from UC with an area under the receiver operator curve (AUROC) of 0.80. Four CpG sites annotated to the NINJ2 gene were found to strongly contribute to the predictive performance. While mass cytometry analyses corroborated the largely similar blood cell composition among patients with PSC-UC, UC and HC, a higher abundance of myeloid cells was observed in UC compared to PSC-UC patients.

Study EGAS00001005832

A novel orthotopic patient-derived xenograft model of radiation-induced glioma following medulloblastoma

Radiation-induced glioma (RIG) is a highly aggressive brain cancer arising as a consequence of radiation therapy. We report a case of RIG that arose in the brain stem following treatment for paediatric medulloblastoma, and the development and characterisation of a matched orthotopic patient-derived xenograft (PDX) model (TK-RIG915). Patient and PDX tumours were analysed using DNA methylation profiling, whole genome sequencing (WGS) and RNA sequencing. While initially thought to be a diffuse intrinsic pontine glioma (DIPG) based on disease location, results from methylation profiling and WGS were not consistent with this diagnosis. Furthermore, clustering analyses based on RNA expression suggested the tumours were distinct from primary DIPG. Additional gene expression analysis demonstrated concordance with a published RIG expression profile. Multiple genetic alterations that enhance PI3K/AKT and Ras/Raf/MEK/ERK signalling were discovered in TK-RIG915 including an activating mutation in PIK3CA, upregulation of PDGFRA and AKT2, inactivating mutations in NF1, and a gain-of-function mutation in PTPN11. Additionally, deletion of CDKN2A/B, increased IDH1 expression, and decreased ARID1A expression were observed. Detection of phospho-S6, -4EBP1 and -ERK via immunohistochemistry confirmed PI3K pathway and ERK activation. Here, we report one of the first PDX models for RIG, which recapitulates the patient disease and is molecularly distinct from primary brain stem glioma. Genetic interrogation of this model has enabled the identification of potential therapeutic vulnerabilities in this currently incurable disease.

Study EGAS00001004709

Small variants in mtDNA Canary Islands - WGS Oxford Nanopore Technologies (ITER)

The mitochondrial genome (mtDNA) is of interest for a range of fields including evolutionary, forensic, and medical genetics. Human mitogenomes can be classified into evolutionary related haplogroups that provide ancestral information and pedigree relationships. Because of this and the advent of high-throughput sequencing (HTS) technology, there is a diversity of bioinformatic tools for haplogroup classification. We present a benchmarking of the 11 most salient tools for human mtDNA classification using empirical whole-genome (WGS) and whole-exome (WES) short-read sequencing data from 36 unrelated donors. Besides, because of its relevance, we also assess the best performing tool in third-generation long noisy read WGS data obtained with nanopore technology for a subset of the donors. We found that, for short-read WGS, most of the tools exhibit high accuracy for haplogroup classification irrespective of the input file used for the analysis. However, for short- read WES, Haplocheck and MixEmt were the most accurate tools. Based on the performance shown for WGS and WES, and the accompanying qualitative assessment, Haplocheck stands out as the most complete tool. For third-generation HTS data, we also showed that Haplocheck was able to accurately retrieve mtDNA haplogroups for all samples assessed, although only after following assembly-based approaches (either based on a referenced-based assembly or a hybrid de novo assembly). Taken together, our results provide guidance for researchers to select the most suitable tool to conduct the mtDNA analyses from HTS data.

Study EGAS00001005677

Single-cell ATAC-seq analysis for COVID19 patients

While SARS-CoV-2 infection causes mild respiratory disease in most individuals, a small group of patients develops severe COVID-19. Dysfunctional innate immune responses have been identified to contribute to differences in COVID-19 severity, but the key regulators are still unknown. Here, we present an integrative single-cell epigenetics, transcriptomic, and genetics analysis of peripheral blood mononuclear cells from hospitalized and convalescent COVID-19 patients. In classical monocytes, we identified 41.3% of significantly up-regulated genes in hospitalized COVID-19 patients potentially induced by differential chromatin accessibility. Sub-clustering and motif-enrichment analyses of monocytes reveal disease condition-specific regulation by transcription factors, such as C/EBPs and SPI1, and their targets, including a long-noncoding RNA LUCAT1, which further regulates interferon responses and is associated with the need for oxygen supply of COVID-19 patients. The interaction between C/EBPs and LUCAT1 was validated through loss-of-function experiments. Finally, we investigated genetic risk variants that exhibit allele-specific open chromatin (ASoC) in promoters/enhancers of COVID-19 patients. Integrating our data with publicly available expression quantitative trait loci and chromosomal interactions indicates that ASoC SNP rs6800484-C is associated with lower expression of CCR2, which may contributeto higher viral loads in lungs and higher risk of COVID-19 hospitalization. Altogether, our study highlights the diverse genetic and epigenetic regulators that contribute to the innate immune responses of different COVID-19 patients.

Study EGAS00001006559

Single-cell RNA-seq analysis for COVID19 patients

While SARS-CoV-2 infection causes mild respiratory disease in most individuals, a small group of patients develops severe COVID-19. Dysfunctional innate immune responses have been identified to contribute to differences in COVID-19 severity, but the key regulators are still unknown. Here, we present an integrative single-cell epigenetics, transcriptomic, and genetics analysis of peripheral blood mononuclear cells from hospitalized and convalescent COVID-19 patients. In classical monocytes, we identified 41.3% of significantly up-regulated genes in hospitalized COVID-19 patients potentially induced by differential chromatin accessibility. Sub-clustering and motif-enrichment analyses of monocytes reveal disease condition-specific regulation by transcription factors, such as C/EBPs and SPI1, and their targets, including a long-noncoding RNA LUCAT1, which further regulates interferon responses and is associated with the need for oxygen supply of COVID-19 patients. The interaction between C/EBPs and LUCAT1 was validated through loss-of-function experiments. Finally, we investigated genetic risk variants that exhibit allele-specific open chromatin (ASoC) in promoters/enhancers of COVID-19 patients. Integrating our data with publicly available expression quantitative trait loci and chromosomal interactions indicates that ASoC SNP rs6800484-C is associated with lower expression of CCR2, which may contributeto higher viral loads in lungs and higher risk of COVID-19 hospitalization. Altogether, our study highlights the diverse genetic and epigenetic regulators that contribute to the innate immune responses of different COVID-19 patients.

Study EGAS00001006560

Atypical B cells and impaired SARS-CoV-2 neutralisation following heterologous vaccination in the elderly

Age is a major risk factor for hospitalization and death after SARS-CoV-2 infection, even in vaccinees. Suboptimal responses to a primary vaccination course have been reported in the elderly, but there is little information regarding the impact of age on responses to booster third doses. Here we show that individuals 70 or older who received a primary two dose schedule with AZD1222 and booster third dose with mRNA vaccine achieved significantly lower neutralizing antibody responses against SARS-CoV-2 spike pseudotyped virus compared to those below 70 at one month post booster. Neither the concentration of serum binding anti spike IgG antibody, nor the frequency of spike- specific B cells showed differences by age grouping. Impaired neutralization potency and breadth post-third dose in the elderly was associated with enrichment of circulating “atypical” spike-specific B cells expressing CD11c and FCRL5. Presence of these atypical B cells was confirmed by scRNA seq. However, when considering individuals who received three doses of mRNA vaccine, we did not observe statistically significant differences in neutralization or enrichment in atypical B cells by age grouping following dose 3. Together, our data reveal impaired neutralising antibody responses in the elderly after heterologous vaccination with AdV followed by mRNA booster dose, but not in homologous three dose mRNA vaccination highlighting that different vaccine formats differentially instruct the memory B cell response.

Study EGAS00001007385

DHODH is a synthetic-lethal vulnerability in BCOR-deficient acute myeloid leukemia

Acute Myeloid Leukemia (AML) remains challenging to treat, especially in cases with mutations in the BCL-6 co- repressor (BCOR), which are associated with poor prognosis and chemo-resistance. In this study, we reveal a synthetic lethal interaction between BCOR and dihydroorotate dehydrogenase (DHODH). We demonstrate that BCOR-deficient cells have a heightened sensitivity to DHODH inhibitors such as brequinar and leflunomide, that are already in clinical use. We confirm that DHODH inhibition selectively induces cell death in BCOR-mutant cells in multiple cellular models, in malignant and non-malignant cells, through chemical and genetic manipulation. Interestingly, we find that the dependency on DHODH does not stem from its role in de novo pyrimidine biosynthesis disruption. Rather, DHODH’s role in the electron transport chain, essential for mitigating reactive oxygen species, may be the physiological vulnerability that pushes BCOR-mutant cells toward cell death when DHODH is inhibited. DHODH inhibitors could be repurposed as targeted therapies for BCOR-mutant tumors, offering a promising strategy for precision medicine in AML and other cancers.

Study EGAS50000001108

6170 results for "fc 26 ps5 coins Buyfc26coins.com is FC 26 coins official site..nqFO"

in 24.97 milliseconds.