Sequence Mutagenesis Research Articles

Regression convolutional neural network models implicate peripheral immune regulatory variants in the predisposition to Alzheimer's disease.

Alzheimer's disease (AD) involves aggregation of amyloid β and tau, neuron loss, cognitive decline, and neuroinflammatory responses. Both resident microglia and peripheral immune cells have been associated with the immune component of AD. However, the relative contribution of resident and peripheral immune cell types to AD predisposition has not been thoroughly explored due to their similarity in gene expression and function. To study the effects of AD-associated variants on cis-regulatory elements, we train convolutional neural network (CNN) regression models that link genome sequence to cell type-specific levels of open chromatin, a proxy for regulatory element activity. We then use in silico mutagenesis of regulatory sequences to predict the relative impact of candidate variants across these cell types. We develop and apply criteria for evaluating our models and refine our models using massively parallel reporter assay (MPRA) data. Our models identify multiple AD-associated variants with a greater predicted impact in peripheral cells relative to microglia or neurons. Our results support their use as models to study the effects of AD-associated variants and even suggest that peripheral immune cells themselves may mediate a component of AD predisposition. We make our library of CNN models and predictions available as a resource for the community to study immune and neurological disorders.

Protein SUMOylation promotes cAMP-independent EPAC1 activation

Protein SUMOylation is a prevalent stress-response posttranslational modification crucial for maintaining cellular homeostasis. Herein, we report that protein SUMOylation modulates cellular signaling mediated by cAMP, an ancient and universal stress-response second messenger. We identify K561 as a primary SUMOylation site in exchange protein directly activated by cAMP (EPAC1) via site-specific mapping of SUMOylation using mass spectrometry. Sequence and site-directed mutagenesis analyses reveal that a functional SUMO-interacting motif in EPAC1 is required for the binding of SUMO-conjugating enzyme UBC9, formation of EPAC1 nuclear condensate, and EPAC1 cellular SUMOylation. Heat shock-induced SUMO modification of EPAC1 promotes Rap1/2 activation in a cAMP-independent manner. Structural modeling and molecular dynamics simulation studies demonstrate that SUMO substituent on K561 of EPAC1 promotes Rap1 interaction by increasing the buried surface area between the SUMOylated receptor and its effector. Our studies identify a functional SUMOylation site in EPAC1 and unveil a novel mechanism in which SUMOylation of EPAC1 leads to its autonomous activation. The findings of SUMOylation-mediated activation of EPAC1 not only provide new insights into our understanding of cellular regulation of EPAC1 but also will open up a new field of experimentation concerning the cross-talk between cAMP/EPAC1 signaling and protein SUMOylation, two major cellular stress response pathways, during cellular homeostasis.

Multiplexed CRISPR-Cas9 mutagenesis of rice PSBS1 noncoding sequences for transgene-free overexpression.

Understanding CRISPR-Cas9's capacity to produce native overexpression (OX) alleles would accelerate agronomic gains achievable by gene editing. To generate OX alleles with increased RNA and protein abundance, we leveraged multiplexed CRISPR-Cas9 mutagenesis of noncoding sequences upstream of the rice PSBS1 gene. We isolated 120 gene-edited alleles with varying non-photochemical quenching (NPQ) capacity in vivo-from knockout to overexpression-using a high-throughput screening pipeline. Overexpression increased OsPsbS1 protein abundance two- to threefold, matching fold changes obtained by transgenesis. Increased PsbS protein abundance enhanced NPQ capacity and water-use efficiency. Across our resolved genetic variation, we identify the role of 5'UTR indels and inversions in driving knockout/knockdown and overexpression phenotypes, respectively. Complex structural variants, such as the 252-kb duplication/inversion generated here, evidence the potential of CRISPR-Cas9 to facilitate significant genomic changes with negligible off-target transcriptomic perturbations. Our results may inform future gene-editing strategies for hypermorphic alleles and have advanced the pursuit of gene-edited, non-transgenic rice plants with accelerated relaxation of photoprotection.

