Whole genome sequencing enhances molecular diagnosis of primary ciliary dyskinesia.

BACKGROUND: Primary ciliary dyskinesia (PCD) is a genetic disorder affecting motile cilia. Most cases are inherited recessively, due to variants in >50 genes that result in abnormal or absent motile cilia. This leads to chronic upper and lower airway disease, subfertility, and laterality defects. Given overlapping clinical features and genetic heterogeneity, diagnosis can be difficult and often occurs late. Of those tested an estimated 30% of genetically screened PCD patients still lack a molecular diagnosis. A molecular diagnosis allows for appropriate clinical management including prediction of phenotypic features correlated to genotype. Here, we aimed to identify how readily a genetic diagnosis could be made using whole genome sequencing (WGS) to facilitate identification of pathogenic variants in known genes as well as novel PCD candidate genes. METHODS: WGS was used to screen for pathogenic variants in eight patients with PCD. RESULTS: 7/8 cases had homozygous or biallelic variants in DNAH5, DNAAF4 or DNAH11 classified as pathogenic or likely pathogenic. Three identified variants were deletions, ranging from 3 to 13 kb, for which WGS identified precise breakpoints, permitting confirmation by Sanger sequencing. WGS yielded identification of a de novo variant in a novel PCD gene TUBB4B. CONCLUSION: Here, WGS uplifted genetic diagnosis of PCD by identifying structural variants and novel modes of inheritance in new candidate genes. WGS could be an important component of the PCD diagnostic toolkit, increasing molecular diagnostic yield from current (70%) levels, and enhancing our understanding of fundamental biology of motile cilia and variants in the noncoding genome.

Rapid evolution of colistin resistance in a bioreactor model of infection of Klebsiella pneumoniae.

Colistin remains an important antibiotic for the therapeutic management of drug-resistant Klebsiella pneumoniae. Despite the numerous reports of colistin resistance in clinical strains, it remains unclear exactly when and how different mutational events arise resulting in reduced colistin susceptibility. Using a bioreactor model of infection, we modelled the emergence of colistin resistance in a susceptible isolate of K. pneumoniae. Genotypic, phenotypic and mathematical analyses of the antibiotic-challenged and un-challenged population indicates that after an initial decline, the population recovers within 24 h due to a small number of "founder cells" which have single point mutations mainly in the regulatory genes encoding crrB and pmrB that when mutated results in up to 100-fold reduction in colistin susceptibility. Our work underlines the rapid development of colistin resistance during treatment or exposure of susceptible K. pneumoniae infections having implications for the use of cationic antimicrobial peptides as a monotherapy.

Control of endothelial cell function and arteriogenesis by MEG3:EZH2 epigenetic regulation of integrin expression.

Epigenetic processes involving long non-coding RNAs regulate endothelial gene expression. However, the underlying regulatory mechanisms causing endothelial dysfunction remain to be elucidated. Enhancer of zeste homolog 2 (EZH2) is an important rheostat of histone H3K27 trimethylation (H3K27me3) that represses endothelial targets, but EZH2 RNA binding capacity and EZH2:RNA functional interactions have not been explored in post-ischemic angiogenesis. We used formaldehyde/UV-assisted crosslinking ligation and sequencing of hybrids and identified a new role for maternally expressed gene 3 (MEG3). MEG3 formed the predominant RNA:RNA hybrid structures in endothelial cells. Moreover, MEG3:EZH2 assists recruitment onto chromatin. By EZH2-chromatin immunoprecipitation, following MEG3 depletion, we demonstrated that MEG3 controls recruitment of EZH2/H3K27me3 onto integrin subunit alpha4 (ITGA4) promoter. Both MEG3 knockdown or EZH2 inhibition (A-395) promoted ITGA4 expression and improved endothelial cell migration and adhesion to fibronectin in vitro. The A-395 inhibitor re-directed MEG3-assisted chromatin remodeling, offering a direct therapeutic benefit by increasing endothelial function and resilience. This approach subsequently increased the expression of ITGA4 in arterioles following ischemic injury in mice, thus promoting arteriogenesis. Our findings show a context-specific role for MEG3 in guiding EZH2 to repress ITGA4. Novel therapeutic strategies could antagonize MEG3:EZH2 interaction for pre-clinical studies.

Ciliopathy patient variants reveal organelle-specific functions for TUBB4B in axonemal microtubules.

