Comprehensive transcriptome analysis of the highly complex Pisum sativum genome using next generation sequencing

Dirigent isoflavene-forming PsPTS2: 3D structure, stereochemical, and kinetic characterization comparison with pterocarpan-forming PsPTS1 homolog in pea

Pea phytoalexins (-)-maackiain and (+)-pisatin have opposite C6a/C11a configurations, but biosynthetically how this occurs is unknown. Pea dirigent-protein (DP) PsPTS2 generates 7,2'-dihydroxy-4',5'-methylenedioxyisoflav-3-ene (DMDIF), and stereoselectivity toward four possible 7,2'-dihydroxy-4',5'-methylenedioxyisoflavan-4-ol (DMDI) stereoisomers was investigated. Stereoisomer configurations were determined using NMR spectroscopy, electronic circular dichroism, and molecular orbital analyses. PsPTS2 efficiently converted cis-(3R,4R)-DMDI into DMDIF 20-fold faster than the trans-(3R,4S)-isomer. The 4R-configured substrate's near β-axial OH orientation significantly enhanced its leaving group abilities in generating A-ring mono-quinone methide (QM), whereas 4S-isomer's α-equatorial-OH was a poorer leaving group. Docking simulations indicated that the 4R-configured β-axial OH was closest to Asp51, whereas 4S-isomer's α-equatorial OH was further away. Neither cis-(3S,4S)- nor trans-(3S,4R)-DMDIs were substrates, even with the former having C3/C4 stereochemistry as in (+)-pisatin. PsPTS2 used cis-(3R,4R)-7,2'-dihydroxy-4'-methoxyisoflavan-4-ol [cis-(3R,4R)-DMI] and C3/C4 stereoisomers to give 2',7-dihydroxy-4'-methoxyisoflav-3-ene (DMIF). DP homologs may exist in licorice (Glycyrrhiza pallidiflora) and tree legume Bolusanthus speciosus, as DMIF occurs in both species. PsPTS1 utilized cis-(3R,4R)-DMDI to give (-)-maackiain 2200-fold more efficiently than with cis-(3R,4R)-DMI to give (-)-medicarpin. PsPTS1 also slowly converted trans-(3S,4R)-DMDI into (+)-maackiain, reflecting the better 4R configured OH leaving group. PsPTS2 and PsPTS1 provisionally provide the means to enable differing C6a and C11a configurations in (+)-pisatin and (-)-maackiain, via identical DP-engendered mono-QM bound intermediate generation, which PsPTS2 either re-aromatizes to give DMDIF or PsPTS1 intramolecularly cyclizes to afford (-)-maackiain. Substrate docking simulations using PsPTS2 and PsPTS1 indicate cis-(3R,4R)-DMDI binds in the anti-configuration in PsPTS2 to afford DMDIF, and the syn-configuration in PsPTS1 to give maackiain.

Transcriptional and Physiological Analysis Reveal New Insights into the Regulation of Fertilization (N, P, K) on the Growth and Synthesis of Medicinal Components of Dendrobium denneanum

Dendrobium denneanum is an important medicinal and ornamental plant. Its ornamental and medicinal values are affected by its vegetative growth conditions and chemical composition accumulation. This study adopted an orthogonal experimental design to treat D. denneanum with nine different levels of nitrogen (N), potassium (K), and phosphorus (P). The morphological indicators of the plant were positively correlated with the nitrogen concentration. The polysaccharide content was the highest at 1500 mg·L^-1 nitrogen and 3000 mg·L^-1 phosphorous and was 26.84% greater than the control. The flavonoid content increased by 36.2% at 500 mg·L^-1 nitrogen, 2000 mg·L^-1 phosphorous, and 300 mg·L^-1 potassium. Principal component score analysis showed that nitrogen had the most significant impact on the various indicators of D. denneanum, followed by phosphorus and potassium. The comprehensive score showed that the T9 treatment (N: 1500 mg·L^-1, P: 3000 mg·L^-1, K: 500 mg·L^-1) had the strongest effect on D. denneanum. Transcriptional analysis showed that compared with the control, the T9 treatment led to 2277 differentially expressed genes (1230 upregulated and 1047 downregulated). This includes fifteen genes enriched in the MAPK signaling pathway, five genes in phenylpropanoid biosynthesis, and two genes in flavonoid biosynthesis. These genes may be involved in regulating plant growth and the biosynthesis of polysaccharides and flavonoids. This study provides guidance for the optimal use of N, P, and K in the cultivation of D. denneanum.

Mendel: From genes to genome

Two hundred years after the birth of Gregor Mendel, it is an appropriate time to reflect on recent developments in the discipline of genetics, particularly advances relating to the prescient friar's model species, the garden pea (Pisum sativum L.). Mendel's study of seven characteristics established the laws of segregation and independent assortment. The genes underlying four of Mendel's loci (A, LE, I, and R) have been characterized at the molecular level for over a decade. However, the three remaining genes, influencing pod color (GP), pod form (V/P), and the position of flowers (FA/FAS), have remained elusive for a variety of reasons, including a lack of detail regarding the loci with which Mendel worked. Here, we discuss potential candidate genes for these characteristics, in light of recent advances in the genetic resources for pea. These advances, including the pea genome sequence and reverse-genetics techniques, have revitalized pea as an excellent model species for physiological-genetic studies. We also discuss the issues that have been raised with Mendel's results, such as the recent controversy regarding the discrete nature of the characters that Mendel chose and the perceived overly-good fit of his segregations to his hypotheses. We also consider the relevance of these controversies to his lasting contribution. Finally, we discuss the use of Mendel's classical results to teach and enthuse future generations of geneticists, not only regarding the core principles of the discipline, but also its history and the role of hypothesis testing.

Genomics Associated Interventions for Heat Stress Tolerance in Cool Season Adapted Grain Legumes

Cool season grain legumes occupy an important place among the agricultural crops and essentially provide multiple benefits including food supply, nutrition security, soil fertility improvement and revenue for farmers all over the world. However, owing to climate change, the average temperature is steadily rising, which negatively affects crop performance and limits their yield. Terminal heat stress that mainly occurred during grain development phases severely harms grain quality and weight in legumes adapted to the cool season, such as lentils, faba beans, chickpeas, field peas, etc. Although, traditional breeding approaches with advanced screening procedures have been employed to identify heat tolerant legume cultivars. Unfortunately, traditional breeding pipelines alone are no longer enough to meet global demands. Genomics-assisted interventions including new-generation sequencing technologies and genotyping platforms have facilitated the development of high-resolution molecular maps, QTL/gene discovery and marker-assisted introgression, thereby improving the efficiency in legumes breeding to develop stress-resilient varieties. Based on the current scenario, we attempted to review the intervention of genomics to decipher different components of tolerance to heat stress and future possibilities of using newly developed genomics-based interventions in cool season adapted grain legumes.

Identification of heat responsive genes in pea stipules and anthers through transcriptional profiling

Field pea (Pisum sativum L.), a cool-season legume crop, is known for poor heat tolerance. Our previous work identified PR11-2 and PR11-90 as heat tolerant and susceptible lines in a recombinant inbred population. CDC Amarillo, a Canadian elite pea variety, was considered as another heat tolerant variety based on its similar field performance as PR11-2. This study aimed to characterize the differential transcription. Plants of these three varieties were stressed for 3 h at 38°C prior to self-pollination, and RNAs from heat stressed anthers and stipules on the same flowering node were extracted and sequenced via the Illumina NovaSeq platform for the characterization of heat responsive genes. In silico results were further validated by qPCR assay. Differentially expressed genes (DEGs) were identified at log2 |fold change (FC)| ≥ 2 between high temperature and control temperature, the three varieties shared 588 DEGs which were up-regulated and 220 genes which were down-regulated in anthers when subjected to heat treatment. In stipules, 879 DEGs (463/416 upregulation/downregulation) were consistent among varieties. The above heat-induced genes of the two plant organs were related to several biological processes i.e., response to heat, protein folding and DNA templated transcription. Ten gene ontology (GO) terms were over-represented in the consistently down-regulated DEGs of the two organs, and these terms were mainly related to cell wall macromolecule metabolism, lipid transport, lipid localization, and lipid metabolic processes. GO enrichment analysis on distinct DEGs of individual pea varieties suggested that heat affected biological processes were dynamic, and variety distinct responses provide insight into molecular mechanisms of heat-tolerance response. Several biological processes, e.g., cellular response to DNA damage stimulus in stipule, electron transport chain in anther that were only observed in heat induced PR11-2 and CDC Amarillo, and their relevance to field pea heat tolerance is worth further validation.

Development and application of the Faba_bean_130K targeted next-generation sequencing SNP genotyping platform based on transcriptome sequencing

KEY MESSAGE

