Genomic Pipeline for Analysis of Mutational Events in Bacteria.

Unveiling the strategies of bacterial adaptation to stress constitute a challenging area of research. The understanding of mechanisms governing emergence of resistance to antimicrobials is of particular importance regarding the increasing threat of antibiotic resistance on public health worldwide. In the last decades, the fast democratization of sequencing technologies along with the development of dedicated bioinformatical tools to process data offered new opportunities to characterize genomic variations underlying bacterial adaptation. Thereby, research teams have now the possibility to dive deeper in the deciphering of bacterial adaptive mechanisms through the identification of specific genetic targets mediating survival to stress. In this chapter, we proposed a step-by-step bioinformatical pipeline enabling the identification of mutational events underlying biocidal stress adaptation associated with antimicrobial resistance development using Escherichia marmotae as an illustrative model.

A novel stoploss mutation CYB5R3 c.906A>G(p.*302Trpext*42) involved in the pathogenesis of hereditary methemoglobinemia.

Recessive congenital methemoglobinemia (RCM) is a hereditary autosomal disorder with an extremely low incidence rate. Here, we report a case of methemoglobinemia type I in a patient with congenital persistent cyanosis. The condition was attributed to a novel compound heterozygous mutation in CYB5R3, characterized by elevated methemoglobin levels (13.4 % of total hemoglobin) and undetectable NADH cytochrome b5 reductase (CYB5R3) activity. Whole-exome sequencing (WES) revealed two heterozygous mutations in CYB5R3: a previously reported pathogenic missense mutation c.611G>A(p.Cys204Tyr) inherited from the father, and a novel stop codon mutation c.906A>G(p.*302Trpext*42) from the mother, the latter mutation assessed as likely pathogenic according to ACMG guidelines. In cells overexpressing the CYB5R3 c.906A>G mutant construct, the CYB5R3 mRNA level was significantly lower than in cells overexpressing the wild-type (WT) CYB5R3 construct. However, there was no significant difference in protein expression levels between the mutant and WT constructs. Notably, an additional protein band of approximately 55 kDa was detected in the mutant cells. Immunofluorescence localization showed that, compared to wild-type CYB5R3, the subcellular localization of the CYB5R3 p.*302Trpext*42 mutant protein did not show significant changes and remained distributed in the endoplasmic reticulum and mitochondria. However, the c.906A>G(p.*302Trpext*42) mutation resulted in increased intracellular reactive oxygen species (ROS) levels and decreased NAD+/NADH ratio, suggesting impaired CYB5R3 function and implicating this novel mutation as likely pathogenic.

SCN1A intronic variants impact on Nav1.1 protein expression and sodium channel function, and associated with epilepsy phenotypic severity.

High-throughput sequencing has identified numerous intronic variants in the SCN1A gene in epilepsy patients. Abnormal mRNA splicing caused by these variants can lead to significant phenotypic differences, but the mechanisms of epileptogenicity and phenotypic differences remain unknown. Two variants, c.4853-1 G>C and c.4853-25 T>A, were identified in intron 25 of SCN1A, which were associated with severe Dravet syndrome (DS) and mild focal epilepsy with febrile seizures plus (FEFS+), respectively. The impact of these variants on protein expression, electrophysiological properties of sodium channels and their correlation with epilepsy severity was investigated through plasmid construction and transfection based on the aberrant spliced mRNA. We found that the expression of truncated mutant proteins was significantly reduced on the cell membrane, and retained in the cytoplasmic endoplasmic reticulum. The mutants caused a decrease in current density, voltage sensitivity, and an increased vulnerability of channel, leading to a partial impairment of sodium channel function. Notably, the expression of DS-related mutant protein on the cell membrane was higher compared to that of FEFS+-related mutant, whereas the sodium channel function impairment caused by DS-related mutant was comparatively milder than that caused by FEFS+-related mutant. Our study suggests that differences in protein expression levels and altered electrophysiological properties of sodium channels play important roles in the manifestation of diverse epileptic phenotypes. The presence of intronic splice site variants may result in severe phenotypes due to the dominant-negative effects, whereas non-canonical splice site variants leading to haploinsufficiency could potentially cause milder phenotypes.

First Brazilian Case Report of Unrelated Patients with Identical ISG15 Mutation.

