Comparison of genetic differentiation at marker loci and quantitative traits

The comparison of the degree of differentiation in neutral marker loci and genes coding quantitative traits with standardized and equivalent measures of genetic differentiation (F_ST and Q_ST, respectively) can provide insights into two important but seldom explored questions in evolutionary genetics: (i) what is the relative importance of random genetic drift and directional natural selection as causes of population differentiation in quantitative traits, and (ii) does the degree of divergence in neutral marker loci predict the degree of divergence in genes coding quantitative traits? Examination of data from 18 independent studies of plants and animals using both standard statistical and meta‐analytical methods revealed a number of interesting points. First, the degree of differentiation in quantitative traits (Q_ST) typically exceeds that observed in neutral marker genes (F_ST), suggesting a prominent role for natural selection in accounting for patterns of quantitative trait differentiation among contemporary populations. Second, the F_ST – Q_ST difference is more pronounced for allozyme markers and morphological traits, than for other kinds of molecular markers and life‐history traits. Third, very few studies reveal situations were Q_ST < F_ST, suggesting that selection pressures, and hence optimal phenotypes, in different populations of the same species are unlikely to be often similar. Fourth, there is a strong correlation between Q_ST and F_ST indices across the different studies for allozyme (r=0.81), microsatellite (r=0.87) and combined (r=0.75) marker data, suggesting that the degree of genetic differentiation in neutral marker loci is closely predictive of the degree of differentiation in loci coding quantitative traits. However, these interpretations are subject to a number of assumptions about the data and methods used to derive the estimates of population differentiation in the two sets of traits.     

Genome scans for the contemporary response to selection in quantitative traits

Genome scans have been an important approach for discovering historical signatures of selection in both model and nonmodel species. An exciting new experimental design for genome scans is to measure the change in allele frequency before and after contemporary selection within a generation, from a single population. The most widely‐used methods, however, have two major limitations: they are based on testing one locus at a time, and they only have power to uncover loci that have evolved under relatively strong selection. On the other hand, complex quantitative traits are common in nature and are caused by several loci of small effect. Selection on a quantitative trait at the phenotypic level is predicted to be accompanied by subtle allele frequency changes in many loci that covary (a polygenic soft sweep), rather than a large, single‐effect allele (a selective sweep). In this issue of Molecular Ecology, Bourret et al. (2014) measure the contemporary response to natural selection across the genome in multiple cohorts of Atlantic salmon during their first year at sea. They introduce a multilocus framework based on groups of markers that covary in their genotypic distribution. While the traditional, single‐locus approach did not find evidence for repeated patterns of selection, the multivariate approach found that a group of covarying SNPs was selected for in different cohorts at one site. Their multilocus framework has potential to be a more fruitful approach for uncovering the genomic basis of adaptation in quantitative traits, although caution should be applied as the framework has yet to be validated with simulated data.     

The genetic differentiation at quantitative trait loci under local adaptation

Most adaptive traits are controlled by large number of genes that may all together be the targets of selection. Adaptation may thus involve multiple but not necessarily substantial allele frequency changes. This has important consequences for the detection of selected loci and implies that a quantitative genetics framework may be more appropriate than the classical 'selective sweep' paradigm. Preferred methods to detect loci involved in local adaptation are based on the detection of 'outlier' values of the allelic differentiation F(ST) . A quantitative genetics framework is adopted here to review theoretical expectations for how allelic differentiation at quantitative trait loci (F(STQ) ) relates to (i), neutral genetic differentiation (F(ST) ) and (ii), phenotypic differentiation (Q(ST) ). We identify cases where results of outlier-based methods are likely to be poor and where differentiation at selected loci conveys little information regarding local adaptation. A first case is when neutral differentiation is high, so that local adaptation does not necessitate increased differentiation. A second case is when local adaptation is reached via an increased covariance of allelic effects rather than via allele frequency changes, which is more likely under high gene flow when the number of loci is high and selection is recent. The comparison of theoretical predictions with observed data from the literature suggests that polygenic local adaptation involving only faint allele frequency changes are very likely in some species such as forest trees and for climate-related traits. Recent methodological improvements that may alleviate the weakness of F(ST) -based detection methods are presented.     