Tubulin, one of the most abundant cytoskeletal building blocks, has numerous isotypes in metazoans encoded by different conserved genes. Whether these distinct isotypes form cell type- and context-specific microtubule structures is poorly understood. Based on a cohort of 12 patients with primary ciliary dyskinesia as well as mouse mutants, we identified and characterized variants in the TUBB4B isotype that specifically perturbed centriole and cilium biogenesis. Distinct TUBB4B variants differentially affected microtubule dynamics and cilia formation in a dominant-negative manner. Structure-function studies revealed that different TUBB4B variants disrupted distinct tubulin interfaces, thereby enabling stratification of patients into three classes of ciliopathic diseases. These findings show that specific tubulin isotypes have distinct and nonredundant subcellular functions and establish a link between tubulinopathies and ciliopathies.

Genetic complexity of diagnostically unresolved Ehlers-Danlos syndrome.

BACKGROUND: The Ehlers-Danlos syndromes (EDS) are heritable disorders of connective tissue (HDCT), reclassified in the 2017 nosology into 13 subtypes. The genetic basis for hypermobile Ehlers-Danlos syndrome (hEDS) remains unknown. METHODS: Whole exome sequencing (WES) was undertaken on 174 EDS patients recruited from a national diagnostic service for complex EDS and a specialist clinic for hEDS. Patients had already undergone expert phenotyping, laboratory investigation and gene sequencing, but were without a genetic diagnosis. Filtered WES data were reviewed for genes underlying Mendelian disorders and loci reported in EDS linkage, transcriptome and genome-wide association studies (GWAS). A genetic burden analysis (Minor Allele Frequency (MAF) <0.05) incorporating 248 Avon Longitudinal Study of Parents and Children (ALSPAC) controls sequenced as part of the UK10K study was undertaken using TASER methodology. RESULTS: Heterozygous pathogenic (P) or likely pathogenic (LP) variants were identified in known EDS and Loeys-Dietz (LDS) genes. Multiple variants of uncertain significance where segregation and functional analysis may enable reclassification were found in genes associated with EDS, LDS, heritable thoracic aortic disease (HTAD), Mendelian disorders with EDS symptomatology and syndromes with EDS-like features. Genetic burden analysis revealed a number of novel loci, although none reached the threshold for genome-wide significance. Variants with biological plausibility were found in genes and pathways not currently associated with EDS or HTAD. CONCLUSIONS: We demonstrate the clinical utility of large panel-based sequencing and WES for patients with complex EDS in distinguishing rare EDS subtypes, LDS and related syndromes. Although many of the P and LP variants reported in this cohort would be identified with current panel testing, they were not at the time of this study, highlighting the use of extended panels and WES as a clinical tool for complex EDS. Our results are consistent with the complex genetic architecture of EDS and suggest a number of novel hEDS and HTAD candidate genes and pathways.

Short-read whole genome sequencing identifies causative variants in most individuals with previously unexplained aniridia.

BACKGROUND: Classic aniridia is a highly penetrant autosomal dominant disorder characterised by congenital absence of the iris, foveal hypoplasia, optic disc anomalies and progressive opacification of the cornea. >90% of cases of classic aniridia are caused by heterozygous, loss-of-function variants affecting the PAX6 locus. METHODS: Short-read whole genome sequencing was performed on 51 (39 affected) individuals from 37 different families who had screened negative for mutations in the PAX6 coding region. RESULTS: Likely causative mutations were identified in 22 out of 37 (59%) families. In 19 out of 22 families, the causative genomic changes have an interpretable deleterious impact on the PAX6 locus. Of these 19 families, 1 has a novel heterozygous PAX6 frameshift variant missed on previous screens, 4 have single nucleotide variants (SNVs) (one novel) affecting essential splice sites of PAX6 5' non-coding exons and 2 have deep intronic SNV (one novel) resulting in gain of a donor splice site. In 12 out of 19, the causative variants are large-scale structural variants; 5 have partial or whole gene deletions of PAX6, 3 have deletions encompassing critical PAX6 cis-regulatory elements, 2 have balanced inversions with disruptive breakpoints within the PAX6 locus and 2 have complex rearrangements disrupting PAX6. The remaining 3 of 22 families have deletions encompassing FOXC1 (a known cause of atypical aniridia). Seven of the causative variants occurred de novo and one cosegregated with familial aniridia. We were unable to establish inheritance status in the remaining probands. No plausibly causative SNVs were identified in PAX6 cis-regulatory elements. CONCLUSION: Whole genome sequencing proves to be an effective diagnostic test in most individuals with previously unexplained aniridia.

