No more non-model species: the promise of next generation sequencing for comparative immunology

The future of amphibian immunology: opportunities and challenges

Historically, amphibians have been essential to our understanding of vertebrate biology and animal development. Because development from egg to tadpole to adult frog can be directly observed, amphibians contributed greatly to our understanding of not only vertebrate animal development but also the development of the immune system. The South African clawed frog (Xenopus laevis) has been key to many of these findings. For example, using Xenopus as a model, the comparative immunology community learned about the contribution of hematopoietic stem cells to development of the immune system and about the diversity of antibodies, B cells, T cells and antigen presenting cells. Amphibians offer many advantages as unique potential model systems to address questions about immune skin interactions, host responses to mycobacteria, the diverse functions of interferons, and immune and mucosal interactions. However, there are also many challenges to advance the research including the lack of specific reagents and well annotated genomes of diverse species. While much is known, many important questions remain. The aim of this short commentary is to look to the future of comparative immunology of amphibians as a group. By identifying some important questions or "information-deficit" areas of research, I hope to pique the interest of younger developing scientists and persuade funding agencies to continue to support comparative immunology studies including those of amphibians.

Transcriptional profiling of Bulinus globosus provides insights into immune gene families in snails supporting the transmission of urogenital schistosomiasis

Schistosomiasis, urogenital and intestinal, afflicts 251 million people worldwide with approximately two-thirds of the patients suffering from the urogenital form of the disease. Freshwater snails of the genus Bulinus (Gastropoda: Planorbidae) serve as obligate intermediate hosts for Schistosoma haematobium, the etiologic agent of human urogenital schistosomiasis. These snails also act as vectors for the transmission of schistosomiasis in livestock and wildlife. Despite their crucial role in human and veterinary medicine, our basic understanding at the molecular level of the entire Bulinus genus, which comprises 37 recognized species, is very limited. In this study, we employed Illumina-based RNA sequencing (RNAseq) to profile the genome-wide transcriptome of Bulinus globosus, one of the most important intermediate hosts for S. haematobium in Africa. A total of 179,221 transcripts (N50 = 1,235) were assembled and the benchmarking universal single-copy orthologs (BUSCO) was estimated to be 97.7%. The analysis revealed a substantial number of transcripts encoding evolutionarily conserved immune-related proteins, particularly C-type lectin (CLECT) domain-containing proteins (n = 316), Toll/Interleukin 1-receptor (TIR)-containing proteins (n = 75), and fibrinogen related domain-containing molecules (FReD) (n = 165). Notably, none of the FReDs are fibrinogen-related proteins (FREPs) (immunoglobulin superfamily (IgSF) + fibrinogen (FBG)). This RNAseq-based transcriptional profile provides new insights into immune capabilities of Bulinus snails, helps provide a framework to explain the complex patterns of compatibility between snails and schistosomes, and improves our overall understanding of comparative immunology.

Unveiling Immunomodulatory Effects of Euglena gracilis in Immunosuppressed Mice: Transcriptome and Pathway Analysis

The immunomodulatory effects of Euglena gracilis (Euglena) and its bioactive component, β-1,3-glucan (paramylon), have been clarified through various studies. However, the detailed mechanisms of the immune regulation remain to be elucidated. This study was designed not only to investigate the immunomodulatory effects but also to determine the genetic mechanisms of Euglena and β-glucan in cyclophosphamide (CCP)-induced immunosuppressed mice. The animals were orally administered saline, Euglena (800 mg/kg B.W.) or β-glucan (400 mg/kg B.W.) for 19 days, and CCP (80 mg/kg B.W.) was subsequently administered to induce immunosuppression in the mice. The mice exhibited significant decreases in body weight, organ weight, and the spleen index. However, there were significant improvements in the spleen weight and the spleen index in CCP-induced mice after the oral administration of Euglena and β-glucan. Transcriptome analysis of the splenocytes revealed immune-related differentially expressed genes (DEGs) regulated in the Euglena- and β-glucantreated groups. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses indicated that pathways related with interleukin (IL)-17 and cAMP play significant roles in regulating T cells, B cells, and inflammatory cytokines. Additionally, Ptgs2, a major inflammatory factor, was exclusively expressed in the Euglena-treated group, suggesting that Euglena's beneficial components, such as carotenoids, could regulate these genes by influencing immune lymphocytes and inflammatory cytokines in CCP-induced mice. This study validated the immunomodulatory effects of Euglena and highlighted its underlying mechanisms, suggesting a positive contribution to the determination of phenotypes associated with immune-related diseases and the research and development of immunotherapies.

Single-Cell RNA Sequencing Reveals Microevolution of the Stickleback Immune System

The risk and severity of pathogen infections in humans, livestock, or wild organisms depend on host immune function, which can vary between closely related host populations or even among individuals. This immune variation can entail between-population differences in immune gene coding sequences, copy number, or expression. In recent years, many studies have focused on population divergence in immunity using whole-tissue transcriptomics. But, whole-tissue transcriptomics cannot distinguish between evolved differences in gene regulation within cells, versus changes in cell composition within the focal tissue. Here, we leverage single-cell transcriptomic approaches to document signatures of microevolution of immune system structure in a natural system, the three-spined stickleback (Gasterosteus aculeatus). We sampled nine adult fish from three populations with variability in resistance to a cestode parasite, Schistocephalus solidus, to create the first comprehensive immune cell atlas for G. aculeatus. Eight broad immune cell types, corresponding to major vertebrate immune cells, were identified. We were also able to document significant variation in both abundance and expression profiles of the individual immune cell types among the three populations of fish. Furthermore, we demonstrate that identified cell type markers can be used to reinterpret traditional transcriptomic data: we reevaluate previously published whole-tissue transcriptome data from a quantitative genetic experimental infection study to gain better resolution relating infection outcomes to inferred cell type variation. Our combined study demonstrates the power of single-cell sequencing to not only document evolutionary phenomena (i.e., microevolution of immune cells) but also increase the power of traditional transcriptomic data sets.

Immune diversity in lophotrochozoans, with a focus on recognition and effector systems

Lophotrochozoa is one of the most species-rich but immunologically poorly explored phyla. Although lack of acquired response in a narrow sense, lophotrochozoans possess various genetic mechanisms that enhance the diversity and specificity of innate immune system. Here, we review the recent advances of comparative immunology studies in lophotrochozoans with focus on immune recognition and effector systems. Haemocytes and coelomocytes are general important yet understudied player. Comparative genomics studies suggest expansion and functional divergence of lophotrochozoan immune reorganization systems is not as "homogeneous and simple" as we thought including the large-scale expansion and molecular divergence of pattern recognition receptors (PRRs) (TLRs, RLRs, lectins, etc.) and signaling adapters (MyD88s etc.), significant domain recombination of immune receptors (RLR, NLRs, lectins, etc.), extensive somatic recombination of fibrinogenrelated proteins (FREPs) in snails. Furthermore, there are repeatedly identified molecular mechanisms that generate immune effector diversity, including high polymorphism of antimicrobial peptides and proteins (AMPs), reactive oxygen and nitrogen species (RONS) and cytokines. Finally, we argue that the next generation omics tools and the recently emerged genome editing technicism will revolutionize our understanding of innate immune system in a comparative immunology perspective.

