Theoretical analysis of mutation hotspots and their DNA sequence context specificity

Insights into the ANKRD11 variants and short-stature phenotype through literature review and ClinVar database search

Ankyrin repeat domain containing-protein 11 (ANKRD11), a transcriptional factor predominantly localized in the cell nucleus, plays a crucial role in the expression regulation of key genes by recruiting chromatin remodelers and interacting with specific transcriptional repressors or activators during numerous biological processes. Its pathogenic variants are strongly linked to the pathogenesis and progression of multisystem disorder known as KBG syndrome. With the widespread application of high-throughput DNA sequencing technologies in clinical medicine, numerous pathogenic variants in the ANKRD11 gene have been reported. Patients with KBG syndrome usually exhibit a broad phenotypic spectrum with a variable degree of severity, even if having identical variants. In addition to distinctive dental, craniofacial and neurodevelopmental abnormalities, patients often present with skeletal anomalies, particularly postnatal short stature. The relationship between ANKRD11 variants and short stature is not well-understood, with limited knowledge regarding its occurrence rate or underlying biological mechanism involved. This review aims to provide an updated analysis of the molecular spectrum associated with ANKRD11 variants, investigate the prevalence of the short stature among patients harboring these variants, evaluate the efficacy of recombinant human growth hormone in treating children with short stature and ANKRD11 variants, and explore the biological mechanisms underlying short stature from both scientific and clinical perspectives. Our investigation indicated that frameshift and nonsense were the most frequent types in 583 pathogenic or likely pathogenic variants identified in the ANKRD11 gene. Among the 245 KBGS patients with height data, approximately 50% displayed short stature. Most patients showed a positive response to rhGH therapy, although the number of patients receiving treatment was limited. ANKRD11 deficiency potentially disrupts longitudinal bone growth by affecting the orderly differentiation of growth plate chondrocytes. Our review offers crucial insights into the association between ANKRD11 variants and short stature and provides valuable guidance for precise clinical diagnosis and treatment of patients with KBG syndrome.

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.

Properties and Mechanisms of Deletions, Insertions, and Substitutions in the Evolutionary History of SARS-CoV-2

SARS-CoV-2 has accumulated many mutations since its emergence in late 2019. Nucleotide substitutions leading to amino acid replacements constitute the primary material for natural selection. Insertions, deletions, and substitutions appear to be critical for coronavirus's macro- and microevolution. Understanding the molecular mechanisms of mutations in the mutational hotspots (positions, loci with recurrent mutations, and nucleotide context) is important for disentangling roles of mutagenesis and selection. In the SARS-CoV-2 genome, deletions and insertions are frequently associated with repetitive sequences, whereas C>U substitutions are often surrounded by nucleotides resembling the APOBEC mutable motifs. We describe various approaches to mutation spectra analyses, including the context features of RNAs that are likely to be involved in the generation of recurrent mutations. We also discuss the interplay between mutations and natural selection as a complex evolutionary trend. The substantial variability and complexity of pipelines for the reconstruction of mutations and the huge number of genomic sequences are major problems for the analyses of mutations in the SARS-CoV-2 genome. As a solution, we advocate for the development of a centralized database of predicted mutations, which needs to be updated on a regular basis.

Increases in computationally predicted deleterious variants of unknown significance and sperm mtDNA copy numbers may be associated with semen quality

BACKGROUND

Mitochondria are essential for sperm motility because they provide the energy required for the movement. Changes in sperm mtDNA, such as point mutations, large-scale deletions, or copy number variations, may interfere with ATP production and reduce sperm motility. However, it is not clear if changes in mtDNA are linked to semen quality.

OBJECTIVES

To explore the association between sperm mitochondrial DNA (mtDNA) changes and semen quality.

MATERIALS AND METHODS

Sixty-five oligo and/or astheno and/or terato patients (O/A/T) patients and 41 controls were recruited from couples undergoing assisted reproduction. Semen and blood samples were collected from the same individual on the day of oocyte retrieval to extract, isolate and purify mtDNA for next-generation sequencing. mtDNA copy numbers were assessed in 64 patient and 39 control sperm DNA samples using quantitative real-time PCR. The 4977 bp deletion was assessed in 20 patient and 20 control sperm DNA samples using polymerase chain reaction.

RESULTS

The mtDNA of patients was more likely to carry pathogenic variants or variants of unknown significance (VUSs) (P = 0.091) with higher heteroplasmy levels (P < 0.05) than that of controls. Interestingly, 33.85% of O/A/T patients (22 out of 65) lacked unique variants in their spermatozoa. but presented an exceptionally high mtDNA copy number (P < 0.0001). Moreover, we observed a decrease in the heteroplasmy level of common mtDNA variants shared by somatic and gamete cells (P < 0.0001) and the emergence of a very large number of de novo mtDNA variants with low-level heteroplasmy in spermatozoa.

DISCUSSION AND CONCLUSION

The increases in the number of computationally predicted deleterious VUS and mtDNA copies in spermatozoa may be associated with semen quality. Exposure to environmental mutation pressure that causes novel mtDNA variants with low-level heteroplasmy may occur during spermatogenesis. Furthermore, when a certain harmful threshold is reached, male germ cells may degrade mtDNA with mutations and replicate the correct mtDNA sequence to maintain the mitochondrial function in spermatozoa.

A new, all-encompassing aetiology of type 1 diabetes

The aetiology of type 1 diabetes (T1D) is considered multifactorial with the contribution of the MHC on chromosome 6 being most important. Multiple factors also contribute to the aetiology of colorectal neoplasia, but the final event causing the change from normal mucosa to polyp and from polyp to cancer is due to a single somatic mutation event. Repeated formation of colorectal neoplasia within an at-risk population results in a predictable, tapering, exponential neoplasia distribution. Critical mutations driving colorectal neoplasia formation occur in mutation-prone DNA. These observations led to three hypotheses related to T1D. First, a single somatic mutation within the MHC of antigen presenting cells results in a change in phenotype from normal to T1D. Second, the distribution of additional autoimmune diseases (AAIDs) among persons with T1D adheres to a predictable, tapering, exponential distribution. And third, critical mutations driving development of T1D occur in mutation-prone DNA. To address the hypotheses in an orderly fashion, a new analytical method called genome-wide aetiology analysis (GWEA) consisting of nine steps is presented. All data required for GWEA of T1D are obtained from peer-reviewed publications or publicly available genome and proteome databases. Critical GWEA steps include AAID distribution among persons with T1D, analysis of at-risk HLA loci for mutation-prone DNA, determination of the role of non-MHC genes on GWAS, and verification of human data by cell culture or animal experiments. GWEA results show that distribution of AAID among persons with T1D adheres to a predictable, tapering, exponential distribution. A single, critical, somatic mutation within the epitope-binding groove of at-risk HLA loci alters HLA-insulin-peptide-T-cell-receptor (TCR) complex binding affinity and creates a new pathway that leads to loss of self-tolerance. The at-risk HLA loci, in particular binding pockets P1, P4 and P9, are encoded by mutation-prone DNA: GC-rich DNA sequence and somatic hypermutation hotspots. All other genes on GWAS can but do not have to amplify the new autoimmune pathway by facilitating DNA mutations, changing peptide binding affinity, reducing signal inhibition or augmenting signal intensity. Animal experiments agree with human studies. In conclusion, T1D is caused by a somatic mutation within the epitope-binding groove of an at-risk HLA gene that affects HLA-insulin-peptide-TCR complex binding affinity and initiates an autoimmune pathway. The nature of the peptide that binds to a mutated epitope-binding groove of an at-risk HLA gene determines the type of autoimmune disease that develops, that is, one at-risk HLA locus, multiple autoimmune diseases. Thus, T1D and AAIDs, and therefore common autoimmune diseases, share a similar somatic mutation-based aetiology.

