Search

The epigenome orchestrates genome accessibility, functionality, and three-dimensional structure. Because epigenetic variation can impact transcription and thus phenotypes, it may contribute to adaptation. Here, we report 1,107 high-quality single-base resolution methylomes and 1,203 transcriptomes from the 1001 Genomes collection of Arabidopsis thaliana. Although the genetic basis of methylation variation is highly complex, geographic origin is a major predictor of genome-wide DNA methylation levels and of altered gene expression caused by epialleles. Comparison to cistrome and epicistrome datasets identifies associations between transcription factor binding sites, methylation, nucleotide variation, and co-expression modules. Physical maps for nine of the most diverse genomes reveal how transposons and other structural variants shape the epigenome, with dramatic effects on immunity genes. The 1001 Epigenomes Project provides a comprehensive resource for understanding how variation in DNA methylation contributes to molecular and non-molecular phenotypes in natural populations of the most studied model plant.

Solanum pennellii is a wild tomato species endemic to Andean regions in South America, where it has evolved to thrive in arid habitats. Because of its extreme stress tolerance and unusual morphology, it is an important donor of germplasm for the cultivated tomato Solanum lycopersicum. Introgression lines (ILs) in which large genomic regions of S. lycopersicum are replaced with the corresponding segments from S. pennellii can show remarkably superior agronomic performance. Here we describe a high-quality genome assembly of the parents of the IL population. By anchoring the S. pennellii genome to the genetic map, we define candidate genes for stress tolerance and provide evidence that transposable elements had a role in the evolution of these traits. Our work paves a path toward further tomato improvement and for deciphering the mechanisms underlying the myriad other agronomic traits that can be improved with S. pennellii germplasm.

Centromeres are critical for cell division, loading CENH3 or CENPA histone variant nucleosomes, directing kinetochore formation and allowing chromosome segregation1,2. Despite their conserved function, centromere size and structure are diverse across species. To understand this centromere paradox3,4, it is necessary to know how centromeric diversity is generated and whether it reflects ancient trans-species variation or, instead, rapid post-speciation divergence. To address these questions, we assembled 346 centromeres from 66 Arabidopsis thaliana and 2 Arabidopsis lyrata accessions, which exhibited a remarkable degree of intra- and inter-species diversity. A. thaliana centromere repeat arrays are embedded in linkage blocks, despite ongoing internal satellite turnover, consistent with roles for unidirectional gene conversion or unequal crossover between sister chromatids in sequence diversification. Additionally, centrophilic ATHILA transposons have recently invaded the satellite arrays. To counter ATHILA invasion, chromosome-specific bursts of satellite homogenization generate higher-order repeats and purge transposons, in line with cycles of repeat evolution. Centromeric sequence changes are even more extreme in comparison between A. thaliana and A. lyrata. Together, our findings identify rapid cycles of transposon invasion and purging through satellite homogenization, which drive centromere evolution and ultimately contribute to speciation.Inter- and intra-species comparison of Arabidopsis centromere variation identifies rapid cycles of transposon invasion and purging through satellite homogenization that drive centromere evolution. [ABSTRACT FROM AUTHOR]

Repeated herbicide applications in agricultural fields exert strong selection on weeds such as blackgrass (Alopecurus myosuroides), which is a major threat for temperate climate cereal crops. This inadvertent selection pressure provides an opportunity for investigating the underlying genetic mechanisms and evolutionary processes of rapid adaptation, which can occur both through mutations in the direct targets of herbicides and through changes in other, often metabolic, pathways, known as non-target-site resistance. How much target-site resistance (TSR) relies on de novo mutations vs. standing variation is important for developing strategies to manage herbicide resistance. We first generated a chromosome-level reference genome for A. myosuroides for population genomic studies of herbicide resistance and genome-wide diversity across Europe in this species. Next, through empirical data in the form of highly accurate long-read amplicons of alleles encoding acetyl-CoA carboxylase (ACCase) and acetolactate synthase (ALS) variants, we showed that most populations with resistance due to TSR mutations—23 out of 27 and six out of nine populations for ACCase and ALS, respectively—contained at least two TSR haplotypes, indicating that soft sweeps are the norm. Finally, through forward-in-time simulations, we inferred that TSR is likely to mainly result from standing genetic variation, with only a minor role for de novo mutations. [ABSTRACT FROM AUTHOR]

