14347 results for ""Pennsylvania+German"+ACS+study+'Socius'+open+access"

in 233.91 milliseconds.

• Genetic insights into the biological mechanisms governing human ovarian ageing

  Reproductive longevity is critical for fertility and impacts healthy ageing in women, yet insights into the underlying biological mechanisms and treatments to preserve it are limited. Here, we identify 290 genetic determinants of ovarian ageing, assessed using normal variation in age at natural menopause (ANM) in ~200,000 women of European ancestry. These common alleles influence clinical extremes of ANM; women in the top 1% of genetic susceptibility have an equivalent risk of premature ovarian insufficiency to those carrying monogenic FMR1 premutations. Identified loci implicate a broad range of DNA damage response (DDR) processes and include loss-of-function variants in key DDR genes. Integration with experimental models demonstrates that these DDR processes act across the life-course to shape the ovarian reserve and its rate of depletion. Furthermore, we demonstrate that experimental manipulation of DDR pathways highlighted by human genetics increase fertility and extend reproductive life in mice. Causal inference analyses using the identified genetic variants indicates that extending reproductive life in women improves bone health and reduces risk of type 2 diabetes, but increases risks of hormone-sensitive cancers. These findings provide insight into the mechanisms governing ovarian ageing, when they act across the life-course, and how they might be targeted by therapeutic approaches to extend fertility and prevent disease.

  Study EGAS00001004947

• Control of Focal Adhesion Kinase Activation by RUNX1-regulated miRNAs in high-risk AML

  We recently described a 16-gene expression signature for improved risk stratification of acute myeloid leukemia (AML) patients called the AML Prognostic Score (APS). A subset of APS high-risk AML patients showed increased levels of focal adhesion kinase (FAK), encoded by the Protein Tyrosine Kinase 2 (PTK2) gene, which was correlated with RUNX1 mutations. RUNX1 mutant cells are more sensitive to PTK2 inhibitors. As we were not able to detect RUNX1 binding sites in the PTK2 promoter, we hypothesized that RUNX1 might regulate micro(mi)RNAs that repress PTK2, such that loss-of-function RUNX1 mutations would result in reduced miRNA expression and derepression of PTK2. Examination of paired RNA-seq and miRNA-seq data from 301 AML cases revealed two miRNAs that positively correlated with RUNX1 expression, contained RUNX1 binding sites in their promoters and were predicted to target PTK2. We show that the hsa-let7a-2-3p and hsa-miR-135a-5p promoters are regulated by RUNX1, and that PTK2 is a direct target of both miRNAs. Even in the absence of RUNX1 mutations, hsa-let7a-2-3p and hsa-miR-135a-5p regulate PTK2 expression, and reduced expression of these two miRNAs sensitizes AML cells to PTK2 inhibition. These data explain how RUNX1 regulates PTK2, and identify potential miRNA biomarkers for targeting AML with PTK2 inhibitors.

  Study EGAS00001006491

• BinDel: software tool for detecting clinically significant microdeletions in low-coverage WGS-based NIPT samples

  Clinically pathogenic chromosomal microdeletions causing genetic disorders such as DiGeorge syndrome are rare genetic aberrations that can cause clinically relevant fetal and childhood developmental deficiencies. Clinical severity of such deficiencies depend on the exact genomic location and genes affected by the fetal chromosomal aberration. Here we present the BinDel, a novel region-aware microdeletion detection software package developed to infer clinically relevant microdeletion risk in low-coverage whole-genome sequencing NIPT data. To test BinDel, we quantified the impact of sequencing coverage, fetal DNA fraction, and region length on microdeletion risk detection accuracy. We also estimated BinDel accuracy on known microdeletion samples and clinically validated aneuploidy samples. BinDel identified each positive control sample as high risk. We also determined that it is critical to take into account that the sample with a detected high microdeletion risk does not have a full chromosome aneuploidy, as the latter can cause erroneous high microdeletion risk findings. We observed that lower sequencing coverage resulted in reduced microdeletion detection accuracy, and higher fetal fractions considerably increased the microdeletion detection accuracy, with coverage becoming increasingly relevant as fetal DNA fraction decreased. In conclusion, we developed an R package-based software tool BinDel for inferring fetal microdeletion risks, which accurately identified all positive control samples with microdeletion or -duplication aberrations as high-risk samples.

  Study EGAS00001006663

• Uncovering tumor intrinsic and extrinsic factors that regulate hepatocellular carcinoma growth using patient derived xenograft assays

  Several major factors limit our understanding of hepatocellular carcinoma (HCC). First, human HCCs are infrequently biopsied for diagnosis and thus are not often biologically interrogated. Second, HCC initiation and progression are strongly influenced by the cirrhotic microenvironment, and the exact contributions of intrinsic and extrinsic tumor factors are unclear. A powerful approach to examine the personalized biology of liver cancers and the influence of host tissues is with patient derived xenograft (PDX) models. In Asia, HCCs from hepatitis B virus patients have been efficiently converted into PDXs, but few parallel efforts from the west have been reported. In a large-scale analysis, we implanted 93 HCCs and 8 cholangiocarcinomas (CCAs) to systematically analyze host factors and to define an optimized platform for PDX development from both surgical and biopsy samples. NSG mice that had undergone partial hepatectomy (PHx) represented the best combination of engraftability, growth, and passageability, but overall rates were low and indicative of a unique intrinsic biology for HCCs in the US. PDX models preserved the histology and genetic features of parental tumors, and ultimately, 9 new models were usable for pre-clinical studies. Intriguingly, HCC PDXs were differentially sensitive to regorafenib and sorafenib and CCA PDXs were also highly sensitive to regorafenib. PDX models functionalize early and advanced stage HCCs and revealed unique biological features of liver cancers from the US.

  Study EGAS00001004020

• Cell-free DNA fragmentation patterns and personalized sequencing reveal circulating tumor DNA in urine and plasma of glioma patients

  Glioma-derived cell-free DNA (cfDNA) is challenging to detect using liquid biopsy as levels in body fluids are low. We determined the glioma-derived DNA fractions in tumor biopsies, in cerebrospinal fluid (CSF), plasma and urine samples, using deep sequencing of personalized capture panels. By sequencing cfDNA across thousands of mutations identified individually in each patient’s matched tumor we detected tumor-derived DNA in plasma (10/12) and urine samples (8/11). The median tumor fraction was 6.4x10-3 in CSF, 3.1x10-5 in plasma and 4.7x10-5 in urine. We identified a shift in the size distribution for mutant cfDNA fragments in these body fluids. Next, we analyzed cfDNA fragment sizes with paired-end shallow whole genome sequencing (WGS) in urine samples from 35 patients with gliomas, 8 individuals with non-malignant brain disorders, and 26 controls (n=69 individuals, 96 samples). cfDNA in urine of glioma patients was significantly more fragmented compared to urine from patients with non-malignant brain disorders (t-test, p=1.7x10-2) and compared to urine of controls (t-test, p=5.2x10-9). The proportion of DNA fragments <60 bp was higher in glioma patients urine and could be used for classification (AUC=0.93). Machine learning models integrating fragment lengths could identify urine samples from glioma patients (AUC=0.97 in cross-validation).

  Study EGAS00001004355

• Loss of LGR4/GPR48 causes severe neonatal salt-wasting due to disrupted WNT signaling altering adrenal zonation

  Disorders of isolated mineralocorticoid deficiency causing potentially life-threatening salt-wasting crisis early in life have been associated with gene variants of aldosterone biosynthesis or resistance, but in some patients no such variants are found. WNT/β-catenin signaling is crucial for differentiation and maintenance of the aldosterone producing adrenal zona glomerulosa (zG). We describe a highly consanguineous family with multiple perinatal deaths or infants presenting at birth with failure to thrive, severe salt-wasting crises associated with isolated hypoaldosteronism, nail anomalies, short stature, and deafness. Whole exome sequencing revealed a homozygous splice variant in the R-SPONDIN receptor LGR4 gene (c.618-1G>C) regulating WNT signaling. The resulting transcripts affected protein function and stability, and resulted in loss of Wnt/β-catenin signaling in vitro. The impact of LGR4 inactivation was analyzed by adrenal cortex specific ablation of Lgr4, using Lgr4Flox/Flox mated with Sf1:Cre mice. Inactivation of Lgr4 within the adrenal cortex in the mouse model caused decreased WNT signaling, aberrant zonation with deficient zG and reduced aldosterone production. Thus, human LGR4 mutations establish a direct link between LGR4 inactivation and decreased canonical WNT signaling with abnormal zG differentiation and endocrine function. Therefore, variants in WNT signaling and its regulators should systematically be considered in familial hyperreninemic hypoaldosteronism.

  Study EGAS00001006808

• Dynamics of sequence and structural cell-free DNA landscapes in small-cell lung cancer

  Patients with small-cell lung cancer (SCLC) have an exceptionally poor prognosis, calling for improved real-time non-invasive biomarkers of therapeutic response. We performed targeted error-correction sequencing on 171 serial plasmas and matched white blood cell (WBC) DNA from 33 patients with metastatic SCLC who received systemic treatment with chemotherapy or immunotherapy-containing regimens. Tumor-derived sequence alterations and plasma aneuploidy were evaluated serially and combined to assess changes in total cell-free tumor load (cfTL). Longitudinal dynamic changes in cfTL were monitored to determine circulating cell-free tumor DNA (ctDNA) molecular response during therapy. Combined tiered analyses of tumor-derived sequence alterations and plasma aneuploidy allowed for assessment of ctDNA molecular response in all patients. Molecular responses captured the therapeutic effect and long-term clinical outcomes in a more accurate and rapid manner compared to radiographic imaging. Patients with sustained molecular responses had longer overall (log-rank p=0.0006) and progression-free (log-rank p<0.0001) survival, with molecular responses detected on average 4 weeks earlier than imaging. ctDNA analyses provide a rapid approach for the assessment of early on-therapy molecular responses and have important implications for the management of patients with SCLC, including the development of improved strategies for real-time tumor burden monitoring.

  Study EGAS00001006831

• cfDNA dataset from extracelular vesicles and paired non-fractionated samples

  Liquid biopsies contain multiple analytes that can be mined to improve the detection and management of cancer. Beyond cell-free DNA (cfDNA), mutations have been detected in DNA associated with extracellular vesicles (EV-DNA). The genome-wide composition and structure of EV-DNA remains poorly characterized, and whether circulating EVs are enriched in tumor signal compared to unfractionated cfDNA is still unclear. Here, using whole genome sequencing from selected lung cancer patients with a relatively high cfDNA tumor content >5%, we determined that the tumor fraction and heterogeneity are comparable between DNA associated with small (<200 nm) EVs and matched plasma cfDNA. DNA in EV fractions, obtained with standardized size-exclusion chromatography, are comprised of short ~150-180 bp fragments and long >1000 bp fragments that are poor in tumor signal. Other fractions only exhibit short fragments with similar tumor DNA content. The composition in bases at the end of EV-DNA fragments, as well as their fragmentation patterns are similar to total plasma cfDNA. Mitochondrial DNA however is relatively enriched in EV fractions. Our results suggests that cfDNA in plasma may be of dual nature, either bound to proteins (including the nucleosome) but also associated to EVs. cfDNA associated to small EV (including exosomes) is however not preferentially enriched in tumor signal.

  Study EGAS00001006848

