Genomic structures and regulation patterns at HPV integration sites in cervical cancer.

Human papillomavirus (HPV) integration has been implicated in transforming HPV infection into cancer, but its genomic consequences have been difficult to study using short-read technologies. To resolve the dysregulation associated with HPV integration, we performed long-read sequencing on 63 cervical cancer genomes. We identified six categories of integration events based on HPV-human genomic structures. Of all HPV integrants, defined as two HPV-human breakpoints bridged by an HPV sequence, 24% contained variable copies of HPV between the breakpoints, a phenomenon we termed heterologous integration. Analysis of DNA methylation within and in proximity to the HPV genome at individual integration events revealed relationships between methylation status of the integrant and its orientation and structure. Dysregulation of the human epigenome and neighboring gene expression in cis with the HPV-integrated allele was observed over megabase-ranges of the genome. By elucidating the structural, epigenetic, and allele-specific impacts of HPV integration, we provide insight into the role of integrated HPV in cervical cancer.

The genome sequence of the Loggerhead sea turtle, Caretta caretta Linnaeus 1758.

We present a genome assembly of Caretta caretta (the Loggerhead sea turtle; Chordata, Testudines, Cheloniidae), generated from genomic data from two unrelated females. The genome sequence is 2.13 gigabases in size. The assembly has a busco completion score of 96.1% and N50 of 130.95 Mb. The majority of the assembly is scaffolded into 28 chromosomal representations with a remaining 2% of the assembly being excluded from these.

Whole-genome and transcriptome analysis of advanced adrenocortical cancer highlights multiple alterations affecting epigenome and DNA repair pathways.

Adrenocortical cancer (ACC) is a rare cancer of the adrenal gland. Several driver mutations have been identified in both primary and metastatic ACCs, but the therapeutic options are still limited. We performed whole-genome and transcriptome sequencing on seven patients with metastatic ACC. Integrative analysis of mutations, RNA expression changes, mutation signature, and homologous recombination deficiency (HRD) analysis was performed. Mutations affecting CTNNB1 and TP53 and frequent loss of heterozygosity (LOH) events were observed in our cohort. Alterations affecting genes involved in cell cycle (RB1, CDKN2A, CDKN2B), DNA repair pathways (MUTYH, BRCA2, ATM, RAD52, MLH1, MSH6), and telomere maintenance (TERF2 and TERT) consisting of somatic and germline mutations, structural variants, and expression outliers were also observed. HRDetect, which aggregates six HRD-associated mutation signatures, identified a subset of cases as HRD. Genomic alterations affecting genes involved in epigenetic regulation were also identified, including structural variants (SWI/SNF genes and histone methyltransferases), and copy gains and concurrent high expression of KDM5A, which may contribute to epigenomic deregulation. Findings from this study highlight HRD and epigenomic pathways as potential therapeutic targets and suggest a subgroup of patients may benefit from a diverse array of molecularly targeted therapies in ACC, a rare disease in urgent need of therapeutic strategies.

Long-read genome sequencing resolves a complex 13q structural variant associated with syndromic anophthalmia.

Microphthalmia, anophthalmia, and coloboma (MAC) are a heterogeneous spectrum of anomalous eye development and degeneration with genetic and environmental etiologies. Structural and copy number variants of chromosome 13 have been implicated in MAC; however, the specific loci involved in disease pathogenesis have not been well-defined. Herein we report a newborn with syndromic degenerative anophthalmia and a complex de novo rearrangement of chromosome 13q. Long-read genome sequencing improved the resolution and clinical interpretation of a duplication-triplication/inversion-duplication (DUP-TRP/INV-DUP) and terminal deletion. Sequence features at the breakpoint junctions suggested microhomology-mediated break-induced replication (MMBIR) of the maternal chromosome as the origin. Comparing this rearrangement to previously reported copy number alterations in 13q, we refine a putative dosage-sensitive critical region for MAC that might provide new insights into its molecular etiology.

The genome of the forest insect pest Pissodes strobi reveals genome expansion and evidence of a Wolbachia endosymbiont.

The highly diverse insect family of true weevils, Curculionidae, includes many agricultural and forest pests. Pissodes strobi, commonly known as the spruce weevil or white pine weevil, is a major pest of spruce and pine forests in North America. Pissodes strobi larvae feed on the apical shoots of young trees, causing stunted growth and can destroy regenerating spruce or pine forests. Here, we describe the nuclear and mitochondrial Pissodes strobi genomes and their annotations, as well as the genome of an apparent Wolbachia endosymbiont. We report a substantial expansion of the weevil nuclear genome, relative to other Curculionidae species, possibly driven by an abundance of class II DNA transposons. The endosymbiont observed belongs to a group (supergroup A) of Wolbachia species that generally form parasitic relationships with their arthropod host.

