Reconstructing an ancestral mammalian immune supercomplex from a marsupial major histocompatibility complex

The Amphibian Major Histocompatibility Complex-A Review and Future Outlook

The major histocompatibility complex (MHC) is a cluster of functionally related genes encoding proteins which, among other functions, mediate immune system activation. While the MHC of many vertebrates has been extensively studied, less is known about the amphibian MHC. This represents an important knowledge gap because amphibians mark the evolutionary transition from an aquatic to a terrestrial lifestyle and often maintain a biphasic lifestyle. Hence, they tend to be exposed to both aquatic and terrestrial pathogen communities, providing opportunities to gain fundamental insights into how the immune system responds to different environmental challenges. Moreover, amphibians are globally threatened by invasive pathogens and the MHC may play a role in combating population decline. In this review, we summarize the current state of knowledge regarding the amphibian MHC and identify the major differences with other vertebrates. We also review how the number of MHC gene copies varies across amphibian groups and how MHC-based variation relates to amphibian ontogeny, behaviour, disease, and phylogeography. We conclude by identifying knowledge gaps and proposing priorities for future research.

Plethora of New Marsupial Genomes Informs Our Knowledge of Marsupial MHC Class II

The major histocompatibility complex (MHC) plays a vital role in the vertebrate immune system due to its role in infection, disease and autoimmunity, or recognition of "self". The marsupial MHC class II genes show divergence from eutherian MHC class II genes and are a unique taxon of therian mammals that give birth to altricial and immunologically naive young providing an opportune study system for investigating evolution of the immune system. Additionally, the MHC in marsupials has been implicated in disease associations, including susceptibility to Chlamydia pecorum infection in koalas. Due to the complexity of the gene family, automated annotation is not possible so here we manually annotate 384 class II MHC genes in 29 marsupial species. We find losses of key components of the marsupial MHC repertoire in the Dasyuromorphia order and the Pseudochiridae family. We perform PGLS analysis to show the gene losses we find are true gene losses and not artifacts of unresolved genome assembly. We investigate the associations between the number of loci and life history traits, including lifespan and reproductive output in lineages of marsupials and hypothesize that gene loss may be linked to the energetic cost and tradeoffs associated with pregnancy and reproduction. We found support for litter size being a significant predictor of the number of DBA and DBB loci, indicating a tradeoff between the energetic requirements of immunity and reproduction. Additionally, we highlight the increased susceptibility of Dasyuridae species to neoplasia and a potential link to MHC gene loss. Finally, these annotations provide a valuable resource to the immunogenetics research community to move forward and further investigate diversity in MHC genes in marsupials.

Extant and extinct bilby genomes combined with Indigenous knowledge improve conservation of a unique Australian marsupial

Ninu (greater bilby, Macrotis lagotis) are desert-dwelling, culturally and ecologically important marsupials. In collaboration with Indigenous rangers and conservation managers, we generated the Ninu chromosome-level genome assembly (3.66 Gbp) and genome sequences for the extinct Yallara (lesser bilby, Macrotis leucura). We developed and tested a scat single-nucleotide polymorphism panel to inform current and future conservation actions, undertake ecological assessments and improve our understanding of Ninu genetic diversity in managed and wild populations. We also assessed the beneficial impact of translocations in the metapopulation (N = 363 Ninu). Resequenced genomes (temperate Ninu, 6; semi-arid Ninu, 6; and Yallara, 4) revealed two major population crashes during global cooling events for both species and differences in Ninu genes involved in anatomical and metabolic pathways. Despite their 45-year captive history, Ninu have fewer long runs of homozygosity than other larger mammals, which may be attributable to their boom-bust life history. Here we investigated the unique Ninu biology using 12 tissue transcriptomes revealing expression of all 115 conserved eutherian chorioallantoic placentation genes in the uterus, an XY1Y2 sex chromosome system and olfactory receptor gene expansions. Together, we demonstrate the holistic value of genomics in improving key conservation actions, understanding unique biological traits and developing tools for Indigenous rangers to monitor remote wild populations.

A pan-cetacean MHC amplicon sequencing panel developed and evaluated in combination with genome assemblies

The major histocompatibility complex (MHC) is a highly polymorphic gene family that is crucial in immunity, and its diversity can be effectively used as a fitness marker for populations. Despite this, MHC remains poorly characterised in non-model species (e.g., cetaceans: whales, dolphins and porpoises) as high gene copy number variation, especially in the fast-evolving class I region, makes analyses of genomic sequences difficult. To date, only small sections of class I and IIa genes have been used to assess functional diversity in cetacean populations. Here, we undertook a systematic characterisation of the MHC class I and IIa regions in available cetacean genomes. We extracted full-length gene sequences to design pan-cetacean primers that amplified the complete exon 2 from MHC class I and IIa genes in one combined sequencing panel. We validated this panel in 19 cetacean species and described 354 alleles for both classes. Furthermore, we identified likely assembly artefacts for many MHC class I assemblies based on the presence of class I genes in the amplicon data compared to missing genes from genomes. Finally, we investigated MHC diversity using the panel in 25 humpback and 30 southern right whales, including four paternity trios for humpback whales. This revealed copy-number variable class I haplotypes in humpback whales, which is likely a common phenomenon across cetaceans. These MHC alleles will form the basis for a cetacean branch of the Immuno-Polymorphism Database (IPD-MHC), a curated resource intended to aid in the systematic compilation of MHC alleles across several species, to support conservation initiatives.

Expanding the evo-devo toolkit: generation of 3D mammary tissue from diverse mammals

The varying pathways of mammary gland development across species and evolutionary history are underexplored, largely due to a lack of model systems. Recent progress in organoid technology holds the promise of enabling in-depth studies of the developmental adaptations that have occurred throughout the evolution of different species, fostering beneficial phenotypes. The practical application of this technology for mammary glands has been mostly confined to rodents and humans. In the current study, we have successfully created next-generation 3D mammary gland organoids from eight eutherian mammals and the first branched organoid of a marsupial mammary gland. Using mammary organoids, we identified a role for ROCK protein in regulating branching morphogenesis, a role that manifests differently in organoids from different mammals. This finding demonstrates the utility of the 3D organoid model for understanding the evolution and adaptations of signaling pathways. These achievements highlight the potential for organoid models to expand our understanding of mammary gland biology and evolution, and their potential utility in studies of lactation or breast cancer.

An Ancestral Major Histocompatibility Complex Organization in Cartilaginous Fish: Reconstructing MHC Origin and Evolution

Cartilaginous fish (sharks, rays, and chimeras) comprise the oldest living jawed vertebrates with a mammalian-like adaptive immune system based on immunoglobulins (Ig), T-cell receptors (TCRs), and the major histocompatibility complex (MHC). Here, we show that the cartilaginous fish "adaptive MHC" is highly regimented and compact, containing (i) a classical MHC class Ia (MHC-Ia) region containing antigen processing (antigen peptide transporters and immunoproteasome) and presenting (MHC-Ia) genes, (ii) an MHC class II (MHC-II) region (with alpha and beta genes) with linkage to beta-2-microglobulin (β2m) and bromodomain-containing 2, (iii) nonclassical MHC class Ib (MHC-Ib) regions with 450 million-year-old lineages, and (iv) a complement C4 associated with the MHC-Ia region. No MHC-Ib genes were found outside of the elasmobranch MHC. Our data suggest that both MHC-I and MHC-II genes arose after the second round of whole-genome duplication (2R) on a human chromosome (huchr) 6 precursor. Further analysis of MHC paralogous regions across early branching taxa from all jawed vertebrate lineages revealed that Ig/TCR genes likely arose on a precursor of the huchr9/12/14 MHC paralog. The β2m gene is linked to the Ig/TCR genes in some vertebrates suggesting that it was present at 1R, perhaps as the donor of C1 domain to the primordial MHC gene. In sum, extant cartilaginous fish exhibit a conserved and prototypical MHC genomic organization with features found in various vertebrates, reflecting the ancestral arrangement for the jawed vertebrates.