Genetic determinants of host- and virus-derived insertions for hepatitis E virus replication

Hepatitis E virus (HEV) is a long-neglected RNA virus and the major causative agent of acute viral hepatitis in humans. Recent data suggest that HEV has a very heterogeneous hypervariable region (HVR), which can tolerate major genomic rearrangements. In this study, we identify insertions of previously undescribed sequence snippets in serum samples of a ribavirin treatment failure patient. These insertions increase viral replication while not affecting sensitivity towards ribavirin in a subgenomic replicon assay. All insertions contain a predicted nuclear localization sequence and alanine scanning mutagenesis of lysine residues in the HVR influences viral replication. Sequential replacement of lysine residues additionally alters intracellular localization in a fluorescence dye-coupled construct. Furthermore, distinct sequence patterns outside the HVR are identified as viral determinants that recapitulate the enhancing effect. In conclusion, patient-derived insertions can increase HEV replication and synergistically acting viral determinants in and outside the HVR are described. These results will help to understand the underlying principles of viral adaptation by viral- and host-sequence snatching during the clinical course of infection.

Constitutive activity of ionotropic glutamate receptors via hydrophobic substitutions in the ligand-binding domain

Neurotransmitter ligands electrically excite neurons by activating ionotropic glutamate receptor (iGluR) ion channels. Knowledge of the iGluR amino acid residues that dominate ligand-induced activation would enable the prediction of function from sequence. We therefore explored the molecular determinants of activity in rat N-methyl-D-aspartate (NMDA)-type iGluRs (NMDA receptors), complex heteromeric iGluRs comprising two glycine-binding GluN1 and two glutamate-binding GluN2 subunits, using amino acid sequence analysis, mutagenesis, and electrophysiology. We find that a broadly conserved aspartate residue controls both ligand potency and channel activity, to the extent that certain substitutions at this position bypass the need for ligand binding in GluN1 subunits, generating NMDA receptors activated solely by glutamate. Furthermore, we identify a homomeric iGluR from the placozoan Trichoplax adhaerens that has utilized native mutations of this crucial residue to evolve into a leak channel that is inhibited by neurotransmitter binding, pointing to a dominant role of this residue throughout the iGluR superfamily.

Pycofitness-Evaluating the fitness landscape of RNA and protein sequences.

The accurate prediction of how mutations change biophysical properties of proteins or RNA is a major goal in computational biology with tremendous impacts on protein design and genetic variant interpretation. Evolutionary approaches such as coevolution can help solving this issue. We present pycofitness, a standalone Python-based software package for the in silico mutagenesis of protein and RNA sequences. It is based on coevolution and, more specifically, on a popular inverse statistical approach, namely direct coupling analysis by pseudo-likelihood maximization. Its efficient implementation and user-friendly command line interface make it an easy-to-use tool even for researchers with no bioinformatics background. To illustrate its strengths, we present three applications in which pycofitness efficiently predicts the deleteriousness of genetic variants and the effect of mutations on protein fitness and thermodynamic stability. https://github.com/KIT-MBS/pycofitness.

Fine-Tuning Gene Expression in Bacteria by Synthetic Promoters.

Promoters are key genetic elements in the initiation and regulation of gene expression. A limited number of natural promoters has been described for the control of gene expression in synthetic biology applications. Therefore, synthetic promoters have been developed to fine-tune the transcription for the desired amount of gene product. Mostly, synthetic promoters are characterized using promoter libraries that are constructed via mutagenesis of promoter sequences. The strength of promoters in the library is determined according to the expression of a reporter gene such as gfp encoding green fluorescent protein. Gene expression can be controlled using inducers. The majority of the studies on gram-negative bacteria are conducted using the expression system of the model organism Escherichia coli while that of the model organism Bacillus subtilis is mostly used in the studies on gram-positive bacteria. Additionally, synthetic promoters for the cyanobacteria, which are phototrophic microorganisms, are evaluated, especially using the model cyanobacterium Synechocystis sp. PCC 6803. Moreover, a variety of algorithms based on machine learning methods were developed to characterize the features of promoter elements. Some of these in silico models were verified using in vitro or in vivo experiments. Identification of novel synthetic promoters with improved features compared to natural ones contributes much to the synthetic biology approaches in terms of fine-tuning gene expression.

