Showing papers in "Heredity in 2014"

Genomic prediction in CIMMYT maize and wheat breeding programs

Abstract: Genomic selection (GS) has been implemented in animal and plant species, and is regarded as a useful tool for accelerating genetic gains. Varying levels of genomic prediction accuracy have been obtained in plants, depending on the prediction problem assessed and on several other factors, such as trait heritability, the relationship between the individuals to be predicted and those used to train the models for prediction, number of markers, sample size and genotype × environment interaction (GE). The main objective of this article is to describe the results of genomic prediction in International Maize and Wheat Improvement Center's (CIMMYT's) maize and wheat breeding programs, from the initial assessment of the predictive ability of different models using pedigree and marker information to the present, when methods for implementing GS in practical global maize and wheat breeding programs are being studied and investigated. Results show that pedigree (population structure) accounts for a sizeable proportion of the prediction accuracy when a global population is the prediction problem to be assessed. However, when the prediction uses unrelated populations to train the prediction equations, prediction accuracy becomes negligible. When genomic prediction includes modeling GE, an increase in prediction accuracy can be achieved by borrowing information from correlated environments. Several questions on how to incorporate GS into CIMMYT's maize and wheat programs remain unanswered and subject to further investigation, for example, prediction within and between related bi-parental crosses. Further research on the quantification of breeding value components for GS in plant breeding populations is required.

Supergenes and their role in evolution.

Abstract: Adaptation is commonly a multidimensional problem, with changes in multiple traits required to match a complex environment. This is epitomized by balanced polymorphisms in which multiple phenotypes co-exist and are maintained in a population by a balance of selective forces. Consideration of such polymorphisms led to the concept of the supergene, where alternative phenotypes in a balanced polymorphism segregate as if controlled by a single genetic locus, resulting from tight genetic linkage between multiple functional loci. Recently, the molecular basis for several supergenes has been resolved. Thus, major chromosomal inversions have been shown to be associated with polymorphisms in butterflies, ants and birds, offering a mechanism for localised reduction in recombination. In several examples of plant self-incompatibility, the functional role of multiple elements within the supergene architecture has been demonstrated, conclusively showing that balanced polymorphism can be maintained at multiple coadapted and tightly linked elements. Despite recent criticism, we argue that the supergene concept remains relevant and is more testable than ever with modern molecular methods.

Toward genomic prediction from whole-genome sequence data: impact of sequencing design on genotype imputation and accuracy of predictions.

Abstract: Genomic prediction from whole-genome sequence data is attractive, as the accuracy of genomic prediction is no longer bounded by extent of linkage disequilibrium between DNA markers and causal mutations affecting the trait, given the causal mutations are in the data set. A cost-effective strategy could be to sequence a small proportion of the population, and impute sequence data to the rest of the reference population. Here, we describe strategies for selecting individuals for sequencing, based on either pedigree relationships or haplotype diversity. Performance of these strategies (number of variants detected and accuracy of imputation) were evaluated in sequence data simulated through a real Belgian Blue cattle pedigree. A strategy (AHAP), which selected a subset of individuals for sequencing that maximized the number of unique haplotypes (from single-nucleotide polymorphism panel data) sequenced gave good performance across a range of variant minor allele frequencies. We then investigated the optimum number of individuals to sequence by fold coverage given a maximum total sequencing effort. At 600 total fold coverage (x 600), the optimum strategy was to sequence 75 individuals at eightfold coverage. Finally, we investigated the accuracy of genomic predictions that could be achieved. The advantage of using imputed sequence data compared with dense SNP array genotypes was highly dependent on the allele frequency spectrum of the causative mutations affecting the trait. When this followed a neutral distribution, the advantage of the imputed sequence data was small; however, when the causal mutations all had low minor allele frequencies, using the sequence data improved the accuracy of genomic prediction by up to 30%.

Entering the second century of maize quantitative genetics.

Abstract: Maize is the most widely grown cereal in the world. In addition to its role in global agriculture, it has also long served as a model organism for genetic research. Maize stands at a genetic crossroads, as it has access to all the tools available for plant genetics but exhibits a genetic architecture more similar to other outcrossing organisms than to self-pollinating crops and model plants. In this review, we summarize recent advances in maize genetics, including the development of powerful populations for genetic mapping and genome-wide association studies (GWAS), and the insights these studies yield on the mechanisms underlying complex maize traits. Most maize traits are controlled by a large number of genes, and linkage analysis of several traits implicates a ‘common gene, rare allele' model of genetic variation where some genes have many individually rare alleles contributing. Most natural alleles exhibit small effect sizes with little-to-no detectable pleiotropy or epistasis. Additionally, many of these genes are locked away in low-recombination regions that encourage the formation of multi-gene blocks that may underlie maize's strong heterotic effect. Domestication left strong marks on the maize genome, and some of the differences in trait architectures may be due to different selective pressures over time. Overall, maize's advantages as a model system make it highly desirable for studying the genetics of outcrossing species, and results from it can provide insight into other such species, including humans.

