Mechanism of a genetic glissando: structural biology of indel mutations

Deep indel mutagenesis reveals the impact of amino acid insertions and deletions on protein stability and function

Amino acid insertions and deletions (indels) are an abundant class of genetic variants. However, compared to substitutions, the effects of indels on protein stability are not well understood. To better understand indels here we analyse new and existing large-scale deep indel mutagenesis (DIM) of structurally diverse proteins. The effects of indels on protein stability vary extensively among and within proteins and are not well predicted by existing computational methods. To address this shortcoming we present INDELi, a series of models that combine experimental or predicted substitution effects and secondary structure information to provide good prediction of the effects of indels on both protein stability and pathogenicity. Moreover, quantifying the effects of indels on protein-protein interactions suggests that insertions can be an important class of gain-of-function variants. Our results provide an overview of the impact of indels on proteins and a method to predict their effects genome-wide.

Comparative genomics and phylogenetic analysis of mitochondrial genomes of Neocinnamomum

BACKGROUND

Neocinnamomum plants are considered a promising feedstock for biodiesel in China, due to the richness in long-chain fatty acids (LCFAs) found in their seeds. However, the mitochondrial genome (mitogenome) of this genus has not yet been systematically described, and the exploration of species relationships within this genus using mitogenome sequences is also an uncharted territory. This has hindered our understanding of mitogenome diversity and the evolutionary relationships within Neocinnamomum.

RESULTS

In this study, a total of 24 individuals representing seven distinct taxa from the genus Neocinnamomum were subjected to Illumina sequencing, and the species N. delavayi was sequenced using Oxford Nanopore sequencing technology. We successfully assembled the mitogenome of N. delavayi, which is 778,066 bp in size and exhibits a single circular structure. The analysis identified 659 dispersed repeats, 211 simple sequence repeats (SSRs), and 30 tandem repeats within the mitogenome. Additionally, 37 homologous fragments, totaling 9929 bp, were found between the mitogenome and the plastid genome (plastome). The codons of 41 protein-coding genes (PCGs) had a preference for ending in A/T, and the codon usage bias of the majority of these genes was influenced by natural selection pressures. Comparative genomic analysis revealed low collinearity and significant gene rearrangements between species. Phylogenetic analysis resulted in the classification of Neocinnamomum into six distinct clades, contradicting previous findings which based on complete plastomes and nuclear ribosomal cistron (nrDNA). In the PCGs of 24 individuals, 86 mutation events were identified, which included three indels and 83 SNPs. Notably, the ccmC gene underwent positive selection in pairwise comparisons of three species pairs. Furthermore, 748 RNA editing sites were predicted within the PCGs of the N. delavayi mitogenome.

CONCLUSIONS

This study enriches our knowledge of the mitogenomes in the family Lauraceae, and provides valuable data and a foundation for genomic evolution research, genetic resource conservation, and molecular breeding in Neocinnamomum.

Converging mechanism of UM171 and KBTBD4 neomorphic cancer mutations

Cancer mutations can create neomorphic protein-protein interactions to drive aberrant function1,2. As a substrate receptor of the CULLIN3-RING E3 ubiquitin ligase complex, KBTBD4 is recurrently mutated in medulloblastoma3, the most common embryonal brain tumour in children4. These mutations impart gain-of-function to KBTBD4 to induce aberrant degradation of the transcriptional corepressor CoREST5. However, their mechanism remains unresolved. Here we establish that KBTBD4 mutations promote CoREST degradation through engaging HDAC1/2 as the direct target of the mutant substrate receptor. Using deep mutational scanning, we chart the mutational landscape of the KBTBD4 cancer hotspot, revealing distinct preferences by which insertions and substitutions can promote gain-of-function and the critical residues involved in the hotspot interaction. Cryo-electron microscopy analysis of two distinct KBTBD4 cancer mutants bound to LSD1-HDAC1-CoREST reveals that a KBTBD4 homodimer asymmetrically engages HDAC1 with two KELCH-repeat β-propeller domains. The interface between HDAC1 and one of the KBTBD4 β-propellers is stabilized by the medulloblastoma mutations, which insert a bulky side chain into the HDAC1 active site pocket. Our structural and mutational analyses inform how this hotspot E3-neosubstrate interface can be chemically modulated. First, we unveil a converging shape-complementarity-based mechanism between gain-of-function E3 mutations and a molecular glue degrader, UM171. Second, we demonstrate that HDAC1/2 inhibitors can block the mutant KBTBD4-HDAC1 interface and proliferation of KBTBD4-mutant medulloblastoma cells. Altogether, our work reveals the structural and mechanistic basis of cancer mutation-driven neomorphic protein-protein interactions.

A complete DNA repair system assembled by two endosymbionts restores heat tolerance of the insect host

DNA repair systems are essential to maintain genome integrity and stability. Some obligate endosymbionts that experience long-term symbiosis with the insect hosts, however, have lost their key components for DNA repair. It is largely unexplored how the bacterial endosymbionts cope with the increased demand for mismatch repairs under heat stresses. Here, we showed that ibpA, a small heat shock protein encoded by Buchnera aphidicola, directly interacted with the cytoskeletal actin to prevent its aggregation in bacteriocytes, thus reinforcing the stability of bacteriocytes. However, the succession of 11 adenines in the promoter of ibpA is extremely prone to mismatching error, e.g., a single adenine deletion, which impairs the induction of ibpA under heat stress. Coinfection with a facultative endosymbiont Serratia symbiotica remarkably reduced the mutagenesis rate in the Buchnera genome and potentially prevented a single adenine deletion in ibpA promoter, thereby alleviating the heat vulnerability of aphid bacteriocytes. Furthermore, Serratia encoded mutH, a conserved core protein of prokaryotic DNA mismatch repair (MMR), accessed to Buchnera cells, which complemented Buchnera mutL and mutS in constituting an active MMR. Our findings imply that a full complement of a prokaryotic MMR system assembled by two bacterial endosymbionts contributes significantly to the thermostability of aphid bacteriocytes in an ibpA-dependent manner, furnishing a distinct molecular link among tripartite symbioses in shaping resilience and adaptation of their insect hosts to occupy other ecological niches.

The mutagenic consequences of defective DNA repair

Multiple separate repair mechanisms safeguard the genome against various types of DNA damage, and their failure can increase the rate of spontaneous mutagenesis. The malfunction of distinct repair mechanisms leads to genomic instability through different mutagenic processes. For example, defective mismatch repair causes high base substitution rates and microsatellite instability, whereas homologous recombination deficiency is characteristically associated with deletions and chromosome instability. This review presents a comprehensive collection of all mutagenic phenotypes associated with the loss of each DNA repair mechanism, drawing on data from a variety of model organisms and mutagenesis assays, and placing greatest emphasis on systematic analyses of human cancer datasets. We describe the latest theories on the mechanism of each mutagenic process, often explained by reliance on an alternative repair pathway or the error-prone replication of unrepaired, damaged DNA. Aided by the concept of mutational signatures, the genomic phenotypes can be used in cancer diagnosis to identify defective DNA repair pathways.

KBTBD4 Cancer Hotspot Mutations Drive Neomorphic Degradation of HDAC1/2 Corepressor Complexes

