3334 results for "*geneti*"

in 90.99 milliseconds.

• Comparison of the diagnostic yield of aCGH and NGS across different neurodevelopmental disorders

  Current recommendations for the genetic diagnosis of neurodevelopmental disorders (NDDs) do not include the use of next generation sequencing (NGS) and are still based on chromosomal microarrays, such as comparative genomic hybridization array (aCGH). The present study compares the diagnostic yield obtained by aCGH and clinical exome sequencing in NDD globally and its spectrum of disorders. To that end, 1412 patients clinically diagnosed with NDDs and studied with aCGH were classified into phenotype categories: Global developmental delay/intellectual disability (GDD/ID); autism spectrum disorder (ASD); and Other NDDs. These categories were further subclassified into subcategories, based on the most frequent accompanying signs and symptoms: isolated forms, forms with epilepsy; forms with micro or macrocephaly and syndromic forms. A subset of 245 patients of the 1412 were subjected to clinical exome sequencing. Diagnostic yield of aCGH and clinical exome sequencing, expressed as the number of solved cases, was compared for each phenotype category and subcategory. Clinical exome sequencing was superior than aCGH for all cases except for isolated ASD, with no cases solved by NGS. Globally, clinical exome sequencing solved 20% of cases (versus 5.7% by aCGH) and the diagnostic yield was highest for all forms of GDD/ID and lowest for Other NDDs (7.1% versus 1.4% by aCGH) and ASD (6.1% versus 3% by aCGH). In the majority of cases, diagnostic yield was higher in the phenotype subcategories than in the mother category. These results support the use of NGS as the first-tier test in the diagnostic algorithm of all NDDs followed by aCGH when necessary.

  Study EGAS00001004949

• Cell-free DNA TAPS for early cancer detection

  Multimodal, genome-wide characterization of epigenetic and genetic information in circulating cell-free DNA (cfDNA) could enable more sensitive early cancer detection. However, due to technological challenges associated with DNA methylation sequencing in low input cfDNA samples, most studies have been limited by DNA damage caused by bisulfite sequencing, or the qualitative nature of enrichment-based sequencing. Recently, we developed TET-assisted Pyridine Borane Sequencing (TAPS), which is a mild, bisulfite-free method for base-resolution direct DNA methylation sequencing. Here we optimized TAPS for cfDNA (cfTAPS) to provide high-quality and high-depth whole-genome cell-free methylomes. We applied cfTAPS to 85 cfDNA samples from patients with hepatocellular carcinoma (HCC) or pancreatic ductal adenocarcinoma (PDAC) and non-cancer controls. From just 10 ng cfDNA (1-3 mL of plasma), we generated the most comprehensive cfDNA methylome to date. We demonstrated that cfTAPS provides multimodal information about cfDNA characteristics, including DNA methylation, tissue of origin, and DNA fragmentation. Integrated analysis of these epigenetic and genetic features enables accurate identification of early HCC and PDAC.

  Study EGAS00001004962

• Genetic heterogeneity and dynamics of transcriptional subtypes in matched primary and recurrent head and neck squamous cell carcinomas

  To get a better understanding of different HNSCC subtypes, whole-exome (matched tumor-normal paits) and whole transcriptome sequencing to to elucidate underlying molecular processes. Here, we present the matched whole-exome sequencing data from 28 patients (normal tissue, tumor tissue, relapse).

  Study EGAS00001005005

• DO NOT USE - Whole genome sequencing of SI-NETs from six patients

  Small intestine neuroendocrine tumor (SI-NET), the most common cancer of the small bowel, often displays a curious multifocal phenotype with several intestinal tumors centered around a regional lymph node metastasis, yet the typical path of evolution of these lesions remains unclear. Here, we determined the complete genome sequences of 24 tumor and 3 adjacent normal tissue samples with their paired normal blood samples (totally 33 whole genomes) from 6 patients with multifocal SI-NETs, allowing elucidation of phylogenetic relationships between multiple intestinal tumors and metastases in individual patients.

  Study EGAS00001005009

• Deep genetic affinity between coastal Pacific and Amazonian natives evidenced by Australasian ancestry

  Analysis of the genetic affinity between coastal Pacific and Amazonian natives (i.e. "Y population" ancestry). Here we analyzed a dataset of newly generated data from four Brazilian indigenous communities genotyped in the Axiom Human Origins array combined with publicly available data from other Native American populations.

  Study EGAS00001005022

• The WID-EC test for the detection and risk prediction of endometrial cancer

  Study EGAS00001005033

• The DNA methylome of cervical cells and risk of ovarian cancer

  Study EGAS00001005045

• Genotype data from 'Evidence of the interplay of genetics and culture in Ethiopia'

  Study EGAS00001005171

• Genetic analysis of HLA and immune escape genes in Diffuse Large B-cell Lymphoma

  Fifty percent of diffuse large B-cell lymphoma (DLBCL) cases lack cell-surface expression of the class-I major histocompatibility complex (MHC-I), thereby escaping immune recognition by cytotoxic T cells. In order to comprehensively identify the mechanisms involved in MHC-I loss, we first performed immunophenotypic analysis of both MHC class-I and -II in 657 cases across the spectrum of B-NHL, revealing that loss of MHC-I, but not MHC-II, is preferentially restricted to DLBCL. We then used whole exome and targeted deep-sequencing to examine genes involved in MHC-I expression in 74 DLBCL samples representative of MHC-I positive and negative cases. We show here that somatic biallelic or monoallelic inactivation of B2M and/or HLA-I is present in 80.9% (34/42) of MHC-I negative tumors. Furthermore, 68.8% (22/32) of MHC-I positive DLBCLs also harbored monoallelic HLA-I genetic alterations (MHC-I positivemono) that lead to allelic imbalance, suggesting allele-specific inactivation. Both MHC-I negative and MHC-I positivemono cases showed a significantly higher mutational burden as well as a higher number of inferred neo-antigens, suggesting co-selection of HLA-I loss and sustained neo-antigen production. Interestingly, the analysis of > 500.000 individuals in two databases revealed homozygosity of germline HLA-I genes in 26-38% of DLBCL patients, a frequency significantly higher than that observed in any other cancer type. In mice, germinal-center B cells lacking HLA-I expression did not progress to lymphoma and were counterselected in the context of oncogene-driven lymphomagenesis, suggesting that additional events are needed to license immune evasion. These results suggest a multi-step process of HLA-I loss including both germ-line and somatic events in DLBCL development, and have direct implications for the pathogenesis and immunotherapeutic targeting of this disease.

  Study EGAS00001005054

• Mutational analysis of an oligoprogressive sarcomatoid hepatocellular carcinoma treated with an immune checkpoint inhibitor.

  In this study, we describe a patient diagnosed with metastatic sarcomatoid hepatocellular carcinoma that showed a remarkable tumor response to Nivolumab, an immune checkpoint inhibitor, followed by an oligoprogression. Molecular analysis of the primary tumor in addition to a bone and brain metastasis was performed and we were able to identify driver genetic alterations indicating a possible role in the tumorigenesis process besides some alterations that may be associated with metastasis.

  Study EGAS00001005064

• Multiple migrations to the Philippines during the last 50,000 years

  A key link to understand human history in Island Southeast Asia is the Philippine archipelago and its poorly investigated genetic diversity. We collected and analyzed the most comprehensive set of population-genomic data for the Philippines, 1028 individuals covering 115 indigenous communities. We demonstrate that the Philippines were populated by at least five waves of human migration. The Cordillerans migrated into the Philippines prior to the arrival of rice agriculture, where some remain to be the least admixed East Asians carrying an ancestry shared by all Austronesian-speaking populations; thereby challenging an exclusive Out-of-Taiwan model of joint farming-language-people dispersal. Altogether, our findings portray Philippines as a crucial gateway, with a multilayered history, that ultimately changed the genetic landscape of the Asia-Pacific region.

  Study EGAS00001005083

• Genetic landscape of inherited retinal dystrophies affected cases in Spain

  Inherited retinal diseases (IRDs), defined by dysfunction or progressive loss of photoreceptors, are disorders characterized by elevated heterogeneity, both at the clinical and genetic levels. Our main goal was to address the genetic landscape of IRD in a Spanish cohort of patients. A retrospective hospital-based cross-sectional study was carried out on IRD affected individuals, referred for genetic testing from all the Spanish autonomous communities. Clinical, demographic and familiar data were collected from each patient, including family pedigree, age of appearance of visual symptoms, presence of any systemic findings and geographical origin. Genetic studies were performed to all families with available DNA using different molecular techniques. Our study provides the general landscape for IRD in Spain and our results have important implications for genetic diagnosis, counselling and new therapeutic strategies to both the Spanish population and other related populations.

  Study EGAS00001005104

• Epigenomic profiling and transcriptomic profiling of LUAD patients tumor tissues and tumor adjacent tissues

  To provide an in-depth understanding of the epigenomic heterogeneity of LUAD, we here investigated the H3K27ac histone modification profiles of tumors and adjacent normal lung tissues from 42 LUAD patients and explored the role of epigenetic alterations in LUAD progression. We also investigated transcriptomic alterations among 36 patients tumor tissues.

  Study EGAS00001005132

• Nasal DNA methylation at three CpG sites predicts childhood allergic disease

  Study EGAS00001005189

• mRNA capture sequencing and RT-qPCR for the detection of pathognomonic, novel and secondary fusion transcripts in formalin-fixed paraffin-embedded tissue: a sarcoma showcase

  RNA sequencing is widely used to study gene abundance patterns, but less frequently to assess structural alterations such as fusion transcripts. In this study, we assess the performance of messenger RNA (mRNA) capture sequencing to identify fusion transcripts in formalin-fixed paraffin-embedded (FFPE) tissue of different types of sarcoma, followed by orthogonal confirmation using RT-qPCR. To validate our workflow, we used the TruSight RNA Pan-Cancer Panel to analyze a first cohort of 6 positive control tumors for which the diagnostic workup demonstrated the presence of a specific chromosomal rearrangement. Fusion transcript calling by FusionCatcher confirmed all known aberrations and enabled the identification of both fusion gene partners and breakpoints. Next, whole transcriptome TruSeq RNA Exome sequencing was applied to a second cohort of 17 patients, diagnosed with fusion gene-negative alveolar rhabdomyosarcoma (ARMS) or undifferentiated round cell sarcoma (URCS), for whom fluorescence in situ hybridization (FISH) did not identify the classical pathognomonic rearrangements. Remarkably, for 6 out of 17 patients, mRNA capture sequencing readily detected a pathognomonic fusion transcript, i.e. PAX3-FOXO1 in 2 ARMS patients, and EWSR1-FLI1, EWSR1-ERG or EWSR1-NFATC2 in 4 URCS patients. These data show that RNA capture sequencing may enhance the detection rate of pathognomonic fusion genes in sarcoma. For the 11 patients with no known pathognomonic fusion transcript, 11 newly identified fusion transcripts were confirmed by RT-qPCR, including COPS3-TOM1L2, NCOA1-DTNB, WWTR1-LINC01986, PLAA-MOB3B, AP1B1-CHEK2 and BRD4-LEUTX fusion transcripts in ARMS patients. Additionally, recurrently detected secondary fusion transcripts in two patients diagnosed with EWSR1-NFATC2-positive sarcoma were confirmed (COPS4-TBC1D9, PICALM-SYTL2, SMG6-VPS53 and UBE2F-ALS2). Our results provide new insights into the underlying genetic causes of these malignancies.

  Study EGAS00001005202