The quantitative genetics of indirect genetic effects: a selective review of modelling issues.

Abstract: Indirect genetic effects (IGE) occur when the genotype of an individual affects the phenotypic trait value of another conspecific individual. IGEs can have profound effects on both the magnitude and the direction of response to selection. Models of inheritance and response to selection in traits subject to IGEs have been developed within two frameworks; a trait-based framework in which IGEs are specified as a direct consequence of individual trait values, and a variance-component framework in which phenotypic variance is decomposed into a direct and an indirect additive genetic component. This work is a selective review of the quantitative genetics of traits affected by IGEs, with a focus on modelling, estimation and interpretation issues. It includes a discussion on variance-component vs trait-based models of IGEs, a review of issues related to the estimation of IGEs from field data, including the estimation of the interaction coefficient Ψ (psi), and a discussion on the relevance of IGEs for response to selection in cases where the strength of interaction varies among pairs of individuals. An investigation of the trait-based model shows that the interaction coefficient Ψ may deviate considerably from the corresponding regression coefficient when feedback occurs. The increasing research effort devoted to IGEs suggests that they are a widespread phenomenon, probably particularly in natural populations and plants. Further work in this field should considerably broaden our understanding of the quantitative genetics of inheritance and response to selection in relation to the social organisation of populations.

Bridging the gap between marker-assisted and genomic selection of heading time and plant height in hybrid wheat

Abstract: Based on data from field trials with a large collection of 135 elite winter wheat inbred lines and 1604 F1 hybrids derived from them, we compared the accuracy of prediction of marker-assisted selection and current genomic selection approaches for the model traits heading time and plant height in a cross-validation approach. For heading time, the high accuracy seen with marker-assisted selection severely dropped with genomic selection approaches RR-BLUP (ridge regression best linear unbiased prediction) and BayesCπ, whereas for plant height, accuracy was low with marker-assisted selection as well as RR-BLUP and BayesCπ. Differences in the linkage disequilibrium structure of the functional and single-nucleotide polymorphism markers relevant for the two traits were identified in a simulation study as a likely explanation for the different trends in accuracies of prediction. A new genomic selection approach, weighted best linear unbiased prediction (W-BLUP), designed to treat the effects of known functional markers more appropriately, proved to increase the accuracy of prediction for both traits and thus closes the gap between marker-assisted and genomic selection.

Male killing Spiroplasma protects Drosophila melanogaster against two parasitoid wasps

Abstract: Maternally transmitted associations between endosymbiotic bacteria and insects are diverse and widespread in nature. Owing to imperfect vertical transmission, many heritable microbes have evolved compensational mechanisms to enhance their persistence in host lineages, such as manipulating host reproduction and conferring fitness benefits to host. Symbiont-mediated defense against natural enemies of hosts is increasingly recognized as an important mechanism by which endosymbionts enhance host fitness. Members of the genus Spiroplasma associated with distantly related Drosophila hosts are known to engage in either reproductive parasitism (i.e., male killing) or defense against natural enemies (the parasitic wasp Leptopilina heterotoma and a nematode). A male-killing strain of Spiroplasma (strain Melanogaster Sex Ratio Organism (MSRO)) co-occurs with Wolbachia (strain wMel) in certain wild populations of the model organism Drosophila melanogaster. We examined the effects of Spiroplasma MSRO and Wolbachia wMel on Drosophila survival against parasitism by two common wasps, Leptopilina heterotoma and Leptopilina boulardi, that differ in their host ranges and host evasion strategies. The results indicate that Spiroplasma MSRO prevents successful development of both wasps, and confers a small, albeit significant, increase in larva-to-adult survival of flies subjected to wasp attacks. We modeled the conditions under which defense can contribute to Spiroplasma persistence. Wolbachia also confers a weak, but significant, survival advantage to flies attacked by L. heterotoma. The host protective effects exhibited by Spiroplasma and Wolbachia are additive and may provide the conditions for such cotransmitted symbionts to become mutualists. Occurrence of Spiroplasma-mediated protection against distinct parasitoids in divergent Drosophila hosts suggests a general protection mechanism.

