Free Standard AU & NZ Shipping For All Book Orders Over $80!
Register      Login
Wildlife Research Wildlife Research Society
Ecology, management and conservation in natural and modified habitats
RESEARCH ARTICLE (Open Access)

Evaluation of genetic markers for the metabarcoding of Australian marsupials from predator scats

Catriona D. Campbell https://orcid.org/0000-0001-7983-1832 A , Anna J. MacDonald https://orcid.org/0000-0003-2972-200X A B and Stephen D. Sarre https://orcid.org/0000-0002-7158-2517 A *
+ Author Affiliations
- Author Affiliations

A Centre for Conservation Ecology and Genomics, Institute for Applied Ecology, University of Canberra, Bruce, ACT 2617, Australia.

B Current Address: Australian Antarctic Division, Kingston, Tas. 7050, Australia.

* Correspondence to: stephen.sarre@canberra.edu.au

Handling Editor: Adam Stow

Wildlife Research 51, WR23134 https://doi.org/10.1071/WR23134
Submitted: 12 October 2023  Accepted: 23 May 2024  Published: 2 July 2024

© 2024 The Author(s) (or their employer(s)). Published by CSIRO Publishing. This is an open access article distributed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License (CC BY-NC-ND)

Abstract

Context

DNA recovered from predator faeces (scats) can be used to determine the presence of fauna and shed light on their life histories and inter-species interactions. DNA metabarcoding, which involves concurrent amplification and sequencing of DNA from multiple taxa, represents an important advance by enabling the simultaneous detection of multiple species from such samples. Although an attractive proposition, metabarcoding requires ‘universally’ applicable genetic markers that can discriminate among a broad range of taxa, while also targeting sequences that are sufficiently short to be amplified from degraded DNA.

Aims

To identify, evaluate, and test metabarcoding DNA markers suitable for the detection of marsupials and other Australian fauna from terrestrial predator scats found in nature.

Methods

We apply a bioinformatic approach using publicly available DNA databases and a locally derived and marker-specific reference-DNA database to evaluate the diagnostic ability and likelihood of amplification of candidate metabarcoding markers for marsupials and other taxa that may be consumed by predators.

Key results

We identify two markers (12SV5 and 16SMam) that are suitable for use and successfully identify marsupial sequences at a high level of resolution. These markers work best in combination because they bring complementary levels of primer specificity and diagnostic ability in detecting multiple prey species as well as the predator. We also show that these samples work well in predator scats sampled from the wild in Tasmania.

Conclusions

These markers provide a useful tool for surveying mammalian predators and their prey and could also be applied to eDNA analyses from other sample types. Improvements to the reference database and further development of markers targeting different taxonomic groups will improve the resolution and usefulness of this approach.

Implications

Metabarcoding of predator scats provides a potent approach to non-invasive wildlife survey that offers the opportunity for the detection of multiple species across all vertebrates.

Keywords: applied ecology, conservation biology, endangered species, faeces, invasive species, molecular ecology, predator–prey interactions, threatened species, vertebrates.

References

Andersen GE, Johnson CN, Barmuta LA, Jones ME (2017) Dietary partitioning of Australia’s two marsupial hypercarnivores, the Tasmanian devil and the spotted-tailed quoll, across their shared distributional range. PLoS ONE 12, e0188529.
| Crossref | Google Scholar |

Applied Biosystems (2009) ‘DNA sequencing by capillary electrophoresis – chemistry guide.’ 2nd edn. (Springer) doi:10.1002/elps.200900218.DNA

Archer M, Bates H, Hand SJ, Evans T, Broome L, McAllan B, Geiser F, Jackson S, Myers T, Gillespie A, Palmer C, Hawke T, Horn AM (2019) The Burramys Project: a conservationist’s reach should exceed history’s grasp, or what is the fossil record for? Philosophical Transactions of the Royal Society B: Biological Sciences 374, 20190221.
| Crossref | Google Scholar | PubMed |

Ashman KR, Watchorn DJ, Whisson DA (2019) Prioritising research efforts for effective species conservation: a review of 145 years of koala research. Mammal Review 49, 189-200.
| Crossref | Google Scholar |

