A MARKER-BASED METHOD FOR INFERENCES ABOUT QUANTITATIVE INHERITANCE IN NATURAL POPULATIONS

A marker-based method for studying quantitative genetic characters in natural populations is presented and evaluated. The method involves regressing quantitative trait similarity on marker-estimated relatedness between individuals. A procedure is first given for estimating the narrow sense heritability and additive genetic correlations among traits, incorporating shared environments. Estimation of the actual variance of relatedness is required for heritability, but not for genetic correlations. The approach is then extended to include isolation by distance of environments, dominance, and shared levels of inbreeding. Investigations of statistical properties show that good estimates do not require great marker polymorphism, but rather require significant variation of actual relatedness; optimal allocation generally favors sampling many individuals at the expense of assaying fewer marker loci; when relatedness declines with physical distance, it is optimal to restrict comparisons to within a certain distance; the power to estimate shared environments and inbreeding effects is reasonable, but estimates of dominance variance may be difficult under certain patterns of relationship; and any linkage of markers to quantitative trait loci does not cause significant problems. This marker-based method makes possible studies with long-lived organisms or with organisms difficult to culture, and opens the possibility that quantitative trait expression in natural environments can be analyzed in an unmanipulative way.     

Identity-by-descent-based heritability analysis in the Northern Finland Birth Cohort

For most complex traits, only a small proportion of heritability is explained by statistically significant associations from genome-wide association studies (GWAS). In order to determine how much heritability can potentially be explained through larger GWAS, several different approaches for estimating total narrow-sense heritability from GWAS data have recently been proposed. These methods include variance components with relatedness estimates from allele-sharing, variance components with relatedness estimates from identity-by-descent (IBD) segments, and regression of phenotypic correlation on relatedness estimates from IBD segments. These methods have not previously been compared on real or simulated data. We analyze the narrow-sense heritability of nine metabolic traits in the Northern Finland Birth Cohort (NFBC) using these methods. We find substantial estimated heritability for several traits, including LDL cholesterol (54 % heritability), HDL cholesterol (46 % heritability), and fasting glucose levels (39 % heritability). Estimates of heritability from the regression-based approach are much lower than variance component estimates in these data, which may be due to the presence of strong population structure. We also investigate the accuracy of the competing approaches using simulated phenotypes based on genotype data from the NFBC. The simulation results substantiate the downward bias of the regression-based approach in the presence of population structure.     

Heritability of phenolics in Quercus laevis inferred using molecular markers.

Studies of quantitative inheritance of phenotypes do not generally encompass the range of environmental conditions to which a population may be exposed in a natural setting and are rarely conducted on long-lived species due to the time required for traditional crossing experiments. We used a marker-based method to estimate relatedness with microsatellite markers in a natural population of a long-lived oak, then used this inferred relatedness to examine quantitative genetic variation in the concentration of foliar phenolics. Estimating heritability using this method requires both significant relatedness and variance in relatedness over distance. However, this population did not show significant variance of relatedness, so only the presence of heritability, and its ranking among traits and environments, could be estimated. Seven foliar phenolics showed a significant relationship between phenotypic similarity and relatedness. The significance of this relationship varied among individual phenolic compounds, as well as by season. Genetic factors appeared to have a more measurable influence on the production of secondary compounds early in the season. After leaf expansion, covariance of relatedness and phenotypic variance appear to become less significant. Therefore heritability may vary seasonally for these traits.     

Consistency between marker- and genealogy-based heritability estimates in an experimental stand of Prosopis alba (Leguminosae)

