Replication slippage involves DNA polymerase pausing and dissociation

Quantifying cell divisions along evolutionary lineages in cancer

Cell division drives somatic evolution but is challenging to quantify. We developed a framework to count cell divisions with DNA replication-related mutations in polyguanine homopolymers. Analyzing 505 samples from 37 patients, we studied the milestones of colorectal cancer evolution. Primary tumors diversify at ~250 divisions from the founder cell, while distant metastasis divergence occurs significantly later, at ~500 divisions. Notably, distant but not lymph node metastases originate from primary tumor regions that have undergone surplus divisions, tying subclonal expansion to metastatic capacity. Then, we analyzed a cohort of 73 multifocal lung cancers and showed that the cell division burden of the tumors' common ancestor distinguishes independent primary tumors from intrapulmonary metastases and correlates with patient survival. In lung cancer too, metastatic capacity is tied to more extensive proliferation. The cell division history of human cancers is easily accessible using our simple framework and contains valuable biological and clinical information.

Fitness landscapes of human microsatellites

Advances in DNA sequencing technology and computation now enable genome-wide scans for natural selection to be conducted on unprecedented scales. By examining patterns of sequence variation among individuals, biologists are identifying genes and variants that affect fitness. Despite this progress, most population genetic methods for characterizing selection assume that variants mutate in a simple manner and at a low rate. Because these assumptions are violated by repetitive sequences, selection remains uncharacterized for an appreciable percentage of the genome. To meet this challenge, we focus on microsatellites, repetitive variants that mutate orders of magnitude faster than single nucleotide variants, can harbor substantial variation, and are known to influence biological function in some cases. We introduce four general models of natural selection that are each characterized by just two parameters, are easily simulated, and are specifically designed for microsatellites. Using a random forests approach to approximate Bayesian computation, we fit these models to carefully chosen microsatellites genotyped in 200 humans from a diverse collection of eight populations. Altogether, we reconstruct detailed fitness landscapes for 43 microsatellites we classify as targets of selection. Microsatellite fitness surfaces are diverse, including a range of selection strengths, contributions from dominance, and variation in the number and size of optimal alleles. Microsatellites that are subject to selection include loci known to cause trinucleotide expansion disorders and modulate gene expression, as well as intergenic loci with no obvious function. The heterogeneity in fitness landscapes we report suggests that genome-scale analyses like those used to assess selection targeting single nucleotide variants run the risk of oversimplifying the evolutionary dynamics of microsatellites. Moreover, our fitness landscapes provide a valuable visualization of the selective dynamics navigated by microsatellites.

Zebrafish: unraveling genetic complexity through duplicated genes

The zebrafish is an invaluable model organism for genetic, developmental, and disease research. Although its high conservation with humans is often cited as justification for its use, the zebrafish harbors oft-ignored genetic characteristics that may provide unique insights into gene structure and function. Zebrafish, along with other teleost fish, underwent an additional round of whole genome duplication after their split from tetrapods-resulting in an abundance of duplicated genes when compared to other vertebrates. These duplicated genes have evolved in distinct ways over the ensuing 350 million years. Thus, each gene within a duplicated gene pair has nuanced differences that create a unique identity. By investigating both members of the gene pair together, we can elucidate the mechanisms that underly protein structure and function and drive the complex interplay within biological systems, such as signal transduction cascades, genetic regulatory networks, and evolution of tissue and organ function. It is crucial to leverage such studies to explore these molecular dynamics, which could have far-reaching implications for both basic science and therapeutic development. Here, we will review the role of gene duplications and the existing models for gene divergence and retention following these events. We will also highlight examples within each of these models where studies comparing duplicated genes in the zebrafish have yielded key insights into protein structure, function, and regulation.

Developmental mosaicism underlying EGFR-mutant lung cancer presenting with multiple primary tumors

Although the development of multiple primary tumors in smokers with lung cancer can be attributed to carcinogen-induced field cancerization, the occurrence of multiple tumors at presentation in individuals with EGFR-mutant lung cancer who lack known environmental exposures remains unexplained. In the present study, we identified ten patients with early stage, resectable, non-small cell lung cancer who presented with multiple, anatomically distinct, EGFR-mutant tumors. We analyzed the phylogenetic relationships among multiple tumors from each patient using whole-exome sequencing (WES) and hypermutable poly(guanine) (poly(G)) repeat genotyping as orthogonal methods for lineage tracing. In four patients, developmental mosaicism, assessed by WES and poly(G) lineage tracing, indicates a common non-germline cell of origin. In two other patients, we identified germline EGFR variants, which confer moderately enhanced signaling when modeled in vitro. Thus, in addition to germline variants, developmental mosaicism defines a distinct mechanism of genetic predisposition to multiple EGFR-mutant primary tumors, with implications for their etiology and clinical management.

We are doing it wrong: Putting homology before phylogeny in cyanobacterial taxonomy

The rapid expansion of whole genome sequencing in bacterial taxonomy has revealed deep evolutionary relationships and speciation signals, but assembly methods often miss true nucleotide diversity in the ribosomal operons. Though it lacks sufficient phylogenetic signal at the species level, the 16S ribosomal RNA gene is still much used in bacterial taxonomy. In cyanobacterial taxonomy, comparisons of 16S-23S Internal Transcribed Spacer (ITS) regions are used to bridge this information gap. Although ITS rRNA region analyses are routinely being used to identify species, researchers often do not identify orthologous operons, which leads to improper comparisons. No method for delineating orthologous operon copies from paralogous ones has been established. A new method for recognizing orthologous ribosomal operons by quantifying the conserved paired nucleotides in a helical domain of the ITS, has been developed. The D1' Index quantifies differences in the ratio of pyrimidines to purines in paired nucleotide sequences of this helix. Comparing 111 operon sequences from 89 strains of Brasilonema, four orthologous operon types were identified. Plotting D1' Index values against the length of helices produced clear separation of orthologs. Most orthologous operons in this study were observed both with and without tRNA genes present. We hypothesize that genomic rearrangement, not gene duplication, is responsible for the variation among orthologs. This new method will allow cyanobacterial taxonomists to utilize ITS rRNA region data more correctly, preventing erroneous taxonomic hypotheses. Moreover, this work could assist genomicists in identifying and preserving evident sequence variability in ribosomal operons, which is an important proxy for evolution in prokaryotes.

Predictive Biomarkers and Resistance Mechanisms of Checkpoint Inhibitors in Malignant Solid Tumors

Predictive biomarkers for immune checkpoint inhibitors (ICIs) in solid tumors such as melanoma, hepatocellular carcinoma (HCC), colorectal cancer (CRC), non-small cell lung cancer (NSCLC), endometrial carcinoma, renal cell carcinoma (RCC), or urothelial carcinoma (UC) include programmed cell death ligand 1 (PD-L1) expression, tumor mutational burden (TMB), defective deoxyribonucleic acid (DNA) mismatch repair (dMMR), microsatellite instability (MSI), and the tumor microenvironment (TME). Over the past decade, several types of ICIs, including cytotoxic T-lymphocyte-associated protein 4 (CTLA-4) inhibitors, anti-programmed cell death 1 (PD-1) antibodies, anti-programmed cell death ligand 1 (PD-L1) antibodies, and anti-lymphocyte activation gene-3 (LAG-3) antibodies have been studied and approved by the Food and Drug Administration (FDA), with ongoing research on others. Recent studies highlight the critical role of the gut microbiome in influencing a positive therapeutic response to ICIs, emphasizing the importance of modeling factors that can maintain a healthy microbiome. However, resistance mechanisms can emerge, such as increased expression of alternative immune checkpoints, T-cell immunoglobulin (Ig), mucin domain-containing protein 3 (TIM-3), LAG-3, impaired antigen presentation, and alterations in the TME. This review aims to synthesize the data regarding the interactions between microbiota and immunotherapy (IT). Understanding these mechanisms is essential for optimizing ICI therapy and developing effective combination strategies.

