3' MACE RNA-sequencing allows for transcriptome profiling in human tissue samples after long-term storage

Transcriptional profiling specifies the pathogen-specific human host response to infectious keratitis

Purpose

Corneal infections are a leading cause of visual impairment and blindness worldwide. Here we applied high-resolution transcriptomic profiling to assess the general and pathogen-specific molecular and cellular mechanisms during human corneal infection.

Methods

Clinical diagnoses of herpes simplex virus (HSV) (n=5) and bacterial/fungal (n=5) keratitis were confirmed by histology. Healthy corneas (n=7) and keratoconus (n=4) samples served as controls. Formalin-fixed, paraffin-embedded (FFPE) human corneal specimens were analyzed using the 3' RNA sequencing method Massive Analysis of cDNA Ends (MACE RNA-seq). The cellular host response was investigated using comprehensive bioinformatic deconvolution (xCell and CYBERSORTx) analyses and by integration with published single cell RNA-seq data of the human cornea.

Results

Our analysis identified 216 and 561 genes, that were specifically overexpressed in viral or bacterial/fungal keratitis, respectively, and allowed to distinguish the two etiologies. The virus-specific host response was driven by adaptive immunity and associated molecular signaling pathways, whereas the bacterial/fungal-specific host response mainly involved innate immunity signaling pathways and cell types. We identified several genes and pathways involved in the host response to infectious keratitis, including CXCL9, CXCR3, and MMP9 for viral, and S100A8/A9, MMP9, and the IL17 pathway for bacterial/fungal keratitis.

Conclusions

High-resolution molecular profiling provides new insights into the human corneal host response to viral and bacterial/fungal infection. Pathogen-specific molecular profiles may provide the foundation for novel diagnostic biomarker and therapeutic approaches that target inflammation-induced damage to corneal host cells with the goal to improve the outcome of infectious keratitis.

A comprehensive assessment of exome capture methods for RNA sequencing of formalin-fixed and paraffin-embedded samples

RNA-Seq analysis of Formalin-Fixed and Paraffin-Embedded (FFPE) samples has emerged as a highly effective approach and is increasingly being used in clinical research and drug development. However, the processing and storage of FFPE samples are known to cause extensive degradation of RNAs, which limits the discovery of gene expression or gene fusion-based biomarkers using RNA sequencing, particularly methods reliant on Poly(A) enrichment. Recently, researchers have developed an exome targeted RNA-Seq methodology that utilizes biotinylated oligonucleotide probes to enrich RNA transcripts of interest, which could overcome these limitations. Nevertheless, the standardization of this experimental framework, including probe designs, sample multiplexing, sequencing read length, and bioinformatic pipelines, remains an essential requirement. In this study, we conducted a comprehensive comparison of three main commercially available exome capture kits and evaluated key experimental parameters, to provide the overview of the advantages and limitations associated with the selection of library preparation protocols and sequencing platforms. The results provide valuable insights into the best practices for obtaining high-quality data from FFPE samples.

Impact of consensus molecular subtypes on survival with and without adjuvant chemotherapy in muscle-invasive urothelial bladder cancer

AIMS

Adjuvant chemotherapy after radical cystectomy can reduce the risk of recurrence and death in advanced muscle-invasive urothelial bladder cancer (MIBC). Molecular subtypes have been shown to be associated with survival. However, their predictive value to guide treatment decisions is controversial and data to use subtypes as guidance for adjuvant chemotherapy is sparse. We aimed to assess survival rates based on MIBC consensus molecular subtypes with and without adjuvant chemotherapy.

METHODS

Gene expression profiles of 143 patients with MIBC undergoing radical cystectomy were determined from formalin-fixed, paraffin-embedded specimen to assign consensus molecular subtypes. Expression of programmed cell death ligand-1 (PD-L1) and immune cell infiltration were determined using multiplex immunofluorescence. Matched-pair analysis was performed to evaluate the effect of adjuvant chemotherapy on overall survival (OS) for molecular subtypes applying Kaplan-Meier and Cox regression survival analyses.

RESULTS

Samples were luminal papillary: 9.1% (n=13), luminal non-specified: 6.3% (n=9), luminal unstable: 4.9% (n=7), stroma-rich: 27.9% (n=40), basal/squamous (Ba/Sq): 48.9% (n=70) and neuroendocrine-like (NE-like): 2.8% (n=4). Ba/Sq tumours had the highest concentration of PD-L1+ tumour and immune cells. Patients with luminal subtypes had better OS than those with NE-like (HR 0.2, 95% CI 0.1 to 0.7, p<0.05) and Ba/Sq (HR 0.5, 95% CI 0.2 to 0.9, p<0.05). No survival benefit with adjuvant chemotherapy was observed for luminal tumours, whereas Ba/Sq had significantly improved survival rates with adjuvant chemotherapy. Retrospective design and sample size are the main limitations.

CONCLUSION

Consensus molecular subtypes can be used to stratify patients with MIBC. Luminal tumours have the best prognosis and less benefit when receiving adjuvant chemotherapy compared with Ba/Sq tumours.

