Elucidation of the Epitranscriptomic RNA Modification Landscape of Chikungunya Virus.

The genomes of positive-sense (+) single-stranded RNA (ssRNA) viruses are believed to be subjected to a wide range of RNA modifications. In this study, we focused on the chikungunya virus (CHIKV) as a model (+) ssRNA virus to study the landscape of viral RNA modification in infected human cells. Among the 32 distinct RNA modifications analysed by mass spectrometry, inosine was found enriched in the genomic CHIKV RNA. However, orthogonal validation by Illumina RNA-seq analyses did not identify any inosine modification along the CHIKV RNA genome. Moreover, CHIKV infection did not alter the expression of ADAR1 isoforms, the enzymes that catalyse the adenosine to inosine conversion. Together, this study highlights the importance of a multidisciplinary approach to assess the presence of RNA modifications in viral RNA genomes.

Endoplasmic Reticulum Stress and Unfolded Protein Response in Vernal Keratoconjunctivitis.

Purpose: Vernal keratoconjunctivitis (VKC) is an ocular allergic disease characterized by a type 2 inflammation, tissue remodeling, and low quality of life for the affected patients. We investigated the involvement of endoplasmic reticulum (ER) stress and unfolded protein response in VKC. Methods: Conjunctival imprints from VKC patients and normal subjects (CTs) were collected, and RNA was isolated, reverse transcribed, and analyzed with the Affymetrix microarray. Differentially expressed genes between VKC patients and CTs were evaluated. Genes related to ER stress, apoptosis, and autophagy were further considered. VKC and CT conjunctival biopsies were analyzed by immunohistochemistry (IHC) with specific antibodies against unfolded protein response (UPR), apoptosis, and inflammation. Conjunctival fibroblast and epithelial cell cultures were exposed to the conditioned medium of activated U937 monocytes and analyzed by quantitative PCR for the expression of UPR, apoptosis, autophagy, and inflammatory markers. Results: ER chaperones HSPA5 (GRP78/BiP) and HYOU1 (GRP170) were upregulated in VKC patients compared to CTs. Genes encoding for ER transmembrane proteins, PKR-like ER kinase (PERK), activating transcription factor 6 (ATF6), ER-associated degradation (ERAD), and autophagy were upregulated, but not those related to apoptosis. Increased positive reactivity of BiP and ATF6 and unchanged expression of apoptosis markers were confirmed by IHC. Cell cultures in stress conditions showed an overexpression of UPR, proinflammatory, apoptosis, and autophagy markers. Conclusions: A significant overexpression of genes encoding for ER stress, UPR, and pro-inflammatory pathway components was reported for VKC. Even though these pathways may lead to ER homeostasis, apoptosis, or inflammation, ER stress in VKC may predominantly contribute to promote inflammation.

Evaluation of Autophagy in Conjunctival Fibroblasts.

Vernal keratoconjunctivitis (VKC) is a serious eye allergy characterized by poorly understood pathogenic mechanisms and a lack of effective treatments. Autophagy, a process involved in both triggering and suppressing immune and inflammatory responses, plays a role in VKC's pathophysiology. Understanding autophagy's involvement in VKC could lead to new treatment possibilities, such as utilizing specific topical substances to induce or inhibit autophagy and prevent severe complications of this eye condition. In our current protocol, we present a robust methodology established in our laboratory for studying autophagy in primary conjunctival fibroblasts. We assess autophagy through techniques like immunocytochemistry, immunoblotting, and qPCR.

Upregulation of Notch Signaling and Cell-Differentiation Inhibitory Transcription Factors in Stable Chronic Obstructive Pulmonary Disease Patients.

