Search Results - EGA European Genome-Phenome Archive

Decoding Human Endogenous Retrovirus Expression in Liver Metastatic Colorectal Cancers: Implications for Diagnosis and Prognosis

The association between Human Endogenous RetroVirus (HERV) expression and cancer can be considered to both reflect the epigenetic deregulation of cancer and facilitate oncogenesis through its promoter/enhancer action on genes. More than 50% of colorectal cancers develop liver metastases and HERV has never been studied in this setting. We collected 400 RNAseq samples from more than 200 patients with primary and liver metastases from published and unpublished databases. We corroborated our findings with tandem mass spectrometry from primary colorectal cancer and single-cell RNAseq from liver metastatic colorectal cancer (lmCRC). We identified global stability in HERV expression between liver metastases and primary colorectal cancers, suggesting an early oncogenic footprint. Therefore, we proposed a list of seventeen HERV loci for lmCRC characterization with ERV-derived peptides but no expression in normal tissues. Among them, we discovered six prognostic markers after surgical resection of lmCRC. Four of them, HERVH_Xp22.32a, HERVH_20p11.23b, HERVH_13q33.3, HERVH_13q31.3, showed poor prognostic value when expressed (log-rank p-value 0.028, 0.0083, 9e-4, 0.05, respectively). HERVH_Xp22.2c (log-rank p-value 0.032) and HERVH_8q21.3b (only in multivariate models) were associated with better prognosis when expressed. These markers showed a cumulative effect on survival when expressed. Some HERVs were associated with decreased cytotoxic immune cells and most of them correlated with cell cycle pathways. Our findings provide new insights into lmCRC classification and suggest novel prognostic biomarkers and potential therapeutic targets.

Study EGAS50000000307

The Dynamic Immune Behavior of Primary and Metastatic Ovarian Carcinoma

Patients with high-grade serous ovarian carcinoma are usually diagnosed at an advanced stage and often develop resistance to conventional therapy. Epithelial ovarian tumors have been shown to be rich in CD8+ lymphocytes, which infiltrate both primary and omental sites. In recent years, several prominent studies have highlighted a positive association between the presence of these immune cells and a good prognosis; however, clinical trials of immune checkpoint inhibitors have so far shown unsatisfactory results in ovarian cancer. In this study, we used an innovative combination of single-cell sequencing and spatial transcriptomics to identify the molecular mechanisms that lead to immunosuppression in primary and metastatic HGSC. Primary tumors consistently showed a more active immune microenvironment than did omental tumors. In addition, while several immune cell subtypes were present in all samples, their relative abundance differed between primary and omental tumors. Finally, we found that primary tumors were mostly populated by dysfunctional CD4 and CD8 T cells in later stages of differentiation, while omental tumors were mostly populated by naïve immune cells; this, in turn, was correlated with expression changes in the interferon α and γ pathways in epithelial cells, showing that cross-communication between the epithelial and immune compartments is important for immune suppression in ovarian cancer. These findings can impact the timing and choice of immune-modulating drugs and should be considered when planning treatment for patients with high-grade serous ovarian carcinoma.

Study EGAS50000000038

Cancer-associated fibroblasts promote drug resistance in ALK-driven lung adenocarcinoma cells by upregulating lipid biosynthesis

Targeted therapy interventions using tyrosine kinase inhibitors (TKIs) provide encouraging treatment responses in ALK-rearranged lung adenocarcinomas, yet resistances occur almost inevitably. Apart from tumor cell-intrinsic resistance mechanisms, accumulating evidence supports a role of cancer-associated fibroblasts (CAFs) in affecting the therapeutic vulnerability of lung cancer cells. Here, we aimed to investigate underlying molecular networks shaping the therapeutic susceptibility of ALK-driven lung adenocarcinoma cells via tumor microenvironmental cues using three-dimensional (3D) spheroid co-culture settings. We show that CAFs promote therapy resistance of lung tumor cells against ALK inhibition by reducing apoptotic cell death and increasing cell proliferation. Using single-cell RNA-sequencing analysis, we show that genes involved in lipogenesis constitute the major transcriptional difference between TKI-treated homo- and heterotypic lung tumor spheroids. CAF-conditioned medium and CAF-secreted factors HGF and NRG1 were both able to promote resistance of 3D-cultured ALK-rearranged lung tumor cells via AKT signaling, which was accompanied by enhanced de novo lipogenesis and supression of lipid peroxidation. Notably, simultaneous targeting of ALK and SREBP-1 was able to overcome the established CAF-driven lipid metabolic-supportive niche of TKI-resistant lung tumor spheroids. Our findings highlight a crucial role of CAFs in mediating ALK-TKI resistance via lipid metabolic reprogramming and suggest new ways to overcome resistance towards molecular directed drugs by targeting vulnerabilities downstream of oncogenic signaling.

Study EGAS50000000135

Diagnostic Genomic Analysis is Prognostic in AYA ALL Patients Treated on a MRD-Stratified Paediatric Protocol

Background: Over the last decade, next generation sequencing has enabled classification of multiple new recurrent genomic drivers of acute lymphoblastic leukemia (ALL). We provide a comprehensive genomic survey of adolescent and young adult (AYA) patients treated on the Australasian Leukaemia and Lymphoma Group (ALLG) ALL06 study. ALL06 clinical trial described the efficacy and tolerability of AYA ALL patients aged 15-39 years treated on a single minimal residual disease (MRD)-stratified platform. In this study we have taken the ALL06 cohort and retrospectively performed comprehensive genomic analyses assessing response of AYA patients with various genomic subtypes treated with a paediatric approach. Methods: Genomic risk classification of B- and T-ALL patients (n=57) enrolled to the study was based on mRNA-Sequencing, MLPA, immunophenotyping and cytogenetics. Results: Using this approach, 36/40 (90%) B- and 13/17 (76.5%) T-ALL patients were classified according to genomic risk. We have linked outcome data with each patient’s genomic risk and found a strong correlation between adverse genomic risk and MRD at the end of consolidation, translating to inferior overall and relapse free survival. Patients with adverse risk genomics who achieved negative MRD status had improved responses compared to those with persistent MRD. Patients with standard risk genomics had excellent responses regardless of their MRD status. Conclusion: Our data argue strongly for incorporation of a genomic risk classification into future ALL treatment paradigms at the time of diagnosis.

Study EGAS50000000752

Ultra-sensitive tumor-informed ctDNA monitoring of treatment response in advanced esophagogastric cancer patients

To explore whether ultra-sensitive circulating tumor DNA (ctDNA) profiling enables earlier prediction of treatment response and detection of disease progression, we applied NeXT Personal, an ultra-sensitive bespoke tumor-informed liquid biopsy platform, to profile tumor samples from the KeyLargo study, a phase II trial in which metastatic esophagogastric cancer (mEGC) patients received capecitabine, oxaliplatin, and pembrolizumab. A total of 24 patients were evaluated, and all were ctDNA-positive at baseline. ctDNA levels varied from 406,067 down to 1.5 parts per million (PPM) with a median limit of detection of 2.03 PPM. ctDNA dynamics were highly correlated with changes in tumor size (ρ = 0.59, p = 7.3×10-9). Lack of early molecular response (50% or greater decrease in ctDNA levels at first available time point after 30 days, C2D1 or C4D1) was associated with worse overall survival (OS) (HR 6.6, 95% CI 1.8-24.1, p = 0.005) and progression-free survival (PFS) (HR 15.4, 95% CI 2.7-87.0, p = 0.002). Lack of molecular clearance of ctDNA was associated with worse OS (HR 6.9, 95% CI 1.5-30.8, p = 0.012) and PFS (HR 19.2, 95% CI 2.4-152.8, p = 0.005). Molecular progression (ctDNA increase) preceded imaging-derived progression by a median lead time of 65 days. These results suggest that ultra-sensitive liquid biopsy approaches could improve treatment decision-making for mEGC patients receiving chemotherapy and immunotherapy.