Optimizing identification of consensus molecular subtypes in muscle-invasive bladder cancer: a comparison of two sequencing methods and gene sets using FFPE specimens

BACKGROUND

Molecular subtypes predict prognosis in muscle-invasive bladder cancer (MIBC) and are explored as predictive markers. To provide a common base for molecular subtyping and facilitate clinical applications, a consensus classification has been developed. However, methods to determine consensus molecular subtypes require validation, particularly when FFPE specimens are used. Here, we aimed to evaluate two gene expression analysis methods on FFPE samples and to compare reduced gene sets to classify tumors into molecular subtypes.

METHODS

RNA was isolated from FFPE blocks of 15 MIBC patients. Massive analysis of 3' cDNA ends (MACE) and the HTG transcriptome panel (HTP) were used to retrieve gene expression. We used normalized, log2-transformed data to call consensus and TCGA subtypes with the consensusMIBC package for R using all available genes, a 68-gene panel (ESSEN1), and a 48-gene panel (ESSEN2).

RESULTS

Fifteen MACE-samples and 14 HTP-samples were available for molecular subtyping. The 14 samples were classified as Ba/Sq in 7 (50%), LumP in 2 (14.3%), LumU in 1 (7.1%), LumNS in 1 (7.1%), stroma-rich in 2 (14.3%) and NE-like in 1 (7.1%) case based on MACE- or HTP-derived transcriptome data. Consensus subtypes were concordant in 71% (10/14) of cases when comparing MACE with HTP data. Four cases with aberrant subtypes had a stroma-rich molecular subtype with either method. The overlap of the molecular consensus subtypes with the reduced ESSEN1 and ESSEN2 panels were 86% and 100%, respectively, with HTP data and 86% with MACE data.

CONCLUSION

Determination of consensus molecular subtypes of MIBC from FFPE samples is feasible using various RNA sequencing methods. Inconsistent classification mainly involves the stroma-rich molecular subtype, which may be the consequence of sample heterogeneity with (stroma)-cell sampling bias and highlights the limitations of bulk RNA-based subclassification. Classification is still reliable when analysis is reduced to selected genes.

The profiles of miR-4510 expression level in breast cancer

MicroRNA that is abnormally produced in breast cells can disrupt biological processes, which can lead to cancer. This study aims to screen differentially expressed genes (DEGs) and ncRNAs (DEncRNAs) in the formalin-fixed paraffin-embedded (FFPE) tissues of breast cancer (BC) as compared with the normal adjacent tissues (NAT), and identify miR-4510 as a novel biomarker of BC. This study looked at differentially expressed genes (DEGs) using MACE-Seq and differentially expressed ncRNAs (DEncRNAs) using the small RNA-Seq. Real-time qPCR was used to determine the level of expression of miR-4510. In this study, MACE-Seq results showed that 26,795 genes, with a p-value < 0.05, were differentially expressed in BC paraffin tissues as compared with NAT. Small RNA-Seq results revealed that 1326 ncRNAs, with a p-value < 0.05, were differentially expressed. We confirmed that miR-4510 was significantly down-expressed (p-value = 0.001) by qRT-PCR in the paraffin tissue of 120 BC patients. Based on eleven computational prediction programs, TP53, TP53INP1, MMP11, and COL1A1 for the miR-4510 were identified as miR-4510 targets. The MACE-seq result showed that the gene of TP53 (p-value = 0.001) and TP53INP1 (p-value = 0.02) was significantly down-regulated, but the gene of MMP11 (p-value = 0.004) and COL1A1 (p-value = 0.0001) was significantly over-expressed in 20 paired specimens of the BC and NAT. We discovered that a single SNP inside the miR-4510 binding site occurred only in BC, in which Guanine (G) changed into Adenine (A). Two SNPs outside the miR-4510 binding site occurred, and Guanine (G) in both BC and NAT was changed into Thymine (T), as compared to the reference sequence (RefSeq). Overall, our results suggested that miR-4510 functions as a tumor suppressor in the BC. Mir-4510 may act as a tumor suppressor, however additional experimental data is needed to corroborate these assumptions and can be exploited as a biomarker for BC.

Local and systemic targets of the MtCLE35-SUNN pathway in the roots of Medicago truncatula

CLE (CLAVATA3/ENDOSPERM SURROUNDING REGION-related) peptides are systemic regulators of legume-rhizobium symbiosis that negatively control the number of nitrogen-fixing nodules. CLE peptides are produced in the root in response to rhizobia inoculation and/or nitrate treatment and are transported to the shoot where they are recognized by the CLV1-like (CLAVATA1-like) receptor kinase. As a result, a shoot-derived signaling pathway is activated that inhibits subsequent nodule development in the root. In Medicago truncatula, MtCLE35 is activated in response to rhizobia and nitrate treatment and the overexpression of this gene systemically inhibits nodulation. The inhibitory effect of MtCLE35 overexpression is dependent on the CLV1-like receptor kinase MtSUNN (SUPER NUMERIC NODULES), suggesting that MtSUNN could be involved in the reception of the MtCLE35 peptide. Yet little is known about the downstream genes regulated by a MtCLE35-activated response in the root. In order to identify genes whose expression levels could be regulated by the MtCLE35-MtSUNN pathway, we performed a MACE-Seq (Massive Analysis of cDNA Ends) transcriptomic analysis of MtCLE35-overexpressing roots. Among upregulated genes, the gene MtSUNN that encodes a putative receptor of MtCLE35 was detected. Moreover, we found that MtSUNN, as well as several other differentially expressed genes, were upregulated locally in MtCLE35-overexpressing roots whereas the MtTML1 and MtTML2 genes were upregulated systemically. Our data suggest that MtCLE35 has both local and systemic effects on target genes in the root.

