Apportionment of global human genetic diversity based on craniometrics and skin color

A number of analyses of classical genetic markers and DNA polymorphisms have shown that the majority of human genetic diversity exists within local populations (approximately 85%), with much less among local populations (approximately 5%) or between major geographic regions or "races" (approximately 10%). Previous analysis of craniometric variation (Relethford [1994] Am J Phys Anthropol 95:53-62) found that between 11-14% of global diversity exists among geographic regions, with the remaining diversity existing within regions. The methods used in this earlier paper are extended to a hierarchical partitioning of genetic diversity in quantitative traits, allowing for assessment of diversity among regions, among local populations within regions, and within local populations. These methods are applied to global data on craniometric variation (57 traits) and skin color. Multivariate analysis of craniometric variation shows results similar to those obtained from genetic markers and DNA polymorphisms: roughly 13% of the total diversity is among regions, 6% among local populations within regions, and 81% within local populations. This distribution is concordant with neutral genetic markers. Skin color shows the opposite pattern, with 88% of total variation among regions, 3% among local populations within regions, and 9% within local populations, a pattern shaped by natural selection. The apportionment of genetic diversity in skin color is atypical, and cannot be used for purposes of classification. If racial groups are based on skin color, it appears unlikely that other genetic and quantitative traits will show the same patterns of variation.     

Comparison of Parameters from Rhythmometric Models with Multiple Components on Hybrid Data

Population multiple components is a statistical tool useful for the analysis of time-dependent hybrid data. With a small number of parameters, it is possible to model and to predict the periodic behavior of a population. In this article, we propose two methods to compare among populations rhythmometric parameters obtained by multiple component analysis. The first is a parametric method based in the usual statistical techniques for comparison of mean vectors in multivariate normal populations. The method, through MANOVA analysis, allows comparison of the MESOR and amplitude-acrophase pair of each component among two or more populations. The second is a nonparametric method, based in bootstrap techniques, to compare parameters from two populations. This test allows one to compare the MESOR, the amplitude, and the acrophase of each fitted component, as well as the global amplitude, orthophase, and bathyphase estimated when all fitted components are harmonics of a fundamental period. The idea is to calculate a confidence interval for the difference of the parameters of interest. If this interval does not contain zero, it can be concluded that the parameters from the two models are different with high probability. An estimation of p-value for the corresponding test can also be calculated. Both methods are illustrated with an example, based on clinical data. The nonparametric test can also be applied to paired data, a special situation of great interest in practice. By the use of similar bootstrap techniques, we illustrate how to construct confidence intervals for any rhythmometric parameter estimated from population multiple components models, including the orthophase, bathyphase, and global amplitude. These tests for comparison of parameters among populations are a needed tool when modeling the nonsinusoidal rhythmic behavior of hybrid data by population multiple component analysis.     

Comparison of parameters from rhythmometric models with multiple components on hybrid data

Population multiple components is a statistical tool useful for the analysis of time-dependent hybrid data. With a small number of parameters, it is possible to model and to predict the periodic behavior of a population. In this article, we propose two methods to compare among populations rhythmometric parameters obtained by multiple component analysis. The first is a parametric method based in the usual statistical techniques for comparison of mean vectors in multivariate normal populations. The method, through MANOVA analysis, allows comparison of the MESOR and amplitude-acrophase pair of each component among two or more populations. The second is a nonparametric method, based in bootstrap techniques, to compare parameters from two populations. This test allows one to compare the MESOR, the amplitude, and the acrophase of each fitted component, as well as the global amplitude, orthophase, and bathyphase estimated when all fitted components are harmonics of a fundamental period. The idea is to calculate a confidence interval for the difference of the parameters of interest. If this interval does not contain zero, it can be concluded that the parameters from the two models are different with high probability. An estimation of p-value for the corresponding test can also be calculated. Both methods are illustrated with an example, based on clinical data. The nonparametric test can also be applied to paired data, a special situation of great interest in practice. By the use of similar bootstrap techniques, we illustrate how to construct confidence intervals for any rhythmometric parameter estimated from population multiple components models, including the orthophase, bathyphase, and global amplitude. These tests for comparison of parameters among populations are a needed tool when modeling the nonsinusoidal rhythmic behavior of hybrid data by population multiple component analysis.     

