VIB-UGent Center for Plant Systems Biology

Advances in modern sequencing technologies allow us to generate sufficient data to analyze hundreds of bacterial genomes from a single machine in a single day. This potential for sequencing massive numbers of genomes calls for fully automated methods to produce high-quality assemblies and variant calls. We introduce Pilon, a fully automated, all-in-one tool for correcting draft assemblies and calling sequence variants of multiple sizes, including very large insertions and deletions. Pilon works with many types of sequence data, but is particularly strong when supplied with paired end data from two Illumina libraries with small e.g., 180 bp and large e.g., 3-5 Kb inserts. Pilon significantly improves draft genome assemblies by correcting bases, fixing mis-assemblies and filling gaps. For both haploid and diploid genomes, Pilon produces more contiguous genomes with fewer errors, enabling identification of more biologically relevant genes. Furthermore, Pilon identifies small variants with high accuracy as compared to state-of-the-art tools and is unique in its ability to accurately identify large sequence variants including duplications and resolve large insertions. Pilon is being used to improve the assemblies of thousands of new genomes and to identify variants from thousands of clinically relevant bacterial strains. Pilon is freely available as open source software.

This paper reports the genome sequence of domesticated tomato, a major crop plant, and a draft sequence for its closest wild relative; comparative genomics reveal very little divergence between the two genomes but some important differences with the potato genome, another important food crop in the genus Solanum. Tomato (Solanum lycopersicum) is a major crop plant and a model system for fruit development. Solanum is one of the largest angiosperm genera1 and includes annual and perennial plants from diverse habitats. Here we present a high-quality genome sequence of domesticated tomato, a draft sequence of its closest wild relative, Solanum pimpinellifolium2, and compare them to each other and to the potato genome (Solanum tuberosum). The two tomato genomes show only 0.6% nucleotide divergence and signs of recent admixture, but show more than 8% divergence from potato, with nine large and several smaller inversions. In contrast to Arabidopsis, but similar to soybean, tomato and potato small RNAs map predominantly to gene-rich chromosomal regions, including gene promoters. The Solanum lineage has experienced two consecutive genome triplications: one that is ancient and shared with rosids, and a more recent one. These triplications set the stage for the neofunctionalization of genes controlling fruit characteristics, such as colour and fleshiness.

During the past decade the pesticidal bacterium Bacillus thuringiensis has been the subject of intensive research. These efforts have yielded considerable data about the complex relationships between the structure, mechanism of action, and genetics of the organism's pesticidal crystal proteins, and a coherent picture of these relationships is beginning to emerge. Other studies have focused on the ecological role of the B. thuringiensis crystal proteins, their performance in agricultural and other natural settings, and the evolution of resistance mechanisms in target pests. Armed with this knowledge base and with the tools of modern biotechnology, researchers are now reporting promising results in engineering more-useful toxins and formulations, in creating transgenic plants that express pesticidal activity, and in constructing integrated management strategies to insure that these products are utilized with maximum efficiency and benefit.

