Development of a knowledge graph framework to ease and empower translational approaches in plant research: a use-case on grain legumes.

While the continuing decline in genotyping and sequencing costs has largely benefited plant research, some key species for meeting the challenges of agriculture remain mostly understudied. As a result, heterogeneous datasets for different traits are available for a significant number of these species. As gene structures and functions are to some extent conserved through evolution, comparative genomics can be used to transfer available knowledge from one species to another. However, such a translational research approach is complex due to the multiplicity of data sources and the non-harmonized description of the data. Here, we provide two pipelines, referred to as structural and functional pipelines, to create a framework for a NoSQL graph-database (Neo4j) to integrate and query heterogeneous data from multiple species. We call this framework Orthology-driven knowledge base framework for translational research (Ortho_KB). The structural pipeline builds bridges across species based on orthology. The functional pipeline integrates biological information, including QTL, and RNA-sequencing datasets, and uses the backbone from the structural pipeline to connect orthologs in the database. Queries can be written using the Neo4j Cypher language and can, for instance, lead to identify genes controlling a common trait across species. To explore the possibilities offered by such a framework, we populated Ortho_KB to obtain OrthoLegKB, an instance dedicated to legumes. The proposed model was evaluated by studying the conservation of a flowering-promoting gene. Through a series of queries, we have demonstrated that our knowledge graph base provides an intuitive and powerful platform to support research and development programmes.

Evidence that a common arbuscular mycorrhizal network alleviates phosphate shortage in interconnected walnut sapling and maize plants.

Under agroforestry practices, inter-specific facilitation between tree rows and cultivated alleys occurs when plants increase the growth of their neighbors especially under nutrient limitation. Owing to a coarse root architecture limiting soil inorganic phosphate (Pi) uptake, walnut trees (Juglans spp.) exhibit dependency on soil-borne symbiotic arbuscular mycorrhizal fungi that extend extra-radical hyphae beyond the root Pi depletion zone. To investigate the benefits of mycorrhizal walnuts in alley cropping, we experimentally simulated an agroforestry system in which walnut rootstocks RX1 (J. regia x J. microcarpa) were connected or not by a common mycelial network (CMN) to maize plants grown under two contrasting Pi levels. Mycorrhizal colonization parameters showed that the inoculum reservoir formed by inoculated walnut donor saplings allowed the mycorrhization of maize recipient roots. Relative to non-mycorrhizal plants and whatever the Pi supply, CMN enabled walnut saplings to access maize Pi fertilization residues according to significant increases in biomass, stem diameter, and expression of JrPHT1;1 and JrPHT1;2, two mycorrhiza-inducible phosphate transporter candidates here identified by phylogenic inference of orthologs. In the lowest Pi supply, stem height, leaf Pi concentration, and biomass of RX1 were significantly higher than in non-mycorrhizal controls, showing that mycorrhizal connections between walnut and maize roots alleviated Pi deficiency in the mycorrhizal RX1 donor plant. Under Pi limitation, maize recipient plants also benefited from mycorrhization relative to controls, as inferred from larger stem diameter and height, biomass, leaf number, N content, and Pi concentration. Mycorrhization-induced Pi uptake generated a higher carbon cost for donor walnut plants than for maize plants by increasing walnut plant photosynthesis to provide the AM fungus with carbon assimilate. Here, we show that CMN alleviates Pi deficiency in co-cultivated walnut and maize plants, and may therefore contribute to limit the use of chemical P fertilizers in agroforestry systems.

Genetic analysis of global faba bean diversity, agronomic traits and selection signatures.

KEY MESSAGE: We identified marker-trait associations for key faba bean agronomic traits and genomic signatures of selection within a global germplasm collection. Faba bean (Vicia faba L.) is a high-protein grain legume crop with great potential for sustainable protein production. However, little is known about the genetics underlying trait diversity. In this study, we used 21,345 high-quality SNP markers to genetically characterize 2678 faba bean genotypes. We performed genome-wide association studies of key agronomic traits using a seven-parent-MAGIC population and detected 238 significant marker-trait associations linked to 12 traits of agronomic importance. Sixty-five of these were stable across multiple environments. Using a non-redundant diversity panel of 685 accessions from 52 countries, we identified three subpopulations differentiated by geographical origin and 33 genomic regions subjected to strong diversifying selection between subpopulations. We found that SNP markers associated with the differentiation of northern and southern accessions explained a significant proportion of agronomic trait variance in the seven-parent-MAGIC population, suggesting that some of these traits were targets of selection during breeding. Our findings point to genomic regions associated with important agronomic traits and selection, facilitating faba bean genomics-based breeding.

