Whole genomes show contrasting trends of population size changes and genomic diversity for an Amazonian endemic passerine over the late quaternary.

The "Amazon tipping point" is a global change scenario resulting in replacement of upland terra-firme forests by large-scale "savannization" of mostly southern and eastern Amazon. Reduced rainfall accompanying the Last Glacial Maximum (LGM) has been proposed to have acted as such a tipping point in the past, with the prediction that terra-firme inhabiting species should have experienced reductions in population size as drier habitats expanded. Here, we use whole-genomes of an Amazonian endemic organism (Scale-backed antbirds - Willisornis spp.) sampled from nine populations across the region to test this historical demography scenario. Populations from southeastern Amazonia and close to the Amazon-Cerrado ecotone exhibited a wide range of demographic patterns, while most of those from northern and western Amazonia experienced uniform expansions between 400 kya and 80-60 kya, with gradual declines toward 20 kya. Southeastern populations of Willisornis were the last to diversify and showed smaller heterozygosity and higher runs of homozygosity values than western and northern populations. These patterns support historical population declines throughout the Amazon that affected more strongly lineages in the southern and eastern areas, where historical "tipping point" conditions existed due to the widespread replacement of humid forest by drier and open vegetation during the LGM.

Genomic scan of a healthcare-associated NDM-1-producing Citrobacter freundii ST18 isolated from a green sea turtle impacted by plastic pollution.

