Evolutionary histories of breast cancer and related clones.

Recent studies have documented frequent evolution of clones carrying common cancer mutations in apparently normal tissues, which are implicated in cancer development1-3. However, our knowledge is still missing with regard to what additional driver events take place in what order, before one or more of these clones in normal tissues ultimately evolve to cancer. Here, using phylogenetic analyses of multiple microdissected samples from both cancer and non-cancer lesions, we show unique evolutionary histories of breast cancers harbouring der(1;16), a common driver alteration found in roughly 20% of breast cancers. The approximate timing of early evolutionary events was estimated from the mutation rate measured in normal epithelial cells. In der(1;16)(+) cancers, the derivative chromosome was acquired from early puberty to late adolescence, followed by the emergence of a common ancestor by the patient's early 30s, from which both cancer and non-cancer clones evolved. Replacing the pre-existing mammary epithelium in the following years, these clones occupied a large area within the premenopausal breast tissues by the time of cancer diagnosis. Evolution of multiple independent cancer founders from the non-cancer ancestors was common, contributing to intratumour heterogeneity. The number of driver events did not correlate with histology, suggesting the role of local microenvironments and/or epigenetic driver events. A similar evolutionary pattern was also observed in another case evolving from an AKT1-mutated founder. Taken together, our findings provide new insight into how breast cancer evolves.

Two Aldehyde Clearance Systems Are Essential to Prevent Lethal Formaldehyde Accumulation in Mice and Humans.

Reactive aldehydes arise as by-products of metabolism and are normally cleared by multiple families of enzymes. We find that mice lacking two aldehyde detoxifying enzymes, mitochondrial ALDH2 and cytoplasmic ADH5, have greatly shortened lifespans and develop leukemia. Hematopoiesis is disrupted profoundly, with a reduction of hematopoietic stem cells and common lymphoid progenitors causing a severely depleted acquired immune system. We show that formaldehyde is a common substrate of ALDH2 and ADH5 and establish methods to quantify elevated blood formaldehyde and formaldehyde-DNA adducts in tissues. Bone-marrow-derived progenitors actively engage DNA repair but also imprint a formaldehyde-driven mutation signature similar to aging-associated human cancer mutation signatures. Furthermore, we identify analogous genetic defects in children causing a previously uncharacterized inherited bone marrow failure and pre-leukemic syndrome. Endogenous formaldehyde clearance alone is therefore critical for hematopoiesis and in limiting mutagenesis in somatic tissues.

Landscape and function of multiple mutations within individual oncogenes.

Sporadic reports have described cancer cases in which multiple driver mutations (MMs) occur in the same oncogene1,2. However, the overall landscape and relevance of MMs remain elusive. Here we carried out a pan-cancer analysis of 60,954 cancer samples, and identified 14 pan-cancer and 6 cancer-type-specific oncogenes in which MMs occur more frequently than expected: 9% of samples with at least one mutation in these genes harboured MMs. In various oncogenes, MMs are preferentially present in cis and show markedly different mutational patterns compared with single mutations in terms of type (missense mutations versus in-frame indels), position and amino-acid substitution, suggesting a cis-acting effect on mutational selection. MMs show an overrepresentation of functionally weak, infrequent mutations, which confer enhanced oncogenicity in combination. Cells with MMs in the PIK3CA and NOTCH1 genes exhibit stronger dependencies on the mutated genes themselves, enhanced downstream signalling activation and/or greater sensitivity to inhibitory drugs than those with single mutations. Together oncogenic MMs are a relatively common driver event, providing the underlying mechanism for clonal selection of suboptimal mutations that are individually rare but collectively account for a substantial proportion of oncogenic mutations.

Germ line DDX41 mutations define a unique subtype of myeloid neoplasms.