Notch signaling is involved in the prevention of cell differentiation and cell fate in various organs, including the lungs. We aimed to determine the transcriptomic and protein expression of Notch receptors, their ligands, and related transcription factors in stable COPD. The expression and localization of Notch receptors, their ligands, and related transcription factors were measured in bronchial biopsies of individuals with stable mild/moderate (MCOPD) (n = 18) or severe/very severe (SCOPD) (n = 16) COPD, control smokers (CSs) (n = 13), and control nonsmokers (CNSs) (n = 11), and in the lung parenchyma of those with MCOPD (n = 13), CSs (n = 10), and CNSs (n = 10) using immunohistochemistry, ELISA tests, and transcriptome analyses. In the bronchial biopsies, Notch4 and HES7 significantly increased in the lamina propria of those with SCOPD compared to those with MCOPD, CSs, and CNSs. In the peripheral lung bronchiolar epithelium, Notch1 significantly increased in those with MCOPD and CSs compared to CNSs. ELISA tests of lung parenchyma homogenates showed significantly increased Notch2 in those with MCOPD compared to CSs and CNSs. Transcriptomic data in lung parenchyma showed increased DLL4 and HES1 mRNA levels in those with MCOPD and CSs compared to CNSs. These data show the increased expression of the Notch pathway in the lungs of those with stable COPD. These alterations may play a role in impairing the regenerative-reparative responses of diseased bronchioles and lung parenchyma.

Comparative analysis of Presence-Absence gene Variations in five hard tick species: impact and functional considerations.

Tick species are vectors of harmful human and animal diseases, and their expansion is raising concerns under the global environmental changes' scenario. Ticks host and transmit bacteria, protozoa and viruses, making the understanding of host-pathogen molecular pathways critical to development of effective disease control strategies. Despite the considerable sizes and repeat contents of tick genomes, individual tick genomics is perhaps the most effective approach to reveal genotypic traits of interest. Presence-Absence gene Variations (PAVs) can contribute to individual differences within species, with dispensable genes carried by subsets of individuals possibly underpinning functional significance at individual or population-levels. We exploited 350 resequencing datasets of Dermacentor silvarum, Haemaphysalis longicornis, Ixodes persulcatus, Rhipicephalus microplus and Rhipicephalus sanguineus hard tick specimens to reveal the extension of PAV and the conservation of dispensable genes among individuals and, comparatively, between species. Overall, we traced 550-3,346 dispensable genes per species and were able to reconstruct 5.3-7 Mb of genomic regions not included in the respective reference genomes, as part of the tick pangenomes. Both dispensable genes and de novo predicted genes indicated that PAVs preferentially impacted mobile genetic elements in these tick species.

FicD genes in invertebrates: A tale of transposons, pathogenic and integrated viruses.

Many gene families are shared across the tree of life between distantly related species because of horizontal gene transfers (HGTs). However, the frequency of HGTs varies strongly between gene families and biotic realms suggesting differential selection pressures and functional bias. One gene family with a wide distribution are FIC-domain containing enzymes (FicDs). FicDs catalyze AMPylation, a post-translational protein modification consisting in the addition of adenosine monophosphate to accessible residues of target proteins. Beside the well-known conservation of FicDs in deuterostomes, we report the presence of a conserved FicD gene ortholog in a large number of protostomes and microbial eukaryotes. We also reported additional FicD gene copies in the genomes of some rotifers, parasitic worms and bivalves. A few dsDNA viruses of these invertebrates, including White spot syndrome virus, Cherax quadricarinatus iridovirus, Ostreid herpesvirus-1 and the beetle nudivirus, carry copies of FicDs, with phylogenetic analysis suggesting a common origin of these FicD copies and the duplicated FicDs of their invertebrate hosts. HGTs and gene duplications possibly mediated by endogenous viruses or genetic mobile elements seem to have contributed to the transfer of AMPylation ability from bacteria and eukaryotes to pathogenic viruses, where this pathway could have been hijacked to promote viral infection.

Viral Diversity in Benthic Abyssal Ecosystems: Ecological and Methodological Considerations.

