Draft assembly and annotation of the Cuban crocodile (Crocodylus rhombifer) genome.

OBJECTIVES: The new data provide an important genomic resource for the Critically Endangered Cuban crocodile (Crocodylus rhombifer). Cuban crocodiles are restricted to the Zapata Swamp in southern Matanzas Province, Cuba, and readily hybridize with the widespread American crocodile (Crocodylus acutus) in areas of sympatry. The reported de novo assembly will contribute to studies of crocodylian evolutionary history and provide a resource for informing Cuban crocodile conservation. DATA DESCRIPTION: The final 2.2 Gb draft genome for C. rhombifer consists of 41,387 scaffolds (contigs: N50 = 104.67 Kb; scaffold: N50-518.55 Kb). Benchmarking Universal Single-Copy Orthologs (BUSCO) identified 92.3% of the 3,354 genes in the vertebrata_odb10 database. Approximately 42% of the genome (960Mbp) comprises repeat elements. We predicted 30,138 unique protein-coding sequences (17,737 unique genes) in the genome assembly. Functional annotation found the top Gene Ontology annotations for Biological Processes, Molecular Function, and Cellular Component were regulation, protein, and intracellular, respectively. This assembly will support future macroevolutionary, conservation, and molecular studies of the Cuban crocodile.

Environmental DNA-based biomonitoring of Cuban Crocodylus and their accompanying vertebrate fauna from Zapata Swamp, Cuba.

Crocodylians globally face considerable challenges, including population decline and extensive habitat modification. Close monitoring of crocodylian populations and their habitats is imperative for the timely detection of population trends, especially in response to management interventions. Here we use eDNA metabarcoding to identify the Critically Endangered Crocodylus rhombifer and the Vulnerable C. acutus, as well as vertebrate community diversity, in Cuba's Zapata Swamp. We tested four different primer sets, including those used previously in Crocodylus population genetic and phylogenetic research, for their efficiency at detecting crocodylian eDNA. We detected C. rhombifer eDNA in 11 out of 15 sampled locations within its historical geographic distribution. We found that data analyses using the VertCOI primers and the mBRAVE bioinformatics pipeline were the most effective molecular marker and pipeline combination for identifying this species from environmental samples. We also identified 55 vertebrate species in environmental samples across the four bioinformatics pipelines- ~ 85% known to be present in the Zapata ecosystem. Among them were eight species previously undetected in the area and eight alien species, including known predators of hatchling crocodiles (e.g., Clarias sp.) and egg predators (e.g., Mus musculus). This study highlights eDNA metabarcoding as a powerful tool for crocodylian biomonitoring within fragile and diverse ecosystems, particularly where fast, non-invasive methods permit detection in economically important areas and will lead to a better understanding of complex human-crocodile interactions and evaluate habitat suitability for potential reintroductions or recovery programs for threatened crocodylian species.

Patterns of cytonuclear discordance and divergence between subspecies of the scarlet macaw (Ara macao) in Central America.

The scarlet macaw, Ara macao, is a neotropical parrot that contains two described subspecies with broadly discrete geographical distributions. One subspecies, A. m. macao, is found from South America north into southwestern Costa Rica, while the second subspecies, A. m. cyanoptera, is found from eastern Costa Rica north into central Mexico. Our previous research using mitochondrial data to examine phylogeographical divergence across the collective range of these two subspecies concluded that they represent distinct evolutionary entities, with minimal contemporary hybridization between them. Here we further examine phylogenetic relationships and patterns of genetic variation between these two subspecies using a dataset of genetic markers derived from their nuclear genomes. Our analyses show clear nuclear divergence between A. m. macao and A. m. cyanoptera in Central America. Collectively however, samples from this region appear genetically more similar to one another than they do to the examined South American (Brazilian) A. m. macao sample. This observation contradicts our previous assessments based on mitochondrial DNA analyses that A. m. macao in Central and South America represent a single phylogeographical group that is evolutionarily distinct from Central American A. m. cyanoptera. Nonetheless, in agreement with our previous findings, ongoing genetic exchange between the two subspecies appears limited. Rather, our analyses indicate that incomplete lineage sorting is the best supported explanation for cytonuclear discordance within these parrots. High-altitude regions in Central America may act as a reproductive barrier, limiting contemporary hybridization between A. m. macao and A. m. cyanoptera. The phylogeographic complexities of scarlet macaw taxa in this region highlight the need for additional evolutionary examinations of these populations.

