Detrimental Mutations Research Articles

PLS‐guided Mutant Recombination to Improve the Stability of Bovine Enterokinases Obtained by Directed Evolution

Activity and thermostability are critical yet challenging to improve simultaneously in enzymes. Using directed evolution, we previously identified bovine enterokinase (EKL) variants with enhanced soluble expression and thermal stability. Partial least squares (PLS) analysis of 321 EKL variants revealed the impact of individual mutations and identified neutral or detrimental mutations in top‐performing variants. Leveraging PLS rankings, we created new variants with fewer mutations and enhanced stability. Most original and PLS‐guided variants exhibited an activity‐stability trade‐off. However, two new triple‐ and quadruple‐mutants improved both activity and stability, surpassing the trade‐off limit. Recombining PLS‐guided mutations likely eliminated neutral or harmful mutations, enhancing stability. MD simulations linked residue‐specific dynamics with stability, pinpointing critical structural regions near aggregation‐prone areas. Our findings validate PLS as a potent strategy to enhance enzyme properties, complementing directed evolution.

On the evolution of chromosomal regions with high gene strand bias in bacteria.

On circular bacterial chromosomes, the majority of genes are coded on the leading strand. This gene strand bias (GSB) ranges from up to 85% in some Bacillota to a little more than 50% in other phyla. The factors determining the extent of the strand bias remain to be found. Here, we report that species in the phylum Gemmatimonadota share a unique chromosome architecture, distinct from neighboring phyla: in a conserved 600-kb region around the terminus of replication, almost all genes were located on the leading strands, while on the remaining part of the chromosome, the strand preference was more balanced. The high strand bias (HSB) region harbors the rRNA clusters, core, and highly expressed genes. Selective pressure for reduction of collisions with DNA replication to minimize detrimental mutations can explain the conservation of essential genes in this region. Repetitive and mobile elements are underrepresented, suggesting reduced recombination frequency by structural isolation from other parts of the chromosome. We propose that the HSB region forms a distinct chromosomal domain. Gemmatimonadota chromosomes evolved mainly by expansion through horizontal gene transfer and duplications outside of the ancient high strand bias region. In support of our hypothesis, we could further identify two Spiroplasma strains on a similar evolutionary path.IMPORTANCEOn bacterial chromosomes, a preferred location of genes on the leading strand has evolved to reduce conflicts between replication and transcription. Despite a vast body of research, the question why bacteria show large differences in their gene strand bias is still not solved. The discovery of "hybrid" chromosomes in different phyla, including Gemmatimonadota, in which a conserved high strand bias is found exclusively in a region at ter, points toward a role of nucleoid structure, additional to replication, in the evolution of strand preferences. A fine-grained structural analysis of the ever-increasing number of available bacterial genomes could help to better understand the forces that shape the sequential and spatial organization of the cell's information content.

Role of vosoritide drug on children's with achondroplasia

Achondroplasia is a genetic condition with an autosomal dominant inheritance pattern. It is caused by detrimental heterozygous FGFR3 gene mutations. Gain-of-function mutations in the FGFR3 gene impair chondrocyte differentiation and proliferation, which has a negative impact on the formation of new bone. Vosoritide, also known as voxzogo, is the first pharmacological intervention to receive approval for the treatment of achondroplasia. This therapeutic agent has been sanctioned for administration via at-home injections, which are to be performed under the supervision of a proficient nurse.

Exploring the effect of disease causing mutations in metal binding sites of human ARSA in metachromatic leukodystrophy.

The arylsulfatase A (ARSA) gene is observed to be deficient in patients with metachromatic leukodystrophy (MLD), a type of lysosomal storage disease. MLD is a severe neurodegenerative disorder characterized by an autosomal recessive inheritance pattern. This study aimed to map the most deleterious mutations at the metal binding sites of ARSA and the amino acids in proximity to the mutated positions. We utilized an array of computational tools, including PredictSNP, MAPP, PhD-SNP, PolyPhen-1, PolyPhen-2, SIFT, SNAP, and ConSurf, to identify the most detrimental mutations potentially implicated in MLD collected from UniProt, ClinVar, and HGMD. Two mutations, D29N and D30H, as being extremely deleterious based on assessments of pathogenicity, conservation, biophysical characteristics, and stability analysis. The D29 and D30 are located at the metal-interacting regions of ARSA and found to undergo post-translational modification, specifically phosphorylation. Henceforth, the in-depth effect of metal binding upon mutation was examined using molecular dynamics simulations (MDS) before and after phosphorylation. The MDS results exhibited high deviation for the D29N and D30H mutations in comparison to the native, and the same was confirmed by significant residue fluctuation and reduced compactness. These structural alterations suggest that such mutations may influence protein functionality, offering potential avenues for personalized therapeutic and providing a basis for potential mutation-specific treatments for severe MLD patients.

