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.

Unveiling invasive insect threats to plant biodiversity: Leveraging eDNA metabarcoding and saturated salt trap solutions for biosurveillance.

The negative global impacts of invasive alien species (IAS) on biodiversity are second only to habitat loss. eDNA metabarcoding allows for a faster and more comprehensive evaluation of community species composition, with a higher taxonomic resolution and less taxonomic expertise required than traditional morphological-based biosurveillance. These advantages have positioned eDNA metabarcoding as the standard method for molecular-based detection of invasive alien species, where fast and accurate detectability allows prompt responses to mitigate their adverse effects. Here, eDNA metabarcoding is used for biosurveillance of invasive alien species regulated by Canada in high-risk areas with four main objectives: i) validate the effectiveness of eDNA metabarcoding of salt trap solutions as a molecular technique for IAS detection, ii) compare detection from DNA extracts obtained from filter quarters versus whole filters, iii) benchmark two different bioinformatic pipelines (MetaWorks and mBRAVE), and iv) compare canopy and ground level trapping. eDNA from up to five IAS (Agrilus planipennis, Daktulosphaira vitifoliae, Lymantria dispar, Popillia japonica, and Trichoferus campestris) were successfully detected across years from 2017 to 2022 in southern Ontario, Canada, with successful morphological validation for all except Lymantria dispar and Trichoferus campestris. Analysis of filter quarters in contrast to whole filters was demonstrated to be insufficient for effective IAS detection in each sample. All IAS were detected in only one filter quarter, suggesting a patchy eDNA distribution on the filter. The MetaWorks and mBRAVE bioinformatics pipelines proved effective in identifying IAS, with MetaWorks yielding a higher success rate when comparing molecular and morphological identifications. Ground-level and canopy-level sampling showed differential IAS recovery rates based on the molecular detection, which also varied per collection year, with all found IAS detected at the canopy level in 2022 while only one (Lymantria dispar) in 2020. The present study ratifies the efficacy and importance of eDNA-based detection in a regulatory context and the utility of adding eDNA metabarcoding of saturated salt trap solutions, a critical tool for IAS detection.

Mitochondrial genome sequencing, mapping, and assembly benchmarking for Culicoides species (Diptera: Ceratopogonidae).

BACKGROUND: Mitochondrial genomes are the most sequenced genomes after bacterial and fungal genomic DNA. However, little information on mitogenomes is available for multiple metazoan taxa, such as Culicoides, a globally distributed, megadiverse genus containing 1,347 species. AIM:  Generating novel mitogenomic information from single Culicoides sonorensis and C. biguttatus specimens, comparing available mitogenome mapping and de novo assembly tools, and identifying the best performing strategy and tools for Culicoides species. RESULTS: We present two novel and fully annotated mitochondrial haplotypes for two Culicoides species, C. sonorensis and C. biguttatus. We also annotated or re-annotated the only available reference mitogenome for C. sonorensis and C. arakawae. All species present a high similarity in mitogenome organization. The general gene arrangement for all Culicoides species was identical to the ancestral insect mitochondrial genome. Only short spacers were found in C. sonorensis (up to 30 bp), contrary to C. biguttatus (up to 114 bp). The mitochondrial genes ATP8, NAD2, NAD6, and LSU rRNA exhibited the highest nucleotide diversity and pairwise interspecific p genetic distance, suggesting that these genes might be suitable and complementary molecular barcodes for Culicoides identification in addition to the commonly utilized COI gene. We observed performance differences between the compared mitogenome generation strategies. The mapping strategy outperformed the de novo assembly strategy, but mapping results were partially biased in the absence of species-specific reference mitogenome. Among the utilized tools, BWA performed best for C. sonorensis while SPAdes, MEGAHIT, and MitoZ were among the best for C. biguttatus. The best-performing mitogenome annotator was MITOS2. Additionally, we were able to recover exogenous mitochondrial DNA from Bos taurus (biting midges host) from a C. biguttatus blood meal sample. CONCLUSIONS: Two novel annotated mitogenome haplotypes for C. sonorensis and C. biguttatus using High-Throughput Sequencing are presented. Current results are useful as the baseline for mitogenome reconstruction of the remaining Culicoides species from single specimens to HTS and genome annotation. Mapping to a species-specific reference mitogenome generated better results for Culicoides mitochondrial genome reconstruction than de novo assembly, while de novo assembly resulted better in the absence of a closely related reference mitogenome. These results have direct implications for molecular-based identification of these vectors of human and zoonotic diseases, setting the basis for using the whole mitochondrial genome as a marker in Culicoides identification.

