Population genetics analysis of the origin of the Oriental fruit fly, Bactrocera dorsalis Hendel (Diptera: Tephritidae), in northern Yunnan Province, China

We examined genetic variation in the Oriental fruit fly, Bactrocera dorsalis (Hendel), using six populations in two regions of Yunnan Province, China, to determine the distribution and likely mechanism for the dispersal of this fly. A 501‐bp portion of the mitochondrial cytochrome oxidase gene was sequenced from a minimum of eight individuals from each population, and 43 haplotypes were observed in the six Bactrocera dorsalis populations. When comparing the genetic diversity of populations in the northern and southern regions, which differ with respect to elevation, climate and plant phenology, we found a significantly greater haplotype diversity in the southern region (permutation test; P < 0.05), suggesting that the northern populations, those at Kunming and Qujing, probably originated from somewhere in the southern region. F_ST and number of pairwise differences revealed a high level of differentiation between the Panxi population and the other populations (permutation test; P < 0.05). Although the difference was marginally insignificant, the Shuitang population seemed to have differentiated from both northern populations. The Mantel test did not detect any isolation due to geographic distance. An amova analysis found that 2.56% of the variance was caused by the Panxi population. Haplotype network analysis showed that none of the six populations had a specific genetic lineage. Together, these analyses suggest that long‐distance dispersal has occurred for this species, and the species most probably took advantage of both a mountain pass and prevailing air currents. The Panxi population was significantly isolated from the others, probably because of its distinguishing habitat features, host plants or the recent reduction of the population size.     

Genetic diversity and structure of natural populations of Plathymenia reticulata (Mimosoideae), a tropical tree from the Brazilian Cerrado

Plathymenia reticulata is a tropical tree native to the Brazilian Cerrado, one of the most important and endangered ecosystems in Brazil. This species presents high-quality wood and potential for recovery of degraded areas. Despite its importance, almost nothing is known about its genetic or ecological features. Random amplified polymorphic DNA (RAPD) markers were used to investigate the genetic diversity and structure of six natural populations of P. reticulata. DNAs from 117 adult individuals were amplified with 10 random primers and Shannon's index and amova were used to evaluate the levels of genetic diversity within and among populations. Through 72 markers, 70.8% of which were polymorphic, it was possible to obtain 117 unique RAPD phenotypes. The levels of genetic variability found in the six populations of P. reticulata were considerable and most of the genetic variation was found between individuals within populations, although pairwise PH(ST) values indicated significant divergence between populations. The among-population component accounted for, respectively, 12.3% and 16% of the genetic variation, according to amova and Shannon's index. These results were compared with other genetic studies on plant species and such a level of differentiation among populations corresponds to that which has usually been observed for outcrossing plants. The importance of maintenance of the P. reticulata populations and implications of the analysis of adult individuals, considering the longevity of this species and the relatively recent Cerrado fragmentation, are discussed.     

Reconstruction of the Sudden Oak Death epidemic in California through microsatellite analysis of the pathogen Phytophthora ramorum

The genetic structure of the clonally reproducing Sudden Oak Death (SOD) pathogen in California was investigated using seven variable microsatellites. A total of 35 multilocus genotypes were identified among 292 samples representative of populations from 14 forest sites and of the nursery trade. amova indicated significant genetic variability both within (44.34%) and among populations (55.66%). Spatial autocorrelation analyses indicated that Moran's index of similarity reached a minimum of 0.1 at 350 m, increased to 0.4 at 1500 m and then decreased to zero at 10 km. These results suggest a bimodal pattern of spread, with medium range dispersal (1500–10 000 m) putatively attributed to the presence of strong winds. Lack of genetic structure was identified for three groups of populations. One group notably included the nurseries' population and two forest populations, both linked to early reports of the pathogen. A neighbour-joining analysis based on pairwise Φ_ST values indicated that the clade inclusive of the nurseries' populations is basal to all California populations. A network analysis identified three common genotypes as the likely founders of the California infestation and proposes a stepwise model for local evolution of novel genotypes. This was supported by the identification in the same locations of novel genotypes and of their 1- or 2-step parents. We hypothesize that the few undifferentiated population groups indicate historical human spread of the pathogen, while the general presence of genetically structured populations indicates that new infestations are currently generated by rare medium or long-range natural movement of the pathogen, followed by local generation of new genotypes.     