Prosopis represents a valuable forest resource in arid and semiarid regions. Management of promising species requires information about genetic parameters, mainly the heritability (h(2)) of quantitative profitable traits. This parameter is traditionally estimated from progeny tests or half-sib analysis conducted in experimental stands. Such an approach estimates h(2) from the ratio of between-family/total phenotypic variance. These analyses are difficult to apply to natural populations of species with a long life cycle, overlapping generations, and a mixed mating system, without genealogical information. A promising alternative is the use of molecular marker information to infer relatedness between individuals and to estimate h(2) from the regression of phenotypic similarity on inferred relatedness. In the current study we compared h(2) of 13 quantitative traits estimated by these two methods in an experimental stand of P. alba, where genealogical information was available. We inferred pairwise relatedness by Ritland's method using six microsatellite loci. Relatedness and heritability estimates from molecular information were highly correlated to the values obtained from genealogical data. Although Ritland's method yields lower h(2) estimates and tends to overestimate genetic correlations between traits, this approach is useful to predict the expected relative gain of different quantitative traits under selection without genealogical information.     

The accuracy of a heritability estimator using molecular information

The heritability of a quantitative trait is a key parameter to quantify the genetic variation present in a population. Although estimates of heritability require accurate information on the genetic relationship among individuals, pedigree data is generally lacking in natural populations. Nowadays, the increasing availability of DNA markers is making possible the estimation of coancestries from neutral molecular information. In 1996, K. Ritland developed an approach to estimate heritability from the regression of the phenotypic similarity on the marker-based coancestry. We carried out simulations to analyze the accuracy of the estimates of heritability obtained by this method using information from a variable number of neutral codominant markers. Because the main application of the estimator is on populations with no family structure, such as natural populations, its accuracy was tested under this scenario. However, the method was also investigated under other scenarios, in order to test the influence of different factors (family structure, assortative mating and phenotypic selection) on the precision. Our results suggest that the main factor causing a directional bias in the estimated heritability is the presence of phenotypic selection, and that very noisy estimates are obtained in the absence of a familiar structure and for small population sizes. The estimated heritabilities from marker-based coancestries showed lower accuracy than the estimated heritabilities from genealogical coancestries. However, a large amount of bias occurred even in the most favourable situation where genealogical coancestries are known. The results also indicate that the molecular markers are more suitable to infer coancestry than inbreeding.     

Using partial genotyping to estimate the genetic and maternal determinants of adaptive traits in a progeny trial of <Emphasis Type="Italic">Fagus sylvatica</Emphasis>

Understanding the determinants of phenotypic variation is critical to evaluate the ability of traits to evolve in a changing environment. In trees, the genetic component of the phenotypic variance is most often estimated based on maternal progeny tests. However, the lack of knowledge about the paternal relatedness hampers the accurate estimation of additive genetic and maternal effects. Here, we investigate how different methods accounting for paternal relatedness allow the estimation of heritability and maternal determinants of adaptive traits in a natural population of Fagus sylvatica L., presenting non-random mating. Twelve potentially adaptive functional traits were measured in 60 maternal families in a nursery. We genotyped a subset of offspring and of all the potentially reproductive adults in the population at 13 microsatellite markers to infer paternal relationships and to estimate average relatedness within and between maternal families. This relatedness information was then used in family and animal models to estimate the components of phenotypic variance. All the studied traits displayed significant genetic variance and moderate heritability. Maternal effects were detected for the diameter increment, stem volume and bud burst. Comparison of family and animal models showed that unbalanced mating system led to only slight departures from maternal family assumptions in the progeny trial. However, neglecting the significant maternal effects led to an overestimation of the heritability. Overall, we highlighted the usefulness of relatedness pattern analyses using polymorphic molecular markers to accurately analyse tree sibling designs.     

Estimating genetic correlations in natural populations

Information on the genetic correlation between traits provides fundamental insight into the constraints on the evolutionary process. Estimates of such correlations are conventionally obtained by raising individuals of known relatedness in artificial environments. However, many species are not readily amenable to controlled breeding programmes, and considerable uncertainty exists over the extent to which estimates derived under benign laboratory conditions reflect the properties of populations in natural settings. Here, non-invasive methods that allow the estimation of genetic correlations from phenotypic measurements derived from individuals of unknown relatedness are introduced. Like the conventional approach, these methods demand large sample sizes in order to yield reasonably precise estimates, and special precautions need to be taken to eliminate bias from shared environmental effects. Provided the sample consists of at least 20% or so relatives, informative estimates of the genetic correlation are obtainable with sample sizes of several hundred individuals, particularly if supplemental information on relatedness is available from polymorphic molecular markers.     

