1
Revista de Biología Tropical, ISSN: 2215-2075, Vol. 71(S4): e57278, diciembre 2023 (Publicado Nov. 01, 2023)
Initial characterization of mitochondrial DNA control region haplotypes
of the Antillean manatee (Trichechus manatus manatus,
Sirenia:Trichechidae) in Guatemala
Grecia Mendez*1; https://orcid.org/0000-0002-8559-6557
Susan Carney2; https://orcid.org/0000-0003-3047-7435
Heidy Amelia Garcia3; https://orcid.org/0000-0002-2392-6768
Ester Quintana-Rizzo4; https://orcid.org/0000-0002-8957-0506
1. Department of Biology, Universidad del Valle, Guatemala; grecia.m.mendez@gmail.com (Correspondence*)
2. Department of Biology, Hood College, Frederick MD, USA; slc239@gmail.com
3. Fundación Defensores de la Naturaleza, Guatemala, Guatemala; hgarcia@defensores.org.gt
4. Department of Biology, Emmanuel College, Boston MA, USA; tetequintana@comcast.net
Received 21-IX-2022. Corrected 13-II-2023. Accepted 11-IV-2023.
ABSTRACT
Introduction: Small populations are at risk of losing genetic variability much faster than large populations; this
subsequently decreases their ability to adapt when facing environmental changes. A small population of the
endangered Antillean manatee (Trichechus manatus manatus) has been identified in Guatemala.
Objective: This study explored the genetic diversity of the Antillean manatee in Guatemala by analysing mito-
chondrial DNA control region haplotypes in the two most important habitats for the species, Bahía La Graciosa,
a coastal bay and Bocas del Polochic, a coastal wetland, both located in the Izabal State.
Methods: Genetic samples were collected using non or minimally invasive sampling techniques: scraping of
epidermal tissue, collection of floating feces, and collection of tissue from carcasses. DNA extractions, DNA
amplification using polymerase chain reaction (PCR), and sequencing of the control D-loop region were used to
process and analyse the samples.
Results: Seven mitochondrial DNA sequences were obtained from 36 samples collected (minimum of four and
maximum of seven individuals). Four haplotypes were identified, A01, A03, A04, and J01. No other Central
American country has reported this number of haplotypes in a manatee population, and it is the first time that
haplotype A01 has been reported for the region. The Guatemalan manatee population comprises at least two
genetic lineages, the Florida/Greater Antilles lineage (haplotypes A01, A03, and A04) and the Mesoamerican
lineage (J01).
Conclusion: Further studies, with the use of nuclear markers, are necessary to understand the population
dynamics between Bahia La Graciosa and Bocas del Polochic to identify the number of management units pres-
ent in the country; also, the degree of relatedness with the Belizean population needs to be established to better
coordinate conservation efforts.
Key words: Non-invasive genetic sampling; endangered species; control D-loop region; Lago de Izabal; Atlantic
coast of Guatemala; conservation management plans.
RESUMEN
Caracterización inicial de haplotipos de la región control de ADN mitocondrial del Manatí Antillano
(Trichechus manatus manatus Sirenia:Trichechidae) en Guatemala
Introducción: Las poblaciones pequeñas corren el riesgo de perder variabilidad genética mucho más rápido que
una población de mayor tamaño; disminuyendo, así mismo, su capacidad de adaptarse ante cambios ambientales.
https://doi.org/10.15517/rev.biol.trop..v71iS4.57278
SUPPLEMENT • MANATEES
2Revista de Biología Tropical, ISSN: 2215-2075 Vol. 71(S4): e57278, diciembre 2023 (Publicado Nov. 01, 2023)
INTRODUCTION
Genetic studies on the West Indian mana-
tee, Trichechus manatus Linnaeus, 1758 have
shown low haplotype and nucleotide diversity,
with 29 haplotypes identified to date, distrib-
uted in three distinct lineages: Cluster I, Cluster
II, and Cluster III (Alvarez-Aleman et al., 2022;
Caballero et al., 2021; Díaz-Ferguson et al.,
2017; Garcia-Rodriguez et al., 1998; Satizabal et
al., 2012; Vianna et al., 2006). Cluster I includes
haplotypes found in Florida and the Greater
Antilles, Cluster II is distributed from the Gulf
of Mexico to the Caribbean coast of South
America, and Cluster III is exclusive to Brazil
and Guyana (Vianna et al., 2006). The Antillean
manatee (T. m. manatus), classified as endan-
gered on the IUCN Red List (Self-Sullivan &
Mignucci-Giannoni, 2008), is one of the two
subspecies of the West Indian manatee and is
found in small population pockets throughout
most of its range.
In Central America, the genetic variability
of the Antillean manatee has been studied in
Belize (Hunter et al., 2010) and Panama (Díaz-
Ferguson et al., 2017); these countries share the
Cluster II lineage (Díaz-Ferguson et al., 2017).
Belize has the largest manatee population in the
western hemisphere with around 1 000 animals
(Hunter et al., 2010) and, through conservation
efforts, could assist in the recovery of the adja-
cent Guatemalan population (Quintana-Rizzo
& Reynolds III, 2010). The most recent estimate
is 150 individuals in Guatemala (Quintana-
Rizzo & Reynolds III, 2010), and threats like
illegal hunting limit the survival of the popu-
lation (Machuca-Coronado & Corona, 2019).
The objective of this study was to explore the
genetic diversity of the Antillean manatee in
Guatemala in two areas that are ecologically
distinct but important for the species, Bocas
del Polochic, an inshore wetland and Bahia La
Graciosa, a coastal bay (Fig. 1A).
Manatee epidermal tissue, feces, and car-
cass skin cuts were collected at the two study
sites between July 2012 and February 2011. To
obtain epidermal tissue, a modified method,
initially developed by Carney et al. (2007) was
employed. A scraper with increased perfora-
tions (from 2 mm to 5 mm wide) and a lighter
pole (4 m of aluminum instead of 2 m of PVC)
were used to scrape the dorsal part of the
manatee without hurting the animal and subse-
quently storing the tissue with 75 % ethanol and
Una pequeña población de la especie en peligro de extinción el manatí antillano (Trichechus manatus manatus)
ha sido identificada en Guatemala.
Objetivo: Este estudio explora la diversidad genética del manatí antillano en Guatemala mediante el análisis de
haplotipos de la región de control del ADN mitocondrial en los dos hábitats más importantes identificados para
la especie, Bahía La Graciosa, un bahía costera y Bocas del Polochic, un humedal costero, ambos localizados en
el departamento de Izabal.
Métodos: Las muestras genéticas se colectaron utilizando técnicas de muestreo no invasivas o mínimamente
invasivas: raspado de tejido epidérmico, recolección de heces y recolección de tejido extraído de cadáveres. Se
usaron extracciones de ADN, amplificación de ADN médiate la reacción en cadena de la polimerasa (PCR) y
secuenciación de la región control D-loop.
Resultados: Se obtuvieron un total de siete secuencias de ADN mitocondrial de 36 muestras recolectadas. Cuatro
haplotipos fueron identificados, A01, A03, A04 y J01. Ningún otro país centroamericano ha reportado esta
cantidad de haplotipos en una población de manatíes y es la primera vez que se reporta el haplotipo A01 para
la región. La población de manatíes guatemaltecos comprende al menos dos linajes genéticos, el linaje Florida/
Antillas Mayores (haplotipos A01, A03 y A04) y el linaje mesoamericano (J01).
Conclusión: Son necesarios más estudios, con el uso de marcadores nucleares, para comprender la dinámica
poblacional entre Bahía La Graciosa y Bocas del Polochic y para poder identificar el número de unidades de
manejo presentes en el país: además, se debe establecer el grado de relación con la población de Belice para coor-
dinar mejor los esfuerzos de conservación.
Palabras clave: Muestro genético no invasivo; especie en peligro de extinción; región control D-loop; Lago de
Izabal; Costa Atlántica de Guatemala; planes de manejo de conservación.
3
Revista de Biología Tropical, ISSN: 2215-2075, Vol. 71(S4): e57278, diciembre 2023 (Publicado Nov. 01, 2023)
1 X TE (pH 8). After a manatee sighting, fecal
samples were collected within a 100 m radius
of the sighting. These were collected using
new disposable plastic bags and latex gloves
(Muschett et al., 2009) and were stored at room
temperature in sterile containers with 95 %
ethanol using a proportion of 1:3 (Muschett et
al., 2009). Lastly, skin cuts from dead carcasses
were collected opportunistically and stored at
-20 °C. Two extraction protocols were used to
treat tissue samples with different degradation
states: the phenol-chloroform method (Mus-
chett et al., 2009) was used for samples with
high DNA concentration and a silica method
(Höss & Pääbo, 1993) was used for low DNA
concentrations. Two extraction protocols were
also used to assay fecal samples: Zhang et al.
(2006) and Marrero et al. (2009). In the Zhang
et al. (2006) protocol, instead of binding the
DNA to a spin column with guanidine thiocya-
nate, the samples were resuspended in 500 µl
of isopropanol at -20 °C. A modification was
used to increase the purity index of the extrac-
tions to 1.6–2.0. This modification consisted of
applying the silica method (Marrero, 2009) to
the resuspended samples.
Mitochondrial control region D-loop DNA
was amplified by PCR using 1X Buffer (Pro-
mega), 4 mM of MgCl2, 150 µM of each dNTPs,
0.3 µM of each primer (CR-4 and CR-5), 1.5 U
of Taq DNA polymerase (Promega) and 1 µl
of extracted DNA with a concentration of 100
ng/µl in a total reaction volume of 25 µl. The
PCR cycling conditions were: 94 °C for 3 min,
followed by 35 cycles of 94 °C for 1 min, 47 °C
for 1 min, 72 °C for 1 min, and a final extension
Fig. 1. (A) Location of the two sampling sites and protected areas: Bocas del Polochic Wildlife Refuge in Lago de Izabal
and Bahía La Graciosa, part of the Punta de Manabique Wildlife Refuge, both located on the Atlantic coast of Guatemala.
(B) Manatee mitochondrial control region D-loop haplotypes found in Guatemala and other parts of the Americas and the
number of samples analyzed per country. Each haplotype is depicted in distinct colors and grouped into its corresponding
cluster, previously defined by Vianna et al. (2006). Data sources: Mexico, Florida and the Dominican Republic (Vianna et
al., 2006), Cuba (Alvarez-Alemán et al., 2022), Puerto Rico (Hunter et al., 2012), Colombia (Caballero et al., 2021), Panama
(Díaz-Ferguson et al., 2017), and Belize (Hunter et al., 2010). (C) Number of genetic haplotypes (HT) and sample size (N)
for Antillean manatee in Guatemala and other Central American countries (Belize and Panama). Since floating fecal samples
cannot be assigned to different manatees, the sample size of Guatemalan manatees in this study could be less than seven,
which is the number of sequences reported but not less than four because four different haplotypes were identified.
4Revista de Biología Tropical, ISSN: 2215-2075 Vol. 71(S4): e57278, diciembre 2023 (Publicado Nov. 01, 2023)
of 72 °C for 1 hr (García-Rodríguez et al.,
1998). The optimum annealing temperature
for fecal samples was 54 °C instead of 47 °C.
PCR products were run on a 1 % agarose gel
and visualized by electrophoresis to confirm
amplification of the expected ~ 410 bp size. The
dNTPs and primers were removed using the
Novagen purification kit. The PCR products
were bi-directionally sequenced (Macrogen
Inc., Korea) to ensure a high confidence level
in each nucleotide.
ChromasPro (2.1.10.1) was used to assem-
ble the forward and reverse strands into con-
tigs. A BLAST search of each contig was used
to identify the top match haplotype, previously
registered and defined by Vianna et al. (2006).
A Clustal Omega alignment was generated
between the contig and the haplotype, where
shared polymorphisms and percentage align-
ment were used to define each sample as its
respective haplotype.
A total of 36 manatee samples were col-
lected. Of these, 29 samples (27 fecal samples
and two tissue samples from carcasses) were
obtained in Bahía La Graciosa, and seven sam-
ples (two fecal samples and five tissue samples
from carcasses) were obtained in Bocas del
Polochic. No epidermal tissue samples were
collected during the study. Noise reduction
was a primary challenge for approaching
manatees since another research was simul-
taneously being carried out in the same boat.
This field experience limited the application of
this specific method in the region. However,
approaching manatees is delicate due to their
elusive behavior.
DNA was amplified from 22 of the samples
collected; all came from Bahía La Graciosa.
Seven samples were successfully sequenced;
one came from tissue samples and six from fecal
samples. The main reason for a lack of sequence
results from tissue samples was contamination;
the preservation, storage, and management of
samples was not optimal. On the other hand,
the limiting factor for processing fecal samples
was degradation; fecal samples degrade very
quickly when exposed to UV radiation.
Double coverage sequences ranged in
length from 390-444 bp. Two sequences had
read lengths of less than 410 bp, the length of
characterized haplotypes in Genbank. Howev-
er, these sequences of 390 and 392 bp in length
were 100 % matches to the J01 haplotype, the
most variable of the haplotype sequences iden-
tified here. To confirm each haplotype, forward
and reverse chromatograms for each manatee
sequence were aligned with the reference hap-
lotype from Genbank (Supplementary Mate-
rial 1). In all cases, there was a 100 % match
between the reference sequence and the chro-
matograms. Only one sequence can be con-
clusively traced back to an individual since it
consisted of a tissue sample from a carcass. The
remaining samples consisted of feces collected
exclusively in the bay area; therefore, there
is a risk that a single individual was sampled
multiple times and/or that these sequences
could only be representative of a few bay area
resident individuals.
Four haplotypes were identified, A01, A04,
A03 and J01 (Fig. 1B and Fig. 1C), with respec-
tive frequencies of 43 %, 14 %, 14 %, and 29 %.
These haplotypes have all been found in Cen-
tral America, except for haplotype A01 (Fig.
1B). All manatee sequences have been reg-
istered in GenBank (OQ587957-OQ587965).
Genetic diversity parameters like haplotype
frequency, haplotype diversity, and nucleotide
diversity could not be determined due to the
aforementioned risk.
This study reports four manatee mito-
chondrial control region D-loop haplotypes for
Guatemala, in a minimum of four individuals
and a maximum of seven manatees. This is
the highest number of haplotypes reported for
Mesoamerica; all other genetically studied pop-
ulations of the region, namely Belize, Panama
and the Caribbean coast of Mexico, have report-
ed three haplotypes (Díaz-Ferguson et al., 2017;
Hunter et al., 2010; Nourisson et al., 2011) (Fig.
1B). This significant finding in a relatively small
sample size strongly indicates that a more com-
prehensive population genetic analysis of the
Guatemalan manatees is warranted. Belize, the
adjacent manatee population to Guatemala, has
5
Revista de Biología Tropical, ISSN: 2215-2075, Vol. 71(S4): e57278, diciembre 2023 (Publicado Nov. 01, 2023)
reported the haplotypes A03, A04, and J01 in
113 individuals; (Fig. 1C) (Hunter et al., 2010).
These three haplotypes have all been identified
in Guatemala, along with the A01 haplotype
that was reported in three out of seven sequenc-
es. Panama has also reported a different set of
three manatee haplotypes: J01, J02, and J03
(Díaz-Ferguson et al., 2017).
Each haplotype identified in Guatemala
was matched with the corresponding cluster
previously named and classified by Vianna et al.
(2006) (Supplementary Material 2). These clus-
ters are lineages inferred using nucleotides and
sequential divergence parameters to construct
phylogenetic relationships. Haplotypes A01,
A03, and A04 belong to Cluster I (Florida and
Greater Antilles) and haplotype J01 belongs to
Cluster II (Mesoamerica) (Vianna et al., 2006).
The presence of Cluster I and Cluster II gives
a bimodal character to the manatee popula-
tion of Guatemala. This pattern has also been
detected in manatees sampled from Mexico,
Belize, Colombia, and Venezuela (Vianna et al.,
2006), indicating that the manatees in Guate-
mala likely share a common ancestry with the
rest of the Mesoamerican manatee population.
The presence of the A01 and the A03 hap-
lotypes in the study suggests that the Guatema-
lan manatees have a genetic relationship with
the North American and/or Greater Antilles
population (Fig. 1B). Haplotype A03 has been
reported in Cuba (Alvarez-Alemán et al., 2022)
and Belize (Hunter et al., 2010). A01 is the only
haplotype reported for the Florida population
and has a broad distribution across the wider
Antilles (Hunter et al., 2012; Vianna et al.,
2006). Florida manatees are known to travel
outside their normal range [e.g., from Florida
to Cuba, Alvarez-Aleman et al. (2010); from
Florida to Mexico, Castelblanco et al. (2021)];
interestingly, the haplotype A01 has not been
reported in the Belizean manatee population,
where 16 times as many samples have been
analyzed to date (Hunter et al., 2010).
A complete mitochondrial diversity study
and an analysis of nuclear DNA markers would
provide a further understanding of the origin
and genetic status of the Guatemalan manatee
population. Nuclear DNA markers would be
fundamental in determining the number of
management units present in the country. They
would aid in delineating population genetic
relationships with neighboring countries, like
Belize, and across the region. Conservation
efforts should focus on determining the distri-
bution of each genetic lineage within the coun-
try as well as the degree of interbreeding and,
prioritize the preservation of each ancestral
lineage and promote habitat/population con-
nectivity where possible.
Ethical statement: the authors declare that
they all agree with this publication and made
significant contributions; that there is no con-
flict of interest of any kind; and that we fol-
lowed all pertinent ethical and legal procedures
and requirements. All financial sources are fully
and clearly stated in the acknowledgments sec-
tion. A signed document has been filed in the
journal archives.
See supplementary material
a04v71s4-SM1
ACKNOWLEDGMENTS
Research was conducted under research
permit 042/2010 and collection permit I-019-
2010 approved by the Guatemalan Govern-
ment, National Council for Protected Areas
(CONAP). We would like to thank our main
sources of funding, FONACON “Fondo Nacio-
nal para la Conservación de la Naturaleza
and Defensores de la Naturaleza. We thank
the entire team of the project “Fortalecimien-
to Institucional para Consolidar el Manejo y
Conservación del Manatí (Trichechus manatus
manatus) en la Costa Atlántica de Guatemala,
Oscar Hugo Machuca, Arnoldo Chaal, and
Marco Tulio Milla who made the fieldwork
possible. To the Department of Biology of the
University, Universidad del Valle de Guatemala,
for making the laboratory part of the study a
6Revista de Biología Tropical, ISSN: 2215-2075 Vol. 71(S4): e57278, diciembre 2023 (Publicado Nov. 01, 2023)
success along with the careful mentorship from
Andres Avalos. We are eternally grateful for
everyones involvement in this project.
REFERENCES
Alvarez-Aleman, A., Hunter, M., Frazer, T. K., Powell, J.
A., Garcia Alfonso, E., & Austin, J. D. (2022). The
first assessment of the genetic diversity and struc-
ture of the endangered West Indian manatee in
Cuba. Genetica, 12, 327–341. https://doi.org/10.1007/
s10709-022-00172-8
Caballero, S., Ortiz-Giral, M. C., Bohorquez, L., Lozano
Mojica, J. D., Caicedo-Herrera, D., Arévalo-Gonzá-
lez, K., & Mignucci-Giannoni, A. A. (2021). Mito-
chondrial genetic diversity, population structure and
detection of Antillean and Amazonian manatees
in Colombia: new areas and new techniques. Fron-
tiers in Genetics, 12, 726916. https://doi.org/10.3389/
fgene.2021.726916
Carney, S. L., Bolen, E. E., Barton, S. L., Scolardi, K.
M., Englund, C. C., Tringali, M. D., & Reynolds
III, J. R. (2007). A minimally invasive method of
field sampling for genetic analyses of the Flori-
da Manatee (Trichechus manatus latirostris). Mari-
ne Mammal Science, 23, 967–975. https://doi.
org/10.1111/j.1748-7692.2007.00155.x
Castelblanco-Martínez, D. N., Alvarez-Aleman, A., Torres,
R., Teague, A. L., Barton, S. L., Rood, K. A., Ramos,
E. A., & Mignucci-Giannoni, A. A. (2021). First
documentation of long-distance travel by a Florida
manatee to the Mexican Caribbean. Ethology Ecology
and Evolution, 34, 545–556. https://doi.org/10.1080/0
3949370.2021.1967457
Díaz-Ferguson, E., Hunter, M., & Guzman, H. M. (2017).
Genetic composition and connectivity of the Anti-
llean Manatee (Trichechus manatus manatus) in Pana-
ma. Aquatic Mammals, 43, 378–383.
García-Rodríguez, A. I., Bowen, M., Domning, D. P.,
Mignucci-Giannoni, A. A., Marmontel, M., Montoya-
Ospina, R. A., Morales-Vela, B., Rudin, M., Bonde, R.
K., & McGuire, P. M. (1998). Phylogeny of West Indian
manatee (Trichechus manatus): how many populations
and how many taxa? Molecular Ecology, 9, 1137–1149.
https://doi.org/10.1046/j.1365-294x.1998.00430.x
Höss, M., & Pääbo, S. (1993). DNA extraction from Pleis-
tocene bones by a silica-based purification method.
Nucleic Acids Research, 21, 3913–3914. https://doi.
org/10.1093/nar/21.16.3913
Hunter, M. E, Auil-Gomez, N. E., Tucker, K. P., Bonde,
R. K., Powell, J., & McGuire, P. M. (2010). Low
genetic variation and evidence of limited disper-
sal in the regionally important Belize manatee.
Animal Conservation, 13, 592–602. https://doi.
org/10.1111/j.1469-1795.2010.00383.x
Hunter, M. E., Mignucci-Giannoni, A. A., Pause Tucker,
K., King, T. L., Bonde, R. K., Gray, B., & McGuire,
P. M. (2012). Puerto Rico and Florida manatees
represent genetically distinct groups. Conservation
Genetics, 13, 1623–1635. https://doi.org/10.1007/
s10592-012-0414-2
Machuca-Coronado, O., & Corona, M. F. (2019). El mana-
tí antillano Trichechus manatus manatus (Sirenia-
Trichechidae) en Guatemala: amenazas y procesos de
conservación. In Kraker C., Calderon A. P., Cabrera A.
A. (Eds). Research Perspectives on the Wild Mammals
of Guatemala (pp. 188–200). Asociación Guatemalte-
ca de Mastozoólogos.
Marrero, P., Fregel, R., Cabrera, V. M., & Nogales, M. (2009).
Extraction of high-quality host DNA from feces and
regurgitated seeds: A useful tool for vertebrate ecolo-
gical studies. Biological Research, 42, 147–151. http://
dx.doi.org/10.4067/S0716-97602009000200002
Muschett, G., Bonacic, C., & Vianna, J. A. (2009). A
noninvasive sampling method for genetic analyses
of the West Indian manatee (Trichechus manatus).
Marine Mammals Science, 25, 955–963. https://doi.
org/10.1111/j.1748-7692.2009.00310.x
Nourisson, C., Morales-Vela, B., Padilla-Saldívar, J., Tucker,
K. P., Clark, A., Olivera-Gomez, L. D., Bonde, R., &
McGuire, P. (2011). Evidence of two genetic clusters
of manatees with low genetic diversity in Mexico and
implications for their conservation. Genetica,139,
833–842. https://doi.org/10.1007/s10709-011-9583-z
Quintana-Rizzo, E., & Reynolds III, J. E. (2010). Regio-
nal Management Plan for the West Indian manatee
(Trichechus manatus) [Technical Report No. 48].
United Nations Environmental Program (CEP) and
Caribbean Environment Programme.
Satizábal, P., Mignucci-Giannoni, A. A., Duchêne, S., Caice-
do-Herrera, D., Perea-Sicchar, C. M., García-Dávila,
C. R.Trujillo, F., & Caballero, S. J. (2012). Phylogeo-
graphy and sex-biased dispersal across riverine mana-
tee populations (Trichechus inunguis and Trichechus
manatus) in South America. PLoS one, 7, e52468.
https://doi.org/10.1371/journal.pone.0052468
Self-Sullivan, C., & Mignucci-Giannoni, A. A. (2008). Tri-
chechus manatus ssp. manatus. The IUCN Red List of
Threatened Species. https://doi.org/10.2305/IUCN.
UK.2008.RLTS.T22105A9359161.en.
Vianna, A. J., Bonde, R. K., Caballero, S., Giraldos, J. P.,
Lima, R. P., Clark, A., Marmontel, M., Morales-Vela,
B., Souzza, M. J., Parr, L., Rodríguez-Lopez, M. A.,
Mignucci-Giannoni, A. A., Powell, J. A., & Santos, F.
R. (2006). Phylogeography, phylogeny and hybridiza-
tion in trichechid sirenians: implication for manatee
conservation. Molecular Ecology, 15, 433–447. https://
doi.org/10.1111/j.1365-294X.2005.02771.x
Zhang, B. W., Li, M., Ma, L. C., & Wei, F. W. (2006). A
widely applicable protocol for DNA isolation from
fecal samples. Biochemical Genetics, 44, 503–512.
https://doi.org/10.1007/s10528-006-9050-1