Avise JC (1986) Mitochondrial DNA and the evolutionary genetics of higher animals. Philosophical Transactions of the Royal Society B: Biological Sciences 312, 325-342.
| Crossref | Google Scholar |

Benson DA, Cavanaugh M, Clark K, Karsch-Mizrachi I, Lipman DJ, Ostell J, Sayers EW (2013) GenBank. Nucleic Acids Research 41, D36-D42.
| Crossref | Google Scholar | PubMed |

Berry O, Sarre SD (2007) Gel-free species identification using melt-curve analysis. Molecular Ecology Notes 7, 1-4.
| Crossref | Google Scholar |

Berry TE, Osterrieder SK, Murray DC, Coghlan ML, Richardson AJ, Grealy AK, Stat M, Bejder L, Bunce M (2017) DNA metabarcoding for diet analysis and biodiversity: a case study using the endangered Australian sea lion (Neophoca cinerea). Ecology and Evolution 7, 5435-5453.
| Crossref | Google Scholar |

Bolger AM, Lohse M, Usadel B (2014) Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114-2120.
| Crossref | Google Scholar | PubMed |

Boyer FF, Mercier CC, Bonin AA, Le Bras Y, Taberlet P, Coissac E (2016) OBITOOLS: a UNIX-inspired software package for DNA metabarcoding. Molecular Ecology Resources 16, 176-182.
| Crossref | Google Scholar | PubMed |

Brown SDJ, Collins RA, Boyer S, Lefort M-C, Malumbres-Olarte J, Vink CJ, Cruickshank RH (2012) Spider: an R package for the analysis of species identity and evolution, with particular reference to DNA barcoding. Molecular Ecology Resources 12, 562-565.
| Crossref | Google Scholar | PubMed |

Burk A, Springer MS (2000) Intergeneric relationships among macropodoidea (Metatheria:Diprotodontia) and the chronicle of kangaroo evolution. Journal of Mammalian Evolution 7, 213-237.
| Crossref | Google Scholar |

Cannon MV, Hester J, Shalkhauser A, Chan ER, Logue K, Small ST, Serre D (2016) In silico assessment of primers for eDNA studies using PrimerTree and application to characterize the biodiversity surrounding the Cuyahoga River. Scientific Reports 6, 22908.
| Crossref | Google Scholar | PubMed |

Coissac E, Riaz T, Puillandre N (2012) Bioinformatic challenges for DNA metabarcoding of plants and animals. Molecular Ecology 21, 1834-1847.
| Crossref | Google Scholar | PubMed |

Collins RA, Cruickshank RH (2013) The seven deadly sins of DNA barcoding. Molecular Ecology Resources 13, 969-975.
| Crossref | Google Scholar | PubMed |

Compson ZG, McClenaghan B, Singer GAC, Fahner NA, Hajibabaei M (2020) Metabarcoding from microbes to mammals: comprehensive bioassessment on a global scale. Frontiers in Ecology and Evolution 8, 581835.
| Crossref | Google Scholar |

Curry CJ, Gibson JF, Shokralla S, Hajibabaei M, Baird DJ (2018) Identifying North American freshwater invertebrates using DNA barcodes: are existing COI sequence libraries fit for purpose? Freshwater Science 37, 178-189.
| Crossref | Google Scholar |

Dasmahapatra KK, Mallet J (2006) Taxonomy: DNA barcodes: recent successes and future prospects. Heredity 97, 254-255.
| Crossref | Google Scholar | PubMed |

Deagle BE, Tollit DJ, Jarman SN, Hindell MA, Trites AW, Gales NJ (2005) Molecular scatology as a tool to study diet: analysis of prey DNA in scats from captive Steller sea lions. Molecular Ecology 14, 1831-1842.
| Crossref | Google Scholar | PubMed |

Deagle BE, Jarman SN, Coissac E, Pompanon F, Taberlet P (2014) DNA metabarcoding and the cytochrome c oxidase subunit I marker: not a perfect match. Biology Letters 10, 20140562.
| Crossref | Google Scholar | PubMed |

Edgar RC (2004) MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Research 32, 1792-1797.
| Crossref | Google Scholar | PubMed |

Ekrem T, Willassen E, Stur E (2007) A comprehensive DNA sequence library is essential for identification with DNA barcodes. Molecular Phylogenetics and Evolution 43, 530-542.
| Crossref | Google Scholar | PubMed |