Genomic consequences of a recent three‐way admixture in supplemented wild brown trout populations revealed by local ancestry tracts

Understanding the evolutionary consequences of human‐mediated introductions of domesticated strains into the wild and their subsequent admixture with natural populations is of major concern in conservation biology. However, the genomic impacts of stocking from distinct sources (locally derived vs. divergent) on the genetic integrity of wild populations remain poorly understood. We designed an approach based on estimating local ancestry along individual chromosomes to provide a detailed picture of genomic admixture in supplemented populations. We used this approach to document admixture consequences in the brown trout Salmo trutta, for which decades of stocking practices have profoundly impacted the genetic make‐up of wild populations. In southern France, small local Mediterranean populations have been subject to successive introductions of domestic strains derived from the Atlantic and Mediterranean lineages. To address the impact of stocking, we evaluate the extent of admixture from both domestic strains within populations, using 75,684 mapped SNPs obtained from double‐digested restriction site‐associated DNA sequencing. Then, the chromosomal ancestry profiles of admixed individuals reveal a wider diversity of hybrid and introgressed genotypes than estimated using classical methods for inferring ancestry and hybrid pedigrees. In addition, the length distribution of introgressed tracts retained different timings of introgression between the two domestic strains. We finally reveal opposite consequences of admixture on the level of polymorphism of the recipient populations between domestic strains. Our study illustrates the potential of using the information contained in the genomic mosaic of ancestry tracts in combination with classical methods based on allele frequencies for analysing multiple‐way admixture with population genomic data.     

NATURAL SELECTION ON QUANTITATIVE TRAITS IN THE BOMBINA HYBRID ZONE

Observations on the means, variances, and covariances of quantitative traits across hybrid zones can give information similar to that from Mendelian markers. In addition, they can identify particular traits through which the cline is maintained. We describe a survey of six traits across the hybrid zone between Bombina bombina and Bombina variegata (Amphibia: Discoglossidae) near Pe??enica in Croatia. We obtained laboratory measuments of the belly pattern, skin thickness, mating call, skeletal form, egg size, and the developmental time of tadpoles. Although offspring from hybrid populations showed no evidence of reduced viability, a third of the F₁ families failed completely, irrespective of the direction of the cross. All traits differed significantly between the taxa. Clines in belly pattern, skin thickness, mating call, and skeletal form were closely concordant with clines in four diagnostic enzyme loci. However, the cline in developmental time was displaced towards bombina, and the cline in egg size was displaced towards variegata. This discordance could be because the traits are not inherited additively or because they are subject to different selection pressures. We favor the latter explanation in the case of developmental time. We show that moderate selection acting directly on a trait suffices to shift its position; rather stronger selection is needed to change its width appreciably. Within hybrid populations, there are significant associations among quantitative traits, and between traits and enzymes. Phenotypic variances also increase in hybrid populations. These observations can be explained by linkage disequilibria among the underlying loci. However, the average magnitude of the covariance between traits is about half that expected from the linkage disequilibria between enzyme loci. The discrepancy is not readily explained by nonadditive gene action. This puzzle is now unresolved and calls for further investigation.     

Efficiency of model-based Bayesian methods for detecting hybrid individuals under different hybridization scenarios and with different numbers of loci

