Smart-seq2 for sensitive full-length transcriptome profiling in single cells

A comprehensive review of the applications of RNA sequencing in celiac disease research

RNA sequencing (RNA-seq) has undergone substantial advancements in recent decades and has emerged as a vital technique for profiling the transcriptome. The transition from bulk sequencing to single-cell and spatial approaches has facilitated the achievement of higher precision at cell resolution. It provides valuable biological knowledge about individual immune cells and aids in the discovery of the molecular mechanisms that contribute to the development of autoimmune diseases. Celiac disease (CeD) is an autoimmune disorder characterized by a strong immune response to gluten consumption. RNA-seq has led to significantly advanced research in multiple fields, particularly in CeD research. It has been instrumental in studies involving comparative transcriptomics, nutritional genomics and wheat research, cancer research in the context of CeD, genetic and noncoding RNA-mediated epigenetic insights, disease monitoring and biomarker discovery, regulation of mitochondrial functions, therapeutic target identification and drug mechanism of action, dietary factors, immune cell profiling and the immune landscape. This review offers a comprehensive examination of recent RNA-seq technology research in the field of CeD, highlighting future challenges and opportunities for its application.

An atlas of small non-coding RNAs in human preimplantation development

Understanding the molecular circuitries that govern early embryogenesis is important, yet our knowledge of these in human preimplantation development remains limited. Small non-coding RNAs (sncRNAs) can regulate gene expression and thus impact blastocyst formation, however, the expression of specific biotypes and their dynamics during preimplantation development remains unknown. Here we identify the abundance of and kinetics of piRNA, rRNA, snoRNA, tRNA, and miRNA from embryonic day (E)3-7 and isolate specific miRNAs and snoRNAs of particular importance in blastocyst formation and pluripotency. These sncRNAs correspond to specific genomic hotspots: an enrichment of the chromosome 19 miRNA cluster (C19MC) in the trophectoderm (TE), and the chromosome 14 miRNA cluster (C14MC) and MEG8-related snoRNAs in the inner cell mass (ICM), which may serve as 'master regulators' of potency and lineage. Additionally, we observe a developmental transition with 21 isomiRs and in tRNA fragment (tRF) codon usage and identify two novel miRNAs. Our analysis provides a comprehensive measure of sncRNA biotypes and their corresponding dynamics throughout human preimplantation development, providing an extensive resource. Better understanding the sncRNA regulatory programmes in human embryogenesis will inform strategies to improve embryo development and outcomes of assisted reproductive technologies. We anticipate broad usage of our data as a resource for studies aimed at understanding embryogenesis, optimising stem cell-based models, assisted reproductive technology, and stem cell biology.

A ligation-independent sequencing method reveals tRNA-derived RNAs with blocked 3' termini

Despite the numerous sequencing methods available, the diversity in RNA size and chemical modification makes it difficult to capture all RNAs in a cell. We developed a method that combines quasi-random priming with template switching to construct sequencing libraries from RNA molecules of any length and with any type of 3' modifications, allowing for the sequencing of virtually all RNA species. Our ligation-independent detection of all types of RNA (LIDAR) is a simple, effective tool to identify and quantify all classes of coding and non-coding RNAs. With LIDAR, we comprehensively characterized the transcriptomes of mouse embryonic stem cells, neural progenitor cells, mouse tissues, and sperm. LIDAR detected a much larger variety of tRNA-derived RNAs (tDRs) compared with traditional ligation-dependent sequencing methods and uncovered tDRs with blocked 3' ends that had previously escaped detection. Therefore, LIDAR can capture all RNAs in a sample and uncover RNA species with potential regulatory functions.

Emerging methods and applications in 3D genomics

Since the advent of Hi-C in 2009, a plethora of high-throughput sequencing methods have emerged to profile the three-dimensional (3D) organization of eukaryotic genomes, igniting the era of 3D genomics. In recent years, the genomic resolution achievable by these approaches has dramatically increased and several single-cell versions of Hi-C have been developed. Moreover, a new repertoire of tools not based on proximity ligation of digested chromatin has emerged, enabling the investigation of the higher-order organization of chromatin in the nucleus. In this review, we summarize the expanding portfolio of 3D genomic technologies, highlighting recent developments and applications from the past three years. Lastly, we present an outlook of where this technology-driven field might be headed.

Maintaining mitochondrial DNA copy number mitigates ROS-induced oocyte decline and female reproductive aging

Oocytes play a crucial role in transmitting maternal mitochondrial DNA (mtDNA), essential for the continuation of species. However, the effects of mitochondrial reactive oxygen species (ROS) on mammalian oocyte maturation and mtDNA maintenance remain unclear. We investigated this by conditionally knocking out the Sod2 gene in primordial follicles, elevating mitochondrial matrix ROS levels from early oocyte stages. Our data indicates that reproductive aging in Sod2 conditional knockout females begins at 6 months, with oxidative stress impairing oocyte quality, particularly affecting OXPHOS complex II and mtDNA-encoded mRNA levels. Despite unchanged mtDNA mutation load, mtDNA copy numbers exhibited significant variations. Strikingly, reducing mtDNA copy numbers by reducing mtSSB protein, crucial for mtDNA replication, accelerated reproductive aging onset to three months, underscoring the critical role of mtDNA copy number maintenance under oxidative stress conditions. This research provides new insights into the relationship among mitochondrial ROS, mtDNA, and reproductive aging, offering potential strategies for delaying aging-related fertility decline.

Private information leakage from single-cell count matrices

The increase in publicly available human single-cell datasets, encompassing millions of cells from many donors, has significantly enhanced our understanding of complex biological processes. However, the accessibility of these datasets raises significant privacy concerns. Due to the inherent noise in single-cell measurements and the scarcity of population-scale single-cell datasets, recent private information quantification studies have focused on bulk gene expression data sharing. To address this gap, we demonstrate that individuals in single-cell gene expression datasets are vulnerable to linking attacks, where attackers can infer their sensitive phenotypic information using publicly available tissue or cell-type-specific expression quantitative trait loci (eQTLs) information. We further develop a method for genotype prediction and genotype-phenotype linking that remains effective without relying on eQTL information. We show that variants from one study can be exploited to uncover private information about individuals in another study.

Single-cell profiling identifies LIN28A mRNA targets in the mouse pluripotent-to-2C-like transition and somatic cell reprogramming

RNA-binding proteins (RBPs) regulate totipotency, pluripotency maintenance, and induction. The intricacies of how they modulate these processes through their interaction with RNAs remain to be elucidated. Here we employed Targets of RBPs Identified By Editing (TRIBE) with single-cell resolution (scTRIBE) to profile the mRNA targets of the key pluripotency regulator LIN28A in mouse embryonic stem cells (ESCs), 2-cell embryo-like cells (2CLCs) and somatic cell reprogramming. LIN28A is known to act by controlling the maturation of the let-7 microRNA but, in addition, it binds to multiple mRNAs and influences their stability and translation efficiency. However, the mRNA targets of LIN28A in 2CLCs and reprogramming are unclear. Through quantitative single-cell analysis of the scTRIBE dataset, we observed a marked increase in the binding of LIN28A to mRNAs of ribosome biogenesis factors and a selected group of totipotency factors in 2CLCs within ESC cultures. Our results suggest that LIN28A extends the half-life of at least some of these mRNAs, providing new insights into its role in the totipotent state. We also uncovered the distinct trajectory-specific LIN28A-mRNA networks in reprogramming, helping explain how LIN28A facilitates the mesenchymal-to-epithelial transition and pluripotency acquisition. Our study not only clarifies the multifunctional role of LIN28A in these processes but also highlights the importance of decoding RNA-protein interactions at the single-cell level.

