Diagnostic evaluation of paediatric autoimmune lymphoproliferative immunodeficiencies (ALPID): a prospective cohort study

BACKGROUND

Lymphoproliferation and autoimmune cytopenias characterise autoimmune lymphoproliferative syndrome. Other conditions sharing these manifestations have been termed autoimmune lymphoproliferative syndrome-like diseases, although they are frequently more severe. The aim of this study was to define the genetic, clinical, and immunological features of these disorders to improve their diagnostic classification.

METHODS

In this prospective cohort study, patients were referred to the Center for Chronic Immunodeficiency in Freiburg, Germany, between Jan 1, 2008 and March 5, 2022. We enrolled patients younger than 18 years with lymphoproliferation and autoimmune cytopenia, lymphoproliferation and at least one additional sign of an inborn error of immunity (SoIEI), bilineage autoimmune cytopenia, or autoimmune cytopenia and at least one additional SoIEI. Autoimmune lymphoproliferative syndrome biomarkers were determined in all patients. Sanger sequencing followed by in-depth genetic studies were recommended for patients with biomarkers indicative of autoimmune lymphoproliferative syndrome, while IEI panels, exome sequencing, or genome sequencing were recommended for patients without such biomarkers. Genetic analyses were done as decided by the treating physician. The study was registered on the German Clinical Trials Register, DRKS00011383, and is ongoing.

FINDINGS

We recruited 431 children referred for autoimmune lymphoproliferative syndrome evaluation, of whom 236 (55%) were included on the basis of lymphoproliferation and autoimmune cytopenia, 148 (34%) on the basis of lymphoproliferation and another SoIEI, 33 (8%) on the basis of autoimmune bicytopenia, and 14 (3%) on the basis of autoimmune cytopenia and another SoIEI. Median age at diagnostic evaluation was 9·8 years (IQR 5·5-13·8), and the cohort comprised 279 (65%) boys and 152 (35%) girls. After biomarker and genetic assessments, autoimmune lymphoproliferative syndrome was diagnosed in 71 (16%) patients. Among the remaining 360 patients, 54 (15%) had mostly autosomal-dominant autoimmune lymphoproliferative immunodeficiencies (AD-ALPID), most commonly affecting JAK-STAT (26 patients), CTLA4-LRBA (14), PI3K (six), RAS (five), or NFκB (three) signalling. 19 (5%) patients had other IEIs, 17 (5%) had non-IEI diagnoses, 79 (22%) were unresolved despite extended genetics (ALPID-U), and 191 (53%) had insufficient genetic workup for diagnosis. 16 (10%) of 161 patients with a final diagnosis had somatic mutations. Alternative classification of patients fulfilling common variable immunodeficiency or Evans syndrome criteria did not increase the proportion of genetic diagnoses.

INTERPRETATION

The ALPID phenotype defined in this study is enriched for patients with genetic diseases treatable with targeted therapies. The term ALPID might be useful to focus diagnostic and therapeutic efforts by triggering extended genetic analysis and consideration of targeted therapies, including in some children currently classified as having common variable immunodeficiency or Evans syndrome.

FUNDING

Deutsche Forschungsgemeinschaft under Germany's Excellence Strategy.

TRANSLATION

For the German translation of the abstract see Supplementary Materials section.

GIGYF1 disruption associates with autism and impaired IGF-1R signaling

Autism spectrum disorder (ASD) represents a group of neurodevelopmental phenotypes with a strong genetic component. An excess of likely gene-disruptive (LGD) mutations in GIGYF1 was implicated in ASD. Here, we report that GIGYF1 is the second-most mutated gene among known ASD high–confidence risk genes. We investigated the inheritance of 46 GIGYF1 LGD variants, including the highly recurrent mutation c.333del:p.L111Rfs*234. Inherited GIGYF1 heterozygous LGD variants were 1.8 times more common than de novo mutations. Among individuals with ASD, cognitive impairments were less likely in those with GIGYF1 LGD variants relative to those with other high-confidence gene mutations. Using a Gigyf1 conditional KO mouse model, we showed that haploinsufficiency in the developing brain led to social impairments without significant cognitive impairments. In contrast, homozygous mice showed more severe social disability as well as cognitive impairments. Gigyf1 deficiency in mice led to a reduction in the number of upper-layer cortical neurons, accompanied by a decrease in proliferation and increase in differentiation of neural progenitor cells. We showed that GIGYF1 regulated the recycling of IGF-1R to the cell surface. KO of GIGYF1 led to a decreased level of IGF-1R on the cell surface, disrupting the IGF-1R/ERK signaling pathway. In summary, our findings show that GIGYF1 is a regulator of IGF-1R recycling. Haploinsufficiency of GIGYF1 was associated with autistic behavior, likely through interference with IGF-1R/ERK signaling pathway.

