Early-Stage Induction of SWI/SNF Mutations during Esophageal Squamous Cell Carcinogenesis

The SWI/SNF chromatin remodeling complex is frequently inactivated by somatic mutations of its various components in various types of cancers, and also by aberrant DNA methylation. However, its somatic mutations and aberrant methylation in esophageal squamous cell carcinomas (ESCCs) have not been fully analyzed. In this study, we aimed to clarify in ESCC, what components of the SWI/SNF complex have somatic mutations and aberrant methylation, and when somatic mutations of the SWI/SNF complex occur. Deep sequencing of components of the SWI/SNF complex using a bench-top next generation sequencer revealed that eight of 92 ESCCs (8.7%) had 11 somatic mutations of 7 genes, ARID1A, ARID2, ATRX, PBRM1, SMARCA4, SMARCAL1, and SMARCC1. The SMARCA4 mutations were located in the Forkhead (85Ser>Leu) and SNF2 family N-terminal (882Glu>Lys) domains. The PBRM1 mutations were located in a bromodomain (80Asn>Ser) and an HMG-box domain (1,377Glu>Lys). For most mutations, their mutant allele frequency was 31–77% (mean 61%) of the fraction of cancer cells in the same samples, indicating that most of the cancer cells in individual ESCC samples had the SWI/SNF mutations on one allele, when present. In addition, a BeadChip array analysis revealed that a component of the SWI/SNF complex, ACTL6B, had aberrant methylation at its promoter CpG island in 18 of 52 ESCCs (34.6%). These results showed that genetic and epigenetic alterations of the SWI/SNF complex are present in ESCCs, and suggested that genetic alterations are induced at an early stage of esophageal squamous cell carcinogenesis.     

Local sequence assembly reveals a high-resolution profile of somatic structural variations in 97 cancer genomes

Genomic structural variations (SVs) are pervasive in many types of cancers. Characterizing their underlying mechanisms and potential molecular consequences is crucial for understanding the basic biology of tumorigenesis. Here, we engineered a local assembly-based algorithm (laSV) that detects SVs with high accuracy from paired-end high-throughput genomic sequencing data and pinpoints their breakpoints at single base-pair resolution. By applying laSV to 97 tumor-normal paired genomic sequencing datasets across six cancer types produced by The Cancer Genome Atlas Research Network, we discovered that non-allelic homologous recombination is the primary mechanism for generating somatic SVs in acute myeloid leukemia. This finding contrasts with results for the other five types of solid tumors, in which non-homologous end joining and microhomology end joining are the predominant mechanisms. We also found that the genes recursively mutated by single nucleotide alterations differed from the genes recursively mutated by SVs, suggesting that these two types of genetic alterations play different roles during cancer progression. We further characterized how the gene structures of the oncogene JAK1 and the tumor suppressors KDM6A and RB1 are affected by somatic SVs and discussed the potential functional implications of intergenic SVs.     

Downregulation of miR-503 Promotes ESCC Cell Proliferation,Migration, and Invasion by Targeting Cyclin D1

Esophageal squamous cell carcinoma (ESCC) is one of the most aggressive cancers in China, but the underlying molecular mechanism of ESCC is still unclear. Involvement of microRNAs has been demonstrated in cancer initiation and progression. Despite the reported function of miR-503 in several human cancers, its detailed anti-oncogenic role and clinical significance in ESCC remain undefined. In this study, we examined miR-503 expression by qPCR and found the downregulation of miR-503 expression in ESCC tissue relative to adjacent normal tissues. Further investigation in the effect of miR-503 on ESCC cell proliferation, migration, and invasion showed that enhanced expression of miR-503 inhibited ESCC aggressive phenotype and overexpression of CCND1 reversed the effect of miR-503-mediated ESCC cell aggressive phenotype. Our study further identified CCND1 as the target gene of miR-503. Thus, miR-503 functions as a tumor suppressor and has an important role in ESCC by targeting CCND1.     

Genomic Analyses Reveal Mutational Signatures and Frequently Altered Genes in Esophageal Squamous Cell Carcinoma

