Mitochondrial genome of the yellow catfish Pelteobagrus fulvidraco and insights into Bagridae phylogenetics

The mitochondrial genome (mitogenome) can provide important information for understanding phylogenetic analysis and molecular evolution. Herein, we amplified the complete mitogenome sequence of Pelteobagrus fulvidraco. The mitogenome was 16,526 bp in length and included 13 protein-coding genes (PCGs), 22 transfer RNA genes, two ribosomal RNA genes and a non-coding control region (D-loop). Both the organization and location of genes in the mitogenome were consistent with those from Siluriformes fishes previously published in GenBank. The phylogenetic relationships based on Bayesian inference (BI) and Maximum likelihood (ML) methods showed that P. fulvidraco has close relationships with Pelteobagrus eupogon and Tachysurus intermedius, suggesting that P. fulvidraco belongs to Tachysurus. This study provides evidence that Tachysurus, Pseudobagrus and Leiocassis do not form monophyly, but that these three genera form a monophyletic group. Our results provide reference for further phylogenetic research of the Bagridae species.     

Genetic Diversity of Tomistoma schlegelii Inferred from mtDNA Markers

The genetic diversity of the endangered crocodile Tomistoma schlegelii was characterized using the protein coding ND 6-tRNA^glu-cyt b and the cytochrome b-control region (cyt b-CR) markers. Concatenate data revealed six haplotypes with an overall haplotype diversity of 0.769 ± 0.039; nucleotide diversity was 0.00535 ± 0.00172. A nearest-neighbor analysis showed that all individuals clustered with four geographic regions (Sumatra, Peninsular Malaysia, Sarawak, and East Kalimantan) and were genetically differentiated. With the exception of the individuals from haplotype H2, which occurred in both Peninsular Malaysia and Sarawak, all other haplotypes were geographically distinct. The H4 lineage, which was found to be the most divergent, clustered exclusively in the basal clade in all phylogenetic trees, and the haplotype network was unconnected at the 95% reconnection limit, suggesting further investigation to establish its possible status as a distinct evolutionary significant unit or a cryptic species.     

Comparative Analysis of Genetic Diversity and Population Structure of Sipunculus nudus as Revealed by Mitochondrial COI Sequences

Genetic diversity and population structure of Sipunculus nudus were evaluated using a 652 base pair fragment of the mitochondrial cytochrome oxidase I gene. The populations were collected from Beihai, Sanya, and Xiamen. A total of 71 polymorphic sites defined 16 distinct haplotypes. The mean haplotype diversity and nucleotide diversity of the three populations were 0.9354 ± 0.0168 and 0.0035 ± 0.0018, respectively. Analysis at the intrapopulation level showed that the Beihai population had the greatest haplotype and nucleotide diversity, followed by the Xiamen and Sanya populations. Analysis of molecular variance showed significant genetic differentiation among the three populations (Fst = 0.0796, P < 0.05). The present results revealed that S. nudus populations had a high level of genetic diversity and distinct population structures.     

Genetic variation of the Chinese longsnout catfish Leiocassis longirostris in the Yangtze River revealed using mitochondrial DNA cytochrome b sequences

The mitochondrial DNA cytochrome b of 132 Leiocassis longirostris collected from 12 localities in the upper to lower reaches of the Yangtze River were amplified and partially sequenced using the PCR technique. The results showed that 27 nucleotide sites were variable along 817 bp length of homologous sequence (3.3%), base substitutions happened mostly at the third codon position. A total of 22 haplotypes were identified, which were characterized with moderate haplotype diversity (h = 0.5417 ± 0.0519), but low nucleotide diversity (π = 0.0019 ± 0.0012). Median-joining network analysis revealed star-shaped patterns with one common central haplotype (H3), whereas mismatch distribution analysis found that the Chinese longsnout catfish fitted a smooth unimodal distribution, which suggested that this species underwent population expansion following bottlenecks and/or they originated from a small number of founding individuals. The time that the total population of Chinese longsnout catfish in the Yangtze River expanded was estimated 169,000–337,000 years before present. The analysis of molecular variance (AMOVA) and phylogenetic reconstructions did not detect significant geographic structure between different river sections, especially between above and below the Gezhouba Dam and the Three Gorges Dam, which suggested that these recently developed dams might have not significantly resulted in population genetic differentiation in the Chinese longsnout catfish.     

Demographic History and Population Structure of Blackfin Flounder (Glyptocephalus stelleri) in Japan Revealed by Mitochondrial Control Region Sequences

