N-Terminal Capping of the αO-Conotoxin Analogue GeX-2 Improves the Serum Stability and Selectivity toward the Human α9α10 Nicotinic Acetylcholine Receptor.

α9α10 nicotinic acetylcholine receptors (nAChRs) are a promising nonopioid analgesic target, with α9α10 nAChR antagonists showing efficacy against chemotherapy-induced hyperalgesia and allodynia. GeX-2, a potent analgesic conotoxin antagonist of α9α10 nAChRs, has limited serum stability. This study improved GeX-2 stability by capping its N-terminal with fatty acids or polyethylene glycol chains, which enhanced its serum stability but eliminated activity at G protein-coupled γ-aminobutyric acid type B (GABAB) receptor-coupled CaV2.2 channels while preserving activity at α9α10 nAChRs. In vivo, α9α10 nAChRs antagonism alone did not alleviate neuropathic pain, highlighting the importance of GABAB receptor-coupled CaV2.2 channels in GeX-2's antinociceptive effects in the chronic constriction injury rat model. The GeX-2 analogue, with an N-terminal methyl group, showed improved activity and selectivity for α9α10 nAChRs, increased serum half-life, and strong analgesic effects in oxaliplatin-induced cold allodynia models. AlphaFold3 and molecular dynamics simulations provided insights into the binding modes and the effects of N-terminal capping, which informed future peptide therapeutic developments.

High-resolution functional mapping of RAD51C by saturation genome editing.

Pathogenic variants in RAD51C confer an elevated risk of breast and ovarian cancer, while individuals homozygous for specific RAD51C alleles may develop Fanconi anemia. Using saturation genome editing (SGE), we functionally assess 9,188 unique variants, including >99.5% of all possible coding sequence single-nucleotide alterations. By computing changes in variant abundance and Gaussian mixture modeling (GMM), we functionally classify 3,094 variants to be disruptive and use clinical truth sets to reveal an accuracy/concordance of variant classification >99.9%. Cell fitness was the primary assay readout allowing us to observe a phenomenon where specific missense variants exhibit distinct depletion kinetics potentially suggesting that they represent hypomorphic alleles. We further explored our exhaustive functional map, revealing critical residues on the RAD51C structure and resolving variants found in cancer-segregating kindred. Furthermore, through interrogation of UK Biobank and a large multi-center ovarian cancer cohort, we find significant associations between SGE-depleted variants and cancer diagnoses.

Impact of essential genes on the success of genome editing experiments generating 3313 new genetically engineered mouse lines.

The International Mouse Phenotyping Consortium (IMPC) systematically produces and phenotypes mouse lines with presumptive null mutations to provide insight into gene function. The IMPC now uses the programmable RNA-guided nuclease Cas9 for its increased capacity and flexibility to efficiently generate null alleles in the C57BL/6N strain. In addition to being a valuable novel and accessible research resource, the production of 3313 knockout mouse lines using comparable protocols provides a rich dataset to analyze experimental and biological variables affecting in vivo gene engineering with Cas9. Mouse line production has two critical steps - generation of founders with the desired allele and germline transmission (GLT) of that allele from founders to offspring. A systematic evaluation of the variables impacting success rates identified gene essentiality as the primary factor influencing successful production of null alleles. Collectively, our findings provide best practice recommendations for using Cas9 to generate alleles in mouse essential genes, many of which are orthologs of genes linked to human disease.

Saturation genome editing of BAP1 functionally classifies somatic and germline variants.

Many variants that we inherit from our parents or acquire de novo or somatically are rare, limiting the precision with which we can associate them with disease. We performed exhaustive saturation genome editing (SGE) of BAP1, the disruption of which is linked to tumorigenesis and altered neurodevelopment. We experimentally characterized 18,108 unique variants, of which 6,196 were found to have abnormal functions, and then used these data to evaluate phenotypic associations in the UK Biobank. We also characterized variants in a large population-ascertained tumor collection, in cancer pedigrees and ClinVar, and explored the behavior of cancer-associated variants compared to that of variants linked to neurodevelopmental phenotypes. Our analyses demonstrated that disruptive germline BAP1 variants were significantly associated with higher circulating levels of the mitogen IGF-1, suggesting a possible pathological mechanism and therapeutic target. Furthermore, we built a variant classifier with >98% sensitivity and specificity and quantify evidence strengths to aid precision variant interpretation.