• WES data (cfTrack study) from NSCLC patients and OC patients

  We also collected samples from 8 NSCLC patients and 4 ovarian cancer patients and. For all 8 NSCLC patients, a tumor biopsy sample, a WBC sample, and three plasma samples were collected. For all 4 ovarian cancer patients, a WBC sample and two serum samples were collected. We collected tumor tissue sample from one ovarian cancer patient (OV4). The cfDNA was extracted from their plasma samples using the QIAamp circulating nucleic acid kit from QIAGEN (Germantown, MD). For serum cfDNA, ampure XP beads size selection was further performed to eliminate gDNA contamination. In brief, 0.5 volume of beads were first added to the cfDNA samples. After incubation, the supernatant was transferred to a new tube and an additional 2.0 volume of beads were added. After 80% ethanol wash, cfDNA was eluted from the beads. FA assays (Agilent Technologies) were performed to rule out the contamination of gDNA in the size selected samples. The cfDNA WES library of all patients and the genomic DNA WES library of the 4 ovarian cancer patients were constructed with the SureSelect XT HS kit from Agilent Technologies (Santa Clara, CA) according to the manufacturer’s protocol. In brief, 10ng of cfDNA was used as input material. After end repair/dA-tailing of cfDNA, the adaptor was ligated. The ligation product was purified with Ampure XP beads (Beckman-Coulter, Atlanta, GA) and the adaptor-ligated library was amplified with index primer in 10-cycle PCR. The amplified library was purified again with Ampure XP beads, and the amount of amplified DNA was measured using the Qubit 1xdsDNA HS assay kit (ThermoFisher, Waltham, MA). 700-1000 ng of DNA sample was hybridized to the capture library and pulled down by streptavidin-coated beads. After washing the beads, the DNA library captured on the beads was re-amplified with 10-cycle PCR. The final libraries were purified by Ampure XP beads. The library concentration was measured by Qubit, and the quality was further examined with Agilent Bioanalyzer before the final step of 2x150bp paired-end sequencing at an average coverage of 200. Whole-exome capture libraries of genomic DNA from the 8 NSCLC patients were constructed via Roche SeqCap EZ Exome V6 (Roche). Enriched exome libraries were sequenced on the Illumina HiSeq 3000 platform (Illumina) to generate 2x100bp paired-end reads at an average coverage of 200.

  Dataset EGAD00001008454

• RRBS data (cfSort study) from noncancer tissue DNA

  The RRBS libraries of the genomic DNA from the 521 tissue samples were constructed following the standard RRBS protocol. 100-200 ng of intact genomic DNA in the volume of 21.5 µl was used as input material. Restriction digestion was done with 2.5 µl 10xCutSmart buffer and 1 µl MspI (NEB) for 18 h at 37 oC and 20 min at 65 oC. 0.5 µl 10xCutSmart buffer, 0.3 µl dACGTP mixture (100 mM dATP, 10 mM dCTP, 10 mM dGTP), 1 µl Klenow (exo-, 5U/µl, NEB) and 2.6 µl RT-PCR water, 0.6 µl 50 mM DTT (ThermoFisher) was added to the mixture for end repair and A-overhang addition with the program 30 oC for 20 min, 37 oC for 1 h and 75 oC for 20 min. Adapter ligation was then performed with 1 µl 10xThermoFisher HC T4 ligase buffer, 0.4 µl 100 mM ATP (ThermoFisher), 0.2 µl 50 mM DTT, 1 µl ThermoFisher HC T4 DNA ligase (30 Weiss Unit/µl), 30 ng home-made duplex UMI adapter with all the cytosines methylated (protocol adopted from Kennedy et al.) at 16 oC for 20 h and 65 oC for 20 min. Bisulfite conversion of the adapter-ligated product was carried out with QIAGEN EpiTect plus DNA bisulfite kit following their protocol for two rounds of conversion. The converted product was purified with Qiagen MinElute spin column and eluted with 20 µl RT-PCR water. PCR amplification was done using the NEBNext Multiplex Oligos for Illumina (2.5 µl of universal and index primer each) and 25 µl KAPA HiFi HotStart Uracil+ ReadyMix (Roche) with the following cycling conditions: 98 oC for 45 s, 9 cycles of 98 oC for 15 s, 60 oC for 30 s and 72 oC for 30 s, followed by a final extension at 72 oC for 5 min. The PCR product was purified with 1x AmpureXP beads and eluted with 30 µl EB buffer. DNA concentration was measured by Qubit 1xdsDNA HS assay. 5% TBE-UREA PAGE and bioanalyzer assay was performed as quality control on each library before sequencing.

  Dataset EGAD00001010880

• Double mutant DNMT3A AML: a unique subtype

  Mutation of DNMT3A, encoding a de novo methyltransferase essential for cytosine methylation, is a common early event in clonal hematopoiesis (CH) and adult acute myeloid leukemia (AML). Spontaneous deamination of methylated cytosines damages DNA, which is repaired by the base excision repair (BER) enzymes MBD4 and TDG. Congenital MBD4-deficiency has been linked to early-onset CH and AML, and is marked by exceedingly high levels of DNA damage and mutation of DNMT3A. Strikingly, wildtype (WT) DNMT3A binds TDG, thereby potentiating its repair activity. Since TDG is the only remaining BER enzyme in MBD4-deficient AML patients capable of repairing methylation damage, we investigated whether mutant DNMT3A negatively affects the repair function of TDG. We found that, whereas WT DNMT3A stimulates TDG function, mutant DNMT3A impairs TDG-mediated repair of DNA damage in vitro. In light of this finding and to extrapolate our observations to the broader AML patient population, we investigate here the genetic profiles and survival outcomes of AML patients with single (SM) versus double mutant (DM) DNMT3A. DM DNMT3A AML patients show a characteristic driver mutation landscape and reduced overall survival when compared to SM DNMT3A AML patients. Importantly, whole-genome sequencing showed a trend for increased DNA damage in primary DM DNMT3A AML samples, especially when DNMT3A mutations are located at the DNMT3A-TDG interaction interface.

  Study EGAS00001007966

• Genome wide variation in the Angolan Namib desert reveals unique Pre-Bantu ancestry

  The populations of the Angolan Namib Desert have been largely neglected in previous surveys of the genomic landscape of southern Africa. Although at present, the Namib is culturally dominated by Southwest Bantu-speaking cattle-herders, the region exhibits an extraordinary ethnographic diversity which includes an array of semi-nomadic peoples whose subsistence strategies fall outside the traditional division between foraging and food production and can thus be referred to as “peripatetic”. Among these small-scale populations are the last speakers of the Kwadi branch of the Khoe-Kwadi language family associated with the introduction of pastoralism into southern Africa (Kwepe), as well as a range of groups whose origins remain enigmatic (Kwisi, Twa and Tjimba). Using genome-wide data from 208 individuals belonging to nine ethnically diverse groups from the Angolan Namib and adjacent areas (Kwepe, Kwisi, Twa, Tjimba, !Xun, Kuvale, Himba, Nyaneka, Ovimbundu) in combination with published data from other regions of Africa, we reconstruct in detail the histories of contact emerging from pre-historic migrations to southern Africa and show that peripatetic groups from southwestern Angola stand out for exhibiting elevated levels of an unique, regionally-specific and highly divergent Pre-Bantu ancestry. These findings highlight the importance of the Namib for understanding the deep genetic structure of Africa.

  Study EGAS00001007011

• Germline whole genome sequencing of patients with Li-Fraumeni syndrome

  Five hundred ng to 1 ug of genomic DNA was submitted to The Centre for Applied Genomics (TCAG) at The Hospital for Sick Children for genomic library preparation and whole genome sequencing. DNA samples were quantified using the Qubit High Sensitivity Assay and purity was assessed using the Nanodrop OD 260/280 ratio. Approximately 500-700 ng of DNA was used as input material for library preparation using the Illumina TruSeq PCR-free DNA Library Prep Kit following the manufacturer’s recommended protocol. In brief, DNA was fragmented to 400 bp on average using sonication on a Covaris LE220 instrument. Fragmented DNA was end-repaired, A-tailed and indexed TruSeq Illumina adapters with overhang-T were added to the DNA. Libraries were then validated on a Bioanalyzer DNA High Sensitivity chip to check for size and absence of primer dimers and quantified by qPCR using Kapa Library Quantification Illumina/ABI Prism Kit protocol (KAPA Biosystems). Validated libraries were pooled in equimolar quantities and paired-end sequenced on an Illumina HiSeq X platform following Illumina’s recommended protocol to generate paired-end reads of 150-bases in length and an average depth of 40X. FASTQ files were aligned to the hg19 reference genome using BWA-mem v0.78. PCR duplicates were marked with Picard MarkDuplicates v1.1.08, and base recalibration and realignment was performed using GATK v2.8.1.

  Study EGAS00001007061

• BMP4 and temozolomide synergize in the majority of patient derived glioblastoma cultures

  The poor prognosis of patients with glioblastoma is partly due to resistance to standard treatment, which combines chemoradiation using temozolomide (TMZ). The concept of cancer stem cells provided new insight into tumor resistance and management, and isolated Glioblastoma stem-like cells (GSCs) which contribute to this therapy resistance. Bone morphogenetic protein 4 (BMP4) is a differentiation factor that stimulates astroglial differentiation of such GSCs, thereby reducing their self-renewal. Exposing GSCs to BMP4 could therefore sensitize them to chemotherapy with TMZ, but previous studies with various treatment paradigms have produced contradictory results in this matter. Therefore, we re-defined the optimal treatment paradigm for the combination of TMZ and BMP4 and subsequently re-assessed the inter-tumoral variability in response to co-treatment with TMZ+BMP4 in vitro. First, the simultaneous treatment was found more effective than sequential therapy. Second, after application of our optimized treatment protocol, 70% of a total of 20 patient-derived glioblastoma cultures showed synergy when the treatments were combined. The synergistic effect was a result of increased cellular apoptosis and not of decreased proliferation. Comparative RNA-sequencing suggested that combination treatment results in decreased MAPK signaling. We conclude that addition of BMP4 to standard therapy with TMZ could enlarge the group of patients with good therapeutic response.

  Study EGAS00001007095

• Structural variant analysis of homologous recombination-deficient genomes

  Homologous recombination (HR) deficiency causes DNA breaks and cytogenetic aberrations. Paradoxically, the types of DNA rearrangements specifically associated with HR-deficient cancers only minimally impact chromosomal structure. Addressing this, we combined a genome graph analysis of short-read whole genome sequencing (WGS) profiles across thousands of tumors with deep linked-read (LR) WGS of 46 BRCA1 or BRCA2 mutant breast cancers to discover a distinct class of HR deficiency-enriched rearrangements called reciprocal pairs. LR WGS showed that reciprocal pairs with identical rearrangement orientations gave rise to one of two distinct chromosomal outcomes, distinguishable only with long molecule data. While one (cis) outcome corresponded to the copy and pasting of a small segment to a distant site, a second (trans) outcome was a quasi-balanced translocation or multi-megabase inversion with substantial (10kb) duplications at each junction. The full spectrum of reciprocal pair outcomes could be explained by an HR-independent replication restart repair mechanism. LR WGS additionally identified single-strand annealing (SSA) as a BRCA2-deficiency specific repair pathway in human cancers. Replication restart- and SSA-associated SVs improved BRCA1- vs. BRCA2- deficiency classification and identified metastatic cancer cases with favorable chemotherapy responses. Our data reveal classes of BRCA1- and BRCA2-deficiency specific rearrangements as drivers of cytogenetic aberrations in HR deficient cells.

  Study EGAS00001007186

• deep-learning-powered tissue deconvolution for cfDNA

  Plasma cell-free DNA (cfDNA) is a noninvasive biomarker for cell death of all organs. Deciphering the tissue origin of cfDNA can reveal abnormal cell death because of diseases, which has great clinical potential in disease detection and monitoring. Despite the great promise, the sensitive and accurate quantification of tissue-derived cfDNA remains challenging to existing methods due to the limited characterization of tissue methylation and the reliance on unsupervised methods. To fully exploit the clinical potential of tissue-derived cfDNA, here we present one of the largest comprehensive and high-resolution methylation atlas based on 521 noncancer tissue samples spanning 29 major types of human tissues. We systematically identified fragment-level tissue-specific methylation patterns and extensively validated them in orthogonal datasets. Based on the rich tissue methylation atlas, we develop the first supervised tissue deconvolution approach, a deep-learning-powered model, cfSort, for sensitive and accurate tissue deconvolution in cfDNA. On the benchmarking data, cfSort showed superior sensitivity and accuracy compared to the existing methods. We further demonstrated the clinical utilities of cfSort with two potential applications: aiding disease diagnosis and monitoring treatment side effects. The tissue-derived cfDNA fraction estimated from cfSort reflected the clinical outcomes of the patients. In summary, the tissue methylation atlas and cfSort enhanced the performance of tissue deconvolution in cfDNA, thus facilitating cfDNA-based disease detection and longitudinal treatment monitoring.

  Study EGAS00001007213

