Swedish population substructure revealed by genome-wide single nucleotide polymorphism data

The use of genome-wide single nucleotide polymorphism (SNP) data has recently proven useful in the study of human population structure. We have studied the internal genetic structure of the Swedish population using more than 350,000 SNPs from 1525 Swedes from all over the country genotyped on the Illumina HumanHap550 array. We have also compared them to 3212 worldwide reference samples, including Finns, northern Germans, British and Russians, based on the more than 29,000 SNPs that overlap between the Illumina and Affymetrix 250K Sty arrays. The Swedes - especially southern Swedes - were genetically close to the Germans and British, while their genetic distance to Finns was substantially longer. The overall structure within Sweden appeared clinal, and the substructure in the southern and middle parts was subtle. In contrast, the northern part of Sweden, Norrland, exhibited pronounced genetic differences both within the area and relative to the rest of the country. These distinctive genetic features of Norrland probably result mainly from isolation by distance and genetic drift caused by low population density. The internal structure within Sweden (F_ST = 0.0005 between provinces) was stronger than that in many Central European populations, although smaller than what has been observed for instance in Finland; importantly, it is of the magnitude that may hamper association studies with a moderate number of markers if cases and controls are not properly matched geographically. Overall, our results underline the potential of genome-wide data in analyzing substructure in populations that might otherwise appear relatively homogeneous, such as the Swedes.     

Dissecting the genetic structure of Korean population using genome-wide SNP arrays

Genome-wide SNP arrays have generated unprecedented quantities of data allow the detection of human evolutionary history and dense genome-wide data also enable the identification of distance ancestry among individuals or ethnic groups. To explain wider aspects of the genetic structure of Koreans and the East Asian population, we analyzed 79 individuals from the Korean HapMap project at 555,352 common single-nucleotide polymorphism loci, and compared this data with the worldwide population groups with the 53 ethnic groups from Human Genome Diversity Panel (HGDP-CEPH). Population differentiation (F_ST), Principal Component Analyses, STRUCTURE and ADMIXTURE are examined. In general, all the individual samples studies here were classified into subset of ethnic groups according to their geographical origins. Korean HapMap individuals were grouped together with East Asian populations from HGDP panel. Recently, a sub-population structure within Korean population has been reported. Our result, however, revealed the genetic homogeneity of Korean population. The ADMIXTURE analysis showed that, overall the Korean populations derive 79 % of their genomic ancestry from southern Asia and have relatively little northern Asian ancestry (21 %). The present work, therefore, provide the evidence that the male-biased southern-to-northern migration influenced not only for the genetic make up of the Y chromosome in the Korean population but also, its autosomal composition.     

中国南北方汉族人群DNA甲基化表观遗传差异研究

目的 族群地域、体貌特征等表型是基因型与环境共同作用的结果。大量基因组学研究表明,汉族人群具有混合特征,内部存在明显的南北遗传差异。本研究旨在探索研究表观基因组在中国南北方汉族人群之间是否存在差异,并筛选差异遗传位点。方法 使用GLINT软件对483份汉族样本的全基因组甲基化芯片数据进行EWAS分析,使用Lasso回归方法筛选位点。使用多元逻辑回归算法构建南北方汉族人群预测模型,通过十折交叉验证的方法评估。结果 筛选出一组南北方汉族之间差异显著的CpG位点,准确性为99.03%,Kappa系数为0.979 6。结论 本研究表明南北方汉族人群之间存在表观遗传差异,本研究为进一步开展不同地域汉族人群之间的表观遗传差异研究奠定了基础。     

Genetic structure of desert ground squirrels over a 20-degree-latitude transect from Oregon through the Baja California peninsula

The genetic structure of populations over a wide geographical area should reflect the demographic and evolutionary processes that have shaped a species across its range. We examined the population genetic structure of antelope ground squirrels (Ammospermophilus leucurus) across the complex of North American deserts from the Great Basin of Oregon to the cape region of the Baja California peninsula. We sampled 73 individuals from 13 major localities over this 2500-km transect, from 43 to 22 degrees north. Our molecular phylogeographical analysis of 555 bp of the mitochondrial cytochrome b gene and 510 bp of the control region revealed great genetic uniformity in a single clade that extends from Oregon to central Baja California. A second distinct clade occupies the southern half of the peninsula. The minimal geographical structure of the northern clade, its low haplotype diversity and the distribution of pairwise differences between haplotypes suggest a rapid northward expansion of the population that must have followed a northward desert habitat shift associated with the most recent Quaternary climate warming and glacial retreat. The higher haplotype diversity within the southern clade and distribution of pairwise differences between haplotypes suggest that the southern clade has a longer, more stable history associated with a southern peninsular refugium. This system, as observed, reflects both historical and contemporary ecological and evolutionary responses to physical environmental gradients within genetically homogeneous populations.     