The Causes for Genomic Instability and How to Try and Reduce Them Through Rational Design of Synthetic DNA

Genetic engineering has revolutionized our ability to manipulate DNA and engineer organisms for various applications. However, this approach can lead to genomic instability, which can result in unwanted effects such as toxicity, mutagenesis, and reduced productivity. To overcome these challenges, smart design of synthetic DNA has emerged as a promising solution. By taking into consideration the intricate relationships between gene expression and cellular metabolism, researchers can design synthetic constructs that minimize metabolic stress on the host cell, reduce mutagenesis, and increase protein yield. In this chapter, we summarize the main challenges of genomic instability in genetic engineering and address the dangers of unknowingly incorporating genomically unstable sequences in synthetic DNA. We also demonstrate the instability of those sequences by the fact that they are selected against conserved sequences in nature. We highlight the benefits of using ESO, a tool for the rational design of DNA for avoiding genetically unstable sequences, and also summarize the main principles and working parameters of the software that allow maximizing its benefits and impact.

The functional mutational landscape of the lacZ gene

The lacZ gene of Escherichia coli encodes β-galactosidase (β-gal), a lactose metabolism enzyme of the lactose operon. Previous chemical modification or site-directed mutagenesis experiments have identified 21 amino acids that are essential for β-gal catalytic activity. We have assembled over 10,000 lacZ mutations from published studies that were collected using a positive selection assay to identify mutations in lacZ that disrupted β-gal function. We analyzed 6,465 independent lacZ mutations that resulted in 2,732 missense mutations that impaired β-gal function. Those mutations affected 492 of the 1,023 lacZ codons, including most of the 21 previously known residues critical for catalytic activity. Most missense mutations occurred near the catalytic site and in regions important for subunit tetramerization. Overall, our work provides a comprehensive and detailed map of the amino acid residues affecting the structure and catalytic activity of the β-gal enzyme.

Fighting the battle against evolution: designing genetically modified organisms for evolutionary stability

Synthetic biology has made significant progress in many areas, but a major challenge that has received limited attention is the evolutionary stability of synthetic constructs made of heterologous genes. The expression of these constructs in microorganisms, that is, production of proteins that are not necessary for the organism, is a metabolic burden, leading to a decrease in relative fitness and make the synthetic constructs unstable over time. This is a significant concern for the synthetic biology community, particularly when it comes to bringing this technology out of the laboratory. In this review, we discuss the issue of evolutionary stability in synthetic biology and review the available tools to address this challenge.

Glycidamide-induced hypermutation in yeast single-stranded DNA reveals a ubiquitous clock-like mutational motif in humans

Mutagens often prefer specific nucleotides or oligonucleotide motifs that can be revealed by studying the hypermutation spectra in single-stranded (ss) DNA. We utilized a yeast model to explore mutagenesis by glycidamide, a simple epoxide formed endogenously in humans from the environmental toxicant acrylamide. Glycidamide caused ssDNA hypermutation in yeast predominantly in cytosines and adenines. The most frequent mutations in adenines occurred in the nAt→nGt trinucleotide motif. Base substitutions A→G in this motif relied on Rev1 translesion polymerase activity. Inactivating Rev1 did not alter the nAt trinucleotide preference, suggesting it may be an intrinsic specificity of the chemical reaction between glycidamide and adenine in the ssDNA. We found this mutational motif enriched in published sequencing data from glycidamide-treated mouse cells and ubiquitous in human cancers. In cancers, this motif was positively correlated with the single base substitution (SBS) smoking-associated SBS4 signature, with the clock-like signatures SBS1, SBS5, and was strongly correlated with smoking history and with age of tumor donors. Clock-like feature of the motif was also revealed in cells of human skin and brain. Given its pervasiveness, we propose that this mutational motif reflects mutagenic lesions to adenines in ssDNA from a potentially broad range of endogenous and exogenous agents.

The roles of sexual selection and sexual conflict in shaping patterns of genome and transcriptome variation

Sexual dimorphism is one of the most prevalent, and often the most extreme, examples of phenotypic variation within species, and arises primarily from genomic variation that is shared between females and males. Many sexual dimorphisms arise through sex differences in gene expression, and sex-biased expression is one way that a single, shared genome can generate multiple, distinct phenotypes. Although many sexual dimorphisms are expected to result from sexual selection, and many studies have invoked the possible role of sexual selection to explain sex-specific traits, the role of sexual selection in the evolution of sexually dimorphic gene expression remains difficult to differentiate from other forms of sex-specific selection. In this Review, we propose a holistic framework for the study of sex-specific selection and transcriptome evolution. We advocate for a comparative approach, across tissues, developmental stages and species, which incorporates an understanding of the molecular mechanisms, including genomic variation and structure, governing gene expression. Such an approach is expected to yield substantial insights into the evolution of genetic variation and have important applications in a variety of fields, including ecology, evolution and behaviour.

Mutation bias and the predictability of evolution

Predicting evolutionary outcomes is an important research goal in a diversity of contexts. The focus of evolutionary forecasting is usually on adaptive processes, and efforts to improve prediction typically focus on selection. However, adaptive processes often rely on new mutations, which can be strongly influenced by predictable biases in mutation. Here, we provide an overview of existing theory and evidence for such mutation-biased adaptation and consider the implications of these results for the problem of prediction, in regard to topics such as the evolution of infectious diseases, resistance to biochemical agents, as well as cancer and other kinds of somatic evolution. We argue that empirical knowledge of mutational biases is likely to improve in the near future, and that this knowledge is readily applicable to the challenges of short-term prediction. This article is part of the theme issue 'Interdisciplinary approaches to predicting evolutionary biology'.

A nonsense mutation in C8orf37 linked with retinitis pigmentosa, early macular degeneration, cataract, and myopia in an arRP family from North India

OBJECTIVE

This study aimed at identifying the underlying genetic defect in a consanguineous autosomal recessive retinitis pigmentosa (arRP) (RP-1175) family having RP with early macular degeneration, cataract, and myopia.

METHODS

Whole-exome sequencing (WES) was performed on the DNA of the proband, and variants observed were validated in the rest of the affected and unaffected family members by Sanger sequencing. Different bioinformatics tools were applied to access the pathogenicity of the observed variant.

