13304 results for ""pdc"+"heritage+language"+census+study"

in 59.29 milliseconds.

• 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

• Patient-derived neuroblastoma model system OHC-NB1

  Accurate disease modeling presents a bottleneck against effective translation of novel neuroblastoma therapies into the clinic. Recently, considerable intra-tumor genetic heterogeneity has been observed in neuroblastoma, which calls for preclinical models that reflect this feature. Flexibility in a testing platform is also desirable to support assessment of different endpoints and tumor cell characteristics. We present the novel neuroblastoma model OHC-NB1, propagated from a bone marrow metastasis from a patient with INRG stage M, MYCN-amplified neuroblastoma at first diagnosis. To model different aspects of the disease, several preclinical models were developed from the bone marrow aspirate: a tissue plastic-adherent cell line, a 3-dimensional (3D) spheroid culture, and cell line-derived xenografts that were subcutaneously passaged in mice. Typical aspects of MYCN-amplified neuroblastoma, such as MYCN amplification in double minutes, 1p deletion, 17q gain, and a diploid karyotype, are present in OHC-NB1 and persist in all culture types. Whole exome sequencing revealed intra-tumor heterogeneity in the bone marrow metastasis reflected in the culture models. The set of OHC-NB1 cultures with a genetically heterogeneous MYCN-amplified neuroblastoma background provide a flexible preclinical testing system that recapitulates important features of metastatic neuroblastoma.

  Study EGAS00001003031

• 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

• A Single-cell Atlas of Progressive TKI Resistance in Chronic Myeloid Leukemia

  High-dimensional scRNA seq and CyTOF were used to identify features predictive of treatment outcome at diagnosis. An optimal TKI response was characterised by erythroid lineage skewing of the CD34+ HSPC compartment. While TKI-resistant CP was reminiscent of a BC-like state, where HSPCs were enriched for inflammation, stemness and quiescence. We designed a machine-learning approach which assessed the prognostic potential of all 32 sub-populations of our scRNA seq dataset. LSCs and NK populations harboured the highest prognostic power at diagnosis. Within the LSCs, erythroid lineage priming and MYC activation were features of an optimal and poor TKI responses, respectively. With the NK cell compartment, apart from a higher abundance of NK cells in TKI optimal responders, a memory-like HLA-DR+ adaptive NK and KLRC1+ NK populations were hallmarks of an optimal and poor TKI responses, respectively. Mechanistically, both cell intrinsic and extrinsic mechanisms contributed towards the higher sensitivity of EPs to TKI therapy. In summary, our high-dimensional single-cell atlas of TKI resistance highlights pivotal transcriptional features associated with TKI-resistance disease, some of which have the potential to be exploited as biomarkers.

  Study EGAS00001005509

• Clonal relationships of ductal carcinoma in situ and recurrent invasive breast cancers defined by genomic analysis

  Ductal Carcinoma In Situ (DCIS) is the most common form of pre-invasive breast cancer and despite treatment a small fraction (5-10%) of DCIS patients present with invasive disease many years later. A fundamental biologic question is whether the invasive disease recurring in the same breast is established by tumor cells in the initial DCIS or represents new unrelated disease. To address this question, we performed genomic analyses on the initial pure DCIS lesion and paired invasive recurrent tumors in 95 patients together with single cell DNA sequencing in a subset of cases. Our data shows that in 75% the invasive recurrence was clonally related to the initial DCIS, suggesting that the tumor cells were not eliminated during the initial treatment with surgery +/- radiotherapy. Surprisingly however, 18% were clonally unrelated to the DCIS, representing new independent lineages, and 7% of cases were ambiguous. Our findings show that although DCIS is often the precursor of invasive recurrence, a significant fraction of invasive recurrences are unrelated to the initial DCIS. This knowledge is essential for accurate risk evaluation of DCIS treatment de-escalation strategies and the identification of predictive biomarkers.

  Study EGAS00001005784

• De novo assembly of 150 Danish genomes reveals rich structural complexity

  Most known genetic variation in human genomes has been called from comparison of short reads to the reference genome, an approach biased against finding complex variation. We sequenced 150 individuals from 50 parent-offspring trios with multiple insert-size libraries to very high coverage. We show that each genome could be independently de novo assembled into a small number of high-quality scaffolds (median N50 > 21 Mb), each of quality comparable to long read assemblies while being very cost-effective. We show that our variant call set from comparing de novo assemblies is far more complete in terms of complex variation than previous studies. Importantly, even the complex 4-5 Mb extended MHC region was assembled and resolved into haplotypes, revealing >700kb novel sequence in this important region of the genome, and major parts of the Y chromosome including some palindromes were assembled with high accuracy. Finally, we show that our variant call-set allows for the genotyping of many more complex variants when used as a reference-panel for imputation into SNP-chip data or into previously resequenced genomes.

  Study EGAS00001002108

• Clonal Decomposition and DNA Replication States Defined by Scaled Single-Cell Genome Sequencing

  Accurate measurement of clonal genotypes, mutational processes, and replication states from individual tumor-cell genomes will facilitate improved understanding of tumor evolution. We have developed DLP+, a scalable single-cell whole-genome sequencing platform implemented using commodity instruments, image-based object recognition, and open source computational methods. Using DLP+, we have generated a resource of 51,926 single-cell genomes and matched cell images from diverse cell types including cell lines, xenografts, and diagnostic samples with limited material. From this resource we have defined variation in mitotic mis-segregation rates across tissue types and genotypes. Analysis of matched genomic and image measurements revealed correlations between cellular morphology and genome ploidy states. Aggregation of cells sharing copy number profiles allowed for calculation of single-nucleotide resolution clonal genotypes and inference of clonal phylogenies and avoided the limitations of bulk deconvolution. Finally, joint analysis over the above features defined clone-specific chromosomal aneuploidy in polyclonal populations.

  Study EGAS00001003190

• B cell receptor repertoire kinetics after SARS-CoV-2 infection and vaccination

  B cells play a central role in the immune response to both SARS-CoV-2 infection and vaccination, but the development of the B cell receptor (BCR) repertoire in both contexts has not been defined nor compared. We analysed serial samples from 171 SARS-CoV-2-infected individuals with a range of disease severities together with 63 vaccine recipients, and found marked differences in the global BCR repertoire after natural infection compared to vaccination. Following infection, the proportion of BCRs bearing IgG1/3 and IgA1 isotypes increased, somatic hypermutation (SHM) was markedly decreased and, in patients with severe disease, expansion of IgM and IgA clones was observed. In contrast, after vaccination the proportion of BCRs bearing IgD/M isotypes increased, SHM was unchanged and expansion of IgG clones was prominent. Infection generated a broad distribution of SARS-CoV-2-specific clones predicted to target the spike protein whilst vaccination produced a more focused response mainly targeting the spike’s receptor-binding domain. These findings offer insights into how different immune exposure to SARS-CoV-2 impacts upon BCR repertoire development, potentially informing vaccine strategies.

  Study EGAS00001005816

• Target gene sequencing for human normal endometrial glands

  To elucidate the timing and mechanism of the clonal expansion of somatic mutations in cancer-associated genes in the normal endometrium, we conducted target sequencing of 112 genes for 1,298 endometrial glands and matched blood samples from 36 women. By collecting endometrial glands from different parts of the endometrium, we showed that multiple glands with the same somatic mutations occupied substantial areas of the endometrium. The 112 genes are as follows: ABCC1, ACRC, ANK3, ARHGAP35, ARID1A, ARID5B, ATCAY, ATM, ATR, BARD1, BCOR, BRCA1, BRCA2, BRD4, BRIP1, CAMTA1, CDC23, CDYL, CFAP54, CHD4, CHEK1, CHEK2, CTCF, CTNNB1, CUX1, DGKA, DISP2, DYNC2H1, EMSY, FAAP24, FAM135B, FAM175A, FAM65C, FANCA, FANCB, FANCC, FANCD2, FANCE, FANCF, FANCG, FANCI, FANCL, FANCM, FAT1, FAT3, FBN2, FBXW7, FGFR2, FRG1, GPR50, HEATR1, HIST1H4B, HNRNPCL1, HOOK3, KIAA1109, KIF26A, KMT2B, KMT2C, KRAS, LAMA2, LRP1B, MLH1, MON2, MRE11A, MSH2, MSH6, MTOR, NBN, PALB2, PHEX, PIK3CA, PIK3R1, PLXNB2, PLXND1, PMS2, POLE, POLR3B, PPP2R1A, PTEN, PTPN13, RAD50, RAD51, RAD51B, RAD51C, RAD51D, RAD52, RAD54B, RAD54L, RICTOR, SACS, SIGLEC9, SLC19A1, SLX4, SPEG, STT3A, TAF1, TAF2, TAS2R31, TFAP2C, TNC, TONSL, TP53, TTC6, UBA7, VNN1, WT1, XIRP2, ZBED6, ZC3H13, ZFHX3, ZFHX4, ZMYM4.

  Study EGAS00001005914

• 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

• Changes in alternative splicing and associated neo-antigens due to therapy

  Recent studies have identified alternative splicing (AS) as a novel source of neoantigens for immunotherapy. Surprisingly little is known about the AS milieu in recurrent glioblastoma (GBM), despite this being the venue for most clinical trials. We profiled 29 primary-recurrent paired human GBM specimens via RNA sequencing and re-analyzed RNA-sequencing data from non-malignant human brain tissues. From these data, we reconstructed the landscape of AS in GBM through recurrence and contrasted that to isoforms observed in non-malignant brain. The AS events we identified were cross-referenced with single-cell GBM atlases to determine cell-type specific splicing patterns. From this we identified novel splicing events in cell-surface proteins that are suitable targets for engineered T-cell therapies. We found recurrent-specific isoforms of mitogen-activated kinase pathway genes which are expressed exclusively by GBM stem-like cells that enhance invasiveness. These studies shed light on the effect of therapy on AS and identify novel targets for emerging immunotherapies.

  Study EGAS00001004524

• Patient-tailored design for selective co-inhibition of leukemic cell subpopulations

  The extensive primary and secondary drug resistance in many cancer types requires rational approaches to design personalized and selective combinatorial therapies that do not only show synergistic effect in overall cancer cell killing but also result in minimal toxic side effects on non-malignant cells. To address the combinatorial explosion in the number of relevant combinations, we implemented a machine learning approach that prioritizes patient-customized drug combinations with a desired synergy-efficacy-toxicity balance by combining single-cell RNA-sequencing with ex vivo single-agent testing in scarce patient-derived primary cells. When applied to two diagnostic and two refractory AML patient cases, each with a different genetic background, our integrated approach accurately predicted patient-specific combinations that were shown to result not only in synergistic cancer cell co-inhibition but were also capable of targeting specific AML cell subpopulations that emerge in differing stages of disease pathogenesis or treatment regimens. Our data-driven approach provides an unbiased means for systematic identification of personalized combinatorial regimens that selectively co-inhibit leukemic cells while avoiding inhibition of non-malignant cells, and highlight the relevance of considering cell heterogeneity for personalized cancer therapy.

  Study EGAS00001004614

• Allele-specific expression of GATA2 due to epigenetic dysregulation in double mutated CEBPA AML

  Transcriptional deregulation is a central event in the development of acute myeloid leukemia (AML). To identify potential disturbances in gene regulation, we conducted an unbiased screen of allele-specific expression (ASE) in 209 AML cases. The gene encoding GATA binding protein 2 (GATA2) displayed ASE more often than any other myeloid or cancer-related gene. GATA2 ASE was strongly associated with CEBPA double mutations (CEBPA DM), with 95% of cases presenting GATA2 ASE. In CEBPA DM AML with GATA2 mutations, the mutated allele was preferentially expressed. We found that GATA2 ASE is a somatic event lost in complete remission, supporting the notion that it plays a role in CEBPA DM AML. Acquisition of GATA2 ASE involved silencing of one allele via promoter methylation, compensated by overactivation of the other allele, thereby preserving expression levels. Notably, promoter methylation was also lost in remission together with GATA2 ASE. In summary, we propose that GATA2 ASE is acquired by epigenetic mechanisms and is a prerequisite for the development of AML with CEBPA DM. This finding constitutes a novel example of an epigenetic hit cooperating with a genetic hit in the pathogenesis of AML.

  Study EGAS00001004684

• Genomics-based personalized oncology in cancer of unknown primary (CUP, project H021)

  To perform a comprehensive genomic characterization of 70 patients suffering from cancer of unknown primary (CUP) we used whole-exome, whole-genome, transcriptome and methylome analysis. We detected a substantial mutational heterogeneity with genes most commonly affected by SNVs, indels and fusions being TP53, TTN, MUC16, ABCA13, COL6A3, KRAS, LRP1B, XIRP2 and CSMD3. The most common fusion involved FGFR2, the most common focal deletion affected CDKN2A. A molecular tumor board recommended genomics-based therapies in 56/70 (80%) patients which were applied in 20/56 (35.7%) cases. Entity predictions based on transcriptome and methylome data could be made in up to 62/70 (88.6%) cases but were conclusive in only 16/48 (33.3%) cases. Germline analysis revealed 6 (likely) pathogenic mutations in 5 patients. Recommended therapies translated into a mean PFS2/1 ratio of 3.61 (median=2.25) with a median PFS1 of 89 days (n=17) compared to a median PFS2 of 182.5 days (n=20). Our data emphasize the clinical benefit of comprehensive genetic approaches in diagnostic and therapeutic management and underline the need for innovative, mechanism-based clinical trials in this heterogeneous group of diseases.

  Study EGAS00001004786

• Accurate detection of tumor-specific fusion genes reveals strongly immunogenic personal neo-antigens

  Fusion genes arising from cancer-associated somatic mutations are a potential rich source for highly immunogenic neo-antigens. However, their exploitation as targets for personalized cancer immunotherapy is currently limited by the lack of computational tools allowing transcriptome-wide identification of unique fusion genes in an accurate and sensitive manner. Here, we present EasyFuse, a computational pipeline, to detect individual and cancer-specific fusion genes in next-generation-sequencing transcriptome data obtained from human cancer samples. Using machine learning, EasyFuse predicts personal fusion genes with high precision and sensitivity and outperforms previously described approaches as qualified by an unprecedented ground-truth dataset of >1500 verification experiments in relevant patient samples. By testing immunogenicity with autologous blood lymphocytes from cancer patients we detected pre-established CD4+ and CD8+ T cell responses for 10 of 21 (48%), and for 1 of 30 (3%) of identified fusion genes, respectively. In conclusion, we demonstrate accurate detection of cancer-specific fusion genes. The high frequency of T cell responses detected in cancer patients support the relevance of private fusion genes as neo-antigens for personalized immunotherapies, especially for tumors with low point mutation burdens.

  Study EGAS00001004877

• 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

• Clonal somatic copy number altered driver events inform drug sensitivity in high-grade serous ovarian cancer

  Chromosomal instability is a major challenge to patient stratification and targeted drug development for high-grade serous ovarian carcinoma (HGSOC). Here we show that somatic copy number alterations (SCNAs) in frequently amplified HGSOC cancer genes significantly correlate with gene expression and methylation status. We identified five prevalent clonal driver SCNAs (chromosomal amplifications encompassing MYC, PIK3CA, CCNE1, KRAS and TERT) from multi-regional HGSOC data and reasoned that their strong selection should prioritise them as key biomarkers for targeted therapies. We used primary HGSOC spheroid models to test interactions between in vitro targeted therapy and SCNAs. MYC chromosomal copy number was associated with in-vitro and clinical response to paclitaxel and in-vitro response to mTORC1/2 inhibition. Activation of the mTOR survival pathway in the context of MYC-amplified HGSOC was statistically associated with increased prevalence of SCNAs in genes from the PI3K pathway. Co-occurrence of amplifications in MYC and genes from the PI3K pathway was independently observed in squamous lung cancer and triple negative breast cancer. These results suggest that identifying co-occurrence of clonal driver SCNA genes could be used to tailor therapeutics for precision medicine.

  Study EGAS00001006200

• MutWP4__CRUK_Grand_Challenge_Mutographs_of_Cancer__Pancreatic_Organoids

  The Mutographs project aims to advance our understanding of the causes of cancer through studies of mutational signatures. Led by Mike Stratton, together with Paul Brennan, Ludmil Alexandrov, Allan Balmain, David Phillips and Peter Campbell, this large-scale international research endeavour was awarded a Cancer Research UK Grand Challenge. Different patterns of somatic mutation are generated by the different environmental, lifestyle and genetic factors that cause cancer, many of them are still unknown. Within Mutographs, Kings College London will characterise the mutational signatures induced by putative human carcinogens in order to identify the origins of mutational signatures found in human cancers. To achieve this human organoid cell cultures will be exposed to a representative catalogue of known or suspected human carcinogens and mutagens and, using whole genome sequencing, the patterns of mutations induced by them will be determined. Somatic mutational signatures will be subsequently extracted by non-negative matrix factorisation methods and correlated with exposure data. Through an enhanced understanding of cancer aetiology, Mutographs unprecedented effort is anticipated to outline modifiable risk factors, lead to new approaches to prevent cancer, and provide opportunities to empower early detection, refine high-risk groups and contribute to further therapeutic development.

  Study EGAS00001003520

• MGP Panel: a comprehensive targeted genomics panel for molecular profiling of multiple myeloma patients

  We designed a comprehensive multiple myeloma (MM) targeted sequencing panel to identify common genomic abnormalities in a single assay and validated it against known standards. The panel comprised 228 genes/exons for mutations, 6 regions for translocations, and 56 regions for copy number abnormalities (CNAs). Toward panel validation, targeted sequencing was conducted on 233 patient samples and further validated using clinical fluorescence in situ hybridization (FISH) (translocations), multiplex ligation probe analysis (MLPA) (CNAs), whole genome sequencing (WGS) (CNAs, mutations, translocations) or droplet digital PCR (ddPCR) of known standards (mutations). Canonical IgH translocations were detected in 43.2% of patients by sequencing, and aligned with FISH except for one patient. CNAs determined by sequencing and MLPA for 22 regions were comparable in 103 samples and concordance between platforms was R2=0.969. VAFs for 74 mutations were compared between sequencing and ddPCR with concordance of R2=0.9849. In summary, we have developed a targeted sequencing panel that is as robust or superior to FISH and WGS. This molecular panel is cost effective, comprehensive, clinically actionable and can be routinely deployed to assist risk stratification at diagnosis or post-treatment to guide sequencing of therapies.

  Study EGAS00001006222

