Single-cell transcriptome profiling highlights the importance of telocyte, kallikrein genes, and alternative splicing in mouse testes aging.

Advancing healthcare for elderly men requires a deeper understanding of testicular aging processes. In this study, we conducted transcriptomic profiling of 43,323 testicular single cells from young and old mice, shedding light on 1032 telocytes-an underexplored testicular cell type in previous research. Our study unveiled 916 age-related differentially expressed genes (age-DEGs), with telocytes emerging as the cell type harboring the highest count of age-DEGs. Of particular interest, four genes (Klk1b21, Klk1b22, Klk1b24, Klk1b27) from the Kallikrein family, specifically expressed in Leydig cells, displayed down-regulation in aged testes. Moreover, cell-type-level splicing analyses unveiled 1838 age-related alternative splicing (AS) events. While we confirmed the presence of more age-DEGs in somatic cells compared to germ cells, unexpectedly, more age-related AS events were identified in germ cells. Further experimental validation highlighted 4930555F03Rik, a non-coding RNA gene exhibiting significant age-related AS changes. Our study represents the first age-related single-cell transcriptomic investigation of testicular telocytes and Kallikrein genes in Leydig cells, as well as the first delineation of cell-type-level AS dynamics during testicular aging in mice.

PUF60 promotes cell cycle and lung cancer progression by regulating alternative splicing of CDC25C.

Alternative splicing (AS) has been implicated in cell cycle regulation and cancer, but the underlying mechanisms are poorly understood. The poly(U)-binding splicing factor 60 (PUF60) is essential for embryonic development and is overexpressed in multiple types of cancer. Here, we report that PUF60 promotes mitotic cell cycle and lung cancer progression by controlling AS of the cell division cycle 25C (CDC25C). Systematic analysis of splicing factors deregulated in lung adenocarcinoma (LUAD) identifies that elevated copy number and expression of PUF60 correlate with poor prognosis. PUF60 depletion inhibits LUAD cell-cycle G2/M transition, cell proliferation, and tumor development. Mechanistically, PUF60 knockdown leads to exon skipping enriched in mitotic cell cycle genes, including CDC25C. Exon 3 skipping in the full-length CDC25C results in nonsense-mediated mRNA decay and a decrease of CDC25C protein, thereby inhibiting cell proliferation. This study establishes PUF60 as a cell cycle regulator and an oncogenic splicing factor in lung cancer.

The unique interplay of mitochondrial oxidative phosphorylation (OXPHOS) and immunity and its potential implication for the sex- and age-related morbidity of severe COVID-19 patients.

Aged male patients are more vulnerable to severe or critical symptoms of COVID-19, but the underlying mechanism remains elusive. In this study, we analyzed previously published scRNA-seq data from a large cohort of COVID-19 patients, castrated and regenerated mice, and bulk RNA-seq of a RNAi library of 400 genes, and revealed that both immunity and OXPHOS displayed cell-type-, sex-, and age-related variation in the severe or critical COVID-19 patients during disease progression, with a more prominent increase in immunity and decrease in OXPHOS in myeloid cells in the males relative to the females (60-69 years old). Male severe or critical patients above 70 years old were an exception in that the compromised negative correlation between OXPHOS and immunity in these patients was associated with its disordered transcriptional regulation. Finally, the expression levels of OXPHOS and androgens were revealed to be positively correlated, and the responses of macrophages to android fluctuation were more striking than other types of detected immune cells in the castrated mice model. Therefore, the interplay of OXPHOS and immunity displayed a cell-type-specific, age-related, and sex-biased pattern, and the underlying potential regulatory role of the hormonal milieu should not be neglected.

Classifying Heart-Sound Signals Based on CNN Trained on MelSpectrum and Log-MelSpectrum Features.

