Evidence of population genetic structure in Ecuadorian Andean bears.

Wildlife conservation in Andean countries is a global priority because of the high levels of biodiversity and endemism. Historically, these countries have had limited resources to monitor wildlife (e.g., through genetic tools) and establish conservation programs. Focusing on the study and emblematic use of a few charismatic species has been a strategic approach to direct efforts for conservation and development planning. Consequently, the Andean bear is a flagship and umbrella species for highly biodiverse Andean countries like Ecuador. The few studies exploring the population genetics of this species have concluded that it has low genetic diversity and few units for conservation as populations appear to be well connected. However, these results might be attributed to ascertainment bias as studies have been performed with heterologous molecular markers. Here, using both mtDNA sequences and species-specific microsatellite markers, we show that Andean bears in Ecuador have population structure. Additionally, we found through the study of three Ecuadorian populations that the species might have a higher genetic diversity than we previously thought. These results could support the revision of research priorities, conservation, and planning strategies to improve connectivity for this species which occurs in crucial biodiversity hotspots.

An Annotated Draft Genome for the Andean Bear, Tremarctos ornatus.

The Andean bear is the only extant member of the Tremarctine subfamily and the only extant ursid species to inhabit South America. Here, we present an annotated de novo assembly of a nuclear genome from a captive-born female Andean bear, Mischief, generated using a combination of short and long DNA and RNA reads. Our final assembly has a length of 2.23 Gb, and a scaffold N50 of 21.12 Mb, contig N50 of 23.5 kb, and BUSCO score of 88%. The Andean bear genome will be a useful resource for exploring the complex phylogenetic history of extinct and extant bear species and for future population genetics studies of Andean bears.

Genetics of the Andean bear (Tremarctos ornatus; Ursidae, Carnivora) in Ecuador: when the Andean Cordilleras are not an Obstacle.

One of the top carnivores in the Andean mountains is the Andean bear (Tremarctos ornatus, Ursidae), the only bear in South America. This is a flagship and key umbrella species in Ecuador because its conservation has a positive impact on the conservation of many other species in the Andes. But to preserve, first one must know the genetic characteristics of a species, among other things. For this, we analyzed six mitochondrial genes and seven nuclear DNA microsatellites of 108 Andean bear specimens sampled throughout Ecuador. We adopted three strategies for analyzing the data: by Province, by Region (north vs south), and by Cordillera. Four main results were obtained. First, the mitochondrial genetic diversity levels were elevated, but there were no differences in genetic diversity by Province or by Cordillera. By Regions, southern Ecuador had higher genetic diversity levels than to northern Ecuador. The genetic diversity for the microsatellites was only medium for the Andean bear at this country. Second, there was clear and significant evidence of female population expansions, for the overall sample, by Province, Region, and Cordillera. This population expansion was determined to have occurred in the time interval of 30,000-20,000 years ago (YA), during the last phase of the Pleistocene. We detected a population decrease to have occurred more recently, within the last 5000 years. It continued until about 300-200 YA when a population increase was again detected. Third, there were, practically, no phylogeographic pattern nor genetic differentiation among Andean bear populations in Ecuador by Province or by Cordillera for either mitochondrial or microsatellite markers. There was a little more genetic differentiation between northern and southern areas. Fourth, there was no trace of significant spatial genetic structure for the Andean bear in Ecuador in agreement with the genetic differentiation analyses. This shows that the Andean Cordilleras in this country did not present an obstacle to the dispersion of this species. Therefore, all of the Andean bear specimens in Ecuador should be treated as a unique Management Unit (MU) for conservation purposes, differently to that determined for other countries as Colombia.

Analysis of the cytochrome c oxidase subunit II (COX2) gene in giant panda, Ailuropoda melanoleuca.

The giant panda, Ailuropoda melanoleuca (Ursidae), has a unique bamboo-based diet; however, this low-energy intake has been sufficient to maintain the metabolic processes of this species since the fourth ice age. As mitochondria are the main sites for energy metabolism in animals, the protein-coding genes involved in mitochondrial respiratory chains, particularly cytochrome c oxidase subunit II (COX2), which is the rate-limiting enzyme in electron transfer, could play an important role in giant panda metabolism. Therefore, the present study aimed to isolate, sequence, and analyze the COX2 DNA from individuals kept at the Giant Panda Protection and Research Center, China, and compare these sequences with those of the other Ursidae family members. Multiple sequence alignment showed that the COX2 gene had three point mutations that defined three haplotypes, with 60% of the sequences corresponding to haplotype I. The neutrality tests revealed that the COX2 gene was conserved throughout evolution, and the maximum likelihood phylogenetic analysis, using homologous sequences from other Ursidae species, showed clustering of the COX2 sequences of giant pandas, suggesting that this gene evolved differently in them.