Plants can defend themselves against a wide array of enemies, from microbes to large animals, yet there is great variability in the effectiveness of such defences, both within and between species. Some of this variation can be explained by conflicting pressures from pathogens with different modes of attack (1). A second explanation comes from an evolutionary 'tug of war', in which pathogens adapt to evade detection, until the plant has evolved new recognition capabilities for pathogen invasion (2-5). If selection is, however, sufficiently strong, susceptible hosts should remain rare. That this is not the case is best explained by costs incurred from constitutive defences in a pest-free environment (6-11). Using a combination of forward genetics and genome-wide association analyses, we demonstrate that allelic diversity at a single locus, ACCELERATED CELL DEATH 6 (ACD6) (12,13), underpins marked pleiotropic differences in both vegetative growth and resistance to microbial infection and herbivory among natural Arabidopsis thaliana strains. A hyperactive ACD6 allele, compared to the reference allele, strongly enhances resistance to a broad range of pathogens from different phyla, but at the same time slows the production of new leaves and greatly reduces the biomass of mature leaves. This allele segregates at intermediate frequency both throughout the worldwide range of A. thaliana and within local populations, consistent with this allele providing substantial fitness benefits despite its marked impact on growth., A survey of A. thaliana accessions collected from the wild revealed extensive environment-dependent variation for leaf initiation rate (Supplementary Table 1). One of the strains, Est-1, which produced leaves more [...]

Little is known about the types of mutations underlying the evolution of species-specific traits. The metal hyperaccumulator Arabidopsis halleri has the rare ability to colonize heavy-metal-polluted soils, and, as an [...]

The gram-negative plant-pathogenic bacterium Xanthomonas campestris pv. vesicatoria is the causative agent of bacterial spot disease in pepper and tomato plants, which leads to economically important yield losses. This pathosystem has become a well-established model for studying bacterial infection strategies. Here, we present the whole-genome sequence of the pepper-pathogenic Xanthomonas campestris pv. vesicatoria strain 85-10, which comprises a 5.17-Mb circular chromosome and four plasmids. The genome has a high G+C content (64.75%) and signatures of extensive genome plasticity. Whole-genome comparisons revealed a gene order similar to both Xanthomonas axonopodis pv. citri and Xanthomonas campestris pv. campestris and a structure completely different from Xanthomonas oryzae pv. oryzae. A total of 548 coding sequences (12.2%) are unique to X. campestris pv. vesicatoria. In addition to a type III secretion system, which is essential for pathogenicity., the genome of strain 85-10 encodes all other types of protein secretion systems described so far in gram-negative bacteria. Remarkably, one of the putative type IV secretion systems encoded on the largest plasmid is similar to the Icm/Dot systems of the human pathogens Legionella pneumophila and Coxiella burnetii. Comparisons with other completely sequenced plant pathogens predicted six novel type III effector proteins and several other virulence factors, including adhesins, cell wall-degrading enzymes, and extracellular polysaccharides.

Summary: Thlaspi arvense (field pennycress) is being domesticated as a winter annual oilseed crop capable of improving ecosystems and intensifying agricultural productivity without increasing land use. It is a selfing diploid with a short life cycle and is amenable to genetic manipulations, making it an accessible field‐based model species for genetics and epigenetics. The availability of a high‐quality reference genome is vital for understanding pennycress physiology and for clarifying its evolutionary history within the Brassicaceae. Here, we present a chromosome‐level genome assembly of var. MN106‐Ref with improved gene annotation and use it to investigate gene structure differences between two accessions (MN108 and Spring32‐10) that are highly amenable to genetic transformation. We describe non‐coding RNAs, pseudogenes and transposable elements, and highlight tissue‐specific expression and methylation patterns. Resequencing of forty wild accessions provided insights into genome‐wide genetic variation, and QTL regions were identified for a seedling colour phenotype. Altogether, these data will serve as a tool for pennycress improvement in general and for translational research across the Brassicaceae. [ABSTRACT FROM AUTHOR]