Major histocompatibility complex genomic investigation of endangered Chinese alligator provides insights into the evolution of tetrapod major histocompatibility complex and survival of critically bottlenecked species

Background

The major histocompatibility complex (MHC) gene family, a vital immune gene family in vertebrates, helps animals defend against pathogens. The polymorphism of MHC genes is important for a species and is considered to be caused by the numerous alleles of MHC antigen-presenting genes. However, the mechanism of this process is unclear due to the lack of data on the MHC structure. The evolutionary trajectories of the tetrapod MHC are also unclear because of insufficient studies on the reptile MHC architecture. Here, we studied the Chinese alligator (Alligator sinensis), which experienced a population bottleneck, but the population increased rapidly in the past 30 years and is proposed to have a unique MHC system to face pathogenic challenges.

Results

We successfully constructed a 2 Mb MHC region using bacterial artificial chromosome (BAC) library and genome data of the Chinese alligator and checked the antigen-presenting genes using transcriptome data and the rapid amplification of cDNA ends (RACE) technique. The MHC architecture reported here uncovers adjacent Class I and Class II subregions and a unique CD1 subregion. This newly added information suggested that the Class I-II structure pattern was more ancient in tetrapods and helped reconstruct the evolution of the MHC region architecture. We also found multiple groups of MHC class I (MHC-I) (12 duplicated loci, belonging to three groups, two of which were novel) and MHC class II (MHC-II) (11 duplicated loci, belonging to two groups) inside the 2 Mb MHC region, and there were three more duplicated MHC-I loci outside it. These highly duplicated antigen-presenting loci had differences in expression, amino acid length of antigen-presenting exons, and splice signal of exon and intron, which together promoted the polymorphism of duplicated genes. Although the MHC antigen-presenting genes were identified as monomorphic or oligomorphic in our previous population study, the loci with high copy numbers and many differences can make up for this loss, presenting another mechanism for polymorphism in antigen presentation. These MHC-I and MHC-IIB loci with low polymorphism for each locus, but high numbers in all, may also contribute to MHC antigen-presenting binding variability in a population.

Conclusion

To summarize, the fine MHC region architecture of reptiles presented in this study completes the evolutionary trajectories of the MHC structure in tetrapods, and these distinctive MHC gene groups in the Chinese alligator may have helped this species to expand rapidly in the past recent years.

Structure and evolution of the squamate major histocompatibility complex as revealed by two Anolis lizard genomes

The major histocompatibility complex (MHC) is an important genomic region for adaptive immunity and has long been studied in ecological and evolutionary contexts, such as disease resistance and mate and kin selection. The MHC has been investigated extensively in mammals and birds but far less so in squamate reptiles, the third major radiation of amniotes. We localized the core MHC genomic region in two squamate species, the green anole (Anolis carolinensis) and brown anole (A. sagrei), and provide the first detailed characterization of the squamate MHC, including the presence and ordering of known MHC genes in these species and comparative assessments of genomic structure and composition in MHC regions. We find that the Anolis MHC, located on chromosome 2 in both species, contains homologs of many previously-identified mammalian MHC genes in a single core MHC region. The repetitive element composition in anole MHC regions was similar to those observed in mammals but had important distinctions, such as higher proportions of DNA transposons. Moreover, longer introns and intergenic regions result in a much larger squamate MHC region (11.7 Mb and 24.6 Mb in the green and brown anole, respectively). Evolutionary analyses of MHC homologs of anoles and other representative amniotes uncovered generally monophyletic relationships between species-specific homologs and a loss of the peptide-binding domain exon 2 in one of two mhc2β gene homologs of each anole species. Signals of diversifying selection in each anole species was evident across codons of mhc1, many of which appear functionally relevant given known structures of this protein from the green anole, chicken, and human. Altogether, our investigation fills a major gap in understanding of amniote MHC diversity and evolution and provides an important foundation for future squamate-specific or vertebrate-wide investigations of the MHC.

Mathematical modeling and stochastic simulations suggest that low-affinity peptides can bisect MHC1-mediated export of high-affinity peptides into "early"- and "late"-phases

The peptide loading complex (PLC) is a multi-protein complex of the endoplasmic reticulum (ER) which optimizes major histocompatibility I (MHC1)-mediated export of intracellular high-affinity peptides. Whilst, the molecular biology of MHC1-mediated export is well supported by empirical data, the stoichiometry, kinetics and spatio-temporal profile of the participating molecular entities are a matter of considerable debate. Here, a low-affinity peptide-driven (LAPD)-model of MHC1-mediated high-affinity peptide export is formulated, implemented, analyzed and simulated. The model is parameterized in terms of the contribution of the shunt reaction to the concentration of exportable MHC1. Theoretical analyses and simulation studies of the model suggest that low-affinity peptides can bisect MHC1-mediated export of high-affinity peptides into time-dependent distinct “early”- and “late”-phases. The net exportable MHC1 (eM1β(t)) is a function of the retrograde (rM1β(t))- and anterograde (aM1β(t))-derived fractions. The “early”-phase is dominated by the contribution of the retrograde/recyclable (rM1β≈61%,aM1β≈39%) pathway to exportable MHC1, is characterized by Tapasin-mediated peptide-editing and is ATP-independent. The “late”-phase on the other hand, is characterized by de novo PLC-assembly, rapid disassembly and a significant contribution of the anterograde pathway to exportable MHC1 (rM1β≈21%,aM1β≈79%). The shunt reaction is rate limiting and may integrate peptide translocation with PLC-assembly/disassembly thereby, regulating peptide export under physiological and pathological (viral infections, dysplastic alterations) conditions.

Major histocompatibility complex genes and locus organization in the Komodo dragon (Varanus komodoensis)

We performed a meta-analysis of the newly assembled Komodo dragon (Varanus komodoensis) genome to characterize the major histocompatibility complex (MHC) of the species. The MHC gene clusters of the Komodo dragon are gene dense, complex, and contain counterparts of many genes of the human MHC. Our analysis identified 20 contigs encompassing ~ 6.9 Mbp of sequence with 223 annotated genes of which many are predicted orthologs to the genes of the human MHC. These MHC contigs range in size from 13.2 kb to 21.5 Mbp, contain an average of one gene per 30 kb, and are thought to occur on at least two chromosomes. Eight contigs, each > 100 kb, could be aligned to the human MHC based on gene content, and these represent gene clusters found in each of the recognized mammalian MHC subregions. The MHC of the Komodo dragon shares organizational features of other non-mammalian taxa. Multiple class Iα and class IIβ genes are indicated, with linkage between classical class I and immunoproteasome genes and between framework class I genes and genes associated with the mammalian class III subregion. These findings are supported in both Komodo genome assemblies and provide new insight into the MHC organization of these unique squamate reptiles.

Specimen-specific drift of densities defines distinct subclasses of extracellular vesicles from human whole saliva

Extracellular vesicles (EVs) in body fluids constitute heterogenous populations, which mirror their diverse parental cells as well as distinct EV-generation pathways. Various methodologies have been proposed to differentiate EVs in order to deepen the current understanding of EV biology. Equilibrium density-gradient centrifugation has often been used to separate EVs based on their buoyant densities; however, the standard conditions used for the method do not necessarily allow all EVs to move to their equilibrium density positions, which complicates the categorization of EVs. Here, by prolonging ultracentrifugation time to 96 h and fractionating EVs both by floating up or spinning down directions, we allowed 111 EV-associated protein markers from the whole saliva of three healthy volunteers to attain equilibrium. Interestingly, the determined buoyant densities of the markers drifted in a specimen-specific manner, and drift patterns differentiated EVs into at least two subclasses. One class carried classical exosomal markers, such as CD63 and CD81, and the other was characterized by the molecules involved in membrane remodeling or vesicle trafficking. Distinct patterns of density drift may represent the differences in generation pathways of EVs.