Large-scale faba bean transcriptome data are available, and the first genotyping platform based on liquid-phase probe targeted capture technology was developed for genetic and molecular breeding studies. Faba bean (Vicia faba L., 2n = 12) is an important food legume crop that is widely grown for multiple uses worldwide. However, no reference genome is currently available due to its very large genome size (approximately 13 Gb) and limited single nucleotide polymorphism (SNP) markers as well as highly efficient genotyping tools have been reported for faba bean. In this study, 16.7 billion clean reads were obtained from transcriptome libraries of flowers and leaves of 102 global faba bean accessions. A total of 243,120 unigenes were de novo assembled and functionally annotated. Moreover, a total of 1,579,411 SNPs were identified and further filtered according to a selection pipeline to develop a high-throughput, flexible, low-cost Faba_bean_130K targeted next-generation sequencing (TNGS) genotyping platform. A set of 69 Chinese faba bean accessions were genotyped with the TNGS genotyping platform, and the average mapping rate of captured reads to reference transcripts was 93.14%, of which 53.23% were located in the targeted regions. The TNGS genotyping results were validated by Sanger sequencing and the average consistency rate reached 93.6%. Comprehensive population genetic analysis was performed on the 69 Chinese faba bean accessions and identified four genetic subgroups correlated with the geographic distribution. This study provides valuable genomic resources and a reliable genotyping tool that could be implemented in genetic and molecular breeding studies to accelerate new cultivar development and improvement in faba bean.

TrancriptomeReconstructoR: data-driven annotation of complex transcriptomes

BACKGROUND

The quality of gene annotation determines the interpretation of results obtained in transcriptomic studies. The growing number of genome sequence information calls for experimental and computational pipelines for de novo transcriptome annotation. Ideally, gene and transcript models should be called from a limited set of key experimental data.

RESULTS

We developed TranscriptomeReconstructoR, an R package which implements a pipeline for automated transcriptome annotation. It relies on integrating features from independent and complementary datasets: (i) full-length RNA-seq for detection of splicing patterns and (ii) high-throughput 5' and 3' tag sequencing data for accurate definition of gene borders. The pipeline can also take a nascent RNA-seq dataset to supplement the called gene model with transient transcripts. We reconstructed de novo the transcriptional landscape of wild type Arabidopsis thaliana seedlings and Saccharomyces cerevisiae cells as a proof-of-principle. A comparison to the existing transcriptome annotations revealed that our gene model is more accurate and comprehensive than the most commonly used community gene models, TAIR10 and Araport11 for A.thaliana and SacCer3 for S.cerevisiae. In particular, we identify multiple transient transcripts missing from the existing annotations. Our new annotations promise to improve the quality of A.thaliana and S.cerevisiae genome research.

CONCLUSIONS

Our proof-of-concept data suggest a cost-efficient strategy for rapid and accurate annotation of complex eukaryotic transcriptomes. We combine the choice of library preparation methods and sequencing platforms with the dedicated computational pipeline implemented in the TranscriptomeReconstructoR package. The pipeline only requires prior knowledge on the reference genomic DNA sequence, but not the transcriptome. The package seamlessly integrates with Bioconductor packages for downstream analysis.

A high coverage reference transcriptome assembly of pea (Pisum sativum L.) mycorrhizal roots

Arbuscular mycorrhiza (AM) is an ancient mutualistic symbiosis formed by 80–90 % of land plant species with the obligatorily biotrophic fungi that belong to the phylum Glomeromycota. This symbiosis is mutually beneficial, as AM fungi feed on plant photosynthesis products, in turn improving the efficiency of nutrient uptake from the environment. The garden pea (Pisum sativum L.), a widely cultivated crop and an important model for genetics, is capable of forming triple symbiotic systems consisting of the plant, AM fungi and nodule bacteria. As transcriptomic and proteomic approaches are being implemented for studying the mutualistic symbioses of pea, a need for a reference transcriptome of genes expressed under these specific conditions for increasing the resolution and the accuracy of other methods arose. Numerous transcriptome assemblies constructed for pea did not include mycorrhizal roots, hence the aim of the study to construct a reference transcriptome assembly of pea mycorrhizal roots. The combined transcriptome of mycorrhizal roots of Pisum sativum cv. Frisson inoculated with Rhizophagus irregularis BEG144 was investigated, and for both the organisms independent transcriptomes were assembled (coverage 177x for pea and 45x for fungus). Genes specific to mycorrhizal roots were found in the assembly, their expression patterns were examined with qPCR on two pea cultivars, Frisson and Finale. The gene expression depended on the inoculation stage and on the pea cultivar. The investigated genes may serve as markers for early stages of inoculation in genetically diverse pea cultivars.

Omics resources and omics-enabled approaches for achieving high productivity and improved quality in pea (Pisum sativum L.)

Pea (Pisum sativum L.), a cool-season legume crop grown in more than 85 countries, is the second most important grain legume and one of the major green vegetables in the world. While pea was historically studied as the genetic model leading to the discovery of the laws of genetics, pea research has lagged behind that of other major legumes in the genomics era, due to its large and complex genome. The evolving climate change and growing population have posed grand challenges to the objective of feeding the world, making it essential to invest research efforts to develop multi-omics resources and advanced breeding tools to support fast and continuous development of improved pea varieties. Recently, the pea researchers have achieved key milestones in omics and molecular breeding. The present review provides an overview of the recent important progress including the development of genetic resource databases, high-throughput genotyping assays, reference genome, genes/QTLs responsible for important traits, transcriptomic, proteomic, and phenomic atlases of various tissues under different conditions. These multi-faceted resources have enabled the successful implementation of various markers for monitoring early-generation populations as in marker-assisted backcrossing breeding programs. The emerging new breeding approaches such as CRISPR, speed breeding, and genomic selection are starting to change the paradigm of pea breeding. Collectively, the rich omics resources and omics-enable breeding approaches will enhance genetic gain in pea breeding and accelerate the release of novel pea varieties to meet the elevating demands on productivity and quality.

Transcriptome analysis of blood for the discovery of sex-related genes in ricefield eel Monopterus albus

The blood acts as a transfer channel for a variety of factors in the whole body. The ricefield eel (Monopterus albus) is a protogynous hermaphrodite vertebrate. Until now, no research has reported an analysis of the blood transcriptome during the process of sexual development in the ricefield eel. In this study, the transcriptome sequencing of blood samples from male and female ricefield eels was completed with a total of 34.70 Gb clean data. The clean data of each sample all reached 5.23 GB, and the percent of the Q30 basic group was 88.62% and above. A total of 106,369 unigenes were obtained after assembly, including 13,296 unigenes with a length of more than 1 kb. Further functional annotation analysis showed that there are 28,522 unigenes that can be annotated. The annotations of genes with differential expression revealed that there were 563 genes with significant differential expression in the blood of male and female ricefield eels, including 91 upregulated genes and 472 downregulated genes. Among which, 14 genes may be closely related to sex differentiation, the qPCR was used to confirmed the expression pattern of those genes and result shown that 11 genes were downregulated and 3 genes were upregulated, consistent with the results of our RNA-Seq analysis. This blood transcript dataset will open future research avenues on ricefield eel sex development and differentiation.

Plant science's next top models

BACKGROUND

Model organisms are at the core of life science research. Notable examples include the mouse as a model for humans, baker's yeast for eukaryotic unicellular life and simple genetics, or the enterobacteria phage λ in virology. Plant research was an exception to this rule, with researchers relying on a variety of non-model plants until the eventual adoption of Arabidopsis thaliana as primary plant model in the 1980s. This proved to be an unprecedented success, and several secondary plant models have since been established. Currently, we are experiencing another wave of expansion in the set of plant models.

SCOPE

Since the 2000s, new model plants have been established to study numerous aspects of plant biology, such as the evolution of land plants, grasses, invasive and parasitic plant life, adaptation to environmental challenges, and the development of morphological diversity. Concurrent with the establishment of new plant models, the advent of the 'omics' era in biology has led to a resurgence of the more complex non-model plants. With this review, we introduce some of the new and fascinating plant models, outline why they are interesting subjects to study, the questions they will help to answer, and the molecular tools that have been established and are available to researchers.

CONCLUSIONS

Understanding the molecular mechanisms underlying all aspects of plant biology can only be achieved with the adoption of a comprehensive set of models, each of which allows the assessment of at least one aspect of plant life. The model plants described here represent a step forward towards our goal to explore and comprehend the diversity of plant form and function. Still, several questions remain unanswered, but the constant development of novel technologies in molecular biology and bioinformatics is already paving the way for the next generation of plant models.

Comprehensive Stress-Based De Novo Transcriptome Assembly and Annotation of Guar (Cyamopsis tetragonoloba (L.) Taub.): An Important Industrial and Forage Crop