Viruses are the most abundant 'biological entities' in the world's oceans. However, technical and methodological constraints limit our understanding of their diversity, particularly in benthic abyssal ecosystems (>4000 m depth). To verify advantages and limitations of analyzing virome DNA subjected either to random amplification or unamplified, we applied shotgun sequencing-by-synthesis to two sample pairs obtained from benthic abyssal sites located in the North-eastern Atlantic Ocean at ca. 4700 m depth. One amplified DNA sample was also subjected to single-molecule long-read sequencing for comparative purposes. Overall, we identified 24,828 viral Operational Taxonomic Units (vOTUs), belonging to 22 viral families. Viral reads were more abundant in the amplified DNA samples (38.5-49.9%) compared to the unamplified ones (4.4-5.8%), with the latter showing a greater viral diversity and 11-16% of dsDNA viruses almost undetectable in the amplified samples. From a procedural point of view, the viromes obtained by direct sequencing (without amplification step) provided a broader overview of both ss and dsDNA viral diversity. Nevertheless, our results suggest that the contextual use of random amplification of the same sample and long-read technology can improve the assessment of viral assemblages by reducing off-target reads.

Impaired autophagy in the lower airways and lung parenchyma in stable COPD.

Background: There is increasing evidence of autophagy activation in COPD, but its role is complex and probably regulated through cell type-specific mechanisms. This study aims to investigate the autophagic process at multiple levels within the respiratory system, using different methods to clarify conflicting results reported so far. Methods: This cross-sectional study was performed on bronchial biopsies and peripheral lung samples obtained from COPD patients (30 and 12 per sample type, respectively) and healthy controls (25 and 22 per sample type, respectively), divided by smoking history. Subjects were matched for age and smoking history. We analysed some of the most important proteins involved in autophagosome formation, such as LC3 and p62, as well as some molecules essential for lysosome function, such as lysosome-associated membrane protein 1 (LAMP1). Immunohistochemistry was used to assess the autophagic process in both sample types. ELISA and transcriptomic analysis were performed on lung samples. Results: We found increased autophagic stimulus in smoking subjects, regardless of respiratory function. This was revealed by immunohistochemistry through a significant increase in LC3 (p<0.01) and LAMP1 (p<0.01) in small airway bronchiolar epithelium, alveolar septa and alveolar macrophages. Similar results were obtained in bronchial biopsy epithelium by evaluating LC3B (p<0.05), also increased in homogenate lung tissue using ELISA (p<0.05). Patients with COPD, unlike the others, showed an increase in p62 by ELISA (p<0.05). No differences were found in transcriptomics analysis. Conclusions: Different techniques, applied at post-transcriptional level, confirm that cigarette smoke stimulates autophagy at multiple levels inside the respiratory system, and that autophagy failure may characterise COPD.

Bone Morphogenic Proteins and Their Antagonists in the Lower Airways of Stable COPD Patients.

BACKGROUND: Bone morphogenic proteins (BMPs) and their antagonists are involved in the tissue development and homeostasis of various organs. OBJECTIVE: To determine transcriptomic and protein expression of BMPs and their antagonists in stable COPD. METHODS: We measured the expression and localization of BMPs and some relevant antagonists in bronchial biopsies of stable mild/moderate COPD (MCOPD) (n = 18), severe/very severe COPD (SCOPD) (n = 16), control smokers (CS) (n = 13), and control non-smokers (CNS) (n = 11), and in lung parenchyma of MCOPD (n = 9), CS (n = 11), and CNS (n = 9) using immunohistochemistry and transcriptome analysis, in vitro after the stimulation of the 16HBE cells. RESULTS: In bronchial biopsies, BMP4 antagonists CRIM1 and chordin were increased in the bronchial epithelium and lamina propria of COPD patients. BMP4 expression was decreased in the bronchial epithelium of SCOPD and MCOPD compared to CNS. Lung transcriptomic data showed non-significant changes between groups. CRIM1 and chordin were significantly decreased in the alveolar macrophages and alveolar septa in COPD patients. External 16HBE treatment with BMP4 protein reduced the bronchial epithelial cell proliferation. CONCLUSIONS: These data show an imbalance between BMP proteins and their antagonists in the lungs of stable COPD. This imbalance may play a role in the remodeling of the airways, altering the regenerative-reparative responses of the diseased bronchioles and lung parenchyma.

Meta-Analysis of Keratoconus Transcriptomic Data Revealed Altered RNA Editing Levels Impacting Keratin Genomic Clusters.

