Future-proofing the koala: Synergising genomic and environmental data for effective species management.

Climatic and evolutionary processes are inextricably linked to conservation. Avoiding extinction in rapidly changing environments often depends upon a species' capacity to adapt in the face of extreme selective pressures. Here, we employed exon capture and high-throughput next-generation sequencing to investigate the mechanisms underlying population structure and adaptive genetic variation in the koala (Phascolarctos cinereus), an iconic Australian marsupial that represents a unique conservation challenge because it is not uniformly threatened across its range. An examination of 250 specimens representing 91 wild source locations revealed that five major genetic clusters currently exist on a continental scale. The initial divergence of these clusters appears to have been concordant with the Mid-Brunhes Transition (~430 to 300 kya), a major climatic reorganisation that increased the amplitude of Pleistocene glacial-interglacial cycles. While signatures of polygenic selection and environmental adaptation were detected, strong evidence for repeated, climate-associated range contractions and demographic bottleneck events suggests that geographically isolated refugia may have played a more significant role in the survival of the koala through the Pleistocene glaciation than in situ adaptation. Consequently, the conservation of genome-wide genetic variation must be aligned with the protection of core koala habitat to increase the resilience of vulnerable populations to accelerating anthropogenic threats. Finally, we propose that the five major genetic clusters identified in this study should be accounted for in future koala conservation efforts (e.g., guiding translocations), as existing management divisions in the states of Queensland and New South Wales do not reflect historic or contemporary population structure.

Phasor histone FLIM-FRET microscopy quantifies spatiotemporal rearrangement of chromatin architecture during the DNA damage response.

To investigate how chromatin architecture is spatiotemporally organized at a double-strand break (DSB) repair locus, we established a biophysical method to quantify chromatin compaction at the nucleosome level during the DNA damage response (DDR). The method is based on phasor image-correlation spectroscopy of histone fluorescence lifetime imaging microscopy (FLIM)-Förster resonance energy transfer (FRET) microscopy data acquired in live cells coexpressing H2B-eGFP and H2B-mCherry. This multiplexed approach generates spatiotemporal maps of nuclear-wide chromatin compaction that, when coupled with laser microirradiation-induced DSBs, quantify the size, stability, and spacing between compact chromatin foci throughout the DDR. Using this technology, we identify that ataxia-telangiectasia mutated (ATM) and RNF8 regulate rapid chromatin decompaction at DSBs and formation of compact chromatin foci surrounding the repair locus. This chromatin architecture serves to demarcate the repair locus from the surrounding nuclear environment and modulate 53BP1 mobility.

First evidence of deviation from Mendelian proportions in a conservation programme.

Classic Mendelian inheritance is the bedrock of population genetics and underpins pedigree-based management of animal populations. However, assumptions of Mendelian inheritance might not be upheld in conservation breeding programmes if early viability selection occurs, even when efforts are made to equalise genetic contributions of breeders. To test this possibility, we investigated deviations from Mendelian proportions in a captive metapopulation of the endangered Tasmanian devil. This marsupial population is ideal for addressing evolutionary questions in conservation due to its large size, range of enclosure types (varying in environmental conditions), good genomic resources (which aid interpretation), and the species' biology. Devil mothers give birth to more offspring than they can nurse in the pouch, providing the potential for intense viability selection amongst embryos. We used data from 140 known sire-dam-offspring triads to isolate within-family selection from population-level mechanisms (such as mate choice or inbreeding), and compared observed offspring genotypes at 123 targeted SNPs to neutral (i.e., Mendelian) expectations. We found lower offspring heterozygosity than expected, and subtle patterns that varied across a gradient of management intensity from zoo-like enclosures to semi-wild environments for some loci. Meiotic drive or maternal-foetal incompatibilities are consistent with our results, although we cannot statistically confirm these mechanisms. We found some evidence that maternal genotype affects annual litter size, suggesting that family-level patterns are driven by differential offspring mortality before birth or during early development. Our results show that deviations from Mendelian inheritance can occur in conservation programmes, despite best-practice management to prevent selection.

From reference genomes to population genomics: comparing three reference-aligned reduced-representation sequencing pipelines in two wildlife species.