Northern gene flow into southeastern East Asians inferred from genome-wide array genotyping

The population history of Southeast (SE) China remains poorly understood due to the sparse sampling of present-day populations and limited modeling with ancient genomic data. We report genome-wide genotyping data from 207 present-day Han Chinese and Hmong-Mien (HM)-speaking She people from Fujian and Taiwan Island, SE China. We coanalyzed 66 Early Neolithic to Iron Age ancient Fujian and Taiwan Island individuals obtained from previously published works to explore the genetic continuity and admixture based on patterns of genetic variations of the high-resolution time transect. We found the genetic differentiation between northern and southern East Asians was defined by a north–south East Asian genetic cline and our studied southern East Asians were clustered in the southern end of this cline. The southeastern coastal modern East Asians are genetically similar to other southern indigenous groups as well as geographically close to Neolithic-to-Iron Age populations, but they also shared excess alleles with post-Neolithic Yellow River ancients, which suggested a southward gene flow on the modern southern coastal gene pool. In addition, we identified one new HM genetic cline in East Asia with the coastal Fujian HM-speaking She localizing at the intersection between HM and Han clines. She people show stronger genetic affinity with southern East Asian indigenous populations, with the main ancestry deriving from groups related to southeastern ancient indigenous rice farmers. The southeastern Han Chinese could be modeled with the primary ancestry deriving from the group related to the Yellow River Basin millet farmers and the remaining from groups related to rice farmers, which was consistent with the northern China origin of modern southeastern Han Chinese and in line with the historically and archaeologically attested southward migrations of Han people and their ancestors. Our estimated north–south admixture time ranges based on the decay of the linkage disequilibrium spanned from the Bronze Age to historic periods, suggesting the recent large-scale population migrations and subsequent admixture participated in the formation of modern Han in SE Asia.     

Randomly Amplified Polymorphic DNA Variation in Populations of Eastern Australian Koalas, Phascolarctos cinereus

Randomly amplified polymorphic DNA (RAPD) variation in populations of the koala, Phascolarctos cinereus, was investigated, revealing significant differences in the level of diversity between southern and northern regions of eastern Australia. Of the 20 polymorphic RAPD markers identified in koalas, 4-7 were polymorphic in southern populations, while 12-17 were polymorphic in northern populations. Analysis of molecular variance revealed a significant difference in the estimated variance between koalas from northern and those from southern regions (P < 0.001), where populations from the north were greater than twice as variable as their southern cousins. The total genetic diversity observed was attributed to regional differences (30.91%), population differences within a region (11.77%), and differences among individuals within a population (57.32%). For the within-region analyses, a large proportion of the genetic diversity was attributable to individual differences within a population, 80.34% for the north and 91.23% for the south. These results demonstrate that RAPD markers are useful for determining population structure among koalas.     

A genome-wide analysis of admixture in Uyghurs and a high-density admixture map for disease-gene discovery

Following up on our previous study, we conducted a genome-wide analysis of admixture for two Uyghur population samples (HGDP-UG and PanAsia-UG), collected from the northern and southern regions of Xinjiang in China, respectively. Both HGDP-UG and PanAsia-UG showed a substantial admixture of East-Asian (EAS) and European (EUR) ancestries, with an empirical estimation of ancestry contribution of 53:47 (EAS:EUR) and 48:52 for HGDP-UG and PanAsia-UG, respectively. The effective admixture time under a model with a single pulse of admixture was estimated as 110 generations and 129 generations, or admixture events occurred about 2200 and 2580 years ago for HGDP-UG and PanAsia-UG, respectively, assuming an average of 20 yr per generation. Despite Uyghurs' earlier history compared to other admixture populations, admixture mapping, holds promise for this population, because of its large size and its mixture of ancestry from different continents. We screened multiple databases and identified a genome-wide single-nucleotide polymorphism panel that can distinguish EAS and EUR ancestry of chromosomal segments in Uyghurs. The panel contains 8150 ancestry-informative markers (AIMs) showing large frequency differences between EAS and EUR populations (F_ST > 0.25, mean F_ST = 0.43) but small frequency differences (7999 AIMs validated) within both populations (F_ST < 0.05, mean F_ST < 0.01). We evaluated the effectiveness of this admixture map for localizing disease genes in two Uyghur populations. To our knowledge, our map constitutes the first practical resource for admixture mapping in Uyghurs, and it will enable studies of diseases showing differences in genetic risk between EUR and EAS populations.     

