Ghost introgression in ricefishes of the genus Adrianichthys in an ancient Wallacean lake.

Because speciation might have been promoted by ancient introgression from an extinct lineage, it is important to detect the existence of 'ghost introgression' in focal taxa and examine its contribution to their diversification. In this study, we examined possible ghost introgression and its contributions to the diversification of ricefishes of the genus Adrianichthys in Lake Poso, an ancient lake on Sulawesi Island, in which some extinctions are known to have occurred. Population-genomic analysis revealed that two extant Adrianichthys species, A. oophorus and A. poptae are reproductively isolated from each other. Comparisons of demographic models demonstrated that introgression from a ghost population, which diverged from the common ancestor of A. oophorus and A. poptae, is essential for reconstructing the demographic history of Adrianichthys. The best model estimated that the divergence of the ghost population greatly predated the divergence between A. oophorus and A. poptae, and that the ghost population secondarily contacted the two extant species within Lake Poso more recently. Genome scans and simulations detected a greatly divergent locus, which cannot be explained without ghost introgression. This locus was also completely segregated between A. oophorus and A. poptae. These findings suggest that variants that came from a ghost population have contributed to the divergence between A. oophorus and A. poptae, but the large time-lag between their divergence and ghost introgression indicates that the contribution of introgression may be restricted.

A new sampling method with zirconium-loaded resin for phosphate oxygen isotope analysis in oligotrophic freshwater systems.

RATIONALE: The phosphate oxygen isotope ratio ( δ 18 O PO 4 ) is a useful technique to trace the sources and biogeochemical cycles of phosphorus (P) in aquatic ecosystems. However, δ 18 O PO 4 has not been widely used in oligotrophic freshwater systems due to technical and methodological difficulties in collecting sufficient phosphate (PO4 ) for the δ 18 O PO 4 analysis, which sometimes requires hundreds of liters of the water sample. In this study, a new approach (PaS-Zir) was developed for the δ 18 O PO 4 analysis in oligotrophic freshwater systems using zirconium (Zr)-loaded (ZrIRC) resin, which has a high affinity for PO4 . METHODS: ZrClO2 was added to Amberlite IRC748 to obtain the ZrIRC resin. The adsorption/desorption experiment using KH2 PO4 with a known value of δ 18 O PO 4 was conducted to determine the adsorption/desorption properties of the resin and the likelihood of isotopic fractionation. By installing mesh bags filled with the resin, the PaS-Zir approach was used in two rivers with low PO4 concentrations (0.2 and 5.3 µmol/L). A conventional sampling method was also performed in the study river with a higher PO4 concentration to validate the efficacy of the PaS-Zir method. RESULTS: The adsorption/desorption experiment demonstrated that the ZrIRC resin possessed a sufficient adsorption capacity (153 µmol/resin-mL) and exhibited little isotopic fractionation during the adsorption/desorption processes. Using the PaS-Zir method, we were able to collect sufficient PO4 samples for the δ 18 O PO 4 analysis from the rivers within at least 4 days of mesh bag installation. The δ 18 O PO 4 values (14.2‰ ± 0.2‰) obtained using the PaS-Zir method were comparable to those obtained using the conventional method (14.0‰ ± 0.03‰). CONCLUSION: We proved that the PaS-Zir method is applicable to oligotrophic freshwater systems and is generally more efficient than the conventional method. In addition, our method is useful for improving the understanding of the P dynamics of oligotrophic ecosystems because of the extremely low concentration of PO4 commonly found in them, which are often prone to P pollution.

A new method for phosphate purification for oxygen isotope ratio analysis in freshwater and soil extracts using solid-phase extraction with zirconium-loaded resin.