BACKGROUND: Recent advances in genomics have greatly increased research opportunities for non-model species. For wildlife, a growing availability of reference genomes means that population genetics is no longer restricted to a small set of anonymous loci. When used in conjunction with a reference genome, reduced-representation sequencing (RRS) provides a cost-effective method for obtaining reliable diversity information for population genetics. Many software tools have been developed to process RRS data, though few studies of non-model species incorporate genome alignment in calling loci. A commonly-used RRS analysis pipeline, Stacks, has this capacity and so it is timely to compare its utility with existing software originally designed for alignment and analysis of whole genome sequencing data. Here we examine population genetic inferences from two species for which reference-aligned reduced-representation data have been collected. Our two study species are a threatened Australian marsupial (Tasmanian devil Sarcophilus harrisii; declining population) and an Arctic-circle migrant bird (pink-footed goose Anser brachyrhynchus; expanding population). Analyses of these data are compared using Stacks versus two widely-used genomics packages, SAMtools and GATK. We also introduce a custom R script to improve the reliability of single nucleotide polymorphism (SNP) calls in all pipelines and conduct population genetic inferences for non-model species with reference genomes. RESULTS: Although we identified orders of magnitude fewer SNPs in our devil dataset than for goose, we found remarkable symmetry between the two species in our assessment of software performance. For both datasets, all three methods were able to delineate population structure, even with varying numbers of loci. For both species, population structure inferences were influenced by the percent of missing data. CONCLUSIONS: For studies of non-model species with a reference genome, we recommend combining Stacks output with further filtering (as included in our R pipeline) for population genetic studies, paying particular attention to potential impact of missing data thresholds. We recognise SAMtools as a viable alternative for researchers more familiar with this software. We caution against the use of GATK in studies with limited computational resources or time.

Development of a SNP-based assay for measuring genetic diversity in the Tasmanian devil insurance population.

BACKGROUND: The Tasmanian devil (Sarcophilus harrisii) has undergone a recent, drastic population decline due to the highly contagious devil facial tumor disease. The tumor is one of only two naturally occurring transmissible cancers and is almost inevitably fatal. In 2006 a disease-free insurance population was established to ensure that the Tasmanian devil is protected from extinction. The insurance program is dependent upon preserving as much wild genetic diversity as possible to maximize the success of subsequent reintroductions to the wild. Accurate genotypic data is vital to the success of the program to ensure that loss of genetic diversity does not occur in captivity. Until recently, microsatellite markers have been used to study devil population genetics, however as genetic diversity is low in the devil and potentially decreasing in the captive population, a more sensitive genotyping assay is required. METHODS: Utilising the devil reference genome and whole genome re-sequencing data, we have identified polymorphic regions for use in a custom genotyping assay. These regions were amplified using PCR and sequenced on the Illumina MiSeq platform to refine a set a markers to genotype the Tasmanian devil insurance population. RESULTS: We have developed a set of single nucleotide polymorphic (SNP) markers, assayed by amplicon sequencing, that provide a high-throughput method for monitoring genetic diversity and assessing familial relationships among devils. To date we have used a total of 267 unique SNPs within both putatively neutral and functional loci to genotype 305 individuals in the Tasmanian devil insurance population. We have used these data to assess genetic diversity in the population as well as resolve the parentage of 21 offspring. CONCLUSIONS: Our molecular data has been incorporated with studbook management practices to provide more accurate pedigree information and to inform breeding recommendations. The assay will continue to be used to monitor the genetic diversity of the insurance population of Tasmanian devils with the aim of reducing inbreeding and maximizing success of reintroductions to the wild.

Lack of genetic diversity across diverse immune genes in an endangered mammal, the Tasmanian devil (Sarcophilus harrisii).