A Genomic Portrait of Haplotype Diversity and Signatures of Selection in Indigenous Southern African Populations

We report a study of genome-wide, dense SNP (∼900K) and copy number polymorphism data of indigenous southern Africans. We demonstrate the genetic contribution to southern and eastern African populations, which involved admixture between indigenous San, Niger-Congo-speaking and populations of Eurasian ancestry. This finding illustrates the need to account for stratification in genome-wide association studies, and that admixture mapping would likely be a successful approach in these populations. We developed a strategy to detect the signature of selection prior to and following putative admixture events. Several genomic regions show an unusual excess of Niger-Kordofanian, and unusual deficiency of both San and Eurasian ancestry, which were considered the footprints of selection after population admixture. Several SNPs with strong allele frequency differences were observed predominantly between the admixed indigenous southern African populations, and their ancestral Eurasian populations. Interestingly, many candidate genes, which were identified within the genomic regions showing signals for selection, were associated with southern African-specific high-risk, mostly communicable diseases, such as malaria, influenza, tuberculosis, and human immunodeficiency virus/AIDs. This observation suggests a potentially important role that these genes might have played in adapting to the environment. Additionally, our analyses of haplotype structure, linkage disequilibrium, recombination, copy number variation and genome-wide admixture highlight, and support the unique position of San relative to both African and non-African populations. This study contributes to a better understanding of population ancestry and selection in south-eastern African populations; and the data and results obtained will support research into the genetic contributions to infectious as well as non-communicable diseases in the region.     

An Overview of the Genetic Structure within the Italian Population from Genome-Wide Data

In spite of the common belief of Europe as reasonably homogeneous at genetic level, advances in high-throughput genotyping technology have resolved several gradients which define different geographical areas with good precision. When Northern and Southern European groups were considered separately, there were clear genetic distinctions. Intra-country genetic differences were also evident, especially in Finland and, to a lesser extent, within other European populations. Here, we present the first analysis using the 125,799 genome-wide Single Nucleotide Polymorphisms (SNPs) data of 1,014 Italians with wide geographical coverage. We showed by using Principal Component analysis and model-based individual ancestry analysis, that the current population of Sardinia can be clearly differentiated genetically from mainland Italy and Sicily, and that a certain degree of genetic differentiation is detectable within the current Italian peninsula population. Pair-wise F_ST statistics Northern and Southern Italy amounts approximately to 0.001 between, and around 0.002 between Northern Italy and Utah residents with Northern and Western European ancestry (CEU). The Italian population also revealed a fine genetic substructure underscoring by the genomic inflation (Sardinia vs. Northern Italy = 3.040 and Northern Italy vs. CEU = 1.427), warning against confounding effects of hidden relatedness and population substructure in association studies.     

Fine-scale north-to-south genetic admixture profile in Shaanxi Han Chinese revealed by genome-wide demographic history reconstruction

The Han Chinese are the world's largest ethnic group residing across China. Shaanxi province in northern China was a pastoral–agricultural interlacing region sensitive to climate change since Neolithic times, which makes it a vital place for studying population dynamics. However, genetic studies of Shaanxi Han are underrepresented due to the lack of high-density sampling and genome-wide data. Here, we genotyped 700 000 single nucleotide polymorphisms (SNPs) in 200 Han individuals from nine populations in Shaanxi and compared with available modern and ancient Eurasian individuals. We revealed a north–south genetic cline in Han Chinese with Shaanxi Han locating at the northern side of the cline. We detected the western Eurasian-related admixture in Shaanxi populations, especially in Guanzhong and Shanbei Han Chinese in proportions of 2%–4.6%. Shaanxi Han were suggested to derive a large part of ancestry (39%–69%) from a lineage that also contributed largely to ancient and present-day Tibetans (85%) as well as southern Han, supporting the common northern China origin of modern Sino-Tibetan-speaking populations and southwestward expansion of millet farmers from the middle-upper Yellow River Basin to the Tibetan Plateau and to southern China. The rest of the ancestry of Shaanxi Han was from a lineage closely related to ancient and present-day Austronesian and Tai-Kadai speaking populations in southern China and Southeast Asia. We also observed a genetic substructure in Shaanxi Han in terms of north–south-related ancestry corresponding well to the latitudes. Maternal mitochondrial DNA and paternal Y-chromosome lineages further demonstrated the aforementioned admixture pattern of Han Chinese in Shaanxi province.     