Whole exome sequencing of low grade serous ovarian carcinoma identifies genomic events associated with clinical outcome.

OBJECTIVES: Low-grade serous ovarian carcinoma (LGSOC) is a distinct, rare, ovarian cancer type characterised by younger patient age and intrinsic chemoresistance. Understanding the molecular landscape is crucial for optimising targeted therapy. METHODS: Genomic data from whole exome sequencing of tumour tissue was analysed in a LGSOC cohort with detailed clinical annotation. RESULTS: 63 cases were analysed and three subgroups identified based on single nucleotide variants: canonical MAPK mutant (cMAPKm: 52%, KRAS/BRAF/NRAS), MAPK-associated gene mutation (MAPK-assoc: 27%) and MAPK wild-type (MAPKwt: 21%). NOTCH pathway disruption occurred across all subgroups. Tumour mutational burden (TMB), mutational signatures and recurrent copy number (CN) changes varied across the cohort with co-occurrence of chromosome 1p loss and 1q gain (CN Chr1pq) a recurrent feature. Low TMB and CN Chr1pq were associated with inferior disease-specific survival (HR 6.43; p < 0.001 and HR 3.29, p = 0.011 respectively). Stepwise genomic classification in relation to outcome resulted in four groups (TMB low; CN Chr1pq; MAPKwt/MAPKassoc; cMAPKm). 5 year disease-specific survival was 46%, 55%, 79% and 100% respectively for these groups. The two most favourable genomic subgroups were enriched for the SBS10b mutational signature, particularly the cMAPKm subgroup. CONCLUSIONS: LGSOC comprises multiple genomic subgroups with distinct clinical and molecular features. Chr1pq CN arm disruption and TMB represent promising methods to identify individuals with poorer prognosis. Further investigation of the molecular basis for these observations is required. MAPKwt cases represent around a fifth of patients. NOTCH inhibitors represent a candidate therapeutic strategy worthy of exploration across these cases.

Distinct histopathological features are associated with molecular subtypes and outcome in low grade serous ovarian carcinoma.

Low grade serous ovarian carcinoma (LGSOC) demonstrates unique clinical and molecular features compared to other ovarian cancer types. The relationship between common histological features of LGSOC and molecular events, such as hormone receptor expression patterns and MAPK gene mutation status, remains poorly understood. Recent data suggest some of these molecular features may be biomarkers of response to recently introduced biologically-targeted therapies, namely endocrine therapy and MEK inhibitors. We utilize a cohort of 63 pathologically-confirmed LGSOC cases with whole exome sequencing and hormone receptor expression data to investigate these relationships. LGSOC cases demonstrated uniformly high oestrogen receptor (ER) expression, but variable progesterone receptor (PR) expression intensity. 60% and 37% of cases demonstrated micropapillary and macropapillary patterns of stromal invasion, respectively. 63% of cases demonstrated desmoplasia, which was significantly associated with advanced disease stage and visible residual disease after cytoreductive surgery. MAPK-mutant cases (KRAS, BRAF, NRAS) more frequently demonstrated macropapillary stromal invasion, while Chr1p loss was associated with desmoplasia and low PR expression. Presence of micropapillary stromal invasion and low PR expression were associated with significantly poorer survival after accounting for stage and residual disease status. Together, these data identify novel relationships between histopathological features and molecularly-defined subgroups in LGSOC.

Toll-like receptor 2 orchestrates a tumor suppressor response in non-small cell lung cancer.

Targeting early-stage lung cancer is vital to improve survival. However, the mechanisms and components of the early tumor suppressor response in lung cancer are not well understood. In this report, we study the role of Toll-like receptor 2 (TLR2), a regulator of oncogene-induced senescence, which is a key tumor suppressor response in premalignancy. Using human lung cancer samples and genetically engineered mouse models, we show that TLR2 is active early in lung tumorigenesis, where it correlates with improved survival and clinical regression. Mechanistically, TLR2 impairs early lung cancer progression via activation of cell intrinsic cell cycle arrest pathways and the proinflammatory senescence-associated secretory phenotype (SASP). The SASP regulates non-cell autonomous anti-tumor responses, such as immune surveillance of premalignant cells, and we observe impaired myeloid cell recruitment to lung tumors after Tlr2 loss. Last, we show that administration of a TLR2 agonist reduces lung tumor growth, highlighting TLR2 as a possible therapeutic target.

