Integrated Microarray and RNAseq Transcriptomic Analysis of Retinal Pigment Epithelium/Choroid in Age-Related Macular Degeneration

Bioinformatics analysis for constructing a cellular senescence-related age-related macular degeneration diagnostic model and identifying relevant disease subtypes to guide treatment

Age-related macular degeneration (AMD) is a condition causing progressive central vision loss. Growing evidence suggests a link between cellular senescence and AMD. However, the exact mechanism by which cellular senescence leads to AMD remains unclear. Employing machine learning, we established an AMD diagnostic model. Through unsupervised clustering, two distinct AMD subtypes were identified. GO, KEGG, and GSVA analyses explored the diverse biological functions associated with the two subtypes. By WGCNA, we constructed a coexpression network of differential genes between the subtypes, revealing the regulatory role of hub genes at the level of transcription factors and miRNAs. We identified 5 genes associated with inflammation for the construction of the AMD diagnostic model. Additionally, we observed that the level of cellular senescence and pathways related to programmed cell death (PCD), such as ferroptosis, necroptosis, and pyroptosis, exhibited higher expression levels in subtype B than A. Immune microenvironments also differed between the subtypes, indicating potentially distinct pathogenic mechanisms and therapeutic targets. In summary, by leveraging cellular senescence-associated gene expression, we developed an AMD diagnostic model. Furthermore, we identified two subtypes with varying expression patterns of senescence genes, revealing their differential roles in programmed cell death, disease progression, and immune microenvironments within AMD.

Roles of transmembrane protein 135 in mitochondrial and peroxisomal functions - implications for age-related retinal disease

Aging is the most significant risk factor for age-related diseases in general, which is true for age-related diseases in the eye including age-related macular degeneration (AMD). Therefore, in order to identify potential therapeutic targets for these diseases, it is crucial to understand the normal aging process and how its mis-regulation could cause age-related diseases at the molecular level. Recently, abnormal lipid metabolism has emerged as one major aspect of age-related symptoms in the retina. Animal models provide excellent means to identify and study factors that regulate lipid metabolism in relation to age-related symptoms. Central to this review is the role of transmembrane protein 135 (TMEM135) in the retina. TMEM135 was identified through the characterization of a mutant mouse strain exhibiting accelerated retinal aging and positional cloning of the responsible mutation within the gene, indicating the crucial role of TMEM135 in regulating the normal aging process in the retina. Over the past decade, the molecular functions of TMEM135 have been explored in various models and tissues, providing insights into the regulation of metabolism, particularly lipid metabolism, through its action in multiple organelles. Studies indicated that TMEM135 is a significant regulator of peroxisomes, mitochondria, and their interaction. Here, we provide an overview of the molecular functions of TMEM135 which is crucial for regulating mitochondria, peroxisomes, and lipids. The review also discusses the age-dependent phenotypes in mice with TMEM135 perturbations, emphasizing the importance of a balanced TMEM135 function for the health of the retina and other tissues including the heart, liver, and adipose tissue. Finally, we explore the potential roles of TMEM135 in human age-related retinal diseases, connecting its functions to the pathobiology of AMD.

A novel feature selection algorithm for identifying hub genes in lung cancer

Lung cancer, a life-threatening disease primarily affecting lung tissue, remains a significant contributor to mortality in both developed and developing nations. Accurate biomarker identification is imperative for effective cancer diagnosis and therapeutic strategies. This study introduces the Voting-Based Enhanced Binary Ebola Optimization Search Algorithm (VBEOSA), an innovative ensemble-based approach combining binary optimization and the Ebola optimization search algorithm. VBEOSA harnesses the collective power of the state-of-the-art classification models through soft voting. Moreover, our research applies VBEOSA to an extensive lung cancer gene expression dataset obtained from TCGA, following essential preprocessing steps including outlier detection and removal, data normalization, and filtration. VBEOSA aids in feature selection, leading to the discovery of key hub genes closely associated with lung cancer, validated through comprehensive protein-protein interaction analysis. Notably, our investigation reveals ten significant hub genes-ADRB2, ACTB, ARRB2, GNGT2, ADRB1, ACTG1, ACACA, ATP5A1, ADCY9, and ADRA1B-each demonstrating substantial involvement in the domain of lung cancer. Furthermore, our pathway analysis sheds light on the prominence of strategic pathways such as salivary secretion and the calcium signaling pathway, providing invaluable insights into the intricate molecular mechanisms underpinning lung cancer. We also utilize the weighted gene co-expression network analysis (WGCNA) method to identify gene modules exhibiting strong correlations with clinical attributes associated with lung cancer. Our findings underscore the efficacy of VBEOSA in feature selection and offer profound insights into the multifaceted molecular landscape of lung cancer. Finally, we are confident that this research has the potential to improve diagnostic capabilities and further enrich our understanding of the disease, thus setting the stage for future advancements in the clinical management of lung cancer. The VBEOSA source codes is publicly available at https://github.com/TEHNAN/VBEOSA-A-Novel-Feature-Selection-Algorithm-for-Identifying-hub-Genes-in-Lung-Cancer .