Extensive ethnolinguistic diversity at the crossroads of North China and South Siberia reflects multiple sources of genetic diversity

North China and South Siberia, populated by Altaic- and Sino-Tibetan-speaking populations, possess extensive ethnolinguistic diversity and serve as the crossroads for the initial peopling of America and western–eastern transcontinental communication. However, the population genetic structure and admixture history of northern East Asians remain poorly understood due to a lack of genome-wide data, especially for Mongolic-speaking people in China. We genotyped genome-wide single nucleotide polymorphisms for 510 individuals from 38 Mongolic, Tungusic, and Sinitic-speaking populations. We first explored the shared alleles and haplotypes within the studied groups. We then merged with 3508 published modern and ancient Eurasian individuals to reconstruct the deep evolutionary and natural selection history of northern East Asians. We identified genetic substructures within Altaic-speaking populations: Western Turkic people harbored more western Eurasian-related ancestry; Northern Mongolic people in Siberia and eastern Tungusic people in Amur River Basin (ARB) possessed a majority of Neolithic ARB related ancestry; Southern Mongolic people in China possessed apparent genetic influence from Neolithic Yellow River Basin (YRB) farmers. Additionally, we found the differentiated admixture history between western and eastern Mongolians and geographically close Northeast Hans: the former received a genetic impact from western Eurasians, and the latter retained the primary Neolithic YRB and ARB ancestry. Moreover, we demonstrated that Kalmyk people from the northern Caucasus Mountains possessed a strong genetic affinity with Neolithic Mongolian Plateau (MP) people, supporting the hypothesis of their eastern Eurasian origin and long-distance migration history. We also illuminated that historical pastoral empires in the MP contributed considerably to the gene pool of northern Mongolic people but rarely to the southern ones. We finally found natural selection signatures in Mongolians associated with alcohol metabolism. Our results demonstrated that the Neolithic ancestral sources from the MP or ARB played an important role in spreading Altaic populations and languages. The observed multisources of genetic diversity contributed significantly to the extensive ethnolinguistic diversity in northern East Asia.     

Structure of allozyme variation in Nordic Silene nutans (Caryophyllaceae): population size, geographical position and immigration history

We investigated allozyme variation in 34 populations of the perennial herb Silene nutans from Sweden and northern Finland, areas that were ice-covered during the last (Weichselian) glaciation. The present geographical structure of genetic variation in S. nutans in Sweden and northern Finland appears to have been mainly shaped by ancient historical processes. Patterns of variation in allele frequencies suggest two major postglacial immigration routes into Sweden, with populations entering the area from both the south and the east and forming a contact zone with admixed populations in central Sweden. While estimates of within-population genetic diversity and allelic richness are significantly correlated with present population size and geographical position (latitude), population size is not correlated with latitude. Low genetic diversity in the northern populations is more likely to have resulted from ancient stochastic events during the process of immigration than from recent population fragmentation. F _IS values are high and increase with latitude. Evidence of recent bottlenecks was detected in several southern Swedish populations: these can be interpreted in terms of population fragmentation as a result of anthropogenic disturbance. Soil pH is uncorrelated with population size and position.  © 2004 The Linnean Society of London, Biological Journal of the Linnean Society, 2004, 81 , 357–371.     

Range expansion and the possibility of an emerging contact zone between two subspecies of Chiffchaff Phylloscopus collybita ssp.