Molecular cloning, overexpression, purification, and sequence analysis of the giant panda (Ailuropoda melanoleuca) ferritin light polypeptide.

The complementary DNA (cDNA) of the giant panda (Ailuropoda melanoleuca) ferritin light polypeptide (FTL) gene was successfully cloned using reverse transcription-polymerase chain reaction technology. We constructed a recombinant expression vector containing FTL cDNA and overexpressed it in Escherichia coli using pET28a plasmids. The expressed protein was then purified by nickel chelate affinity chromatography. The cloned cDNA fragment was 580 bp long and contained an open reading frame of 525 bp. The deduced protein sequence was composed of 175 amino acids and had an estimated molecular weight of 19.90 kDa, with an isoelectric point of 5.53. Topology prediction revealed one N-glycosylation site, two casein kinase II phosphorylation sites, one N-myristoylation site, two protein kinase C phosphorylation sites, and one cell attachment sequence. Alignment indicated that the nucleotide and deduced amino acid sequences are highly conserved across several mammals, including Homo sapiens, Cavia porcellus, Equus caballus, and Felis catus, among others. The FTL gene was readily expressed in E. coli, which gave rise to the accumulation of a polypeptide of the expected size (25.50 kDa, including an N-terminal polyhistidine tag).

Synergistic roles of climate warming and human occupation in Patagonian megafaunal extinctions during the Last Deglaciation.

The causes of Late Pleistocene megafaunal extinctions (60,000 to 11,650 years ago, hereafter 60 to 11.65 ka) remain contentious, with major phases coinciding with both human arrival and climate change around the world. The Americas provide a unique opportunity to disentangle these factors as human colonization took place over a narrow time frame (~15 to 14.6 ka) but during contrasting temperature trends across each continent. Unfortunately, limited data sets in South America have so far precluded detailed comparison. We analyze genetic and radiocarbon data from 89 and 71 Patagonian megafaunal bones, respectively, more than doubling the high-quality Pleistocene megafaunal radiocarbon data sets from the region. We identify a narrow megafaunal extinction phase 12,280 ± 110 years ago, some 1 to 3 thousand years after initial human presence in the area. Although humans arrived immediately prior to a cold phase, the Antarctic Cold Reversal stadial, megafaunal extinctions did not occur until the stadial finished and the subsequent warming phase commenced some 1 to 3 thousand years later. The increased resolution provided by the Patagonian material reveals that the sequence of climate and extinction events in North and South America were temporally inverted, but in both cases, megafaunal extinctions did not occur until human presence and climate warming coincided. Overall, metapopulation processes involving subpopulation connectivity on a continental scale appear to have been critical for megafaunal species survival of both climate change and human impacts.

Ancient mitochondrial DNA reveals convergent evolution of giant short-faced bears (Tremarctinae) in North and South America.

The Tremarctinae are a subfamily of bears endemic to the New World, including two of the largest terrestrial mammalian carnivores that have ever lived: the giant, short-faced bears Arctodus simus from North America and Arctotherium angustidens from South America (greater than or equal to 1000 kg). Arctotherium angustidens became extinct during the Early Pleistocene, whereas Arctodus simus went extinct at the very end of the Pleistocene. The only living tremarctine is the spectacled bear (Tremarctos ornatus), a largely herbivorous bear that is today only found in South America. The relationships among the spectacled bears (Tremarctos), South American short-faced bears (Arctotherium) and North American short-faced bears (Arctodus) remain uncertain. In this study, we sequenced a mitochondrial genome from an Arctotherium femur preserved in a Chilean cave. Our molecular phylogenetic analyses revealed that the South American short-faced bears were more closely related to the extant South American spectacled bear than to the North American short-faced bears. This result suggests striking convergent evolution of giant forms in the two groups of short-faced bears (Arctodus and Arctotherium), potentially as an adaptation to dominate competition for megafaunal carcasses.

Identification of conserved microRNAs in peripheral blood from giant panda: expression of mammary gland-related microRNAs during late pregnancy and early lactation.

