The complexity of expressed kappa light chains in egg-laying mammals

Best genome sequencing strategies for annotation of complex immune gene families in wildlife

BACKGROUND

The biodiversity crisis and increasing impact of wildlife disease on animal and human health provides impetus for studying immune genes in wildlife. Despite the recent boom in genomes for wildlife species, immune genes are poorly annotated in nonmodel species owing to their high level of polymorphism and complex genomic organisation. Our research over the past decade and a half on Tasmanian devils and koalas highlights the importance of genomics and accurate immune annotations to investigate disease in wildlife. Given this, we have increasingly been asked the minimum levels of genome quality required to effectively annotate immune genes in order to study immunogenetic diversity. Here we set out to answer this question by manually annotating immune genes in 5 marsupial genomes and 1 monotreme genome to determine the impact of sequencing data type, assembly quality, and automated annotation on accurate immune annotation.

RESULTS

Genome quality is directly linked to our ability to annotate complex immune gene families, with long reads and scaffolding technologies required to reassemble immune gene clusters and elucidate evolution, organisation, and true gene content of the immune repertoire. Draft-quality genomes generated from short reads with HiC or 10× Chromium linked reads were unable to achieve this. Despite mammalian BUSCOv5 scores of up to 94.1% amongst the 6 genomes, automated annotation pipelines incorrectly annotated up to 59% of manually annotated immune genes regardless of assembly quality or method of automated annotation.

CONCLUSIONS

Our results demonstrate that long reads and scaffolding technologies, alongside manual annotation, are required to accurately study the immune gene repertoire of wildlife species.

Marsupial and monotreme serum immunoglobulin binding by proteins A, G and L and anti-kangaroo antibody

Serological studies are often conducted to examine exposure to infectious agents in wildlife populations. However, specific immunological reagents for wildlife species are seldom available and can limit the study of infectious diseases in these animals. This study examined the ability of four commercially available immunoglobulin-binding reagents to bind serum immunoglobulins from 17 species within the Marsupialia and Monotremata. Serum samples were assessed for binding, using immunoblots and ELISAs (Enzyme-linked immunosorbent assays), to three microbially-derived proteins - staphylococcal protein A, streptococcal protein G and peptostreptococcal protein L. Additionally, an anti-kangaroo antibody was included for comparison. The inter- and intra-familial binding patterns of the reagents to serum immunoglobulins varied and evolutionary distance between animal species was not an accurate predictor of the ability of reagents to bind immunoglobulins. Results from this study can be used to inform the selection of appropriate immunological reagents in future serological studies in these clades.

A comparative overview of immunoglobulin genes and the generation of their diversity in tetrapods

In the past several decades, immunoglobulin (Ig) genes have been extensively characterized in many tetrapod species. This review focuses on the expressed Ig isotypes and the diversity of Ig genes in mammals, birds, reptiles, and amphibians. With regard to heavy chains, five Ig isotypes - IgM, IgD, IgG, IgA, and IgE - have been reported in mammals. Among these isotypes, IgM, IgD, and IgA (or its analog, IgX) are also found in non-mammalian tetrapods. Birds, reptiles, and amphibians express IgY, which is considered the precursor of IgG and IgE. Some species have developed unique isotypes of Ig, such as IgO in the platypus, IgF in Xenopus, and IgY (ΔFc) in ducks and turtles. The κ and λ light chains are both utilized in tetrapods, but the usage frequencies of κ and λ chains differ greatly among species. The diversity of Ig genes depends on several factors, including the germline repertoire and recombinatorial and post-recombinatorial diversity, and different species have evolved distinct mechanisms to generate antibody diversity.

Immunome database for marsupials and monotremes

BACKGROUND

To understand the evolutionary origins of our own immune system, we need to characterise the immune system of our distant relatives, the marsupials and monotremes. The recent sequencing of the genomes of two marsupials (opossum and tammar wallaby) and a monotreme (platypus) provides an opportunity to characterise the immune gene repertoires of these model organisms. This was required as many genes involved in immunity evolve rapidly and fail to be detected by automated gene annotation pipelines.

DESCRIPTION

