Fragile Sites Research Articles

Chromosomal rearrangements associated with SMC5/6 deficiency in DNA replication.

Completion of DNA replication before chromosome segregation is essential for the stable maintenance of the genome. Under replication stress, DNA synthesis may persist beyond S phase, especially in genomic regions that are difficult to proceed with the replication processes. Incomplete replication in mitosis emerges as non-disjoined segment in mitotic chromosomes leading to anaphase bridges. The resulting chromosome rearrangements are not well characterized, however. Here, we report that incomplete replication due to SMC5/6 deficiency impairs sister chromatid disjunction at difficult-to-replicate regions, including common fragile sites. These non-disjoined regions manifest as cytologically defined symmetric gaps, causing anaphase bridges. These bridges break at the gaps, leading to telomere loss, micronucleation, and fragmentation. Subsequently, fusions between telomere-deficient chromosomes generate complex chromosomal rearrangements, including dicentric chromosomes, suggesting the occurrence of breakage-fusion-bridge cycle. Additionally, chromosomes in micronuclei were pulverized, indicative of chromothripsis. Our findings suggest that incomplete replication facilitates complex chromosomal rearrangements, which may contribute to genomic instability in human cancers.

C9orf72 expansion creates the unstable folate-sensitive fragile site FRA9A

Abstract The hyper-unstable Chr9p21 locus, harbouring the interferon gene cluster, oncogenes and C9orf72, is linked to multiple diseases. C9orf72 (GGGGCC)n expansions (C9orf72Exp) are associated with incompletely penetrant amyotrophic lateral sclerosis, frontotemporal dementia and autoimmune disorders. C9orf72Exp patients display hyperactive cGAS-STING-linked interferon immune and DNA damage responses, but the source of immuno-stimulatory or damaged DNA is unknown. Here, we show C9orf72Exp in pre-symptomatic and ALS-FTD patient cells and brains cause the folate-sensitive chromosomal fragile site, FRA9A. FRA9A centers on > 33kb of C9orf72 as highly-compacted chromatin embedded in an 8.2Mb fragility zone spanning 9p21, encompassing 46 genes, making FRA9A one of the largest fragile sites. C9orf72Exp cells show chromosomal instability, heightened global- and Chr9p-enriched sister-chromatid exchanges, truncated-Chr9s, acentric-Chr9s and Chr9-containing micronuclei, providing endogenous sources of damaged and immunostimulatory DNA. Cells from one C9orf72Exp patient contained highly-rearranged FRA9A-expressing Chr9 with Chr9-wide dysregulated gene expression. Somatic C9orf72Exp repeat instability and chromosomal fragility are sensitive to folate-deficiency. Age-dependent repeat instability, chromosomal fragility, and chromosomal instability can be transferred to CNS and peripheral tissues of transgenic C9orf72Exp mice, implicating C9orf72Exp as the source. Our results highlight unappreciated effects of C9orf72 expansions that trigger vitamin-sensitive chromosome fragility, adding structural variations to the disease-enriched 9p21 locus, and likely elsewhere.

Replication stress induces POLQ-mediated structural variant formation throughout common fragile sites after entry into mitosis

Genomic structural variants (SVs) greatly impact human health, but much is unknown about the mechanisms that generate the largest class of nonrecurrent alterations. Common fragile sites (CFSs) are unstable loci that provide a model for SV formation, especially large deletions, under replication stress. We study SV junction formation as it occurs in human cell lines by applying error-minimized capture sequencing to CFS DNA harvested after low-dose aphidicolin treatment. SV junctions form throughout CFS genes at a 5-fold higher rate after cells pass from G2 into M-phase. Neither SV formation nor CFS expression depend on mitotic DNA synthesis (MiDAS), an error-prone form of replication active at CFSs. Instead, analysis of tens of thousands of de novo SV junctions combined with DNA repair pathway inhibition reveal a primary role for DNA polymerase theta (POLQ)-mediated end-joining (TMEJ). We propose an important role for mitotic TMEJ in nonrecurrent SV formation genome wide.