Global gene expression responses of Atlantic salmon skin to Moritella viscosa

Moritella viscosa is a Gram-negative pathogen that causes large, chronic ulcers, known as winter-ulcer disease, in the skin of several fish species including Atlantic salmon. We used a bath challenge approach to profile the transcriptome responses of M. viscosa-infected Atlantic salmon skin at the lesion (Mv-At) and away from the lesion (Mv-Aw) sites. M. viscosa infection was confirmed through RNA-based qPCR assays. RNA-Seq identified 5212 and 2911 transcripts differentially expressed in the Mv-At compared to no-infection control and Mv-Aw groups, respectively. Also, there were 563 differentially expressed transcripts when comparing the Mv-Aw to control samples. Our results suggest that M. viscosa caused massive and strong, but largely infection site-focused, transcriptome dysregulations in Atlantic salmon skin, and its effects beyond the skin lesion site were comparably subtle. The M. viscosa-induced transcripts of Atlantic salmon were mainly involved in innate and adaptive immune response-related pathways, whereas the suppressed transcripts by this pathogen were largely connected to developmental and cellular processes. As validated by qPCR, M. viscosa dysregulated transcripts encoding receptors, signal transducers, transcription factors and immune effectors playing roles in TLR- and IFN-dependent pathways as well as immunoregulation, antigen presentation and T-cell development. This study broadened the current understanding of molecular pathways underlying M. viscosa-triggered responses of Atlantic salmon, and identified biomarkers that may assist to diagnose and combat this pathogen.

Computational Methods for the Discovery and Annotation of Viral Integrations

The transfer of genetic material between viruses and eukaryotic cells is pervasive. Somatic integrations of DNA viruses and retroviruses have been linked to persistent viral infection and genotoxic effects. Integrations into germline cells, referred to as Endogenous Viral Elements (EVEs), can be co-opted for host functions. Besides DNA viruses and retroviruses, EVEs can also derive from nonretroviral RNA viruses, which have often been observed in piRNA clusters. Here, we describe a bioinformatic framework to annotate EVEs in a genome assembly, study their widespread occurrence and polymorphism and identify sample-specific viral integrations using whole genome sequencing data.

A high-quality genome and comparison of short- versus long-read transcriptome of the palaearctic duck Aythya fuligula (tufted duck)

BACKGROUND

The tufted duck is a non-model organism that experiences high mortality in highly pathogenic avian influenza outbreaks. It belongs to the same bird family (Anatidae) as the mallard, one of the best-studied natural hosts of low-pathogenic avian influenza viruses. Studies in non-model bird species are crucial to disentangle the role of the host response in avian influenza virus infection in the natural reservoir. Such endeavour requires a high-quality genome assembly and transcriptome.

FINDINGS

This study presents the first high-quality, chromosome-level reference genome assembly of the tufted duck using the Vertebrate Genomes Project pipeline. We sequenced RNA (complementary DNA) from brain, ileum, lung, ovary, spleen, and testis using Illumina short-read and Pacific Biosciences long-read sequencing platforms, which were used for annotation. We found 34 autosomes plus Z and W sex chromosomes in the curated genome assembly, with 99.6% of the sequence assigned to chromosomes. Functional annotation revealed 14,099 protein-coding genes that generate 111,934 transcripts, which implies a mean of 7.9 isoforms per gene. We also identified 246 small RNA families.

CONCLUSIONS

This annotated genome contributes to continuing research into the host response in avian influenza virus infections in a natural reservoir. Our findings from a comparison between short-read and long-read reference transcriptomics contribute to a deeper understanding of these competing options. In this study, both technologies complemented each other. We expect this annotation to be a foundation for further comparative and evolutionary genomic studies, including many waterfowl relatives with differing susceptibilities to avian influenza viruses.

Between-population differences in constitutive and infection-induced gene expression in threespine stickleback

Vertebrate immunity is a complex system consisting of a mix of constitutive and inducible defences. Furthermore, host immunity is subject to selective pressure from a range of parasites and pathogens which can produce variation in these defences across populations. As populations evolve immune responses to parasites, they may adapt via a combination of (1) constitutive differences, (2) shared inducible responses, or (3) divergent inducible responses. Here, we leverage a powerful natural host‐parasite model system (Gasterosteus aculeatus and Schistochephalus solidus) to tease apart the relative contributions of these three types of adaptations to among‐population divergence in response to parasites. Gene expression analyses revealed limited evidence of significant divergence in constitutive expression of immune defence, and strong signatures of conserved inducible responses to the parasite. Furthermore, our results highlight a handful of immune‐related genes which show divergent inducible responses which may contribute disproportionately to functional differences in infection success or failure. In addition to investigating variation in evolutionary adaptation to parasite selection, we also leverage this unique data set to improve understanding of cellular mechanisms underlying a putative resistance phenotype (fibrosis). Combined, our results provide a case study in evolutionary immunology showing that a very small number of genes may contribute to genotype differences in infection response.

Immune Gene Expression Covaries with Gut Microbiome Composition in Stickleback

Commensal microbial communities have immense effects on their vertebrate hosts, contributing to a number of physiological functions, as well as host fitness. In particular, host immunity is strongly linked to microbiota composition through poorly understood bi-directional links. Gene expression may be a potential mediator of these links between microbial communities and host function. However, few studies have investigated connections between microbiota composition and expression of host immune genes in complex systems. Here, we leverage a large study of laboratory-raised fish from the species Gasterosteus aculeatus (three-spined stickleback) to document correlations between gene expression and microbiome composition. First, we examined correlations between microbiome alpha diversity and gene expression. Our results demonstrate robust positive associations between microbial alpha diversity and expression of host immune genes. Next, we examined correlations between host gene expression and abundance of microbial taxa. We identified 15 microbial families that were highly correlated with host gene expression. These families were all tightly correlated with host expression of immune genes and processes, falling into one of three categories—those positively correlated, negatively correlated, and neutrally related to immune processes. Furthermore, we highlight several important immune processes that are commonly associated with the abundance of these taxa, including both macrophage and B cell functions. Further functional characterization of microbial taxa will help disentangle the mechanisms of the correlations described here. In sum, our study supports prevailing hypotheses of intimate links between host immunity and gut microbiome composition.

Examining adaptive evolution of immune activity: opportunities provided by gastropods in the age of 'omics'

Parasites threaten all free-living organisms, including molluscs. Understanding the evolution of immune defence traits in natural host populations is crucial for predicting their long-term performance under continuous infection risk. Adaptive trait evolution requires that traits are subject to selection (i.e. contribute to organismal fitness) and that they are heritable. Despite broad interest in the evolutionary ecology of immune activity in animals, the understanding of selection on and evolutionary potential of immune defence traits is far from comprehensive. For instance, empirical observations are only rarely in line with theoretical predictions of immune activity being subject to stabilizing selection. This discrepancy may be because ecoimmunological studies can typically cover only a fraction of the complexity of an animal immune system. Similarly, molecular immunology/immunogenetics studies provide a mechanistic understanding of immunity, but neglect variation that arises from natural genetic differences among individuals and from environmental conditions. Here, we review the current literature on natural selection on and evolutionary potential of immune traits in animals, signal how merging ecological immunology and genomics will strengthen evolutionary ecological research on immunity, and indicate research opportunities for molluscan gastropods for which well-established ecological understanding and/or 'immune-omics' resources are already available. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