Accuracy of genomic selection models in a large population of open-pollinated families in white spruce.

Abstract: Genomic selection (GS) is of interest in breeding because of its potential for predicting the genetic value of individuals and increasing genetic gains per unit of time. To date, very few studies have reported empirical results of GS potential in the context of large population sizes and long breeding cycles such as for boreal trees. In this study, we assessed the effectiveness of marker-aided selection in an undomesticated white spruce (Picea glauca (Moench) Voss) population of large effective size using a GS approach. A discovery population of 1694 trees representative of 214 open-pollinated families from 43 natural populations was phenotyped for 12 wood and growth traits and genotyped for 6385 single-nucleotide polymorphisms (SNPs) mined in 2660 gene sequences. GS models were built to predict estimated breeding values using all the available SNPs or SNP subsets of the largest absolute effects, and they were validated using various cross-validation schemes. The accuracy of genomic estimated breeding values (GEBVs) varied from 0.327 to 0.435 when the training and the validation data sets shared half-sibs that were on average 90% of the accuracies achieved through traditionally estimated breeding values. The trend was also the same for validation across sites. As expected, the accuracy of GEBVs obtained after cross-validation with individuals of unknown relatedness was lower with about half of the accuracy achieved when half-sibs were present. We showed that with the marker densities used in the current study, predictions with low to moderate accuracy could be obtained within a large undomesticated population of related individuals, potentially resulting in larger gains per unit of time with GS than with the traditional approach.

The evolution of genomic imprinting: theories, predictions and empirical tests

Abstract: The epigenetic phenomenon of genomic imprinting has motivated the development of numerous theories for its evolutionary origins and genomic distribution In this review, we examine the three theories that have best withstood theoretical and empirical scrutiny These are: Haig and colleagues' kinship theory; Day and Bonduriansky's sexual antagonism theory; and Wolf and Hager's maternal-offspring coadaptation theory These theories have fundamentally different perspectives on the adaptive significance of imprinting The kinship theory views imprinting as a mechanism to change gene dosage, with imprinting evolving because of the differential effect that gene dosage has on the fitness of matrilineal and patrilineal relatives The sexual antagonism and maternal-offspring coadaptation theories view genomic imprinting as a mechanism to modify the resemblance of an individual to its two parents, with imprinting evolving to increase the probability of expressing the fitter of the two alleles at a locus In an effort to stimulate further empirical work on the topic, we carefully detail the logic and assumptions of all three theories, clarify the specific predictions of each and suggest tests to discriminate between these alternative theories for why particular genes are imprinted

Comparing G: multivariate analysis of genetic variation in multiple populations.

Abstract: The additive genetic variance–covariance matrix (G) summarizes the multivariate genetic relationships among a set of traits. The geometry of G describes the distribution of multivariate genetic variance, and generates genetic constraints that bias the direction of evolution. Determining if and how the multivariate genetic variance evolves has been limited by a number of analytical challenges in comparing G-matrices. Current methods for the comparison of G typically share several drawbacks: metrics that lack a direct relationship to evolutionary theory, the inability to be applied in conjunction with complex experimental designs, difficulties with determining statistical confidence in inferred differences and an inherently pair-wise focus. Here, we present a cohesive and general analytical framework for the comparative analysis of G that addresses these issues, and that incorporates and extends current methods with a strong geometrical basis. We describe the application of random skewers, common subspace analysis, the 4th-order genetic covariance tensor and the decomposition of the multivariate breeders equation, all within a Bayesian framework. We illustrate these methods using data from an artificial selection experiment on eight traits in Drosophila serrata, where a multi-generational pedigree was available to estimate G in each of six populations. One method, the tensor, elegantly captures all of the variation in genetic variance among populations, and allows the identification of the trait combinations that differ most in genetic variance. The tensor approach is likely to be the most generally applicable method to the comparison of G-matrices from any sampling or experimental design.

Genome-wide association study of handedness excludes simple genetic models

Abstract: Handedness is a human behavioural phenotype that appears to be congenital, and is often assumed to be inherited, but for which the developmental origin and underlying causation(s) have been elusive. Models of the genetic basis of variation in handedness have been proposed that fit different features of the observed resemblance between relatives, but none has been decisively tested or a corresponding causative locus identified. In this study, we applied data from well-characterised individuals studied at the London Twin Research Unit. Analysis of genome-wide SNP data from 3940 twins failed to identify any locus associated with handedness at a genome-wide level of significance. The most straightforward interpretation of our analyses is that they exclude the simplest formulations of the ‘right-shift' model of Annett and the ‘dextral/chance' model of McManus, although more complex modifications of those models are still compatible with our observations. For polygenic effects, our study is inadequately powered to reliably detect alleles with effect sizes corresponding to an odds ratio of 1.2, but should have good power to detect effects at an odds ratio of 2 or more.