RESULTS

A nonsense mutation i.e., c.555G > A (p.Trp185Ter) in C8orf37 in homozygous form, has been identified that segregated with the disease in the affected members. c.555G > A was absent in unaffected family members and in 107 ethnically matched controls, therefore ruling out its possibility of being a polymorphism.

CONCLUSIONS

Present study identifies a nonsense mutation (c.555G > A) at codon 185 in C8orf37 linked with arRP, early macular degeneration, posterior subcapsular cataract, and myopia. The identical mutation has previously been reported in a Pakistani family with isolated RP and in a Chinese family with RP and macular degeneration. This variable expressivity of the identified mutation c.555G > A in C8orf37 in the analyzed Indian family may be attributed to the presence of the modifier alleles. Also, Trp185 might be a mutation hotspot in Asian arRP patients and in the future, p.Trp185Ter in C8orf37 may be tested during initial screening in arRP cases especially belonging to a similar population.

Cancer driver mutations: predictions and reality

Cancer cells accumulate many genetic alterations throughout their lifetime, but only a few of them drive cancer progression, termed driver mutations. Driver mutations may vary between cancer types and patients, can remain latent for a long time and become drivers at particular cancer stages, or may drive oncogenesis only in conjunction with other mutations. The high mutational, biochemical, and histological tumor heterogeneity makes driver mutation identification very challenging. In this review we summarize recent efforts to identify driver mutations in cancer and annotate their effects. We underline the success of computational methods to predict driver mutations in finding novel cancer biomarkers, including in circulating tumor DNA (ctDNA). We also report on the boundaries of their applicability in clinical research.

Improved NGS-based detection of microsatellite instability using tumor-only data

Microsatellite instability (MSI) is a molecular signature of mismatch repair deficiency (dMMR), a predictive marker of immune checkpoint inhibitor therapy response. Despite its recognized pan-cancer value, most methods only support detection of this signature in colorectal cancer. In addition to the tissue-specific differences that impact the sensitivity of MSI detection in other tissues, the performance of most methods is also affected by patient ethnicity, tumor content, and other sample-specific properties. These limitations are particularly important when only tumor samples are available and restrict the performance and adoption of MSI testing. Here we introduce MSIdetect, a novel solution for NGS-based MSI detection. MSIdetect models the impact of indel burden and tumor content on read coverage at a set of homopolymer regions that we found are minimally impacted by sample-specific factors. We validated MSIdetect in 139 Formalin-Fixed Paraffin-Embedded (FFPE) clinical samples from colorectal and endometrial cancer as well as other more challenging tumor types, such as glioma or sebaceous adenoma or carcinoma. Based on analysis of these samples, MSIdetect displays 100% specificity and 96.3% sensitivity. Limit of detection analysis supports that MSIdetect is sensitive even in samples with relatively low tumor content and limited microsatellite instability. Finally, the results obtained using MSIdetect in tumor-only data correlate well (R=0.988) with what is obtained using tumor-normal matched pairs, demonstrating that the solution addresses the challenges posed by MSI detection from tumor-only data. The accuracy of MSI detection by MSIdetect in different cancer types coupled with the flexibility afforded by NGS-based testing will support the adoption of MSI testing in the clinical setting and increase the number of patients identified that are likely to benefit from immune checkpoint inhibitor therapy.

Bulky Adducts in Clustered DNA Lesions: Causes of Resistance to the NER System

The nucleotide excision repair (NER) system removes a wide range of bulky DNA lesions that cause significant distortions of the regular double helix structure. These lesions, mainly bulky covalent DNA adducts, are induced by ultraviolet and ionizing radiation or the interaction between exogenous/endogenous chemically active substances and nitrogenous DNA bases. As the number of DNA lesions increases, e.g., due to intensive chemotherapy and combination therapy of various diseases or DNA repair impairment, clustered lesions containing bulky adducts may occur. Clustered lesions are two or more lesions located within one or two turns of the DNA helix. Despite the fact that repair of single DNA lesions by the NER system in eukaryotic cells has been studied quite thoroughly, the repair mechanism of these lesions in clusters remains obscure. Identification of the structural features of the DNA regions containing irreparable clustered lesions is of considerable interest, in particular due to a relationship between the efficiency of some antitumor drugs and the activity of cellular repair systems. In this review, we analyzed data on the induction of clustered lesions containing bulky adducts, the potential biological significance of these lesions, and methods for quantification of DNA lesions and considered the causes for the inhibition of NER-catalyzed excision of clustered bulky lesions.

Analyses of mutational patterns induced by formaldehyde and acetaldehyde reveal similarity to a common mutational signature

Formaldehyde and acetaldehyde are reactive small molecules produced endogenously in cells as well as being environmental contaminants. Both of these small aldehydes are classified as human carcinogens, since they are known to damage DNA and exposure is linked to cancer incidence. However, the mutagenic properties of formaldehyde and acetaldehyde remain incompletely understood, at least in part because they are relatively weak mutagens. Here, we use a highly sensitive yeast genetic reporter system featuring controlled generation of long single-stranded DNA regions to show that both small aldehydes induced mutational patterns characterized by predominantly C/G → A/T, C/G → T/A, and T/A → C/G substitutions, each in similar proportions. We observed an excess of C/G → A/T transversions when compared to mock-treated controls. Many of these C/G → A/T transversions occurred at TC/GA motifs. Interestingly, the formaldehyde mutational pattern resembles single base substitution signature 40 from the Catalog of Somatic Mutations in Cancer. Single base substitution signature 40 is a mutational signature of unknown etiology. We also noted that acetaldehyde treatment caused an excess of deletion events longer than 4 bases while formaldehyde did not. This latter result could be another distinguishing feature between the mutational patterns of these simple aldehydes. These findings shed new light on the characteristics of 2 important, commonly occurring mutagens.

Effect of Silver Nanoparticles on the In Vitro Regeneration, Biochemical, Genetic, and Phenotype Variation in Adventitious Shoots Produced from Leaf Explants in Chrysanthemum

