Identification of natural killer cell receptor clusters in the platypus genome reveals an expansion of C-type lectin genes

A Genomic and Transcriptomic Analysis of the C-Type Lectin Gene Family Reveals Highly Expanded and Diversified Repertoires in Bivalves

C-type lectins belong to a widely conserved family of lectins characterized in Metazoa. They show important functional diversity and immune implications, mainly as pathogen recognition receptors. In this work, C-type lectin-like proteins (CTLs) of a set of metazoan species were analyzed, revealing an important expansion in bivalve mollusks, which contrasted with the reduced repertoires of other mollusks, such as cephalopods. Orthology relationships demonstrated that these expanded repertoires consisted of CTL subfamilies conserved within Mollusca or Bivalvia and of lineage-specific subfamilies with orthology only between closely related species. Transcriptomic analyses revealed the importance of the bivalve subfamilies in mucosal immunity, as they were mainly expressed in the digestive gland and gills and modulated with specific stimuli. CTL domain-containing proteins that had additional domains (CTLDcps) were also studied, revealing interesting gene families with different conservation degrees of the CTL domain across orthologs from different taxa. Unique bivalve CTLDcps with specific domain architectures were revealed, corresponding to uncharacterized bivalve proteins with putative immune function according to their transcriptomic modulation, which could constitute interesting targets for functional characterization.

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.

Unraveling the LRC Evolution in Mammals: IGSF1 and A1BG Provide the Keys

Receptors of the leukocyte receptor cluster (LRC) play a range of important functions in the human immune system. However, the evolution of the LRC remains poorly understood, even in m\ammals not to mention nonmammalian vertebrates. We conducted a comprehensive bioinformatics analysis of the LRC-related genes in the publicly available genomes of six species that represent eutherian, marsupial, and monotreme lineages of mammals. As a result, the LRCs of African elephant and armadillo were characterized, two new genes, IGSF1 and A1BG, were attributed to the LRC of eutherian mammals, the LRC gene content was substantially extended in the short-tailed opossum and Tasmanian devil and, finally, four LRC genes were identified in the platypus genome. These findings have for the first time provided a solid basis for inference of the LRC phylogeny across mammals. Our analysis suggests that the mammalian LRC family likely derived from two ancestral genes, which evolved in a lineage-specific manner by expansion/contraction, extensive exon shuffling, and sequence divergence. The striking structural and functional diversity of eutherian LRC molecules appears largely lineage specific. The only family member retained in all the three mammalian lineages is a collagen-binding receptor OSCAR. Strong sequence conservation of a transmembrane domain known to associate with FcRγ suggests an adaptive role of this domain subtype in the LRC evolution.

Comparative genomics of the NKG2D ligand gene family

NKG2D ligands (NKG2DLs) are a group of stress-inducible major histocompatibility complex (MHC) class I-like molecules that act as a danger signal alerting the immune system to the presence of abnormal cells. In mammals, two families of NKG2DL genes have been identified: the MIC gene family encoded in the MHC region and the ULBP gene family encoded outside the MHC region in most species. Some mammals have a third family of NKG2DL-like class I genes which we named MILL (MHC class I-like located near the leukocyte receptor complex). Despite the fact that MILL genes are more closely related to MIC genes than ULBP genes are to MIC genes, MILL molecules do not function as NKG2DLs, and their function remains unknown. With the progress of mammalian genome projects, information on the MIC, ULBP, and MILL gene families became available in many mammalian species. Here, we summarize such information and discuss the origin and evolution of the NKG2DL gene family from the viewpoint of host-pathogen coevolution.

Identification, characterisation and expression analysis of natural killer receptor genes in Chlamydia pecorum infected koalas (Phascolarctos cinereus)

BACKGROUND

Koalas (Phascolarctos cinereus), an iconic Australian marsupial, are being heavily impacted by the spread of Chlamydia pecorum, an obligate intracellular bacterial pathogen. Koalas vary in their response to this pathogen, with some showing no symptoms, while others suffer severe symptoms leading to infertility, blindness or death. Little is known about the pathology of this disease and the immune response against it in this host. Studies have demonstrated that natural killer (NK) cells, key components of the innate immune system, are involved in the immune response to chlamydial infections in humans. These cells can directly lyse cells infected by intracellular pathogens and their ability to recognise these infected cells is mediated through NK receptors on their surface. These are encoded in two regions of the genome, the leukocyte receptor complex (LRC) and the natural killer complex (NKC). These two families evolve rapidly and different repertoires of genes, which have evolved by gene duplication, are seen in different species.