The forage crop Guar (Cyamopsis tetragonoloba (L.) Taub.) has the ability to endure heat, drought, and mild salinity. A complete image on its genic architecture will promote our understanding about gene expression networks and different tolerance mechanisms at the molecular level. Therefore, whole mRNA sequence approach on the Guar plant was conducted to provide a snapshot of the mRNA information in the cell under salinity, heat, and drought stresses to be integrated with previous transcriptomic studies. RNA-Seq technology was employed to perform a 2 × 100 paired-end sequencing using an Illumina HiSeq 2500 platform for the transcriptome of leaves of C. tetragonoloba under normal, heat, drought, and salinity conditions. Trinity was used to achieve a de novo assembly followed by gene annotation, functional classification, metabolic pathway analysis, and identification of SSR markers. A total of 218.2 million paired-end raw reads (~44 Gbp) were generated. Of those, 193.5M paired-end reads of high quality were used to reconstruct a total of 161,058 transcripts (~266 Mbp) with N50 of 2552 bp and 61,508 putative genes. There were 6463 proteins having >90% full-length coverage against the Swiss-Prot database and 94% complete orthologs against Embryophyta. Approximately, 62.87% of transcripts were blasted, 50.46% mapped, and 43.50% annotated. A total of 4715 InterProScan families, 3441 domains, 74 repeats, and 490 sites were detected. Biological processes, molecular functions, and cellular components comprised 64.12%, 25.42%, and 10.4%, respectively. The transcriptome was associated with 985 enzymes and 156 KEGG pathways. A total of 27,066 SSRs were gained with an average frequency of one SSR/9.825 kb in the assembled transcripts. This resulting data will be helpful for the advanced analysis of Guar to multi-stress tolerance.

Development of PCR markers specific to Dasypyrum villosum genome based on transcriptome data and their application in breeding Triticum aestivum-D. villosum#4 alien chromosome lines

Background

Dasypyrum villosum is an important wild species of wheat (Triticum aestivum L.) and harbors many desirable genes that can be used to improve various traits of wheat. Compared with other D. villosum accessions, D. villosum#4 still remains less studied. In particular, chromosomes of D. villosum#4 except 6V#4 have not been introduced into wheat by addition or substitution and translocation, which is an essential step to identify and apply the alien desired genes. RNA-seq technology can generate large amounts of transcriptome sequences and accelerate the development of chromosome-specific molecular markers and assisted selection of alien chromosome line.

Results

We obtained the transcriptome of D. villosum#4 via a high-throughput sequencing technique, and then developed 76 markers specific to each chromosome arm of D. villosum#4 based on the bioinformatic analysis of the transcriptome data. The D. villosum#4 sequences containing the specific DNA markers were expected to be involved in different genes, among which most had functions in metabolic processes. Consequently, we mapped these newly developed molecular markers to the homologous chromosome of barley and obtained the chromosome localization of these markers on barley genome. Then we analyzed the collinearity of these markers among D. villosum, wheat, and barley. In succession, we identified six types of T. aestivum-D. villosum#4 alien chromosome lines which had one or more than one D. villosum#4 chromosome in the cross and backcross BC₃F₅ populations between T. durum–D. villosum#4 amphidiploid TH3 and wheat cv. Wan7107 by employing the selected specific markers, some of which were further confirmed to be translocation or addition lines by genomic in situ hybridization (GISH).

Conclusion

Seventy-six PCR markers specific to chromosomes of D. villosum#4 based on transcriptome data were developed in the current study and their collinearity among D. villosum, wheat, and barley were carried out. Six types of Triticum aestivum-D. villosum#4 alien chromosome lines were identified by using 12 developed markers and some of which were further confirmed by GISH. These novel T. aestivum-D. villosum#4 chromosome lines have great potential to be used for the introduction of desirable genes from D. villosum#4 into wheat by chromosomal translocation to breed new wheat varieties.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5630-4) contains supplementary material, which is available to authorized users.

TIR1 auxin receptors are implicated in the differential response to 4-Cl-IAA and IAA in developing pea fruit

The auxins indole-3-acetic acid (IAA) and 4-chloroindole-3-acetic acid (4-Cl-IAA) occur naturally in pea (Pisum sativum); however, only 4-Cl-IAA mimics the presence of seeds in stimulating pericarp growth. To examine if this differential auxin effect is mediated through TIR1/AFB auxin receptors, pea TIR1 and AFB2 homologs were functionally characterized in Arabidopsis, and receptor expression, and auxin distribution and action were profiled in developing pea fruits. PsTIR1a, PsTIR1b, and PsAFB2 restored the auxin-sensitive root growth response to the mutant Arabidopsis seedlings Attir1-10 and/or Attir1-10 afb2-3. Expression of PsTIR1 or AtTIR1 in Attir1-10 afb2-3 mutants also restored the greater root inhibitory response of 4-Cl-IAA compared to that of IAA, implicating TIR1 receptors in this response. The ability of 4-Cl-IAA to stimulate a stronger DR5::GUS auxin response than IAA at the same concentration in pea pericarps was associated with its ability to enrich the auxin-receptor transcript pool with PsTIR1a and PsAFB2 by decreasing the transcript abundance of PsTIR1b (mimicking results in pericarps with developing seeds). Therefore, the markedly different effect of IAA and 4-Cl-IAA on pea fruit growth may at least partially involve TIR1/AFB receptors and the differential modulation of their population, resulting in specific Aux/IAA protein degradation that leads to an auxin-specific tissue response.

Changes in gene expression during germination reveal pea genotypes with either "quiescence" or "escape" mechanisms of waterlogging tolerance

Waterlogging causes germination failure in pea (Pisum sativum L.). Three genotypes (BARI Motorshuti-3, Natore local-2 [NL-2], and Kaspa) contrasting in ability to germinate in waterlogged soil were exposed to different durations of waterlogging. Whole genome RNAseq was employed to capture differentially expressing genes. The ability to germinate in waterlogged soil was associated with testa colour and testa membrane integrity as confirmed by electrical conductivity measurements. Genotypes Kaspa and NL-2 displayed different mechanisms of tolerance. In Kaspa, an energy conserving strategy was indicated by a strong upregulation of tyrosine protein kinsase and down regulation of linoleate 9S-lipoxygenase 5, a fat metabolism gene. In contrast, a faster energy utilization strategy was suggested in NL-2 by the marked upregulation of a subtilase family protein and peroxisomal adenine nucleotide carrier 2, a fat metabolizing gene. Waterlogging susceptibility in germinating seeds of genotype BARI Motorshuti-3 was linked to upregulation of a kunitz-type trypsin/protease inhibitor that blocks protein metabolism and may lead to excessive lipid metabolism and the membrane leakage associated with waterlogging damage. Pathway analyses based on gene ontologies showed seed storage protein metabolism as upregulated in tolerant genotypes and downregulated in the sensitive genotype. Understanding the tolerance mechanism provides a platform to breed for adaptation to waterlogging stress at germination in pea.

Identification of Stipules reduced, a leaf morphology gene in pea (Pisum sativum)

Pea (Pisum sativum) is one of relatively few genetically amenable plant species with compound leaves. Pea leaves have a variety of specialized organs: leaflets, tendrils, pulvini and stipules, which enable the identification of mutations that transform or affect distinct parts of the leaf. Characterization of these mutations offers insights into the development and evolution of novel leaf traits. The previously characterized morphological gene Cochleata, conferring stipule identity, was known to interact with Stipules reduced (St), which conditions stipule size in pea, but the St gene remained unknown. Here we analysed Fast Neutron irradiated pea mutants by restriction site associated DNA sequencing. We identified St as a gene encoding a C2H2 zinc finger transcription factor that is regulated by Cochleata. St regulates both cell division and cell expansion in the stipule. Our approach shows how systematic genome-wide screens can be used successfully for the analysis of traits in species for which whole genome sequences are not available.

Pea Marker Database (PMD) - A new online database combining known pea (Pisum sativum L.) gene-based markers

Pea (Pisum sativum L.) is the oldest model object of plant genetics and one of the most agriculturally important legumes in the world. Since the pea genome has not been sequenced yet, identification of genes responsible for mutant phenotypes or desirable agricultural traits is usually performed via genetic mapping followed by candidate gene search. Such mapping is best carried out using gene-based molecular markers, as it opens the possibility for exploiting genome synteny between pea and its close relative Medicago truncatula Gaertn., possessing sequenced and annotated genome. In the last 5 years, a large number of pea gene-based molecular markers have been designed and mapped owing to the rapid evolution of "next-generation sequencing" technologies. However, the access to the complete set of markers designed worldwide is limited because the data are not uniformed and therefore hard to use. The Pea Marker Database was designed to combine the information about pea markers in a form of user-friendly and practical online tool. Version 1 (PMD1) comprises information about 2484 genic markers, including their locations in linkage groups, the sequences of corresponding pea transcripts and the names of related genes in M. truncatula. Version 2 (PMD2) is an updated version comprising 15944 pea markers in the same format with several advanced features. To test the performance of the PMD, fine mapping of pea symbiotic genes Sym13 and Sym27 in linkage groups VII and V, respectively, was carried out. The results of mapping allowed us to propose the Sen1 gene (a homologue of SEN1 gene of Lotus japonicus (Regel) K. Larsen) as the best candidate gene for Sym13, and to narrow the list of possible candidate genes for Sym27 to ten, thus proving PMD to be useful for pea gene mapping and cloning. All information contained in PMD1 and PMD2 is available at www.peamarker.arriam.ru.