The Chiffchaff Phylloscopus collybita is represented in Sweden by two different subspecies; the northern well-established abietinus and the southern recently established collybita which has expanded its range northward during the past two decades. At present, an area approximately 500 km wide separates the two subspecies. In order to document differences between the northern and southern populations we compared morphology, vocalisation, habitat choice, and neutral genetic variation in mitochondrial (mt) DNA and at four microsatellite loci of 30 male Chiffchaffs from each subspecies. Our results show significant differences in several morphological traits and in song. Playback experiments revealed a significant difference in aggressive behaviour depending on which population-specific song that was played to the birds. Mitochondrial DNA was geographically structured with ∼90% of the birds carrying a mt haplotype matching their sample population. No allelic differences at the microsatellite loci were found between populations. Our data demonstrate a substantial differentiation between the northern and southern populations despite gene flow, clearly separating them into the subspecies abietinus and collybita .     

Contemporary genetic structure of Xantus's Hummingbird (Basilinna xantusii) in the Baja California peninsula

Genetic structure and phylogeographic patterns of natural populations are of great importance in assessing the conservation status of species. These population properties can be estimated using molecular markers of either mitochondrial DNA (mtDNA) or nuclear DNA (nDNA) to understand the historical, ecological and dispersal patterns that influence genetic exchange within and between populations. Basilinna xantusii is a sexually dimorphic hummingbird endemic to the Baja California Peninsula (BCP). It comprises three ancestral mitochondrial lineages linked to vicariant events, late Pleistocene climate changes and the geographical distribution of oases. This study aimed to determine and understand the current population genetic structure of this hummingbird. The genotypes of 16 microsatellite loci from 100 individuals collected across the geographical range of this species were compared with mtDNA sequences previously published. Cluster analyses identified five populations, two with almost no genetic admixture in the northern part of the BCP and three others with varying levels of admixed ancestry across the BCP. In San José de Magdalena, at the northernmost end of the range of Xantus's Hummingbird, 40% of individuals collected belong to one genetic cluster and the remaining 60% to another, both genetic clusters showing very little admixed ancestry. We hypothesize that, despite being in sympatry, these individuals do not interbreed, unlike the other populations where individuals showed ancestry coefficients of the other genetic groups. The philopatric behaviour of males and the long-range dispersal capacity of females probably determine the observed genetic differentiation pattern. The mito-nuclear discordance detected could be due to the molecular markers used and to female-biased dispersal. Gene flow is asymmetric in this species, being greater from north to south than vice versa, which is probably related to differences in the seasonality of precipitation across the BCP and to urbanization of the oases.     

Genetic population structure in horse chestnut (Aesculus hippocastanum L.): effects of human‐mediated expansion across Europe

The geographical distribution of existing populations of horse chestnut (Aesculus hippocastanum L.) in Europe is determined by past demographic events during the Quaternary. In the present study we evaluate the imprints that northward expansions originated from common ancestry at southern Europe may have left on the present patterns of genetic variation for horse chestnut across the continent. Genetic diversity and levels of population structure in a European south–north gradient, ranging from the Balkans to the Scandinavian Peninsula, were determined with Amplified Fragment Length Polymorphism (AFLP) markers in 159 loci. A family of rarefaction techniques for the estimation of gene diversity was used to exclude potential confounding effects as a result of the unequal sample sizes. The results indicate that northern populations are not more genetically depleted than southern populations, thus suggesting that diversity for this species is not correlated with latitudinal distribution. Detailed hypotheses based on prediction models for different historical events associated with human‐mediated spread of cultivation are examined for a better understanding of the current genetic patterns of regional differentiation.     

Geographic structure of genetic variation in the widespread woodland grass Milium effusum L. A comparison between two regions with contrasting history and geomorphology.

Allozyme variation in the forest grass Milium effusum L. was studied in 21-23 populations within each of two equally sized densely sampled areas in northern and southern Sweden. In addition, 25 populations from other parts of Eurasia were studied for comparison. The structure of variation was analysed with both diversity statistics and measures based on allelic richness at a standardised sample size. The species was found to be highly variable, but no clear geographic patterns in the distribution of alleles or in overall genetic differentiation were found, either within the two regions or within the whole sample. Thus, no inferences about the direction of postglacial migration could be made. Obviously, migration and gene flow must have taken place in a manner capable of randomising the distribution of alleles. However, there were clear differences in levels and structuring of the variation between the two regions. Levels of variation, both in terms of genetic diversity and allelic richness, were lower in northern Sweden as compared with southern Sweden. In contrast, different measures of geographic structure all showed higher levels of population differentiation in the northern region. This is interpreted as due to different geomorphological conditions in the two regions, creating a relatively continuous habitat and gene flow in the southern region as compared with the northern region where the species, although common, is confined to narrow and mutually isolated corridors in the landscape.     