Transcriptional Profiling Identifies Prognostic Gene Signatures for Conjunctival Extranodal Marginal Zone Lymphoma

This study characterizes the transcriptional profile and the cellular tumor microenvironment of conjunctival extranodal marginal zone lymphoma (EMZL) and identifies prognostically relevant biomarkers. Ten formalin-fixed and paraffin-embedded conjunctival EMZL and eight healthy conjunctival specimens were analyzed by Massive Analysis of cDNA Ends (MACE) RNA sequencing. The 3417 upregulated genes in conjunctival EMZL were involved in processes such as B cell proliferation and Rac protein signaling, whereas the 1188 downregulated genes contributed most significantly to oxidative phosphorylation and UV protection. The tumor microenvironment, as determined by deconvolution analysis, was mainly composed of multiple B cell subtypes which reflects the tumor's B cell lineage. However, several T cell types, including T helper 2 cells and regulatory T cells, as well as innate immune cell types, such as anti-inflammatory macrophages and plasmacytoid dendritic cells, were also strongly enriched in conjunctival EMZL. A 13-biomarker prognostic panel, including S100A8 and S100A9, classified ocular and extraocular tumor recurrence, exceeded prognostic accuracy of Ann Arbor and American Joint Committee on Cancer (AJCC) staging, and demonstrated prognostic value for patient survival in 21 different cancer types in a database of 12,332 tumor patients. These findings may lead to new options of targeted therapy and may improve prognostic prediction for conjunctival EMZL.

Web-based gene expression analysis-paving the way to decode healthy and diseased ocular tissue

BACKGROUND

Gene expression analysis using RNA sequencing has helped to improve the understanding of many diseases. Databases, such as the Gene Expression Omnibus database of the National Center for Biotechnology Information provide RNA sequencing raw data from various diseased tissue types but their analysis requires advanced bioinformatics skills. Therefore, specific ocular databases provide the transcriptional profiles of different ocular tissues and in addition enable intuitive web-based data analysis.

OBJECTIVE

The aim of this narrative review is to provide an overview of ocular transcriptome databases and to compare them with the Human Eye Transcriptome Atlas newly established in Freiburg.

METHODS

PubMed literature search.

RESULTS

A total of nine ocular transcriptome databases focusing on different aspects were identified. The iSyTE and Express platforms specialize in gene expression during lens and retinal development in mice, whereas retina.tigem.it, Eye in a Disk, and Spectacle focus on selected ocular tissues such as the retina. Spectacle, UCSC Cell Browser and Single Cell Portal allow intuitive exploration of single cell RNA sequencing data derived from retinal, choroid, cornea, iris, trabecular meshwork and sclera specimens. The microarray profiles of a variety of healthy ocular tissues are included in the Ocular Tissue Database. The Human Eye Transcriptome Atlas provides the largest collection of different ocular tissue types, contains the highest number of ocular diseases and is characterized by a high level of quality achieved by methodological consistency.

CONCLUSION

Ocular transcriptome databases provide comprehensive and intuitive insights into the transcriptional profiles of a variety of healthy and diseased ocular tissues. Thus, they improve our understanding of the underlying molecular mediators, support hypothesis generation and help in the search for new diagnostic and therapeutic targets for various ocular diseases.

[Next-generation sequencing in ophthalmology]

New methods for basic research are finding their way into ophthalmology and expand the options for research, diagnostics and treatment. This article centers on the study of gene activity in cells and tissues by next-generation sequencing (NGS). Following a brief introduction to the basic principles of NGS, this article focuses on transcriptome analysis by RNA sequencing using examples from ophthalmology. The RNA sequencing provides a comprehensive and unbiased overview of gene activity in cells and tissues and thus forms an important foundation for generating new testable scientific hypotheses. It thus contributes to the in-depth characterization of pathological changes and supports the development of new diagnostic and therapeutic approaches.

Transcriptional Profiling Provides New Insights into Organ Culture-Induced Changes in Human Donor Corneas