To understand the evolution of developmental processes, nonmodel organisms in the nematodes, insects, and vertebrates are compared with established model systems. Often, these comparisons suffer from the inability to apply sophisticated technologies to these nonmodel species. In the nematode Pristionchus pacificus, cellular and genetic analyses are used to compare vulva development to that of Caenorhabditis elegans. However, substantial changes in gene function between P. pacificus and C. elegans limit the use of candidate gene approaches in studying P. pacificus mutations. To facilitate map-based cloning of mutations in P. pacificus, we constructed a BAC-based genetic linkage map. A BAC library of 13,440 clones was generated and completely end sequenced. By comparing BAC end and EST sequences between the 'wild-type' strain P. pacificus var. California and the polymorphic strain P. pacificus var. Washington, 133 single-stranded conformational polymorphisms were identified. These markers were tested on a meiotic mapping panel of 46 randomly picked [F.sub.2] animals after a cross of the two strains, providing' the first genetic linkage map of P. pacificus. A mapping strategy using two selected markers per chromosome was devised and the efficiency of this approach was illustrated by the mapping of the Ppa-unc-1/Twitchin gene.

Abstract Background The guppy, Poecilia reticulata, is a well-known model organism for studying inheritance and variation of male ornamental traits as well as adaptation to different river habitats. However, genomic resources for studying this important model were not previously widely available. Results With the aim of generating molecular markers for genetic mapping of the guppy, cDNA libraries were constructed from embryos and different adult organs to generate expressed sequence tags (ESTs). About 18,000 ESTs were annotated according to BLASTN and BLASTX results and the sequence information from the 3' UTRs was exploited to generate PCR primers for re-sequencing of genomic DNA from different wild type strains. By comparison of EST-linked genomic sequences from at least four different ecotypes, about 1,700 polymorphisms were identified, representing about 400 distinct genes. Two interconnected MySQL databases were built to organize the ESTs and markers, respectively. A robust phylogeny of the guppy was reconstructed, based on 10 different nuclear genes. Conclusion Our EST and marker databases provide useful tools for genetic mapping and phylogenetic studies of the guppy.

Amaranthus tuberculatus , Amaranthus hybridus , and Amaranthus palmeri are agronomically important weed species. Here, we present the most contiguous draft assemblies of these three species to date. We utilized a combination of Pacific Biosciences long-read sequencing and chromatin contact mapping information to assemble and order sequences of A. palmeri to near-chromosome-level resolution, with scaffold N50 of 20.1 Mb. To resolve the issues of heterozygosity and coassembly of alleles in diploid species, we adapted the trio binning approach to produce haplotype assemblies of A. tuberculatus and A. hybridus. This approach resulted in an improved assembly of A. tuberculatus , and the first genome assembly for A. hybridus , with contig N50s of 2.58 and 2.26 Mb, respectively. Species-specific transcriptomes and information from related species were used to predict transcripts within each assembly. Syntenic comparisons of these species and Amaranthus hypochondriacus identified sites of genomic rearrangement, including duplication and translocation, whereas genetic map construction within A. tuberculatus highlighted the need for further ordering of the A. hybridus and A. tuberculatus contigs. These multiple reference genomes will accelerate genomic studies in these species to further our understanding of weedy evolution within Amaranthus. [ABSTRACT FROM AUTHOR]

The handheld Oxford Nanopore MinION sequencer generates ultra-long reads with minimal cost and time requirements, which makes sequencing genomes at the bench feasible. Here, we sequence the gold standard Arabidopsis thaliana genome (KBS-Mac-74 accession) on the bench with the MinION sequencer, and assemble the genome using typical consumer computing hardware (4 Cores, 16 Gb RAM) into chromosome arms (62 contigs with an N50 length of 12.3 Mb). We validate the contiguity and quality of the assembly with two independent single-molecule technologies, Bionano optical genome maps and Pacific Biosciences Sequel sequencing. The new A. thaliana KBS-Mac-74 genome enables resolution of a quantitative trait locus that had previously been recalcitrant to a Sanger-based BAC sequencing approach. In summary, we demonstrate that even when the purpose is to understand complex structural variation at a single region of the genome, complete genome assembly is becoming the simplest way to achieve this goal. [ABSTRACT FROM AUTHOR]

