A yeast living ancestor reveals the origin of genomic introgressions.

Genome introgressions drive evolution across the animal1, plant2 and fungal3 kingdoms. Introgressions initiate from archaic admixtures followed by repeated backcrossing to one parental species. However, how introgressions arise in reproductively isolated species, such as yeast4, has remained unclear. Here we identify a clonal descendant of the ancestral yeast hybrid that founded the extant Saccharomyces cerevisiae Alpechin lineage5, which carries abundant Saccharomyces paradoxus introgressions. We show that this clonal descendant, hereafter defined as a 'living ancestor', retained the ancestral genome structure of the first-generation hybrid with contiguous S. cerevisiae and S. paradoxus subgenomes. The ancestral first-generation hybrid underwent catastrophic genomic instability through more than a hundred mitotic recombination events, mainly manifesting as homozygous genome blocks generated by loss of heterozygosity. These homozygous sequence blocks rescue hybrid fertility by restoring meiotic recombination and are the direct origins of the introgressions present in the Alpechin lineage. We suggest a plausible route for introgression evolution through the reconstruction of extinct stages and propose that genome instability allows hybrids to overcome reproductive isolation and enables introgressions to emerge.

ScRAPdb: an integrated pan-omics database for the Saccharomyces cerevisiae reference assembly panel.