Novel and unique properties of nanomaterials, which are not apparent in larger-size forms of the same material, encourage the undertaking of studies exploring the multifaced effects of nanomaterials on plants. The results of such studies are not only scientifically relevant but, additionally, can be implemented to plant production and/or breeding. This study aimed to verify the applicability of silver nanoparticles (AgNPs) as a mutagen in chrysanthemum breeding. Chrysanthemum × grandiflorum (Ramat.) Kitam. ‘Lilac Wonder’ and ‘Richmond’ leaf explants were cultured on the modified MS medium supplemented with 0.6 mg·L⁻¹ 6-benzylaminopurine (BAP) and 2 mg·L⁻¹ indole-3-acetic acid (IAA) and treated with AgNPs (spherical; 20 nm in diameter size; 0, 50, and 100 mg·L⁻¹). AgNPs strongly suppressed the capability of leaf explants to form adventitious shoots and the efficiency of shoot regeneration. The content of primary and secondary metabolites (chlorophyll a, chlorophyll b, total chlorophylls, carotenoids, anthocyanins, phenolic compounds) and the activity of enzymatic antioxidants (superoxide dismutase and guaiacol peroxide) in leaf explants varied depending on the AgNPs treatment and age of culture. Phenotype variations of ex vitro cultivated chrysanthemums, covering the color and pigment content in the inflorescence, were detected in one 50 mg·L⁻¹ AgNPs-derived and five 100 mg·L⁻¹ AgNPs-derived ‘Lilac Wonder’ plants and were manifested as the color change from pink to burgundy-gold. However, no changes in inflorescence color/shape were found among AgNPs-treated ‘Richmond’ chrysanthemums. On the other hand, the stem height, number of leaves, and chlorophyll content in leaves varied depending on the AgNPs treatment and the cultivar analyzed. A significant effect of AgNPs on the genetic variation occurrence was found. A nearly two-fold higher share of polymorphic products, in both cultivars studied, was generated by RAPD markers than by SCoTs. To conclude, protocols using leaf explant treatment with AgNPs can be used as a novel breeding technique in chrysanthemum. However, the individual cultivars may differ in biochemical response, the efficiency of in vitro regeneration, genetic variation, and frequency of induced mutations in flowering plants.

Biocontainment of Genetically Engineered Algae

Algae (including eukaryotic microalgae and cyanobacteria) have been genetically engineered to convert light and carbon dioxide to many industrially and commercially relevant chemicals including biofuels, materials, and nutritional products. At industrial scale, genetically engineered algae may be cultivated outdoors in open ponds or in closed photobioreactors. In either case, industry would need to address a potential risk of the release of the engineered algae into the natural environment, resulting in potential negative impacts to the environment. Genetic biocontainment strategies are therefore under development to reduce the probability that these engineered bacteria can survive outside of the laboratory or industrial setting. These include active strategies that aim to kill the escaped cells by expression of toxic proteins, and passive strategies that use knockouts of native genes to reduce fitness outside of the controlled environment of labs and industrial cultivation systems. Several biocontainment strategies have demonstrated escape frequencies below detection limits. However, they have typically done so in carefully controlled experiments which may fail to capture mechanisms of escape that may arise in the more complex natural environment. The selection of biocontainment strategies that can effectively kill cells outside the lab, while maintaining maximum productivity inside the lab and without the need for relatively expensive chemicals will benefit from further attention.

In silico evolution of nucleic acid-binding proteins from a nonfunctional scaffold

Directed evolution emulates the process of natural selection to produce proteins with improved or altered functions. These approaches have proven to be very powerful but are technically challenging and particularly time and resource intensive. To bypass these limitations, we constructed a system to perform the entire process of directed evolution in silico. We employed iterative computational cycles of mutation and evaluation to predict mutations that confer high-affinity binding activities for DNA and RNA to an initial de novo designed protein with no inherent function. Beneficial mutations revealed modes of nucleic acid recognition not previously observed in natural proteins, highlighting the ability of computational directed evolution to access new molecular functions. Furthermore, the process by which new functions were obtained closely resembles natural evolution and can provide insights into the contributions of mutation rate, population size and selective pressure on functionalization of macromolecules in nature.

Mutations as Levy flights

Data from a long time evolution experiment with Escherichia Coli and from a large study on copy number variations in subjects with European ancestry are analyzed in order to argue that mutations can be described as Levy flights in the mutation space. These Levy flights have at least two components: random single-base substitutions and large DNA rearrangements. From the data, we get estimations for the time rates of both events and the size distribution function of large rearrangements.

Short loop functional commonality identified in leukaemia proteome highlights crucial protein sub-networks

Direct drug targeting of mutated proteins in cancer is not always possible and efficacy can be nullified by compensating protein-protein interactions (PPIs). Here, we establish an in silico pipeline to identify specific PPI sub-networks containing mutated proteins as potential targets, which we apply to mutation data of four different leukaemias. Our method is based on extracting cyclic interactions of a small number of proteins topologically and functionally linked in the Protein-Protein Interaction Network (PPIN), which we call short loop network motifs (SLM). We uncover a new property of PPINs named 'short loop commonality' to measure indirect PPIs occurring via common SLM interactions. This detects 'modules' of PPI networks enriched with annotated biological functions of proteins containing mutation hotspots, exemplified by FLT3 and other receptor tyrosine kinase proteins. We further identify functional dependency or mutual exclusivity of short loop commonality pairs in large-scale cellular CRISPR-Cas9 knockout screening data. Our pipeline provides a new strategy for identifying new therapeutic targets for drug discovery.

Use of signals of positive and negative selection to distinguish cancer genes and passenger genes

A major goal of cancer genomics is to identify all genes that play critical roles in carcinogenesis. Most approaches focused on genes positively selected for mutations that drive carcinogenesis and neglected the role of negative selection. Some studies have actually concluded that negative selection has no role in cancer evolution. We have re-examined the role of negative selection in tumor evolution through the analysis of the patterns of somatic mutations affecting the coding sequences of human genes. Our analyses have confirmed that tumor suppressor genes are positively selected for inactivating mutations, oncogenes, however, were found to display signals of both negative selection for inactivating mutations and positive selection for activating mutations. Significantly, we have identified numerous human genes that show signs of strong negative selection during tumor evolution, suggesting that their functional integrity is essential for the growth and survival of tumor cells.

Unravelling roles of error-prone DNA polymerases in shaping cancer genomes

Mutagenesis is a key hallmark and enabling characteristic of cancer cells, yet the diverse underlying mutagenic mechanisms that shape cancer genomes are not understood. This review will consider the emerging challenge of determining how DNA damage response pathways-both tolerance and repair-act upon specific forms of DNA damage to generate mutations characteristic of tumors. DNA polymerases are typically the ultimate mutagenic effectors of DNA repair pathways. Therefore, understanding the contributions of DNA polymerases is critical to develop a more comprehensive picture of mutagenic mechanisms in tumors. Selection of an appropriate DNA polymerase-whether error-free or error-prone-for a particular DNA template is critical to the maintenance of genome stability. We review different modes of DNA polymerase dysregulation including mutation, polymorphism, and over-expression of the polymerases themselves or their associated activators. Based upon recent findings connecting DNA polymerases with specific mechanisms of mutagenesis, we propose that compensation for DNA repair defects by error-prone polymerases may be a general paradigm molding the mutational landscape of cancer cells. Notably, we demonstrate that correlation of error-prone polymerase expression with mutation burden in a subset of patient tumors from The Cancer Genome Atlas can identify mechanistic hypotheses for further testing. We contrast experimental approaches from broad, genome-wide strategies to approaches with a narrower focus on a few hundred base pairs of DNA. In addition, we consider recent developments in computational annotation of patient tumor data to identify patterns of mutagenesis. Finally, we discuss the innovations and future experiments that will develop a more comprehensive portrait of mutagenic mechanisms in human tumors.