• ATAC-Seq of CD4 T cell subsets

  Blood samples were taken from 4 age-matched healthy volunteers (two males and two females), and naïve CD4 T cells were isolated using the naïve CD4+ T Cell Isolation Kit II, human (Miltenyi Biotec). Two hundred thousand cells were plated in a 96-well plate, activated with Dynabeads Human T-activator CD3/CD28 (bead-to-cell ratio of 1:1, as recommended by the manufacturer) (ThermoFisher Scientific) and polarized towards CD4 T cell subtypes (Teff, Th1, Th2, rTh17, pTh17, and Treg cells) with specific cocktails of recombinant cytokines and neutralizing antibodies (see details in the material and methods) for each subset. At day 4 post-stimulation, forty thousand cells were used to perform ATAC-Seq on CD4 T cell subtypes. Again, crude nuclei of live CD4+ T cell subsets were treated with Tagment DNA buffer and Tagment DNA Enzyme (Nextera DNA Library Prep Kit, Illumina), and then the DNA was purified by MinElute PCR Purification Kit (Qiagen). Transposed DNA fragments were amplified using specific adapters followed by purification with MinElute PCR Purification Kit (Qiagen). Fragments from 240-360pb were selected in the PippinHT system (Sage Science). The quality of the library and its DNA concentration were assessed by Bioanalyzer instruments (Agilent Technologies) and ultimately submitted for sequencing using Illumina HiSeq 2500 sequencer, V4 chemistry.

  Study EGAS00001007345

• Machine Learning Signal Enrichment for Ultrasensitive Plasma Tumor Burden Monitoring

  In solid tumor oncology, circulating tumor DNA (ctDNA) is poised to transform care through accurate assessment of minimal residual disease (MRD) and therapeutic response monitoring. To overcome the sparsity of ctDNA fragments in low tumor fraction (TF) settings and increase MRD sensitivity, we previously leveraged genome-wide mutational integration through plasma whole genome sequencing (WGS). We now introduce MRD-EDGE, a composite machine learning-guided WGS ctDNA single nucleotide variant (SNV) and copy number variant (CNV) detection platform designed to increase signal enrichment. MRD-EDGESNV uses deep learning and a ctDNA-specific feature space to increase SNV signal-to-noise enrichment in WGS by 300X compared to previous WGS error suppression. MRD-EDGECNV also reduces the degree of aneuploidy needed for ultrasensitive CNV detection through WGS from 1 Gb to 200 Mb, thereby expanding its applicability to a wider range of solid tumors. We harness the improved performance to identify MRD following surgery in multiple cancer types, track changes in tumor burden in response to neoadjuvant immunotherapy in non-small cell lung cancer (NSCLC) and demonstrate ctDNA shedding in precancerous colorectal adenomas. Finally, the radical signal-to-noise enrichment in MRD-EDGESNV enables de novo mutation calling in melanoma and small cell lung cancer (SCLC) without matched tumor, yielding clinically informative TF monitoring for patients on immune checkpoint inhibition (ICI).

  Study EGAS00001007451

• Machine learning guided signal enrichment for ultrasensitive plasma tumor burden monitoring

  In solid tumor oncology, circulating tumor DNA (ctDNA) is poised to transform care through accurate assessment of minimal residual disease (MRD) and therapeutic response monitoring. To overcome the sparsity of ctDNA fragments in low tumor fraction (TF) settings and increase MRD sensitivity, we previously leveraged genome-wide mutational integration through plasma whole genome sequencing (WGS). We now introduce MRD-EDGE, a composite machine learning-guided WGS ctDNA single nucleotide variant (SNV) and copy number variant (CNV) detection platform designed to increase signal enrichment. MRD-EDGESNV uses deep learning and a ctDNA-specific feature space to increase SNV signal-to-noise enrichment in WGS by 300X compared to previous WGS error suppression. MRD-EDGECNV also reduces the degree of aneuploidy needed for ultrasensitive CNV detection through WGS from 1 Gb to 200 Mb, thereby expanding its applicability to a wider range of solid tumors. We harness the improved performance to identify MRD following surgery in multiple cancer types, track changes in tumor burden in response to neoadjuvant immunotherapy in non-small cell lung cancer (NSCLC) and demonstrate ctDNA shedding in precancerous colorectal adenomas. Finally, the radical signal-to-noise enrichment in MRD-EDGESNV enables de novo mutation calling in melanoma and small cell lung cancer (SCLC) without matched tumor, yielding clinically informative TF monitoring for patients on immune checkpoint inhibition (ICI).

  Study EGAS00001007545

• The double-hit signature identifies double-hit diffuse large B-cell lymphoma with genetic events cryptic to FISH​

  High-grade B-cell lymphomas with MYC and BCL2 and/or BCL6 rearrangements (HGBL-DH/THs) include a group of diffuse large B-cell lymphomas (DLBCLs) with inferior outcomes after standard chemoimmunotherapy. We recently described a gene expression signature that identifies 27% of germinal center B-cell DLBCLs (GCB-DLBCLs) as having a double-hit-like expression pattern (DHITsig) and inferior outcomes; however, only half of these cases have both MYC and BCL2 translocations identifiable using standard breakapart fluorescence in situ hybridization (FISH). Here, 20 DHITsig+ GCB-DLBCLs apparently lacking MYC and/or BCL2 rearrangements underwent whole-genome sequencing. This revealed 6 tumors with MYC or BCL2 rearrangements that were cryptic to breakapart FISH. Copy-number analysis identified 3 tumors with MYC and 6 tumors with MIR17HG gains or amplifications, both of which may contribute to dysregulation of MYC and its downstream pathways. Focal deletions of the PVT1 promoter were observed exclusively among DHITsig+ tumors lacking MYC translocations; this may also contribute to MYC overexpression. These results highlight that FISH fails to identify all HGBL-DH/THs, while revealing a range of other genetic mechanisms potentially underlying MYC dysregulation in DHITsig+ DLBCL, suggesting that gene expression profiling is more sensitive for identifying the biology underlying poor outcomes in GCB-DLBCL.

  Study EGAS00001004285

• Bulk-tissue RNA-sequencing paired nuclear and cytoplasmic fractions of anterior prefrontal cortex, cerebellar cortex and putamen tissues from post-mortem neuropathologically-confirmed control individuals.

  Gaining insight into the genetic regulation of gene expression in human brain is key to the interpretation of genome-wide association studies for major neurological and neuropsychiatric diseases. Expression quantitative trait loci (eQTL) analyses have largely been used to achieve this, providing valuable insights into the genetic regulation of steady-state RNA in human brain, but not distinguishing between molecular processes regulating transcription and stability. RNA quantification within cellular fractions can disentangle these processes in cell types and tissues which are challenging to model in vitro. We investigated the underlying molecular processes driving the genetic regulation of gene expression specific to a cellular fraction using allele-specific expression (ASE). Applying ASE analysis to genomic and transcriptomic data from paired nuclear and cytoplasmic fractions of anterior prefrontal cortex, cerebellar cortex and putamen tissues from 4 post-mortem neuropathologically-confirmed control human brains, we demonstrate that a significant proportion of genetic regulation of gene expression occurs post-transcriptionally in the cytoplasm, with genes undergoing this form of regulation more likely to be synaptic. These findings have implications for understanding the structure of gene expression regulation in human brain, and importantly the interpretation of rapidly growing single-nucleus brain RNA-sequencing and eQTL datasets, where cytoplasm-specific regulatory events could be missed.

  Study EGAS00001006380

• Characterization of individual foci of multicentric/multifocal breast cancer using targeted next generation sequencing

  Multifocality or multicentricity in breast cancer may be defined as the presence of two or more tumor foci within a single quadrant of the breast or within different quadrants of the same breast, respectively. This original classification of the breast cancer as multicentric or multifocal was based on the assumption that cancers arising in the same quadrant were more likely to arise from the same ductal structures than those occurring in separate areas of the breast. The problem with these definitions is that the ?quadrants? of the breast are arbitrary external designations, as no internal boundaries do exist. This project will therefore focus both on synchronous multifocal and multicentric tumors. The incidence of multifocal and multicentric breast cancers was reported to be between 13 and 75% depending on the definition used, the extent of the pathologic sampling of the breast and whether in situ disease is considered evidence of multicentricity (1). Although this incidence is variable, those figures show that it is a frequent phenomenon. Multiple (multifocal/multicentric) breast carcinomas, especially when occurring in the same breast, represent a real challenge for both pathologists and clinicians in terms of identifying the cellular origin and the best therapeutic management of the cancer. Multifocality or multicentricity has been associated with a number of more aggressive features including an increased rate of regional lymph node metastases and adverse patient outcome when compared with unifocal tumors (2-3), and a possible increased risk of local recurrence following breast conserving surgery (4). For the moment, the literature is divided on whether there is a corresponding impact on survival outcomes. Today, the current convention to stage and to treat multifocal and multicentric tumors is the classical tumor-node-metastasis (TNM) staging guidelines with which tumor size is assessed by the largest tumor focus without taking other foci of disease into consideration. If some papers, as the recent one from Lynch and colleagues, support the current staging convention (3), others, however, as Boyages et al. suggested that aggregate size and not the size of the largest lesion should be considered in order to refine the prognostic assessment of those tumors (5). On the top of that, the question whether multifocal/multicentric carcinomas are due to the spread of a single carcinoma throughout the breast or is due to multiple carcinomas arising simultaneously has been a matter of debate. Some studies suggested that multifocal breast cancer may result from either intramammary spread from a single primary tumor or multiple synchronous primary tumors; whereas others suggest that multiple breast carcinomas always arise from the same clone (6-8). Recently, Pietri and colleagues analyzed the biological characterization of a series of 113 multifocal/multicentric breast cancers (8) which were diagnosed over a 5-year period. The expression of estrogen (ER) and progesterone (PgR) receptors, Ki-67 proliferative index, expression of HER2 and tumor grading were prospectively determined in each tumor focus, and mismatches among foci were recorded. Mismatches in ER status were present in 5 (4.4%) cases and PgR in 18 (15.9%) cases. Mismatches in tumor grading were present in 21 cases (18.6%), proliferative index (Ki-67) in 17 (15%) cases and HER2 status in 11 (9.7%) cases. Interestingly, this heterogeneity among foci has led to 14 (12.4%) patients receiving different adjuvant treatments compared with what would have been indicated if we had only taken into account the biologic status of the primary tumor. This study therefore showed that differences in biological characteristics of multifocal/multicentric lesions play a crucial role in the adjuvant treatment decision making process. In this study, we will concentrate on a larger series of patients with multifocal invasive ductal breast cancer lesions. We aim at: 1. Evaluating the incidence of multifocality according to the different breast cancer molecular subtypes (ER-/HER2-, HER2+, ER+/HER2-). 2. Evaluating the incidence of multifocality in patients with hereditary breast cancer disease (presence of germline BRCA1 or BRCA2 mutations). Moreover, we would like to investigate if multifocal lesions with BRCA1 or BRCA2 mutations exhibit a characteristic combination of substitution mutation signatures and a distinctive profile of deletions as demonstrated recently by Nik-Zainal and colleagues (9). 3. Correlating multifocality with clinical information in order to define its influence on patients? survival (DFS and OS). 4. Carrying high coverage targeted gene sequencing of driver cancer genes and genes whose mutation is of therapeutic importance in order to compare clinically-relevant genetic differences between several multifocal breast cancer lesions. 5. Evaluating the impact of the distance between the different lesions on the clinical outcome but also on the genetic differences. 6. Comparing gene expression patterns between several multifocal breast cancer lesions and correlate them with the results of the targeted genes screen. 7. Characterizing the genomic and transcriptomic status of cancer related genes in metastatic lesions (local recurrence, positive lymph node or distant metastatic sites) from the same multifocal invasive ductal breast cancer patients in order to evaluate the consequence of genomic and transcriptomic heterogeneity of multifocal lesions on metastatic lesions. Multiple (multifocal/multicentric) breast carcinomas, especially when occurring in the same breast, represent a real challenge for both pathologists and clinicians in terms of identifying the cellular origin and the best therapeutic choice. This project has the potential to identify genetic/transcriptomic differences existing between several lesions constituting multifocal breast cancers, which in the routine clinical practice are usually considered to be homogeneous among them. We foresee validating significant results in a larger series of patients and this, in turn, could have a remarkable impact on the treatment and clinical management of multifocal breast cancers. Indeed, we hope to provide some evidence whether or not each focus matters in multifocal and multicentric breast cancer to define the adequate therapeutic approach, especially in the context of targeted therapies. The work to be done at Sanger will be target gene screen pooling of 1400 samples.

  Dataset EGAD00001000624