Analysis of Evolutionary Expression Relationships of Tea Plant Genes Based on the Identification of RCI2 Gene Members in Tea Plant (Camellia sinensis)

Plants in nature frequently encounter abiotic stresses, leading to constant mutation and evolution of plant genes as part of the adaptation process. Rare cold-inducible (RCI2) genes are a specific class of genes closely associated with abiotic stress. In our study, we identified 11 RCI2 genes in tea plants. Our evolutionary analysis revealed a high similarity between CsRCI2A and CsRCI2B with AtRCI2A and AtRCI2B. Furthermore, we found that CsRCI2B, CsRCI2D, CsRCI2I, and CsRCI2J were expressed in both drought and low-temperature conditions, as well as in different tissues. The transcript abundance of these genes varied significantly across different treatment times and tissues. However, these genes had fewer mutation sites, indicating that they underwent amino acid mutations during the evolutionary process after abiotic stress, while maintaining a stable gene structure without excessive mutation sites. Our results suggest that most of the genes in tea plants may degrade under the influence of the external environment. However, core genes, which do not undergo detrimental mutations, play a crucial role in tea plants after experiencing abiotic stress. This confirms that RCI2 genes in tea plants selectively undergo deleterious and beneficial mutations in response to abiotic stress, ultimately impacting gene expression.

CSIG-14. THE IMPACT OF CORTICOTROPHIN RELEASING HORMONE STRESS NEURO-HORMONES ON THE VIABILITY AND PROLIFERATION OF BRAIN TUMOURS

Abstract The majority of brain tumour sub-types are of glial origin. These include high-grade gliomas (HGGs), which lack effective targeted therapies and are thus invariably lethal. We have shown that stress leads to proliferation of brain glia. Such stress-induced glial cell division could provide opportunities for detrimental genetic mutations which increase the risk of oncogenesis. Identifying stress-induced glial plasticity could uncover molecular targets that directly influence glioma viability. Corticotropin-releasing hormone (CRH) stress neurohormones, including urocortin I-III (UCN I-III) are key mediators of stress signalling in the brain including neuroplastic processes. The aims of this project were to compare the expression CRH peptides and receptors in different types of brain tumours and their impact on cells viability and proliferation. Five different brain tumour cell lines were used (two biopsy-derived-in-house-established HGGs – UP007 and SEBTA023 – and three commercially available Medulloblastomas, MBs – UW402, CHLA-01-MED and CHLA-01R-MED). The gene expression profiles of the CRH family were accessed by qRT-PCR (n = 5) and the protein expression by Immunocytochemistry and flow cytometry (FC, n = 3). The effects of targeted drug modulation were assessed through viability (MTS, n = 3) and proliferation assays (Cyquant NF cell proliferation, n = 3). Statistical significance of qRT-PCRs and FC were calculated through a Kruskal-Wallis test, and for viability and proliferation assays through a two-way ANOVA approach. Values of P < 0.05 were accepted as significant. The qRT-PCRs revealed a higher expression of stress neuro-hormones mRNA in MBs compared to HGGs. Immunocytochemistry confirmed the protein expression of CRH and downstream receptors in all cell lines, while FC revealed a higher expression of CRH-Rs in HGGs than MBs. The pharmacological modulation of CRH-Rs induced contrasting effects on cellular viability and proliferation, suggesting patient-specific CRH-CRH-R pathway. This raises the prospect of repurposing available CRH-R ligands for targeted treatment of susceptible tumour types.

Reanalysis of mtDNA mutations of human primordial germ cells (PGCs) reveals NUMT contamination and suggests that selection in PGCs may be positive

