Population genetic structure in the endemic cyprinid fish Microphysogobio alticorpus in Taiwan: Evidence for a new phylogeographical area

The Formosa Bank is a significant zoogeographical barrier for the freshwater fish in Southern Taiwan. Here, we present population genetic structure analyses on Microphysogobio alticorpus populations in Taiwan, and biogeographic pattern to infer the relationship between historic dispersal geological dynamics. A total of 24 haplotypes were genotyped for mitochondrial CR + cyt b region in 134 specimens collected from 9 populations. Relatively high levels of haplotype diversity (h = 0.896) and low levels of nucleotide diversity (π = 0.005) were detected in M. alticorpus. Two major phylogenetic haplotype groups, Clade N1+N2 and S, were revealed and estimated to diverge about 0.121 myr (million years) ago. Haplotype network and population genetic analyses indicated significant genetic structure (F_ST = 0.876), largely concordant with the geographical location of the populations. According to SAMOVA analysis, M. alticorpus populations were divided into five units that can be explained by episodes of dispersal and population expansion followed by population fragmentation and restricted gene flow. The results of the topological pattern and network suggested that the Kaoping River was the first habitat colonized by the ancestor of M. alticorpus, before the northward dispersal.     

Comparative phylogeography of two codistributed endemic cyprinids in southeastern Taiwan

Two cyprinid fishes, Spinibarbus hollandi and Onychostoma alticorpus, are endemic to southeastern Taiwan. This study examined the phylogeography of these two codistributed primary freshwater fishes using mitochondrial DNA cytochrome b sequences (1140 bp) to search for general patterns in the effect of historical changes in southeastern Taiwan. In total, 135 specimens belonging to these two species were collected from five populations. These two codistributed species revealed similar genetic variation patterns. The genetic variation in both species was very low, and the geographical distribution of the genetic variation corresponded neither to the drainage structure nor to the geographical distances between the samples. The results of a statistical dispersal-vicariance analysis suggested that the ancestral populations of these two species were distributed in southern Taiwan before their dispersal. Our study suggests that the initial colonization occurred in the Kaoping River followed by eastern and northward dispersal. Our results also indicate that the Central Range in Taiwan did not act as a barrier to the dispersal of S. hollandi or O. alticorpus.     

Population genetic structure in wild and hatchery populations of white cloud mountain minnow (Tanichthys albonubes): Recommendations for conservation

We used mitochondrial (mt) cytochrome b gene (cyt b) to compare the genetic variability in three hatchery broodstocks of white cloud mountain minnow with the variability in six wild populations sampled in two river drainages. A total of 43 haplotypes in 102 specimens were observed, with no haplotype shared between wild and hatchery populations. The nucleotide diversity of the wild samples (0.048) was significantly higher than that of the hatchery ones (0.007), but the haplotype diversity was almost similar between them. Two major phylogenetic haplotype groups were revealed and estimated to diverge about 6.531 myr (million years) ago. Significant genetic differentiation was revealed between wild and hatchery populations as well as among nine sampled populations, suggesting at chance effect during the founding process for the hatchery population and a subsequent genetic drift. According to the network, the connection between wild and hatchery populations indicates that present hatchery populations originated from single wild population. We suggested that two regions (Pearl River system and Lu River) identified by reciprocal mtDNA monophyly and SAMOVA should be regarded as three different ESUs and two different MUs in South China, respectively.     

Phylogeography of lateral plate dimorphism in the freshwater type of ninespine sticklebacks, genus Pungitius

The freshwater type of ninespine sticklebacks, genus Pungitius, is widely distributed in northern Japan and reproductively isolated from other genetically divergent types endemic to small regions in Japan. This type expresses dimorphism in its lateral plate morphology: complete and partial row morphs. The two morphs show a parapatric distribution in Japan. To clarify the process involving the distribution of these two morphs, we examined their phylogeography based on restriction fragment length polymorphism of an entire mitochondrial DNA (mtDNA). The survey was carried out with seven restriction enzymes on the populations of the freshwater type collected from 41 localities across the distribution range in Japan, and 6 further Pungitius populations from the Okhotsk Sea coast of Russia were appended. An unweighted pair-group method with arithmetic averages (UPGMA) tree among 54 mtDNA haplotypes resolved eight clustering groups that differed in sequence divergence by approximately 1.3%–2.1%. Two of the eight groups were found only in Russia. mtDNA phylogenies constructed by neighbor-joining and Wagner parsimony methods suggested that the haplotypes of each plate morph were polyphyletic. The geographic distribution pattern of these groups suggests that they should be classified into two broad categories, one with extensive distribution and the other with localized distribution of the constituent haplotypes within a group. The former groups were found mainly in the populations with the completely plated morph and the latter groups with the partially plated morph. It is supposed that twice dispersals of dimorphic or complete plated ancestors and genetic differentiation during the interglacial played an important role in the formation of the present distribution of the two morphs in Japan. Received: March 28, 2000 / Revised: November 3, 2000 / Accepted: January 16, 2001     

