MHC class II diversity of koala (Phascolarctos cinereus) populations across their range

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.

Immunohistochemical characterization of the immune cell response during chlamydial infection in the male and female koala (Phascolarctos cinereus) reproductive tract

Chlamydiosis is one of the main causes of the progressive decline of koala populations in eastern Australia. While histologic, immunologic, and molecular studies have provided insights into the basic function of the koala immune system, the in situ immune cell signatures during chlamydial infection of the reproductive tract in koalas have not been investigated. Thirty-two female koalas and 47 males presented to wildlife hospitals with clinical signs suggestive of Chlamydia infection were euthanized with the entire reproductive tract collected for histology; immunohistochemistry (IHC) for T-cell (CD3ε, CD4, and CD8α), B-cell (CD79b), and human leukocyte antigen (HLA)-DR markers; and quantitative real-time polymerase chain reaction (rtPCR) for Chlamydia pecorum. T-cells, B-cells, and HLA-DR-positive cells were observed in both the lower and upper reproductive tracts of male and female koalas with a statistically significant associations between the degree of the inflammatory reaction; the number of CD3, CD4, CD79b, and HLA-DR positive cells; and the PCR load. CD4-positive cells were negatively associated with the severity of the gross lesions. The distribution of immune cells was also variable according to the location within the genital tract in both male and female koalas. These preliminary results represent a step forward towards further exploring mechanisms behind chlamydial infection immunopathogenesis, thus providing valuable information about the immune response and infectious diseases in free-ranging koalas.

Geographic patterns of koala retrovirus genetic diversity, endogenization, and subtype distributions

Koala retrovirus (KoRV) subtype A (KoRV-A) is currently in transition from exogenous virus to endogenous viral element, providing an ideal system to elucidate retroviral-host coevolution. We characterized KoRV geography using fecal DNA from 192 samples across 20 populations throughout the koala's range. We reveal an abrupt change in KoRV genetics and incidence at the Victoria/New South Wales state border. In northern koalas, pol gene copies were ubiquitously present at above five per cell, consistent with endogenous KoRV. In southern koalas, pol copies were detected in only 25.8% of koalas and always at copy numbers below one, while the env gene was detected in all animals and in a majority at copy numbers above one per cell. These results suggest that southern koalas carry partial endogenous KoRV-like sequences. Deep sequencing of the env hypervariable region revealed three putatively endogenous KoRV-A sequences in northern koalas and a single, distinct sequence present in all southern koalas. Among northern populations, env sequence diversity decreased with distance from the equator, suggesting infectious KoRV-A invaded the koala genome in northern Australia and then spread south. The exogenous KoRV subtypes (B to K), two novel subtypes, and intermediate subtypes were detected in all northern koala populations but were strikingly absent from all southern animals tested. Apart from KoRV subtype D, these exogenous subtypes were generally locally prevalent but geographically restricted, producing KoRV genetic differentiation among northern populations. This suggests that sporadic evolution and local transmission of the exogenous subtypes have occurred within northern Australia, but this has not extended into animals within southern Australia.

Differences in constitutive innate immunity between divergent Australian marsupials

Understanding immunity in wildlife populations is important from both One Health and conservation perspectives. The constitutive innate immune system is the first line of defence against pathogens, and comparisons among taxa can test the impact of evolution and life history on immune function. We investigated serum bacterial killing ability (BKA) of five marsupial species that employ varying life history strategies, demonstrated to influence immunity in other vertebrates. The brushtail possum and eastern grey kangaroo had the greatest BKA, while ringtail possums and koalas had the least. These differences were independent of social structure, captivity status and phylogeny, but were associated with diet and body size. Sex and disease status had no effect on BKA in koalas, however potential for differences between wild and captive koalas warrants further investigation. The current study has provided a foundation for future investigations into how adaptive and innate immunity interact in marsupials from an eco-evolutionary perspective.

Transcriptomic and genomic variants between koala populations reveals underlying genetic components to disorders in a bottlenecked population

