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.

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.

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.

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.

Two interbreeding populations of Saccharomyces cerevisiae strains coexist in cachaça fermentations from Brazil.

In this study, the phylogenetic relationships between cachaça strains of Saccharomyces cerevisiae isolated from different geographical areas in Brazil were obtained on the basis of sequences of one mitochondrial (COX2) and three nuclear (EGT2, CAT8, and BRE5) genes. This analysis allowed us to demonstrate that different types of strains coexist in cachaça fermentations: wine strains, exhibiting alleles related or identical to those present in European wine strains; native strains, containing alleles similar to those found in strains isolated from traditional fermentations from Latin America, North America, Malaysian, Japan, or West Africa; and their intraspecific hybrids or 'mestizo' strains, heterozygous for both types of alleles. Wine strains and hybrids with high proportions of wine-type alleles predominate in southern and southeastern Brazil, where cachaça production coexists with winemaking. The high frequency of 'wine-type' alleles in these regions is probably due to the arrival of wine immigrant strains introduced from Europe in the nearby wineries due to the winemaking practices. However, in north and northeastern states, regions less suited or not suited for vine growing and winemaking, wine-type alleles are much less frequent because 'mestizo' strains with intermediate or higher proportions of 'native-type' alleles are predominant.

Genetic variation and population structure in the endangered Hermann's tortoise: the roles of geography and human-mediated processes.

The Hermann's tortoise (Testudo hermanni) is an endangered land tortoise distributed in disjoint populations across Mediterranean Europe. We investigated its genetic variation by typing 1 mitochondrial locus and 9 nuclear microsatellites in approximately 300 individuals from 22 localities. Our goal was to understand the relative impact of natural and human-mediated processes in shaping the genetic structure and to identify the genetic priorities for the conservation of this species. We found that 1) all geographic areas are highly differentiated, mainly as a function of their distance but with a clear genetic discontinuity (F st values larger than 0.4) between the Eastern and the Western subspecies; 2) the contact zone between subspecies is located farthest to the west than previously believed, and it probably coincides with the delta of the largest Italian river; 3) extinction events due to climatic conditions in the Upper Palaeolithic and subsequent human-mediated translocations in the Neolithic possibly explain the unexpected similarity among Spain, Sicily, and Corsica. For conservation purposes, the large majority of genetic pools appears native although hybridization among subspecies, related to extensive 20th century trade of tortoises across Europe, is observed in Spain and some Italian samples. Most populations do not seem at immediate risk of low genetic variation, except the French population, which has very low nuclear genetic diversity (heterozygosity = 0.25) and where 50 out of 51 sampled animals shared the same mitochondrial sequence. In general, restocking and reintroduction plans should carefully consider the genetic background of the individuals.

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.

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.

Nuclear markers reveal a complex introgression pattern among marine turtle species on the Brazilian coast.

Surprisingly, a high frequency of interspecific sea turtle hybrids has been previously recorded in a nesting site along a short stretch of the Brazilian coast. Mitochondrial DNA data indicated that as much as 43% of the females identified as Eretmochelys imbricata are hybrids in this area (Bahia State of Brazil). It is a remarkable find, because most of the nesting sites surveyed worldwide, including some in northern Brazil, presents no hybrids, and rare Caribbean sites present no more than 2% of hybrids. Thus, a detailed understanding of the hybridization process is needed to evaluate natural or anthropogenic causes of this regional phenomenon in Brazil, which could be an important factor affecting the conservation of this population. We analysed a set of 12 nuclear markers to investigate the pattern of hybridization involving three species of sea turtles: hawksbill (E. imbricata), loggerhead (Caretta caretta) and olive ridley (Lepidochelys olivacea). Our data indicate that most of the individuals in the crossings L. olivacea × E. imbricata and L. olivacea × C. caretta are F1 hybrids, whereas C. caretta × E. imbricata crossings present F1 and backcrosses with both parental species. In addition, the C. caretta × E. imbricata hybridization seems to be gender and species biased, and we also found one individual with evidence of multispecies hybridization among C. caretta × E. imbricata × Chelonia mydas. The overall results also indicate that hybridization in this area is a recent phenomenon, spanning at least two generations or ~40 years.

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.

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.

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.

Geographic patterns of genetic variation in a broadly distributed marine vertebrate: new insights into loggerhead turtle stock structure from expanded mitochondrial DNA sequences.

Previous genetic studies have demonstrated that natal homing shapes the stock structure of marine turtle nesting populations. However, widespread sharing of common haplotypes based on short segments of the mitochondrial control region often limits resolution of the demographic connectivity of populations. Recent studies employing longer control region sequences to resolve haplotype sharing have focused on regional assessments of genetic structure and phylogeography. Here we synthesize available control region sequences for loggerhead turtles from the Mediterranean Sea, Atlantic, and western Indian Ocean basins. These data represent six of the nine globally significant regional management units (RMUs) for the species and include novel sequence data from Brazil, Cape Verde, South Africa and Oman. Genetic tests of differentiation among 42 rookeries represented by short sequences (380 bp haplotypes from 3,486 samples) and 40 rookeries represented by long sequences (∼800 bp haplotypes from 3,434 samples) supported the distinction of the six RMUs analyzed as well as recognition of at least 18 demographically independent management units (MUs) with respect to female natal homing. A total of 59 haplotypes were resolved. These haplotypes belonged to two highly divergent global lineages, with haplogroup I represented primarily by CC-A1, CC-A4, and CC-A11 variants and haplogroup II represented by CC-A2 and derived variants. Geographic distribution patterns of haplogroup II haplotypes and the nested position of CC-A11.6 from Oman among the Atlantic haplotypes invoke recent colonization of the Indian Ocean from the Atlantic for both global lineages. The haplotypes we confirmed for western Indian Ocean RMUs allow reinterpretation of previous mixed stock analysis and further suggest that contemporary migratory connectivity between the Indian and Atlantic Oceans occurs on a broader scale than previously hypothesized. This study represents a valuable model for conducting comprehensive international cooperative data management and research in marine ecology.