PRKDC-Mediated NHEJ May Play a Crucial Role in Aneuploidy of Chromosome 8-Driven Progression of Ovarian Cancer.

High malignancy is a prominent characteristic of epithelial ovarian cancer (EOC), emphasizing the necessity for further elucidation of the potential mechanisms underlying cancer progression. Aneuploidy and copy number variation (CNV) partially contribute to the heightened malignancy observed in EOC; however, the precise features of aneuploidy and their underlying molecular patterns, as well as the relationship between CNV and aneuploidy in EOC, remain unclear. In this study, we employed single-cell sequencing data along with The Cancer Genome Atlas (TCGA) to investigate aneuploidy and CNV in EOC. The technique of fluorescence in situ hybridization (FISH) was employed using specific probes. The copy number variation within the genomic region of chromosome 8 (42754568-47889815) was assessed and utilized as a representative measure for the ploidy status of individual cells in chromosome 8. Differential expression analysis was performed between different subgroups based on chromosome 8 ploidy. Gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), protein-protein interaction (PPI), and hub-gene analyses were subsequently utilized to identify crucial genes involved. By classifying enriched tumor cells into distinct subtypes based on chromosome 8 ploidy combined with TCGA data integration, we identified key genes driving chromosome 8 aneuploidy in EOC, revealing that PRKDC gene involvement through the mediated non-homologous end-joining pathway may play a pivotal role in disease progression. Further validation through analysis of the GEO and TCGA database and survival assessment, considering both mRNA expression levels and CNV status of PRKDC, has confirmed its involvement in the progression of EOC. Further functional analysis revealed an upregulation of PRKDC in both ovarian EOC cells and tissues, with its expression showing a significant correlation with the extent of copy number variation (CNV) on chromosome 8. Taken together, CNV amplification and aneuploidy of chromosome 8 are important characteristics of EOC. PRKDC and the mediated NHEJ pathway may play a crucial role in driving aneuploidy on chromosome 8 during the progression of EOC.

Unveiling gene regulatory networks during cellular state transitions without linkage across time points.

Time-stamped cross-sectional data, which lack linkage across time points, are commonly generated in single-cell transcriptional profiling. Many previous methods for inferring gene regulatory networks (GRNs) driving cell-state transitions relied on constructing single-cell temporal ordering. Introducing COSLIR (COvariance restricted Sparse LInear Regression), we presented a direct approach to reconstructing GRNs that govern cell-state transitions, utilizing only the first and second moments of samples between two consecutive time points. Simulations validated COSLIR's perfect accuracy in the oracle case and demonstrated its robust performance in real-world scenarios. When applied to single-cell RT-PCR and RNAseq datasets in developmental biology, COSLIR competed favorably with existing methods. Notably, its running time remained nearly independent of the number of cells. Therefore, COSLIR emerges as a promising addition to GRN reconstruction methods under cell-state transitions, bypassing the single-cell temporal ordering to enhance accuracy and efficiency in single-cell transcriptional profiling.

Single-Cell Chromatin Accessibility Analysis Reveals the Epigenetic Basis and Signature Transcription Factors for the Molecular Subtypes of Colorectal Cancers.

Colorectal cancer is a highly heterogeneous disease, with well-characterized subtypes based on genome, DNA methylome, and transcriptome signatures. To chart the epigenetic landscape of colorectal cancers, we generated a high-quality single-cell chromatin accessibility atlas of epithelial cells for 29 patients. Abnormal chromatin states acquired in adenomas were largely retained in colorectal cancers, which were tightly accompanied by opposite changes of DNA methylation. Unsupervised analysis on malignant cells revealed two epigenetic subtypes, exactly matching the iCMS classification, and key iCMS-specific transcription factors (TFs) were identified, including HNF4A and PPARA for iCMS2 tumors and FOXA3 and MAFK for iCMS3 tumors. Notably, subtype-specific TFs bind to distinct target gene sets and contribute to both interpatient similarities and diversities for both chromatin accessibilities and RNA expressions. Moreover, we identified CpG-island methylator phenotypes and pinpointed chromatin state signatures and TF regulators for the CIMP-high subtype. Our work systematically revealed the epigenetic basis of the well-known iCMS and CIMP classifications of colorectal cancers. SIGNIFICANCE: Our work revealed the epigenetic basis of the well-known iCMS and CIMP classifications of colorectal cancers. Moreover, interpatient minor similarities and major diversities of chromatin accessibility signatures of TF target genes can faithfully explain the corresponding interpatient minor similarities and major diversities of RNA expression signatures of colorectal cancers, respectively. This article is featured in Selected Articles from This Issue, p. 897.

