The chromatin insulator CTCF regulates HPV18 transcript splicing and differentiation-dependent late gene expression.

The ubiquitous host protein, CCCTC-binding factor (CTCF), is an essential regulator of cellular transcription and functions to maintain epigenetic boundaries, stabilise chromatin loops and regulate splicing of alternative exons. We have previously demonstrated that CTCF binds to the E2 open reading frame (ORF) of human papillomavirus (HPV) 18 and functions to repress viral oncogene expression in undifferentiated keratinocytes by co-ordinating an epigenetically repressed chromatin loop within HPV episomes. Keratinocyte differentiation disrupts CTCF-dependent chromatin looping of HPV18 episomes promoting induction of enhanced viral oncogene expression. To further characterise CTCF function in HPV transcription control we utilised direct, long-read Nanopore RNA-sequencing which provides information on the structure and abundance of full-length transcripts. Nanopore analysis of primary human keratinocytes containing HPV18 episomes before and after synchronous differentiation allowed quantification of viral transcript species, including the identification of low abundance novel transcripts. Comparison of transcripts produced in wild type HPV18 genome-containing cells to those identified in CTCF-binding deficient genome-containing cells identifies CTCF as a key regulator of differentiation-dependent late promoter activation, required for efficient E1^E4 and L1 protein expression. Furthermore, our data show that CTCF binding at the E2 ORF promotes usage of the downstream weak splice donor (SD) sites SD3165 and SD3284, to the dominant E4 splice acceptor site at nucleotide 3434. These findings demonstrate that in the HPV life cycle both early and late virus transcription programmes are facilitated by recruitment of CTCF to the E2 ORF.

SARS-CoV-2 Testing in the Community: Testing Positive Samples with the TaqMan SARS-CoV-2 Mutation Panel To Find Variants in Real Time.

Genome sequencing is a powerful tool for identifying SARS-CoV-2 variant lineages; however, there can be limitations due to sequence dropout when used to identify specific key mutations. Recently, ThermoFisher Scientific has developed genotyping assays to help bridge the gap between testing capacity and sequencing capability to generate real-time genotyping results based on specific variants. Over a 6-week period during the months of April and May 2021, we set out to assess the ThermoFisher TaqMan mutation panel genotyping assay, initially for three mutations of concern and then for an additional two mutations of concern, against SARS-CoV-2-positive clinical samples and the corresponding COVID-19 Genomics UK Consortium (COG-UK) sequencing data. We demonstrate that genotyping is a powerful in-depth technique for identifying specific mutations, is an excellent complement to genome sequencing, and has real clinical health value potential, allowing laboratories to report and take action on variants of concern much more quickly.

The molecular landscape of well differentiated retroperitoneal liposarcoma.

Well differentiated liposarcoma (WD-LPS) is a relatively rare tumour, with fewer than 50 cases occurring per year in the UK. These tumours are both chemotherapy- and radiotherapy-resistant and present a significant treatment challenge requiring radical surgery. Little is known of the molecular landscape of these tumours and no current targets for molecular therapy exist. We aimed to carry out a comprehensive molecular characterisation of WD-LPS via whole genome sequencing, RNA sequencing, and methylation array analysis. A recurrent mutation within exon 1 of FOXD4L3 was observed (chr9:70,918,189A>T; c.322A>T; p.Lys108Ter). Recurrent mutations were also observed in Wnt signalling, immunity, DNA repair, and hypoxia-associated genes. Recurrent amplification of HGMA2 was observed, although this was in fact part of a general amplification of the region around this gene. Recurrent gene fusions in HGMA2, SDHA, TSPAN31, and MDM2 were also observed as well as consistent rearrangements between chromosome 6 and chromosome 12. Our study has demonstrated a recurrent mutation within FOXD4L3, which shows evidence of interaction with the PAX pathway to promote tumourigenesis. © 2021 The Authors. The Journal of Pathology published by John Wiley & Sons, Ltd. on behalf of The Pathological Society of Great Britain and Ireland.

Disruption of CTCF-YY1-dependent looping of the human papillomavirus genome activates differentiation-induced viral oncogene transcription.

