Single Nucleotide Mutations Research Articles

PatU3 plays a central role in coordinating cell division and differentiation in pattern formation of filamentous cyanobacterium Nostoc sp. PCC 7120.

Spatial periodic signal for cell differentiation in some multicellular organisms is generated according to Turing's principle for pattern formation. How a dividing cell responds to the signal of differentiation is addressed with the filamentous cyanobacterium Nostoc sp. PCC 7120, which forms the patterned distribution of heterocysts. We show that differentiation of a dividing cell was delayed until its division was completed and only one daughter cell became heterocyst. A mutant of patU3, which encodes an inhibitor of heterocyst formation, showed no such delay and formed heterocyst pairs from the daughter cells of cell division or dumbbell-shaped heterocysts from the cells undergoing cytokinesis. The patA mutant, which forms heterocysts only at the filament ends, restored intercalary heterocysts by a single nucleotide mutation of patU3, and double mutants of patU3/patA and patU3/hetF had the phenotypes of the patU3 mutant. We provide evidence that HetF, which can degrade PatU3, is recruited to cell divisome through its C-terminal domain. A HetF mutant with its N-terminal peptidase domain but lacking the C-terminal domain could not prevent the formation of heterocyst pairs, suggesting that the divisome recruitment of HetF is needed to sequester HetF for the delay of differentiation in dividing cells. Our study demonstrates that PatU3 plays a key role in cell-division coupled control of differentiation.

A novel missense compound heterozygous variant in TLR1 gene is associated with susceptibility to rheumatoid arthritis - structural perspective and functional annotations.

Besides human leukocyte antigen (HLA-DRB1) locus, more than 100 loci across the genome have been identified and linked with the onset, expression and/or progression of rheumatoid arthritis (RA). However, there are still grey areas in our understanding of the key genetic contributors of the disease, particularly in familial cases. In the present study, we have performed the whole exome sequencing (WES) of RA patients from two consanguineous families of Pakistan in a quest to identify novel, high-impact, RA-susceptibility genetic variants. Through stepwise filtering, around 17,000 variants (common in the affected members) were recognized, out of which 2651 were predicted to be deleterious. Of these, 196 had direct relevance to RA. When selected for homozygous recessive mode of inheritance, two novel pathogenic variants (c.1324T>C, p.Cys442→Arg442; c.2036T>C, p.Ile679→Thr679) in the TLR1 gene displayed the role of compound heterozygosity in modulating the phenotypic expression and penetrance of RA. The structural and functional consequences of the TLR1 missense single nucleotide mutations (Cys442→Arg442; Ile679→Thr679) were evaluated through molecular dynamic simulation (MDS) studies. Analysis showed domain's rigidification, conferring stability to mutant TLR1-TIR/TIRAP-TIR complex with concomitant increase in molecular interactions with pro-inflammatory cytokines, compared to the wild-type conformation. Gene co-expression network analysis highlighted interlinked partnering genes along with interleukin-6 production of TLR1 (corrected p-value 2.98e-4) and acetylcholine receptor activity of CHRNG (corrected p-value 6.12e-2) as highly enriched associated functions. The results, validated through case-control study subjects, suggested that the variants identified through WES and confirmed through Sanger sequencing and MDS are the novel disease variants and are likely to confer RA-susceptibility, independently and/or in a family-specific context. Key Points • Exploration of population based/ethno-specific big data is imperative to identify novel causal variants of RA. • Two new deleterious missense mutations in mutational hotspot exon 4 of TLR1 gene have been identified in Pakistani RA patients. • MD simulation data provides evidence for domain's rigidification, conferring stability to mutant TLR1-TIR/TIRAP-TIR complex, with concomitant increase in production of pro-inflammatory cytokines, thus adding to the onset/erosive outcome of RA.

A single amino acid substitution in MdLAZY1A dominantly impairs shoot gravitropism in Malus.