The giant panda (Ailuropoda melanoleuca) is one of the world's most endangered mammals, and it has evolved several unusual biological and behavioral traits. During puberty, pregnancy, lactation, and involution, the mammary gland undergoes profound morphological and functional changes. A large number of microRNAs (miRNAs) have been identified to be involved in mammary gland development and lactation. In this study, we identified 202 conserved mature miRNAs, corresponding to 147 pre-miRNAs, in giant panda peripheral blood using a small RNA-sequencing approach. In addition, 27 miRNA families and 29 miRNA clusters were identified. We analyzed the arm selection preference of pre-miRNAs and found that: 1) most giant panda pre-miRNAs generated one-strand miRNAs, and the 5p-arm only miRNAs have a higher expression level than 3p-arm only miRNAs; 2) there were more 5p-arm dominant miRNAs than 3p-arm dominant miRNAs; and 3) 5p-arm dominant miRNAs have a larger fold change within miRNA pairs than 3p-arm dominant miRNAs. Expression of 12 lactation-related miRNAs was detected across late pregnancy and early lactation stages by qPCR, and seven miRNAs were identified as clustered in one significant model. Most of these clustered miRNAs exhibited inhibitory roles in proliferation and differentiation of mammary epithelial cells. Functional analysis highlighted important roles of the seven as signed miRNAs in mammary development and metabolic changes, including blood vessel morphogenesis, macromolecule biosynthesis, cell cycle regulation, and protein transport.

Estimation of the growth curve and heritability of the growth rate for giant panda (Ailuropoda melanoleuca) cubs.

Giant panda cubs have a low survival rate during the newborn and early growth stages. However, the growth and developmental parameters of giant panda cubs during the early lactation stage (from birth to 6 months) are not well known. We examined the growth and development of giant panda cubs by the Chapman growth curve model and estimated the heritability of the maximum growth rate at the early lactation stage. We found that 83 giant panda cubs reached their maximum growth rate at approximately 75-120 days after birth. The body weight of cubs at 75 days was 4285.99 g. Furthermore, we estimated that the heritability of the maximum growth rate was moderate (h(2) = 0.38). Our study describes the growth and development of giant panda cubs at the early lactation stage and provides valuable growth benchmarks. We anticipate that our results will be a starting point for more detailed research on increasing the survival rate of giant panda cubs. Feeding programs for giant panda cubs need further improvement.

Molecular cloning, characterization, and bioactivity analysis of interleukin 18 in giant panda (Ailuropoda melanoleuca).

Interleukin 18 (IL-18), as a member of IL-1 superfamily, is an important pleiotropic cytokine that modulates Th1 immune responses. In this report, we cloned and identified a homolog of IL-18 in giant panda (Ailuropoda melanoleuca) (designated as AmIL-18) from peripheral blood mononuclear cells stimulated with lipopolysaccharide. The open readin g frame of AmIL-18 cDNA is 579 bp encoding a deduced protein of 192 amino acids. AmIL-18 gDNA fragments contained 5 exons and 4 introns. The amino acid sequence of AmIL-18 shared 23.9 to 87.0% identity with other species. To evaluate the effects of AmIL-18 on the immune response, we expressed the recombinant AmIL-18 in Escherichia coli BL21 (DE3). The fusion protein PET-AmIL-18 was purified by nickel affinity column chromatography and verified by sodium dodecyl sulfate polyacrylamide gel electrophoresis and Western blot analysis. The biological function of purified PET-AmIL-18 was determined on mouse splenocytes by quantitative real-time polymerase chain reaction. INF-γ and other cytokines were increased when stimulated by PET-AmIL-18, particularly when combined with recombinant human interleukin 12, while a Th2-type cytokine, interleukin-4, was strikingly suppressed. These results will provide information for the potential use of recombinant proteins to manipulate the immune response in giant pandas and facilitate the study to protect this treasured species.

Comparative analysis and molecular characterization of genomic sequences and proteins of FABP4 and FABP5 from the giant panda (Ailuropoda melanoleuca).

