scDMV: a zero-one inflated beta mixture model for DNA methylation variability with scBS-seq data.

MOTIVATION: The utilization of single-cell bisulfite sequencing (scBS-seq) methods allows for precise analysis of DNA methylation patterns at the individual cell level, enabling the identification of rare populations, revealing cell-specific epigenetic changes, and improving differential methylation analysis. Nonetheless, the presence of sparse data and an overabundance of zeros and ones, attributed to limited sequencing depth and coverage, frequently results in reduced precision accuracy during the process of differential methylation detection using scBS-seq. Consequently, there is a pressing demand for an innovative differential methylation analysis approach that effectively tackles these data characteristics and enhances recognition accuracy. RESULTS: We propose a novel beta mixture approach called scDMV for analyzing methylation differences in single-cell bisulfite sequencing data, which effectively handles excess zeros and ones and accommodates low-input sequencing. Our extensive simulation studies demonstrate that the scDMV approach outperforms several alternative methods in terms of sensitivity, precision, and controlling the false positive rate. Moreover, in real data applications, we observe that scDMV exhibits higher precision and sensitivity in identifying differentially methylated regions, even with low-input samples. In addition, scDMV reveals important information for GO enrichment analysis with single-cell whole-genome sequencing data that are often overlooked by other methods. AVAILABILITY AND IMPLEMENTATION: The scDMV method, along with a comprehensive tutorial, can be accessed as an R package on the following GitHub repository: https://github.com/PLX-m/scDMV.

Endogenous mutant Huntingtin alters the corticogenesis via lowering Golgi recruiting ARF1 in cortical organoid.

Pathogenic mutant huntingtin (mHTT) infiltrates the adult Huntington's disease (HD) brain and impairs fetal corticogenesis. However, most HD animal models rarely recapitulate neuroanatomical alterations in adult HD and developing brains. Thus, the human cortical organoid (hCO) is an alternative approach to decode mHTT pathogenesis precisely during human corticogenesis. Here, we replicated the altered corticogenesis in the HD fetal brain using HD patient-derived hCOs. Our HD-hCOs had pathological phenotypes, including deficient junctional complexes in the neural tubes, delayed postmitotic neuronal maturation, dysregulated fate specification of cortical neuron subtypes, and abnormalities in early HD subcortical projections during corticogenesis, revealing a causal link between impaired progenitor cells and chaotic cortical neuronal layering in the HD brain. We identified novel long, oriented, and enriched polyQ assemblies of HTTs that hold large flat Golgi stacks and scaffold clathrin+ vesicles in the neural tubes of hCOs. Flat Golgi stacks conjugated polyQ assemblies by ADP-ribosylation factor 1 (ARF1). Inhibiting ARF1 activation with Brefeldin A (BFA) disassociated polyQ assemblies from Golgi. PolyQ assembles with mHTT scaffolded fewer ARF1 and formed shorter polyQ assembles with fewer and shorter Golgi and clathrin vesicles in neural tubes of HD-hCOs compared with those in hCOs. Inhibiting the activation of ARF1 by BFA in healthy hCOs replicated impaired junctional complexes in the neural tubes. Together, endogenous polyQ assemblies with mHTT reduced the Golgi recruiting ARF1 in the neuroepithelium, impaired the Golgi structure and activities, and altered the corticogenesis in HD-hCO.

Mendelian randomization analysis reveals a causal effect of Streptococcus salivarius on diabetic retinopathy through regulating host fasting glucose.