Plant architecture is 1 of the most important factors that determines crop yield potential and productivity. In apple (Malus domestica), genetic improvement of tree architecture has been challenging due to a long juvenile phase and growth as complex trees composed of a distinct scion and a rootstock. To better understand the genetic control of apple tree architecture, the dominant weeping growth phenotype was investigated. We report the identification of MdLAZY1A (MD13G1122400) as the genetic determinant underpinning the Weeping (W) locus that largely controls weeping growth in Malus. MdLAZY1A is 1 of the 4 paralogs in apple that are most closely related to AtLAZY1 involved in gravitropism in Arabidopsis (Arabidopsis thaliana). The weeping allele (MdLAZY1A-W) contains a single nucleotide mutation c.584T>C that leads to a leucine to proline (L195P) substitution within a predicted transmembrane domain that colocalizes with Region III, 1 of the 5 conserved regions in LAZY1-like proteins. Subcellular localization revealed that MdLAZY1A localizes to the plasma membrane and nucleus in plant cells. Overexpressing the weeping allele in apple cultivar Royal Gala (RG) with standard growth habit impaired its gravitropic response and altered the growth to weeping-like. Suppressing the standard allele (MdLAZY1A-S) by RNA interference (RNAi) in RG similarly changed the branch growth direction to downward. Overall, the L195P mutation in MdLAZY1A is genetically causal for weeping growth, underscoring not only the crucial roles of residue L195 and Region III in MdLAZY1A-mediated gravitropic response but also a potential DNA base editing target for tree architecture improvement in Malus and other crops.

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.

Competitive aptamer switch for modulating ligand binding affinity

Aptamers are short, single-stranded DNA or RNA molecules that selectively bind to a target molecule. Aptamer-complement duplex (ACD) is often used to design molecular switches capable of producing a detectable signal or triggering a structural change upon aptamer binding to a target. However, such aptamer switch generally faces an increased dissociation constant (Kd) due to the energy barrier of the complementary duplex. We reported a competitive hybridization mechanism to modulate the binding affinity of an ACD to a target adenosine. Using the computation-guided design, we calculated the aptamer folding energy for the duplex length from 11-nt to 15-nt, and experimentally measured increased apparent Kd values resulted from these extended duplexes. Using a set of strands to compete with the ACD hybridization, it reduced the aptamer folding energy to facilitate aptamer switches with decreased apparent Kd values ranging from over 400 μM without a competing strand to ∼30 μM with a competing strand. This competitive aptamer switch was also found sensitive to single-nucleotide mutations of a competing strand. Our work provides an approach to modulate the binding affinity and the sensitivity of aptamer-complement duplexes, which could be useful in the nucleic acids-based sensing and nanomedicine.

Experimental estimates of germline mutation rate in eukaryotes: a phylogenetic meta-analysis.

Mutation is the ultimate source of all genetic variation, and over the last 10 years the ready availability of whole-genome sequencing has permitted direct estimation of mutation rate for many non-model species across the tree of life. In this meta-analysis, we make a comprehensive search of the literature for mutation rate estimates in eukaryotes, identifying 140 mutation accumulation (MA) and parent-offspring (PO) sequencing studies covering 134 species. Based on these data, we revisit differences in the single-nucleotide mutation (SNM) rate between different phylogenetic lineages and update the known relationships between mutation rate and generation time, genome size, and nucleotide diversity-while accounting for phylogenetic nonindependence. We do not find a significant difference between MA and PO in estimated mutation rates, but we confirm that mammal and plant lineages have higher mutation rates than arthropods and that unicellular eukaryotes have the lowest mutation rates. We find that mutation rates are higher in species with longer generation times and larger genome sizes, even when accounting for phylogenetic relationships. Moreover, although nucleotide diversity is positively correlated with mutation rate, the gradient of the relationship is significantly less than one (on a logarithmic scale), consistent with higher mutation rates in populations with smaller effective size. For the 29 species for which data are available, we find that indel mutation rates are positively correlated with nucleotide mutation rates and that short deletions are generally more common than short insertions. Nevertheless, despite recent progress, no estimates of either SNM or indel mutation rates are available for the majority of deeply branching eukaryotic lineages-or even for most animal phyla. Even among charismatic megafauna, experimental mutation rate estimates remain unknown for amphibia and scarce for reptiles and fish.