Germ line DDX41 variants have been implicated in late-onset myeloid neoplasms (MNs). Despite an increasing number of publications, many important features of DDX41-mutated MNs remain to be elucidated. Here we performed a comprehensive characterization of DDX41-mutated MNs, enrolling a total of 346 patients with DDX41 pathogenic/likely-pathogenic (P/LP) germ line variants and/or somatic mutations from 9082 MN patients, together with 525 first-degree relatives of DDX41-mutated and wild-type (WT) patients. P/LP DDX41 germ line variants explained ∼80% of known germ line predisposition to MNs in adults. These risk variants were 10-fold more enriched in Japanese MN cases (n = 4461) compared with the general population of Japan (n = 20 238). This enrichment of DDX41 risk alleles was much more prominent in male than female (20.7 vs 5.0). P/LP DDX41 variants conferred a large risk of developing MNs, which was negligible until 40 years of age but rapidly increased to 49% by 90 years of age. Patients with myelodysplastic syndromes (MDS) along with a DDX41-mutation rapidly progressed to acute myeloid leukemia (AML), which was however, confined to those having truncating variants. Comutation patterns at diagnosis and at progression to AML were substantially different between DDX41-mutated and WT cases, in which none of the comutations affected clinical outcomes. Even TP53 mutations made no exceptions and their dismal effect, including multihit allelic status, on survival was almost completely mitigated by the presence of DDX41 mutations. Finally, outcomes were not affected by the conventional risk stratifications including the revised/molecular International Prognostic Scoring System. Our findings establish that MDS with DDX41-mutation defines a unique subtype of MNs that is distinct from other MNs.

Age-related remodelling of oesophageal epithelia by mutated cancer drivers.

Clonal expansion in aged normal tissues has been implicated in the development of cancer. However, the chronology and risk dependence of the expansion are poorly understood. Here we intensively sequence 682 micro-scale oesophageal samples and show, in physiologically normal oesophageal epithelia, the progressive age-related expansion of clones that carry mutations in driver genes (predominantly NOTCH1), which is substantially accelerated by alcohol consumption and by smoking. Driver-mutated clones emerge multifocally from early childhood and increase their number and size with ageing, and ultimately replace almost the entire oesophageal epithelium in the extremely elderly. Compared with mutations in oesophageal cancer, there is a marked overrepresentation of NOTCH1 and PPM1D mutations in physiologically normal oesophageal epithelia; these mutations can be acquired before late adolescence (as early as early infancy) and significantly increase in number with heavy smoking and drinking. The remodelling of the oesophageal epithelium by driver-mutated clones is an inevitable consequence of normal ageing, which-depending on lifestyle risks-may affect cancer development.

Cutibacterium acnes invades prostate epithelial cells to induce BRCAness as a possible pathogen of prostate cancer.

BACKGROUND: Abundant evidence suggests that chronic inflammation is linked to prostate cancer and that infection is a possible cause of prostate cancer. METHODS: To identify microbiota or pathogens associated with prostate cancer, we investigated the transcriptomes of 20 human prostate cancer tissues. We performed de novo assembly of nonhuman sequences from RNA-seq data. RESULTS: We identified four bacteria as candidate microbiota in the prostate, including Moraxella osloensis, Uncultured chroococcidiopsis, Cutibacterium acnes, and Micrococcus luteus. Among these, C. acnes was detected in 19 of 20 prostate cancer tissue samples by immunohistochemistry. We then analyzed the gene expression profiles of prostate epithelial cells infected in vitro with C. acnes and found significant changes in homologous recombination (HR) and the Fanconi anemia pathway. Notably, electron microscopy demonstrated that C. acnes invaded prostate epithelial cells and localized in perinuclear vesicles, whereas analysis of γH2AX foci and HR assays demonstrated impaired HR repair. In particular, BRCA2 was significantly downregulated in C. acnes-infected cells. CONCLUSIONS: These findings suggest that C. acnes infection in the prostate could lead to HR deficiency (BRCAness) which promotes DNA double-strand breaks, thereby increasing the risk of cancer development.

Whole-genome landscape of adult T-cell leukemia/lymphoma.

