Small Molecule Targets TMED9 and Promotes Lysosomal Degradation to Reverse Proteinopathy.

Intracellular accumulation of misfolded proteins causes toxic proteinopathies, diseases without targeted therapies. Mucin 1 kidney disease (MKD) results from a frameshift mutation in the MUC1 gene (MUC1-fs). Here, we show that MKD is a toxic proteinopathy. Intracellular MUC1-fs accumulation activated the ATF6 unfolded protein response (UPR) branch. We identified BRD4780, a small molecule that clears MUC1-fs from patient cells, from kidneys of knockin mice and from patient kidney organoids. MUC1-fs is trapped in TMED9 cargo receptor-containing vesicles of the early secretory pathway. BRD4780 binds TMED9, releases MUC1-fs, and re-routes it for lysosomal degradation, an effect phenocopied by TMED9 deletion. Our findings reveal BRD4780 as a promising lead for the treatment of MKD and other toxic proteinopathies. Generally, we elucidate a novel mechanism for the entrapment of misfolded proteins by cargo receptors and a strategy for their release and anterograde trafficking to the lysosome.

Tracking Cancer Evolution Reveals Constrained Routes to Metastases: TRACERx Renal.

Clear-cell renal cell carcinoma (ccRCC) exhibits a broad range of metastatic phenotypes that have not been systematically studied to date. Here, we analyzed 575 primary and 335 metastatic biopsies across 100 patients with metastatic ccRCC, including two cases sampledat post-mortem. Metastatic competence was afforded by chromosome complexity, and we identify 9p loss as a highly selected event driving metastasis and ccRCC-related mortality (p = 0.0014). Distinct patterns of metastatic dissemination were observed, including rapid progression to multiple tissue sites seeded by primary tumors of monoclonal structure. By contrast, we observed attenuated progression in cases characterized by high primary tumor heterogeneity, with metastatic competence acquired gradually and initial progression to solitary metastasis. Finally, we observed early divergence of primitive ancestral clones and protracted latency of up to two decades as a feature of pancreatic metastases.

An interdisciplinary consensus on the management of brain metastases in patients with renal cell carcinoma.

Brain metastases are a challenging manifestation of renal cell carcinoma. We have a limited understanding of brain metastasis tumor and immune biology, drivers of resistance to systemic treatment, and their overall poor prognosis. Current data support a multimodal treatment strategy with radiation treatment and/or surgery. Nonetheless, the optimal approach for the management of brain metastases from renal cell carcinoma remains unclear. To improve patient care, the authors sought to standardize practical management strategies. They performed an unstructured literature review and elaborated on the current management strategies through an international group of experts from different disciplines assembled via the network of the International Kidney Cancer Coalition. Experts from different disciplines were administered a survey to answer questions related to current challenges and unmet patient needs. On the basis of the integrated approach of literature review and survey study results, the authors built algorithms for the management of single and multiple brain metastases in patients with renal cell carcinoma. The literature review, consensus statements, and algorithms presented in this report can serve as a framework guiding treatment decisions for patients. CA Cancer J Clin. 2022;72:454-489.

Bipolar impact and phasing of Heinrich-type climate variability.

During the last ice age, the Laurentide Ice Sheet exhibited extreme iceberg discharge events that are recorded in North Atlantic sediments1. These Heinrich events have far-reaching climate impacts, including widespread disruptions to hydrological and biogeochemical cycles2-4. They occurred during Heinrich stadials-cold periods with strongly weakened Atlantic overturning circulation5-7. Heinrich-type variability is not distinctive in Greenland water isotope ratios, a well-dated site temperature proxy8, complicating efforts to assess their regional climate impact and phasing against Antarctic climate change. Here we show that Heinrich events have no detectable temperature impact on Greenland and cooling occurs at the onset of several Heinrich stadials, and that both types of Heinrich variability have a distinct imprint on Antarctic climate. Antarctic ice cores show accelerated warming that is synchronous with increases in methane during Heinrich events, suggesting an atmospheric teleconnection9, despite the absence of a Greenland climate signal. Greenland ice-core nitrogen stable isotope ratios, a sensitive temperature proxy, indicate an abrupt cooling of about three degrees Celsius at the onset of Heinrich Stadial 1 (17.8 thousand years before present, where present is defined as 1950). Antarctic warming lags this cooling by 133 ± 93 years, consistent with an oceanic teleconnection. Paradoxically, proximal sites are less affected by Heinrich events than remote sites, suggesting spatially complex event dynamics.