• Novel optineurin frameshift insertion causing familial frontotemporal dementia and parkinsonism without amyotrophic lateral sclerosis

  Frontotemporal Dementia (FTD) is a common cause of Young Onset Dementia and has diverse clinical manifestations involving behaviour, executive function, language and motor function, including parkinsonism. Up to 50% of FTD patients report a positive family history, supporting a strong genetic basis, particularly in cases with both FTD and amyotrophic lateral sclerosis (FTD-ALS). Mutations in three genes are associated with the majority of familial FTD (fFTD) cases - microtubule associated protein tau gene (MAPT), progranulin gene (PGRN), and hexanucleotide repeat expansions in chromosome 9 open reading frame 72 (C9orf72), whilst mutations in other genes such as optineurin (OPTN) have rarely been reported. Mutations in OPTN have been reported mostly in familial and sporadic cases of ALS, or in rare cases of FTD-ALS, but not in association with pure or predominant FTD and/or parkinsonian phenotype. Here, we report for the first time, a family from the Philippines with 4 affected members harbouring the same OPTN frameshift insertion, presenting with FTD-related phenotypes, including one sibling with predominant parkinsonism resembling corticobasal syndrome. Notably, none of the affected members showed any evidence of motor neuron disease or ALS at the time of writing, both clinically and on electrophysiological testing, expanding the phenotypic spectrum of OPTN mutations. Close follow-up of mutation carriers for the development of new clinical features and wider investigation of additional family members with further genetic analyses will be conducted to investigate the possibility of other genetic modifiers in this family which could explain phenotypic heterogeneity.

  Study EGAS00001005220

• Targeted sequencing of follicular lymphoma tumour samples from the UK's Haematological Malignancy Research Network

  Follicular lymphoma displays morphological and clinical heterogeneity, underpinned by shared molecular pathogenesis. Mutations in epigenetic regulators alongside t(14;18) have been identified as disease initiating events. Additional mutational diversity suggests disease heterogeneity. Whether this resolves into mechanistically distinct subsets remains an open question. Targeted sequencing was applied to an unselected population-based FL cohort (n=548) diagnosed in a catchment population of ~4 million (14 centers), with full clinical follow-up (n=538) that included 96 DLBCL transformations. We investigated whether molecular subclusters of FL can be identified, and whether mutational data provide predictive information relating to transformation.

  Study EGAS00001005238

• Radiation-induced gliomas represent H3-/IDH-wild type pediatric gliomas with recurrent PDGFRA amplification and loss of CDKN2A/B

  Long-term complications such as radiation-induced second malignancies occur in a subset of patients following radiation-therapy, particularly relevant in pediatric patients due to the long follow-up period in case of survival. Radiation-induced gliomas (RIGs) have been reported in patients after treatment with cranial irradiation for various primary malignancies such as acute lymphoblastic leukemia (ALL) and medulloblastoma (MB). We perform comprehensive (epi-) genetic and expression profiling of RIGs arising after cranial irradiation for MB (n=23) and ALL (n=9). Our study reveals a unifying molecular signature for the majority of RIGs, with recurrent PDGFRA amplification and loss of CDKN2A/B and an absence of somatic hotspot mutations in genes encoding histone 3 variants or IDH1/2, uncovering diagnostic markers and potentially actionable targets.

  Study EGAS00001005243

• Molecular programs of melanoma brain metastases (MBM)

  We performed paired-end whole transcriptome sequencing of fresh frozen or intraoperative brain metastases of therapeutically treated stage IV melanoma patients. In addition, we analysed the epigenetic landscape of tumors usind Illumina 850k bead chip technology.

  Study EGAS00001005976

• Multimodal single-cell and bulk glioma analyses

  Glioma intratumoral heterogeneity enables adaptation to challenging microenvironments and contributes to therapeutic resistance. We integrated 914 single-cell DNA methylomes, 55,284 single-cell transcriptomes, and bulk multi-omic profiles across 11 adult IDH-mutant or IDH-wild-type gliomas to delineate sources of intratumoral heterogeneity. We show that local DNA methylation disorder associates with cell-to-cell DNA methylation differences, is elevated in more aggressive tumors, links with transcriptional disruption, and is altered in environmental stress response. Glioma cells under in vitro hypoxic and irradiation stress increased local DNA methylation disorder and shifted cell states. We identified a positive association between genetic and epigenetic instability that was supported in bulk longitudinally collected DNA methylation data. Increased DNA methylation disorder associated with accelerated disease progression, and recurrently selected DNA methylation changes were enriched for environmental stress response pathways. Our work identifies an epigenetically facilitated adaptive stress response process and highlights the importance of epigenetic heterogeneity in shaping therapeutic outcomes.

  Study EGAS00001005300

• COVID19 Host Genetic Initiative

  The COVID-19 host genetics initiative is a bottom-up collaborative effort that has three main goals: 1. Provide an environment to foster the sharing of resources to facilitate COVID-19 host genetics research (e.g. protocols, questionnaires). 2. Organize analytical activities across studies to identify genetic determinants of COVID-19 susceptibility and severity. 3. Provide a platform to share the results from such activities, as well as the individual-level data where possible, to benefit the broader scientific community. Here, we release the harmonized individual-level data for the subset of individuals who participated in COVID-19 HGI. Please refer to this data dictionary: https://docs.google.com/spreadsheets/d/1whu_MjShXxbdoAJMCYxZ5C-vWOcRivwfWglTcBRNZgQ/edit#gid=470747846 Please cite both Nakanishi, T el al. J Clin Invest 2021 (PMID: 34597274) and COVID-19 Host Genetics Initiative, Nature 2021 (PMID: 34237774) when you use these data for your publications.

  Study EGAS00001005304

• Bulk-tissue RNA-sequencing of anterior cingulate cortex samples derived from Lewy body disease patients

  Parkinson’s disease (PD), Parkinson’s disease with dementia (PDD) and dementia with Lewy bodies (DLB) are three clinically, genetically and neuropathologically overlapping neurodegenerative diseases collectively known as the Lewy body diseases (LBDs). A variety of molecular mechanisms have been implicated in PD pathogenesis, but the mechanisms underlying PDD and DLB remain largely unknown, a knowledge gap that presents an impediment to the discovery of disease-modifying therapies. Transcriptomic profiling can contribute to addressing this gap, but remains limited in the LBDs. Here, we applied paired bulk-tissue and single-nucleus RNA-sequencing to anterior cingulate cortex samples derived from 28 individuals, including healthy controls, PD, PDD and DLB cases (n = 7 per group), to transcriptomically profile the LBDs. Using this approach, we (i) found transcriptional alterations in multiple cell types across the LBDs; (ii) discovered evidence for widespread dysregulation of RNA splicing, particularly in PDD and DLB; (iii) identified potential splicing factors, with links to other dementia-related neurodegenerative diseases, coordinating this dysregulation; and (iv) identified transcriptomic commonalities and distinctions between the LBDs that inform understanding of the relationships between these three clinical disorders. Together, these findings have important implications for the design of RNA-targeted therapies for these diseases and highlight a potential molecular “window” of therapeutic opportunity between the initial onset of PD and subsequent development of Lewy body dementia.

  Study EGAS00001005305

• Evolutionary predictability of genetic versus non genetic resistance to anticancer drugs in melanoma

  Therapy resistance arises from heterogeneous drug-tolerant persister cells or minimal residual disease (MRD) through genetic and nongenetic mechanisms. A key question is whether specific molecular features of the MRD ecosystem determine which of these two distinct trajectories will eventually prevail. We show that, in melanoma exposed to MAPK therapeutics, emergence of a transient neural crest stem cell (NCSC) population in MRD concurs with the development of nongenetic resistance. This increase relies on a GDNF-dependent signaling cascade, which activates the AKT survival pathway in a focal adhesion kinase (FAK)-dependent manner. Ablation of the NCSC population through FAK inhibition delays relapse in patient-derived tumor xenografts. Strikingly, all tumors that ultimately escape this treatment exhibit resistance-conferring genetic alterations and increased sensitivity to ERK inhibition. These findings identify an approach that abrogates the nongenetic resistance trajectory in melanoma and demonstrate that the cellular composition of MRD deterministically imposes distinct drug resistance evolutionary paths.

  Study EGAS00001005314

• Mutational bias in spermatogonia impacts the anatomy of regulatory sites in the human genome

  Mutation in the germline is the ultimate source of genetic variation, but little is known about the influence of germline chromatin structure on mutational processes. Using ATAC-seq, we profile the open chromatin landscape of human spermatogonia, the most proliferative cell-type of the germline, identifying transcription factor binding sites (TFBSs) and PRDM9-binding sites, a subset of which will initiate meiotic recombination. We observe an increase in rare structural variant (SV) breakpoints at PRDM9-bound sites, implicating meiotic recombination in the generation of structural variation. Many germline TFBSs, such as NRF, are also associated with increased rates of SV breakpoints, apparently independent of recombination. Singleton short insertions (>=5 bp) are highly enriched at TFBSs, particularly at sites bound by testis active TFs, and their rates correlate with those of structural variant breakpoints. Short insertions often duplicate the TFBS motif, leading to clustering of motif sites near regulatory regions in this male-driven evolutionary process. Increased mutation loads at germline TFBSs disproportionately affect neural enhancers with activity in spermatogonia, potentially altering neurodevelopmental regulatory architecture. Local chromatin structure in spermatogonia is thus pervasive in shaping both evolution and disease.

  Study EGAS00001005366