The demographic history and population genetic structure of the blackfin flounder (Glyptocephalus stelleri) along coastal regions of Japan were investigated. Genetic variation in DNA sequences was examined from the first hypervariable region of the mitochondrial DNA control region. A high level of haplotypic diversity (h = 0.99 ± 0.004) was detected, indicating a high level of intrapopulation genetic diversity. The starburst structure of the minimum spanning tree suggested a very recent origin for most haplotypes. The demographic history of G. stelleri was examined using neutrality tests and mismatch distribution analysis, which also indicated a Pleistocene population expansion at about 124,100–413,400 years ago. Hierarchical molecular variance analysis and conventional population Fst comparisons revealed no significant genetic differentiation throughout the range examined.     

Genetic variation at mtDNA and microsatellite loci in Chinese longsnout catfish (<Emphasis Type="Italic">Leiocassis longirostris</Emphasis>)

Genetic variability and population structure of the Chinese longsnout catfish Leiocassis longirostris Günther in the Yangtze River was examined with mitochondrial control region sequences and nuclear microsatellite markers. A 705-bp segment of the mitochondrial DNA control region was sequenced from 132 samples, which identified a total of 61 haplotypes. The Chinese longsnout catfish in the Yangtze River was characterized with high haplotype diversity (h = 0.9770 ± 0.0041) but low nucleotide diversity (π = 0.0081 ± 0.0043). Median-joining network analysis revealed a star-shaped pattern and mismatch distribution analysis found a smooth unimodal distribution, which suggested that this species in the Yangtze River underwent a population expansion following bottlenecks and/or they originated from a small size of founding population. It was estimated that the possible time of population expansion was 139,000–435,000 years before present, a time period in the middle Pleistocene. The analysis of molecular variance and phylogenetic reconstructions did not detect significant geographic structure between different river sections. This pattern of genetic variation was further evidenced with nuclear microsatellite markers. The genetic differentiation between above and below the Gezhouba Dam and Three Gorges Dam is very small at mitochondrial and nuclear levels, which suggested that these recently developed dams might have not significantly resulted in population genetic fragmentation in the Chinese longsnout catfish. However, the potential exacerbation of genetic structuring by the dams should not be overlooked in the future.     

Analyses of genetic population structure of two ecologically important mangrove tree species, Bruguiera gymnorrhiza and Kandelia obovata from different river basins of Iriomote Island of the Ryukyu Archipelago, Japan

We analyzed nuclear and chloroplast microsatellite makers to assess genetic diversity and examine genetic structure of two mangrove tree species, Bruguiera gymnorrhiza and Kandelia obovata recovered from nine major river basins of Iriomote Island of the Ryukyu Archipelago, Japan. The average number of alleles per nuclear locus per population was 2.6 in B. gymnorrhiza and 1.7 in K. obovata. Bayesian clustering analysis using InStruct identified two genetic clusters in B. gymnorrhiza and three clusters in K. obovata. Chloroplast microsatellites revealed two dominant haplotypes from B. gymnorrhiza and three haplotypes from K. obovata. The overall result suggests low genetic diversity for both species. AMOVA for nuclear microsatellites showed 9.3?% of population variation in B. gymnorrhiza. Although genetic differentiation between several populations was significant in this species, F _ST suggested low to moderate level of differentiations between most of the populations. Distribution of genetic clusters and chloroplast haplotypes also suggested weak differentiations among B. gymnorrhiza populations. In K. obovata, population variation was, however, relatively high (27.8?%). The high differentiation between K. obovata populations was also suggested from the F _ST and genetic clusters from nuclear microsatellites, and chloroplast haplotypes. A significant correlation between chloroplast genetic distances and coastline distances as well as haplotype distributions for both species suggest that propagules from each species mostly disperse to the neighboring river basins. While significant F _IS and higher percentage of admixed clusters from nuclear microsatellites suggested inbreeding, continual gene exchange by propagule dispersal among populations especially among neighboring populations was suggested for both species from maternally inherited chloroplast microsatellites analyses.     

Karyotypic diversification in Hypostomus Lacépède, 1803 (Siluriformes, Loricariidae): biogeographical and phylogenetic perspectives