A cooperative PNPase-Hfq-RNA carrier complex facilitates bacterial riboregulation.

Polynucleotide phosphorylase (PNPase) is an ancient exoribonuclease conserved in the course of evolution and is found in species as diverse as bacteria and humans. Paradoxically, Escherichia coli PNPase can act not only as an RNA degrading enzyme but also by an unknown mechanism as a chaperone for small regulatory RNAs (sRNAs), with pleiotropic consequences for gene regulation. We present structures of the ternary assembly formed by PNPase, the RNA chaperone Hfq, and sRNA and show that this complex boosts sRNA stability in vitro. Comparison of structures for PNPase in RNA carrier and degradation modes reveals how the RNA is rerouted away from the active site through interactions with Hfq and the KH and S1 domains. Together, these data explain how PNPase is repurposed to protect sRNAs from cellular ribonucleases such as RNase E and could aid RNA presentation to facilitate regulatory actions on target genes.

MicroRNA-142 Is Critical for the Homeostasis and Function of Type 1 Innate Lymphoid Cells.

Natural killer (NK) cells are cytotoxic type 1 innate lymphoid cells (ILCs) that defend against viruses and mediate anti-tumor responses, yet mechanisms controlling their development and function remain incompletely understood. We hypothesized that the abundantly expressed microRNA-142 (miR-142) is a critical regulator of type 1 ILC biology. Interleukin-15 (IL-15) signaling induced miR-142 expression, whereas global and ILC-specific miR-142-deficient mice exhibited a cell-intrinsic loss of NK cells. Death of NK cells resulted from diminished IL-15 receptor signaling within miR-142-deficient mice, likely via reduced suppressor of cytokine signaling-1 (Socs1) regulation by miR-142-5p. ILCs persisting in Mir142-/- mice demonstrated increased expression of the miR-142-3p target αV integrin, which supported their survival. Global miR-142-deficient mice exhibited an expansion of ILC1-like cells concurrent with increased transforming growth factor-ß (TGF-ß) signaling. Further, miR-142-deficient mice had reduced NK-cell-dependent function and increased susceptibility to murine cytomegalovirus (MCMV) infection. Thus, miR-142 critically integrates environmental cues for proper type 1 ILC homeostasis and defense against viral infection.

GTSF1 accelerates target RNA cleavage by PIWI-clade Argonaute proteins.

Argonaute proteins use nucleic acid guides to find and bind specific DNA or RNA target sequences. Argonaute proteins have diverse biological functions and many retain their ancestral endoribonuclease activity, cleaving the phosphodiester bond between target nucleotides t10 and t11. In animals, the PIWI proteins-a specialized class of Argonaute proteins-use 21-35 nucleotide PIWI-interacting RNAs (piRNAs) to direct transposon silencing, protect the germline genome, and regulate gene expression during gametogenesis1. The piRNA pathway is required for fertility in one or both sexes of nearly all animals. Both piRNA production and function require RNA cleavage catalysed by PIWI proteins. Spermatogenesis in mice and other placental mammals requires three distinct, developmentally regulated PIWI proteins: MIWI (PIWIL1), MILI (PIWIL2) and MIWI22-4 (PIWIL4). The piRNA-guided endoribonuclease activities of MIWI and MILI are essential for the production of functional sperm5,6. piRNA-directed silencing in mice and insects also requires GTSF1, a PIWI-associated protein of unknown function7-12. Here we report that GTSF1 potentiates the weak, intrinsic, piRNA-directed RNA cleavage activities of PIWI proteins, transforming them into efficient endoribonucleases. GTSF1 is thus an example of an auxiliary protein that potentiates the catalytic activity of an Argonaute protein.