High levels of population genetic differentiation in the American crocodile (Crocodylus acutus).

The American crocodile (Crocodylus acutus) is a widely distributed species across coastal and brackish areas of the Neotropical region of the Americas and the Greater Antilles. Available information on patterns of genetic differentiation in C. acutus shows a complex structuring influenced by interspecific interactions (mainly hybridization) and anthropogenic actions (mostly historical hunting, recent poaching, habitat loss and fragmentation, and unintentional translocation of individuals). In this study, we used data on mitochondrial DNA control region and 11 nuclear polymorphic microsatellite loci to assess the degree of population structure of C. acutus in South America, North America, Central America and the Greater Antilles. We used traditional genetic differentiation indices, Bayesian clustering and multivariate methods to create a more comprehensive picture of the genetic relationships within the species across its range. Analyses of mtDNA and microsatellite loci show evidence of a strong population genetic structure in the American crocodile, with unique populations in each sampling locality. Our results support previous findings showing large degrees of genetic differentiation between the continental and the Greater Antillean C. acutus. We report three new haplotypes unique to Venezuela, which are considerably less distant from the Central and North American haplotypes than to the Greater Antillean ones. Our findings reveal genetic population differentiation between Cuban and Jamaican C. acutus and offer the first evidence of strong genetic differentiation among the populations of Greater Antillean C. acutus.

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.

The Chytrid Fungus, Batrachochytrium dendrobatidis, is Widespread Among Cuban Amphibians.

The fungus Batrachochytrium dendrobatidis (Bd) is a generalist amphibian pathogen responsible for chytridiomycosis. It was documented for the first time in Cuba in 2007, the apparent cause of the decline in one species of toad. In a recent survey, Bd was reported only for the highlands of Central Cuba. In the present study, we reexamined the geographic distribution and level of impact of Bd in Cuba by conducting an island-wide sampling in 10 localities and collecting skin swabs from 18 species and 28 environmental samples. We report detection of Bd in 60% of sampled sites and in 58% of sampled taxa. We show that Bd is associated with riparian, arboreal and terrestrial species, and it was estimated to occur in approximately 30% of the aquatic habitats we sampled. In addition, we confirmed that a dying individual of the species Eleutherodactylus casparii was severely infected with Bd. We also rise concern about the endanger toad Peltophryne longinasus and about three species of endemic riparian frogs that were not detected during our surveys. This study demonstrates that this pathogen is widespread throughout Cuba and provides relevant evidence to advance our understanding of its detection in amphibians and the aquatic environment in Cuba and about the occurrence of Bd in species with different ecologies. We provide valuable baseline information for Bd risk assessment and decision-making processes to mitigate its negative impact on Cuban amphibians.

Genetic evidence supports a distinct lineage of American crocodile (Crocodylus acutus) in the Greater Antilles.

Four species of true crocodile (genus Crocodylus) have been described from the Americas. Three of these crocodile species exhibit non-overlapping distributions-Crocodylus intermedius in South America, C. moreletii along the Caribbean coast of Mesoamerica, and C. rhombifer confined to Cuba. The fourth, C. acutus, is narrowly sympatric with each of the other three species. In this study, we sampled 113 crocodiles across Crocodylus populations in Cuba, as well as exemplar populations in Belize and Florida (USA), and sequenced three regions of the mitochondrial genome (D-loop, cytochrome b, cytochrome oxidase I; 3,626 base pair long dataset) that overlapped with published data previously collected from Colombia, Jamaica, and the Cayman Islands. Phylogenetic analyses of these data revealed two, paraphyletic lineages of C. acutus. One lineage, found in the continental Americas, is the sister taxon to C. intermedius, while the Greater Antillean lineage is most closely related to C. rhombifer. In addition to the paraphyly of the two C. acutus lineages, we recovered a 5.4% estimate of Tamura-Nei genetic divergence between the Antillean and continental clades. The reconstructed paraphyly, distinct phylogenetic affinities and high genetic divergence between Antillean and continental C. acutus populations are consistent with interspecific differentiation within the genus and suggest that the current taxon recognized as C. acutus is more likely a complex of cryptic species warranting a reassessment of current taxonomy. Moreover, the inclusion, for the first time, of samples from the western population of the American crocodile in Cuba revealed evidence for continental mtDNA haplotypes in the Antilles, suggesting this area may constitute a transition zone between distinct lineages of C. acutus. Further study using nuclear character data is warranted to more fully characterize this cryptic diversity, resolve taxonomic uncertainty, and inform conservation planning in this system.