Adult T-cell leukemia/lymphoma (ATL) is an aggressive neoplasm immunophenotypically resembling regulatory T cells, associated with human T-cell leukemia virus type-1. Here, we performed whole-genome sequencing (WGS) of 150 ATL cases to reveal the overarching landscape of genetic alterations in ATL. We discovered frequent (33%) loss-of-function alterations preferentially targeting the CIC long isoform, which were overlooked by previous exome-centric studies of various cancer types. Long but not short isoform-specific inactivation of Cic selectively increased CD4+CD25+Foxp3+ T cells in vivo. We also found recurrent (13%) 3'-truncations of REL, which induce transcriptional upregulation and generate gain-of-function proteins. More importantly, REL truncations are also common in diffuse large B-cell lymphoma, especially in germinal center B-cell-like subtype (12%). In the non-coding genome, we identified recurrent mutations in regulatory elements, particularly splice sites, of several driver genes. In addition, we characterized the different mutational processes operative in clustered hypermutation sites within and outside immunoglobulin/T-cell receptor genes and identified the mutational enrichment at the binding sites of host and viral transcription factors, suggesting their activities in ATL. By combining the analyses for coding and noncoding mutations, structural variations, and copy number alterations, we discovered 56 recurrently altered driver genes, including 11 novel ones. Finally, ATL cases were classified into 2 molecular groups with distinct clinical and genetic characteristics based on the driver alteration profile. Our findings not only help to improve diagnostic and therapeutic strategies in ATL, but also provide insights into T-cell biology and have implications for genome-wide cancer driver discovery.

Prognostic significance of chronic kidney disease and impaired renal function in Japanese patients with COVID-19.

BACKGROUND: Renal impairment is a predictor of coronavirus disease (COVID-19) severity. No studies have compared COVID-19 outcomes in patients with chronic kidney disease (CKD) and patients with impaired renal function without a prior diagnosis of CKD. This study aimed to identify the impact of pre-existing impaired renal function without CKD on COVID-19 outcomes. METHODS: This retrospective study included 3,637 patients with COVID-19 classified into three groups by CKD history and estimated glomerular filtration rate (eGFR) on referral: Group 1 (n = 2,460), normal renal function without a CKD history; Group 2 (n = 905), impaired renal function without a CKD history; and Group 3 (n = 272), history of CKD. We compared the clinical characteristics of these groups and assessed the effect of CKD and impaired renal function on critical outcomes (requirement for respiratory support with high-flow oxygen devices, invasive mechanical ventilation, or extracorporeal membrane oxygen, and death during hospitalization) using multivariable logistic regression. RESULTS: The prevalence of comorbidities (hypertension, diabetes, and cardiovascular disease) and incidence of inflammatory responses (white blood counts, and C-reactive protein, procalcitonin, and D-dimer levels) and complications (bacterial infection and heart failure) were higher in Groups 2 and 3 than that in Group 1. The incidence of critical outcomes was 10.8%, 17.7%, and 26.8% in Groups 1, 2, and 3, respectively. The mortality rate and the rate of requiring IMV support was lowest in Group 1 and highest in Group 3. Compared with Group 1, the risk of critical outcomes was higher in Group 2 (adjusted odds ratio [aOR]: 1.32, 95% confidence interval [CI]: 1.03-1.70, P = 0.030) and Group 3 (aOR: 1.94, 95% CI: 1.36-2.78, P < 0.001). Additionally, the eGFR was significantly associated with critical outcomes in Groups 2 (odds ratio [OR]: 2.89, 95% CI: 1.64-4.98, P < 0.001) and 3 (OR: 1.87, 95% CI: 1.08-3.23, P = 0.025) only. CONCLUSIONS: Clinicians should consider pre-existing CKD and impaired renal function at the time of COVID-19 diagnosis for the management of COVID-19.

Differential Gene Regulatory Network Analysis between Azacitidine-Sensitive and -Resistant Cell Lines.

Azacitidine, a DNA methylation inhibitor, is employed for the treatment of acute myeloid leukemia (AML). However, drug resistance remains a major challenge for effective azacitidine chemotherapy, though several studies have attempted to uncover the mechanisms of azacitidine resistance. With the aim to identify the mechanisms underlying acquired azacitidine resistance in cancer cell lines, we developed a computational strategy that can identify differentially regulated gene networks between drug-sensitive and -resistant cell lines by extending the existing method, differentially coexpressed gene sets (DiffCoEx). The technique specifically focuses on cell line-specific gene network analysis. We applied our method to gene networks specific to azacitidine sensitivity and identified differentially regulated gene networks between azacitidine-sensitive and -resistant cell lines. The molecular interplay between the metallothionein gene family, C19orf33, ELF3, GRB7, IL18, NRN1, and RBM47 were identified as differentially regulated gene network in drug resistant cell lines. The biological mechanisms associated with azacitidine and AML for the markers in the identified networks were verified through the literature. Our results suggest that controlling the identified genes (e.g., the metallothionein gene family) and "cellular response"-related pathways ("cellular response to zinc ion", "cellular response to copper ion", and "cellular response to cadmium ion", where the enriched functional-related genes are MT2A, MT1F, MT1G, and MT1E) may provide crucial clues to address azacitidine resistance in patients with AML. We expect that our strategy will be a useful tool to uncover patient-specific molecular interplay that provides crucial clues for precision medicine in not only gastric cancer but also complex diseases.

