Epigenome-Wide Study Identifies Epigenetic Outliers in Normal Mucosa of Patients with Colorectal Cancer.

Nongenetic predisposition to colorectal cancer continues to be difficult to measure precisely, hampering efforts in targeted prevention and screening. Epigenetic changes in the normal mucosa of patients with colorectal cancer can serve as a tool in predicting colorectal cancer outcomes. We identified epigenetic changes affecting the normal mucosa of patients with colorectal cancer. DNA methylation profiling on normal colon mucosa from 77 patients with colorectal cancer and 68 controls identified a distinct subgroup of normally-appearing mucosa with markedly disrupted DNA methylation at a large number of CpGs, termed as "Outlier Methylation Phenotype" (OMP) and are present in 15 of 77 patients with cancer versus 0 of 68 controls (P < 0.001). Similar findings were also seen in publicly available datasets. Comparison of normal colon mucosa transcription profiles of patients with OMP cancer with those of patients with non-OMP cancer indicates genes whose promoters are hypermethylated in the OMP patients are also transcriptionally downregulated, and that many of the genes most affected are involved in interactions between epithelial cells, the mucus layer, and the microbiome. Analysis of 16S rRNA profiles suggests that normal colon mucosa of OMPs are enriched in bacterial genera associated with colorectal cancer risk, advanced tumor stage, chronic intestinal inflammation, malignant transformation, nosocomial infections, and KRAS mutations. In conclusion, our study identifies an epigenetically distinct OMP group in the normal mucosa of patients with colorectal cancer that is characterized by a disrupted methylome, altered gene expression, and microbial dysbiosis. Prospective studies are needed to determine whether OMP could serve as a biomarker for an elevated epigenetic risk for colorectal cancer development. PREVENTION RELEVANCE: Our study identifies an epigenetically distinct OMP group in the normal mucosa of patients with colorectal cancer that is characterized by a disrupted methylome, altered gene expression, and microbial dysbiosis. Identification of OMPs in healthy controls and patients with colorectal cancer will lead to prevention and better prognosis, respectively.

DNA metabarcoding reveals that coyotes in New York City consume wide variety of native prey species and human food.

Carnivores are currently colonizing cities where they were previously absent. These urban environments are novel ecosystems characterized by habitat degradation and fragmentation, availability of human food, and different prey assemblages than surrounding areas. Coyotes (Canis latrans) established a breeding population in New York City (NYC) over the last few decades, but their ecology within NYC is poorly understood. In this study, we used non-invasive scat sampling and DNA metabarcoding to profile vertebrate, invertebrate, and plant dietary items with the goal to compare the diets of urban coyotes to those inhabiting non-urban areas. We found that both urban and non-urban coyotes consumed a variety of plants and animals as well as human food. Raccoons (Procyon lotor) were an important food item for coyotes within and outside NYC. In contrast, white-tailed deer (Odocoileus virginianus) were mainly eaten by coyotes inhabiting non-urban areas. Domestic chicken (Gallus gallus) was the human food item found in most scats from both urban and non-urban coyotes. Domestic cats (Felis catus) were consumed by urban coyotes but were detected in only a small proportion of the scats (<5%), which differs markedly from high rates of cat depredation in some other cities. In addition, we compared our genetic metabarcoding analysis to a morphological analysis of the same scat samples. We found that the detection similarity between the two methods was low and it varied depending on the type of diet item.

Cascading effects of habitat loss on ectoparasite-associated bacterial microbiomes.

Suitable habitat fragment size, isolation, and distance from a source are important variables influencing community composition of plants and animals, but the role of these environmental factors in determining composition and variation of host-associated microbial communities is poorly known. In parasite-associated microbial communities, it is hypothesized that evolution and ecology of an arthropod parasite will influence its microbiome more than broader environmental factors, but this hypothesis has not been extensively tested. To examine the influence of the broader environment on the parasite microbiome, we applied high-throughput sequencing of the V4 region of 16S rRNA to characterize the microbiome of 222 obligate ectoparasitic bat flies (Streblidae and Nycteribiidae) collected from 155 bats (representing six species) from ten habitat fragments in the Atlantic Forest of Brazil. Parasite species identity is the strongest driver of microbiome composition. To a lesser extent, reduction in habitat fragment area, but not isolation, is associated with an increase in connectance and betweenness centrality of bacterial association networks driven by changes in the diversity of the parasite community. Controlling for the parasite community, bacterial network topology covaries with habitat patch area and exhibits parasite-species specific responses to environmental change. Taken together, habitat loss may have cascading consequences for communities of interacting macro- and microorgansims.

