Alternative haplotypes of antigen processing genes in zebrafish diverged early in vertebrate evolution

An Ancestral Major Histocompatibility Complex Organization in Cartilaginous Fish: Reconstructing MHC Origin and Evolution

Cartilaginous fish (sharks, rays, and chimeras) comprise the oldest living jawed vertebrates with a mammalian-like adaptive immune system based on immunoglobulins (Ig), T-cell receptors (TCRs), and the major histocompatibility complex (MHC). Here, we show that the cartilaginous fish "adaptive MHC" is highly regimented and compact, containing (i) a classical MHC class Ia (MHC-Ia) region containing antigen processing (antigen peptide transporters and immunoproteasome) and presenting (MHC-Ia) genes, (ii) an MHC class II (MHC-II) region (with alpha and beta genes) with linkage to beta-2-microglobulin (β2m) and bromodomain-containing 2, (iii) nonclassical MHC class Ib (MHC-Ib) regions with 450 million-year-old lineages, and (iv) a complement C4 associated with the MHC-Ia region. No MHC-Ib genes were found outside of the elasmobranch MHC. Our data suggest that both MHC-I and MHC-II genes arose after the second round of whole-genome duplication (2R) on a human chromosome (huchr) 6 precursor. Further analysis of MHC paralogous regions across early branching taxa from all jawed vertebrate lineages revealed that Ig/TCR genes likely arose on a precursor of the huchr9/12/14 MHC paralog. The β2m gene is linked to the Ig/TCR genes in some vertebrates suggesting that it was present at 1R, perhaps as the donor of C1 domain to the primordial MHC gene. In sum, extant cartilaginous fish exhibit a conserved and prototypical MHC genomic organization with features found in various vertebrates, reflecting the ancestral arrangement for the jawed vertebrates.

Ancient duplication, coevolution, and selection at the MHC class IIA and IIB genes of birds

Introduction

The Major Histocompatibility Complex (MHC) of vertebrates is a dynamically evolving multigene family primarily responsible for recognizing non-self peptide antigens and triggering a pathogen-specific adaptive immune response. In birds, the MHC was previously thought to evolve via concerted evolution with high degree of gene homogenization and the rapid loss of orthology. However, the discovery of two ancient avian MHC-IIB gene lineages (DAB1 and DAB2) originating before the radiation of extant birds indicated that despite the action of concerted evolution, orthology may be detectable for long evolutionary periods.

Methods

Here, we take advantage of rapidly accumulating digital genomic resources to search for the signal of an ancient duplication at the avian MHC-IIA genes, as well as to compare phylogenetic distribution and selection between MHC-IIA and IIB gene lineages.

Results

The analysis of MHC sequences from over 230 species representing ca. 70 bird families revealed the presence of two ancient MHC-IIA gene lineages (DAA1 and DAA2) and showed that their phylogenetic distribution matches exactly the distribution of DAB1 and DAB2 lineages, suggesting tight coevolution. The early post-duplication divergence of DAA1 and DAA2 was driven by positive selection fixing radical amino acid differences within the membrane-proximal domain and, most probably, being functionally related to the interactions between α2 and β2 chains of the MHC-II heterodimer. We detected no evidence for an overall (gene-wide) relaxation or intensification of selection at either DAA1/DAB1 or DAA2/DAB2, but codon-specific differences in selection signature were found at the peptide-binding sites between the two gene lineages, perhaps implying specialization to different pathogen regimes.

Discussion

Our results suggest that specific pairing of MHC-II α and β chains may have an adaptive significance, a conclusion that advances knowledge on the macroevolution of the avian MHC-II and opens exciting novel directions for future research.

Debaryomyces hansenii supplementation in low fish meal diets promotes growth, modulates microbiota and enhances intestinal condition in juvenile marine fish

BACKGROUND

The development of a sustainable business model with social acceptance, makes necessary to develop new strategies to guarantee the growth, health, and well-being of farmed animals. Debaryomyces hansenii is a yeast species that can be used as a probiotic in aquaculture due to its capacity to i) promote cell proliferation and differentiation, ii) have immunostimulatory effects, iii) modulate gut microbiota, and/or iv) enhance the digestive function. To provide inside into the effects of D. hansenii on juveniles of gilthead seabream (Sparus aurata) condition, we integrated the evaluation of the main key performance indicators coupled with the integrative analysis of the intestine condition, through histological and microbiota state, and its transcriptomic profiling.

RESULTS

After 70 days of a nutritional trial in which a diet with low levels of fishmeal (7%) was supplemented with 1.1% of D. hansenii (17.2 × 10⁵ CFU), an increase of ca. 12% in somatic growth was observed together with an improvement in feed conversion in fish fed a yeast-supplemented diet. In terms of intestinal condition, this probiotic modulated gut microbiota without affecting the intestine cell organization, whereas an increase in the staining intensity of mucins rich in carboxylated and weakly sulphated glycoconjugates coupled with changes in the affinity for certain lectins were noted in goblet cells. Changes in microbiota were characterized by the reduction in abundance of several groups of Proteobacteria, especially those characterized as opportunistic groups. The microarrays-based transcriptomic analysis found 232 differential expressed genes in the anterior-mid intestine of S. aurata, that were mostly related to metabolic, antioxidant, immune, and symbiotic processes.

CONCLUSIONS

Dietary administration of D. hansenii enhanced somatic growth and improved feed efficiency parameters, results that were coupled to an improvement of intestinal condition as histochemical and transcriptomic tools indicated. This probiotic yeast stimulated host-microbiota interactions without altering the intestinal cell organization nor generating dysbiosis, which demonstrated its safety as a feed additive. At the transcriptomic level, D. hansenii promoted metabolic pathways, mainly protein-related, sphingolipid, and thymidylate pathways, in addition to enhance antioxidant-related intestinal mechanisms, and to regulate sentinel immune processes, potentiating the defensive capacity meanwhile maintaining the homeostatic status of the intestine.