Diabetic retinopathy (DR) is one of leading causes of vision loss in adults with increasing prevalence worldwide. Increasing evidence has emphasized the importance of gut microbiome in the aetiology and development of DR. However, the causal relationship between gut microbes and DR remains largely unknown. To investigate the causal associations of DR with gut microbes and DR risk factors, we employed two-sample Mendelian Randomization (MR) analyses to estimate the causal effects of 207 gut microbes on DR outcomes. Inputs for MR included Genome-wide Association Study (GWAS) summary statistics of 207 taxa of gut microbes (the Dutch Microbiome Project) and 21 risk factors for DR. The GWAS summary statistics data of DR was from the FinnGen Research Project. Data analysis was performed in May 2023. We identified eight bacterial taxa that exhibited significant causal associations with DR (FDR < 0.05). Among them, genus Collinsella and species Collinsella aerofaciens were associated with increased risk of DR, while the species Bacteroides faecis, Burkholderiales bacterium_1_1_47, Ruminococcus torques, Streptococcus salivarius, genus Burkholderiales_noname and family Burkholderiales_noname showed protective effects against DR. Notably, we found that the causal effect of species Streptococcus salivarius on DR was mediated through the level of host fasting glucose, a well-established risk factor for DR. Our results reveal that specific gut microbes may be causally linked to DR via mediating host metabolic risk factors, highlighting potential novel therapeutic or preventive targets for DR.

Death associated protein like 1 acts as a novel tumor suppressor in melanoma by increasing the stability of P21 protein.

Melanoma is a primary malignant tumor with high lethality, which occurs in the skin and eye tissues, while the molecular mechanisms of melanomagenesis remain largely unknown. Here, we show that death-associated protein-like 1 (DAPL1) expression is lower in melanoma tissues than in paracancerous tissues or nevus tissues, and Uveal melanoma patients with lower DAPL1 expression have a poorer survival rate than those with higher expression of DAPL1. Overexpression of DAPL1 inhibits proliferation of cultured melanoma cells, whereas knockdown of DAPL1 increases cell proliferation. Tumor transplantation experiment results also demonstrate that DAPL1 inhibits tumorigenesis of melanoma cells both in subretinal and subcutaneous tissues of nude mice in vivo. Finally, DAPL1 inhibits proliferation of melanoma cells by increasing the protein level of P21 via decreasing the ubiquitin mediated degradation of P21 and promoting its stability. Conversely, knockdown of P21 neutralizes the effects of inhibition of DAPL1 on melanoma cell proliferation and enhances the severity of melanoma tumorigenesis. These results suggest that DAPL1 is a novel melanoma tumor suppressor gene and thus a potential therapeutic target for melanoma.

Protocol for detecting rare and common genetic associations in whole-exome sequencing studies using MAGICpipeline.

Whole-exome sequencing (WES) is a major approach to uncovering gene-disease associations and pinpointing effector genes. Here, we present a protocol for estimating genetic associations of rare and common variants in large-scale case-control WES studies using MAGICpipeline, an open-access analysis pipeline. We describe steps for assessing gene-based rare-variant association analyses by incorporating multiple variant pathogenic annotations and statistical techniques. We then detail procedures for identifying disease-related modules and hub genes using weighted correlation network analysis, a systems biology approach. For complete details on the use and execution of this protocol, please refer to Su et al. (2023).1.

Histone 3 lysine 9 acetylation-specific reprogramming regulates esophageal squamous cell carcinoma progression and metastasis.

Dysregulation of histone acetylation is widely implicated in tumorigenesis, yet its specific roles in the progression and metastasis of esophageal squamous cell carcinoma (ESCC) remain unclear. Here, we profiled the genome-wide landscapes of H3K9ac for paired adjacent normal (Nor), primary ESCC (EC) and metastatic lymph node (LNC) esophageal tissues from three ESCC patients. Compared to H3K27ac, we identified a distinct epigenetic reprogramming specific to H3K9ac in EC and LNC samples relative to Nor samples. This H3K9ac-related reprogramming contributed to the transcriptomic aberration of targeting genes, which were functionally associated with tumorigenesis and metastasis. Notably, genes with gained H3K9ac signals in both primary and metastatic lymph node samples (common-gained gene) were significantly enriched in oncogenes. Single-cell RNA-seq analysis further revealed that the corresponding top 15 common-gained genes preferred to be enriched in mesenchymal cells with high metastatic potential. Additionally, in vitro experiment demonstrated that the removal of H3K9ac from the common-gained gene MSI1 significantly downregulated its transcription, resulting in deficiencies in ESCC cell proliferation and migration. Together, our findings revealed the distinct characteristics of H3K9ac in esophageal squamous cell carcinogenesis and metastasis, and highlighted the potential therapeutic avenue for intervening ESCC through epigenetic modulation via H3K9ac.