You are more than what you eat: potentially adaptive enrichment of microbiome functions across bat dietary niches.

BACKGROUND: Animals evolved in a microbial world, and their gut microbial symbionts have played a role in their ecological diversification. While many recent studies report patterns of phylosymbiosis between hosts and their gut bacteria, fewer studies examine the potentially adaptive functional contributions of these microbes to the dietary habits of their hosts. In this study, we examined predicted metabolic pathways in the gut bacteria of more than 500 individual bats belonging to 60 species and compare the enrichment of these functions across hosts with distinct dietary ecologies. RESULTS: We found that predicted microbiome functions were differentially enriched across hosts with different diets. Using a machine-learning approach, we also found that inferred microbiome functions could be used to predict specialized host diets with reasonable accuracy. We detected a relationship between both host phylogeny and diet with respect to microbiome functional repertoires. Because many predicted functions could potentially fill nutritional gaps for bats with specialized diets, we considered pathways discriminating dietary niches as traits of the host and fit them to comparative phylogenetic models of evolution. Our results suggest that some, but not all, predicted microbiome functions may evolve toward adaptive optima and thus be visible to the forces of natural selection operating on hosts over evolutionary time. CONCLUSIONS: Our results suggest that bats with specialized diets may partially rely on their gut microbes to fulfill or augment critical nutritional pathways, including essential amino acid synthesis, fatty acid biosynthesis, and the generation of cofactors and vitamins essential for proper nutrition. Our work adds to a growing body of literature suggesting that animal microbiomes are structured by a combination of ecological and evolutionary processes and sets the stage for future metagenomic and metabolic characterization of the bat microbiome to explore links between bacterial metabolism and host nutrition.

Molecular diet analysis of neotropical bats based on fecal DNA metabarcoding.

Bat communities in the Neotropics are some of the most speciose assemblages of mammals on Earth, with regions supporting more than 100 sympatric species with diverse feeding ecologies. Because bats are small, nocturnal, and volant, it is difficult to directly observe their feeding habits, which has resulted in their classification into broadly defined dietary guilds (e.g., insectivores, carnivores, and frugivores). Apart from these broad guilds, we lack detailed dietary information for many species and therefore have only a limited understanding of interaction networks linking bats and their diet items. In this study, we used DNA metabarcoding of plants, arthropods, and vertebrates to investigate the diets of 25 bat species from the tropical dry forests of Lamanai, Belize. Our results report some of the first detection of diet items for the focal bat taxa, adding rich and novel natural history information to the field of bat ecology. This study represents a comprehensive first effort to apply DNA metabarcoding to bat diets at Lamanai and provides a useful methodological framework for future studies testing hypotheses about coexistence and niche differentiation in the context of modern high-throughput molecular data.

Population genetic structure and habitat connectivity for jaguar (Panthera onca) conservation in Central Belize.

