Targeted single-cell proteomic analysis identifies new liquid biopsy biomarkers associated with multiple myeloma.

Multiple myeloma (MM) is accompanied by alterations to the normal plasma cell (PC) proteome, leading to changes to the tumor microenvironment and disease progression. There is a great need for understanding the consequences that lead to MM progression and for the discovery of new biomarkers that can aid clinical diagnostics and serve as targets for therapeutics. This study demonstrates the applicability of utilizing the single-cell high-definition liquid biopsy assay (HDSCA) and imaging mass cytometry to characterize the proteomic profile of myeloma. In our study, we analyzed ~87,000 cells from seven patient samples (bone marrow and peripheral blood) across the myeloma disease spectrum and utilized our multiplexed panel to characterize the expression of clinical markers for PC classification, additional potential therapeutic targets, and the tumor microenvironment cells. Our analysis showed BCMA, ICAM3 (CD50), CD221, and CS1 (SLAMF7) as the most abundantly expressed markers on PCs across all myeloma stages, with BCMA, ICAM3, and CD221 having significantly higher expression levels on disease versus precursor PCs. Additionally, we identify significantly elevated levels of expression for CD74, MUM1, CD229, CD44, IGLL5, Cyclin D1, UBA52, and CD317 on PCs from overt disease conditions compared to those from precursor states.

Single cell clonotypic and transcriptional evolution of multiple myeloma precursor disease.

Multiple myeloma remains an incurable disease, and the cellular and molecular evolution from precursor conditions, including monoclonal gammopathy of undetermined significance and smoldering multiple myeloma, is incompletely understood. Here, we combine single-cell RNA and B cell receptor sequencing from fifty-two patients with myeloma precursors in comparison with myeloma and normal donors. Our comprehensive analysis reveals early genomic drivers of malignant transformation, distinct transcriptional features, and divergent clonal expansion in hyperdiploid versus non-hyperdiploid samples. Additionally, we observe intra-patient heterogeneity with potential therapeutic implications and identify distinct patterns of evolution from myeloma precursor disease to myeloma. We also demonstrate distinctive characteristics of the microenvironment associated with specific genomic changes in myeloma cells. These findings add to our knowledge about myeloma precursor disease progression, providing valuable insights into patient risk stratification, biomarker discovery, and possible clinical applications.

Noninvasive genomic profiling of somatic mutations in oral cavity cancers.

OBJECTIVES: Somatic mutations may predict prognosis, therapeutic response, or cancer progression. We evaluated targeted sequencing of oral rinse samples (ORS) for non-invasive mutational profiling of oral squamous cell carcinomas (OSCC). MATERIALS AND METHODS: A custom hybrid capture panel targeting 42 frequently mutated genes in OSCC was used to identify DNA sequence variants in matched ORS and fresh-frozen tumors from 120 newly-diagnosed patients. Receiver operating characteristic (ROC) curves determined the optimal variant allele fraction (VAF) cutoff for variant discrimination in ORS. Behavioral, clinical, and analytical factors were evaluated for impacts on assay performance. RESULTS: Half of tumors involved oral tongue (50 %), and a majority were T1-T2 tumor stage (55 %). Median depth of sequencing coverage was 260X for OSCC and 1,563X for ORS. Frequencies of single nucleotide variants (SNVs) at highly mutated genes (including TP53, FAT1, HRAS, NOTCH1, CDKN2A, CASP8, NFE2L2, and PIK3CA) in OSCC were highly correlated with TCGA data (R = 0.96, p = 2.5E-22). An ROC curve with area-under-the-curve (AUC) of 0.80 showed that, at an optimal VAF cutoff of 0.10 %, ORS provided 76 % sensitivity, 96 % specificity, but precision of only 2.6E-4. At this VAF cutoff, 206 of 270 SNVs in OSCC were detected in matched ORS. Sensitivity varied by patient, T stage and target gene. Neither downsampled ORS as matched control nor a naïve Bayesian classifier adjusting for sequencing bias appreciably improved assay performance. CONCLUSION: Targeted sequencing of ORS provides moderate assay performance for noninvasive detection of SNVs in OSCC. Our findings strongly rationalize further clinical and laboratory optimization of this assay, including strategies to improve precision.

Intratumoral Heterogeneity and Clonal Evolution Induced by HPV Integration.

