The origins of limnetic forms and cryptic divergence in Gnathopogon fishes (Cyprinidae) in Japan

The cyprinid species of the genus Gnathopogon, exhibiting flexible morphological and ecological variation, include limnetic life forms. We examined the origin of the limnetic forms and the population divergence of the Japanese Gnathopogon species, using molecular phylogenetic and phylogeographic analyses. A Bayesian phylogenetic inference approach based on mtDNA cytochrome b sequence data revealed three major lineages in G. elongatus. One of them formed a monophyletic group with the limnetic species G. caerulescens, which is endemic to an ancient lake, Lake Biwa. The divergence of the G. caerulescens lineage was estimated to date back to the early Pleistocene. This precedes the formation of the extensive pelagic environment in the present Lake Biwa. However, the recent genetic divergence of G. caerulescens was inferred to originate in the present Lake Biwa in the late Pleistocene. Another lacustrine population in the Mikata Lakes was shown to belong to a different lineage from G. caerulescens. The majority of the population possessed unique, but non-monophyletic, haplotypes, suggesting a short evolutionary history. One of the cryptic lineages of G. elongatus discovered in the Ina Valley, the lower area of Lake Suwa, might be related to the extinct lacustrine subspecies G. elongatus suwae, which has been replaced by introduced congeners. The previous and ongoing introductions of Gnathopogon fishes would have produced genetic disturbance to the indigenous populations.     

Annual variations in length–weight relationship of Gymnogobius isaza (Tanaka, 1916) from Lake Biwa,Japan

The study describes length–weight relationships (LWRs) of Gymnogobius isaza (Tanaka, 1916) collected from Lake Biwa, Japan, between 1975 and 2002. Annual values of the allometric exponent of LWRs ranged from 2.47 to 3.42. This study is one of the few to report annual variations in LWR.     

Trophic position estimates of formalin-fixed samples with nitrogen isotopic compositions of amino acids: an application to gobiid fish (Isaza) in Lake Biwa, Japan

To evaluate the influence of formalin-fixation on the nitrogen isotopic composition of individual amino acids, we conducted a formalin-fixation experiment over 62 weeks using muscle tissues of largemouth bass (Micropterus salmoides) and freshwater minnow (Zacco platypus) collected in Lake Biwa. The results indicated that the formalin-fixation does not affect the nitrogen isotopic composition of amino acids including glutamic acid and phenylalanine from muscle protein of the fish. We applied this technique to formalin-fixed gobiid fish “Isaza” (Gymnogobius isaza, Tanaka) specimens that had been collected from Lake Biwa through the twentieth century. The pattern of the δ¹⁵N values for these amino acids remains relatively constant, even though each amino acid exhibits a gradual increase by ~3 ‰ from 1916 to 1992. With a formula proposed by Chikaraishi et al. (Limnol Oceanogr Method 7:740–750, 2009), we estimated the temporal variations in the trophic position of Isaza in Lake Biwa. The trophic position of Isaza remained quite constant (3.2–3.3) across a major eutrophication period in 1960–1980. The estimated trophic position of Isaza fish is consistent with the stomach content analysis in the previous reports. The constant trophic position suggests that the eutrophication did not apparently affect the trophic position of Isaza fish, although Isaza’s dietary preference changed from zooplankton to gammarids around 1970.     

Multiple colonizations of Lake Biwa by Sarcocheilichthys fishes and their population history