Genetic consequences of habitat fragmentation during a range expansion

Abstract: We investigate the effect of habitat fragmentation on the genetic diversity of a species experiencing a range expansion. These two evolutionary processes have not been studied yet, at the same time, owing to the difficulties of deriving analytic results for non-equilibrium models. Here we provide a description of their interaction by using extensive spatial and temporal coalescent simulations and we suggest guidelines for a proper genetic sampling to detect fragmentation. To model habitat fragmentation, we simulated a two-dimensional lattice of demes partitioned into groups (patches) by adding barriers to dispersal. After letting a population expand on this grid, we sampled lineages from the lattice at several scales and studied their coalescent history. We find that in order to detect fragmentation, one needs to extensively sample at a local level rather than at a landscape level. This is because the gene genealogy of a scattered sample is less sensitive to the presence of genetic barriers. Considering the effect of temporal changes of fragmentation intensities, we find that at least 10, but often >100, generations are needed to affect local genetic diversity and population structure. This result explains why recent habitat fragmentation does not always lead to detectable signatures in the genetic structure of populations. Finally, as expected, long-distance dispersal increases local genetic diversity and decreases levels of population differentiation, efficiently counteracting the effects of fragmentation.

The evolution of quantitative traits in complex environments.

Abstract: Species inhabit complex environments and respond to selection imposed by numerous abiotic and biotic conditions that vary in both space and time. Environmental heterogeneity strongly influences trait evolution and patterns of adaptive population differentiation. For example, heterogeneity can favor local adaptation, or can promote the evolution of plastic genotypes that alter their phenotypes based on the conditions they encounter. Different abiotic and biotic agents of selection can act synergistically to either accelerate or constrain trait evolution. The environmental context has profound effects on quantitative genetic parameters. For instance, heritabilities measured in controlled conditions often exceed those measured in the field; thus, laboratory experiments could overestimate the potential for a population to respond to selection. Nevertheless, most studies of the genetic basis of ecologically relevant traits are conducted in simplified laboratory environments, which do not reflect the complexity of nature. Here, we advocate for manipulative field experiments in the native ranges of plant species that differ in mating system, life-history strategy and growth form. Field studies are vital to evaluate the roles of disparate agents of selection, to elucidate the targets of selection and to develop a nuanced perspective on the evolution of quantitative traits. Quantitative genetics field studies will also shed light on the potential for natural populations to adapt to novel climates in highly fragmented landscapes. Drawing from our experience with the ecological model system Boechera (Brassicaceae), we discuss advancements possible through dedicated field studies, highlight future research directions and examine the challenges associated with field studies.

Genome-wide association study identifies novel loci associated with resistance to bovine tuberculosis.

Abstract: Tuberculosis (TB) caused by Mycobacterium bovis is a re-emerging disease of livestock that is of major economic importance worldwide, as well as being a zoonotic risk. There is significant heritability for host resistance to bovine TB (bTB) in dairy cattle. To identify resistance loci for bTB, we undertook a genome-wide association study in female Holstein–Friesian cattle with 592 cases and 559 age-matched controls from case herds. Cases and controls were categorised into distinct phenotypes: skin test and lesion positive vs skin test negative on multiple occasions, respectively. These animals were genotyped with the Illumina BovineHD 700K BeadChip. Genome-wide rapid association using linear and logistic mixed models and regression (GRAMMAR), regional heritability mapping (RHM) and haplotype-sharing analysis identified two novel resistance loci that attained chromosome-wise significance, protein tyrosine phosphatase receptor T (PTPRT; P=4.8 × 10−7) and myosin IIIB (MYO3B; P=5.4 × 10−6). We estimated that 21% of the phenotypic variance in TB resistance could be explained by all of the informative single-nucleotide polymorphisms, of which the region encompassing the PTPRT gene accounted for 6.2% of the variance and a further 3.6% was associated with a putative copy number variant in MYO3B. The results from this study add to our understanding of variation in host control of infection and suggest that genetic marker-based selection for resistance to bTB has the potential to make a significant contribution to bTB control.

Genome-wide signatures of population bottlenecks and diversifying selection in European wolves