Towards Stem/Progenitor Cell-Based Therapies for Retinal Degeneration.

Retinal degeneration (RD) is a leading cause of blindness worldwide and includes conditions such as retinitis pigmentosa (RP), age-related macular degeneration (AMD), and Stargardt's disease (STGD). These diseases result in the permanent loss of vision due to the progressive and irreversible degeneration of retinal cells, including photoreceptors (PR) and the retinal pigment epithelium (RPE). The adult human retina has limited abilities to regenerate and repair itself, making it challenging to achieve complete self-replenishment and functional repair of retinal cells. Currently, there is no effective clinical treatment for RD. Stem cell therapy, which involves transplanting exogenous stem cells such as retinal progenitor cells (RPCs), embryonic stem cells (ESCs), induced pluripotent stem cells (iPSCs), and mesenchymal stem cells (MSCs), or activating endogenous stem cells like Müller Glia (MG) cells, holds great promise for regenerating and repairing retinal cells in the treatment of RD. Several preclinical and clinical studies have shown the potential of stem cell-based therapies for RD. However, the clinical translation of these therapies for the reconstruction of substantial vision still faces significant challenges. This review provides a comprehensive overview of stem/progenitor cell-based therapy strategies for RD, summarizes recent advances in preclinical studies and clinical trials, and highlights the major challenges in using stem/progenitor cell-based therapies for RD.

Screening Mutations of the Monogenic Syndromic High Myopia by Whole Exome Sequencing From MAGIC Project.

Purpose: This observational study aimed to identify mutations in monogenic syndromic high myopia (msHM) using data from reported samples (n = 9370) of the Myopia Associated Genetics and Intervention Consortium (MAGIC) project. Methods: The targeted panel containing 298 msHM-related genes was constructed and screening of clinically actionable variants was performed based on whole exome sequencing. Capillary sequencing was used to verify the identified gene mutations in the probands and perform segregation analysis with their relatives. Results: A total of 381 candidate variants in 84 genes and 85 eye diseases were found to contribute to msHM in 3.6% (335/9370) of patients with HM. Among them, the 22 genes with the most variations accounted for 62.7% of the diagnostic cases. In the genotype-phenotype association analysis, 60% (201/335) of suspected msHM cases were recalled and 25 patients (12.4%) received a definitive genetic diagnosis. Pathogenic variants were distributed in 18 msHM-related diseases, mainly involving retinal dystrophy genes (e.g. TRPM1, CACNA1F, and FZD4), connective tissue disease genes (e.g. FBN1 and COL2A1), corneal or lens development genes (HSF4, GJA8, and MIP), and other genes (TEK). The msHM gene mutation types were allocated to four categories: nonsense mutations (36%), missense mutations (36%), frameshift mutations (20%), and splice site mutations (8%). Conclusions: This study highlights the importance of thorough molecular subtyping of msHM to provide appropriate genetic counselling and multispecialty care for children and adolescents with HM.

Multimodal analysis of cfDNA methylomes for early detecting esophageal squamous cell carcinoma and precancerous lesions.

Detecting early-stage esophageal squamous cell carcinoma (ESCC) and precancerous lesions is critical for improving survival. Here, we conduct whole-genome bisulfite sequencing (WGBS) on 460 cfDNA samples from patients with non-metastatic ESCC or precancerous lesions and matched healthy controls. We develop an expanded multimodal analysis (EMMA) framework to simultaneously identify cfDNA methylation, copy number variants (CNVs), and fragmentation markers in cfDNA WGBS data. cfDNA methylation markers are the earliest and most sensitive, detectable in 70% of ESCCs and 50% of precancerous lesions, and associated with molecular subtypes and tumor microenvironments. CNVs and fragmentation features show high specificity but are linked to late-stage disease. EMMA significantly improves detection rates, increasing AUCs from 0.90 to 0.99, and detects 87% of ESCCs and 62% of precancerous lesions with >95% specificity in validation cohorts. Our findings demonstrate the potential of multimodal analysis of cfDNA methylome for early detection and monitoring of molecular characteristics in ESCC.

