Finding the sources of missing heritability in a yeast cross

For many traits, including susceptibility to common diseases in humans, causal loci uncovered by genetic mapping studies explain only a minority of the heritable contribution to trait variation. Multiple explanations for this "missing heritability" have been proposed. Here we use a large cross between two yeast strains to accurately estimate different sources of heritable variation for 46 quantitative traits and to detect underlying loci with high statistical power. We find that the detected loci explain nearly the entire additive contribution to heritable variation for the traits studied. We also show that the contribution to heritability of gene-gene interactions varies among traits, from near zero to 50%. Detected two-locus interactions explain only a minority of this contribution. These results substantially advance our understanding of the missing heritability problem and have important implications for future studies of complex and quantitative traits

Natural Allelic Variation Defines a Role for ATMYC1: Trichome Cell Fate Determination

The molecular nature of biological variation is not well understood. Indeed, many questions persist regarding the types of molecular changes and the classes of genes that underlie morphological variation within and among species. Here we have taken a candidate gene approach based on previous mapping results to identify the gene and ultimately a polymorphism that underlies a trichome density QTL in Arabidopsis thaliana. Our results show that natural allelic variation in the transcription factor ATMYC1 alters trichome density in A. thaliana; this is the first reported function for ATMYC1. Using site-directed mutagenesis and yeast two-hybrid experiments, we demonstrate that a single amino acid replacement in ATMYC1, discovered in four ecotypes, eliminates known protein–protein interactions in the trichome initiation pathway. Additionally, in a broad screen for molecular variation at ATMYC1, including 72 A. thaliana ecotypes, a high-frequency block of variation was detected that results in >10% amino acid replacement within one of the eight exons of the gene. This sequence variation harbors a strong signal of divergent selection but has no measurable effect on trichome density. Homologs of ATMYC1 are pleiotropic, however, so this block of variation may be the result of natural selection having acted on another trait, while maintaining the trichome density role of the gene. These results show that ATMYC1 is an important source of variation for epidermal traits in A. thaliana and indicate that the transcription factors that make up the TTG1 genetic pathway generally may be important sources of epidermal variation in plants

High-throughput SNP genotyping in the highly heterozygous genome of Eucalyptus: assay success, polymorphism and transferability across species

<p>Abstract</p> <p>Background</p> <p>High-throughput SNP genotyping has become an essential requirement for molecular breeding and population genomics studies in plant species. Large scale SNP developments have been reported for several mainstream crops. A growing interest now exists to expand the speed and resolution of genetic analysis to outbred species with highly heterozygous genomes. When nucleotide diversity is high, a refined diagnosis of the target SNP sequence context is needed to convert queried SNPs into high-quality genotypes using the Golden Gate Genotyping Technology (GGGT). This issue becomes exacerbated when attempting to transfer SNPs across species, a scarcely explored topic in plants, and likely to become significant for population genomics and inter specific breeding applications in less domesticated and less funded plant genera.</p> <p>Results</p> <p>We have successfully developed the first set of 768 SNPs assayed by the GGGT for the highly heterozygous genome of <it>Eucalyptus </it>from a mixed Sanger/454 database with 1,164,695 ESTs and the preliminary 4.5X draft genome sequence for <it>E. grandis</it>. A systematic assessment of <it>in silico </it>SNP filtering requirements showed that stringent constraints on the SNP surrounding sequences have a significant impact on SNP genotyping performance and polymorphism. SNP assay success was high for the 288 SNPs selected with more rigorous <it>in silico </it>constraints; 93% of them provided high quality genotype calls and 71% of them were polymorphic in a diverse panel of 96 individuals of five different species.</p> <p>SNP reliability was high across nine <it>Eucalyptus </it>species belonging to three sections within subgenus Symphomyrtus and still satisfactory across species of two additional subgenera, although polymorphism declined as phylogenetic distance increased.</p> <p>Conclusions</p> <p>This study indicates that the GGGT performs well both within and across species of <it>Eucalyptus </it>notwithstanding its nucleotide diversity ≥2%. The development of a much larger array of informative SNPs across multiple <it>Eucalyptus </it>species is feasible, although strongly dependent on having a representative and sufficiently deep collection of sequences from many individuals of each target species. A higher density SNP platform will be instrumental to undertake genome-wide phylogenetic and population genomics studies and to implement molecular breeding by Genomic Selection in <it>Eucalyptus</it>.</p

Single-nucleotide polymorphism discovery by high-throughput sequencing in sorghum

<p>Abstract</p> <p>Background</p> <p>Eight diverse sorghum (<it>Sorghum bicolor </it>L. Moench) accessions were subjected to short-read genome sequencing to characterize the distribution of single-nucleotide polymorphisms (SNPs). Two strategies were used for DNA library preparation. Missing SNP genotype data were imputed by local haplotype comparison. The effect of library type and genomic diversity on SNP discovery and imputation are evaluated.</p> <p>Results</p> <p>Alignment of eight genome equivalents (6 Gb) to the public reference genome revealed 283,000 SNPs at ≥82% confirmation probability. Sequencing from libraries constructed to limit sequencing to start at defined restriction sites led to genotyping 10-fold more SNPs in all 8 accessions, and correctly imputing 11% more missing data, than from semirandom libraries. The SNP yield advantage of the reduced-representation method was less than expected, since up to one fifth of reads started at noncanonical restriction sites and up to one third of restriction sites predicted <it>in silico </it>to yield unique alignments were not sampled at near-saturation. For imputation accuracy, the availability of a genomically similar accession in the germplasm panel was more important than panel size or sequencing coverage.</p> <p>Conclusions</p> <p>A sequence quantity of 3 million 50-base reads per accession using a <it>Bsr</it>FI library would conservatively provide satisfactory genotyping of 96,000 sorghum SNPs. For most reliable SNP-genotype imputation in shallowly sequenced genomes, germplasm panels should consist of pairs or groups of genomically similar entries. These results may help in designing strategies for economical genotyping-by-sequencing of large numbers of plant accessions.</p

