Mutational signatures of colorectal cancers according to distinct computational workflows.

Tumor mutational signatures have gained prominence in cancer research, yet the lack of standardized methods hinders reproducibility and robustness. Leveraging colorectal cancer (CRC) as a model, we explored the influence of computational parameters on mutational signature analyses across 230 CRC cell lines and 152 CRC patients. Results were validated in three independent datasets: 483 endometrial cancer patients stratified by mismatch repair (MMR) status, 35 lung cancer patients by smoking status and 12 patient-derived organoids (PDOs) annotated for colibactin exposure. Assessing various bioinformatic tools, reference datasets and input data sizes including whole genome sequencing, whole exome sequencing and a pan-cancer gene panel, we demonstrated significant variability in the results. We report that the use of distinct algorithms and references led to statistically different results, highlighting how arbitrary choices may induce variability in the mutational signature contributions. Furthermore, we found a differential contribution of mutational signatures between coding and intergenic regions and defined the minimum number of somatic variants required for reliable mutational signature assignment. To facilitate the identification of the most suitable workflows, we developed Comparative Mutational Signature analysis on Coding and Extragenic Regions (CoMSCER), a bioinformatic tool which allows researchers to easily perform comparative mutational signature analysis by coupling the results from several tools and public reference datasets and to assess mutational signature contributions in coding and non-coding genomic regions. In conclusion, our study provides a comparative framework to elucidate the impact of distinct computational workflows on mutational signatures.

Emergence of MET hyper-amplification at progression to MET and BRAF inhibition in colorectal cancer.

BACKGROUND: Combined MET and BRAF inhibition showed clinical benefit in a patient with rectal cancer carrying BRAFV600E and MET amplification. However after 4 months, acquired resistance emerged and the patient deceased shortly after disease progression. The mechanism of resistance to this drug combination is unknown. METHODS: We analysed plasma circulating tumour DNA obtained at progression by exome sequencing and digital PCR. MET gene and mRNA in situ hybridisation analyses in two bioptic specimens obtained at progression were used to confirm the plasma data. RESULTS: We identified in plasma MET gene hyper-amplification as a potential mechanism underlying therapy resistance. Increased MET gene copy and transcript levels were detected in liver and lymph node metastatic biopsies. Finally, transduction of MET in BRAF mutant colorectal cancer cells conferred refractoriness to BRAF and MET inhibition. CONCLUSIONS: We identified in a rectal cancer patient MET gene hyper-amplification as mechanism of resistance to dual BRAF and MET inhibition.

Small noncoding RNAs in cells transformed by human T-cell leukemia virus type 1: a role for a tRNA fragment as a primer for reverse transcriptase.

UNLABELLED: The present study employed mass sequencing of small RNA libraries to identify the repertoire of small noncoding RNAs expressed in normal CD4(+) T cells compared to cells transformed with human T-cell leukemia virus type 1 (HTLV-1), the causative agent of adult T-cell leukemia/lymphoma (ATLL). The results revealed distinct patterns of microRNA expression in HTLV-1-infected CD4(+) T-cell lines with respect to their normal counterparts. In addition, a search for virus-encoded microRNAs yielded 2 sequences that originated from the plus strand of the HTLV-1 genome. Several sequences derived from tRNAs were expressed at substantial levels in both uninfected and infected cells. One of the most abundant tRNA fragments (tRF-3019) was derived from the 3' end of tRNA-proline. tRF-3019 exhibited perfect sequence complementarity to the primer binding site of HTLV-1. The results of an in vitro reverse transcriptase assay verified that tRF-3019 was capable of priming HTLV-1 reverse transcriptase. Both tRNA-proline and tRF-3019 were detected in virus particles isolated from HTLV-1-infected cells. These findings suggest that tRF-3019 may play an important role in priming HTLV-1 reverse transcription and could thus represent a novel target to control HTLV-1 infection. IMPORTANCE: Small noncoding RNAs, a growing family of regulatory RNAs that includes microRNAs and tRNA fragments, have recently emerged as key players in many biological processes, including viral infection and cancer. In the present study, we employed mass sequencing to identify the repertoire of small noncoding RNAs in normal T cells compared to T cells transformed with human T-cell leukemia virus type 1 (HTLV-1), a retrovirus that causes adult T-cell leukemia/lymphoma. The results revealed a distinct pattern of microRNA expression in HTLV-1-infected cells and a tRNA fragment (tRF-3019) that was packaged into virions and capable of priming HTLV-1 reverse transcription, a key event in the retroviral life cycle. These findings indicate tRF-3019 could represent a novel target for therapies aimed at controlling HTLV-1 infection.