BACKGROUND: Connectivity among jaguar (Panthera onca) populations will ensure natural gene flow and the long-term survival of the species throughout its range. Jaguar conservation efforts have focused primarily on connecting suitable habitat in a broad-scale. Accelerated habitat reduction, human-wildlife conflict, limited funding, and the complexity of jaguar behaviour have proven challenging to maintain connectivity between populations effectively. Here, we used non-invasive genetic sampling and individual-based conservation genetic analyses to assess genetic diversity and levels of genetic connectivity between individuals in the Cockscomb Basin Wildlife Sanctuary and the Maya Forest Corridor. We used expert knowledge and scientific literature to develop models of landscape permeability based on circuit theory with fine-scale landscape features as ecosystem types, distance to human settlements and roads to predict the most probable jaguar movement across central Belize. RESULTS: We used 12 highly polymorphic microsatellite loci to identify 50 individual jaguars. We detected high levels of genetic diversity across loci (HE = 0.61, HO = 0.55, and NA = 9.33). Using Bayesian clustering and multivariate models to assess gene flow and genetic structure, we identified one single group of jaguars (K = 1). We identified critical areas for jaguar movement that fall outside the boundaries of current protected areas in central Belize. We detected two main areas of high landscape permeability in a stretch of approximately 18 km between Sittee River Forest Reserve and Manatee Forest Reserve that may increase functional connectivity and facilitate jaguar dispersal from and to Cockscomb Basin Wildlife Sanctuary. Our analysis provides important insights on fine-scale genetic and landscape connectivity of jaguars in central Belize, an area of conservation concern. CONCLUSIONS: The results of our study demonstrate high levels of relatively recent gene flow for jaguars between two study sites in central Belize. Our landscape analysis detected corridors of expected jaguar movement between the Cockscomb Basin Wildlife Sanctuary and the Maya Forest Corridor. We highlight the importance of maintaining already established corridors and consolidating new areas that further promote jaguar movement across suitable habitat beyond the boundaries of currently protected areas. Continued conservation efforts within identified corridors will further maintain and increase genetic connectivity in central Belize.

Gut microbiota and their putative metabolic functions in fragmented Bengal tiger population of Nepal.

Bengal tigers (Panthera tigris tigris) serve a pivotal role as an apex predator in forest ecosystems. To increase our knowledge on factors impacting the viability and health of this endangered species, we studied the gut microbiota in 32 individual Bengal tigers from three geographically separated areas (Chitwan National Park (CNP), Bardia National Park (BNP) and Suklaphanta Wildlife Reserve (SWR)) in Nepal, using noninvasive genetic sampling methods. Gut microbiota influence the immune system, impact various physiological functions, and modulates metabolic reactions, that ultimately impact the host health, behavior and development. Across the tiger populations in Nepal, we found significant differences in the composition of microbial communities based on their geographic locations. Specifically, we detected significant differences between CNP and the other two protected areas (CNP vs BNP: pseudo t = 1.944, P = 0.006; CNP vs SWR: pseudo t = 1.9942, P = 0.0071), but no differences between BNP and SWR. This mirrors what has been found for tiger gene flow in the same populations, suggesting gut microbiota composition and host gene flow may be linked. Furthermore, predictive metagenome functional content analysis (PICRUSt) revealed a higher functional enrichment and diversity for significant gut microbiota in the Chitwan tiger population and the lowest enrichment and diversity in Suklaphanta. The CNP tiger population contained higher proportions of microbiota that are associated with predicted functions relevant for metabolism of amino acid, lipid, xenobiotics biodegradation, terpenoides and polyketides than the SWR population. We conclude the tiger population structure, gut microbiota profile and associated functional metabolic categories are correlated, with geographically most separated CNP and SWR tiger population having the most distinct and different host genotype and microbiota profiles. Our work dramatically expands the understanding of tiger microbiota in wild populations and provides a valuable case study on how to investigate genetic diversity at different hierarchical levels, including hosts as well as their microbial communities.

Does faecal matter reflect location? An initial assessment of isotopic variability between consumed prey remains and faecal matter for wild jaguars.

Faecal isotopic analysis may complement other non-invasive wildlife survey tools for monitoring landscape use by carnivores, such as motion-detecting cameras and non-invasive genetic sampling. We analysed carbon, nitrogen, and strontium isotopes in faecal matter produced by jaguars (Panthera onca) as well as bones from consumed prey at the Mountain Pine Ridge Forest Reserve (MPR) in Belize, Central America. The MPR is ideally suited for a spatial isotope study as vegetation and geology both vary considerably. The isotopic composition of faecal matter should reflect the habitat and geology where consumed prey lived. We used bone from consumed prey recovered from jaguar scats as a proxy for diet. Faecal matter and bone showed comparable spatial isotopic trends, suggesting that the isotopic composition of jaguar faeces can be used to detect foraging in different habitats (pine forest versus broadleaf forest) or on different geologies (Mesozoic carbonates; Palaeozoic granite, contact metamorphics, and metasediments). This result is reassuring as bones are not always present in carnivore scats. Studying landscape use by cryptic and wide-ranging carnivore species like jaguars remains challenging. Isotopic analysis of faecal matter complements the existing array of non-invasive spatial monitoring tools.

