Using VarScan 2 for Germline Variant Calling and Somatic Mutation Detection

Integrative molecular and spatial analysis reveals evolutionary dynamics and tumor-immune interplay of in situ and invasive acral melanoma

In acral melanoma (AM), progression from in situ (AMis) to invasive AM (iAM) leads to significantly reduced survival. However, evolutionary dynamics during this process remain elusive. Here, we report integrative molecular and spatial characterization of 147 AMs using genomics, bulk and single-cell transcriptomics, and spatial transcriptomics and proteomics. Vertical invasion from AMis to iAM displays an early and monoclonal seeding pattern. The subsequent regional expansion of iAM exhibits two distinct patterns, clonal expansion and subclonal diversification. Notably, molecular subtyping reveals an aggressive iAM subset featured with subclonal diversification, increased epithelial-mesenchymal transition (EMT), and spatial enrichment of APOE+/CD163+ macrophages. In vitro and ex vivo experiments further demonstrate that APOE+CD163+ macrophages promote tumor EMT via IGF1-IGF1R interaction. Adnexal involvement can predict AMis with higher invasive potential whereas APOE and CD163 serve as prognostic biomarkers for iAM. Altogether, our results provide implications for the early detection and treatment of AM.

Mutations in Podospora anserina MCM1 and VelC Trigger Spontaneous Development of Barren Fruiting Bodies

The ascomycete Podospora anserina is a heterothallic filamentous fungus found mainly on herbivore dung. It is commonly used in laboratories as a model system, and its complete life cycle lasting eight days is well mastered in vitro. The main objective of our team is to understand better the global process of fruiting body development, named perithecia, induced normally in this species by fertilization. Three allelic mutants, named pfd3, pfd9, and pfd23 (for "promoting fruiting body development") obtained by UV mutagenesis, were selected in view of their abilities to promote barren perithecium development without fertilization. By complete genome sequencing of pfd3 and pfd9, and mutant complementation, we identified point mutations in the mcm1 gene as responsible for spontaneous perithecium development. MCM1 proteins are MADS box transcription factors that control diverse developmental processes in plants, metazoans, and fungi. We also identified using the same methods a mutation in the VelC gene as responsible for spontaneous perithecium development in the vacua mutant. The VelC protein belongs to the velvet family of regulators involved in the control of development and secondary metabolite production. A key role of MCM1 and VelC in coordinating the development of P. anserina perithecia with gamete formation and fertilization is highlighted.

Multiple DNA repair pathways prevent acetaldehyde-induced mutagenesis in yeast

Acetaldehyde is the primary metabolite of alcohol and is present in many environmental sources including tobacco smoke. Acetaldehyde is genotoxic, whereby it can form DNA adducts and lead to mutagenesis. Individuals with defects in acetaldehyde clearance pathways have increased susceptibility to alcohol-associated cancers. Moreover, a mutation signature specific to acetaldehyde exposure is widespread in alcohol and smoking-associated cancers. However, the pathways that repair acetaldehyde-induced DNA damage and thus prevent mutagenesis are vaguely understood. Here, we used Saccharomyces cerevisiae to systematically delete genes in each of the major DNA repair pathways to identify those that alter acetaldehyde-induced mutagenesis. We found that deletion of the nucleotide excision repair (NER) genes, RAD1 or RAD14, led to an increase in mutagenesis upon acetaldehyde exposure. Acetaldehyde-induced mutations were dependent on translesion synthesis as well as DNA inter-strand crosslink (ICL) repair in Δrad1 strains. Moreover, whole genome sequencing of the mutated isolates demonstrated an increase in C→A changes coupled with an enrichment of gCn→A changes in the acetaldehyde-treated Δrad1 isolates. The gCn→A mutation signature has been shown to be diagnostic of acetaldehyde exposure in yeast and in human cancers. We also demonstrated that the deletion of the two DNA-protein crosslink (DPC) repair proteases, WSS1 and DDI1, also led to increased acetaldehyde-induced mutagenesis. Defects in base excision repair (BER) led to a mild increase in mutagenesis, while defects in mismatch repair (MMR), homologous recombination repair (HR) and post replicative repair pathways did not impact mutagenesis upon acetaldehyde exposure. Our results in yeast were further corroborated upon analysis of whole exome sequenced liver cancers, wherein, tumors with defects in ERCC1 and ERCC4 (NER), FANCD2 (ICL repair) or SPRTN (DPC repair) carried a higher gCn→A mutation load than tumors with no deleterious mutations in these genes. Our findings demonstrate that multiple DNA repair pathways protect against acetaldehyde-induced mutagenesis.

Combined DNA Analysis from Stool and Blood Samples Improves Tumor Tracking and Assessment of Clonal Heterogeneity in Localized Rectal Cancer Patients

Objectives: In this study, stool samples were evaluated for tumor mutation analysis via a targeted next generation sequencing (NGS) approach in a small patient cohort suffering from localized rectal cancer. Introduction: Colorectal cancer (CRC) causes the second highest cancer-related death rate worldwide. Thus, improvements in disease assessment and monitoring that may facilitate treatment allocation and allow organ-sparing "watch-and-wait" treatment strategies are highly relevant for a significant number of CRC patients. Methods: Stool-based results were compared with mutation profiles derived from liquid biopsies and the gold standard procedure of tumor biopsy from the same patients. A workflow was established that enables the detection of de-novo tumor mutations in stool samples of CRC patients via ultra-sensitive cell-free tumor DNA target enrichment. Results: Notably, only a 19% overall concordance was found in mutational profiles across the compared sample specimens of stool, tumor, and liquid biopsies. Conclusion: Based on these results, the analysis of stool and liquid biopsy samples can provide important additional information on tumor heterogeneity and potentially on the assessment of minimal residual disease and clonal tumor evolution.

Overexpression of MTHFD2 represents an inflamed tumor microenvironment and precisely predicts the molecular subtype and immunotherapy response of bladder cancer

Introduction

Methylenetetrahydrofolate dehydrogenase 2 (MTHFD2), whose aberrant expression is common in cancers, has recently been identified as a potential regulator of immune response. However, its immune-related role in bladder cancer (BLCA) and its association with immunotherapy efficacy remain unclear.

Methods

RNA sequencing data from The Cancer Genome Atlas (TCGA) was applied to analyze the immunological roles and prognostic value of MTHFD2 in pan-cancers. The association of MTHFD2 with several immunological features of tumor microenvironment (TME), including cancer-immunity cycle, immune cells infiltration, immune checkpoints expression, and T cell inflamed score was analyzed in TCGA-BLCA cohort. The predictors of cancer treatments effectiveness, including the expression and mutation of certain genes, molecular subtypes, and several signatures were evaluated as well. These results were validated by another independent cohort (GSE48075). Finally, the predictive value of MTHFD2 for TME and immunotherapy efficacy were validated using immunohistochemistry assay and RNA sequencing data from IMvigor210 cohort, respectively.

Results

MTHFD2 was found to be positively associated with several immunological features of an inflamed tumor microenvironment (TME) in various cancers and could predict BLCA patients' prognosis. In BLCA, high expression of MTHFD2 was observed to be positively related with the cancer-immunity cycle, the infiltration of several immune cells, and the expression of immunoregulators and T-cell inflamed scores, indicating a positive correlation with the inflamed TME. Moreover, patients with high MTHFD2 expression were more likely to be basal-like subtypes and respond to BLCA treatments, including immunotherapy, chemotherapy, and target therapy. The clinical data of the IMvigor210 cohort confirmed the higher response rates and better survival benefits of immunotherapy in high-MTHFD2-expression patients.

Conclusion

Collectively, high MTHFD2 predicts an inflamed TME, a basal-like subtype, and a better response to various therapeutic strategies, especially the ICB therapy, in bladder cancer.

Base-Excision Repair Mutational Signature in Two Sebaceous Carcinomas of the Eyelid