An additional proofreader contributes to DNA replication fidelity in mycobacteria

The removal of mis-incorporated nucleotides by proofreading activity ensures DNA replication fidelity. Whereas the ε-exonuclease DnaQ is a well-established proofreader in the model organism Escherichia coli, it has been shown that proofreading in a majority of bacteria relies on the polymerase and histidinol phosphatase (PHP) domain of replicative polymerase, despite the presence of a DnaQ homolog that is structurally and functionally distinct from E. coli DnaQ. However, the biological functions of this type of noncanonical DnaQ remain unclear. Here, we provide independent evidence that noncanonical DnaQ functions as an additional proofreader for mycobacteria. Using the mutation accumulation assay in combination with whole-genome sequencing, we showed that depletion of DnaQ in Mycolicibacterium smegmatis leads to an increased mutation rate, resulting in AT-biased mutagenesis and increased insertions/deletions in the homopolymer tract. Our results showed that mycobacterial DnaQ binds to the β clamp and functions synergistically with the PHP domain proofreader to correct replication errors. Furthermore, the loss of dnaQ results in replication fork dysfunction, leading to attenuated growth and increased mutagenesis on subinhibitory fluoroquinolones potentially due to increased vulnerability to fork collapse. By analyzing the sequence polymorphism of dnaQ in clinical isolates of Mycobacterium tuberculosis (Mtb), we demonstrated that a naturally evolved DnaQ variant prevalent in Mtb lineage 4.3 may enable hypermutability and is associated with drug resistance. These results establish a coproofreading model and suggest a division of labor between DnaQ and PHP domain proofreader. This study also provides real-world evidence that a mutator-driven evolutionary pathway may exist during the adaptation of Mtb.

Quantifying Replication Slippage Error in Cryptosporidium Metabarcoding Studies

Genetic variation in Cryptosporidium, a common protozoan gut parasite in humans, is often based on marker genes containing trinucleotide repeats, which differentiate subtypes and track outbreaks. However, repeat regions have high replication slippage rates, making it difficult to discern biological diversity from error. Here, we synthesized Cryptosporidium DNA in clonal plasmid vectors, amplified them in different mock community ratios, and sequenced them using next-generation sequencing to determine the rate of replication slippage with dada2. Our results indicate that slippage rates increase with the length of the repeat region and can contribute to error rates of up to 20%.

AMPLON: Amplifying DNA with Multiarm Priming and Looping Optimization of Nucleic Acid

Nucleic acid amplification, the bedrock of biotechnology and molecular diagnostics, surges in applications-especially isothermal approaches-heightening the demand for advanced and precisely engineered methods. Here, a novel approach for amplifying DNA with multiarm priming and looping optimization of nucleic acid (AMPLON) is presented. AMPLON relies on a novel polymeric material with unique set of multiarm polyethylene glycol-DNA primers for efficient DNA amplification under isothermal conditions. Each arm carries single-stranded DNA complementing the sense or antisense sequence of the target DNA. The amplification reaction begins with antisense arms binding to the target DNA, forming a template for sense-carrying arms to direct multiarm large DNA amplicon synthesis through successive DNA looping and unlooping steps. Using human immunodeficiency virus type 1 (HIV-1) as a model clinical target, AMPLON exhibits high sensitivity, detecting target concentrations as low as 100 copies mL-1. Compared to a quantitative real-time polymerase chain reaction assay using sensitive primers, AMPLON reliably identifies HIV-1 RNA in plasma samples (n = 20) with a significant agreement rate of 95%. With its ability to achieve highly specific and sensitive target amplification within 30 min, AMPLON holds immense potential to transform the field of nucleic acid research and unleashing new possibilities in medicine and biotechnology.

Sequencing and characterizing short tandem repeats in the human genome

Short tandem repeats (STRs) are highly polymorphic sequences throughout the human genome that are composed of repeated copies of a 1-6-bp motif. Over 1 million variable STR loci are known, some of which regulate gene expression and influence complex traits, such as height. Moreover, variants in at least 60 STR loci cause genetic disorders, including Huntington disease and fragile X syndrome. Accurately identifying and genotyping STR variants is challenging, in particular mapping short reads to repetitive regions and inferring expanded repeat lengths. Recent advances in sequencing technology and computational tools for STR genotyping from sequencing data promise to help overcome this challenge and solve genetically unresolved cases and the 'missing heritability' of polygenic traits. Here, we compare STR genotyping methods, analytical tools and their applications to understand the effect of STR variation on health and disease. We identify emergent opportunities to refine genotyping and quality-control approaches as well as to integrate STRs into variant-calling workflows and large cohort analyses.

Effects of trimer repeats on Psidium guajava L. gene expression and prospection of functional microsatellite markers

Most research on trinucleotide repeats (TRs) focuses on human diseases, with few on the impact of TR expansions on plant gene expression. This work investigates TRs' effect on global gene expression in Psidium guajava L., a plant species with widespread distribution and significant relevance in the food, pharmacology, and economics sectors. We analyzed TR-containing coding sequences in 1,107 transcripts from 2,256 genes across root, shoot, young leaf, old leaf, and flower bud tissues of the Brazilian guava cultivars Cortibel RM and Paluma. Structural analysis revealed TR sequences with small repeat numbers (5-9) starting with cytosine or guanine or containing these bases. Functional annotation indicated TR-containing genes' involvement in cellular structures and processes (especially cell membranes and signal recognition), stress response, and resistance. Gene expression analysis showed significant variation, with a subset of highly expressed genes in both cultivars. Differential expression highlighted numerous down-regulated genes in Cortibel RM tissues, but not in Paluma, suggesting interplay between tissues and cultivars. Among 72 differentially expressed genes with TRs, 24 form miRNAs, 13 encode transcription factors, and 11 are associated with transposable elements. In addition, a set of 20 SSR-annotated, transcribed, and differentially expressed genes with TRs was selected as phenotypic markers for Psidium guajava and, potentially for closely related species as well.

Unravelling Chlamydia trachomatis diversity in Amhara, Ethiopia: MLVA-ompA sequencing as a molecular typing tool for trachoma

Trachoma is the leading infectious cause of blindness worldwide and is now largely confined to around 40 low- and middle-income countries. It is caused by Chlamydia trachomatis (Ct), a contagious intracellular bacterium. The World Health Organization recommends mass drug administration (MDA) with azithromycin for treatment and control of ocular Ct infections, alongside improving facial cleanliness and environmental conditions to reduce transmission. To understand the molecular epidemiology of trachoma, especially in the context of MDA and transmission dynamics, the identification of Ct genotypes could be useful. While many studies have used the Ct major outer membrane protein gene (ompA) for genotyping, it has limitations. Our study applies a typing system novel to trachoma, Multiple Loci Variable Number Tandem Repeat Analysis combined with ompA (MLVA-ompA). Ocular swabs were collected post-MDA from four trachoma-endemic zones in Ethiopia between 2011-2017. DNA from 300 children with high Ct polymerase chain reaction (PCR) loads was typed using MLVA-ompA, utilizing 3 variable number tandem repeat (VNTR) loci within the Ct genome. Results show that MLVA-ompA exhibited high discriminatory power (0.981) surpassing the recommended threshold for epidemiological studies. We identified 87 MLVA-ompA variants across 26 districts. No significant associations were found between variants and clinical signs or chlamydial load. Notably, overall Ct diversity significantly decreased after additional MDA rounds, with a higher proportion of serovar A post-MDA. Despite challenges in sequencing one VNTR locus (CT1299), MLVA-ompA demonstrated cost-effectiveness and efficiency relative to whole genome sequencing, providing valuable information for trachoma control programs on local epidemiology. The findings suggest the potential of MLVA-ompA as a reliable tool for typing ocular Ct and understanding transmission dynamics, aiding in the development of targeted interventions for trachoma control.