Immune gene variation associated with chromosome-scale differences among individual zebrafish genomes

Immune genes have evolved to maintain exceptional diversity, offering robust defense against pathogens. We performed genomic assembly to examine immune gene variation in zebrafish. Gene pathway analysis identified immune genes as significantly enriched among genes with evidence of positive selection. A large subset of genes was absent from analysis of coding sequences due to apparent lack of reads, prompting us to examine genes overlapping zero coverage regions (ZCRs), defined as 2 kb stretches without mapped reads. Immune genes were identified as highly enriched within ZCRs, including over 60% of major histocompatibility complex (MHC) genes and NOD-like receptor (NLR) genes, mediators of direct and indirect pathogen recognition. This variation was most highly concentrated throughout one arm of chromosome 4 carrying a large cluster of NLR genes, associated with large-scale structural variation covering more than half of the chromosome. Our genomic assemblies uncovered alternative haplotypes and distinct complements of immune genes among individual zebrafish, including the MHC Class II locus on chromosome 8 and the NLR gene cluster on chromosome 4. While previous studies have shown marked variation in NLR genes between vertebrate species, our study highlights extensive variation in NLR gene regions between individuals of the same species. Taken together, these findings provide evidence of immune gene variation on a scale previously unknown in other vertebrate species and raise questions about potential impact on immune function.

Characteristics and expression profiles of MHC class Ⅰ molecules in Carassius auratus

Major histocompatibility complex class Ⅰ (MHC Ⅰ) molecules play a vital role in adaptive immune systems in vertebrates by presenting antigens to effector T cells. Understanding the expression profiling of MHC Ⅰ molecules in fish is essential for improving our knowledge of the relationship between microbial infection and adaptive immunity. In this study, we conducted a comprehensive analysis of MHC Ⅰ gene characteristics in Carassius auratus, an important freshwater aquaculture fish in China that is susceptible to Cyprinid herpesvirus 2 (CyHV-2) infection. We identified approximately 20 MHC Ⅰ genes discussed, including U, Z, and L lineage genes. However, only U and Z lineage proteins were identified in the kidney of Carassius auratus using high pH reversed-phase chromatography and mass spectrometry. The L lineage proteins were either not expressed or present at an extremely low level in the kidneys of Carassius auratus. We also used targeted proteomics to analyze changes in protein MHC Ⅰ molecules abundance in healthy and CyHV-2-infected Carassius auratus. We observed that five MHC Ⅰ molecules were upregulated, and Caau-UFA was downregulated in the diseased group. This study is the first to reveal the expression of MHC Ⅰ molecules at a large scale in Cyprinids, which enhances our understanding of fish adaptive immune systems.

A highly diverse set of novel immunoglobulin-like transcript (NILT) genes in zebrafish indicates a wide range of functions with complex relationships to mammalian receptors

Multiple novel immunoglobulin-like transcripts (NILTs) have been identified from salmon, trout, and carp. NILTs typically encode activating or inhibitory transmembrane receptors with extracellular immunoglobulin (Ig) domains. Although predicted to provide immune recognition in ray-finned fish, we currently lack a definitive framework of NILT diversity, thereby limiting our predictions for their evolutionary origin and function. In order to better understand the diversity of NILTs and their possible roles in immune function, we identified five NILT loci in the Atlantic salmon (Salmo salar) genome, defined 86 NILT Ig domains within a 3-Mbp region of zebrafish (Danio rerio) chromosome 1, and described 41 NILT Ig domains as part of an alternative haplotype for this same genomic region. We then identified transcripts encoded by 43 different NILT genes which reflect an unprecedented diversity of Ig domain sequences and combinations for a family of non-recombining receptors within a single species. Zebrafish NILTs include a sole putative activating receptor but extensive inhibitory and secreted forms as well as membrane-bound forms with no known signaling motifs. These results reveal a higher level of genetic complexity, interindividual variation, and sequence diversity for NILTs than previously described, suggesting that this gene family likely plays multiple roles in host immunity.

Some thoughts about what non-mammalian jawed vertebrates are telling us about antigen processing and peptide loading of MHC molecules

The major histocompatibility complex (MHC) of mammals encodes highly polymorphic classical class I and class II molecules with crucial roles in immune responses, as well as various nonclassical molecules encoded by the MHC and elsewhere in the genome that have a variety of functions. These MHC molecules are supported by antigen processing and peptide loading pathways which are well-understood in mammals. This review considers what has been learned about the MHC, MHC molecules and the supporting pathways in non-mammalian jawed vertebrates. From the initial understanding from work with the chicken MHC, a great deal of diversity in the structure and function has been found. Are there underlying principles?

•

The genomic organisation of the MHC varies enormously across jawed vertebrates.

•

Total numbers of MHC genes vary among vertebrates, with only a few classical MHC genes.

•

Some nonclassical MHC and classical pathway genes appear earlier than others.

•

Obvious co-evolution within MHC pathways occurs in some species, but not others.

•

The promiscuity of interactions may correlate with differences in genomic organisation.

Antigen Presentation and Autophagy in Teleost Adaptive Immunity