METHODS

In this study we aimed to characterise genes belonging to the NK receptor clusters in the koala by searching available koala transcriptomes using a combination of search methods. We developed a qPCR assay to quantify relative expression of four genes, two encoded within the NK receptor cluster (CLEC1B, CLEC4E) and two known to play a role in NK response to Chalmydia in humans (NCR3, PRF1).

RESULTS

We found that the NK receptor repertoire of the koala closely resembles that of the Tasmanian devil, with minimal genes in the NKC, but with lineage specific expansions in the LRC. Additional genes important for NK cell activity, NCR3 and PRF1, were also identified and characterised. In a preliminary study to investigate whether these genes are involved in the koala immune response to infection by its chlamydial pathogen, C. pecorum, we investigated the expression of four genes in koalas with active chlamydia infection, those with past infection and those without infection using qPCR. This analysis revealed that one of these four, CLEC4E, may be upregulated in response to chlamydia infection.

CONCLUSION

We have characterised genes of the NKC and LRC in koalas and have discovered evidence that one of these genes may be upregulated in koalas with chlamydia, suggesting that these receptors may play a role in the immune response of koalas to chlamydia infection.

Major histocompatibility complex genomics and human disease

Over several decades, various forms of genomic analysis of the human major histocompatibility complex (MHC) have been extremely successful in picking up many disease associations. This is to be expected, as the MHC region is one of the most gene-dense and polymorphic stretches of human DNA. It also encodes proteins critical to immunity, including several controlling antigen processing and presentation. Single-nucleotide polymorphism genotyping and human leukocyte antigen (HLA) imputation now permit the screening of large sample sets, a technique further facilitated by high-throughput sequencing. These methods promise to yield more precise contributions of MHC variants to disease. However, interpretation of MHC-disease associations in terms of the functions of variants has been problematic. Most studies confirm the paramount importance of class I and class II molecules, which are key to resistance to infection. Infection is likely driving the extreme variation of these genes across the human population, but this has been difficult to demonstrate. In contrast, many associations with autoimmune conditions have been shown to be specific to certain class I and class II alleles. Interestingly, conditions other than infections and autoimmunity are also associated with the MHC, including some cancers and neuropathies. These associations could be indirect, owing, for example, to the infectious history of a particular individual and selective pressures operating at the population level.

The hammerhead ribozyme: structure, catalysis, and gene regulation

The hammerhead ribozyme has long been considered a prototype for understanding RNA catalysis, but discrepancies between the earlier crystal structures of a minimal hammerhead self-cleaving motif and various biochemical investigations frustrated attempt to understand hammerhead ribozyme catalysis in terms of structure. With the discovery that a tertiary contact distal from the ribozyme's active site greatly enhances its catalytic prowess, and the emergence of new corresponding crystal structures of full-length hammerhead ribozymes, a unified understanding of catalysis in terms of the structure is now possible. A mechanism in which the invariant residue G12 functions as a general base, and the 2'-OH moiety of the invariant G8, itself forming a tertiary base pair with the invariant C3, is the general acid, appears consistent with both the crystal structure and biochemical experimental results. Originally discovered in the context of plant satellite RNA viruses, the hammerhead more recently has been found embedded in the 3'-untranslated region of mature mammalian mRNAs, suggesting additional biological roles in genetic regulation.

Identification of natural killer cell receptor genes in the genome of the marsupial Tasmanian devil (Sarcophilus harrisii)