Releasates of riboflavin/UV-treated platelets: Microvesicles suppress cytokine-mediated endothelial cell migration/proliferation.

BACKGROUND: Accelerated development of the platelet (PLT) storage lesion upon pathogen inactivation (PI) is associated with the release of proteins from granules and platelet microvesicles (PMVs). Whether PI treatments alter the interaction between PLT factors and the vessel endothelium is of interest in understanding the risk profile of these technologies. STUDY DESIGN AND METHODS: In a pool-and-split study, one platelet concentrate (PC) was treated with riboflavin/UV (RF/UV) light, while the other one was kept as an untreated control. Releasates and PMV-depleted releasates were prepared by differential centrifugation steps on days 0, 1, 5, and 7 of storage. Cytokine/chemokine release following PI treatment was analyzed by an antibody array, and results were verified by the enzyme-linked immunosorbent assay. PMVs were enumerated by CD41 labeling and flow cytometry. Wound scratch assays were performed using cultured Ea.hy926 cells exposed to the differently prepared releasates. Effects of releasates on the phosphorylation levels of kinases ERK and p38 expressed by endothelial cells were analyzed by immunoblot. RESULTS: Cytokine/chemokine assays identified a 2-fold increase in epidermal growth factor released from PCs treated with RF/UV light compared with control. PMV count increased ~100-fold following PI treatment. Unmodified releasates and PMV-depleted releasates displayed different contributions to the kinetics of endothelial cell wound closure. This observation was associated with an increased ERK versus unaltered p38 activation in the endothelial cells. CONCLUSION: This study identified an inhibitory impact of PMVs on endothelial cell migration/proliferation upon stimulation by released cytokines and PMVs from PLTs treated with RF/UV light for endothelial cell wound closure.

Analysis of Ugandan cervical carcinomas identifies human papillomavirus clade-specific epigenome and transcriptome landscapes.

Cervical cancer is the most common cancer affecting sub-Saharan African women and is prevalent among HIV-positive (HIV+) individuals. No comprehensive profiling of cancer genomes, transcriptomes or epigenomes has been performed in this population thus far. We characterized 118 tumors from Ugandan patients, of whom 72 were HIV+, and performed extended mutation analysis on an additional 89 tumors. We detected human papillomavirus (HPV)-clade-specific differences in tumor DNA methylation, promoter- and enhancer-associated histone marks, gene expression and pathway dysregulation. Changes in histone modification at HPV integration events were correlated with upregulation of nearby genes and endogenous retroviruses.

Endogenous Retrovirus Transcript Levels Are Associated with Immunogenic Signatures in Multiple Metastatic Cancer Types.

Next-generation sequencing of solid tumors has revealed variable signatures of immunogenicity across tumors, but underlying molecular characteristics driving such variation are not fully understood. Although expression of endogenous retrovirus (ERV)-containing transcripts can provide a source of tumor-specific neoantigen in some cancer models, associations between ERV levels and immunogenicity across different types of metastatic cancer are not well established. We performed bioinformatics analysis of genomic, transcriptomic, and clinical data across an integrated cohort of 199 patients with metastatic breast, colorectal, and pancreatic ductal adenocarcinoma tumors. Within each cancer type, we identified a subgroup of viral mimicry tumors in which increased ERV levels were coupled with transcriptional signatures of autonomous antiviral response and immunogenicity. In addition, viral mimicry colorectal and pancreatic tumors showed increased expression of DNA demethylation gene TET2 Taken together, these data demonstrate the existence of an ERV-associated viral mimicry phenotype across three distinct metastatic cancer types, while indicating links between ERV abundance, epigenetic dysregulation, and immunogenicity.

Pan-cancer analysis of advanced patient tumors reveals interactions between therapy and genomic landscapes.