Early childhood epilepsies: epidemiology, classification, aetiology, and socio-economic determinants

Epilepsies of early childhood are frequently resistant to therapy and often associated with cognitive and behavioural comorbidity. Aetiology focused precision medicine, notably gene-based therapies, may prevent seizures and comorbidities. Epidemiological data utilizing modern diagnostic techniques including whole genome sequencing and neuroimaging can inform diagnostic strategies and therapeutic trials. We present a 3-year, multicentre prospective cohort study, involving all children under 3 years of age in Scotland presenting with epilepsies. We used two independent sources for case identification: clinical reporting and EEG record review. Capture-recapture methodology was then used to improve the accuracy of incidence estimates. Socio-demographic and clinical details were obtained at presentation, and 24 months later. Children were extensively investigated for aetiology. Whole genome sequencing was offered for all patients with drug-resistant epilepsy for whom no aetiology could yet be identified. Multivariate logistic regression modelling was used to determine associations between clinical features, aetiology, and outcome. Three hundred and ninety children were recruited over 3 years. The adjusted incidence of epilepsies presenting in the first 3 years of life was 239 per 100 000 live births [95% confidence interval (CI) 216-263]. There was a socio-economic gradient to incidence, with a significantly higher incidence in the most deprived quintile (301 per 100 000 live births, 95% CI 251-357) compared with the least deprived quintile (182 per 100 000 live births, 95% CI 139-233), χ2 odds ratio = 1.7 (95% CI 1.3-2.2). The relationship between deprivation and incidence was only observed in the group without identified aetiology, suggesting that populations living in higher deprivation areas have greater multifactorial risk for epilepsy. Aetiology was determined in 54% of children, and epilepsy syndrome was classified in 54%. Thirty-one per cent had an identified genetic cause for their epilepsy. We present novel data on the aetiological spectrum of the most commonly presenting epilepsies of early childhood. Twenty-four months after presentation, 36% of children had drug-resistant epilepsy (DRE), and 49% had global developmental delay (GDD). Identification of an aetiology was the strongest determinant of both DRE and GDD. Aetiology was determined in 82% of those with DRE, and 75% of those with GDD. In young children with epilepsy, genetic testing should be prioritized as it has the highest yield of any investigation and is most likely to inform precision therapy and prognosis. Epilepsies in early childhood are 30% more common than previously reported. Epilepsies of undetermined aetiology present more frequently in deprived communities. This likely reflects increased multifactorial risk within these populations.

Towards realizing the vision of precision medicine: AI based prediction of clinical drug response

Accurate and individualized prediction of response to therapies is central to precision medicine. However, because of the generally complex and multifaceted nature of clinical drug response, realizing this vision is highly challenging, requiring integrating different data types from the same individual into one prediction model. We used the anti-epileptic drug brivaracetam as a case study and combine a hybrid data/knowledge-driven feature extraction with machine learning to systematically integrate clinical and genetic data from a clinical discovery dataset (n = 235 patients). We constructed a model that successfully predicts clinical drug response [area under the curve (AUC) = 0.76] and show that even with limited sample size, integrating high-dimensional genetics data with clinical data can inform drug response prediction. After further validation on data collected from an independently conducted clinical study (AUC = 0.75), we extensively explore our model to gain insights into the determinants of drug response, and identify various clinical and genetic characteristics predisposing to poor response. Finally, we assess the potential impact of our model on clinical trial design and demonstrate that, by enriching for probable responders, significant reductions in clinical study sizes may be achieved. To our knowledge, our model represents the first retrospectively validated machine learning model linking drug mechanism of action and the genetic, clinical and demographic background in epilepsy patients to clinical drug response. Hence, it provides a blueprint for how machine learning-based multimodal data integration can act as a driver in achieving the goals of precision medicine in fields such as neurology.

Early-onset autoimmunity associated with SOCS1 haploinsufficiency

Autoimmunity can occur when a checkpoint of self-tolerance fails. The study of familial autoimmune diseases can reveal pathophysiological mechanisms involved in more common autoimmune diseases. Here, by whole-exome/genome sequencing we identify heterozygous, autosomal-dominant, germline loss-of-function mutations in the SOCS1 gene in ten patients from five unrelated families with early onset autoimmune manifestations. The intracellular protein SOCS1 is known to downregulate cytokine signaling by inhibiting the JAK-STAT pathway. Accordingly, patient-derived lymphocytes exhibit increased STAT activation in vitro in response to interferon-γ, IL-2 and IL-4 that is reverted by the JAK1/JAK2 inhibitor ruxolitinib. This effect is associated with a series of in vitro and in vivo immune abnormalities consistent with lymphocyte hyperactivity. Hence, SOCS1 haploinsufficiency causes a dominantly inherited predisposition to early onset autoimmune diseases related to cytokine hypersensitivity of immune cells.