The human papillomavirus (HPV) genome is integrated into host DNA in most HPV-positive cancers, but the consequences for chromosomal integrity are unknown. Continuous long-read sequencing of oropharyngeal cancers and cancer cell lines identified a previously undescribed form of structural variation, "heterocateny," characterized by diverse, interrelated, and repetitive patterns of concatemerized virus and host DNA segments within a cancer. Unique breakpoints shared across structural variants facilitated stepwise reconstruction of their evolution from a common molecular ancestor. This analysis revealed that virus and virus-host concatemers are unstable and, upon insertion into and excision from chromosomes, facilitate capture, amplification, and recombination of host DNA and chromosomal rearrangements. Evidence of heterocateny was detected in extrachromosomal and intrachromosomal DNA. These findings indicate that heterocateny is driven by the dynamic, aberrant replication and recombination of an oncogenic DNA virus, thereby extending known consequences of HPV integration to include promotion of intratumoral heterogeneity and clonal evolution. SIGNIFICANCE: Long-read sequencing of HPV-positive cancers revealed "heterocateny," a previously unreported form of genomic structural variation characterized by heterogeneous, interrelated, and repetitive genomic rearrangements within a tumor. Heterocateny is driven by unstable concatemerized HPV genomes, which facilitate capture, rearrangement, and amplification of host DNA, and promotes intratumoral heterogeneity and clonal evolution. See related commentary by McBride and White, p. 814. This article is highlighted in the In This Issue feature, p. 799.

Enrichment-Free Single-Cell Detection and Morphogenomic Profiling of Myeloma Patient Samples to Delineate Circulating Rare Plasma Cell Clones.

Multiple myeloma is an incurable malignancy that initiates from a bone marrow resident clonal plasma cell and acquires successive mutational changes and genomic alterations, eventually resulting in tumor burden accumulation and end-organ damage. It has been recently recognized that myeloma secondary genomic events result in extensive sub-clonal heterogeneity both in localized bone marrow areas and circulating peripheral blood plasma cells. Rare genomic subclones, including myeloma initiating cells, could be the drivers of disease progression and recurrence. Additionally, evaluation of rare myeloma cells in blood for disease monitoring has numerous advantages over invasive bone marrow biopsies. To this end, an unbiased method for detecting rare cells and delineating their genomic makeup enables disease detection and monitoring in conditions with low abundant cancer cells. In this study, we applied an enrichment-free four-plex (CD138, CD56, CD45, DAPI) immunofluorescence assay and single-cell DNA sequencing for morphogenomic characterization of plasma cells to detect and delineate common and rare plasma cells and discriminate between normal and malignant plasma cells in paired blood and bone marrow aspirates from five patients with newly diagnosed myeloma (N = 4) and monoclonal gammopathy of undetermined significance (n = 1). Morphological analysis confirms CD138+CD56+ cells in the peripheral blood carry genomic alterations that are clonally identical to those in the bone marrow. A subset of altered CD138+CD56- cells are also found in the peripheral blood consistent with the known variability in CD56 expression as a marker of plasma cell malignancy. Bone marrow tumor clinical cytogenetics is highly correlated with the single-cell copy number alterations of the liquid biopsy rare cells. A subset of rare cells harbors genetic alterations not detected by standard clinical diagnostic methods of random localized bone marrow biopsies. This enrichment-free morphogenomic approach detects and characterizes rare cell populations derived from the liquid biopsies that are consistent with clinical diagnosis and have the potential to extend our understanding of subclonality at the single-cell level in this disease. Assay validation in larger patient cohorts has the potential to offer liquid biopsy for disease monitoring with similar or improved disease detection as traditional blind bone marrow biopsies.

Diverse tumorigenic consequences of human papillomavirus integration in primary oropharyngeal cancers.

Human papillomavirus (HPV) causes 5% of all cancers and frequently integrates into host chromosomes. The HPV oncoproteins E6 and E7 are necessary but insufficient for cancer formation, indicating that additional secondary genetic events are required. Here, we investigate potential oncogenic impacts of virus integration. Analysis of 105 HPV-positive oropharyngeal cancers by whole-genome sequencing detects virus integration in 77%, revealing five statistically significant sites of recurrent integration near genes that regulate epithelial stem cell maintenance (i.e., SOX2, TP63, FGFR, MYC) and immune evasion (i.e., CD274). Genomic copy number hyperamplification is enriched 16-fold near HPV integrants, and the extent of focal host genomic instability increases with their local density. The frequency of genes expressed at extreme outlier levels is increased 86-fold within ±150 kb of integrants. Across 95% of tumors with integration, host gene transcription is disrupted via intragenic integrants, chimeric transcription, outlier expression, gene breaking, and/or de novo expression of noncoding or imprinted genes. We conclude that virus integration can contribute to carcinogenesis in a large majority of HPV-positive oropharyngeal cancers by inducing extensive disruption of host genome structure and gene expression.

