DNA mutation motifs in the genes associated with inherited diseases

K-mer-based Approaches to Bridging Pangenomics and Population Genetics

Many commonly studied species now have more than one chromosome-scale genome assembly, revealing a large amount of genetic diversity previously missed by approaches that map short reads to a single reference. However, many species still lack multiple reference genomes and correctly aligning references to build pangenomes can be challenging for many species, limiting our ability to study this missing genomic variation in population genetics. Here, we argue that k-mers are a very useful but underutilized tool for bridging the reference-focused paradigms of population genetics with the reference-free paradigms of pangenomics. We review current literature on the uses of k-mers for performing three core components of most population genetics analyses: identifying, measuring, and explaining patterns of genetic variation. We also demonstrate how different k-mer-based measures of genetic variation behave in population genetic simulations according to the choice of k, depth of sequencing coverage, and degree of data compression. Overall, we find that k-mer-based measures of genetic diversity scale consistently with pairwise nucleotide diversity (π) up to values of about π=0.025 (R2=0.97) for neutrally evolving populations. For populations with even more variation, using shorter k-mers will maintain the scalability up to at least π=0.1. Furthermore, in our simulated populations, k-mer dissimilarity values can be reliably approximated from counting bloom filters, highlighting a potential avenue to decreasing the memory burden of k-mer-based genomic dissimilarity analyses. For future studies, there is a great opportunity to further develop methods to identifying selected loci using k-mers.

An Update on Neuroaging on Earth and in Spaceflight

Over 400 articles on the pathophysiology of brain aging, neuroaging, and neurodegeneration were reviewed, with a focus on epigenetic mechanisms and numerous non-coding RNAs. In particular, this review the accent is on microRNAs, the discovery of whose pivotal role in gene regulation was recognized by the 2024 Nobel Prize in Physiology or Medicine. Aging is not a gradual process that can be easily modeled and described. Instead, multiple temporal processes occur during aging, and they can lead to mosaic changes that are not uniform in pace. The rate of change depends on a combination of external and internal factors and can be boosted in accelerated aging. The rate can decrease in decelerated aging due to individual structural and functional reserves created by cognitive, physical training, or pharmacological interventions. Neuroaging can be caused by genetic changes, epigenetic modifications, oxidative stress, inflammation, lifestyle, and environmental factors, which are especially noticeable in space environments where adaptive changes can trigger aging-like processes. Numerous candidate molecular biomarkers specific to neuroaging need to be validated to develop diagnostics and countermeasures.

Reappraisal of the Concept of Accelerated Aging in Neurodegeneration and Beyond

BACKGROUND

Genetic and epigenetic changes, oxidative stress and inflammation influence the rate of aging, which diseases, lifestyle and environmental factors can further accelerate. In accelerated aging (AA), the biological age exceeds the chronological age.

OBJECTIVE

The objective of this study is to reappraise the AA concept critically, considering its weaknesses and limitations.

METHODS

We reviewed more than 300 recent articles dealing with the physiology of brain aging and neurodegeneration pathophysiology.

RESULTS

(1) Application of the AA concept to individual organs outside the brain is challenging as organs of different systems age at different rates. (2) There is a need to consider the deceleration of aging due to the potential use of the individual structure-functional reserves. The latter can be restored by pharmacological and/or cognitive therapy, environment, etc. (3) The AA concept lacks both standardised terminology and methodology. (4) Changes in specific molecular biomarkers (MBM) reflect aging-related processes; however, numerous MBM candidates should be validated to consolidate the AA theory. (5) The exact nature of many potential causal factors, biological outcomes and interactions between the former and the latter remain largely unclear.

CONCLUSIONS

Although AA is commonly recognised as a perspective theory, it still suffers from a number of gaps and limitations that assume the necessity for an updated AA concept.

Whole exome sequencing in Alopecia Areata identifies rare variants in KRT82

