Use of high throughput sequencing to observe genome dynamics at a single cell level

Exploring the Role of Clustered Mutations in Carcinogenesis and Their Potential Clinical Implications in Cancer

Abnormal cell proliferation and growth leading to cancer primarily result from cumulative genome mutations. Single gene mutations alone do not fully explain cancer onset and progression; instead, clustered mutations-simultaneous occurrences of multiple mutations-are considered to be pivotal in cancer development and advancement. These mutations can affect different genes and pathways, resulting in cells undergoing malignant transformation with multiple functional abnormalities. Clustered mutations influence cancer growth rates, metastatic potential, and drug treatment sensitivity. This summary highlights the various types and characteristics of clustered mutations to understand their associations with carcinogenesis and discusses their potential clinical significance in cancer. As a unique mutation type, clustered mutations may involve genomic instability, DNA repair mechanism defects, and environmental exposures, potentially correlating with responsiveness to immunotherapy. Understanding the characteristics and underlying processes of clustered mutations enhances our comprehension of carcinogenesis and cancer progression, providing new diagnostic and therapeutic approaches for cancer.

A diet of oxidative stress-adapted bacteria improves stress resistance and lifespan in C. elegans via p38-MAPK

Organisms across taxa face stresses including variable temperature, redox imbalance, and xenobiotics. Successfully responding to stress and restoring homeostasis are crucial for survival. Aging is associated with a decreased stress response and alterations in the microbiome, which contribute to disease development. Animals and their microbiota share their environment; however, microbes have short generation time and can rapidly evolve and potentially affect host physiology during stress. Here, we leverage Caenorhabditis elegans and its simplified bacterial diet to demonstrate how microbial adaptation to oxidative stress affects the host's lifespan and stress response. We find that worms fed stress-evolved bacteria exhibit enhanced stress resistance and an extended lifespan. Through comprehensive genetic and metabolic analysis, we find that iron in stress-evolved bacteria enhances worm stress resistance and lifespan via activation of the mitogen-activated protein kinase pathway. In conclusion, our study provides evidence that understanding microbial stress-mediated adaptations could be used to slow aging and alleviate age-related health decline.

Identification and characterization of non-coding RNA networks in infected macrophages revealing the pathogenesis of F. nucleatum-associated diseases

BACKGROUND

F. nucleatum, as an important periodontal pathogen, is not only closely associated with the development of periodontitis, but also implicated in systemic diseases. Macrophages may act as an important mediator in the pathogenic process of F. nucleatum infection. As non-coding RNAs (ncRNAs) have attracted extensive attention as important epigenetic regulatory mechanisms recently, we focus on the competing endogenous RNA (ceRNA) regulatory networks to elucidate the pathogenesis of F. nucleatum-associated diseases.

RESULTS

We screen abnormally expressed mRNAs, miRNAs, lncRNAs and circRNAs in macrophages after F. nucleatum infection via the whole transcriptome sequencing technology, including 375 mRNAs, 5 miRNAs, 64 lncRNAs, and 180 circRNAs. The accuracy of RNA-seq and microRNA-seq result was further verified by qRT-PCR analysis. GO and KEGG analysis show that the differentially expressed genes were mainly involved in MAPK pathway, Toll-like receptor pathway, NF-κB pathway and apoptosis. KEGG disease analysis reveals that they were closely involved in immune system diseases, cardiovascular disease, cancers, inflammatory bowel disease (IBD) et al. We constructed the underlying lncRNA/circRNA-miRNA-mRNA networks to understand their interaction based on the correlation analysis between the differentially expressed RNAs, and then screen the core non-coding RNAs. In which, AKT2 is controlled by hsa_circ_0078617, hsa_circ_0069227, hsa_circ_0084089, lncRNA NUP210, lncRNA ABCB9, lncRNA DIXDC1, lncRNA ATXN1 and lncRNA XLOC_237387 through miR-150-5p; hsa_circ_0001165, hsa_circ_0008460, hsa_circ_0001118, lncRNA XLOC_237387 and lncRNA ATXN1 were identified as the ceRNAs of hsa-miR-146a-3p and thereby indirectly modulating the expression of MITF.

CONCLUSIONS

Our data identified promising candidate ncRNAs responsible for regulating immune response in the F. nucleatum-associated diseases, offering new insights regarding the pathogenic mechanism of this pathogen.

RNA-sequencing analysis of bisphenol A biodegradation by white-rot fungus Phanerochaete sordida YK-624

