Goat domestication was critical for agriculture and civilization, but its underlying genetic changes and selection regimes remain unclear. Here, we analyze the genomes of worldwide domestic goats, wild caprid species, and historical remains, providing evidence of an ancient introgression event from a West Caucasian tur-like species to the ancestor of domestic goats. One introgressed locus with a strong signature of selection harbors the MUC6 gene, which encodes a gastrointestinally secreted mucin. Experiments revealed that the nearly fixed introgressed haplotype confers enhanced immune resistance to gastrointestinal pathogens. Another locus with a strong signal of selection may be related to behavior. The selected alleles at these two loci emerged in domestic goats at least 7200 and 8100 years ago, respectively, and increased to high frequencies concurrent with the expansion of the ubiquitous modern mitochondrial haplogroup A. Tracking these archaeologically cryptic evolutionary transformations provides new insights into the mechanisms of animal domestication.

Mesenchymal stem cell (MSC)-mediated immunomodulation has been harnessed for the treatment of human diseases, but its underlying mechanism has not been fully understood. Dead cells, including apoptotic cells have immunomodulatory properties. It has been repeatedly reported that the proportion of nonviable MSCs in a MSC therapeutic preparation varied from 5~50% in the ongoing clinical trials. It is conceivable that the nonviable cells in a MSC therapeutic preparation may play a role in the therapeutic effects of MSCs. We found that the MSC therapeutic preparation in the present study had about 5% dead MSCs (DMSCs), characterized by apoptotic cells. Namely, 1 × 106 MSCs in the preparation contained about 5 × 104 DMSCs. We found that the treatment with even 5 × 104 DMSCs alone had the equal therapeutic effects as with 1 × 106 MSCs. This protective effect of the dead MSCs alone was confirmed in four mouse models, including concanavalin A (ConA)- and carbon tetrachloride (CCl4)-induced acute liver injury, LPS-induced lung injury and spinal cord injury. We also found that the infused MSCs died by apoptosis in vivo. Furthermore, the therapeutic effect was attributed to the elevated level of phosphatidylserine (PS) upon the injection of MSCs or DMSCs. The direct administration of PS liposomes (PSLs) mimic apoptotic cell fragments also exerted the protective effects as MSCs and DMSCs. The Mer tyrosine kinase (MerTK) deficiency or the knockout of chemokine receptor C-C motif chemokine receptor 2 (CCR2) reversed these protective effects of MSCs or DMSCs. These results revealed that DMSCs alone in the therapeutic stem cell preparation or the apoptotic cells induced in vivo may exert the same immunomodulatory property as the "living MSCs preparation" through releasing PS, which was further recognized by MerTK and participated in modulating immune cells.

Tumor vasculature often exhibits disorder and inefficiency. Vascular normalization offers potential for alleviating hypoxia and optimizing drug delivery in tumors. However, identifying effective agents is hindered by a lack of robust screening. We aimed to establish a comprehensive method using the zebrafish functional xenograft vasculature platform (zFXVP) to visualize and quantify tumor vasculature changes. Employing zFXVP, we systematically screened compounds, identifying PF-502 as a robust vascular normalization agent. Mechanistic studies showed PF-502 induces endothelial cell-cycle arrest, streamlines vasculature, and activates Notch1 signaling, enhancing stability and hemodynamics. In murine models, PF-502 exhibited pronounced vascular normalization and improved drug delivery at a sub-maximum tolerated dose. These findings highlight zFXVP's utility and suggest PF-502 as a promising adjunctive for vascular normalization in clinical settings.

Leukemogenesis is proposed to be a multistep process by which normal hematopoietic stem and progenitor cells are transformed into full-blown leukemic cells, the details of which are not fully understood. Here, we performed serial single-cell transcriptome analyses of preleukemic and leukemic cells (PLCs) and constructed the cellular and molecular transformation trajectory in a Myc-driven acute myeloid leukemia (AML) model in mice, which represented the transformation course in patients. We found that the Myc targets were gradually up-regulated along the trajectory. Among them were splicing factors, which showed stage-specific prognosis for AML patients. Furthermore, we dissected the detailed gene network of a tipping point for hematopoietic stem and progenitor cells (HSPCs) to generate initiating PLCs, which was characterized by dramatically increased splicing factors and unusual RNA velocity. In the late stage, PLCs acquired explosive heterogeneity through RNA alternative splicing. Among them, the Hsp90aa1hi subpopulation was conserved in both human and mouse AML and associated with poor prognosis. Exon 4 skipping of Tmem134 was identified in these cells. While the exon skipping product Tmem134β promoted the cell cycle, full-length Tmem134α delayed tumorigenesis. Our study emphasized the critical roles of RNA splicing in the full process of leukemogenesis.