Abstract: Genomic resources developed for domesticated species provide powerful tools for studying the evolutionary history of their wild relatives. Here we use 61K single-nucleotide polymorphisms (SNPs) evenly spaced throughout the canine nuclear genome to analyse evolutionary relationships among the three largest European populations of grey wolves in comparison with other populations worldwide, and investigate genome-wide effects of demographic bottlenecks and signatures of selection. European wolves have a discontinuous range, with large and connected populations in Eastern Europe and relatively smaller, isolated populations in Italy and the Iberian Peninsula. Our results suggest a continuous decline in wolf numbers in Europe since the Late Pleistocene, and long-term isolation and bottlenecks in the Italian and Iberian populations following their divergence from the Eastern European population. The Italian and Iberian populations have low genetic variability and high linkage disequilibrium, but relatively few autozygous segments across the genome. This last characteristic clearly distinguishes them from populations that underwent recent drastic demographic declines or founder events, and implies long-term bottlenecks in these two populations. Although genetic drift due to spatial isolation and bottlenecks seems to be a major evolutionary force diversifying the European populations, we detected 35 loci that are putatively under diversifying selection. Two of these loci flank the canine platelet-derived growth factor gene, which affects bone growth and may influence differences in body size between wolf populations. This study demonstrates the power of population genomics for identifying genetic signals of demographic bottlenecks and detecting signatures of directional selection in bottlenecked populations, despite their low background variability.

Estimating genome-wide heterozygosity: effects of demographic history and marker type

Abstract: Heterozygosity–fitness correlations (HFCs) are often used to link individual genetic variation to differences in fitness. However, most studies examining HFCs find weak or no correlations. Here, we derive broad theoretical predictions about how many loci are needed to adequately measure genomic heterozygosity assuming different levels of identity disequilibrium (ID), a proxy for inbreeding. We then evaluate the expected ability to detect HFCs using an empirical data set of 200 microsatellites and 412 single nucleotide polymorphisms (SNPs) genotyped in two populations of bighorn sheep (Ovis canadensis), with different demographic histories. In both populations, heterozygosity was significantly correlated across marker types, although the strength of the correlation was weaker in a native population compared with one founded via translocation and later supplemented with additional individuals. Despite being bi-allelic, SNPs had similar correlations to genome-wide heterozygosity as microsatellites in both populations. For both marker types, this association became stronger and less variable as more markers were considered. Both populations had significant levels of ID; however, estimates were an order of magnitude lower in the native population. As with heterozygosity, SNPs performed similarly to microsatellites, and precision and accuracy of the estimates of ID increased as more loci were considered. Although dependent on the demographic history of the population considered, these results illustrate that genome-wide heterozygosity, and therefore HFCs, are best measured by a large number of markers, a feat now more realistically accomplished with SNPs than microsatellites.

Using next-generation RNA sequencing to identify imprinted genes

Abstract: Genomic imprinting is manifested as differential allelic expression (DAE) depending on the parent-of-origin. The most direct way to identify imprinted genes is to directly score the DAE in a context where one can identify which parent transmitted each allele. Because many genes display DAE, simply scoring DAE in an individual is not sufficient to identify imprinted genes. In this paper, we outline many technical aspects of a scheme for identification of imprinted genes that makes use of RNA sequencing (RNA-seq) from tissues isolated from F1 offspring derived from the pair of reciprocal crosses. Ideally, the parental lines are from two inbred strains that are not closely related to each other. Aspects of tissue purity, RNA extraction, library preparation and bioinformatic inference of imprinting are all covered. These methods have already been applied in a number of organisms, and one of the most striking results is the evolutionary fluidity with which novel imprinted genes are gained and lost within genomes. The general methodology is also applicable to a wide range of other biological problems that require quantification of allele-specific expression using RNA-seq, such as cis-regulation of gene expression, X chromosome inactivation and random monoallelic expression.

Evidence for multiple sex-determining loci in Tasmanian Atlantic salmon (Salmo salar).

Abstract: Phenotypic sex in salmonids is determined primarily by a genetic male heterogametic system; yet, sex reversal can be accomplished via hormonal treatment. In Tasmanian Atlantic salmon aquaculture, to overcome problems associated with early sexual maturation in males, sex-reversed females are crossed with normal females to produce all female stock. However, phenotypic distinction of sex-reversed females (neo-males) from true males is problematic. We set out to identify genetic markers that could make this distinction. Microsatellite markers from chromosome 2 (Ssa02), to which the sex-determining locus (SEX) has been mapped in two Scottish Atlantic salmon families, did not predict sex in a pilot study of seven families. A TaqMan 64 SNP genome-wide scan suggested SEX was on Ssa06 in these families, and this was confirmed by microsatellite markers. A survey of 58 families in total representing 38 male lineages in the SALTAS breeding program found that 34 of the families had SEX on Ssa02, in 22 of the families SEX was on Ssa06, and two of the families had a third SEX locus, on Ssa03. A PCR test using primers designed from the recently published sdY gene is consistent with Tasmanian Atlantic salmon having a single sex-determining gene that may be located on at least three linkage groups.