Infectious diseases are a burden for aquaculture. Antigen processing and presentation (APP) to the immune effector cells that fight pathogens is key in the adaptive immune response. At the core of the adaptive immunity that appeared in lower vertebrates during evolution are the variable genes encoding the major histocompatibility complex (MHC). MHC class I molecules mainly present peptides processed in the cytosol by the proteasome and transported to the cell surface of all cells through secretory compartments. Professional antigen-presenting cells (pAPC) also express MHC class II molecules, which normally present peptides processed from exogenous antigens through lysosomal pathways. Autophagy is an intracellular self-degradation process that is conserved in all eukaryotes and is induced by starvation to contribute to cellular homeostasis. Self-digestion during autophagy mainly occurs by the fusion of autophagosomes, which engulf portions of cytosol and fuse with lysosomes (macroautophagy) or assisted by chaperones (chaperone-mediated autophagy, CMA) that deliver proteins to lysosomes. Thus, during self-degradation, antigens can be processed to be presented by the MHC to immune effector cells, thus, linking autophagy to APP. This review is focused on the essential components of the APP that are conserved in teleost fish and the increasing evidence related to the modulation of APP and autophagy during pathogen infection.

Coevolution between MHC class I and Antigen Processing Genes in salamanders

Proteins encoded by antigen-processing genes (APGs) provide major histocompatibility complex (MHC) class I (MHC-I) with antigenic peptides. In mammals, polymorphic multigenic MHC-I family is served by monomorphic APGs, whereas in certain nonmammalian species both MHC-I and APGs are polymorphic and coevolve within stable haplotypes. Coevolution was suggested as an ancestral gnathostome feature, presumably enabling only a single highly expressed classical MHC-I gene. In this view coevolution, while optimizing some aspects of adaptive immunity, would also limit its flexibility by preventing the expansion of classical MHC-I into a multigene family. However, some nonmammalian taxa, such as salamanders, have multiple highly expressed MHC-I genes, suggesting either that coevolution is relaxed or that it does not prevent the establishment of multigene MHC-I. To distinguish between these two alternatives, we use salamanders (30 species from 16 genera representing six families) to test, within a comparative framework, a major prediction of the coevolution hypothesis: the positive correlation between MHC-I and APG diversity. We found that MHC-I diversity explained both within-individual and species-wide diversity of two APGs, TAP1 and TAP2, supporting their coevolution with MHC-I, whereas no consistent effect was detected for the other three APGs (PSMB8, PSMB9, and TAPBP). Our results imply that although coevolution occurs in salamanders, it does not preclude the expansion of the MHC-I gene family. Contrary to the previous suggestions, nonmammalian vertebrates thus may be able to accommodate diverse selection pressures with flexibility granted by rapid expansion or contraction of the MHC-I family, while retaining the benefits of coevolution between MHC-I and TAPs.

Major histocompatibility complex genes and locus organization in the Komodo dragon (Varanus komodoensis)

We performed a meta-analysis of the newly assembled Komodo dragon (Varanus komodoensis) genome to characterize the major histocompatibility complex (MHC) of the species. The MHC gene clusters of the Komodo dragon are gene dense, complex, and contain counterparts of many genes of the human MHC. Our analysis identified 20 contigs encompassing ~ 6.9 Mbp of sequence with 223 annotated genes of which many are predicted orthologs to the genes of the human MHC. These MHC contigs range in size from 13.2 kb to 21.5 Mbp, contain an average of one gene per 30 kb, and are thought to occur on at least two chromosomes. Eight contigs, each > 100 kb, could be aligned to the human MHC based on gene content, and these represent gene clusters found in each of the recognized mammalian MHC subregions. The MHC of the Komodo dragon shares organizational features of other non-mammalian taxa. Multiple class Iα and class IIβ genes are indicated, with linkage between classical class I and immunoproteasome genes and between framework class I genes and genes associated with the mammalian class III subregion. These findings are supported in both Komodo genome assemblies and provide new insight into the MHC organization of these unique squamate reptiles.

The Laboratory Domestication of Zebrafish: From Diverse Populations to Inbred Substrains

We know from human genetic studies that practically all aspects of biology are strongly influenced by the genetic background, as reflected in the advent of “personalized medicine.” Yet, with few exceptions, this is not taken into account when using laboratory populations as animal model systems for research in these fields. Laboratory strains of zebrafish (Danio rerio) are widely used for research in vertebrate developmental biology, behavior, and physiology, for modeling diseases, and for testing pharmaceutic compounds in vivo. However, all of these strains are derived from artificial bottleneck events and therefore are likely to represent only a fraction of the genetic diversity present within the species. Here, we use restriction site-associated DNA sequencing to genetically characterize wild populations of zebrafish from India, Nepal, and Bangladesh, and to compare them to previously published data on four common laboratory strains. We measured nucleotide diversity, heterozygosity, and allele frequency spectra, and find that wild zebrafish are much more diverse than laboratory strains. Further, in wild zebrafish, there is a clear signal of GC-biased gene conversion that is missing in laboratory strains. We also find that zebrafish populations in Nepal and Bangladesh are most distinct from all other strains studied, making them an attractive subject for future studies of zebrafish population genetics and molecular ecology. Finally, isolates of the same strains kept in different laboratories show a pattern of ongoing differentiation into genetically distinct substrains. Together, our findings broaden the basis for future genetic, physiological, pharmaceutic, and evolutionary studies in Danio rerio.

Characterization of novel zebrafish MHC class I U lineage genes and their haplotype

Major histocompatibility complex (MHC) genes are essential for distinguishing between individuals in all jawed vertebrates. Although MHC class I (mhc1) genes in zebrafish comprise distinct haplotypes, not all members of the mhc1 gene family have been fully characterized. In this study, we report the identification of two novel U lineage genes isolated from the WIK strain of zebrafish. These new mhc1 genes, named una and uoa, are located in tandem on chromosome 19 with >70% homology to previously isolated U genes. Sequencing of their neighboring genes revealed that una and uoa form a unique haplotype different from the previously known U lineage haplotypes. Additionally, we determined the expression profiles of U, Z, and L genes in three different tissues. These findings collectively suggest that mhc1 U lineage genes and their haplotypes in zebrafish are more divergent than previously considered, and their expression patterns vary significantly among different tissues.