Corneal transplantation is one of the most common forms of tissue transplantation worldwide. Donor corneal tissue used in transplantation is provided by eye banks, which store the tissue in culture medium after procurement. To date, the effects of cell culture on human corneal tissue have not been fully elucidated. Using the 3' RNA sequencing method for massive analysis of cDNA ends (MACE), we show that cultivation of corneal tissue leads to significant changes in a variety of molecular processes in human corneal tissue that go well beyond aspects of previously known culture effects. Functionally grouped network analysis revealed nine major groups of biological processes that were affected by corneal organ culture, among them keratinization, hypoxia, and angiogenesis, with genes from each group being affected by culture time. A cell type deconvolution analysis revealed significant modulations of the corneal immune cell profile in a time dependent manner. The results suggest that current culture conditions should be further refined and that prolonged cultivation may be detrimental. Recently, we showed that MACE enables transcriptional profiling of formalin-fixed and paraffin-embedded (FFPE) conjunctival tissue with high accuracy even after more than 10 years of storage. Here we demonstrate that MACE provides comparable results for native and FFPE corneal tissue, confirming that the technology is suitable for transcriptome analysis of a wide range of archived diseased corneal samples stored in histological archives. Finally, our data underscore the feasibility of bioinformatics cell-type enrichment analysis in bulk RNA-seq data to profile immune cell composition in fixed and archived corneal tissue samples, for which RNA-seq analysis of individual cells is often not possible.

Transcriptional and Distributional Profiling of Microglia in Retinal Angiomatous Proliferation

Macular neovascularization type 3, formerly known as retinal angiomatous proliferation (RAP), is a hallmark of age-related macular degeneration and is associated with an accumulation of myeloid cells, such as microglia (MG) and infiltrating blood-derived macrophages (MAC). However, the contribution of MG and MAC to the myeloid cell pool at RAP sites and their exact functions remain unknown. In this study, we combined a microglia-specific reporter mouse line with a mouse model for RAP to identify the contribution of MG and MAC to myeloid cell accumulation at RAP and determined the transcriptional profile of MG using RNA sequencing. We found that MG are the most abundant myeloid cell population around RAP, whereas MAC are rarely, if ever, associated with late stages of RAP. RNA sequencing of RAP-associated MG showed that differentially expressed genes mainly contribute to immune-associated processes, including chemotaxis and migration in early RAP and proliferative capacity in late RAP, which was confirmed by immunohistochemistry. Interestingly, MG upregulated only a few angiomodulatory factors, suggesting a rather low angiogenic potential. In summary, we showed that MG are the dominant myeloid cell population at RAP sites. Moreover, MG significantly altered their transcriptional profile during RAP formation, activating immune-associated processes and exhibiting enhanced proliferation, however, without showing substantial upregulation of angiomodulatory factors.

In-Depth Molecular Profiling Specifies Human Retinal Microglia Identity

Microglia are the tissue-resident macrophages of the retina and brain, being critically involved in organ development, tissue homeostasis, and response to cellular damage. Until now, little is known about the molecular signature of human retinal microglia and how it differs from the one of brain microglia and peripheral monocytes. In addition, it is not yet clear to what extent murine retinal microglia resemble those of humans, which represents an important prerequisite for translational research. The present study applies fluorescence-activated cell sorting to isolate human retinal microglia from enucleated eyes and compares their transcriptional profile with the one of whole retinal tissue, human brain microglia as well as classical, intermediate and non-classical monocytes. Finally, human retinal microglia are compared to murine retinal microglia, isolated from Cx3cr1^GFP/+ mice. Whereas human retinal microglia exhibited a high grade of similarity in comparison to their counterparts in the brain, several enriched genes were identified in retinal microglia when compared to whole retinal tissue, as well as classical, intermediate, and non-classical monocytes. In relation to whole retina sequencing, several risk genes associated with age-related macular degeneration (AMD) and diabetic retinopathy (DR) were preferentially expressed in retinal microglia, indicating their potential pathophysiological involvement. Although a high degree of similarity was observed between human and murine retinal microglia, several species-specific genes were identified, which should be kept in mind when employing mouse models to investigate retinal microglia biology. In summary, this study provides detailed insights into the molecular profile of human retinal microglia, identifies a plethora of tissue-specific and species-specific genes in comparison to human brain microglia and murine retinal microglia, and thus highlights the significance of retinal microglia in human retinal diseases and for translational research approaches.

[Web-based gene expression analysis-paving the way to decode healthy and diseased ocular tissue]

BACKGROUND

Gene expression analysis using RNA sequencing has helped to improve the understanding of many diseases. Databases, such as the Gene Expression Omnibus database of the National Center for Biotechnology Information provide RNA sequencing raw data from various diseased tissue types but their analysis requires advanced bioinformatics skills. Therefore, specific ocular databases provide the transcriptional profiles of different ocular tissues and in addition enable intuitive web-based data analysis.

OBJECTIVE

The aim of this narrative review is to provide an overview of ocular transcriptome databases and to compare them with the Human Eye Transcriptome Atlas newly established in Freiburg.

METHODS

PubMed literature search.

RESULTS

A total of nine ocular transcriptome databases focusing on different aspects were identified. The iSyTE and Express platforms specialize in gene expression during lens and retinal development in mice, whereas retina.tigem.it, Eye in a Disk, and Spectacle focus on selected ocular tissues such as the retina. Spectacle, UCSC Cell Browser and Single Cell Portal allow intuitive exploration of single cell RNA sequencing data derived from retinal, choroid, cornea, iris, trabecular meshwork and sclera specimens. The microarray profiles of a variety of healthy ocular tissues are included in the Ocular Tissue Database. The Human Eye Transcriptome Atlas provides the largest collection of different ocular tissue types, contains the highest number of ocular diseases and is characterized by a high level of quality achieved by methodological consistency.