Transcriptome profiling of Lymnaea stagnalis (Gastropoda) for ecoimmunological research

BACKGROUND

Host immune function can contribute to numerous ecological/evolutionary processes. Ecoimmunological studies, however, typically use one/few phenotypic immune assays and thus do not consider the complexity of the immune system. Therefore, "omics" resources that allow quantifying immune activity across multiple pathways are needed for ecoimmunological models. We applied short-read based RNAseq (Illumina NextSeq 500, PE-81) to characterise transcriptome profiles of Lymnaea stagnalis (Gastropoda), a multipurpose model snail species. We used a genetically diverse snail stock and exposed individuals to immune elicitors (injury, bacterial/trematode pathogens) and changes in environmental conditions that can alter immune activity (temperature, food availability).

RESULTS

Immune defence factors identified in the de novo assembly covered elements broadly described in other gastropods. For instance, pathogen-recognition receptors (PRR) and lectins activate Toll-like receptor (TLR) pathway and cytokines that regulate cellular and humoral defences. Surprisingly, only modest diversity of antimicrobial peptides and fibrinogen related proteins were detected when compared with other taxa. Additionally, multiple defence factors that may contribute to the phenotypic immune assays used to quantify antibacterial activity and phenoloxidase (PO)/melanisation-type reaction in this species were found. Experimental treatments revealed factors from non-self recognition (lectins) and signalling (TLR pathway, cytokines) to effectors (e.g., antibacterial proteins, PO enzymes) whose transcription depended on immune stimuli and environmental conditions, as well as components of snail physiology/metabolism that may drive these effects. Interestingly, the transcription of many factors (e.g., PRR, lectins, cytokines, PO enzymes, antibacterial proteins) showed high among-individual variation.

CONCLUSIONS

Our results indicate several uniform aspects of gastropod immunity, but also apparent differences between L. stagnalis and some previously examined taxa. Interestingly, in addition to immune defence factors that responded to immune elicitors and changes in environmental conditions, many factors showed high among-individual variation across experimental snails. We propose that such factors are highly important to be included in future ecoimmunological studies because they may be the key determinants of differences in parasite resistance among individuals both within and between natural snail populations.

Physiology and immunity of the invasive giant African snail, Achatina (Lissachatina) fulica, intermediate host of Angiostrongylus cantonensis

As one of the most successful invasive land snail species, Achatina (Lissachatina) fulica Bowdich, 1822 has achieved wide global distribution, particularly in (sub)tropical regions, with further dispersal likely due to climate change. This species of giant African snails (up to 17 cm shell length) is a pest that has extensive negative impact on agriculture and can serve as vector for several parasites, including Angiostrongylus cantonensis, a nematode parasite that causes (human) eosinophilic meningitis, an emergent disease. Investigation showed that A. cantonensis infection negatively impacts the metabolism of A. fulica by depleting polysaccharide stores of the intermediate host, compromising the energy balance of the snail. A review of the literature indicates that A. fulica possesses potent innate type immune defenses to counter infection, including phagocytic hemocytes capable of deploying reactive oxygen species and lectins for non-self recognition, a serine protease-dependent coagulation response (not observed in other taxa of gastropods), as well as antimicrobial proteins including achacin, an antimicrobial protein. A recent chromosome level genome assembly will facilitate progressively detailed characterization of these immune features of A. fulica. We strongly encourage further immunological studies of A. fulica, ranging from organismal level to molecular biology to gain better understanding of the A. fulica internal defense response to nematode pathogens like A. cantonensis and the contribution of immune function to the invasiveness of (snail) species. Characterization of immunity of A. fulica, representing the understudied Stylommatophora (panpulmonate landsnails) will also broaden the comparative immunology of Gastropoda.

Snail defence responses to parasite infection: The Lymnaea stagnalis-Trichobilharzia szidati model

Lymnaea stagnalis is a common freshwater gastropod. Importantly, the snail serves as the intermediate host for more than one hundred species of digenetic trematodes, including the avian schistosome Trichobilharzia szidati, a causative agent of cercarial dermatitis in humans. Infection of L. stagnalis by T. szidati initiates a dynamic confrontation between the host and the parasite that culminates in immunocompatibility ensuring survival and development of larvae. Unfortunately, the molecular mechanisms determining this immunocompatibility remain poorly characterised. By employing a variety of immune elicitors, including chemical compounds, PAMPs and bacteria, research in the last two decades has elucidated some of the molecular processes that regulate the snail internal defence response such as haemocyte signalling pathways. These discoveries provide a framework for future studies of molecular interactions between T. szidati and L. stagnalis to help elucidate factors and mechanisms enabling transmission of schistosome parasites. Moreover, support from recently available next generation sequence data and CRISPR-enabled functional genomics should further enable L. stagnalis as an important model for comparative immunology and contribute to a more comprehensive understanding of immune functions in gastropod molluscs.

Transcriptome-wide analysis of immune responses in Eriocheir sinensis hemocytes after challenge with different microbial derivatives

Knowledge about how Eriocheir sinensis interacts with microorganisms in its ambient environment is still lacking. Using RNA-Seq, we determined the most conserved genes and pathways compared with other animals and detected highly-induced immune genes in E. sinensis hemocytes post-in vivo challenge with different microbial derivatives. In total, 33.2 million high-quality reads were generated and assembled into 177,679 contigs. Completeness assessment and functional annotation were performed. Lipopolysaccharide, peptidoglycan, and β-1, 3-glucan stimulation induced 373, 173, and 108 differentially expressed (DE) transcripts, respectively. GO terms such as 'G-protein-coupled receptor binding', 'negative regulation of mitogen-activated protein kinase activity', and 'positive regulation of blood circulation' were enriched in the DE transcripts. Quantitative real-time PCR validated the data for selected genes. Our data contribute to understanding the immune defense mechanism in E. sinensis and the development of the innate immune system, thereby providing insights into disease control and prevention in aquaculture.

Transcriptomic analysis of oligochaete immune responses to myxosporeans infection: Branchiura sowerbyi infected with Myxobolus cultus

