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.

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.

Helping koalas battle disease - Recent advances in Chlamydia and koala retrovirus (KoRV) disease understanding and treatment in koalas

The iconic Australian marsupial, the koala (Phascolarctos cinereus), has suffered dramatic population declines as a result of habitat loss and fragmentation, disease, vehicle collision mortality, dog attacks, bushfires and climate change. In 2012, koalas were officially declared vulnerable by the Australian government and listed as a threatened species. In response, research into diseases affecting koalas has expanded rapidly. The two major pathogens affecting koalas are Chlamydia pecorum, leading to chlamydial disease and koala retrovirus (KoRV). In the last eight years, these pathogens and their diseases have received focused study regarding their sources, genetics, prevalence, disease presentation and transmission. This has led to vast improvements in pathogen detection and treatment, including the ongoing development of vaccines for each as a management and control strategy. This review will summarize and highlight the important advances made in understanding and combating C. pecorum and KoRV in koalas, since they were declared a threatened species. With complementary advances having also been made from the koala genome sequence and in our understanding of the koala immune system, we are primed to make a significant positive impact on koala health into the future.

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.

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.

Host stress physiology and Trypanosoma haemoparasite infection influence innate immunity in the woylie (Bettongia penicillata)

Understanding immune function is critical to conserving wildlife in view of infectious disease threats, particularly in threatened species vulnerable to stress, immunocompromise and infection. However, few studies examine stress, immune function and infection in wildlife. We used a flow cytometry protocol developed for human infants to assess phagocytosis, a key component of innate immunity, in a critically endangered marsupial, the woylie (Bettongia penicillata). The effects of stress physiology and Trypanosoma infection on phagocytosis were investigated. Blood and faecal samples were collected from woylies in a captive facility over three months. Trypanosoma status was determined using PCR. Faecal cortisol metabolites (FCM) were quantified by enzyme-immunoassay. Mean phagocytosis measured was >90%. An interaction between sex and FCM influenced the percentage of phagocytosing leukocytes, possibly reflecting the influence of sex hormones and glucocorticoids. An interaction between Trypanosoma status and FCM influenced phagocytosis index, suggesting that stress physiology and infection status influence innate immunity.

Identification and analysis of divergent immune gene families within the Tasmanian devil genome

Background

The Tasmanian devil (Sarcophilus harrisii) is being threatened with extinction in the wild by a disease known as devil facial tumour disease (DFTD). In order to prevent the spread of this disease a thorough understanding of the Tasmanian devil immune system and its response to the disease is required. In 2011 and 2012 two genome sequencing projects of the Tasmania devil were released. This has provided us with the raw data required to begin to investigate the Tasmanian devil immunome in depth. In this study we characterise immune gene families of the Tasmanian devil. We focus on immunoglobulins, T cell receptors and cytokine families.

Results

We identify and describe 119 cytokines including 40 interleukins, 39 chemokines, 8 interferons, 18 tumour necrosis family cytokines and 14 additional cytokines. Constant regions for immunoglobulins and T cell receptors were also identified. The repertoire of genes in these families was similar to the opossum, however devil specific duplications were seen and orthologs to eutherian genes not previously identified in any marsupial were also identified.

Conclusions

By using multiple data sources as well as targeted search methods, highly divergent genes across the Tasmanian devil immune system were identified and characterised. This understanding will allow for the development of devil specific assays and reagents and allow for future studies into the immune response of the Tasmanian devil immune system to DFTD.

Electronic supplementary material

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

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSION

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

SNP marker discovery in koala TLR genes

Toll-like receptors (TLRs) play a crucial role in the early defence against invading pathogens, yet our understanding of TLRs in marsupial immunity is limited. Here, we describe the characterisation of nine TLRs from a koala immune tissue transcriptome and one TLR from a draft sequence of the koala genome and the subsequent development of an assay to study genetic diversity in these genes. We surveyed genetic diversity in 20 koalas from New South Wales, Australia and showed that one gene, TLR10 is monomorphic, while the other nine TLR genes have between two and 12 alleles. 40 SNPs (16 non-synonymous) were identified across the ten TLR genes. These markers provide a springboard to future studies on innate immunity in the koala, a species under threat from two major infectious diseases.

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.

