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.

The Biological Networks Gene Ontology tool (BiNGO) is an open-source Java tool to determine which Gene Ontology (GO) terms are significantly overrepresented in a set of genes. BiNGO can be used either on a list of genes, pasted as text, or interactively on subgraphs of biological networks visualized in Cytoscape. BiNGO maps the predominant functional themes of the tested gene set on the GO hierarchy, and takes advantage of Cytoscape's versatile visualization environment to produce an intuitive and customizable visual representation of the results.

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.

Diatoms, a type of microscopic marine and freshwater alga, dominate the oceans and are responsible for about a fifth of the primary productivity on Earth. The complete genome sequence of Phaeodactylum tricornutum is reported in this issue, the second diatom to be sequenced. Comparisons with Thalassiosira pseudonana, the first, reveal that hundreds of diatom genes have been acquired by gene transfer from bacteria — or vice versa. Gene transfer appears to have been common during diatom evolution, creating unorthodox combinations of genes — including some from plants and animals — likely to play major roles in nutrient management and environmental signalling. Diatoms are photosynthetic secondary endosymbionts found throughout marine and freshwater environments, and are believed to be responsible for around one-fifth of the primary productivity on Earth1,2. The genome sequence of the marine centric diatom Thalassiosira pseudonana was recently reported, revealing a wealth of information about diatom biology3,4,5. Here we report the complete genome sequence of the pennate diatom Phaeodactylum tricornutum and compare it with that of T. pseudonana to clarify evolutionary origins, functional significance and ubiquity of these features throughout diatoms. In spite of the fact that the pennate and centric lineages have only been diverging for 90 million years, their genome structures are dramatically different and a substantial fraction of genes (∼40%) are not shared by these representatives of the two lineages. Analysis of molecular divergence compared with yeasts and metazoans reveals rapid rates of gene diversification in diatoms. Contributing factors include selective gene family expansions, differential losses and gains of genes and introns, and differential mobilization of transposable elements. Most significantly, we document the presence of hundreds of genes from bacteria. More than 300 of these gene transfers are found in both diatoms, attesting to their ancient origins, and many are likely to provide novel possibilities for metabolite management and for perception of environmental signals. These findings go a long way towards explaining the incredible diversity and success of the diatoms in contemporary oceans.

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.

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 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.

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.

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.

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.

Eucalypts are the world's most widely planted hardwood trees. Their outstanding diversity, adaptability and growth have made them a global renewable resource of fibre and energy. We sequenced and assembled >94% of the 640-megabase genome of Eucalyptus grandis. Of 36,376 predicted protein-coding genes, 34% occur in tandem duplications, the largest proportion thus far in plant genomes. Eucalyptus also shows the highest diversity of genes for specialized metabolites such as terpenes that act as chemical defence and provide unique pharmaceutical oils. Genome sequencing of the E. grandis sister species E. globulus and a set of inbred E. grandis tree genomes reveals dynamic genome evolution and hotspots of inbreeding depression. The E. grandis genome is the first reference for the eudicot order Myrtales and is placed here sister to the eurosids. This resource expands our understanding of the unique biology of large woody perennials and provides a powerful tool to accelerate comparative biology, breeding and biotechnology.

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.

Being continuously exposed to variable environmental conditions, plants produce phytohormones to react quickly and specifically to these changes. The phytohormone ethylene is produced in response to multiple stresses. While the role of ethylene in defense responses to pathogens is widely recognized, recent studies in arabidopsis and crop species highlight an emerging key role for ethylene in the regulation of organ growth and yield under abiotic stress. Molecular connections between ethylene and growth-regulatory pathways have been uncovered, and altering the expression of ethylene response factors (ERFs) provides a new strategy for targeted ethylene-response engineering. Crops with optimized ethylene responses show improved growth in the field, opening new windows for future crop improvement. This review focuses on how ethylene regulates shoot growth, with an emphasis on leaves.

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.