Accurate detection of offspring resulting from hybridization between individuals of distinct populations has a range of applications in conservation and population genetics. We assessed the hybrid identification efficiency of two methods (implemented in the STRUCTURE and NEWHYBRIDS programs) which are tailored to identifying hybrid individuals but use different approaches. Simulated first- and second-generation hybrids were used to assess the performance of these two methods in detecting recent hybridization under scenarios with different levels of genetic divergence and varying numbers of loci. Despite the different approaches of the methods, the hybrid detection efficiency was generally similar and neither of the two methods outperformed the other in all scenarios assessed. Interestingly, hybrid detection efficiency was only minimally affected by whether reference population allele frequency information was included or not. In terms of genotyping effort, efficient detection of F1 hybrid individuals requires the use of 12 or 24 loci with pairwise F(ST) between hybridizing parental populations of 0.21 or 0.12, respectively. While achievable, these locus numbers are nevertheless higher than the number of loci currently commonly applied in population genetic studies. The method of STRUCTURE seemed to be less sensitive to the proportion of hybrids included in the sample, while NEWHYBRIDS seemed to perform slightly better when individuals from both backcross and F1 hybrid classes were present in the sample. However, separating backcrosses from purebred parental individuals requires a considerable genotyping effort (at least 48 loci), even when divergence between parental populations is high.     

汉族和藏族的皮肤反射系数分析

人类的肤色是在遗传和环境因素的共同作用下形成的。最近在人类复杂性状基因定位中,统计学和遗传学研究方法的发展,使得肤色相关基因有可能利用这些方法来甄别。而且,有了Photovolt ColorWalk色度计等便携的新式光谱反射系数测量工具,研究者可以方便准确地测量大量人群的皮肤反射系数作为遗传学研究用的人类表型性状。我们衣物Photovolt ColorWark测试了372人的一个汉族群体和274人的一个藏族群体的上臂内侧不受阳光照射的皮肤反射系数,以建立反射系数测量的基准数据。我们调查了年龄、性别、居住地纬度、种族等不同因素对皮肤反射系数的影响,也调查了皮肤反射系数的正态分布。在这些研究结果的基础上我们设计了皮肤反射系数遗传学研究的策略。     

Cytoplasmic microtubules in transformed mouse x nontransformed human cell hybrids: correlation with in vitro growth.

Patterns of cytoplasmic microtubules in somatic cell hybrids between transformed mouse cells and nontransformed human skin fibroblasts were examined using antitubulin antibodies as an immunofluorescent probe. Nontransformed cells have been shown to exhibit an extensive cytoplasmic microtubule complex (CMTC), while in transformed cells, this complex is greatly diminished. The hybrid populations contained both types of cells. In addition, they contained cells with previously undescribed intermediate CMTC phenotypes. The percentage of each phenotype present in hybrid populations was determined for sixteen hybrid clones. Seven clones were found which appeared transformed on the basis of their CMTC pattern. The others were comprised of various proportions of all the cell types described. Repeated quantitation of the proportions of these types in the hybrid populations showed them to be stable with time in culture. Growth in vitro of the hybrid clones was assayed by determining their saturation densities, their plating efficiencies on plastic and their colony-forming abilities in soft agar. In vitro growth of a cell population was found to be directly proportional to the percentage of cells in the population which showed the greatly diminished CMTC pattern which has been described for transformed cells. This is strong evidence that a greatly reduced CMTC is associated with transformed behavior, especially the increased capacity of transformed cells for in vitro growth.     

Demographic history shapes genomic ancestry in hybrid zones

Demographic factors such as migration rate and population size can impede or facilitate speciation. In hybrid zones, reproductive boundaries between species are tested and demography mediates the opportunity for admixture between lineages that are partially isolated. Genomic ancestry is a powerful tool for revealing the history of admixed populations, but models and methods based on local ancestry are rarely applied to structured hybrid zones. To understand the effects of demography on ancestry in hybrids zones, we performed individual‐based simulations under a stepping‐stone model, treating migration rate, deme size, and hybrid zone age as parameters. We find that the number of ancestry junctions (the transition points between genomic regions with different ancestries) and heterogenicity (the genomic proportion heterozygous for ancestry) are often closely connected to demographic history. Reducing deme size reduces junction number and heterogenicity. Elevating migration rate increases heterogenicity, but migration affects junction number in more complex ways. We highlight the junction frequency spectrum as a novel and informative summary of ancestry that responds to demographic history. A substantial proportion of junctions are expected to fix when migration is limited or deme size is small, changing the shape of the spectrum. Our findings suggest that genomic patterns of ancestry could be used to infer demographic history in hybrid zones.     