Fatty acid binding proteins (FABPs) are a family of small, highly conserved cytoplasmic proteins that bind long-chain fatty acids and other hydrophobic ligands. In this study, cDNA and genomic sequences of FABP4 and FABP5 were cloned successfully from the giant panda (Ailuropoda melanoleuca) using reverse transcription polymerase chain reaction (RT-PCR) technology and touchdown-PCR. The cDNAs of FABP4 and FABP5 cloned from the giant panda were 400 and 413 bp in length, containing an open reading frame of 399 and 408 bp, encoding 132 and 135 amino acids, respectively. The genomic sequences of FABP4 and FABP5 were 3976 and 3962 bp, respectively, which each contained four exons and three introns. Sequence alignment indicated a high degree of homology with reported FABP sequences of other mammals at both the amino acid and DNA levels. Topology prediction revealed seven protein kinase C phosphorylation sites, two casein kinase II phosphorylation sites, two N-myristoylation sites, and one cytosolic fatty acid-binding protein signature in the FABP4 protein, and three N-glycosylation sites, three protein kinase C phosphorylation sites, one casein kinase II phosphorylation site, one N-myristoylation site, one amidation site, and one cytosolic fatty acid-binding protein signature in the FABP5 protein. The FABP4 and FABP5 genes were overexpressed in Escherichia coli BL21 and they produced the expected 16.8- and 17.0-kDa polypeptides. The results obtained in this study provide information for further in-depth research of this system, which has great value of both theoretical and practical significance.

Overexpression, purification, and pharmacologic evaluation of anticancer activity of ribosomal protein L24 from the giant panda (Ailuropoda melanoleuca).

The ribosomal protein L24 (RPL24) belongs to the L24E family of ribosomal proteins and is located in the cytoplasm. The purpose of this study was to investigate the structure and anti-cancer function of RPL24 of the giant panda (Ailuropoda melanoleuca). The complementary DNA of RPL24 was cloned successfully using reverse transcription-polymerase chain reaction technology. We constructed a recombinant expression vector containing RPL24 complementary DNA and overexpressed it in Escherichia coli using pET28a plasmids. The expression product obtained was purified using Ni-chelating affinity chromatography. The results indicated that the length of the fragment cloned is 509 bp, and it contains an open-reading frame of 474 bp encoding 157 amino acids. Primary structure analysis revealed that the molecular weight of the putative RPL24 protein is 17.78 kDa with a theoretical isoelectric point of 11.86. The RPL24 gene is readily expressed in E. coli, and the RPL24 fused with the N-terminal histidine-tagged protein to give rise to the accumulation of an expected 23.51-kDa polypeptide. The inhibitory rate in mice treated with 0.1 mg/mL RPL24, the highest of 3 doses administered, can reach 67.662%, which may be comparable to the response to mannatide. The histology of organs with tumors showed that the tissues in the RPL24 group displayed a looser arrangement compared with that in the control group. Furthermore, no obvious damage was apparent in other organs, such as heart, lung, and kidney. The data showed that the recombinant RPL24 had time and dose dependency on the cell growth inhibition rate. Human laryngeal carcinoma Hep-2 cells treated with 0.3125-10 µg/mL RPL24 for 24 h displayed significant cell growth inhibition (P < 0.05; N = 6) in assays using 3-(4,5-dimethylthiazol- 2-yl)-2,5-diphenyltetrazolium bromide compared with that in control (untreated) cells. By contrast, human hepatoma Hep G-2 cells displayed no significant change (P > 0.05; N = 6) from control (untreated) cells. RPL24 has time and dose dependency on Hep-2 cell growth inhibition. The data indicate that the effect at low concentrations is better than that at high concentrations, and the concentration of 0.625 µg/mL provides the best rate of growth inhibition. Further research is ongoing to determine the bioactive principles of recombinant RPL24 protein that are responsible for its anticancer activity.

Coexpression of interleukin-6 and -2 from giant panda in Escherichia coli and the biological activity of the fusion protein.

To construct a fusion cytokine protein with more and stronger bioactivities to enhance the immunity of the cytokine alone, we expressed interleukin (IL)-6/(IL)-2 from giant panda (Ailuropoda melanoleuca) in Escherichia coli as a 59.4-kDa fusion protein. Subsequently, the inclusion bodies were solubilized with 8 M urea and applied onto a Ni-nitrilotriacetic acid column. The final production of IL-6/IL-2 reached 6 mg/L in soluble form, and the purified final product was >96% pure. In Western blot assays, the recombinant IL-6/IL-2 was recognized by polyclonal antibodies against IL-6 and IL-2 of giant panda. The results demonstrated that the protein mixture contained correctly folded IL-2 and IL-6 proteins. A 3-(4,5-dimethylthiazol-2-yl) 2,5-diphenyltetrazolium bromide assay demonstrated that IL-6/IL-2 can promote lymphocyte proliferation and differentiation. These data suggest that the fusion protein could be used to develop a novel immunoadjuvant to enhance the immunity of animals against infectious diseases.