SDH mutations, as potential predictor of chemotherapy prognosis in small cell lung cancer patients

PurposeSmall cell lung cancer (SCLC) is an aggressive and rapidly progressive malignant tumor characterized by a poor prognosis. Chemotherapy remains the primary treatment in clinical practice; however, reliable biomarkers for predicting chemotherapy outcomes are scarce.MethodsIn this study, 78 SCLC patients were stratified into “good” or “poor” prognosis cohorts based on their overall survival (OS) following surgery and chemotherapeutic treatment. Next-generation sequencing was employed to analyze the mutation status of 315 tumorigenesis-associated genes in tumor tissues obtained from the patients. The random forest (RF) method, validated by the support vector machine (SVM), was utilized to identify single nucleotide mutations (SNVs) with predictive power. To verify the prognosis effect of SNVs, samples from the cbioportal database were utilized.ResultsThe SVM and RF methods confirmed that 20 genes positively contributed to prognosis prediction, displaying an area under the validation curve with a value of 0.89. In the corresponding OS analysis, all patients with SDH, STAT3 and PDCD1LG2 mutations were in the poor prognosis cohort (15/15, 100%). Analysis of public databases further confirms that SDH mutations are significantly associated with worse OS.ConclusionOur results provide a potential stratification of chemotherapy prognosis in SCLC patients, and have certain guiding significance for subsequent precise targeted therapy.

Mimicking the last step of gene elongation: hints from the bacterial hisF gene

Gene elongation consists in an in-tandem duplication of a gene and divergence and fusion of the two copies, resulting in a gene constituted by two divergent paralogous modules. Many present-day proteins show internal repeats of amino acid sequences, generated by gene elongation events; however, gene elongation is still a poorly studied evolutionary molecular mechanism. The most documented case is that of the histidine biosynthetic genes hisA and hisF, which derive from the gene elongation of an ancestral gene half the size of the extant ones.The aim of this work was to experimentally simulate the possible last step of the gene elongation event occurred during hisF gene evolution under selective pressure conditions.Azospirillum brasilense hisF gene, carrying a single nucleotide mutation that generates a stop codon between the two halves of the gene, was used to transform the histidine-auxotrophic Escherichia coli strain FB182 (hisF892). The transformed strain was subjected to selective pressure (i.e., low concentration/absence of histidine in the growth medium) and the obtained mutants were characterized.The restoration of prototrophy was strongly dependent on the time of incubation and on the strength of the selective pressure. The mutations involved the introduced stop codon with a single base substitution and none of the mutants restored the wild-type codon. Possible correlations between the different mutations and i) E. coli codon usage, ii) three-dimensional structures of the mutated HisF proteins, and iii) growth ability of the mutants were investigated. On the contrary, when the experiment was repeated by mutating a more conserved codon, only a synonymous substitution was obtained.Thus, experiments performed in this study allowed to mimic a possible gene elongation event occurred during the evolution of hisF gene, evidencing the ability of bacterial cells to modify their genome in short times under selective conditions.

Abstract PR005: Multi-level functional genomics reveals molecular and cellular oncogenicity of patient-based 3’ untranslated region mutations