An Exploratory Data Analysis from Ovine and Bovine RNA-Seq Identifies Pathways and Key Genes Related to Cervical Dilatation

The present study developed a review and exploration of data in public and already validated repositories. The main objective was to identify the pathways involved in ruminants' cervical dilatation, which are conserved between cattle and sheep in the follicular and luteal phases of the reproductive cycle. In cattle, 1961 genes were more differentially expressed in the follicular phase and 1560 in the luteal phase. An amount of 24 genes were considered exclusively expressed from these. A total of 18 genes were in the follicular phase and 6 genes were in the luteal phase. In sheep, 2126 genes were more differentially expressed in the follicular phase and 2469 genes were more differentially expressed in the luteal phase. Hoxb genes were identified in both species and are correlated with the PI3K/Akt pathway. PI3K/Akt was also found in both cattle and sheep, appearing prominently in the follicular and luteal phases of both species. Our analyses have pointed out that the PI3K/Akt pathway and the Hoxb genes appear in prominence in modulating mechanisms that involve estrus alterations in the cervix. PI3K/Akt appears to be an important pathway in the cervical relaxation process.

Circular RNA expression and the competitive endogenous RNA network in pathological, age-related macular degeneration events: A cross-platform normalization study

Age-related macular degeneration (AMD) causes irreversible blindness in people aged over 50 worldwide. The dysfunction of the retinal pigment epithelium is the primary cause of atrophic AMD. In the current study, we used the ComBat and Training Distribution Matching method to integrate data obtained from the Gene Expression Omnibus database. We analyzed the integrated sequencing data by the Gene Set Enrichment Analysis. Peroxisome and tumor necrosis factor-α (TNF-α) signaling and nuclear factor kappa B (NF-κB) were among the top 10 pathways, and thus we selected them to construct AMD cell models to identify differentially expressed circular RNAs (circRNAs). We then constructed a competing endogenous RNA network, which is related to differentially expressed circRNAs. This network included seven circRNAs, 15 microRNAs, and 82 mRNAs. The Kyoto Encyclopedia of Genes and Genomes analysis of mRNAs in this network showed that the hypoxia-inducible factor-1 (HIF-1) signaling pathway was a common downstream event. The results of the current study may provide insights into the pathological processes of atrophic AMD.

Analysis of Wild Type and Variant B Cystatin C Interactome in Retinal Pigment Epithelium Cells Reveals Variant B Interacting Mitochondrial Proteins

Cystatin C, a secreted cysteine protease inhibitor, is abundantly expressed in retinal pigment epithelium (RPE) cells. A mutation in the protein's leader sequence, corresponding to formation of an alternate variant B protein, has been linked with an increased risk for both age-related macular degeneration (AMD) and Alzheimer's disease (AD). Variant B cystatin C displays intracellular mistrafficking with partial mitochondrial association. We hypothesized that variant B cystatin C interacts with mitochondrial proteins and impacts mitochondrial function. We sought to determine how the interactome of the disease-related variant B cystatin C differs from that of the wild-type (WT) form. For this purpose, we expressed cystatin C Halo-tag fusion constructs in RPE cells to pull down proteins interacting with either the WT or variant B form, followed by identification and quantification by mass spectrometry. We identified a total of 28 interacting proteins, of which 8 were exclusively pulled down by variant B cystatin C. These included 18 kDa translocator protein (TSPO) and cytochrome B5 type B, both of which are localized to the mitochondrial outer membrane. Variant B cystatin C expression also affected RPE mitochondrial function with increased membrane potential and susceptibility to damage-induced ROS production. The findings help us to understand how variant B cystatin C differs functionally from the WT form and provide leads to RPE processes adversely affected by the variant B genotype.

Data mining and network analysis reveals C-X-C chemokine receptor type 5 is involved in the pathophysiology of age-related macular degeneration