Rapid adaptation and speciation have not been well documented for organisms in Lake Biwa, Japan, the oldest ancient lake in East Asia. To examine these processes, we focused on the divergence of Sarcocheilichthys gudgeons and conducted phylogeographic and population genetic analyses using mitochondrial DNA sequences and microsatellite data. No remarkable genetic differentiation was found between two species, S. biwaensis and S. variegatus microoculus, both endemic to Lake Biwa, whereas this species group, including S. v. variegatus, was genetically divided into three major regional groups: the Lake Biwa, Ise Bay, and western groups. Divergent mtDNA haplotypes were included within the Lake Biwa and western groups, strongly suggesting secondary contacts among allopatrically isolated populations. Dating for mtDNA phylogeny using a geological constraint suggested the colonization of Lake Biwa by multiple Sarcocheilichthys lineages that diverged from each other in the early–middle Pleistocene. Coalescent-based population analyses indicated that the local populations colonized the rocky bottom habitat in Lake Biwa from other habitats after the Last Glacial Maximum, likely reflecting past environmental changes in the lake, including the disappearance of rocky areas during the glacial climate. Divergent morphological adaptation in Sarcocheilichthys associated with substrate type may have rapidly proceeded along with environmental changes.     

Range expansion and lineage admixture of the Japanese evergreen tree Machilus thunbergii in central Japan

We investigated the range expansion histories of Machilus thunbergii populations in the Kinki region of central Japan on the basis of nuclear microsatellite data. In the Kinki region, M. thunbergii is typically found in the coastal area, with some fragmented populations inland, around Lake Biwa. Phylogenetic and Bayesian clustering analysis (STRUCTURE analysis) revealed that the inland populations have different genetic components between the west and east sides of Lake Biwa. The population located on the north side of the lake has an admixture of the two genetically differentiated lineages, contributing to an increase in the genetic diversity of the population. Populations around Lake Biwa had lost rare alleles and the F value obtained from STRUCTURE analysis was lower in the coastal populations than in the lake populations. These results suggest that populations around Lake Biwa experienced a bottleneck due to a founder effect during the initial migration to the lake and that glacial refugia of M. thunbergii in the Kinki region existed along the coast.     

Long-Term Changes in the Diet of Gymnogobius isaza from Lake Biwa,Japan: Effects of Body Size and Environmental Prey Availability

Body size and environmental prey availability are both key factors determining feeding habits of gape-limited fish predators. However, our understanding of their interactive or relative effects is still limited. In this study, we performed quantitative dietary analysis of different body sizes of goby (Gymnogobius isaza) specimens collected from Lake Biwa between 1962 and 2004. First, we report that the diet was composed mainly of zooplankton (cladocerans and copepods) before the 1980s, and thereafter, shifted to zoobenthos (gammarids). This foraging shift coincided with, and thus can be linked to, known historical events in the lake at that time: decrease in zooplankton abundance with the alleviation of eutrophication, increase in fish body size resulting from fish population collapse, and increase in gammarid abundance due to reduced fish predation pressure. Supporting this view, our data analyses revealed how the long-term changes in the diet composition would be co-mediated by changes in fish body size and environmental prey availability. Specifically, while zoobenthos abundance strongly affected the fish diet composition, larger (smaller) fish preferred zoobenthos (zooplankton). Furthermore, the body size effects were stronger than those of prey availability. These results provide the best long-term evidence that fish feeding habits vary over decades with its body size and prey community due to anthropogenic disturbances.     

Is the Relationship between Body Size and Trophic Niche Position Time-Invariant in a Predatory Fish? First Stable Isotope Evidence

Characterizing relationships between individual body size and trophic niche position is essential for understanding how population and food-web dynamics are mediated by size-dependent trophic interactions. However, whether (and how) intraspecific size-trophic relationships (i.e., trophic ontogeny pattern at the population level) vary with time remains poorly understood. Using archival specimens of a freshwater predatory fish Gymnogobius isaza (Tanaka 1916) from Lake Biwa, Japan, we assembled a long-term (>40 years) time-series of the size-dependence of trophic niche position by examining nitrogen stable isotope ratios (δ ¹⁵N) of the fish specimens. The size-dependence of trophic niche position was defined as the slope of the relationship between δ ¹⁵N and log body size. Our analyses showed that the slope was significantly positive in about 60% of years and null in other years, changing through time. This is the first quantitative (i.e., stable isotope) evidence of long-term variability in the size-trophic relationship in a predatory fish. This finding had implications for the fish trophic dynamics, despite that about 60% of the yearly values were not statistically different from the long-term average. We proposed hypotheses for the underlying mechanism of the time-varying size-trophic relationship.     