The Myxozoa are endoparasites characterized by a two-host life cycle that typically involves invertebrates and vertebrates as definitive and intermediate hosts, respectively. However, little is known about invertebrate-myxosporean interactions, particularly about patterns of host immune defense. We used RNA-sequencing to identify genes that are possibly involved in the immune responses of the oligochaete Branchiura sowerbyi naturally infected with Myxobolus cultus. De novo assembly of the B. sowerbyi transcriptome yielded 119,031 unigenes, with an average length of 896 bp and an N50 length of 1754 bp. Comparative transcriptome analysis revealed 4059 differentially expressed genes (DEGs) between M. cultus-infected and uninfected B. sowerbyi groups, including 3802 upregulated genes and 257 downregulated genes. Among the B. sowerbyi immune factors implicated in the responses to M. cultus infection, DEGs related to lectins, ubiquitin-mediated proteolysis, phagocytosis, oxidative-antioxidative responses, proteases, and protease inhibitors were upregulated. The expression of some immune-related molecules such as calmodulin, heat shock proteins, antimicrobial peptides, lysenin, and serum amyoid A protein were also significantly upregulated. The expression patterns of 14 immune-related DEGs identified by RNA-seq were validated by quantitative real-time polymerase chain reaction. This study is the first attempt to characterize the B. sowerbyi transcriptome and identify immune-related molecules possibly associated with M. cultus infection. It is also the first report of invertebrate host-myxosporean interactions at the transcriptomic level. Our results will facilitate the elucidation of adaptive evolution mechanisms of myxosporean parasites in the definitive host and the genetic basis for differences in resistance of invertebrate hosts of different genotypes to a myxosporean species.

De novo assembly, characterization of tissue-specific transcriptomes and identification of immune related genes from the scallop Argopecten purpuratus

The scallop Argopecten purpuratus is one of the most economically important cultured mollusks on the coasts from Chile and Peru but its production has declined, in part, due to the emergence of mass mortality events of unknown origin. Driven by this scenario, increasing progress has been made in recent years in the comprehension of immune response mechanisms in this species. However, it is still not entirely understood how different mucosal interfaces participate and cooperate with the immune competent cells, the hemocytes, in the immune defense. Thus, in this work we aimed to characterize the transcriptome of three tissues with immune relevance from A. purpuratus by next-generation sequencing and de novo transcriptome assembly. For this, 18 cDNA libraries were constructed from digestive gland, gills and hemocytes tissues of scallops from different immune conditions and sequenced by the Illumina HiSeq4000 platform. A total of 967.964.884 raw reads were obtained and 967.432.652 clean reads were generated. The clean reads were de novo assembled into 46.601 high quality contigs and 32.299 (69.31%) contigs were subsequently annotated. In addition, three de novo specific assemblies were performed from clean reads obtained from each tissue cDNA libraries for their comparison. Gene ontology (GO) and KEGG analyses revealed that annotated sequences from digestive gland, gills and hemocytes could be classified into both general and specific subcategory terms and known biological pathways, respectively, according to the tissue nature. Finally, several immune related candidate genes were identified, and the differential expression of tissue-specific genes was established, suggesting they could display specific roles in the host defense. The data presented in this study provide the first insight into the tissue specific transcriptome profiles of A. purpuratus, which should be considered for further research on the interplay between the hemocytes and mucosal immune responses.

Transcriptome analysis and discovery of genes involved in immune pathways in Solen strictus (Gould, 1861) under Vibrio anguillarum

The Solen strictus (Gould, 1861) has been recognized as an important marine economic bivalve in Eastern and Southeast Asia. To gain a better understanding of the S. strictus immune system and its related genes in response to bacterial infections, we performed a comparative gene transcription analysis from S. strictus with Vibrio anguillarum through RNA-Seq technology, meanwhile the differentially expressed genes (DEGs) were investigated. After assembly, a total of 195,774 transcripts with an average length of 996 bp were obtained. Total 153,038 unigenes were annotated in the nr, Swiss-Prot, KEGG, COG, KOG, GO and Pfam databases, and 56,597 unigenes (36.98%) were annotated in at least one database. After bacterial challenge, there were 1588 significant differentially expressed genes (DEGs) between the challenged and control groups, including 999 up-regulated and 589 down-regulated genes. All the DEGs were classified into three gene ontology categories, and allocated to 225 KEGG pathways. Immune-related genes were detected from immune system pathways among the top 20 enriched pathways, such as Toll-like receptor signaling, RIG-I-like receptor signaling and NOD-like receptor signaling pathway. In addition, 56,079 potential simple sequence repeats (SSRs) and 1,031,521 candidate single nucleotide polymorphisms (SNPs) were detected and identified in the S. strictus transcriptome. Results of the present study will provide valuable theoretical resources for future genetic and genomic research on S. strictus. The research results will be helpful for improving the efficiency and quality of artificial breeding, establishing genetic linkage map, and enhancing health management for this species.

Profiling of the TCRβ repertoire in non-model species using high-throughput sequencing

In recent years, immune repertoire profiling with high-throughput sequencing (HTS) has advanced our understanding of adaptive immunity. However, fast progress in the field applied mostly to human and mouse research, with only few studies devoted to other model vertebrates. We present the first in-depth characterization of the T-cell receptor (TCR) repertoire in a non-model mammal (bank vole, Myodes glareolus), widely used in ecological and evolutionary research. We used RNA from spleens, 5′RACE and HTS to describe V and J segments of TCRβ, qualitatively characterize preferential V–J segment usage and CDR3 length distribution. Overall orthology to murine genes was preserved, with 11 J and 37 V genes found in voles (although 3 V genes lacked a close orthologue). Further, we implemented unique molecular identifiers for quantitative analysis of CDR3 repertoire with stringent error correction. A conservative, lower bound estimation of the TCRβ repertoire was similar to that found for mice (1.7–2.3 × 10⁵ clonotypes). We hope that by providing an easy-to-follow molecular protocol and on-line bioinformatics tools that do not require reference sequences (AmpliTCR and AmpliCDR3), we will encourage HTS immune repertoire profiling in other non-model vertebrates, thus opening new research avenues in e.g. comparative immunology, ecology and evolutionary biology.

Comparative immunological study of the snail Physella acuta (Hygrophila, Pulmonata) reveals shared and unique aspects of gastropod immunobiology

The freshwater snail Physella acuta was selected to expand the perspective of comparative snail immunology. Analysis of Physella acuta, belonging to the Physidae, taxonomic sister family to Planorbidae, affords family-level comparison of immune features characterized from Biomphalaria glabrata, the model snail often used to interpret general gastropod immunity. To capture constitutive and induced immune sequences, transcriptomes of an individual Physella acuta snail, 12 h post injection with bacteria (Gram -/+) and one sham-exposed snail were recorded with 454 pyrosequencing. Assembly yielded a combined reference transcriptome containing 24,288 transcripts. Additionally, genomic Illumina reads were obtained (∼15-fold coverage). Recovery of transcripts for two macin-like antimicrobial peptides (AMPs), 12 aplysianins, four LBP/BPIs and three physalysins indicated that Physella acuta shares a similar organization of antimicrobial defenses with Biomphalaria glabrata, contrasting a modest AMP arsenal with a diverse set of antimicrobial proteins. The lack of predicted transmembrane domains in all seven Physella acuta PGRP transcripts supports the notion that gastropods do not employ cell-bound PGRP receptors, different from ecdysozoan invertebrates yet similar to mammals (vertebrate deuterostomes). The well-documented sequence diversification by Biomphalaria glabrata FREPs (immune lectins comprising immunoglobulin superfamily domains and fibrinogen domains), resulting from somatic mutations of a large FREP gene family is hypothesized to be unique to Planorbidae; Physella acuta revealed just two bonafide FREP genes and these were not diversified. Furthermore, the flatworm parasite Echinostoma paraensei, confirmed here to infect both snail species, did not evoke from Physella acuta the abundant expression of FREP proteins at 2, 4 and 8 days post exposure that was previously observed from Biomphalaria glabrata. The Physella acuta reference transcriptome also revealed 24 unique transcripts encoding proteins consisting of a single fibrinogen-related domain (FReDs), with a short N-terminal sequence encoding either a signal peptide, transmembrane domain or no predicted features. The Physella acuta FReDs are candidate immune genes based on implication of similar sequences in immunity of bivalve molluscs. Overall, comparative analysis of snails of sister families elucidated the potential for taxon-specific immune features and investigation of strategically selected species will provide a more comprehensive view of gastropod immunity.