Recent evolutionary history of the bluethroat (Luscinia svecica) across Eurasia

We analysed mitochondrial DNA (mtDNA) sequences from 154 bluethroats (Luscinia svecica) sampled at 21 sites throughout much of their Eurasian range. A previously reported, single base-pair mtDNA difference between L. s. svecica and L. s. namnetum was inconsistent upon expanded geographical sampling. A significant FST value (0.29) and an isolation-by-distance effect show the existence of geographical differentiation. Phylogenetic analysis of haplotypes revealed northern and southern groups, although lineage sorting is incomplete. There was no geographical structure to the haplotype tree within groups, and currently recognized subspecies were not supported. A minimum evolution tree based on pairwise mtDNA genetic distances among average samples showed the same two broadly distributed northern and southern groups. These groups abut in the centre of the latitudinal range, and were possibly isolated by forest that developed and spread westward over the last 15 000 years. Pairwise FST values averaged 0.16 in the southern group, 0.04 in the northern group, and 0.42 between groups. Mismatch distributions suggested population growth in each group, with that in the south being more recent. In the northern group, the geographical pattern in tau suggested northward and eastward expansion. Analysis of nucleotide diversity suggested westward expansion in the southern group. The northern group had higher nucleotide diversity than the southern group, consistent with a larger current population size in the north. Given the significant FST, incompletely sorted haplotype tree, and broadly patterned minimum evolution tree, L. svevica appears to represent a species at an intermediate stage of differentiation between panmixia and reciprocal monophyly.     

Latitudinal divergence in a widespread amphibian: Contrasting patterns of neutral and adaptive genomic variation

Stochastic effects from demographic processes and selection are expected to shape the distribution of genetic variation in spatially heterogeneous environments. As the amount of genetic variation is central for long‐term persistence of populations, understanding how these processes affect variation over large‐scale geographical gradients is pivotal. We investigated the distribution of neutral and putatively adaptive genetic variation, and reconstructed demographic history in the moor frog (Rana arvalis) using 136 individuals from 15 populations along a 1,700‐km latitudinal gradient from northern Germany to northern Sweden. Using double digest restriction‐site associated DNA sequencing we obtained 27,590 single nucleotide polymorphisms (SNPs), and identified differentiation outliers and SNPs associated with growing season length. The populations grouped into a southern and a northern cluster, representing two phylogeographical lineages from different post‐glacial colonization routes. Hybrid index estimation and demographic model selection showed strong support for a southern and northern lineage and evidence of gene flow between regions located on each side of a contact zone. However, patterns of past gene flow over the contact zone differed between neutral and putatively adaptive SNPs. While neutral nucleotide diversity was higher along the southern than the northern part of the gradient, nucleotide diversity in differentiation outliers showed the opposite pattern, suggesting differences in the relative strength of selection and drift along the gradient. Variation associated with growing season length decreased with latitude along the southern part of the gradient, but not along the northern part where variation was lower, suggesting stronger climate‐mediated selection in the north. Outlier SNPs included loci involved in immunity and developmental processes.     

Phylogenetic divergence in leatherside chub (Gila copei) inferred from mitochondrial cytochrome b sequences

We examined intra-specific phylogenetic relationships in leatherside chub, Gila copei. The complete mitochondrial (mt) cytochrome b gene (1140 bp) was sequenced for 30 individuals from 10 populations that span the geographical distribution of this species. Traditional phylogenetic analyses revealed two deeply divergent and evolutionarily distinct mtDNA clades that are geographically separated in northern and southern drainage basins. Interpopulation sequence variation between clades ranged from 7.7 to 8.1%. The northern clade was genetically more similar and phylogenetically more closely related to the selected out-group Lepidomeda m. mollispinus than to the southern clade, suggesting that the taxonomy of this species may require revision. Sequence variation among populations within clades ranged from 0 to 0.3% in the north and from 0 to 0.7% in the south. Statistical parsimony was used to construct phylogenetic networks of haplotypes within clades. Nested clade analysis revealed that geographical fragmentation has played an important role in genetic structuring within northern and southern clades.