Alopecia areata is a complex genetic disease that results in hair loss due to the autoimmune-mediated attack of the hair follicle. We previously defined a role for both rare and common variants in our earlier GWAS and linkage studies. Here, we identify rare variants contributing to Alopecia Areata using a whole exome sequencing and gene-level burden analyses approach on 849 Alopecia Areata patients compared to 15,640 controls. KRT82 is identified as an Alopecia Areata risk gene with rare damaging variants in 51 heterozygous Alopecia Areata individuals (6.01%), achieving genome-wide significance (p = 2.18E−07). KRT82 encodes a hair-specific type II keratin that is exclusively expressed in the hair shaft cuticle during anagen phase, and its expression is decreased in Alopecia Areata patient skin and hair follicles. Finally, we find that cases with an identified damaging KRT82 variant and reduced KRT82 expression have elevated perifollicular CD8 infiltrates. In this work, we utilize whole exome sequencing to successfully identify a significant Alopecia Areata disease-relevant gene, KRT82, and reveal a proposed mechanism for rare variant predisposition leading to disrupted hair shaft integrity.

Common variants have been discovered to be associated with Alopecia Areata; however, rare variants have been less well studied. Here, the authors use whole-exome sequencing to identify associated rare variants in the hair keratin gene KRT82. Further, they find that individuals with Alopecia Areata have reduced expression of KRT82 in the skin and hair follicle.

Design of Oligonucleotides for Allele-Specific Amplification Based on PCR and Isothermal Techniques

Single-nucleotide variations have been associated to various genetic diseases, variations on drug efficiency, and differences in cancer prognostics. The detection of these changes in nucleic acid sequences from patient samples is particularly useful for accurate diagnosis, therapeutics, and disease management. A reliable allele-specific amplification is still an important challenge for molecular-based diagnostic technologies. In the last years, allele-specific primers have been designed for promoting the enrichment of certain variants, based on a higher stability of primer/template duplexes. Also, several methods are based on the addition of a blocking oligonucleotide that prevent the amplification of a specific variant, enabling that other DNA variants can be observed. In this context, genotyping methods based on isothermal amplification techniques are increasing, especially those assays aimed to be deployed at point-of-care applications. The correct selection of target sequences is crucial for reaching the required analytical performances, in terms of reaction time, amplification yield, and selectivity. The present chapter describes the design criteria for the selection of primers and blockers for relevant PCR approaches and novel isothermal strategies. Several successful examples are provided in order to highlight the main design restrictions and the potential to be extended to other applications.

Bending of Canonical and G/T Mismatched DNAs

Mismatched base pairs alter the flexibility and intrinsic curvature of DNA. The role of such DNA features is not fully understood in the mismatch repair pathway. MutS/DNA complexes exhibit DNA bending, PHE intercalation, and changes of base-pair parameters near the mismatch. Recently, we have shown that base-pair opening in the absence of MutS can discriminate mismatches from canonical base pairs better than DNA bending. However, DNA bending in the absence of MutS was found to be rather challenging to describe correctly. Here, we present a computational study on the DNA bending of canonical and G/T mismatched DNAs. Five types of geometric parameters covering template-based bending toward the experimental DNA structure, global, and local geometry parameters were employed in biased molecular dynamics in the absence of MutS. None of these parameters showed higher discrimination than the base-pair opening. Only roll could induce a sharply localized bending of DNA as observed in the experimental MutS/DNA structure. Further, we demonstrated that the intercalation of benzene mimicking PHE decreases the energetic cost of DNA bending without any effect on mismatch discrimination.

A recurrent SHANK3 frameshift variant in Autism Spectrum Disorder

Autism Spectrum Disorder (ASD) is genetically complex with ~100 copy number variants and genes involved. To try to establish more definitive genotype and phenotype correlations in ASD, we searched genome sequence data, and the literature, for recurrent predicted damaging sequence-level variants affecting single genes. We identified 18 individuals from 16 unrelated families carrying a heterozygous guanine duplication (c.3679dup; p.Ala1227Glyfs*69) occurring within a string of 8 guanines (genomic location [hg38]g.50,721,512dup) affecting SHANK3, a prototypical ASD gene (0.08% of ASD-affected individuals carried the predicted p.Ala1227Glyfs*69 frameshift variant). Most probands carried de novo mutations, but five individuals in three families inherited it through somatic mosaicism. We scrutinized the phenotype of p.Ala1227Glyfs*69 carriers, and while everyone (17/17) formally tested for ASD carried a diagnosis, there was the variable expression of core ASD features both within and between families. Defining such recurrent mutational mechanisms underlying an ASD outcome is important for genetic counseling and early intervention.