NLRC5 regulates expression of MHC-I and provides a target for anti-tumor immunity in transmissible cancers

Purpose

Downregulation of MHC class I (MHC-I) is a common immune evasion strategy of many cancers. Similarly, two allogeneic clonal transmissible cancers have killed thousands of wild Tasmanian devils (Sarcophilus harrisii) and also modulate MHC-I expression to evade anti-cancer and allograft responses. IFNG treatment restores MHC-I expression on devil facial tumor (DFT) cells but is insufficient to control tumor growth. Transcriptional co-activator NLRC5 is a master regulator of MHC-I in humans and mice but its role in transmissible cancers remains unknown. In this study, we explored the regulation and role of MHC-I in these unique genetically mis-matched tumors.

Methods

We used transcriptome and flow cytometric analyses to determine how MHC-I shapes allogeneic and anti-tumor responses. Cell lines that overexpress NLRC5 to drive antigen presentation, and B2M-knockout cell lines incapable of presenting antigen on MHC-I were used to probe the role of MHC-I in rare cases of tumor regressions.

Results

Transcriptomic results suggest that NLRC5 plays a major role in MHC-I regulation in devils. NLRC5 was shown to drive the expression of many components of the antigen presentation pathway but did not upregulate PDL1. Serum from devils with tumor regressions showed strong binding to IFNG-treated and NLRC5 cell lines; antibody binding to IFNG-treated and NRLC5 transgenic tumor cells was diminished or absent following B2M knockout.

Conclusion

MHC-I could be identified as a target for anti-tumor and allogeneic immunity. Consequently, NLRC5 could be a promising target for immunotherapy and vaccines to protect devils from transmissible cancers and inform development of transplant and cancer therapies for humans.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00432-021-03601-x.

Annotation of immune genes in the extinct thylacine (Thylacinus cynocephalus)

Advances in genome sequencing technology have enabled genomes of extinct species to be sequenced. However, given the fragmented nature of these genome assemblies, it is not clear whether it is possible to comprehensively annotate highly variable and repetitive genes such as those involved in immunity. As such, immune genes have only been investigated in a handful of extinct genomes, mainly in human lineages. In 2018 the genome of the thylacine (Thylacinus cynocephalus), a carnivorous marsupial from Tasmania that went extinct in 1936, was sequenced. Here we attempt to characterise the immune repertoire of the thylacine and determine similarity to its closest relative with a genome available, the Tasmanian devil (Sarcophilus harrisii), as well as other marsupials. Members from all major immune gene families were identified. However, variable regions could not be characterised, and complex families such as the major histocompatibility complex (MHC) were highly fragmented and located across multiple small scaffolds. As such, at a gene level we were unable to reconstruct full-length coding sequences for the majority of thylacine immune genes. Despite this, we identified genes encoding functionally important receptors and immune effector molecules, which suggests the functional capacity of the thylacine immune system was similar to other mammals. However, the high number of partial immune gene sequences identified limits our ability to reconstruct an accurate picture of the thylacine immune repertoire.

A chromosome-scale reference genome and genome-wide genetic variations elucidate adaptation in yak

Yak is an important livestock animal for the people indigenous to the harsh, oxygen‐limited Qinghai‐Tibetan Plateau and Hindu Kush ranges of the Himalayas. The yak genome was sequenced in 2012, but its assembly was fragmented because of the inherent limitations of the Illumina sequencing technology used to analyse it. An accurate and complete reference genome is essential for the study of genetic variations in this species. Long‐read sequences are more complete than their short‐read counterparts and have been successfully applied towards high‐quality genome assembly for various species. In this study, we present a high‐quality chromosome‐scale yak genome assembly (BosGru_PB_v1.0) constructed with long‐read sequencing and chromatin interaction technologies. Compared to an existing yak genome assembly (BosGru_v2.0), BosGru_PB_v1.0 shows substantially improved chromosome sequence continuity, reduced repetitive structure ambiguity, and gene model completeness. To characterize genetic variation in yak, we generated de novo genome assemblies based on Illumina short reads for seven recognized domestic yak breeds in Tibet and Sichuan and one wild yak from Hoh Xil. We compared these eight assemblies to the BosGru_PB_v1.0 genome, obtained a comprehensive map of yak genetic diversity at the whole‐genome level, and identified several protein‐coding genes absent from the BosGru_PB_v1.0 assembly. Despite the genetic bottleneck experienced by wild yak, their diversity was nonetheless higher than that of domestic yak. Here, we identified breed‐specific sequences and genes by whole‐genome alignment, which may facilitate yak breed identification.

Sequence and chromatin determinants of transcription factor binding and the establishment of cell type-specific binding patterns

Transcription factors (TFs) selectively bind distinct sets of sites in different cell types. Such cell type-specific binding specificity is expected to result from interplay between the TF's intrinsic sequence preferences, cooperative interactions with other regulatory proteins, and cell type-specific chromatin landscapes. Cell type-specific TF binding events are highly correlated with patterns of chromatin accessibility and active histone modifications in the same cell type. However, since concurrent chromatin may itself be a consequence of TF binding, chromatin landscapes measured prior to TF activation provide more useful insights into how cell type-specific TF binding events became established in the first place. Here, we review the various sequence and chromatin determinants of cell type-specific TF binding specificity. We identify the current challenges and opportunities associated with computational approaches to characterizing, imputing, and predicting cell type-specific TF binding patterns. We further focus on studies that characterize TF binding in dynamic regulatory settings, and we discuss how these studies are leading to a more complex and nuanced understanding of dynamic protein-DNA binding activities. We propose that TF binding activities at individual sites can be viewed along a two-dimensional continuum of local sequence and chromatin context. Under this view, cell type-specific TF binding activities may result from either strongly favorable sequence features or strongly favorable chromatin context.

The major histocompatibility complex of Old World camelids: Class I and class I-related genes

The genomic structure of the Major Histocompatibility Complex (MHC) region and variation in selected MHC class I related genes in Old World camels, Camelus bactrianus and Camelus dromedaries were studied. The overall genomic organization of the camel MHC region follows a general pattern observed in other mammalian species and individual MHC loci appear to be well conserved. Selected MHC class I genes B-67 and BL3-7 exhibited unexpectedly low variability, even when compared to other camel MHC class I related genes MR1 and MICA. Interspecific SNP and allele sharing are relatively common, and frequencies of heterozygotes are usually low. Such a low variation in a genomic region generally considered as one of the most polymorphic in vertebrate genomes is unusual. Evolutionary relationships between MHC class I related genes and their counterparts from other species seem to be rather complex. Often, they do not follow the general evolutionary history of the species concerned. Close evolutionary relationships of individual MHC class I loci between camels, humans and dogs were observed. Based on the results of this study and on our data on MHC class II genes, the extent and the pattern of polymorphism of the MHC region of Old World camelids differed from most mammalian groups studied so far. Camels thus seem to be an important model for our understanding of the role of genetic diversity in immune functions, especially in the context of unique features of their immunoglobulin and T-cell receptor genes.

MHC Class II variability in bare-nosed wombats (Vombatus ursinus)