SSBP1 mutations cause mtDNA depletion underlying a complex optic atrophy disorder

Inherited optic neuropathies include complex phenotypes, mostly driven by mitochondrial dysfunction. We report an optic atrophy spectrum disorder, including retinal macular dystrophy and kidney insufficiency leading to transplantation, associated with mitochondrial DNA (mtDNA) depletion without accumulation of multiple deletions. By whole-exome sequencing, we identified mutations affecting the mitochondrial single-strand binding protein (SSBP1) in 4 families with dominant and 1 with recessive inheritance. We show that SSBP1 mutations in patient-derived fibroblasts variably affect the amount of SSBP1 protein and alter multimer formation, but not the binding to ssDNA. SSBP1 mutations impaired mtDNA, nucleoids, and 7S-DNA amounts as well as mtDNA replication, affecting replisome machinery. The variable mtDNA depletion in cells was reflected in severity of mitochondrial dysfunction, including respiratory efficiency, OXPHOS subunits, and complex amount and assembly. mtDNA depletion and cytochrome c oxidase-negative cells were found ex vivo in biopsies of affected tissues, such as kidney and skeletal muscle. Reduced efficiency of mtDNA replication was also reproduced in vitro, confirming the pathogenic mechanism. Furthermore, ssbp1 suppression in zebrafish induced signs of nephropathy and reduced optic nerve size, the latter phenotype complemented by WT mRNA but not by SSBP1 mutant transcripts. This previously unrecognized disease of mtDNA maintenance implicates SSBP1 mutations as a cause of human pathology.

Pan-cancer analysis of whole genomes

Cancer is driven by genetic change, and the advent of massively parallel sequencing has enabled systematic documentation of this variation at the whole-genome scale1-3. Here we report the integrative analysis of 2,658 whole-cancer genomes and their matching normal tissues across 38 tumour types from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). We describe the generation of the PCAWG resource, facilitated by international data sharing using compute clouds. On average, cancer genomes contained 4-5 driver mutations when combining coding and non-coding genomic elements; however, in around 5% of cases no drivers were identified, suggesting that cancer driver discovery is not yet complete. Chromothripsis, in which many clustered structural variants arise in a single catastrophic event, is frequently an early event in tumour evolution; in acral melanoma, for example, these events precede most somatic point mutations and affect several cancer-associated genes simultaneously. Cancers with abnormal telomere maintenance often originate from tissues with low replicative activity and show several mechanisms of preventing telomere attrition to critical levels. Common and rare germline variants affect patterns of somatic mutation, including point mutations, structural variants and somatic retrotransposition. A collection of papers from the PCAWG Consortium describes non-coding mutations that drive cancer beyond those in the TERT promoter4; identifies new signatures of mutational processes that cause base substitutions, small insertions and deletions and structural variation5,6; analyses timings and patterns of tumour evolution7; describes the diverse transcriptional consequences of somatic mutation on splicing, expression levels, fusion genes and promoter activity8,9; and evaluates a range of more-specialized features of cancer genomes8,10-18.

Partial Loss of USP9X Function Leads to a Male Neurodevelopmental and Behavioral Disorder Converging on Transforming Growth Factor β Signaling

BACKGROUND

The X-chromosome gene USP9X encodes a deubiquitylating enzyme that has been associated with neurodevelopmental disorders primarily in female subjects. USP9X escapes X inactivation, and in female subjects de novo heterozygous copy number loss or truncating mutations cause haploinsufficiency culminating in a recognizable syndrome with intellectual disability and signature brain and congenital abnormalities. In contrast, the involvement of USP9X in male neurodevelopmental disorders remains tentative.

METHODS

We used clinically recommended guidelines to collect and interrogate the pathogenicity of 44 USP9X variants associated with neurodevelopmental disorders in males. Functional studies in patient-derived cell lines and mice were used to determine mechanisms of pathology.

RESULTS

Twelve missense variants showed strong evidence of pathogenicity. We define a characteristic phenotype of the central nervous system (white matter disturbances, thin corpus callosum, and widened ventricles); global delay with significant alteration of speech, language, and behavior; hypotonia; joint hypermobility; visual system defects; and other common congenital and dysmorphic features. Comparison of in silico and phenotypical features align additional variants of unknown significance with likely pathogenicity. In support of partial loss-of-function mechanisms, using patient-derived cell lines, we show loss of only specific USP9X substrates that regulate neurodevelopmental signaling pathways and a united defect in transforming growth factor β signaling. In addition, we find correlates of the male phenotype in Usp9x brain-specific knockout mice, and further resolve loss of hippocampal-dependent learning and memory.