• Discovery of a highly widespread bacteriophage family and its associations to metabolic syndrome gut microbiomes

  There is significant interest in altering the course of cardiometabolic disease development via the gut microbiome. Nevertheless, the highly abundant phage members of the complex gut ecosystem -which impact gut bacteria- remain understudied. Here, we characterized gut virome changes associated with metabolic syndrome (MetS), a highly prevalent clinical condition preceding cardiometabolic disease, in 196 participants with a combination of whole genome shotgun and virus like particle sequencing. MetS gut virome populations exhibited decreased richness and diversity, but larger inter-individual variation. These populations were enriched in phages infecting Bacteroidaceae and depleted in those infecting Ruminococcaeae. Differential abundance analysis identified eighteen viral clusters (VCs) as significantly associated with either MetS or healthy viromes. Among these are a MetS-associated Roseburia VC that is related to healthy control-associated Faecalibacterium and Oscillibacter VCs. Further analysis of these VCs revealed the Candidatus Heliusviridae, a highly widespread gut phage lineage found in 90+% of the participants. The identification of the temperate Ca. Heliusviridae provides a novel starting point to a better understanding of the effect that phages have on their bacterial hosts and the role that this plays in MetS

  Study EGAS00001006260

• Clonal relationships of ductal carcinoma in situ and recurrent invasive breast cancers defined by genomic analysis

  Ductal Carcinoma In Situ (DCIS) is the most common form of pre-invasive breast cancer and despite treatment a small fraction (5-10%) of DCIS patients present with invasive disease many years later. A fundamental biologic question is whether the invasive disease recurring in the same breast is established by tumor cells in the initial DCIS or represents new unrelated disease. To address this question, we performed genomic analyses on the initial pure DCIS lesion and paired invasive recurrent tumors in 95 patients together with single cell DNA sequencing in a subset of cases. Our data shows that in 75% the invasive recurrence was clonally related to the initial DCIS, suggesting that the tumor cells were not eliminated during the initial treatment with surgery +/- radiotherapy. Surprisingly however, 18% were clonally unrelated to the DCIS, representing new independent lineages, and 7% of cases were ambiguous. Our findings show that although DCIS is often the precursor of invasive recurrence, a significant fraction of invasive recurrences are unrelated to the initial DCIS. This knowledge is essential for accurate risk evaluation of DCIS treatment de-escalation strategies and the identification of predictive biomarkers.

  Study EGAS00001006306

• Pharmacogenomic analysis of patient-derived tumor cells in gynecologic cancers

  Gynecologic malignancy is one of the leading causes of mortality in female adults worldwide. Comprehensive genomic analysis has revealed a list of molecular aberrations that are essential to tumorigenesis, progression, and metastasis of gynecologic tumors. However, targeting such alterations has frequently led to treatment failures due to underlying genomic complexity and activation of various tumor cell survival pathway molecules. A compilation of molecular characterization of tumors with pharmacological drug response is the next step towards clinical application of patient-tailored treatment regimens. Toward this goal, we have established a library of 139 gynecologic tumors including cervical, endometrial, and epithelial ovarian cancers (EOCs) and uterine sarcomas that were genomically and/or pharmacologically annotated and explored dynamic pharmacogenomic associations against 37 molecularly targeted drugs. We discovered lineage-specific drug sensitivities based on subcategorization of gynecologic tumors and identified TP53 mutation as a molecular determinant that elicit therapeutic response to poly (ADP-Ribose) polymerase (PARP) inhibitor. We further identified transcriptome expression of inhibitor of DNA biding 2 (ID2) as a potential predictive biomarker for treatment response to olaparib. Together, our results demonstrate potential utility of rapid drug screening combined with genomic profiling for precision treatment of gynecologic cancers.

  Study EGAS00001003556

• Treatment-mediated selection of lethal prostate cancer clones defined by copy number architectures

  Despite initial responses to hormone treatment, metastatic prostate cancer invariably evolves to a lethal state. To characterize the intra-patient relationships of metastases that evade treatment, we performed genome- wide copy number profiling and bespoke approaches targeting the androgen receptor (AR) on 142 metastatic regions from 10 organs harvested post-mortem from nine men who died from prostate cancer. We identified diverse and patient-unique alterations clustering around the AR in metastases from every patient with evidence of independent acquisition of related genomic changes within an individual and, in some patients, the co-existence of AR-neutral clones. Using the genomic boundaries of pan-autosome copy number change, we confirmed a common clone of origin across metastases and diagnostic biopsies; and identified in individual patients, clusters of metastases occupied by dominant clones with diverged autosomal copy number alterations. Autosome-defined clusters were characterized by cluster-specific AR gene architectures that in two index cases were topologically more congruent than by chance (p-values 0.03, 3.07x10-8). Integration with anatomical site suggested patterns of spread and points of genomic divergence. Copy number boundaries identified treatment-selected clones with putatively distinct lethal trajectories.

  Study EGAS00001006598

• DNA methylation landscape of prostate cancer

  Metastases from primary prostate cancers to rare locations, such as the brain, are becoming more common due to longer life expectancy thanks to improved treatments. Epigenetic dysregulation is a feature of primary prostate cancer, and distinct DNA methylation profiles are associated with the mutually exclusive SPOP mutant or TMPRSS2-ERG fusion genetic backgrounds. Using a cohort of 42 patients with prostate cancer brain metastases (PCBM), with matched primary samples available for 17 patients, we carried out DNA methylation analysis using Illumina EPIC array. We examined the epigenetic distinction between primary PCa and PCBM, the association between epigenetic alterations and mutational background, and particular epigenetic alterations that may be associated with PCBM in particular. Multiregion sampling of PCBM revealed epigenetic heterogeneity in metastases. We found aberrant methylation in PCBM is associated with mutational background, PRC2 complex activity and upregulation of EZH2, an effect that is particularly pronounced in SPOP mutant PCBM. While we detected a hypermethylator phenotype in PCBM, we also observed hypomethylation occurring at the promoters of genes involved in neuroactive ligand-receptor interactions or the calcium signalling pathway, such as GABRB3 and GRM5, a change not observed in non-brain prostate cancer metastases.

  Study EGAS00001006670

• Pre-neoplastic liver colonization by 11p15.5 altered mosaic cells in young children with hepatoblastoma

  In children, liver tumors are a heterogeneous group representing 1-2% of all cancers with the most common diagnosis being hepatoblastoma (HB). While HBs are predominantly sporadic, around 15% of cases develop as part of predisposition syndromes such as familial adenomatous polyposis (APC mutation) or Beckwith-Wiedemann (BWS; 11p15.5 locus altered). Here, we identified a mosaic genetic alteration of 11p15.5 locus in the liver in 17% of HB patients without a clinical diagnosis of BWS syndrome. These mosaic 11p15.5 alterations were never found in children with other types of liver tumor (hepatocellular carcinoma or adenoma, fibrolamellar carcinoma). Bulk transcriptomic analysis revealed the abundance of extracellular matrix, angiogenesis, and progenitor gene signatures in mosaic livers. In situ hybridization assay RNAscope allowed for the detection of 11p15.5 altered hepatocytes mosaicism in liver FFPE sections of HB patients presenting IGF2 overexpression and H19 downregulation together with an alteration of the zonal architecture of the liver. Spatial transcriptomics analysis revealed an overexpression of the IGF2 pathway, lipogenesis, and platelet activation in 11p15.5 altered hepatocytes. These data provide insights for 11p15.5 mosaicism detection as well as functional consequences of these alterations during the early steps of carcinogenesis.

  Study EGAS00001006692

• Targeted therapy of advanced parathyroid carcinoma guided by genomic and transcriptomic profiling (hipo_021)

  Parathyroid carcinoma (PC) is an orphan malignancy with a high risk of recurrence after surgery. Tumor-directed systemic treatments for PC are not established. We used whole-genome and RNA sequencing in four patients with advanced PC to identify molecular alterations that could guide clinical management. In two cases, the genomic and transcriptomic profiles provided targets for experimental therapies that resulted in biochemical response and prolonged disease stabilization: (i) immune checkpoint inhibition with pembrolizumab based on high tumor mutational burden and a single-base substitution signature associated with APOBEC (apolipoprotein B mRNA editing enzyme, catalytic polypeptide-like) overactivation; (ii) multi-receptor tyrosine kinase inhibition with lenvatinib due to overexpression of FGFR1 and RET and, later in the course of disease, PARP inhibition with olaparib prompted by signs of defective homologous recombination DNA repair. In addition, our data provided new insights into the molecular landscape of PC with respect to the genome-wide footprints of specific mutational processes and pathogenic germline alterations. These data underscore the potential of comprehensive molecular analyses to improve care for patients with ultra-rare cancers based on insight into disease biology.

  Study EGAS00001006747

• Case Report: Rare IKZF1 gene fusions identified in neonate with congenital KMT2A-rearranged Acute Lymphoblastic Leukemia

  Here we describe a rare case of congenital KMT2Ar ALL presenting with co-occurring IKZF1 gene fusions and a predictably aggressive disease trajectory. We report for the first time, the novel IKZF1::TUT1 and KDM2A::IKZF1 gene fusions. Rearrangements in-volving KMT2A are commonly retained in relapsed infant ALL, however, in this case the KMT2A::AFF1 gene fusion did not appear to be the lesion driving leukemic relapse. Instead, our data suggest that relapse was driven by IKZF1::TUT1. This gene fusion remained in all samples investigated, including the on-blinatumomab therapy sample taken immediately prior to relapse. Conversely, the KMT2A::AFF1 gene fusion was only detected in the diagnosis and refractory post-induction samples highlighting a key role for IKZF1::TUT1 in disease pathogenesis. Intriguingly, both IKZF1 gene fusions are predicted to be out-of-frame, however, our data demonstrate the IKZF1 gene is still expressed. This is not unprecedented and it has previously been observed that out-of-frame fusions can cause transcriptional activation/repression of genes involved in the fusions leading to increases or decreases of their expression and the as-sociated functional outcomes.

  Study EGAS00001006947

• HLA-DR is absent in primitive macrophages through epigenetic silencing of master regulator CIITA

  The earliest macrophages are generated during embryonic development from erythro-myeloid progenitors (EMPs) via primitive haematopoiesis. This process is poorly understood in humans and in the mouse is thought to be spatially restricted to the yolk sac. Human fetal placental macrophages, Hofbauer cells (HBC), arise during the primitive haematopoietic wave at ~18 days post conception. Here we identify a population of placental erythro-myeloid progenitors (PEMPs) in the early human placenta that give rise to HBC. PEMP are fetal CD34+CD43+ progenitors found exclusively at early gestational timepoints. Transcriptomic analyses reveal that PEMP have conserved features of primitive yolk sac EMPs, including the lack of HLF expression. Using in vitro single-cell culture experiments we show that PEMP generate HBC-like cells that completely lack HLA-DR expression  and demonstrate that this is a conserved feature of all macrophages generated by primitive haematopoiesis in humans. Finally, we demonstrate that the class II transactivator, CTIIA, the master regulator HLA-DR expression, is epigenetically silenced during primitive haematopoiesis, leading to HLA-DRneg macrophages. These findings demonstrate that HBC are derived locally from PEMP and establish the human placenta as an additional site of primitive haematopoiesis.

  Study EGAS00001006981

• RNAseq

  Familial Mediterranean Fever (FMF) is a prototypical periodic fever syndrome caused by genetic variation in MEFV. While it is known that especially IL-1β responses are dysregulated in FMF, its innate immune landscape has not been comprehensively described. Therefore, we extensively characterized the function of monocytes and neutrophils in patients with FMF in between disease attacks, as well as their circulating inflammatory proteome. We found that monocyte IL-1 and IL-6 production was enhanced following a range of stimulations, in concordance with alterations in the plasma inflammatory proteome. In contrast, neutrophil function in FMF was normal. Additionally, ATAC-seq and RNA-seq analyses revealed important epigenetic and transcriptional reprogramming arguing for monocyte function dysregulation. Interestingly, chromatin-accessibility was down-regulated in genomic regions related to cellular responses. This argues that the primary immune dysregulation in monocytes due to MEFV mutations leads to a chronic inflammatory profile that is subsequently associated with counterregulatory epigenetic and transcriptional changes reminiscent of tolerance (rather than trained immunity) induction. These data increase our understanding of the innate immune changes in FMF, that can be used in the design of novel approaches to manage chronic inflammation of these patients.

  Study EGAS00001007165

• Single-cell profiling of co-cultures of GSCCs and macrophages

  Glioblastoma (GBM) remains the most malignant primary brain tumor, with a median survival rarely exceeding 2 years. Tumor heterogeneity and an immunosuppressive microenvironment are key factors contributing to the poor response rates of current therapeutic approaches. GBM-associated macrophages (GAMs) often exhibit immunosuppressive features that promote tumor progression. However, their dynamic interactions with GBM tumor cells remain poorly understood. Here, we used patient-derived GBM stem cell cultures and combined single-cell RNA sequencing of GAM-GBM co-cultures and real-time in vivo monitoring of GAM-GBM interactions in orthotopic zebrafish xenograft models to provide insight into the cellular, molecular, and spatial heterogeneity. Our analyses revealed substantial heterogeneity across GBM patients in GBM-induced GAM polarization and the ability to attract and activate GAMs – features that correlated with patient survival. Differential gene expression analysis, immunohistochemistry on original tumor samples, and knock-out experiments in zebrafish subsequently identified LGALS1 as a primary regulator of immunosuppression. Overall, our work highlights that GAM-GBM interactions can be studied in a clinically relevant way using co-cultures and avatar models, while offering new opportunities to identify promising immune-modulating targets.

  Study EGAS00001007482

• Massively parallel functional dissection of schizophrenia associated non-coding genetic variants

  Schizophrenia (SCZ) is a severe mental disorder affecting 1% of the world population. SCZ is characterized by an underlying genetic architecture that is highly polygenic. Genome wide association studies have identified thousands of genetic variants that are statistically linked to the disease. However, the translation of these associations into insights on the pathomechanisms has been challenging because the causal genetic variants, their molecular function, and their target genes remain largely unknown. To address these questions, we combined induced pluripotent stem cell technology with a massively parallel variant annotation pipeline (MVAP) to functionally characterize 35,000 SCZ associated non-coding genetic variants. This approach identified a set of 620 (1.7%) single nucleotide polymorphisms as functional on the molecular level in a highly cell type and condition specific fashion. Subsequent multi-modal integration of epigenomic data combined with CRISPR screening in human neurons enabled us to systematically translate SCZ variant associations into target genes, biological processes, and ultimately alterations of neuronal physiology. These results provide a new high-resolution map of functional variant-gene combinations and offer comprehensive biological insights into the developmental context and stimulus dependent molecular processes modulated by SCZ genetic variation beyond statistical association.

  Study EGAS00001007542

• Immune and clinicopathological features predict HER2-positive breast cancer prognosis in the neoadjuvant NeoALTTO and CALGB 40601 trials

  The identification of prognostic drivers in patients receiving neoadjuvant therapy is crucial for treatment optimization in HER2-positive breast cancer, with the immune microenvironment being a key factor. Here, we investigated the complexity of B and T cell receptor (BCR and TCR) repertoires, reflecting antigen-specific immune response, in the context of two phase III trials, namely NeoALTTO and CALGB 40601, evaluating neoadjuvant paclitaxel with trastuzumab and/or lapatinib in women affected by HER2-positive breast cancer. BCR/TCR features were heterogeneous, and BCR measures describing clonal expansion demonstrated independent prognostic value. A model developed in NeoALTTO and validated in CALGB 40601 predicted event-free survival (EFS) by integrating hormone receptor and clinical nodal status, treatment response (breast pCR)breast pathological complete response (pCR) as well as tumor-infiltrating lymphocytes and BCR evenness, measuring the equality of BCR clonal distribution. The HER2-EveNT score derived from the model allows the identification of patients with 5-year EFS > 90%, and, in those not achieving pCR, of a subgroup of immune-enriched tumors with an excellent outcome despite residual disease.

  Study EGAS00001007563