cDNA, genomic sequence cloning and overexpression of giant panda (Ailuropoda melanoleuca) mitochondrial ATP synthase ATP5G1.

The ATP5G1 gene is one of the three genes that encode mitochondrial ATP synthase subunit c of the proton channel. We cloned the cDNA and determined the genomic sequence of the ATP5G1 gene from the giant panda (Ailuropoda melanoleuca) using RT-PCR technology and touchdown-PCR, respectively. The cloned cDNA fragment contains an open reading frame of 411 bp encoding 136 amino acids; the length of the genomic sequence is of 1838 bp, containing three exons and two introns. Alignment analysis revealed that the nucleotide sequence and the deduced protein sequence are highly conserved compared to Homo sapiens, Mus musculus, Rattus norvegicus, Bos taurus, and Sus scrofa. The homologies for nucleotide sequences of the giant panda ATP5G1 to those of these species are 93.92, 92.21, 92.46, 93.67, and 92.46%, respectively, and the homologies for amino acid sequences are 90.44, 95.59, 93.38, 94.12, and 91.91%, respectively. Topology prediction showed that there is one protein kinase C phosphorylation site, one casein kinase II phosphorylation site, five N-myristoylation sites, and one ATP synthase c subunit signature in the ATP5G1 protein of the giant panda. The cDNA of ATP5G1 was transfected into Escherichia coli, and the ATP5G1 fused with the N-terminally GST-tagged protein gave rise to accumulation of an expected 40-kDa polypeptide, which had the characteristics of the predicted protein.

cDNA, genomic sequence cloning and overexpression of ribosomal protein gene L9 (rpL9) of the giant panda (Ailuropoda melanoleuca).

The ribosomal protein L9 (RPL9), a component of the large subunit of the ribosome, has an unusual structure, comprising two compact globular domains connected by an α-helix; it interacts with 23 S rRNA. To obtain information about rpL9 of Ailuropoda melanoleuca (the giant panda) we designed primers based on the known mammalian nucleotide sequence. RT-PCR and PCR strategies were employed to isolate cDNA and the rpL9 gene from A. melanoleuca; these were sequenced and analyzed. We overexpressed cDNA of the rpL9 gene in Escherichia coli BL21. The cloned cDNA fragment was 627 bp in length, containing an open reading frame of 579 bp. The deduced protein is composed of 192 amino acids, with an estimated molecular mass of 21.86 kDa and an isoelectric point of 10.36. The length of the genomic sequence is 3807 bp, including six exons and five introns. Based on alignment analysis, rpL9 has high similarity among species; we found 85% agreement of DNA and amino acid sequences with the other species that have been analyzed. Based on topology predictions, there are two N-glycosylation sites, five protein kinase C phosphorylation sites, one casein kinase II phosphorylation site, two tyrosine kinase phosphorylation sites, three N-myristoylation sites, one amidation site, and one ribosomal protein L6 signature 2 in the L9 protein of A. melanoleuca. The rpL9 gene can be readily expressed in E. coli; it fuses with the N-terminal GST-tagged protein, giving rise to the accumulation of an expected 26.51-kDa polypeptide, which is in good agreement with the predicted molecular weight. This expression product could be used for purification and further study of its function.

cDNA, genomic sequence cloning and overexpression of the ribosomal protein S13 gene in the giant panda (Ailuropoda melanoleuca).

The cDNA and the genomic sequence of ribosomal protein S13 (RPS13) of the giant panda (Ailuropoda melanoleuca) was cloned using reverse transcription-polymerase chain reaction (RT-PCR) and touchdown-PCR, respectively. These two sequences were sequenced and analyzed, and the cDNA of the RPS13 gene was overexpressed in Escherichia coli BL21. We compared the nucleotide sequences of the coding region and the amino acid sequences with those of seven other mammalian species retrieved from GenBank. The cDNA fragment of the RPS13 cloned from the giant panda is 496 bp in size, containing an open-reading frame of 456 bp, encoding 151 amino acids. The length of the genomic sequence is 2277 bp, with five exons and four introns. The coding sequence shows a high degree of homology to those of Homo sapiens, Bos taurus, Canis lupus familiaris, Macaca mulatta, Mus musculus, Rattus norvegicus, and Pan troglodytes; the degree of homology was 91.23, 94.30, 94.74, 92.11, 87.94, 87.72, and 91.45%, respectively. The homologies for the deduced amino acid sequences reached as high as 99%. Primary structure analysis revealed that the molecular weight of the putative RPS13 protein is 17.22325 kDa, with a theoretical pI of 10.42. Based on topology prediction, there is one protein kinase C phosphorylation site, one casein kinase II phosphorylation site, two N-myristoylation sites, and one ribosomal protein S15 signature in the RPS13 protein of the giant panda. The RPS13 gene can be expressed in E. coli and the RPS13 protein fused with the N-terminally GST-tagged form, which gave rise to the addition of an expected 43-kDa polypeptide.