Introduction: Keratoconus (KC) is an ocular disorder with a multifactorial origin. Transcriptomic analyses (RNA-seq) revealed deregulations of coding (mRNA) and non-coding RNAs (ncRNAs) in KC, suggesting that mRNA-ncRNA co-regulations can promote the onset of KC. The present study investigates the modulation of RNA editing mediated by the adenosine deaminase acting on dsRNA (ADAR) enzyme in KC. Materials: The level of ADAR-mediated RNA editing in KC and healthy corneas were determined by two indexes in two different sequencing datasets. REDIportal was used to localize known editing sites, whereas new putative sites were de novo identified in the most extended dataset only and their possible impact was evaluated. Western Blot analysis was used to measure the level of ADAR1 in the cornea from independent samples. Results: KC was characterized by a statistically significant lower RNA-editing level compared to controls, resulting in a lower editing frequency, and less edited bases. The distribution of the editing sites along the human genome showed considerable differences between groups, particularly relevant in the chromosome 12 regions encoding for Keratin type II cluster. A total of 32 recoding sites were characterized, 17 representing novel sites. JUP, KRT17, KRT76, and KRT79 were edited with higher frequencies in KC than in controls, whereas BLCAP, COG3, KRT1, KRT75, and RRNAD1 were less edited. Both gene expression and protein levels of ADAR1 appeared not regulated between diseased and controls. Conclusions: Our findings demonstrated an altered RNA-editing in KC possibly linked to the peculiar cellular conditions. The functional implications should be further investigated.

Tracing RNA viruses associated with Nudibranchia gastropods.

Background: Nudibranchia is an under-studied taxonomic group of gastropods, including more than 3,000 species with colourful and extravagant body shapes and peculiar predatory and defensive strategies. Although symbiosis with bacteria has been reported, no data are available for the nudibranch microbiome nor regarding viruses possibly associated with these geographically widespread species. Methods: Based on 47 available RNA sequencing datasets including more than two billion reads of 35 nudibranch species, a meta-transcriptome assembly was constructed. Taxonomic searches with DIAMOND, RNA-dependent-RNA-polymerase identification with palmscan and viral hallmark genes identification by VirSorter2 in combination with CheckV were applied to identify genuine viral genomes, which were then annotated using CAT. Results: A total of 20 viral genomes were identified as bona fide viruses, among 552 putative viral contigs resembling both RNA viruses of the Negarnaviricota, Pisuviricota, Kitrinoviricota phyla and actively transcribing DNA viruses of the Cossaviricota and Nucleocytoviricota phyla. The 20 commonly identified viruses showed similarity with RNA viruses identified in other RNA-seq experiments and can be putatively associated with bacteria, plant and arthropod hosts by co-occurence analysis. The RNA samples having the highest viral abundances showed a heterogenous and mostly sample-specific distribution of the identified viruses, suggesting that nudibranchs possess diversified and mostly unknown viral communities.

ADAR-Editing during Ostreid Herpesvirus 1 Infection in Crassostrea gigas: Facts and Limitations.

Ostreid herpesvirus-1 (OsHV-1) RNAs are enzymatically modified by A-to-I conversions during the infection of Crassostrea gigas. The increase of ADAR1 expression and hyper-editing activity parallel to OsHV-1 RNAs suggests a functional connection between dsRNA editing and antiviral responses. We analyzed 87 RNA-seq data sets from immuno-primed, resistant, and susceptible oysters exposed to OsHV-1 to compare the ADAR hyper-editing levels on host and viral transcripts and trace hyper-editing on the oyster genes. Host RNAs were more hyper-edited than viral RNAs, despite the increased editing of viral RNAs in late infection phases. A set of genes, representing ∼0.5% of the oyster transcriptome and including several tripartite motif-containing sequences, were constantly hyper-edited. Conversely, we identified genes involved in antiviral response, miRNA maturation, and epigenetic regulation that were hyper-edited in specific conditions only. Despite technical and biological bottlenecks that hamper the understanding of the bivalve "RNA editome," available tools and technologies can be adapted to bivalve mollusks. IMPORTANCE Ostreid herpesvirus-1 (OsHV-1) is a harmful pathogen of bivalve species, such as oysters. However, knowledge is lacking about host-virus interactions at the molecular level, hampering the possibility of a correct management of viral outbreaks and related massive mortalities. Notably, OsHV-1 transcripts are massively modified by host RNA editing enzyme during infection, resulting in multiple A-to-I variations along RNAs assuming double-strand conformations. The impact of these modifications on host transcripts is, however, not completely clear. Analyzing RNA-seq data of oysters infected with OsHV-1, we revealed that ∼0.5% of the oyster transcriptome is always enzymatically modified by ADAR, whereas genes involved in antiviral response, miRNA maturation, and epigenetic regulation were hyper-edited in specific conditions only. Despite our results, relevant technical bottlenecks impair an accurate quantification of RNA editing events, making necessary an approach specifically dedicated to the progressive understanding of oyster "RNA editome."