Introducing synthetic thermostable RNase inhibitors to single-cell RNA-seq

Single-cell RNA-sequencing (scRNAseq) is revolutionizing biomedicine, propelled by advances in methodology, ease of use, and cost reduction of library preparation. Over the past decade, there have been remarkable technical improvements in most aspects of single-cell transcriptomics. Yet, little to no progress has been made in advancing RNase inhibition despite maintained RNA integrity being critical during cell collection, storage, and cDNA library generation. Here, we demonstrate that a synthetic thermostable RNase inhibitor (SEQURNA) yields single-cell libraries of equal or superior quality compared to ubiquitously used protein-based recombinant RNase inhibitors (RRIs). Importantly, the synthetic RNase inhibitor provides additional unique improvements in reproducibility and throughput, enables new experimental workflows including retained RNase inhibition throughout heat cycles, and can reduce the need for dry-ice transports. In summary, replacing RRIs represents a substantial advancement in the field of single-cell transcriptomics.

Intercellular extrachromosomal DNA copy-number heterogeneity drives neuroblastoma cell state diversity

Neuroblastoma exhibits significant inter- and intra-tumor genetic heterogeneity and varying clinical outcomes. Extrachromosomal DNAs (ecDNAs) may drive this heterogeneity by independently segregating during cell division, leading to rapid oncogene amplification. While ecDNA-mediated oncogene amplification is linked to poor prognosis in various cancers, the effects of ecDNA copy-number heterogeneity on intermediate phenotypes are poorly understood. Here, we leverage DNA and RNA sequencing from the same single cells in cell lines and neuroblastoma patients to investigate these effects. By analyzing ecDNA amplicon structures, we reveal extensive intercellular ecDNA copy-number heterogeneity. We also provide direct evidence of how this heterogeneity influences the expression of cargo genes, including MYCN and its downstream targets, and the overall transcriptional state of neuroblastoma cells. Our findings highlight the role of ecDNA copy number in promoting rapid adaptability of cellular states within tumors, underscoring the need for ecDNA-specific treatment strategies to address tumor formation and adaptation.

Fasting shapes chromatin architecture through an mTOR/RNA Pol I axis

Chromatin architecture is a fundamental mediator of genome function. Fasting is a major environmental cue across the animal kingdom, yet how it impacts three-dimensional (3D) genome organization is unknown. Here we show that fasting induces an intestine-specific, reversible and large-scale spatial reorganization of chromatin in Caenorhabditis elegans. This fasting-induced 3D genome reorganization requires inhibition of the nutrient-sensing mTOR pathway, acting through the regulation of RNA Pol I, but not Pol II nor Pol III, and is accompanied by remodelling of the nucleolus. By uncoupling the 3D genome configuration from the animal's nutritional status, we find that the expression of metabolic and stress-related genes increases when the spatial reorganization of chromatin occurs, showing that the 3D genome might support the transcriptional response in fasted animals. Our work documents a large-scale chromatin reorganization triggered by fasting and reveals that mTOR and RNA Pol I shape genome architecture in response to nutrients.

Deciphering breast cancer dynamics: insights from single-cell and spatial profiling in the multi-omics era

As one of the most common tumors in women, the pathogenesis and tumor heterogeneity of breast cancer have long been the focal point of research, with the emergence of tumor metastasis and drug resistance posing persistent clinical challenges. The emergence of single-cell sequencing (SCS) technology has introduced novel approaches for gaining comprehensive insights into the biological behavior of malignant tumors. SCS is a high-throughput technology that has rapidly developed in the past decade, providing high-throughput molecular insights at the individual cell level. Furthermore, the advent of multitemporal point sampling and spatial omics also greatly enhances our understanding of cellular dynamics at both temporal and spatial levels. The paper provides a comprehensive overview of the historical development of SCS, and highlights the most recent advancements in utilizing SCS and spatial omics for breast cancer research. The findings from these studies will serve as valuable references for future advancements in basic research, clinical diagnosis, and treatment of breast cancer.

Claustrum and dorsal endopiriform cortex complex cell-identity is determined by Nurr1 and regulates hallucinogenic-like states in mice

The Claustrum/dorsal endopiriform cortex complex (CLA) is an enigmatic brain region with extensive glutamatergic projections to multiple cortical areas. The transcription factor Nurr1 is highly expressed in the CLA, but its role in this region is not understood. By using conditional gene-targeted mice, we show that Nurr1 is a crucial regulator of CLA neuron identity. Although CLA neurons remain intact in the absence of Nurr1, the distinctive gene expression pattern in the CLA is abolished. CLA has been hypothesized to control hallucinations, but little is known of how the CLA responds to hallucinogens. After the deletion of Nurr1 in the CLA, both hallucinogen receptor expression and signaling are lost. Furthermore, functional ultrasound and Neuropixel electrophysiological recordings revealed that the hallucinogenic-receptor agonists' effects on functional connectivity between prefrontal and sensorimotor cortices are altered in Nurr1-ablated mice. Our findings suggest that Nurr1-targeted strategies provide additional avenues for functional studies of the CLA.

The non-canonical bivalent gene Wfdc15a controls spermatogenic protease and immune homeostasis

Male infertility can be caused by chromosomal abnormalities, mutations and epigenetic defects. Epigenetic modifiers pre-program hundreds of spermatogenic genes in spermatogonial stem cells (SSCs) for expression later in spermatids, but it remains mostly unclear whether and how those genes are involved in fertility. Here, we report that Wfdc15a, a WFDC family protease inhibitor pre-programmed by KMT2B, is essential for spermatogenesis. We found that Wfdc15a is a non-canonical bivalent gene carrying both H3K4me3 and facultative H3K9me3 in SSCs, but is later activated along with the loss of H3K9me3 and acquisition of H3K27ac during meiosis. We show that WFDC15A deficiency causes defective spermiogenesis at the beginning of spermatid elongation. Notably, depletion of WFDC15A causes substantial disturbance of the testicular protease-antiprotease network and leads to an orchitis-like inflammatory response associated with TNFα expression in round spermatids. Together, our results reveal a unique epigenetic program regulating innate immunity crucial for fertility.

NACC2, a molecular effector of miR-132 regulation at the interface between adult neurogenesis and Alzheimer's disease