WWOX tuning of oleic acid signaling orchestrates immunosuppressive macrophage polarization and sensitizes hepatocellular carcinoma to immunotherapy

BackgroundImmune checkpoint inhibitors (ICIs) are therapeutically effective for hepatocellular carcinoma (HCC) but are individually selective. This study examined the role of specific common fragile sites (CFSs) related gene in HCC...

A phenome-wide association study of methylated GC-rich repeats identifies a GCC repeat expansion in AFF3 associated with intellectual disability.

GC-rich tandem repeat expansions (TREs) are often associated with DNA methylation, gene silencing and folate-sensitive fragile sites, and underlie several congenital and late-onset disorders. Through a combination of DNA-methylation profiling and tandem repeat genotyping, we identified 24 methylated TREs and investigated their effects on human traits using phenome-wide association studies in 168,641 individuals from the UK Biobank, identifying 156 significant TRE-trait associations involving 17 different TREs. Of these, a GCC expansion in the promoter of AFF3 was associated with a 2.4-fold reduced probability of completing secondary education, an effect size comparable to several recurrent pathogenic microdeletions. In a cohort of 6,371 probands with neurodevelopmental problems of suspected genetic etiology, we observed a significant enrichment of AFF3 expansions compared with controls. With a population prevalence that is at least fivefold higher than the TRE that causes fragile X syndrome, AFF3 expansions represent a major cause of neurodevelopmental delay.

Polyubiquitinated PCNA triggers SLX4-mediated break-induced replication in alternative lengthening of telomeres (ALT)cancer cells.

Replication stresses are the major source of break-induced replication (BIR). Here, we show that in alternative lengthening of telomeres (ALT) cells, replication stress-induced polyubiquitinated proliferating cell nuclear antigen (PCNA) (polyUb-PCNA) triggers BIR at telomeres and the common fragile site (CFS). Consistently, depleting RAD18, a PCNA ubiquitinating enzyme, reduces the occurrence of ALT-associated promyelocytic leukemia (PML) bodies (APBs) and mitotic DNA synthesis at telomeres and CFS, both of which are mediated by BIR. In contrast, inhibiting ubiquitin-specific protease 1 (USP1), an Ub-PCNA deubiquitinating enzyme, results in an increase in the above phenotypes in a RAD18- and UBE2N (the PCNA polyubiquitinating enzyme)-dependent manner. Furthermore, deficiency of ATAD5, which facilitates USP1 activity and unloads PCNAs, augments recombination-associated phenotypes. Mechanistically, telomeric polyUb-PCNA accumulates SLX4, a nuclease scaffold, at telomeres through its ubiquitin-binding domain and increases telomere damage. Consistently, APB increase induced by Ub-PCNA depends on SLX4 and structure-specific endonucleases. Taken together, our results identified the polyUb-PCNA-SLX4 axis as a trigger for directing BIR.

High Prevalence of Chromosomal Rearrangements and LINE Retrotranspositions Detected in Formalin-Fixed, Paraffin-Embedded Colorectal Cancer Tissue

Structural variants (SVs) caused by chromosomal rearrangements in common fragile sites (CFS) or LINE retrotranspositions are highly prevalent in colorectal cancer (CRC). However, methodology for targeted detection of these SVs is lacking. We here report the use of formalin-fixed paraffin-embedded targeted locus capture (FFPE-TLC) sequencing as a novel technology for targeted detection of tumor-specific SVs. Analysis of 29 FFPE colorectal tumor and 8 matched normal samples revealed tumor-specific SVs in 24 patients (83%) with a median of 2 SVs per patient (range 1-21). 104 SVs were found in the CFS-associated genes MACROD2, PRKN, FHIT and WWOX in 18 patients (62%) and 39 SVs caused by three LINE transposable elements were found in 15 patients (52%). Tumor-specificity of SVs was independently verified by droplet digital PCR of tumor tissue DNA and their applicability as plasma circulating tumor DNA biomarkers was demonstrated. We conclude that FFPE-TLC sequencing enables the detection of tumor-specific SVs caused by chromosomal rearrangements and LINE retrotranspositions in FFPE tissue. Therefore, FFPE-TLC sequencing facilitates the investigation of the biological and clinical impact of SVs using FFPE material from (retrospective) cohorts of cancer patients and has potential clinical applicability to detect SV biomarkers in the routine molecular diagnostics setting.