• The patterns and dynamics of genomic instability in metastatic pancreatic cancer

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

  Study EGAS00000000064

• Nanopore sequencing enables allelic phasing of FLG loss-of-function variants, intragenic copy number variation and methylation status in atopic dermatitis and ichthyosis vulgaris

  Loss-of-function (LoF) variants in the FLG gene are associated with ichthyosis vulgaris (IV) and atopic dermatitis (AD). IV and AD patients with two FLG LoF variants are reported to have a more severe disease course. However, subsets of compound heterozygous patients display phenotypic heterogeneity despite a prediction of a severe phenotype from their genotype. The underlying genetic mechanism of this phenotypic variability is unknown. One possibility is that two LoF FLG variants could be located in cis (on the same allele), thus altering the predicted gene dosage. To investigate the allelic contribution of compound heterozygous FLG LoF mutations to IV phenotypic variability, we employ long reads from Oxford Nanopore Technologies, utilising two methodologies to phase the entire coding region of FLG into parental alleles: an amplicon-sequencing approach and an adaptive-sampling strategy. We investigated the allelic phasing of FLG LoF variants and intragenic-repeat copy-number variations (CNV) of 21 subjects with known compound heterozygous genotypes and mixed IV severity. Our results suggest that IV severity is not linked to the presence of cis variants in our IV patient cohort. Nanopore reads enabled detection and accurate phasing of CNVs and LoF variants into parental alleles, providing a framework for future in-depth genetic analysis. These results highlight the usefulness of utilising long reads for interrogating the FLG gene and other highly repetitive genes for genetic stratification related to clinical phenotypes.

  Study EGAS50000000166

• The tissue origin of highly elevated cell-free DNA in patients with and without cancer

  Forty years ago, the first measurements of cell free DNA (cfDNA) from cancer patients revealed that its total concentration was often elevated compared to healthy controls. The basis for this increase, and in particular the source of the extra cfDNA in such patients, has never been determined. To address this issue, we assessed methylation patterns of cfDNA from 57 individuals, 34 of whom had unequivocally elevated cfDNA concentrations (i.e., >15 ng/mL of plasma). We found that the excess cfDNA in cancer patients did not come from either neoplastic cells or surrounding cells from the tumor's tissue of origin. Instead, 80% of the cfDNA was derived from neutrophils, B cells or T cells, with neutrophils accounting for ~2/3. B cells contributed more than T cells despite the great excess of T cells over B cells in the blood. To see if this result was generalizable to other clinical situations in which elevated cfDNA concentrations have been described, we evaluated nine patients before and one day following surgery. After surgical resection, cfDNA concentration increased by an average of 16.1-fold, as expected. The majority of this increase (55%) was contributed by leukocytes, though damaged hepatocytes contributed a substantial minority (31% of increase). These data suggest that the elevated concentrations of cfDNA in cancer patients, as with patients after surgery, arises largely from dying leukocytes—not from the neoplastic cells or the tissue of origin of the cancer.

  Study EGAS00001005400

• Characterizing the Role of the Immune Microenvironment in Multiple Myeloma Progression at a Single Cell Level

  The plasma cell disease multiple myeloma (MM) is highly dependent on its bone marrow microenvironment (BMME). Early alterations within the bone marrow microenvironment that contribute to the progression of multiple myeloma (MM) from its precursor stages could hold the key to identifying novel therapeutic approaches, yet the intrinsic variability in cellular populations between patients together with differences in sample processing and analysis methods have made it difficult to identify consistent changes between datasets. Here, we used single-cell RNA sequencing of BM cells from precursor stages, monoclonal gammopathy of unknown significance (MGUS), smoldering MM (SMM) and newly diagnosed MM and analyzed our data in combination with a previously published dataset that used a similar patient population and sample processing. Despite vast interpatient heterogeneity, some alterations were seen consistently in both datasets. We identified changes in immune cell populations as disease progressed that were characterized by a substantial decrease in memory and naïve CD4 T cells and an increase in CD8+ effector T cells and T-regulatory cells. These alterations were further accompanied by an enrichment of non-clonal memory B cells and an increase of CD14 and CD16 monocytes in MM compared to its precursor stages. These results provide crucial information of the immune changes associated with the progression to clinical MM and can help to develop immune based strategies for patient stratification and early therapeutic intervention. 

  Study phs002756

• Integrated Epigenetic Maps of Human Embryonic, Extraembryonic and Adult Cells

  Working cooperatively with other Mapping Centers and the Data Coordination Center (EDACC) funded by this Roadmap mechanism we have comprehensively mapped the epigenomes of select human cells with significant relevance to complex human disease. Our effort has been focused on cells relevant to human health and complex disease including cells from the blood, brain, breast, skin, male germ cells, extraembryonic tissues and trophoblast, and U.S. Government-approved lines of human embryonic stem cells. We have incorporated high quality, homogeneous cells from males and females, and two predominant racial groups, and biological replicates of each cell type. The resulting comprehensive maps include up to 6 histone modifications selected for their opposing roles in regulating active and inactive chromatin (H3K4me1, H3K4me3, H3K9me3, H3K9ac, H3K27ac, H3K27me3, H3K36me3), DNA methylation, miRNA and gene expression, and for select cases whole genome shotgun sequence. Our cell and data production pipeline incorporates verification and data validation with independent methods, and operates under a model motivated by increased data production and decrease cost. Our group in conjunction with the other REMC teams, the EDACC, ENCODE, Epigenetics of Human Health and Disease centers and the NIH Roadmap program are developing methods, tools and reference epigenome maps for the research community that will make the promise of epigenetics in understand and treating human complex disease a reality. Our reference epigenomes will enable new disciplines including human population epigenetics, comparative epigenomics, neuroepigenetics, and therapeutic epigenetics for tissue regeneration and reversal of disease.

  Study phs000791

• Cell-Free, Methylated DNA in Blood Samples Reveals Tissue-Specific, Cellular Damage from Radiation Treatment

  Radiation therapy is an effective cancer treatment, although damage to surrounding healthy tissues can also occur. Cell-type specific DNA methylation patterns can be used to decode the cellular origins of cell-free DNA fragments, allowing for minimally-invasive monitoring of tissue damage. To evaluate whether changes in cell-free DNA methylation can indicate damages to tissues in patients treated with radiation, we collected serum samples from 15 breast cancer patients at three timepoints during their standard-of-care radiation therapy after surgery. A baseline sample was taken for each patient before onset of radiation therapy and a second End-Of-Treatment (EOT) sample was taken 30 minutes after the last treatment fraction. Finally, a recovery sample was taken one month after completion of radiation therapy. From serum samples, cell-free DNA was isolated and bisulfite capture-sequencing DNA methylation data were generated. Bisulfite capture-sequencing DNA methylation data were also generated from buffy coat and tissue-specific endothelial gDNA samples as well as from serum/plasma samples of healthy controls. The analysis of serum samples from breast cancer patients undergoing radiation treatment revealed distinct tissue-specific epithelial and endothelial responses to radiation across multiple organs. In conclusion, as a proof of concept we show that cell-type specific methylation signatures can be applied to detect cellular injury from radiation treatment using minimally invasive cell-free DNA in blood samples.

  Study phs003290

• NIMH (National Institute of Mental Health) De Novo Mutation Identification in Taiwanese Schizophrenia Trios

  The substantial reproductive impact of schizophrenia, for which affected individuals have fewer than half as many offspring as unaffected individuals do, implies that mutations of largest effect will frequently be de novo mutations. Ascertaining exome sequence variation in father-mother-offspring trios allows such mutations to be identified and distinguished from the far-larger amount of rare variation that is inherited by each individual. The pursuit of this approach in a large, well-powered cohort of trios can also provide lessons that inform the development of such gene discovery strategies more generally in human genetics. Schizophrenia trios from the Taiwanese population are being collected by Dr. Ming Tsuang (PI, UC San Diego, California) and investigators in Taiwan (PI, Dr. Hai Gwo Hwu; both funded by NIMH grant 1R01MH085560; Expanding Rapid Ascertainment Networks of Schizophrenia Families in Taiwan). A total of 3800 trios are anticipated to be collected by May 2013. This represents a highly homogenous national sample from the same ancestral population. DNA samples will be obtained from the NIMH Repository, Rutgers University Cell and DNA Repository (described below) and stored at the Broad Institute. Genetic and data analyses will be performed at the Broad Institute. We propose to sequence the whole exome of trios by hybrid capture and Illumina next generation sequencing and perform targeted genotyping and validation of variants (SNPs, indels and CNVs) using several molecular methods, to include emulsion-based PCR and Sanger sequencing.

  Study phs001196

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

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

  Study phs001491

• RNA-Sequencing of B-Lymphoblastic Leukemia with Glucocorticoids and PI3K Delta Inhibition

  Mononuclear cells were isolated from bone marrow or peripheral blood of 7 pediatric subjects with B-lymphoblastic leukemia at the time of diagnosis. To determine the effect of PI3K delta inhibition on glucocorticoid-induced gene regulation, B-lymphoblastic leukemia cells were treated with glucocorticoids +/- the PI3K delta inhibitor, idelalisib, for 24 hours prior to RNA extraction and sequencing. Biological replicates were performed for 4 of the 7 subjects; 2 of these subjects had specimens demonstrating an additive response to the combination of glucocorticoids and idelalisib in viability testing in vitro, and 2 subjects' specimens demonstrated a synergistic response. All 7 subjects were analyzed together and the 4 specimens with biological replicates were also evaluated in an additional analysis. An average of 40 million reads per specimen and treatment condition were obtained. We compared gene expression in each treatment condition to vehicle-treated cells and expression in glucocorticoid or idelalisib only treated cells to combination treated cells. Using a high dose of dexamethasone (25 or 50 nM) and an equivalent dose of prednisolone, we found that the addition of idelalisib (500 nM) enhanced prednisolone-induced gene regulation but not dexamethasone-induced gene regulation. Surprisingly, the specimens with additive response to combination therapy demonstrated enhancement of prednisolone-induced gene regulation with the addition of idelalisib, while the specimens with synergistic responses did not. Raw sequencing data (FASTQ) and count tables will be available through dbGaP.

  Study phs003085

• NINDS Family-Based Whole-Genome Sequencing to Find Modifiers of Age of Onset in Huntington's Disease

  The goal of our studies is to identify genetic modifiers of neurodegeneration in Huntington's disease (HD). HD is caused by expansion of CAG repeats in the huntingtin (Htt) gene, with longer stretches often leading to more rapid disease onset and progression. Yet, for a given number of repeats, the age of symptom onset can be variable, differing by up to decades. Thus, the age of onset of motor symptoms in HD is only partly explained by the length of the CAG expansion. Available evidence suggests that genetic modifiers contribute to the variation in HD onset. Identifying genetic modifiers is important because they may provide critical insights into HD pathogenesis and reveal key pathways that could be targeted by novel HD therapeutics. This is important since there are no disease-modifying therapies for HD, and mHtt is an unattractive small-molecule drug target. We recruited 21 HD families with varying characteristics of disease progression and age of onset and obtained medical histories, clinical records and DNA samples that were subjected to whole-genome sequencing (WGS). These WGS data describe families of 104 subjects, including HD patients and their unaffected family members. These individuals were selected based on individual clinical histories and family structures that best fit our criteria for expressing potential genetic modifiers. We are testing the hypothesis that novel, rare genetic variants contribute to HD and those genetic modifiers can be identified by WGS.

  Study phs001071

• Systems genetics in human endothelial cells identifies non-coding variants modifying enhancers, expression, and complex disease traits

  Identification of causal variants and mechanisms underlying complex disease traits in humans requires strategies to locate and fine-map functional regulatory variants in disease-relevant cell types. To discover functional regulatory variants in primary aortic endothelial cells (ECs) from humans, we collected genetic, transcriptomic, and four epigenomic phenotypes in a population of 153 human donors representing individuals of both sexes and three major ancestries. We found thousands of EC eQTLs that were not present in GTEx at all ranges of effect sizes, indicating novel functional variants not observable from tissue data. We performed several epigenetic assays across 53 donors' EC to enable molecular QTL mapping, which included chromatin accessibility, histone modification, and transcription factor binding of proteins ERG and NF-kB/p65 in two activation states. We discovered over 3000 regulatory elements where cis-regulatory variants associated to significant differences in epigenetic marks. Co-localization with our eQTLs, TF motif mutations, and 3D confirmation capture data in ECs enabled our discovery of hundreds of high-confidence functional regulatory elements that perturb endothelial molecular phenotypes. Furthermore, this set of variants is enriched at GWAS loci for Coronary Artery disease and other disease traits, with some SNPs demonstrating pleiotropic effects. Together, our dataset and analytical approach provide the genetics and vascular biology communities with specific variants affecting EC biology and serves as a proof-of-principle for how to detect functional regulatory variants in human cells with relevance to complex traits.

  Study phs002057