BACKGROUND: ISG15 deficiency is a mixed syndrome of Mendelian susceptibility to mycobacterial infections (MSMD), a rare inherited condition characterized primarily by recurrent infections from low-virulence mycobacteria and monogenic type I interferonopathy. OBJECTIVE: To characterize the laboratory and molecular features of two patients from different families affected by the same ISG15 variant. METHODS: We began with clinical characterization and investigation, assessed IL-12/IFN-γ production, performed genetic characterization through WES and Sanger sequencing, conducted an in silico molecular analysis of the genetic ISG15 variant's protein impact, and utilized RNAseq for transcriptome analysis to understand pathway impacts on ISG15-deficient subjects from unrelated families. RESULTS: A mutation in the ISG15 gene was identified, affecting two patients treated in different hospitals and cities in Brazil (Fortaleza and Sao Paulo), who are also members of unrelated families. Both patients showed low IFN-γ production when stimulated with BCG or BCG + IL-12. ISG15 deficiency presented with two distinct clinical phenotypes: infectious and neurological. It was identified that both patients are homozygous for the variant (c.83 T > A). Furthermore, it was observed that the mutant protein p.L28Q results in an unstable protein with increased flexibility (ΔΔG: -2.400 kcal/mol). Transcriptome analysis revealed 1321 differentially expressed genes, with significant upregulation in interferon pathways, showing higher expression in patients compared to controls. CONCLUSION: This study describes the first reported cases in Brazil of two unrelated patients with the same ISG15 mutation c.83 T > A, exhibiting infectious features such as mycobacterial infections and systemic candidiasis, neurological findings, and skin lesions, without adverse reactions to the BCG vaccine. CLINICAL IMPLICATIONS: Reporting ISG15 gene mutations in Brazilian patients enhances understanding of genetic susceptibilities, guiding effective diagnostics and treatment. Identifying high-risk individuals aids clinical practices, genetic counseling, and influences public health policies. We have identified the first case in Brazil of the same ISG15 variant c.83 T > A that was identified in two unrelated patients with distinct clinical phenotypes, infectious and neurological.

Phylogenetic analysis, genetic diversity, and epidemiology of pigeon paramyxovirus type 1 in China.

Pigeon paramyxovirus type 1 (PPMV-1) poses significant economic challenges to the pigeon industry in China. However, information about the prevalence, genetic diversity, and epidemiology of PPMV-1 in China is still lacking. In this study, we isolated six strains of PPMV-1 from Hubei and Zhejiang provinces in 2022. All six isolates were found to belong to subgenotype VI.2.1.1.2.2. Five of them were identified as mesogenic and one as lentogenic. Multiple mutations were observed in the F and HN proteins of these isolates. Comprehensive analysis of global PPMV-1 strains highlighted the dominance of genotype VI, showing that VI.2.1.1.2.2 has been the dominant subgenotype since 2011. We also identified 36 host-specific amino acid substitutions that are unique to PPMV-1 in comparison to chicken-origin NDVs. The data reported here contribute to our understanding of the epidemiology, genetic diversity, and prevalence of PPMV-1 and serve as a valuable reference for the prevention and control of PPMV-1.

Differential modulation of mutant CALR and JAK2 V617F-driven oncogenesis by HLA genotype in myeloproliferative neoplasms.

It has been demonstrated previously that human leukocyte antigen class I (HLA-I) and class II (HLA-II) alleles may modulate JAK2 V617F and CALR mutation (CALRmut)-associated oncogenesis in myeloproliferative neoplasms (MPNs). However, the role of immunogenetic factors in MPNs remains underexplored. We aimed to investigate the potential involvement of HLA genes in CALRmut+ MPNs. High-resolution genotyping of HLA-I and -II loci was conducted in 42 CALRmut+ and 158 JAK2 V617F+ MPN patients and 1,083 healthy controls. A global analysis of the diversity of HLA-I genotypes revealed no significant differences between CALRmut+ patients and controls. However, one HLA-I allele (C*06:02) showed an inverse correlation with presence of CALR mutation. A meta-analysis across independent cohorts and healthy individuals from the 1000 Genomes Project confirmed an inverse correlation between the presentation capabilities of the HLA-I loci for JAK2 V617F and CALRmut-derived peptides in both patients and healthy individuals. scRNA-Seq analysis revealed low expression of TAP1 and CIITA genes in CALRmut+ hematopoietic stem and progenitor cells. In conclusion, the HLA-I genotype differentially restricts JAK2 V617F and CALRmut-driven oncogenesis potentially explaining the mutual exclusivity of the two mutations and differences in their presentation latency. These findings have practical implications for the development of neoantigen-based vaccines in MPNs.

Case Series: ATRX Variants in Four Patients with Metastatic Pheochromocytoma.