Innovative Tools for DNA Topology Probing in Human Cells Reveal a Build-Up of Positive Supercoils Following Replication Stress at Telomeres and at the FRA3B Fragile Site.

Linear unconstrained DNA cannot harbor supercoils since these supercoils can diffuse and be eliminated by free rotation of the DNA strands at the end of the molecule. Mammalian telomeres, despite constituting the ends of linear chromosomes, can hold supercoils and be subjected to topological stress. While negative supercoiling was previously observed, thus proving the existence of telomeric topological constraints, positive supercoils were never probed due to the lack of an appropriate tool. Indeed, the few tools available currently could only investigate unwound (Trioxsalen) or overwound (GapR) DNA topology (variations in twist) but not the variations in writhe (supercoils and plectonemes). To address this question, we have designed innovative tools aimed at analyzing both positive and negative DNA writhe in cells. Using them, we could observe the build-up of positive supercoils following replication stress and inhibition of Topoisomerase 2 on telomeres. TRF2 depletion caused both telomere relaxation and an increase in positive supercoils while the inhibition of Histone Deacetylase I and II by TSA only caused telomere relaxation. Moving outside telomeres, we also observed a build-up of positive supercoils on the FRA3B fragile site following replication stress, suggesting a topological model of DNA fragility for this site.

Evaluation of the gene fusion landscape in early onset sporadic rectal cancer reveals association with chromatin architecture and genome stability.

Gene fusions represent a distinct class of structural variants identified frequently in cancer genomes across cancer types. Several gene fusions exhibit gain of oncogenic function and thus have been the focus of development of efficient targeted therapies. However, investigation of fusion landscape in early-onset sporadic rectal cancer, a poorly studied colorectal cancer subtype prevalent in developing countries, has not been performed. Here, we present a comprehensive landscape of gene fusions in EOSRC and CRC using patient derived tumor samples and data from The Cancer Genome Atlas, respectively. Gene Ontology analysis revealed enrichment of unique biological process terms associated with 5'- and 3'- fusion partner genes. Extensive network analysis highlighted genes exhibiting significant promiscuity in fusion formation and their association with chromosome fragile sites. Investigation of fusion formation in the context of global chromatin architecture unraveled a novel mode of gene activation that arose from fusion between genes located in orthogonal chromatin compartments. The study provides novel evidence linking fusions to genome stability and architecture and unearthed a hitherto unidentified mode of gene activation in cancer.

TDP1 phosphorylation by CDK1 in mitosis promotes MUS81-dependent repair of trapped Top1-DNA covalent complexes

Topoisomerase 1 (Top1) controls DNA topology, relieves DNA supercoiling during replication and transcription, and is critical for mitotic progression to the G1 phase. Tyrosyl-DNA phosphodiesterase 1 (TDP1) mediates the removal of trapped Top1-DNA covalent complexes (Top1cc). Here, we identify CDK1-dependent phosphorylation of TDP1 at residue S61 during mitosis. A TDP1 variant defective for S61 phosphorylation (TDP1-S61A) is trapped on the mitotic chromosomes, triggering DNA damage and mitotic defects. Moreover, we show that Top1cc repair in mitosis occurs via a MUS81-dependent DNA repair mechanism. Replication stress induced by camptothecin or aphidicolin leads to TDP1-S61A enrichment at common fragile sites, which over-stimulates MUS81-dependent chromatid breaks, anaphase bridges, and micronuclei, ultimately culminating in the formation of 53BP1 nuclear bodies during G1 phase. Our findings provide new insights into the cell cycle-dependent regulation of TDP1 dynamics for the repair of trapped Top1-DNA covalent complexes during mitosis that prevents genomic instability following replication stress.

A post-irradiation-induced replication stress promotes RET proto-oncogene breakage.