Comparison of subcutaneous versus intranasal immunization of male koalas (Phascolarctos cinereus) for induction of mucosal and systemic immunity against Chlamydia pecorum

Chlamydia pecorum infections are debilitating in the koala, contributing significantly to morbidity and mortality, with current antibiotic treatments having minimal success and adversely affecting gut microflora. This, combined with the sometimes-asymptomatic nature of the infection, suggests that an efficacious anti-chlamydial vaccine is required to control chlamydial infections in the koala. To date vaccination studies have focused primarily on female koalas, however, given the physiological differences between male and female reproductive tracts, we tested the efficacy of a vaccine in 12 captive male koalas. We evaluated the potential of both subcutaneous and intranasal vaccine delivery to elicit mucosal immunity in male koalas. Our results showed that both intranasal and subcutaneous delivery of a vaccine consisting of C. pecorum major outer membrane protein (MOMP) and the adjuvant immunostimulating complex (ISC) induced significant immune responses in male koalas. Subcutaneous immunization elicited stronger cell-mediated responses in peripheral blood lymphocytes (PBL), and greater plasma antibody levels whereas the intranasal immunization elicited stronger humoral responses in urogenital tract (UGT) secretions. This is the first time a Chlamydia vaccine has been tested in the male koala and the first assessment of a mucosal vaccination route in this species. Our results suggest that vaccination of male koalas can elicit mucosal immunity and could contribute to the long-term survivability of wild populations of the koala.

IL12Rβ1: the cytokine receptor that we used to know

Human IL12RB1 encodes IL12Rβ1, a type I transmembrane receptor that is an essential component of the IL12- and IL23-signaling complex. IL12RB1 is well-established as being a promoter of delayed type hypersensitivity (DTH), the immunological reaction that limits tuberculosis. However, recent data demonstrate that in addition to promoting DTH, IL12RB1 also promotes autoimmunity. The contradictory roles of IL12RB1 in human health raises the question, what are the factors governing IL12RB1 function in a given individual, and how is inter-individual variability in IL12RB1 function introduced? Here we review recent data that demonstrate individual variability in IL12RB1 function is introduced at the epigenetic, genomic polymorphism, and mRNA splicing levels. Where and how these differences contribute to disease susceptibility and outcome are also reviewed. Collectively, recent data support a model wherein IL12RB1 sequence variability - whether introduced at the genomic or post-transcriptional level - contributes to disease, and that human IL12RB1 is not as simple a gene as we once believed.

Interleukin 17A is an immune marker for chlamydial disease severity and pathogenesis in the koala (Phascolarctos cinereus)

The koala (Phascolarctos cinereus) is an iconic Australian marsupial species that is facing many threats to its survival. Chlamydia pecorum infections are a significant contributor to this ongoing decline. A major limiting factor in our ability to manage and control chlamydial disease in koalas is a limited understanding of the koala's cell-mediated immune response to infections by this bacterial pathogen. To identify immunological markers associated with chlamydial infection and disease in koalas, we used koala-specific Quantitative Real Time PCR (qrtPCR) assays to profile the cytokine responses of Peripheral Blood Mononuclear Cells (PBMCs) collected from 41 koalas with different stages of chlamydial disease. Target cytokines included the principal Th1 (Interferon gamma; IFNγ), Th2 (Interleukin 10; IL10), and pro-inflammatory cytokines (Tumor Necrosis Factor alpha; TNFα). A novel koala-specific IL17A qrtPCR assay was also developed as part of this study to quantitate the gene expression of this Th17 cytokine in koalas. A statistically significant higher IL17A gene expression was observed in animals with current chlamydial disease compared to animals with asymptomatic chlamydial infection. A modest up-regulation of pro-inflammatory cytokines, such as TNFα and IFNγ, was also observed in these animals with signs of current chlamydial disease. IL10 gene expression was not evident in the majority of animals from both groups. Future longitudinal studies are now required to confirm the role played by cytokines in pathology and/or protection against C. pecorum infection in the koala.