With the gradual maturity of sequencing technology, many microbiome studies have published, driving the emergence and advance of related analysis tools. R language is the widely used platform for microbiome data analysis for powerful functions. However, tens of thousands of R packages and numerous similar analysis tools have brought major challenges for many researchers to explore microbiome data. How to choose suitable, efficient, convenient, and easy-to-learn tools from the numerous R packages has become a problem for many microbiome researchers. We have organized 324 common R packages for microbiome analysis and classified them according to application categories (diversity, difference, biomarker, correlation and network, functional prediction, and others), which could help researchers quickly find relevant R packages for microbiome analysis. Furthermore, we systematically sorted the integrated R packages (phyloseq, microbiome, MicrobiomeAnalystR, Animalcules, microeco, and amplicon) for microbiome analysis, and summarized the advantages and limitations, which will help researchers choose the appropriate tools. Finally, we thoroughly reviewed the R packages for microbiome analysis, summarized most of the common analysis content in the microbiome, and formed the most suitable pipeline for microbiome analysis. This paper is accompanied by hundreds of examples with 10,000 lines codes in GitHub, which can help beginners to learn, also help analysts compare and test different tools. This paper systematically sorts the application of R in microbiome, providing an important theoretical basis and practical reference for the development of better microbiome tools in the future. All the code is available at GitHub github.com/taowenmicro/EasyMicrobiomeR.

Age-dependent loss of skeletal muscle mass and function is a feature of sarcopenia, and increases the risk of many aging-related metabolic diseases. Here, we report phenotypic and single-nucleus transcriptomic analyses of non-human primate skeletal muscle aging. A higher transcriptional fluctuation was observed in myonuclei relative to other interstitial cell types, indicating a higher susceptibility of skeletal muscle fiber to aging. We found a downregulation of FOXO3 in aged primate skeletal muscle, and identified FOXO3 as a hub transcription factor maintaining skeletal muscle homeostasis. Through the establishment of a complementary experimental pipeline based on a human pluripotent stem cell-derived myotube model, we revealed that silence of FOXO3 accelerates human myotube senescence, whereas genetic activation of endogenous FOXO3 alleviates human myotube aging. Altogether, based on a combination of monkey skeletal muscle and human myotube aging research models, we unraveled the pivotal role of the FOXO3 in safeguarding primate skeletal muscle from aging, providing a comprehensive resource for the development of clinical diagnosis and targeted therapeutic interventions against human skeletal muscle aging and the onset of sarcopenia along with aging-related disorders.

Hair loss affects millions of people at some time in their life, and safe and efficient treatments for hair loss are a significant unmet medical need. We report that topical delivery of quercetin (Que) stimulates resting hair follicles to grow with rapid follicular keratinocyte proliferation and replenishes perifollicular microvasculature in mice. We construct dynamic single-cell transcriptome landscape over the course of hair regrowth and find that Que treatment stimulates the differentiation trajectory in the hair follicles and induces an angiogenic signature in dermal endothelial cells by activating HIF-1α in endothelial cells. Skin administration of a HIF-1α agonist partially recapitulates the pro-angiogenesis and hair-growing effects of Que. Together, these findings provide a molecular understanding for the efficacy of Que in hair regrowth, which underscores the translational potential of targeting the hair follicle niche as a strategy for regenerative medicine, and suggest a route of pharmacological intervention that may promote hair regrowth.

Emerging evidence suggests that intron-detaining transcripts (IDTs) are a nucleus-detained and polyadenylated mRNA pool for cell to quickly and effectively respond to environmental stimuli and stress. However, the underlying mechanisms of detained intron (DI) splicing are still largely unknown. Here, we suggest that post-transcriptional DI splicing is paused at the B^act state, an active spliceosome but not catalytically primed, which depends on Smad Nuclear Interacting Protein 1 (SNIP1) and RNPS1 (a serine-rich RNA binding protein) interaction. RNPS1 and B^act components preferentially dock at DIs and the RNPS1 docking is sufficient to trigger spliceosome pausing. Haploinsufficiency of Snip1 attenuates neurodegeneration and globally rescues IDT accumulation caused by a previously reported mutant U2 snRNA, a basal spliceosomal component. Snip1 conditional knockout in the cerebellum decreases DI splicing efficiency and causes neurodegeneration. Therefore, we suggest that SNIP1 and RNPS1 form a molecular brake to promote spliceosome pausing, and that its misregulation contributes to neurodegeneration.