Coadaptation and conflict, misconception and muddle, in the evolution of genomic imprinting

Abstract: Coadaptation and conflict, misconception and muddle, in the evolution of genomic imprinting

The specification of imprints in mammals

Abstract: At the heart of genomic imprinting in mammals are imprinting control regions (ICRs), which are the discrete genetic elements that confer imprinted monoallelic expression to several genes in imprinted gene clusters. A characteristic of the known ICRs is that they acquire different epigenetic states, exemplified by differences in DNA methylation, in the sperm and egg, and these imprint marks remain on the sperm- and oocyte-derived alleles into the next generation as a lifelong memory of parental origin. Although there has been much focus on gametic marking of ICRs as the point of imprint specification, recent mechanistic studies and genome-wide DNA methylation profiling do not support the existence of a specific imprinting machinery in germ cells. Rather, ICRs are part of more widespread methylation events that occur during gametogenesis. Instead, a decisive component in the specification of imprints is the choice of which sites of gamete-derived methylation to maintain in the zygote and preimplantation embryo at a time when much of the remainder of the genome is being demethylated. Among the factors involved in this selection, the zinc-finger protein Zfp57 can be regarded as an imprint-specific, sequence-specific DNA binding factor responsible for maintaining methylation at most ICRs. The recent insights into the balance of gametic and zygotic contributions to imprint specification should help understand mechanistic opportunities and constraints on the evolution of imprinting in mammals.

Competition as a source of constraint on life history evolution in natural populations

Abstract: Competition among individuals is central to our understanding of ecology and population dynamics. However, it could also have major implications for the evolution of resource-dependent life history traits (for example, growth, fecundity) that are important determinants of fitness in natural populations. This is because when competition occurs, the phenotype of each individual will be causally influenced by the phenotypes, and so the genotypes, of competitors. Theory tells us that indirect genetic effects arising from competitive interactions will give rise to the phenomenon of ‘evolutionary environmental deterioration', and act as a source of evolutionary constraint on resource-dependent traits under natural selection. However, just how important this constraint is remains an unanswered question. This article seeks to stimulate empirical research in this area, first highlighting some patterns emerging from life history studies that are consistent with a competition-based model of evolutionary constraint, before describing several quantitative modelling strategies that could be usefully applied. A recurrent theme is that rigorous quantification of a competition's impact on life history evolution will require an understanding of the causal pathways and behavioural processes by which genetic (co)variance structures arise. Knowledge of the G-matrix among life history traits is not, in and of itself, sufficient to identify the constraints caused by competition.

Post-GWAS: where next? More samples, more SNPs or more biology?

Abstract: The power of genome-wide association studies (GWAS) rests on several foundations: (i) there is a significant amount of additive genetic variation, (ii) individual causal polymorphisms often have sizable effects and (iii) they segregate at moderate-to-intermediate frequencies, or will be effectively 'tagged' by polymorphisms that do. Each of these assumptions has recently been questioned. (i) Why should genetic variation appear additive given that the underlying molecular networks are highly nonlinear? (ii) A new generation of relatedness-based analyses directs us back to the nearly infinitesimal model for effect sizes that quantitative genetics was long based upon. (iii) Larger effect causal polymorphisms are often low frequency, as selection might lead us to expect. Here, we review these issues and other findings that appear to question many of the foundations of the optimism GWAS prompted. We then present a roadmap emerging as one possible future for quantitative genetics. We argue that in future GWAS should move beyond purely statistical grounds. One promising approach is to build upon the combination of population genetic models and molecular biological knowledge. This combined treatment, however, requires fitting experimental data to models that are very complex, as well as accurate capturing of the uncertainty of resulting inference. This problem can be resolved through Bayesian analysis and tools such as approximate Bayesian computation-a method growing in popularity in population genetic analysis. We show a case example of anterior-posterior segmentation in Drosophila, and argue that similar approaches will be helpful as a GWAS augmentation, in human and agricultural research.

Assessing patterns of hybridization between North Atlantic eels using diagnostic single-nucleotide polymorphisms