Cancer mutations can create neomorphic protein-protein interactions to drive aberrant function 1 . As a substrate receptor of the CULLIN3-RBX1 E3 ubiquitin ligase complex, KBTBD4 is recurrently mutated in medulloblastoma (MB) 2 , the most common embryonal brain tumor in children, and pineoblastoma 3 . These mutations impart gain-of-function to KBTBD4 to induce aberrant degradation of the transcriptional corepressor CoREST 4 . However, their mechanism of action remains unresolved. Here, we elucidate the mechanistic basis by which KBTBD4 mutations promote CoREST degradation through engaging HDAC1/2, the direct neomorphic target of the substrate receptor. Using deep mutational scanning, we systematically map the mutational landscape of the KBTBD4 cancer hotspot, revealing distinct preferences by which insertions and substitutions can promote gain-of-function and the critical residues involved in the hotspot interaction. Cryo-electron microscopy (cryo-EM) analysis of two distinct KBTBD4 cancer mutants bound to LSD1-HDAC1-CoREST reveals that a KBTBD4 homodimer asymmetrically engages HDAC1 with two KELCH-repeat propeller domains. The interface between HDAC1 and one of the KBTBD4 propellers is stabilized by the MB mutations, which directly insert a bulky side chain into the active site pocket of HDAC1. Our structural and mutational analyses inform how this hotspot E3-neo-substrate interface can be chemically modulated. First, our results unveil a converging shape complementarity-based mechanism between gain-of-function E3 mutations and a molecular glue degrader, UM171. Second, we demonstrate that HDAC1/2 inhibitors can block the mutant KBTBD4-HDAC1 interface, the aberrant degradation of CoREST, and the growth of KBTBD4-mutant MB models. Altogether, our work reveals the structural and mechanistic basis of cancer mutation-driven neomorphic protein-protein interactions and pharmacological strategies to modulate their action for therapeutic applications.

New whole-genome alignment tools are needed for tapping into plant diversity

Genome alignment is one of the most foundational methods for genome sequence studies. With rapid advances in sequencing and assembly technologies, these newly assembled genomes present challenges for alignment tools to meet the increased complexity and scale. Plant genome alignment is technologically challenging because of frequent whole-genome duplications (WGDs) as well as chromosome rearrangements and fractionation, high nucleotide diversity, widespread structural variation, and high transposable element (TE) activity causing large proportions of repeat elements. We summarize classical pairwise and multiple genome alignment (MGA) methods, and highlight techniques that are widely used or are being developed by the plant research community. We also outline the remaining challenges for precise genome alignment and the interpretation of alignment results in plants.

Galidesivir Triphosphate Promotes Stalling of Dengue-2 Virus Polymerase Immediately Prior to Incorporation

Millions of people are infected by the dengue and Zika viruses each year, resulting in significant morbidity and mortality. Galidesivir is an adenosine nucleoside analog that can attenuate flavivirus replication in cell-based assays and animal models of infection. Galidesivir is converted to the triphosphorylated form by host kinases and subsequently incorporated into viral RNA by viral RNA polymerases. This has been proposed to lead to the delayed termination of RNA synthesis. Here, we report direct in vitro testing of the effects of Galidesivir triphosphate on dengue-2 and Zika virus polymerase activity. Galidesivir triphosphate was chemically synthesized, and inhibition of RNA synthesis followed using a dinucleotide-primed assay with a homopolymeric poly(U) template. Galidesivir triphosphate was equipotent against dengue-2 and Zika polymerases, with IC50 values of 42 ± 12 μM and 47 ± 5 μM, respectively, at an ATP concentration of 20 μM. RNA primer extension assays show that the dengue-2 polymerase stalls while attempting to add a Galidesivir nucleotide to the nascent RNA chain, evidenced by the accumulation of RNA products truncated immediately upstream of Galidesivir incorporation sites. Nevertheless, Galidesivir is incorporated at isolated sites with low efficiency, leading to the subsequent synthesis of full-length RNA with no evidence of delayed chain termination. The incorporation of Galidesivir at consecutive sites is strongly disfavored, highlighting the potential for modulation of inhibitory effects of nucleoside analogs by the template sequence. Our results suggest that attenuation of dengue replication by Galidesivir may not derive from the early termination of RNA synthesis following Galidesivir incorporation.

TIVAN-indel: a computational framework for annotating and predicting non-coding regulatory small insertions and deletions

MOTIVATION

Small insertion and deletion (sindel) of human genome has an important implication for human disease. One important mechanism for non-coding sindel (nc-sindel) to have an impact on human diseases and phenotypes is through the regulation of gene expression. Nevertheless, current sequencing experiments may lack statistical power and resolution to pinpoint the functional sindel due to lower minor allele frequency or small effect size. As an alternative strategy, a supervised machine learning method can identify the otherwise masked functional sindels by predicting their regulatory potential directly. However, computational methods for annotating and predicting the regulatory sindels, especially in the non-coding regions, are underdeveloped.

RESULTS

By leveraging labeled nc-sindels identified by cis-expression quantitative trait loci analyses across 44 tissues in Genotype-Tissue Expression (GTEx), and a compilation of both generic functional annotations and large-scale epigenomic profiles, we develop TIssue-specific Variant Annotation for Non-coding indel (TIVAN-indel), which is a supervised computational framework for predicting non-coding regulatory sindels. As a result, we demonstrate that TIVAN-indel achieves the best prediction performance in both with-tissue prediction and cross-tissue prediction. As an independent evaluation, we train TIVAN-indel from the 'Whole Blood' tissue in GTEx and test the model using 15 immune cell types from an independent study named Database of Immune Cell Expression. Lastly, we perform an enrichment analysis for both true and predicted sindels in key regulatory regions such as chromatin interactions, open chromatin regions and histone modification sites, and find biologically meaningful enrichment patterns.

AVAILABILITY AND IMPLEMENTATION

https://github.com/lichen-lab/TIVAN-indel.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Insertions and deletions in protein evolution and engineering

Protein evolution or engineering studies are traditionally focused on amino acid substitutions and the way these contribute to fitness. Meanwhile, the insertion and deletion of amino acids is often overlooked, despite being one of the most common sources of genetic variation. Recent methodological advances and successful engineering stories have demonstrated that the time is ripe for greater emphasis on these mutations and their understudied effects. This review highlights the evolutionary importance and biotechnological relevance of insertions and deletions (indels). We provide a comprehensive overview of approaches that can be employed to include indels in random, (semi)-rational or computational protein engineering pipelines. Furthermore, we discuss the tolerance to indels at the structural level, address how domain indels can link the function of unrelated proteins, and feature studies that illustrate the surprising and intriguing potential of frameshift mutations.

Two New Mutations in the CEL Gene Causing Diabetes and Hereditary Pancreatitis: How to Correctly Identify MODY8 Cases

CONTEXT

Maturity onset diabetes of the young, type 8 (MODY8) is associated with mutations in the CEL gene, which encodes the digestive enzyme carboxyl ester lipase. Several diabetes cases and families have in recent years been attributed to mutations in CEL without any functional or clinical evidence provided.

OBJECTIVE

To facilitate correct MODY8 diagnostics, we screened 2 cohorts of diabetes patients and delineated the phenotype.

METHODS

Young, lean Swedish and Finnish patients with a diagnosis of type 2 diabetes (352 cases, 406 controls) were screened for mutations in the CEL gene. We also screened 58 Czech MODY cases who had tested negative for common MODY genes. For CEL mutation-positive subjects, family history was recorded, and clinical investigations and pancreatic imaging performed.

RESULTS