• Male Infertility: Genetics of Spermatogenic Failure

  Spermatogenesis is a complex biological process that requires the coordination of thousands of genes. Perhaps because of the large number of genes involved, spermatogenic failure occurs frequently and affects approximately 1% of men. While environmental and genetic factors likely contribute to this disorder, it is thought that the majority of cases have an underlying genetic basis. The two most common genetic causes are deletions on the Y chromosome and cytogenetic abnormalities (/e.g./ Klinefelter syndrome, XXY), which each account for roughly 10-15% of cases of complete spermatogenic failure. Point mutations in several other genes have also been linked to spermatogenic failure, but their collective prevalence is very low. Thus, for approximately 70% of men with spermatogenic failure, the genetic cause remains unknown. We hypothesize that a disproportionate number of these remaining genetic variants reside on the sex chromosomes because they are hemizygous in males and contain a disproportionate number of genes expressed in the testis. To identify genes that are required for spermatogenesis, we have performed capture-based targeted sequencing in 301 men with azoospermia (complete absence of sperm) and 300 fertile controls. Our sequencing is focused on coding genes and conserved, non-coding sequences of the X and Y chromosomes. In addition, we have targeted ~497 autosomal genes that display exclusive or predominant expression in the testis. The total sequence space targeted is 21.3 Mb, and we require >70% of targets reach an average coverage of 20X.

  Study phs001023

• Enrichment of Lung Microbiome with Supraglottic Taxa which is Associated with Increased Pulmonary Inflammation

  Background: The lung microbiome of healthy individuals frequently harbors oral organisms. Despite evidence that micro-aspiration is commonly associated with smoking-related lung diseases, the effects of lung microbiome enrichment with upper airway taxa on inflammation has not been studied. We hypothesize that the presence of oral microorganisms in the lung microbiome is associated with enhanced pulmonary inflammation. Methods: We sampled bronchoalveolar lavage (BAL) from the lower airways of 29 asymptomatic subjects (9 never-smokers, 14 former-smokers and 6 current-smokers). We quantified, amplified, and sequenced 16S rRNA genes from BAL samples by qPCR and 454 sequencing. Pulmonary inflammation was assessed by exhaled nitric oxide (eNO), BAL lymphocytes and neutrophils. Results: BAL had lower total 16S than supraglottic samples and higher than saline background. Bacterial communities in the lower airway clustered in two distinct groups that we designated as pneumotypes. The rRNA gene concentration and microbial community of the first pneumotype was similar to that of the saline background. The second pneumotype had higher rRNA gene concentration and higher relative abundance of supraglottic-characteristic taxa (SCT), such as Veillonella and Prevotella, and we called it pneumotypeSCT. Smoking had no effect on pneumotype allocation, alpha or beta diversity. PneumotypeSCT was associated with higher BAL lymphocyte-count (p = 0.007), BAL neutrophil-count (p = 0.034) and eNO (p = 0.022). Conclusion: A pneumotype with high relative abundance of supraglottic-characteristic taxa is associated with enhanced subclinical lung inflammation.

  Study phs000633

• Rapid Whole Genome Sequencing for Genetic Disease Diagnosis in Neonatal Intensive Care Units

  Monogenic diseases are frequent causes of neonatal morbidity and mortality, and disease presentations are often undifferentiated at birth. More than 3,500 monogenic diseases have been characterized, but clinical testing is available for only some of them and many feature clinical and genetic heterogeneity. As such, an immense unmet need exists for improved molecular diagnosis in infants. Because disease progression is extremely rapid, albeit heterogeneous, in newborns, molecular diagnoses must occur quickly to be relevant for clinical decision-making. We describe 50-hour differential diagnosis of genetic disorders by whole-genome sequencing (WGS) that features automated bioinformatic analysis and is intended to be a prototype for use in neonatal intensive care units. Retrospective 50-hour WGS identified known molecular diagnoses in two children. Prospective WGS disclosed potential molecular diagnosis of a severe GJB2-related skin disease in one neonate; BRAT1-related lethal neonatal rigidity and multifocal seizure syndrome in another infant, identified BCL9L as a novel, recessive visceral heterotaxy gene (HTX6) in a pedigree, and ruled out known candidate genes in one infants. Sequencing of parents or affected siblings expedited the identification of disease gene in prospective cases. Thus, rapid WGS can potentially broaden and foreshorten differential diagnosis, resulting in fewer empirical treatments and faster progression to genetic and prognostic counseling. Reprinted from Saunders et. al, Rapid Whole-Genome Sequencing for Genetic Disease Diagnosis in Neonatal Intensive Care Units. Sci. Transl. Med. 4, 154ra135 (2012; PMID: 23035047) with permission from AAAS.

  Study phs000564

• A super-enhancer associated with CD47 links pro-inflammatory signaling to CD47 upregulation in breast cancer

  CD47 is a cell surface molecule that inhibits phagocytosis of cells that express it by binding to its receptor, SIRPα, on macrophages and other immune cells. CD47 is expressed at different levels in normal cells, however, in cancer cells, CD47 transcript and protein expression is aberrantly increased. Here we sought to uncover the regulators of CD47 transcription, including active enhancers that increase its aberrant expression in cancer cells, in order to reveal mechanisms by which different neoplastic cells generate this dominant 'don't eat me' signal. Enhancers are genomic regions, often referred to as "switches", that can turn on or off the transcription of target genes. Recently the discovery of super-enhancers (SEs) has given more insight into the regulatory architecture of key genes that are highly expressed in a specific cell type, during a particular developmental stage or in disease. By analyzing the CD47 regulatory genomic landscape, we discovered: i) A distinct super-enhancer (SE) is associated with CD47 upregulation in breast cancer cells ii) Disruption of CD47 SEs by using the BRD4 inhibitor JQ1 robustly reduces CD47 gene expression; and iii) The TNF-NFKB1 signaling pathway is directly involved in the regulation of CD47 by interacting with a distal downstream constituent enhancer located within a CD47-associated SE specific to breast cancer. Our results describe a novel mechanism that cancer cells have evolved to drive CD47 overexpression to escape immune surveillance.

  Study phs001264

• Indexing Genes Impacted by Copy Number Variation in Developmental Disorders

  A challenge in clinical genomics is to predict whether copy number variation (CNV) affecting a gene or multiple genes will manifest as disease. Increasing recognition of gene dosage effects in neurodevelopmental disorders prompted us to develop a computational approach based on critical-exon (highly expressed in brain, highly conserved) examination for potential etiologic effects. Using a large CNV dataset, our updated analyses revealed significant (P < 1.64 x 10-15) enrichment of critical-exons within rare CNVs in cases compared to controls. Separately, we used a weighted gene co-expression network analysis (WGCNA) to construct an unbiased protein module from prenatal and adult tissues and found it significantly enriched for critical exons in prenatal (P < 1.15 x 10-50, OR = 2.11) and adult (P < 6.03 x 10-18, OR = 1.55) development. WGCNA yielded 1,206 proteins for which we prioritized the corresponding genes as likely to have a role in neurodevelopmental disorders. We compared the gene lists obtained from critical-exon and WGCNA analysis and found 438 candidate genes associated with CNVs annotated as pathogenic, or variants of uncertain significance (VOUS), from among 10,619 developmental delay cases. We identified genes containing CNVs previously considered to be VOUS, to be new candidate genes for neurodevelopmental disorders (GIT1, MVB12B and PPP1R9A) demonstrating the utility of this strategy to index the clinical effects of CNVs. Reprinted from Uddin et al, 2016, with permission from Scientific Reports.

  Study phs001154

• Genomic Profiling of Peripheral T-cell Lymphoma

  Peripheral T-cell lymphomas (PTCLs) are malignant and highly aggressive hematologic tumors arising from mature post-thymic T cells. The diagnosis of PTCL includes diverse lymphoma subgroups, altogether accounting for about 15% of all non-Hodgkin lymphomas. Despite recent advances in uncovering the mutational landscape of PTCLs, the prognoses of most PTCLs remain poor. To overcome the limits of exome sequencing and identify new genomic abnormalities in PTCL, we performed RNAseq and whole genome sequencing. From our analyses, we identified a recurrent in-frame deletion (VAV1 Δ778-786) generated by a focal deletion-driven alternative splicing mechanism as well as novel VAV1 gene fusions (VAV1-THAP4, VAV1-MYO1F, and VAV1-S100A7). These genetic lesions result in increased activation of VAV1 catalytic-dependent (MAPK, JNK) and non-catalytic-dependent (nuclear factor of activated T cells, NFAT) VAV1 effector pathways. These results support an oncogenic role for VAV1 signaling in the pathogenesis of PTCL (Abate et al., 2017). Recently, we have also identified a novel chimeric oncogene, FYN-TRAF3IP2, in the RNA-seq and whole genome sequencing analyses. FYN-TRAF3IP2 leads to aberrant NF-κβ activity downstream of T-cell receptor signaling and induces NF-κβ-driven T-cell transformation when expressed in hematopoietic progenitors in mice. This novel gene fusion supports the oncogenic driver role of NF-κβ signaling in the pathogenesis of PTCL (Kim et al., in review).

  Study phs001962

• Genomic Analysis of Pediatric Low Grade Gliomas

  Pediatric low-grade gliomas (PLGGs) are among the most common solid tumors in children but, apart from mutations or duplications in the BRAF kinase in specific subclasses, few genetic driver events are known. Diffuse PLGGs compose a set of uncommon subtypes that exhibit invasive growth and are therefore especially challenging clinically. These tumors are particularly poorly understood. We performed high-resolution copy-number analysis of 44 diffuse PLGGs to identify recurrent alterations. Diffuse PLGGs exhibited fewer such alterations than adult low-grade gliomas, but we identified several significantly recurrent events. The most significant event, 8q13.1 gains, was observed in 28% of diffuse astrocytoma WHO grade II (DA2) and resulted in partial duplication of the transcription factor MYBL1 with truncation of its C-terminal negative-regulatory domain. A similar recurrent deletion-truncation breakpoint was identified in two angiocentric gliomas in the related gene MYB on 6q23.3. Whole genome sequencing of a MYBL1-rearranged diffuse astrocytoma grade II demonstrated MYBL1 tandem duplication and few other events. Two truncated MYBL1 transcripts identified in this tumor induced anchorage-independent growth when expressed in 3T3 cells and tumor formation in nude mice. Truncated transcripts were also expressed in two additional tumors with MYBL1 partial duplication. Our results define clinically relevant molecular subclasses of diffuse PLGGs and highlight a potential role for the MYB family in the biology of low-grade gliomas. "Reprinted from www.pnas.org/cgi/doi/10.1073/pnas.1300252110 with permission from PNAS."

  Study phs000614

• A New Pipeline to Predict and Confirm Tumor Neoantigens Predicts Better Response to Checkpoint Blockade

  Mutations that drive oncogenesis in cancer can generate neoantigens that may be recognized by the immune system. Identification of these neoantigens remains challenging due to the complexity of the major histocompatibility complex (MHC) antigen and T-cell receptor interaction. Here we describe the development of a systematic approach to efficiently identify and validate immunogenic neoantigens. Whole exome sequencing of tissue from a melanoma patient was used to identify nonsynonymous mutations, followed by MHC binding prediction and identification of tumor clonal architecture. The top 18 putative class I neoantigens were selected for immunogenicity testing via a novel in vitro pipeline in HLA-A201 healthy donor blood. Naive CD8 T cells from donors were stimulated with allogeneic dendritic cells pulsed with peptide pools and then with individual peptides. The presence of antigen-specific T cells was determined via functional assays. We identified one putative neoantigen that expanded T cells specific to the mutant form of the peptide and validated this pipeline in a subset of patients with bladder tumors treated with durvalumab (N = 5). Within this cohort, the top predicted neoantigens from all patients were immunogenic in vitro. Finally, we looked at overall survival in the whole durvalumab-treated bladder cohort (N = 37) by stratifying patients by tertile measure of tumor mutation burden (TMB) or neoantigen load. Patients with higher neoantigen and TMB load tended to show better overall survival.

  Study phs002269