Multiple Paternity in a Reintroduced Population of the Orinoco Crocodile (Crocodylus intermedius) at the El Frío Biological Station, Venezuela.

The success of a reintroduction program is determined by the ability of individuals to reproduce and thrive. Hence, an understanding of the mating system and breeding strategies of reintroduced species can be critical to the success, evaluation and effective management of reintroduction programs. As one of the most threatened crocodile species in the world, the Orinoco crocodile (Crocodylus intermedius) has been reduced to only a few wild populations in the Llanos of Venezuela and Colombia. One of these populations was founded by reintroduction at Caño Macanillal and La Ramera lagoon within the El Frío Biological Station, Venezuela. Twenty egg clutches of C. intermedius were collected at the El Frío Biological Station for incubation in the lab and release of juveniles after one year. Analyzing 17 polymorphic microsatellite loci from 335 hatchlings we found multiple paternity in C. intermedius, with half of the 20 clutches fathered by two or three males. Sixteen mothers and 14 fathers were inferred by reconstruction of multilocus parental genotypes. Our findings showed skewed paternal contributions to multiple-sired clutches in four of the clutches (40%), leading to an overall unequal contribution of offspring among fathers with six of the 14 inferred males fathering 90% of the total offspring, and three of those six males fathering more than 70% of the total offspring. Our results provide the first evidence of multiple paternity occurring in the Orinoco crocodile and confirm the success of reintroduction efforts of this critically endangered species in the El Frío Biological Station, Venezuela.

The complete mitochondrial genome structure of the jaguar (Panthera onca).

The jaguar (Panthera onca) is the largest felid in the Western hemisphere, and the only member of the Panthera genus in the New World. The jaguar inhabits most countries within Central and South America, and is considered near threatened by the International Union for the Conservation of Nature. This study represents the first sequence of the entire jaguar mitogenome, which was the only Panthera mitogenome that had not been sequenced. The jaguar mitogenome is 17,049 bases and possesses the same molecular structure as other felid mitogenomes. Bayesian inference (BI) and maximum likelihood (ML) were used to determine the phylogenetic placement of the jaguar within the Panthera genus. Both BI and ML analyses revealed the jaguar to be sister to the tiger/leopard/snow leopard clade.

Integrating molecular, phenotypic and environmental data to elucidate patterns of crocodile hybridization in Belize.

The genus Crocodylus comprises 12 currently recognized species, many of which can be difficult to differentiate phenotypically. Interspecific hybridization among crocodiles is known to occur in captivity and has been documented between some species in the wild. The identification of hybrid individuals is of importance for management and monitoring of crocodilians, many of which are Convention on International Trade in Endangered Species (CITES) listed. In this study, both mitochondrial and nuclear DNA markers were evaluated for their use in confirming a suspected hybrid zone between American crocodile (Crocodylus acutus) and Morelet's crocodile (Crocodylus moreletii) populations in southern Belize where individuals and nests exhibiting atypical phenotypic features had previously been observed. Patterns observed in both phenotypic and molecular data indicate possible behavioural and ecological characteristics associated with hybridization events. The results of the combined analyses found that the majority of suspected hybrid samples represent crosses between female C. acutus and male C. moreletii. Phenotypic data could statistically identify hybrids, although morphological overlap between hybrids and C. moreletii reduced reliability of identification based solely on field characters. Ecologically, C. acutus was exclusively found in saline waters, whereas hybrids and C. moreletii were largely absent in these conditions. A hypothesized correlation between unidirectional hybridization and destruction of C. acutus breeding habitats warrants additional research.

The interplay of homing and dispersal in green turtles: a focus on the southwestern atlantic.