High gene flow across large geographic scales reduces extinction risk for a highly specialised coral feeding butterflyfish

The vulnerability of ecologically specialised species to environmental fluctuations has been well documented. However, population genetic structure can influence vulnerability to environmental change and recent studies have indicated that specialised species may have lower genetic diversity and greater population structuring compared to their generalist counterparts. To examine whether there were differences in population genetic structure between a dietary specialist (Chaetodon trifascialis) and a dietary generalist (Chaetodon lunulatus) we compared the demographic history and levels of gene flow of two related coral-feeding butterflyfishes. Using allele frequencies of ≥11 microsatellite loci and >350 bases of mitochondrial control region sequence our analyses of C. trifascialis and C. lunulatus from five locations across the Pacific Ocean revealed contrasting demographic histories and levels of genetic structure. Heterozygosity excess tests, neutrality tests and mismatch distributions were all highly significant in the dietary specialist C. trifascialis (all P < 0.01), suggesting genetic bottlenecks have occurred in all locations. In contrast, we found little evidence of genetic bottlenecks for the dietary generalist C. lunulatus. High gene flow and low genetic structuring was detected among locations for C. trifascialis (amova: R(ST) = 0.0027, P = 0.371; Φ(ST) = 0.068, P < 0.0001). Contrary to our expectations, a greater level of genetic structuring between locations was detected for C. lunulatus (amova: R(ST) = 0.0277, Φ(ST) = 0.166, both P < 0.0001). These results suggest that dietary specialisation may affect demographic history through reductions in population size following resource declines, without affecting population structure through reductions in gene flow in the same way that habitat specialisation appears to. Although C. trifascialis is highly vulnerable to coral loss, the high gene flow detected here suggests populations will be able to recover from local declines through the migration of individuals.     

Gene flow, dispersal, and nested clade analysis among populations of the stonefly Peltoperla tarteri in the southern Appalachians

We examined the genetic structure and phylogeography of populations of the stonefly Peltoperla tarteri in the Southern Appalachians to determine the extent and likely mechanism for dispersal of this stream insect. A 454-base-pair (bp) portion of the mitochondrial control region was sequenced from a minimum of 20 individuals from eight populations. Pairwise FST and exact tests showed high levels of differentiation among almost all populations except those on the same stream. amova analysis detected significant genetic differentiation between streams within drainages (phi(SD) = 0.14, P < 0.001), and there was a slight positive correlation between aquatic distance and genetic distance (r = 0.295, P = 0.03). According to nested clade analysis, the present day pattern of genetic variation in P. tarteri is the result of a historical range expansion coupled with restricted gene flow with isolation by distance. Together, these analyses suggest that adult dispersal is limited and that movement by larvae is the primary dispersal mechanism for P. tarteri.     

A Monte Carlo Permutation Test for Random Mating Using Genome Sequences

Testing for random mating of a population is important in population genetics, because deviations from randomness of mating may indicate inbreeding, population stratification, natural selection, or sampling bias. However, current methods use only observed numbers of genotypes and alleles, and do not take advantage of the fact that the advent of sequencing technology provides an opportunity to investigate this topic in unprecedented detail. To address this opportunity, a novel statistical test for random mating is required in population genomics studies for which large sequencing datasets are generally available. Here, we propose a Monte-Carlo-based-permutation test (MCP) as an approach to detect random mating. Computer simulations used to evaluate the performance of the permutation test indicate that its type I error is well controlled and that its statistical power is greater than that of the commonly used chi-square test (CHI). Our simulation study shows the power of our test is greater for datasets characterized by lower levels of migration between subpopulations. In addition, test power increases with increasing recombination rate, sample size, and divergence time of subpopulations. For populations exhibiting limited migration and having average levels of population divergence, the statistical power approaches 1 for sequences longer than 1Mbp and for samples of 400 individuals or more. Taken together, our results suggest that our permutation test is a valuable tool to detect random mating of populations, especially in population genomics studies.     