The Tasmanian devil (Sarcophilus harrisii) is threatened with extinction due to the spread of devil facial tumour disease. Polymorphisms in immune genes can provide adaptive potential to resist diseases. Previous studies in diversity at immune loci in wild species have almost exclusively focused on genes of the major histocompatibility complex (MHC); however, these genes only account for a fraction of immune gene diversity. Devils lack diversity at functionally important immunity loci, including MHC and Toll-like receptor genes. Whether there are polymorphisms at devil immune genes outside these two families is unknown. Here, we identify polymorphisms in a wide range of key immune genes, and develop assays to type single nucleotide polymorphisms (SNPs) within a subset of these genes. A total of 167 immune genes were examined, including cytokines, chemokines and natural killer cell receptors. Using genome-level data from ten devils, SNPs within coding regions, introns and 10 kb flanking genes of interest were identified. We found low polymorphism across 167 immune genes examined bioinformatically using whole-genome data. From this data, we developed long amplicon assays to target nine genes. These amplicons were sequenced in 29-220 devils and found to contain 78 SNPs, including eight SNPS within exons. Despite the extreme paucity of genetic diversity within these genes, signatures of balancing selection were exhibited by one chemokine gene, suggesting that remaining diversity may hold adaptive potential. The low functional diversity may leave devils highly vulnerable to infectious disease, and therefore, monitoring and preserving remaining diversity will be critical for the long-term management of this species. Examining genetic variation in diverse immune genes should be a priority for threatened wildlife species. This study can act as a model for broad-scale immunogenetic diversity analysis in threatened species.

The CARD9 Gene in Koalas (Phascolarctos cinereus): Does It Play a Role in the Cryptococcus-Koala Interaction?

Cryptococcus is a genus of fungal pathogens that can infect and cause disease in a range of host species and is particularly prominent in koalas (Phascolarctos cinerus). Like other host species, koalas display a range of outcomes upon exposure to environmental Cryptococcus, from external nasal colonization to asymptomatic invasive infection and, in rare cases, severe clinical disease resulting in death. Host factors contributing to these varied outcomes are poorly understood. Due to their close relationship with eucalypt trees (a key environmental niche for Cryptococcus gattii) and suspected continual exposure to the pathogen, koalas provide a unique opportunity to examine host susceptibility in natural infections. Caspase recruitment domain-containing protein 9 (CARD9) is a key intracellular signaling protein in the fungal innate immune response. Humans with mutations in CARD9 succumb to several different severe and chronic fungal infections. This study is the first to sequence and explore CARD9 variation in multiple koalas using Sanger sequencing. Four CARD9 exons were successfully sequenced in 22 koalas from a New South Wales, Australia population. We found minimal variation between koalas across all four exons, an observation that was also made when CARD9 sequences were compared between koalas and six other species, including humans and mice. Ten single-nucleotide polymorphisms (SNP) were identified in this study and explored in the context of cryptococcal exposure outcomes. While we did not find any significant association with variation in cryptococcal outcomes, we found a high degree of conservation between species at several SNP loci that requires further investigation. The findings from this study lay the groundwork for further investigations of CARD9 and Cryptococcus both in koalas and other species, and highlight several considerations for future studies.

A novel multi-variate immunological approach, reveals immune variation associated with environmental conditions, and co-infection in the koala (Phascolarctos cinereus).

External signs of disease are frequently used as indicators of disease susceptibility. However, immune profiling can be a more effective indicator to understand how host responses to infection may be shaped by host, pathogen and environmental factors. To better inform wildlife health assessment and research directions, we investigated the utility of a novel multivariate immunophenotyping approach examining innate and adaptive immune responses in differing climatic, pathogen co-infection and demographic contexts across two koala (Phascolarctos cinereus) populations in New South Wales: the Liverpool Plains (LP), and Southern Highlands to South-west Sydney (SHSWS). Relative to the comparatively healthy SHSWS, the LP had greater and more variable innate immune gene expression (IL-1ß, IL-6), and KoRV transcription. During extreme heat and drought, koalas from the LP displayed upregulation of a stress pathway gene and reduced adaptive immune genes expression, haematocrit and plasma protein, suggesting the possibility of environmental impacts through multiple pathways. In those koalas, KoRV transcription status, Chlamydia pecorum infection loads, and visible urogenital inflammation were not associated with immune variation, suggesting that immune markers were more sensitive indicators of real-time impacts than observed disease outcomes.

Expanding the known distribution of phascolartid gammaherpesvirus 1 in koalas to populations across Queensland and New South Wales.