Few reports have highlighted the rare presence of somatic ATRX variants in clinically aggressive, metastatic pheochromocytoma/paraganglioma (PCC/PGL); however, none have addressed detailed clinical presentation (including biochemistry and imaging) and management of these patients. Here, we address these clinical features and management based on four PCC patients with somatic ATRX variants from our National Institutes of Health PCC/PGL cohort. A total of 192 patients underwent exome sequencing (germline, somatic, or both), and four males were found to have somatic ATRX variants (with additional somatic VHL and FH oncogenic variants in patients 2 and 4, respectively). Per-lesion and per-patient comparisons were performed among functional imaging scans performed at the NIH. Biochemical phenotype and response to systemic treatment were evaluated. This mini-series supports prior studies showing aggressive/metastatic PCC in patients with somatic ATRX variants, as all developed widespread metastatic disease. All four PCC patients presented with noradrenergic biochemical phenotype, and some with significant elevation in 3-methoxytyramine. 18F-FDOPA PET/CT was found to be the superior functional imaging modality, with 100% lesion detection rate when compared to that of 68Ga-DOTATATE, 18F-FDG, 18F-FDA, and 123I-MIBG scans. While patients did not respond to chemotherapy or tyrosine kinase inhibitors, they responded to targeted radiotherapy using high-specific-activity 131I-MIBG (Azedra®) or 177Lu-DOTATATE (Lutathera®).

Deep Sequencing and Phenotyping in an Australian Tuberous Sclerosis Complex "No Mutations Identified" Cohort.

Tuberous sclerosis complex (TSC) is a variable multisystem disorder. The "no mutations identified" (NMI) group are reportedly phenotypically milder than those with an identified molecular cause, and often have mosaic or intronic variants not detected by standard sequencing methods. METHODS: We describe the phenotypes in an Australian TSC NMI group (n = 18) and a molecular testing strategy implementable in a diagnostic laboratory. Massively parallel sequencing (MPS) of the whole genomic regions of TSC1 and TSC2 was performed using DNA extracted from multiple tissue samples per participant. RESULTS: Our study showed that the phenotype in TSC NMI individuals can be similar to those with heterozygous, particularly TSC1, variants. Although neurodevelopmental outcomes can be less severe, the number of organ systems involved was similar to the non-mosaic groups. A diagnostic yield of 72% (13/18) was achieved, with the majority (10/13) being mosaic variants and the remainder heterozygous variants missed on previous testing. CONCLUSION: Testing DNA from multiple tissue samples allowed for validation of otherwise discarded low-level mosaic variants and detection of mosaic variants by MPS without excessive cost or the need for specialised techniques. Implementing this approach in a diagnostic setting is viable and allows optimal clinical care of patients with NMI TSC.

Distinct landscape and clinical implications of therapy-related clonal hematopoiesis.

Therapy-related clonal hematopoiesis (t-CH) is defined as clonal hematopoiesis detected in individuals previously treated with chemotherapy and/or radiation therapy. With the increased use of genetic analysis in oncological care, the detection of t-CH among cancer patients is becoming increasingly common. t-CH arises through the selective bottleneck imposed by chemotherapies and potentially through direct mutagenesis from chemotherapies, resulting in a distinct mutational landscape enriched with mutations in DNA damage-response pathway genes such as TP53, PPM1D, and CHEK2. Emerging evidence sheds light on the mechanisms of t-CH development and potential strategies to mitigate its emergence. Due to its unique characteristics that predominantly affect cancer patients, t-CH has clinical implications distinct from those of CH in the general population. This Review discusses the potential mechanisms of t-CH development, its mutational landscape, mutant-drug relationships, and its clinical significance. We highlight the distinct nature of t-CH and call for intensified research in this field.

Clonal hematopoiesis and atherosclerosis.

Clonal hematopoiesis of indeterminate potential (CHIP) has emerged as a previously unrecognized, potent, age-related, and common risk factor for atherosclerosis. Somatic mutations in certain known leukemia driver genes give rise to clones of mutant cells in peripheral blood. The increased risk of developing hematologic malignancy does not, on its own, explain excess mortality in individuals with CHIP. Cardiovascular disease accounts for much of this gap. Experimental evidence supports the causality of certain CHIP mutations in accelerated atherosclerosis. CHIP due to mutations in different driver genes varies in their promotion of atherosclerotic events and in the region of augmented atherosclerotic involvement. For example, CHIP due to mutations in DNMT3a appears less atherogenic than CHIP that arises from TET2 or JAK2, forms of CHIP that incite inflammation. The recognition of certain CHIP mutations as promoters of atherosclerotic risk has opened new insights into understanding of the pathophysiology of this disease. The accentuated cardiovascular risk and involvement of distinct pathways of various forms of CHIP also inform novel approaches to allocation of targeted therapies, affording a step toward personalized medicine.