High quantitative and no molecular differentiation of a freshwater snail (Galba truncatula) between temporary and permanent water habitats

We investigate the variation in quantitative and molecular traits in the freshwater snail Galba truncatula, from permanent and temporary water habitats. Using a common garden experiment, we measured 20 quantitative traits and molecular variation using seven microsatellites in 17 populations belonging to these two habitats. We estimated trait means in each habitat. We also estimated the distributions of overall genetic quantitative variation (QST), and of molecular variation (FST), within and between habitats. Overall, we observed a lack of association between molecular and quantitative variance. Among habitats, we found QST>FST, an indication of selection for different optima. Individuals from temporary water habitat matured older, at a larger size and were less fecund than individuals from permanent water habitat. We discuss these findings in the light of several theories for life-history traits evolution.     

Technical note: quantitative measures of iris color using high resolution photographs

Our understanding of the genetic architecture of iris color is still limited. This is partly related to difficulties associated with obtaining quantitative measurements of eye color. Here we introduce a new automated method for measuring iris color using high resolution photographs. This method extracts color measurements in the CIE 1976 L*a*b* (CIELAB) color space from a 256 by 256 pixel square sampled from the 9:00 meridian of the iris. Color is defined across three dimensions: L* (the lightness coordinate), a* (the red-green coordinate), and b* (the blue-yellow coordinate). We applied this method to a sample of individuals of diverse ancestry (East Asian, European and South Asian) that was genotyped for the HERC2 rs12913832 polymorphism, which is strongly associated with blue eye color. We identified substantial variation in the CIELAB color space, not only in the European sample, but also in the East Asian and South Asian samples. As expected, rs12913832 was significantly associated with quantitative iris color measurements in subjects of European ancestry. However, this SNP was also strongly associated with iris color in the South Asian sample, although there were no participants with blue irides in this sample. The usefulness of this method is not restricted only to the study of iris pigmentation. High-resolution pictures of the iris will also make it possible to study the genetic variation involved in iris textural patterns, which show substantial heritability in human populations.     

The distribution of the effects of genes affecting quantitative traits in livestock

Meta-analysis of information from quantitative trait loci (QTL) mapping experiments was used to derive distributions of the effects of genes affecting quantitative traits. The two limitations of such information, that QTL effects as reported include experimental error, and that mapping experiments can only detect QTL above a certain size, were accounted for. Data from pig and dairy mapping experiments were used. Gamma distributions of QTL effects were fitted with maximum likelihood. The derived distributions were moderately leptokurtic, consistent with many genes of small effect and few of large effect. Seventeen percent and 35% of the leading QTL explained 90% of the genetic variance for the dairy and pig distributions respectively. The number of segregating genes affecting a quantitative trait in dairy populations was predicted assuming genes affecting a quantitative trait were neutral with respect to fitness. Between 50 and 100 genes were predicted, depending on the effective population size assumed. As data for the analysis included no QTL of small effect, the ability to estimate the number of QTL of small effect must inevitably be weak. It may be that there are more QTL of small effect than predicted by our gamma distributions. Nevertheless, the distributions have important implications for QTL mapping experiments and Marker Assisted Selection (MAS). Powerful mapping experiments, able to detect QTL of 0.1σ_p, will be required to detect enough QTL to explain 90% the genetic variance for a quantitative trait.     

Mapping quantitative trait loci regulating chicken body composition traits