Specific inactivation of two immunomodulatory SIGLEC genes during human evolution.

Sialic acid-recognizing Ig-like lectins (Siglecs) are signaling receptors that modulate immune responses, and are targeted for interactions by certain pathogens. We describe two primate Siglecs that were rendered nonfunctional by single genetic events during hominin evolution after our common ancestor with the chimpanzee. SIGLEC13 was deleted by an Alu-mediated recombination event, and a single base pair deletion disrupted the ORF of SIGLEC17. Siglec-13 is expressed on chimpanzee monocytes, innate immune cells that react to bacteria. The human SIGLEC17P pseudogene mRNA is still expressed at high levels in human natural killer cells, which bridge innate and adaptive immune responses. As both resulting pseudogenes are homozygous in all human populations, we resurrected the originally encoded proteins and examined their functions. Chimpanzee Siglec-13 and the resurrected human Siglec-17 recruit a signaling adapter and bind sialic acids. Expression of either Siglec in innate immune cells alters inflammatory cytokine secretion in response to Toll-like receptor-4 stimulation. Both Siglecs can also be engaged by two potentially lethal sialylated bacterial pathogens of newborns and infants, agents with a potential impact on reproductive fitness. Neanderthal and Denisovan genomes show human-like sequences at both loci, corroborating estimates that the initial pseudogenization events occurred in the common ancestral population of these hominins. Both loci also show limited polymorphic diversity, suggesting selection forces predating the origin of modern humans. Taken together, these data suggest that genetic elimination of Siglec-13 and/or Siglec-17 represents signatures of infectious and/or other inflammatory selective processes contributing to population restrictions during hominin origins.

Comparative genomics revealed key molecular targets to rapidly convert a reference rifamycin-producing bacterial strain into an overproducer by genetic engineering.

Rifamycins are mainstay agents in treatment of many widespread diseases, but how an improved rifamycin producer can be created is still incompletely understood. Here, we describe a comparative genomic approach to investigate the mutational patterns introduced by the classical mutate-and-screen method in the genome of an improved rifamycin producer. Comparing the genome of the rifamycin B overproducer Amycolatopsis mediterranei HP-130 with those of the reference strains A. mediterranei S699 and U32, we identified 250 variations, affecting 227 coding sequences (CDS), 109 of which were HP-130-specific since they were absent in both S699 and U32. Mutational and transcriptional patterns indicated a series of genomic manipulations that not only proved the causative effect of mutB2 (coding for methylmalonyl-CoA mutase large subunit) and argS2 (coding for arginyl tRNA synthetase) mutations on the overproduction of rifamycin, but also constituted a rational strategy to genetically engineer a reference strain into an overproducer.

Preclinical efficacy of carfilzomib in BRAF-mutant colorectal cancer models.

Serine/threonine-protein kinase B-raf (BRAF) mutations are found in 8-15% of colorectal cancer patients and identify a subset of tumors with poor outcome in the metastatic setting. We have previously reported that BRAF-mutant human cells display a high rate of protein production, causing proteotoxic stress, and are selectively sensitive to the proteasome inhibitors bortezomib and carfilzomib. In this work, we tested whether carfilzomib could restrain the growth of BRAF-mutant colorectal tumors not only by targeting cancer cells directly, but also by promoting an immune-mediated antitumor response. In human and mouse colorectal cancer cells, carfilzomib triggered robust endoplasmic reticulum stress and autophagy, followed by the emission of immunogenic-damage-associated molecules. Intravenous administration of carfilzomib delayed the growth of BRAF-mutant murine tumors and mobilized the danger-signal proteins calreticulin and high mobility group box 1 (HMGB1). Analyses of drug-treated samples revealed increased intratumor recruitment of activated cytotoxic T cells and natural killers, concomitant with the downregulation of forkhead box protein P3 (Foxp3)+ T-cell surface glycoprotein CD4 (CD4)+ T cells, indicating that carfilzomib promotes reshaping of the immune microenvironment of BRAF-mutant murine colorectal tumors. These results will inform the design of clinical trials in BRAF-mutant colorectal cancer patients.