Personalized medicine aims to develop tailored treatments for individual patients based on specific mutations present in the affected organ. This approach has proven paramount in cancer treatment, as each tumor carries distinct driver mutations that respond to targeted drugs and, in some cases, may confer resistance to other therapies. Particularly for rare conditions, personalized medicine has the potential to revolutionize treatment strategies. Rare cancers often lack extensive datasets of molecular and pathological information, large-scale trials for novel therapies, and established treatment guidelines. Consequently, surgery is frequently the only viable option for many rare tumors, when feasible, as traditional multimodal approaches employed for more common cancers often play a limited role. Sebaceous carcinoma of the eyelid is an exceptionally rare cancer affecting the eye's adnexal tissues, most frequently reported in Asia, but whose prevalence is significantly increasing even in Europe and the US. The sole established curative treatment is surgical excision, which can lead to significant disfigurement. In cases of metastatic sebaceous carcinoma, validated drug options are currently lacking. In this project, we set out to characterize the mutational landscape of two sebaceous carcinomas of the eyelid following surgical excision. Utilizing available bioinformatics tools, we demonstrated our ability to identify common features promptly and accurately in both tumors. These features included a Base-Excision Repair mutational signature, a notably high tumor mutational burden, and key driver mutations in somatic tissues. These findings had not been previously reported in similar studies. This report underscores how, in the case of rare tumors, it is possible to comprehensively characterize the mutational landscape of each individual case, potentially opening doors to targeted therapeutic options.

Identifying rare genetic variants in 21 highly multiplex autism families: the role of diagnosis and autistic traits

Autism is a highly heritable, heterogeneous, neurodevelopmental condition. Large-scale genetic studies, predominantly focussing on simplex families and clinical diagnoses of autism have identified hundreds of genes associated with autism. Yet, the contribution of these classes of genes to multiplex families and autistic traits still warrants investigation. Here, we conducted whole-genome sequencing of 21 highly multiplex autism families, with at least three autistic individuals in each family, to prioritise genes associated with autism. Using a combination of both autistic traits and clinical diagnosis of autism, we identify rare variants in genes associated with autism, and related neurodevelopmental conditions in multiple families. We identify a modest excess of these variants in autistic individuals compared to individuals without an autism diagnosis. Finally, we identify a convergence of the genes identified in molecular pathways related to development and neurogenesis. In sum, our analysis provides initial evidence to demonstrate the value of integrating autism diagnosis and autistic traits to prioritise genes.

Validation of genetic variants from NGS data using deep convolutional neural networks

Accurate somatic variant calling from next-generation sequencing data is one most important tasks in personalised cancer therapy. The sophistication of the available technologies is ever-increasing, yet, manual candidate refinement is still a necessary step in state-of-the-art processing pipelines. This limits reproducibility and introduces a bottleneck with respect to scalability. We demonstrate that the validation of genetic variants can be improved using a machine learning approach resting on a Convolutional Neural Network, trained using existing human annotation. In contrast to existing approaches, we introduce a way in which contextual data from sequencing tracks can be included into the automated assessment. A rigorous evaluation shows that the resulting model is robust and performs on par with trained researchers following published standard operating procedure.

High resolution genomes of multiple Xiphophorus species provide new insights into microevolution, hybrid incompatibility, and epistasis

Because of diverged adaptative phenotypes, fish species of the genus Xiphophorus have contributed to a wide range of research for a century. Existing Xiphophorus genome assemblies are not at the chromosomal level and are prone to sequence gaps, thus hindering advancement of the intra- and inter-species differences for evolutionary, comparative, and translational biomedical studies. Herein, we assembled high-quality chromosome-level genome assemblies for three distantly related Xiphophorus species, namely, X. maculatus, X. couchianus, and X. hellerii Our overall goal is to precisely assess microevolutionary processes in the clade to ascertain molecular events that led to the divergence of the Xiphophorus species and to progress understanding of genetic incompatibility to disease. In particular, we measured intra- and inter-species divergence and assessed gene expression dysregulation in reciprocal interspecies hybrids among the three species. We found expanded gene families and positively selected genes associated with live bearing, a special mode of reproduction. We also found positively selected gene families are significantly enriched in nonpolymorphic transposable elements, suggesting the dispersal of these nonpolymorphic transposable elements has accompanied the evolution of the genes, possibly by incorporating new regulatory elements in support of the Britten-Davidson hypothesis. We characterized inter-specific polymorphisms, structural variants, and polymorphic transposable element insertions and assessed their association to interspecies hybridization-induced gene expression dysregulation related to specific disease states in humans.

Newly identified sex chromosomes in the Sphagnum (peat moss) genome alter carbon sequestration and ecosystem dynamics

Peatlands are crucial sinks for atmospheric carbon but are critically threatened due to warming climates. Sphagnum (peat moss) species are keystone members of peatland communities where they actively engineer hyperacidic conditions, which improves their competitive advantage and accelerates ecosystem-level carbon sequestration. To dissect the molecular and physiological sources of this unique biology, we generated chromosome-scale genomes of two Sphagnum species: S. divinum and S. angustifolium. Sphagnum genomes show no gene colinearity with any other reference genome to date, demonstrating that Sphagnum represents an unsampled lineage of land plant evolution. The genomes also revealed an average recombination rate an order of magnitude higher than vascular land plants and short putative U/V sex chromosomes. These newly described sex chromosomes interact with autosomal loci that significantly impact growth across diverse pH conditions. This discovery demonstrates that the ability of Sphagnum to sequester carbon in acidic peat bogs is mediated by interactions between sex, autosomes and environment.

DNA damage and somatic mutations in mammalian cells after irradiation with a nail polish dryer

Ultraviolet A light is commonly emitted by UV-nail polish dryers with recent reports suggesting that long-term use may increase the risk for developing skin cancer. However, no experimental evaluation has been conducted to reveal the effect of radiation emitted by UV-nail polish dryers on mammalian cells. Here, we show that irradiation by a UV-nail polish dryer causes high levels of reactive oxygen species, consistent with 8-oxo-7,8-dihydroguanine damage and mitochondrial dysfunction. Analysis of somatic mutations reveals a dose-dependent increase of C:G>A:T substitutions in irradiated samples with mutagenic patterns similar to mutational signatures previously attributed to reactive oxygen species. In summary, this study demonstrates that radiation emitted by UV-nail polish dryers can both damage DNA and permanently engrave mutations on the genomes of primary mouse embryonic fibroblasts, human foreskin fibroblasts, and human epidermal keratinocytes.

Temporal and spatial stability of the EM/PM molecular subtypes in adult diffuse glioma

Detailed characterizations of genomic alterations have not identified subtype-specific vulnerabilities in adult gliomas. Mapping gliomas into developmental programs may uncover new vulnerabilities that are not strictly related to genomic alterations. After identifying conserved gene modules co-expressed with EGFR or PDGFRA (EM or PM), we recently proposed an EM/PM classification scheme for adult gliomas in a histological subtype- and grade-independent manner. By using cohorts of bulk samples, paired primary and recurrent samples, multi-region samples from the same glioma, single-cell RNA-seq samples, and clinical samples, we here demonstrate the temporal and spatial stability of the EM and PM subtypes. The EM and PM subtypes, which progress in a subtype-specific mode, are robustly maintained in paired longitudinal samples. Elevated activities of cell proliferation, genomic instability and microenvironment, rather than subtype switching, mark recurrent gliomas. Within individual gliomas, the EM/PM subtype was preserved across regions and single cells. Malignant cells in the EM and PM gliomas were correlated to neural stem cell and oligodendrocyte progenitor cell compartment, respectively. Thus, while genetic makeup may change during progression and/or within different tumor areas, adult gliomas evolve within a neurodevelopmental framework of the EM and PM molecular subtypes. The dysregulated developmental pathways embedded in these molecular subtypes may contain subtype-specific vulnerabilities.

Oral cancer prediction by noninvasive genetic screening

Oral squamous cell carcinomas (OSCCs) develop in genetically altered epithelium in the mucosal lining, also coined as fields, which are mostly not visible but occasionally present as white oral leukoplakia (OL) lesions. We developed a noninvasive genetic assay using next-generation sequencing (NGS) on brushed cells to detect the presence of genetically altered fields, including those that are not macroscopically visible. The assay demonstrated high accuracy in OL patients when brush samples were compared with biopsies as gold standard. In a cohort of Fanconi anemia patients, detection of mutations in prospectively collected oral brushes predicted oral cancer also when visible abnormalities were absent. We further provide insight in the molecular landscape of OL with frequent changes of TP53, FAT1 and NOTCH1. NGS analysis of noninvasively collected samples offers a highly accurate method to detect genetically altered fields in the oral cavity, and predicts development of OSCC in high-risk individuals. Noninvasive genetic screening can be employed to screen high-risk populations for cancer and precancer, map the extension of OL lesions beyond what is visible, map the oral cavity for precancerous changes even when visible abnormalities are absent, test accuracy of promising imaging modalities, monitor interventions and determine genetic progression as well as the natural history of the disease in the human patient.

