Search

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

The resolution of the gar genome affords an opportunity to examine the diversification and functional specialization of immune effector molecules at a distant and potentially informative point in phylogenetic development. Although innate immunity is effected by a particularly large number of different families of molecules, the focus here is to provide detailed characterization of several families of innate receptors that are encoded in large multigene families, for which orthologous forms can be identified in other species of bony fish but not in other vertebrate groups as well as those for which orthologs are present in other vertebrate species. The results indicate that although teleost fish and the gar, as a holostean reference species, share gene families thought previously to be restricted to the teleost fish, the manner in which the members of the multigene families of innate immune receptors have undergone diversification is different in these two major phylogenetic radiations. It appears that both the total genome duplication and different patterns of genetic selection have influenced the derivation and stabilization of innate immune genes in a substantial manner during the course of vertebrate evolution.

Antigen processing and presentation genes found within the MHC are among the most highly polymorphic genes of vertebrate genomes, providing populations with diverse immune responses to a wide array of pathogens. Here, we describe transcriptome, exome, and whole-genome sequencing of clonal zebrafish, uncovering the most extensive diversity within the antigen processing and presentation genes of any species yet examined. Our CG2 clonal zebrafish assembly provides genomic context within a remarkably divergent haplotype of the core MHC region on chromosome 19 for six expressed genes not found in the zebrafish reference genome: mhc1uga, proteasome-β 9b (psmb9b), psmb8f, and previously unknown genes psmb13b, tap2d, and tap2e We identify ancient lineages for Psmb13 within a proteasome branch previously thought to be monomorphic and provide evidence of substantial lineage diversity within each of three major trifurcations of catalytic-type proteasome subunits in vertebrates: Psmb5/Psmb8/Psmb11, Psmb6/Psmb9/Psmb12, and Psmb7/Psmb10/Psmb13. Strikingly, nearby tap2 and MHC class I genes also retain ancient sequence lineages, indicating that alternative lineages may have been preserved throughout the entire MHC pathway since early diversification of the adaptive immune system ∼500 Mya. Furthermore, polymorphisms within the three MHC pathway steps (antigen cleavage, transport, and presentation) are each predicted to alter peptide specificity. Lastly, comparative analysis shows that antigen processing gene diversity is far more extensive than previously realized (with ancient coelacanth psmb8 lineages, shark psmb13, and tap2t and psmb10 outside the teleost MHC), implying distinct immune functions and conserved roles in shaping MHC pathway evolution throughout vertebrates.

Teleost genomes encode multiple multigene families of immunoglobulin domain-containing innate immune receptors (IIIRs) with unknown function and no clear mammalian orthologs. However, the genomic organization of IIIR gene clusters and the structure and signaling motifs of the proteins they encode are similar to those of mammalian innate immune receptor families such as the killer cell immunoglobulin-like receptors (KIRs), leukocyte immunoglobulin-like receptors (LILRs), Fc receptors, triggering receptors expressed on myeloid cells (TREMs) and CD300s. Teleost IIIRs include novel immune-type receptors (NITRs); diverse immunoglobulin domain containing proteins (DICPs); polymeric immunoglobulin receptor-like proteins (PIGRLs); novel immunoglobulin-like transcripts (NILTs) and leukocyte immune-type receptors (LITRs). The accumulation of genomic sequence data has revealed that IIIR gene clusters in zebrafish display haplotypic and gene content variation. This intraspecific genetic variation, as well as significant interspecific variation, frequently confounds the identification of definitive orthologous IIIR sequences between teleost species. Nevertheless, by defining which teleost lineages encode (and do not encode) different IIIR families, predictions can be made about the presence (or absence) of specific IIIR families in each teleost lineage. It is anticipated that further investigations into available genomic resources and the sequencing of a variety of multiple teleost genomes will identify additional IIIR families and permit the modeling of the evolutionary origins of IIIRs.

Bony fish encode multiple multi-gene families of membrane receptors that are comprised of immunoglobulin (Ig) domains and are predicted to function in innate immunity. One of these families, the diverse immunoglobulin (Ig) domain-containing protein (DICP) genes, maps to three chromosomal loci in zebrafish. Most DICPs possess one or two Ig ectodomains and include membrane-bound and secreted forms. Membrane-bound DICPs include putative inhibitory and activating receptors. Recombinant DICP Ig domains bind lipids with varying specificity, a characteristic shared with mammalian CD300 and TREM family members. Numerous DICP transcripts amplified from different lines of zebrafish did not match the zebrafish reference genome sequence suggesting polymorphic and haplotypic variation. The expression of DICPs in three different lines of zebrafish has been characterized employing PCR-based strategies. Certain DICPs exhibit restricted expression in adult tissues whereas others are expressed ubiquitously. Transcripts of a subset of DICPs can be detected during embryonic development suggesting roles in embryonic immunity or other developmental processes. Transcripts representing eleven previously uncharacterized DICP sequences were identified. The assignment of two of these sequences to an unplaced genomic scaffold resulted in the identification of an alternative DICP haplotype that is linked to a MHC class I Z lineage haplotype on zebrafish chromosome 3. The linkage of DICP and MHC Class I genes also is observable in the genomes of the related grass carp (Ctenopharyngodon idellus) and common carp (Cyprinus carpio) suggesting that this is a shared character with the last common Cyprinidae ancestor.

Nat. Genet. 48, 427–437 (2016); published online 7 March 2016; corrected after print 25 April 2016 As we intended, other researchers have been able to use the draft spotted gar genome sequence available from the Broad Institute website since December 2011, the assembly LepOcu1 publicly available from NCBI since 13 January 2012 under accession code GCA000242695.