Advanced and metastatic tumors with complex treatment histories drive cancer mortality. Here we describe the POG570 cohort, a comprehensive whole-genome, transcriptome and clinical dataset, amenable for exploration of the impacts of therapies on genomic landscapes. Previous exposure to DNA-damaging chemotherapies and mutations affecting DNA repair genes, including POLQ and genes encoding Polζ, were associated with genome-wide, therapy-induced mutagenesis. Exposure to platinum therapies coincided with signatures SBS31 and DSB5 and, when combined with DNA synthesis inhibitors, signature SBS17b. Alterations in ESR1, EGFR, CTNNB1, FGFR1, VEGFA and DPYD were consistent with drug resistance and sensitivity. Recurrent noncoding events were found in regulatory region hotspots of genes including TERT, PLEKHS1, AP2A1 and ADGRG6. Mutation burden and immune signatures corresponded with overall survival and response to immunotherapy. Our data offer a rich resource for investigation of advanced cancers and interpretation of whole-genome and transcriptome sequencing in the context of a cancer clinic.

Altered Gene Expression along the Glycolysis-Cholesterol Synthesis Axis Is Associated with Outcome in Pancreatic Cancer.

PURPOSE: Identification of clinically actionable molecular subtypes of pancreatic ductal adenocarcinoma (PDAC) is key to improving patient outcome. Intertumoral metabolic heterogeneity contributes to cancer survival and the balance between distinct metabolic pathways may influence PDAC outcome. We hypothesized that PDAC can be stratified into prognostic metabolic subgroups based on alterations in the expression of genes involved in glycolysis and cholesterol synthesis. EXPERIMENTAL DESIGN: We performed bioinformatics analysis of genomic, transcriptomic, and clinical data in an integrated cohort of 325 resectable and nonresectable PDAC. The resectable datasets included retrospective The Cancer Genome Atlas (TCGA) and the International Cancer Genome Consortium (ICGC) cohorts. The nonresectable PDAC cohort studies included prospective COMPASS, PanGen, and BC Cancer Personalized OncoGenomics program (POG). RESULTS: On the basis of the median normalized expression of glycolytic and cholesterogenic genes, four subgroups were identified: quiescent, glycolytic, cholesterogenic, and mixed. Glycolytic tumors were associated with the shortest median survival in resectable (log-rank test P = 0.018) and metastatic settings (log-rank test P = 0.027). Patients with cholesterogenic tumors had the longest median survival. KRAS and MYC-amplified tumors had higher expression of glycolytic genes than tumors with normal or lost copies of the oncogenes (Wilcoxon rank sum test P = 0.015). Glycolytic tumors had the lowest expression of mitochondrial pyruvate carriers MPC1 and MPC2. Glycolytic and cholesterogenic gene expression correlated with the expression of prognostic PDAC subtype classifier genes. CONCLUSIONS: Metabolic classification specific to glycolytic and cholesterogenic pathways provides novel biological insight into previously established PDAC subtypes and may help develop personalized therapies targeting unique tumor metabolic profiles.See related commentary by Mehla and Singh, p. 6.

Text-mining clinically relevant cancer biomarkers for curation into the CIViC database.

BACKGROUND: Precision oncology involves analysis of individual cancer samples to understand the genes and pathways involved in the development and progression of a cancer. To improve patient care, knowledge of diagnostic, prognostic, predisposing, and drug response markers is essential. Several knowledgebases have been created by different groups to collate evidence for these associations. These include the open-access Clinical Interpretation of Variants in Cancer (CIViC) knowledgebase. These databases rely on time-consuming manual curation from skilled experts who read and interpret the relevant biomedical literature. METHODS: To aid in this curation and provide the greatest coverage for these databases, particularly CIViC, we propose the use of text mining approaches to extract these clinically relevant biomarkers from all available published literature. To this end, a group of cancer genomics experts annotated sentences that discussed biomarkers with their clinical associations and achieved good inter-annotator agreement. We then used a supervised learning approach to construct the CIViCmine knowledgebase. RESULTS: We extracted 121,589 relevant sentences from PubMed abstracts and PubMed Central Open Access full-text papers. CIViCmine contains over 87,412 biomarkers associated with 8035 genes, 337 drugs, and 572 cancer types, representing 25,818 abstracts and 39,795 full-text publications. CONCLUSIONS: Through integration with CIVIC, we provide a prioritized list of curatable clinically relevant cancer biomarkers as well as a resource that is valuable to other knowledgebases and precision cancer analysts in general. All data is publically available and distributed with a Creative Commons Zero license. The CIViCmine knowledgebase is available at http://bionlp.bcgsc.ca/civicmine/.

The Genome of the Steller Sea Lion (Eumetopias jubatus).

