Global Analysis of Genetic, Epigenetic and Transcriptional Polymorphisms in Arabidopsis thaliana Using Whole Genome Tiling Arrays

Whole genome tiling arrays provide a high resolution platform for profiling of genetic, epigenetic, and gene expression polymorphisms. In this study we surveyed natural genomic variation in cytosine methylation among Arabidopsis thaliana wild accessions Columbia (Col) and Vancouver (Van) by comparing hybridization intensity difference between genomic DNA digested with either methylation-sensitive (HpaII) or -insensitive (MspI) restriction enzyme. Single Feature Polymorphisms (SFPs) were assayed on a full set of 1,683,620 unique features of Arabidopsis Tiling Array 1.0F (Affymetrix), while constitutive and polymorphic CG methylation were assayed on a subset of 54,519 features, which contain a 5′CCGG3′ restriction site. 138,552 SFPs (1% FDR) were identified across enzyme treatments, which preferentially accumulated in pericentromeric regions. Our study also demonstrates that at least 8% of all analyzed CCGG sites were constitutively methylated across the two strains, while about 10% of all analyzed CCGG sites were differentially methylated between the two strains. Within euchromatin arms, both constitutive and polymorphic CG methylation accumulated in central regions of genes but under-represented toward the 5′ and 3′ ends of the coding sequences. Nevertheless, polymorphic methylation occurred much more frequently in gene ends than constitutive methylation. Inheritance of methylation polymorphisms in reciprocal F1 hybrids was predominantly additive, with F1 plants generally showing levels of methylation intermediate between the parents. By comparing gene expression profiles, using matched tissue samples, we found that magnitude of methylation polymorphism immediately upstream or downstream of the gene was inversely correlated with the degree of expression variation for that gene. In contrast, methylation polymorphism within genic region showed weak positive correlation with expression variation. Our results demonstrated extensive genetic and epigenetic polymorphisms between Arabidopsis accessions and suggested a possible relationship between natural CG methylation variation and gene expression variation

Network Analysis Identifies ELF3 as a QTL for the Shade Avoidance Response in Arabidopsis

Quantitative Trait Loci (QTL) analyses in immortal populations are a powerful method for exploring the genetic mechanisms that control interactions of organisms with their environment. However, QTL analyses frequently do not culminate in the identification of a causal gene due to the large chromosomal regions often underlying QTLs. A reasonable approach to inform the process of causal gene identification is to incorporate additional genome-wide information, which is becoming increasingly accessible. In this work, we perform QTL analysis of the shade avoidance response in the Bayreuth-0 (Bay-0, CS954) x Shahdara (Sha, CS929) recombinant inbred line population of Arabidopsis. We take advantage of the complex pleiotropic nature of this trait to perform network analysis using co-expression, eQTL and functional classification from publicly available datasets to help us find good candidate genes for our strongest QTL, SAR2. This novel network analysis detected EARLY FLOWERING 3 (ELF3; AT2G25930) as the most likely candidate gene affecting the shade avoidance response in our population. Further genetic and transgenic experiments confirmed ELF3 as the causative gene for SAR2. The Bay-0 and Sha alleles of ELF3 differentially regulate developmental time and circadian clock period length in Arabidopsis, and the extent of this regulation is dependent on the light environment. This is the first time that ELF3 has been implicated in the shade avoidance response and that different natural alleles of this gene are shown to have phenotypic effects. In summary, we show that development of networks to inform candidate gene identification for QTLs is a promising technique that can significantly accelerate the process of QTL cloning

Molecular evolution of Adh and LEAFY and the phylogenetic utility of their introns in Pyrus (Rosaceae)

<p>Abstract</p> <p>Background</p> <p>The genus <it>Pyrus </it>belongs to the tribe Pyreae (the former subfamily Maloideae) of the family Rosaceae, and includes one of the most important commercial fruit crops, pear. The phylogeny of <it>Pyrus </it>has not been definitively reconstructed. In our previous efforts, the internal transcribed spacer region (ITS) revealed a poorly resolved phylogeny due to non-concerted evolution of nrDNA arrays. Therefore, introns of low copy nuclear genes (LCNG) are explored here for improved resolution. However, paralogs and lineage sorting are still two challenges for applying LCNGs in phylogenetic studies, and at least two independent nuclear loci should be compared. In this work the second intron of <it>LEAFY </it>and the alcohol dehydrogenase gene (<it>Adh</it>) were selected to investigate their molecular evolution and phylogenetic utility.</p> <p>Results</p> <p>DNA sequence analyses revealed a complex ortholog and paralog structure of <it>Adh </it>genes in <it>Pyrus </it>and <it>Malus</it>, the pears and apples. Comparisons between sequences from RT-PCR and genomic PCR indicate that some <it>Adh </it>homologs are putatively nonfunctional. A partial region of <it>Adh1 </it>was sequenced for 18 <it>Pyrus </it>species and three subparalogs representing <it>Adh1-1 </it>were identified. These led to poorly resolved phylogenies due to low sequence divergence and the inclusion of putative recombinants. For the second intron of <it>LEAFY</it>, multiple inparalogs were discovered for both <it>LFY1int2 </it>and <it>LFY2int2</it>. <it>LFY1int2 </it>is inadequate for phylogenetic analysis due to lineage sorting of two inparalogs. <it>LFY2int2-N</it>, however, showed a relatively high sequence divergence and led to the best-resolved phylogeny. This study documents the coexistence of outparalogs and inparalogs, and lineage sorting of these paralogs and orthologous copies. It reveals putative recombinants that can lead to incorrect phylogenetic inferences, and presents an improved phylogenetic resolution of <it>Pyrus </it>using <it>LFY2int2-N</it>.</p> <p>Conclusions</p> <p>Our study represents the first phylogenetic analyses based on LCNGs in <it>Pyrus</it>. Ancient and recent duplications lead to a complex structure of <it>Adh </it>outparalogs and inparalogs in <it>Pyrus </it>and <it>Malus</it>, resulting in neofunctionalization, nonfunctionalization and possible subfunctionalization. Among all investigated orthologs, <it>LFY2int2-N </it>is the best nuclear marker for phylogenetic reconstruction of <it>Pyrus </it>due to suitable sequence divergence and the absence of lineage sorting.</p

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