Bisphenol A (BPA) is a representative example of an endocrine-disrupting chemical. It is one of the most produced chemical substances in the world, but it causes harmful effects in organisms, such that the effective degradation of BPA is critical. The white-rot fungus Phanerochaete sordida YK-624 has been shown to effectively degrade BPA under ligninolytic and non-ligninolytic conditions. However, it is still unclear what kinds of enzymes are involved in BPA degradation. To explore the mechanism of BPA degradation, the present study analysed the functional genes of P. sordida YK-624 using RNA-sequencing (RNA-Seq). Oxidation-reduction process and metabolic pathway were enriched under ligninolytic and non-ligninolytic conditions by Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. It is suggested that BPA might be used as a carbon source by P. sordida YK-624. Lignin peroxidase and cytochrome P450 were detected in upregulated differentially expressed genes (DEGs). The lignin-degrading enzyme lignin peroxidase and the intracellular cytochrome P450 system were involved in BPA degradation by P. sordida YK-624, respectively. Furthermore, quantitative real-time PCR (qPCR) was used to validate the reliability of the RNA-Seq results.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03298-w.

Damage-Induced Mutation Clustering in Gram-Positive Bacteria: Preliminary Data

The phenomenon of a nonrandom distribution of mutations in a genome has been observed for many years. In fact, recent findings have indicated the presence of mutation clusters in different biological systems, including chemically treated yeast, transgenic mice, and human cancer cells. Until now, an asymmetrical distribution of mutations was only described in a single bacterial species. Here, we used ethyl methanesulfonate mutagenesis and a whole-genome sequencing approach to determine if this phenomenon is universal and not confined to Gram-negative bacteria. The Gram-positive bacterium Bacillus subtilis was selected for ethyl methanesulfonate treatment, followed by the next-generation sequencing of several mutagenized B. subtilis genomes. A nonrandom distribution of mutations was observed. This pilot study with a limited number of sequenced clones may indicate not only the universality of the phenomenon of mutation clusters but also the effectiveness of the use of a whole-genome sequencing approach in studying this phenomenon.

Pervasive lesion segregation shapes cancer genome evolution

Cancers arise through the acquisition of oncogenic mutations and grow by clonal expansion1,2. Here we reveal that most mutagenic DNA lesions are not resolved into a mutated DNA base pair within a single cell cycle. Instead, DNA lesions segregate, unrepaired, into daughter cells for multiple cell generations, resulting in the chromosome-scale phasing of subsequent mutations. We characterize this process in mutagen-induced mouse liver tumours and show that DNA replication across persisting lesions can produce multiple alternative alleles in successive cell divisions, thereby generating both multiallelic and combinatorial genetic diversity. The phasing of lesions enables accurate measurement of strand-biased repair processes, quantification of oncogenic selection and fine mapping of sister-chromatid-exchange events. Finally, we demonstrate that lesion segregation is a unifying property of exogenous mutagens, including UV light and chemotherapy agents in human cells and tumours, which has profound implications for the evolution and adaptation of cancer genomes.

Hypermutation in single-stranded DNA

Regions of genomic DNA can become single-stranded in the course of normal replication and transcription as well as during DNA repair. Abnormal repair and replication intermediates can contain large stretches of persistent single-stranded DNA, which is extremely vulnerable to DNA damaging agents and hypermutation. Since such single-stranded DNA spans only a fraction of the genome at a given instance, hypermutation in these regions leads to tightly-spaced mutation clusters. This phenomenon of hypermutation in single-stranded DNA has been documented in several experimental models as well as in cancer genomes. Recently, hypermutated single-stranded RNA viral genomes also have been documented. Moreover, indications of hypermutation in single-stranded DNA may also be found in the human germline. This review will summarize key current knowledge and the recent developments in understanding the diverse mechanisms and sources of ssDNA hypermutation.

FtsA G50E mutant suppresses the essential requirement for FtsK during bacterial cell division in Escherichia coli

In Escherichia coli, the N-terminal domain of the essential protein FtsK (FtsK_N) is proposed to modulate septum formation through the formation of dynamic and essential protein interactions with both the Z-ring and late-stage division machinery. Using genomic mutagenesis, complementation analysis, and in vitro pull-down assays, we aimed to identify protein interaction partners of FtsK essential to its function during division. Here, we identified the cytoplasmic Z-ring membrane anchoring protein FtsA as a direct protein-protein interaction partner of FtsK. Random genomic mutagenesis of an ftsK temperature-sensitive strain of E. coli revealed an FtsA point mutation (G50E) that is able to fully restore normal cell growth and morphology, and further targeted site-directed mutagenesis of FtsA revealed several other point mutations capable of fully suppressing the essential requirement for functional FtsK. Together, this provides insight into a potential novel co-complex formed between these components during division and suggests FtsA may directly impact FtsK function.