The complex life cycle of oncogenic human papillomavirus (HPV) initiates in undifferentiated basal epithelial keratinocytes where expression of the E6 and E7 oncogenes is restricted. Upon epithelial differentiation, E6/E7 transcription is increased through unknown mechanisms to drive cellular proliferation required to support virus replication. We report that the chromatin-organising CCCTC-binding factor (CTCF) promotes the formation of a chromatin loop in the HPV genome that epigenetically represses viral enhancer activity controlling E6/E7 expression. CTCF-dependent looping is dependent on the expression of the CTCF-associated Yin Yang 1 (YY1) transcription factor and polycomb repressor complex (PRC) recruitment, resulting in trimethylation of histone H3 at lysine 27. We show that viral oncogene up-regulation during cellular differentiation results from YY1 down-regulation, disruption of viral genome looping, and a loss of epigenetic repression of viral enhancer activity. Our data therefore reveal a key role for CTCF-YY1-dependent looping in the HPV life cycle and identify a regulatory mechanism that could be disrupted in HPV carcinogenesis.

LARP1 isoform expression in human cancer cell lines.

LARP1 is an oncogenic RNA-binding protein required for ribosome biogenesis and cancer cell survival. From published in vitro studies, there is disparity over which of two different LARP1 protein isoforms (termed the long LI-LARP1 and short SI-LARP1) is the canonical. Here, after conducting a series of biochemical and cellular assays, we conclude that LI-LARP1 (NM_033551.3 > NP_056130.2) is the dominantly expressed form. We observe that SI-LARP1 (NM_015315.5> NP_056130.2) is epigenetically repressed and that this repression is evolutionarily conserved in all but a small subclade of mammalian species. As with other LARP family members, there are multiple potential LARP1 mRNA isoforms that appear to be censored within the nucleus. The capacity of the cell to modulate splicing and expression of these apparently 'redundant' mRNAs hints at contextually specific mechanisms of LARP1 expression.

Targeted deep sequencing of urothelial bladder cancers and associated urinary DNA: a 23-gene panel with utility for non-invasive diagnosis and risk stratification.

OBJECTIVES: To develop a focused panel of somatic mutations (SMs) present in the majority of urothelial bladder cancers (UBCs), to investigate the diagnostic and prognostic utility of this panel, and to compare the identification of SMs in urinary cell-pellet (cp)DNA and cell-free (cf)DNA as part of the development of a non-invasive clinical assay. PATIENTS AND METHODS: A panel of SMs was validated by targeted deep-sequencing of tumour DNA from 956 patients with UBC. In addition, amplicon and capture-based targeted sequencing measured mutant allele frequencies (MAFs) of SMs in 314 urine cpDNAs and 153 urine cfDNAs. The association of SMs with grade, stage, and clinical outcomes was investigated by univariate and multivariate Cox models. Concordance between SMs detected in tumour tissue and cpDNA and cfDNA was assessed. RESULTS: The panel comprised SMs in 23 genes: TERT (promoter), FGFR3, PIK3CA, TP53, ERCC2, RHOB, ERBB2, HRAS, RXRA, ELF3, CDKN1A, KRAS, KDM6A, AKT1, FBXW7, ERBB3, SF3B1, CTNNB1, BRAF, C3orf70, CREBBP, CDKN2A, and NRAS; 93.5-98.3% of UBCs of all grades and stages harboured ≥1 SM (mean: 2.5 SMs/tumour). RAS mutations were associated with better overall survival (P = 0.04). Mutations in RXRA, RHOB and TERT (promoter) were associated with shorter time to recurrence (P < 0.05). MAFs in urinary cfDNA and cpDNA were highly correlated; using a capture-based approach, >94% of tumour SMs were detected in both cpDNA and cfDNA. CONCLUSIONS: SMs are reliably detected in urinary cpDNA and cfDNA. The technical capability to identify very low MAFs is essential to reliably detect UBC, regardless of the use of cpDNA or cfDNA. This 23-gene panel shows promise for the non-invasive diagnosis and risk stratification of UBC.