CONCLUSIONS

Our data demonstrate the involvement of USP9X variants in a distinctive neurodevelopmental and behavioral syndrome in male subjects and identify plausible mechanisms of pathogenesis centered on disrupted transforming growth factor β signaling and hippocampal function.

Whole-Genome and Epigenomic Landscapes of Etiologically Distinct Subtypes of Cholangiocarcinoma

Cholangiocarcinoma (CCA) is a hepatobiliary malignancy exhibiting high incidence in countries with endemic liver-fluke infection. We analyzed 489 CCAs from 10 countries, combining whole-genome (71 cases), targeted/exome, copy-number, gene expression, and DNA methylation information. Integrative clustering defined 4 CCA clusters-fluke-positive CCAs (clusters 1/2) are enriched in ERBB2 amplifications and TP53 mutations; conversely, fluke-negative CCAs (clusters 3/4) exhibit high copy-number alterations and PD-1/PD-L2 expression, or epigenetic mutations (IDH1/2, BAP1) and FGFR/PRKA-related gene rearrangements. Whole-genome analysis highlighted FGFR2 3' untranslated region deletion as a mechanism of FGFR2 upregulation. Integration of noncoding promoter mutations with protein-DNA binding profiles demonstrates pervasive modulation of H3K27me3-associated sites in CCA. Clusters 1 and 4 exhibit distinct DNA hypermethylation patterns targeting either CpG islands or shores-mutation signature and subclonality analysis suggests that these reflect different mutational pathways. Our results exemplify how genetics, epigenetics, and environmental carcinogens can interplay across different geographies to generate distinct molecular subtypes of cancer.Significance: Integrated whole-genome and epigenomic analysis of CCA on an international scale identifies new CCA driver genes, noncoding promoter mutations, and structural variants. CCA molecular landscapes differ radically by etiology, underscoring how distinct cancer subtypes in the same organ may arise through different extrinsic and intrinsic carcinogenic processes. Cancer Discov; 7(10); 1116-35. ©2017 AACR.This article is highlighted in the In This Issue feature, p. 1047.

Genomic and proteomic characterization of ARID1A chromatin remodeller in ampullary tumors

AT rich interactive domain 1A (ARID1A) is one of the most commonly mutated genes in a broad variety of tumors. The mechanisms that involve ARID1A in ampullary cancer progression remains elusive. Here, we evaluated the frequency of ARID1A and KRAS mutations in ampullary adenomas and adenocarcinomas and in duodenal adenocarcinomas from two cohorts of patients from Singapore and Romania, correlated with clinical and pathological tumor features, and assessed the functional role of ARID1A. In the ampullary adenocarcinomas, the frequency of KRAS and ARID1A mutations was 34.7% and 8.2% respectively, with a loss or reduction of ARID1A protein in 17.2% of the cases. ARID1A mutational status was significantly correlated with ARID1A protein expression level (P=0.023). There was a significant difference in frequency of ARID1A mutation between Romania and Singapore (2.7% versus 25%, P=0.04), suggestive of different etiologies. One somatic mutation was detected in the ampullary adenoma group. In vitro studies indicated the tumor suppressive role of ARID1A. Our results warrant further investigation of this chromatin remodeller as a potential early biomarker of the disease, as well as identification of therapeutic targets in ARID1A mutated ampullary cancers.

JAK-STAT and G-protein-coupled receptor signaling pathways are frequently altered in epitheliotropic intestinal T-cell lymphoma

Epitheliotropic intestinal T-cell lymphoma (EITL, also known as type II enteropathy-associated T-cell lymphoma) is an aggressive intestinal disease with poor prognosis and its molecular alterations have not been comprehensively characterized. We aimed to identify actionable easy-to-screen alterations that would allow better diagnostics and/or treatment of this deadly disease. By performing whole-exome sequencing of four EITL tumor-normal pairs, followed by amplicon deep sequencing of 42 tumor samples, frequent alterations of the JAK-STAT and G-protein-coupled receptor (GPCR) signaling pathways were discovered in a large portion of samples. Specifically, STAT5B was mutated in a remarkable 63% of cases, JAK3 in 35% and GNAI2 in 24%, with the majority occurring at known activating hotspots in key functional domains. Moreover, STAT5B locus carried copy-neutral loss of heterozygosity resulting in the duplication of the mutant copy, suggesting the importance of mutant STAT5B dosage for the development of EITL. Dysregulation of the JAK-STAT and GPCR pathways was also supported by gene expression profiling and further verified in patient tumor samples. In vitro overexpression of GNAI2 mutants led to the upregulation of pERK1/2, a member of MEK-ERK pathway. Notably, inhibitors of both JAK-STAT and MEK-ERK pathways effectively reduced viability of patient-derived primary EITL cells, indicating potential therapeutic strategies for this neoplasm with no effective treatment currently available.