MGS-Fast: Metagenomic shotgun data fast annotation using microbial gene catalogs.

BACKGROUND: Current methods used for annotating metagenomics shotgun sequencing (MGS) data rely on a computationally intensive and low-stringency approach of mapping each read to a generic database of proteins or reference microbial genomes. RESULTS: We developed MGS-Fast, an analysis approach for shotgun whole-genome metagenomic data utilizing Bowtie2 DNA-DNA alignment of reads that is an alternative to using the integrated catalog of reference genes database of well-annotated genes compiled from human microbiome data. This method is rapid and provides high-stringency matches (>90% DNA sequence identity) of the metagenomics reads to genes with annotated functions. We demonstrate the use of this method with data from a study of liver disease and synthetic reads, and Human Microbiome Project shotgun data, to detect differentially abundant Kyoto Encyclopedia of Genes and Genomes gene functions in these experiments. This rapid annotation method is freely available as a Galaxy workflow within a Docker image. CONCLUSIONS: MGS-Fast can confidently transfer functional annotations from gene databases to metagenomic reads, with speed and accuracy.

miCloud: A Plug-n-Play, Extensible, On-Premises Bioinformatics Cloud for Seamless Execution of Complex Next-Generation Sequencing Data Analysis Pipelines.

The availability of low-cost small-factor sequencers, such as the Illumina MiSeq, MiniSeq, or iSeq, have paved the way for democratizing genomics sequencing, providing researchers in minority universities with access to the technology that was previously only affordable by institutions with large core facilities. However, these instruments are not bundled with software for performing bioinformatics data analysis, and the data analysis can be the main bottleneck for independent laboratories or even small clinical facilities that consider adopting genomic sequencing for medical applications. To address this issue, we have developed miCloud, a bioinformatics platform that enables genomic data analysis through a fully featured data analysis cloud, which seamlessly integrates with genome sequencers over the local network. The miCloud can be easily deployed without any prior bioinformatics expertise on any computing environment, from a laboratory computer workstation to a university computer cluster. Our platform not only provides access to a set of preconfigured RNA-Seq and CHIP-Seq bioinformatics pipelines, but also enables users to develop or install new preconfigured tools from the large selection available on open-source online Docker container repositories. The miCloud built-in analysis pipelines are also integrated with the Visual Omics Explorer framework (Kim et al., 2016), which provides rich interactive visualizations and publication-ready graphics from the next-generation sequencing data. Ultimately, the miCloud demonstrates a bioinformatics approach that can be adopted in the field for standardizing genomic data analysis, similarly to the way molecular biology sample preparation kits have standardized laboratory operations.

Species, sex and geo-location identification of seized tiger (Panthera tigris tigris) parts in Nepal-A molecular forensic approach.