Different milk diets have substantial effects on the jejunal mucosal immune system of pre-weaning calves, as demonstrated by whole transcriptome sequencing

There is increasing evidence that nutrition during early mammalian life has a strong influence on health and performance in later life. However, there are conflicting data concerning the appropriate milk diet. This discrepancy particularly applies to ruminants, a group of mammals that switch from monogastric status to rumination during weaning. Little is known regarding how the whole genome expression pattern in the juvenile ruminant gut is affected by alternative milk diets. Thus, we performed a next-generation-sequencing-based holistic whole transcriptome analysis of the jejunum in male pre-weaned German Holstein calves fed diets with restricted or unlimited access to milk during the first 8 weeks of life. Both groups were provided hay and concentrate ad libitum. The analysis of jejunal mucosa samples collected 80 days after birth and four weeks after the end of the feeding regimes revealed 275 differentially expressed loci. While the differentially expressed loci comprised 67 genes encoding proteins relevant to metabolism or metabolic adaptation, the most distinct difference between the two groups was the consistently lower activation of the immune system in calves that experienced restricted milk access compared to calves fed milk ad libitum. In conclusion, different early life milk diets had significant prolonged effects on the intestinal immune system.

Transcriptome analysis and discovery of genes involved in immune pathways from coelomocytes of Onchidium struma after bacterial challenge

Onchidium struma widely distributes in subtidal and low-tidal zones, which is considered to be an economical species with rich nutrition, a valuable biomonitor for heavy metal pollution and a representative species for evolution from ocean to land. However, there is limited genetic information available for O. struma development. This study compared transcriptomic profiles of coelomocytes from normal and bacteria infected O. struma by Illumina-based paired-end sequencing to explore the molecular immune mechanism of O. struma against bacterial infection. After assembly, a total of 92,450 unigenes with an average length of 1019 bp were obtained. Approximately 34,964 (37.82%) unigenes were annotated in the Nr NCBI database and 40.1% of unigenes were similar with that of Aplysia californica. Among them, 7609 unigenes were classified into three Gene Ontology (GO) categories: biological process (3250 unigenes, 42.7%), cellular component (2,281, 30.0%) and molecular function (2078 unigenes, 27.3%). A total of 22,776 unigenes were aligned to the Clusters of Orthologous Groups (COG) of proteins and classified into 25 functional categories. Following bacterial infection, 10,623 differently expressed unigenes (DEGs) were identified, including 7644 up-regulated and 2979 down-regulated unigenes. Further KEGG analysis annotated 11,681 DEGs to 42 pathways, and 11 pathways were identified to be related with diseases and immune system. To our knowledge, it was first time to analyze transcriptome profiles of O. struma. Results of the present study will provide valuable theoretical resources for future genetic and genomic research on O. struma. The research results will be helpful for improving the efficiency and quality of artificial breeding, establishing genetic linkage map, and enhancing health management for this species.

Evolutionary dynamics of an expressed MHC class IIβ locus in the Ranidae (Anura) uncovered by genome walking and high-throughput amplicon sequencing

The Major Histocompatibility Complex (MHC) is a genomic region encoding immune loci that are important and frequently used markers in studies of adaptive genetic variation and disease resistance. Given the primary role of infectious diseases in contributing to global amphibian declines, we characterized the hypervariable exon 2 and flanking introns of the MHC Class IIβ chain for 17 species of frogs in the Ranidae, a speciose and cosmopolitan family facing widespread pathogen infections and declines. We find high levels of genetic variation concentrated in the Peptide Binding Region (PBR) of the exon. Ten codons are under positive selection, nine of which are located in the mammal-defined PBR. We hypothesize that the tenth codon (residue 21) is an amphibian-specific PBR site that may be important in disease resistance. Trans-species and trans-generic polymorphisms are evident from exon-based genealogies, and co-phylogenetic analyses between intron, exon and mitochondrial based reconstructions reveal incongruent topologies, likely due to different locus histories. We developed two sets of barcoded adapters that reliably amplify a single and likely functional locus in all screened species using both 454 and Illumina based sequencing methods. These primers provide a resource for multiplexing and directly sequencing hundreds of samples in a single sequencing run, avoiding the labour and chimeric sequences associated with cloning, and enabling MHC population genetic analyses. Although the primers are currently limited to the 17 species we tested, these sequences and protocols provide a useful genetic resource and can serve as a starting point for future disease, adaptation and conservation studies across a range of anuran taxa.

Integrated multi-omic analyses in Biomphalaria-Schistosoma dialogue reveal the immunobiological significance of FREP-SmPoMuc interaction

The fresh water snail Biomphalaria glabrata is one of the vectors of the trematode pathogen Schistosoma mansoni, which is one of the agents responsible of human schistosomiasis. In this host-parasite interaction, co-evolutionary dynamic results into an infectivity mosaic known as compatibility polymorphism. Integrative approaches including large scale molecular approaches have been conducted in recent years to improve our understanding of the mechanisms underlying compatibility. This review presents the combination of integrated Multi-Omic approaches leading to the discovery of two repertoires of polymorphic and/or diversified interacting molecules: the parasite antigens S. mansoni polymorphic mucins (SmPoMucs) and the B. glabrata immune receptors fibrinogen-related proteins (FREPs). We argue that their interactions may be major components for defining the compatible/incompatible status of a specific snail/schistosome combination.

Genomic and transcriptomic approaches to study immunology in cyprinids: What is next?

Accelerated by the introduction of Next-Generation Sequencing (NGS), a number of genomes of cyprinid fish species have been drafted, leading to a highly valuable collective resource of comparative genome information on cyprinids (Cyprinidae). In addition, NGS-based transcriptome analyses of different developmental stages, organs, or cell types, increasingly contribute to the understanding of complex physiological processes, including immune responses. Cyprinids are a highly interesting family because they comprise one of the most-diversified families of teleosts and because of their variation in ploidy level, with diploid, triploid, tetraploid, hexaploid and sometimes even octoploid species. The wealth of data obtained from NGS technologies provides both challenges and opportunities for immunological research, which will be discussed here. Correct interpretation of ploidy effects on immune responses requires knowledge of the degree of functional divergence between duplicated genes, which can differ even between closely-related cyprinid fish species. We summarize NGS-based progress in analysing immune responses and discuss the importance of respecting the presence of (multiple) duplicated gene sequences when performing transcriptome analyses for detailed understanding of complex physiological processes. Progressively, advances in NGS technology are providing workable methods to further elucidate the implications of gene duplication events and functional divergence of duplicates genes and proteins involved in immune responses in cyprinids. We conclude with discussing how future applications of NGS technologies and analysis methods could enhance immunological research and understanding.

