- Short Report
- Open Access
Development of microsatellite markers for Hyacinth macaw (Anodorhynchus hyacinthinus) and their cross-amplification in other parrot species
BMC Research Notesvolume 8, Article number: 736 (2015)
Hyacinth macaw Anodorhynchus hyacinthinus is the largest parrot of the world and is considered vulnerable to extinction due to its habitat loss and illegal trade associated to the international pet market demand. Genetic studies on this species are still incipient to generate a consistent characterization of the population dynamics and to develop appropriate conservation strategies. In this sense, microsatellite markers may support the detection of a population genetic structure for this bird species. However, at this time, none Hyacinth macaw species-specific primers for microsatellite loci have been so far established. This study aimed to develop and characterize polymorphic microsatellite markers for A. hyacinthinus and to check for their cross-amplification in other parrot species.
Sequences containing repeated dinucleotide motifs were prospected and optimized from a genomic library that was enriched for microsatellites using magnetic beads. The analyses of 43–57 samples from wild individuals of three distinct Brazilian subpopulations led to the characterization of five polymorphic microsatellite loci. Allele richness per locus ranged from two to 12. Three loci exhibited observed heterozygosity values higher than 50 %, but the overall average value among all loci was close to 45 %. In addition, successful primer cross-amplification was verified in seven other investigated species of Neotropical parrots.
The newly developed markers have shown to be potentially useful for in situ and ex situ population studies to support future conservation actions of Hyacinth macaw and other parrots.
Hyacinth macaw Anodorhynchus hyacinthinus (Latham, 1790; Psittaciformes, Psittacidae) is the largest parrot of the world and has great importance in conservation biology due to its flagship species status [1, 2]. In Brazil, the species is mainly distributed in three non-continuous regions that correspond to different biomes, exhibiting dissimilar faunal and floral compositions. In these areas, Hyacinth macaw has different preference for food items and nesting sites [3, 4]. Moreover, their populations are sensitive to the continuous habitat loss and have also been directly affected by the illegal trade, which led this species to be considered as vulnerable to extinction . Genetic studies have gained prominence to aid in the conservation of a wide-range endangered vertebrates [6–9], including bird species . Among the genetic markers, microsatellite loci have shown to be useful in variable contexts, notably to estimate the influence of genetic components in population studies , to trace the geographic origin of unknown individuals  or to determine the kinship of siblings , among other applications. Currently, there is a very limited source of microsatellite loci characterized in Neotropical parrots, despite of the species richness of this group . Available studies on Hyacinth macaw have used nonspecific microsatellite markers to assess the genetic variability of this species , their population genetic structure [15, 16] and to identify the probable geographic origin of a rescued individual from captivity . Despite one pair of primers that was developed for A. hyacinthinus (Davis S, unpublished data), all these studies used heterologous primers characterized for other parrot species such as Ara ararauna [17, 18], Amazona guildingii [19, 20], and Psittacus erithacus . However, these other birds do not represent sister taxa of the genus Anodorhynchus and are not closely interrelated [22, 23]. In this sense, the use of heterologous microsatellite primers generally leads to a decrease in the polymorphism level as the phylogenetic distance among species increases [24–26]. In fact, for Hyacinth macaw, the number of alleles and observed heterozygosity for all analysed heterologous loci were considerably lower when compared to those species from which the primers were originally described . In addition, a lower genetic variability has been observed in threatened parrots when compared to other non-threatened species . So far, there are no species-specific primers for microsatellite loci available for Hyacinth macaw described in literature. Thus, considering the vulnerable status of this species, a proper assessment of the genetic background of extant populations requires the analysis of a considerable number of loci. Consequently, the development of microsatellite specific primers for A. hyacinthinus could enlarge the rol of available molecular tools for the accurate assessment of its genetic diversity and population structure. Thereby, in this study we report the isolation and characterization of microsatellite loci for Hyacinth macaw and their cross-amplification in other species of Neotropical parrots.
Total genomic DNA from 57 individuals of A. hyacinthinus from distinct Brazilian regions (Table 1) was isolated following the extraction protocol using proteinase K and phenol–chloroform . A DNA sample of one animal was used on a digestion assay and fragments between 700 and 1200 bp were selected by digestion with the restriction enzyme AfaI (Invitrogen) and further ligated to adapters with known sequences (Rsa21 and Rsa25). Then these fragments were enriched throughout hybridization to probes containing dinucleotide repeats (GT and CT) chemically linked to biotin and selected using magnetic beads covered with streptavidin. Subsequently, the recovered fragments were inserted into a cloning vector pGEM-T Easy (Promega), and electroporation was used to transform competent cells of Escherichia coli (XL1-Blue). Then, plasmid extraction was performed using the alkaline lysis method , and T7 and SP6 primers were used for sequencing the library inserts with an ABIPrism 3500xL Genetic Analyzer (Applied Biosystems). The protocol used for microsatellite enrichment technique was described in Billotte et al.  with some modifications. Generated nucleotide sequences were screened for microsatellites regions in the program SSRIT . Primer3 plus  and OligoAnalyzer 3.1 (Integrated DNA Technologies) softwares were used for the design of the species-specific primers and to verify the occurrence of hairpins, self-dimers, and heterodimers. PCR conditions were as follows: initial denaturation at 95 °C for 10 min, followed by 35 cycles at 95 °C for 1 min, 45–65 °C for 40 s (nine annealing temperatures were applied: 45, 53, 55, 57, 59, 60, 61, 63, and 65 °C), 72 °C for 40 s, and a final extension at 72 °C for 10 min. Each PCR contained 20–50 ng of DNA, 0.1 μL of TaqDNA polymerase (5 U/μL, Invitrogen), 1 μL of each primer (10 mM), 0.4 μL of MgCl2 (25 mM), 1.2 μL of PCR buffer (10×), 1 μL of dNTPs (2.5 mM), and autoclaved Milli-Q water to complete 12 μL. Approximately 3 μL of the products were electrophoresed through 1.5 % agarose gels to verify amplification success.
For the genetic analyses, we collected approximately 0.1 mL of peripheral blood from the brachial veins of the lower wing of chicks found in natural nests, using disposable syringes and needles. The nests were accessed by an alpinism technique adapted to trees. Subsequently, the blood samples were transferred into microtubes containing approximately 0.5 mL of absolute ethanol and kept at air temperature. In addition, several research partners provided samples that were collected in previous field campaigns (details in Table 1). This work was approved by the federal government authoritative (IBAMA; permission Number: 36590-1). Microsatellite markers that yielded a clear amplification band on agarose gels were then tested for polymorphism level using 43–57 samples of A. hyacinthinus belonging to three potentially sub-populations of this species (Table 1). The polymorphism level was determined by a genotyping method throughout PCR amplification of the DNA samples with a sequence-specific forward primer composed by a M13 (−21) tail with distinct fluorescent dye labels (HEX or FAM, Applied Biosystems) , following a previously established protocol . Fragments genotyping was performed in ABIPrism 3500xL or ABIPrism 3130xL Genetic Analyzers (Applied Biosystems) in order to detect the alleles size. Each reaction contained 0.1 μL of TaqDNA polymerase (5 U/μL, Invitrogen), 0.1 μL of the forward primer with a M13 tail (5′CACGACGTTGTAAAACGAC-3′) (10 μM) , 0.3 μL of the reverse primer (10 μM), 0.2 μL of fluorescent dye label (10 mM), 0.4 μL of MgCl2 (25 mM), 1 μL of dNTPs (4 mM), 1.2 μL of PCR buffer (10×), 20–50 ng of DNA, and autoclaved Milli-Q water to complete 12 μL. PCR conditions were as follows: initial denaturation at 95 °C for 10 min, followed by 35 cycles at 95 °C for 1 min, 57–60 °C for 40 s (as detailed in Table 2), and 72 °C for 40 s, followed by 10 M13 tail cycles at 95 °C for 1 min, 55 °C for 40 s, 72 °C for 40 s, and a final extension step at 72 °C for 10 min. Approximately 2 μL of the products were electrophoresed through 1.5 % agarose gels to evaluate amplification success and to estimate their concentration. An aliquot of 1.0 μL of each amplified product was mixed with 0.2 μL of a molecular marker GeneScan ™ ROX ™-500 STANDARD (2 fmol, Applied Biosystem) and Hi-Di™ formamide (Applied Biosystems) to complete 10 μL of reaction. The mixture was analyzed in automatic sequencer ABIPrism 3500 (Applied Biosystems). The peaks and the sizes of the fragments for each allele were obtained by GeneMarker v2.6.3. The Genepop 4.2 program  was used to check the number of alleles per locus, linkage disequilibrium between loci, expected heterozygosity (He), observed heterozygosity (Ho), inbreeding coefficient (Fis) and Hardy–Weinberg equilibrium. The software Cervus was used to estimate the PIC (polymorphic information content) of each locus and the Micro-checker program  was used to assess the presence of null alleles, genotyping errors, and allele dropout. Additionally, we tested the transferability of the described primer sets in other seven Neotropical parrot species (Amazona guildingii, A. ochrocephala, Ara severus, A. macao, A. ararauna, A. chloropterus, and Anodorhynchus leari). For this purpose, PCR conditions were as follows: an initial denaturation step at 95 °C for 10 min, followed by 10 cycles at 95 °C for 1 min, 60 °C touchdown decreasing 0.5 °C at every cycle, and 72 °C for 40 s, followed by 25 cycles at 95 °C for 1 min, 55 °C for 40 s, 72 °C for 40 s and a final extension at 72 °C for 10 min. Each PCR contained 20–50 ng of DNA, 0.2 μL of TaqDNA polymerase (5 U/μL, Invitrogen), 1 μL of each primer (10 mM), 0.4 μL of MgCl2 (25 mM), 1.2 μL of PCR buffer (10×), 1 μL of dNTPs (2.5 mM), and autoclaved Milli-Q water to complete 12 μL. An aliquot of 3 μL of each product was electrophoresed in 1.5 % agarose gels to evaluate amplification success.
From a library composed by a total of 96 bacterial colonies, 52 clones (54.17 %) were recovered and sequenced. The majority of the generated sequences, however, were not suitable for primers design or possessed a low number of repeat motifs and therefore had to be discarded. From this filtering step, only six clones could be used for the primers design (11.53 %).
As outcome, with the exception of the AnH23 locus that was found to be monomorphic, five isolated microsatellite loci were polymorphic in A. hyacinthinus (Table 2). The number of alleles per locus varied from two (AnH6 and AnH33) to 12 (AnH34), the observed heterozygosity ranged from 0.175 (AnH10) to 0.721 (AnH34), with an average of 0.440 for all loci, and the expected heterozygosity per locus ranged from 0.193 (AnH10) to 0.818 (AnH34) with an average of 0.509 between all loci (Table 2). All analyzed polymorphic loci have shown no evidence of linkage disequilibrium (p > 0.05) and no locus was found in Hardy–Weinberg disequilibrium (p > 0.01). In the analysis carried out in the Micro-checker program , the presence of null alleles was not detected. In addition, among the five polymorphic loci, two (AnH6 and AnH33) were considered reasonably informative (PIC between 0.25 and 0.5) and two (AnH17 and AnH34) were highly informative (PIC > 0.5) (Table 2) . Furthermore, the designed primers for Hyacinth macaw generated amplification products for other seven species of Neotropical parrots. The primer sets AnH6, AnH33, and AnH17 led to amplification results for two, three, and five species, respectively, whereas both AnH10 and AnH34 primer sets were effective for six species (Table 3). Although the AnH23 locus was monomorphic in A. hyachintinus, primers for this microsatellite also generated amplification products for all the other tested parrots and their polymorphism should be further investigated. Curiously, primers for the AnH6 locus did not lead to amplification results in A. leari that is the phylogenetically closest species to Hyacinth macaw, probably due to specific mutations in the primers annealing regions. The described primers in this work have proven to be functional and can serve as important tools to determine the variability and the genetic structure of populations of Hyacinth macaw and other Neotropical parrots and may assist in in situ and ex situ conservation plans.
Availability of supporting data
The microsatellites data supporting the results of this article are available in GenBank at NCBI (http://www.ncbi.nlm.nih.gov/genbank/). Accession Numbers KP860340 to KP860345.
Sick H. Ornitologia Brasileira. 3rd ed. Rio de Janeiro: Editora Nova Fronteira; 1997.
Guedes NMR. Management and conservation of the large macaws in the wild. Ornitol Neotrop. 2004;15:279–83.
Guedes NMR. Reproductive biology of Hyacinth Macaw (Anodorhynchus hyacinthinus) in the Pantanal-MS, Brazil. Dissertation. Piracicaba: University of São Paulo; 1993.
Presti FT, Oliveira-Marques AR, Silva GF, Miyaki CY, Guedes NM. Notas sobre alguns aspectos da biologia da arara-azul (Anodorhynchus hyacinthinus) (Psittaciformes: Psittacidae) na região do Carajás, Pará. Atual Ornitol. 2009;151:4–7.
BirdLife International. Species factsheet: Anodorhynchus hyacinthinus. 2015. http://www.birdlife.org. Accessed 18 March 2015.
Martinez-Cruz B, Godoy JA, Negro JJ. Population genetics after fragmentation: the case of the endangered Spanish imperial eagle (Aquila adalberti). Mol Ecol. 2004;13:2243–55.
Mockford SW, Herman TB, Snyder M, Wright JM. Conservation genetics of Blanding’s turtle and its application in the identification of evolutionarily significant units. Conserv Genet. 2007;8(1):209–19.
Chapman DD, Pinhal D, Shivji MS. Tracking the fin trade: genetic stock identification in western Atlantic scalloped hammerhead sharks Sphyrna lewini. ESR. 2009;9:221–8.
Haag T, Santos AS, Sana DA, Morato RG, Cullen L Jr, Crawshaw PG Jr, De Angelo C, Di Bitetti MS, Salzano FM, Eizirik E. The effect of habitat fragmentation on the genetic structure of a top predator: loss of diversity and high differentiation among remnant populations of Atlantic Forest jaguars (Panthera onca). Mol Ecol. 2010;19(22):4906–21.
Oliveira EJ, Padua JG, Zucchi MI, Vencovsky R, Vieira MLC. Origin, evolution and genome distribution of microsatellites. Genet Mol Biol. 2006;29(2):294–307.
Manel S, Berthier P, Luikart G. Detecting wildlife poaching: identifying the origin of individuals with bayesian assignment tests and multilocus genotypes. Conserv Biol. 2002;16(6):650–9.
Berger-Wolf TY, Sheikh SI, DasGupta B, Ashley MV, Caballero IC, Chaovalitwongse W, Putrevu SL. Reconstructing sibling relationships in wild populations. Bioinformatics. 2007;23(13):i49–56.
Presti FT, Wasko A. A review of microsatellite markers and their application on genetic diversity studies in parrots. Open J Genet. 2014;4:69–77.
Presti FT, Oliveira-Marques AR, Caparroz R, Biondo C, Miyaki CY. Comparative analysis of microsatellite variability in five macaw species (Psittaciformes, Psittacidae): application for conservation. Genet Mol Biol. 2011;34:348–52.
Faria PJ, Guedes NMR, Yamashita C, Martuscelli P, Miyaki CY. Genetic variation and population structure of the endangered Hyacinth Macaw (Anodorhynchus hyacinthinus): implications for conservation. Biodivers Conserv. 2008;17:765–79.
Presti FT, Guedes NMR, Antas PTZ, Miyaki CY. Population genetic structure in hyacinth macaws (Anodorhynchus hyacinthinus) and identification of the probable origin of confiscated individuals. J Hered. 2015;106(S1):491–502.
Caparroz R, Miyaki CY, Baker AJ. Characterization of microsatellite loci in the Blue and gold Macaw, Ara ararauna (Psittaciformes: Aves). Mol Ecol Notes. 2003;10:1046–8.
Gebhardt KJ, Waits LP. Cross-species amplification and optimization of microsatellite markers for use in six neotropical parrots. Mol Ecol Res. 2008;4:835–9.
Russello MA, Calcagnotto D, DeSalle R, Amato G. Characterization of microsatellite loci in the endangered St. Vicent Parrot, Amazona guildingii. Mol Ecol Notes. 2001;1:13–5.
Russello MA, Lin K, Amato G, Caccone A. Additional microsatellite loci for endangered St. Vicent Parrot, Amazona guildingii. Conserv Genet. 2005;6(4):643–5.
Taylor TD, Parkin DT. Characterization of 12 microsatellite primer pairs for the African grey parrot, Psittacus erithacus and their conservation across the Psittaciformes. Mol Ecol Notes. 2007;7:163–7.
Tavares ES, Baker AJ, Pereira SL, Miyaki CY. Phylogenetic relationships and historical biogeography of neotropical parrots (Psittaciformes: Psittacidae: Arini) inferred from mitochondrial and nuclear DNA sequences. Syst Biol. 2006;55:454–70.
Tokita M, Kiyoshi T, Armstrong KN. Evolution of craniofacial novelty in parrots through developmental modularity and heterochromy. Evol Dev. 2007;9:590–601.
Galbusera P, van Dongen S, Matthysen E. Cross species amplification of microsatellite primers in passerine birds. Conserv Genet. 2000;1:163–8.
Zane L, Bargelloni L, Patarnello T. Strategies for microsatellite isolation: a review. Mol Ecol. 2002;11:1–16.
Primmer CR, Painter JN, Koskinen MT, Palo JU, Merila J. Factors affecting avian cross-species microsatellite amplification. J Avian Biol. 2005;36:348–60.
Bruford MW, Hanotte O, Brookfield JFY, Burke T. Single locus and multilocus DNA fingerprinting. In: Hoelzel CAR, editor. Molecular genetic analyses of populations: a pratical approach. New York: Oxford University Press; 1992. p. 225–69.
Birboim HC, Doly J. A rapid alkaline extraction procedure for screening recombinant plasmid DNA. Nucleic Acids. 1979;7(6):1513–23.
Billotte N, Lagoda PJL, Risterucci AM, Baurens FC. Microsatellite-enriched libraries: applied methodology for the development of SSR markers in tropical crops. Fruits. 1999;54:277–88.
Temnykh S, Lukashova A, Cartinhour S, DeClerck G, Lipovich L, McCouch S. Computational and experimental analysis of microsatellites in rice (Oryza sativa L.): frequency, length variation, transposon associations, and genetic marker potential. Genome Res. 2001;11(8):1441–52.
Untergasser A, Nijveen H, Rao X, Bisseling T, Geurts R, Leunissen JAM. Primer3Plus, an enhanced web interface to Primer3. Nucleic Acids Res. 2007;35(suppl 2):w71–4.
Schuelke M. An economic method for the fluorescent labeling of PCR fragments. Nat Biotechnol. 2000;18(2):233–4.
Boutin-Ganache I, Raposo M, Raymond M, Descepper CF. M13-tailed primers improve the readability and usability of microsatellite analyses performed with two different allele sizing methods. Biotechniques. 2001;31:1–3.
Raymond M, Rousset F. GENEPOP (version 1.2): population genetic software for exact tests and ecumenicism. J Hered. 1995;86:248–9.
Brookfield JFY. A simple new method for estimating null allele frequency from heterozygote deficiency. Mol Ecol. 1996;5(3):453–5.
Botstein D, White RL, Skolnick M, Davis RW. Construction of a genetic linkage map in man using restriction fragment length polymorphisms. Am J Hum Genet. 1980;32:314–31.
HS performed the construction of microsatellite-enriched genomic library, optimized PCR reactions, microsatellite genotyping, and data analysis and also drafted the manuscript. FP helped in genotyping, in the discussion of the results and revised drafts of the manuscript. AW contributed with materials, equipment and revised the manuscript. DP coordinated the study, helped in the discussion of the results and critically edited and revised all drafts of the manuscript. All authors read and approved the final manuscript.
This work had a financial support of VALE S.A., São Paulo Research Foundation (FAPESP), Coordenadoria de Aperfeiçoamento de Pessoal de Nível Superior (CAPES). The Instituto Chico Mendes de Proteção à Biodiversidade and Instituto Arara Azul supported field campaings and data collection. The authors thanks Laboratório de Genética e Evolução Animal from University of São Paulo (USP) by samples donation and Grace Ferreira Da Silva, Luiz Pereira Rodrigues (Mogno Meio Ambiente), Dezivaldo Ribeiro Nascimento, Fernando, Patrick Karassawa, Ludmila Conrado, Talita Aleixo Roberto de Almeida, Mayla Barbirato, Frederico Drumond Martins and Humberto de Souza Baldan for helping in field work.
The authors declare that they have no competing interests.