Assumption-free estimation of heritability from genome-wide identity-by-descent sharing between full siblings

The study of continuously varying, quantitative traits is important in evolutionary biology, agriculture, and medicine. Variation in such traits is attributable to many, possibly interacting, genes whose expression may be sensitive to the environment, which makes their dissection into underlying causative factors difficult. An important population parameter for quantitative traits is heritability, the proportion of total variance that is due to genetic factors. Response to artificial and natural selection and the degree of resemblance between relatives are all a function of this parameter. Following the classic paper by R. A. Fisher in 1918, the estimation of additive and dominance genetic variance and heritability in populations is based upon the expected proportion of genes shared between different types of relatives, and explicit, often controversial and untestable models of genetic and non-genetic causes of family resemblance. With genome-wide coverage of genetic markers it is now possible to estimate such parameters solely within families using the actual degree of identity-by-descent sharing between relatives. Using genome scans on 4,401 quasi-independent sib pairs of which 3,375 pairs had phenotypes, we estimated the heritability of height from empirical genome-wide identity-by-descent sharing, which varied from 0.374 to 0.617 (mean 0.498, standard deviation 0.036). The variance in identity-by-descent sharing per chromosome and per genome was consistent with theory. The maximum likelihood estimate of the heritability for height was 0.80 with no evidence for non-genetic causes of sib resemblance, consistent with results from independent twin and family studies but using an entirely separate source of information. Our application shows that it is feasible to estimate genetic variance solely from within-family segregation and provides an independent validation of previously untestable assumptions. Given sufficient data, our new paradigm will allow the estimation of genetic variation for disease susceptibility and quantitative traits that is free from confounding with non-genetic factors and will allow partitioning of genetic variation into additive and non-additive components.     

Estimating quantitative genetic parameters in wild populations: a comparison of pedigree and genomic approaches

The estimation of quantitative genetic parameters in wild populations is generally limited by the accuracy and completeness of the available pedigree information. Using relatedness at genomewide markers can potentially remove this limitation and lead to less biased and more precise estimates. We estimated heritability, maternal genetic effects and genetic correlations for body size traits in an unmanaged long‐term study population of Soay sheep on St Kilda using three increasingly complete and accurate estimates of relatedness: (i) Pedigree 1, using observation‐derived maternal links and microsatellite‐derived paternal links; (ii) Pedigree 2, using SNP‐derived assignment of both maternity and paternity; and (iii) whole‐genome relatedness at 37 037 autosomal SNPs. In initial analyses, heritability estimates were strikingly similar for all three methods, while standard errors were systematically lower in analyses based on Pedigree 2 and genomic relatedness. Genetic correlations were generally strong, differed little between the three estimates of relatedness and the standard errors declined only very slightly with improved relatedness information. When partitioning maternal effects into separate genetic and environmental components, maternal genetic effects found in juvenile traits increased substantially across the three relatedness estimates. Heritability declined compared to parallel models where only a maternal environment effect was fitted, suggesting that maternal genetic effects are confounded with direct genetic effects and that more accurate estimates of relatedness were better able to separate maternal genetic effects from direct genetic effects. We found that the heritability captured by SNP markers asymptoted at about half the SNPs available, suggesting that denser marker panels are not necessarily required for precise and unbiased heritability estimates. Finally, we present guidelines for the use of genomic relatedness in future quantitative genetics studies in natural populations.     

Quantifying the genetic component of phenotypic variation in unpedigreed wild plants: tailoring genomic scan for within-population use