Structural investigation of pathogenic RFC1 AAGGG pentanucleotide repeats reveals a role of G-quadruplex in dysregulated gene expression in CANVAS

An expansion of AAGGG pentanucleotide repeats in the replication factor C subunit 1 (RFC1) gene is the genetic cause of cerebellar ataxia, neuropathy, and vestibular areflexia syndrome (CANVAS), and it also links to several other neurodegenerative diseases including the Parkinson's disease. However, the pathogenic mechanism of RFC1 AAGGG repeat expansion remains enigmatic. Here, we report that the pathogenic RFC1 AAGGG repeats form DNA and RNA parallel G-quadruplex (G4) structures that play a role in impairing biological processes. We determine the first high-resolution nuclear magnetic resonance (NMR) structure of a bimolecular parallel G4 formed by d(AAGGG)2AA and reveal how AAGGG repeats fold into a higher-order structure composed of three G-tetrad layers, and further demonstrate the formation of intramolecular G4s in longer DNA and RNA repeats. The pathogenic AAGGG repeats, but not the nonpathogenic AAAAG repeats, form G4 structures to stall DNA replication and reduce gene expression via impairing the translation process in a repeat-length-dependent manner. Our results provide an unprecedented structural basis for understanding the pathogenic mechanism of AAGGG repeat expansion associated with CANVAS. In addition, the high-resolution structures resolved in this study will facilitate rational design of small-molecule ligands and helicases targeting G4s formed by AAGGG repeats for therapeutic interventions.

Guide RNA scaffold variants enabled easy cloning of large gRNA cluster for multiplexed gene editing

Cas9 protein-mediated gene editing has revolutionized genetic manipulation in most organisms. There are many cases where multiplexed gene editing is needed. Cas9 is capable of multiplex gene editing when expressed with multiple guide RNAs. Conventional cloning methods for multiplexed gene editing vector is not efficient due to repeated use of a single-guide RNA scaffold and inefficient ligation. In this study, we conducted structure-guided mutagenesis and random mutagenesis on the original sgRNA scaffold and identified a large number of functional sgRNA scaffold variants. With these scaffold variants and different tRNAs, fusion polymerase chain reaction protocol was developed to rapidly synthesize spacer-scaffold-tRNA-spacer units with up to 9 targets. In conjunction with golden gate cloning, gene editing vectors with up to 24 target sites were efficiently cloned in one-step cloning. One such gene editing vector targeting 12 genes in tomato were tested in stable transformation and 10 out of the 12 genes were found mutated in a single transgenic line. To facilitate the application of multiplexed gene editing using these scaffold variants and tRNAs from different species, a webserver was created to generate primer sets and provide template sequences for the synthesis of large sgRNA expression units based on the user-supplied target sequences and species.

New Discoveries on Protein Recruitment and Regulation during the Early Stages of the DNA Damage Response Pathways

Maintaining genomic stability and properly repairing damaged DNA is essential to staying healthy and preserving cellular homeostasis. The five major pathways involved in repairing eukaryotic DNA include base excision repair (BER), nucleotide excision repair (NER), mismatch repair (MMR), non-homologous end joining (NHEJ), and homologous recombination (HR). When these pathways do not properly repair damaged DNA, genomic stability is compromised and can contribute to diseases such as cancer. It is essential that the causes of DNA damage and the consequent repair pathways are fully understood, yet the initial recruitment and regulation of DNA damage response proteins remains unclear. In this review, the causes of DNA damage, the various mechanisms of DNA damage repair, and the current research regarding the early steps of each major pathway were investigated.

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.

Ultra-high static magnetic field induces a change in the spectrum but not frequency of DNA spontaneous mutations in Arabidopsis thaliana

Biological effects of magnetic fields have been extensively studied in plants, microorganisms and animals, and applications of magnetic fields in regulation of plant growth and phytoprotection is a promising field in sustainable agriculture. However, the effect of magnetic fields especially ultra-high static magnetic field (UHSMF) on genomic stability is largely unclear. Here, we investigated the mutagenicity of 24.5, 30.5 and 33.0 T UHSMFs with the gradient of 150, 95 and 0 T/m, respectively, via whole genome sequencing. Our results showed that 1 h exposure of Arabidopsis dried seeds to UHSMFs has no significant effect on the average rate of DNA mutations including single nucleotide variations and InDels (insertions and deletions) in comparison with the control, but 33.0 T and 24.5 T treatments lead to a significant change in the rate of nucleotide transitions and InDels longer than 3 bp, respectively, suggesting that both strength and gradient of UHSMF impact molecular spectrum of DNA mutations. We also found that the decreased transition rate in UHSMF groups is correlated with the upstream flanking sequences of G and C mutation sites. Furthermore, the germination rate of seeds exposed to 24.5 T SMF with -150 T/m gradient showed a significant decrease at 24 hours after sowing. Overall, our data lay a basis for precisely assessing the potential risk of UHSMF on DNA stability, and for elucidating molecular mechanism underlying gradient SMF-regulated biological processes in the future.

A deleterious frameshift insertion mutation in the ZNF142 gene leads to intellectual developmental disorder with impaired speech in three affected siblings: Clinical features and literature review

BACKGROUND

ZNF142 gene is a protein-coding gene encoding Zinc Finger Protein 142. ZNF proteins are a vast group of cellular effectors with a wide range of functions such as signal transduction, transcriptional regulation, meiotic recombination, DNA repair, development, and cell migration. Mutations in the ZNF142 gene are related to neurodevelopmental disorder with impaired speech and hyperkinetic movements (NEDISHM). This study on a family with three affected siblings identified a pathogenic frameshift insertion variant. In addition, we conducted a review of the literature on previously reported ZNF142 gene variants and their clinical manifestations.

MATERIALS AND METHODS

Three affected siblings with severe intellectual developmental disabilities and speech impairments, their parents, and other sibs in the family were included. The patients were studied by the whole exome sequencing. Sanger sequencing, co-segregation analysis, and in silico analysis were carried out to verify candidate variant. The identified variant was interpreted based on the ACMG guideline.

RESULTS

We identified a frameshift insertion variant in the ZNF142 gene, NM_001379659.1: c.3755dup (NP_001366588.1:p.Arg1253ThrfsTer15), that was related to the clinical features of three patients. The identified variant was found to be pathogenic.

CONCLUSION

The current study findings expand the existing knowledge of the variant on the ZNF142 gene implicated in the neurodevelopmental disorder, intellectual disability, and impaired speech and it presents a detailed clinical feature associated with related conditions. The data have implications for genetic diagnosis and counseling in families with the same disorders.

Uncovering the Genetic and Molecular Features of Huntington's Disease in Northern Colombia

