Population Genetic Structure of Transcaucasian Mole Vole (<Emphasis Type="Italic">Ellobius lutescens</Emphasis>) Along Zagros Mountains,Iran

The Zagros and Alborz mountainous ridges can be regarded as one of the most interesting and known physical barriers responsible for the vicariance event. Based on the probable effect of Zagros Mountains on the rodent population vicariance, a research on Caucasian mole vole phylogeography, population genetic structure and diversity was designed along the mentioned mountainous areas. To this end, a total of 38 tissue samples were collected from the northern parts of the study area to the southern parts. Obtained mitochondrial cytb (1041 bp) sequences were used in this phylogenetic analysis. The phylogenetic analysis was based on the TRN+I evolutionary model and gaining Bayesian phylogenetic tree with maximum verification. By using median joining logic, the relationships between different acquired haplotypes were analyzed. It was shown that the Caucasian mole vole population had been disjointed (based on posterior probability of 1 and 100 bootstraps) along the Zagros mountainous ridges, especially in both geographical extremes located in the northern and southern parts of the mountainous ridges. Meanwhile from the 38 analyzed sequences, 17 haplotypes were obtained, of which 10 haplotypes were unique. The mutational steps between haplotypes were assessed by generating statistical parsimony haplotype networks, which yielded 36 mutational steps between the northern and southern populations. Based on neutrality tests and analyzing their power under sudden population expansions, it was found that this event happened around the northern and southern populations. Genetic distance of two percent between the northern and southern populations indicated the existence of local adaptations by these two groups, which can be regarded as evolutionary units.     

An ancient mitochondrial polymorphism in Adalis bipunctata linked to a sex-ratio-distorting bacterium

Sex-ratio-distorting microbes are common parasites of arthropods. Although the reasons they have invaded and spread though populations are well understood, their subsequent dynamics within those populations are virtually unknown. We have found that different strains of a male-killing Rickettsia bacterium infecting the beetle Adalia bipunctata are associated with distinct mitochondrial haplotypes, which is expected as both the mitochondria and the bacteria are maternally transmitted. These mitochondrial haplotypes shared a common ancestor >2 million years ago, and their overall diversity is significantly greater than expected under neutrality from comparisons with a nuclear gene. Furthermore, a variety of statistical tests show strong deviations from neutrality in mitochondrial but not in nuclear genes. We therefore conclude that natural selection is probably maintaining a polymorphism of different Rickettsia strains in this species. Despite the age of the different mitochondrial haplotypes, there is very little genetic diversity within them. Furthermore, there is considerable variation in mitochondrial haplotype and bacterial strain frequency between populations, despite it being thought that this species has fairly low levels of population structure. We conclude that the fitness of these male killers may be negatively frequency dependent or different strains may be favored in different populations. These hypotheses await experimental confirmation.     

Demographic expansion of parasitic nematodes of livestock based on mitochondrial DNA regions that conflict with the infinite-sites model

Mitochondrial ND4 sequences of populations of four species of parasitic nematodes of livestock were subjected to demographic analyses. Deviation from selective neutrality was detectable using the frequency spectrum of segregating sites and highly negative neutrality statistics. However, the mitochondrial data sets do not comply with the infinite-sites model that underlies these tests, and as a consequence, it was not established whether these features are solely a result of population expansion, or whether aspects of the molecular evolution of these mitochondrial regions are also involved. Coalescent analyses based on Fu's Fs neutrality test, which incorporated estimates of rate heterogeneity, the transition-transversion ratio and nucleotide bias, as well as analyses that are fairly robust to deviations from the infinite-sites model supported population expansion. Also analyses that do not depend on the infinite-sites model suggested historical population expansion of these nematodes. The very similar time since expansion, the absence of signatures of positive selection in ND4 and the logical association with human demography imply that selective sweeps of mitochondrial variants are less probable, and that expansion is the most likely scenario for the parasitic nematodes of livestock. The methods used to characterize the expansion have different assumptions and emphasize different aspects of expansions. The resulting restrictions on the interpretation of expansions are outlined.     

Identification and characterization of t haplotypes in wild mice populations using molecular markers