IFI207, a young and fast-evolving protein, controls retroviral replication via the STING pathway.

Mammalian AIM-2-like receptor (ALR) proteins bind nucleic acids and initiate production of type I interferons or inflammasome assembly, thereby contributing to host innate immunity. In mice, the Alr locus is highly polymorphic at the sequence and copy number level, and we show here that it is one of the most dynamic regions of the genome. One rapidly evolving gene within this region, Ifi207, was introduced to the Mus genome by gene conversion or an unequal recombination event a few million years ago. Ifi207 has a large, distinctive repeat region that differs in sequence and length among Mus species and even closely related inbred Mus musculus strains. We show that IFI207 controls murine leukemia virus (MLV) infection in vivo and that it plays a role in the STING-mediated response to cGAMP, dsDNA, DMXXA, and MLV. IFI207 binds to STING, and inclusion of its repeat region appears to stabilize STING protein. The Alr locus and Ifi207 provide a clear example of the evolutionary innovation of gene function, possibly as a result of host-pathogen co-evolution.IMPORTANCEThe Red Queen hypothesis predicts that the arms race between pathogens and the host may accelerate evolution of both sides, and therefore causes higher diversity in virulence factors and immune-related proteins, respectively . The Alr gene family in mice has undergone rapid evolution in the last few million years and includes the creation of two novel members, MndaL and Ifi207. Ifi207, in particular, became highly divergent, with significant genetic changes between highly related inbred mice. IFI207 protein acts in the STING pathway and contributes to anti-retroviral resistance via a novel mechanism. The data show that under the pressure of host-pathogen coevolution in a dynamic locus, gene conversion and recombination between gene family members creates new genes with novel and essential functions that play diverse roles in biological processes.

Variant ranking pipeline for complex familial disorders.

Identifying genetic susceptibility factors for complex disorders remains a challenging task. To analyze collections of small and large pedigrees where genetic heterogeneity is likely, but biological commonalities are plausible, we have developed a weights-based pipeline to prioritize variants and genes. The Weights-based vAriant Ranking in Pedigrees (WARP) pipeline prioritizes variants using 5 weights: disease incidence rate, number of cases in a family, genome fraction shared amongst cases in a family, allele frequency and variant deleteriousness. Weights, except for the population allele frequency weight, are normalized between 0 and 1. Weights are combined multiplicatively to produce family-specific-variant weights that are then averaged across all families in which the variant is observed to generate a multifamily weight. Sorting multifamily weights in descending order creates a ranked list of variants and genes for further investigation. WARP was validated using familial melanoma sequence data from the European Genome-phenome Archive. The pipeline identified variation in known germline melanoma genes POT1, MITF and BAP1 in 4 out of 13 families (31%). Analysis of the other 9 families identified several interesting genes, some of which might have a role in melanoma. WARP provides an approach to identify disease predisposing genes in studies with small and large pedigrees.

Intratumoral antigen signaling traps CD8+ T cells to confine exhaustion to the tumor site.

Immunotherapy advances have been hindered by difficulties in tracking the behaviors of lymphocytes after antigen signaling. Here, we assessed the behavior of T cells active within tumors through the development of the antigen receptor signaling reporter (AgRSR) mouse, fate-mapping lymphocytes responding to antigens at specific times and locations. Contrary to reports describing the ready egress of T cells out of the tumor, we find that intratumoral antigen signaling traps CD8+ T cells in the tumor. These clonal populations expand and become increasingly exhausted over time. By contrast, antigen-signaled regulatory T cell (Treg) clonal populations readily recirculate out of the tumor. Consequently, intratumoral antigen signaling acts as a gatekeeper to compartmentalize CD8+ T cell responses, even within the same clonotype, thus enabling exhausted T cells to remain confined to a specific tumor tissue site.

NUDCD3 deficiency disrupts V(D)J recombination to cause SCID and Omenn syndrome.