Towards the introducing of resistance to powdery mildew from Lycopersicon hirsutum into L. esculentum

Genes of resistance to Oidium lycopersicum from Lycopersicon hirsutum LA 1775 were introduced to L. esculentum. Breeding procedures were based on a one-way programme up to the F2 generation and then four different methods were adopted to obtain F4 and BC4 populations. Screening tests among those hybrid populations were performed in a greenhouse and showed segregation for resistance to powdery mildew due to different genetic backgrounds of the families derived from four breeding methods that changed the status of the gene/genes responsible for resistance to powdery mildew. F4 and BC4 populations varied in relation to morphological traits (fruit size and weight, seed and fruit productivity, number of locules). There was a significant progress in breeding in comparison to L. hirsutum regarding fruit size and weight, and the number of locules. Values of two other traits: seed and fruit productivity, that are correlated with self- and cross-compatibility, were low and similar to L. hirsutum. Therefore, another one or two backcrosses will probably improve seed and fruit productivity.     

Mating patterns in a hybrid zone of fire-bellied toads (Bombina): inferences from adult and full-sib genotypes

We present two novel methods to infer mating patterns from genetic data. They differ from existing statistical methods of parentage inference in that they apply to populations that deviate from Hardy-Weinberg and linkage equilibrium, and so are suited for the study of assortative mating in hybrid zones. The core data set consists of genotypes at several loci for a number of full-sib clutches of unknown parentage. Our inference is based throughout on estimates of allelic associations within and across loci, such as heterozygote deficit and pairwise linkage disequilibrium. In the first method, the most likely parents of a given clutch are determined from the genotypic distribution of the associated adult population, given an explicit model of nonrandom mating. This leads to estimates of the strength of assortment. The second approach is based solely on the offspring genotypes and relies on the fact that a linear relation exists between associations among the offspring and those in the population of breeding pairs. We apply both methods to a sample from the hybrid zone between the fire-bellied toads Bombina bombina and B. variegata (Anura: Disco glossidae) in Croatia. Consistently, both approaches provide no evidence for a departure from random mating, despite adequate statistical power. Instead, B. variegata-like individuals among the adults contributed disproportionately to the offspring cohort, consistent with their preference for the type of breeding habitat in which this study was conducted.     

Skin pigmentation,biogeographical ancestry and admixture mapping

Ancestry informative markers (AIMs) are genetic loci showing alleles with large frequency differences between populations. AIMs can be used to estimate biogeographical ancestry at the level of the population, subgroup (e.g. cases and controls) and individual. Ancestry estimates at both the subgroup and individual level can be directly instructive regarding the genetics of the phenotypes that differ qualitatively or in frequency between populations. These estimates can provide a compelling foundation for the use of admixture mapping (AM) methods to identify the genes underlying these traits. We present details of a panel of 34 AIMs and demonstrate how such studies can proceed, by using skin pigmentation as a model phenotype. We have genotyped these markers in two population samples with primarily African ancestry, viz. African Americans from Washington D.C. and an African Caribbean sample from Britain, and in a sample of European Americans from Pennsylvania. In the two African population samples, we observed significant correlations between estimates of individual ancestry and skin pigmentation as measured by reflectometry (R(2)=0.21, P<0.0001 for the African-American sample and R(2)=0.16, P<0.0001 for the British African-Caribbean sample). These correlations confirm the validity of the ancestry estimates and also indicate the high level of population structure related to admixture, a level that characterizes these populations and that is detectable by using other tests to identify genetic structure. We have also applied two methods of admixture mapping to test for the effects of three candidate genes (TYR, OCA2, MC1R) on pigmentation. We show that TYR and OCA2 have measurable effects on skin pigmentation differences between the west African and west European parental populations. This work indicates that it is possible to estimate the individual ancestry of a person based on DNA analysis with a reasonable number of well-defined genetic markers. The implications and applications of ancestry estimates in biomedical research are discussed.     