HotSpot3D web server: an integrated resource for mutation analysis in protein 3D structures

MOTIVATION

HotSpot3D is a widely used software for identifying mutation hotspots on the 3D structures of proteins. To further assist users, we developed a new HotSpot3D web server to make this software more versatile, convenient and interactive.

RESULTS

The HotSpot3D web server performs data pre-processing, clustering, visualization and log-viewing on one stop. Users can interactively explore each cluster and easily re-visualize the mutational clusters within browsers. We also provide a database that allows users to search and visualize proximal mutations from 33 cancers in the Cancer Genome Atlas.

AVAILABILITY AND IMPLEMENTATION

http://niulab.scgrid.cn/HotSpot3D/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Quasi-species nature and differential gene expression of severe acute respiratory syndrome coronavirus 2 and phylogenetic analysis of a novel Iranian strain

A novel coronavirus related to severe acute respiratory syndrome virus, (SARS-CoV-2) is the causal agent of the COVID-19 pandemic. Despite the genetic mutations across the SARS-CoV-2 genome being recently investigated, its transcriptomic genetic polymorphisms at inter-host level and the viral gene expression level based on each Open Reading Frame (ORF) remains unclear. Using available High Throughput Sequencing (HTS) data and based on SARS-CoV-2 infected human transcriptomic data, this study presents a high-resolution map of SARS-CoV-2 single nucleotide polymorphism (SNP) hotspots in a viral population at inter-host level. Four throat swab samples from COVID-19 infected patients were pooled, with RNA-Seq read retrieved from SRA NCBI to detect 21 SNPs and a replacement across the SARS-CoV-2 genomic population. Twenty-two RNA modification sites on viral transcripts were identified that may cause inter-host genetic diversity of this virus. In addition, the canonical genomic RNAs of N ORF showed higher expression in transcriptomic data and reverse transcriptase quantitative PCR compared to other SARS-CoV-2 ORFs, indicating the importance of this ORF in virus replication or other major functions in virus cycle. Phylogenetic and ancestral sequence analyses based on the entire genome revealed that SARS-CoV-2 is possibly derived from a recombination event between SARS-CoV and Bat SARS-like CoV. Ancestor analysis of the isolates from different locations including Iran suggest shared Chinese ancestry. These results propose the importance of potential inter-host level genetic variations to the evolution of SARS-COV-2, and the formation of viral quasi-species. The RNA modifications discovered in this study may cause amino acid sequence changes in polyprotein, spike protein, product of ORF8 and nucleocapsid (N) protein, suggesting further insights to understanding the functional impacts of mutations in the life cycle and pathogenicity of SARS-CoV-2.

On the complexity measures of mutation hotspots in human TP53 protein

The role of sequence complexity in 23 051 somatic missense mutations including 73 well-known mutation hotspots across 22 major cancers was studied in human TP53 proteins. A role for sequence complexity in TP53 protein mutations is suggested since (i) the mutation rate significantly increases in low amino acid pair bias complexity; (ii) probability distribution complexity increases following single point substitution mutations and strikingly increases after mutation at the mutation hotspots including six detectable hotspot mutations (R175, G245, R248, R249, R273, and R282); and (iii) the degree of increase in distribution complexity is significantly correlated with the frequency of missense mutations (r = -0.5758, P < 0.0001) across 20 major types of solid tumors. These results are consistent with the hypothesis that amino acid pair bias and distribution probability may be used as novel measures for protein sequence complexity, and the degree of complexity is related to its susceptibility to mutation, as such, it may be used as a predictor for modeling protein mutations in human cancers.

Landscape of drug-resistance mutations in kinase regulatory hotspots

More than 48 kinase inhibitors (KIs) have been approved by Food and Drug Administration. However, drug-resistance (DR) eventually occurs, and secondary mutations have been found in the previously targeted primary-mutated cancer cells. Cancer and drug research communities recognize the importance of the kinase domain (KD) mutations for kinasopathies. So far, a systematic investigation of kinase mutations on DR hotspots has not been done yet. In this study, we systematically investigated four types of representative mutation hotspots (gatekeeper, G-loop, αC-helix and A-loop) associated with DR in 538 human protein kinases using large-scale cancer data sets (TCGA, ICGC, COSMIC and GDSC). Our results revealed 358 kinases harboring 3318 mutations that covered 702 drug resistance hotspot residues. Among them, 197 kinases had multiple genetic variants on each residue. We further computationally assessed and validated the epidermal growth factor receptor mutations on protein structure and drug-binding efficacy. This is the first study to provide a landscape view of DR-associated mutation hotspots in kinase's secondary structures, and its knowledge will help the development of effective next-generation KIs for better precision medicine.

Shoot Tip Cryopreservation of Lamprocapnos spectabilis (L.) Fukuhara Using Different Approaches and Evaluation of Stability on the Molecular, Biochemical, and Plant Architecture Levels

The aim of this study is to optimize and evaluate the effectiveness of vitrification, droplet-vitrification, and encapsulation-vitrification techniques in the cryopreservation of Lamprocapnos spectabilis (L.) Fukuhara 'Gold Heart', a popular medicinal and ornamental plant species. In vitro-derived shoot tips were used in the experiments. All three techniques were based on explant dehydration with plant vitrification solution 3 (PVS3; 50% glycerol and 50% sucrose) for 0, 30, 60, 90, 120, 150, or 180 min. The recovered microshoots were subjected to morphometric, biochemical, and molecular analyses (RAPD, ISSR, SCoT). The highest recovery level was reported with the encapsulation-vitrification protocol based on 150 min dehydration (73.1%), while the vitrification technique was the least effective (maximum 25.8% recovery). Explants cryopreserved with the encapsulation-vitrification technique produced the highest mean number of shoots (4.9); moreover, this technique was optimal in terms of rooting efficiency. The highest fresh weight of shoots, on the other hand, was found with the vitrification protocol based on a 30-min PVS3 treatment. The concentrations of chlorophyll a and b were lower in all cryopreservation-derived plants, compared to the untreated control. On the other hand, short dehydration and cryopreservation of non-encapsulated explants stimulated the synthesis of anthocyanins. A small genetic variation in 5% of all samples analyzed was detected by RAPD and ISSR marker systems. Only plants recovered from the encapsulation-vitrification protocol had no DNA sequence alternations.

Characterization of DNA methylation-based markers for human body fluid identification in forensics: a critical review

Body fluid identification in crime scene investigations aids in reconstruction of crime scenes. Several studies have identified and reported differentially methylated sites (DMSs) and regions (DMRs) which differ between forensically relevant tissues (tDMRs) and body fluids. Diverse factors affect methylation patterns such as the environment, diets, lifestyle, disease, ethnicity, genetic variation, amongst others. Thus, it is important to analyse the stability of markers employed for forensic identification. Furthermore, even though epigenetic modifications are described as stable and heritable, epigenetic inheritance of potential markers for body fluid identification needs to be assessed in the long term. Here, we discuss the current status of reported DNA methylation-based markers and their verification studies. Such thorough investigation is crucial to develop a stable panel of DNA methylation-based markers for accurate body fluid identification.