As a unicellular eukaryote, the budding yeast Saccharomyces cerevisiae strikes a unique balance between biological complexity and experimental tractability, serving as a long-standing classic model for both basic and applied studies. Recently, S. cerevisiae further emerged as a leading system for studying natural diversity of genome evolution and its associated functional implication at population scales. Having high-quality comparative and functional genomics data are critical for such efforts. Here, we exhaustively expanded the telomere-to-telomere (T2T) S. cerevisiae reference assembly panel (ScRAP) that we previously constructed for 142 strains to cover high-quality genome assemblies and annotations of 264 S. cerevisiae strains from diverse geographical and ecological niches and also 33 outgroup strains from all the other Saccharomyces species complex. We created a dedicated online database, ScRAPdb (https://www.evomicslab.org/db/ScRAPdb/), to host this expanded pangenome collection. Furthermore, ScRAPdb also integrates an array of population-scale pan-omics atlases (pantranscriptome, panproteome and panphenome) and extensive data exploration toolkits for intuitive genomics analyses. All curated data and downstream analysis results can be easily downloaded from ScRAPdb. We expect ScRAPdb to become a highly valuable platform for the yeast community and beyond, leading to a pan-omics understanding of the global genetic and phenotypic diversity.

Genome instability footprint under rapamycin and hydroxyurea treatments.

The mutational processes dictating the accumulation of mutations in genomes are shaped by genetic background, environment and their interactions. Accurate quantification of mutation rates and spectra under drugs has important implications in disease treatment. Here, we used whole-genome sequencing and time-resolved growth phenotyping of yeast mutation accumulation lines to give a detailed view of the mutagenic effects of rapamycin and hydroxyurea on the genome and cell growth. Mutation rates depended on the genetic backgrounds but were only marginally affected by rapamycin. As a remarkable exception, rapamycin treatment was associated with frequent chromosome XII amplifications, which compensated for rapamycin induced rDNA repeat contraction on this chromosome and served to maintain rDNA content homeostasis and fitness. In hydroxyurea, a wide range of mutation rates were elevated regardless of the genetic backgrounds, with a particularly high occurrence of aneuploidy that associated with dramatic fitness loss. Hydroxyurea also induced a high T-to-G and low C-to-A transversion rate that reversed the common G/C-to-A/T bias in yeast and gave rise to a broad range of structural variants, including mtDNA deletions. The hydroxyurea mutation footprint was consistent with the activation of error-prone DNA polymerase activities and non-homologues end joining repair pathways. Taken together, our study provides an in-depth view of mutation rates and signatures in rapamycin and hydroxyurea and their impact on cell fitness, which brings insights for assessing their chronic effects on genome integrity.

RecombineX: A generalized computational framework for automatic high-throughput gamete genotyping and tetrad-based recombination analysis.

Meiotic recombination is an essential biological process that ensures faithful chromosome segregation and promotes parental allele shuffling. Tetrad analysis is a powerful approach to quantify the genetic makeups and recombination landscapes of meiotic products. Here we present RecombineX (https://github.com/yjx1217/RecombineX), a generalized computational framework that automates the full workflow of marker identification, gamete genotyping, and tetrad-based recombination profiling based on any organism or genetic background with batch processing capability. Aside from conventional reference-based analysis, RecombineX can also perform analysis based on parental genome assemblies, which facilitates analyzing meiotic recombination landscapes in their native genomic contexts. Additional features such as copy number variation profiling and missing genotype inference further enhance downstream analysis. RecombineX also includes a dedicate module for simulating the genomes and reads of recombinant tetrads, which enables fine-tuned simulation-based hypothesis testing. This simulation module revealed the power and accuracy of RecombineX even when analyzing tetrads with very low sequencing depths (e.g., 1-2X). Tetrad sequencing data from the budding yeast Saccharomyces cerevisiae and green alga Chlamydomonas reinhardtii were further used to demonstrate the accuracy and robustness of RecombineX for organisms with both small and large genomes, manifesting RecombineX as an all-around one stop solution for future tetrad analysis. Interestingly, our re-analysis of the budding yeast tetrad sequencing data with RecombineX and Oxford Nanopore sequencing revealed two unusual structural rearrangement events that were not noticed before, which exemplify the occasional genome instability triggered by meiosis.

GetPrimers: A generalized PCR-based genetic targeting primer designer enabling easy and standardized targeted gene modification across multiple systems.

Genetic targeting (e.g., gene knockout and tagging) based on polymerase chain reaction (PCR) is a simple yet powerful approach for studying gene functions. Although originally developed in classic budding and fission yeast models, the same principle applies to other eukaryotic systems with efficient homologous recombination. One-step PCR-based genetic targeting is conventionally used but the sizes of the homologous arms that it generates for recombination-mediated genetic targeting are usually limited. Alternatively, gene targeting can also be performed via fusion PCR, which can create homologous arms that are orders of magnitude larger, therefore substantially increasing the efficiency of recombination-mediated genetic targeting. Here, we present GetPrimers (https://www.evomicslab.org/app/getprimers/), a generalized computational framework and web tool to assist automatic targeting and verification primer design for both one-step PCR-based and fusion PCR-based genetic targeting experiments. Moreover, GetPrimers by design runs for any given genetic background of any species with full genome scalability. Therefore, GetPrimers is capable of empowering high-throughput functional genomic assays at multipopulation and multispecies levels. Comprehensive experimental validations have been performed for targeting and verification primers designed by GetPrimers across multiple organism systems and experimental setups. We anticipate GetPrimers to become a highly useful and popular tool to facilitate easy and standardized gene modification across multiple systems.

Vasorin (VASN) overexpression promotes pulmonary metastasis and resistance to adjuvant chemotherapy in patients with locally advanced rectal cancer.

BACKGROUND: LARC patients commonly receive adjuvant therapy, however, hidden micrometastases still limit the improvement of OS. This study aims to investigate the impact of VASN in rectal cancer with pulmonary metastasis and understand the underlying molecular mechanisms to guide adjuvant chemotherapy selection. METHODS: Sequencing data from rectal cancer patients with pulmonary metastasis from Sun Yat-sen University Cancer Center (SYSUCC) and publicly available data were meticulously analyzed. The functional role of VASN in pulmonary metastasis was validated in vivo and in vitro. Coimmunoprecipitation (co-IP), immunofluorescence, and rescue experiments were conducted to unravel potential molecular mechanisms of VASN. Moreover, VASN expression levels in tumor samples were examined and analyzed for their correlations with pulmonary metastasis status, tumor stage, adjuvant chemotherapy benefit, and survival outcome. RESULTS: Our study revealed a significant association between high VASN expression and pulmonary metastasis in LARC patients. Experiments in vitro and in vivo demonstrated that VASN could promote the cell proliferation, metastasis, and drug resistance of colorectal cancer. Mechanistically, VASN interacts with the NOTCH1 protein, leading to concurrent activation of the NOTCH and MAPK pathways. Clinically, pulmonary metastasis and advanced tumor stage were observed in 90% of VASN-positive patients and 53.5% of VASN-high patients, respectively, and VASN-high patients had a lower five-year survival rate than VASN-low patients (26.7% vs. 83.7%). Moreover, the Cox analysis and OS analysis indicated that VASN was an independent prognostic factor for OS (HR = 7.4, P value < 0.001) and a predictor of adjuvant therapy efficacy in rectal cancer. CONCLUSIONS: Our study highlights the role of VASN in decreasing drug sensitivity and activating the NOTCH and MAPK pathways, which leads to tumorigenesis and pulmonary metastasis. Both experimental and clinical data support that rectal cancer patients with VASN overexpression detected in biopsies have a higher risk of pulmonary metastasis and adjuvant chemotherapy resistance.

Genome evolution across 1,011 Saccharomyces cerevisiae isolates.

Large-scale population genomic surveys are essential to explore the phenotypic diversity of natural populations. Here we report the whole-genome sequencing and phenotyping of 1,011 Saccharomyces cerevisiae isolates, which together provide an accurate evolutionary picture of the genomic variants that shape the species-wide phenotypic landscape of this yeast. Genomic analyses support a single 'out-of-China' origin for this species, followed by several independent domestication events. Although domesticated isolates exhibit high variation in ploidy, aneuploidy and genome content, genome evolution in wild isolates is mainly driven by the accumulation of single nucleotide polymorphisms. A common feature is the extensive loss of heterozygosity, which represents an essential source of inter-individual variation in this mainly asexual species. Most of the single nucleotide polymorphisms, including experimentally identified functional polymorphisms, are present at very low frequencies. The largest numbers of variants identified by genome-wide association are copy-number changes, which have a greater phenotypic effect than do single nucleotide polymorphisms. This resource will guide future population genomics and genotype-phenotype studies in this classic model system.

Intragenic repeat expansion in the cell wall protein gene HPF1 controls yeast chronological aging.

Aging varies among individuals due to both genetics and environment, but the underlying molecular mechanisms remain largely unknown. Using a highly recombined Saccharomyces cerevisiae population, we found 30 distinct quantitative trait loci (QTLs) that control chronological life span (CLS) in calorie-rich and calorie-restricted environments and under rapamycin exposure. Calorie restriction and rapamycin extended life span in virtually all genotypes but through different genetic variants. We tracked the two major QTLs to the cell wall glycoprotein genes FLO11 and HPF1 We found that massive expansion of intragenic tandem repeats within the N-terminal domain of HPF1 was sufficient to cause pronounced life span shortening. Life span impairment by HPF1 was buffered by rapamycin but not by calorie restriction. The HPF1 repeat expansion shifted yeast cells from a sedentary to a buoyant state, thereby increasing their exposure to surrounding oxygen. The higher oxygenation altered methionine, lipid, and purine metabolism, and inhibited quiescence, which explains the life span shortening. We conclude that fast-evolving intragenic repeat expansions can fundamentally change the relationship between cells and their environment with profound effects on cellular lifestyle and longevity.

Molecular asymmetry in the cephalochordate embryo revealed by single-blastomere transcriptome profiling.

Studies in various animals have shown that asymmetrically localized maternal transcripts play important roles in axial patterning and cell fate specification in early embryos. However, comprehensive analyses of the maternal transcriptomes with spatial information are scarce and limited to a handful of model organisms. In cephalochordates (amphioxus), an early branching chordate group, maternal transcripts of germline determinants form a compact granule that is inherited by a single blastomere during cleavage stages. Further blastomere separation experiments suggest that other transcripts associated with the granule are likely responsible for organizing the posterior structure in amphioxus; however, the identities of these determinants remain unknown. In this study, we used high-throughput RNA sequencing of separated blastomeres to examine asymmetrically localized transcripts in two-cell and eight-cell stage embryos of the amphioxus Branchiostoma floridae. We identified 111 and 391 differentially enriched transcripts at the 2-cell stage and the 8-cell stage, respectively, and used in situ hybridization to validate the spatial distribution patterns for a subset of these transcripts. The identified transcripts could be categorized into two major groups: (1) vegetal tier/germ granule-enriched and (2) animal tier/anterior-enriched transcripts. Using zebrafish as a surrogate model system, we showed that overexpression of one animal tier/anterior-localized amphioxus transcript, zfp665, causes a dorsalization/anteriorization phenotype in zebrafish embryos by downregulating the expression of the ventral gene, eve1, suggesting a potential function of zfp665 in early axial patterning. Our results provide a global transcriptomic blueprint for early-stage amphioxus embryos. This dataset represents a rich platform to guide future characterization of molecular players in early amphioxus development and to elucidate conservation and divergence of developmental programs during chordate evolution.

Natural variants suppress mutations in hundreds of essential genes.

The consequence of a mutation can be influenced by the context in which it operates. For example, loss of gene function may be tolerated in one genetic background, and lethal in another. The extent to which mutant phenotypes are malleable, the architecture of modifiers and the identities of causal genes remain largely unknown. Here, we measure the fitness effects of ~ 1,100 temperature-sensitive alleles of yeast essential genes in the context of variation from ten different natural genetic backgrounds and map the modifiers for 19 combinations. Altogether, fitness defects for 149 of the 580 tested genes (26%) could be suppressed by genetic variation in at least one yeast strain. Suppression was generally driven by gain-of-function of a single, strong modifier gene, and involved both genes encoding complex or pathway partners suppressing specific temperature-sensitive alleles, as well as general modifiers altering the effect of many alleles. The emerging frequency of suppression and range of possible mechanisms suggest that a substantial fraction of monogenic diseases could be managed by modulating other gene products.

Human RAP1 specifically protects telomeres of senescent cells from DNA damage.

Repressor/activator protein 1 (RAP1) is a highly evolutionarily conserved protein found at telomeres. Although yeast Rap1 is a key telomere capping protein preventing non-homologous end joining (NHEJ) and consequently telomere fusions, its role at mammalian telomeres in vivo is still controversial. Here, we demonstrate that RAP1 is required to protect telomeres in replicative senescent human cells. Downregulation of RAP1 in these cells, but not in young or dividing pre-senescent cells, leads to telomere uncapping and fusions. The anti-fusion effect of RAP1 was further explored in a HeLa cell line where RAP1 expression was depleted through an inducible CRISPR/Cas9 strategy. Depletion of RAP1 in these cells gives rise to telomere fusions only when telomerase is inhibited. We further show that the fusions triggered by RAP1 loss are dependent upon DNA ligase IV. We conclude that human RAP1 is specifically involved in protecting critically short telomeres. This has important implications for the functions of telomeres in senescent cells.

Hemichordate genomes and deuterostome origins.

Acorn worms, also known as enteropneust (literally, 'gut-breathing') hemichordates, are marine invertebrates that share features with echinoderms and chordates. Together, these three phyla comprise the deuterostomes. Here we report the draft genome sequences of two acorn worms, Saccoglossus kowalevskii and Ptychodera flava. By comparing them with diverse bilaterian genomes, we identify shared traits that were probably inherited from the last common deuterostome ancestor, and then explore evolutionary trajectories leading from this ancestor to hemichordates, echinoderms and chordates. The hemichordate genomes exhibit extensive conserved synteny with amphioxus and other bilaterians, and deeply conserved non-coding sequences that are candidates for conserved gene-regulatory elements. Notably, hemichordates possess a deuterostome-specific genomic cluster of four ordered transcription factor genes, the expression of which is associated with the development of pharyngeal 'gill' slits, the foremost morphological innovation of early deuterostomes, and is probably central to their filter-feeding lifestyle. Comparative analysis reveals numerous deuterostome-specific gene novelties, including genes found in deuterostomes and marine microbes, but not other animals. The putative functions of these genes can be linked to physiological, metabolic and developmental specializations of the filter-feeding ancestor.

Shared Molecular Targets Confer Resistance over Short and Long Evolutionary Timescales.

Pre-existing and de novo genetic variants can both drive adaptation to environmental changes, but their relative contributions and interplay remain poorly understood. Here we investigated the evolutionary dynamics in drug-treated yeast populations with different levels of pre-existing variation by experimental evolution coupled with time-resolved sequencing and phenotyping. We found a doubling of pre-existing variation alone boosts the adaptation by 64.1% and 51.5% in hydroxyurea and rapamycin, respectively. The causative pre-existing and de novo variants were selected on shared targets: RNR4 in hydroxyurea and TOR1, TOR2 in rapamycin. Interestingly, the pre-existing and de novo TOR variants map to different functional domains and act via distinct mechanisms. The pre-existing TOR variants from two domesticated strains exhibited opposite rapamycin resistance effects, reflecting lineage-specific functional divergence. This study provides a dynamic view on how pre-existing and de novo variants interactively drive adaptation and deepens our understanding of clonally evolving populations.

simuG: a general-purpose genome simulator.

SUMMARY: Simulated genomes with pre-defined and random genomic variants can be very useful for benchmarking genomic and bioinformatics analyses. Here we introduce simuG, a lightweight tool for simulating the full-spectrum of genomic variants (single nucleotide polymorphisms, Insertions/Deletions, copy number variants, inversions and translocations) for any organisms (including human). The simplicity and versatility of simuG make it a unique general-purpose genome simulator for a wide-range of simulation-based applications. AVAILABILITY AND IMPLEMENTATION: Code in Perl along with user manual and testing data is available at https://github.com/yjx1217/simuG. This software is free for use under the MIT license. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

A critical time window for the analgesic effect of central histamine in the partial sciatic ligation model of neuropathic pain.

BACKGROUND: It is known that histamine participates in pain modulation. However, the effect of central histamine on neuropathic pain is not fully understood. Here, we report a critical time window for the analgesic effect of central histamine in the partial sciatic nerve ligation model of neuropathic pain. METHODS: Neuropathic pain was induced by partial sciatic nerve ligation (PSL) in rats, wild-type (C57BL/6J) mice and HDC(-/-) (histidine decarboxylase gene knockout) and IL-1R(-/-) (interleukin-1 receptor gene knockout) mice. Histidine, a precursor of histamine that can increase the central histamine levels, was administered intraperitoneally (i.p.). Histidine decarboxylase (HDC) enzyme inhibitor α-fluoromethylhistidine was administered intracerebroventricularly (i.c.v.). Histamine H1 receptor antagonist mepyramine and H2 receptor antagonist cimetidine were given intrathecally (i.t.) and intracisternally (i.c.). Withdrawal thresholds to tactile and heat stimuli were measured with a set of von Frey hairs and infrared laser, respectively. Immunohistochemistry and Western blot were carried out to evaluate the morphology of microglia and IL-1ß production, respectively. RESULTS: Histidine (100 mg/kg, i.p.) administered throughout days 0-3, 0-7, or 0-14 postoperatively (PO) alleviated mechanical allodynia and thermal hyperalgesia in the hindpaw following PSL in rats. Intrathecal histamine reversed PSL-induced thermal hyperalgesia in a dose-dependent manner and intracisternal histamine alleviated both mechanical allodynia and thermal hyperalgesia. Moreover, α-fluoromethylhistidine (i.c.v.) abrogated the analgesic effect of histidine. However, histidine treatment initiated later than the first postoperative day (treatment periods included days 2-3, 4-7, and 8-14 PO) did not show an analgesic effect. In addition, histidine treatment initiated immediately, but not 3 days after PSL, inhibited microglial activation and IL-1ß upregulation in the lumbar spinal cord, in parallel with its effects on behavioral hypersensitivity. Moreover, the inhibitory effects on pain hypersensitivity and spinal microglial activation were absent in HDC(-/-) mice and IL-1R(-/-) mice. H1 receptor antagonist mepyramine (200 ng/rat i.t. or i.c.), but not H2 receptor antagonist cimetidine (200, 500 ng/rat i.t. or 500 ng/rat i.c.), blocked the effects of histidine on pain behavior and spinal microglia. CONCLUSIONS: These results demonstrate that central histamine is analgesic within a critical time window in the PSL model of neuropathic pain via histamine H1 receptors. This effect may partly relate to the inhibition of microglial activation and IL-1ß production in the spinal cord following nerve injury.

Ploidy influences the functional attributes of de novo lager yeast hybrids.

The genomes of hybrid organisms, such as lager yeast (Saccharomyces cerevisiae × Saccharomyces eubayanus), contain orthologous genes, the functionality and effect of which may differ depending on their origin and copy number. How the parental subgenomes in lager yeast contribute to important phenotypic traits such as fermentation performance, aroma production, and stress tolerance remains poorly understood. Here, three de novo lager yeast hybrids with different ploidy levels (allodiploid, allotriploid, and allotetraploid) were generated through hybridization techniques without genetic modification. The hybrids were characterized in fermentations of both high gravity wort (15 °P) and very high gravity wort (25 °P), which were monitored for aroma compound and sugar concentrations. The hybrid strains with higher DNA content performed better during fermentation and produced higher concentrations of flavor-active esters in both worts. The hybrid strains also outperformed both the parent strains. Genome sequencing revealed that several genes related to the formation of flavor-active esters (ATF1, ATF2¸ EHT1, EEB1, and BAT1) were present in higher copy numbers in the higher ploidy hybrid strains. A direct relationship between gene copy number and transcript level was also observed. The measured ester concentrations and transcript levels also suggest that the functionality of the S. cerevisiae- and S. eubayanus-derived gene products differs. The results contribute to our understanding of the complex molecular mechanisms that determine phenotypes in lager yeast hybrids and are expected to facilitate targeted strain development through interspecific hybridization.

ERK1 and ERK2 present functional redundancy in tetrapods despite higher evolution rate of ERK1.

BACKGROUND: The Ras/Raf/MEK/ERK signaling pathway is involved in essential cell processes and it is abnormally activated in ~30 % of cancers and cognitive disorders. Two ERK isoforms have been described, ERK1 and ERK2; ERK2 being regarded by many as essential due to the embryonic lethality of ERK2 knock-out mice, whereas mice lacking ERK1 are viable and fertile. The controversial question of why we have two ERKs and whether they have differential functions or display functional redundancy has not yet been resolved. RESULTS: To investigate this question we used a novel approach based on comparing the evolution of ERK isoforms' sequences and protein expression across vertebrates. We gathered and cloned erk1 and erk2 coding sequences and we examined protein expression of isoforms in brain extracts in all major clades of vertebrate evolution. For the first time, we measured each isoforms' relative protein level in phylogenetically distant animals using anti-phospho antibodies targeting active ERKs. We demonstrate that squamates (lizards, snakes and geckos), despite having both genes, do not express ERK2 protein whereas other tetrapods either do not express ERK1 protein or have lost the erk1 gene. To demonstrate the unexpected squamates' lack of ERK2 expression, we targeted each ERK isoform in lizard primary fibroblasts by specific siRNA-mediated knockdown. We also found that undetectable expression of ERK2 in lizard is compensated by a greater strength of lizard's erk1 promoter. Finally, phylogenetic analysis revealed that ERK1 amino acids sequences evolve faster than ERK2's likely due to genomic factors, including a large difference in gene size, rather than from functional differences since amino acids essential for function are kept invariant. CONCLUSIONS: ERK isoforms appeared by a single gene duplication at the onset of vertebrate evolution at least 400 Mya. Our results demonstrate that tetrapods can live by expressing either one or both ERK isoforms, supporting the notion that ERK1/2 act interchangeably. Substrate recognition sites and catalytic cleft are nearly invariant in all vertebrate ERKs further suggesting functional redundancy. We suggest that future ERK research should shift towards understanding the role and regulation of total ERK quantity, especially in light of newly described erk2 gene amplification identified in tumors.

NanoTrans: an integrated computational framework for comprehensive transcriptome analysis with nanopore direct RNA sequencing.

Nanopore direct RNA sequencing (DRS) provides the direct access to native RNA strands with full-length information, shedding light on rich qualitative and quantitative properties of gene expression profiles. Here with NanoTrans, we present an integrated computational framework that comprehensively covers all major DRS-based application scopes, including isoform clustering and quantification, poly(A) tail length estimation, RNA modification profiling, and fusion gene detection. In addition to its merit in providing such a streamlined one-stop solution, NanoTrans also shines in its workflow-orientated modular design, batch processing capability, all-in-one tabular and graphic report output, as well as automatic installation and configuration supports. Finally, by applying NanoTrans to real DRS datasets of yeast, Arabidopsis, as well as human embryonic kidney and cancer cell lines, we further demonstrate its utility, effectiveness, and efficacy across a wide range of DRS-based application settings.

Full-spectral genome analysis of natural killer/T cell lymphoma highlights impacts of genome instability in driving its progression.

BACKGROUND: Natural killer/T cell lymphoma (NKTCL) is a clinically and genetically heterogeneous disease with poor prognosis. Genome sequencing and mutation characterization provides a powerful approach for patient stratification, treatment target discovery, and etiology identification. However, previous studies mostly concentrated on base-level mutations in primary NKTCL, whereas the large-scale genomic alterations in NKTCL and the mutational landscapes in relapsed/refractory NKTCL remain largely unexplored. METHODS: Here, we assembled whole-genome sequencing and whole-exome sequencing data from 163 patients with primary or relapsed/refractory NKTCL and compared their somatic mutational landscapes at both nucleotide and structure levels. RESULTS: Our study not only confirmed previously reported common NKTCL mutational targets like STAT3, TP53, and DDX3X but also unveiled several novel high-frequency mutational targets such as PRDM9, DST, and RBMX. In terms of the overall mutational landscape, we observed striking differences between primary and relapsed/refractory NKTCL patient groups, with the latter exhibits higher levels of tumor mutation burden, copy number variants (CNVs), and structural variants (SVs), indicating a strong signal of genomic instability. Complex structural rearrangements such as chromothripsis and focal amplification are also significantly enriched in relapsed/refractory NKTCL patients, exerting a substantial impact on prognosis. Accordingly, we devised a novel molecular subtyping system (i.e., C0-C4) with distinct prognosis by integrating potential driver mutations at both nucleotide and structural levels, which further provides an informative guidance for novel treatments that target these specific driver mutations and genome instability as a whole. CONCLUSIONS: The striking differences underlying the mutational landscapes between the primary and relapsed/refractory NKTCL patients highlight the importance of genomic instability in driving the progression of NKTCL. Our newly proposed molecular subtyping system is valuable in assisting patient stratification and novel treatment design towards a better prognosis in the age of precision medicine.

Alternations in inflammatory macrophage niche drive phenotypic and functional plasticity of Kupffer cells.

Inflammatory signals lead to recruitment of circulating monocytes and induce their differentiation into pro-inflammatory macrophages. Therefore, whether blocking inflammatory monocytes can mitigate disease progression is being actively evaluated. Here, we employ multiple lineage-tracing models and show that monocyte-derived macrophages (mo-mac) are the major population of immunosuppressive, liver metastasis-associated macrophages (LMAM), while the proportion of Kupffer cells (KC) as liver-resident macrophages is diminished in metastatic nodules. Paradoxically, genetic ablation of mo-macs results in only a marginal decrease in LMAMs. Using a proliferation-recording system and a KC-tracing model in a monocyte-deficient background, we find that LMAMs can be replenished either via increased local macrophage proliferation or by promoting KC infiltration. In the latter regard, KCs undergo transient proliferation and exhibit substantial phenotypic and functional alterations through epigenetic reprogramming following the vacating of macrophage niches by monocyte depletion. Our data thus suggest that a simultaneous blockade of monocyte recruitment and macrophage proliferation may effectively target immunosuppressive myelopoiesis and reprogram the microenvironment towards an immunostimulatory state.