CONCLUSION

Ocular transcriptome databases provide comprehensive and intuitive insights into the transcriptional profiles of a variety of healthy and diseased ocular tissues. Thus, they improve our understanding of the underlying molecular mediators, support hypothesis generation and help in the search for new diagnostic and therapeutic targets for various ocular diseases.

Characterization of the Cellular Microenvironment and Novel Specific Biomarkers in Pterygia Using RNA Sequencing

With a worldwide prevalence of ~12%, pterygium is a common degenerative and environmentally triggered ocular surface disorder characterized by wing-shaped growth of conjunctival tissue onto the cornea that can lead to blindness if left untreated. This study characterizes the transcriptional profile and the cellular microenvironment of conjunctival pterygia and identifies novel pterygia-specific biomarkers. Formalin-fixed and paraffin-embedded pterygia as well as healthy conjunctival specimens were analyzed using MACE RNA sequencing (n = 8 each) and immunohistochemistry (pterygia n = 7, control n = 3). According to the bioinformatic cell type enrichment analysis using xCell, the cellular microenvironment of pterygia was characterized by an enrichment of myofibroblasts, T-lymphocytes and various antigen-presenting cells, including dendritic cells and macrophages. Differentially expressed genes that were increased in pterygia compared to control tissue were mainly involved in autophagy (including DCN, TMBIM6), cellular response to stress (including TPT1, DDX5) as well as fibroblast proliferation and epithelial to mesenchymal transition (including CTNNB1, TGFBR1, and FN1). Immunohistochemical analysis confirmed a significantly increased FN1 stromal immunoreactivity in pterygia when compared to control tissue. In addition, a variety of factors involved in apoptosis were significantly downregulated in pterygia, including LCN2, CTSD, and NISCH. Furthermore, 450 pterygia-specific biomarkers were identified by including transcriptional data of different ocular surface pathologies serving as controls (training group), which were then validated using transcriptional data of cultured human pterygium cells. Among the most pterygia-specific factors were transcripts such as AHNAK, RTN4, TPT1, FSTL1, and SPARC. Immunohistochemical validation of SPARC revealed a significantly increased stromal immunoreactivity in pterygia when compared to controls, most notably in vessels and intravascular vessel wall-adherent mononuclear cells. Taken together, the present study provides new insights into the cellular microenvironment and the transcriptional profile of pterygia, identifies new and specific biomarkers and in addition to fibrosis-related genes, uncovers autophagy, stress response and apoptosis modulation as pterygium-associated processes. These findings expand our understanding of the pathophysiology of pterygia, provide new diagnostic tools, and may enable new targeted therapeutic options for this common and sight-threatening ocular surface disease.

Comparative transcriptome analysis of human and murine choroidal neovascularization identifies fibroblast growth factor inducible-14 as phylogenetically conserved mediator of neovascular age-related macular degeneration

BACKGROUND

Visual outcome of patients with neovascular age-related macular degeneration has significantly improved during the last years following the introduction of anti-vascular endothelial growth factor (VEGF) therapy. However, about one third of patients show persistent exudation and decreasing visual acuity despite recurrent anti-VEGF treatment, which implies a role of other, still unknown proangiogenic mediators.

METHODS

The present study applied transcriptional profiling of human and mouse (C57BL/6J wildtype) choroidal neovascularization (CNV) membranes each with reference to healthy control tissue to identify yet unrecognized mediators of CNV formation. Key factors were further investigated by immunohistochemistry as well as by intravitreal inhibition experiments and multiplex protein assays in the laser-induced CNV mouse model.

FINDINGS

Transcriptional profiles of CNV membranes were characterized by enhanced activation of blood vessel development, cytoskeletal organization, and cytokine production, with angiogenesis and wound healing processes predominating in humans and activation of immune processes in mice. Besides several species-specific factors, 95 phylogenetically conserved CNV-associated genes were detected, among which fibroblast growth factor inducible-14 (FN14), a member of the tumor necrosis factor (TNF) receptor family, was identified as a key player of CNV formation. Blocking the pathway by intravitreal injection of a FN14 decoy receptor modulated the cytokine profile - most notably IL-6 - and led to a significant reduction of CNV size in vivo.

INTERPRETATION

This study characterizes the transcriptome of human and mouse CNV membranes in an unprejudiced manner and identifies FN14 as a phylogenetically conserved mediator of CNV formation and a promising new therapeutic target for neovascular AMD.

FUNDING

This study was funded by the Helmut Ecker Foundation and the Volker Homann Foundation.

Immunosenescence in Choroidal Neovascularization (CNV)-Transcriptional Profiling of Naïve and CNV-Associated Retinal Myeloid Cells during Aging