Phylogeography of Asian weaver ants, Oecophylla smaragdina

The phylogeography of Oecophylla smaragdina was studied using the mitochondrial cytochrome b gene (Cytb), cytochrome oxidase subunit I (COI), and nuclear long-wavelength opsin gene (LW Rh). Weaver ants were collected from 35 localities and from one to nine colonies per locality. Neighbor-joining trees inferred from 647 bp of Cytb and 1,026 bp of COI using Oecophylla longinoda as an outgroup indicated that the haplotypes of O. smaragdina were clearly separated into seven groups: group 1 of India excluding West Bengal; group 2 of Bengal, Indochinese Peninsula, Malay Peninsula and Greater Sunda Islands, including Lombok and Sumbawa; group 3 of the Philippines; group 4 of Flores; group 5 of Sulawesi; group 6 of Halmahera; and group 7 of New Guinea and Australia. This grouping was also supported by a strict consensus tree derived from maximum parsimony and maximum likelihood trees. In addition, two haplotypes of LW Rh were found in O. smaragdina: one in group 2 and another in all the other groups. Comparison to haplotypes in other hymenopteran species suggests that group 2 is younger than other groups of O. smaragdina. The clustering of the seven groups was coincident with geological evidence of the distribution of continents, islands, and seas during glacial periods.     

Phylogeography of the endangered tideland snail Batillaria zonalis in the Japanese and Ryukyu Islands

The phylogeography of the endangered tideland snail Batillaria zonalis in the Japanese and Ryukyu Islands was analyzed on the basis of nucleotide sequences of a mitochondrial gene for cytochrome oxidase c subunit I (COI). Extremely low genetic diversity was found in populations at both the northern and southern boundaries of the geographic distribution of this species in Japan, i.e., Sendai Bay and Iriomotejima Island, respectively, which might be attributed to the population bottleneck due to historical environmental variations and/or the recent foundation of populations in the marginal part of the inhabitable range. Most populations contained unique rare haplotypes, and significant genetic differentiation on the whole was shown, while no clear geographic genetic structure was detected between the Japanese and Ryukyu Islands or over the distribution area of B. zonalis in Japan, with the exception of significant genetic divergences in Ago Bay in the central part of Honshu and the southern part of Okinawajima Island.     

Phylogeography of the cyprinid Squalius aradensis and implications for conservation of the endemic freshwater fauna of southern Portugal

The Iberian cyprinid fauna, characterized by the presence of numerous endemic species, has suffered from significant habitat degradation. The critically endangered Squalius aradensis is restricted to small drainages of southern Portugal, habitats that typically exhibit a characteristic Mediterranean-type heterogeneous hydrological system throughout the year, including alternation of flooding events during winter and complete drought in large river sections during summer. To assess the effect of historical and recent processes on genetic diversity in S. aradensis we examined within- and among-population variability in cytochrome b and six polymorphic microsatellite loci. Estimates of genetic diversity in time and space through the combined use of traditional Phi-/F-statistics, phylogenetic trees, ordination methods and nested clade analysis indicated significant and congruent structuring among populations. Data suggest that the Arade drainage represent the evolutionary centre of the species, with subsequent allopatric fragmentation across drainages. Factors other than isolation by distance strongly affected the within-drainage genetic differentiation observed in these Mediterranean-type drainages, including recent population expansion from a bottleneck event and restricted gene flow imposed by a long-term barrier (brackish water area). Significant correlation was found between S. aradensis allelic diversity and upstream drainage area. The relevance of findings for conservation issues is discussed in relation to local intermittent hydrological conditions, the highly restricted distribution and the critically endangered status of the species.     