Koala populations across the east coast of Australia are under threat of extinction with little known about the presence or distribution of a potential pathogen, phascolartid gammaherpesvirus 1 (PhaHV-1) across these threatened populations. Co-infections with PhaHV-1 and Chlamydia pecorum may be common and there is currently a limited understanding of the impact of these co-infections on koala health. To address these knowledge gaps, archived clinical and field-collected koala samples were examined by quantitative polymerase chain reaction to determine the distribution of PhaHV-1 in previously untested populations across New South Wales and Queensland. We detected PhaHV-1 in all regions surveyed with differences in detection rate between clinical samples from rescued koalas (26%) and field-collected samples from free-living koalas (8%). This may reflect increased viral shedding in koalas that have been admitted into care. We have corroborated previous work indicating greater detection of PhaHV-1 with increasing age in koalas and an association between PhaHV-1 and C. pecorum detection. Our work highlights the need for continued surveillance of PhaHV-1 in koala populations to inform management interventions, and targeted research to understand the pathogenesis of PhaHV-1 and determine the impact of infection and co-infection with C. pecorum.

Real-time analysis of metabolic activity within Lactobacillus acidophilus by phasor fluorescence lifetime imaging microscopy of NADH.

Nicotinamide adenine dinucleotide (NADH) is an endogenous fluorescent molecule commonly used as a metabolic biomarker. Fluorescence lifetime imaging microscopy (FLIM) is a method in which the fluorescence decay is measured at each pixel of an image. While the fluorescence spectrum of free and protein-bound NADH is very similar, free and protein-bound NADH display very different decay profiles. Therefore, FLIM can provide a way to distinguish free/bound NADH at the level of single bacteria within biological samples. The phasor technique is a graphical method to analyse the entire image and to produce a histogram of pixels with different decay profile. In this study, NADH fluorescence decay profiles within Lactobacillus acidophilus samples treated using different protocols indicated discernible variations. Clear distinctions between fluorescence decay profiles of NADH in samples of artificially heightened metabolic activity in comparison to those of samples lacking an accessible carbon source were obtained.

Development of diagnostic and point of care assays for a gammaherpesvirus infecting koalas.

The recent listing of koala populations as endangered across much of their range has highlighted the need for better management interventions. Disease is a key threat to koala populations but currently there is no information across the threatened populations on the distribution or impact of a gammaherpesvirus, phascolarctid gammaherpesvirus 1 (PhaHV-1). PhaHV-1 is known to infect koalas in southern populations which are, at present, not threatened. Current testing for PhaHV-1 involves lengthy laboratory techniques that do not permit quantification of viral load. In order to better understand distribution, prevalence and impacts of PhaHV-1 infections across koala populations, diagnostic and rapid point of care tests are required. We have developed two novel assays, a qPCR assay and an isothermal assay, that will enable researchers, clinicians and wildlife managers to reliably and rapidly test for PhaHV-1 in koalas. The ability to rapidly diagnose and quantify viral load will aid quarantine practices, inform translocation management and guide research into the clinical significance and impacts of PhaHV-1 infection in koalas.

NADH distribution in live progenitor stem cells by phasor-fluorescence lifetime image microscopy.

NADH is a naturally fluorescent metabolite associated with cellular respiration. Exploiting the different fluorescence lifetime of free and bound NADH has the potential to quantify the relative amount of bound and free NADH, enhancing understanding of cellular processes including apoptosis, cancer pathology, and enzyme kinetics. We use the phasor-fluorescence lifetime image microscopy approach to spatially map NADH in both the free and bound forms of live undifferentiated and differentiated myoblast cells. The phasor approach graphically depicts the change in lifetime at a pixel level without the requirement for fitting the decay. Comparison of the spatial distribution of NADH in the nucleus of cells induced to differentiate through serum starvation and undifferentiated cells show differing distributions of bound and free NADH. Undifferentiated cells displayed a short lifetime indicative of free NADH in the nucleus and a longer lifetime attributed to the presence of bound NADH outside of the nucleus. Differentiating cells displayed redistribution of free NADH with decreased relative concentration of free NADH within the nucleus whereas the majority of NADH was found in the cytoplasm.

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.

Characterization of reproductive gene diversity in the endangered Tasmanian devil.