Historical hunting pressures on koalas in the southern part of their range in Australia have led to a marked genetic bottleneck when compared with their northern counterparts. There are a range of suspected genetic disorders such as testicular abnormalities, oxalate nephrosis and microcephaly reported at higher prevalence in these genetically restricted southern animals. This paper reports analysis of differential expression of genes from RNAseq of lymph nodes, SNPs present in genes and the fixation index (population differentiation due to genetic structure) of these SNPs from two populations, one in south east Queensland, representative of the northern genotype and one in the Mount Lofty Ranges South Australia, representative of the southern genotype. SNPs that differ between these two populations were significantly enriched in genes associated with brain diseases. Genes which were differentially expressed between the two populations included many associated with brain development or disease, and in addition a number associated with testicular development, including the androgen receptor. Finally, one of the 8 genes both differentially expressed and with a statistical difference in SNP frequency between populations was SLC26A6 (solute carrier family 26 member 6), an anion transporter that was upregulated in SA koalas and is associated with oxalate transport and calcium oxalate uroliths in humans. Together the differences in SNPs and gene expression described in this paper suggest an underlying genetic basis for several disorders commonly seen in southern Australian koalas, supporting the need for further research into the genetic basis of these conditions, and highlighting that genetic selection in managed populations may need to be considered in the future.

Incidence of sarcoptic mange in bare-nosed wombats (Vombatus ursinus)

Sarcoptic mange is a threat to bare-nosed wombats (Vombatus ursinus) in Australia and a major factor contributing to the decline in populations of this species. It is caused by a mite, Sarcoptes scabiei, and if untreated can lead to severe irritation, blindness, starvation and eventual death. Modes of transmission likely occur through direct contact between wombats and indirectly from wombat burrows. Our study aimed to estimate and compare the numbers of bare-nosed wombats at three study sites, and to determine how rainfall, temperature and humidity influenced the incidence of sarcoptic mange in the populations. The abundances of macropods and rabbits were also compared with wombat abundance and sarcoptic mange prevalence at the three study sites. Across the study sites, 1655 bare-nosed wombats were observed. Sarcoptic mange prevalence was between 7.0% and 40.7%, depending on site, season and year. Sarcoptic mange prevalence correlated with yearly rainfall, with a higher incidence of sarcoptic mange occurring in higher-rainfall years. Higher numbers of macropods in the study areas also correlated with higher incidences of sarcoptic mange in the wombat populations.

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.

Prevalence of cryptococcal antigenemia and nasal colonization in a free-ranging koala population

Cryptococcosis, caused by environmental fungi in the Cryptococcus neoformans and Cryptococcus gattii species complexes, affects a variety of hosts, including koalas (Phascolarctos cinereus). Cryptococcal antigenemia and nasal colonization are well characterized in captive koalas, but free-ranging populations have not been studied systematically. Free-ranging koalas (181) from the Liverpool Plains region of New South Wales, Australia, were tested for cryptococcal antigenemia (lateral flow immunoassay) and nasal colonization (bird seed agar culture). Results were related to environmental and individual koala characteristics. Eucalypt trees (14) were also randomly tested for the presence of Cryptococcus spp. by bird seed agar culture. In sum, 5.5% (10/181) and 6.6% (12/181) of koalas were positive for antigenemia and nasal colonization, respectively, on at least one occasion. And 64.3% (9/14) of eucalypts were culture-positive for Cryptococcus spp. URA5 restriction fragment length polymorphism analysis identified most isolates as C. gattii VGI, while C. neoformans VNI was only found in one koala and one tree. Colonized koalas were significantly more likely to test positive for antigenemia. No associations between antigenemia or colonization, and external environmental characteristics (the relative abundance of Eucalyptus camaldulensis and season), or individual koala characteristics (body condition, sex, and age), could be established, suggesting that antigenemia and colonization are random outcomes of host-pathogen-environment interactions. The relationship between positive antigenemia status and a relatively high abundance of E. camaldulensis requires further investigation. This study characterizes cryptococcosis in a free-ranging koala population, expands the ecological niche of the C. gattii/C. neoformans species complexes and highlights free-ranging koalas as important sentinels for this disease.

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.