Axolotl (Ambystoma mexicanum) is an excellent model for investigating regeneration, the interaction between regenerative and developmental processes, comparative genomics, and evolution. The brain, which serves as the material basis of consciousness, learning, memory, and behavior, is the most complex and advanced organ in axolotl. The modulation of transcription factors is a crucial aspect in determining the function of diverse regions within the brain. There is, however, no comprehensive understanding of the gene regulatory network of axolotl brain regions. Here, we utilized single-cell ATAC sequencing to generate the chromatin accessibility landscapes of 81,199 cells from the olfactory bulb, telencephalon, diencephalon and mesencephalon, hypothalamus and pituitary, and the rhombencephalon. Based on these data, we identified key transcription factors specific to distinct cell types and compared cell type functions across brain regions. Our results provide a foundation for comprehensive analysis of gene regulatory programs, which are valuable for future studies of axolotl brain development, regeneration, and evolution, as well as on the mechanisms underlying cell-type diversity in vertebrate brains.

Early consumption of starter feed promotes rumen development in lambs. We examined rumen development in lambs fed starter feed for 5 weeks using histological and biochemical analyses and by performing high-throughput sequencing in rumen tissues. Additionally, rumen contents of starter feed-fed lambs were compared to those of breast milk-fed controls. Our physiological and biochemical findings revealed that early starter consumption facilitated rumen development, changed the pattern of ruminal fermentation, and increased the amylase and carboxymethylcellulase activities of rumen micro-organisms. RNA-seq analysis revealed 225 differentially expressed genes between the rumens of breast milk- and starter feed-fed lambs. These DEGs were involved in many metabolic pathways, particularly lipid and carbohydrate metabolism, and included HMGCL and HMGCS2. Sequencing analysis of 16S rRNA genes revealed that ruminal bacterial communities were more diverse in breast milk-than in starter feed-fed lambs, and each group had a distinct microbiota. We conclude that early starter feeding is beneficial to rumen development and physiological function in lambs. The underlying mechanism may involve the stimulation of ruminal ketogenesis and butanoate metabolism via HMGCL and HMGCS2 combined with changes in the fermentation type induced by ruminal microbiota. Overall, this study provides insights into the molecular mechanisms of rumen development in sheep.

Chromosome copy-number variations are a hallmark of cancer. Among them, the prevalent chromosome 17p deletions are associated with poor prognosis and can promote tumorigenesis more than TP53 loss. Here, we use multiple functional genetic strategies and identify a new 17p tumor suppressor gene (TSG), plant homeodomain finger protein 23 (PHF23). Its deficiency impairs B-cell differentiation and promotes immature B-lymphoblastic malignancy. Mechanistically, we demonstrate that PHF23, an H3K4me3 reader, directly binds the SIN3-HDAC complex through its N-terminus and represses its deacetylation activity on H3K27ac. Thus, the PHF23-SIN3-HDAC (PSH) complex coordinates these two major active histone markers for the activation of downstream TSGs and differentiation-related genes. Furthermore, dysregulation of the PSH complex is essential for the development and maintenance of PHF23-deficient and 17p-deleted tumors. Hence, our study reveals a novel epigenetic regulatory mechanism that contributes to the pathology of 17p-deleted cancers and suggests a susceptibility in this disease. SIGNIFICANCE: We identify PHF23, encoding an H3K4me3 reader, as a new TSG on chromosome 17p, which is frequently deleted in human cancers. Mechanistically, PHF23 forms a previously unreported histone-modifying complex, the PSH complex, which regulates gene activation through a synergistic link between H3K4me3 and H3K27ac.This article is highlighted in the In This Issue feature, p. 1.