Simultaneous de novo calling and phasing of genetic variants at chromosome-scale using NanoStrand-seq.

The successful accomplishment of the first telomere-to-telomere human genome assembly, T2T-CHM13, marked a milestone in achieving completeness of the human reference genome. The upcoming era of genome study will focus on fully phased diploid genome assembly, with an emphasis on genetic differences between individual haplotypes. Most existing sequencing approaches only achieved localized haplotype phasing and relied on additional pedigree information for further whole-chromosome scale phasing. The short-read-based Strand-seq method is able to directly phase single nucleotide polymorphisms (SNPs) at whole-chromosome scale but falls short when it comes to phasing structural variations (SVs). To shed light on this issue, we developed a Nanopore sequencing platform-based Strand-seq approach, which we named NanoStrand-seq. This method allowed for de novo SNP calling with high precision (99.52%) and acheived a superior phasing accuracy (0.02% Hamming error rate) at whole-chromosome scale, a level of performance comparable to Strand-seq for haplotype phasing of the GM12878 genome. Importantly, we demonstrated that NanoStrand-seq can efficiently resolve the MHC locus, a highly polymorphic genomic region. Moreover, NanoStrand-seq enabled independent direct calling and phasing of deletions and insertions at whole-chromosome level; when applied to long genomic regions of SNP homozygosity, it outperformed the strategy that combined Strand-seq with bulk long-read sequencing. Finally, we showed that, like Strand-seq, NanoStrand-seq was also applicable to primary cultured cells. Together, here we provided a novel methodology that enabled interrogation of a full spectrum of haplotype-resolved SNPs and SVs at whole-chromosome scale, with broad applications for species with diploid or even potentially polypoid genomes.

A human fetal cerebellar map of the late second trimester reveals developmental molecular characteristics and abnormality in trisomy 21.

Our understanding of human fetal cerebellum development during the late second trimester, a critical period for the generation of astrocytes, oligodendrocytes, and unipolar brush cells (UBCs), remains limited. Here, we performed single-cell RNA sequencing (scRNA-seq) in human fetal cerebellum samples from gestational weeks (GWs) 18-25. We find that proliferating UBC progenitors distribute in the subventricular zone of the rhombic lip (RLSVZ) near white matter (WM), forming a layer structure. We also delineate two trajectories from astrogenic radial glia (ARGs) to Bergmann glial progenitors (BGPs) and recognize oligodendrogenic radial glia (ORGs) as one source of primitive oligodendrocyte progenitor cells (PriOPCs). Additionally, our scRNA-seq analysis of the trisomy 21 fetal cerebellum at this stage reveals abnormal upregulated genes in pathways such as the cell adhesion pathway and focal adhesion pathway, which potentially promote neuronal differentiation. Overall, our research provides valuable insights into normal and abnormal development of the human fetal cerebellum.

Mapping crossover events of mouse meiotic recombination by restriction fragment ligation-based Refresh-seq.