Human papillomavirus insertions identify the PIM family of serine/threonine kinases as targetable driver genes in head and neck squamous cell carcinoma.

Human papillomavirus (HPV) insertions in cancer genomes have been linked to various forms of focal genomic instability and altered expression of neighboring genes. Here we tested the hypothesis that investigation of HPV insertions in a head and neck cancer squamous cell carcinoma (HNSCC) cell line would identify targetable driver genes contributing to oncogenesis of other HNSCC. In the cell line UPCI:SCC090 HPV16 integration amplified the PIM1 serine/threonine kinase gene ~16-fold, thereby increasing transcript and protein levels. We used genetic and pharmacological approaches to inhibit PIM kinases in this and other HNSCC cell lines. Knockdown of PIM1 transcripts by transfected short hairpin RNAs reduced UPCI:SCC090 viability. CRISPR/Cas9-mediated mutagenesis of PIM1 caused cell cycle arrest and apoptosis. Pharmacological inhibition of PIM family kinases decreased growth of UPCI:SCC090 and additional HNSCC cell lines in vitro and a xenograft UPCI:SCC090 model in vivo. Based on established interactions between intracellular signaling pathways and relatively high levels of gene expression in almost all HNSCC, we also evaluated combinations of PIM kinase and epidermal growth factor receptor (EGFR) inhibitors. Dual inhibition of these pathways resulted in supra-additive cell death. These data support clinical testing of PIM inhibitors alone or in combination in HNSCC.

Identification of Rare Variants Predisposing to Thyroid Cancer.

Background: Familial non-medullary thyroid cancer (NMTC) accounts for a relatively small proportion of thyroid cancer cases, but it displays strong genetic predisposition. So far, only a few NMTC susceptible genes and low-penetrance variants contributing to NMTC have been described. This study aimed to identify rare germline variants that may predispose individuals to NMTC by sequencing a cohort of 17 NMTC families. Methods: Whole-genome sequencing and genome-wide linkage analysis were performed in 17 NMTC families. MendelScan and BasePlayer were applied to screen germline variants followed by customized filtering. The remaining candidate variants were subsequently validated by Sanger sequencing. A panel of 277 known cancer predisposition genes was also screened in these families. Results: A total of 41 rare coding candidate variants in 40 genes identified by whole-genome sequencing are reported, including 24 missense, five frameshift, five splice change, and seven nonsense variants. Sanger sequencing confirmed all 41 rare variants and proved their co-segregation with NMTC in the extended pedigrees. In silico functional analysis of the candidate genes using Ingenuity Pathway Analysis showed that cancer was the top category of "Diseases and Disorders." Additionally, a targeted search displayed six variants in known cancer predisposition genes, including one frameshift variant and five missense variants. Conclusions: The data identify rare germline variants that may play important roles in NMTC predisposition. It is proposed that in future research including functional characterization, these variants and genes be considered primary candidates for thyroid cancer predisposition.

Human papillomavirus and the landscape of secondary genetic alterations in oral cancers.