Ionizing radiation generates genomic instability by promoting the accumulation of chromosomal rearrangements. The oncogenic translocation RET/PTC1 is present in more than 70% of radiation-induced thyroid cancers. Both RET and CCDC6, the genes implicated in RET/PTC1, are found within common fragile sites - chromosomal regions prone to DNA breakage during slight replication stress. Given that irradiated cells become more susceptible to genomic destabilization due to the accumulation of replication-stress-related double-strand breaks (DSBs), we explored whether RET and CCDC6 exhibit DNA breakage under replicative stress several days post-irradiation of thyroid cells. We analyzed the dynamic of DNA replication in human thyroid epithelial cells (HThy-ori-3.1) 4 days post a 5-Gy exposure using molecular DNA combing. The DNA replication schedule was evaluated through replication-timing experiments. We implemented a ChIP-qPCR assay to determine whether the RET and CCDC6 genes break following irradiation. Our study indicates that replicative stress, occurring several days post-irradiation in thyroid cells, primarily causes DSBs in the RET gene. We discovered that both the RET and CCDC6 genes undergo late replication in thyroid cells. However, only RET's replication rate is notably delayed after irradiation. The findings suggest that post-irradiation in the RET gene causes a breakage in the replication fork, which could potentially invade another genomic area, including CCDC6. As a result, this could greatly contribute to the high prevalence of chromosomal RET/PTC rearrangements seen in patients exposed to external radiation.

The long noncoding RNA MANCR to promote and stabilize triple negative breast cancer.

e13150 Background: Cancer remains a leading cause of death worldwide, despite many advances in cancer treatment over the past decade. Triple Negative Breast Cancer (TNBC) is among the highest with a poor prognosis. Early diagnosis and modern treatment regimens for breast cancer (BC) have led to a high disease-free survival rate. Still, survival decreases significantly for patients with metastatic disease. Long non-coding RNAs are among a recent class of epigenetic regulators that function in the nucleus to support the stability of the cells and maintain the fidelity of chromatin interactions. Our laboratory discovered the long noncoding RNA MANCR (LINC00704) as being upregulated in human BC tumors and this is indicative of a worse survival rate than in patients with low MANCR-expressing tumors. MANCR is aberrantly expressed in TNBC cells, and these cells are highly dependent on MANCR to retain their tumorigenic characteristics. Methods: Our functional in vitro studies in MDA-MB-231 TNBC cells used GapmerRs to knockdown MANCR. For in vivo studies, TNBC cells were injected into the 4th mammary fat pad of mice, allowed initial tumor formation, and then treated the animals with 2 nmol/g of a negative control or MANCR targeting GapmeR. We also identified genome interaction sites at single nucleotide resolution by chromatin isolation by RNA purification (ChIRP)-seq to determine the mechanism of MANCR activity in TNBC cells. Results: We found that MANCR knockdown promotes DNA damage and decreases cell proliferation, migration, anchorage-independent colony formation, transwell invasion, and overall survival. Additionally, targeting MANCR in animals with existing tumors drastically inhibited tumor growth and the end-point tumor mass over time. The MANCR GapmeR treatment also inhibited the ability of TNBC cells to circulate and disseminate to distant organs. After performing ChIRP-seq in the MDA-MB-231 cells, we identified 988 genome-wide binding sites that exhibit MANCR interactions, of which 79% are intergenic and 21% are genic regions. Furthermore, 127 MANCR ChIRP peaks were found to overlap with fragile sites in the genome, indicating MANCR provides stability to these sites. Conclusions: These data suggest that targeting MANCR has therapeutic potential for patients with “MANCR-high” TNBC tumors by disrupting genome stability. Indeed, our in vivo studies demonstrate that “MANCR-high” TNBC tumors require MANCR to rapidly grow and promote disease progression. Lastly, many of the MANCR-chromatin interactions identified were found in intergenic regions and overlap with fragile sites within the genome. Collectively these data strongly indicate that MANCR stabilizes the TNBC genome, and disrupting genome stability by targeting MANCR has therapeutic potential.

Genome-wide identification of replication fork stalling/pausing sites and the interplay between RNA Pol II transcription and DNA replication progression