Robust Genetic Analysis of the X-Linked Anophthalmic (Ie) Mouse.

Anophthalmia (missing eye) describes a failure of early embryonic ocular development. Mutations in a relatively small set of genes account for 75% of bilateral anophthalmia cases, yet 25% of families currently are left without a molecular diagnosis. Here, we report our experimental work that aimed to uncover the developmental and genetic basis of the anophthalmia characterising the X-linked Ie (eye-ear reduction) X-ray-induced allele in mouse that was first identified in 1947. Histological analysis of the embryonic phenotype showed failure of normal eye development after the optic vesicle stage with particularly severe malformation of the ventral retina. Linkage analysis mapped this mutation to a ~6 Mb region on the X chromosome. Short- and long-read whole-genome sequencing (WGS) of affected and unaffected male littermates confirmed the Ie linkage but identified no plausible causative variants or structural rearrangements. These analyses did reduce the critical candidate interval and revealed evidence of multiple variants within the ancestral DNA, although none were found that altered coding sequences or that were unique to Ie. To investigate early embryonic events at a genetic level, we then generated mouse ES cells derived from male Ie embryos and wild type littermates. RNA-seq and accessible chromatin sequencing (ATAC-seq) data generated from cultured optic vesicle organoids did not reveal any large differences in gene expression or accessibility of putative cis-regulatory elements between Ie and wild type. However, an unbiased TF-footprinting analysis of accessible chromatin regions did provide evidence of a genome-wide reduction in binding of transcription factors associated with ventral eye development in Ie, and evidence of an increase in binding of the Zic-family of transcription factors, including Zic3, which is located within the Ie-refined critical interval. We conclude that the refined Ie critical region at chrX: 56,145,000-58,385,000 contains multiple genetic variants that may be linked to altered cis regulation but does not contain a convincing causative mutation. Changes in the binding of key transcription factors to chromatin causing altered gene expression during development, possibly through a subtle mis-regulation of Zic3, presents a plausible cause for the anophthalmia phenotype observed in Ie, but further work is required to determine the precise causative allele and its genetic mechanism.

Gorham-Stout case report: a multi-omic analysis reveals recurrent fusions as new potential drivers of the disease.

BACKGROUND: Gorham-Stout disease is a rare condition characterized by vascular proliferation and the massive destruction of bone tissue. With less than 400 cases in the literature of Gorham-Stout syndrome, we performed a unique study combining whole-genome sequencing and RNA-Seq to probe the genomic features and differentially expressed pathways of a presented case, revealing new possible drivers and biomarkers of the disease. CASE PRESENTATION: We present a case report of a white 45-year-old female patient with marked bone loss of the left humerus associated with vascular proliferation, diagnosed with Gorham-Stout disease. The analysis of whole-genome sequencing showed a dominance of large structural DNA rearrangements. Particularly, rearrangements in chromosomes seven, twelve, and twenty could contribute to the development of the disease, especially a gene fusion involving ATG101 that could affect macroautophagy. The study of RNA-sequencing data from the patient uncovered the PI3K/AKT/mTOR pathway as the most affected signaling cascade in the Gorham-Stout lesional tissue. Furthermore, M2 macrophage infiltration was detected using immunohistochemical staining and confirmed by deconvolution of the RNA-seq expression data. CONCLUSIONS: The way that DNA and RNA aberrations lead to Gorham-Stout disease is poorly understood due to the limited number of studies focusing on this rare disease. Our study provides the first glimpse into this facet of the disease, exposing new possible therapeutic targets and facilitating the clinicopathological diagnosis of Gorham-Stout disease.

Multiomic Characterization of High-Grade Serous Ovarian Carcinoma Enables High-Resolution Patient Stratification.

