Characterisation of MHC class I genes in the koala

Using bioinformatics to investigate functional diversity: a case study of MHC diversity in koalas

Conservation genomics can greatly improve conservation outcomes of threatened populations, including those impacted by disease. Understanding diversity within immune gene families, including the major histocompatibility complex (MHC) and toll-like receptors (TLR), is important due to the role they play in disease resilience and susceptibility. With recent advancements in sequencing technologies and bioinformatic tools, the cost of generating high-quality sequence data has significantly decreased and made it possible to investigate diversity across entire gene families in large numbers of individuals compared to investigating only a few genes or a few populations previously. Here, we use the koala as a case study for investigating functional diversity across populations. We utilised previous target enrichment data and 438 whole genomes to firstly, determine the level of sequencing depth required to investigate MHC diversity and, secondly, determine the current level of diversity in MHC genes in koala populations. We determined for low complexity, conserved genes such as TLR genes 10 × sequencing depth is sufficient to reliably genotype more than 90% of variants, whereas for complex genes such as the MHC greater than 20 × and preferably 30 × sequencing depth is required. We used whole genome data to identify 270 biallelic SNPs across 24 MHC genes as well as copy number variation (CNV) within class I and class II genes and conduct supertype analysis. Overall, we have provided a bioinformatic workflow for investigating variation in a complex immune gene family from whole genome sequencing data and determined current levels of diversity within koala MHC genes.

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.

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.

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.

The Koala Immune Response to Chlamydial Infection and Vaccine Development-Advancing Our Immunological Understanding

Simple Summary

Chlamydia is a major pathogen of the Australian marsupial, the koala (Phascolarctos cinereus). One approach to improving this situation is to develop a vaccine. Human Chlamydia research suggests that an effective anti-chlamydial response will involve a balance between a cell-mediated Th1 response and a humoral Th2 responses, involving systemic IgG and mucosal IgA. Characterization of koalas with chlamydial disease suggests that increased expression for similar immunological pathways and monitoring of koalas’ post-vaccination can be successful and subsequently lead to improved vaccines. These findings offer optimism that a chlamydial vaccine for wider distribution to koalas is not far off.

Abstract

Chlamydia is a significant pathogen for many species, including the much-loved Australian marsupial, the koala (Phascolarctos cinereus). To combat this situation, focused research has gone into the development and refinement of a chlamydial vaccine for koalas. The foundation of this process has involved characterising the immune response of koalas to both natural chlamydial infection as well as vaccination. From parallels in human and mouse research, it is well-established that an effective anti-chlamydial response will involve a balance of cell-mediated Th1 responses involving interferon-gamma (IFN-γ), humoral Th2 responses involving systemic IgG and mucosal IgA, and inflammatory Th17 responses involving interleukin 17 (IL-17) and neutrophils. Characterisation of koalas with chlamydial disease has shown increased expression within all three of these major immunological pathways and monitoring of koalas’ post-vaccination has detected further enhancements to these key pathways. These findings offer optimism that a chlamydial vaccine for wider distribution to koalas is not far off. Recent advances in marsupial genetic knowledge and general nucleic acid assay technology have moved koala immunological research a step closer to other mammalian research systems. However, koala-specific reagents to directly assay cytokine levels and cell-surface markers are still needed to progress our understanding of koala immunology.

High-throughput immunogenetic typing of koalas suggests possible link between MHC alleles and cancers

Characterizing the allelic diversity within major histocompatibility complex (MHC) genes is an important way of determining the potential genetic resilience of a population to infectious and ecological pressures. For the koala (Phascolarctos cinereus), endemic diseases, anthropogenic factors and climate change are all placing increased pressure on this vulnerable marsupial. To increase the ability of researchers to study MHC genetics in koalas, this study developed and tested a high-throughput immunogenetic profiling methodology for targeting MHC class I UA and UC genes and MHC class II DAB, DBB, DCB and DMB genes in a population of 82 captive koalas. This approach was validated by comparing the determined allelic profiles from 36 koala family units (18 dam-sire-joey units and 18 parent-joey pairs), finding 96% overall congruence within family profiles. Cancers are a significant cause of morbidity in koalas and the risk factors remain undetermined. Our analysis of this captive population revealed several novel MHC alleles, including a potential link between the DBB*03 allele and a risk of developing cancer. This method offers a reliable, high-throughput protocol for expanded study into koala immunogenetics.

Marsupial chromosomics: bridging the gap between genomes and chromosomes