BackgroundDNA replication progression can be affected by the presence of physical barriers like the RNA polymerases, leading to replication stress and DNA damage. Nonetheless, we do not know how transcription influences overall DNA replication progression.ResultsTo characterize sites where DNA replication forks stall and pause, we establish a genome-wide approach to identify them. This approach uses multiple timepoints during S-phase to identify replication fork/stalling hotspots as replication progresses through the genome. These sites are typically associated with increased DNA damage, overlapped with fragile sites and with breakpoints of rearrangements identified in cancers but do not overlap with replication origins. Overlaying these sites with a genome-wide analysis of RNA polymerase II transcription, we find that replication fork stalling/pausing sites inside genes are directly related to transcription progression and activity. Indeed, we find that slowing down transcription elongation slows down directly replication progression through genes. This indicates that transcription and replication can coexist over the same regions. Importantly, rearrangements found in cancers overlapping transcription-replication collision sites are detected in non-transformed cells and increase following treatment with ATM and ATR inhibitors. At the same time, we find instances where transcription activity favors replication progression because it reduces histone density.ConclusionsAltogether, our findings highlight how transcription and replication overlap during S-phase, with both positive and negative consequences for replication fork progression and genome stability by the coexistence of these two processes.

A Unifying Hypothesis for the Genome Dynamics Proposed to Underlie Neuropsychiatric Phenotypes.

The sheer number of gene variants and the extent of the observed clinical and molecular heterogeneity recorded in neuropsychiatric disorders (NPDs) could be due to the magnified downstream effects initiated by a smaller group of genomic higher-order alterations in response to endogenous or environmental stress. Chromosomal common fragile sites (CFS) are functionally linked with microRNAs, gene copy number variants (CNVs), sub-microscopic deletions and duplications of DNA, rare single-nucleotide variants (SNVs/SNPs), and small insertions/deletions (indels), as well as chromosomal translocations, gene duplications, altered methylation, microRNA and L1 transposon activity, and 3-D chromosomal topology characteristics. These genomic structural features have been linked with various NPDs in mostly isolated reports and have usually only been viewed as areas harboring potential candidate genes of interest. The suggestion to use a higher level entry point (the 'fragilome' and associated features) activated by a central mechanism ('stress') for studying NPD genetics has the potential to unify the existing vast number of different observations in this field. This approach may explain the continuum of gene findings distributed between affected and unaffected individuals, the clustering of NPD phenotypes and overlapping comorbidities, the extensive clinical and molecular heterogeneity, and the association with certain other medical disorders.

Chronic replication stress invokes mitochondria dysfunction via impaired parkin activity

Replication stress is a major contributor to tumorigenesis because it provides a source of chromosomal rearrangements via recombination events. PARK2, which encodes parkin, a regulator of mitochondrial homeostasis, is located on one of the common fragile sites that are prone to rearrangement by replication stress, indicating that replication stress may potentially impact mitochondrial homeostasis. Here, we show that chronic low-dose replication stress causes a fixed reduction in parkin expression, which is associated with mitochondrial dysfunction, indicated by an increase in mtROS. Consistent with the major role of parkin in mitophagy, reduction in parkin protein expression was associated with a slight decrease in mitophagy and changes in mitochondrial morphology. In contrast, cells expressing ectopic PARK2 gene does not show mtROS increases and changes in mitochondrial morphology even after exposure to chronic replication stress, suggesting that intrinsic fragility at PARK2 loci associated with parkin reduction is responsible for mitochondrial dysfunction caused by chronic replication stress. As endogenous replication stress and mitochondrial dysfunction are both involved in multiple pathophysiology, our data support the therapeutic development of recovery of parkin expression in human healthcare.

An application of logistic regression model to the study of constitutional imbalances in human chromosomal fragile sites

An application of logistic regression model to the study of constitutional imbalances in human chromosomal fragile sites

Replication of [AT/TA]25 Microsatellite Sequences by Human DNA Polymerase δ Holoenzymes Is Dependent on dNTP and RPA Levels.