Molecular Evolution of Antigen-Processing Genes in Salamanders: Do They Coevolve with MHC Class I Genes?

Proteins encoded by antigen-processing genes (APGs) prepare antigens for presentation by the major histocompatibility complex class I (MHC I) molecules. Coevolution between APGs and MHC I genes has been proposed as the ancestral gnathostome condition. The hypothesis predicts a single highly expressed MHC I gene and tight linkage between APGs and MHC I. In addition, APGs should evolve under positive selection, a consequence of the adaptive evolution in MHC I. The presence of multiple highly expressed MHC I genes in some teleosts, birds, and urodeles appears incompatible with the coevolution hypothesis. Here, we use urodele amphibians to test two key expectations derived from the coevolution hypothesis: 1) the linkage between APGs and MHC I was studied in Lissotriton newts and 2) the evidence for adaptive evolution in APGs was assessed using 42 urodele species comprising 21 genera from seven families. We demonstrated that five APGs (PSMB8, PSMB9, TAP1, TAP2, and TAPBP) are tightly linked (<0.5 cM) to MHC I. Although all APGs showed some codons under episodic positive selection, we did not find a pervasive signal of positive selection expected under the coevolution hypothesis. Gene duplications, putative gene losses, and divergent allelic lineages detected in some APGs demonstrate considerable evolutionary dynamics of APGs in salamanders. Overall, our results indicate that if coevolution between APGs and MHC I occurred in urodeles, it would be more complex than envisaged in the original formulation of the hypothesis.

MHC class I evolution; from Northern pike to salmonids

Background

Salmonids are of major importance both as farmed and wild animals. With the changing environment comes changes in pathogenic pressures so understanding the immune system of all salmonid species is of essence. Major histocompatibility complex (MHC) genes are key players in the adaptive immune system signalling infection to responding T-cells populations. Classical MHC class I (MHCI) genes, defined by high polymorphism, broad expression patterns and peptide binding ability, have a key role in inducing immunity. In salmonids, the fourth whole genome duplication that occurred 94 million years ago has provided salmonids with duplicate MHCI regions, while Northern Pike, a basal sister clade to salmonids, represent a species which has not experienced this whole genome duplication.

Results

Comparing the gene organization and evolution of MHC class I gene sequences in Northern pike versus salmonids displays a complex picture of how many of these genes evolved. Regional salmonid Ia and Ib Z lineage gene duplicates are not orthologs to the Northern pike Z lineage sequences. Instead, salmonids have experienced unique gene duplications in both duplicate regions as well as in the Salmo and Oncorhynchus branch. Species-specific gene duplications are even more pronounced for some L lineage genes.

Conclusions

Although both Northern pike as well as salmonids have expanded their U and Z lineage genes, these gene duplications occurred separately in pike and in salmonids. However, the similarity between these duplications suggest the transposable machinery was present in a common ancestor. The salmonid MHCIa and MHCIb regions were formed during the 94 MYA since the split from pike and before the Oncorhynchus and Salmo branch separated. As seen in tetrapods, the non-classical U lineage genes are diversified duplicates of their classical counterpart. One MHCI lineage, the L lineage, experienced massive species-specific gene duplications after Oncorhynchus and Salmo split approximately 25 MYA. Based on what we currently know about L lineage genes, this large variation in number of L lineage genes also signals a large functional diversity in salmonids.

Revisiting the Teleost Thymus: Current Knowledge and Future Perspectives

Simple Summary

The thymus is the immune organ producing T lymphocytes that are essential to create immunity after encountering pathogens or vaccination. This review summarizes the thymus localization and histological studies, cell composition, and function in teleost fishes. We also describe how seasonal changes, photoperiod, water temperature fluctuations, and hormones can affect thymus development in fish species. Overall, the information helps identify future studies needed to understand thymus function in fish species and the immune system’s evolutionary origins. Since fish are exposed to pathogens, especially under aquaculture conditions, knowledge about the fish thymus and T lymphocyte can also help improve fish farming protocols, considering intrinsic and environmental conditions that can contribute to achieving the best vaccine responsiveness for disease resistance.

Abstract

The thymus in vertebrates plays a critical role in producing functionally competent T-lymphocytes. Phylogenetically, the thymus emerges early during evolution in jawed cartilaginous fish, and it is usually a bilateral organ placed subcutaneously at the dorsal commissure of the operculum. In this review, we summarize the current understanding of the thymus localization, histology studies, cell composition, and function in teleost fishes. Furthermore, we consider environmental factors that affect thymus development, such as seasonal changes, photoperiod, water temperature fluctuations and hormones. Further analysis of the thymus cell distribution and function will help us understand how key stages for developing functional T cells occur in fish, and how thymus dynamics can be modulated by external factors like photoperiod. Overall, the information presented here helps identify the knowledge gaps and future steps needed for a better understanding of the immunobiology of fish thymus.

Selective Stimulation of Duplicated Atlantic Salmon MHC Pathway Genes by Interferon-Gamma