A multi-centre clinico-genetic analysis of the VPS35 gene in Parkinson disease indicates reduced penetrance for disease-associated variants.

BACKGROUND: Two recent studies identified a mutation (p.Asp620Asn) in the vacuolar protein sorting 35 gene as a cause for an autosomal dominant form of Parkinson disease . Although additional missense variants were described, their pathogenic role yet remains inconclusive. METHODS AND RESULTS: We performed the largest multi-center study to ascertain the frequency and pathogenicity of the reported vacuolar protein sorting 35 gene variants in more than 15,000 individuals worldwide. p.Asp620Asn was detected in 5 familial and 2 sporadic PD cases and not in healthy controls, p.Leu774Met in 6 cases and 1 control, p.Gly51Ser in 3 cases and 2 controls. Overall analyses did not reveal any significant increased risk for p.Leu774Met and p.Gly51Ser in our cohort. CONCLUSIONS: Our study apart from identifying the p.Asp620Asn variant in familial cases also identified it in idiopathic Parkinson disease cases, and thus provides genetic evidence for a role of p.Asp620Asn in Parkinson disease in different populations worldwide.

Patient Derived Organoids Confirm That PI3K/AKT Signalling Is an Escape Pathway for Radioresistance and a Target for Therapy in Rectal Cancer.

Objectives: Partial or total resistance to preoperative chemoradiotherapy occurs in more than half of locally advanced rectal cancer patients. Several novel or repurposed drugs have been trialled to improve cancer cell sensitivity to radiotherapy, with limited success. We aimed to understand the mechanisms of resistance to chemoradiotherapy in rectal cancer using patient derived organoid models. Design: To understand the mechanisms underlying this resistance, we compared the pre-treatment transcriptomes of patient-derived organoids (PDO) with measured radiotherapy sensitivity to identify biological pathways involved in radiation resistance coupled with single cell sequencing, genome wide CRISPR-Cas9 and targeted drug screens. Results: RNA sequencing enrichment analysis revealed upregulation of PI3K/AKT/mTOR and epithelial mesenchymal transition pathway genes in radioresistant PDOs. Single-cell sequencing of pre & post-irradiation PDOs showed mTORC1 and PI3K/AKT upregulation, which was confirmed by a genome-wide CRSIPR-Cas9 knockout screen using irradiated colorectal cancer (CRC) cell lines. We then tested the efficiency of dual PI3K/mTOR inhibitors in improving cancer cell sensitivity to radiotherapy. After irradiation, significant AKT phosphorylation was detected (p=0.027) which was abrogated with dual PI3K/mTOR inhibitors and lead to significant radiosensitisation of the HCT116 cell line and radiation resistant PDO lines. Conclusions: The PI3K/AKT/mTOR pathway upregulation contributes to radioresistance and its targeted pharmacological inhibition leads to significant radiosensitisation in CRC organoids, making it a potential target for clinical trials.

Recurrent SARS-CoV-2 mutations in immunodeficient patients.

Long-term severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections in immunodeficient patients are an important source of variation for the virus but are understudied. Many case studies have been published which describe one or a small number of long-term infected individuals but no study has combined these sequences into a cohesive dataset. This work aims to rectify this and study the genomics of this patient group through a combination of literature searches as well as identifying new case series directly from the COVID-19 Genomics UK (COG-UK) dataset. The spike gene receptor-binding domain and N-terminal domain (NTD) were identified as mutation hotspots. Numerous mutations associated with variants of concern were observed to emerge recurrently. Additionally a mutation in the envelope gene, T30I was determined to be the second most frequent recurrently occurring mutation arising in persistent infections. A high proportion of recurrent mutations in immunodeficient individuals are associated with ACE2 affinity, immune escape, or viral packaging optimisation. There is an apparent selective pressure for mutations that aid cell-cell transmission within the host or persistence which are often different from mutations that aid inter-host transmission, although the fact that multiple recurrent de novo mutations are considered defining for variants of concern strongly indicates that this potential source of novel variants should not be discounted.

Genetic modifiers in rare disorders: the case of fragile X syndrome.