Two cases (1 Swedish and 1 Czech) with germline mutation in CEL were identified. Clinical and radiological investigations of these 2 probands and their families revealed dominantly inherited insulin-dependent diabetes, pancreatic exocrine dysfunction, and atrophic pancreas with lipomatosis and cysts. Notably, hereditary pancreatitis was the predominant phenotype in 1 pedigree. Both families carried single-base pair deletions in the proximal part of the CEL variable number of tandem repeat (VNTR) region in exon 11. The mutations are predicted to lead to aberrant protein tails that make the CEL protein susceptible to aggregation.

CONCLUSION

The diagnosis of MODY8 requires a pancreatic exocrine phenotype and a deletion in the CEL VNTR in addition to dominantly inherited diabetes. CEL screening may be warranted also in families with hereditary pancreatitis of unknown genetic etiology.

Ferrocene as a potential electrochemical reporting surrogate of abasic sites in DNA

We have evaluated the possibility of replacing abasic sites with ferrocene for enzymatic synthesis of canonical and modified DNA.

Meiosis and beyond - understanding the mechanistic and evolutionary processes shaping the germline genome

The separation of germ cell populations from the soma is part of the evolutionary transition to multicellularity. Only genetic information present in the germ cells will be inherited by future generations, and any molecular processes affecting the germline genome are therefore likely to be passed on. Despite its prevalence across taxonomic kingdoms, we are only starting to understand details of the underlying micro-evolutionary processes occurring at the germline genome level. These include segregation, recombination, mutation and selection and can occur at any stage during germline differentiation and mitotic germline proliferation to meiosis and post-meiotic gamete maturation. Selection acting on germ cells at any stage from the diploid germ cell to the haploid gametes may cause significant deviations from Mendelian inheritance and may be more widespread than previously assumed. The mechanisms that affect and potentially alter the genomic sequence and allele frequencies in the germline are pivotal to our understanding of heritability. With the rise of new sequencing technologies, we are now able to address some of these unanswered questions. In this review, we comment on the most recent developments in this field and identify current gaps in our knowledge.

Variable Spontaneous Mutation and Loss of Heterozygosity among Heterozygous Genomes in Yeast

Mutation and recombination are the primary sources of genetic variation. To better understand the evolution of genetic variation, it is crucial to comprehensively investigate the processes involving mutation accumulation and recombination. In this study, we performed mutation accumulation experiments on four heterozygous diploid yeast species in the Saccharomycodaceae family to determine spontaneous mutation rates, mutation spectra, and losses of heterozygosity (LOH). We observed substantial variation in mutation rates and mutation spectra. We also observed high LOH rates (1.65-11.07×10-6 events per heterozygous site per cell division). Biases in spontaneous mutation and LOH together with selection ultimately shape the variable genome-wide nucleotide landscape in yeast species.

Early developmental asymmetries in cell lineage trees in living individuals

Mosaic mutations can be used to track cell lineages in humans. We used cell cloning to analyze embryonic cell lineages in two living individuals and a postmortem human specimen. Of 10 reconstructed postzygotic divisions, none resulted in balanced contributions of daughter lineages to tissues. In both living individuals, one of two lineages from the first cleavage was dominant across tissues, with 90% frequency in blood. We propose that the efficiency of DNA repair contributes to lineage imbalance. Allocation of lineages in postmortem brain correlated with anterior-posterior axis, associating lineage history with cell fate choices in embryos. We establish a minimally invasive framework for defining cell lineages in any living individual, which paves the way for studying their relevance in health and disease.

Mononucleotide repeat expansions with non-natural polymerase substrates

Replicative strand slippage is a biological phenomenon, ubiquitous among different organisms. However, slippage events are also relevant to non-natural replication models utilizing synthetic polymerase substrates. Strand slippage may notably affect the outcome of the primer extension reaction with repetitive templates in the presence of non-natural nucleoside triphosphates. In the current paper, we studied the ability of Taq, Vent (exo-), and Deep Vent (exo-) polymerases to produce truncated, full size, or expanded modified strands utilizing non-natural 2′-deoxyuridine nucleotide analogues and different variants of the homopolymer template. Our data suggest that the slippage of the primer strand is dependent on the duplex fluttering, incorporation efficiency for a particular polymerase-dNTP pair, rate of non-templated base addition, and presence of competing nucleotides.

High-fidelity DNA ligation enforces accurate Okazaki fragment maturation during DNA replication

DNA ligase 1 (LIG1, Cdc9 in yeast) finalizes eukaryotic nuclear DNA replication by sealing Okazaki fragments using DNA end-joining reactions that strongly discriminate against incorrectly paired DNA substrates. Whether intrinsic ligation fidelity contributes to the accuracy of replication of the nuclear genome is unknown. Here, we show that an engineered low-fidelity LIG1Cdc9 variant confers a novel mutator phenotype in yeast typified by the accumulation of single base insertion mutations in homonucleotide runs. The rate at which these additions are generated increases upon concomitant inactivation of DNA mismatch repair, or by inactivation of the Fen1Rad27 Okazaki fragment maturation (OFM) nuclease. Biochemical and structural data establish that LIG1Cdc9 normally avoids erroneous ligation of DNA polymerase slippage products, and this protection is compromised by mutation of a LIG1Cdc9 high-fidelity metal binding site. Collectively, our data indicate that high-fidelity DNA ligation is required to prevent insertion mutations, and that this may be particularly critical following strand displacement synthesis during the completion of OFM.

Evolutionary history of cotranscriptional editing in the paramyxoviral phosphoprotein gene

The phosphoprotein gene of the paramyxoviruses encodes multiple protein products. The P, V, and W proteins are generated by transcriptional slippage. This process results in the insertion of non-templated guanosine nucleosides into the mRNA at a conserved edit site. The P protein is an essential component of the viral RNA polymerase and is encoded by a faithful copy of the gene in the majority of paramyxoviruses. However, in some cases, the non-essential V protein is encoded by default and guanosines must be inserted into the mRNA in order to encode P. The number of guanosines inserted into the P gene can be described by a probability distribution, which varies between viruses. In this article, we review the nature of these distributions, which can be inferred from mRNA sequencing data, and reconstruct the evolutionary history of cotranscriptional editing in the paramyxovirus family. Our model suggests that, throughout known history of the family, the system has switched from a P default to a V default mode four times; complete loss of the editing system has occurred twice, the canonical zinc finger domain of the V protein has been deleted or heavily mutated a further two times, and the W protein has independently evolved a novel function three times. Finally, we review the physical mechanisms of cotranscriptional editing via slippage of the viral RNA polymerase.

In vitro evolution of antibody affinity via insertional scanning mutagenesis of an entire antibody variable region

We report a systematic combinatorial exploration of affinity enhancement of antibodies by insertions and deletions (InDels). Transposon-based introduction of InDels via the method TRIAD (transposition-based random insertion and deletion mutagenesis) was used to generate large libraries with random in-frame InDels across the entire single-chain variable fragment gene that were further recombined and screened by ribosome display. Knowledge of potential insertion points from TRIAD libraries formed the basis of exploration of length and sequence diversity of novel insertions by insertional-scanning mutagenesis (InScaM). An overall 256-fold affinity improvement of an anti-IL-13 antibody BAK1 as a result of InDel mutagenesis and combination with known point mutations validates this approach, and suggests that the results of this InDel mutagenesis and conventional exploration of point mutations can synergize to generate antibodies with higher affinity.

Single-Cell Sequencing of Mouse Thymocytes Reveals Mutational Landscape Shaped by Replication Errors, Mismatch Repair, and H3K36me3