• Whole genome and transcriptome sequencing of lung cancer patients and cell lines at Genentech

  Version 1 Whole genome sequencing was applied to tumor and adjacent normal lung tissue in an individual non-small-cell lung cancer patient. We present an analysis of somatic changes identified throughout the tumor genome, including single-nucleotide variants, copy number variants, and large-scale chromosomal rearrangements. Over 50,000 high-confidence single-nucleotide variants were identified, revealing evidence of substantial smoking-related DNA damage as well as distinct mutational pressures within the tumor resulting in uneven distribution of somatic mutations across the genome. Version 2 Lung cancer is a highly heterogeneous disease in terms of both underlying genetic lesions and response to therapeutic treatments. We performed deep whole genome sequencing and transcriptome sequencing on 19 lung cancer cell lines and 3 lung tumor/normal pairs. Overall, our data show that cell line models exhibit similar mutation spectra to human tumor samples. Smoker and never-smoker cancer samples exhibit distinguishable patterns of mutations. A number of epigenetic regulators are frequently altered by mutations or copy number changes. A systematic survey of splice-site mutations identified over 100 splice site mutations associated with cancer specific aberrant splicing, including mutations in several known cancer-related genes. Differential usages of splice isoforms were also studied. Taken together, these data present a comprehensive genomic landscape of a large number of lung cancer samples and further demonstrate that cancer specific alternative splicing is a widespread phenomenon that has potential utility as therapeutic biomarkers.

  Study phs000299

• Tumor Mutation Burden, Expressed Neoantigen and Immune Microenvironment in Diffuse Gliomas

  A consistent correlation between tumor mutation burden (TMB) and tumor immune microenvironment has not been observed in gliomas as in other cancers. Driver germline and somatic mutations, TMB, neoantigen, and immune cell signatures were analyzed using whole exome sequencing (WES) and transcriptome sequencing of tumor and WES of matched germline DNA in a cohort of 66 glioma samples (44 IDH-mutant and 22 IDH-wildtype). Fourteen samples revealed a hypermutator phenotype (HMP). Eight pathogenic (P) or likely pathogenic (LP) germline variants were detected in 9 (19%) patients. Six of these 8 genes were DNA damage repair genes. P/LP germline variants were found in 22% of IDH-mutant gliomas and 12.5%of IDH-wildtype gliomas (p= 0.7). TMB was correlated with expressed neoantigen but showed an inverse correlation with immune score (R =−0.46, p= 0.03) in IDH-wildtype tumors and no correlation in IDH-mutant tumors. Antigen Processing and Presentation (APP) score correlates with immune score. The APP is higher in non-hypermutator phenotype (NHMP) samples versus HMP samples in the IDH-wildtype gliomas, but higher in HMP samples versus NHMP samples in the IDH-mutant gliomas. Conclusion: TMB was inversely correlated with the immune score in the IDH-wildtype gliomas, and showed no correlation in the IDH-mutant tumors. The APP was correlated with the immune score, and may be further investigated as a biomarker for response to immunotherapy in gliomas. Studies of germline variants in a large glioma cohort are warranted.

  Study phs002653

• Genetic Basis for Clinical Response to CTLA-4 Blockade in Melanoma

  Immune checkpoint inhibitors are effective cancer treatments, but molecular determinants of clinical benefit are unknown. Ipilimumab and tremelimumab are antibodies against cytotoxic T-lymphocyte antigen 4 (CTLA-4). Anti-CTLA-4 treatment prolongs overall survival in patients with melanoma. CTLA-4 blockade activates T cells and enables them to destroy tumor cells. We obtained tumor tissue from patients with melanoma who were treated with ipilimumab or tremelimumab. Whole-exome sequencing was performed on tumors and matched blood samples. Somatic mutations and candidate neoantigens generated from these mutations were characterized. Neoantigen peptides were tested for the ability to activate lymphocytes from ipilimumab-treated patients. Malignant melanoma exomes from 64 patients treated with CTLA-4 blockade were characterized with the use of massively parallel sequencing. A discovery set consisted of 11 patients who derived a long-term clinical benefit and 14 patients who derived a minimal benefit or no benefit. Mutational load was associated with the degree of clinical benefit (P = 0.01) but alone was not sufficient to predict benefit. Using genomewide somatic neoepitope analysis and patient-specific HLA typing, we identified candidate tumor neoantigens for each patient. We elucidated a neo-antigen landscape that is specifically present in tumors with a strong response to CTLA-4 blockade. We validated this signature in a second set of 39 patients with melanoma who were treated with anti-CTLA-4 antibodies. Predicted neoantigens activated T cells from the patients treated with ipilimumab.

  Study phs001041

• University of Texas PDX Development and Trial Center Grant

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

  Study phs001980

• Extent and Significance of Bacterial DNA Integrations in the Human Cancer Genome

  The integration of exogenous DNA into the human genome can cause somatic mutations associated with oncogenesis. For example, the insertion of HPV DNA into human chromosomes is the single most important event leading to tumorigenesis in cervical cancer. It is also now preventable with vaccines against HPV. In contrast to viral DNA integrations, the instances and repercussions of bacterial DNA integration into the somatic human genome are less clear. Yet bacterial DNA integrations (BDI) found to be associated with tumorigenesis could also be prevented using therapeutics like vaccines that limit exposure to the bacteria. Previous analysis of publicly available sequencing data from the Cancer Genome Atlas (TCGA) showed BDIs from Pseudomonas spp. rRNA into the 5'-UTR of four proto-oncogenes as well as in Ig (the immunoglobulin gene) of gastric cancer samples. The original samples from TCGA were not available to validate these findings. Therefore, a different cohort of gastric cancer patients was identified in order to establish the presence BDIs. Whole exome sequencing was completed on seven tumor samples and six adjacent gastric samples, as well as one intestinal metaplasia sample and one non-atrophic gastritis sample. Four of those samples were also investigated with transcriptomics and whole genome sequencing. Given that 15-20% of cancers worldwide are linked to infections, more work is needed to understand the role of infectious agents in oncogenesis and the mechanisms underlying that role.

  Study phs001936

• Small Intestine Neuroendocrine Tumors (Carcinoid Tumors)

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

  Study phs000579

• Exome Sequencing of Esophageal Adenocarcinoma

  The incidence of esophageal adenocarcinoma (EAC) has risen 600% over the last 30 years. When coupled to the five-year survival rate of only 15% for this disease, identification of new therapeutic targets for EAC is of great importance. Here, we analyze the mutational spectra identified from the whole genome sequencing of 16 EACs and complete exome sequencing of 149 EAC tumors and matched germline samples. We identify a novel mutation signature marked by high prevalence of A to C transversions at AA*(N) trinucleotides. Notably, the development of EAC is preceded by pathologic intestinal metaplasia of the lower esophagus, Barrett's esophagus, itself a response to injury from gastric-esophageal reflux disease (GERD). Thus, we hypothesize to represent these AA*(N) mutations to a novel signature of GERD-induced mutagenesis. Analysis of the exome data identified 26 genes subject to statistically significant recurrent mutation. Of these 26 genes, only TP53, CDKN2A, SMAD4, and PIK3CA had been previously implicated in EAC. Novel mutations include those targeting chromatin modifying factors and novel candidate contributors to EAC: SPG20, TLR4, ELMO1 and DOCK2. Of these, ELMO1 and DOCK2 are notably direct dimerization partners that act to activate Rac1 to enhance cellular invasiveness, thus generating new hypotheses about the potential for the Rac1 pathway to be a novel target for therapy of EAC. "Reprinted from 'Unraveling the mutational complexity of esophageal adenocarcinoma', with permission from Nature Genetics."

  Study phs000598

• Region-Specific Neural Stem Cell Lineages Revealed by Single-Cell RNA-Sequences from Human Embryonic Stem Cells

  During development of the human brain, multiple cell types with diverse regional identities are generated. Here we report a system to generate early human brain forebrain and mid/hindbrain cell types from human embryonic stem cells (hESCs), and infer and experimentally confirm a lineage tree for the generation of these types based on single-cell RNA-Seq analysis. We engineered SOX2Cit/+ and DCXCit/Y hESC lines to target progenitors and neurons throughout neural differentiation for single-cell transcriptomic profiling, then identified discrete cell types consisting of both rostral (cortical) and caudal (mid/hindbrain) identities. Direct comparison of the cell types were made to primary tissues using gene expression atlases and fetal human brain single-cell gene expression data, and this established that the cell types resembled early human brain cell types, including preplate cells. From the single-cell transcriptomic data a Bayesian algorithm generated a unified lineage tree, and predicted novel regulatory transcription factors. The lineage tree highlighted a prominent bifurcation between cortical and mid/hindbrain cell types, confirmed by clonal analysis experiments. We demonstrated that cell types from either branch could preferentially be generated by manipulation of the canonical Wnt/beta-catenin pathway. In summary, we present an experimentally validated lineage tree that encompasses multiple brain regions, and our work sheds light on the molecular regulation of region-specific neural lineages during human brain development.

  Study phs001205

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

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

  Study phs000942

• Targeted Linked-Read DNA-seq Analysis of Castration-Resistant Prostate Cancers

  Androgen receptor (AR) inhibition is standard of care for advanced prostate cancer (PC). However, its efficacy is limited by progression to castration-resistant PC (CRPC), usually due to AR re-activation via mechanisms that include AR amplification and structural rearrangement. These two classes of AR alterations often co-occur in CRPC tumors, but it is unclear whether this reflects intercellular or intracellular heterogeneity of AR. It is important to resolve the mechanisms of AR genomic rearrangements to develop new therapies and predictive biomarkers for PC. We developed a targeted linked-read DNA sequencing assay to interrogate alterations in AR and a panel of genes known to be amplified or mutated in CRPC. Linked-read DNA-sequencing libraries were generated using the 10X Genomics Chromium platform and subjected to hybridization enrichment using a custom-designed Agilent SureSelect bait panel. Post-capture DNA-seq libraries were pooled and sequenced on an Illumina HiSeq 2500 using 2x150 bp settings. This targeted linked-read DNA-seq assay was used to interrogate a set of 6 benchmarking samples harboring known AR gene rearrangements and other alterations in AR, as well a panel of 6 patient-derived xenograft CRPC tumors and 23 clinical CRPC tumors. In this cohort, we identified a subset of samples that developed complex, multiply-rearranged AR gene structures in conjunction with very high AR copy number. The data that will be made available through dbGaP are raw FASTQ files.

  Study phs003343

• Center for Sub-Cellular Genomics

  A cell is a highly complex system with distributed molecular physiologies in structured sub-cellular compartments whose interplay with the nuclear genome determine the functional characteristics of the cell. A classic example of distributed genomic processes is found in neurons. Learning and memory require modulation of individual synapses through RNA localization, localized translation, and localized metabolites, such as those from dendritic mitochondria. Dendrites of neurons integrate distributed synaptic signals into both electrical and nuclear transcriptional response. Dysfunction of these distributed genomic functions in neurons can result in a broad spectrum of neuropsychiatric diseases such as bipolar and depressive disorders, autism, among others. The Center for Sub-Cellular Genomics' (CSG) mission is to develop new technologies to enable the dissection of complex genomic interactions within a single cell. New technologies include methods to target DNA and RNA sampling at sub-cellular resolution in live cells; high-throughput isolation and microfluidics chemistry for sub-cellular fractions; sub-cellular resolution mass-spectrometry; and computational tools for analysis of sub-cellular genomics data. The new technologies developed by CSG will enable the broader genomics community to associate sub-cellular processes to cell function and dysfunction. These technologies will be utilized in the center to generate a comprehensive atlas of genomic signatures of sub-cellular compartments in the mammalian neuron system. The technologies and neuron-focused data obtained here will be relevant to neurodegenerative and neuropsychiatric diseases including schizophrenia, depression, autism spectrum disorder, traumatic brain injury, and Alzheimer's.

  Study phs002120