The giant diploid faba genome unlocks variation in a global protein crop.

Increasing the proportion of locally produced plant protein in currently meat-rich diets could substantially reduce greenhouse gas emissions and loss of biodiversity1. However, plant protein production is hampered by the lack of a cool-season legume equivalent to soybean in agronomic value2. Faba bean (Vicia faba L.) has a high yield potential and is well suited for cultivation in temperate regions, but genomic resources are scarce. Here, we report a high-quality chromosome-scale assembly of the faba bean genome and show that it has expanded to a massive 13 Gb in size through an imbalance between the rates of amplification and elimination of retrotransposons and satellite repeats. Genes and recombination events are evenly dispersed across chromosomes and the gene space is remarkably compact considering the genome size, although with substantial copy number variation driven by tandem duplication. Demonstrating practical application of the genome sequence, we develop a targeted genotyping assay and use high-resolution genome-wide association analysis to dissect the genetic basis of seed size and hilum colour. The resources presented constitute a genomics-based breeding platform for faba bean, enabling breeders and geneticists to accelerate the improvement of sustainable protein production across the Mediterranean, subtropical and northern temperate agroecological zones.

Assembly of the 81.6 Mb centromere of pea chromosome 6 elucidates the structure and evolution of metapolycentric chromosomes.

Centromeres in the legume genera Pisum and Lathyrus exhibit unique morphological characteristics, including extended primary constrictions and multiple separate domains of centromeric chromatin. These so-called metapolycentromeres resemble an intermediate form between monocentric and holocentric types, and therefore provide a great opportunity for studying the transitions between different types of centromere organizations. However, because of the exceedingly large and highly repetitive nature of metapolycentromeres, highly contiguous assemblies needed for these studies are lacking. Here, we report on the assembly and analysis of a 177.6 Mb region of pea (Pisum sativum) chromosome 6, including the 81.6 Mb centromere region (CEN6) and adjacent chromosome arms. Genes, DNA methylation profiles, and most of the repeats were uniformly distributed within the centromere, and their densities in CEN6 and chromosome arms were similar. The exception was an accumulation of satellite DNA in CEN6, where it formed multiple arrays up to 2 Mb in length. Centromeric chromatin, characterized by the presence of the CENH3 protein, was predominantly associated with arrays of three different satellite repeats; however, five other satellites present in CEN6 lacked CENH3. The presence of CENH3 chromatin was found to determine the spatial distribution of the respective satellites during the cell cycle. Finally, oligo-FISH painting experiments, performed using probes specifically designed to label the genomic regions corresponding to CEN6 in Pisum, Lathyrus, and Vicia species, revealed that metapolycentromeres evolved via the expansion of centromeric chromatin into neighboring chromosomal regions and the accumulation of novel satellite repeats. However, in some of these species, centromere evolution also involved chromosomal translocations and centromere repositioning.

A major-effect genetic locus, ApRVII, controlling resistance against both adapted and non-adapted aphid biotypes in pea.

KEY MESSAGE: A genome-wide association study for pea resistance against a pea-adapted biotype and a non-adapted biotype of the aphid, Acyrthosiphon pisum, identified a genomic region conferring resistance to both biotypes. In a context of reduced insecticide use, the development of cultivars resistant to insect pests is crucial for an integrated pest management. Pea (Pisum sativum) is a crop of major importance among cultivated legumes, for the supply of dietary proteins and nitrogen in low-input cropping systems. However, yields of the pea crop have become unstable due to plant parasites. The pea aphid (Acyrthosiphon pisum) is an insect pest species forming a complex of biotypes, each one adapted to feed on one or a few related legume species. This study aimed to identify resistance to A. pisum and the underlying genetic determinism by examining a collection of 240 pea genotypes. The collection was screened against a pea-adapted biotype and a non-adapted biotype of A. pisum to characterize their resistant phenotype. Partial resistance was observed in some pea genotypes exposed to the pea-adapted biotype. Many pea genotypes were completely resistant to non-adapted biotype, but some exhibited partial susceptibility. A genome-wide association study, using pea exome-capture sequencing data, enabled the identification of the major-effect quantitative trait locus ApRVII on the chromosome 7. ApRVII includes linkage disequilibrium blocks significantly associated with resistance to one or both of the two aphid biotypes studied. Finally, we identified candidate genes underlying ApRVII that are potentially involved in plant-aphid interactions and marker haplotypes linked with aphid resistance. This study sets the ground for the functional characterization of molecular pathways involved in pea defence to the aphids but also is a step forward for breeding aphid-resistant cultivars.

Proteomics of developing pea seeds reveals a complex antioxidant network underlying the response to sulfur deficiency and water stress.