DNA mismatch repair (MMR) corrects replication errors and is recruited by the histone mark H3K36me3, enriched in exons of transcriptionally active genes. To dissect in vivo the mutational landscape shaped by these processes, we employed single-cell exome sequencing on T cells of wild-type and MMR-deficient (Mlh1-/-) mice. Within active genes, we uncovered a spatial bias in MMR efficiency: 3' exons, often H3K36me3-enriched, acquire significantly fewer MMR-dependent mutations compared with 5' exons. Huwe1 and Mcm7 genes, both active during lymphocyte development, stood out as mutational hotspots in MMR-deficient cells, demonstrating their intrinsic vulnerability to replication error in this cell type. Both genes are H3K36me3-enriched, which can explain MMR-mediated elimination of replication errors in wild-type cells. Thus, H3K36me3 can boost MMR in transcriptionally active regions, both locally and globally. This offers an attractive concept of thrifty MMR targeting, where critical genes in each cell type enjoy preferential shielding against de novo mutations.

RNAIndel: discovering somatic coding indels from tumor RNA-Seq data

MOTIVATION

Reliable identification of expressed somatic insertions/deletions (indels) is an unmet need due to artifacts generated in PCR-based RNA-Seq library preparation and the lack of normal RNA-Seq data, presenting analytical challenges for discovery of somatic indels in tumor transcriptome.

RESULTS

We present RNAIndel, a tool for predicting somatic, germline and artifact indels from tumor RNA-Seq data. RNAIndel leverages features derived from indel sequence context and biological effect in a machine-learning framework. Except for tumor samples with microsatellite instability, RNAIndel robustly predicts 88-100% of somatic indels in five diverse test datasets of pediatric and adult cancers, even recovering subclonal (VAF range 0.01-0.15) driver indels missed by targeted deep-sequencing, outperforming the current best-practice for RNA-Seq variant calling which had 57% sensitivity but with 14 times more false positives.

AVAILABILITY AND IMPLEMENTATION

RNAIndel is freely available at https://github.com/stjude/RNAIndel.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Transcription-coupled repair and mismatch repair contribute towards preserving genome integrity at mononucleotide repeat tracts

The mechanisms that underpin how insertions or deletions (indels) become fixed in DNA have primarily been ascribed to replication-related and/or double-strand break (DSB)-related processes. Here, we introduce a method to evaluate indels, orientating them relative to gene transcription. In so doing, we reveal a number of surprising findings: First, there is a transcriptional strand asymmetry in the distribution of mononucleotide repeat tracts in the reference human genome. Second, there is a strong transcriptional strand asymmetry of indels across 2,575 whole genome sequenced human cancers. We suggest that this is due to the activity of transcription-coupled nucleotide excision repair (TC-NER). Furthermore, TC-NER interacts with mismatch repair (MMR) under physiological conditions to produce strand bias. Finally, we show how insertions and deletions differ in their dependencies on these repair pathways. Our analytical approach reveals insights into the contribution of DNA repair towards indel mutagenesis in human cells.

Frameshifting preserves key physicochemical properties of proteins

Frameshifts in protein coding sequences are widely perceived as resulting in either nonfunctional or even deleterious protein products. Indeed, frameshifts typically lead to markedly altered protein sequences and premature stop codons. By analyzing complete proteomes from all three domains of life, we demonstrate that, in contrast, several key physicochemical properties of protein sequences exhibit significant robustness against +1 and -1 frameshifts. In particular, we show that hydrophobicity profiles of many protein sequences remain largely invariant upon frameshifting. For example, over 2,900 human proteins exhibit a Pearson's correlation coefficient R between the hydrophobicity profiles of the original and the +1-frameshifted variants greater than 0.7, despite an average sequence identity between the two of only 6.5% in this group. We observe a similar effect for protein sequence profiles of affinity for certain nucleobases as well as protein sequence profiles of intrinsic disorder. Finally, analysis of significance and optimality demonstrates that frameshift stability is embedded in the structure of the universal genetic code and may have contributed to shaping it. Our results suggest that frameshifting may be a powerful evolutionary mechanism for creating new proteins with vastly different sequences, yet similar physicochemical properties to the proteins from which they originate.

Proofreading of single nucleotide insertion/deletion replication errors analyzed by MALDI-TOF mass spectrometry assay

Small nucleotide insertion/deletion (indel) errors are one of the common replication errors in DNA synthesis. The most frequent occurrence of indel error was thought to be due to repeated sequences being prone to slippage during DNA replication. Proofreading and DNA mismatch repair are important factors in indel error correction to maintain the high fidelity of genetic information transactions. We employed a MALDI-TOF mass spectrometry (MS) analysis to measure the efficiency of Klenow polymerase (KF) proofreading of indel errors. Herein, a non-labeled and non-radio-isotopic oligonucleotide primer is annealed to a template DNA forming a single nucleotide indel error and was proofread by KF in the presence of a combination of different deoxyribonucleotide triphosphates and/or dideoxyribonucleotide triphosphates. The proofreading products were identified by the KF modified mass change of the primer. We examined proofreading of DNAs containing indel errors at various positions of the primer-template junction. We found that indel errors located 1-5-nucleotides (nt) from the primer terminus can be proofread efficiently, while insertion/deletions at 6-nt from the 3' end are partially corrected and extended. Indels located 7-9-nt from the primer terminus escape proofreading and are elongated by polymerase. The possible underlying mechanisms of these observations are discussed in the context of the polymerase and primer-template junction interactions via a structure analysis.

Mode and Tempo of Microsatellite Length Change in a Malaria Parasite Mutation Accumulation Experiment

Malaria parasites have small extremely AT-rich genomes: microsatellite repeats (1–9 bp) comprise 11% of the genome and genetic variation in natural populations is dominated by repeat changes in microsatellites rather than point mutations. This experiment was designed to quantify microsatellite mutation patterns in Plasmodium falciparum. We established 31 parasite cultures derived from a single parasite cell and maintained these for 114–267 days with frequent reductions to a single cell, so parasites accumulated mutations during ∼13,207 cell divisions. We Illumina sequenced the genomes of both progenitor and end-point mutation accumulation (MA) parasite lines in duplicate to validate stringent calling parameters. Microsatellite calls were 99.89% (GATK), 99.99% (freeBayes), and 99.96% (HipSTR) concordant in duplicate sequence runs from independent sequence libraries, whereas introduction of microsatellite mutations into the reference genome revealed a low false negative calling rate (0.68%). We observed 98 microsatellite mutations. We highlight several conclusions: microsatellite mutation rates (3.12 × 10⁻⁷ to 2.16 × 10⁻⁸/cell division) are associated with both repeat number and repeat motif like other organisms studied. However, 41% of changes resulted from loss or gain of more than one repeat: this was particularly true for long repeat arrays. Unlike other eukaryotes, we found no insertions or deletions that were not associated with repeats or homology regions. Overall, microsatellite mutation rates are among the lowest recorded and comparable to those in another AT-rich protozoan (Dictyostelium). However, a single infection (>10¹¹ parasites) will still contain over 2.16 × 10³ to 3.12 × 10⁴ independent mutations at any single microsatellite locus.

Elucidation of de novo small insertion/deletion biology with parent-of-origin phasing