Hypostomus is a species-rich genus of fish with unclear systematics and phylogenetic relationships. Ten species of Hypostomus (H. albopunctatus, H. ancistroides, H. cochliodon, H. commersoni, H. faveolus, H. hermanni, H. aff. paulinus, H. regani, H. strigaticeps and H. topavae) were cytogenetically analyzed through Giemsa staining and silver nitrate impregnation, and the obtained data were correlated to the available biogeographical and phylogenetic analyses for the genus. Although the silver stained nucleolar organizer regions (AgNORs) were found to vary significantly among the species, the diploid numbers could be correlated to the distribution of the species on northern and southern South America river basins. Species with the lower diploid numbers (2n = 64) were associated to northern hydrographic basins and showed a single AgNORs bearing pair. Diploid numbers of 66–68 chromosomes and of 70–84 chromosomes were correlated to two major clades within Hypostomus and southern hydrographic basins, and showed AgNORs varying on number and position. Our results show that cytogenetic data can be correlated to the phylogeny and biogeography of the genus, helping to clarify its complex evolutionary history.     

Low genetic diversity and evidence of population structure among subspecies of Nerodia harteri, a threatened water snake endemic to Texas

Nerodia harteri is a threatened small-bodied water snake that occupies one of the most restricted ranges of any snake within the continental United States. It is found closely associated with rivers and tributaries in the Colorado and Brazos river basins, which flow through north-central Texas. Nerodia harteri has been at the center of debate owing to conflicts between conservation efforts and the construction of dams that change or destroy its preferred habitat. Additionally, its taxonomic status has also been under contention with some authors recognizing two subspecies, the Brazos water snake (N. h. harteri) and the Concho water snake (N. h. paucimaculata), whereas other authors consider each separate species. Despite its relatively recent discovery during the 1940s, N. harteri has been the subject of several ecological studies, yet no population genetic assessment of either subspecies has been performed to date. We first evaluated the phylogenetic placement of both subspecies among other North American Nerodia using partial sequence data from the mitochondrial gene cytochrome-b. We then tested for population subdivision among four rivers encompassing the range of N. harteri and tested for the presence of admixture between river basins using mitochondrial sequence data (920?bp of cyt-b) and five cross-species amplified microsatellite loci. We found low mitochondrial haplotype diversity represented by two unique haplotypes in each river basin, which were separated by no more than four nucleotide changes. Nuclear loci showed low genetic diversity and population structuring within and among river basins. We did not find conclusive evidence of admixture between basins, and we support the presence of two separate evolutionarily significant units and two separate management units corresponding to each major river basin. Given increasing natural and anthropogenic threats, we recommend continued ecological and genetic monitoring of both subspecies.     

Comparative Study on the Genetic Diversity of GHR Gene in Tibetan Cattle and Holstein Cows

Due to the phenotype-based artificial selection in domestic cattle, the underlying functional genes may be indirectly selected and show decreasing diversity in theory. The growth hormone receptor (GHR) gene has been widely proposed to significantly associate with critical economic traits in cattle. In the present study, we comparatively studied the genetic diversity of GHR in Tibetan cattle (a traditional unselected breed, n = 93) and Chinese Holstein cow (the intensively selected breed, n = 94). The Tibetan yak (n = 38) was also included as an outgroup breed. A total of 21 variants were detected by sequencing 1279 bp genomic fragments encompassing the largest exon 9. Twelve haplotypes (H1～H12) constructed by 15 coding SNPs were presented as a star-like network profile, in which haplotype H2 was located at the central position and almost occupied by Tibetan yaks. Furthermore, H2 was also identical to the formerly reported sequence specific to African cattle. Only haplotype H5 was simultaneously shared by all three breeds. Tibetan cattle showed higher nucleotide diversity (0.00215 ± 0.00015) and haplotype diversity (0.678 ± 0.026) than Holstein cow. Conclusively, we found Tibetan cattle have retained relatively high genetic variation of GHR. The predominant presence of African cattle specific H2 in the outgroup yak breed would highlight its ancestral relationship, which may be used as one informative molecular marker in the phylogenetic studies.     

Genetic diversity and population structure of the endangered ripon barbel,Barbus altianalis (Boulenger, 1900) in Lake Victoria catchment,Kenya based on mitochondrial DNA sequences