We have developed a database of immune genes from the tammar wallaby, red-necked wallaby, northern brown bandicoot, brush-tail possum, opossum, echidna and platypus. The resource contains 2,235 newly identified sequences and 3,197 sequences which had been described previously. This comprehensive dataset was built from a variety of sources, including EST projects and expert-curated gene predictions generated through a variety of methods including chained-BLAST and sensitive HMMER searches. To facilitate systems-based research we have grouped sequences based on broad Gene Ontology categories as well as by specific functional immune groups. Sequences can be extracted by keyword, gene name, protein domain and organism name. Users can also search the database using BLAST.

CONCLUSION

The Immunome Database for Marsupials and Monotremes (IDMM) is a comprehensive database of all known marsupial and monotreme immune genes. It provides a single point of reference for genomic and transcriptomic datasets. Data from other marsupial and monotreme species will be added to the database as it become available. This resource will be utilized by marsupial and monotreme immunologists as well as researchers interested in the evolution of mammalian immunity.

Those other mammals: the immunoglobulins and T cell receptors of marsupials and monotremes

This review summarizes analyses of marsupial and monotreme immunoglobulin and T cell receptor genetics and expression published over the past decade. Analyses of recently completed whole genome sequences from the opossum and the platypus have yielded insight into the evolution of the common antigen receptor systems, as well as discovery of novel receptors that appear to have been lost in eutherian mammals. These species are also useful for investigation of the development of the immune system in organisms notable for giving birth to highly altricial young, as well as the evolution of maternal immunity through comparison of oviparous and viviparous mammals.

Cloning and structural analysis of two highly divergent IgA isotypes, IgA1 and IgA2 from the duck billed platypus, Ornithorhynchus anatinus

To trace the emergence of modern IgA isotypes during vertebrate evolution we have studied the immunoglobulin repertoire of a model monotreme, the platypus. Two highly divergent IgA-like isotypes (IgA1 and IgA2) were identified and their primary structures were determined from full-length cDNAs. A comparative analysis of the amino acid sequences for IgA from various animal species showed that the two platypus IgA isotypes form a branch clearly separated from their eutherian (placental) counterparts. However, they still conform to the general structure of eutherian IgA, with a hinge region and three constant domains. This indicates that the deletion of the second domain and the formation of a hinge region in IgA did occur very early during mammalian evolution, more than 166 million years ago. The two IgA isotypes in platypus differ in primary structure and appear to have arisen from a very early gene duplication, possibly preceding the metatherian eutherian split. Interestingly, one of these isotypes, IgA1, appears to be expressed in only the platypus, but is present in the echidna based on Southern blot analysis. The platypus may require a more effective mucosal immunity, with two highly divergent IgA forms, than the terrestrial echidna, due to its lifestyle, where it is exposed to pathogens both on land and in the water.

Immunoglobulins, antibody repertoire and B cell development

Swine share with most placental mammals the same five antibody isotypes and same two light chain types. Loci encoding lambda, kappa and Ig heavy chains appear to be organized as they are in other mammals. Swine differ from rodents and primates, but are similar to rabbits in using a single VH family (VH3) to encode their variable heavy chain domain, but not the family used by cattle, another artiodactyl. Distinct from other hoofed mammals and rodents, Ckappa:Clambda usage resembles the 1:1 ratio seen in primates. Since IgG subclasses diversified after speciation, same name subclass homologs do not exist among swine and other mammals unless very closely related. Swine possess six putative IgG subclasses that appear to have diversified by gene duplication and exon shuffle while retaining motifs that can bind to FcgammaRs, FcRn, C1q, protein A and protein G. The epithelial chorial placenta of swine and the precosial nature of their offspring have made piglets excellent models for studies on fetal antibody repertoire development and on the postnatal role of gut colonization, maternal colostrum and neonatal infection on the development of adaptive immunity during the "critical window" of immunological development. This chapter traces the study of the humoral immune system of this species through its various eras of discovery and compiles the results in tables and figures that should be a useful reference for educators and investigators.

The oldest platypus and its bearing on divergence timing of the platypus and echidna clades