Topokaryotyping demonstrates single cell variability and stress dependent variations in nuclear envelope associated domains

Analysis of large-scale interphase genome positioning with reference to a nuclear landmark has recently been studied using sequencing-based single cell approaches. However, these approaches are dependent upon technically challenging, time consuming and costly high throughput sequencing technologies, requiring specialized bioinformatics tools and expertise. Here, we propose a novel, affordable and robust microscopy-based single cell approach, termed Topokaryotyping, to analyze and reconstruct the interphase positioning of genomic loci relative to a given nuclear landmark, detectable as banding pattern on mitotic chromosomes. This is accomplished by proximity-dependent histone labeling, where biotin ligase BirA fused to nuclear envelope marker Emerin was coexpressed together with Biotin Acceptor Peptide (BAP)-histone fusion followed by (i) biotin labeling, (ii) generation of mitotic spreads, (iii) detection of the biotin label on mitotic chromosomes and (iv) their identification by karyotyping. Using Topokaryotyping, we identified both cooperativity and stochasticity in the positioning of emerin-associated chromatin domains in individual cells. Furthermore, the chromosome-banding pattern showed dynamic changes in emerin-associated domains upon physical and radiological stress. In summary, Topokaryotyping is a sensitive and reliable technique to quantitatively analyze spatial positioning of genomic regions interacting with a given nuclear landmark at the single cell level in various experimental conditions.

Real-time dynamics of mutagenesis reveal the chronology of DNA repair and damage tolerance responses in single cells

Evolutionary processes are driven by diverse molecular mechanisms that act in the creation and prevention of mutations. It remains unclear how these mechanisms are regulated because limitations of existing mutation assays have precluded measuring how mutation rates vary over time in single cells. Toward this goal, I detected nascent DNA mismatches as a proxy for mutagenesis and simultaneously followed gene expression dynamics in single Escherichia coli cells using microfluidics. This general microscopy-based approach revealed the real-time dynamics of mutagenesis in response to DNA alkylation damage and antibiotic treatments. It also enabled relating the creation of DNA mismatches to the chronology of the underlying molecular processes. By avoiding population averaging, I discovered cell-to-cell variation in mutagenesis that correlated with heterogeneity in the expression of alternative responses to DNA damage. Pulses of mutagenesis are shown to arise from transient DNA repair deficiency. Constitutive expression of DNA repair pathways and induction of damage tolerance by the SOS response compensate for delays in the activation of inducible DNA repair mechanisms, together providing robustness against the toxic and mutagenic effects of DNA alkylation damage.

A streptomycin resistance marker in H. parasuis based on site-directed mutations in rpsL gene to perform unmarked in-frame mutations and to verify natural transformation

Haemophilus parasuis is a member of the family Pasteurellaceae and a major causative agent of Glässer’s disease. This bacterium is normally a benign swine commensal but may become a deadly pathogen upon penetration into multiple tissues, contributing to severe lesions in swine. We have established a successive natural transformation-based markerless mutation system in this species. However, the two-step mutation system requires screening of natural competent cells, and cannot delete genes which regulate natural competence per se. In this study, we successfully obtained streptomycin-resistant derivatives from H. parasuis wild type strain SC1401 by using ethyl methane sulfonate (EMS, CH₃SO₂OC₂H₅). Upon sequencing and site-directed mutations, we uncovered that the EMS-induced point mutation in rpsL at codon 43rd (AAA → AGA; K43R) or at 88th (AAA → AGA; K88R) confers a much higher streptomycin resistance than clinical isolates. We have applied the streptomycin resistance marker as a positive selection marker to perform homologous recombination through conjugation and successfully generated a double unmarked in-frame targeted mutant 1401D88△tfox△arcA. Combined with a natural transformation-based knockout system and this genetic technique, multiple deletion mutants or attenuated strains of H. parasuis can be easily constructed. Moreover, the mutant genetic marker rpsL and streptomycin resistant phenotypes can serve as an effective tool to select naturally competent strains, and to verify natural transformation quantitatively.

Annotation and characterization of Cd-responsive metal transporter genes in rapeseed (Brassica napus)

