Understanding quantitative genetic variation

Sex-ratio evolution in nuclear-cytoplasmic gynodioecy when restoration is a threshold trait

Gynodioecious plant species, which have populations consisting of female and hermaphrodite individuals, usually have complex sex determination involving cytoplasmic male sterility (CMS) alleles interacting with nuclear restorers of fertility. In response to recent evidence, we present a model of sex-ratio evolution in which restoration of male fertility is a threshold trait. We find that females are maintained at low frequencies for all biologically relevant parameter values. Furthermore, this model predicts periodically high female frequencies (>50%) under conditions of lower female seed fecundity advantages (compensation, x = 5%) and pleiotropic fitness effects associated with restorers of fertility (costs of restoration, y = 20%) than in other models. This model explains the maintenance of females in species that have previously experienced invasions of CMS alleles and the evolution of multiple restorers. Sensitivity of the model to small changes in cost and compensation values and to initial conditions may explain why populations of the same species vary widely for sex ratio.

Adaptation in a spider mite population after long-term evolution on a single host plant

Evolution in a single environment is expected to erode genetic variability, thereby precluding adaptation to novel environments. To test this, a large population of spider mites kept on cucumber for approximately 300 generations was used to establish populations on novel host plants (tomato or pepper), and changes in traits associated to adaptation were measured after 15 generations. Using a half-sib design, we investigated whether trait changes were related to genetic variation in the base population. Juvenile survival and fecundity exhibited genetic variation and increased in experimental populations on novel hosts. Conversely, no variation was detected for host choice and developmental time and these traits did not evolve. Longevity remained unchanged on novel hosts despite the presence of genetic variation, suggesting weak selection for this trait. Hence, patterns of evolutionary changes generally matched those of genetic variation, and changes in some traits were not hindered by long-term evolution in a constant environment.

Genome partitioning of genetic variation for height from 11,214 sibling pairs

Height has been used for more than a century as a model by which to understand quantitative genetic variation in humans. We report that the entire genome appears to contribute to its additive genetic variance. We used genotypes and phenotypes of 11,214 sibling pairs from three countries to partition additive genetic variance across the genome. Using genome scans to estimate the proportion of the genomes of each chromosome from siblings that were identical by descent, we estimated the heritability of height contributed by each of the 22 autosomes and the X chromosome. We show that additive genetic variance is spread across multiple chromosomes and that at least six chromosomes (i.e., 3, 4, 8, 15, 17, and 18) are responsible for the observed variation. Indeed, the data are not inconsistent with a uniform spread of trait loci throughout the genome. Our estimate of the variance explained by a chromosome is correlated with the number of times suggestive or significant linkage with height has been reported for that chromosome. Variance due to dominance was not significant but was difficult to assess because of the high sampling correlation between additive and dominance components. Results were consistent with the absence of any large between-chromosome epistatic effects. Notwithstanding the proposed architecture of complex traits that involves widespread gene-gene and gene-environment interactions, our results suggest that variation in height in humans can be explained by many loci distributed over all autosomes, with an additive mode of gene action.

Maintenance of genetic variation in sexual ornaments: a review of the mechanisms

Female preferences for elaborate male sexual traits have been documented in a number of species in which males contribute only genes to the next generation. In such systems, mate choice has been hypothesised to benefit females genetically. For the genetic benefits to be possible there must be additive genetic variation (V A) for sexual ornaments, such that highly ornamented males can pass fitter genes on to the progeny of choosy females. Here, I review the mechanisms that can contribute to the maintenance of this variation. The variation may be limited to sexual ornaments, resulting in Fisherian benefits in terms of the increased reproductive success of male progeny produced by choosy females. Alternatively, ornaments may capture V A in other life-history traits. In the latter case, "good genes" benefits may apply in terms of improved performance of the progeny of either sex. Some mechanisms, however, such as negative pleiotropy, sexually antagonistic variation or overdominance, can maintain V A in ornaments and other life-history traits with little variation in total fitness, leaving little room for any genetic benefits of mate choice. Distinguishing between these mechanisms has consequences not only for the theory of sexual selection, but also for evolution of sex and for biological conservation. I discuss how the traditional ways of testing for genetic benefits can usefully be supplemented by tests detecting benefits resulting from specific mechanisms maintaining V A in sexual ornaments.

Prediction of individual genetic risk to disease from genome-wide association studies

Empirical studies suggest that the effect sizes of individual causal risk alleles underlying complex genetic diseases are small, with most genotype relative risks in the range of 1.1-2.0. Although the increased risk of disease for a carrier is small for any single locus, knowledge of multiple-risk alleles throughout the genome could allow the identification of individuals that are at high risk. In this study, we investigate the number and effect size of risk loci that underlie complex disease constrained by the disease parameters of prevalence and heritability. Then we quantify the value of prediction of genetic risk to disease using a range of realistic combinations of the number, size, and distribution of risk effects that underlie complex diseases. We propose an approach to assess the genetic risk of a disease in healthy individuals, based on dense genome-wide SNP panels. We test this approach using simulation. When the number of loci contributing to the disease is >50, a large case-control study is needed to identify a set of risk loci for use in predicting the disease risk of healthy people not included in the case-control study. For diseases controlled by 1000 loci of mean relative risk of only 1.04, a case-control study with 10,000 cases and controls can lead to selection of approximately 75 loci that explain >50% of the genetic variance. The 5% of people with the highest predicted risk are three to seven times more likely to suffer the disease than the population average, depending on heritability and disease prevalence. Whether an individual with known genetic risk develops the disease depends on known and unknown environmental factors.