METTL17 is an Fe-S cluster checkpoint for mitochondrial translation

Friedreich's ataxia (FA) is a debilitating, multisystemic disease caused by the depletion of frataxin (FXN), a mitochondrial iron-sulfur (Fe-S) cluster biogenesis factor. To understand the cellular pathogenesis of FA, we performed quantitative proteomics in FXN-deficient human cells. Nearly every annotated Fe-S cluster-containing protein was depleted, indicating that as a rule, cluster binding confers stability to Fe-S proteins. We also observed depletion of a small mitoribosomal assembly factor METTL17 and evidence of impaired mitochondrial translation. Using comparative sequence analysis, mutagenesis, biochemistry, and cryoelectron microscopy, we show that METTL17 binds to the mitoribosomal small subunit during late assembly and harbors a previously unrecognized [Fe4S4]2+ cluster required for its stability. METTL17 overexpression rescued the mitochondrial translation and bioenergetic defects, but not the cellular growth, of FXN-depleted cells. These findings suggest that METTL17 acts as an Fe-S cluster checkpoint, promoting translation of Fe-S cluster-rich oxidative phosphorylation (OXPHOS) proteins only when Fe-S cofactors are replete.

Functional characterization of a TerC family protein of Riemerella anatipestifer in manganese detoxification and virulence.

Manganese (Mn) is an essential element for bacteria, but the overload of manganese is toxic. In a previous study, we showed that the cation diffusion facilitator protein MetA and the resistance-nodulation-division efflux pump MetB are responsible for Mn efflux in the bacterial pathogen Riemerella anatipestifer CH-1. However, whether this bacterium encodes additional manganese efflux proteins is unclear. In this study, we show that R. anatipestifer CH-1 encodes a tellurium resistance C (TerC) family protein with low similarity to other characterized TerC family proteins. Compared to the wild type (WT), the terC mutant of R. anatipestifer CH-1 (∆terC) is sensitive to Mn(II) intoxication. The ability of TerC to export manganese is higher than that of MetB but lower than that of MetA. Consistently, terC deletion (∆terC) led to intracellular accumulation of Mn2+ under excess manganese conditions. Further study showed that ∆terC was more sensitive than the WT to the oxidant hypoclorite but not to hydrogen peroxide. Mutagenesis studies showed that the mutant at amino acid sites of Glu116 (E116), Asp122 (D122), Glu245 (E245) Asp248 (D248), and Asp254 (D254) may be involved in the ability of TerC to export manganese. The transcription of terC was upregulated under excess manganese and downregulated under iron-limited conditions. However, this was not dependent on the manganese metabolism regulator MetR. In contrast to a strain lacking the manganese efflux pump MetA or MetB, the terC mutant is attenuated in virulence in a duckling model of infection due to increased sensitivity to duck serum. Finally, comparative analysis showed that homologs of TerC are distributed across the bacterial kingdom, suggesting that TerC exerts a conserved manganese efflux function.IMPORTANCERiemerella anatipestifer is a notorious bacterial pathogen of ducks and other birds. In R. anatipestifer, the genes involved in manganese efflux have not been completely identified, although MetA and MetB have been identified as two manganese exporters. Additionally, the function of TerC family proteins in manganese efflux is controversial. Here, we demonstrated that a TerC family protein helps prevent Mn(II) intoxication in R. anatipestifer and that the ability of TerC to export manganese is intermediate compared to that of MetA and MetB. Sequence analysis and mutagenesis studies showed that the conserved key amino sites of TerC are Glu116, Asp122, Glu245, Asp248, and Asp254. The transcription of terC was regulated by manganese excess and iron limitation. Finally, we show that TerC plays a role in the virulence of R. anatipestifer due to the increased sensitivity to duck serum, rather than the increased sensitivity to manganese. Taken together, these results expand our understanding of manganese efflux and the pathogenic mechanisms of R. anatipestifer.