Inborn errors of T cell development present a pediatric emergency in which timely curative therapy is informed by molecular diagnosis. In 11 affected patients across four consanguineous kindreds, we detected homozygosity for a single deleterious missense variant in the gene NudC domain-containing 3 (NUDCD3). Two infants had severe combined immunodeficiency with the complete absence of T and B cells (T -B- SCID), whereas nine showed classical features of Omenn syndrome (OS). Restricted antigen receptor gene usage by residual T lymphocytes suggested impaired V(D)J recombination. Patient cells showed reduced expression of NUDCD3 protein and diminished ability to support RAG-mediated recombination in vitro, which was associated with pathologic sequestration of RAG1 in the nucleoli. Although impaired V(D)J recombination in a mouse model bearing the homologous variant led to milder immunologic abnormalities, NUDCD3 is absolutely required for healthy T and B cell development in humans.

Tabernanthalog and ibogainalog inhibit the α7 and α9α10 nicotinic acetylcholine receptors via different mechanisms and with higher potency than the GABAA receptor and CaV2.2 channel.

In this study, we have investigated the pharmacological activity and structural interaction of two novel psychoplastogens, tabernanthalog (TBG) and ibogainalog (IBG) at heterologously-expressed rat (r) and human (h) nicotinic acetylcholine receptors (nAChRs), the rα1ß2γ2L γ-aminobutyric acid type A receptor (GABAAR), and the human voltage-gated N-type calcium channel (CaV2.2 channel). Both compounds inhibited the nAChRs with the following receptor selectivity: α9α10 > α7 > α3ß2 â‰… α3ß4, indicating that ß2/ß4 subunits are relatively less important for their activity. The potencies of TBG and IBG were comparable at hα7 and hα9α10 subtypes, and comparable to their rat counterparts. TBG- and IBG-induced inhibition of rα7 was ACh concentration-independent and voltage-dependent, whereas rα9α10 inhibition was ACh concentration-dependent and voltage-independent, suggesting that they interact with the α7 ion channel pore and α9α10 orthosteric ligand binding site, respectively. These results were supported by molecular docking studies showing that at the α7 model TBG forms stable interactions with luminal rings at 9', 13', and 16', whereas IBG mostly interacts with the extracellular-transmembrane junction. In the α9α10 model, however, these compounds interacted with several residues from the principal (+) and complementary (-) sides in the transmitter binding site. Ibogaminalog (DM506) also interacted with a non-luminal site at α7, and one α9α10 orthosteric site. TBG and IBG inhibited the GABAAR and CaV2.2 channels with 10 to 30-fold lower potencies. In sum, we show that TBG and IBG inhibit the α7 and α9α10 nAChRs by noncompetitive and competitive mechanisms, respectively, and with higher potency than the GABAAR and CaV2.2 channel.

A gut microbial signature for combination immune checkpoint blockade across cancer types.

Immune checkpoint blockade (ICB) targeting programmed cell death protein 1 (PD-1) and cytotoxic T lymphocyte protein 4 (CTLA-4) can induce remarkable, yet unpredictable, responses across a variety of cancers. Studies suggest that there is a relationship between a cancer patient's gut microbiota composition and clinical response to ICB; however, defining microbiome-based biomarkers that generalize across cohorts has been challenging. This may relate to previous efforts quantifying microbiota to species (or higher taxonomic rank) abundances, whereas microbial functions are often strain specific. Here, we performed deep shotgun metagenomic sequencing of baseline fecal samples from a unique, richly annotated phase 2 trial cohort of patients with diverse rare cancers treated with combination ICB (n = 106 discovery cohort). We demonstrate that strain-resolved microbial abundances improve machine learning predictions of ICB response and 12-month progression-free survival relative to models built using species-rank quantifications or comprehensive pretreatment clinical factors. Through a meta-analysis of gut metagenomes from a further six comparable studies (n = 364 validation cohort), we found cross-cancer (and cross-country) validity of strain-response signatures, but only when the training and test cohorts used concordant ICB regimens (anti-PD-1 monotherapy or combination anti-PD-1 plus anti-CTLA-4). This suggests that future development of gut microbiome diagnostics or therapeutics should be tailored according to ICB treatment regimen rather than according to cancer type.

Rational Design of Potent α-Conotoxin PeIA Analogues with Non-Natural Amino Acids for the Inhibition of Human α9α10 Nicotinic Acetylcholine Receptors.