The mechanisms underlying de novo insertion/deletion (indel) genesis, such as polymerase slippage, have been hypothesized but not well characterized in the human genome. We implemented two methodological improvements, which were leveraged to dissect indel mutagenesis. We assigned de novo variants to parent-of-origin (i.e., phasing) with low-coverage long-read whole-genome sequencing, achieving better phasing compared to short-read sequencing (medians of 84% and 23%, respectively). We then wrote an application programming interface to classify indels into three subtypes according to sequence context. Across three cohorts with different phasing methods (N_trios  = 540, all cohorts), we observed that one de novo indel subtype, change in copy count (CCC), was significantly correlated with father's (p = 7.1 × 10^-4 ) but not mother's (p = .45) age at conception. We replicated this effect in three cohorts without de novo phasing (p_paternal  = 1.9 × 10^-9 , p_maternal  = .61; N_trios  = 3,391, all cohorts). Although this is consistent with polymerase slippage during spermatogenesis, the percentage of variance explained by paternal age was low, and we did not observe an association with replication timing. These results suggest that spermatogenesis-specific events have a minor role in CCC indel mutagenesis, one not observed for other indel subtypes nor for maternal age in general. These results have implications for indel modeling in evolution and disease.

Causes and Effects of Loss of Classical Nonhomologous End Joining Pathway in Parasitic Eukaryotes

We report frequent losses of components of the classical nonhomologous end joining pathway (C-NHEJ), one of the main eukaryotic tools for end joining repair of DNA double-strand breaks, in several lineages of parasitic protists. Moreover, we have identified a single lineage among trypanosomatid flagellates that has lost Ku70 and Ku80, the core C-NHEJ components, and accumulated numerous insertions in many protein-coding genes. We propose a correlation between these two phenomena and discuss the possible impact of the C-NHEJ loss on genome evolution and transition to the parasitic lifestyle.IMPORTANCE Parasites tend to evolve small and compact genomes, generally endowed with a high mutation rate, compared with those of their free-living relatives. However, the mechanisms by which they achieve these features, independently in unrelated lineages, remain largely unknown. We argue that the loss of the classical nonhomologous end joining pathway components may be one of the crucial steps responsible for characteristic features of parasite genomes.

Phylogenetic Signal of Indels and the Neoavian Radiation

The early radiation of Neoaves has been hypothesized to be an intractable “hard polytomy”. We explore the fundamental properties of insertion/deletion alleles (indels), an under-utilized form of genomic data with the potential to help solve this. We scored >5 million indels from >7000 pan-genomic intronic and ultraconserved element (UCE) loci in 48 representatives of all neoavian orders. We found that intronic and UCE indels exhibited less homoplasy than nucleotide (nt) data. Gene trees estimated using indel data were less resolved than those estimated using nt data. Nevertheless, Accurate Species TRee Algorithm (ASTRAL) species trees estimated using indels were generally similar to nt-based ASTRAL trees, albeit with lower support. However, the power of indel gene trees became clear when we combined them with nt gene trees, including a striking result for UCEs. The individual UCE indel and nt ASTRAL trees were incongruent with each other and with the intron ASTRAL trees; however, the combined indel+nt ASTRAL tree was much more congruent with the intronic trees. Finally, combining indel and nt data for both introns and UCEs provided sufficient power to reduce the scope of the polytomy that was previously proposed for several supraordinal lineages of Neoaves.

Pathogenicity and functional impact of non-frameshifting insertion/deletion variation in the human genome

Differentiation between phenotypically neutral and disease-causing genetic variation remains an open and relevant problem. Among different types of variation, non-frameshifting insertions and deletions (indels) represent an understudied group with widespread phenotypic consequences. To address this challenge, we present a machine learning method, MutPred-Indel, that predicts pathogenicity and identifies types of functional residues impacted by non-frameshifting insertion/deletion variation. The model shows good predictive performance as well as the ability to identify impacted structural and functional residues including secondary structure, intrinsic disorder, metal and macromolecular binding, post-translational modifications, allosteric sites, and catalytic residues. We identify structural and functional mechanisms impacted preferentially by germline variation from the Human Gene Mutation Database, recurrent somatic variation from COSMIC in the context of different cancers, as well as de novo variants from families with autism spectrum disorder. Further, the distributions of pathogenicity prediction scores generated by MutPred-Indel are shown to differentiate highly recurrent from non-recurrent somatic variation. Collectively, we present a framework to facilitate the interrogation of both pathogenicity and the functional effects of non-frameshifting insertion/deletion variants. The MutPred-Indel webserver is available at http://mutpred.mutdb.org/.

R-loop formation by dCas9 is mutagenic in Saccharomyces cerevisiae

Cas9 binds and cleaves specific DNA sequences by inducing the formation of an R-loop between the guide RNA and its genomic target. While targeting of active Cas9 to a genomic locus is highly mutagenic because Cas9 creates DNA double strand breaks, targeting of dead Cas9 (dCas9) is presumed not to be mutagenic, as dCas9 lacks DNA endonuclease activity. Here, we show that dCas9 targeting induces mutations in yeast, particularly when targeted to the non-transcribed strand of a gene. dCas9-induced mutations cluster near the guide RNA target region and are comprised of single nucleotide substitutions, small insertions and deletions, and even complex mutations, depending upon the particular guide RNA target. We show that many of these mutations are a consequence of cytosine deamination events occurring on the non-target strand of the dCas9-induced R-loop, while others are associated with homopolymer instability or translesion DNA synthesis. Targeting of dCas9 by a mismatch-containing guide RNA also increases CAN1 mutation frequency, particularly in an ung1Δ mutant strain, suggesting that dCas9 induces mutations through similar mechanisms at off-target sites. These findings indicate that DNA binding by dCas9 is mutagenic in yeast, likely because dCas9 induces the formation of an R-loop at its target site.

Patterns of Nucleotide Deletion and Insertion Inferred from Bacterial Pseudogenes

Pseudogenes are a paradigm of neutral evolution and their study has the potential to reveal intrinsic mutational biases. However, this potential is mitigated by the fact that pseudogenes are quickly purged from bacterial genomes. Here, we assembled a large set of pseudogenes from genomes experiencing reductive evolution as well as functional references for which we could establish reliable phylogenetic relationships. Using this unique dataset, we identified 857 independent insertion and deletion mutations and discover a pervasive bias towards deletions, but not insertions, with sizes multiples of 3 nt. We further show that selective constraints for the preservation of gene frame are unlikely to account for the observed mutational bias and propose that a mechanistic bias in alternative end-joining repair, a recombination-independent double strand break DNA repair mechanism, is responsible for the accumulation of 3n deletions.

A Fluorescent G-Quadruplex Sensor for Chemical RNA Copying

Non-enzymatic RNA replication may have been one of the processes involved in the appearance of life on Earth. Attempts to recreate this process in a laboratory setting have not been successful thus far, highlighting a critical need for finding prebiotic conditions that increase the rate and the yield. Now a highly parallel assay for template directed RNA synthesis is presented that relies on the intrinsic fluorescence of a 2-aminopurine modified G-quadruplex. The application of the assay to examine the combined influence of multiple variables including pH, divalent metal concentrations and ribonucleotide concentrations on the copying of RNA sequences is demonstrated. The assay enables a direct survey of physical and chemical conditions, potentially prebiotic, which could enable the chemical replication of RNA.

Muver, a computational framework for accurately calling accumulated mutations

BACKGROUND

Identification of mutations from next-generation sequencing data typically requires a balance between sensitivity and accuracy. This is particularly true of DNA insertions and deletions (indels), that can impart significant phenotypic consequences on cells but are harder to call than substitution mutations from whole genome mutation accumulation experiments. To overcome these difficulties, we present muver, a computational framework that integrates established bioinformatics tools with novel analytical methods to generate mutation calls with the extremely low false positive rates and high sensitivity required for accurate mutation rate determination and comparison.