Tiger (Panthera tigris) populations are in danger across their entire range due to habitat loss, poaching and the demand for tiger parts. The Bengal tiger (Panthera tigris tigris) is an endangered apex predator with a population size estimated to be less than 200 in Nepal. In spite of strict wildlife protection laws, illegal trade of tiger parts is increasing; and Nepal has become one of the major sources and transit routes for poached wildlife parts. Identification of wildlife parts is often challenging for law enforcement officials due to inadequate training and lack of available tools. Here, we describe a molecular forensic approach to gain insight into illegally trafficked tiger parts seized across Nepal. We created Nepal's first comprehensive reference genetic database of wild tigers through the Nepal Tiger Genome Project (2011-2013). This database has nuclear DNA microsatellite genotype and sex profiles, including geo-spatial information, of over 60% (n = 120) of the wild tigers of Nepal. We analyzed 15 putative cases of confiscated poached tiger parts and all were confirmed to be of tiger. Ten samples were identified as male and five were female. We determined probable geo-source location for 9 of the 14 samples with 6-8 nuclear DNA microsatellite loci using inferences from four different statistical assignment methods. Six samples were assigned to Bardia National Park and one of these was an exact match to a female tiger previously profiled in our fecal DNA reference database. Two tiger samples were assigned to Shuklaphanta Wildlife Reserve and one to Chitwan National Park. We are unable to definitively assign five tiger samples which could be offspring dispersers or might have come from tiger population outside of Nepal. Our study revealed that the western region, particularly Bardia National Park, is a poaching hotspot for illegal tiger trade in Nepal. We present feasibility of using molecular forensic based evidence to incriminate criminals in a court of law in the fight against wildlife crime.

Ecological correlates of the spatial co-occurrence of sympatric mammalian carnivores worldwide.

The composition of local mammalian carnivore communities has far-reaching effects on terrestrial ecosystems worldwide. To better understand how carnivore communities are structured, we analysed camera trap data for 108 087 trap days across 12 countries spanning five continents. We estimate local probabilities of co-occurrence among 768 species pairs from the order Carnivora and evaluate how shared ecological traits correlate with probabilities of co-occurrence. Within individual study areas, species pairs co-occurred more frequently than expected at random. Co-occurrence probabilities were greatest for species pairs that shared ecological traits including similar body size, temporal activity pattern and diet. However, co-occurrence decreased as compared to other species pairs when the pair included a large-bodied carnivore. Our results suggest that a combination of shared traits and top-down regulation by large carnivores shape local carnivore communities globally.

Comparing Microbiome Sampling Methods in a Wild Mammal: Fecal and Intestinal Samples Record Different Signals of Host Ecology, Evolution.

The gut microbiome is a community of host-associated symbiotic microbes that fulfills multiple key roles in host metabolism, immune function, and tissue development. Given the ability of the microbiome to impact host fitness, there is increasing interest in studying the microbiome of wild animals to better understand these communities in the context of host ecology and evolution. Human microbiome research protocols are well established, but wildlife microbiome research is still a developing field. Currently, there is no standardized set of best practices guiding the collection of microbiome samples from wildlife. Gut microflora are typically sampled either by fecal collection, rectal swabbing, or by destructively sampling the intestinal contents of the host animal. Studies rarely include more than one sampling technique and no comparison of these methods currently exists for a wild mammal. Although some studies have hypothesized that the fecal microbiome is a nested subset of the intestinal microbiome, this hypothesis has not been formally tested. To address these issues, we examined guano (feces) and distal intestinal mucosa from 19 species of free-ranging bats from Lamanai, Belize, using 16S rRNA amplicon sequencing to compare microbial communities across sample types. We found that the diversity and composition of intestine and guano samples differed substantially. In addition, we conclude that signatures of host evolution are retained by studying gut microbiomes based on mucosal tissue samples, but not fecal samples. Conversely, fecal samples retained more signal of host diet than intestinal samples. These results suggest that fecal and intestinal sampling methods are not interchangeable, and that these two microbiotas record different information about the host from which they are isolated.

Assessment of genetic diversity, population structure, and gene flow of tigers (Panthera tigris tigris) across Nepal's Terai Arc Landscape.