As part of a population genetics survey of the hybrid zone between mouse subspecies Mus musculus domesticus and M. m. musculus, we identified and characterized the t haplotypes in 1068 mice from 186 different populations in a 2500 km2 area in central Jutland. On the basis of two t-specific PCR markers, 130 mice possessed this haplotype. The allele frequencies at six microsatellites on the third and fourth chromosomal inversions of the t region were sufficiently different between t-bearing and non-t-bearing mice, and linkage disequilibria sufficiently marked on the t haplotype, to be able to reconstitute the genotype of most t haplotypes. A total of three frequent and 15 rarer haplotypes were identified. These haplotypes resemble each other more than they resemble a panel of known haplotypes from a wide range of geographical regions, except for tw73, which was also extracted from Jutland. The patterns of variation at the microsatellite loci suggest that the Jutland haplotypes were derived from a small number of haplotypes, followed by recombination between complementing haplotypes. Further evidence of recombination came from complementation tests that we performed, showing the lack of concordance between the degrees of complementation and of molecular resemblance between haplotypes. This study shows that it is possible to characterize the presence and variation of t haplotypes by a population genetics approach using simple molecular markers. However recombination between t haplotypes has occurred frequently enough to obscure the links between this variation and the biological properties of distortion and lethality of the haplotypes that originally colonized Jutland.     

MILDLY DELETERIOUS MUTATIONS IN AVIAN MITOCHONDRIAL DNA: EVIDENCE FROM NEUTRALITY TESTS

To determine whether mildly deleterious mutations (MDMs) are present in nonrecombining genomes such as avian mitochondrial DNA (mtDNA), I analyzed molecular data from 14 studies using the neutrality tests of Tajima (1989a) and McDonald and Kreitman (1991). The presence of MDMs in mtDNA is inferred from trends observed across species in estimates of heterozygosity (θ and π) and by comparisons of polymorphism and divergence using the neutrality index (NI). Assuming neutrality, θ equals π and NI equals one. In this study, however, θ is greater than π more often than expected by chance, which reflects an excess of low-frequency alleles, and NI values presented here and elsewhere are consistently greater than one, which suggests an excess of nonsynonymous mutations within species (polymorphism) relative to between species (divergence). These observations suggest that, within species, there is an excess of rare haplotypes and that these haplotypes are carrying MDMs. The excess rare haplotypes may need to be accounted for when estimating population genetic parameters that assume strict neutrality.     

Phylogeography of the genus Dasiphora (Rosaceae) in the Qinghai‐Tibetan Plateau: divergence blurred by expansion

Plateau uprisings and climatic oscillations are considered to have caused extensive allopatric divergences that account for the rich species diversity of the Qinghai‐Tibetan Plateau (QTP). However, secondary contact during range shifts in the Quaternary glacial cycles or inter‐uplift stages may have restored the gene flow between species and so counteracted these divergences, particularly in rapidly‐adapting dominant elements. We tested this hypothesis by determining the phylogeographical history of Dasiphora (Rosaceae), a genus of two species that are widely distributed on the QTP and co‐exist in numerous localities. We sequenced two chloroplast DNA fragments (rbcL, trnT‐L) for 559 individuals from 87 populations. Bayesian methods were used to identify phylogenetic relationships and to estimate divergence times. Demographic histories were inferred using neutrality tests, mismatch distribution analysis, and coalescent simulation. A total of 112 haplotypes that clustered into three major groups were identified. The formation of these groups and their subgroups was dated to between the Pliocene and the late Pleistocene. In addition, we found that some groups underwent multiple extensive expansions. Species‐specific haplotypes were identified for each species, although these haplotypes phylogenetically intermixed. These results suggest that recent plateau uplifts and climatic oscillations might have caused the deep divergences observed within this genus. However, later range expansions probably blurred these divergences and possible species boundaries. Our results shed new light on the complex evolutionary history of the QTP alpine plants. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 111 , 777–788.     

A General Approach for Haplotype Phasing across the Full Spectrum of Relatedness