The study of adaptive genetic variation in natural populations is central to evolutionary biology. Quantitative genetics methods, however, are hardly applicable to long-lived organisms, and current knowledge on adaptive genetic variation in wild plants mostly refers to annuals and short-lived perennials. Studies on long-lived species are essential to explore possible life-history correlates of genetic variation, selection, and trait heritability. In this paper, we propose a method based on molecular markers to quantify the genetic basis of individual phenotypic differences in wild plants under natural conditions. Rather than focusing on inferring individual relatedness to estimate the heritability of phenotypic traits, we directly estimate the proportion of observed phenotypic variance that is statistically accounted for by genotypic differences between individuals. This is achieved by (i) identifying loci that are correlated across individuals with the phenotypic trait of interest by means of an amplified fragment length polymorphism (AFLP)-based explorative genomic scan, and (ii) fitting multiple regression and linear random effect models to estimate the effects of genotype, environment and genotype × environment on phenotypes. We apply this method to estimate genotypic and environmental effects on cumulative maternal fecundity in a wild population of the long-lived Viola cazorlensis monitored for 20 years. Results show that between 56–63% (depending on estimation method) of phenotypic variance in fecundity is accounted for by genotypic differences in 11 AFLP loci that are significantly related to fecundity. Genotype × environment effects accounted for 38% of fecundity variance, which may help to explain the unexpectedly high levels of genetic variance for fecundity found.     

How does selfing affect the genetic variance of quantitative traits? An updated meta‐analysis on empirical results in angiosperm species

Most theoretical works predict that selfing should reduce the level of additive genetic variance available for quantitative traits within natural populations. Despite a growing number of quantitative genetic studies undertaken during the last two decades, this prediction is still not well supported empirically. To resolve this issue and confirm or reject theoretical predictions, we reviewed quantitative trait heritability estimates from natural plant populations with different rates of self‐fertilization and carried out a meta‐analysis. In accordance with models of polygenic traits under stabilizing selection, we found that the fraction of additive genetic variance is negatively correlated with the selfing rate. Although the mating system explains a moderate fraction of the variance, the mean reduction of narrow‐sense heritability values between strictly allogamous and predominantly selfing populations is strong, around 60%. Because some nonadditive components of genetic variance become selectable under inbreeding, we determine whether self‐fertilization affects the relative contribution of these components to genetic variance by comparing narrow‐sense heritability estimates from outcrossing populations with broad‐sense heritability estimated in autogamous populations. Results suggest that these nonadditive components of variance may restore some genetic variance in predominantly selfing populations; it remains, however, uncertain how these nonadditive components will contribute to adaptation.     

Marker-based quantitative genetics in the wild?: the heritability and genetic correlation of chemical defenses in eucalyptus

Marker-based methods for estimating heritability and genetic correlation in the wild have attracted interest because traditional methods may be impractical or introduce bias via G x E effects, mating system variation, and sampling effects. However, they have not been widely used, especially in plants. A regression-based approach, which uses a continuous measure of genetic relatedness, promises to be particularly appropriate for use in plants with mixed-mating systems and overlapping generations. Using this method, we found significant narrow-sense heritability of foliar defense chemicals in a natural population of Eucalyptus melliodora. We also demonstrated a genetic basis for the phenotypic correlation underlying an ecological example of conditioned flavor aversion involving different biosynthetic pathways. Our results revealed that heritability estimates depend on the spatial scale of the analysis in a way that offers insight into the distribution of genetic and environmental variance. This study is the first to successfully use a marker-based method to measure quantitative genetic parameters in a tree. We suggest that this method will prove to be a useful tool in other studies and offer some recommendations for future applications of the method.     

Accuracy of dominant markers for estimation of relatedness and heritability in an experimental stand of Prosopis alba (Leguminosae)