Huntington's disease (HD) is a genetic disorder caused by a CAG trinucleotide expansion in the huntingtin (HTT) gene. Juan de Acosta, Atlántico, a city located on the Caribbean coast of Colombia, is home to the world's second-largest HD pedigree. Here, we include 291 descendants of this pedigree with at least one family member with HD. Blood samples were collected, and genomic DNA was extracted. We quantified the HTT CAG expansion using an amplicon sequencing protocol. The genetic heterogeneity was measured as the ratio of the mosaicism allele's read peak and the slippage ratio of the allele's read peak from our sequence data. The statistical and bioinformatic analyses were performed with a significance threshold of p < 0.05. We found that the average HTT CAG repeat length in all participants was 21.91 (SD = 8.92). Of the 291 participants, 33 (11.3%, 18 females) had a positive molecular diagnosis for HD. Most affected individuals were adults, and the most common primary and secondary alleles were 17/7 (CAG/CCG) and 17/10 (CAG/CCG), respectively. The mosaicism increased with age in the participants with HD, while the slippage analyses revealed differences by the HD allele type only for the secondary allele. The slippage tended to increase with the HTT CAG repeat length in the participants with HD, but the increase was not statistically significant. This study analyzed the genetic and molecular features of 291 participants, including 33 with HD. We found that the mosaicism increased with age in the participants with HD, particularly for the secondary allele. The most common haplotype was 17/7_17/10. The slippage for the secondary allele varied by the HD allele type, but there was no significant difference in the slippage by sex. Our findings offer valuable insights into HD and could have implications for future research and clinical management.

Poly (ADP-ribose) polymerase 1 and neurodegenerative diseases: Past, present, and future

Poly (ADP-ribose) polymerase 1 (PARP1) is a first responder that recognizes DNA damage and facilitates its repair. Neurodegenerative diseases, characterized by progressive neuron loss driven by various risk factors, including DNA damage, have increasingly shed light on the pivotal involvement of PARP1. During the early phases of neurodegenerative diseases, PARP1 experiences controlled activation to swiftly address mild DNA damage, thereby contributing to maintain brain homeostasis. However, in late stages, exacerbated PARP1 activation precipitated by severe DNA damage exacerbates the disease condition. Consequently, inhibition of PARP1 overactivation emerges as a promising therapeutic approach for neurodegenerative diseases. In this review, we comprehensively synthesize and explore the multifaceted role of PARP1 in neurodegenerative diseases, with a particular emphasis on its over-activation in the aggregation of misfolded proteins, dysfunction of the autophagy-lysosome pathway, mitochondrial dysfunction, neuroinflammation, and blood-brain barrier (BBB) injury. Additionally, we encapsulate the therapeutic applications and limitations intrinsic of PARP1 inhibitors, mainly including limited specificity, intricate pathway dynamics, constrained clinical translation, and the heterogeneity of patient cohorts. We also explore and discuss the potential synergistic implementation of these inhibitors alongside other agents targeting DNA damage cascades within neurodegenerative diseases. Simultaneously, we propose several recommendations for the utilization of PARP1 inhibitors within the realm of neurodegenerative disorders, encompassing factors like the disease-specific roles of PARP1, combinatorial therapeutic strategies, and personalized medical interventions. Lastly, the encompassing review presents a forward-looking perspective along with strategic recommendations that could guide future research endeavors in this field.

Non-B DNA structures as a booster of genome instability

Non-canonical secondary structures (NCSs) are alternative nucleic acid structures that differ from the canonical B-DNA conformation. NCSs often occur in repetitive DNA sequences and can adopt different conformations depending on the sequence. The majority of these structures form in the context of physiological processes, such as transcription-associated R-loops, G4s, as well as hairpins and slipped-strand DNA, whose formation can be dependent on DNA replication. It is therefore not surprising that NCSs play important roles in the regulation of key biological processes. In the last years, increasing published data have supported their biological role thanks to genome-wide studies and the development of bioinformatic prediction tools. Data have also highlighted the pathological role of these secondary structures. Indeed, the alteration or stabilization of NCSs can cause the impairment of transcription and DNA replication, modification in chromatin structure and DNA damage. These events lead to a wide range of recombination events, deletions, mutations and chromosomal aberrations, well-known hallmarks of genome instability which are strongly associated with human diseases. In this review, we summarize molecular processes through which NCSs trigger genome instability, with a focus on G-quadruplex, i-motif, R-loop, Z-DNA, hairpin, cruciform and multi-stranded structures known as triplexes.

Germline mutations and developmental mosaicism underlying EGFR-mutant lung cancer

While the development of multiple primary tumors in smokers with lung cancer can be attributed to carcinogen-induced field cancerization, the occurrence of multiple primary tumors in individuals with EGFR-mutant lung cancer who lack known environmental exposures remains unexplained. We identified ten patients with early-stage, resectable non-small cell lung cancer who presented with multiple anatomically distinct EGFR-mutant tumors. We analyzed the phylogenetic relationships among multiple tumors from each patient using whole exome sequencing (WES) and hypermutable poly-guanine (poly-G) repeat genotyping, as orthogonal methods for lineage tracing. In two patients, we identified germline EGFR variants, which confer moderately enhanced signaling when modeled in vitro. In four other patients, developmental mosaicism is supported by the poly-G lineage tracing and WES, indicating a common non-germline cell-of-origin. Thus, developmental mosaicism and germline variants define two distinct mechanisms of genetic predisposition to multiple EGFR-mutant primary tumors, with implications for understanding their etiology and clinical management.

The Mitochondrial Genome of the Holoparasitic Plant Thonningia sanguinea Provides Insights into the Evolution of the Multichromosomal Structure

Multichromosomal mitochondrial genome (mitogenome) structures have repeatedly evolved in many lineages of angiosperms. However, the underlying mechanism remains unclear. The mitogenomes of three genera of Balanophoraceae, namely Lophophytum, Ombrophytum, and Rhopalocnemis, have already been sequenced and assembled, all showing a highly multichromosomal structure, albeit with different genome and chromosome sizes. It is expected that characterization of additional lineages of this family may expand the knowledge of mitogenome diversity and provide insights into the evolution of the plant mitogenome structure and size. Here, we assembled and characterized the mitogenome of Thonningia sanguinea, which, together with Balanophora, forms a clade sister to the clade comprising Lophophytum, Ombrophytum, and Rhopalocnemis. The mitogenome of T. sanguinea possesses a multichromosomal structure of 18 circular chromosomes of 8.7-19.2 kb, with a total size of 246,247 bp. There are very limited shared regions and poor chromosomal correspondence between T. sanguinea and other Balanophoraceae species, suggesting frequent rearrangements and rapid sequence turnover. Numerous medium- and small-sized repeats were identified in the T. sanguinea mitogenome; however, no repeat-mediated recombination was detected, which was verified by Illumina reads mapping and PCR experiments. Intraspecific mitogenome variations in T. sanguinea are mostly insertions and deletions, some of which can lead to degradation of perfect repeats in one or two accessions. Based on the mitogenome features of T. sanguinea, we propose a mechanism to explain the evolution of a multichromosomal mitogenome from a master circle, which involves mutation in organellar DNA replication, recombination and repair genes, decrease of recombination, and repeat degradation.

High-resolution mapping reveals the mechanism and contribution of genome insertions and deletions to RNA virus evolution

RNA viruses rapidly adapt to selective conditions due to the high intrinsic mutation rates of their RNA-dependent RNA polymerases (RdRps). Insertions and deletions (indels) in viral genomes are major contributors to both deleterious mutational load and evolutionary novelty, but remain understudied. To characterize the mechanistic details of their formation and evolutionary dynamics during infection, we developed a hybrid experimental-bioinformatic approach. This approach, called MultiMatch, extracts insertions and deletions from ultradeep sequencing experiments, including those occurring at extremely low frequencies, allowing us to map their genomic distribution and quantify the rates at which they occur. Mapping indel mutations in adapting poliovirus and dengue virus populations, we determine the rates of indel generation and identify mechanistic and functional constraints shaping indel diversity. Using poliovirus RdRp variants of distinct fidelity and genome recombination rates, we demonstrate tradeoffs between fidelity and Indel generation. Additionally, we show that maintaining translation frame and viral RNA structures constrain the Indel landscape and that, due to these significant fitness effects, Indels exert a significant deleterious load on adapting viral populations. Conversely, we uncover positively selected Indels that modulate RNA structure, generate protein variants, and produce defective interfering genomes in viral populations. Together, our analyses establish the kinetic and mechanistic tradeoffs between misincorporation, recombination, and Indel rates and reveal functional principles defining the central role of Indels in virus evolution, emergence, and the regulation of viral infection.