Background: Widespread and more frequently occurring drought conditions are a consequence of global warming and increase the demand for tolerant crop varieties to feed the growing world population. A better understanding of the molecular mechanisms underlying the water deficit response of crops will enable targeted breeding strategies to develop robust cultivars. Results: In the present study, the transcriptional response of maize (Zea mays L.) primary roots to low water potentials was monitored by RNA sequencing (RNA-Seq) experiments. After 6 h and 24 h of mild (-0.2 MPa) and severe (-0.8 MPa) water deficit conditions, the primary root transcriptomes of seedlings grown under water deficit and control conditions were compared. The number of responsive genes was dependent on and increased with intensification of water deficit treatment. After short-term mild and severe water deficit 249 and 3,000 genes were differentially expressed, respectively. After a 24 h treatment the number of affected genes increased to 7,267 and 12,838 for mild and severe water deficit, respectively, including more than 80% of the short-term responsive genes. About half of the differentially expressed genes were up-regulated and maximal fold-changes increased with treatment intensity to more than 300-fold. A consensus set of 53 genes was differentially regulated independently of the nature of deficit treatment. Characterization revealed an overrepresentation of the Gene Ontology (GO) categories “oxidoreductase activity” and “heme binding” among regulated genes connecting the water deficit response to ROS metabolism. Conclusion: This study gives a comprehensive insight in water deficit responsive genes in young maize primary roots and provides a set of candidate genes that merit further genetic analyses in the future. [ABSTRACT FROM AUTHOR]

Color intensity is a relevant feature for red wine quality, which depends on the accumulation of athocyanins during berry ripening and their interaction with other compunds present in the wine. While grape color is a varietal feature, emerging spontaneous somatic variation for this trait can be selected for cultivar improvement. Here we studied the genomic origin of a darker berry clone (VN21) selected in 'Tempranillo Tinto' (TT) cultivar, which enables the production of wines with higher color intensity. Using a diploid genome assembly of TT produced following a trio binning approach from PacBio and Nanopore sequencing, and genome re-sequencing data of VN21 compared to other TT clones, we identified a 10 Mb deletion in chromosome 11 that likely affected only the L1 meristem cell layer of VN21. An RNA-seq analysis identified a general down-regulation of genes within the chimeric hemizygous segment in the berry skin of VN21. Down-regulated genes where also enriched in wax biosynthesis functional category genes, including one CER1 and one ABCG32 homologs located in the hemizygous segment. SEM images showed that wax accumulation is impaired in the berry cuticle of VN21, which likely leads to the shiny color of VN21 berries. Candidate loss of function polymorphisms remaing hemizygous in VN21 after the chimeric deletion were also detected from TT haplotype comparison. Our findings show that very large hemizygous deletions can stabilize as periclinal chimeras in grapevine clones, giving rise to pleiotropic mutant phenotypes that can be exploited for cultivar innovation. [ABSTRACT FROM AUTHOR]

Unlike the situation in animals, the final morphology of the plant body is highly modulated by the environment. During Arabidopsis development, intrinsic factors provide the framework for basic patterning processes. CLASS III HOMEODOMAIN LEUCINE ZIPPER (HD-ZIPIII) transcription factors are involved in embryo, shoot and root patterning. During vegetative growth HD-ZIPIII proteins control several polarity set-up processes such as in leaves and the vascular system. We have identified several direct target genes of the HD-ZIPIII transcription factor REVOLUTA (REV) using a chromatin immunoprecipitation/DNA sequencing (ChIP-Seq) approach. This analysis revealed that REV acts upstream of auxin biosynthesis and affects directly the expression of several class II HD-ZIP transcription factors that have been shown to act in the shade-avoidance response pathway. We show that, as well as involvement in basic patterning, HD-ZIPIII transcription factors have a critical role in the control of the elongation growth that is induced when plants experience shade. Leaf polarity is established by the opposed actions of HD-ZIPIII and KANADI transcription factors. Finally, our study reveals that the module that consists of HD-ZIPIII/KANADI transcription factors controls shade growth antagonistically and that this antagonism is manifested in the opposed regulation of shared target genes. [ABSTRACT FROM AUTHOR]

Cohesin is a protein complex that forms a ring around sister chromatids thus holding them together. The ring is composed of three proteins: Smc1, Smc3 and Scc1. The roles of three additional proteins that associate with the ring, Scc3, Pds5 and Wpl1, are not well understood. It has been proposed that these three factors form a complex that stabilizes the ring and prevents it from opening. This activity promotes sister chromatid cohesion but at the same time poses an obstacle for the initial entrapment of sister DNAs. This hindrance to cohesion establishment is overcome during DNA replication via acetylation of the Smc3 subunit by the Eco1 acetyltransferase. However, the full mechanistic consequences of Smc3 acetylation remain unknown. In the current work, we test the requirement of Scc3 and Pds5 for the stable association of cohesin with DNA. We investigated the consequences of Scc3 and Pds5 depletion in vivo using degron tagging in budding yeast. The previously described DHFR--based N-terminal degron as well as a novel Eco1-derived C-terminal degron were employed in our study. Scc3 and Pds5 associate with cohesin complexes independently of each other and require the Scc1 "core" ubunit for their association with chromosomes. Contrary to previous data for Scc1 downregulation, depletion of either Scc3 or Pds5 had a strong effect on sister chromatid cohesion but not on cohesin binding to DNA. Quantity, stability and genome-wide distribution of cohesin complexes remained mostly unchanged after the depletion of Scc3 and Pds5. Our findings are inconsistent with a previously proposed model that Scc3 and Pds5 are cohesin maintenance factors required for cohesin ring stability or for maintaining its association with DNA. We propose that Scc3 and Pds5 specifically function during cohesion establishment in S phase. [ABSTRACT FROM AUTHOR]