Our previous studies found that the C-X-C motif chemokine receptor 5 (CXCR5) loss leads to retinal pigment epithelium (RPE) dysfunction and AMD pathogenesis. The current study aimed to characterize the G protein-coupled receptor (GPCR) structure of CXCR5 and analyze its interactions with AMD-related risk genes. The sequence alignments, homology model of CXCR5 and structural assessment analysis were performed. Data and text mining were then performed to identify AMD-related risk genes and their interaction with CXCR5 using statistical and mathematical algorithms. Sequence alignment and phylogenetic tree analysis revealed that human CXCR5 was highly similar (85.4839%) to the rabbit. The least similarity (33.871%) was found to be in zebrafish compared to the other species. The CXCR5 model structural assessment and secondary structure analysis exhibited an excellent model. Network analysis revealed that IL10, TNF, ICAM1, CXCL1, CXCL8, APP, TLR4, SELL, C3, IL17A and CCR2 were the most connected genes CXCR5. These findings suggest that CXCR5 signaling may regulate the biological function of RPE and modulate AMD pathophysiology via GPCR signaling and interacting with identified AMD risk genes. In summary, the data presented here provide novel and crucial insights into the molecular mechanisms of CXCR5 involvement in AMD.Communicated by Ramaswamy H. Sarma.

Intercellular communication analysis of the human retinal pigment epithelial and choroidal cells predicts pathways associated with aging, cellular senescence and age-related macular degeneration

The retinal pigment epithelium (RPE) and the choroid are ocular tissues with fundamental roles in supporting neuroretinal function. The pathogenesis of age-related macular degeneration (AMD), a leading cause of irreversible blindness for which aging is the highest risk factor is closely linked with progressive impairment of various functions of these tissues. Cellular senescence, marked by cell cycle arrest and secretion of proinflammatory factors, is known to be associated with aging and has been proposed as a potential driver of AMD. Here, we investigated the role played by intercellular communication in the RPE/choroid within the context of aging, senescence and AMD. We inferred cell–cell interactions in the RPE/choroid by applying CellChat and scDiffCom on a publicly available scRNA-seq dataset from three human donors with and without AMD. We identified age-regulated ligand and receptor genes by using limma on a separate publicly available bulk microarray dataset providing RPE/choroid samples at multiple time points. Cellular senescence was investigated by assigning a score to each cell and each sample of these scRNA-seq and microarray datasets, respectively, based on the expression of key signature genes determined by a previous senescence meta-analysis. We identified VEGF-, BMP-and tenascin-mediated pathways supporting some of the strongest cell–cell interactions between RPE cells, fibroblasts and choroidal endothelial cells and as strong intercellular communication pathways related to both aging and senescence. Their signaling strength was enhanced between subpopulations of cells having high senescence scores. Predominant ligands of these pathways were upregulated with age whereas predominant receptors were downregulated. Globally, we also observed that cells from AMD samples presented slightly bigger senescence scores than normal cells and that the senescence score positively correlated with age in bulk samples (R = 0.26, value of p < 0.01). Hence, our analysis provides novel information on RPE/choroid intercellular communication that gives insights into the connection between aging, senescence and AMD.

Absence of Gut Microbiota Is Associated with RPE/Choroid Transcriptomic Changes Related to Age-Related Macular Degeneration Pathobiology and Decreased Choroidal Neovascularization

Studies have begun to reveal significant connections between the gut microbiome and various retinal diseases, including age-related macular degeneration (AMD). As critical supporting tissues of the retina, the retinal pigment epithelium (RPE) and underlying choroid play a critical role in retinal homeostasis and degeneration. However, the relationship between the microbiome and RPE/choroid remains poorly understood, particularly in animal models of AMD. In order to better elucidate this role, we performed high-throughput RNA sequencing of RPE/choroid tissue in germ-free (GF) and specific pathogen-free (SPF) mice. Furthermore, utilizing a specialized laser-induced choroidal neovascularization (CNV) model that we developed, we compared CNV size and inflammatory response between GF and SPF mice. After correction of raw data, 660 differentially expressed genes (DEGs) were identified, including those involved in angiogenesis regulation, scavenger and cytokine receptor activity, and inflammatory response-all of which have been implicated in AMD pathogenesis. Among lasered mice, the GF group showed significantly decreased CNV lesion size and microglial infiltration around CNV compared to the SPF group. Together, these findings provide evidence for a potential gut-RPE/choroidal axis as well as a correlation with neovascular features of AMD.

Transcriptomic and proteomic retinal pigment epithelium signatures of age-related macular degeneration