The generation of new neurons at the hippocampal neurogenic niche, known as adult hippocampal neurogenesis (AHN), and its impairment, have been implicated in Alzheimer's disease (AD). MicroRNA-132 (miR-132), the most consistently downregulated microRNA (miRNA) in AD, was recently identified as a potent regulator of AHN, exerting multilayered proneurogenic effects in adult neural stem cells (NSCs) and their progeny. Supplementing miR-132 in AD mouse brain restores AHN and relevant memory deficits, yet the exact mechanisms involved are still unknown. Here, we identify NACC2 as a novel miR-132 target implicated in both AHN and AD. miR-132 deficiency in mouse hippocampus induces Nacc2 expression and inflammatory signaling in adult NSCs. We show that miR-132-dependent regulation of NACC2 is involved in the initial stages of human NSC differentiation towards astrocytes and neurons. Later, NACC2 function in astrocytic maturation becomes uncoupled from miR-132. We demonstrate that NACC2 is present in reactive astrocytes surrounding amyloid plaques in mouse and human AD hippocampus, and that there is an anticorrelation between miR-132 and NACC2 levels in AD and upon induction of inflammation. Unraveling the molecular mechanisms by which miR-132 regulates neurogenesis and cellular reactivity in AD, will provide valuable insights towards its possible application as a therapeutic target.

Efficient genome replication in influenza A virus requires NS2 and sequence beyond the canonical promoter

Influenza A virus encodes promoters in both the sense and antisense orientations. These support the generation of new genomes, antigenomes, and mRNA transcripts. Using minimal replication assays-transfections with viral polymerase, nucleoprotein, and a genomic template-the influenza promoter sequences were identified as 13nt at the 5' end of the viral genomic RNA (U13) and 12nt at the 3' end (U12). Other than the fourth 3' nucleotide, the U12 and U13 sequences are identical between all eight RNA molecules that comprise the segmented influenza genome. Despite possessing identical promoters, individual segments can exhibit different transcriptional dynamics during infection. However flu promoter sequences were defined in experiments without influenza NS2, a protein which modulates transcription and replication differentially between genomic segments. This suggests that the identity of the "complete" promoter may depend on NS2. Here we assess how internal sequences of two genomic segments, HA and PB1, may contribute to NS2-dependent replication as well as map such interactions down to individual nucleotides in PB1. We find that the expression of NS2 significantly alters sequence requirements for efficient replication beyond the identical U12 and U13 sequence, providing a mechanism for the divergent replication and transcription dynamics across the influenza A virus genome.

An integrated approach identifies the molecular underpinnings of murine anterior visceral endoderm migration

The anterior visceral endoderm (AVE) differs from the surrounding visceral endoderm (VE) in its migratory behavior and ability to restrict primitive streak formation to the opposite side of the mouse embryo. To characterize the molecular bases for the unique properties of the AVE, we combined single-cell RNA sequencing of the VE prior to and during AVE migration with phosphoproteomics, high-resolution live-imaging, and short-term lineage labeling and intervention. This identified the transient nature of the AVE with attenuation of "anteriorizing" gene expression as cells migrate and the emergence of heterogeneities in transcriptional states relative to the AVE's position. Using cell communication analysis, we identified the requirement of semaphorin signaling for normal AVE migration. Lattice light-sheet microscopy showed that Sema6D mutants have abnormalities in basal projections and migration speed. These findings point to a tight coupling between transcriptional state and position of the AVE and identify molecular controllers of AVE migration.

AMAR-seq: Automated Multimodal Sequencing of DNA Methylation, Chromatin Accessibility, and RNA Expression with Single-Cell Resolution

Parallel single-cell multimodal sequencing is the most intuitive and precise tool for cellular status research. In this study, we propose AMAR-seq to automate methylation, chromatin accessibility, and RNA expression coanalysis with single-cell precision. We validated the accuracy and robustness of AMAR-seq in comparison with standard single-omics methods. The high gene detection rate and genome coverage of AMAR-seq enabled us to establish a genome-wide gene expression regulatory atlas and triple-omics landscape with single base resolution and implement single-cell copy number variation analysis. Applying AMAR-seq to investigate the process of mouse embryonic stem cell differentiation, we revealed the dynamic coupling of the epigenome and transcriptome, which may contribute to unraveling the molecular mechanisms of early embryonic development. Collectively, we propose AMAR-seq for the in-depth and accurate establishment of single-cell multiomics regulatory patterns in a cost-effective, efficient, and automated manner, paving the way for insightful dissection of complex life processes.

Transcriptional re-programming of liver-resident iNKT cells into T-regulatory type-1-like liver iNKT cells involves extensive gene de-methylation

Unlike conventional CD4+ T cells, which are phenotypically and functionally plastic, invariant NKT (iNKT) cells generally exist in a terminally differentiated state. Naïve CD4+ T cells can acquire alternative epigenetic states in response to different cues, but it remains unclear whether peripheral iNKT cells are epigenetically stable or malleable. Repetitive encounters of liver-resident iNKT cells (LiNKTs) with alpha-galactosylceramide (αGalCer)/CD1d-coated nanoparticles (NPs) can trigger their differentiation into a LiNKT cell subset expressing a T regulatory type 1 (TR1)-like (LiNKTR1) transcriptional signature. Here we dissect the epigenetic underpinnings of the LiNKT-LiNKTR1 conversion as compared to those underlying the peptide-major histocompatibility complex (pMHC)-NP-induced T-follicular helper (TFH)-to-TR1 transdifferentiation process. We show that gene upregulation during the LINKT-to-LiNKTR1 cell conversion is associated with demethylation of gene bodies, inter-genic regions, promoters and distal gene regulatory elements, in the absence of major changes in chromatin exposure or deposition of expression-promoting histone marks. In contrast, the naïve CD4+ T cell-to-TFH differentiation process involves extensive remodeling of the chromatin and the acquisition of a broad repertoire of epigenetic modifications that are then largely inherited by TFH cell-derived TR1 cell progeny. These observations indicate that LiNKT cells are epigenetically malleable and particularly susceptible to gene de-methylation.

State of the Art in Sub-Phenotyping Midbrain Dopamine Neurons

Dopaminergic neurons in the ventral tegmental area (VTA) and the substantia nigra pars compacta (SNpc) comprise around 75% of all dopaminergic neurons in the human brain. While both groups of dopaminergic neurons are in close proximity in the midbrain and partially overlap, development, function, and impairments in these two classes of neurons are highly diverse. The molecular and cellular mechanisms underlying these differences are not yet fully understood, but research over the past decade has highlighted the need to differentiate between these two classes of dopaminergic neurons during their development and in the mature brain. This differentiation is crucial not only for understanding fundamental circuitry formation in the brain but also for developing therapies targeted to specific dopaminergic neuron classes without affecting others. In this review, we summarize the state of the art in our understanding of the differences between the dopaminergic neurons of the VTA and the SNpc, such as anatomy, structure, morphology, output and input, electrophysiology, development, and disorders, and discuss the current technologies and methods available for studying these two classes of dopaminergic neurons, highlighting their advantages, limitations, and the necessary improvements required to achieve more-precise therapeutic interventions.

Transcriptome analysis reveals that defects in cell cycle regulation contribute to preimplantation embryo arrest