Study EGAS50000001224

Delineating pediatric brain tumor progression using single-nuclei sequencing

Pediatric brain cancers are the most common solid tumors in children with varying survival rates depending on molecular features of cancer cells and associated microenvironment. Several studies have used single-cell RNA sequencing to define transcriptional programs active in cancer cells and highlighted supportive roles of non-cancerous cells embedded in these tumors. Although application of single-cell RNA sequencing technology has expanded our understanding of intratumoral heterogeneity and underlying transcriptional mechanisms in brain tumors at diagnosis, mechanisms of disease remains elusive. Therefore, the purpose of this proposal is to characterize cellular and transcriptional states of malignant and non-malignant compartments at primary diagnosis and at relapse. In this study, we performed single-nuclei sequencing of 42 patient samples taken at various disease stages (initial diagnosis, progressive disease, and relapsed disease) from 20 patients with low-grade glioma or Atypical Teratoid Rhabdoid Tumor (ATRT). Bulk whole genome and RNA-seq data for these samples are available from the Children’s Brain Tumor Network and the Gabriella Miller Kids First Data Resource Center. Findings from our study will help define the mechanistic shifts of tumor and non-cancerous cells over course of disease progression. Collectively with the single-cell Pediatric Cancer Atlas funded, this project lays the foundation of an important data set focused on lethal relapsed cancers that can continue to be expanded to identify cancer-intrinsic and extrinsic microenvironmental factors supporting recurrence of pediatric brain tumors.

Study EGAS50000001288

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

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

Study EGAS50000001181

Genomic landscapes of endometrioid and mucinous ovarian cancers and morphologically similar tumor types

While endometrioid and mucinous ovarian cancers represent nearly a quarter of ovarian cancers, their molecular characteristics and pathologic origins are poorly understood. To characterize the genomic and epigenomic alterations of these ovarian cancer subtypes and evaluate links to morphologically similar tumors from other sites, we performed sequence, copy number, mutation signature, rearrangement, and methylation analyses from tumor samples and matched normal tissues of 133 patients. These analyses included patients with ovarian endometrioid (n=44), ovarian mucinous (n=43), uterine endometrial (n=15), and gastrointestinal mucinous cancers (n=31). In addition to identifying genes previously known to be involved in these tumors, we identified alterations in RAD51C, NOTCH4, SMARCA1/4 and JAK1 in ovarian endometrioid, ESR1 in uterine endometrioid, and SMARCA4 in ovarian mucinous cancers. Whole genome sequencing revealed rearrangements involving PTEN, NF1, and NF2 in ovarian endometrioid tumors, and NF1 and MED1 in ovarian mucinous cancers. The number of alterations and affected genes in ovarian and uterine endometrioid cancers were similar as were genome-wide methylation profiles, suggesting that these tumors may have a similar tissue of origin. By contrast, methylation patterns in ovarian mucinous cancers were more similar to other gastrointestinal mucinous cancers, including mucinous cancers of the stomach, colon, and pancreas. These analyses provide insights into the genomic landscapes and origins of mucinous and endometrioid ovarian carcinomas, providing new avenues for early clinical intervention and management of patients with these cancers.

Study EGAS50000000523

Systematic Growth Factor Profiling Platform for 3D Tumor Models Reveals EstradiolResponsive Cellular Mechanisms of Immunotherapy Resistance

Epithelial ovarian cancer (EOC) exhibits significant molecular heterogeneity, which is thought to be influenced by the tumor microenvironment (TME). However, the specific factors driving this heterogeneity remain poorly understood. In this study, we integrated growth factoromics and genomic data to identify TME factors contributing to EOC heterogeneity and immunosuppression. Our analysis revealed distinct subgroups of EOC samples with differential responsiveness to specific growth factors, particularly estradiol and R-spondin. Notably, samples with high responsiveness to estradiol and Wnt (C1 cluster) were enriched in the Mesenchymal (MES) molecular subtype, while samples with high responsiveness to R-spondin (C2 cluster) were associated with the Differentiated (DIF) and Immunoreactive (IMR) subtypes. We discovered that the G protein-coupled estrogen receptor (GPER), rather than the classical estrogen receptor alpha (ERα), was the key mediator of estrogen responsiveness and a determinant of molecular subtype classification. Single-cell sequencing analysis identified subsets of estrogen-responsive malignant cell (MAL.PDCD5) and cancer-associated fibroblast (FB.TNFSF10), which were associated with immunosuppressive patterns in the TCGA cohort. Importantly, the interaction between these cell populations significantly increased upon estrogen treatment, suggesting a potential mechanism by which estrogen signaling may contribute to the development of an immunosuppressive TME. Our findings highlight the critical role of estrogen signaling in shaping the molecular heterogeneity and immunosuppressive TME of EOC and provide novel insights into potential therapeutic targets for overcoming immunosuppression in this malignancy.

Study EGAS50000000422

Tumor-reactive heterotypic CD8 T cell clusters from clinical samples

Emerging evidence suggests a correlation between CD8+ T cell-tumor cell proximity and anti-tumor immune response. However, it remains unknown whether these cells exist as functional clusters that can be isolated from clinical samples. By conventional and imaging flow cytometry, we show here that from 21/21 human melanoma metastases, we could isolate heterotypic clusters, comprising CD8+ T cells interacting with one or more tumor cells and/or antigen-presenting cells (APCs). Single-cell RNA-sequencing revealed that T cells from clusters were enriched for gene signatures associated with tumor reactivity and exhaustion. Clustered T cells exhibited increased TCR clonality indicative of expansion, while TCR-matched T cells showed more exhaustion and co-modulation when conjugated to APCs than to tumor cells. T cells that were expanded from clusters ex vivo exerted on average 9-fold increased killing activity towards autologous melanomas, which was accompanied by enhanced cytokine production. Upon adoptive cell transfer into mice, T cells from clusters showed improved patient-derived melanoma control, which was associated with increased T cell infiltration and activation. Together, these results demonstrate that tumor-reactive CD8+ T cells are enriched in functional clusters with tumor cells and/or APCs and that they can be isolated and expanded from clinical samples. Typically excluded by single cell gating in flow cytometry, these distinct heterotypic T cell clusters serve as a valuable source to decipher functional tumor-immune cell interactions, while they may also be therapeutically explored.

Study EGAS50000000785

The circulating cell-free DNA landscape in sepsis is dominated by impaired liver clearance

Circulating cell-free DNA (cfDNA) is a promising molecular biomarker. However, its utility in severe infection remains poorly understood. Here, we isolated cfDNA from sepsis patients and controls, demonstrating a 41-fold increase in the amount of cfDNA in circulation during disease. We used sequencing to reconstruct cfDNA methylomes, fragmentation profiles, and nucleosome footprints across 58 samples. We observed no difference in cfDNA composition between patients and controls, challenging the idea that cfDNA increases due to higher immune cell death during sepsis. Instead, we suggest that liver dysfunction prevents efficient clearance of cfDNA during disease. This was supported by fragmentation and end-motif patterns, both of which showed evidence of cfDNA being exposed to circulating nucleases for a prolonged period, proportionally to the extent of liver dysfunction. Variation in cfDNA of megakaryocyte-erythroid progenitor origin was also a significant contributor in sepsis, increasing over time. Moreover, we showed that cfDNA retains nucleosome footprints with cell type-specific gene activity information. We developed a novel approach to study nucleosome phasing that successfully recovers tissue-specific signatures. By combining this with single-cell data, we demonstrated that sepsis patients with liver dysfunction have higher amounts of cfDNA derived from Kupffer cells and the liver parenchyma. In conclusion, we present the first high-throughput multi-modal study of cfDNA during sepsis, which will serve as a reference point for future studies on the role of this biomarker in critical illness.

Study EGAS50000001033

Relapse CHL study

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

Study EGAS00001008222