DeepGraFT: A novel semantic segmentation auxiliary ROI-based deep learning framework for effective fundus tessellation classification.

Fundus tessellation (FT) is a prevalent clinical feature associated with myopia and has implications in the development of myopic maculopathy, which causes irreversible visual impairment. Accurate classification of FT in color fundus photo can help predict the disease progression and prognosis. However, the lack of precise detection and classification tools has created an unmet medical need, underscoring the importance of exploring the clinical utility of FT. Thus, to address this gap, we introduce an automatic FT grading system (called DeepGraFT) using classification-and-segmentation co-decision models by deep learning. ConvNeXt, utilizing transfer learning from pretrained ImageNet weights, was employed for the classification algorithm, aligning with a region of interest based on the ETDRS grading system to boost performance. A segmentation model was developed to detect FT exits, complementing the classification for improved grading accuracy. The training set of DeepGraFT was from our in-house cohort (MAGIC), and the validation sets consisted of the rest part of in-house cohort and an independent public cohort (UK Biobank). DeepGraFT demonstrated a high performance in the training stage and achieved an impressive accuracy in validation phase (in-house cohort: 86.85 %; public cohort: 81.50 %). Furthermore, our findings demonstrated that DeepGraFT surpasses machine learning-based classification models in FT classification, achieving a 5.57 % increase in accuracy. Ablation analysis revealed that the introduced modules significantly enhanced classification effectiveness and elevated accuracy from 79.85 % to 86.85 %. Further analysis using the results provided by DeepGraFT unveiled a significant negative association between FT and spherical equivalent (SE) in the UK Biobank cohort. In conclusion, DeepGraFT accentuates potential benefits of the deep learning model in automating the grading of FT and allows for potential utility as a clinical-decision support tool for predicting progression of pathological myopia.

Exome-wide association study identifies KDELR3 mutations in extreme myopia.

Extreme myopia (EM), defined as a spherical equivalent (SE) ≤ -10.00 diopters (D), is one of the leading causes of sight impairment. Known EM-associated variants only explain limited risk and are inadequate for clinical decision-making. To discover risk genes, we performed a whole-exome sequencing (WES) on 449 EM individuals and 9606 controls. We find a significant excess of rare protein-truncating variants (PTVs) in EM cases, enriched in the retrograde vesicle-mediated transport pathway. Employing single-cell RNA-sequencing (scRNA-seq) and a single-cell polygenic burden score (scPBS), we pinpointed PI16 + /SFRP4+ fibroblasts as the most relevant cell type. We observed that KDELR3 is highly expressed in scleral fibroblast and involved in scleral extracellular matrix (ECM) organization. The zebrafish model revealed that kdelr3 downregulation leads to elongated ocular axial length and increased lens diameter. Together, our study provides insight into the genetics of EM in humans and highlights KDELR3's role in EM pathogenesis.

Augmentation of scleral glycolysis promotes myopia through histone lactylation.

Myopia is characterized of maladaptive increases in scleral fibroblast-to-myofibroblast transdifferentiation (FMT). Scleral hypoxia is a significant factor contributing to myopia, but how hypoxia induces myopia is poorly understood. Here, we showed that myopia in mice and guinea pigs was associated with hypoxia-induced increases in key glycolytic enzymes expression and lactate levels in the sclera. Promotion of scleral glycolysis or lactate production induced FMT and myopia; conversely, suppression of glycolysis or lactate production eliminated or inhibited FMT and myopia. Mechanistically, increasing scleral glycolysis-lactate levels promoted FMT and myopia via H3K18la, and this promoted Notch1 expression. Genetic analyses identified a significant enrichment of two genes encoding glycolytic enzymes, ENO2 and TPI1. Moreover, increasing sugar intake in guinea pigs not only induced myopia but also enhanced the response to myopia induction via the scleral glycolysis-lactate-histone lactylation pathway. Collectively, we suggest that scleral glycolysis contributes to myopia by promoting FMT via lactate-induced histone lactylation.