The resilience of the mitochondrial genome (mtDNA) to a high mutational pressure depends, in part, on negative purifying selection in the germline. A paradigm in the field has been that such selection, at least in part, takes place in primordial germ cells (PGCs). Specifically, Floros et al. (Nature Cell Biology 20: 144–51) reported an increase in the synonymity of mtDNA mutations (a sign of purifying selection) between early-stage and late-stage PGCs. We re-analyzed Floros’ et al. data and determined that their mutational dataset was significantly contaminated with single nucleotide variants (SNVs) derived from a nuclear sequence of mtDNA origin (NUMT) located on chromosome 5. Contamination was caused by co-amplification of the NUMT sequence by cross-specific PCR primers. Importantly, when we removed NUMT-derived SNVs, the evidence of purifying selection was abolished. In addition to bulk PGCs, Floros et al. reported the analysis of single-cell late-stage PGCs, which were amplified with different sets of PCR primers that cannot amplify the NUMT sequence. Accordingly, there were no NUMT-derived SNVs among single PGC mutations. Interestingly, single PGC mutations show adecreaseof synonymity with increased intracellular mutant fraction. More specifically, nonsynonymous mutations show faster intracellular genetic drift towards higher mutant fraction than synonymous ones. This pattern is incompatible with predominantly negative selection. This suggests that germline selection of mtDNA mutations is a complex phenomenon and that the part of this process that takes place in PGCs may be predominantly positive. However counterintuitive, positive germline selection of detrimental mtDNA mutations has been reported previously andpotentially may be evolutionarily advantageous.

Genomic estimates of mutation and substitution rates contradict the evolutionary speed hypothesis of the latitudinal diversity gradient.

The latitudinal diversity gradient (LDG) refers to a decrease in biodiversity from the equator to the poles. The evolutionary speed hypothesis, backed by the metabolic theory of ecology, asserts that nucleotide mutation and substitution rates per site per year are higher and thereby speciation rates are higher at higher temperatures, generating the LDG. However, prior empirical investigations of the relationship between the temperature and mutation or substitution rate were based on a few genes and the results were mixed. We here revisit this relationship using genomic data. No significant correlation between the temperature and mutation rate is found in 13 prokaryotes or in 107 eukaryotes. An analysis of 234 diverse trios of bacterial taxa indicates that the synonymous substitution rate is not significantly associated with the growth temperature. The same data, however, reveal a significant negative association between the nonsynonymous substitution rate and temperature, which is explainable by a larger fraction of detrimental nonsynonymous mutations at higher temperatures due to a stronger demand for protein stability. We conclude that the evolutionary speed hypothesis of the LDG is unsupported by genomic data and advise that future mechanistic studies of the LDG should focus on other hypotheses.

Cold-Driven Hemoglobin Evolution in Antarctic Notothenioid Fishes Prior to Hemoglobin Gene Loss in White-Blooded Icefishes.

Expression of multiple hemoglobin isoforms with differing physiochemical properties likely helps species adapt to different environmental and physiological conditions. Antarctic notothenioid fishes inhabit the icy Southern Ocean and display fewer hemoglobin isoforms, each with less affinity for oxygen than temperate relatives. Reduced hemoglobin multiplicity was proposed to result from relaxed selective pressure in the cold, thermally stable, and highly oxygenated Antarctic waters. These conditions also permitted the survival and diversification of white-blooded icefishes, the only vertebrates living without hemoglobin. To understand hemoglobin evolution during adaptation to freezing water, we analyzed hemoglobin genes from 36 notothenioid genome assemblies. Results showed that adaptation to frigid conditions shaped hemoglobin gene evolution by episodic diversifying selection concomitant with cold adaptation and by pervasive evolution in Antarctic notothenioids compared to temperate relatives, likely a continuing adaptation to Antarctic conditions. Analysis of hemoglobin gene expression in adult hematopoietic organs in various temperate and Antarctic species further revealed a switch in hemoglobin gene expression underlying hemoglobin multiplicity reduction in Antarctic fish, leading to a single hemoglobin isoform in adult plunderfishes and dragonfishes, the sister groups to icefishes. The predicted high hemoglobin multiplicity in Antarctic fish embryos based on transcriptomic data, however, raises questions about the molecular bases and physiological implications of diverse hemoglobin isoforms in embryos compared to adults. This analysis supports the hypothesis that the last common icefish ancestor was vulnerable to detrimental mutations affecting the single ancestral expressed alpha- and beta-globin gene pair, potentially predisposing their subsequent loss.