Current understanding of spatial ecology is insufficient in many threatened marine species, failing to provide a solid basis for conservation and management. To address this issue for globally endangered green turtles, we investigated their population distribution by sequencing a mitochondrial control region segment from the Rocas Atoll courtship area (n = 30 males) and four feeding grounds (FGs) in Brazil (n = 397), and compared our findings to published data (n (nesting) = 1205; n (feeding) = 1587). At Rocas Atoll, the first Atlantic courtship area sequenced to date, we found males were differentiated from local juveniles but not from nesting females. In combination with tag data, this indicates possible male philopatry. The most common haplotypes detected at the study sites were CMA-08 and CMA-05, and significant temporal variation was not revealed. Although feeding grounds were differentiated overall, intra-regional structure was less pronounced. Ascension was the primary natal source of the study FGs, with Surinam and Trindade as secondary sources. The study clarified the primary connectivity between Trindade and Brazil. Possible linkages to African populations were considered, but there was insufficient resolution to conclusively determine this connection. The distribution of FG haplotype lineages was nonrandom and indicative of regional clustering. The study investigated impacts of population size, geographic distance, ocean currents, and juvenile natal homing on connectivity, addressed calls for increased genetic sampling in the southwestern Atlantic, and provided data important for conservation of globally endangered green turtles.

Population structure and gene flow of the yellow anaconda (Eunectes notaeus) in northern Argentina.

Yellow anacondas (Eunectes notaeus) are large, semiaquatic boid snakes found in wetland systems in South America. These snakes are commercially harvested under a sustainable management plan in Argentina, so information regarding population structuring can be helpful for determination of management units. We evaluated genetic structure and migration using partial sequences from the mitochondrial control region and mitochondrial genes cyt-b and ND4 for 183 samples collected within northern Argentina. A group of landscape features and environmental variables including several treatments of temperature and precipitation were explored as potential drivers of observed genetic patterns. We found significant population structure between most putative population comparisons and bidirectional but asymmetric migration in several cases. The configuration of rivers and wetlands was found to be significantly associated with yellow anaconda population structure (IBD), and important for gene flow, although genetic distances were not significantly correlated with the environmental variables used here. More in-depth analyses of environmental data may be needed to fully understand the importance of environmental conditions on population structure and migration. These analyses indicate that our putative populations are demographically distinct and should be treated as such in Argentina's management plan for the harvesting of yellow anacondas.

Conservation genetics of harvested river turtles, Podocnemis expansa and Podocnemis unifilis, in the Peruvian Amazon: all roads lead to Iquitos.

We present a mtDNA analysis of Podocnemis expansa (n = 81) and Podocnemis unifilis (n = 228) turtles traded in Peru to evaluate the potential origin of these animals. In particular, we were interested in the relationship between samples reported in the Iquitos markets (IMs) and a Pacaya Samiria Natural Reserve (PSNR) where illegal hunting is presumed. Our mtDNA data showed that, for both species, all haplotypes found within the PSNR were observed in the IM, and that these markets also displayed haplotypes not documented in the reserve. This suggests that the IMs are recipients of Podocnemis turtles from within and outside the PSNR. The fact that most of the haplotype diversity observed in the markets was not found within the PSNR strongly suggests that Podocnemis genetic diversity is exploited in areas where conservation actions are limited. Hence, we recommend expanding Podocnemis conservation efforts outside of protected areas.

Health assessment of the ex situ population of St Vincent parrots (Amazona guildingii) in St Vincent and the Grenadines.

The St Vincent Amazon parrot (Amazona guildingii), listed as a CITES appendix I animal, is endemic only to the mainland of St Vincent and the Grenadines (SVG), Lesser Antilles. A health assessment that included physical examination, hematology, blood parasite evaluation, plasma biochemical analysis, plasma protein electrophoresis, level of exposure to selected infectious agents, and fecal parasites was performed in July 2002 on 37 (51%) of the parrots in the captive population in SVG. Clinical abnormalities noted in the 37 parrots included poor feathers, lipomas, abnormal choanal papillae, obesity, leg ulcers, respiratory abnormalities, cardiac abnormalities, seizure activity, old fractures, missing digits, skin disease associated with mites (Knemidokoptes species), oral granuloma, and a thin, friable beak. Only 7 of the birds were clinically normal on physical examination. Results of hematologic testing, plasma biochemical analysis, and plasma electrophoresis were not statistically different between female and male parrots. No blood parasites were found in any of the 32 samples examined. None of the 36 parrots evaluated had antibodies to the 12 infectious agents tested. Of the 21 fecal samples available, Capillaria species was detected in 1 bird. Findings from this study, in addition to nutritional, genetic, and husbandry evaluations, have been used to make recommendations to the Forestry Department of SVG for improvements in husbandry and veterinary care of this ex situ population.