Key Points
Lymphotoxin axis is not essential for formation of O-MALT in ectotherms and birds.
Vertebrate O-MALT structures are enriched in neuroactive ligand/receptor genes.
Mammalian PPs and LNs are enriched in genes involved in olfactory transduction.
Abstract
The adaptive immune system of all jawed vertebrates relies on the presence of B and T cell lymphocytes that aggregate in specific body sites to form primary and secondary lymphoid structures. Secondary lymphoid organs include organized MALT (O-MALT) such as the tonsils and Peyer patches. O-MALT became progressively organized during vertebrate evolution, and the TNF superfamily of genes has been identified as essential for the formation and maintenance of O-MALT and other secondary and tertiary lymphoid structures in mammals. Yet, the molecular drivers of O-MALT structures found in ectotherms and birds remain essentially unknown. In this study, we provide evidence that TNFSFs, such as lymphotoxins, are likely not a universal mechanism to maintain O-MALT structures in adulthood of teleost fish, sarcopterygian fish, or birds. Although a role for TNFSF2 (TNF-α) cannot be ruled out, transcriptomics suggest that maintenance of O-MALT in nonmammalian vertebrates relies on expression of diverse genes with shared biological functions in neuronal signaling. Importantly, we identify that expression of many genes with olfactory function is a unique feature of mammalian Peyer patches but not the O-MALT of birds or ectotherms. These results provide a new view of O-MALT evolution in vertebrates and indicate that different genes with shared biological functions may have driven the formation of these lymphoid structures by a process of convergent evolution.
Introduction
The development and organization of lymphoid tissues was a vital step in the evolution of the vertebrate immune system (1–3). Lymphoid organs can be classified into primary and secondary lymphoid organs (SLOs). Primary lymphoid organs are sites of lymphogenesis where lymphocyte progenitor cells differentiate into B and T lymphocytes. These sites include the bone marrow and the thymus in mammals, the thymus and bursa of Fabricius in birds and the thymus and head kidney in teleosts. SLOs are sites where immune cells interact with each other and immune responses are activated against foreign Ags. SLOs, such as the spleen, lymph nodes (LNs), and MALT, are therefore the birthplace of effective adaptive immune responses. SLO formation is genetically programmed, whereas tertiary lymphoid structures (TLS) are considered ectopic lymphoid accumulations that appear during adulthood in response to environmental stimuli (4).
MALT are immune inductive sites located at mucosal barriers that provide increased protection at areas of high pathogen encounter, allowing for efficient Ag trapping and rapid activation of the adaptive immune response. MALT includes diffuse MALT and organized MALT (O-MALT). Diffuse MALT consists of immune cells scattered throughout the epithelium acting as sentinels against invading pathogens and is present in all vertebrates from agnathans to mammals. O-MALT are composed of clusters of lymphocytes and can be found in both ectotherms and endotherms. Similar to other SLOs, the O-MALT of endotherms is segregated into B and T cell zones, it contains follicular dendritic cells (FDCs) and the germinal center (GC) reaction allows for affinity maturation of the adaptive immune response via selection of high-affinity B cell clones. In ectotherms, however, little to modest affinity maturation can be detected, and SLOs do not have well-defined B and T cell zones (except for the spleen).
The evolutionary origins of O-MALT have long been a subject of debate among evolutionary immunologists. O-MALT was thought to have emerged in anuran amphibians as primitive lymphoid aggregates (LAs), but in 2015, primitive O-MALT structures were found in the gut and nasopharyngeal tissue of African lungfish (Sarcopterygii) (5) revealing that O-MALT is an innovation that predates the emergence of tetrapods. Lungfish LAs are thought to be SLOs because they were present in similar anatomical locations in all animals examined, and tertiary lymphoid organs (TLOs) coined as inducible LAs appeared in response to infection (5). Lungfish LAs share features with previously described LA in ectotherms because no compartmentalization into B and T cell zones or GCs were identified. Importantly, lungfish LAs are mostly composed of T cells and, to a lesser extent, B cells. Additionally, no evidence for somatic hypermutation was evident in lungfish Las, and therefore, the function of these structures remains enigmatic.
Although bona fide LAs are not present in teleost fish, a unique O-MALT–like structure known as interbranchial lymphoid tissue (ILT) has been reported in the gill arch of salmonids. The presence of ILT in Atlantic salmon is intriguing because no other lymphoid structures appear to be present in association with other teleost mucosal barriers such as the gastrointestinal tract. Similar to lungfish LAs and in contrast to mammalian SLOs, salmonid ILT mostly consists of diverse T cell clusters, shows high expression of CCL19, no expression of RAG-1, no B and T cell areas, and no GC formation (5–10). This is not surprising because GC reactions as defined in endotherms do not appear to occur in teleosts. Although salmon ILT is not present in yolk-sac larvae and first appears in juveniles (10), it is thought to be an SLO (10) and not a TLO because its localization and structure is similar in all individuals. In support, during infection, ILT decreases in size (8) rather than displaying the classical emergence of lymphocyte clusters at sites of inflammation associated with TLOs.
O-MALT evolution is also complex within endotherms. Birds, such as turkeys and chickens, do not possess LNs but do possess a pharyngeal tonsil, cecal tonsils (CT), and other O-MALT structures in their gastrointestinal tract that have a high level of organization, with FDCs that form part of the GCs (11–13). In birds, CT, similar to Peyer patches (PPs), develop during embryogenesis and therefore are present at birth (14). Despite the apparent similarity with mammalian SLOs, birds continue to surprise evolutionary immunologists because of the extensive reductionism of the immune gene repertoire in their genomes (12).
TNF superfamily (TNFSF) members play broad biological roles in cell proliferation, differentiation, inflammation, regulation of affinity maturation, and cell death (15–17). Additionally, TNFSF members are known for their importance in lymphoid tissue organogenesis and maintenance. Using several mouse knockout models (18–21), eight TNFSF members appear to be required for O-MALT formation, organization, maintenance, and function in mammals (15). Interestingly, TNFSF2 is found in all gnathostomes and some invertebrates (22), and lymphotoxins arose as a result of a gene the duplication process, as evidenced by the tandem arrangement in the human and mouse genomes and their high amino acid identity (23). Whereas knockout models of lymphotoxins and lymphotoxin receptors unequivocally indicate absence of organized SLOs (18, 24–26), TNFSF2/TNFRSF1A/1B knockout mice show very diverse phenotypes with respect to PP organogenesis (20, 21, 26–29).
In support to the TNFSF role in lymphoid development, several studies have suggested a progressive expansion of TNFSFs during vertebrate evolution, potentially explaining the progressive organization of lymphoid structures found from bony fish to mammals (5, 30). However, there are also lines of evidence against the TNFSF hypothesis. First, amphibians express TNFSF1 (LTA) and TNFSF3 (LTB); yet, their SLOs do not have GCs, and only modest levels of affinity maturation are observed (2, 31). Second, birds lack some TNFSFs that are vital for lymphocytic organization in mammals, such as TNFSF1 and TNFSF3 (13, 32) as well as other TNFSFs, such as TNFSF13 (APRIL) and TNFSF12 (TWEAK) (17). Thus, it appears that the function of lymphotoxins in amphibians is not the same as in mammals and that birds must use other mechanisms to generate and maintain O-MALT structures. Combined, these lines of evidence suggest that the molecular drivers of lymphocyte organization may be different in different vertebrate groups.
We hypothesize, in this study, that given the diversity of O-MALT structures and their different degrees of lymphocytic organization observed in vertebrates, these structures arose by convergent evolution creating lymphocytic aggregations of similar form but potentially different functions. As a consequence, whereas TNFSF may be essential for mammalian O-MALT formation and maintenance, alternative molecular drivers may be responsible for these processes in other vertebrate groups. To test this hypothesis, we use a comparative phylogenetic approach by performing RNA sequencing (RNA-Seq) of different O-MALT structures obtained from mammals, birds, sarcopterygian fish (lungfish), and bony fish, as well as in-depth analyses of TNFSF expression in the O-MALT in these four vertebrate groups. Our results support the notion that the TNFSF hypothesis likely does not explain the diversity of O-MALT structures present in vertebrates and provide a new model in which the molecular drivers of O-MALT formation may require molecules involved in neuronal signaling.
Materials and Methods
Animals
Juvenile Protopterus dolloi (slender lungfish) (0.5 kg) were obtained from ExoticFishShop.com (https://exoticfishshop.net/am. Each RNA-Seq library was prepared from tissues from one individual. All animal studies were reviewed and approved by the Office of Animal Care Compliance at the University of New Mexico (16-200384-MC, mouse protocol number 16-200497-HSC) and the United States Department of Agriculture, Beltsville Agriculture Research Center turkey protocol (number 17-008).
Tissue sampling
Lungfish nasal LAs were dissected as explained in (5
Histology
For light microscopy, tissue samples from trout ILT, lungfish LA, turkey CT, and mouse PPs were fixed in 4% paraformaldehyde overnight, transferred to 70% ethanol, and embedded in paraffin (n = 3, except for turkey, n = 2). Samples were sectioned at a thickness of 5 μm, dewaxed in xylene, and stained using H&E for general morphological analysis. A total of six sections per sample were observed, and images were acquired with a Nikon Eclipse Ti-S Inverted Microscope and NIS-Elements Advanced Research Software (Version 4.20.02).
RNA-Seq and assembly
33). The quality of the paired-end reads from our Illumina run and the SRAs were assessed using FastQC (34) and poor-quality reads were trimmed out using Trimmomattic set to default parameters (35). The trimmed reads were then assembled into de novo transcriptomes using Trinity (36, 37). The success of the assembly was assessed by realigning our raw fastq reads to the corresponding transcriptome using BWA (38) and samtools (39).
Data mining and TNFSF phylogenetic analysis
Published genomes for the representative species listed in Supplemental Table I were searched for TNFSF members in NCBI and Ensembl (40). To search for gene expression patterns in O-MALT from mouse, turkey, lungfish, and rainbow trout, TNFSF and TNFRSF protein sequences were downloaded from NCBI and used as queries for TBLASTN searches in our de novo–assembled transcriptomes (Supplemental Table II). A summary of transcriptome quality metrics is shown in Supplemental Table II. The resulting nucleotide sequences from these searches were used as queries for BLASTX searches in NCBI using hits with an E-value lower than or equal to 1 × 10−5. Only the top hit for each search was used unless the top hit was an uncharacterized sequence, in which case the second hit, if characterized, was used. To ensure the detection of all distant TNF homologs, profile hidden Markov modeling (HMM) was implemented. To do this, our de novo transcriptomes were translated into protein databases with the longest open reading frame for each sequence using TransDecoder-5.0.0 (37) set with a minimum length of 100 aa. Raw HMM profiles for TNF (PF00229) and TNFR (PF00020) were downloaded from Pfam (41) and used in HMMER (http://hmmer.org/) to search our translated transcriptomes. Sequence alignments for TNFSF1, TNFSF2, and TNFRSF3 for all vertebrate classes were conducted using MAFFT, a multiple sequence alignment program (42). Selected TNF sequences were used to construct phylogenetic trees for ligands and receptors, respectively. Neighbor-joining phylogenetic trees were constructed using the Poisson correction with a bootstrap value of 1000 in MEGA X, as previously explained in (43) (Supplemental Figs. 1, 2).
To compare all the genes expressed in the generated transcriptomes, we first used DIAMOND (44) with the UniProt database (45). Resulting UniProt accession numbers from each transcriptome were then compared in Venny2.1 (46), and the common and unique accession numbers were used for gene ontology and KEGG pathway analysis using the DAVID bioinformatics database (47). Scatterplots identifying significant KEGG pathways (p < 0.05) were created in R (48).
Data availability
The datasets generated and/or analyzed during the current study are available in the NCBI under https://www.ncbi.nlm.nih.gov/bioproject/PRJNA486850 (see Supplemental Table II). Lungfish TNF sequences were submitted to NCBI GenBank (African lungfish TNFSF accession numbers MK935171-MK935184, African lungfish TNFRSF accession numbers MK965520-MK965537, and South American TNFSF accession numbers MN536217-MN536233.
Results
Histological analysis of MALTs
As previously reported, histological examination of O-MALT in mice, turkey, African lungfish, and rainbow trout revealed the presence of poorly organized O-MALT in ectotherms, whereas highly organized O-MALT structures can be found in mice and turkeys (Fig. 1). Lungfish LAs have previously been reported to have a diameter between 300 and 350 μm (5). ILT was first discovered in salmon as lymphocytic accumulations at the base of each gill arch as well as a trailing edge to the distal end of the gill filament (9), and we also identified these structures in adult rainbow trout gill. Turkey CT and mouse PP have similar mean diameters of ∼300 μm and, as previously reported, both of these structures showed a high degree of histological organization with defined compartmentalization into B and T cell zones (49, 50). Lungfish LA, as previously described (5), showed no compartmentalization, and it was composed of random clustering of lymphocytes with no distinct zones. Although we did not attempt to identify B and T cell zones in trout ILT in this study, these are not present in salmon ILT (9) or any other bony fish SLOs studied to date but are present in the spleen of cartilaginous fish (51).
Phylogenetic representation of O-MALT formation in bony jawed vertebrates (right) and light microscopy images of H&E-stained sections from the different O-MALT structures used in this study.
Analysis of TNFSF and TNFRSF in vertebrate genomes
Bioinformatic analyses were performed among the major vertebrate classes to identify all TNFSF and TNFRSF genes with a focus on TNFSFs previously described as key factors in lymphoid tissue formation in mammals (Fig. 2A). Tacchi et al. (5) reported the presence of 18 TNF ligands and 27 TNFR in humans, 13 TNF ligands and nine TNFR in teleost, 14 TNF ligands and 15 TNFR in the coelacanth, and 13 TNF ligands and 14 TNFR in African lungfish. To revisit the TNFSF theory, we expanded our searches to include newly available genomes/transcriptomes for all vertebrate classes (Supplemental Table I). BLAST searches revealed the presence of 14 TNFSF ligands and 25 TNFSF receptors in teleost genomes. BALM was found in the newly sequenced teleost genomes, along with receptors TNFRSF1B, TNFRSF4, TNFRSF6B, TNFRSF7 (CD27), TNFRSF8, TNFRSF9, TNFRSF10B, TNFRSF11A, TNFRSF11B, TNFRSF12A, TNFRSF13B, TNFRSF14, TNFRSF18, TNFRSF19, and TNFRSF25 (Tables I, II). In coelacanth, BLAST searches did not reveal the presence of any new ligands but identified 10 more TNFR: TNFRSF1B, TNFRSF4, TNFRSF7, TNFRSF10B, TNFRSF11B, TNFRSF13B, TNFRSF17, TNFRSF18, TNFRSF25, and TNFRSF27. Searching all new African lungfish transcriptomes revealed the expression of ligands TNFSF9, TNFSF11, and BALM as well as receptors TNFRSF1B, TNFRSF8, TNFRSF12A, TNFRSF13C, and EDAR. In South American lungfish, our transcriptomes showed the expression of 10 ligands and 18 receptors. Amphibian and reptile genomes showed the expression of 12 ligands and 21 receptors and 16 ligands and 24 receptors, respectively. The presence of 12 ligands and 27 receptors was detected in the bird genomes searched. As previously reported, TNFSF1, TNFSF3, and TNFRSF3 (LTBR) were absent from bird genomes (13) (Tables I, II). Finally, our searches confirmed all previously reported TNFSF molecules in human.
Expression of TNFs vital to secondary lymphoid tissue formation in vertebrates. (A) Expression of TNFs in previously published genomes and transcriptomes. (B) Expression of TNFs in O-MALT sequenced from mouse PP, turkey CT, African lungfish LA, and trout ILT. *, an amino acid percentage identity with human ≥30%; +, amino an acid percentage identity with human between 30 and 20%; −, amino acid percentage identity with human of <20%. The red boxes denote important TNFs for secondary lymphoid tissue development that were absent in searched genomes.
Because of the low sequence similarity between TNFSF molecules (52, 53), to identify distant TNFSF homologs we performed additional structural searches using HMM using the TNF homology domain (THD). HMM analyses revealed the presence of three more TNFR (TNFRSF4, TNFRSF8, and TNFRSF13C) in South American lungfish not found through BLAST searches. In African lungfish, HMM searches identified the presence of one additional ligand, TNFSF9, and three additional receptors, TNFRSF EDAR, TNFRSF8, and TNFRSF9. HMM searches in trout, turkey, and mouse transcriptomes did not result in the identification of any novel TNFSFs or TNFRSFs (Tables I, II).
When focusing on the TNFs vital for lymphocyte organization in mice, we found that four of these molecules are evolutionarily conserved from bony fish to mammals, including TNFSF5 (CD40L), TNFRSF1A, TNFRSF5 (CD40), and TNFRSF11A. TNFSF1 and TNFSF3 showed complex evolutionary histories because they are present in the coelacanth but appear to be absent in Australian lungfish. The majority of tetrapods have TNFSF1 and TNFSF3 but they have been lost in Aves (54–56). Our analysis indicates that TNFRSF3 was lost in early tetrapods because it is present in bony fish and lungfish but not in amphibians. TNFRSF3 genes are found in mammals, suggesting deletion of TNFRSF3 in the amphibian lineage (Fig. 2A).
Phylogenetic analysis showed the homology between TNFSFs and TNFRSFs in bony vertebrates. TNFSF ligands mostly grouped within their respective family clade. In particular, within their clade, TNFSF ligands 10, 11, 5, 1, 2, 8, and 12 had bootstrap values higher than 50%. Hence, for these families, the lower sequence divergence may reflect a conservation of function. As expected, in lungfish, TNFSF ligands appeared more closely related to tetrapod TNFSF ligands compared with those of bony fish. TNFSF1 and TNFSF2 ligands shared a bootstrap value of 51%, with no clear separation among their members (Supplemental Fig. 1). As seen for some TNFSF, the phylogenetic tree for TNFRSFs showed that some receptor families are well conserved with bootstrap values higher than 50%. Specifically, the members of TNFRSF19, 19L, 27, EDAR, 16, 8, 21 and 11A families all appeared well conserved, with no clear separation among teleosts and other vertebrates. The only exception was salmon TNFRSF11A, which was more closely related to salmon and coelacanth TNFRSF5 than to other TNFRSF11A molecules. Interestingly, TNFRSF11B and 6 shared a bootstrap value of 75%, and their members appeared mixed within the two clades. Of particular relevance, we observed that in the majority of the TNFRSF clades, lungfish sequences appeared more closely related to those of teleosts than tetrapod counterparts (Supplemental Fig. 2).
TNFSF2, as predicted, was found in all scanned genomes/transcriptomes (Fig. 2). As expected, the amino acid alignment of TNFSF2 molecules from representative vertebrate species showed high similarity among them as well as with the TNFSF1 family. Specifically, several conserved amino acids were found in the 10 β-strands domains (A–H) present in human, with the highest degree of conservation observed in domain C (Fig. 3). Interestingly, the amino acid sequence of domain C in birds (Gallus gallus) showed several amino acid substitutions that deviate the TNFSF2 sequence considerably from that of mammals, amphibians, and lungfish. The corresponding receptors (TNFRSF1A and TNFRSF1B) were also found in all vertebrate groups, except for amphibians, in which surprisingly, TNFRSF1A is present but TNFRSF1B is missing (Table II).
TNFSF2 amino acid alignment. Boxes indicate amino acids that make up the β strand of the typical jelly roll fold for TNF family members. Highlighted in black are conserved amino acids in the THD. Black shading denotes the most conserved amino acids and gray denotes partially conserved amino acids. The black arrows are used to show the entire conserved domain when it did not fit on the same line. Bold letters denote a distinct subunit that makes up TNFSF2.
These results provide a revised view of the complexity of TNFSF evolution in vertebrates and provide further support for the notion that SLOs emerged in each class of vertebrates in novel ways that are not always dependent on TNFSF members such as the lymphotoxin axis.
Expression of TNFSF and TNFRSF in ectotherm and endotherm O-MALT structures
Because of the lack of published lungfish genomes or transcriptomes obtained from O-MALT, we performed RNA-Seq on O-MALT tissues from mice, turkeys, African lungfish, and rainbow trout (Supplemental Table II). Searching these new transcriptomes indicated a similar expression pattern of TNFSFs as previously shown in Fig. 2A and 2B. Interestingly, we observed the expression of all TNFs vital to lymphoid tissue formation in mammals in the LA of the lungfish. Because a lungfish genome is not available, orthology assignment between lungfish and human TNFSF molecules can only be predicted based on the percentage of amino acid identity (57). Percentage amino acid identity analysis for TNFSF ligand and receptors present in both African lungfish and humans revealed that 17 out of the 32 molecules shared >30% sequence identity supporting orthology (58). Of the eight critical TNFs identified in mammalian studies (15), lungfish TNFSF1, TNFSF5, TNFSF11, and TNFRSF11 have >30% sequence identity with their human counterparts (Table III). The lower sequence identity of the remaining four members raises the question of whether they are involved in lymphoid organ development and organization. To gain some insights into this question, we performed amino acid sequence alignments and motif analyses of vertebrate TNFSF1 and TNFRSF3 molecules. Amino acid alignments for TNFSF1 showed the presence of the conserved THD in all jawed vertebrates. Additionally, this alignment confirmed that teleost TNF-new (TNFN) does not contain the conserved THD for TNFSF1, and therefore, it is not an ortholog of mammalian TNFSF1 (22) (Fig. 4). Similar to tetrapods, lungfish TNFSF1 contained the conserved THD that makes up the typical jelly roll conformation of TNFs (52) and residues vital for binding TNFRSF3 were conserved in mammalian and African lungfish TNFSF1 sequences. In support of this, the percentage identity between human and African lungfish TNFSF1 was above 30% (Table III). Combined, this analysis suggests that lungfish TNFSF1 is an ortholog of mammalian TNFSF1. Although we did not perform amino acid sequence alignments for TNFSF3, the amino acid identity between lungfish and human TNFSF3 is ∼27% (Table III), a value not sufficient to ascertain that the lungfish molecule carries out similar functions to the mammalian counterpart. Phylogenetic analysis showed that all vertebrate TNFSF3 molecules form one clade, whereas TNFSF1 molecules appear to be intermixed with TNFSF2 in a single clade, indicating that they are closely related.
TNFSF1 amino acid alignment. Boxes indicate amino acids that make up the β strand of the typical jelly roll fold for TNF family members. Highlighted in black are conserved amino acids in the THD. The black arrows are used to show the entire conserved domain when it did not fit on the same line. The black arrowheads denoted important amino acids that bind to TNFRSF3. Each bold letter denotes a distinct subunit that makes up TNFSF1.
TNFR do not contain a THD, but instead, are composed of short cysteine-rich domains (CRDs), which are essential for the interaction between TNSRSF3 and the TNFSF1 and TNFSF3 heterodimer (52). Amino acid alignment showed that CRDs are conserved in all the vertebrate TNFRSF3 molecules analyzed (Fig. 5A). The intracellular domain of TNFRSF3 contains both conventional and unconventional TRAF binding domains in coelacanth and tetrapods (59). However, teleost TNFRSF3 molecules, as well as African and South American lungfish TNFRSF3, lack these TRAF binding domains. Instead, teleost and South American lungfish have a bony fish–specific conserved TRAF binding motif, whereas the African lungfish does not (Fig. 5B). The lack of a TRAF binding motif suggests that African lungfish TNFRSF3 may not signal the same way as its mammalian counterpart and therefore may have an alternative function. In support, percentage identity analysis of TNFRSF3 sequences aligned in Fig. 5A and 5B show high sequence identity (66.98%) between the human and mouse sequences, whereas the percentage identity between human and African lungfish molecules is only 20.9%, indicating lack of orthology (Fig. 5C). This result supports the phylogenetic analysis findings that showed two separate clades for TNFRSF3, a clade containing teleost, lungfishes, and amphibian molecules and a second clade containing human, mouse, and coelacanth molecules (Supplemental Fig. 2).
Comparison of TNFRSF3 amino acid sequences among jawed vertebrates. (A) Extracellular domain of TNFRSF3, with triangles indicating extracellular CRD and residues required for ligand binding in mammals. (B) Alignment of TNFRSF3 intracellular domain. First box indicates unconventional TRAF binding motif, second box indicates conventional TRAF binding motif, and third box indicates a reported bony fish–specific conserved TRAF binding motif. (C) Amino acid sequence identity of TNFRSF3 molecules aligned in (A) and (B), darker shades of gray indicate higher percentage identity. The amino acids shaded in black, dark gray, and light gray denote conserved amino acids going from most conserved to least conserved across vertebrate lineage, respectively. The “v” above certain columns indicates the conserved cysteines that denote the TNFRSF motif.
As expected, we also noted the absence of TNFSF1, TNFSF3, or TNFRSF3 in turkey CT, as well as absence of TNFSF1 and TNFSF3 expression in rainbow trout, suggesting that alternative molecular mechanisms other than the lymphotoxin axis drive mucosal lymphoid tissue formation in birds and teleost fish. TNFSF2 transcripts were detected in all the O-MALT transcriptomes generated in this study, as well as the receptors TNFRSF1A and TNFRSF1B. Therefore, in the absence of the lymphotoxin axis, TNFSF2 may be driving the clustering of lymphocytes in nonmammalian species.
Unbiased search of molecular drivers of O-MALT structure in vertebrates
To elucidate previously unidentified molecular drivers and biological pathways involved in O-MALT formation, we compared the transcriptomes generated from four different O-MALT tissues from each representative vertebrate species. We obtained a total of 280,740 transcripts in the mouse PP transcriptome, 164,478 transcripts in the turkey CT transcriptome, 162,353 transcripts in the lungfish LA transcriptome, and 286,901 transcripts in the trout ILT transcriptome. We first removed genes expressed in negative-control tissue (nonlymphoid tissues) transcriptomes from each of the four species. This resulted in a total of 6483 transcripts for the mouse PP, 27,172 transcripts for turkey CT, 6772 transcripts for the lungfish LA, and 13,214 transcripts for the trout ILT. As illustrated by the Venn diagram, we observed few transcripts shared among all four O-MALT transcriptomes (Fig. 6). Turkey CT and trout ILT had the most transcripts in common with 3096 (6.9%), whereas mouse PP and lungfish LA had the lowest percentage of genes in common, with only 263 (0.6%). Turkey CT had the highest number of unique genes, with 20,520 (46%), whereas mouse PP and lungfish LA had a similar percentage of unique genes, with ∼4000 (8.7–8.8%) each, and trout ILT had 8300 (18.6%) (Fig. 6).
Venn diagram of genes shared between different O-MALT transcriptomes generated in this study. Analysis was performed after removing genes found in control (non–O-MALT) tissues for each corresponding species.
We next performed gene ontology and KEGG pathway analysis of the respective gene lists that were unique or shared among O-MALTs (Fig. 7). Interestingly, genes unique to mouse PPs were significantly enriched in genes belonging to olfactory transduction (Table IV). These included olfactory receptor (OR) family 8, subfamily B, member 4; OR family 2, subfamily H, member 1; and cyclic nucleotide gated channel β 1 (Table IV). To confirm these findings, we searched previously published SRAs from mouse inguinal LNs. We found that mouse LNs are also enriched in genes belonging to the olfactory transduction pathway, including OR families and cyclic nucleotide gated channels. Interestingly, 14 genes that were unique to lungfish LA were also enriched in the olfactory transduction KEGG pathway (Table IV). Additionally, we found a total of 97 genes that were unique to turkey CT and 20 genes that were unique to trout ILT that were enriched in the neuroactive ligand/receptor pathway. These included glycine receptors, cholinergic receptors, purinergic receptors, thyroid-stimulating hormone receptors, among others (Table V). When looking at genes shared between lungfish LA and turkey CT we found that they shared 12 genes that were enriched in the neuroactive ligand/receptor interaction pathway. Combined, our findings highlight little overlap in the transcriptome of different vertebrate O-MALT structures and unveil unique biological processes that likely drive lymphocyte aggregation in each vertebrate group, such as olfactory-related genes in mammals and sarcopterygian fish and neuronal-derived signals in birds and ectotherms.
Significantly enriched KEGG pathways from O-MALT transcriptomic data after removing genes expressed in corresponding negative-control (non–O-MALT) transcriptomes for each species.
Discussion
The organization of lymphocytes and other immune cells in discrete lymphoid structures is one of the hallmarks of the immune system of endotherms. This organization is believed to increase cell–Ag interactions, optimize Ag presentation and T cell stimulation, and lead to efficient selection of high-affinity B cell clones during the maturation of the immune response. SLOs, such as spleen, LNs, and O-MALT, develop during embryogenesis. In contrast, TLSs develop after birth and require chronic inflammation or microbial signals to develop. TLS include the isolated lymphoid follicles in the small intestine, the inducible bronchial-associated lymphoid tissue, and the TLS commonly associated with tumors (5, 60–62).
TNFSFs are vital molecules for the process of lymphoid tissue formation and organization in mammals (63). TNFSF molecules and signaling pathways govern both the formation of SLOs, such as the spleen, LNs, and PPs, as well as FDC, GC, and TLS (61–65). The current model for mammalian LN and PP formation relies on lymphoid tissue inducer (LTi) cells and an initial neuronal-derived signal (66). It has been proposed that retinoic acid released by the vagus nerve induces expression of CXCL13 by mesenchymal cells (66). However, evidence for the neuronal contributions to the initiation of this process are limited because only one previous study introduced this concept and follow-up experimental support is currently lacking. Precursor LTi cells start to cluster and signal through TNFSF11 and TNFRSF11, which initiates the expression of TNFSF1 and TNFSF3 on the precursor LTi cells that then become mature LTi cells. Mature LTi cells also rely on the lymphotoxin axis to initiate the expression of chemokines, adhesion molecules, and cytokines that facilitate the attraction and retention of more hematopoietic cells causing the growth and maintenance of the LN (63, 67–69). Thus, TNFRSF3-mediated signaling in response to the LTα1β2 binding is the main pathway for promoting mammalian lymphoid tissue development. TNFSF1-deficient, TNFSF3-deficient, and TNFRSF3-deficient mice all have similar defects during secondary and tertiary lymphoid organogenesis (70–72), whereas TNFSF2 deficiency results in diverse phenotypes, ranging from altered SLO morphology to no changes in SLO appearance, depending on the mouse model used (20, 21, 26–29). Specific differences in LN and PP organogenesis have, nevertheless, been reported (63). For instance, mice deficient in TNFRSF1A, a receptor for LTα1β2, lack the formation of PPs but still retain LN development (21). Additionally, unlike LNs, PP formation requires a population of CD11c+ dendritic-like cells that express LTβ and the receptor tyrosine kinase RET. Interestingly, RET was previously described as a neuroregulator (73–75), and innate lymphoid cells (ILCs) ILC3 in the gut express RET and respond to glial-derived neurotrophic factors (76). Additionally, PPs but not LNs are absent or significantly reduced in CXCL13- or CXCR5-deficient animals (77) as well as sharpin-deficient animals (78).
Despite the clear involvement of TNFSF molecules in mammalian SLO and TLS biology, it is still unclear whether TNFSFs also govern O-MALT formation and maintenance in ectotherms and whether these structures serve similar immunological functions to those described in mammals. The evolution of TNFSF and TNFRSF is very complex, and it involved small-scale duplications as well as large, genome-wide duplications (53). The coevolution of TNFSF and TNFRSF families is characterized by functional convergence, as evidenced by the fact that multiple ligands share the same receptor and that different ligand-receptor interactions can result in the same biological outcome (53). Of interest, Drosophila has one TNF ligand (Eiger) and one TNFR ligand (Wengen) (79), and Drosophila TNFSF molecules are highly expressed in the nervous system (80). We previously reported the expansion in TNFSF genes in African lungfish compared with other ectotherms and proposed that this expansion may explain the presence of O-MALT in lungfish (5). Since that study, new genomes and transcriptomes have become available, allowing us to revisit the TNFSF phylogeny. Additional structural searches also allowed us to refine our approaches that were previously limited to BLAST searches. Thus, in this study we report that TNFSF does not appear to have expanded in sarcopterygian fish because lower numbers of overall TNFSF genes were found in lungfish compared with bony fish. However, expansion of TNFs within teleosts may not be surprising, given the multiple genome duplication events that have occurred over bony fish evolution (81). P. dolloi is also a tetraploid species, and therefore, similar events could have contributed to TNFSF expansions in this lungfish species. Because of the lack of functional experiments at this point, it is unclear whether the TNFSF genes identified in ectotherms encode molecules with homologous functions to their mammalian counterparts, and further work is needed to clarify these questions.
Given the importance of the lymphotoxin axis that was revealed in murine studies, we focused on lymphotoxin sequence analyses in this study. TNFSF1 and TNFSF3 appear to be absent in teleost fish despite the fact that salmonids have ILT in their gills, a structure predominantly composed of T cells that lack the canonical B and T cell areas and GC formation found in endotherms (7, 9). In accord, previous efforts to discern teleost lymphotoxins showed that TNFRSF3 is present in teleosts but TNFSF1 and TNFSF3 are not. Teleost fish have a novel TNF, TNFN, that at first was considered to be the candidate homolog for mammalian lymphotoxins, but a recent study and our amino acid alignment showed that the structure and behavior of TNFN did not relate to either TNFSF1 or TNFSF3 (Fig. 4) (30, 59). Rather, TNFSF14 was found to bind to TNFRSF3, showing that, originally, TNFRSF3 was the receptor for TNFSF14. Recently, single-cell RNA-Seq analysis of RAG−/− zebrafish revealed the expression of ILC1-, ILC2-, and ILC3-type markers in zebrafish (82). Specifically, an ILC3-like cell subset found in zebrafish shares some transcriptional similarities to the mammalian LTi cells and exclusively expresses tnfb (TNFSF2), the ortholog to human LTα and a marker of LTi cells. Interestingly, we detected a partial contig for fish type I TNFSF2 in trout ILT by Illumina sequencing, and therefore, it is possible that TNFSF2 signaling alone is sufficient to drive lymphocytic aggregation in nonmammalian species. Because bony fish have two types of TNFSF2 (named I and II) (83), further studies are therefore required, for instance, to ascertain a role for each TNFSF2 gene in teleost ILT formation.
Amino acid sequence alignment of TNFRSF3 from jawed vertebrates showed that, whereas the extracellular domains are highly conserved, the transmembrane domain as well as the intracellular domain are not. Thus, the percentage identity at the amino acid level between Protopterus sp. TNFRSF3 and human TNFRSF3 is ∼20%, a low value that suggests that these two molecules are not orthologs. This finding was further supported by the TNFRSF tree reported in this study (Supplemental Fig. 2). Functionally, our results suggest that TNFRSF3 signaling, as described in mammalian studies, may not occur in lungfish and further studies should investigate the biological function of lungfish TNFRSF3. Similarly, both amphibian and bony fish TNFRSF3 only shared a 20% amino acid identity with human TNFRSF3, indicating that ectotherm TNFRSF3 molecules may lack the signaling capabilities of their mammalian counterparts. Our findings, therefore, indicate that molecular signaling pathways other than the lymphotoxin pathway must be responsible for the formation and maintenance of teleost ILT and lungfish nasal O-MALT structures.
Although birds do not have LNs, they have O-MALT structures, such as CT and pharyngeal tonsils, that are structurally very similar to those found in mammals with FDCs and GCs (84). Mammalian studies have emphasized the importance of signaling through TNFRSF3 for lymphoid tissue formation, but the lack of the lymphotoxin axis in birds has not impeded their ability to form and maintain SLOs with a mammalian-like level of compartmentalization (54). The loss of lymphotoxins in birds is not a rare event. Birds have the unique ability to compensate for losses in numerous immune genes/loci, including the Ig kappa L chain, peroxiredoxin, and other TNFSF members (15, 85–87). Thus, in agreement with previously published studies (32, 54), we found that bird genomes do not contain TNFSF1, TNFSF3, or TNFRSF3. Furthermore, our comparative transcriptomic data using turkey CT suggests that these structures may be generated by very different pathways in birds and mammals because no common gene expression signatures were found. Because we performed RNA-Seq analysis on O-MALT from adult animals, the gene expression patterns likely captured maintenance rather than organogenesis of this structure. Future studies should evaluate gene expression profiling during embryogenesis of bird SLOs to gain a deeper insight into the genes involved in their early organogenesis.
The interactions between the nervous system and the immune system are complex and bidirectional (88, 89). Along with TNFSFs, neurons have been proposed to produce metabolites such as retinoic acid that could contribute to SLO formation. In addition, a subset of ILC3s are LTi-like cells characterized by the expression of retinoic acid–related orphan receptor γ T (RORγt), which are crucial initiators of SLO organogenesis (90). As mentioned earlier, ILC3 also expresses RET and responds to neurotrophic factors, providing additional evidence that neuronal signals may regulate SLO formation. Yet, this hypothesis has not been explored in other vertebrate groups, and experimental evidence is limited at this point. Our findings support a conserved role for neuronal-derived signals in SLO maintenance in adulthood. We found that neuroactive ligand and receptor pathways were enriched in trout, lungfish, turkey, and mouse SLO transcriptomic datasets. Among the genes identified in this pathway, we identified cholinergic receptors, dopamine receptors, γ-aminobutyric acid type A receptors, and cholinergic nicotinic receptors. Nicotinic acetylcholine receptors are present within SLOs, which are highly innervated with cholinergic fibers and have been shown to influence lymphocyte development (91). We found expression of five different γ-aminobutyric acid type A receptor α genes in turkey CT but not in other SLO transcriptomes examined. These receptors have been shown to affect a wide variety of immune processes, such as cytokine secretion, cell proliferation, phagocytic activity, and chemotaxis (92). Interestingly, we observed little overlap in the specific genes within this pathway expressed in each SLO, suggesting that neuronal functions rather than specific suites of genes guide SLO maintenance in different species. It is important to note that not all animals sampled in this study were the same age. Whereas African lungfish and trout tissues were obtained from juvenile/preadult stages, mouse and turkey samples were collected from adults. Thus, we cannot rule out that the gene expression patterns found in each SLO depend on the developmental stage as well as immunological and metabolic states at which each animal was sampled. Despite these caveats, we propose, in this study, that neuronal signals may guide the aggregation of T cells and B cells and the maintenance of SLO structures across vertebrates. Further studies are warranted to test this new paradigm.
OR genes are fast evolving genes and the largest gene family in mammalian genomes. In mice, 4% of all genes are olfactory genes, whereas in humans, olfactory genes make 1.4% of the genome (93). Although this expansion of OR genes in mammals is thought to reflect the importance of olfaction in the survival of this vertebrate group, ORs also play vital roles outside of the olfactory system and have widespread expression patterns (94, 95). For instance, OR ectopic expression is thought to be vital for cell-to-cell recognition, cell cycle, and migration in many cell types (96, 97). Additionally, odorant receptors regulate immune cell functions such as T cell migration and T cell accumulation in LNs (98) or chemotaxis of lung macrophages (99). The present study identified that both mammalian PPs, mammalian LNs, and lungfish LAs were significantly enriched in genes belonging to the olfactory transduction pathway. Because the lungfish LA we sampled was a nasal LA, we cannot rule out contamination of olfactory gene expression from connective tissue surrounding the LA, but contamination issues could not explain OR expression in mouse PPs or LNs. This observation is supported by the fact that ectopic expression of specific OR genes is high in mammalian lymphoid tissues. For instance, OR family 52, subfamily N, member 4 and OR family 56, subfamily B, member 1 expression is highest in human PPs out of 18 tissues measured (100) and OR 2A4 expression is highest in human LNs out of 33 tissues analyzed. Importantly, LTi cells express G protein–coupled receptor 183 (GPR183). This receptor is required for the formation of TLS such as cryptopatches and isolated lymphoid follicles (90). In mice that are deficient in G protein αi2 subunit, PPs regress and colitis is observed, emphasizing the importance of G proteins in the formation and retention of SLO and TLS (101). Expression of OR genes in different subsets of immune cells within SLOs and TLS may allow for fast and specific mechanisms of cell-to-cell communication. Given the complexity, size, and evolutionary rates of the OR gene family in mammals, further work is warranted to define the nature of OR-mediated interactions within mammalian lymphoid structures.
In summary, our findings provide new insights into the molecules that drive the complexity and diversity of lymphocytic aggregates found at mucosal tissues of vertebrates. Exhaustive TNFSF phylogenetic studies indicate that, as suggested by others, this superfamily is not universally used by all vertebrates to achieve O-MALT formation and maintenance. We propose a new paradigm, in which genes with functions in neuroactive ligand and receptor interactions as well as olfactory-related genes may drive O-MALT maintenance in a group-specific manner. Although we do not provide functional evidence for this theory, our transcriptomic analyses suggest that neuronal inputs other than retinoic acid production are required for vertebrate SLO maintenance. Additionally, ectopic expression of olfactory genes appears to be a unique innovation of mammalian SLOs. Cooperation between olfactory genes and TNFSF members, particularly TNFSF2, may also occur and be critical for the structure and function of O-MALT. Further studies are warranted to confirm these new theories in this study proposed.
Disclosures
The authors have no financial conflicts of interest.
Acknowledgments
We thank Dr. Martin F. Flajnik for critical review of the manuscript and Dr. Lijing Bu for help with bioinformatics analysis. Dr. Judy Cannon provided mouse tissue samples, Dr. Katina Krasnec provided turkey tissue samples, and Dr. M. Forstner provided Australian lungfish tissue samples. We thank Dr. Amemiya for kindly sharing lungfish RNA-Seq datasets.
Footnotes
This work was supported by the National Institutes of Health/National Science Foundation Award 1456940 (to I.S.) and National Institutes of Health/Center for Evolutionary and Theoretical Immunology/Centers of Biomedical Research Excellence Grant P20GM103452.
The datasets presented in Supplemental Table II have been submitted to the National Center for Biotechnology Information (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA4868500) under accession number PRJNA4868500. Lungfish TNF sequences have been submitted to the National Center for Biotechnology Information’s GenBank (https://www.ncbi.nlm.nih.gov/genbank/) under accession numbers MK935171–MK935184 (African lungfish TNF super family), MK965520–MK965537 (African lungfish TNFR super family), and MN536217– MN536233 (South American TNF superfamily).
The online version of this article contains supplemental material.
Abbreviations used in this article:
- CRD
- cysteine-rich domain
- CT
- cecal tonsil
- FDC
- follicular dendritic cell
- GC
- germinal center
- HMM
- hidden Markov modeling
- ILC
- innate lymphoid cell
- ILT
- interbranchial lymphoid tissue
- LA
- lymphoid aggregate
- LN
- lymph node
- LTi
- lymphoid tissue inducer
- NCBI
- National Center for Biotechnology Information
- O-MALT
- organized MALT
- OR
- olfactory receptor
- PP
- Peyer patch
- RNA-Seq
- RNA sequencing
- SLO
- secondary lymphoid organ
- SRA
- Sequence Read Archive
- THD
- TNF homology domain
- TLO
- tertiary lymphoid organ
- TLS
- tertiary lymphoid structure
- TNFN
- TNF-new
- TNFSF
- TNF superfamily.
- Received August 30, 2019.
- Accepted February 26, 2020.
- Copyright © 2020 by The American Association of Immunologists, Inc.