De novo variants are a common cause of genetic hearing loss

PURPOSE

De novo variants (DNVs) are a well-recognized cause of genetic disorders. The contribution of DNVs to hearing loss (HL) is poorly characterized. We aimed to evaluate the rate of DNVs in HL-associated genes and assess their contribution to HL.

METHODS

Targeted genomic enrichment and massively parallel sequencing were used for molecular testing of all exons and flanking intronic sequences of known HL-associated genes, with no exclusions on the basis of type of HL or clinical features. Segregation analysis was performed, and previous reports of DNVs in PubMed and ClinVar were reviewed to characterize the rate, distribution, and spectrum of DNVs in HL.

RESULTS

DNVs were detected in 10% (24/238) of trios for whom segregation analysis was performed. Overall, DNVs were causative in at least ∼1% of probands for whom a genetic diagnosis was resolved, with marked variability based on inheritance mode and phenotype. DNVs of MITF were most common (21% of DNVs), followed by GATA3 (13%), STRC (13%), and ACTG1 (8%). Review of reported DNVs revealed gene-specific variability in contribution of DNV to the mutational spectrum of HL-associated genes.

CONCLUSION

DNVs are a relatively common cause of genetic HL and must be considered in all cases of sporadic HL.

Development of a high-throughput SNP array for sea cucumber (Apostichopus japonicus) and its application in genomic selection with MCP regularized deep neural networks

High-throughput single nucleotide polymorphism (SNP) genotyping assays are powerful tools for genetic studies and genomic breeding applications for many species. Though large numbers of SNPs have been identified in sea cucumber (Apostichopus japonicus), but, as yet, no high-throughput genotyping platform is available for this species. In this study, we designed and developed a high-throughput 24 K SNP genotyping array named HaishenSNP24K for A. japonicus, based on the multi-objective-local optimization (MOLO) algorithm and HD-Marker genotyping method. The SNP array exhibited a relatively high genotyping call rate (> 96%), genotyping accuracy (>95%) and exhibited highly polymorphic in sea cucumber populations. In addition, we also assessed its application in genomic selection (GS). Deep neural networks (DNN) that can capture the complicated interactions of genes have been proposed as a promising tool in GS for SNP-based genomic prediction of complex traits in animal breeding. To overcome the problem of over-fitting when using the HaishenSNP24K array as high-dimensional DNN input, we developed minmax concave penalty (MCP) regularization for sparse deep neural networks (DNN-MCP) that finds an optimal sparse structure of a DNN by minimizing the square error subject to the non-convex penalty MCP on the parameters (weights and biases). Compared to two linear models, namely RR-GBLUP and Bayes B, and the nonlinear model DNN, DNN-MCP has greatly improved the genomic prediction ability for three quantitative traits (e.g., wet weight, dry weight and survival time) in the sea cucumber population. To the best of our knowledge, this is the first work to develop a high-throughput SNP array for A. japonicus and a new model DNN-MCP for genomic prediction of complex traits in GS. The present results provide evidence that supports the HaishenSNP24K array with DNN-MCP will be valuable for genetic studies and molecular breeding in A. japonicus.

Circulating Tumor DNA Characteristics Based on Next Generation Sequencing and Its Correlation With Clinical Parameters in Patients With Lymphoma

Background

Lymphoma is a heterogeneous group of tumors in terms of morphological subtypes, molecular alterations, and management. However, data on circulating tumor DNA (ctDNA) mutated genes are limited. The purpose of this study was to investigate the features of the ctDNA mutated genes, the prognosis, and the association between the ctDNA mutated genes and the clinical parameters in lymphoma.

Methods

Differences in the ctDNA between the mutated genes and the prognosis of 59 patients with Hodgkin’s lymphoma (HL) (10.2%), germinal center B-cell–like lymphoma (GCB) (28.8%), nongerminal center B-cell–like lymphoma (non-GCB) (50.8%), and marginal zone lymphoma (MZL) (10.2%) were analyzed by next generation sequencing (NGS) targeting 121 lymphoma-relevant genes.

Results

Genetic alterations were identified in the ctDNA samples with a median of 6 variants per sample. The genetic variation of the ctDNA in the plasma was found to be significantly correlated with the clinical indices in lymphoma. The genetic heterogeneity of different lymphoma subtypes was clearly observed in the ctDNAs from HL, GCB, non-GCB, and MZL, confirming that distinct molecular mechanisms are involved in the pathogenesis of different lymphomas.

Conclusion

Our findings suggest that NGS-based ctDNA mutation analysis reveals genetic heterogeneity across lymphoma subtypes, with potential implications for discovering therapeutic targets, exploring genomic evolution, and developing risk-adaptive therapies.

Integrative genomic analysis reveals low T-cell infiltration as the primary feature of tobacco use in HPV-positive oropharyngeal cancer

Although tobacco use is an independent adverse prognostic feature in HPV(+) oropharyngeal squamous cell carcinoma (OPSCC), the biologic features associated with tobacco use have not been systematically investigated. We characterized genomic and immunologic features associated with tobacco use through whole exome sequencing, mRNA hybridization, and immunohistochemical staining in 47 HPV(+) OPSCC tumors. Low expression of transcripts in a T cell-inflamed gene expression profile (TGEP) was associated with tobacco use at diagnosis and lower overall and disease-free survival. Tobacco use was associated with an increased proportion of T > C substitutions and a lower proportion of expected mutational signatures, but not with increases in mutational burden or recurrent oncogenic mutations. Our findings suggest that rather than increased mutational burden, tobacco's primary and clinically relevant association in HPV(+) OPSCC is immunosuppression of the tumor immune microenvironment. Quantitative assays of T cell infiltration merit further study as prognostic markers in HPV(+) OPSCC.

TGS1 mediates 2,2,7-trimethyl guanosine capping of the human telomerase RNA to direct telomerase dependent telomere maintenance

Pathways that direct the selection of the telomerase-dependent or recombination-based, alternative lengthening of telomere (ALT) maintenance pathway in cancer cells are poorly understood. Using human lung cancer cells and tumor organoids we show that formation of the 2,2,7-trimethylguanosine (TMG) cap structure at the human telomerase RNA 5′ end by the Trimethylguanosine Synthase 1 (TGS1) is central for recruiting telomerase to telomeres and engaging Cajal bodies in telomere maintenance. TGS1 depletion or inhibition by the natural nucleoside sinefungin impairs telomerase recruitment to telomeres leading to Exonuclease 1 mediated generation of telomere 3′ end protrusions that engage in RAD51-dependent, homology directed recombination and the activation of key features of the ALT pathway. This indicates a critical role for 2,2,7-TMG capping of the RNA component of human telomerase (hTR) in enforcing telomerase-dependent telomere maintenance to restrict the formation of telomeric substrates conductive to ALT. Our work introduces a targetable pathway of telomere maintenance that holds relevance for telomere-related diseases such as cancer and aging.

Telomerase protects chromosome ends in stem cells and cancer cells. Here the authors show that Trimethylguaonsine Synthase 1 (TGS-1) – dependent trimethylguanosine capping of the RNA component of the human telomerase complex has an important role in directing telomere dependent telomere maintenance and suppressing the ALT pathway in cancer cells.

Comparison of Next-Generation Sequencing and Polymerase Chain Reaction for Personalized Treatment-Related Genomic Status in Patients with Metastatic Colorectal Cancer

Personalized treatments based on the genetic profiles of tumors can simultaneously optimize efficacy and minimize toxicity, which is beneficial for improving patient outcomes. This study aimed to integrate gene alterations associated with predictive and prognostic outcomes in patients with metastatic colorectal cancer (mCRC) with polymerase chain reaction (PCR) and in-house next-generation sequencing (NGS) to detect KRAS, NRAS, and BRAF mutations. In the present study, 41 patients with mCRC were assessed between August 2017 and June 2019 at a single institution. The overall concordance between NGS and PCR results for detecting KRAS, NRAS, and BRAF mutations was considerably high (87.8–92.7%), with only 15 discrepant results between PCR and NGS. Our companion diagnostic test analyzes KRAS, NRAS, and BRAF as a panel of CRC molecular targets; therefore, it has the advantages of requiring fewer specimens and being more time and cost efficient than conventional testing for separate analyses, allowing for the simultaneous analysis of multiple genes.

