Hotspots of Human Mutation

The long and short of hyperdivergent regions

The increasing prevalence of genome sequencing and assembly has uncovered evidence of hyperdivergent genomic regions - loci with excess genetic diversity - in species across the tree of life. Hyperdivergent regions are often enriched for genes that mediate environmental responses, such as immunity, parasitism, and sensory perception. Especially in self-fertilizing species where the majority of the genome is homozygous, the existence of hyperdivergent regions might imply the historical action of evolutionary forces such as introgression and/or balancing selection. We anticipate that the application of new sequencing technologies, broader taxonomic sampling, and evolutionary modeling of hyperdivergent regions will provide insights into the mechanisms that generate and maintain genetic diversity within and between species.

The theory of massively repeated evolution and full identifications of cancer-driving nucleotides (CDNs)

Tumorigenesis, like most complex genetic traits, is driven by the joint actions of many mutations. At the nucleotide level, such mutations are cancer-driving nucleotides (CDNs). The full sets of CDNs are necessary, and perhaps even sufficient, for the understanding and treatment of each cancer patient. Currently, only a small fraction of CDNs is known as most mutations accrued in tumors are not drivers. We now develop the theory of CDNs on the basis that cancer evolution is massively repeated in millions of individuals. Hence, any advantageous mutation should recur frequently and, conversely, any mutation that does not is either a passenger or deleterious mutation. In the TCGA cancer database (sample size n=300-1000), point mutations may recur in i out of n patients. This study explores a wide range of mutation characteristics to determine the limit of recurrences (i*) driven solely by neutral evolution. Since no neutral mutation can reach i*=3, all mutations recurring at i≥3 are CDNs. The theory shows the feasibility of identifying almost all CDNs if n increases to 100,000 for each cancer type. At present, only <10% of CDNs have been identified. When the full sets of CDNs are identified, the evolutionary mechanism of tumorigenesis in each case can be known and, importantly, gene targeted therapy will be far more effective in treatment and robust against drug resistance.

Gastric Organoid, a Promising Modeling for Gastric Stem Cell Homeostasis and Therapeutic Application

The elucidation of the pathophysiology underlying various diseases necessitates the development of research platforms that faithfully mimic in vivo conditions. Traditional model systems such as two-dimensional cell cultures and animal models have proven inadequate in capturing the complexities of human disease modeling. However, recent strides in organoid culture systems have opened up new avenues for comprehending gastric stem cell homeostasis and associated diseases, notably gastric cancer. Given the significance of gastric cancer, a thorough understanding of its pathophysiology and molecular underpinnings is imperative. To this end, the utilization of patient-derived organoid libraries emerges as a remarkable platform, as it faithfully mirrors patient-specific characteristics, including mutation profiles and drug sensitivities. Furthermore, genetic manipulation of gastric organoids facilitates the exploration of molecular mechanisms underlying gastric cancer development. This review provides a comprehensive overview of recent advancements in various adult stem cell-derived gastric organoid models and their diverse applications.

Detection and analysis of complex structural variation in human genomes across populations and in brains of donors with psychiatric disorders

Complex structural variations (cxSVs) are often overlooked in genome analyses due to detection challenges. We developed ARC-SV, a probabilistic and machine-learning-based method that enables accurate detection and reconstruction of cxSVs from standard datasets. By applying ARC-SV across 4,262 genomes representing all continental populations, we identified cxSVs as a significant source of natural human genetic variation. Rare cxSVs have a propensity to occur in neural genes and loci that underwent rapid human-specific evolution, including those regulating corticogenesis. By performing single-nucleus multiomics in postmortem brains, we discovered cxSVs associated with differential gene expression and chromatin accessibility across various brain regions and cell types. Additionally, cxSVs detected in brains of psychiatric cases are enriched for linkage with psychiatric GWAS risk alleles detected in the same brains. Furthermore, our analysis revealed significantly decreased brain-region- and cell-type-specific expression of cxSV genes, specifically for psychiatric cases, implicating cxSVs in the molecular etiology of major neuropsychiatric disorders.

A Hypermutable Region in the DISP2 Gene Links to Natural Selection and Late-Onset Neurocognitive Disorders in Humans

(CCG) short tandem repeats (STRs) are predominantly enriched in genic regions, mutation hotspots for C to T truncating substitutions, and involved in various neurological and neurodevelopmental disorders. However, intact blocks of this class of STRs are widely overlooked with respect to their link with natural selection. The human neuron-specific gene, DISP2 (dispatched RND transporter family member 2), contains a (CCG) repeat in its 5' untranslated region. Here, we sequenced this STR in a sample of 448 Iranian individuals, consisting of late-onset neurocognitive disorder (NCD) (N = 203) and controls (N = 245). We found that the region spanning the (CCG) repeat was highly mutated, resulting in several flanking (CCG) residues. However, an 8-repeat of the (CCG) repeat was predominantly abundant (frequency = 0.92) across the two groups. While the overall distribution of genotypes was not different between the two groups (p > 0.05), we detected four genotypes in the NCD group only (2% of the NCD genotypes, Mid-p = 0.02), consisting of extreme short alleles, 5- and 6-repeats, that were not detected in the control group. The patients harboring those genotypes received the diagnoses of probable Alzheimer's disease and vascular dementia. We also found six genotypes in the control group only (2.5% of the control genotypes, Mid-p = 0.01) that consisted of the 8-repeat and extreme long alleles, 9- and 10-repeats, of which the 10-repeat was not detected in the NCD group. The (CCG) repeat specifically expanded in primates. In conclusion, we report an indication of natural selection at a novel hypermutable region in the human genome and divergent alleles and genotypes in late-onset NhCDs and controls. These findings reinforce the hypothesis that a collection of rare alleles and genotypes in a number of genes may unambiguously contribute to the cognition impairment component of late-onset NCDs.

The human SMAD9 (GCC) repeat links to natural selection and late-onset neurocognitive disorders

INTRODUCTION

Whereas (GCC)-repeats are overrepresented in genic regions, and mutation hotspots, they are largely unexplored with regard to their link with natural selection. Across numerous primate species and tissues, SMAD9 (SMAD Family Member 9) reaches highest level of expression in the human brain. This gene contains a (GCC)-repeat in the interval between + 1 and + 60 of the transcription start site, which is in the high-ranking (GCC)-repeats with respect to length.

METHODS

Here we sequenced this (GCC)-repeat in 396 Iranian individuals, consisting of late-onset neurocognitive disorder (NCD) (N = 181) and controls (N = 215).

RESULTS

We detected two predominantly abundant alleles of 7 and 9 repeats, forming 96.2% of the allele pool. The (GCC)7/(GCC)9 ratio was in the reverse order in the NCD group versus controls (p = 0.005), resulting from excess of (GCC)7 in the NCD group (p = 0.003) and (GCC)9 in the controls (p = 0.01). Five genotypes, predominantly consisting of (GCC)7 and lacking (GCC)9 were detected in the NCD group only (p = 0.008). The patients harboring those genotypes received the diagnoses of Alzheimer's disease (AD) and vascular dementia (VD). Five genotypes consisting of (GCC)9 and lacking (GCC)7 were detected in the control group only (p = 0.002). The group-specific genotypes formed approximately 4% of the genotype pool in the human samples studied.

CONCLUSION

We propose natural selection and a novel locus for late-onset AD and VD at the SMAD9 (GCC)-repeat in humans.

Towards Reliable Detection of Introgression in the Presence of Among-Species Rate Variation