α-Conotoxins (α-CTxs) are structurally related peptides that antagonize nicotinic acetylcholine receptors (nAChRs), which may serve as new alternatives to opioid-based treatment for pain-related conditions. The non-natural amino acid analogues of α-CTxs have been demonstrated with improved potency compared to the native peptide. In this study, we chemically synthesized Dab/Dap-substituted analogues of α-CTx PeIA and evaluated their activity at heterologously expressed human α9α10 nAChRs. PeIA[S4Dap, S9Dap] had the most potent half-maximal inhibitory concentration (IC50) of 0.93 nM. Molecular dynamic simulations suggested that the side chain amino group of Dap4 formed additional hydrogen bonds with S168 and D169 of the receptor and Dap9 formed an extra hydrogen bond interaction with Q34, which is distinctive to PeIA. Overall, our findings provide new insights into further development of more potent analogues of α-CTxs, and PeIA[S4Dap, S9Dap] has potential as a drug candidate for the treatment of chronic neuropathic pain.

NFIL3 contributes to cytotoxic T lymphocyte-mediated killing.

Cytotoxic T lymphocytes (CTLs) are key effectors of the adaptive immune system that recognize and eliminate virally infected and cancerous cells. In naive CD8+ T cells, T-cell receptor (TCR) engagement drives a number of transcriptional, translational and proliferation changes over the course of hours and days leading to differentiation into CTLs. To gain a better insight into this mechanism, we compared the transcriptional profiles of naive CD8+ T cells to those of activated CTLs. To find new regulators of CTL function, we performed a selective clustered regularly interspaced short palindromic repeats (CRISPR) screen on upregulated genes and identified nuclear factor IL-3 (NFIL3) as a potential regulator of cytotoxicity. Although NFIL3 has established roles in several immune cells including natural killer, Treg, dendritic and CD4+ T cells, its function in CD8+ CTLs is less well understood. Using CRISPR/Cas9 editing, we found that removing NFIL3 in CTLs resulted in a marked decrease in cytotoxicity. We found that in CTLs lacking NFIL3 TCR-induced extracellular signal-regulated kinase phosphorylation, immune synapse formation and granule release were all intact while cytotoxicity was functionally impaired in vitro. Strikingly, NFIL3 controls the production of cytolytic proteins as well as effector cytokines. Thus, NFIL3 plays a cell intrinsic role in modulating cytolytic mechanisms in CTLs.

Targeted action to increase inclusion at the Wellcome Sanger Institute.

The Sanger Excellence Fellowship has been established to increase the representation of researchers with Black-heritage backgrounds at a leading research centre in the UK.

Dual Antagonism of α9α10 nAChR and GABAB Receptor-Coupled CaV2.2 Channels by an Analgesic αO-Conotoxin Analogue.

Pain severely affects the physical and mental health of patients. The need to develop nonopioid analgesic drugs to meet medical demands is urgent. In this study, we designed a truncated analogue of αO-conotoxin, named GeX-2, based on disulfide-bond deletion and sequence truncation. GeX-2 retained the potency of its parent peptide at the human α9α10 nAChR and exhibited potent inhibitory activity at CaV2.2 channels via activation of the GABAB receptor (GABABR). Importantly, GeX-2 significantly alleviated pain in the rat model of chronic constriction injury. The dual inhibition of GeX-2 at both α9α10 nAChRs and CaV2.2 channels is speculated to synergistically mediate the potent analgesic effects. Results from site-directed mutagenesis assay and computational modeling suggest that GeX-2 preferentially interacts with the α10(+)α10(-) binding site of α9α10 nAChR and favorably binds to the top region of the GABABR2 subunit. The study offers vital insights into the molecular action mechanism of GeX-2, demonstrating its potential as a novel nonopioid analgesic.

Modeling methyl-sensitive transcription factor motifs with an expanded epigenetic alphabet.

BACKGROUND: Transcription factors bind DNA in specific sequence contexts. In addition to distinguishing one nucleobase from another, some transcription factors can distinguish between unmodified and modified bases. Current models of transcription factor binding tend not to take DNA modifications into account, while the recent few that do often have limitations. This makes a comprehensive and accurate profiling of transcription factor affinities difficult. RESULTS: Here, we develop methods to identify transcription factor binding sites in modified DNA. Our models expand the standard A/C/G/T DNA alphabet to include cytosine modifications. We develop Cytomod to create modified genomic sequences and we also enhance the MEME Suite, adding the capacity to handle custom alphabets. We adapt the well-established position weight matrix (PWM) model of transcription factor binding affinity to this expanded DNA alphabet. Using these methods, we identify modification-sensitive transcription factor binding motifs. We confirm established binding preferences, such as the preference of ZFP57 and C/EBPß for methylated motifs and the preference of c-Myc for unmethylated E-box motifs. CONCLUSIONS: Using known binding preferences to tune model parameters, we discover novel modified motifs for a wide array of transcription factors. Finally, we validate our binding preference predictions for OCT4 using cleavage under targets and release using nuclease (CUT&RUN) experiments across conventional, methylation-, and hydroxymethylation-enriched sequences. Our approach readily extends to other DNA modifications. As more genome-wide single-base resolution modification data becomes available, we expect that our method will yield insights into altered transcription factor binding affinities across many different modifications.