The testis is pivotal for male reproduction, and its progressive functional decline in aging is associated with infertility. However, the regulatory mechanism underlying primate testicular aging remains largely elusive. Here, we resolve the aging-related cellular and molecular alterations of primate testicular aging by establishing a single-nucleus transcriptomic atlas. Gene-expression patterns along the spermatogenesis trajectory revealed molecular programs associated with attrition of spermatogonial stem cell reservoir, disturbed meiosis and impaired spermiogenesis along the sequential continuum. Remarkably, Sertoli cell was identified as the cell type most susceptible to aging, given its deeply perturbed age-associated transcriptional profiles. Concomitantly, downregulation of the transcription factor Wilms’ Tumor 1 (WT1), essential for Sertoli cell homeostasis, was associated with accelerated cellular senescence, disrupted tight junctions, and a compromised cell identity signature, which altogether may help create a hostile microenvironment for spermatogenesis. Collectively, our study depicts in-depth transcriptomic traits of non-human primate (NHP) testicular aging at single-cell resolution, providing potential diagnostic biomarkers and targets for therapeutic interventions against testicular aging and age-related male reproductive diseases.

Human pluripotent stem cells provide an inexhaustible model to study human embryogenesis in vitro. Recent studies have provided diverse models to generate human blastoids by self-organization of different pluripotent stem cells or somatic reprogramming intermediates. However, whether blastoids can be generated from other cell types or whether they can recapitulate postimplantation development in vitro is unknown. Here, we develop a strategy to generate human blastoids from heterogeneous intermediates with epiblast, trophectoderm, and primitive endoderm signatures of the primed-to-naïve conversion process, which resemble natural blastocysts in morphological architecture, composition of cell lineages, transcriptome, and lineage differentiation potential. In addition, these blastoids reflect many features of human peri-implantation and pregastrulation development when further cultured in an in vitro 3D culture system. In summary, our study provides an alternative strategy to generate human blastoids and offers insights into human early embryogenesis by modeling peri- and postimplantation development in vitro.

Mammals exhibit limited heart regeneration ability, which can lead to heart failure after myocardial infarction. In contrast, zebrafish exhibit remarkable cardiac regeneration capacity. Several cell types and signaling pathways have been reported to participate in this process. However, a comprehensive analysis of how different cells and signals interact and coordinate to regulate cardiac regeneration is unavailable. We collected major cardiac cell types from zebrafish and performed high-precision single-cell transcriptome analyses during both development and post-injury regeneration. We revealed the cellular heterogeneity as well as the molecular progress of cardiomyocytes during these processes, and identified a subtype of atrial cardiomyocyte exhibiting a stem-like state which may transdifferentiate into ventricular cardiomyocytes during regeneration. Furthermore, we identified a regeneration-induced cell (RIC) population in the epicardium-derived cells (EPDC), and demonstrated Angiopoietin 4 (Angpt4) as a specific regulator of heart regeneration. angpt4 expression is specifically and transiently activated in RIC, which initiates a signaling cascade from EPDC to endocardium through the Tie2-MAPK pathway, and further induces activation of cathepsin K in cardiomyocytes through RA signaling. Loss of angpt4 leads to defects in scar tissue resolution and cardiomyocyte proliferation, while overexpression of angpt4 accelerates regeneration. Furthermore, we found that ANGPT4 could enhance proliferation of neonatal rat cardiomyocytes, and promote cardiac repair in mice after myocardial infarction, indicating that the function of Angpt4 is conserved in mammals. Our study provides a mechanistic understanding of heart regeneration at single-cell precision, identifies Angpt4 as a key regulator of cardiomyocyte proliferation and regeneration, and offers a novel therapeutic target for improved recovery after human heart injuries.