Elbrecht V, Leese F (2016) PrimerMiner: an R package for development and in silico validation of DNA metabarcoding primers. Methods in Ecology and Evolution 8, 622-626.
| Crossref | Google Scholar |

Fancourt BA, Hawkins CE, Nicol SC (2018) Mechanisms of climate-change-induced species decline: spatial, temporal and long-term variation in the diet of an endangered marsupial carnivore, the eastern quoll. Wildlife Research 45, 737-750.
| Crossref | Google Scholar |

Ficetola GF, Taberlet P (2023) Towards exhaustive community ecology via DNA metabarcoding. Molecular Ecology 32, 6320-6329.
| Crossref | Google Scholar |

Frankham GJ, Handasyde KA, Eldridge MDB (2012) Novel insights into the phylogenetic relationships of the endangered marsupial genus Potorous. Molecular Phylogenetics and Evolution 64, 592-602.
| Crossref | Google Scholar | PubMed |

Goldberg CS, Strickler KM, Fremier AK (2018) Degradation and dispersion limit environmental DNA detection of rare amphibians in wetlands: increasing efficacy of sampling designs. Science of The Total Environment 633, 695-703.
| Crossref | Google Scholar | PubMed |

Hajibabaei M, Singer GAC, Hebert PDN, Hickey DA (2007) DNA barcoding: how it complements taxonomy, molecular phylogenetics and population genetics. Trends in Genetics 23, 167-172.
| Crossref | Google Scholar | PubMed |

Hebert PDN, Ratnasingham S, de Waard JR (2003a) Barcoding animal life: cytochrome c oxidase subunit 1 divergences among closely related species. Proceedings of the Royal Society of London. Series B: Biological Sciences 270, S96-S99.
| Crossref | Google Scholar |

Hebert PDN, Cywinska A, Ball SL, deWaard JR (2003b) Biological identifications through DNA barcodes. Proceedings of the Royal Society of London. Series B: Biological Sciences 270, 313-321.
| Crossref | Google Scholar |

Hebert PDN, Penton EH, Burns JM, Janzen DH, Hallwachs W (2004) Ten species in one: DNA barcoding reveals cryptic species in the neotropical skipper butterfly Astraptes fulgerator. Proceedings of the National Academy of Sciences of the United States of America 101, 14812-14817.
| Crossref | Google Scholar |

Ji Y, Ashton L, Pedley SM, Edwards DP, Tang Y, Nakamura A, Kitching R, Dolman PM, Woodcock P, Edwards FA, Larsen TH, Hsu WW, Benedick S, Hamer KC, Wilcove DS, Bruce C, Wang X, Levi T, Lott M, Emerson BC, Yu DW (2013) Reliable, verifiable and efficient monitoring of biodiversity via metabarcoding. Ecology Letters 16, 1245-1257.
| Crossref | Google Scholar | PubMed |

Jones ME, Barmuta LA (1998) Diet overlap and relative abundance of sympatric dasyurid carnivores: a hypothesis of competition. Journal of Animal Ecology 67, 410-421.
| Crossref | Google Scholar |

Katoh K, Misawa K, Kuma K, Miyata T (2002) MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Research 30, 3059-3066.
| Crossref | Google Scholar | PubMed |

Keck F, Blackman RC, Bossart R, Brantschen J, Couton M, Hurlemann S, Kirschner D, Locher N, Zhang H, Altermatt F (2022) Meta-analysis shows both congruence and complementarity of DNA and eDNA metabarcoding to traditional methods for biological community assessment. Molecular Ecology 31, 1820-1835.
| Crossref | Google Scholar | PubMed |

Keck F, Couton M, Altermatt F (2023) Navigating the seven challenges of taxonomic reference databases in metabarcoding analyses. Molecular Ecology Resources 23, 742-755.
| Crossref | Google Scholar | PubMed |

Kimura M (1980) A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. Journal of Molecular Evolution 16, 111-120.
| Crossref | Google Scholar | PubMed |

Kitano T, Umetsu K, Tian W, Osawa M (2007) Two universal primer sets for species identification among vertebrates. International Journal of Legal Medicine 121, 423-427.
| Crossref | Google Scholar | PubMed |