A de novo missense mutation in MPP2 confers an increased risk of Vogt–Koyanagi–Harada disease as shown by trio-based whole-exome sequencing

Vogt–Koyanagi–Harada (VKH) disease is a leading cause of blindness in young and middle-aged people. However, the etiology of VKH disease remains unclear. Here, we performed the first trio-based whole-exome sequencing study, which enrolled 25 VKH patients and 50 controls, followed by a study of 2081 VKH patients from a Han Chinese population to uncover detrimental mutations. A total of 15 de novo mutations in VKH patients were identified, with one of the most important being the membrane palmitoylated protein 2 (MPP2) p.K315N (MPP2-N315) mutation. The MPP2-N315 mutation was highly deleterious according to bioinformatic predictions. Additionally, this mutation appears rare, being absent from the 1000 Genome Project and Genome Aggregation Database, and it is highly conserved in 10 species, including humans and mice. Subsequent studies showed that pathological phenotypes and retinal vascular leakage were aggravated in MPP2-N315 mutation knock-in or MPP2-N315 adeno-associated virus-treated mice with experimental autoimmune uveitis (EAU). In vitro, we used clustered regularly interspaced short palindromic repeats (CRISPR‒Cas9) gene editing technology to delete intrinsic MPP2 before overexpressing wild-type MPP2 or MPP2-N315. Levels of cytokines, such as IL-1β, IL-17E, and vascular endothelial growth factor A, were increased, and barrier function was destroyed in the MPP2-N315 mutant ARPE19 cells. Mechanistically, the MPP2-N315 mutation had a stronger ability to directly bind to ANXA2 than MPP2-K315, as shown by LC‒MS/MS and Co-IP, and resulted in activation of the ERK3/IL-17E pathway. Overall, our results demonstrated that the MPP2-K315N mutation may increase susceptibility to VKH disease.

Computational Identification of Most Deleterious Missense Mutations in Human PD-1 Gene.

Traditional cancer treatment approaches are often hindered by the presence of toxic side effects and the high rate of relapse observed in treated organs. In contrast, novel immunotherapeutic strategies targeting immune checkpoint inhibitors, particularly PD-1, have demonstrated promising results with minimal adverse effects. However, the emergence of immunotherapeutic-resistant tumors, predominantly caused by intrinsic mutations, poses a significant obstacle to successful treatment outcomes. Consequently, the primary objective of this study was to screen for the most detrimental missense mutations in the PD-1 gene associated with immunotherapeutic resistance. To achieve this aim, a comprehensive screening process utilizing 20 web servers, incorporating both sequence- and structure-based methodologies, was undertaken. Through meticulous analysis and mutual disease association sorting, four specific missense mutations were successfully identified. These mutations, namely, R38C, D61V, R94C, and D117V, emerged as the leading contributors to genetic cancer progression and immunotherapeutic resistance against PD-1 blockers. The findings presented in this study are supported by multiple lines of evidence. A thorough examination of protein topology, structural alignment, docking interactions with PD-L1, and protein flexibility collectively confirmed the pathogenic nature of these sorted mutations. By considering these various aspects, we have gained a comprehensive understanding of the underlying mechanisms driving immunotherapeutic resistance. In conclusion, the comprehensive screening process undertaken in this study has successfully identified R38C, D61V, R94C, and D117V as the primary mutations contributing to genetic cancer progression and immunotherapeutic resistance against PD-1 blockers. The integration of protein topology analysis, structural alignment, docking studies with PD-L1, and assessment of protein flexibility have collectively provided robust evidence to support the pathogenic significance of these mutations.

Sunflower Hybrids and Inbred Lines Adopt Different Physiological Strategies and Proteome Responses to Cope with Water Deficit.

Sunflower is a hybrid crop that is considered moderately drought-tolerant and adapted to new cropping systems required for the agro-ecological transition. Here, we studied the impact of hybridity status (hybrids vs. inbred lines) on the responses to drought at the molecular and eco-physiological level exploiting publicly available datasets. Eco-physiological traits and leaf proteomes were measured in eight inbred lines and their sixteen hybrids grown in the high-throughput phenotyping platform Phenotoul-Heliaphen. Hybrids and parental lines showed different growth strategies: hybrids grew faster in the absence of water constraint and arrested their growth more abruptly than inbred lines when subjected to water deficit. We identified 471 differentially accumulated proteins, of which 256 were regulated by drought. The amplitude of up- and downregulations was greater in hybrids than in inbred lines. Our results show that hybrids respond more strongly to water deficit at the molecular and eco-physiological levels. Because of presence/absence polymorphism, hybrids potentially contain more genes than their parental inbred lines. We propose that detrimental homozygous mutations and the lower number of genes in inbred lines lead to a constitutive defense mechanism that may explain the lower growth of inbred lines under well-watered conditions and their lower reactivity to water deficit.