Genetic diversity and population structure of the Mediterranean sesame core collection with use of genome-wide SNPs developed by double digest RAD-Seq

The Mediterranean sesame core collection contains agro-morphologically superior sesame accessions from geographically diverse regions in four continents. In the present investigation, the genetic diversity and population structure of this collection was analyzed with 5292 high-quality SNPs discovered by double-digest restriction site associated DNA (ddRAD) sequencing, a cost-effective and flexible next-generation sequencing method. The genetic distance between pairs of accessions varied from 0.023 to 0.524. The gene diversity was higher in accessions from Asia than from America, Africa, and Europe. The highest genetic differentiation was observed between accessions collected from America and Europe. Structure analysis showed the presence of three subpopulations among the sesame accessions, and only six accessions were placed in an admixture group. Phylogenetic tree and principal coordinate analysis clustered the accessions based on their countries of origin. However, no clear division was evident among the sesame accessions with regard to their continental locations. This result was supported by an AMOVA analysis, which revealed a genetic variation among continental groups of 5.53% of the total variation. The large number of SNPs clearly indicated that the Mediterranean sesame core collection is a highly diverse genetic resource. The collection can be exploited by breeders to select appropriate accessions that will provide high genetic gain in sesame improvement programs. The high-quality SNP data generated here should also be used in genome-wide association studies to explore qualitative trait loci and SNPs related to economically and agronomically important traits in sesame.

High-resolution mapping of DNA polymerase fidelity using nucleotide imbalances and next-generation sequencing

DNA polymerase fidelity is affected by both intrinsic properties and environmental conditions. Current strategies for measuring DNA polymerase error rate in vitro are constrained by low error subtype sensitivity, poor scalability, and lack of flexibility in types of sequence contexts that can be tested. We have developed the Magnification via Nucleotide Imbalance Fidelity (MagNIFi) assay, a scalable next-generation sequencing assay that uses a biased deoxynucleotide pool to quantitatively shift error rates into a range where errors are frequent and hence measurement is robust, while still allowing for accurate mapping to error rates under typical conditions. This assay is compatible with a wide range of fidelity-modulating conditions, and enables high-throughput analysis of sequence context effects on base substitution and single nucleotide deletion fidelity using a built-in template library. We validate this assay by comparing to previously established fidelity metrics, and use it to investigate neighboring sequence-mediated effects on fidelity for several DNA polymerases. Through these demonstrations, we establish the MagNIFi assay for robust, high-throughput analysis of DNA polymerase fidelity.

Bending of DNA duplexes with mutation motifs

Mutations can be induced by environmental factors but also arise spontaneously during DNA replication or due to deamination of methylated cytosines at CpG dinucleotides. Sites where mutations occur with higher frequency than would be expected by chance are termed hotspots while sites that contain mutations rarely are termed coldspots. Mutations are permanently scanned and repaired by repair systems. Among them, the mismatch repair targets base pair mismatches, which are discriminated from canonical base pairs by probing altered elasticity of DNA. Using biased molecular dynamics simulations, we investigated the elasticity of coldspots and hotspots motifs detected in human genes associated with inherited disorders, and also of motifs with Czech population hotspots and de novo mutations. Main attention was paid to mutations leading to G/T and A+/C pairs. We observed that hotspots without CpG/CpHpG sequences are less flexible than coldspots, which indicates that flexible sequences are more effectively repaired. In contrary, hotspots with CpG/CpHpG sequences exhibited increased flexibility as coldspots. Their mutability is more likely related to spontaneous deamination of methylated cytosines leading to C > T mutations, which are primarily targeted by base excision repair. We corroborated conclusions based on computer simulations by measuring melting curves of hotspots and coldspots containing G/T mismatch.

Mutational signatures and mutable motifs in cancer genomes

Cancer is a genetic disorder, meaning that a plethora of different mutations, whether somatic or germ line, underlie the etiology of the 'Emperor of Maladies'. Point mutations, chromosomal rearrangements and copy number changes, whether they have occurred spontaneously in predisposed individuals or have been induced by intrinsic or extrinsic (environmental) mutagens, lead to the activation of oncogenes and inactivation of tumor suppressor genes, thereby promoting malignancy. This scenario has now been recognized and experimentally confirmed in a wide range of different contexts. Over the past decade, a surge in available sequencing technologies has allowed the sequencing of whole genomes from liquid malignancies and solid tumors belonging to different types and stages of cancer, giving birth to the new field of cancer genomics. One of the most striking discoveries has been that cancer genomes are highly enriched with mutations of specific kinds. It has been suggested that these mutations can be classified into 'families' based on their mutational signatures. A mutational signature may be regarded as a type of base substitution (e.g. C:G to T:A) within a particular context of neighboring nucleotide sequence (the bases upstream and/or downstream of the mutation). These mutational signatures, supplemented by mutable motifs (a wider mutational context), promise to help us to understand the nature of the mutational processes that operate during tumor evolution because they represent the footprints of interactions between DNA, mutagens and the enzymes of the repair/replication/modification pathways.

Exploring background mutational processes to decipher cancer genetic heterogeneity

Much remains unknown about the progression and heterogeneity of mutational processes in different cancers and their diagnostic and clinical potential. A growing body of evidence supports mutation rate dependence on the local DNA sequence context for various types of mutations. We propose several tools for the analysis of cancer context-dependent mutations, which are implemented in an online computational framework MutaGene. The framework explores DNA context-dependent mutational patterns and underlying somatic cancer mutagenesis, analyzes mutational profiles of cancer samples, identifies the combinations of underlying mutagenic processes including those related to infidelity of DNA replication and repair machinery, and various other endogenous and exogenous mutagenic factors. As a result, the combination of mutagenic processes can be identified in any query sample with subsequent comparison to mutational profiles derived from malignant and benign samples. In addition, mutagen or cancer-specific mutational background models are applied to calculate expected DNA and protein site mutability to decouple relative contributions of mutagenesis and selection in carcinogenesis, thus elucidating the site-specific driving events in cancer. MutaGene is freely available at https://www.ncbi.nlm.nih.gov/projects/mutagene/.

Genetic and molecular findings of 38 Iranian patients with chronic granulomatous disease caused by p47-phox defect