Patients with preimplantation embryo arrest (PREMBA) often experience assisted reproductive failure primarily due to the lack of transferable embryos, and the molecular mechanisms underlying PREMBA remain unclear. In our study, the embryos from five women with recurrent preimplantation embryo arrest and three women with tubal factor infertility were used for single-embryo transcriptome sequencing. Meanwhile, the transcriptomes of normal human preimplantation embryos obtained from GSE36552 were utilized to perform a comparative analysis with the transcriptomes of PREMBA embryos. Our results showed dysregulation of the cell cycle phase transition might be a potential pathogenic factor contributing to PREMBA. Through integrated analysis of the differentially expressed genes (DEGs) and weighted gene co-expression network analysis (WGCNA), we identified a number of hub genes using the protein-protein interaction network. The top 5 hub genes were as follows: CCNB2, BUB1B, CDC25A, CCNB3, and PLK3. The expression of hub genes was validated in PREMBA embryos and donated embryos using RT-qPCR. The knockdown of Ccnb2 in mouse zygotes led to an increase in embryo fragmentation, a rise in apoptosis, and a reduction in blastocyst formation. Furthermore, silencing the expression of CDC25A in HEK293T cells resulted in a decrease in cell proliferation and an increase in apoptosis, providing further support for our findings. Our findings could predict the development outcomes of preimplantation embryos and be used as potential therapeutic targets to prevent recurrent failures of IVF/ICSI attempts.

Germline loss in C. elegans enhances longevity by disrupting adhesion between niche and stem cells

Ageing and fertility are intertwined. Germline loss extends the lifespan in various organisms, termed gonadal longevity. However, the original longevity signal from the somatic gonad remains poorly understood. Here, we focused on the interaction between germline stem cells (GSCs) and their niche, the distal tip cells (DTCs), to explore the barely known longevity signal from the somatic gonad in C. elegans. We found that removing germline disrupts the cell adhesions between GSC and DTC, causing a significant transcriptomic change in DTC through hmp-2/β-catenin and two GATA transcription factors, elt-3 and pqm-1 in this niche cell. Inhibiting elt-3 and pqm-1 in DTC suppresses gonadal longevity. Moreover, we further identified the TGF-β ligand, tig-2, as the cytokine from DTC upon the loss of germline, which evokes the downstream gonadal longevity signalling throughout the body. Our findings thus reveal the source of the longevity signalling in response to germline removal, highlighting the stem cell niche as a critical signalling hub in ageing.

The heterogeneity of erythroid cells: insight at the single-cell transcriptome level

Erythroid cells, the most prevalent cell type in blood, are one of the earliest products and permeate through the entire process of hematopoietic development in the human body, the oxygen-transporting function of which is crucial for maintaining overall health and life support. Previous investigations into erythrocyte differentiation and development have primarily focused on population-level analyses, lacking the single-cell perspective essential for comprehending the intricate pathways of erythroid maturation, differentiation, and the encompassing cellular heterogeneity. The continuous optimization of single-cell transcriptome sequencing technology, or single-cell RNA sequencing (scRNA-seq), provides a powerful tool for life sciences research, which has a particular superiority in the identification of unprecedented cell subgroups, the analyzing of cellular heterogeneity, and the transcriptomic characteristics of individual cells. Over the past decade, remarkable strides have been taken in the realm of single-cell RNA sequencing technology, profoundly enhancing our understanding of erythroid cells. In this review, we systematically summarize the recent developments in single-cell transcriptome sequencing technology and emphasize their substantial impact on the study of erythroid cells, highlighting their contributions, including the exploration of functional heterogeneity within erythroid populations, the identification of novel erythrocyte subgroups, the tracking of different erythroid lineages, and the unveiling of mechanisms governing erythroid fate decisions. These findings not only invigorate erythroid cell research but also offer new perspectives on the management of diseases related to erythroid cells.

Single-Cell RNA-Sequencing: Opening New Horizons for Breast Cancer Research

Breast cancer is the most prevalent malignant tumor among women with high heterogeneity. Traditional techniques frequently struggle to comprehensively capture the intricacy and variety of cellular states and interactions within breast cancer. As global precision medicine rapidly advances, single-cell RNA sequencing (scRNA-seq) has become a highly effective technique, revolutionizing breast cancer research by offering unprecedented insights into the cellular heterogeneity and complexity of breast cancer. This cutting-edge technology facilitates the analysis of gene expression profiles at the single-cell level, uncovering diverse cell types and states within the tumor microenvironment. By dissecting the cellular composition and transcriptional signatures of breast cancer cells, scRNA-seq provides new perspectives for understanding the mechanisms behind tumor therapy, drug resistance and metastasis in breast cancer. In this review, we summarized the working principle and workflow of scRNA-seq and emphasized the major applications and discoveries of scRNA-seq in breast cancer research, highlighting its impact on our comprehension of breast cancer biology and its potential for guiding personalized treatment strategies.

Application and research progress of single cell sequencing technology in leukemia

Leukemia is a malignant tumor with high heterogeneity and a complex evolutionary process. It is difficult to resolve the heterogeneity and clonal evolution of leukemia cells by applying traditional bulk sequencing techniques, thus preventing a deep understanding of the mechanisms of leukemia development and the identification of potential therapeutic targets. However, with the development and application of single-cell sequencing technology, it is now possible to investigate the gene expression profile, mutations, and epigenetic features of leukemia at the single-cell level, thus providing a new perspective for leukemia research. In this article, we review the recent applications and advances of single-cell sequencing technology in leukemia research, discuss its potential for enhancing our understanding of the mechanisms of leukemia development, discovering therapeutic targets and personalized treatment, and provide reference guidelines for the significance of this technology in clinical research.

Human Liver Organoids to Predict the Outcome of Kasai Portoenterostomy

BACKGROUND

Kasai portoenterostomy (KPE) remains the primary intervention for biliary atresia (BA), but its outcomes are highly variable. Reliable prognostic biomarkers remain elusive, complicating the management and prediction of postoperative progression.

METHOD

Liver biopsies from BA patients taken at and after KPE (post-KPE) were used to generate organoids for RNA-sequencing analysis. Control organoids were derived from non-BA livers. Differential gene expression and enrichment analyses were performed to assess post-KPE transcriptomic changes between native liver survivors (NLS) and patients who eventually became liver transplant recipients (LTR).

RESULTS

Organoid datasets: 70 from liver biopsies at KPE (10 patients), 112 from post-KPE livers (13 livers; 12 patients), and 47 from control livers (9 patients). At KPE, BA organoids displayed mainly hepatocyte expression, a trait notably reduced in control organoids. Similarly, post-KPE organoids from NLS revealed a significant decrease in hepatocyte expression features and an overall increase in cholangiocyte expression features. A similar hepatocyte-to-cholangiocyte expression transition was evidenced in paired liver organoids (at- and post-KPE) generated from an NLS. In contrast, post-KPE organoids from LTR maintained a high level of hepatocyte expression features.

CONCLUSION

Our study demonstrated that an elevated expression of hepatocyte features in KPE organoids may indicate aberrant cholangiocyte development in BA livers. In contrast, a post-KPE hepatocyte-to-cholangiocyte expression transition in NLS may imply effective biliary recovery. The lack of this transition in LTR organoids indicates ongoing disease progression, highlighting the potential for organoid-based transcriptomic profiling to inform KPE success and guide BA management.