De novo transcriptome assembly reveals high transcriptional complexity in Pisum sativum axillary buds and shows rapid changes in expression of diurnally regulated genes

BACKGROUND

The decision for a bud to grow into a branch is a key regulatory process affecting plant architecture. In order to study molecular processes regulating axillary bud outgrowth in the model plant garden pea (Pisum sativum), we sequenced the axillary bud transcriptome and performed de novo transcriptome assembly.

RESULTS

We assembled a pea axillary bud transcriptome into 81,774 transcripts comprised of 194,067 isoforms. This new pea transcriptome resource is both comprehensive and representative, as shown by comparison to other available pea sequence resources. Over half of the transcriptome could be annotated based on sequence homology to Arabidopsis thaliana proteins, while almost one quarter of the isoforms were identified as putative long non-coding RNAs (lncRNAs). This transcriptome will be useful in studies of pea buds because it includes genes expressed specifically in buds which are not represented in other transcriptome studies. We also investigated the impact of a short time collection series on gene expression. Differential gene expression analysis identified 142 transcripts changing within the short 170 min time frame that the buds were harvested within. Thirty-three of these transcripts are implicated in diurnal fluctuations in other flowering plants, while the remaining transcripts include 31 putative lncRNA. Further investigation of the differentially expressed transcripts found an enrichment of genes involved in post-transcriptional regulation, including RNA processing and modification, as well as genes involved in fatty acid biosynthesis and oxidative phosphorylation.

CONCLUSIONS

We have sequenced and assembled a high quality pea bud transcriptome containing both coding and non-coding RNA transcripts that will be useful for further studies into axillary bud outgrowth. Over the short sample collection time frame of just 170 min, we identified differentially expressed coding and non-coding RNA, some of which are implicated in diurnal regulation, highlighting the utility of our transcriptome resource in identifying gene expression changes and informing future experimental designs.

Proteomics offers insight to the mechanism behind Pisum sativum L. response to pea seed-borne mosaic virus (PSbMV)

Pea seed-borne mosaic virus (PSbMV) significantly reduces yields in a broad spectra of legumes. The eukaryotic translation initiation factor has been shown to confer resistance to this pathogen, thus implying that translation and proteome dynamics play a role in resistance. This study presents the results of a proteome-wide analysis of Pisum sativum L. response to PSbMV infection. LC-MS profiling of two contrasting pea cultivars, resistant (B99) and susceptible (Raman) to PSbMV infection, detected >2300 proteins, 116 of which responded to PSbMV ten and/or twenty days post-inoculation. These differentially abundant proteins are involved in number of processes that have previously been reported in the plant-pathogen response, including protein and amino acid metabolism, stress signaling, redox homeostasis, carbohydrate metabolism, and lipid metabolism. We complemented our proteome-wide analysis work with targeted analyses of free amino acids and selected small molecules, fatty acid profiling, and enzyme activity assays. Data from these additional experiments support our findings and validate the biological relevance of the observed proteome changes. We found surprising similarities in the resistant and susceptible cultivars, which implies that a seemingly unaffected plant, with no detectable levels of PSbMV, actively suppresses viral replication.

BIOLOGICAL SIGNIFICANCE

Plant resistance to PSbMV is connected to translation initiation factors, yet the processes involved are still poorly understood at the proteome level. To the best of our knowledge, this is the first survey of the global proteomic response to PSbMV in plants. The combination of label-free LC-MS profiling and two contrasting cultivars (resistant and susceptible) provided highly sensitive snapshots of protein abundance in response to PSbMV infection. PSbMV is a member of the largest family of plant viruses and our results are in accordance with previously characterized potyvirus-responsive proteomes. Hence, the results of this study can further extend our knowledge about these pathogens. We also show that even though no viral replication is detected in the PSbMV-resistant cultivar B99, it is still significantly affected by PSbMV inoculation.

A Combined Comparative Transcriptomic, Metabolomic, and Anatomical Analyses of Two Key Domestication Traits: Pod Dehiscence and Seed Dormancy in Pea (Pisum sp.)

The origin of the agriculture was one of the turning points in human history, and a central part of this was the evolution of new plant forms, domesticated crops. Seed dispersal and germination are two key traits which have been selected to facilitate cultivation and harvesting of crops. The objective of this study was to analyze anatomical structure of seed coat and pod, identify metabolic compounds associated with water-impermeable seed coat and differentially expressed genes involved in pea seed dormancy and pod dehiscence. Comparative anatomical, metabolomics, and transcriptomic analyses were carried out on wild dormant, dehiscent Pisum elatius (JI64, VIR320) and cultivated, indehiscent Pisum sativum non-dormant (JI92, Cameor) and recombinant inbred lines (RILs). Considerable differences were found in texture of testa surface, length of macrosclereids, and seed coat thickness. Histochemical and biochemical analyses indicated genotype related variation in composition and heterogeneity of seed coat cell walls within macrosclereids. Liquid chromatography-electrospray ionization/mass spectrometry and Laser desorption/ionization-mass spectrometry of separated seed coats revealed significantly higher contents of proanthocyanidins (dimer and trimer of gallocatechin), quercetin, and myricetin rhamnosides and hydroxylated fatty acids in dormant compared to non-dormant genotypes. Bulk Segregant Analysis coupled to high throughput RNA sequencing resulted in identification of 770 and 148 differentially expressed genes between dormant and non-dormant seeds or dehiscent and indehiscent pods, respectively. The expression of 14 selected dormancy-related genes was studied by qRT-PCR. Of these, expression pattern of four genes: porin (MACE-S082), peroxisomal membrane PEX14-like protein (MACE-S108), 4-coumarate CoA ligase (MACE-S131), and UDP-glucosyl transferase (MACE-S139) was in agreement in all four genotypes with Massive analysis of cDNA Ends (MACE) data. In case of pod dehiscence, the analysis of two candidate genes (SHATTERING and SHATTERPROOF) and three out of 20 MACE identified genes (MACE-P004, MACE-P013, MACE-P015) showed down-expression in dorsal and ventral pod suture of indehiscent genotypes. Moreover, MACE-P015, the homolog of peptidoglycan-binding domain or proline-rich extensin-like protein mapped correctly to predicted Dpo1 locus on PsLGIII. This integrated analysis of the seed coat in wild and cultivated pea provides new insight as well as raises new questions associated with domestication and seed dormancy and pod dehiscence.

Generation and Characterisation of a Reference Transcriptome for Lentil (Lens culinaris Medik.)

RNA-Seq using second-generation sequencing technologies permits generation of a reference unigene set for a given species, in the absence of a well-annotated genome sequence, supporting functional genomics studies, gene characterisation and detailed expression analysis for specific morphophysiological or environmental stress response traits. A reference unigene set for lentil has been developed, consisting of 58,986 contigs and scaffolds with an N50 length of 1719 bp. Comparison to gene complements from related species, reference protein databases, previously published lentil transcriptomes and a draft genome sequence validated the current dataset in terms of degree of completeness and utility. A large proportion (98%) of unigenes were expressed in more than one tissue, at varying levels. Candidate genes associated with mechanisms of tolerance to both boron toxicity and time of flowering were identified, which can eventually be used for the development of gene-based markers. This study has provided a comprehensive, assembled and annotated reference gene set for lentil that can be used for multiple applications, permitting identification of genes for pathway-specific expression analysis, genetic modification approaches, development of resources for genotypic analysis, and assistance in the annotation of a future lentil genome sequence.

Reconstruction of a composite comparative map composed of ten legume genomes

The Fabaceae (legume family) is the third largest and the second of agricultural importance among flowering plant groups. In this study, we report the reconstruction of a composite comparative map composed of ten legume genomes, including seven species from the galegoid clade (Medicago truncatula, Medicago sativa, Lens culinaris, Pisum sativum, Lotus japonicus, Cicer arietinum, Vicia faba) and three species from the phaseoloid clade (Vigna radiata, Phaseolus vulgaris, Glycine max). To accomplish this comparison, a total of 209 cross-species gene-derived markers were employed. The comparative analysis resulted in a single extensive genetic/genomic network composed of 93 chromosomes or linkage groups, from which 110 synteny blocks and other evolutionary events (e.g., 13 inversions) were identified. This comparative map also allowed us to deduce several large scale evolutionary events, such as chromosome fusion/fission, with which might explain differences in chromosome numbers among compared species or between the two clades. As a result, useful properties of cross-species genic markers were re-verified as an efficient tool for cross-species translation of genomic information, and similar approaches, combined with a high throughput bioinformatic marker design program, should be effective for applying the knowledge of trait-associated genes to other important crop species for breeding purposes. Here, we provide a basic comparative framework for the ten legume species, and expect to be usefully applied towards the crop improvement in legume breeding.

Identification of LATE BLOOMER2 as a CYCLING DOF FACTOR Homolog Reveals Conserved and Divergent Features of the Flowering Response to Photoperiod in Pea