In higher plants, heavy metal transporters are responsible for metal uptake, translocation and homeostasis. These metals include essential metals such as zinc (Zn) or manganese (Mn) and non-essential metals like cadmium (Cd) or lead (Pb). Although a few heavy metal transporters have been well identified in model plants (e.g. Arabidopsis and rice), little is known about their functionality in rapeseed (Brassica napus). B. napus is an important oil crop ranking the third largest sources of vegetable oil over the world. Importantly, B. napus has long been considered as a desirable candidate for phytoremediation owning to its massive dry weight productivity and moderate to high Cd accumulation. In this study, 270 metal transporter genes (MTGs) from B. napus genome were identified and annotated using bioinformatics and high-throughput sequencing. Most of the MTGs (74.8%, 202/270) were validated by RNA-sequencing (RNA-seq) the seedling libraries. Based on the sequence identity, nine superfamilies including YSL, OPT, NRAMP, COPT, ZIP, CDF/MTP, HMA, MRP and PDR have been classified. RNA-sequencing profiled 202 non-redundant MTGs from B. napus seedlings, of which, 108 MTGs were differentially expressed and 62 genes were significantly induced under Cd stress. These differentially expressed genes (DEGs) are dispersed in the rapeseed genome. Some of the genes were well confirmed by qRT-PCR. Analysis of the genomic distribution of MTGs on B. napus chromosomes revealed that their evolutional expansion was probably through localized allele duplications.

Clusters of Multiple Mutations: Incidence and Molecular Mechanisms

It has been long understood that mutation distribution is not completely random across genomic space and in time. Indeed, recent surprising discoveries identified multiple simultaneous mutations occurring in tiny regions within chromosomes while the rest of the genome remains relatively mutation-free. Mechanistic elucidation of these phenomena, called mutation showers, mutation clusters, or kataegis, in parallel with findings of abundant clustered mutagenesis in cancer genomes, is ongoing. So far, the combination of factors most important for clustered mutagenesis is the induction of DNA lesions within unusually long and persistent single-strand DNA intermediates. In addition to being a fascinating phenomenon, clustered mutagenesis also became an indispensable tool for identifying a previously unrecognized major source of mutation in cancer, APOBEC cytidine deaminases. Future research on clustered mutagenesis may shed light onto important mechanistic details of genome maintenance, with potentially profound implications for human health.

RNA-seq based transcriptomic analysis of single bacterial cells

Gene-expression heterogeneity among individual cells determines the fate of a bacterial population. Here we report the first bacterial single-cell RNA sequencing (RNA-seq), BaSiC RNA-seq, a method integrating RNA isolation, cDNA synthesis and amplification, and RNA-seq analysis of the whole transcriptome of single cyanobacterium Synechocystis sp. PCC 6803 cells which typically contain approximately 5-7 femtogram total RNA per cell. We applied the method to 3 Synechocystis single cells at 24 h and 3 single cells at 72 h after nitrogen-starvation stress treatment, as well as their bulk-cell controls under the same conditions, to determine the heterogeneity upon environmental stress. With 82-98% and 31-48% of all putative Synechocystis genes identified in single cells of 24 and 72 h, respectively, the results demonstrated that the method could achieve good identification of the transcripts in single bacterial cells. In addition, the preliminary results from nitrogen-starved cells also showed a possible increasing gene-expression heterogeneity from 24 h to 72 h after nitrogen starvation stress. Moreover, preliminary analysis of single-cell transcriptomic datasets revealed that genes from the "Mobile elements" functional category have the most significant increase of gene-expression heterogeneity upon stress, which was further confirmed by single-cell RT-qPCR analysis of gene expression in 24 randomly selected cells.

A comprehensive survey of the mutagenic impact of common cancer cytotoxics

BACKGROUND

Genomic mutations caused by cytotoxic agents used in cancer chemotherapy may cause secondary malignancies as well as contribute to the evolution of treatment-resistant tumour cells. The stable diploid genome of the chicken DT40 lymphoblast cell line, an established DNA repair model system, is well suited to accurately assay genomic mutations.

RESULTS

We use whole genome sequencing of multiple DT40 clones to determine the mutagenic effect of eight common cytotoxics used for the treatment of millions of patients worldwide. We determine the spontaneous mutagenesis rate at 2.3 × 10(-10) per base per cell division and find that cisplatin, cyclophosphamide and etoposide induce extra base substitutions with distinct spectra. After four cycles of exposure, cisplatin induces 0.8 mutations per Mb, equivalent to the median mutational burden in common leukaemias. Cisplatin-induced mutations, including short insertions and deletions, are mainly located at sites of putative intrastrand crosslinks. We find two of the newly defined cisplatin-specific mutation types as causes of the reversion of BRCA2 mutations in emerging cisplatin-resistant tumours or cell clones. Gemcitabine, 5-fluorouracil, hydroxyurea, doxorubicin and paclitaxel have no measurable mutagenic effect. The cisplatin-induced mutation spectrum shows good correlation with cancer mutation signatures attributed to smoking and other sources of guanine-directed base damage.