The intelligent classification of heart-sound signals can assist clinicians in the rapid diagnosis of cardiovascular diseases. Mel-frequency cepstral coefficients (MelSpectrums) and log Mel-frequency cepstral coefficients (Log-MelSpectrums) based on a short-time Fourier transform (STFT) can represent the temporal and spectral structures of original heart-sound signals. Recently, various systems based on convolutional neural networks (CNNs) trained on the MelSpectrum and Log-MelSpectrum of segmental heart-sound frames that outperform systems using handcrafted features have been presented and classified heart-sound signals accurately. However, there is no a priori evidence of the best input representation for classifying heart sounds when using CNN models. Therefore, in this study, the MelSpectrum and Log-MelSpectrum features of heart-sound signals combined with a mathematical model of cardiac-sound acquisition were analysed theoretically. Both the experimental results and theoretical analysis demonstrated that the Log-MelSpectrum features can reduce the classification difference between domains and improve the performance of CNNs for heart-sound classification.

Epigenetic regulation of human-specific gene expression in the prefrontal cortex.

BACKGROUND: Changes in gene expression levels during brain development are thought to have played an important role in the evolution of human cognition. With the advent of high-throughput sequencing technologies, changes in brain developmental expression patterns, as well as human-specific brain gene expression, have been characterized. However, interpreting the origin of evolutionarily advanced cognition in human brains requires a deeper understanding of the regulation of gene expression, including the epigenomic context, along the primate genome. Here, we used chromatin immunoprecipitation sequencing (ChIP-seq) to measure the genome-wide profiles of histone H3 lysine 4 trimethylation (H3K4me3) and histone H3 lysine 27 acetylation (H3K27ac), both of which are associated with transcriptional activation in the prefrontal cortex of humans, chimpanzees, and rhesus macaques. RESULTS: We found a discrete functional association, in which H3K4me3HP gain was significantly associated with myelination assembly and signaling transmission, while H3K4me3HP loss played a vital role in synaptic activity. Moreover, H3K27acHP gain was enriched in interneuron and oligodendrocyte markers, and H3K27acHP loss was enriched in CA1 pyramidal neuron markers. Using strand-specific RNA sequencing (ssRNA-seq), we first demonstrated that approximately 7 and 2% of human-specific expressed genes were epigenetically marked by H3K4me3HP and H3K27acHP, respectively, providing robust support for causal involvement of histones in gene expression. We also revealed the co-activation role of epigenetic modification and transcription factors in human-specific transcriptome evolution. Mechanistically, histone-modifying enzymes at least partially contribute to an epigenetic disturbance among primates, especially for the H3K27ac epigenomic marker. In line with this, peaks enriched in the macaque lineage were found to be driven by upregulated acetyl enzymes. CONCLUSIONS: Our results comprehensively elucidated a causal species-specific gene-histone-enzyme landscape in the prefrontal cortex and highlighted the regulatory interaction that drove transcriptional activation.

Epididymis cell atlas in a patient with a sex development disorder and a novel NR5A1 gene mutation.

This study aims to characterize the cell atlas of the epididymis derived from a 46,XY disorders of sex development (DSD) patient with a novel heterozygous mutation of the nuclear receptor subfamily 5 group A member 1 (NR5A1) gene. Next-generation sequencing found a heterozygous c.124C>G mutation in NR5A1 that resulted in a p.Q42E missense mutation in the conserved DNA-binding domain of NR5A1. The patient demonstrated feminization of external genitalia and Tanner stage 1 breast development. The surgical procedure revealed a morphologically normal epididymis and vas deferens but a dysplastic testis. Microfluidic-based single-cell RNA sequencing (scRNA-seq) analysis found that the fibroblast cells were significantly increased (approximately 46.5%), whereas the number of main epididymal epithelial cells (approximately 9.2%), such as principal cells and basal cells, was dramatically decreased. Bioinformatics analysis of cell-cell communications and gene regulatory networks at the single-cell level inferred that epididymal epithelial cell loss and fibroblast occupation are associated with the epithelial-to-mesenchymal transition (EMT) process. The present study provides a cell atlas of the epididymis of a patient with 46,XY DSD and serves as an important resource for understanding the pathophysiology of DSD.