Kumar S, Stecher G, Tamura K (2016) MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Molecular Biology and Evolution 33, 1870-1874.
| Crossref | Google Scholar |

Lambert DM, Baker A, Huynen L, Haddrath O, Hebert PDN, Millar CD (2005) Is a large-scale DNA-based inventory of ancient life possible? Journal of Heredity 96, 279-284.
| Crossref | Google Scholar | PubMed |

Lindner MF, Gonçalves LT, Bianchi FM, Ferrari A, Cavalleri A (2023) Tiny insects, big troubles: a review of BOLD’s COI database for Thysanoptera (Insecta). Bulletin of Entomological Research 113, 1-13.
| Crossref | Google Scholar |

MacDonald AJ, Sarre SD (2015) Species assignment from trace DNA sequences: an in silico assessment of the test used to survey for foxes in Tasmania. Journal of Applied Ecology 52, 1649-1655.
| Crossref | Google Scholar |

MacDonald AJ, Sarre SD (2017) A framework for developing and validating taxon-specific primers for specimen identification from environmental DNA. Molecular Ecology Resources 17, 708-720.
| Crossref | Google Scholar | PubMed |

MacDonald AJ, Sarre SD, FitzSimmons NN, Aitken N (2011) Determining microsatellite genotyping reliability and mutation detection ability: an approach using small-pool PCR from sperm DNA. Molecular Genetics and Genomics 285, 1-18.
| Crossref | Google Scholar | PubMed |

Macqueen P, Seddon JM, Austin JJ, Hamilton S, Goldizen AW (2010) Phylogenetics of the pademelons (Macropodidae: Thylogale) and historical biogeography of the Australo-Papuan region. Molecular Phylogenetics and Evolution 57, 1134-1148.
| Crossref | Google Scholar | PubMed |

May-Collado LJ, Kilpatrick CW, Agnarsson I (2015) Mammals from ‘down under’: a multi-gene species-level phylogeny of marsupial mammals (Mammalia, Metatheria). PeerJ 3, e805.
| Crossref | Google Scholar | PubMed |

McInnes JC, Alderman R, Deagle BE, Lea M-A, Raymond B, Jarman SN (2017) Optimised scat collection protocols for dietary DNA metabarcoding in vertebrates. Methods in Ecology and Evolution 8, 192-202.
| Crossref | Google Scholar |

McLennan EA, Wise P, Lee AV, Grueber CE, Belov K, Hogg CJ (2022) DNA metabarcoding reveals a broad dietary range for Tasmanian devils introduced to a naive ecosystem. Ecology and Evolution 12, e8936.
| Crossref | Google Scholar | PubMed |

Mena JL, Yagui H, Tejeda V, Bonifaz E, Bellemain E, Valentini A, Tobler MW, Sánchez-Vendizú P, Lyet A (2021) Environmental DNA metabarcoding as a useful tool for evaluating terrestrial mammal diversity in tropical forests. Ecological Applications 31, e02335.
| Crossref | Google Scholar | PubMed |

Meredith RW, Westerman M, Springer MS (2009) A phylogeny of Diprotodontia (Marsupialia) based on sequences for five nuclear genes. Molecular Phylogenetics and Evolution 51, 554-571.
| Crossref | Google Scholar | PubMed |

Meusnier I, Singer GAC, Landry J-F, Hickey DA, Hebert PDN, Hajibabaei M (2008) A universal DNA mini-barcode for biodiversity analysis. BMC Genomics 9, 214.
| Crossref | Google Scholar | PubMed |

Mitchell KJ, Pratt RC, Watson LN, Gibb GC, Llamas B, Kasper M, Edson J, Hopwood B, Male D, Armstrong KN, Meyer M, Hofreiter M, Austin J, Donnellan SC, Lee MSY, Phillips MJ, Cooper A (2014) Molecular phylogeny, biogeography, and habitat preference evolution of marsupials. Molecular Biology and Evolution 31, 2322-2330.
| Crossref | Google Scholar | PubMed |

Modave E (2017) Identification, distribution and diet of Tasmanian predators inferred by scat DNA. Doctoral thesis, University of Canberrra, ACT, Australia.