Differentiation of tundra/taiga and boreal coniferous forest wolves: genetics, coat colour and association with migratory caribou

The grey wolf has one of the largest historic distributions of any terrestrial mammal and can disperse over great distances across imposing topographic barriers. As a result, geographical distance and physical obstacles to dispersal may not be consequential factors in the evolutionary divergence of wolf populations. However, recent studies suggest ecological features can constrain gene flow. We tested whether wolf-prey associations in uninterrupted tundra and forested regions of Canada explained differences in migratory behaviour, genetics, and coat colour of wolves. Satellite-telemetry data demonstrated that tundra wolves (n = 19) migrate annually with caribou (n = 19) from denning areas in the tundra to wintering areas south of the treeline. In contrast, nearby boreal coniferous forest wolves are territorial and associated year round with resident prey. Spatially explicit analysis of 14 autosomal microsatellite loci (n = 404 individuals) found two genetic clusters corresponding to tundra vs. boreal coniferous forest wolves. A sex bias in gene flow was inferred based on higher levels of mtDNA divergence (F(ST) = 0.282, 0.028 and 0.033; P < 0.0001 for mitochondrial, nuclear autosomal and Y-chromosome markers, respectively). Phenotypic differentiation was substantial as 93% of wolves from tundra populations exhibited light colouration whereas only 38% of boreal coniferous forest wolves did (chi(2) = 64.52, P < 0.0001). The sharp boundary representing this discontinuity was the southern limit of the caribou migration. These findings show that substantial genetic and phenotypic differentiation in highly mobile mammals can be caused by prey-habitat specialization rather than distance or topographic barriers. The presence of a distinct wolf ecotype in the tundra of North America highlights the need to preserve migratory populations.

Reconciling strong stabilizing selection with the maintenance of genetic variation in a natural population of black field crickets (Teleogryllus commodus)

Genetic variation in single traits, including those closely related to fitness, is pervasive and generally high. By contrast, theory predicts that several forms of selection, including stabilizing selection, will eliminate genetic variation. Stabilizing selection in natural populations tends to be stronger than that assumed in theoretical models of the maintenance of genetic variation. The widespread presence of genetic variation in the presence of strong stabilizing selection is a persistent problem in evolutionary genetics that currently has no compelling explanation. The recent insight that stabilizing selection often acts most strongly on trait combinations via correlational selection may reconcile this problem. Here we show that for a set of male call properties in the cricket Teleogryllus commodus, the pattern of multivariate stabilizing sexual selection is closely associated with the degree of additive genetic variance. The multivariate trait combinations experiencing the strongest stabilizing selection harbored very little genetic variation while combinations under weak selection contained most of the genetic variation. Our experiment provides empirical support for the prediction that a small number of trait combinations experiencing strong stabilizing selection will have reduced genetic variance and that genetically independent trait combinations experiencing weak selection can simultaneously harbor much higher levels of genetic variance.

Robust associations of four new chromosome regions from genome-wide analyses of type 1 diabetes