Abstract: The two North Atlantic eel species, the European eel (Anguilla anguilla) and the American eel (Anguilla rostrata), spawn in partial sympatry in the Sargasso Sea, providing ample opportunity to interbreed. In this study, we used a RAD (Restriction site Associated DNA) sequencing approach to identify species-specific diagnostic single-nucleotide polymorphisms (SNPs) and design a low-density array that combined with screening of a diagnostic mitochondrial DNA marker. Eels from Iceland (N=159) and from the neighboring Faroe Islands (N=29) were genotyped, along with 94 larvae (49 European and 45 American eel) collected in the Sargasso Sea. Our SNP survey showed that the majority of Icelandic eels are pure European eels but there is also an important contribution of individuals of admixed ancestry (10.7%). Although most of the hybrids were identified as F1 hybrids from European eel female × American eel male crosses, backcrosses were also detected, including a first-generation backcross (F1 hybrid × pure European eel) and three individuals identified as second-generation backcrosses originating from American eel × F1 hybrid backcrosses interbreeding with pure European eels. In comparison, no hybrids were observed in the Faroe Islands, the closest bodies of land to Iceland. It is possible that hybrids show an intermediate migratory behaviour between the two parental species that ultimately brings hybrid larvae to the shores of Iceland, situated roughly halfway between the Sargasso Sea and Europe. Only two hybrids were observed among Sargasso Sea larvae, both backcrosses, but no F1 hybrids, that points to temporal variation in the occurrence of hybridization.

Experimental evolution and the dynamics of genomic mutation rate modifiers.

Abstract: Because genes that affect mutation rates are themselves subject to mutation, mutation rates can be influenced by natural selection and other evolutionary forces. The population genetics of mutation rate modifier alleles has been a subject of theoretical interest for many decades. Here, we review experimental contributions to our understanding of mutation rate modifier dynamics. Numerous evolution experiments have shown that mutator alleles (modifiers that elevate the genomic mutation rate) can readily rise to high frequencies via genetic hitchhiking in non-recombining microbial populations. Whereas these results certainly provide an explanatory framework for observations of sporadically high mutation rates in pathogenic microbes and in cancer lineages, it is nonetheless true that most natural populations have very low mutation rates. This raises the interesting question of how mutator hitchhiking is suppressed or its phenotypic effect reversed in natural populations. Very little experimental work has addressed this question; with this in mind, we identify some promising areas for future experimental investigation.

Quantitative genetic study of the adaptive process

Abstract: The additive genetic variance with respect to absolute fitness, VA(W), divided by mean absolute fitness, , sets the rate of ongoing adaptation. Fisher's key insight yielding this quantitative prediction of adaptive evolution, known as the Fundamental Theorem of Natural Selection, is well appreciated by evolutionists. Nevertheless, extremely scant information about VA(W) is available for natural populations. Consequently, the capacity for fitness increase via natural selection is unknown. Particularly in the current context of rapid environmental change, which is likely to reduce fitness directly and, consequently, the size and persistence of populations, the urgency of advancing understanding of immediate adaptive capacity is extreme. We here explore reasons for the dearth of empirical information about VA(W), despite its theoretical renown and critical evolutionary role. Of these reasons, we suggest that expectations that VA(W) is negligible, in general, together with severe statistical challenges of estimating it, may largely account for the limited empirical emphasis on it. To develop insight into the dynamics of VA(W) in a changing environment, we have conducted individual-based genetically explicit simulations. We show that, as optimizing selection on a trait changes steadily over generations, VA(W) can grow considerably, supporting more rapid adaptation than would the VA(W) of the base population. We call for direct evaluation of VA(W) and in support of prediction of rates adaptive evolution, and we advocate for the use of aster modeling as a rigorous basis for achieving this goal.

Relatedness severely impacts accuracy of marker-assisted selection for disease resistance in hybrid wheat.

Abstract: The accuracy of genomic selection depends on the relatedness between the members of the set in which marker effects are estimated based on evaluation data and the types for which performance is predicted. Here, we investigate the impact of relatedness on the performance of marker-assisted selection for fungal disease resistance in hybrid wheat. A large and diverse mapping population of 1739 elite European winter wheat inbred lines and hybrids was evaluated for powdery mildew, leaf rust and stripe rust resistance in multi-location field trials and fingerprinted with 9 k and 90 k SNP arrays. Comparison of the accuracies of prediction achieved with data sets from the two marker arrays revealed a crucial role for a sufficiently high marker density in genome-wide association mapping. Cross-validation studies using test sets with varying degrees of relationship to the corresponding estimation sets revealed that close relatedness leads to a substantial increase in the proportion of total genotypic variance explained by the identified QTL and consequently to an overoptimistic judgment of the precision of marker-assisted selection.

