A computational pipeline for functional gene discovery

Many computational pipelines exist for the detection of differentially expressed genes. However, computational pipelines for functional gene detection rarely exist. We developed a new computational pipeline for functional gene identification from transcriptome profiling data. Key features of the pipeline include batch effect correction, clustering optimization by gap statistics, gene ontology analysis of clustered genes, and literature analysis for functional gene discovery. By leveraging this pipeline on RNA-seq datasets from two mouse retinal development studies, we identified 7 candidate genes involved in the formation of the photoreceptor outer segment. The expression of top three candidate genes (Pde8b, Laptm4b, and Nr1h4) in the outer segment of the developing mouse retina were experimentally validated by immunohistochemical analysis. This computational pipeline can accurately predict novel functional gene for a specific biological process, e.g., development of the outer segment and synapses of the photoreceptor cells in the mouse retina. This pipeline can also be useful to discover functional genes for other biological processes and in other organs and tissues.

Loss of CEP162 function at the primary cilium delays ciliogenesis and causes retinal ciliopathy in humans

Ciliopathies often comprise retinal degeneration since the photoreceptor outer segment is an adapted primary cilium. CEP162 is a distal end centriolar protein required for proper transition zone assembly during ciliogenesis and whose loss causes ciliopathy in zebrafish. CEP162 has so far not been implicated in human disease. Here, we identified a homozygous CEP162 frameshift variant, c.1935dupA (p.(E646R*5)), in retinitis pigmentosa patients from two unrelated Moroccan families, likely representing a founder allele. We found that even though mRNA levels were reduced, the truncated CEP162-E646R*5 protein was expressed and localized to the mitotic spindle during mitosis, but not at the basal body of the cilium. In CEP162 knockdown cells, expression of the truncated CEP162-E646R*5 protein is unable to restore ciliation indicating its loss of function at the cilium. In patient fibroblasts, cilia overcome the absence of CEP162 from the primary cilium by delaying ciliogenesis through the persistence of CP110 at the mother centriole. The patient fibroblasts are ultimately able to extend some abnormally long cilia that are missing key transition zone components. Defective transition zone formation likely disproportionately affects the long-living ciliary outer segment of photoreceptors resulting in retinal dystrophy. CEP162 is expressed in human retina, and we show that wild-type CEP162, but not truncated CEP162-E646R*5, specifically localizes to the distal end of centrioles of mouse photoreceptor cilia. Together, our genetic, cell-based, and in vivo modeling establish that CEP162 deficiency causes retinal ciliopathy in humans.

A Novel GUCA1A Variant Associated with Cone Dystrophy Alters cGMP Signaling in Photoreceptors by Strongly Interacting with and Hyperactivating Retinal Guanylate Cyclase

Guanylate cyclase-activating protein 1 (GCAP1), encoded by the GUCA1A gene, is a neuronal calcium sensor protein involved in shaping the photoresponse kinetics in cones and rods. GCAP1 accelerates or slows the cGMP synthesis operated by retinal guanylate cyclase (GC) based on the light-dependent levels of intracellular Ca2+, thereby ensuring a timely regulation of the phototransduction cascade. We found a novel variant of GUCA1A in a patient affected by autosomal dominant cone dystrophy (adCOD), leading to the Asn104His (N104H) amino acid substitution at the protein level. While biochemical analysis of the recombinant protein showed impaired Ca2+ sensitivity of the variant, structural properties investigated by circular dichroism and limited proteolysis excluded major structural rearrangements induced by the mutation. Analytical gel filtration profiles and dynamic light scattering were compatible with a dimeric protein both in the presence of Mg2+ alone and Mg2+ and Ca2+. Enzymatic assays showed that N104H-GCAP1 strongly interacts with the GC, with an affinity that doubles that of the WT. The doubled IC50 value of the novel variant (520 nM for N104H vs. 260 nM for the WT) is compatible with a constitutive activity of GC at physiological levels of Ca2+. The structural region at the interface with the GC may acquire enhanced flexibility under high Ca2+ conditions, as suggested by 2 μs molecular dynamics simulations. The altered interaction with GC would cause hyper-activity of the enzyme at both low and high Ca2+ levels, which would ultimately lead to toxic accumulation of cGMP and Ca2+ in the photoreceptor outer segment, thus triggering cell death.

Deletion of prominin-1 in mice results in disrupted photoreceptor outer segment protein homeostasis

AIM: To illustrate the underlying mechanism how prominin-1 (also known as Prom1) mutation contribute to progressive photoreceptor degeneration. METHODS: A CRISPR-mediated Prom1 knockout (Prom1-KO) mice model in the C57BL/6 was generated and the photoreceptor degeneration phenotypes by means of structural and functional tests were demonstrated. Immunohistochemistry and immunoblot analysis were performed to reveal the localization and quantity of related outer segment (OS) proteins. RESULTS: The Prom1-KO mice developed the photoreceptor degeneration phenotype including the decreased outer nuclear layer (ONL) thickness and compromised electroretinogram amplitude. Immunohistochemistry analysis revealed impaired trafficking of photoreceptor OS proteins. Immunoblot data demonstrated decreased photoreceptor OS proteins. CONCLUSION: Prom1 deprivation causes progressive photoreceptor degeneration. Prom1 is essential for maintaining normal trafficking and normal quantity of photoreceptor OS proteins. The new light is shed on the pathogenic mechanism underlying photoreceptor degeneration caused by Prom1 mutation.