Esophageal squamous cell carcinoma (ESCC) is one of the most common cancers worldwide and the fourth most lethal cancer in China. However, although genomic studies have identified some mutations associated with ESCC, we know little of the mutational processes responsible. To identify genome-wide mutational signatures, we performed either whole-genome sequencing (WGS) or whole-exome sequencing (WES) on 104 ESCC individuals and combined our data with those of 88 previously reported samples. An APOBEC-mediated mutational signature in 47% of 192 tumors suggests that APOBEC-catalyzed deamination provides a source of DNA damage in ESCC. Moreover, PIK3CA hotspot mutations (c.1624G>A [p.Glu542Lys] and c.1633G>A [p.Glu545Lys]) were enriched in APOBEC-signature tumors, and no smoking-associated signature was observed in ESCC. In the samples analyzed by WGS, we identified focal (<100 kb) amplifications of CBX4 and CBX8. In our combined cohort, we identified frequent inactivating mutations in AJUBA, ZNF750, and PTCH1 and the chromatin-remodeling genes CREBBP and BAP1, in addition to known mutations. Functional analyses suggest roles for several genes (CBX4, CBX8, AJUBA, and ZNF750) in ESCC. Notably, high activity of hedgehog signaling and the PI3K pathway in approximately 60% of 104 ESCC tumors indicates that therapies targeting these pathways might be particularly promising strategies for ESCC. Collectively, our data provide comprehensive insights into the mutational signatures of ESCC and identify markers for early diagnosis and potential therapeutic targets.     

Evolutionary Genomics of Structural Variation in Asian Rice (Oryza sativa) Domestication

Structural variants (SVs) are a largely unstudied feature of plant genome evolution, despite the fact that SVs contribute substantially to phenotypes. In this study, we discovered SVs across a population sample of 347 high-coverage, resequenced genomes of Asian rice (Oryza sativa) and its wild ancestor (O. rufipogon). In addition to this short-read data set, we also inferred SVs from whole-genome assemblies and long-read data. Comparisons among data sets revealed different features of genome variability. For example, genome alignment identified a large (∼4.3 Mb) inversion in indica rice varieties relative to japonica varieties, and long-read analyses suggest that ∼9% of genes from the outgroup (O. longistaminata) are hemizygous. We focused, however, on the resequencing sample to investigate the population genomics of SVs. Clustering analyses with SVs recapitulated the rice cultivar groups that were also inferred from SNPs. However, the site-frequency spectrum of each SV type—which included inversions, duplications, deletions, translocations, and mobile element insertions—was skewed toward lower frequency variants than synonymous SNPs, suggesting that SVs may be predominantly deleterious. Among transposable elements, SINE and mariner insertions were found at especially low frequency. We also used SVs to study domestication by contrasting between rice and O. rufipogon. Cultivated genomes contained ∼25% more derived SVs and mobile element insertions than O. rufipogon, indicating that SVs contribute to the cost of domestication in rice. Peaks of SV divergence were enriched for known domestication genes, but we also detected hundreds of genes gained and lost during domestication, some of which were enriched for traits of agronomic interest.     

Genetic Aberrations in Imatinib-Resistant Dermatofibrosarcoma Protuberans Revealed by Whole Genome Sequencing

Dermatofibrosarcoma protuberans (DFSP) is a very rare soft tissue sarcoma. DFSP often reveals a specific chromosome translocation, t(17;22)(q22;q13), which results in the fusion of collagen 1 alpha 1 (COL1A1) gene and platelet-derived growth factor-B (PDGFB) gene. The COL1A1-PDGFB fusion protein activates the PDGFB receptor and resultant constitutive activation of PDGFR receptor is essential in the pathogenesis of DFSP. Thus, blocking PDGFR receptor activation with imatinib has shown promising activity in the treatment of advanced and metastatic DFSP. Despite the success with targeted agents in cancers, acquired drug resistance eventually occurs. Here, we tried to identify potential drug resistance mechanisms against imatinib in a 46-year old female with DFSP who initially responded well to imatinib but suffered rapid disease progression. We performed whole-genome sequencing of both pre-treatment and post-treatment tumor tissue to identify the mutational events associated with imatinib resistance. No significant copy number alterations, insertion, and deletions were identified during imatinib treatment. Of note, we identified newly emerged 8 non-synonymous somatic mutations of the genes (ACAP2, CARD10, KIAA0556, PAAQR7, PPP1R39, SAFB2, STARD9, and ZFYVE9) in the imatinib-resistant tumor tissue. This study revealed diverse possible candidate mechanisms by which imatinib resistance to PDGFRB inhibition may arise in DFSP, and highlights the usefulness of whole-genome sequencing in identifying drug resistance mechanisms and in pursuing genome-directed, personalized anti-cancer therapy.     