Human papillomavirus (HPV) is a necessary but insufficient cause of a subset of oral squamous cell carcinomas (OSCCs) that is increasing markedly in frequency. To identify contributory, secondary genetic alterations in these cancers, we used comprehensive genomics methods to compare 149 HPV-positive and 335 HPV-negative OSCC tumor/normal pairs. Different behavioral risk factors underlying the two OSCC types were reflected in distinctive genomic mutational signatures. In HPV-positive OSCCs, the signatures of APOBEC cytosine deaminase editing, associated with anti-viral immunity, were strongly linked to overall mutational burden. In contrast, in HPV-negative OSCCs, T>C substitutions in the sequence context 5'-ATN-3' correlated with tobacco exposure. Universal expression of HPV E6*1 and E7 oncogenes was a sine qua non of HPV-positive OSCCs. Significant enrichment of somatic mutations was confirmed or newly identified in PIK3CA, KMT2D, FGFR3, FBXW7, DDX3X, PTEN, TRAF3, RB1, CYLD, RIPK4, ZNF750, EP300, CASZ1, TAF5, RBL1, IFNGR1, and NFKBIA Of these, many affect host pathways already targeted by HPV oncoproteins, including the p53 and pRB pathways, or disrupt host defenses against viral infections, including interferon (IFN) and nuclear factor kappa B signaling. Frequent copy number changes were associated with concordant changes in gene expression. Chr 11q (including CCND1) and 14q (including DICER1 and AKT1) were recurrently lost in HPV-positive OSCCs, in contrast to their gains in HPV-negative OSCCs. High-ranking variant allele fractions implicated ZNF750, PIK3CA, and EP300 mutations as candidate driver events in HPV-positive cancers. We conclude that virus-host interactions cooperatively shape the unique genetic features of these cancers, distinguishing them from their HPV-negative counterparts.

Metaplastic breast cancer in a patient with neurofibromatosis type 1 and somatic loss of heterozygosity.

Metaplastic breast carcinoma (MBC) is rare and has a poor prognosis. Here we describe genetic analysis of a 41-yr-old female patient with MBC and neurofibromatosis type I (NF1). She initially presented with pT3N1a, grade 3 MBC, but lung metastases were discovered subsequently. To identify the molecular cause of her NF1, we screened for germline mutations disrupting NF1 or SPRED1, revealing a heterozygous germline single-nucleotide variant (SNV) in exon 21 of NF1 at c.2709G>A, Chr 17: 29556342. By report, this variant disrupts pre-mRNA splicing of NF1 transcripts. No pathogenic mutations were identified in SPRED1 A potential association between MBC and NF1 was reported in eight previous cases, but none underwent detailed genomics analysis. To identify additional candidate germline variants potentially predisposing to MBC, we conducted targeted exome sequencing of 279 established cancer-causing genes in a control blood sample, disclosing four rare SNVs. Analysis of her breast tumor showed markedly altered variant allelic fractions (VAFs) for two (50%) of them, revealing somatic loss of heterozygosity (LOH) at germline SNVs. Of these, only the VAF of the pathogenic SNV in NF1 was increased in the tumor. Tumor sequencing demonstrated five somatic mutations altering TP53, BRCA1, and other genes potentially contributing to cancer formation. Because somatic LOH at certain germline SNVs can enhance their impacts, we conclude that increased allelic imbalance of the pathogenic SNV in NF1 likely contributed to tumorigenesis. Our results highlight a need to assess predisposing genetic factors and LOH that can cause rare, aggressive diseases such as MBC in NF1.

HPV integration hijacks and multimerizes a cellular enhancer to generate a viral-cellular super-enhancer that drives high viral oncogene expression.

Integration of human papillomavirus (HPV) genomes into cellular chromatin is common in HPV-associated cancers. Integration is random, and each site is unique depending on how and where the virus integrates. We recently showed that tandemly integrated HPV16 could result in the formation of a super-enhancer-like element that drives transcription of the viral oncogenes. Here, we characterize the chromatin landscape and genomic architecture of this integration locus to elucidate the mechanisms that promoted de novo super-enhancer formation. Using next-generation sequencing and molecular combing/fiber-FISH, we show that ~26 copies of HPV16 are integrated into an intergenic region of chromosome 2p23.2, interspersed with 25 kb of amplified, flanking cellular DNA. This interspersed, co-amplified viral-host pattern is frequent in HPV-associated cancers and here we designate it as Type III integration. An abundant viral-cellular fusion transcript encoding the viral E6/E7 oncogenes is expressed from the integration locus and the chromatin encompassing both the viral enhancer and a region in the adjacent amplified cellular sequences is strongly enriched in the super-enhancer markers H3K27ac and Brd4. Notably, the peak in the amplified cellular sequence corresponds to an epithelial-cell-type specific enhancer. Thus, HPV16 integration generated a super-enhancer-like element composed of tandem interspersed copies of the viral upstream regulatory region and a cellular enhancer, to drive high levels of oncogene expression.

In Vivo Genome Editing Restores Dystrophin Expression and Cardiac Function in Dystrophic Mice.