Noninvasive prenatal molecular karyotyping from maternal plasma

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

Study EGAS00001000439

Neoadjuvant combination PD-L1 plus CTLA-4 blockade in patients with cisplatin-ineligible operable urothelial carcinoma

Immune checkpoint therapy (ICT) is being tested in the neoadjuvant setting for patients with localized urothelial carcinoma (UC), with one study reporting data in cisplatin-ineligible patients who received anti PD-L1 monotherapy. The study reported that patients with bulky tumors, a known high-risk feature defined as greater than clinical T2 disease, had fewer responses, with pathologic complete response (pCR) rate of 17%. Here, we report on the first pilot combination neoadjuvant trial (NCT02812420) with anti-PD-L1 (durvalumab) plus anti-CTLA-4 (tremelimumab) in cisplatin-ineligible patients, with all tumors identified as having high risk features (N=28). Primary endpoint was safety and we observed 6 of 28 patients (21%) with grade ≥3 immune-related adverse events, consisting of asymptomatic laboratory abnormalities (N=4), hepatitis and colitis (N=2). We also observed pCR of 37.5% and downstaging to pT1 or less in 58% of patients who completed surgery (N=24). In summary, we provide initial safety, efficacy and biomarker data with neoadjuvant combination anti-PD-L1 plus anti-CTLA-4, which warrants further development for patients with localized UC, especially cisplatin-ineligible patients with high-risk features who do not currently have an established standard-of-care neoadjuvant treatment. Publication: Gao et al. Neoadjuvant PD-L1 plus CTLA-4 blockade in patients with cisplatin-ineligible operable high-risk urothelial carcinoma. Nature Medicine volume 26, pages1845–1851(2020)

Study EGAS00001004074

Pseudodiastrophic dysplasia expands the known phenotypic spectrum of defects in proteoglycan biosynthesis

Pseudodiastrophic dysplasia (PDD) is a severe skeletal dysplasia associated with prenatal manifestation and early lethality. Clinically, PDD is classified as a ‘dysplasia with multiple joint dislocations’ however, the molecular aetiology of the disorder is currently unknown. In order to identify the genetic defects underlying PDD, we performed whole exome sequencing (WES) on 3 patients from 2 unrelated families, clinically diagnosed with PDD. WES resulted in the identification of bi-allelic variants in the established skeletal dysplasia genes, B3GAT3 (Family 1) and CANT1 (Family 2). Mutations in these genes have previously been reported to cause ‘multiple joint dislocations, short stature, and craniofacial dysmorphism with or without congenital heart defects’ (B3GAT3) and Desbuquois dysplasia 1 (CANT1); disorders in the same nosological group as PDD. Follow-up of the B3GAT3 variants demonstrated significantly reduced B3GAT3/GlcAT-I expression. Downstream in vitro functional analysis revealed abolished biosynthesis of glycosaminoglycan side chains on proteoglycans. Functional evaluation of the CANT1 variant showed impaired nucleotidase activity, which results in inhibition of glycosaminoglycan synthesis through accumulation of uridine diphosphate. For the families described in this study, the PDD phenotype was caused by mutations in the known skeletal dysplasia genes B3GAT3 and CANT1, demonstrating the advantage of genomic analyses in delineating the molecular diagnosis of skeletal dysplasias. This finding expands the phenotypic spectrum of B3GAT3- and CANT1- related skeletal dysplasias to include PDD, and highlights the significant phenotypic overlap of conditions within the proteoglycan biosynthesis pathway.

Study EGAS00001004145

Targeted sequencing of candidate regions on chromosome 22q predisposing to multiple schwannomas

Constitutional LZTR1 or SMARCB1 pathogenic variants have been found in ~86% of familial and ~40% of sporadic schwannomatosis cases. Hence, we performed massively parallel sequencing of the entire LZTR1, SMARCB1 and NF2 genomic loci in 35 individuals with schwannomas negative for constitutional first-hit pathogenic variants in the LZTR1/SMARCB1/NF2 coding sequences, however with 22q deletion and/or a different NF2 pathogenic variant in each tumor, including six cases with only one tumor available. Furthermore, we verified whether any other LZTR1/SMARCB1/NF2 (likely) pathogenic variant(s) could be found in 16 cases carrying a SMARCB1 constitutional variant in the 3’ untranslated region (3’-UTR) c.*17C>T, c.*70C>T or c.*82C>T. As no additional variants were found, functional studies were performed to clarify the effect of these 3’-UTR variants on the transcript. The 3’-UTR variants c.*17C>T and c.*82C>T showed pathogenicity by negatively affecting SMARCB1 transcript level. Two novel deep intronic SMARCB1 variants, c.500+883T>G and c.500+887G>A, resulting in out-of-frame missplicing of intron 4, were identified in two unrelated individuals. Further resequencing of the entire repeat - masked genomics sequences of chromosome 22q in individuals negative for pathogenic variants in the SMARCB1/LZTR1/NF2 coding- and non-coding regions revealed five potential schwannomatosis-predisposing candidate genes, i.e. MYO18B, NEFH, SGSM1, SGSM3 and SBF1, pending further verification.

Study EGAS00001005680

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

Evolution of the African pygmy phenotype

The evolutionary history of the human ‘pygmy’ phenotype (small body size), a characteristic of African and Southeast Asian rainforest hunter-gatherers, is largely unknown. Here we use a genome-wide admixture mapping analysis to identify 16 genomic regions that are significantly associated with the pygmy phenotype in the Batwa, a rainforest hunter-gatherer population from Uganda (East Central Africa). The identified genomic regions have multiple attributes that provide supporting evidence of genuine association with the pygmy phenotype, including enrichments for SNPs previously associated with stature variation in Europeans and for genes with growth hormone receptor and regulation functions. To test adaptive evolutionary hypotheses we computed the haplotype-based iHS statistic and the level of population differentiation (FST) between the Batwa and their agricultural neighbors, the Bakiga, for each SNP. Both |iHS| and FST values were significantly higher for SNPs within the Batwa pygmy phenotype-associated regions than the remainder of the genome, a signature of polygenic adaptation. In contrast, when we expanded our analysis to include Baka rainforest hunter-gatherers from Cameroon and Gabon (West Central Africa) and Nzebi and Nzime neighboring agriculturalists, we did not observe elevated |iHS| or FST values in these genomic regions. Together, these results suggest adaptive and at least partially convergent origins of the pygmy phenotype even within Africa, supporting the hypothesis that small body size confers a selective advantage for tropical rainforest hunter-gatherers but raising questions about the antiquity of this behavior.

Study EGAS00001000908

Illumina HumanCoreExome genotyping data from the TEENAGE (TEENs of Attica: Genes and Environment) study.

The TEENAGE study is a cross-sectional study comprising healthy Greek adolescents. The study was conducted according to the guidelines specified in the Declaration of Helsinki and all procedures involving human subjects were approved by the Institutional Review Board of Harokopio University and the Greek Ministry of Education, Lifelong Learning and Religious Affairs. The TEENAGE study target population comprised adolescent students aged 13 to 15 years attending the first three classes of public secondary schools located in the wider Athens area of Attica. The regions across Attica from which students were recruited were physically similar and all urban in nature. Sampling of regions and schools selected to participate was random. Students were enrolled from the wider urban region to ensure that they were coming from families with varying educational background and socio-economic status. Moreover, in order to exclude any vast economic disparities, all participating schools were exclusively public schools, which are free of charge. A total of twenty-four schools participated in the study out of 285 junior secondary public schools in Attica. Anthropometric, medical and dietary information was obtained. DNA samples were genotyped using Illumina HumanCoreExome BeadChips (Illumina, San Diego, CA, USA) at the Wellcome Trust Sanger Institute, Hinxton, UK. Genotypes were called using GenCall (Illumina Genome Studio) followed by zCall and quality control (QC) was performed in two stages (pre- and post-zCall).

Study EGAS00001001001

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

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