Studies of gene diversity are used to investigate population dynamics, including immunological fitness. Aside from the selection of an appropriate gene target, an important factor that underpins these studies is the ability to recover viable DNA samples from native animals that are protected, threatened or difficult to sample or locate such as the bare-nosed wombat (Vombatus ursinus). In this study, we used genomic DNA extracted from muscle tissue samples and also identified the optimal method to extract DNA from fresh wombat scat samples to enable further analyses to be performed using non-invasive techniques. The DNA was probed via the polymerase chain reaction using previously targeted marsupial Major Histocompatibility Complex (MHC) gene primers. These genes are highly variable and associated with binding and presentation of pathogens in the immune system. Twenty-three unique MHC Class II DAB V. ursinus gene sequences were translated to 21 unique predicted peptide sequences from 34 individual tissue or scat samples. Vombatus ursinus MHC Class II DAB gene and peptide sequences were most similar to other marsupial DNA and peptide sequences. Further analysis also indicated the likelihood of MHC Class II DAB family membership through motif identification. Additional sampling is required to assess the full level of diversity of MHC Class II DAB genes among V. ursinus populations; however, this study is the first to identify MHC genes in a wombat and will advance immunological and disease studies of the species.

Immunogenetics of marsupial B-cells

Marsupials and eutherians are mammals that differ in their physiological traits, predominately their reproductive and developmental strategies; eutherians give birth to well-developed young, while marsupials are born highly altricial after a much shorter gestation. These developmental traits also result in differences in the development of the immune system of eutherian and marsupial species. In eutherians, B-cells are the key to humoral immunity as they are found in multiple lymphoid organs and have the unique ability to mediate the production of antigen-specific antibodies in the presence of extracellular pathogens. The development of B-cells in marsupials has been reported and hypothesised to be similar to that of eutherians, except that haematopoiesis occurs in the liver, postpartum, until the bone marrow fully matures. In eutherians, specific genes are linked to specific stages in B-cell development, maturation, and differentiation processes, and have been identified including immunoglobulins (heavy and light chains), cluster of differentiation markers (CD10, 19, 34 and CD79α/β), signal transduction molecules (BTK, Lyn and Syk) and transcriptional regulators (EBF1, E2A, and Pax5). This review aims to discuss the known similarities and differences between marsupial and eutherian B-cells, in regards to their genetic presence, homology, and developmental stages, as well as to highlight the areas requiring further investigation. By enhancing our understanding of the genes that are involved with B-cells in the marsupial lineage, it will, in turn, aid our understanding of the marsupial immune system and support the development of specific immunological reagents for research and wildlife conservation purposes.

Origin and evolution of the major histocompatibility complex class I region in eutherian mammals

Major histocompatibility complex (MHC) genes in vertebrates are vital in defending against pathogenic infections. To gain new insights into the evolution of MHC Class I (MHCI) genes and test competing hypotheses on the origin of the MHCI region in eutherian mammals, we studied available genome assemblies of nine species in Afrotheria, Xenarthra, and Laurasiatheria, and successfully characterized the MHCI region in six species. The following numbers of putatively functional genes were detected: in the elephant, four, one, and eight in the extended class I region, and κ and β duplication blocks, respectively; in the tenrec, one in the κ duplication block; and in the four bat species, one or two in the β duplication block. Our results indicate that MHCI genes in the κ and β duplication blocks may have originated in the common ancestor of eutherian mammals. In the elephant, tenrec, and all four bats, some MHCI genes occurred outside the MHCI region, suggesting that eutherians may have a more complex MHCI genomic organization than previously thought. Bat-specific three- or five-amino-acid insertions were detected in the MHCI α1 domain in all four bats studied, suggesting that pathogen defense in bats relies on MHCIs having a wider peptide-binding groove, as previously assayed by a bat MHCI gene with a three-amino-acid insertion showing a larger peptide repertoire than in other mammals. Our study adds to knowledge on the diversity of eutherian MHCI genes, which may have been shaped in a taxon-specific manner.

The Marine Mammal Class II Major Histocompatibility Complex Organization

Sirenians share with cetaceans and pinnipeds several convergent traits selected for the aquatic lifestyle. Living in water poses new challenges not only for locomotion and feeding but also for combating new pathogens, which may render the immune system one of the best tools aquatic mammals have for dealing with aquatic microbial threats. So far, only cetaceans have had their class II Major Histocompatibility Complex (MHC) organization characterized, despite the importance of MHC genes for adaptive immune responses. This study aims to characterize the organization of the marine mammal class II MHC using publicly available genomes. We located class II sequences in the genomes of one sirenian, four pinnipeds and eight cetaceans using NCBI-BLAST and reannotated the sequences using local BLAST search with exon and intron libraries. Scaffolds containing class II sequences were compared using dotplot analysis and introns were used for phylogenetic analysis. The manatee class II region shares overall synteny with other mammals, however most DR loci were translocated from the canonical location, past the extended class II region. Detailed analysis of the genomes of closely related taxa revealed that this presumed translocation is shared with all other living afrotherians. Other presumptive chromosome rearrangements in Afrotheria are the deletion of DQ loci in Afrosoricida and deletion of DP in E. telfairi. Pinnipeds share the main features of dog MHC: lack of a functional pair of DPA/DPB genes and inverted DRB locus between DQ and DO subregions. All cetaceans share the Cetartiodactyla inversion separating class II genes into two subregions: class IIa, with DR and DQ genes, and class IIb, with non-classic genes and a DRB pseudogene. These results point to three distinct and unheralded class II MHC structures in marine mammals: one canonical organization but lacking DP genes in pinnipeds; one bearing an inversion separating IIa and IIb subregions lacking DP genes found in cetaceans; and one with a translocation separating the most diverse class II gene from the MHC found in afrotherians and presumptive functional DR, DQ, and DP genes. Future functional research will reveal how these aquatic mammals cope with pathogen pressures with these divergent MHC organizations.

Evolution of Copy Number at the MHC Varies across the Avian Tree of Life

The evolution of the major histocompatibility complex (MHC) is shaped by frequent gene duplications and deletions, which generate extensive variation in the number of loci (gene copies) between different taxa. Here, we collected estimates of copy number at the MHC for over 250 bird species from 68 families. We found contrasting patterns of copy number evolution between MHC class I and class IIB, which encode receptors for intra- and extracellular pathogens, respectively. Across the avian evolutionary tree, there was evidence of accelerated evolution and stabilizing selection acting on copy number at class I, while copy number at class IIB was primarily influenced by fluctuating selection and drift. Reconstruction of MHC copy number variation showed ancestrally low numbers of MHC loci in nonpasserines and evolution toward larger numbers of loci in passerines. Different passerine lineages had the highest duplication rates for MHC class I (Sylvioidea) and class IIB (Muscicapoidea and Passeroidea). We also found support for the correlated evolution of MHC copy number and life-history traits such as lifespan and migratory behavior. These results suggest that MHC copy number evolution in birds has been driven by life histories and differences in exposure to intra- and extracellular pathogens.

Disentangling the mechanisms of mate choice in a captive koala population

Successful captive breeding programs are crucial to the long-term survival of many threatened species. However, pair incompatibility (breeding failure) limits sustainability of many captive populations. Understanding whether the drivers of this incompatibility are behavioral, genetic, or a combination of both, is crucial to improving breeding programs. We used 28 years of pairing data from the San Diego Zoo koala colony, plus genetic analyses using both major histocompatibility complex (MHC)-linked and non-MHC-linked microsatellite markers, to show that both genetic and non-genetic factors can influence mating success. Male age was reconfirmed to be a contributing factor to the likelihood of a koala pair copulating. This trend could also be related to a pair's age difference, which was highly correlated with male age in our dataset. Familiarity was reconfirmed to increase the probability of a successful copulation. Our data provided evidence that females select mates based on MHC and genome-wide similarity. Male heterozygosity at MHC class II loci was associated with both pre- and post-copulatory female choice. Genome-wide similarity, and similarity at the MHC class II DAB locus, were also associated with female choice at the post-copulatory level. Finally, certain MHC-linked alleles were associated with either increased or decreased mating success. We predict that utilizing a variety of behavioral and MHC-dependent mate choice mechanisms improves female fitness through increased reproductive success. This study highlights the complexity of mate choice mechanisms in a species, and the importance of ascertaining mate choice mechanisms to improve the success of captive breeding programs.