RATIONALE: Phosphate (PO4 ) oxygen isotope (δ18 OPO4 ) analysis is increasingly applied to elucidate phosphorus cycling. Due to its usefulness, analytical methods continue to be developed and improved to increase processing efficiency and applicability to various sample types. A new pretreatment procedure to obtain clean Ag3 PO4 using solid-phase extraction (SPE) with zirconium-loaded resin (ZrME), which can selectively adsorb PO4 , is presented and evaluated here. METHODS: Our method comprises (1) PO4 concentration, (2) PO4 separation by SPE, (3) cation removal, (4) Cl- removal, and (5) formation of Ag3 PO4 . The method was tested by comparing the resulting δ18 OPO4 of KH2 PO4 reagent, soil extracts (NaHCO3 , NaOH, and HCl), freshwater, and seawater with data obtained using a conventional pretreatment method. RESULTS: PO4 recovery of our method ranged from 79.2% to 97.8% for KH2 PO4 , soil extracts, and freshwater. Although the recovery rate indicated incomplete desorption of PO4 from the ZrME columns, our method produced high-purity Ag3 PO4 and accurate δ18 OPO4 values (i.e., consistent with those obtained using conventional pretreatment methods). However, for seawater, the PO4 recovery was low (1.1%), probably due to the high concentrations of F- and SO4 2- which interfere with PO4 adsorption on the columns. Experiments indicate that the ZrME columns could be regenerated and used repeatedly at least three times. CONCLUSIONS: We demonstrated the utility of ZrME for purification of PO4 from freshwater and soil extracts for δ18 OPO4 analysis. Multiple samples could be processed in three days using this method, increasing sample throughput and potentially facilitating more widespread use of δ18 OPO4 analysis to deepen our understanding of phosphorus cycling in natural environments.

Resource partitioning is not coupled with assortative mating in sympatrically divergent ricefish in a Wallacean ancient lake.

Sympatric speciation is considered to be difficult without the coupling between ecological traits that allow resource partitioning and reproductive traits that allow assortative mating. Such "magic traits" are known to be involved in most of the compelling examples of sympatric speciation. In this study, we report a possible case of sympatric speciation without magic traits. Three species of ricefish (genus Oryzias) are suggested to have diverged sympatrically within Lake Poso, an ancient lake in Sulawesi. An analysis of genome-wide single-nucleotide polymorphisms showed that these three species are reproductively isolated from each other throughout the lake. Stable isotope analyses revealed that the three species use different food resources, which reflect differences in their feeding morphologies (gill rakers and digestive tracts) and feeding sites. Field and laboratory observations showed that O. nebulosus and O. orthognathus share a mating habitat of cobbles, where they scatter fertilized eggs, whereas this site is never used by O. nigrimas, indicating that assortative mating is partly achieved by spatial isolation. The small, less-adhesive eggs of O. nebulosus and O. orthognathus probably reflect their adaptation to spawning on cobble beaches. Laboratory mating experiments showed strong prezygotic isolation between O. nebulosus and O. orthognathus, which is achieved by strong species recognition presumably by both sexes based on species-specific mating dances and nuptial coloration. In summary, the assortative mating of O. nebulosus and O. orthognathus is probably not coupled to resource partitioning. We discussed how sympatric speciation among these species might have been achieved even without magic traits.

Identification of Phosphorus Sources in a Watershed Using a Phosphate Oxygen Isoscape Approach.

Identifying nonpoint phosphorus (P) sources in a watershed is essential for addressing cultural eutrophication and for proposing best-management solutions. The oxygen isotope ratio of phosphate (δ18OPO4) can shed light on P sources and P cycling in ecosystems. This is the first assessment of the δ18OPO4 distribution in a whole catchment, namely, the Yasu River Watershed in Japan. The observed δ18OPO4 values in the river water varied spatially from 10.3‰ to 17.6‰. To identify P sources in the watershed, we used an isoscape approach involving a multiple-linear-regression model based on land use and lithological types. We constructed two isoscape models, one using data only from the whole watershed and the other using data from the small tributaries. The model results explain 69% and 96% of the spatial variation in the river water δ18OPO4. The lower R2 value for the whole watershed model is attributed to the relatively large travel time for P in the main stream of the lower catchment that can result in cumulative biological P recycling. Isoscape maps and a correlation analysis reveal the relative importance of P loading from paddy fields and bedrock. This work demonstrates the utility of δ18OPO4 isoscape models for assessing nonpoint P sources in watershed ecosystems.

Using food network unfolding to evaluate food-web complexity in terms of biodiversity: theory and applications.

Food-web complexity often hinders disentangling functionally relevant aspects of food-web structure and its relationships to biodiversity. Here, we present a theoretical framework to evaluate food-web complexity in terms of biodiversity. Food network unfolding is a theoretical method to transform a complex food web into a linear food chain based on ecosystem processes. Based on this method, we can define three biodiversity indices, horizontal diversity (DH ), vertical diversity (DV ) and range diversity (DR ), which are associated with the species diversity within each trophic level, diversity of trophic levels, and diversity in resource use, respectively. These indices are related to Shannon's diversity index (H'), where H' = DH  + DV  - DR . Application of the framework to three riverine macroinvertebrate communities revealed that D indices, calculated from biomass and stable isotope features, captured well the anthropogenic, seasonal, or other within-site changes in food-web structures that could not be captured with H' alone.