Profiling ATM regulated genes in Drosophila at physiological condition and after ionizing radiation.

BACKGROUND: ATM (ataxia-telangiectasia mutated) protein kinase is highly conserved in metazoan, and plays a critical role at DNA damage response, oxidative stress, metabolic stress, immunity, RNA biogenesis etc. Systemic profiling of ATM regulated genes, including protein-coding genes, miRNAs, and long non-coding RNAs, will greatly improve our understanding of ATM functions and its regulation.  RESULTS: 1) differentially expressed protein-coding genes, miRNAs, and long non-coding RNAs in atm mutated flies were identified at physiological condition and after X-ray irradiation. 2) functions of differentially expressed genes in atm mutated flies, regardless of protein-coding genes or non-coding RNAs, are closely related with metabolic process, immune response, DNA damage response or oxidative stress. 3) these phenomena are persistent after irradiation. 4) there is a cross-talk regulation towards miRNAs by ATM, E2f1, and p53 during development and after irradiation. 5) knock-out flies or knock-down flies of most irradiation-induced miRNAs were sensitive to ionizing radiation. CONCLUSIONS: We provide a valuable resource of protein-coding genes, miRNAs, and long non-coding RNAs, for understanding ATM functions and regulations. Our work provides the new evidence of inter-dependence among ATM-E2F1-p53 for the regulation of miRNAs.

Accurate tumor clonal structures require single-cell analysis.

Tumor clonal structure is closely related to future progression, which has been mainly investigated as mutation abundance clustering in bulk samples. With relatively limited studies at single-cell resolution, a systematic comparison of the two approaches is still lacking. Here, using bulk and single-cell mutational data from the liver and colorectal cancers, we checked whether co-mutations determined by single-cell analysis had corresponding bulk variant allele frequency (VAF) peaks. While bulk analysis suggested the absence of subclonal peaks and, possibly, neutral evolution in some cases, the single-cell analysis identified coexisting subclones. The overlaps of bulk VAF ranges for co-mutations from different subclones made it difficult to separate them. Complex subclonal structures and dynamic evolution could be hidden under the seemingly clonal neutral pattern at the bulk level, suggesting single-cell analysis is necessary to avoid underestimation of tumor heterogeneity.

BET bromodomain inhibitor JQ1 regulates spermatid development by changing chromatin conformation in mouse spermatogenesis.

As a BET bromodomain inhibitor, JQ1 has been proven have efficacy against a number of different cancers. In terms of male reproduction, JQ1 may be used as a new type of contraceptive, since JQ1 treatment in male mice could lead to germ cell defects and a decrease of sperm motility, moreover, this effect is reversible. However, the mechanism of JQ1 acting on gene regulation in spermatogenesis remains unclear. Here, we performed single-cell RNA sequencing (scRNA-seq) on mouse testes treated with JQ1 or vehicle control to determine the transcriptional regulatory function of JQ1 in spermatogenesis at the single cell resolution. We confirmed that JQ1 treatment could increase the numbers of somatic cells and spermatocytes and decrease the numbers of spermatid cells. Gene Ontology (GO) analysis demonstrated that differentially expressed genes which were down-regulated after JQ1 injection were mainly enriched in "DNA conformation change" biological process in early developmental germ cells and "spermatid development" biological process in spermatid cells. ATAC-seq data further confirmed that JQ1 injection could change the open state of chromatin. In addition, JQ1 could change the numbers of accessible meiotic DNA double-stranded break sites and the types of transcription factor motif that functioned in pachytene spermatocytes and round spermatids. The multi-omics analysis revealed that JQ1 had the ability to regulate gene transcription by changing chromatin conformation in mouse spermatogenesis, which would potentiate the availability of JQ1 in male contraceptive.

