Evaluating the Phylogenetic Status of the Extinct Japanese Otter on the Basis of Mitochondrial Genome Analysis

Low but highly geographically structured genomic diversity of East Asian Eurasian otters and its conservation implications

Populations of Eurasian otters Lutra lutra, one of the most widely distributed apex predators in Eurasia, have been depleted mainly since the 1950s. However, a lack of information about their genomic diversity and how they are organized geographically in East Asia severely impedes our ability to monitor and conserve them in particular management units. Here, we re-sequenced and analyzed 20 otter genomes spanning continental East Asia, including a population at Kinmen, a small island off the Fujian coast, China. The otters form three genetic clusters (one of L. l. lutra in the north and two of L. l. chinensis in the south), which have diverged in the Holocene. These three clusters should be recognized as three conservation management units to monitor and manage independently. The heterozygosity of the East Asian otters is as low as that of the threatened carnivores sequenced. Historical effective population size trajectories inferred from genomic variations suggest that their low genomic diversity could be partially attributed to changes in the climate since the mid-Pleistocene and anthropogenic intervention since the Holocene. However, no evidence of genetic erosion, mutation load, or high level of inbreeding was detected in the presumably isolated Kinmen Island population. Any future in situ conservation efforts should consider this information for the conservation management units.

Geologically calibrated mammalian tree and its correlation with global events, including the emergence of humans

A robust timetree for Mammalia was constructed using the time calibration function of BEAST v1.10.4 and MEGA 11. The analysis involved the application of times of the most recent common ancestors, including a total of 19 mammalian fossil calibration ages following Benton et al. (Palaeontologia Electronica, 2015, 1-106) for their minimum ages. Additionally, fossil calibration ages for Gorilla, Pan, and a geologic event calibration age for otters were incorporated. Using these calibration ages, I constructed a geologically calibrated tree that estimates the age of the Homo and Pan splitting to be 5.69 Ma. The tree carries several significant implications. First, after the initial rifting at 120 Ma, the Atlantic Ocean expanded by over 500 km around Chron 34 (84 Ma), and vicariant speciation between Afrotheria (Africa) and Xenarthra (South America) appears to have commenced around 70 Ma. Additionally, ordinal level differentiations began immediately following the K-Pg boundary (66.0 Ma), supporting previous hypothesis that mammalian radiation rapidly filled ecological niches left vacant by non-avian dinosaurs. I constructed a diagram depicting the relationship between base substitution rate and age using an additional function in BEAST v1.10.4. The diagram reveals an exponential increase in the base substitution rate approaching recent times. This increased base substitution rate during the Neogene period may have contributed to the expansion of biodiversity, including the extensive adaptive radiation that led to the evolution of Homo sapiens. One significant driving factor behind this radiation could be attributed to the emergence and proliferation of C4 grasses since 20 Ma. These grasses have played a role in increasing carbon fixation, reducing atmospheric CO2 concentration, inducing global cooling, and initiating Quaternary glacial-interglacial cycles, thereby causing significant climatic changes.

Genomics Reveals Complex Population History and Unexpected Diversity of Eurasian Otters (Lutra lutra) in Britain Relative to Genetic Methods

Conservation genetic analyses of many endangered species have been based on genotyping of microsatellite loci and sequencing of short fragments of mtDNA. The increase in power and resolution afforded by whole genome approaches may challenge conclusions made on limited numbers of loci and maternally inherited haploid markers. Here, we provide a matched comparison of whole genome sequencing versus microsatellite and control region (CR) genotyping for Eurasian otters (Lutra lutra). Previous work identified four genetically differentiated "stronghold" populations of otter in Britain, derived from regional populations that survived the population crash of the 1950s-1980s. Using whole genome resequencing data from 45 samples from across the British stronghold populations, we confirmed some aspects of population structure derived from previous marker-driven studies. Importantly, we showed that genomic signals of the population crash bottlenecks matched evidence from otter population surveys. Unexpectedly, two strongly divergent mitochondrial lineages were identified that were undetectable using CR fragments, and otters in the east of England were genetically distinct and surprisingly variable. We hypothesize that this previously unsuspected variability may derive from past releases of Eurasian otters from other, non-British source populations in England around the time of the population bottleneck. Our work highlights that even reasonably well-studied species may harbor genetic surprises, if studied using modern high-throughput sequencing methods.

Analysis of the Complete Mitochondrial Genome of Pteronura brasiliensis and Lontra canadensis