Chromothripsis is a common mechanism driving genomic rearrangements in primary and metastatic colorectal cancer

Background

Structural rearrangements form a major class of somatic variation in cancer genomes. Local chromosome shattering, termed chromothripsis, is a mechanism proposed to be the cause of clustered chromosomal rearrangements and was recently described to occur in a small percentage of tumors. The significance of these clusters for tumor development or metastatic spread is largely unclear.

Results

We used genome-wide long mate-pair sequencing and SNP array profiling to reveal that chromothripsis is a widespread phenomenon in primary colorectal cancer and metastases. We find large and small chromothripsis events in nearly every colorectal tumor sample and show that several breakpoints of chromothripsis clusters and isolated rearrangements affect cancer genes, including NOTCH2, EXO1 and MLL3. We complemented the structural variation studies by sequencing the coding regions of a cancer exome in all colorectal tumor samples and found somatic mutations in 24 genes, including APC, KRAS, SMAD4 and PIK3CA. A pairwise comparison of somatic variations in primary and metastatic samples indicated that many chromothripsis clusters, isolated rearrangements and point mutations are exclusively present in either the primary tumor or the metastasis and may affect cancer genes in a lesion-specific manner.

Conclusions

We conclude that chromothripsis is a prevalent mechanism driving structural rearrangements in colorectal cancer and show that a complex interplay between point mutations, simple copy number changes and chromothripsis events drive colorectal tumor development and metastasis.     

CCGD-ESCC: A Comprehensive Database for Genetic Variants Associated with Esophageal Squamous Cell Carcinoma in Chinese Population

Esophageal squamous-cell carcinoma (ESCC) is one of the most lethal malignancies in the world and occurs at particularly higher frequency in China. While several genome-wide association studies (GWAS) of germline variants and whole-genome or whole-exome sequencing studies of somatic mutations in ESCC have been published, there is no comprehensive database publically available for this cancer. Here, we developed the Chinese Cancer Genomic Database-Esophageal Squamous Cell Carcinoma (CCGD-ESCC) database, which contains the associations of 69,593 single nucleotide polymorphisms (SNPs) with ESCC risk in 2022 cases and 2039 controls, survival time of 1006 ESCC patients (survival GWAS) and gene expression (expression quantitative trait loci, eQTL) in 94 ESCC patients. Moreover, this database also provides the associations between 8833 somatic mutations and survival time in 675 ESCC patients. Our user-friendly database is a resource useful for biologists and oncologists not only in identifying the associations of genetic variants or somatic mutations with the development and progression of ESCC but also in studying the underlying mechanisms for tumorigenesis of the cancer. CCGD-ESCC is freely accessible at http://db.cbi.pku.edu.cn/ccgd/ESCCdb.     

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.     

An Integrative Analysis to Identify Driver Genes in Esophageal Squamous Cell Carcinoma

