Insights into ADAR gene complement, expression patterns, and RNA editing landscape in Chlamys farreri.

Adenosine Deaminases Acting on RNA (ADARs) are evolutionarily conserved enzymes known to convert adenosine to inosine in double-stranded RNAs and participate in host-virus interactions. Conducting a meta-analysis of available transcriptome data, we identified and characterised eight ADAR transcripts in Chlamys farreri, a farmed marine scallop susceptible to Acute viral necrosis virus (AVNV) infections and mortality outbreaks. Accordingly, we identified six ADAR genes in the Zhikong scallop genome, revised previous gene annotations, and traced alternative splicing variants. In detail, each ADAR gene encodes a unique combination of functional domains, always including the Adenosine deaminase domain, RNA binding domains and, in one case, two copies of a Z-DNA binding domain. After phylogenetic analysis, five C. farreri ADARs clustered in the ADAR1 clade along with sequences from diverse animal phyla. Gene expression analysis indicated CF051320 as the most expressed ADAR, especially in the eye and male gonad. The other four ADAR1 genes and one ADAR2 gene exhibited variable expression levels, with CF105370 and CF051320 significantly increasing during early scallop development. ADAR-mediated single-base editing, evaluated across adult C. farreri tissues and developmental stages, was mainly detectable in intergenic regions (83 % and 85 %, respectively). Overall, the expression patterns of the six ADAR genes together with the editing and hyper-editing values computed on scallops RNA-seq samples support the adaptive value of ADAR1-mediated editing, particularly in the pre-settling larval stages.

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.

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.

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.

Parallel analysis of miRNAs and mRNAs suggests distinct regulatory networks in Crassostrea gigas infected by Ostreid herpesvirus 1.

BACKGROUND: Since 2008, the aquaculture production of Crassostrea gigas was heavily affected by mass mortalities associated to Ostreid herpesvirus 1 (OsHV-1) microvariants worldwide. Transcriptomic studies revealed the major antiviral pathways of the oyster immune response while other findings suggested that also small non-coding RNAs (sncRNA) such as microRNAs might act as key regulators of the oyster response against OsHV-1. To explore the explicit connection between small non-coding and protein-coding transcripts, we performed paired whole transcriptome analysis of sncRNA and messenger RNA (mRNA) in six oysters selected for different intensities of OsHV-1 infection. RESULTS: The mRNA profiles of the naturally infected oysters were mostly governed by the transcriptional activity of OsHV-1, with several differentially expressed genes mapping to the interferon, toll, apoptosis, and pro-PO pathways. In contrast, miRNA profiles suggested more complex regulatory mechanisms, with 15 differentially expressed miRNAs (DE-miRNA) pointing to a possible modulation of the host response during OsHV-1 infection. We predicted 68 interactions between DE-miRNAs and oyster 3'-UTRs, but only few of them involved antiviral genes. The sncRNA reads assigned to OsHV-1 rather resembled mRNA degradation products, suggesting the absence of genuine viral miRNAs. CONCLUSIONS: We provided data describing the miRNAome during OsHV-1 infection in C. gigas. This information can be used to understand the role of miRNAs in healthy and diseased oysters, to identify new targets for functional studies and, eventually to disentangle cause and effect relationships during viral infections in marine mollusks.

A-to-I editing of Malacoherpesviridae RNAs supports the antiviral role of ADAR1 in mollusks.