Recent colonization by a coastal plant of inland habitats at an ancient freshwater lake,Lake Biwa: multilocus sequencing and a demographic history of Lathyrus japonicus (Fabaceae)

Ancient lakes have been recognized as “long‐term isolated islands” in terrestrial ecosystems. Lake Biwa, one of the few ancient lakes that formed around 4 million years ago, harbors many coastal species that commonly inhabit seashores. The beach pea, Lathyrus japonicus, is a typical coastal species of this freshwater lake, where morphological, physiological, and genetic differentiations have been reported between Biwa and coastal populations. Whether Biwa populations were isolated for long periods throughout Pleistocene climatic oscillations and subsequent range shifts is unclear. We assessed population genetic structure and demography of beach pea in this ancient freshwater lake using the sequences of eight nuclear loci. The results of STRUCTURE analyses showed evidence of admixture between Biwa and coastal populations, reflecting recent gene flow. The estimated demographic parameters implemented by the isolation with migration model (IM model) revealed a recent divergence (postglacial period) of Biwa populations, with some gene flow from Biwa to coastal populations. In addition, Biwa populations were significantly smaller in size than the ancestral or coastal populations. Our study suggests that a Holocene thermal maximum, when transgression could allow seeds from coastal plants to access Lake Biwa, was involved in the origin of the Biwa populations and their genetic divergence. Thus, coastal populations might have migrated to Lake Biwa relatively recently. Our study concluded that ancestral migrants in Lake Biwa were derived from small founding populations and accelerated genetic isolation of Biwa populations during short‐term isolation.     

Diversification and Demography of the Oriental Garden Lizard (Calotes versicolor) on Hainan Island and the Adjacent Mainland

The Oriental garden lizard (Calotes versicolor) is one of the few non-gekkonid lizards that are geographically widespread in the tropics. We investigated its population dynamics on Hainan Island and the adjacent mainland of China and Vietnam, focusing on the impact of cyclic upheaval and submergence of land bridges during the Pleistocene. Our Bayesian phylogenetic analysis reveals two mitochondrial lineages, A and B, which are estimated to have coalesced about 0.26 million years ago (95% credibility interval: 0.05–0.61 million years ago). Lineage A contains individuals mainly from central and southern Wuzhi Mountain on Hainan Island, whereas lineage B mainly comprises individuals from other sites on the island plus the adjacent mainland. The estimated coalescence times within lineages A (0.05 million years ago) and B (0.13 million years ago) fall within a period of cyclical land-bridge formation and disappearance in the Pleistocene. A spatial analysis of molecular variance identified two distinct population groupings: I, primarily containing lineage A, and II, mainly consisting of lineage B. However, haplotypes from lineages A and B occur sympatrically, suggesting that gene flow is ongoing. Neither Wuzhi Mountain nor Qiongzhou Strait and Gulf of Tonkin act as barriers to gene flow among C. versicolor populations. Analyses of the data using mismatch distributions and extended Bayesian skyline plots provide evidence of a relatively stable population size through time for Group I, and moderate population expansions and contractions during the end of the Pleistocene for Group II. We conclude that the phylogeographical patterns of C. versicolor are the combined product of Pleistocene sea-level oscillations and nonphysical barriers to gene flow.     

Long‐term isolation of the coastal plant Calystegia soldanella (Convolvulaceae) in ancient freshwater Lake Biwa,Japan