Study EGAS00001006596

Whole exome sequencing of an invasive diffuse intrinsic pontine glioma sampled from different sites

The lateral ventricle subventricular zone (SVZ) is a frequent site of high-grade glioma (HGG) spread and tumor invasion of the SVZ is a predictor of a worse clinical prognosis. A range of HGG types invade the stem cell niche of the SVZ, including both adult glioblastoma and pediatric high-grade gliomas such as diffuse intrinsic pontine glioma (DIPG). The cellular and molecular mechanisms mediating this frequent invasion of the SVZ are poorly understood. Here we demonstrate that neural precursor cell (NPC):glioma cell communication underpins the propensity of glioma to colonize the SVZ. SVZ NPCs secrete chemoattractant signals toward which glioma cells home. Biochemical and proteomic analyses of factors secreted by SVZ NPCs revealed a discrete list of candidate proteins. Necessity and sufficiency testing implicated the neurite outgrowth-promoting factor pleiotrophin, along with three required binding partners (secreted protein acidic and rich in cysteine (SPARC), SPARC-like protein 1, and heat shock protein 90B), as key mediators of this chemoattractant effect. Pleiotrophin protein expression is strongly enriched in the SVZ, and knockdown of pleiotrophin expression starkly reduced glioma invasion of the SVZ in the adult murine brain. Pleiotrophin, in complex with the three binding partners, activated the Rho/ROCK pathway in DIPG cells, and inhibition of ROCK resulted in decreased DIPG invasion toward SVZ neural precursor cell-secreted factors. These findings demonstrate a pathogenic role for neural precursor cell – glioma interactions and potential therapeutic targets to limit glioma invasion.

Study EGAS00001002326

Platelet RNAseq data for SLFN14 K219N patients

Pathogenic missense variants in SLFN14, which encodes an RNA endoribonuclease protein that regulates rRNA degradation, are known to cause inherited thrombocytopenia with impaired platelet aggregation and ATP secretion. Despite rather mild laboratory defects, these patients display an obvious bleeding phenotype. The function of SLFN14 in megakaryocyte (MK) and platelet biology is unknown. This study aims to characterize the platelet transcriptome in patients with a SLFN14 K219N variant and model the disease in the immortalized megakaryocyte cell line imMKCL. Total platelet RNA was sequenced for two patients and 19 healthy controls. Differential gene expression analysis yielded a total of 2999 and 2888 significantly (|log2FC|>1, FDR<0.05) up- and downregulated genes, respectively. Remarkably, these downregulated genes were not enriched for any biological pathway while upregulated genes were enriched for pathways involved in (mitochondrial) translation and transcription with a significant upregulation of 134 ribosomal protein genes (RPG). Heterozygous and homozygous SLFN14 K219N imMKCL showed a defect in MK differentiation that aggravated during cell passaging of undifferentiated cells and in proplatelet formation. SLFN14 defective platelets and MK showed signs of rRNA degradation while this was absent in undifferentiated imMKCL cells. Upregulation of RPG through increased mTOR signaling in SLFN14 K219N MK seems to be a compensatory response for rRNA degradation. Indeed, mTOR inhibition with rapamycin resulted in further enhanced rRNA degradation in SLFN14 K219N MK. Taken together, our study indicates dysregulation of mTOR-regulated ribosomal biogenesis as the disease mechanism for SLFN14-related thrombocytopenia.

Study EGAS00001006339

Plasma-Seq of patients with metastatic prostate cancer

Study 1 In this study we analysed patients with metastatic prostate cancer to scan their tumor genomes noninvasively in plasma DNA. We wanted to make whole-genome analysis from plasma DNA amenable to clinical routine applications and developed an approach based on a benchtop high-throughput platform, i.e. Illuminas MiSeq instrument. We performed whole-genome sequencing from plasma at a shallow sequencing depth to establish a genome-wide copy number profile of the tumor at low costs within 2 days. The genome-wide profiling in the plasma of our patients revealed multiple copy number aberrations including those previously reported in prostate tumors, such as losses in 8p and gains in 8q. High-level copy number gains in the AR locus were observed in patients with CRPC but not with CSPC disease. We identified the TMPRSS2-ERG rearrangement associated 3-Mbp deletion on chromosome 21 and found corresponding fusion plasma fragments in these cases. In an index case multiregional sequencing of the primary tumor identified different copy number changes in each sector, suggesting multifocal disease. Our plasma analyses of this index case, performed 13 years after resection of the primary tumor, revealed novel chromosomal rearrangements, which were stable in serial plasma analyses over a 9 months period, which is consistent with the presence of one metastatic clone. Our approach provides specific genomic signatures within 2 days which may therefore serve as “liquid biopsy”.

Study EGAS00001000451

RB1 Gene Inactivation by Chromothripsis in Human Retinoblastoma

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

Study EGAS00001000598

Evolutionary origin and methylation status of human intronic CpG islands that are not present in mouse

Imprinting of the human RB1 gene is due to the presence of a differentially methylated CpG island (CGI) in intron 2, which is part of a retrogene derived from the PPP1R26 gene on chromosome 9. The murine Rb1 gene does not have this retrogene and is not imprinted. We have investigated whether the RB1/Rb1 locus is unique with respect to these differences. For this we have compared the CGIs from human and mouse by in silico analyses and found out that there are at least 2033 human intronic CGIs that are not present in the mouse. Among these CGIs, 104 show sequences similarities elsewhere in the human genome, which suggests that they arose by retrotransposition. We could narrow down the time points when most of these CGIs appeared during evolution. Their methylation status was analysed in two monocyte methylome datasets from whole genome bisulfite sequencing and in 18 published methylomes. Five CGIs appear to be differentially methylated, which suggests that they are subject of genomic imprinting. They are located in introns of the following genes: RB1, ASRGL1, PARP11, PDXDC1 and MYO1D. Four of these CGIs are parts of retrogenes. Our study supports the notion that imprinting builds on host defence mechanisms by which the genome protects itself against retrotransposed and foreign DNA and that the epigenetic fate of the retrotransposed DNA depends on selective forces at the integration site.

Study EGAS00001000719

Open_Targets_020_Epigenomes_of_Cell_Lines

There is currently a drive to establish cell based assay systems of greater human biological and disease relevance through the use of well characterised transformed cell lines, primary cells and complex cellular models (e.g. co-culture, 3D models). However, although the field is gaining valuable experience in running more non-standard & complex cell assays for target validation and compound pharmacology studies, there is the lack of a systematic approach to determine if this expansion in cell assay models is reflected in increased human biological and disease relevance. The increasing wealth of publically available transcriptomic, and epigenome (ENCODE and Epigenome Roadmap) data represents an ideal reference mechanism for determining the relationship between cell types used for target & compound studies to primary human cells and tissues from both healthy volunteers & patients. The CTTV020 epigenomes of cell line project aims to generate epigenetic and transcriptomic profiles of cell lines and compare these with existing and newly generated reference data sets from human tissue and cell types. The aim is to identify assay systems which will provide greater confidence in translating target biology and compound pharmacology to patients. Multiple cell types commonly used within research have been grouped according to biology. Examples include erythroid, lung epithelial, hepatocyte cell types and immortalised models of monocyte / macrophage biology. This data is part of a pre-publication release. For information on the proper use of pre-publication data shared by the Wellcome Trust Sanger Institute please see http://www.sanger.ac.uk/datasharing/

Study EGAS00001003138

Tissue-specific mutation accumulation in human adult stem cells during life