Paired miRNA and RNA sequencing provides a first insight into molecular defense mechanisms of Scapharca broughtonii during ostreid herpesvirus-1 infection.

Ostreid herpesvirus 1 (OsHV-1) infection caused mortalities with relevant economic losses in bivalve aquaculture industry worldwide. Initially described as an oyster pathogen, OsHV-1 can infect other bivalve species, like the blood clam Scapharca broughtonii. However, at present, little is known about the molecular interactions during OsHV-1 infection in the blood clam. We produced paired miRNA and total RNA-seq data to investigate the blood clam transcriptional changes from 0 to 72 h after experimental infection with OsHV-1. High-throughput miRNA sequencing of 24 libraries revealed 580 conserved and 270 new blood clam miRNAs, whereas no genuine miRNA was identified for OsHV-1. Total 88-203 differently expressed miRNAs were identified per time point, mostly up-regulated and mainly targeting metabolic pathways. Most of the blood clam mRNAs, in contrast, were down-regulated up to 60 h post-injection, with the trend analysis revealing the activation of immune genes only when comparing the early and latest stage of infection. Taken together, paired short and long RNA data suggested a miRNA-mediated down-regulation of host metabolic and energetic processes as a possible antiviral strategy during early infection stages, whereas antiviral pathways appeared upregulated only at late infection.

Tracing the SARS-CoV-2 infection on the ocular surface: Overview and preliminary corneoscleral transcriptome sequencing.

COVID-19's impact on the ocular surface has already been recognized, however the molecular mechanisms induced by the infection on the ocular surface are still unclear. The aim of this paper is to provide a first overview of the transcriptional perturbations caused by SARS-CoV-2 on the ocular surface by analyzing gene expression profile of corneoscleral ring samples from post-mortem SARS-CoV-2 positive donors (PD). The presence of SARS-CoV-2 on the ocular surface, in tears and corneal tissues has rarely been detected in infected individuals in both the presence and the absence of ocular manifestations. In this preliminary study, 6 human corneoscleral tissues of 3 PD and two tissues from a negative donor (CTRL) were obtained at the local eye bank. The presence of genomic and sub-genomic SARS-CoV-2 RNAs was assessed by qRT-PCR, while transcriptome analysis (RNA-sequencing) was performed by Illumina. Principal Component Analysis (PCA), search for differentially expressed genes (DEGs) and Gene Ontology (GO)-enrichment analysis were performed. Three samples from PD were found positive for SARS-CoV-2 genomic RNA, although the absence of sub-genomic RNAs indicated an inactive virus. PCA analysis grouped 3 different clusters, one including CTRL, and the other two including, respectively, PD with undetected SARS-CoV-2 (PD-SARS-neg) and PD with detected SARS-CoV-2 (PD-SARS-pos). The DEGs in common with the 2 PD clusters included several genes associable to the interferon pathway, such as ADAMTS4, RSAD2, MMP1, IL6, ISG15 and proinflammatory cytokines. Among the down-regulated genes we found AQP5. GO analysis revealed 77 GO terms over-represented in PD-SARS-neg vs. CTRL, and 17 GO terms in PD-SARS-pos vs. CTRL. The presence of SARS-CoV-2 RNA and RNA-sequencing reads in ocular surface tissues supports the possibility that the eye acts as an entry route. The modulation of early responsive genes, together with several ISGs suggests a potential protective responsiveness of the ocular tissues to SARS-CoV-2.

A glimpse on metazoan ZNFX1 helicases, ancient players of antiviral innate immunity.