Lake Biwa is an ancient freshwater lake that was formed approximately 4 Mya and harbours many coastal plants that commonly inhabit the seashore. We used chloroplast DNA haplotype analysis using two spacer sequences and simple sequence repeat (SSR) analysis using eight nuclear microsatellite markers to detect genomic signatures indicating long‐term isolation of inland populations of Calystegia soldanella in Lake Biwa from coastal populations. We used 348 samples from 63 populations for haplotype analysis and 478 samples from 27 populations for SSR analysis covering the inland and coastal distribution of the species. We detected seven haplotypes, and the distribution pattern of these haplotypes was geographically highly structured between Lake Biwa and the coast. Nuclear SSR analysis also supported genetic differentiation between Lake Biwa and coastal populations (analyses of molecular variance, 43%), and the grouping of Lake Biwa and coastal populations by a Neighbour‐joining tree. In addition, genetic diversity of the inland populations (mean H_E = 0.153) was significantly lower than that of coastal populations (mean H_E = 0.328). These results suggested that inland populations at Lake Biwa have been isolated from coastal populations for a very long time. The inland populations most likely experienced a bottleneck effect, resulting in sufficient in situ genetic divergence to clearly distinguish them from coastal populations. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2011, 102 , 51–66.     

Introgression of Oncorhynchus masou subsp. (Biwa salmon) genome into lake-run O. m. ishikawae (Amago salmon) introduced into Lake Biwa, Japan

An investigation to determine the extent of introgression between Biwa salmon Oncorhynchus masou subsp. and lake-run Amago salmon O. m. ishikawae from Lake Biwa was conducted using mtDNA sequence and amplified fragment length polymorphism (AFLP) analysis. Hatchery-reared Amago salmon from the Samegai Trout Farm, the main source of stock released into the inflowing rivers of Lake Biwa, were also investigated. The extent of mtDNA introgression between these two subspecies was quite low: an anomalous haplotype was detected in just 1 of 55 Biwa salmon and 2 of 63 lake-run Amago salmon. Many of the mtDNA haplotypes found in lake-run Amago salmon were common to those in hatchery-reared Amago salmon. In the AFLP analysis, the nuclear genome of Amago salmon was hardly detected in Biwa salmon, but the Biwa salmon genome comprised 13.4% of the total amplified fragments from lake-run Amago salmon, being detected in different proportions in 28 of the 63 investigated individuals of the latter subspecies. These results suggest that introgression has occurred between Biwa salmon and lake-run Amago salmon descending from hatchery-reared Amago salmon from the Samegai Trout Farm.     

Phylogeny and historical demography of endemic fishes in Lake Biwa: the ancient lake as a promoter of evolution and diversification of freshwater fishes in western Japan

To elucidate the origins of the endemic fish of Lake Biwa, an ancient lake in Japan, and the role of the lake in the diversification of freshwater fish in western Japan, we established a molecular phylogenetic framework with an absolute time scale and inferred the historical demography of a large set of fish species in and around the lake. We used mtDNA sequences obtained from a total of 190 specimens, including 11 endemic species of Lake Biwa and their related species, for phylogenetic analyses with divergence time estimations and from a total of 2319 specimens of 42 species (including 14 endemics) occurring in the lake for population genetic analyses. Phylogenetic analysis suggested that some of the endemic species diverged from their closest relatives earlier (1.3–13.0 Ma) than the period in which the present environmental characteristics of the lake started to develop (ca. 0.4 Ma), whereas others diverged more recently (after 0.4 Ma). In contrast, historical demographic parameters suggested that almost all species, including endemic and nonendemic ones, expanded their populations after the development of the present lake environment. In phylogeographic analyses, common or very close haplotypes of some species were obtained from Lake Biwa and other regions of western Japan. The phylogenetic and historical demographic evidence suggests that there was a time lag between phylogenetic divergence and population establishment and that phenotypic adaptation of some endemic species to the limnetic environment occurred much later than the divergences of those endemic lineages. Population structure and phylogeographic patterns suggest that Lake Biwa has functioned not only as the center of adaptive evolution but also as a reservoir for fish diversity in western Japan.     