CONCLUSION

This study provides support for the use of cell line mutagenesis assays to validate or predict the mutagenic effect of environmental and iatrogenic exposures. Our results suggest genetic reversion due to cisplatin-induced mutations as a distinct mechanism for developing resistance.

New technologies for DNA analysis--a review of the READNA Project

The REvolutionary Approaches and Devices for Nucleic Acid analysis (READNA) project received funding from the European Commission for 41/2 years. The objectives of the project revolved around technological developments in nucleic acid analysis. The project partners have discovered, created and developed a huge body of insights into nucleic acid analysis, ranging from improvements and implementation of current technologies to the most promising sequencing technologies that constitute a 3(rd) and 4(th) generation of sequencing methods with nanopores and in situ sequencing, respectively.

Genomic approaches to DNA repair and mutagenesis

DNA damage is a constant threat to cells, causing cytotoxicity as well as inducing genetic alterations. The steady-state abundance of DNA lesions in a cell is minimized by a variety of DNA repair mechanisms, including DNA strand break repair, mismatch repair, nucleotide excision repair, base excision repair, and ribonucleotide excision repair. The efficiencies and mechanisms by which these pathways remove damage from chromosomes have been primarily characterized by investigating the processing of lesions at defined genomic loci, among bulk genomic DNA, on episomal DNA constructs, or using in vitro substrates. However, the structure of a chromosome is heterogeneous, consisting of heavily protein-bound heterochromatic regions, open regulatory regions, actively transcribed genes, and even areas of transient single stranded DNA. Consequently, DNA repair pathways function in a much more diverse set of chromosomal contexts than can be readily assessed using previous methods. Recent efforts to develop whole genome maps of DNA damage, repair processes, and even mutations promise to greatly expand our understanding of DNA repair and mutagenesis. Here we review the current efforts to utilize whole genome maps of DNA damage and mutation to understand how different chromosomal contexts affect DNA excision repair pathways.

Identifying structural variation in haploid microbial genomes from short-read resequencing data using breseq

Background

Mutations that alter chromosomal structure play critical roles in evolution and disease, including in the origin of new lifestyles and pathogenic traits in microbes. Large-scale rearrangements in genomes are often mediated by recombination events involving new or existing copies of mobile genetic elements, recently duplicated genes, or other repetitive sequences. Most current software programs for predicting structural variation from short-read DNA resequencing data are intended primarily for use on human genomes. They typically disregard information in reads mapping to repeat sequences, and significant post-processing and manual examination of their output is often required to rule out false-positive predictions and precisely describe mutational events.

Results

We have implemented an algorithm for identifying structural variation from DNA resequencing data as part of the breseq computational pipeline for predicting mutations in haploid microbial genomes. Our method evaluates the support for new sequence junctions present in a clonal sample from split-read alignments to a reference genome, including matches to repeat sequences. Then, it uses a statistical model of read coverage evenness to accept or reject these predictions. Finally, breseq combines predictions of new junctions and deleted chromosomal regions to output biologically relevant descriptions of mutations and their effects on genes. We demonstrate the performance of breseq on simulated Escherichia coli genomes with deletions generating unique breakpoint sequences, new insertions of mobile genetic elements, and deletions mediated by mobile elements. Then, we reanalyze data from an E. coli K-12 mutation accumulation evolution experiment in which structural variation was not previously identified. Transposon insertions and large-scale chromosomal changes detected by breseq account for ~25% of spontaneous mutations in this strain. In all cases, we find that breseq is able to reliably predict structural variation with modest read-depth coverage of the reference genome (>40-fold).

Conclusions

Using breseq to predict structural variation should be useful for studies of microbial epidemiology, experimental evolution, synthetic biology, and genetics when a reference genome for a closely related strain is available. In these cases, breseq can discover mutations that may be responsible for important or unintended changes in genomes that might otherwise go undetected.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1039) contains supplementary material, which is available to authorized users.

Benchmarking undedicated cloud computing providers for analysis of genomic datasets