Koala immunology and infectious diseases: How much can the koala bear?

Infectious diseases are contributing to the decline of the iconic Australian marsupial, the koala (Phascolarctos cinereus). Infections with the obligate intracellular bacteria, Chlamydia pecorum, cause debilitating ocular and urogenital-tract disease while the koala-retrovirus (KoRV) has been implicated in host immunosuppression and exacerbation of chlamydial pathogenesis. Although histological studies have provided insight into the basic architecture of koala immune tissues, our understanding of the koala immune response to infectious disease has been limited, until recently, by a lack of species-specific immune reagents. Recent advances in the characterisation of key immune genes have focused on advancing our understanding of the immune response to Chlamydia infection, revealing commonalities in disease pathologies and immunity between koalas and other hosts and paving the way for the development of a koala Chlamydia vaccine. This review summarises these recent findings and highlights key aspects of the koala immune system requiring further attention with particular regard to their most prominent infectious diseases.

Weird Animals, Sex, and Genome Evolution

Making my career in Australia exposed me to the tyranny of distance, but it gave me opportunities to study our unique native fauna. Distantly related animal species present genetic variation that we can use to explore the most fundamental biological structures and processes. I have compared chromosomes and genomes of kangaroos and platypus, tiger snakes and emus, devils (Tasmanian) and dragons (lizards). I particularly love the challenges posed by sex chromosomes, which, apart from determining sex, provide stunning examples of epigenetic control and break all the evolutionary rules that we currently understand. Here I describe some of those amazing animals and the insights on genome structure, function, and evolution they have afforded us. I also describe my sometimes-random walk in science and the factors and people who influenced my direction. Being a woman in science is still not easy, and I hope others will find encouragement and empathy in my story.

Generalists and Specialists: A New View of How MHC Class I Molecules Fight Infectious Pathogens

In comparison with the major histocompatibility complexes (MHCs) of typical mammals, the chicken MHC is simple and compact with a single dominantly expressed class I molecule that can determine the immune response. In addition to providing useful information for the poultry industry and allowing insights into the evolution of the adaptive immune system, the simplicity of the chicken MHC has allowed the discovery of phenomena that are more difficult to discern in the more complicated mammalian systems. This review discusses the new concept that poorly expressed promiscuous class I alleles act as generalists to protect against a wide variety of infectious pathogens, while highly expressed fastidious class I alleles can act as specialists to protect against new and dangerous pathogens.

A broad overview of classical MHC I expression and bound peptides reveals an inverse correlation between repertoire breadth and cell-surface expression in some chicken and human alleles.

Several chicken class I alleles with wide peptide-binding repertoires (promiscuity) are associated with resistance to a variety of common diseases.

Conversely, a narrow peptide-binding repertoire (fastidiousness) in some human HLA-B alleles is associated with resistance to HIV progression.

Cell-surface expression of some classical class I alleles depends on the regulation of translocation to the cell surface rather than of transcription or translation. MHC translocation is influenced by peptide translocation in chickens and by tapasin interaction in humans.

Transcriptome sequencing of the long-nosed bandicoot (Perameles nasuta) reveals conservation and innovation of immune genes in the marsupial order Peramelemorphia

Bandicoots are omnivorous marsupials of the order Peramelemorphia. Conservation concerns and their unique biological characteristics suggest peramelomorphs are worthy research subjects, but knowledge of their genetics and immunology has lagged behind that of other high-profile marsupials. Here, we characterise the transcriptome of the long-nose bandicoot (Perameles nasuta), the first high-throughput data set from any peramelomorph. We investigate the immune gene repertoire of the bandicoot, with a focus on key immune gene families, and compare to previously characterised marsupial and mammalian species. We find that the immune gene complement in bandicoot is often conserved with respect to other marsupials; however, the diversity of expressed transcripts in several key families, such as major histocompatibility complex, T cell receptor μ and natural killer cell receptors, appears greater in the bandicoot than other Australian marsupials, including devil and koala. This transcriptome is an important first step for future studies of bandicoots and the bilby, allowing for population level analysis and construction of bandicoot-specific immunological reagents and assays. Such studies will be critical to understanding the immunology and physiology of Peramelemorphia and to inform the conservation of these unique marsupials.

Characterisation of MHC class I genes in the koala

Koala (Phascolarctos cinereus) populations are on the decline across the majority of Australia's mainland. Two major diseases threatening the long-term survival of affected koala populations are caused by obligate intracellular pathogens: Chlamydia and koala retrovirus (KoRV). To improve our understanding of the koala immune system, we characterised their major histocompatibility complex (MHC) class I genes, which are centrally involved in presenting foreign peptides derived from intracellular pathogens to cytotoxic T cells. A total of 11 class I genes were identified in the koala genome. Three genes, Phci-UA, UB and UC, showed relatively high genetic variability and were expressed in all 12 examined tissues, whereas the other eight genes had tissue-specific expression and limited polymorphism. Evidence of diversifying selection was detected in Phci-UA and UC, while gene conversion may have played a role in creating new alleles at Phci-UB. We propose that Phci-UA, UB and UC are likely classical MHC genes of koalas, and further research is needed to understand their role in koala chlamydial and KoRV infections.

Insights into the ancestral organisation of the mammalian MHC class II region from the genome of the pteropid bat, Pteropus alecto

BACKGROUND

Bats are an extremely successful group of mammals and possess a variety of unique characteristics, including their ability to co-exist with a diverse range of pathogens. The major histocompatibility complex (MHC) is the most gene dense and polymorphic region of the genome and MHC class II (MHC-II) molecules play a vital role in the presentation of antigens derived from extracellular pathogens and activation of the adaptive immune response. Characterisation of the MHC-II region of bats is crucial for understanding the evolution of the MHC and of the role of pathogens in shaping the immune system.

RESULTS

Here we describe the relatively contracted MHC-II region of the Australian black flying-fox (Pteropus alecto), providing the first detailed insight into the MHC-II region of any species of bat. Twelve MHC-II genes, including one locus (DRB2) located outside the class II region, were identified on a single scaffold in the bat genome. The presence of a class II locus outside the MHC-II region is atypical and provides evidence for an ancient class II duplication block. Two non-classical loci, DO and DM and two classical, DQ and DR loci, were identified in P. alecto. A putative classical, DPB pseudogene was also identified. The bat's antigen processing cluster, though contracted, remains highly conserved, thus supporting its importance in antigen presentation and disease resistance.

CONCLUSIONS

This detailed characterisation of the bat MHC-II region helps to fill a phylogenetic gap in the evolution of the mammalian class II region and is a stepping stone towards better understanding of the immune responses in bats to viral, bacterial, fungal and parasitic infections.

Evolutionary analysis of two complement C4 genes: Ancient duplication and conservation during jawed vertebrate evolution

The complement C4 is a thioester-containing protein, and a histidine (H) residue catalyzes the cleavage of the thioester to allow covalent binding to carbohydrates on target cells. Some mammalian and teleost species possess an additional isotype where the catalytic H is replaced by an aspartic acid (D), which binds preferentially to proteins. We found the two C4 isotypes in many other jawed vertebrates, including sharks and birds/reptiles. Phylogenetic analysis suggested that C4 gene duplication occurred in the early days of the jawed vertebrate evolution. The D-type C4 of bony fish except for mammals formed a cluster, termed D-lineage. The D-lineage genes were located in a syntenic region outside MHC, and evolved conservatively. Mammals lost the D-lineage before speciation, but D-type C4 was regenerated by recent gene duplication in some mammalian species or groups. Dual C4 molecules with different substrate specificities would have contributed to development of the antibody-dependent classical pathway.