Divergence with gene flow within the recent chipmunk radiation ( Tamias )

Abstract: Increasing data have supported the importance of divergence with gene flow (DGF) in the generation of biological diversity. In such cases, lineage divergence occurs on a shorter timescale than does the completion of reproductive isolation. Although it is critical to explore the mechanisms driving divergence and preventing homogenization by hybridization, it is equally important to document cases of DGF in nature. Here we synthesize data that have accumulated over the last dozen or so years on DGF in the chipmunk (Tamias) radiation with new data that quantify very high rates of mitochondrial DNA (mtDNA) introgression among para- and sympatric species in the T. quadrivittatus group in the central and southern Rocky Mountains. These new data (188 cytochrome b sequences) bring the total number of sequences up to 1871; roughly 16% (298) of the chipmunks we have sequenced exhibit introgressed mtDNA. This includes ongoing introgression between subspecies and between both closely related and distantly related taxa. In addition, we have identified several taxa that are apparently fixed for ancient introgressions and in which there is no evidence of ongoing introgression. A recurrent observation is that these introgressions occur between ecologically and morphologically diverged, sometimes non-sister taxa that engage in well-documented niche partitioning. Thus, the chipmunk radiation in western North America represents an excellent mammalian example of speciation in the face of recurrent gene flow among lineages and where biogeography, habitat differentiation and mating systems suggest important roles for both ecological and sexual selection.

Comparing the intersex genetic correlation for fitness across novel environments in the fruit fly, Drosophila serrata

Abstract: Sexually antagonistic genetic variation can pose limits to the independent evolution and adaptation of the sexes. The extent of sexually antagonistic variation is reflected in the intersex genetic correlation for fitness (rwFM). Previous estimates of this correlation have been mostly limited to populations in environments to which they are already well adapted, making it difficult to gauge the importance of sexually antagonistic genetic variance during the early stages of adaptation, such as that occurring following abrupt environmental change or upon the colonization of new habitat. Here we assayed male and female lifetime fitness in a population of Drosophila serrata in four novel laboratory environments. We found that rwFM varied significantly across environments, with point estimates ranging from positive to negative values of considerable magnitude. We also found that the variability among estimates was because, at least in part, of significant differences among environments in the genetic variances of both male and female fitness, with no evidence of any significant changes in the intersex covariance itself, although standard errors of these estimates were large. Our results illustrate the unpredictable nature of rwFM in novel environments and suggest that, although sexually antagonistic genetic variance can be pronounced in some novel environments, it may have little effect in constraining the early stages of adaptation in others.

Non-conflict theories for the evolution of genomic imprinting

Abstract: Theories focused on kinship and the genetic conflict it induces are widely considered to be the primary explanations for the evolution of genomic imprinting. However, there have appeared many competing ideas that do not involve kinship/conflict. These ideas are often overlooked because kinship/conflict is entrenched in the literature, especially outside evolutionary biology. Here we provide a critical overview of these non-conflict theories, providing an accessible perspective into this literature. We suggest that some of these alternative hypotheses may, in fact, provide tenable explanations of the evolution of imprinting for at least some loci.

Genetic diversity and classification of Oryza sativa with emphasis on Chinese rice germplasm.

Abstract: Despite extensive studies on cultivated rice, the genetic structure and subdivision of this crop remain unclear at both global and local scales. Using 84 nuclear simple sequence repeat markers, we genotyped a panel of 153 global rice cultivars covering all previously recognized groups and 826 cultivars representing the diversity of Chinese rice germplasm. On the basis of model-based grouping, neighbour-joining tree and principal coordinate analysis, we confirmed the widely accepted five major groups of rice cultivars (indica, aus, aromatic, temperate japonica and tropical japonica), and demonstrated that rayada rice was unique in genealogy and should be treated as a new (the sixth) major group of rice germplasm. With reference to the global classification of rice cultivars, we identified three major groups (indica, temperate japonica and tropical japonica) in Chinese rice germplasm and showed that Chinese temperate japonica contained higher diversity than that of global samples, whereas Chinese indica and tropical japonica maintained slightly lower diversity than that present in the global samples. Particularly, we observed that all seasonal, drought-tolerant and endosperm types occurred within each of three major groups of Chinese cultivars, which does not support previous claims that seasonal differentiation exists in Indica and drought-tolerant differentiation is present in Japonica. It is most likely that differentiation of cultivar types arose multiple times stemming from artificial selection for adaptation to local environments.