Transcriptome-wide gene expression outlier analysis pinpoints therapeutic vulnerabilities in colorectal cancer.

Multiple strategies are continuously being explored to expand the drug target repertoire in solid tumors. We devised a novel computational workflow for transcriptome-wide gene expression outlier analysis that allows the systematic identification of both overexpression and underexpression events in cancer cells. Here, it was applied to expression values obtained through RNA sequencing in 226 colorectal cancer (CRC) cell lines that were also characterized by whole-exome sequencing and microarray-based DNA methylation profiling. We found cell models displaying an abnormally high or low expression level for 3533 and 965 genes, respectively. Gene expression abnormalities that have been previously associated with clinically relevant features of CRC cell lines were confirmed. Moreover, by integrating multi-omics data, we identified both genetic and epigenetic alternations underlying outlier expression values. Importantly, our atlas of CRC gene expression outliers can guide the discovery of novel drug targets and biomarkers. As a proof of concept, we found that CRC cell lines lacking expression of the MTAP gene are sensitive to treatment with a PRMT5-MTA inhibitor (MRTX1719). Finally, other tumor types may also benefit from this approach.

Case report: Preclinical efficacy of NEDD8 and proteasome inhibitors in patient-derived models of signet ring high-grade mucinous colorectal cancer from a Lynch syndrome patient.

High-grade mucinous colorectal cancer (HGM CRC) is particularly aggressive, prone to metastasis and treatment resistance, frequently accompanied by "signet ring" cancer cells. A sizeable fraction of HGM CRCs (20-40%) arises in the context of the Lynch Syndrome, an autosomal hereditary syndrome that predisposes to microsatellite instable (MSI) CRC. Development of patient-derived preclinical models for this challenging subtype of colorectal cancer represents an unmet need in oncology. We describe here successful propagation of preclinical models from a case of early-onset, MSI-positive metastatic colorectal cancer in a male Lynch syndrome patient, refractory to standard care (FOLFOX6, FOLFIRI-Panitumumab) and, surprisingly, also to immunotherapy. Surgical material from a debulking operation was implanted in NOD/SCID mice, successfully yielding one patient-derived xenograft (PDX). PDX explants were subsequently used to generate 2D and 3D cell cultures. Histologically, all models resembled the tumor of origin, displaying a high-grade mucinous phenotype with signet ring cells. For preclinical exploration of alternative treatments, in light of recent findings, we considered inhibition of the proteasome by bortezomib and of the related NEDD8 pathway by pevonedistat. Indeed, sensitivity to bortezomib was observed in mucinous adenocarcinoma of the lung, and we previously found that HGM CRC is preferentially sensitive to pevonedistat in models with low or absent expression of cadherin 17 (CDH17), a differentiation marker. We therefore performed IHC on the tumor and models, and observed no CDH17 expression, suggesting sensitivity to pevonedistat. Both bortezomib and pevonedistat showed strong activity on 2D cells at 72 hours and on 3D organoids at 7 days, thus providing valid options for in vivo testing. Accordingly, three PDX cohorts were treated for four weeks, respectively with vehicle, bortezomib and pevonedistat. Both drugs significantly reduced tumor growth, as compared to the vehicle group. Interestingly, while bortezomib was more effective in vitro, pevonedistat was more effective in vivo. Drug efficacy was further substantiated by a reduction of cellularity and of Ki67-positive cells in the treated tumors. These results highlight proteasome and NEDD8 inhibition as potentially effective therapeutic approaches against Lynch syndrome-associated HGM CRC, also when the disease is refractory to all available treatment options.

Genetic and pharmacological modulation of DNA mismatch repair heterogeneous tumors promotes immune surveillance.

Patients affected by colorectal cancer (CRC) with DNA mismatch repair deficiency (MMRd), often respond to immune checkpoint blockade therapies, while those with mismatch repair-proficient (MMRp) tumors generally do not. Interestingly, a subset of MMRp CRCs contains variable fractions of MMRd cells, but it is unknown how their presence impacts immune surveillance. We asked whether modulation of the MMRd fraction in MMR heterogeneous tumors acts as an endogenous cancer vaccine by promoting immune surveillance. To test this hypothesis, we use isogenic MMRp (Mlh1+/+) and MMRd (Mlh1-/-) mouse CRC cells. MMRp/MMRd cells mixed at different ratios are injected in immunocompetent mice and tumor rejection is observed when at least 50% of cells are MMRd. To enrich the MMRd fraction, MMRp/MMRd tumors are treated with 6-thioguanine, which leads to tumor rejection. These results suggest that genetic and pharmacological modulation of the DNA mismatch repair machinery potentiate the immunogenicity of MMR heterogeneous tumors.