A major bottleneck in biological discovery is now emerging at the computational level. Cloud computing offers a dynamic means whereby small and medium-sized laboratories can rapidly adjust their computational capacity. We benchmarked two established cloud computing services, Amazon Web Services Elastic MapReduce (EMR) on Amazon EC2 instances and Google Compute Engine (GCE), using publicly available genomic datasets (E.coli CC102 strain and a Han Chinese male genome) and a standard bioinformatic pipeline on a Hadoop-based platform. Wall-clock time for complete assembly differed by 52.9% (95% CI: 27.5-78.2) for E.coli and 53.5% (95% CI: 34.4-72.6) for human genome, with GCE being more efficient than EMR. The cost of running this experiment on EMR and GCE differed significantly, with the costs on EMR being 257.3% (95% CI: 211.5-303.1) and 173.9% (95% CI: 134.6-213.1) more expensive for E.coli and human assemblies respectively. Thus, GCE was found to outperform EMR both in terms of cost and wall-clock time. Our findings confirm that cloud computing is an efficient and potentially cost-effective alternative for analysis of large genomic datasets. In addition to releasing our cost-effectiveness comparison, we present available ready-to-use scripts for establishing Hadoop instances with Ganglia monitoring on EC2 or GCE.

A new workflow for whole-genome sequencing of single human cells

Unbiased amplification of the whole-genome amplification (WGA) of single cells is crucial to study cancer evolution and genetic heterogeneity, but is challenging due to the high complexity of the human genome. Here, we present a new workflow combining an efficient adapter-linker PCR-based WGA method with second-generation sequencing. This approach allows comparison of single cells at base pair resolution. Amplification recovered up to 74% of the human genome. Copy-number variants and loss of heterozygosity detected in single cell genomes showed concordance of up to 99% to pooled genomic DNA. Allele frequencies of mutations could be determined accurately due to an allele dropout rate of only 2%, clearly demonstrating the low bias of our PCR-based WGA approach. Sequencing with paired-end reads allowed genome-wide analysis of structural variants. By direct comparison to other WGA methods, we further endorse its suitability to analyze genetic heterogeneity.

Quantum biology at the cellular level--elements of the research program

Quantum biology is emerging as a new field at the intersection between fundamental physics and biology, promising novel insights into the nature and origin of biological order. We discuss several elements of QBCL (quantum biology at cellular level) - a research program designed to extend the reach of quantum concepts to higher than molecular levels of biological organization. We propose a new general way to address the issue of environmentally induced decoherence and macroscopic superpositions in biological systems, emphasizing the 'basis-dependent' nature of these concepts. We introduce the notion of 'formal superposition' and distinguish it from that of Schroedinger's cat (i.e., a superposition of macroscopically distinct states). Whereas the latter notion presents a genuine foundational problem, the former one contradicts neither common sense nor observation, and may be used to describe cellular 'decision-making' and adaptation. We stress that the interpretation of the notion of 'formal superposition' should involve non-classical correlations between molecular events in a cell. Further, we describe how better understanding of the physics of Life can shed new light on the mechanism driving evolutionary adaptation (viz., 'Basis-Dependent Selection', BDS). Experimental tests of BDS and the potential role of synthetic biology in closing the 'evolvability mechanism' loophole are also discussed.

Direct DNA Lesion Reversal and Excision Repair in Escherichia coli

Cellular DNA is constantly challenged by various endogenous and exogenous genotoxic factors that inevitably lead to DNA damage: structural and chemical modifications of primary DNA sequence. These DNA lesions are either cytotoxic, because they block DNA replication and transcription, or mutagenic due to the miscoding nature of the DNA modifications, or both, and are believed to contribute to cell lethality and mutagenesis. Studies on DNA repair in Escherichia coli spearheaded formulation of principal strategies to counteract DNA damage and mutagenesis, such as: direct lesion reversal, DNA excision repair, mismatch and recombinational repair and genotoxic stress signalling pathways. These DNA repair pathways are universal among cellular organisms. Mechanistic principles used for each repair strategies are fundamentally different. Direct lesion reversal removes DNA damage without need for excision and de novo DNA synthesis, whereas DNA excision repair that includes pathways such as base excision, nucleotide excision, alternative excision and mismatch repair, proceeds through phosphodiester bond breakage, de novo DNA synthesis and ligation. Cell signalling systems, such as adaptive and oxidative stress responses, although not DNA repair pathways per se, are nevertheless essential to counteract DNA damage and mutagenesis. The present review focuses on the nature of DNA damage, direct lesion reversal, DNA excision repair pathways and adaptive and oxidative stress responses in E. coli.

Genotyping in the cloud with Crossbow

