Structure of the mitochondrial control region of the Eurasian otter (Lutra lutra; Carnivora, Mustelidae): patterns of genetic heterogeneity and implications for conservation of the species in Italy

Characterisation of the Complete Mitochondrial Genome of Critically Endangered Mustela lutreola (Carnivora: Mustelidae) and Its Phylogenetic and Conservation Implications

In this paper, a complete mitochondrial genome of the critically endangered European mink Mustela lutreola L., 1761 is reported. The mitogenome was 16,504 bp in length and encoded the typical 13 protein-coding genes, two ribosomal RNA genes and 22 transfer RNA genes, and harboured a putative control region. The A+T content of the entire genome was 60.06% (A > T > C > G), and the AT-skew and GC-skew were 0.093 and −0.308, respectively. The encoding-strand identity of genes and their order were consistent with a collinear gene order characteristic for vertebrate mitogenomes. The start codons of all protein-coding genes were the typical ATN. In eight cases, they were ended by complete stop codons, while five had incomplete termination codons (TA or T). All tRNAs had a typical cloverleaf secondary structure, except tRNA^Ser(AGC) and tRNA^Lys, which lacked the DHU stem and had reduced DHU loop, respectively. Both rRNAs were capable of folding into complex secondary structures, containing unmatched base pairs. Eighty-one single nucleotide variants (substitutions and indels) were identified. Comparative interspecies analyses confirmed the close phylogenetic relationship of the European mink to the so-called ferret group, clustering the European polecat, the steppe polecat and the black-footed ferret. The obtained results are expected to provide useful molecular data, informing and supporting effective conservation measures to save M. lutreola.

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.

Complete Mitochondrial Genome of Three Species of the Genus Microtus (Arvicolinae, Rodentia)

The 65 species of the genus Microtus have unusual sex-related genetic features and a high rate of karyotype variation. However, only nine complete mitogenomes for these species are currently available. We describe the complete mitogenome sequences of three Microtus, which vary in length from 16,295 bp to 16,331 bp, contain 13 protein-coding genes (PCGs), two ribosomal RNA genes, 22 transfer RNA genes and a control region. The length of the 13 PCGs and the coded proteins is the same in all three species, and the start and stop codons are conserved. The non-coding regions include the L-strand origin of replication, with the same sequence of 35 bp, and the control region, which varies between 896 bp and 930 bp in length. The control region includes three domains (Domains I, II and III) with extended termination-associated sequences (ETAS-1 and ETAS-2) in Domain I. Domain II and Domain III include five (CSB-B, C, D, E and F) and three (CSB-1, CSB-2, and CSB-3) conserved sequence blocks, respectively. Phylogenetic reconstructions using the mitochondrial genomes of all the available Microtus species and one representative species from another genus of the Arvicolinae subfamily reproduced the established phylogenetic relationships for all the Arvicolinae genera that were analyzed.

The complete mitochondrial genome of Talpa aquitania (Talpidae; Insectivora), a mole species endemic to northern Spain and southern France

The complete mitogenome sequence of Talpa aquitania, a recently described Talpa species, was assembled using whole-genome sequencing data. It varies in length from 16,776 to 16,846 bp, contains 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes, one origin of L-strand replication, and a control region. In the control region, which varied from 1320 to 1390 bp, we identified the extended termination-associated sequence (ETAS-1 and ETAS-2) and the conserved sequence blocks (CSB-1, 2, 3, B, C, D, E, F). In addition, this region includes a 10 bp tandem repeat DNA sequence, with a variable number of repeats that suggest the existence of heteroplasmy. Phylogeny reconstructions based on Maximum Likelihood, Neighbor-joining and Bayesian inference analyses yielded phylogenies with similar topologies demonstrating that T. aquitania and T. occidentalis are sister species.

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.

Mitochondrial heteroplasmy in vertebrates using ChIP-sequencing data

BACKGROUND