Describing the diversity of Ag specific receptors in vertebrates: Contribution of repertoire deep sequencing

During the last decades, gene and cDNA cloning identified TCR and Ig genes across vertebrates; genome sequencing of TCR and Ig loci in many species revealed the different organizations selected during evolution under the pressure of generating diverse repertoires of Ag receptors. By detecting clonotypes over a wide range of frequency, deep sequencing of Ig and TCR transcripts provides a new way to compare the structure of expressed repertoires in species of various sizes, at different stages of development, with different physiologies, and displaying multiple adaptations to the environment. In this review, we provide a short overview of the technologies currently used to produce global description of immune repertoires, describe how they have already been used in comparative immunology, and we discuss the future potential of such approaches. The development of these methodologies in new species holds promise for new discoveries concerning particular adaptations. As an example, understanding the development of adaptive immunity across metamorphosis in frogs has been made possible by such approaches. Repertoire sequencing is now widely used, not only in basic research but also in the context of immunotherapy and vaccination. Analysis of fish responses to pathogens and vaccines has already benefited from these methods. Finally, we also discuss potential advances based on repertoire sequencing of multigene families of immune sensors and effectors in invertebrates.

Omics and cytokine discovery in fish: Presenting the Yellowtail kingfish (Seriola lalandi) as a case study

A continued programme of research is essential to overcome production bottlenecks in any aquacultured fish species. Since the introduction of genetic and molecular techniques, the quality of immune research undertaken in fish has greatly improved. Thousands of species specific cytokine genes have been discovered, which can be used to conduct more sensitive studies to understand how fish physiology is affected by aquaculture environments or disease. Newly available transcriptomic technologies, make it increasingly easier to study the immunogenetics of farmed species for which little data exists. This paper reviews how the application of transcriptomic procedures such as RNA Sequencing (RNA-Seq) can advance fish research. As a case study, we present some preliminary findings using RNA-Seq to identify cytokine related genes in Seriola lalandi. These will allow in-depth investigations to understand the immune responses of these fish in response to environmental change or disease and help in the development of therapeutic approaches.

Comparative immunogenomics of molluscs

Comparative immunology, studying both vertebrates and invertebrates, provided the earliest descriptions of phagocytosis as a general immune mechanism. However, the large scale of animal diversity challenges all-inclusive investigations and the field of immunology has developed by mostly emphasizing study of a few vertebrate species. In addressing the lack of comprehensive understanding of animal immunity, especially that of invertebrates, comparative immunology helps toward management of invertebrates that are food sources, agricultural pests, pathogens, or transmit diseases, and helps interpret the evolution of animal immunity. Initial studies showed that the Mollusca (second largest animal phylum), and invertebrates in general, possess innate defenses but lack the lymphocytic immune system that characterizes vertebrate immunology. Recognizing the reality of both common and taxon-specific immune features, and applying up-to-date cell and molecular research capabilities, in-depth studies of a select number of bivalve and gastropod species continue to reveal novel aspects of molluscan immunity. The genomics era heralded a new stage of comparative immunology; large-scale efforts yielded an initial set of full molluscan genome sequences that is available for analyses of full complements of immune genes and regulatory sequences. Next-generation sequencing (NGS), due to lower cost and effort required, allows individual researchers to generate large sequence datasets for growing numbers of molluscs. RNAseq provides expression profiles that enable discovery of immune genes and genome sequences reveal distribution and diversity of immune factors across molluscan phylogeny. Although computational de novo sequence assembly will benefit from continued development and automated annotation may require some experimental validation, NGS is a powerful tool for comparative immunology, especially increasing coverage of the extensive molluscan diversity. To date, immunogenomics revealed new levels of complexity of molluscan defense by indicating sequence heterogeneity in individual snails and bivalves, and members of expanded immune gene families are expressed differentially to generate pathogen-specific defense responses.

De novo transcriptome assembly and developmental mode specific gene expression of Pygospio elegans

Species with multiple different larval developmental modes are interesting models for the study of mechanisms underlying developmental mode transitions and life history evolution. Pygospio elegans, a small, tube-dwelling polychaete worm commonly found in estuarine and marine habitats around the northern hemisphere, is one species with variable developmental modes. To provide new genomic resources for studying P. elegans and to address the differences in gene expression between individuals producing offspring with different larval developmental modes, we performed whole transcriptome Illumina RNA sequencing of adult worms from two populations and prepared a de novo assembly of the P. elegans transcriptome. The transcriptome comprises 66,233 unigenes, of which 33,807 contain predicted coding sequences, 26,448 have at least one functional annotation, and 3,076 are classified as putative long non-coding RNAs. We found more than 8,000 unigenes significantly differentially expressed between adult worms from populations producing either planktonic or benthic larvae. This comprehensive transcriptome resource for P. elegans adds to the available genomic data for annelids and can be used to uncover mechanisms allowing developmental variation in this and potentially other marine invertebrate species.

Transcriptome Analysis and Identification of Differentially Expressed Transcripts of Immune-Related Genes in Spleen of Gosling and Adult Goose

The goose (Anser cygnoides), having high nutritional value, high-quality feathers and high economic benefit, is an economically important poultry species. However, the molecular mechanisms underlying the higher susceptibility to pathogens in goslings than in adult geese remains poorly understood. In this study, the histological sections of spleen tissue from a two-week-old gosling and an adult goose, respectively, were subjected to comparative analysis. The spleen of gosling was mainly composed of mesenchyma, accompanied by scattered lymphocytes, whereas the spleen parenchyma was well developed in the adult goose. To investigate goose immune-related genes, we performed deep transcriptome and gene expression analyses of the spleen samples using paired-end sequencing technology (Illumina). In total, 50,390 unigenes were assembled using Trinity software and TGICL software. Moreover, these assembled unigenes were annotated with gene descriptions and gene ontology (GO) analysis was performed. Through Kyoto encyclopedia of genes and genomes (KEGG) analysis, we investigated 558 important immune-relevant unigenes and 23 predicted cytokines. In addition, 22 immune-related genes with differential expression between gosling and adult goose were identified, among which the three genes showing largest differences in expression were immunoglobulin alpha heavy chain (IgH), mannan-binding lectin serine protease 1 isoform X1 (MASP1) and C-X-C chemokine receptor type 4 (CXCR4). Finally, of these 22 differentially expressed immune-related genes, seven genes, including tumor necrosis factor receptor superfamily member 13B (TNFRSF13B), C-C motif chemokine 4-like (CCL4), CXCR4, interleukin 2 receptor alpha (IL2RA), MHC class I heavy chain (MHCIα), transporter of antigen processing 2 (TAP2) IgH, were confirmed by quantitative real-time PCR (qRT-PCR). The expression levels of all the candidate unigenes were up-regulated in adult geese other than that of TNFRSF13B. The comparative analysis of the spleen transcriptomes of gosling and adult goose may promote better understanding of immune molecular development in goose.