The molecular pathways responsible for the flowering response to photoperiod have been extensively studied in Arabidopsis thaliana and cereals but remain poorly understood in other major plant groups. Here, we describe a dominant mutant at the LATE BLOOMER2 (LATE2) locus in pea (Pisum sativum) that is late-flowering with a reduced response to photoperiod. LATE2 acts downstream of light signaling and the circadian clock to control expression of the main photoperiod-regulated FT gene, FTb2, implying that it plays a primary role in photoperiod measurement. Mapping identified the CYCLING DOF FACTOR gene CDFc1 as a strong candidate for LATE2, and the late2-1D mutant was found to carry a missense mutation in CDFc1 that impairs its capacity to bind to the blue-light photoreceptor FKF1 in yeast two-hybrid assays and delays flowering in Arabidopsis when overexpressed. Arabidopsis CDF genes are important negative regulators of CONSTANS (CO) transcription, but we found no effect of LATE2 on the transcription of pea CO-LIKE genes, nor on genes in any other families previously implicated in the activation of FT in Arabidopsis. Our results reveal an important component of the pea photoperiod response pathway and support the view that regulation of FTb2 expression by photoperiod occurs via a CO-independent mechanism.

Assessing zygosity in progeny of transgenic plants: current methods and perspectives

Homozygosity is highly desirable in transgenic plants research to ensure the stable integration and inheritance of transgene(s). Simple, reliable and high-throughput techniques to detect the zygosity of transgenic events in plants are invaluable tools for biotechnology and plant breeding companies. Currently, a number of basic techniques are being used to determine the zygosity of transgenic plants in T1 generation. For successful application of any technique, precision and simplicity of approach combined with the power of resolution are important parameters. On the basis of simplicity, resolution and cost involved, the available techniques have been classified into three major classes which are conventional methods, current methods and next generation methods. Conventional methods include antibiotic marker-based selection and the highly labor intensive Southern blot analysis. In contrast, methods such as real time PCR, TAIL PCR and competitive PCR are not only cost effective but rapid as well. Moreover, methods such as NGS, digital PCR and loop-mediated isothermal amplification also provide a cost effective, fast and not so labor intensive substitute of current methods. In this review, we have attempted to compare and contrast all the available efficient methods to distinguish homozygous plants in progeny of transgenics. This review also provides information of various techniques available for determining zygosity in plants so as to permit researchers to make informed choices of techniques that best suit their analyses. More importantly, detection and subsequent selection of homozygous individuals is central for facilitating the movement of transgenic plants from the laboratory to the field.

OEP40, a Regulated Glucose-permeable β-Barrel Solute Channel in the Chloroplast Outer Envelope Membrane

Chloroplasts and mitochondria are unique endosymbiotic cellular organelles surrounded by two membranes. Essential metabolic networking between these compartments and their hosting cells requires the exchange of a large number of biochemical pathway intermediates in a directed and coordinated fashion across their inner and outer envelope membranes. Here, we describe the identification and functional characterization of a highly specific, regulated solute channel in the outer envelope of chloroplasts, named OEP40. Loss of OEP40 function in Arabidopsis thaliana results in early flowering under cold temperature. The reconstituted recombinant OEP40 protein forms a high conductance β-barrel ion channel with subconductant states in planar lipid bilayers. The OEP40 channel is slightly cation-selective PK+/PCl- ≈ 4:1 and rectifying (i⃗/i⃖ ≅ 2) with a slope conductance of Ḡmax ≅ 690 picosiemens. The OEP40 channel has a restriction zone diameter of ≅1.4 nm and is permeable for glucose, glucose 1-phosphate and glucose 6-phosphate, but not for maltose. Moreover, channel properties are regulated by trehalose 6-phosphate, which cannot permeate. Altogether, our results indicate that OEP40 is a "glucose-gate" in the outer envelope membrane of chloroplasts, facilitating selective metabolite exchange between chloroplasts and the surrounding cell.

SNP discovery and genetic mapping using genotyping by sequencing of whole genome genomic DNA from a pea RIL population

BACKGROUND

Progress in genetics and breeding in pea still suffers from the limited availability of molecular resources. SNP markers that can be identified through affordable sequencing processes, without the need for prior genome reduction or a reference genome to assemble sequencing data would allow the discovery and genetic mapping of thousands of molecular markers. Such an approach could significantly speed up genetic studies and marker assisted breeding for non-model species.

RESULTS

A total of 419,024 SNPs were discovered using HiSeq whole genome sequencing of four pea lines, followed by direct identification of SNP markers without assembly using the discoSnp tool. Subsequent filtering led to the identification of 131,850 highly designable SNPs, polymorphic between at least two of the four pea lines. A subset of 64,754 SNPs was called and genotyped by short read sequencing on a subpopulation of 48 RILs from the cross 'Baccara' x 'PI180693'. This data was used to construct a WGGBS-derived pea genetic map comprising 64,263 markers. This map is collinear with previous pea consensus maps and therefore with the Medicago truncatula genome. Sequencing of four additional pea lines showed that 33 % to 64 % of the mapped SNPs, depending on the pairs of lines considered, are polymorphic and can therefore be useful in other crosses. The subsequent genotyping of a subset of 1000 SNPs, chosen for their mapping positions using a KASP™ assay, showed that almost all generated SNPs are highly designable and that most (95 %) deliver highly qualitative genotyping results. Using rather low sequencing coverages in SNP discovery and in SNP inferring did not hinder the identification of hundreds of thousands of high quality SNPs.

CONCLUSIONS

The development and optimization of appropriate tools in SNP discovery and genetic mapping have allowed us to make available a massive new genomic resource in pea. It will be useful for both fine mapping within chosen QTL confidence intervals and marker assisted breeding for important traits in pea improvement.

Reversible Burst of Transcriptional Changes during Induction of Crassulacean Acid Metabolism in Talinum triangulare

Drought tolerance is a key factor for agriculture in the 21st century as it is a major determinant of plant survival in natural ecosystems as well as crop productivity. Plants have evolved a range of mechanisms to cope with drought, including a specialized type of photosynthesis termed Crassulacean acid metabolism (CAM). CAM is associated with stomatal closure during the day as atmospheric CO2 is assimilated primarily during the night, thus reducing transpirational water loss. The tropical herbaceous perennial species Talinum triangulare is capable of transitioning, in a facultative, reversible manner, from C3 photosynthesis to weakly expressed CAM in response to drought stress. The transcriptional regulation of this transition has been studied. Combining mRNA-Seq with targeted metabolite measurements, we found highly elevated levels of CAM-cycle enzyme transcripts and their metabolic products in T. triangulare leaves upon water deprivation. The carbohydrate metabolism is rewired to reduce the use of reserves for growth to support the CAM-cycle and the synthesis of compatible solutes. This large-scale expression dataset of drought-induced CAM demonstrates transcriptional regulation of the C3-CAM transition. We identified candidate transcription factors to mediate this photosynthetic plasticity, which may contribute in the future to the design of more drought-tolerant crops via engineered CAM.

KNOTTED1-LIKE HOMEOBOX 3: a new regulator of symbiotic nodule development

KNOX transcription factors (TFs) regulate different aspects of plant development essentially through their effects on phytohormone metabolism. In particular, KNOX TF SHOOTMERISTEMLESS activates the cytokinin biosynthesis ISOPENTENYL TRANSFERASE (IPT) genes in the shoot apical meristem. However, the role of KNOX TFs in symbiotic nodule development and their possible effects on phytohormone metabolism during nodulation have not been studied to date. Cytokinin is a well-known regulator of nodule development, playing the key role in the regulation of cell division during nodule primordium formation. Recently, the activation of IPT genes was shown to take place during nodulation. Therefore, it was hypothesized that KNOX TFs may regulate nodule development and activate cytokinin biosynthesis upon nodulation. This study analysed the expression of different KNOX genes in Medicago truncatula Gaertn. and Pisum sativum L. Among them, the KNOX3 gene was upregulated in response to rhizobial inoculation in both species. pKNOX3::GUS activity was observed in developing nodule primordium. KNOX3 ectopic expression caused the formation of nodule-like structures on transgenic root without bacterial inoculation, a phenotype similar to one described previously for legumes with constitutive activation of the cytokinin receptor. Furthermore, in transgenic roots with MtKNOX3 knockdown, downregulation of A-type cytokinin response genes was found, as well as the MtIPT3 and LONELYGUY2 (MtLOG2) gene being involved in cytokinin activation. Taken together, these findings suggest that KNOX3 gene is involved in symbiotic nodule development and may regulate cytokinin biosynthesis/activation upon nodule development in legume plants.

Comparative Transcriptome Analysis between the Fungal Plant Pathogens Sclerotinia sclerotiorum and S. trifoliorum Using RNA Sequencing