Mitochondrial DNA reveals a mosaic pattern of phylogeographical structure in Atlantic and Mediterranean populations of anchovy (Engraulis encrasicolus)

This study extends the geographic coverage of a previous study of mitochondrial DNA restriction fragment length polymorphism in European anchovy. Both studies together include 24 samples representing 17 localities extending from the Black Sea, through the Mediterranean Sea to the eastern Atlantic as far south as Dakar, Senegal. Eighty-eight haplotypes define two clades (A and B) separated by 3.2% sequence divergence. Clade A has a star-like genealogy indicative of a recent population expansion. Clade B has a more complex genealogy, consisting of several haplotypes at intermediate frequencies. The distributions of these clades consist of a mosaic with abrupt changes between some areas and gradients between other areas. Clade A predominates the Black and Aegean seas, but is present throughout the Mediterranean. Unexpectedly, new data show that clade A is also at a high frequency in the Atlantic, from Portugal to at least Senegal. Overall, the level of genetic differentiation among populations is high (F(ST)=0.148, p<0.0001), with the greatest differences between basins. AMOVA reveals four main geographical groups: Atlantic, central Mediterranean, Aegean Sea, and Black Sea. Mismatch distribution clearly indicates historical bottleneck and population expansion for clade A, while for clade B such evidence is equivocal. This difference may reflect a range expansion for both clades, but with higher gene flow (Nm values) between demes for clade A. Both contemporary and historical processes are important in shaping the complex genetic population structure of European anchovy.     

Phylogenetic relationships and genetic diversity of badgers from the Korean Peninsula: Implications for the taxonomic status of the Korean badger

Accurate taxonomic classification of wildlife species is crucial for guiding biological research and for developing effective management and conservation programs. The taxonomic status of Eurasian badgers from South Korea remains poorly resolved. Here we assessed the phylogenetic relationships and genetic variation of Eurasia badgers using partial mitochondrial fragments to elucidate the evolutionary history and taxonomic status of badgers from the Korean Peninsula. Forty-eight unique haplotypes from 125 individuals were observed. Phylogenetic reconstructions and reduced median networks indicate that Eurasian badgers consisted of four geographic clades (Japan, Eastern Eurasia, Western Eurasia, and Caucasus) with a relatively weak split observed within Eastern Eurasia. Estimated divergence time between the Japanese and Eastern Eurasian clades, including the Korean population, was 467,100 years (69,200–1,085,500 years). The results of this study support the hypothesis that the Japanese badger migrated from the Eurasian continent over the Korea-Japan land bridge and that the Korean Peninsula was an important refugia during the Pleistocene. Our study confirmed that the South Korean badger, Meles meles, belongs to the Eastern Eurasian clade. Based on these results and those of previous studies, we recommend that the scientific name of the Korean badger be changed from M. meles to Meles leucurus (Asian badger).     

Genetic population structure of <Emphasis Type="Italic">Lefua echigonia</Emphasis> inferred from allozymic and mitochondrial cytochrome <Emphasis Type="Italic">b</Emphasis> variations

We examined the genetic population structure of Lefua echigonia (Japanese name, hotoke-dojo) using polymerase chain reaction restriction-fragment length polymorphisms (PCR-RFLPs) of the mitochondrial cytochrome b gene and two allozymic loci. The phylogenetic relationships of L. echigonia and those among L. echigonia, Lefua sp. (nagare hotoke-dojo), and L. nikkonis (ezo hotoke-dojo) were also investigated based on the nucleotide sequences of the cytochrome b gene. PCR-RFLP analysis revealed 18 mitotypes in L. echigonia, 2 in Lefua sp., and 1 in L. nikkonis. Phylogenetic trees based on the cytochrome b sequences indicated that the 18 mitotypes in L. echigonia were divided into five distinct groups (South-Kanto, North-Kanto, Tohoku, Echigo, and Tokai-Kinki clades) that differed by 8.5–15.3%, reflecting region-specific geographic distributions. The distributions of alleles in two allozymic loci roughly corresponded to those of the mitotype groups. The divergence times of the five groups were estimated to be about 3.4–7.7 million years ago by applying a general rate for mitochondrial DNA, suggesting that the divergence among them might have occurred in the late Tertiary. It can be inferred that the regional differentiation of each group was mainly due to geographic isolation and that this has been maintained, because the boundaries among the groups corresponded to geological features. The trees also supported the existence of three taxa, L. echigonia, Lefua sp., and L. nikkonis. We concluded that Lefua sp. was distinguished from other species in Lefua by morphological and ecological characters and also by genetic divergences of the cytochrome b gene. Our study also demonstrated the superior efficacy and simplicity of PCR-RFLP analysis as a method for detecting genetic variation in L. echigonia.     