BackgroundFew driver genes have been well established in esophageal squamous cell carcinoma (ESCC). Identification of the genomic aberrations that contribute to changes in gene expression profiles can be used to predict driver genes.MethodsWe searched for driver genes in ESCC by integrative analysis of gene expression microarray profiles and copy number data. To narrow down candidate genes, we performed survival analysis on expression data and tested the genetic vulnerability of each genes using public RNAi screening data. We confirmed the results by performing RNAi experiments and evaluating the clinical relevance of candidate genes in an independent ESCC cohort.ResultsWe found 10 significantly recurrent copy number alterations accompanying gene expression changes, including loci 11q13.2, 7p11.2, 3q26.33, and 17q12, which harbored CCND1, EGFR, SOX2, and ERBB2, respectively. Analysis of survival data and RNAi screening data suggested that GRB7, located on 17q12, was a driver gene in ESCC. In ESCC cell lines harboring 17q12 amplification, knockdown of GRB7 reduced the proliferation, migration, and invasion capacities of cells. Moreover, siRNA targeting GRB7 had a synergistic inhibitory effect when combined with trastuzumab, an anti-ERBB2 antibody. Survival analysis of the independent cohort also showed that high GRB7 expression was associated with poor prognosis in ESCC.ConclusionOur integrative analysis provided important insights into ESCC pathogenesis. We identified GRB7 as a novel ESCC driver gene and potential new therapeutic target.     

Genetic Variants in Epidermal Growth Factor Receptor Pathway Genes and Risk of Esophageal Squamous Cell Carcinoma and Gastric Cancer in a Chinese Population

The epidermal growth factor receptor (EGFR) signaling pathway regulates cell proliferation, differentiation, and survival, and is frequently dysregulated in esophageal and gastric cancers. Few studies have comprehensively examined the association between germline genetic variants in the EGFR pathway and risk of esophageal and gastric cancers. Based on a genome-wide association study in a Han Chinese population, we examined 3443 SNPs in 127 genes in the EGFR pathway for 1942 esophageal squamous cell carcinomas (ESCCs), 1758 gastric cancers (GCs), and 2111 controls. SNP-level analyses were conducted using logistic regression models. We applied the resampling-based adaptive rank truncated product approach to determine the gene- and pathway-level associations. The EGFR pathway was significantly associated with GC risk (P = 2.16×10⁻³). Gene-level analyses found 10 genes to be associated with GC, including FYN, MAPK8, MAP2K4, GNAI3, MAP2K1, TLN1, PRLR, PLCG2, RPS6KB2, and PIK3R3 (P<0.05). For ESCC, we did not observe a significant pathway-level association (P = 0.72), but gene-level analyses suggested associations between GNAI3, CHRNE, PAK4, WASL, and ITCH, and ESCC (P<0.05). Our data suggest an association between specific genes in the EGFR signaling pathway and risk of GC and ESCC. Further studies are warranted to validate these associations and to investigate underlying mechanisms.     

Genetic and molecular analysis of synaptic vesicle recycling in Drosophila

Following exocytosis, one of the major presynaptic events is replenishing synaptic vesicles (SVs) to ensure the possibility of continuous synaptic transmission. The nerve terminal is thought to recycle SVs through clathrin-mediated endocytosis and by a clathrin-independent pathway called ‘kiss and run’. This review highlights the use of the genetic model organism, the fruit fly (Drosophila melanogaster), in dissecting the molecular mechanisms of clathrin-mediated endocytosis in recycling SVs at neuromuscular junctions (NMJs). Analyses of endocytotic mutants in Drosophila indicate that clathrin-mediated endocytosis may be essential for SV recycling, including a putative fast recycling mechanism uncovered recently. Further, a rather complex picture begins to emerge suggesting that clathrin-mediated endocytosis involves several sequential steps mediated by a large number of proteins. Finally, these studies also reveal that SV proteins may be selectively retrieved into nascent SVs by clathrin accessory proteins and defects in protein retrieval have significant impacts on synaptic transmission. Following the completion of the Drosophila Genome Project and the development of gene targeting and RNAi approaches, genetic studies in Drosophila have become increasingly efficient. Hence, Drosophila is expected to continue to serve as an important model organism for studies of SV recycling.     

Werner syndrome protein suppresses the formation of large deletions during the replication of human telomeric sequences