Quantitative evaluation of hybridization and the impact on biodiversity conservation

Anthropogenic hybridization is an increasing conservation threat worldwide. In South Africa, recent hybridization is threatening numerous ungulate taxa. For example, the genetic integrity of the near‐threatened bontebok (Damaliscus pygargus pygargus) is threatened by hybridization with the more common blesbok (D. p. phillipsi). Identifying nonadmixed parental and admixed individuals is challenging based on the morphological traits alone; however, molecular analyses may allow for accurate detection. Once hybrids are identified, population simulation software may assist in determining the optimal conservation management strategy, although quantitative evaluation of hybrid management is rarely performed. In this study, our objectives were to describe species‐wide and localized rates of hybridization in nearly 3,000 individuals based on 12 microsatellite loci, quantify the accuracy of hybrid assignment software (STRUCTURE and NEWHYBRIDS), and determine an optimal threshold of bontebok ancestry for management purposes. According to multiple methods, we identified 2,051 bontebok, 657 hybrids, and 29 blesbok. More than two‐thirds of locations contained at least some hybrid individuals, with populations varying in the degree of introgression. HYBRIDLAB was used to simulate four generations of coexistence between bontebok and blesbok, and to optimize a threshold of ancestry, where most hybrids will be detected and removed, and the fewest nonadmixed bontebok individuals misclassified as hybrids. Overall, a threshold Q‐value (admixture coefficient) of 0.90 would remove 94% of hybrid animals, while a threshold of 0.95 would remove 98% of hybrid animals but also 8% of nonadmixed bontebok. To this end, a threshold of 0.90 was identified as optimal and has since been implemented in formal policy by a provincial nature conservation agency. Due to widespread hybridization, effective conservation plans should be established and enforced to conserve native populations that are genetically unique.     

Natural hybridization between black spruce and red spruce

Using species-specific random amplified polymorphic DNA (RAPD) markers and morphological characters, natural hybridization between the closely related black spruce Picea mariana (Mill.) B.S.P. and red spruce P. rubens Sarg. was evaluated in natural populations from north-eastern North America. Sampling included populations from both areas of allopatry and also 14 populations from part of the area of sympatry located in the province of Québec and covering several thousands of square kilometres. Classification results from RAPD species-specific markers and from a discriminant function based on morphology were compared. Molecular analysis of the allopatric populations indicated a small amount of interspecific gene leakage with no asymmetric directionality to introgression. A high occurrence of hybrid/introgressant individuals was observed within sympatric populations, suggesting weak reproductive isolation. As expected, the detection of such individuals was more efficient using molecular markers than with morphological traits. The hybrid zone appeared extensive with variable species structure and, in some stands, the main component composed of hybrid/introgressant trees. Implications for the genecology and genetic management of these species are discussed.     

Some problems of the study and representation of Siberian flora biodiversity in connection with its conservation

The present knowledge of flora of Siberia is quite insufficient, which is a considerable obstacle to the detection of its biodiversity for conserving the gene pool. Planning conservation measures should be focused not only on species, but also on intraspecific taxa and their phylogenetic uniqueness. In the absence of genetic data, available morphological and geographical methods should be used, so that, when it is impossible to establish the actual origin and relationship of taxa, the existing morphological diverstity is at least represented, as far as it is known to be based on genetic diversity. Phenetic maps can be of much use in the study of intraspecific morphological diversity. The estimation of geographic variability and morphological diversity, as well as the evaluation of territories, can be based on such maps. To represent adequately the biodiversity existing within poorly studied, presumably hybrid plant groups, until actual origin and relationship are known, they should probably be forcedly and provisionally considered as a special type of hybrid complexes, representing the unclear present day taxonomic situation. Such complexes could include populations and individuals with morphological characters of two or more different species, until the systematic position of such populations and individuals is further explored. Until the actual taxonomic status and relationship of the components are established, they could be regarded as subspecies permitted by the Code, or as certain recorded morphological deviations from the type, without assigning any taxonomic status to them--depending on the available data on variability and distribution. In the future, the resulting provisional information on morphological diversity would help to concentrate the efforts of biologists, in possession of the newest methods, on the most important objects, and serve as the scientific base for effective measures aimed at the conservation and management of the vast gene pool of the Siberian flora.     