Fragile sites are unstable genomic regions that are prone to breakage during stressed DNA replication. Several common fragile sites (CFS) contain A+T-rich regions including perfect [AT/TA] microsatellite repeats that may collapse into hairpins when in single-stranded DNA (ssDNA) form and coincide with chromosomal hotspots for breakage and rearrangements. While many factors contribute to CFS instability, evidence exists for replication stalling within [AT/TA] microsatellite repeats. Currently, it is unknown how stress causes replication stalling within [AT/TA] microsatellite repeats. To investigate this, we utilized FRET to characterize the structures of [AT/TA]25 sequences and also reconstituted lagging strand replication to characterize the progression of pol δ holoenzymes through A+T-rich sequences. The results indicate that [AT/TA]25 sequences adopt hairpins that are unwound by the major ssDNA-binding complex, RPA, and the progression of pol δ holoenzymes through A+T-rich sequences saturated with RPA is dependent on the template sequence and dNTP concentration. Importantly, the effects of RPA on the replication of [AT/TA]25 sequences are dependent on dNTP concentration, whereas the effects of RPA on the replication of A+T-rich, nonstructure-forming sequences are independent of dNTP concentration. Collectively, these results reveal complexities in lagging strand replication and provide novel insights into how [AT/TA] microsatellite repeats contribute to genome instability.

Abstract 401: The impact of DNA structural variants on PRKN gene function and colorectal cancer pathogenesis

