Plant-RRBS, a bisulfite and next-generation sequencing-based methylome profiling method enriching for coverage of cytosine positions

Background: Cytosine methylation in plant genomes is important for the regulation of gene transcription and transposon activity. Genome-wide methylomes are studied upon mutation of the DNA methyltransferases, adaptation to environmental stresses or during development. However, from basic biology to breeding programs, there is a need to monitor multiple samples to determine transgenerational methylation inheritance or differential cytosine methylation. Methylome data obtained by sodium hydrogen sulfite (bisulfite)-conversion and next-generation sequencing (NGS) provide genome- wide information on cytosine methylation. However, a profiling method that detects cytosine methylation state dispersed over the genome would allow high-throughput analysis of multiple plant samples with distinct epigenetic signatures. We use specific restriction endonucleases to enrich for cytosine coverage in a bisulfite and NGS-based profiling method, which was compared to whole-genome bisulfite sequencing of the same plant material. Methods: We established an effective methylome profiling method in plants, termed plant-reduced representation bisulfite sequencing (plant-RRBS), using optimized double restriction endonuclease digestion, fragment end repair, adapter ligation, followed by bisulfite conversion, PCR amplification and NGS. We report a performant laboratory protocol and a straightforward bioinformatics data analysis pipeline for plant-RRBS, applicable for any reference-sequenced plant species. Results: As a proof of concept, methylome profiling was performed using an Oryza sativa ssp. indica pure breeding line and a derived epigenetically altered line (epiline). Plant-RRBS detects methylation levels at tens of millions of cytosine positions deduced from bisulfite conversion in multiple samples. To evaluate the method, the coverage of cytosine positions, the intra-line similarity and the differential cytosine methylation levels between the pure breeding line and the epiline were determined. Plant-RRBS reproducibly covers commonly up to one fourth of the cytosine positions in the rice genome when using MspI-DpnII within a group of five biological replicates of a line. The method predominantly detects cytosine methylation in putative promoter regions and not-annotated regions in rice. Conclusions: Plant-RRBS offers high-throughput and broad, genome- dispersed methylation detection by effective read number generation obtained from reproducibly covered genome fractions using optimized endonuclease combinations, facilitating comparative analyses of multi-sample studies for cytosine methylation and transgenerational stability in experimental material and plant breeding populations

Natural history of Arabidopsis thaliana and oomycete symbioses

Molecular ecology of plant–microbe interactions has immediate significance for filling a gap in knowledge between the laboratory discipline of molecular biology and the largely theoretical discipline of evolutionary ecology. Somewhere in between lies conservation biology, aimed at protection of habitats and the diversity of species housed within them. A seemingly insignificant wildflower called Arabidopsis thaliana has an important contribution to make in this endeavour. It has already transformed botanical research with deepening understanding of molecular processes within the species and across the Plant Kingdom; and has begun to revolutionize plant breeding by providing an invaluable catalogue of gene sequences that can be used to design the most precise molecular markers attainable for marker-assisted selection of valued traits. This review describes how A. thaliana and two of its natural biotrophic parasites could be seminal as a model for exploring the biogeography and molecular ecology of plant–microbe interactions, and specifically, for testing hypotheses proposed from the geographic mosaic theory of co-evolution

A genomic analysis of disease-resistance genes encoding nucleotide binding sites in Sorghum bicolor

A large set of candidate nucleotide-binding site (NBS)-encoding genes related to disease resistance was identified in the sorghum (Sorghum bicolor) genome. These resistance (R) genes were characterized based on their structural diversity, physical chromosomal location and phylogenetic relationships. Based on their N-terminal motifs and leucine-rich repeats (LRR), 50 non-regular NBS genes and 224 regular NBS genes were identified in 274 candidate NBS genes. The regular NBS genes were classified into ten types: CNL, CN, CNLX, CNX, CNXL, CXN, NX, N, NL and NLX. The vast majority (97%) of NBS genes occurred in gene clusters, indicating extensive gene duplication in the evolution of S. bicolor NBS genes. Analysis of the S. bicolor NBS phylogenetic tree revealed two major clades. Most NBS genes were located at the distal tip of the long arms of the ten sorghum chromosomes, a pattern significantly different from rice and Arabidopsis, the NBS genes of which have a random chromosomal distribution

PeanutMap: an online genome database for comparative molecular maps of peanut

BACKGROUND: Molecular maps have been developed for many species, and are of particular importance for varietal development and comparative genomics. However, despite the existence of multiple sets of linkage maps, databases of these data are lacking for many species, including peanut. DESCRIPTION: PeanutMap provides a web-based interface for viewing specific linkage groups of a map set. PeanutMap can display and compare multiple maps of a set based upon marker or trait correspondences, which is particularly important as cultivated peanut is a disomic tetraploid. The database can also compare linkage groups among multiple map sets, allowing identification of corresponding linkage groups from results of different research projects. Data from the two published peanut genome map sets, and also from three maps sets of phenotypic traits are present in the database. Data from PeanutMap have been incorporated into the Legume Information System website to allow peanut map data to be used for cross-species comparisons. CONCLUSION: The utility of the database is expected to increase as several SSR-based maps are being developed currently, and expanded efforts for comparative mapping of legumes are underway. Optimal use of these data will benefit from the development of tools to facilitate comparative analysis

Arabidopsis thaliana encodes a bacterial-type heterodimeric isopropylmalate isomerase involved in both Leu biosynthesis and the Met chain elongation pathway of glucosinolate formation