Genome scans were conducted on an F₂ resource population derived from intercross of the White Plymouth Rock with the Silkies Fowl to detect QTL affecting chicken body composition traits. The population was genotyped with 129 microsatellite markers and phenotyped for 12 body composition traits on 238 F₂ individuals from 15 full-sib families. In total, 21 genome-wide QTL were found to be responsible for 11 traits, including two newly studied traits of proventriculus weight and shank girth. Three QTL were genome-wide significant: at 499 c m on GGA1 (explained 3.6% of phenotypic variance, P  < 0.01) and 51 c m on GGA5 (explained 3.3% of phenotypic variance, P  < 0.05) for the shank & claw weight and 502 c m on GGA1 (explained 1.4% of phenotypic variance, P  < 0.05) for wing weight. The QTL on GGA1 seemed to have pleiotropic effects, also affecting gizzard weight at 490 c m , shank girth at 489 c m and intestine length at 481 c m . It is suggested that further efforts be made to understand the possible pleiotropic effects of the QTL on GGA1 and that on GGA5 for two shank-related traits.     

Mapping of quantitative trait loci for clinical–chemical traits in swine

Clinical–chemical traits are diagnostic parameters essential for characterization of health and disease in veterinary practice. The traits show significant variability and are under genetic control, but little is known about the fundamental genetic architecture of this variability, especially in swine. We have identified QTL for alkaline phosphatase (ALP), lactate (LAC), bilirubin (BIL), creatinine (CRE) and ionized sodium (Na⁺), potassium (K⁺) and calcium (Ca⁺⁺) from the serum of 139 F₂ pigs from a Meishan/Pietrain family before and after challenge with Sarcocystis miescheriana , a protozoan parasite of muscle. After infection, the pigs passed through three stages representing acute disease, subclinical disease and chronic disease. Forty-two QTL influencing clinical–chemical traits during these different stages were identified on 15 chromosomes. Eleven of the QTL were significant on a genome-wide level; 31 QTL were chromosome-wide significant. QTL showed specific health/disease patterns with respect to the baseline values of the traits as well as the values obtained through the different stages of disease. QTL influencing different traits at different times were found primarily on chromosomes 1, 3, 7 and 14. The most prominent QTL for the investigated clinical–chemical traits mapped to SSC3 and 7. Baseline traits of ALP, LAC, BIL, Ca⁺⁺ and K⁺ were influenced by QTL regions on SSC3, 6, 7, 8 and 13. Single QTL explained up to 21.7% of F₂ phenotypic variance. Our analysis confirms that variation of clinical–chemical traits is associated with multiple chromosomal regions.     

Tests for major genes affecting quantitative traits in wild radish,Raphanus raphanistrum

The number of gene loci coding for quantitative traits is an important issue in genetics. However, there are still very few empirical data on this point, especially in natural populations. I tested for major gene effects on ten quantitative traits in wild radish, using an indirect method based on the patterns of family means and within and between family variances for traits. This method should reveal whether a single locus is responsible for most of the variation in a trait. Eight of the traits measured were morphological dimensions of leaves and flowers; no strong evidence for major gene effects on these traits was found. In contrast, evidence for major gene effects was found in the other two traits, emergence time and flowering time.     

Mapping of quantitative trait loci for carcass traits in a Japanese Black (Wagyu) cattle population

To detect quantitative trait loci (QTL) that influence economically important traits in a purebred Japanese Black cattle population, we performed a preliminary genome-wide scan using 187 microsatellite markers across a paternal half-sib family composed of 258 offspring. We located six QTL at the 1% chromosome-wise level on bovine chromosomes (BTA) 4, 6, 13, 14 and 21. A second screen of these six QTL regions using 138 additional paternal offspring half-sib from the same sire, provided further support for five QTL: carcass weight on BTA14 (22-39 cM), one for rib thickness on BTA6 (27-58 cM) and three for beef marbling score (BMS) on BTA4 (59-67 cM), BTA6 (68-89 cM) and BTA21 (75-84 cM). The location of QTL for subcutaneous fat thickness on BTA13 was not supported by the second screen (P > 0.05). We determined that the combined contribution of the three QTLs for BMS was 10.1% of the total variance. The combined phenotypic average of these three Q was significantly different (P < 0.001) from those of other allele combinations. Analysis of additional half-sib families will be necessary to confirm these QTL.     

甘蔗品种资源的表型遗传多样性