The estimation of genetic components of phenotypic variance is based on the resemblance between relatives. In natural populations of most forest tree species without genealogical information, a possible alternative approach is the use of relatedness estimates obtained indirectly from molecular marker data. Heritability (h ²) is then estimated from the covariance of estimated relatedness and phenotypic resemblance. In a stand of Prosopis alba planted in 1991 in Argentina, relatedness was estimated for all individual pairs of trees by means of the information proceeding from 128 dominant markers (57 AFLPs and 71 ISSRs) and compared with estimates obtained from six microsatellite loci previously studied. We empirically compared the accuracy of different relatedness estimators based on dominant markers proposed by Lynch and Milligan (Mol Ecol 3:91–99, 1994), Hardy (Mol Ecol 12:1577–1588, 2003), Wang (Mol Ecol 13:3169–3178, 2004), and Ritland (Mol Ecol 14:3157–3165, 2005). Heritabilities of 13 quantitative traits were then estimated from the regression of pairwise phenotypic distances on pairwise relatedness according to Ritland (Genet Res 67:175–185, 1996a). Relatedness inferred from molecular markers was in all cases significantly correlated with actual relatedness, although Ritland's estimator showed the highest bias but the lowest variance. Dominant marker-based h ² estimates were evidently downwards biased and showed poor correlation with those based on family records. In conclusion, the use of dominant molecular markers evidently produces much greater underestimates of h ² than from using co-dominant ones, attributable to the lower accuracy in the indirect estimation of relatedness coefficient. Many traits with enough genetic variability as to respond readily to selection would remain undetected; only those with very high heritability would show significant h ² estimates.     

Intra-population comparison of vegetative and floral trait heritabilities estimated from molecular markers in wild Aquilegia populations

Measuring heritable genetic variation is important for understanding patterns of trait evolution in wild populations, and yet studies of quantitative genetic parameters estimated directly in the field are limited by logistic constraints, such as the difficulties of inferring relatedness among individuals in the wild. Marker-based approaches have received attention because they can potentially be applied directly to wild populations. For long-lived, self-compatible plant species where pedigrees are inadequate, the regression-based method proposed by Ritland has the appeal of estimating heritabilities from marker-based estimates of relatedness. The method has been difficult to implement in some plant populations, however, because it requires significant variance in relatedness across the population. Here, we show that the method can be readily applied to compare the ability of different traits to respond to selection, within populations. For several taxa of the perennial herb genus Aquilegia, we estimated heritabilities of floral and vegetative traits and, combined with estimates of natural selection, compared the ability to respond to selection of both types of traits under current conditions. The intra-population comparisons showed that vegetative traits have a higher potential for evolution, because although they are as heritable as floral traits, selection on them is stronger. These patterns of potential evolution are consistent with macroevolutionary trends in the European lineage of the genus.     

RAD‐sequencing for estimating genomic relatedness matrix‐based heritability in the wild: A case study in roe deer

Estimating the evolutionary potential of quantitative traits and reliably predicting responses to selection in wild populations are important challenges in evolutionary biology. The genomic revolution has opened up opportunities for measuring relatedness among individuals with precision, enabling pedigree‐free estimation of trait heritabilities in wild populations. However, until now, most quantitative genetic studies based on a genomic relatedness matrix (GRM) have focused on long‐term monitored populations for which traditional pedigrees were also available, and have often had access to knowledge of genome sequence and variability. Here, we investigated the potential of RAD‐sequencing for estimating heritability in a free‐ranging roe deer (Capreolous capreolus) population for which no prior genomic resources were available. We propose a step‐by‐step analytical framework to optimize the quality and quantity of the genomic data and explore the impact of the single nucleotide polymorphism (SNP) calling and filtering processes on the GRM structure and GRM‐based heritability estimates. As expected, our results show that sequence coverage strongly affects the number of recovered loci, the genotyping error rate and the amount of missing data. Ultimately, this had little effect on heritability estimates and their standard errors, provided that the GRM was built from a minimum number of loci (above 7,000). Genomic relatedness matrix‐based heritability estimates thus appear robust to a moderate level of genotyping errors in the SNP data set. We also showed that quality filters, such as the removal of low‐frequency variants, affect the relatedness structure of the GRM, generating lower h² estimates. Our work illustrates the huge potential of RAD‐sequencing for estimating GRM‐based heritability in virtually any natural population.     

Marker-inferred relatedness as a tool for detecting heritability in nature