The last steps of the Leu biosynthetic pathway and the Met chain elongation cycle for glucosinolate formation share identical reaction types suggesting a close evolutionary relationship of these pathways. Both pathways involve the condensation of acetyl-CoA and a 2-oxo acid, isomerization of the resulting 2-malate derivative to form a 3-malate derivative, the oxidation-decarboxylation of the 3-malate derivative to give an elongated 2-oxo acid, and transamination to generate the corresponding amino acid. We have now analyzed the genes encoding the isomerization reaction, the second step of this sequence, in Arabidopsis thaliana. One gene encodes the large subunit and three encode small subunits of this enzyme, referred to as isopropylmalate isomerase (IPMI) with respect to the Leu pathway. Metabolic profiling of large subunit mutants revealed accumulation of intermediates of both Leu biosynthesis and Met chain elongation, and an altered composition of aliphatic glucosinolates demonstrating the function of this gene in both pathways. In contrast, the small subunits appear to be specialized to either Leu biosynthesis or Met chain elongation. Green fluorescent protein tagging experiments confirms the import of one of the IPMI small subunits into the chloroplast, the localization of the Met chain elongation pathway in these organelles. These results suggest the presence of different heterodimeric IPMIs in Arabidopsis chloroplasts with distinct substrate specificities for Leu or glucosinolate metabolism determined by the nature of the different small subunit

Dynamics of Vulmar/VulMITE group of transposable elements in Chenopodiaceae subfamily Betoideae

Transposable elements are important factors driving plant genome evolution. Upon their mobilization, novel insertion polymorphisms are being created. We investigated differences in copy number and insertion polymorphism of a group of Mariner-like transposable elements Vulmar and related VulMITE miniature inverted-repeat transposable elements (MITEs) in species representing subfamily Betoideae. Insertion sites of these elements were identified using a modified transposon display protocol, allowing amplification of longer fragments representing regions flanking insertion sites. Subsequently, a subset of TD fragments was converted into insertion site-based polymorphism (ISBP) markers. The investigated group of transposable elements was the most abundant in accessions representing the section Beta, showing intraspecific insertion polymorphisms likely resulting from their recent activity. In contrast, no unique insertions were observed for species of the genus Beta section Corollinae, while a set of section-specific insertions was observed in the genus Patellifolia, however, only two of them were polymorphic between P. procumbens and P. webbiana. We hypothesize that Vulmar and VulMITE elements were inactivated in the section Corollinae, while they remained active in the section Beta and the genus Patellifolia. The ISBP markers generally confirmed the insertion patterns observed with TD markers, including presence of distinct subsets of TE insertions specific to Beta and Patellifolia

Patterns of polymorphism and linkage disequilibrium in cultivated barley

We carried out a genome-wide analysis of polymorphism (4,596 SNP loci across 190 elite cultivated accessions) chosen to represent the available genetic variation in current elite North West European and North American barley germplasm. Population sub-structure, patterns of diversity and linkage disequilibrium varied considerably across the seven barley chromosomes. Gene-rich and rarely recombining haplotype blocks that may represent up to 60% of the physical length of barley chromosomes extended across the ‘genetic centromeres’. By positioning 2,132 bi-parentally mapped SNP markers with minimum allele frequencies higher than 0.10 by association mapping, 87.3% were located to within 5 cM of their original genetic map position. We show that at this current marker density genetically diverse populations of relatively small size are sufficient to fine map simple traits, providing they are not strongly stratified within the sample, fall outside the genetic centromeres and population sub-structure is effectively controlled in the analysis. Our results have important implications for association mapping, positional cloning, physical mapping and practical plant breeding in barley and other major world cereals including wheat and rye that exhibit comparable genome and genetic features

Assessing the functional coherence of modules found in multiple-evidence networks from Arabidopsis

<p>Abstract</p> <p>Background</p> <p>Combining multiple evidence-types from different information sources has the potential to reveal new relationships in biological systems. The integrated information can be represented as a relationship network, and clustering the network can suggest possible functional modules. The value of such modules for gaining insight into the underlying biological processes depends on their functional coherence. The challenges that we wish to address are to define and quantify the functional coherence of modules in relationship networks, so that they can be used to infer function of as yet unannotated proteins, to discover previously unknown roles of proteins in diseases as well as for better understanding of the regulation and interrelationship between different elements of complex biological systems.</p> <p>Results</p> <p>We have defined the functional coherence of modules with respect to the Gene Ontology (GO) by considering two complementary aspects: (i) the fragmentation of the GO functional categories into the different modules and (ii) the most representative functions of the modules. We have proposed a set of metrics to evaluate these two aspects and demonstrated their utility in <it>Arabidopsis thaliana</it>. We selected 2355 proteins for which experimentally established protein-protein interaction (PPI) data were available. From these we have constructed five relationship networks, four based on single types of data: PPI, co-expression, co-occurrence of protein names in scientific literature abstracts and sequence similarity and a fifth one combining these four evidence types. The ability of these networks to suggest biologically meaningful grouping of proteins was explored by applying Markov clustering and then by measuring the functional coherence of the clusters.</p> <p>Conclusions</p> <p>Relationship networks integrating multiple evidence-types are biologically informative and allow more proteins to be assigned to a putative functional module. Using additional evidence types concentrates the functional annotations in a smaller number of modules without unduly compromising their consistency. These results indicate that integration of more data sources improves the ability to uncover functional association between proteins, both by allowing more proteins to be linked and producing a network where modular structure more closely reflects the hierarchy in the gene ontology.</p