PATTERNS OF MATING IN WILD SUNFLOWER HYBRID ZONES

Theory predicts that homoploid hybrid speciation will be facilitated by selfing, yet most well-documented hybrid species are outcrossers. One possible explanation for this puzzle is that conditions in hybrid populations may favor selfing, even in otherwise outcrossing species. For example, in self-incompatible plants, mixtures of self and interspecific pollen often induce selfing. Here, we examine patterns of mating in three hybrid zones and four “pure” populations of Helianthus annuus and H. petiolaris, wild, self-incompatible sunflower species that are thought to have parented three homoploid hybrid species. Fourteen to 16 maternal families from each pure population and 44–46 maternal families from each hybrid zone were analyzed for seven polymorphic isozyme loci. Maximum-likelihood (ML) methods were used to estimate multilocus outcrossing rates (T_m) and hybridization frequencies for each maternal family, each phenotypic group within each hybrid zone (annuus-like, hybrid, and petiolaris-like), and each population. As predicted for self-incompatible species, all four parental populations have outcrossing rate ML estimates of 1.0. Within the hybrid zones, outcrossing rates were lowest in the H. annuus–like fraction of the population (0.73, 0.72, and 0.74 in the three hybrid zones, respectively), largely intermediate in the H. petiolaris–like group (0.94, 0.90, and 0.94), and highest in the hybrid group (0.97, 0.93, and 0.97). Although outcrossing rates are lower in hybrid zones than in pure populations, it is unlikely that the observed decrease facilitates hybrid speciation because outcrossing rates in the critical hybrid fraction of the population do not differ significantly from 1.0. Dividing the outcrossed pollen pool into intraspecific and interspecific components revealed that maternal plants are largely fertilized by conspecific pollen, confirming an important role for pollen competition as a reproductive barrier. Highly sterile hybrid plants do not appear to discriminate between parental species pollen, but hybrids with higher fertility tend to be fertilized by pollen from the parental group they resemble genetically. Thus, gametic selection leads to substantial assortative mating in these hybrid zones.     

PATTERNS OF VARIATION AND LINKAGE DISEQUILIBRIUM IN A FIELD CRICKET HYBRID ZONE

The distribution of multilocus genotypes found within a natural hybrid zone is determined by the sample of genotypes present when the hybrid zone first formed, by subsequent patterns of genetic exchange between the hybridizing taxa, and by drift and selection within each of the hybrid zone populations. We have used anonymous nuclear DNA restriction fragment polymorphisms (RFLPs) to characterize the array of multilocus genotypes present within a well-studied hybrid zone between two eastern North American field crickets, Gryllus pennsylvanicus and Gryllus firmus. These crickets hybridize along a zone of contact that extends from New England to Virginia. Previous studies have shown that both premating and postmating barriers exist between the two cricket species, but the absence of diagnostic morphological and allozyme markers has made it difficult to assess the consequences of these barriers for genetic exchange. Analyses based on four diagnostic anonymous nuclear markers indicate that hybrid zone populations in Connecticut contain few F₁ hybrids, and that nonrandom associations persist among nuclear gene markers, between nuclear and cytoplasmic markers, and between molecular markers and morphology. Field cricket populations within the hybrid zone are not “hybrid swarms” but consist primarily of crickets that are very much like one or the other of the parental species. Despite ample opportunity for genetic exchange and evidence for introgression at some loci, the two species remain quite distinct. Such a pattern appears to be characteristic of many natural hybrid zones.     

Skin color variation in eastern Nepal