• SG Peranakan Project: Genetic admixture in the culturally unique Peranakan Chinese population in Southeast Asia

  The Peranakan Chinese are culturally unique descendants of immigrants from China who settled in the Malay Archipelago ~300-500 years ago. Today, among large communities in Southeast Asia, the Peranakans have preserved Chinese traditions with strong influence from the local indigenous Malays. Yet, whether or to what extent genetic admixture co-occurred with the cultural mixture has been a topic of ongoing debate. We performed whole-genome sequencing (WGS) on 177 Singapore (SG) Peranakans and analyzed the data jointly with WGS data of Asian and European populations. We estimated that Peranakan Chinese inherited ~5.62% (95% confidence interval [CI]: 4.75-6.46%) Malay ancestry, much higher than that in SG Chinese (1.08%, 0.69-1.53%), southern Chinese (0.86%, 0.57-1.31%), and northern Chinese (0.25%, 0.18-0.33%). A sex-biased admixture history, in which the Malay ancestry was contributed primarily by females, was supported by X chromosomal variants, and mitochondrial (MT) and Y haplogroups. Finally, we identified an ancient admixture event shared by Peranakan Chinese and SG Chinese ~1,612 (95% CI: 1,345-1,923) years ago, coinciding with the settlement history of Han Chinese in southern China, apart from the recent admixture event with Malays unique to Peranakan Chinese ~190 (159-213) years ago. These findings greatly advance our understanding of the dispersal history of Chinese and their interaction with indigenous populations in Southeast Asia.

  Study EGAS00001005379

• Analysis of RAD51C promoter methylation using targeted bisulfite sequencing (amplicon sequencing) in ovarian cancer pre-clinical models and patient samples.

  In high-grade serous ovarian carcinoma (HGSC), deleterious mutations in the DNA repair gene RAD51C are established drivers of defective homologous recombination and are emerging biomarkers of PARP inhibitor (PARPi) sensitivity. RAD51C promoter methylation (meRAD51C) is detected at similar frequencies to mutations, yet its effects on PARPi responses remain unresolved. In this study, three HGSC patient-derived xenograft (PDX) models with methylation at most or all examined CpG sites in the RAD51C promoter show responses to PARPi. Both complete and heterogeneous methylation patterns were associated with RAD51C gene silencing and homologous recombination deficiency (HRD). PDX models lost meRAD51C following treatment with PARPi rucaparib or niraparib, where a single unmethylated copy of RAD51C was sufficient to drive PARPi resistance. Genomic copy number profiling of one of the PDX models using SNP arrays revealed that this resistance was acquired independently in two genetically distinct lineages. In a cohort of 11 patients with RAD51C-methylated HGSC, various patterns of meRAD51C were associated with genomic 'scarring', indicative of HRD history, but exhibited no clear correlations with clinical outcome. Differences in methylation stability under treatment pressure were also observed between patients, where one HGSC was found to maintain meRAD51C after 6 lines of therapy (4 platinum-based), whilst another HGSC sample was found to have heterozygous meRAD51C and elevated RAD51C gene expression (relative to homozygous meRAD51C controls) after only neo-adjuvant chemotherapy. As meRAD51C loss in a single gene copy was sufficient to cause PARPi resistance in PDX, methylation zygosity should be carefully assessed in previously treated patients when considering PARPi therapy.

  Study EGAS00001005395

• Germline RUNX1 Variation and Predisposition to Childhood Acute Lymphoblastic

  Somatic mutations of RUNX1, which encodes the myeloid and lymphoid transcriptional factor RUNX1, are common in both B- and T- acute lymphoid leukemia (ALL) and are associated with poor prognosis of T-ALL. However, there has been no comprehensive investigation of the pattern or prevalence of RUNX1 germline mutation in both B- and T-ALL. Here we report germline RUNX1 variants in 1.23% of B-ALL and 2.11% of T-ALL, identifying 31 unique variants in 62 B-ALL and 18 unique variants in 26 T-ALL children. The majority of frameshift and nonsense variants affected RUNX1 function in transcriptional regulation, hematopoiesis, and cellular proliferation. We identified JAK3 as the most frequent somatic mutation in T-ALL with RUNX1 variants. These results not only identify RUNX1 as a leukemia predisposition gene but also further underline the importance of germline genetic variants to the development of ALL

  Study EGAS00001005403

• Interferon lambda 4 impairs viral antigen presentation and attenuates T cell responses

  Genetic variants of the interferon lambda (IFNL) gene locus are strongly associated with spontaneous and IFN treatment induced clearance of hepatitis C virus (HCV) infections. Individuals with the ancestral IFNL4-dG allele are not able to clear HCV in the acute phase and develop chronic hepatitis C (CHC) with more than 90% probability. Paradoxically, the IFNL4-dG allele encodes a fully functional IFNλ4 protein with antiviral activity against HCV and coronavirus. Here we describe a hitherto unknown effect of IFNλ4 on HCV antigen presentation. Only minor amounts of IFNλ4 are secreted, because the protein is largely retained in the endoplasmic reticulum (ER) where it induces ER stress. Stressed cells were significantly weaker activators of HCV specific CD8+ T cells. This was not due to reduced MHC I surface presentation or extracellular IFNλ4 effects, since T cell responses could be restored by exogenous loading of MHC with HCV antigens. Rather, IFNλ4 induced ER stress impaired HCV antigen processing and/or loading onto the MHC I complex. Our results provide a potential explanation for the IFNλ4 – HCV paradox.

  Study EGAS00001005396

• Genomic characterization of co-existing biliary tract intraepithelial neoplasia and carcinoma lesions reveals distinct evolutionary paths of gallbladder cancer

  Gallbladder carcinoma is the most common cancer of the biliary tract with dismal survival largely due to delayed diagnosis. Biliary tract intraepithelial neoplasia (BilIN) is the common benign tumor that is suspected to be precancerous lesions. However, the genetic and evolutionary relationships between BilIN and carcinoma remain unclear. Here we performed whole-exome sequencing of coexisting low-grade BilIN (adenoma), high-grade BilIN, and carcinoma lesions, and normal tissues from the same patients.

  Study EGAS00001005402

• Philippine Ayta possess the highest level of Denisovan ancestry in the world

  Multiple lines of evidence show that modern humans interbred with archaic Denisovans. Here, we report an account of shared demographic history between Australasians and Denisovans distinctively in Island Southeast Asia. Our analyses are based on ~2.3 million genotypes from 118 ethnic groups of the Philippines including 25 diverse self-identified Negrito populations, along with high-coverage genomes of AustraloPapuans and Ayta Magbukon Negritos. We show that Ayta Magbukon possess the highest level of Denisovan ancestry in the world – ~30-40% greater than that of Australians and Papuans – consistent with an independent admixture event into Negritos from Denisovans. Together with the recently described Homo luzonensis, we suggest that there were multiple archaic species that inhabited the Philippines prior to the arrival of modern humans, and that these archaic groups may have been genetically related. Altogether, our findings unveil a complex intertwined history of modern and archaic humans in the Asia-Pacific region, where distinct Islander Denisovan populations differentially admixed with incoming Australasians across multiple locations and at various points in time.

  Study EGAS00001005407

• Philippine Ayta possess the highest level of Denisovan ancestry in the world

  Multiple lines of evidence show that modern humans interbred with archaic Denisovans. Here, we report an account of shared demographic history between Australasians and Denisovans distinctively in Island Southeast Asia. Our analyses are based on ~2.3 million genotypes from 118 ethnic groups of the Philippines including 25 diverse self-identified Negrito populations, along with high-coverage genomes of AustraloPapuans and Ayta Magbukon Negritos. We show that Ayta Magbukon possess the highest level of Denisovan ancestry in the world – ~30-40% greater than that of Australians and Papuans – consistent with an independent admixture event into Negritos from Denisovans. Together with the recently described Homo luzonensis, we suggest that there were multiple archaic species that inhabited the Philippines prior to the arrival of modern humans, and that these archaic groups may have been genetically related. Altogether, our findings unveil a complex intertwined history of modern and archaic humans in the Asia-Pacific region, where distinct Islander Denisovan populations differentially admixed with incoming Australasians across multiple locations and at various points in time.

  Study EGAS00001005408

• Whole genome sequencing (bam files) of 5 samples of myxofibrosarcoma and 5 matched pairs

  Myxofibrosarcoma (MFS) is a rare subtype of sarcomas in the elderly, whose genetic basis is poorly understood. To elucidate it, the whole genome sequence was performed.

  Study EGAS00001005444

• Genetics causes of male infertility in 185 patient-parent trios from Netherlands and UK

  In this, whole exome sequencing was performed on 185 patients suffering from Azoospermia or severe oligozoospermia and on their respective fertile parents in order to investigate de novo and inherited genetic variants likely to be the cause of infertility in these men. After bioinformatic analysis and comprehensive variant interpretation several de novo point mutations and copy number variations were identified as well as multiple homozygous and compound heterozygous mutations in recessive genes that are likely influencing the infertility phenotypes observed in the patients.

  Study EGAS00001005417

• Feasibility of Ultra-Rapid Exome Sequencing in Critically Ill Infants and Children With Suspected Monogenic Conditions in the Australian Public Health Care System

  Critically ill infants and children with suspected monogenic conditions underwent ultra-rapid whole exome genetic testing. A molecular diagnosis was established in 51% of the patients. This study suggests feasibility of ultra-rapid genomic testing in critically ill pediatric patients.

  Study EGAS00001005430

• Spatial concordance of DNA methylation classification in diffuse glioma

  Background: Intratumoral heterogeneity is a hallmark of diffuse gliomas. DNA methylation profiling is an emerging approach in the clinical classification of brain tumors. The goal of this study is to investigate the effects of intratumoral heterogeneity on classification confidence. Methods: We used neuronavigation to acquire 133 image-guided and spatially-separated stereotactic biopsy samples from 16 adult patients with a diffuse glioma (7 IDH-wildtype and 2 IDH-mutant glioblastoma, 6 diffuse astrocytoma, IDH-mutant and 1 oligodendroglioma, IDH-mutant and 1p19q codeleted), which we characterized using DNA methylation arrays. Samples were obtained from regions with and without abnormalities on contrast enhanced T1 weighted and fluid-attenuated inversion recovery MRI. Methylation profiles were analyzed to devise a three-dimensional reconstruction of (epi)genetic heterogeneity. Tumor purity was assessed from clonal methylation sites. Results: Molecular aberrations indicated that tumor was found outside imaging abnormalities, underlining the infiltrative nature of this tumor and the limitations of current routine imaging modalities. We demonstrate that tumor purity is highly variable between samples and explains a substantial part of apparent epigenetic spatial heterogeneity. We observed that DNA methylation subtypes are often, but not always, conserved in space taking tumor purity and prediction accuracy into account. Conclusion: Our results underscore the infiltrative nature of diffuse gliomas and suggest that DNA methylation subtypes are relatively concordant in this tumor type, although some heterogeneity exists.

  Study EGAS00001005434