Importance of base-pair opening for mismatch recognition

Mismatch repair is a highly conserved cellular pathway responsible for repairing mismatched dsDNA. Errors are detected by the MutS enzyme, which most likely senses altered mechanical property of damaged dsDNA rather than a specific molecular pattern. While the curved shape of dsDNA in crystallographic MutS/DNA structures suggests the role of DNA bending, the theoretical support is not fully convincing. Here, we present a computational study focused on a base-pair opening into the minor groove, a specific base-pair motion observed upon interaction with MutS. Propensities for the opening were evaluated in terms of two base-pair parameters: Opening and Shear. We tested all possible base pairs in anti/anti, anti/syn and syn/anti orientations and found clear discrimination between mismatches and canonical base-pairs only for the opening into the minor groove. Besides, the discrimination gap was also confirmed in hotspot and coldspot sequences, indicating that the opening could play a more significant role in the mismatch recognition than previously recognized. Our findings can be helpful for a better understanding of sequence-dependent mutability. Further, detailed structural characterization of mismatches can serve for designing anti-cancer drugs targeting mismatched base pairs.

Density Functional Theory in the Prediction of Mutagenicity: A Perspective

As a field, computational toxicology is concerned with using in silico models to predict and understand the origins of toxicity. It is fast, relatively inexpensive, and avoids the ethical conundrum of using animals in scientific experimentation. In this perspective, we discuss the importance of computational models in toxicology, with a specific focus on the different model types that can be used in predictive toxicological approaches toward mutagenicity (SARs and QSARs). We then focus on how quantum chemical methods, such as density functional theory (DFT), have previously been used in the prediction of mutagenicity. It is then discussed how DFT allows for the development of new chemical descriptors that focus on capturing the steric and energetic effects that influence toxicological reactions. We hope to demonstrate the role that DFT plays in understanding the fundamental, intrinsic chemistry of toxicological reactions in predictive toxicology.

Bending of DNA duplexes with mutation motifs

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

Review of Different Sequence Motif Finding Algorithms

The DNA motif discovery is a primary step in many systems for studying gene function. Motif discovery plays a vital role in identification of Transcription Factor Binding Sites (TFBSs) that help in learning the mechanisms for regulation of gene expression. Over the past decades, different algorithms were used to design fast and accurate motif discovery tools. These algorithms are generally classified into consensus or probabilistic approaches that many of them are time-consuming and easily trapped in a local optimum. Nature-inspired algorithms and many of combinatorial algorithms are recently proposed to overcome these problems. This paper presents a general classification of motif discovery algorithms with new sub-categories that facilitate building a successful motif discovery algorithm. It also presents a summary of comparison between them.

Simulating DNA Chip Design Using All-Electronic Graphene-Based Substrates

In this paper, we present a theoretical investigation of an all-electronic biochip based on graphene to detect DNA including a full dynamical treatment for the environment. Our proposed device design is based on the changes in the electronic transport properties of graphene interacting with DNA strands under the effect of the solvent. To investigate these systems, we applied a hybrid methodology, combining quantum and classical mechanics (QM/MM) coupled to non-equilibrium Green's functions, allowing for the calculations of electronic transport. Our results show that the proposed device has high sensitivity towards the presence of DNA, and, combined with the presence of a specific DNA probe in the form of a single-strand, it presents good selectivity towards specific nucleotide sequences.

iPS-Cell Technology and the Problem of Genetic Instability-Can It Ever Be Safe for Clinical Use?