BACKGROUND: Carbapenemase-producing Citrobacter freundii has been reported as a leading cause of healthcare-associated infections. Particularly, C. freundii belonging to the sequence type (ST) 18 is considered to be an emerging nosocomial clone. OBJECTIVES: To report the genomic background and phylogenomic analysis of a multidrug-resistant NDM-1-producing C. freundii ST18 (strain CF135931) isolated from an endangered green sea turtle affected by plastic pollution in Brazil. METHODS: Genomic DNA was extracted and sequenced using the Illumina NextSeq platform. De novo assembly was performed by CLC Workbench, and in silico analysis accomplished by bioinformatics tools. For phylogenomic analysis, publicly available C. freundii (txid:546) genome assemblies were retrieved from the NCBI database. RESULTS: The genome size was calculated at 5 290 351 bp, comprising 5263 total genes, 4 rRNAs, 77 tRNAs, 11ncRNAs, and 176 pseudogenes. The strain belonged to C. freundii ST18, whereas resistome analysis predicted genes encoding resistance to ß-lactams (blaNDM-1, blaOXA-1, blaCMY-117, and blaTEM-1C), aminoglycosides (aph(3'')-Ib, aadA16, aph(3')-VI, aac(6')-Ib-cr, and aph(6)-Id), quinolones (aac(6')-Ib-cr), macrolides (mph(A) and erm(B)), sulphonamides (sul1 and sul2), tetracyclines (tetA and tetD), and trimethoprim (dfrA27). The phylogenomic analysis revealed that CF135931 strain is closely related to international human-associated ST18 clones producing NDM-1. CONCLUSION: Genomic surveillance efforts are necessary for robust monitoring of the emergence of drug-resistant strains and WHO critical priority pathogens within a One Health framework. In this regard, this draft genome and associated data can improve understanding of dissemination dynamics of nosocomial clones of carbapenemase-producing C. freundii beyond hospital walls. In fact, the emergence of NDM-1-producing C. freundii of global ST18 in wildlife deserves considerable attention.

The genome sequence of the Loggerhead sea turtle, Caretta caretta Linnaeus 1758.

We present a genome assembly of Caretta caretta (the Loggerhead sea turtle; Chordata, Testudines, Cheloniidae), generated from genomic data from two unrelated females. The genome sequence is 2.13 gigabases in size. The assembly has a busco completion score of 96.1% and N50 of 130.95 Mb. The majority of the assembly is scaffolded into 28 chromosomal representations with a remaining 2% of the assembly being excluded from these.

Metabarcoding of fecal pellets in wild muskox populations reveals negative relationships between microbiome and diet alpha diversity.

Microbiome diversity and diet composition concomitantly influence species health, fitness, immunity, and digestion. In environments where diet varies spatially and temporally, microbiome plasticity may promote rapid host adaptation to available resources. For northern ungulates in particular, metabarcoding of noninvasively collected fecal pellets presents unprecedented insights into their diverse ecological requirements and niches by clarifying the interrelationships of microbiomes, key to deriving nutrients, in context of altered forage availability in changing climates. Muskoxen (Ovibos moschatus) are Arctic-adapted species that experience fluctuating qualities and quantities of vegetation. Geography and seasonality have been noted to influence microbiome composition and diversity in muskoxen, yet it is unclear how their microbiomes intersect with diet. Following observations from other species, we hypothesized increasing diet diversity would result in higher microbiome diversity in muskoxen. We assessed diet composition in muskoxen using three common plant metabarcoding markers and explored correlations with microbiome data. Patterns of dietary diversity and composition were not fully concordant among the markers used, yet all reflected the primary consumption of willows and sedges. Individuals with similar diets had more similar microbiomes, yet in contrast to most literature, yielded negative relationships between microbiome and diet alpha diversity. This negative correlation may reflect the unique capacities of muskoxen to survive solely on high-fiber Arctic forage and provide insight into their resiliency to exploit changing dietary resources in a rapidly warming Arctic altering vegetation diversity.

Tracing Eastern Wolf Origins From Whole-Genome Data in Context of Extensive Hybridization.

Southeastern Canada is inhabited by an amalgam of hybridizing wolf-like canids, raising fundamental questions regarding their taxonomy, origins, and timing of hybridization events. Eastern wolves (Canis lycaon), specifically, have been the subject of significant controversy, being viewed as either a distinct taxonomic entity of conservation concern or a recent hybrid of coyotes (C. latrans) and grey wolves (C. lupus). Mitochondrial DNA analyses show some evidence of eastern wolves being North American evolved canids. In contrast, nuclear genome studies indicate eastern wolves are best described as a hybrid entity, but with unclear timing of hybridization events. To test hypotheses related to these competing findings we sequenced whole genomes of 25 individuals, representative of extant Canadian wolf-like canid types of known origin and levels of contemporary hybridization. Here we present data describing eastern wolves as a distinct taxonomic entity that evolved separately from grey wolves for the past ∼67,000 years with an admixture event with coyotes ∼37,000 years ago. We show that Great Lakes wolves originated as a product of admixture between grey wolves and eastern wolves after the last glaciation (∼8,000 years ago) while eastern coyotes originated as a product of admixture between "western" coyotes and eastern wolves during the last century. Eastern wolf nuclear genomes appear shaped by historical and contemporary gene flow with grey wolves and coyotes, yet evolutionary uniqueness remains among eastern wolves currently inhabiting a restricted range in southeastern Canada.

Evidence of backcross inviability and mitochondrial DNA paternal leakage in sea turtle hybrids.

Hybridization is known to be part of many species' evolutionary history. Sea turtles have a fascinating hybridization system in which species separated by as much as 43 million years are still capable of hybridizing. Indeed, the largest nesting populations in Brazil of loggerheads (Caretta caretta) and hawksbills (Eretmochelys imbricata) have a high incidence of hybrids between these two species. A third species, olive ridleys (Lepidochelys olivacea), is also known to hybridize although at a smaller scale. Here, we used restriction site-associated DNA sequencing (RAD-Seq) markers, mitogenomes, and satellite-telemetry to investigate the patterns of hybridization and introgression in the Brazilian sea turtle population and their relationship with the migratory behaviours between feeding and nesting aggregations. We also explicitly test if the mixing of two divergent genomes in sea turtle hybrids causes mitochondrial paternal leakage. We developed a new species-specific PCR-assay capable of detecting mitochondrial DNA (mtDNA) inheritance from both parental species and performed ultra-deep sequencing to estimate the abundance of each mtDNA type. Our results show that all adult hybrids are first generation (F1) and most display a loggerhead migratory behaviour. We detected paternal leakage in F1 hybrids and different proportions of mitochondria from maternal and paternal species. Although previous studies showed no significant fitness decrease in hatchlings, our results support genetically-related hybrid breakdown possibly caused by cytonuclear incompatibility. Further research on hybrids from other populations in addition to Brazil and between different species will show if backcross inviability and mitochondrial paternal leakage is observed across sea turtle species.

Global phylogeography of ridley sea turtles (Lepidochelys spp.): evolution, demography, connectivity, and conservation.

Globally distributed marine taxa are well suited for investigations of biogeographic impacts on genetic diversity, connectivity, and population demography. The sea turtle genus Lepidochelys includes the wide-ranging and abundant olive ridley (L. olivacea), and the geographically restricted and 'Critically Endangered' Kemp's ridley (L. kempii). To investigate their historical biogeography, we analyzed a large dataset of mitochondrial DNA (mtDNA) sequences from olive (n = 943) and Kemp's (n = 287) ridleys, and genotyped 15 nuclear microsatellite loci in a global sample of olive ridleys (n = 285). We found that the ridley species split ~ 7.5 million years ago, before the Panama Isthmus closure. The most ancient mitochondrial olive ridley lineage, located in the Indian Ocean, was dated to ~ 2.2 Mya. Both mitochondrial and nuclear markers revealed significant structure for olive ridleys between Atlantic (ATL), East Pacific (EP), and Indo-West Pacific (IWP) areas. However, the divergence of mtDNA clades was very recent (< 1 Mya) with low within- clade diversity, supporting a recurrent extinction-recolonization model for these ocean regions. All data showed that ATL and IWP groups were more closely related than those in the EP, with mtDNA data supporting recent recolonization of the ATL from the IWP. Individual olive ridley dispersal between the ATL, EP, and IN/IWP could be interpreted as more male- than female-biased, and genetic diversity was lowest in the Atlantic Ocean. All populations showed signs of recent expansion, and estimated time frames were concordant with their recent colonization history. Investigating species abundance and distribution changes over time is central to evolutionary biology, and this study provides a historical biogeographic context for marine vertebrate conservation and management. Supplementary Information: The online version contains supplementary material available at 10.1007/s10592-022-01465-3.

Divergence and hybridization in sea turtles: Inferences from genome data show evidence of ancient gene flow between species.

Reconstructing past events of hybridization and population size changes are required to understand speciation mechanisms and current patterns of genetic diversity, and ultimately contribute to species' conservation. Sea turtles are ancient species currently facing anthropogenic threats including climate change, fisheries, and illegal hunting. Five of the seven extant sea turtle species are known to currently hybridize, especially along the Brazilian coast where some populations can have ~32%-42% of hybrids. Although frequently observed today, it is not clear what role hybridization plays in the evolutionary diversification of this group of reptiles. In this study, we generated whole genome resequencing data of the five globally distributed sea turtle species to estimate a calibrated phylogeny and the population size dynamics, and to understand the role of hybridization in shaping the genomes of these ancient species. Our results reveal discordant species divergence dates between mitochondrial and nuclear genomes, with a high frequency of conflicting trees throughout the nuclear genome suggesting that some sea turtle species frequently hybridized in the past. The reconstruction of the species' demography showed a general decline in effective population sizes with no signs of recovery, except for the leatherback sea turtle. Furthermore, we discuss the influence of reference bias in our estimates. We show long-lasting ancestral gene flow events within Chelonioidea that continued for millions of years after initial divergence. Speciation with gene flow is a common pattern in marine species, and it raises questions whether current hybridization events should be considered as a part of these species' evolutionary history or a conservation issue.

Author Correction: Population-specific signatures of intra-individual mitochondrial DNA heteroplasmy and their potential evolutionary advantages.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

New Genetic Insights About Hybridization and Population Structure of Hawksbill and Loggerhead Turtles From Brazil.

An extremely high incidence of hybridization among sea turtles is found along the Brazilian coast. This atypical phenomenon and its impact on sea turtle conservation can be elucidated through research focused on the evolutionary history of sea turtles. We assessed high-quality multilocus haplotypes of 143 samples of the 5 species of sea turtles that occur along the Brazilian coast to investigate the hybridization process and the population structure of hawksbill (Eretmochelys imbricata) and loggerhead turtles (Caretta caretta). The multilocus data were initially used to characterize interspecific hybrids. Introgression (F2 hybrids) was only confirmed in hatchlings of F1 hybrid females (hawksbill × loggerhead), indicating that introgression was either previously overestimated and F2 hybrids may not survive to adulthood, or the first-generation hybrid females nesting in Brazil were born as recent as few decades ago. Phylogenetic analyses using nuclear markers recovered the mtDNA-based Indo-Pacific and Atlantic lineages for hawksbill turtles, demonstrating a deep genetic divergence dating from the early Pliocene. In addition, loggerhead turtles that share a common feeding area and belong to distinct Indo-Pacific and Atlantic mtDNA clades present no clear genetic differentiation at the nuclear level. Finally, our results indicate that hawksbill and loggerhead rookeries along the Brazilian coast are likely connected by male-mediated gene flow.

Combined use of mitochondrial and nuclear genetic markers further reveal immature marine turtle hybrids along the South Western Atlantic.

Marine turtle hybridization is usually sporadic and involves reports of only a few individuals; however, Brazilian populations have high hybridization rates. Here we investigated the presence of hybrids in morphologically identified immature hawksbills (Eretmochelys imbricata) along the South Western Atlantic (SWA). We sequenced one mitochondrial (D-Loop) and three nuclear DNA (RAG1, RAG2, and CMOS) markers to better understand the patterns and characteristics of hybrids. We identified 22 hybrids (n = 270), 11 of them at the extreme South of the SWA. Uruguay had the highest hybrid frequency in the SWA (~37.5%) followed by southern Brazil with 30%. These are common areas for loggerheads (Caretta caretta) but uncommon for hawksbills, and these hybrids may be adopting the behavior of loggerheads. By analyzing nuclear markers, we can infer that 50% of the sampled hybrids are first generation (F1) and 36% are the result of backcrosses between hybrids and pure E. imbricata (> F1). We also report for the first time immature E. imbricata x Lepidochelys olivacea hybrids at the Brazilian coast. Considering the high frequency of hybrids in the SWA, continuous monitoring should be performed to assess the fitness, genetic integrity, and extent of changes in the gene pools of involved populations.

High prevalence of subclinical frog virus 3 infection in freshwater turtles of Ontario, Canada.

Ranaviruses have been associated with chelonian mortality. In Canada, the first two cases of ranavirus were detected in turtles in 2018 in Ontario, although a subsequent survey of its prevalence failed to detect additional positive cases. To confirm the prevalence of ranavirus in turtles in Ontario, we used a more sensitive method to investigate if lower level persistent infection was present in the population. Here we report results via a combination of qPCR, PCR, Sanger sequencing and genome sequencing from turtles from across Ontario, with no clinical signs of illness. We found 2 positives with high viral load and 5 positives with low viral load. Histopathology found subtle histological changes. DNA sequences identified two types of frog virus 3 (FV3), and genome sequencing identified a ranavirus similar to wild-type FV3. Our results show that the virus has been present in Ontario's turtles as subclinical infections.

Complete mitochondrial genome of the Florida manatee (Trichechus manatus latirotris, Sirenia).

The Florida manatee (Trichechus manatus latirostris) is an endangered subspecies of the West Indian manatee (T. manatus), which inhabits inland and marine waters of southeastern United States. In this study, we assembled the mitochondrial genome (mtDNA) of the Florida manatee from whole genome shotgun reads. As a result, we show that the currently annotated T. manatus mtDNA belongs to a different species, the Amazonian manatee (T. inunguis). The newly assembled Florida manatee mtDNA is 16,881 bp in length, with 13 protein-coding genes, two ribosomal RNAs (rRNAs), 22 transfer RNAs (tRNAs) and one non-coding control region (D-loop). Phylogenetic analysis based on the control region indicates the newly assembled mtDNA is haplotype A01, characteristic of T. m. latirostris, while the current mtDNA associated with the Florida manatee genome assembly has a Ti02 haplotype that is found in Amazonian manatees and hybrids.

Population-specific signatures of intra-individual mitochondrial DNA heteroplasmy and their potential evolutionary advantages.

Heteroplasmy is the existence of more than one mitochondrial DNA (mtDNA) variant within a cell. The evolutionary mechanisms of heteroplasmy are not fully understood, despite being a very common phenomenon. Here we combined heteroplasmy measurements using high throughput sequencing on green turtles (Chelonia mydas) with simulations to understand how heteroplasmy modulates population diversity across generations and under different demographic scenarios. We found heteroplasmy to be widespread in all individuals analysed, with consistent signal in individuals across time and tissue. Significant shifts in haplotype composition were found from mother to offspring, signalling the effect of the cellular bottleneck during oogenesis as included in the model. Our model of mtDNA inheritance indicated that heteroplasmy favoured the increase of population diversity through time and buffered against population bottlenecks, thus indicating the importance of this phenomenon in species with reduced population sizes and frequent population bottlenecks like marine turtles. Individuals with recent haplotypes showed higher levels of heteroplasmy than the individuals with ancient haplotypes, suggesting a potential advantage of maintaining established copies when new mutations arise. We recommend using heteroplasmy through high throughput sequencing in marine turtles, as well as other wildlife populations, for diversity assessment, population genetics, and mixed stock analysis.

Geography, seasonality, and host-associated population structure influence the fecal microbiome of a genetically depauparate Arctic mammal.

The Canadian Arctic is an extreme environment with low floral and faunal diversity characterized by major seasonal shifts in temperature, moisture, and daylight. Muskoxen (Ovibos moschatus) are one of few large herbivores able to survive this harsh environment. Microbiome research of the gastrointestinal tract may hold clues as to how muskoxen exist in the Arctic, but also how this species may respond to rapid environmental changes. In this study, we investigated the effects of season (spring/summer/winter), year (2007-2016), and host genetic structure on population-level microbiome variation in muskoxen from the Canadian Arctic. We utilized 16S rRNA gene sequencing to characterize the fecal microbial communities of 78 male muskoxen encompassing two population genetic clusters. These clusters are defined by Arctic Mainland and Island populations, including the following: (a) two mainland sampling locations of the Northwest Territories and Nunavut and (b) four locations of Victoria Island. Between these geographic populations, we found that differences in the microbiome reflected host-associated genetic cluster with evidence of migration. Within populations, seasonality influenced bacterial diversity with no significant differences between years of sampling. We found evidence of pathogenic bacteria, with significantly higher presence in mainland samples. Our findings demonstrate the effects of seasonality and the role of host population-level structure in driving fecal microbiome differences in a large Arctic mammal.

Frog Virus 3 Genomes Reveal Prevalent Recombination between Ranavirus Lineages and Their Origins in Canada.

Ranaviruses are pathogens associated with the decline of amphibian populations across much of their distribution. In North America, frog virus 3 (FV3) is a widely distributed pathogen with wild populations of amphibians harboring different lineages and putative recombinants between FV3 and common midwife toad virus (CMTV). These recombinants have higher pathogenicity, and CMTV-derived genes associated with virulence are reported in wild strains in Canada. However, while FV3 is linked to amphibian die-offs in North America, CMTVs have been reported only in commercial frog farms in North America. We sequenced complete genomes of 18 FV3 isolates from three amphibian species to characterize genetic diversity of the lineages in Canada and infer possible recombinant regions. The 18 FV3 isolates displayed different signals of recombination, varying from none to interspersed recombination with previously isolated CMTV-like viruses. In general, most recombination breakpoints were located within open reading frames (ORFs), generating new ORFs and proteins that were a mixture between FV3 and CMTV. A combined spatial and temporal phylogeny suggests the presence of the FV3 lineage in Canada is relatively contemporary (<100 years), corroborating the hypothesis that both CMTV- and FV3-like viruses spread to North America when the international commercial amphibian trade started. Our results highlight the importance of pathogen surveillance and viral dynamics using full genomes to more clearly understand the mechanisms of disease origin and spread.IMPORTANCE Amphibian populations are declining worldwide, and these declines have been linked to a number of anthropogenic factors, including disease. Among the pathogens associated with amphibian mortality, ranaviruses have caused massive die-offs across continents. In North America, frog virus 3 (FV3) is a widespread ranavirus that can infect wild and captive amphibians. In this study, we sequenced full FV3 genomes isolated from frogs in Canada. We report widespread recombination between FV3 and common midwife toad virus (CMTV). Phylogenies indicate a recent origin for FV3 in Canada, possibly as a result of international amphibian trade.

Low intraspecific variation of Frog virus 3 with evidence for novel FV3-like isolates in central and northwestern Canada.

Frog virus 3 (FV3) and FV3-like ranaviruses can infect a variety of cold-blooded aquatic species and present a primary threat to amphibians across the globe. Previous studies of FV3-like viruses have largely investigated higher-level phylogenetic distinctions of these pathogens via portions of the conserved major capsid protein (MCP), and the putative virulence gene vIF-2α. Few studies, however, have investigated the spatial distribution of FV3 variants at the population level3-data that can be used to further understand the spatial epidemiology of this disease. In this study, we sequenced the MCP and vIF-2α of 127 FV3-positive amphibians sampled from Canadian water bodies in Ontario, northeastern Alberta, and southern Northwest Territories to explore whether intraspecific genetic variation exists within FV3. There was a lack of variation at the 2 markers across these regions, suggesting that there is a lack of FV3 sequence diversity in Canada, which may hint at a single source of infection that has spread. However, an undocumented variant termed Wood Buffalo ranavirus (WBRV) was detected in samples from 3 sites in Alberta and Northwest Territories that clustered within the FV3-like lineage with 99.3% sequence homology for MCP. For vIF-2α, all sequences were the expected truncated variant except for 6 samples in Ontario. These latter sequences were suggestive of recombination with common midwife toad virus (CMTV). The lack of variation suggests that higher-resolution genome analyses will be required to further explore the spatial spread and intraspecific variation of the disease.

Demography and selection shape transcriptomic divergence in field crickets.

Gene flow, demography, and selection can result in similar patterns of genomic variation and disentangling their effects is key to understanding speciation. Here, we assess transcriptomic variation to unravel the evolutionary history of Gryllus rubens and Gryllus texensis, cryptic field cricket species with highly divergent mating behavior. We infer their demographic history and screen their transcriptomes for footprints of selection in the context of the inferred demography. We find strong support for a long history of bidirectional gene flow, which ceased during the late Pleistocene, and a bottleneck in G. rubens consistent with a peripatric origin of this species. Importantly, the demographic history has likely strongly shaped patterns of genetic differentiation (empirical FST distribution). Concordantly, FST -based selection detection uncovers a large number of outliers, likely comprising many false positives, echoing recent theoretical insights. Alternative genetic signatures of positive selection, informed by the demographic history of the sibling species, highlighted a smaller set of loci; many of these are candidates for controlling variation in mating behavior. Our results underscore the importance of demography in shaping overall patterns of genetic divergence and highlight that examining both demography and selection facilitates a more complete understanding of genetic divergence during speciation.

DNA sampling from eggshells and microsatellite genotyping in rare tropical birds: Case study on Brazilian Merganser

Abstract This study shows that sampling maternal DNA from hatched and abandoned eggshells is a viable noninvasive strategy for studying the genetics of rare or endangered tropical birds, as exemplified here by the Brazilian Merganser (Mergus octosetaceus). Eighteen microsatellites were isolated from enriched libraries and nine heterologous loci from related species were tested. Seven loci were amplified successfully, with five of them being polymorphic. These loci exhibited amplicons ranging from 110 to 254 bp for 132 samples, with 60 from eggshells and 72 from blood or muscle samples. The number of alleles for M. octosetaceus ranged from one to six (mean = 3.71), which is low compared to M. merganser (1-15 alleles), a 'least concern' species. Genetic diversity did not differ significantly between noninvasive and invasive samples (Z(u) = 0.31, p = 0.37). Thus, noninvasive sampling, as demonstrated here with eggshells, provides an efficient means to assess genetic diversity in tropical birds without the need to capture and handle them.

DNA sampling from eggshells and microsatellite genotyping in rare tropical birds: Case study on Brazilian Merganser.

This study shows that sampling maternal DNA from hatched and abandoned eggshells is a viable noninvasive strategy for studying the genetics of rare or endangered tropical birds, as exemplified here by the Brazilian Merganser (Mergus octosetaceus). Eighteen microsatellites were isolated from enriched libraries and nine heterologous loci from related species were tested. Seven loci were amplified successfully, with five of them being polymorphic. These loci exhibited amplicons ranging from 110 to 254 bp for 132 samples, with 60 from eggshells and 72 from blood or muscle samples. The number of alleles for M. octosetaceus ranged from one to six (mean = 3.71), which is low compared to M. merganser (1-15 alleles), a 'least concern' species. Genetic diversity did not differ significantly between noninvasive and invasive samples (Z(u) = 0.31, p = 0.37). Thus, noninvasive sampling, as demonstrated here with eggshells, provides an efficient means to assess genetic diversity in tropical birds without the need to capture and handle them.