Werner syndrome (WS) is a disorder characterized by features of premature aging and increased cancer that is caused by loss of the RecQ helicase WRN. Telomeres consisting of duplex TTAGGG repeats in humans protect chromosome ends and sustain cellular proliferation. WRN prevents the loss of telomeres replicated from the G-rich strand, which can form secondary G-quadruplex (G4) structures. Here, we dissected WRN roles in the replication of telomeric sequences by examining factors inherent to telomeric repeats, such as G4 DNA, independently from other factors at chromosome ends that can also impede replication. For this we used the supF shuttle vector (SV) mutagenesis assay. We demonstrate that SVs with [TTAGGG]₆ sequences are stably replicated in human cells, and that the repeats suppress the frequency of large deletions despite G4 folding potential. WRN depletion increased the supF mutant frequency for both the telomeric and non-telomeric SVs, compared with the control cells, but this increase was much greater (27-fold) for telomeric SVs. The higher SV mutant frequencies in WRN-deficient cells were primarily due to an increase in large sequence deletions and rearrangements. However, WRN depletion caused a more dramatic increase in deletions and rearrangements arising within the telomeric SV (70-fold), compared with non-telomeric SV (8-fold). Our results indicate that WRN prevents large deletions and rearrangements during replication, and that this role is particularly important in templates with telomeric sequence. This provides a possible explanation for increased telomere loss in WS cells.     

Transposon Domestication versus Mutualism in Ciliate Genome Rearrangements

Ciliated protists rearrange their genomes dramatically during nuclear development via chromosome fragmentation and DNA deletion to produce a trimmer and highly reorganized somatic genome. The deleted portion of the genome includes potentially active transposons or transposon-like sequences that reside in the germline. Three independent studies recently showed that transposase proteins of the DDE/DDD superfamily are indispensible for DNA processing in three distantly related ciliates. In the spirotrich Oxytricha trifallax, high copy-number germline-limited transposons mediate their own excision from the somatic genome but also contribute to programmed genome rearrangement through a remarkable transposon mutualism with the host. By contrast, the genomes of two oligohymenophorean ciliates, Tetrahymena thermophila and Paramecium tetraurelia, encode homologous PiggyBac-like transposases as single-copy genes in both their germline and somatic genomes. These domesticated transposases are essential for deletion of thousands of different internal sequences in these species. This review contrasts the events underlying somatic genome reduction in three different ciliates and considers their evolutionary origins and the relationships among their distinct mechanisms for genome remodeling.     

Homology-mediated end-capping as a primary step of sister chromatid fusion in the breakage-fusion-bridge cycles

Breakage-fusion-bridge (BFB) cycle is a series of chromosome breaks and duplications that could lead to the increased copy number of a genomic segment (gene amplification). A critical step of BFB cycles leading to gene amplification is a palindromic fusion of sister chromatids following the rupture of a dicentric chromosome during mitosis. It is currently unknown how sister chromatid fusion is produced from a mitotic break. To delineate the process, we took an integrated genomic, cytogenetic and molecular approach for the recurrent MCL1 amplicon at chromosome 1 in human tumor cells. A newly developed next-generation sequencing-based approach identified a cluster of palindromic fusions within the amplicon at ∼50-kb intervals, indicating a series of breaks and fusions by BFB cycles. The physical location of the amplicon (at the end of a broken chromosome) further indicated BFB cycles as underlying processes. Three palindromic fusions were mediated by the homologies between two nearby inverted Alu repeats, whereas the other two fusions exhibited microhomology-mediated events. Such breakpoint sequences indicate that homology-mediated fold-back capping of broken ends followed by DNA replication is an underlying mechanism of sister chromatid fusion. Our results elucidate nucleotide-level events during BFB cycles and end processing for naturally occurring mitotic breaks.     

Cryptic and Complex Chromosomal Aberrations in Early-Onset Neuropsychiatric Disorders