Differential Responses of Two Ecologically Similar Case-Bearing Caddisfly Species to a Fish Chemical Cue: Implications for a Coexistence Mechanism.

The mechanisms for the coexistence of multiple species occupying the same ecological niche are often puzzling. Predator effects on competitively superior species is one possible mechanism. In this study, we tested whether the presence of size-selective predators (fishes) acts as a mechanism for the coexistence of two species of case-bearing caddisfly larvae, Perissoneura paradoxa and Psilotreta kisoensis (Odontoceridae, Trichoptera). The larvae of these two species have similar ecological and life history traits except their body size, and they have been found to coexist only in habitats shared with predatory fishes. Experiments on intra and interspecific competition revealed that the larger Pe. paradoxa always outcompeted the smaller Ps. kisoensis in the absence of predatory fishes, suggesting that Pe. paradoxa performed intra-guild predation on Ps. kisoensis. We also conducted experiments to examine how strongly each of these species responded in terms of case repair with/without a predator chemical cue after their cases were partly dismantled. Perissoneura paradoxa exhibited a stronger case repair response in the presence of a predator chemical cue than that exhibited by Ps. kisoensis, suggesting that Pe. paradoxa is more vulnerable to fish predation, probably because their body size is in the preferred prey range of fishes. We suggest that the presence of predators works in the favor of smaller, subordinate species through size-selective predator effects, enabling these two competitive species to coexist in the same habitat.

Genomic architecture of habitat-related divergence and signature of directional selection in the body shapes of Gnathopogon fishes.

Evolution of ecomorphologically relevant traits such as body shapes is important to colonize and persist in a novel environment. Habitat-related adaptive divergence of these traits is therefore common among animals. We studied the genomic architecture of habitat-related divergence in the body shape of Gnathopogon fishes, a novel example of lake-stream ecomorphological divergence, and tested for the action of directional selection on body shape differentiation. Compared to stream-dwelling Gnathopogon elongatus, the sister species Gnathopogon caerulescens, exclusively inhabiting a large ancient lake, had an elongated body, increased proportion of the caudal region and small head, which would be advantageous in the limnetic environment. Using an F2 interspecific cross between the two Gnathopogon species (195 individuals), quantitative trait locus (QTL) analysis with geometric morphometric quantification of body shape and restriction-site associated DNA sequencing-derived markers (1622 loci) identified 26 significant QTLs associated with the interspecific differences of body shape-related traits. These QTLs had small to moderate effects, supporting polygenic inheritance of the body shape-related traits. Each QTL was mostly located on different genomic regions, while colocalized QTLs were detected for some ecomorphologically relevant traits that are proxy of body and caudal peduncle depths, suggesting different degree of modularity among traits. The directions of the body shape QTLs were mostly consistent with the interspecific difference, and QTL sign test suggested a genetic signature of directional selection in the body shape divergence. Thus, we successfully elucidated the genomic architecture underlying the adaptive changes of the quantitative and complex morphological trait in a novel system.

Stable nitrogen isotopic composition of amino acids reveals food web structure in stream ecosystems.

The stable N isotopic composition of individual amino acids (SIAA) has recently been used to estimate trophic positions (TPs) of animals in several simple food chain systems. However, it is unknown whether the SIAA is applicable to more complex food web analysis. In this study we measured the SIAA of stream macroinvertebrates, fishes, and their potential food sources (periphyton and leaf litter of terrestrial C3 plants) collected from upper and lower sites in two streams having contrasting riparian landscapes. The stable N isotope ratios of glutamic acid and phenylalanine confirmed that for primary producers (periphyton and C3 litter) the TP was 1, and for primary consumers (e.g., mayfly and caddisfly larvae) it was 2. We built a two-source mixing model to estimate the relative contributions of aquatic and terrestrial sources to secondary and higher consumers (e.g., stonefly larva and fishes) prior to the TP calculation. The estimated TPs (2.3-3.5) roughly corresponded to their omnivorous and carnivorous feeding habits, respectively. We found that the SIAA method offers substantial advantages over traditional bulk method for food web analysis because it defines the food web structure based on the metabolic pathway of amino groups, and can be used to estimate food web structure under conditions where the bulk method cannot be used. Our result provides evidence that the SIAA method is applicable to the analysis of complex food webs, where heterogeneous resources are mixed.