The Wellcome Trust Case Control Consortium (WTCCC) primary genome-wide association (GWA) scan on seven diseases, including the multifactorial autoimmune disease type 1 diabetes (T1D), shows associations at P < 5 x 10(-7) between T1D and six chromosome regions: 12q24, 12q13, 16p13, 18p11, 12p13 and 4q27. Here, we attempted to validate these and six other top findings in 4,000 individuals with T1D, 5,000 controls and 2,997 family trios independent of the WTCCC study. We confirmed unequivocally the associations of 12q24, 12q13, 16p13 and 18p11 (P(follow-up) Collapse

Key Words Collapse

MESH Headings

• Adolescent
• Case-Control Studies
• Chromosome Mapping
• Diabetes Mellitus, Type 1/genetics
• Genetic Predisposition to Disease
• Genome, Human
• Humans
• Polymorphism, Single Nucleotide

Collapse

Grants

• 061859 Wellcome Trust
• 089989 Wellcome Trust
• G0000934 Medical Research Council
• Wellcome Trust
• 061858 Wellcome Trust

Collapse

Affiliation(s)

John A Todd Juvenile Diabetes Research Foundation/Wellcome Trust Diabetes and Inflammation Laboratory, Department of Medical Genetics, University of Cambridge, Addenbrooke's Hospital, Cambridge CB2 0XY, UK.
Neil M Walker
Jason D Cooper
Deborah J Smyth
Kate Downes
Vincent Plagnol
Rebecca Bailey
Sergey Nejentsev
Sarah F Field
Felicity Payne
Christopher E Lowe
Jeffrey S Szeszko
Jason P Hafler
Lauren Zeitels
Jennie H M Yang
Adrian Vella
Sarah Nutland
Helen E Stevens
Helen Schuilenburg
Gillian Coleman
Meeta Maisuria
William Meadows
Luc J Smink
Barry Healy
Oliver S Burren
Alex A C Lam
Nigel R Ovington
James Allen
Ellen Adlem
Hin-Tak Leung
Chris Wallace
Joanna M M Howson
Cristian Guja
Constantin Ionescu-Tîrgovişte
Matthew J Simmonds
Joanne M Heward
Stephen C L Gough
David B Dunger
Linda S Wicker
David G Clayton

Collapse

Sexually antagonistic genetic variation for fitness in red deer

Evolutionary theory predicts the depletion of genetic variation in natural populations as a result of the effects of selection, but genetic variation is nevertheless abundant for many traits that are under directional or stabilizing selection. Evolutionary geneticists commonly try to explain this paradox with mechanisms that lead to a balance between mutation and selection. However, theoretical predictions of equilibrium genetic variance under mutation-selection balance are usually lower than the observed values, and the reason for this is unknown. The potential role of sexually antagonistic selection in maintaining genetic variation has received little attention in this debate, surprisingly given its potential ubiquity in dioecious organisms. At fitness-related loci, a given genotype may be selected in opposite directions in the two sexes. Such sexually antagonistic selection will reduce the otherwise-expected positive genetic correlation between male and female fitness. Both theory and experimental data suggest that males and females of the same species may have divergent genetic optima, but supporting data from wild populations are still scarce. Here we present evidence for sexually antagonistic fitness variation in a natural population, using data from a long-term study of red deer (Cervus elaphus). We show that male red deer with relatively high fitness fathered, on average, daughters with relatively low fitness. This was due to a negative genetic correlation between estimates of fitness in males and females. In particular, we show that selection favours males that carry low breeding values for female fitness. Our results demonstrate that sexually antagonistic selection can lead to a trade-off between the optimal genotypes for males and females; this mechanism will have profound effects on the operation of selection and the maintenance of genetic variation in natural populations.

Knockdown resistance to heat stress and slow recovery from chill coma are genetically associated in a quantitative trait locus region of chromosome 2 inDrosophila melanogaster

In insects, two ecologically relevant traits of thermal adaptation are knockdown resistance to high temperature (KRHT) and chill-coma recovery (CCR). Chromosome 2 of Drosophila melanogaster was tested for quantitative trait loci (QTL) affecting both CCR and KRHT in backcrosses between homosequential lines that are fixed for the standard (noninverted) sequence of this autosome. These lines were obtained by artificial selection on KRHT and subsequent inbreeding from a stock that was derived from a single wild population. Heat-induced expression of the 70KD heat-shock protein (Hsp70) was also examined for variation between the lines. Composite interval mapping was performed for each trait on each reciprocal backcross, identifying one QTL region in the middle of chromosome 2 for both KRHT and CCR. The largest estimates of additive effects were found in pericentromeric regions of chromosome 2, accounting for 10-14% (CCR) and 10-17% (KRHT) of the phenotypic variance in BC populations. No QTL was found in the region of the heat-shock factor (hsf) gene. However, the two parental lines have diverged in the heat-induced Hsp70 expression. Distribution of KRHT QTL on chromosome 2 was similar between this study based on crosses between lines selected from a single wild population and previous work based on crosses between selection lines from different continents. Colocalized QTL showed a trade-off association between CCR and KRHT, which should be the result of either multiple, tightly linked trait-specific genes or a single gene with pleiotropic effects on the traits. We discuss candidate loci contained within the QTL regions.

Selection responses of means and inbreeding depression for female fecundity in Drosophila melanogaster suggest contributions from intermediate-frequency alleles to quantitative trait variation

The extent to which quantitative trait variability is caused by rare alleles maintained by mutation, versus intermediate-frequency alleles maintained by balancing selection, is an unsolved problem of evolutionary genetics. We describe the results of an experiment to examine the effects of selection on the mean and extent of inbreeding depression for early female fecundity in Drosophila melanogaster. Theory predicts that rare, partially recessive deleterious alleles should cause a much larger change in the effect of inbreeding than in the mean of the outbred population, with the change in inbreeding effect having an opposite sign to the change in mean. The present experiment fails to support this prediction, suggesting that intermediate-frequency alleles contribute substantially to genetic variation in early fecundity.

Complexity of type 2 diabetes mellitus data sets emerging from nutrigenomic research: a case for dimensionality reduction?

Nutrigenomics promises personalized nutrition and an improvement in preventing, delaying, and reducing the symptoms of chronic diseases such as diabetes. Nutritional genomics is the study of how foods affect the expression of genetic information in an individual and how an individual's genetic makeup affects the metabolism and response to nutrients and other bioactive components in food. The path to those promises has significant challenges, from experimental designs that include analysis of genetic heterogeneity to the complexities of food and environmental factors. One of the more significant complications in developing the knowledge base and potential applications is how to analyze high-dimensional datasets of genetic, nutrient, metabolomic (clinical), and other variables influencing health and disease processes. Type 2 diabetes mellitus (T2DM) is used as an illustration of the challenges in studying complex phenotypes with nutrigenomics concepts and approaches.

Linkage and segregation analysis of black and brindle coat color in domestic dogs

Mutations of pigment type switching have provided basic insight into melanocortin physiology and evolutionary adaptation. In all vertebrates that have been studied to date, two key genes, Agouti and Melanocortin 1 receptor (Mc1r), encode a ligand-receptor system that controls the switch between synthesis of red-yellow pheomelanin vs. black-brown eumelanin. However, in domestic dogs, historical studies based on pedigree and segregation analysis have suggested that the pigment type-switching system is more complicated and fundamentally different from other mammals. Using a genomewide linkage scan on a Labrador x greyhound cross segregating for black, yellow, and brindle coat colors, we demonstrate that pigment type switching is controlled by an additional gene, the K locus. Our results reveal three alleles with a dominance order of black (K(B)) > brindle (k(br)) > yellow (k(y)), whose genetic map position on dog chromosome 16 is distinct from the predicted location of other pigmentation genes. Interaction studies reveal that Mc1r is epistatic to variation at Agouti or K and that the epistatic relationship between Agouti and K depends on the alleles being tested. These findings suggest a molecular model for a new component of the melanocortin signaling pathway and reveal how coat-color patterns and pigmentary diversity have been shaped by recent selection.

The genetic basis of zinc tolerance in the metallophyte Arabidopsis halleri ssp. halleri (Brassicaceae): an analysis of quantitative trait loci

The species Arabidopsis halleri, an emerging model for the study of heavy metal tolerance and accumulation in plants, has evolved a high level of constitutive zinc tolerance. Mapping of quantitative trait loci (QTL) was used to investigate the genetic architecture of zinc tolerance in this species. A first-generation backcross progeny of A. halleri ssp. halleri from a highly contaminated industrial site and its nontolerant relative A. lyrata ssp. petraea was produced and used for QTL mapping of zinc tolerance. A genetic map covering most of the A. halleri genome was constructed using 85 markers. Among these markers, 65 were anchored in A. thaliana and revealed high synteny with other Arabidopsis genomes. Three QTL of comparable magnitude on three different linkage groups were identified. At all QTL positions zinc tolerance was enhanced by A. halleri alleles, indicating directional selection for higher zinc tolerance in this species. The two-LOD support intervals associated with these QTL cover 24, 4, and 13 cM. The importance of each of these three regions is emphasized by their colocalization with HMA4, MTP1-A, and MTP1-B, respectively, three genes well known to be involved in metal homeostasis and tolerance in plants.

Joint estimates of quantitative trait locus effect and frequency using synthetic recombinant populations of Drosophila melanogaster

We develop and implement a strategy to map QTL in two synthetic populations of Drosophila melanogaster each initiated with eight inbred founder strains. These recombinant populations allow simultaneous estimates of QTL location, effect, and frequency. Five X-linked QTL influencing bristle number were resolved to intervals of approximately 1.3 cM. We confirm previous observations of bristle number QTL distal to 4A at the tip of the chromosome and identify two novel QTL in 7F-8C, an interval that does not include any classic bristle number candidate genes. If QTL at the tip of the X are biallelic they appear to be intermediate in frequency, although there is evidence that these QTL may reside in multiallelic haplotypes. Conversely, the two QTL mapping to the middle of the X chromosome are likely rare: in each case the minor allele is observed in only 1 of the 16 founders. Assuming additivity and biallelism we estimate that identified QTL contribute 1.0 and 8.7%, respectively, to total phenotypic variation in male abdominal and sternopleural bristle number in nature. Models that seek to explain the maintenance of genetic variation make different predictions about the population frequency of QTL alleles. Thus, mapping QTL in eight-way recombinant populations can distinguish between these models.

THE GENETIC ARCHITECTURE OF IMMUNE DEFENSE AND REPRODUCTION IN MALE BOMBUS TERRESTRIS BUMBLEBEES

Understanding the architecture of genetic variation, that is the number, effect, location, and interaction, of genes responsible for phenotypic variability in nature is important for the understanding of microevolutionary processes. In this study, we have used a quantitative trait loci (QTL) approach to uncover the genetic architecture of fitness-relevant traits associated with reproduction and immune defense in male Bombus terrestris bumblebees. Three male reproductive investment traits, the number and length of the produced sperm and the size of the accessory glands, were studied. Two branches of the insect innate immune system, the activation of the Phenoloxidase-cascade and the hemolymph's antibacterial activity, were investigated. We found that variation in most of the studied traits is based on a network of minor QTLs and epistatic interactions. Unexpectedly, there was no evidence for phenotypic or genetic trade-offs between the presumably costly investment in immune defense and reproductive effort in this population for the measured traits. In fact, we found a positive correlation, both, in phenotype and genetic architecture for the number of produced sperm and antibacterial activity against an insect pathogen. A major finding for all traits analyzed was that the epistatic interactions accounted for a major proportion of the explained phenotypic variance. Especially for traits involved in immune defense, this pattern highlights the possible role of parasites in the evolution and maintenance of recombination and sexual reproduction.

Prediction of phenotype and gene expression for combinations of mutations

Molecular interactions provide paths for information flows. Genetic interactions reveal active information flows and reflect their functional consequences. We integrated these complementary data types to model the transcription network controlling cell differentiation in yeast. Genetic interactions were inferred from linear decomposition of gene expression data and were used to direct the construction of a molecular interaction network mediating these genetic effects. This network included both known and novel regulatory influences, and predicted genetic interactions. For corresponding combinations of mutations, the network model predicted quantitative gene expression profiles and precise phenotypic effects. Multiple predictions were tested and verified.

Distributions of epistasis in microbes fit predictions from a fitness landscape model

How do the fitness effects of several mutations combine? Despite its simplicity, this question is central to the understanding of multilocus evolution. Epistasis (the interaction between alleles at different loci), especially epistasis for fitness traits such as reproduction and survival, influences evolutionary predictions "almost whenever multilocus genetics matters". Yet very few models have sought to predict epistasis, and none has been empirically tested. Here we show that the distribution of epistasis can be predicted from the distribution of single mutation effects, based on a simple fitness landscape model. We show that this prediction closely matches the empirical measures of epistasis that have been obtained for Escherichia coli and the RNA virus vesicular stomatitis virus. Our results suggest that a simple fitness landscape model may be sufficient to quantitatively capture the complex nature of gene interactions. This model may offer a simple and widely applicable alternative to complex metabolic network models, in particular for making evolutionary predictions.

Limits to evolution at range margins: when and why does adaptation fail?

What stops populations expanding into new territory beyond the edge of a range margin? Recent models addressing this problem have brought together population genetics and population ecology, and some have included interactions among species at range edges. Here, we review these models of adaptation at environmental or parapatric margins, and discuss the contrasting effects of migration in either swamping local adaptation, or supplying the genetic variation that is necessary for adaptation to continue. We illustrate how studying adaptation at range margins (both with and without hybridization) can provide insight into the genetic and ecological factors that limit evolution more generally, especially in response to current rates of environmental change.

The effect of non-additive genetic interactions on selection in multi-locus genetic models

Additive genetic variance might usually be expected to decrease in a finite population because of genetic drift. However, both theoretical and empirical studies have shown that the additive genetic variance of a population could, in some cases, actually increase owing to the action of genetic drift in presence of non-additive effects. We used Monte-Carlo simulations to address a less-well-studied issue: the effects of directional truncation selection on a trait affected by non-additive genetic variation. We investigated the effects on genetic variance and the response to selection. We compared two different genetic models, representing various numbers of loci. We found that the additive genetic variance could also increase in the case of truncation selection, when dominance and epistasis was present. Additive-by-additive epistatic effects generally gave a higher increase in additive variance compared to dominance. However, the magnitude of the increase differed depending on the particular model and on the number of loci.

Variation and constraint in plant evolution and development

The goal of this short review is to consider the interrelated phenomena of phenotypic variation and genetic constraint with respect to plant diversity. The unique aspects of plants, including sessile habit, modular growth and diverse developmental programs expressed at the phytomer level, merit a specific examination of the genetic basis of their phenotypic variation, and how they experience and escape genetic constraint. Numerous QTL studies with wild and domesticated plants reveal that most phenotypic traits are polygenic but vary in the number and effect of the loci contributing, from a few loci of large effects to many with small effects. Further, somatic mutations, developmental plasticity and epigenetic variation, especially gene methylation, can contribute to increases in phenotypic variation. The flip side of these processes, genetic constraint, can similarly be the result of many factors, including pleiotropy, canalization and genetic redundancy. Genetic constraint is not only a mechanism to prevent change, however, it can also serve to direct evolution along certain paths. Ultimately, genetic constraint often comes full circle and is released through events such as hybridization, genome duplication and epigenetic remodeling. We are just beginning to understand how these processes can operate simultaneously during the evolution of ecologically important traits in plants.

Genetic architecture of leaf ecophysiological traits in Helianthus

We investigated quantitative trait loci (QTLs) for several leaf chemistry traits in early-generation hybrids between Helianthus annuus and Helianthus petiolaris, the parental species of the ancient diploid hybrid sunflower species Helianthus anomalus, Helianthus deserticola, and Helianthus paradoxus. We grew individuals of a second-generation backcross (BC(2)) toward H. petiolaris under optimum conditions in a glass house experiment. Trait values were measured once for each individual. In addition, genotypic data previously determined for each individual were employed for composite interval mapping of QTLs. We detected QTLs for leaf carbon concentration, leaf nitrogen concentration, leaf nitrogen per unit area, and photosynthetic nitrogen use efficiency. Leaf carbon isotope discrimination (delta(13)C) and leaf nitrogen isotopic composition (delta(15)N) were analyzed, but no significant QTLs were found for these traits. Interestingly, two neighboring loci explained a relatively large percentage of the variation in leaf nitrogen per unit area. This was notable because leaf nitrogen has been shown to strongly affect the fitness of early-generation sunflower hybrids in the H. anomalus habitat, and QTLs of large effect are expected to respond relatively quickly to selection. We speculate that the genetic architecture underlying leaf nitrogen may have facilitated the colonization of active desert sand dunes by H. anomalus.

Accurate, precise modeling of cell proliferation kinetics from time-lapse imaging and automated image analysis of agar yeast culture arrays

Background

Genome-wide mutant strain collections have increased demand for high throughput cellular phenotyping (HTCP). For example, investigators use HTCP to investigate interactions between gene deletion mutations and additional chemical or genetic perturbations by assessing differences in cell proliferation among the collection of 5000 S. cerevisiae gene deletion strains. Such studies have thus far been predominantly qualitative, using agar cell arrays to subjectively score growth differences. Quantitative systems level analysis of gene interactions would be enabled by more precise HTCP methods, such as kinetic analysis of cell proliferation in liquid culture by optical density. However, requirements for processing liquid cultures make them relatively cumbersome and low throughput compared to agar. To improve HTCP performance and advance capabilities for quantifying interactions, YeastXtract software was developed for automated analysis of cell array images.

Results

YeastXtract software was developed for kinetic growth curve analysis of spotted agar cultures. The accuracy and precision for image analysis of agar culture arrays was comparable to OD measurements of liquid cultures. Using YeastXtract, image intensity vs. biomass of spot cultures was linearly correlated over two orders of magnitude. Thus cell proliferation could be measured over about seven generations, including four to five generations of relatively constant exponential phase growth. Spot area normalization reduced the variation in measurements of total growth efficiency. A growth model, based on the logistic function, increased precision and accuracy of maximum specific rate measurements, compared to empirical methods. The logistic function model was also more robust against data sparseness, meaning that less data was required to obtain accurate, precise, quantitative growth phenotypes.

Conclusion

Microbial cultures spotted onto agar media are widely used for genotype-phenotype analysis, however quantitative HTCP methods capable of measuring kinetic growth rates have not been available previously. YeastXtract provides objective, automated, quantitative, image analysis of agar cell culture arrays. Fitting the resulting data to a logistic equation-based growth model yields robust, accurate growth rate information. These methods allow the incorporation of imaging and automated image analysis of cell arrays, grown on solid agar media, into HTCP-driven experimental approaches, such as global, quantitative analysis of gene interaction networks.

On the synthesis and interpretation of consistent but weak gene-disease associations in the era of genome-wide association studies

Emerging technologies are allowing researchers to study hundreds of thousands of genetic variants simultaneously as risk factors for common complex diseases. Both theoretical considerations and empirical evidence suggest that specific genetic variants causally associated with common diseases will have small effects (risk ratios mostly <2.0). However, the combination of even a few small effects (e.g. effects of fewer than 20 common genetic variants) could account for a sizeable population attributable fraction of common diseases and shed important light on disease pathogenesis and environmental determinants. Nevertheless, the inauguration of genome-wide association studies only magnifies the challenge of differentiating between the expected, true weak associations from the numerous spurious effects caused by misclassification, confounding and significance-chasing biases. Standards are urgently needed for presenting and interpreting cumulative evidence on gene-disease associations, especially for consistent but weak associations. Criteria for synthesis of the evidence should include sound methods for study conduct and analysis, biological plausibility, experimental evidence and adequate replication in large-scale, collaborative studies. Efforts by the Human Genome Epidemiology Network (HuGENet) are currently ongoing to streamline and operationalize these criteria for data on genetic associations with common diseases.

Congenic strains confirm aerobic running capacity quantitative trait loci on rat chromosome 16 and identify possible intermediate phenotypes

We previously identified two inbred rat strains divergent for treadmill aerobic running capacity (ARC), the low-performing Copenhagen (COP) and the high-performing DA rats, and used an F(2)(COPxDA) population to identify ARC quantitative trait loci (QTLs) on rat chromosome 16 (RNO16) and the proximal portion of rat chromosome 3 (RNO3). Two congenic rat strains were bred to further investigate these ARC QTLs by introgressing RNO16 and the proximal portion of RNO3 from DA rats into the genetic background of COP rats and were named COP.DA(chr 16) and COP.DA(chr 3), respectively. COP.DA(chr 16) rats had significantly greater ARC compared with COP rats (696.7 +/- 38.2 m vs. 571.9 +/- 27.5 m, P = 0.03). COP.DA(chr 3) rats had increased, although not significant, ARC compared with COP rats (643.6 +/- 40.9 m vs. 571.9 +/- 27.5 m). COP.DA(chr 16) rats had significantly greater subcutaneous abdominal fat, as well as decreased fasting triglyceride levels, compared with COP rats (P < 0.05), indicating that genes responsible for strain differences in fat metabolism are also located on RNO16. While this colocalization of QTLs may be coincidental, it is also possible that these differences in energy balance may be associated with the superior running performance of COP.DA(chr 16) consomic rats.

Genetics of global gene expression

A new field of genetic analysis of global gene expression has emerged in recent years, driven by the realization that traditional techniques of linkage and association analysis can be applied to thousands of transcript levels measured by microarrays. Genetic dissection of transcript abundance has shed light on the architecture of quantitative traits, provided a new approach for connecting DNA sequence variation with phenotypic variation, and improved our understanding of transcriptional regulation and regulatory variation.

Development of a linkage map and mapping of phenotypic polymorphisms in a free-living population of Soay sheep (Ovis aries)

An understanding of the determinants of trait variation and the selective forces acting on it in natural populations would give insights into the process of evolution. The combination of long-term studies of individuals living in the wild and better genomic resources for nonmodel organisms makes achieving this goal feasible. This article reports the development of a complete linkage map in a pedigree of free-living Soay sheep on St. Kilda and its application to mapping the loci responsible for three morphological polymorphisms for which the maintenance of variation demands explanation. The map was derived from 251 microsatellite and four allozyme markers and covers 3350 cM (approximately 90% of the sheep genome) at approximately 15-cM intervals. Marker order was consistent with the published sheep map with the exception of one region on chromosome 1 and one on chromosome 12. Coat color maps to chromosome 2 where a strong candidate gene, tyrosinase-related protein 1 (TYRP1), has also been mapped. Coat pattern maps to chromosome 13, close to the candidate locus Agouti. Horn type maps to chromosome 10, a location similar to that previously identified in domestic sheep. These findings represent an advance in the dissection of the genetic diversity in the wild and provide the foundation for QTL analyses in the study population.

Cheating can stabilize cooperation in mutualisms

Mutualisms present a challenge for evolutionary theory. How is cooperation maintained in the face of selection for selfishness and cheating? Both theory and data suggest that partner choice, where one species preferentially directs aid to the more cooperative members of the other species, is central to cooperation in many mutualisms. However, the theory has only so far considered the evolutionary effects of partner choice on one of the species in a mutualism in isolation. Here, we investigate the co-evolution of cooperation and choice in a choosy host and its symbiont. Our model reveals that even though choice and cooperation may be initially selected, it will often be unstable. This is because choice reduces variation in the symbiont and, therefore, tends to remove the selective incentive for its own maintenance (a scenario paralleled in the lek paradox in female choice and policing in within-species cooperation). However, we also show that when variability is reintroduced into symbionts each generation, in the form of less cooperative individuals, choice is maintained. This suggests that the presence of cheaters and cheater species in many mutualisms is central to the maintenance of partner choice and, paradoxically, cooperation itself.

Dissecting the complex genetic basis of mate choice

The genetic analysis of mate choice is fraught with difficulties. Males produce complex signals and displays that can consist of a combination of acoustic, visual, chemical and behavioural phenotypes. Furthermore, female preferences for these male traits are notoriously difficult to quantify. During mate choice, genes not only affect the phenotypes of the individual they are in, but can influence the expression of traits in other individuals. How can genetic analyses be conducted to encompass this complexity? Tighter integration of classical quantitative genetic approaches with modern genomic technologies promises to advance our understanding of the complex genetic basis of mate choice.

Case-control studies in the genomic era: a clinician's guide

The goal of case-control association studies is to find genetic variants in the human genome that influence common traits. The Human Genome and HapMap projects have added fresh impetus to this goal by cataloguing the raw genetic data behind human DNA variation. Studies that associate these genetic variants with phenotype improve both molecular diagnostics and drug discovery and offer clinicians important opportunities to improve care of patients. In this review I focus on case-control studies, which are the most widely used design and expected to be the most powerful. I also address the problem of case-control non-replication, which is widespread despite enormous effort and use of resources. Important causes of non-replication include inadequate statistical power to detect small and moderate effects, phenotype heterogeneity, population stratification, publication bias, and multiple comparison testing.

Artificial Selection on Male Longevity Influences Age‐Dependent Reproductive Effort in the Black Field Cricket Teleogryllus commodus

Although the trade-off between reproductive effort and longevity is central to both sexual selection and evolutionary theories of aging, there has been little synthesis between these fields. Here, we selected directly on adult longevity of male field crickets Teleogryllus commodus and measured the correlated responses of age-dependent male reproductive effort, female lifetime fecundity, and several other life-history traits. Male longevity responded significantly to five generations of divergent selection. Males from downward-selected lines commenced calling sooner and reached their peak calling effort at a younger age. They called more per night and, despite living less than half as long, called more overall than males selected for increased longevity. Females from the downward-selected lines lived significantly shorter lives than females from the upward-selected lines but still produced the same number of offspring. Nymph survival, development time, and body size and weight at eclosion did not show significant correlated response to selection on male longevity, despite evidence for substantial genetic variation in each of these traits. Collectively, our findings directly support the antagonistic pleiotropy model of aging and suggest an important role for sexual selection in the aging process.

Genetic architecture of the dog: sexual size dimorphism and functional morphology

Purebred dogs are a valuable resource for genetic analysis of quantitative traits. Quantitative traits are complex, controlled by many genes that are contained within regions of the genome known as quantitative trait loci (QTL). The genetic architecture of quantitative traits is defined by the characteristics of these genes: their number, the magnitude of their effects, their positions in the genome and their interactions with each other. QTL analysis is a valuable tool for exploring genetic architecture, and highlighting regions of the genome that contribute to the variation of a trait within a population.

How to be early flowering: an evolutionary perspective

In wild and cultivated annual plant species, flowering time is an important life-history trait that coordinates the life cycle with local environmental conditions. Extensive studies on the genetic basis of flowering time in the model species Arabidopsis thaliana have revealed a complex genetic network that can detect environmental and internal signals. Based on this knowledge and on known pleiotropic effects associated with flowering time genes, we suggest that a natural shift towards an early-flowering life cycle might involve only particular functional regions in a limited number of genes. Our predictions are supported by genetic theories of adaptation and by recent data about genes associated with natural variation. We analyse the extent to which these predictions can also apply to crop species.

Natural genetic variation in whole-genome expression in Arabidopsis thaliana: the impact of physiological QTL introgression

A long-standing and fundamental question in biology is how genes influence complex phenotypes. Combining near-isogenic line mapping with genome expression profiling offers a unique opportunity for exploring the functional relationship between genotype and phenotype and for generating candidate genes for future study. We used a whole-genome microarray produced with ink-jet technology to measure the relative expression level of over 21,500 genes from an Arabidopsis thaliana near-isogenic line (NIL) and its recurrent parent. The NIL material contained two introgressions (bottom of chromosome II and top of chromosome III) of the Cvi-1 ecotype in a Ler-2 ecotype genome background. Each introgression 'captures' a Cvi allele of a physiological quantitative trait loci (QTL) that our previous studies have shown increases transpiration and reduces water-use efficiency at the whole-plant level. We used a mixed model anova framework for assessing sources of expression variability and for evaluating statistical significance in our array experiment. We discovered 25 differentially expressed genes in the introgression at a false-discovery rate (FDR) cut-off of 0.20 and identified new candidate genes for both QTL regions. Several differentially expressed genes were confirmed with QRT-PCR (quantitative reverse transcription-polymerase chain reaction) assays. In contrast, we found no statistically significant differentially expressed genes outside of the QTL introgressions after controlling for multiple tests. We discuss these results in the context of candidate genes, cloning QTL, and phenotypic evolution.

Additive genetic variation and the distribution of QTN effects among sites

Early quantitative genetic theories emphasize the distribution of gene effects at individual loci or the distribution of mutant effects in maintaining quantitative genetic variation, but infrequently stress the distribution of gene effects among loci. In this study, we analyse the effects of the distribution of quantitative trait nucleotides (QTN) effects among sites under artificial and stabilizing selection. Wright's formula is applied to describing the density distribution of allele frequencies of multiple diallelic QTN at the equilibrium of mutation-drift-selection. Our results demonstrate that the distribution of QTN effects among sites can affect additive genetic variation in terms of total additive variance, average gene diversity, per-class contribution of QTN effects and per-QTN contribution. When the distribution of QTN effects among sites is changed from L-shaped to bell-shaped or to be a flatter, both the total additive variance and the average gene diversity are changed. Per-class and per-QTN contributions exhibit different distribution patterns. The L-shaped distribution indicates the predominant role of the aggregative effects from the QTN of small finite effects. The bell-shaped or flatter distributions indicate the predominance of the QTN of intermediate and large effects. These predictions highlight the significance of the distribution of QTN effects among sites in interpreting the maintenance of quantitative genetic variation at the fine genome scale.

Rangewide molecular structuring in the Utah sucker (Catostomus ardens)

The Utah sucker (Catostomus ardens) is endemic to the Bonneville Basin and the upper Snake River drainage in western North America, and is thought to hybridize with the federally endangered June sucker (Chasmistes liorus mictus) in Utah Lake (Bonneville Basin). Here we describe the discovery of a major subdivision in Utah suckers (4.5% mitochondrial sequence divergence) between the ancient Snake River drainage and the Bonneville Basin. This boundary has not previously been recognized in Utah suckers based on morphologic variation, but has been recently described in two endemic cyprinids in the region. Populations in valleys east of the Wasatch Mountains in Utah clustered with the Snake River populations, suggesting that these valleys may have had an ancient hydrologic connection to the Snake River. We also found evidence of population isolation within the Bonneville Basin, corresponding to two Pleistocene sub-basins of the ancient Lake Bonneville. In contrast, we found no molecular evidence for deep divergence between Utah suckers and June suckers in Utah Lake or for a history of hybridization between divergent lineages in that population, although we recognize that demographic events may have obscured this signal. These findings suggest that the morphological differences between Utah and June suckers in Utah Lake may be the result of strong, and relatively recent, ecological selection. In summary, morphological and molecular characters seem to vary along different axes in different portions of the range of this taxon, providing an interesting system for studying the contributions of neutral and adaptive variation to species diversity.

Natural variation in MAM within and between populations of Arabidopsis lyrata determines glucosinolate phenotype

The genetic variation that underlies the glucosinolate phenotype of Arabidopsis lyrata ssp. petraea was investigated between and within populations. A candidate glucosinolate biosynthetic locus (MAM, containing methylthioalkylmalate synthase genes) was mapped in A. lyrata to a location on linkage group 6 corresponding to the homologous location for MAM in A. thaliana. In A. thaliana MAM is responsible for side chain elongation in aliphatic glucosinolates, and the MAM phenotype can be characterized by the ratios of long- to short-chain glucosinolates. A quantitative trait loci (QTL) analysis of glucosinolate ratios in an A. lyrata interpopulation cross found one QTL at MAM. Additional QTL were identified for total indolic glucosinolates and for the ratio of aliphatic to indolic glucosinolates. MAM was then used as the candidate gene for a within-population cosegregation analysis in a natural A. lyrata population from Germany. Extensive variation in microsatellite markers at MAM was found and this variation cosegregated with the same glucosinolate ratios as in the QTL study. The combined results indicate that both between- and within-population genetic variation in the MAM region determines phenotypic variation in glucosinolate side chains in A. lyrata.

Studying phenotypic evolution using multivariate quantitative genetics

Quantitative genetics provides a powerful framework for studying phenotypic evolution and the evolution of adaptive genetic variation. Central to the approach is G, the matrix of additive genetic variances and covariances. G summarizes the genetic basis of the traits and can be used to predict the phenotypic response to multivariate selection or to drift. Recent analytical and computational advances have improved both the power and the accessibility of the necessary multivariate statistics. It is now possible to study the relationships between G and other evolutionary parameters, such as those describing the mutational input, the shape and orientation of the adaptive landscape, and the phenotypic divergence among populations. At the same time, we are moving towards a greater understanding of how the genetic variation summarized by G evolves. Computer simulations of the evolution of G, innovations in matrix comparison methods, and rapid development of powerful molecular genetic tools have all opened the way for dissecting the interaction between allelic variation and evolutionary process. Here I discuss some current uses of G, problems with the application of these approaches, and identify avenues for future research.

Sexual reproduction selects for robustness and negative epistasis in artificial gene networks

The mutational deterministic hypothesis for the origin and maintenance of sexual reproduction posits that sex enhances the ability of natural selection to purge deleterious mutations after recombination brings them together into single genomes. This explanation requires negative epistasis, a type of genetic interaction where mutations are more harmful in combination than expected from their separate effects. The conceptual appeal of the mutational deterministic hypothesis has been offset by our inability to identify the mechanistic and evolutionary bases of negative epistasis. Here we show that negative epistasis can evolve as a consequence of sexual reproduction itself. Using an artificial gene network model, we find that recombination between gene networks imposes selection for genetic robustness, and that negative epistasis evolves as a by-product of this selection. Our results suggest that sexual reproduction selects for conditions that favour its own maintenance, a case of evolution forging its own path.