Marsupials have unique features that make them particularly interesting to study, and sequencing of marsupial genomes is helping to understand their evolution. A decade ago, it was a huge feat to sequence the first marsupial genome. Now, the advances in sequencing technology have made the sequencing of many more marsupial genomes possible. However, the DNA sequence is only one component of the structures it is packaged into: chromosomes. Knowing the arrangement of the DNA sequence on each chromosome is essential for a genome assembly to be used to its full potential. The importance of combining sequence information with cytogenetics has previously been demonstrated for rapidly evolving regions of the genome, such as the sex chromosomes, as well as for reconstructing the ancestral marsupial karyotype and understanding the chromosome rearrangements involved in the Tasmanian devil facial tumour disease. Despite the recent advances in sequencing technology assisting in genome assembly, physical anchoring of the sequence to chromosomes is required to achieve a chromosome-level assembly. Once chromosome-level assemblies are achieved for more marsupials, we will be able to investigate changes in the packaging and interactions between chromosomes to gain an understanding of the role genome architecture has played during marsupial evolution.

Koala immunogenetics and chlamydial strain type are more directly involved in chlamydial disease progression in koalas from two south east Queensland koala populations than koala retrovirus subtypes

Chlamydial disease control is increasingly utilised as a management tool to stabilise declining koala populations, and yet we have a limited understanding of the factors that contribute to disease progression. To examine the impact of host and pathogen genetics, we selected two geographically separated south east Queensland koala populations, differentially affected by chlamydial disease, and analysed koala major histocompatibility complex (MHC) genes, circulating strains of Chlamydia pecorum and koala retrovirus (KoRV) subtypes in longitudinally sampled, well-defined clinical groups. We found that koala immunogenetics and chlamydial genotypes differed between the populations. Disease progression was associated with specific MHC alleles, and we identified two putative susceptibility (DCb 03, DBb 04) and protective (DAb 10, UC 01:01) variants. Chlamydial genotypes belonging to both Multi-Locus Sequence Typing sequence type (ST) 69 and ompA genotype F were associated with disease progression, whereas ST 281 was associated with the absence of disease. We also detected different ompA genotypes, but not different STs, when long-term infections were monitored over time. By comparison, KoRV profiles were not significantly associated with disease progression. These findings suggest that chlamydial genotypes vary in pathogenicity and that koala immunogenetics and chlamydial strains are more directly involved in disease progression than KoRV subtypes.

Therapeutic vaccination of koalas harbouring endogenous koala retrovirus (KoRV) improves antibody responses and reduces circulating viral load

The long-term survival of the koala is under serious threat from multiple factors, including infectious disease agents such as Chlamydia and koala retrovirus (KoRV). KoRV is present in both exogenous and endogenous forms, depending on the geographical location of the population. In the northern half of Australia, it is present as an endogenous infection in all koalas, making a case for an urgent need to develop a therapeutic vaccine that might prevent KoRV-associated pathologies in these koalas. To this end, we determined the therapeutic effects of vaccinating koalas harbouring endogenous KoRV with a recombinant KoRV Env protein combined with a Tri-adjuvant. We found that vaccination led to a significant increase in circulating anti-KoRV IgG levels, as well as increase in neutralising antibodies. Our study also showed that post-vaccination antibodies were able to recognize epitopes on the Env protein that were unrecognised pre-vaccination, as well as resulting in an increase in the recognition of the previously recognised epitopes. The vaccine also induced antibodies that were cross-reactive against multiple KoRV-subtypes. Finally, we found a complete clearance of KoRV-A in plasma from koalas that had detectable levels of KoRV-A pre-vaccination. Similarly, there was a significant reduction in the expression of KoRV-B viral RNA levels post-vaccination. Collectively, this study showed that koalas harbouring endogenous KoRV can benefit from prophylactic vaccination against KoRV using a recombinant KoRV-A Env protein and that the mechanism of this protection might be through the boosting of natural anti-KoRV antibodies and expanding the breadth of the recognised epitopes.

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.

Adaptation and conservation insights from the koala genome

The koala, the only extant species of the marsupial family Phascolarctidae, is classified as 'vulnerable' due to habitat loss and widespread disease. We sequenced the koala genome, producing a complete and contiguous marsupial reference genome, including centromeres. We reveal that the koala's ability to detoxify eucalypt foliage may be due to expansions within a cytochrome P450 gene family, and its ability to smell, taste and moderate ingestion of plant secondary metabolites may be due to expansions in the vomeronasal and taste receptors. We characterized novel lactation proteins that protect young in the pouch and annotated immune genes important for response to chlamydial disease. Historical demography showed a substantial population crash coincident with the decline of Australian megafauna, while contemporary populations had biogeographic boundaries and increased inbreeding in populations affected by historic translocations. We identified genetically diverse populations that require habitat corridors and instituting of translocation programs to aid the koala's survival in the wild.