UBTF-internal tandem duplication as a novel poor prognostic factor in pediatric acute myeloid leukemia.

The prognosis of pediatric acute myeloid leukemia (AML) has improved via stratification therapy. However, relapse or death occurs in 30%-40% of cases. Novel genetic factors for pediatric AML need to be elucidated to improve prognosis. We detected recurrent internal tandem duplication in upstream binding transcription factor (UBTF-ITD) in 1.2% (6/503) of Japanese pediatric patients with de novo AML. No UBTF-ITD was detected in 175 adult patients with AML or in 65 cell lines that included 15 AML, 39 acute lymphoblastic leukemia, five chronic myeloid leukemia, and six neuroblastoma cell lines. All UBTF-ITDs were found in exon 13 and shared a duplicated region. UBTF-ITD was more frequently detected in patients with trisomy 8, FLT3-ITD, WT1 mutation, and/or high PRDM16 expression (trisomy 8, 3/6; FLT3-ITD, 5/6; WT1 mutation, 2/6; and high PRDM16 expression, 6/6). Gene expression patterns of patients with UBTF-ITD were similar to those of patients with NUP98::NSD1 or FUS::ERG. Survival analysis of the AML-05 cohort revealed that patients with UBTF-ITD had worse outcomes than those without UBTF-ITD (3-year event-free survival, 20% vs. 55%; 3-year overall survival, 40% vs. 74%). Moreover, among the 27 patients with trisomy 8, all three patients with UBTF -ITD had a poor prognosis resulting in early events (relapse or non-complete remission) within 1 year. Our findings suggest that UBTF-ITD may be a novel and significant prognostic factor for pediatric patients with AML.

TP53 and RB1 alterations characterize poor prognostic subgroups in pediatric acute myeloid leukemia.

Pediatric acute myeloid leukemia (AML) is a poor prognostic subtype of pediatric leukemia. However, the detailed characteristics of many genetic abnormalities are yet to be established in this disease. Although TP53 and RB1 are established as representative tumor suppressor genes in various cancers, alterations of these two genes, especially RB1, have not been characterized in pediatric AML. We performed next-generation sequencing in 328 pediatric AML patients from the Japanese AML-05 trial to ascertain TP53 and RB1 alterations, and their prognostic implications. We identified seven patients with TP53 alterations (2.1%) and six patients with RB1 alterations (1.8%). These alterations were found in only patients without RUNX1::RUNX1T1, CBFB::MYH11, or KMT2A rearrangements. TP53 and RB1 were frequently co-deleted with their neighboring genes PRPF8 and ELF1, respectively. Patients with TP53 alterations had significantly lower 5-year overall survival (OS; 14.3% vs. 71.4%, p < 0.001) and lower 5-year event-free survival (EFS; 0% vs. 56.3%, p < 0.001); similarly, patients with RB1 had significantly lower 5-year OS (0% vs. 71.8%, p < 0.001) and lower 5-year EFS (0% vs. 56.0%, p < 0.001) when compared to patients without these alterations. In gene expression analyses, oxidative phosphorylation, glycolysis, and protein secretion were upregulated in patients with TP53 and/or RB1 alterations. Additionally, Kaplan-Meier analysis revealed that high expressions of SLC2A5, KCNAB2, and CD300LF were related to poor OS of non-core-binding factor AML patients (p < 0.001, p = 0.001, and p = 0.021, respectively). This study will contribute to the development of risk-stratified therapy and precision medicine in pediatric AML.

Convergent genomic diversity and novel BCAA metabolism in intrahepatic cholangiocarcinoma.