RT-based Sanger sequencing of RNAs containing complex RNA repetitive elements

Although next-generation sequencing (NGS) technologies have revolutionized our ability to sequence DNA with high-throughput, the chain termination-based Sanger sequencing method remains a widely used approach for DNA sequence analysis due to its simplicity, low cost and high accuracy. In particular, high accuracy makes Sanger sequencing the "gold standard" for sequence validation in basic research and clinical applications. During the early days of Sanger sequencing development, reverse transcriptase (RT)-based RNA sequencing was also explored and showed great promise, but the approach did not acquire popularity over time due to the limited processivity and low template unwinding capability of Avian Myeloblastosis Virus (AMV) RT, and other RT enzymes available at the time. RNA molecules have complex features, often containing repetitive sequences and stable secondary or tertiary structures. While these features are required for RNA biological function, they represent strong obstacles for retroviral RTs. Repetitive sequences and stable structures cause reverse transcription errors and premature primer extension stops, making chain termination-based methods unfeasible. MarathonRT is an ultra-processive RT encoded group II intron that can copy RNA molecules of any sequence and structure in a single cycle, making it an ideal RT enzyme for Sanger RNA sequencing. In this chapter, we upgrade the Sanger RNA sequencing method by replacing AMV RT with MarathonRT, providing a simple, robust method for direct RNA sequence analysis. The guidance for troubleshooting and further optimization are also provided.

Evolutionary Study of Protein Short Tandem Repeats in Protein Families

Tandem repeats in proteins are patterns of residues repeated directly adjacent to each other. The evolution of these repeats can be assessed by using groups of homologous sequences, which can help pointing to events of unit duplication or deletion. High pressure in a protein family for variation of a given type of repeat might point to their function. Here, we propose the analysis of protein families to calculate protein short tandem repeats (pSTRs) in each protein sequence and assess their variability within the family in terms of number of units. To facilitate this analysis, we developed the pSTR tool, a method to analyze the evolution of protein short tandem repeats in a given protein family by pairwise comparisons between evolutionarily related protein sequences. We evaluated pSTR unit number variation in protein families of 12 complete metazoan proteomes. We hypothesize that families with more dynamic ensembles of repeats could reflect particular roles of these repeats in processes that require more adaptability.

Slipped-strand mispairing within a polycytidine tract in transcriptional regulator mga leads to M protein phase variation and Mga length polymorphism in Group A Streptococcus

The M protein, a major virulence factor of Group A Streptococcus (GAS), is regulated by the multigene regulator Mga. An unexplained phenomena frequently occurring with in vitro genetic manipulation or culturing of M1T1 GAS strains is the loss of M protein production. This study was aimed at elucidating the basis for the loss of M protein production. The majority of M protein-negative (M^-) variants had one C deletion at a tract of 8 cytidines starting at base 1,571 of the M1 mga gene, which is designated as c.1571C[8]. The C deletion led to a c.1571C[7] mga variant that has an open reading frame shift and encodes a Mga-M protein fusion protein. Transformation with a plasmid containing wild-type mga restored the production of the M protein in the c.1571C[7] mga variant. Isolates producing M protein (M⁺) were recovered following growth of the c.1571C[7] M protein-negative variant subcutaneously in mice. The majority of the recovered isolates with reestablished M protein production had reverted back from c.1571C[7] to c.1571C[8] tract and some M⁺ isolates lost another C in the c.1571C[7] tract, leading to a c.1571C[6] variant that encodes a functional Mga with 13 extra amino acid residues at the C-terminus compared with wild-type Mga. The nonfunctional c.1571C[7] and functional c.1571C[6] variants are present in M1, M12, M14, and M23 strains in NCBI genome databases, and a G-to-A nonsense mutation at base 1,657 of M12 c.1574C[7] mga leads to a functional c.1574C[7]/1657A mga variant and is common in clinical M12 isolates. The numbers of the C repeats in this polycytidine tract and the polymorphism at base 1,657 lead to polymorphism in the size of Mga among clinical isolates. These findings demonstrate the slipped-strand mispairing within the c.1574C[8] tract of mga as a reversible switch controlling M protein production phase variation in multiple GAS common M types.

Origin of minicircular mitochondrial genomes in red algae

Eukaryotic organelle genomes are generally of conserved size and gene content within phylogenetic groups. However, significant variation in genome structure may occur. Here, we report that the Stylonematophyceae red algae contain multipartite circular mitochondrial genomes (i.e., minicircles) which encode one or two genes bounded by a specific cassette and a conserved constant region. These minicircles are visualized using fluorescence microscope and scanning electron microscope, proving the circularity. Mitochondrial gene sets are reduced in these highly divergent mitogenomes. Newly generated chromosome-level nuclear genome assembly of Rhodosorus marinus reveals that most mitochondrial ribosomal subunit genes are transferred to the nuclear genome. Hetero-concatemers that resulted from recombination between minicircles and unique gene inventory that is responsible for mitochondrial genome stability may explain how the transition from typical mitochondrial genome to minicircles occurs. Our results offer inspiration on minicircular organelle genome formation and highlight an extreme case of mitochondrial gene inventory reduction.

HTLV-1 bZIP factor impairs DNA mismatch repair system

Adult T-cell leukemia (ATL) is a peripheral T-cell malignancy caused by human T-cell leukemia virus type 1 (HTLV-1). Microsatellite instability (MSI) has been observed in ATL cells. Although MSI results from impaired mismatch repair (MMR) pathway, no null mutations in the genes encoding MMR factors are detectable in ATL cells. Thus, it is unclear whether or not impairment of MMR causes the MSI in ATL cells. HTLV-1 bZIP factor (HBZ) protein interacts with numerous host transcription factors and significantly contributes to disease pathogenesis and progression. Here we investigated the effect of HBZ on MMR in normal cells. The ectopic expression of HBZ in MMR-proficient cells induced MSI, and also suppressed the expression of several MMR factors. We then hypothesized that the HBZ compromises MMR by interfering with a transcription factor, nuclear respiratory factor 1 (NRF-1), and identified the consensus NRF-1 binding site at the promoter of the gene encoding MutS homologue 2 (MSH2), an essential MMR factor. The luciferase reporter assay revealed that NRF-1 overexpression enhanced MSH2 promoter activity, while co-expression of HBZ reversed this enhancement. These results supported the idea that HBZ suppresses the transcription of MSH2 by inhibiting NRF-1. Our data demonstrate that HBZ causes impaired MMR, and may imply a novel oncogenesis driven by HTLV-1.

Shadow enhancers mediate trade-offs between transcriptional noise and fidelity