Monotremes have left a poor fossil record, and paleontology has been virtually mute during two decades of discussion about molecular clock estimates of the timing of divergence between the platypus and echidna clades. We describe evidence from high-resolution x-ray computed tomography indicating that Teinolophos, an Early Cretaceous fossil from Australia's Flat Rocks locality (121-112.5 Ma), lies within the crown clade Monotremata, as a basal platypus. Strict molecular clock estimates of the divergence between platypus and echidnas range from 17 to 80 Ma, but Teinolophos suggests that the two monotreme clades were already distinct in the Early Cretaceous, and that their divergence may predate even the oldest strict molecular estimates by at least 50%. We generated relaxed molecular clock models using three different data sets, but only one yielded a date overlapping with the age of Teinolophos. Morphology suggests that Teinolophos is a platypus in both phylogenetic and ecological aspects, and tends to contradict the popular view of rapid Cenozoic monotreme diversification. Whereas the monotreme fossil record is still sparse and open to interpretation, the new data are consistent with much slower ecological, morphological, and taxonomic diversification rates for monotremes than in their sister taxon, the therian mammals. This alternative view of a deep geological history for monotremes suggests that rate heterogeneities may have affected mammalian evolution in such a way as to defeat strict molecular clock models and to challenge even relaxed molecular clock models when applied to mammalian history at a deep temporal scale.

Analysis of the expression of immunoglobulins throughout lactation suggests two periods of immune transfer in the tammar wallaby (Macropus eugenii)

Marsupial young are born in an under-developed state without mature immune responses. Prior to the maturation of an immune system, marsupial young are heavily reliant upon immune factors secreted in the milk to defend them against potential microbial pathogens in the environment. In this study, we identified and characterized the immunoglobulin heavy chain constant regions, light chains, polymeric Ig receptor (pIgR), J chain, neonatal Fc receptor (alpha chain) (FcRn) and the chemokine CCL28 from the model marsupial species, the tammar wallaby (Macropus eugenii). Low levels of conservation were seen in motifs in C alpha and C gamma associated with receptor binding and or transcytosis, and this may have potential implications for functionality. We evaluated the expression of immunoglobulin genes in the tammar mammary gland throughout lactation and found that two periods of increased expression of immunoglobulin genes occur. These two periods coincide with the birth of the young, and with its first emergence from the pouch. This increased expression may represent a strategy for maternal immunological protection of the pouch young.

TCR gamma chain diversity in the spleen of the duckbill platypus (Ornithorhynchus anatinus)

TCR gamma (TRG) chain diversity in splenic gammadelta T cells was determined for an egg-laying mammal (or monotreme), the duckbill platypus. Three distinct V subgroups were found in the expressed TRG chains and these three subgroups are members of a clade not found so far in eutherian mammals or birds. Each subgroup contains approximately five V gene segments, and their overall divergence is much less than is found in eutherians and birds, consistent with their recent evolution from an ancestral V gene segment. The platypus TRG locus also contains three C region genes and many of the residues involved in TCR function, such as interactions with CD3, were conserved in the monotreme C regions. All non-eutherian mammals (monotremes and marsupials) lacked the second cysteine residue necessary to form the intradomain disulfide bond in the C region, a loss apparently due to independent mutations in marsupials and monotremes. Monotreme TRGC regions also had among the most variation in the length of the connecting peptide region described for any species due to repeated motifs.

Status and applications of genomic resources for the gray, short-tailed opossum, Monodelphis domestica, an American marsupial model for comparative biology

Owing to its small size, favourable reproductive characteristics, and simple husbandry, the gray, short-tailed opossum, Monodelphis domestica, has become the most widely distributed and intensively utilised laboratory-bred research marsupial in the world today. This article provides an overview of the current state and future projections of genomic resources for this species and discusses the potential impact of this growing resource base on active research areas that use M. domestica as a model system. The resources discussed include: fully arrayed, bacterial artificial chromosome (BAC) libraries; an expanding linkage map; developing full-genome BAC-contig and chromosomal fluorescence in situ hybridisation maps; public websites providing access to the M. domestica whole-genome-shotgun sequence trace database and the whole-genome sequence assembly; and a new project underway to create an expressed-sequence database and microchip expression arrays for functional genomics applications. Major research areas discussed span a variety of genetic, evolutionary, physiologic, reproductive, developmental, and behavioural topics, including: comparative immunogenetics; genomic imprinting; reproductive biology; neurobiology; photobiology and carcinogenesis; genetics of lipoprotein metabolism; developmental and behavioural endocrinology; sexual differentiation and development; embryonic and fetal development; meiotic recombination; genome evolution; molecular evolution and phylogenetics; and more.

The Platypus Is in Its Place: Nuclear Genes and Indels Confirm the Sister Group Relation of Monotremes and Therians