Microfluidics: reframing biological enquiry.

The underlying physical properties of microfluidic tools have led to new biological insights through the development of microsystems that can manipulate, mimic and measure biology at a resolution that has not been possible with macroscale tools. Microsystems readily handle sub-microlitre volumes, precisely route predictable laminar fluid flows and match both perturbations and measurements to the length scales and timescales of biological systems. The advent of fabrication techniques that do not require highly specialized engineering facilities is fuelling the broad dissemination of microfluidic systems and their adaptation to specific biological questions. We describe how our understanding of molecular and cell biology is being and will continue to be advanced by precision microfluidic approaches and posit that microfluidic tools - in conjunction with advanced imaging, bioinformatics and molecular biology approaches - will transform biology into a precision science.

Structural basis for piRNA targeting.

PIWI proteins use PIWI-interacting RNAs (piRNAs) to identify and silence transposable elements and thereby maintain genome integrity between metazoan generations1. The targeting of transposable elements by PIWI has been compared to mRNA target recognition by Argonaute proteins2,3, which use microRNA (miRNA) guides, but the extent to which piRNAs resemble miRNAs is not known. Here we present cryo-electron microscopy structures of a PIWI-piRNA complex from the sponge Ephydatia fluviatilis with and without target RNAs, and a biochemical analysis of target recognition. Mirroring Argonaute, PIWI identifies targets using the piRNA seed region. However, PIWI creates a much weaker seed so that stable target association requires further piRNA-target pairing, making piRNAs less promiscuous than miRNAs. Beyond the seed, the structure of PIWI facilitates piRNA-target pairing in a manner that is tolerant of mismatches, leading to long-lived PIWI-piRNA-target interactions that may accumulate on transposable-element transcripts. PIWI ensures targeting fidelity by physically blocking the propagation of piRNA-target interactions in the absence of faithful seed pairing, and by requiring an extended piRNA-target duplex to reach an endonucleolytically active conformation. PIWI proteins thereby minimize off-targeting cellular mRNAs while defending against evolving genomic threats.

Host JAK-STAT activity is a target of parasitoid wasp virulence strategies.

Innate immune responses that allow hosts to survive infection depend on the action of multiple conserved signaling pathways. Pathogens and parasites in turn have evolved virulence factors to target these immune signaling pathways in an attempt to overcome host immunity. Consequently, the interactions between host immune molecules and pathogen virulence factors play an important role in determining the outcome of an infection. The immune responses of Drosophila melanogaster provide a valuable model to understand immune signaling and host-pathogen interactions. Flies are commonly infected by parasitoid wasps and mount a coordinated cellular immune response following infection. This response is characterized by the production of specialized blood cells called lamellocytes that form a tight capsule around wasp eggs in the host hemocoel. The conserved JAK-STAT signaling pathway has been implicated in lamellocyte proliferation and is required for successful encapsulation of wasp eggs. Here we show that activity of Stat92E, the D. melanogaster STAT ortholog, is induced in immune tissues following parasitoid infection. Virulent wasp species are able to suppress Stat92E activity during infection, suggesting they target JAK-STAT pathway activation as a virulence strategy. Furthermore, two wasp species (Leptopilina guineaensis and Ganaspis xanthopoda) suppress phenotypes associated with a gain-of-function mutation in hopscotch, the D. melanogaster JAK ortholog, indicating that they inhibit the activity of the core signaling components of the JAK-STAT pathway. Our data suggest that parasitoid wasp virulence factors block JAK-STAT signaling to overcome fly immune defenses.

Specimen collection is essential for modern science.

Natural history museums are vital repositories of specimens, samples and data that inform about the natural world; this Formal Comment revisits a Perspective that advocated for the adoption of compassionate collection practices, querying whether it will ever be possible to completely do away with whole animal specimen collection.