Clonal hematopoiesis and hematological malignancy.

Clonal hematopoiesis (CH), the expansion of hematopoietic stem cells and their progeny driven by somatic mutations in leukemia-associated genes, is a common phenomenon that rises in prevalence with advancing age to affect most people older than 70 years. CH remains subclinical in most carriers, but, in a minority, it progresses to a myeloid neoplasm, such as acute myeloid leukemia, myelodysplastic syndrome, or myeloproliferative neoplasm. Over the last decade, advances in our understanding of CH, its molecular landscape, and the risks associated with different driver gene mutations have culminated in recent developments that allow for a more precise estimation of myeloid neoplasia risk in CH carriers. In turn, this is leading to the development of translational and clinical programs to intercept and prevent CH from developing into myeloid neoplasia. Here, we give an overview of the spectrum of CH driver mutations, what is known about their pathophysiology, and how this informs the risk of incident myeloid malignancy.

Molecular chaperones: Guardians of tumor suppressor stability and function.

The term 'tumor suppressor' describes a widely diverse set of genes that are generally involved in the suppression of metastasis, but lead to tumorigenesis upon loss-of-function mutations. Despite the protein products of tumor suppressors exhibiting drastically different structures and functions, many share a common regulatory mechanism-they are molecular chaperone 'clients'. Clients of molecular chaperones depend on an intracellular network of chaperones and co-chaperones to maintain stability. Mutations of tumor suppressors that disrupt proper chaperoning prevent the cell from maintaining sufficient protein levels for physiological function. This review discusses the role of the molecular chaperones Hsp70 and Hsp90 in maintaining the stability and functional integrity of tumor suppressors. The contribution of cochaperones prefoldin, HOP, Aha1, p23, FNIP1/2 and Tsc1 as well as the chaperonin TRiC to tumor suppressor stability is also discussed. Genes implicated in renal cell carcinoma development-VHL, TSC1/2, and FLCN-will be used as examples to explore this concept, as well as how pathogenic mutations of tumor suppressors cause disease by disrupting protein chaperoning, maturation, and function.

Flagellar motility is mutagenic.

Flagella are highly complex rotary molecular machines that enable bacteria to not only migrate to optimal environments but also to promote range expansion, competitiveness, virulence, and antibiotic survival. Flagellar motility is an energy-demanding process, where the sum of its production (biosynthesis) and operation (rotation) costs has been estimated to total ~10% of the entire energy budget of an Escherichia coli cell. The acquisition of such a costly adaptation process is expected to secure short-term benefits by increasing competitiveness and survival, as well as long-term evolutionary fitness gains. While the role of flagellar motility in bacterial survival has been widely reported, its direct influence on the rate of evolution remains unclear. We show here that both production and operation costs contribute to elevated mutation rates. Our findings suggest that flagellar movement may be an important player in tuning the rate of bacterial evolution.

PPP3CB overexpression mediates EGFR TKI resistance in lung tumors via calcineurin/MEK/ERK signaling.

Despite initial high response rates to first-line EGFR TKI, all non-small-cell lung cancer (NSCLC) with EGFR-activating mutation will ultimately develop resistance to treatment. Identification of resistance mechanisms is critical to adapt treatment and improve patient outcomes. Here, we show that a PPP3CB transcript that encodes full-length catalytic subunit 2B of calcineurin accumulates in EGFR-mutant NSCLC cells with acquired resistance against different EGFR TKIs and in post-progression biopsies of NSCLC patients treated with EGFR TKIs. Neutralization of PPP3CB by siRNA or inactivation of calcineurin by cyclosporin A induces apoptosis in resistant cells treated with EGFR TKIs. Mechanistically, EGFR TKIs increase the cytosolic level of calcium and trigger activation of a calcineurin/MEK/ERK pathway that prevents apoptosis. Combining EGFR, calcineurin, and MEK inhibitors overcomes resistance to EGFR TKI in both in vitro and in vivo models. Our results identify PPP3CB overexpression as a new mechanism of acquired resistance to EGFR TKIs, and provide a promising therapeutic approach for NSCLC patients that progress under TKI treatment.

AVM: A Manually Curated Database of Aerosol-transmitted Virus Mutations, Human Diseases, and Drugs.

