Research Article |
Corresponding author: Kenji Saitoh ( ksaitoh@affrc.go.jp ) Academic editor: Szabolcs Lengyel
© 2017 Kenji Saitoh, Noriyasu Suzuki, Masumi Ozaki, Kazuhiro Ishii, Tetsuya Sado, Takahiro Morosawa, Takatoshi Tsunagawa, Masaru Tsuchiya.
This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Citation:
Saitoh K, Suzuki N, Ozaki M, Ishii K, Sado T, Morosawa T, Tsunagawa T, Tsuchiya M (2017) Natural habitats uncovered? – Genetic structure of known and newly found localities of the endangered bitterling Pseudorhodeus tanago (Cyprinidae). Nature Conservation 17: 19-33. https://doi.org/10.3897/natureconservation.17.10939
|
Overuse of natural resources by humans is a major threat to biodiversity. Overuse often involves species of economic or esthetic value, and fish are a typical example for a group that is exploited both for economic reasons (for human consumption) and for esthetic reasons (e.g. by aquarists). Pseudorhodeus tanago (Tanaka, 1909) (formerly known as Tanakia tanago) is a small colorful but legally protected (fishing, keeping and transfer are banned) bitterling fish distributed around Tokyo, Japan. Whereas it is critically endangered and more and more habitat loss has occurred, at least four stocks have been newly found during the last decade. To explore whether emergence of these newly found habitats is a consequence of incomplete survey, we genotyped mitochondrial cytochrome b sequence of P. tanago from 17 localities and an illegal home aquarium. Populations known by the past extensive survey (13 localities) showed geographically structured population genetic characteristics. Population-specific haplotypes were common indicating past divergence and bottleneck events. Four (north, {center + west}, south_1, south_2) or five (north, center, west, south_1, south_2) geographic groups were detectable as for these known localities. On the other hand, newly found stocks were polymorphic and showed identical haplotypes from distant known localities. If we assume historical basis of distribution and genetic characteristics of these newly found stocks, it must be a series of unlikely geological events and haplotype sorting. We discuss potential issues posed by these questionable stocks.
Bottleneck, bucket biology, conservation genetics, fish dumping, home aquarium, phylogeography, poaching, satoyama, Tanakia tanago
Overuse of organisms by hunting or fishing for trade or esthetic purposes is one of the biggest threats to biodiversity. Controlling these activities is fundamental for ensuring the persistence of endangered organisms, particularly those that have traits attractive to humans. Alerts to poaching are thus necessary on those organisms legally protected. Not only tusks of elephants and rhinos but beautifully colored bodies of fish as well have attracted violators. The red list of Japanese brackish and freshwater fish (
Pseudorhodeus tanago (Tanaka, 1909) (formerly known as Tanakia tanago transferred to the new genus,
Records of 45 localities of P. tanago habitats were established as part of an extensive survey (
From the habitat characteristics of P. tanago, confined and scattered among headwaters of fine branching dales in agricultural landscape on hill terrains, we assume geographically structured population genetic characteristics of this species.
Whereas more and more local populations have been lost recently, P. tanago were newly found in a few localities in the last decade (‘newly found’ locality or stock hereafter) (Fig.
Pseudorhodeus tanago from 13 known localities collected from 1993 through 2013 (55 individuals), four newly found stocks collected from 2010 through 2014 (22 individuals), and three individuals seized from an illegal home aquarium were materials of this research (Table
Past (a) and present (b) geographic distribution of Pseudorhodeus tanago. Dots indicate records of presence of the bitterling. Light blue painted areas contain one to three localities from which specimens of the present report came. Exact places of these localities are kept secret by authorities to protect from poaching. Numbers on map b stand for localities of samples listed in Table
# | Locality | Collecting date, year | Number of individuals |
---|---|---|---|
1 | Onjuku, Chiba | Dec. 19, 2012 | 4 |
2 | Katsuura 1, Chiba | Jun. 14, 2012 | 8 |
3 | Katsuura 2, Chiba | Jun. 14, 2012 | 4 |
4 | Katsuura 3, Chiba | Aug. 22, 2012 | 4 |
5 | Isumi A, Chiba | Jun. 21, 2012 | 4 |
6 | Isumi B, Chiba | Jun. 26, 2012 | 8 |
7 | Mobara, Chiba | Oct. 26, 2012 | 4 |
8 | Nagara, Chiba* | Jan. 28, 2013 | 4 |
9 | Namegawa, Saitama | 2003 | 5 |
10 | Small pond near Yokohama* | 2003 year class | 3 |
11 | South-east Tochigi | 2006 | 3 |
12 | North Tochigi | 1993-1994 | 1 |
13 | Handa natural habitat conservation area, Tochigi | 1993 | 3 |
14 | Central Chiba** | May 16, 2010 | 3 |
15 | North Chiba** | Jun. 3, 2010, Jul. 25, 2011 | 2 |
16 | Ibaraki A** | Nov. 6, 2014 | 16 |
17 | Ibaraki B** | 2014 | 1 |
18 | Seized from a home aquarium | 2014 | 3 |
Extraction of DNA was done from fin clips with QuickGene DNA Tissue kit on QuickGene-810 (Kurabo, Neyagawa, Japan). PCR primers were L14695 on the L-strand (AATTYTTGCTCRGACTCTAACC) and H15910 on the H-strand (GATCTTCGGATTACAAGACCGAT) which worked for amplifying a 1219 bp mitochondrial DNA fragment encompassing the whole cytochrome b gene and flanking tRNA partial sequences. PCR reaction mixture of 12.5 μL contained 1μL of template DNA, 0.96 μL of dNTP mix (2.5 nmol each), 1.2 μL of 10x ExTaq buffer, 0.06 μL (0.3 U) of ExTaq (Takara, Shiga, Japan), 1 μL of primers (5 pmol each), 7.28 μL of Milli-Q grade water. PCR reaction was of touchdown profile (
Indices of nucleotide divergence were calculated with Arlequin v.3.5 (
Sequencing Pseudorhodeus tanago mitochondrial cytochrome b revealed 10 haplotypes (Table
# \ Haplotype | Hap01 | Hap02 | Hap03 | Hap04 | Hap05 | Hap06 | Hap07 | Hap08 | Hap09 | Hap10 |
---|---|---|---|---|---|---|---|---|---|---|
1 | 4 | |||||||||
2 | 8 | |||||||||
3 | 4 | |||||||||
4 | 4 | |||||||||
5 | 2 | 2 | ||||||||
6 | 7 | 1 | ||||||||
7 | 4 | |||||||||
8* | 4 | |||||||||
9 | 5 | |||||||||
10* | 3 | |||||||||
11 | 3 | |||||||||
12 | 1 | |||||||||
13 | 3 | |||||||||
14** | 1 | 2 | ||||||||
15** | 1 | 1 | ||||||||
16** | 8 | 5 | 3 | |||||||
17** | 1 | |||||||||
18 | 1 | 2 |
Parsimonious haplotype network of mitochondrial DNA sequences of Pseudorhodeus tanago (Hap01–10). Each line connecting haplotypes indicates one nucleotide difference. Small open circles stand for missing haplotypes. Locality number (#) and number of individuals observed (parentheses) indicate haplotype sharing. Blue solid lines encircle haplotype groups (north, center, west and south). Blue broken line indicates sub-grouping within the south haplotype group exhibited by SAMOVA analysis (right, south_1; left, south_2).
Haplotypes and variable sites of Pseudorhodeus tanago mitochondrial cytochrome b.
Nucleotide site * | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
4 | 4 | 4 | 4 | 4 | 4 | 4 | 4 | 4 | 4 | 4 | 5 | 5 | 5 | 5 | 5 | 5 | 5 | 5 | |
4 | 4 | 5 | 5 | 5 | 5 | 8 | 9 | 9 | 9 | 9 | 0 | 1 | 2 | 3 | 4 | 4 | 4 | 5 | |
1 | 7 | 2 | 2 | 2 | 8 | 0 | 3 | 5 | 6 | 9 | 3 | 0 | 2 | 0 | 1 | 3 | 4 | 1 | |
2 | 2 | 0 | 6 | 9 | 0 | 7 | 7 | 5 | 7 | 4 | 0 | 0 | 8 | 1 | 5 | 2 | 7 | 3 | |
Haplotype | ** | ** | |||||||||||||||||
Hap01 | T | G | A | C | G | G | G | C | G | G | G | T | A | A | T | G | C | A | G |
Hap02 | T | G | A | C | G | G | A | C | G | G | G | C | A | A | T | G | C | A | G |
Hap03 | T | G | A | C | G | G | A | C | G | G | G | C | A | G | T | G | C | A | G |
Hap04 | T | G | A | C | G | G | A | C | G | G | G | C | A | A | T | G | T | A | G |
Hap05 | T | G | A | C | G | G | A | C | G | G | G | C | A | A | T | A | T | A | G |
Hap06 | T | A | A | T | G | G | A | C | G | C | A | C | A | A | T | G | T | A | G |
Hap07 | T | A | A | T | G | G | A | C | G | G | A | C | A | A | C | G | T | A | G |
Hap08 | T | A | A | T | G | A | A | T | G | G | A | C | A | A | T | G | T | G | G |
Hap09 | C | A | G | T | A | G | A | C | A | A | A | C | G | A | T | G | T | A | A |
Hap10 | C | A | G | T | A | G | A | C | A | A | A | C | G | A | T | G | T | G | A |
A haplotype network coincides with geographically structured population genetic characteristics of this species (Fig.
Geographic grouping of localities by AMOVA (SAMOVA) under assumptions of four or five geographic groups was similar to the haplotype grouping based on the haplotype network without geographic information (Table
Geographic grouping of known localities was largely coincident with river connectivity (Fig.
On the other hand, placement of specimens from newly found localities onto the haplotype network obscured the geographic population genetic structure, though SAMOVA analysis apportioned these localities to either of the South or Center geographic groups according to their haplotype composition. Coverage of variation among geographic groups reduced from 74.43–80.5% to 72.76–76.18% (Table
Contribution to variation in different patterns of population/stock grouping.
Grouping | Number of groups | Among groups | Among localities within groups | Within localities | Variance among groups | |||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
Sum of squares | % of variation | Fct | P(Fct) | Sum of squares | % of variation | Fst | P(Fst) | Sum of squares | % of variation | Fsc | P(Fsc) | |||
Known localities only | ||||||||||||||
North + center + west + south_1 (#1 + #3 + #6) + south_2 (#2 + #4 + #5) | 5 | 110,494 | 80,50 | 0,805 | low* | 15,562 | 13,35 | 0,939 | low | 8,125 | 6,15 | 0,685 | low | 2,532 |
North + {center + west} + south_1 (#1 + #3 + #6) + south_2 (#2 + #4 + #5) | 4 | 101,475 | 74,43 | 0,744 | low | 24,582 | 19,52 | 0,940 | low | 8,125 | 6,04 | 0,764 | low | 2,382 |
All localities (known + newly found) | ||||||||||||||
North + {center + #14} + west + south_3 (#2) + south4 (others) | 5 | 121,061 | 76,18 | 0,762 | low | 22,598 | 11,29 | 0,875 | low | 23,562 | 12,52 | 0,474 | low | 2,389 |
North + {center + #14} + west + south | 4 | 105,356 | 72,76 | 0,728 | low | 38,302 | 16,08 | 0,888 | low | 23,562 | 11,16 | 0,590 | low | 2,560 |
Genetic (πxy, below diagonal) and geographic (km, above diagonal) distance among collecting localities.
# | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | 13 | 14 | 15 | 16 | 17 |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
1 | 3.45 | 4.41 | 5.10 | 15.18 | 6.54 | 32.34 | 30.18 | 131.78 | 77.30 | 137.40 | 186.42 | 192.10 | 44.87 | 61.76 | 100.22 | 97.24 | |
2 | 4.000 | 1.00 | 2.41 | 15.54 | 9.22 | 32.37 | 29.88 | 129.87 | 74.61 | 137.21 | 186.22 | 192.17 | 44.35 | 61.48 | 101.46 | 96.42 | |
3 | 0.000 | 4.000 | 1.79 | 16.07 | 10.19 | 32.75 | 30.17 | 129.62 | 74.09 | 137.49 | 186.49 | 192.51 | 44.57 | 61.75 | 102.10 | 96.52 | |
4 | 2.000 | 2.000 | 2.000 | 17.85 | 11.39 | 34.54 | 31.94 | 130.99 | 75.05 | 139.25 | 188.25 | 194.29 | 46.31 | 63.51 | 103.84 | 98.20 | |
5 | 2.333 | 1.333 | 2.333 | 0.333 | 11.27 | 17.19 | 15.30 | 119.45 | 68.91 | 122.25 | 171.28 | 176.92 | 30.06 | 46.69 | 86.17 | 82.61 | |
6 | 0.000 | 3.250 | 0.000 | 1.500 | 1.708 | 28.08 | 26.52 | 130.38 | 78.05 | 132.94 | 181.94 | 187.25 | 41.26 | 57.63 | 94.19 | 93.82 | |
7 | 7.000 | 5.000 | 7.000 | 5.000 | 4.333 | 5.250 | 3.81 | 105.55 | 61.36 | 105.06 | 154.09 | 159.82 | 13.64 | 29.55 | 71.49 | 65.94 | |
8 | 7.000 | 5.000 | 7.000 | 5.000 | 4.333 | 5.250 | 0.000 | 105.40 | 59.23 | 107.34 | 156.35 | 162.36 | 14.76 | 31.62 | 75.08 | 67.31 | |
9 | 7.000 | 5.000 | 7.000 | 5.000 | 4.333 | 5.500 | 2.000 | 2.000 | 62.86 | 70.58 | 98.52 | 112.62 | 92.31 | 82.41 | 109.34 | 51.58 | |
10* | 9.000 | 7.000 | 9.000 | 7.000 | 6.333 | 7.500 | 4.000 | 4.000 | 4.000 | 103.40 | 146.69 | 140.01 | 52.65 | 55.93 | 106.07 | 59.01 | |
11 | 7.000 | 5.000 | 7.000 | 5.000 | 4.333 | 5.250 | 0.000 | 0.000 | 2.000 | 4.000 | 49.02 | 56.15 | 93.02 | 75.74 | 61.16 | 45.70 | |
12 | 13.000 | 11.000 | 13.000 | 11.000 | 10.333 | 11.375 | 7.000 | 7.000 | 8.000 | 10.000 | 7.000 | 15.79 | 141.99 | 124.74 | 103.09 | 92.74 | |
13 | 14.000 | 12.000 | 14.000 | 12.000 | 11.333 | 12.375 | 8.000 | 8.000 | 9.000 | 9.000 | 8.000 | 1.000 | 148.38 | 130.99 | 103.10 | 101.6 | |
14** | 3.667 | 2.333 | 3.667 | 1.667 | 1.000 | 2.333 | 0.000 | 0.000 | 1.333 | 3.333 | 0.000 | 6.667 | 7.667 | 17.40 | 65.84 | 52.56 | |
15** | 0.000 | 2.000 | 0.000 | 0.000 | 0.333 | -0.250 | 5.000 | 5.000 | 5.000 | 7.000 | 5.000 | 11.000 | 12.000 | 1.667 | 52.55 | 37.24 | |
16** | 3.000 | 1.000 | 3.000 | 1.000 | 0.333 | 2.250 | 4.000 | 4.000 | 4.000 | 6.000 | 4.000 | 10.000 | 11.000 | 1.333 | 1.000 | 58.96 | |
17** | 0.492 | 2.117 | 0.492 | 0.492 | 0.638 | 0.195 | 5.117 | 5.117 | 5.117 | 7.117 | 5.117 | 11.117 | 12.117 | 1.908 | -0.508 | 1.117 |
Geographic grouping of localities (encircled by gray solid lines) and river connectivity in 16th century (blue solid lines). Gray broken line encircles localities with the South geographic group including newly found localities (#15–17). Another newly found locality (#14) could be a member of the Center geographic group, but its haplotype composition is not typical among localities of the group. Arrows a and b denote historical overflow or stream capture terrains between river basins. Blue broken line indicates ancient river flow, and a blue arrow denotes change of the river flow at 500 ka.
The low sequence diversity of Pseudorhodeus tanago represented as monotypy in most known localities (11/13) conformed to their habitat characteristics (Table
Isolation on a wider and geological scale is also responsible for divergence among geographic groups. The North geographic group contained haplotypes distant from others (Hap01 and Hap10), and ancestors of the North and other geographic groups thus diverged first. The central part of the entire geographic range of P. tanago encircled the flat plain that underwent repeated marine transgressions (Fig.
Polymorphism and haplotype composition among newly found stocks, on the other hand, posed a question about the characteristics of genetic and geographic population structure of P. tanago. If we assume historical basis of distribution and genetic characteristics of these newly found stocks, it must be a series of unlikely geological events and haplotype sorting.
Including these newly found stocks, P. tanago as a whole showed a bipartite population structure in which some localities show genetic variation while others rarely do so (Table
There is no information such as confession or witness of violators about unauthorized releasing of P. tanago into the newly found localities. Genetic polymorphism of newly found localities may arrange stocks of these localities as conservation target of higher priority, if habitats are under risks of extinction. Recent decline of habitats indicates mitigation activities have not worked well (
Natural distribution in these newly found localities is, however, questionable because of unlikely association between haplotype/population relationships and river connectivity. Because bitterling fishes have attractive coloration to humans, there have been reports on unauthorized intentional releasing activities possibly implicated in home aquaria (
We dedicate this paper to Toshihiro Ishinabe (deceased) (Kannonzaki Nature Museum) and his contributions to conservation of Pseudorhodeus tanago over two decades carrying out habitat monitoring and ex situ breeding. KS, NS and KI conceived this research. KS analyzed data and drafted the manuscript. NS, MO, KI, TS, TM, TT and MT conducted field survey. All authors read, improved and approved the manuscript. This research was financially supported by a Protective Breeding Program of Ministry of the Environment, Japan. Permission for research use of P. tanago was provided by Agency for Cultural Affairs and Ministry of the Environment, Japan to NS and KI. We have no financial conflict of interests about contents of this publication. Specimens were collected under a condition in which the authorities expect non-disclosure of exact places of habitats except for a protected and publicly opened locality (#13). NS and KI are employees of two of these authorities that do not disclose some localities. Toshihiro Ishinabe, Hitoshi Kubota (Tochigi Prefectural Government), Kouki Kano (Ibaraki University), Masaki Miya (Natural History Museum and Institute, Chiba), Naoyuki Suguro (Kanagawa Prefectural Fisheries Technology Center) and Kazuhiro Umezawa (Saitama Prefectural Government) helped us collecting specimens. Maya Suzuki (Kanto Regional Environment Office, Ministry of the Environment) helped and coordinated field survey and sample collection. Takahiro Ohtsuka (National Research Institute of Fisheries Science) helped us for laboratory works.