The purpose of this paper is to provide skin color data on several endogamous groups of eastern Nepal and to demonstrate genetic microdifferentiation in skin color. Skin reflectance measures were taken at the upper inner arm and forehead sites, using the British DSL Model 99 Reflectance Spectrophotometer fitted with blue, green, and red filters. Measurements on 484 males representing six endogamous groups (Jirel, Sunwar, Sherpa, Tamang, Brahman, Chetri) were utilized. After adjusting for group-specific age effects, multivariate likelihood ratio permutation tests were used to assess ethnic differences in skin reflectance means and covariance matrices. Ethnic group membership had a highly significant effect on skin color at the upper inner arm site and at the forehead site. Differential tanning responses among groups were also detected and may represent the influence of genotype-environment interaction on reflectance traits. Mahalanobis distance analysis revealed patterns of microdifferentiation that are inconsistent with oral historical accounts. Little support was found for the ethnohistorical belief that Jirels are a hybrid group descended from Sunwars and Sherpas, even though this scenario is supported by linguistic, anthropometric, and dermatoglyphic data. We suggest that while reflectometric studies of skin pigmentation may be useful in assessing the macrodifferentiation of human populations, the use of skin color in differentiation studies at the local population level needs to be carefully evaluated.     

Genetic structure of natural populations of Castanea sativa in Turkey: evidence of a hybrid zone

This study points out the evidence of a hybrid zone between two groups of genetically differentiated populations of chestnut (Castanea sativa Mill.) in Turkey. Genetic structure, gene flow and introgression levels, based on 16 allozyme loci, were investigated on 34 population samples spanning the entire C. sativa distribution area in this country. The occurrence of the hybrid zone, located in the Bithynian region, was inferred in a group of populations showing the following genetic characteristics: (i) enhanced genetic variability and intermediate allelic frequencies between those of the Western and Eastern groups of populations; (ii) sharp and concordant changes in allele frequencies; (iii) decreased gene flow with the Western and Eastern populations. Starting from the cline width estimated to be 324 km, strength of selection was evaluated from the gene flow distance, as indicated from the degree of genetic structuring outside the hybrid zone. Evolutionary processes shaping the observed genetic differentiation and introgression are discussed on the basis of palynological data, palaeoclimatic events and evidence of hybridization found in other plant and animal species in the same region.     

Estimating ancestry and heterozygosity of hybrids using molecular markers

ABSTRACT: BACKGROUND: Hybridization, genetic mixture of distinct populations, gives rise to myriad recombinant genotypes. Characterizing the genomic composition of hybrids is critical for studies of hybrid zone dynamics, inheritance of traits, and consequences of hybridization for evolution and conservation. Hybrid genomes are often summarized either by an estimate of the proportion of alleles coming from each ancestral population or classification into discrete categories like F1, F2, backcross, or merely ``hybrid' vs. ``pure'. In most cases, it is not realistic to classify individuals into the restricted set of classes produced in the first two generations of admixture. However, the continuous ancestry index misses an important dimension of the genotype. Joint consideration of ancestry together with interclass heterozygosity (proportion of loci with alleles from both ancestral populations) captures all of the information in the discrete classification without the unrealistic assumption that only two generations of admixture have transpired. Here I describe a maximum likelihood method for joint estimation of ancestry and interclass heterozygosity. RESULTS: I present two worked examples illustrating the value of the approach for describing variation among hybrid populations and evaluating the validity of the assumption underlying discrete classification. Naively classifying natural hybrids into the standard six line cross categories can be misleading, and false classification can be a serious problem for datasets with few molecular markers. CONCLUSIONS: My analysis underscores previous work showing that many (50 or more) ancestry informative markers are needed to avoid erroneous classification. However, valuable inferences can be obtained by focusing directly on distributions of ancestry and heterozygosity. Estimating and visualizing the joint distribution of ancestry and interclass heterozygosity is an effective way to compare the genetic structure of hybrid populations and these estimates can be used in classic quantitative genetic methods for assessing additive, dominant, and epistatic genetic effects on hybrid phenotypes and fitness. The methods are implemented in a freely available package ``HIest' for the R statistical software (http://cran.r-project.org/web/packages/HIest/index.html).