The Arabidopsis SWEET1 and SWEET2 uniporters recognize similar substrates while differing in subcellular localization

Sugars Will Eventually be Exported Transporters (SWEETs) are central for sugar allocation in plants. The SWEET family has approximately 20 homologs in most plant genomes, and despite extensive research on their structures and molecular functions, it is still unclear how diverse SWEETs recognize different substrates. Previous work using SweetTrac1, a biosensor constructed by the intramolecular fusion of a conformation-sensitive fluorescent protein in the plasma membrane transporter SWEET1 from Arabidopsis thaliana, identified common features in the transporter's substrates. Here, we report SweetTrac2, a new biosensor based on the Arabidopsis vacuole membrane transporter SWEET2, and use it to explore the substrate specificity of this second protein. Our results show that SWEET1 and SWEET2 recognize similar substrates but some with different affinities. Sequence comparison and mutagenesis analysis support the conclusion that the differences in affinity depend on nonspecific interactions involving previously uncharacterized residues in the substrate-binding pocket. Furthermore, SweetTrac2 can be an effective tool for monitoring sugar transport at vacuolar membranes that would be otherwise challenging to study.

Functional and structural characterization of Streptococcus pneumoniae pyruvate kinase involved in fosfomycin resistance

Glycolysis is the primary metabolic pathway in the strictly fermentative Streptococcus pneumoniae, which is a major human pathogen associated with antibiotic resistance. Pyruvate kinase (PYK) is the last enzyme in this pathway that catalyzes the production of pyruvate from phosphoenolpyruvate (PEP) and plays a crucial role in controlling carbon flux; however, while S.pneumoniae PYK (SpPYK) is indispensable for growth, surprisingly little is known about its functional properties. Here, we report that compromising mutations in SpPYK confers resistance to the antibiotic fosfomycin, which inhibits the peptidoglycan synthesis enzyme MurA, implying a direct link between PYK and cell wall biogenesis. The crystal structures of SpPYK in the apo and ligand-bound states reveal key interactions that contribute to its conformational change as well as residues responsible for the recognition of PEP and the allosteric activator fructose 1,6-bisphosphate (FBP). Strikingly, FBP binding was observed at a location distinct from previously reported PYK effector binding sites. Furthermore, we show that SpPYK could be engineered to become more responsive to glucose 6-phosphate instead of FBP by sequence and structure-guided mutagenesis of the effector binding site. Together, our work sheds light on the regulatory mechanism of SpPYK and lays the groundwork for antibiotic development that targets this essential enzyme.

Functional differences between voltage-gated proton channels from flowering and non-flowering plants.

pH regulation is essential for plant growth and development, as well as for the response to abiotic stress and immunity. Some of the constituents of the pH regulation network have been identified, but little is known about the contribution of proton channels to plant physiology. The existence of plant Hv proteins was first proposed more than a decade ago based on sequence homology with proton channels identified in algae, but whether these proteins can function as proton channels has not been demonstrated. We have identified Hv proton channels from flowering plants that gate with a unique modality, as they can be turned on by an electrical stimulus if they have been first exposed to a mechanical stimulus in a process that we call priming. Homolog channels from non-flowering plants do not share this characteristic, as they can be turned on by the electrical stimulus alone like their counterparts in the animal kingdom. Using sequence homology analysis, mutagenesis, and electrophysiological measurements, combined with AI-based structural modeling, we investigated the mechanism of the new gating modality. We found several molecular features of the S4 transmembrane segment of the channel that are responsible for mechanical priming. The activity of human Hv1 is known to be potentiated by membrane stretch through a mechanism poorly understood. Here, we propose a simple relationship between priming of Hvs from flowering plants and stretch-mediated potentiation of human Hv1.

Random Mutagenesis by PCR.