Demographic population structure of black howler monkeys in fragmented and continuous forest in Chiapas,Mexico: Implications for conservation

For wild primates, demography studies are increasingly recognized as necessary for assessing the viability of vulnerable populations experiencing rapid environmental change. In particular, anthropogenic changes such as habitat loss and fragmentation can cause ecological and behavioral changes in small, isolated populations, which may, over time, alter population density and demographic structure (age/sex classes and group composition) in fragment populations relative to continuous forest populations. We compared our study population of Endangered black howler monkeys (Alouatta pigra) in 34 forest fragments around Palenque National Park (PNP), Mexico (62 groups, 407 individuals), to the adjacent population in PNP, protected primary forest (21 groups, 134 individuals), and to previous research on black howlers in fragments in our study area (18 groups, 115 individuals). We used χ² and Mann–Whitney U tests to address the questions: (a) what is the current black howler demographic population structure in unprotected forest fragments around PNP? (b) How does it compare to PNP's stable, continuous population? (c) How has it changed over time? Compared to the PNP population, the fragment populations showed higher density, a significantly lower proportion of multimale groups, and significantly fewer adult males per group. The population's age/sex structure in the fragmented landscape has been stable over the last 17 years, but differed in a higher proportion of multifemale groups, higher density, and higher patch occupancy in the present. In the context of conservation, some of our results may be positive as they indicate possible population growth over time. However, long-term scarcity of adult males in fragments and associated effects on population demographic structure might be cause for concern, in that it may affect gene flow and genetic diversity. The scarcity of adult males might stem from males experiencing increased mortality while dispersing in the fragmented landscape, whereas females might be becoming more philopatric in fragments.     

Genetic variation and population differentiation in the lichen-forming ascomycete Xanthoria parietina on the island Storfosna, central Norway

Genetic diversity and fine-scale population structure in the lichen-forming ascomycete Xanthoria parietina was investigated using sequence variation in part of the intergenic spacer (IGS) and the complete internal transcribed spacer (ITS) regions of the nuclear ribosomal DNA. Sampling included 213 and 225 individuals, respectively, from seven populations in two different habitats, bark and rock, on the island Storfosna off the central west coast of Norway. Both markers revealed significant variation and a total of 10 IGS and 16 ITS haplotypes were found. There were no signs of significant positive spatial autocorrelation at any spatial size class down to 10% of transect length, nor did we find significant deviations from neutrality or signs of historical population expansion. Analysis of molecular variance (amova) indicated that most of the genetic variance observed was within populations, but when populations were grouped according to habitat, more than a quarter of the variance was explained among groups. Pairwise comparisons of populations (F(ST), exact tests of population differentiation) revealed significant differentiation between populations in different habitats (on bark or rock), but not between populations in the same habitat. Haplotype networks show that internal and presumably old haplotypes are shared between habitats, whereas terminal haplotypes tend to be unique to a habitat, mostly bark. We interpret the observed pattern to mean that there is no evidence of restricted gene flow between populations in the same habitat at the present spatial scale (interpopulation distances one or a few kilometres). On the other hand, differentiation between habitats is considerable, which we attribute to restricted gene flow between habitats (habitat isolation). Evidence suggests that the observed differentiation did not evolve locally. Estimates of divergence time between populations in the respective habitats indicate that an ancestral population started to diverge at least 34,000 years ago but probably much further back in time.     

Genetic diversity and differentiation of guanaco populations from Argentina inferred from microsatellite data