Aerosol-transmitted viruses possess strong infectivity and can spread over long distances, earning the difficult-to-control title. They cause various human diseases and pose serious threats to human health. Mutations can increase the transmissibility and virulence of the strains, reducing the protection provided by vaccines and weakening the efficacy of antiviral drugs. In this study, we established a manually curated database (termed AVM) to store information on aerosol-transmitted viral mutations (VMs). The current version of the AVM contains 42,041 VMs (including 2613 immune escape mutations), 45 clinical information datasets, and 407 drugs/antibodies/vaccines. Additionally, we recorded 88 human diseases associated with viruses and found that the same virus can target multiple organs in the body, leading to diverse diseases. Furthermore, the AVM database offers a straightforward user interface for browsing, retrieving, and downloading information. This database is a comprehensive resource that can provide timely and valuable information on the transmission, treatment, and diseases caused by aerosol-transmitted viruses (http://www.bio-bigdata.center/AVM).

A general multi-scale description of metastable adaptive motion across fitness valleys.

We consider a stochastic individual-based model of adaptive dynamics on a finite trait graph G = ( V , E ) . The evolution is driven by a linear birth rate, a density dependent logistic death rate and the possibility of mutations along the directed edges in E. We study the limit of small mutation rates for a simultaneously diverging population size. Closing the gap between Bovier et al. (Ann Appl Probab 29(6):3541-358, 2019) and Coquille et al. (Electron J Probab 26:1-37, 2021) we give a precise description of transitions between evolutionary stable conditions (ESC), where multiple mutations are needed to cross a valley in the fitness landscape. The system shows a metastable behaviour on several divergent time scales, corresponding to the widths of these fitness valleys. We develop the framework of a meta graph that is constituted of ESCs and possible metastable transitions between them. This allows for a concise description of the multi-scale jump chain arising from concatenating several jumps. Finally, for each of the various time scales, we prove the convergence of the population process to a Markov jump process visiting only ESCs of sufficiently high stability.

CFTR Exon 10 deleterious mutations in patients with congenital bilateral absence of vas deferens in a cohort of Pakistani patients.

Congenital bilateral absence of vas deferens (CBAVD) is a urological syndrome of Wolffian ducts and is responsible for male infertility and obstructive azoospermia. This study is designed to explore the integrity of exon 10 of CFTR and its role in male infertility in a cohort of CBVAD patients in Pakistan. Genomic DNA was extracted from 17 male patients with CBAVD having clinical symptoms, and 10 healthy controls via phenol-chloroform method. Exon 10 of the CFTR gene was amplified, using PCR with specific primers and DNA screening was done by Sanger sequencing. Sequencing results were analyzed using freeware Serial Cloner, SnapGene, BioEdit and FinchTV. Furthermore, bioinformatics tools were used to analyze the mutations and their impact on the protein function and stability. We have identified 4 mutations on exon 10 of CFTR in 6 out of 17 patients. Two of the mutations were missense variants V456A, K464E, and the other two were silent mutations G437G, S431S. The identified variant V456A was present in 4 of the studied patients. Whereas, the presence of K464E in our patients further weighs on the crucial importance for its strategic location to influence the gene function at post-transcriptional and protein level. Furthermore, Polyphen-2 and SIFT analyze the mutations as harmful and deleterious. The recurrence of V456A and tactically conserved locality of K464E are evidence of their potential role in CBAVD patients and in male infertility. The data can contribute in developing genetic testing and treatment of CBAVD.

Isoform-specific C-terminal phosphorylation drives autoinhibition of Casein kinase 1.

Casein kinase 1δ (CK1δ) controls essential biological processes including circadian rhythms and wingless-related integration site (Wnt) signaling, but how its activity is regulated is not well understood. CK1δ is inhibited by autophosphorylation of its intrinsically disordered C-terminal tail. Two CK1 splice variants, δ1 and δ2, are known to have very different effects on circadian rhythms. These variants differ only in the last 16 residues of the tail, referred to as the extreme C termini (XCT), but with marked changes in potential phosphorylation sites. Here, we test whether the XCT of these variants have different effects in autoinhibition of the kinase. Using NMR and hydrogen/deuterium exchange mass spectrometry, we show that the δ1 XCT is preferentially phosphorylated by the kinase and the δ1 tail makes more extensive interactions across the kinase domain. Mutation of δ1-specific XCT phosphorylation sites increases kinase activity both in vitro and in cells and leads to changes in the circadian period, similar to what is reported in vivo. Mechanistically, loss of the phosphorylation sites in XCT disrupts tail interaction with the kinase domain. δ1 autoinhibition relies on conserved anion-binding sites around the CK1 active site, demonstrating a common mode of product inhibition of CK1δ. These findings demonstrate how a phosphorylation cycle controls the activity of this essential kinase.