Methods employed in genome-wide association studies are not feasible ways to explore genotype-phenotype associations in rare disorders due to limited statistical power. An alternative approach is to examine relationships among specific single nucleotide polymorphisms (SNPs), selected a priori, and behavioural characteristics. Here, we adopt this strategy to examine relationships between three SNPs (5-HTTLPR, MAOA, COMT) and specific clinically-relevant behaviours that are phenotypic of fragile X syndrome (FXS) but vary in severity and frequency across individuals. Sixty-four males with FXS participated in the current study. Data from standardised informant measures of challenging behaviour (defined as physical aggression, property destruction, stereotyped behaviour, and self-injury), autism symptomatology, attention-deficit-hyperactivity-disorder characteristics, repetitive behaviour and mood/interest and pleasure were compared between each SNP genotype. No association was observed between behavioural characteristics and either 5-HTTLPR (serotonin) or MAOA (monoamine oxidase) genotypes. However, compared to the COMT (dopamine) AG and GG genotypes, the AA genotype was associated with greater interest and pleasure in the environment, and with reduced risk for property destruction, stereotyped behaviour and compulsive behaviour. The results suggest that common genetic variation in the COMT genotype affecting dopamine levels in the brain may contribute to the variability of challenging and repetitive behaviours and interest and pleasure in this population. This study identifies a role for additional genetic risk in understanding the neural and genetic mechanisms contributing to phenotypic variability in neurodevelopmental disorders, and highlights the merit of investigating SNPs that are selected a priori on a theoretical basis in rare populations.

Complete response to neoadjuvant chemoradiotherapy in rectal cancer is associated with RAS/AKT mutations and high tumour mutational burden.

BACKGROUND: Pathological complete response (pathCR) in rectal cancer is beneficial, as up to 75% of patients do not experience regrowth of the primary tumour, but it is poorly understood. We hypothesised that the changes seen in the pre-treatment biopsies of pathCR but not seen in residual tumour after chemoradiotherapy were the determinants of responsiveness. METHODS: Two groups of patients with either complete response (pathCR group, N = 24) or no response (poor response group, N = 24) were retrieved. Pre-treatment biopsies of cancers from these patients underwent high read depth amplicon sequencing for a targeted panel, exome sequencing, methylation profiling and immunohistochemistry for DNA repair pathway proteins. RESULTS: Twenty four patients who underwent pathCR and twenty-four who underwent poor response underwent molecular characterisation. Patients in the pathCR group had significantly higher tumour mutational burden and neoantigen load, frequent copy number alterations but fewer structural variants and enrichment for driver mutations in the PI3K/AKT/mTOR signalling pathway. There were no significant differences in tumour heterogeneity as measured by MATH score. Methylation analysis demonstrated enrichment for hypomethyation in the PI3K/AKT/mTOR signalling pathway. DISCUSSION: The phenomenon of pathCR in rectal cancer may be related to immunovisibility caused by a high tumour mutational burden phenotype. Potential therapy resistance mechanisms involve the PI3K/AKT/mTOR signalling pathway, but tumour heterogeneity does not seem to play a role in resistance.

Clinical-grade whole-genome sequencing and 3' transcriptome analysis of colorectal cancer patients.

BACKGROUND: Clinical-grade whole-genome sequencing (cWGS) has the potential to become the standard of care within the clinic because of its breadth of coverage and lack of bias towards certain regions of the genome. Colorectal cancer presents a difficult treatment paradigm, with over 40% of patients presenting at diagnosis with metastatic disease. We hypothesised that cWGS coupled with 3' transcriptome analysis would give new insights into colorectal cancer. METHODS: Patients underwent PCR-free whole-genome sequencing and alignment and variant calling using a standardised pipeline to output SNVs, indels, SVs and CNAs. Additional insights into the mutational signatures and tumour biology were gained by the use of 3' RNA-seq. RESULTS: Fifty-four patients were studied in total. Driver analysis identified the Wnt pathway gene APC as the only consistently mutated driver in colorectal cancer. Alterations in the PI3K/mTOR pathways were seen as previously observed in CRC. Multiple private CNAs, SVs and gene fusions were unique to individual tumours. Approximately 30% of patients had a tumour mutational burden of > 10 mutations/Mb of DNA, suggesting suitability for immunotherapy. CONCLUSIONS: Clinical whole-genome sequencing offers a potential avenue for the identification of private genomic variation that may confer sensitivity to targeted agents and offer patients new options for targeted therapies.