Structural variation (SV) is a significant component of the genetic etiology of both neurodevelopmental and psychiatric disorders; however, routine guidelines for clinical genetic screening have been established only in the former category. Genome-wide chromosomal microarray (CMA) can detect genomic imbalances such as copy-number variants (CNVs), but balanced chromosomal abnormalities (BCAs) still require karyotyping for clinical detection. Moreover, submicroscopic BCAs and subarray threshold CNVs are intractable, or cryptic, to both CMA and karyotyping. Here, we performed whole-genome sequencing using large-insert jumping libraries to delineate both cytogenetically visible and cryptic SVs in a single test among 30 clinically referred youth representing a range of severe neuropsychiatric conditions. We detected 96 SVs per person on average that passed filtering criteria above our highest-confidence resolution (6,305 bp) and an additional 111 SVs per genome below this resolution. These SVs rearranged 3.8 Mb of genomic sequence and resulted in 42 putative loss-of-function (LoF) or gain-of-function mutations per person. We estimate that 80% of the LoF variants were cryptic to clinical CMA. We found myriad complex and cryptic rearrangements, including a “paired” duplication (360 kb, 169 kb) that flanks a 5.25 Mb inversion that appears in 7 additional cases from clinical CNV data among 47,562 individuals. Following convergent genomic profiling of these independent clinical CNV data, we interpreted three SVs to be of potential clinical significance. These data indicate that sequence-based delineation of the full SV mutational spectrum warrants exploration in youth referred for neuropsychiatric evaluation and clinical diagnostic SV screening more broadly.     

A recurrent deletion of DPY19L2 causes infertility in man by blocking sperm head elongation and acrosome formation

An increasing number of couples require medical assistance to achieve a pregnancy, and more than 2% of the births in Western countries now result from assisted reproductive technologies. To identify genetic variants responsible for male infertility, we performed a whole-genome SNP scan on patients presenting with total globozoospermia, a primary infertility phenotype characterized by the presence of 100% round acrosomeless spermatozoa in the ejaculate. This strategy allowed us to identify in most patients (15/20) a 200 kb homozygous deletion encompassing only DPY19L2, which is highly expressed in the testis. Although there was no known function for DPY19L2 in humans, previous work indicated that its ortholog in C. elegans is involved in cell polarity. In man, the DPY19L2 region has been described as a copy-number variant (CNV) found to be duplicated and heterozygously deleted in healthy individuals. We show here that the breakpoints of the deletions are located on a highly homologous 28 kb low copy repeat (LCR) sequence present on each side of DPY19L2, indicating that the identified deletions were probably produced by nonallelic homologous recombination (NAHR) between these two regions. We demonstrate that patients with globozoospermia have a homozygous deletion of DPY19L2, thus indicating that DPY19L2 is necessary in men for sperm head elongation and acrosome formation. A molecular diagnosis can now be proposed to affected men; the presence of the deletion confirms the diagnosis of globozoospermia and assigns a poor prognosis for the success of in vitro fertilization.     

Oral Microbiota and Risk for Esophageal Squamous Cell Carcinoma in a High-Risk Area of China

Poor oral health has been linked with an increased risk of esophageal squamous cell carcinoma (ESCC). We investigated whether alteration of oral microbiota is associated with ESCC risk. Fasting saliva samples were collected from 87 incident and histopathologicallly diagnosed ESCC cases, 63 subjects with dysplasia and 85 healthy controls. All subjects were also interviewed with a questionnaire. V3–V4 region of 16S rRNA was amplified and sequenced by 454-pyrosequencing platform. Carriage of each genus was compared by means of multivariate-adjusted odds ratios derived from logistic regression model. Relative abundance was compared using Metastats method. Beta diversity was estimated using Unifrac and weighted Unifrac distances. Principal coordinate analysis (PCoA) was applied to ordinate dissimilarity matrices. Multinomial logistic regression was used to compare the coordinates between different groups. ESCC subjects had an overall decreased microbial diversity compared to control and dysplasia subjects (P<0.001). Decreased carriage of genera Lautropia, Bulleidia, Catonella, Corynebacterium, Moryella, Peptococcus and Cardiobacterium were found in ESCC subjects compared to non-ESCC subjects. Multinomial logistic regression analyses on PCoA coordinates also revealed that ESCC subjects had significantly different levels for several coordinates compared to non-ESCC subjects. In conclusion, we observed a correlation between altered salivary bacterial microbiota and ESCC risk. The results of our study on the saliva microbiome are of particular interest as it reflects the shift in microbial communities. Further studies are warranted to verify this finding, and if being verified, to explore the underlying mechanisms.