Enhancers are stretches of regulatory DNA that bind transcription factors (TFs) and regulate the expression of a target gene. Shadow enhancers are two or more enhancers that regulate the same target gene in space and time and are associated with most animal developmental genes. These multi-enhancer systems can drive more consistent transcription than single enhancer systems. Nevertheless, it remains unclear why shadow enhancer TF binding sites are distributed across multiple enhancers rather than within a single large enhancer. Here, we use a computational approach to study systems with varying numbers of TF binding sites and enhancers. We employ chemical reaction networks with stochastic dynamics to determine the trends in transcriptional noise and fidelity, two key performance objectives of enhancers. This reveals that while additive shadow enhancers do not differ in noise and fidelity from their single enhancer counterparts, sub- and superadditive shadow enhancers have noise and fidelity trade-offs not available to single enhancers. We also use our computational approach to compare the duplication and splitting of a single enhancer as mechanisms for the generation of shadow enhancers and find that the duplication of enhancers can decrease noise and increase fidelity, although at the metabolic cost of increased RNA production. A saturation mechanism for enhancer interactions similarly improves on both of these metrics. Taken together, this work highlights that shadow enhancer systems may exist for several reasons: genetic drift or the tuning of key functions of enhancers, including transcription fidelity, noise and output.

DNA mutilation: A telltale sign of cancer inception

DNA damage is a discrepancy in its chemical structure precipitated by a multitude of factors. Most DNA damages can be repaired efficiently through diverse restorative mechanisms subjective to the type of damage. DNA-damaging agents elicit a medley of cellular retorts like cell cycle arrest, followed by DNA repair mechanisms or apoptosis. An unrepaired DNA damage in a nonreplicating cell does not generally engender mutations but a similar scenario in replicating cell routes to permanent modification of genetic material shrugging to carcinogenesis. DNA mutilation can be allied to disarray in bases, debasement of backbone, or crosslinks. Base damages or backbone damages like single-strand and double-strand DNA breaks are usually produced by reactive oxygen species and ionizing radiation. This substantial DNA damage has broadly been considered to be caused by various exogenous and endogenous agents with variable rates of causality and decrees of risk, sourcing toward cancer or other diseases, necessitating furtherance in diagnostics at sequential points. The purpose of this article is to review in detail the various types of DNA damages, their contributory factors, and recent developments in their identification.

The motif composition of variable number tandem repeats impacts gene expression

Understanding the impact of DNA variation on human traits is a fundamental question in human genetics. Variable number tandem repeats (VNTRs) make up ∼3% of the human genome but are often excluded from association analysis owing to poor read mappability or divergent repeat content. Although methods exist to estimate VNTR length from short-read data, it is known that VNTRs vary in both length and repeat (motif) composition. Here, we use a repeat-pangenome graph (RPGG) constructed on 35 haplotype-resolved assemblies to detect variation in both VNTR length and repeat composition. We align population-scale data from the Genotype-Tissue Expression (GTEx) Consortium to examine how variations in sequence composition may be linked to expression, including cases independent of overall VNTR length. We find that 9422 out of 39,125 VNTRs are associated with nearby gene expression through motif variations, of which only 23.4% are accessible from length. Fine-mapping identifies 174 genes to be likely driven by variation in certain VNTR motifs and not overall length. We highlight two genes, CACNA1C and RNF213, that have expression associated with motif variation, showing the utility of RPGG analysis as a new approach for trait association in multiallelic and highly variable loci.

Challenges and Therapeutic Opportunities in the dMMR/MSI-H Colorectal Cancer Landscape

About 5 to 15% of all colorectal cancers harbor mismatch repair deficient/microsatellite instability-high status (dMMR/MSI-H) that associates with high tumor mutation burden and increased immunogenicity. As a result, and in contrast to other colorectal cancer phenotypes, a significant subset of dMMR/MSI-H cancer patients strongly benefit from immunotherapy. Yet, a large proportion of these tumors remain unresponsive to any immuno-modulating treatment. For this reason, current efforts are focused on the characterization of resistance mechanisms and the identification of predictive biomarkers to guide therapeutic decision-making. Here, we provide an overview on the new advances related to the diagnosis and definition of dMMR/MSI-H status and focus on the distinct clinical, functional, and molecular cues that associate with dMMR/MSI-H colorectal cancer. We review the development of novel predictive factors of response or resistance to immunotherapy and their potential application in the clinical setting. Finally, we discuss current and emerging strategies applied to the treatment of localized and metastatic dMMR/MSI-H colorectal tumors in the neoadjuvant and adjuvant setting.

Microsatellite Instability and Aberrant Pre-mRNA Splicing: How Intimate Is It?

Cancers that belong to the microsatellite instability (MSI) class can account for up to 15% of all cancers of the digestive tract. These cancers are characterized by inactivation, through the mutation or epigenetic silencing of one or several genes from the DNA MisMatch Repair (MMR) machinery, including MLH1, MLH3, MSH2, MSH3, MSH6, PMS1, PMS2 and Exo1. The unrepaired DNA replication errors turn into mutations at several thousand sites that contain repetitive sequences, mainly mono- or dinucleotides, and some of them are related to Lynch syndrome, a predisposition condition linked to a germline mutation in one of these genes. In addition, some mutations shortening the microsatellite (MS) stretch could occur in the 3'-intronic regions, i.e., in the ATM (ATM serine/threonine kinase), MRE11 (MRE11 homolog) or the HSP110 (Heat shock protein family H) genes. In these three cases, aberrant pre-mRNA splicing was observed, and it was characterized by the occurrence of selective exon skipping in mature mRNAs. Because both the ATM and MRE11 genes, which as act as players in the MNR (MRE11/NBS1 (Nibrin)/RAD50 (RAD50 double strand break repair protein) DNA damage repair system, participate in double strand breaks (DSB) repair, their frequent splicing alterations in MSI cancers lead to impaired activity. This reveals the existence of a functional link between the MMR/DSB repair systems and the pre-mRNA splicing machinery, the diverted function of which is the consequence of mutations in the MS sequences.

Tandem repeats in giant archaeal Borg elements undergo rapid evolution and create new intrinsically disordered regions in proteins

Borgs are huge, linear extrachromosomal elements associated with anaerobic methane-oxidizing archaea. Striking features of Borg genomes are pervasive tandem direct repeat (TR) regions. Here, we present six new Borg genomes and investigate the characteristics of TRs in all ten complete Borg genomes. We find that TR regions are rapidly evolving, recently formed, arise independently, and are virtually absent in host Methanoperedens genomes. Flanking partial repeats and A-enriched character constrain the TR formation mechanism. TRs can be in intergenic regions, where they might serve as regulatory RNAs, or in open reading frames (ORFs). TRs in ORFs are under very strong selective pressure, leading to perfect amino acid TRs (aaTRs) that are commonly intrinsically disordered regions. Proteins with aaTRs are often extracellular or membrane proteins, and functionally similar or homologous proteins often have aaTRs composed of the same amino acids. We propose that Borg aaTR-proteins functionally diversify Methanoperedens and all TRs are crucial for specific Borg-host associations and possibly cospeciation.

Spelling out the roles of individual nucleoporins in nuclear export of mRNA

The Nuclear Pore Complex (NPC) represents a critical passage through the nuclear envelope for nuclear import and export that impacts nearly every cellular process at some level. Recent technological advances in the form of Auxin Inducible Degron (AID) strategies and Single-Point Edge-Excitation sub-Diffraction (SPEED) microscopy have enabled us to provide new insight into the distinct functions and roles of nuclear basket nucleoporins (Nups) upon nuclear docking and export for mRNAs. In this paper, we provide a review of our recent findings as well as an assessment of new techniques, updated models, and future perspectives in the studies of mRNA's nuclear export.

Tail Wags Dog’s SINE: Retropositional Mechanisms of Can SINE Depend on Its A-Tail Structure

Simple Summary