A RAD-based linkage map and comparative genomics in the gudgeons (genus Gnathopogon, Cyprinidae).

BACKGROUND: The construction of linkage maps is a first step in exploring the genetic basis for adaptive phenotypic divergence in closely related species by quantitative trait locus (QTL) analysis. Linkage maps are also useful for comparative genomics in non-model organisms. Advances in genomics technologies make it more feasible than ever to study the genetics of adaptation in natural populations. Restriction-site associated DNA (RAD) sequencing in next-generation sequencers facilitates the development of many genetic markers and genotyping. We aimed to construct a linkage map of the gudgeons of the genus Gnathopogon (Cyprinidae) for comparative genomics with the zebrafish Danio rerio (a member of the same family as gudgeons) and for the future QTL analysis of the genetic architecture underlying adaptive phenotypic evolution of Gnathopogon. RESULTS: We constructed the first genetic linkage map of Gnathopogon using a 198 F2 interspecific cross between two closely related species in Japan: river-dwelling Gnathopogon elongatus and lake-dwelling Gnathopogon caerulescens. Based on 1,622 RAD-tag markers, a linkage map spanning 1,390.9 cM with 25 linkage groups and an average marker interval of 0.87 cM was constructed. We also identified a region involving female-specific transmission ratio distortion (TRD). Synteny and collinearity were extensively conserved between Gnathopogon and zebrafish. CONCLUSIONS: The dense SNP-based linkage map presented here provides a basis for future QTL analysis. It will also be useful for transferring genomic information from a "traditional" model fish species, zebrafish, to screen candidate genes underlying ecologically important traits of the gudgeons.

Genotypic composition and the relationship between genotypic composition and geographical proximity of the cyanobacterium Microcystis aeruginosa in western Japan.

Microcystis aeruginosa is one of the bloom-forming harmful algae in freshwater ecosystems. We genetically characterized Microcystis populations during bloom-forming periods in various reservoirs, lakes, and ponds in Japan during 2009. Using phylogenetic analysis, we evaluated the relationship between current genotype expansions and geographic location within western Japan and intraspecific variation. Microcystis aeruginosa colonies were isolated at 15 sites and were analyzed by sequencing the 16S-23S internal transcribed spacer (ITS) region of the ribosomal operon, and the potential to produce toxins was assessed by PCR-based detection of the microcystin synthetase gene mcyG. In total, 171 colonies were separated into 41 genotypes. The highest genotypic composition was detected in the south basin of Lake Biwa and the lowest in Lagoon Iba. Cluster analysis indicated no obvious association between genotypic composition and geographic distance. Thus, clear genetic differentiation accompanied by geographic origins was not found in western Japan. The resulting neighbor-joining tree revealed 3 clusters, 2 of which contained strains that showed both nonamplification and amplification of the mcyG gene.

Development of nine markers and characterization of the microsatellite loci in the endangered Gymnogobius isaza (Gobiidae).

Gymnogobius isaza is a freshwater goby endemic to Lake Biwa, Japan. They experienced a drastic demographic bottleneck in the 1950s and 1980s and slightly recovered thereafter, but the population size is still very small. To reveal dynamics of genetic diversity of G. isaza, we developed nine microsatellite markers based on the sequence data of a related goby Chaenogobius annularis. Nine SSR (Simple Sequence Repeats) markers were successfully amplified for raw and formalin-fixed fish samples. The number of alleles and expected heterozygosities ranged from one to 10 and from 0.06 to 0.84, respectively, for the current samples, while one to 12 and 0.09 to 0.83 for historical samples. The markers described here will be useful for investigating the genetic diversity and gene flow and for conservation of G. isaza.

Community capability building for environmental conservation in Lake Biwa (Japan) through an adaptive and abductive approach.

In the south basin of Lake Biwa, Shiga, Japan, overgrown aquatic weeds (submerged macrophytes) impede cruising boats and cause unpleasant odors and undesirable waste when washed ashore. To address this socio-ecological problem, Shiga Prefectural Government implemented a public program to remove overgrown weeds and compost them ashore to conserve the lake environment, while coastal inhabitants and occasional volunteers remove weeds from the beaches to maintain the quality of the living environment. However, these effects are limited because of disjointed social networks. We applied an adaptive and abductive approach to develop community capability to jointly address this problem by sharing academic knowledge with local actors and empowering them. The initial multifaceted reviews, including interviews and postal questionnaire surveys, revealed that the agro-economic value of composted weeds declined in historical and socio-psychological contexts and that most of the unengaged public relied on local governments to address environmental problems. These findings were synthesized and assessed with workshop participants, including local inhabitants, governmental agents, businesspeople, social entrepreneurs, and research experts, to unearth the best solution. The workshops resulted in the development of an e-point system, called Biwa Point, to promote and acknowledge voluntary environmental conservation activities, including beach cleaning. It may contribute to enhancing the socio-ecological capability of communities. Additionally, ethical issues, such as publication of inconvenient truths, undesired interpretation by the researchers, and social constraints in research methods, arose through our research practice. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s42532-021-00078-3.