• Molecular Characterization of Hemimegalencephaly

  Somatic mutations, in which a fraction of the cells in the body have a deleterious mutation, are well recognized in cancer but only recently appreciated in neurological disease. Given a somatic mutation in a population of progenitor cells, all daughter cells inherit the mutation and are able to express the resultant phenotype as a function of the differentiation program. We recently identified the first de novo somatic mutations in the developing brain in the condition hemimegalencephaly (HME), a catastrophic focal epilepsy condition associated with a malformation of cerebral cortical development (MCD). HME is one of the most severe MCD syndromes, characterized by massive hamartomatous overgrowth of either of the two cerebral hemispheres. Cerebral hemispherectomy is a frequent treatment for refractory epilepsy, allowing the sampling of diseased tissue. By comparing DNA from a diseased brain with DNA from blood/saliva, we identified de novo somatic mutations in PIK3CA, AKT3 or MTOR, part of the mTOR pathway. Mutations were present in 8-40% of sequenced alleles in various brain regions sampled during surgery and in some codons known to activate proteins.• Population informationPatient samples were collected from multiple centers. The population is Caucasian and Hispanic• Molecular technologies employedWhole-exome sequencing• Principal findings of the studyWe were the first to describe mTOR-related and non-mTOR mutations causing hemimegalencephaly, and we described two-hit mutational models to explain the genetic pathogenesis of hemimegalencephaly

  Study phs002156

• Trisomy 21 Dosage Compensation Map (T21DoCoMap)

  Excess gene dosage from human chromosome 21 (HSA21), due to trisomy or translocation of HSA21 material, causes Down syndrome (DS). We derived trisomy 21 (T21) induced pluripotent stem cells (iPSCs) alongside otherwise isogenic euploid controls from mosaic DS fibroblasts (GM00260, first described in Lubiniecki et al., PMID: 225951) of the NIGMS genetic cell repository, to enable analysis of whole-chromosomal and individual gene dosage imbalance in the context of in vitro models of DS development. Virtual karyotyping (by Infinium CytoSNP-850k), and phased genotyping (by linked-read whole-genome sequencing) of the resulting iPSC lines provide a deeply-characterized genomic platform for dissection of HSA21 gene dosage. We apply this highly unique resource to allelic analysis of mRNA sequencing of iPSCs and differentiated neural lineages, and demonstrate its utility in the context of inducible HSA21 dosage compensation. To this end, we equipped one distinct HSA21 copy with a doxycycline-inducible XIST transgene, which endogenously silences one X in 46,XX females. Standard and allele-specific RNA-seq analysis confirms mono-allelic T21 silencing to be virtually complete and irreversible, further supported by DNA methylation data (Infinium methylEPIC). This full-length HSA21 map of phased single-nucleotide and insertion/deletion variants enables allele-specific applications for gene-level dissection T21 dosage effects. Because our XIST remains inducible in post-mitotic T21 neurons and astrocytes, we perform single-nuclei RNA-seq to demonstrate XIST efficiently represses genes even after terminal differentiation.

  Study phs002397

• Machine learning to detect the SINEs of cancer

  We previously described an approach called RealSeqS to evaluate aneuploidy in plasma cell-free DNA (cfDNA) through the amplification of ~350,000 repeated elements with a single primer. We hypothesized that an unbiased evaluation of the large amount of sequencing data obtained with RealSeqS might reveal other differences between plasma samples from patients with and without cancer. This hypothesis was tested through the development of a novel machine-learning approach called Alu Profile Learning Using Sequencing (A-PLUS) and its application to samples from 5108 individuals, 2037 with cancer and the remainder without cancer. Samples from cancer patients and controls were pre-specified into four cohorts used for: 1) model training, 2) analyte integration and threshold determination, 3) validation, and 4) reproducibility. A-PLUS alone provided a sensitivity of 40.5% across 11 different cancer types in the Validation Cohort, at a specificity of 98.5%. Combining A-PLUS with aneuploidy and 8 common protein biomarkers detected 51% of 1167 cancers at 98.9% specificity. We found that part of the power of A-PLUS could be ascribed to a single feature – the global reduction of AluS sub-family elements in the circulating DNA of cancer patients. We confirmed this reduction through the analysis of another independent dataset obtained with a very different approach (whole genome sequencing). The evaluation of Alu elements therefore has the potential to enhance the performance of several methods designed for the earlier detection of cancer.

  Study EGAS00001007169

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

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

  Study phs002999

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

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

  Study EGAS50000000605

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

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

  Study EGAS50000000307

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

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

  Study EGAS50000000135