Phylogeography and mitochondrial diversity of extirpated brown bear (Ursus arctos) populations in the contiguous United States and Mexico.

The fossil record indicates that the brown bear (Ursus arctos) colonized North America from Asia over 50 000 years ago. The species historically occupied the western United States and northern Mexico but has been extirpated from over 99% of this range in the last two centuries. To evaluate colonization hypotheses, subspecific classifications, and historical patterns and levels of genetic diversity in this region, we sequenced 229 nucleotides of the mitochondrial DNA control region in 108 museum specimens. The work was set in a global context by synthesizing all previous brown bear control region sequences from around the world. In mid-latitude North America a single moderately diverse clade is observed, represented by 23 haplotypes with up to 3.5% divergence. Only eight of 23 haplotypes (35%) are observed in the extensively sampled extant populations suggesting a substantial loss of genetic variability. The restriction of all haplotypes from mid-latitude North America to a single clade suggests that this region was founded by bears with a similar maternal ancestry. However, the levels and distributions of diversity also suggest that the colonizing population was not a small founder event, and that expansion occurred long enough ago for local mutations to accrue. Our data are consistent with recent genetic evidence that brown bears were south of the ice prior to the last glacial maximum. There is no support for previous subspecies designations, although bears of the southwestern United States may have had a distinctive, but recent, pattern of ancestry.

Microsatellite analysis of the spectacled bear (Tremarctos ornatus) across its range distribution.

DNA samples of the spectacled bear (Tremarctos ornatus) from five Andean countries, Venezuela, Colombia, Ecuador, Peru and Bolivia, were analyzed for nine microsatellite loci. Seven of them were polymorphic, which led us to investigate several population-genetic parameters. Private alleles and significant differences in gene frequencies were found among the populations studied, which demonstrated the extent of genetic differentiation among the spectacled bear populations. The levels of gene diversity measured with these microsatellites were rather modest in this species. Hardy-Weinberg disequilibrium was especially found for the overall and the Ecuadorian samples, and might be due to the Wahl-und effect or consanguinity. Significant genetic heterogeneity was mainly observed among the Colombian and the Ecuadorian populations. Markov chain Monte Carlo simulations clearly showed that two different gene pools were present, one present in the Venezuelan-Colombian bears and other in the Ecuadorian ones.

Molecular population genetic analysis of the spectacled bear (Tremarctos ornatus) in the northern Andean area.

Eighty-two Andean bear samples obtained in three South American Andean countries were analyzed using five hypervariable microsatellite markers. Neither the overall sample nor the samples coming from each of the countries analyzed were in Hardy-Weinberg equilibrium. This is attributable to the Wahlund effect caused by a fragmentation of the bear populations. The genetic variability levels found were in general low for this kind of molecular markers (H = 0.38 for the global population). The situation of the Ecuadorian bear population seems to be endangered (H = 0.24). This is the lowest level recorded for any bear population. The genetic heterogeneity among the three populations was large and highly significant (F(ST) = 0.39; R(ST) = 0.32-0.49). Likewise, the gene flow estimates were remarkably low among these populations (Nm = 0.2-0.3). This means that the populations are strongly isolated. Different methods were applied to determine the effective numbers in these populations. A significant spatial structure of the genotypes was analyzed with different procedures. In all cases, an isolation-by-distance structure was detected. This could be a consequence of the original South America colonization. There was no evidence for a recent genetic bottleneck. This could mean that the low heterozygosity and the population fragmentation could be explained by ancient events related to the bear colonization, with the arrival of humans beings in the Americas, and/or with the glacial maximum, 16.000-30.000 years ago.