Non-invasive Genetic Hair Sampling of a Population of
←
→
Page content transcription
If your browser does not render page correctly, please read the page content below
Non-invasive Genetic Hair Sampling of a Population of Bridled Nailtail Wallaby (Onychogalea fraenata) on Avocet Nature Refuge Leanne Henry (Simpson) - 220055558 Master of Scientific Studies – Thesis SCI695 University of New England Supervisor: Stuart Cairns Co-supervisor: Jane Hughes (Griffith University)
Non-invasive Genetic Hair Sampling of a Population of Bridled Nailtail Wallaby (Onychogalea fraenata) on Avocet Nature Refuge Leanne Henry (Simpson) - 220055558 Master of Scientific Studies – Thesis SCI695 University of New England Supervisor: Stuart Cairns Co-supervisor: Jane Hughes (Griffith University)
Acknowledgements I would first like to thank the Spooner family for acknowledging and embracing the natural values on their cattle property near Springsure and providing a home for one Australia’s most endangered species. Hugo Spooner provided endless assistance and support throughout the project and for that I am very grateful. Many thanks must also go to my primary supervisor Stuart Cairns for all of his guidance, analytical and statistical knowledge, as well as field support. To my co-supervisor Jane Hughes and research assistant Kathryn Real from Griffith University, I cannot thank you enough for your direction, assistance and patience with all things genetic. To the Norman Wettenhall Foundation, the Save the Bilby Fund and the Dept. of National Parks, Recreation, Sport and Racing's, Queensland Parks and Wildlife Service, all my thanks go to them for their generous financial contributions and their tremendous support of this project. A scientific purposes permit (WITK09715411) was granted through the Department of Environment and Resource Management (now Dept. of Environment and Heritage Protection) and ethics approval was granted from the University of New England Ethics Approval Committee (AEC11/053). To Janelle Lowry from DEHP, with her dedication for all things nailtail and tremendous support in all aspects of the project, I am unable to thank you enough for your advice and endless assistance. Also to Alan Horsup, from DEHP who I would like to sincerely thank for sharing his knowledge of hair sampling and for all of his support and assistance throughout the project. Many thanks go to Tina Janssen and her family for their hospitality and assistance in the initial stages of this project. This project received an enormous amount of support from organisations such as the Fitzroy Basin Association, in particular Graeme Armstrong who provided an enormous amount of assistance and direction, and Lyndal Rolfe, the Bridled Nailtail Trust, namely Fiachra Kearney, the Conservation and Wildlife Management Branch of the SSAA, especially Mark Woods and Glenys-Julie Harris, Murray Haseler from Bush Heritage and Wild Mob Volunteers Australia. I would like to sincerely thank them all for providing a significant amount of both advice and field support which was boundless throughout the project. This project would not have been possible without the help of the many volunteers who gave their support and, at times, a significant amount of sweat to help me complete the field component of this project. Many thanks go to Janelle Lowry, Alan Horsup, Bennett Henry, Colleen and Ken King, Lorelle Campbell, Graham Armstrong, Lyndy Marshall, Lyndal Rolfe, Stuart Cairns, Gerhard Koertner, Shellie Cash, Mark Woods and Glenys-Julie Harris. i
Many people from the QPWS provided a substantial amount of support and reassurance throughout the project. In particular, I would like to thank Andrew Dinwoodie, Barry Nolan, John Augusteyn, Graham Hemson, Rhonda Melzer, Peter Moore and Leigh Harris and special thanks goes to my work supervisor, Scott Brook for his support and patience over the last 12 months. Thanks also goes to all of the wonderful people from the Griffith University Molecular Ecology laboratory, for providing endless assistance and opening my eyes to the world of genetics; in particular Kathryn Real, Dan Schmidt and Joel Huey. Many thanks go to Susan Nuske from the University of Queensland for her guidance and advice and to Lauren Young and Kerry Walsh from Central Queensland University for their assistance in the initial stages and during the analysis stage of this project. To my very supportive husband-to-be, Bennett Henry I cannot thank you enough for your patience, hard work and endless support throughout the project. To all of my friends and colleagues who offered advice, kept me sane and gave me an excuse to take a break, I whole heartedly thank you. Videos of this hair sampling study on Avocet NR can be viewed on the YouTube channel ‘Nailtail1’. http://www.youtube.com/user/Nailtail1 ii
Abstract Sampling wild animal populations is of fundamental importance particularly when managing rare and threatened species. This study investigates the use of non-invasive genetic hair sampling to census a reintroduced population of bridled nailtail wallaby (Onychogalea fraenata) on Avocet Nature Refuge (NR), and compares the technique to traditional cage trap sampling. Cage trapping was conducted using 146 metal traps at Avocet NR in June 2011 over four nights. Hairs for genetic sampling were collected from 180 hair traps (78 runway and 102 triangle) on Avocet NR in October 2011 over 18 nights. Bridled nailtail wallaby hair was macroscopically identified before DNA was extracted with chelex for genotyping. Tissue samples from cage trapped individuals and hair samples were extracted and genotyped on five highly polymorphic microsatellite loci. Error rates were calculated and incorporated into the results as extracted hair samples produce low quality and small quantities of DNA. Cage trapping resulted in a total of 20 captures and the identification of 10 individual bridled nailtail wallabies. Hair trapping resulted in a total of 422 captures and identified 66 individual wallabies. The two hair trapping methods used, showed no significant difference in their ability to capture a hair sample. With the amplification of only 71 samples at four or more loci, only 66 different bridled nailtail wallabies were able to be confidently identified (PI 0.999). The poor amplification of loci was most likely due to the large fragment length (>150bps) of the loci B29, B87 and B151. Using two different methods, the effective population size (Ne) of bridled nailtail wallabies was calculated at 29 (18-50) and 35 (24-61). The initial cost to set up a cage trapping survey at Avocet NR is $70,676 when compared to $46,613 for a hair trap survey. Ongoing surveys for cage trapping cost $28,546 and for hair trapping the cost is $40,857, this includes $12,000 for genetic analysis. Hair trapping is much more effective in capturing bridled nailtail wallaby samples from the population than the currently used cage trapping method. A population estimate based on the capture-mark-recapture (CMR) model was unable to be conducted for either trap survey method due to the small sample size of identifiable individuals and the low number of recaptures. Using the current technique, cage trapping is not suitable for estimating population size at Avocet NR. To improve sample size for non-invasive hair samples, it is recommended to investigate different extraction methods and establish new highly variable markers with a fragment length
Table of Contents Acknowledgements ....................................................................................................................................i Abstract.................................................................................................................................................... iii List of Figures ........................................................................................................................................... v List of Tables ............................................................................................................................................ v 1. Introduction ......................................................................................................................................... 1 1.1 Bridled nailtail wallaby .................................................................................................................... 2 1.2 Non-invasive genetic hair sampling ................................................................................................. 4 2. Methods ............................................................................................................................................... 7 2.1 Study site ........................................................................................................................................ 7 2.2 Cage trapping.................................................................................................................................. 7 2.3 Hair trapping .................................................................................................................................. 8 2.4 Hair analysis.................................................................................................................................. 10 2.5 Genetic analysis ............................................................................................................................ 11 2.5.1 DNA Extraction .................................................................................................................... 11 2.5.2 Microsatellite genotyping........................................................................................................ 11 2.5.3 Statistical analysis ................................................................................................................... 12 2.6 Cost effectiveness ......................................................................................................................... 14 3. Results ................................................................................................................................................ 15 3.1 Cage Trapping .............................................................................................................................. 15 3.2 Hair trapping ................................................................................................................................ 15 3.3 Genetic analysis ............................................................................................................................ 17 3.4 Cost comparison ........................................................................................................................... 19 4. Discussion .......................................................................................................................................... 21 4.1 General observations .................................................................................................................... 21 4.2 Trap method ................................................................................................................................. 21 4.2.1 Animal ethics ......................................................................................................................... 24 4.3 Genetic Analysis ........................................................................................................................... 25 4.4 Cost effectiveness ......................................................................................................................... 28 5. Recommendations .............................................................................................................................. 31 5.1 Hair trapping ................................................................................................................................ 31 5.2 Genetic analysis ............................................................................................................................ 31 References .............................................................................................................................................. 34 Appendices............................................................................................................................................. 38 iv
List of Figures Fig. 1: Bridled nailtail wallaby (Onychogalea fraenata) 3 Fig. 2: Former and current distribution of the bridled nailtail wallaby 4 Fig. 3: Metal cage trap used to capture bridled nailtail wallabies 8 Fig. 4: Two types of hair traps used in this study: triangle hair trap and runway hair trap 9 Fig. 5: The number of bridled nailtail wallaby captures using cage traps 15 Fig. 6: The total number of hair trap sites that collected bridled nailtail wallaby hair 16 Fig. 7: The average number of individual samples collected from each runway and 17 triangle trap Fig. 8: The COLONY output showing the relationships between the individuals 19 based on genetic data Fig. 9: Cost associated with the initial setup and ongoing survey for hair and cage trap 20 methods List of Tables Table 1: Genetic analysis of hair samples including error rates 18 v
1. Introduction Sampling of wild populations is of primary importance in the management of many rare and threatened species across the world. Estimating the population size of a species enables managers to monitor for positive or negative trends that may result from particular management actions. One of the more commonly used population size estimation methods is known as capture-mark-recapture (CMR) (Attiwill & Wilson 2006). The aim of the CMR method is to capture individuals from the population, mark them and release them back into the population so they can be recaptured on subsequent occasions (Amstrup et al. 2005). This method is commonly used for estimating the population size of many small to medium-sized mammal species (Minta et al. 1989; Nichols 1992; Attiwill & Wilson 2006; Augusteyn et al. 2010; Ruibal et al. 2010). The accuracy of the CMR method relies on all individuals having an equal chance of being captured and then recaptured on subsequent sampling occasions (White et al. 1982; Amstrup et al. 2005). Using cage or foot-hold traps are effective ways to capture a proportion of the population so long as the species is abundant with high tradability (Attiwill & Wilson 2006). Invasive sampling using traps can, however, change the behaviour of individuals. Attractants (food rewards) may make some individuals ‘trap happy’, while the experience of being trapped and then handled can make others ‘trap shy’ (Attiwill & Wilson 2006). Live-trap types can also impact directly on the health of the individual or their young by causing injury or death, or by increasing the chance of predation (Johnson 1997; Lemckert et al. 2006). In more recent times, non-invasive genetic sampling methods have been developed to better estimate the population size, sex ratio, distribution and genetic variability of species where invasive or other non-invasive methods such as the mark-resight method are not feasible. Traditional capture methods such as cage or foot-hold trapping have proved ineffective in the population size estimation of several species; including bears (Mowat & Strobeck 2000; Romain- Bondi et al. 2004), otters (Depue & Ben-David 2007) and wombats (Hoyle et al. 1995). Cage traps were used as the capture method to estimate the population size of free-ranging bridled nailtail wallabies (Onychogalea fraenata) from 1993-1996 (Fisher et al. 2000) and since 2007 (Augusteyn et al. 2010). When applied to CMR models, this capture method produces inconsistencies when estimating population size. When two trapping events from the same year were compared, the known minimum population size based on marked individuals was higher than the population estimate using the CMR model (Augusteyn et al. 2010). Mark-resight data from distance estimation methods was used prior to 2002 but became unreliable owing to reduced visibility caused by encroaching introduced grasses (Lundie-Jenkins & Lowry 2005a). Non-invasive Genetic Hair Sampling Page 1
The use of genetic material in association with wildlife management provides opportunities for assessing the genetic variability and structure within a population, and how it interacts with other populations (immigration/emigration). It also provides information on how individuals interact with each other (breeding patterns), what areas of habitat they use, and enables the best combinations of individuals to be identified for use in breeding and translocation programs (Lindenmayer & Burgman 2005; Attiwill & Wilson 2006). This information assists managers to manage rare and threatened species in the wild or in captivity, more effectively. Genetic structure is directly related to species’ survival as the greater the variability within a population the greater the chance that individuals will be equipped to deal with stochastic events (Attiwill & Wilson 2006). Traditional methods of obtaining DNA from animals involve tissue or blood samples being collected from confined individuals (Henry & Russello 2011). With improvements in DNA extraction and polymerase chain reaction (PCR) methods (Taberlet et al. 1996), as well as improvements in the development of analytical methods and user-friendly analysis software, very small quantities of DNA can now be utilised to genotype individuals (Henry & Russello 2011). These advances have paved the way for the use of non-invasively collected genetic material from hair, scats, feathers and urine, to determine the size and structure of wild animal populations (Dreher et al. 2009; Henry & Russello 2011). This study examines the use of a non-invasive hair sampling method on a small population of bridled nailtail wallabies (Onychogalea fraenata) on Avocet Nature Refuge (NR) in central Queensland. The study aimed to compare two types of hair trap, to genetically identify individuals from remotely collected hair samples, and from that data, to describe the genetic structure and estimate the size of the population. 1.1 Bridled nailtail wallaby The bridled nailtail wallaby (O. fraenata) is one of Australia’s medium sized macropods most at risk of extinction due to habitat loss, competition with domestic livestock and predation from introduced animals (Lundie-Jenkins & Lowry 2005a; Van Dyck & Strahan 2008) (Fig. 1). This species is listed as ‘Endangered’ under Commonwealth legislation, the Environmental Protection and Biodiversity Conservation Act 1999, and Queensland legislation, the Nature Conservation Act 1992. The bridled nailtail wallaby is one of three related Australian wallabies possessing a horny spur at the end of its tail; the other two being the crescent nailtail wallaby (O. lunata) which is declared extinct and the northern nailtail wallaby (O. unguifera) which is common across northern Australia (Van Dyck & Strahan 2008). Non-invasive Genetic Hair Sampling Page 2
Fig. 1. Bridled nailtail wallaby (Onychogalea fraenata) on Avocet Nature Refuge in central Queensland. The bridled nailtail wallaby weighs between 5-8 kg for males and 4-6 kg for females; reaching sexual maturity after only 270 days and 136 days, respectively (Hendrikz & Johnson 1999; Lundie-Jenkins & Lowry 2005a). Females exhibit embryonic diapause and have comparatively high fecundity with the ability to give birth to three young per year under favourable conditions (Johnson 1997). The bridled nailtail wallaby lives for 6-8 years in the wild and were previously known to inhabit areas of acacia shrubland and grassy woodland within the Brigalow Belt bioregion of semi-arid eastern Australia (Van Dyck & Strahan 2008) (Fig. 2). This medium-sized wallaby was first described in 1840 by John Gould (Flannery et al. 1990). The bridled nailtail wallaby received its name from the distinctive white ‘bridle’ running from the centre of the neck down behind the forearm and the horny spur on the tip of its tail (Van Dyck & Strahan 2008). In the 1800s the bridled nailtail wallaby was widespread through the semi-arid areas of central Queensland to northern Victoria (Gordon & Lawrie 1980). Its numbers had declined in most areas by the early 1900s, with the last recorded sighting in 1937 (Pople et al. 2001). The bridled nailtail wallaby was presumed extinct until being rediscovered on two properties in central Queensland near Dingo in 1973 (Gordon & Lawrie 1980; Pople et al. 2001) where they inhabit brigalow (Acacia harpophylla), poplar box (Eucalyptus populnea) and rosewood (Acacia rhodoxylon) communities (Fisher 2000). Non-invasive Genetic Hair Sampling Page 3
Fig. 2. Former and current distributions of the bridled nailtail wallaby throughout eastern Australia (modified from Van Dyck & Strahan 2008). Since the rediscovery of the bridled nailtail wallaby, captive breeding programs have been developed and translocations have taken place, to establish new populations (Lundie-Jenkins & Lowry 2005a). There are currently three free-range populations in Queensland; the remnant population at Taunton National Park (Scientific) (NP) and the translocated populations at Idalia National Park (NP), and Avocet Nature Refuge (NR) (DERM 2011a). There is also another population within an enclosure at the Australian Wildlife Conservancy, Scotia Sanctuary in south- west New South Wales with trial translocations occurring outside of the enclosure since 2010 (AWC 2012) (Fig. 2). Despite every effort to preserve the bridled nailtail wallaby, it appears that the Taunton NP population has suffered a substantial decline from approximately 1400 individuals in 1991 (Lundie-Jenkins & Lowry 2005a), to possibly fewer than 150 in 2011 (Augusteyn et al. 2010; DERM 2011a). CMR population estimates using cage traps show the Idalia NP population to also be fewer than 150, while the Avocet NR population was estimated at fewer than 50 individuals in 2008 (DERM 2011a; Kingsley et al. 2012). Several genetic studies have been conducted on bridled nailtail wallaby populations in an attempt to better understand this species and enable managers to improve genetic diversity in the translocated populations through captive breeding (Sigg 2004; Lundie-Jenkins & Lowry 2005b). Genetic analysis of stored tissue samples was undertaken in 2008 on all populations of bridled nailtail wallabies in an attempt to select additional individuals that could contribute new alleles to the captive breeding program and ultimately populations at the translocation sites (Seddon 2008). Non-invasive Genetic Hair Sampling Page 4
1.2 Non-invasive genetic hair sampling Non-invasive sampling using collected genetic material has been used successfully to estimate population size, occupancy, abundance, sex ratio and genetic diversity (Woods et al. 1999; Bremner-Harrison et al. 2006; Downey et al. 2007). This method has been applied to numerous mammal species, including felids (Downey et al. 2007), bears (Woods et al. 1999; Mowat & Strobeck 2000; Romain-Bondi et al. 2004; Gervasi et al. 2008), wombats (Sloane et al. 2000; Banks et al. 2003; Walker et al. 2006), quolls (Ruibal et al. 2010) and otters (Depue & Ben-David 2007; Hajkova et al. 2009). This form of population monitoring suits species that have large home ranges, small population sizes or are cryptic in nature (Taberlet & Luikart 1999; Downey et al. 2007). The use of conventional CMR sampling has proved unreliable for several species because they are restricted by low sample size and low detection probabilities which violates CMR assumptions (Williams et al. 2009). Non-invasive genetic sampling involves the collection of a DNA sample, which is usually obtained from either hair, feathers, urine or scats. Hair has been collected from wombats for DNA analysis using double-sided tape across burrow entrances (Sloane et al. 2000) and from bears collected from a barbed wire enclosure (Woods et al. 1999). In order to identify individuals from these samples, DNA is extracted and then genotyped (Piggott & Taylor 2003; Waits & Paetkau 2005; Attiwill & Wilson 2006). A pilot study was conducted in 2010 on collecting and genotyping hair from bridled naitlatil wallabies, outlining the best hair collection methods and error rates from DNA analysis (Nuske 2010). This form of non-invasive sampling enables researchers to obtain DNA from an individual with minimal to no interaction with the species involved (Henry & Russello 2011). This improves the power of CMR models by increasing the likelihood of obtaining information from cryptic species as well as ensuring that samples conform to CMR assumptions (Bremner- Harrison et al. 2006; Downey et al. 2007; Dreher et al. 2009). Based on the genetic information gathered from hair sampling an effective population size is able to be calculated. The effective population size (Ne) is a measure of genetic diversity, being the number of individuals that have equal opportunity to contribute genes to the next generation (Lindenmayer & Burgman 2005). Knowledge of Ne is particularly important when managing small populations of isolated species since genetic diversity is, in most cases, directly related to their ability to persist (Lindenmayer & Burgman 2005; Janecka et al. 2008). Non-invasive Genetic Hair Sampling Page 5
The estimation of the effective and actual population size is crucial for establishing a goal set for endangered species management (Williams et al. 2009). The drawback of using non- invasive genetic methods is the prevalence of genotyping errors that may result in the over or under estimation of a population (Creel et al. 2003). Several methods have been established to account for such errors and conservative methods are recommended when genotyping low quality and small quantities of DNA as found in hair and scat samples (Taberlet & Luikart 1999; Pompanon et al. 2005). These methods include repeat PCRs, calculating error rates based on probability of identity (PI) (Pompanon et al. 2005), using highly polymorphic loci and the use of dinucleotide repeat motifs and shorter fragment lengths
2. Methods 2.1 Study site Avocet Nature Refuge (NR) is an 1140 ha portion of a 4500 ha cattle grazing property located approximately 37 km south of Emerald and 35 km north of Springsure in central Queensland (S23o51’12.5” E148o09’54.5”) (H. Spooner, pers. com. 2012) (Fig. 2). Rainfall for this site was gathered from Hugo Spooner from the Avocet house approximately 4.5 km north-west of Avocet NR. A Nature Refuge in Queensland is established and gazetted by a legally binding, perpetual agreement attached to the land title under the Nature Conservation Act 1992 (McCosker 1998; DERM 2011b). Avocet NR was gazetted at the request of the land holder, Hugo Spooner, in 1998 to protect an area of endangered native brigalow woodland. Avocet NR receives on average 625 mm of rain annually and contains a mixture of Brigalow, poplar box, silver-leafed ironbark (Eucalyptus melanophloia), narrow-leafed ironbark (Eucalyptus crebra) and lancewood (Acacia shirleyi) communities, located within the brigalow belt bioregion (McCosker 1998; DERM 2012). Avocet NR adjoins the Goonderoo Bush Heritage Reserve which was gazetted a Nature Refuge in 2003 to also protect natural brigalow woodland communities (BHA 2012) (Appendix 1). Avocet NR was recognised by the Environmental Protection Agency (EPA) (now Department of Environment and Heritage Protection) in 1998 as a suitable site for the reintroduction of a population of bridled nailtail wallabies (J. Lowry, pers. com. 2012). With support and assistance from the Spooner family, as well as Emerald Shire Council (now part of the Central Highlands Regional Council) and groups such as the Sporting Shooters Association of Australia (SSAA), Bridled Nailtail Wallaby Trust and the local community, a total of 166 bridled nailtail wallabies were introduced over 14 releases, onto Avocet NR between 2001 to 2005 (MacLeod et al. 2010; Kingsley et al. 2012). Since then, other partners such as the Fitzroy Basin Association and Minerva Mine have contributed to the recovery of this species. 2.2 Cage trapping Cage trapping of bridled nailtail wallabies was undertaken on Avocet NR by the Queensland Parks and Wildlife Service (QPWS) from the 14-18 June 2011. In the three months preceding cage trapping, the Avocet area received a total of 304 mm of rain (H. Spooner, pers. com. 2012). Since this area is a cattle grazing property and the attractant used for bridled nailtail wallabies is also attractive to cattle, the cattle were removed from the cage trapping area prior to the commencement of the survey. Trapping consisted of setting 146 (380 x 380 x 760 mm) treadle triggered, wire cage traps lined with marine carpet (Fig. 3) and placed between 100-150 m apart Non-invasive Genetic Hair Sampling Page 7
(Appendix 2). Cage traps were brought from Taunton NP and laid out by a team of two people over two days. This included the time it took to conduct minor cage trap repairs (J. Lowry, pers. com. 2011). Prior to trapping, traps were pre-fed with premium lucerne hay, in a ‘locked open’ position, on three occasions over a period of one week. During trapping, traps were set about an hour before dusk and checked twice each night for four consecutive nights (eight occasions). Traps were then closed at dawn to prevent capture during the day (J. Lowry, pers. com. 2011). Any lucerne outside of the traps was removed at each visit to increase the chance of bridled nailtail wallabies entering the traps and triggering the treadle. Each trapped bridled nailtail wallaby was placed in a hessian bag and processed. This involved taking a 3-mm diameter sterile biopsy punch of tissue from the ear for DNA analysis and inserting a passive implant transponder (microchip) between the scapular for ‘first-time captures’ and recording standard body measurements (Jackson 2003; Nolan 2012). See Appendix 3 for datasheet used. Once this was completed, the animals were immediately released at the site of capture. Collection of the cage traps post-trapping involves extensive washing and disinfecting (Virkon) of the traps to prevent the spread of disease and parasites (Chapman et al. 2005; Nolan 2012). Fig 3. Metal wire cage trap lined with carpet used in the capture of bridled nailtail wallabies, baited with lucerne. Non-invasive Genetic Hair Sampling Page 8
2.3 Hair trapping Two types of hair traps were used to obtain hair samples from bridled nailtail wallabies: triangle and runway traps (Fig. 4). The types of hair traps used in this study were based on a pilot study of non-invasive hair trapping of bridled nailtail wallabies conducted at Taunton NP (Nuske 2010). Triangle hair traps consisted of 8-mm diameter rope tied around three to four trees or metal stakes in the shape of a triangle. Double-sided tape (4970 TESA tape) was attached to the undersides of the two to three, 60-80 cm long, sections of the triangle. Runway hair traps consisted of two stakes or trees located either side of an animal pad about 30-40 cm apart with double-sided tape strung between them. All hair traps were set at a height of 20-25 cm. A triangle or runway trap was selected depending on their appropriateness to the site (i.e. the presence of three small trees close together or well used runways). A total of 180 hair traps (78 runways and 102 triangles) were set along 30 km of track on Avocet NR (Appendix 4). Hair traps were placed every 50-100 m in areas of previous bridled nailtail wallaby captures or sightings and every 100-300 m apart in areas where bridled nailtail wallabies had not previously been sighted or where the habitat was deemed to be of less suitable for the bridled nailtail wallaby (Fisher 2000). Hair traps were installed during the period from the end of September to the beginning of October 2011. Selection of sites and installation of ropes and stakes took three full days for 4-5 people. A small amount of attractant (premium lucerne and sweet potato) was placed in the centre of each triangle and on both sides and through each runway (Fig. 4). During pre-feeding, Fig. 4. The two trap types used for hair sampling: triangle hair trap with attractant (left) and a runway hair trap with attractant (right). Non-invasive Genetic Hair Sampling Page 9
hair traps had a 5 cm piece of double-sided tape attached to allow bridled nailtail wallabies to habituate to the smell of the tape. The traps were baited twice over a one-week period before all trap sites were set with hair tape. TESA tape was used as it has proven to be very effective in the capture of wombat hair (Banks et al. 2003; Walker et al. 2006). Hair traps were set for a total of 18 consecutive nights during the period 7-25 October 2011. All traps were checked on four occasions in this period, every 3-5 days. All hair traps with hair, had that section of tape carefully removed. The sample was placed flat on grease-proof paper (Glad Bake), folded and stored in a labelled envelope. Fresh tape was placed across the section removed and fresh attractant was placed at the site. Fresh attractant was also placed at the hair traps which had collected no hair. All hair traps were cleared and removed by the end of October 2011. Envelopes containing hair samples were kept dry in a refrigerator on site (Avocet NR) until they could be transported to the QPWS laboratory in Rockhampton. The samples were then kept dry in a container at room temperature until the DNA extraction could be completed. 2.4 Hair analysis Collected hairs were distinguished from those of other species by visual identification of hair colour, pattern and thickness based on the Hair Identification Key for Australian Mammals (Triggs & Brunner 2002). Using the species list from the Avocet NR evaluation document (McCosker 1998), a table was created to assist with quick on-site identification of hair (Appendix 5). In order to assist in the confirmation of hair identification, ten motion-sensing infra-red (IR) cameras (Bushnell Trail Cameras) were set up at hair traps prior to, and during the survey where attractant was being taken. Hair samples were identified by their appearance and the length and diameter of guard hairs. Guard hairs are the larger or coarser hairs that form the main part of the coat (Triggs & Brunner 2002). Bridled nailtail wallaby guard hairs were
2.5 Genetic analysis 2.5.1 DNA Extraction Tissue samples were collected using 3-mm diameter sterile ear biopsy punches during the CMR cage trapping survey on Avocet NR in June 2011. They were placed in 1.5 ml Eppendorf tubes containing 90% ethanol and stored in a refrigerator at 4oC before being transported on ice to the Griffith University Molecular Ecology laboratory. Half of the tissue was digested with 200 µl of CTAB buffer and 8 µl of proteinase-K and heated to 55oC for approximately five hours. DNA extraction was carried out using a modified CTAB extraction method (Doyle and Doyle 1987). The DNA extracts were then stored in a refrigerator at 4oC. A study conducted in 2010 on the Taunton NP bridled nailtail population confirmed that DNA from a single bridled nailtail wallaby hair could be successfully genotyped when compared to DNA extracted from tissue samples (Nuske 2010). In the laboratory, hairs with a clear follicle (white skin at the base of the hair) were selected for DNA extraction. Given that more than one individual may have passed under the same piece of tape during the same sampling period, the use of pooled hair samples was not considered appropriate (Sloane et al. 2000; Roon et al. 2005; Walker et al. 2006; Depue & Ben-David 2007). Approximately 10 mm was cut from each hair and placed follicle end down in a 1.5 ml Eppendorf tube containing either 200 µl or 100 µl of 5% chelex solution, as described in Sloane et al. (2000). Scissors and forceps were sterilised between the handling of individual hairs. Three replicate individual hairs were collected from each sample (clump of hair on the tape). The second or third replicate was genotyped only in the event that the first or second hair failed to amplify. Follicles in chelex solution were boiled for 10 min then stored in a freezer at -20oC. Extractions were performed in the QPWS laboratory 1-3 months after hair collection. 2.5.2 Microsatellite genotyping Tissue and hair samples were transported to the Griffith University Molecular Ecology laboratory and stored at 4oC prior to genotyping. Five highly polymorphic microsatellite loci (B90, B87, B29, B151, and YM148) that have previously been used to identify bridled nailtail wallabies were used for the microsatellite genotyping of extracted hair samples (Pope et al. 1996; Sigg et al. 2005). Polymerase Chain Reactions (PCRs) were conducted for both the tissue and hair samples. The PCR for extracted tissue samples contained a 10 µl total solution of 0.5 µl of DNA extract, 1x reaction buffer, 1.5 mM of MgCl2, 0.1 µM of forward primer, 0.4 µM of reverse primer, 0.4 µM of fluorescent tag primer (Real et al. 2009), 0.2 mM of dNTPs and 0.3 u of Red Taq (Astral). Non-invasive Genetic Hair Sampling Page 11
PCRs for hair samples extracted in 200 µl of chelex contained a total solution of 12.5 µl consisting of 8.825 µl of DNA extract, 1x reaction buffer, 1.5 mM of MgCl2, 0.1 µM of forward primer, 0.4 µM of reverse primer, 0.4 µM of fluorescent tag primer (Real et al. 2009), 0.2 mM of dNTPs and 0.3 u of Red Taq (Astral). Hair samples extracted in 100 µl of chelex contained at total solution of 10 µl consisting of 4 µl of DNA extract, 1x reaction buffer, 1.5 mM of MgCl2, 0.1 µM of forward primer, 0.4 µM of reverse primer, 0.4 µM of fluorescent tag primer (Real et al. 2009), 0.2 mM of dNTPs and 0.3 u of Red Taq (Astral). Separate PCRs were conducted for each locus. Tissue and hair sample amplifications for B90 and YM148 microsatellite loci were carried out at an initial denaturing of 94oC for 5 min before 35 cycles of 94oC for 30 s of denaturing, 55oC for 30 s of annealing and 72oC for 30 s of extension before a final 72oC for 10 min of extension. Tissue sample PCR conditions for B29, B87 and B151 were the same as for B90 and YM148 except for a different annealing temperature of 58oC. Owing to the poor amplification of B29, B87 and B151 for hair extracts, PCR conditions were altered for these loci to an initial denaturation of 94oC for 5 min, before 40 cycles of 94oC for 30 s of denaturing, 58oC for 90 s of annealing and 72oC for 90 s of extension before a final 72oC for 10 min. All loci were pooled for each tissue sample, while for hair samples they were combined into two pools for genotyping comprising of YM148, B90 and B29, and B87 and B151 to reduce any further dilution of the amplified fragments. The microsatellite fragments were run on an ABI 3130 genetic analyser. Scoring of microsatellite genotypes was undertaken using Genemapper 4.0 (ABI software). 2.5.3 Statistical analysis The effectiveness of the two types of hair traps was compared using a contingency table analysis using the Chi-square statistic (Van Emden 2008). The success of the two types of traps in relation to the number of samples collected from each was compared using a single factor analysis of variance (ANOVA). A correction factor was calculated for analysing the success of the trap types (runways and triangles) in collecting samples as the total number of each trap type was not equal across the survey area (102 triangles and 78 runways). Corrected No. = 50% x No. successful traps of successful traps Actual % The size range of each locus was obtained from Sigg (2004). In a genetic data study conducted on bridled nailtail wallabies in 2008, the Avocet NR population contained an average of 7.4 alleles (Na) per locus based on the five polymorphic loci used in this study (Seddon 2008). Additional individuals were introduced after the collection of the samples used to calculate the 2008 Na. In this study microsatellite data were analysed with the MS Excel-based program Non-invasive Genetic Hair Sampling Page 12
GenAlEx 6 (Peakall & Smouse 2006) to determine the number of alleles (Na) at each locus and the number of matches to distinguish unique genotypes at two and 4-5 loci. The program Micro- checker was used to locate and explain any inconsistencies at each locus and also outlined possible reasons for loci being out of Hardy-Weinberg equilibrium (HWE) (Van Oosterhout et al. 2004). A Hardy-Weinberg Exact Test was conducted on samples that amplified at five loci using the program GENEPOP (Raymond & Rousset 1995; Rousset 2008) to establish if each locus was in HWE, and to also evaluate the inbreeding coefficient (FIS). HWE is based on the main assumptions that the population is large, mating randomly and there is no selection taking place (Ridley 2004). The genotypic frequencies of a natural population typically conform to Hardy- Weinberg expectations (Attiwill & Wilson 2006). If a population is out of HWE, it suggests that something, such as non-random mating, may be occurring in the population (Ridley 2004). Inbreeding occurs when naturally random-mating populations are so small that mating with close relatives results in offspring receiving identical copies of a gene from both parents. This produces a high proportion of homozygotes within the population (Ridley 2004; Attiwill & Wilson 2006). The inbreeding coefficient for an individual relative to the subpopulation (FIS), is the proportion of genes in an individual that are identical because they are derived from a common ancestor (Lindenmayer & Burgman 2005). Samples containing low quality and small quantities of DNA, such as in this study, can be prone to genotyping errors. These are primarily ‘allelic dropout’, which results in a heterozygote being scored as a homozygote, or ‘false alleles’, where an additional allele is produced to show three alleles or a false heterozygote (Taberlet et al. 1996; Pompanon et al. 2005). In order to reduce the likelihood of the occurrence of these errors, any genotypes that differed at only one locus by the scoring of a heterozygote, were assumed to be the same genotype (Taberlet & Luikart 1999). Duplicates of the same samples were genotyped and compared to produce a mean genotyping error rate at each locus (Pompanon et al. 2005). The successful amplification rate was based on the number of samples, out of the 422 collected that amplified at each locus. For individuals that could be genotyped at four or more loci, the program Cervus 3.0 (Kalinowski et al. 2007) was used to establish the exclusion probability (PI) for individual and sibling identity. This is effectively the probability that two individuals or siblings will not share the same multi-locus genotype by chance. The program Cervus 3.0 is able to accommodate simple genotyping errors such as null alleles when assigning parentage and exclusion probabilities (Kalinowski et al. 2007). The aim of selecting primers should be to ensure that the multi-locus Non-invasive Genetic Hair Sampling Page 13
genotype has an exclusion probability (PI) as close as possible to one; this meaning that there is a zero probability that they will share the same genotype by chance. Two programs were used to estimate the effective population size (Ne) based on the unique genotypes that amplified at four or more loci. COLONY which uses “the maximum likelihood method to assign individuals in a sample into full-sib families nested within half-sib families (colonies) using these individuals' multi-locus genotypes without parental information” (Jones & Wang 2009). ONeSAMP which uses “summary statistics and approximate Bayesian computation to estimate Ne from a single sample of microsatellite data” (Tallmon et al. 2008). The effective population size (Ne) is the number of individuals that would experience a loss of genetic variability at the same rate as a randomly mating hypothetical population (Lindenmayer & Burgman 2005). Different programs, such as those that assess effective population size (Ne) and probability of identity have the ability to account for different levels of missing data. For this reason, analysis was able to be carried out for samples that amplified on four loci, together with those that amplified on all five loci (Kalinowski et al. 2007; Tallmon et al. 2008; Jones & Wang 2009). 2.6 Cost effectiveness The comparison of the cost effectiveness of cage trapping versus hair trapping takes into account set up costs for establishing and conducting a new survey, as well as ongoing costs for conducting a survey of the Avocet NR bridled nailtail wallaby population. This analysis assumes that both cage and hair trapping methods would be undertaken by paid Queensland Government employees. This comparison does not include costs associated with the initial planning phase of establishing either method at a new location as it is likely that they would incur a similar investment. The costs presented in this comparison take into account the cost of materials, fuel, staff wages ($45/hour), travel allowance ($100/night/staff member) and the cost of data analysis and genetic analysis (Appendix 6). Non-invasive Genetic Hair Sampling Page 14
3. Results Over the three months prior to cage trapping in June 2011, Avocet received a total of 304 mm of rainfall. No rainfall was recorded during the four nights of cage trapping. In the three months leading up to and including the first hair sampling period in October 2011, Avocet received a total of 16 mm of rainfall, predominantly in July. Before the end of the second hair sampling period, a hail storm went through the sampling area and 18 mm of rainfall was recorded. Prior to the third collection of samples, another 8 mm fell over Avocet. No rainfall was recorded between the third and final sampling period. 3.1 Cage Trapping During the cage trapping survey there was a total of 20 captures of bridled nailtail wallabies, comprising 10 individuals (Fig. 5). Of the 146 sites of which cage traps were set, 15 captured bridled nailtail wallabies in close proximity to Ron’s dam (Appendix 2). Four bridled nailtail wallabies were captured once while the other individuals were recaptured between 1-5 times over the eight trapping occasions (J. Lowry, pers. com. 2011). The number of new individuals declined over the survey period, while the total number of captures increased until the final trap night, when only two individuals were trapped (Fig. 5). 3.2 Hair trapping A total of 438 hair samples were visually identified as bridled nailtail wallaby from hair trap sites. Of these, 422 were genetically identified as bridled nailtail wallaby and were collected from 60 of the 180 hair trap sites (Appendix 4). All bridled nailtail wallaby hair samples were collected within Regional Ecosystems (REs) dominated by brigalow or poplar box at Avocet NR (Appendix 1). 8 New individuals 7 Total captures 6 No. BNTW 5 4 3 2 1 0 14-Jun 15-Jun 16-Jun 17-Jun Trap night Fig. 5. The number of bridled nailtail wallaby (BNTW) captures and new individuals in cage traps at Avocet NR in June 2011. Each trap night contains two trapping occasions. Non-invasive Genetic Hair Sampling Page 15
Using the corrected totals, runway traps were slightly more successful than triangle traps in sampling period A, while triangle traps were marginally more successful for sampling periods B, C and D (Fig. 6). An average of 45 traps (42-52) produced at least one bridled nailtail wallaby sample over the four sampling periods. There was no clear trend in the number of samples collected, with totals consistently high from the beginning to the end of the survey. For none of the four sampling occasions was there a statistical difference between the preference of the two trap types when collecting bridled nailtail wallaby samples (X2, 0.05) (Fig 6). Overall, triangle traps produced a significantly higher number of individual bridled nailtail wallaby hair samples per site than did runway traps (F1,58 = 11.46; P
8 Runway Triangle 7 Average no. samples per site 6 5 4 3 2 1 0 A B C D Sampling period Fig. 7. The average number of individual bridled nailtail wallaby hair samples collected from each runway and triangle trap for each of the four sampling periods. Bars show the minimum and maximum number of samples collected. 3.3 Genetic analysis Of the 438 hair samples identified using the hair ID key, 422 samples amplified on at least one locus, 223 amplified at the two loci B90 and YM148, and 71 amplified at four or more loci. Of the 223 samples that amplified at the B90 and YM148 loci, 93 were unique genotypes. Of the 71 that amplified on at least four loci, 66 were unique genotypes. Based on the samples that amplified at 4-5 loci (which produced the most reliable result), the minimum number of bridled nailtail wallabies known to be alive on Avocet NR is at least 66 individuals. None of the five loci were in Hardy-Weinberg Equilibrium (HWE). However, locus YM148 was close at 0.049. The FIS value for B90 was negative, indicating an excess of heterozygotes relative to HWE predictions. All other loci were positive, indicating an excess in homozygosity (Table 1). Micro-checker indicated homozygote excess at three of the five loci (B29, B87 and B151). This is most likely to result from null alleles which may be caused by allelic dropout or mutations at the primer site preventing the large allele from amplifying (Marucco et al. 2011). Given that locus B90 and YM148 showed only small deviations from expected frequencies, it is possible that additional sampling would reveal this population to be in HWE (Van Oosterhout et al. 2004). Non-invasive Genetic Hair Sampling Page 17
Table 1. The size range of the microsatellite fragment in base pairs (bps) (Sigg 2004), number of alleles at each locus (Na), significance value of a HWE test, inbreeding coefficient (FIS), error rate and total successful amplification at each locus for all samples. Locus Size range Na Significance FIS Error rate Successful (bp) estimate amplification B90 90-121 7 0.0126 -0.1020 0.088 0.782 YM148 96-106 9 0.0490 0.1857 0.083 0.889 B29 228-272 7 0.0001 0.2472 0* 0.154 B87 159-179 5 0.0000 0.4228 0.056 0.445 B151 206-246 12 0.0002 0.3222 0.079 0.232 * Insufficient data to complete error rate The total number of alleles detected across all five loci was 40 for the Avocet population. The error rate for each locus was calculated based on replicated PCR runs and is given in Table 1. The successful amplification of locus B90 and YM148 was much higher than that of the other three loci whose amplification was less than 45% (Table 1). Due to the low rate of amplification for B29, repeated PCRs failed to produce enough data to enable an error rate to be determined. The size range of each locus appears to be correlated with the error rate. Loci with shorter fragment lengths had higher amplification success (Table 1). The unique genotypes from the 66 individuals that amplified at four or more loci have a combined exclusion probability (identity) of 0.999 and a combined exclusion probability (sibling identity) of 0.984 (Kalinowski et al. 2007). It is, therefore, highly unlikely that any two individuals will share the same genotype by chance when using the five highly polymorphic microsatellite loci used in this study. The program COLONY detected a large proportion of half siblings suggesting that individuals in the population are closely related (Fig. 8). COLONY calculated the effective population size (Ne) of the bridled nailtail wallaby population on Avocet NR to be 29 with 95% creditable limits of 18 to 50 (Jones & Wang 2009). The effective population size was also calculated using the program ONeSAMP, which estimated Ne at 35, 95% creditable limits of 24 to 61 (Tallmon et al. 2008). Non-invasive Genetic Hair Sampling Page 18
Fig. 8. Graph produced by COLONY showing the relationship between individuals based on the 66 genotypes that amplified at four to five loci. 3.4 Cost comparison The cost of initially setting up and conducting a CMR survey on bridled nailtail wallabies on Avocet NR was $70,676 for cage trapping and $46,613 for a hair trapping census. The cost of conducting additional surveys would be $28,546 for cage trapping and $40,857 for hair trapping (Appendix 6). Cage trapping is therefore, approximately 34% more expensive to set up than hair trapping, while subsequent surveys for hair sampling cost 30% more due to additional costs associated with genetic analysis (Fig. 9) (Appendix 6). Non-invasive Genetic Hair Sampling Page 19
$45,000 Initial Costs Cage Trap $45,000 Ongoing Costs Cage Trap $40,000 Hair Trap $40,000 Hair Trap $35,000 $35,000 $30,000 $30,000 $25,000 $25,000 $20,000 $20,000 $15,000 $15,000 $10,000 $10,000 $5,000 $5,000 $0 $0 Materials Staff Vehicles Analysis Materials Staff Vehicles Analysis Fig. 9. The initial set up cost to undertake a cage or hair trapping survey assuming all new equipment was purchased. Ongoing costs for additional hair and cage trap surveys including maintenance and small replacement costs. Cost itemised by materials, staff wages, vehicle costs, and analysis which includes modelling and contracted genetic processing. The cost of the materials (primarily cage traps) for setting up a cage trapping survey, at over $43,492, is significantly higher than for hair trapping at $3,600. The ongoing analysis costs of conducting a hair survey are significantly higher due to the $12,000 required to genotype 400 hair samples (Appendix 6). The costs associated with staff time and vehicle costs are relatively similar for both survey methods (Fig. 9). Non-invasive Genetic Hair Sampling Page 20
4. Discussion 4.1 General observations In June 2012, Avocet NR had received a total 304 mm of rainfall in the three months leading up to cage trapping which resulted in much of the area maintaining good grass and herb coverage prior to pre-feeding. With a high proportion of feed available bridled nailtail wallabies may have been less likely to enter the cage traps (Barnett & Dutton 1995). For the hair trap survey in October 2011, Avocet NR was extremely dry with little obvious fresh pasture growth during pre- feeding and over the first sampling period. Rainfall in the first and second week saw the survey area get progressively greener, with fresh grass and forb growth across previously bare areas. The amount of fresh pasture growth appeared to have no influence on the number of successful hair sites during the 18 night hair trapping survey which remained consistently between 42-52 sites (Fig. 6). The number of individual bridled nailtail wallabies visiting each site in response to pasture growth remains untested due the small number of individuals that were able to be identified genetically. Bridled nailtail wallaby samples were not obtained from any of the cage or hair traps around the release site near the northern dam (Alice’s dam) which is surrounded by mature brigalow with little understorey (Appendices 1 & 2). Fisher (2000) found that bridled nailtail wallabies from the remnant population at Taunton NP did not inhabit open grassy woodlands, such as silver-leafed ironbark communities, preferring a vegetation structure that consists of more woody stems and a dense canopy cover at 0.5 m high (predominantly brigalow regrowth). The single outlier hair trap that collected numerous bridled nailtail wallaby hair samples was approximately 1.5 km away from all the other successful hair traps. This hair trap was situated within a poplar box, brigalow mix community (RE 11.10.12) with some low scrubs on the edge of a 20-30 m wide cleared fire break that adjoins a cleared paddock (Appendix 1). All the other hair traps around this block (eastern and northern side) were placed under small sparse regrowth within the cleared firebreak. 4.2 Trap method The cage trapping event on Avocet NR conducted in June 2011, appears to have provided an inaccurate count of the number of bridled nailtail wallabies, trapping a total of 10 individuals. This is significantly lower than the 66 individuals identified through genotyping of hair trapped samples in October 2011. The low cage capture rate may have been influenced by factors such as trap preparation. As a result of the long wet season and the use of the cage traps from other sites, pre-feeding was only conducted for one week prior to trapping instead of the minimum two weeks normally undertaken (J. Lowry, pers. com. 2011). Unlike the other two bridled nailtail Non-invasive Genetic Hair Sampling Page 21
wallaby populations in Qld, the Avocet NR population is infrequently surveyed and animals may have become unfamiliar with cage traps and attractants, possibly resulting in the low trapping rates (Barnett & Dutton 1995). If this were the case, then a longer pre-feeding time would need to be employed to ensure the likelihood that trapping individuals is consistent with that of other bridled nailtail wallaby populations. With the likelihood of a small population size of bridled nailtail wallabies on Avocet NR, rectifying these inconsistencies may improve the capture rate. However, cage trapping is still unlikely to provide an accurate estimate of population size owing to the small sample size obtained from live trapping (Henry & Russello 2011). For small population sizes (0.30 are required for accurate population estimates (Marucco et al. 2011), with recaptures occurring over multiple occasions for all individuals (Minta et al. 1989). Non-invasive genetic sampling has the advantage of obtaining recaptures from a single sampling session (Petit & Valiere 2006). Capture heterogeneity can be a concern when using trapping data in models such as CMR (Minta et al. 1989; Dreher et al. 2009; Ebert et al. 2010). Overall, there appears to be no substantial evidence to suggest a difference in sex ratio when cage trapping bridled nailtail wallabies (Kingsley et al. 2012), although this should be tested for hair sampling by sexing DNA samples from hair trapped individuals. The use of triangle and runway hair traps across Avocet NR proved very effective. This technique collected 422 bridled nailtail wallaby hair samples on the four sampling occasions (18 nights), compared to cage trapping returning only 20 captures over eight occasions (four nights). Triangles produced more hair samples per trap than did runway traps. It is unclear whether these additional samples provided additional information on the number of individuals or population size of bridled nailtail wallabies due to the low percentage of identifiable individuals. To minimise stress and to reduce any major impact live trapping may have on the population, cage trapping is restricted to four consecutive nights (J. Lowry, pers. com. 2011). The ability to leave hair traps out for an extended period of time significantly increases the sample size for use in a CMR population size estimate (Banks et al. 2003; Depue & Ben-David 2007; Henry & Russello 2011). Hair trapping also enables a much larger area to be surveyed given that traps do not have to be checked at a particular time of day or on a daily basis (Mowat & Strobeck 2000; Stenglein et al. 2011). When cage trapping bridled nailtail wallabies, the number of captures generally declines over the trapping period with many individuals only trapped on one occasion. Cage trapping the remnant Taunton NP population in 2012 saw 68% of bridled nailtail wallabies trapped only once, over 11 occasions (A. Dinwoodie, pers. com. 2012). Where species show an avoidance to traps, Non-invasive Genetic Hair Sampling Page 22
You can also read