Genetic structure of the endangered Irrawaddy dolphin (Orcaella brevirostris) in the Gulf of Thailand

The Irrawaddy dolphin (Orcaella brevirostris) is an endangered, small cetacean species which is widely distributed in rivers, estuaries, and coastal waters throughout the tropical and subtropical Indo-Pacific. Despite the extensive distribution of this species, little is known of individual movements or genetic exchange among regions in Thailand. Here, we evaluate the genetic diversity and genetic structure of O. brevirostris in the eastern, northern and western Gulf of Thailand, and Andaman Sea. Although phylogenetic relationships and network analysis based on 15 haplotypes obtained from 32 individuals reveal no obvious divergence, significant genetic differentiation in mitochondrial DNA (overall F_ST = 0.226, P < 0.001; Φ_ST = 0.252, P < 0.001) is apparent among regions. Of 18 tested microsatellite loci, 10 are polymorphic and successfully characterized in 28 individuals, revealing significant genetic differentiation (overall F_ST = 0.077, P < 0.05) among the four sampling sites. Structure analysis reveals two inferred genetic clusters. Additionally, Mantel analysis demonstrates individual-by-individual genetic distances and geographic distances follow an isolation-by-distance model. We speculate that the significant genetic structure of O. brevirostris in Thailand is associated with a combination of geographical distribution patterns, environmental and anthropogenic factors, and local adaptations.     

Phylogeography of the temperate tree species Quercus acutissima in China: Inferences from chloroplast DNA variations

Quercus acutissima is one of the most widespread temperate deciduous tree species in China. To study its phylogeographical pattern and demographic history, three chloroplast DNA fragments (atpB-rbcL, psbA-trnH and trnS-trnG) from 401 individuals representing 30 populations were sequenced. A total of 19 haplotypes were identified, and these showed a weak phylogeographical structure (N_ST = 0.689 > G_ST = 0.630, P > 0.05) at the species level. The Q. acutissima population harboured a high level of genetic diversity (H_T = 0.791), and the genetic variation mainly resided among populations (59.54%). The unimodal mismatch distribution and significantly negative Fu's F_S value indicate that the Q. acutissima population experienced rapid range expansion, which probably occurred between 0.37 and 0.12 Ma. Molecular phylogeography and ecological niche modelling (ENM) data suggest the existence of multiple localized glacial refugia in central China (e.g., the Qinling, Dabashan and Dabieshan mountain ranges) and southwestern China (Yunnan-Guizhou Plateau and its adjacent regions) during the Quaternary glaciations. Our study showed that geographical heterogeneity and climate changes may have shaped the genetic structure and phylogeographical pattern of this tree species in China.     

Phylogeography and intraspecific genetic variation of prochilodontid fishes endemic to rivers of northern South America

Phylogeography and intraspecific genetic variation were studied in prochilodontids endemic to the Orinoco, Essequibo and Amazon River basins of northern South America. Portions of two protein‐encoding mitochondrial (mt) DNA genes, ND4 and COI, were examined using single‐strand conformational polymorphisms (SSCPs) and nucleotide sequencing. Phylogeographic analysis indicated that the geographically widespread Prochilodus rubrotaeniatus is paraphyletic, with individuals from the Orinoco sharing most recent common ancestry with co‐occurring Prochilodus mariae . A second Prochilodus rubrotaeniatus clade was composed of haplotypes found in the Rio Cuyuní(Essequibo Basin) and tributaries of the Rio Negro (Amazon Basin). Intraspecific genetic analysis suggested that a complex set of processes have influenced patterns of genetic variation in prochilodontid lineages. Prochilodus rubrotaeniatus is monomorphic at both loci in the Rio Negro and probably recently colonized this basin from the Rio Essequibo. Only two of 55 P. mariae exhibited variant haplotypes, and both had resulted from non‐synonymous changes in the ND4 region. These observations were counter to neutral expectation and consistent with the action of natural selection on the mitochondrion. Overall, these analyses implicate vicariance, demography and selection for driving diversification of prochilodontids in northern South America.     