Genotype data from 14 microsatellite markers were used to assess the genetic diversity and differentiation of four guanaco populations from Argentine Patagonia. These animals were recently captured in the wild and maintained in semi-captivity for fibre production. Considerable genetic diversity in these populations was suggested by the finding of a total of 162 alleles, an average mean number of alleles per locus ranging from 6.50 to 8.19, and H(e) values ranging from 0.66 to 0.74. Assessment of population differentiation showed moderate but significant values of F(ST)=0.071 (P=0.000) and R(ST)=0.083 (P=0.000). An amova test showed that the genetic variation among populations was 5.6% while within populations it was 94.4%. A number of 6.6 migrants per generation may support these results. Unambiguous individual assignment to original populations was obtained for the Pilcaniyeu, Las Heras and La Esperanza populations. The erroneous assignment of 18.75% Rio Mayo individuals to the Las Heras population can be explained by the low genetic differentiation found between these two populations. Thirty-nine of 56 loci per population combinations were in Hardy--Weinberg disequilibrium because of guanaco heterozygote deficiency, which may be explained by population subdivision. The high level of genetic diversity of the guanacos analysed here indicates that the Patagonian guanaco constitutes an important genetic resource for conservation or economic utilization programmes.     

种群微卫星DNA分析中样本量对各种遗传多样性度量指标的影响

我们以中国飞蝗种群的微卫星遗传分析数据为例 ,评估了取样对种群遗传多样性指标的影响 ,结果显示 :样本大小与所观测到的每位点等位基因数、平均等位基因数及基因丰富度指数均呈显著正相关 ,而与期望杂合度无显著相关 ;微卫星位点多态性的高低直接影响所观测到的种群基因丰富度及其检测所需的样本量 ;对大多数种群遗传和分子生态学研究而言 ,30 - 5 0个个体是微卫星DNA分析所需要的最小样本量。基因丰富度经过稀疏法或多次随机抽样法校正后 ,可适用于瓶颈效应等种群历史数量变动的检测。另外 ,在研究中 ,还应避免采集时间的不同及样本的性比构成所可能造成的对种群遗传结构的影响     

Genetic differentiation in relation to marine landscape in a broadcast-spawning bivalve mollusc (Placopecten magellanicus)

Marine bivalves are sessile or sedentary as adults but have planktonic larvae which can potentially disperse over large distances. Consequently larval transport is expected to play a prominent role in facilitating gene flow and determining population structure. The sea scallop (Placopecten magellanicus) is a dioecious species with high fecundity, broadcast spawning and a c. 30-day planktonic larval stage, yet it forms discrete populations or 'beds' which have significantly different dynamics and characteristics. We analysed variation at six microsatellite loci in 12 locations throughout the geographic range of the species from Newfoundland, Canada, to New Jersey, USA. Significant differentiation was present and the maximum pairwise theta value, between one of the Newfoundland samples in the north and a sample from the southern portion of the range, was high at 0.061. Other proximate pairs of samples had no detectable genetic differentiation. Mantel tests indicated a significant isolation by distance, but only when one of the populations was excluded. A landscape genetic approach was used to detect areas of low gene flow using a joint analysis of spatial and genetic information. The two major putative barriers inferred by Monmonier's algorithm were then used to define regions for an analysis of molecular variance (amova). That analysis showed a significant but low percentage (1.2%) of the variation to be partitioned among regions, negligible variation among populations within regions, and the majority of the variance distributed between individuals within populations. Prominent currents were concordant with the demarcation of the regions, while a novel approach of using particle tracking software to mimic scallop larval dispersal was employed to interpret within-region genetic patterns.     

Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study