Pea is a legume crop producing protein-rich seeds and is increasingly in demand for human consumption and animal feed. The aim of this study was to explore the proteome of developing pea seeds at three key stages covering embryogenesis, the transition to seed-filling, and the beginning of storage-protein synthesis, and to investigate how the proteome was influenced by S deficiency and water stress, applied either separately or combined. Of the 3184 proteins quantified by shotgun proteomics, 2473 accumulated at particular stages, thus providing insights into the proteome dynamics at these stages. Differential analyses in response to the stresses and inference of a protein network using the whole proteomics dataset identified a cluster of antioxidant proteins (including a glutathione S-transferase, a methionine sulfoxide reductase, and a thioredoxin) possibly involved in maintaining redox homeostasis during early seed development and preventing cellular damage under stress conditions. Integration of the proteomics data with previously obtained transcriptomics data at the transition to seed-filling revealed the transcriptional events associated with the accumulation of the stress-regulated antioxidant proteins. This transcriptional defense response involves genes of sulfate homeostasis and assimilation, thus providing candidates for targeted studies aimed at dissecting the signaling cascade linking S metabolism to antioxidant processes in developing seeds.

A reference genome for pea provides insight into legume genome evolution.

We report the first annotated chromosome-level reference genome assembly for pea, Gregor Mendel's original genetic model. Phylogenetics and paleogenomics show genomic rearrangements across legumes and suggest a major role for repetitive elements in pea genome evolution. Compared to other sequenced Leguminosae genomes, the pea genome shows intense gene dynamics, most likely associated with genome size expansion when the Fabeae diverged from its sister tribes. During Pisum evolution, translocation and transposition differentially occurred across lineages. This reference sequence will accelerate our understanding of the molecular basis of agronomically important traits and support crop improvement.

Water stress combined with sulfur deficiency in pea affects yield components but mitigates the effect of deficiency on seed globulin composition.

Water stress and sulfur (S) deficiency are two constraints increasingly faced by crops due to climate change and low-input agricultural practices. To investigate their interaction in the grain legume pea (Pisum sativum), sulfate was depleted at the mid-vegetative stage and a moderate 9-d water stress period was imposed during the early reproductive phase. The combination of the stresses impeded reproductive processes in a synergistic manner, reducing seed weight and seed number, and inducing seed abortion, which highlighted the paramount importance of sulfur for maintaining seed yield components under water stress. On the other hand, the moderate water stress mitigated the negative effect of sulfur deficiency on the accumulation of S-rich globulins (11S) in seeds, probably due to a lower seed sink strength for nitrogen, enabling a readjustment of the ratio of S-poor (7S) to 11S globulins. Transcriptome analysis of developing seeds at the end of the combined stress period indicated that similar biological processes were regulated in response to sulfur deficiency and to the combined stress, but that the extent of the transcriptional regulation was greater under sulfur deficiency. Seeds from plants subjected to the combined stresses showed a specific up-regulation of a set of transcription factor and SUMO ligase genes, indicating the establishment of unique regulatory processes when sulfur deficiency is combined with water stress.

Whole-genome landscape of Medicago truncatula symbiotic genes.

Advances in deciphering the functional architecture of eukaryotic genomes have been facilitated by recent breakthroughs in sequencing technologies, enabling a more comprehensive representation of genes and repeat elements in genome sequence assemblies, as well as more sensitive and tissue-specific analyses of gene expression. Here we show that PacBio sequencing has led to a substantially improved genome assembly of Medicago truncatula A17, a legume model species notable for endosymbiosis studies1, and has enabled the identification of genome rearrangements between genotypes at a near-base-pair resolution. Annotation of the new M. truncatula genome sequence has allowed for a thorough analysis of transposable elements and their dynamics, as well as the identification of new players involved in symbiotic nodule development, in particular 1,037 upregulated long non-coding RNAs (lncRNAs). We have also discovered that a substantial proportion (~35% and 38%, respectively) of the genes upregulated in nodules or expressed in the nodule differentiation zone colocalize in genomic clusters (270 and 211, respectively), here termed symbiotic islands. These islands contain numerous expressed lncRNA genes and display differentially both DNA methylation and histone marks. Epigenetic regulations and lncRNAs are therefore attractive candidate elements for the orchestration of symbiotic gene expression in the M. truncatula genome.

Development of two major resources for pea genomics: the GenoPea 13.2K SNP Array and a high-density, high-resolution consensus genetic map.

