Outgroup, alignment and modelling improvements indicate that two TNFSF13-like genes existed in the vertebrate ancestor

Sharks Provide Evidence for a Highly Complex TNFSF Repertoire in the Jawed Vertebrate Ancestor

Cytokines of the TNF superfamily (TNFSF) control many immunological processes and are implicated in the etiology of many immune disorders and diseases. Despite their obvious biological importance, the TNFSF repertoires of many species remain poorly characterized. In this study, we perform detailed bioinformatic, phylogenetic, and syntenic analyses of five cartilaginous fish genomes to identify their TNFSF repertoires. Strikingly, we find that shark genomes harbor ∼30 TNFSF genes, more than any other vertebrate examined to date and substantially more than humans. This is due to better retention of the ancestral jawed vertebrate TNFSF repertoire than any other jawed vertebrate lineage, combined with lineage-specific gene family expansions. All human TNFSFs appear in shark genomes, except for lymphotoxin-α (LTA; TNFSF1) and TNF (TNFSF2), and CD70 (TNFSF7) and 4-1BBL (TNFSF9), which diverged by tandem duplications early in tetrapod and mammalian evolution, respectively. Although lacking one-to-one LTA and TNF orthologs, sharks have evolved lineage-specific clusters of LTA/TNF co-orthologs. Other key findings include the presence of two BAFF (TNFSF13B) genes along with orthologs of APRIL (TNFSF13) and BALM (TNFSF13C) in sharks, and that all cartilaginous fish genomes harbor an ∼400-million-year-old cluster of multiple FASLG (TNFSF6) orthologs. Finally, sharks have retained seven ancestral jawed vertebrate TNFSF genes lost in humans. Taken together, our data indicate that the jawed vertebrate ancestor possessed a much larger and diverse TNFSF repertoire than previously hypothesized and oppose the idea that the cartilaginous fish immune system is “primitive” compared with that of mammals.

Exploration of the Nurse Shark (Ginglymostoma cirratum) Plasma Immunoproteome Using High-Resolution LC-MS/MS

Many animals of scientific importance lack species-specific reagents (e.g., monoclonal antibodies) for in-depth studies of immune proteins. Mass spectrometry (MS)-based proteomics has emerged as a useful method for monitoring changes in protein abundance and modifications in non-model species. It can be used to quantify hundreds of candidate immune molecules simultaneously without the generation of new reagents. Here, we used MS-based proteomics to identify and quantify candidate immune proteins in the plasma of the nurse shark (Ginglymostoma cirratum), a cartilaginous fish and representative of the most basal extant vertebrate lineage with an immunoglobulin-based immune system. Mass spectrometry-based LC-MS/MS was performed on the blood plasma of nurse sharks immunized with human serum albumin (n=4) or sham immunized (n=1), and sampled at days 0 (baseline control), 1, 2, 3, 5, 7, 14, 21, 28, 25, 42 and 49. An antigen-specific antibody response was experimentally confirmed post-immunization. To provide a high-quality reference to identify proteins, we assembled and annotated a multi-tissue de novo transcriptome integrating long- and short-read sequence data. This comprised 62,682 contigs containing open reading frames (ORFs) with a length >80 amino acids. Using this transcriptome, we reliably identified 626 plasma proteins which were broadly categorized into coagulation, immune, and metabolic functional groups. To assess the feasibility of performing LC-MS/MS proteomics in nurse shark in the absence of species-specific protein annotations, we compared the results to an alternative strategy, mapping peptides to proteins predicted in the genome assembly of a related species, the whale shark (Rhincodon typus). This approach reliably identified 297 proteins, indicating that useful data on the plasma proteome may be obtained in many instances despite the absence of a species-specific reference protein database. Among the plasma proteins defined against the nurse shark transcriptome, fifteen showed consistent changes in abundance across the immunized shark individuals, indicating a role in the immune response. These included alpha-2-macroglobulin (A2M) and a novel protein yet to be characterized in diverse vertebrate lineages. Overall, this study enhances genetic and protein-level resources for nurse shark research and vastly improves our understanding of the elasmobranch plasma proteome, including its remodelling following immune stimulation.

The whale shark genome reveals patterns of vertebrate gene family evolution