RESULTS

Muver uses statistical comparison of ancestral and descendant allelic frequencies to identify variant loci and assigns genotypes with models that include per-sample assessments of sequencing errors by mutation type and repeat context. Muver identifies maximally parsimonious mutation pathways that connect these genotypes, differentiating potential allelic conversion events and delineating ambiguities in mutation location, type, and size. Benchmarking with a human gold standard father-son pair demonstrates muver's sensitivity and low false positive rates. In DNA mismatch repair (MMR) deficient Saccharomyces cerevisiae, muver detects multi-base deletions in homopolymers longer than the replicative polymerase footprint at rates greater than predicted for sequential single-base deletions, implying a novel multi-repeat-unit slippage mechanism.

CONCLUSIONS

Benchmarking results demonstrate the high accuracy and sensitivity achieved with muver, particularly for indels, relative to available tools. Applied to an MMR-deficient Saccharomyces cerevisiae system, muver mutation calls facilitate mechanistic insights into DNA replication fidelity.

Site-Specific Self-Catalyzed DNA Depurination: A Biological Mechanism That Leads to Mutations and Creates Sequence Diversity

Self-catalyzed DNA depurination is a sequence-specific physiological mechanism mediated by spontaneous extrusion of a stem-loop catalytic intermediate. Hydrolysis of the 5'G residue of the 5'GA/TGG loop and of the first 5'A residue of the 5'GAGA loop, together with particular first stem base pairs, specifies their hydrolysis without involving protein, cofactor, or cation. As such, this mechanism is the only known DNA catalytic activity exploited by nature. The consensus sequences for self-depurination of such G- and A-loop residues occur in all genomes examined across the phyla, averaging one site every 2,000-4,000 base pairs. Because apurinic sites are subject to error-prone repair, leading to substitution and short frameshift mutations, they are both a source of genome damage and a means for creating sequence diversity. Their marked overrepresentation in genomes, and largely unchanging density from the lowest to the highest organisms, indicate their selection over the course of evolution. The mutagenicity at such sites in many human genes is associated with loss of function of key proteins responsible for diverse diseases.

Error-prone replication bypass of the imidazole ring-opened formamidopyrimidine deoxyguanosine adduct

Addition of hydroxyl radicals to the C8 position of 2'-deoxyguanosine generates an 8-hydroxyguanyl radical that can be converted into either 8-oxo-7,8-dihydro-2'-deoxyguanosine or N-(2-deoxy-d-pentofuranosyl)-N-(2,6-diamino-4-hydroxy-5-formamidopyrimidine) (Fapy-dG). The Fapy-dG adduct can adopt different conformations and in particular, can exist in an unnatural α anomeric configuration in addition to canonical β configuration. Previous studies reported that in 5'-TGN-3' sequences, Fapy-dG predominantly induced G → T transversions in both mammalian cells and Escherichia coli, suggesting that mutations could be formed either via insertion of a dA opposite the 5' dT due to primer/template misalignment or as result of direct miscoding. To address this question, single-stranded vectors containing a site-specific Fapy-dG adduct were generated to vary the identity of the 5' nucleotide. Following vector replication in primate cells (COS7), complex mutation spectra were observed that included ∼3-5% G → T transversions and ∼14-21% G → A transitions. There was no correlation apparent between the identity of the 5' nucleotide and spectra of mutations. When conditions for vector preparation were modified to favor the β anomer, frequencies of both G → T and G → A substitutions were significantly reduced. Mutation frequencies in wild-type E. coli and a mutant deficient in damage-inducible DNA polymerases were significantly lower than detected in COS7 and spectra were dominated by deletions. Thus, mutagenic bypass of Fapy-dG can proceed via mechanisms that are different from the previously proposed primer/template misalignment or direct misinsertions of dA or dT opposite to the β anomer of Fapy-dG. Environ. Mol. Mutagen. 58:182-189, 2017. © 2017 Wiley Periodicals, Inc.

Targeted next generation sequencing of the entire vitamin D receptor gene reveals polymorphisms correlated with vitamin D deficiency among older Filipino women with and without fragility fracture

This study aimed to discover genetic variants in the entire 101 kB vitamin D receptor (VDR) gene for vitamin D deficiency in a group of postmenopausal Filipino women using targeted next generation sequencing (TNGS) approach in a case-control study design. A total of 50 women with and without osteoporotic fracture seen at the Philippine Orthopedic Center were included. Blood samples were collected for determination of serum vitamin D, calcium, phosphorus, glucose, blood urea nitrogen, creatinine, aspartate aminotransferase, alanine aminotransferase and as primary source for targeted VDR gene sequencing using the Ion Torrent Personal Genome Machine. The variant calling was based on the GATK best practice workflow and annotated using Annovar tool. A total of 1496 unique variants in the whole 101-kb VDR gene were identified. Novel sequence variations not registered in the dbSNP database were found among cases and controls at a rate of 23.1% and 16.6% of total discovered variants, respectively. One disease-associated enhancer showed statistically significant association to low serum 25-hydroxy vitamin D levels (Pearson chi-square P-value=0.009). The transcription factor binding site prediction program PROMO predicted the disruption of three transcription factor binding sites in this enhancer region. These findings show the power of TNGS in identifying sequence variations in a very large gene and the surprising results obtained in this study greatly expand the catalog of known VDR sequence variants that may represent an important clue in the emergence of vitamin D deficiency. Such information will also provide the additional guidance necessary toward a personalized nutritional advice to reach sufficient vitamin D status.

Simple sequence repeat variations expedite phage divergence: Mechanisms of indels and gene mutations

Phages are the most abundant biological entities and influence prokaryotic communities on Earth. Comparing closely related genomes sheds light on molecular events shaping phage evolution. Simple sequence repeat (SSR) variations impart over half of the genomic changes between T7M and T3, indicating an important role of SSRs in accelerating phage genetic divergence. Differences in coding and noncoding regions of phages infecting different hosts, coliphages T7M and T3, Yersinia phage ϕYeO3-12, and Salmonella phage ϕSG-JL2, frequently arise from SSR variations. Such variations modify noncoding and coding regions; the latter efficiently changes multiple amino acids, thereby hastening protein evolution. Four classes of events are found to drive SSR variations: insertion/deletion of SSR units, expansion/contraction of SSRs without alteration of genome length, changes of repeat motifs, and generation/loss of repeats. The categorization demonstrates the ways SSRs mutate in genomes during phage evolution. Indels are common constituents of genome variations and human diseases, yet, how they occur without preexisting repeat sequence is less understood. Non-repeat-unit-based misalignment-elongation (NRUBME) is proposed to be one mechanism for indels without adjacent repeats. NRUBME or consecutive NRUBME may also change repeat motifs or generate new repeats. NRUBME invoking a non-Watson-Crick base pair explains insertions that initiate mononucleotide repeats. Furthermore, NRUBME successfully interprets many inexplicable human di- to tetranucleotide repeat generations. This study provides the first evidence of SSR variations expediting phage divergence, and enables insights into the events and mechanisms of genome evolution. NRUBME allows us to emulate natural evolution to design indels for various applications.

Genome-wide analysis of the specificity and mechanisms of replication infidelity driven by imbalanced dNTP pools