A Computational Pipeline for Functional Gene Discovery

Many computational pipelines exist for the detection of differentially expressed genes. However, computational pipelines for functional gene detection are rarely exist. We developed a new computational pipeline for functional gene identification from transcriptome profiling data. Key features of the pipeline include clustering optimization by gap statistics, gene ontology analysis for each cluster, and literature analysis for functional gene discovery. By leveraging this pipeline on RNA-seq datasets of mouse retinal development studies, we identified 14 candidate genes involved in the formation of the photoreceptor outer segment. The expression of top three candidate genes (Pde8b, Laptm4b, and Nr1h4) in the outer segment of the developing mouse retina were experimentally validated by immunohistochemical analysis. This computational pipeline can accurately predict novel functional gene for a specific biological process, e.g., the outer segment development of the photoreceptor cells in the mouse retina. This pipeline is also applicable to functional gene discovery for any other biological processes and in any other organs and tissues.

Acute RhoA/Rho Kinase Inhibition is Sufficient to Restore Phagocytic Capacity to Retinal Pigment Epithelium Lacking the Engulfment Receptor MerTK

The diurnal phagocytosis of spent photoreceptor outer segment fragments (POS) by retinal pigment epithelial (RPE) cells is essential for visual function. POS internalization by RPE cells requires the assembly of F-actin phagocytic cups beneath surface-tethered POS and Mer tyrosine kinase (MerTK) signaling. The activation of the Rho family GTPase Rac1 is necessary for phagocytic cup formation, and Rac1 is activated normally in MerTK-deficient RPE. We show here that mutant RPE lacking MerTK and wild-type RPE deprived of MerTK ligand both fail to form phagocytic cups regardless of Rac1 activation. However, in wild-type RPE in vivo, a decrease in RhoA activity coincides with the daily phagocytosis burst, while RhoA activity in MerTK-deficient RPE is constant. Elevating RhoA activity blocks phagocytic cup formation and phagocytosis by wild-type RPE. Conversely, inhibiting RhoA effector Rho kinases (ROCKs) rescues both F-actin assembly and POS internalization of primary RPE if MerTK or its ligand are lacking. Most strikingly, acute ROCK inhibition is sufficient to induce the formation and acidification of endogenous POS phagosomes by MerTK-deficient RPE ex vivo. Altogether, RhoA pathway inactivation is a necessary and sufficient downstream effect of MerTK phagocytic signaling such that the acute manipulation of cytosolic ROCK activity suffices to restore phagocytic capacity to MerTK-deficient RPE.

Deficiency in Lyst function leads to accumulation of secreted proteases and predisposition to mechanic stress-induced retinal detachment

Chediak–Higashi syndrome , caused by mutations in the Lys osome T rafficking Regulator ( Lyst ) gene, is a recessive hypopigmentation disorder characterized by albinism, neuropathies, neurodegeneration , and defective immune response s , with e nlargement of lysosomes and lysosome-related organelles. Although recent studies have suggested that Lyst mutations impair the regulation of sizes of lysosome and lysosome-related organelle , the underlying pathogenic mechanism of Chediak–Higashi syndrome is still unclear. Here we show striking evidence that deficiency in LYST protein function leads to accumulation of photoreceptor outer segment phagosomes in retinal pigment epitheli al cells , and reduce s adhesion between photoreceptor outer segment and retinal pigment epithelial cells in a mouse model of Chediak–Higashi syndrome . In addition, we observe elevated levels of cathepsin s , matrix metallopeptidase ( MMP ) 3 and oxidative stress markers in the retinal pigment epitheli um of Lyst mutant s . Previous reports showed that impaired degradation of photoreceptor outer segment phagosomes causes elevated oxidative stress , which could consequently lead to increase s of cysteine cathepsins and MMPs in the extracellular matrix. T aken together , we conclude that the loss of LYST function causes accumulation of phagosomes in the retinal pigment epitheli um and elevation of several extracellular matrix -remodeling proteases through oxidative stress , which may, in turn, reduce retina l adhesion. Our work reveals previously unreported pathogenic events in the retinal pigment epitheli um caused by Lyst deficiency , which may place Chediak–Higashi syndrome patients at increased risk for retinal detachment . The same pathogenic events may be conserved in other professional phagocytic cells, such as macrophages in the immune system, contributing to overall Chediak–Higashi syndrome pathology.

The Bardet-Biedl protein Bbs1 controls photoreceptor outer segment protein and lipid composition

Primary cilia are key sensory organelles whose dysfunction leads to ciliopathy disorders such as Bardet-Biedl syndrome (BBS). Retinal degeneration is common in ciliopathies, since the outer segments (OSs) of photoreceptors are highly specialized primary cilia. BBS1, encoded by the most commonly mutated BBS-associated gene, is part of the BBSome protein complex. Using a new bbs1 zebrafish mutant, we show that retinal development and photoreceptor differentiation are unaffected by Bbs1-loss, supported by an initially unaffected transcriptome. Quantitative proteomics and lipidomics on isolated OSs show that Bbs1 is required for BBSome-entry into OSs and that Bbs1-loss leads to accumulation of membrane-associated proteins in OSs, with enrichment in proteins involved in lipid homeostasis. Disruption of the tightly regulated OS lipid composition with increased OS cholesterol content are paralleled by early functional visual deficits, which precede progressive OS morphological anomalies. Our findings identify a new role for Bbs1/BBSome in OS lipid homeostasis and suggest a new pathomechanism underlying retinal degeneration in BBS.