Many existing cohorts contain a range of relatedness between genotyped individuals, either by design or by chance. Haplotype estimation in such cohorts is a central step in many downstream analyses. Using genotypes from six cohorts from isolated populations and two cohorts from non-isolated populations, we have investigated the performance of different phasing methods designed for nominally ‘unrelated’ individuals. We find that SHAPEIT2 produces much lower switch error rates in all cohorts compared to other methods, including those designed specifically for isolated populations. In particular, when large amounts of IBD sharing is present, SHAPEIT2 infers close to perfect haplotypes. Based on these results we have developed a general strategy for phasing cohorts with any level of implicit or explicit relatedness between individuals. First SHAPEIT2 is run ignoring all explicit family information. We then apply a novel HMM method (duoHMM) to combine the SHAPEIT2 haplotypes with any family information to infer the inheritance pattern of each meiosis at all sites across each chromosome. This allows the correction of switch errors, detection of recombination events and genotyping errors. We show that the method detects numbers of recombination events that align very well with expectations based on genetic maps, and that it infers far fewer spurious recombination events than Merlin. The method can also detect genotyping errors and infer recombination events in otherwise uninformative families, such as trios and duos. The detected recombination events can be used in association scans for recombination phenotypes. The method provides a simple and unified approach to haplotype estimation, that will be of interest to researchers in the fields of human, animal and plant genetics.     

Comparisons of site- and haplotype-frequency methods for detecting positive selection

In this report, we compare the differences between various site- and haplotype-frequency tests in their power to detect positive selection by doing computer simulations. Our results are the following. 1) Although haplotype-frequency tests that are conditional on the number of haplotypes (K) were developed for nonrecombining haplotypes, these tests are insensitive to recombination. Such tests, including the Ewens-Watterson (EW) test, can therefore be applied to recombining haplotypes. 2) Tests conditional on the number of segregating sites (S) become overly conservative in the presence of recombination. 3) The EW test is usually the most powerful test during the sweep phase, especially when the local recombination rate is high. 4) The "extended haplotype homozygosity" test relies heavily on the prior knowledge of the target of selection. With that knowledge, it is the most powerful test, whereas in the absence of this prior information, the test has little power. We also study the sensitivities of the haplotype-frequency tests to background selection and various demographic forces. We find that these tests are sensitive to some forces other than positive selection. To alleviate the problem of low specificity, compound tests, such as the DH test (Zeng et al. 2006), may be a solution. In the companion paper (Zeng K, Shi S, Wu C-I, in preparation), we use the EW test to devise 2 compound tests, which are more powerful in detecting positive selection than DH, but are also relatively insensitive to demography.     

Genetic structure and demographic history of Lymantria dispar (Linnaeus, 1758) (Lepidoptera: Erebidae) in its area of origin and adjacent areas

We analyzed the population genetic structure and demographic history of 20 Lymantria dispar populations from Far East Asia using microsatellite loci and mitochondrial genes. In the microsatellite analysis, the genetic distances based on pairwise F_ST values ranged from 0.0087 to 0.1171. A NeighborNet network based on pairwise F_ST genetic distances showed that the 20 regional populations were divided into five groups. Bayesian clustering analysis (K = 3) demonstrated the same groupings. The populations in the Korean Peninsula and adjacent regions, in particular, showed a mixed genetic pattern. In the mitochondrial genetic analysis based on 98 haplotypes, the median‐joining network exhibited a star shape that was focused on three high‐frequency haplotypes (Haplotype 1: central Korea and adjacent regions, Group 1; Haplotype 37: southern Korea, Group 2; and Haplotype 90: Hokkaido area, Group 3) connected by low‐frequency haplotypes. The mismatch distribution dividing the three groups was unimodal. In the neutral test, Tajima's D and Fu's FS tests were negative. We can thus infer that the Far East Asian populations of L. dispar underwent a sudden population expansion. Based on the age expansion parameter, the expansion time was inferred to be approximately 53,652 years before present (ybp) for Group 1, approximately 65,043 ybp for Group 2, and approximately 76,086 ybp for Group 3. We propose that the mixed genetic pattern of the inland populations of Far East Asia is due to these expansions and that the inland populations of the region should be treated as valid subspecies that are distinguishable from other subspecies by genetic traits.     