Lipophagy, the selective engulfment of lipid droplets (LDs) by autophagosomes for lysosomal degradation, is critical to lipid and energy homeostasis. Here we show that the lipid transfer protein ORP8 is located on LDs and mediates the encapsulation of LDs by autophagosomal membranes. This function of ORP8 is independent of its lipid transporter activity and is achieved through direct interaction with phagophore-anchored LC3/GABARAPs. Upon lipophagy induction, ORP8 has increased localization on LDs and is phosphorylated by AMPK, thereby enhancing its affinity for LC3/GABARAPs. Deletion of ORP8 or interruption of ORP8-LC3/GABARAP interaction results in accumulation of LDs and increased intracellular triglyceride. Overexpression of ORP8 alleviates LD and triglyceride deposition in the liver of ob/ob mice, and Osbpl8^-/- mice exhibit liver lipid clearance defects. Our results suggest that ORP8 is a lipophagy receptor that plays a key role in cellular lipid metabolism.

As an important part of tumor microenvironment, neutrophils are poorly understood due to their spatiotemporal heterogeneity in tumorigenesis. Here we defined, at single-cell resolution, CD44⁻CXCR2⁻ neutrophils as tumor-specific neutrophils (tsNeus) in both mouse and human gastric cancer (GC). We uncovered a Hippo regulon in neutrophils with unique YAP signature genes (e.g., ICAM1, CD14, EGR1) distinct from those identified in epithelial and/or cancer cells. Importantly, knockout of YAP/TAZ in neutrophils impaired their differentiation into CD54⁺ tsNeus and reduced their antitumor activity, leading to accelerated GC progression. Moreover, the relative amounts of CD54⁺ tsNeus were found to be negatively associated with GC progression and positively associated with patient survival. Interestingly, GC patients receiving neoadjuvant chemotherapy had increased numbers of CD54⁺ tsNeus. Furthermore, pharmacologically enhancing YAP activity selectively activated neutrophils to suppress refractory GC, with no significant inflammation-related side effects. Thus, our work characterized tumor-specific neutrophils in GC and revealed an essential role of YAP/TAZ-CD54 axis in tsNeus, opening a new possibility to develop neutrophil-based antitumor therapeutics.

Approximately 140 million people worldwide are homozygous carriers of APOE4 (ϵ4), a strong genetic risk factor for late onset familial and sporadic Alzheimer’s disease (AD), 91% of whom will develop AD at earlier age than heterozygous carriers and noncarriers. Susceptibility to AD could be reduced by targeted editing of APOE4, but a technical basis for controlling the off-target effects of base editors is necessary to develop low-risk personalized gene therapies. Here, we first screened eight cytosine base editor variants at four injection stages (from 1- to 8-cell stage), and found that FNLS-YE1 variant in 8-cell embryos achieved the comparable base conversion rate (up to 100%) with the lowest bystander effects. In particular, 80% of AD-susceptible ϵ4 allele copies were converted to the AD-neutral ϵ3 allele in human ϵ4-carrying embryos. Stringent control measures combined with targeted deep sequencing, whole genome sequencing, and RNA sequencing showed no DNA or RNA off-target events in FNLS-YE1-treated human embryos or their derived stem cells. Furthermore, base editing with FNLS-YE1 showed no effects on embryo development to the blastocyst stage. Finally, we also demonstrated FNLS-YE1 could introduce known protective variants in human embryos to potentially reduce human susceptivity to systemic lupus erythematosus and familial hypercholesterolemia. Our study therefore suggests that base editing with FNLS-YE1 can efficiently and safely introduce known preventive variants in 8-cell human embryos, a potential approach for reducing human susceptibility to AD or other genetic diseases.

The clustered regularly interspaced short palindromic repeats (CRISPR)-Cas9 system has been widely used for genome engineering and transcriptional regulation in many different organisms. Current CRISPR-activation (CRISPRa) platforms often require multiple components because of inefficient transcriptional activation. Here, we fused different phase-separation proteins to dCas9-VPR (dCas9-VP64-P65-RTA) and observed robust increases in transcriptional activation efficiency. Notably, human NUP98 (nucleoporin 98) and FUS (fused in sarcoma) IDR domains were best at enhancing dCas9-VPR activity, with dCas9-VPR-FUS IDR (VPRF) outperforming the other CRISPRa systems tested in this study in both activation efficiency and system simplicity. dCas9-VPRF overcomes the target strand bias and widens gRNA designing windows without affecting the off-target effect of dCas9-VPR. These findings demonstrate the feasibility of using phase-separation proteins to assist in the regulation of gene expression and support the broad appeal of the dCas9-VPRF system in basic and clinical applications.