Within the mammalian immune system, natural killer (NK) cells contribute to the first line of defence against infectious agents and tumours. Their activity is regulated, in part, by cell surface NK cell receptors. NK receptors can be divided into two unrelated, but functionally analogous superfamilies based on the structure of their extracellular ligand-binding domains. Receptors belonging to the C-type lectin superfamily are predominantly encoded in the natural killer complex (NKC), while receptors belonging to the immunoglobulin superfamily are predominantly encoded in the leukocyte receptor complex (LRC). Natural killer cell receptors are emerging as a rapidly evolving gene family which can display significant intra- and interspecific variation. To date, most studies have focused on eutherian mammals, with significantly less known about the evolution of these receptors in marsupials. Here, we describe the identification of 43 immunoglobulin domain-containing LRC genes in the genome of the Tasmanian devil (Sarcophilus harrisii), the largest remaining marsupial carnivore and only the second marsupial species to be studied. We also identify orthologs of NKC genes KLRK1, CD69, CLEC4E, CLEC1B, CLEC1A and an ortholog of an opossum NKC receptor. Characterisation of these regions in a second, distantly related marsupial provides new insights into the dynamic evolutionary histories of these receptors in mammals. Understanding the functional role of these genes is also important for the development of therapeutic agents against Devil Facial Tumour Disease, a contagious cancer that threatens the Tasmanian devil with extinction.

Immunogenetics of the NKG2D ligand gene family

NKG2D ligands (NKG2DLs) are a group of major histocompatibility complex (MHC) class I-like molecules, the expression of which is induced by cellular stresses such as infection, tumorigenesis, heat shock, tissue damage, and DNA damage. They act as a molecular danger signal alerting the immune system for infected or neoplastic cells. Mammals have two families of NKG2DL genes: the MHC-encoded MIC gene family and the ULBP gene family encoded outside the MHC region in most mammals. Rodents such as mice and rats lack the MIC family of ligands. Interestingly, some mammals have NKG2DL-like molecules named MILL that are phylogenetically related to MIC, but do not function as NKG2DLs. In this paper, we review our current knowledge of the MIC, ULBP, and MILL gene families in representative mammalian species and discuss the origin and evolution of the NKG2DL gene family.

The immune gene repertoire of an important viral reservoir, the Australian black flying fox

Background

Bats are the natural reservoir host for a range of emerging and re-emerging viruses, including SARS-like coronaviruses, Ebola viruses, henipaviruses and Rabies viruses. However, the mechanisms responsible for the control of viral replication in bats are not understood and there is little information available on any aspect of antiviral immunity in bats. Massively parallel sequencing of the bat transcriptome provides the opportunity for rapid gene discovery. Although the genomes of one megabat and one microbat have now been sequenced to low coverage, no transcriptomic datasets have been reported from any bat species. In this study, we describe the immune transcriptome of the Australian flying fox, Pteropus alecto, providing an important resource for identification of genes involved in a range of activities including antiviral immunity.

Results

Towards understanding the adaptations that have allowed bats to coexist with viruses, we have de novo assembled transcriptome sequence from immune tissues and stimulated cells from P. alecto. We identified about 18,600 genes involved in a broad range of activities with the most highly expressed genes involved in cell growth and maintenance, enzyme activity, cellular components and metabolism and energy pathways. 3.5% of the bat transcribed genes corresponded to immune genes and a total of about 500 immune genes were identified, providing an overview of both innate and adaptive immunity. A small proportion of transcripts found no match with annotated sequences in any of the public databases and may represent bat-specific transcripts.

Conclusions

This study represents the first reported bat transcriptome dataset and provides a survey of expressed bat genes that complement existing bat genomic data. In addition, these data provide insight into genes relevant to the antiviral responses of bats, and form a basis for examining the roles of these molecules in immune response to viral infection.

Evolution of the C-type lectin-like receptor genes of the DECTIN-1 cluster in the NK gene complex

Pattern recognition receptors are crucial in initiating and shaping innate and adaptive immune responses and often belong to families of structurally and evolutionarily related proteins. The human C-type lectin-like receptors encoded in the DECTIN-1 cluster within the NK gene complex contain prominent receptors with pattern recognition function, such as DECTIN-1 and LOX-1. All members of this cluster share significant homology and are considered to have arisen from subsequent gene duplications. Recent developments in sequencing and the availability of comprehensive sequence data comprising many species showed that the receptors of the DECTIN-1 cluster are not only homologous to each other but also highly conserved between species. Even in Caenorhabditis elegans, genes displaying homology to the mammalian C-type lectin-like receptors have been detected. In this paper, we conduct a comprehensive phylogenetic survey and give an up-to-date overview of the currently available data on the evolutionary emergence of the DECTIN-1 cluster genes.

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.