Analyses of non-coding somatic drivers in 2,658 cancer whole genomes.

The discovery of drivers of cancer has traditionally focused on protein-coding genes1-4. Here we present analyses of driver point mutations and structural variants in non-coding regions across 2,658 genomes from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium5 of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). For point mutations, we developed a statistically rigorous strategy for combining significance levels from multiple methods of driver discovery that overcomes the limitations of individual methods. For structural variants, we present two methods of driver discovery, and identify regions that are significantly affected by recurrent breakpoints and recurrent somatic juxtapositions. Our analyses confirm previously reported drivers6,7, raise doubts about others and identify novel candidates, including point mutations in the 5' region of TP53, in the 3' untranslated regions of NFKBIZ and TOB1, focal deletions in BRD4 and rearrangements in the loci of AKR1C genes. We show that although point mutations and structural variants that drive cancer are less frequent in non-coding genes and regulatory sequences than in protein-coding genes, additional examples of these drivers will be found as more cancer genomes become available.

Assessing Target Specificity of the Small Molecule Inhibitor MARIMASTAT to Snake Venom Toxins: A Novel Application of Thermal Proteome Profiling.

New treatments that circumvent the pitfalls of traditional antivenom therapies are critical to address the problem of snakebite globally. Numerous snake venom toxin inhibitors have shown promising cross-species neutralization of medically significant venom toxins in vivo and in vitro. The development of high-throughput approaches for the screening of such inhibitors could accelerate their identification, testing, and implementation and thus holds exciting potential for improving the treatments and outcomes of snakebite envenomation worldwide. Energetics-based proteomic approaches, including thermal proteome profiling and proteome integral solubility alteration (PISA) assays, represent "deep proteomics" methods for high throughput, proteome-wide identification of drug targets and ligands. In the following study, we apply thermal proteome profiling and PISA methods to characterize the interactions between venom toxin proteoforms in Crotalus atrox (Western Diamondback Rattlesnake) and the snake venom metalloprotease (SVMP) inhibitor marimastat. We investigate its venom proteome-wide effects and characterize its interactions with specific SVMP proteoforms, as well as its potential targeting of non-SVMP venom toxin families. We also compare the performance of PISA thermal window and soluble supernatant with insoluble precipitate using two inhibitor concentrations, providing the first demonstration of the utility of a sensitive high-throughput PISA-based approach to assess the direct targets of small molecule inhibitors for snake venom.

Engineering of glycoside hydrolase family 7 cellobiohydrolases directed by natural diversity screening.

Protein engineering and screening of processive fungal cellobiohydrolases (CBHs) remain challenging due to limited expression hosts, synergy-dependency, and recalcitrant substrates. In particular, glycoside hydrolase family 7 (GH7) CBHs are critically important for the bioeconomy and typically difficult to engineer. Here, we target the discovery of highly active natural GH7 CBHs and engineering of variants with improved activity. Using experimentally assayed activities of genome mined CBHs, we applied sequence and structural alignments to top performers to identify key point mutations linked to improved activity. From ∼1500 known GH7 sequences, an evolutionarily diverse subset of 57 GH7 CBH genes was expressed in Trichoderma reesei and screened using a multiplexed activity screening assay. Ten catalytically enhanced natural variants were identified, produced, purified, and tested for efficacy using industrially relevant conditions and substrates. Three key amino acids in CBHs with performance comparable or superior to Penicillium funiculosum Cel7A were identified and combinatorially engineered into P. funiculosum cel7a, expressed in T. reesei, and assayed on lignocellulosic biomass. The top performer generated using this combined approach of natural diversity genome mining, experimental assays, and computational modeling produced a 41% increase in conversion extent over native P. funiculosum Cel7A, a 55% increase over the current industrial standard T. reesei Cel7A, and 10% improvement over Aspergillus oryzae Cel7C, the best natural GH7 CBH previously identified in our laboratory.