High-throughput dynamic analysis of differentially expressed genes in splenic dendritic cells from mice infected with Schistosoma japonicum

Dendritic cells are the initiation and key point of immune response and play a role in immune regulation. So we explored the mechanisms involved in immune regulation of dendritic cells (DCs) against schistosomiasis using mice infected with Schistosoma japonicum. Splenic DCs from normal mice and mice with acute and chronic S. japonicum infection were sorted by flow cytometry. The numbers and functions of differentially expressed genes (DEGs) in DCs were determined by high-throughput analysis. All DEGs with transcription-level fold changes of ≥2 were selected and matched to corresponding genes in databases. Annotations and cluster analysis of DEGs were performed to compare differences between groups. Six important DEGs about immune regulation-CD86, TLR2, DC-SIGN, Capase3, PD-L2, and IL-7r were selected, and their transcription levels at different stages of schistosomisis were validated by qPCR. The Venn diagram of DEGs implied some genes are functional at all stages during S. japonicum infection, while others are only involved at certain stages. GO and KEGG pathway annotations indicated that these DEGs mainly belong to biological regulation, regulation of biological process, regulation of cellular process, antigen processing and presentation, cell adhesion molecules, cytokine-cytokine receptor interaction and Toll-like receptor signaling. Cluster analysis revealed immune regulation existed in splenic DCs. The results above indicated that the mechanisms underlying immune regulation to S. japonicum infection in mice are very complex. The present high-throughput dynamic analysis of DEGs in splenic DCs provides valuable insights into the molecular mechanisms underlying immune regulation in S. japonicum infection.

The UT family of MHC class I loci unique to non-eutherian mammals has limited polymorphism and tissue specific patterns of expression in the opossum

BACKGROUND

The Major Histocompatibility Complex (MHC) class I family of genes encode for molecules that have well-conserved structures, but have evolved to perform diverse functions. The availability of the gray, short-tailed opossum, Monodelphis domestica whole genome sequence has allowed for analysis of MHC class I gene content in this marsupial. Utilization of a novel method to search for MHC related domain structures revealed a previously unknown family of MHC class I-related genes. These genes, named UT1-17, are clustered on chromosome 1 in the opossum, unlinked to the MHC region. UT genes are only found in marsupial and monotreme genomes, consistent with being ancient in mammals yet lost in eutherian mammals. This study investigates the expression and polymorphism of the UT loci in the opossum to gain insight into their possible function.

RESULTS

Of the 17 opossum UT genes, most have restricted tissue transcription patterns, with the thymus and skin being the most common sites. Full-length structure of 11 UT transcripts revealed genes varying between five and eight exons, typical for class I family members. There were only two alternative splice variants found. The UT genes also have limited polymorphism and little evidence of positive selection. One locus, UT8, was chosen for further analysis due to its conservation amongst marsupials and generic characteristics. UT8 transcription is limited to developing αβ thymocytes, and is absent from mature αβ T cells in peripheral lymphoid tissues.

CONCLUSION

The overall characteristics and features of UT genes including low polymorphism and restricted tissue expression make it likely that the molecules encoded by UT genes perform roles other than antigenic peptide presentation.

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.

Co-evolution with chicken class I genes

The concept of co-evolution (or co-adaptation) has a long history, but application at molecular levels (e.g., 'supergenes' in genetics) is more recent, with a consensus definition still developing. One interesting example is the chicken major histocompatibility complex (MHC). In contrast to typical mammals that have many class I and class I-like genes, only two classical class I genes, two CD1 genes and some non-classical Rfp-Y genes are known in chicken, and all are found on the microchromosome that bears the MHC. Rarity of recombination between the closely linked and polymorphic genes encoding classical class I and TAPs allows co-evolution, leading to a single dominantly expressed class I molecule in each MHC haplotype, with strong functional consequences in terms of resistance to infectious pathogens. Chicken tapasin is highly polymorphic, but co-evolution with TAP and class I genes remains unclear. T-cell receptors, natural killer (NK) cell receptors, and CD8 co-receptor genes are found on non-MHC chromosomes, with some evidence for co-evolution of surface residues and number of genes along the avian and mammalian lineages. Over even longer periods, co-evolution has been invoked to explain how the adaptive immune system of jawed vertebrates arose from closely linked receptor, ligand, and antigen-processing genes in the primordial MHC.

Evolution and comparative analysis of the bat MHC-I region

Bats are natural hosts to numerous viruses and have ancient origins, having diverged from other eutherian mammals early in evolution. These characteristics place them in an important position to provide insights into the evolution of the mammalian immune system and antiviral immunity. We describe the first detailed partial map of a bat (Pteropus alecto) MHC-I region with comparative analysis of the MHC-I region and genes. The bat MHC-I region is highly condensed, yet relatively conserved in organisation, and is unusual in that MHC-I genes are present within only one of the three highly conserved class I duplication blocks. We hypothesise that MHC-I genes first originated in the β duplication block, and subsequently duplicated in a step-wise manner across the MHC-I region during mammalian evolution. Furthermore, bat MHC-I genes contain unique insertions within their peptide-binding grooves potentially affecting the peptide repertoire presented to T cells, which may have implications for the ability of bats to control infection without overt disease.

miR-628, a microRNA that is induced by Toll-like receptor stimulation, regulates porcine innate immune responses

Mammalian innate and acquired immune responses involve a coordinated, sequential, and self limiting sequence of events controlled by positive and negative regulatory mechanism. MicroRNAs have been implicated as a negative regulator for diverse biological events including immune responses. However, the involvement of miRNAs in regulating the immune responses is just beginning to be explored. Here, we characterized the expression profiling of 375 microRNAs in porcine monocytes induced by lipopolysaccharide (LPS), and result shows that several of them are endotoxin-responsive genes. Through promoter analysis, the miR-628 was found to be a NF-κB dependent gene. Importantly, miR-628 was predicted to base-pair with sequences in the 3'-UTR of the myeloid differentiation protein 88 (MyD88) gene. And we found that the UTR inhibit expression of a linked reporter gene coding a key adapter molecule downstream of Toll-like receptors (TLRs), resulting in suppressing of the TLR signaling. Therefore, we not only propose a role of miR-628 in control of the TLR signaling through a negative feedback regulation loop involving down-regulation of MyD88 protein levels, but results may also contribute to rational target selection orchestrating the inflammatory responses.

Marsupials and monotremes possess a novel family of MHC class I genes that is lost from the eutherian lineage

Background

Major histocompatibility complex (MHC) class I genes are found in the genomes of all jawed vertebrates. The evolution of this gene family is closely tied to the evolution of the vertebrate genome. Family members are frequently found in four paralogous regions, which were formed in two rounds of genome duplication in the early vertebrates, but in some species class Is have been subject to additional duplication or translocation, creating additional clusters. The gene family is traditionally grouped into two subtypes: classical MHC class I genes that are usually MHC-linked, highly polymorphic, expressed in a broad range of tissues and present endogenously-derived peptides to cytotoxic T-cells; and non-classical MHC class I genes generally have lower polymorphism, may have tissue-specific expression and have evolved to perform immune-related or non-immune functions. As immune genes can evolve rapidly and are subject to different selection pressure, we hypothesised that there may be divergent, as yet unannotated or uncharacterised class I genes.

Results