The role of interspecific hybridization has recently seen increasing attention, especially in the context of diversification dynamics. Genomic research has now made it abundantly clear that both hybridization and introgression-the exchange of genetic material through hybridization and backcrossing-are far more common than previously thought. Besides cases of ongoing or recent genetic exchange between taxa, an increasing number of studies report "ancient introgression"- referring to results of hybridization that took place in the distant past. However, it is not clear whether commonly used methods for the detection of introgression are applicable to such old systems, given that most of these methods were originally developed for analyses at the level of populations and recently diverged species, affected by recent or ongoing genetic exchange. In particular, the assumption of constant evolutionary rates, which is implicit in many commonly used approaches, is more likely to be violated as evolutionary divergence increases. To test the limitations of introgression detection methods when being applied to old systems, we simulated thousands of genomic datasets under a wide range of settings, with varying degrees of among-species rate variation and introgression. Using these simulated datasets, we showed that some commonly applied statistical methods, including the D-statistic and certain tests based on sets of local phylogenetic trees, can produce false-positive signals of introgression between divergent taxa that have different rates of evolution. These misleading signals are caused by the presence of homoplasies occurring at different rates in different lineages. To distinguish between the patterns caused by rate variation and genuine introgression, we developed a new test that is based on the expected clustering of introgressed sites along the genome and implemented this test in the program Dsuite.

Genome organization and stability in mammalian pre-implantation development

A largely stable genome is required for normal development, even as genetic change is an integral aspect of reproduction, genetic adaptation and evolution. Recent studies highlight a critical window of mammalian development with intrinsic DNA replication stress and genome instability in the first cell divisions after fertilization. Patterns of DNA replication and genome stability are established very early in mammals, alongside patterns of nuclear organization, and before the emergence of gene expression patterns, and prior to cell specification and germline formation. The study of DNA replication and genome stability in the mammalian embryo provides a unique cellular system due to the resetting of the epigenome to a totipotent state, and the de novo establishment of the patterns of nuclear organization, gene expression, DNA methylation, histone modifications and DNA replication. Studies on DNA replication and genome stability in the early mammalian embryo is relevant for understanding both normal and disease-causing genetic variation, and to uncover basic principles of genome regulation.

A transposase-derived gene required for human brain development

DNA transposable elements and transposase-derived genes are present in most living organisms, including vertebrates, but their function is largely unknown. PiggyBac Transposable Element Derived 5 (PGBD5) is an evolutionarily conserved vertebrate DNA transposase-derived gene with retained nuclease activity in human cells. Vertebrate brain development is known to be associated with prominent neuronal cell death and DNA breaks, but their causes and functions are not well understood. Here, we show that PGBD5 contributes to normal brain development in mice and humans, where its deficiency causes disorder of intellectual disability, movement, and seizures. In mice, Pgbd5 is required for the developmental induction of post-mitotic DNA breaks and recurrent somatic genome rearrangements. In the brain cortex, loss of Pgbd5 leads to aberrant differentiation and gene expression of distinct neuronal populations, including specific types of glutamatergic neurons, which explains the features of PGBD5 deficiency in humans. Thus, PGBD5 might be a transposase-derived enzyme required for brain development in mammals.

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

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

Machine learning enables pan-cancer identification of mutational hotspots at persistent CTCF binding sites

CCCTC-binding factor (CTCF) is an insulator protein that binds to a highly conserved DNA motif and facilitates regulation of three-dimensional (3D) nuclear architecture and transcription. CTCF binding sites (CTCF-BSs) reside in non-coding DNA and are frequently mutated in cancer. Our previous study identified a small subclass of CTCF-BSs that are resistant to CTCF knock down, termed persistent CTCF binding sites (P-CTCF-BSs). P-CTCF-BSs show high binding conservation and potentially regulate cell-type constitutive 3D chromatin architecture. Here, using ICGC sequencing data we made the striking observation that P-CTCF-BSs display a highly elevated mutation rate in breast and prostate cancer when compared to all CTCF-BSs. To address whether P-CTCF-BS mutations are also enriched in other cell-types, we developed CTCF-INSITE-a tool utilising machine learning to predict persistence based on genetic and epigenetic features of experimentally-determined P-CTCF-BSs. Notably, predicted P-CTCF-BSs also show a significantly elevated mutational burden in all 12 cancer-types tested. Enrichment was even stronger for P-CTCF-BS mutations with predicted functional impact to CTCF binding and chromatin looping. Using in vitro binding assays we validated that P-CTCF-BS cancer mutations, predicted to be disruptive, indeed reduced CTCF binding. Together this study reveals a new subclass of cancer specific CTCF-BS DNA mutations and provides insights into their importance in genome organization in a pan-cancer setting.

Effects of PCNA Stability on the Formation of Mutations

The fidelity of replication, especially in the presence of DNA damage, is essential for the proper function of cells. Mutations that inactivate genes involved in DNA damage repair or bypass are enriched in several types of cancer cells. Thus, it is important to further our understanding of the mechanisms governing replication fidelity. PCNA is a ring-shaped complex that encircles DNA at the front of the replication fork, at the double-stranded/single-stranded DNA junction. It serves as a processivity factor for the different DNA replication polymerases, allowing them to replicate longer stretches of DNA by physically tethering them to the DNA and preventing their detachment. In addition, PCNA also regulates and coordinates different DNA damage bypass pathways meant to allow DNA replication in the presence of DNA damage. Due to its essentiality and the numerous functions it has in the cell, much is still unclear about PCNA. Here, we utilize PCNA mutants that lower the stability of the PCNA complex on the chromatin, and thus tend to disassociate and fall from the DNA. Using these mutants, we show that PCNA's physical presence on the DNA can prevent DNA misalignment at repetitive sequences, leading to increased mutation formation. We also show that PCNA-interacting proteins play an important role in strengthening the ring's stability on the chromatin. Such repetitive sequence-induced mutations are common in several human diseases and it is important to study their formation and the mechanisms guarding against them.

Clinical characteristics of BRAT1-related disease: a systematic literature review

BACKGROUND

BRAT1 (BRCA1-associated ataxia telangiectasia mutated activator 1) is involved in many important biological processes, including DNA damage response and maintenance of mitochondrial homeostasis. Dysfunctional BRAT1 causes variable clinical phenotypes, which hinders BRAT1-related disease from recognition and diagnosis.

METHODS

Preferred Reporting Items for Systematic Reviews and Meta-Analyses statement was the guideline for this systematic review. MEDLINE was searched by terms ("BAAT1" and "BRAT1") from inception until June 21, 2022.

RESULTS

Twenty-eight studies, screened out of 49 records, were included for data extraction. The data from fifty patients with mutated BRAT1 were collected. There are 3 high relevant phenotypes, 4 medium relevant phenotypes and 3 low relevant phenotypes. Eye-related abnormal features were most frequently reported: 27 abnormal features were observed. Thirty-nine kinds of pathogenic nucleotide change in BRAT1 were reported. Top three common mutations of BRAT1 were c.638_639insA (16 cases), c.1395G > A (5 cases) and c.294dupA (4 cases). Homozygous mutations in BRAT1 presented a more severe phenotype than those who are compound heterozygotes.

CONCLUSIONS

This is the first comprehensive systematic review to present quantitative data about clinical characteristics of BRAT1-related disease, which helps doctors to recognize and diagnose it easier.

The Engineered Drug 3'UTRMYC1-18 Degrades the c-MYC-STAT5A/B-PD-L1 Complex In Vivo to Inhibit Metastatic Triple-Negative Breast Cancer