High-definition mapping of retroviral integration sites identifies active regulatory elements in human multipotent hematopoietic progenitors.

Integration of retroviral vectors in the human genome follows nonrandom patterns that favor insertional deregulation of gene expression and increase the risk of their use in clinical gene therapy. The molecular basis of retroviral target site selection is still poorly understood. We used deep sequencing technology to build genomewide, high-definition maps of > 60 000 integration sites of Moloney murine leukemia virus (MLV)- and HIV-based retroviral vectors in the genome of human CD34(+) multipotent hematopoietic progenitor cells (HPCs) and used gene expression profiling, chromatin immunoprecipitation, and bioinformatics to associate integration to genetic and epigenetic features of the HPC genome. Clusters of recurrent MLV integrations identify regulatory elements (alternative promoters, enhancers, evolutionarily conserved noncoding regions) within or around protein-coding genes and microRNAs with crucial functions in HPC growth and differentiation, bearing epigenetic marks of active or poised transcription (H3K4me1, H3K4me2, H3K4me3, H3K9Ac, Pol II) and specialized chromatin configurations (H2A.Z). Overall, we mapped 3500 high-frequency integration clusters, which represent a new resource for the identification of transcriptionally active regulatory elements. High-definition MLV integration maps provide a rational basis for predicting genotoxic risks in gene therapy and a new tool for genomewide identification of promoters and regulatory elements controlling hematopoietic stem and progenitor cell functions.

Comparative genomics and transcriptional profiles of Saccharopolyspora erythraea NRRL 2338 and a classically improved erythromycin over-producing strain.

BACKGROUND: The molecular mechanisms altered by the traditional mutation and screening approach during the improvement of antibiotic-producing microorganisms are still poorly understood although this information is essential to design rational strategies for industrial strain improvement. In this study, we applied comparative genomics to identify all genetic changes occurring during the development of an erythromycin overproducer obtained using the traditional mutate-and- screen method. RESULTS: Compared with the parental Saccharopolyspora erythraea NRRL 2338, the genome of the overproducing strain presents 117 deletion, 78 insertion and 12 transposition sites, with 71 insertion/deletion sites mapping within coding sequences (CDSs) and generating frame-shift mutations. Single nucleotide variations are present in 144 CDSs. Overall, the genomic variations affect 227 proteins of the overproducing strain and a considerable number of mutations alter genes of key enzymes in the central carbon and nitrogen metabolism and in the biosynthesis of secondary metabolites, resulting in the redirection of common precursors toward erythromycin biosynthesis. Interestingly, several mutations inactivate genes coding for proteins that play fundamental roles in basic transcription and translation machineries including the transcription anti-termination factor NusB and the transcription elongation factor Efp. These mutations, along with those affecting genes coding for pleiotropic or pathway-specific regulators, affect global expression profile as demonstrated by a comparative analysis of the parental and overproducer expression profiles. Genomic data, finally, suggest that the mutate-and-screen process might have been accelerated by mutations in DNA repair genes. CONCLUSIONS: This study helps to clarify the mechanisms underlying antibiotic overproduction providing valuable information about new possible molecular targets for rationale strain improvement.

Targeting the DNA Damage Response Pathways and Replication Stress in Colorectal Cancer.

PURPOSE: Genomic instability is a hallmark of cancer and targeting DNA damage response (DDR) is emerging as a promising therapeutic strategy in different solid tumors. The effectiveness of targeting DDR in colorectal cancer has not been extensively explored. EXPERIMENTAL DESIGN: We challenged 112 cell models recapitulating the genomic landscape of metastatic colorectal cancer with ATM, ATR, CHK1, WEE1, and DNA-PK inhibitors, in parallel with chemotherapeutic agents. We focused then on ATR inhibitors (ATRi) and, to identify putative biomarkers of response and resistance, we analyzed at multiple levels colorectal cancer models highly sensitive or resistant to these drugs. RESULTS: We found that around 30% of colorectal cancers, including those carrying KRAS and BRAF mutations and unresponsive to targeted agents, are sensitive to at least one DDR inhibitor. By investigating potential biomarkers of response to ATRi, we found that ATRi-sensitive cells displayed reduced phospho-RPA32 foci at basal level, while ATRi-resistant cells showed increased RAD51 foci formation in response to replication stress. Lack of ATM and RAD51C expression was associated with ATRi sensitivity. Analysis of mutational signatures and HRDetect score identified a subgroup of ATRi-sensitive models. Organoids derived from patients with metastatic colorectal cancer recapitulated findings obtained in cell lines. CONCLUSIONS: In conclusion, a subset of colorectal cancers refractory to current therapies could benefit from inhibitors of DDR pathways and replication stress. A composite biomarker involving phospho-RPA32 and RAD51 foci, lack of ATM and RAD51C expression, as well as analysis of mutational signatures could be used to identify colorectal cancers likely to respond to ATRi.