Single-cell whole-genome sequencing methods have undergone great improvements over the past decade. However, allele dropout, which means the inability to detect both alleles simultaneously in an individual diploid cell, largely restricts the application of these methods particularly for medical applications. Here, we develop a new single-cell whole-genome sequencing method based on third-generation sequencing (TGS) platform named Refresh-seq (restriction fragment ligation-based genome amplification and TGS). It is based on restriction endonuclease cutting and ligation strategy in which two alleles in an individual cell can be cut into equal fragments and tend to be amplified simultaneously. As a new single-cell long-read genome sequencing method, Refresh-seq features much lower allele dropout rate compared with SMOOTH-seq. Furthermore, we apply Refresh-seq to 688 sperm cells and 272 female haploid cells (secondary polar bodies and parthenogenetic oocytes) from F1 hybrid mice. We acquire high-resolution genetic map of mouse meiosis recombination at low sequencing depth and reveal the sexual dimorphism in meiotic crossovers. We also phase the structure variations (deletions and insertions) in sperm cells and female haploid cells with high precision. Refresh-seq shows great performance in screening aneuploid sperm cells and oocytes due to the low allele dropout rate and has great potential for medical applications such as preimplantation genetic diagnosis.

Single-cell analysis reveals insights into epithelial abnormalities in ovarian endometriosis.

Ovarian endometriosis is characterized by the growth of endometrial tissue within the ovary, causing infertility and chronic pain. However, its pathophysiology remains unclear. Utilizing high-precision single-cell RNA sequencing, we profile the normal, eutopic, and ectopic endometrium from 34 individuals across proliferative and secretory phases. We observe an increased proportion of ciliated cells in both eutopic and ectopic endometrium, characterized by a diminished expression of estrogen sulfotransferase, which likely confers apoptosis resistance. After translocating to ectopic lesions, endometrial epithelium upregulates nicotinamide N-methyltransferase expression that inhibits apoptosis by promoting deacetylation and subsequent nuclear exclusion of transcription factor forkhead box protein O1, thereby leading to the downregulation of the apoptotic gene BIM. Moreover, epithelial cells in ectopic lesions elevate HLA class II complex expression, which stimulates CD4+ T cells and consequently contributes to chronic inflammation. Altogether, our study provides a comprehensive atlas of ovarian endometriosis and highlights potential therapeutic targets for modulating apoptosis and inflammation.

Integrative analysis of transcriptome, DNA methylome and chromatin accessibility reveals candidate therapeutic targets in hypertrophic cardiomyopathy.

Hypertrophic cardiomyopathy (HCM) is the most common inherited heart disease and is characterized by primary left ventricular hypertrophy usually caused by mutations in sarcomere genes. The mechanism underlying cardiac remodeling in HCM remains incompletely understood. An investigation of HCM through integrative analysis at multi-omics levels will be helpful for treating HCM. DNA methylation and chromatin accessibility, as well as gene expression, were assessed by nucleosome occupancy and methylome sequencing (NOMe-seq) and RNA-seq, respectively, using the cardiac tissues of HCM patients. Compared with those of the controls, the transcriptome, DNA methylome and chromatin accessibility of the HCM myocardium showed multifaceted differences. At the transcriptome level, HCM hearts returned to the fetal gene program through decreased sarcomeric and metabolic gene expression and increased extracellular matrix gene expression. In the DNA methylome, hypermethylated and hypomethylated differentially methylated regions (DMRs) were identified in HCM. At the chromatin accessibility level, HCM hearts showed changes in different genome elements. Several transcription factors (TFs), including SP1 and EGR1, exhibited a fetal-like pattern of binding motifs in nucleosome-depleted regions (NDRs) in HCM. In particular, the inhibition of SP1 or EGR1 in an HCM mouse model harboring sarcomere mutations markedly alleviated the HCM phenotype of the mutant mice and reversed fetal gene reprogramming. Overall, this study not only provides a high precision multi-omics map of HCM heart tissue but also sheds light on the therapeutic strategy by intervening the fetal gene reprogramming in HCM.

A single-cell characterised signature integrating heterogeneity and microenvironment of lung adenocarcinoma for prognostic stratification.