c-MYC is overexpressed in 70% of human cancers, including triple-negative breast cancer (TNBC), yet there is no clinically approved drug that directly targets it. Here, we engineered the mRNA-stabilizing poly U sequences within the 3'UTR of c-MYC to specifically destabilize and promote the degradation of c-MYC transcripts. Interestingly, the engineered derivative outcompetes the endogenous overexpressed c-MYC mRNA, leading to reduced c-MYC mRNA and protein levels. The iron oxide nanocages (IO-nanocages) complexed with MYC-destabilizing constructs inhibited primary and metastatic tumors in mice bearing TNBC and significantly prolonged survival by degrading the c-MYC-STAT5A/B-PD-L1 complexes that drive c-MYC-positive TNBC. Taken together, we have described a novel therapy for c-MYC-driven TNBC and uncovered c-MYC-STAT5A/B-PD-L1 interaction as the target.

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

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

A molecular and cellular perspective on human brain evolution and tempo

The evolution of the modern human brain was accompanied by distinct molecular and cellular specializations, which underpin our diverse cognitive abilities but also increase our susceptibility to neurological diseases. These features, some specific to humans and others shared with related species, manifest during different stages of brain development. In this multi-stage process, neural stem cells proliferate to produce a large and diverse progenitor pool, giving rise to excitatory or inhibitory neurons that integrate into circuits during further maturation. This process unfolds over varying time scales across species and has progressively become slower in the human lineage, with differences in tempo correlating with differences in brain size, cell number and diversity, and connectivity. Here we introduce the terms 'bradychrony' and 'tachycrony' to describe slowed and accelerated developmental tempos, respectively. We review how recent technical advances across disciplines, including advanced engineering of in vitro models, functional comparative genetics and high-throughput single-cell profiling, are leading to a deeper understanding of how specializations of the human brain arise during bradychronic neurodevelopment. Emerging insights point to a central role for genetics, gene-regulatory networks, cellular innovations and developmental tempo, which together contribute to the establishment of human specializations during various stages of neurodevelopment and at different points in evolution.

Genomic Mosaicism of the Brain: Origin, Impact, and Utility

Genomic mosaicism describes the phenomenon where some but not all cells within a tissue harbor unique genetic mutations. Traditionally, research focused on the impact of genomic mosaicism on clinical phenotype-motivated by its involvement in cancers and overgrowth syndromes. More recently, we increasingly shifted towards the plethora of neutral mosaic variants that can act as recorders of cellular lineage and environmental exposures. Here, we summarize the current state of the field of genomic mosaicism research with a special emphasis on our current understanding of this phenomenon in brain development and homeostasis. Although the field of genomic mosaicism has a rich history, technological advances in the last decade have changed our approaches and greatly improved our knowledge. We will provide current definitions and an overview of contemporary detection approaches for genomic mosaicism. Finally, we will discuss the impact and utility of genomic mosaicism.

Founder mutations and rare disease in the Arab world

Founder mutations are disease-causing variants that occur frequently in geographically or culturally isolated groups whose shared ancestor(s) carried the pathogenic variant. While some disease alleles may vanish from the genetic pool due to natural selection, variants with weaker effects may survive for a long time, thereby enhancing the prevalence of some rare diseases. These are predominantly autosomal recessive diseases but can also be autosomal dominant traits with late-onset or mild phenotypes. Cultural practices, such as endogamy and consanguinity, in these isolated groups lead to higher prevalence of such rare diseases compared to the rest of the population and worldwide. In this Perspective, we define population isolates and the underlying genetic mechanisms for accumulating founder mutations. We also discuss the current and potential scientific, clinical and public-health implications of studying founder mutations in population isolates around the world, with a particular focus on the Arab population.

HATCHet2: clone- and haplotype-specific copy number inference from bulk tumor sequencing data

Bulk DNA sequencing of multiple samples from the same tumor is becoming common, yet most methods to infer copy-number aberrations (CNAs) from this data analyze individual samples independently. We introduce HATCHet2, an algorithm to identify haplotype- and clone-specific CNAs simultaneously from multiple bulk samples. HATCHet2 extends the earlier HATCHet method by improving identification of focal CNAs and introducing a novel statistic, the minor haplotype B-allele frequency (mhBAF), that enables identification of mirrored-subclonal CNAs. We demonstrate HATCHet2's improved accuracy using simulations and a single-cell sequencing dataset. HATCHet2 analysis of 10 prostate cancer patients reveals previously unreported mirrored-subclonal CNAs affecting cancer genes.

Assisting the analysis of insertions and deletions using regional allele frequencies