The complete mitochondrial genome of an 11,450-year-old aurochsen (Bos primigenius) from Central Italy.

BACKGROUND: Bos primigenius, the aurochs, is the wild ancestor of modern cattle breeds and was formerly widespread across Eurasia and northern Africa. After a progressive decline, the species became extinct in 1627. The origin of modern taurine breeds in Europe is debated. Archaeological and early genetic evidence point to a single Near Eastern origin and a subsequent spread during the diffusion of herding and farming. More recent genetic data are instead compatible with local domestication events or at least some level of local introgression from the aurochs. Here we present the analysis of the complete mitochondrial genome of a pre-Neolithic Italian aurochs. RESULTS: In this study, we applied a combined strategy employing both multiplex PCR amplifications and 454 pyrosequencing technology to sequence the complete mitochondrial genome of an 11,450-year-old aurochs specimen from Central Italy. Phylogenetic analysis of the aurochs mtDNA genome supports the conclusions from previous studies of short mtDNA fragments--namely that Italian aurochsen were genetically very similar to modern cattle breeds, but highly divergent from the North-Central European aurochsen. CONCLUSIONS: Complete mitochondrial genome sequences are now available for several modern cattle and two pre-Neolithic mtDNA genomes from very different geographic areas. These data suggest that previously identified sub-groups within the widespread modern cattle mitochondrial T clade are polyphyletic, and they support the hypothesis that modern European breeds have multiple geographic origins.

Complete mitochondrial genome sequence of the Tyrolean Iceman.

The Tyrolean Iceman was a witness to the Neolithic-Copper Age transition in Central Europe 5350-5100 years ago, and his mummified corpse was recovered from an Alpine glacier on the Austro-Italian border in 1991 [1]. Using a mixed sequencing procedure based on PCR amplification and 454 sequencing of pooled amplification products, we have retrieved the first complete mitochondrial-genome sequence of a prehistoric European. We have then compared it with 115 related extant lineages from mitochondrial haplogroup K. We found that the Iceman belonged to a branch of mitochondrial haplogroup K1 that has not yet been identified in modern European populations. This is the oldest complete Homo sapiens mtDNA genome generated to date. The results point to the potential significance of complete-ancient-mtDNA studies in addressing questions concerning the genetic history of human populations that the phylogeography of modern lineages is unable to tackle.

An Ariadne's thread to the identification and annotation of noncoding RNAs in eukaryotes.

Non-protein coding RNAs (ncRNAs) have emerged as a vast and heterogeneous portion of eukaryotic transcriptomes. Several ncRNA families, either short (<200 nucleotides, nt) or long (>200 nt), have been described and implicated in a variety of biological processes, from translation to gene expression regulation and nuclear trafficking. Most probably, other families are still to be discovered. Computational methods for ncRNA research require different approaches from the ones normally used in the prediction of protein-coding genes. Indeed, primary sequence alone is often insufficient to infer ncRNA functionality, whereas secondary structure and local conservation of portions of the transcript could provide useful information for both the prediction and the functional annotation of ncRNAs. Here we present an overview of computational methods and bioinformatics resources currently available for studying ncRNA genes, introducing the common themes as well as the different approaches required for long and short ncRNA identification and annotation.

Complete genome sequence of a serotype 11A, ST62 Streptococcus pneumoniae invasive isolate.