Crossbow is a scalable, portable, and automatic cloud computing tool for identifying SNPs from high-coverage, short-read resequencing data. It is built on Apache Hadoop, an implementation of the MapReduce software framework. Hadoop allows Crossbow to distribute read alignment and SNP calling subtasks over a cluster of commodity computers. Two robust tools, Bowtie and SOAPsnp, implement the fundamental alignment and variant calling operations respectively, and have demonstrated capabilities within Crossbow of analyzing approximately one billion short reads per hour on a commodity Hadoop cluster with 320 cores. Through protocol examples, this unit will demonstrate the use of Crossbow for identifying variations in three different operating modes: on a Hadoop cluster, on a single computer, and on the Amazon Elastic MapReduce cloud computing service.

Developing insights into the mechanisms of evolution of bacterial pathogens from whole-genome sequences

Evolution of bacterial pathogen populations has been detected in a variety of ways including phenotypic tests, such as metabolic activity, reaction to antisera and drug resistance and genotypic tests that measure variation in chromosome structure, repetitive loci and individual gene sequences. While informative, these methods only capture a small subset of the total variation and, therefore, have limited resolution. Advances in sequencing technologies have made it feasible to capture whole-genome sequence variation for each sample under study, providing the potential to detect all changes at all positions in the genome from single nucleotide changes to large-scale insertions and deletions. In this review, we focus on recent work that has applied this powerful new approach and summarize some of the advances that this has brought in our understanding of the details of how bacterial pathogens evolve.

Two mechanisms produce mutation hotspots at DNA breaks in Escherichia coli

Mutation hotspots and showers occur across phylogeny and profoundly influence genome evolution, yet the mechanisms that produce hotspots remain obscure. We report that DNA double-strand breaks (DSBs) provoke mutation hotspots via stress-induced mutation in Escherichia coli. With tet reporters placed 2 kb to 2 Mb (half the genome) away from an I-SceI site, RpoS/DinB-dependent mutations occur maximally within the first 2 kb and decrease logarithmically to ∼60 kb. A weak mutation tail extends to 1 Mb. Hotspotting occurs independently of I-site/tet-reporter-pair position in the genome, upstream and downstream in the replication path. RecD, which allows RecBCD DSB-exonuclease activity, is required for strong local but not long-distance hotspotting, indicating that double-strand resection and gap-filling synthesis underlie local hotspotting, and newly illuminating DSB resection in vivo. Hotspotting near DSBs opens the possibility that specific genomic regions could be targeted for mutagenesis, and could also promote concerted evolution (coincident mutations) within genes/gene clusters, an important issue in the evolution of protein functions.

A single nucleotide polymorphism in ycdC alters tRNA synthetase expression and results in hypersecretion in Escherichia coli

The most important approach to the development of platform organisms for recombinant protein production relies on random mutagenesis and phenotypic selection. Complex phenotypes, including those associated with significantly elevated expression and secretion of heterologous proteins, are the result of multiple genomic mutations. Using next generation sequencing, a parent and derivative hypersecreter strain (B41) of Escherichia coli were sequenced with an average coverage of 52.8X and 55X, respectively. A new base-pair calling program, revealed a single nucleotide polymorphism in the B41 genome at position 1,074,787, resulting in translation termination near the N-terminus of a transcriptional regulator protein, RutR, coded by the ycdC gene. We verified the hypersecretion phenotype in a ycdC::Tn5 mutant and observed a 3.4-fold increase in active hemolysin secretion, consistent with the increase observed in B41 strain. mRNA expression profiling showed decreased expression of tRNA-synthetases and some amino acid transporters in the ycdC::Tn5 mutant. This study demonstrates the power of next generation sequencing to characterize mutants leading to successful metabolic engineering strategies for strain improvement.

Escherichia coli frameshift mutation rate depends on the chromosomal context but not on the GATC content near the mutation site

Different studies have suggested that mutation rate varies at different positions in the genome. In this work we analyzed if the chromosomal context and/or the presence of GATC sites can affect the frameshift mutation rate in the Escherichia coli genome. We show that in a mismatch repair deficient background, a condition where the mutation rate reflects the fidelity of the DNA polymerization process, the frameshift mutation rate could vary up to four times among different chromosomal contexts. Furthermore, the mismatch repair efficiency could vary up to eight times when compared at different chromosomal locations, indicating that detection and/or repair of frameshift events also depends on the chromosomal context. Also, GATC sequences have been proved to be essential for the correct functioning of the E. coli mismatch repair system. Using bacteriophage heteroduplexes molecules it has been shown that GATC influence the mismatch repair efficiency in a distance- and number-dependent manner, being almost nonfunctional when GATC sequences are located at 1 kb or more from the mutation site. Interestingly, we found that in E. coli genomic DNA the mismatch repair system can efficiently function even if the nearest GATC sequence is located more than 2 kb away from the mutation site. The results presented in this work show that even though frameshift mutations can be efficiently generated and/or repaired anywhere in the genome, these processes can be modulated by the chromosomal context that surrounds the mutation site.

