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.

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.

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.

Group Lasso Regularized Deep Learning for Cancer Prognosis from Multi-Omics and Clinical Features.

Accurate prognosis of patients with cancer is important for the stratification of patients, the optimization of treatment strategies, and the design of clinical trials. Both clinical features and molecular data can be used for this purpose, for instance, to predict the survival of patients censored at specific time points. Multi-omics data, including genome-wide gene expression, methylation, protein expression, copy number alteration, and somatic mutation data, are becoming increasingly common in cancer studies. To harness the rich information in multi-omics data, we developed GDP (Group lass regularized Deep learning for cancer Prognosis), a computational tool for survival prediction using both clinical and multi-omics data. GDP integrated a deep learning framework and Cox proportional hazard model (CPH) together, and applied group lasso regularization to incorporate gene-level group prior knowledge into the model training process. We evaluated its performance in both simulated and real data from The Cancer Genome Atlas (TCGA) project. In simulated data, our results supported the importance of group prior information in the regularization of the model. Compared to the standard lasso regularization, we showed that group lasso achieved higher prediction accuracy when the group prior knowledge was provided. We also found that GDP performed better than CPH for complex survival data. Furthermore, analysis on real data demonstrated that GDP performed favorably against other methods in several cancers with large-scale omics data sets, such as glioblastoma multiforme, kidney renal clear cell carcinoma, and bladder urothelial carcinoma. In summary, we demonstrated that GDP is a powerful tool for prognosis of patients with cancer, especially when large-scale molecular features are available.