Approximately 850 bp of the mitochondrial control region was used to assess the genetic diversity, population structure and demographic expansion of the endangered cyprinid Barbus altianalis, a species known to be potamodramous in the Lake Victoria drainage system. The 196 samples taken from the four main rivers draining the Lake Victoria catchment (Nzoia, Yala, Nyando and Sondu–Miriu) yielded 49 mitochondrial DNA haplotypes; 83.7% thereof were private haplotypes restricted to particular rivers. The overall mean haplotype diversity was high (0.93663 ± 0.008) and ranged between 0.566 (Sondu – Miriu) and 0.944 (Nzoia). The overall mean nucleotide diversity was low (0.01322 ± 0.00141), ranging from 0.0342 (Sondu – Miriu) to 0.0267 (Nzoia). Population differentiation tests revealed strong and highly significant (P ≤ 0.001) segregation of populations in the four river basins. F_ST values among the four river‐based populations ranged from 0.05202 to 0.44352. The samples formed two main haplotype networks based on a 95% parsimony criterion, each exhibiting a strong signature of past population expansion. The smaller network was restricted to the River Nzoia, whereas the larger network contained representatives from all four rivers; within this the central haplotypes were found in more than one river, whereas the peripheral haplotypes tended to be river‐specific. The degree of population differentiation and the number of river‐specific haplotypes are too high to be explained by recent anthropogenic impacts alone and suggest that the species has probably existed in the Lake Victoria catchment as two populations: the now ‘extinct’ migratory population and the extant river restricted non‐migratory populations.     

Population genetics of Ageratina adenophora using inter-simple sequence repeat (ISSR) molecular markers in China

Abstract

Understanding distribution and diversity of invasive weeds is essential for the development of efficient control measures against it. In the present study, inter-simple sequence repeat (ISSR) markers were used to assess the biogeographic relationships among populations of the invasive Crofton weed (Ageratina adenophora (Spreng.)) during 2004–2006 in China. A total of 100 ISSR primers with di-, tri-, tetra- and penta-nucleotide repeats were screened, from which 20 polymorphic and informative primers were selected. Amplification of the 20 primers generated a total of 479 polymorphic bands among the 64 weed populations, and a high level of genetic diversity (H _E = 0.1541 ± 0.0193) was detected in A. adenophora. Neighbor-joining (NJ) cluster analysis based on genetic distances among populations grouped the populations according to their geographical origin, i.e. (1) populations of southwestern Guizhou, (2) populations of Liangshan city in Sichuan, (3) populations of western Guizhou, (4) Guangxi populations plus Chongqing populations, (5) populations of southern Yunnan, and (6) populations of Yangtze River Valleys in Sichuan plus populations of western Yunnan. A significant positive correlation between geographical and genetic distance was detected by the Mantel test (r = 0.183, p = 0.0012). Based on the divergence relationships revealed by ISSR markers, it was assumed that A. adenophora mainly dispersed through wind and water in China.     

The Phylogeographical Pattern and Conservation of the Chinese Cobra (Naja atra) across Its Range Based on Mitochondrial Control Region Sequences

The vulnerable Chinese cobra (Naja atra) ranges from southeastern China south of the Yangtze River to northern Vietnam and Laos. Large mountain ranges and water bodies may influence the pattern of genetic diversity of this species. We sequenced the mitochondrial DNA control region (1029 bp) using 285 individuals collected from 23 localities across the species'' range and obtained 18 sequences unique to Taiwan from GenBank for phylogenetic and population analysis. Two distinct clades were identified, one including haplotypes from the two westernmost localities (Hekou and Miyi) and the other including haplotypes from all sampling sites except Miyi. A strong population structure was found (Φst = 0.76, P<0.0001) with high haplotype diversity (h = 1.00) and low nucleotide diversity (π = 0.0049). The Luoxiao and Nanling Mountains act as historical geographical barriers limiting gene exchange. In the haplotype network there were two “star” clusters. Haplotypes from populations east of the Luoxiao Mountains were represented within one cluster and haplotypes from populations west of the mountain range within the other, with haplotypes from populations south of the Nanling Mountains in between. Lineage sorting between mainland and island populations is incomplete. It remains unknown as to how much adaptive differentiation there is between population groups or within each group. We caution against long-distance transfers within any group, especially when environmental differences are apparent.     

长江中下游湖泊短颌鲚线粒体控制区遗传多样性分析

