Does natural selection promote population divergence? A comparative analysis of population structure using amplified fragment length polymorphism markers and quantitative traits

Divergent natural selection is considered an important force in plant evolution leading to phenotypic differentiation between populations exploiting different environments. Extending an earlier greenhouse study of population differentiation in the selfing annual plant Senecio vulgaris, we estimated the degree of population divergence in several quantitative traits related to growth and life history and compared these estimates with those based on presumably neutral molecular markers (amplified fragment length polymorphisms; AFLPs). This approach allowed us to disentangle the effects of divergent selection from that of other evolutionary forces (e.g. genetic drift). Five populations were examined from each of two habitat types (ruderal and agricultural habitats). We found a high proportion of total genetic variance to be among populations, both for AFLP markers (phiST = 0.49) and for quantitative traits (range of QST: 0.26-0.77). There was a strong correlation between molecular and quantitative genetic differentiation between pairs of populations (Mantel's r = 0.59). However, estimates of population differentiation in several quantitative traits exceeded the neutral expectation (estimated from AFLP data), suggesting that divergent selection contributed to phenotypic differentiation, especially between populations from ruderal and agricultural habitats. Estimates of within-population variation in AFLP markers and quantitative genetic were poorly correlated, indicating that molecular marker data may be of limited value to predict the evolutionary potential of populations of S. vulgaris.     

The effect of SNP discovery method and sample size on estimation of population genetic data for Chinese and Indian rhesus macaques (<Emphasis Type="Italic">Macaca mulatta</Emphasis>)

This study was designed to address issues regarding sample size and marker location that have arisen from the discovery of SNPs in the genomes of poorly characterized primate species and the application of these markers to the study of primate population genetics. We predict the effect of discovery sample size on the probability of discovering both rare and common SNPs and then compare this prediction with the proportion of common and rare SNPs discovered when different numbers of individuals are sequenced. Second, we examine the effect of genomic region on estimates of common population genetic data, comparing markers from both coding and non-coding regions of the rhesus macaque genome and the population genetic data calculated from these markers, to measure the degree and direction of bias introduced by SNPs located in coding versus non-coding regions of the genome. We found that both discovery sample size and genomic region surveyed affect SNP marker attributes and population genetic estimates, even when these are calculated from an expanded data set containing more individuals than the original discovery data set. Although none of the SNP detection methods or genomic regions tested in this study was completely uninformative, these results show that each has a different kind of genetic variation that is suitable for different purposes, and each introduces specific types of bias. Given that each SNP marker has an individual evolutionary history, we calculated that the most complete and unbiased representation of the genetic diversity present in the individual can be obtained by incorporating at least 10 individuals into the discovery sample set, to ensure the discovery of both common and rare polymorphisms.     

Biodiversity of 52 chicken populations assessed by microsatellite typing of DNA pools

In a project on the biodiversity of chickens funded by the European Commission (EC), eight laboratories collaborated to assess the genetic variation within and between 52 populations from a wide range of chicken types. Twenty-two di-nucleotide microsatellite markers were used to genotype DNA pools of 50 birds from each population. The polymorphism measures for the average, the least polymorphic population (inbred C line) and the most polymorphic population (Gallus gallus spadiceus) were, respectively, as follows: number of alleles per locus, per population: 3.5, 1.3 and 5.2; average gene diversity across markers: 0.47, 0.05 and 0.64; and proportion of polymorphic markers: 0.91, 0.25 and 1.0. These were in good agreement with the breeding history of the populations. For instance, unselected populations were found to be more polymorphic than selected breeds such as layers. Thus DNA pools are effective in the preliminary assessment of genetic variation of populations and markers. Mean genetic distance indicates the extent to which a given population shares its genetic diversity with that of the whole tested gene pool and is a useful criterion for conservation of diversity. The distribution of population-specific (private) alleles and the amount of genetic variation shared among populations supports the hypothesis that the red jungle fowl is the main progenitor of the domesticated chicken.     

籼粳杂种双单倍体的配子选择

对典型的灿与粳稻杂种,窄叶青８号／京系１７Ｆ１花药进行培养获得的１３２个双单倍体的形态特性、同工酶与ＲＦＬＰ标记的分离与重组进行了考察分析,研究是否存在配子选择问题。结果表明：（１）对４个重要数量性状和６个涉及籼、粳特征的形态指数进行考察所获数据均为连续分布,并呈正态曲线;（２）用８种同工酶对５２个ＤＨ系分析结果表明,只有２种同工酶显著偏离期望的１：１比率,而灿与粳的总基因型比率相近;（３）应用１６７个ＲＦＬＰ标记对１３２个ＤＨ系进行的分析发现,有３６％标记发生偏分离,但偏籼与偏粳的比率相近,两个亲本基因组在ＤＨ群体中所占比率相同（各５０％）,各种基因组成呈正态分布。综上所述,本研究虽观察到一些轻微偏分离现象,但籼粳基因基本上随机分离与重组,等位基因总频率未偏离１：１比率。     

Quality Control of Genotypes Using Heritability Estimates of Gene Content at the Marker

Quality control filtering of single-nucleotide polymorphisms (SNPs) is a key step when analyzing genomic data. Here we present a practical method to identify low-quality SNPs, meaning markers whose genotypes are wrongly assigned for a large proportion of individuals, by estimating the heritability of gene content at each marker, where gene content is the number of copies of a particular reference allele in a genotype of an animal (0, 1, or 2). If there is no mutation at the marker, gene content has an additive heritability of 1 by construction. The method uses restricted maximum likelihood (REML) to estimate heritability of gene content at each SNP and also builds a likelihood-ratio test statistic to test for zero error variance in genotyping. As a by-product, estimates of the allele frequencies of markers at the base population are obtained. Using simulated data with 10% permutation error (4% actual error) in genotyping, the method had a specificity of 0.96 (4% of correct markers are rejected) and a sensitivity of 0.99 (1% of wrong markers are accepted) if markers with heritability lower than 0.975 are discarded. Checking of Mendelian errors resulted in a lower sensitivity (0.84) for the same simulation. The proposed method is further illustrated with a real data set with genotypes from 3534 animals genotyped for 50,433 markers from the Illumina PorcineSNP60 chip and a pedigree of 6473 individuals; those markers underwent very little quality control. A total of 4099 markers with P-values lower than 0.01 were discarded based on our method, with associated estimates of heritability as low as 0.12. Contrary to other techniques, our method uses all information in the population simultaneously, can be used in any population with markers and pedigree recordings, and is simple to implement using standard software for REML estimation. Scripts for its use are provided.     

Sequence-tagged-sites (STSs) of cDNA clones in Cryptomeria japonica and their evaluation as molecular markers in conifers

 We have generated 66 sequence-tagged-site (STS) markers from cDNA clones of Cryptomeria japonica, and 60% of them have already been mapped into C. japonica linkage groups. All of the STS markers showed a single fragment following polymerase chain reaction (PCR) amplification. We investigated by polymorphism of these STS markers in a mapped F₂ population and 15 plus trees by means of a restriction endonuclease analysis. Polymorphism levels were 10.6% and 22.7% in the F₂ population and the 15 plus trees, respectively. PCR amplification levels of the 66 STS markers in 14 conifer species varied depending on their genetic relationship with C. japonica. Taxodium, which is closely related to C. japonica, had the most amplifications (31.82%), followed by Sequoiadendron giganteum, which is of the same family. The average proportion of PCR amplifications in each family gradually declined in the following order: from Taxodiaceae to Cuppresaceae, Sciadopityaceae, Pinaceae, and Taxaceae. These results are in general agreement with a molecular phylogenetic relationship based on chloroplast DNA. The 66 STS markers will be useful as on anchor point for genome mapping and population genetics, and some of them will also be useful when studying other conifers. Received: 8 September 1996 / Accepted: 8 November 1996     

Environmental and other factors influencing fertility in Sri Lanka

The study utilizes the complete tabulations from the 1971 census and data from the Registrar General's and Meteorological departments for the same year. It applies a multiple regression analysis of the total fertility rate for each district for 1971 on a series of environmental variables. These include the proportion of the rural population by district to the all-island rural population; the district mean annual rainfall and the proportion of the district population employed in the major rural occupations to the total employed. To these 3 environmental variables were added 2 socioeconomic variables: the proportion of the population aged 10+ years and who are literate, and infant mortality. Results indicate that the 5 independent variables are clearly intercorrelated. The multiple regression analysis shows that the 5 variables together account for 76.4% of the total variation in district total fertility rate. It is argued that this and other studies undertaken previously provide useful pointers to the type of variable to be considered in any policy aimed at population control and indicate where major efforts should be directed.     

Comparing the distribution of RAPD and RFLP markers in a high density linkage map of sugar beet.

The distribution of RAPD markers was compared with that of RFLP markers in a high density linkage map of sugar beet. The same mapping population of 161 F2 individuals was used to generate all the marker data. The total map comprises 160 RAPD and 248 RFLP markers covering 508 cM. Both the RAPD and the RFLP markers show a high degree of clustering over the nine linkage groups. The pattern is compatible with a strong distal localization of recombination in the sugar beet. It leads generally to one major cluster of markers in the centre of each linkage group. In regions of high marker density, dominant RAPD markers present in either linkage phase and codominant RFLP markers are subclustered relative to each other. This phenomenon is shown to be attributable to: (i) effects of the mapping procedure when dominant and codominant data are combined, (ii) effects of the mapping procedure when dominant data in both linkage phases are combined, and (iii) genuine differences in the way RAPD and RFLP markers are recruited.     

A genome-scan method to identify selected loci appropriate for both dominant and codominant markers: a Bayesian perspective

Identifying loci under natural selection from genomic surveys is of great interest in different research areas. Commonly used methods to separate neutral effects from adaptive effects are based on locus-specific population differentiation coefficients to identify outliers. Here we extend such an approach to estimate directly the probability that each locus is subject to selection using a Bayesian method. We also extend it to allow the use of dominant markers like AFLPs. It has been shown that this model is robust to complex demographic scenarios for neutral genetic differentiation. Here we show that the inclusion of isolated populations that underwent a strong bottleneck can lead to a high rate of false positives. Nevertheless, we demonstrate that it is possible to avoid them by carefully choosing the populations that should be included in the analysis. We analyze two previously published data sets: a human data set of codominant markers and a Littorina saxatilis data set of dominant markers. We also perform a detailed sensitivity study to compare the power of the method using amplified fragment length polymorphism (AFLP), SNP, and microsatellite markers. The method has been implemented in a new software available at our website (http://www-leca.ujf-grenoble.fr/logiciels.htm).     

Comparative microsatellite linkage analysis and genetic structure of two populations of F6 lines derived from Lycopersicon pimpinellifolium and L. cheesmanii

A population of recombinant inbred lines (RILs) has several advantages over its F₂ population counterpart with respect to quantitative trait loci (QTLs) and genomic studies. The objective of the investigation reported here was the comparative characterization by simple sequence repeat (SSR) and sequence characterized amplified region (SCAR) markers of two populations of F₆ lines derived from Lycopersicon pimpinellifolium (P population, consisting of 142 lines) and L. cheesmanii (C population, consisting of 115 lines) and sharing the female parent, L. esculentum var. cerasiforme. Almost the same percentage of polymorphic markers was found for each population although a different set of markers was involved. The proportion of SSR primer pairs (93 in total) that resulted in polymorphism for the main band was larger (55–56%) than for SCAR ones (13–16%). The C population showed the largest proportion of markers with zygotic and gametic segregation distortion, which is in agreement with the larger genetic distance reported between L. esculentum and L. cheesmanii than with the former and L. pimpinellifolium. Zygotic distortion corresponded primarily to an excess of heterozygotes in both populations, suggesting that the increment of homozygosity was the main factor limiting viability/self-fertility of the lines. Despite both populations sharing the female parent, P alleles were slightly favored in the P population while E alleles were the most frequently fixed in the C population. A linkage map for each population was obtained, with the average distances between consecutive markers being 3.8 cM or 3.4 cM depending on the population. Discrepancy between the maps for the location of only four markers on chromosomes 3, 6 and 10 was observed. Two possible causes of this discrepancy were investigated and can not be discarded: (1) the presence of duplicated markers and (2) segregation distortion caused by the selective advantage of gametes carrying one of the two alleles. This marker characterization of both populations will continue and will enable the comparative QTLs and candidate gene analysis of complex traits towards a more efficient utilization of genetic resources and breeding strategies.     

Integration of trinucleotide microsatellites into a linkage map of Citrus

 We report the successful assignment of the first seven microsatellite markers to the Citrus RFLP and isozyme map. A total of 14 microsatellite primer pairs were developed and tested for amplification and product-length polymorphism within a population of plants previously used for linkage-map construction. In each case, the successfully assigned microsatellite mapped to the termini of a different linkage group indicating a widespread distribution throughout the genome. Analysis of allele segregation revealed that two of nine microsatellites displayed a significant deviation from expected ratios (P>0.5). This was compared with other marker types within Citrus and a similar proportion of skewed loci was also found to be present. The analysis of two markers was complicated by the non-amplification of an inherited null allele within the mapping population. The successful integration of microsatellites into the genetic map of Citrus demonstrates the utility of this marker type for genetic analysis within wide intergeneric plant crosses. Received: 16 September 1996 / Accepted: 18 October 1996     

Statistical confidence for likelihood-based paternity inference in natural populations

Paternity inference using highly polymorphic codominant markers is becoming common in the study of natural populations. However, multiple males are often found to be genetically compatible with each offspring tested, even when the probability of excluding an unrelated male is high. While various methods exist for evaluating the likelihood of paternity of each nonexcluded male, interpreting these likelihoods has hitherto been difficult, and no method takes account of the incomplete sampling and error-prone genetic data typical of large-scale studies of natural systems. We derive likelihood ratios for paternity inference with codominant markers taking account of typing error, and define a statistic Δ for resolving paternity. Using allele frequencies from the study population in question, a simulation program generates criteria for Δ that permit assignment of paternity to the most likely male with a known level of statistical confidence. The simulation takes account of the number of candidate males, the proportion of males that are sampled and gaps and errors in genetic data. We explore the potentially confounding effect of relatives and show that the method is robust to their presence under commonly encountered conditions. The method is demonstrated using genetic data from the intensively studied red deer ( Cervus elaphus ) population on the island of Rum, Scotland. The Windows-based computer program, CERVUS , described in this study is available from the authors. CERVUS can be used to calculate allele frequencies, run simulations and perform parentage analysis using data from all types of codominant markers.     

AFLP markers closely linked to a major gene essential for nucellar embryony (apomixis<Emphasis Type="Italic">)</Emphasis> in <Emphasis Type="Italic">Citrus maxima</Emphasis> × <Emphasis Type="Italic">Poncirus trifoliata</Emphasis>

Some citrus varieties express a form of apomixis termed nucellar embryony in which the adventive embryos develop from nucellus tissue surrounding the embryo sac. This trait results in many seeds containing multiple embryos (polyembryony). Inheritance of the frequency of polyembryony was studied in 88 progeny from a cross of Citrus maxima (monoembryonic) × Poncirus trifoliata (polyembryonic). The frequency of polyembryonic seed produced by each progeny was determined by scoring 100–500 seeds for the number of seedlings to emerge from each seed. Two groups of eight individuals from each extreme of the population were chosen for bulked segregant analysis with amplified fragment length polymorphism markers amplified with 256 primer combinations. Candidate markers identified in the bulks as linked to the trait were tested on the 32 individuals used to create the bulks and then on the remaining plants in the population. Five candidate markers tightly linked to polyembryony in P. trifoliata were identified. Specific marker alleles were present in nearly all progeny that produced polyembryonic seed, and alternate alleles were present in nearly all progeny that produced only monoembryonic seed. The region defined by these markers very likely contains a gene that is essential for the production of polyembryonic seeds by apomixis, but also shows segregation distortion. The proportion of polyembryonic seeds varied widely among the hybrid progeny, probably due to other genes. Scoring 119 progeny of a P. trifoliata selfed population for the closely linked markers and the proportion of polyembryonic seeds confirmed close linkage between these markers and polyembryony.     

Rapid Assessment of Genetic Ancestry in Populations of Unknown Origin by Genome-Wide Genotyping of Pooled Samples

As we move forward from the current generation of genome-wide association (GWA) studies, additional cohorts of different ancestries will be studied to increase power, fine map association signals, and generalize association results to additional populations. Knowledge of genetic ancestry as well as population substructure will become increasingly important for GWA studies in populations of unknown ancestry. Here we propose genotyping pooled DNA samples using genome-wide SNP arrays as a viable option to efficiently and inexpensively estimate admixture proportion and identify ancestry informative markers (AIMs) in populations of unknown origin. We constructed DNA pools from African American, Native Hawaiian, Latina, and Jamaican samples and genotyped them using the Affymetrix 6.0 array. Aided by individual genotype data from the African American cohort, we established quality control filters to remove poorly performing SNPs and estimated allele frequencies for the remaining SNPs in each panel. We then applied a regression-based method to estimate the proportion of admixture in each cohort using the allele frequencies estimated from pooling and populations from the International HapMap Consortium as reference panels, and identified AIMs unique to each population. In this study, we demonstrated that genotyping pooled DNA samples yields estimates of admixture proportion that are both consistent with our knowledge of population history and similar to those obtained by genotyping known AIMs. Furthermore, through validation by individual genotyping, we demonstrated that pooling is quite effective for identifying SNPs with large allele frequency differences (i.e., AIMs) and that these AIMs are able to differentiate two closely related populations (HapMap JPT and CHB).     

combi.pl: a computer program to combine data sets with inconsistent microsatellite marker allele size information

Combining two data sets with allele information from overlapping microsatellite markers is often desirable, particularly in population genetic studies where a substantial body of published data exists. When genotyping is performed in different laboratories, allele size calling may not be presumed to be consistent. Our approach solves this problem by assigning allele sizes across studies using maximum-likelihood theory. Using data overlaps in samples and markers, allele shifts between two studies are calculated for each overlapping marker and a single file containing allele frequencies of consistent alleles is produced. The program (combi.pl) is written in PERL and available at http://data40.uni-tz.gwdg.de/~htaeube.     

Restriction fragment length polymorphisms of type I collagen locus 2 (COL1A2) in two communities of African ancestry and other mixed populations of northwestern Ecuador

Restriction fragment length polymorphisms are good anthropological markers for discriminating geographically distinct populations at both the allele and the haplotype level. Two communities of African ancestry and ladinos, mestizos, and mulattoes living in the Esmeraldas province in northwestern Ecuador were analyzed for three RFLPs (EcoRI, RsaI, and MspI) of the COL1A2 gene. Also, the same markers were studied in a population sample from Spain to compare the allele and haplotype frequencies of the Esmeraldas populations with those of their representative European parental population. Data for the native American and sub-Saharan African founder components were available from the literature. No significant levels of differentiation between the two African Ecuadoran communities emerged from either the frequency analysis of each single marker and all three RFLP markers together or from the AMOVA. The ladinos and mestizos also showed a rather similar distribution of allele and haplotype frequencies, confirming that the two ethnic terms do not correspond to genetically different populations. The comparison with the supposed founding European, sub-Saharan African, and native American populations indicated a large presence of African genes in the gene pool of both communities, with a higher proportion of the Amerindian component in Viche than in Rio Cayapas. The present findings confirm the previous genetic admixture estimates based on nuclear and mitochondrial DNA markers and the demographic data.     

鳞翅目昆虫基因组中微卫星DNA的特征以及对其分离的影响

本文根据我们对鳞翅目昆虫棉铃虫和松毛虫以及其它动物 (筏蜘蛛、朱、鳕鱼和飞蝗 )的微卫星富集性基因组DNA文库的筛选和分析结果 ,结合其它实验室已发表的资料 ,对鳞翅目昆虫基因组中微卫星DNA的丰度和结构特点进行了较为系统的分析。结果表明 :与其它类群相比 ,尽管鳞翅目昆虫物种间存在差异 ,但其基因组中存在明显偏多的侧翼序列重复的、以多拷贝形式存在的微卫星位点 ,且其中相当一部分以基因家族的形式存在。微卫星DNA家族通常可以在序列分析阶段被识别出来 ,但很多多拷贝位点只有通过一系列后续分析才能被检查出来。这应是鳞翅目昆虫中微卫星位点的优化率相对偏低的主要原因。棉铃虫和松毛虫基因组中三相重复微卫星丰度相对较高 ,从而从某种程度上补偿了这些物种微卫星分离过程中因丰度低、多拷贝位点比例高所带来的困难。棉铃虫微卫星DNA家族侧翼序列中多聚T/A序列的存在表明 ,逆转录转座或逆转录侵染可能是在基因组中形成多拷贝微卫星位点和微卫星DNA家族的重要机制之一     

Detection of agriculturally important QTLs in chickens and analysis of the factors affecting genotyping strategy

Three single cross populations were generated in order to analyze factors affecting the ability to detect true linkage with minimum false positive or false negative associations, and to detect associations between markers and quantitative traits. The three populations are: (1) a broiler x broiler cross of a single sire and 34 dams, resulting in 266 progeny; (2) a broiler x broiler cross of a single sire and 41 dams resulting in 360 progeny; and (3) a broiler x layer cross of a single sire with 56 dams resulting in 1180 progeny. Based on these three resource populations we show that: a) gradient selective genotyping was more effective than the random selective genotyping; b) selective genotyping was significant at a selected proportion less than 62% of the cumulative truncation point; c) as few as 10% of selected individuals (5% of each of the two tails) were sufficient to show significant association between markers and phenotypes; d) a gradient slices approach was more powerful than using replicates of the extreme groups; and e) in resource populations resulting from crosses between lines of different backgrounds, most of the microsatellite markers used are polymorphic. We also used simulation to test factors affecting power to detect true associations between markers and traits that are hard to detect in experimental resource populations. Using defined populations in the simulation, we concluded that the following guidelines provide reliable detection of linked QTLs: 1) the resource population size should be larger than 100; 2) a QTL effect larger than 0.4 SD is detectable with a reasonable number of markers (>100) and resource population size (>200 subjects); 3) the DNA pool from each tail of the trait distribution should contain at least 10% of the resource family; 4) each of the two DNA pools should include more than 35 individuals. Some of these guidelines that were deduced from the simulation analysis have been confirmed in the experimental part of this study.     

Development of microsatellite markers and construction of genetic map in rice blast pathogen Magnaporthe grisea

Magnaporthe grisea is the most destructive fungal pathogen of rice and a model organism for studying plant-pathogen interaction. Molecular markers and genetic maps are useful tools for genetic studies. In this study, based on the released genome sequence data of M. grisea, we investigated 446 simple sequence repeat (SSR) loci and developed 313 SSR markers, which showed polymorphisms among nine isolates from rice (including a laboratory strain 2539). The number of alleles of each marker ranged 2-9 with an average of 3.3. The polymorphic information content (PIC) of each marker ranged 0.20-0.89 with an average of 0.53. Using a population derived from a cross between isolates Guy11 and 2539, we constructed a genetic map of M. grisea consisting of 176 SSR markers. The map covers a total length of 1247 cM, equivalent to a physical length of about 35.0 Mb or 93% of the genome, with an average distance of 7.1cM between adjacent markers. A web-based database of the SSR markers and the genetic map was established (http://ibi.zju.edu.cn/pgl/MGM/index.html).