With fewer than 200 tigers (Panthera tigris tigris) left in Nepal, that are generally confined to five protected areas across the Terai Arc Landscape, genetic studies are needed to provide crucial information on diversity and connectivity for devising an effective country-wide tiger conservation strategy. As part of the Nepal Tiger Genome Project, we studied landscape change, genetic variation, population structure, and gene flow of tigers across the Terai Arc Landscape by conducting Nepal's first comprehensive and systematic scat-based, non-invasive genetic survey. Of the 770 scat samples collected opportunistically from five protected areas and six presumed corridors, 412 were tiger (57%). Out of ten microsatellite loci, we retain eight markers that were used in identifying 78 individual tigers. We used this dataset to examine population structure, genetic variation, contemporary gene flow, and potential population bottlenecks of tigers in Nepal. We detected three genetic clusters consistent with three demographic sub-populations and found moderate levels of genetic variation (He = 0.61, AR = 3.51) and genetic differentiation (FST = 0.14) across the landscape. We detected 3-7 migrants, confirming the potential for dispersal-mediated gene flow across the landscape. We found evidence of a bottleneck signature likely caused by large-scale land-use change documented in the last two centuries in the Terai forest. Securing tiger habitat including functional forest corridors is essential to enhance gene flow across the landscape and ensure long-term tiger survival. This requires cooperation among multiple stakeholders and careful conservation planning to prevent detrimental effects of anthropogenic activities on tigers.

Range-Wide Snow Leopard Phylogeography Supports Three Subspecies.

The snow leopard, Panthera uncia, is an elusive high-altitude specialist that inhabits vast, inaccessible habitat across Asia. We conducted the first range-wide genetic assessment of snow leopards based on noninvasive scat surveys. Thirty-three microsatellites were genotyped and a total of 683 bp of mitochondrial DNA sequenced in 70 individuals. Snow leopards exhibited low genetic diversity at microsatellites (AN = 5.8, HO = 0.433, HE = 0.568), virtually no mtDNA variation, and underwent a bottleneck in the Holocene (∼8000 years ago) coinciding with increased temperatures, precipitation, and upward treeline shift in the Tibetan Plateau. Multiple analyses supported 3 primary genetic clusters: (1) Northern (the Altai region), (2) Central (core Himalaya and Tibetan Plateau), and (3) Western (Tian Shan, Pamir, trans-Himalaya regions). Accordingly, we recognize 3 subspecies, Panthera uncia irbis (Northern group), Panthera uncia uncia (Western group), and Panthera uncia uncioides (Central group) based upon genetic distinctness, low levels of admixture, unambiguous population assignment, and geographic separation. The patterns of variation were consistent with desert-basin "barrier effects" of the Gobi isolating the northern subspecies (Mongolia), and the trans-Himalaya dividing the central (Qinghai, Tibet, Bhutan, and Nepal) and western subspecies (India, Pakistan, Tajikistan, and Kyrgyzstan). Hierarchical Bayesian clustering analysis revealed additional subdivision into a minimum of 6 proposed management units: western Mongolia, southern Mongolia, Tian Shan, Pamir-Himalaya, Tibet-Himalaya, and Qinghai, with spatial autocorrelation suggesting potential connectivity by dispersing individuals up to ∼400 km. We provide a foundation for global conservation of snow leopard subspecies, and set the stage for in-depth landscape genetics and genomic studies.

Genetic Diversity and Population Structure of Mesoamerican Jaguars (Panthera onca): Implications for Conservation and Management.

Mesoamerican jaguars (Panthera onca) have been extirpated from over 77% of their historic range, inhabiting fragmented landscapes at potentially reduced population sizes. Maintaining and restoring genetic diversity and connectivity across human-altered landscapes has become a major conservation priority; nonetheless large-scale genetic monitoring of natural populations is rare. This is the first regional conservation genetic study of jaguars to primarily use fecal samples collected in the wild across five Mesoamerican countries: Belize, Costa Rica, Guatemala, Honduras, and Mexico. We genotyped 445 jaguar fecal samples and examined patterns of genetic diversity and connectivity among 115 individual jaguars using data from 12 microsatellite loci. Overall, moderate levels of genetic variation were detected (NA = 4.50 ± 1.05, AR = 3.43 ± 0.22, HE = 0.59 ± 0.04), with Mexico having the lowest genetic diversity, followed by Honduras, Guatemala, Belize, and Costa Rica. Population-based gene flow measures (FST = 0.09 to 0.15, Dest = 0.09 to 0.21), principal component analysis, and Bayesian clustering applied in a hierarchical framework revealed significant genetic structure in Mesoamerican jaguars, roughly grouping individuals into four genetic clusters with varying levels of admixture. Gene flow was highest among Selva Maya jaguars (northern Guatemala and central Belize), whereas genetic differentiation among all other sampling sites was moderate. Genetic subdivision was most pronounced between Selva Maya and Honduran jaguars, suggesting limited jaguar movement between these close geographic regions and ultimately refuting the hypothesis of contemporary panmixia. To maintain a critical linkage for jaguars dispersing through the Mesoamerican landscape and ensure long-term viability of this near threatened species, we recommend continued management and maintenance of jaguar corridors. The baseline genetic data provided by this study underscores the importance of understanding levels of genetic diversity and connectivity to making informed management and conservation decisions with the goal to maintain functional connectivity across the region.