P. brasiliensis and L. canadensis are two otter species, which successfully occupied semi-aquatic habitats and diverged from other Mustelidae. Herein, the full-length mitochondrial genome sequences were constructed for these two otter species for the first time. Comparative mitochondrial genome, selection pressure, and phylogenetic independent contrasts (PICs) analyses were conducted to determine the structure and evolutionary characteristics of their mitochondrial genomes. Phylogenetic analyses were also conducted to confirm these two otter species' phylogenetic position. The results demonstrated that the mitochondrial genome structure of P. brasiliensis and L. canadensis were consistent across Mustelidae. However, selection pressure analyses demonstrated that the evolutionary rates of mitochondrial genome protein-coding genes (PCGs) ND1, ND4, and ND4L were higher in otters than in terrestrial Mustelidae, whereas the evolutionary rates of ND2, ND6, and COX1 were lower in otters. Additionally, PIC analysis demonstrated that the evolutionary rates of ND2, ND4, and ND4L markedly correlated with a niche type. Phylogenetic analysis showed that P. brasiliensis is situated at the base of the evolutionary tree of otters, and then L. canadensis diverged from it. This study suggests a divergent evolutionary pattern of Mustelidae mitochondrial genome PCGs, prompting the otters' adaptation to semi-aquatic habitats.

On the occurrence of Eurasian Otter Lutra lutra (Carnivora: Mustelidae) in Neeru stream of Chenab catchment, Jammu & Kashmir, India

This communication reports the first photographic record of Eurasian Otter Lutra lutra  in a hill stream in Jammu & Kashmir, putting an end to doubts over its presence in the upper Chenab catchment. Three individuals were photographed during a limited camera trap survey conducted in Neeru stream, a left bank tributary of river Chenab during mid-October 2020. We argue that rapid human population influx, infrastructure expansion, and pollution have altered the hydro morphology of Neeru stream, affecting the otter population. This observation calls for more intensive otter surveys in the nearby smaller basins of Neeru, Kalnai, & Sewa and other large tributaries of Chenab River, combining occupancy surveys with camera traps for improved conservation and management of the species in the region.

Taxonomic status of otter species in Nakai-Nam Theun National Park, Lao PDR, based on DNA evidence

Otter populations are threatened by habitat loss, pollution, conflicts with humans, and illegal wildlife trade to meet the demand for pets, for their fur, and for parts used in traditional medicines. Baseline information on the distribution, population genetic diversity, and connectivity is crucial to inform conservation management decisions; however, reliable data from otter populations in Southeast Asia remain scarce. In this study, we conducted baseline otter fecal DNA surveys based on mitochondrial DNA (mtDNA) to identify species, assess the occurrence, and map the spatial distribution of genetic diversity and evolutionary relationships of otter populations using 1700 bp Cytochrome B - Control Region and mitogenome from Nakai-Nam Theun National Park in the Annamite Mountains of Lao PDR. Of the total 56 samples identified to species, the majority (87.5%) was of the widely distributed Eurasian otter with three haplotypes (Lutra lutra; LLLA01-LLLA03), with a calculated haplotype diversity of 0.600 and a nucleotide diversity of 0.00141 based on mitogenome. The second species was the Asian small-clawed otter with only one haplotype detected (Aonyx cinereus; ACLA01). All Eurasian otter haplotypes were newly characterized and clustered within the strongly supported South-Southeast-North Asian clade of Lutra lutra. Compared with the European clade, the high mtDNA diversity of Lutra lutra in Nakai-Nam Theun National Park potentially reflects long-term demographic stability and lesser degree of population bottleneck during the last glacial maxima (LGM, ~21,000 years ago). The single haplotype detected in Asian small-clawed otters had not been detected in previous genetic studies. Our research is the first otter-specific noninvasive genetic study in Lao PDR and provides baseline insights into the otter population diversity in a regional priority site for biodiversity conservation.

Genetic and viability assessment of a reintroduced Eurasian otter Lutra lutra population on the River Ticino, Italy