The identification of genetically homogeneous groups of individuals is a long standing issue in population genetics. A recent Bayesian algorithm implemented in the software STRUCTURE allows the identification of such groups. However, the ability of this algorithm to detect the true number of clusters (K) in a sample of individuals when patterns of dispersal among populations are not homogeneous has not been tested. The goal of this study is to carry out such tests, using various dispersal scenarios from data generated with an individual-based model. We found that in most cases the estimated 'log probability of data' does not provide a correct estimation of the number of clusters, K. However, using an ad hoc statistic DeltaK based on the rate of change in the log probability of data between successive K values, we found that STRUCTURE accurately detects the uppermost hierarchical level of structure for the scenarios we tested. As might be expected, the results are sensitive to the type of genetic marker used (AFLP vs. microsatellite), the number of loci scored, the number of populations sampled, and the number of individuals typed in each sample.     

DNA from bird-dispersed seed and wind-disseminated pollen provides insights into postglacial colonization and population genetic structure of whitebark pine (Pinus albicaulis)

Uniparentally inherited mitochondrial (mt)DNA and chloroplast (cp)DNA microsatellites (cpSSRs) were used to examine population genetic structure and biogeographic patterns of bird-dispersed seed and wind-disseminated pollen of whitebark pine (Pinus albicaulis Engelm.). Sampling was conducted from 41 populations throughout the range of the species. Analyses provide evidence for an ancestral haplotype and two derived mtDNA haplotypes with distinct regional distributions. An abrupt contact zone between mtDNA haplotypes in the Cascade Range suggests postglacial biogeographic movements. Among three cpSSR loci, 42 haplotypes were detected within 28 cpSSR sample populations that were aggregated into six regions. Analysis of molecular variance (amova) was used to determine the hierarchical genetic structure of cpSSRs. amova and population pairwise comparisons (FST ) of cpSSR, and geographical distribution of mtDNA haplotypes provide insights into historical changes in biogeography. The genetic data suggest that whitebark pine has been intimately tied to climatic change and associated glaciation, which has led to range movements facilitated by seed dispersal by Clark's nutcracker (Nucifraga columbiana Wilson). The two hypotheses proposed to explain the genetic structure are: (i) a northward expansion into Canada and the northern Cascades in the early Holocene; and (ii) historical gene flow between Idaho and the Oregon Cascades when more continuous habitat existed in Central Oregon during the late Pleistocene. Genetic structure and insights gained from historical seed movements provide a basis on which to develop recovery plans for a species that is at risk from multiple threats.     

Mitochondrial DNA inherited variants are associated with successful aging and longevity in humans.

Mitochondrial DNA (mtDNA) is characterized by high variability, maternal inheritance, and absence of recombination. Studies of human populations have revealed ancestral associated polymorphisms whose combination defines groups of mtDNA types (haplogroups) that are currently used to reconstruct human evolution lineages. We used such inherited mtDNA markers to compare mtDNA population pools between a sample of individuals selected for successful aging and longevity (212 subjects older than 100 years and in good clinical condition) and a sample of 275 younger individuals (median age 38 years) carefully matched as to sex and geographic origin (northern and southern Italy). All nine haplogroups that are typical of Europeans were found in both samples, but male centenarians emerged in northern Italy as a particular sample: 1) mtDNA haplogroup frequency distribution was different between centenarians and younger individuals (P=0.017 by permutation tests); and 2) the frequency of the J haplogroup was notably higher in centenarians than in younger individuals (P=0.0052 by Fisher exact test). Since haplogroups are defined on the basis of inherited variants, these data show that mtDNA inherited variability could play a role in successful aging and longevity.     

Microsatellite analysis reveals interpopulation differentiation and gene flow in the endangered tree Changiostyrax dolichocarpa (Styracaceae) with fragmented distribution in central China