The absolute and relative concentrations of the four dNTPs are key determinants of DNA replication fidelity, yet the consequences of altered dNTP pools on replication fidelity have not previously been investigated on a genome-wide scale. Here, we use deep sequencing to determine the types, rates and locations of uncorrected replication errors that accumulate in the nuclear genome of a mismatch repair-deficient diploid yeast strain with elevated dCTP and dTTP concentrations. These imbalanced dNTP pools promote replication errors in specific DNA sequence motifs suggesting increased misinsertion and increased mismatch extension at the expense of proofreading. Interestingly, substitution rates are similar for leading and lagging strand replication, but are higher in regions replicated late in S phase. Remarkably, the rate of single base deletions is preferentially increased in coding sequences and in short rather than long mononucleotides runs. Based on DNA sequence motifs, we propose two distinct mechanisms for generating single base deletions in vivo. Collectively, the results indicate that elevated dCTP and dTTP pools increase mismatch formation and decrease error correction across the nuclear genome, and most strongly increases mutation rates in coding and late replicating sequences.

Identification and characterization of a -1 reading frameshift in the heavy chain constant region of an IgG1 recombinant monoclonal antibody produced in CHO cells

Frameshifts lead to complete alteration of the intended amino acid sequences, and therefore may affect the biological activities of protein therapeutics and pose potential immunogenicity risks. We report here the identification and characterization of a novel -1 frameshift variant in a recombinant IgG1 therapeutic monoclonal antibody (mAb) produced in Chinese hamster ovary cells during the cell line selection studies. The variant was initially observed as an atypical post-monomer fragment peak in size exclusion chromatography. Characterization of the fragment peak using intact and reduced liquid chromatography-mass spectrometry (LC-MS) analyses determined that the fragment consisted of a normal light chain disulfide-linked to an aberrant 26 kDa fragment that could not be assigned to any HC fragment even after considering common modifications. Further analysis using LC-MS/MS peptide mapping revealed that the aberrant fragment contained the expected HC amino acid sequence (1-232) followed by a 20-mer novel sequence corresponding to expression of heavy chain DNA sequence in the -1 reading frame. Examination of the DNA sequence around the frameshift initiation site revealed that a mononucleotide repeat GGGGGG located in the IgG1 HC constant region was most likely the structural root cause of the frameshift. Rapid identification of the frameshift allowed us to avoid use of a problematic cell line containing the frameshift as the production cell line. The frameshift reported here may be observed in other mAb products and the hypothesis-driven analytical approaches employed here may be valuable for rapid identification and characterization of frameshift variants in other recombinant proteins.

A small indel mutation in an anthocyanin transporter causes variegated colouration of peach flowers

The ornamental peach cultivar 'Hongbaihuatao (HBH)' can simultaneously bear pink, red, and variegated flowers on a single tree. Anthocyanin content in pink flowers is extremely low, being only 10% that of a red flower. Surprisingly, the expression of anthocyanin structural and potential regulatory genes in white flowers was not significantly lower than that in both pink and red flowers. However, proteomic analysis revealed a GST encoded by a gene-regulator involved in anthocyanin transport (Riant)-which is expressed in the red flower, but almost undetectable in the variegated flower. The Riant gene contains an insertion-deletion (indel) polymorphism in exon 3. In white flowers, the Riant gene is interrupted by a 2-bp insertion in the last exon, which causes a frameshift and a premature stop codon. In contrast, both pink and red flowers that arise from bud sports are heterozygous for the Riant locus, with one functional allele due to the 2-bp deletion or a novel 1-bp insertion. Southern blot analysis indicated that the Riant gene occurs in a single copy in the peach genome and it is not interrupted by a transposon. The function of the Riant gene was confirmed by its ectopic expression in the Arabidopsis tt19 mutant, where it complements the anthocyanin phenotype, but not the proanthocyanidin pigmentation in seed coat. Collectively,these results indicate that a small indel mutation in the Riant gene, which is not the result of a transposon insertion or excision, causes variegated colouration of peach flowers.

Genome-wide insertion-deletion (InDel) marker discovery and genotyping for genomics-assisted breeding applications in chickpea

We developed 21,499 genome-wide insertion-deletion (InDel) markers (2- to 54-bp in silico fragment length polymorphism) by comparing the genomic sequences of four (desi, kabuli and wild C. reticulatum) chickpea [Cicer arietinum (L.)] accessions. InDel markers showing 2- to 6-bp fragment length polymorphism among accessions were abundant (76.8%) in the chickpea genome. The physically mapped 7,643 and 13,856 markers on eight chromosomes and unanchored scaffolds, respectively, were structurally and functionally annotated. The 4,506 coding (23% large-effect frameshift mutations) and regulatory InDel markers were identified from 3,228 genes (representing 11.7% of total 27,571 desi genes), suggesting their functional relevance for trait association/genetic mapping. High amplification (97%) and intra-specific polymorphic (60-83%) potential and wider genetic diversity (15-89%) were detected by genome-wide 6,254 InDel markers among desi, kabuli and wild accessions using even a simpler cost-effective agarose gel-based assay. This signifies added advantages of this user-friendly genetic marker system for manifold large-scale genotyping applications in laboratories with limited infrastructure and resources. Utilizing 6,254 InDel markers-based high-density (inter-marker distance: 0.212 cM) inter-specific genetic linkage map (ICC 4958 × ICC 17160) of chickpea as a reference, three major genomic regions harboring six flowering and maturity time robust QTLs (16.4-27.5% phenotypic variation explained, 8.1-11.5 logarithm of odds) were identified. Integration of genetic and physical maps at these target QTL intervals mapped on three chromosomes delineated five InDel markers-containing candidate genes tightly linked to the QTLs governing flowering and maturity time in chickpea. Taken together, our study demonstrated the practical utility of developing and high-throughput genotyping of such beneficial InDel markers at a genome-wide scale to expedite genomics-assisted breeding applications in chickpea.

The complete mitochondrial genome sequence of the green microalga Lobosphaera (Parietochloris) incisa reveals a new type of palindromic repetitive repeat

Background

Lobosphaera incisa, formerly known as Myrmecia incisa and then Parietochloris incisa, is an oleaginous unicellular green alga belonging to the class Trebouxiophyceae (Chlorophyta). It is the richest known plant source of arachidonic acid, an ω-6 poly-unsaturated fatty acid valued by the pharmaceutical and baby-food industries. It is therefore an organism of high biotechnological interest, and we recently reported the sequence of its chloroplast genome.

Results

We now report the complete sequence of the mitochondrial genome of L. incisa from high-throughput Illumina short-read sequencing. The circular chromosome of 69,997 bp is predicted to encode a total of 64 genes, some harboring specific self-splicing group I and group II introns. Overall, the gene content is highly similar to that of the mitochondrial genomes of other Trebouxiophyceae, with 34 protein-coding, 3 rRNA, and 27 tRNA genes. Genes are distributed in two clusters located on different DNA strands, a bipartite arrangement that suggests expression from two divergent promoters yielding polycistronic primary transcripts. The L. incisa mitochondrial genome contains families of intergenic dispersed DNA repeat sequences that are not shared with other known mitochondrial genomes of Trebouxiophyceae. The most peculiar feature of the genome is a repetitive palindromic repeat, the LIMP (L. Incisa Mitochondrial Palindrome), found 19 times in the genome. It is formed by repetitions of an AACCA pentanucleotide, followed by an invariant 7-nt loop and a complementary repeat of the TGGTT motif. Analysis of the genome sequencing reads indicates that the LIMP can be a substrate for large-scale genomic rearrangements. We speculate that LIMPs can act as origins of replication. Deep sequencing of the L. incisa transcriptome also suggests that the LIMPs with long stems are sites of transcript processing. The genome also contains five copies of a related palindromic repeat, the HyLIMP, with a 10-nt motif related to that of the LIMP.