• Discovery of immunotherapy targets for pediatric solid and brain tumors by exon-level expression

  Immunotherapy with chimeric antigen receptor T cells for pediatric solid and brain tumors is constrained by available targetable antigens. Cancer-specific exons present a promising reservoir of targets; however, these have not been explored and validated systematically in a pan-cancer fashion. To identify cancer specific exon targets, we analyze 1,532 RNA-seq datasets from 16 types of pediatric solid and brain tumors for comparison with normal tissues using a newly developed workflow. We find 2,933 exons in 157 genes encoding proteins of the surfaceome or matrisome with high cancer specificity either at the gene (n=148) or the alternatively spliced isoform (n=9) level. Expression of selected alternatively spliced targets, including the EDB domain of fibronectin 1, and gene targets, such as COL11A1, are validated in pediatric patient derived xenograft tumors. We generate T cells expressing chimeric antigen receptors specific for the EDB domain or COL11A1 and demonstrate that these have antitumor activity. The full target list, explorable via an interactive web portal (https://cseminer.stjude.org/), provides a rich resource for developing immunotherapy of pediatric solid and brain tumors using gene or AS targets with high expression specificity in cancer.

  Study EGAS00001007766

• M116 Proteome Extracellular Vesicle Profiling

  For extracellular vesicles (ECV) extraction from plasma, 1M ammonium acetate was added to precipitate ECVs on ice for 45 min. Then, 100 mM ammonium acetate was added to the mixture, and ECVs were precipitated by centrifugation at 20,000g for 30 min. ECVs were washed with 50 mM ammonium bicarbonate (Sigma-Aldrich, Cat. 1066-33-7). Then, 600 µl of 1% ammonium deoxycholate (Sigma Aldrich, Cat. K2755-1MG) were added. The concentration of protein in each sample was measured using a bicinchoninic acid assay (BCA assay). Ammonium bicarbonate was used to dilute 500 µg of protein into a final volume of 500 µl. Next, dithiothreitol (Sigma-Aldrich, Cat. D9779) was added to obtain a final concentration of 20 mM, followed by iodoacetamide (Sigma-Aldrich, Cat. I6125) to a final concentration of 40 mM. Next, trypsin (Roche, Cat. RTRYP-RO) was added to the sample in a 1:25 protein ratio and incubated at 37°C overnight. Next day, formic acid (ThermoFisher Scientific, Cat. 28905) in a final concentration of 0.1% was added and extraction of proteins was done using Empore™ Solid Phase Extraction Cartridges (3M), following manufacturer’s instructions. The eluted samples were then centrifuged for 90 min. using a speed vacuum centrifuge (Thermo, RC1010), followed by snap freezing in liquid nitrogen. Then, the samples were kept in a freeze dryer (LyoDry Compact Benchtop, MechaTech) overnight. Next, the samples were reconstituted in 30 µl of 0.1% formic acid (FA) and an o-Phthaladehyde (Oparil) assay was performed to determine the concentration of each sample. After that, the sample was prepared in a concentration of 0.5 µg/µl using 0.1% FA and alcohol dehydrogenase (ADH). The samples were prepared in glass mass spectrometry vials for proteomic analysis using a Waters Synapt G2Si High-Definition Mass Spectrometry (Waters Corporation) operated by the MassLynx 4.1., 110 min. running time with 2 µl of an injection containing 1 µg of peptide. Quality controls were also run along with samples to guarantee consistency. Pooled quality controls were made from all samples, in which the samples were run at the beginning, middle and end of the mass spectrometry run. Samples were randomized before running the experiment. The proteomic data was then imported into Progenesis software 4.2 (Nonlinear Dynamic, UK) to identify and quantify peptides and proteins.

  Dataset EGAD50000001679

• Genome Diversity in Africa Project: Benin (2021-02-16)

  Globally, human populations show structured genetic diversity as a result of geographical dispersion, selection and drift. Understanding this genetic variation can provide insights into the evolutionary processes that shape both human adaptation and variation in disease. Populations from SSA have the highest levels of genetic diversity. This characteristic, in addition to historical genetic admixture, can lead to complexities in the design of studies assessing the genetic determinants of disease and human variation. However, such studies of African populations are also likely to provide new opportunities to discover novel disease susceptibility loci and variants and refine gene-disease association signals. A systematic assessment of genetic diversity within SSA would facilitate genomic epidemiological studies in the region. The Genome Diversity in Africa Project (GDAP) aims to produce a comprehensive catalogue of human genetic variation in SSA, including single nucleotide polymorphisms (SNPs), structural variants, and haplotypes. This resource will make a substantial contribution to understanding patterns of genetic diversity within and among populations in SSA, as well as providing a global resource to help design, implement and interpret genomic studies in SSA populations and studies comprising globally diverse populations, complementing existing genomic resources. Specifically, we plan to carry out high depth whole genome sequencing of up to 2000 individuals across Africa (25 individuals from each ethnolinguistic group). Our scientific objectives are to: 1) develop a resource that provides a comprehensive catalogue of genetic variation in populations from SSA accessible to the global scientific community; 2) characterise population genetic diversity, structure, gene flow and admixture across SSA; 3) develop a cost-efficient, next-generation genotype array for diverse populations across SSA; and 4) facilitate whole genome-sequencing association studies of complex traits and diseases by developing a reference panel for imputation and resource for enhancing fine-mapping disease susceptibility loci. These scientific objectives will be supported by cross-cutting operational activities, including network and management of the consortium, research ethics, and research capacity building in statistical genetics and bioinformatics This data is part of a pre-publication release. For information on the proper use of pre-publication data shared by the Wellcome Trust Sanger Institute (including details of any publication moratoria), please see http://www.sanger.ac.uk/datasharing/ . This dataset contains all the data available for this study on 2021-02-16.

  Dataset EGAD00001006970

• Genome‐wide postnatal changes in immunity following fetal inflammatory response

  The fetal inflammatory response (FIR) increases the risk of perinatal brain injury, particularly in extremely low gestational age newborns (ELGANs, < 28 weeks of gestation). One of the mechanisms contributing to such a risk is a postnatal intermittent or sustained systemic inflammation (ISSI) following FIR. The link between prenatal and postnatal systemic inflammation is supported by the presence of well‐established inflammatory biomarkers in the umbilical cord and peripheral blood. However, the extent of molecular changes contributing to this association is unknown. Using RNA sequencing and mass spectrometry proteomics, we profiled the transcriptome and proteome of archived neonatal dried blood spot (DBS) specimens from 21 ELGANs. Comparing FIR‐affected and unaffected ELGANs, we identified 782 gene and 27 protein expression changes of 50% magnitude or more, and an experiment‐wide significance level below 5% false discovery rate. These expression changes confirm the robust postnatal activation of the innate immune system in FIR‐affected ELGANs and reveal for the first time an impairment of their adaptive immunity. In turn, the altered pathways provide clues about the molecular mechanisms triggering ISSI after FIR, and the onset of perinatal brain injury.

  Study EGAS00001003635

• Pan-neuroblastoma analysis reveals age- and signature-associated driver alterations

  Neuroblastoma is a pediatric malignancy with heterogeneous clinical outcomes depending on age and other factors. To better understand neuroblastoma pathogenesis, here we analyze whole-genome, whole-exome and/or transcriptome sequencing data from 702 neuroblastoma samples, half of which are new data. Over 95% of samples harbor at least one recurrent driver alteration and most aberrations, including MYCN, ATRX, and TERT alterations, differ in frequency by age. MYCN alterations occur at median 2.3 years of age, TERT at 3.8 years, and ATRX at 5.6 years, suggesting age-specific susceptibility to genetic alterations. Mutational signature analysis shows that reactive oxygen-species (ROS)-induced mutations are the most common cause of single nucleotide variant (SNV) drivers in neuroblastoma, including most ALK and Ras pathway SNVs. ROS-induced mutagenesis is both an early event, as it causes truncal variants, and continuous as it also induces relapse-specific mutations. ROS-induced mutations are more abundant in neuroblastomas with MYCN amplification, 17q gains, and increased expression of mitochondrial ribosome genes. We observe recurrent FGFR1 N546K or internal tandem duplication variants in 6 patients and ALK N-terminal structural alterations in 5 samples, identifying additional patients amenable to precision therapy.

  Study EGAS00001003931

• SINGLE-CELL RNA SEQUENCING Single-cell RNA sequencing was performed on 13 ‘mild-moderate’ and 10 ‘critical’ COVID19 PBMC samples

  Summary Epidemiological and clinical reports have indicated that SARS-CoV-2 virulence hinges upon the triggering of an aberrant host immune response, more so than on direct virus-induced cellular damage. To elucidate the immunopathology underlying COVID-19 severity, we performed cytokine and multiplex immune profiling in mild-moderate and critically ill COVID-19 patients. Hypercytokinemia in COVID-19 differed from the IFN-γ-driven cytokine storm in macrophage activation syndrome, and was more pronounced in critical versus mild-moderate COVID-19. Systems modelling of cytokine levels paired with deep-immune profiling showed that classical monocytes drive this hyper-inflammatory phenotype and that a reduction in T-lymphocytes correlates with disease severity, with CD8+ cells being disproportionately affected. Expression of antigen presenting machinery was also reduced in critical disease. Furthermore, we found that neutrophils contributed to disease severity and local tissue damage by amplification of hypercytokinemia and the formation of neutrophil extracellular traps. Together these findings suggest a myeloid-driven immunopathology, in which hyperactivated neutrophils and an ineffective adaptive immune system act as mediators of COVID-19 disease severity.

  Study EGAS00001005039

• RNA-sequencing of tumors from 45 patients with recurrent or metastatic gastric cancer treated with immune checkpoint inhibitors

  Genomic profiling can provide prognostic and predictive information to guide clinical care. Biomarkers that reliably predict patient response to chemotherapy and immune checkpoint inhibition in gastric cancer are lacking. In this work, we used our machine learning algorithm NTriPath to identify a gastric-cancer specific 32-gene signature. Using unsupervised clustering on expression levels of these 32 genes in tumors from 567 patients, we identified four molecular subtypes that were prognostic for survival. We then built a support vector machine with linear kernel and the binary classifier to generate a risk score that is prognostic for 5-year overall survival and validated the risk score using three independent datasets. We also found that the molecular subtypes predicted response to adjuvant 5-fluorouracil and platinum after gastrectomy and immune checkpoint inhibitors in patients with metastatic or recurrent disease. The 32-gene signature is a promising prognostic and predictive biomarker to guide the clinical care of gastric cancer patients and should be validated in a prospective manner. This EGA archive includes both St.Mary hospital cohort samples (n=28) and Yonsei immunotherapy patient samples (n=17).

  Study EGAS00001005588

• HTAN MCL Pre-Cancer Atlas Pilot Project - Targeted Sequencing Development Study

  The HTAN-MCL Pre-Cancer Atlas Pilot Project (PCAPP) is the result of a collaboration between the seven members of the MCL consortium. Across four organ types, PCAPP's goal is to collect and profile pre-malignant lesions for gene expression, DNA mutations, single-cell gene expression and immune-environment. Most PML are small in size and only available come from formalin fixed paraffin embedded archived tissue. The primary goal of PCAPP is to 1) understand the logistical challenges of PML specimen collection, 2) document technical limitations of the assays that are specific to the PML and 3) overcome them to support the generation of a more comprehensive Pre-Cancer Atlas in the future. The current upload provides RNA and DNA sequencing from participants with DCIS who were studied at the University of San Diego and the University of Vermont. Description of the overall study: A. Background/Significance One of the critical barriers to developing new approaches for cancer detection and prevention is the lack of understanding of the key molecular and cellular changes that cause cancer initiation and progression. Unlike the extensive work that has been done profiling advanced stage tumors, few studies have comprehensively profiled the molecular alterations found in precancerous tissues. Premalignant lesions are currently characterized by histologic changes that precede the development of invasive carcinoma1,2.These lesions can often be identified in regions surrounding an invasive tumor or in biopsies taken from patients undergoing diagnostic evaluation for suspicion of cancer. Currently, limited metrics exist to identify lesions that will likely progress to carcinoma and require intervention from those that will naturally regress or remain stable3,4. Characterization of the molecular alterations in premalignant lesions and the corresponding changes in the microenvironment would hasten the development of biomarkers for early detection and risk stratification as well as suggest preventive interventions to reverse or delay the development of cancer. Our pilot study will establish the feasibility of transcriptomic, genomic and immune profiling of FFPE premalignant lesions from multiple organ sites, collected and profiled with uniform SOPs across multiple institutions within the MCL consortium. We will characterize the molecular alterations in precancerous lesions and the corresponding microenvironment in four major organ sites, in order to uncover the molecular and cellular determinants of premalignancy, and establish standardized sequencing and immunohistochemistry protocols on FFPE precancerous tissue. We will also evaluate the technical feasibility of single nuclei sequencing of small FFPE pre-cancer lesions. Successful completion of the proposed pilot study will set the stage for expansion and development of a comprehensive Pre-Cancer Atlas (PCA) as part of the NCI's moonshot.B. Specific Aims Aim 1: Collect premalignant lesions (PML) and their associated microenvironment via LCM from FFPE tissue across four organ sites (breast, lung, pancreas & prostate). Aim 2: Perform bulk RNA and DNA seq on premalignant FFPE samples (and flash frozen tissue where available) and compare the genomic/transcriptomic alterations within and across organ sites. C. ApproachAim 1: Collect premalignant lesions (PML) and their associated microenvironment via LCM from FFPE tissue across four organ sites (breast, lung, pancreas & prostate). MethodsI. Patient Population/Sample Collection: Overview of the sites collecting PML tissue from the respective organs is provided in Table 1 and a full description of the biospecimens to be obtained is described in detail for each organ type below. Table 1. Breakdown of cohort by tissue type and collection site.Organ siteBreastLungPancreasProstateType of PMLDCISAAH, Squamous Dysplasia/CISIPMNsPINCollection of PatientsUCSF/UCSDUVMBU*/UCLAVanderbilt/MoffittMDACC*JHUStanford*# of Patients201920 (10 of each type)20 (10 of each type)242020Total patients per Organ39402440Note: single nuclei/cell RNA-Seq will be performed on 4-5 FFPE samples from each of the organ types 1. DCIS lesions from breast tissue: DCIS lesions will be collected from 39 patients (20 from UCSF/UCSD & 19 from UVM) with primary low or high-grade DCIS diagnosed from a breast core biopsy. Subsequent resected lumpectomy or mastectomy tissues will be prospectively sampled in the vicinity of the prior biopsy site using multiple approaches: 1) Live cells (heterogeneous mix) will be obtained as a cell scrape slurry from the lesion surface or by fine needle aspirate (FNA); 2) For a subset of specimens where size is sufficient, a block of breast tissue with DCIS will be fresh-frozen; 3) The remainder of the specimen will be taken for routine formalin-fixation and paraffin-embedding (FFPE). The FFPE sample will be annotated to identify the matched FFPE tissue block adjacent to the fresh-frozen sample and will be sectioned for use in bulk and single nuclei sequencing . We will dissect DCIS, adjacent normal and when available, associated carcinoma. In addition, when possible, normal tissue will be collected from a tissue block lacking lesions as well as collection of blood. A subset of patients (n = 5 | FFPE, flash frozen and fresh) will be sent to the Broad Institute for single nuclei/cell sequencing.2. AAH and squamous dysplastic/CIS lesions from airway and lung tissue: For squamous cell lung cancer, we will collect endobronchial biopsies from abnormal airway regions identified on autofluroscence bronchoscopy or identify PMLs in the margins of resected lung tissue. We will study 20 patients (5 each from BU/UCLA/Vanderbilt/Moffitt) with pre-invasive squamous lesions (moderate-severe dysplasia or carcinoma in situ (CIS)) identified on pathologic examination. LCM of the premalignant region and adjacent normal epithelium will be performed as well as the invasive tumor for those collected from the resection margin (n=5 from UCLA). On a subset of lesions collected at bronchoscopy (n=5), we will collect additional biopsies that will be flash frozen and fresh for single nuclei and cell sequencing, respectively, performed at the Broad Institute. In parallel to the work at the Broad, BU will perform single cell RNA-seq on these freshly cell sorted tissues (n = 5). Blood will be collected on all patients for genomic studies. For lung adenocarcinoma, we will collect resected FFPE lung tissues from 20 patients (10 from UCLA and 10 from Vanderbilt/Moffitt) with early stage lung adenocarcinoma that harbor atypical adenomatous (AAH) premalignant lesions in the resection margin. We will LCM multiple AAH regions (3-5 per patient) as well as adjacent regions of normal epithelium and invasive adenocarcinoma. In addition, blood will be collected on all patients for genomic studies. 3. IPMNs from pancreatic tissue: For pancreatic cancer PML, we will collect low and high grade lesions from 24 patients representing macroscopic Intraductal Papillary Mucinous Neoplasms (IPMN) (n=24) from surgically resected specimens along with blood samples. Archival FFPE specimens of microscopic PanIN lesions, occurring multi-focally adjacent to invasive PDAC, and archival IPMN lesions (with or without associated invasive cancer), along with the adjacent normal tissue, will undergo LCM and utilized for bulk DNA and RNA sequencing. If matched frozen tissues are available for a subset of these FFPE samples, we will bank for comparison of profiles. Because IPMNs are macroscopic lesions, they provide an opportunity for obtaining the samples fresh and therefore can be used for single cell sequencing (in contrast to PanINs). Therefore, 5 freshly obtained IPMNs will be used for the single cell RNA sequencing studies performed at both the Broad Institute and MDACC, and the matched FFPE and/or frozen sections from these lesions (obtained from the adjacent PML) will be sent to Broad Institute as a pilot to assess "single nuclei" RNA sequencing.4. PINs from prostate tissue: For prostate cancer PML, there will be 40 samples of Prostatic Intraepithelial Neoplasia (PIN) collected between the Stanford and JHU sites (20 cases per site). At the Stanford site, 20 prostate specimens detected by PSA screening who have/will undergo surgery (radical prostatectomy) for clinically localized disease will make up the final cohort. The age range of the participants would be 40-75, and we anticipate that 18 will be Caucasian, 1 Asian and 1 Latino or African American based on the practice demographics practice at Stanford. Clinical and MRI data will also be collected for these samples. We will collect low grade (e.g. Gleason score of 6/Grade group 1; n=10) and high grade (Gleason score 4+3=7 or higher/Grade group 3 or higher; n=10) PINs from FFPE samples that have prostate carcinoma. In addition to obtaining LCM archival samples of low and high grade PIN, we will also obtain normal prostatic epithelial from the peripheral, central and transition zones as well as multiple samples of prostate carcinoma in order to obtain the spectrum of Gleason grades in the carcinoma as needed. LCM samples will be used for bulk DNA and RNA sequencing. In addition, single cells will be dissected from FFPE samples to prepare single cell RNA seq libraries using techniques developed at Stanford, and FFPE tissue will be sent to the Broad for single nuclei sequencing. When available, flash frozen and fresh samples from these prostates will be archived and prepared for single nuclei and cell sequencing, respectively, at the Broad Institute and at Stanford (single cell only). JHU will also capture 10 cases (5 grade group 1 and 5 grade group 2) of high grade PIN, normal and invasive adenocarcinoma using frozen sections from fresh frozen tissues. When possible these will be from the same patients as the FFPE samples. Since frozen sections can be quite challenging to morphologically determine high grade PIN from normal epithelium, for these samples we will perform a number of additional tissue-based characterizations. These will include a multicolor combined basal cells (p63 and CK903) and PIN/carcinoma markers (AMACR) referred to in the cocktail as "PIN4", c-MYC (referred to as MYC) protein5, by IHC and mRNA by in situ hybridization (AM De Marzo, Q Zheng unpublished observations), telomere length by in situ hybridization6 and the 5'ETS/45S rRNA7. For these slides, the whole slides will be scanned with a Hammamatsu Nanozoomer with a 40x objective and regions of interest will be annotated as a guide for LCM.II. Laser-capture microdissection (LCM): FFPE tissue blocks will be sectioned at 7μm thickness and serial sections will be stained with H&E. LCM will be performed utilizing standard LCM systems, such as Leica LMD7000 and ArcturusXT at each site. Regions of premalignancy will be dissected and RNA/DNA will be extracted from microdissected cells using the Qiagen All Prep DNA/RNA FFPE Kit. Aim 2: Perform bulk RNA and DNA seq on premalignant FFPE samples and compare the genomic/transcriptomic alterations within and across organ sites. Rationale: There have been limited studies characterizing the genomic and transcriptomic landscape of premalignant lesions associated with breast, pancreatic, lung or prostate cancers. Characterizing the molecular determinants of premalignant disease that are unique and shared across multiple organs will enable new candidate biomarkers for early detection and novel therapeutic strategies for early intervention. MethodsBulk RNA-seq of LCM FFPE tissue: All participating sites will perform bulk RNA-seq in accordance with SOPs developed at BU. In brief, total RNA will be isolated from LCM'd lesion and associated microenvironment tissue using the Qiagen All Prep DNA/RNA FFPE Kit and quality will be assessed with the Agilent Bioanalyzer. Libraries will be generated with the Illumina TruSeq Access kit (for FFPE samples). They will be sequenced on the Illumina HiSeq2500 with 75base-pair paired-end reads. Quality of FASTQ files will be assessed with FastQC. Reads will be aligned to the human genome with STAR and gene-level and isoform-level expression will be quantified with RSEM. Splice junction saturation, transcript integrity, and biotype distributions will be calculated for each sample with RSeQC. DESeq2 and EdgeR will be used to identify associations between gene expression profiles and clinical variables while controlling for confounding covariates. BU will serve as an RNA-seq Core to assess reproducibility of FFPE RNA-seq methods across sites. We will perform RNA-seq according to the SOP listed above on a subset of samples for each organ type (total n ~ 20). Bulk seq of DNA from FFPE tissue: All participating sites will perform targeted or whole exome-seq (WES) in accordance with SOPs. In brief, DNA from laser captured material will be isolated using the Qiagen All Prep DNA/RNA FFPE Kit and undergo stringent quality control to ensure high quality input material for genomic profiling. Purified DNA (ideally 100-200 ng) will be used for library preparation and amplification, followed by next generation sequencing using standard protocols distributed by CDMG. Exome-seq methods are considered standardized, thus we will not need a DNA-seq Core to assess reproducibility across sites. We anticipate local centers will use Illumina paired end reads, following the following general approach. 1) DNA library preparation: Paired-end libraries will be prepared following the manufacturer's protocols (Illumina and Agilent), fragmented to 150-200 bp 2) Capture of targeted exome: Whole exome capture will be carried out using the protocol for Agilent's SureSelect Human All Exon kit. Purified capture products will be amplified using the SureSelect GA PCR primers (Agilent) for 12 cycles. 3) Sequencing will be carried out for the captured libraries using at least 100 bp paired-end reads. To achieve high level sensitivity and accuracy for detecting all the mutations in the whole exome, each sample will be sequenced at 200X mean depth. 4) Read mapping and alignment and variant analysis: Sequence short reads will be aligned to a reference genome (NCBI human genome assembly build 38) using BWA-MEM. Local realignment of aligned reads will be performed using Genome Analysis Toolkit (GATK).Data QC: To ensure scientific rigor and consistency among sites in RNA and DNA processing we will include a preliminary analysis of steps in processing and analysis. Protocols for extraction of high quality RNA and DNA from formalin fixed paraffin embedded (FFPE) tissues, which will be used extensively in these studies continue to improve and may have variable implementation among the sites participating in this study. To evaluate consistency of preliminary steps in processing and downstream analyses, we will initially distribute slides from one large FFPE fixed cancer of origin from prostate, breast, lung and pancreatic cancer. Analysis of these samples will allow us to review the DNA and RNA characteristics (yield, purity and strand length) among sites. Downstream analysis of these same samples will also allow us to compare among sites the consistency of variant calls among centers. We will be able to identify if there are some times of calls (such as small insertion deletions) that are more variable among centers versus other types of calls (such as relative gene expression or single base pair substitutions) that we expect to be less variable and to characterize the reliability of findings across sites. We are also including a 5% blind duplicate analysis of RNA sequencing. Samples will be analysed by the participating genomics cores without knowledge of the phenotype. RNA seq and CNA analyses are normalized for batch effects. We will also compare the observed sex to the self-reported sex as based on RNA profiles and exome sequencing of X chromosome genes as another check for processing accuracy and sample management. D. References 1. Wacholder, S. Precursors in Cancer Epidemiology: Aligning Definition and Function. Cancer Epidemiol. Prev. Biomark. 22, 521-527 (2013). PMID: 23549395.2. Berman, J. J. Precancer: The Beginning and the End of Cancer. (Jones & Bartlett Learning, 2011).3. Nasiell, K., Nasiell, M. & Vaćlavinková, V. Behavior of moderate cervical dysplasia during long-term follow-up. Obstet. Gynecol. 61, 609-614 (1983). PMID: 6835614.4. Merrick, D. T. et al. Persistence of Bronchial Dysplasia Is Associated with Development of Invasive Squamous Cell Carcinoma. Cancer Prev. Res. (Phila. Pa.) 9, 96-104 (2016). PMID: 26542061.5. Gurel, B. et al. Nuclear MYC protein overexpression is an early alteration in human prostate carcinogenesis. Mod. Pathol. Off. J. U. S. Can. Acad. Pathol. Inc 21, 1156-1167 (2008). PMID: 18567993.6. Meeker, A. K. et al. Telomere shortening is an early somatic DNA alteration in human prostate tumorigenesis. Cancer Res. 62, 6405-6409 (2002). PMID: 12438224.7. Guner, G. et al. Novel Assay to Detect RNA Polymerase I Activity In Vivo. Mol. Cancer Res. MCR 15, 577-584 (2017). PMID: 28119429.

  Study phs002225