Genomic features of renal cell carcinoma developed during end-stage renal disease and dialysis.

Patients with end-stage renal disease (ESRD) or receiving dialysis have a much higher risk for renal cell carcinoma (RCC), but carcinogenic mechanisms and genomic features remain little explored and undefined. This study's goal was to identify the genomic features of ESRD RCC and characterize them for associations with tumor histology and dialysis exposure. In this study, we obtained 33 RCCs, with various histological subtypes, that developed in ESRD patients receiving dialysis and performed whole-genome sequencing and transcriptome analyses. Driver events, copy-number alteration (CNA) analysis and mutational signature profiling were performed using an analysis pipeline that integrated data from germline and somatic SNVs, Indels and structural variants as well as CNAs, while transcriptome data were analyzed for differentially expressed genes and through gene set enrichment analysis. ESRD related clear cell RCCs' driver genes and mutations mirrored those in sporadic ccRCCs. Longer dialysis periods significantly correlated with a rare mutational signature SBS23, whose etiology is unknown, and increased mitochondrial copy number. All acquired cystic disease (ACD)-RCCs, which developed specifically in ESRD patients, showed chromosome 16q amplification. Gene expression analysis suggests similarity between certain ACD-RCCs and papillary RCCs and in TCGA papillary RCCs with chromosome 16 gain identified enrichment for genes related to DNA repair, as well as pathways related to reactive oxygen species, oxidative phosphorylation and targets of Myc. This analysis suggests that ESRD or dialysis could induce types of cellular stress that impact some specific types of genomic damage leading to oncogenesis.

Variant spectrum of von Hippel-Lindau disease and its genomic heterogeneity in Japan.

Von Hippel-Lindau (VHL) disease is an autosomal dominant, inherited syndrome with variants in the VHL gene, causing predisposition to multi-organ neoplasms with vessel abnormality. Germline variants in VHL can be detected in 80-90% of patients clinically diagnosed with VHL disease. Here, we summarize the results of genetic tests for 206 Japanese VHL families, and elucidate the molecular mechanisms of VHL disease, especially in variant-negative unsolved cases. Of the 206 families, genetic diagnosis was positive in 175 families (85%), including 134 families (65%) diagnosed by exon sequencing (15 novel variants) and 41 (20%) diagnosed by multiplex ligation-dependent probe amplification (MLPA) (one novel variant). The deleterious variants were significantly enriched in VHL disease Type 1. Interestingly, five synonymous or non-synonymous variants within exon 2 caused exon 2 skipping, which is the first report of exon 2 skipping caused by several missense variants. Whole genome and target deep sequencing analysis were performed for 22 unsolved cases with no variant identified and found three cases with VHL mosaicism (variant allele frequency: 2.5-22%), one with mobile element insertion in the VHL promoter region, and two with a pathogenic variant of BAP1 or SDHB. The variants associated with VHL disease are heterogeneous, and for more accuracy of the genetic diagnosis of VHL disease, comprehensive genome and DNA/RNA analyses are required to detect VHL mosaicism, complicated structure variants and other related gene variants.

Candidate pathway analysis of surfactant proteins identifies CTSH and SFTA2 that influences lung cancer risk.