The fungal plant pathogens Sclerotinia sclerotiorum and S. trifoliorum are morphologically similar, but differ considerably in host range. In an effort to elucidate mechanisms of the host range difference, transcriptomes of the 2 species at vegetative growth stage were compared to gain further insight into commonality and uniqueness in gene expression and pathogenic mechanisms of the 2 closely related pathogens. A total of 23133 and 21043 unique transcripts were obtained from S. sclerotiorum and S. trifoliorum, respectively. Approximately 43% of the transcripts were genes with known functions for both species. Among 1411 orthologous contigs, about 10% (147) were more highly (>3-fold) expressed in S. trifoliorum than in S. sclerotiorum, and about 12% (173) of the orthologs were more highly (>3-fold) expressed in S. sclerotiorum than in S. trifoliorum. The expression levels of genes on the supercontig 30 have the highest correlation coefficient value between the 2 species. Twenty-seven contigs were found to be new and unique for S. trifoliorum. Additionally, differences in expressed genes involved in pathogenesis like oxalate biosynthesis and endopolygalacturonases were detected between the 2 species. The analyses of the transcriptomes not only discovered similarities and uniqueness in gene expression between the 2 closely related species, providing additional information for annotation the S. sclerotiorum genome, but also provided foundation for comparing the transcriptomes with host-infecting transcriptomes.

Genomic Tools in Pea Breeding Programs: Status and Perspectives

Pea (Pisum sativum L.) is an annual cool-season legume and one of the oldest domesticated crops. Dry pea seeds contain 22-25% protein, complex starch and fiber constituents, and a rich array of vitamins, minerals, and phytochemicals which make them a valuable source for human consumption and livestock feed. Dry pea ranks third to common bean and chickpea as the most widely grown pulse in the world with more than 11 million tons produced in 2013. Pea breeding has achieved great success since the time of Mendel's experiments in the mid-1800s. However, several traits still require significant improvement for better yield stability in a larger growing area. Key breeding objectives in pea include improving biotic and abiotic stress resistance and enhancing yield components and seed quality. Taking advantage of the diversity present in the pea genepool, many mapping populations have been constructed in the last decades and efforts have been deployed to identify loci involved in the control of target traits and further introgress them into elite breeding materials. Pea now benefits from next-generation sequencing and high-throughput genotyping technologies that are paving the way for genome-wide association studies and genomic selection approaches. This review covers the significant development and deployment of genomic tools for pea breeding in recent years. Future prospects are discussed especially in light of current progress toward deciphering the pea genome.

De Novo Assembly of the Pea (Pisum sativum L.) Nodule Transcriptome

The large size and complexity of the garden pea (Pisum sativum L.) genome hamper its sequencing and the discovery of pea gene resources. Although transcriptome sequencing provides extensive information about expressed genes, some tissue-specific transcripts can only be identified from particular organs under appropriate conditions. In this study, we performed RNA sequencing of polyadenylated transcripts from young pea nodules and root tips on an Illumina GAIIx system, followed by de novo transcriptome assembly using the Trinity program. We obtained more than 58,000 and 37,000 contigs from "Nodules" and "Root Tips" assemblies, respectively. The quality of the assemblies was assessed by comparison with pea expressed sequence tags and transcriptome sequencing project data available from NCBI website. The "Nodules" assembly was compared with the "Root Tips" assembly and with pea transcriptome sequencing data from projects indicating tissue specificity. As a result, approximately 13,000 nodule-specific contigs were found and annotated by alignment to known plant protein-coding sequences and by Gene Ontology searching. Of these, 581 sequences were found to possess full CDSs and could thus be considered as novel nodule-specific transcripts of pea. The information about pea nodule-specific gene sequences can be applied for gene-based markers creation, polymorphism studies, and real-time PCR.

De novo transcriptome analysis of Medicago falcata reveals novel insights about the mechanisms underlying abiotic stress-responsive pathway

Background

The entire world is facing a deteriorating environment. Understanding the mechanisms underlying plant responses to external abiotic stresses is important for breeding stress-tolerant crops and herbages. Phytohormones play critical regulatory roles in plants in the response to external and internal cues to regulate growth and development. Medicago falcata is one of the stress-tolerant candidate leguminous species and is able to fix atmospheric nitrogen. This ability allows leguminous plants to grow in nitrogen deficient soils.

Methods

We performed Illumina sequencing of cDNA prepared from abiotic stress treated M. falcata. Sequencedreads were assembled to provide a transcriptome resource. Transcripts were annotated using BLASTsearches against the NCBI non-redundant database and gene ontology definitions were assigned. Acomparison among the three abiotic stress treated samples was carried out. The expression of transcriptswas confirmed with qRT-PCR.

Results

We present an abiotic stress-responsive M. falcata transcriptome using next-generation sequencing data from samples grown under standard, dehydration, high salinity, and cold conditions. We combined reads from all samples and de novo assembled 98,515 transcripts to build the M. falcata gene index. A comprehensive analysis of the transcriptome revealed abiotic stress-responsive mechanisms underlying the metabolism and core signalling components of major phytohormones. We identified nod factor signalling pathways during early symbiotic nodulation that are modified by abiotic stresses. Additionally, a global comparison of homology between the M. falcata and M. truncatula transcriptomes, along with five other leguminous species, revealed a high level of global sequence conservation within the family.

Conclusions