RATIONALE: Duchenne muscular dystrophy is a severe inherited form of muscular dystrophy caused by mutations in the reading frame of the dystrophin gene disrupting its protein expression. Dystrophic cardiomyopathy is a leading cause of death in Duchenne muscular dystrophy patients, and currently no effective treatment exists to halt its progression. Recent advancement in genome editing technologies offers a promising therapeutic approach in restoring dystrophin protein expression. However, the impact of this approach on Duchenne muscular dystrophy cardiac function has yet to be evaluated. Therefore, we assessed the therapeutic efficacy of CRISPR (clustered regularly interspaced short palindromic repeats)-mediated genome editing on dystrophin expression and cardiac function in mdx/Utr+/- mice after a single systemic delivery of recombinant adeno-associated virus. OBJECTIVE: To examine the efficiency and physiological impact of CRISPR-mediated genome editing on cardiac dystrophin expression and function in dystrophic mice. METHODS AND RESULTS: Here, we packaged SaCas9 (clustered regularly interspaced short palindromic repeat-associated 9 from Staphylococcus aureus) and guide RNA constructs into an adeno-associated virus vector and systemically delivered them to mdx/Utr+/- neonates. We showed that CRIPSR-mediated genome editing efficiently excised the mutant exon 23 in dystrophic mice, and immunofluorescence data supported the restoration of dystrophin protein expression in dystrophic cardiac muscles to a level approaching 40%. Moreover, there was a noted restoration in the architecture of cardiac muscle fibers and a reduction in the extent of fibrosis in dystrophin-deficient hearts. The contractility of cardiac papillary muscles was also restored in CRISPR-edited cardiac muscles compared with untreated controls. Furthermore, our targeted deep sequencing results confirmed that our adeno-associated virus-CRISPR/Cas9 strategy was very efficient in deleting the ≈23 kb of intervening genomic sequences. CONCLUSIONS: This study provides evidence for using CRISPR-based genome editing as a potential therapeutic approach for restoring dystrophic cardiomyopathy structurally and functionally.

Dynamic silencing of somatic L1 retrotransposon insertions reflects the developmental and cellular contexts of their genomic integration.

BACKGROUND: The ongoing mobilization of mammalian transposable elements (TEs) contributes to natural genetic variation. To survey the epigenetic control and expression of reporter genes inserted by L1 retrotransposition in diverse cellular and genomic contexts, we engineered highly sensitive, real-time L1 retrotransposon reporter constructs. RESULTS: Here we describe different patterns of expression and epigenetic controls of newly inserted sequences retrotransposed by L1 in various somatic cells and tissues including cultured human cancer cells, mouse embryonic stem cells, and tissues of pseudofounder transgenic mice and their progeny. In cancer cell lines, the newly inserted sequences typically underwent rapid transcriptional gene silencing, but they lacked cytosine methylation even after many cell divisions. L1 reporter expression was reversible and oscillated frequently. Silenced or variegated reporter expression was strongly and uniformly reactivated by treatment with inhibitors of histone deacetylation, revealing the mechanism for their silencing. By contrast, de novo integrants retrotransposed by L1 in pluripotent mouse embryonic stem (ES) cells underwent rapid silencing by dense cytosine methylation. Similarly, de novo cytosine methylation also was identified at new integrants when studied in several distinct somatic tissues of adult founder mice. Pre-existing L1 elements in cultured human cancer cells were stably silenced by dense cytosine methylation, whereas their transcription modestly increased when cytosine methylation was experimentally reduced in cells lacking DNA methyltransferases DNMT1 and DNMT3b. As a control, reporter genes mobilized by piggyBac (PB), a DNA transposon, revealed relatively stable and robust expression without apparent silencing in both cultured cancer cells and ES cells. CONCLUSIONS: We hypothesize that the de novo methylation marks at newly inserted sequences retrotransposed by L1 in early pre-implantation development are maintained or re-established in adult somatic tissues. By contrast, histone deacetylation reversibly silences L1 reporter insertions that had mobilized at later timepoints in somatic development and differentiation, e.g., in cancer cell lines. We conclude that the cellular contexts of L1 retrotransposition can determine expression or silencing of newly integrated sequences. We propose a model whereby reporter expression from somatic TE insertions reflects the timing, molecular mechanism, epigenetic controls and the genomic, cellular and developmental contexts of their integration.

HPV positive neuroendocrine cervical cancer cells are dependent on Myc but not E6/E7 viral oncogenes.