为探明长江中下游不同湖泊中短颌鲚(Coilia brachygnathus)遗传多样性水平和遗传分化程度,以洞庭湖、长湖、巢湖3个地理群体作为研究对象,采用线粒体控制区序列为分子标记,分别应用软件Dna SP 5.0、Arlequin3.1.1、MEGA5.0和Network 5.1进行了遗传参数统计和单倍型间分子变异分析(AMOVA),构建邻接系统树及单倍型网络图。对长江中下游短颌鲚野生群体的遗传多样性和遗传结构进行分析。结果显示,用来分析的1 236 bp D-loop区序列中共90个变异位点,54个简约信息位点。长江中下游3个地理群体中共发现58个单倍型,单倍型多样性(h)范围0.949~0.982,核苷酸多样性范围0.004 99~0.006 21,说明长江中下游3个湖泊短颌鲚地理群体具有较高的遗传多样性水平。3个短颌鲚地理群体遗传分化指数(Fst)为0.265 95,呈现出中等程度的分化水平,主要表现在巢湖群体与其他群体之间处于中等程度分化水平。依据遗传距离构建系统发育树及单倍型网络图也出现相类似结果。     

Diversity of fish species in suichuan river and shushui river and conservation value,China

Based on the data collected from the Suichuan and Shushui Rivers, both tributaries of the Ganjiang River between April and July of 2015, the fish species identified were classified into 68 and 46 species, 14 and 12 families in the Suichuan River and Shushui River respectively. Cyprinidae is the most common family that accounts for 52.9% and 58.7% in the total number of fish species in the Suichuan River and Shushui River, respectively. The dominant species were Pseudohemiculter dispar, Squalidus argentatus, Silurus asotus and Leptobotia elongate for the Suichuan River and Squalidus argentatus, Acrossocheilus parallens, Pseudohemiculter dispar, Silurus asotus, Pelteobagrus fulvidraco, Acrossocheilus fasciatus and Zacco platypus for the Shushui River. The diversity of fish species showed more abundant and diverse for the Suichuan River. The composition of ecotype of fish indicated the rich diversities of ecotype in both the Suichuan and Shushui Rivers. The tributary (Suichuan River) and main stream of the Ganjiang River demonstrated a highly fauna similarity and the fish resources indicated its significance to maintain the fish diversity in the middle of the Ganjiang River and its benefits to the existence of the fish species in mountain streams. Dam construction, sand excavation and heavy metal pollution are the most significant threat to fish diversity and ecosystem functioning in Ganjiang River basin. In order to protect fish diversity and fisheries more effectively, relevant laws should be strengthened and conservation areas should be established for the survival of freshwater fish species.     

Common polymorphisms of calpain-10 and the risk of polycystic ovary syndrome in Tunisian population: a case–control study

Recent studies have suggested that calpain-10 (CAPN10) gene polymorphisms play a role in the susceptibility to polycystic ovary syndrome (PCOS). The aim of the present study was to evaluate the possible association between three single nucleotide polymorphisms (SNPs) in CAPN10 gene: UCSNP-43 (rs3792267), UCSNP-19 (rs3842570), and UCSNP-63 (rs5030952) and PCOS in Tunisian cases and control women. Study subjects included 127 women with PCOS (mean age 29.8 ± 4.7 year) and 150 healthy women (mean age 33.5 ± 5.6 year). CAPN10 genotyping was carried-out by direct PCR and PCR–RFLP. Linkage disequilibrium pattern in the genomic region explored was determined by HAPLOVIEW 4.2 while reconstruction of haplotypes was done using PHASE 2.1. The phylogenetic distribution of haplotypes in the population was determined by ARLEQUIN 2.000. Six haplotypes were observed. None of SNPs associated with PCOS or its components while the haplotype H4 associated with the phenotype PCOS-obese (P < 0.025). Moreover the pair of haplotypes H1/H4 strongly associated with high blood-pressure (OR = 14.4, P < 0.012). This work confirms the association of CAPN10 gene with metabolic components in PCOS and highlights the role of haplotypes as strong and efficient genetic markers.     

Genetic diversity and structure of an endemic and critically endangered stream river salamander (Caudata: Ambystoma leorae) in Mexico

Small or isolated populations are highly susceptible to stochastic events. They are prone and vulnerable to random demographic or environmental fluctuations that could lead to extinction due to the loss of alleles through genetic drift and increased inbreeding. We studied Ambystoma leorae an endemic and critically threatened species. We analyzed the genetic diversity and structure, effective population size, presence of bottlenecks and inbreeding coefficient of 96 individuals based on nine microsatellite loci. We found high levels of genetic diversity expressed as heterozygosity (H_o = 0.804, H_e = 0.613, H_e* = 0.626 and H_Nei = 0.622). The population presents few alleles (4–9 per locus) and genotypes (3–14 per locus) compared with other mole salamanders species. We identified three genetically differentiated subpopulations with a significant level of genetic structure (F_ST = 0.021, R_ST = 0.044 y D_est = 0.010, 95 % CI). We also detected a reduction signal in population size and evidence of a genetic bottleneck (M = 0.367). The effective population size is small (Ne = 45.2), but similar to another mole salamanders with restricted distributions or with recently fragmented habitat. The inbreeding coefficient levels detected are low (F_IS = ?0.619–0.102) as is gene flow. Despite, high levels of genetic diversity A. leorae is critically endangered because it is a small isolated population.     