Cumulative number of cell divisions as a meaningful timescale for adaptive laboratory evolution of Escherichia coli

Adaptive laboratory evolution (ALE) under controlled conditions has become a valuable approach for the study of the genetic and biochemical basis for microbial adaptation under a given selection pressure. Conventionally, the timescale in ALE experiments has been set in terms of number of generations. As mutations are believed to occur primarily during cell division in growing cultures, the cumulative number of cell divisions (CCD) would be an alternative way to set the timescale for ALE. Here we show that in short-term ALE (up to 40–50 days), Escherichia coli, under growth rate selection pressure, was found to undergo approximately 10^11.2 total cumulative cell divisions in the population to produce a new stable growth phenotype that results from 2 to 8 mutations. Continuous exposure to a low level of the mutagen N-methyl-N′-nitro-N-nitrosoguanidine was found to accelerate this timescale and led to a superior growth rate phenotype with a much larger number of mutations as determined with whole-genome sequencing. These results would be useful for the fundamental kinetics of the ALE process in designing ALE experiments and provide a basis for its quantitative description.

Damage-induced localized hypermutability

Genome instability continuously presents perils of cancer, genetic disease and death of a cell or an organism. At the same time, it provides for genome plasticity that is essential for development and evolution. We address here the genome instability confined to a small fraction of DNA adjacent to free DNA ends at uncapped telomeres and double-strand breaks. We found that budding yeast cells can tolerate nearly 20 kilobase regions of subtelomeric single-strand DNA that contain multiple UV-damaged nucleotides. During restoration to the double-strand state, multiple mutations are generated by error-prone translesion synthesis. Genome-wide sequencing demonstrated that multiple regions of damage-induced localized hypermutability can be tolerated, which leads to the simultaneous appearance of multiple mutation clusters in the genomes of UV- irradiated cells. High multiplicity and density of mutations suggest that this novel form of genome instability may play significant roles in generating new alleles for evolutionary selection as well as in the incidence of cancer and genetic disease.

Most "dark matter" transcripts are associated with known genes

Short-read RNA sequencing in mouse and human tissues shows that most transcripts are encoded within or nearby known genes and that most of the genome is not transcribed.

A series of reports over the last few years have indicated that a much larger portion of the mammalian genome is transcribed than can be accounted for by currently annotated genes, but the quantity and nature of these additional transcripts remains unclear. Here, we have used data from single- and paired-end RNA-Seq and tiling arrays to assess the quantity and composition of transcripts in PolyA+ RNA from human and mouse tissues. Relative to tiling arrays, RNA-Seq identifies many fewer transcribed regions (“seqfrags”) outside known exons and ncRNAs. Most nonexonic seqfrags are in introns, raising the possibility that they are fragments of pre-mRNAs. The chromosomal locations of the majority of intergenic seqfrags in RNA-Seq data are near known genes, consistent with alternative cleavage and polyadenylation site usage, promoter- and terminator-associated transcripts, or new alternative exons; indeed, reads that bridge splice sites identified 4,544 new exons, affecting 3,554 genes. Most of the remaining seqfrags correspond to either single reads that display characteristics of random sampling from a low-level background or several thousand small transcripts (median length = 111 bp) present at higher levels, which also tend to display sequence conservation and originate from regions with open chromatin. We conclude that, while there are bona fide new intergenic transcripts, their number and abundance is generally low in comparison to known exons, and the genome is not as pervasively transcribed as previously reported.

The human genome was sequenced a decade ago, but its exact gene composition remains a subject of debate. The number of protein-coding genes is much lower than initially expected, and the number of distinct transcripts is much larger than the number of protein-coding genes. Moreover, the proportion of the genome that is transcribed in any given cell type remains an open question: results from “tiling” microarray analyses suggest that transcription is pervasive and that most of the genome is transcribed, whereas new deep sequencing-based methods suggest that most transcripts originate from known genes. We have addressed this discrepancy by comparing samples from the same tissues using both technologies. Our analyses indicate that RNA sequencing appears more reliable for transcripts with low expression levels, that most transcripts correspond to known genes or are near known genes, and that many transcripts may represent new exons or aberrant products of the transcription process. We also identify several thousand small transcripts that map outside known genes; their sequences are often conserved and are often encoded in regions of open chromatin. We propose that most of these transcripts may be by-products of the activity of enhancers, which associate with promoters as part of their role as long-range gene regulatory sites. Overall, however, we find that most of the genome is not appreciably transcribed.