Application of a novel method of sensitive genome searching of available vertebrate genome sequences revealed a new, extensive sub-family of divergent MHC class I genes, denoted as UT, which has not previously been characterized. These class I genes are found in both American and Australian marsupials, and in monotremes, at an evolutionary chromosomal breakpoint, but are not present in non-mammalian genomes and have been lost from the eutherian lineage. We show that UT family members are expressed in the thymus of the gray short-tailed opossum and in other immune tissues of several Australian marsupials. Structural homology modelling shows that the proteins encoded by this family are predicted to have an open, though short, antigen-binding groove.

Conclusions

We have identified a novel sub-family of putatively non-classical MHC class I genes that are specific to marsupials and monotremes. This family was present in the ancestral mammal and is found in extant marsupials and monotremes, but has been lost from the eutherian lineage. The function of this family is as yet unknown, however, their predicted structure may be consistent with presentation of antigens to T-cells.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1745-4) contains supplementary material, which is available to authorized users.

Evolution of V genes from the TRV loci of mammals

Information concerning the evolution of T lymphocyte receptors (TCR) can be deciphered from that part of the molecule that recognizes antigen presented by major histocompatibility complex (MHC), namely the variable (V) regions. The genes that code for these variable regions are found within the TCR loci. Here, we describe a study of the evolutionary origin of V genes that code for the α and β chains of the TCR loci of mammals. In particular, we demonstrate that most of the 35 TRAV and 25 TRBV conserved genes found in Primates are also found in other Eutheria, while in Marsupials, Monotremes, and Reptiles, these genes diversified in a different manner. We also show that in mammals, all TRAV genes are derived from five ancestral genes, while all TRBV genes originate from four such genes. In Reptiles, the five TRAV and three out of the four TRBV ancestral genes exist, as well as other V genes not found in mammals. We also studied the TRGV and TRDV loci from all mammals, and we show a relationship of the TRDV to the TRAV locus throughout evolutionary time.

Characterisation of major histocompatibility complex class I genes at the fetal-maternal interface of marsupials

Major histocompatibility complex class I molecules (MHC-I) are expressed at the cell surface and are responsible for the presentation of self and non-self antigen repertoires to the immune system. Eutherian mammals express both classical and non-classical MHC-I molecules in the placenta, the latter of which are thought to modulate the maternal immune response during pregnancy. Marsupials last shared a common ancestor with eutherian mammals such as humans and mice over 160 million years ago. Since, like eutherians, they have an intra-uterine development dependent on a placenta, albeit a short-lived and less invasive one, they provide an opportunity to investigate the evolution of MHC-I expression at the fetal-maternal interface. We have characterised MHC-I mRNA expression in reproductive tissues of the tammar wallaby (Macropus eugenii) from the time of placental attachment to day 25 of the 26.5 day pregnancy. Putative classical MHC-I genes were expressed in the choriovitelline placenta, fetus, and gravid endometrium throughout the whole of this period. The MHC-I classical sequences were phylogenetically most similar to the Maeu-UC (50/100 clones) and Maeu-UA genes (7/100 clones). Expression of three non-classical MHC-I genes (Maeu-UD, Maeu-UK and Maeu-UM) were also present in placental samples. The results suggest that expression of classical and non-classical MHC-I genes in extant marsupial and eutherian mammals may have been necessary for the evolution of the ancestral therian placenta and survival of the mammalian fetus at the maternal-fetal interface.

Major Histocompatibility Complex Genes Map to Two Chromosomes in an Evolutionarily Ancient Reptile, the Tuatara Sphenodon punctatus

Major histocompatibility complex (MHC) genes are a central component of the vertebrate immune system and usually exist in a single genomic region. However, considerable differences in MHC organization and size exist between different vertebrate lineages. Reptiles occupy a key evolutionary position for understanding how variation in MHC structure evolved in vertebrates, but information on the structure of the MHC region in reptiles is limited. In this study, we investigate the organization and cytogenetic location of MHC genes in the tuatara (Sphenodon punctatus), the sole extant representative of the early-diverging reptilian order Rhynchocephalia. Sequencing and mapping of 12 clones containing class I and II MHC genes from a bacterial artificial chromosome library indicated that the core MHC region is located on chromosome 13q. However, duplication and translocation of MHC genes outside of the core region was evident, because additional class I MHC genes were located on chromosome 4p. We found a total of seven class I sequences and 11 class II β sequences, with evidence for duplication and pseudogenization of genes within the tuatara lineage. The tuatara MHC is characterized by high repeat content and low gene density compared with other species and we found no antigen processing or MHC framework genes on the MHC gene-containing clones. Our findings indicate substantial differences in MHC organization in tuatara compared with mammalian and avian MHCs and highlight the dynamic nature of the MHC. Further sequencing and annotation of tuatara and other reptile MHCs will determine if the tuatara MHC is representative of nonavian reptiles in general.

Immunomics of the koala (Phascolarctos cinereus)

The study of the koala transcriptome has the potential to advance our understanding of its immunome--immunological reaction of a given host to foreign antigens--and to help combat infectious diseases (e.g., chlamydiosis) that impede ongoing conservation efforts. We used Illumina sequencing of cDNA to characterize genes expressed in two different koala tissues of immunological importance, blood and spleen. We generated nearly 600 million raw sequence reads, and about 285 million of these were subsequently assembled and condensed into ~70,000 subcomponents that represent putative transcripts. We annotated ~16% of these subcomponents and identified those related to infection and the immune response, including Toll-like receptors (TLRs), RIG-I-like receptors (RLRs), major histocompatibility complex (MHC) genes, and koala retrovirus (KoRV). Using phylogenetic analyses, we identified 29 koala genes in these target categories and report their concordance with currently accepted gene groups. By mapping multiple sequencing reads to transcripts, we identified 56 putative SNPs in genes of interest. The distribution of these SNPs indicates that MHC genes (34 SNPs) are more diverse than KoRV (12 SNPs), TLRs (8 SNPs), or RLRs (2 SNPs). Our sequence data also indicate that KoRV sequences are highly expressed in the transcriptome. Our efforts have produced full-length sequences for potentially important immune genes in koala, which should serve as targets for future investigations that aim to conserve koala populations.

The opossum MHC genomic region revisited

The gray short-tailed opossum Monodelphis domestica is one of the few marsupial species for which a high quality whole genome sequence is available and the major histocompatibility complex (MHC) region has been annotated. Previous analyses revealed only a single locus within the opossum MHC region, designated Modo-UA1, with the features expected for encoding a functionally classical class I α-chain. Nine other class I genes found within the MHC are highly divergent and have features usually associated with non-classical roles. The original annotation, however, was based on an early version of the opossum genome assembly. More recent analyses of allelic variation in individual opossums revealed too many Modo-UA1 sequences per individual to be accounted for by a single MHC class I locus found in the genome assembly. A reanalysis of a later generation assembly, MonDom5, revealed the presence of two additional loci, now designated Modo-UA3 and UA4, in a region that was expanded and more complete than in the earlier assembly. Modo-UA1, UA3, and UA4 are all transcribed, although Modo-UA4 transcripts are rarer. Modo-UA4 is also relatively non-polymorphic. Evidence presented support the accuracy of the later assembly and the existence of three related class I genes in the opossum, making opossums more typical of mammals and most tetrapods by having multiple apparent classical MHC class I loci.

Comparative genome analyses reveal distinct structure in the saltwater crocodile MHC