Interrater reliability of interictal EEG waveforms in Lennox-Gastaut Syndrome.

OBJECTIVE: Identification of EEG waveforms is critical for diagnosing Lennox-Gastaut Syndrome (LGS) but is complicated by the progressive nature of the disease. Here, we assess the interrater reliability (IRR) among pediatric epileptologists for classifying EEG waveforms associated with LGS. METHODS: A novel automated algorithm was used to objectively identify epochs of EEG with transient high power, which were termed events of interest (EOIs). The algorithm was applied to EEG from 20 LGS subjects and 20 healthy controls during NREM sleep, and 1350 EOIs were identified. Three raters independently reviewed the EOIs within isolated 15-second EEG segments in a randomized, blinded fashion. For each EOI, the raters assigned a waveform label (spike and slow wave, generalized paroxysmal fast activity, seizure, spindle, vertex, muscle, artifact, nothing, or other) and indicated the perceived subject type (LGS or control). RESULTS: Labeling of subject type had 85% accuracy across all EOIs and an IRR of κ =0.790, suggesting that brief segments of EEG containing high-power waveforms can be reliably classified as pathological or normal. Waveform labels were less consistent, with κ =0.558, and the results were highly variable for different categories of waveforms. Label mismatches typically occurred when one reviewer selected "nothing," suggesting that reviewers had different thresholds for applying named labels. SIGNIFICANCE: Classification of EEG waveforms associated with LGS has weak IRR, due in part to varying thresholds applied during visual review. Computational methods to objectively define EEG biomarkers of LGS may improve IRR and aid clinical decision-making.

Saturation genome editing of DDX3X clarifies pathogenicity of germline and somatic variation.

Loss-of-function of DDX3X is a leading cause of neurodevelopmental disorders (NDD) in females. DDX3X is also a somatically mutated cancer driver gene proposed to have tumour promoting and suppressing effects. We perform saturation genome editing of DDX3X, testing in vitro the functional impact of 12,776 nucleotide variants. We identify 3432 functionally abnormal variants, in three distinct classes. We train a machine learning classifier to identify functionally abnormal variants of NDD-relevance. This classifier has at least 97% sensitivity and 99% specificity to detect variants pathogenic for NDD, substantially out-performing in silico predictors, and resolving up to 93% of variants of uncertain significance. Moreover, functionally-abnormal variants can account for almost all of the excess nonsynonymous DDX3X somatic mutations seen in DDX3X-driven cancers. Systematic maps of variant effects generated in experimentally tractable cell types have the potential to transform clinical interpretation of both germline and somatic disease-associated variation.

Knockout or inhibition of USP30 protects dopaminergic neurons in a Parkinson's disease mouse model.

Mutations in SNCA, the gene encoding α-synuclein (αSyn), cause familial Parkinson's disease (PD) and aberrant αSyn is a key pathological hallmark of idiopathic PD. This α-synucleinopathy leads to mitochondrial dysfunction, which may drive dopaminergic neurodegeneration. PARKIN and PINK1, mutated in autosomal recessive PD, regulate the preferential autophagic clearance of dysfunctional mitochondria ("mitophagy") by inducing ubiquitylation of mitochondrial proteins, a process counteracted by deubiquitylation via USP30. Here we show that loss of USP30 in Usp30 knockout mice protects against behavioral deficits and leads to increased mitophagy, decreased phospho-S129 αSyn, and attenuation of SN dopaminergic neuronal loss induced by αSyn. These observations were recapitulated with a potent, selective, brain-penetrant USP30 inhibitor, MTX115325, with good drug-like properties. These data strongly support further study of USP30 inhibition as a potential disease-modifying therapy for PD.