Computational screening of pathogenic missense nsSNPs in heme oxygenase 1 (HMOX1) gene and their structural and functional consequences

Heme Oxygenase 1 (HMOX1) is a cytoprotective enzyme, exhibiting the highest activity in the spleen, catalyzing the heme ring breakdown into products of biological significance- biliverdin, CO, and Fe2+. In vascular cells, HMOX1 possesses strong anti-apoptotic, antioxidant, anti-proliferative, anti-inflammatory, and immunomodulatory actions. The majority of these activities are crucial for the prevention of atherogenesis. Single amino acid substitutions in proteins generated by missense non-synonymous single nucleotide polymorphism (nsSNPs) in the protein-encoding regions of genes are potent enough to cause significant medical challenges due to the alteration of protein structure and function. The current study aimed at characterizing and analyzing high-risk nsSNPs associated with the human HMOX1 gene. Preliminary screening of the total available 288 missense SNPs was performed through the lens of deleteriousness and stability prediction tools. Finally, a total of seven nsSNPs (Y58D, A131T, Y134H, F166S, F167S, R183S and M186V) were found to be most deleterious by all tools that are present at highly conserved positions. Molecular dynamics simulations (MDS) analysis explained the mutational effects on the dynamic action of the wild-type and mutant proteins. In a nutshell, R183S (rs749644285) was identified as a highly detrimental mutation that could significantly render the enzymatic activity of HMOX1. The finding of this computational analysis might help subject the experimental confirmatory analysis to characterize the role of nsSNPs in HMOX1. Communicated by Ramaswamy H. Sarma

Identifying the Molecular Drivers of Pathogenic Aldehyde Dehydrogenase Missense Mutations in Cancer and Non-Cancer Diseases.

Human aldehyde dehydrogenases (ALDHs) comprising 19 isoenzymes play a vital role on both endogenous and exogenous aldehyde metabolism. This NAD(P)-dependent catalytic process relies on the intact structural and functional activity of the cofactor binding, substrate interaction, and the oligomerization of ALDHs. Disruptions on the activity of ALDHs, however, could result in the accumulation of cytotoxic aldehydes, which have been linked with a wide range of diseases, including both cancers as well as neurological and developmental disorders. In our previous works, we have successfully characterised the structure-function relationships of the missense variants of other proteins. We, therefore, applied a similar analysis pipeline to identify potential molecular drivers of pathogenic ALDH missense mutations. Variants data were first carefully curated and labelled as cancer-risk, non-cancer diseases, and benign. We then leveraged various computational biophysical methods to describe the changes caused by missense mutations, informing a bias of detrimental mutations with destabilising effects. Cooperating with these insights, several machine learning approaches were further utilised to investigate the combination of features, revealing the necessity of the conservation of ALDHs. Our work aims to provide important biological perspectives on pathogenic consequences of missense mutations of ALDHs, which could be invaluable resources in the development of cancer treatment.

High BM plasma S100A8/A9 is associated with a perturbed microenvironmentand poor prognosis in myelodysplastic syndromes.

S100A8/A9 is a proinflammatory protein and plays an essential role in the pathogenesis of myelodysplastic syndromes (MDS) via the S100A8/A9-Toll-like receptors axis. While S100A8/A9 levels have been used as biomarkers in many inflammatory diseases, their clinical relevance has not been conclusively resolved in MDS. To address this, we used an enzyme-linked immunosorbent assay to quantify S100A8/A9 heterodimers in bone marrow (BM) plasma from 215 MDS patients and compared S100A8/A9 levels across patients with various disease risks and genotypes. S100A8/A9 levels correlated with ASXL1 variant allele frequencies significantly. Moreover, mutant ASXL1 with concurrent RUNX1, STAG2, ZRSR2, or EZH2 mutations was associated with higher S100A8/A9 levels. We further showed that higher S100A8/A9 independently predicted inferior leukemia-free survival and overall survival in MDS patients, irrespective of age, Revised International Prognostic Scoring System subgroups, and known detrimental mutations. Lastly, through deep-sequenced transcriptomic analysis, we demonstrated that higher S100A8/A9 in the BM intimated a perturbed microenvironment with enhanced myeloid-derived suppressor cell-mediated tumor immune escape signal, altered metabolism, and impairment in the functions and quantities of CD8+ T cells and NK cells. S100A8/A9 in the BM microenvironment may be a potential biomarker in the prognostication of MDS and target for novel therapy.