Mutation profiling of circulating tumor DNA identifies distinct mutation patterns in non-Hodgkin lymphoma

OBJECTIVE

Circulating tumor DNA (ctDNA) is emerging as a versatile biomarker for noninvasive genotyping and response monitoring in specific B-cell lymphomas; however, few studies have been conducted to explore ctDNA-based mutation profiling across non-Hodgkin lymphomas (NHLs) and genomic changes after initiation of chemotherapy.

METHODS

A targeted sequencing of 362 genes was performed to detect the mutation profiles in paired blood and tissue samples from 42 NHL patients. Genomic alterations were explored in 11 diffuse large B-cell lymphoma (DLBCL) patients using paired blood samples collected pre- and post-R-CHOP chemotherapy.

RESULTS

The frequencies of PIM1, MYD88, MYC, ZNF292, JAK, and MAF mutations were higher in aggressive than in indolent B-cell lymphoma and NK/T subtypes. Tumor mutation burden in blood samples was higher in aggressive than in indolent B-cell lymphomas and higher in patients who progressed than in those who responded to treatments. Our data also revealed significant enhance of concordance index through integrating mutated genes that were significantly associated with prognosis into International Prognostic Index-based prognostic model. Moreover, acquisition of mutations such as PCLO_p.L1220Tfs*3 was associated with resistance to R-CHOP in DLBCL patients.

CONCLUSIONS

Our findings illustrated distinct mutation patterns across various NHL subtypes and suggested the association of genomic alterations in ctDNA with treatment outcomes.

Circulating Cell-Free DNA Reflects the Clonal Evolution of Breast Cancer Tumors

Simple Summary

Liquid biopsy of cell-free DNA (cfDNA) is proposed as potential method for the early detection of breast cancer (BC) metastases and following the clonal evolution of BC. Though the use of liquid biopsy is widely discussed, only few studies have demonstrated such usage so far. The aim of this study was to evaluate how accurately cfDNA resembles the genetic profile of tumor DNA and how liquid biopsy reflects the clonal evolution of 18 Eastern-Finnish BC cases with locoregional or distant metastases. Although notable discordance between the sequenced cfDNA and tumor DNA was observed, our results show that liquid biopsy reflects the heterogeneity and clonal evolution of BC and may help to identify potential driver variants and therapeutic targets that are not detected with the sequencing of tumor DNA. This information may be used to detect metastatic BC earlier and to support decision-making in clinical practice.

Abstract

Liquid biopsy of cell-free DNA (cfDNA) is proposed as a potential method for the early detection of breast cancer (BC) metastases and following the clonal evolution of BC. Though the use of liquid biopsy is a widely discussed topic in the field, only a few studies have demonstrated such usage so far. We sequenced the DNA of matched primary tumor and metastatic sites together with the matched cfDNA samples from 18 Eastern Finnish BC patients and investigated how well cfDNA reflected the clonal evolution of BC interpreted from tumor DNA. On average, liquid biopsy detected 56.2 ± 7.2% of the somatic variants that were present either in the matched primary tumor or metastatic sites. Despite the high discordance observed between matched samples, liquid biopsy was found to reflect the clonal evolution of BC and identify novel driver variants and therapeutic targets absent from the tumor DNA. Tumor-specific somatic variants were detected in cfDNA at the time of diagnosis and 8.4 ± 2.4 months prior to detection of locoregional recurrence or distant metastases. Our results demonstrate that the sequencing of cfDNA may be used for the early detection of locoregional and distant BC metastases. Observed discordance between tumor DNA sequencing and liquid biopsy supports the parallel sequencing of cfDNA and tumor DNA to yield the most comprehensive overview for the genetic landscape of BC.

Salt stress responses and SNP-based phylogenetic analysis of Thai rice cultivars

Genetic diversity is important for developing salt-tolerant rice (Oryza sativa L.) cultivars. Certain Thai rice accessions display salt tolerance at the adult or reproductive stage, but their response to salinity at the seedling stage is unknown. In this study, a total of 10 rice cultivars/line, including eight Thai cultivars and standard salt-tolerant cultivar and susceptible line, were screened using a hydroponic system to identify salt-tolerant genotypes at the seedling stage. Different morphophysiological indicators were used to classify tolerant and susceptible genotypes. Phylogenetic analyses were performed to determine the evolutionary relationships between the cultivars. Results showed that 'Lai Mahk', 'Jao Khao', 'Luang Pratahn', and 'Ma Gawk' exhibited salt stress tolerance comparable with the standard salt-tolerance check 'Pokkali'. Whole-exome single-nucleotide polymorphism (SNP)-based phylogenetic analysis showed that the Thai rice cultivars were monophyletic and distantly related to Pokkali and IR29. Lai Mahk and Luang Pratahn were found closely related when using the whole-exome SNPs for the analysis. This is also the case for the analysis of SNPs from 164 salt-tolerance genes and transcription regulatory genes. The tolerant cultivars shared the same haplotype for 16 genes. Overall, the findings of this study identified four salt-stress-tolerant Thai rice cultivars, which could be used in rice breeding programs for salinity tolerance.

GCAT|Panel, a comprehensive structural variant haplotype map of the Iberian population from high-coverage whole-genome sequencing

The combined analysis of haplotype panels with phenotype clinical cohorts is a common approach to explore the genetic architecture of human diseases. However, genetic studies are mainly based on single nucleotide variants (SNVs) and small insertions and deletions (indels). Here, we contribute to fill this gap by generating a dense haplotype map focused on the identification, characterization, and phasing of structural variants (SVs). By integrating multiple variant identification methods and Logistic Regression Models (LRMs), we present a catalogue of 35 431 441 variants, including 89 178 SVs (≥50 bp), 30 325 064 SNVs and 5 017 199 indels, across 785 Illumina high coverage (30x) whole-genomes from the Iberian GCAT Cohort, containing a median of 3.52M SNVs, 606 336 indels and 6393 SVs per individual. The haplotype panel is able to impute up to 14 360 728 SNVs/indels and 23 179 SVs, showing a 2.7-fold increase for SVs compared with available genetic variation panels. The value of this panel for SVs analysis is shown through an imputed rare Alu element located in a new locus associated with Mononeuritis of lower limb, a rare neuromuscular disease. This study represents the first deep characterization of genetic variation within the Iberian population and the first operational haplotype panel to systematically include the SVs into genome-wide genetic studies.

Patient-specific TBX5-G125R Variant Induces Profound Transcriptional Deregulation and Atrial Dysfunction

Background: The pathogenic missense variant p.G125R in TBX5 causes Holt-Oram syndrome (HOS; hand-heart syndrome) and early onset of atrial fibrillation. Revealing how an altered key developmental transcription factor modulates cardiac physiology in vivo will provide unique insights into the mechanisms underlying atrial fibrillation in these patients. Methods: We analyzed electrocardiograms (ECGs) of an extended family pedigree of HOS patients. Next, we introduced the TBX5-p.G125R variant in the mouse genome (Tbx5^G125R) and performed electrophysiological analyses (ECG, optical mapping, patch clamp, intracellular calcium measurements), transcriptomics (single nuclei and tissue RNA sequencing) and epigenetic profiling (ATAC-sequencing, H3K27ac CUT&RUN-sequencing). Results: We discovered high incidence of atrial extra systoles and atrioventricular conduction disturbances in HOS patients. Tbx5^G125R/+ mice were morphologically unaffected and displayed variable RR intervals, atrial extra systoles and susceptibility to atrial fibrillation, reminiscent of TBX5-p.G125R patients. Atrial conduction velocity was not affected but systolic and diastolic intracellular calcium concentrations were decreased and action potentials prolonged in isolated cardiomyocytes of Tbx5^G125R/+ mice compared to controls. Transcriptional profiling of atria revealed most profound transcriptional changes in cardiomyocytes versus other cell types, and identified over a thousand coding and non-coding transcripts that were differentially expressed. Epigenetic profiling uncovered thousands of TBX5-p.G125R sensitive putative regulatory elements (including enhancers) that gained accessibility in atrial cardiomyocytes. The majority of sites with increased accessibility were occupied by Tbx5. The small group of sites with reduced accessibility was enriched for DNA binding motifs of members of the SP- and KLF families of transcription factors. These data show that Tbx5-p.G125R induces changes in regulatory element activity, altered transcriptional regulation and changed cardiomyocyte behavior, possibly caused by altered DNA binding and cooperativity properties. Conclusions: Our data reveal how a disease-causing missense variant in TBX5 induces profound changes in the atrial transcriptional regulatory network and epigenetic state in vivo, leading to arrhythmia reminiscent of those seen in human TBX5-p.G125R variant carriers.