• Whole exome sequencing (bam files) of 55 samples of myxofibrosarcoma and 44 matched pairs

  Myxofibrosarcoma (MFS) is a rare subtype of sarcomas in the elderly, whose genetic basis is poorly understood. To elucidate it, the whole exome sequence was performed.

  Study EGAS00001005442

• Targeted-capture sequencing (bam files) of 81 samples of myxofibrosarcoma and 44 matched pairs

  Myxofibrosarcoma (MFS) is a rare subtype of sarcomas in the elderly, whose genetic basis is poorly understood. To elucidate it, the Targeted-capture sequencing was performed.

  Study EGAS00001005443

• Next Generation Sequencing platform for targeted Personalized Therapy of Leukemia - Acute myeloid leukemia

  The advent of high-throughput next generation sequencing (NGS) technologies that are revolutionizing genomics and transcriptomics by providing a single base resolution tool for a unified deep analysis of diseases complexity allows a fast and cost-efficient fine-scale assessment of the genetic variability hidden within cohorts of patients affected by the same leukemia. That being so, by potentially highlighting inter-individual differences that may play a role in the differential success of diverse therapeutic interventions, they promise to be crucial for selecting the most appropriate medical treatments. The project aims at the identification of novel prognostic biomarkers for acute myeloid leukemia (AML) and studied the molecular differences between aneuploid and euploid AML.

  Study EGAS00001005464

• The PEMDAC phase 2 study of pembrolizumab and entinostat in patients with metastatic uveal melanoma

  Patients with metastatic uveal melanoma (UM) have an abysmal prognosis. Preclinical studies have suggested that epigenetic therapy could enhance immunogenicity of cancer cells. Here we test if epigenetic therapy would enhance PD-1 immunotherapy in patients with metastatic UM. We report the results of the PEMDAC phase 2 clinical trial (n=29; NCT02697630) where the HDAC inhibitor entinostat was combined with the PD-1 inhibitor pembrolizumab in patients with metastatic UM. The primary endpoint was objective response rate (ORR), and was met with an ORR of 14%. The clinical benefit rate at 18 weeks was 28%, median progression free survival was 2.1 months and the median overall survival was 13.4 months. Toxicities were manageable, and there were no treatment-related deaths. Extensive genomics studies were performed using DNA/RNA and single cell sequencing and flow cytometry. Objective responses and prolonged survival were seen in three patients with BAP1 wildtype tumors, and in one patient with an iris melanoma that exhibited a UV signature. Longer survival also correlated with low baseline ctDNA levels or LDH. In conclusion, HDAC inhibition and anti-PD1 immunotherapy results in durable responses in a subset of patients with metastatic UM. Further exploration of combined immunotherapy and epigenetic therapy in metastatic UM is warranted.

  Study EGAS00001005478

• An epigenomics time course analysis of covid19 patients from Quebec, Canada

  Investigating the cellular and epigenetic changes at play in severely ill hospitalized COVID-19 patients at admission and how it evolves with prognosis. Single-cell RNA-seq analysis done on longitudinal samples highlighted changes in the monocytic compartment, which we explored using whole-genome epigenomics methods (ATACseq abd WGBS).

  Study EGAS00001005468

• Epigenetic encoding, heritability and plasticity of glioma transcriptional cell states

  Study EGAS00001005472

• Role of cohesin/CTCF in human monocyte differentiation

  Here we study genome architecture in a postmitotic differentiation setting, the differentiation of human blood monocytes (MO). We profiled and compared epigenetic, transcriptome and 3D conformation landscapes during MO differentiation (either into dendritic cells or macrophages) across the genome, with and without siRNA-mediated depletion of RD21 or CTCF, and detected numerous architectural changes, ranging from higher level compartments

  Study EGAS00001005508

• Methylation clocks - individual colon, small intestine and endometrial crypts

  Study EGAS00001005514

• Multi-omics integration reveals only minor long-term molecular and functional sequelae in immune cells of individuals recovered from COVID-19

  The majority of COVID-19 patients experience mild to moderate disease course and recover within a few weeks. An increasing number of studies characterized the long-term changes in the specific anti-SARS-CoV-2 immune responses, but how COVID-19 shapes the innate and heterologous adaptive immune system after recovery is less well known. To comprehensively investigate the post-SARS-CoV-2 infection sequelae on the immune system, we performed a multi-omics study by integrating single-cell RNA-sequencing, single-cell ATAC-sequencing, genome-wide DNA methylation profiling, and functional validation experiments in convalescent COVID-19 and healthy individuals. We showed that immune responses generally recover without major sequelae after COVID-19. However, subtle differences persist at the transcriptomic level in monocytes, with downregulation of the interferon pathway, while DNA methylation also displays moderate changes in convalescent COVID-19 individuals. However, these differences did not affect the cytokine production capacity of PBMCs upon different bacterial, viral, and fungal stimuli, although baseline release of IL-1RA and IFN?? was higher in convalescent individuals. In conclusion, despite minor differences in epigenetic and transcriptional programs, the immune system of convalescent COVID-19 patients largely recovers to the homeostatic level of healthy individuals.

  Study EGAS00001005529