PURPOSE: High-grade serous ovarian carcinoma (HGSOC) is the most common ovarian cancer type; most patients experience disease recurrence that accumulates chemoresistance, leading to treatment failure. Genomic and transcriptomic features have been associated with differential outcome and treatment response. However, the relationship between events at the gene sequence, copy number, and gene-expression levels remains poorly defined. EXPERIMENTAL DESIGN: We perform multiomic characterization of a large HGSOC cohort (n = 362) with detailed clinical annotation to interrogate the relationship between patient subgroups defined by specific molecular events. RESULTS: BRCA2-mutant (BRCA2m) and EMSY-overexpressing cases demonstrated prolonged survival [multivariable hazard ratios (HR) 0.40 and 0.51] and significantly higher first- and second-line chemotherapy response rate. CCNE1-gained (CCNE1g) cases demonstrated underrepresentation of FIGO stage IV cases, with shorter survival but no significant difference in treatment response. We demonstrate marked overlap between the TCGA- and Tothill-derived subtypes. IMR/C2 cases displayed higher BRCA1/2m frequency (25.5%, 32.5%) and significantly greater immune cell infiltration, whereas PRO/C5 cases had the highest CCNE1g rate (23.9%, 22.2%) and were uniformly low in immune cell infiltration. The survival benefit for cases with aberrations in homologous recombination repair (HRR) genes was apparent across all transcriptomic subtypes (HR range, 0.48-0.68). There was significant co-occurrence of RB loss and HRR gene aberrations; RB loss was further associated with favorable survival within HRR-aberrant cases (multivariable HR, 0.50). CONCLUSIONS: These data paint a high-resolution picture of the molecular landscape in HGSOC, better defining patients who may benefit most from specific molecular therapeutics and highlighting those for whom novel treatment strategies are needed to improve outcomes.

Ancient DNA at the edge of the world: Continental immigration and the persistence of Neolithic male lineages in Bronze Age Orkney.

Orkney was a major cultural center during the Neolithic, 3800 to 2500 BC. Farming flourished, permanent stone settlements and chambered tombs were constructed, and long-range contacts were sustained. From ∼3200 BC, the number, density, and extravagance of settlements increased, and new ceremonial monuments and ceramic styles, possibly originating in Orkney, spread across Britain and Ireland. By ∼2800 BC, this phenomenon was waning, although Neolithic traditions persisted to at least 2500 BC. Unlike elsewhere in Britain, there is little material evidence to suggest a Beaker presence, suggesting that Orkney may have developed along an insular trajectory during the second millennium BC. We tested this by comparing new genomic evidence from 22 Bronze Age and 3 Iron Age burials in northwest Orkney with Neolithic burials from across the archipelago. We identified signals of inward migration on a scale unsuspected from the archaeological record: As elsewhere in Bronze Age Britain, much of the population displayed significant genome-wide ancestry deriving ultimately from the Pontic-Caspian Steppe. However, uniquely in northern and central Europe, most of the male lineages were inherited from the local Neolithic. This suggests that some male descendants of Neolithic Orkney may have remained distinct well into the Bronze Age, although there are signs that this had dwindled by the Iron Age. Furthermore, although the majority of mitochondrial DNA lineages evidently arrived afresh with the Bronze Age, we also find evidence for continuity in the female line of descent from Mesolithic Britain into the Bronze Age and even to the present day.

In Vivo Modeling of Patient Genetic Heterogeneity Identifies New Ways to Target Cholangiocarcinoma.

Intrahepatic cholangiocarcinoma (ICC) is an aggressive malignancy of the bile ducts within the liver characterized by high levels of genetic heterogeneity. In the context of such genetic variability, determining which oncogenic mutations drive ICC growth has been difficult, and developing modes of patient stratification and targeted therapies remains challenging. Here we model the interactions between rare mutations with more common driver genes and combine in silico analysis of patient data with highly multiplexed in vivo CRISPR-spCas9 screens to perform a functional in vivo study into the role genetic heterogeneity plays in driving ICC. Novel tumor suppressors were uncovered, which, when lost, cooperate with the RAS oncoprotein to drive ICC growth. Focusing on a set of driver mutations that interact with KRAS to initiate aggressive, sarcomatoid-type ICC revealed that tumor growth relies on Wnt and PI3K signaling. Pharmacologic coinhibition of Wnt and PI3K in vivo impeded ICC growth regardless of mutational profile. Therefore, Wnt and PI3K activity should be considered as a signature by which patients can be stratified for treatment independent of tumor genotype, and inhibitors of these pathways should be levied to treat ICC. SIGNIFICANCE: This work shows that, despite significant genetic heterogeneity, intrahepatic cholangiocarcinoma relies on a limited number of signaling pathways to grow, suggesting common therapeutic vulnerabilities across patients.