Immunosenescence is considered a possible factor in the development of age-related macular degeneration and choroidal neovascularization (CNV). However, age-related changes of myeloid cells (MCs), such as microglia and macrophages, in the healthy retina or during CNV formation are ill-defined. In this study, Cx3cr1-positive MCs were isolated by fluorescence-activated cell sorting from six-week (young) and two-year-old (old) Cx3cr1^GFP/+ mice, both during physiological aging and laser-induced CNV development. High-throughput RNA-sequencing was performed to define the age-dependent transcriptional differences in MCs during physiological aging and CNV development, complemented by immunohistochemical characterization and the quantification of MCs, as well as CNV size measurements. These analyses revealed that myeloid cells change their transcriptional profile during both aging and CNV development. In the steady state, senescent MCs demonstrated an upregulation of factors contributing to cell proliferation and chemotaxis, such as Cxcl13 and Cxcl14, as well as the downregulation of microglial signature genes. During CNV formation, aged myeloid cells revealed a significant upregulation of angiogenic factors such as Arg1 and Lrg1 concomitant with significantly enlarged CNV and an increased accumulation of MCs in aged mice in comparison to young mice. Future studies need to clarify whether this observation is an epiphenomenon or a causal relationship to determine the role of immunosenescence in CNV formation.

Application of the 3' mRNA-Seq using unique molecular identifiers in highly degraded RNA derived from formalin-fixed, paraffin-embedded tissue

Background

Archival formalin-fixed, paraffin-embedded (FFPE) tissue samples with clinical and histological data are a singularly valuable resource for developing new molecular biomarkers. However, transcriptome analysis remains challenging with standard mRNA-seq methods as FFPE derived-RNA samples are often highly modified and fragmented. The recently developed 3′ mRNA-seq method sequences the 3′ region of mRNA using unique molecular identifiers (UMI), thus generating gene expression data with minimal PCR bias. In this study, we evaluated the performance of 3′ mRNA-Seq using Lexogen QuantSeq 3′ mRNA-Seq Library Prep Kit FWD with UMI, comparing with TruSeq Stranded mRNA-Seq and RNA Exome Capture kit. The fresh-frozen (FF) and FFPE tissues yielded nucleotide sizes range from 13 to > 70% of DV200 values; input amounts ranged from 1 ng to 100 ng for validation.

Results

The total mapped reads of QuantSeq 3′ mRNA-Seq to the reference genome ranged from 99 to 74% across all samples. After PCR bias correction, 3 to 56% of total sequenced reads were retained. QuantSeq 3′ mRNA-Seq data showed highly reproducible data across replicates in Universal Human Reference RNA (UHR, R > 0.94) at input amounts from 1 ng to 100 ng, and FF and FFPE paired samples (R = 0.92) at 10 ng. Severely degraded FFPE RNA with ≤30% of DV200 value showed good concordance (R > 0.87) with 100 ng input. A moderate correlation was observed when directly comparing QuantSeq 3′ mRNA-Seq data with TruSeq Stranded mRNA-Seq (R = 0.78) and RNA Exome Capture data (R > 0.67).

Conclusion

In this study, QuantSeq 3′ mRNA-Seq with PCR bias correction using UMI is shown to be a suitable method for gene quantification in both FF and FFPE RNAs. 3′ mRNA-Seq with UMI may be applied to severely degraded RNA from FFPE tissues generating high-quality sequencing data.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08068-1.

Cardiac SARS-CoV-2 infection is associated with pro-inflammatory transcriptomic alterations within the heart

Aims

Cardiac involvement in COVID-19 is associated with adverse outcome. However, it is unclear whether cell-specific consequences are associated with cardiac SARS-CoV-2 infection. Therefore, we investigated heart tissue utilizing in situ hybridization, immunohistochemistry, and RNA-sequencing in consecutive autopsy cases to quantify virus load and characterize cardiac involvement in COVID-19.

Methods and results

In this study, 95 SARS-CoV-2-positive autopsy cases were included. A relevant SARS-CoV-2 virus load in the cardiac tissue was detected in 41/95 deceased (43%). Massive analysis of cDNA ends (MACE)-RNA-sequencing was performed to identify molecular pathomechanisms caused by the infection of the heart. A signature matrix was generated based on the single-cell dataset ‘Heart Cell Atlas’ and used for digital cytometry on the MACE-RNA-sequencing data. Thus, immune cell fractions were estimated and revealed no difference in immune cell numbers in cases with and without cardiac infection. This result was confirmed by quantitative immunohistological diagnosis. MACE-RNA-sequencing revealed 19 differentially expressed genes (DEGs) with a q-value <0.05 (e.g. up: IFI44L, IFT3, TRIM25; down: NPPB, MB, MYPN). The upregulated DEGs were linked to interferon pathways and originate predominantly from endothelial cells. In contrast, the downregulated DEGs originate predominately from cardiomyocytes. Immunofluorescent staining showed viral protein in cells positive for the endothelial marker ICAM1 but rarely in cardiomyocytes. The Gene Ontology (GO) term analysis revealed that downregulated GO terms were linked to cardiomyocyte structure, whereas upregulated GO terms were linked to anti-virus immune response.

Conclusion

This study reveals that cardiac infection induced transcriptomic alterations mainly linked to immune response and destruction of cardiomyocytes. While endothelial cells are primarily targeted by the virus, we suggest cardiomyocyte destruction by paracrine effects. Increased pro-inflammatory gene expression was detected in SARS-CoV-2-infected cardiac tissue but no increased SARS-CoV-2 associated immune cell infiltration was observed.