Abstract It is clear from the high mortality rate in advanced cancers that we need a more complete understanding of the many sources of oncogenic dysregulation in tumors. The importance of the non-coding genome in disease has recently begun to be revealed through expanded use of whole genome sequencing. In particular, 3’ untranslated region (3’UTR) somatic mutations represent an important but largely unexplored avenue of alternative oncogenic gene dysregulation. Individual instances of 3’UTR-mediated oncogenicity are known, including oncogenic RNA binding proteins and microRNAs. However, a comprehensive, high-throughput study of patient-based 3’UTR mutations and their effect on post-transcriptional gene regulation in cancer has yet to be undertaken. To determine the significance of 3’UTR mutations in advanced disease, we identify 3’UTR somatic variants across 185 metastatic castration-resistant prostate tumors, discovering 14,497 single-nucleotide mutations, which are enriched in oncogenic pathways and 3’UTR regulatory elements. We develop two complementary massively parallel reporter assays (MPRAs) to measure how thousands of these patient-based mutations affect two distinct levels of post-transcriptional gene regulation: mRNA translation and stability. These MPRAs identify hundreds of functional variants that allow us to define three determinants of mutation significance: sequence conservation, known regulatory elements, and RNA structure. Furthermore, we demonstrate the clinical relevance of these mutations, observing that CRISPR-Cas9 base editing of distinct patient 3’UTR mutations into endogenous cellular loci increases oncogenic mRNA translation and cellular stress resistance. Finally, we go back to patients, illustrating that those harboring oncogenic 3’UTR mutations discovered by our methods display particularly poor prognosis. This work represents an unprecedented view of the extent to which disease-relevant 3’UTR mutations affect mRNA stability, translation, and cancer progression, uncovering principles of regulatory functionality and potential therapeutic targets in previously unexplored regulatory regions. Citation Format: Samantha L. Schuster, Sonali Arora, Cynthia L. Wladyka, Lukas Corey, Bethany L. Stackhouse, Lori Kollath, Eva Corey, Lawrence D. True, Dave Young, Patrick J. Paddison, Andrew C. Hsieh. Multi-level functional genomics reveals molecular and cellular oncogenicity of patient-based 3’ untranslated region mutations [abstract]. In: Proceedings of the AACR Special Conference: Advances in Prostate Cancer Research; 2023 Mar 15-18; Denver, Colorado. Philadelphia (PA): AACR; Cancer Res 2023;83(11 Suppl):Abstract nr PR005.

Multiplexed and accurate quantification strategy for miRNA based on specific terminal-mediated PCR with equivalent amplification.

MicroRNAs (miRNAs) are recognized as potential biomarkers for the early diagnosis and prognosis of different diseases. Multiplexed and accurate miRNA quantification methods with equivalent detection efficiency are particularly crucial due to their complex biological functions and lack of a unified internal reference gene. Here, a unique multiplexed miRNA detection method, named Specific Terminal-Mediated miRNA PCR (STEM-Mi-PCR), was developed. It mainly includes a linear reverse transcription step using tailored-designed target specific capture primers, followed by an exponential amplification process using two universal primers to execute the multiplex assay. For proof of concept, four miRNAs were used as models to develop a multiplexed detection assay within one tube simultaneously and then evaluate the performance of the established STEM-Mi-PCR. The sensitivity of the 4-plexed assay was approximately 100 aM with an equivalent amplification efficiency (95.67±8.58%), and had no cross-reactivity each other with high specificity. Quantification of different miRNAs in twenty patients' tissues shown variation from approximately pM to fM concentration level, demonstrating the possibility of practical application of the established method. Moreover, this method was extraordinarily capable of single nucleotide mutation discrimination in different let-7 family members with no more than 0.7% nonspecific detection signal. Hence, the STEM-Mi-PCR we proposed here paves an easy and promising way for miRNA profiling in future clinical applications.

Electrical Detection Assay Based on Programmable Nucleic Acid Probe for Efficient Single-Nucleotide Polymorphism Identification.

The large-scale pandemic and fast evolution of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants have triggered an urgent need for an efficient and sensitive on-site nucleic acid testing method with single-nucleotide polymorphism (SNP) identification capability. Here, we report a multiplexed electrical detection assay based on a paperclip-shaped nucleic acid probe (PNprobe) functionalized field-effect transistor (FET) biosensor for highly sensitive and specific detection and discrimination of SARS-CoV-2 variants. The three-stem structure of the PNprobe significantly amplifies the thermodynamic stability difference between variant RNAs that differ in a single-nucleotide mutation. With the assistance of combinatorial FET detection channels, the assay realizes simultaneously the detection and identification of key mutations of seven SARS-CoV-2 variants, including nucleotide substitutions and deletions at single-nucleotide resolution within 15 min. For 70 simulated throat swab samples, the multiplexed electrical detection assay shows an identification accuracy of 97.1% for the discrimination of SARS-CoV-2 variants. Our designed multiplexed electrical detection assay with SNP identification capability provides an efficient tool to achieve scalable pandemic screening.