Gradual accumulation of mutations in human adult stem cells during life is associated with various age-related diseases, including cancer. The number of stem cell divisions throughout life is believed to be a major determinant for mutation accumulation and could explain the extreme variation of cancer incidence across different organs. Yet, mutation patterns and rates of healthy adult stem cells remain unknown. Here, we determined genome-wide mutation patterns in primary adult stem cells of the small intestine, colon and liver of human donors with ages ranging from 3 to 87 years. We find that the number of mutations increases linearly with age up to several thousand mutations per cell at 87 years of age, while mutation spectra remain constant throughout life. Small intestine and colon stem cells have a 2-fold higher mutation rate per year compared with liver stem cells. These differences could be exclusively attributed to the mutagenic action of spontaneous deamination of cytosine residues and may reflect the high stem cell division rate in these tissues. The genomic distribution of somatic mutations is non-random and predominantly associated with DNA replication dynamics in the small intestine and colon, and with transcription in the liver. These results indicate that a stable balance between various mutagenic and DNA repair processes is maintained throughout life and that the activity of these processes in adult stem cells varies between tissues.

Study EGAS00001001682

Cistrome-partitioning reveals convergence of somatic mutations and risk-variants on master transcription regulators in primary prostate tumors

Thousands of noncoding somatic Single Nucleotide Variants (SNVs) of unknown function are reported in tumors. Partitioning the genome according to cistromes, reveals the enrichment of somatic SNVs in prostate tumors as opposed to adjacent normal tissue cistromes of master transcription regulators, including AR, FOXA1 and HOXB13. This parallels enrichment of prostate cancer genetic predispositions over these transcription regulators’ tumor cistromes, exemplified at the 8q24 locus harboring both risk-variants and somatic SNVs in cis-regulatory elements, upregulating MYC expression and altering the binding of transcription regulators to DNA. However, Massively-Parallel Reporter Assays reveal that few SNVs can alter the transactivation potential of individual CREs. Instead, SNVs accumulate, similarly to inherited riskvariants, in cistromes of master transcription regulators required for prostate cancer development. Difficulties in inferring the biological significance of noncoding mutations have limited their inclusion in precision genomics medicine pipelines. Most attempts to delineate a role for noncoding mutations relied on detecting evidence for positive selection within individual CREs, such as reported for the TERT gene promoter. By considering the enrichment of noncoding mutations in cistromes as opposed to individual CREs, we reveal their specificity towards master transcription regulators that promote prostate cancer development, a feature shared with inherited risk-variants. Overall, our work provides a blueprint for the functional interpretation of noncoding mutations in genomic tests relying on defining cis-regulatory units according to cistrome-partitioning to identify cancer driver transcription regulators.

Study EGAS00001003928

Characterization of HCV-specific CD4 T cells during DAA-Therapy

Background: Chronic HCV-infection is characterized by a severe impairment of HCV-specific CD4 T cell help that is driven by chronic antigen stimulation. We aimed to study the fate of HCV-specific CD4 T cells after viral elimination. Methods: HCV-specific CD4 T cell responses were longitudinally analyzed using MHC class II tetramer-technology, multicolor flow cytometry and RNA sequencing in a cohort of chronically HCV-infected patients undergoing therapy with direct-acting antivirals. In addition, HCV-specific neutralizing antibodies and CXCL13 levels were analyzed. Results: We observed that the frequency of HCV-specific CD4 T cells increased within two weeks after initiation of DAA therapy. Multicolor flow cytometry revealed a downregulation of exhaustion and activation markers and an upregulation of memory-associated markers. While cells with a Th1 phenotype were the predominant subset at baseline, cells with phenotypic and transcriptional characteristics of follicular T helper cells increasingly shaped the circulating HCV-specific CD4 T cell repertoire, suggesting antigen-independent survival of this subset. These changes were accompanied by a decline of HCV-specific neutralizing antibodies and the germinal center activity. Conclusion: We identified a population of HCV-specific CD4 T cells with a follicular T helper cell signature that is maintained after therapy-induced elimination of persistent infection and may constitute an important target population for vaccination efforts to prevent re-infection and immunotherapeutic approaches for persistent viral infections.

Study EGAS00001003950

Molecular Attributes Underlying Central Nervous System and Systemic Relapse in Diffuse Large B-cell Lymphoma

Central nervous system (CNS) relapse in diffuse large B-cell lymphoma remains a significant conundrum. While clinical risk scores and biomarkers allow the identification of patients with higher risk of CNS relapse, their predictive ability is low. Herein, we report findings from a comprehensive study examining gene expression changes, mutations and clonal evolution. Firstly, considering diagnostic samples from patients with CNS relapse, systemic relapse or no relapse, we found a striking enrichment of gene expression changes reflecting protein biosynthesis in patients who experienced CNS or systemic relapse. Furthermore, diagnostic specimens of relapsing cases were characterized by low immune infiltrates. The most discriminative gene feature in our analysis was SLAMF1, high expression of which was associated with favourable outcomes. Next, we assessed whether specific gene mutations were associated with subsequent CNS or systemic relapse. Activated B-cell-like genes such as PIM1, CD79B and PRDM1 were more frequently mutated in cases with CNS relapse, whereas TP53, BCL2 and FOXO1 mutations were associated with systemic relapse. Lastly, we performed tumor evolution analysis in five patients, identifying significant clonal divergence that accompanied CNS relapse. In summary, we identified both common biology underlying both CNS and systemic relapse, independent of cell-of-origin (COO), and specific COO-enriched gene mutations that conferred a higher risk of either CNS or systemic relapse. Moreover, we identified a pattern of clonal evolution that can explain difficulties in both predicting and treating CNS relapse.

Study EGAS00001004057

UK10K_RARE_HYPERCHOL

In the UK10K project we propose a series of complementary genetic approaches to find new low frequency/rare variants contributing to disease phenotypes. These will be based on obtaining the genome wide sequence of 4000 samples from the TwinsUK and ALSPAC cohorts (at 6x sequence coverage), and the exome sequence (protein coding regions and related conserved sequence) of 6000 samples selected for extreme phenotypes. Our studies will focus primarily on cardiovascular-related quantitative traits, obesity and related metabolic traits, neurodevelopmental disorders and a limited number of extreme clinical phenotypes that will provide proof-of-concept for future familial trait sequencing. We will analyse directly quantitative traits in the cohorts and the selected traits in the extreme samples, and also use imputation down to 0.1% allele frequency to extend the analyses to further sample sets with genome wide genotype data. In each case we will investigate indels and larger structural variants as well as SNPs, and use statistical methods that combine rare variants in a locus or pathway as well as single-variant approaches.Familial Hypercholesterolemia is a condition where the affected person has a consistently high level of LDL which can lead to early clogging of the coronary arteries. All patients selected for this study will have been found not to carry the common APOB and PCSK9 mutations, and to have no detectable LDLR mutations by testing for 18 common mutations. For further information with regard to this cohort please contact Steve Humphries(steve.humphries@ucl.ac.uk).

Study EGAS00001000129

Integration of genomics and histology reveals diagnosis and effective therapy of refractory cancer of unknown primary with PDL1 amplification (H021)

Identification of the tissue of origin in cancer of unknown primary (CUP) poses a diagnostic challenge and is critical for directing site-specific therapy. Currently, clinical decision making in patients with CUP primarily relies on histopathology and clinical features. Comprehensive molecular profiling has the potential to contribute to diagnostic categorization and, most importantly, guide CUP therapy through identification of actionable lesions. We here report the case of an advanced-stage malignancy initially mimicking poorly differentiated soft-tissue sarcoma that did not respond to multi-agent chemotherapy. Molecular profiling within a clinical whole-exome and transcriptome sequencing program revealed a heterozygous, highly amplified KRAS G12S mutation, compound-heterozygous TP53 mutation/deletion, high mutational load, and focal amplification of chromosomes 9p (including PDL1 [CD274] and JAK2) and 10p (including GATA3). Integrated analysis of molecular data and conventional histopathology suggested a diagnosis of triple-negative breast cancer (TNBC) and provided a rationale for immune checkpoint inhibitor therapy with pembrolizumab, which resulted in rapid clinical improvement and a lasting partial remission. Analysis of 157 TNBC samples from The Cancer Genome Atlas revealed focal, high-level PDL1 amplification coinciding with excessive PDL1 mRNA expression in 24% of cases. Collectively, these results illustrate the diagnostic utility of multidimensional tumor profiling in cases with non-descript histology and immune phenotype, demonstrate the predictive power of genomic PDL1 amplification for immune checkpoint inhibition, and suggest a targeted therapeutic strategy in chromosome 9p24.1/PDL1-amplified cancers.H021