Modave E, MacDonald AJ, Sarre SD (2017) A single mini-barcode test to screen for Australian mammalian predators from environmental samples. GigaScience 6, 1-13.
| Crossref | Google Scholar | PubMed |

Mugnai F, Costantini F, Chenuil A, Leduc M, Gutiérrez Ortega JM, Meglécz E (2023) Be positive: customized reference databases and new, local barcodes balance false taxonomic assignments in metabarcoding studies. PeerJ 7, e14616.
| Crossref | Google Scholar | PubMed |

Munro NT, McIntyre S, Macdonald B, Cunningham SA, Gordon IJ, Cunningham RB, Manning AD (2019) Returning a lost process by reintroducing a locally extinct digging marsupial. PeerJ 2019, 1-24.
| Crossref | Google Scholar |

Murray DC, Haile J, Dortch J, White NE, Haouchar D, Bellgard MI, Allcock RJ, Prideaux GJ, Bunce M (2013) Scrapheap challenge: a novel bulk-bone metabarcoding method to investigate ancient DNA in faunal assemblages. Scientific Reports 3, 3371.
| Crossref | Google Scholar |

Osborne MJ, Christidis L, Norman JA (2002) Molecular phylogenetics of the Diprotodontia (kangaroos, wombats, koala, possums, and allies). Molecular Phylogenetics and Evolution 25, 219-228.
| Crossref | Google Scholar |

Pääbo S, Poinar H, Serre D, Jaenicke-Després V, Hebler J, Rohland N, Kuch M, Krause J, Vigilant L, Hofreiter M (2004) Genetic analyses from ancient DNA. Annual Review of Genetics 38, 645-679.
| Crossref | Google Scholar | PubMed |

Pompanon F, Deagle BE, Symondson WOC, Brown DS, Jarman SN, Taberlet P (2012) Who is eating what: diet assessment using next generation sequencing. Molecular Ecology 21, 1931-1950.
| Crossref | Google Scholar | PubMed |

Porter TM, Hajibabaei M (2018) Scaling up: a guide to high-throughput genomic approaches for biodiversity analysis. Molecular Ecology 27 313-338.
| Crossref | Google Scholar |

Ramsey DSL, MacDonald AJ, Quasim S, Barclay C, Sarre SD (2015) An examination of the accuracy of a sequential PCR and sequencing test used to detect the incursion of an invasive species: the case of the red fox in Tasmania. Journal of Applied Ecology 52, 562-570.
| Crossref | Google Scholar |

Ramsey DSL, Barclay C, Campbell CD, Dewar E, MacDonald AJ, Modave E, Quasim S, Sarre SD (2018) Detecting rare carnivores using scats: Implications for monitoring a fox incursion into Tasmania. Ecology and Evolution 8, 732-743.
| Crossref | Google Scholar | PubMed |