• Enhanced Adjuvanticity of a Personal Neoantigen Vaccine Generates Potent Neoantigen-Specific Immunity

  Personalized neoantigen-targeting vaccines can induce durable and diverse T cell responses and have been associated with improved outcomes. To augment immunogenicity, we tested a synthetic long peptide vaccine formulated with the adjuvants Montanide and poly-ICLC in 10 patients with melanoma. The vaccines, co-administered locally with ipilimumab and systemically with nivolumab, generated de novo ex vivo T cell responses against the majority of immunizing neoepitopes in all 9 fully vaccinated patients, and ex vivo CD8+ T cell responses in 6 of 9. Hundreds of circulating and intratumoral T cell receptor (TCR) clonotypes emerged following vaccination that were distinct from those arising after PD-1 inhibition. By linking vaccine-specificity and cellular phenotype of TCRs in post-treatment tumors at single cell resolution, we demonstrate remodeling of the intratumoral T cell repertoire following vaccination. These observations demonstrate the impact of multi-pronged immune adjuvanticity on the T cell repertoire in the context of neoantigen-targeting vaccines. We are submitting WES data from whole blood samples and tumor FFPE samples, and RNA-seq data from tumor FFPE samples. We are also submitting RNA-seq data from purified tumor samples. scRNA-seq data will also be submitted for skin samples (pre and post neoantigen vaccine dose #3), and tumor samples before starting treatment, before starting neoantigen vaccine Dose #1, and before neoantigen vaccine dose #4. Bulk TCR-seq data will be submitted for blood samples.

  Study phs003919

• HuBMAP: A Spatially Resolved Molecular Atlas of Human Endothelium

  The human circulatory system plays a vital role in the proper functioning of every organ in the human body. Because vascular related diseases, such as atherosclerosis, do not occur homogeneously throughout the vascular tree, differences in the endothelium of various organs must exist. There is even evidence of heterogeneity of adjacent endothelial cells. We are studying normal endothelium throughout the body and across individuals of different ages to characterize these differences. The first step will be to procure normal tissue samples from a wide range of organs. From disaggregated samples of these tissues, we will perform single cell combinatorial indexing (sci-) to generate chromatin accessibility (sci-ATAC-seq) and transcriptome (sci-RNA-seq) data from millions of single cells across hundreds of samples.  Molecular heterogeneity might correlate with multiple facets of the anatomical distribution of endothelium, including vascular vs. lymphatic, arteries vs. veins, capillaries vs. larger vessels, differences between organs, organ zonation, cell-to-cell differences, etc. Based on heterogeneous markers defined from deep profiling of disaggregated cells, we will employ seqFISH to spatially relate the distinct endothelial subpopulations arising from the single nuclei data to histologic preparations. For example, this will allow us to identify subpopulations with particular markers as arterial, capillary, venous, or lymphatic. The end product will include multiparameter images in which the complete gene expression and chromatin accessibility profile can be interactively visualized for every cell.

  Study phs002267

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

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

  Study EGAS50000000523

• CRISPR screening in human trophoblast stem cells reveals both shared and distinct aspects of human and mouse placental development

  The placenta serves as the interface between the mother and fetus, facilitating the exchange of gases and nutrients between their separate blood circulation systems. Trophoblasts in the placenta play a central role in this process. Our current understanding of mammalian trophoblast development relies largely on mouse models. However, given the diversification of mammalian placentas, findings from the mouse placenta cannot be readily extrapolated to other mammalian species, including humans. To fill this knowledge gap, we performed CRISPR knockout (KO) screening in human trophoblast stem cells (hTSCs). We targeted genes essential for mouse placental development and identified more than 100 genes as critical regulators in both human hTSCs and mouse placentas. Among them, we further characterized in detail two transcription factors, DLX3 and GCM1, and revealed their essential roles in hTSC differentiation. Moreover, a gene function-based comparison between human and mouse trophoblast subtypes suggests that their relationship may differ significantly from previous assumptions based on tissue localization or cellular function. Notably, our data reveal that hTSCs may not be analogous to mouse TSCs or the extraembryonic ectoderm (ExE) in which in vivo TSCs reside. Instead, hTSCs may be analogous to progenitor cells in the mouse ectoplacental cone and chorion. This finding is consistent with the absence of ExE-like structures during human placental development. Our data not only deepen our understanding of human trophoblast development but also facilitate cross-species comparison of mammalian placentas.

  Study JGAS000659

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

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

  Study JGAS000646

• Homozygous inactivation of CHEK2 is linked to a familial case of multiple primary lung cancer with accompanying cancers in other organs

  In clinical practice, there are a number of cancer patients with clear family histories, but not belonging to existing familial malignancies. Recent advances in genomic technologies have enhanced the possibility of identifying causative genes in such cases. Two siblings, an elder sister and a younger brother, were found to have multiple primary lung cancers at the age of 60. The former subsequently developed breast cancer, and had a history of myoma uteri. The latter had initially developed prostate cancer at the age of 59, and had a history of colon cancer. SNP genotyping revealed that approximately 10% of the genomes were homozygous in both patients. Exome sequencing revealed non-synonymous mutations in four genes in the runs of homozygosity: CHEK2, FCGRT, INPP5J, and SFI1. Evolutionary conservation of primary protein structures suggested the functional importance of the CHEK2 mutation, R474C. This mutation altered the tertiary structure of CHK2 by disrupting the salt bridge between R474 and E394. No such structural changes were observed with the other three mutated genes. Subsequent cell-based transfection analysis revealed that CHK2 R474C was unstable and scarcely activated. We concluded that CHEK2 is a causative gene in this case of familial cancer. Although homozygous inactivation of CHEK2 in mice led to cancers in multiple organs, accumulation of additional human cases is needed to establish its pathogenic role in humans.

  Study JGAS000057

• DNA methylation at HBV integrants and flanking host genomes

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

  Study JGAS000015

• Tumor-reactive heterotypic CD8 T cell clusters from clinical samples

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

  Study EGAS50000000785

• Noninvasive prenatal molecular karyotyping from maternal plasma

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

  Study EGAS00001000439

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

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

  Study EGAS00001005680

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

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

  Study EGAS00001005737

• Evolution of the African pygmy phenotype

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

  Study EGAS00001000908

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

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

  Study EGAS00001002326

• Open_Targets_020_Epigenomes_of_Cell_Lines

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

  Study EGAS00001003138

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

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

  Study EGAS00001001682

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

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

  Study EGAS00001003928

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

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

  Study EGAS00001001846

• Persistent STAG2 mutation in recurrent pediatric glioblastoma

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

  Study EGAS00001004340

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

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

  Study EGAS00001000476

• Massively parallel nanowell-based single-cell gene expression profiling

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

  Study EGAS00001002320

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

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

  Study EGAS00001002423

• Identification_of_synthetic_lethal_genes_by_CRISPR_Cas9_library

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

  Study EGAS00001002117

• The evolutionary history of human colitis-associated colorectal cancer

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

  Study EGAS00001003028

• Whole genome sequencing with linked reads of pediatric glioblastoma samples

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

  Study EGAS00001003432

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

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

  Study EGAS00001003547

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

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

  Study EGAS00001003661

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

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

  Study EGAS00001004282

• UK10K NEURO ASD GALLAGHER

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

  Study EGAS00001000112

• Identification of a Fumarate-Hydratase Deficient-like RCC Subtype with Convergent Pathway Alterations

  The development of molecularly targeted therapies builds upon the identification of tumor clades with shared molecular abnormalities. Hereditary leiomyomatosis and renal cell carcinoma (HLRCC) is a familial syndrome resulting from germline mutation of the two-hit tumor suppressor gene, Fumarate Hydratase (FH). HLRCC renal tumors present often in young adults, exhibit unique histological features and are particularly aggressive. HLRCC is uncommon and poorly understood, and whether there are sporadic renal tumors with similar characteristics, is unclear. Notably, FH can also be inactivated somatically in RCC, and should be included among FH-deficient (FHD) kidney tumors. We identified a tumor clade without FH mutations but with similar morphological and molecular characteristics as well as aggressive clinical behavior, which we name FHDL (FHD-like). We show that FHD/FHDL tumors represent a unique subtype of CIMP (CpG Island Methylator Phenotype) with the highest levels of overall genome methylation and CDKN2A silencing. FHD/FHDL tumors are characterized by convergent activation of the Hippo pathway as determined by frequent and mutually exclusive mutations in NF2 and PTPN14. AKR1B10, an NRF2 target, was upregulated in FHD/FHDL tumors as well as GATA3, which suggests a shared origin in the renal medulla and is consistent with topological studies. There is loss of 5-hydroxymethylcytosine in both FHD and FHDL tumors. These data define a group of particularly aggressive renal cancers with shared molecular features and potential targets for therapeutic intervention.

  Study EGAS00001002646

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

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

  Study EGAS00001002678

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

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

  Study EGAS00001004301

• Genomic Alterations in Gingivo-buccal Cancer: ICGC-India Project_YR01

  As a part of the ICGC, India has undertaken genomic studies on gingivobuccal cancer of the oral cavity, which is the most prevalent form of cancer among men in India. There are various known environmental (life-style) correlates of this cancer, the most important of which are tobacco chewing and HPV infection. Paired DNA samples – isolated from the tumor tissue and from the blood of fifty patients – have been analyzed to catalog germline and somatic mutations. Association between each observed genomic alteration and exposure to environmental risk factors is being explored.Detailed clinical characterization of the patients, collection of data on demographic and environmental exposures, and isolation of DNA samples from blood and tumor tissues collected from each patient are being done at the Advanced Centre for Research, Treatment and Education on Cancer, Mumbai. Exome capture and deep resequencing are being performed at the National Institute of Biomedical Genomics, Kalyani. We have analyzed the exomes of paired blood and tumor DNA samples using two orthogonal platforms – exome capture by Roche-Nimblegen SeqCap Ez kit followed by sequencing on Roche-454 GS FLX, and exome capture by Illumina TrueSeq Exome Enrichment kit followed by sequencing on Illumina HiSeq-2000. Each exome is being sequenced at a mean depth of 30x. In addition, we have used Illumina Omni Quad SNP-chips to generate genotype data at 1.14 million genomic markers. The average concordance between SNP-chip and sequence data is 99.8%.

  Study EGAS00001000249

• Whole Tumor spatial heterogeneity in metastatic melanoma

  Heterogeneous inter-tumoral responses, sustained periods of apparent clinical benefit despite lack of objective response to various therapies, and even spontaneous remission are well known within a subpopulation of advanced melanoma patients. The molecular and cellular dynamics facilitating long-term survival remain poorly defined, particularly in the current era ofexposure to multiple potentially active therapies. We studied an exceptional case of long-term survival in a patient with non-responding metastatic melanoma in order to characterize the clonal and microenvironmental factors.We performed immunogenomic analysis of 3 metachronous tumors, and of 67 intratumor sub-regions of one metastatic lesion, usingwhole exome sequencing, RNA-sequencing, immunohistochemistry, and T cell receptor sequencing. Longitudinal analysesidentified mutations in severalknownresistance genes affecting distinctmetastases. Intra-tumoral immune signatures revealed spatially-distinct pockets of immune activation and suppression throughout the metastasis. A specific T cell VBeta CDR3 rearrangement was identified as both dominant and recurrent not only across multiple spatial points within a single tumor mass, but also across metachronous tumors spanning the patient’s disease course. Immunophenotyping of this T cell population with single-cell sequencingsuggestedrepeatedT-cell priming events leading to the persistence of both activated and exhausted T cells bearing this same TCR-Beta at any given time.Our findings highlight genomic and immune heterogeneityin metastatic disease, extending to the intra-tumoral leveland provide a rationale for integrative spatiotemporal analysis of metastatic lesions in the context of immune checkpoint blockade.

  Study EGAS00001003292

• Exploiting evolutionary steering in cancer therapy

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

  Study EGAS00001003200

• Probabilistic cell type assignment of single-cell transcriptomic data reveals spatiotemporal microenvironment dynamics in human cancers

  Single-cell RNA sequencing (scRNA-seq) has transformed biomedical research, enabling decomposition of complex tissues into disaggregated, functionally distinct cell types. For many applications, investigators wish to identify cell types with known marker genes. Typically, such cell type assignments are performed through unsupervised clustering followed by manual annotation based on these marker genes, or via "mapping" procedures to existing data. However, the manual interpretation required in the former case scales poorly to large datasets, which are also often prone to batch effects, while existing data for purified cell types must be available for the latter. Furthermore, unsupervised clustering can be error-prone, leading to under- and over- clustering of the cell types of interest. To overcome these issues we present CellAssign, a probabilistic model that leverages prior knowledge of cell type marker genes to annotate scRNA-seq data into pre-defined and de novo cell types. CellAssign automates the process of assigning cells in a highly scalable manner across large datasets while simultaneously controlling for batch and patient effects. We demonstrate the analytical advantages of CellAssign through extensive simulations and exemplify real-world utility to profile the spatial dynamics of high-grade serous ovarian cancer and the temporal dynamics of follicular lymphoma. Our analysis reveals subclonal malignant phenotypes and points towards an evolutionary interplay between immune and cancer cell populations with cancer cells escaping immune recognition.

  Study EGAS00001003452

• Long-read trio sequencing of unsolved patients with intellectual disability

  Long-read sequencing (LRS) has the potential to comprehensively identify all medically-relevant genome variation, including variation commonly missed by short read sequencing (SRS) approaches. To determine this potential, we performed LRS up to 40x coverage using the Pacific Biosciences Sequel instrument for 5 patient-parent trios with intellectual disability, whose aetiology remained unresolved after SRS exomes and genomes. A first systematic assessment of LRS coverage showed that ~35Mb of the genome was only accessible by LRS and not SRS. Structural variant (SV) calling yielded on average 28,292 SV calls per individual, totalling 12.9 MB of sequence. Trio-based analyses showed concordance for up to 95% of these SV calls across the genome, and 80% of the LRS SV calls were not identified by SRS. De novo mutation analysis did not identify any de novo SVs, confirming these are rare events. Because of high sequence coverage, we were also able to call single nucleotide substitutions. On average, we identified 3 million substitutions per genome, with a Mendelian inheritance concordance of up to 97%. Of these, approximately 100,000 were located in the ~35Mb of the genome that was only captured by LRS. Moreover, these variants affected the coding sequence of 64 genes, including 32 known to cause Mendelian disorders. Although no clinical diagnoses were made in these 5 patients, our data show the potential added value of LRS compared to SRS for identifying medically relevant genome variation.

  Study EGAS00001004319

• The Medical Genome Reference Bank: a whole genome data resource of 4,000 healthy elderly individuals.

  Allele frequency data from human reference populations is of increasing value for filtering and assignment of pathogenicity to genetic variants. Aged and healthy populations are more likely to be selectively depleted of pathogenic alleles, and therefore particularly suitable as a reference populations for the major diseases of clinical and public health importance. However, reference studies of the healthy elderly have remained under represented in human genetics. We have developed the Medical Genome Reference Bank (MGRB), a large scale comprehensive whole genome dataset of confirmed healthy elderly individuals, to provide a publicly accessible resource for health related research, and for clinical genetics. It also represents a useful resource for studying the genetics of healthy aging. The MGRB comprises 4,000 healthy, older individuals with no reported history of cancer, cardiovascular disease or dementia, recruited from two Australian community based cohorts. DNA derived from blood samples will be subject to whole genome sequencing. The MGRB will measure genome wide genetic variation in 4,000 individuals, mostly of European decent, aged 60 to 95 years (mean age 75 years). The MGRB has committed to a policy of data sharing, employing a hierarchical data management system to maintain participant privacy and confidentiality, whilst maximizing research and clinical usage of the database. The MGRB will represent a dataset of international significance, broadly accessible to the clinical and genetic research community.

  Study EGAS00001003511

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

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

  Study EGAS00001001146

• Initial whole genome sequencing of plasma cell neoplasms in First Responders exposed to the World Trade Center attack of September 11, 2001

  The World Trade Center (WTC) attack of September 11, 2001 created an unprecedented environmental exposure to known and suspected carcinogens. High incidence of multiple myeloma (MM) and precursor conditions has been reported among first responders to the WTC disaster. To expand on our prior screening studies, and to characterize the genomic impact of the exposure to known and potential carcinogens in the WTC debris, we were motivated to perform whole genome sequencing (WGS) of WTC first responders and recovery workers who developed a plasma cell disorder after the attack. PATIENTS AND METHODS: We performed WGS of 9 CD138-positive bone marrow mononuclear samples from patients who were diagnosed with plasma cell disorders after the WTC disaster. RESULTS: No significant differences were observed in comparing the post-WTC driver and mutational signatures landscape with 110 previously published WGS from 56 patients with MM and the CoMMpass WGS cohort (n=752). Leveraging constant activity of the single base substitution mutational signatures 1 and 5 over time, we estimated that tumor-initiating chromosomal gains were windowed to both pre- and post-WTC exposure. CONCLUSIONS: Although limitations in sample size preclude any definitive conclusions, our findings suggest that the observed increased incidence of plasma cell neoplasms in this population is due to complex and heterogeneous effects of the WTC exposure that may have initiated or contributed to progression of malignancy.

  Study EGAS00001004467

• Targets of MEK inhibition in DIPG

  The survival of children with diffuse intrinsic pontine glioma (DIPG) remains dismal, with new treatments desperately needed. In the era of precision medicine, targeted therapies represent an exciting treatment opportunity, yet resistance can rapidly emerge, playing an important role in treatment failure. In a prospective biopsy-stratified clinical trial, we combined detailed molecular profiling (methylation BeadArray, exome, RNAseq, phospho-proteomics) linked to drug screening in newly-established patient-derived models of DIPG in vitro and in vivo. We identified a high degree of in vitro sensitivity to the MEK inhibitor trametinib in samples which harboured genetic alterations targeting the MAPK pathway, including the non-canonical BRAF_G469V mutation, and those affecting PIK3R1 and NF1. However, treatment of PDX models and of a patient with trametinib at relapse failed to elicit a significant response. We generated trametinib-resistant clones in the BRAF_G469V model through continuous drug exposure, and identified acquired mutations in MEK1/2 (MEK1_K57N, MEK1_I141S and MEK2_I115N) with sustained pathway up-regulation. These cells showed the hallmarks of mesenchymal transition, and expression signatures overlapping with inherently trametinib-insensitive primary patient-derived cells that predicted an observed sensitivity to dasatinib. Combinations of trametinib with dasatinib and the downstream ERK inhibitor ulixertinib showed highly synergistic effects in vitro. These data highlight the MAPK pathway as a therapeutic target in DIPG, and show the importance of parallel resistance modelling and rational combinatorial treatments likely to be required for meaningful clinical translation.

  Study EGAS00001004495

• The clonal and mutational evolution spectrum of primary triple negative breast cancers

  Primary triple negative breast cancers (TNBC) represent approximately 16% of all breast cancers and are a tumour type defined by exclusion,for which comprehensive landscapes of somatic mutation have not been determined. Here we show in 104 early TNBC cases, that at the time of diagnosis these cancers exhibit a wide and continuous spectrum of genomic evolution, with some exhibiting only a handful of somatic aberrations in a few pathways, whereas others contain hundreds of somatic events and multiple pathways implicated. Integration with matched whole transcriptome sequence data revealed that only 36% of mutations are expressed. By examining single nucleotide variant (SNV) allelic abundance derived from deep re-sequencing (median > 20,000 fold) measurements in 2414 somatic mutations, we determine for the first time in an epithelial tumour, the relative abundance of clonal genotypes among cases in the population. We show that TNBC vary widely and continuously in their clonal frequencies at the time of diagnosis, with basal subtype TNBC exhibiting more variation than non-basal TNBC. Although p53 and somatic mutations appear clonally dominant compared with other pathways, in some tumours their clonal frequencies are incompatible with founder status. Mutations in cytoskeletal and cell shape/motility proteins occurred at lower clonal frequencies, suggesting they occurred later during tumour progression. Taken together our results show that future attempts to dissect the biology and therapeutic responses of TNBC will require the determination of individual tumour clonal genotypes.

  Study EGAS00001000132