Polymorphic simple sequence repeat (SSR) markers were used to investigate the impact of habitat fragmentation on the population structure and gene flow of Changiostyrax dolichocarpa, a critically endangered tree in central China. Intrapopulation genetic diversity, population structure and gene flow in the five extant populations of this species were analysed by eight SSR markers. Intrapopulation genetic diversity results suggest that C. dolichocarpa remnants maintained a relatively high degree of genetic diversity despite severe fragmentation. Low genetic differentiation among populations was found based on Wright's F(ST) and amova analysis. Both the F(ST)-based estimate and private allele method revealed high historical gene flow among the remnant populations. Recent immigrants, detected by assignment tests, tend to decrease from the grandparent generation to the current generation. The potentially highly restricted current gene flow among fragments may render the fragmented populations of C. dolichocarpa at a higher risk of local extinction several generations after fragmentation. Both in situ and ex situ conservation management for the remnant populations of C. dolichocarpa are therefore urgently needed to rescue remaining genetic diversity.     

Hierarchical structure of the Sicilian goats revealed by Bayesian analyses of microsatellite information

Genetic structure and relationship amongst the main goat populations in Sicily (Girgentana, Derivata di Siria, Maltese and Messinese) were analysed using information from 19 microsatellite markers genotyped on 173 individuals. A posterior Bayesian approach implemented in the program STRUCTURE revealed a hierarchical structure with two clusters at the first level (Girgentana vs. Messinese, Derivata di Siria and Maltese), explaining 4.8% of variation (amova Ф_ST estimate). Seven clusters nested within these first two clusters (further differentiations of Girgentana, Derivata di Siria and Maltese), explaining 8.5% of variation (amova Ф_SC estimate). The analyses and methods applied in this study indicate their power to detect subtle population structure.     

Microsatellite markers to assess the influence of population size, isolation and demographic change on the genetic structure of the UK butterfly Polyommatus bellargus

Five microsatellite DNA markers were isolated and used to quantify population genetic structure among a subset of UK populations of the Adonis blue (Polyommatus bellargus Rottemburg). Specifically, whether population size, degree of isolation or history of bottlenecking in 1976-1978 can explain current patterns of genetic variation. The butterfly is at its northern range limit in the UK, where it exists as a highly fragmented metapopulation on isolated pockets of calcareous grassland. Most populations were affected by a severe bottleneck in the late 1970s, when a drought caused the host plant (Hippocrepis comosa) to wilt. Mantel tests and spatial autocorrelation analysis indicated a significant effect of isolation by distance among the UK populations, a relationship that broke down at greater geographical scales (> 23.85 km), probably because of large areas of unsuitable habitat presenting barriers to gene flow. Similarly, amova revealed that variation among geographical regions was almost double that observed within regions. Larger populations were found to support significantly higher levels of genetic diversity, suggesting that small populations may lose genetic diversity through drift. If, as in other butterfly species, low genetic diversity increases the probability of population extinction, then these populations are likely to be under threat. Neither isolation nor a history of bottlenecks appeared to influence genetic diversity. The results indicate that adequate population size a crucial factor in the conservation of genetic diversity in P. bellargus in the UK.     

Post-glacial range fragmentation is responsible for the current distribution of Potentilla matsumurae Th. Wolf (Rosaceae) in the Japanese archipelago

Aim We aimed to elucidate how the current geographic distribution of alpine plants in the Japanese archipelago was shaped during Quaternary climatic oscillations, using Potentilla matsumurae as a case study. According to previous phylogeographic studies, post‐glacial range fragmentation (vicariance scenario) and stepwise migration (dispersal scenario) are both possible. We thus aimed to assess which scenario is more probable for the distribution changes of alpine plants in the Japanese archipelago. Location The alpine zone in the Japanese archipelago. Methods Using amplified fragment length polymorphism we determined the genotype of 161 individuals of P. matsumurae from 22 populations. Relationships among individuals and populations were examined using principal coordinates analysis and a neighbour‐joining (NJ) tree, respectively. To examine the genetic population structure, we performed analysis of molecular variance (amova ) and structure analysis. Results Differentiation between central Honshu and northern Japan was not very strong based on the principal coordinates analysis among individuals, the NJ tree of populations (59% bootstrap support), or amova (12% of genetic variation). Moreover, structure analysis did not detect clear geographic differentiation across populations. Although the populations in central Honshu were structured geographically (Mantel test: r = 0.45, P < 0.005; NJ tree), those in northern Japan did not exhibit geographic structure regardless of geographic distance (Mantel test: r = 0.26, P = 0.03; NJ tree). Population relationships in the NJ tree did not always reflect the geographic location. Main conclusions The current geographic structure of P. matsumurae could not be explained by stepwise migration. This suggests that a single continuous distribution during the last glacial period was later fragmented, perhaps by recovering forest, during the post‐glacial period, resulting in the current distribution and phylogeographic structure of P. matsumurae. Our data support the vicariance scenario.     