Enzymatic Cleavage-Mediated Extension Stalling Enables Accurate Recognition and Quantification of Locus-Specific Uracil Modification in DNA.

Chemical modifications in DNA have profound influences on the structures and functions of DNA. Uracil, a naturally occurring DNA modification, can originate from the deamination of cytosine or arise from misincorporation of dUTP into DNA during DNA replication. Uracil in DNA will imperil genomic stability due to their potential in producing detrimental mutations. An in-depth understanding of the functions of uracil modification requires the accurate determination of its site as well as content in genomes. Herein, we characterized that a new member of the uracil-DNA glycosylase (UDG) family enzyme (UdgX-H109S) could selectively cleave both uracil-containing single-stranded DNA (ssDNA) and double-stranded DNA (dsDNA). Based on this unique property of UdgX-H109S, we developed an enzymatic cleavage-mediated extension stalling (ECES) method for the locus-specific detection and quantification of uracil in genomic DNA. In the ECES method, UdgX-H109S specifically recognizes and cleaves the N-glycosidic bond of uracil from dsDNA and generates an apurinic/apyrimidinic (AP) site, which could be broken by APE1 to form a one-nucleotide gap. The specific cleavage by UdgX-H109S is then evaluated and quantified by qPCR. With the developed ECES approach, we demonstrated that the level of uracil at position Chr4:50566961 in genomic DNA of breast cancer tissues was significantly decreased. Collectively, the ECES method has been proved to be accurate and reproducible in the locus-specific quantification of uracil in genomic DNA from biological and clinical samples.

Predicting the effects of rare genetic variants on oncogenic signaling pathways: A computational analysis of HRAS protein function.

The HRAS gene plays a crucial role in regulating essential cellular processes for life, and this gene's misregulation is linked to the development of various types of cancers. Nonsynonymous single nucleotide polymorphisms (nsSNPs) within the coding region of HRAS can cause detrimental mutations that disrupt wild-type protein function. In the current investigation, we have employed in-silico methodologies to anticipate the consequences of infrequent genetic variations on the functional properties of the HRAS protein. We have discovered a total of 50 nsSNPs, of which 23 were located in the exon region of the HRAS gene and denoting that they were expected to cause harm or be deleterious. Out of these 23, 10 nsSNPs ([G60V], [G60D], [R123P], [D38H], [I46T], [G115R], [R123G], [P11OL], [A59L], and [G13R]) were identified as having the most delterious effect based on results of SIFT analysis and PolyPhen2 scores ranging from 0.53 to 69. The DDG values -3.21 kcal/mol to 0.87 kcal/mol represent the free energy change associated with protein stability upon mutation. Interestingly, we identified that the three mutations (Y4C, T58I, and Y12E) were found to improve the structural stability of the protein. We performed molecular dynamics (MD) simulations to investigate the structural and dynamic effects of HRAS mutations. Our results showed that the stable model of HRAS had a significantly lower energy value of -18756 kj/mol compared to the initial model of -108915 kj/mol. The RMSD value for the wild-type complex was 4.40 Å, and the binding energies for the G60V, G60D, and D38H mutants were -107.09 kcal/mol, -109.42 kcal/mol, and -107.18 kcal/mol, respectively as compared to wild-type HRAS protein had -105.85 kcal/mol. The result of our investigation presents convincing corroboration for the potential functional significance of nsSNPs in augmenting HRAS expression and adding to the activation of malignant oncogenic signalling pathways.

Single-cell individual full-length mtDNA sequencing by iMiGseq uncovers unexpected heteroplasmy shifts in mtDNA editing.