Abundant copy-number loss of CYCLOPS and STOP genes in gastric adenocarcinoma

BACKGROUND

Gastric cancer, a leading cause of cancer death worldwide, has been little studied compared with other cancers that impose similar health burdens. Our goal is to assess genomic copy-number loss and the possible functional consequences and therapeutic implications thereof across a large series of gastric adenocarcinomas.

METHODS

We used high-density single-nucleotide polymorphism microarrays to determine patterns of copy-number loss and allelic imbalance in 74 gastric adenocarcinomas. We investigated whether suppressor of tumorigenesis and/or proliferation (STOP) genes are associated with genomic copy-number loss. We also analyzed the extent to which copy-number loss affects Copy-number alterations Yielding Cancer Liabilities Owing to Partial losS (CYCLOPS) genes-genes that may be attractive targets for therapeutic inhibition when partially deleted.

RESULTS

The proportion of the genome subject to copy-number loss varies considerably from tumor to tumor, with a median of 5.5 %, and a mean of 12 % (range 0-58.5 %). On average, 91 STOP genes were subject to copy-number loss per tumor (median 35, range 0-452), and STOP genes tended to have lower copy-number compared with the rest of the genes. Furthermore, on average, 1.6 CYCLOPS genes per tumor were both subject to copy-number loss and downregulated, and 51.4 % of the tumors had at least one such gene.

CONCLUSIONS

The enrichment of STOP genes in regions of copy-number loss indicates that their deletion may contribute to gastric carcinogenesis. Furthermore, the presence of several deleted and downregulated CYCLOPS genes in some tumors suggests potential therapeutic targets in these tumors.

Mutational landscapes of tongue carcinoma reveal recurrent mutations in genes of therapeutic and prognostic relevance

Background

Carcinoma of the oral tongue (OTSCC) is the most common malignancy of the oral cavity, characterized by frequent recurrence and poor survival. The last three decades has witnessed a change in the OTSCC epidemiological profile, with increasing incidence in younger patients, females and never-smokers. Here, we sought to characterize the OTSCC genomic landscape and to determine factors that may delineate the genetic basis of this disease, inform prognosis and identify targets for therapeutic intervention.

Methods

Seventy-eight cases were subjected to whole-exome (n = 18) and targeted deep sequencing (n = 60).

Results

While the most common mutation was in TP53, the OTSCC genetic landscape differed from previously described cohorts of patients with head and neck tumors: OTSCCs demonstrated frequent mutations in DST and RNF213, while alterations in CDKN2A and NOTCH1 were significantly less frequent. Despite a lack of previously reported NOTCH1 mutations, integrated analysis showed enrichments of alterations affecting Notch signaling in OTSCC. Importantly, these Notch pathway alterations were prognostic on multivariate analyses. A high proportion of OTSCCs also presented with alterations in drug targetable and chromatin remodeling genes. Patients harboring mutations in actionable pathways were more likely to succumb from recurrent disease compared with those who did not, suggesting that the former should be considered for treatment with targeted compounds in future trials.

Conclusions

Our study defines the Asian OTSCC mutational landscape, highlighting the key role of Notch signaling in oral tongue tumorigenesis. We also observed somatic mutations in multiple therapeutically relevant genes, which may represent candidate drug targets in this highly lethal tumor type.

Electronic supplementary material

The online version of this article (doi:10.1186/s13073-015-0219-2) contains supplementary material, which is available to authorized users.

Exome-wide Sequencing Shows Low Mutation Rates and Identifies Novel Mutated Genes in Seminomas

BACKGROUND

Testicular germ cell tumors are the most common cancer diagnosed in young men, and seminomas are the most common type of these cancers. There have been no exome-wide examinations of genes mutated in seminomas or of overall rates of nonsilent somatic mutations in these tumors.

OBJECTIVE

The objective was to analyze somatic mutations in seminomas to determine which genes are affected and to determine rates of nonsilent mutations.

DESIGN, SETTING, AND PARTICIPANTS

Eight seminomas and matched normal samples were surgically obtained from eight patients.

INTERVENTION

DNA was extracted from tissue samples and exome sequenced on massively parallel Illumina DNA sequencers. Single-nucleotide polymorphism chip-based copy number analysis was also performed to assess copy number alterations.

OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS

The DNA sequencing read data were analyzed to detect somatic mutations including single-nucleotide substitutions and short insertions and deletions. The detected mutations were validated by independent sequencing and further checked for subclonality.

RESULTS AND LIMITATIONS