A prenylated dsRNA sensor protects against severe COVID-19.

Inherited genetic factors can influence the severity of COVID-19, but the molecular explanation underpinning a genetic association is often unclear. Intracellular antiviral defenses can inhibit the replication of viruses and reduce disease severity. To better understand the antiviral defenses relevant to COVID-19, we used interferon-stimulated gene (ISG) expression screening to reveal that 2'-5'-oligoadenylate synthetase 1 (OAS1), through ribonuclease L, potently inhibits severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). We show that a common splice-acceptor single-nucleotide polymorphism (Rs10774671) governs whether patients express prenylated OAS1 isoforms that are membrane-associated and sense-specific regions of SARS-CoV-2 RNAs or if they only express cytosolic, nonprenylated OAS1 that does not efficiently detect SARS-CoV-2. In hospitalized patients, expression of prenylated OAS1 was associated with protection from severe COVID-19, suggesting that this antiviral defense is a major component of a protective antiviral response.

Integrated molecular characterisation of endometrioid ovarian carcinoma identifies opportunities for stratification.

Endometrioid ovarian carcinoma (EnOC) is an under-investigated ovarian cancer type. Recent studies have described disease subtypes defined by genomics and hormone receptor expression patterns; here, we determine the relationship between these subtyping layers to define the molecular landscape of EnOC with high granularity and identify therapeutic vulnerabilities in high-risk cases. Whole exome sequencing data were integrated with progesterone and oestrogen receptor (PR and ER) expression-defined subtypes in 90 EnOC cases following robust pathological assessment, revealing dominant clinical and molecular features in the resulting integrated subtypes. We demonstrate significant correlation between subtyping approaches: PR-high (PR + /ER + , PR + /ER-) cases were predominantly CTNNB1-mutant (73.2% vs 18.4%, P < 0.001), while PR-low (PR-/ER + , PR-/ER-) cases displayed higher TP53 mutation frequency (38.8% vs 7.3%, P = 0.001), greater genomic complexity (P = 0.007) and more frequent copy number alterations (P = 0.001). PR-high EnOC patients experience favourable disease-specific survival independent of clinicopathological and genomic features (HR = 0.16, 95% CI 0.04-0.71). TP53 mutation further delineates the outcome of patients with PR-low tumours (HR = 2.56, 95% CI 1.14-5.75). A simple, routinely applicable, classification algorithm utilising immunohistochemistry for PR and p53 recapitulated these subtypes and their survival profiles. The genomic profile of high-risk EnOC subtypes suggests that inhibitors of the MAPK and PI3K-AKT pathways, alongside PARP inhibitors, represent promising candidate agents for improving patient survival. Patients with PR-low TP53-mutant EnOC have the greatest unmet clinical need, while PR-high tumours-which are typically CTNNB1-mutant and TP53 wild-type-experience excellent survival and may represent candidates for trials investigating de-escalation of adjuvant chemotherapy to agents such as endocrine therapy.

Structural Variants at the BRCA1/2 Loci are a Common Source of Homologous Repair Deficiency in High-grade Serous Ovarian Carcinoma.

PURPOSE: The abundance and effects of structural variation at BRCA1/2 in tumors are not well understood. In particular, the impact of these events on homologous recombination repair deficiency (HRD) has yet to be demonstrated. EXPERIMENTAL DESIGN: Exploiting a large collection of whole-genome sequencing data from high-grade serous ovarian carcinoma (N = 205) together with matched RNA sequencing for the majority of tumors (N = 150), we have comprehensively characterized mutation and expression at BRCA1/2. RESULTS: In addition to the known spectrum of short somatic mutations (SSM), we discovered that multi-megabase structural variants (SV) were a frequent, unappreciated source of BRCA1/2 disruption in these tumors, and we found a genome-wide enrichment for large deletions at the BRCA1/2 loci across the cohort. These SVs independently affected a substantial proportion of patients (16%) in addition to those affected by SSMs (24%), conferring HRD and impacting patient survival. We also detail compound deficiencies involving SSMs and SVs at both loci, demonstrating that the strongest risk of HRD emerges from combined SVs at both BRCA1 and BRCA2 in the absence of SSMs. Furthermore, these SVs are abundant and disruptive in other cancer types. CONCLUSIONS: These results extend our understanding of the mutational landscape underlying HRD, increase the number of patients predicted to benefit from therapies exploiting HRD, and suggest there is currently untapped potential in SV detection for patient stratification.