Abstract Background: Cancer is caused by somatic DNA alterations, such as single nucleotide variants, somatic copy number alterations and structural variants (SVs). PRKN is among the genes most frequently affected by SVs in colorectal cancer (CRC), with 37% of primary cancers and 56% of metastatic lesions having focal deletions in PRKN. PRKN encodes the E3 ubiquitin ligase Parkin, which regulates the biological activity and degradation of substrates via post-translational modifications. The genomic location of PRKN is within a common fragile site. Therefore, it remains unclear whether the high frequency of SVs in PRKN are passenger events due to genomic instability or driver events of CRC pathogenesis. Aim: Our study aims to investigate the impact of SVs on the functionality of Parkin in CRC pathogenesis. Methods: We used publicly available DNA deep whole genome sequencing (WGS) and tumor-matched RNA sequencing data (n=394) to estimate the impact of SVs in PRKN on tumor biology. For this, the transcriptome-wide p-value distribution generated by differential expression analysis was used as a measure of biological impact. The HCT116-MLH1-repaired cell line and patient-derived adenoma organoids were genetically modified using CRISPR/Cas9 to mimic the SVs observed in patient material and analyze their impact. RNA-sequencing, whole proteome and ubiquitin remnant-specific proteomics were used to assess the effect of SVs on gene function and identify altered biological processes. Results: Analysis of WGS and RNA sequencing of patient material showed that SVs in PRKN have a large impact on CRC biology, comparable to the impact of mutations in TP53 (impact score 0.379 vs 0.332). The most pronounced effects are seen in the ubiquitin-like, RING0 and RING1 domains. These domains are important for recognition and/or binding to ubiquitin, E2 ubiquitin-conjugating enzymes and substrates. In the HCT116-MLH1-repaired CRC cell line, hemizygous introduction of an SV in PRKN had minimal effects on genome-wide gene expression. In contrast, homozygous deletion altered the expression of nearly 2000 genes. Analyses of ubiquitin motif remnants showed differential ubiquitination of proteins involved in DNA repair following the introduction of an SV in PRKN in the CRC cell line. These proteins include TOP2A (Log2 fold change (L2FC) = -3.69, P = 5.71−5), RAD18 (L2FC = -4.55, P = 1.35−4) and Spartan (SPRTN; L2FC = -2.57, P = 7.43−4). Conclusion: SVs in PRKN have a high impact on tumor biology and result in aberrant gene expression and post-translational modifications of proteins related to DNA repair. Further research will expand to the patient-derived adenoma organoids and focus on understanding if these recurrent somatic alterations contribute to genomic instability and the development of CRC pathogenesis. Citation Format: Elisabeth J. van Bree, Soufyan Lakbir, Carmen Rubio-Alarcón, Renske de Wit, Anne Bolijn, Pien Delis-van Diemen, Marianne Tijssen, Ellen Stelloo, Harma Feitsma, Gerrit A. Meijer, Beatriz Carvalho, Sanne Abeln, Remond J. Fijneman. The impact of DNA structural variants on PRKN gene function and colorectal cancer pathogenesis [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 401.

Mapping the landscape of HPV integration and characterising virus and host genome interactions in HPV-positive oropharyngeal squamous cell carcinoma.

Human papillomavirus (HPV) integration into the host genome is an important factor in HPV(+)OPSCC carcinogenesis, in conjunction with HPV oncoproteins E6/E7. However, a well-studied investigation about virus-host interaction still needs to be completed. Our objective is to characterise HPV integration to investigate potential mechanisms of tumourigenesis independent of E6/E7 oncoproteins. High-throughput viral integration detection was performed on 109 HPV(+)OPSCC tumours with relevant clinicopathological information. Of these tumours, 38 tumours underwent targeted gene sequencing, 29 underwent whole exome sequencing and 26 underwent RNA sequencing. HPV integration was detected in 94% of tumours (with a mean integration count of 337). Tumours occurring at the tonsil/oropharyngeal wall that exhibit higher PD-L1 expression demonstrated increased integration sites (p=.024). HPV exhibited a propensity for integration at genomic sites located within specific fragile sites (FRA19A) or genes associated with functional roles such as cell proliferation and differentiation (PTEN, AR), immune evasion (CD274) and glycoprotein biosynthesis process (FUT8). The viral oncogenes E2, E4, E6 and E7 tended to remain intact. HPV fragments displayed enrichment within host copy number variation (CNV) regions. However, insertions into genes related to altered homologous recombination repair were infrequent. Genes with integration had distinct expression levels. Fifty-nine genes whose expression level was affected by viral integration were identified, for example, EPHB1, which was reported to be involved in cellular protein metabolic process. HPV can promote oncogenesis through recurrent integration into functional host genome regions, leading to subsequent genomic aberrations and gene expression disruption. This study characterises viral integrations and virus-host interactions, enhancing our understanding of HPV-related carcinogenesis mechanisms.

Cell-Free Genic Extrachromosomal Circular DNA Profiles of DNase Knockouts Associated with Systemic Lupus Erythematosus and Relation with Common Fragile Sites.

Cell-free extrachromosomal circular DNA (cf-eccDNA) has been proposed as a promising early biomarker for disease diagnosis, progression and drug response. Its established biomarker features are changes in the number and length distribution of cf-eccDNA. Another novel promising biomarker is a set of eccDNA excised from a panel of genes specific to a condition compared to a control. Deficiencies in two endonucleases that specifically target DNA, Dnase1 and Dnase1l3, are associated with systemic lupus erythematosus (SLE). To study the genic eccDNA profiles in the case of their deficiencies, we mapped sequenced eccDNA data from plasma, liver and buffy coat from Dnase1 and Dnase1l3 knockouts (KOs), and wild type controls in mouse. Next, we performed an eccDNA differential analysis between KO and control groups using our DifCir algorithm. We found a specific genic cf-eccDNA fingerprint of the Dnase1l3 group compared to the wild type controls involving 131 genes; 26% of them were associated with human chromosomal fragile sites (CFSs) and with a statistically significant enrichment of CFS-associated genes. We found six genes in common with the genic cf-eccDNA profile of SLE patients with DNASE1L3 deficiency, namely Rorb, Mvb12b, Osbpl10, Fto, Tnik and Arhgap10; all of them were specific and present in all human plasma samples, and none of them were associated with CFSs. A not so distinctive genic cf-eccDNA difference involving only seven genes was observed in the case of the Dnase1 group compared to the wild type. In tissue-liver and buffy coat-we did not detect the same genic eccDNA difference observed in the plasma samples. These results point to a specific role of a set of genic eccDNA in plasma from DNase KOs, as well as a relation with CFS genes, confirming the promise of the genic cf-eccDNA in studying diseases and the need for further research on the relationship between eccDNA and CFSs.