The ontogeny and dynamics of mtDNA heteroplasmy remain unclear due to limitations of current mtDNA sequencing methods. We developed individual Mitochondrial Genome sequencing (iMiGseq) of full-length mtDNA for ultra-sensitive variant detection, complete haplotyping, and unbiased evaluation of heteroplasmy levels, all at the individual mtDNA molecule level. iMiGseq uncovered unappreciated levels of heteroplasmic variants in single cells well below the conventional NGS detection limit and provided accurate quantitation of heteroplasmy level. iMiGseq resolved the complete haplotype of individual mtDNA in single oocytes and revealed genetic linkage of de novo mutations. iMiGseq detected sequential acquisition of detrimental mutations, including large deletions, in defective mtDNA in NARP/Leigh syndrome patient-derived induced pluripotent stem cells. iMiGseq identified unintended heteroplasmy shifts in mitoTALEN editing, while showing no appreciable level of unintended mutations in DdCBE-mediated mtDNA base editing. Therefore, iMiGseq could not only help elucidate the mitochondrial etiology of diseases, but also evaluate the safety of various mtDNA editing strategies.

Cryo-EM structures of human SPCA1a reveal the mechanism of Ca2+/Mn2+ transport into the Golgi apparatus.

Secretory pathway Ca2+/Mn2+ ATPase 1 (SPCA1) actively transports cytosolic Ca2+ and Mn2+ into the Golgi lumen, playing a crucial role in cellular calcium and manganese homeostasis. Detrimental mutations of the ATP2C1 gene encoding SPCA1 cause Hailey-Hailey disease. Here, using nanobody/megabody technologies, we determined cryo-electron microscopy structures of human SPCA1a in the ATP and Ca2+/Mn2+-bound (E1-ATP) state and the metal-free phosphorylated (E2P) state at 3.1- to 3.3-Å resolutions. The structures revealed that Ca2+ and Mn2+ share the same metal ion-binding pocket with similar but notably different coordination geometries in the transmembrane domain, corresponding to the second Ca2+-binding site in sarco/endoplasmic reticulum Ca2+-ATPase (SERCA). In the E1-ATP to E2P transition, SPCA1a undergoes similar domain rearrangements to those of SERCA. Meanwhile, SPCA1a shows larger conformational and positional flexibility of the second and sixth transmembrane helices, possibly explaining its wider metal ion specificity. These structural findings illuminate the unique mechanisms of SPCA1a-mediated Ca2+/Mn2+ transport.

Cryo-EM structures of the TolQ membrane protein in complex with inner membrane protein partner TolR from the opportunistic pathogen Acinetobacter baumannii provide information on structure and stoichiometry.

TolQ-TolR-TolA is an energy-transducing inner-membrane sub-complex conserved across Gram-negative bacteria. It maintains connection with periplasmic TolB-Pal sub-complex in an energy-dependent manner through TolA protein spanning the periplasm. Utilizing proton motive force (PMF), Tol-Pal complex contributes to the integrity of outer membrane and the whole cell envelope. TolQ remains the only part of Tol-Pal complex with unknown structure up to now. To address this, we determined 3.02 and 3.34 Å cryo-electron microscopy (cryo-EM) structures of the integral inner-membrane protein TolQ and in complex with TolR, respectively from Acinetobacter baumannii. A. baumannii is one of the most dreaded and clinically important Gram-negative pathogens causing a broad-spectrum nosocomial infection mainly among debilitated patients admitted to intensive care units. Our structural analysis revealed that TolQ forms a pentamer with a channel-like hydrophobic trans-membrane pore slightly shifted above the membrane opening into a large negatively charged cytoplasmic chamber which is in common with the closest TolQ-TolR paralogs ExbB-ExbD and MotA-MotB. We validated observed structural features of TolQ by prior functional data. According to our analysis the most functionally detrimental point mutations mapped on structure are localized inside the channel. Obtained structures suggest that two trans-membrane domains of TolR dimer fit inside TolQ channel in TolQ-TolR complex forming a proton-flow pathway for PMF utilization. Cryo-EM structures solved in this work were funded in whole or in part with Refeyn Travel Grant and U.S. federal funds from the National Institute of Allergy and Infectious Diseases, National Institutes of Health, Department of Health and Human Services, under contract No. HHSN272201700060C (Center for Structural Genomics of Infectious Diseases (CSGID); http://csgid.org).