Pulmonary surfactant is a lipoprotein synthesized and secreted by alveolar type II cells in lung. We evaluated the associations between 200,139 single nucleotide polymorphisms (SNPs) of 40 surfactant-related genes and lung cancer risk using genotyped data from two independent lung cancer genome-wide association studies. Discovery data included 18,082 cases and 13,780 controls of European ancestry. Replication data included 1,914 cases and 3,065 controls of European descent. Using multivariate logistic regression, we found novel SNPs in surfactant-related genes CTSH [rs34577742 C > T, odds ratio (OR) = 0.90, 95% confidence interval (CI) = 0.89-0.93, P = 7.64 × 10-9] and SFTA2 (rs3095153 G > A, OR = 1.16, 95% CI = 1.10-1.21, P = 1.27 × 10-9) associated with overall lung cancer in the discovery data and validated in an independent replication data-CTSH (rs34577742 C > T, OR = 0.88, 95% CI = 0.80-0.96, P = 5.76 × 10-3) and SFTA2 (rs3095153 G > A, OR = 1.14, 95% CI = 1.01-1.28, P = 3.25 × 10-2). Among ever smokers, we found SNPs in CTSH (rs34577742 C > T, OR = 0.89, 95% CI = 0.85-0.92, P = 1.94 × 10-7) and SFTA2 (rs3095152 G > A, OR = 1.20, 95% CI = 1.14-1.27, P = 4.25 × 10-11) associated with overall lung cancer in the discovery data and validated in the replication data-CTSH (rs34577742 C > T, OR = 0.88, 95% CI = 0.79-0.97, P = 1.64 × 10-2) and SFTA2 (rs3095152 G > A, OR = 1.15, 95% CI = 1.01-1.30, P = 3.81 × 10-2). Subsequent transcriptome-wide association study using expression weights from a lung expression quantitative trait loci study revealed genes most strongly associated with lung cancer are CTSH (PTWAS = 2.44 × 10-4) and SFTA2 (PTWAS = 2.32 × 10-6).

Snake venom gene expression is coordinated by novel regulatory architecture and the integration of multiple co-opted vertebrate pathways.

Understanding how regulatory mechanisms evolve is critical for understanding the processes that give rise to novel phenotypes. Snake venom systems represent a valuable and tractable model for testing hypotheses related to the evolution of novel regulatory networks, yet the regulatory mechanisms underlying venom production remain poorly understood. Here, we use functional genomics approaches to investigate venom regulatory architecture in the prairie rattlesnake and identify cis-regulatory sequences (enhancers and promoters), trans-regulatory transcription factors, and integrated signaling cascades involved in the regulation of snake venom genes. We find evidence that two conserved vertebrate pathways, the extracellular signal-regulated kinase and unfolded protein response pathways, were co-opted to regulate snake venom. In one large venom gene family (snake venom serine proteases), this co-option was likely facilitated by the activity of transposable elements. Patterns of snake venom gene enhancer conservation, in some cases spanning 50 million yr of lineage divergence, highlight early origins and subsequent lineage-specific adaptations that have accompanied the evolution of venom regulatory architecture. We also identify features of chromatin structure involved in venom regulation, including topologically associated domains and CTCF loops that underscore the potential importance of novel chromatin structure to coevolve when duplicated genes evolve new regulatory control. Our findings provide a model for understanding how novel regulatory systems may evolve through a combination of genomic processes, including tandem duplication of genes and regulatory sequences, cis-regulatory sequence seeding by transposable elements, and diverse transcriptional regulatory proteins controlled by a co-opted regulatory cascade.

Bayesian semiparametric model for sequential treatment decisions with informative timing.

We develop a Bayesian semiparametric model for the impact of dynamic treatment rules on survival among patients diagnosed with pediatric acute myeloid leukemia (AML). The data consist of a subset of patients enrolled in a phase III clinical trial in which patients move through a sequence of four treatment courses. At each course, they undergo treatment that may or may not include anthracyclines (ACT). While ACT is known to be effective at treating AML, it is also cardiotoxic and can lead to early death for some patients. Our task is to estimate the potential survival probability under hypothetical dynamic ACT treatment strategies, but there are several impediments. First, since ACT is not randomized, its effect on survival is confounded over time. Second, subjects initiate the next course depending on when they recover from the previous course, making timing potentially informative of subsequent treatment and survival. Third, patients may die or drop out before ever completing the full treatment sequence. We develop a generative Bayesian semiparametric model based on Gamma Process priors to address these complexities. At each treatment course, the model captures subjects' transition to subsequent treatment or death in continuous time. G-computation is used to compute a posterior over potential survival probability that is adjusted for time-varying confounding. Using our approach, we estimate the efficacy of hypothetical treatment rules that dynamically modify ACT based on evolving cardiac function.