Single-cell RNA-sequencing of zebrafish hair cells reveals novel genes potentially involved in hearing loss.

Hair cells play key roles in hearing and balance, and hair cell loss would result in hearing loss or vestibular dysfunction. Cellular and molecular research in hair cell biology provides us a better understanding of hearing and deafness. Zebrafish, owing to their hair cell-enriched organs, have been widely applied in hair cell-related research worldwide. Similar to mammals, zebrafish have inner ear hair cells. In addition, they also have lateral line neuromast hair cells. These different types of hair cells vary in morphology and function. However, systematic analysis of their molecular characteristics remains lacking. In this study, we analyzed the GFP+ cells isolated from Tg(Brn3c:mGFP) larvae with GFP expression in all hair cells using single-cell RNA-sequencing (scRNA-seq). Three subtypes of hair cells, namely macula hair cell (MHC), crista hair cell (CHC), and neuromast hair cell (NHC), were characterized and validated by whole-mount in situ hybridization analysis of marker genes. The hair cell scRNA-seq data revealed hair cell-specific genes, including hearing loss genes that have been identified in humans and novel genes potentially involved in hair cell formation and function. Two novel genes were discovered to specifically function in NHCs and MHCs, corresponding to their specific expression in NHCs and MHCs. This study allows us to understand the specific genes in hair cell subpopulations of zebrafish, which will shed light on the genetics of both human vestibular and cochlear hair cell function.

Single Cell Transcriptome Sequencing of Zebrafish Testis Revealed Novel Spermatogenesis Marker Genes and Stronger Leydig-Germ Cell Paracrine Interactions.

Spermatogenesis in testis is an important process for sexual reproduction, and worldwide about 10-15 percent of couples suffer from infertility. It is of importance to study spermatogenesis at single cell level in both of human and model organisms. Currently, single-cell RNA sequencing technologies (scRNA-seq) had been extensively applied to the study of cellular components and its gene regulations in the testes of different species, including human, monkey, mouse, and fly, but not in zebrafish. Zebrafish was a widely used model organism in biology and had been extensively used for the study of spermatogenesis in the previous studies. Therefore, it is also important to profile the transcriptome of zebrafish testis at single cell level. In this study, the transcriptomes of 14, 315 single cells from adult male zebrafish testes were profiled by scRNA-seq, and 10 cell populations were revealed, including Leydig cell, Sertoli cell, spermatogonia cell (SPG), spermatocyte, and spermatids. Notably, thousands of cell-type specific novel marker genes were identified, including sumo3b for SPG, krt18a.1 for Sertoli cells, larp1b and edrf1 for spermatids, which were also validated by RNA in situ hybridization experiments. Interestingly, through Ligand-Receptor (LR) analyses, zebrafish Leydig cells demonstrated stronger paracrine influence on germ cells than Sertoli cells. Overall, this study could be an important resource for the study of spermatogenesis in zebrafish and might also facilitate the study of the genes associated with human infertility through using zebrafish as a model organism.

The repertoire of testicular extracellular vesicle cargoes and their involvement in inter-compartmental communication associated with spermatogenesis.

BACKGROUND: Spermatogenesis is regulated by a complex network of intercellular communication processes. Extracellular vesicles (EVs) are one of the important mediators in intercellular communication. Previous reports have demonstrated the involvement of EVs from the epididymis and prostate in sperm maturation and function. However, the presence of EVs in the testis and their potential involvement in spermatogenesis has not been explored. Here, we have established a testis dissociation protocol that allows the isolation and characterization of testicular EVs. RESULTS: We show that testicular EVs are specifically and efficiently taken up by somatic cells and germ cells, including the spermatozoa in the interstitial space and the seminiferous tubule compartments. We profiled the proteome of testicular EVs and probed the cell types that release them, revealing the potential contributions from the Leydig cells and testicular macrophages. Moreover, we sequenced the small RNA cargoes of testicular EVs and identified sets of small non-coding RNAs that were overlooked in the testis transcriptome. Selected miRNA candidates in testicular EVs were found in sperm RNA payload and demonstrated specific resistance towards ribonuclease A independent of the vesicle membrane. Small molecule inhibition of EV secretion perturbed spermatogenesis via inter-compartmental communication. CONCLUSIONS: Together, our study provides a valuable resource on the repertoire of cargoes carried by testicular EVs and uncovers a physiological function of testicular EVs in inter-compartmental communication associated to spermatogenesis.