BACKGROUND: Adenosine deaminase enzymes of the ADAR family are conserved in metazoans. They convert adenine into inosine in dsRNAs and thus alter both structural properties and the coding potential of their substrates. Acting on exogenous dsRNAs, ADAR1 exerts a pro- or anti-viral role in vertebrates and Drosophila. RESULTS: We traced 4 ADAR homologs in 14 lophotrochozoan genomes and we classified them into ADAD, ADAR1 or ADAR2, based on phylogenetic and structural analyses of the enzymatic domain. Using RNA-seq and quantitative real time PCR we demonstrated the upregulation of one ADAR1 homolog in the bivalve Crassostrea gigas and in the gastropod Haliotis diversicolor supertexta during Ostreid herpesvirus-1 or Haliotid herpesvirus-1 infection. Accordingly, we demonstrated an extensive ADAR-mediated editing of viral RNAs. Single nucleotide variation (SNV) profiles obtained by pairing RNA- and DNA-seq data from the viral infected individuals resulted to be mostly compatible with ADAR-mediated A-to-I editing (up to 97%). SNVs occurred at low frequency in genomic hotspots, denoted by the overlapping of viral genes encoded on opposite DNA strands. The SNV sites and their upstream neighbor nucleotide indicated the targeting of selected adenosines. The analysis of viral sequences suggested that, under the pressure of the ADAR editing, the two Malacoherpesviridae genomes have evolved to reduce the number of deamination targets. CONCLUSIONS: We report, for the first time, evidence of an extensive editing of Malacoherpesviridae RNAs attributable to host ADAR1 enzymes. The analysis of base neighbor preferences, structural features and expression profiles of molluscan ADAR1 supports the conservation of the enzyme function among metazoans and further suggested that ADAR1 exerts an antiviral role in mollusks.

Expansion and loss events characterized the occurrence of MIF-like genes in bivalves.

Macrophage migration inhibitory factor (MIF) dynamically connects innate and adaptive immune systems in vertebrate animals, allowing highly orchestrated systemic responses to various insults. The occurrence of MIF-like genes in non-vertebrate organisms suggests its origin from an ancestral metazoan gene, whose function is still a matter of debate. In the present work, by analyzing available genomic and transcriptomic data from bivalve mollusks, we identified 137 MIF-like sequences, which were classified into three types, based on phylogeny and conservation of key residues: MIF, D-DT, and the lineage-specific type MDL. Comparative genomics revealed syntenic conservation of homologous genes at the family level, the loss of D-DT in the Ostreidae family as well as the expansion of MIF-like genes in the Mytilidae family, possibly underpinning the neofunctionalization of duplicated gene copies. In M. galloprovincialis, MIF and one D-DT were mostly expressed in haemocytes and mantle rim of untreated animals, while D-DT paralogs often showed very limited expression, suggesting an accessory role or their persistence as relict genes.

Identification of a newly described OsHV-1 µvar from the North Adriatic Sea (Italy).

The surveillance activities for abnormal bivalve mortality events in Italy include the diagnosis of ostreid herpesvirus type 1 (OsHV-1) in symptomatic oysters. OsHV-1-positive oysters (Crassostrea gigas) were used as a source for in vivo virus propagation and a virus-rich sample was selected to perform shotgun sequencing based on Illumina technology. Starting from this unpurified supernatant sample from gills and mantle, we generated 3.5 million reads (2×300 bp) and de novo assembled the whole genome of an Italian OsHV-1 microvariant (OsHV-1-PT). The OsHV-1-PT genome encodes 125 putative ORFs, 7 of which had not previously been predicted in other sequenced Malacoherpesviridae. Overall, OsHV-1-PT displays typical microvariant OsHV-1 genome features, while few polymorphisms (0.08 %) determine its uniqueness. As little is known about the genetic determinants of OsHV-1 virulence, comparing complete OsHV-1 genomes supports a better understanding of the virus pathogenicity and provides new insights into virus-host interactions.

Dual transcriptomic analysis of Ostreid herpesvirus 1 infected Scapharca broughtonii with an emphasis on viral anti-apoptosis activities and host oxidative bursts.