Intraspecific structure of sable <Emphasis Type="Italic">Martes zibellina</Emphasis> L. Inferred from nucleotide variation of the mitochondrial DNA cytochrome <Emphasis Type="Italic">b</Emphasis> gene

A fragment of the mitochondrial DNA (mtDNA) cytochrome b gene was sequenced in sable from Magadan oblast, Khabarovsk krai, and Kamchatka. Using phylogenetic analysis, the presence of two clusters (A and BC), with the divergence value of 1.4%, was demonstrated. Analysis of the cytochrome b gene median networks indicated that split of the ancestral population took place in early Pleistocene (about one Myr ago), while expansion of its more young phylogenetic group A occurred in late Pleistocene, about 120000 years ago.     

Juvenile migration of the exclusively pelagic cyprinid,Gnathopogon caerulescens (Honmoroko) in Lake Biwa,Central Japan

Migration of wild and cultivated juvenile honmoroko Gnathopogon caerulescens of from the spawning and nursery areas in Lake Biwa were investigated, both in the Ibanaiko Lagoon and its outlet to Daido River, using beam‐trawl surveys in 2013 and 2014. The study demonstrated migration of G. caerulescens from a nursery lagoon toward Lake Biwa after the juvenile stage. These findings appear to be the first direct evidence for migration of an exclusively pelagic cyprinid species from a littoral nursery to a pelagic adult habitat in a large deep lake.     

Disturbance of the indigenous gene pool of the threatened brook lamprey Lethenteron sp. S by intraspecific introgression and habitat fragmentation

Artificial introduction and habitat fragmentation affect the indigenous gene pools of fluvial animals. To investigate the effect of human activities on the genetic population structure of vulnerable brook lamprey Lethenteron sp. S in a single river system, samples from 12 tributaries of the Jinzu River, Japan, were analyzed using mitochondrial and microsatellite DNA markers. Exogenous Lake Biwa (Japan) haplotypes and alleles were detected in lampreys from several Jinzu River tributaries. Since Lake Biwa is the source of the commercial ayu fish Plecoglossus altivelis that is introduced in the Jinzu River, the exogenous Lake Biwa lamprey genotypes in the Jinzu River probably originated from the Lake Biwa lampreys that were unintentionally introduced along with the ayu fish. Bayesian admixture and mitochondrial DNA analyses revealed various genetic disturbance phases of the exogenous genotypes in the Jinzu River, such as the six indigenous populations, four admixed populations with low frequencies (average admixture proportion = 0.02–0.04; exogenous haplotype proportion <0.01), one introgressed population (0.71 and 0.57) and one population almost displaced by exogenous genotypes (0.93 and 0.96). Samples from three tributaries with weirs were genetically differentiated from the others by using pairwise F _ST and Bayesian analyses; the results suggested isolation by the weirs. Reduced mitochondrial DNA diversity was observed in 1 of the 3 samples probably due to reduced population size. These findings indicate that the indigenous lamprey populations in the Jinzu River are seriously affected by introgression with exogenous genotypes via unintentional introduction and habitat fragmentation by weirs.     

A re‐examination of the molecular phylogeny and biogeography of the genus Schizothorax (Teleostei: Cyprinidae) through enhanced sampling,with emphasis on the species in the Yunnan–Guizhou Plateau,China

The mitochondrial DNA cytochrome b sequences of 36 Schizothorax species from 51 localities in the Yunnan–Guizhou Plateau (YGP) and its adjacent areas were analysed. Maximum parsimony, Maximum likelihood and Bayesian phylogenetic analyses were performed to examine the relationships of Schizothorax species. A hypothesis of the phylogenetic relationships of the species is given. A relaxed molecular clock based on Bayesian evolutionary analysis was used to tentatively calculate the divergence times of Schizothorax. Samples from the YGP were tentatively grouped into three geographically distributed clades: the Tsangpo‐Irrawaddy, the Mekong‐Salween and the Trans‐Jinsha River (including Jinsha, Red, Nanpan and Beipan Rivers). Calibration of the molecular clock revealed that two geological periods, the late Miocene about 10 million years before present (Myr BP) and the Pliocene (4.0 Myr BP), were important times in the vicariant speciation of Schizothorax. The phylogenetic history of the species is congruent with events caused by the uplift of the Tibetan Plateau and the YGP. The divergence of Schizothorax species in YGP began in the Pliocene. Our phylogenetic trees did not support the hypothesis that the paleo Jinsha River was drained through the Yangtze River‐Jianchuan Lake‐Erhai Lake to the Red River. Schizothorax in the Beipan River were derived from the Jinsha River.     