Nanopore Sequencing and Hi-C Based De Novo Assembly of Trachidermus fasciatus Genome.

Trachidermus fasciatus is a roughskin sculpin fish widespread across the coastal areas of East Asia. Due to environmental destruction and overfishing, the population of this species is under threat. In order to protect this endangered species, it is important to have the genome sequenced. Reference genomes are essential for studying population genetics, domestic farming, and genetic resource protection. However, currently, no reference genome is available for Trachidermus fasciatus, and this has greatly hindered the research on this species. In this study, we integrated nanopore long-read sequencing, Illumina short-read sequencing, and Hi-C methods to thoroughly assemble the Trachidermus fasciatus genome. Our results provided a chromosome-level high-quality genome assembly with a predicted genome size of 542.6 Mbp (2n = 40) and a scaffold N50 of 24.9 Mbp. The BUSCO value for genome assembly completeness was higher than 96%, and the single-base accuracy was 99.997%. Based on EVM-StringTie genome annotation, a total of 19,147 protein-coding genes were identified, including 35,093 mRNA transcripts. In addition, a novel gene-finding strategy named RNR was introduced, and in total, 51 (82) novel genes (transcripts) were identified. Lastly, we present here the first reference genome for Trachidermus fasciatus; this sequence is expected to greatly facilitate future research on this species.

Deep Learning Methods for Heart Sounds Classification: A Systematic Review.

The automated classification of heart sounds plays a significant role in the diagnosis of cardiovascular diseases (CVDs). With the recent introduction of medical big data and artificial intelligence technology, there has been an increased focus on the development of deep learning approaches for heart sound classification. However, despite significant achievements in this field, there are still limitations due to insufficient data, inefficient training, and the unavailability of effective models. With the aim of improving the accuracy of heart sounds classification, an in-depth systematic review and an analysis of existing deep learning methods were performed in the present study, with an emphasis on the convolutional neural network (CNN) and recurrent neural network (RNN) methods developed over the last five years. This paper also discusses the challenges and expected future trends in the application of deep learning to heart sounds classification with the objective of providing an essential reference for further study.

Spatio-temporal landscape of mouse epididymal cells and specific mitochondria-rich segments defined by large-scale single-cell RNA-seq.

Spermatozoa acquire their fertilizing ability and forward motility during epididymal transit, suggesting the importance of the epididymis. Although the cell atlas of the epididymis was reported recently, the heterogeneity of the cells and the gene expression profile in the epididymal tube are still largely unknown. Considering single-cell RNA sequencing results, we thoroughly studied the cell composition, spatio-temporal differences in differentially expressed genes (DEGs) in epididymal segments and mitochondria throughout the epididymis with sufficient cell numbers. In total, 40,623 cells were detected and further clustered into 8 identified cell populations. Focused analyses revealed the subpopulations of principal cells, basal cells, clear/narrow cells, and halo/T cells. Notably, two subtypes of principal cells, the Prc7 and Prc8 subpopulations were enriched as stereocilia-like cells according to GO analysis. Further analysis demonstrated the spatially specific pattern of the DEGs in each cell cluster. Unexpectedly, the abundance of mitochondria and mitochondrial transcription (MT) was found to be higher in the corpus and cauda epididymis than in the caput epididymis by scRNA-seq, immunostaining, and qPCR validation. In addition, the spatio-temporal profile of the DEGs from the P42 and P56 epididymis, including transiting spermatozoa, was depicted. Overall, our study presented the single-cell transcriptome atlas of the mouse epididymis and revealed the novel distribution pattern of mitochondria and key genes that may be linked to sperm functionalities in the first wave and subsequent wave of sperm, providing a roadmap to be emulated in efforts to achieve sperm maturation regulation in the epididymis.