Multivariate QST-FST comparisons: a neutrality test for the evolution of the g matrix in structured populations

Neutrality tests in quantitative genetics provide a statistical framework for the detection of selection on polygenic traits in wild populations. However, the existing method based on comparisons of divergence at neutral markers and quantitative traits (Q(st)-F(st)) suffers from several limitations that hinder a clear interpretation of the results with typical empirical designs. In this article, we propose a multivariate extension of this neutrality test based on empirical estimates of the among-populations (D) and within-populations (G) covariance matrices by MANOVA. A simple pattern is expected under neutrality: D = 2F(st)/(1 - F(st))G, so that neutrality implies both proportionality of the two matrices and a specific value of the proportionality coefficient. This pattern is tested using Flury's framework for matrix comparison [common principal-component (CPC) analysis], a well-known tool in G matrix evolution studies. We show the importance of using a Bartlett adjustment of the test for the small sample sizes typically found in empirical studies. We propose a dual test: (i) that the proportionality coefficient is not different from its neutral expectation [2F(st)/(1 - F(st))] and (ii) that the MANOVA estimates of mean square matrices between and among populations are proportional. These two tests combined provide a more stringent test for neutrality than the classic Q(st)-F(st) comparison and avoid several statistical problems. Extensive simulations of realistic empirical designs suggest that these tests correctly detect the expected pattern under neutrality and have enough power to efficiently detect mild to strong selection (homogeneous, heterogeneous, or mixed) when it is occurring on a set of traits. This method also provides a rigorous and quantitative framework for disentangling the effects of different selection regimes and of drift on the evolution of the G matrix. We discuss practical requirements for the proper application of our test in empirical studies and potential extensions.     

Directional selection and the site-frequency spectrum.

In this article we explore statistical properties of the maximum-likelihood estimates (MLEs) of the selection and mutation parameters in a Poisson random field population genetics model of directional selection at DNA sites. We derive the asymptotic variances and covariance of the MLEs and explore the power of the likelihood ratio tests (LRT) of neutrality for varying levels of mutation and selection as well as the robustness of the LRT to deviations from the assumption of free recombination among sites. We also discuss the coverage of confidence intervals on the basis of two standard-likelihood methods. We find that the LRT has high power to detect deviations from neutrality and that the maximum-likelihood estimation performs very well when the ancestral states of all mutations in the sample are known. When the ancestral states are not known, the test has high power to detect deviations from neutrality for negative selection but not for positive selection. We also find that the LRT is not robust to deviations from the assumption of independence among sites.     

Perspective: detecting adaptive molecular polymorphism: lessons from the MHC

Abstract. In the 1960s, when population geneticists first began to collect data on the amount of genetic variation in natural populations, balancing selection was invoked as a possible explanation for how such high levels of molecular variation are maintained. However, the predictions of the neutral theory of molecular evolution have since become the standard by which cases of balancing selection may be inferred. Here we review the evidence for balancing selection acting on the major histocompatibility complex (MHC) of vertebrates, a genetic system that defies many of the predictions of neutrality. We apply many widely used tests of neutrality to MHC data as a benchmark for assessing the power of these tests. These tests can be categorized as detecting selection in the current generation, over the history of populations, or over the histories of species. We find that selection is not detectable in MHC datasets in every generation, population, or every evolutionary lineage. This suggests either that selection on the MHC is heterogeneous or that many of the current neutrality tests lack sufficient power to detect the selection consistently. Additionally, we identify a potential inference problem associated with several tests of neutrality. We demonstrate that the signals of selection may be generated in a relatively short period of microevolutionary time, yet these signals may take exceptionally long periods of time to be erased in the absence of selection. This is especially true for the neutrality test based on the ratio of nonsynonymous to synonymous substitutions. Inference of the nature of the selection events that create such signals should be approached with caution. However, a combination of tests on different time scales may overcome such problems.     

Evolutionary correlation between control region sequence and restriction polymorphisms in the mitochondrial genome of a large Senegalese Mandenka sample