为了提高甘蔗品种资源的利用效率, 为甘蔗遗传育种亲本的选择、杂交组合配制及核心种质构建提供理论指导, 我们利用17个质量性状和5个数量性状分析了来自13个国家共20个地区的1,160份甘蔗品种资源的遗传变异、遗传结构和遗传距离等。结果表明: 5个数量性状在不同来源地品种群体之间的变异系数值(CV)变化较大, 说明不同来源地品种群体的数量遗传变异有较大差异, 其中来自海南的品种其群体遗传变异最丰富。质量性状的Shannon-Wiener多样性指数表明, 来自美国的品种群体遗传多样性最高, 其次为中国台湾, 第三为澳大利亚, 说明上述3个地区的甘蔗种质创新比较活跃, 在遗传育种中使用了更多表型性状多样化的亲本。不同来源地品种其群体的遗传分化系数(Gst)和基因流(Nm)显示甘蔗品种群体表型性状的遗传变异主要来自来源地内部, 且不同来源地品种群体之间存在较大的基因交流。遗传距离和UPGMA聚类分析结果表明, 各来源地品种群体之间遗传距离在0.0261–0.2945之间, 其中以福建和广东的品种最为相似, 其次为古巴和美国、广西和云南、澳大利亚和菲律宾、江西和四川、巴西和法国, 说明上述地区在杂交亲本的选择上比较相近。鉴于此, 在遗传育种中应加大利用具有丰富遗传多样性的品种材料并尽量避免选择同一组的品种相互杂交, 同时对于与其他来源地品种群体遗传距离较远的墨西哥品种群体在亲本选配时应给予更多关注。     

Molecular and quantitative genetic differentiation across Europe in yellow dung flies

Relating geographic variation in quantitative traits to underlying population structure is crucial for understanding processes driving population differentiation, isolation and ultimately speciation. Our study represents a comprehensive population genetic survey of the yellow dung fly Scathophaga stercoraria, an important model organism for evolutionary and ecological studies, over a broad geographic scale across Europe (10 populations from the Swiss Alps to Iceland). We simultaneously assessed differentiation in five quantitative traits (body size, development time, growth rate, proportion of diapausing individuals and duration of diapause), to compare differentiation in neutral marker loci (F(ST)) to that of quantitative traits (Q(ST)). Despite long distances and uninhabitable areas between sampled populations, population structuring was very low but significant (F(ST) = 0.007, 13 microsatellite markers; F(ST) = 0.012, three allozyme markers; F(ST) = 0.007, markers combined). However, only two populations (Iceland and Sweden) showed significant allelic differentiation to all other populations. We estimated high levels of gene flow [effective number of migrants (Nm) = 6.2], there was no isolation by distance, and no indication of past genetic bottlenecks (i.e. founder events) and associated loss of genetic diversity in any northern or island population. In contrast to the low population structure, quantitative traits were strongly genetically differentiated among populations, following latitudinal clines, suggesting that selection is responsible for life history differentiation in yellow dung flies across Europe.     

藏东木本植物群落功能性状分布与环境的关系

研究植物群落功能性状间的相关关系及其对环境变化的响应,能够有效揭示植物功能性状的权衡模式及其对环境的适应策略。藏东昌都地区位于横断山脉西北部,复杂气候地貌孕育了丰富的植物资源,是青藏高原森林灌丛生态系统主要组分和国际生物多样性保护的热点地区。以藏东森林灌丛群落优势木本植物为研究对象,在大量野外调查基础上,采用相关分析、主成分分析、线性回归和方差分析等方法,研究了该区域植物功能性状间的相关关系、功能性状对环境变量的响应规律以及功能性状的变异来源。结果表明：(1)藏东木本植物表现出适应高寒环境的性状权衡模式,即：比叶面积、叶体积较小而叶干物质含量较大,叶磷含量和叶钾含量协同变化;(2)海拔和气候变量共同驱动着藏东木本植物功能性状的变化,并且藏东木本植物倾向于采取“高投入—慢回报”提高御寒能力的保守型适应策略;(3)海拔是影响藏东植物功能性状变异最显著的环境变量,种间变异在藏东植物群落功能性状随环境变化中起主要作用。研究结果揭示了藏东木本植物功能性状的权衡模式及其对高寒环境的适应策略,有助于加深对藏东自然植物资源分布规律和生态功能的认识,为区域生态系统功能和生物多样性保护提供科学依据。     