Systematic profiling of early regulators during tissue regeneration using zebrafish model.

Great progresses have been made in comprehension of tissue regeneration process. However, one of the central questions in regeneration research remains to be deciphered is what factors initiate regenerative process. In present study, we focused on systematic profiling of early regulators in tissue regeneration via high-throughput screening on zebrafish caudal fin model. Firstly, 53 GO-annotated regeneration-related genes, which were specifically activated upon fin amputation, were identified according to the transcriptomic analysis. Moreover, qRT-PCR analysis of a couple of randomly selected genes from the aforementioned gene list validated our sequencing results. These studies confirmed the reliability of transcriptome sequencing analysis. Fibroblast growth factor 20a (fgf20a) is a key initial factor in the regeneration of zebrafish. Through a gene expression correlation analysis, we discovered a collection of 70 genes correlating with fgf20a, whose expression increased promptly at 2 days post amputation (dpa) and went down to the basal level until the completion of fin regeneration. In addition, two genes, socs3b and nppc, were chosen to investigate their functions during the fin regeneration. Inhibition of either of those genes significantly delayed the regenerative process. Taken together, we provided a simple and effective time-saving strategy that may serve as a tool for identifying early regulators in regeneration and identified 71 genes as early regulators of fin regeneration.

The slc4a2b gene is required for hair cell development in zebrafish.

Hair cells (HCs) function as important sensory receptors that can detect movement in their immediate environment. HCs in the inner ear can sense acoustic signals, while in aquatic vertebrates HCs can also detect movements, vibrations, and pressure gradients in the surrounding water. Many genes are responsible for the development of HCs, and developmental defects in HCs can lead to hearing loss and other sensory dysfunctions. Here, we found that the solute carrier family 4, member 2b (slc4a2b) gene, which is a member of the anion-exchange family, is expressed in the otic vesicles and lateral line neuromasts in developing zebrafish embryos. An in silico analysis showed that the slc4a2b is evolutionarily conserved, and we found that loss of function of slc4a2b resulted in a decreased number of HCs in zebrafish neuromasts due to increased HC apoptosis. Taken together, we conclude that slc4a2b plays a critical role in the development of HCs in zebrafish.

Transcriptome Sequencing Unravels Potential Biomarkers at Different Stages of Cerebral Ischemic Stroke.

Ischemic stroke, which accounts for 87% of all strokes, constitutes the leading cause of morbidity and mortality in China. Although the genetics and epigenetics of stroke have been extensively investigated, few studies have examined their relationships at different stages of stroke. This study assessed the characteristics of transcriptome changes at different stages of ischemic stroke using a mouse model of transient middle cerebral artery occlusion (tMCAO) and bioinformatics analyses. Cerebral cortex tissues from tMCAO mice at days 1, 3, 7, 14, and 28 were removed for RNA-Seq and small RNA-Seq library construction, sequencing, and bioinformatics analysis. We identified differentially expressed (DE) genes and miRNAs and revealed an association of the up-regulated or down-regulated DEmiRNAs with the correspondingly altered DEgene targets at each time point. In addition, different biological pathways were activated at different time points; thus, three groups of miRNAs were verified that may represent potential clinical biomarkers corresponding to days 1, 3, and 7 after ischemic stroke. Notably, this represents the first functional association of some of these miRNAs with stroke, e.g., miR-2137, miR-874-5p, and miR-5099. Together, our findings lay the foundation for the transition from a single-point, single-drug stroke treatment approach to multiple-time-point multi-drug combination therapies.