Comprehensive Molecular Landscape of Cetuximab Resistance in Head and Neck Cancer Cell Lines

Cetuximab is the sole anti-EGFR monoclonal antibody that is FDA approved to treat head and neck squamous cell carcinoma (HNSCC). However, no predictive biomarkers of cetuximab response are known for HNSCC. Herein, we address the molecular mechanisms underlying cetuximab resistance in an in vitro model. We established a cetuximab resistant model (FaDu), using increased cetuximab concentrations for more than eight months. The resistance and parental cells were evaluated for cell viability and functional assays. Protein expression was analyzed by Western blot and human cell surface panel by lyoplate. The mutational profile and copy number alterations (CNA) were analyzed using whole-exome sequencing (WES) and the NanoString platform. FaDu resistant clones exhibited at least two-fold higher IC50 compared to the parental cell line. WES showed relevant mutations in several cancer-related genes, and the comparative mRNA expression analysis showed 36 differentially expressed genes associated with EGFR tyrosine kinase inhibitors resistance, RAS, MAPK, and mTOR signaling. Importantly, we observed that overexpression of KRAS, RhoA, and CD44 was associated with cetuximab resistance. Protein analysis revealed EGFR phosphorylation inhibition and mTOR increase in resistant cells. Moreover, the resistant cell line demonstrated an aggressive phenotype with a significant increase in adhesion, the number of colonies, and migration rates. Overall, we identified several molecular alterations in the cetuximab resistant cell line that may constitute novel biomarkers of cetuximab response such as mTOR and RhoA overexpression. These findings indicate new strategies to overcome anti-EGFR resistance in HNSCC.

Identification and Validation of a Novel Six-Gene Expression Signature for Predicting Hepatocellular Carcinoma Prognosis

Background

Hepatocellular carcinoma (HCC) is a highly lethal disease. Effective prognostic tools to guide clinical decision-making for HCC patients are lacking.

Objective

We aimed to establish a robust prognostic model based on differentially expressed genes (DEGs) in HCC.

Methods

Using datasets from The Cancer Genome Atlas (TCGA), the Gene Expression Omnibus (GEO), and the International Genome Consortium (ICGC), DEGs between HCC tissues and adjacent normal tissues were identified. Using TCGA dataset as the training cohort, we applied the least absolute shrinkage and selection operator (LASSO) algorithm and multivariate Cox regression analyses to identify a multi-gene expression signature. Proportional hazard assumptions and multicollinearity among covariates were evaluated while building the model. The ICGC cohort was used for validation. The Pearson test was used to evaluate the correlation between tumor mutational burden and risk score. Through single-sample gene set enrichment analysis, we investigated the role of signature genes in the HCC microenvironment.

Results

A total of 274 DEGs were identified, and a six-DEG prognostic model was developed. Patients were stratified into low- or high-risk groups based on risk scoring by the model. Kaplan-Meier analysis revealed significant differences in overall survival and progression-free interval. Through univariate and multivariate Cox analyses, the model proved to be an independent prognostic factor compared to other clinic-pathological parameters. Time-dependent receiver operating characteristic curve analysis revealed satisfactory prediction of overall survival, but not progression-free interval. Functional enrichment analysis showed that cancer-related pathways were enriched, while immune infiltration analyses differed between the two risk groups. The risk score did not correlate with levels of PD-1, PD-L1, CTLA4, or tumor mutational burden.

Conclusions

We propose a six-gene expression signature that could help to determine HCC patient prognosis. These genes may serve as biomarkers in HCC and support personalized disease management.

Role and Regulation of Poly-3-Hydroxybutyrate in Nitrogen Fixation in Azorhizobium caulinodans

An Azorhizobium caulinodans phaC mutant (OPS0865) unable to make poly-3-hydroxybutyrate (PHB), grows poorly on many carbon sources and cannot fix nitrogen in laboratory culture. However, when inoculated onto its host plant, Sesbania rostrata, the phaC mutant consistently fixed nitrogen. Upon reisolation from S. rostrata root nodules, a suppressor strain (OPS0921) was isolated that has significantly improved growth on a variety of carbon sources and also fixes nitrogen in laboratory culture. The suppressor retains the original mutation and is unable to synthesize PHB. Genome sequencing revealed a suppressor transition mutation, G to A (position 357,354), 13 bases upstream of the ATG start codon of phaR in its putative ribosome binding site (RBS). PhaR is the global regulator of PHB synthesis but also has other roles in regulation within the cell. In comparison with the wild type, translation from the phaR native RBS is increased approximately sixfold in the phaC mutant background, suggesting that the level of PhaR is controlled by PHB. Translation from the phaR mutated RBS (RBS*) of the suppressor mutant strain (OPS0921) is locked at a low basal rate and unaffected by the phaC mutation, suggesting that RBS* renders the level of PhaR insensitive to regulation by PHB. In the original phaC mutant (OPS0865), the lack of nitrogen fixation and poor growth on many carbon sources is likely to be due to increased levels of PhaR causing dysregulation of its complex regulon, because PHB formation, per se, is not required for effective nitrogen fixation in A. caulinodans.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Clinical Characteristics of Patients and Whole Genome Sequencing-Based Surveillance of Escherichia coli Community-Onset Bloodstream Infections at a Non-tertiary Hospital in CHINA

Background:Escherichia coli is the most common pathogens in patients with community-onset blood stream infections (COBSI). Knowledge of the epidemiology of this disease is crucial to improve allocation of health resources, formulate isolation strategies that prevent transmission, and guide empirical antibiotic therapy.

Methods: This retrospective observational study examined patients with E. coli COBSI (EC-COBSI) at a non-tertiary hospital in China. Whole-genome sequencing and analysis of the isolates was performed. The relationships of clinical variables with antimicrobial resistance and the genetic background of the isolates were examined.

Results: There were 148 isolates in patients with EC-COBSI. All isolates were susceptible to ceftazidime/avibactam, carbapenems, and tigecycline; 35.1% were positive for extended spectrum β-lactamase (ESBL+); and bla_CTX–M–14 was the most common ESBL gene. Patients with ESBL- isolates were more likely to receive appropriate empiric treatment than those with ESBL+ isolates (61.5% vs. 91.4%, p < 0.001), but these two groups had similar mortality rates. The overall 30-day mortality rate was 9.5%. Phylogenetic analysis showed that the isolates were diverse, and that the main sequence types (STs) were ST95, ST131, and ST69. Intra-abdominal infection was the primary source of disease, and isolates from these patients had lower frequencies of virulence genes.

Conclusion: The mortality rate of patients with EC-COBSI was unrelated to ESBL status of the isolates. Most isolates had low resistance to most of the tested antimicrobial agents. The isolates were diverse, and multiple strains were related. Prevention and control of EC-COBSI should target prevention of patient colonization and the living environment.

Enterovirus A71 Induces Neurological Diseases and Dynamic Variants in Oral Infection of Human SCARB2-Transgenic Weaned Mice