• Resource for Genetic Epidemiology Research on Adult Health and Aging (GERA)

  The Resource for Genetic Epidemiology Research on Aging (GERA) Cohort was created by a RC2 "Grand Opportunity" grant that was awarded to the Kaiser Permanente Research Program on Genes, Environment, and Health (RPGEH) and the UCSF Institute for Human Genetics (AG036607; Schaefer/Risch, PIs). The RC2 project enabled genome-wide SNP genotyping (GWAS) to be conducted on a cohort of over 100,000 adults who are members of the Kaiser Permanente Medical Care Plan, Northern California Region (KPNC), and participating in its RPGEH. The purpose of the RPGEH is to facilitate research on the genetic and environmental factors that affect health and disease by linking together clinical data from electronic health records, survey data on demographic and behavioral factors, and environmental data from various sources, with genetic data derived from biospecimens collected from participants. At the time of the award of the RC2 project in late 2009, the RPGEH had established a cohort of about 140,000 individuals who had answered a detailed survey, provided saliva samples for extraction of DNA, and given broad consent for the use of their data in studies of health and disease. To maximize the diversity of the resulting sample, the GERA cohort was formed by including all racial and ethnic minority participants with saliva samples (N = 20,925; 19%); the remaining participants were drawn sequentially and randomly from white non-Hispanic participants (89,341; 81%). A total of 110,266 participant samples were included to ensure that at least 100,000 were successfully assayed. The resulting GERA cohort is 42% male, 58% female, and ranges in age from 18 to over 100 years old with an average age of 63 years at the time of the RPGEH survey (2007). The sample is ethnically diverse, generally well-educated with above average income. Approximately 69% of the participants are married or living with a partner. Length of membership in KPNC averages 23.5 years. UCSF and RPGEH investigators worked with the genomics company Affymetrix to design four custom microarrays for genotyping each of the four major race-ethnicity groups included in the GERA Cohort, described in detail in Hoffmann et al., 2011a and 2011b. Following genotyping and quality control procedures, and after removal of invalid, discordant, or withdrawn samples, about 103,000 participants were successfully genotyped. The resulting genotypic data were linked to survey data and data abstracted from the electronic medical records. As described below, all RPGEH participants were mailed new consent forms with explicit discussion of the placement of data in the NIH-maintained dbGaP. About 77% of participants returned completed consent forms, resulting in a final sample size of 78,486 participants in the GERA Cohort with data for deposit into dbGaP. Origins of the RPGEH GERA Cohort The goal in creating the RPGEH GERA cohort was to create a large, multiethnic, and comprehensive population-based resource for research into the genetic and environmental basis of common age-related diseases and their treatment, and factors influencing healthy aging and longevity. The GERA Cohort consists of a diverse cohort of more than 100,000 adults who are members of the Kaiser Permanente Medical Care Plan, Northern California Region (KPNC), and participating in its Research Program on Genes, Environment and Health (RPGEH). KPNC is an integrated health care delivery system with a population of about 3.3 million people in northern California. The membership of KPNC is representative of the general population in the 14 county area in which facilities are located, although the membership is underrepresented for the extremes of income at both ends of the spectrum. The RPGEH utilizes the longitudinal electronic health records (EHR) of KPNC to obtain clinical, laboratory, imaging and pharmacy information on all cohort members, to which personal demographic, behavioral and health characteristics have been added through member surveys. The GERA Cohort comprises a subsample of the RPGEH participant cohort, and was created through the RC2 award from the NIA, NIMH, and NIH Common Fund as described above. GERA Study Design The GERA Cohort is a subsample, as described above, of the longitudinal cohort enrolled in the Kaiser Permanente RPGEH. The RPGEH cohort includes about 400,000 survey participants of whom about 200,000 have provided broad consent and a sample of saliva or blood for use in studies of genetic and environmental factors in health and disease. The GERA Cohort was developed from a mailed survey sent to all adult members of KPNC who had been members for two years or more in 2007. All survey respondents were contacted and asked to complete a consent form; those who completed consent forms were asked to provide a saliva sample. Additional male participants were added to the RPGEH through inclusion of the Northern California sample of the California Men's Health Study (CMHS) cohort of about 40,000 men from KPNC, ages 45-69 years old at the time of the CMHS survey in 2002-2003. The CMHS participants contributed about 15,400 saliva samples to the RPGEH and were eligible for inclusion in the GERA Cohort. CMHS participants were included according to the same sampling design as for the RPGEH cohort as a whole. Specifically, all minority participants were selected for inclusion in order to maximize representation of minorities in the GERA Cohort, and Non-Hispanic White participants were selected at random to complete the sample of 110,266 GERA Cohort participants. GERA Genotypic Data High-density genotyping was conducted at UCSF using custom designed Affymetrix Axiom arrays, as described in Hoffmann et al. (2011a; 2011b). To maximize genome-wide coverage of common and less common variants, four specific arrays were designed for individuals of Non-Hispanic White (EUR), East Asian (EAS), African-American (AFR), and Latino (LAT) race/ethnicity. There was broad overlap among the SNPs on the arrays, which were designed using a hybrid greedy imputation algorithm (Hoffmann et al., 2011b) applied to genotype information validated by Affymetrix from the 1000 Genomes Project. However, in order to capture low frequency variants specific to particular race-ethnicity groups, SNP content varies between arrays. A more detailed description of the process of genotyping and results is included in Genotyping of DNA Samples. Description of the analyses of population structure and development of principal components for adjustment of population structure is included in Population Structure Analysis. GERA Phenotypic Data RPGEH and CMHS Survey Data. The sources of data on demographic and behavioral factors deposited in dbGaP for the GERA Cohort are the RPGEH and CMHS surveys. Data on common demographic factors such as gender, race/ethnicity, marital status, and education and on behavioral factors such as smoking, alcohol consumption, and body mass index, have been cleaned, edited, reconciled between the two surveys, and compiled into summary indices, where appropriate, for deposition into dbGaP. A more complete description of the survey variables is included in Survey Variables Documentation. Please note that the terms of use of the GERA Cohort Data, as specified in the Data Use Certification (DUC), prohibit the use of survey variables as outcomes in analyses. For example, a genome-wide association study (GWAS) of education or smoking is not permitted as specified by the DUC. Only health conditions can be used as outcome variables in analyses. Health Conditions derived from Kaiser Permanente Electronic Medical Records. Data on the occurrence of health conditions in participants in the GERA Cohort have been derived from summarizing ICD-9 coded diagnoses in Kaiser Permanente's electronic medical records. An algorithm that aggregates specific ICD-9 codes into appropriate diagnostic groups for selected conditions is applied to outpatient and inpatient databases; see Disease and Conditions Definitions Documentation for details. The criterion for including a condition as "present" for a participant is the occurrence of two or more diagnoses within a diagnostic category occurring on separate days. Two or more is used as the criterion in order to reduce false positives due to mistakes or rule-out diagnoses. When compared with validated disease registries, the criterion of 2+ diagnoses yields high specificity and good sensitivity. ICD-9 codes in the electronic records are specified in several ways. For outpatient visits occurring during the period 1995 to 2006, diagnoses were assigned by the treating physician who endorsed specific diagnoses on an optically scanned list that varied by specialty. Beginning in 2006 with the advent of an integrated, fully electronic medical record, outpatient diagnoses are made by physicians/ providers using a pull down menu. Discharge diagnoses from inpatient stays are specified by physicians and coded by specially trained coders. Databases of ICD-9 codes for diagnoses assigned at outpatient visits, or as one of the discharge diagnoses following inpatient stays, are complete and available for all KPNC members dating back to 1995. Although the average length of KPNC membership among GERA cohort members is 23.5 years in 2007, not all have been members since 1995, so the history for some conditions, such as those that are not chronic or recurrent, may not be complete for all cohort members. The year of first membership in KPNC is included as a variable in the list of survey variables, enabling investigators to estimate the number of years of observation of each Cohort member. RPGEH Access and Collaborations Website and Procedures The RPGEH maintains a web portal for inquiries and applications for collaboration and access to data. The url is: https://rpgehportal.kaiser.org/. RPGEH has an application process and an Access Review Committee that reviews applications for collaboration and use. For more details, please contact RPGEH through the website.

  Study phs000674

• The Vaginal Microbiome: Disease, Genetics and the Environment

  The vagina is an interactive interface with the environment, and as such is covered by a protective epithelial surface. This surface is colonized by a biofilm consisting of bacteria and other microorganisms, which through a variety of mechanisms serve to protect the host from invasion by pathogens. Alterations in the normal vaginal microflora; i.e., dysbioses, particularly those associated with bacterial vaginosis, the most common disorder of the female reproductive tract, contribute to risk of sexually transmitted infections, infection by HIV, and are associated with increased risk of adverse pregnancy outcomes including but not limited to spontaneous preterm birth. Physiologic alterations in the host (e.g., menopause and pregnancy), which are beginning to be investigated as potential selective conditions for change in the "normal" vaginal flora, and their impact on disease susceptibility and transmission, remain to be more definitively elucidated. The effects of chronically abnormal physiologic states (e.g., diabetes mellitus) on normal vaginal flora have not been well described or studied. Finally, an almost unexplored area of inquiry is the genetic contribution, including race/ethnicity, to the establishment and maintenance of a "normal" vaginal flora, under normal and physiologically altered circumstances. Our research will shed light on how the vaginal microbiome contributes to adverse obstetrical outcomes and sexually transmitted infections in diverse populations. This project addresses the following questions: First: Do the genes of the host contribute to the composition of the vaginal microbiome? We hypothesize that a woman's genetic composition significantly affects the ability of certain commensal, parasitic and pathogenic microbes to colonize and/or infect the genital tract. Thus, we will compare and quantify the microbial populations inhabiting the vaginas of monozygotic and dizygotic twins from the Mid Atlantic Twin Registry, and, in the process, address the question of whether there is a relationship between the microbiomes of the vagina, mouth and GI tract. Second: What changes in the vaginal microbiome are associated with common physiological perturbations or non-infectious pathological states of the host?We hypothesize that "altered" physiologic (pregnancy, menopause) and pathologic (chronic disease, hysterectomy) conditions, or environmental "exposures" (exogenous hormones, antibiotics, chronic immunosuppressant, smoking; douching) can predictably alter the vaginal microenvironment. These alterations will lead to changes in microbial populations within the vagina. Changes in the microbial populations may have impacts, positive or more likely negative, on the spontaneous and future well-being of the affected individual. We are characterizing the effects of these "altered" physiologic and pathologic conditions, and environmental exposures, on the composition of the vaginal microbiome. Finally, we are examining the impact of the important physiological changes incurred during pregnancy on the composition of the microbiome of the female reproductive tract. Third: What changes in the vaginal microbiome are associated with relevant infectious diseases and conditions? We are testing the hypothesis that infectious diseases predictably alter the vaginal microbiome, and that these changes have an impact on the disease susceptibility, process, and outcome. A predilection for bacterial vaginosis, vaginitis, HIV infection, or other sexually transmitted diseases, is likely associated with a women's vaginal microbial composition. Thus, we will characterize samples from women with a variety of these infectious conditions to determine the contribution of their microbiomes to the disease process and susceptibility. We are addressing these questions using a combination of high throughput 'nextgen' sequencing technologies. Early studies were performed using the Roche 454 FLX and the first generation Illumina Genome Analyzer II instruments and more recently employing Illumina MiSeq and HiSeq 2500 and 4000 platforms, installed in the Genomics Core of the Nucleic Acids Research Facilities at VCU. Thus, segments of the 16S rRNA genes are amplified from the complex samples taken from various target sites in and around the vagina of each study participant. These segments are sequenced and subjected to taxonomic classification protocols to identify and quantify the bacterial taxa present in each sample. Additionally, total DNA isolated from these samples are subjected to shotgun sequence analysis to empirically reconstruct the metabolic potential of these microbiomes. Finally, specific bacterial clones are completely sequenced and analyzed to associate the unavoidable strain and isolate diversity with the clinical phenotypes presented.

  Study phs000256

• Human Microbiome Project Demonstration Study of Cutaneous Microbiome in Psoriasis

  Psoriasis, a highly prevalent disease of humans of unknown cause, is a chronic inflammatory disorder primarily involving skin, with distinctive clinical characteristics. With the newly developed tools that facilitate microbiome research, it now is possible to assess whether the cutaneous microbiome plays a role in the pathogenesis of this disorder. Preliminary data from our studies suggest that the cutaneous microbiome in psoriasis is complex and possibly different from normal. To deal with this complexity, we propose to examine the cutaneous microbiome in relation to psoriasis with explorations at several taxonomic and informatic levels. Our overall objective is to examine how changes in the normal cutaneous microbiome contribute to the pathogenesis of psoriasis. Since causality is complex and often difficult to prove, and beyond the scope of this RFP, our overall hypothesis is that there are alterations in the cutaneous microbiome in areas of skin affected by psoriasis in comparison with the range observed in clinically unaffected areas, or in healthy persons. We also hypothesize that the characteristics of the microbiome may affect clinical responses to the immunomodulatory agents used to treat psoriasis. An alternative hypothesis is that effective treatment of psoriasis with systemic immunomodulatory agents will not substantially affect the disordered microbial ecosystem. Such observations would provide evidence for the roles of the microbiota in this disorder. Since an important consideration in microbiome research is the optimal level (e.g. phylum, genus, species, strain, gene) at which to examine a scientific question, and we are not yet certain what are the optimal levels for psoriasis, this also will be examined. Our studies of psoriasis should allow development of both approaches and tools that will have general utility for Microbiome research. To test our hypothesis, we propose the following specific aims: 1) To understand the cutaneous microbiome species composition overlaying psoriatic lesions; 2)To investigate differences in metagenome content for psoriatic lesions compared to normal skin; 3) To identify differences in the transcriptional profiles of the microbiome and the host between normal skin and psoriatic lesions using high-throughput sequencing; and 4) To estimate the effects of systemic immunomodulatory therapy for psoriasis on microbiome composition. In total, these studies should help us understand the role of the microbiome in psoriasis pathogenesis. We sought to characterize and compare the cutaneous microbiota of psoriatic lesions (lesion), unaffected contralateral skin from psoriatic patients (normal), and similar skin loci in matched healthy controls (control) in order to discern patterns that govern skin colonization and their relationship to clinical diagnosis. Using high-throughput 16S rRNA sequencing, we assayed the cutaneous bacterial communities of 51 matched triplets and characterized these samples using community data analysis techniques.

  Study phs000251