There are currently no treatments for geographic atrophy, the advanced form of age-related macular degeneration. Hence, innovative studies are needed to model this condition and prevent or delay its progression. Induced pluripotent stem cells generated from patients with geographic atrophy and healthy individuals were differentiated to retinal pigment epithelium. Integrating transcriptional profiles of 127,659 retinal pigment epithelium cells generated from 43 individuals with geographic atrophy and 36 controls with genotype data, we identify 445 expression quantitative trait loci in cis that are asssociated with disease status and specific to retinal pigment epithelium subpopulations. Transcriptomics and proteomics approaches identify molecular pathways significantly upregulated in geographic atrophy, including in mitochondrial functions, metabolic pathways and extracellular cellular matrix reorganization. Five significant protein quantitative trait loci that regulate protein expression in the retinal pigment epithelium and in geographic atrophy are identified - two of which share variants with cis- expression quantitative trait loci, including proteins involved in mitochondrial biology and neurodegeneration. Investigation of mitochondrial metabolism confirms mitochondrial dysfunction as a core constitutive difference of the retinal pigment epithelium from patients with geographic atrophy. This study uncovers important differences in retinal pigment epithelium homeostasis associated with geographic atrophy.

A mutation in transmembrane protein 135 impairs lipid metabolism in mouse eyecups

Aging is a significant factor in the development of age-related diseases but how aging disrupts cellular homeostasis to cause age-related retinal disease is unknown. Here, we further our studies on transmembrane protein 135 (Tmem135), a gene involved in retinal aging, by examining the transcriptomic profiles of wild-type, heterozygous and homozygous Tmem135 mutant posterior eyecup samples through RNA sequencing (RNA-Seq). We found significant gene expression changes in both heterozygous and homozygous Tmem135 mutant mouse eyecups that correlate with visual function deficits. Further analysis revealed that expression of many genes involved in lipid metabolism are changed due to the Tmem135 mutation. Consistent with these changes, we found increased lipid accumulation in mutant Tmem135 eyecup samples. Since mutant Tmem135 mice have similar ocular pathologies as human age-related macular degeneration (AMD) eyes, we compared our homozygous Tmem135 mutant eyecup RNA-Seq dataset with transcriptomic datasets of human AMD donor eyes. We found similar changes in genes involved in lipid metabolism between the homozygous Tmem135 mutant eyecups and AMD donor eyes. Our study suggests that the Tmem135 mutation affects lipid metabolism as similarly observed in human AMD eyes, thus Tmem135 mutant mice can serve as a good model for the role of dysregulated lipid metabolism in AMD.

Long Non-coding RNAs Gabarapl2 and Chrnb2 Positively Regulate Inflammatory Signaling in a Mouse Model of Dry Eye

Purpose: To elucidate the expression profile and the potential role of long non-coding ribonucleic acids (RNAs) (lncRNAs) in a dry eye disease (DED) model.

Methods: A DED model was established in C57BL/6J mice with 0.2% benzalkonium chloride (BAC) twice a day for 14 days. The differentially expressed lncRNAs were detected by RNA-seq technology (Gene Expression Omnibus, GEO GSE186450) and the aberrantly expressed lncRNAs were further verified by RT-qPCR. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were conducted to predicate the related candidate genes and potential pathological pathways. Cells from a human corneal epithelial cell line (HCECs) were cultured under hyperosmolarity. The regulation of inflammatory factors by silencing potential targeted lncRNAs was verified in vitro in HCECs.

Results: In our study, a significant increase in corneal fluorescence staining and a reduction in tear production were observed in DED mice at all follow-ups compared with the controls, and the differences were increasing over time. In total, 2,649 upregulated and 704 downregulated lncRNAs were identified in DED mice. We selected six aberrantly expressed and most abundant lncRNAs and performed RT-qPCR using the samples for RNA-seq. Chrnb2, Gabarapl2, and Usp31 were thereby confirmed as the most significantly altered lncRNAs. Pathway analysis revealed that the neuroactive ligand–receptor interaction signaling pathway was the most enriched, followed by the calcium signaling pathway and cytokine–cytokine receptor interaction. Following treatment of Gabarapl2 siRNA and Chrnb2 siRNA, tumor necrosis factor-α (TNF-α), interleukin (IL)-1β, and IL-6 were significantly downregulated in the HCECs.

Conclusion: Our study suggests that Chrnb2 and Gabarapl2 may be involved in the inflammation response by regulating TNF-α, IL-1β, and IL-6 in DED. These candidate lncRNAs may be both potential biomarkers and therapeutic targets for DED.