Study EGAS00001001846

Persistent STAG2 mutation in recurrent pediatric glioblastoma

Similar to their adult counterparts, the prognosis for pediatric patients with high-grade gliomas remains poor. At time of recurrence, treatment options are limited and remain without consensus. This report describes the genetic findings, obtained from whole exome sequencing of a pediatric patient with glioblastoma who underwent multiple surgical resections and treatment with standard chemoradiation, as well as a novel recombinant poliovirus vaccine therapy. Strikingly, despite the variety of treatments, there was persistence of a tumor clone, characterized by a deleterious STAG2 mutation, whose deficiency in pre-clinical studies can cause aneuploidy and aberrant mitotic progression, but remains understudied in the clinical setting. There was near elimination of an EGFR mutated and amplified tumor clone after gross total resection, standard chemoradiation, and poliovirus therapy, followed by the emergence of a persistently STAG2 mutated clone, with rare mutations in PTPN11 and BRAF, the latter composed of a novel deleterious mutation previously not reported in pediatric glioblastoma (p.D594G). This was accompanied by a mutation signature shift towards one characterized by increased DNA damage repair defects, consistent with the known underlying STAG2 deficiency. As such, this case represents a novel report following the clinical and genetic progression of a STAG2 mutated glioblastoma, including treatment with a novel and emerging immunotherapy. Although STAG2 deficiency comprises only a small subset of gliomas, this case adds clinical evidence to existing pre-clinical data supporting a role for STAG2 mutations in tumor formation and resistance to standard therapies.

Study EGAS00001004340

Negligible impact on missing heritability of autoimmune-locus rare coding-region variants

Genome wide association studies (GWAS) have identified common variants of modest effect size at hundreds of loci for common autoimmune diseases - however a substantial fraction of heritability remains unexplained, and to which rare variants may contribute. To discover rare variants and test them for association with a phenotype, the majority of studies re-sequence a small initial sample size and then genotype the discovered variants in a larger sample set. This approach will fail to analyse a large fraction of the rare variants present in the entire sample set. Here we perform simultaneous amplicon sequencing-based variant discovery and genotyping for coding exons of 25 GWAS risk genes in 41,911 white-European origin UK subjects comprising 24,892 subjects with six autoimmune disease phenotypes and 17,019 controls, and show that rare coding-region variants at known loci play a negligible role in common autoimmune disease susceptibility. These results do not support the rare variant - synthetic genome-wide association hypothesis. Many known autoimmune disease risk loci contain multiple independent common and low frequency variant association signals, and so genes in these loci are a priori stronger candidates to harbor rare coding-region variants than other genes. Our data suggest that the missing heritability for common autoimmune diseases may not lie in the rare coding-region variant portion of the allelic spectrum, but perhaps as others have proposed in very many common variant signals of weak effect.

Study EGAS00001000476

Massively parallel nanowell-based single-cell gene expression profiling

Technological advances have enabled transcriptome characterization of cell types at the single-cell level providing new biological insights. New methods that enable simple yet high-throughput single-cell expression profiling are highly desirable. Here we report a novel nanowell-based single-cell RNA sequencing system, ICELL8, which enables processing of thousands of cells per sample. The system employs a 5184-nanowell-containing microchip to capture ~1300 single cells and process them. Each nanowell contains preprinted oligonucleotides encoding poly-d(T), a unique well barcode, and a unique molecular identifier. The ICELL8 system uses imaging software to identify nanowells containing viable single cells and only wells with single cells are processed into sequencing libraries. Here, we report the performance and utility of ICELL8 using samples of increasing complexity from cultured cells to mouse solid tissue samples. Our assessment of the system to discriminate between mixed human and mouse cells showed that ICELL8 has a low cell multiplet rate (< 3%) and low cross-cell contamination. We characterized single-cell transcriptomes of more than a thousand cultured human and mouse cells as well as 468 mouse pancreatic islets cells. We were able to identify distinct cell types in pancreatic islets, including alpha, beta, delta and gamma cells. Overall, ICELL8 provides efficient and cost-effective single-cell expression profiling of thousands of cells, allowing researchers to decipher single-cell transcriptomes within complex biological samples.

Study EGAS00001002320

Genome of the Netherlands