Delta SARS-CoV-2 s2m Structure, Dynamics, and Entropy: Consequences of the G15U Mutation.

Bioinformatic analysis of the Delta SARS-CoV-2 genome reveals a single nucleotide mutation (G15U) in the stem-loop II motif (s2m) relative to ancestral SARS-CoV-2. Despite sequence similarity, unexpected differences between SARS-CoV-2 and Delta SARS-CoV-2 s2m homodimerization experiments require the discovery of unknown structural and thermodynamic changes necessary to rationalize the data. Using our reported SARS-CoV-2 s2m model, we induced the G15U substitution and performed 3.5 microseconds of unbiased molecular dynamics simulation at 283 and 310 K. The resultant Delta s2m adopted a secondary structure consistent with our reported NMR data, resulting in significant deviations in the tertiary structure and dynamics from our SARS-CoV-2 s2m model. First, we find differences in the overall three-dimensional structure, where the characteristic 90° L-shaped kink of the SARS-CoV-2 s2m did not form in the Delta s2m resulting in a "linear" hairpin with limited bending dynamics. Delta s2m helical parameters are calculated to align closely with A-form RNA, effectively eliminating a hinge point to form the L-shape kink by correcting an upper stem defect in SARS-CoV-2 induced by a noncanonical and dynamic G:A base pair. Ultimately, the shape difference rationalizes the migration differences in reported electrophoresis experiments. Second, increased fluctuation of the Delta s2m palindromic sequence, within the terminal loop, compared to SARS-CoV-2 s2m results in an estimated increase of entropy of 6.8 kcal/mol at 310 K relative to the SARS-CoV-2 s2m. The entropic difference offers a unique perspective on why the Delta s2m homodimerizes less spontaneously, forming fewer kissing dimers and extended duplexes compared to SARS-CoV-2. In this work, both the L-shape reduction and palindromic entropic penalty provides an explanation of our reported in vitro electrophoresis homodimerization results. Ultimately, the structural, dynamical, and entropic differences between the SARS-CoV-2 s2m and Delta s2m serve to establish a foundation for future studies of the s2m function in the viral lifecycle.

Genetic Diversity among Rose Rosette Virus Isolates: A Roadmap towards Studies of Gene Function and Pathogenicity.

The phylogenetic relationships of ninety-five rose rosette virus (RRV) isolates with full-length genomic sequences were analyzed. These isolates were recovered mostly from commercial roses that are vegetatively propagated rather than grown from seed. First, the genome segments were concatenated, and the maximum likelihood (ML) tree shows that the branches arrange independent of their geographic origination. There were six major groups of isolates, with 54 isolates in group 6 and distributed in two subgroups. An analysis of nucleotide diversity across the concatenated isolates showed lower genetic differences among RNAs encoding the core proteins required for encapsidation than the latter genome segments. Recombination breakpoints were identified near the junctions of several genome segments, suggesting that the genetic exchange of segments contributes to differences among isolates. The ML analysis of individual RNA segments revealed different relationship patterns among isolates, which supports the notion of genome reassortment. We tracked the branch positions of two newly sequenced isolates to highlight how genome segments relate to segments of other isolates. RNA6 has an interesting pattern of single-nucleotide mutations that appear to influence amino acid changes in the protein products derived from ORF6a and ORF6b. The P6a proteins were typically 61 residues, although three isolates encoded P6a proteins truncated to 29 residues, and four proteins extended 76-94 residues. Homologous P5 and P7 proteins appear to be evolving independently. These results suggest greater diversity among RRV isolates than previously recognized.