Mitochondrial heteroplasmy, the presence of more than one mitochondrial DNA (mtDNA) variant in a cell or individual, is not as uncommon as previously thought. It is mostly due to the high mutation rate of the mtDNA and limited repair mechanisms present in the mitochondrion. Motivated by mitochondrial diseases, much focus has been placed into studying this phenomenon in human samples and in medical contexts. To place these results in an evolutionary context and to explore general principles of heteroplasmy, we describe an integrated cross-species evaluation of heteroplasmy in mammals that exploits previously reported NGS data. Focusing on ChIP-seq experiments, we developed a novel approach to detect heteroplasmy from the concomitant mitochondrial DNA fraction sequenced in these experiments.

RESULTS

We first demonstrate that the sequencing coverage of mtDNA in ChIP-seq experiments is sufficient for heteroplasmy detection. We then describe a novel detection method for accurate detection of heteroplasmies, which also accounts for the error rate of NGS technology. Applying this method to 79 individuals from 16 species resulted in 107 heteroplasmic positions present in a total of 45 individuals. Further analysis revealed that the majority of detected heteroplasmies occur in intergenic regions.

CONCLUSION

In addition to documenting the prevalence of mtDNA in ChIP-seq data, the results of our mitochondrial heteroplasmy detection method suggest that mitochondrial heteroplasmies identified across vertebrates share similar characteristics as found for human heteroplasmies. Although largely consistent with previous studies in individual vertebrates, our integrated cross-species analysis provides valuable insights into the evolutionary dynamics of mitochondrial heteroplasmy.

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

The Japanese otter lived throughout four main Japanese islands, but it has not been observed in the wild since 1979 and was declared extinct in 2012. Although recent taxonomic and molecular phylogenetic studies suggest that it should be treated as an independent species, International Union for Conservation of Nature Red List considers it as subspecies of Lutra lutra. Therefore, the taxonomic status of this species needs to be resolved. Here we determined the complete mitochondrial genome of two Japanese otters caught in Kanagawa and Kochi prefectures and five Eurasian otters (L. lutra). We reconstructed a molecular phylogenetic tree to estimate the phylogenetic position of the Japanese otter in Lutrinae using the Japanese otters and the other 11 Lutrinae species on the basis of ND5 (692 bp) and cytochrome b (1,140 bp) sequences. We observed that the two Japanese otters had close relationships with Eurasian otters, forming a monophyletic group (100% bootstrap probability). To elucidate detailed phylogenetic relationships among the Japanese and Eurasian otters, we reconstructed a maximum likelihood tree according to mitochondrial genome sequences (14,740 bp). The Japanese otter (JO1) collected in Kanagawa was deeply nested in the Eurasian otter clade, whereas the Japanese otter (JO2) collected in Kochi formed a distinct independent lineage in the Lutra clade. The estimated molecular divergences time for the ancestral lineages of the Japanese otters was 0.10 Ma (95%: 0.06-0.16 Ma) and 1.27 Ma (95%: 0.98-1.59 Ma) for JO1 and JO2 lineages, respectively. Thus, JO1 was identified as a member of L. lutra; JO2 represented the old Japanese otter lineage, which may be a distinct new species or subspecies of Lutra. We suggest that the ancestral population of the JO2 lineage migrated to Japan via the land bridge that existed between western Japanese islands and Asian continent at 1.27 Ma.

Molecular characteristics and evolution of the mitochondrial control region in three genera (Hipposideridae: Hipposideros Aselliscus and Coelops) of leaf-nosed bats

The mitochondrial control region (CR) was sequenced for three genera of Hipposideridae to give a detailed overview of its features. The CR of leaf-nosed bats (1288-1560 bp) was divided into three domains like that of other mammals. In addition to the common conserved blocks (ETAS1, ETAS2, F-B boxes, CSB1, CSB2, and CSB3) found in all species, a CSB1-like element was also detected in the conserved sequence blocks (CSB). Repeated motifs were examined in the ETAS of Aselliscus stoliczkanus (26 bp) and Hipposideros bicolor (80 bp) and were present in the CSB of all individuals (6, 8, 16, and 20 bp). Phylogenetic reconstructions using the CR sequences indicated that the phylogenetic relationships among Hipposideros species were consistent with the results of other molecular and phenetic analyses. Aselliscus and Coelops had a closer relationship. But the central domain could not be used for phylogenetic analyses at family and genus levels due to its high conservation.