Dynamic human retinal pigment epithelium (RPE) and choroid architecture based on single-cell transcriptomic landscape analysis.

The retinal pigment epithelium (RPE) and choroid are located behind the human retina and have multiple functions in the human visual system. Knowledge of the RPE and choroid cells and their gene expression profiles are fundamental for understanding retinal disease mechanisms and therapeutic strategies. Here, we sequenced the RNA of about 0.3 million single cells from human RPE and choroids across two regions and seven ages, revealing regional and age differences within the human RPE and choroid. Cell-cell interactions highlight the broad connectivity networks between the RPE and different choroid cell types. Moreover, the transcription factors and their target genes change during aging. The coding of somatic variations increases during aging in the human RPE and choroid at the single-cell level. Moreover, we identified ELN as a candidate for improving RPE degeneration and choroidal structure during aging. The mapping of the molecular architecture of the human RPE and choroid improves our understanding of the human vision support system and offers potential insights into the intervention targets for retinal diseases.

Novel truncating variants in CTNNB1 cause familial exudative vitreoretinopathy.

BACKGROUND: Familial exudative vitreoretinopathy (FEVR) is an inheritable blinding disorder with clinical and genetic heterogeneity. Heterozygous variants in the CTNNB1 gene have been reported to cause FEVR. However, the pathogenic basis of CTNNB1-associated FEVR has not been fully explored. METHODS: Whole-exome sequencing was performed on the genomic DNA of probands. Dual-luciferase reporter assay, western blotting and co-immunoprecipitation were used to characterise the impacts of variants. Quantitative real-time PCR, EdU (5-ethynyl-2'-deoxyuridine) incorporation assay and immunocytochemistry were performed on the primary human retinal microvascular endothelial cells (HRECs) to investigate the effect of CTNNB1 depletion on the downstream genes involved in Norrin/ß-catenin signalling, cell proliferation and junctional integrity, respectively. Transendothelial electrical resistance assay was applied to measure endothelial permeability. Heterozygous endothelial-specific Ctnnb1-knockout mouse mice were generated to verify FEVR-like phenotypes in the retina. RESULTS: We identified two novel heterozygous variants (p.Leu103Ter and p.Val199LeufsTer11) and one previously reported heterozygous variant (p.His369ThrfsTer2) in the CTNNB1 gene. These variants caused truncation and degradation of ß-catenin that reduced Norrin/ß-catenin signalling activity. Additionally, knockdown (KD) of CTNNB1 in HRECs led to diminished mRNA levels of Norrin/ß-catenin targeted genes, reduced cell proliferation and compromised junctional integrity. The Cre-mediated heterozygous deletion of Ctnnb1 in mouse endothelial cells (ECs) resulted in FEVR-like phenotypes. Moreover, LiCl treatment partially rescued the defects in CTNNB1-KD HRECs and EC-specific Ctnnb1 heterozygous knockout mice. CONCLUSION: Our findings reinforced the current pathogenesis of Norrin/ß-catenin for FEVR and expanded the causative variant spectrum of CTNNB1 for the prenatal diagnosis and genetic counselling of FEVR.

Deep Sc-RNA sequencing decoding the molecular dynamic architecture of the human retina.

The human retina serves as a light detector and signals transmission tissue. Advanced insights into retinal disease mechanisms and therapeutic strategies require a deep understanding of healthy retina molecular events. Here, we sequenced the mRNA of over 0.6 million single cells from human retinas across six regions at nine different ages. Sixty cell sub-types have been identified from the human mature retinas with unique markers. We revealed regional and age differences of gene expression profiles within the human retina. Cell-cell interaction analysis indicated a rich synaptic connection within the retinal cells. Gene expression regulon analysis revealed the specific expression of transcription factors and their regulated genes in human retina cell types. Some of the gene's expression, such as DKK3, are elevated in aged retinas. A further functional investigation suggested that over expression of DKK3 could impact mitochondrial stability. Overall, decoding the molecular dynamic architecture of the human retina improves our understanding of the vision system.

A correlation-based feature analysis of physical examination indicators can help predict the overall underlying health status using machine learning.

As a systematic investigation of the correlations between physical examination indicators (PEIs) is lacking, most PEIs are currently independently used for disease warning. This results in the general physical examination having limited diagnostic values. Here, we systematically analyzed the correlations in 221 PEIs between healthy and 34 unhealthy statuses in 803,614 individuals in China. Specifically, the study population included 711,928 healthy participants, 51,341 patients with hypertension, 12,878 patients with diabetes, and 34,997 patients with other unhealthy statuses. We found rich relevance between PEIs in the healthy physical status (7662 significant correlations, 31.5%). However, in the disease conditions, the PEI correlations changed. We focused on the difference in PEIs between healthy and 35 unhealthy physical statuses and found 1239 significant PEI differences, suggesting that they could be candidate disease markers. Finally, we established machine learning algorithms to predict health status using 15-16% of the PEIs through feature extraction, reaching a 66-99% accurate prediction, depending on the physical status. This new reference of the PEI correlation provides rich information for chronic disease diagnosis. The developed machine learning algorithms can fundamentally affect the practice of general physical examinations.

A genetic variant in IL-6 lowering its expression is protective for critical patients with COVID-19.

Critical coronavirus disease 2019 (COVID-19) is associated with high mortality and potential genetic factors have been reported to be involved in the development of critical COVID-19. We performed a genome-wide association study to identify the genetic factors responsible for developing critical COVID-19. 632 critical patients with COVID-19 and 3021 healthy controls from the Chinese population were recruited. First, we identified a genome-wide significant difference of IL-6 rs2069837 (p = 9.73 × 10-15, OR = 0.41) between 437 critical patients with COVID-19 and 2551 normal controls in the discovery cohort. When replicated these findings in a set of 195 patients with critical COVID-19 and 470 healthy controls, we detected significant association of rs2069837 with COVID-19 (p = 8.89 × 10-3, OR = 0.67). This variant surpassed the formal threshold for genome-wide significance (combined p = 4.64 × 10-16, OR = 0.49). Further analysis revealed that there was a significantly stronger expression of IL-6 in the serum from patients with critical COVID-19 than in that from patients with asymptomatic COVID-19. An in vitro assay showed that the A to G allele changes in rs2069837 within IL-6 obviously decreased the luciferase expression activity. When analyzing the effect of this variant on the IL-6 in the serum based on the rs2069837 genotype, we found that the A to G variation in rs2069837 decreased the expression of IL-6, especially in the male. Overall, we identified a genetic variant in IL-6 that protects against critical conditions with COVID-19 though decreasing IL-6 expression in the serum.

Sulfated galactoglucan impedes xenografted lung cancer cell growth by blocking angiogenesis via binding BMPRs.

Evidences propound tumor growth may be impeded by blocking angiogenesis. Before we showed that sulfated glucan or arabinogalactan might bind to BMP2 or its receptors to inhibit angiogenesis. Whether sulfated galactoglucan can target both BMPRIA and BMPRII to impede angiogenesis and tumor cells growth is still vague. Here, we prepare galactoglucan and its sulfated derivatives Sul-CDA-0.05. The sulfate groups substituted are at the C-6 of 1, 4-linked α-Glcp and 1, 4-linked α-Galp backbone and at the C-6 of branch chain T-linked α-Glcp. Sul-CDA-0.05 can inhibit angiogenesis in vitro and in vivo. Indeed, Sul-CDA-0.05 impedes xenografted A549 lung tumor cells growth. Mechanism study demonstrates that this polysaccharide may target both BMPRIA and BMPRII to block BMP/Smad/Id1 signaling and attenuate VEGF and its transcription factor. Our evidences suggest that Sul-CDA-0.05 may be a new drug candidate for anti-lung cancer therapy by targeting both BMPRIA and BMPRII.

Integrated analysis on transcriptome and behaviors defines HTT repeat-dependent network modules in Huntington's disease.

Huntington's disease (HD) is caused by a CAG repeat expansion in the huntingtin (HTT) gene. Knock-in mice carrying a CAG repeat-expanded Htt will develop HD phenotypes. Previous studies suggested dysregulated molecular networks in a CAG length genotype- and the age-dependent manner in brain tissues from knock-in mice carrying expanded Htt CAG repeats. Furthermore, a large-scale phenome analysis defined a behavioral signature for HD genotype in knock-in mice carrying expanded Htt CAG repeats. However, an integrated analysis correlating phenotype features with genotypes (CAG repeat expansions) was not conducted previously. In this study, we revealed the landscape of the behavioral features and gene expression correlations based on 445 mRNA samples and 445 microRNA samples, together with behavioral features (396 PhenoCube behaviors and 111 NeuroCube behaviors) in Htt CAG-knock-in mice. We identified 37 behavioral features that were significantly associated with CAG repeat length including the number of steps and hind limb stand duration. The behavioral features were associated with several gene coexpression groups involved in neuronal dysfunctions, which were also supported by the single-cell RNA sequencing data in the striatum and the spatial gene expression in the brain. We also identified 15 chemicals with significant responses for genes with enriched behavioral features, most of them are agonist or antagonist for dopamine receptors and serotonin receptors used for neurology/psychiatry. Our study provides further evidence that abnormal neuronal signal transduction in the striatum plays an important role in causing HD-related phenotypic behaviors and provided rich information for the further pharmacotherapeutic intervention possibility for HD.

Age-related macular degeneration: Epidemiology, genetics, pathophysiology, diagnosis, and targeted therapy.

Age-related macular degeneration (AMD) is a complex eye disorder and is the leading cause of incurable blindness worldwide in the elderly. Clinically, AMD initially affects the central area of retina known as the macula and it is classified as early stage to late stage (advanced AMD). The advanced AMD is classified into the nonexudative or atrophic form (dry AMD) and the exudative or neovascular form (wet AMD). More severe vision loss is typically associated with the wet form. Multiple genetic factors, lipid metabolism, oxidative stress and aging, play a role in the etiology of AMD. Dysregulation in genetic to AMD is established to 46%-71% of disease contribution, with CFH and ARMS2/HTRA1 to be the two most notable risk loci among the 103 identified AMD associated loci so far. Chronic cigarette smoking is the most proven consistently risk living habits for AMD. Deep learning algorithm has been developed based on image recognition to distinguish wet AMD and normal macula with high accuracy. Currently, anti-vascular endothelial growth factor (VEGF) therapy is highly effective at treating wet AMD. Several new generation AMD drugs and iPSC-derived RPE cell therapy are in the clinical trial stage and are promising to improve AMD treatment in the near future.

Single-cell transcriptome reveals diversity of Müller cells with different metabolic-mitochondrial signatures in normal and degenerated macula.

Müller cell is the most abundant glial cell in mammalian retina, supporting the functions of photoreceptors and other retinal neurons via maintaining environmental homeostasis. In response to injury and/or neuronal degeneration, Müller cells undergo morphological and functional alternations, known as reactive gliosis documented in multiple retinal diseases, including age-related macular degeneration (AMD), retinitis pigmentosa, diabetic retinopathy, and traumatic retinal detachment. But the functional consequences of Müller glia cell reactivation or even the regulatory networks of the retinal gliosis are still controversial. In this study, we reveal different subpopulations of Müller cells with distinct metabolic-mitochondrial signatures by integrating single cell transcriptomic data from Early AMD patients and healthy donors. Our results show that a portion of Müller cells exhibits low mitochondrial DNA (mtDNA) expressions, reduced protein synthesis, impaired homeostatic regulation, decreased proliferative ability but enhanced proangiogenic function. Interestingly, the major alternation of Müller cells in Early AMD retina is the change of subpopulation abundance, rather than generation of new subcluster. Transcription factor enrichment analysis further highlights the key regulators of metabolic-mitochondrial states of Müller glias in Early AMD patients especially. Our study demonstrates new characteristics of retinal gliosis associated with Early AMD and suggests the possibility to prevent degeneration by intervening mitochondrial functions of Müller cells.

Whole-Exome Sequencing Reveals Novel TSPAN12 Variants in Autosomal Dominant Familial Exudative Vitreoretinopathy.

Background: Familial exudative vitreoretinopathy (FEVR), a group of rare inherited retinal vascular disorders, is the major cause of vision loss in juveniles. At present, the diagnosis of FEVR remains difficult due to its clinical and genetic heterogeneities. Aims: To identify the causative genetic variants in two unrelated FEVR-affected families: one Indian family and one Chinese Han family. Materials and Methods: Five affected patients from two families were recruited for this study. Whole-exome sequencing was applied to the probands, and Sanger sequencing was performed for validation. Stringent whole-exome sequence data analyses were performed to evaluate all of the identified pathogenic variants. Results: Two novel variants in the TSPAN12 gene, were identified: a missense variant c.437 T > G (p.Leu146Arg); and a nonsense variant c.477 C > A (p.Cys159*). Both variants cosegregated with the disease in the investigated FEVR-affected families. Additionally, both variants inactivated the ability of TSPAN12 protein to enhance Norrin/ß-catenin signaling. Conclusion: This study expands the mutational spectrum of TSPAN12 for FEVR.

Whole-Exome Sequencing Identified DLG1 as a Candidate Gene for Familial Exudative Vitreoretinopathy.

Purpose: Familial exudative vitreoretinopathy (FEVR) is a blinding retinal vascular disease. Clinically, FEVR is characterized by incomplete vascularization of the peripheral retina and pathological neovascularization. Only about 50% of FEVR cases can be explained by known FEVR disease gene variations. This study aimed to identify novel genes associated with the FEVR phenotype and explore their pathogenic mechanisms. Materials and Methods: Exome sequencing analyses were conducted on one Chinese family with FEVR whose affected members did not exhibit pathogenic variants in the known FEVR genes (verified using Sanger sequencing analysis). Functions of the affected proteins were evaluated using reporter assays. Western blot analysis was used to detect mutant protein expression and the genes' pathogenic mechanisms. Results: A rare novel heterozygous variant in DLG1 (c.1792A>G; p.S598G) was identified. The amino acid residues surrounding the identified variant are highly conserved among vertebrates. A luciferase reporter assay revealed that the mutant DLG1 protein DLG1-S598G lost its ability to activate Wnt signaling. Moreover, a knockdown (KD) of DLG1 in human primary retinal endothelial cells impaired tube formation. Mechanistically, DLG1 KD led to a reduction in phosphorylated VEGFR2, an essential receptor for the angiogenic potency that signals the vascular endothelial growth factor molecule. Conclusions: The data reported here demonstrate that DLG1 is a novel candidate gene for FEVR.

Untargeted metabolomics for uncovering plasma biological markers of wet age-related macular degeneration.

Wet age-related macular degeneration (wAMD) causes central vision loss and represents a major health problem in elderly people. Here we have used untargeted metabolomics using UHPLC-MS to profile plasma from 127 patients with wAMD (67 choroidal neovascularization (CNV) and 60 polypoidal choroidal vasculopathy (PCV)) and 50 controls. A total of 545 biochemicals were detected. Among them, 17 metabolites presented difference between patients with wAMD and controls. Most of them were oxidized lipids (N=6, 35.29%). Comparing to controls, 28 and 18 differential metabolites were identified in patients with CNV and PCV, respectively. Two metabolites, hyodeoxycholic acid and L-tryptophanamide, were differently distributed between PCV and CNV. We first investigated the genetic association with metabolites in wet AMD (CFH rs800292 and HTRA1 rs10490924). We identified six differential metabolites between the GG and AA genotypes of CFH rs800292, five differential metabolites between the GG and AA genotypes of HTRA1 rs10490924, and four differential metabolites between the GG and GA genotypes of rs10490924. We selected four metabolites (cyclamic acid, hyodeoxycholic acid, L-tryptophanamide and O-phosphorylethanolamine) for in vitro experiments. Among them, cyclamic acid reduced the activity, inhibited the proliferation, increased the apoptosis and necrosis in human retinal pigment epithelial cells (HRPECs). L-tryptophanamide affected the proliferation, apoptosis and necrosis in HRPECs, and promoted the tube formation and migration in primary human retinal endothelial cells (HRECs). Hyodeoxycholic acid and O-phosphorylethanolamine inhibited the tube formation and migration in HRECs. The results suggested that differential metabolites have certain effects on wAMD pathogenesis-related HRPECs and HRECs.

Dynamic blood single-cell immune responses in patients with COVID-19.

The 2019 coronavirus disease (COVID-19) outbreak caused by the SARS-CoV-2 virus is an ongoing global health emergency. However, the virus' pathogenesis remains unclear, and there is no cure for the disease. We investigated the dynamic changes of blood immune response in patients with COVID-19 at different stages by using 5' gene expression, T cell receptor (TCR), and B cell receptors (BCR) V(D)J transcriptome analysis at a single-cell resolution. We obtained single-cell mRNA sequencing (scRNA-seq) data of 341,420 peripheral blood mononuclear cells (PBMCs) and 185,430 clonotypic T cells and 28,802 clonotypic B cells from 25 samples of 16 patients with COVID-19 for dynamic studies. In addition, we used three control samples. We found expansion of dendritic cells (DCs), CD14+ monocytes, and megakaryocytes progenitor cells (MP)/platelets and a reduction of naïve CD4+ T lymphocytes in patients with COVID-19, along with a significant decrease of CD8+ T lymphocytes, and natural killer cells (NKs) in patients in critical condition. The type I interferon (IFN-I), mitogen-activated protein kinase (MAPK), and ferroptosis pathways were activated while the disease was active, and recovered gradually after patient conditions improved. Consistent with this finding, the mRNA level of IFN-I signal-induced gene IFI27 was significantly increased in patients with COVID-19 compared with that of the controls in a validation cohort that included 38 patients and 35 controls. The concentration of interferon-α (IFN-α) in the serum of patients with COVID-19 increased significantly compared with that of the controls in an additional cohort of 215 patients with COVID-19 and 106 controls, further suggesting the important role of the IFN-I pathway in the immune response of COVID-19. TCR and BCR sequences analyses indicated that patients with COVID-19 developed specific immune responses against SARS-CoV-2 antigens. Our study reveals a dynamic landscape of human blood immune responses to SARS-CoV-2 infection, providing clues for therapeutic potentials in treating COVID-19.

Catenin α 1 mutations cause familial exudative vitreoretinopathy by overactivating Norrin/ß-catenin signaling.

Familial exudative vitreoretinopathy (FEVR) is a severe retinal vascular disease that causes blindness. FEVR has been linked to mutations in several genes associated with inactivation of the Norrin/ß-catenin signaling pathway, but these account for only approximately 50% of cases. We report that mutations in α-catenin (CTNNA1) cause FEVR by overactivating the ß-catenin pathway and disrupting cell adherens junctions. We identified 3 heterozygous mutations in CTNNA1 (p.F72S, p.R376Cfs*27, and p.P893L) by exome sequencing and further demonstrated that FEVR-associated mutations led to overactivation of Norrin/ß-catenin signaling as a result of impaired protein interactions within the cadherin-catenin complex. The clinical features of FEVR were reproduced in mice lacking Ctnna1 in vascular endothelial cells (ECs) or with overactivated ß-catenin signaling by an EC-specific gain-of-function allele of Ctnnb1. In isolated mouse lung ECs, both CTNNA1-P893L and F72S mutants failed to rescue either the disrupted F-actin arrangement or the VE-cadherin and CTNNB1 distribution. Moreover, we discovered that compound heterozygous Ctnna1 F72S and a deletion allele could cause a similar phenotype. Furthermore, in a FEVR family, we identified a mutation of LRP5, which activates Norrin/ß-catenin signaling, and the corresponding knockin mice exhibited a partial FEVR-like phenotype. Our study demonstrates that the precise regulation of ß-catenin activation is critical for retinal vascular development and provides new insights into the pathogenesis of FEVR.

Scleral HIF-1α is a prominent regulatory candidate for genetic and environmental interactions in human myopia pathogenesis.

BACKGROUND: Myopia is a good model for understanding the interaction between genetics and environmental stimuli. Here we dissect the biological processes affecting myopia progression. METHODS: Human Genetic Analyses: (1) gene set analysis (GSA) of new genome wide association study (GWAS) data for 593 individuals with high myopia (refraction ≤ -6 diopters [D]); (2) over-representation analysis (ORA) of 196 genes with de novo mutations, identified by whole genome sequencing of 45 high-myopia trio families, and (3) ORA of 284 previously reported myopia risk genes. Contributions of the enriched signaling pathways in mediating the genetic and environmental interactions during myopia development were investigated in vivo and in vitro. RESULTS: All three genetic analyses showed significant enrichment of four KEGG signaling pathways, including amphetamine addiction, extracellular matrix (ECM) receptor interaction, neuroactive ligand-receptor interaction, and regulation of actin cytoskeleton pathways. In individuals with extremely high myopia (refraction ≤ -10 D), the GSA of GWAS data revealed significant enrichment of the HIF-1α signaling pathway. Using human scleral fibroblasts, silencing the key nodal genes within protein-protein interaction networks for the enriched pathways antagonized the hypoxia-induced increase in myofibroblast transdifferentiation. In mice, scleral HIF-1α downregulation led to hyperopia, whereas upregulation resulted in myopia. In human subjects, near work, a risk factor for myopia, significantly decreased choroidal blood perfusion, which might cause scleral hypoxia. INTERPRETATION: Our study implicated the HIF-1α signaling pathway in promoting human myopia through mediating interactions between genetic and environmental factors. FUNDING: National Natural Science Foundation of China grants; Natural Science Foundation of Zhejiang Province.

A comprehensive map of disease networks and molecular drug discoveries for glaucoma.

Glaucoma is the leading cause of irreversible blindness worldwide. The molecular etiology of glaucoma is complex and unclear. At present, there are few drugs available for glaucoma treatment. The aim of the present study was to perform a systematic analysis of glaucoma candidate drugs/chemicals based on glaucoma genes, including genetic factors and differentially expressed (DE) genes. In total, 401 genes from the genetic databases and 1656 genes from the DE gene analysis were included in further analyses. In terms of glaucoma-related genetic factors, 54 pathways were significantly enriched (FDR < 0.05), and 96 pathways for DE genes were significantly enriched (FDR < 0.05). A search of the PheWAS database for diseases associated with glaucoma-related genes returned 1,289 diseases, and a search for diseases associated with DE glaucoma-related genes returned 1,356 diseases. Cardiovascular diseases, neurodegenerative diseases, cancer, and ophthalmic diseases were highly related to glaucoma genes. A search of the DGIdb, KEGG, and CLUE databases revealed a set of drugs/chemicals targeting glaucoma genes. A subsequent analysis of the electronic medical records (EMRs) of 136,128 patients treated in Sichuan Provincial People's Hospital for candidate drug usage and the onset of glaucoma revealed nine candidate drugs. Among these drugs, individuals treated with nicardipine had the lowest incidence of glaucoma. Taken together with the information from the drug databases, the 40 most likely candidate drugs for glaucoma treatment were highlighted. Based on these findings, we concluded that the molecular mechanism of glaucoma is complex and may be a reflection of systemic diseases. A set of ready-to-use candidate drugs targeting glaucoma genes may be developed for glaucoma clinical drug treatments. Our results provide a systematic interpretation of glaucoma genes, interactions with other systemic diseases, and candidate drugs/chemicals.

A pectin-like polysaccharide from Polygala tenuifolia inhibits pancreatic cancer cell growth in vitro and in vivo by inducing apoptosis and suppressing autophagy.

Pancreatic ductal adenocarcinoma (PDAC) is one of the most lethal malignant carcinomas, which is characterized by apoptosis- and autophagy-dependent tumorigenic growth. Autophagy constitutes a stress adaptation that suppresses apoptosis. To explore new leading compound against PDAC, a pectin-like polysaccharide named RP02-1, was purified from roots of Polygala tenuifolia. Bioactivity test showed that RP02-1 might inhibit pancreatic cancer cells growth in vitro and in vivo. RP02-1 could inhibit pancreatic cancer cell (AsPC-1 and BxPC-3) proliferation, migration and colony formation. Mechanism study suggested that RP02-1 induced pancreatic cancer cells apoptosis, which was detected by Bcl-2 down-regulation, Bax up-regulation and conversion from Caspase 3 to Cleaved Caspase 3. Interestingly, autophagy was suppressed by RP02-1 treatment concentration-dependently through affenuatingBeclin-1, ATG5 and LC3B expression in BxPC-3 cells. In addition, RP02-1 could inhibit autophagy induced by Pennogenin 3-O-beta-chacotrioside. However, RP02-1 had almost no toxicity both in vitro and in vivo. The above results suggested that RP02-1 might be a potential leading compound for new drug candidate development for human PDAC treatment via inducing apoptosis and against autophagy.