Genomic comparisons reveal biogeographic and anthropogenic impacts in the koala (Phascolarctos cinereus): a dietary-specialist species distributed across heterogeneous environments

The Australian koala is an iconic marsupial with highly specific dietary requirements distributed across heterogeneous environments, over a large geographic range. The distribution and genetic structure of koala populations has been heavily influenced by human actions, specifically habitat modification, hunting and translocation of koalas. There is currently limited information on population diversity and gene flow at a species-wide scale, or with consideration to the potential impacts of local adaptation. Using species-wide sampling across heterogeneous environments, and high-density genome-wide markers (SNPs and PAVs), we show that most koala populations display levels of diversity comparable to other outbred species, except for those populations impacted by population reductions. Genetic clustering analysis and phylogenetic reconstruction reveals a lack of support for current taxonomic classification of three koala subspecies, with only a single evolutionary significant unit supported. Furthermore, ~70% of genetic variance is accounted for at the individual level. The Sydney Basin region is highlighted as a unique reservoir of genetic diversity, having higher diversity levels (i.e., Blue Mountains region; AvHecorr=0.20, PL% = 68.6). Broad-scale population differentiation is primarily driven by an isolation by distance genetic structure model (49% of genetic variance), with clinal local adaptation corresponding to habitat bioregions. Signatures of selection were detected between bioregions, with no single region returning evidence of strong selection. The results of this study show that although the koala is widely considered to be a dietary-specialist species, this apparent specialisation has not limited the koala's ability to maintain gene flow and adapt across divergent environments as long as the required food source is available.

Demographic, environmental and genetic determinants of mating success in captive koalas (Phascolarctos cinereus)

Many factors have been shown to affect mating behavior. For instance, genes of the major histocompatibility complex (MHC) are known to influence mate choice in a wide variety of vertebrate species. The genetic management of captive populations can be confounded if intrinsic mate choice reduces or eliminates reproductive success between carefully chosen breeding pairs. For example, the San Diego Zoo koala colony only has a 45% copulation rate for matched individuals. Herein, we investigated determinants of koala mating success using breeding records (1984-2010) and genotypes for 52 individuals at four MHC markers. We quantified MHC diversity according to functional amino acids, heterozygosity, and the probability of producing a heterozygous offspring. We then used categorical analysis and logistic regression to investigate both copulation and parturition success. In addition, we also examined age, day length, and average pairwise kinship. Our post-hoc power analysis indicates that at a power level of 1-β = 0.8, we should have been able to detect strong MHC preferences. However, we did not find a significant MHC effect on either copulation or parturition success with one exception: pairs with lower or no production of a joey had significantly lower MHC functional amino acid diversity in the categorical analysis. In contrast, day length and dam age (or age difference of the pair) consistently had an effect on mating success. These findings may be leveraged to improve the success of attempted pairs, conserve resources, and facilitate genetic management.

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.

Transcriptome analyses of immune tissues from three Japanese frogs (genus Rana ) reveals their utility in characterizing major histocompatibility complex class II

Background

In Japan and East Asia, endemic frogs appear to be tolerant or not susceptible to chytridiomycosis, a deadly amphibian disease caused by the chytrid fungus Batrachochytridium dendrobatidis (Bd). Japanese frogs may have evolved mechanisms of immune resistance to pathogens such as Bd. This study characterizes immune genes expressed in various tissues of healthy Japanese Rana frogs.

Results

We generated transcriptome data sets of skin, spleen and blood from three adult Japanese Ranidae frogs (Japanese brown frog Rana japonica, the montane brown frog Rana ornativentris, and Tago’s brown frog Rana tagoi tagoi) as well as whole body of R. japonica and R. ornativentris tadpoles. From this, we identified tissue- and stage-specific differentially expressed genes; in particular, the spleen was most enriched for immune-related genes. A specific immune gene, major histocompatibility complex class IIB (MHC-IIB), was further characterized due to its role in pathogen recognition. We identified a total of 33 MHC-IIB variants from the three focal species (n = 7 individuals each), which displayed evolutionary signatures related to increased MHC variation, including balancing selection. Our supertyping analyses of MHC-IIB variants from Japanese frogs and previously studied frog species identified potential physiochemical properties of MHC-II that may be important for recognizing and binding chytrid-related antigens.