Induction of cellular immune responses rely on Major histocompatibility complex (MHC) molecules presenting pathogenic peptides to T cells. Peptide processing, transport, loading and editing is a constitutive process in most cell types, but is accelerated upon infection. Recently, an unexpected complexity in the number of functional genes involved in MHC class I peptide cleavage, peptide transport, peptide loading and editing was found in teleosts, originating from the second and third whole genome duplication events. Salmonids have expanded upon this with functional duplicates also from a fourth unique salmonid whole genome duplication. However, little is known about how individual gene duplicates respond in the context of stimulation. Here we set out to investigate how interferon gamma (IFNg) regulates the transcription of immune genes in Atlantic salmon with particular focus on gene duplicates and MHC pathways. We identified a range of response patterns in Atlantic salmon gene duplicates, with upregulation of all duplicates for some genes, like interferon regulatory factor 1 (IRF1) and interferon induced protein 44-like (IFI44.L), but only induction of one or a few duplicates of other genes, such as TAPBP and ERAP2. A master regulator turned out to be the IRF1 and not the enhanceosome as seen in mammals. If IRF1 also collaborates with CIITA and possibly NLRC5 in regulating IFNg induction of MHCI and MHCII expression in Atlantic salmon, as in zebrafish, remains to be established. Altogether, our results show the importance of deciphering between gene duplicates, as they often respond very differently to stimulation and may have different biological functions.

Cellular diversity of the regenerating caudal fin

Zebrafish faithfully regenerate their caudal fin after amputation. During this process, both differentiated cells and resident progenitors migrate to the wound site and undergo lineage-restricted, programmed cellular state transitions to populate the new regenerate. Until now, systematic characterizations of cells comprising the new regenerate and molecular definitions of their state transitions have been lacking. We hereby characterize the dynamics of gene regulatory programs during fin regeneration by creating single-cell transcriptome maps of both preinjury and regenerating fin tissues at 1/2/4 days post-amputation. We consistently identified epithelial, mesenchymal, and hematopoietic populations across all stages. We found common and cell type-specific cell cycle programs associated with proliferation. In addition to defining the processes of epithelial replenishment and mesenchymal differentiation, we also identified molecular signatures that could better distinguish epithelial and mesenchymal subpopulations in fish. The insights for natural cell state transitions during regeneration point to new directions for studying this regeneration model.

Linking Virus Discovery to Immune Responses Visualized during Zebrafish Infections

The discovery of new viruses currently outpaces our capacity for experimental examination of infection biology. To better couple virus discovery with immunology, we genetically modified zebrafish to visually report on virus infections. After generating a strain that expresses green fluorescent protein (GFP) under an interferon-stimulated gene promoter, we repeatedly observed transgenic larvae spontaneously expressing GFP days after hatching. RNA sequencing comparisons of co-housed GFP-positive and GFP-negative zebrafish revealed a naturally occurring picornavirus that induced a canonical interferon-mediated response and hundreds of antiviral defense genes not observed following immunostimulatory treatments or experimental infections with other viruses. Among the many genes induced by picornavirus infection was a large set encoding guanosine triphosphatase (GTPase) of immunity-associated proteins (GIMAPs). The GIMAP gene family is massively expanded in fish genomes and may also play a crucial role in antiviral responses in mammals, including humans. We subsequently detected zebrafish picornavirus in publicly available sequencing data from seemingly asymptomatic zebrafish in many research institutes and found that it altered gene expression in a previous study of zebrafish development. Experiments revealed a horizontal mode of virus transmission, highlighting a system for studying the spread of picornavirus infections within and between individuals. Our study describes a naturally occurring picornavirus that elicits strong antiviral responses in zebrafish and provides new strategies for simultaneously discovering viruses and their impact on vertebrate hosts.

Transcriptome Ortholog Alignment Sequence Tools (TOAST) for phylogenomic dataset assembly

Background

Advances in next-generation sequencing technologies have reduced the cost of whole transcriptome analyses, allowing characterization of non-model species at unprecedented levels. The rapid pace of transcriptomic sequencing has driven the public accumulation of a wealth of data for phylogenomic analyses, however lack of tools aimed towards phylogeneticists to efficiently identify orthologous sequences currently hinders effective harnessing of this resource.

Results

We introduce TOAST, an open source R software package that can utilize the ortholog searches based on the software Benchmarking Universal Single-Copy Orthologs (BUSCO) to assemble multiple sequence alignments of orthologous loci from transcriptomes for any group of organisms. By streamlining search, query, and alignment, TOAST automates the generation of locus and concatenated alignments, and also presents a series of outputs from which users can not only explore missing data patterns across their alignments, but also reassemble alignments based on user-defined acceptable missing data levels for a given research question.

Conclusions

TOAST provides a comprehensive set of tools for assembly of sequence alignments of orthologs for comparative transcriptomic and phylogenomic studies. This software empowers easy assembly of public and novel sequences for any target database of candidate orthologs, and fills a critically needed niche for tools that enable quantification and testing of the impact of missing data. As open-source software, TOAST is fully customizable for integration into existing or novel custom informatic pipelines for phylogenomic inference.

Software, a detailed manual, and example data files are available through github carolinafishes.github.io

An Ancient, MHC-Linked, Nonclassical Class I Lineage in Cartilaginous Fish

Cartilaginous fishes, or chondrichthyans, are the oldest jawed vertebrates that have an adaptive immune system based on the MHC and Ig superfamily-based AgR. In this basal group of jawed vertebrates, we identified a third nonclassical MHC class I lineage (UDA), which is present in all species analyzed within the two major cartilaginous subclasses, Holocephali (chimaeras) and Elasmobranchii (sharks, skates, and rays). The deduced amino acid sequences of UDA have eight out of nine typically invariant residues that bind to the N and C termini of bound peptide found in most vertebrae classical class I (UAA); additionally, the other predicted 28 peptide-binding residues are perfectly conserved in all elasmobranch UDA sequences. UDA is distinct from UAA in its differential tissue distribution and its lower expression levels and is mono- or oligomorphic unlike the highly polymorphic UAA UDA has a low copy number in elasmobranchs but is multicopy in the holocephalan spotted ratfish (Hydrolagus colliei). Using a nurse shark (Ginglymostoma cirratum) family, we found that UDA is MHC linked but separable by recombination from the tightly linked cluster of UAA, TAP, and LMP genes, the so-called class I region found in most nonmammalian vertebrates. UDA has predicted structural features that are similar to certain nonclassical class I genes in other vertebrates, and, unlike polymorpic classical class I, we anticipate that it may bind to a conserved set of specialized peptides.