Ratnasingham S, Hebert PDN (2007) BOLD: the barcode of life data system (http://www.barcodinglife.org). Molecular Ecology Notes 7, 355-364.
| Crossref | Google Scholar | PubMed |

Reside AE, Briscoe NJ, Dickman CR, Greenville AC, Hradsky BA, Kark S, Kearney MR, Kutt AS, Nimmo DG, Pavey CR, Read JL, Ritchie EG, Roshier D, Skroblin A, Stone Z, West M, Fisher DO (2019) Persistence through tough times: fixed and shifting refuges in threatened species conservation. Biodiversity and Conservation 28, 1303-1330.
| Crossref | Google Scholar |

Riaz T, Shehzad W, Viari A, Pompanon F, Taberlet P, Coissac E (2011) EcoPrimers: inference of new DNA barcode markers from whole genome sequence analysis. Nucleic Acids Research 39, e145.
| Crossref | Google Scholar | PubMed |

Sarre SD, MacDonald AJ, Barclay C, Saunders GR, Ramsey DSL (2013) Foxes are now widespread in Tasmania: DNA detection defines the distribution of this rare but invasive carnivore. Journal of Applied Ecology 50, 459-468.
| Crossref | Google Scholar |

Schnell IB, Bohmann K, Gilbert MTP (2015) Tag jumps illuminated – reducing sequence-to-sample misidentifications in metabarcoding studies. Molecular Ecology Resources 15, 1289-1303.
| Crossref | Google Scholar |

Seymour M, Durance I, Cosby BJ, Ransom-Jones E, Deiner K, Ormerod SJ, Colbourne JK, Wilgar G, Carvalho GR, de Bruyn M, Edwards F, Emmett BA, Bik HM, Creer S (2018) Acidity promotes degradation of multi-species environmental DNA in lotic mesocosms. Communications Biology 1, 4.
| Crossref | Google Scholar |

Springer MS, Westerman M, Kirsch JAW (1994) Relationships among orders and families of marsupials based on 12S ribosomal DNA sequences and the timing of the marsupial radiation. Journal of Mammalian Evolution 2, 85-115.
| Crossref | Google Scholar |

Springer MS, Hollar LJ, Burk A (1995) Compensatory substitutions and the evolution of the mitochondrial 12S rRNA gene in mammals. Molecular Biology and Evolution 12, 1138-1150.
| Google Scholar |

Staats M, Arulandhu AJ, Gravendeel B, Holst-Jensen A, Scholtens I, Peelen T, Prins TW, Kok E (2016) Advances in DNA metabarcoding for food and wildlife forensic species identification. Analytical and Bioanalytical Chemistry 408, 4615-4630.
| Crossref | Google Scholar | PubMed |

Taberlet PP, Coissac E, Pompanon F, Brochmann C, Willerslev E (2012) Towards next-generation biodiversity assessment using DNA metabarcoding. Molecular Ecology 21, 2045-2050.
| Crossref | Google Scholar | PubMed |

Tamura K, Nei M (1993) Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Molecular Biology and Evolution 10, 512-526.
| Crossref | Google Scholar | PubMed |

Taylor PG (1996) Reproducibility of ancient DNA sequences from extinct Pleistocene fauna. Molecular Biology and Evolution 13, 283-285.
| Crossref | Google Scholar |

Westerman M, Burk A, Amrine-Madsen HM, Prideaux GJ, Case JA, Springer MS (2002) Molecular evidence for the last survivor of an ancient kangaroo lineage. Journal of Mammalian Evolution 9, 209-223.
| Google Scholar |

Westerman M, Loke S, Springer MS (2004) Molecular phylogenetic relationships of two extinct potoroid marsupials, Potorous platyops and Caloprymnus campestris (Potoroinae: Marsupialia). Molecular Phylogenetics and Evolution 31, 476-485.
| Crossref | Google Scholar | PubMed |

Westerman M, Krajewski C, Kear BP, Meehan L, Meredith RW, Emerling CA, Springer MS (2016) Phylogenetic relationships of dasyuromorphian marsupials revisited. Zoological Journal of the Linnean Society 176, 686-701.
| Crossref | Google Scholar |

Wilkes Walburn J, Rourke ML, Furlan E, DiBattista JD, Broadhurst MK, Fowler AM, Hughes JM, Fielder S (2022) Robust environmental DNA assay development and validation: a case study with two vulnerable Australian fish. Aquatic Conservation: Marine and Freshwater Ecosystems 32, 1225-1231.
| Crossref | Google Scholar |

Yamamoto S, Masuda R, Sato Y, Sado T, Araki H, Kondoh M, Minamoto T, Miya M (2017) Environmental DNA metabarcoding reveals local fish communities in a species-rich coastal sea. Scientific Reports 7, 40368.
| Crossref | Google Scholar | PubMed |

Yang L, Tan Z, Wang D, Xue L, Guan M-X, Huang T, Li R (2014) Species identification through mitochondrial rRNA genetic analysis. Scientific Reports 4, 4089.
| Crossref | Google Scholar | PubMed |

Zeale MRK, Butlin RK, Barker GLA, Lees DC, Jones G (2011) Taxon-specific PCR for DNA barcoding arthropod prey in bat faeces. Molecular Ecology Resources 11, 236-244.
| Crossref | Google Scholar | PubMed |

Zwickl DJ (2006) Genetic algorithm approaches for the phylogenetic analysis of large biological sequence datasets under the maximum likelihood criterion. PhD dissertation, University of Texas, USA. Available at http://www.zo.utexas.edu/faculty/antisense/garli/Garli.html