BACKGROUND: Driver alterations may represent novel candidates for driver gene-guided therapy; however, intrahepatic cholangiocarcinoma (ICC) with multiple genomic aberrations makes them intractable. Therefore, the pathogenesis and metabolic changes of ICC need to be understood to develop new treatment strategies. We aimed to unravel the evolution of ICC and identify ICC-specific metabolic characteristics to investigate the metabolic pathway associated with ICC development using multiregional sampling to encompass the intra- and inter-tumoral heterogeneity. METHODS: We performed the genomic, transcriptomic, proteomic and metabolomic analysis of 39-77 ICC tumour samples and eleven normal samples. Further, we analysed their cell proliferation and viability. RESULTS: We demonstrated that intra-tumoral heterogeneity of ICCs with distinct driver genes per case exhibited neutral evolution, regardless of their tumour stage. Upregulation of BCAT1 and BCAT2 indicated the involvement of 'Val Leu Ile degradation pathway'. ICCs exhibit the accumulation of ubiquitous metabolites, such as branched-chain amino acids including valine, leucine, and isoleucine, to negatively affect cancer prognosis. We revealed that this metabolic pathway was almost ubiquitously altered in all cases with genomic diversity and might play important roles in tumour progression and overall survival. CONCLUSIONS: We propose a novel ICC onco-metabolic pathway that could enable the development of new therapeutic interventions.

Divergent lncRNA MYMLR regulates MYC by eliciting DNA looping and promoter-enhancer interaction.

Long non-coding RNAs (lncRNAs) function in a wide range of processes by diverse mechanisms, though their roles in regulation of oncogenes and/or tumor suppressors remain rather elusive. We performed a global search for lncRNAs affecting MYC activity using a systems biology-based approach with a K supercomputer and the GIMLET algorism based on local distance correlations. Consequently, MYMLR was identified and experimentally shown to maintain MYC transcriptional activity and cell cycle progression despite the low levels of expression. A proteomic search for MYMLR-binding proteins identified PCBP2, while it was also found that MYMLR places a 557-kb upstream enhancer region in the proximity of the MYC promoter in cooperation with PCBP2. These findings implicate a crucial role for MYMLR in regulation of the archetypical oncogene MYC and warrant future studies regarding the involvement of low copy number lncRNAs in regulation of other crucial oncogenes and tumor suppressor genes.

Bone Marrow Failure and Immunodeficiency Associated with Human RAD50 Variants.

PURPOSE: The MRE11-RAD50-NBN (MRN) complex plays a key role in recognizing and signaling DNA double-strand breaks. Pathogenic variants in NBN and MRE11 give rise to the autosomal-recessive diseases, Nijmegen breakage syndrome (NBS) and ataxia telangiectasia-like disorder, respectively. The clinical consequences of pathogenic variants in RAD50 are incompletely understood. We aimed to characterize a newly identified RAD50 deficiency/NBS-like disorder (NBSLD) patient with bone marrow failure and immunodeficiency. METHODS: We report on a girl with microcephaly, mental retardation, bird-like face, short stature, bone marrow failure and B-cell immunodeficiency. We searched for candidate gene by whole-exome sequencing and analyzed the cellular phenotype of patient-derived fibroblasts using immunoblotting, radiation sensitivity assays and lentiviral complementation experiments. RESULTS: Compound heterozygosity for two variants in the RAD50 gene (p.Arg83His and p.Glu485Ter) was identified in this patient. The expression of RAD50 protein and MRN complex formation was maintained in the cells derived from this patient. DNA damage-induced activation of the ATM kinase was markedly decreased, which was restored by the expression of wild-type (WT) RAD50. Radiosensitivity appeared inconspicuous in the patient-derived cell line as assessed by colony formation assay. The RAD50R83H missense substitution did not rescue the mitotic defect in complementation experiments using RAD50-deficient fibroblasts, whereas RAD50WT did. The RAD50E485X nonsense variant was associated with in-frame skipping of exon 10 (p.Glu485_545del). CONCLUSION: These findings indicate important roles of RAD50 in human bone marrow and immune cells. RAD50 deficiency/NBSLD can manifest as a distinct inborn error of immunity characterized by bone marrow failure and B-cell immunodeficiency.

Comprehensive analysis of indels in whole-genome microsatellite regions and microsatellite instability across 21 cancer types.