BACKGROUND: Streptococcus pneumoniae is an important human pathogen representing a major cause of morbidity and mortality worldwide. We sequenced the genome of a serotype 11A, ST62 S. pneumoniae invasive isolate (AP200), that was erythromycin-resistant due to the presence of the erm(TR) determinant, and carried out analysis of the genome organization and comparison with other pneumococcal genomes. RESULTS: The genome sequence of S. pneumoniae AP200 is 2,130,580 base pair in length. The genome carries 2216 coding sequences (CDS), 56 tRNA, and 12 rRNA genes. Of the CDSs, 72.9% have a predicted biological known function. AP200 contains the pilus islet 2 and, although its phenotype corresponds to serotype 11A, it contains an 11D capsular locus. Chromosomal rearrangements resulting from a large inversion across the replication axis, and horizontal gene transfer events were observed. The chromosomal inversion is likely implicated in the rebalance of the chromosomal architecture affected by the insertions of two large exogenous elements, the erm(TR)-carrying Tn1806 and a functional prophage designated φSpn_200. Tn1806 is 52,457 bp in size and comprises 49 ORFs. Comparative analysis of Tn1806 revealed the presence of a similar genetic element or part of it in related species such as Streptococcus pyogenes and also in the anaerobic species Finegoldia magna, Anaerococcus prevotii and Clostridium difficile. The genome of φSpn_200 is 35,989 bp in size and is organized in 47 ORFs grouped into five functional modules. Prophages similar to φSpn_200 were found in pneumococci and in other streptococcal species, showing a high degree of exchange of functional modules. φSpn_200 viral particles have morphologic characteristics typical of the Siphoviridae family and are capable of infecting a pneumococcal recipient strain. CONCLUSIONS: The sequence of S. pneumoniae AP200 chromosome revealed a dynamic genome, characterized by chromosomal rearrangements and horizontal gene transfers. The overall diversity of AP200 is driven mainly by the presence of the exogenous elements Tn1806 and φSpn_200 that show large gene exchanges with other genetic elements of different bacterial species. These genetic elements likely provide AP200 with additional genes, such as those conferring antibiotic-resistance, promoting its adaptation to the environment.

A Subset of Colorectal Cancers with Cross-Sensitivity to Olaparib and Oxaliplatin.

PURPOSE: Defects in the homologous recombination (HR) repair pathway are of clinical interest due to sensitivity of HR-deficient cells to PARP inhibitors. We were interested in defining PARP vulnerability in patients with metastatic colorectal cancer (mCRC) carrying KRAS and BRAF mutations who display poor prognosis, have limited therapeutic options, and represent an unmet clinical need. EXPERIMENTAL DESIGN: We tested colorectal cancer cell lines, patient-derived organoids (PDO), and patient-derived xenografts (PDX) enriched for KRAS and BRAF mutations for sensitivity to the PARP inhibitor olaparib, and the chemotherapeutic agents oxaliplatin and 5-fluorouracil (5-FU). Genomic profiles and DNA repair proficiency of colorectal cancer models were compared with pharmacologic response. RESULTS: Thirteen of 99 (around 13%) colorectal cancer cell lines were highly sensitive to clinically active concentrations of olaparib and displayed functional deficiency in HR. Response to PARP blockade was positively correlated with sensitivity to oxaliplatin in colorectal cancer cell lines as well as patient-derived organoids. Treatment of PDXs with olaparib impaired tumor growth and maintenance therapy with PARP blockade after initial oxaliplatin response delayed disease progression in mice. CONCLUSIONS: These results indicate that a colorectal cancer subset characterized by poor prognosis and limited therapeutic options is vulnerable to PARP inhibition and suggest that PDO-based drug-screening assays can be used to identify patients with colorectal cancer likely to benefit from olaparib. As patients with mCRC almost invariably receive therapies based on oxaliplatin, "maintenance" treatment with PARP inhibitors warrants further clinical investigation.

Phenotypes and gene expression profiles of Saccharopolyspora erythraea rifampicin-resistant (rif) mutants affected in erythromycin production.