Bioinformatics analyses on the immune status of renal transplant patients, a systemic research of renal transplantation

BACKGROUND

Kidney transplantation is the most effective treatment for end-stage renal disease. Allograft rejections severely affect survivals of allograft kidneys and recipients.

METHODS

Using bioinformatics approaches, the present study was designed to investigate immune status in renal transplant recipients. Fifteen datasets from Gene Expression Omnibus (GEO) were collected and analysed. Analysis of gene enrichment and protein-protein interactions were also used.

RESULTS

There were 40 differentially expressed genes (DEGs) identified in chronic rejection group when compared with stable recipients, which were enriched in allograft rejection module. There were 135 DEGs identified in acute rejection patients, compared with stable recipients, in which most genes were enriched in allograft rejection and immune deficiency. There were 288 DEGs identified in stable recipients when compared to healthy subjects. Most genes were related to chemokine signalling pathway. In integrated comparisons, expressions of MHC molecules and immunoglobulins were increased in both acute and chronic rejection; expressions of LILRB and MAP 4 K1 were increased in acute rejection patients, but not in stable recipients. There were no overlapping DEGs in blood samples of transplant recipients.

CONCLUSION

By performing bioinformatics analysis on the immune status of kidney transplant patients, the present study reports several DEGs in the renal biopsy of transplant recipients, which are requested to be validated in clinical practice.

The round goby genome provides insights into mechanisms that may facilitate biological invasions

Background

The invasive benthic round goby (Neogobius melanostomus) is the most successful temperate invasive fish and has spread in aquatic ecosystems on both sides of the Atlantic. Invasive species constitute powerful in situ experimental systems to study fast adaptation and directional selection on short ecological timescales and present promising case studies to understand factors involved the impressive ability of some species to colonize novel environments. We seize the unique opportunity presented by the round goby invasion to study genomic substrates potentially involved in colonization success.

Results

We report a highly contiguous long-read-based genome and analyze gene families that we hypothesize to relate to the ability of these fish to deal with novel environments. The analyses provide novel insights from the large evolutionary scale to the small species-specific scale. We describe expansions in specific cytochrome P450 enzymes, a remarkably diverse innate immune system, an ancient duplication in red light vision accompanied by red skin fluorescence, evolutionary patterns of epigenetic regulators, and the presence of osmoregulatory genes that may have contributed to the round goby’s capacity to invade cold and salty waters. A recurring theme across all analyzed gene families is gene expansions.

Conclusions

The expanded innate immune system of round goby may potentially contribute to its ability to colonize novel areas. Since other gene families also feature copy number expansions in the round goby, and since other Gobiidae also feature fascinating environmental adaptations and are excellent colonizers, further long-read genome approaches across the goby family may reveal whether gene copy number expansions are more generally related to the ability to conquer new habitats in Gobiidae or in fish.

Electronic supplementary material

The online version of this article (10.1186/s12915-019-0731-8) contains supplementary material, which is available to authorized users.

Differential functional patterns of memory CD4+ and CD8+ T-cells from volunteers immunized with Ty21a typhoid vaccine observed using a recombinant Escherichia coli system expressing S. Typhi proteins

It is widely accepted that CD4⁺ and CD8⁺ T-cells play a significant role in protection against Salmonella enterica serovar Typhi (S. Typhi), the causative agent of the typhoid fever. However, the antigen specificity of these T-cells remains largely unknown. Previously, we demonstrated the feasibility of using a recombinant Escherichia coli (E. coli) expression system to uncover the antigen specificity of CD4⁺ and CD8⁺ T cells. Here, we expanded these studies to include the evaluation of 12 additional S. Typhi proteins: 4 outer membrane proteins (OmpH, OmpL, OmpR, OmpX), 3 Vi-polysaccharide biosynthesis proteins (TviA, TviB, TviE), 3 cold shock proteins (CspA, CspB, CspC), and 2 conserved hypothetical proteins (Chp 1 and Chp2), all selected based on the bioinformatic analyses of the content of putative T-cell epitopes. CD4⁺ and CD8⁺ T cells from 15 adult volunteers, obtained before and 42 days after immunization with oral live attenuated Ty21a vaccine, were assessed for their functionality (i.e., production of cytokines and cytotoxic expression markers in response to stimulation with selected antigens) as measured by flow cytometry. Although volunteers differed on their T-cell antigen specificity, we observed T-cell immune responses against all S. Typhi proteins evaluated. These responses included 9 proteins, OmpH, OmpR, TviA, TviE, CspA, CspB, CspC, Chp 1 and Chp 2, which have not been previously reported to elicit T-cell responses. Interestingly, we also observed that, regardless of the protein, the functional patterns of the memory T-cells were different between CD4⁺ and CD8⁺ T cells. In sum, these studies demonstrated the feasibility of using bioinformatic analysis and the E. coli expressing system described here to uncover novel immunogenic T-cell proteins that could serve as potential targets for the production of protein-based vaccines.

Inferring the "Primordial Immune Complex": Origins of MHC Class I and Antigen Receptors Revealed by Comparative Genomics