Conclusions

This is one of the first studies to generate transcriptomic resources for Japanese frogs, and contributes to further understanding the immunogenetic factors associated with resistance to infectious diseases in amphibians such as chytridiomycosis. Notably, MHC-IIB supertyping analyses identified unique functional properties of specific MHC-IIB alleles that may partially contribute to Bd resistance, and such properties provide a springboard for future experimental validation.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4404-0) contains supplementary material, which is available to authorized users.

Limited polymorphism of the functional MHC class II B gene in the black-spotted frog (Pelophylax nigromaculatus) identified by locus-specific genotyping

Amphibians can be more vulnerable to environmental changes and diseases than other species because of their complex life cycle and physiological requirements. Therefore, understanding the adaptation of amphibians to environmental changes is crucial for their conservation. Major histocompatibility complex (MHC) presents an excellent tool for the investigation of adaptive variations and the assessment of adaptive potential because it can be under strong diversifying selection. Here, we isolated the MHC class II B (MHCIIB) gene from cDNA sequences of the black-spotted frog (Pelophylax nigromaculatus), a widespread amphibian species in China, and designed locus-specific primers to characterize adaptive variability of this amphibian. Ten alleles were identified from 67 individual frogs of three populations and no more than two alleles were present in each individual animal. Furthermore, none of the sequences had indels or/and stop codons, which is in good agreement with locus-specific amplification of a functional gene. However, we found low polymorphism at both nucleotide and amino acid levels, even in the antigen-binding region. Purifying selection acting at this locus was supported by the findings that the d_N/d_S ratio across all alleles was below 1 and that negatively selected sites were detected by different methods. Allele frequency distributions were significantly different among geographic populations, indicating that physiographic factors may have strong effect on the genetic structure of the black-spotted frog. This study revealed limited polymorphism of three adjacent black-spotted frog populations at the functional MHCIIB locus, which may be attributed to region-specific differences. The locus-specific genotyping technique developed in this study would provide a foundation for future studies on adaptive divergence among different frog populations.

Genetic wealth, population health: Major histocompatibility complex variation in captive and wild ring-tailed lemurs (Lemur catta)

Across species, diversity at the major histocompatibility complex (MHC) is critical to individual disease resistance and, hence, to population health; however, MHC diversity can be reduced in small, fragmented, or isolated populations. Given the need for comparative studies of functional genetic diversity, we investigated whether MHC diversity differs between populations which are open, that is experiencing gene flow, versus populations which are closed, that is isolated from other populations. Using the endangered ring-tailed lemur (Lemur catta) as a model, we compared two populations under long-term study: a relatively "open," wild population (n = 180) derived from Bezà Mahafaly Special Reserve, Madagascar (2003-2013) and a "closed," captive population (n = 121) derived from the Duke Lemur Center (DLC, 1980-2013) and from the Indianapolis and Cincinnati Zoos (2012). For all animals, we assessed MHC-DRB diversity and, across populations, we compared the number of unique MHC-DRB alleles and their distributions. Wild individuals possessed more MHC-DRB alleles than did captive individuals, and overall, the wild population had more unique MHC-DRB alleles that were more evenly distributed than did the captive population. Despite management efforts to maintain or increase genetic diversity in the DLC population, MHC diversity remained static from 1980 to 2010. Since 2010, however, captive-breeding efforts resulted in the MHC diversity of offspring increasing to a level commensurate with that found in wild individuals. Therefore, loss of genetic diversity in lemurs, owing to small founder populations or reduced gene flow, can be mitigated by managed breeding efforts. Quantifying MHC diversity within individuals and between populations is the necessary first step to identifying potential improvements to captive management and conservation plans.

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.

Koala retrovirus genotyping analyses reveal a low prevalence of KoRV-A in Victorian koalas and an association with clinical disease

PURPOSE