The plant Arabidopsis thaliana occurs naturally in many different habitats throughout Eurasia. As a foundation for identifying genetic variation contributing to adaptation to diverse environments, a 1001 Genomes Project to sequence geographically diverse A. thaliana strains has been initiated. Here we present the first phase of this project, based on population-scale sequencing of 80 strains drawn from eight regions throughout the species' native range. We describe the majority of common small-scale polymorphisms as well as many larger insertions and deletions in the A. thaliana pan-genome, their effects on gene function, and the patterns of local and global linkage among these variants. The action of processes other than spontaneous mutation is identified by comparing the spectrum of mutations that have accumulated since A. thaliana diverged from its closest relative 10 million years ago with the spectrum observed in the laboratory. Recent species-wide selective sweeps are rare, and potentially deleterious mutations are more common in marginal populations. [ABSTRACT FROM AUTHOR]

Here, we report on the complete genome sequence of the hyperthermophilic Crenarchaeum Thermoproteus tenax (strain Kra1, DSM 2078T) a type strain of the crenarchaeotal order Thermoproteales. Its circular 1.84-megabase genome harbors no extrachromosomal elements and 2,051 open reading frames are identified, covering 90.6% of the complete sequence, which represents a high coding density. Derived from the gene content, T. tenax is a representative member of the Crenarchaeota. The organism is strictly anaerobic and sulfur-dependent with optimal growth at 86°C and pH 5.6. One particular feature is the great metabolic versatility, which is not accompanied by a distinct increase of genome size or information density as compared to other Crenarchaeota. T. tenax is able to grow chemolithoautotrophically (CO2/H2) as well as chemoorganoheterotrophically in presence of various organic substrates. All pathways for synthesizing the 20 proteinogenic amino acids are present. In addition, two presumably complete gene sets for NADH:quinone oxidoreductase (complex I) were identified in the genome and there is evidence that either NADH or reduced ferredoxin might serve as electron donor. Beside the typical archaeal A0A1-ATP synthase, a membrane-bound pyrophosphatase is found, which might contribute to energy conservation. Surprisingly, all genes required for dissimilatory sulfate reduction are present, which is confirmed by growth experiments. Mentionable is furthermore, the presence of two proteins (ParA family ATPase, actin-like protein) that might be involved in cell division in Thermoproteales, where the ESCRT system is absent, and of genes involved in genetic competence (DprA, ComF) that is so far unique within Archaea. [ABSTRACT FROM AUTHOR]

The S locus, a single polymorphic locus, is responsible for self-incompatibility (SI) in the Brassicaceae family and many related plant families. Despite its importance, our knowledge of S-locus evolution is largely restricted to the causal genes encoding the S-locus receptor kinase (SRK) receptor and S-locus cysteine-rich protein (SCR) ligand of the SI system. Here, we present high- quality sequences of the genomic region of six S-locus haplotypes: Arabidopsis (Arabidopsis thaliana; one haplotype), Arabidopsis lyrata (four haptotypes), and Capsella rubella (one haplotype). We compared these with reference S-locus haplotypes of the self- compatible Arabidopsis and its SI congener A. lyrata. We subsequently reconstructed the likely genomic organization of the S locus in the most recent common ancestor of Arabidopsis and Capsella. As previously reported, the two SI-determining genes, SCR and SRK, showed a pattern of coevolution. In addition, consistent with previous studies, we found that duplication, gene conversion, and positive selection have been important factors in the evolution of these two genes and appear to contribute to the generation of new recognition specificities. Intriguingly, the inactive pseudo-S-locus haplotype in the self-compatible species C. rubella is likely to be an old S-locus haplotype that only very recently became fixed when C. rubella split off from its SI ancestor, Capsella grandiflora. [ABSTRACT FROM AUTHOR]