LEVEL OF EVIDENCE

Level III.

Transcriptomic profiling of the oocyte-cumulus-granulosa cell complex from estrogen receptor β knockout mice

OBJECTIVE

To study the role of estrogen receptor β in follicle development and maturation and the response to gonadotropin stimulation aiming at superovulation.

DESIGN

Experimental study and transcriptomic analysis.

SETTING

Karolinka Institutet, medical university.

ANIMAL(S)

Healthy wild-type (WT) and estrogen receptor β knockout (Esr2-KO) female mice undergoing superovulation at 4 weeks, 7 weeks, and 6 months of age.

INTERVENTION(S)

Not applicable.

MAIN OUTCOME MEASURE(S)

Oocyte yield after superovulation, transcriptomic profiling of cumulus-granulosa cell complexes and oocytes, and immunohistochemical analyses.

RESULT(S)

Superovulation of estrogen receptor β (ERβ) knockout mice resulted in reduced oocyte yield at 6 months of age compared with WT mice, but younger mice had similar yields. RNA-seq analysis of cumulus cells from superovulated WT and Esr2-KO mice identified genes and pathways associated with among others adhesion, proliferation, Wnt-signaling, and placed ERβ in bipotential granulosa cell cluster. Loss of ERβ increased expression of the other estrogen receptors Esr1 and Gper1.

CONCLUSION(S)

Our results show that ERβ has an important role in regulating ovulation in response to exogenous gonadotropins in 6-month-old mice, but not in younger mice. Our transcriptomic and immunohistochemical observations suggest a dysregulation of the granulosa cell communication and lack of tight coordination between granulosa cell replication and antrum expansion. A significant upregulation of other estrogen receptors may support a compensatory mechanism sustaining fertility during younger age in Esr2-KO mice.

Transcriptomic analysis of the spatiotemporal axis of oogenesis and fertilization in C. elegans

Caenorhabditis elegans hermaphrodite presents a unique model to study the formation of oocytes. However, the size of the model animal and difficulties in retrieval of specific stages of the germline have obviated closer systematic studies of this process throughout the years. Here, we present a transcriptomic level analysis into the oogenesis of C. elegans hermaphrodites. We dissected a hermaphrodite gonad into seven sections corresponding to the mitotic distal region, the pachytene region, the diplotene region, the early diakinesis region and the 3 most proximal oocytes, and deeply sequenced the transcriptome of each of them along with that of the fertilized egg using a single-cell RNA-seq (scRNA-seq) protocol. We identified specific gene expression events as well as gene splicing events in finer detail along the gonad and provided novel insights into underlying mechanisms of the oogenesis process. Furthermore, through careful review of relevant research literature coupled with patterns observed in our analysis, we delineate transcripts that may serve functions in the interactions between the germline and cells of the somatic gonad. These results expand our knowledge of the transcriptomic space of the C. elegans germline and lay a foundation on which future studies of the germline can be based upon.

The transgenerational effects of maternal low-protein diet during lactation on offspring

Environmental factors such as diet and lifestyle can influence the health of both mothers and offspring. However, its transgenerational transmission and underlying mechanisms remain largely unknown. Here, using a maternal lactation-period low-protein diet (LPD) mouse model, we show that maternal LPD during lactation causes decreased survival and stunted growth, significantly reduces ovulation and litter size, and alters the gut microbiome in the female LPD-F1 offspring. The transcriptome of LPD-F1 metaphase II (MII) oocytes shows that differentially expressed genes are enriched in female pregnancy and multiple metabolic processes. Moreover, maternal LPD causes early stunted growth and impairs metabolic health, which is transmitted over two generations. The methylome alteration of LPD-F1 oocytes can be partly transmitted to the F2 oocytes. Together, our results reveal that LPD during lactation transgenerationally affects offspring health, probably via oocyte epigenetic changes.

Discovery of a monoclonal, high-affinity CD8+ T-cell clone following natural hepatitis C virus infection

CD8+ T cells recognizing their cognate antigen are typically recruited as a polyclonal population consisting of multiple clonotypes with varying T-cell receptor (TCR) affinity to the target peptide-major histocompatibility complex (pMHC) complex. Advances in single-cell sequencing have increased accessibility toward identifying TCRs with matched antigens. Here we present the discovery of a monoclonal CD8+ T-cell population with specificity for a hepatitis C virus (HCV)-derived human leukocyte antigen (HLA) class I epitope (HLA-B*07:02 GPRLGVRAT) which was isolated directly ex vivo from an individual with an episode of acutely resolved HCV infection. This population was absent before infection and underwent expansion and stable maintenance for at least 2 years after infection as measured by HLA-multimer staining. Furthermore, the monoclonal clonotype was characterized by an unusually long dissociation time (half-life = 794 s and koff = 5.73 × 10-4) for its target antigen when compared with previously published results. A comparison with related populations of HCV-specific populations derived from the same individual and a second individual suggested that high-affinity TCR-pMHC interactions may be inherent to epitope identity and shape the phenotype of responses which has implications for rational TCR selection and design in the age of personalized immunotherapies.

Single-cell RNA sequencing technology in human spermatogenesis: Progresses and perspectives

Spermatogenesis, a key part of the spermiation process, is regulated by a combination of key cells, such as primordial germ cells, spermatogonial stem cells, and somatic cells, such as Sertoli cells. Abnormal spermatogenesis can lead to azoospermia, testicular tumors, and other diseases related to male infertility. The application of single-cell RNA sequencing (scRNA-seq) technology in male reproduction is gradually increasing with its unique insight into deep mining and analysis. The data cover different periods of neonatal, prepubertal, pubertal, and adult stages. Different types of male infertility diseases including obstructive and non-obstructive azoospermia (NOA), Klinefelter Syndrome (KS), Sertoli Cell Only Syndrome (SCOS), and testicular tumors are also covered. We briefly review the principles and application of scRNA-seq and summarize the research results and application directions in spermatogenesis in different periods and pathological states. Moreover, we discuss the challenges of applying this technology in male reproduction and the prospects of combining it with other technologies.

Population Diversity at the Single-Cell Level

Population-scale single-cell genomics is a transformative approach for unraveling the intricate links between genetic and cellular variation. This approach is facilitated by cutting-edge experimental methodologies, including the development of high-throughput single-cell multiomics and advances in multiplexed environmental and genetic perturbations. Examining the effects of natural or synthetic genetic variants across cellular contexts provides insights into the mutual influence of genetics and the environment in shaping cellular heterogeneity. The development of computational methodologies further enables detailed quantitative analysis of molecular variation, offering an opportunity to examine the respective roles of stochastic, intercellular, and interindividual variation. Future opportunities lie in leveraging long-read sequencing, refining disease-relevant cellular models, and embracing predictive and generative machine learning models. These advancements hold the potential for a deeper understanding of the genetic architecture of human molecular traits, which in turn has important implications for understanding the genetic causes of human disease.

Cohesin prevents cross-domain gene coactivation