Koala retrovirus (KoRV) is undergoing endogenization into the genome of koalas in Australia, providing an opportunity to assess the effect of retrovirus infection on the health of a population. The prevalence of KoRV in north-eastern Australia (Queensland and New South Wales) is 100 %, whereas previous preliminary investigations in south-eastern Australia (Victoria) suggested KoRV is present at a lower prevalence, although the values have varied widely. Here, we describe a large study of free-ranging koalas in Victoria to estimate the prevalence of KoRV and assess the clinical significance of KoRV infection in wild koalas.

METHODOLOGY

Blood or spleen samples from 648 koalas where tested for KoRV provirus, and subsequently genotyped, using PCRs to detect the pol and env genes respectively. Clinical data was also recorded where possible and analysed in comparison to infection status.

RESULTS

The prevalence of KoRV was 24.7 % (160/648). KoRV-A was detected in 141/160 cases, but KoRV-B, a genotype associated with neoplasia in captive koalas, was not detected. The genotype in 19 cases could not be determined. Genomic differences between KoRV in Victoria and type strains may have impacted genotyping. Factors associated with KoRV infection, based on multivariable analysis, were low body condition score, region sampled, and 'wet bottom' (a staining of the fur around the rump associated with chronic urinary incontinence). Koalas with wet bottom were nearly twice as likely to have KoRV provirus detected than those without wet bottom (odds ratio=1.90, 95 % confidence interval 1.21, 2.98).

CONCLUSION

Our findings have important implications for the conservation of this iconic species, particularly regarding translocation potential of Victorian koalas.

Historical gene flow and profound spatial genetic structure among golden pheasant populations suggested by multi-locus analysis

Major histocompatibility complex (MHC) is a good marker system for geographical genetics since they are functional genes in the immune system that are likely to affect the fitness of the individual, and the survival and evolutionary potential of a population in a changing environment. Golden pheasant (Chrysolophus pictus) is a wild Phasianidae distributed in central and north China. In this study, we used a locus-specific genotyping technique for MHC IIB genes of golden pheasant. Combining with microsatellites (simple sequence repeat, SSR) and mitochondrial DNA (mtDNA) D-loop region, we investigated the demographic history and illuminate genetic structure of this bird in detail. SYR (south of Yangtze river) - NYR (north of Yangtze river) lineages, separated by Yangtze River, were defined in genetic structure of MHC IIB. NYR was supposed as refuge during glacial period, suggested by diversity parameters and more ancient alleles in this region. Based on this hypothesis, there was gene flow from NYR to SYR, which was proved by three pieces of evidence: (1) distinct demographic histories of SYR (kept stable) and NYR (experienced expansion); (2) specific affiliation of LC in genetic structure of SSR and MHC genes; (3) significant gene flow from NYR to SYR. Moreover, we also found balancing selection by combination of three Grouping A2's regions (SC, QL and North) into one in Grouping B4 (NYR) and no pattern of isolation by distance (IBD) found in MHC IIB, whereas for SSR we found a relatively strong and significant IBD. Several mechanisms in the evolution of MHC IIB genes, including recombination, historically positive selection, trans-species evolution and concerted evolution, were shown by molecular and phylogenetic analysis. Overall these results suggest the Yangtze River was inferred to be a geological barrier for this avian and NYR might experience population expansion, which invaded into a neighboring region. This study contributes to the understanding of the effects of geographic features on contemporary patterns of genetic variation in the golden pheasant in China, and helps us to define the adaptive unite (AU) for this avian.

Major Histocompatibility Complex Class II in the red-tailed phascogale (Phascogale calura)

