Abstract
Advances in high-throughput sequencing (HTS) have made it a powerful resource for the conservation of threatened species, providing information at both population and individual levels to inform management decisions. In South America, however, the application of HTS in conservation has been limited, primarily due to challenges in funding and access to advanced genomic equipment and analytical expertise. Darwin’s fox (Lycalopex fulvipes), endemic to Chile’s Valdivian Temperate Rainforest, is the most endangered canid in South America with a small and declining population estimated at less than 1000 mature individuals. Despite its endangered status, significant knowledge gaps remain. Here we highlight the potential of HTS to address these challenges, such as clarifying its taxonomy, demographic history, geographic distribution, population structure, genetic diversity, and pathogen exposure. Integrating molecular data into conservation planning will be pivotal in ensuring the long-term survival of Darwin’s fox by identifying priorities for targeted management interventions, highlighting areas of critical habitat for conservation, and guiding genetic rescue efforts to enhance genetic diversity and resilience.
1 Introduction
The Chilean endemic Darwin’s fox (Lycalopex fulvipes) is the most endangered canid in South America (). This small, solitary, omnivorous species is obligate to forest habitats and primarily confined to the dense understory of the Valdivian Temperate Rainforest in southern Chile, which is recognized as a biodiversity hotspot threatened by unsustainable commercial logging and large-scale deforestation (). Darwin’s fox populations persist in native forest remnants within the Nahuelbuta mountain range, where fewer than 100 mature individuals remain (), on Chiloé Island, home to fewer than 500 mature individuals (), as well as in the Valdivian coastal range (Vilà et al., 2004; ), and Gorbea () dominated by agricultural land and some remaining native forest. These populations exhibit slight ecological, behavioural and phenotypic differences. Mainland foxes primarily inhabit dense forests and are predominantly nocturnal, while Chiloé foxes are more habitat-flexible, exhibiting coastal foraging and more diurnal activity. Additionally, Chiloé individuals are slightly smaller on average ().
Additional unpublished sightings raise the question about the true distribution and the existence of undiscovered populations. Beyond habitat loss and human-induced disturbances, Darwin’s fox faces significant threats from feral and free-ranging domestic dogs which attack them (), disrupt their behaviour (), and expose them to pathogens, posing the risk of disease spillover (; ).
Until the mid 1990’s, Darwin’s fox was considered a subspecies of the South American grey fox (Lycalopex griseus). However, the use of mitochondrial DNA (mtDNA) markers led to its classification as a distinct species (Yahnke et al., 1996), though further studies disagree on its phylogenetic position within the genus Lycalopex (; ; ).
Despite its Endangered status on the IUCN Red List of Threatened Species () and in Chilean legislation (DS 151/2007 MINSEGPRES), there has been no comprehensive assessment of population structure or consistent evaluation of intraspecific variation, and genetic monitoring for the species remains absent. However, studies have suggested that genetic diversity among Darwin’s foxes on Chiloé Island is lower compared to their mainland counterparts in Nahuelbuta (Yahnke et al., 1996; Vilà et al., 2004; ; ). Additionally, Darwin’s fox exhibits extremely low genome-wide heterozygosity, with a significant portion of its autosomal genome characterised by extensive runs of homozygosity (ROH) ().
In September 2023, the “Plan for the Recovery, Conservation and Management of Darwin’s fox” () was submitted to the Chilean Ministry of Environment. This recovery plan underscores the importance of closing key knowledge gaps about Darwin’s fox, which are critical for shaping effective conservation management strategies and actions. Molecular data are expected to play a pivotal role in their success, as genetic markers can provide insights at the population or individual level that are otherwise difficult to obtain. Some major knowledge gaps that still need to be addressed include: (i) What is the phylogenetic position of Darwin’s fox within the genus Lycalopex? (ii) What is the evolutionary and phylogeographic history of the species? (iii) What is the current distribution of Darwin’s foxes? (iv) To what extent are remnant populations connected? (v) How extensive and widespread is inbreeding? (vi) What pathogens are Darwin’s foxes exposed to? While some of these questions can be answered using traditional methods (e.g. camera trap surveillance, parasite egg counts from faecal samples, non-invasive sample screening), many can only be adequately addressed through molecular approaches, particularly through high-throughput sequencing (HTS) techniques.
Here we focus on how HTS approaches can help to address critical, immediate, and conservation-relevant issues. Other topics that can be studied using the same or similar techniques fall outside the scope of this review.
2 Use of high-throughput sequencing in defining strategies for Darwin’s fox conservation management
The main HTS technologies are provided by Illumina, Pacific BioScience (PacBio), and Oxford Nanopore Technologies (ONT). Illumina platforms generate short, high-accuracy sequences ranging from 50 to 300 base pairs (bp) in length, either as single or paired-end reads. These sequences are applicable to a wide range of experimental designs, from whole genome sequencing (WGS) to metagenomics. In contrast, PacBio platforms produce long reads with an average length of 20 kilobases (kb), which are advantageous for resolving complex genomic regions and detecting structural variants. ONT can produce even longer reads, with some kits capable of generating sequences exceeding 50 kb, but has lower base-calling accuracy compared with PacBio. Additionally, the portability of certain ONT devices makes this technology suitable for field-based, on-site sequencing.
Three main HTS approaches can be followed: sequencing of the whole genome (whole-genome sequencing, WGS), sequencing only parts of the genome (reduced representation approach, RRA), or sequencing environmental or invertebrate-derived DNA (eDNA/iDNA).
WGS (Table 1) provides complete genetic information of a specimen and thus unravels its complete genetic landscape, including genetic diversity, inbreeding levels, evolutionary and demographic history, and even gene-environment associations via whole genome bisulfite sequencing (WGBS), which is useful to detect epigenetically modified (methylated) sites. To date, only two Darwin’s foxes have been sequenced through WGS (), and a chromosome level assembly has not yet been generated.
RRAs such as RNA sequencing (RNA-seq), Restriction Site Associated DNA Sequencing (RAD-seq), Targeted Capture and SNP Arrays (Table 1) retrieve sequence information from a subset of the genome, utilising methods designed to obtain/target specific regions of the genome, which enables comparison amongst samples. These methods allow the cost-effective study of genetic diversity within and between populations.
Environmental genomics (eDNA/iDNA) utilizes genetic material shed by organisms into their surroundings, such as water, soil, or air, or from blood consumed by invertebrates (“invertebrate-derived DNA”) (). These methods enable non-invasive species detection, biodiversity assessment, and population monitoring.
Table 1
| Technique | Advantages | Disadvantages | Sample type | Phylogenetics | Hybridization | Population structure | Genetic diversity | Historical demography | Genetic load | Example cases |
|---|---|---|---|---|---|---|---|---|---|---|
| Whole-Genome Sequencing (WGS) Encompasses the entire genome, including coding and non-coding regions, regulatory sequences, repetitive elements, and structural variation. | • Provides the highest resolution for identifying genetic variation and population structure. • Detects most types of variants including rare and novel mutations. • Enables precise detection of ROH for assessing inbreeding. • Facilitates integration with other omics data, e.g. transcriptomics, epigenomics. | • High cost compared to targeted or RRA methods, especially for population-level studies. • Generates massive amounts of data, requiring significant computational resources for storage, processing, and analysis. • Degraded DNA may result in lower coverage or increased error rates. | • Fresh tissue; provides high-quality, high-quantity DNA. • Non-invasive; e.g. hair with roots, feathers with quill, faeces. May contain contaminants that use up sequencing real estate. • Ancient and museum samples; e.g. bones, teeth, skin, and preserved tissue. Deeper sequencing required. | yes | yes | yes | yes | yes | yes | () WGS addressing taxonomy, hybridization, genetic diversity, and historical demography of canids in South America. () WGS to investigate the effect of migration, inbreeding and genetic load on juvenile survival in arctic foxes. |
| Restriction Site Associated DNA Sequencing (RAD-seq) Sequences regions adjacent to restriction enzyme cutting sites, reducing genome complexity. Multiple variants exist, differentiated by number and use of restriction enzymes (e.g. ddRAD, 3RAD). | • No reference genome needed, making it suitable for non-model organisms. • Low cost per sample. Cost-effective for large-scale studies compared to whole-genome sequencing. • Resolution (i.e. SNP-density) can be fine-tuned with choice of restriction enzymes. | • Incomplete genome coverage (SNP detection limited to regions near restriction sites). • Underestimates genetic diversity. • Restriction enzyme cut site dependent. • Susceptible to allele/locus dropout, relevant for studies of multiple taxa. | • Fresh tissue; provides high-quality, high-quantity DNA. • Non-invasive; e.g. hair with roots, feathers with quill. Works for low-to-medium coverage RAD-seq. • Hybridization and historical demography analyses limited to relatively recent events. | yes | yes | yes | yes | yes | no | (von Holdt, 2022) RAD-seq to quantify red wolf ancestry in coyotes (hybridization), determine population structure, and genetic diversity. () RAD-seq to study sea turtle hybridization, develop a hybrid-index and evaluate reproductive output of hybrids. |
| RNA Sequencing (RNA-seq) Sequences complementary DNA (cDNA) synthesized from RNA, providing data on transcribed sequences and their abundance; includes coding and non-coding RNA. | • No reference genome needed, benefits from annotated assembly. • Reveals gene expression, functional variation, transcript isoforms, mutations, and regulatory elements. • Provides detailed information about gene expression patterns specific to tissue/cell type or ‘treatment’. | • Costly for low abundance transcripts. • Limited to transcribed sequences. • RNA is unstable and prone to rapid degradation, making it challenging for studies of rare and elusive species. • Tissue and situational specificity does not represent the whole organism or population. | • Fresh tissue; provides high-quality RNA from specific tissues (e.g., blood, liver, brain). Tissue must match the study’s focus. • Requires immediate stabilization with RNAlater, flash freezing in liquid nitrogen, or storage at −80°C to prevent RNA degradation. | yes | yes | no | yes | no | yes | () RNA-seq for variant detection without a reference genome, to identify population structure and loci under selection in white-footed mice. () RNA-seq to detect adaptive evolution of immune-related genes of wolve´s blood transcriptome. |
| Targeted Capture Selectively enriches target region(s). Short DNA/RNA probes complementary to target regions are used to hybridize with target loci, while off-target DNA is washed away. | • No reference genome required. Probes can be generated from RRA approaches (e.g. RAD-seq loci) or related taxa. • Effective for degraded, fragmented or contaminated DNA. • High consistency among samples; sequence variation does not lead to allele/locus dropout. | • Generally requires a priori knowledge of target loci for probe design. • Can be costly for small projects (due to cost of probes). | • Fresh tissue; provides high-quality, high-quantity DNA. • Non-invasive; e.g. hair with roots, feathers with quill, faeces. • Ancient and museum samples; e.g. bones, teeth, skin. • Formalin fixed paraffin embedded. • eDNA/iDNA | yes | yes | yes | yes | no | yes | () Cross-species capture to generate sequence data to design a SNP-based monitoring tool for non-invasively collected samples. () Reconstructed ancestral sequences used to design probes for capture of divergent taxa without available reference. |
| SNP Arrays Detection of predefined SNPs by hybridising DNA to a microarray chip containing allele-specific oligonucleotide probes. | • High reproducibility and consistency across samples and studies. • Arrays designed for model species can be utilized for studies of related non-model taxa. • Highly specific and cost effective. • No reference genome needed. | • Limited to predefined variants, missing novel or rare SNPs and structural variants. • SNP selection for array design requires a priori information. • Severely degraded samples may fail, resulting in allele/locus dropout. | • Fresh tissue; provides high-quality, high-quantity DNA. • Non-invasive; e.g. hair with roots, feathers with quill, faeces. | yes | yes | yes | yes | no | no | (vonHoldt et al., 2013) Utilization of dog SNP array to develop SNP-panel for species identification and detection of hybridization. (vonHoldt et al., 2011) Uses dog SNP array to study evolutionary relationship among wolf-like canids. |
| Environmental DNA (eDNA), Invertebrate DNA (iDNA) Utilises DNA from environ-mental sources (eDNA), or from invertebrates (iDNA) that act as “DNA collectors” from vertebrates. Follows either a metagenomics or metabarcoding approach. | • Eliminates the need for direct sampling of organisms, making it ideal for endangered or elusive species. • Potentially detects all species in the DNA pool from the extracted sample, enabling ecosystem-wide monitoring. • Low cost for surveillance compared to non-HTS methods (e.g. camera trapping). | • DNA is often degraded and present in small amounts. • Risk of contamination during collection and processing. • Relies on reference databases, which may be incomplete for non-model species. • Provides less detailed genetic information than direct sampling. | • Environmental; e.g. water, soil, sediment, surfaces. Requires DNA extraction tailored to inhibitors (e.g. humic acids for soil/sediment). • Invertebrate derived; e.g. leeches, mosquitoes, ticks. May require pooling of samples from a locality to obtain sufficient material. | yes | no | no | no | no | no | () Combines eDNA with hybrid capture to identify mammal species from water samples in Namibia and Tanzania. () Utilization of eDNA from water and soil samples for indirect monitoring of maned wolf’s distribution in Argentina. |
Brief summary of high-throughput sequencing techniques and their potential usage for conservation related analyses.
3 Taxonomic uncertainty
Poor taxonomic assessment can lead to species misidentification, misallocation of resources, and ineffective protection measures, ultimately hindering conservation efforts (). While the species status of Darwin’s fox is undisputed, its precise phylogenetic placement within the genus Lycalopex remains unresolved. The rapid divergence of Lycalopex taxa began 1.3 million years ago (Mya), and started between 0.7 and 0.27 Mya for Darwin’s fox (Yahnke et al., 1996; ; ; ), complicating phylogenetic reconstruction. This is primarily due to the retention of ancestral polymorphisms (i.e. incomplete lineage sorting) and hybridization among various Lycalopex species (; ; ; ). Recent studies examining taxonomic relationships within the Lycalopex genus using various genetic markers (Yahnke et al., 1996; Vilà et al., 2004; ; ; ; ), and whole genomes () disagree on species relationships within Lycalopex (Figure 1).
Figure 1
Darwin’s fox is sympatric with two other Lycalopex species (L. culpaeus and L. griseus) in parts of its range, but it remains unknown whether hybridization occurs among these taxa. From a conservation perspective, hybridization can be a double-edged sword: it may threaten endangered species by diluting their gene pool (genetic swamping) or, conversely, increase genetic diversity and adaptability, thereby enhancing resilience (
WGS is likely the most effective method for reconstructing the Lycalopex species tree, as it enables a comprehensive evaluation of phylogenetic incongruences by sampling across both coding and non-coding regions, detecting rare variants, and (potentially) incorporating structural variants in analyses (
Broad sampling across the full geographic distribution of all Lycalopex species is advisable, as introgression may be geographically localized or restricted to specific lineages. Ideally, samples with uncertain provenance or heritage, such as those from zoos, should be avoided to prevent confounding results.
4 Uncertain distribution, abundance and connectivity
The extent of Darwin’s fox geographic distribution is unknown (Figure 1), as are population numbers, their sizes, and potential connectivity among populations, rendering effective conservation measures difficult. Underestimating the species’ range risks neglecting populations that could serve as critical genetic reservoirs, leading to a loss of genetic diversity. Maintaining that diversity, however, is essential for the species’ adaptive potential and long-term survival. The oldest documented and the most thoroughly studied population, first recorded by Darwin in 1840, is located on Chiloé Island, which represents the southernmost edge of the known distribution range of the species. The northernmost population resides in the Nahuelbuta area (
Detecting and surveying an elusive species in dense rainforests using traditional methods can be costly, logistically challenging and time-intensive. Non-invasive approaches, such as use of eDNA (
Promising areas for eDNA/iDNA surveillance include regions predicted as suitable habitat by niche modelling (
5 Genetic diversity and population structure
Surveying Darwin’s foxes’ intraspecific genetic variation is essential for evaluating population structure, genetic differentiation, isolation times, and detecting bottlenecks or signs of genomic erosion. It helps determine whether geographical distances or barriers contribute to genetic divergence among populations. Preserving remaining genetic diversity is critical, as signs of inbreeding are already present in Darwin’s foxes (
Early research on Darwin’s fox genetic variability focused on mtDNA control-region sequences, revealing that foxes on Chiloé Island shared the same haplotype, while mainland foxes had distinct haplotypes, suggesting differentiation between these populations (Yahnke et al., 1996; Vilà et al., 2004;
Although these pioneering studies have provided valuable insights, their limited sample and/or marker numbers limit the generalizability of the findings regarding genetic variability. Future research should aim for broader sampling across and within populations to better capture intraspecific variation. WGS offers the most comprehensive data on genetic diversity (e.g. SNPs, indels, runs of homozygosity, and structural variants), enabling detailed analyses of population structure, demography, connectivity, kinship, divergence times, and more (
An important factor in selecting a HTS approach is its ability to assess genetic load, which is vital for understanding inbreeding depression and the population dynamics of deleterious alleles. Genetic load arises from the accumulation of harmful variants that reduce fitness by increasing expression of recessive deleterious alleles and potentially fixing them through genetic drift; expression of these harmful variants can negatively affect health, adaptability, and reproduction (
6 Pathogen community and adaptive immune system
Diseases are an important, yet often underestimated factor influencing species demography, particularly when new pathogens are introduced into naive populations. Dogs roaming close to and within protected areas of Darwin’s fox distribution range, unvaccinated and untreated for parasites (
Several bacterial pathogens have been detected in Darwin’s foxes, including Toxoplasma gondii, Leptospira sp., Mycoplasma haemocanis (
HTS is revolutionising pathogen detection by enabling broad-spectrum analysis, integrating data from hosts, vectors, and environmental samples to provide a comprehensive understanding of pathogen transmission (
Genomic regions like the Major Histocompatibility Complex (MHC), or Dog Leukocyte Antigen (DLA) in canids, are essential for adaptive immunity and pathogen response (Yuhki et al., 2007). Reduced diversity at these loci can increase susceptibility to disease, while introgressive hybridization with other canids may enhance variation and resilience to pathogens. The uncharacterized diversity of DLA genes in Darwin’s foxes raises uncertainty about whether populations exhibit reduced variation at these loci, which is detrimental for developing strategies to maintain functional diversity and enhance the species’ resilience to disease (
7 Discussion and future perspectives
High-throughput sequencing has emerged as an increasingly valuable tool for conservationists, gaining prominence in wildlife management due to its enhanced accessibility and effectiveness. It generates highly informative data that can support critical decision-making, as evidenced by efforts to conserve species on the brink of extinction, such as the Iberian lynx (
The integration of genomic data into conservation decision-making is indispensable in shaping effective conservation policies, allocation of resources and designing management plans, like the one proposed to Chilean authorities by the N.G.O.
Finally, as a flagship species for the Valdivian Temperate Rainforest, efforts to protect Darwin’s foxes would also contribute to the conservation of the whole biodiversity of this unique ecosystem, which harbours many threatened and emblematic species, such as the Pudu (Pudu puda), Southern River Otter (Lontra provocax), Darwin’s frog (Rhinoderma darwinii), Monito del monte (Dromiciops gliroides), Alerce tree (Fitzroya cupressoides) and the Long-nosed shrew opossum (Rhyncholestes raphanurus).
Statements
Author contributions
CV: Conceptualization, Visualization, Writing – original draft, Writing – review & editing. JG: Writing – original draft, Writing – review & editing. JF: Supervision, Writing – review & editing. DF: Conceptualization, Supervision, Writing – original draft, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by the Leibniz Institut für Zoo und Wildtierforschung (IZW), as well as through the scholarship grant provided by the Agencia Nacional de Investigación y Desarrollo (ANID) (grant n° 62210037) and the Deutscher Akademischer Austauschdienst (DAAD) (grant n° 91825331). Open Access funding was provided by the IZW's Institutional partnership with Frontiers within the Open Access Publishing Framework Agreement.
Acknowledgments
We would like to thank Dr. Juliana Vianna, Dr. Francisco Fonturbel, Felipe Osorio, Juan José Saez, and Eduardo Pizarro for their valuable feedback on early drafts of the manuscript.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcosc.2025.1512531/full#supplementary-material
References
1
AbramsJ. F.HörigL. A.BrozovicR.AxtnerJ.Crampton-PlattA.MohamedA.et al. (2019). Shifting up a gear with iDNA: From mammal detection events to standardised surveys. J. Appl. Ecol.56, 1637–1648. doi: 10.1111/jpe.2019.56.issue-7
2
Acosta-JamettG.ContrerasS.MuñozP.BriceñoC.ChirgwinC.HernándezF. (2018). Description of gastrointestinal parasitism through coprologic survey in Darwin’s fox, Lycalopex fulvipes (Martin 1837), and kodkod, Leopardus guigna (Molina 1782), in Chiloé island, Chile. Gayana82, 160–165. doi: 10.4067/S0717-65382018000200160
3
Acosta-JamettG.SurotD.CortésM.MarambioV.ValenzuelaC.VallverduA.et al. (2015). Epidemiology of canine distemper and canine parvovirus in domestic dogs in urban and rural areas of the Araucanía region in Chile. Vet. Microbiol.178, 260–264. doi: 10.1016/j.vetmic.2015.05.012
4
Alvarez-EstapeM.FontsereC.Serres-ArmeroA.KudernaL. F. K.DobryninP.GuidaraH.et al. (2022). Insights from the rescue and breeding management of Cuvier’s gazelle (Gazella cuvieri) through whole-genome sequencing. Evolutionary Appl.15, 351–364. doi: 10.1111/eva.13336
5
AmavetP.PolettaG.OdettiL.MarcóM. V. P.SiroskiP. (2023). Detection of the maned wolf, a cryptic and vulnerable species, through environmental DNA studies. J. Nat. Conserv.74, 126439. doi: 10.1016/j.jnc.2023.126439
6
AndresK. J.LodgeD. M.SethiS. A.AndrésJ. (2023). Detecting and analysing intraspecific genetic variation with eDNA: From population genetics to species abundance. Mol. Ecol.32, 4118–4132. doi: 10.1111/mec.v32.15
7
AndrewsK. R.GoodJ. M.MillerM. R.LuikartG.HohenloheP. A. (2016). Harnessing the power of RADseq for ecological and evolutionary genomics. Nat. Rev. Genet.17, 81–92. doi: 10.1038/nrg.2015.28
8
ArantesL. S.FerreiraL. C. L.DrillerM.FilhoF. P. M.MazzoniC.SantosF. R. (2020). Genomic evidence of recent hybridization between sea turtles at Abrolhos Archipelago and its association to low reproductive output. Sci. Rep.10, 12847. doi: 10.1038/s41598-020-69613-8
9
Balagué-DobónL.CáceresA.GonzálezJ. R. (2022). Fully exploiting SNP arrays: a systematic review on the tools to extract underlying genomic structure. Brief Bioinform.23, bbac043. doi: 10.1093/bib/bbac043
10
BassD.ChristisonK. W.StentifordG. D.CookL. S. J.HartikainenH. (2023). Environmental DNA/RNA for pathogen and parasite detection, surveillance, and ecology. Trends Parasitol.39, 285–304. doi: 10.1016/j.pt.2022.12.010
11
BellD. A.RobinsonZ. L.FunkW. C.FitzpatrickS. W.AllendorfF. W.TallmonD. A.et al. (2019). The exciting potential and remaining uncertainties of genetic rescue. Trends Ecol. Evol.34, 1070–1079. doi: 10.1016/j.tree.2019.06.006
12
BengK. C.CorlettR. T. (2020). Applications of environmental DNA (eDNA) in ecology and conservation: opportunities, challenges and prospects. Biodivers Conserv.29, 2089–2121. doi: 10.1007/s10531-020-01980-0
13
BlanchongJ. A.RobinsonS. J.SamuelM. D.FosterJ. T. (2016). Application of genetics and genomics to wildlife epidemiology: Genetics and Wildlife Epidemiology. J. Wildl Manage80, 593–608. doi: 10.1002/jwmg.v80.4
14
CabelloJ. E.DávilaJ. A. (2014). Isolation and characterization of microsatellite loci in Darwin’s fox (Lycalopex fulvipes) and cross-amplification in other canid species. Conserv. Genet. Resour6, 759–761. doi: 10.1007/s12686-014-0208-6
15
CabelloJ.EsperónF.NapolitanoC.HidalgoE.DávilaJ. A.MillánJ. (2013). Molecular identification of a novel gammaherpesvirus in the endangered Darwin’s fox (Lycalopex fulvipes). J. Gen. Virol.94, 2745–2749. doi: 10.1099/vir.0.057851-0
16
CairnsK. M.ShannonL. M.Koler-MatznickJ.BallardJ. W. O.BoykoA. R. (2018). Elucidating biogeographical patterns in Australian native canids using genome wide SNPs. PLoS One13, e0198754. doi: 10.1371/journal.pone.0198754
17
CarvalhoC. S.de OliveiraM. E.Rodriguez-CastroK. G.SaranholiB. H.GalettiP. M.Jr. (2022). Efficiency of eDNA and iDNA in assessing vertebrate diversity and its abundance. Mol. Ecol. Resour.22, 1262–1273. doi: 10.1111/1755-0998.13543
18
ChavezD. E.GronauI.HainsT.DikowR. B.FrandsenP. B.FigueiróH. V.et al. (2022). Comparative genomics uncovers the evolutionary history, demography, and molecular adaptations of South American canids. Proc. Natl. Acad. Sci. U.S.A.119, e2205986119. doi: 10.1073/pnas.2205986119
19
ChemisquyM. A.PrevostiF. J.MartínezP.RaimondiV.StomJ. E. C.Acosta-JamettG.et al. (2019). How many species of grey foxes (Canidae, Carnivora) are there in southern South America? Mastozool Neotrop26, 81–97. doi: 10.31687/saremMN.19.26.1.0.16
20
Chiloé SilvestreO. N. G. (2023). Plan de Recuperación, Conservación y Gestión del Zorro Chilote O de Darwin (Lycalopex fulvipes). Chile: Ministerio del Medio Ambiente. Available online at: https://consultasciudadanas.mma.gob.cl/storage/consultation/YBOBEjvOGhxEd2VNvPpBgKcV71lrPOHRlwHc5eIa.pdf.
21
CleavelandS.MlengeyaT.KaareM.HaydonD.LemboT.LaurensonM. K.et al. (2007). The conservation relevance of epidemiological research into carnivore viral diseases in the serengeti. Conserv. Biol.21, 612–622. doi: 10.1111/j.1523-1739.2007.00701.x
22
CockerillC. A.HasselgrenM.DussexN.DalénL.von SethJ.AngerbjörnA.et al. (2022). Genomic consequences of fragmentation in the endangered Fennoscandian arctic fox (Vulpes lagopus). Genes13 (11), 2124. doi: 10.3390/genes13112124
23
ComizzoliP. (2017). Biobanking and fertility preservation for rare and endangered species. Anim. Reprod.14, 30–33. doi: 10.21451/1984-3143-AR889
24
D’elíaOrtloffSanchezGuinez (2013). A new geographic record of the endangered Darwin’s fox Lycalopex fulvipes (Carnivora: Canidae): filling the distributional gap. Rev. Chil Anat. 86 (4), 485–488. doi: 10.4067/S0716-078X2013000400010
25
Di CataldoS.Hidalgo-HermosoE.SacristánI.CevidanesA.NapolitanoC.HernándezC. V.et al. (2020). Hemoplasmas are endemic and cause asymptomatic infection in the endangered Darwin’s fox (Lycalopex fulvipes). Appl. Environ. Microbiol.86 (12), e00779-20. doi: 10.1128/AEM.00779-20
26
EscobarL. E.QiaoH.CabelloJ.PetersonA. T. (2018). Ecological niche modeling re-examined: A case study with the Darwin’s fox. Ecol. Evol.8, 4757–4770. doi: 10.1002/ece3.2018.8.issue-10
27
EscodaL.HawlitschekO.González-EstebanJ.CastresanaJ. (2022). Methodological challenges in the genomic analysis of an endangered mammal population with low genetic diversity. Sci. Rep.12, 21390. doi: 10.1038/s41598-022-25619-y
28
FariasA. A.SepúlvedaM. A.Silva-RodríguezE. A.EgurenA.GonzálezD.JordánN. I.et al. (2014). A new population of Darwin’s fox (Lycalopex fulvipes) in the Valdivian Coastal Range. Rev. Chil Hist Nat.87, 3. doi: 10.1186/0717-6317-87-3
29
FavariniM. O.SimãoT. L. L.MacedoG. S.GarcezF. S.OliveiraL. R.Cárdenas-AlayzaS.et al. (2022). Complex evolutionary history of the South American fox genus lycalopex (Mammalia, Carnivora, Canidae) inferred from multiple mitochondrial and nuclear markers. Diversity14, 642. doi: 10.3390/d14080642
30
FörsterD. W.BullJ. K.LenzD.AutenriethM.PaijmansJ. L. A.KrausR. H. S.et al. (2018). Targeted resequencing of coding DNA sequences for SNP discovery in nonmodel species. Mol. Ecol. Resour18, 1356–1373. doi: 10.1111/1755-0998.12924
31
GarcezF. S.TchaickaL.LemosF. G.KasperC. B.DalponteJ. C.EizirikE. (2024). Phylogeographic analyses of an endemic Neotropical fox (Lycalopex vetulus) reveal evidence of hybridization with a different canid species (L. gymnocercus). J. Hered. 115 (4), 399–410. doi: 10.1093/jhered/esae012
32
HarrisS. E.Munshi-SouthJ.ObergfellC.O’NeillR. (2013). Signatures of rapid evolution in urban and rural transcriptomes of white-footed mice (Peromyscus leucopus) in the New York metropolitan area. PLoS One8, e74938. doi: 10.1371/journal.pone.0074938
33
HasselgrenM.DussexN.von SethJ.AngerbjörnA.OlsenR.-A.DalénL.et al. (2021). Genomic and fitness consequences of inbreeding in an endangered carnivore. Mol. Ecol.30, 2790–2799. doi: 10.1111/mec.v30.12
34
Hidalgo-HermosoE.CabelloJ.VegaC.Kroeger-GómezH.Moreira-ArceD.NapolitanoC.et al. (2020). An eight-year survey for canine distemper virus indicates lack of exposure in the endangered Darwin’s fox (Lycalopex fulvipes). J. Wildl Dis.56, 482–485. doi: 10.7589/2019-08-195
35
Hidalgo-HermosoE.CabelloJ.VerasayJ.Moreira-ArceD.HidalgoM.AbalosP.et al. (2022). Serosurvey for selected parasitic and bacterial pathogens in Darwin’s fox (Lycalopex fulvipes): Not only dog diseases are a threat. J. Wildl Dis.58, 76–85. doi: 10.7589/JWD-D-21-00024
36
Howard-McCombeJ.JamiesonA.CarmagniniA.RussoI.-R. M.GhazaliM.CampbellR.et al. (2023). Genetic swamping of the critically endangered Scottish wildcat was recent and accelerated by disease. Curr. Biol.33, 4761–9.e5. doi: 10.1016/j.cub.2023.10.026
37
JiménezJ. E. (2007). Ecology of a coastal population of the critically endangered Darwin’s fox (Pseudalopex fulvipes) on Chiloé Island, southern Chile. J. Zool271, 63–77. doi: 10.1111/j.1469-7998.2006.00218.x
38
KimM.BaeJ.OhB.RhimH.YangM.-S.YangS.et al. (2023). Surveillance of wild animals carrying infectious agents based on high-throughput screening platform in the Republic of Korea. BMC Vet. Res.19, 158. doi: 10.1186/s12917-023-03714-0
39
Kleinman-RuizD.Lucena-PerezM.VillanuevaB.FernándezJ.SaveljevA. P.RatkiewiczM.et al. (2022). Purging of deleterious burden in the endangered Iberian lynx. Proc. Natl. Acad. Sci.119, e2110614119. doi: 10.1073/pnas.2110614119
40
Kleinman-RuizD.SorianoL.Casas-MarceM.SzychtaC.SánchezI.FernándezJ.et al. (2019). Genetic evaluation of the Iberian lynx ex situ conservation programme. Heredity123, 647–661. doi: 10.1038/s41437-019-0217-z
41
LiuG.ZhangH.SunG.ZhaoC.ShangS.GaoX.et al. (2017). Characterization of the peripheral blood transcriptome and adaptive evolution of the MHC I and TLR gene families in the wolf (Canis lupus). BMC Genomics18, 584. doi: 10.1186/s12864-017-3983-0
42
MedelR. G.JiménezJ. E.JaksićF. M.YáñezJ.ArmestoJ. J. (1990). Discovery of a continental population of the rare Darwin’s fox, Dusicyon fulvipes (Martin, 1837) in Chile. Biol. Conserv.51, 71–77. doi: 10.1016/0006-3207(90)90033-L
43
MichelA. L.Van HeerdenH.CrossleyB. M.Al DahoukS.PrasseD.RuttenV. (2021). Pathogen detection and disease diagnosis in wildlife: challenges and opportunities: -EN- -FR- Détection des agents pathogènes et diagnostic des maladies dans la faune sauvage: difficultés spécifiques et perspectives -ES- Detección de patógenos y diagnóstico de enfermedades en la fauna silvestre: dificultades y oportunidades. Rev. Sci. Tech40, 105–118. doi: 10.20506/rst.issue.40.1.3205
44
MolinaC.CastilloA.SamaniegoH. (2018). Evaluación del nicho ambiental de Lycalopex fulvipes (zorro de Darwin) y la incidencia del cambio climático sobre su distribución geográfica. Gayana82, 65–78. doi: 10.4067/S0717-65382018000100065
45
Moreira-ArceD.VergaraP. M.BoutinS.CarrascoG.BrionesR.SotoG. E.et al. (2016). Mesocarnivores respond to fine-grain habitat structure in a mosaic landscape comprised by commercial forest plantations in southern Chile. For Ecol. Manage369, 135–143. doi: 10.1016/j.foreco.2016.03.024
46
Morrison WRI. I. I.LohrJ. L.DuchenP.WilchesR.TrujilloD.MairM.et al. (2009). The impact of taxonomic change on conservation: Does it kill, can it save, or is it just irrelevant? Biol. Conserv.142, 3201–3206. doi: 10.1016/j.biocon.2009.07.019
47
NapolitanoC.ClavijoC.Rojas-BonziV.MiñoC. I.González-MayaJ. F.BouN.et al. (2024). Understanding the conservation-genetics gap in Latin America: challenges and opportunities to integrate genetics into conservation practices. Front. Genet.15, 1425531. doi: 10.3389/fgene.2024.1425531
48
OnoratoD. P.CunninghamM. W.LotzM.CriffieldM.ShindleD.JohnsonA.et al. (2024). Multi-generational benefits of genetic rescue. Sci. Rep.14, 17519. doi: 10.1038/s41598-024-67033-6
49
PaijmansJ. L. A.BarlowA.HennebergerK.FickelJ.HofreiterM.FoersterD. W. G. (2020). Ancestral mitogenome capture of the Southeast Asian banded linsang. PloS One15, e0234385. doi: 10.1371/journal.pone.0234385
50
PeriniF. A.RussoC. A. M.SchragoC. G. (2010). The evolution of South American endemic canids: a history of rapid diversification and morphological parallelism. J. Evol. Biol.23, 311–322. doi: 10.1111/j.1420-9101.2009.01901.x
51
PerryG. H.MelstedP.MarioniJ. C.WangY.BainerR.PickrellJ. K.et al. (2012). Comparative RNA sequencing reveals substantial genetic variation in endangered primates. Genome Res.22, 602–610. doi: 10.1101/gr.130468.111
52
PizarroE. J.Julio-KalajžićB.Sallaberry-PincheiraN.MuñozV.González-AcuñaD.CabelloJ.et al. (2023). Species delimitation and intraspecific diversification in recently diverged South American foxes. Mamm Res. 69, 71–87. doi: 10.1007/s13364-023-00717-y
53
PlasilM.FutasJ.JelinekA.BurgerP. A.HorinP. (2022). Comparative genomics of the major histocompatibility complex (MHC) of felids. Front. Genet.13, 829891. doi: 10.3389/fgene.2022.829891
54
QuinnC. B.Preckler-QuisquaterS.BuchalskiM. R.SacksB. N. (2024). Whole genomes inform genetic rescue strategy for Montane red foxes in north America. Mol. Biol. Evol.41, msae193. doi: 10.1093/molbev/msae193
55
RakotoariveloA. R.RambudaT.TaronU. H.StalderG.O’DonoghueP.RobovskýJ.et al. (2024). Complex patterns of gene flow and convergence in the evolutionary history of the spiral-horned antelopes (Tragelaphini). Mol. Phylogenet Evol.198, 108131. doi: 10.1016/j.ympev.2024.108131
56
RobinsonJ.KyriazisC. C.YuanS. C.LohmuellerK. E. (2023). Deleterious variation in natural populations and implications for conservation genetics. Annu. Rev. Anim. Biosci.11, 93–114. doi: 10.1146/annurev-animal-080522-093311
57
RobinsonJ. A.RäikkönenJ.VucetichL. M.VucetichJ. A.PetersonR. O.LohmuellerK. E.et al. (2019). Genomic signatures of extensive inbreeding in Isle Royale wolves, a population on the threshold of extinction. Sci. Adv.5, eaau0757. doi: 10.1126/sciadv.aau0757
58
SeeberP. A.McEwenG. K.LöberU.FörsterD. W.EastM. L.MelzheimerJ.et al. (2019). Terrestrial mammal surveillance using hybridization capture of environmental DNA from African waterholes. Mol. Ecol. Resour19, 1486–1496. doi: 10.1111/1755-0998.13069
59
Sillero-Zubiri E. byC.HoffmannM.MacdonaldD. W. (2004). Canids: Foxes, Wolves, Jackals and Dogs. Sillero-ZubiriCHoffmannMMacdonaldDW. (IUCN, Gland, Switzerland and Cambridge, UK).
60
Silva-RodríguezE.FariasA.Moreira-ArceD.CabelloJHidalgo-HermosoE.LucheriniM.et al. (2016). Lycalopex fulvipes (errata version published in 2016). The IUCN Red List of Threatened Species 2016: e.T41586A107263066. doi: 10.2305/IUCN.UK.2016-1.RLTS.T41586A85370871.en
61
Silva-RodríguezE. A.OvandoE.GonzálezD.ZambranoB.SepúlvedaM. A.SvenssonG. L.et al. (2018). Large-scale assessment of the presence of Darwin’s fox across its newly discovered range. Mamm Biol.92, 45–53. doi: 10.1016/j.mambio.2018.04.003
62
SommerS. (2005). The importance of immune gene variability (MHC) in evolutionary ecology and conservation. Front. Zool2, 16. doi: 10.1186/1742-9994-2-16
63
StangeM.BarrettR. D. H.HendryA. P. (2021). The importance of genomic variation for biodiversity, ecosystems and people. Nat. Rev. Genet.22, 89–105. doi: 10.1038/s41576-020-00288-7
64
TchaickaL.de FreitasT. R. O.BagerA.VidalS. L.LucheriniM.IriarteA.et al. (2016). Molecular assessment of the phylogeny and biogeography of a recently diversified endemic group of South American canids (Mammalia: Carnivora: Canidae). Genet. Mol. Biol.39, 442–451. doi: 10.1590/1678-4685-GMB-2015-0189
65
Tomasco IH.GiorelloF. M.BoullosaN.FeijooM.LanzoneC.LessaE. P. (2022). The contribution of incomplete lineage sorting and introgression to the evolutionary history of the fast-evolving genus Ctenomys (Rodentia, Ctenomyidae). Mol. Phylogenet Evol.176, 107593. doi: 10.1016/j.ympev.2022.107593
66
Van BormS.BelákS.FreimanisG.FusaroA.GranbergF.HöperD.et al. (2015). “Next-generation sequencing in veterinary medicine: how can the massive amount of information arising from high-throughput technologies improve diagnosis, control, and management of infectious diseases?,” in Veterinary Infection Biology: Molecular Diagnostics and High-Throughput Strategies, vol. 1247 . Eds. CunhaM. V.InácioJ. (Springer New York, New York, NY), 415–436.
67
VilàC.LeonardJ. A.IriarteA.O’BrienS. J.JohnsonW. E.WayneR. K. (2004). Detecting the vanishing populations of the highly endangered Darwin’s fox, Pseudalopex fulvipes. Anim. Conserv. Forum7, 147–153. doi: 10.1017/S1367943004001271
68
von HoldtK. (2022). “Conclusion: towards a southern sociology,” in Critical Engagement with Public Sociology, vol. p . Eds. BezuidenhoutA.SonwabileM.von HoldtK. (Bristol University Press), 265–279.
69
vonHoldtB. M.PollingerJ. P.EarlD. A.KnowlesJ. C.BoykoA. R.ParkerH.et al. (2011). A genome-wide perspective on the evolutionary history of enigmatic wolf-like canids. Genome Res.21, 1294–1305. doi: 10.1101/gr.116301.110
70
vonHoldtB. M.PollingerJ. P.EarlD. A.ParkerH. G.OstranderE. A.WayneR. K. (2013). Identification of recent hybridization between gray wolves and domesticated dogs by SNP genotyping. Mamm Genome24, 80–88. doi: 10.1007/s00335-012-9432-0
71
WangD.CoscoyL.ZylberbergM.AvilaP. C.BousheyH. A.GanemD.et al. (2002). Microarray-based detection and genotyping of viral pathogens. Proc. Natl. Acad. Sci. U.S.A.99, 15687–15692. doi: 10.1073/pnas.242579699
72
WiselyS. M.RyderO. A.SantymireR. M.EngelhardtJ. F.NovakB. J. (2015). A road map for 21st century genetic restoration: gene pool enrichment of the black-footed ferret. J. Heredity106, 581–592. doi: 10.1093/jhered/esv041
73
WrightB. R.FarquharsonK. A.McLennanE. A.BelovK.HoggC. J.GrueberC. E. (2020). A demonstration of conservation genomics for threatened species management. Mol. Ecol. Resour20, 1526–1541. doi: 10.1111/1755-0998.13211
74
YahnkeC. J.JohnsonW. E.GeffenE.SmithD.HertelF.RoyM. S.et al. (1996). Darwin’s fox: A distinct endangered species in a vanishing habitat. Conserv. Biol.10, 366–375. doi: 10.1046/j.1523-1739.1996.10020366.x
75
YuhkiN.BeckT.StephensR.NeelamB.O’BrienS. J. (2007). Comparative genomic structure of human, dog, and cat MHC: HLA, DLA, and FLA. J. Hered98, 390–399. doi: 10.1093/jhered/esm056
Summary
Keywords
conservation genomics, high-throughput sequencing, Lycalopex fulvipes, genetic diversity, Chile
Citation
Valenzuela-Turner C, Grau JH, Fickel J and Förster DW (2025) Preserving Darwin’s fox: genomic tools for the conservation of South America’s most endangered canid. Front. Conserv. Sci. 6:1512531. doi: 10.3389/fcosc.2025.1512531
Received
16 October 2024
Accepted
06 January 2025
Published
23 January 2025
Volume
6 - 2025
Edited by
Juan Pablo Jaramillo-Correa, National Autonomous University of Mexico, Mexico
Reviewed by
Jaime Gasca-Pineda, Posdoctoral Fellow, Mexico
Gustavo P. Lorenzana, Universidad de la Sierra, Mexico
Updates

Check for updates
Copyright
© 2025 Valenzuela-Turner, Grau, Fickel and Förster.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Cristóbal Valenzuela-Turner, valenzuela@izw-berlin.de
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.