A limited role for gene duplications in the evolution of platypus venom

Gene duplication followed by adaptive selection is believed to be the primary driver of venom evolution. However, to date, no studies have evaluated the importance of gene duplications for venom evolution using a genomic approach. The availability of a sequenced genome and a venom gland transcriptome for the enigmatic platypus provides a unique opportunity to explore the role that gene duplication plays in venom evolution. Here, we identify gene duplication events and correlate them with expressed transcripts in an in-season venom gland. Gene duplicates (1,508) were identified. These duplicated pairs (421), including genes that have undergone multiple rounds of gene duplications, were expressed in the venom gland. The majority of these genes are involved in metabolism and protein synthesis not toxin functions. Twelve secretory genes including serine proteases, metalloproteinases, and protease inhibitors likely to produce symptoms of envenomation such as vasodilation and pain were detected. Only 16 of 107 platypus genes with high similarity to known toxins evolved through gene duplication. Platypus venom C-type natriuretic peptides and nerve growth factor do not possess lineage-specific gene duplicates. Extensive duplications, believed to increase the potency of toxic content and promote toxin diversification, were not found. This is the first study to take a genome-wide approach in order to examine the impact of gene duplication on venom evolution. Our findings support the idea that adaptive selection acts on gene duplicates to drive the independent evolution and functional diversification of similar venom genes in venomous species. However, gene duplications alone do not explain the "venome" of the platypus. Other mechanisms, such as alternative splicing and mutation, may be important in venom innovation.

The phylogenetic origins of natural killer receptors and recognition: relationships, possibilities, and realities

Natural killer (NK) cells affect a form of innate immunity that recognizes and eliminates cells that are infected with certain viruses or have undergone malignant transformation. In mammals, this recognition can be mediated through immunoglobulin- (Ig) and/or lectin-type NK receptors (NKRs). NKR genes in mammals range from minimally polymorphic single-copy genes to complex multigene families that exhibit high levels of haplotypic complexity and exhibit significant interspecific variation. Certain single-copy NKR genes that are present in one mammal are present as expanded multigene families in other mammals. These observations highlight NKRs as one of the most rapidly evolving eukaryotic gene families and likely reflect the influence of pathogens, especially viruses, on their evolution. Although well characterized in human and mice, cytotoxic cells that are functionally similar to NK cells have been identified in species ranging from birds to reptiles, amphibians and fish. Although numerous receptors have been identified in non-mammalian vertebrates that share structural relationships with mammalian NKRs, functionally defining these lower vertebrate molecules as NKRs is confounded by methodological and interpretive complexities. Nevertheless, several lines of evidence suggest that NK-type function or its equivalent has sustained a long evolutionary history throughout vertebrate species.

Comparative genomic analysis of mammalian NKG2D ligand family genes provides insights into their origin and evolution

NKG2D is a major activating receptor of natural killer cells. Its ligands are major histocompatibility complex (MHC) class I-like molecules whose expression is induced by cellular stresses such as infections and tumorigenesis. Humans have two families of NKG2D ligands (NKG2DL): MHC class I-related chains (MIC) encoded in the MHC and UL16-binding proteins (ULBP) encoded outside the MHC. By contrast, mice have only the latter family of ligands; instead, they have non-MHC-encoded MILL molecules that are closely related to MIC, but do not function as NKG2DL. To gain insights into the origin and evolution of MIC, ULBP, and MILL gene families, we conducted comparative genomic analysis of NKG2DL family genes in five mammalian species. In the opossum MHC, we identified a ULBP-like gene adjacent to a previously described MIC-like gene, suggesting that ULBP genes were originally encoded in the MHC. The opossum genome also contained a transcribed MILL-like gene in a region syntenic to the rodent regions encoding MILL molecules. These observations indicate that MIC-, ULBP-, and MILL-like genes emerged before the divergence of placental and marsupial mammals. Comparison of the human, cattle, rat, mouse, and opossum genomes indicates that after emigration from the MHC, ULBP genes underwent extensive duplications in each species. In mice, some of the ULBP genes appear to have been translocated telomerically on the same chromosome, forming a major cluster of existent NKG2DL genes.