Validation of an Effective Protocol for Culicoides Latreille (Diptera: Ceratopogonidae) Detection Using eDNA Metabarcoding.

eDNA metabarcoding is an effective molecular-based identification method for the biosurveillance of flighted insects. An eDNA surveillance approach maintains specimens for secondary morphological identification useful for regulatory applications. This study identified Culicoides species using eDNA metabarcoding and compared these results to morphological identifications of trapped specimens. Insects were collected using ultraviolet (UV) lighted fan traps containing a saturated salt (NaCl) solution from two locations in Guelph, Ontario, Canada. There were forty-two Culicoides specimens collected in total. Molecular identification detected four species, C. biguttatus, C. stellifer, C. obsoletus, and C. mulrennani. Using morphological identification, two out of these four taxonomic ranks were confirmed at the species level (C. biguttatus and C. stellifer) and one was confirmed at the subgenus level (Avaritia [C. obsoletus]). No molecular detection of Culicoides species occurred in traps with an abundance of less than three individuals per taxon. The inconsistency in identifying Culicoides specimens to the species level punctuates the need for curated DNA reference libraries for Culicoides. In conclusion, the saturated salt (NaCl) solution preserved the Culicoides' morphological characteristics and the eDNA.

Optimization and validation of a cost-effective protocol for biosurveillance of invasive alien species.

Environmental DNA (eDNA) metabarcoding has revolutionized biodiversity monitoring and invasive pest biosurveillance programs. The introduction of insect pests considered invasive alien species (IAS) into a non-native range poses a threat to native plant health. The early detection of IAS can allow for prompt actions by regulating authorities, thereby mitigating their impacts. In the present study, we optimized and validated a fast and cost-effective eDNA metabarcoding protocol for biosurveillance of IAS and characterization of insect and microorganism diversity. Forty-eight traps were placed, following the CFIA's annual forest insect trapping survey, at four locations in southern Ontario that are high risk for forest IAS. We collected insects and eDNA samples using Lindgren funnel traps that contained a saturated salt (NaCl) solution in the collection jar. Using cytochrome c oxidase I (COI) as a molecular marker, a modified Illumina protocol effectively identified 2,535 Barcode Index Numbers (BINs). BINs were distributed among 57 Orders and 304 Families, with the vast majority being arthropods. Two IAS (Agrilus planipennis and Lymantria dispar) are regulated by the Canadian Food Inspection Agency (CFIA) as plant health pests, are known to occur in the study area, and were identified through eDNA in collected traps. Similarly, using 16S ribosomal RNA and nuclear ribosomal internal transcribed spacer (ITS), five bacterial and three fungal genera, which contain species of regulatory concern across several Canadian jurisdictions, were recovered from all sampling locations. Our study results reaffirm the effectiveness and importance of integrating eDNA metabarcoding as part of identification protocols in biosurveillance programs.

Biosurveillance for invasive insect pest species using an environmental DNA metabarcoding approach and a high salt trap collection fluid.