One of the components of NADPH oxidase is p47-phox, encoded by NCF1 gene. This study aims to find new genetic changes and clinical features in 38 Iranian patients with autosomal recessive chronic granulomatous disease (AR-CGD) caused by NCF1 gene defect. Patients who had abnormal NBT and DHR-1,2,3 assay with loss of p47-phox in Western blotting were included in this study. After recording demographic and clinical data, PCR amplification was performed followed by direct sequencing for all exons and exon-intron boundaries. The most common form of CGD in Iran was AR-CGD due to consanguinity marriages. Among patients with AR-CGD, NCF1 deficiency was found to be more common than other forms. Cutaneous involvements (53%), pulmonary infections (50%) and lymphadenopathy (29%) were more prevalent than other clinical manifestations of CGD. Mutation analysis of NCF1 gene identified five different mutations. Homozygous delta GT deletion (c.75_76delGT) was the most frequent mutation and was detected in more than 63% of families. Six families had a nonsense mutation in exon 7 (c.579G > A). Two novel mutations were found in exon 4 in two families, including a missense mutation (c.328C > T) and a nine-nucleotide deletion (c.331_339delTGTCCCCAC). Genetic detection of these mutations may result in early diagnosis and prevention of possible complications of the disease. This could be useful for timely decision-making for haematopoietic stem cell transplantation and for carrier detection as well as prenatal diagnosis of next children in the affected families. Our findings might help to predict outcomes, raise awareness and help effective treatment in these patients.

Electron injection from mitochondrial transcription factor A to DNA associated with thymine dimer photo repair

Electron transfer through π-stacked arrays of double-stranded DNA contributes to the redox chemistry of bases, including guanine oxidation and thymine-thymine dimer repair by photolyase. 5-Bromouracil is an attractive photoreactive thymine analogue that can be used to investigate electron transfer in DNA, and is a useful probe for protein-DNA interaction analysis. In the present study using ^BrU we found that UV irradiation facilitated electron injection from mitochondrial transcription factor A into DNA. We also observed that this electron injection could lead to repair of a thymine-thymine dimer.

Spatial Vulnerabilities of the Escherichia coli Genome to Spontaneous Mutations Revealed with Improved Duplex Sequencing

Investigation of spontaneous mutations by next-generation sequencing technology has attracted extensive attention lately due to the fundamental roles of spontaneous mutations in evolution and pathological processes. However, these studies only focused on the mutations accumulated through many generations during long-term (possibly be years of) culturing, but not the freshly generated mutations that occur at very low frequencies. In this study, we established a molecularly barcoded deep sequencing strategy to detect low abundant spontaneous mutations in genomes of bacteria cell cultures. Genome-wide spontaneous mutations in 15 Escherichia coli cell culture samples were defined with a high confidence (P < 0.01). We also developed a hotspot-calling approach based on the run-length encoding algorithm to find the genomic regions that are vulnerable to the spontaneous mutations. The hotspots for the mutations appeared to be highly conserved across the bacteria samples. Further biological annotation of these regions indicated that most of the spontaneous mutations were located at the repeat domains or nonfunctional domains of the genomes, suggesting the existence of mechanisms that could somehow prevent the occurrence of mutations in crucial genic areas. This study provides a more faithful picture of mutation occurrence and spectra in a single expansion process without long-term culturing.

The Presence of Genotoxic and/or Pro-inflammatory Bacterial Genes in Gut Metagenomic Databases and Their Possible Link With Inflammatory Bowel Diseases

Background: The human gut microbiota is a dynamic community of microorganisms that mediate important biochemical processes. Differences in the gut microbial composition have been associated with inflammatory bowel diseases (IBD) and other intestinal disorders. In this study, we quantified and compared the frequencies of eight genotoxic and/or pro-inflammatory bacterial genes found in metagenomic Whole Genome Sequences (mWGSs) of samples from individuals with IBD vs. a cohort of healthy human subjects. Methods: The eight selected gene sequences were clbN, clbB, cif, cnf-1, usp, tcpC from Escherichia coli, gelE from Enterococcus faecalis and murB from Akkermansia muciniphila. We also included the sequences for the conserved murB genes from E. coli and E. faecalis as markers for the presence of Enterobacteriaceae or Enterococci in the samples. The gene sequences were chosen based on their previously reported ability to disrupt normal cellular processes to either promote inflammation or to cause DNA damage in cultured cells or animal models, which could be linked to a role in IBD. The selected sequences were searched in three different mWGS datasets accessed through the Human Microbiome Project (HMP): a healthy cohort (N = 251), a Crohn's disease cohort (N = 60) and an ulcerative colitis cohort (N = 17). Results: Firstly, the sequences for the murB housekeeping genes from Enterobacteriaceae and Enterococci were more frequently found in the IBD cohorts (32% E. coli in IBD vs. 12% in healthy; 13% E. faecalis in IBD vs. 3% in healthy) than in the healthy cohort, confirming earlier reports of a higher presence of both of these taxa in IBD. For some of the sequences in our study, especially usp and gelE, their frequency was even more sharply increased in the IBD cohorts than in the healthy cohort, suggesting an association with IBD that is not easily explained by the increased presence of E. coli or E. faecalis in those samples. Conclusion: Our results suggest a significant association between the presence of some of these genotoxic or pro-inflammatory gene sequences and IBDs. In addition, these results illustrate the power and limitations of the HMP database in the detection of possible clinical correlations for individual bacterial genes.

Molecular Genetic Characterization of Mutagenesis Using a Highly Sensitive Single-Stranded DNA Reporter System in Budding Yeast

Mutations are permanent alterations to the coding content of DNA. They are starting material for the Darwinian evolution of species by natural selection, which has yielded an amazing diversity of life on Earth. Mutations can also be the fundamental basis of serious human maladies, most notably cancers. In this chapter, I describe a highly sensitive reporter system for the molecular genetic analysis of mutagenesis, featuring controlled generation of long stretches of single-stranded DNA in budding yeast cells. This system is ~100- to ~1000-fold more susceptible to mutation than conventional double-stranded DNA reporters, and is well suited for generating large mutational datasets to investigate the properties of mutagens.

DNA mutation motifs in the genes associated with inherited diseases

Mutations in human genes can be responsible for inherited genetic disorders and cancer. Mutations can arise due to environmental factors or spontaneously. It has been shown that certain DNA sequences are more prone to mutate. These sites are termed hotspots and exhibit a higher mutation frequency than expected by chance. In contrast, DNA sequences with lower mutation frequencies than expected by chance are termed coldspots. Mutation hotspots are usually derived from a mutation spectrum, which reflects particular population where an effect of a common ancestor plays a role. To detect coldspots/hotspots unaffected by population bias, we analysed the presence of germline mutations obtained from HGMD database in the 5-nucleotide segments repeatedly occurring in genes associated with common inherited disorders, in particular, the PAH, LDLR, CFTR, F8, and F9 genes. Statistically significant sequences (mutational motifs) rarely associated with mutations (coldspots) and frequently associated with mutations (hotspots) exhibited characteristic sequence patterns, e.g. coldspots contained purine tract while hotspots showed alternating purine-pyrimidine bases, often with the presence of CpG dinucleotide. Using molecular dynamics simulations and free energy calculations, we analysed the global bending properties of two selected coldspots and two hotspots with a G/T mismatch. We observed that the coldspots were inherently more flexible than the hotspots. We assume that this property might be critical for effective mismatch repair as DNA with a mutation recognized by MutSα protein is noticeably bent.

Fluorescence-Based Reporters for Detection of Mutagenesis in E. coli