The use of induced Pluripotent Stem Cells (iPSC) as a source of autologous tissues shows great promise in regenerative medicine. Nevertheless, several major challenges remain to be addressed before iPSC-derived cells can be used in therapy, and experience of their clinical use is extremely limited. In this review, the factors affecting the safe translation of iPSC to the clinic are considered, together with an account of efforts being made to overcome these issues. The review draws upon experiences with pluripotent stem-cell therapeutics, including clinical trials involving human embryonic stem cells and the widely transplanted mesenchymal stem cells. The discussion covers concerns relating to: (i) the reprogramming process; (ii) the detection and removal of incompletely differentiated and pluripotent cells from the resulting medicinal products; and (iii) genomic and epigenetic changes, and the evolutionary and selective processes occurring during culture expansion, associated with production of iPSC-therapeutics. In addition, (iv) methods for the practical culture-at-scale and standardization required for routine clinical use are considered. Finally, (v) the potential of iPSC in the treatment of human disease is evaluated in the light of what is known about the reprogramming process, the behavior of cells in culture, and the performance of iPSC in pre-clinical studies.

Role of Inosine⁻Uracil Base Pairs in the Canonical RNA Duplexes

Adenosine to inosine (A–I) editing is the most common modification of double-stranded RNA (dsRNA). This change is mediated by adenosine deaminases acting on RNA (ADARs) enzymes with a preference of U>A>C>G for 5′ neighbor and G>C=A>U or G>C>U=A for 3′ neighbor. A–I editing occurs most frequently in the non-coding regions containing repetitive elements such as ALUs. It leads to disruption of RNA duplex structure, which prevents induction of innate immune response. We employed standard and biased molecular dynamics (MD) simulations to analyze the behavior of RNA duplexes with single and tandem inosine–uracil (I–U) base pairs in different sequence context. Our analysis showed that the I–U pairs induce changes in base pair and base pair step parameters and have different dynamics when compared with standard canonical base pairs. In particular, the first I–U pair from tandem I–U/I–U systems exhibited increased dynamics depending on its neighboring 5′ base. We discovered that UII sequence, which is frequently edited, has lower flexibility compared with other sequences (AII, GII, CII), hence it only modestly disrupts dsRNA. This might indicate that the UAA motifs in ALUs do not have to be sufficiently effective in preventing immune signaling.

Truncating Variants in NAA15 Are Associated with Variable Levels of Intellectual Disability, Autism Spectrum Disorder, and Congenital Anomalies

N-alpha-acetylation is a common co-translational protein modification that is essential for normal cell function in humans. We previously identified the genetic basis of an X-linked infantile lethal Mendelian disorder involving a c.109T>C (p.Ser37Pro) missense variant in NAA10, which encodes the catalytic subunit of the N-terminal acetyltransferase A (NatA) complex. The auxiliary subunit of the NatA complex, NAA15, is the dimeric binding partner for NAA10. Through a genotype-first approach with whole-exome or genome sequencing (WES/WGS) and targeted sequencing analysis, we identified and phenotypically characterized 38 individuals from 33 unrelated families with 25 different de novo or inherited, dominantly acting likely gene disrupting (LGD) variants in NAA15. Clinical features of affected individuals with LGD variants in NAA15 include variable levels of intellectual disability, delayed speech and motor milestones, and autism spectrum disorder. Additionally, mild craniofacial dysmorphology, congenital cardiac anomalies, and seizures are present in some subjects. RNA analysis in cell lines from two individuals showed degradation of the transcripts with LGD variants, probably as a result of nonsense-mediated decay. Functional assays in yeast confirmed a deleterious effect for two of the LGD variants in NAA15. Further supporting a mechanism of haploinsufficiency, individuals with copy-number variant (CNV) deletions involving NAA15 and surrounding genes can present with mild intellectual disability, mild dysmorphic features, motor delays, and decreased growth. We propose that defects in NatA-mediated N-terminal acetylation (NTA) lead to variable levels of neurodevelopmental disorders in humans, supporting the importance of the NatA complex in normal human development.