The human zinc finger NFX1-type containing 1 (ZNFX1) is an interferon-stimulated protein associated to the outer mitochondrial membrane, able to bind dsRNAs and interact with MAVS proteins, promoting type I IFN response in the early stage of viral infection. An N-terminal Armadillo (ARM)-type fold and a large helicase core (P-loop) and zinc fingers confer RNA-binding and ATPase activities to ZNFX1. We studied the phylogenetic distribution of metazoan ZNFX1s, ZNFX1 gene expression trends and genomic and protein signatures during viral infection of invertebrates. Based on 221 ZNFX1 sequences, we obtained a polyphyletic tree with a taxonomy-consistent branching at the phylum-level only. In metazoan genomes, ZNFX1 genes were found either in single copy, with up to some tens of exons in vertebrates, or in multiple copies, with one or a few exons and one of them sometimes encompassing most of the coding sequence, in invertebrates like sponges, sea urchins and mollusks. Structural analyses of selected ZNFX1 proteins showed high conservation of the helicase region (P-loop), an overall conserved region and domain architecture, an ARM-fold mostly traceable, and the presence of intrinsically disordered regions of varying length and position. The remarkable over-expression of ZNFX1 in bivalve and gastropod mollusks infected with dsDNA viruses underscores the antiviral role of ZNFX1, whereas nothing similar was found in virus-infected nematodes and corals. Whether the functional diversification reported in the C. elegans ZNFX1 occurs in other metazoan proteins remains to be established.

Viral Decoys: The Only Two Herpesviruses Infecting Invertebrates Evolved Different Transcriptional Strategies to Deflect Post-Transcriptional Editing.

The highly versatile group of Herpesviruses cause disease in a wide range of hosts. In invertebrates, only two herpesviruses are known: the malacoherpesviruses HaHV-1 and OsHV-1 infecting gastropods and bivalves, respectively. To understand viral transcript architecture and diversity we first reconstructed full-length viral genomes of HaHV-1 infecting Haliotis diversicolor supertexta and OsHV-1 infecting Scapharca broughtonii by DNA-seq. We then used RNA-seq over the time-course of experimental infections to establish viral transcriptional dynamics, followed by PacBio long-read sequencing of full-length transcripts to untangle viral transcript architectures at two selected time points. Despite similarities in genome structure, in the number of genes and in the diverse transcriptomic architectures, we measured a ten-fold higher transcript variability in HaHV-1, with more extended antisense gene transcription. Transcriptional dynamics also appeared different, both in timing and expression trends. Both viruses were heavily affected by post-transcriptional modifications performed by ADAR1 affecting sense-antisense gene pairs forming dsRNAs. However, OsHV-1 concentrated these modifications in a few genomic hotspots, whereas HaHV-1 diluted ADAR1 impact by elongated and polycistronic transcripts distributed over its whole genome. These transcriptional strategies might thus provide alternative potential roles for sense-antisense transcription in viral transcriptomes to evade the host's immune response in different virus-host combinations.

The Evolution and Diversity of Interleukin-17 Highlight an Expansion in Marine Invertebrates and Its Conserved Role in Mucosal Immunity.

The interleukin-17 (IL-17) family consists of proinflammatory cytokines conserved during evolution. A comparative genomics approach was applied to examine IL-17 throughout evolution from poriferans to higher vertebrates. Cnidaria was highlighted as the most ancient diverged phylum, and several evolutionary patterns were revealed. Large expansions of the IL-17 repertoire were observed in marine molluscs and echinoderm species. We further studied this expansion in filter-fed Mytilus galloprovincialis, which is a bivalve with a highly effective innate immune system supported by a variable pangenome. We recovered 379 unique IL-17 sequences and 96 receptors from individual genomes that were classified into 23 and 6 isoforms after phylogenetic analyses. Mussel IL-17 isoforms were conserved among individuals and shared between closely related Mytilidae species. Certain isoforms were specifically implicated in the response to a waterborne infection with Vibrio splendidus in mussel gills. The involvement of IL-17 in mucosal immune responses could be conserved in higher vertebrates from these ancestral lineages.