The Steller sea lion is the largest member of the Otariidae family and is found in the coastal waters of the northern Pacific Rim. Here, we present the Steller sea lion genome, determined through DNA sequencing approaches that utilized microfluidic partitioning library construction, as well as nanopore technologies. These methods constructed a highly contiguous assembly with a scaffold N50 length of over 14 megabases, a contig N50 length of over 242 kilobases and a total length of 2.404 gigabases. As a measure of completeness, 95.1% of 4104 highly conserved mammalian genes were found to be complete within the assembly. Further annotation identified 19,668 protein coding genes. The assembled genome sequence and underlying sequence data can be found at the National Center for Biotechnology Information (NCBI) under the BioProject accession number PRJNA475770.

Phagocytosis of Aspergillus fumigatus by Human Bronchial Epithelial Cells Is Mediated by the Arp2/3 Complex and WIPF2.

Aspergillus fumigatus is an opportunistic fungal pathogen capable of causing severe infection in humans. One of the limitations in our understanding of how A. fumigatus causes infection concerns the initial stages of infection, notably the initial interaction between inhaled spores or conidia and the human airway. Using publicly-available datasets, we identified the Arp2/3 complex and the WAS-Interacting Protein Family Member 2 WIPF2 as being potentially responsible for internalization of conidia by airway epithelial cells. Using a cell culture model, we demonstrate that RNAi-mediated knockdown of WIPF2 significantly reduces internalization of conidia into airway epithelial cells. Furthermore, we demonstrate that inhibition of Arp2/3 by a small molecule inhibitor causes similar effects. Using super-resolution fluorescence microscopy, we demonstrate that WIPF2 is transiently localized to the site of bound conidia. Overall, we demonstrate the active role of the Arp2/3 complex and WIPF2 in mediating the internalization of A. fumigatus conidia into human airway epithelial cells.

Genome-wide discovery of somatic coding and noncoding mutations in pediatric endemic and sporadic Burkitt lymphoma.

Although generally curable with intensive chemotherapy in resource-rich settings, Burkitt lymphoma (BL) remains a deadly disease in older patients and in sub-Saharan Africa. Epstein-Barr virus (EBV) positivity is a feature in more than 90% of cases in malaria-endemic regions, and up to 30% elsewhere. However, the molecular features of BL have not been comprehensively evaluated when taking into account tumor EBV status or geographic origin. Through an integrative analysis of whole-genome and transcriptome data, we show a striking genome-wide increase in aberrant somatic hypermutation in EBV-positive tumors, supporting a link between EBV and activation-induced cytidine deaminase (AICDA) activity. In addition to identifying novel candidate BL genes such as SIN3A, USP7, and CHD8, we demonstrate that EBV-positive tumors had significantly fewer driver mutations, especially among genes with roles in apoptosis. We also found immunoglobulin variable region genes that were disproportionally used to encode clonal B-cell receptors (BCRs) in the tumors. These include IGHV4-34, known to produce autoreactive antibodies, and IGKV3-20, a feature described in other B-cell malignancies but not yet in BL. Our results suggest that tumor EBV status defines a specific BL phenotype irrespective of geographic origin, with particular molecular properties and distinct pathogenic mechanisms. The novel mutation patterns identified here imply rational use of DNA-damaging chemotherapy in some patients with BL and targeted agents such as the CDK4/6 inhibitor palbociclib in others, whereas the importance of BCR signaling in BL strengthens the potential benefit of inhibitors for PI3K, Syk, and Src family kinases among these patients.

Transcriptomic and proteomic host response to Aspergillus fumigatus conidia in an air-liquid interface model of human bronchial epithelium.