Comparative analyses suggest that the MHC was derived from a prevertebrate "primordial immune complex" (PIC). PIC duplicated twice in the well-studied two rounds of genome-wide duplications (2R) early in vertebrate evolution, generating four MHC paralogous regions (predominantly on human chromosomes [chr] 1, 6, 9, 19). Examining chiefly the amphibian Xenopus laevis, but also other vertebrates, we identified their MHC paralogues and mapped MHC class I, AgR, and "framework" genes. Most class I genes mapped to MHC paralogues, but a cluster of Xenopus MHC class Ib genes (xnc), which previously was mapped outside of the MHC paralogues, was surrounded by genes syntenic to mammalian CD1 genes, a region previously proposed as an MHC paralogue on human chr 1. Thus, this gene block is instead the result of a translocation that we call the translocated part of the MHC paralogous region (MHCtrans) Analyses of Xenopus class I genes, as well as MHCtrans, suggest that class I arose at 1R on the chr 6/19 ancestor. Of great interest are nonrearranging AgR-like genes mapping to three MHC paralogues; thus, PIC clearly contained several AgR precursor loci, predating MHC class I/II. However, all rearranging AgR genes were found on paralogues derived from the chr 19 precursor, suggesting that invasion of a variable (V) exon by the RAG transposon occurred after 2R. We propose models for the evolutionary history of MHC/TCR/Ig and speculate on the dichotomy between the jawless (lamprey and hagfish) and jawed vertebrate adaptive immune systems, as we found genes related to variable lymphocyte receptors also map to MHC paralogues.

Major Histocompatibility Complex (MHC) Genes and Disease Resistance in Fish

Fascinating about classical major histocompatibility complex (MHC) molecules is their polymorphism. The present study is a review and discussion of the fish MHC situation. The basic pattern of MHC variation in fish is similar to mammals, with MHC class I versus class II, and polymorphic classical versus nonpolymorphic nonclassical. However, in many or all teleost fishes, important differences with mammalian or human MHC were observed: (1) The allelic/haplotype diversification levels of classical MHC class I tend to be much higher than in mammals and involve structural positions within but also outside the peptide binding groove; (2) Teleost fish classical MHC class I and class II loci are not linked. The present article summarizes previous studies that performed quantitative trait loci (QTL) analysis for mapping differences in teleost fish disease resistance, and discusses them from MHC point of view. Overall, those QTL studies suggest the possible importance of genomic regions including classical MHC class II and nonclassical MHC class I genes, whereas similar observations were not made for the genomic regions with the highly diversified classical MHC class I alleles. It must be concluded that despite decades of knowing MHC polymorphism in jawed vertebrate species including fish, firm conclusions (as opposed to appealing hypotheses) on the reasons for MHC polymorphism cannot be made, and that the types of polymorphism observed in fish may not be explained by disease-resistance models alone.

Distribution of ancient α1 and α2 domain lineages between two classical MHC class I genes and their alleles in grass carp

Major histocompatibility complex (MHC) class I molecules play a crucial role in the immune response by binding and presenting pathogen-derived peptides to specific CD8⁺ T cells. From cDNA of 20 individuals of wild grass carp (Ctenopharyngodon idellus), we could amplify one or two alleles each of classical MHC class I genes Ctid-UAA and Ctid-UBA. In total, 27 and 22 unique alleles of Ctid-UAA and Ctid-UBA were found. The leader, α1, transmembrane and cytoplasmic regions distinguish between Ctid-UAA and Ctid-UBA, and their encoded α1 domain sequences belong to the ancient lineages α1-V and α1-II, respectively, which separated several hundred million years ago. However, Ctid-UAA and Ctid-UBA share allelic lineage variation in their α2 and α3 sequences, in a pattern suggestive of past interlocus recombination events that transferred α2+α3 fragments. The allelic Ctid-UAA and Ctid-UBA variation involves ancient variation between domain lineages α2-I and α2-II, which in the present study was dated back to before the ancestral separation of teleost fish and spotted gar (> 300 million years ago). This is the first report with compelling evidence that recombination events combining different ancient α1 and α2 domain lineages had a major impact on the allelic variation of two different classical MHC class I genes within the same species.

Origin and evolution of the specialized forms of proteasomes involved in antigen presentation

Proteasomes are a multi-subunit protease complex that produces peptides bound by major histocompatibility complex (MHC) class I molecules. Phylogenetic studies indicate that two specialized forms of proteasomes, immunoproteasomes and thymoproteasomes, and the proteasome activator PA28αβ emerged in a common ancestor of jawed vertebrates which acquired adaptive immunity based on the MHC, T cell receptors, and B cell receptors ~ 500 million years ago. Comparative genomics studies now provide strong evidence that the genes coding for the immunoproteasome subunits emerged by genome-wide duplication. On the other hand, the gene encoding the thymoproteasome subunit β5t emerged by tandem duplication from the gene coding for the β5 subunit. Strikingly, birds lack immunoproteasomes, thymoproteasomes, and the proteasome activator PA28αβ, raising an interesting question of whether they have evolved any compensatory mechanisms.

A cold-blooded view of adaptive immunity

The adaptive immune system arose 500 million years ago in ectothermic (cold-blooded) vertebrates. Classically, the adaptive immune system has been defined by the presence of lymphocytes expressing recombination-activating gene (RAG)-dependent antigen receptors and the MHC. These features are found in all jawed vertebrates, including cartilaginous and bony fish, amphibians and reptiles and are most likely also found in the oldest class of jawed vertebrates, the extinct placoderms. However, with the discovery of an adaptive immune system in jawless fish based on an entirely different set of antigen receptors - the variable lymphocyte receptors - the divergence of T and B cells, and perhaps innate-like lymphocytes, goes back to the origin of all vertebrates. This Review explores how recent developments in comparative immunology have furthered our understanding of the origins and function of the adaptive immune system.

Mechanisms of Fish Macrophage Antimicrobial Immunity