A Comparative Analysis of Genetic Diversity and Structure in Jaguars (Panthera onca), Pumas (Puma concolor), and Ocelots (Leopardus pardalis) in Fragmented Landscapes of a Critical Mesoamerican Linkage Zone.

With increasing anthropogenic impact and landscape change, terrestrial carnivore populations are becoming more fragmented. Thus, it is crucial to genetically monitor wild carnivores and quantify changes in genetic diversity and gene flow in response to these threats. This study combined the use of scat detector dogs and molecular scatology to conduct the first genetic study on wild populations of multiple Neotropical felids coexisting across a fragmented landscape in Belize, Central America. We analyzed data from 14 polymorphic microsatellite loci in 1053 scat samples collected from wild jaguars (Panthera onca), pumas (Puma concolor), and ocelots (Leopardus pardalis). We assessed levels of genetic diversity, defined potential genetic clusters, and examined gene flow for the three target species on a countrywide scale using a combination of individual- and population-based analyses. Wild felids in Belize showed moderate levels of genetic variation, with jaguars having the lowest diversity estimates (HE = 0.57 ± 0.02; AR = 3.36 ± 0.09), followed by pumas (HE = 0.57 ± 0.08; AR = 4.20 ± 0.16), and ocelots (HE = 0.63 ± 0.03; AR = 4.16 ± 0.08). We observed low to moderate levels of genetic differentiation for all three target species, with jaguars showing the lowest degree of genetic subdivision across the country, followed by ocelots and pumas. Although levels of genetic diversity and gene flow were still fairly high, we detected evidence of fine-scale genetic subdivision, indicating that levels of genetic connectivity for wild felids in Belize are likely to decrease if habitat loss and fragmentation continue at the current rate. Our study demonstrates the value of understanding fine-scale patterns of gene flow in multiple co-occurring felid species of conservation concern, which is vital for wildlife movement corridor planning and prioritizing future conservation and management efforts within human-impacted landscapes.

Noninvasive individual and species identification of jaguars (Panthera onca), pumas (Puma concolor) and ocelots (Leopardus pardalis) in Belize, Central America using cross-species microsatellites and faecal DNA.

There is a great need to develop efficient, noninvasive genetic sampling methods to study wild populations of multiple, co-occurring, threatened felids. This is especially important for molecular scatology studies occurring in challenging tropical environments where DNA degrades quickly and the quality of faecal samples varies greatly. We optimized 14 polymorphic microsatellite loci for jaguars (Panthera onca), pumas (Puma concolor) and ocelots (Leopardus pardalis) and assessed their utility for cross-species amplification. Additionally, we tested their reliability for species and individual identification using DNA from faeces of wild felids detected by a scat detector dog across Belize in Central America. All microsatellite loci were successfully amplified in the three target species, were polymorphic with average expected heterozygosities of HE = 0.60 ± 0.18 (SD) for jaguars, HE = 0.65 ± 0.21 (SD) for pumas and HE = 0.70 ± 0.13 (SD) for ocelots and had an overall PCR amplification success of 61%. We used this nuclear DNA primer set to successfully identify species and individuals from 49% of 1053 field-collected scat samples. This set of optimized microsatellite multiplexes represents a powerful tool for future efforts to conduct noninvasive studies on multiple, wild Neotropical felids.