The contrast between the disruption of genome topology after cohesin loss and the lack of downstream gene expression changes instigates intense debates regarding the structure-function relationship between genome and gene regulation. Here, by analyzing transcriptome and chromatin accessibility at the single-cell level, we discover that, instead of dictating population-wide gene expression levels, cohesin supplies a general function to neutralize stochastic coexpression tendencies of cis-linked genes in single cells. Notably, cohesin loss induces widespread gene coactivation and chromatin co-opening tens of million bases apart in cis. Spatial genome and protein imaging reveals that cohesin prevents gene co-bursting along the chromosome and blocks spatial mixing of transcriptional hubs. Single-molecule imaging shows that cohesin confines the exploration of diverse enhancer and core promoter binding transcriptional regulators. Together, these results support that cohesin arranges nuclear topology to control gene coexpression in single cells.

Dynamics and timing of diversification events of ciliated eukaryotes from a large phylogenomic perspective

Ciliophora, an exceptionally diverse lineage of unicellular eukaryotes, exhibits a remarkable range of species richness across classes in the ciliate Tree of Life. In this study, we have acquired transcriptome and genome data from 40 representative species in seven ciliate classes. Utilizing 247 genes and 105 taxa, we devised a comprehensive phylogenomic tree for Ciliophora, encompassing over 60 % of orders and constituting the most extensive dataset of ciliate species to date. We established a robust phylogenetic framework that encompasses ambiguous taxa and the major classes within the phylum. Our findings support the monophyly of each of two subphyla (Postciliodesmatophora and Intramacronucleata), along with three subclades (Protocruzia, CONTHREEP, and SAPML) nested within Intramacronucleata, and elucidate evolutionary positions among the major classes within the phylum. Drawing on the robust ciliate Tree of Life and three constraints, we estimated the radiation of Ciliophora around 1175 Ma during the middle of the Proterozoic Eon, and most of the ciliate classes diverged from their sister lineage during the latter half of this period. Additionally, based on the time-calibrated tree and species richness pattern, we investigated net diversification rates of Ciliophora and its classes. The global net diversification rate for Ciliophora was estimated at 0.004979 species/Ma. Heterogeneity in net diversification rates was evident at the class level, with faster rates observed in Oligohymenophorea and Spirotrichea than other classes within the subclades CONTHREEP and SAPML, respectively. Notably, our analysis suggests that variations in net diversification rates, rather than clade ages, appear to contribute to the differences in species richness in Ciliophora at the class level.

Degenerative and regenerative peripheral processes are associated with persistent painful chemotherapy-induced neuropathies in males and females

This study investigated the time course of gene expression changes during the progression of persistent painful neuropathy caused by paclitaxel (PTX) in male and female mouse hindpaws and dorsal root ganglia (DRG). Bulk RNA-seq was used to examine these gene expression changes at 1, 16, and 31 days post-last PTX. At these time points, differentially expressed genes (DEGs) were predominantly related to the reduction or increase in epithelial, skin, bone, and muscle development and to angiogenesis, myelination, axonogenesis, and neurogenesis. These processes are accompanied by the regulation of DEGs related to the cytoskeleton, extracellular matrix organization, and cellular energy production. This gene plasticity during the progression of persistent painful neuropathy could be interpreted as a biological process linked to tissue regeneration/degeneration. In contrast, gene plasticity related to immune processes was minimal at 1-31 days after PTX. It was also noted that despite similarities in biological processes and pain chronicity between males and females, specific DEGs differed dramatically according to sex. The main conclusions of this study are that gene expression plasticity in hindpaw and DRG during PTX neuropathy progression similar to tissue regeneration and degeneration, minimally affects immune system processes and is heavily sex-dependent at the individual gene level.

Diversity of microglial transcriptional responses during opioid exposure and neuropathic pain

ABSTRACT

Microglia take on an altered morphology during chronic opioid treatment. This morphological change is broadly used to identify the activated microglial state associated with opioid side effects, including tolerance and opioid-induced hyperalgesia (OIH). Microglia display similar morphological responses in the spinal cord after peripheral nerve injury (PNI). Consistent with this observation, functional studies have suggested that microglia activated by opioids or PNI engage common molecular mechanisms to induce hypersensitivity. In this article, we conducted deep RNA sequencing (RNA-seq) and morphological analysis of spinal cord microglia in male mice to comprehensively interrogate transcriptional states and mechanistic commonality between multiple models of OIH and PNI. After PNI, we identify an early proliferative transcriptional event across models that precedes the upregulation of histological markers of microglial activation. However, we found no proliferative transcriptional response associated with opioid-induced microglial activation, consistent with histological data, indicating that the number of microglia remains stable during morphine treatment, whereas their morphological response differs from PNI models. Collectively, these results establish the diversity of pain-associated microglial transcriptomic responses and point towards the targeting of distinct insult-specific microglial responses to treat OIH, PNI, or other central nervous system pathologies.

Mining alternative splicing patterns in scRNA-seq data using scASfind

Single-cell RNA-seq (scRNA-seq) is widely used for transcriptome profiling, but most analyses focus on gene-level events, with less attention devoted to alternative splicing. Here, we present scASfind, a novel computational method to allow for quantitative analysis of cell type-specific splicing events using full-length scRNA-seq data. ScASfind utilizes an efficient data structure to store the percent spliced-in value for each splicing event. This makes it possible to exhaustively search for patterns among all differential splicing events, allowing us to identify marker events, mutually exclusive events, and events involving large blocks of exons that are specific to one or more cell types.

scCross: a deep generative model for unifying single-cell multi-omics with seamless integration, cross-modal generation, and in silico exploration

Single-cell multi-omics data reveal complex cellular states, providing significant insights into cellular dynamics and disease. Yet, integration of multi-omics data presents challenges. Some modalities have not reached the robustness or clarity of established transcriptomics. Coupled with data scarcity for less established modalities and integration intricacies, these challenges limit our ability to maximize single-cell omics benefits. We introduce scCross, a tool leveraging variational autoencoders, generative adversarial networks, and the mutual nearest neighbors (MNN) technique for modality alignment. By enabling single-cell cross-modal data generation, multi-omics data simulation, and in silico cellular perturbations, scCross enhances the utility of single-cell multi-omics studies.

Compound models and Pearson residuals for single-cell RNA-seq data without UMIs

Recent work employed Pearson residuals from Poisson or negative binomial models to normalize UMI data. To extend this approach to non-UMI data, we model the additional amplification step with a compound distribution: we assume that sequenced RNA molecules follow a negative binomial distribution, and are then replicated following an amplification distribution. We show how this model leads to compound Pearson residuals, which yield meaningful gene selection and embeddings of Smart-seq2 datasets. Further, we suggest that amplification distributions across several sequencing protocols can be described by a broken power law. The resulting compound model captures previously unexplained overdispersion and zero-inflation patterns in non-UMI data.

Spatial transcriptomic profiling of isolated microregions in tissue sections utilizing laser-induced forward transfer