Microsatellites are repeats of 1- to 6-bp units, and approximately 10 million microsatellites have been identified across the human genome. Microsatellites are vulnerable to DNA mismatch errors and have thus been used to detect cancers with mismatch repair deficiency. To reveal the mutational landscape of microsatellite repeat regions at the genome level, we analyzed approximately 20.1 billion microsatellites in 2717 whole genomes of pan-cancer samples across 21 tissue types. First, we developed a new insertion and deletion caller (MIMcall) that takes into consideration the error patterns of different types of microsatellites. Among the 2717 pan-cancer samples, our analysis identified 31 samples, including colorectal, uterus, and stomach cancers, with a higher proportion of mutated microsatellite (≥0.03), which we defined as microsatellite instability (MSI) cancers of genome-wide level. Next, we found 20 highly mutated microsatellites that can be used to detect MSI cancers with high sensitivity. Third, we found that replication timing and DNA shape were significantly associated with mutation rates of microsatellites. Last, analysis of mutations in mismatch repair genes showed that somatic SNVs and short indels had larger functional impacts than germline mutations and structural variations. Our analysis provides a comprehensive picture of mutations in the microsatellite regions and reveals possible causes of mutations, as well as provides a useful marker set for MSI detection.

Genetic analysis of low-grade adenosquamous carcinoma of the breast progressing to high-grade metaplastic carcinoma.

PURPOSE: Low-grade adenosquamous carcinoma (LGASC) is a rare type of metaplastic carcinoma of the breast (MBC) with an indolent clinical course. A few LGASC cases with high-grade transformation have been reported; however, the genetics underlying malignant progression of LGASC remain unclear. METHODS: We performed whole-genome sequencing analysis on five MBCs from four patients, including one case with matching primary LGASC and a lymph node metastatic tumor consisting of high-grade MBC with a predominant metaplastic squamous cell carcinoma component (MSC) that progressed from LGASC and three cases of independent de novo MSC. RESULTS: Unlike de novo MSC, LGASC and its associated MSC showed no TP53 mutation and tended to contain fewer structural variants than de novo MSC. Both LGASC and its associated MSC harbored the common GNAS c.C2530T:p.Arg844Cys mutation, which was more frequently detected in the cancer cell fraction of MSC. MSC associated with LGASC showed additional pathogenic deletions of multiple tumor-suppressor genes, such as KMT2D and BTG1. Copy number analysis revealed potential 18q loss of heterozygosity in both LGASC and associated MSC. The frequency of SMAD4::DCC fusion due to deletions increased with progression to MSC; however, chimeric proteins were not detected. SMAD4 protein expression was already decreased at the LGASC stage due to unknown mechanisms. CONCLUSION: Not only LGASC but also its associated high-grade MBC may be genetically different from de novo high-grade MBC. Progression from LGASC to high-grade MBC may involve the concentration of driver mutations caused by clonal selection and inactivation of tumor-suppressor genes.

Role of the Orphan Transporter SLC35E1 in the Nuclear Egress of Herpes Simplex Virus 1.

This study developed a system consisting of two rounds of screening cellular proteins involved in the nuclear egress of herpes simplex virus 1 (HSV-1). Using this system, we first screened cellular proteins that interacted with the HSV-1 nuclear egress complex (NEC) consisting of UL34 and UL31 in HSV-1-infected cells, which are critical for the nuclear egress of HSV-1, by tandem affinity purification coupled with mass spectrometry-based proteomics technology. Next, we performed CRISPR/Cas9-based screening of live HSV-1-infected reporter cells under fluorescence microscopy using single guide RNAs targeting the cellular proteins identified in the first proteomic screening to detect the mislocalization of the lamin-associated protein emerin, which is a phenotype for defects in HSV-1 nuclear egress. This study focused on a cellular orphan transporter SLC35E1, one of the cellular proteins identified by the screening system. Knockout of SLC35E1 reduced HSV-1 replication and induced membranous invaginations containing perinuclear enveloped virions (PEVs) adjacent to the nuclear membrane (NM), aberrant accumulation of PEVs in the perinuclear space between the inner and outer NMs and the invagination structures, and mislocalization of the NEC. These effects were similar to those of previously reported mutation(s) in HSV-1 proteins and depletion of cellular proteins that are important for HSV-1 de-envelopment, one of the steps required for HSV-1 nuclear egress. Our newly established screening system enabled us to identify a novel cellular protein required for efficient HSV-1 de-envelopment. IMPORTANCE The identification of cellular protein(s) that interact with viral effector proteins and function in important viral procedures is necessary for enhancing our understanding of the mechanics of various viral processes. In this study, we established a new system consisting of interactome screening for the herpes simplex virus 1 (HSV-1) nuclear egress complex (NEC), followed by loss-of-function screening to target the identified putative NEC-interacting cellular proteins to detect a defect in HSV-1 nuclear egress. This newly established system identified SLC35E1, an orphan transporter, as a novel cellular protein required for efficient HSV-1 de-envelopment, providing an insight into the mechanisms involved in this viral procedure.