Conclusions

The mitochondrial genome of L. incisa encodes a unique type of repetitive palindromic repeat sequence, the LIMP, which can mediate genome rearrangements and play a role in mitochondrial gene expression. Experimental studies are needed to confirm and further characterize the functional role(s) of the LIMP.

DNA polymerase from temperate phage Bam35 is endowed with processive polymerization and abasic sites translesion synthesis capacity

DNA polymerases (DNAPs) responsible for genome replication are highly faithful enzymes that nonetheless cannot deal with damaged DNA. In contrast, translesion synthesis (TLS) DNAPs are suitable for replicating modified template bases, although resulting in very low-fidelity products. Here we report the biochemical characterization of the temperate bacteriophage Bam35 DNA polymerase (B35DNAP), which belongs to the protein-primed subgroup of family B DNAPs, along with phage Φ29 and other viral and mobile element polymerases. B35DNAP is a highly faithful DNAP that can couple strand displacement to processive DNA synthesis. These properties allow it to perform multiple displacement amplification of plasmid DNA with a very low error rate. Despite its fidelity and proofreading activity, B35DNAP was able to successfully perform abasic site TLS without template realignment and inserting preferably an A opposite the abasic site (A rule). Moreover, deletion of the TPR2 subdomain, required for processivity, impaired primer extension beyond the abasic site. Taken together, these findings suggest that B35DNAP may perform faithful and processive genome replication in vivo and, when required, TLS of abasic sites.

Differences in genome-wide repeat sequence instability conferred by proofreading and mismatch repair defects

Mutation rates are used to calibrate molecular clocks and to link genetic variants with human disease. However, mutation rates are not uniform across each eukaryotic genome. Rates for insertion/deletion (indel) mutations have been found to vary widely when examined in vitro and at specific loci in vivo. Here, we report the genome-wide rates of formation and repair of indels made during replication of yeast nuclear DNA. Using over 6000 indels accumulated in four mismatch repair (MMR) defective strains, and statistical corrections for false negatives, we find that indel rates increase by 100 000-fold with increasing homonucleotide run length, representing the greatest effect on replication fidelity of any known genomic parameter. Nonetheless, long genomic homopolymer runs are overrepresented relative to random chance, implying positive selection. Proofreading defects in the replicative polymerases selectively increase indel rates in short repetitive tracts, likely reflecting the distance over which Pols δ and ϵ interact with duplex DNA upstream of the polymerase active site. In contrast, MMR defects hugely increase indel mutagenesis in long repetitive sequences. Because repetitive sequences are not uniformly distributed among genomic functional elements, the quantitatively different consequences on genome-wide repeat sequence instability conferred by defects in proofreading and MMR have important biological implications.

Noncanonical secondary structure stabilizes mitochondrial tRNA(Ser(UCN)) by reducing the entropic cost of tertiary folding

Mammalian mitochondrial tRNA(Ser(UCN)) (mt-tRNA(Ser)) and pyrrolysine tRNA (tRNA(Pyl)) fold to near-canonical three-dimensional structures despite having noncanonical secondary structures with shortened interhelical loops that disrupt the conserved tRNA tertiary interaction network. How these noncanonical tRNAs compensate for their loss of tertiary interactions remains unclear. Furthermore, in human mt-tRNA(Ser), lengthening the variable loop by the 7472insC mutation reduces mt-tRNA(Ser) concentration in vivo through poorly understood mechanisms and is strongly associated with diseases such as deafness and epilepsy. Using simulations of the TOPRNA coarse-grained model, we show that increased topological constraints encoded by the unique secondary structure of wild-type mt-tRNA(Ser) decrease the entropic cost of folding by ∼2.5 kcal/mol compared to canonical tRNA, offsetting its loss of tertiary interactions. Further simulations show that the pathogenic 7472insC mutation disrupts topological constraints and hence destabilizes the mutant mt-tRNA(Ser) by ∼0.6 kcal/mol relative to wild-type. UV melting experiments confirm that insertion mutations lower mt-tRNA(Ser) melting temperature by 6-9 °C and increase the folding free energy by 0.8-1.7 kcal/mol in a largely sequence- and salt-independent manner, in quantitative agreement with our simulation predictions. Our results show that topological constraints provide a quantitative framework for describing key aspects of RNA folding behavior and also provide the first evidence of a pathogenic mutation that is due to disruption of topological constraints.

Fidelity consequences of the impaired interaction between DNA polymerase epsilon and the GINS complex

DNA polymerase epsilon interacts with the CMG (Cdc45-MCM-GINS) complex by Dpb2p, the non-catalytic subunit of DNA polymerase epsilon. It is postulated that CMG is responsible for targeting of Pol ɛ to the leading strand. We isolated a mutator dpb2-100 allele which encodes the mutant form of Dpb2p. We showed previously that Dpb2-100p has impaired interactions with Pol2p, the catalytic subunit of Pol ɛ. Here, we present that Dpb2-100p has strongly impaired interaction with the Psf1 and Psf3 subunits of the GINS complex. Our in vitro results suggest that while dpb2-100 does not alter Pol ɛ's biochemical properties including catalytic efficiency, processivity or proofreading activity - it moderately decreases the fidelity of DNA synthesis. As the in vitro results did not explain the strong in vivo mutator effect of the dpb2-100 allele we analyzed the mutation spectrum in vivo. The analysis of the mutation rates in the dpb2-100 mutant indicated an increased participation of the error-prone DNA polymerase zeta in replication. However, even in the absence of Pol ζ activity the presence of the dpb2-100 allele was mutagenic, indicating that a significant part of mutagenesis is Pol ζ-independent. A strong synergistic mutator effect observed for transversions in the triple mutant dpb2-100 pol2-4 rev3Δ as compared to pol2-4 rev3Δ and dpb2-100 rev3Δ suggests that in the presence of the dpb2-100 allele the number of replication errors is enhanced. We hypothesize that in the dpb2-100 strain, where the interaction between Pol ɛ and GINS is weakened, the access of Pol δ to the leading strand may be increased. The increased participation of Pol δ on the leading strand in the dpb2-100 mutant may explain the synergistic mutator effect observed in the dpb2-100 pol3-5DV double mutant.

The pattern of DNA cleavage intensity around indels

Indels (insertions and deletions) are the second most common form of genetic variations in the eukaryotic genomes and are responsible for a multitude of genetic diseases. Despite its significance, detailed molecular mechanisms for indel generation are still unclear. Here we examined 2,656,597 small human and mouse germline indels, 16,742 human somatic indels, 10,599 large human insertions, and 5,822 large chimpanzee insertions and systematically analyzed the patterns of DNA cleavage intensities in the 200 base pair regions surrounding these indels. Our results show that DNA cleavage intensities close to the start and end points of indels are significantly lower than other regions, for both small human germline and somatic indels and also for mouse small indels. Compared to small indels, the patterns of DNA cleavage intensity around large indels are more complex, and there are two low intensity regions near each end of the indels that are approximately 13 bp apart from each other. Detailed analyses of a subset of indels show that there is slight difference in cleavage intensity distribution between insertion indels and deletion indels that could be contributed by their respective enrichment of different repetitive elements. These results will provide new insight into indel generation mechanisms.