Remarkable spatial memory in a migratory cardinalfish.

The ability to orient and navigate within a certain environment is essential for all animals, and spatial memory enables animals to remember the locations of such markers as predators, home, and food. Here we report that the migratory marine cardinalfish Apogon notatus has the potential to retain long-term spatial memory comparable to that of other animals. Female A. notatus establish a small territory on a shallow boulder bottom to pair and spawn with males. We carried out field research in two consecutive breeding seasons on territory settlement by individually marked females. Females maintained a territory at the same site throughout one breeding season. After overwintering in deep water, many of them (82.1%) returned to their breeding ground next spring and most occupied the same site as in the previous season, with only a 0.56 m shift on average. Our results suggest that female A. notatus have long-distance homing ability to pinpoint the exact location of their previous territory, and retain spatial memory for as long as 6 months.

Analysis of biomagnification of persistent organic pollutants in the aquatic food web of the Mekong Delta, South Vietnam using stable carbon and nitrogen isotopes.

The present study elucidated the biomagnification profiles of persistent organic pollutants (POPs) through a tropical aquatic food web of Vietnam based on trophic characterization using stable nitrogen analysis. Various biological samples collected from the main stream of the Mekong Delta were provided for the analysis for both POPs, and stable nitrogen and carbon isotope ratios. Of the POPs analyzed, dichlorodiphenyltrichloroethane and its metabolites (DDTs) were the predominant contaminants with concentrations ranging from 0.058 to 12 ng/g wet weight, followed by polychlorinated biphenyls (PCBs) at 0.017-8.9 ng/g, chlordane compounds (CHLs) at 0.0043-0.76 ng/g, tris-4-chlorophenyl methane (TCPMe) at N.D.-0.26 ng/g, hexachlorocyclohexane isomers (HCHs) at N.D.-0.20 ng/g and hexachlorobenzene (HCB) at 0.0021-0.096 ng/g. Significant positive increases of concentrations in DDTs, CHLs, and TCPMe against the stable nitrogen ratio (delta(15)N) were detected, while, concentrations of HCHs and HCB showed no significant increase. The slopes of the regression equations between the log-transformed concentrations of these POPs and delta(15)N were used as indices of biomagnification. The slopes of the POPs for which positive biomagnification was detected ranged from 0.149 to 0.177 on a wet weight basis. The slopes of DDTs and CHLs were less than those reported for a marine food web of the Arctic Ocean, indicating that less biomagnification had occurred in the tropical food web. Of the isomers of CHLs, unlike the studies of the Arctic Ocean, oxychlordane did not undergo significant biomagnification through the food web of the Mekong Delta. This difference is considered to be due to a lack of marine mammals, which might metabolize cis- and trans-chlordane to oxychlordane, in the Mekong Delta ecosystem. The biomagnification profile of TCPMe is reported for the first time in the present study.

Biomagnification of trace elements in the aquatic food web in the Mekong Delta, South Vietnam using stable carbon and nitrogen isotope analysis.

In the present study, we report the concentrations of 21 trace elements (V, Cr, Mn, Co, Cu, Zn, As, Se, Rb, Sr, Mo, Ag, Cd, Sn, Sb, Cs, Ba, Hg, Tl, Pb, and Bi), as well as the results of the analysis of stable carbon and nitrogen isotopes, of the various biota that make up the food web in the main stream of the Mekong Delta near Can Tho, South Vietnam. A significant trophic level-dependent increase was found in concentrations of Se, Rb, and Hg with increasing delta(15)N, indicating that an overall biomagnification of these elements occurred. However, the increase of Hg concentration per trophic level was lower than in previous studies. In contrast, the concentration of Mn showed an opposite trend in the food web of the Mekong Delta. In addition to these overall trends, the present study revealed that the biomagnification profiles of trace metals differ between crustaceans and fishes; concentrations of Mn, Cu, Zn, As, Sr, Mo, Ag, Cd, Sb, Cs, Ba, Tl, and Pb were significantly higher in crustaceans, whereas fishes showed higher concentrations of Cr, Rb, and Hg (trophic level determined by delta(15)N). The differences in the biomagnification profile between the major taxa might be attributed to differences in metal accumulation and in detoxification abilities such as possessing a metal-binding protein, e.g., metallothionein (MT).