Using conditional cell reprogramming, we generated a stable cell culture of an extremely rare and aggressive neuroendocrine cervical cancer. The cultured cells contained HPV-16, formed colonies in soft agar and rapidly produced tumors in immunodeficient mice. The HPV-16 genome was integrated adjacent to the Myc gene, both of which were amplified 40-fold. Analysis of RNA transcripts detected fusion of the HPV/Myc genes, arising from apparent microhomologous recombination. Spectral karyotyping (SKY) and fluorescent-in-situ hybridization (FISH) demonstrated coordinate localization and translocation of the amplified Myc and HPV genes on chromosomes 8 and 21. Similar to the primary tumor, tumor cell cultures expressed very high levels of the Myc protein and, in contrast to all other HPV-positive cervical cancer cell lines, they harbored a gain-of-function mutation in p53 (R273C). Unexpectedly, viral oncogene knockdown had no effect on the growth of the cells, but it did inhibit the proliferation of a conventional HPV-16 positive cervical cancer cell line. Knockdown of Myc, but not the mutant p53, significantly inhibited tumor cell proliferation. On the basis of these data, we propose that the primary driver of transformation in this aggressive cervical cancer is not HPV oncogene expression but rather the overexpression of Myc.

Merkel Cell Polyomavirus Exhibits Dominant Control of the Tumor Genome and Transcriptome in Virus-Associated Merkel Cell Carcinoma.

Merkel cell polyomavirus is the primary etiological agent of the aggressive skin cancer Merkel cell carcinoma (MCC). Recent studies have revealed that UV radiation is the primary mechanism for somatic mutagenesis in nonviral forms of MCC. Here, we analyze the whole transcriptomes and genomes of primary MCC tumors. Our study reveals that virus-associated tumors have minimally altered genomes compared to non-virus-associated tumors, which are dominated by UV-mediated mutations. Although virus-associated tumors contain relatively small mutation burdens, they exhibit a distinct mutation signature with observable transcriptionally biased kataegic events. In addition, viral integration sites overlap focal genome amplifications in virus-associated tumors, suggesting a potential mechanism for these events. Collectively, our studies indicate that Merkel cell polyomavirus is capable of hijacking cellular processes and driving tumorigenesis to the same severity as tens of thousands of somatic genome alterations. IMPORTANCE: A variety of mutagenic processes that shape the evolution of tumors are critical determinants of disease outcome. Here, we sequenced the entire genome of virus-positive and virus-negative primary Merkel cell carcinomas (MCCs), revealing distinct mutation spectra and corresponding expression profiles. Our studies highlight the strong effect that Merkel cell polyomavirus has on the divergent development of viral MCC compared to the somatic alterations that typically drive nonviral tumorigenesis. A more comprehensive understanding of the distinct mutagenic processes operative in viral and nonviral MCCs has implications for the effective treatment of these tumors.

A germline mutation in SRRM2, a splicing factor gene, is implicated in papillary thyroid carcinoma predisposition.

Papillary thyroid carcinoma (PTC) displays strong but so far largely uncharacterized heritability. Here we studied genetic predisposition in a family with six affected individuals. We genotyped all available family members and conducted whole exome sequencing of blood DNA from two affected individuals. Haplotype analysis and other genetic criteria narrowed our list of candidates to a germline variant in the serine/arginine repetitive matrix 2 gene (SRRM2). This heterozygous variant, c.1037C > T (Ser346Phe or S346F; rs149019598) cosegregated with PTC in the family. It was not found in 138 other PTC families. It was found in 7/1,170 sporadic PTC cases and in 0/1,404 controls (p = 0.004). The encoded protein SRRM2 (also called SRm300) is part of the RNA splicing machinery. To evaluate the possibility that the S346F missense mutation affects alternative splicing, we compared RNA-Seq data in leukocytes from three mutation carriers and three controls. Significant differences in alternative splicing were identified for 1,642 exons, of which a subset of 7 exons was verified experimentally. The results confirmed a higher ratio of inclusion of exons in mutation carriers. These data suggest that the S346F mutation in SRRM2 predisposes to PTC by affecting alternative splicing of unidentified downstream target genes.

Multiple functional variants in long-range enhancer elements contribute to the risk of SNP rs965513 in thyroid cancer.