We present here the first comparative analysis at the population level between Restriction Fragment Length Polymorphism (RFLP) and control region sequence polymorphism in a large and homogeneous Senegalese Mandenka sample. Eleven RFLP haplotypes and 60 different sequences are found in 119 individuals, revealing that a very high level of mtDNA diversity can be maintained in a small population. A sequence neighbor- joining tree and an analysis of molecular variance show that sequences associated with a given restriction haplotype are evolutionarily highly correlated: sequencing generally leads to the subtyping of RFLP haplotypes. Evolutionary relationships among RFLP haplotypes inferred from restriction site differences are in good agreement with those inferred from sequence data. A single difference is observed and is likely due to a single restriction homoplasy having occurred in the control region. Selective neutrality tests on both RFLP and sequence data accept the hypotheses of mtDNA neutrality and population equilibrium. The deep coalescence times (exceeding 50,000 yr) of sequences associated with the two most frequent restriction haplotypes confirm that the Niokolo Mandenka population has not passed through a recent bottleneck and that gene flow is maintained among West African populations despite ethnic differences.      

秦岭地区赤胸梳爪步甲种群的遗传多样性和扩张历史分析（英文）

赤胸梳爪步Dolichus halensis (Schaller)（鞘翅目：步甲科）是重要的捕食性天敌昆虫,在我国分布广泛。为揭示其种群遗传多样性和扩张机制, 本研究以秦岭地区为中心,以线粒体Cox1 tRNALeu Cox2基因片段为分子标记,对来自于24个采集点共191个个体进行了检测分析。在长度为1 601 bp的碱基中共检测到45个变异位点,定义了53个单倍型,单倍型多样性高（H_d=0.796）,而核苷酸多样性较低（P_i=0.0033）。系统发育分析结果表明该地区该物种存在两大进化枝。分子变异分析（AMOVA）表明86.61%的变异来源于种群内。SAMOVA和PERMUT分析结果一致,表明秦岭地区分布的赤胸梳爪步甲种群不存在明显的谱系地理结构。中性检验和错配分布分析的结果一致,表明该物种在秦岭地区曾经发生过种群扩张。综上,认为赤胸疏爪步甲种群经历过冰期后的扩散。     

Spinocerebellar ataxia type 8: molecular genetic comparisons and haplotype analysis of 37 families with ataxia