Phylogeny of salmonids (salmoniformes: Salmonidae) and its molecular dating: Analysis of mtDNA data

Phylogenetic relationships among 41 species of salmonid fish and some aspects of their diversification-time history were studied using the GenBank and original mtDNA data. The position of the root of the Salmonidae phylogenetic tree was uncertain. Among the possible variants, the most reasonable seems to be that in which thymallins are grouped into the same clade as coregonins and the lineage of salmonins occupied a basal position relative to this clade. The genera of Salmoninae formed two distinct clades, i.e., (Brachymystax, Hucho) and (Salmo, Parahucho, (Salvelinus, (Parasalmo, Oncorhynchus)). Furthermore, the genera Parasalmo and Oncorhynchus were reciprocally monophyletic. The congruence of Salmonidae phylogenetic trees obtained using different types of phylogenetic markers is discussed. According to Bayesian dating, ancestral lineages of salmonids and their sister esocoids diverged about 106 million years ago. Sometime after, probably 100–70 million years ago, the salmonid-specific whole genome duplication took place. The divergence of salmonid lineages on the genus level occurred much later, within the time interval of 42–20 million years ago. The main wave of the diversification of salmonids at the species level occurred during the last 12 million years. The possible effect of genome duplication on the Salmonidae diversification pattern is discussed.     

Sequence evolution of mitochondrial tRNA genes and deep-branch animal phylogenetics

Mitochondrial DNA sequences are often used to construct molecular phylogenetic trees among closely related animals. In order to examine the usefulness of mtDNA sequences for deep-branch phylogenetics, genes in previously reported mtDNA sequences were analyzed among several animals that diverged 20–600 million years ago. Unambiguous alignment was achieved for stem-forming regions of mitochondrial tRNA genes by virtue of their conservative secondary structures. Sequences derived from stem parts of the mitochondrial tRNA genes appeared to accumulate much variation linearly for a long period of time: nearly 100 Myr for transition differences and more than 350 Myr for transversion differences. This characteristic could be attributed, in part, to the structural variability of mitochondrial tRNAs, which have fewer restrictions on their tertiary structure than do nonmitochondrial tRNAs. The tRNA sequence data served to reconstruct a well-established phylogeny of the animals with 100% bootstrap probabilities by both maximum parsimony and neighbor joining methods. By contrast, mitochondrial protein genes coding for cytochrome b and cytochrome oxidase subunit I did not reconstruct the established phylogeny or did so only weakly, although a variety of fractions of the protein gene sequences were subjected to tree-building. This discouraging phylogenetic performance of mitochondrial protein genes, especially with respect to branches originating over 300 Myr ago, was not simply due to high randomness in the data. It may have been due to the relative susceptibility of the protein genes to natural selection as compared with the stem parts of mitochondrial tRNA genes. On the basis of these results, it is proposed that mitochondrial tRNA genes may be useful in resolving deep branches in animal phylogenies with divergences that occurred some hundreds of Myr ago. For this purpose, we designed a set of primers with which mtDNA fragments encompassing clustered tRNA genes were successfully amplified from various vertebrates by the polymerase chain reaction.Abbreviations AA stem amino acid-acceptor stem - AC stem anticodon stem - COI cytochrome oxidase subunit I - cytb cytochrome b - D stem dihydrouridine stem - MP maximum parsimony - mtDNA mitochondrial DNA - Myr million years - NJ neighbor joining - PCR polymerase chain reaction - Ti transition - T stem tC stem - Tv transversion Correspondence to: Y. Kumazawa     

Phylogenetic relationships of four species of floating gobies (Gymnogobius) as inferred from partial mitochondrial cytochrome b gene sequences

 To clarify the intra- and interspecific relationships of four Gymnogobius species, G. urotaenia, G. isaza, Gymnogobius sp. 1 (sumiukigori), and Gymnogobius sp. 2 (shimaukigori), partial mitochondrial cytochrome b gene sequences of 639 bp were obtained for a total of 31 specimens from Japan and Korea, plus 2 outgroup specimens. Twenty-nine haplotypes were identified in the ingroup, with a total of 122 variable sites (19.1%). The individuals regarded as the same species morphologically were monophyletic genetically. Sequence differences between amphidromous individuals of three species distributed in both Japan and Korea were relatively small (0.16–1.25%). The largest intraspecific sequence difference was observed between individuals of G. urotaenia from Lake Biwa and those from other localities (1.25–2.19%). Interspecific sequence differences ranged from 4.07% to 13.46%; neighbor-joining, maximum-parsimony, and maximum-likelihood methods indicated that Gymnogobius sp. 2 diverged first, followed by G. isaza, with Gymnogobius sp. 1 and G. urotaenia being monophyletic. The estimated divergence time of each species, based on estimated divergence rates for mitochondrial protein-coding genes already reported (0.8–2.8%/my), suggested that speciation occurred mainly in the Pliocene (possibly Miocene), with G. isaza (a Lake Biwa endemic) diverging significantly earlier (probably Lake Kouga stage) than estimated in previous studies. In contrast, according to the previous hypothesis, the substitution rates were highly overestimated to about 12–20%/my. Received: November 19, 2001 / Revised: May 20, 2002 / Accepted: June 18, 2002 Acknowledgments Special thanks are due to Katsutoshi Watanabe and Yuji Yamazaki for valuable advice and providing fish samples. We also thank Kouji Nakayama, Toshiyuki Ohkawa, Motoomi Yamaguchi, and members of the Laboratory of Marine Biology, Fukui Prefectural University, and members of the Laboratory of Molecular Marine Biology, Ocean Research Institute, University of Tokyo, and the Laboratory of Marine Stock Enhancement, University of Kyoto, for their help and advice during the present study. Correspondence to:Shigeo Harada     

PCR-RFLP analysis of mitochondrial DNA cytochrome b gene among Haruan (Channa striatus) in Malaysia

Haruan (Channa striatus) is in great demand in the Malaysian domestic fish market. In the present study, mtDNA cyt b was used to investigate genetic variation of C. striatus among populations in Peninsular Malaysia. The overall population of C. striatus demonstrated a high level of haplotype diversity (h) and a low-to-moderate level of nucleotide diversity (π). Analysis of molecular variance (AMOVA) results showed a significantly different genetic differentiation among 6 populations (F_ST = 0.37566, P = 0.01). Gene flow (Nm) was high and ranged from 0.32469 to infinity (∞). No significant relationship between genetic distance and geographic distance was detected. A UPGMA tree based on the distance matrix of net interpopulation nucleotide divergence (d_A) and haplotype network of mtDNA cyt b revealed that C. striatus is divided into 2 major clades. The neutrality and mismatch distribution tests for all populations suggested that C. striatus in the study areas had undergone population expansion. The estimated time of population expansion in the mtDNA cyt b of C. striatus populations occurred 0.72-6.19 million years ago. Genetic diversity of mtDNA cyt b and population structure among Haruan populations in Peninsular Malaysia will be useful in fisheries management for standardization for Good Agriculture Practices (GAP) in fish-farming technology, as well as providing the basis for Good Manufacturing Practices (GMP).