Profiling gene expression while preserving cell locations aids in the comprehensive understanding of cell fates in multicellular organisms. However, simple and flexible isolation of microregions of interest (mROIs) for spatial transcriptomics is still challenging. We present a laser-induced forward transfer (LIFT)-based method combined with a full-length mRNA-sequencing protocol (LIFT-seq) for profiling region-specific tissues. LIFT-seq demonstrated that mROIs from two adjacent sections could reliably and sensitively detect and display gene expression. In addition, LIFT-seq can identify region-specific mROIs in the mouse cortex and hippocampus. Finally, LIFT-seq identified marker genes in different layers of the cortex with very similar expression patterns. These genes were then validated using in situ hybridization (ISH) results. Therefore, LIFT-seq will be a valuable and efficient technique for profiling the spatial transcriptome in various tissues.

Single-cell omics: experimental workflow, data analyses and applications

Cells are the fundamental units of biological systems and exhibit unique development trajectories and molecular features. Our exploration of how the genomes orchestrate the formation and maintenance of each cell, and control the cellular phenotypes of various organismsis, is both captivating and intricate. Since the inception of the first single-cell RNA technology, technologies related to single-cell sequencing have experienced rapid advancements in recent years. These technologies have expanded horizontally to include single-cell genome, epigenome, proteome, and metabolome, while vertically, they have progressed to integrate multiple omics data and incorporate additional information such as spatial scRNA-seq and CRISPR screening. Single-cell omics represent a groundbreaking advancement in the biomedical field, offering profound insights into the understanding of complex diseases, including cancers. Here, we comprehensively summarize recent advances in single-cell omics technologies, with a specific focus on the methodology section. This overview aims to guide researchers in selecting appropriate methods for single-cell sequencing and related data analysis.

FIND-seq: high-throughput nucleic acid cytometry for rare single-cell transcriptomics

Rare cells have an important role in development and disease, and methods for isolating and studying cell subsets are therefore an essential part of biology research. Such methods traditionally rely on labeled antibodies targeted to cell surface proteins, but large public databases and sophisticated computational approaches increasingly define cell subsets on the basis of genomic, epigenomic and transcriptomic sequencing data. Methods for isolating cells on the basis of nucleic acid sequences powerfully complement these approaches by providing experimental access to cell subsets discovered in cell atlases, as well as those that cannot be otherwise isolated, including cells infected with pathogens, with specific DNA mutations or with unique transcriptional or splicing signatures. We recently developed a nucleic acid cytometry platform called 'focused interrogation of cells by nucleic acid detection and sequencing' (FIND-seq), capable of isolating rare cells on the basis of RNA or DNA markers, followed by bulk or single-cell transcriptomic analysis. This platform has previously been used to characterize the splicing-dependent activation of the transcription factor XBP1 in astrocytes and HIV persistence in memory CD4 T cells from people on long-term antiretroviral therapy. Here, we outline the molecular and microfluidic steps involved in performing FIND-seq, including protocol updates that allow detection and whole transcriptome sequencing of rare HIV-infected cells that harbor genetically intact virus genomes. FIND-seq requires knowledge of microfluidics, optics and molecular biology. We expect that FIND-seq, and this comprehensive protocol, will enable mechanistic studies of rare HIV+ cells, as well as other cell subsets that were previously difficult to recover and sequence.

Pan-cancer mapping of single CD8+ T cell profiles reveals a TCF1:CXCR6 axis regulating CD28 co-stimulation and anti-tumor immunity

CD8+ T cells must persist and function in diverse tumor microenvironments to exert their effects. Thus, understanding common underlying expression programs could better inform the next generation of immunotherapies. We apply a generalizable matrix factorization algorithm that recovers both shared and context-specific expression programs from diverse datasets to a single-cell RNA sequencing (scRNA-seq) compendium of 33,161 CD8+ T cells from 132 patients with seven human cancers. Our meta-single-cell analyses uncover a pan-cancer T cell dysfunction program that predicts clinical non-response to checkpoint blockade in melanoma and highlights CXCR6 as a pan-cancer marker of chronically activated T cells. Cxcr6 is trans-activated by AP-1 and repressed by TCF1. Using mouse models, we show that Cxcr6 deletion in CD8+ T cells increases apoptosis of PD1+TIM3+ cells, dampens CD28 signaling, and compromises tumor growth control. Our study uncovers a TCF1:CXCR6 axis that counterbalances PD1-mediated suppression of CD8+ cell responses and is essential for effective anti-tumor immunity.

TET3 Contributes to Exercise-Induced Functional Axon Regeneration and Visual Restoration

Axons have intrinsically poor regenerative capacity in the mature central nervous system (CNS), leading to permanent neurological impairments in individuals. There is growing evidence that exercise is a powerful physiological intervention that can obviously enhance cell rejuvenate capacity, but its molecular mechanisms that mediate the axonal regenerative benefits remain largely unclear. Using the eye as the CNS model, here it is first indicated that placing mice in an exercise stimulation environment induced DNA methylation patterns and transcriptomes of retinal ganglion cell, promoted axon regeneration after injury, and reversed vision loss in aged mice. These beneficial effects are dependent on the DNA demethylases TET3-mediated epigenetic effects, which increased the expression of genes associated with the regenerative growth programs, such as STAT3, Wnt5a, Klf6. Exercise training also shows with the improved mitochondrial and metabolic dysfunction in retinas and optic nerves via TET3. Collectively, these results suggested that the increased regenerative capacity induced by enhancing physical activity is mediated through epigenetic reprogramming in mouse model of optic nerve injury and in aged mouse. Understanding the molecular mechanism underlying exercise-dependent neuronal plasticity led to the identification of novel targets for ameliorating pathologies associated with etiologically diverse diseases.

Single-cell total-RNA profiling unveils regulatory hubs of transcription factors

Recent development of RNA velocity uses master equations to establish the kinetics of the life cycle of RNAs from unspliced RNA to spliced RNA (i.e., mature RNA) to degradation. To feed this kinetic analysis, simultaneous measurement of unspliced RNA and spliced RNA in single cells is greatly desired. However, the majority of single-cell RNA-seq chemistry primarily captures mature RNA species to measure gene expressions. Here, we develop a one-step total-RNA chemistry-based single-cell RNA-seq method: snapTotal-seq. We benchmark this method with multiple single-cell RNA-seq assays in their performance in kinetic analysis of cell cycle by RNA velocity. Next, with LASSO regression between transcription factors, we identify the critical regulatory hubs mediating the cell cycle dynamics. We also apply snapTotal-seq to profile the oncogene-induced senescence and identify the key regulatory hubs governing the entry of senescence. Furthermore, from the comparative analysis of unspliced RNA and spliced RNA, we identify a significant portion of genes whose expression changes occur in spliced RNA but not to the same degree in unspliced RNA, indicating these gene expression changes are mainly controlled by post-transcriptional regulation. Overall, we demonstrate that snapTotal-seq can provide enriched information about gene regulation, especially during the transition between cell states.

In situ Patch-seq analysis of microglia reveals a lack of stress genes as found in FACS-isolated microglia