RoDiCE: robust differential protein co-expression analysis for cancer complexome.

MOTIVATION: The full spectrum of abnormalities in cancer-associated protein complexes remains largely unknown. Comparing the co-expression structure of each protein complex between tumor and healthy cells may provide insights regarding cancer-specific protein dysfunction. However, the technical limitations of mass spectrometry-based proteomics, including contamination with biological protein variants, causes noise that leads to non-negligible over- (or under-) estimating co-expression. RESULTS: We propose a robust algorithm for identifying protein complex aberrations in cancer based on differential protein co-expression testing. Our method based on a copula is sufficient for improving identification accuracy with noisy data compared to conventional linear correlation-based approaches. As an application, we use large-scale proteomic data from renal cancer to show that important protein complexes, regulatory signaling pathways and drug targets can be identified. The proposed approach surpasses traditional linear correlations to provide insights into higher-order differential co-expression structures. AVAILABILITY AND IMPLEMENTATION: https://github.com/ymatts/RoDiCE. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Identification of bacteriophage genome sequences with representation learning.

MOTIVATION: Bacteriophages/phages are the viruses that infect and replicate within bacteria and archaea, and rich in human body. To investigate the relationship between phages and microbial communities, the identification of phages from metagenome sequences is the first step. Currently, there are two main methods for identifying phages: database-based (alignment-based) methods and alignment-free methods. Database-based methods typically use a large number of sequences as references; alignment-free methods usually learn the features of the sequences with machine learning and deep learning models. RESULTS: We propose INHERIT which uses a deep representation learning model to integrate both database-based and alignment-free methods, combining the strengths of both. Pre-training is used as an alternative way of acquiring knowledge representations from existing databases, while the BERT-style deep learning framework retains the advantage of alignment-free methods. We compare INHERIT with four existing methods on a third-party benchmark dataset. Our experiments show that INHERIT achieves a better performance with the F1-score of 0.9932. In addition, we find that pre-training two species separately helps the non-alignment deep learning model make more accurate predictions. AVAILABILITY AND IMPLEMENTATION: The codes of INHERIT are now available in: https://github.com/Celestial-Bai/INHERIT. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Frequent genetic alterations in immune checkpoint-related genes in intravascular large B-cell lymphoma.

Intravascular large B-cell lymphoma (IVLBCL) is a unique type of extranodal lymphoma characterized by selective growth of tumor cells in small vessels without lymphadenopathy. Greater understanding of the molecular pathogenesis of IVLBCL is hampered by the paucity of lymphoma cells in biopsy specimens, creating a limitation in obtaining sufficient tumor materials. To uncover the genetic landscape of IVLBCL, we performed whole-exome sequencing (WES) of 21 patients with IVLBCL using plasma-derived cell-free DNA (cfDNA) (n = 18), patient-derived xenograft tumors (n = 4), and tumor DNA from bone marrow (BM) mononuclear cells (n = 2). The concentration of cfDNA in IVLBCL was significantly higher than that in diffuse large B-cell lymphoma (DLBCL) (P < .0001) and healthy donors (P = .0053), allowing us to perform WES; most mutations detected in BM tumor DNA were successfully captured in cfDNA and xenograft. IVLBCL showed a high frequency of genetic lesions characteristic of activated B-cell-type DLBCL, with the former showing conspicuously higher frequencies (compared with nodal DLBCL) of mutations in MYD88 (57%), CD79B (67%), SETD1B (57%), and HLA-B (57%). We also found that 8 IVLBCL (38%) harbored rearrangements of programmed cell death 1 ligand 1 and 2 (PD-L1/PD-L2) involving the 3' untranslated region; such rearrangements are implicated in immune evasion via PD-L1/PD-L2 overexpression. Our data demonstrate the utility of cfDNA and imply important roles for immune evasion in IVLBCL pathogenesis and PD-1/PD-L1/PD-L2 blockade in therapeutics for IVLBCL.