Mutagenesis in model organisms following exposure to chemicals is used as an indicator of genotoxicity. Mutagenesis assays are also used to study mechanisms of DNA homeostasis. This chapter focuses on detection of mutagenesis in prokaryotes, which boils down to two approaches: reporter inactivation (forward mutation assay) and reversion of an inactivating mutation (reversion mutation assay). Both methods are labor intensive, involving visual screening, quantification of colonies on solid media, or determining a Poisson distribution in liquid culture. Here, we present two reversion reporters for in vivo mutagenesis that produce a quantitative output, and thus have the potential to greatly reduce the amount of test chemical and labor involved in these assays. This output is obtained by coupling a TEM β lactamase-based reversion assay with GFP fluorescence, either by placing the two genes on the same plasmid or by fusing them translationally and interrupting the N-terminus of the chimeric ORF with a stop codon. We also describe a reporter aimed at facilitating the monitoring of continuous mutagenesis in mutator strains. This reporter couples two reversion markers, allowing the temporal separation of mutation events in time, thus providing information about the dynamics of mutagenesis in mutator strains. Here, we describe these reporter systems, provide protocols for use, and demonstrate their key functional features using error-prone Pol I mutagenesis as a source of mutations.

A high-throughput assay for quantitative measurement of PCR errors

The accuracy with which DNA polymerase can replicate a template DNA sequence is an extremely important property that can vary by an order of magnitude from one enzyme to another. The rate of nucleotide misincorporation is shaped by multiple factors, including PCR conditions and proofreading capabilities, and proper assessment of polymerase error rate is essential for a wide range of sensitive PCR-based assays. In this paper, we describe a method for studying polymerase errors with exceptional resolution, which combines unique molecular identifier tagging and high-throughput sequencing. Our protocol is less laborious than commonly-used methods, and is also scalable, robust and accurate. In a series of nine PCR assays, we have measured a range of polymerase accuracies that is in line with previous observations. However, we were also able to comprehensively describe individual errors introduced by each polymerase after either 20 PCR cycles or a linear amplification, revealing specific substitution preferences and the diversity of PCR error frequency profiles. We also demonstrate that the detected high-frequency PCR errors are highly recurrent and that the position in the template sequence and polymerase-specific substitution preferences are among the major factors influencing the observed PCR error rate.

Chronic Lymphocytic Leukemia with Mutated IGHV4-34 Receptors: Shared and Distinct Immunogenetic Features and Clinical Outcomes

Purpose: We sought to investigate whether B cell receptor immunoglobulin (BcR IG) stereotypy is associated with particular clinicobiological features among chronic lymphocytic leukemia (CLL) patients expressing mutated BcR IG (M-CLL) encoded by the IGHV4-34 gene, and also ascertain whether these associations could refine prognostication.Experimental Design: In a series of 19,907 CLL cases with available immunogenetic information, we identified 339 IGHV4-34-expressing cases assigned to one of the four largest stereotyped M-CLL subsets, namely subsets #4, #16, #29 and #201, and investigated in detail their clinicobiological characteristics and disease outcomes.Results: We identified shared and subset-specific patterns of somatic hypermutation (SHM) among patients assigned to these subsets. The greatest similarity was observed between subsets #4 and #16, both including IgG-switched cases (IgG-CLL). In contrast, the least similarity was detected between subsets #16 and #201, the latter concerning IgM/D-expressing CLL. Significant differences between subsets also involved disease stage at diagnosis and the presence of specific genomic aberrations. IgG subsets #4 and #16 emerged as particularly indolent with a significantly (P < 0.05) longer time-to-first-treatment (TTFT; median TTFT: not yet reached) compared with the IgM/D subsets #29 and #201 (median TTFT: 11 and 12 years, respectively).Conclusions: Our findings support the notion that BcR IG stereotypy further refines prognostication in CLL, superseding the immunogenetic distinction based solely on SHM load. In addition, the observed distinct genetic aberration landscapes and clinical heterogeneity suggest that not all M-CLL cases are equal, prompting further research into the underlying biological background with the ultimate aim of tailored patient management. Clin Cancer Res; 23(17); 5292-301. ©2017 AACR.

Genetic Diversity, Population Structure, and Linkage Disequilibrium of a Core Collection of Ziziphus jujuba Assessed with Genome-wide SNPs Developed by Genotyping-by-sequencing and SSR Markers

Chinese jujube (Ziziphus jujuba Mill) is an economically important fruit species native to China with high nutritious and medicinal value. Genotyping-by-sequencing was used to detect and genotype single nucleotide polymorphisms (SNPs) in a core collection of 150 Chinese jujube accessions and further to characterize their genetic diversity, population structure, and linkage disequilibrium (LD). A total of 4,680 high-quality SNPs were identified, of which 38 sets of tri-allelic SNPs were detected. The average polymorphism information content (PIC) values based on bi-allelic SNPs and tri-allelic SNPs were 0.27 and 0.38, respectively. STRUCTURE and principal coordinate analyses based on SNPs revealed that the 150 accessions could be clustered into two groups. However, neighbor-joining trees indicated the accessions should be grouped into three major clusters. Our data confirm that the resolving power for genetic diversity was similar for the SSRs and SNPs. In contrast, regarding population structure, the resolving power was higher for SSRs than for SNPs. The LD pattern in Chinese jujube was investigated for the first time. We observed a relatively rapid LD decay with a short range (∼10 kb) for all pseudo-chromosomes and for individual pseudo-chromosomes. Our findings provide important information for future genome-wide association analyses and marker-assisted selective breeding of Chinese jujube.

Population Heterogeneity in Mutation Rate Increases the Frequency of Higher-Order Mutants and Reduces Long-Term Mutational Load

Mutation rate is a crucial evolutionary parameter that has typically been treated as a constant in population genetic analyses. However, the propensity to mutate is likely to vary among co-existing individuals within a population, due to genetic polymorphisms, heterogeneous environmental influences, and random physiological fluctuations. We review the evidence for mutation rate heterogeneity and explore its consequences by extending classic population genetic models to allow an arbitrary distribution of mutation rate among individuals, either with or without inheritance. With this general new framework, we rigorously establish the effects of heterogeneity at various evolutionary timescales. In a single generation, variation of mutation rate about the mean increases the probability of producing zero or many simultaneous mutations on a genome. Over multiple generations of mutation and selection, heterogeneity accelerates the appearance of both deleterious and beneficial multi-point mutants. At mutation-selection balance, higher-order mutant frequencies are likewise boosted, while lower-order mutants exhibit subtler effects; nonetheless, population mean fitness is always enhanced. We quantify the dependencies on moments of the mutation rate distribution and selection coefficients, and clarify the role of mutation rate inheritance. While typical methods of estimating mutation rate will recover only the population mean, analyses assuming mutation rate is fixed to this mean could underestimate the potential for multi-locus adaptation, including medically relevant evolution in pathogenic and cancerous populations. We discuss the potential to empirically parameterize mutation rate distributions, which have to date hardly been quantified.