猪科的系统发育与系统地理分化

猪科动物是最为昌盛的哺乳动物之一,但它们的分类及系统发育关系却尚待解决。为此,用PCR直接测序法测定了一个红河猪(Potamochoerusporcus)、一个马来西亚野猪(Susbarbatus)以及数头欧亚野猪(Susscorfa)线粒体细胞色素b全序列1140bp。结合从GenBank中获得的其他猪科物种胞色素b序列,用邻接法和最大简约法构建了猪科物种的系统发育关系,其结果显示苏拉威西鹿猪与其他猪科物种为姐妹群的关系。非洲物种疣猪和河猪聚为一枝,而欧亚猪属物种聚为另一支。各物种在系统发育关系中的位置与它们的地理分布对应。猪属物种可以被划分为3个种组:分布于菲律宾群岛、苏拉威西岛及其邻近岛屿的猪属物种形成一个种组,包括苏拉威西岛疣猪、菲律宾疣猪和印度尼西亚疣猪;欧亚野猪种组,包括欧亚野猪和姬猪;须猪和爪哇疣猪形成一个种组。     

黑嘴鸥种群的遗传分化与迁徙模式

用分子生物学方法研究了黑嘴鸥(Larus saundersi)种群的遗传结构和分化.对我国黑嘴鸥现存的3个地方种群50个样本的线粒体DNA控制区部分片段进行了序列测定,结合已公布的的两个韩国样本的序列,比对后所获得的550 bp的序列中,共发现49个变异位点,定义了37种线粒体单元型,其中仅1种为我国3个地方种群间共享单元型,多数为各地方种群内特有.我国3个地方种群内的单元型多样性和核苷酸多态性都很高;4个地方种群合并成一个大种群时, 其线粒体单元型多样性与核苷酸多样性仍很高,分别为0.974%±0.012%和0.510%±0.042%.进一步分析表明,4个地方种群间核苷酸歧异度在0.433%-0.585%之间,种群分化程度FST在-0.03176-0.48063之间,基因流Nm在1.03-33.79之间,平均遗传距离为0.0051±0.0011.UPGMA系统发生树及最大简约网络图表明, 黑嘴鸥的演化关系呈现星状辐射,单元型之间呈现一种混杂的分布格局,地方种群没有形成明显的系统地理结构.同时运用彩色旗标技术,初步掌握了4个地方种群迁徙规律.对比分子研究结果,验证了不同地方种群间的遗传分化与其迁徙模式基本一致,同时认为黑嘴鸥不同地方繁殖种群从越冬地返回繁殖地的过程中会发生交叉.最后,鉴于黑嘴鸥不同地方种群间的遗传分化不显著,其迁徙行为具有差异性,建议在保护和遗传管理中应分别对待.     

Phylogeography and genetic structure of disjunct Salix arbutifolia populations in Japan

Disjunct geographic distributions of boreal plant species at the southern edges of their ranges are expected to result in low genetic diversity and high genetic differentiation in the disjunct populations. This prediction was tested in a riparian willow, Salix arbutifolia, distributed in the northeastern Eurasian continent and the Sakhalin, Hokkaido, and Honshu Islands, using chloroplast DNA haplotypes and nuclear microsatellite genotypes. Hokkaido and Honshu populations shared a chloroplast haplotype identical to a closely related species, S. cardiophylla. This haplotype was divergent from haplotypes in the Eurasian continent (Primorsky) and the Sakhalin Island. In the nuclear microsatellites, most Hokkaido populations were genetically closer to Primorsky populations than to Sakhalin populations in spite of the geographical vicinity between Sakhalin and Hokkaido. The unexpected genetic divergence between Sakhalin and Hokkaido implies a complicated history of migration and colonization. The most peripheral populations in Honshu had the lowest genetic diversity and were most differentiated from the others. Thus, low genetic diversity and high genetic differentiation at the range periphery suggest substantial effects of genetic drift on genetic structure in the disjunct populations of Salix arbutifolia at the southern edge of its range.     

Migratory behaviour shapes spatial genetic structure of cyprinid fishes within the Lake Malawi catchment