This paper presents a perspective of how inferred relatedness, based on genetic marker data such as microsatellites or amplified fragment length polymorphisms (AFLPs), can be used to demonstrate quantitative genetic variation in natural populations. Variation at two levels is considered: among pairs of individuals within populations, and among pairs of subpopulations within a population. In the former, inferred pairwise relatedness, combined with trait measures, allow estimates of heritability 'in the wild'. In the latter, estimates of QST are obtained, in the absence of known heritabilities, via estimates of pairwise FST. Estimators of relatedness based on the 'Kronecker operator' are given. Both methods require actual variation of relationship, a rarely studied aspect of population structure, and not necessarily present. Some conditions for appropriate population structures in the wild are identified, in part through a review of recent studies.     

Pedigree-free animal models: the relatedness matrix reloaded

Animal models typically require a known genetic pedigree to estimate quantitative genetic parameters. Here we test whether animal models can alternatively be based on estimates of relatedness derived entirely from molecular marker data. Our case study is the morphology of a wild bird population, for which we report estimates of the genetic variance-covariance matrices (G) of six morphological traits using three methods: the traditional animal model; a molecular marker-based approach to estimate heritability based on Ritland's pairwise regression method; and a new approach using a molecular genealogy arranged in a relatedness matrix (R) to replace the pedigree in an animal model. Using the traditional animal model, we found significant genetic variance for all six traits and positive genetic covariance among traits. The pairwise regression method did not return reliable estimates of quantitative genetic parameters in this population, with estimates of genetic variance and covariance typically being very small or negative. In contrast, we found mixed evidence for the use of the pedigree-free animal model. Similar to the pairwise regression method, the pedigree-free approach performed poorly when the full-rank R matrix based on the molecular genealogy was employed. However, performance improved substantially when we reduced the dimensionality of the R matrix in order to maximize the signal to noise ratio. Using reduced-rank R matrices generated estimates of genetic variance that were much closer to those from the traditional model. Nevertheless, this method was less reliable at estimating covariances, which were often estimated to be negative. Taken together, these results suggest that pedigree-free animal models can recover quantitative genetic information, although the signal remains relatively weak. It remains to be determined whether this problem can be overcome by the use of a more powerful battery of molecular markers and improved methods for reconstructing genealogies.     

Heritability of Wing Morphology in a Natural Population of Drosophila gouveai

The natural and laboratory heritabilities of a series of parameters related to wing size and shape were estimated in a population of Drosophila gouveai (repleta group) under field and laboratory conditions. A morphometric analysis was done using 17 wing parameters related to wing landmark positions obtained using the method of the best adjustment of an ellipse to the wing edge. Three parameters (thetaA, thetaC and thetaD) showed highly significant heritability in the wild (average 0.61), whereas only wing size (W(SI)) had significant heritability in the laboratory (0.71). The additive genetic variance of most parameters was greater in the wild than in the laboratory. These results showed that some parameters possessed a substantial genetic additive component in their phenotypic variance, and that morphometric parameters of D. gouveai wings are appropriate quantitative markers for assessing morphological differentiation among populations.     

HOW MUCH HERITABLE VARIATION CAN BE MAINTAINED IN FINITE POPULATIONS BY MUTATION–SELECTION BALANCE?

The joint effects of stabilizing selection, mutation, recombination, and random drift on the genetic variability of a polygenic character in a finite population are investigated. A simulation study is performed to test the validity of various analytical predictions on the equilibrium genetic variance. A new formula for the expected equilibrium variance is derived that approximates the observed equilibrium variance very closely for all parameter combinations we have tested. The computer model simulates the continuum-of-alleles model of Crow and Kimura. However, it is completely stochastic in the sense that it models evolution as a Markov process and does not use any deterministic evolution equations. The theoretical results are compared with heritability estimates from laboratory and natural populations. Heritabilities ranging from 20% to 50%, as observed even in lab populations under a constant environment, can only be explained by a mutation-selection balance if the phenotypic character is neutral or the number of genes contributing to the trait is sufficiently high, typically several hundred, or if there are a few highly variable loci that influence quantitative traits.