The transition from primary colorectal cancer to isolated peritoneal malignancy is associated with an increased tumour mutational burden.

Colorectal Peritoneal metastases (CPM) develop in 15% of colorectal cancers. Cytoreductive surgery and heated intraperitoneal chemotherapy (CRS & HIPEC) is the current standard of care in selected patients with limited resectable CPM. Despite selection using known prognostic factors survival is varied and morbidity and mortality are relatively high. There is a need to improve patient selection and a paucity of research concerning the biology of isolated CPM. We aimed to determine the biology associated with transition from primary CRC to CPM and of patients with CPM not responding to treatment with CRS & HIPEC, to identify those suitable for treatment with CRS & HIPEC and to identify targets for existing repurposed or novel treatment strategies. A cohort of patients with CPM treated with CRS & HIPEC was recruited and divided according to prognosis. Molecular profiling of the transcriptome (n = 25), epigenome (n = 24) and genome (n = 21) of CPM and matched primary CRC was performed. CPM were characterised by frequent Wnt/ ß catenin negative regulator mutations, TET2 mutations, mismatch repair mutations and high tumour mutational burden. Here we show the molecular features associated with CPM development and associated with not responding to CRS & HIPEC. Potential applications include improving patient selection for treatment with CRS & HIPEC and in future research into novel and personalised treatments targeting the molecular features identified here.

Rapid, highly accurate and cost-effective open-source simultaneous complete HLA typing and phasing of class I and II alleles using nanopore sequencing.

Accurate rapid genotyping of the genes within the HLA region presents many difficulties because of the complexity of this region. Here we present the results of our proof of concept nanopore-based long read polymerase chain reaction (PCR) solution for HLA genotyping. For 15 HLA anthropology-based samples and 13 NHS Blood and Transplant derived samples 40 ng of genomic DNA underwent long-range PCR for class I and II HLA alleles. Pooled PCR products were sequenced on the Oxford Nanopore MinIoON R9.4.1 flow cell. Sequenced reads had HLA genotype assigned with HLA-LA. Called genotypes were compared with reference derived from a combination of short-read next-generation sequencing, Sanger sequence and/or single-site polymorphism (SSP) typing. For concordance, accuracy was 100%, 98.4%, 97.5% and 95.1% for the first, second, third and fourth fields, respectively, to four field accuracy where it was available, otherwise three field in 28 samples for class I calls and 17 samples for class II calls. Phasing of maternal and paternal alleles, as well as phasing based identification of runs of homozygosity, was shown successfully. Time for assay run was 8 hours and the reconstruction of HLA typing data was 15 minutes. Assay cost was £55 ($80USD)/sample. We have developed a rapid and cost-effective long-range PCR and nanopore sequencing-based assay that can genotype the genes within HLA region to up to four field accuracy, identify runs of homozygosity in HLA, reconstruct maternal and paternal haplotypes and can be scaled from multi-sample runs to a single sample.

Treatment of COVID-19 with remdesivir in the absence of humoral immunity: a case report.

The response to the coronavirus disease 2019 (COVID-19) pandemic has been hampered by lack of an effective severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) antiviral therapy. Here we report the use of remdesivir in a patient with COVID-19 and the prototypic genetic antibody deficiency X-linked agammaglobulinaemia (XLA). Despite evidence of complement activation and a robust T cell response, the patient developed persistent SARS-CoV-2 pneumonitis, without progressing to multi-organ involvement. This unusual clinical course is consistent with a contribution of antibodies to both viral clearance and progression to severe disease. In the absence of these confounders, we take an experimental medicine approach to examine the in vivo utility of remdesivir. Over two independent courses of treatment, we observe a temporally correlated clinical and virological response, leading to clinical resolution and viral clearance, with no evidence of acquired drug resistance. We therefore provide evidence for the antiviral efficacy of remdesivir in vivo, and its potential benefit in selected patients.