Enterovirus A71 (EV-A71) and many members of the Picornaviridae family are neurotropic pathogens of global concern. These viruses are primarily transmitted through the fecal-oral route, and thus suitable animal models of oral infection are needed to investigate viral pathogenesis. An animal model of oral infection was developed using transgenic mice expressing human SCARB2 (hSCARB2 Tg), murine-adapted EV-A71/MP4 virus, and EV-A71/MP4 virus with an engineered nanoluciferase gene that allows imaging of viral replication and spread in infected mice. Next-generation sequencing of EV-A71 genomes in the tissues and organs of infected mice was also performed. Oral inoculation of EV-A71/MP4 or nanoluciferase-carrying MP4 virus stably induced neurological symptoms and death in infected 21-day-old weaned mice. In vivo bioluminescence imaging of infected mice and tissue immunostaining of viral antigens indicated that orally inoculated virus can spread to the central nervous system (CNS) and other tissues. Next-generating sequencing further identified diverse mutations in viral genomes that can potentially contribute to viral pathogenesis. This study presents an EV-A71 oral infection murine model that efficiently infects weaned mice and allows tracking of viral spread, features that can facilitate research into viral pathogenesis and neuroinvasion via the natural route of infection. IMPORTANCE Enterovirus A71 (EV-A71), a positive-strand RNA virus of the Picornaviridae, poses a persistent global public health problem. EV-A71 is primarily transmitted through the fecal-oral route, and thus suitable animal models of oral infection are needed to investigate viral pathogenesis. We present an animal model of EV-A71 infection that enables the natural route of oral infection in weaned and nonimmunocompromised 21-day-old hSCARB2 transgenic mice. Our results demonstrate that severe disease and death could be stably induced, and viral invasion of the CNS could be replicated in this model, similar to severe real-world EV-A71 infections. We also developed a nanoluciferase-containing EV-A71 virus that can be used with this animal model to track viral spread after oral infection in real time. Such a model offers several advantages over existing animal models and can facilitate future research into viral spread, tissue tropism, and viral pathogenesis, all pressing issues that remain unaddressed for EV-A71 infections.

A POLE Splice Site Deletion Detected in a Patient with Biclonal CLL and Prostate Cancer: A Case Report

Chronic lymphocytic leukemia (CLL) is considered a clonal B cell malignancy. Sporadically, CLL cases with multiple productive heavy and light-chain rearrangements were detected, thus leading to a bi- or oligoclonal CLL disease with leukemic cells originating either from different B cells or otherwise descending from secondary immunoglobulin rearrangement events. This suggests a potential role of clonal hematopoiesis or germline predisposition in these cases. During the screening of 75 CLL cases for kappa and lambda light-chain rearrangements, we could detect a single case with CLL cells expressing two distinct kappa and lambda light chains paired with two separate immunoglobulin heavy-chain variable regions. Furthermore, this patient also developed a prostate carcinoma. Targeted genome sequencing of highly purified light-chain specific CLL clones from this patient and from the prostate carcinoma revealed the presence of a rare germline polymorphism in the POLE gene. Hence, our data suggest that the detected SNP may predispose for cancer, particularly for CLL.

Whole Genome Sequencing of Methicillin-Resistant Staphylococcus epidermidis Clinical Isolates Reveals Variable Composite SCCmec ACME among Different STs in a Tertiary Care Hospital in Oman

Staphylococcus epidermidis has been recently recognized as an emerging nosocomial pathogen. There are concerns over the increasing virulence potential of this commensal due to the capabilities of transferring mobile genetic elements to Staphylococcus aureus through staphylococcal chromosomal cassette (SCCmec) and the closely related arginine catabolic mobile element (ACME) and the copper and mercury resistance island (COMER). The potential pathogenicity of S. epidermidis, particularly from blood stream infections, has been poorly investigated. In this study, 24 S. epidermidis isolated from blood stream infections from Oman were investigated using whole genome sequence analysis. Core genome phylogenetic trees revealed one third of the isolates belong to the multidrug resistance ST-2. Genomic analysis unraveled a common occurrence of SCCmec type IV and ACME element predominantly type I arranged in a composite island. The genetic composition of ACME was highly variable among isolates of same or different STs. The COMER-like island was absent in all of our isolates. Reduced copper susceptibility was observed among isolates of ST-2 and ACME type I, followed by ACME type V. In conclusion, in this work, we identify a prevalent occurrence of highly variable ACME elements in different hospital STs of S. epidermidis in Oman, thus strongly suggesting the hypothesis that ACME types evolved from closely related STs.

Bam-readcount - rapid generation of basepair-resolution sequence metrics

Bam-readcount is a utility for generating low-level information about sequencing data at specific nucleotide positions. Originally designed to help filter genomic mutation calls, the metrics it outputs are useful as input for variant detection tools and for resolving ambiguity between variant callers1,2. In addition, it has found broad applicability in diverse fields including tumor evolution, single-cell genomics, climate change ecology, and tracking community spread of SARS-CoV-2.3-6.

Contributions of cis- and trans-Regulatory Evolution to Transcriptomic Divergence across Populations in the Drosophila mojavensis Larval Brain

Natural selection on gene expression was originally predicted to result primarily in cis- rather than trans-regulatory evolution, due to the expectation of reduced pleiotropy. Despite this, numerous studies have ascribed recent evolutionary divergence in gene expression predominantly to trans-regulation. Performing RNA-seq on single isofemale lines from genetically distinct populations of the cactophilic fly Drosophila mojavensis and their F₁ hybrids, we recapitulated this pattern in both larval brains and whole bodies. However, we demonstrate that improving the measurement of brain expression divergence between populations by using seven additional genotypes considerably reduces the estimate of trans-regulatory contributions to expression evolution. We argue that the finding of trans-regulatory predominance can result from biases due to environmental variation in expression or other sources of noise, and that cis-regulation is likely a greater contributor to transcriptional evolution across D. mojavensis populations. Lastly, we merge these lines of data to identify several previously hypothesized and intriguing novel candidate genes, and suggest that the integration of regulatory and population-level transcriptomic data can provide useful filters for the identification of potentially adaptive genes.

Characterization of a RAD51C-silenced high-grade serous ovarian cancer model during development of PARP inhibitor resistance

Acquired PARP inhibitor (PARPi) resistance in BRCA1- or BRCA2-mutant ovarian cancer often results from secondary mutations that restore expression of functional protein. RAD51C is a less commonly studied ovarian cancer susceptibility gene whose promoter is sometimes methylated, leading to homologous recombination (HR) deficiency and PARPi sensitivity. For this study, the PARPi-sensitive patient-derived ovarian cancer xenograft PH039, which lacks HR gene mutations but harbors RAD51C promoter methylation, was selected for PARPi resistance by cyclical niraparib treatment in vivo. PH039 acquired PARPi resistance by the third treatment cycle and grew through subsequent treatment with either niraparib or rucaparib. Transcriptional profiling throughout the course of resistance development showed widespread pathway level changes along with a marked increase in RAD51C mRNA, which reflected loss of RAD51C promoter methylation. Analysis of ovarian cancer samples from the ARIEL2 Part 1 clinical trial of rucaparib monotherapy likewise indicated an association between loss of RAD51C methylation prior to on-study biopsy and limited response. Interestingly, the PARPi resistant PH039 model remained platinum sensitive. Collectively, these results not only indicate that PARPi treatment pressure can reverse RAD51C methylation and restore RAD51C expression, but also provide a model for studying the clinical observation that PARPi and platinum sensitivity are sometimes dissociated.

Adenine Base Editing Reduces Misfolded Protein Accumulation and Toxicity in Alpha-1 Antitrypsin Deficient Patient iPSC-Hepatocytes

Alpha-1 antitrypsin deficiency (AATD) is most commonly caused by the Z mutation, a single-base substitution that leads to AAT protein misfolding and associated liver and lung disease. In this study, we apply adenine base editors to correct the Z mutation in patient induced pluripotent stem cells (iPSCs) and iPSC-derived hepatocytes (iHeps). We demonstrate that correction of the Z mutation in patient iPSCs reduces aberrant AAT accumulation and increases its secretion. Adenine base editing (ABE) of differentiated iHeps decreases ER stress in edited cells, as demonstrated by single-cell RNA sequencing. We find ABE to be highly efficient in iPSCs and do not identify off-target genomic mutations by whole-genome sequencing. These results reveal the feasibility and utility of base editing to correct the Z mutation in AATD patient cells.

Immunogenomic profiling and pathological response results from a clinical trial of docetaxel and carboplatin in triple-negative breast cancer

PURPOSE

Patients with triple-negative breast cancer (TNBC) who do not achieve pathological complete response (pCR) following neoadjuvant chemotherapy have a high risk of recurrence and death. Molecular characterization may identify patients unlikely to achieve pCR. This neoadjuvant trial was conducted to determine the pCR rate with docetaxel and carboplatin and to identify molecular alterations and/or immune gene signatures predicting pCR.

EXPERIMENTAL DESIGN

Patients with clinical stages II/III TNBC received 6 cycles of docetaxel and carboplatin. The primary objective was to determine if neoadjuvant docetaxel and carboplatin would increase the pCR rate in TNBC compared to historical expectations. We performed whole-exome sequencing (WES) and immune profiling on pre-treatment tumor samples to identify alterations that may predict pCR. Thirteen matching on-treatment samples were also analyzed to assess changes in molecular profiles.