With the increase in global trade and warming patterns, the movement, introduction, and establishment of non-native insect species has increased. A rapid and effective early detection biosurveillance program to identify species of concern is needed to reduce future impacts and costs associated with introduced non-native species. One of the challenges facing insect surveillance trapping methods is the sheer volume of individual specimens in the collections. Although molecular identification methods are improving, they currently have limitations (e.g., destructive processing of specimens) and a protocol addressing these limitations can support regulatory applications that need morphological evidence to corroborate molecular data.The novel protocol presented here uses a metabarcoding approach to amplify environmental DNA from a saturated salt solution trap fluid, which retains trap specimens for downstream morphological identifications. The use of a saturated salt solution to preserve specimens in traps addresses issues with the high evaporation rate of ethanol in traps, and public safety concerns with other fluid preservation options with unattended traps in public settings.Using a metabarcoding approach, a 407-nucleotide segment of the cytochrome c oxidase subunit 1 (COI) animal barcode region was successfully amplified from Lindgren funnel trap collection fluids. These traps were placed in forested areas to survey for wood-boring beetles of regulatory concern. Our results displayed successful amplification of target taxa, including the molecular identification of the Japanese Beetle Popillia japonica, a species regulated in Canada. A second species, Anisandrus maiche, recently introduced to North America, was identified in every trap. The genus Lymantria, which contains numerous species of concern to North American woodlands, was also detected. Also, there were six other species identified of interest due to their potential impacts on native and crop flora and fauna.Our results show how this protocol can be used as an efficient method for the surveillance of insects using a trap with a saturated salt solution and eDNA metabarcoding to detect species of regulatory concern.

Phylogenomics reveals novel relationships among Neotropical crocodiles (Crocodylus spp.).

Extant species in the order Crocodylia are remnants of an ancient lineage of large-bodied archosaur reptiles. Despite decades of systematic studies, phylogenetic relationships among members of the genus Crocodylus (true crocodiles) in the Neotropics are poorly understood. Here we estimated phylogenomic relationships among the four extant Crocodylus species in the Americas. Species-tree reconstructions using genotypic data from 17,538 SNPs collected for 33 individuals spanning six Crocodylus species (four ingroup and two outgroup) revealed novel relationships for all Neotropical species. For the first time, C. acutus, the American crocodile, was recovered as monophyletic when individuals from Antillean and continental populations were analyzed together. Our results also contradict previous inferences based on mitochondrial DNA data and a limited number of nuclear markers by robustly grouping Morelet's crocodile (C. moreletii) as the sister species to C. acutus, suggesting a novel phylogeographic hypothesis for the group. The present study punctuates the importance of using nuclear genome-wide information and representative sampling for resolving phylogenetic relationships, especially in broadly distributed species and those with complex evolutionary histories.

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.

Founded: Genetic Reconstruction of Lineage Diversity and Kinship Informs Ex situ Conservation of Cuban Amazon Parrots (Amazona leucocephala).

Captive breeding is a widespread conservation strategy, yet such programs rarely include empirical genetic data for assessing management assumptions and meeting conservation goals. Cuban Amazon parrots (Amazona leucocephala) are considered vulnerable, and multiple on-island captive populations have been established from wild-caught and confiscated individuals of unknown ancestry. Here, we used mitochondrial haplotypic and nuclear genotypic data at 9 microsatellite loci to quantify the extent and distribution of genetic variation within and among captive populations in Zapata Swamp and Managua, Cuba, and to estimate kinship among breeders (n = 88). Using Bayesian clustering analysis, we detected 2 distinct clusters within the Zapata population, one of which was shared with Managua. Individuals from the cluster unique to Zapata possessed mitochondrial haplotypes with affinities to Cuban subspecies (A. l. leucocephala, A. l. palmarum); the shared cluster was similar, but also included haplotypes closely related to the subspecies restricted to Cayman Brac (A. l. hesterna). Overall mean kinship was low within each captive population (-0.026 to -0.012), with 19 and 11 recommended breeding pairs in Zapata and Managua, respectively, ranked according to mean kinship and informed by molecular sexing. Our results highlight the importance of understanding population history within ex situ management programs, while providing genetic information to directly inform Cuban parrot conservation.

Caracterización morfológica y evaluación de resistencia a Fusarium oxysporum en especies silvestres del género Solanum sección Lycopersicon / Morphological characterization and resistance evaluation to Fusarium oxysporum in wild species Solanum genera Lycopersicon section