The ark shell, Scapharca (Anadara) broughtonii, is an economically important marine shellfish species in Northwestern Pacific. Mass mortalities of ark shell adults related to Ostreid herpesvirus-1 (OsHV-1) infection have occurred frequently since 2012. However, due to the lack of transcriptomic resource of ark shells, the molecular mechanisms underpinning the virus-host interaction remains largely undetermined. In the present study, we resolved the dual transcriptome changes of OsHV-1 infected ark shell with Illumina sequencing. A total of 44 M sequence reads were generated, of which 67,119 reads were mapped to the OsHV-1 genome. De novo assembly of host reads resulted in 276,997 unigenes. 74,529 (26.90%), 47,653 (17.20%) and 19, 611 (7.07%) unigenes were annotated into GO, KOG and KEGG database, respectively. According to RSEM expression values, we identified 2998 differentially expressed genes (DEGs) between control and challenged groups, which included 2065 up-regulated unigenes and 933 down-regulated unigenes. Further analysis of functional pathways indicated that OsHV-1 could inhibit host cell apoptosis mainly by the up-regulation of inhibitor of apoptosis protein (IAP), and thus facilitating its successful replication. While host hemoglobins could induce oxidative burst by suppressing its peroxidase activity, and thus defense against OsHV-1 infection. Although we reported a narrow expression of the OsHV-1 genome compared to Crassostrea gigas infection, we highlighted several common viral genes highly expressed in the two hosts, suggesting an important functional role. This study offers insights into the pathogenesis mechanisms of OsHV-1 infection in bivalve mollusks of the Arcidae family.

Induced expression of cathelicidins in trout (Oncorhynchus mykiss) challenged with four different bacterial pathogens.

Cathelicidins are an important family of antimicrobial peptide effectors of innate immunity in vertebrates. Two members of this group, CATH-1 and CATH-2, have been identified and characterized in teleosts (ray-finned fish). In this study, we investigated the expression of these genes in different tissues of rainbow trout challenged with 4 different inactivated pathogens. By using qPCR, we detected a strong induction of both cath-1 and cath-2 genes within 24 hours after intraperitoneal inoculation with Lactococcus garvieae, Yersinia ruckeri, Aeromonas salmonicida, or Flavobacterium psychrophilum cells. Up to 700-fold induction of cath-2 was observed in the spleen of animals challenged with Y. ruckeri. Moreover, we found differences in the intensity and timing of gene up-regulation in the analyzed tissues. The overall results highlight the importance of cathelicidins in the immune response mechanisms of salmonids.

Mortality occurrence and pathogen detection in Crassostrea gigas and Mytilus galloprovincialis close-growing in shallow waters (Goro lagoon, Italy).

The complex interactions occurring between farmed bivalves and their potential pathogens in the circumstances of global climate changes are current matter of study, owing to the recurrent production breakdowns reported in Europe and other regions of the world. In the frame of Project FP7-KBBE-2010-4 BIVALIFE, we investigated the occurrence of mortality and potential pathogens during the Spring-Summer transition in Crassostrea gigas and Mytilus galloprovincialis cohabiting in the shallow waters of one northern Italian lagoon (Sacca di Goro, Adriatic Sea) and regarded as susceptible and resistant species, respectively. In 2011, limited bivalve mortality was detected in the open-field trial performed with 6-12 month old spat whereas subsequent trials with 2-3 month old spat produced almost complete (2012) and considerable (2013) oyster mortality. Macroscopical examination and histology excluded the presence of notifiable pathogens but, in the sampling preceding the massive oyster spat mortality of 2012, a µdeleted variant of OsHV-1 DNA was found in wide-ranging amounts in all analyzed oysters in conjunction with substantial levels of Vibrio splendidus and Vibrio aestuarianus. The large oyster spat mortality with borderline OsHV-1 positivity recorded in 2013 supports the multi-factorial etiology of the syndrome. This is the first report of a OsHV-1 (under a form interpreted as the variant µVar) in the Goro lagoon. Transcriptional host footprints are under investigation to better understand the bivalve response to environmental factors, included viral and bacterial pathogens, in relation to the observed mortalities.

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.

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.

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.