Effects of resource-dependent cannibalism on population size distribution and individual life history in a case-bearing caddisfly.

Resource availability often determines the intensity of cannibalism, which has a considerable effect on population size distribution and individual life history. Larvae of the caddisfly Psilotreta kisoensis build portable cases from sedimentary sands and often display cannibalism. For this species, the availability of preferable case material is a critical factor that affects larval fitness, and material is locally variable depending on the underlying geology. In this study, we investigated how sand quality as a case material determines cannibalism frequency among larvae and, in turn, how the differential cannibalism frequency affects the body-size distribution and voltinism. Rearing experiments within a cohort revealed that a bimodal size distribution developed regardless of material quality. However, as the preferable material became abundant, the proportion of larger to smaller individuals increased. Consecutive experiments suggested that smaller larvae were more frequently cannibalized by larger ones and excluded from the population when preferable smooth material was abundant. This frequent cannibalism resulted in a bimodal size distribution with a significantly higher proportion of larger compared to smaller individuals. The size-dependent cannibalism was significantly suppressed when the larvae were raised in an environment with a scarcity of the preferable case material. This is probably because larvae cannot enjoy the benefit of rapid growth by cannibalism due to the difficulties in enlarging their case. At low cannibalism the growth of smaller individuals was stunted, and this was probably due to risk of cannibalism by larger individuals. This growth reduction in small individuals led to a bimodal size-distribution but with a lower proportion of larger to smaller individuals compared to at high cannibalism. A field study in two streams showed a similar size distribution of larvae as was found in the rearing experiment. The bimodal ratio has consequences for life history, since a size-bimodal population causes a cohort splitting: only larvae that were fully grown at 1 year had a univoltine life cycle, whereas larvae with a stunted growth continued their larval life for another year (semivoltine). This study suggests that availability of preferable case building material is an important factor that affects cannibalism, which in turn affects larval population size structure and cohort splitting.

Correction: Effects of resource-dependent cannibalism on population size distribution and individual life history in a case-bearing caddisfly.

[This corrects the article DOI: 10.1371/journal.pone.0191925.].

Spatial genetic structure and body size divergence in endangered Gymnogobius isaza in ancient Lake Biwa.

Gymnogobius isaza is a freshwater goby endemic to ancient Lake Biwa, the largest lake in Japan. The species is now listed as 'Critically Endangered' in the Red Data Book of Japan. Nevertheless, it remains subject to fishing without any specific management strategies. Previous studies using mitochondrial DNA markers showed that this fish species has two cryptic lineages. However, little is known about spatial genetic structure and ecological differences across the broad lakescape. In this study, we collected fish samples at nine locations along the lakeshore during the breeding season and tested for the presence of spatial heterogeneity in the lineage's composition while measuring body size as the most fundamental biological trait. The results showed that the major lineage dominated all the sampling locations whereas the minor lineage consisted of only 11% (16/143) of samples. Furthermore, although their spatial distributions overlapped (i.e. the two lineages may be well mixed), we found it possible that the minor lineage may have a potentially narrower distribution than the major lineage. In addition, we found that the two lineages differ in body size; specifically, the minor lineage is smaller in size. From the viewpoint of genetic diversity conservation and sustainable resource use, this fish should be managed as two genetic stocks and spatial and/or body size-based fishery management is desirable, with particular attention to the minor (smaller sized) lineage.

Integrated trophic position decreases in more diverse communities of stream food webs.

The relationship between biodiversity and ecosystem functioning is an important theme in environmental sciences. We propose a new index for configuration of the biomass pyramid in an ecosystem, named integrated trophic position (iTP). The iTP is defined as a sum of trophic positions (i.e. the average number of steps involved in biomass transfer) of all the animals in a food web integrated by their individual biomass. The observed iTP for stream macroinvertebrates ranged from 2.39 to 2.79 and was negatively correlated with the species density and the Shannon-Wiener diversity index of the local community. The results indicate a lower efficiency of biomass transfer in more diverse communities, which may be explained by the variance in edibility hypothesis and/or the trophic omnivory hypothesis. We found a negative effect of biodiversity on ecosystem functioning.