Viral S protein histochemistry reveals few potential SARS-CoV-2 entry sites in human ocular tissues

Despite the reported low expression of the primary SARS-CoV-2 receptor ACE2 in distinct ocular tissues, some clinical evidence suggests that SARS-CoV-2 can infect the eye. In this study, we explored potential entry sites for SARS-CoV-2 by viral S protein histochemistry on various ocular tissues and compared the staining patterns with RNA and protein expression of TMPRSS2 and ACE2. Potential viral entry sites were investigated by histochemistry using tagged recombinant viral S protein on 52 ocular tissue samples including specimens of the cornea, conjunctiva, lid margin, lacrimal gland tissue, retina, choroid, and RPE. In addition, ACE2 and TMPRSS2 immunohistochemistry were performed on the same ocular tissue, each with distinct antibodies binding to different epitopes. Lung tissue samples were used as positive controls. Finally, bulk RNA sequencing (RNA-Seq) was used to determine the expression of ACE2 and its auxiliary factors in the tissues mentioned above. S protein histochemistry revealed a positive staining in lung tissue but absent staining in the cornea, the conjunctiva, eye lid samples, the lacrimal glands, the retina and the optic nerve which was supported by hardly any immunoreactivity for ACE2 and TMPRSS2 and scarce ACE2 and TMPRSS2 RNA expression. Negligible staining with antibodies targeting ACE2 or TMPRSS2 was seen in the main and accessory lacrimal glands. In contrast, ocular staining (S protein, ACE2, TMPRSS2) was distinctly present in pigmented cells of the RPE and choroid, as well as in the ciliary body and the iris stroma. S protein histochemistry revealed hardly any SARS-CoV-2 entry sites in all ocular tissues examined. Similarly, no significant ACE2 or TMPRSS2 expression was found in extra- and intraocular tissue. While this study suggest a rather low risk of ocular infection with SARS-CoV-2, it should be noted, that potential viral entry sites may increase in response to inflammation or in certain disease states.

Time- and Stimulus-Dependent Characteristics of Innate Immune Cells in Organ-Cultured Human Corneal Tissue

PURPOSE

The pattern of immune cells infiltrating the corneal stroma has been extensively studied in mice, but data on human tissue have been far less elaborate. To further characterize the number and differentiation state of resident immune cells in organ-cultured human corneal tissue, we employed a comprehensive bioinformatic deconvolution (xCell) of bulk RNA-sequencing (RNA-seq) data, immunohistochemistry (IHC), and flow cytometry (FC).

METHODS

A transcriptome-based analysis of immune cell types in human corneal samples was performed. The results were validated by IHC, focusing on the identification of pro-inflammatory (M1) and regulatory (M2) macrophages. A protocol was established to identify these 2 different macrophage populations in human corneal tissue by means of FC. Subsequently, corneal samples in organ culture were differentially stimulated by IL-10, IL-4 &amp; IL-13, or LPS and macrophage populations were evaluated regarding their response to these stimuli. Furthermore, cell survival was analyzed in correlation with time in organ culture.

RESULTS

xCell-based mathematical deconvolution of bulk RNA-seq data revealed the presence of CD8 T cells, Th17 cells, dendritic cells, and macrophages as the predominant immune cell types in organ-cultured human corneal tissue. Furthermore, RNA-seq allowed the detection of different macrophage marker genes in corneal samples, including PTPRC (CD45), ITGAM (CD11b), CD14, and CD74. Our RNA-seq data showed no evidence of a relevant presence of monocytes in human corneal tissue. The presence of different macrophage subtypes was confirmed by IHC. The disintegration and subsequent FC analysis of human corneal samples showed the presence of both M1 (HLA-DR+, CD282+, CD86+, and CD284+) and M2 (CD163+ and CD206+) macrophage subtypes. Furthermore, we found that the total number of macrophages in corneal samples decreased more than the total cell count with increasing tissue culture time. Treatment with IL-10 led to higher total cell counts per cornea and to an increased expression of the M2 marker CD163 (p < 0.05) while expression levels of various M1 macrophage markers were not significantly reduced by interleukin treatment.

CONCLUSIONS

Regarding different macrophage populations, untreated human corneas showed more M1 than M2 macrophages. With increasing organ culture time, these macrophages decreased. In terms of cell dynamics, adding interleukins to the organ culture medium influenced the phenotype of macrophages within the cornea as detected by FC. Modifying the immunomodulatory properties of human grafts appears a promising approach to further reduce the risk of graft rejection in patients. In this context, treatment with interleukins was more effective in upregulating M2 macrophages than in suppressing M1 macrophages in corneal tissue.

Mapping the origin and fate of myeloid cells in distinct compartments of the eye by single-cell profiling