Chondrichthyes (cartilaginous fishes) are fundamental for understanding vertebrate evolution, yet their genomes are understudied. We report long-read sequencing of the whale shark genome to generate the best gapless chondrichthyan genome assembly yet with higher contig contiguity than all other cartilaginous fish genomes, and studied vertebrate genomic evolution of ancestral gene families, immunity, and gigantism. We found a major increase in gene families at the origin of gnathostomes (jawed vertebrates) independent of their genome duplication. We studied vertebrate pathogen recognition receptors (PRRs), which are key in initiating innate immune defense, and found diverse patterns of gene family evolution, demonstrating that adaptive immunity in gnathostomes did not fully displace germline-encoded PRR innovation. We also discovered a new toll-like receptor (TLR29) and three NOD1 copies in the whale shark. We found chondrichthyan and giant vertebrate genomes had decreased substitution rates compared to other vertebrates, but gene family expansion rates varied among vertebrate giants, suggesting substitution and expansion rates of gene families are decoupled in vertebrate genomes. Finally, we found gene families that shifted in expansion rate in vertebrate giants were enriched for human cancer-related genes, consistent with gigantism requiring adaptations to suppress cancer.

Evidence for sponges as sister to all other animals from partitioned phylogenomics with mixture models and recoding

Resolving the relationships between the major lineages in the animal tree of life is necessary to understand the origin and evolution of key animal traits. Sponges, characterized by their simple body plan, were traditionally considered the sister group of all other animal lineages, implying a gradual increase in animal complexity from unicellularity to complex multicellularity. However, the availability of genomic data has sparked tremendous controversy as some phylogenomic studies support comb jellies taking this position, requiring secondary loss or independent origins of complex traits. Here we show that incorporating site-heterogeneous mixture models and recoding into partitioned phylogenomics alleviates systematic errors that hamper commonly-applied phylogenetic models. Testing on real datasets, we show a great improvement in model-fit that attenuates branching artefacts induced by systematic error. We reanalyse key datasets and show that partitioned phylogenomics does not support comb jellies as sister to other animals at either the supermatrix or partition-specific level.

Tumor Necrosis Factor Superfamily: Ancestral Functions and Remodeling in Early Vertebrate Evolution

The evolution of the tumor necrosis factor superfamily (TNFSF) in early vertebrates is inferred by comparing the TNFSF genes found in humans and nine fishes: three agnathans, two chondrichthyans, three actinopterygians, and the sarcopterygian Latimeria chalumnae. By combining phylogenetic and synteny analyses, the TNFSF sequences detected are classified into five clusters of genes and 24 orthology groups. A model for their evolution since the origin of vertebrates is proposed. Fifteen TNFSF genes emerged from just three progenitors due to the whole-genome duplications (WGDs) that occurred before the agnathan/gnathostome split. Later, gnathostomes not only kept most of the genes emerged in the WGDs but soon added several tandem duplicates. More recently, complex, lineage-specific patterns of duplications and losses occurred in different gnathostome lineages. In agnathan species only seven to eight TNFSF genes are detected, because this lineage soon lost six of the genes emerged in the ancestral WGDs and additional losses in both hagfishes and lampreys later occurred. The orthologs of many of these lost genes are, in mammals, ligands of death-domain-containing TNFSF receptors, indicating that the extrinsic apoptotic pathway became simplified in the agnathan lineage. From the patterns of emergence of these genes, it is deduced that both the regulation of apoptosis and the control of the NF-κB pathway that depends in modern mammals on TNFSF members emerged before the ancestral vertebrate WGDs.

Ancient BCMA-like Genes Herald B Cell Regulation in Lampreys

The TNF superfamily ligands BAFF and APRIL interact with three receptors, BAFFR, BCMA, and TACI, to play discrete and crucial roles in regulating B cell selection and homeostasis in mammals. The interactions between these ligands and receptors are both specific and redundant: BAFFR binds BAFF, whereas BCMA and TACI bind to either BAFF or APRIL. In a previous phylogenetic inquiry, we identified and characterized a BAFF-like gene in lampreys, which, with hagfish, are the only extant jawless vertebrates, both of which have B-like and T-like lymphocytes. To gain insight into lymphocyte regulation in jawless vertebrates, in this study we identified two BCMA-like genes in lampreys, BCMAL1 and BCMAL2, which were found to be preferentially expressed by B-like lymphocytes. In vitro analyses indicated that the lamprey BAFF-like protein can bind to a BCMA-like receptor Ig fusion protein and to both BCMAL1- and BCMAL2-transfected cells. Discriminating regulatory roles for the two BCMA-like molecules are suggested by their differential expression before and after activation of the B-like lymphocytes in lampreys. Our composite results imply that BAFF-based mechanisms for B cell regulation evolved before the divergence of jawed and jawless vertebrates.

Discovery of All Three Types in Cartilaginous Fishes Enables Phylogenetic Resolution of the Origins and Evolution of Interferons