The rate of nonsynonymous somatic mutations averaged 0.31 mutations/Mb. We detected nonsilent somatic mutations in 96 genes that were not previously known to be mutated in seminomas, of which some may be driver mutations. Many of the mutations appear to have been present in subclonal populations. In addition, two genes, KIT and KRAS, were affected in two tumors each with mutations that were previously observed in other cancers and are presumably oncogenic.

CONCLUSIONS

Our study, the first report on exome sequencing of seminomas, detected somatic mutations in 96 new genes, several of which may be targetable drivers. Furthermore, our results show that seminoma mutation rates are five times higher than previously thought, but are nevertheless low compared to other common cancers. Similar low rates are seen in other cancers that also have excellent rates of remission achieved with chemotherapy.

PATIENT SUMMARY

We examined the DNA sequences of seminomas, the most common type of testicular germ cell cancer. Our study identified 96 new genes in which mutations occurred during seminoma development, some of which might contribute to cancer development or progression. The study also showed that the rates of DNA mutations during seminoma development are higher than previously thought, but still lower than for other common solid-organ cancers. Such low rates are also observed among other cancers that, like seminomas, show excellent rates of disease remission after chemotherapy.

Abstract 5168: Gastric adenocarcinomas show pervasive loss of heterozygosity and enrichment for “STOP” genes in regions of hemizygous deletion

Purpose: Gastric cancer, the second leading cause of cancer death worldwide, has been little studied compared with other cancers that impose similar burdens on public health. Our goal is to assess loss of heterozygosity (LOH) and hemizygous deletion across a large series of gastric adenocarcinomas using the most sensitive and comprehensive method currently available.

Experimental design: We used high-density single-nucleotide-polymorphism microarrays that assay genotype and copy number at 906,600 sites and analyzed the results to determine patterns of LOH and hemizygous deletion in 77 gastric adenocarcinomas. We investigated whether “STOP” genes - genes that tend to impede proliferation - are associated with hemizygous deletion. We also analyzed the extent to which hemizygous deletion affected CYCLOPS (Copy number alterations Yielding Cancer Liabilities Owing to Partial losS) genes - genes that may be attractive targets for therapeutic inhibition when hemizygously deleted.

Results: We found that LOH is pervasive: on average 27% of each tumor genome is subject to LOH, and > 98% of single nucleotide polymorphisms assayed were subject to LOH in at least 10% of tumors. Furthermore, 26% of LOH was due to hemizygous deletion, and STOP genes were significantly associated with regions of frequent hemizygous deletion; on average, 61 STOP genes were hemizygously deleted per tumor. Furthermore, on average, 4.65 CYCLOPS genes were hemizygously deleted per tumor, and 54.5% of the tumors had at least one CYCLOPS gene deleted.

Conclusions: These findings suggest that hemizygous deletion of anti-proliferative STOP genes contributes substantially to gastric carcinogenesis. Furthermore, the presence of several hemizygously deleted CYLOPS genes in some tumors suggests potential therapeutic targets in these tumors.

Citation Format: Ioana Cutcutache, Alice Yingting Wu, John R. McPherson, Zhengdeng Lei, Niantao Deng, Wai Keong Wong, Khee Chee Soo, Weng Hoong Chan, London Lucien Ooi, Roy E. Welsch, Patrick Tan, Steven G. Rozen. Gastric adenocarcinomas show pervasive loss of heterozygosity and enrichment for “STOP” genes in regions of hemizygous deletion. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 5168. doi:10.1158/1538-7445.AM2014-5168

Whole-genome reconstruction and mutational signatures in gastric cancer

Background

Gastric cancer is the second highest cause of global cancer mortality. To explore the complete repertoire of somatic alterations in gastric cancer, we combined massively parallel short read and DNA paired-end tag sequencing to present the first whole-genome analysis of two gastric adenocarcinomas, one with chromosomal instability and the other with microsatellite instability.

Results

Integrative analysis and de novo assemblies revealed the architecture of a wild-type KRAS amplification, a common driver event in gastric cancer. We discovered three distinct mutational signatures in gastric cancer - against a genome-wide backdrop of oxidative and microsatellite instability-related mutational signatures, we identified the first exome-specific mutational signature. Further characterization of the impact of these signatures by combining sequencing data from 40 complete gastric cancer exomes and targeted screening of an additional 94 independent gastric tumors uncovered ACVR2A, RPL22 and LMAN1 as recurrently mutated genes in microsatellite instability-positive gastric cancer and PAPPA as a recurrently mutated gene in TP53 wild-type gastric cancer.

Conclusions

These results highlight how whole-genome cancer sequencing can uncover information relevant to tissue-specific carcinogenesis that would otherwise be missed from exome-sequencing data.

Primer3--new capabilities and interfaces