Validation of epigenetic markers to identify colitis associated cancer: Results of module 1 of the ENDCAP-C study.

BACKGROUND: Chronic inflammation caused by ulcerative colitis (UC) causes a pro-neoplastic drive in the inflamed colon, leading to a markedly greater risk of invasive malignancy compared to the general population. Despite surveillance protocols, 50% of cases proceed to cancer before neoplasia is detected. The Enhanced Neoplasia Detection and Cancer Prevention in Chronic Colitis (ENDCaP-C) trial is an observational multi-centre test accuracy study to ascertain the role of molecular markers in improving the detection of dysplasia. We aimed to validate previously identified biomarkers of neoplasia in a retrospective cohort and create predictive models for later validation in a prospective cohort. METHODS: A retrospective analysis using bisulphite pyrosequencing of an 11 marker panel (SFRP1, SFRP2, SRP4, SRP5, WIF1, TUBB6, SOX7, APC1A, APC2, MINT1, RUNX3) in samples from 35 patients with cancer, 78 with dysplasia and 343 without neoplasia undergoing surveillance for UC associated neoplasia across 6 medical centres. Predictive models for UC associated cancer/dysplasia were created in the setting of neoplastic and non-neoplastic mucosa. FINDINGS: For neoplastic mucosa a five marker panel (SFRP2, SFRP4, WIF1, APC1A, APC2) was accurate in detecting pre-cancerous and invasive neoplasia (AUC = 0.83; 95% CI: 0.79, 0.88), and dysplasia (AUC = 0.88; (0.84, 0.91). For non-neoplastic mucosa a four marker panel (APC1A, SFRP4, SFRP5, SOX7) had modest accuracy (AUC = 0.68; 95% CI: 0.62,0.73) in predicting associated bowel neoplasia through the methylation signature of distant non-neoplastic colonic mucosa. INTERPRETATION: This multiplex methylation marker panel is accurate in the detection of ulcerative colitis associated dysplasia and neoplasia and is currently being validated in a prospective clinical trial. FUNDING: The ENDCAP-C study was funded by the National Institute for Health Research Efficacy and Mechanism Evaluation (EME) Programme (11/100/29).

Evaluation of the interaction between LRRK2 and PARK16 loci in determining risk of Parkinson's disease: analysis of a large multicenter study.

A recent study MacLeod et al. has shown that an interaction between variants at the LRRK2 and PARK16 loci influences risk of development of Parkinson's disease (PD). Our study examines the proposed interaction between LRRK2 and PARK16 variants in modifying PD risk using a large multicenter series of PD patients (7715) and controls (8261) from sites participating in the Genetic Epidemiology of Parkinson's Disease Consortium. Our data does not support a strong direct interaction between LRRK2 and PARK16 variants; however, given the role of retromer and lysosomal pathways in PD, further studies are warranted.

CCNF mutations in amyotrophic lateral sclerosis and frontotemporal dementia.

Amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD) are overlapping, fatal neurodegenerative disorders in which the molecular and pathogenic basis remains poorly understood. Ubiquitinated protein aggregates, of which TDP-43 is a major component, are a characteristic pathological feature of most ALS and FTD patients. Here we use genome-wide linkage analysis in a large ALS/FTD kindred to identify a novel disease locus on chromosome 16p13.3. Whole-exome sequencing identified a CCNF missense mutation at this locus. Interrogation of international cohorts identified additional novel CCNF variants in familial and sporadic ALS and FTD. Enrichment of rare protein-altering CCNF variants was evident in a large sporadic ALS replication cohort. CCNF encodes cyclin F, a component of an E3 ubiquitin-protein ligase complex (SCF(Cyclin F)). Expression of mutant CCNF in neuronal cells caused abnormal ubiquitination and accumulation of ubiquitinated proteins, including TDP-43 and a SCF(Cyclin F) substrate. This implicates common mechanisms, linked to protein homeostasis, underlying neuronal degeneration.