The [A] allele of SNP rs965513 in 9q22 has been consistently shown to be highly associated with increased papillary thyroid cancer (PTC) risk with an odds ratio of ∼1.8 as determined by genome-wide association studies, yet the molecular mechanisms remain poorly understood. Previously, we noted that the expression of two genes in the region, forkhead box E1 (FOXE1) and PTC susceptibility candidate 2 (PTCSC2), is regulated by rs965513 in unaffected thyroid tissue, but the underlying mechanisms were not elucidated. Here, we fine-mapped the 9q22 region in PTC and controls and detected an ∼33-kb linkage disequilibrium block (containing the lead SNP rs965513) that significantly associates with PTC risk. Chromatin characteristics and regulatory element signatures in this block disclosed at least three regulatory elements functioning as enhancers. These enhancers harbor at least four SNPs (rs7864322, rs12352658, rs7847449, and rs10759944) that serve as functional variants. The variant genotypes are associated with differential enhancer activities and/or transcription factor binding activities. Using the chromosome conformation capture methodology, long-range looping interactions of these elements with the promoter region shared by FOXE1 and PTCSC2 in a human papillary thyroid carcinoma cell line (KTC-1) and unaffected thyroid tissue were found. Our results suggest that multiple variants coinherited with the lead SNP and located in long-range enhancers are involved in the transcriptional regulation of FOXE1 and PTCSC2 expression. These results explain the mechanism by which the risk allele of rs965513 predisposes to thyroid cancer.

The methyltransferase Setdb2 mediates virus-induced susceptibility to bacterial superinfection.

Immune responses are tightly regulated to ensure efficient pathogen clearance while avoiding tissue damage. Here we report that Setdb2 was the only protein lysine methyltransferase induced during infection with influenza virus. Setdb2 expression depended on signaling via type I interferons, and Setdb2 repressed expression of the gene encoding the neutrophil attractant CXCL1 and other genes that are targets of the transcription factor NF-κB. This coincided with occupancy by Setdb2 at the Cxcl1 promoter, which in the absence of Setdb2 displayed diminished trimethylation of histone H3 Lys9 (H3K9me3). Mice with a hypomorphic gene-trap construct of Setdb2 exhibited increased infiltration of neutrophils during sterile lung inflammation and were less sensitive to bacterial superinfection after infection with influenza virus. This suggested that a Setdb2-mediated regulatory crosstalk between the type I interferons and NF-κB pathways represents an important mechanism for virus-induced susceptibility to bacterial superinfection.

An antisense promoter in mouse L1 retrotransposon open reading frame-1 initiates expression of diverse fusion transcripts and limits retrotransposition.

Between 6 and 30% of human and mouse transcripts are initiated from transposable elements. However, the promoters driving such transcriptional activity are mostly unknown. We experimentally characterized an antisense (AS) promoter in mouse L1 retrotransposons for the first time, oriented antiparallel to the coding strand of L1 open reading frame-1. We found that AS transcription is mediated by RNA polymerase II. Rapid amplification of cDNA ends cloning mapped transcription start sites adjacent to the AS promoter. We identified >100 novel fusion transcripts, of which many were conserved across divergent mouse lineages, suggesting conservation of potential functions. To evaluate whether AS L1 transcription could regulate L1 retrotransposition, we replaced portions of native open reading frame-1 in donor elements by synonymously recoded sequences. The resulting L1 elements lacked AS promoter activity and retrotransposed more frequently than endogenous L1s. Overexpression of AS L1 transcripts also reduced L1 retrotransposition. This suppression of retrotransposition was largely independent of Dicer. Our experiments shed new light on how AS fusion transcripts are initiated from endogenous L1 elements across the mouse genome. Such AS transcription can contribute substantially both to natural transcriptional variation and to endogenous regulation of L1 retrotransposition.

Genome-wide analysis of HPV integration in human cancers reveals recurrent, focal genomic instability.

Genomic instability is a hallmark of human cancers, including the 5% caused by human papillomavirus (HPV). Here we report a striking association between HPV integration and adjacent host genomic structural variation in human cancer cell lines and primary tumors. Whole-genome sequencing revealed HPV integrants flanking and bridging extensive host genomic amplifications and rearrangements, including deletions, inversions, and chromosomal translocations. We present a model of "looping" by which HPV integrant-mediated DNA replication and recombination may result in viral-host DNA concatemers, frequently disrupting genes involved in oncogenesis and amplifying HPV oncogenes E6 and E7. Our high-resolution results shed new light on a catastrophic process, distinct from chromothripsis and other mutational processes, by which HPV directly promotes genomic instability.