The role of fibrinogen-related proteins in the gastropod immune response

Fibrinogen-related proteins or FREPs constitute a large family of molecules, defined by the presence of a fibrinogen-related domain (FReD). These molecules are found in all animals and are diverse in both form and function. Here, we review the current understanding of gastropod FREPs, which are characterized by the presence of a fibrinogen domain connected to one or two immunoglobulin superfamily domains by way of a short interceding region. We present a historical perspective on the discovery of FREPs in gastropods followed by a summary of advances made in the nearly two decades of research focused on the characterization of FREPs in Biomphalaria glabrata (BgFREPs). Topics covered include BgFREP genomic architecture, predicted structure and known functions, structural comparisons between BgFREPs, and evidence of somatic diversification. Also examined are the expression patterns of BgFREPs during snail development and immunological challenges. Recent functional characterization of the role BgFREPs play in the defence response against digenean trematodes is also presented, as well as new data investigating the nucleotide-level genomic conservation of FREPs among Pulmonate gastropods. Finally, we identify areas in need of further research. These include confirming and identifying the specific binding targets of BgFREPs and elucidating how they later engage snail haemocytes to elicit an immunological response, precise mechanisms and importance of BgFREP diversification, characterizing the tissue expression patterns of BgFREPs, as well as addressing whether gastropod FREPs retain immunological importance in alternative snail-trematode associations or more broadly in snail-pathogen interactions.

The early stages of the immune response of the European abalone Haliotis tuberculata to a Vibrio harveyi infection

Vibrio harveyi is a marine bacterial pathogen responsible for episodic abalone mortalities in France, Japan and Australia. In the European abalone, V. harveyi invades the circulatory system in a few hours after exposure and is lethal after 2 days of infection. In this study, we investigated the responses of European abalone immune cells over the first 24 h of infection. Results revealed an initial induction of immune gene expression including Rel/NF-kB, Mpeg and Clathrin. It is rapidly followed by a significant immuno-suppression characterized by reduced cellular hemocyte parameters, immune response gene expressions and enzymatic activities. Interestingly, Ferritin was overexpressed after 24 h of infection suggesting that abalone attempt to counter V. harveyi infection using soluble effectors. Immune function alteration was positively correlated with V. harveyi concentration. This study provides the evidence that V. harveyi has a hemolytic activity and an immuno-suppressive effect in the European abalone.

The use of -omic tools in the study of disease processes in marine bivalve mollusks

Our understanding of disease processes and host-pathogen interactions in model species has benefited greatly from the application of medium and high-throughput genomic, metagenomic, epigenomic, transcriptomic, and proteomic analyses. The rate at which new, low-cost, high-throughput -omic technologies are being developed has also led to an expansion in the number of studies aimed at gaining a better understanding of disease processes in bivalves. This review provides a catalogue of the genetic and -omic tools available for bivalve species and examples of how -omics has contributed to the advancement of marine bivalve disease research, with a special focus in the areas of immunity, bivalve-pathogen interactions, mechanisms of disease resistance and pathogen virulence, and disease diagnosis. The analysis of bivalve genomes and transcriptomes has revealed that many immune and stress-related gene families are expanded in the bivalve taxa examined thus far. In addition, the analysis of proteomes confirms that responses to infection are influenced by epigenetic, post-transcriptional, and post-translational modifications. The few studies performed in bivalves show that epigenetic modifications are non-random, suggesting a role for epigenetics in regulating the interactions between bivalves and their environments. Despite the progress -omic tools have enabled in the field of marine bivalve disease processes, there is much more work to be done. To date, only three bivalve genomes have been sequenced completely, with assembly status at different levels of completion. Transcriptome datasets are relatively easy and inexpensive to generate, but their interpretation will benefit greatly from high quality genome assemblies and improved data analysis pipelines. Finally, metagenomic, epigenomic, proteomic, and metabolomic studies focused on bivalve disease processes are currently limited but their expansion should be facilitated as more transcriptome datasets and complete genome sequences become available for marine bivalve species.

De Novo transcriptome analysis of Oncomelania hupensis after molluscicide treatment by next-generation sequencing: implications for biology and future snail interventions

The freshwater snail Oncomelania hupensis is the only intermediate host of Schistosoma japonicum, which causes schistosomiasis. This disease is endemic in the Far East, especially in mainland China. Because niclosamide is the only molluscicide recommended by the World Health Organization, 50% wettable powder of niclosamide ethanolamine salt (WPN), the only chemical molluscicide available in China, has been widely used as the main snail control method for over two decades. Recently, a novel molluscicide derived from niclosamide, the salt of quinoid-2',5-dichloro-4'-nitro-salicylanilide (Liu Dai Shui Yang An, LDS), has been developed and proven to have the same molluscicidal effect as WPN, with lower cost and significantly lower toxicity to fish than WPN. The mechanism by which these molluscicides cause snail death is not known. Here, we report the next-generation transcriptome sequencing of O. hupensis; 145,008,667 clean reads were generated and assembled into 254,286 unigenes. Using GO and KEGG databases, 14,860 unigenes were assigned GO annotations and 4,686 unigenes were mapped to 250 KEGG pathways. Many sequences involved in key processes associated with biological regulation and innate immunity have been identified. After the snails were exposed to LDS and WPN, 254 unigenes showed significant differential expression. These genes were shown to be involved in cell structure defects and the inhibition of neurohumoral transmission and energy metabolism, which may cause snail death. Gene expression patterns differed after exposure to LDS and WPN, and these differences must be elucidated by the identification and annotation of these unknown unigenes. We believe that this first large-scale transcriptome dataset for O. hupensis will provide an opportunity for the in-depth analysis of this biomedically important freshwater snail at the molecular level and accelerate studies of the O. hupensis genome. The data elucidating the molluscicidal mechanism will be of great benefit in future snail control efforts.

Diversity of animal immune receptors and the origins of recognition complexity in the deuterostomes

Invertebrate animals are characterized by extraordinary diversity in terms of body plan, life history and life span. The past impression that invertebrate immune responses are controlled by relatively simple innate systems is increasingly contradicted by genomic analyses that reveal significant evolutionary novelty and complexity. One accessible measure of this complexity is the multiplicity of genes encoding homologs of pattern recognition receptors. These multigene families vary significantly in size, and their sequence character suggests that they vary in function. At the same time, certain aspects of downstream signaling appear to be conserved. Here, we analyze five major classes of immune recognition receptors from newly available animal genome sequences. These include the Toll-like receptors (TLR), Nod-like receptors (NLR), SRCR domain scavenger receptors, peptidoglycan recognition proteins (PGRP), and Gram negative binding proteins (GNBP). We discuss innate immune complexity in the invertebrate deuterostomes, which was first recognized in sea urchins, within the wider context of emerging genomic information across animal phyla.

Lectin-like molecules in transcriptome of Littorina littorea hemocytes