On the River Ticino in northern Italy, a small number of captive Eurasian otters Lutra lutra, belonging to the European breeding programme for self-sustaining captive populations, were reintroduced in 1997, after the species had been declared locally extinct in the 1980s. We surveyed for otter signs in 2008, 2010, 2016–2017 and 2018, confirming the presence of what is probably a small population. To assess the abundance and viability of the population, we genotyped fresh spraints collected during the last two surveys, using 11 microsatellite markers, and modelled the population trend using Vortex. A minimum of six individuals were identified from 25 faecal samples. The analysis of mitochondrial DNA determined that the reintroduced otters share a transversion that is characteristic of the Asiatic subspecies Lutra lutra barang, confirming the contribution of the Asiatic subspecies to the genetic pool of the captive-bred founder population. Population size was consistent with the release of three pairs of otters and all models implied that the number of founders was too small to ensure the long-term survival of the population. Stochastic factors are therefore likely to threaten the success of this reintroduction.

Evolutionary history of Carnivora (Mammalia, Laurasiatheria) inferred from mitochondrial genomes

The order Carnivora, which currently includes 296 species classified into 16 families, is distributed across all continents. The phylogeny and the timing of diversification of members of the order are still a matter of debate. Here, complete mitochondrial genomes were analysed to reconstruct the phylogenetic relationships and to estimate divergence times among species of Carnivora. We assembled 51 new mitogenomes from 13 families, and aligned them with available mitogenomes by selecting only those showing more than 1% of nucleotide divergence and excluding those suspected to be of low-quality or from misidentified taxa. Our final alignment included 220 taxa representing 2,442 mitogenomes. Our analyses led to a robust resolution of suprafamilial and intrafamilial relationships. We identified 21 fossil calibration points to estimate a molecular timescale for carnivorans. According to our divergence time estimates, crown carnivorans appeared during or just after the Early Eocene Climatic Optimum; all major groups of Caniformia (Cynoidea/Arctoidea; Ursidae; Musteloidea/Pinnipedia) diverged from each other during the Eocene, while all major groups of Feliformia (Nandiniidae; Feloidea; Viverroidea) diversified more recently during the Oligocene, with a basal divergence of Nandinia at the Eocene/Oligocene transition; intrafamilial divergences occurred during the Miocene, except for the Procyonidae, as Potos separated from other genera during the Oligocene.

Aquatic Adaptation and Depleted Diversity: A Deep Dive into the Genomes of the Sea Otter and Giant Otter

Despite its recent invasion into the marine realm, the sea otter (Enhydra lutris) has evolved a suite of adaptations for life in cold coastal waters, including limb modifications and dense insulating fur. This uniquely dense coat led to the near-extinction of sea otters during the 18th-20th century fur trade and an extreme population bottleneck. We used the de novo genome of the southern sea otter (E. l. nereis) to reconstruct its evolutionary history, identify genes influencing aquatic adaptation, and detect signals of population bottlenecks. We compared the genome of the southern sea otter with the tropical freshwater-living giant otter (Pteronura brasiliensis) to assess common and divergent genomic trends between otter species, and with the closely related northern sea otter (E. l. kenyoni) to uncover population-level trends. We found signals of positive selection in genes related to aquatic adaptations, particularly limb development and polygenic selection on genes related to hair follicle development. We found extensive pseudogenization of olfactory receptor genes in both the sea otter and giant otter lineages, consistent with patterns of sensory gene loss in other aquatic mammals. At the population level, the southern sea otter and the northern sea otter showed extremely low genomic diversity, signals of recent inbreeding, and demographic histories marked by population declines. These declines may predate the fur trade and appear to have resulted in an increase in putatively deleterious variants that could impact the future recovery of the sea otter.

What is the taxonomic status of East Asian otter species based on molecular evidence?: focus on the position of the Japanese otter holotype specimen from museum

The Japanese otter (Lutra nippon), once inhabited in most islands of Japan, is now considered as an extinct species. Although the Japanese otter is regarded as a distinct species from the Eurasian otter (L. lutra), its phylogeny and taxonomic status are based on limited information on morphological and genetic data, and thus further clarification is required. Here, we assessed the phylogenetic relationship among the genus Lutra and taxonomic status of L. nippon by using the complete sequences of cytochrome b gene of its holotype. The present phylogenic trees supported that the genus Lutra specimens largely formed monophyletic group, with L. sumatrana as a basal to other Lutra species. Within Lutra species, L. nippon was distantly related with L. lutra. The European otter population of L. l. lutra were clustered together with its subspecies, L. l. chinensis rather than the same subspecies, Korean otter population. The discrepancy between the genetic data and traditional taxonomy justifies the necessity of reexamination of the current subspecific classification system of Eurasian otters. Level of genetic divergence between the holotype of L. nippon and L. lutra was two to three-fold lower than those among the other sister species of the Lutrinae. Based on the level of divergence between the L. nippon and L. lutra, and insufficient evidence of morphological difference between them, it is suggested that designation of Japanese otter as a separate species from L. lutra will be reconsidered.