The major histocompatibility complex (MHC) is a dynamic genome region with an essential role in the adaptive immunity of vertebrates, especially antigen presentation. The MHC is generally divided into subregions (classes I, II and III) containing genes of similar function across species, but with different gene number and organisation. Crocodylia (crocodilians) are widely distributed and represent an evolutionary distinct group among higher vertebrates, but the genomic organisation of MHC within this lineage has been largely unexplored. Here, we studied the MHC region of the saltwater crocodile (Crocodylus porosus) and compared it with that of other taxa. We characterised genomic clusters encompassing MHC class I and class II genes in the saltwater crocodile based on sequencing of bacterial artificial chromosomes. Six gene clusters spanning ∼452 kb were identified to contain nine MHC class I genes, six MHC class II genes, three TAP genes, and a TRIM gene. These MHC class I and class II genes were in separate scaffold regions and were greater in length (2-6 times longer) than their counterparts in well-studied fowl B loci, suggesting that the compaction of avian MHC occurred after the crocodilian-avian split. Comparative analyses between the saltwater crocodile MHC and that from the alligator and gharial showed large syntenic areas (>80% identity) with similar gene order. Comparisons with other vertebrates showed that the saltwater crocodile had MHC class I genes located along with TAP, consistent with birds studied. Linkage between MHC class I and TRIM39 observed in the saltwater crocodile resembled MHC in eutherians compared, but absent in avian MHC, suggesting that the saltwater crocodile MHC appears to have gene organisation intermediate between these two lineages. These observations suggest that the structure of the saltwater crocodile MHC, and other crocodilians, can help determine the MHC that was present in the ancestors of archosaurs.

Characterisation of non-classical MHC class I genes in the Tasmanian devil (Sarcophilus harrisii)

The Tasmanian devil (Sarcophilus harrisii) is a carnivorous marsupial that is under threat of extinction due to an unusual transmissible disease called Devil Facial Tumour Disease (DFTD). Previous studies on the classical MHC genes have provided important insights into immune responses in this endangered species; however, so far, very little is known about the non-classical MHC genes of this species, which can also play significant roles in the immune system. Here, we report characterisation of five non-classical class I genes in the Tasmanian devil, including Saha-UD, -UK, -UM, -MR1 and -CD1. Saha-UD has been isolated previously and is known to have low genetic polymorphism, though its categorisation as classical or non-classical gene has remained undetermined. In this study, we observed tissue-specific expression of Saha-UD, suggesting that it is more characteristic of a non-classical gene. Restricted tissue expression patterns were also observed for other genes, with an exception of Saha-MR1 being ubiquitously expressed in all examined tissues. Saha-UK, -UM and -MR1 were found to be genetically monomorphic, while four alleles were found at Saha-CD1 with signs of positive selection detected within the α1 domain. Among the four Saha-CD1 alleles, one predominant allele (Saha-CD1*01) showed a high allele frequency of 0.906 in the Tasmanian devil population, resulting in a low heterozygosity (0.188) at this locus. Alternative splicing takes place in Saha-CD1, giving rise to a full-length transcript and a splice variant lacking intact antigen-binding, β2m-binding, transmembrane and cytoplasmic domains.

Identification of MHCII variants associated with chlamydial disease in the koala (Phascolarctos cinereus)

Chlamydiosis, the most common infectious disease in koalas, can cause chronic urogenital tract fibrosis and infertility. High titres of serum immunoglobulin G against 10 kDa and 60 kDa chlamydial heat-shock proteins (c-hsp10 and c-hsp60) are associated with fibrous occlusion of the koala uterus and uterine tube. Murine and human studies have identified associations between specific major histocompatibility complex class II (MHCII) alleles or genotypes, and higher c-hsp 60 antibody levels or chlamydia-associated disease and infertility. In this study, we characterised partial MHCII DAB and DBB genes in female koalas (n = 94) from a single geographic population, and investigated associations among antibody responses to c-hsp60 quantified by ELISA, susceptibility to chlamydial infection, or age. The identification of three candidate MHCII variants provides additional support for the functional role of MHCII in the koala, and will inform more focused future studies. This is the first study to investigate an association between MHC genes with chlamydial pathogenesis in a non-model, free-ranging species.

Selection and trans-species polymorphism of major histocompatibility complex class II genes in the order Crocodylia

Major Histocompatibility Complex (MHC) class II genes encode for molecules that aid in the presentation of antigens to helper T cells. MHC characterisation within and between major vertebrate taxa has shed light on the evolutionary mechanisms shaping the diversity within this genomic region, though little characterisation has been performed within the Order Crocodylia. Here we investigate the extent and effect of selective pressures and trans-species polymorphism on MHC class II α and β evolution among 20 extant species of Crocodylia. Selection detection analyses showed that diversifying selection influenced MHC class II β diversity, whilst diversity within MHC class II α is the result of strong purifying selection. Comparison of translated sequences between species revealed the presence of twelve trans-species polymorphisms, some of which appear to be specific to the genera Crocodylus and Caiman. Phylogenetic reconstruction clustered MHC class II α sequences into two major clades representing the families Crocodilidae and Alligatoridae. However, no further subdivision within these clades was evident and, based on the observation that most MHC class II α sequences shared the same trans-species polymorphisms, it is possible that they correspond to the same gene lineage across species. In contrast, phylogenetic analyses of MHC class II β sequences showed a mixture of subclades containing sequences from Crocodilidae and/or Alligatoridae, illustrating orthologous relationships among those genes. Interestingly, two of the subclades containing sequences from both Crocodilidae and Alligatoridae shared specific trans-species polymorphisms, suggesting that they may belong to ancient lineages pre-dating the divergence of these two families from the common ancestor 85-90 million years ago. The results presented herein provide an immunogenetic resource that may be used to further assess MHC diversity and functionality in Crocodylia.

A study of immunomodulatory genes responses to macrophages of Schistosoma japonicum infection during different stages by microarray analysis

Macrophages initiate, modulate, and also serve as final effector cells in immune responses during the course of schistosomal infections. In this study, we investigated the gene expression profile and functional changes of macrophages in immune responses against the Schistosoma japonicum by microarray analysis. Hierarchical clustering analysis demonstrated that a significant switch in gene transformation associated with a type-1 response and linked with a type-2 cytokine phenotype occurs between 4.5 and 8 weeks post-infection. Moreover, the gene profiles at 3 later time-points following egg challenge were similar in complexity and magnitude. The data also showed that there were mostly inhibition of gene expression related TLR, IFN, MHC and TNFrsf at the switch between 4.5 and 8 weeks post-infection, It is suggested that these immunomodulatory genes may be down-regulated in defense against S. japonicum eggs and granuloma pathology. The induction of alternatively activated macrophage (AAMϕ) was important for dampening the inflammation in hepatic granulomas and contributing to a decrease in cytotoxicity. The gene expressions involved in repair/remodeling during liver fibrosis were also observed after egg production. Understanding the immune mechanisms associated with parasitic resistance, pathology of parasite infection, and parasite growth will provide useful insight on host-schistosome interactions and for the control of schistosomiasis.

HLA-DO as the optimizer of epitope selection for MHC class II antigen presentation

Processing of antigens for presentation to helper T cells by MHC class II involves HLA-DM (DM) and HLA-DO (DO) accessory molecules. A mechanistic understanding of DO in this process has been missing. The leading model on its function proposes that DO inhibits the effects of DM. To directly study DO functions, we designed a recombinant soluble DO and expressed it in insect cells. The kinetics of binding and dissociation of several peptides to HLA-DR1 (DR1) molecules in the presence of DM and DO were measured. We found that DO reduced binding of DR1 to some peptides, and enhanced the binding of some other peptides to DR1. Interestingly, these enhancing and reducing effects were observed in the presence, or absence, of DM. We found that peptides that were negatively affected by DO were DM-sensitive, whereas peptides that were enhanced by DO were DM-resistant. The positive and negative effects of DO could only be measured on binding kinetics as peptide dissociation kinetics were not affected by DO. Using Surface Plasmon Resonance, we demonstrate direct binding of DO to a peptide-receptive, but not a closed conformation of DR1. We propose that DO imposes another layer of control on epitope selection during antigen processing.