• Genome Variation among HIV-Resistant People with Hemophilia

  Rare individuals are highly resistant to infection with human immunodeficiency virus (HIV). Studies of candidate genes resulted in the discovery of a 32bp deletion in the CC-chemokine receptor 5 gene (CCR5Δ32), which rendered this critical co-receptor for primary HIV infection to be non-functional. Pharmacologic and vaccine-induced blockade of CCR5 is being pursued to treat and prevent HIV infection and other conditions. The allele frequency of CCR5Δ32 among persons of European ancestry is approximately 10%. CCR5Δ32/Δ32 homozygotes are almost totally resistant to HIV infection. People with severe hemophilia A require frequent replacement with clotting Factor VIII (FVIII) to control hemorrhage. Prior to the discovery of HIV in 1984 and licensure of recombinant FVIII in the late 1980s, people with severe hemophilia A were treated with plasma-derived FVIII and thus were intensively exposed to HIV. Only 5% of such patients were not infected with HIV. Of these, approximately 1/3 were CCR5Δ32/Δ32 homozygotes. The remaining 2/3 of these people who were highly resistant to HIV remain unexplained. This project seeks to discover genome variations among people who are highly resistant to HIV infection. Such variation is likely to serve as a target for reducing the morbidity and incidence of HIV.

  Study phs000445

• Analysis of the Whole Transcriptomes of Human Testis with Complete Spermatogenesis and of Human Testes with Sertoli Cell-Only Syndrome

  Sertoli cell-only syndrome is severe form of human male infertility in which most seminiferous tubules appear to lack all spermatogenic cells, including spermatogonial stem cells (SSCs). However, a few small tubule segments of some patients have active spermatogenesis and, thus, functional stem cell niches and SSCs. Normally SSCs replicate, migrate and refill adjacent empty niches, but this does not appear to occur in SCO syndrome. We hypothesized that this failure occurs because most niches are dysfunctional. As Sertoli cells are essential to formation of these niches, we used RNAseq to compare the transcriptomes of human testes with qualitatively normal (complete) spermatogenesis (n=4) with the transcriptomes of human testes with SCO syndrome (n=7). We then focused our analysis on the expression of transcripts that bioinformatic analyses identified as Sertoli cell signature transcripts. Results show that Sertoli cells in SCO testes express abnormally low levels of GDNF, FGF8 and BMP4, all of which are important regulators of mouse SSCs and/or progenitor spermatogonia. Sertoli cells in SCO testes express significantly reduced levels of transcripts for proteins that polarize the Sertoli cell plasma membrane and regulate the trafficking of cell adhesion and gap junction proteins in and out of that plasma membrane.

  Study phs001777

• Whole Transcriptome Sequencing of Human Tumor Cells and Hematopoietic Stem and Progenitor Cells During Aging and Bone Marrow Disorders

  Whole transcriptome sequencing analyses were performed on human hematopoietic stem and progenitor cells during aging and bone marrow disorders. For this purpose, normal and diseased bone marrow or peripheral blood-derived mononuclear cells were FACS-purified into hematopoietic stem (CD34+CD38-Lin-) or hematopoietic progenitor (CD34+CD38+Lin-) populations. Primary human samples from healthy individuals were used along with comparative analyses in some pre-malignant and malignant conditions, including myelodysplastic syndrome (MDS), acute myeloid leukemia (AML) and plasma cell dyscrasias (e.g., multiple myeloma). Studies were also performed on total mononuclear cell fractions from patients with newly-diagnosed or relapsed multiple myeloma, or plasma cell leukemia. Briefly, RNA was extracted from all samples and quality was assessed using a Bioanalyzer instrument. Only samples with RNA integrity (RIN) values >7 were subjected to library preparation for RNA-sequencing. The SMARTer v4 cDNA amplification kit was used along with library construction by the NEB Ultra DNA kit for Illumina. Paired-end sequencing was performed on Illumina NextSeq instruments (Scripps Next Generation Sequencing Core), with all samples run over two separate flow-cells (to generate additional sequencing reads); >95M reads were generated in total for each sample. The RNA-sequencing data are available as FASTQ files.

  Study phs002291

• Germ Cell and Associated Heme Malignancies Evolve from a Common Shared Precursor

  Germ cell tumors (GCTs) are the most common cancer in men between the ages of 15-40. While most patients are cured, those with disease arising in the mediastinum have distinctly poor outcomes. One in every 17 patients with primary mediastinal non-seminomatous GCTs develop an incurable hematologic malignancy and prior data intriguingly suggests a clonal relationship exists between hematologic malignancies and GCTs in these cases. To date however, the precise clonal relationship between GCTs and the diverse additional somatic malignancies arising in such individuals has not been determined. Here, we traced the clonal evolution and characterized the genetic features of each neoplasm from a cohort of fifteen patients with GCTs and associated hematologic malignancies. We discovered that GCTs and hematologic malignancies developing in such individuals evolved from a common shared precursor, nearly all of which harbored allelically imbalanced TP53 and/or RAS pathway mutations. Hematologic malignancies arising in this setting genetically resembled mediastinal GCTs rather than de novo myeloid neoplasms. Our findings argue that this scenario represents a unique clinical syndrome, distinct from de novo GCTs or hematologic malignancies, initiated by an ancestral precursor which gives rise to the parallel evolution of GCTs and blood cancers in these patients. Reprinted from PMID: 32897884, with permission from JCI.

  Study phs002231

• TNF-induced dynamic regulation of mRNA stabilome in rheumatoid arthritis fibroblast-like synoviocytes

  During rheumatoid arthritis (RA), TNF activates fibroblast-like synoviocytes (FLS) inducing in a temporal order a constellation of genes, which perpetuate synovial inflammation. Although the molecular mechanisms regulating TNF-induced transcription are well characterized, little is known about the impact of mRNA stability on gene expression and the impact of TNF on decay rates of mRNA transcripts in FLS. To address these issues we performed RNA sequencing and genome-wide analysis of the mRNA stabilome in RA FLS. We found that TNF induces a biphasic gene expression program: initially, the inducible transcriptome consists primarily of unstable transcripts but progressively switches and becomes dominated by very stable transcripts. This temporal switch is due to: a) TNF-induced prolonged stabilization of previously unstable transcripts that enables progressive transcript accumulation over days and b) sustained expression and late induction of very stable transcripts. TNF- induced mRNA stabilization in RA FLS occurs during the late phase of TNF response, is MAPK-dependent, and involves several genes with pathogenic potential such as IL6, CXCL1, CXCL3, CXCL8/IL8, CCL2, and PTGS2. These results provide the first insights into genome-wide regulation of mRNA stability in RA FLS and highlight the potential contribution of dynamic regulation of the mRNA stabilome by TNF to chronic synovitis.

  Study phs001371

• Identification and Targeting of Inflammatory Macrophage-Fibroblast Crosstalk in Rheumatoid Arthritis

  Macrophages tailor their function according to the signals found in tissue microenvironments, assuming a wide spectrum of phenotypes. A detailed understanding of macrophage phenotypes in human tissues is limited. Using single-cell RNA sequencing, we defined distinct macrophage subsets in the joints of patients with the autoimmune disease rheumatoid arthritis (RA), which affects ~1% of the population. The subset we refer to as HBEGF+ inflammatory macrophages is enriched in RA tissues and is shaped by resident fibroblasts and the cytokine tumor necrosis factor (TNF). These macrophages promoted fibroblast invasiveness in an epidermal growth factor receptor-dependent manner, indicating that intercellular cross-talk in this inflamed setting reshapes both cell types and contributes to fibroblast-mediated joint destruction. In an ex vivo synovial tissue assay, most medications used to treat RA patients targeted HBEGF+ inflammatory macrophages; however, in some cases, redirecting them into a state is not expected to resolve inflammation. These data highlight how advances in our understanding of chronically inflamed human tissues and the effects of medications therein can be achieved by studies on local macrophage phenotypes and intercellular interactions. Reprinted from Kuo, Ding et al., Science Translational Medicine 2019, PMID: 31068444, with permission from American Association for the Advancement of Science.

  Study phs001340

• Genomic Sequencing of Cervical Cancers

  Cervical cancer is responsible for 10-15% of cancer related deaths in women worldwide. The etiological role of infection with high-risk human papilloma viruses (HPV) in carcinomas of the cervix is well established. In general, the development of cervical carcinomas follows a progression from persistent HPV infection through precancerous lesions to invasive cancer. Previous studies have implicated somatic mutations in PIK3CA, PTEN, TP53, STK11 and KRAS as well as chromosome-arm level copy number alterations in the pathogenesis of cervical carcinomas. Here, we report whole exome sequencing analysis of 118 cervical carcinoma-normal paired samples from patients in Norway and Mexico, as well as transcriptome sequencing of 80 cases and whole genome sequencing of 13 tumor-normal pairs. Novel somatic mutations include recurrent E322K substitutions in the MAPK1 gene encoding the ERK2 kinase and inactivating mutations in the HLA-B gene. In addition, recurrent somatic mutations in FBXW7, EP300, and NFE2L2 are novel in the context of primary cervical carcinomas. Analysis of HPV integration sites revealed recurrent integration into the RAD51B locus as well as co-occurrence of HPV genome integration and copy number gains within several genomic loci. These findings shed new light on the pathogenesis of cervical carcinomas and suggest potential novel therapeutic targets.

  Study phs000600

• Single Cell and Tissue Level Functional Genomics Analysis of Astrocyte-Related Mechanisms in Taupathy

  This data is the first part of a larger single cell and tissue level functional genomics analysis of dementias with tau pathology.  Cellular resolution will clarify the role of glial genes and pathways in neurodegenerative disease. This includes understanding the role of microglia and astrocytes in neurodegenerative disorders characterized by Tau protein pathology. To achieve a refined view of microglial neuroinflammatory signaling in tau-associated neurodegeneration, we performed a functional genomic analysis that integrated time-course gene expression data from mouse models expressing human mutant MAPT (the tau gene) with single nuclear sequencing data obtained from patients with frontotemporal dementia with Tau protein pathology (Pick’s disease). This included microglia-specific gene expression data obtained from single nuclear sequencing of postmortem frozen precentral gyrus samples collected from 16 total individuals, including 8 subjects with clinical diagnosis of behavior variant frontotemporal dementia and neuropathological diagnosis of Pick’s disease (Tau protein pathology) and controls (n=8). Our findings present detailed molecular profiles of disease-associated microglial transitions over the course of tau-associated neurodegeneration. Furthermore, we demonstrate overlap between microglia profiles obtained from precentral gyrus of patients with Pick’s disease with those up-regulated during the neurodegenerative stage in mice expressing human MAPT mutations.

  Study phs002197

• Whole genome sequencing of core-binding factor leukemia

  Pediatric de novo acute myeloid leukemia (AML) is a heterogeneous disease that can be divided into clinically distinct subtypes based on the presence of specific chromosomal abnormalities or gene alterations. One of the best characterized subtypes of AML involves leukemias with alterations of the core-binding factor (CBF)-complex, which comprises the FAB subtypes M2 and M4Eo and associates with a favorable outcome. Patients with the AML M2 subtype harbor a translocation between chromosomes 8 and 21 [t(8;21)] that yields the chimeric fusion gene RUNX1(AML1)-RUNX1T1(ETO), while patients with AML M4Eo express the chimeric fusion gene CBFβ-SMMHC(MYH11) as a result of an inversion/translocation event of chromosome 16 [inv(16)/t(16;16)]. In an effort to define the total complement of genetic changes in CBF-leukemia, we performed paired-end whole genome sequencing (WGS) on diagnostic leukemia blasts and matched germ line samples from 17 pediatric CBF-leukemia patients using the Illumina platform. Somatic alterations, including single nucleotide variations (SNVs) and structural variations (SVs), including insertions, deletions, inversions, and inter- and intra-chromosomal rearrangements, were detected using complementary analysis pipelines (Bambino, CREST and CONSERTING). Recurrent screening of identified mutations will be performed in a cohort of approximately 94 cases of CBF-leukemias.

  Study phs000414

• International Multi-Center ADHD Genetics Project

  The goal of the project is to complete a 600,000 tag SNP genome-wide association scan of 958 parent-child trios from the International Multisite ADHD Genetics (IMAGE) project, in order to assess the association of SNP markers with ADHD, analyze quantitative ADHD phenotypes, complete copy number analyses, assess parent of origin effects and season of birth effects, and test for epistasis among apparently uncorrelated genes. Acquiring DNA Samples All consent forms stipulate that the samples can only be used by researchers who have been approved by the National Institute of Mental Health (NIMH), National Institutes of Health. All consent forms, except those used at the Zürich site (N=141 subjects), explicitly indicate that the samples may be used by researchers from commercial enterprises seeking to benefit financially from the analysis of the samples. The Zürich consent does not prohibit such use. The Zürich consent form also included an "opt out" that allowed the subjects to indicate that they did not want their samples stored at the NIMH repository or used by researchers external to the project. No subjects enrolled in the project opted out. Consent groups and participant set ADHD (ADHD): 2758 (924 trios)

  Study phs000016

• Cross-Species Single-Cell Analysis of Pancreatic Ductal Adenocarcinoma Reveals Antigen-Presenting Cancer-Associated Fibroblasts

  Cancer-associated fibroblasts (CAFs) are major players in the progression and drug resistance of pancreatic ductal adenocarcinoma (PDAC). CAFs constitute a diverse cell population consisting of several recently described subtypes, although the extent of CAF heterogeneity has remained undefined. Here we employ single-cell RNA-sequencing to thoroughly characterize the neoplastic and tumor microenvironment content of human PDAC tumors. Six human PDAC tumor specimens from six patients were collected, and processed for single-cell RNA-sequencing analysis. Adjacent-normal pancreas tissue was also collected from two of the patients. Tumor samples were digested, and fluorescence-activated cell sorting was used to isolate viable cells. For one tumor sample, viable, CD45-negative, CD31-negative, and EpCAM-negative cells were also isolated to enrich for CAFs. The 10X Chromium platform was then used to isolate single cells for RNA-sequencing analysis. This work has demonstrated the differences in immune cell populations between adjacent-normal and tumor tissues, and identified subpopulations of epithelial cells and CAFs present in PDAC tumors. This high-throughput analysis is a resource to better understand the cell populations present in PDAC, and may ultimately aid in the development of more effective therapies for this deadly malignancy.

  Study phs001840

• Integrated Genomic Analysis of Periampullary Tumors at The Human Genome Sequencing Center, Baylor College of Medicine (HGSC-BCM)

  Adenocarcinomas arising in the complex environment of the ampulla of Vater constitute a histopathological heterogenous group, presumably originating from the different epithelial cellular constituents present at the site: pancreas, bile duct, and intestinal duodenum. These tumors have been described in many different ways: intra-ampullary, periampullary, intra-ampullary papillary-tubular neoplasm, ampullary-ductal, periampullary-duodenal, and ampullary-not otherwise specified. These varied classifications reflect the difficulty in classifying these tumors into specific groups. Only the tumors clearly localized in the bile duct or duodenum are identified as distal cholangiocarcinomas (CAC) or duodenal adenocarcinomas (DUOAC). The current classification is based on macroscopic features that may distinguish the epithelium of origin, microscopic features, clinicopathological criteria, histopathology and expression of differential markers. This classification is subjective and prone to inter-observer variability and significantly impacts on treatment selection and therapeutic development. In order to define subtypes of periampullary cancer with clinical relevance, we performed whole exome sequencing and copy number analysis of 160 cancers arising in the periampullary region, 62 of these clearly arising from either the bile duct (n = 44), or the duodenum (n = 18) and 98 periampullary cancers (AMPAC) where the epithelium of origin could not be clearly defined.

  Study phs000895

• Trisomy21: Risk Factors for Chromosome Nondisjunction (T21NDJ)

  The overall goal of the project is to identify genetic risk factors associated with chromosome 21 nondisjunction in the oocyte. The dataset derives from multi-site collection of live birth probands with Down syndrome due to standard trisomy 21 (T21) and their biological parents. The type of nondisjunction error (NDJ)(maternal or paternal) in all cases has been determined to be maternal in origin based on the chromosome 21 variants contributed from parent to proband. The Center of Inherited Disease with support from NICHD has conducted genome-wide genotyping using the Illuminia Human OmniExpress Plus Exome array on approximately 800 women who have been identified through their offspring with DS and have been characterized as having a maternal meiois I (MI) or meiosis II (MII) nondisjunction error. Genotypes from biological fathers of the offspring with DS can be used with the data on mothers to better define the type of nondisjunction error (MI or MII) and to refine the chromosome 21 recombination profile. We provide the type of nondisjunction error, knowing that this will be updated based on the new panel. We do not provide the recombination profile, as this can be best defined using the new comprehensive set of SNPs in the OmniExpress panel.

  Study phs000718