MLH1 Loss in Primary Prostate Cancer.

Among mismatch repair-deficient (MMRd) prostate cancers, loss of MLH1 is relatively uncommon and few cases have been reported in detail. Here, we describe the molecular features of two cases of primary prostate cancer with MLH1 loss detected by immunohistochemistry, and in one case, confirmed via transcriptomic profiling. Both cases were microsatellite stable on standard polymerase chain reaction (PCR)-based microsatellite instability (MSI) testing, but showed evidence of MSI on a newer PCR-based long mononucleotide repeat (LMR) assay and by next-generation sequencing. Germline testing was negative for Lynch syndrome-associated mutations in both cases. Targeted or whole-exome tumor sequencing using multiple commercial/academic platforms (Foundation, Tempus, JHU, and UW-OncoPlex) showed modestly elevated, though variable, tumor mutation burden estimates (2.3-10 mutations/Mb) consistent with MMRd, but without identifiable pathogenic single-nucleotide or indel mutations in MLH1. Copy-number analysis confirmed biallelic MLH1 loss in one case and monoallelic MLH1 loss in the second case, without evidence of MLH1 promoter hypermethylation in either. The second patient was treated with single-agent pembrolizumab and demonstrated a short-lived prostate-specific antigen response. These cases highlight the challenges in identifying MLH1-deficient prostate cancers using standard MSI testing and commercial sequencing panels, and support the utility of immunohistochemical assays and LMR- or sequencing-based MSI testing for detection of MMRd prostate cancers.

PO-01-123 THE BIOPHYSICAL ACTION OF MICRORNA IS ESSENTIAL TO MAINTAIN THE NORMAL FUNCTION AND ELECTROPHYSIOLOGY OF THE HEART

MicroRNAs (miRs) regulate most (if not all) biological events via a broadly-recoginized RNA interference (RNAi) mechnism. miR1, encoded by miR1-1 and miR1-2 genes, is the predominant miRs of the heart and plays critical roles in heart development and cardiac diseases. Previous studies of miRs have focused on the canonical RNAi mechanism. Recently, we discovered a novel biophysical action of miR and found that miR1 directly binds to ion channels physically modulates the electrophysiology of cardiomyocytes. There is an arrhythmia-associated human single nucleotide polymorphism (hSNP)–hSNP14A/G (rs776480338) on mature miR1, in which the 14th nucleotide of miR1 is mutated from A to G. We found that hSNP14A/G specifically disrupts this biophysical action while maintaining its RNAi function. To define the specific physiological significance of the biophysical action of miR1 in regulating the function of the heart. By using CRISPR/Cas9 technology, we created single nucleotide mutation of hSNP14A/G on mouse miR1-1 and miR1-2 genes and have successfully developed 14G-mutated homozygous transgenic mice (14G-Homo, miR1-114G/G;miR1-214G/G), in which the biophysical function of miR1 is specifically removed while the RNAi function is maintained. The function of whole heart was evaluated by echocardiagram and surface ECG. Cellular electrophysiology and contraction of isolated cardiomyocytes were studied by patch clamping and IonOptix. 14G-Homo mice had higher mortality than wildtype (WT) mice. Compared to WT hearts, adult 14G-homo heart (2-month age) had significantly lower ejection fraction and fractional shortening with thinner wall of left ventricle, which became more severe at the age of 6 months; some mice developed severe heart failure. Surface ECG showed that 14G-Homo hearts had slower heart rate and the prolonged QRS interval and QT interval. Optial mapping revealed a slower conduction velocity in 14G-Homo heart with a high inducibility of ventricular arrhythmias. Patch clamping of isolated ventricular cardiomyocytes showed that, compared to WT cells, 14G-Homo cardiomyocytes had a prolonged duration to action potention with DADs. Compared to WT cells, 14G-Homo myocytes had bigger IK1 and IL-Ca2+ and showed faster Ca2+ transients. A specific delection of miR’s biophysical fuction leads to arrhythmia and heart failure. Beyond the RNAi, miRs are ion channel modulators; the biophysical action of miRs is deeply involved in maintaining the homeostasis of the heart.