Random mutagenesis of DNA sequences has the advantage of generating DNA sequences with novel properties, either directly in the case of aptamers or through subsequent transcription/translation of the mutated sequence in the case of proteins. In both cases no prior structural or mechanistic knowledge of the molecule is required. For sequences greater than 100bp, one of the easiest methods to introduce the mutations is to use Error-prone PCR (EP-PCR) as discussed in this chapter. When coupled with an appropriate selection or high throughput screening methodology, PCR-based random mutagenesis can provide a powerful tool for modern molecular biologists.

Nanomolar biosensor for detection of phenylacetic acid and L-phenylalanine

Phenylacetic acid (PAA) plays an important role in the plant hormone biology, possess a substantial antimicrobial activity and widely used building block for several pharmaceutical drugs. Technological advancement for PAA and other small molecule investigation in vivo can help to yield interesting insights into their involvement in biological processes. Transcription factor-based whole-cell biosensors have recently become useful analytical tools for detection and monitoring of various molecules. Compared to conventional analytical methods, the whole-cell biosensor offers cost-effective, rapid and high-throughput solution. In this paper, we identify and characterize PAA-inducible gene expression system CnPaaX/PpaaA2 from Cupriavidus necator H16. Nucleotide sequence analysis and mutagenesis reveal a conservative region with palindromic sequence motifs upstream to paaA2, which are likely used for PaaX-mediated transcriptional regulation. Using inducible system CnPaaX/PpaaA2 we develop C. necator-based whole-cell biosensor that exhibits a dose-dependent response to the exogenous PAA. Moreover, the biosensor is highly specific to the PAA compared to naturally occurring hydroxylated forms of this compound. Importantly, it exhibits a limit of measurement as low as 9.54 nM demonstrating the highest level of sensitivity reported so far. By combining genetic elements of PAA-inducible system with phenylpyruvate decarboxylase gene ipdC, we develop biosensor suitable for determining of L-phenylalanine.

CRISPR base editing of cis-regulatory elements enables the perturbation of neurodegeneration-linked genes

CRISPR technology has demonstrated broad utility for controlling target gene expression; however, there remains a need for strategies capable of modulating expression via the precise editing of non-coding regulatory elements. Here, we demonstrate that CRISPR base editors, a class of gene-modifying proteins capable of creating single-base substitutions in DNA, can be used to perturb gene expression via their targeted mutagenesis of cis-acting sequences. Using the promoter region of the human huntingtin (HTT) gene as an initial target, we show that editing of the binding site for the transcription factor NF-κB led to a marked reduction in HTT gene expression in base-edited cell populations. We found that these gene perturbations were persistent and specific, as a transcriptome-wide RNA analysis revealed minimal off-target effects resulting from the action of the base editor protein. We further demonstrate that this base-editing platform could influence gene expression invivo as its delivery to a mouse model of Huntington's disease led to a potent decrease in HTT mRNA in striatal neurons. Finally, to illustrate the applicability of this concept, we target the amyloid precursor protein, showing that multiplex editing of its promoter region significantly perturbed its expression. These findings demonstrate the potential for base editors to regulate target gene expression.

Overcoming the limitations of cytokines to improve cancer therapy.

Cytokines are pleiotropic soluble proteins used by immune cells to orchestrate a coordinated response against pathogens and malignancies. In cancer immunotherapy, cytokine-based drugs can be developed potentiating pro-inflammatory cytokines or blocking immunosuppressive cytokines. However, the complexity of the mechanisms of action of cytokines requires the use of biotechnological strategies to minimize systemic toxicity, while potentiating the antitumor response. Sequence mutagenesis, fusion proteins and gene therapy strategies are employed to enhance the half-life in circulation, target the desired bioactivity to the tumor microenvironment, and to optimize the therapeutic window of cytokines. In this review, we provide an overview of the different strategies currently being pursued in pre-clinical and clinical studies to make the most of cytokines for cancer immunotherapy.

Identification of a Catalytic Active but Non-Aggregating MDM2 RING Domain Variant