Interindividual variation at genes known to play a role in reproduction may impact reproductive fitness. The Tasmanian devil is an endangered Australian marsupial with low genetic diversity. Recent work has shown concerning declines in productivity in both wild and captive populations over time. Understanding whether functional diversity exists at reproductive genes in the Tasmanian devil is a key first step in identifying genes that may influence productivity. We characterized single nucleotide polymorphisms (SNPs) at 214 genes involved in reproduction in 37 Tasmanian devils. Twenty genes contained nonsynonymous substitutions, with genes involved in embryogenesis, fertilization and hormonal regulation of reproduction displaying greater numbers of nonsynonymous SNPs than synonymous SNPs. Two genes, ADAMTS9 and NANOG, showed putative signatures of balancing selection indicating that natural selection is maintaining diversity at these genes despite the species exhibiting low overall levels of genetic diversity. We will use this information in future to examine the interplay between reproductive gene variation and reproductive fitness in Tasmanian devil populations.

A demonstration of conservation genomics for threatened species management.

As species extinction rates increase, genomics provides a powerful tool to support intensive management of threatened species. We use the Tasmanian devil (Sarcophilus harrisii) to demonstrate how conservation genomics can be implemented in threatened species management. We conducted whole genome sequencing (WGS) of 25 individuals from the captive breeding programme and reduced-representation sequencing (RRS) of 98 founders of the same programme. A subset of the WGS samples was also sequenced by RRS, allowing us to directly compare genome-wide heterozygosity with estimates from RRS data. We found good congruence in interindividual variation and gene-ontology classifications between the two data sets, indicating that our RRS data reflect the genome well. We also attempted genome-wide association studies with both data sets (regarding breeding success), but the genomic data suffered from small sample size, while the RRS data suffered from lack of precision, highlighting a key trade-off in the design of conservation genomic research. Nevertheless, we identified a number of candidate genes that may be associated with variation in breeding success. Individual heterozygosity, as measured by WGS or RRS, was not associated with breeding success in captivity but was negatively associated with litter sizes of breeding females in the RRS data set. Our findings enable conservation managers to have confidence in RRS data while understanding its limitations, and provide avenues for further investigation into which processes underlie variation in breeding success in captive Tasmanian devils. We caution, however, that deep functional insights using RRS may be impaired by a lack of precision, especially when marker density is low.

Too much of a good thing? Finding the most informative genetic data set to answer conservation questions.

Molecular markers are a useful tool allowing conservation and population managers to shed light on genetic processes affecting threatened populations. However, as technological advancements in molecular techniques continue to evolve, conservationists are frequently faced with new genetic markers, each with nuanced variation in their characteristics as well as advantages and disadvantages for informing various questions. We used a well-studied population of Tasmanian devils (Sarcophilus harrisii) from Maria Island, Tasmania, to illustrate the issues associated with combining multiple genetic data sets and to help answer a question posed by many population managers: which data set will provide the most precise and accurate estimates of the population processes we are trying to measure? We analysed individual heterozygosity (as internal relatedness, IR) of 96 individuals, calculated using four genetic marker types (putatively neutral microsatellites, major histocompatibility complex-linked microsatellites, reduced representation sequencing, and candidate region resequencing). We found no correlation in IR values across marker types, suggesting that various genetic markers reflect different aspects of genomic diversity. In addition, some marker types were more informative than others for conservation decision-making. Reduced representation sequencing provided the highest precision (lowest error) for estimating population-level genetic diversity, and most closely reflected genome-wide heterozygosity both theoretically and empirically. Within the conservation context, our results highlight important considerations when choosing a molecular technique for wildlife genetics.

Variants in the host genome may inhibit tumour growth in devil facial tumours: evidence from genome-wide association.

Devil facial tumour disease (DFTD) has decimated wild populations of Tasmanian devils (Sarcophilus harrisii) due to its ability to avoid immune detection and pass from host to host by biting. A small number of devils have been observed to spontaneously recover from the disease which is otherwise fatal. We have sequenced the genomes of these rare cases and compared them to the genomes of devils who succumbed to the disease. Genome-wide association, based on this limited sampling, highlighted two key genomic regions potentially associated with ability to survive DFTD. Following targeted genotyping in additional samples, both of these loci remain significantly different between cases and controls, with the PAX3 locus retaining significance at the 0.001 level, though genome-wide significance was not achieved. We propose that PAX3 may be involved in a regulatory pathway that influences the slowing of tumour growth and may allow more time for an immune response to be mounted in animals with regressed tumours. This provides an intriguing hypothesis for further research and could provide a novel route of treatment for this devastating disease.