Evaluation of a Novel FKS1 R1354H Mutation Associated with Caspofungin Resistance in Candida auris Using the CRISPR-Cas9 System.

Outbreaks of invasive infections, with high mortality rates, caused by multidrug-resistant Candida auris have been reported worldwide. Although hotspot mutations in FKS1 are an established cause of echinocandin resistance, the actual contribution of these mutations to echinocandin resistance remains unknown. Here, we sequenced the FKS1 gene of a caspofungin-resistant clinical isolate (clade I) and identified a novel resistance mutation (G4061A inducing R1354H). We applied the clustered regularly interspaced short palindromic repeats (CRISPR)-Cas9 system to generate a recovered strain (H1354R) in which only this single nucleotide mutation was reverted to its wild-type sequence. We also generated mutant strains with only the R1354H mutation introduced into C. auris wild-type strains (clade I and II) and analyzed their antifungal susceptibility. Compared to their parental strains, the R1354H mutants exhibited a 4- to 16-fold increase in caspofungin minimum inhibitory concentration (MIC) while the H1354R reverted strain exhibited a 4-fold decrease in caspofungin MIC. In a mouse model of disseminated candidiasis, the in vivo therapeutic effect of caspofungin was more closely related to the FKS1 R1354H mutation and the virulence of the strain than its in vitro MIC. The CRISPR-Cas9 system could thus aid in elucidating the mechanism underlying drug resistance in C. auris.

A natural single nucleotide mutation in the small regulatory RNA ArcZ of Dickeya solani switches off the antimicrobial activities against yeast and bacteria.

The necrotrophic plant pathogenic bacterium Dickeya solani emerged in the potato agrosystem in Europe. All isolated strains of D. solani contain several large polyketide synthase/non-ribosomal peptide synthetase (PKS/NRPS) gene clusters. Analogy with genes described in other bacteria suggests that the clusters ooc and zms are involved in the production of secondary metabolites of the oocydin and zeamine families, respectively. A third cluster named sol was recently shown to produce an antifungal molecule. In this study, we constructed mutants impaired in each of the three secondary metabolite clusters sol, ooc, and zms to compare first the phenotype of the D. solani wild-type strain D s0432-1 with its associated mutants. We demonstrated the antimicrobial functions of these three PKS/NRPS clusters against bacteria, yeasts or fungi. The cluster sol, conserved in several other Dickeya species, produces a secondary metabolite inhibiting yeasts. Phenotyping and comparative genomics of different D. solani wild-type isolates revealed that the small regulatory RNA ArcZ plays a major role in the control of the clusters sol and zms. A single-point mutation, conserved in some Dickeya wild-type strains, including the D. solani type strain IPO 2222, impairs the ArcZ function by affecting its processing into an active form.

Abstract 2553: Elucidation of the downstream signaling of PIK3CA genomic alterations in serous ovarian cancer