• Pancreatic islets PISA RNA-seq samples

  GWAS have identified more than 700 genetic signals associated with type 2 diabetes (T2D). To gain insight into the underlying molecular mechanisms, we created the Translational human pancreatic Islet Genotype tissue-Expression Resource (TIGER), aggregating >500 human islet RNA-seq and genotyping datasets. We imputed genotypes using 4 reference panels and meta-analyzed cohorts to improve coverage of expression quantitative trait loci (eQTL) and developed a method to combine allele-specific expression across samples (cASE). We identified >1 million islet eQTLs (56% novel), of which 53 colocalize with T2D signals (60% novel). Among them, a low-frequency allele that reduces T2D risk by half increases CCND2 expression. We identified 8 novel cASE colocalizations, among which an SLC30A8 T2D associated variant. We make all the data available through the open-access TIGER portal (http://tiger.bsc.es), which represents a comprehensive human islet genomic data resource to elucidate how genetic variation affects islet function and translate this into therapeutic insight and precision medicine for T2D.

  Study EGAS00001005535

• Investigation of the keratinocytic gene expression pattern in Hidradenitis suppurativa

  The underlying pathogenetic factors generating the innate immune signal necessary for T cell activation, initiation and chronification of Hidradenitis suppurativa (HS, also known as Acne inversa) are still poorly understood. Emerging evidence suggests that defective keratinocyte function critically contributes to HS disease development and progression. To elucidate the role of keratinocytes in HS lesion formation, we compared the transcriptomes of isolated lesional and perilesional HS epidermis by RNA sequencing.

  Study EGAS00001005544

• Genomic characterization of retinoblastoma (targeted sequencing)

  Retinoblastoma is the most frequent intraocular malignancy in children, originating from a maturing cone precursor in the developing retina. Little is known on the molecular basis underlying the biological and clinical behavior of this cancer. Here, using multi-omics data, we demonstrate the existence of two retinoblastoma subtypes. Subtype 1, of earlier onset, includes most of the heritable forms. It harbors few genetic alterations other than the initiating RB1 inactivation and corresponds to differentiated tumors expressing mature cone markers. By contrast, subtype 2 tumors harbor frequent recurrent genetic alterations including MYCN-amplification. They express markers of less differentiated cone together with neuronal/ganglion cell markers with marked inter- and intra-tumor heterogeneity. The cone dedifferentiation in subtype 2 is associated with stemness features including low immune and interferon response, E2F and MYC/MYCN activation and a higher propensity for metastasis. The recognition of these two subtypes, one maintaining a cone-differentiated state, and the other, more aggressive, associated with cone dedifferentiation and expression of neuronal markers, opens up important biological and clinical perspectives for retinoblastomas. In this study, we performed targeted sequencing of the exonic regions of RB1, BCOR and ARID1A on 37 samples from the series of 102 retinoblastomas. Among them, 23 samples were not subjected to whole-exome sequencing, and 14 samples have the whole-exome sequencing data (deposited in EGAS00001005248).

  Study EGAS00001005550

• Pediatric study using genome sequencing

  This study used whole genome sequencing to assess genome variants from DNA from pediatric individuals seen in genetics clinic who were thought to have genetic conditions.

  Study EGAS00001005553

• DLBCL NGS Genomic Datasets of non-China cohort from Phoenix Clinical Trial

  Study biopsies (n=773; 189 cases as China cohort and 584 cases as non-China cohort) underwent whole exome sequencing or targeted gene panel DNA sequencing and RNA sequencing, allowing the identification of somatic mutations. Tumors were assigned to DLBCL genetic subtypes as the MCD (n=110), BN2 (n=47) or N1 (n=28) using the LymphGen algorithm based on their somatically mutated genes. In younger patients (age≤60), ibrutinib plus R-CHOP was associated with a 100% event free survival (EFS) at 3-years in the MCD and N1 groups, whereas R-CHOP alone was associated with a 3-year EFS of 42.9% (p=0.01) and 50% (p=0.016) in the MCD and N1 groups, respectively. This analysis showed a benefit of ibrutinib with R-CHOP chemotherapy in particular genetic subtypes, providing mechanistic support for its survival benefit in younger patients with non-GCB DLBCL. The EGA collects RNA-Seq and NuGene targeted sequencing datasets of 584 cases for non-China cohort.

  Study EGAS00001005554

• Mitochondrial DNA sequencing of iPSC, parental cells, and iPSC derived cardiomyocytes

  The aim of this work was to analyze to which extent mitochondrial DNA (mtDNA) variants in human induced pluripotent stem cells (iPSC) originate from their parental cells or from de novo mutagenesis, and whether dynamics in variant heteroplasmy levels during prolonged expansion culture of iPSC are caused by inter- or intra-cellular selection or genetic drift. Hence, mtDNA of 26 early passage iPSC clones from, in total, 10 donors and the corresponding parental cell populations (here endothelial cells) were sequenced. Furthermore, 7 iPSC clones were expanded up to passage 50 and the change in variant heteroplasmy levels was monitored. Lastly, the effect of differentiation on the variant heteroplasmy levels was investigated during directed cardiomyocyte differentiation of 4 iPSC clones.

  Study EGAS00001005560

• Whole genome sequencing of a breast cancer cohort with known functional homologous recombination status

  Optimal patient selection for PARP inhibitors and double strand break (DSB) inducing chemotherapy is of great clinical importance. In principle, the use of these therapies can be extended beyond germline BRCA1/2 mutated cancers. Although several methods for detection of BRCA-like tumors have been developed, a golden standard is lacking. We have developed the RECAP test, a functional HR assay exploiting the formation of RAD51 foci in proliferating cells after ex vivo irradiation of fresh primary breast cancer (BrC) tissue. The RECAP test seems robust and has an advantage over genetic tests, due to its functional character for exploring the HR phenotype. Another type of HRD test is the BRCAness classifier, which is based on specific genomic patterns derived from copy number data of BRCA1/2 mutated BrCs that also occur in sporadic cancers. The BRCAness classifier is a robust test that is easily applicable in the clinic and has also been validated to predict in vivo response to high dose chemotherapy in clinical trials. More recently, a new HRD test based on genomic signatures which can be extracted from whole genome sequencing has been published, HRDetect. Studies comparing these different HRD tests within the same cohort of patients do not yet exist. Here, we have the unique opportunity to compare three HRD tests within the same cohort. We aim to compare three different HRD tests (RECAP, BRCAness classifier and HRDetect) within the same cohort of primary BrCs (n=41). As a secondary aim, we will analyze additional genetic alterations in HRD tumors (both BRCA-deficient and BRCA-proficient HRD tumors).

  Study EGAS00001005572

• Genetics and therapeutic responses to TIL therapy of pancreatic cancer PDX models

  Pancreatic cancer is the 7th leading cause of cancer-related deaths worldwide. Checkpoint immunotherapy has not yet showed encouraging results in pancreatic cancer possibly owing to a poor immunogenicity and/or an immune suppressive microenvironment. The aim of this study was to develop patient-derived xenografts (PDX) models and to assess if autologous tumor-infiltrating lymphocytes (TILs) would have anti-tumoral activity in pancreatic cancer. Tumor biopsies from 29 patients were subcutaneously transplanted into NOG mice. Tumor growths were confirmed in 11 out of 29 transplantations. The PDX tumors histologically resembled their original biopsies, but since stromal cells in the PDX model tumors were from mouse, their gene expression differed from the original biopsies. Immune checkpoint ligands other than PD-L1 were expressed in pancreatic cancers but PD-L1 was rarely expressed. One of the three tumors that did express PD-L1 was an adenosquamous cancer and another had a mismatch repair deficiency. TILs were expanded from six tumors and were injected into NOG or human interleukin-2 transgenic-NOG (hIL2-NOG) mice carrying PDX tumors. Regression of tumors could be verified in hIL2-NOG mice in 3 of the 6 PDX models treated with autologous TILs, including the adenosquamous PDX model. In conclusion, PDX models of pancreatic cancer can be used to learn more about tumor characteristics and biomarkers, and to evaluate responses to ACT and combination therapies. The major benefit of the model is that modifications of T cells can be tested in an autologous humanized mouse model to gain preclinical data to support the initiation of a clinical trial.

  Study EGAS00001005596

• IVF Whole genome prediction

  Preimplantation Genetic Testing (PGT) of in vitro fertilized (IVF) embryos has been proposed as a method to reduce transmission of common disease, although comprehensive genetic assessment remains unavailable. Here we use a combination of molecular and statistical techniques to reliably and accurately infer inherited genome sequence in 110 embryos. We combine polygenic risk scores (PRS) and rare variation to model susceptibility across 12 common conditions.

  Study EGAS00001005619

• 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

• SudanMitoSeq: Sudanese mitochondrial sequencing

  In various contexts, mitochondrial function or dysfunction can be linked to mitochondrial genome variations. The use of mitochondrial genetics thus promises personalized diagnostics and treatments. In order to devise specific precision medicine approaches based on mitochondrial genetic variation and to test them within clinical studies, control data is indispensable. Such control data needs to comprehensively cover genetic variation commonly observed and thus expected. In the context of this study we assessed whether current, world-wide mitochondrial data sufficiently represents the region of North and East Africa, that is, whether current reference data is ready for precision medicine in this region. Towards this, we sequenced mitochondrial genomes of 159 Sudanese individuals provided as part of this EGA study and analyzed them together with various other, publicly available data concerning mitochondrial variation. For details, please refer to the publication.

  Study EGAS00001005669

• Germline biallelic mutation affecting the transcription factor Helios causes pleiotropic defects of immunity

  We studied a patient with recurrent respiratory infections and hypogammaglobulinemia and identified a germline homozygous missense mutation in IKZF2 encoding Helios (p.Ile325Val). We show that HeliosI325V retains DNA-binding and dimerization properties, but loses interaction with several partners, including epigenetic remodelers HDAC1, HDAC3 and the ATAC complex. Single-cell RNA-sequencing of peripheral blood mononuclear cells revealed gene expression signatures indicative of a shift towards pro-inflammatory, effector-like status in the patient’s T cells. We observed an upregulation of the anti-inflammatory gene, TSC22D3, encoding the glucocorticoid-induced leucine zipper (GILZ), whose expression was reported to increase upon IL-2 deprivation. We validated this finding in primary T cells, where we observed a pronounced reduction in IL-2 production upon their stimulation, together with a defect in their ability to proliferate, rendering them more anergic. Collectively, we identify a novel germline-encoded inborn error of immunity and define a role for Helios in conventional T cells, whereby interactions with specific binding partners is necessary to mediate the transcriptional programs that enable T-cell homeostasis in health and disease.

  Study EGAS00001005675

• 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

• 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

• Pediatric HGG WES and RNA-Seq

  Paediatric high-grade glioma (pHGG) is the leading cause of paediatric brain tumour death. Here we performed RNA and whole exome sequencing on pHGG samples to analyse transcriptomic and genetic alterations associated with this devastating disease, and to understand the relationship between them.

  Study EGAS00001005687

• Human genomic and phenotypic synthetic data for the study of rare diseases

  The purpose of this dataset is to facilitate development of technical implementations for rare disease data integration, analysis, discovery, and federated access. This synthetic dataset includes clinical and genomic data from 6 rare disease cases. It consists of 18 whole genomes (6 index cases with their parents) which have genetic background based on public human data sequenced in the context of the Illumina Platinum initiative (Eberle, MA et al. (2017)) and made available by the HapMap project (https://www.genome.gov/10001688/international-hapmap-project). In each of the cases, real causative variants correlating with the phenotypic data provided were spiked-in. The cases included in this synthetic dataset correspond to the following type of disorders: CASE 1- Congenital myasthenic syndrome (Autosomal Dominant -de novo variant) CASE 2- Macular dystrophy (Autosomal Dominant) CASE 3- Muscular dystrophy (Autosomal Recessive-compound heterozygous variants) CASE 4- Mitochondrial disorder (Autosomal Recessive-consanguineous case - homozygous variant) CASE 5- Breast cancer (Autosomal Dominant) CASE 6- Similar as case 1 for patient matchmaking tests: Congenital myasthenic syndrome (Autosomal Dominant-de novo variant) For each case you will be able to download the following data: clinical information (phenopackets per individual and pedigree per family), raw genomic data (FASTQ and BAMs) and processed genomic data (vcfs). When using the data, the following should be acknowledged: the RD-Connect GPAP (https://platform.rd-connect.eu/), EC H2020 project EJP-RD (grant # 825575), EC H2020 project B1MG (grant # 951724) and Generalitat de Catalunya VEIS project (grant # 001-P-001647).

  Study EGAS00001005702

• Genetic variability in exon 1 of the glucocorticoid receptor gene NR3C1 is associated with postoperative complications

  Patients having major surgery experience inflammation post-operatively, which may contribute to post-operative morbidity. Endogenous glucocorticoids (GCs) are an essential part of the stress response, but this response varies between individuals, which may in turn affect clinical outcome in general, and post-operative inflammation, specifically. Exon 1 of the NR3C1 gene, encoding the glucocorticoid receptor, contains an established region of differential regulation. DNA methylation patterns in this region have been found to differ between individuals. We investigated the methylation status and genotype in the Cytosine-phosphate-Guanine (CpG) island in exon 1 of NR3C1 in 24 patients who underwent major abdominal surgery and correlated it with post-operative complications. DNA was extracted from peripheral blood leukocytes and underwent targeted bisulfite sequencing of the CpG island. Complications were graded according to the Clavien-Dindo classification and 14 out of 24 had post-operative complications. Multifactorial and partial least square methods were used to analyse the data. A homogenous demethylated pattern was observed in all patients and no single CpG methylation was associated with post-operative complications. Four single nucleotide polymorphisms (SNPs) were significantly correlated with higher Clavien-Dindo scores. Genetic variability in the chr5:143,402,505-143,405,805 region of exon 1 of the glucocorticoid receptor gene NR3C1, but not DNA methylation, was associated with more severe post-operative complications in patients having major abdominal surgery. This indicates that the response to GCs may be of clinical importance for inflammatory conditions.

  Study EGAS00001005737

• Proteogenomics of chronic lymphocytic leukemia

  Cancer heterogeneity at the proteome level may explain differences in therapy response and prognosis beyond the currently established genomic and transcriptomic based diagnostics. The relevance of proteomics for disease classifications remains to be established in clinically heterogeneous cancer entities such as chronic lymphocytic leukemia (CLL). Here, we characterized the proteome and transcriptome in-depth alongside genetic and ex-vivo drug response profiling in a clinically well annotated CLL discovery cohort (n= 68). Unsupervised clustering of the proteome data revealed six subgroups. Five of these proteomic groups were associated with genetic features, while one group was only detectable at the proteome level. This new group was characterized by accelerated disease progression, high spliceosomal protein abundances associated with aberrant splicing, and low B cell receptor signaling protein abundances (ASB-CLL). We developed classifiers to identify ASB-CLL based on its characteristic proteome or splicing signature in two independent cohorts (n= 165, n= 169) and confirmed that ASB-CLL comprises about 20 % of CLL patients. The inferior overall survival observed in ASB-CLL was independent of both TP53- and IGHV mutation status. Our multi-omics analysis refines the classification of CLL and highlights the potential of proteomics to improve cancer patient stratification beyond genetic and transcriptomic profiling. This is the transcriptomics data of the discovery cohort.

  Study EGAS00001005746

• Genome-wide Rare Variant Score Associates With Morphological Subtypes of Autism Spectrum Disorder

  Defining different genetic subtypes of Autism Spectrum Disorder (ASD) can help predict developmental outcomes. Based on minor physical and major congenital anomalies, we categorized 325 Canadian children with ASD into dysmorphic and non-dysmorphic subgroupings. We developed a method to calculate a genome-wide rare variant score (GRVS) from whole-genome sequencing (WGS) data for each ASD proband. GRVS is a sum of the number of variants in morphology-associated coding and non-coding regions, weighted by their effect sizes. Probands with dysmorphic ASD had a significantly higher GRVS compared to those with non-dysmorphic ASD (P= 0.027). Using the polygenic transmission disequilibrium test, we observed an over-transmission of ASD-associated common variants in non-dysmorphic ASD probands (P= 2.9×10-3). We replicated our findings using 442 ASD probands that had accompanying morphology data in the Simons Simplex Collection. Our findings provide support for an alternative genomic classification of ASD subgroups using morphology data, which may inform intervention protocols.

  Study EGAS00001005753

• Searching for genetic modulators of the phenotypic heterogeneity in Brugada Syndrome

  Study EGAS00001005848

• Identification of Germline Monoallelic Mutations in IKZF2 in Patients with Immune Dysregulation.

  Helios, encoded by IKZF2, is a member of the Ikaros family of transcription factors with pivotal roles in T-follicular helper, NK- and T-regulatory cell physiology. Somatic IKZF2 mutations are frequently found in lymphoid malignancies. Although germline mutations in IKZF1 and IKZF3, encoding Ikaros and Aiolos, have recently been identified in patients with phenotypically similar immunodeficiency syndromes, the effect of germline mutations in IKZF2 on human hematopoiesis and immunity remains enigmatic. We identified germline IKZF2 mutations (one nonsense (p.R291X)- and 4 distinct missense variants) in six patients with systemic lupus erythematosus, immune thrombocytopenia or EBV-associated hemophagocytic lymphohistiocytosis. Patients exhibited hypogammaglobulinemia, decreased number of T-follicular helper and NKcells. Single-cell RNA sequencing of PBMCs from the patient carrying the R291X variant revealed upregulation of pro-inflammatory genes associated with T-cell receptor activation and T-cell exhaustion. Functional assays revealed the inability of HeliosR291X to homodimerize and bind target DNA as dimers. Moreover, proteomic analysis by proximity-dependent Biotin Identification revealed aberrant interaction of 3/5 Helios mutants with core components of the NuRD complex conveying HELIOS mediated epigenetic and transcriptional dysregulation.

  Study EGAS00001005874

• Genome-wide genotype data for 1,433 ni-Vanuatu

  Here, we generated genome-wide genotype data for 1,433 contemporary ni-Vanuatu and assessed their fine-scale genetic structure, in order to address central questions about the settlement of western Remote Oceania.

  Study EGAS00001005910

• Combination therapies to inhibit LCK tyrosine kinase and mTOR signaling in T-cell Acute Lymphoblastic Leukemia

  Relapse and refractory T cell acute lymphoblastic leukemia (T-ALL) has a poor prognosis and new combination therapies are sorely needed. This work uncovered potent drug synergies between AKT/mTORC1 inhibitors and the general tyrosine kinase-inhibitor, dasatinib. The combination effectively curbed T-ALL growth in human cell lines and primary human T-ALL samples, and was well tolerated and effective in suppressing leukemia growth in patient-derived xenografts grown in mice. Mechanistically, dasatinib inhibited phosphorylation and activation of the lymphocyte-specific protein tyrosine kinase (LCK) to blunt the T-cell receptor (TCR) signaling pathway and when complexed with mTORC1 inhibition, induced potent T-ALL cell killing through reducing MCL-1 protein expression. Analysis of responses across T-ALL subtypes and association with specific genetic mutations or expression patterns failed to identify any obvious biomarkers or correlations in predicting therapy responses, although there appeared to be a trend toward higher TCR/LCK pathway activity in samples that exhibited synergistic therapy responses.

  Study EGAS00001005945

• The Proteogenomic Subtypes of Acute Myeloid Leukemia

  Acute myeloid leukemia (AML) is an aggressive blood cancer with poor prognosis. We report a comprehensive proteogenomic analysis of bone-marrow biopsies from 252 uniformly treated AML patients to elucidate the molecular pathophysiology of AML in order to inform future diagnostic and therapeutic approaches. In addition to in-depth quantitative proteomics, our analysis includes cytogenetic profiling and DNA/RNA sequencing. We identify five proteomic AML subtypes, each reflecting specific biological features spanning genomic boundaries. Two of these proteomic subtypes correlate with patient outcome, but none are exclusively associated with specific genomic aberrations. Remarkably, one subtype (Mito-AML), which is only captured in the proteome, is characterized by high expression of mitochondrial proteins and confers poor outcome, with reduced remission rate and shorter overall survival upon treatment with intensive induction chemotherapy. Functional analyses reveal that Mito-AML is metabolically wired towards stronger complex I dependent respiration and is more responsive to treatment with the BCL2 inhibitor venetoclax.

  Study EGAS00001005950

• An Isogenic Cell Line Panel for Sequence-based Screening of Targeted Anti-cancer Drugs

  We describe the creation and characterization of an isogenic cell line panel representing common cancer pathways, with multiple features optimized for high-throughput screening. More than 1,800 cell lines from three normal human cells were generated using CRISPR-technologies. Surprisingly, we discovered most of these lines did not result in complete gene inactivation, despite integration of sgRNA at the desired genomic site. However, a subset of the lines harbored true, biallelic disruptions of the targeted tumor suppressor gene, yielding a final panel of 100 well-characterize lines covering 19 pathways frequently subject to loss of function in cancers. This panel included genetic markers optimized for sequence-based ratiometric assays for drug-based screening assays. To illustrate the potential utility of this panel, we developed a multiplexed high-throughput screen that identified Wee1 inhibitor MK-1775 as a selective growth inhibitor of cells with inactivation of TP53. These cell lines and screening approach should prove useful for researchers studying a variety of cellular and biochemical phenomena

  Study EGAS00001005974

• Germline TP53 mutations undergo copy number gain years prior to tumor diagnosis

  Li-Fraumeni syndrome (LFS) is a hereditary cancer predisposition syndrome associated with germline TP53 mutations. We performed whole-genome sequence (WGS) analysis of tumors from 22 TP53 mutation carriers. We observed somatic mutations affecting Wnt, PI3K/AKT signaling, epigenetic modifiers and homologous recombination genes as well as mutational signatures associated with prior chemotherapy. We identified near-ubiquitous early loss of heterozygosity of TP53, with gain of the mutant allele. This occurred earlier in these tumors compared to tumors with somatic TP53 mutations, suggesting the timing of this mark may distinguish germline from somatic TP53 mutations. Phylogenetic trees of tumor evolution, reconstructed from bulk and multi-region WGS, revealed that LFS tumors exhibit comparatively limited heterogeneity. Overall, our study delineates early copy number gains of mutant TP53 as a characteristic mutational process in LFS tumorigenesis, likely arising very early in life or in utero.years prior to tumor diagnosis.

  Study EGAS00001005982

• Host genotyping data from Dutch adult bacterial meningitis patients and linked bacterial genome sequences

  Streptococcus pneumoniae is a common nasopharyngeal colonizer, but can also cause life-threatening invasive diseases such as empyema, bacteremia and meningitis. Genetic variation of host and pathogen is known to play a role in invasive pneumococcal disease, though to what extent is unknown. This study includes 1146 samples of host genotyping data (genotyped and imputed) from Illumina Omni arrays. Samples were collected from adults (>16 yrs) patients with CSF confirmed bacterial meningitis in the Netherlands between 2006 and 2015. Metadata includes patient outcome, species of bacteria, and for 467 samples a link to an ENA run with the associated bacterial genome (S. pneumoniae only).

  Study EGAS00001005993

• Multimodal immunogenomic biomarker analysis of tumors from pediatric patients enrolled to a phase 1-2 study of single-agent atezolizumab

  We report herein an extensive exploratory biomarker analysis of refractory tumors taken from pediatric patients prior to receiving atezolizumab monotherapy in the phase 1-2 iMATRIX-atezolizumab trial (NCT02541604). A high percentage of CD8+ T cells and elevated protein levels of programmed cell death ligand 1 (PD-L1) were associated with progression-free survival (PFS). T-cell receptor (TCR) sequencing revealed that diverse infiltrating TCR repertoire at baseline was prognostic. We found no associations between panel-based tumor mutation burden (TMB) or specific genetic aberrations with PFS in this study. Through a pan-cancer gene co-expression network analysis, we developed a novel tumor-agnostic Pediatric Cytotoxicity and Antigen Presentation (PedCAP) signature that was associated with improved PFS in the iMATRIX-atezo study. Our study highlights features of immune response in pediatric cancers when treated with immune checkpoint inhibitors and provides a multi-biomarker pediatric immunogram framework to guide prospective clinical trials in pediatric cancers.

  Study EGAS00001006004

• Human genome-wide variations in the Massim region

  The Massim, a cultural region encompassing the southeastern tip of mainland Papua New Guinea (PNG) and nearby islands, is famous for the Kula trading network, in which different valuables circulate in different directions among the islands. To explore Massim genetic history, we generated genome-wide data from across the region and found variable levels of Papuan-related (indigenous) ancestry, including a novel ancestry associated with Rossel Island, and Austronesian-related ancestry that arrived later. We find genetic evidence for different patterns of contact across PNG, the Massim, and the Bismarck and Solomon Archipelagoes for Austronesian-related vs. Papuan-related ancestry, with more geographic restriction for the latter. Moreover, Kula-practicing groups share more genetic similarity than do other groups, and this likely predates the origin of Kula, suggesting that high between-group contact facilitated the formation of Kula. Our study provides the first comprehensive genome-wide assessment of Massim inhabitants and new insights into the fascinating Kula system.

  Study EGAS00001006010

• Drug-induced epigenomic plasticity reprograms circadian rhythm regulation to drive prostate cancer towards androgen-independence (RNA-seq)

  Understanding how prostate cancer cells adapt to AR-targeted interventions is critical for identifying novel drug targets to improve the clinical management of treatment-resistant disease. Our study revealed an enzalutamide-induced epigenetic plasticity towards pro-survival signaling, and uncovered circadian regulator ARNTL as an acquired vulnerability after AR inhibition, presenting a novel clinical lead for therapeutic development.

  Study EGAS00001006016

• Drug-induced epigenomic plasticity reprograms circadian rhythm regulation to drive prostate cancer towards androgen-independence (ChIP-seq)

  Understanding how prostate cancer cells adapt to AR-targeted interventions is critical for identifying novel drug targets to improve the clinical management of treatment-resistant disease. Our study revealed an enzalutamide-induced epigenetic plasticity towards pro-survival signaling, and uncovered circadian regulator ARNTL as an acquired vulnerability after AR inhibition, presenting a novel clinical lead for therapeutic development.

  Study EGAS00001006017

• Papuan Y chromosome Diversity Panel

  Island Southeast Asia and Oceania host one of the world’s richest assemblages of human phenotypic, linguistic and cultural diversity. Despite this, the region’s male genetic lineages are globally among the last to remain unresolved. We compiled ~9.7 Mb of Y chromosome sequence from a diverse sample of over 380 men from this region, including 152 first reported here. The granularity of this dataset allows us to fully resolve and date the regional Y chromosome phylogeny. This new high-resolution tree confirms two main population bursts: multiple rapid diversifications following the region’s initial settlement ~50 kya, and extensive expansions <6 kya. Notably, ~40-25 kya the deep rooting local lineages of C-M130, M-P256, and S-B254 show almost no further branching events in Island Southeast Asia, New Guinea and Australia, matching a similar pause in diversification seen in maternal mitochondrial DNA lineages. The main local lineages start diversifying ~25 kya, at the time of the Last Glacial Maximum. This improved Y chromosome topology highlights localized events with important historical implications, including pre-Holocene contact between Mainland and Island Southeast Asia, potential interactions between Australia and the Papuan world, and a sustained period of diversification following the flooding of the ancient Sunda and Sahul continents as the insular landscape observed today formed. The high-resolution phylogeny of the Y chromosome presented here thus enables a detailed exploration of

  Study EGAS00001006025

• CIRdb: Array genotype data

  The current inhabitants of the Canary Islands have a unique genetic makeup in the European diversity landscape due to the existence of African footprints from recent admixture events, especially of North African components (>20%). The underrepresentation of non-Europeans in genetic studies and the sizable North African ancestry, which is nearly absent from all existing catalogs of worldwide genetic diversity, justify the need to develop CIRdb, a population-specific reference catalog of natural genetic variation in the Canary Islanders. Based on array genotyping of the selected unrelated donors and comparisons against available datasets from European, sub-Saharan, and North African populations, we illustrate the intermediate genetic differentiation of Canary Islanders between Europeans and North Africans and the existence of within-population differences that are likely driven by genetic isolation. In addition, here we describe the overall design and the methods that are being implemented to develop CIRdb.

  Study EGAS00001006050

• Insights from genome-wide data from Thailand and Laos

  We provide the first detailed insights into the genetic profiles of Thai/Lao ethnolinguistic groups, which should be helpful for reconstructing human genetic history in MSEA and selecting populations for participation in ongoing whole genome sequence and biomedical studies

  Study EGAS00001006053

• Analysis of recent shared ancestry in a familial cohort identifies coding and noncoding autism spectrum disorder variants

  Autism spectrum disorder (ASD) is a collection of neurodevelopmental disorders characterized by deficits in social communication and restricted, repetitive patterns of behavior or interests. ASD is highly heritable, but genetically and phenotypically heterogeneous, reducing the power to identify causative genes. We performed whole genome sequencing (WGS) in an ASD cohort of 68 individuals from 22 families enriched for recent shared ancestry. We identified an average of 3.07 million variants per genome, of which an average of 112,512 were rare. We mapped runs of homozygosity (ROHs) in affected individuals and found an average genomic homozygosity of 9.65%, consistent with expectations for multiple generations of consanguineous unions. We identified potentially pathogenic rare exonic or splice site variants in 12 known (including KMT2C, SCN1A, SPTBN1, SYNE1, ZNF292) and 12 candidate (including CHD5, GRB10, PPP1R13B) ASD genes. Furthermore, we annotated noncoding variants in ROHs with brain-specific regulatory elements and identified putative disease-causing variants within brain-specific promoters and enhancers for 5 known ASD and neurodevelopmental disease genes (ACTG1, AUTS2, CTNND2, CNTNAP4, SPTBN4). We also identified copy number variants in two known ASD and neurodevelopmental disease loci in two affected individuals. In total we identified potentially etiological variants in known ASD or neurodevelopmental disease genes for ~61% (14/23) of affected individuals. We combined WGS with homozygosity mapping and regulatory element annotations to identify candidate ASD variants. Our analyses add to the growing number of ASD genes and variants and emphasize the importance of leveraging recent shared ancestry to map disease variants in complex neurodevelopmental disorders.

  Study EGAS00001006058

• Chromatin profiles classify castration-resistant prostate cancers suggesting therapeutic targets

  Untreated prostate cancers rely on androgen receptor (AR) signaling for growth and survival, forming the basis for the initial efficacy of androgen deprivation therapy (ADT). Yet, the disease can relapse and progress to a lethal stage termed castration-resistant prostate cancer (CRPC). Reactivation of AR signaling represents the most common driver of CRPC growth and next generation AR signaling inhibitors (ARSIs) are now used in combination with ADT as a first line therapy. However, ARSIs can result in selective pressure generating AR-independent tumors. The transition from AR-dependence frequently accompanies a change in phenotype resembling developmental trans-differentiation or ‘lineage plasticity’. Neuroendocrine prostate cancer, which lacks a defined pathologic classification, is the most studied type of lineage plasticity. However, most AR-null tumors do not exhibit neuroendocrine features and are classified as ‘double-negative prostate cancer’, the drivers of which are poorly defined. Lineage plasticity studies in CRPC are limited by the lack of genetically defined patient-derived models that recapitulate the disease spectrum. To address this, we developed a biobank of organoids generated from patient biopsies to study the landscape of metastatic CRPC and allow for functional validation assays. Proteins called transcription factors (TFs) are drivers of tumor lineage plasticity. To identify the key TFs that drive the growth of AR-independent tumors, we integrated epigenetic and transcriptomic data generated from CRPC models. By presenting a map of the chromatin accessibility and transcriptomic landscape of CRPC using a diverse biobank of organoids, PDXs and cell lines that recapitulate the heterogeneity of metastatic prostate cancer, we validate CRPC-AR and CRPC-NE subtypes and reveal two subtypes of AR-negative/low samples as well as their respective masterTFs. Additional analysis revealed a model in which YAP, TAZ, TEAD and AP-1 function together and drive oncogenic growth in CRPC-SCL samples. Overall, we show here how an approach to stratify CRPC patients into four subtypes using their transcriptomic signatures can potentially inform upon appropriate clinical decisions.

  Study EGAS00001006059

• Submission 6 - study_title 1

  Submission 6 - description 1

  Study EGAS50000000420

• High-throughput telomere length measurement at nucleotide resolution using the PacBio high fidelity sequencing platform

  Telomeres are specialized nucleoprotein structures at the ends of linear chromosomes that are essential for the stability as well as the complete replication of the human genome. The progressive shortening of steady-state telomere length in normal human somatic cells is a promising biomarker for age-associated diseases. However, there remain substantial challenges in quantifying telomere length in clinics and research laboratories due to the lack of accurate and high-throughput telomere length measurement method. Here, we describe a novel workflow to capture and analyze telomeres accurately at the single-nucleotide resolution. By using newly designed telobaits that are complementary to the single-stranded G-rich 3’ telomeric overhangs, we can specifically and efficiently capture full-length telomeres and subject them to single-molecule real-time (SMRT) sequencing. We show that the telomere length can be measured at nucleotide resolution using the PacBio high fidelity (HiFi) sequencing platform, which provides the necessary accuracy and throughput to uncover the steady-state telomere length distribution in human culture cell lines as well as clinical patient samples. This new method is a valuable tool not only in population-wide epidemiology studies but also for the longitudinal assessment of telomere attrition in individuals. Our results also reveal the extreme heterogeneity of telomeric variant sequences (TVSs) that is dispersed throughout the telomere repeat region. Importantly, the unique distribution pattern and sequence of TVSs in individuals may be used as a fingerprint for sample identification. In addition, the presence of TVSs disrupts the continuity of the canonical (5’-TTAGGG-3’)n telomere repeats, which could affect the binding of shelterin complexes at the chromosomal ends. Therefore, TVSs may act as a form of genetic polymorphisms that impairs telomere protection. Together with the absolute telomere length in each individual, TVSs may have profound implications in human aging and diseases.

  Study EGAS00001006103

• Differentially methylated CpGs in response to growth hormone administration in children with idiopathic short stature

  Background: Recombinant human growth hormone (rhGH) has shown a great growth-promoting potential in children with idiopathic short stature (ISS). However, the response to rhGH differs across individuals, largely due to genetic and epigenetic heterogeneity. Since epigenetic marks on the methylome can be dynamically influenced by GH, we performed a comprehensive pharmaco-epigenomics analysis of DNA methylation changes associated with long-term rhGH administration in children with ISS. Results: We measured DNA methylation profiles before and after GH treatment (with a duration of ~18 months in average) on 47 healthy children using customized methylC-seq capture sequencing. Their changes were compared and associated with changes in plasma IGF1 by adjusting sex, age, treatment duration, and estimated blood proportions. We observed a considerable inter-individual heterogeneity of DNA methylation changes responding to GH treatment. We identified 267 response-associated differentially methylated cytosines (DMC) that were enriched in promoter regions, CpG islands and blood cell-type specific regulatory elements. Furthermore, the genes associated with these DMCs were enriched in the biology process of "cell development", "neuron differentiation" and "developmental growth", and in the TGF-beta signaling pathway, PPAR Alpha pathway, endoderm differentiation pathway, adipocytokine signaling pathway as well as PI3K-Akt signaling pathway, and cAMP signaling pathway. Conclusion: Our study provides a first insight in DNA methylation changes associated with rhGH administration, which may help understand mechanisms of epigenetic regulation on GH-responsive genes.

  Study EGAS00001006115

• Genomic History of the Solomon Islands

  The aim of this study was to provide for a better understanding of the genetic relationships among and between populations of Near and Remote Oceania, with a focus on the initial settlement of the Solomon Islands and their role in the first dispersal of people from Near into Remote Oceania. Irina Pugach, Ana T Duggan, D Andrew Merriwether, Françoise R Friedlaender, Jonathan S Friedlaender, Mark Stoneking, "The Gateway from Near into Remote Oceania: New Insights from Genome-Wide Data", Molecular Biology and Evolution, Volume 35, Issue 4, April 2018, Pages 871–886

  Study EGAS00001006116

• Conserved features of TERT promoter duplications reveal an activation mechanism that mimics hotspot mutations in cancer

  Mutations in the TERT promoter are the single most common non-coding mutation in cancer and represent the genetic underpinnings of tumor cell immortality. Beyond the two most common point mutations, G228A and G250A, which selectively recruit the ETS factor GABP to activate TERT, the significance of other variants in the TERT promoter are unknown. We identified duplications of wildtype sequence within the core promoter region of TERT in 7 different cancer types that have strikingly similar features including size, insertion position, and inclusion of one of the native ETS motifs. Each duplication activates the TERT promoter to a similar level as G228A and G250A and is critically dependent on the insertion site. The GABP tetramer binds to the TERT duplicated promoter sequence by virtue of the native ETS motif and its duplicated version with precise spacing, and it is necessary for the transcriptional activation by all duplications tested. Spatiotemporal analysis in a multifocal glioblastoma shows the duplication is clonal and its activation of TERT is readily detectable at the single cell level and in bulk tumor tissue. We conclude that recurrent TERT promoter duplications of the native ETS sequence are functionally and mechanistically equivalent to the hotspot mutations that confer tumor cell immortality. The shared mechanism of these divergent somatic genetic alterations suggests a strong selective pressure for recruitment of the GABP tetramer to activate TERT.

  Study EGAS00001006118

• Loss of Epigenetic Barrier is Required for Enhancer Hijacking-Mediated Oncogenic Transcription

  Structural variants (SVs) involving enhancer hijacking can disrupt chromatin topologies to cause oncogene activation in cancer genomes, yet the molecular determinants for the transcriptional output of enhancer hijacking remain largely unknown. We developed a multimodal approach to integrate genome sequencing, chromosome conformation, and sequence-based deep learning for quantitative analysis of transcriptional effects and structural reorganization imposed by SVs in leukemic genomes. We identified candidate pathogenic SVs including recurrent t(5;14) translocations that cause the hijacking of BCL11B enhancers for oncogenic activation of TLX3-dependent transcriptional programs. By engineering patient-associated t(5;14) in isogenic leukemia cells, we uncovered an uncharacterized mechanism whereby DNA methylation serves as an epigenetic barrier to enhancer hijacking and loss of epigenetic barrier is a molecular determinant for the transcriptional output of pathogenic SVs. Hence, leveraging the epigenetic barriers of SV-mediated oncogenic programs may provide new opportunities to reprogram gene regulation as epigenetic therapies in human disease.

  Study EGAS00001006140

• Integrated single-cell profiling dissects cell-state-specific enhancer landscapes of human tumor infiltrating T cells.

  Study EGAS00001006141

• Mapping genetic effects on cell type-specific chromatin accessibility using single nucleus ATAC-seq

  Gene regulation is highly cell type-specific and understanding the function of non-coding genetic variants associated with complex traits requires molecular phenotyping at cell type resolution. In this study we performed single nucleus ATAC-seq (snATAC-seq) and genotyping in peripheral blood mononuclear cells from 13 individuals. We mapped chromatin accessibility QTLs (caQTLs) in each immune cell type and sub-type in a subset of 10 samples of European genetic ancestry.

  Study EGAS00001006184

• Enhancer profiling to identify identifies epigenetic markers of endocrine resistance in metastatic castration resistant prostate cancer patients

  Study EGAS00001006161

• Molecular landscape of blastic plasmacytoid dendritic cell neoplasm

  Blastic plasmacytoid dendritic cell neoplasm (BPDCN) is an aggressive malignancy assumed to originate from plasmacytoid dendritic cells (pDCs), which mostly affects the skin, bone marrow and lymph nodes and sequentially other organ systems. We applied paired WES/RNA-seq combined with genome-wide copy-number analysis to characterize 47 BPDCN patients regarding mutational drivers, cytogenetic aberrations and gene-expression profiles.

  Study EGAS00001006166

• Heat selection enables highly scalable methylome profiling in cell-free DNA for noninvasive monitoring of cancer patients

  Genome-wide analysis of cell-free DNA (cfDNA) methylation profile has been recognized as a promising approach for sensitive and specific detection of many cancers. However, scaling such genome-wide assays for clinical translation is impractical due to the high cost of whole genome bisulfite sequencing. We have shown that the small fraction of GC-rich genome is highly enriched in CpG sites and disproportionately harbors the majority of cancer-specific methylation signature. Here, we report on the simple but effective Heat enrichment of CpG-rich regions for Bisulfite Sequencing (Heatrich-BS) platform that allows for focused methylation profiling in these highly informative regions. Our novel method and bioinformatics algorithm enable accurate tumor burden estimation with high sensitivity and quantitative tracking of colorectal cancer patient’s response to treatment, at much reduced sequencing cost suitable for frequent monitoring. We also show, for the first time, tumor epigenetic subtyping from cfDNA using Heatrich-BS, which could enable patient stratification from non-invasive liquid biopsy. As such, Heatrich-BS holds great potential for highly scalable screening and regular monitoring of cancer using liquid biopsy.

  Study EGAS00001006198

• Somatic L1 retrotranspositions in normal human cells

  Genomic mutations accumulate in normal cells during lifetime; however, somatic mosaicism generated by mobilization of transposable elements is not well understood. Here, we explored whole-genome sequences of single-cell expanded clones to illustrate the mutational landscape produced by L1. In addition, we analyzed promoter DNA methylation status of retrotransposition-competent L1 and read-through RNA expression profiles to investigate the epigenetic regulation of L1 activity.

  Study EGAS00001006213

• Submission 58 - study_title 1

  Submission 58 - description 1

  Study EGAS50000000476

• Extensive patient-to-patient single nuclei transcriptome heterogeneity in pheochromocytomas and paragangliomas

  Pheochromocytomas (PC) and paragangliomas (PG) are rare neuroendocrine tumors of varied genetic makeup, associated with high cardiovascular morbidity and a variable risk of malignancy. We focused on PCPG tumors with germline SDHB and RET mutations, representing distinct prognostic groups with worse or better prognoses, respectively. We applied single-nuclei RNA sequencing (snRNA-seq) on tissue samples from 11 patients.

  Study EGAS00001006230

• Salivary Gland Cancer TSO500 dataset

  Salivary gland cancer (SGC) is a rare cancer for which systemic treatment options are limited. Therefore, it is important to characterize its genetic landscape in search for actionable aberrations. This research aimed to identify these actionable aberrations by NGS-based analysis of DNA (single and multiple nucleotide variants, copy number variants, microsatellite instability and tumor mutational burden) in a large cohort of SGC patients. DNA was extracted from archival tissue of 119 patients with various SGC subtypes and sequenced using a targeted NGS panel encompassing 523 cancer related genes (TruSight Oncology 500, TSO500). The 119 patients included 46 AdCC patients, 43 SDC patients, 15 MEC patients, 9 AciCC patients and 6 patients with other subtypes (1 secretory carcinoma, 1 polymorphic adenocarcinoma (PAC), 1 adenocarcinoma NOS, 1 myoepithelial carcinoma, 1 epithelial/myoepithelial carcinoma and 1 mixed PAC/myoepithelial carcinoma). Median exon coverage ranged from 80-904 (median 308) unique reads and the median percentage of exon coverage with at least 100 unique reads was median 96.6% in these 119 samples.

  Study EGAS00001006232

• Submission 56 - study_title 1

  Submission 56 - description 1

  Study EGAS50000000472

• Genetic Determinants of Mannose-binding Lectin Activity Predispose to Thromboembolic Complications in Critical COVID-19

  Study EGAS00001006266

• Submission 25 - study_title 1

  Submission 25 - description 1

  Study EGAS50000000486

• Genetic and epigenetic variation at regulatory regions contribute to cancer evolution under endocrine treatment

  Comprehensive profiling of hormone-dependent breast cancer (HDBC) has identified hundreds of protein-coding alterations contributing to cancer initiation, but only a handful have been linked to endocrine therapy resistance, potentially contributing to 40% of relapses. If other mechanisms underlie the evolution of HDBC under adjuvant therapy is currently unknown. In this work, we employ integrative functional genomics to dissect the contribution of cis-regulatory elements (CREs) to cancer evolution by focusing on 12 megabases of non-coding DNA, including clonal enhancers, gene promoters, and boundaries of topologically associating domains. Massive parallel perturbation in vitro reveals context-dependent roles for many of these CREs, with a specific impact on dormancy entrance and endocrine therapy resistance9. Profiling of CRE somatic alterations in a unique, longitudinal cohort of patients treated with endocrine therapies identifies non-coding changes involved in therapy resistance. Overall, our data uncover actionable transient transcriptional programs critical for dormant persister cells and unveil new regulatory nodes driving evolutionary trajectories towards disease progression.

  Study EGAS00001006340

• Presence of rare potential pathogenic variants in subjects under 65 years old with very severe or fatal COVID‑19

  Rare variants affecting host defense against pathogens could be involved in COVID‑19 severity and may help explain fatal outcomes in young and middle‑aged patients. Our aim was to report the presence of rare genetic variants in certain genes, by using whole exome sequencing, in a selected group of COVID‑19 patients under 65 years who required intubation or resulting in death (n = 44). To this end, different etiopathogenic mechanisms were explored using gene prioritization‑based analysis in which genes involved in immune response, immunodeficiencies or blood coagulation were studied. We detected 44 different variants of interest, in 29 different patients (66%). Some of these variants were previously described as pathogenic and were located in genes mainly involved in immune response. A network analysis, including the 42 genes with candidate variants, showed three main components, consisting of 25 highly interconnected genes related to immune response and two additional networks composed by genes enriched in carbohydrate metabolism and in DNA metabolism and repair processes. In conclusion, we have detected candidate variants that may potentially influence COVID‑19 outcome in our cohort of patients. Further studies are needed to confirm the ultimate role of the genetic variants described in the present study on COVID‑19 severity.

  Study EGAS00001006372