The Genome of the Netherlands (GoNL) Project characterizes DNA sequence variation, common and rare, for SNVs, short insertions and deletions (indels) and larger deletions in 769 individuals of Dutch ancestry selected from five biobanks under the auspices of the Dutch hub of the Biobanking and Biomolecular Research Infrastructure (BBMRI-NL). The samples come from a representative sample of 250 trio-families from all provinces in the Netherlands. The parent-offspring trios include adult individuals ranging in age from 19 to 87 years (mean=53 years; SD=16 years) from birth cohorts 1910-1994. Sequencing was done on blood-derived DNA from uncultured cells and accomplished coverage was 14-15x. Samples where contributed by LifeLines (http://lifelines.nl/lifelines-research/general), The Leiden Longevity Study (http://www.healthy-ageing.nl; http://www.langleven.net), The Netherlands Twin Registry (NTR: http://www.tweelingenregister. org), The Rotterdam studies, (http://www.erasmus-epidemiology.nl/rotterdamstudy) and the Genetic Research in Isolated Populations program (http://www.epib.nl/research/geneticepi/research.html#gip). The sequencing was carried out in collaboration with the Beijing Institute for Genomics (BGI). The analysis was done by a consortium lead by UMCG, LUMC, Erasmus MC, VU university and UMCU, see http://www.nlgenome.nl. Funding for the project was provided by the Netherlands Organization for Scientific Research under award number 184021007, dated July 9, 2009 and made available as a Rainbow Project of the Biobanking and Biomolecular Research Infrastructure Netherlands (BBMRI-NL).

Study EGAS00001000644

Copy number profiling of putative cancer stem from pleural effusion aspirates from breast cancer patients

A common symptom during late stage breast cancer disease is pleural effusion, which is related to poor prognosis. In pleural effusions malignant cells can be detected, which indicates a metastatic spread from the primary tumor site. Beside the usefulness for investigation of the metastatic process, it has been a source to isolate cells with putative cancer stem cell (CSC) properties. For this study, pleural effusion aspirate from 25 metastatic breast cancer patients were processed to establish in vitro culture. Patient-derived aspirates were cultured under sphere forming conditions and isolated primary cultures were sorted for cancer stem cell populations ALDH1+ and CD44+CD24-/low and sphere forming efficiency of CSC and non-CSC subpopulation was determined. Unsorted cells had a higher potential to form spheres compared to CSC subpopulations. Subsequently, DNA was isolated from pleural effusions before and after cell sorting, to compare their copy number profiles generated by low-coverage whole genome sequencing to DNA from FFPE primary tumors or bone metastasis. Unbalanced profiles could be detected in four pleural samples before sorting, in two ALDH1+ sorted samples, in one CD44+CD24-/low and in all control samples. Additionally, spheroid culture from one pleural sample could be passaged on and was characterized. Therefore, pleural effusions are an appropriate media for enrichment of putative CSC. However, low-coverage sequencing is difficult due to low quantities of abnormal cells compared to normal cells.

Study EGAS00001002343

Cisplatin increases sensitivity to FGFR inhibition in patient-derived xenograft models of lung squamous cell carcinoma

Lung squamous cell carcinoma (SqCC) is a molecularly complex and genomically unstable disease. No targeted therapy is currently approved for lung SqCC, although potential oncogenic drivers of SqCC have been identified, including amplification of the fibroblast growth factor receptor 1 (FGFR1). Reports from a recently completed clinical trial indicate low response rates in patients treated with FGFR tyrosine kinase inhibitors, suggesting inadequacy of FGFR1 amplification as a biomarker of response, or the need for combination treatment. We aimed to develop accurate models of lung SqCC and determine improved targeted therapies for these tumors. We show that detection of FGFR1 mRNA by RNA in situ hybridisation is a better predictor of response to FGFR inhibition than FGFR1 gene amplification using clinically relevant patient-derived xenograft (PDX) models of lung SqCC. FGFR1-overexpressing tumors were observed in all histological subtypes of non-small cell lung cancers (NSCLC) as assessed on a tissue microarray, indicating a broader range of tumors that may respond to FGFR inhibitors. In FGFR1-overexpressing PDX tumors, we observed increased differentiation and reduced proliferation following FGFR inhibition. Combination therapy with cisplatin was able to increase tumor cell death, and dramatically prolonged animal survival compared to single agent treatment. Our data suggest that FGFR tyrosine kinase inhibitors can benefit NSCLC patients with FGFR1-overexpressing tumors and provides a rationale for clinical trials combining cisplatin with FGFR inhibitors.

Study EGAS00001002423

Identification_of_synthetic_lethal_genes_by_CRISPR_Cas9_library

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

Study EGAS00001002117

The evolutionary history of human colitis-associated colorectal cancer

Inflammatory bowel disease confers an increased lifetime risk of developing colorectal cancer (CRC), and colitis-associated CRC (CA-CRC) is molecularly distinct from sporadic CRC (S-CRC). Here we have dissected the evolutionary history of CA-CRC using multi-region sequencing. Exome sequencing was performed on fresh-frozen multiple regions of carcinoma, adjacent non-cancerous mucosa and blood from 12 CA-CRC patients (n=55 exomes), and key variants were validated with orthogonal methods. Genome-wide copy number profiling was performed using SNP arrays and low-pass whole genome sequencing on archival non-dysplastic mucosa (n=9), low-grade dysplasia (LGD, n=30), high-grade dysplasia (HGD, n=13), mixed LGD/HGD (n=7) and CA-CRC (n=19). Phylogenetic trees were reconstructed, and evolutionary analysis used to reveal the temporal sequence of events leading to CA-CRC. 10/12 tumors were microsatellite stable with a median mutation burden of 3.0 single nucleotide alterations (SNAs) per Mb, around 20% higher than S-CRC (2.5 SNAs/Mb), and consistent with elevated aging-associated mutational processes. Non-dysplastic mucosa had considerable mutation burden (median 47 SNAs), including mutations shared with the neighboring CA-CRC, indicating a pre-cancer mutational field. CA-CRCs were often near-triploid (42%) or near-tetraploid (21%) and phylogenetic analysis revealed that copy number alterations (CNAs) began to accrue in non-dysplastic bowel, but the LGD/HGD transition often involved a punctuated catastrophic CNA increase.

Study EGAS00001003028

Whole genome sequencing with linked reads of pediatric glioblastoma samples

Pediatric glioblastoma (pGBM) is a lethal cancer with no effective therapies. To understand mechanisms of tumor evolution in this cancer, we performed whole genome sequencing with linked reads on longitudinally resected pGBM samples. Our analyses showed that all diagnostic and recurrent samples were collections of genetically diverse subclones. Clonal composition rapidly evolved at recurrence, with less than 8% of non-synonymous single nucleotide variants being shared in diagnostic-recurrent pairs. In order to track the origins of the mutational events we observed in pGBM, we generated whole genome datasets for two patients and their parents. These trios showed that genetic variants could be (i) somatic, (ii) inherited from a healthy parent, or (iii) arose de novo in the germlines of pGBM patients. Analysis of variant allele frequencies supported a model of tumor growth involving slow-cycling cancer stem cells that give rise to fast-proliferating progenitor-like cells and to non-dividing cells. Interestingly, radiation and anti-mitotic chemotherapeutics did not increase overall tumor burden upon recurrence. These findings support an important role for slow-cycling stem cell populations in contributing to recurrences, since slow-cycling cell populations are expected to be less prone to genotoxic stress induced by these treatments and to accumulate few mutations. Our results highlight the need for new targeted treatments that account for the complex functional hierarchies and genomic heterogeneity of pGBM.

Study EGAS00001003432

Germline variants collaborate with somatic mutations and initiate and/or drive disease in primary myelodysplastic syndrome (MDS) and therapy-related myeloid neoplasms (t-MN)

Therapy-related myeloid neoplasms (T-MN) are poorly characterized secondary hematological malignancies following chemotherapy/radiotherapy exposure. We compared the clinical and mutational characteristics of T-MN (n=129) and primary myelodysplastic syndrome (P-MDS, n=108) patients. Although the somatic mutation frequency was similar between T-MN and P-MDS patients (93% in both groups), the pattern was distinct. TP53 mutations were more frequent in T-MN (29.5% versus 7%), while spliceosomal complex mutations were more common in P-MDS (56.5% versus 25.6%).In contrast to P-MDS, the ring sideroblasts (RS) phenotype was not associated with better survival in T-MN, most probably due to genetic association with TP53 mutations. SF3B1 was mutated in 96% of P-MDS with ≥15% RS, but in only 32% T-MN. TP53 mutations were detected in 92% T-MN with ≥15% RS and SF3B1 wild-type cases. Interestingly, T-MN and P-MDS patients with “Very low” or “Low” Revised International Prognostic Scoring System (IPSS-R) showed similar biological and clinical characteristics. In a Cox-regression analysis, TP53 mutation was a poor prognostic factor in T-MN, independent of IPSS-R cytogenetic, disease-modifying therapy and NRAS mutation.Our data have direct implications for T-MN management and provide evidence that, in addition to conventional disease parameters, mutational analysis should be incorporated in T-MN risk stratification.

Study EGAS00001003547

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

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

Study EGAS00001003661

Genetic and epigenetic characterization of adenoid cystic carcinoma

Translocation events are frequent in cancer and may create chimeric fusions or regulatory rearrangements that drive oncogene overexpression. Although regulatory rearrangements are increasingly recognized in hematopoietic and even solid tumors, the underlying mechanisms remain obscure. Here we identify super-enhancer translocations that drive overexpression of the oncogenic transcription factor MYB as a unifying theme in adenoid cystic carcinoma (ACC). Whole genome sequencing data and chromatin state maps for 13 primary tumors and xenografts reveal distinct chromosomal rearrangements that juxtapose super-enhancer clusters to the MYB locus. Chromosome conformation capture confirms that the translocated enhancers physically interact with the MYB promoter. Remarkably, the MYB product itself binds to the translocated enhancers, thus creating a positive feedback loop that sustains its own expression. MYB also binds a large number of active enhancers that drive opposing regulatory programs in alternate cell lineages in ACC. MYB cooperates with the transcription factor TP63 in the myoepithelial component of the tumors, but promotes a Notch program in the luminal epithelial component. Bromodomain inhibitors, which disrupt enhancer function, slow tumor growth in ACC primagraft models in vivo, but are ineffective against high grade tumors that harbor activating mutations in the Notch pathway. Thus, our study identifies super-enhancer translocations as a unifying feature of ACC, and provides insight into the mechanism by which sustained MYB overexpression drives alternate cell fates in this disease.Processed ChIP-seq data is available on GEO under accession number GSE76465.

Study EGAS00001001457

A no-stop mutation in MAGEB4 is a possible cause of rare X-linked azoospermia and oligozoospermia in a consanguineous Turkish family

Purpose: The purpose of this study was to identify mutations that cause non-syndromic male infertility using whole exome sequencing of family cases.Methods: We recruited a consanguineous Turkish family comprising nine siblings with male triplets; two of the triplets were infertile as well as one younger infertile brother. Whole exome sequencing (WES) performed on two azoospermic brothers, identified a mutation in the MAGEB4 gene which was confirmed via Sanger sequencing and then screened for on control groups and unrelated infertile subjects. The effect of the mutation on mRNA and protein levels were tested after in vitro cell transfection. Structural features of MAGEB4 were predicted throughout the conserved MAGE domain.Results: The novel single base substitution (c.1041A>T) in the X-linked MAGEB4 gene was identified as a no-stop mutation. The mutation is predicted to add 24 amino acids to the C-terminus of MAGEB4. Our functional studies were unable to detect any effect either on mRNA stability, intracellular localization of the protein or the ability to homo / heterodimerize with other MAGE proteins. We thus hypothesize that these additional amino acids may affect the proper protein interactions with MAGEB4 partners. Conclusion: The whole exome analysis of a consanguineous Turkish family revealed MAGEB4 as a possible new X-linked cause of inherited male infertility. This study provides the first clue to the physiological function of a MAGE protein.

Study EGAS00001002272

Successful BRAF/MEK-inhibition in a young patient with BRAF V600E-mutated extrapancreatic acinar cell carcinoma (HIPO-021)

Pancreatic acinar cell carcinoma (PAC) is a rare disease with a poor prognosis. Treatment options for metastatic PAC are limited and often follow chemotherapeutic regimens for pancreatic ductal adenocarcinoma. Although recurrent genomic alterations, such as BRAF fusions and defects in genes involved in homologous recombination DNA repair, have been described in PAC, data on the clinical efficacy of molecularly guided, targeted treatment are scarce. Here we describe the case of a 27-yr-old patient with BRAF V600E-mutated PAC who was successfully treated with a combination of BRAF and MEK inhibitors. The patient presented to our clinic with abdominal pain and weight loss. Imaging showed extensive retroperitoneal disease as well as mediastinal lymphadenopathy. Because of elevated alpha-fetoprotein (AFP) levels and inconclusive histologic findings, a germ cell tumor was suspected; however, PEI chemotherapy was unsuccessful. A repeat biopsy yielded the diagnosis of PAC and treatment with FOLFIRINOX was initiated. Comprehensive molecular profiling within the MASTER (Molecularly Aided Stratification for Tumor Eradication Research) precision oncology program revealed a somatic BRAF V600E mutation and a germline PALB2 stop-gain mutation. Therapy was therefore switched to BRAF/MEK inhibition, resulting in almost complete remission and disease control for 12 mo and a remarkable improvement in the patient's general condition. These results indicate that BRAF alterations are a valid therapeutic target in PAC that should be routinely assessed in this patient population.

Study EGAS00001004282

UK10K COHORT TWINSUK

The UK10K project proposes a series of complementary genetic approaches to find new low-frequency/rare variants contributing to disease phenotypes. These will be based on obtaining the genome-wide sequence of 4000 samples from the TwinsUK and ALSPAC cohorts (at 6x sequence coverage), and the exome sequence (protein-coding regions and related conserved sequence) of 6000 samples selected for extreme phenotypes. Our studies will focus primarily on cardiovascular-related quantitative traits, obesity and related metabolic traits, neurodevelopmental disorders and a limited number of extreme clinical phenotypes that will provide proof-of-concept for future familial trait sequencing. We will directly analyse quantitative traits in the cohorts and the selected traits in the extreme samples, and also use imputation down to 0.1% allele frequency to extend the analyses to further sample sets with genome wide genotype data. In each case we will investigate indels and larger structural variants as well as SNPs, and use statistical methods that combine rare variants in a locus or pathway as well as single-variant approaches. The TwinsUK samples will be part of the cohort study and will undergo whole genome sequencing. For further information with regard to this cohort please contact Brent Richards (brent.richards@kcl.ac.uk) or Nicole Soranzo (ns6@sanger.ac.uk).

Study EGAS00001000108

UK10K NEURO ASD BIONED

In the UK10K project we propose a series of complementary genetic approaches to find new low frequency/rare variants contributing to disease phenotypes. These will be based on obtaining the genome wide sequence of 4000 samples from the TwinsUK and ALSPAC cohorts (at 6x sequence coverage), and the exome sequence (protein coding regions and related conserved sequence) of 6000 samples selected for extreme phenotypes. Our studies will focus primarily on cardiovascular-related quantitative traits, obesity and related metabolic traits, neurodevelopmental disorders and a limited number of extreme clinical phenotypes that will provide proof-of-concept for future familial trait sequencing. We will analyse directly quantitative traits in the cohorts and the selected traits in the extreme samples, and also use imputation down to 0.1% allele frequency to extend the analyses to further sample sets with genome wide genotype data. In each case we will investigate indels and larger structural variants as well as SNPs, and use statistical methods that combine rare variants in a locus or pathway as well as single-variant approaches. The BioNED (Biomarkers for Childhood onset neuropsychiatric disorders) study has been carrying out detailed phenotypic assessments evaluating children with an autism spectrum disorder. These assessments included ADI-R, ADOS, neuropsychology, EEG etc. There are 56 DNA samples from this study (25 extracted from blood). For further information with regard to this cohort please contact Patrick Bolton (patrick.bolton@kcl.ac.uk).

Study EGAS00001000111

UK10K NEURO ASD GALLAGHER

In the UK10K project we propose a series of complementary genetic approaches to find new low frequency/rare variants contributing to disease phenotypes. These will be based on obtaining the genome wide sequence of 4000 samples from the TwinsUK and ALSPAC cohorts (at 6x sequence coverage), and the exome sequence (protein coding regions and related conserved sequence) of 6000 samples selected for extreme phenotypes. Our studies will focus primarily on cardiovascular-related quantitative traits, obesity and related metabolic traits, neurodevelopmental disorders and a limited number of extreme clinical phenotypes that will provide proof-of-concept for future familial trait sequencing. We will analyse directly quantitative traits in the cohorts and the selected traits in the extreme samples, and also use imputation down to 0.1% allele frequency to extend the analyses to further sample sets with genome wide genotype data. In each case we will investigate indels and larger structural variants as well as SNPs, and use statistical methods that combine rare variants in a locus or pathway as well as single-variant approaches. This is an Irish sample set of individuals with ASD (approximately 50% with comorbid intellectual disability). Individuals have been diagnosed with ADI/ ADOS, measures of cognition/ adaptive function. They represent a more severe, narrowly defined cohort of ASD subjects. Family histories are available for some with measures of broader phenotype. For further information on this cohort please contact Nadia Bolshakova (bolshakn@tcd.ie).

Study EGAS00001000112

Performance comparison of three DNA extraction kits on human whole-exome data from formalin-fixed paraffin-embedded normal and tumor samples

We generated 42 human whole-exome sequencing data sets from fresh-frozen (FF) and FFPE samples. These samples include normal and tumor tissues from two different organs (liver and colon), that we extracted with three different FFPE extraction kits (QIAamp DNA FFPE Tissue kit and GeneRead DNA FFPE kit from Qiagen, Maxwell\textsuperscript{TM} RSC DNA FFPE Kit from Promega). Variant calling analysis shows a very high rate of concordance between matched FF / FFPE pairs and equivalent performance for the three kits we analyzed. We find a significant variation in the difference of total number of variants called between FF and FFPE samples for the three different FFPE DNA extraction kits. Coverage analysis shows that FFPE samples have less good indicators than FF samples, yet the coverage quality remains above accepted thresholds. We detect limited but significant variations in coverage indicator values between the three FFPE extraction kits. Globally, the GeneRead and QIAamp kits have better variant calling and coverage indicators than the Maxwell kit on the samples used in this study, although this kit performs better on some indicators and has advantages in terms of practical usage. Taken together, our results confirm the potential of FFPE samples analysis for clinical genomic studies, but also indicate that the choice of a FFPE DNA extraction kit should be done with careful testing and analysis beforehand in order to maximize the accuracy of the results.

Study EGAS00001002631

23014 results for "ega"

in 130.13 milliseconds.