Polymerase chain reaction (PCR) is a basic molecular biology technique with a multiplicity of uses, including deoxyribonucleic acid cloning and sequencing, functional analysis of genes, diagnosis of diseases, genotyping and discovery of genetic variants. Reliable primer design is crucial for successful PCR, and for over a decade, the open-source Primer3 software has been widely used for primer design, often in high-throughput genomics applications. It has also been incorporated into numerous publicly available software packages and web services. During this period, we have greatly expanded Primer3’s functionality. In this article, we describe Primer3’s current capabilities, emphasizing recent improvements. The most notable enhancements incorporate more accurate thermodynamic models in the primer design process, both to improve melting temperature prediction and to reduce the likelihood that primers will form hairpins or dimers. Additional enhancements include more precise control of primer placement—a change motivated partly by opportunities to use whole-genome sequences to improve primer specificity. We also added features to increase ease of use, including the ability to save and re-use parameter settings and the ability to require that individual primers not be used in more than one primer pair. We have made the core code more modular and provided cleaner programming interfaces to further ease integration with other software. These improvements position Primer3 for continued use with genome-scale data in the decade ahead.

Primer-BLAST: a tool to design target-specific primers for polymerase chain reaction

BACKGROUND

Choosing appropriate primers is probably the single most important factor affecting the polymerase chain reaction (PCR). Specific amplification of the intended target requires that primers do not have matches to other targets in certain orientations and within certain distances that allow undesired amplification. The process of designing specific primers typically involves two stages. First, the primers flanking regions of interest are generated either manually or using software tools; then they are searched against an appropriate nucleotide sequence database using tools such as BLAST to examine the potential targets. However, the latter is not an easy process as one needs to examine many details between primers and targets, such as the number and the positions of matched bases, the primer orientations and distance between forward and reverse primers. The complexity of such analysis usually makes this a time-consuming and very difficult task for users, especially when the primers have a large number of hits. Furthermore, although the BLAST program has been widely used for primer target detection, it is in fact not an ideal tool for this purpose as BLAST is a local alignment algorithm and does not necessarily return complete match information over the entire primer range.

RESULTS

We present a new software tool called Primer-BLAST to alleviate the difficulty in designing target-specific primers. This tool combines BLAST with a global alignment algorithm to ensure a full primer-target alignment and is sensitive enough to detect targets that have a significant number of mismatches to primers. Primer-BLAST allows users to design new target-specific primers in one step as well as to check the specificity of pre-existing primers. Primer-BLAST also supports placing primers based on exon/intron locations and excluding single nucleotide polymorphism (SNP) sites in primers.

CONCLUSIONS

We describe a robust and fully implemented general purpose primer design tool that designs target-specific PCR primers. Primer-BLAST offers flexible options to adjust the specificity threshold and other primer properties. This tool is publicly available at http://www.ncbi.nlm.nih.gov/tools/primer-blast.

Janus kinase 3-activating mutations identified in natural killer/T-cell lymphoma

The molecular pathogenesis of natural killer/T-cell lymphoma (NKTCL) is not well understood. We conducted whole-exome sequencing and identified Janus kinase 3 (JAK3) somatic-activating mutations (A572V and A573V) in 2 of 4 patients with NKTCLs. Further validation of the prevalence of JAK3 mutations was determined by Sanger sequencing and high-resolution melt (HRM) analysis in an additional 61 cases. In total, 23 of 65 (35.4%) cases harbored JAK3 mutations. Functional characterization of the JAK3 mutations support its involvement in cytokine-independent JAK/STAT constitutive activation leading to increased cell growth. Moreover, treatment of both JAK3-mutant and wild-type NKTCL cell lines with a novel pan-JAK inhibitor, CP-690550, resulted in dose-dependent reduction of phosphorylated STAT5, reduced cell viability, and increased apoptosis. Hence, targeting the deregulated JAK/STAT pathway could be a promising therapy for patients with NKTCLs.

SIGNIFICANCE

Gene mutations causing NKTCL have not been fully identified. Through exome sequencing, we identified activating mutations of JAK3 that may play a significant role in the pathogenesis of NKTCLs. Our findings have important implications for the management of patients with NKTCLs.

Exome sequencing of liver fluke-associated cholangiocarcinoma