Genetic mechanisms of critical illness in COVID-19.

Host-mediated lung inflammation is present1, and drives mortality2, in the critical illness caused by coronavirus disease 2019 (COVID-19). Host genetic variants associated with critical illness may identify mechanistic targets for therapeutic development3. Here we report the results of the GenOMICC (Genetics Of Mortality In Critical Care) genome-wide association study in 2,244 critically ill patients with COVID-19 from 208 UK intensive care units. We have identified and replicated the following new genome-wide significant associations: on chromosome 12q24.13 (rs10735079, P = 1.65 × 10-8) in a gene cluster that encodes antiviral restriction enzyme activators (OAS1, OAS2 and OAS3); on chromosome 19p13.2 (rs74956615, P = 2.3 × 10-8) near the gene that encodes tyrosine kinase 2 (TYK2); on chromosome 19p13.3 (rs2109069, P = 3.98 ×  10-12) within the gene that encodes dipeptidyl peptidase 9 (DPP9); and on chromosome 21q22.1 (rs2236757, P = 4.99 × 10-8) in the interferon receptor gene IFNAR2. We identified potential targets for repurposing of licensed medications: using Mendelian randomization, we found evidence that low expression of IFNAR2, or high expression of TYK2, are associated with life-threatening disease; and transcriptome-wide association in lung tissue revealed that high expression of the monocyte-macrophage chemotactic receptor CCR2 is associated with severe COVID-19. Our results identify robust genetic signals relating to key host antiviral defence mechanisms and mediators of inflammatory organ damage in COVID-19. Both mechanisms may be amenable to targeted treatment with existing drugs. However, large-scale randomized clinical trials will be essential before any change to clinical practice.

Molecular stratification of endometrioid ovarian carcinoma predicts clinical outcome.

Endometrioid ovarian carcinoma (EnOC) demonstrates substantial clinical and molecular heterogeneity. Here, we report whole exome sequencing of 112 EnOC cases following rigorous pathological assessment. We detect a high frequency of mutation in CTNNB1 (43%), PIK3CA (43%), ARID1A (36%), PTEN (29%), KRAS (26%), TP53 (26%) and SOX8 (19%), a recurrently-mutated gene previously unreported in EnOC. POLE and mismatch repair protein-encoding genes were mutated at lower frequency (6%, 18%) with significant co-occurrence. A molecular taxonomy is constructed, identifying clinically distinct EnOC subtypes: cases with TP53 mutation demonstrate greater genomic complexity, are commonly FIGO stage III/IV at diagnosis (48%), are frequently incompletely debulked (44%) and demonstrate inferior survival; conversely, cases with CTNNB1 mutation, which is mutually exclusive with TP53 mutation, demonstrate low genomic complexity and excellent clinical outcome, and are predominantly stage I/II at diagnosis (89%) and completely resected (87%). Moreover, we identify the WNT, MAPK/RAS and PI3K pathways as good candidate targets for molecular therapeutics in EnOC.

Increased ultra-rare variant load in an isolated Scottish population impacts exonic and regulatory regions.

Human population isolates provide a snapshot of the impact of historical demographic processes on population genetics. Such data facilitate studies of the functional impact of rare sequence variants on biomedical phenotypes, as strong genetic drift can result in higher frequencies of variants that are otherwise rare. We present the first whole genome sequencing (WGS) study of the VIKING cohort, a representative collection of samples from the isolated Shetland population in northern Scotland, and explore how its genetic characteristics compare to a mainland Scottish population. Our analyses reveal the strong contributions played by the founder effect and genetic drift in shaping genomic variation in the VIKING cohort. About one tenth of all high-quality variants discovered are unique to the VIKING cohort or are seen at frequencies at least ten fold higher than in more cosmopolitan control populations. Multiple lines of evidence also suggest relaxation of purifying selection during the evolutionary history of the Shetland isolate. We demonstrate enrichment of ultra-rare VIKING variants in exonic regions and for the first time we also show that ultra-rare variants are enriched within regulatory regions, particularly promoters, suggesting that gene expression patterns may diverge relatively rapidly in human isolates.