Accurate estimation of population allele frequency (AF) is crucial for gene discovery and genetic diagnostics. However, determining AF for frameshift-inducing small insertions and deletions (indels) faces challenges due to discrepancies in mapping and variant calling methods. Here, we propose an innovative approach to assess indel AF. We developed CRAFTS-indels (Calculating Regional Allele Frequency Targeting Small indels), an algorithm that combines AF of distinct indels within a given region and provides "regional AF" (rAF). We tested and validated CRAFTS-indels using three independent datasets: gnomAD v2 (n=125,748 samples), an internal dataset (IGM; n=39,367), and the UK BioBank (UKBB; n=469,835). By comparing rAF against standard AF, we identified rare indels with rAF exceeding standard AF (sAF≤10-4 and rAF>10-4) as "rAF-hi" indels. Notably, a high percentage of rare indels were "rAF-hi", with a higher proportion in gnomAD v2 (11-20%) and IGM (11-22%) compared to the UKBB (5-9% depending on the CRAFTS-indels' parameters). Analysis of the overlap of regions based on their rAF with low complexity regions and with ClinVar classification supported the pertinence of rAF. Using the internal dataset, we illustrated the utility of CRAFTS-indel in the analysis of de novo variants and the potential negative impact of rAF-hi indels in gene discovery. In summary, annotation of indels with cohort specific rAF can be used to handle some of the limitations of current annotation pipelines and facilitate detection of novel gene disease associations. CRAFTS-indels offers a user-friendly approach to providing rAF annotation. It can be integrated into public databases such as gnomAD, UKBB and used by ClinVar to revise indel classifications.

CardioHotspots: a database of mutational hotspots for cardiac disorders

Mutational hotspots are DNA regions with an abnormally high frequency of genetic variants. Identifying whether a variant is located in a mutational hotspot is critical for determining the variant's role in disorder predisposition, development, and treatment response. Despite their significance, current databases on mutational hotspots are limited to the oncology domain. However, identifying mutational hotspots is critical for any disorder in which genetics plays a role. This is true for the world's leading cause of death: cardiac disorders. In this work, we present CardioHotspots, a literature-based database of manually curated hotspots for cardiac diseases. This is the only database we know of that provides high-quality and easily accessible information about hotspots associated with cardiac disorders. CardioHotspots is publicly accessible via a web-based platform (https://genomics-hub.pros.dsic.upv.es:3099/). Database URL: https://genomics-hub.pros.dsic.upv.es:3099/.

The Segregation of p.Arg68Ter-CLDN14 Mutation in a Syrian Deaf Family, Phenotypic Variations, and Comparative Analysis with the GJB2 Gene

Hearing impairment, a rare inherited condition, is notably prevalent in populations with high rates of consanguinity. The most common form observed globally is autosomal recessive non-syndromic hearing loss. Despite its prevalence, this genetic disorder is characterized by a substantial genetic diversity, making diagnosis and screening challenging. The emergence of advanced next-generation sequencing (NGS) technologies has significantly advanced the discovery of genes and variants linked to various conditions, such as hearing loss. In this study, our objective was to identify the specific variant causing hearing loss in a family from Syria using clinical exome sequencing. The proband in the family exhibited profound deafness as shown by pure-tone audiometry results. The analysis of the different variants obtained by NGS revealed the presence of a nonsense mutation within the CLDN14 gene. Through Sanger sequencing, we verified that this variant segregates with the disease and was not present in the control population. Moreover, we conducted a comprehensive review of all reported deafness-related CLDN14 mutations and their associated phenotypes. Furthermore, we endeavored to carry out a comparative analysis between the CLDN14 and GJB2 genes, with the objective of identifying potential factors that could explain the notable discrepancy in mutation frequency between these two genes.

Extended Enrichment for Ultrasensitive Detection of Low-Frequency Mutations by Long Blocker Displacement Amplification

Detecting low-frequency DNA mutations hotspots cluster is critical for cancer diagnosis but remains challenging. Quantitative PCR (qPCR) is constrained by sensitivity, and allele-specific PCR is restricted by throughput. Here we develop a long blocker displacement amplification (LBDA) coupled with qPCR for ultrasensitive and multiplexed variants detection. By designing long blocker oligos to perfectly match wildtype sequences while mispairing with mutants, long blockers enable 14-44 nt enrichment regions which is 2-fold longer than normal BDA in the experiments. For wild template with a specific nucleotide, LBDA can detect different mutation types down to 0.5 % variant allele frequency (VAF) in one reaction, with median enrichment fold of 1,000 on 21 mutant DNA templates compared to the wild type. We applied LBDA-qPCR to detect KRAS and NRAS mutation hotspots, utilizing a single plex assay capable of covering 81 mutations and tested in synthetic templates and colorectal cancer tissue samples. Moreover, the mutation types were verified through Sanger sequencing, demonstrating concordance with results obtained from next generation sequencing. Overall, LBDA-qPCR provides a simple yet ultrasensitive approach for multiplexed detection of low VAF mutations hotspots, presenting a powerful tool for cancer diagnosis and monitoring.

Deletion in a regulatory region is associated with underexpression of miR-148b‑3p in patients with prostate cancer

Prostate cancer (PCa) is the leading cause of cancer-related death in men. This pathology is complex and heterogeneous; therefore, elucidating the molecular mechanisms that lead to its origin and progression is imperative. MicroRNAs (miRNAs or miRs) are part of the epigenetic machinery that regulates the expression of human genes, therefore, mutations in the genes that encode them can lead to a dysregulation in their expression, which directly impacts their target genes, which could be oncogenes or tumor suppressor genes. In PCa several dysregulated expression levels of miRNAs are associated with perturbed cellular processes. A differential expression of miRNAs such as miR-145-5p and miR-148-3p has been observed in PCa, possibly due to mutations in regions near the miRNAs. However, the molecular mechanisms that lead to the dysregulation of these miRNAs still need to be clarified. Therefore, the present study aimed to analyze the expression of miRNAs and their relationship with mutations in patients with and without PCa. In total, 71 patients were analyzed: 41 of whom had PCa (CAP group) and 30 with benign pathology (BPD group). Underexpression was observed in miR-145-5p and miR-148b-3p in PCa patients (P=0.03 and P=0.001, respectively). In miR-145-5p, no mutations related to its expression were identified. For miR-148b-3p, a set of mutations were identified in the chr12:54337042/54337043 region, which were grouped into the mutation named DelsAAG. Although this mutation's abnormal allele is related to PCa (P=0.017), a statistically significant difference was observed in the expression of miR-148b-3p between carriers and non-carriers of the mutated allele, identifying a mechanism likely to be involved in the miR-148b-3p dysregulation.

Positive Selection Drives cis-regulatory Evolution Across the Threespine Stickleback Y Chromosome

Allele-specific gene expression evolves rapidly on heteromorphic sex chromosomes. Over time, the accumulation of mutations on the Y chromosome leads to widespread loss of gametolog expression, relative to the X chromosome. It remains unclear if expression evolution on degrading Y chromosomes is primarily driven by mutations that accumulate through processes of selective interference, or if positive selection can also favor the down-regulation of coding regions on the Y chromosome that contain deleterious mutations. Identifying the relative rates of cis-regulatory sequence evolution across Y chromosomes has been challenging due to the limited number of reference assemblies. The threespine stickleback (Gasterosteus aculeatus) Y chromosome is an excellent model to identify how regulatory mutations accumulate on Y chromosomes due to its intermediate state of divergence from the X chromosome. A large number of Y-linked gametologs still exist across 3 differently aged evolutionary strata to test these hypotheses. We found that putative enhancer regions on the Y chromosome exhibited elevated substitution rates and decreased polymorphism when compared to nonfunctional sites, like intergenic regions and synonymous sites. This suggests that many cis-regulatory regions are under positive selection on the Y chromosome. This divergence was correlated with X-biased gametolog expression, indicating the loss of expression from the Y chromosome may be favored by selection. Our findings provide evidence that Y-linked cis-regulatory regions exhibit signs of positive selection quickly after the suppression of recombination and allow comparisons with recent theoretical models that suggest the rapid divergence of regulatory regions may be favored to mask deleterious mutations on the Y chromosome.

RP1 Dominant p.Ser740* Pathogenic Variant in 20 Knowingly Unrelated Families Affected by Rod-Cone Dystrophy: Potential Founder Effect in Western Sicily

Background and Objectives. Retinitis pigmentosa (RP) is the most common inherited rod-cone dystrophy (RCD), resulting in nyctalopia, progressive visual field, and visual acuity decay in the late stages. The autosomal dominant form (ADRP) accounts for about 20% of RPs. Among the over 30 genes found to date related to ADRP, RP1 pathogenic variants have been identified in 5-10% of cases. In a cohort of RCD patients from the Palermo province on the island of Sicily, we identified a prevalent nonsense variant in RP1, which was associated with ADRP. The objective of our study was to analyse the clinical and molecular data of this patient cohort and to evaluate the potential presence of a founder effect. Materials and Methods. From 2005 to January 2023, 84 probands originating from Western Sicily (Italy) with a diagnosis of RCD or RP and their relatives underwent deep phenotyping, which was performed in various Italian clinical institutions. Molecular characterisation of patients and familial segregation of pathogenic variants were carried out in different laboratories using Sanger and/or next-generation sequencing (NGS). Results. Among 84 probands with RCD/RP, we found 28 heterozygotes for the RP1 variant c.2219C>G, p.Ser740* ((NM_006269.2)*, which was therefore significantly prevalent in this patient cohort. After a careful interview process, we ascertained that some of these patients shared the same pedigree. Therefore, we were ultimately able to define 20 independent family groups with no traceable consanguinity. Lastly, analysis of clinical data showed, in our patients, that the p.Ser740* nonsense variant was often associated with a late-onset and relatively mild phenotype. Conclusions. The high prevalence of the p.Ser740* variant in ADRP patients from Western Sicily suggests the presence of a founder effect, which has useful implications for the molecular diagnosis of RCD in patients coming from this Italian region. This variant can be primarily searched for in RP-affected subjects displaying compatible modes of transmission and phenotypes, with an advantage in terms of the required costs and time for analysis. Moreover, given its high prevalence, the RP1 p.Ser740* variant could represent a potential candidate for the development of therapeutic strategies based on gene editing or translational read-through therapy for suppression of nonsense variants.

DNA polymerase ε and δ variants drive mutagenesis in polypurine tracts in human tumors

Alterations in the exonuclease domain of DNA polymerase ε cause ultramutated cancers. These cancers accumulate AGA>ATA transversions; however, their genomic features beyond the trinucleotide motifs are obscure. We analyze the extended DNA context of ultramutation using whole-exome sequencing data from 524 endometrial and 395 colorectal tumors. We find that G>T transversions in POLE-mutant tumors predominantly affect sequences containing at least six consecutive purines, with a striking preference for certain positions within polypurine tracts. Using this signature, we develop a machine-learning classifier to identify tumors with hitherto unknown POLE drivers and validate two drivers, POLE-E978G and POLE-S461L, by functional assays in yeast. Unlike other pathogenic variants, the E978G substitution affects the polymerase domain of Pol ε. We further show that tumors with POLD1 drivers share the extended signature of POLE ultramutation. These findings expand the understanding of ultramutation mechanisms and highlight peculiar mutagenic properties of polypurine tracts in the human genome.

Methods of Detection and Mechanisms of Origin of Complex Structural Genome Variations

Based on classical karyotyping, structural genome variations (SVs) have generally been considered to be either "simple" (with one or two breakpoints) or "complex" (with more than two breakpoints). Studying the breakpoints of SVs at nucleotide resolution revealed additional, subtle structural variations, such that even "simple" SVs turned out to be "complex." Genome-wide sequencing methods, such as fosmid and paired-end mapping, short-read and long-read whole genome sequencing, and single-molecule optical mapping, also indicated that the number of SVs per individual was considerably larger than expected from karyotyping and high-resolution chromosomal array-based studies. Interestingly, SVs were detected in studies of cohorts of individuals without clinical phenotypes. The common denominator of all SVs appears to be a failure to accurately repair DNA double-strand breaks (DSBs) or to halt cell cycle progression if DSBs persist. This review discusses the various DSB response mechanisms during the mitotic cell cycle and during meiosis and their regulation. Emphasis is given to the molecular mechanisms involved in the formation of translocations, deletions, duplications, and inversions during or shortly after meiosis I. Recently, CRISPR-Cas9 studies have provided unexpected insights into the formation of translocations and chromothripsis by both breakage-fusion-bridge and micronucleus-dependent mechanisms.

A new, all-encompassing aetiology of type 1 diabetes

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

Induction of Translational Readthrough on Protein Tyrosine Phosphatases Targeted by Premature Termination Codon Mutations in Human Disease

Nonsense mutations generating premature termination codons (PTCs) in various genes are frequently associated with somatic cancer and hereditary human diseases since PTCs commonly generate truncated proteins with defective or altered function. Induced translational readthrough during protein biosynthesis facilitates the incorporation of an amino acid at the position of a PTC, allowing the synthesis of a complete protein. This may evade the pathological effect of the PTC mutation and provide new therapeutic opportunities. Several protein tyrosine phosphatases (PTPs) genes are targeted by PTC in human disease, the tumor suppressor PTEN being the more prominent paradigm. Here, using PTEN and laforin as examples, two PTPs from the dual-specificity phosphatase subfamily, we describe methodologies to analyze in silico the distribution and frequency of pathogenic PTC in PTP genes. We also summarize laboratory protocols and technical notes to study the induced translational readthrough reconstitution of the synthesis of PTP targeted by PTC in association with disease in cellular models.

Small Indels in the Androgen Receptor Gene: Phenotype Implications and Mechanisms of Mutagenesis

CONTEXT

Despite high abundance of small indels in human genomes, their precise roles and underlying mechanisms of mutagenesis in Mendelian disorders require further investigation.

OBJECTIVE

To profile the distribution, functional implications, and mechanisms of small indels in the androgen receptor (AR) gene in individuals with androgen insensitivity syndrome (AIS).

METHODS

We conducted a systematic review of previously reported indels within the coding region of the AR gene, including 3 novel indels. Distribution throughout the AR coding region was examined and compared with genomic population data. Additionally, we assessed their impact on the AIS phenotype and investigated potential mechanisms driving their occurrence.

RESULTS

A total of 82 indels in AIS were included. Notably, all frameshift indels exhibited complete AIS. The distribution of indels across the AR gene showed a predominance in the N-terminal domain, most leading to frameshift mutations. Small deletions accounted for 59.7%. Most indels occurred in nonrepetitive sequences, with 15.8% situated within triplet regions. Gene burden analysis demonstrated significant enrichment of frameshift indels in AIS compared with controls (P < .00001), and deletions were overrepresented in AIS (P < .00001).

CONCLUSION

Our findings underscore a robust genotype-phenotype relationship regarding small indels in the AR gene in AIS, with a vast majority presenting complete AIS. Triplet regions and homopolymeric runs emerged as prone loci for small indels within the AR. Most were frameshift indels, with polymerase slippage potentially explaining half of AR indel occurrences. Complex frameshift indels exhibited association with palindromic runs. These discoveries advance understanding of the genetic basis of AIS and shed light on potential mechanisms underlying pathogenic small indel events.

Neoantigen-targeted TCR-engineered T cell immunotherapy: current advances and challenges

Adoptive cell therapy using T cell receptor-engineered T cells (TCR-T) is a promising approach for cancer therapy with an expectation of no significant side effects. In the human body, mature T cells are armed with an incredible diversity of T cell receptors (TCRs) that theoretically react to the variety of random mutations generated by tumor cells. The outcomes, however, of current clinical trials using TCR-T cell therapies are not very successful especially involving solid tumors. The therapy still faces numerous challenges in the efficient screening of tumor-specific antigens and their cognate TCRs. In this review, we first introduce TCR structure-based antigen recognition and signaling, then describe recent advances in neoantigens and their specific TCR screening technologies, and finally summarize ongoing clinical trials of TCR-T therapies against neoantigens. More importantly, we also present the current challenges of TCR-T cell-based immunotherapies, e.g., the safety of viral vectors, the mismatch of T cell receptor, the impediment of suppressive tumor microenvironment. Finally, we highlight new insights and directions for personalized TCR-T therapy.

Mutational hotspots lead to robust but suboptimal adaptive outcomes in certain environments

The observed mutational spectrum of adaptive outcomes can be constrained by many factors. For example, mutational biases can narrow the observed spectrum by increasing the rate of mutation at isolated sites in the genome. In contrast, complex environments can shift the observed spectrum by defining fitness consequences of mutational routes. We investigate the impact of different nutrient environments on the evolution of motility in Pseudomonas fluorescens Pf0-2x (an engineered non-motile derivative of Pf0-1) in the presence and absence of a strong mutational hotspot. Previous work has shown that this mutational hotspot can be built and broken via six silent mutations, which provide rapid access to a mutation that rescues swimming motility and confers the strongest swimming phenotype in specific environments. Here, we evolved a hotspot and non-hotspot variant strain of Pf0-2x for motility under nutrient-rich (LB) and nutrient-limiting (M9) environmental conditions. We observed the hotspot strain consistently evolved faster across all environmental conditions and its mutational spectrum was robust to environmental differences. However, the non-hotspot strain had a distinct mutational spectrum that changed depending on the nutrient environment. Interestingly, while alternative adaptive mutations in nutrient-rich environments were equal to, or less effective than, the hotspot mutation, the majority of these mutations in nutrient-limited conditions produced superior swimmers. Our competition experiments mirrored these findings, underscoring the role of environment in defining both the mutational spectrum and the associated phenotype strength. This indicates that while mutational hotspots working in concert with natural selection can speed up access to robust adaptive mutations (which can provide a competitive advantage in evolving populations), they can limit exploration of the mutational landscape, restricting access to potentially stronger phenotypes in specific environments.

Localizing unmapped sequences with families to validate the Telomere-to-Telomere assembly and identify new hotspots for genetic diversity

Although it is ubiquitous in genomics, the current human reference genome (GRCh38) is incomplete: It is missing large sections of heterochromatic sequence, and as a singular, linear reference genome, it does not represent the full spectrum of human genetic diversity. To characterize gaps in GRCh38 and human genetic diversity, we developed an algorithm for sequence location approximation using nuclear families (ASLAN) to identify the region of origin of reads that do not align to GRCh38. Using unmapped reads and variant calls from whole-genome sequences (WGSs), ASLAN uses a maximum likelihood model to identify the most likely region of the genome that a subsequence belongs to given the distribution of the subsequence in the unmapped reads and phasings of families. Validating ASLAN on synthetic data and on reads from the alternative haplotypes in the decoy genome, ASLAN localizes >90% of 100-bp sequences with >92% accuracy and ∼1 Mb of resolution. We then ran ASLAN on 100-mers from unmapped reads from WGS from more than 700 families, and compared ASLAN localizations to alignment of the 100-mers to the recently released T2T-CHM13 assembly. We found that many unmapped reads in GRCh38 originate from telomeres and centromeres that are gaps in GRCh38. ASLAN localizations are in high concordance with T2T-CHM13 alignments, except in the centromeres of the acrocentric chromosomes. Comparing ASLAN localizations and T2T-CHM13 alignments, we identified sequences missing from T2T-CHM13 or sequences with high divergence from their aligned region in T2T-CHM13, highlighting new hotspots for genetic diversity.

Human LINE-1 retrotransposons: impacts on the genome and regulation by host factors

Genome sequencing revealed that nearly half of the human genome is comprised of transposable elements. Although most of these elements have been rendered inactive due to mutations, full-length intact long interspersed element-1 (LINE-1 or L1) copies retain the ability to mobilize through RNA intermediates by a so-called "copy-and-paste" mechanism, termed retrotransposition. L1 is the only known autonomous mobile genetic element in the genome, and its retrotransposition contributes to inter- or intra-individual genetic variation within the human population. However, L1 retrotransposition also poses a threat to genome integrity due to gene disruption and chromosomal instability. Moreover, recent studies suggest that aberrant L1 expression can impact human health by causing diseases such as cancer and chronic inflammation that might lead to autoimmune disorders. To counteract these adverse effects, the host cells have evolved multiple layers of defense mechanisms at the epigenetic, RNA and protein levels. Intriguingly, several host factors have also been reported to facilitate L1 retrotransposition, suggesting that there is competition between negative and positive regulation of L1 by host factors. Here, we summarize the known host proteins that regulate L1 activity at different stages of the replication cycle and discuss how these factors modulate disease-associated phenotypes caused by L1.

Sequence dependencies and mutation rates of localized mutational processes in cancer

BACKGROUND

Cancer mutations accumulate through replication errors and DNA damage coupled with incomplete repair. Individual mutational processes often show nucleotide sequence and functional region preferences. As a result, some sequence contexts mutate at much higher rates than others, with additional variation found between functional regions. Mutational hotspots, with recurrent mutations across cancer samples, represent genomic positions with elevated mutation rates, often caused by highly localized mutational processes.

METHODS

We count the 11-mer genomic sequences across the genome, and using the PCAWG set of 2583 pan-cancer whole genomes, we associate 11-mers with mutational signatures, hotspots of single nucleotide variants, and specific genomic regions. We evaluate the mutation rates of individual and combined sets of 11-mers and derive mutational sequence motifs.

RESULTS

We show that hotspots generally identify highly mutable sequence contexts. Using these, we show that some mutational signatures are enriched in hotspot sequence contexts, corresponding to well-defined sequence preferences for the underlying localized mutational processes. This includes signature 17b (of unknown etiology) and signatures 62 (POLE deficiency), 7a (UV), and 72 (linked to lymphomas). In some cases, the mutation rate and sequence preference increase further when focusing on certain genomic regions, such as signature 62 in transcribed regions, where the mutation rate is increased up to 9-folds over cancer type and mutational signature average.

CONCLUSIONS

We summarize our findings in a catalog of localized mutational processes, their sequence preferences, and their estimated mutation rates.

Structural underpinnings of mutation rate variations in the human genome

Single nucleotide mutation rates have critical implications for human evolution and genetic diseases. Importantly, the rates vary substantially across the genome and the principles underlying such variations remain poorly understood. A recent model explained much of this variation by considering higher-order nucleotide interactions in the 7-mer sequence context around mutated nucleotides. This model's success implicates a connection between DNA shape and mutation rates. DNA shape, i.e. structural properties like helical twist and tilt, is known to capture interactions between nucleotides within a local context. Thus, we hypothesized that changes in DNA shape features at and around mutated positions can explain mutation rate variations in the human genome. Indeed, DNA shape-based models of mutation rates showed similar or improved performance over current nucleotide sequence-based models. These models accurately characterized mutation hotspots in the human genome and revealed the shape features whose interactions underlie mutation rate variations. DNA shape also impacts mutation rates within putative functional regions like transcription factor binding sites where we find a strong association between DNA shape and position-specific mutation rates. This work demonstrates the structural underpinnings of nucleotide mutations in the human genome and lays the groundwork for future models of genetic variations to incorporate DNA shape.

HATCHet2: clone- and haplotype-specific copy number inference from bulk tumor sequencing data

Multi-region DNA sequencing of primary tumors and metastases from individual patients helps identify somatic aberrations driving cancer development. However, most methods to infer copy-number aberrations (CNAs) analyze individual samples. We introduce HATCHet2 to identify haplotype- and clone-specific CNAs simultaneously from multiple bulk samples. HATCHet2 introduces a novel statistic, the mirrored haplotype B-allele frequency (mhBAF), to identify mirrored-subclonal CNAs having different numbers of copies of parental haplotypes in different tumor clones. HATCHet2 also has high accuracy in identifying focal CNAs and extends the earlier HATCHet method in several directions. We demonstrate HATCHet2's improved accuracy using simulations and a single-cell sequencing dataset. HATCHet2 analysis of 50 prostate cancer samples from 10 patients reveals previously-unreported mirrored-subclonal CNAs affecting cancer genes.

De novo retinoic acid receptor beta (RARB) variant associated with microphthalmia and dystonia

BACKGROUND

Definition of the individual genotypes that cause a Mendelian phenotype is of great importance both to clinical diagnostics and disease characterization. Heterozygous de novo gain-of-function missense variants in RARB are associated with syndromic microphthalmia 12 (MCOPS12), a developmental disorder characterized by eye malformations and variable involvement of other organs. A subset of patients were described with poorly delineated movement disorders. Additionally, RARB bi-allelic loss-of-function variants, inherited from asymptomatic heterozygous carrier parents, have been found in a recessive family with four MCOPS12-affected members.

PATIENT/METHODS

We used trio whole-exome sequencing to explore the molecular basis of disease in an individual with congenital eye abnormality and movement disorder. All patients with reported RARB variants were reviewed.

RESULTS

We report on identification of a heterozygous de novo RARB nonsense variant in a girl with microphthalmia and progressive generalized dystonia. Public database entries indicate that the de novo variant is recurrently present in clinically affected subjects but a literature report has not yet been available.

CONCLUSIONS

We provide the first detailed evidence for a role of dominant RARB truncating alterations in congenital eye-brain disease, expanding the spectrum of MCOPS12-associated mutations. Considered together with the published family with bi-allelic variants, the data suggest manifestation and non-manifestation of disease in relation to almost identical RARB loss-of-function variations, an apparent paradox that is seen in a growing number of human genetic conditions associated with both recessive and dominant inheritance patterns.

Preleukemic Fusion Genes Induced via Ionizing Radiation

Although the prevalence of leukemia is increasing, the agents responsible for this increase are not definitely known. While ionizing radiation (IR) was classified as a group one carcinogen by the IARC, the IR-induced cancers, including leukemia, are indistinguishable from those that are caused by other factors, so the risk estimation relies on epidemiological data. Several epidemiological studies on atomic bomb survivors and persons undergoing IR exposure during medical investigations or radiotherapy showed an association between radiation and leukemia. IR is also known to induce chromosomal translocations. Specific chromosomal translocations resulting in preleukemic fusion genes (PFGs) are generally accepted to be the first hit in the onset of many leukemias. Several studies indicated that incidence of PFGs in healthy newborns is up to 100-times higher than childhood leukemia with the same chromosomal aberrations. Because of this fact, it has been suggested that PFGs are not able to induce leukemia alone, but secondary mutations are necessary. PFGs also have to occur in specific cell populations of hematopoetic stem cells with higher leukemogenic potential. In this review, we describe the connection between IR, PFGs, and cancer, focusing on recurrent PFGs where an association with IR has been established.

DNA repair mechanisms that promote insertion-deletion events during immunoglobulin gene diversification

Insertions and deletions (indels) are low-frequency deleterious genomic DNA alterations. Despite their rarity, indels are common, and insertions leading to long complementarity-determining region 3 (CDR3) are vital for antigen-binding functions in broadly neutralizing and polyreactive antibodies targeting viruses. Because of challenges in detecting indels, the mechanism that generates indels during immunoglobulin diversification processes remains poorly understood. We carried out ultra-deep profiling of indels and systematically dissected the underlying mechanisms using passenger-immunoglobulin mouse models. We found that activation-induced cytidine deaminase-dependent ±1-base pair (bp) indels are the most prevalent indel events, biasing deleterious outcomes, whereas longer in-frame indels, especially insertions that can extend the CDR3 length, are rare outcomes. The ±1-bp indels are channeled by base excision repair, but longer indels require additional DNA-processing factors. Ectopic expression of a DNA exonuclease or perturbation of the balance of DNA polymerases can increase the frequency of longer indels, thus paving the way for models that can generate antibodies with long CDR3. Our study reveals the mechanisms that generate beneficial and deleterious indels during the process of antibody somatic hypermutation and has implications in understanding the detrimental genomic alterations in various conditions, including tumorigenesis.

Pangenome obtained by long-read sequencing of 11 genomes reveal hidden functional structural variants in pigs

Long-read sequencing (LRS) facilitates both the genome assembly and the discovery of structural variants (SVs). Here, we built a graph-based pig pangenome by incorporating 11 LRS genomes with an average of 94.01% BUSCO completeness score, revealing 206-Mb novel sequences. We discovered 183,352 nonredundant SVs (63% novel), representing 12.12% of the reference genome. By genotyping SVs in an additional 196 short-read sequencing samples, we identified thousands of population stratified SVs. Particularly, we detected 7,568 Tibetan specific SVs, some of which demonstrate significant population differentiation between Tibetan and low-altitude pigs, which might be associated with the high-altitude hypoxia adaptation in Tibetan pigs. Further integrating functional genomic data, the most promising candidate genes within the SVs that might contribute to the high-altitude hypoxia adaptation were discovered. Overall, our study generates a benchmark pangenome resource for illustrating the important roles of SVs in adaptive evolution, domestication, and genetic improvement of agronomic traits in pigs.

HotSPOT: A Computational Tool to Design Targeted Sequencing Panels to Assess Early Photocarcinogenesis

Mutations found in skin are acquired in specific patterns, clustering around mutation-prone genomic locations. The most mutation-prone genomic areas, mutation hotspots, first induce the growth of small cell clones in healthy skin. Mutations accumulate over time, and clones with driver mutations may give rise to skin cancer. Early mutation accumulation is a crucial first step in photocarcinogenesis. Therefore, a sufficient understanding of the process may help predict disease onset and identify avenues for skin cancer prevention. Early epidermal mutation profiles are typically established using high-depth targeted next-generation sequencing. However, there is currently a lack of tools for designing custom panels to capture mutation-enriched genomic regions efficiently. To address this issue, we created a computational algorithm that implements a pseudo-exhaustive approach to identify the best genomic areas to target. We benchmarked the current algorithm in three independent mutation datasets of human epidermal samples. Compared to the sequencing panel designs originally used in these publications, the mutation capture efficacy (number of mutations/base pairs sequenced) of our designed panel improved 9.6-12.1-fold. Mutation burden in the chronically sun-exposed and intermittently sun-exposed normal epidermis was measured within genomic regions identified by hotSPOT based on cutaneous squamous cell carcinoma (cSCC) mutation patterns. We found a significant increase in mutation capture efficacy and mutation burden in cSCC hotspots in chronically sun-exposed vs. intermittently sun-exposed epidermis (p < 0.0001). Our results show that our hotSPOT web application provides a publicly available resource for researchers to design custom panels, enabling efficient detection of somatic mutations in clinically normal tissues and other similar targeted sequencing studies. Moreover, hotSPOT also enables the comparison of mutation burden between normal tissues and cancer.

Hereditary Diffuse Gastric Cancer

Hereditary diffuse gastric cancer (HDGC) is a dominantly inherited cancer syndrome characterized by a high incidence of diffuse gastric cancer (DGC) and lobular breast cancer (LBC). HDGC is caused by germline mutations in 2 genes involved in the epithelial adherens junction complex, CDH1 and CTNNA1. We discuss the genetics of HDGC and the variability of its clinical phenotype, in particular the variable penetrance of advanced DGC and LBC, both within and between families. We review the pathology of the disease, the mechanism of tumor initiation, and its natural history. Finally, we describe current best practice for the clinical management of HDGC, including emerging genetic testing criteria for the identification of new families, methods for endoscopic surveillance, the complications associated with prophylactic surgery, postoperative quality of life, and the emerging field of HDGC chemoprevention.

A Comprehensive Investigation of Genomic Variants in Prostate Cancer Reveals 30 Putative Regulatory Variants

Prostate cancer (PC) is the most frequently diagnosed non-skin cancer in the world. Previous studies have shown that genomic alterations represent the most common mechanism for molecular alterations responsible for the development and progression of PC. This highlights the importance of identifying functional genomic variants for early detection in high-risk PC individuals. Great efforts have been made to identify common protein-coding genetic variations; however, the impact of non-coding variations, including regulatory genetic variants, is not well understood. Identification of these variants and the underlying target genes will be a key step in improving the detection and treatment of PC. To gain an understanding of the functional impact of genetic variants, and in particular, regulatory variants in PC, we developed an integrative pipeline (AGV) that uses whole genome/exome sequences, GWAS SNPs, chromosome conformation capture data, and ChIP-Seq signals to investigate the potential impact of genomic variants on the underlying target genes in PC. We identified 646 putative regulatory variants, of which 30 significantly altered the expression of at least one protein-coding gene. Our analysis of chromatin interactions data (Hi-C) revealed that the 30 putative regulatory variants could affect 131 coding and non-coding genes. Interestingly, our study identified the 131 protein-coding genes that are involved in disease-related pathways, including Reactome and MSigDB, for most of which targeted treatment options are currently available. Notably, our analysis revealed several non-coding RNAs, including RP11-136K7.2 and RAMP2-AS1, as potential enhancer elements of the protein-coding genes CDH12 and EZH1, respectively. Our results provide a comprehensive map of genomic variants in PC and reveal their potential contribution to prostate cancer progression and development.

Clinical and genetic characteristics of Tunisian children with infantile nephropathic cystinosis

BACKGROUND

Nephropathic cystinosis is an autosomal recessive disease caused by a mutation in the CTNS gene which encodes cystinosin, a lysosomal cystine transporter. The spectrum of mutations in the CTNS gene is not well defined in the North African population. Here, we investigated twelve patients with nephropathic cystinosis belonging to eight Tunisian families in order to analyze the clinical and genetic characteristics of Tunisian children with infantile nephropathic cystinosis.

METHODS

Clinical data were collected retrospectively. Molecular analysis of the CTNS gene was performed by Sanger sequencing.

RESULTS

We describe a new splicing mutation c.971-1G > C in the homozygous state in 6/12 patients which seems to be a founder mutation. The reported deletion of 23nt c.771_793 Del (p.Gly258Serfs*30) was detected in a homozygous state in one patient and in a heterozygous compound state with the c.971-1G > C mutation in 3/12 patients. Two of 12 patients have a deletion of exons 4 and 5 of the CTNS gene. None of our patients had the most common 57-kb deletion.

CONCLUSIONS

The mutational spectrum in the Tunisian population is different from previously described populations. Thus, a molecular diagnostic strategy must be implemented in Tunisia, by targeting as a priority the common mutations described in this country. Such a strategy will allow a cost-effective diagnosis confirmation as well as early administration of treatment with oral cysteamine. A higher resolution version of the Graphical abstract is available as Supplementary information.

Chromatin accessibility: methods, mechanisms, and biological insights

Access to DNA is a prerequisite to the execution of essential cellular processes that include transcription, replication, chromosomal segregation, and DNA repair. How the proteins that regulate these processes function in the context of chromatin and its dynamic architectures is an intensive field of study. Over the past decade, genome-wide assays and new imaging approaches have enabled a greater understanding of how access to the genome is regulated by nucleosomes and associated proteins. Additional mechanisms that may control DNA accessibility in vivo include chromatin compaction and phase separation - processes that are beginning to be understood. Here, we review the ongoing development of accessibility measurements, we summarize the different molecular and structural mechanisms that shape the accessibility landscape, and we detail the many important biological functions that are linked to chromatin accessibility.

A (GCC) repeat in SBF1 reveals a novel biological phenomenon in human and links to late onset neurocognitive disorder

The human SBF1 (SET binding factor 1) gene, alternatively known as MTMR5, is predominantly expressed in the brain, and its epigenetic dysregulation is linked to late-onset neurocognitive disorders (NCDs), such as Alzheimer’s disease. This gene contains a (GCC)-repeat at the interval between + 1 and + 60 of the transcription start site (SBF1-202 ENST00000380817.8). We sequenced the SBF1 (GCC)-repeat in a sample of 542 Iranian individuals, consisting of late-onset NCDs (N = 260) and controls (N = 282). While multiple alleles were detected at this locus, the 8 and 9 repeats were predominantly abundant, forming > 95% of the allele pool across the two groups. Among a number of anomalies, the allele distribution was significantly different in the NCD group versus controls (Fisher’s exact p = 0.006), primarily as a result of enrichment of the 8-repeat in the former. The genotype distribution departed from the Hardy–Weinberg principle in both groups (p < 0.001), and was significantly different between the two groups (Fisher’s exact p = 0.001). We detected significantly low frequency of the 8/9 genotype in both groups, higher frequency of this genotype in the NCD group, and reverse order of 8/8 versus 9/9 genotypes in the NCD group versus controls. Biased heterozygous/heterozygous ratios were also detected for the 6/8 versus 6/9 genotypes (in favor of 6/8) across the human samples studied (Fisher’s exact p = 0.0001). Bioinformatics studies revealed that the number of (GCC)-repeats may change the RNA secondary structure and interaction sites at least across human exon 1. This STR was specifically expanded beyond 2-repeats in primates. In conclusion, we report indication of a novel biological phenomenon, in which there is selection against certain heterozygous genotypes at a STR locus in human. We also report different allele and genotype distribution at this STR locus in late-onset NCD versus controls. In view of the location of this STR in the 5′ untranslated region, RNA/RNA or RNA/DNA heterodimer formation of the involved genotypes and alternative RNA processing and/or translation should be considered.

Individual Genetic Heterogeneity

Genetic variation has been widely covered in literature, however, not from the perspective of an individual in any species. Here, a synthesis of genetic concepts and variations relevant for individual genetic constitution is provided. All the different levels of genetic information and variation are covered, ranging from whether an organism is unmixed or hybrid, has variations in genome, chromosomes, and more locally in DNA regions, to epigenetic variants or alterations in selfish genetic elements. Genetic constitution and heterogeneity of microbiota are highly relevant for health and wellbeing of an individual. Mutation rates vary widely for variation types, e.g., due to the sequence context. Genetic information guides numerous aspects in organisms. Types of inheritance, whether Mendelian or non-Mendelian, zygosity, sexual reproduction, and sex determination are covered. Functions of DNA and functional effects of variations are introduced, along with mechanism that reduce and modulate functional effects, including TARAR countermeasures and intraindividual genetic conflict. TARAR countermeasures for tolerance, avoidance, repair, attenuation, and resistance are essential for life, integrity of genetic information, and gene expression. The genetic composition, effects of variations, and their expression are considered also in diseases and personalized medicine. The text synthesizes knowledge and insight on individual genetic heterogeneity and organizes and systematizes the central concepts.

Maintenance of quantitative genetic variance in complex, multi-trait phenotypes: The contribution of rare, large effect variants in two Drosophila species

The interaction of evolutionary processes to determine quantitative genetic variation has implications for contemporary and future phenotypic evolution, as well as for our ability to detect causal genetic variants. While theoretical studies have provided robust predictions to discriminate among competing models, empirical assessment of these has been limited. In particular, theory highlights the importance of pleiotropy in resolving observations of selection and mutation, but empirical investigations have typically been limited to few traits. Here, we applied high-dimensional Bayesian Sparse Factor Genetic modeling to gene expression datasets in 2 species, Drosophila melanogaster and Drosophila serrata, to explore the distributions of genetic variance across high-dimensional phenotypic space. Surprisingly, most of the heritable trait covariation was due to few lines (genotypes) with extreme [>3 interquartile ranges (IQR) from the median] values. Intriguingly, while genotypes extreme for a multivariate factor also tended to have a higher proportion of individual traits that were extreme, we also observed genotypes that were extreme for multivariate factors but not for any individual trait. We observed other consistent differences between heritable multivariate factors with outlier lines vs those factors without extreme values, including differences in gene functions. We use these observations to identify further data required to advance our understanding of the evolutionary dynamics and nature of standing genetic variation for quantitative traits.