GCPII modulates oxidative stress and prostate cancer susceptibility through changes in methylation of RASSF1, BNIP3, GSTP1 and Ec-SOD

Glutamate carboxypeptidase II (GCPII) haplotypes were found to influence susceptibility to prostate cancer. In the current study, we have elucidated the impact of these haplotypes on the expression of PSMA, BNIP3, Ec-SOD, GSTP1 and RASSF1 genes to understand the epigenetic basis of oxidative stress and prostate cancer risk. Expression analysis was carried out by RT-PCR. Bisulphite treated DNA was subjected to MS-PCR and COBRA for epigenetic studies. Plasma MDA and glutathione levels were measured. In prostate cancer, upregulation of BNIP3 (204.4 ± 23.77 vs. 143.9 ± 16.42 %, p = 0.03); and downregulation of Ec-SOD (105.8 ± 13.69 vs. 176.3 ± 21.1 %, p = 0.027) and RASSF1A (16.67 ± 16.0 vs. 90.8 ± 8.5 %, p = 0.0048) was observed. Hypomethylation of BNIP3 (31.25 ± 16.19 vs. 45.70 ± 2.42 %, p < 0.0001), hypermethylation of Ec-SOD (71.4 ± 6.75 vs. 10.0 ± 3.78 %, p < 0.0001) and RASSF1 (76.25 ± 12.53 vs. 30.0 ± 8.82 %, p = 0.0077) was observed in prostate cancer. The gene expression signature of PSMA, BNIP3, Ec-SOD, GSTP1, clearly demarcated cases and controls (AUC = 0.89 in the ROC curve). D191V variant of GCPII showed positive association with oxidative stress and inverse association with Ec-SOD expression. H475Y variant showed positive association with Ec-SOD expression and inverse association with oxidative stress. R190W variant was found to reduce oxidative stress by increasing glutathione levels. GCPII genetic variants contribute to increased oxidative stress and prostate cancer risk by modulating the CpG island methylation of Ec-SOD.     

Low genetic diversity and limited gene flow in a dominant mangrove tree species (Rhizophora stylosa) at its northern biogeographical limit across the chain of three Sakishima islands of the Japanese archipelago as revealed by chloroplast and nuclear SSR analysis

We examined the genetic diversity, population structure and gene flow in a dominant mangrove tree (Rhizophora stylosa) at its northern biogeographical limit in Sakishima islands of the Japanese archipelago. Simple sequence repeat (SSR) markers from chloroplast (cpSSR) and nuclear DNA were used to analyze 16 populations recovered from various river basins across the chain of three Sakishima islands—Iriomote, Ishigaki and Miyako. The average number of alleles (1.7–2.7) and observed heterozygosities (0.031–0.216) at nuclear SSR and haploid diversity (0.000–0.489) at cpSSR across the populations suggested low genetic diversity in R. stylosa in Sakishima islands. cpSSR analysis identified two haplotypes, and Bayesian clustering analysis (nuclear SSR) revealed two genetic clusters. Analysis of molecular variance (nuclear SSR) showed significant population differentiations. Pairwise tests consistently revealed significant differentiation between most of the population pairs; however, the degrees of differentiations are generally correspondent to the relative geographical distances as suggested from pairwise F _ST and cpSSR genetic distances. Moreover, Mantel tests showed some signals of correlations between genetic distances (nuclear and chloroplast) and geographical distances. These results suggest that combined contribution of gene flow via pollen and propagule dispersal in R. stylosa mostly occurred between neighboring river basins. The appearances of two cpSSR haplotypes (maternal lineages) as well as two nuclear genetic clusters (putative ancestral lineages) at various river basins support the hypothesis that present-day R. stylosa populations across the Sakishima islands were established from few identical founders; however, significant differentiations among various river basins most likely resulted from the limited gene flow and high inbreeding.