We applied the patch-seq technique to harvest transcripts from individual microglial cells from cortex, hippocampus and corpus callosum of acute brain slices from adult mice. After recording membrane currents with the patch-clamp technique, the cytoplasm was collected via the pipette and underwent adapted SMART-seq2 preparation with subsequent sequencing. On average, 4138 genes were detected in 113 cells from hippocampus, corpus callosum and cortex, including microglia markers such as Tmem119, P2ry12 and Siglec-H. Comparing our dataset to previously published single cell mRNA sequencing data from FACS-isolated microglia indicated that two clusters of cells were absent in our patch-seq dataset. Pathway analysis of marker genes in FACS-specific clusters revealed association with microglial activation and stress response. This indicates that under normal conditions microglia in situ lack transcripts associated with a stress-response, and that the microglia-isolation procedure by mechanical dissociation and FACS triggers the expression of genes related to activation and stress.

Understanding testicular single cell transcriptional atlas: from developmental complications to male infertility

Spermatogenesis is a multi-step biological process where mitotically active diploid (2n) spermatogonia differentiate into haploid (n) spermatozoa via regulated meiotic programming. The alarming rise in male infertility has become a global concern during the past decade thereby demanding an extensive profiling of testicular gene expression. Advancements in Next-Generation Sequencing (NGS) technologies have revolutionized our empathy towards complex biological events including spermatogenesis. However, despite multiple attempts made in the past to reveal the testicular transcriptional signature(s) either with bulk tissues or at the single-cell, level, comprehensive reviews on testicular transcriptomics and associated disorders are limited. Notably, technologies explicating the genome-wide gene expression patterns during various stages of spermatogenic progression provide the dynamic molecular landscape of testicular transcription. Our review discusses the advantages of single-cell RNA-sequencing (Sc-RNA-seq) over bulk RNA-seq concerning testicular tissues. Additionally, we highlight the cellular heterogeneity, spatial transcriptomics, dynamic gene expression and cell-to-cell interactions with distinct cell populations within the testes including germ cells (Gc), Sertoli cells (Sc), Peritubular cells (PTc), Leydig cells (Lc), etc. Furthermore, we provide a summary of key finding of single-cell transcriptomic studies that have shed light on developmental mechanisms implicated in testicular disorders and male infertility. These insights emphasize the pivotal roles of Sc-RNA-seq in advancing our knowledge regarding testicular transcriptional landscape and may serve as a potential resource to formulate future clinical interventions for male reproductive health.

Depletion of placental brain-derived neurotrophic factor (BDNF) is attributed to premature ovarian insufficiency (POI) in mice offspring

INTRODUCTION

Premature ovarian insufficiency (POI) is one of the causes of female infertility. Unexplained POI is increasingly affecting women in their reproductive years. However, the etiology of POI is diverse and remains elusive. We and others have shown that brain-derived neurotrophic factor (BDNF) plays an important role in adult ovarian function. Here, we report on a novel role of BDNF in the Developmental Origins of POI.

METHODS

Placental BDNF knockout mice were created using CRISPR/CAS9. Homozygous knockout (cKO(HO)) mice didn't survive, while heterozygous knockout (cKO(HE)) mice did. BDNF reduction in cKO(HE) mice was confirmed via immunohistochemistry and Western blots. Ovaries were collected from cKO(HE) mice at various ages, analyzing ovarian metrics, FSH expression, and litter sizes. In one-month-old mice, oocyte numbers were assessed using super-ovulation, and oocyte gene expression was analyzed with smart RNAseq. Ovaries of P7 mice were studied with SEM, and gene expression was confirmed with RT-qPCR. Alkaline phosphatase staining at E11.5 and immunofluorescence for cyclinD1 assessed germ cell number and cell proliferation.

RESULTS

cKO(HE) mice had decreased ovarian function and litter size in adulthood. They were insensitive to ovulation induction drugs manifested by lower oocyte release after superovulation in one-month-old cKO(HE) mice. The transcriptome and SEM results indicate that mitochondria-mediated cell death or aging might occur in cKO(HE) ovaries. Decreased placental BDNF led to diminished primordial germ cell proliferation at E11.5 and ovarian reserve which may underlie POI in adulthood.

CONCLUSION

The current results showed decreased placental BDNF diminished primordial germ cell proliferation in female fetuses during pregnancy and POI in adulthood. Our findings can provide insights into understanding the underlying mechanisms of POI.

DDX20 is required for cell-cycle reentry of prospermatogonia and establishment of spermatogonial stem cell pool during testicular development in mice

RNA-binding proteins (RBPs), as key regulators of mRNA fate, are abundantly expressed in the testis. However, RBPs associated with human male infertility remain largely unknown. Through bioinformatic analyses, we identified 62 such RBPs, including an evolutionarily conserved RBP, DEAD-box helicase 20 (DDX20). Male germ-cell-specific inactivation of Ddx20 at E15.5 caused T1-propsermatogonia (T1-ProSG) to fail to reenter cell cycle during the first week of testicular development in mice. Consequently, neither the foundational spermatogonial stem cell (SSC) pool nor progenitor spermatogonia were ever formed in the knockout testes. Mechanistically, DDX20 functions to control the translation of its target mRNAs, many of which encode cell-cycle-related regulators, by interacting with key components of the translational machinery in prospermatogonia. Our data demonstrate a previously unreported function of DDX20 as a translational regulator of critical cell-cycle-related genes, which is essential for cell-cycle reentry of T1-ProSG and formation of the SSC pool.

Beaconet: A Reference-Free Method for Integrating Multiple Batches of Single-Cell Transcriptomic Data in Original Molecular Space

Integrating multiple single-cell datasets is essential for the comprehensive understanding of cell heterogeneity. Batch effect is the undesired systematic variations among technologies or experimental laboratories that distort biological signals and hinder the integration of single-cell datasets. However, existing methods typically rely on a selected dataset as a reference, leading to inconsistent integration performance using different references, or embed cells into uninterpretable low-dimensional feature space. To overcome these limitations, a reference-free method, Beaconet, for integrating multiple single-cell transcriptomic datasets in original molecular space by aligning the global distribution of each batch using an adversarial correction network is presented. Through extensive comparisons with 13 state-of-the-art methods, it is demonstrated that Beaconet can effectively remove batch effect while preserving biological variations and is superior to existing unsupervised methods using all possible references in overall performance. Furthermore, Beaconet performs integration in the original molecular feature space, enabling the characterization of cell types and downstream differential expression analysis directly using integrated data with gene-expression features. Additionally, when applying to large-scale atlas data integration, Beaconet shows notable advantages in both time- and space-efficiencies. In summary, Beaconet serves as an effective and efficient batch effect removal tool that can facilitate the integration of single-cell datasets in a reference-free and molecular feature-preserved mode.

New insights into oocyte cytoplasmic lattice-associated proteins

Oocyte maturation and preimplantation embryo development are critical to successful pregnancy outcomes and the correct establishment and maintenance of genomic imprinting. Thanks to novel technologies and omics studies in human patients and mouse models, the importance of the proteins associated with the cytoplasmic lattices (CPLs), highly abundant structures found in the cytoplasm of mammalian oocytes and preimplantation embryos, in the maternal to zygotic transition is becoming increasingly evident. This review highlights the recent discoveries on the role of these proteins in protein storage and other oocyte cytoplasmic processes, epigenetic reprogramming, and zygotic genome activation (ZGA). A better comprehension of these events may significantly improve clinical diagnosis and pave the way for targeted interventions aiming to correct or mitigate female fertility issues and genomic imprinting disorders.