BACKGROUND: The high heterogeneity of tumour and the complexity of tumour microenvironment (TME) greatly impacted the tumour development and the prognosis of cancer in the era of immunotherapy. In this study, we aimed to portray the single cell-characterised landscape of lung adenocarcinoma (LUAD), and develop an integrated signature incorporating both tumour heterogeneity and TME for prognosis stratification. METHODS: Single-cell tagged reverse transcription sequencing (STRT-seq) was performed on tumour tissues and matched normal tissues from 14 patients with LUAD for immune landscape depiction and candidate key genes selection for signature construction. Kaplan-Meier survival analyses and in-vitro cell experiments were conducted to confirm the gene functions. The transcriptomic profile of 1949 patients from 11 independent cohorts including nine public datasets and two in-house cohorts were obtained for validation. FINDINGS: We selected 11 key genes closely related to cell-to-cell interaction, tumour development, T cell phenotype transformation, and Ma/Mo cell distribution, including HLA-DPB1, FAM83A, ITGB4, OAS1, FHL2, S100P, FSCN1, SFTPD, SPP1, DBH-AS1, CST3, and established an integrated 11-gene signature, stratifying patients to High-Score or Low-Score group for better or worse prognosis. Moreover, the prognostically-predictive potency of the signature was validated by 11 independent cohorts, and the immunotherapeutic predictive potency was also validated by our in-house cohort treated by immunotherapy. Additionally, the in-vitro cell experiments and drug sensitivity prediction further confirmed the gene function and generalizability of this signature across the entire RNA profile spectrum. INTERPRETATION: This single cell-characterised 11-gene signature might offer insights for prognosis stratification and potential guidance for treatment selection. FUNDING: Support for the study was provided by National key research and development project (2022YFC2505004, 2022YFC2505000 to Z.W. and J.W.), Beijing Natural Science Foundation (7242114 to J.X.), National Natural Science Foundation of China of China (82102886 to J.X., 81871889 and 82072586 to Z.W.), Beijing Nova Program (20220484119 to J.X.), NSFC general program (82272796 to J.W.), NSFC special program (82241229 to J.W.), CAMS Innovation Fund for Medical Sciences (2021-1-I2M-012, 2022-I2M-1-009 to Z.W. and J.W.), Beijing Natural Science Foundation (7212084 to Z.W.), CAMS Key lab of translational research on lung cancer (2018PT31035 to J.W.), Aiyou Foundation (KY201701 to J.W.). Medical Oncology Key Foundation of Cancer Hospital Chinese Academy of Medical Sciences (CICAMS-MOCP2022003 to J.X.).

DEAD-box helicase 17 (DDX17) protects cardiac function by promoting mitochondrial homeostasis in heart failure.

DEAD-box helicase 17 (DDX17) is a typical member of the DEAD-box family with transcriptional cofactor activity. Although DDX17 is abundantly expressed in the myocardium, its role in heart is not fully understood. We generated cardiomyocyte-specific Ddx17-knockout mice (Ddx17-cKO), cardiomyocyte-specific Ddx17 transgenic mice (Ddx17-Tg), and various models of cardiomyocyte injury and heart failure (HF). DDX17 is downregulated in the myocardium of mouse models of heart failure and cardiomyocyte injury. Cardiomyocyte-specific knockout of Ddx17 promotes autophagic flux blockage and cardiomyocyte apoptosis, leading to progressive cardiac dysfunction, maladaptive remodeling and progression to heart failure. Restoration of DDX17 expression in cardiomyocytes protects cardiac function under pathological conditions. Further studies showed that DDX17 can bind to the transcriptional repressor B-cell lymphoma 6 (BCL6) and inhibit the expression of dynamin-related protein 1 (DRP1). When DDX17 expression is reduced, transcriptional repression of BCL6 is attenuated, leading to increased DRP1 expression and mitochondrial fission, which in turn leads to impaired mitochondrial homeostasis and heart failure. We also investigated the correlation of DDX17 expression with cardiac function and DRP1 expression in myocardial biopsy samples from patients with heart failure. These findings suggest that DDX17 protects cardiac function by promoting mitochondrial homeostasis through the BCL6-DRP1 pathway in heart failure.