As a key regulator of the tumour suppressor protein p53, MDM2 is involved in various types of cancer and has thus been an attractive drug target. So far, small molecule design has primarily focussed on the N-terminal p53-binding domain although on-target toxicity effects have been reported. Targeting the catalytic RING domain of MDM2 resembles an alternative approach to drug MDM2 with the idea to prevent MDM2-mediated ubiquitination of p53 while retaining MDM2′s ability to bind p53. The design of RING inhibitors has been limited by the extensive aggregation tendency of the RING domain, making it challenging to undertake co-crystallization attempts with potential inhibitors. Here we compare the purification profiles of the MDM2 RING domain from several species and show that the MDM2 RING domain of other species than human is much less prone to aggregate although the overall structure of the RING domain is conserved. Through sequence comparison and mutagenesis analyses, we identify a single point mutation, G443T, which greatly enhances the dimeric fraction of human MDM2 RING domain during purification. Neither does the mutation alter the structure of the RING domain, nor does it affect E2(UbcH5B)–Ub binding and activity. Hence, MDM2-G443T facilitates studies involving binding partners that would be hampered by the low solubility of the wild-type RING domain. Furthermore, it will be valuable for the development of MDM2 RING inhibitors.

From the notebook to recombinant protein production in Escherichia coli: Design of expression vectors and gene cloning.

Research in recombinant protein expression in microorganism hosts spans half a century. The field has evolved from mostly trial-and-error approaches to more rational strategies, including careful design of the expression vectors and the coding sequence for the protein of interest. It is important to reflect on many aspects about vector construction, such as codon usage, integration site, coding sequence mutagenesis and many others. In this chapter, we overview methods and considerations to generate a suitable construct and anticipate possible experimental roadblocks.

Genome-wide functional analysis of hot pepper immune receptors reveals an autonomous NLR clade in seed plants.

Summary Plants possess hundreds of intracellular immune receptors encoding nucleotide‐binding domain leucine‐rich repeat (NLR) proteins. Full‐length NLRs or a specific domain of NLRs often induce plant cell death in the absence of pathogen infection.In this study we used genome‐wide transient expression analysis to identify a group of NLRs (ANLs; ancient and autonomous NLRs) carrying autoactive coiled‐coil (CCA) domains in pepper (Capsicum annuum). CCA‐mediated cell death mimics hypersensitive cell death triggered by the interaction between NLRs and pathogen effectors.Sequence alignment and mutagenesis analyses revealed that the intact α1 helix of CCAs is critical for both CCA‐ and ANL‐mediated cell death. Cell death induced by CCAs does not require NRG1/ADR1 or NRC type helper NLRs, suggesting ANLs may function as singleton NLRs. We also found that CCAs localize to the plasma membrane, as demonstrated for Arabidopsis singleton NLR ZAR1. Extended studies revealed that autoactive CCAs are well conserved in other Solanaceae plants as well as in rice, a monocot plant. Further phylogenetic analyses revealed that ANLs are present in all tested seed plants (spermatophytes).Our study not only uncovers the autonomous NLR clade in plants but also provides powerful resources for dissecting the underlying molecular mechanism of NLR‐mediated cell death in plants.

Genome of Tetraploid Fish Schizothorax o'connori Provides Insights into Early Re-diploidization and High-Altitude Adaptation.

SummaryWhole-genome duplications (WGDs) of Schizothoracinae are believed to have played a significant role in speciation and environmental adaptation on the Qinghai-Tibet Plateau (QTP). Here, we present a genome for Schizothorax o'connori, a QTP endemic fish and showed the species as a young tetraploid with a recent WGD later than ∼1.23 mya. We exhibited that massive insertions between duplicated genomes caused by transposon bursts could induce mutagenesis in adjacent sequences and alter the expression of neighboring genes, representing an early re-diploidization process in a polyploid genome after WGD. Meanwhile, we found that many genes involved in DNA repair and folate transport/metabolism experienced natural selection and might contribute to the environmental adaptation of this species. Therefore, the S. o'connori genome could serve as a young tetraploid model for investigations of early re-diploidization in polyploid genomes and offers an invaluable genetic resource for environmental adaptation studies of the endemic fish of the QTP.