Similar to the brain, the eye is considered an immune‐privileged organ where tissue‐resident macrophages provide the major immune cell constituents. However, little is known about spatially restricted macrophage subsets within different eye compartments with regard to their origin, function, and fate during health and disease. Here, we combined single‐cell analysis, fate mapping, parabiosis, and computational modeling to comprehensively examine myeloid subsets in distinct parts of the eye during homeostasis. This approach allowed us to identify myeloid subsets displaying diverse transcriptional states. During choroidal neovascularization, a typical hallmark of neovascular age‐related macular degeneration (AMD), we recognized disease‐specific macrophage subpopulations with distinct molecular signatures. Our results highlight the heterogeneity of myeloid subsets and their dynamics in the eye that provide new insights into the innate immune system in this organ which may offer new therapeutic targets for ophthalmological diseases.

Secreted Phosphoprotein 1 Expression in Retinal Mononuclear Phagocytes Links Murine to Human Choroidal Neovascularization

Age-related macular degeneration (AMD) represents the most common cause of blindness in the elderly in the Western world. An impairment of the outer blood-retina barrier and a localized inflammatory microenvironment cause sprouting of choroidal neovascular membranes (CNV) in neovascular AMD that are in intimate contact with surrounding myeloid cells, such as retinal microglia, and ultimately lead to visual impairment. The discovery of novel target molecules to interfere with angiogenesis and inflammation is vital for future treatment approaches in AMD patients. To explore the transcriptional profile and the function of retinal microglia at sites of CNV, we performed a comprehensive RNA-seq analysis of retinal microglia in the mouse model of laser-induced choroidal neovascularization (mCNV). Here, we identified the angiogenic factor Osteopontin (Opn), also known as "secreted phosphoprotein 1" (Spp1), as one of the most highly expressed genes in retinal microglia in the course of CNV formation. We confirmed the presence of SPP1 at the lesion site in recruited retinal microglia in Cx3cr1 ^CreER:Rosa26-tdTomato reporter mice by confocal microscopy and in whole retinal tissue lysates by ELISA highlighting a massive local production of SPP1. Inhibition of SPP1 by intravitreal injection of an anti-SPP1 antibody significantly increased the lesion size compared to IgG-treated control eyes. In line with our results in rodents, we found an increased SPP1 mRNA expression in surgically extracted human choroidal neovascular (hCNV) membranes by the quantitative RNA-seq approach of massive analysis of cDNA ends (MACE). Numerous IBA1⁺SPP1⁺ myeloid cells were detected in human CNV membranes. Taken together, these results highlight the importance of SPP1 in the formation of CNV and potentially offer new opportunities for therapeutic intervention by modulating the SPP1 pathway.

MACE RNA sequencing analysis of conjunctival squamous cell carcinoma and papilloma using formalin-fixed paraffin-embedded tumor tissue

Recent advances in the field of biomedical research allow for elucidation of the transcriptional signature of rare tumors such as conjunctival squamous cell carcinoma (SCC). In this study we compare its expression profile to conjunctival papilloma (Pap) and healthy conjunctival tissue (Ctrl) and develop a classification tool to differentiate these entities. Seven conjunctival SCC, seven Pap and ten Ctrl were formalin-fixed and paraffin-embedded (FFPE) and analyzed using Massive Analysis of cDNA Ends (MACE) RNA sequencing. Differentially expressed genes (DEG) and gene ontology (GO) clusters were explored and the abundance of involved cell types was quantified by xCell. Finally, a classification model was developed to distinguish SCC from Pap and Ctrl. Among the most prominent DEG in SCC a plethora of keratins were upregulated when compared to Pap and Ctrl. xCell analysis revealed an enrichment of immune cells, including activated dendritic cells and T-helper type 1 cells (Th1), in SCC when compared to Ctrl. The generated classification model could reliably discriminate between the three entities according to the expression pattern of 30 factors. This study provides a transcriptome-wide gene expression profile of rare conjunctival SCC. The analysis identifies distinct keratins, as well as dendritic and Th1 cells as important mediators in SCC. Finally, the provided gene expression classifier may become an aid to the conventional histological classification of conjunctival tumors in uncertain cases.

Transcriptional characterization of conjunctival melanoma identifies the cellular tumor microenvironment and prognostic gene signatures

This study characterizes the transcriptome and the cellular tumor microenvironment (TME) of conjunctival melanoma (CM) and identifies prognostically relevant biomarkers. 12 formalin-fixed and paraffin-embedded CM were analyzed by MACE RNA sequencing, including six cases each with good or poor clinical outcome, the latter being defined by local recurrence and/or systemic metastases. Eight healthy conjunctival specimens served as controls. The TME of CM, as determined by bioinformatic cell type enrichment analysis, was characterized by the enrichment of melanocytes, pericytes and especially various immune cell types, such as plasmacytoid dendritic cells, natural killer T cells, B cells and mast cells. Differentially expressed genes between CM and control were mainly involved in inhibition of apoptosis, proteolysis and response to growth factors. POU3F3, BIRC5 and 7 were among the top expressed genes associated with inhibition of apoptosis. 20 genes, among them CENPK, INHA, USP33, CASP3, SNORA73B, AAR2, SNRNP48 and GPN1, were identified as prognostically relevant factors reaching high classification accuracy (area under the curve: 1.0). The present study provides new insights into the TME and the transcriptional profile of CM and additionally identifies new prognostic biomarkers. These results add new diagnostic tools and may lead to new options of targeted therapy for CM.