New insights into the role of MHC diversity in devil facial tumour disease.

BACKGROUND: Devil facial tumour disease (DFTD) is a fatal contagious cancer that has decimated Tasmanian devil populations. The tumour has spread without invoking immune responses, possibly due to low levels of Major Histocompatibility Complex (MHC) diversity in Tasmanian devils. Animals from a region in north-western Tasmania have lower infection rates than those in the east of the state. This area is a genetic transition zone between sub-populations, with individuals from north-western Tasmania displaying greater diversity than eastern devils at MHC genes, primarily through MHC class I gene copy number variation. Here we test the hypothesis that animals that remain healthy and tumour free show predictable differences at MHC loci compared to animals that develop the disease. METHODOLOGY/PRINCIPAL FINDINGS: We compared MHC class I sequences in 29 healthy and 22 diseased Tasmanian devils from West Pencil Pine, a population in north-western Tasmania exhibiting reduced disease impacts of DFTD. Amplified alleles were assigned to four loci, Saha-UA, Saha-UB, Saha-UC and Saha-UD based on recently obtained genomic sequence data. Copy number variation (caused by a deletion) at Saha-UA was confirmed using a PCR assay. No association between the frequency of this deletion and disease status was identified. All individuals had alleles at Saha-UD, disproving theories of disease susceptibility relating to copy number variation at this locus. Genetic variation between the two sub-groups (healthy and diseased) was also compared using eight MHC-linked microsatellite markers. No significant differences were identified in allele frequency, however differences were noted in the genotype frequencies of two microsatellites located near non-antigen presenting genes within the MHC. CONCLUSIONS/SIGNIFICANCE: We did not find predictable differences in MHC class I copy number variation to account for differences in susceptibility to DFTD. Genotypic data was equivocal but indentified genomic areas for further study.

Phasor-FLIM analysis of NADH distribution and localization in the nucleus of live progenitor myoblast cells.

Analysis of the cellular distributions of coenzymes including NADH may aid in understanding a cells metabolic status. We altered serum concentration (0, 2, and 10%) to induce living myoblast cells to undergo the early stages of differentiation. Through microscopy and phasor-FLIM, we spatially mapped and identified variations in the distribution of free and bound NADH. Undifferentiated cells displayed abundant free NADH within the nucleus along with specific regions of more bound NADH. Complete serum starvation dramatically increased the fraction of bound NADH in the nucleus, indicating heightened requirement for transcriptional processes. In comparison, cells exposed to 2% serum exhibited intermediate free nuclear NADH fraction. Overall our results suggest an order of events in which a cell metabolic status alters significantly during the early stages of serum induced differentiation.

Attitudes toward consumption and conservation of tigers in China.

A heated debate has recently emerged between tiger farmers and conservationists about the potential consequences of lifting the ban on trade in farmed tiger products in China. This debate has caused unfounded speculation about the extent of the potential market for tiger products. To fill this knowledge gap, we surveyed 1880 residents from a total of six Chinese cities to understand Urban Chinese tiger consumption behavior, knowledge of trade issues and attitudes towards tiger conservation. We found that 43% of respondents had consumed some product alleged to contain tiger parts. Within this user-group, 71% said that they preferred wild products over farmed ones. The two predominant products used were tiger bone plasters (38%) and tiger bone wine (6.4%). 88% of respondents knew that it was illegal to buy or sell tiger products, and 93% agreed that a ban in trade of tiger parts was necessary to conserve wild tigers. These results indicate that while Urban Chinese people are generally supportive of tiger conservation, there is a huge residual demand for tiger products that could resurge if the ban on trade in tiger parts is lifted in China. We suspect that the current supply of the market is predominantly met by fakes or substitutes branded as tiger medicines, but not listing tiger as an ingredient. We suggest that the Traditional Chinese Medicine community should consider re-branding these products as bone-healing medicines in order to reduce the residual demand for real tiger parts over the long-term. The lifting of the current ban on trade in farmed tiger parts may cause a surge in demand for wild tiger parts that consumers say are better. Because of the low input costs associated with poaching, wild-sourced parts would consistently undercut the prices of farmed tigers that could easily be laundered on a legal market. We therefore recommend that the Chinese authorities maintain the ban on trade in tiger parts, and work to improve the enforcement of the existing ban.