The common periwinkle Littorina littorea was introduced in the list of models for comparative immunobiology as a representative of phylogenetically important taxon Caenogastropoda. Using Illumina sequencing technology, we de novo assembled the transcriptome of Littorina littorea hemocytes from 182 million mRNA-Seq pair-end 100 bp reads into a total of 15,526 contigs clustered in 4472 unigenes. The transcriptome profile was analyzed for presence of carbohydrate-binding molecules in a variety of architectural contexts. Hemocytes' repertoire of lectin-like proteins bearing conserved carbohydrate-recognition domains (CRDs) is highly diversified, including 11 of 15 lectin families earlier described in animals, as well as the novel members of lectin family found for the first time in mollusc species. The new molluscan lineage-specific domain combinations were confirmed by cloning and sequencing, including the fuco-lectin related molecules (FLReMs) composed of N-terminal region with no sequence homology to any known protein, a middle Fucolectin Tachylectin-4 Pentaxrin (FTP) domain, and a C-terminal epidermal growth factor (EGF) repeat region. The repertoire of lectin-like molecules is discussed in terms of their potential participation in the receptor phase of immune response. In total, immune-associated functions may be attributed to 70 transcripts belonging to 6 lectin families. These lectin-like genes show low overlap between species of invertebrates, suggesting relatively rapid evolution of immune-associated genes in the group. The repertoire provides valuable candidates for further characterization of the gene functions in mollusc immunity.

Fibrinogen-Related Proteins (FREPs) in Mollusks

Anti-parasite responses of the snail Biomphalaria glabrata involve antigen-reactive plasma lectins termed fibrinogen-related proteins (FREPs) comprising a C-terminal fibrinogen (FBG) domain and one or two upstream immunoglobulin domains. FREPs are highly polymorphic; they derive from several gene families with multiple loci and alleles that are diversified by exon loss, alternative splicing, and random somatic mutation (gene conversion and point mutations). Individual B. glabrata snails have dynamically distinct FREP sequence repertoires. The immune relevance of B. glabrata FREPs is indicated by FREP binding to polymorphic antigens of (snail-specific) digenean parasites and altered resistance of B. glabrata to digeneans following RNAi knockdown of FREPs. The compatibility polymorphism hypothesis proposes that FREP mutation increases the range of germline-encoded immune recognition in B. glabrata to counter antigenically-varied parasites. Somatic mutation may result from sequence exchange among tandemly arranged FREP genes in the genome, and analysis of sequence variants also suggests involvement of cytidine deaminase-like activity or epigenetic regulation. Without current indications of selection or retention of effective sequence variants toward immunological memory, FREP diversification is thought to afford B. glabrata immunity that is anticipatory but not adaptive. More remains to be learned about this system; other mollusks elaborate diversified lectins consisting of single FBG domains, and bona fide FREPs were reported from additional gastropod species, but these may not be diversified. Future comparative immunological studies and gene discovery driven by next-generation sequencing will further clarify taxonomic distribution of FREP diversification and the underlying mutator mechanisms as a component of immune function in mollusks.

A family of variable immunoglobulin and lectin domain containing molecules in the snail Biomphalaria glabrata

Technical limitations have hindered comprehensive studies of highly variable immune response molecules that are thought to have evolved due to pathogen-mediated selection such as fibrinogen-related proteins (FREPs) from Biomphalaria glabrata. FREPs combine upstream immunoglobulin superfamily (IgSF) domains with a C-terminal fibrinogen-related domain (FreD) and participate in reactions against trematode parasites. From RNAseq data we assembled a de novo reference transcriptome of B. glabrata to investigate the diversity of FREP transcripts. This study increased over two fold the number of bonafide FREP subfamilies and revealed important sequence diversity within FREP12 subfamily. We also report the discovery of related molecules that feature one or two IgSF domains associated with different C-terminal lectin domains, named C-type lectin-related proteins (CREPs) and Galectin-related protein (GREP). Together, the highly similar FREPs, CREPs and GREP were designated VIgL (Variable Immunoglobulin and Lectin domain containing molecules).

Dual analysis of host and pathogen transcriptomes in ostreid herpesvirus 1-positive Crassostrea gigas

Ostreid herpesvirus type 1 (OsHV-1) has become a problematic infective agent for the Pacific oyster Crassostrea gigas. In particular, the OsHV-1 μVar subtype has been associated with severe mortality episodes in oyster spat and juvenile oysters in France and other regions of the world. Factors enhancing the infectivity of the virus and its interactions with susceptible and resistant bivalve hosts are still to be understood, and only few studies have explored the expression of oyster or viral genes during productive infections. In this work, we have performed a dual RNA sequencing analysis on an oyster sample with a high viral load. High sequence coverage allowed us to thoroughly explore the OsHV-1 transcriptome and identify the activated molecular pathways in C. gigas. The identification of several highly induced and defence-related oyster transcripts supports the crucial role played by the innate immune system against the virus and opportunistic microbes possibly contributing to subsequent spat mortality.

Tracing evolutionary relicts of positive selection on eight malaria-related immune genes in mammals

Plasmodium-induced malaria widely infects primates and other mammals. Multiple past studies have revealed that positive selection could be the main evolutionary force triggering the genetic diversity of anti-malaria resistance-associated genes in human or primates. However, researchers focused most of their attention on the infra-generic and intra-specific genome evolution rather than analyzing the complete evolutionary history of mammals. Here we extend previous research by testing the evolutionary link of natural selection on eight candidate genes associated with malaria resistance in mammals. Three of the eight genes were detected to be affected by recombination, including TNF-α, iNOS and DARC. Positive selection was detected in the rest five immunogenes multiple times in different ancestral lineages of extant species throughout the mammalian evolution. Signals of positive selection were exposed in four malaria-related immunogenes in primates: CCL2, IL-10, HO1 and CD36. However, selection signals of G6PD have only been detected in non-primate eutherians. Significantly higher evolutionary rates and more radical amino acid replacement were also detected in primate CD36, suggesting its functional divergence from other eutherians. Prevalent positive selection throughout the evolutionary trajectory of mammalian malaria-related genes supports the arms race evolutionary hypothesis of host genetic response of mammalian immunogenes to infectious pathogens.

Functional Impact of RNA editing and ADARs on regulation of gene expression: perspectives from deep sequencing studies

Cells regulate gene expression at multiple levels leading to a balance between robustness and complexity within their proteome. One core molecular step contributing to this important balance during metazoan gene expression is RNA editing, such as the co-transcriptional recoding of RNA transcripts catalyzed by the adenosine deaminse acting on RNA (ADAR) family of enzymes. Understanding of the adenosine-to-inosine RNA editing process has been broadened considerably by the next generation sequencing (NGS) technology, which allows for in-depth demarcation of an RNA editome at nucleotide resolution. However, critical issues remain unresolved with regard to how RNA editing cooperates with other transcript-associated events to underpin regulated gene expression. Here we review the growing body of evidence, provided by recent NGS-based studies, that links RNA editing to other mechanisms of post-transcriptional RNA processing and gene expression regulation including alternative splicing, transcript stability and localization, and the biogenesis and function of microRNAs (miRNAs). We also discuss the possibility that systematic integration of NGS data may be employed to establish the rules of an “RNA editing code”, which may give us new insights into the functional consequences of RNA editing.