Con el objetivo de evaluar la variación morfológica y la resistencia a Fusarium oxysporum en la colección del banco de germoplasma de la Universidad Nacional de Loja (UNL), 146 accesiones de cuatro especies silvestres fueron seleccionadas (Solanum pimpinellifolium, S. neorickii, S. habrochaites, S. lycopersicum var. cerasiforme); y, 20 accesiones de tomate cultivado (S. lycopersicum). Un total de 723 plántulas se transplantaron al campo y se evaluaron morfológicamente durante todo su ciclo biológico, con un descriptor que incluyó 20 caracteres cuantitativos y 25 cualitativos. Se estimó la variabilidad morfológica y se detectaron diferencias fenotípicas relacionadas con el fruto, características vegetativas de las plantas y los componentes del rendimiento. Las variables morfológicas mostraron que S. habrochaites es la especie más diferenciada con relación a las restantes. La podredumbre vascular causada por F. oxysporum es responsable de enormes pérdidas en la producción de tomate en el mundo; sin embargo, las especies silvestres que originaron las variedades cultivadas poseen genes de interés agronómico, utilizados por décadas para la generación de variedades resistentes. Por ello, entre 64 y 228 individuos por especie, se evaluaron por tres procedimientos de inoculación. Todas las especies fueron significativamente diferentes en los niveles de resistencia, tolerancia o susceptibilidad. S. neorickii fue más resistente y tolerante, seguida de S. pimpinellifolium y S. lycopersicum var. cerasiforme, mientras que todas las variedades cultivadas (S. lycopersicum) fueron susceptibles. Tales circunstancias permitieron determinar sin lugar a dudas, que las plantas sobrevivientes y tolerantes poseen un genotipo particular que determina esa característica.

In order to evaluate the morphological variation and Fusarium oxysporum resistance in the National University of Loja (UNL) genebank collection, 146 wild accessions of four species were selected (Solanum pimpinellifolium, S. neorickii, S. habrochaites , S. lycopersicum var. cerasiforme) and 20 accessions of cultivated tomato (S. lycopersicum). For a total of 723 seedlings were transplanted to the field and evaluated morphologically throughout its life cycle, with a descriptor that included 20 quantitative and 25 qualitative characters. Morphological variability was estimated were detected phenotypic differences related the fruit, plant vegetative characteristics and yield components. Morphological variables showed that S. habrochaites is the most differentiated with respect to the others. Vascular rot caused by F. oxysporum is responsible for huge losses in tomato production in the world, but the wild species that gave rise to the cultivated varieties have agronomic interest genes used for decades to generate resistant varieties. Therefore, between 64 and 228 individuals per species were evaluated for three inoculation procedures. All species were significantly different between levels of resistance, susceptibility or tolerance. S. neorickii was more resistant and tolerant, followed by S. pimpinellifolium and S. lycopersicum var. cerasiforme, while all cultivars (S. lycopersicum) were susceptible. Such circumstances allowed undoubtedly determine that survivors and tolerant plants have a particular genotype that determines the characteristic.

Evolutionary history of Cuban crocodiles Crocodylus rhombifer and Crocodylus acutus inferred from multilocus markers.

Among crocodilians, Crocodylus rhombifer is one of the world's most endangered species with the smallest natural distribution. In Cuba, this endemic species coexists with the American crocodile (Crocodylus acutus). Hybridization between these two species is well known in captivity and might occur in the wild, but has never been demonstrated genetically. Here, we combined molecular data with environmental, geographic, and fossil data to infer the evolutionary history of Crocodylus in the Cuban Archipelago, and to evaluate genealogical support for species boundaries. We analyzed seven microsatellite loci plus DNA sequence data from nuclear (RAG-1) and mitochondrial (cytochrome b and cytochrome oxidase I) genes from 89 wild-caught individuals in Cuba, Grand Cayman Island, Jamaica, and Central America, and two samples from zoo collections. Microsatellites showed evidence of introgression, suggesting potential hybridization among Cuban groups. In Cuba, C. acutus contained one mitochondrial DNA (mtDNA) haplotype, whereas C. rhombifer contained two haplotypes. MtDNA data showed that C. acutus is paraphyletic with respect to C. rhombifer, revealing 1% sequence divergence between species within Cuba vs. 8% divergence between Cuban forms and mainland C. acutus. We suggest that hybridization has been a historical as well as a current phenomenon between C. acutus and C. rhombifer. These findings suggest that long-term conservation of crocodiles in Cuba will require identification of genetically pure and hybrid individuals, and a decrease in anthropogenic activities. We also recommend more extensive morphological and genetic analyses of Cuban population to establish clear boundaries of the hybrid zone between C. acutus and C. rhombifer.