Aspergillus fumigatus (A. fumigatus) is a wide-spread fungus that is a potent allergen in hypersensitive individuals but also an opportunistic pathogen in immunocompromised patients. It reproduces asexually by releasing airborne conidiospores (conidia). Upon inhalation, fungal conidia are capable of reaching the airway epithelial cells (AECs) in bronchial and alveolar tissues. Previous studies have predominantly used submerged monolayer cultures for studying this host-pathogen interaction; however, these cultures do not recapitulate the mucocililary differentiation phenotype of the in vivo epithelium in the respiratory tract. Thus, the aim of this study was to use well-differentiated primary human bronchial epithelial cells (HBECs) grown at the air-liquid interface (ALI) to determine their transcriptomic and proteomic responses following interaction with A. fumigatus conidia. We visualized conidial interaction with HBECs using confocal laser scanning microscopy (CLSM), and applied NanoString nCounter and shotgun proteomics to assess gene expression changes in the human cells upon interaction with A. fumigatus conidia. Western blot analysis was used to assess the expression of top three differentially expressed proteins, CALR, SET and NUCB2. CLSM showed that, unlike submerged monolayer cultures, well-differentiated ALI cultures of primary HBECs were estimated to internalize less than 1% of bound conidia. Nevertheless, transcriptomic and proteomic analyses revealed numerous differentially expressed host genes; these were enriched for pathways including apoptosis/autophagy, translation, unfolded protein response and cell cycle (up-regulated); complement and coagulation pathways, iron homeostasis, nonsense mediated decay and rRNA binding (down-regulated). CALR and SET were confirmed to be up-regulated in ALI cultures of primary HBECs upon exposure to A. fumigatus via western blot analysis. Therefore, using transcriptomics and proteomics approaches, ALI models recapitulating the bronchial epithelial barrier in the conductive zone of the respiratory tract can provide novel insights to the molecular response of bronchial epithelial cells upon exposure to A. fumigatus conidia.

Genome-wide discovery of somatic regulatory variants in diffuse large B-cell lymphoma.

Diffuse large B-cell lymphoma (DLBCL) is an aggressive cancer originating from mature B-cells. Prognosis is strongly associated with molecular subgroup, although the driver mutations that distinguish the two main subgroups remain poorly defined. Through an integrative analysis of whole genomes, exomes, and transcriptomes, we have uncovered genes and non-coding loci that are commonly mutated in DLBCL. Our analysis has identified novel cis-regulatory sites, and implicates recurrent mutations in the 3' UTR of NFKBIZ as a novel mechanism of oncogene deregulation and NF-κB pathway activation in the activated B-cell (ABC) subgroup. Small amplifications associated with over-expression of FCGR2B (the Fcγ receptor protein IIB), primarily in the germinal centre B-cell (GCB) subgroup, correlate with poor patient outcomes suggestive of a novel oncogene. These results expand the list of subgroup driver mutations that may facilitate implementation of improved diagnostic assays and could offer new avenues for the development of targeted therapeutics.

Interactions of Aspergillus fumigatus Conidia with Airway Epithelial Cells: A Critical Review.

Aspergillus fumigatus is an environmental filamentous fungus that also acts as an opportunistic pathogen able to cause a variety of symptoms, from an allergic response to a life-threatening disseminated fungal infection. The infectious agents are inhaled conidia whose first point of contact is most likely to be an airway epithelial cell (AEC). The interaction between epithelial cells and conidia is multifaceted and complex, and has implications for later steps in pathogenesis. Increasing evidence has demonstrated a key role for the airway epithelium in the response to respiratory pathogens, particularly at early stages of infection; therefore, elucidating the early stages of interaction of conidia with AECs is essential to understand the establishment of infection in cohorts of at-risk patients. Here, we present a comprehensive review of the early interactions between A. fumigatus and AECs, including bronchial and alveolar epithelial cells. We describe mechanisms of adhesion, internalization of conidia by AECs, the immune response of AECs, as well as the role of fungal virulence factors, and patterns of fungal gene expression characteristic of early infection. A clear understanding of the mechanisms involved in the early establishment of infection by A. fumigatus could point to novel targets for therapy and prophylaxis.

Systems Biology Approaches for Host-Fungal Interactions: An Expanding Multi-Omics Frontier.

Opportunistic fungal infections are an increasing threat for global health, and for immunocompromised patients in particular. These infections are characterized by interaction between fungal pathogen and host cells. The exact mechanisms and the attendant variability in host and fungal pathogen interaction remain to be fully elucidated. The field of systems biology aims to characterize a biological system, and utilize this knowledge to predict the system's response to stimuli such as fungal exposures. A multi-omics approach, for example, combining data from genomics, proteomics, metabolomics, would allow a more comprehensive and pan-optic "two systems" biology of both the host and the fungal pathogen. In this review and literature analysis, we present highly specialized and nascent methods for analysis of multiple -omes of biological systems, in addition to emerging single-molecule visualization techniques that may assist in determining biological relevance of multi-omics data. We provide an overview of computational methods for modeling of gene regulatory networks, including some that have been applied towards the study of an interacting host and pathogen. In sum, comprehensive characterizations of host-fungal pathogen systems are now possible, and utilization of these cutting-edge multi-omics strategies may yield advances in better understanding of both host biology and fungal pathogens at a systems scale.