A genome scan for quantitative trait loci affecting three ear traits in a White Duroc × Chinese Erhualian resource population

Chinese Erhualian pigs have larger and floppier ears compared with White Duroc pigs (small, half- or fully-pricked ears). To identify quantitative trait loci (QTL) for ear weight and area as well as erectness, a genome-wide scan with 194 microsatellites was performed in a White Duroc × Chinese Erhualian resource population (>1000 F₂ animals). Twenty-three genome-wide significant QTL and 12 suggestive QTL were identified. All QTL for ear erectness and size detected in two previous studies, bar two on SSC6 and 9, were confirmed here. The 1% genome-wide significant QTL at 70 cM on SSC5 and at 58 cM on SSC7 have profound and pleiotropic effects on the three ear traits, with Erhualian alleles increasing weight and area but decreasing erectness. Notably, the 95% confidence interval of the QTL for weight and area on SSC7 spanned only 3 cM. New QTL reaching 1% genome-wide significance were found on SSC8 (at 37 cM) for all three ear traits, on SSC4 and 16 for weight and area, and on SSCX for area. Unexpectedly, Erhualian alleles at these loci were associated with lighter and smaller or erect ear. Some new suggestive QTL were also found on other chromosome regions. Almost all the QTL for weight and area had essentially additive effects, while the QTL for erectness on SSC2, 5 and 7 showed not only additive effects but also partial dominance effects of Erhualian alleles. The two most significant QTL on SSC7 and SSC5 could be promising targets for fine mapping and identification of the causative mutations.     

Genetic and molecular dissection of quantitative traits in rice

Recent progress in the generation of a molecular genetic map and markers for rice has made possible a new phase of mapping individual genes associated with complex traits. This type of analysis is often referred to as quantitative trait locus (QTL) analysis. Increasing numbers of QTL analyses are providing enormous amounts of information about QTLs, such as the numbers of loci involved, their chromosomal locations and gene effects. Clarification of genetic bases of complex traits has a big impact not only on fundamental research on rice plant development, but it also has practical benefits for rice breeding. In this review, we summarize recent progress of QTL analysis of several complex traits in rice. A strategy for positional cloning of genes at QTLs is also discussed.     

RST Calc: a collection of computer programs for calculating estimates of genetic differentiation from microsatellite data and determining their significance

R _ST, an analogue of F _ST, provides a convenient approach for estimating levels of genetic differentiation from microsatellite data. This paper examines current approaches for calculating estimates of R _ST and suggests a weighting scheme based on the transformation of allele sizes at loci across data sets. Combined within an analysis of variance framework this scheme yields an estimator of R _ST analogous to the θ estimator of F _ST. Software for the IBM-PC is described which carries out such calculations and assesses the significance of R _ST or Nm estimates using bootstrap and permutation tests.     

Population structure in the barn swallow,Hirundo rustica: a comparison between neutral DNA markers and quantitative traits

Local adaptation to variable environments can generate clinal variation in morphological traits. Alternatively, similar patterns of clinal variation may be generated simply as a result of genetic drift/migration balance. Teasing apart these different processes is a continuing focus in evolutionary ecology. We compare genetic differentiation at molecular loci and quantitative traits to analyse the effect of these different processes in a morphological latitudinal cline of the barn swallow, Hirundo rustica, breeding across Europe. The results obtained show no structuring at neutral microsatellite loci, which contrasts with positive structuring at five quantitative morphometric traits. This supports the hypothesis that the observed morphometric cline in barn swallows is the result of selection acting in a spatially heterogeneous environment. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 99 , 306–314.