Interferons orchestrate host antiviral responses in jawed vertebrates. They are categorized into three classes; IFN1 and IFN3 are the primary antiviral cytokine lineages, while IFN2 responds to a broader variety of pathogens. The evolutionary relationships within and between these three classes have proven difficult to resolve. Here, we reassess interferon evolution, considering key phylogenetic pitfalls including taxon sampling, alignment quality, model adequacy, and outgroup choice. We reveal that cartilaginous fishes, and hence the jawed vertebrate ancestor, possess(ed) orthologs of all three interferon classes. We show that IFN3 groups sister to IFN1, resolve the origins of the human IFN3 lineages, and find that intronless IFN3s emerged at least three times. IFN2 genes are highly conserved, except for IFN-γ-rel, which we confirm resulted from a teleost-specific duplication. Our analyses show that IFN1 phylogeny is highly sensitive to phylogenetic error. By accounting for this, we describe a new backbone IFN1 phylogeny that implies several IFN1 genes existed in the jawed vertebrate ancestor. One of these is represented by the intronless IFN1s of tetrapods, including mammalian-like repertoires of reptile IFN1s and a subset of amphibian IFN1s, in addition to newly-identified intron-containing shark IFN1 genes. IFN-f, previously only found in teleosts, likely represents another ancestral jawed vertebrate IFN1 family member, suggesting the current classification of fish IFN1s into two groups based on the number of cysteines may need revision. The providence of the remaining fish IFN1s and the coelacanth IFN1s proved difficult to resolve, but they may also be ancestral jawed vertebrate IFN1 lineages. Finally, a large group of amphibian-specific IFN1s falls sister to all other IFN1s and was likely also present in the jawed vertebrate ancestor. Our results verify that intronless IFN1s have evolved multiple times in amphibians and indicate that no one-to-one orthology exists between mammal and reptile IFN1s. Our data also imply that diversification of the multiple IFN1s present in the jawed vertebrate ancestor has occurred through a rapid birth-death process, consistent with functional maintenance over a 450-million-year host-pathogen arms race. In summary, this study reveals a new model of interferon evolution important to our understanding of jawed vertebrate antiviral immunity.

Phylotranscriptomics suggests the jawed vertebrate ancestor could generate diverse helper and regulatory T cell subsets

BACKGROUND

The cartilaginous fishes diverged from other jawed vertebrates ~ 450 million years ago (mya). Despite this key evolutionary position, the only high-quality cartilaginous fish genome available is for the elephant shark (Callorhinchus milii), a chimaera whose ancestors split from the elasmobranch lineage ~ 420 mya. Initial analysis of this resource led to proposals that key components of the cartilaginous fish adaptive immune system, most notably their array of T cell subsets, was primitive compared to mammals. This proposal is at odds with the robust, antigen-specific antibody responses reported in elasmobranchs following immunization. To explore this discrepancy, we generated a multi-tissue transcriptome for small-spotted catshark (Scyliorhinus canicula), a tractable elasmobranch model for functional studies. We searched this, and other newly available sequence datasets, for CD4+ T cell subset-defining genes, aiming to confirm the presence or absence of each subset in cartilaginous fishes.

RESULTS

We generated a new transcriptome based on a normalised, multi-tissue RNA pool, aiming to maximise representation of tissue-specific and lowly expressed genes. We utilized multiple transcriptomic datasets and assembly variants in phylogenetic reconstructions to unambiguously identify several T cell subset-specific molecules in cartilaginous fishes for the first time, including interleukins, interleukin receptors, and key transcription factors. Our results reveal the inability of standard phylogenetic reconstruction approaches to capture the site-specific evolutionary processes of fast-evolving immune genes but show that site-heterogeneous mixture models can adequately do so.

CONCLUSIONS

Our analyses reveal that cartilaginous fishes are capable of producing a range of CD4+ T cell subsets comparable to that of mammals. Further, that the key molecules required for the differentiation and functioning of these subsets existed in the jawed vertebrate ancestor. Additionally, we highlight the importance of considering phylogenetic diversity and, where possible, utilizing multiple datasets for individual species, to accurately infer gene presence or absence at higher taxonomic levels.

Evolutionary history of the T cell receptor complex as revealed by small-spotted catshark (Scyliorhinus canicula)

In every jawed vertebrate species studied so far, the T cell receptor (TCR) complex is composed of two different TCR chains (α/β or γ/δ) and a number of CD3 subunits responsible for transmitting signals into the T cell. In this study, we characterised all of the TCR and CD3 genes of small-spotted catshark (Scyliorhinus canicula) and analysed their expression in a broad range of tissues. While the TCR complex is highly conserved across jawed vertebrates, we identified a number of differences in catshark, most notably the presence of two copies of both TCRβ and CD3γδ, and the absence of a functionally-important proline rich region from CD3ε. We also demonstrate that TCRβ has duplicated independently multiple times in jawed vertebrate evolution, bringing additional diversity to the TCR complex. This study reveals new insights about the evolutionary history of the TCR complex and raises new avenues for future exploration.