The genomes of higher organisms including humans are invaded by millions of repetitive elements (transposons), which can sometimes be deleterious or beneficial for hosts. Many aspects of the mechanisms underlying the expansion of transposons in the genomes remain unclear. Short retrotransposons (SINEs) are one of the most abundant classes of genomic repeats. Their amplification relies on two major processes: transcription and reverse transcription. Here, short retrotransposons of dogs and other canids called Can SINE were analyzed. Their amplification was extraordinarily active in the wolf and, particularly, dog breeds relative to other canids. We also studied a variation of their transcription mechanism involving the polyadenylation of transcripts. An analysis of specific signals involved in this process allowed us to conclude that Can SINEs could alternate amplification with and without polyadenylation in their evolution. Understanding the mechanisms of transposon replication can shed light on the mechanisms of genome function.

Abstract

SINEs, non-autonomous short retrotransposons, are widespread in mammalian genomes. Their transcripts are generated by RNA polymerase III (pol III). Transcripts of certain SINEs can be polyadenylated, which requires polyadenylation and pol III termination signals in their sequences. Our sequence analysis divided Can SINEs in canids into four subfamilies, older a1 and a2 and younger b1 and b2. Can_b2 and to a lesser extent Can_b1 remained retrotranspositionally active, while the amplification of Can_a1 and Can_a2 ceased long ago. An extraordinarily high Can amplification was revealed in different dog breeds. Functional polyadenylation signals were analyzed in Can subfamilies, particularly in fractions of recently amplified, i.e., active copies. The transcription of various Can constructs transfected into HeLa cells proposed AATAAA and (TC)_n as functional polyadenylation signals. Our analysis indicates that older Can subfamilies (a1, a2, and b1) with an active transcription terminator were amplified by the T⁺ mechanism (with polyadenylation of pol III transcripts). In the currently active Can_b2 subfamily, the amplification mechanisms with (T⁺) and without the polyadenylation of pol III transcripts (T⁻) irregularly alternate. The active transcription terminator tends to shorten, which renders it nonfunctional and favors a switch to the T⁻ retrotransposition. The activity of a truncated terminator is occasionally restored by its elongation, which rehabilitates the T⁺ retrotransposition for a particular SINE copy.

Mismatch repair is a double-edged sword in the battle against microsatellite instability

Roughly 3% of the human genome consists of microsatellites or tracts of short tandem repeats (STRs). These STRs are often unstable, undergoing high-frequency expansions (increases) or contractions (decreases) in the number of repeat units. Some microsatellite instability (MSI) is seen at multiple STRs within a single cell and is associated with certain types of cancer. A second form of MSI is characterised by expansion of a single gene-specific STR and such expansions are responsible for a group of 40+ human genetic disorders known as the repeat expansion diseases (REDs). While the mismatch repair (MMR) pathway prevents genome-wide MSI, emerging evidence suggests that some MMR factors are directly involved in generating expansions in the REDs. Thus, MMR suppresses some forms of expansion while some MMR factors promote expansion in other contexts. This review will cover what is known about the paradoxical effect of MMR on microsatellite expansion in mammalian cells.

Non-canonical DNA/RNA structures associated with the pathogenesis of Fragile X-associated tremor/ataxia syndrome and Fragile X syndrome

Fragile X-associated tremor/ataxia syndrome (FXTAS) and fragile X syndrome (FXS) are primary examples of fragile X-related disorders (FXDs) caused by abnormal expansion of CGG repeats above a certain threshold in the 5'-untranslated region of the fragile X mental retardation (FMR1) gene. Both diseases have distinct clinical manifestations and molecular pathogenesis. FXTAS is a late-adult-onset neurodegenerative disorder caused by a premutation (PM) allele (CGG expansion of 55-200 repeats), resulting in FMR1 gene hyperexpression. On the other hand, FXS is a neurodevelopmental disorder that results from a full mutation (FM) allele (CGG expansions of ≥200 repeats) leading to heterochromatization and transcriptional silencing of the FMR1 gene. The main challenge is to determine how CGG repeat expansion affects the fundamentally distinct nature of FMR1 expression in FM and PM ranges. Abnormal CGG repeat expansions form a variety of non-canonical DNA and RNA structures that can disrupt various cellular processes and cause distinct effects in PM and FM alleles. Here, we review these structures and how they are related to underlying mutations and disease pathology in FXS and FXTAS. Finally, as new CGG expansions within the genome have been identified, it will be interesting to determine their implications in disease pathology and treatment.

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.

Increased Frequency of Indels in Hypervariable Regions of SARS-CoV-2 Proteins—A Possible Signature of Adaptive Selection

Most attention in the surveillance of evolving SARS-CoV-2 genome has been centered on nucleotide substitutions in the spike glycoprotein. We show that, as the pandemic extends into its second year, the numbers and ratio of genomes with in-frame insertions and deletions (indels) increases significantly, especially among the variants of concern (VOCs). Monitoring of the SARS-CoV-2 genome evolution shows that co-occurrence (i.e., highly correlated presence) of indels, especially deletions on spike N-terminal domain and non-structural protein 6 (NSP6) is a shared feature in several VOCs such as Alpha, Beta, Delta, and Omicron. Indels distribution is correlated with spike mutations associated with immune escape and growth in the number of genomes with indels coincides with the increasing population resistance due to vaccination and previous infections. Indels occur most frequently in the spike, but also in other proteins, especially those involved in interactions with the host immune system. We also showed that indels concentrate in regions of individual SARS-CoV-2 proteins known as hypervariable regions (HVRs) that are mostly located in specific loop regions. Structural analysis suggests that indels remodel viral proteins’ surfaces at common epitopes and interaction interfaces, affecting the virus’ interactions with host proteins. We hypothesize that the increased frequency of indels, the non-random distribution of them and their independent co-occurrence in several VOCs is another mechanism of response to elevated global population immunity.

Group II intron and repeat-rich red algal mitochondrial genomes demonstrate the dynamic recent history of autocatalytic RNAs

Background

Group II introns are mobile genetic elements that can insert at specific target sequences, however, their origins are often challenging to reconstruct because of rapid sequence decay following invasion and spread into different sites. To advance understanding of group II intron spread, we studied the intron-rich mitochondrial genome (mitogenome) in the unicellular red alga, Porphyridium.

Results

Analysis of mitogenomes in three closely related species in this genus revealed they were 3–6-fold larger in size (56–132 kbp) than in other red algae, that have genomes of size 21–43 kbp. This discrepancy is explained by two factors, group II intron invasion and expansion of repeated sequences in large intergenic regions. Phylogenetic analysis demonstrates that many mitogenome group II intron families are specific to Porphyridium, whereas others are closely related to sequences in fungi and in the red alga-derived plastids of stramenopiles. Network analysis of intron-encoded proteins (IEPs) shows a clear link between plastid and mitochondrial IEPs in distantly related species, with both groups associated with prokaryotic sequences.

Conclusion

Our analysis of group II introns in Porphyridium mitogenomes demonstrates the dynamic nature of group II intron evolution, strongly supports the lateral movement of group II introns among diverse eukaryotes, and reveals their ability to proliferate, once integrated in mitochondrial DNA.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01200-3.

OUP accepted manuscript

Conflicts between transcription and replication machinery are a potent source of replication stress and genome instability; however, no technique currently exists to identify endogenous genomic locations prone to transcription–replication interactions. Here, we report a novel method to identify genomic loci prone to transcription–replication interactions termed transcription–replication immunoprecipitation on nascent DNA sequencing, TRIPn-Seq. TRIPn-Seq employs the sequential immunoprecipitation of RNA polymerase 2 phosphorylated at serine 5 (RNAP2s5) followed by enrichment of nascent DNA previously labeled with bromodeoxyuridine. Using TRIPn-Seq, we mapped 1009 unique transcription–replication interactions (TRIs) in mouse primary B cells characterized by a bimodal pattern of RNAP2s5, bidirectional transcription, an enrichment of RNA:DNA hybrids, and a high probability of forming G-quadruplexes. TRIs are highly enriched at transcription start sites and map to early replicating regions. TRIs exhibit enhanced Replication Protein A association and TRI-associated genes exhibit higher replication fork termination than control transcription start sites, two marks of replication stress. TRIs colocalize with double-strand DNA breaks, are enriched for deletions, and accumulate mutations in tumors. We propose that replication stress at TRIs induces mutations potentially contributing to age-related disease, as well as tumor formation and development.