Opisthorchis viverrini-related cholangiocarcinoma (CCA), a fatal bile duct cancer, is a major public health concern in areas endemic for this parasite. We report here whole-exome sequencing of eight O. viverrini-related tumors and matched normal tissue. We identified and validated 206 somatic mutations in 187 genes using Sanger sequencing and selected 15 genes for mutation prevalence screening in an additional 46 individuals with CCA (cases). In addition to the known cancer-related genes TP53 (mutated in 44.4% of cases), KRAS (16.7%) and SMAD4 (16.7%), we identified somatic mutations in 10 newly implicated genes in 14.8-3.7% of cases. These included inactivating mutations in MLL3 (in 14.8% of cases), ROBO2 (9.3%), RNF43 (9.3%) and PEG3 (5.6%), and activating mutations in the GNAS oncogene (9.3%). These genes have functions that can be broadly grouped into three biological classes: (i) deactivation of histone modifiers, (ii) activation of G protein signaling and (iii) loss of genome stability. This study provides insight into the mutational landscape contributing to O. viverrini-related CCA.

First somatic mutation of E2F1 in a critical DNA binding residue discovered in well-differentiated papillary mesothelioma of the peritoneum

Background

Well differentiated papillary mesothelioma of the peritoneum (WDPMP) is a rare variant of epithelial mesothelioma of low malignancy potential, usually found in women with no history of asbestos exposure. In this study, we perform the first exome sequencing of WDPMP.

Results

WDPMP exome sequencing reveals the first somatic mutation of E2F1, R166H, to be identified in human cancer. The location is in the evolutionarily conserved DNA binding domain and computationally predicted to be mutated in the critical contact point between E2F1 and its DNA target. We show that the R166H mutation abrogates E2F1's DNA binding ability and is associated with reduced activation of E2F1 downstream target genes. Mutant E2F1 proteins are also observed in higher quantities when compared with wild-type E2F1 protein levels and the mutant protein's resistance to degradation was found to be the cause of its accumulation within mutant over-expressing cells. Cells over-expressing wild-type E2F1 show decreased proliferation compared to mutant over-expressing cells, but cell proliferation rates of mutant over-expressing cells were comparable to cells over-expressing the empty vector.

Conclusions

The R166H mutation in E2F1 is shown to have a deleterious effect on its DNA binding ability as well as increasing its stability and subsequent accumulation in R166H mutant cells. Based on the results, two compatible theories can be formed: R166H mutation appears to allow for protein over-expression while minimizing the apoptotic consequence and the R166H mutation may behave similarly to SV40 large T antigen, inhibiting tumor suppressive functions of retinoblastoma protein 1.

Assessing matched normal and tumor pairs in next-generation sequencing studies

Next generation sequencing technology has revolutionized the study of cancers. Through matched normal-tumor pairs, it is now possible to identify genome-wide germline and somatic mutations. The generation and analysis of the data requires rigorous quality checks and filtering, and the current analytical pipeline is constantly undergoing improvements. We noted however that in analyzing matched pairs, there is an implicit assumption that the sequenced data are matched, without any quality check such as those implemented in association studies. There are serious implications in this assumption as identification of germline and rare somatic variants depend on the normal sample being the matched pair. Using a genetics concept on measuring relatedness between individuals, we demonstrate that the matchedness of tumor pairs can be quantified and should be included as part of a quality protocol in analysis of sequenced data. Despite the mutation changes in cancer samples, matched tumor-normal pairs are still relatively similar in sequence compared to non-matched pairs. We demonstrate that the approach can be used to assess the mutation landscape between individuals.

Genetic and structural variation in the gastric cancer kinome revealed through targeted deep sequencing

Genetic alterations in kinases have been linked to multiple human pathologies. To explore the landscape of kinase genetic variation in gastric cancer (GC), we used targeted, paired-end deep sequencing to analyze 532 protein and phosphoinositide kinases in 14 GC cell lines. We identified 10,604 single-nucleotide variants (SNV) in kinase exons including greater than 300 novel nonsynonymous SNVs. Family-wise analysis of the nonsynonymous SNVs revealed a significant enrichment in mitogen-activated protein kinase (MAPK)-related genes (P < 0.01), suggesting a preferential involvement of this kinase family in GC. A potential antioncogenic role for MAP2K4, a gene exhibiting recurrent alterations in 2 lines, was functionally supported by siRNA knockdown and overexpression studies in wild-type and MAP2K4 variant lines. The deep sequencing data also revealed novel, large-scale structural rearrangement events involving kinases including gene fusions involving CDK12 and the ERBB2 receptor tyrosine kinase in MKN7 cells. Integrating SNVs and copy number alterations, we identified Hs746T as a cell line exhibiting both splice-site mutations and genomic amplification of MET, resulting in MET protein overexpression. When applied to primary GCs, we identified somatic mutations in 8 kinases, 4 of which were recurrently altered in both primary tumors and cell lines (MAP3K6, STK31, FER, and CDKL5). These results demonstrate that how targeted deep sequencing approaches can deliver unprecedented multilevel characterization of a medically and pharmacologically relevant gene family. The catalog of kinome genetic variants assembled here may broaden our knowledge on kinases and provide useful information on genetic alterations in GC.