JASPAR (http://jaspar.genereg.net/) is an open-access database containing manually curated, non-redundant transcription factor (TF) binding profiles for TFs across six taxonomic groups. In this 9th release, we expanded the CORE collection with 341 new profiles (148 for plants, 101 for vertebrates, 85 for urochordates, and 7 for insects), which corresponds to a 19% expansion over the previous release. We added 298 new profiles to the Unvalidated collection when no orthogonal evidence was found in the literature. All the profiles were clustered to provide familial binding profiles for each taxonomic group. Moreover, we revised the structural classification of DNA binding domains to consider plant-specific TFs. This release introduces word clouds to represent the scientific knowledge associated with each TF. We updated the genome tracks of TFBSs predicted with JASPAR profiles in eight organisms; the human and mouse TFBS predictions can be visualized as native tracks in the UCSC Genome Browser. Finally, we provide a new tool to perform JASPAR TFBS enrichment analysis in user-provided genomic regions. All the data is accessible through the JASPAR website, its associated RESTful API, the R/Bioconductor data package, and a new Python package, pyJASPAR, that facilitates serverless access to the data.

JASPAR (http://jaspar.genereg.net) is an open-access database of curated, non-redundant transcription factor (TF)-binding profiles stored as position frequency matrices (PFMs) and TF flexible models (TFFMs) for TFs across multiple species in six taxonomic groups. In the 2018 release of JASPAR, the CORE collection has been expanded with 322 new PFMs (60 for vertebrates and 262 for plants) and 33 PFMs were updated (24 for vertebrates, 8 for plants and 1 for insects). These new profiles represent a 30% expansion compared to the 2016 release. In addition, we have introduced 316 TFFMs (95 for vertebrates, 218 for plants and 3 for insects). This release incorporates clusters of similar PFMs in each taxon and each TF class per taxon. The JASPAR 2018 CORE vertebrate collection of PFMs was used to predict TF-binding sites in the human genome. The predictions are made available to the scientific community through a UCSC Genome Browser track data hub. Finally, this update comes with a new web framework with an interactive and responsive user-interface, along with new features. All the underlying data can be retrieved programmatically using a RESTful API and through the JASPAR 2018 R/Bioconductor package.

Conifers have dominated forests for more than 200 million years and are of huge ecological and economic importance. Here we present the draft assembly of the 20-gigabase genome of Norway spruce (Picea abies), the first available for any gymnosperm. The number of well-supported genes (28,354) is similar to the >100 times smaller genome of Arabidopsis thaliana, and there is no evidence of a recent whole-genome duplication in the gymnosperm lineage. Instead, the large genome size seems to result from the slow and steady accumulation of a diverse set of long-terminal repeat transposable elements, possibly owing to the lack of an efficient elimination mechanism. Comparative sequencing of Pinus sylvestris, Abies sibirica, Juniperus communis, Taxus baccata and Gnetum gnemon reveals that the transposable element diversity is shared among extant conifers. Expression of 24-nucleotide small RNAs, previously implicated in transposable element silencing, is tissue-specific and much lower than in other plants. We further identify numerous long (>10,000 base pairs) introns, gene-like fragments, uncharacterized long non-coding RNAs and short RNAs. This opens up new genomic avenues for conifer forestry and breeding. The draft genome of the Norway spruce (P. abies) is presented; this is the first gymnosperm genome to be sequenced and reveals a large genome size (20 Gb) resulting from the accumulation of transposable elements, and comparative sequencing of five other gymnosperm genomes provides insights into conifer genome evolution. The first draft gymnosperm genome, that of a Norway spruce (Picea abies), is published this week by the Spruce Genome Project consortium. The genome is from a tree originally collected in 1959 in eastern Jämtland, central Sweden. At 20 gigabases, the genome is more than a hundred times larger than that of the model plant species Arabidopsis, but the two contain a similar number of genes. The large genome size is the result of an accumulation of transposable elements. Comparative sequencing of five further gymnosperm genomes suggests that transposable element diversity is shared among extant conifers. The sequence data are available for public access from the ConGenIE website ( http://congenie.org/ ).

Sequencing of Medicago truncatula, a model organism of legume biology, shows that genome duplications had a role in the evolution of endosymbiotic nitrogen fixation. Legumes are unusual among plants in that they can carry out endosymbiotic nitrogen fixation with rhizobial bacteria. The genome of Medicago truncatula (also known as barrel medic or barrel clover), a well-established model for the study of legume biology, has now been sequenced. Genome analysis shows that M. truncatula has undergone several rounds of whole-genome duplication, and that the duplication that took place approximately 58 million years ago played an important part in the evolution of endosymbiotic nitrogen fixation. Legumes (Fabaceae or Leguminosae) are unique among cultivated plants for their ability to carry out endosymbiotic nitrogen fixation with rhizobial bacteria, a process that takes place in a specialized structure known as the nodule. Legumes belong to one of the two main groups of eurosids, the Fabidae, which includes most species capable of endosymbiotic nitrogen fixation1. Legumes comprise several evolutionary lineages derived from a common ancestor 60 million years ago (Myr ago). Papilionoids are the largest clade, dating nearly to the origin of legumes and containing most cultivated species2. Medicago truncatula is a long-established model for the study of legume biology. Here we describe the draft sequence of the M. truncatula euchromatin based on a recently completed BAC assembly supplemented with Illumina shotgun sequence, together capturing ∼94% of all M. truncatula genes. A whole-genome duplication (WGD) approximately 58 Myr ago had a major role in shaping the M. truncatula genome and thereby contributed to the evolution of endosymbiotic nitrogen fixation. Subsequent to the WGD, the M. truncatula genome experienced higher levels of rearrangement than two other sequenced legumes, Glycine max and Lotus japonicus. M. truncatula is a close relative of alfalfa (Medicago sativa), a widely cultivated crop with limited genomics tools and complex autotetraploid genetics. As such, the M. truncatula genome sequence provides significant opportunities to expand alfalfa’s genomic toolbox.

The filamentous fungus Aspergillus niger is widely exploited by the fermentation industry for the production of enzymes and organic acids, particularly citric acid. We sequenced the 33.9-megabase genome of A. niger CBS 513.88, the ancestor of currently used enzyme production strains. A high level of synteny was observed with other aspergilli sequenced. Strong function predictions were made for 6,506 of the 14,165 open reading frames identified. A detailed description of the components of the protein secretion pathway was made and striking differences in the hydrolytic enzyme spectra of aspergilli were observed. A reconstructed metabolic network comprising 1,069 unique reactions illustrates the versatile metabolism of A. niger. Noteworthy is the large number of major facilitator superfamily transporters and fungal zinc binuclear cluster transcription factors, and the presence of putative gene clusters for fumonisin and ochratoxin A synthesis.

Jasmonates (JAs) trigger an important transcriptional reprogramming of plant cells to modulate both basal development and stress responses. In spite of the importance of transcriptional regulation, only one transcription factor (TF), the Arabidopsis thaliana basic helix-loop-helix MYC2, has been described so far as a direct target of JAZ repressors. By means of yeast two-hybrid screening and tandem affinity purification strategies, we identified two previously unknown targets of JAZ repressors, the TFs MYC3 and MYC4, phylogenetically closely related to MYC2. We show that MYC3 and MYC4 interact in vitro and in vivo with JAZ repressors and also form homo- and heterodimers with MYC2 and among themselves. They both are nuclear proteins that bind DNA with sequence specificity similar to that of MYC2. Loss-of-function mutations in any of these two TFs impair full responsiveness to JA and enhance the JA insensitivity of myc2 mutants. Moreover, the triple mutant myc2 myc3 myc4 is as impaired as coi1-1 in the activation of several, but not all, JA-mediated responses such as the defense against bacterial pathogens and insect herbivory. Our results show that MYC3 and MYC4 are activators of JA-regulated programs that act additively with MYC2 to regulate specifically different subsets of the JA-dependent transcriptional response.

The spider mite Tetranychus urticae is a cosmopolitan agricultural pest with an extensive host plant range and an extreme record of pesticide resistance. Here we present the completely sequenced and annotated spider mite genome, representing the first complete chelicerate genome. At 90 megabases T. urticae has the smallest sequenced arthropod genome. Compared with other arthropods, the spider mite genome shows unique changes in the hormonal environment and organization of the Hox complex, and also reveals evolutionary innovation of silk production. We find strong signatures of polyphagy and detoxification in gene families associated with feeding on different hosts and in new gene families acquired by lateral gene transfer. Deep transcriptome analysis of mites feeding on different plants shows how this pest responds to a changing host environment. The T. urticae genome thus offers new insights into arthropod evolution and plant–herbivore interactions, and provides unique opportunities for developing novel plant protection strategies. The genome of the spider mite Tetranychus urticae is sequenced, providing insights into its polyphagous feeding, silk production, hormonal repertoire and reduced Hox cluster. The spider mite (Tetranychus urticae) is a common agricultural pest that feeds on a wide range of hosts — including maize (corn), soya, tomatoes and peppers — and is notoriously resistant to pesticides. Its genome has now been sequenced and analysed, providing insights into its hormonal repertoire and the evolution of silk production. Transcriptome analysis of mites feeding on different plants reveals how this pest defends itself in a changing host environment and gives pointers to possible non-pesticide plant-protection strategies. The genome encodes 17 fibroin genes, and physical tests of spider-mite silk show it to be a natural nanomaterial with fibres that are more than 100 times thinner than those produced by silk spiders.

BACKGROUND: Worldwide, grapes and their derived products have a large market. The cultivated grape species Vitis vinifera has potential to become a model for fruit trees genetics. Like many plant species, it is highly heterozygous, which is an additional challenge to modern whole genome shotgun sequencing. In this paper a high quality draft genome sequence of a cultivated clone of V. vinifera Pinot Noir is presented. PRINCIPAL FINDINGS: We estimate the genome size of V. vinifera to be 504.6 Mb. Genomic sequences corresponding to 477.1 Mb were assembled in 2,093 metacontigs and 435.1 Mb were anchored to the 19 linkage groups (LGs). The number of predicted genes is 29,585, of which 96.1% were assigned to LGs. This assembly of the grape genome provides candidate genes implicated in traits relevant to grapevine cultivation, such as those influencing wine quality, via secondary metabolites, and those connected with the extreme susceptibility of grape to pathogens. Single nucleotide polymorphism (SNP) distribution was consistent with a diffuse haplotype structure across the genome. Of around 2,000,000 SNPs, 1,751,176 were mapped to chromosomes and one or more of them were identified in 86.7% of anchored genes. The relative age of grape duplicated genes was estimated and this made possible to reveal a relatively recent Vitis-specific large scale duplication event concerning at least 10 chromosomes (duplication not reported before). CONCLUSIONS: Sanger shotgun sequencing and highly efficient sequencing by synthesis (SBS), together with dedicated assembly programs, resolved a complex heterozygous genome. A consensus sequence of the genome and a set of mapped marker loci were generated. Homologous chromosomes of Pinot Noir differ by 11.2% of their DNA (hemizygous DNA plus chromosomal gaps). SNP markers are offered as a tool with the potential of introducing a new era in the molecular breeding of grape.

The crystal proteins of Bacillus thuringiensis have been extensively studied because of their pesticidal properties and their high natural levels of production. The increasingly rapid characterization of new crystal protein genes, triggered by an effort to discover proteins with new pesticidal properties, has resulted in a variety of sequences and activities that no longer fit the original nomenclature system proposed in 1989. Bacillus thuringiensis pesticidal crystal protein (Cry and Cyt) nomenclature was initially based on insecticidal activity for the primary ranking criterion. Many exceptions to this systematic arrangement have become apparent, however, making the nomenclature system inconsistent. Additionally, the original nomenclature, with four activity-based primary ranks for 13 genes, did not anticipate the current 73 holotype sequences that form many more than the original four subgroups. A new nomenclature, based on hierarchical clustering using amino acid sequence identity, is proposed. Roman numerals have been exchanged for Arabic numerals in the primary rank (e.g., Cry1Aa) to better accommodate the large number of expected new sequences. In this proposal, 133 crystal proteins comprising 24 primary ranks are systematically arranged.

Disease suppressive soils typically develop after a disease outbreak due to the subsequent assembly of protective microbiota in the rhizosphere. The role of the plant immune system in the assemblage of a protective rhizosphere microbiome is largely unknown. In this study, we demonstrate that Arabidopsis thaliana specifically promotes three bacterial species in the rhizosphere upon foliar defense activation by the downy mildew pathogen Hyaloperonospora arabidopsidis. The promoted bacteria were isolated and found to interact synergistically in biofilm formation in vitro. Although separately these bacteria did not affect the plant significantly, together they induced systemic resistance against downy mildew and promoted growth of the plant. Moreover, we show that the soil-mediated legacy of a primary population of downy mildew infected plants confers enhanced protection against this pathogen in a second population of plants growing in the same soil. Together our results indicate that plants can adjust their root microbiome upon pathogen infection and specifically recruit a group of disease resistance-inducing and growth-promoting beneficial microbes, therewith potentially maximizing the chance of survival of their offspring that will grow in the same soil.

The genome of Ectocarpus, a model organism for brown algae, has been sequenced. Brown algae are complex photosynthetic organisms that have adapted to life in rocky coastal environments. Genome analysis sheds light on this adaptation and reveals an extended set of light-harvesting and pigment biosynthesis genes and novel metabolic processes such as halide metabolism. Comparative genomic analyses highlight the likely importance of a family of receptor kinases and related molecules in the evolution of multicellularity in plants, animals and brown algae. The genome of Ectocarpus siliculosis, a model for the study of brown algae, has been sequenced. These seaweeds are complex photosynthetic organisms that have adapted to rocky coastal environments. Genome analysis sheds light on this adaptation, revealing an extended set of light-harvesting and pigment biosynthesis genes, and new metabolic processes such as halide metabolism. Comparative analyses are also significant with respect to the evolution of multicellularity in plants, animals and brown algae. Brown algae (Phaeophyceae) are complex photosynthetic organisms with a very different evolutionary history to green plants, to which they are only distantly related1. These seaweeds are the dominant species in rocky coastal ecosystems and they exhibit many interesting adaptations to these, often harsh, environments. Brown algae are also one of only a small number of eukaryotic lineages that have evolved complex multicellularity (Fig. 1). We report the 214 million base pair (Mbp) genome sequence of the filamentous seaweed Ectocarpus siliculosus (Dillwyn) Lyngbye, a model organism for brown algae2,3,4,5, closely related to the kelps6,7 (Fig. 1). Genome features such as the presence of an extended set of light-harvesting and pigment biosynthesis genes and new metabolic processes such as halide metabolism help explain the ability of this organism to cope with the highly variable tidal environment. The evolution of multicellularity in this lineage is correlated with the presence of a rich array of signal transduction genes. Of particular interest is the presence of a family of receptor kinases, as the independent evolution of related molecules has been linked with the emergence of multicellularity in both the animal and green plant lineages. The Ectocarpus genome sequence represents an important step towards developing this organism as a model species, providing the possibility to combine genomic and genetic2 approaches to explore these and other4,5 aspects of brown algal biology further.

A new generation lignocellulose biorefinery uses heterogeneous catalysis for the high-yield production of a handful of chemicals from wood.

The differential distribution of the plant signaling molecule auxin is required for many aspects of plant development. Local auxin maxima and gradients arise as a result of local auxin metabolism and, predominantly, from directional cell-to-cell transport. In this primer, we discuss how the coordinated activity of several auxin influx and efflux systems, which transport auxin across the plasma membrane, mediates directional auxin flow. This activity crucially contributes to the correct setting of developmental cues in embryogenesis, organogenesis, vascular tissue formation and directional growth in response to environmental stimuli.

Rust fungi are some of the most devastating pathogens of crop plants. They are obligate biotrophs, which extract nutrients only from living plant tissues and cannot grow apart from their hosts. Their lifestyle has slowed the dissection of molecular mechanisms underlying host invasion and avoidance or suppression of plant innate immunity. We sequenced the 101-Mb genome of Melampsora larici-populina, the causal agent of poplar leaf rust, and the 89-Mb genome of Puccinia graminis f. sp. tritici, the causal agent of wheat and barley stem rust. We then compared the 16,399 predicted proteins of M. larici-populina with the 17,773 predicted proteins of P. graminis f. sp tritici. Genomic features related to their obligate biotrophic lifestyle include expanded lineage-specific gene families, a large repertoire of effector-like small secreted proteins, impaired nitrogen and sulfur assimilation pathways, and expanded families of amino acid and oligopeptide membrane transporters. The dramatic up-regulation of transcripts coding for small secreted proteins, secreted hydrolytic enzymes, and transporters in planta suggests that they play a role in host infection and nutrient acquisition. Some of these genomic hallmarks are mirrored in the genomes of other microbial eukaryotes that have independently evolved to infect plants, indicating convergent adaptation to a biotrophic existence inside plant cells.

BACKGROUND: Phosphorylation is the most frequent post-translational modification made to proteins and may regulate protein activity as either a molecular digital switch or a rheostat. Despite the cornucopia of high-throughput (HTP) phosphoproteomic data in the last decade, it remains unclear how many proteins are phosphorylated and how many phosphorylation sites (p-sites) can exist in total within a eukaryotic proteome. We present the first reliable estimates of the total number of phosphoproteins and p-sites for four eukaryotes (human, mouse, Arabidopsis, and yeast). RESULTS: In all, 187 HTP phosphoproteomic datasets were filtered, compiled, and studied along with two low-throughput (LTP) compendia. Estimates of the number of phosphoproteins and p-sites were inferred by two methods: Capture-Recapture, and fitting the saturation curve of cumulative redundant vs. cumulative non-redundant phosphoproteins/p-sites. Estimates were also adjusted for different levels of noise within the individual datasets and other confounding factors. We estimate that in total, 13 000, 11 000, and 3000 phosphoproteins and 230 000, 156 000, and 40 000 p-sites exist in human, mouse, and yeast, respectively, whereas estimates for Arabidopsis were not as reliable. CONCLUSIONS: Most of the phosphoproteins have been discovered for human, mouse, and yeast, while the dataset for Arabidopsis is still far from complete. The datasets for p-sites are not as close to saturation as those for phosphoproteins. Integration of the LTP data suggests that current HTP phosphoproteomics appears to be capable of capturing 70 % to 95 % of total phosphoproteins, but only 40 % to 60 % of total p-sites.

Mix of metabolites tunes root microbiota Uncharacterized biosynthetic genes in plant genomes suggest that plants make a plethora of specialized metabolites. Huang et al. reconstructed three biosynthetic networks from the small mustard plant Arabidopsis thaliana . Promiscuous acyltransferases and dehydrogenases contributed to metabolite diversification. The plant may use these specialized metabolites to modulate the microbiota surrounding its roots. Disruption of the pathways and intervention with purified compounds caused changes in the root microbiota. Science , this issue p. eaau6389

Polyploidy has been hypothesized to be both an evolutionary dead-end and a source for evolutionary innovation and species diversification. Although polyploid organisms, especially plants, abound, the apparent nonrandom long-term establishment of genome duplications suggests a link with environmental conditions. Whole-genome duplications seem to correlate with periods of extinction or global change, while polyploids often thrive in harsh or disturbed environments. Evidence is also accumulating that biotic interactions, for instance, with pathogens or mutualists, affect polyploids differently than nonpolyploids. Here, we review recent findings and insights on the effect of both abiotic and biotic stress on polyploids versus nonpolyploids and propose that stress response in general is an important and even determining factor in the establishment and success of polyploidy.