We reported elsewhere that an untranslated CTG expansion causes the dominantly inherited neurodegenerative disorder spinocerebellar ataxia type 8 (SCA8). SCA8 shows a complex inheritance pattern with extremes of incomplete penetrance, in which often only one or two affected individuals are found in a given family. SCA8 expansions have also been found in control chromosomes, indicating that separate genetic or environmental factors increase disease penetrance among SCA8-expansion-carrying patients with ataxia. We describe the molecular genetic features and disease penetrance of 37 different families with SCA8 ataxia from the United States, Canada, Japan, and Mexico. Haplotype analysis using 17 STR markers spanning an approximately 1-Mb region was performed on the families with ataxia, on a group of expansion carriers in the general population, and on psychiatric patients, to clarify the genetic basis of the reduced penetrance and to investigate whether CTG expansions among different populations share a common ancestral background. Two major ancestrally related haplotypes (A and A') were found among white families with ataxia, normal controls, and patients with major psychosis, indicating a common ancestral origin of both pathogenic and nonpathogenic SCA8 expansions among whites. Two additional and distinct haplotypes were found among a group of Japanese families with ataxia (haplotype B) and a Mexican family with ataxia (haplotype C). Our finding that SCA8 expansions on three independently arising haplotypes are found among patients with ataxia and cosegregate with ataxia when multiple family members are affected further supports the direct role of the CTG expansion in disease pathogenesis.     

Testing for recombination in a short nuclear DNA sequence of the European meadow grasshopper, Chorthippus parallelus

Single-copy nuclear DNA sequences have high potential as a source of genetic markers for population analyses. However, the difficulties that arise when haplotypes that are the product of recombinational rearrangements are present require additional consideration. Two statistical methods for identifying potential recombinants by detecting anomalies in the distribution of variable sites along sequences were used to screen sequences from a single-copy nuclear DNA fragment, cpnl-1, of the European meadow grasshopper (Chorthippus parallelus). Five of the 71 haplotypes in the cpnl-1 data set showed nonrandom distribution of polymorphic sites using both methods. The second method pinpointed an additional four haplotypes. Estimates of the rate of recombination in the entire data set were obtained using standard methods. It is concluded that cpnl-1 haplotypes have been involved in recombination or gene conversion events at a rate more than twice the mutation rate. This confirms that recombination and gene conversion are significant factors in the generation of haplotype variation in nuclear gene sequences. The cpnl-1 haplotypes identified by the tests were present only in populations that have had recent contact; the Balkan and Turkish refugial populations and their post-glacial colonies to the north. This is discussed in relation to the phylogenetic inferences drawn from the same data in a previous report.     

Recombination sites in the HLA class II region are haplotype dependent

We have analyzed DNA sequence polymorphisms of DQ alpha and DQ beta chains from three haplotypes from the DRw52 family: DR5 DQw1 (FPA, GM3106), DRw6 DQw1 (CB6B, 10w9060), and DRw6 DQw3 (AMALA, 10w9064). The results indicate that the DR5 DQw1 and DRw6 DQw1 haplotypes have arisen by recombination between the DR beta 1 and DQ alpha loci. This contrasts with our previous analysis of DR4 DQ"Wa", DR3 DQ"Wa", and DR7 DQw3 haplotypes, all of which appear to have arisen by virtue of recombination between DQ alpha and DQ beta. Thus, there appear to be at least two different sites where recombination has occurred within the DR and DQ subregions. These differing patterns of recombination were interpreted in the context of the three major family groups of class II haplotypes, the DRw53, DRw52, and DR1/2 haplotype families. The data indicate that haplotypes from these family groups tend to undergo recombination at different locations. We propose that these differences in site of recombination are a reflection of differences in the molecular organization of the haplotypes belonging to each family group.     

Recombination and the frequency spectrum in Drosophila melanogaster and Drosophila simulans

Most "tests of neutrality" assess whether particular data sets depart from the predictions of a standard neutral model with no recombination. For Drosophila, where nuclear polymorphism data routinely show evidence of genetic exchange, the assumption of no recombination is often unrealistic. In addition, while conservative, this assumption is made at the cost of a great loss in power. Perhaps as a result, tests of the frequency spectrum based on zero recombination suggest an adequate fit of Drosophila polymorphism data to the predictions of the standard neutral model. Here, we analyze the frequency spectrum of a large number of loci in Drosophila melanogaster and D. simulans using two summary statistics. We use an estimate of the population recombination rate based on a laboratory estimate of the rate of crossing over per physical length and an estimate of the species' effective population size. In contrast to previous studies, we find that roughly half of the loci depart from the predictions of the standard neutral model. The extent of the departure depends on the exact recombination rate, but the global pattern that emerges is robust. Interestingly, these departures from neutral expectations are not unidirectional. The large variance in outcomes may be due to a complex demographic history and inconsistent sampling, or to the pervasive action of natural selection.     

An approach to mapping haplotype-specific recombination sites in human MHC class III

Studies of the major histocompatibility complex (MHC) in mouse indicate that the recombination sites are not randomly distributed and their occurrence is haplotype-dependent. No data concerning haplotype-specific recombination sites in human are available due to the low number of informative families. To investigate haplotype-specific recombination sites in human MHC, we here describe an approach based on identification of recombinant haplotypes derived from one conserved haplotype at the population level. The recombination sites were mapped by comparing polymorphic markers between the recombinant and assumed original haplotypes. We tested this approach on the extended haplotype HLA A3; B47; Bf ^* F; C4A ^* 1; C4B ^* Q0; DR7, which is most suitable for this analysis. First, it carries a number of rare markers, and second, the haplotype, albeit rare in the general population, is frequent in patients with 21-hydroxylase (21OH) defect. We observed recombinants derived from this haplotype in patients with 21OH defect. All these haplotypes had the centromeric part (from Bf to DR) identical to the original haplotype, but they differed in HLA A and B. We therefore assumed that they underwent recombinations in the segment that separates the Bf and HLA B genes. Polymorphic markers indicated that all break points mapped to two segments near the TNF locus. This approach makes possible the mapping of preferential recombination sites in different haplotypes.