Single nucleotide polymorphism (SNP) arrays represent important genotyping tools for innovative strategies in both basic research and applied breeding. Pea is an important food, feed and sustainable crop with a large (about 4.45 Gbp) but not yet available genome sequence. In the present study, 12 pea recombinant inbred line populations were genotyped using the newly developed GenoPea 13.2K SNP Array. Individual and consensus genetic maps were built providing insights into the structure and organization of the pea genome. Largely collinear genetic maps of 3918-8503 SNPs were obtained from all mapping populations, and only two of these exhibited putative chromosomal rearrangement signatures. Similar distortion patterns in different populations were noted. A total of 12 802 transcript-derived SNP markers placed on a 15 079-marker high-density, high-resolution consensus map allowed the identification of ohnologue-rich regions within the pea genome and the localization of local duplicates. Dense syntenic networks with sequenced legume genomes were further established, paving the way for the identification of the molecular bases of important agronomic traits segregating in the mapping populations. The information gained on the structure and organization of the genome from this research will undoubtedly contribute to the understanding of the evolution of the pea genome and to its assembly. The GenoPea 13.2K SNP Array and individual and consensus genetic maps are valuable genomic tools for plant scientists to strengthen pea as a model for genetics and physiology and enhance breeding.

Full-length de novo assembly of RNA-seq data in pea (Pisum sativum L.) provides a gene expression atlas and gives insights into root nodulation in this species.

Next-generation sequencing technologies allow an almost exhaustive survey of the transcriptome, even in species with no available genome sequence. To produce a Unigene set representing most of the expressed genes of pea, 20 cDNA libraries produced from various plant tissues harvested at various developmental stages from plants grown under contrasting nitrogen conditions were sequenced. Around one billion reads and 100 Gb of sequence were de novo assembled. Following several steps of redundancy reduction, 46 099 contigs with N50 length of 1667 nt were identified. These constitute the 'Caméor' Unigene set. The high depth of sequencing allowed identification of rare transcripts and detected expression for approximately 80% of contigs in each library. The Unigene set is now available online (http://bios.dijon.inra.fr/FATAL/cgi/pscam.cgi), allowing (i) searches for pea orthologs of candidate genes based on gene sequences from other species, or based on annotation, (ii) determination of transcript expression patterns using various metrics, (iii) identification of uncharacterized genes with interesting patterns of expression, and (iv) comparison of gene ontology pathways between tissues. This resource has allowed identification of the pea orthologs of major nodulation genes characterized in recent years in model species, as a major step towards deciphering unresolved pea nodulation phenotypes. In addition to a remarkable conservation of the early transcriptome nodulation apparatus between pea and Medicago truncatula, some specific features were highlighted. The resource provides a reference for the pea exome, and will facilitate transcriptome and proteome approaches as well as SNP discovery in pea.

Linkage to the mating-type locus across the genus Microbotryum: insights into nonrecombining chromosomes.

Parallels have been drawn between the evolution of nonrecombining regions in fungal mating-type chromosomes and animal and plant sex chromosomes, particularly regarding the stages of recombination cessation forming evolutionary strata of allelic divergence. Currently, evidence and explanations for recombination cessation in fungi are sparse, and the presence of evolutionary strata has been examined in a minimal number of fungal taxa. Here, the basidiomycete genus Microbotryum was used to determine the history of recombination cessation for loci on the mating-type chromosomes. Ancestry of linkage with mating type for 13 loci was assessed across 20 species by a phylogenetic method. No locus was found to exhibit trans-specific polymorphism for alternate alleles as old as the mating pheromone receptor, indicating that ages of linkage to mating type varied among the loci. The ordering of loci in the ancestry of linkage to mating type does not agree with their previously proposed assignments to evolutionary strata. This study suggests that processes capable of influencing divergence between alternate alleles may act at loci in the nonrecombining regions (e.g., gene conversion) and encourages further work to dissect the evolutionary processes acting upon genomic regions that determine mating compatibility.

Genome expansion and gene loss in powdery mildew fungi reveal tradeoffs in extreme parasitism.

Powdery mildews are phytopathogens whose growth and reproduction are entirely dependent on living plant cells. The molecular basis of this life-style, obligate biotrophy, remains unknown. We present the genome analysis of barley powdery mildew, Blumeria graminis f.sp. hordei (Blumeria), as well as a comparison with the analysis of two powdery mildews pathogenic on dicotyledonous plants. These genomes display massive retrotransposon proliferation, genome-size expansion, and gene losses. The missing genes encode enzymes of primary and secondary metabolism, carbohydrate-active enzymes, and transporters, probably reflecting their redundancy in an exclusively biotrophic life-style. Among the 248 candidate effectors of pathogenesis identified in the Blumeria genome, very few (less than 10) define a core set conserved in all three mildews, suggesting that most effectors represent species-specific adaptations.