Overcrowding conditions and temperatures shifts regularly manifest in large-scale infections of farmed fish, resulting in economic losses for the global aquaculture industries. Increased understanding of the functional mechanisms of fish antimicrobial host defenses is an important step forward in prevention of pathogen-induced morbidity and mortality in aquaculture setting. Like other vertebrates, macrophage-lineage cells are integral to fish immune responses and for this reason, much of the recent fish immunology research has focused on fish macrophage biology. These studies have revealed notable similarities as well as striking differences in the molecular strategies by which fish and higher vertebrates control their respective macrophage polarization and functionality. In this review, we address the current understanding of the biological mechanisms of teleost macrophage functional heterogeneity and immunity, focusing on the key cytokine regulators that control fish macrophage development and their antimicrobial armamentarium.

Whole genome duplications have provided teleosts with many roads to peptide loaded MHC class I molecules

Background

In sharks, chickens, rats, frogs, medaka and zebrafish there is haplotypic variation in MHC class I and closely linked genes involved in antigen processing, peptide translocation and peptide loading. At least in chicken, such MHCIa haplotypes of MHCIa, TAP2 and Tapasin are shown to influence the repertoire of pathogen epitopes being presented to CD8+ T-cells with subsequent effect on cell-mediated immune responses.

Results

Examining MHCI haplotype variation in Atlantic salmon using transcriptome and genome resources we found little evidence for polymorphism in antigen processing genes closely linked to the classical MHCIa genes. Looking at other genes involved in MHCI assembly and antigen processing we found retention of functional gene duplicates originating from the second vertebrate genome duplication event providing cyprinids, salmonids, and neoteleosts with the potential of several different peptide-loading complexes. One of these gene duplications has also been retained in the tetrapod lineage with orthologs in frogs, birds and opossum.

Conclusion

We postulate that the unique salmonid whole genome duplication (SGD) is responsible for eliminating haplotypic content in the paralog MHCIa regions possibly due to frequent recombination and reorganization events at early stages after the SGD. In return, multiple rounds of whole genome duplications has provided Atlantic salmon, other teleosts and even lower vertebrates with alternative peptide loading complexes. How this affects antigen presentation remains to be established.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1138-9) contains supplementary material, which is available to authorized users.

PSMB5 plays a dual role in cancer development and immunosuppression

Tumor progression and metastasis are dependent on the intrinsic properties of tumor cells and the influence of microenvironment including the immune system. It would be important to identify target drug that can inhibit cancer cell and activate immune cells. Proteasome β subunits (PSMB) family, one component of the ubiquitin-proteasome system, has been demonstrated to play an important role in tumor cells and immune cells. Therefore, we used a bioinformatics approach to examine the potential role of PSMB family. Analysis of breast TCGA and METABRIC database revealed that high expression of PSMB5 was observed in breast cancer tissue and that high expression of PSMB5 predicted worse survival. In addition, high expression of PSMB5 was observed in M2 macrophages. Based on our bioinformatics analysis, we hypothesized that PSMB5 contained immunosuppressive and oncogenic characteristics. To study the effects of PSMB5 on the cancer cell and macrophage in vitro, we silenced PSMB5 expression with shRNA in THP-1 monocytes and MDA-MB-231 cells respectively. Knockdown of PSMB5 promoted human THP-1 monocyte differentiation into M1 macrophage. On the other hand, knockdown PSMB5 gene expression inhibited MDA-MB-231 cell growth and migration by colony formation assay and boyden chamber. Collectively, our data demonstrated that delivery of PSMB5 shRNA suppressed cell growth and activated defensive M1 macrophages in vitro. Furthermore, lentiviral delivery of PSMB5 shRNA significantly decreased tumor growth in a subcutaneous mouse model. In conclusion, our bioinformatics study and functional experiments revealed that PSMB5 served as novel cancer therapeutic targets. These results also demonstrated a novel translational approach to improve cancer immunotherapy.

Use of a novel antigen expressing system to study the Salmonella enterica serovar Typhi protein recognition by T cells

Salmonella enterica serovar Typhi (S. Typhi), the causative agent of the typhoid fever, is a pathogen of great public health importance. Typhoid vaccines have the potential to be cost-effective measures towards combating this disease, yet the antigens triggering host protective immune responses are largely unknown. Given the key role of cellular-mediated immunity in S. Typhi protection, it is crucial to identify S. Typhi proteins involved in T-cell responses. Here, cells from individuals immunized with Ty21a typhoid vaccine were collected before and after immunization and used as effectors. We also used an innovative antigen expressing system based on the infection of B-cells with recombinant Escherichia coli (E. coli) expressing one of four S. Typhi gene products (i.e., SifA, OmpC, FliC, GroEL) as targets. Using flow cytometry, we found that the pattern of response to specific S. Typhi proteins was variable. Some individuals responded to all four proteins while others responded to only one or two proteins. We next evaluated whether T-cells responding to recombinant E. coli also possess the ability to respond to purified proteins. We observed that CD4⁺ cell responses, but not CD8⁺ cell responses, to recombinant E. coli were significantly associated with the responses to purified proteins. Thus, our results demonstrate the feasibility of using an E. coli expressing system to uncover the antigen specificity of T-cells and highlight its applicability to vaccine studies. These results also emphasize the importance of selecting the stimuli appropriately when evaluating CD4⁺ and CD8⁺ cell responses.

Salmonella enterica serovar Typhi (S. Typhi) is the causative agent of the life-threatening typhoid fever that affects 11.9–20.6 million individuals annually in low-income and middle-income countries. The T-cells, CD4⁺ and CD8⁺ T cells, play a significant role in protection against S. Typhi infection. Yet, the antigens triggering host protective immune responses recognized by these cells are largely unknown. To address this shortcoming, in this study we used an E. coli expression system methodology for identifying immunogenic proteins of S. Typhi. We found that although the pattern of response to individual S. Typhi proteins was variable among the typhoid vaccinees, the E. coli expressing system uncovered the antigen specificity of T-cells, and highlight its applicability to vaccine studies.