Small cell lung cancer (SCLC) is notorious for its early and frequent metastases, which contribute to it as a recalcitrant malignancy. To understand the molecular mechanisms underlying SCLC metastasis, we generated SCLC mouse models with orthotopically transplanted genome-edited lung organoids and performed multiomics analyses. We found that a deficiency of KMT2C, a histone H3 lysine 4 methyltransferase frequently mutated in extensive-stage SCLC, promoted multiple-organ metastases in mice. Metastatic and KMT2C-deficient SCLC displayed both histone and DNA hypomethylation. Mechanistically, KMT2C directly regulated the expression of DNMT3A, a de novo DNA methyltransferase, through histone methylation. Forced DNMT3A expression restrained metastasis of KMT2C-deficient SCLC through repressing metastasis-promoting MEIS/HOX genes. Further, S-(5'-adenosyl)-L-methionine, the common cofactor of histone and DNA methyltransferases, inhibited SCLC metastasis. Thus, our study revealed a concerted epigenetic reprogramming of KMT2C- and DNMT3A-mediated histone and DNA hypomethylation underlying SCLC metastasis, which suggested a potential epigenetic therapeutic vulnerability.

Endometrial cancer (EC) is the most common female reproductive tract cancer and its incidence has been continuously increasing in recent years. The underlying mechanisms of EC tumorigenesis remain unclear, and efficient target therapies are lacking, for both of which feasible endometrial cancer animal models are essential but currently limited. Here, an organoid and genome editing-based strategy to generate primary, orthotopic, and driver-defined ECs in mice is reported. These models faithfully recapitulate the molecular and pathohistological characteristics of human diseases. The authors names these models and similar models for other cancers as organoid-initiated precision cancer models (OPCMs). Importantly, this approach can conveniently introduce any driver mutation or a combination of driver mutations. Using these models,it is shown that the mutations in Pik3ca and Pik3r1 cooperate with Pten loss to promote endometrial adenocarcinoma in mice. In contrast, the Kras G12D mutati led to endometrial squamous cell carcinoma. Then, tumor organoids are derived from these mouse EC models and performed high-throughput drug screening and validation. The results reveal distinct vulnerabilities of ECs with different mutations. Taken together, this study develops a multiplexing approach to model EC in mice and demonstrates its value for understanding the pathology of and exploring the potential treatments for this malignancy.

The genetic information coded in DNA leads to trait innovation via a gene regulatory network (GRN) in development. Here, we developed a conserved non-coding element interpretation method to integrate multi-omics data into gene regulatory network (CNEReg) to investigate the ruminant multi-chambered stomach innovation. We generated paired expression and chromatin accessibility data during rumen and esophagus development in sheep, and revealed 1601 active ruminant-specific conserved non-coding elements (active-RSCNEs). To interpret the function of these active-RSCNEs, we defined toolkit transcription factors (TTFs) and modeled their regulation on rumen-specific genes via batteries of active-RSCNEs during development. Our developmental GRN revealed 18 TTFs and 313 active-RSCNEs regulating 7 rumen functional modules. Notably, 6 TTFs (OTX1, SOX21, HOXC8, SOX2, TP63, and PPARG), as well as 16 active-RSCNEs, functionally distinguished the rumen from the esophagus. Our study provides a systematic approach to understanding how gene regulation evolves and shapes complex traits by putting evo-devo concepts into practice with developmental multi-omics data.

Gastric cancer (GC) is one of the most frequent and lethal malignancies in the world. However, our understanding of the mechanisms underlying its initiation and progression is limited. Here, we generate a series of primary GC models in mice with genome-edited gastric organoids, which elucidate the genetic drivers for sequential transformation from dysplasia to well-differentiated and poorly differentiated GC. Further, we find that the orthotopic GC, but not the subcutaneous GC even with the same genetic drivers, display remote metastasis, suggesting critical roles of the microenvironment in GC metastasis. Through single-cell RNA-seq analyses and functional studies, we show that the interaction between fibronectin 1 on stomach-specific macrophages and integrin a6β4 on GC cells promotes remote metastases. Taken together, our studies propose a strategy to model GC and dissect the genetic and microenvironmental factors driving the full-range gastric tumorigenesis.

Epigenetic regulations play crucial roles in the pathogenesis of metabolic-associated fatty liver disease; therefore, elucidating the biological functions of differential miRNAs helps us to understand the pathogenesis. Herein, we discovered miR-337-3p was decreased in patients with NAFLD from Gene Expression Omnibus dataset, which was replicated in various cell and mouse models with lipid disorders. Subsequently, overexpression of miR-337-3p in vivo could ameliorate hepatic lipid accumulation, reduce fasting blood glucose, and improve insulin resistance. Meanwhile, we determined miR-337-3p might influence multiple genes involved in glycolipid metabolism through mass spectrometry detection, bioinformatics analysis, and experimental verification. Finally, we selected HMGCR as a representative example to investigate the molecular mechanism of miR-337-3p regulating these genes, where the seed region of miR-337-3p bound to 3'UTR of HMGCR to inhibit HMGCR translation. In conclusion, we discovered a new function of miR-337-3p in glycolipid metabolism and that might be a new therapeutic target of MAFLD.

An accurate prediction of the intracranial infiltration tendency and drug response of individual glioblastoma (GBM) cells is essential for personalized prognosis and treatment for this disease. However, the clinical utility of mouse patient-derived orthotopic xenograft (PDOX) models remains limited given current technical constraints, including difficulty in generating sufficient sample numbers from small tissue samples and a long latency period for results. To overcome these issues, we established zebrafish GBM xenografts of diverse origin, which can tolerate intracranial engraftment and maintain their unique histological features. Subsequent single-cell RNA-sequencing (scRNA-seq) analysis confirmed significant transcriptional identity to that of invading GBM microtumors observed in the proportionally larger brains of model animals and humans. Endothelial scRNA-seq confirmed that the zebrafish blood-brain barrier is homologous to the mammalian blood-brain barrier. Finally, we established a rapid and efficient zebrafish PDOX (zPDOX) model, which can predict long-term outcomes of GBM patients within 20 days. The zPDOX model provides a novel avenue for precision medicine of GBM, especially for the evaluation of intracranial infiltration tendency and prediction of individual drug sensitivity.

DNA methylation-related genes, including TET2, IDH2, and DNMT3A are highly frequently mutated in angioimmunoblastic T-cell lymphoma (AITL), an aggressive malignancy of T follicular helper (Tfh) cells associated with aberrant immune features. It has been shown that TET2 loss cooperates with RHOAG17V to promote AITL in mice but the functional role of DNMT3A mutations in AITL remains unclear. Here, we report that DNMT3AR882H, the most common mutation of DNMT3A in AITL, accelerates the development of Tet2-/-; RHOAG17V AITL in mice, indicated by the expansion of malignant Tfh cells and aberrant B cells, skin rash, and significantly shortened disease-free survival. To understand the underlying cellular and molecular mechanisms, we performed single-cell transcriptome analyses of lymph nodes of mice transplanted with Tet2-/-, Tet2-/-; RHOAG17V or DNMT3AR882H; Tet2-/-; RHOAG17V hematopoietic stem and progenitor cells. These single-cell landscapes reveal that DNMT3A mutation further activates Tfh cells and leads to rapid and terminal differentiation of B cells, probably through enhancing the interacting PD1/PD-L1, ICOS/ICOSL, CD28/CD86, and ICAM1/ITGAL pairs. Our study establishes the functional roles of DNMT3A mutation in AITL and sheds light on the molecular mechanisms of this disease.

Pain and anxiety comorbidities are a common health problem, but the neural mechanisms underlying comorbidity remain unclear. We propose that comorbidity implies that similar brain regions and neural circuits, with the lateral septum (LS) as a major candidate, process pain and anxiety. From results of behavioral and neurophysiological experiments combined with selective LS manipulation in mice, we find that LS GABAergic neurons were critical for both pain and anxiety. Selective activation of LS GABAergic neurons induced hyperalgesia and anxiety-like behaviors. In contrast, selective inhibition of LS GABAergic neurons reduced nocifensive withdrawal responses and anxiety-like behaviors. This was found in two mouse models, one for chronic inflammatory pain (induced by complete Freund's adjuvant) and one for anxiety (induced by chronic restraint stress). Additionally, using TetTag chemogenetics to functionally mark LS neurons, we found that activation of LS neurons by acute pain stimulation could induce anxiety-like behaviors and vice versa. Furthermore, we show that LS GABAergic projection to the lateral hypothalamus (LH) plays an important role in the regulation of pain and anxiety comorbidities. Our study revealed that LS GABAergic neurons, and especially the LSGABAergic-LH circuit, are a critical to the modulation of pain and anxiety comorbidities.

Geese are extremely adept in utilizing plant-derived roughage within their diet. However, the intestinal microbiome of geese remains limited, especially the dietary effect on microbial diversity. Cassava foliage was widely used in animal feed, but little information is available for geese. In this study, the geese were fed with control diet (CK), experimental diet supplemented with 5% cassava foliage (CF5) or 10% (CF10) for 42 days, respectively. The cecal samples were collected after animals were killed. High-throughput sequencing technology was used to investigate the microbial diversity in the caecum of geese with different dietary supplements. Taxonomic analysis indicated that the predominant phyla were distinct with different dietary treatments. The phyla Firmicutes (51.4%), Bacteroidetes (29.55%) and Proteobacteria (7.90%) were dominant in the CK group, but Bacteroidetes (65.19% and 67.29%,) Firmicutes (18.01% and 17.39%), Proteobacteria (8.72% and 10.18%), Synergistete (2.51% and 1.76%) and Spirochaetes (2.60% and 1.46%) were dominant in CF5 and CF10 groups. The abundance of Firmicutes was negatively correlated with the supplementation of cassava foliage. However, the abundance of Bacteroidetes and Proteobacteria were positively correlated with the supplementation of cassava foliage. Our study also revealed that the microbial communities were significantly different at genus levels. Genes related to nutrient and energy metabolism, immunity and signal transduction pathways were primarily enriched by the microbiome.

Preweaned rumen development is vital for animal health and efficient fermentation. In this study, we integrated ruminal transcriptomic and metagenomic data to explore the dynamics of rumen functions, microbial colonization, and their functional interactions during the first 8 weeks of life in goats.

Immunophenotyping of the tumor microenvironment (TME) is essential for enhancing immunotherapy efficacy. However, strategies for characterizing the TME exhibit significant heterogeneity. Here, we show that endoplasmic reticular oxidoreductase-1α (ERO1A) mediates an immune-suppressive TME and attenuates the response to PD-1 blockade. Ablation of ERO1A in tumor cells substantially incites anti-tumor T cell immunity and promotes the efficacy of aPD-1 in therapeutic models. Single-cell RNA-sequencing analyses confirm that ERO1A correlates with immunosuppression and dysfunction of CD8+ T cells along anti-PD-1 treatment. In human lung cancer, high ERO1A expression is associated with a higher risk of recurrence following neoadjuvant immunotherapy. Mechanistically, ERO1A ablation impairs the balance between IRE1α and PERK signaling activities and induces lethal unfolded protein responses in tumor cells undergoing endoplasmic reticulum stress, thereby enhancing anti-tumor immunity via immunogenic cell death. These findings reveal how tumor ERO1A induces immunosuppression, highlighting its potential as a therapeutic target for cancer immunotherapy.

Poly (ADP-ribose) polymerase inhibitors (PARPi) are selectively active in ovarian cancer (OC) with homologous recombination (HR) deficiency (HRD) caused by mutations in BRCA1/2 and other DNA repair pathway members. We sought molecular targeted therapy that induce HRD in HR-proficient cells to induce synthetic lethality with PARPi and extend the utility of PARPi. Here, we demonstrate that lysine-specific demethylase 1 (LSD1) is an important regulator for OC. Importantly, genetic depletion or pharmacological inhibition of LSD1 induces HRD and sensitizes HR-proficient OC cells to PARPi in vitro and in multiple in vivo models. Mechanistically, LSD1 inhibition directly impairs transcription of BRCA1/2 and RAD51, three genes essential for HR, dependently of its canonical demethylase function. Collectively, our work indicates combination with LSD1 inhibitor could greatly expand the utility of PARPi to patients with HR-proficient tumor, warranting assessment in human clinical trials.

Deletions of chromosome 17p, where TP53 gene locates, are the most frequent chromosome alterations in human cancers and associated with poor outcomes in patients. Our previous work suggested that there were p53-independent mechanisms involved in chromosome 17p deletions-driven cancers. Here, we report that altered arachidonate metabolism, due to the deficiency of mouse Alox8 on chromosome 11B3 (homologous to human ALOX15B on chromosome 17p), contributes to the B cell malignancy. While the metabolites produced from lipoxygenase pathway reduced, chromosome 11B3 deletions or Alox8 loss, lead to upregulating its paralleling cyclooxygenase pathway, indicated by the increased levels of oncometabolite prostaglandin E2. Ectopic PGE2 prevented the apoptosis and differentiation of pre-B cells. Further studies revealed that Alox8 deficiency dramatically and specifically induced Cox-2(Ptgs2) gene expression. Repressing Cox-2 by its shRNAs impaired the tumorigenesis driven by Alox8 loss. And, in turn, tumor cells with Alox8 or 11B3 loss were sensitive to the COX-2 inhibitor celecoxib. This correlation between COX-2 upregulation and chromosome 17p deletions was consistent in human B-cell lymphomas. Hence, our studies reveal that the arachidonate metabolism abnormality with unbalanced ALOX and COX pathways underlies human cancers with 17p deletions and suggest new susceptibility for this disease.