i
Conservation Genetics and Reproduction in Three Australian Marsupial Species
© J. Gould
Emily J. Miller 2008
A dissertation presented to the University of New South Wales in fulfillment of requirement for the degree of Doctorate of Philosophy in Biological Sciences
ii
“…on [the two largest] these islands are large numbers of cats, which are creatures of miraculous form, as big as a hare; the Head is similar to [that] of a Civet cat, the forepaws are very short, about a finger long. Whereon they have five small Nails, or small fingers, as an ape’s fore-paw, and the two hind legs are at least half an ell long, they run on the flat of the joint of the leg, so that they are not quick in running. The tail is very long, the same as a Meerkat [lemur]; if they are going to eat they sit on their hind legs and take the food with their fore-paws and eat exactly the same as squirrels or apes do.”
- Description of the first sighting of a tammar wallaby in the Abrolhos Islands, Western Australia by Francisco Pelsaert from the translation of Heeres (1899) (originally published 1648)
“DNA…doesn’t lie.” - Mark D. B. Eldridge
For my Family
Image on front cover: Macropus eugenii (Tammar wallaby), reproduced from John Gould, The Mammals of Australia, 1863.
iii
Declaration of Originality
I hereby declare that this submission is my own work and to the best of my knowledge it contains no materials previously published or written by another person, or substantial proportions of material which have been accepted for the award of any other degree or diploma at UNSW or any other educational institution, except where due acknowledgement is made in the thesis. Any contribution made to the research by others, with whom I have worked at UNSW or elsewhere, is explicitly acknowledged in the thesis. I also declare that the intellectual content of this thesis is the product of my own work, except to the extent that assistance from others in the project's design and conception or in style, presentation and linguistic expression is acknowledged.
Signed: ……………………………………… Date: …………………………….. Emily J. Miller
iv
Copyright Statement
I hereby grant the University of New South Wales or its agents the right to archive and to make available my thesis or dissertation in whole or part in the University libraries in all forms of media, now or here after known, subject to the provisions of the Copyright Act 1968. I retain all proprietary rights, such as patent rights. I also retain the right to use in future works (such as articles or books) all or part of this thesis or dissertation.
I also authorise University Microfilms to use the 350 word abstract of my thesis in Dissertation Abstract International (this is applicable to doctoral theses only). I have either used no substantial portions of copyright material in my thesis or I have obtained permission to use copyright material; where permission has not been granted I have applied/will apply for a partial restriction of the digital copy of my thesis or dissertation.
Signed: ……………………………………… Date: …………………………….. Emily J. Miller
v
Authenticity Statement
I certify that the Library deposit digital copy is a direct equivalent of the final officially approved version of my thesis. No emendation of content has occurred and if there are any minor variations in formatting, they are the result of the conversion to digital format.
Signed: ……………………………………… Date: …………………………….. Emily J. Miller
vi
Preface
This thesis consists of five stand-alone papers (Chapters two to six) that are being prepared for submission to international journals of high standing. Each chapter is therefore self-contained and some repetition occurs. To prevent unnecessary duplication a single reference list is provided at the end of the thesis formatted in the style for Conservation Biology.
This thesis is a compilation of my own work, with guidance from my principal supervisors, Catherine Herbert and Mark Eldridge. The contributions of co-authors are detailed below. Chapter 1 2 3 4 5 6 7 Conception EJM EJM, EJM, EJM, EJM, EJM, EJM CAH, CAH, CAH, CAH, CAH, MDBE MDBE MDBE MDBE RJ Sample/data NA EJM, EJM, EJM, EJM, EJM, NA Collection CAH, CAH, CAH, CAH, CAH, NT, NT, MDBE, MDBE, JC, BM, BM, KM, JC, JN, MC, MC, NT, JN MW, HR, NM PO, TF NN, BJ NM, KM Analysis NA EJM EJM EJM EJM EJM NA Writing EJM EJM EJM EJM EJM EJM EJM EJM – Emily J. Miller; CAH – Catherine A. Herbert; MDBE – Mark D. B. Eldridge; RJ – Robert Johnson; NT – Neil Thomas; BM – Brian MacMahon; MC – Martin Clarke; HR – Howard Robinson; NN – Nicole Noakes; NM – Nicola Marlow; KM – Keith Morris; PO – Peter Orell; BJ – Brent Johnson; JC – James Cook; JN – Jan Nedved; MW – Michelle Wilson; TF – Terry Fletcher.
vii
Abstract
Many Australian marsupial species require active population management to ensure their survival in the wild. Such management should be based on a sound understanding of species biology. This thesis explores how knowledge of reproduction and genetics can be applied to the management of three Australian marsupial species faced with contrasting management scenarios.
The ‘vulnerable’ greater bilby is the sole remaining desert bandicoot in Australia. They are a secretive, solitary species whose mating system is unclear. This research examined temporal changes in genetic diversity within two captive breeding programs utilising different management strategies. Using seven microsatellite loci, this study found the regular translocation of new individuals into the population maintained genetic diversity. Parentage analysis revealed the bilby to have a promiscious mating system. Sires and non-sires could not be distinguished by morphological traits.
The tammar wallaby is a polygynous, solitary species that is threatened on mainland Australia, but overabundant on some offshore islands. The population genetics of tammars from the Abrolhos Islands in Western Australia were examined using nine autosomal and four Y-linked microsatellite loci, and mitochondrial DNA. There was a relationship between island size, population size and genetic diversity. The Abrolhos populations have significantly lower genetic diversity and are more inbred than mainland tammars and all sampled populations were significantly differentiated. The Abrolhos and mainland populations should be treated as separate Management Units.
The eastern grey kangaroo is a gregarious, polygynous species that is often locally overabundant. To determine traits influencing male reproductive success, behavioural, morphological, physiological and genetic data were examined and showed dominance status, body size and testosterone concentrations were important factors. Sires were also significantly more heterozygous and genetically dissimilar to females, than non- sires.
As body condition influences individual fitness, and management decisions; five body condition indices (BCI) calculated from morphological data were validated using
viii serum biochemistry and haematology in two kangaroo populations with contrasting body condition. Blood parameters were found to be more reliable indicators of condition, questioning the credibility of BCIs currently used in management.
These studies demonstrate the importance of reproductive and genetic data in assisting wildlife management, regardless of a species conservation status.
ix
Acknowledgements
A career in biology was unexpected and arose from a friend (Joanna Gurung) insisting I attend a biology lecture conducted by Professor David Briscoe, as it was apparently unlike anything else. Eight years later the journey has been incredible and I have experienced so much personal and professional growth, travelled to amazing places and worked with fantastic people. However, none of this would have been possible without the love, support and encouragement from my supervisors, family, friends and colleagues.
I am indebted to my supervisors Dr Cath Herbert, Dr Mark Eldridge and Professor Des Cooper. I am privileged to have had the opportunity to work with you all. I am grateful to Cath for her friendship and giving me the opportunity to do a PhD and sparking my initial curiosity in the biology of marsupials. Cath’s passion, enthusiasm and sheer brilliance have enabled me to learn and experience incredible opportunities. I solely credit Cath for my introduction to the art of ‘fondling testicles’ – all in the name of science of course! Without Cath this PhD would not exist.
Mark Eldridge inspired my passion for genetics, as well as broadening my interests/hobbies to include other taxa. From the very first time we did laboratory work together and discovered sneaky copulations occurring in tammar wallabies, I have never looked back. Mark has been a pillar of support for me during my PhD through continual encouragement and teaching me to appreciate the eccentricities that come with being a biologist. Mark has provided me very useful insights at times of crisis such as “there is something fundamentally wrong with the universe” and “there is nothing worse than a bit of integer overflow”.
Professor Des Cooper has enlightened me over the years learning the intricacies of cricket, and cricket players and possesses inspiring wisdom and knowledge. This PhD involved much fieldwork. James Cook and Jan Nedved dedicated a large amount of their time and energy into capturing kangaroos and wallabies so I could ‘fondle their testicles’, as James puts it. James provided much support, particularly when dealing with spiders. A big thankyou to all the volunteers who have assisted with fieldwork.
x
To Lee Ann Rollins, Joachim Elenz, and the Elenz clan (Mischa, Adam and Klara) who have become dear friends and a second family to me. Thank you all for keeping me sane, particularly Lee who always help me keep things in perspective and company in the lab. The Rollins/Elenz family has treated me as one of their own and words cannot express my gratitude and how much you have all come to mean to me.
To my dear friends, Richard (Dr Dick) Lane, Ryan (no name) Moore, (Ranger) Dan Tilbury, without your love, support and shenanigans I would not have stayed “sane”. I know you will share the same excitement and relief upon my completion. I love you all! A big thanks to Paul Herbert for all the stirring and motivation, the nicknames I will never live down and the endless hours devoted to my fish tank to keep me diverted. Maria Cardoso and Mick Freeman who have keep me grounded and provided a good laugh whenever needed. Tiffanie Nelson provided great support and friendship. I think we were separated at birth in a previous life; Michael Whitehead who has changed me forever, always encouraged my warped sense of humour, and nurtured my nerdiness. To my many other friends that were pivotal in making this thesis happen: Enhua Lee, Alex James, Fiona Thomson, Celine Frere, The entire PBSG Clan, Matt Fahey, Illara Clyde, Kellie Wilson, and Boyd the Barman for the best Guinness poured in Sydney.
A special thanks to the Molecular Ecology and Evolutionary Facility Members – Jennifer Sinclair, Michael Whitehead, Celine Frere, Maria Cardoso, Gianluca Maio, Merel Dalebout, Jackie Chan, Steve Hamilton – and the founder Professor Bill Sherwin. Bill has taught me that the world can be explained by mathematical equations; thank you Clare for letting me play with the MEEF technician; Jen and Clare for the useful manuscript comments; the ‘Samuels Crew’ and all those creatures in the dungeon; Kris Carlyon, Snoop Lothian, Susi Zajitschek, Angela Moles, Russell Bonduriansky, Rob Brooks, Jonathon Russell, Nelika Hughes; Angela Higgins and Jeremy Shearman (Ramaciotti Centre, UNSW) for processing my GeneScan and sequencing; David Warton for statistical advice (Chapter 5).
This research would not be possible without the animals and staff at Waratah Park Earth Sanctuary and Australia Walkabout Wildlife Park; or the many wonderful people
xi with whom we have collaborated with from the Department of Environment and Conservation, Western Australia – Keith Morris, Neil Thomas, Nicole Noakes, Brent Johnson, Peter Orell, Nicola Marlow, Brian MacMahon, Martin Clarke, Howard Robinson; the Taronga Zoo staff and vets, in particular Dr Robert Johnson, Dr Larry Vogelnest and Dr Kimberly Vinette-Herin; Mark Adams (SA Museum) for access to the wild WA bilby samples. Karina Acevedo-Whitehouse for access to the IR macro; Cameron Wood and the staff at the Royal North Shore Hospital; Jay Cox and Captain Crankypants (The Rat Patrol) for sailing us around the Abrolhos Islands.
I would like to acknowledge the ARC Kangaroo Genome Centre of Excellence and the Fertility Management of Koalas and Kangaroos Contraception Program who provided funding for this research. Some funding was also provided by the State Trustees M. A. Ingram Trust and The Royal Zoological Society of New South Wales. All experimental work in the Chapter 4 was approved by the Department of Environment and Conservation (WA) Animal Ethics Committee under the approval numbers 10/2005 and 40/2007. All experimental work in Chapters 5 and 6 was carried out was approved by the University of New South Wales Animal Ethics Committee under the approval number 05/121.
Last but not least, a big thank you to my family. Dad, Mum, Ang, Ben, Michelle, Lachlan, Alison, Nat, Hannah, Summer, Nana, Jon, Andrew and Melinda. You have all constantly been there to listen to me and offer advice. Had it not been for you support, I would not have undertaken such a conquest (and finished!). I love you all very much and am grateful that I am part of the family.
xii
Conference Presentations
Miller, E. J., Eldridge, M. D. B., Morris, K, and Herbert, C. A (2008) Swimming tammars? Relationships amongst the tammar wallaby (Macropus eugenii) populations in the Houtman Abrolhos Archipelago, Western Australia. Oral presentation at the Australian Mammal Society Conference, Darwin, Northern Territory, Australia.
Herbert, C. A., Morris, K., Orell, P. Miller, E. J., Eldridge, M. D. B., and Renfree, M. (2008) Living on the edge: reproductive ecology of tammar wallabies inhabiting the Abrolhos Archipelago, Western Australia. Australian Mammal Society Conference, Darwin, Northern Territory, Australia.
Miller, E. J., Eldridge, M. D. B., and Herbert, C. A. (2007) Morphological, socio-endocrine and genetic traits that influence male reproductive success in the eastern grey kangaroo (Macropus giganteus). Oral presentation at the 6th International Zoo and Wildlife Research Conference of Behaviour, Physiology and Genetics, Berlin, Germany
Miller, E. J., Eldridge, M. D. B., Cooper, D. W. and Herbert, C. A. (2007) Sex in marsupials. Oral Presentation at the Australian Museum Research Forum, Sydney, New South Wales, Australia.
Miller, E. J., Eldridge, M. D. B., Herbert, C. A. et al (2007) Bigger is not better: male reproductive success in a captive breeding program for the greater bilby (Macrotis lagotis) Oral presentation at the Genetics Society for Australia Conference, Sydney, New South Wales, Australia.
Miller, E. J., Eldridge, M. D. B., Herbert, C. A. et al (2006) Traits that influence male reproductive success in the greater bilby (Macrotis lagotis): implications for management. Oral presentation at the Australasian Wildlife Management Society Conference, Auckland, New Zealand.
Miller, E. J., Eldridge, M. D. B., Cooper, D. W. and Herbert, C. A. (2006) Morphological and genetic traits potentially influencing male reproductive success in the eastern grey kangaroo (Macropus giganteus): implications for management. Poster presentation at the Australian Mammal Society Conference, Melbourne, Victoria, Australia.
xiii
Table of Contents
Declaration of Originality ...... iii Copyright Statement ...... iv Authenticity Statement...... v Preface...... vi Abstract ...... vii Acknowledgements...... ix Conference Presentations...... xii Table of Contents...... xiii List of Tables ...... xviii List of Figures ...... xxi List of Plates...... xxiv Chapter 1 General Introduction ...... 25 1.1 Conservation Biology...... 25 1.2 Conservation Genetics...... 27 1.3 Mating Systems and Reproduction...... 30 1.4 Wildlife Management...... 32 1.5 Study Species ...... 36 1.5.1 The Greater Bilby (Macrotis lagotis Reid 1837)...... 36 1.5.2 The Tammar Wallaby (Macropus eugenii Desmarest 1817) ...... 37 1.5.3 The Eastern Grey Kangaroo (Macropus giganteus Shaw 1790) ...... 40 1.6 Study Aims ...... 42 Chapter 2 Bilbies behind bars: the impact of captive management on genetic diversity in a threatened species ...... 45 2.1 Introduction ...... 45 2.2 Materials and Methods ...... 49 2.2.1 Study populations ...... 49 2.2.2 Microsatellite genotyping ...... 51 2.2.3 Genetic analyses ...... 52 2.3 Results ...... 54 2.3.1 Genetic diversity...... 54
xiv
2.3.2 Temporal changes in genetic diversity ...... 55 2.3.3 Genetic diversity in comparison to wild populations ...... 58 2.3.4 Genetic diversity in comparison to studbook estimates...... 59 2.4 Discussion ...... 61 2.4.1 Genetic diversity in the captive bilby populations ...... 61 2.4.2 Comparison to wild populations...... 63 2.4.3 Comparison to studbook estimates ...... 64 2.4.4 Implications for management ...... 65 2.4.5 Conclusions...... 66 Chapter 3 The genetic mating system, male reproductive success and selection on male traits in the Greater Bilby (Macrotis lagotis) ...... 68 3.1 Introduction ...... 68 3.2 Materials and Methods ...... 72 3.2.1 Study population, data and sample collection ...... 72 3.2.2 Microsatellite genotyping ...... 74 3.2.3 Parentage analysis...... 74 3.2.4 Male reproductive success and selection analysis...... 76 3.3 Results ...... 77 3.3.1 Parentage assignment...... 77 3.3.2 Morphological traits and male reproductive success...... 78 3.3.3 Selection analysis...... 80 3.4 Discussion ...... 83 3.4.1 Mating system...... 83 3.4.2 Male reproductive success...... 85 3.4.3 Selection on male traits...... 86 3.4.4 Implications for conservation and management ...... 86 3.4.5 Conclusions...... 87 Chapter 4 Swimming tammars? Genetics of three island populations of tammar wallabies (Macropus eugenii) in the Houtman Abrolhos Acrhipelago, Western Australia ..89 4.1 Introduction ...... 89 4.2 Materials and Methods ...... 96
xv
4.2.1 Study populations ...... 96 4.2.2 Sample collection and DNA extraction ...... 97 4.2.3 Microsatellite amplification and screening...... 98 4.2.4 Mitochondrial DNA amplification and screening ...... 98 4.2.5 Estimates of autosomal and Y-linked microsatellite diversity ...... 99 4.2.6 Estimates of mtDNA haplotypic diversity...... 100 4.2.7 Population structure and gene flow ...... 101 4.2.8 Phylogenetic analysis...... 101 4.2.9 Computer simulations...... 102 4.3 Results ...... 103 4.3.1 Autosomal and Y-linked microsatellite diversity ...... 103 4.3.2 mtDNA diversity ...... 106 4.3.3 Population structure and gene flow ...... 110 4.3.4 Phylogenetics...... 111 4.3.5 Computer simulations...... 113 4.4 Discussion ...... 114 4.4.1 Genetic diversity...... 114 4.4.2 Population structure and gene flow ...... 119 4.4.3 North Island founder and population genetics...... 120 4.4.4 Phylogenetics and conservation significance ...... 121 4.4.5 Conclusions...... 123 Chapter 5 Dominance, body size and internal relatedness influence male reproductive success in eastern grey kangaroos (Macropus giganteus)...... 124 5.1 Introduction ...... 124 5.2 Materials and Methods ...... 129 5.2.1 Study populations ...... 129 5.2.2 Animal capture and handling...... 130 5.2.3 Behavioural dominance ...... 132 5.2.4 Genetic assignment of paternity ...... 132 5.2.5 Paternity analysis ...... 133 5.2.6 Testosterone radioimmunoassay...... 133 5.2.7 Data analysis...... 134
xvi
5.3 Results ...... 136 5.3.1 Male dominance...... 136 5.3.2 Characterisation of microsatellite loci and paternity assignment ...... 139 5.3.3 Traits that influence male reproductive success ...... 139 5.3.4 Relatedness, heterozygosity and reproductive success...... 141 5.4 Discussion ...... 142 5.4.1 Male dominance...... 142 5.4.2 Male reproductive success, dominance, body size and testosterone ...... 143 5.4.3 Male reproductive success and genetic traits...... 144 5.4.4 Implications for management ...... 146 5.4.5 Conclusions...... 147 Chapter 6 Validation of body condition indices using serum biochemistry and haematology in a Macropodid species ...... 148 6.1 Introduction ...... 148 6.2 Materials and Methods ...... 151 6.2.1 Study populations ...... 151 6.2.2 Sample collection...... 152 6.2.3 Biochemical and haematological analyses ...... 154 6.2.4 Kidney fat index ...... 155 6.2.5 Data Analysis...... 156 6.3 Results ...... 157 6.3.1 Serum biochemistry and haematology...... 158 6.3.2 Validation of assumptions for OLS ...... 161 6.3.3 Body condition indices ...... 164 6.3.4 Validation of body condition indices using blood parameters ...... 167 6.4 Discussion ...... 167 6.4.1 Serum biochemistry and haematology...... 168 6.4.2 Body condition indices ...... 170 6.4.3 Relationship between blood parameters and body condition indices...... 170 6.4.4 Conclusions...... 172
xvii
Chapter 7 General Discussion: Marsupial genetics, reproduction and the implications for wildlife management ...... 173 References ...... 182 Appendix 2 ...... 218 Appendix 3 ...... 220 Appendix 4 ...... 223 Appendix 5 ...... 225 Appendix 6 ...... 227 Appendix 7 ...... 228 Appendix 8 ...... 229 Appendix 9 ...... 232 Photo Gallery...... 236
xviii
List of Tables
Table 1-1 Comparison of the biology of the three Australian marsupials species selected for this study, the greater bilby (Macrotis lagotis), tammar wallaby (Macropus eugenii) and eastern grey kangaroo (M. giganteus)...... 44 Table 2-1 Summary of the overall genetic diversity calculated across all years from seven microsatellite loci for two captive bred populations of the greater bilby (Macrotis lagotis), Return to Dryandra (Dryandra) and Peron Captive Breeding Centre (Peron), Western Australia (mean ± se)...... 54
Table 2-2 Pair-wise genetic differentiation (FST) partitioned over time for two captive bred populations of greater bilbies, Macrotis lagotis, (a) Return to Dryandra (Dryandra) and (b) Peron Captive Breeding Centre (Peron), Western Australia...... 58 Table 2-3 Genetic diversity for wild Northern Territory and Queensland populations, in comparison to the captive Western Australian (WA) populations of greater bilbies (Macrotis lagotis)...... 59 Table 2-4 Comparison of the genetic diversity in the captive greater bilby (Macrotis agotis) populations with other wild marsupial populations of varying conservation status...... 64 Table 3-1 Male morphological traits (mean ± se) measured in the greater bilby (Macrotis lagotis) to examine the strength of selection...... 81 Table 3-2 Standardised linear gradients ( ) and matrix of quadratic and correlational selection gradients ( ) in the greater bilby (Macrotis lagotis)...... 81 Table 3-3 The M matrix of eigenvalues from the canonical analysis of presented in Table 3-2...... 82 Table 4-1 Summary of autosomal microsatellite diversity indices for the East Wallabi Island (EWI), West Wallabi Island (WWI) and North Island (NI) tammar wallaby (Macropus eugenii) populations, and compared to a mainland Western Australian population...... 104
Table 4-2 Genetic differentiation (FST) amongst East Wallabi Island (EWI), West Wallabi Island (WWI) and North Island (NI) tammar wallaby (Macropus eugenii) populations, and the mainland Western Australian population...... 104
xix
Table 4-3 Characteristics of the four Y-linked microsatellite loci in the Abrolhos tammar wallaby (Macropus eugenii) populations (East Wallabi Island, EWI n = 20; West Wallabi Island, WWI n = 16, North Island, NI n = 19)...... 105 Table 4-4 Summary of mtDNA control region (642 bp) diversity indices for East Wallabi Island (EWI), West Wallabi Island (WWI) and North Island (NI) tammar wallaby (Macropus eugenii) populations...... 107 Table 5-1 Number of eastern grey kangaroos (Macropus giganteus) observed (sampled and genotyped) in each of the three semi-free ranging populations...... 130 Table 5-2 Summary of the dominance rank (I) for each male, and hierarchy linearity (h) for three semi-free ranging captive populations of eastern grey kangaroos (Macropus giganteus)...... 137 Table 5-3 Summary of the size range of the 10 morphological traits measured in male eastern grey kangaroos (Macropus giganteus) across all individuals and the mean (± se) for alpha and lower ranked males...... 138 Table 5-4 Summary of the mean allelic diversity (AD), range of alleles, mean and
range of observed (Ho) and expected (He) heterozygosity for the three populations of eastern grey kangaroos (Macropus giganteus) used in this study...... 139 Table 6-1 Summary of the body condition indices used to compare two populations of eastern grey kangaroos (Macropus giganteus) representing a population in ‘good’ condition and a population in ‘poor’ condition....154 Table 6-2 Summary of the main causes of fluctuation in the selected serum biochemistry and haematology parameters used in this study. The arrows denote an increase ( ) or decrease ( ) in the circulating levels of analytes ...... 155 Table 6-3 Summary of the mean (± se) and reference range of the serum biochemistry and haematology parameters for two eastern grey kangaroo (Macropus giganteus) populations. Population A = ‘good’ condition, and Population B = ‘poor’ condition in (a) females, and (b) males. ‘Pooled’ represents the mean calculated from pooled population data from A and B
xx
and therefore depicts the mean values across the full spectrum of body conditions in this species...... 159 Table 6-4 Summary of morphological measurements (mean ± se) for females and males in two populations of eastern grey kangaroos (Macropus giganteus) (Population A = ‘good’ condition, and Population B = ‘poor’ condition). * denotes a significant difference between populations...... 161
xxi
List of Figures
Figure 2-1 Distribution of the greater bilby (Macrotis lagotis) in Australia (black = present distribution; mid-grey = historic distribution at European settlement; pale grey = Late-Holocene sub-fossil). The locations of the two captive breeding programs referred to in this study are indicated by black circles...... 48 Figure 2-2 Genetic diversity over time for the greater bilby (Macrotis lagotis), captive breeding programs, Return to Dryandra (Dryandra, ) and Peron Captive Breeding Centre (Peron, ), Western Australia. (a) Mean allelic
diversity (AD) and (b) mean heterozygosity (He)...... 56
Figure 2-3 Inbreeding coefficients (FIS) over time for two captive bred populations of greater bilbies (Macrotis lagotis); Return to Dryandra (Dryandra, ) and Peron Captive Breeding Centre (Peron, ), Western Australia. The dashed arrows indicate when new individuals were translocated into each colony (Dryandra below graph, Peron above graph) and the number represents the number of individuals...... 57 Figure 2-4 Comparison of the studbook (dashed line) and genotypic (solid line) estimates of (a) genetic diversity, and (b) inbreeding coefficients between 1999 and 2007 for the Peron Captive Breeding Centre (Peron) colony, Western Australia...... 60 Figure 3-1 Percent (%) males siring offspring (± se) between 2000 and 2004 in a semi free-ranging captive greater bilby (Macrotis lagotis) population. ..79 Figure 3-2 Relationship between male body weight (g) and the number of offspring sired between 2000 and 2004 in a semi free-ranging captive greater bilby (Macrotis lagotis) population...... 79 Figure 3-3 Relationship between body weight (g) of male greater bilbies (Macrotis lagotis) and the females with which they sired offspring in a semi free- ranging captive population...... 80 Figure 3-4 Visualisation of the fitness surface on the two major axes of nonlinear
selection, m1 and m4. Reproductive success of male greater bilbies (Macrotis lagotis) was standardised to a mean of one...... 82 Figure 4-1 Present and former distribution of the tammar wallaby (Macropus eugenii) (black = present distribution; grey = historic; (A) indicates the
xxii
location the Houtman Abrolhos Archipelago and (B) the Tutanning Nature Reserve, Western Australia...... 93 Figure 4-2 Geographic relationship of East Wallabi, West Wallabi, and North Islands, within the Wallabi Group of the Houtman Abrolhos Archipelago, Western Australia (Hesperian 2007)...... 97 Figure 4-3 Distribution of Y-linked microsatellite haplotypes across the three Abrolhos Island tammar wallaby (Macropus eugenii) populations, East Wallabi Island (black), West Wallabi Island (shaded) and North Island (white)...... 105 Figure 4-4 Distribution of mtDNA control region haplotypes across the three Abrolhos Island tammar wallaby (Macropus eugenii) populations, East Wallabi Island (black), West Wallabi Island (shaded) and North Island (white)...... 106 Figure 4-5 Variable sites within the examined 595 base pair segment of mtDNA control region from 102 tammar wallabies (Macropus eugenii)...... 109 Figure 4-6 Posterior probabilities from the STRUCTURE analysis indicated the data were structured into two clusters. (a) The log probability data, Ln P(D), as a function of K; (b) The rate of change in the log probability of data, K, as a function of K...... 110 Figure 4-7 The level of admixture among tammar wallabies (Macropus eugenii) in the three Wallabi Group Islands. Each vertical bar represents a single individual...... 111 Figure 4-8 Maximum likelihood (ML) tree of relationship amongst Abrolhos, mainland WA and SA tammar wallaby (Macropus eugenii) mtDNA control region sequences (595 bp). Robustness is indicated by bootstrap values ( 50%, above branches) and Bayesian posterior probabilities ( 0.50, below branches)...... 112 Figure 4-9 The predicted rate of loss of genetic diversity (allelic diversity, AD and
heterozygosity, He) for (a) East Wallabi Island and (b) West Wallabi Island populations of tammar wallabies (Macropus eugenii) in the
Abrolhos Archipelago (Ne = 40 ( ); Ne = 60 ( ); Ne = 100 ( ))...... 113 Figure 5-1 Relationship between male dominance rank and body weight (kg) in eastern grey kangaroos (Macropus giganteus) in three semi-free ranging
xxiii
captive populations (Population A ( ) n = 6; Population B ( ) n = 9; and Population C ( ) n = 6). Dominance rank (highest rank = 1 to the lowest rank = 9)...... 138 Figure 5-2 Relationship between male dominance rank and serum testosterone concentration (nmol/L) in three semi-free ranging captive populations of eastern grey kangaroos (Macropus giganteus) (a) Population A, (b) Population B, and (c) Population C. Dominance rank (highest rank = 1 to the lowest rank = 9)...... 140 Figure 5-3 Percent (%) offspring sired by each male eastern grey kangaroo (Macropus giganteus) of varying dominance rank in three semi-free ranging captive populations (Population A (black) n = 6; Population B (shaded) n = 9; Population C (white) n = 6). Dominance rank (highest rank = 1 to the lowest rank = 9)...... 141 Figure 5-4 The mean (± se) difference in individual heterozygosity between sires ( ) and non-sires ( ) for internal relatedness (IR), standardised
heterozygosity (SH) and mean d2 in the eastern grey kangaroo (Macropus giganteus)...... 142 Figure 6-1 Graphs of the ln-transformed data showing the relationship between body weight (kg) and leg length (cm) for (a) female (n = 31) and (b) male (n = 20) eastern grey kangaroos (Macropus giganteus) from two populations...... 162 Figure 6-2 Graphs of the pooled population ln-transformed data showing the relationship between body weight (kg) and pes length (cm) for (a) female (n = 31) and (b) male (n = 20) eastern grey kangaroos (Macropus giganteus) from two populations...... 163 Figure 6-3 Individual differences in the relationship between ln-transformed body weight (kg) and leg length (cm) for (a) females, and (b) males in Population A (‘good’ condition; ) and Population B (‘poor’ condition; )...... 165 Figure 6-4 Individual differences in the relationship between ln-transformed body weight (kg) and pes length (cm) for (a) females, and (b) males in Population A (‘good’ condition; ) and Population B (‘poor’ condition; )...... 166
xxiv
List of Plates
Plate 1-1 The greater bilby (Macrotis lagotis)...... 37 Plate 1-2 A female tammar wallaby (Macropus eugenii) with a pouch young...... 40 Plate 1-3 A female eastern grey kangaroo (Macropus giganteus) with a joey...... 42 Plate 2-1 The Dryandra woodlands, Western Australia, that surrounds the Return to Dryandra captive breeding facility and serves as a reintroduction site for the locally extinct greater bilby (Macrotis lagotis)...... 50 Plate 2-2 Predator proof fencing surrounding the Return to Dryandra captive breeding facility for threatened species in Western Australia...... 51 Plate 3-1 Collecting morphological measurements from a male greater bilby (Macrotis lagotis), (a) head length (mm) and (b) testis length (mm)...... 73 Plate 3-2 Inspecting the pouch of a female greater bilby (Macrotis lagotis) that had twins present...... 73 Plate 4-1 The harsh environment on (a) West Wallabi, (b) East Wallabi, and (c) North Islands in the Houtman Abrolhos Archipelago, Western Australia...... 96 Plate 4-2 Evidence of kinked tails in the Abrolhos Island populations of tammar wallabies, Macropus eugenii (a) categorised subjectively from right to left, a major kink, a minor kink and a straight tail, (b) close up of a tail with a more prominent kink...... 116 Plate 4-3 Photographs of sexually mature male tammar wallabies (Macropus eugenii) exhibiting (a) unilateral cryptorchism on North Island, Western Australia, and (b) a male with normal testicles from Kangaroo Island, South Australia...... 117 Plate 5-1 Male eastern grey kangaroos (Macropus giganteus)...... 129 Plate 5-2 Collecting a range of morphological measurements (cm) from male eastern grey kangaroos (Macropus giganteus). (a) head length, (b) forearm length, (c) leg length, (d) pes length, (e) tail length, and (f) testis size...... 131 Plate 6-1 Collecting a blood sample from the lateral tail vein of a male eastern grey kangaroo (Macropus giganteus)...... 154
Chapter 1 General introduction 25
Chapter 1 General Introduction
Many wild animal populations currently require some form of active management to ensure their survival. Human modifications to the environment, including habitat destruction and fragmentation, modification of faunal assemblages through the removal or introduction of predators and competitors, modification of floral communities (e.g. conversion of shrubland to pasture), alteration to fire regimes and water sources, and the process of urbanisation have resulted in dramatic changes to the distribution and abundance of endemic species (Caughley 1987). Since the onset of European settlement in Australia, marsupial populations have fluctuated greatly. In most cases there has been a substantial reduction in the abundance and range of small to medium sized mammals, with many becoming extinct or being eliminated from the mainland (Short & Smith 1994). But, some marsupial populations have benefited and are now considered overabundant at a local or regional scale (Herbert 2007; Kinnear et al. 2002). In some cases these species are secure throughout their range, for example the five large kangaroo species (DEH 2007), while in other cases there are geographic and temporal variations in the species conservation status, for example koalas and some small wallaby species (Herbert 2007). Effective management of species at both ends of this continuum, i.e. from critically endangered to overabundant, requires a thorough understanding of many aspects of species biology.
1.1 Conservation Biology
The goal of conservation biology is to conserve global biodiversity. Conservation biology has been described as a ‘crisis discipline’, and encompasses a broad range of fields including ecology, demography, genetics, wildlife and resource biology (Soulé 1985). The current losses of biodiversity and rates of extinctions are unprecedented (Myers & Knoll 2001) and the role of conservation biology has become crucial in the 21st Century. Of the 5 487 mammals species worldwide, 1 141 (21%) are threatened and since 1600, 76 species of mammal have gone extinct worldwide (IUCN 2008). At
Chapter 1 General introduction 26 least half of these extinctions have occurred in Australia in the last 200 years following European settlement of the continent (IUCN 2008; Short & Smith 1994), and approximately 22% of all Australian species are now in decline (Ceballos & Ehrlich 2002). Deterministic and stochastic threats have been identified as the two main types of threats contributing to extinctions (Caughley 1994). Deterministic threats include habitat loss, introduced species, overexploitation, pollution, disease and climate change (Vié et al. 2008). Stochastic threats are random changes in demographic, environmental, catastrophic and genetic factors (Frankham et al. 2002; IUCN 2006). For many species it is still difficult to quantify threatening processes as their biology is poorly known.
There is evidence that conservation efforts can bring species back from the brink of extinction (Brooke et al. 2008). The conservation status of the black-footed ferret (Mustela nigripes) changed from ‘Extinct in the Wild’ to ‘Endangered’ after a successful reintroduction by the United States Fish and Wildlife Service into eight western states and Mexico between 1991 and 2008 (Vié et al. 2008). Similarly, after successful reintroductions in Mongolia since the early 1990s, the Prezwalski’s horse (Equus prezwalskii) has been downgraded from ‘Extinct in the Wild’ in 1996 to ‘Critically Endangered’ (Vié et al. 2008). Within Australia, several species have had their conservation status downgraded due to intensive conservation efforts including the western quoll (Dasyurus geoffroii), Lumholtz’s tree kangaroo (Dendrolagus lumholtzi), tammar wallaby (Macropus eugenii), western brush wallaby (M. irma), burrowing bettong (Bettongia lesueur), Shark-Bay mouse (Pseudomys fieldi) and the greater stick-nest rat (Leporillus conditor) (IUCN 2008).
Threatened species are typically the focus of conservation efforts and taxa that are common receive little attention (Gaston & Fuller 2007). Common species are fundamental to ecosystem structure and function (Hobbs & Mooney 1998) and since their distribution is typically widespread, their impacts often cover large geographical areas (Gaston 2000). Modifications to the landscape have led to many common species suffering considerable reductions in population size (Hobbs & Mooney 1998), but without the threat of impending extinction (Gaston & Fuller 2007). Other common species have adapted to and exploited the modifications, increasing their distribution
Chapter 1 General introduction 27 and/or abundance within their range (Calaby & Grigg 1989; Caughley 1987). Overabundant or expanding populations can have a number of negative impacts on other native species that are rarer and/or less adaptable (Garrott et al. 1993), including reducing natural diversity by monopolising resources, introducing or spreading infectious diseases and parasites, changing the species composition or relative abundance of sympatric species, and even causing local extinctions (Garrott et al. 1993; Noss 1990; Temple 1990). Therefore managing overabundant species is important for the conservation of threatened species.
1.2 Conservation Genetics
Genetics has become an intrinsic component of conservation biology (Frankel 1974), and an instrumental tool in wildlife management and conservation (Hedrick 2001) for both wild and captive populations. Genetics was applied more generally to conservation in the early 1980s after the establishment of the genetic principles for conservation biology (Frankel & Soulé 1981; Schonewald-Cox et al. 1983; Soulé & Wilcox 1980). Subsequently conservation genetics emerged as a prominent discipline and is defined as
“the application of genetics to preserve species as dynamic entities capable of coping with environmental change. It encompasses genetic management of small populations, resolution of taxonomic uncertainties, defining management units within species and the use of molecular analyses in forensics and understanding species’ biology” (Frankham et al. 2002, pp 1).
Techniques in molecular biology have advanced rapidly during the past few decades (Sunnucks 2000). The first method for measuring genetic diversity was discovered in 1966 using protein electrophoresis (Harris 1966; Lewontin & Hubby 1966). Subsequently, several techniques have been devised for examining genetic diversity, such as DNA fingerprinting (Jeffreys et al. 1985a, b), microsatellites (Jarne & Logoda 1996), and single nucleotide polymorphisms (Brooks 2003). The development of the polymerase chain reaction (PCR) revolutionised genetics as it allowed the
Chapter 1 General introduction 28 amplification of small amounts of template DNA, broadening the scope of questions able to be addressed in conservation biology (Bartlett & Stirling 2003). Continual advances in technology and the development of powerful analysis methods (Luikart & England 1999) has led to the sequencing of mitochondrial and nuclear genes, as well as the recent advent of sequencing entire genomes (Sunnucks 2000). These advances are allowing scientists to tap into data that would be otherwise unobtainable (Sunnucks 2000).
The IUCN Red List has recognised that the threat status of 15% of mammalian taxa could not be determined as they are ‘data deficient’ (IUCN 2008). Conservation genetics has been widely applied in wildlife management as it provides unique insights into species biology and makes available information for prioritising conservation efforts. Conservation genetics is applicable at three levels: that of the species, population and the individual (Frankham 1995a), and all levels are important to conserve (Frankham et al. 2002) as they are fundamental units in ecology and evolution. On a species level, genetics enables the clarification and identification of taxonomic uncertainties, which is fundamental to conserving a species. For example, until recently, only one species of tuatara (Sphenodon genus) was recognised in New Zealand, which represented the sole extant genus of the entire reptilian Order. Molecular analyses have revealed this species should be divided into three distinct taxonomic groups and has had significant implications for conservation planning (Daugherty et al. 1990). The increased sensitivity of technology and methodology has facilitated the detection of introgression (Mallet 2005; Rhymer & Simberloff 1996), which has conservation implications through the loss of the genetic integrity of a species, for example recovery efforts for the highly threatened Ethiopian wolf (Canis simiensis) have struggled to find pure bred animals for captive breeding due to introgression with domestic dogs (Gottelli et al. 1994). Phylogeography has also become an increasingly valuable technique for the identification of a species origin, patterns of geographic dispersal and subsequent colonisation(s) (Edwards et al. 2007). For example, Y-linked microsatellite haplotypes have been used to elucidate the global dispersal patterns of human populations (Kayser et al. 2001).
Chapter 1 General introduction 29
The determination of genetic relationships on a population level is essential for obtaining information on the genetic distinctness of populations, distinct lineages, identification of separate management units, populations of high conservation priority, delineating species mating systems (Frankham et al. 2002; Hughes 1998), as well as the demographic and historical structure of a population (Escudero et al. 2003). For example, the origin of Ireland’s fauna has been the subject of much debate but with the use of mitochondrial DNA, Finnegan et al. (2008) were able to confirm the Irish red squirrel (Sciurus vulgaris) population arose from a colonisation event on the island, as well as delineate regional genetic structure of the species and identify management units for conservation. Genetics is also important for obtaining information on the impact of genetic stochasticity and its effect on random genetic drift, inbreeding, changes in allele frequencies and genetic diversity (Frankham 1995a; Shaffer 1981). Such information is required for the genetic management of species, particularly threatened, captive bred populations to avoid any further genetic deterioration in addition to what may have already occurred in the declining wild population.
The identification of ‘units’ for conservation has become increasingly popular in conservation biology, but also hotly debated (Crandall et al. 2000; Fraser & Bernatchez 2001; Moritz 1994; Moritz 1999; Ryder 1986; Waples 1991), however all aim to conserve the adaptive genetic variance within species (Fraser & Bernatchez 2001). Two prominent concepts for units of conservation are Evolutionary Significant Units (ESUs) and Management Units (MUs) within ESUs (Hedrick et al. 2001c). An ESU is a group of individuals with a similar genetic composition that demonstrate deep evolutionary divergence from other groups within the same species (Frankham et al. 2002; Moritz 1994). The original definitions of ESUs incorporated two main elements: reproductive isolation (i.e. genetic distinctness) and ecological distinctness (Ryder 1986; Waples 1991). More recently, Moritz (1994) shifted the emphasis to genetic distinctness, suggesting ESUs should be recognised as ‘reciprocally monophyletic for mtDNA alleles and show significant divergence of allele frequencies at nuclear loci’ (Moritz 1994). MUs focus on identifying interactions between closely related populations within a species (Moritz 1994; Moritz 1999) and represent the ecological component of ESUs. Populations that have diverged in allele frequency but do not yet
Chapter 1 General introduction 30 show reciprocal monophyly for mtDNA are also significant for conservation, and are referred to as MUs (Moritz 1994; Moritz 1999).
Much research focuses on the individual level, examining individual fitness traits, identifying source populations for individuals, paternity assignment, reproductive success, and mate choice (Charmantier & Sheldon 2006; Mays & Hill 2004). Identifying source populations of individuals is essential for their translocation or reintroduction into other populations to avoid inbreeding and introgression (Frankham 1995a; Frankham et al. 2002). Accurate parentage data is essential for wildlife managers when making decisions regarding the pairing of individuals for mating and the transfer of individuals among wild and captive populations. The widespread application of genetics has also revealed inconsistencies in the presumed relationship between social organisation, parentage and mating system (Ambs et al. 1999; Coltman et al. 1999a; Issac 2005; Worthington Wilmer et al. 1999), transforming how mating systems are understood.
1.3 Mating Systems and Reproduction
Crucial for effective management is a knowledge of a species mating system and parentage to determine effective population sizes, detect inbreeding, genetic differentiation, and genetic diversity, all of which affect population growth, survival and evolutionary potential (Ralls & Ballou 1986; Reed & Frankham 2003). Parentage cannot be determined using behavioural data alone as genetic analyses have revealed inconsistencies in the presumed relationship based on observational data (Ambs et al. 1999; Coltman et al. 1999a; Worthington Wilmer et al. 1999). Mating systems describe the social and genetic relationships of individuals in relation to reproduction (Emlen & Oring 1977). The diversity in mammalian mating systems is a product of the reproductive strategies of individuals and their environments (Clutton-Brock 1989; Emlen & Oring 1977; Reynolds 1996). The classic mating system classifications are monogamy, polygyny, polyandry and promiscuity, though the precise definition of each can differ even within the field of evolutionary biology (Andersson 1994).
Chapter 1 General introduction 31
Typically, monogamy is defined as the continuing bond and exclusive relationship between a male and a female. Polygyny occurs when a single male mates with several females, and polyandry is the reverse, whereby a single female mates with several males (Wittenberger 1981). In a promiscuous mating system, both males and females mate with more than one partner and there is no long term relationship (Wittenberger 1979). Polygyny commonly occurs in species where males do not provide parental care (i.e. most mammals). As a result, males are able to invest their energy into competing for resources and mates (Emlen & Oring 1977). The ability to monopolise resources is a key factor in the intensity of sexual selection. Polygyny is predicted to occur when there is spatial clumping of females in such a way that enables males to defend multiple females from being accessed by other males (Clutton-Brock 1989; Emlen & Oring 1977) and is characterised by greater variance in male reproductive success than in females (Hoogland & Foltz 1982). Polyandrous mating systems are characterised by females pairing with multiple males during a single breeding season, however, true polyandry has been documented in very few species (Wittenberger 1979). In contrast, promiscuity generally occurs when males are unable to successfully monopolise access to females. It often occurs in species where female groups are unstable, males do not provide parental care, and male home ranges are distributed throughout the home ranges of several females (Clutton-Brock 1989). Males are often unable to defend territories and therefore cannot monopolise access to females (Ramsay & Stirling 1986).
According to sexual selection theory (Darwin 1859), males ought to compete for access to receptive females, the most limiting resource for male reproductive success (Trivers 1972). There is evidence in several taxa demonstrating that a larger body size is advantageous when competing for access to receptive females (Clinchy et al. 2004; Fisher & Cockburn 2005; Fisher & Lara 1999). The male-biased sexual size dimorphism that exists in many mammalian taxa is often attributed to competition between individuals for reproductive opportunities, encouraging the evolution of secondary sexual traits (Andersson 1994; Birkhead 2000). Female preference for particular traits could potentially increase the strength of selection on male traits. For example, if females prefer larger males because their size indicates quality, strength, or
Chapter 1 General introduction 32 fighting ability (Andersson 1994; Bercovitch et al. 2003; Coltman et al. 2001; Schulte- Hostedde & Millar 2002), then we would expect strong selection for body size.
Mating system and parentage studies are important for understanding ecological and evolutionary forces influencing individual fitness, and thus reproductive strategies, and are strongly linked with effective population size (Ne). Ne is defined as:
“the number of individuals that would give rise to the calculated inbreeding coefficient, loss of heterozygosity or variance in allele frequency if they behaved in the manner of an idealised population” (Frankham et al. 2002, pp 189).
In mating systems where there is high reproductive skew and some individuals are excluded from breeding (Nunney 1991, 1993), Ne can become very low (Frankham 1995b; Nunney & Elam 1994). In populations of conservation concern, this is of particular importance as Ne determines the rate of loss of genetic diversity and inbreeding which contributes to the risk of extinction (Frankham 1995b; Frankham et al. 2002).
1.4 Wildlife Management
Many wildlife populations need to be actively managed to stem the loss of global biodiversity and extinctions. Due to limited resources available for conservation and management, the application of genetics is an important tool for prioritising conservation efforts. Effective conservation requires the identification of key threatening processes contributing to a species decline, and the implementation of management practices to modify their effects (Caughley & Gunn 1996). The action taken depends on the processes involved and there is often no single solution. Management actions may include increasing the suitability and area of habitat, captive breeding, reintroductions and/or translocations to artificially improve gene flow and re- establishment of populations, cross-fostering and predator control (Lindenmayer & Burgman 2005). Conservation genetics can contribute to each step of the decision-
Chapter 1 General introduction 33 making process through resolving taxonomic uncertainties so that managers are confident of the status of populations, as well as identifying the relationships among the populations they aim to conserve. Decisions regarding the future of wild and captive populations also rely on knowledge of genetic diversity, levels of inbreeding, and the degree of differentiation among populations (Frankham et al. 2002).
Usually the causes of species decline are multifaceted therefore management strategies need to take a similar approach. For example in Western Australia (WA), predation by the introduced European red fox (Vulpes vulpes) has contributed to the decline of many species. In response, the largest wildlife recovery program ever undertaken in Australia, the Western Shield Program, was established in 1996 by Western Australia’s Department of Environment and Conservation (DEC). The main aims of this program are to maximise the recovery of sustainable populations of threatened species through predator control, and reintroduce threatened species into the former habitats where they have become locally extinct (CALM 1999). The program continues to monitor the recovery of these populations, increase awareness though education and public relations, and makes the best use of new and existing research to enhance its recovery programs (CALM 1999). Captive breeding programs were established for 14 threatened species in WA to increase the number of individuals available for wild reintroductions (Mawson 2004a). In response to predator control, population numbers for many species have increased and species have expanded beyond the areas undergoing predator control (Kinnear et al. 2002), which has been instrumental in the delisting of several species. For example the status of the brush-tailed bettong (B. penicillata) has been changed from ‘Endangered’ to ‘Lower Risk’, and both the western quoll (D. geoffroii) and numbat (Myrmecobius fasciatus) have been de-listed from ‘Endangered’ to ‘Vulnerable’ (Orell 2004). Many species have also been reintroduced into areas where they had become locally extinct (Mawson 2004a).
In contrast to the broad pattern of mammal declines across Australia, some species have experienced an increase in abundance and/or have expanded their range in response to human-induced alterations to the landscape (Hobbs & Mooney 1998). ‘Overabundance’ is a value-laden term, which is usually applied to a population when they exceed the carrying capacity of the environment, have undesired effects on the
Chapter 1 General introduction 34 ecosystem and/or other species, or interfere with humans and livestock (Caughley 1987; Herbert 2007). Animals at such high densities pose a risk to their own welfare through a reduction in available resources which leads to reductions in fecundity, disease resistance, and survival. Often welfare considerations play an important role in management decisions, for example culling animals deemed to be in poor condition (DEH 2007). There are several options for managing overabundant populations including doing nothing, translocation, commercial and non-commercial harvesting, restricting dispersal though fencing, reducing fertility (fertility control) through contraception or sterilisation, and euthanasia (Caughley 1987; Coulson 1998; Herbert 2007). The appropriate management strategy employed is usually governed by socio- economic and political factors (Coulson 1998), particularly when these populations are located in close proximity to nature reserves and urban areas (Adderton Herbert 2004), or involve iconic species, for example the koala (Herbert 2007).
The potential evolutionary consequences of harvesting wildlife have received little attention until recently. There is now a growing body of evidence suggesting that practices such as size-selective harvesting (or trophy hunting) can have negative ecological and evolutionary consequences at a population level, for example changes in selective pressures on morphological traits (Allendorf et al. 2008). As males are often larger than females, size-selective harvesting is usually sex-biased towards males. Coltman et al. (2003) have shown that size-selective harvesting in big-horn rams (Ovis canadensis) has lead to a significant reduction in body weight and horn size of males in one population that has been intensively studied for 30 years. In addition to potential phenotypic changes in males, sex-specific culling can influence population dynamics by reducing female fecundity, disrupting gene flow between populations and altering fine scale genetic structure (Allendorf et al. 2008; Harris et al. 2002). There is debate about the extent to which male sex-biased harvesting could limit population size by reducing female fecundity when trophy males are removed (Fairall 1985), with some scientists suggesting that this is unlikely if young males can fertilise females in the absence of trophy males (Mysterud et al. 2002). But, sex-biased harvesting is thought to have led to reproductive collapse in some species, for example several ungulates (Fairall 1985; Milner-Gulland et al. 2003; Mysterud et al. 2002; Solberg et al. 2002) and elephants (Milner-Gulland & Mace 1991). In a long-term (24 year) study
Chapter 1 General introduction 35 of a red deer (Cervus elaphus) population on the Isle of Rum, removal of hunting pressure led to an increase in the number of breeding females, which was associated with a decrease in the levels of polygyny and a decline in female population genetic structure (Nussey et al. 2006). These examples show the extent to which hunting (or the removal of the pressures exerted by hunting) can influence population genetic structure with long term evolutionary consequences. Unfortunately there are very few long-term studies on harvested mammalian populations, so there is very little empirical data on the consequences of selective harvesting. The first step towards developing a greater understanding of the genetic consequences of harvesting is to collect data on the population genetic substructure, mating system and the extent to which targeted traits are heritable in the target populations (Allendorf et al. 2008).
Fertility control is becoming increasingly attractive to control overabundant mammals, particularly highly valued native mammals (Garrott et al. 1993), as it is advocated as a more humane, ethical method than culling. Methods of fertility control include surgical sterilisation, steroidal and non-steroidal contraception, immunocontraception and dopamine agonists that inhibit lactation. The advantages and disadvantages of each methods are reviewed elsewhere (Adderton Herbert 2004; Cooper & Herbert 2001) and will not be dealt with in this thesis. The management of overabundant native species requires a balance between overcoming the negative impacts of the species, yet maintaining an adequate population base with sufficient genetic diversity to ensure species survival within the area. From a genetic perspective, the use of a reversible fertility control is beneficial as it would enable genetic diversity to be maintained in populations (Adderton Herbert 2004; Herbert 2007). This may be an especially important issue for some of the smaller macropodid species that are threatened on a national scale, but have temporal or localised population eruptions such as the black footed rock wallabies (Petrogale lateralis) in WA and koalas (P. cinereus) in Victoria (Herbert 2007). A recent study by Tanaka et al. (submitted) has concluded that regular contraception, in which all females are allowed to reproduce at some stage during their life, will not affect effective population size (Ne in the genetic sense).
Chapter 1 General introduction 36
1.5 Study Species
In order to explore these issues in the context of wildlife management, reproductive and genetic data was collected and applied to three Australian marsupial species: the greater bilby (Macrotis lagotis), tammar wallaby (Macropus eugenii) and the eastern grey kangaroo (M. giganteus). Each species has a different life history strategy and has been affected by anthropogenic changes to the landscape in varying ways, resulting in contrasting current management objectives. The differences between the three study species examined in this thesis are summarised in Table 1-1.
1.5.1 The Greater Bilby (Macrotis lagotis Reid 1837)
The greater bilby is an iconic species for conservation and the sole remaining species of desert bandicoot in Australia (Johnson 2002). Bilbies are distinguished by their soft, silky hair, long rabbit-like ears, and long furry black and white tail (Plate 1-1). They are an omnivorous, nocturnal marsupial that live solitarily or in pairs (Johnson 2002). Bilbies dig burrows up to two metres deep, and individuals may have up to a dozen burrows they use within their home range (Moseby & O'Donnell 2003). Bilbies exhibit sexual dimorphism with males being, on average, one to two times larger than females (males: 800 – 2500g; females: 600 – 1100g). They are capable of breeding all year round, producing one to three offspring per litter, and up to four litters per year under ideal conditions. Young permanently exit the pouch at approximately 75 days, before becoming independent at 90 days of age, and males do not provide parental care (Southgate et al. 2000). The bilby mating system is not well documented due to their secretive nature and the difficulties associated with wild observations. The mating system of the bilby is thought to be polygynous or promiscuous (Johnson & Johnson 1983; Lee & Cockburn 1985; Moritz et al. 1997).
Bilbies have experienced an 80% reduction in their range throughout the arid and semi-arid zones of Australia (Southgate 1990) and are now confined to north-western Central Australia (from the Tanami Desert, south to Warburton, WA,) with some isolated satellite populations in Queensland (QLD) (Southgate 1990). In addition to habitat loss, bilbies have suffered a significant decline in numbers due to competition with introduced species and predation (Moritz et al. 1997), and the species is listed as
Chapter 1 General introduction 37
‘Vulnerable’ nationally (IUCN 2006). As a conservation response, a National Recovery Plan has been established (Pavey 2006) that aims to implement predator control for European red foxes (V. vulpes) and feral cats (Felis catus) in areas where bilbies are in decline, maintain husbandry and captive breeding programs with a coordinated management of populations; reintroduce captive-bred animals to areas of their former range and continue refinement of methodology for monitoring wild populations (Pavey 2006).
© E. Miller
Plate 1-1 The greater bilby (Macrotis lagotis).
1.5.2 The Tammar Wallaby (Macropus eugenii Desmarest 1817)
The tammar wallaby is one of the smallest Macropus species, with a dark grey-brown upper coat, a pale buff-grey coat beneath, and reddish arms, feet and flanks (Plate 1-2). Most individuals also have a faint white cheek stripe (Smith & Hinds 2002). Tammar wallabies rest during the day in low, dense vegetation and move into open grassy areas for feeding during dusk and the evening (Smith & Hinds 2002). They are a solitary species whose home range overlaps with several other individuals while they feed in the same area, but no territorial behaviour has been documented (Smith & Hinds
Chapter 1 General introduction 38
2002). The tammar wallaby is sexually dimorphic with males being one to two times larger than females (males: 6 – 10kg; females: 4 – 6kg). Tammar wallabies have a polygynous mating system, in which females usually mate with more than one male (Miller et al. in press; Rudd 1994; Smith & Hinds 2002). Breeding is highly seasonal and synchronised, with a large proportion of births occurring in late January and early February (Tyndale-Biscoe & Renfree 1987). Females exhibit a post-partum oestrus within an hour or two of giving birth and all mating activity usually concludes within four hours of parturition (Renfree et al. 1989; Rudd 1994). The resultant conceptus grows to a 100-cell blastocyst before entering embryonic diapause (Tyndale-Biscoe & Renfree 1987). The loss or removal of a pouch young during the breeding season (the period of decreasing day length from the summer solstice to the winter solstice) stimulates the quiescent blastocyst to reactivate. Birth and a new post-partum oestrus then occurs approximately 26.5 ± 0.4 days later (Renfree et al. 1989). Loss or removal of a young during the non-breeding season (the period of increasing day length between the winter and summer solstice) does not result in reactivation of the quiescent blastocyst (Tyndale-Biscoe & Renfree 1987).
Tammar wallabies have a fragmented distribution with populations inhabiting mainland south-western Australia, Garden Island, two islands in the Recherche Archipelago, and three islands in the Abrolhos Archipelago (WA) and Kangaroo Island, South Australia (SA). The SA mainland, Flinders Island and St Peters Island populations are now listed as ‘Extinct in the Wild’ (Poole et al. 1991). Two subspecies of tammar wallaby are commonly recognised in Australia, M. e. derbianus (WA) and M. e. decres (Kangaroo Island, SA), and despite a growing body of evidence indicating morphological and genetic distinctions, they are generally grouped together and referred to as M. eugenii (DEH 2004). Declines in the tammar wallaby have been attributed to land clearing for agriculture, predation by introduced species and stochastic environmental events such as bushfires (DEH 2004).
The tammar wallaby poses interesting issues for wildlife management. They have suffered a dramatic decline in numbers on the mainland and some islands, but thrive on several other offshore islands. Tammar wallabies are listed as ‘Low Risk (near
Chapter 1 General introduction 39 threatened)’ on a national scale (IUCN 2006), however there is a marked difference between the mainland and island populations. Tammar wallabies from New Zealand that are believed to be descendents from the original mainland SA population (Taylor & Cooper 1999) are currently being bred in captivity for reintroduction to mainland SA (DEH 2004; Taylor & Cooper 1999). On mainland WA, the tammar wallaby was de-listed from the State Threatened Fauna List as a result of their recovery under the Western Shield fox baiting program (Kinnear et al. 2002; Possingham et al. 2003). They are now in sufficiently high densities at a localised scale in the wheat belt of south-western WA that these populations are being used as a source for translocations into other areas of their former range (Mawson 2004b). On several offshore islands they are considered overabundant, for example on Kangaroo Island (SA) tammar wallabies are at such high densities, thousands are culled annually to reduce competition with domestic stock for pasture (Wright & Stott 1999). Despite the tammar wallaby being a model species for a wide range of marsupial research (Hinds et al. 1990; Smith & Hinds 2002; Tyndale-Biscoe & Renfree 1987), surprisingly little is known about their evolutionary ecology.
Chapter 1 General introduction 40
© E. Miller
Plate 1-2 A female tammar wallaby (Macropus eugenii) with a pouch young.
1.5.3 The Eastern Grey Kangaroo (Macropus giganteus Shaw 1790)
The eastern grey kangaroo (Plate 1-3) is a crepuscular, polygynous species with a complex social organisation as they are gregarious, territorial and hierarchical (Jarman 1983). Eastern grey kangaroos exhibit pronounced sexual size dimorphism with males being up to four times larger than females at first oestrus (males: 19 – 85kg; females: 13 – 42kg), and males compete intensively for access to resources (Coulson 2008; Ganslosser 1989; Jarman & Southwell 1986). They are capable of breeding all year round, although there is a peak in the summer months, October to March (Poole 2002; Poole & Pilton 1964). This peak in breeding is associated with elevated levels of testosterone concentrations in male kangaroos, but with no significant seasonal variation in testis size (Nave 2002). Females are polyoestrous and monovular, but unlike the majority of species in the family Macropodidae, the gestation period is eight to nine days shorter than the oestrous cycle, and they do not exhibit post-partum
Chapter 1 General introduction 41 oestrus. Instead females return to oestrus 11 days following the loss of pouch young, or during the second half of lactation when the pouch young is greater than 160 days old (Poole & Pilton 1964). Eastern grey kangaroos live in small groups of less than six individuals that consist mainly of mothers and daughters, as young males generally disperse as young adults. At dawn and dusk, several groups whose home-ranges overlap, aggregate where there are concentrated resources, for example a preferred feeding area, forming a mob (Jarman & Southwell 1986).
Eastern grey kangaroos have a wide and almost continuous distribution through the eastern states of Australia, from QLD to Tasmania, from the inland plains to the coast (Poole 2002). They occupy a range of habitats including woodland, shrubland, open forest, and semi-arid mallee and mulga scrubs (Poole 2002) and have increased in number throughout their range in response to landscape modifications such as the conversion of forest and bush land to agriculture and the increase in permanent watering points (Calaby & Grigg 1989; Newsome 1965; Oliver 1986).
Management of kangaroos is a contentious issue in Australia as they are an iconic species. The proximity of many eastern grey kangaroo populations to urban areas makes management operations subject to public scrutiny as they are highly visible (Adderton Herbert 2004). Current management strategies include translocation, fertility control, commercial and non-commercial harvesting, fencing of areas, and euthanasia (Adderton Herbert 2004; Garrott et al. 1993). The formal policy framework for kangaroo management in Australia stipulates that culling should take place to alleviate or prevent the suffering of individuals (DEH 2007). The health and condition of the population is monitored, both prior to the cull to assess its necessity, and after the cull to monitor its success. Commercial harvesting of kangaroos is biased towards the largest males within the population (DEH 2007). Tenhumberg et al. (2004) attempted to model the impact of harvesting on kangaroos and found that size-selective harvesting could result in significantly smaller kangaroos for a given age when the entire population is subject to harvesting. Dispersal of individuals from non-harvested populations into the area may help to mitigate the genetic effects of harvesting, as well as maintain population size and structure (Tenhumberg et al. 2004) which is essential for maintaining their evolutionary potential (Ralls & Ballou 1986; Reed & Frankham
Chapter 1 General introduction 42
2003). But, we currently do not have an understanding of many of the variables necessary for input into these models, such as knowledge of the genetic mating system, the genetic substructure of the population and the heritability of the selected traits.
© E. Miller Plate 1-3 A female eastern grey kangaroo (Macropus giganteus) with a joey.
1.6 Study Aims
This thesis reports research on conservation genetics and reproduction in three Australian marsupial species. The broad aim was to use conservation genetic techniques to explore aspects of the three selected species life history (e.g. mating systems, Chapters 2 and 5), population structure (e.g. delineation of management units, Chapter 4) and the response of populations to management (e.g. captive breeding strategies, Chapter 3). In each case, the data were used to test current conservation genetics theories and to demonstrate how an understanding of population genetics is
Chapter 1 General introduction 43 critical for effective population management. In addition, the final research chapter (Chapter 6) sought to validate five commonly used body condition indices. Such indices are often employed to evaluate population health in the context of management and are an important tool in the field of evolutionary ecology by allowing us to objectively measure the quality of individuals. The specific aims of each chapter were:
(i) To investigate the impact of contrasting management strategies on the temporal patterns of genetic diversity for two captive breeding programs for a threatened species, the greater bilby (Chapter 2). (ii) To elucidate the genetic mating system, patterns of male reproductive success and the strength of sexual selection operating on male morphological traits in an elusive threatened species, the greater bilby (Chapter 3). (iii) To examine the population genetics and patterns of gene flow among three island populations of tammar wallabies and identify whether they are of particular genetic or conservation significance (Chapter 4). (iv) To identify whether particular behavioural, morphological, physiological, or genetic traits influence male reproductive success in the eastern grey kangaroo (Chapter 5). (v) To validate five commonly used body condition indices using serum biochemistry and haematology using the eastern grey kangaroo as a model species (Chapter 6).
Chapter 1 General introduction 44
Table 1-1 Comparison of the biology of the three Australian marsupials species selected for this study, the greater bilby (Macrotis lagotis), tammar wallaby (Macropus eugenii) and eastern grey kangaroo (M. giganteus). See section 1.5 for references for the below information.
Species Greater bilby Tammar wallaby Eastern grey kangaroo Conservation Vulnerable Low risk (near threatened) Least concern status
Distribution Once widespread, now confined to Remnant in WA mainland, several offshore Wide almost continuous between inland north-west Central Australia. islands in WA and SA. Extinct on SA and east coast of Australia Satellite populations in QLD mainland and Flinders Island
Habitat Semi-arid to arid Dense vegetation for daytime shelter; open Semi-mallee scrub through to woodland grassy area for feeding to forest; open grassy area for feeding
Active Nocturnal Crepuscular/nocturnal Crepuscular
Social Solitary Live in large mixed sex groups/solitary, not Gregarious, live in large mixed sex organisation territorial, hierarchical groups, territorial, hierarchical
Mating system Uncertain, polygynous or Polygynous Polygynous promiscuous
Sexual Yes Yes Yes - pronounced dimorphism Breeding All year round Highly seasonal, late Jan – early Feb with All year round, births peak in summer secondary peak late Feb – early March (Oct – Mar)
No. offspring 1 -3/litter 1 1
Management Captive breeding and Varied as some populations are threatened Overabundant local populations reintroductions, predator control and others are overabundant. Translocations, throughout range. Several management captive breeding and reintroductions, techniques including culling and fertility predator control and culling control
Chapter 2 Genetic diversity in contrasting captive breeding programs 45
Chapter 2 Bilbies behind bars: the impact of captive management on genetic diversity in a threatened species
2.1 Introduction
It is estimated that 2000-3000 terrestrial vertebrate species will require captive breeding over the next 200 years to save them from extinction (Frankham et al. 2002; Soulé et al. 1986; Tudge 1995). Already, many species have been preserved in captivity following their extinction in the wild, for example Prezwalski’s horse (Equus prezwalskii), the black-footed ferret (Mustela nigripes) and the Arabian oryx (Oryx leucoryx) (Frankham et al. 2002). For many other species, captive breeding programs have been established as a form of insurance against extinction in the wild. Captive populations provide a source of individuals for reintroductions and should be able to establish a self-sustaining wild population with high reproductive fitness and ample genetic diversity (Frankham et al. 2002). The International Union for Conservation of Nature (IUCN) recommends that captive populations be founded before wild populations drop below 1000 individuals. The advantages of establishing a captive population at this stage include using wild individuals with low inbreeding levels, minimising the impact of removing individuals from the wild and allowing sufficient time to develop suitable husbandry techniques (IUCN 2006).
The maintenance of genetic diversity is an important goal in captive breeding for several reasons. First, genetic diversity supplies the raw material for adaptive change (Darwin 1859; Frankham et al. 2002). Maintaining the adaptive potential of any
Chapter 2 Genetic diversity in contrasting captive breeding programs 46 population, wild or captive, ensures the population has the greatest chances of coping with environmental change. Second, higher levels of heterozygosity have been associated with higher fitness in plants, vertebrates and invertebrates (Reed & Frankham 2003). Third, inbreeding increases homozygosity and can lead to an accumulation of deleterious alleles and reduced fitness i.e. inbreeding depression (Frankham et al. 2002). Often captive breeding programs are founded with a small number of individuals, which can result in a loss of genetic diversity, inbreeding and reduced fitness. In turn, this will reduce a population’s adaptive potential thus lowering the likelihood of a population’s long-term survival (Ralls & Ballou 1986; Reed & Frankham 2003). The amount of genetic diversity lost over time is determined not so much by the number of individuals present in a population (N, census size), but by the genetically effective population size (Ne) (Wright 1969). Therefore Ne is an important consideration for management strategies.
Several management strategies have been recommended to retain maximum levels of genetic diversity and minimise levels of inbreeding in captive populations, such as the ‘maximum avoidance of inbreeding’ (MAI) and ‘minimising kinship’ (MK) strategies.
MAI involves the equalisation of family sizes and the doubling of Ne (Frankham et al. 2002). MK involves managing pedigreed populations so that the matings are made between individuals that are the most distantly related in the population. The latter method has been shown to be more effective in maintaining genetic diversity than alternative techniques such as random choice of parents or MAI, and also helps reduce unequal representation of genotypes from founder individuals (Montgomery et al. 1997). Despite this, no significant difference in reproductive fitness has been detected between MAI and MK techniques when experimentally tested using 40 replicate populations of the vinegar fly (Drosophila melanogaster) managed over four generations (Montgomery et al. 1997). This study compares the genetic diversity of two captive breeding programs of the greater bilby (Macrotis lagotis) that have implemented different management strategies, unmanipulated mating and minimising kinship.
Introducing novel genotypes into a population (‘genetic rescue’) through translocation and exchange of individuals between different captive breeding programs can also help
Chapter 2 Genetic diversity in contrasting captive breeding programs 47 alleviate the effects of inbreeding and loss of heterozygosity (Bryant et al. 1999; Spielman & Frankham 1992). ‘Genetic rescue’ is based on the assertion that immigrants into a population can introduce new genetic variation and increase the fitness in populations that are experiencing inbreeding depression (Darwin 1883; Ingvarsson & Whitlock 2000; Tallmon et al. 2004; Whitlock et al. 2000; Wright 1931, 1940). This has been tested in recent studies in both natural and experimental populations (Tallmon et al. 2004). For example, a reduction in the proportion of stillborns, increase in genetic diversity, recruitment and population growth in a population of adders (Vipera berus) has been attributed to the introduction of immigrants (Madsen et al. 1999). The improved fitness is thought to be largely due to heterosis in the offspring that arises from matings between immigrants and local individuals. However, if the immigrant is highly genetically divergent from the local population, this can lead to outbreeding depression and can result in a subsequent reduction in fitness (Tallmon et al. 2004).
Over the past 200 years, more mammals have become extinct in Australia than any other country (Short & Smith 1994). Within Australia, captive programs for many species have been established with the aim of reintroduction into the wild, for example the western barred bandicoot (Perameles bougainville), Gilbert’s potoroo (Potorous gilbertii), banded hare-wallaby (Lagostrophus fasciatus), mallee fowl (Leipoa ocellata), greater stick nest rat (Leporillus conditor) and the greater bilby (Mawson 2004a). The greater bilby is Australia’s only remaining species of desert bandicoot and an iconic species for conservation. They are burrowing, nocturnal marsupials, belonging to the subfamily Thylacomyinae (Johnson 2002). During the past century bilby populations have been in rapid decline. Prior to European settlement, bilbies were widespread throughout the arid and semi-arid zones of Australia, but now are confined to just 20% of their former range (Southgate 1990). It has been suggested that competition with introduced herbivores (cattle, sheep and rabbits) and predation by introduced red foxes (Vulpes vulpes) and feral cats (Felis catus) have reduced populations as well as contributed to the failure of some reintroduction attempts (Moritz et al. 1997). The greater bilby is listed as ‘Vulnerable’ to extinction by the IUCN (IUCN 2006).
Chapter 2 Genetic diversity in contrasting captive breeding programs 48
As part of the national recovery plan for the greater bilby, several captive breeding programs have been established with the management goal of retaining at least 90% of current genetic variation for 100 years (Pavey 2006). Two such programs have been established in Western Australia (Figure 2-1): Return to Dryandra (hereafter referred to as Dryandra) and Peron Captive Breeding Centre (hereafter referred to as Peron). Both colonies were established within 12 months of one another, each founded with seven individuals, and the same sex ratio (3 females: 4 males). Dryandra is located in the Dryandra Woodlands, an area where bilbies have become locally extinct. Management has simulated the natural environment and mating was unmanipulated under semi-free ranging conditions. In contrast, Peron maintained a pedigree with the aim of minimising the mean kinship between individuals to ensure the maximum retention of genetic diversity. Both programs were developed by Western Australia’s Department of Environment and Conservation (DEC) and form a component of the largest wildlife recovery program ever undertaken in Australia, Western Shield, which aims to expand predator control and then reintroduce native animals to their former habitats (CALM 1999).
Francis Peron National Park
Dryandra
Figure 2-1 Distribution of the greater bilby (Macrotis lagotis) in Australia (black = present distribution; mid-grey = historic distribution at European settlement; pale grey = Late-Holocene sub-fossil). The locations of the two captive breeding programs referred to in this study are indicated by black circles (after DEC 2004).
Chapter 2 Genetic diversity in contrasting captive breeding programs 49
To date, only two studies have examined the genetic variation of current and historically subdivided bilby populations (Moritz et al. 1997; Southgate & Adams 1993). There is no knowledge of how founder numbers and management strategies impact genetic diversity in captive breeding programs over time in a threatened marsupial. Such information is important for decisions influencing conservation and future management. The Dryandra and Peron captive colonies provide a unique opportunity to determine the impact of different management strategies on genetic diversity in this species, as well as being immediately applicable to captive breeding of any species in the world.
The aims of this study were to: (1) measure and compare genetic diversity within and between each captive colony overall, and compare these to wild populations; (2) monitor temporal changes in genetic diversity; and (3) assess the reliability of studbook estimates of genetic diversity and inbreeding compared to those calculated from microsatellite data.
2.2 Materials and Methods
2.2.1 Study populations
The Dryandra (32° 46’S, 116° 58’E) population was established in 1998. The breeding enclosure consists of 20 hectares (ha) of natural vegetation (Plate 2-1) surrounded by a 2.5m electrified fence (Plate 2-2). This is divided into two 10ha sections by a conventional fence, and one of these enclosures provides the basis of this study. The founder population consisted of five individuals of diverse origins from Kanyana (WA), and two wild caught individuals from north-western WA. Kanyana is a wildlife rehabilitation centre for wild bilbies that has a small breeding colony of individuals with a known pedigree. The founder sex ratio at Dryandra was three females to four males. Animals were semi-free ranging and mating was unmanipulated within their enclosure. The Peron population was established in 1997, and is located on 1050km2 of the Peron Peninsula (25° 42’S, 113° 32’E), which is severed from the mainland at its isthmus by a 3.4km electric fence. Animals were housed in small mesh-covered pens to accommodate breeding pens and large outdoor pens to accommodate small family groups. The founder population consisted of four individuals from Kanyana and three
Chapter 2 Genetic diversity in contrasting captive breeding programs 50 wild caught individuals, and all founders originated from north-western WA. The founder sex ratio was three females to four males and there was a known pedigree based on the breeding of the founder individuals. Animals were selectively bred to minimise kinship. Both populations had additional individuals translocated into the colony over time from Kanyana, or the wild, since establishment, as well as animals periodically being removed for reintroduction into the wild. In both populations, all individuals were identified using subcutaneous PIT tags (Destron, USA). Tissue samples were collected from individuals upon first capture as an independent juvenile or upon introduction into the colony and were stored in dimethyl sulfoxide (DMSO) prior to DNA extraction.
© E. Miller
Plate 2-1 The Dryandra woodlands, Western Australia, that surrounds the Return to Dryandra captive breeding facility and serves as a reintroduction site for the locally extinct greater bilby (Macrotis lagotis).
Chapter 2 Genetic diversity in contrasting captive breeding programs 51
© E. Miller
Plate 2-2 Predator proof fencing surrounding the Return to Dryandra captive breeding facility for threatened species in Western Australia.
2.2.2 Microsatellite genotyping
DNA was extracted from 2mm ear biopsies (Dryandra n = 216; Peron n = 266) using a standard salting-out method (Sunnucks & Hales 1996; Appendix 1). To enable a comparison between the captive populations and WA wild-caught individuals, frozen blood homogenates were obtained from the South Australian Museum from six wild caught bilbies from northern WA (lodged by R. Southgate). DNA was extracted from the blood homogenates using a FlexiGene DNA kit (Qiagen, Australia). Each individual was screened at nine polymorphic microsatellite loci (Bil02, Bil16, Bil17 Bil22, Bil41, Bil55, Bil56, Bil63 and Bil66) previously characterised from the greater
Chapter 2 Genetic diversity in contrasting captive breeding programs 52 bilby (Moritz et al. 1997). Genotyping was carried out using multiplexed PCRs, whereby fluorescently labeled primers enabled the simultaneous amplification of many targets of interest, using several pairs of primers in a single reaction. Two multiplex combinations were used (i) Bil55, Bil22, Bil16, Bil41, Bil17 and (ii) Bil02, Bil56, Bil63, Bil66. All loci were amplified in a “touchdown” PCR, whereby the initial annealing temperature was decreased at 1ºC increments and run for a total of ten cycles each. Ten micro-litre reactions were performed using a Multiplex PCR kit (Qiagen, Australia), 0.2μM of each primer and 60 – 80ng genomic DNA. Amplifications were carried out in a PTC-220 thermocycler (MJ Research, USA) using an initial HotStarTaq activation step at 95ºC for 15 min, a total of 50 cycles of 94ºC for 30 s, 60, 59, 58, 57, 56 and 55ºC for 1 min 30 s (10 cycles each), 72ºC for 1 min 30 s, and a final extension at 72ºC for 10 min. The PCR products were analysed in a 48 capillary AB 3730 DNA Analyser (Applied Biosystems, USA). The DNA fragments were sized and quantified using GeneMapper 3.7 (Applied Biosystems, USA).
2.2.3 Genetic analyses
The genotype file was checked for duplicate entries, scoring errors due to null alleles, short allele dominance and stutter bands using the program MICRO-CHECKER (van Oosterhout et al. 2004). Locus independence and Hardy-Weinberg equilibrium tests were conducted using GenePop 3.4 (Raymond & Rousset 2003) using a Markov chain method (1000 iterations). The statistical significance levels were corrected for multiple comparisons using sequential Bonferroni adjustments (Rice 1989). Genetic diversity was estimated for each population by calculating allelic diversity (AD), observed (Ho) and expected (He) heterozygosities using GENALEX 6.0 (Peakall & Smouse 2006). A Wilcoxon signed rank test was used to detect significant differences in diversity between Dryandra and Peron. To assess temporal changes in genetic diversity within each population, the dataset for each colony was further divided into annual ‘sub- populations’. Each sub-population represents every independent individual present in the population in a 12 month period. Due to the small number of individuals present in the first two years for both populations (1998 and 1999 for Dryandra (98/99) and 1997 and 1998 for Peron (97/98)), these years were pooled and examined as a single ‘sub- population’. A One-Way ANOVA using Tukey’s Post-Hoc Multiple Comparisons was used to test for differences in AD, Ho and He between years within each colony.
Chapter 2 Genetic diversity in contrasting captive breeding programs 53
The F-statistics, inbreeding coefficients (FIS) and genetic differentiation (FST), were calculated using Weir & Cockerham (1984) estimators from the microsatellite data using FSTAT 2.9.3.2 (Goudet 2001). The significance levels of FIS and FST were determined after 10 000 permutations. Population bottlenecks were tested for using BOTTLENECK 1.2.02 (Cornuet & Luikart 1996), assuming a two-phase model (TPM) with a 95% single-step mutations and 5% multiple step mutations, and a variance among multiple steps of 12. A Wilcoxon test was used as suggested for relatively low number of loci (Piry et al. 1999).
Microsatellite variation (AD and He) of both captive populations was compared to published data for the same microsatellite loci for wild M. lagotis populations from the Northern Territory (NT; n = 19) and Queensland (QLD; n = 27) (Moritz et al. 1997) using a Wilcoxon signed rank test in SPSS 15.0. To enable this comparison, Bil66, which was out of Hardy-Weinberg equilibrium for this study, was included in the analysis. Note Bil63 was also out of Hardy-Weinberg equilibrium for the wild QLD population (Moritz et al. 1997).
Population data for Peron was collected by DEC (WA) and maintained in SPARKS (ISIS 1992). Two pedigree-based statistics were generated using PM2000 (Pollak et al. 2002): (i) mean inbreeding coefficients (F) were calculated using SPARKS (ISIS 1992) and corroborated using PM2000; and (ii) retained gene diversity, that is, the proportion of heterozygotes expected in the descendant population under Hardy- Weinberg equilibrium assumptions relative to the founding population. Retained gene diversity is calculated by allocating two hypothetical alleles to each founder and running a gene-drop simulation analysis (1000 iterations) to determine what proportion of the founding genetic diversity is present in the current population according to Mendelian inheritance (Lacy & Ballou 2002). The studbook estimates were compared to the genotypic estimates of diversity using a Wilcoxon signed rank test. All statistical analyses were conducted in SPSS 15.0.
Chapter 2 Genetic diversity in contrasting captive breeding programs 54
2.3 Results
2.3.1 Genetic diversity
Eight of the nine loci amplified successfully (Bil17 failed to amplify). MICROCHECKER found five duplicate entries of individuals that had been sampled twice, and so the duplicates were removed from the dataset. All loci except Bil66 were in Hardy-Weinberg equilibrium. In both captive populations Bil66 showed a significant excess of homozygotes (p > 0.05) and evidence of null alleles was detected by MICROCHECKER. There was no evidence of allelic dropout in any locus. Bil17 and Bil66 were excluded from further analyses. The two captive bilby populations were polymorphic for all other loci. For both populations, across the seven loci, values of Ho varied between 0.394 and 0.823 and for He between 0.483 and 0.848 (Table 2-1). There were four to nine alleles per locus. One private allele was identified in the Dryandra population and nine were identified in the Peron population. The annual population allele frequency data for Dryandra and Peron are presented in Appendix 2 and 3, respectively. The overall allele frequency data comparing Dryandra and Peron are presented in Appendix 4.
Table 2-1 Summary of the overall genetic diversity calculated across all years from seven microsatellite loci for two captive bred populations of the greater bilby (Macrotis lagotis), Return to Dryandra (Dryandra) and Peron Captive Breeding Centre (Peron), Western Australia (mean ± se).
Parameter Dryandra Peron n 216 266 AD 6.4 (± 0.6) 7.6 (± 0.6) Range AD 4 – 8 5 – 9
Mean Ho 0.681 (± 0.047) 0.664 (± 0.055)
Range Ho 0.463 – 0.823 0.394 – 0.763
Mean He 0.709 (± 0.041) 0.758 (± 0.030)
Range He 0.483 – 0.813 0.628 – 0.848
Sample size (n), allelic diversity (AD), observed heterozygosity (Ho) and expected heterozygosity (He).
Chapter 2 Genetic diversity in contrasting captive breeding programs 55
2.3.2 Temporal changes in genetic diversity
Changes in allelic diversity (AD) and heterozygosity (He) over time for each population are shown in Figures 2-2 (a) and (b), respectively. The Dryandra colony showed an initial decline in AD, followed by a slight increase then plateau. The Peron colony showed an initial increase in allelic diversity then remained fairly constant until 2004, when a slight reduction in AD is evident. Dryandra has significantly lower AD than Peron (WRS z = -2.226, p = 0.026). The level of He in both populations has remained relatively constant since establishment, with a gradual increase in He evident in the Peron colony, but no significant difference to Dryandra (WRS z = -1.823, p = 0.068). In both colonies there were no significant differences between years for AD
(Dryandra: F6, 42 = 0.169, p = 0.984; Peron: F9, 60 = 1.062, p = 0.403) or He (Dryandra:
F6, 42 = 0.171, p = 0.983; Peron: F9, 60 = 1.404, p = 0.207).
The levels of inbreeding (FIS) over time are shown in Figure 2-3. The Dryandra population became increasingly outbred three years after establishment (mean FIS = 0.080; range = 0.016 – 0.123). This trend continued until 2003 when the levels of inbreeding began to increase. A significant FIS was detected in 2005 (p = 0.001). The positive values for Peron indicate more homozygotes being present in the population
(mean FIS = 0.136; range = -0.035 – 0.075), and FIS was significant in 2003, 2004 and 2005 (p = 0.001 each year). The number of individuals translocated into each population each year is indicated on this graph.
Chapter 2 Genetic diversity in contrasting captive breeding programs 56
(a) 10.0
9.0 )
AD 8.0
7.0
6.0 Allelic diversity ( diversity Allelic 5.0
4.0 97/98 98/99 2000 2001 2002 2003 2004 2005 2006 Year
(b)
1.0
) 0.9 e H 0.8
0.7
0.6
Mean heterozygosity ( 0.5
0.4 97/98 98/99 2000 2001 2002 2003 2004 2005 2006 Year
Figure 2-2 Genetic diversity over time for the greater bilby (Macrotis lagotis), captive breeding programs, Return to Dryandra (Dryandra, ) and Peron Captive Breeding Centre (Peron, ), Western Australia. (a) Mean allelic diversity (AD) and (b) mean heterozygosity (He).
Chapter 2 Genetic diversity in contrasting captive breeding programs 57
0.30 12 3 2 2 2 ) IS F 0.10
97/98 98/99 2000 2001 2002 2003 2004 2005 2006 -0.10 Inbreeding( coefficient
15 16 8 9 13 3 -0.30 Year
Figure 2-3 Inbreeding coefficients (FIS) over time for two captive bred populations of greater bilbies (Macrotis lagotis); Return to Dryandra (Dryandra, ) and Peron Captive Breeding Centre (Peron, ), Western Australia. The dashed arrows indicate when new individuals were translocated into each colony (Dryandra below graph, Peron above graph) and the number represents the number of individuals.
The overall genetic differentiation between Dryandra and Peron was low, but significant (FST = 0.056). There was also significant genetic differentiation within each captive colony over time (Dryandra, Table 2-2(a); Peron Table 2-2(b)). Dryandra showed no evidence of a genetic bottleneck in any year, whereas Peron showed evidence of a genetic bottleneck only in 2006 (p = 0.019).
Chapter 2 Genetic diversity in contrasting captive breeding programs 58
Table 2-2 Pair-wise genetic differentiation (FST) partitioned over time for two captive bred populations of greater bilbies, Macrotis lagotis, (a) Return to Dryandra (Dryandra) and (b) Peron Captive Breeding Centre (Peron), Western Australia.
(a) Pop 2000 2001 2002 2003 2004 2005 98/99 0.002 0.009 0.013 0.008 0.011 0.013 2000 0.015 0.026* 0.018* 0.021* 0.020* 2001 0.001 0.007 0.011* 0.022* 2002 -0.001 0.004 0.017* 2003 -0.004 0.003 2004 -0.002 Significant differentiation indicated by *
(b) Pop 1999 2000 2001 2002 2003 2004 2005 2006 97/98 -0.007 0.007 -0.009 0.015 0.054* 0.067* 0.075* 0.083* 1999 -0.006 -0.005 0.019* 0.048* 0.055* 0.066* 0.067* 2000 0.000 0.030* 0.055* 0.058* 0.073* 0.070* 2001 0.016* 0.047* 0.060* 0.070* 0.075* 2002 0.004 0.017* 0.027* 0.047* 2003 0.004 0.022* 0.052* 2004 0.004 0.023* 2005 0.003 Significant differentiation indicated by *
2.3.3 Genetic diversity in comparison to wild populations
Dryandra had significantly lower levels of He (WRS z = -2.240, p = 0.025) and AD (WRS z = -1.761, p = 0.033) than the wild NT population. There were no significant differences in genetic diversity between Peron and NT, or between both captive populations and the wild QLD population (Table 2-3). Due to the small sample size of wild WA individuals (n = 6), the data could not be analysed statistically. The wild WA bilbies shared 100% of their alleles with the captive individuals: 78.9% were shared with Dryandra individuals and 100% were shared with the Peron population. Three
Chapter 2 Genetic diversity in contrasting captive breeding programs 59 alleles were present in the wild WA and Peron individuals that were not present in the Dryandra population.
Table 2-3 Genetic diversity for wild Northern Territory and Queensland populations, in comparison to the captive Western Australian (WA) populations of greater bilbies (Macrotis lagotis).
Population n Mean AD Mean He Northern Territory* 19 8.9 0.813 Queensland* 27 7.7 0.654 Dryandra (WA) 218 6.4 0.709 Peron (WA) 266 7.6 0.758
Sample size (n), allelic diversity (AD) and expected heterozygosity (He). *published data from Moritz et al. (1997).
2.3.4 Genetic diversity in comparison to studbook estimates
Overall the studbook estimates of genetic diversity for Peron were significantly higher than those estimated from the genotypic data (WRS z = -2.192 p = 0.028), and show the reverse trend (Figure 2-4(a)). The studbook estimates suggest that diversity slowly decreased over time, whereas the genotypic data shows a gradual increase. With higher estimates of genetic diversity, the studbook estimates of inbreeding were significantly lower than those calculated from the genotypic data (WRS z = -2.666 p = 0.008; Figure 2–4(b)), but indicated a similar trend in inbreeding over time.
Chapter 2 Genetic diversity in contrasting captive breeding programs 60
(a)
1
0.9 ) e H 0.8
0.7
0.6
Heterozygosity ( Heterozygosity 0.5
0.4 1999 2000 2001 2002 2003 2004 2005 2006 2007 Year
(b)
0.3 )
IS 0.2 F 0.1
0 1999 2000 2001 2002 2003 2004 2005 2006 2007 -0.1
Inbreeding ( coefficient -0.2
-0.3 Year
Figure 2-4 Comparison of the studbook (dashed line) and genotypic (solid line) estimates of (a) genetic diversity, and (b) inbreeding coefficients between 1999 and 2007 for the Peron Captive Breeding Centre (Peron) colony, Western Australia.
Chapter 2 Genetic diversity in contrasting captive breeding programs 61
2.4 Discussion
Maximising and maintaining genetic diversity is an important goal for captive breeding programs to ensure long-term population sustainability. This study found that (i) genetic diversity was maintained over time in both populations despite the different management strategies implemented; (ii) the supplementation of new individuals into the colonies provided a new source of genetic variation and aided the maintenance of genetic diversity; (iii) Peron had similar levels of genetic diversity to wild NT bilbies, but Dryandra had significantly lower diversity, and (iv) the studbook estimates of genetic diversity were significantly higher, and inbreeding significantly lower than those calculated from the genotypic data.
2.4.1 Genetic diversity in the captive bilby populations
The levels of allelic diversity and heterozygosity remained relatively constant in both captive breeding colonies, although Dryandra had significantly lower allelic diversity than Peron. In captive populations, genetic variation arises from the contribution of founding individuals, the introduction of immigrants and mutation (Ballou 1984). In small, closed populations mutation rates are likely to be too low to generate variation (Ballou 1984; Hedrick 2000). In addition, allelic diversity is strongly linked to population size (Frankham 1996) and small closed populations are expected to lose diversity due to random genetic drift (Lacy 1989). Consequently, we would expect the erosion of genetic diversity over time due to genetic drift in these captive bilby populations. However, we observed that genetic diversity was maintained over time in both the Dryandra and Peron captive breeding colonies, though at different levels.
Both captive bilby colonies started with similar levels of allelic diversity and in the year following establishment received a similar number of new individuals, yet after 12 months Peron had significantly higher allelic diversity than Dryandra. This difference possibly arose because those individuals introduced into Peron contributed more novel diversity than those introduced into Dryandra. Founding populations with unrelated individuals is important to minimise the loss of genetic diversity and inbreeding at the foundation for long-term conservation (Frankham et al. 2002; Gautschi et al. 2003). The strategy of minimising kinship through the use of a pedigree
Chapter 2 Genetic diversity in contrasting captive breeding programs 62 at Peron may have also contributed to the higher levels of diversity retained through initially equalising the founder contribution, that is, limiting the reproduction and genetic representation of the founder individuals (Ballou & Lacy 1995; Loebel et al. 1992). Equalisation of founder contributions is thought to diminish genetic drift by enlarging effective population size, thus resulting in a higher retention of allelic diversity (Lacy 1989). However, it is difficult to determine which management strategy has been more effective at maintaining diversity since both populations continued to translocate individuals into the population over time.
The translocation of new individuals into the populations may have buffered the effects of inbreeding and genetic differentiation within each colony. Between 1999 and 2005, 57 new individuals were translocated into the Dryandra colony. There was no significant differentiation between the founding and established population in 2005. Dryandra showed some evidence of significant (but low level) inbreeding in 2005, but there was no evidence of a genetic bottleneck. In contrast, Peron translocated fewer individuals (n = 21) between 1999 and 2004, experienced higher levels of significant inbreeding (2003 – 2005) and showed evidence of the beginnings of a genetic bottleneck by 2006. In addition, by 2003, Peron had significantly differentiated from its founding population. Similarly, a study of the bridled nailtail wallaby (Onychogalea fraenata) found captive-bred individuals and their wild-born offspring were significantly differentiated from the wild remnant populations within four generations as a result of rapid genetic drift arising from a small number of founders (n = 7) and subsequent loss of allelic diversity (Sigg 2006). They also showed that when individuals were removed from the captive breeding program for reintroduction into the wild, allelic diversity, but not heterozygosity decreased (Sigg 2006). If individuals are being removed from a population genetic diversity will decline if new genetic stock is not added to the population (Frankham 1996; Nei 1987; Sigg 2006).
There are substantial benefits for small, partially inbred populations when even a single unrelated animal is translocated into the population, for example increases in reproductive fitness (Bryant et al. 1999; Spielman & Frankham 1992). The strategy of transferring individuals among captive breeding facilities has proven to be beneficial for other species, for example, the bearded vulture, Gypaetus barbatus (Gautschi et al.
Chapter 2 Genetic diversity in contrasting captive breeding programs 63
2003), and a single immigrant increased heterozygosity and led to a rapid spread of new alleles in a re-founded population of the Scandinavian wolf, Canis lupus (Ingvarsson 2003; Vilà et al. 2003). Exchanging individuals with other facilities is an effective strategy for maintaining genetic diversity, but populations should continue to be monitored for genetic changes. There is evidence from salmon that suggests that breeding programs designed to supplement wild fisheries can reduce the fitness of natural populations through a reduced effective population size, mutation accumulation, and genetic adaptation to captivity (Araki et al. 2007; Goodman 2005; Heath et al. 2003; Wang & Ryman 2001).
2.4.2 Comparison to wild populations
In comparison with wild bilby populations, Peron had similar levels of diversity, but Dryandra had significantly lower heterozygosity and allelic diversity than the wild NT population but not wild QLD bilbies. When comparing the wild WA bilby genotypes with the captive populations, the results showed that Dryandra had not captured as much wild genetic variation as it potentially could have. A strategy to improve the representation of wild bilby diversity in Dryandra would be to source more wild caught individuals and/or exchange more individuals with Peron. However, the levels of diversity in both captive populations were within the range of values observed in other marsupial taxa including other ‘Vulnerable’ taxa (Table 2-4). As most captive populations are small, and/or founded with small numbers, it is not surprising that evidence suggests captive populations generally have lower genetic diversity than wild populations (Jiang et al. 2005; Kubota et al. 2008; Neveua et al. 1998). In the context of other captive and wild species, the genetic diversity in these captive populations of the greater bilby is reasonable.
Chapter 2 Genetic diversity in contrasting captive breeding programs 64
Table 2-4 Comparison of the genetic diversity in the captive greater bilby (Macrotis lagotis) populations with other wild marsupial populations of varying conservation status (after Bowyer et al. (2002)).
IUCN category No. populations AD He Least concern 8 12.0 – 5.2 0.86 – 0.66 Near threatened 4 11.1 – 5.3 0.85 – 0.60 Vulnerable* 3 8.9 – 6.0 0.86 – 0.65 Endangered 2 11.6 – 6.0 0.83 – 0.72 Critically endangered 1 1.8 0.27
Allelic diversity (AD) and heterozygosity (He). * Dryandra mean He = 0.709, AD = 6.4;
Peron He = 0.758, AD = 7.6.
2.4.3 Comparison to studbook estimates
In the Peron population there was a significant difference between studbook and genotypic estimates of genetic diversity and inbreeding. The studbook calculations based on pedigree data significantly overestimated genetic diversity, and consequently significantly underestimated levels of inbreeding. The studbook estimates also showed a declining trend in genetic diversity over time, but the reverse trend was apparent in the genotypic data. The Peron management program aimed to retain at least 90% of the genetic diversity found in the wild or from the source population. The estimates of genetic diversity from the studbook indicated that the Peron population was on target initially, but as time progressed, the diversity in the population decreased by 18.8% compared to the target of losing no more than 10%. In contrast, the genotypic data showed a significantly lower level of diversity in the population, but with genetic diversity increasing over time. As expected, the reverse was true for inbreeding. The studbook calculation for inbreeding was consistently lower than that estimated from the microsatellite data.
A difficulty that faces most captive managed populations is that the relationship amongst the founders is often unknown, and they are assumed to be unrelated. This assumption of unrelatedness between founder individuals will affect studbook calculations (Nielsen et al. 2007). If the founders are related, the studbook calculations
Chapter 2 Genetic diversity in contrasting captive breeding programs 65 are likely to overestimate genetic diversity, and consequently underestimate inbreeding in the population, as shown in this study. There is conflicting evidence in the literature regarding the accuracy of studbook estimates compared to genotypic estimates. While some studies have found them to be inaccurate (Haig et al. 1994; Jones et al. 2002; Morin & Ryder 1991), others have found studbook and genotypic estimates to be congruent and therefore suitable as a tool for genetic management (Nielsen et al. 2007; Wisely et al. 2003). As well as incorrect assumptions regarding the interrelationship of founders, these discrepancies can arise from fundamental differences in how the measures are calculated, inaccurate or incomplete studbook data, genotyping data, or uninformative genetic loci. These results highlight the importance of ensuring that founding individuals are unrelated prior to commencement and validation of studbook estimates of diversity with genotypic data as they form the basis for the genetic management of many threatened populations.
2.4.4 Implications for management
Species loss and the threat of extinction is a worldwide problem. Many species require captive breeding, and utilise studbook estimates of diversity to gauge the ‘genetic health’ of the captive population. Molecular genetics is a tool that has facilitated the evaluation of the effectiveness of management decisions and their impact on genetic diversity. In small, closed populations, founder contributions and new immigrants into a population are important for providing a new source of genetic variation. However, where possible, founders should be tested for relatedness. This will not only facilitate maximisation of genetic diversity, but also satisfy an underlying assumption of studbook diversity calculations which may improve the accuracy. In the case of the greater bilby, recruiting additional unrelated founder individuals and increasing the gene flow between captive bilby populations may assist in maintaining genetic variation in both the captive and the wild reintroduced populations. This will reduce the effects of genetic drift, genetic adaptation to captivity, inbreeding and genetic differentiation among populations. An additional benefit of exchanging individuals among populations could be further increasing the effective population size, thus maintaining the evolutionary potential of the greater bilby.
Chapter 2 Genetic diversity in contrasting captive breeding programs 66
The translocation of individuals between populations, zoos and wildlife institutions can provide ‘genetic rescue’ and be an effective conservation tool to mitigate the effects of small founder numbers, unequal founder contributions, inbreeding, genetic differentiation and genetic drift, but wildlife managers need to be aware of the potential problems that can occur. In conjunction with genetic management, breeding programs should utilise biological data to maximise genetic diversity. For example, a species mating system can influence the genetic structure of a population. If the mating system is polygynous and there is a large reproductive skew with only a few males contributing to the gene pool, an effective strategy would be to manipulate mating patterns by using a specific number of sires (Oyama et al. 2007). In a promiscuous mating system, more males participate in breeding, lowering the variance in male reproductive success which in turn would increase the effective population size and slow rates of inbreeding. An increase in the number of sires should lead to lower levels of relatedness within the population. Knowledge of the bilby mating system will help to further elucidate the relative influence of captive breeding strategies, genetic rescue and species mating systems on the maintenance of genetic diversity within captive breeding programs.
2.4.5 Conclusions
In summary, this study found the levels of genetic diversity were maintained over time in both captive breeding programs for the threatened greater bilby. Since both colonies translocated new individuals into the populations regularly, it was difficult to determine which strategy, if any, was more effective in maintaining genetic diversity. The introduction of new individuals helped mitigate the risk of genetic erosion, inbreeding and genetic differentiation that is expected to occur in small, closed populations founded with a small number of individuals. Although Dryandra had lower levels of diversity than the wild NT bilby population, the overall level of genetic variability in the captive bilby populations were comparable to that of wild bilby populations and other marsupial taxa. The studbook estimates of genetic diversity were overestimated in comparison to those calculated by the microsatellite data, and thus, levels of inbreeding were underestimated. The use of studbook genetic estimates should not be solely relied upon to evaluate the genetic health of the captive bilby population. Captive breeding programs should validate their studbook estimates with
Chapter 2 Genetic diversity in contrasting captive breeding programs 67 population specific genotypic data. This study highlights the importance of replenishing captive populations with new stock, especially post-animal removal for reintroductions.
Chapter 3 Genetic mating system, reproductive success and selection 68
Chapter 3 The genetic mating system, male reproductive success and selection on male traits in the Greater Bilby (Macrotis lagotis)
3.1 Introduction
Captive breeding programs are likely to become an increasingly important component of conservation strategies for terrestrial vertebrates in the future. Several species have already been preserved in captivity following their extinction in the wild, for example the California condor (Gymnogyps californianus) and Père David’s deer (Elaphurus davidianus) and it has been estimated that between 2000 and 3000 terrestrial vertebrate species will require captive breeding to prevent their extinction in the next 200 years (Frankham et al. 2002; Soulé et al. 1986; Tudge 1995). Captive breeding programs also serve as a form of insurance against wild extinction, providing individuals for wild reintroductions. Ideally, individuals selected for reintroduction should be physically healthy with a known high reproductive output and abundant genetic variation (Frankham et al. 2002), but this information is often not available because of both practical and financial constraints. In the absence of this information, individuals may be selected for release based on their body size as there is evidence in several species that males with a larger body size often have enhanced reproductive success (Andersson 1994; Birkhead 2000). Fundamental to successfully meeting this criterion is an understanding of the target species biology, including knowledge of their behaviour, reproduction, social organisation and mating system.
Chapter 3 Genetic mating system, reproductive success and selection 69
Knowledge of the mating system and parentage is crucial for the effective management of captive breeding programs in order to determine effective population sizes, detect inbreeding, genetic differentiation, and genetic diversity, all of which affect population growth and survival (Ralls & Ballou 1986; Reed & Frankham 2003). Parentage cannot be reliably determined using behavioural data alone. Advances in genetics have revealed inconsistencies in the presumed relationship between social organisation, parentage and mating system (Ambs et al. 1999; Coltman et al. 1999a; Worthington Wilmer et al. 1999), revolutionising how mating systems are understood. For example, around 95% of avian mating systems are classified as socially monogamous (Schwagmeyer & Ketterson 1999), however, genetic data has shown that extra-pair fertilisations occur on average in 13% (range 0 – 76%) of species (Westneat & Sherman 1997), leading to a better understanding of the evolution of sexual selection, mate choice and parental care. This study uses genetic markers to examine the mating system, male reproductive success and strength of selection on male morphological traits in a cryptic, threatened species that is the subject of a captive breeding and reintroduction program.
The diversity in mammalian mating systems is a product of the reproductive strategies of individuals and their environments, for example how many individuals males or females mate with and whether pair bonds form (Clutton-Brock 1989; Emlen & Oring 1977; Reynolds 1996). The classic mating system classifications are monogamy, polygyny, polyandry and promiscuity, though the precise definition of each can differ even within the field of evolutionary biology (Andersson 1994). Typically, monogamy is defined as the long-standing bond and exclusive relationship between a male and a female. Polygyny occurs when a single male mates with several females, and polyandry is the opposite, whereby a single female mates with several males (Wittenberger 1981). In a promiscuous mating system both males and females mate with more than one partner and no long term relationship forms (Wittenberger 1979). Polygyny commonly occurs in species where males do not provide parental care (i.e. most mammals). Consequently, males are able to invest their energy into competing for resources and mates (Emlen & Oring 1977). The ability to monopolise resources is a key factor in the intensity of sexual selection. Polygyny is expected to occur when there is spatial clumping of females in such a way that enables males to defend
Chapter 3 Genetic mating system, reproductive success and selection 70 multiple females from being accessed by other males (Clutton-Brock 1989; Emlen & Oring 1977) and is characterised by greater variance in male reproductive success than in females (Hoogland & Foltz 1982). Polyandrous mating systems are characterised by females pairing with multiple males during a single breeding season, but true polyandry has only been documented in a few mammalian species (Wittenberger 1979). In contrast, promiscuity generally occurs when males are unable to successfully monopolise access to females. It often occurs in species where female groups are unstable, males provide no parental care, and male home ranges are distributed throughout the home ranges of several females (Clutton-Brock 1989). Males are often unable to defend territories and therefore cannot monopolise access to females (Ramsay & Stirling 1986).
According to sexual selection theory, males ought to compete for access to receptive females, the most limiting resource for male reproductive success (Trivers 1972). There is evidence in numerous mammalian taxa demonstrating that a larger body size is advantageous when competing for access to receptive females (Clinchy et al. 2004; Fisher & Cockburn 2005; Fisher & Lara 1999). The male-biased sexual size dimorphism that exists in many mammalian taxa is often attributed to competition between individuals for reproductive opportunities, encouraging the evolution of secondary sexual traits (Andersson 1994; Birkhead 2000). Female preference for particular traits could potentially increase the strength of selection on male traits. For example, if females prefer larger males because their size indicated quality, strength or fighting ability (Andersson 1994; Bercovitch et al. 2003; Coltman et al. 2001; Schulte- Hostedde & Millar 2002), then we would expect strong selection for body size. Variation in body size may also influence male reproductive success through size- assortative mating. In several mammals, a positive correlation between mother’s body weight and litter size has been detected, for example possums (Julien-Laferriere & Atramentowicz 1990), opossums (Hossler et al. 1994), rodents (Kaufman & Kaufman 1987; McClure 1981; Myers & Master 1983; Svendsen 1964) and squirrels (Neuhaus 2000; Risch et al. 2007). Given the higher fecundity of larger females, positive size- assortative mating theoretically should differentially increase the fitness of larger males. Under this scenario, large male body size is not important per se, but rather male body size relative to female body size.
Chapter 3 Genetic mating system, reproductive success and selection 71
The greater bilby (Macrotis lagotis) is an iconic species for conservation in Australia. They are the sole remaining species of desert bandicoot and are listed as ‘Vulnerable’ to extinction (IUCN 2006). Bilbies are now restricted to 20% of their former distribution in the arid and semi-arid zones of Australia (Southgate 1990) due to habitat loss, competition with introduced species and predation (Moritz et al. 1997). Bilbies are distinguished by their soft, silky hair, long rabbit-like ears and long furry black and white tail. They are an omnivorous, nocturnal marsupial that live solitary or in pairs. They dig burrows up to two metres deep, and individuals may have up to a dozen burrows they use within their home range (mean: males = 316 ± 128ha; females = 18 ± 4ha) (Moseby & O'Donnell 2003). Bilbies are sexually dimorphic with males being larger than females (males: 800 – 2500g; females: 600 – 1100g). They are capable of breeding all year round, producing one to three offspring per litter, and up to four litters per year under ideal conditions. Males do not provide parental care (Southgate et al. 2000).
The bilby mating system is not well documented due to their secretive nature and the difficulties associated with wild observations. They are thought to be polygynous or promiscuous (Lee & Cockburn 1985). Johnson & Johnson (1983) monitored the behaviour of three adult males, two adult females and two female young in captivity (18m x 12m pen). They observed a rigid dominance hierarchy, maintained with little destructive aggression. Dominant males maintained access to all the well used burrows in the enclosure and chased subordinates out of the burrow (Johnson & Johnson 1983). Moritz et al. (1997) used genetic data to perform a parentage exclusion analysis for a colony of wild Queensland bilbies (n = 17). The partial pedigree constructed was suggestive of bilbies being strongly polygynous as one male mated with three females to sire seven of the eight offspring (Moritz et al. 1997).
Several bilby captive breeding programs have been established across Australia with the aim of wild reintroduction. Return to Dryandra (RTD) is a program located within the Dryandra Woodland, Western Australia (WA), an area where bilbies have become locally extinct. This program was developed by Western Australia’s Department of Environment and Conservation (DEC) and forms a component of the Western Shield
Chapter 3 Genetic mating system, reproductive success and selection 72 program, which aims to expand predator control and then reintroduce native animals to their former habitats (CALM 1999). Within the RTD population, the male selection criteria for reintroduction back into the wild is based the assumption that large males are monopolising mating opportunities within the large breeding enclosure. The removal of reproductively dominant males from the captive colony should theoretically aid the maintenance of genetic diversity both within the colony and within the reintroduced population, by enabling other captive males to gain fertilisations, and by releasing a healthy individual with a high reproductive potential. Therefore it is vital to understand the breeding biology and population genetics of the greater bilby to determine if this is an appropriate strategy for effective management of captive breeding and reintroduction programs.
This study aims to (i) clarify the mating system of the greater bilby using microsatellite markers; (ii) examine the variance in male reproductive success and determine whether a small number of males monopolise paternity; (iii) examine whether morphological traits are associated with male reproductive success; and (iv) determine if there is strong sexual selection acting on male traits to enhance their reproductive success. Most mating system research to date has been conducted in eutherian mammals, and marsupials have received little attention. This study makes an important contribution to the future conservation and management of the greater bilby, and adds to the growing body of scientific literature about threatened species.
3.2 Materials and Methods
3.2.1 Study population, data and sample collection
The RTD breeding facility is located in south-west WA (32º 46’S, 116º 58’E) and consists of two 10ha enclosed areas of natural vegetation, surrounded by a 2.5m electrified fence. This study was based on samples collected from one of the 10ha enclosures between 1999 and 2005. Food was supplied to the animals in nine small and four large feed hoppers in the enclosure. Animals were trapped quarterly between 1999 and 2005 using small cage traps covered with a Hessian sack and baited in the afternoon with a mixture of peanut butter, oats and anchovies. The traps were checked at 2200 and the following morning at sunrise. All animals were individually identified
Chapter 3 Genetic mating system, reproductive success and selection 73 using PIT tags (Destron, USA) inserted under the skin. A small ear biopsy (2mm) was collected from all individuals when they were first trapped. Each individual was weighed, and head (Plate 3-1(a)) and pes (foot) length were measured using vernier callipers. Additional measurements for scrotum length and breadth were also collected from males (Plate 3-1(b)), and the reproductive status was checked for females (Plate 3-2). Once the animals were processed, each individual was returned to their point of capture and released.
(a) (b)
© E. Miller © E. Miller
Plate 3-1 Collecting morphological measurements from a male greater bilby (Macrotis lagotis), (a) head length (mm) and (b) testis length (mm).
© E. Miller Plate 3-2 Inspecting the pouch of a female greater bilby (Macrotis lagotis) that had twins present.
Chapter 3 Genetic mating system, reproductive success and selection 74
3.2.2 Microsatellite genotyping
DNA was successfully extracted from small tissue biopsies from 216 individuals using a salting out method (Sunnucks & Hales 1996; Appendix 1). Each individual was screened at nine polymorphic microsatellite loci (Bil02, Bil16, Bil17, Bil22, Bil41, Bil55, Bil56, Bil63 and Bil66) previously characterised from the greater bilby (Moritz et al. 1997). Genotyping was carried out using multiplex PCR, whereby fluorescently labeled primers enabled the simultaneous amplification of many targets of interest, using several pairs of primers in a single reaction. Two multiplex combinations were used; (i) Bil55, Bil22, Bil16, Bil41, Bil17; and (ii) Bil02, Bil56, Bil63, Bil66. All loci were best amplified in a “touchdown” PCR, whereby the initial annealing temperature was decreased at 1ºC increments and run for a total of ten cycles each. Ten micro-litre reactions were performed using a Multiplex PCR kit (Qiagen, Australia), 0.2μM of each primer and 60 – 80ng genomic DNA. Amplifications were carried out in a PTC- 220 thermocycler (MJ Research, USA) using an initial HotStarTaq activation step at 95ºC for 15 min, a total of 50 cycles of 94ºC for 30 s, 60, 59, 58, 57, 56 and 55ºC for 1 min 30 s (10 cycles each), 72ºC for 1 min 30 s, and final extension at 72ºC for 10 min. The PCR products were analysed in a 48 capillary AB 3730 DNA Analyser (Applied Biosystems, USA). The DNA fragments were sized and quantified using GeneMapper 3.7 (Applied Biosystems, USA).
3.2.3 Parentage analysis
Locus independence and Hardy-Weinberg equilibrium tests were conducted using GenePop 3.4 (Raymond & Rousset 2003) using a Markov chain method (1000 iterations). The statistical significance levels were corrected for multiple comparisons using sequential Bonferroni adjustments (Rice 1989). Using MICROCHECKER, the genotype file was checked for duplicate entries, scoring errors due to null alleles, and stutter bands (van Oosterhout et al. 2004). Null alleles can result in parent-offspring mismatches, therefore loci with a high frequency of null alleles were excluded from paternity analysis.
Paternity and maternity were examined using CERVUS 2.0 (Marshall et al. 1998), a program that assigns paternity using a likelihood based approach. CERVUS calculates
Chapter 3 Genetic mating system, reproductive success and selection 75 a likelihood ratio (LOD) score for each candidate parent and assigns paternity to the most likely individual at a given statistical confidence, taking into account scoring errors, missing data and the proportion of candidate parents sampled. The program conducts simulations based on these input parameters to determine the significance levels of the LOD scores and estimate the expected rate of the successful parentage assignment in the population. Paternity may be assigned at 95% confidence (strict) criterion or 80% confidence (relaxed) criterion (Marshall et al. 1998). Offspring that were assigned at both the 80 and 95% confidence level were included in subsequent analyses.
To be able to assign an individual as a putative parent rather than an offspring, knowledge of individual’s age was required for the parentage analysis. Age class was calculated for all individuals using head length and body weight (Southgate 2005), and date of birth was estimated based on their age class. All individuals that were born in captivity were considered offspring (n = 145). To narrow down the number of candidate parents per offspring, the dataset was divided into annual ‘populations’. Each population represents every independent individual present in the population in a 12 month period. The assigned parentages were cross referenced with the candidate parent’s age and location. Of the total number of individuals relevant to this study (n = 232), 93% were genotyped successfully (n = 216). The proportion of loci sampled was 0.982, and the typing error was estimated as 0.006. Ten thousand iterations were performed for this simulation. Paternity was assigned to individuals with the highest LOD scores at the 95% and 80% statistical confidence levels. The paternity results were examined to assess whether males were successful in siring offspring in multiple years. Maternity was examined for evidence of matings with the same or different males within and between years. New individuals (total n = 57) were translocated into the population between 1999 and 2004. To examine whether females preferred new or ‘novel’ males introduced into the population, the proportion of sires that were immigrants was calculated.
Chapter 3 Genetic mating system, reproductive success and selection 76
3.2.4 Male reproductive success and selection analysis
Data Reduction using a Principal Components Analysis (PCA) was applied to determine the most suitable morphological variables for further analysis. A One-Way ANOVA using Tukey’s Post-Hoc Multiple Comparisons was used to test the relative importance of the various traits for male reproductive success by dividing the males into two subsets; sires and non-sires. A linear regression was used to establish whether male body weight could predict the number of offspring sired. Additionally, we tested for positive assortative mating between males and females with respect to body size. The body weight of successful males was regressed on the body weight of females with whom maternity was assigned. All statistical analyses were carried out using SPSS 14.0.
To assess whether variance in male reproductive success was arising from variance in fitness components (i.e., natural selection) (Arnold & Wade 1984), a selection analysis was performed. The strength of linear and non-linear selection upon male traits was calculated using a multiple regression based method (Lande & Arnold 1983) that accounts for the effects of correlation among traits. Briefly, male reproductive success was standardised so that the mean was equal to one, and the four fitness components, body weight (g), head length (mm), pes length (mm) and testicular volume (mm3), were transformed into Z scores for each male (n = 147). Testicular volume was calculated using the formula for an oblate spheroid, V = ( /6).B2.L, where V = volume, B = breadth of testis and L = length of testis (Williamson et al. 1990). A multiple linear regression was fitted to the standardised traits and measure of reproductive success to calculate the vector of linear selection gradients ( ). A quadratic regression model was then applied to estimate the matrix of nonlinear selection gradients ( ). As nonlinear selection is often underestimated, a canonical analysis was performed to find the major axis of the response surface, the M matrix (Blows & Brooks 2003; Phillips & Arnold 1989). The appropriate gradients have been doubled in this analysis as it has been suggested that quadratic regression coefficients
( ii) obtained from statistical packages (including SPSS 14.0) need to be doubled to obtain the estimated quadratic selection gradients (Lande & Arnold 1983; Stinchcombe et al. 2008). The strength of the nonlinear selection along each of the eigenvectors (mi) of is given by their eigenvalues ( i) using PopTools 3.0.3 (Hood
Chapter 3 Genetic mating system, reproductive success and selection 77
2008). Linear selection gradients of the eigenvectors (mi) of are given by . The major axis of selection was defined as the eigenvectors (mi) that display a significant level of linear and/or nonlinear selection. Selection on the major axis of the response surface was visualised using a spline three-dimensional surface plot in STATISTICA 7.0.
3.3 Results
3.3.1 Parentage assignment
MICROCHECKER found three duplicate individuals that were removed from the dataset. Hardy-Weinberg equilibrium was rejected for Bil66, as there was a significant (p > 0.05) excess of homozygotes and evidence of null alleles. There was no evidence of allele dropout for all loci. One locus (Bil17) failed to amplify successfully. Both Bil66 and Bil17 were excluded from further analyses. All remaining loci were polymorphic with the total number of alleles ranging from three to eight (mean =
6.4 ± 0.6) and the mean heterozygosity (He) estimates varying from 0.327 to 0.810 (mean = 0.709 ± 0.041). These data indicate there is an adequate level of genetic diversity to determine paternity. Allele frequency data are presented in Appendix 2.
Paternity was confidently assigned to 55 individuals at 80% statistical confidence, and 31 (56%) of these at the 95% confidence level. Paternity was not monopolised by a single or small group of males. Females did not show a preference towards immigrant males as there was no discrepancy between the proportion of immigrant males and their representation as sires within the population (19.6%). In any given year, 59.2 ± 9.3% of males in the population did not sire any offspring (Figure 3-1). Of the reproductively successful males, the majority sired one offspring (69.8 ± 3.6%), and 28.3 ± 1.8% sired multiple offspring (two or three) across one or two years, but not necessarily consecutively. One male sired six offspring in 2003. Maternity was confidently assigned to 55 individuals at 80% confidence and 25 (45%) of these at the 95% confidence level. Both paternity and maternity was confidently assigned to 20 offspring with 80% confidence. Only three of the 20 females had paternity assigned for multiple offspring (two each). Females did not repeatedly mate with the male that
Chapter 3 Genetic mating system, reproductive success and selection 78 sired their previous offspring i.e. we did not find any indication for mate fidelity. However, this needs to be interpreted with caution due to the small sample size (n = 6).
3.3.2 Morphological traits and male reproductive success
The PCA indicated that body weight, head length, scrotum length and width accounted for 82% of the variation in male reproductive success and were included in further analysis. Pes length was excluded as it contributed little to the model. Contrary to our predictions, sires and non-sires could not be distinguished from one another based on their morphological traits. There was no significant difference between sires and non- sires in body weight (F4, 143 = 1.135, p = 0.342), head length (F4, 143 = 0.590, p =
0.670), scrotum width (F4, 143 = 0.680, p = 0.607) or scrotum length (F4, 143 = 0.663, p = 0.619). The regression analyses further supported these results as there was no relationship between male body weight and the number of offspring sired in any given 2 year (r = 0.066; F1, 18 = 0.079; p = 0.782; Figure 3-2). In addition to the lack of evidence that morphological traits influence male reproductive success, there was no evidence of non-random (assortative) mating between males and females with respect to body weight. The regression revealed the body weight of females was not significantly correlated with that of the males who fathered their offspring (r2 = 0.003;
F1, 18 = 0.050, p = 0.826; Figure 3-3), suggesting random mating likely occurs in bilbies.
Chapter 3 Genetic mating system, reproductive success and selection 79
100
80
60
40
Percentage (%) of males 20
0 0123456 No. offspring sired
Figure 3-1 Percent (%) males siring offspring (± se) between 2000 and 2004 in a semi free-ranging captive greater bilby (Macrotis lagotis) population.
7
6
5
4
3
No. offspring sired sired offspring No. 2
1
0 0 500 1000 1500 2000 2500 3000 Body weight (g)
Figure 3-2 Relationship between male body weight (g) and the number of offspring sired between 2000 and 2004 in a semi free-ranging captive greater bilby (Macrotis lagotis) population.
Chapter 3 Genetic mating system, reproductive success and selection 80
1500 )
1000 Female weight (g weight Female
500 500 1000 1500 2000 Male weight (g)
Figure 3-3 Relationship between body weight (g) of male greater bilbies (Macrotis lagotis) and the females with which they sired offspring in a semi free-ranging captive population.
3.3.3 Selection analysis
To test whether selection was acting on male traits (summarised in Table 3-1) to enhance reproductive success, male reproductive success was standardised as a measure of fitness in a selection analysis. There was no evidence for overall strong linear or nonlinear selection acting on the male morphological traits measured in this study. This result was qualitatively the same, regardless of whether canonical rotation was applied to the data (Tables 3-2 and 3-3). However, there was no significant trend for non-linear selection on some aspects of the data. The canonical rotation of the matrix by symmetric Eigenanalysis of returned four eigenvectors (mi) representing the major axis of the response surface. The M matrix and their associated eigenvectors
(Table 3-3) show eigenvectors m1 and m4 to be under significant nonlinear selection.
Body weight and pes length had the heaviest loading for m1, and pes length had for m4.
The major axis of nonlinear selection, m1 and m4, indicates that there does not appear to be a single optimum phenotype associated with fitness (reproductive success) in
Chapter 3 Genetic mating system, reproductive success and selection 81 male bilbies. Males of various weights, with varying pes length had high fitness (Figure 3-4). Nevertheless, heavier males with a medium pes had higher reproductive success in this study. Body weight and pes length are correlated traits, and although the selection analysis is designed to cope with correlational data, this may explain why the overall selection was not significant.
Table 3-1 Male morphological traits (mean ± se) measured in the greater bilby (Macrotis lagotis) to examine the strength of selection.
Trait n Mean se Range Body weight (g) 147 1213.96 33.47 499 – 2483 Head length (mm) 147 104.88 0.82 85 – 138 Pes length (mm) 147 97.17 2.68 79 – 111 Testicular volume (mm3) 147 30.05 0.10 3 – 17
Table 3-2 Standardised linear gradients ( ) and matrix of quadratic and correlational selection gradients ( ) in the greater bilby (Macrotis lagotis). The multiple regressions describing selection were not significant (linear F4, 143 = 0.034, p = 0.875; nonlinear F4,
143 = 1.044, p = 0.387).