The minicircular and extremely heteroplasmic mitogenome of the holoparasitic plant Rhopalocnemis phalloides

The plastid and nuclear genomes of parasitic plants exhibit deeply altered architectures,^1-13 whereas the few examined mitogenomes range from deeply altered to conventional.^14-20 To provide further insight on mitogenome evolution in parasitic plants, we report the highly modified mitogenome of Rhopalocnemis phalloides, a holoparasite in Balanophoraceae. Its mitogenome is uniquely arranged in 21 minicircular chromosomes that vary in size from 4,949 to 7,861 bp, with a total length of only 130,713 bp. All chromosomes share an identical 896 bp conserved region, with a large stem-loop that acts as the origin of replication, flanked on each side by hypervariable and semi-conserved regions. Similar minicircular structures with shared and unique regions have been observed in parasitic animals and free-living protists,^21-24 suggesting convergent structural evolution. Southern blots confirm both the minicircular structure and the replication origin of the mitochondrial chromosomes. PacBio reads provide evidence for chromosome recombination and rolling-circle replication for the R. phalloides mitogenome. Despite its small size, the mitogenome harbors a typical set of genes and introns within the unique regions of each chromosome, yet introns are the smallest among seed plants and ferns. The mitogenome also exhibits extreme heteroplasmy, predominantly involving short indels and more complex variants, many of which cause potential loss-of-function mutations for some gene copies. All heteroplasmic variants are transcribed, and functional and nonfunctional protein-coding variants are spliced and RNA edited. Our findings offer a unique perspective into how mitogenomes of parasitic plants can be deeply altered and shed light on plant mitogenome replication.

The Complexity of the Ovine and Caprine Keratin-Associated Protein Genes

Sheep (Ovis aries) and goats (Capra hircus) have, for more than a millennia, been a source of fibres for human use, be it for use in clothing and furnishings, for insulation, for decorative and ceremonial purposes, or for combinations thereof. While use of these natural fibres has in some respects been superseded by the use of synthetic and plant-based fibres, increased accounting for the carbon and water footprint of these fibres is creating a re-emergence of interest in fibres derived from sheep and goats. The keratin-associated proteins (KAPs) are structural components of wool and hair fibres, where they form a matrix that cross-links with the keratin intermediate filaments (KIFs), the other main structural component of the fibres. Since the first report of a complete KAP protein sequence in the late 1960s, considerable effort has been made to identify the KAP proteins and their genes in mammals, and to ascertain how these genes and proteins control fibre growth and characteristics. This effort is ongoing, with more and more being understood about the structure and function of the genes. This review consolidates that knowledge and suggests future directions for research to further our understanding.

Next-Generation Sequencing Somatic and Germline Assay Troubleshooting Guide Derived From Proficiency Testing Data

CONTEXT.—

Next-generation sequencing-based assays are increasingly used in clinical molecular laboratories to detect somatic variants in solid tumors and hematologic malignancies and to detect constitutional variants. Proficiency testing data are potential sources of information about challenges in performing these assays.

OBJECTIVE.—

To examine the most common sources of unacceptable results from the College of American Pathologists Next-Generation Sequencing Bioinformatics, Hematological Malignancies, Solid Tumor, and Germline surveys, and provide recommendations on how to avoid these pitfalls and improve performance.

DESIGN.—

The College of American Pathologists next-generation sequencing somatic and germline proficiency testing survey results from 2016 to 2019 were analyzed to identify the most common causes of unacceptable results.

RESULTS.—

On somatic and germline proficiency testing surveys, 95.9% (18 815/19 623) and 97.8% (33 890/34 641) of all variants were correctly identified, respectively. The most common causes of unacceptable results related to sequencing were false-negative errors in genomic regions that were difficult to sequence because of high GC content. False-positive errors occurred in the context of homopolymers and pseudogenes. Recurrent errors in variant annotation were seen for dinucleotide and duplication variants and included unacceptable transcript selection and outdated variant nomenclature. A small percentage of preanalytic or postanalytic errors were attributed to specimen swaps and transcription errors.

CONCLUSIONS.—

Laboratories demonstrate overall excellent performance for detecting variants in both somatic and germline proficiency testing surveys. Proficiency testing survey results highlight infrequent, but recurrent, analytic and nonanalytic challenges in performing next- generation sequencing-based assays and point to remedies to help laboratories improve performance.

DNA damage in cancer development: special implications in viral oncogenesis

DNA lesions arise from a combination of physiological/metabolic sources and exogenous environmental influences. When left unrepaired, these alterations accumulate in the cells and can give rise to mutations that change the function of important proteins (i.e. tumor suppressors, oncoproteins), or cause chromosomal rearrangements (i.e. gene fusions) that also result in the deregulation of key cellular molecules. Progressive acquisition of such genetic changes promotes uncontrolled cell proliferation and evasion of cell death, and hence plays a key role in carcinogenesis. Another less-studied consequence of DNA damage accumulating in the host genome is the integration of oncogenic DNA viruses such as Human papillomavirus, Merkel cell polyomavirus, and Hepatitis B virus. This critical step of viral-induced carcinogenesis is thought to be particularly facilitated by DNA breaks in both viral and host genomes. Therefore, the impact of DNA damage on carcinogenesis is magnified in the case of such oncoviruses via the additional effect of increasing integration frequency. In this review, we briefly present the various endogenous and exogenous factors that cause different types of DNA damage. Next, we discuss the contribution of these lesions in cancer development. Finally, we examine the amplified effect of DNA damage in viral-induced oncogenesis and summarize the limited data existing in the literature related to DNA damage-induced viral integration. To conclude, additional research is needed to assess the DNA damage pathways involved in the transition from viral infection to cancer. Discovering that a certain DNA damaging agent increases the likelihood of viral integration will enable the development of prophylactic and therapeutic strategies designed specifically to prevent such integration, with an ultimate goal of reducing or eliminating these viral-induced malignancies.

Nonspecific Synthesis in the Reactions of Isothermal Nucleic Acid Amplification

The review focuses on the main factors involved in the formation of nonspecific products in isothermal nucleic acid amplification, such as mispriming, ab initio DNA synthesis, and additional activities of DNA polymerases, and discusses approaches to prevent formation of such nonspecific products in LAMP, RPA, NASBA, RCA, SDA, LSDA, NDA, and EXPAR.

Profiling variable-number tandem repeat variation across populations using repeat-pangenome graphs

Variable number tandem repeats (VNTRs) are composed of consecutive repetitive DNA with hypervariable repeat count and composition. They include protein coding sequences and associations with clinical disorders. It has been difficult to incorporate VNTR analysis in disease studies that use short-read sequencing because the traditional approach of mapping to the human reference is less effective for repetitive and divergent sequences. In this work, we solve VNTR mapping for short reads with a repeat-pangenome graph (RPGG), a data structure that encodes both the population diversity and repeat structure of VNTR loci from multiple haplotype-resolved assemblies. We develop software to build a RPGG, and use the RPGG to estimate VNTR composition with short reads. We use this to discover VNTRs with length stratified by continental population, and expression quantitative trait loci, indicating that RPGG analysis of VNTRs will be critical for future studies of diversity and disease.