Abstract Genomic aberrations of PIK3CA, which encodes the catalytic subunit p110α of the class IA phosphatidylinositol 3-kinase (PI3K), is frequently observed in cancers. In ovarian cancer of the serous subtype (the most common histotype of the disease), PIK3CA amplification is dominant over single nucleotide mutations in multiple patient cohorts, suggesting the possibility that PIK3CA amplification is an alternative to PIK3CA “hot spot” mutation for serous ovarian tumorigenesis. However, whether PIK3CA amplification and mutation exert the same functional consequences on serous ovarian tumorigenesis remains uncertain. This study aims to explore and compare the functional roles between the two types of PIK3CA gene aberrations on tumorigenicity, cellular signaling, and therapeutic responses in serous ovarian cancer. We confirmed that both PIK3CA amplification and E545K mutation could induce multiple tumorigenic phenotypes in CRISPR/Cas9-mediated E545K knock-in or wild-type PIK3CA overexpressing serous ovarian cancer cell lines. Signaling pathway analyses showed that E545K knock-in cells, but not PIK3CA-overexpressing cells, had the MAPK signaling pathway activated and were sensitive to MAPK pathway inhibition. An induction of AKT phosphorylation and PDK1 phosphorylation was observed in both E545K knock-in cells and wild-type PIK3CA-overexpressing cells. Accordingly, these cells displayed higher sensitivity towards p110α or AKT inhibitors. Intriguingly, in line with previous report which demonstrated YAP as AKT substrate, the activation of AKT associated with an increase in YAP S127 phosphorylation in the knock-in or PIK3CA-overexpressing cells. AKT inhibition suppressed YAP phosphorylation. Corroborating with the role of S127 in nuclear exclusion of YAP, YAP accumulated in the cytoplasm of these cells, leading to reduced transcriptional activity of YAP. While the functional significance of YAP cytoplasmic accumulation is yet to be investigated, our results reveal the common and distinct downstream pathways between PIK3CA amplification and mutation. These findings may be applicable to p110α-targeting therapies in serous ovarian cancer. Citation Format: Shibo Zhang, Hei Ip Hong, Victor CY Mak, Yuan Zhou, Panpan Wang, Lydia Wai Ting Cheung. Elucidation of the downstream signaling of PIK3CA genomic alterations in serous ovarian cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 2553.

Evolution of the SARS-CoV-2 Mutational Spectrum.

SARS-CoV-2 evolves rapidly in part because of its high mutation rate. Here, we examine whether this mutational process itself has changed during viral evolution. To do this, we quantify the relative rates of different types of single-nucleotide mutations at 4-fold degenerate sites in the viral genome across millions of human SARS-CoV-2 sequences. We find clear shifts in the relative rates of several types of mutations during SARS-CoV-2 evolution. The most striking trend is a roughly 2-fold decrease in the relative rate of G→T mutations in Omicron versus early clades, as was recently noted by Ruis et al. (2022. Mutational spectra distinguish SARS-CoV-2 replication niches. bioRxiv, doi:10.1101/2022.09.27.509649). There is also a decrease in the relative rate of C→T mutations in Delta, and other subtle changes in the mutation spectrum along the phylogeny. We speculate that these changes in the mutation spectrum could arise from viral mutations that affect genome replication, packaging, and antagonization of host innate-immune factors, although environmental factors could also play a role. Interestingly, the mutation spectrum of Omicron is more similar than that of earlier SARS-CoV-2 clades to the spectrum that shaped the long-term evolution of sarbecoviruses. Overall, our work shows that the mutation process is itself a dynamic variable during SARS-CoV-2 evolution and suggests that human SARS-CoV-2 may be trending toward a mutation spectrum more similar to that of other animal sarbecoviruses.

Epistasis decreases with increasing antibiotic pressure but not temperature.

Predicting mutational effects is essential for the control of antibiotic resistance (ABR). Predictions are difficult when there are strong genotype-by-environment (G × E), gene-by-gene (G × G or epistatic) or gene-by-gene-by-environment (G × G × E) interactions. We quantified G × G × E effects in Escherichia coli across environmental gradients. We created intergenic fitness landscapes using gene knock-outs and single-nucleotide ABR mutations previously identified to vary in the extent of G × E effects in our environments of interest. Then, we measured competitive fitness across a complete combinatorial set of temperature and antibiotic dosage gradients. In this way, we assessed the predictability of 15 fitness landscapes across 12 different but related environments. We found G × G interactions and rugged fitness landscapes in the absence of antibiotic, but as antibiotic concentration increased, the fitness effects of ABR genotypes quickly overshadowed those of gene knock-outs, and the landscapes became smoother. Our work reiterates that some single mutants, like those conferring resistance or susceptibility to antibiotics, have consistent effects across genetic backgrounds in stressful environments. Thus, although epistasis may reduce the predictability of evolution in benign environments, evolution may be more predictable in adverse environments. This article is part of the theme issue 'Interdisciplinary approaches to predicting evolutionary biology'.