Ancestry and divergence of subtropical montane forest isolates: molecular biogeography of the genus Abies (Pinaceae) in southern México and Guatemala

The genus Abies has a complex history in southern México and Guatemala. In this region, four closely related species, Abies flinckii , A. guatemalensis , A. hickelii , and A. religiosa , are distributed in fragmented and isolated montane populations. Range-wide genetic variation was investigated across species using cytoplasmic DNA markers with contrasted inheritance. Variation at two maternally inherited mitochondrial DNA markers was low. All species shared two of the nine mitotypes detected, while the remaining seven mitochondrial DNA types were restricted to a few isolated stands. Mitochondrial genetic differentiation across taxa was high ( G _ST = 0.933), it was not related to the taxonomic identity ( amova ; P  > 0.05) of the populations, and it was not phylogeographically structured ( G _ST ≈  N _ST). In contrast, variation at three paternally inherited chloroplast DNA microsatellites was high. Chloroplast genetic differentiation was lower ( G _ST = 0.402; R _ST = 0.547) than for mitochondrial DNA, but it was significantly related to taxonomy ( amova ; P  < 0.001), and exhibited a significant phylogeographical structure ( G _ST <  R _ST). Different analyses of population structure indicated that A. flinckii was the most divergent taxon, while the remaining three species formed a relatively homogeneous group. However, a small number of the populations of these three taxa, all located at the limits of their respective ranges or in the Transverse Volcanic Belt, diverged from this main cluster. These trends suggest that the Mesoamerican Abies share a recent common ancestor and that their divergence and speciation is mainly driven by genetic drift and isolation during the warm interglacial periods.     

Association analysis of the RGS4 gene in Han Chinese and Scottish populations with schizophrenia

We investigated the RGS4 as a susceptibility gene for schizophrenia in Chinese Han (184 trios and 138 sibling pairs, a total of 322 families) and Scottish (580 cases and 620 controls) populations using both a family trio and case-control design. Both the samples had statistical power greater than 70% to detect a heterozygote genotype relative risk of >1.2 for frequent RGS4-risk alleles. We genotyped four single nucleotide polymorphisms (SNPs) which have previously been associated with schizophrenia as either individually or part of haplotypes. Allele frequencies and linkage disequilibrium between the SNPs was similar in the two populations. In the Chinese sample, no individual SNPs or any of their haplotypes were associated with schizophrenia. In the Scottish population, one SNP (SNP7) was significantly over-represented in the cases compared with the controls (0.44 vs. 0.38; A allele; chi(2) 7.08, P = 0.011 after correction for correlation between markers by permutation testing). One two-marker haplotype, composed of alleles T and A of SNP4 and SNP7, respectively, showed individual significance after correction by permutation testing (chi(2) 6.8; P = 0.04). None of the full four-marker haplotypes showed association, including the G-G-G-G haplotype previously associated with schizophrenia in more than one sample and the A-T-A-A haplotype. Thus, our data do not directly replicate previous associations of RGS4, but association with SNP 7 in the Scottish population provides some support for a role in schizophrenia susceptibility. We cannot conclusively exclude RGS4, as associated haplotypes are likely to be surrogates for unknown causative alleles, whose relationship with overlying haplotypes may differ between the population groups. Differences in the association seen across the two populations could result from methodological factors such as diagnostic differences but most likely result from ethnic differences in haplotype structures within RGS4.