RESULTS

Fifty-eight of 127 (45.7%) patients achieved pCR. There was a non-significant trend toward higher mutation burden for patients with residual cancer burden (RCB) 0/I versus RCB II/III (median 80 versus 68 variants, p 0.88). TP53 was the most frequently mutated gene, observed in 85.7% of tumors. EGFR, RB1, RAD51AP2, SDK2, L1CAM, KPRP, PCDHA1, CACNA1S, CFAP58, COL22A1, and COL4A5 mutations were observed almost exclusively in pre-treatment samples from patients who achieved pCR. Seven mutations in PCDHA1 were observed in pre-treatment samples from patients who did not achieve pCR. Several immune gene signatures including IDO1, PD-L1, interferon gamma signaling, CTLA4, cytotoxicity, tumor inflammation signature, inflammatory chemokines, cytotoxic cells, lymphoid, PD-L2, exhausted CD8, Tregs, and immunoproteasome were upregulated in pre-treatment samples from patients who achieved pCR.

CONCLUSION

Neoadjuvant docetaxel and carboplatin resulted in a pCR of 45.7%. WES and immune profiling differentiated patients with and without pCR.

TRIAL REGISTRATION

Clinical trial information: NCT02124902, Registered 24 April 2014 & NCT02547987, Registered 10 September 2015.

AID Contributes to Accelerated Disease Progression in the TCL1 Mouse Transplant Model for CLL

Simple Summary

Cancers, such as chronic lymphocytic leukemia, frequently acquire consecutive somatic mutations in the genome, which contribute to disease progression and treatment resistance. Activation-induced deaminase is an enzyme responsible for generating the highly diverse B cell repertoire but it can also induce substantial collateral damage within the genome of cells. Hence, it is important to assess whether AID contributes to cancer mutations and to the course of disease. This research shows that AID contributes to the acquisition of somatic cancer-specific mutations in a mouse model for chronic lymphocytic leukemia reflected in prolonged overall survival of leukemic mice lacking AID expression. These data should initiate future studies to assess the effect of AID inhibition on the occurrence of drug resistance.

Abstract

Adaptive somatic mutations conferring treatment resistance and accelerated disease progression is still a major problem in cancer therapy. Additionally in CLL, patients receiving novel, efficient drugs frequently become treatment refractory and eventually relapse. Activation-induced deaminase (AID) is a cytosine deaminase that catalyzes somatic hypermutation of genomic DNA at the immunoglobulin locus in activated B cells. As considerable off-target mutations by AID have been discerned in chronic lymphocytic leukemia, it is essential to investigate to which extent these mutations contribute to disease progression to estimate whether AID inhibition could counteract drug resistance mechanisms. In this study, we examined the TCL1 mouse model for CLL on an AID pro- and deficient background by comparing disease development and mutational landscapes. We provide evidence that AID contributes to the acquisition of somatic cancer-specific mutations also in the TCL1 model and accelerates CLL development particularly in the transplant setting. We conclude that AID is directly determining the fitness of the CLL clone, which prompts further studies to assess the effect of AID inhibition on the occurrence of drug resistance.

The bone microenvironment increases phenotypic plasticity of ER+ breast cancer cells

Estrogen receptor-positive (ER+) breast cancer exhibits a strong bone tropism in metastasis. How the bone microenvironment (BME) impacts ER signaling and endocrine therapy remains poorly understood. Here, we discover that the osteogenic niche transiently and reversibly reduces ER expression and activities specifically in bone micrometastases (BMMs), leading to endocrine resistance. As BMMs progress, the ER reduction and endocrine resistance may partially recover in cancer cells away from the osteogenic niche, creating phenotypic heterogeneity in macrometastases. Using multiple approaches, including an evolving barcoding strategy, we demonstrated that this process is independent of clonal selection, and represents an EZH2-mediated epigenomic reprogramming. EZH2 drives ER+ BMMs toward a basal and stem-like state. EZH2 inhibition reverses endocrine resistance. These data exemplify how epigenomic adaptation to BME promotes phenotypic plasticity of metastatic seeds, fosters intra-metastatic heterogeneity, and alters therapeutic responses. Our study provides insights into the clinical enigma of ER+ metastatic recurrences despite endocrine therapies.

Dissecting Clonal Hematopoiesis in Tissues of Classical Hodgkin Lymphoma Patients

Clonal hematopoiesis predisposes to hematological malignancies. However, clonal hematopoiesis is understudied in classical Hodgkin lymphoma (cHL), a mature B-cell neoplasm exhibiting the most abundant microenvironment. We analyzed clonal hematopoiesis in 40 cHL cases by sequencing microdissected tumor cells and matched normal cells from blood and/or lymph nodes. Five patients had blood and/or tissue clonal hematopoiesis. In three of five patients (all failing first-line therapy), clonal hematopoiesis spread through the tissue microenvironment extensively, and featured mutant DNMT3A^R882H , KRAS^G60D and DNMT3A^R882H +TET2^Q1274 * in 33%, 92% and 60% of non-neoplastic cells, respectively. In the latter case, DNMT3A/TET2-mutant clonal hematopoiesis seeded the neoplastic clone, which was infected by the Epstein-Barr virus and showed almost no other somatic mutations exome-wide. In the former case, DNMT3A-mutant clonal hematopoiesis did not originate the neoplastic clone despite dominating the blood and B-cell lineage (~94% leukocytes; ~96% mature blood B cells), yet led to NPM1-mutated acute myeloid leukemia 6 years after therapy for cHL. Our results expand to cHL the spectrum of hematologic malignancies associated with clonal hematopoiesis.

Transitioning the Molecular Tumor Board from Proof of Concept to Clinical Routine: A German Single-Center Analysis

Molecular precision oncology faces two major challenges: first, to identify relevant and actionable molecular variants in a rapidly changing field and second, to provide access to a broad patient population. Here, we report a four-year experience of the Molecular Tumor Board (MTB) of the Comprehensive Cancer Center Freiburg (Germany) including workflows and process optimizations. This retrospective single-center study includes data on 488 patients enrolled in the MTB from February 2015 through December 2018. Recommendations include individual molecular diagnostics, molecular stratified therapies, assessment of treatment adherence and patient outcomes including overall survival. The majority of MTB patients presented with stage IV oncologic malignancies (90.6%) and underwent an average of 2.1 previous lines of therapy. Individual diagnostic recommendations were given to 487 patients (99.8%). A treatment recommendation was given in 264 of all cases (54.1%) which included a molecularly matched treatment in 212 patients (43.4%). The 264 treatment recommendations were implemented in 76 patients (28.8%). Stable disease was observed in 19 patients (25.0%), 17 had partial response (22.4%) and five showed a complete remission (6.6%). An objective response was achieved in 28.9% of cases with implemented recommendations and for 4.5% of the total population (22 of 488 patients). By optimizing the MTB workflow, case-discussions per session increased significantly while treatment adherence and outcome remained stable over time. Our data demonstrate the feasibility and effectiveness of molecular-guided personalized therapy for cancer patients in a clinical routine setting showing a low but robust and durable disease control rate over time.

Genomic divergence and differential gene expression between crustacean ecotypes across a marine thermal gradient

Environmental gradients between marine biogeographical provinces separate distinct faunal communities. However, the absence of absolute dispersal barriers allows numerous species to occur on both sides of such boundaries. While the regional populations of such widespread species are often morphologically indistinguishable from each other, genetic evidence suggests that they represent unique ecotypes, and likely even cryptic species, that may be uniquely adapted to their local environment. Here, we explored genomic divergence in four sympatric southern African decapod crustaceans whose ranges span the boundary between the cool-temperate west coast (south-eastern Atlantic) and the warm-temperate south coast (south-western Indian Ocean) near the southern tip of the African continent. Using genome-wide data, we found that all four species comprise distinct west coast and south coast ecotypes, with molecular dating suggesting divergence during the Pleistocene. Transcriptomic data from the hepatopancreas of twelve specimens of one of these species, the mudprawn Upogebia africana, which were exposed to either 10 °C or 20 °C, showed a clear difference in gene expression profiles between the west- and south coast ecotypes. This difference was particularly clear at 10 °C, where individuals from the south coast experienced a 'transcriptomic shock'. This low temperature is more typical of the west coast during upwelling events, and the physiological stress experienced by the south coast ecotype under such conditions may explain its absence from that coastline. Our results shed new light on the processes involved in driving genomic divergence and incipient speciation along coastlines with porous dispersal barriers.