BACKGROUND: There is evidence from previous works that bacterial secondary metabolism may be stimulated by genetic manipulation of RNA polymerase (RNAP). In this study we have used rifampicin selection as a strategy to genetically improve the erythromycin producer Saccharopolyspora erythraea. RESULTS: Spontaneous rifampicin-resistant (rif) mutants were isolated from the parental strain NRRL2338 and two rif mutations mapping within rpoB, S444F and Q426R, were characterized. With respect to the parental strain, S444F mutants exhibited higher respiratory performance and up to four-fold higher final erythromycin yields; in contrast, Q426R mutants were slow-growing, developmental-defective and severely impaired in erythromycin production. DNA microarray analysis demonstrated that these rif mutations deeply changed the transcriptional profile of S. erythraea. The expression of genes coding for key enzymes of carbon (and energy) and nitrogen central metabolism was dramatically altered in turn affecting the flux of metabolites through erythromycin feeder pathways. In particular, the valine catabolic pathway that supplies propionyl-CoA for biosynthesis of the erythromycin precursor 6-deoxyerythronolide B was strongly up-regulated in the S444F mutants, while the expression of the biosynthetic gene cluster of erythromycin (ery) was not significantly affected. In contrast, the ery cluster was down-regulated (<2-fold) in the Q426R mutants. These strains also exhibited an impressive stimulation of the nitrogen regulon, which may contribute to lower erythromycin yields as erythromycin production was strongly inhibited by ammonium. CONCLUSION: Rifampicin selection is a simple and reliable tool to investigate novel links between primary and secondary metabolism and morphological differentiation in S. erythraea and to improve erythromycin production. At the same time genome-wide analysis of expression profiles using DNA microarrays allowed information to be gained about the mechanisms underlying the stimulatory/inhibitory effects of the rif mutations on erythromycin production.

A Genomic Analysis Workflow for Colorectal Cancer Precision Oncology.

BACKGROUND: The diagnosis of colorectal cancer (CRC) is routinely accomplished through histopathologic examination. Prognostic information and treatment decisions are mainly determined by TNM classification, first defined in 1968. In the last decade, patient-specific CRC genomic landscapes were shown to provide important prognostic and predictive information. Therefore, there is a need for developing next generation sequencing (NGS) and bioinformatic workflows that can be routinely used for the assessment of prognostic and predictive biomarkers. MATERIALS AND METHODS: To foster the application of genomics in the clinical management of CRCs, the IDEA workflow has been built to easily adapt to the availability of patient specimens and the clinical question that is being asked. Initially, IDEA deploys ad-hoc NGS assays to interrogate predefined genomic target sequences (from 600 kb to 30 Mb) with optimal detection sensitivity. Next, sequencing data are processed through an integrated bioinformatic pipeline to assess single nucleotide variants, insertions and deletions, gene copy-number alterations, and chromosomal rearrangements. The overall results are gathered into a user-friendly report. RESULTS: We provide evidence that IDEA is capable of identifying clinically relevant molecular alterations. When optimized to analyze circulating tumor DNA, IDEA can be used to monitor response and relapse in the blood of patients with metastatic CRC receiving targeted agents. IDEA detected primary and secondary resistance mechanisms to ERBB2 blockade including sub-clonal RAS and BRAF mutations. CONCLUSIONS: The IDEA workflow provides a flexible platform to integrate NGS and bioinformatic tools for refined diagnosis and management of patients with advanced CRC.

Evolving neoantigen profiles in colorectal cancers with DNA repair defects.

BACKGROUND: Neoantigens that arise as a consequence of tumor-specific mutations can be recognized by T lymphocytes leading to effective immune surveillance. In colorectal cancer (CRC) and other tumor types, a high number of neoantigens is associated with patient response to immune therapies. The molecular processes governing the generation of neoantigens and their turnover in cancer cells are poorly understood. We exploited CRC as a model system to understand how alterations in DNA repair pathways modulate neoantigen profiles over time. METHODS: We performed whole exome sequencing (WES) and RNA sequencing (RNAseq) in CRC cell lines, in vitro and in vivo, and in CRC patient-derived xenografts (PDXs) to track longitudinally genomic profiles, clonal evolution, mutational signatures, and predicted neoantigens. RESULTS: The majority of CRC models showed remarkably stable mutational and neoantigen profiles; however, those carrying defects in DNA repair genes continuously diversified. Rapidly evolving and evolutionary stable CRCs displayed characteristic genomic signatures and transcriptional profiles. Downregulation of molecules implicated in antigen presentation occurred selectively in highly mutated and rapidly evolving CRC. CONCLUSIONS: These results indicate that CRCs carrying alterations in DNA repair pathways display dynamic neoantigen patterns that fluctuate over time. We define CRC subsets characterized by slow and fast evolvability and link this phenotype to downregulation of antigen-presenting cellular mechanisms. Longitudinal monitoring of the neoantigen landscape could be relevant in the context of precision medicine.