M. falcata is shown to be a model candidate for studying abiotic stress-responsive mechanisms in legumes. This global gene expression analysis provides new insights into the biochemical and molecular mechanisms involved in the acclimation to abiotic stresses. Our data provides many gene candidates that might be used for herbage and crop breeding. Additionally, FalcataBase (http://bioinformatics.cau.edu.cn/falcata/) was built for storing these data.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2019-x) contains supplementary material, which is available to authorized users.

Pre-fractionation strategies to resolve pea (Pisum sativum) sub-proteomes

Legumes are important crop plants and pea (Pisum sativum L.) has been investigated as a model with respect to several physiological aspects. The sequencing of the pea genome has not been completed. Therefore, proteomic approaches are currently limited. Nevertheless, the increasing numbers of available EST-databases as well as the high homology of the pea and medicago genome (Medicago truncatula Gaertner) allow the successful identification of proteins. Due to the un-sequenced pea genome, pre-fractionation approaches have been used in pea proteomic surveys in the past. Aside from a number of selective proteome studies on crude extracts and the chloroplast, few studies have targeted other components such as the pea secretome, an important sub-proteome of interest due to its role in abiotic and biotic stress processes. The secretome itself can be further divided into different sub-proteomes (plasma membrane, apoplast, cell wall proteins). Cell fractionation in combination with different gel-electrophoresis, chromatography methods and protein identification by mass spectrometry are important partners to gain insight into pea sub-proteomes, post-translational modifications and protein functions. Overall, pea proteomics needs to link numerous existing physiological and biochemical data to gain further insight into adaptation processes, which play important roles in field applications. Future developments and directions in pea proteomics are discussed.

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.

Deciphering common and specific transcriptional immune responses in pea towards the oomycete pathogens Aphanomyces euteiches and Phytophthora pisi

BACKGROUND

Root rot caused by Aphanomyces euteiches is one of the most destructive pea diseases while a distantly related species P. pisi has been recently described as the agent of pea and faba bean root rot. These two oomycete pathogens with different pathogenicity factor repertories have both evolved specific mechanisms to infect pea. However, little is known about the genes and mechanisms of defence against these pathogens in pea. In the present study, the transcriptomic response of pea to these two pathogens was investigated at two time points during early phase of infection using a Medicago truncatula microarray.

RESULTS

Of the 37,976 genes analysed, 574 and 817 were differentially expressed in response to A. euteiches at 6 hpi and 20 hpi, respectively, while 544 and 611 genes were differentially regulated against P. pisi at 6 hpi and 20 hpi, respectively. Differentially expressed genes associated with plant immunity responses were involved in cell wall reinforcement, hormonal signalling and phenylpropanoid metabolism. Activation of cell wall modification, regulation of jasmonic acid biosynthesis and induction of ethylene signalling pathway were among the common transcriptional responses to both of these oomycetes. However, induction of chalcone synthesis and the auxin pathway were specific transcriptional changes against A. euteiches.

CONCLUSIONS

Our results demonstrate a global view of differentially expressed pea genes during compatible interactions with P. pisi and A. euteiches at an early phase of infection. The results suggest that distinct signalling pathways are triggered in pea by these two pathogens that lead to common and specific immune mechanisms in response to these two oomycetes. The generated knowledge may eventually be used in breeding pea varieties with resistance against root rot disease.

De novo assembly and characterisation of the field pea transcriptome using RNA-Seq

BACKGROUND

Field pea (Pisum sativum L.) is a cool-season grain legume that is cultivated world-wide for both human consumption and stock-feed purposes. Enhancement of genetic and genomic resources for field pea will permit improved understanding of the control of traits relevant to crop productivity and quality. Advances in second-generation sequencing and associated bioinformatics analysis now provide unprecedented opportunities for the development of such resources. The objective of this study was to perform transcriptome sequencing and characterisation from two genotypes of field pea that differ in terms of seed and plant morphological characteristics.

RESULTS

Transcriptome sequencing was performed with RNA templates from multiple tissues of the field pea genotypes Kaspa and Parafield. Tissue samples were collected at various growth stages, and a total of 23 cDNA libraries were sequenced using Illumina high-throughput sequencing platforms. A total of 407 and 352 million paired-end reads from the Kaspa and Parafield transcriptomes, respectively were assembled into 129,282 and 149,272 contigs, which were filtered on the basis of known gene annotations, presence of open reading frames (ORFs), reciprocal matches and degree of coverage. Totals of 126,335 contigs from Kaspa and 145,730 from Parafield were subsequently selected as the reference set. Reciprocal sequence analysis revealed that c. 87% of contigs were expressed in both cultivars, while a small proportion were unique to each genotype. Reads from different libraries were aligned to the genotype-specific assemblies in order to identify and characterise expression of contigs on a tissue-specific basis, of which 87% were expressed in more than one tissue, while others showed distinct expression patterns in specific tissues, providing unique transcriptome signatures.

CONCLUSION

This study provided a comprehensive assembled and annotated transcriptome set for field pea that can be used for development of genetic markers, in order to assess genetic diversity, construct linkage maps, perform trait-dissection and implement whole-genome selection strategies in varietal improvement programs, as well to identify target genes for genetic modification approaches on the basis of annotation and expression analysis. In addition, the reference field pea transcriptome will prove highly valuable for comparative genomics studies and construction of a finalised genome sequence.

Deep sequencing of the Mexican avocado transcriptome, an ancient angiosperm with a high content of fatty acids

BACKGROUND

Avocado (Persea americana) is an economically important tropical fruit considered to be a good source of fatty acids. Despite its importance, the molecular and cellular characterization of biochemical and developmental processes in avocado is limited due to the lack of transcriptome and genomic information.

RESULTS

The transcriptomes of seeds, roots, stems, leaves, aerial buds and flowers were determined using different sequencing platforms. Additionally, the transcriptomes of three different stages of fruit ripening (pre-climacteric, climacteric and post-climacteric) were also analyzed. The analysis of the RNAseqatlas presented here reveals strong differences in gene expression patterns between different organs, especially between root and flower, but also reveals similarities among the gene expression patterns in other organs, such as stem, leaves and aerial buds (vegetative organs) or seed and fruit (storage organs). Important regulators, functional categories, and differentially expressed genes involved in avocado fruit ripening were identified. Additionally, to demonstrate the utility of the avocado gene expression atlas, we investigated the expression patterns of genes implicated in fatty acid metabolism and fruit ripening.

CONCLUSIONS

A description of transcriptomic changes occurring during fruit ripening was obtained in Mexican avocado, contributing to a dynamic view of the expression patterns of genes involved in fatty acid biosynthesis and the fruit ripening process.

The Medicago sativa gene index 1.2: a web-accessible gene expression atlas for investigating expression differences between Medicago sativa subspecies

Background

Alfalfa (Medicago sativa L.) is the primary forage legume crop species in the United States and plays essential economic and ecological roles in agricultural systems across the country. Modern alfalfa is the result of hybridization between tetraploid M. sativa ssp. sativa and M. sativa ssp. falcata. Due to its large and complex genome, there are few genomic resources available for alfalfa improvement.

Results

A de novo transcriptome assembly from two alfalfa subspecies, M. sativa ssp. sativa (B47) and M. sativa ssp. falcata (F56) was developed using Illumina RNA-seq technology. Transcripts from roots, nitrogen-fixing root nodules, leaves, flowers, elongating stem internodes, and post-elongation stem internodes were assembled into the Medicago sativa Gene Index 1.2 (MSGI 1.2) representing 112,626 unique transcript sequences. Nodule-specific and transcripts involved in cell wall biosynthesis were identified. Statistical analyses identified 20,447 transcripts differentially expressed between the two subspecies. Pair-wise comparisons of each tissue combination identified 58,932 sequences differentially expressed in B47 and 69,143 sequences differentially expressed in F56. Comparing transcript abundance in floral tissues of B47 and F56 identified expression differences in sequences involved in anthocyanin and carotenoid synthesis, which determine flower pigmentation. Single nucleotide polymorphisms (SNPs) unique to each M. sativa subspecies (110,241) were identified.

Conclusions

The Medicago sativa Gene Index 1.2 increases the expressed sequence data available for alfalfa by ninefold and can be expanded as additional experiments are performed. The MSGI 1.2 transcriptome sequences, annotations, expression profiles, and SNPs were assembled into the Alfalfa Gene Index and Expression Database (AGED) at http://plantgrn.noble.org/AGED/, a publicly available genomic resource for alfalfa improvement and legume research.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1718-7) contains supplementary material, which is available to authorized users.

The structure of plant photosystem I super-complex at 2.8 Å resolution

Most life forms on Earth are supported by solar energy harnessed by oxygenic photosynthesis. In eukaryotes, photosynthesis is achieved by large membrane-embedded super-complexes, containing reaction centers and connected antennae. Here, we report the structure of the higher plant PSI-LHCI super-complex determined at 2.8 Å resolution. The structure includes 16 subunits and more than 200 prosthetic groups, which are mostly light harvesting pigments. The complete structures of the four LhcA subunits of LHCI include 52 chlorophyll a and 9 chlorophyll b molecules, as well as 10 carotenoids and 4 lipids. The structure of PSI-LHCI includes detailed protein pigments and pigment-pigment interactions, essential for the mechanism of excitation energy transfer and its modulation in one of nature's most efficient photochemical machines.

De novo assembly of a genome-wide transcriptome map of Vicia faba (L.) for transfer cell research

Vicia faba (L.) is an important cool-season grain legume species used widely in agriculture but also in plant physiology research, particularly as an experimental model to study transfer cell (TC) development. TCs are specialized nutrient transport cells in plants, characterized by invaginated wall ingrowths with amplified plasma membrane surface area enriched with transporter proteins that facilitate nutrient transfer. Many TCs are formed by trans-differentiation from differentiated cells at apoplasmic/symplasmic boundaries in nutrient transport. Adaxial epidermal cells of isolated cotyledons can be induced to form functional TCs, thus providing a valuable experimental system to investigate genetic regulation of TC trans-differentiation. The genome of V. faba is exceedingly large (ca. 13 Gb), however, and limited genomic information is available for this species. To provide a resource for future transcript profiling of epidermal TC differentiation, we have undertaken de novo assembly of a genome-wide transcriptome map for V. faba. Illumina paired-end sequencing of total RNA pooled from different tissues and different stages, including isolated cotyledons induced to form epidermal TCs, generated 69.5 M reads, of which 65.8 M were used for assembly following trimming and quality control. Assembly using a De-Bruijn graph-based approach generated 21,297 contigs, of which 80.6% were successfully annotated against GO terms. The assembly was validated against known V. faba cDNAs held in GenBank, including transcripts previously identified as being specifically expressed in epidermal cells across TC trans-differentiation. This genome-wide transcriptome map therefore provides a valuable tool for future transcript profiling of epidermal TC trans-differentiation, and also enriches the genetic resources available for this important legume crop species.

Pea VEGETATIVE2 Is an FD Homolog That Is Essential for Flowering and Compound Inflorescence Development

As knowledge of the gene networks regulating inflorescence development in Arabidopsis thaliana improves, the current challenge is to characterize this system in different groups of crop species with different inflorescence architecture. Pea (Pisum sativum) has served as a model for development of the compound raceme, characteristic of many legume species, and in this study, we characterize the pea VEGETATIVE2 (VEG2) locus, showing that it is critical for regulation of flowering and inflorescence development and identifying it as a homolog of the bZIP transcription factor FD. Through detailed phenotypic characterizations of veg2 mutants, expression analyses, and the use of protein-protein interaction assays, we find that VEG2 has important roles during each stage of development of the pea compound inflorescence. Our results suggest that VEG2 acts in conjunction with multiple FLOWERING LOCUS T (FT) proteins to regulate expression of downstream target genes, including TERMINAL FLOWER1, LEAFY, and MADS box homologs, and to facilitate cross-regulation within the FT gene family. These findings further extend our understanding of the mechanisms underlying compound inflorescence development in pea and may have wider implications for future manipulation of inflorescence architecture in related legume crop species.

Oep23 forms an ion channel in the chloroplast outer envelope

BACKGROUND

Metabolite, ion and protein translocation into chloroplasts occurs across two membranes, the inner and the outer envelope. Solute and metabolite channels fulfill very important functions in integrating the organelles into the metabolic network of the cell. However so far only a few have been identified. Here we describe the identification and the characterization of the outer envelope protein of 23 kDa, Oep23 from garden pea.

RESULTS

Oep23 is found in the entire plant lineage from green algae to flowering plants. It is expressed in all organs and developmental states tested so far. The reconstituted recombinant protein Oep23 from pea forms a high conductance ion channel with a maximal conductance in the fully open state of 466 ± 14pS at a holding potential of +100 mV (in 250 mM KCl). The Oep23 channel is cation selective (PK+ : PCl- = 15 : 1) with a voltage dependent open probability of maximal Vmem = 0 mV.

CONCLUSION

The data indicate that the Oep23 activity represents a single channel unit and does not assemble into a multiple pore complex like bacterial type porins or mitochondrial voltage dependent anion channel. Thus, Oep23 represents a new member of ion channels in the outer envelope of chloroplasts involved in solute exchange.

Next-generation sequencing is a robust strategy for the high-throughput detection of zygosity in transgenic maize

Simple and reliable, high-throughput techniques to detect the zygosity of transgenic events in plants are valuable for biotechnology and plant breeding companies seeking robust genotyping data for the assessment of new lines and the monitoring of breeding programs. We show that next-generation sequencing (NGS) applied to short PCR products spanning the transgene integration site provides accurate zygosity data that are more robust and reliable than those generated by PCR-based methods. The NGS reads covered the 5' border of the transgenic events (incorporating part of the transgene and the flanking genomic DNA), or the genomic sequences flanking the unfilled transgene integration site at the wild-type locus. We compared the NGS method to competitive real-time PCR with transgene-specific and wild-type-specific primer/probe pairs, one pair matching the 5' genomic flanking sequence and 5' part of the transgene and the other matching the unfilled transgene integration site. Although both NGS and real-time PCR provided useful zygosity data, the NGS technique was favorable because it needed fewer optimization steps. It also provided statistically more-reliable evidence for the presence of each allele because each product was often covered by more than 100 reads. The NGS method is also more suitable for the genotyping of large panels of plants because up to 80 million reads can be produced in one sequencing run. Our novel method is therefore ideal for the rapid and accurate genotyping of large numbers of samples.

FAX1, a novel membrane protein mediating plastid fatty acid export

Fatty acid synthesis in plants occurs in plastids, and thus, export for subsequent acyl editing and lipid assembly in the cytosol and endoplasmatic reticulum is required. Yet, the transport mechanism for plastid fatty acids still remains enigmatic. We isolated FAX1 (fatty acid export 1), a novel protein, which inserts into the chloroplast inner envelope by α-helical membrane-spanning domains. Detailed phenotypic and ultrastructural analyses of FAX1 mutants in Arabidopsis thaliana showed that FAX1 function is crucial for biomass production, male fertility and synthesis of fatty acid-derived compounds such as lipids, ketone waxes, or pollen cell wall material. Determination of lipid, fatty acid, and wax contents by mass spectrometry revealed that endoplasmatic reticulum (ER)-derived lipids decreased when FAX1 was missing, but levels of several plastid-produced species increased. FAX1 over-expressing lines showed the opposite behavior, including a pronounced increase of triacyglycerol oils in flowers and leaves. Furthermore, the cuticular layer of stems from fax1 knockout lines was specifically reduced in C29 ketone wax compounds. Differential gene expression in FAX1 mutants as determined by DNA microarray analysis confirmed phenotypes and metabolic imbalances. Since in yeast FAX1 could complement for fatty acid transport, we concluded that FAX1 mediates fatty acid export from plastids. In vertebrates, FAX1 relatives are structurally related, mitochondrial membrane proteins of so-far unknown function. Therefore, this protein family might represent a powerful tool not only to increase lipid/biofuel production in plants but also to explore novel transport systems involved in vertebrate fatty acid and lipid metabolism.

The novel protein FAX1 mediates the export of free fatty acids across the inner membrane of chloroplasts so that they can be processed in other plant cell organelles to generate oils, waxes, and other lipids.

Fatty acid synthesis in plants occurs in chloroplasts—the organelle more commonly known for conducting photosynthesis. For subsequent lipid assembly to be possible in the endoplasmatic reticulum (ER), export of these fatty acids across the chloroplast envelope membranes is required. The mechanism of this transport until now has not been known. We isolated FAX1 (fatty acid export 1), a novel membrane protein in chloroplast inner envelopes. FAX1 function is crucial for biomass production, male fertility, and the synthesis of fatty acid-derived compounds like lipids, waxes, or cell wall material of pollen grains. Whereas ER-derived lipids decreased when FAX1 was missing, levels of plastid-produced lipids increased. FAX1 over-expressing mutants showed the opposite behavior, including an increase of triacyglycerol oils. Because FAX1 could complement for fatty acid transport in yeast, we concluded that FAX1 mediates the export of free fatty acids from chloroplasts. In vertebrates, FAX1 relatives are structurally related proteins of so-far unknown function in mitochondria. This protein family may thus represent a powerful tool not only to increase lipid oil and biofuel production in plants but also to explore novel transport systems in animals.

Comparative Transcriptomic Analyses of Vegetable and Grain Pea (Pisum sativum L.) Seed Development

Understanding the molecular mechanisms regulating pea seed developmental process is extremely important for pea breeding. In this study, we used high-throughput RNA-Seq and bioinformatics analyses to examine the changes in gene expression during seed development in vegetable pea and grain pea, and compare the gene expression profiles of these two pea types. RNA-Seq generated 18.7 G of raw data, which were then de novo assembled into 77,273 unigenes with a mean length of 930 bp. Our results illustrate that transcriptional control during pea seed development is a highly coordinated process. There were 459 and 801 genes differentially expressed at early and late seed maturation stages between vegetable pea and grain pea, respectively. Soluble sugar and starch metabolism related genes were significantly activated during the development of pea seeds coinciding with the onset of accumulation of sugar and starch in the seeds. A comparative analysis of genes involved in sugar and starch biosynthesis in vegetable pea (high seed soluble sugar and low starch) and grain pea (high seed starch and low soluble sugar) revealed that differential expression of related genes at late development stages results in a negative correlation between soluble sugar and starch biosynthetic flux in vegetable and grain pea seeds. RNA-Seq data was validated by using real-time quantitative RT-PCR analysis for 30 randomly selected genes. To our knowledge, this work represents the first report of seed development transcriptomics in pea. The obtained results provide a foundation to support future efforts to unravel the underlying mechanisms that control the developmental biology of pea seeds, and serve as a valuable resource for improving pea breeding.

Plastidic type I signal peptidase 1 is a redox-dependent thylakoidal processing peptidase

Thylakoids are the photosynthetic membranes in chloroplasts and cyanobacteria. The aqueous phase inside the thylakoid known as the thylakoid lumen plays an essential role in the photosynthetic electron transport. The presence and significance of thiol-disulfide exchange in this compartment have been recognized but remain poorly understood. All proteins found free in the thylakoid lumen and some proteins associated to the thylakoid membrane require an N-terminal targeting signal, which is removed in the lumen by a membrane-bound serine protease called thylakoidal processing peptidase (TPP). TPP is homologous to Escherichia coli type I signal peptidase (SPI) called LepB. Genetic data indicate that plastidic SPI 1 (Plsp1) is the main TPP in Arabidopsis thaliana (Arabidopsis) although biochemical evidence had been lacking. Here we demonstrate catalytic activity of bacterially produced Arabidopsis Plsp1. Recombinant Plsp1 showed processing activity against various TPP substrates at a level comparable to that of LepB. Plsp1 and LepB were also similar in the pH optima, sensitivity to arylomycin variants and a preference for the residue at -3 to the cleavage site within a substrate. Plsp1 orthologs found in angiosperms contain two unique Cys residues located in the lumen. Results of processing assays suggested that these residues were redox active and formation of a disulfide bond between them was necessary for the activity of recombinant Arabidopsis Plsp1. Furthermore, Plsp1 in Arabidopsis and pea thylakoids migrated faster under non-reducing conditions than under reducing conditions on SDS-PAGE. These results underpin the notion that Plsp1 is a redox-dependent signal peptidase in the thylakoid lumen.

Integrative structural annotation of de novo RNA-Seq provides an accurate reference gene set of the enormous genome of the onion (Allium cepa L.)

The onion (Allium cepa L.) is one of the most widely cultivated and consumed vegetable crops in the world. Although a considerable amount of onion transcriptome data has been deposited into public databases, the sequences of the protein-coding genes are not accurate enough to be used, owing to non-coding sequences intermixed with the coding sequences. We generated a high-quality, annotated onion transcriptome from de novo sequence assembly and intensive structural annotation using the integrated structural gene annotation pipeline (ISGAP), which identified 54,165 protein-coding genes among 165,179 assembled transcripts totalling 203.0 Mb by eliminating the intron sequences. ISGAP performed reliable annotation, recognizing accurate gene structures based on reference proteins, and ab initio gene models of the assembled transcripts. Integrative functional annotation and gene-based SNP analysis revealed a whole biological repertoire of genes and transcriptomic variation in the onion. The method developed in this study provides a powerful tool for the construction of reference gene sets for organisms based solely on de novo transcriptome data. Furthermore, the reference genes and their variation described here for the onion represent essential tools for molecular breeding and gene cloning in Allium spp.