A Genome-Wide Association Study Identifies Variants Underlying the Arabidopsis thaliana Shade Avoidance Response

Shade avoidance is an ecologically and molecularly well-understood set of plant developmental responses that occur when the ratio of red to far-red light (R∶FR) is reduced as a result of foliar shade. Here, a genome-wide association study (GWAS) in Arabidopsis thaliana was used to identify variants underlying one of these responses: increased hypocotyl elongation. Four hypocotyl phenotypes were included in the study, including height in high R∶FR conditions (simulated sun), height in low R∶FR conditions (simulated shade), and two different indices of the response of height to low R∶FR. GWAS results showed that variation in these traits is controlled by many loci of small to moderate effect. A known PHYC variant contributing to hypocotyl height variation was identified and lists of significantly associated genes were enriched in a priori candidates, suggesting that this GWAS was capable of generating meaningful results. Using metadata such as expression data, GO terms, and other annotation, we were also able to identify variants in candidate de novo genes. Patterns of significance among our four phenotypes allowed us to categorize associations into three groups: those that affected hypocotyl height without influencing shade avoidance, those that affected shade avoidance in a height-dependent fashion, and those that exerted specific control over shade avoidance. This grouping allowed for the development of explicit hypotheses about the genetics underlying shade avoidance variation. Additionally, the response to shade did not exhibit any marked geographic distribution, suggesting that variation in low R∶FR–induced hypocotyl elongation may represent a response to local conditions

Adaptive Value of Phenological Traits in Stressful Environments: Predictions Based on Seed Production and Laboratory Natural Selection

Phenological traits often show variation within and among natural populations of annual plants. Nevertheless, the adaptive value of post-anthesis traits is seldom tested. In this study, we estimated the adaptive values of pre- and post-anthesis traits in two stressful environments (water stress and interspecific competition), using the selfing annual species Arabidopsis thaliana. By estimating seed production and by performing laboratory natural selection (LNS), we assessed the strength and nature (directional, disruptive and stabilizing) of selection acting on phenological traits in A. thaliana under the two tested stress conditions, each with four intensities. Both the type of stress and its intensity affected the strength and nature of selection, as did genetic constraints among phenological traits. Under water stress, both experimental approaches demonstrated directional selection for a shorter life cycle, although bolting time imposes a genetic constraint on the length of the interval between bolting and anthesis. Under interspecific competition, results from the two experimental approaches showed discrepancies. Estimation of seed production predicted directional selection toward early pre-anthesis traits and long post-anthesis periods. In contrast, the LNS approach suggested neutrality for all phenological traits. This study opens questions on adaptation in complex natural environment where many selective pressures act simultaneously

Genome-wide association mapping reveals a rich genetic architecture of complex traits in Oryza sativa

Asian rice, Oryza sativa is a cultivated, inbreeding species that feeds over half of the world's population. Understanding the genetic basis of diverse physiological, developmental, and morphological traits provides the basis for improving yield, quality and sustainability of rice. Here we show the results of a genome-wide association study based on genotyping 44,100 SNP variants across 413 diverse accessions of O. sativa collected from 82 countries that were systematically phenotyped for 34 traits. Using cross-population-based mapping strategies, we identified dozens of common variants influencing numerous complex traits. Significant heterogeneity was observed in the genetic architecture associated with subpopulation structure and response to environment. This work establishes an open-source translational research platform for genome-wide association studies in rice that directly links molecular variation in genes and metabolic pathways with the germplasm resources needed to accelerate varietal development and crop improvement

Genome-wide association mapping identifies a new arsenate reductase enzyme critical for limiting arsenic accumulation in plants

Inorganic arsenic is a carcinogen, and its ingestion through foods such as rice presents a significant risk to human health. Plants chemically reduce arsenate to arsenite. Using genome-wide association (GWA) mapping of loci controlling natural variation in arsenic accumulation in Arabidopsis thaliana allowed us to identify the arsenate reductase required for this reduction, which we named High Arsenic Content 1 (HAC1). Complementation verified the identity of HAC1, and expression in Escherichia coli lacking a functional arsenate reductase confirmed the arsenate reductase activity of HAC1. The HAC1 protein accumulates in the epidermis, the outer cell layer of the root, and also in the pericycle cells surrounding the central vascular tissue. Plants lacking HAC1 lose their ability to efflux arsenite from roots, leading to both increased transport of arsenic into the central vascular tissue and on into the shoot. HAC1 therefore functions to reduce arsenate to arsenite in the outer cell layer of the root, facilitating efflux of arsenic as arsenite back into the soil to limit both its accumulation in the root and transport to the shoot. Arsenate reduction by HAC1 in the pericycle may play a role in limiting arsenic loading into the xylem. Loss of HAC1-encoded arsenic reduction leads to a significant increase in arsenic accumulation in shoots, causing an increased sensitivity to arsenate toxicity. We also confirmed the previous observation that the ACR2 arsenate reductase in A. thaliana plays no detectable role in arsenic metabolism. Furthermore, ACR2 does not interact epistatically with HAC1, since arsenic metabolism in the acr2 hac1 double mutant is disrupted in an identical manner to that described for the hac1 single mutant. Our identification of HAC1 and its associated natural variation provides an important new resource for the development of low arsenic-containing food such as rice