Diversity in major histocompatibility complex (MHC) genes can be correlated with the level of immunological fitness of an individual or group of individuals. This study tested published primer sets designed to amplify fragments of the MHC Class II DAB and DBB genes to amplify the equivalent gene fragments in red-tailed phascogales (Phascogale calura). Seventeen genomic DNA samples extracted from phascogale muscle tissue were used to amplify the initial DAB and DBB fragments; however, only DAB PCR proved successful. The fragments were 172 bp in length between the primers and had a high level of identity to other known marsupial MHC Class II DAB gene sequences (89–98%), including those of the koala (Phascolarctos cinereus), Tasmanian devil (Sarcophilus harrisii), common brushtail possum (Trichosurus vulpecula) and several wallaby species. Multiple sequence alignment revealed limited variability of MHC Class II genes between the individuals, but eight individual sequences in total. Genomic DNA was subsequently extracted from three fresh red-tailed phascogale scat samples and DAB fragments successfully amplified. The technique will allow for red-tailed phascogales to be sampled non-invasively in the wild and to determine the level of MHC diversity among individuals in the population.

Antibody and Cytokine Responses of Koalas (Phascolarctos cinereus) Vaccinated with Recombinant Chlamydial Major Outer Membrane Protein (MOMP) with Two Different Adjuvants

Developing a vaccine against Chlamydia is key to combating widespread mortalities and morbidities associated with this infection in koalas (Phascolarctos cinereus). In previous studies, we have shown that two or three doses of a Recombinant Major Outer Membrane Protein (rMOMP) antigen-based vaccine, combined with immune stimulating complex (ISC) adjuvant, results in strong cellular and humoral immune responses in koalas. We have also separately evaluated a single dose vaccine, utilising a tri-adjuvant formula that comprises polyphosphazine based poly I: C and host defense peptides, with the same antigen. This formulation also produced strong cellular and humoral immune responses in captive koalas. In this current study, we directly compared the host immune responses of two sub-groups of wild Chlamydia negative koalas in one population vaccinated with the rMOMP protein antigen and adjuvanted with either the ISC or tri-adjuvant formula. Overall, both adjuvants produced strong Chlamydia-specific cellular (IFN-γ and IL-17A) responses in circulating PBMCs as well as MOMP-specific and functional, in vitro neutralising antibodies. While the immune responses were similar, there were adjuvant-specific immune differences between the two adjuvants, particularly in relation to the specificity of the MOMP epitope antibody responses.

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.

A transcriptome resource for the koala (Phascolarctos cinereus): insights into koala retrovirus transcription and sequence diversity

Background

The koala, Phascolarctos cinereus, is a biologically unique and evolutionarily distinct Australian arboreal marsupial. The goal of this study was to sequence the transcriptome from several tissues of two geographically separate koalas, and to create the first comprehensive catalog of annotated transcripts for this species, enabling detailed analysis of the unique attributes of this threatened native marsupial, including infection by the koala retrovirus.

Results

RNA-Seq data was generated from a range of tissues from one male and one female koala and assembled de novo into transcripts using Velvet-Oases. Transcript abundance in each tissue was estimated. Transcripts were searched for likely protein-coding regions and a non-redundant set of 117,563 putative protein sequences was produced. In similarity searches there were 84,907 (72%) sequences that aligned to at least one sequence in the NCBI nr protein database. The best alignments were to sequences from other marsupials. After applying a reciprocal best hit requirement of koala sequences to those from tammar wallaby, Tasmanian devil and the gray short-tailed opossum, we estimate that our transcriptome dataset represents approximately 15,000 koala genes. The marsupial alignment information was used to look for potential gene duplications and we report evidence for copy number expansion of the alpha amylase gene, and of an aldehyde reductase gene.

Koala retrovirus (KoRV) transcripts were detected in the transcriptomes. These were analysed in detail and the structure of the spliced envelope gene transcript was determined. There was appreciable sequence diversity within KoRV, with 233 sites in the KoRV genome showing small insertions/deletions or single nucleotide polymorphisms. Both koalas had sequences from the KoRV-A subtype, but the male koala transcriptome has, in addition, sequences more closely related to the KoRV-B subtype. This is the first report of a KoRV-B-like sequence in a wild population.

Conclusions

This transcriptomic dataset is a useful resource for molecular genetic studies of the koala, for evolutionary genetic studies of marsupials, for validation and annotation of the koala genome sequence, and for investigation of koala retrovirus. Annotated transcripts can be browsed and queried at http://koalagenome.org.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-786) contains supplementary material, which is available to authorized users.

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.