Identification of novel haplotypes and interpretation of gene flow of mitochondrial DNA control region of Eurasian otter (Lutra lutra) for the effective conservation

The number and distribution of Eurasian otters have declined during twentieth century due to human activity and water pollution. The global conservation status of Eurasian otter is presently ‘Near Threatened (NT)’ and strictly protected by being listed on the international legislation and conventions. A number of studies using the mitochondrial DNA (mtDNA) control region (CR) have been conducted in order to effectively apply conservation and reintroduction programs, especially in Europe. However, aside from Europe, there have been few studies concerning genetic diversity and phylogeny of Eurasian otters. Therefore, in this study, we sequenced partial mtDNA CR sequences (232 bp) from five South Korean Eurasian otters and analyzed 27 otters originating from parts of northeast Asia (South Korea, China, Japan and Russia (Sakhalin)), and Europe. Out of 232 bp partial mtDNA CR sequences, 13 polymorphic sites (5.6%) were identified and 4 novel mtDNA CR haplotypes (Lut16–19) were discovered from 12 Eurasian otters originating from northeast Asian region. In this study, a comprehensive analysis of genetic diversity and population structure of Eurasian otter between Europe and northeast Asia continents were conducted. Of these, different past demographic histories in Pleistocene period might have largely impacted the genetic structure of each population differently. In addition, low degree of gene flow, isolation by distance (IBD) pattern from geographically wide distanced dataset and analysis of molecular variance (AMOVA) also represented distinct genetic characteristics of Eurasian otter between Europe and northeast Asia.

Characterization of the complete mitochondrial genome of Ortleppascaris sinensis (Nematoda: Heterocheilidae) and comparative mitogenomic analysis of eighteen Ascaridida nematodes

Ortleppascaris sinensis (Nematoda: Ascaridida) is a dominant intestinal nematode of the captive Chinese alligator. However, the epidemiology, molecular ecology and population genetics of this parasite remain largely unexplored. In this study, the complete mitochondrial (mt) genome sequence of O. sinensis was first determined using a polymerase chain reaction (PCR)-based primer-walking strategy, and this is also the first sequencing of the complete mitochondrial genome of a member of the genus Ortleppascaris. The circular mitochondrial genome (13,828 bp) of O. sinensis contained 12 protein-coding, 22 transfer RNA and 2 ribosomal RNA genes, but lacked the ATP synthetase subunit 8 gene. Finally, phylogenetic analysis of mtDNAs indicated that the genus Ortleppascaris should be attributed to the family Heterocheilidae. It is necessary to sequence more mtNDAs of Ortleppascaris nematodes in the future to test and confirm our conclusion. The complete mitochondrial genome sequence of O. sinensis reported here should contribute to molecular diagnosis, epidemiological investigations and ecological studies of O. sinensis and other related Ascaridida nematodes.

Geographic variation of craniodental morphology of the Eurasian otter (Lutra lutra) in East Asia

Craniodental morphology of the Eurasian otter (Lutra lutra) in the Korean Peninsula, Japanese islands and Kinmen Island (Taiwan) was studied using geometric morphometrics to identify the skull variations between the populations. Forty adult skulls were examined (29 specimens from the Korean Peninsula, six from Shikoku, Honshu and Hokkaido of Japan, and five from Kinmen Island). Images of the dorsal and ventral views of the skull and the right lateral view of the mandible were analyzed. Specimens from the Korean Peninsula were larger than those from the Japanese islands and Kinmen Island. However, no correlation was observed between the shape variations in the three populations and the centroid size of the skull. The Mann-Whitney U-test showed that relative warps (RWs) RW1, RW2 and RW4 of the dorsal view and RW2 of the ventral view of the skull differed significantly between the populations. Some craniodental differences between the populations were seen in the dorsal and ventral views of the skull, mostly at the snout and parietal regions. The MANOVA test revealed significant differences between the specimens from the Japanese islands and Korean Peninsula and between the specimens from the Korean Peninsula and Kinmen Island. RWs plots showed an overlap of all three populations. In conclusion, the comparisons of the three examined populations revealed significant differences in their craniodental morphology.