• NIDDK IBD Genetics Consortium Repository Immunochip

  The Immunochip is a custom Illumina Infinium chip comprising 196,524 SNPs and small indels selected primarily based on GWAS analysis of 12 autoimmune and inflammatory diseases. The chip has two purposes: fine mapping of 289 established associations corresponding to 187 distinct loci, and deep replication of suggestive, but not yet proven, associations. Fine-mapping regions were defined as 0.2cM centered on GWAS hit SNPs, and all SNPs and short indels in these regions from the 1000 Genomes Project low coverage pilot CEU samples, as well as variants discovered in resequencing experiments conducted by groups collaborating in the chip design were selected for inclusion. Replication of autoimmune and inflammatory GWAS (including Crohn's disease and ulcerative colitis) contributed the bulk of the remaining SNP lists. Approximately 25,000 SNPs were included as replication of unrelated diseases as part of the WTCCC2 project, which serve as useful null SNPs for these analyses. In total, approximately 240,000 SNPs were selected for inclusion, with an assay design success rate of ~80%. Reprinted from Jostins et al. (2012), with permission from Springer Nature (see https://www.nature.com/reprints/permission-requests.html). Jostins, L., et al. (2012) Host-microbe interactions have shaped the genetic architecture of inflammatory bowel disease. Nature, 491, 119-24. PMID: 23128233.

  Study phs001721

• Germline hypomorphic CARD11 mutations in severe atopic disease

  We described five patients from four unrelated families who had severe atopic dermatitis and susceptibility to mild skin viral infections. Some patients also exhibited other allergic symptoms such as asthma and food allergy. We performed targeted sequencing for affected patients A-I of family A and B-I of family B, whole exome sequencing for D-I, D-II.1 and healthy mother D-II.3 of family D, and whole genome sequencing for C-I and healthy mother C-II.2 of family C. We identified three different novel damaging missense mutations of CARD11 (GeneID:84433) in patients A-I (L194P), B-I (R975W), and C-I (E57D). In patients D-I and D-II.1, a dup183_196 (p.K196_V197insMKEERDSYNDELVK) mutation was found in the coil-coil domain of the protein. These mutations act as dominant negatives impeding the scaffolding function of the wild type protein, which leads to downstream NF-kB and mTORC1 signaling defects after TCR stimulation. Sanger sequencing later identified the corresponding CARD11 mutations in the affected family members (C-II.1, D-II.2) that are not included in the dbGaP submission. Hypomorphic CARD11 dominant negative mutants are novel causes of severe atopic disease.

  Study phs001369

• Discovery of Non-ETS Gene Fusions in Human Prostate Cancer using Next Generation RNA Sequencing

  Half of prostate cancers harbor gene fusions between TMPRSS2 and members of the ETS transcription factor family. To date little is known about the presence of non-ETS fusion events in prostate cancer. We employed next-generation transcriptome sequencing (RNA-Seq) in order to explore the whole transcriptome of 25 human prostate cancer samples for the presence of chimeric fusion transcripts. We generated more than 1 billion sequence reads and used a novel computational approach (FusionSeq) in order to identify novel gene fusion candidates with high confidence. In total, we discovered and characterized seven new cancer-specific gene fusions, two involving the ETS genes ETV1 and ERG, and five involving non-ETS genes such as CDKN1A (p21), CD9 and IKBKB (IKK-beta), genes known to exhibit key biological roles in cellular homeostasis or assumed to be critical in tumorigenesis of other tumor entities, as well as the oncogene PIGU and the tumor suppressor gene RSRC2. The novel gene fusions are found to be of low frequency but interestingly, the non-ETS fusions were all present in prostate cancer harboring the TMPRSS2-ERG gene fusion. Future work will focus on determining if the ETS rearrangements in prostate cancer are associated or directly predispose to a rearrangement prone phenotype.

  Study phs000310

• Somatic Mutation Profile by Next Generation Sequencing in HER2+ Breast Cancer

  The HER pathway is the driving force behind 30% of human breast cancers. It is important to understand how targeted therapies block different cellular pathways, and mechanisms of escape from this blockage. Therapies directed at HER2 establish a successful treatment paradigm, but de novo and acquired resistance exist. The HER signaling system is a complex network with four receptors and eleven ligands, a phosphorylation signaling cascade, and many transcription factors, all complicated by both positive and negative feedback circuits. Analysis of genomes, exomes and transcriptomes by next generation sequencing is aimed at uncovering the genetic factors responsible for patient responses to HER2-directed therapies. We are sequencing HER2-overexpressing cancers, in order to identify potential somatic changes that may better select patients who will benefit from therapy, to determine new targets that may overcome resistance, and to improve outcomes with known current HER2-targeted therapies. Whole exome capture sequencing will determine somatic mutation profiles in HER2-overexpressing tumors, to comprehensively characterize the somatic alterations, with the goal of identifying those patients most likely to respond, as well as discovering new targets that may overcome resistance to HER2-directed therapy. We augment the whole exome data with RNA-seq data to determine expression levels of somatic mutations we discover.

  Study phs000770

• Targeted MitoExome Sequencing of Mitochondrial OXPHOS Diseases (Massachusetts General Hospital)

  Human disorders of mitochondrial oxidative phosphorylation (OXPHOS) represent a devastating collection of inherited diseases. These disorders impact at least 1:5000 live births, and are characterized by multi-organ system involvement. They are characterized by remarkable locus heterogeneity, with mutations in the mtDNA as well as in over 77 nuclear genes identified to date. It is estimated that additional genes may be mutated in these disorders. To discover the genetic causes of mitochondrial OXPHOS diseases, we performed targeted, deep sequencing of the entire mitochondrial genome (mtDNA) and the coding exons of over 1000 nuclear genes encoding the mitochondrial proteome. We applied this 'MitoExome' sequencing to 124 unrelated patients with a wide range of OXPHOS disease presentations from the Massachusetts General Hospital Mitochondrial Disorders Clinic. The 2.3Mb targeted region was captured by hybrid selection and Illumina sequenced with paired 76bp reads. The total set of 1605 targeted nuclear genes included 1013 genes with strong evidence of mitochondrial localization from the MitoCarta database, 377 genes with weaker evidence of mitochondrial localization from the MitoP2 database and other sources, and 215 genes known to cause other inborn errors of metabolism. Approximately 88% of targeted bases were well-covered (>20X), with mean 200X coverage per targeted base.

  Study phs000339

• Rapid Evolution of a Skin Lightening Allele in Southern African KhoeSan

  The genetic architecture and polygenicity of skin pigmentation is explored in KhoeSan communities, including Nama individuals whose phenotype and genotype data are provided in this dataset. By pairing genotypes and quantitative spectrophotometric skin pigmentation phenotypes, we show that skin pigmentation is highly heritable in KhoeSan, yet known pigmentation loci only explains a small fraction of the phenotypic variance. Using genome-wide association analyses, we identified both canonical and non-canonical skin pigmentation loci. We show that this phenotype is more polygenic and complex than previously characterized. (Martin et al., Cell 2017. PMID: 29195075) Following up on the top associated signal with large effect size in the gene SLC24A5, we demonstrate that the canonical Eurasian nonsynonymous allele was introduced into KhoeSan via a recent migration ~2 kya and was under extremely strong selection. The derived allele was present at high frequency despite controlling for gene flow. With high-throughput sequences in the captured SLC24A5 region, we show that the most common derived haplotype is identical amongst Europeans, eastern African and KhoeSan. Using 4-population demographic simulations with selection, we show that the allele was introduced into the KhoeSan only 2,000 ya via a back-to-Africa migration and then experienced a selective sweep.

  Study phs001753

• Generation of a Cellular Atlas of Human Adipose Tissue

  White adipose tissue (WAT), once regarded as morphologically and functionally bland, is now recognized to be dynamic, plastic, heterogenous, and involved in a wide array of biological processes including energy homeostasis, glucose and lipid handling, blood pressure control, and host defense.  High fat feeding and other metabolic stressors cause dramatic changes in adipose morphology, physiology, and cellular composition, and alterations in adiposity are associated with insulin resistance, dyslipidemia, and Type 2 diabetes (T2D). Here we provide detailed cellular atlases of human and murine subcutaneous and visceral white fat at single-cell resolution across a range of body weight. We identify subpopulations of adipocytes, adipose stem and progenitor cells (ASPCs), vascular, and immune cells and demonstrate commonalities and differences across species and dietary conditions. We link specific cell types to increased risk of metabolic disease, and we provide an initial blueprint for a comprehensive set of interactions between individual cell types in the adipose niche in leanness and obesity. We also utilize bulk RNA sequencing data of enzymatically digested human adipocytes to scale these findings to a larger cohort. These data comprise an extensive resource for the exploration of genes, traits, and cell types in the function of WAT across species, depots, and nutritional conditions. 

  Study phs002766

• A Chromatin Accessibility Atlas of the Developing Human Telencephalon

  Gene expression differs between cell types and regions within complex tissues such as the developing brain. To discover regulatory elements underlying this specificity, we generated genome-wide maps of chromatin accessibility in nine anatomically-defined regions of the developing human telencephalon. Additionally, we defined the histone modification landscape of the prefrontal cortex and generated chromatin accessibility maps of its upper and deep layers. We predicted a subset of open chromatin regions (18%) that are most likely to be active enhancers, many of which are dynamic with 26% differing between early and late mid-gestation and 28% present in only one brain region. These predicted regulatory elements (pREs) are enriched proximal to genes with expression differences across developmental stages, regions, and cortical laminae; they harbor distinct sequence motifs that suggest potential upstream regulators. We leveraged this atlas to predict and validate novel regulatory elements of genes that control cortex laminar identity and genes associated with autism spectrum disorder (ASD). These include enhancers proximal to FEZF2 and BCL11A that were validated in mouse, and an enhancer of ASD gene SLC6A1 containing two functional de novo mutations in individuals with ASD whose enhancer function we validated via CRISPRa. We also include data for transcription factor binding in mid-fetal cortex tissues, using CUT&RUN methods.

  Study phs002033

• Joint-Specific TF Regulation in RA

  The mechanisms responsible for the distribution and severity of joint involvement in rheumatoid arthritis (RA) are not known. To explore whether site-specific fibroblast-like synoviocyte (FLS) biology might be associated with location-specific synovitis and explain the predilection for hand (wrist/metacarpal phalangeal joints) involvement in RA, we generated bulk, whole-genome transcriptomic and chromatin accessibility data from FLS. Samples were matched in age, sex, and were homogeneous in that all had significant joint damage with longstanding disease requiring surgical intervention. Analysis revealed joint-specific patterns of FLS phenotypes, with proliferative, migratory, proinflammatory, and matrix-degrading characteristics observed in resting FLS derived from the hand joints compared with hip or knee. TNF-stimulation amplified these differences, with greater enrichment of proinflammatory and proliferative genes in hand FLS compared with hip and knee FLS. Hand FLS also had the greatest expression of markers associated with an ‘activated' state relative to the ‘resting' state, with the greatest cytokine and MMP expression in TNF-stimulated hand FLS. Predicted differences in proliferation and migration were biologically validated with hand FLS exhibiting greater migration and cell growth than hip or knee FLS. Distinctive joint-specific FLS biology associated with a more aggressive inflammatory response might contribute to the distribution and severity of joint involvement in RA.

  Study phs003633

• Tumor detection by analysis of both symmetric- and hemi-methylation of plasma cell-free DNA

  Aberrant DNA methylation plays a critical role in tumorigenesis. While DNA methylation has been used for cancer detection and classification, DNA hemi-methylation, a novel epigenetic mark has not been analyzed extensively in cancer epigenomes. Here we report a strand-specific (ss) sequencing method (MeDIP-Seq) for plasma cell free (cf) DNA (sscf-MeDIP-Seq), which can analyze both symmetrically methylated DNA regions (DMRs) as well as hemi-methylated regions (DHMRs). Using the sscf-MeDIP-Seq method, we analyzed plasma cfDNA methylomes of 271 samples from subjects with liver cancer and brain cancer and from individuals without cancer (controls). Among them, 215 samples were chosen as the discovery cohort for the identification of DMRs and DHMRs specific in each subject group and for the training of machine learning models of multi-cancer detection (MCD) using DMRs, DHMRs and DMRs+DHMRs as inputs. These models were then used to predict the 56 samples in the validation cohort. We found that models trained with DMRs+DHMRs as inputs in general outperformed models trained with DMRs or DHMRs alone, with AUROC being 0.971, 0.981, and 0.99 in predicting control, liver and brain cancer samples in the validation cohort, respectively, by the DMR+DHMR-trained models.

  Study phs003462

• Melanoma Genome Sequencing Project

  Melanoma is the most prevalent cause of skin cancer morbidity and mortality. In order to characterize the full range of somatic mutations that may drive the growth of melanoma, we are sequencing tumor and normal DNA from a set of roughly 150 melanomas. For the majority of samples (approximately 90% of the cases), mutations in protein coding genes will be assessed in the exonic DNA of tumor-normal pairs using hybrid capture and paired-end DNA sequencing. For a few of DNA samples (approximately 10% of the cases), the entire genomes will be analyzed to assess the possible contribution of complex structural rearrangements contributing to oncogenesis. The sequencing data are supplemented by copy number profiling on the same tumors using high-density SNP arrays. Integration of these approaches will enable the unbiased and comprehensive characterization of both known and novel recurrent DNA alterations that arise in melanoma.Versions 1-4 contain whole exome and bulk RNA sequencing data for patients with melanoma.Version 5 - We compared somatic genomic profiles from matched pre-treatment and post-resistance tumor biopsies in patients with metastatic melanoma who exhibited heterogeneous responses to immune checkpoint inhibitors (ICIs). This version includes whole exome sequencing (WES) data from 14 patients.

  Study phs000452

• Childhood Cancer Data Initiative (CCDI): Pediatric In Vivo Testing Program –Sarcoma, Kidney, and Liver Cancers

  The goal was to complete whole exome sequencing (WES) and RNA sequencing (RNA-Seq) on patient derived xenograft samples. Cancer in children is rare with approximately 15,700 new cases diagnosed annually in children 21 years or younger in the U.S. Through use of multimodality therapy (surgery, radiation therapy, and aggressive chemotherapy), 70% of patients will be 'cured' of their disease, and 5-year Event-Free Survival (EFS) exceeds 80%. Consequently, the number of patients that can be enrolled in phase I/II clinical trials is small, and most patients will have been extensively treated, hence drug/radiation resistant. Thus, preclinical studies that accurately translate into effective clinical therapy are an essential component of pediatric drug development. Our group has contributed to studies in the Pediatric Preclinical Testing Program (PPTP) and Pediatric Preclinical Testing Consortium (PPTC) that have led to clinical studies through Children's Oncology Group (COG). Of importance, we have developed and characterized over 330 Patient Derived Xenografts (PDX), developed from tumors both at diagnosis and relapse, that can be used to facilitate pediatric drug development as directed by FDA under the Research to Accelerate Cures and Equity for Children Act (RACE for Children Act).

  Study phs003160

• National Cancer Institute (NCI) Early Onset Malignancies Initiative (EOMI): Molecular profiling of Breast, Colon, Kidney, Liver, Multiple Myeloma, and Prostate among Racially and Ethnically Diverse Populations

  The Early Onset Malignancies Initiative (EOMI) is focused on collecting specimens from racially and ethnically diverse population groups with six different types of cancer (breast, prostate, colon, liver, kidney, and multiple myeloma). This initiative is a partnership between the Center for Research Strategy, the Center for Cancer Genomics and the Division of Cancer (DCP). DCP's NCI Community Oncology Research Program (NCORP) is a critical partner. The acquisition of the cancer and matching normal biospecimens in this initiative are from newly diagnosed early onset cancer patients who identify as Non-Hispanic Black, Hispanic, American Indian or Alaska Native, or Non-Hispanic White. The Minority/Underserved NCORPs sites will collect the samples for this project. By carrying out this initiative in the NCORP's national network of investigators, cancer care providers, and research institutions, the EOMI will harness the community-based partnerships to improve participation among these racial/ethnic groups. The NCORP investigators, as participants in NCI's clinical trials network, will annotate the biologic samples with clinical information that will enhance the ability to answer questions about molecular determinants of early disease, treatment response, and/or prognosis. The EOMI will be a step forward in addressing racial/ethnic disparities and in understanding cancer in underserved populations.

  Study phs001952

• Single-Cell Genomic and Transcriptomic Analysis of the Aging Human Brain

  We studied human brain aging using single-cell whole genome sequencing and single-nucleus RNA sequencing of the pre-frontal cortex from 19 donors ranging in age from 0.4 years to 104 years. All donors were free of neuropathology and none received a clinical neurological diagnosis prior to death. snRNA-seq libraries were prepared in batches that always included old and young donors to mitigate batch effects in the preparations. Neurons for single-cell WGS were collected from the same tissue region as the snRNA-seq for each donor but represent independent cells. Multiplexed error-robust fluorescent in-situ hybridization (MERFISH) data was also collected from pre-frontal cortex samples of a subset of the 19 donors. Transcriptomic results identified infant-specific cell clusters, and common-down regulation of cell-essential genes involved in ribosomes, transport, homeostasis and metabolism during aging, across all cell types. scWGS identified two-age associated mutational signatures, as well as gene-length and expression-level dependent somatic mutation rates in neurons. Raw sequencing data in fastq format for single-nucleus RNA-seq and single-cell DNA-seq, as well as cell and gene count matrix, transcript metadata, and cell metadata for MERFISH will be available in dbGaP.

  Study phs003445

• Examination of Engineered LINE-1 Integration Events in HeLa Cells

  Long Interspsesed Element-1 (LINE1 or L1) encodes two proteins (ORF1p and ORF2p) required for its mobilization (i.e. retrotransposition) in the human genome. Human ORF2p, a 150kDa protein, has endonuclease (EN) and reverse transcriptase (RT) activities. L1 EN is classified as apurinic/apyrimidin-like endonuclease (APE) and is responsible for initiating target-site primer reverse transcription (TPRT) during the L1 retrotransposition cycle. During canonical TPRT, the L1 EN activity makes a single-strand endonucleolytic nick at a double-strand DNA target sequence in genomic DNA. The endonucleolytic nick is typically made at an L1 EN degnerate consensus cleavage sequence (5'-TTTT/A-3') exposing a free 3' hydroxy group that can be used as a primer by the ORF2p L1 RT activity to reverse transcribe the associated L1 RNA. Here, we characterized over 21,000 de novo engineered L1 retrotransposition events in HeLa cells to determine if L1 exhibits overt integration preferences in the human genome. Bru-seq of HeLa cells to determine actively transcribed regions of the genome was also performed. This Bru-seq data was used to determine if engineered L1 integration events preferentially integrate into transcribed regions of the genome, and to determine which strand during transcription may present itself as preferential during L1 retrotransposition.

  Study phs001669

• Transcriptome sequencing of human colon organoid after co-cultivation with Bifidobacterium longum

  The large intestine has a dense milieu of indigenous bacteria, generating a complex ecosystem with crosstalk between individual bacteria and host cells. In vitro host cell modeling and bacterial interactions at the anaerobic interphase have elucidated the crosstalk molecular basis. Although classical cell lines derived from patients with colorectal cancer including Caco-2 cells are used, whether they adequately mimic normal colonic epithelial physiology is unclear. To address this, we performed transcriptome profiling of Caco-2 and Monolayer cells derived from healthy Human Colonic Organoid (MHCO) cultured hemi-anaerobically. Coculture with the anaerobic gut bacteria, Bifidobacterium longum subsp. longum differentiated the probiotic effects of test cells from those of physiologically normal intestinal and colorectal cancer cells. We cataloged non- or overlapping gene signatures where gene profiles of Caco-2 cells represented absorptive cells in the small intestinal epithelium, and MHCO cells showed complete colonic epithelium signature, including stem/progenitor, goblet, and enteroendocrine cells colonocytes. Characteristic gene expression changes related to lipid metabolism, inflammation, and cell-cell adhesion were observed in cocultured live Bifidobacterium longum and Caco-2 or MHCO cells. B. longum-stimulated MHCO cells exhibited barrier-enhancing characteristics, as demonstrated in clinical trials. Our data represent a valuable resource for understanding gut microbe and host cell communication.

  Study JGAS000740

• Bruno et al.: Interferon gamma rebalances immunopathological signatures in Chronic Granulomatous Disease through metabolic rewiring

  Chronic granulomatous disease (CGD) is a primary immunodeficiency characterized by recurrent life-threatening infections and hyperinflammatory complications. It is caused by mutations in the NADPH oxidase complex and the consequent loss of reactive oxygen species (ROS) production. Recombinant human interferon gamma (rIFN-γ) prophylaxis reduces the risk of severe infections, but the mechanisms behind its efficacy in CGD are still an open question, as it does not restore NADPH oxidase-dependent ROS production. Here, we show that myeloid innate immune cells of CGD patients are transcriptionally and functionally reprogrammed to a hyperactive inflammatory status, displaying an impaired in vitro induction of trained immunity. CGD monocytes have reduced intracellular amino acids concentrations and profound functional metabolic defects, both at the level of glycolysis and mitochondrial respiration. Ex vivo and in vivo treatment with IFN-γ restored these metabolic defects and reduced excessive IL-1β and IL-6 production in response to fungal stimuli in CGD monocytes. These data suggest that prophylactic rIFNγ modulates the metabolic status of innate immune cells in CGD. These data shed light on the effects of NADPH-derived ROS deficiency to the metabolic programs of immune cells and pose the basis for targeting this immunometabolic axis, potentially beyond CGD, with IFN-γ immunotherapy.

  Study EGAS00001005463

• Comprehensive gene analysis of colorectal cancer cases

  Colorectal cancer is the third most common morbidity and the fourth leading cause of cancer death in the world.To characterize the mutation profiles of colorectal cancer with liver metastasis,we performed both whole-exosome sequencing and RNA sequencing for 12 cases of colorectal cancer with liver metastasis (group A) and 16 cases of colorectal cancer without liver metastasis (group B). As a result of driver mutation search by Mutsig pipeline, PTEN was identified as a significant mutant gene in group A alone.ADAMTS10,NELL1 and RXFP3 were identified as candidate genes for metastasis by MuTect,SID and VarScan pipeline.Furthermore, in all samples, the copy number of chromosomes 1p,4,7,8,13,15, 17p,18,20 and 22q were significantly abnormal. Chromosomes 1p,4,15,18,17p,22q were deleted and chromosomes 7, 8,13 and 20 were amplified by CNV.Pathway analysed revealed Wnt/Ca2+ and PCP pathway was involved in metastasis of cl cancer.9 kinds of fusion genes were identified in deFuse pipeline. Among them, ADAP1 fusion genes were identified in 2 samples, and ADAP1 fusion genes were further identified in 3 samples by our in-house algorithm. It is suggested that new fusion genes present in such multiple samples may be involved in the development of colorectal cancer.

  Study JGAS000128

• Rapid detection of expanded short tandem repeats in personal genomics using hybrid sequencing.

  Long expansions of short tandem repeats (STRs); i.e., DNA repeats of 2-6 nucleotides, are associated with some genetic diseases. Cost-efficient high throughput sequencing can quickly produce billions of short reads that would be useful for uncovering disease-associated STRs. However, enumerating STRs in short reads remains largely unexplored because of the difficulty in elucidating STRs much longer than 100 bp, the typical length of short reads. We propose ab initio procedures for sensing and locating long STRs promptly by utilizing the frequency distribution of all STRs and paired-end read information. We validated the reproducibility of this method using biological replicates and used it to locate an STR associated with a brain disease (SCA31). Subsequently, we sequenced this STR site in eleven SCA31 samples using SMRT? sequencing (Pacific Biosciences), determined 2.3-3.1 kb sequences at nucleotide resolution, and revealed that (TGGAA)- and (TAAAATAGAA)-repeat expansions determined the instability of the repeat expansions associated with SCA31. Our method could also identify common STRs, (AAAG)- and (AAAAG)-repeat expansions, which are remarkably expanded at four positions in an SCA31 sample. This is the first proposed method for rapidly finding disease-associated long STRs in personal genomes using hybrid sequencing of short and long reads.

  Study JGAS000002

• BHD-associated kidney cancer

  Background Birt-Hogg-Dub? (BHD) syndrome is a hereditary kidney cancer syndrome, which predisposes patients to develop kidney cancer, cutaneous fibrofolliculomas and pulmonary cysts. The responsible gene, FLCN is a tumor suppressor for kidney cancer which plays an important role for energy homeostasis through the regulation of mitochondrial oxidative metabolism. However, how FLCN-deficiency leads to renal tumorigenesis is unclear. Results Whole-exome sequencing analysis of BHD-associated kidney cancer revealed that copy number variations (CNV) of BHD-associated kidney cancer are considerably different from those already reported in sporadic cases. Somatic mutation analysis revealed BHD-associated kidney cancer has quite few common mutations in specific genes; however, chromatin remodeling genes in total were frequently mutated in BHD-associated kidney cancer (17/29 tumors, 59%). Metabolite analysis revealed that glycolysis, pentose phosphate pathway and lipid metabolism were upregulated whereas the rate of beta oxydation and TCA cycle is decreased in BHD-associated kidney cancer compared to normal adjacent kidney. Conclusions BHD-associated kidney cancer displays a unique molecular characteristics which is totally different from that of sporadic kidney cancer, providing the mechanistic insight into tumorigenesis under FLCN deficiency as well as the foundation for the development of novel therapeutics for kidney cancer.

  Study JGAS000115

• Genomic and Genetic Analysis of Brain Tumors and Analysis of Their Clinicopathological Significance

  The treatment outcome of glioma patients remains disappointing, and more studies are needed to improve therapeutic strategies. In particular, it is crucial to understand the mechanisms underlying resistance to treatment and malignant transformation, as well as the intratumoral heterogeneity and interaction of tumor cells and host immunity. Recent advances in omics analysis techniques, such as next-generation sequencing (NGS), have enabled us to obtain omics data more easily and accelerated our ability to integrate this information and precisely characterize cases. In our project, we have already collected more than 100 pairs of human glioma/normal tissue samples, of which some are multiple specimens collected from initial and recurrent tumors, and some include multiple specimens collected from different regions within a single tumor. We used our in-house algorithms, composed of alignment and variant calling programs, to analyze omics data from multiple platforms, such as whole-exome sequencing (WES), RNA sequencing (RNA-seq), and methylation array, with our original cohort of glioma patients. Using several informatics tools, the information about genetic and epigenetic status and expression signatures were analyzed. Those data are expected to provide key information for investigating the mechanisms discussed above. Through this project, we hope to improve strategies for personalized therapy based on comprehensive omics data.

  Study JGAS000004

• The impact of CCR8+ regulatory T cells on cytotoxic T cell function in human lung cancer

  Regulatory T cells (Tregs) suppress the host immune response and maintain immune homeostasis. Tregs also promote cancer progression and are involved in resistance to immune checkpoint inhibitor treatments. Recent studies identified selective CCR8 expression on tumor-infiltrating Tregs; CCR8+ Tregs have been indicated as a possible new target of cancer immunotherapy. Here, we investigated the features of CCR8+ Tregs in lung cancer patients. CCR8+ Tregs were highly activated and infiltration of CCR8+ Tregs in tumors was associated with poor prognosis in lung cancer patients. We also investigated their immune suppressive function, especially the influence on cytotoxic T lymphocyte cell function. The Cancer Genome Analysis revealed that CD8 T cell activities were suppressed in high CCR8-expressing tumors. Additionally, depletion of CCR8+ cells enhanced CD8 T cell function in an ex vivo culture of lung tumor-infiltrating cells. Moreover, CCR8+ Tregs, but not CCR8��� Tregs, induced from human PBMCs markedly suppressed CD8 T cell cytotoxicity. Finally, we demonstrated the therapeutic effect of targeting CCR8 in a murine model of lung cancer. These findings reveal the significance of CCR8+ Tregs for immunosuppression in lung cancer, especially via cytotoxic T lymphocyte cell suppression, and suggest the potential value of CCR8-targeted therapy for cancer treatment.

  Study JGAS000454

• Microbiomic and immunogenic biomarkers of adjuvant chemotherapy efficacy in stage III colorectal cancer

  The standard treatment for stage III colorectal cancer (CRC) is surgery, followed by adjuvant chemotherapy (AC). However, the efficacy of AC is limited, and approximately 30% of patients experience recurrence. Therefore, easily assessable molecular biomarkers that can predict response to AC are needed. Here, we aimed to identify predictive biomarkers of response to AC using tumor transcriptomic data. Tumor specimens from 253 patients who underwent surgery for stage III CRC between 1997 and 2019 were analyzed using RNA sequencing. Transcriptional subtyping and analyses of gut microbiome, immune cell fractions, and T-cell receptor (TCR) were performed to identify factors associated with recurrence-free survival (RFS) in patients treated with or without AC. Of 253 patients, 118 (46.6%) received AC. Twenty phyla and 799 species of the microbiomes were identified, among which Fusobacteria were identified in 247 patients (97.6%). Multivariable Cox proportional hazards regression analysis of patients who received AC showed that intratumoral Fusobacterium and M1 macrophage abundance and diversity of TCR alpha and delta (TRAD) were independently associated with RFS. Fusobacterium and M1 macrophage abundance, and TRAD diversity were associated with the response to AC, suggesting that these factors could serve as biomarkers for personalized treatment in stage III CRC.

  Study JGAS000875

• Shedding light over COVID-19 susceptibility and severity

  “One of the most remarkable features of SARS-CoV-2 infection is the variation in consequence ranging from asymptomatic to life-threatening” state the authors of the majestic Nature paper released in accelerated pre-view on July 2021. This variation in consequences represent not only an issue of burning medical relevance, but also a scientific and intellectual challenge for scientists around the globe. Researchers are known to be easily triggered by a challenge, and they recently demonstrated to be even better at double their efforts to focus on an urgent societal need. In this work led and coordinated by the COVID-19 Host Genetics Initiative (COVID-19 HGI), researchers gathered data from 49,562 cases and two million controls from 46 studies performed worldwide. The majority of the data has been deposited at the EGA, where 29 datasets have been collected under the study number EGAS00001005304. Data gathered from diverse sources conferred a huge robustness to the results. Genetic data from different ancestry origin were included in the analyses, even though data from patients with European origin was highly predominant (approx. 75%). Separated analysis of the different ancestries gave rise to relevant differences in the results, highlighting how important it is to take this aspect in consideration for future studies. Data from the individual studies was compiled into three categories according to COVID-19 severity for a meta-analysis: (a) critically ill cases (b) moderate or severe cases (c) all cases. They were paired with non-COVID-19-reported controls. All three analysis resulted in a total of 13 genetic loci associated with COVID-19 severity. Out of those, 9 loci associated specifically with critical illness, while 6 (a few overlapped) associated with moderated illness. Some genomic position matched previous findings. Interestingly, working with such an enormous amount of data allowed the depth to discriminate association with disease susceptibility vs severity, an issue of large medical relevance. Several loci mapped into biological relevant genes, like TYK2, PPP1R15A and FOXP4, previously implicated in immune response and/or lung functions. Researchers participating in the COVID-19 Host Genetics Initiative have worked in an unprecedented time-tight collaboration, prioritizing data sharing and collective results to tackle the pandemic affecting the world. The scientific community won’t miss any of the lessons to learn from this, and the EGA is ready to keep working to enable safe data sharing and fruitful collaborations.

  Blog covid-19-susceptibility-and-severity

• Recurrent somatic JAK-STAT mutations within a novel RUNX1-mutated pedigree

  The acquisition of somatic mutations is an emerging field of investigation in familial leukemia. Currently, genetic profiles in familial MDS/AML are considered analogous to sporadic disease, although the patterns of clonal evolution within families are poorly defined. We performed whole exome profiling of tumour samples from a novel RUNX1 mutated family, to determine the stepwise evolution of MDS/AML across 4 young siblings. Three siblings developed monocytic AML/RAEB2 at 5 years of age, with hepatosplenomegaly and somatic mutations upregulating JAK-STAT signalling, the latter are typically detected in <5% of sporadic MDS/AML. Two siblings acquired the canonical JAK2 V617F mutation, while another acquired a unique missense mutation of SH2B3, a negative regulator of JAK2. Notably, 2/3 siblings demonstrated dosage amplification of these mutations due to acquired uniparental disomy of chromosomes 9p and 12q (encompassing JAK2 and SH2B3, respectively). All 4 siblings were heterozygous for the 46/1 JAK2 haplotype associated with predisposition to sporadic V617F myeloproliferative disorders, which likely influenced the acquisition of JAK2 mutations. Our findings provide further evidence that relatives with shared germline mutations may acquire somatic mutations in a non-random manner leading to convergent patterns of disease evolution.

  Study EGAS00001001862

• Comparison of EGF and PDGF driven glioblastomas.

  Glioblastoma multiforme (GBM) is the most common and aggressive type of primary brain tumor. Epidermal growth-factor (EGF) and platelet-derived growth-factor (PDGF) receptors are frequently amplified and/or possess gain-of-function mutations in these tumors. Despite years of research into tyrosine-kinase inhibitors, clinical trials have shown disappointing efficacy. One of the underlying reasons is intra-tumor heterogeneity. To assess the effect of intra-tumor heterogeneity on gene expression, we derived a general approach to map single-cell expression profiles to sequentially-acquired mutations identified from exome sequencing. Using 288 single cells, we constructed high-resolution phylogenies of EGF-driven and PDGF-driven GBMs, modeling transcriptional kinetics during tumor evolution. A deletion in a PDGF-receptor domain involved in dimerization correlated to an up-regulation of the PI3K/AKT pathway in a proneural GBM, and ectopic expression of this mutant PDGFRA in vitro enhanced proliferation and invasiveness. Descending the phylogenetic tree of the PDGFR-driven tumor corresponded to a progressive induction of an oligodendrocyte progenitor-like cell type that expresses growth and angiogenic factors. In contrast, phylogenetic analysis of an EGFR-amplified tumor showed a gradual up-regulation of pro-invasive genes.

  Study EGAS00001001900

• The Genomic Landscape of Childhood and Adolescent Melanoma

  Despite remarkable advances in the genomic characterization of adult melanoma, the molecular pathogenesis of pediatric melanoma remains largely unknown. We analyzed 15 conventional melanomas (CM), 3 melanomas arising in congenital nevi (CNM), and 5 spitzoid melanomas (SM), using various platforms, including whole genome or exome sequencing, molecular inversion probe assay, and/or targeted sequencing. CM demonstrated a high burden of somatic single nucleotide variations (SNV), with each case containing a TERT promoter (TERT-p) mutation, 13/15 containing an activating BRAF V600 mutation, and >80% of the identified SNV consistent with UV damage. By contrast, the 3 CNM contained an activating NRAS Q61 mutation and no TERT-p mutations. SM were characterized by chromosomal rearrangements resulting in activated kinase signaling in 40%, and an absence of TERT-p mutations except for the one SM that succumbed to hematogenous metastasis. We conclude that pediatric CM has a very similar UV-induced mutational spectrum to that found in the adult counterpart, emphasizing the need to promote sun protection practices in early life and to improve access to therapeutic agents being explored in adults in young patients. By contrast, the pathogenesis of CNM appears to be distinct. TERT-p mutations may identify the rare subset of spitzoid melanocytic lesions prone to disseminate.

  Study EGAS00001000901

• LySeqST: A targeted sequencing assay for robust genomic classification of diffuse large B-cell lymphoma

  Diffuse large B-cell lymphoma (DLBCL) is a biologically heterogeneous disease. Two genomic classification systems, LymphGen and DLBclass, are capable of reproducibly classifying single tumors into subtypes with prognostic relevance in the context of standard-of-care R-CHOP chemoimmunotherapy. Preliminary data from subgroup analyses of clinical trials suggest that the distinct biological features are targetable with drugs that exploit subtype-specific vulnerabilities, and large clinical trials are being designed to test this hypothesis. Although whole exome sequencing (WES) remains the gold standard for each of these classification systems, many laboratories may opt for smaller targeted sequencing panels that reduce costs associated with sequencing, data storage, and processing. Here, we describe the design and validation of a targeted sequencing panel (LySeqST) that captures the genomic features required specifically for LymphGen classification. We perform in silico and real-world validation vs. WES data and demonstrate that LySeqST can be used to accurately classify DLBCL tumors. In an unselected population-based cohort, we determine the real-world proportions and validate the prognostic relevance of LymphGen subtypes, confirming that only tumors expressing the dark zone gene expression signature are consistently associated with inferior outcomes. Our results support the use of LySeqST for accurate genomic classification of DLBCL.

  Study EGAS50000001601

• The evolutionary steps from primary to metastatic prostate cancer are largely uncharted, and the ability to use DNA present in body fluids as correlates of aggregate metastatic status is under-examined. We reconstructed phylogenies in ten prostate cancer patients with fatal disease using deep targeted sequencing of the prostate, adjacent and distant organs, as well as plasma, serum, and cerebrospinal fluid at various time points. A total of 163 samples are studied.

  The evolutionary steps from primary tumor to metastasis in prostate cancer are largely uncharted, and the ability to use serum, plasma, and cerebrospinal fluid as a correlate of aggregate metastatic tumor genomic status has not been tested. We used deep targeted sequencing to reconstruct tumor evolution in ten prostate cancer patients with fatal disease encompassing examination of the prostate and adjacent and distant organs, as well as plasma, serum, and cerebrospinal fluid at various time points. We show that there is substantial evolution from a common ancestor within the prostate that results in branching to multiple lineages which form an intermixed multi-clonal primary tumor mass. After the occurrence of key driver aberrations, one of these lineages will metastasize to multiple sites in a sequential fashion. These metastatic sites are then susceptible to being populated by cells from other intra-prostatic lineages or from other metastases. Genomic representation of metastases in body fluids is not uniform. Cerebrospinal fluid analysis can detect lineages not detected in circulating DNA, suggesting possible clinical utility.

  Study EGAS00001003848

• Pancreatic, Small-intestinal and Pulmonary Neuroendocrine Tumors

  Purpose: The primary origin of neuroendocrine tumor metastases can be difficult to determine by histopathology alone, but is critical for therapeutic decision making. DNA methylation-based profiling is now routinely used in the diagnostic workup of brain tumors. This has been enabled by the availability of cost-efficient array-based platforms. We have extended these efforts to augment histopathological diagnosis in neuroendocrine tumors. Experimental Design and Results: We compiled data of 69 small-intestinal, pulmonary, and pancreatic neuroendocrine tumors. These data were used to build a ridge regression calibrated random forest classification algorithm (NEN-ID) that predicts the origin of tumor samples with high accuracy (> 95%). The model was validated during 3x3 nested cross validation and tested in a local (n=26) and external (n=172) cohort. In addition, we show that our diagnostic approach is robust across a range of possible confounding experimental parameters such as tumor purity and array quality. A software infrastructure and online user interface was built to make the model available to the scientific community. Conclusions: This DNA methylation-based prediction model can be used in the workup for patients with neuroendocrine tumors of unknown primary. To facilitate validation and clinical implementation, we provide a user-friendly, publicly available web-based version of NEN-ID.

  Study EGAS00001004878

• Emergence of oncofetal plasticity is ubiquitous in early colorectal cancers

  Metastasis formation is classically considered a late-stage event in colorectal cancer (CRC) evolution. Yet, the time and patterning by which metastatic competence is acquired remain poorly understood. Here, we show that metastasis-associated oncofetal cell states already emerge at the earliest stages of CRC, concurrent with invasive front formation. However, while necessary for metastasis, we detect them ubiquitously among early non-metastatic cancers, highlighting additional bottlenecks like immune evasion. To understand how oncofetal cells first emerge, we generated multiregional organoid models that reflect successive tumor progression stages within individual early-stage CRCs. Whole-genome sequencing and growth factor-dependency assays exclude tumor cell-intrinsic acquired traits. In contrast, single-cell spatial atlases of the tumor-microenvironment before and after malignant transformation revealed stereotypic patterning of fibroblast subtypes resembling normal tissue architecture, resulting in distinct regional microenvironments. At the onset of malignant growth into the submucosa, the first cancer-associated fibroblasts (CAFs) to appear strongly resemble submucosal trophocytes and colocalize with oncofetal cell states at invasive fronts. Functionally, fibroblast-organoid co-cultures confirm that these trophocyte-like CAFs induce plastic transitioning to oncofetal states. Thus, interactions between tumor and submucosal fibroblasts directly following malignant transformation dictate the timing and location at which oncofetal plasticity first occurs during CRC progression.

  Study EGAS50000001532

• Sclerosing epithelioid fibrosarcoma case report

  Sclerosing epithelioid fibrosarcoma (SEF) is a rare and aggressive soft tissue sarcomathought to originate in fibroblasts of the tissues surrounding tendons, ligaments and muscles. Minimally responsive to conventional cytotoxic chemotherapies, greater than 50% of SEF patients experience relapse and/or metastatic disease. SEF is most commonly discovered in middle-aged and elderly adults, but also rarely in children. A common gene fusion occurring between the EWSR1 and CREB3L1 genes has been observed in 80-90% of SEF cases. We describe here the youngest SEF patient reported to date (a 3-year-old Caucasian male) who presented with numerous bony and lung metastases. Additionally, we perform a comprehensive literature review of all SEF-related articles published since the disease was first characterized. The patient described in this report experienced persistent disease progression despite being heavily treated with multiple surgeries, radiation, numerous chemotherapies and targeted therapeutics. Primary tumor, metastatic and relapse sites were sequenced by whole genome, whole exome, and deep transcriptome next generation sequencing with comparison to a patient-matched normal blood sample. Consistent across all sequencing analyses was the disease-defining EWSR1-CREB3L1 fusion as a single feature consensus. We provide an analysis of our genomic findings and a discussion as to potential therapeutic strategies for SEF.

  Study EGAS00001005214

• Hematopoietic Tet2 inactivation enhances the response to checkpoint blockade immunotherapy

  Somatic mutations inactivating TET2 are among the most common drivers of clonal hematopoiesis (CH). While TET2 inactivation is associated with monocyte-derived inflammation and improved chimeric antigen-receptor-T cell function, its impact on immunotherapy response is unknown. In murine models, hematopoietic Tet2 mutation enhanced the immune checkpoint blockade (ICB) response via the combined presence of phagocytes, CD4 and CD8 T cells. This effect was lost with myeloid- or T-cell restricted Tet2 inactivation, or in mice with 20% Tet2-mutant hematopoiesis. Mechanistically, in Tet2-mutant tumor-infiltrating leukocytes (TILs), ICB preferentially restricted cell states linked to tumor progression while inducing anti-tumor states. Tet2-mutant monocytes activated costimulatory programs, while Tet2-mutant T cells showed enhanced T cell memory signatures, alongside decreased exhaustion and regulatory phenotypes. This murine data was clinically relevant, since tumors from colorectal cancer and melanoma patients with TET2 driver mutation-CH (TET2-CH) showed enhanced immune infiltration, inflammation, and T cell activation. In melanoma patients treated with ICB, TET2-CH was associated with 6-fold greater odds of clinical benefit. Collectively, we demonstrate that hematopoietic Tet2 inactivation primes leukocytes for anti-tumor states associated with immunotherapy response and provides a potential biomarker for personalized therapy.

  Study EGAS50000001191

• IDENTIFICATION AND TARGETED MANAGEMENT OF A PATIENT WITH A NEURODEGENERATIVE DISORDER CAUSED BY BIALLELIC MUTATIONS IN SLC5A6

  We describe a sibling pair displaying an early infantile-onset, progressive neurodegenerative phenotype, with symptoms of developmental delay and epileptic encephalopathy developing from 12 to 14 months of age. Using whole exome sequencing, compound heterozygous variants were identified in SLC5A6, which encodes the sodium-dependent multivitamin transporter (SMVT) protein. SMVT is an important transporter of the B-group vitamins biotin, pantothenate and lipoate. The protein is ubiquitously expressed and has major roles in vitamin uptake in the digestive system, as well as transport of these vitamins across the blood-brain barrier. Pathogenicity of the identified variants was demonstrated by impaired biotin uptake of mutant SMVT. Identification of this vitamin transporter as the genetic basis of this disorder guided targeted therapeutic intervention, resulting clinically in improvement of the patient’s neurocognitive and neuromotor function. This is the second report of biallelic mutations in SLC5A6 leading to a neurodegenerative disorder due to impaired biotin, pantothenate and lipoate uptake. The genetic and phenotypic overlap of these cases confirms mutations in SLC5A6 as the genetic cause of this disease phenotype. Recognition of the genetic disorder caused by SLC5A6 mutations is essential for early diagnosis and to facilitate timely intervention by triple vitamin (biotin, pantothenate and lipoate) replacement therapy.

  Study EGAS00001003861

• Patient-derived pediatric brain tumor organoids faithfully recapitulate primary tumors

  Extensive molecular analysis by many groups has revealed a heterogeneous landscape of embryonal brain tumors, including medulloblastoma (MB) and atypical teratoid/rhabdoid tumor (ATRT). MB, the most common malignant pediatric brain tumor, is comprised of four molecularly distinct subgroups: Wingless (WNT), Sonic Hedgehog (SHH), Group 3 (G3), and Group 4 (G4). Among these, MB-SHH harboring MYCN amplification and TP53 mutation and MB-G3 with MYC overexpression or amplification constitute the most aggressive tumors, have the poorest clinical outcomes, and lack effective therapeutic options for recurrent disease. Here we describe the establishment of tumor organoid (TO) models from patient-derived orthotopic xenografts (PDOX) of MYCN-amplified and TP53-mutant MB-SHH, MYC-amplified MB-G3 and two ATRT which all grew in mice when reimplanted into the cortices of naïve recipient animals (TOX). Methylation, bulk and single-cell RNA sequencing, and whole genome sequencing demonstrated that TOs and TOXs faithfully recapitulated the molecular landscape and heterogeneity of PDOXs and primary patient tumors compared to established cell lines. The development of TOs will facilitate in vitro high-throughput drug or CRISPR screens without the need for use tumor-bearing mice. This will also accelerate the identification and validation of new vulnerabilities and in vivo preclinical testing toward clinical trials.

  Study EGAS00001008305

• Characterizing microbiome-directed fibre snacks in gnotobiotic mice and humans

  Knowledge of the interrelationships between what we eat and the configurations of our gut microbial communities is providing important insights into how food components that are not directly metabolized by human enzymes are linked to our physiology and health status. Changing food preferences brought about by Westernization that have deleterious health effects1,2, plus rapid population expansion, ongoing challenges to sustainable agriculture, and other forces contributing to increased food insecurity, are catalyzing efforts to identify more nutritious and affordable foods3. The gut microbial community is complex, dynamic, and exhibits considerable intra- and interpersonal variation in its composition and functions. The massive number of potential interactions between its components makes it challenging to define the mechanisms by which food ingredients affect community properties. There is also a paucity of information about the ‘bioactive’ ingredients of foods that influence the fitness and expressed functions of community members. Here, plant fibres, from different sustainable sources and targeting distinct features of obese human gut microbiomes in gnotobiotic mice, were formulated into snack prototypes and used to supplement controlled diets consumed by overweight and obese adults; the results revealed fibre-specific changes in their microbiomes that were linked to changes in their plasma proteomes indicative of altered physiologic state.

  Study EGAS00001005268

• Tumor-derived exosomes modulate PD-L1 expression in monocytes

  In chronic lymphocytic leukemia (CLL), monocytes and macrophages are skewed toward protumorigenic phenotypes, including the release of tumor-supportive cytokines and the expression of immunosuppressive molecules such as programmed cell death 1 ligand 1 (PD-L1). To understand the mechanism driving protumorigenic skewing in CLL, we evaluated the role of tumor cell–derived exosomes in the cross-talk with monocytes. We carried out RNA sequencing and proteome analyses of CLL-derived exosomes and identified noncoding Y RNA hY4 as a highly abundant RNA species that is enriched in exosomes from plasma of CLL patients compared with healthy donor samples. Transfer of CLL-derived exosomes or hY4 alone to monocytes resulted in key CLL-associated phenotypes, including the release of cytokines, such as C-C motif chemokine ligand 2 (CCL2), CCL4, and interleukin-6, and the expression of PD-L1. These responses were abolished in Toll-like receptor 7 (TLR7)–deficient monocytes, suggesting exosomal hY4 as a driver of TLR7 signaling. Pharmacologic inhibition of endosomal TLR signaling resulted in a substantially reduced activation of monocytes in vitro and attenuated CLL development in vivo. Our results indicate that exosome-mediated transfer of noncoding RNAs to monocytes contributes to cancer-related inflammation and concurrent immune escape via PD-L1 expression.

  Study EGAS00001002377

• The Spatial Heterogeneity in Multiple Myeloma - from the Subclonal Architecture to the Immune Microenvironment (partly hipo_K08K)

  Multiple myeloma (MM) cells show pronounced heterogeneity not only within a given patient but also between spatially separated tumor loci in the bone marrow. Understanding this spatial heterogeneity is emerging as a critical challenge for the successful treatment of the disease. Yet, our understanding of spatial differences in the subclonal architecture, molecular signatures and interactions with the tumor microenvironment remains very limited. To address this shortcoming, we performed bulk and single-cell multi-region sequencing, including random bone marrow samples from the iliac crest and paired imaging-guided focal lesion specimens from 15 newly diagnosed MM patients. We found a median of 5 subclones per patient and unique subclones in focal lesions. Central features of spatial heterogeneity included a consistent down-regulation of the chemoattractant cytokines CXCL7 and CXCL12 in focal lesion tumor cells as well as a significant depletion of macrophages in the focal lesion microenvironment. In contrast, a site-specific expansion of T-cell clones was not detected in focal lesions. In conclusion, our results demonstrate the relevance of considering spatial heterogeneity with potential implications for immunotherapies and molecular studies analyzing the role of gene signatures and MM-microenvironment interactions in disease progression.

  Study EGAS00001006090

• UK10K NEURO ABERDEEN

  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 sample set comprises cases of schizophrenia with additional cognitive measurements, collected in Aberdeen, Scotland.For further information on the Aberdeen cohort please contact David St Clair (d.stclair@abdn.ac.uk).

  Study EGAS00001000109

• Ethiopia_Genome_Project__high_coverage_

  The Ethiopian area stands among the most ancient ones ever occupied by human populations and their ancestors. Particularly, according to archaeological evidences, it is possible to trace back the presence of Hominids up to at least 3 million years ago. Furthermore, the present day human populations show a great cultural, linguistic and historic diversity which makes them essential candidate to investigate a considerable part of the African variability. Following the typing of 300 Ethiopian samples on Illumina Omni 1M (see Human Variability in Ethiopia project, previously approved by the Genotyping committee) we now have a clearer idea on which populations living in the area include the most of the diversity.This project therefore aims to sequence the whole genome of 300 individuals at low (4-8x) depth belonging to the six most representative populations of the Ethiopian area to produce a unique catalogue of variants peculiar of the North East Africa. Furthermore 6 samples (one from each population) will also be sequenced at high (30x) depth to ensure full coverage of the diversity spectrum.The retrieved variants will be of great help in evaluating the demographic dynamics of those populations as well as shedding light on the migrations out of Africa.

  Study EGAS00001000237

• Exome_trios_in_patients_with_gastroschisis

  Gastroschisis (MIM 230750) is a herniation of the intestines through a defect of the abdominal wall lateral to the umbilicus (usually on the right side), and it is not covered by a membrane [Ledbetter, 2012]. Gastroschisis is a congenital anomaly with increasing incidence, easy prenatal diagnosis and extremely variable postnatal outcomes. On the basis of clinical manifestations, epidemiologic charateristics, and the presence and type of additional malformations, gastroschisis could be considered a heterogeneous condition with no gene/s discovered yet. This congenital anomaly affects approximately 1-3 infancts per 10,000 live births [Calzolari et al.1995;Parker et al.,2010] Current knowledge about causative mutations/variants. To date, no single gene has been linked to gastroschisis. Some publications have tried to link this malformation to variants in genes (such as AEBP1 (adipocyte enhancer binding protein) gene [Feldkamp et al,. 2012] or the VEGF-NOS3 pathway [Lammer et al., 2008]. Previously, a Scribble mutant mouse model (circletail) was reported to exhibit gastroschisis, however recent studies demonstrated that the Scribble knockout fetus exhibits exomphalos phenotype of gastroschisis [Carnagham et al., 2013]. This data is part of a pre-publication release. For information on the proper use of pre-publication data shared by the Wellcome Trust Sanger Institute (including details of any publication moratoria), please see http://www.sanger.ac.uk/datasharing/

  Study EGAS00001002664

• Allele Balance Bias Identifies Systematic Genotyping Errors and False Disease Associations

  In recent years, Next Generation Sequencing (NGS) has become a cornerstone of clinical genetics and diagnostics. Many clinical applications require high precision, especially if rare events such as somatic mutations in cancer or genetic variants causing rare diseases need to be identified. Although random sequencing errors can be modeled statistically and deep sequencing minimizes their impact, systematic errors remain a problem even at high depth of coverage. Understanding their source is crucial to increase precision of clinical NGS applications. In this work, we studied the relation between recurrent biases in allele balance (AB), systematic errors and false positive variant calls across a large cohort of human samples analyzed by whole exome sequencing (WES). We have modeled the allele balance distribution for biallelic genotypes in 987 WES samples in order to identify positions recurrently deviating significantly from the expectation, a phenomenon we termed allele balance bias (ABB). Furthermore, we have developed a genotype callability score based on ABB for all positions of the human exome, which detects false positive variant calls that passed state-of-the-art filters. Finally, we demonstrate the use of ABB for detection of false associations proposed by rare variant association studies (RVAS).

  Study EGAS00001003027

• Exome sequencing identifies mutation of the ribosome in T-cell acute lymphoblastic leukemia

  Acute lymphoblastic leukemia (ALL) is the most common pediatric malignancy and responds well to therapy in children, but not in adults. We used exome sequencing of 67 T-cell ALL (T- ALL) samples to gain insight in the mutational spectrum and age related differences of these leukemias. We detected protein-altering mutations in 508 genes, with an average of 8.2 mutations in pediatric and 21.0 in adult T-ALL. Based on stringent filtering, we predict 15 of these genes to be drivers, including 7 genes that were not previously implicated in T-ALL. We identify CNOT3 as a new tumor suppressor, as it is mutated in 12.2 % of adult T-ALLs and its knock-down causes de novo tumors in a sensitized drosophila eye tumor model. In addition, we identify mutations in RPL5 and RPL10, proteins of the 60S ribosomal complex2, in 18.2 % of pediatric T-ALLs. Mutation of arginine 98 in the X-linked RPL10 gene was most frequent. Yeast cells expressing the Rpl10 R98S mutant showed a clear ribosome biogenesis defect, excluding that this mutation is a passenger mutation. Our data provide insights in the mutational landscape of pediatric and adult T-ALL and identify the ribosome as a potential oncogenic factor in cancer.

  Study EGAS00001000296

• The Genomic Landscape of Core-Binding Factor Acute Myeloid Leukemias

  Pediatric de novo acute myeloid leukemia (AML) is a heterogeneous disease that can be divided into clinically distinct subtypes based on the presence of specific chromosomal abnormalities or gene alterations. One of the best characterized subtypes of AML involves leukemias with alterations of the core-binding factor (CBF)-complex, which comprises the FAB subtypes M2 and M4Eo and associates with a favorable outcome. Patients with the AML M2 subtype harbor a translocation between chromosomes 8 and 21 [t(8;21)] that yields the chimeric fusion gene RUNX1(AML1)-RUNX1T1(ETO), while patients with AML M4Eo express the chimeric fusion gene CBFβ-SMMHC(MYH11) as a result of an inversion/translocation event of chromosome 16 [inv(16)/t(16;16)]. In an effort to define the total complement of genetic changes in CBF-leukemia, we performed paired-end whole genome sequencing (WGS) on diagnostic leukemia blasts and matched germ line samples from 17 pediatric CBF-leukemia patients using the Illumina platform. Somatic alterations, including single nucleotide variations (SNVs) and structural variations (SVs), including insertions, deletions, inversions, and inter- and intra-chromosomal rearrangements, were detected using complementary analysis pipelines (Bambino, CREST and CONSERTING). Recurrent screening of identified mutations will be performed in a cohort of approximately 94 cases of CBF-leukemias.

  Study EGAS00001000349