Morphological data supports monotremes as the sister group of Theria (extant marsupials + eutherians), but phylogenetic analyses of 12 mitochondrial protein-coding genes have strongly supported the grouping of monotremes with marsupials: the Marsupionta hypothesis. Various nuclear genes tend to support Theria, but a comprehensive study of long concatenated sequences and broad taxon sampling is lacking. We therefore determined sequences from six nuclear genes and obtained additional sequences from the databases to create two large and independent nuclear data sets. One (data set I) emphasized taxon sampling and comprised five genes, with a concatenated length of 2,793 bp, from 21 species (two monotremes, six marsupials, nine placentals, and four outgroups). The other (data set II) emphasized gene sampling and comprised eight genes and three proteins, with a concatenated length of 10,773 bp or 3,669 amino acids, from five taxa (a monotreme, a marsupial, a rodent, human, and chicken). Both data sets were analyzed by parsimony, minimum evolution, maximum likelihood, and Bayesian methods using various models and data partitions. Data set I gave bootstrap support values for Theria between 55% and 100%, while support for Marsupionta was at most 12.3%. Taking base compositional bias into account generally increased the support for Theria. Data set II exclusively supported Theria, with the highest possible values and significantly rejected Marsupionta. Independent phylogenetic evidence in support of Theria was obtained from two single amino acid deletions and one insertion, while no supporting insertions and deletions were found for Marsupionta. On the basis of our data sets, the time of divergence between Monotremata and Theria was estimated at 231-217 MYA and between Marsupialia and Eutheria at 193-186 MYA. The morphological evidence for a basal position of Monotremata, well separated from Theria, is thus fully supported by the available molecular data from nuclear genes.

High variability in complementarity-determining regions compensates for a low number of V gene families in the λ light chain locus of the platypus

Based on the analysis of a panel of variable (V) region sequences from the Australian duck-billed platypus and the Australian short beaked echidna, the monotremes were found to express a highly diversified Vlambda repertoire. High variability was observed both in sequence and in length of all three CDR regions. However, all monotreme sequences were found to form a separate branch on a distance tree, and the monotremes appear to express only two Vlambda gene families. The appearance of all Vlambda gene segments in one branch on the distance tree gives further support for the notion that deletions of entire V region clans or families, followed by successive rounds of gene duplications may be a relatively common phenomenon during vertebrate evolution. Four different constant region sequences were also identified and a preferential use of certain J segments to each constant region was observed. A more detailed picture of the locus was obtained by analysis of genomic DNA by Southern blot and PCR. The organization of the lambda locus involves multiple V and several constant region genes with one or several joining segments positioned upstream of each constant region, similar to the organization in mouse and man. An mRNA frequency analysis shows that the lambda light chain accounts for more than 90% of the light chain transcripts in the spleen. The abundance and the high variability indicate that light chain diversity at the lambda locus contributes significantly to the antigen-binding repertoire in monotremes. A high lambda to kappa light chain ratio also indicates that variability in the CDR regions is more important for the repertoire size than the total number of V gene families.

Unusually similar patterns of antibody V segment diversity in distantly related marsupials

A pattern of coevolution between the V gene segments of Ig H and L chains has been noted previously by several investigators. Species with restricted germline V(H) diversity tend to have limited germline V(L) diversity, whereas species with high levels of germline V(H) diversity have more diverse V(L) gene segments. Evidence for a limited pool of V(H) but diverse V(L) gene segments in a South American opossum, Monodelphis domestica, is consistent with this marsupial being an exception to the pattern. To determine whether M. domestica is unique or the norm for marsupials, the V(H) and V(L) of an Australian possum, Trichosurus vulpecula, were characterized. The Ig repertoire in T. vulpecula is also derived from a restricted V(H) pool but a diverse V(L) pool. The V(L) gene segments of T. vulpecula are highly complex and contain lineages that predate the separation of marsupials and placental mammals. Thus, neither marsupial follows a pattern of coevolution of V(H) and V(L) gene segments observed in other mammals. Rather, marsupial V(H) and V(L) complexity appears to be evolving divergently, retaining diversity in V(L) perhaps to compensate for limited V(H) diversity. There is a high degree of similarity between the V(H) and V(L) in M. domestica and T. vulpecula, with the majority of V(L) families being shared between both species. All marsupial V(H) sequences isolated so far form a common clade of closely related sequences, and in contrast to the V(L) genes, the V(H) likely underwent a major loss of diversity early in marsupial evolution.