• Genetic differences among freshwater fish populations are dependent on life‐history characteristics of the species, including the range of adult dispersal and the extent of homing to natal breeding grounds. However, the effects of variation in such characteristics on population genetic connectivity are rarely studied comparatively among closely related species.
• We studied population genetic structure within three congeneric cyprinid species from the Lake Malawi catchment that differ substantially in life‐history traits and conservation status, using a combination of microsatellite and mitochondrial DNA markers. Mpasa (Opsaridium microlepis) is a large (70 cm total length) migratory species that spawns in rivers, but as an adult is exclusively known from the main lake body. Sanjika (Opsaridium microcephalum), is a medium size (30 cm total length) species that exists in lake breeding, river‐lake migratory and apparently landlocked populations. Dwarf sanjika (Opsaridium tweddleorum) is a small non‐migratory species (15 cm total length) that persists in small tributaries surrounding the main lake and adjoining rivers.
• The results revealed striking differences among the three species in spatial genetic structuring. The river‐lake migratory mpasa showed only weak yet significant population genetic structure within the main Lake Malawi catchment, suggesting that there is no strong natal homing. The habitat‐generalist sanjika showed only weak spatial genetic differentiation at microsatellite loci within the Lake Malawi catchment, but moderate structure in mitochondrial DNA, potentially reflecting male‐biased dispersal. The river‐restricted dwarf sanjika showed strong genetic structure in both microsatellite and mitochondrial DNA, suggesting strictly limited dispersal at both adult and juvenile stages.
• We conclude that contrasting migration life histories have resulted in dramatically different patterns of population genetic structure among these congeneric species. The observed patterns demonstrate how divergent life‐history evolution may strongly influence broader patterns of population genetic connectivity in freshwater fish, with consequences for management and conservation. Specifically the results suggesting gene flow among Lake Malawi populations of mpasa, an IUCN red‐listed ‘Endangered’ species endemic to the lake catchment, imply that conservation initiatives operating at both local and catchment scales are needed to reverse local population decline.

     

Phylogeography of the frog Leptodactylus validus (Amphibia: Anura): Patterns and timing of colonization events in the Lesser Antilles

The frog Leptodactylus validus occurs in northern South America, Trinidad and Tobago, and the southern Lesser Antilles (Grenada and St. Vincent). Mitochondrial DNA sequences were used to perform a nested clade phylogeographic analysis (NCPA), to date colonization events, and to analyze colonization patterns using on a relaxed molecular clock and coalescent simulations. L. validus originated on the mainland and first colonized Trinidad with subsequent independent colonizations of Tobago and the Lesser Antilles from Trinidad. The NCPA suggests a historical vicariant event between populations in Trinidad and Tobago from those in the Lesser Antilles. The colonization of Trinidad occurred 1 million years ago (mya) and the colonization of the Lesser Antillean islands occurred 0.4 mya. The coalescent approach supported the scenario where L. validus dispersed from Trinidad to St. Vincent and from there to Grenada, a dispersal event that could have been mediated by human introduction as recent as 1600 years ago.     

Phylogeography and genetic diversity of the red squirrel (Sciurus vulgaris) in China: Implications for the species' postglacial expansion history

The geographical distribution of Sciurus vulgaris spans much of the Palearctic from western Europe and the UK eastward to the pacific coast of East Asia. S. vulgaris occurs in China only in the far northwest and the northeast. Understanding of the species’ postglacial expansion history in East Asia has been limited by the paucity of molecular data. In this study, we used partial D-loop and cytochrome b gene sequences to assess mitochondrial DNA variation in S. vulgaris in China. Our objectives were to (1) determine phylogeographical patterns of S. vulgaris in China; (2) understand the species’ postglacial expansion history in this region; and (3) quantify genetic diversity levels within S. vulgaris populations in China. We identified a supported phylogenetic group in S. vulgaris from China, and found no tendency for haplotypes to cluster by geographic region. Our analysis of S. vulgaris from China and other regions supports the hypothesis that the Calabria region of southern Italy is a glacial refugium for the species. We tentatively propose a postglacial expansion pattern for the squirrels: migrating from Calabria via Central and Eastern Europe to Russia and from there to China, and firstly to the northwest and then to northeast in China. We found high levels of genetic diversity in S. vulgaris populations across China as a whole, and discussed its influential factors.