UV-exposure, endogenous DNA damage, and DNA replication errors shape the spectra of genome changes in human skin

Human skin is continuously exposed to environmental DNA damage leading to the accumulation of somatic mutations over the lifetime of an individual. Mutagenesis in human skin cells can be also caused by endogenous DNA damage and by DNA replication errors. The contributions of these processes to the somatic mutation load in the skin of healthy humans has so far not been accurately assessed because the low numbers of mutations from current sequencing methodologies preclude the distinction between sequencing errors and true somatic genome changes. In this work, we sequenced genomes of single cell-derived clonal lineages obtained from primary skin cells of a large cohort of healthy individuals across a wide range of ages. We report here the range of mutation load and a comprehensive view of the various somatic genome changes that accumulate in skin cells. We demonstrate that UV-induced base substitutions, insertions and deletions are prominent even in sun-shielded skin. In addition, we detect accumulation of mutations due to spontaneous deamination of methylated cytosines as well as insertions and deletions characteristic of DNA replication errors in these cells. The endogenously induced somatic mutations and indels also demonstrate a linear increase with age, while UV-induced mutation load is age-independent. Finally, we show that DNA replication stalling at common fragile sites are potent sources of gross chromosomal rearrangements in human cells. Thus, somatic mutations in skin of healthy individuals reflect the interplay of environmental and endogenous factors in facilitating genome instability and carcinogenesis.

Skin forms the first barrier against a variety of environmental toxins and DNA damaging agents. Additionally, DNA of skin cells suffer from endogenous damage and errors during replication. Altogether, these lesions cause a variety of genome changes resulting in disease including cancer. However, the accurate measurement of the range and complete spectrum of genome changes in healthy skin was missing due to technical or biological limitations of prior studies. We present here accurate measurements of the various types of somatic genome changes that we found in skin fibroblasts and melanocytes from 21 donors ranging in ages from 25 to 79 years, which allowed to distinguish age related from age independent changes. Our cohort contains both White and African American donors, allowing an estimation of the impacts of skin color on mutagenesis. As a result, we revealed the complete spectrum and determined the range of somatic genome changes and their etiologies in healthy human skin fibroblasts and melanocytes and highlighted molecular mechanisms underlying these changes. Therefore, our study introduces a base line for defining disease levels of genome instability in skin.

Highly accurate long-read HiFi sequencing data for five complex genomes

The PacBio^® HiFi sequencing method yields highly accurate long-read sequencing datasets with read lengths averaging 10–25 kb and accuracies greater than 99.5%. These accurate long reads can be used to improve results for complex applications such as single nucleotide and structural variant detection, genome assembly, assembly of difficult polyploid or highly repetitive genomes, and assembly of metagenomes. Currently, there is a need for sample data sets to both evaluate the benefits of these long accurate reads as well as for development of bioinformatic tools including genome assemblers, variant callers, and haplotyping algorithms. We present deep coverage HiFi datasets for five complex samples including the two inbred model genomes Mus musculus and Zea mays, as well as two complex genomes, octoploid Fragaria × ananassa and the diploid anuran Rana muscosa. Additionally, we release sequence data from a mock metagenome community. The datasets reported here can be used without restriction to develop new algorithms and explore complex genome structure and evolution. Data were generated on the PacBio Sequel II System.

Measurement(s)	DNA • genome • Metagenome
Technology Type(s)	DNA sequencing • PacBio Sequel System
Factor Type(s)	organism that had its genome sequenced
Sample Characteristic - Organism	Mus musculus • Rana muscosa • Fragaria x ananassa • Zea mays

Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.12855527

Best practices for variant calling in clinical sequencing

Next-generation sequencing technologies have enabled a dramatic expansion of clinical genetic testing both for inherited conditions and diseases such as cancer. Accurate variant calling in NGS data is a critical step upon which virtually all downstream analysis and interpretation processes rely. Just as NGS technologies have evolved considerably over the past 10 years, so too have the software tools and approaches for detecting sequence variants in clinical samples. In this review, I discuss the current best practices for variant calling in clinical sequencing studies, with a particular emphasis on trio sequencing for inherited disorders and somatic mutation detection in cancer patients. I describe the relative strengths and weaknesses of panel, exome, and whole-genome sequencing for variant detection. Recommended tools and strategies for calling variants of different classes are also provided, along with guidance on variant review, validation, and benchmarking to ensure optimal performance. Although NGS technologies are continually evolving, and new capabilities (such as long-read single-molecule sequencing) are emerging, the “best practice” principles in this review should be relevant to clinical variant calling in the long term.

Next Generation Sequencing and Bioinformatics Analysis of Family Genetic Inheritance

Mendelian and complex genetic trait diseases continue to burden and affect society both socially and economically. The lack of effective tests has hampered diagnosis thus, the affected lack proper prognosis. Mendelian diseases are caused by genetic mutations in a singular gene while complex trait diseases are caused by the accumulation of mutations in either linked or unlinked genomic regions. Significant advances have been made in identifying novel diseases associated mutations especially with the introduction of next generation and third generation sequencing. Regardless, some diseases are still without diagnosis as most tests rely on SNP genotyping panels developed from population based genetic analyses. Analysis of family genetic inheritance using whole genomes, whole exomes or a panel of genes has been shown to be effective in identifying disease-causing mutations. In this review, we discuss next generation and third generation sequencing platforms, bioinformatic tools and genetic resources commonly used to analyze family based genomic data with a focus on identifying inherited or novel disease-causing mutations. Additionally, we also highlight the analytical, ethical and regulatory challenges associated with analyzing personal genomes which constitute the data used for family genetic inheritance.

Cathepsin D deficiency in mammary epithelium transiently stalls breast cancer by interference with mTORC1 signaling

Cathepsin D (CTSD) is a lysosomal protease and a marker of poor prognosis in breast cancer. However, the cells responsible for this association and the function of CTSD in cancer are still incompletely understood. By using a conditional CTSD knockout mouse crossed to the transgenic MMTV-PyMT breast cancer model we demonstrate that CTSD deficiency in the mammary epithelium, but not in myeloid cells, blocked tumor development in a cell-autonomous manner. We show that lack of CTSD impaired mechanistic Target of Rapamycin Complex 1 (mTORC1) signaling and induced reversible cellular quiescence. In line, CTSD-deficient tumors started to grow with a two-month delay and quiescent Ctsd^-/- tumor cells re-started proliferation upon long-term culture. This was accompanied by rewiring of oncogenic gene expression and signaling pathways, while mTORC1 signaling remained permanently disabled in CTSD-deficient cells. Together, these studies reveal a tumor cell-autonomous effect of CTSD deficiency, and establish a pivotal role of this protease in the cellular response to oncogenic stimuli.

The lysosomal aspartic protease Cathepsin D (CTSD) is associated with breast cancer progression. Here the authors show that selective inactivation of CTSD in mammary epithelium delays tumor onset due to impaired mTORC1 signaling, but resumes malignant growth due to compensatory oncogenic pathways

Rapid and economical drug resistance profiling with Nanopore MinION for clinical specimens with low bacillary burden of Mycobacterium tuberculosis

Objective

We designed and tested a Nanopore sequencing panel for direct tuberculosis drug resistance profiling. The panel targeted 10 resistance-associated loci. We assessed the feasibility of amplifying and sequencing these loci from 23 clinical specimens with low bacillary burden.

Results

At least 8 loci were successfully amplified from the majority for predicting first- and second-line drug resistance (14/23, 60.87%), and the 12 specimens yielding all 10 targets were sequenced with Nanopore MinION and Illumina MiSeq. MinION sequencing data was corrected by Nanopolish and recurrent variants were filtered. A total of 67,082 bases across all consensus sequences were analyzed, with 67,019 bases called by both MinION and MiSeq as wildtype. For the 41 single nucleotide variants (SNVs) called by MiSeq with 100% variant allelic frequency (VAF), 39 (95.1%) were called by MinION. For the 22 mixed bases called by MiSeq, a SNV with the highest VAF (70%) was called by MinION. With short assay time, reasonable reagent cost as well as continuously improving sequencing chemistry and signal correction pipelines, this Nanopore method can be a viable option for direct tuberculosis drug resistance profiling in the near future.