The Medicago genome provides insight into the evolution of rhizobial symbioses

Fluorescence in situ hybridization and optical mapping to correct scaffold arrangement in the tomato genome

The order and orientation (arrangement) of all 91 sequenced scaffolds in the 12 pseudomolecules of the recently published tomato (Solanum lycopersicum, 2n = 2x = 24) genome sequence were positioned based on marker order in a high-density linkage map. Here, we report the arrangement of these scaffolds determined by two independent physical methods, bacterial artificial chromosome–fluorescence in situ hybridization (BAC-FISH) and optical mapping. By localizing BACs at the ends of scaffolds to spreads of tomato synaptonemal complexes (pachytene chromosomes), we showed that 45 scaffolds, representing one-third of the tomato genome, were arranged differently than predicted by the linkage map. These scaffolds occur mostly in pericentric heterochromatin where 77% of the tomato genome is located and where linkage mapping is less accurate due to reduced crossing over. Although useful for only part of the genome, optical mapping results were in complete agreement with scaffold arrangement by FISH but often disagreed with scaffold arrangement based on the linkage map. The scaffold arrangement based on FISH and optical mapping changes the positions of hundreds of markers in the linkage map, especially in heterochromatin. These results suggest that similar errors exist in pseudomolecules from other large genomes that have been assembled using only linkage maps to predict scaffold arrangement, and these errors can be corrected using FISH and/or optical mapping. Of note, BAC-FISH also permits estimates of the sizes of gaps between scaffolds, and unanchored BACs are often visualized by FISH in gaps between scaffolds and thus represent starting points for filling these gaps.

Genome variations account for different response to three mineral elements between Medicago truncatula ecotypes Jemalong A17 and R108

BACKGROUND

Resequencing can be used to identify genome variations underpinning many morphological and physiological phenotypes. Legume model plant Medicago truncatula ecotypes Jemalong A17 (J. A17) and R108 differ in their responses to mineral toxicity of aluminum and sodium, and mineral deficiency of iron in growth medium. The difference may result from their genome variations, but no experimental evidence supports this hypothesis.

RESULTS

A total of 12,750 structure variations, 135,045 short insertions/deletions and 764,154 single nucleotide polymorphisms were identified by resequencing the genome of R108. The suppressed expression of MtAACT that encodes a putative aluminum-induced citrate efflux transporter by deletion of partial sequence of the second intron may account for the less aluminum-induced citrate exudation and greater accumulation of aluminum in roots of R108 than in roots of J. A17, thus rendering R108 more sensitive to aluminum toxicity. The higher expression-level of MtZpt2-1 encoding a TFIIIA-related transcription factor in J. A17 than R108 under conditions of salt stress can be explained by the greater number of stress-responsive elements in its promoter sequence, thus conferring J. A17 more tolerant to salt stress than R108 plants by activating the expression of downstream stress-responsive genes. YSLs (Yellow Stripe-Likes) are involved in long-distance transport of iron in plants. We found that an YSL gene was deleted in the genome of R108 plants, thus rendering R108 less tolerance to iron deficiency than J. A17 plants.

CONCLUSIONS

The deletion or change in several genes may account for the different responses of M. truncatula ecotypes J. A17 and R108 to mineral toxicity of aluminum and sodium as well as iron deficiency. Uncovering genome variations by resequencing is an effective method to identify different traits between species/ecotypes that are genetically related. These findings demonstrate that analyses of genome variations by resequencing can shed important light on differences in responses of M. truncatula ecotypes to abiotic stress in general and mineral stress in particular.

Synthetic biology approaches to engineering the nitrogen symbiosis in cereals

Nitrogen is abundant in the earth's atmosphere but, unlike carbon, cannot be directly assimilated by plants. The limitation this places on plant productivity has been circumvented in contemporary agriculture through the production and application of chemical fertilizers. The chemical reduction of nitrogen for this purpose consumes large amounts of energy and the reactive nitrogen released into the environment as a result of fertilizer application leads to greenhouse gas emissions, as well as widespread eutrophication of aquatic ecosystems. The environmental impacts are intensified by injudicious use of fertilizers in many parts of the world. Simultaneously, limitations in the production and supply of chemical fertilizers in other regions are leading to low agricultural productivity and malnutrition. Nitrogen can be directly fixed from the atmosphere by some bacteria and Archaea, which possess the enzyme nitrogenase. Some plant species, most notably legumes, have evolved close symbiotic associations with nitrogen-fixing bacteria. Engineering cereal crops with the capability to fix their own nitrogen could one day address the problems created by the over- and under-use of nitrogen fertilizers in agriculture. This could be achieved either by expression of a functional nitrogenase enzyme in the cells of the cereal crop or through transferring the capability to form a symbiotic association with nitrogen-fixing bacteria. While potentially transformative, these biotechnological approaches are challenging; however, with recent advances in synthetic biology they are viable long-term goals. This review discusses the possibility of these biotechnological solutions to the nitrogen problem, focusing on engineering the nitrogen symbiosis in cereals.

New evidence of ancestral polyploidy in the Genistoid legume Lupinus angustifolius L. (narrow-leafed lupin)

This is the first clear evidence of duplication and/or triplication of large chromosomal regions in a genome of a Genistoid legume, the most basal clade of Papilionoid legumes. Lupinus angustifolius L. (narrow-leafed lupin) is the most widely cultivated species of Genistoid legume, grown for its high-protein grain. As a member of this most basal clade of Papilionoid legumes, L. angustifolius serves as a useful model for exploring legume genome evolution. Here, we report an improved reference genetic map of L. angustifolius comprising 1207 loci, including 299 newly developed Diversity Arrays Technology markers and 54 new gene-based PCR markers. A comparison between the L. angustifolius and Medicago truncatula genomes was performed using 394 sequence-tagged site markers acting as bridging points between the two genomes. The improved L. angustifolius genetic map, the updated M. truncatula genome assembly and the increased number of bridging points between the genomes together substantially enhanced the resolution of synteny and chromosomal colinearity between these genomes compared to previous reports. While a high degree of syntenic fragmentation was observed that was consistent with the large evolutionary distance between the L. angustifolius and M. truncatula genomes, there were striking examples of conserved colinearity of loci between these genomes. Compelling evidence was found of large-scale duplication and/or triplication in the L. angustifolius genome, consistent with one or more ancestral polyploidy events.

Patterns of divergence of a large family of nodule cysteine-rich peptides in accessions of Medicago truncatula

The nodule cysteine-rich (NCR) groups of defensin-like (DEFL) genes are one of the largest gene families expressed in the nodules of some legume plants. They have only been observed in the inverted repeat loss clade (IRLC) of legumes, which includes the model legume Medicago truncatula. NCRs are reported to play an important role in plant-microbe interactions. To understand their diversity we analyzed their expression and sequence polymorphisms among four accessions of M. truncatula. A significant expression and nucleotide variation was observed among the genes. We then used 26 accessions to estimate the selection pressures shaping evolution among the accessions by calculating the nucleotide diversity at non-synonymous and synonymous sites in the coding region. The mature peptides of the orthologous NCRs had signatures of both purifying and diversifying selection pressures, unlike the seed DEFLs, which predominantly exhibited purifying selection. The expression, sequence variation and apparent diversifying selection in NCRs within the Medicago species indicates rapid and recent evolution, and suggests that this family of genes is actively evolving to adapt to different environments and is acquiring new functions.

An improved genome release (version Mt4.0) for the model legume Medicago truncatula

BACKGROUND

Medicago truncatula, a close relative of alfalfa, is a preeminent model for studying nitrogen fixation, symbiosis, and legume genomics. The Medicago sequencing project began in 2003 with the goal to decipher sequences originated from the euchromatic portion of the genome. The initial sequencing approach was based on a BAC tiling path, culminating in a BAC-based assembly (Mt3.5) as well as an in-depth analysis of the genome published in 2011.

RESULTS

Here we describe a further improved and refined version of the M. truncatula genome (Mt4.0) based on de novo whole genome shotgun assembly of a majority of Illumina and 454 reads using ALLPATHS-LG. The ALLPATHS-LG scaffolds were anchored onto the pseudomolecules on the basis of alignments to both the optical map and the genotyping-by-sequencing (GBS) map. The Mt4.0 pseudomolecules encompass ~360 Mb of actual sequences spanning 390 Mb of which ~330 Mb align perfectly with the optical map, presenting a drastic improvement over the BAC-based Mt3.5 which only contained 70% sequences (~250 Mb) of the current version. Most of the sequences and genes that previously resided on the unanchored portion of Mt3.5 have now been incorporated into the Mt4.0 pseudomolecules, with the exception of ~28 Mb of unplaced sequences. With regard to gene annotation, the genome has been re-annotated through our gene prediction pipeline, which integrates EST, RNA-seq, protein and gene prediction evidences. A total of 50,894 genes (31,661 high confidence and 19,233 low confidence) are included in Mt4.0 which overlapped with ~82% of the gene loci annotated in Mt3.5. Of the remaining genes, 14% of the Mt3.5 genes have been deprecated to an "unsupported" status and 4% are absent from the Mt4.0 predictions.

CONCLUSIONS

Mt4.0 and its associated resources, such as genome browsers, BLAST-able datasets and gene information pages, can be found on the JCVI Medicago web site (http://www.jcvi.org/medicago). The assembly and annotation has been deposited in GenBank (BioProject: PRJNA10791). The heavily curated chromosomal sequences and associated gene models of Medicago will serve as a better reference for legume biology and comparative genomics.

An improved genome release (version Mt4.0) for the model legume Medicago truncatula

Background

Medicago truncatula, a close relative of alfalfa, is a preeminent model for studying nitrogen fixation, symbiosis, and legume genomics. The Medicago sequencing project began in 2003 with the goal to decipher sequences originated from the euchromatic portion of the genome. The initial sequencing approach was based on a BAC tiling path, culminating in a BAC-based assembly (Mt3.5) as well as an in-depth analysis of the genome published in 2011.

Results

Here we describe a further improved and refined version of the M. truncatula genome (Mt4.0) based on de novo whole genome shotgun assembly of a majority of Illumina and 454 reads using ALLPATHS-LG. The ALLPATHS-LG scaffolds were anchored onto the pseudomolecules on the basis of alignments to both the optical map and the genotyping-by-sequencing (GBS) map. The Mt4.0 pseudomolecules encompass ~360 Mb of actual sequences spanning 390 Mb of which ~330 Mb align perfectly with the optical map, presenting a drastic improvement over the BAC-based Mt3.5 which only contained 70% sequences (~250 Mb) of the current version. Most of the sequences and genes that previously resided on the unanchored portion of Mt3.5 have now been incorporated into the Mt4.0 pseudomolecules, with the exception of ~28 Mb of unplaced sequences. With regard to gene annotation, the genome has been re-annotated through our gene prediction pipeline, which integrates EST, RNA-seq, protein and gene prediction evidences. A total of 50,894 genes (31,661 high confidence and 19,233 low confidence) are included in Mt4.0 which overlapped with ~82% of the gene loci annotated in Mt3.5. Of the remaining genes, 14% of the Mt3.5 genes have been deprecated to an “unsupported” status and 4% are absent from the Mt4.0 predictions.

Conclusions

Mt4.0 and its associated resources, such as genome browsers, BLAST-able datasets and gene information pages, can be found on the JCVI Medicago web site (http://www.jcvi.org/medicago). The assembly and annotation has been deposited in GenBank (BioProject: PRJNA10791). The heavily curated chromosomal sequences and associated gene models of Medicago will serve as a better reference for legume biology and comparative genomics.

An improved genome release (version Mt4.0) for the model legume Medicago truncatula

BACKGROUND

Medicago truncatula, a close relative of alfalfa, is a preeminent model for studying nitrogen fixation, symbiosis, and legume genomics. The Medicago sequencing project began in 2003 with the goal to decipher sequences originated from the euchromatic portion of the genome. The initial sequencing approach was based on a BAC tiling path, culminating in a BAC-based assembly (Mt3.5) as well as an in-depth analysis of the genome published in 2011.

RESULTS

Here we describe a further improved and refined version of the M. truncatula genome (Mt4.0) based on de novo whole genome shotgun assembly of a majority of Illumina and 454 reads using ALLPATHS-LG. The ALLPATHS-LG scaffolds were anchored onto the pseudomolecules on the basis of alignments to both the optical map and the genotyping-by-sequencing (GBS) map. The Mt4.0 pseudomolecules encompass ~360 Mb of actual sequences spanning 390 Mb of which ~330 Mb align perfectly with the optical map, presenting a drastic improvement over the BAC-based Mt3.5 which only contained 70% sequences (~250 Mb) of the current version. Most of the sequences and genes that previously resided on the unanchored portion of Mt3.5 have now been incorporated into the Mt4.0 pseudomolecules, with the exception of ~28 Mb of unplaced sequences. With regard to gene annotation, the genome has been re-annotated through our gene prediction pipeline, which integrates EST, RNA-seq, protein and gene prediction evidences. A total of 50,894 genes (31,661 high confidence and 19,233 low confidence) are included in Mt4.0 which overlapped with ~82% of the gene loci annotated in Mt3.5. Of the remaining genes, 14% of the Mt3.5 genes have been deprecated to an "unsupported" status and 4% are absent from the Mt4.0 predictions.

CONCLUSIONS

Mt4.0 and its associated resources, such as genome browsers, BLAST-able datasets and gene information pages, can be found on the JCVI Medicago web site (http://www.jcvi.org/medicago). The assembly and annotation has been deposited in GenBank (BioProject: PRJNA10791). The heavily curated chromosomal sequences and associated gene models of Medicago will serve as a better reference for legume biology and comparative genomics.

Genetic analysis of tolerance to the root lesion nematode Pratylenchus neglectus in the legume Medicago littoralis

BACKGROUND

The nematode Pratylenchus neglectus has a wide host range and is able to feed on the root systems of cereals, oilseeds, grain and pasture legumes. Under the Mediterranean low rainfall environments of Australia, annual Medicago pasture legumes are used in rotation with cereals to fix atmospheric nitrogen and improve soil parameters. Considerable efforts are being made in breeding programs to improve resistance and tolerance to Pratylenchus neglectus in the major crops wheat and barley, which makes it vital to develop appropriate selection tools in medics.

RESULTS

A strong source of tolerance to root damage by the root lesion nematode (RLN) Pratylenchus neglectus had previously been identified in line RH-1 (strand medic, M. littoralis). Using RH-1, we have developed a single seed descent (SSD) population of 138 lines by crossing it to the intolerant cultivar Herald. After inoculation, RLN-associated root damage clearly segregated in the population. Genetic analysis was performed by constructing a genetic map using simple sequence repeat (SSR) and gene-based SNP markers. A highly significant quantitative trait locus (QTL), QPnTolMl.1, was identified explaining 49% of the phenotypic variation in the SSD population. All SSRs and gene-based markers in the QTL region were derived from chromosome 1 of the sequenced genome of the closely related species M. truncatula. Gene-based markers were validated in advanced breeding lines derived from the RH-1 parent and also a second RLN tolerance source, RH-2 (M. truncatula ssp. tricycla). Comparative analysis to sequenced legume genomes showed that the physical QTL interval exists as a synteny block in Lotus japonicus, common bean, soybean and chickpea. Furthermore, using the sequenced genome information of M. truncatula, the QTL interval contains 55 genes out of which five are discussed as potential candidate genes responsible for the mapped tolerance.

CONCLUSION

The closely linked set of SNP-based PCR markers is directly applicable to select for two different sources of RLN tolerance in breeding programs. Moreover, genome sequence information has allowed proposing candidate genes for further functional analysis and nominates QPnTolMl.1 as a target locus for RLN tolerance in economically important grain legumes, e.g. chickpea.

Cross-family translational genomics of abiotic stress-responsive genes between Arabidopsis and Medicago truncatula

Cross-species translation of genomic information may play a pivotal role in applying biological knowledge gained from relatively simple model system to other less studied, but related, genomes. The information of abiotic stress (ABS)-responsive genes in Arabidopsis was identified and translated into the legume model system, Medicago truncatula. Various data resources, such as TAIR/AtGI DB, expression profiles and literatures, were used to build a genome-wide list of ABS genes. tBlastX/BlastP similarity search tools and manual inspection of alignments were used to identify orthologous genes between the two genomes. A total of 1,377 genes were finally collected and classified into 18 functional criteria of gene ontology (GO). The data analysis according to the expression cues showed that there was substantial level of interaction among three major types (i.e., drought, salinity and cold stress) of abiotic stresses. In an attempt to translate the ABS genes between these two species, genomic locations for each gene were mapped using an in-house-developed comparative analysis platform. The comparative analysis revealed that fragmental colinearity, represented by only 37 synteny blocks, existed between Arabidopsis and M. truncatula. Based on the combination of E-value and alignment remarks, estimated translation rate was 60.2% for this cross-family translation. As a prelude of the functional comparative genomic approaches, in-silico gene network/interactome analyses were conducted to predict key components in the ABS responses, and one of the sub-networks was integrated with corresponding comparative map. The results demonstrated that core members of the sub-network were well aligned with previously reported ABS regulatory networks. Taken together, the results indicate that network-based integrative approaches of comparative and functional genomics are important to interpret and translate genomic information for complex traits such as abiotic stresses.

Transcriptome analysis of salt tolerant common bean (Phaseolus vulgaris L.) under saline conditions

Salinity is one of the important abiotic stress factors that limit crop production. Common bean, Phaseolus vulgaris L., a major protein source in developing countries, is highly affected by soil salinity and the information on genes that play a role in salt tolerance is scarce. We aimed to identify differentially expressed genes (DEGs) and related pathways by comprehensive analysis of transcriptomes of both root and leaf tissues of the tolerant genotype grown under saline and control conditions in hydroponic system. We have generated a total of 158 million high-quality reads which were assembled into 83,774 all-unigenes with a mean length of 813 bp and N50 of 1,449 bp. Among the all-unigenes, 58,171 were assigned with Nr annotations after homology analyses. It was revealed that 6,422 and 4,555 all-unigenes were differentially expressed upon salt stress in leaf and root tissues respectively. Validation of the RNA-seq quantifications (RPKM values) was performed by qRT-PCR (Quantitative Reverse Transcription PCR) analyses. Enrichment analyses of DEGs based on GO and KEGG databases have shown that both leaf and root tissues regulate energy metabolism, transmembrane transport activity, and secondary metabolites to cope with salinity. A total of 2,678 putative common bean transcription factors were identified and classified under 59 transcription factor families; among them 441 were salt responsive. The data generated in this study will help in understanding the fundamentals of salt tolerance in common bean and will provide resources for functional genomic studies.

A paradigm for endosymbiotic life: cell differentiation of Rhizobium bacteria provoked by host plant factors

Symbiosis between Rhizobium bacteria and legumes leads to the formation of the root nodule. The endosymbiotic bacteria reside in polyploid host cells as membrane-surrounded vesicles where they reduce atmospheric nitrogen to support plant growth by supplying ammonia in exchange for carbon sources and energy. The morphology and physiology of endosymbionts, despite their common function, are highly divergent in different hosts. In galegoid plants, the endosymbionts are terminally differentiated, uncultivable polyploid cells, with remarkably elongated and even branched Y-shaped cells. Bacteroid differentiation is controlled by host peptides, many of which have antibacterial activity and require the bacterial function of BacA. Although the precise and combined action of several hundred host peptides and BacA has yet to be discovered, similarities, especially to certain insect-bacterium symbioses involving likewise host peptides for manipulation of endosymbionts, suggest convergent evolution. Rhizobium-legume symbiosis provides a rich source of information for understanding host-controlled endosymbiotic life in eukaryotic cells.

Two distinct plastid genome configurations and unprecedented intraspecies length variation in the accD coding region in Medicago truncatula

We fully sequenced four and partially sequenced six additional plastid genomes of the model legume Medicago truncatula. Three accessions, Jemalong 2HA, Borung and Paraggio, belong to ssp. truncatula, and R108 to ssp. tricycla. We report here that the R108 ptDNA has a ∼45-kb inversion compared with the ptDNA in ssp. truncatula, mediated by a short, imperfect repeat. DNA gel blot analyses of seven additional ssp. tricycla accessions detected only one of the two alternative genome arrangements, represented by three and four accessions each. Furthermore, we found a variable number of repeats in the essential accD and ycf1 coding regions. The repeats within accD are recombinationally active, yielding variable-length insertions and deletions in the central part of the coding region. The length of ACCD was distinct in each of the 10 sequenced ecotypes, ranging between 650 and 796 amino acids. The repeats in the ycf1 coding region are also recombinationally active, yielding short indels in 10 regions of the reading frames. Thus, the plastid genome variability we report here could be linked to repeat-mediated genome rearrangements. However, the rate of recombination was sufficiently low, so that no heterogeneity of ptDNA could be observed in populations maintained by single-seed descent.

The WRKY Transcription Factor Genes in Lotus japonicus

WRKY transcription factor genes play critical roles in plant growth and development, as well as stress responses. WRKY genes have been examined in various higher plants, but they have not been characterized in Lotus japonicus. The recent release of the L. japonicus whole genome sequence provides an opportunity for a genome wide analysis of WRKY genes in this species. In this study, we identified 61 WRKY genes in the L. japonicus genome. Based on the WRKY protein structure, L. japonicus WRKY (LjWRKY) genes can be classified into three groups (I-III). Investigations of gene copy number and gene clusters indicate that only one gene duplication event occurred on chromosome 4 and no clustered genes were detected on chromosomes 3 or 6. Researchers previously believed that group II and III WRKY domains were derived from the C-terminal WRKY domain of group I. Our results suggest that some WRKY genes in group II originated from the N-terminal domain of group I WRKY genes. Additional evidence to support this hypothesis was obtained by Medicago truncatula WRKY (MtWRKY) protein motif analysis. We found that LjWRKY and MtWRKY group III genes are under purifying selection, suggesting that WRKY genes will become increasingly structured and functionally conserved.

Growth conditions determine the DNF2 requirement for symbiosis

Rhizobia and legumes are able to interact in a symbiotic way leading to the development of root nodules. Within nodules, rhizobia fix nitrogen for the benefit of the plant. These interactions are efficient because spectacularly high densities of nitrogen fixing rhizobia are maintained in the plant cells. DNF2, a Medicago truncatula gene has been described as required for nitrogen fixation, bacteroid's persistence and to prevent defense-like reactions in the nodules. This manuscript shows that a Rhizobium mutant unable to differentiate is not sufficient to trigger defense-like reactions in this organ. Furthermore, we show that the requirement of DNF2 for effective symbiosis can be overcome by permissive growth conditions. The dnf2 knockout mutants grown in vitro on agarose or Phytagel as gelling agents are able to produce nodules fixing nitrogen with the same efficiency as the wild-type. However, when agarose medium is supplemented with the plant defense elicitor ulvan, the dnf2 mutant recovers the fix- phenotype. Together, our data show that plant growth conditions impact the gene requirement for symbiotic nitrogen fixation and suggest that they influence the symbiotic suppression of defense reactions in nodules.

BPG3 is a novel chloroplast protein that involves the greening of leaves and related to brassinosteroid signaling

Brassinosteroids are plant steroid hormones that regulate plant organs and chloroplast development. The detailed molecular mechanism for plant development by BR signaling is yet to be revealed, and many points regarding the relationship between BR signaling and chloroplast development remain unknown. We identify here the dominant mutant Brz-insensitive-pale green3-1D (bpg3-1D) from the Arabidopsis FOX lines that show reduced sensitivity to the chlorophyll accumulation promoted by the BR biosynthesis inhibitor, Brassinazole (Brz), in the light. BPG3 encodes a novel chloroplast protein that is evolutionally conserved in bacteria, algae, and higher plants. The expression of BPG3 was induced by light and Brz. The inhibition of electron transport in photosystem II of the chloroplasts was detected in bpg3-1D. These results suggest that BPG3 played an important role in regulating photosynthesis in the chloroplast under BR signaling.

High-density genome-wide association mapping implicates an F-box encoding gene in Medicago truncatula resistance to Aphanomyces euteiches

• The use of quantitative disease resistance (QDR) is a promising strategy for promoting durable resistance to plant pathogens, but genes involved in QDR are largely unknown. To identify genetic components and accelerate improvement of QDR in legumes to the root pathogen Aphanomyces euteiches, we took advantage of both the recently generated massive genomic data for Medicago truncatula and natural variation of this model legume. • A high-density (≈5.1 million single nucleotide polymorphisms (SNPs)) genome-wide association study (GWAS) was performed with both in vitro and glasshouse phenotyping data collected for 179 lines. • GWAS identified several candidate genes and pinpointed two independent major loci on the top of chromosome 3 that were detected in both phenotyping methods. Candidate SNPs in the most significant locus (σ(A)²= 23%) were in the promoter and coding regions of an F-box protein coding gene. Subsequent qRT-PCR and bioinformatic analyses performed on 20 lines demonstrated that resistance is associated with mutations directly affecting the interaction domain of the F-box protein rather than gene expression. • These results refine the position of previously identified QTL to specific candidate genes, suggest potential molecular mechanisms, and identify new loci explaining QDR against A. euteiches.

An integrated analysis of plant and bacterial gene expression in symbiotic root nodules using laser-capture microdissection coupled to RNA sequencing

Rhizobium-induced root nodules are specialized organs for symbiotic nitrogen fixation. Indeterminate-type nodules are formed from an apical meristem and exhibit a spatial zonation which corresponds to successive developmental stages. To get a dynamic and integrated view of plant and bacterial gene expression associated with nodule development, we used a sensitive and comprehensive approach based upon oriented high-depth RNA sequencing coupled to laser microdissection of nodule regions. This study, focused on the association between the model legume Medicago truncatula and its symbiont Sinorhizobium meliloti, led to the production of 942 million sequencing read pairs that were unambiguously mapped on plant and bacterial genomes. Bioinformatic and statistical analyses enabled in-depth comparison, at a whole-genome level, of gene expression in specific nodule zones. Previously characterized symbiotic genes displayed the expected spatial pattern of expression, thus validating the robustness of our approach. We illustrate the use of this resource by examining gene expression associated with three essential elements of nodule development, namely meristem activity, cell differentiation and selected signaling processes related to bacterial Nod factors and redox status. We found that transcription factor genes essential for the control of the root apical meristem were also expressed in the nodule meristem, while the plant mRNAs most enriched in nodules compared with roots were mostly associated with zones comprising both plant and bacterial partners. The data, accessible on a dedicated website, represent a rich resource for microbiologists and plant biologists to address a variety of questions of both fundamental and applied interest.

Integrated physical, genetic and genome map of chickpea (Cicer arietinum L.)

Physical map of chickpea was developed for the reference chickpea genotype (ICC 4958) using bacterial artificial chromosome (BAC) libraries targeting 71,094 clones (~12× coverage). High information content fingerprinting (HICF) of these clones gave high-quality fingerprinting data for 67,483 clones, and 1,174 contigs comprising 46,112 clones and 3,256 singletons were defined. In brief, 574 Mb genome size was assembled in 1,174 contigs with an average of 0.49 Mb per contig and 3,256 singletons represent 407 Mb genome. The physical map was linked with two genetic maps with the help of 245 BAC-end sequence (BES)-derived simple sequence repeat (SSR) markers. This allowed locating some of the BACs in the vicinity of some important quantitative trait loci (QTLs) for drought tolerance and reistance to Fusarium wilt and Ascochyta blight. In addition, fingerprinted contig (FPC) assembly was also integrated with the draft genome sequence of chickpea. As a result, ~965 BACs including 163 minimum tilling path (MTP) clones could be mapped on eight pseudo-molecules of chickpea forming 491 hypothetical contigs representing 54,013,992 bp (~54 Mb) of the draft genome. Comprehensive analysis of markers in abiotic and biotic stress tolerance QTL regions led to identification of 654, 306 and 23 genes in drought tolerance "QTL-hotspot" region, Ascochyta blight resistance QTL region and Fusarium wilt resistance QTL region, respectively. Integrated physical, genetic and genome map should provide a foundation for cloning and isolation of QTLs/genes for molecular dissection of traits as well as markers for molecular breeding for chickpea improvement.

Population differentiation for germination and early seedling root growth traits under saline conditions in the annual legume Medicago truncatula (Fabaceae)

PREMISE OF THE STUDY

Seedling establishment and survival are highly sensitive to soil salinity and plants that evolved in saline environments are likely to express traits that increase fitness in those environments. Such traits are of ecological interest and they may have practical value for improving salt tolerance in cultivated species. We examined responses to soil salinity and tested potential mechanisms of salt tolerance in Medicago truncatula, using genotypes that originated from natural populations occurring on saline and nonsaline soils.

METHODS

Germination and seedling responses were quantified and compared between saline and nonsaline origin genotypes. Germination treatments included a range of sodium chloride (NaCl) concentrations in both offspring and parental environments. Seedling treatments included NaCl, abscisic acid (ABA), and potassium chloride (KCl).

KEY RESULTS

Saline origin genotypes displayed greater salinity tolerance for germination and seedling traits relative to nonsaline origin genotypes. We observed population specific differences for the effects of salinity on time to germination and for the impact of parental environment on germination rates. ABA and NaCl treatments had similar negative effects on root growth, although relative sensitivities differed, with saline population less sensitive to NaCl and more sensitive to ABA compared to their nonsaline counterparts.

CONCLUSIONS

We report population differentiation for germination and seedling growth traits under saline conditions among populations derived from saline and nonsaline environments. These observations are consistent with a syndrome of adaptations for salinity tolerance during early plant development, including traits that are common among saline environments and those that are idiosyncratic to local populations.

Genome sequence of the date palm Phoenix dactylifera L

Date palm (Phoenix dactylifera L.) is a cultivated woody plant species with agricultural and economic importance. Here we report a genome assembly for an elite variety (Khalas), which is 605.4 Mb in size and covers >90% of the genome (~671 Mb) and >96% of its genes (~41,660 genes). Genomic sequence analysis demonstrates that P. dactylifera experienced a clear genome-wide duplication after either ancient whole genome duplications or massive segmental duplications. Genetic diversity analysis indicates that its stress resistance and sugar metabolism-related genes tend to be enriched in the chromosomal regions where the density of single-nucleotide polymorphisms is relatively low. Using transcriptomic data, we also illustrate the date palm’s unique sugar metabolism that underlies fruit development and ripening. Our large-scale genomic and transcriptomic data pave the way for further genomic studies not only on P. dactylifera but also other Arecaceae plants.

The date palm is one of the most economically important plants of the palm family. Here, the authors present a high-quality genome assembly of the date palm Phoenix dactylifera, and reveal insights into the unique sugar metabolism underlying fruit ripening.

Genomic survey, gene expression analysis and structural modeling suggest diverse roles of DNA methyltransferases in legumes

DNA methylation plays a crucial role in development through inheritable gene silencing. Plants possess three types of DNA methyltransferases (MTases), namely Methyltransferase (MET), Chromomethylase (CMT) and Domains Rearranged Methyltransferase (DRM), which maintain methylation at CG, CHG and CHH sites. DNA MTases have not been studied in legumes so far. Here, we report the identification and analysis of putative DNA MTases in five legumes, including chickpea, soybean, pigeonpea, Medicago and Lotus. MTases in legumes could be classified in known MET, CMT, DRM and DNA nucleotide methyltransferases (DNMT2) subfamilies based on their domain organization. First three MTases represent DNA MTases, whereas DNMT2 represents a transfer RNA (tRNA) MTase. Structural comparison of all the MTases in plants with known MTases in mammalian and plant systems have been reported to assign structural features in context of biological functions of these proteins. The structure analysis clearly specified regions crucial for protein-protein interactions and regions important for nucleosome binding in various domains of CMT and MET proteins. In addition, structural model of DRM suggested that circular permutation of motifs does not have any effect on overall structure of DNA methyltransferase domain. These results provide valuable insights into role of various domains in molecular recognition and should facilitate mechanistic understanding of their function in mediating specific methylation patterns. Further, the comprehensive gene expression analyses of MTases in legumes provided evidence of their role in various developmental processes throughout the plant life cycle and response to various abiotic stresses. Overall, our study will be very helpful in establishing the specific functions of DNA MTases in legumes.

Dynamic recruitment of amino acid transporters to the insect/symbiont interface

Symbiosis is well known to influence bacterial symbiont genome evolution and has recently been shown to shape eukaryotic host genomes. Intriguing patterns of host genome evolution, including remarkable numbers of gene duplications, have been observed in the pea aphid, a sap-feeding insect that relies on a bacterial endosymbiont for amino acid provisioning. Previously, we proposed that gene duplication has been important for the evolution of symbiosis based on aphid-specific gene duplication in amino acid transporters (AATs), with some paralogs highly expressed in the cells housing symbionts (bacteriocytes). Here, we use a comparative approach to test the role of gene duplication in enabling recruitment of AATs to bacteriocytes. Using genomic and transcriptomic data, we annotate AATs from sap-feeding and non sap-feeding insects and find that, like aphids, AAT gene families have undergone independent large-scale gene duplications in three of four additional sap-feeding insects. RNA-seq differential expression data indicate that, like aphids, the sap-feeding citrus mealybug possesses several lineage-specific bacteriocyte-enriched paralogs. Further, differential expression data combined with quantitative PCR support independent evolution of bacteriocyte enrichment in sap-feeding insect AATs. Although these data indicate that gene duplication is not necessary to initiate host/symbiont amino acid exchange, they support a role for gene duplication in enabling AATs to mediate novel host/symbiont interactions broadly in the sap-feeding suborder Sternorrhyncha. In combination with recent studies on other symbiotic systems, gene duplication is emerging as a general pattern in host genome evolution.

Response of Medicago truncatula seedlings to colonization by Salmonella enterica and Escherichia coli O157:H7

Disease outbreaks due to the consumption of legume seedlings contaminated with human enteric bacterial pathogens like Escherichia coli O157:H7 and Salmonella enterica are reported every year. Besides contaminations occurring during food processing, pathogens present on the surface or interior of plant tissues are also responsible for such outbreaks. In the present study, surface and internal colonization of Medicago truncatula, a close relative of alfalfa, by Salmonella enterica and Escherichia coli O157:H7 were observed even with inoculum levels as low as two bacteria per plant. Furthermore, expression analyses revealed that approximately 30% of Medicago truncatula genes were commonly regulated in response to both of these enteric pathogens. This study highlights that very low inoculum doses trigger responses from the host plant and that both of these human enteric pathogens may in part use similar mechanisms to colonize legume seedlings.

Host plant peptides elicit a transcriptional response to control the Sinorhizobium meliloti cell cycle during symbiosis

The α-proteobacterium Sinorhizobium meliloti establishes a chronic intracellular infection during the symbiosis with its legume hosts. Within specialized host cells, S. meliloti differentiates into highly polyploid, enlarged nitrogen-fixing bacteroids. This differentiation is driven by host cells through the production of defensin-like peptides called "nodule-specific cysteine-rich" (NCR) peptides. Recent research has shown that synthesized NCR peptides exhibit antimicrobial activity at high concentrations but cause bacterial endoreduplication at sublethal concentrations. We leveraged synchronized S. meliloti populations to determine how treatment with a sublethal NCR peptide affects the cell cycle and physiology of bacteria at the molecular level. We found that at sublethal levels a representative NCR peptide specifically blocks cell division and antagonizes Z-ring function. Gene-expression profiling revealed that the cell division block was produced, in part, through the substantial transcriptional response elicited by sublethal NCR treatment that affected ∼15% of the genome. Expression of critical cell-cycle regulators, including ctrA, and cell division genes, including genes required for Z-ring function, were greatly attenuated in NCR-treated cells. In addition, our experiments identified important symbiosis functions and stress responses that are induced by sublethal levels of NCR peptides and other antimicrobial peptides. Several of these stress-response pathways also are found in related α-proteobacterial pathogens and might be used by S. meliloti to sense host cues during infection. Our data suggest a model in which, in addition to provoking stress responses, NCR peptides target intracellular regulatory pathways to drive S. meliloti endoreduplication and differentiation during symbiosis.

Genome assembly and annotation for red clover (Trifolium pratense; Fabaceae)

PREMISE OF THE STUDY

Red clover (Trifolium pratense) is an important forage plant from the legume family with great importance in agronomy and livestock nourishment. Nevertheless, assembling its medium-sized genome presents a challenge, given current hardware and software possibilities. Next-generation sequencing technologies enable us to generate large amounts of sequence data at low cost. In this study, the genome assembly and red clover genome features are presented.

METHODS

First, assembly software was assessed using data sets from a closely related species to find the best possible combination of assembler plus error correction program to assemble the red clover genome. The newly sequenced genome was characterized by repetitive content, number of protein-coding and nonprotein-coding genes, and gene families and functions. Genome features were also compared with those of other sequenced plant species.

KEY RESULTS

Abyss with Echo correction was used for de novo assembly of the red clover genome. The presented assembly comprises ∼314.6 Mbp. In contrast to leguminous species with comparable genome sizes, the genome of T. pratense contains a larger repetitive portion and more abundant retrotransposons and DNA transposons. Overall, 47 398 protein-coding genes were annotated from 64 761 predicted genes. Comparative analysis revealed several gene families that are characteristic for T. pratense. Resistance genes, leghemoglobins, and nodule-specific cystein-rich peptides were identified and compared with other sequenced species.

CONCLUSIONS

The presented red clover genomic data constitute a resource for improvement through molecular breeding and for comparison to other sequenced plant species.

A nodule-specific lipid transfer protein AsE246 participates in transport of plant-synthesized lipids to symbiosome membrane and is essential for nodule organogenesis in Chinese milk vetch

Rhizobia in legume root nodules fix nitrogen in symbiosomes, organelle-like structures in which a membrane from the host plant surrounds the symbiotic bacteria. However, the components that transport plant-synthesized lipids to the symbiosome membrane remain unknown. This study identified and functionally characterized the Chinese milk vetch (Astragalus sinicus) lipid transfer protein AsE246, which is specifically expressed in nodules. It was found that AsE246 can bind lipids in vitro. More importantly, AsE246 can bind the plant-synthesized membrane lipid digalactosyldiacylglycerol in vivo. Immunofluorescence and immunoelectron microscopy showed that AsE246 and digalactosyldiacylglycerol localize in the symbiosome membrane and are present in infection threads. Overexpression of AsE246 resulted in increased nodule numbers; knockdown of AsE246 resulted in reduced nodule numbers, decreased lipids contents in nodules, diminished nitrogen fixation activity, and abnormal development of symbiosomes. AsE246 knockdown also resulted in fewer infection threads, nodule primordia, and nodules, while AsE246 overexpression resulted in more infection threads and nodule primordia, suggesting that AsE246 affects nodule organogenesis associated with infection thread formation. Taken together, these results indicate that AsE246 contributes to lipids transport to the symbiosome membrane, and this transport is required for effective legume-rhizobium symbiosis.

An efficient reverse genetics platform in the model legume Medicago truncatula

Medicago truncatula is one of the model species for legume studies. In an effort to develop legume genetics resources, > 21 700 Tnt1 retrotransposon insertion lines have been generated. To facilitate fast-growing needs in functional genomics, two reverse genetics approaches have been established: web-based database searching and PCR-based reverse screening. More than 840 genes have been reverse screened using the PCR-based approach over the past 6 yr to identify mutants in these genes. Overall, c. 84% (705 genes) success rate was achieved in identifying mutants with at least one Tnt1 insertion, of which c. 50% (358 genes) had three or more alleles. To demonstrate the utility of the two reverse genetics platforms, two mutant alleles were isolated for each of the two floral homeotic MADS-box genes, MtPISTILATA and MtAGAMOUS. Molecular and genetic analyses indicate that Tnt1 insertions in exons of both genes are responsible for the defects in floral organ development. In summary, we have developed two efficient reverse genetics platforms to facilitate functional characterization of M. truncatula genes.

Evolution of a symbiotic receptor through gene duplications in the legume-rhizobium mutualism

The symbiosis between legumes and nitrogen-fixing rhizobia co-opted pre-existing endomycorrhizal features. In particular, both symbionts release lipo-chitooligosaccharides (LCOs) that are recognized by LysM-type receptor kinases. We investigated the evolutionary history of rhizobial LCO receptor genes MtLYK3-LjNFR1 to gain insight into the evolutionary origin of the rhizobial symbiosis. We performed a phylogenetic analysis integrating gene copies from nonlegumes and legumes, including the non-nodulating, phylogenetically basal legume Cercis chinensis. Signatures of differentiation between copies were investigated through patterns of molecular evolution. We show that two rounds of duplication preceded the evolution of the rhizobial symbiosis in legumes. Molecular evolution patterns indicate that the resulting three paralogous gene copies experienced different selective constraints. In particular, one copy maintained the ancestral function, and another specialized into perception of rhizobial LCOs. It has been suggested that legume LCO receptors evolved from a putative ancestral defense-related chitin receptor through the acquisition of two kinase motifs. However, the phylogenetic analysis shows that these domains are actually ancestral, suggesting that this scenario is unlikely. Our study underlines the evolutionary significance of gene duplication and subsequent neofunctionalization in MtLYK3-LjNFR1 genes. We hypothesize that their ancestor was more likely a mycorrhizal LCO receptor, than a defense-related receptor kinase.

Overexpression of Medicago SVP genes causes floral defects and delayed flowering in Arabidopsis but only affects floral development in Medicago

The MADS-domain transcription factor SHORT VEGETATIVE PHASE plays a key role as a repressor of the transition to flowering and as a regulator of early floral development in Arabidopsis thaliana (Arabidopsis). However, no flowering-time repressors have been functionally identified in the model legume Medicago truncatula (Medicago). In this study, phylogenetic analysis of two closely-related MtSVP-like sequences, MtSVP1 and MtSVP2, showed that their predicted proteins clustered together within the eudicot SVP clade. To determine if the MtSVP-like genes have a role in flowering, they were functionally characterized in Medicago and Arabidopsis. Transcripts of both MtSVP genes were abundant and broadly expressed in vegetative tissues but were detected at much lower levels in flowers in Medicago. Over-expression of the MtSVP genes in Arabidopsis resulted in delayed flowering and flowers with many abnormal phenotypes such as leafy sepals, changes to floral organ number and longer pedicels than the wild type. By contrast, in transgenic Medicago, over-expression of MtSVP1 resulted in alterations to flower development, but did not alter flowering time, suggesting that MtSVP1 may not function to repress the transition to flowering in Medicago.

Comprehensive transcriptome assembly of Chickpea (Cicer arietinum L.) using sanger and next generation sequencing platforms: development and applications

A comprehensive transcriptome assembly of chickpea has been developed using 134.95 million Illumina single-end reads, 7.12 million single-end FLX/454 reads and 139,214 Sanger expressed sequence tags (ESTs) from >17 genotypes. This hybrid transcriptome assembly, referred to as Cicer arietinumTranscriptome Assembly version 2 (CaTA v2, available at http://data.comparative-legumes.org/transcriptomes/cicar/lista_cicar-201201), comprising 46,369 transcript assembly contigs (TACs) has an N50 length of 1,726 bp and a maximum contig size of 15,644 bp. Putative functions were determined for 32,869 (70.8%) of the TACs and gene ontology assignments were determined for 21,471 (46.3%). The new transcriptome assembly was compared with the previously available chickpea transcriptome assemblies as well as to the chickpea genome. Comparative analysis of CaTA v2 against transcriptomes of three legumes - Medicago, soybean and common bean, resulted in 27,771 TACs common to all three legumes indicating strong conservation of genes across legumes. CaTA v2 was also used for identification of simple sequence repeats (SSRs) and intron spanning regions (ISRs) for developing molecular markers. ISRs were identified by aligning TACs to the Medicago genome, and their putative mapping positions at chromosomal level were identified using transcript map of chickpea. Primer pairs were designed for 4,990 ISRs, each representing a single contig for which predicted positions are inferred and distributed across eight linkage groups. A subset of randomly selected ISRs representing all eight chickpea linkage groups were validated on five chickpea genotypes and showed 20% polymorphism with average polymorphic information content (PIC) of 0.27. In summary, the hybrid transcriptome assembly developed and novel markers identified can be used for a variety of applications such as gene discovery, marker-trait association, diversity analysis etc., to advance genetics research and breeding applications in chickpea and other related legumes.

Genomic signature of adaptation to climate in Medicago truncatula

Local adaptation and adaptive clines are pervasive in natural plant populations, yet the effects of these types of adaptation on genomic diversity are not well understood. With a data set of 202 accessions of Medicago truncatula genotyped at almost 2 million single nucleotide polymorphisms, we used mixed linear models to identify candidate loci responsible for adaptation to three climatic gradients-annual mean temperature (AMT), precipitation in the wettest month (PWM), and isothermality (ITH)-representing the major axes of climate variation across the species' range. Loci with the strongest association to these climate gradients tagged genome regions with high sequence similarity to genes with functional roles in thermal tolerance, drought tolerance, or resistance to herbivores of pathogens. Genotypes at these candidate loci also predicted the performance of an independent sample of plant accessions grown in climate-controlled conditions. Compared to a genome-wide sample of randomly drawn reference SNPs, candidates for two climate gradients, AMT and PWM, were significantly enriched for genic regions, and genome segments flanking genic AMT and PWM candidates harbored less nucleotide diversity, elevated differentiation between haplotypes carrying alternate alleles, and an overrepresentation of the most common haplotypes. These patterns of diversity are consistent with a history of soft selective sweeps acting on loci underlying adaptation to climate, but not with a history of long-term balancing selection.

Keel petal incision: a simple and efficient method for genetic crossing in Medicago truncatula

BACKGROUND

Genetic crossing is an essential tool in both forward and reverse genetic approaches to understand the biological functions of genes. For Medicago truncatula (barrel medic) various crossing techniques have been used which differ in the methods used to dissect the female parent's unopened flower bud to remove immature anthers for prevention of self-pollination. Previously described methods including front, side or back incision methods may damage the flower bud, impeding successful fertilization and/or seed development because they may allow pollen to dislodge and floral organs to desiccate after crossing, all of which diminish the success rates of crossing.

RESULTS

We report the keel petal incision method for genetic crossing in M. truncatula ecotype R108 and demonstrate successful crosses with two other M. truncatula ecotypes, A17 and A20. In the method presented here, an incision is made along the central line of the keel petal from the bottom 1/3rd of the female parent's flower bud to its distal end. This allows easy removal of anthers from the flower bud and access for cross-pollination. After pollination, the stigma and the deposited pollen from the male donor are covered by the keel petal, wing petals and standard petal, forming a natural pouch. The pouch prevents dislodging of deposited pollen from the stigma and protects the internal floral organs from drying out, without using cling-film or water-containing chambers to maintain a humid environment. The keel petal incision method showed an approximate 80% success rate in the M. truncatula R108 ecotype and also in other ecotypes including Jemalong A17 and A20.

CONCLUSIONS

Our keel petal incision protocol shows marked improvement over existing methods with respect to the ease of crossing and the percentage of successful crosses. Developed for the M. truncatula R108 ecotype, the protocol has been demonstrated with A17 and A20 ecotypes and is expected to work with other ecotypes. Investigators of varying experience have achieved genetic crosses in M. truncatula using this method.

Genome-wide analysis and expression profiling suggest diverse roles of GH3 genes during development and abiotic stress responses in legumes

Growth hormone auxin regulates various cellular processes by altering the expression of diverse genes in plants. Among various auxin-responsive genes, GH3 genes maintain endogenous auxin homeostasis by conjugating excess of auxin with amino acids. GH3 genes have been characterized in many plant species, but not in legumes. In the present work, we identified members of GH3 gene family and analyzed their chromosomal distribution, gene structure, gene duplication and phylogenetic analysis in different legumes, including chickpea, soybean, Medicago, and Lotus. A comprehensive expression analysis in different vegetative and reproductive tissues/stages revealed that many of GH3 genes were expressed in a tissue-specific manner. Notably, chickpea CaGH3-3, soybean GmGH3-8 and -25, and Lotus LjGH3-4, -5, -9 and -18 genes were up-regulated in root, indicating their putative role in root development. In addition, chickpea CaGH3-1 and -7, and Medicago MtGH3-7, -8, and -9 were found to be highly induced under drought and/or salt stresses, suggesting their role in abiotic stress responses. We also observed the examples of differential expression pattern of duplicated GH3 genes in soybean, indicating their functional diversification. Furthermore, analyses of three-dimensional structures, active site residues and ligand preferences provided molecular insights into function of GH3 genes in legumes. The analysis presented here would help in investigation of precise function of GH3 genes in legumes during development and stress conditions.

Regulation of Compound Leaf Development

Leaf morphology is one of the most variable, yet inheritable, traits in the plant kingdom. How plants develop a variety of forms and shapes is a major biological question. Here, we discuss some recent progress in understanding the development of compound or dissected leaves in model species, such as tomato (Solanum lycopersicum), Cardamine hirsuta and Medicago truncatula, with an emphasis on recent discoveries in legumes. We also discuss progress in gene regulations and hormonal actions in compound leaf development. These studies facilitate our understanding of the underlying regulatory mechanisms and put forward a prospective in compound leaf studies.

The genome of the recently domesticated crop plant sugar beet (Beta vulgaris)

Sugar beet (Beta vulgaris ssp. vulgaris) is an important crop of temperate climates which provides nearly 30% of the world's annual sugar production and is a source for bioethanol and animal feed. The species belongs to the order of Caryophylalles, is diploid with 2n = 18 chromosomes, has an estimated genome size of 714-758 megabases and shares an ancient genome triplication with other eudicot plants. Leafy beets have been cultivated since Roman times, but sugar beet is one of the most recently domesticated crops. It arose in the late eighteenth century when lines accumulating sugar in the storage root were selected from crosses made with chard and fodder beet. Here we present a reference genome sequence for sugar beet as the first non-rosid, non-asterid eudicot genome, advancing comparative genomics and phylogenetic reconstructions. The genome sequence comprises 567 megabases, of which 85% could be assigned to chromosomes. The assembly covers a large proportion of the repetitive sequence content that was estimated to be 63%. We predicted 27,421 protein-coding genes supported by transcript data and annotated them on the basis of sequence homology. Phylogenetic analyses provided evidence for the separation of Caryophyllales before the split of asterids and rosids, and revealed lineage-specific gene family expansions and losses. We sequenced spinach (Spinacia oleracea), another Caryophyllales species, and validated features that separate this clade from rosids and asterids. Intraspecific genomic variation was analysed based on the genome sequences of sea beet (Beta vulgaris ssp. maritima; progenitor of all beet crops) and four additional sugar beet accessions. We identified seven million variant positions in the reference genome, and also large regions of low variability, indicating artificial selection. The sugar beet genome sequence enables the identification of genes affecting agronomically relevant traits, supports molecular breeding and maximizes the plant's potential in energy biotechnology.

Structure and expression of barley starch phosphorylase genes

The function of starch phosphorylase has long been debated on the regulation of starch metabolism during the growth and development of plants. In this study, we isolated starch phosphorylase genes (Pho1 and Pho2) from barley, characterized their gene and protein structures, predicated their promoter's cis-elements and analyzed expression patterns. Multiple alignments of these genes showed that (1) both Pho1 and Pho2 genes possess 15 exons and 14 introns in all but three of the species analyzed, Aegilops tauschii (for Pho1 which contains 16 exons and 15 introns), potato (for Pho1b which contains 14 exons and 13 introns), and Triticum uraru (for Pho2 which contains 15 exons and 14 introns); (2) the exon-intron junctions of Pho1 and Pho2 flanking the ligand-binding sites are more conservative than the other regions. Analysis of protein sequences revealed that Pho1 and Pho2 were highly homologous except for two regions, the N terminal domain and the L78 insertion region. The results of real-time quantitative PCR (RT-qPCR) indicated that Pho2 is mainly expressed in germinating seeds, and the expression of Pho1 is similar to that of starch synthesis genes during seed development in barley. Microarray-based analysis indicated that the accumulation of Pho1 or Pho2 transcripts exhibited uniform pattern both in various tissues and various stages of seed development among species of barley, rice, and Arabidopsis. Pho1 of barley was significantly down-regulated under cold and drought treatments, and up-regulated under stem rust infection. Pho2 exhibited similar expression to Pho1 in barley. However, significant difference in expression was not detected for either Pho1 or Pho2 under any of the investigated abiotic stresses. In Arabidopsis, significant down-regulation was detected for Pho1 (PHS1) under abscisic acid (ABA) and for Pho2 (PHS2) under cold, salt, and ABA. Our results provide valuable information to genetically manipulate phosphorylase genes and to further elucidate their regulatory mechanism in the starch biosynthetic pathway.

The peptide-encoding CEP1 gene modulates lateral root and nodule numbers in Medicago truncatula

The role of MtCEP1, a member of the CEP (C-terminally encoded peptide) signaling peptide family, was examined in Medicago truncatula root development. MtCEP1 was expressed in root tips, vascular tissue, and young lateral organs, and was up-regulated by low nitrogen levels and, independently, by elevated CO2. Overexpressing MtCEP1 or applying MtCEP1 peptide to roots elicited developmental phenotypes: inhibition of lateral root formation, enhancement of nodulation, and the induction of periodic circumferential root swellings, which arose from cortical, epidermal, and pericycle cell divisions and featured an additional cortical cell layer. MtCEP peptide addition to other legume species induced similar phenotypes. The enhancement of nodulation by MtCEP1 is partially tolerant to high nitrate, which normally strongly suppresses nodulation. These nodules develop faster, are larger, and fix more nitrogen in the absence and presence of inhibiting nitrate levels. At 25mM nitrate, nodules formed on pre-existing swelling sites induced by MtCEP1 overexpression. RNA interference-mediated silencing of several MtCEP genes revealed a negative correlation between transcript levels of MtCEP1 and MtCEP2 with the number of lateral roots. MtCEP1 peptide-dependent phenotypes were abolished or attenuated by altering or deleting key residues in its 15 amino acid domain. RNA-Seq analysis revealed that 89 and 116 genes were significantly up- and down-regulated, respectively, by MtCEP1 overexpression, including transcription factors WRKY, bZIP, ERF, and MYB, homologues of LOB29, SUPERROOT2, and BABY BOOM. Taken together, the data suggest that the MtCEP1 peptide modulates lateral root and nodule development in M. truncatula.

Structural and functional characterization of annexin 1 from Medicago truncatula

Annexins are calcium- and membrane-binding proteins that have been shown to have diverse properties such as actin, integrin and GTP binding, both in animals and plants. Recently, Medicago truncatula annexin 1 (AnnMt1) has been suggested to participate in nodulation (Nod factor signaling) and mycorrhization in legume plants. In this report we demonstrate for the first time that recombinant AnnMt1 (rec-AnnMt1) mediates membrane permeabilization to cations with conductance ranging from 16 pS to 329 pS. In agreement with other structurally determined annexins, homology modeling of AnnMt1 suggests that most of the functional determinants are found on the convex surface of the modeled structure. In conclusion, we propose a potential constitutive role of AnnMt1 in Nod factor signaling as a non-specific ion channel.

Legume adaptation to sulfur deficiency revealed by comparing nutrient allocation and seed traits in Medicago truncatula

Reductions in sulfur dioxide emissions and the use of sulfur-free mineral fertilizers are decreasing soil sulfur levels and threaten the adequate fertilization of most crops. To provide knowledge regarding legume adaptation to sulfur restriction, we subjected Medicago truncatula, a model legume species, to sulfur deficiency at various developmental stages, and compared the yield, nutrient allocation and seed traits. This comparative analysis revealed that sulfur deficiency at the mid-vegetative stage decreased yield and altered the allocation of nitrogen and carbon to seeds, leading to reduced levels of major oligosaccharides in mature seeds, whose germination was dramatically affected. In contrast, during the reproductive period, sulfur deficiency had little influence on yield and nutrient allocation, but the seeds germinated slowly and were characterized by low levels of a biotinylated protein, a putative indicator of germination vigor that has not been previously related to sulfur nutrition. Significantly, plants deprived of sulfur at an intermediary stage (flowering) adapted well by remobilizing nutrients from source organs to seeds, ensuring adequate quantities of carbon and nitrogen in seeds. This efficient remobilization of photosynthates may be explained by vacuolar sulfate efflux to maintain leaf metabolism throughout reproductive growth, as suggested by transcript and metabolite profiling. The seeds from these plants, deprived of sulfur at the floral transition, contained normal levels of major oligosaccharides but their germination was delayed, consistent with low levels of sucrose and the glycolytic enzymes required to restart seed metabolism during imbibition. Overall, our findings provide an integrative view of the legume response to sulfur deficiency.

The phytochrome gene family in soybean and a dominant negative effect of a soybean PHYA transgene on endogenous Arabidopsis PHYA

KEY MESSAGE

The evolutionary origin of the phytochrome genes in soybean was analyzed. The expression profiles of PHYA paralogs were characterized. The heterologous expression of GmPHYA1 in Arabidopsis resulted in longer hypocotyls. The phytochromes (PHY) are a small family of red/far-red light photoreceptors which regulate a number of important developmental responses in plants. So far, the members of the PHY gene family in soybean (Glycine max) remain unclear and an understanding of each member's physiological functions is limited. Our present in silico analysis revealed that the soybean genome harbors four PHYA, two PHYB and two PHYE, totally four pairs of eight PHY loci. The phylogenetic analysis suggested that the four PHY paralogous pairs originated from the latest round of genome duplication (~13 million years ago) and the four copies of PHYA were remnants of the two rounds of genome duplication (~58 and ~13 million years ago). A possible evolutionary history of PHYA homologs in the three legume species (soybean, Medicago truncatula, and Lotus japonicus) was proposed and the fate of duplicate soybean PHYA genes following polyploidization was discussed. The expression profiles of a soybean PHYA paralogous pair (GmPHYA1 and GmPHYA2) showed that the transcript abundance was highest in the aerial organs of young plants. The physiological role of GmPHYA1 was explored by observing the de-etiolation phenotype of transgenic Arabidopsis plants constitutively expressing GmPHYA1. The GmPHYA1 protein interfered with the function of endogenous PHYA with respect to de-etiolation in a dominant negative manner when exogenously expressed in Arabidopsis. The elucidation of the PHY gene family members in soybean provide us with a general description and understanding of the photoreceptor gene family in this important crop plant.

Patterns of Selection in Plant Genomes

Plants show a wide range of variation in mating system, ploidy level, and demographic history, allowing for unique opportunities to investigate the evolutionary and genetic factors affecting genome-wide patterns of positive and negative selection. In this review, we highlight recent progress in our understanding of the extent and nature of selection on plant genomes. We discuss differences in selection as they relate to variation in demography, recombination, mating system, and ploidy. We focus on the population genetic consequences of these factors and argue that, although variation in the magnitude of purifying selection is well documented, quantifying rates of positive selection and disentangling the relative importance of recombination, demography, and ploidy are ongoing challenges. Large-scale comparative studies that examine the relative and joint importance of these processes, combined with explicit models of population history and selection, are key and feasible goals for future work.

A tandem array of CBF/DREB1 genes is located in a major freezing tolerance QTL region on Medicago truncatula chromosome 6

BACKGROUND

Freezing provokes severe yield losses to different fall-sown annual legumes. Understanding the molecular bases of freezing tolerance is of great interest for breeding programs. Medicago truncatula Gaertn. is an annual temperate forage legume that has been chosen as a model species for agronomically and economically important legume crops. The present study aimed to identify positional candidate genes for a major freezing tolerance quantitative trait locus that was previously mapped to M. truncatula chromosome 6 (Mt-FTQTL6) using the LR3 population derived from a cross between the freezing-tolerant accession F83005-5 and the freezing-sensitive accession DZA045-5.

RESULTS

The confidence interval of Mt-FTQTL6 was narrowed down to the region comprised between markers MTIC153 and NT6054 using recombinant F7 and F8 lines. A bacterial-artificial chromosome (BAC) clone contig map was constructed in an attempt to close the residual assembly gap existing therein. Twenty positional candidate genes including twelve C-repeat binding factor (CBF)/dehydration-responsive element binding factor 1 (DREB1) genes were identified from BAC-derived sequences and whole-genome shotgun sequences (WGS). CBF/DREB1 genes are organized in a tandem array within an approximately 296-Kb region. Eleven CBF/DREB1 genes were isolated and sequenced from F83005-5 and DZA045-5 which revealed high polymorphism among these accessions. Unique features characterizing CBF/DREB1 genes from M. truncatula, such as alternative splicing and large tandem duplication, are elucidated for the first time.

CONCLUSIONS

Overall, twenty genes were identified as potential candidates to explain Mt-FTQTL6 effect. Their future functional characterization will uncover the gene(s) involved in freezing tolerance difference observed between F83005-5 and DZA045-5. Knowledge transfer for breeding improvement of crop legumes is expected. Furthermore, CBF/DREB1 related data will certainly have a large impact on research studies targeting this group of transcriptional activators in M. truncatula and other legume species.

Detecting small plant peptides using SPADA (Small Peptide Alignment Discovery Application)

BACKGROUND

Small peptides encoded as one- or two-exon genes in plants have recently been shown to affect multiple aspects of plant development, reproduction and defense responses. However, popular similarity search tools and gene prediction techniques generally fail to identify most members belonging to this class of genes. This is largely due to the high sequence divergence among family members and the limited availability of experimentally verified small peptides to use as training sets for homology search and ab initio prediction. Consequently, there is an urgent need for both experimental and computational studies in order to further advance the accurate prediction of small peptides.

RESULTS

We present here a homology-based gene prediction program to accurately predict small peptides at the genome level. Given a high-quality profile alignment, SPADA identifies and annotates nearly all family members in tested genomes with better performance than all general-purpose gene prediction programs surveyed. We find numerous mis-annotations in the current Arabidopsis thaliana and Medicago truncatula genome databases using SPADA, most of which have RNA-Seq expression support. We also show that SPADA works well on other classes of small secreted peptides in plants (e.g., self-incompatibility protein homologues) as well as non-secreted peptides outside the plant kingdom (e.g., the alpha-amanitin toxin gene family in the mushroom, Amanita bisporigera).

CONCLUSIONS

SPADA is a free software tool that accurately identifies and predicts the gene structure for short peptides with one or two exons. SPADA is able to incorporate information from profile alignments into the model prediction process and makes use of it to score different candidate models. SPADA achieves high sensitivity and specificity in predicting small plant peptides such as the cysteine-rich peptide families. A systematic application of SPADA to other classes of small peptides by research communities will greatly improve the genome annotation of different protein families in public genome databases.

PlantTFcat: an online plant transcription factor and transcriptional regulator categorization and analysis tool

Background

Plants regulate intrinsic gene expression through transcription factors (TFs), transcriptional regulators (TRs), chromatin regulators (CRs), and the basal transcription machinery. An understanding of plant gene regulatory mechanisms at a systems level requires the identification of these regulatory elements on a genomic scale.

Results

Here, we present PlantTFcat, a high-performance web-based analysis tool that is designed to identify and categorize plant TF/TR/CR genes from genome-scale protein and nucleic acid sequences by systematically analyzing InterProScan domain patterns in protein sequences. The comprehensive prediction logics that are included in PlantTFcat are based on relationships between gene families and conserved domains from 108 published plant TF/TR/CR families. These prediction logics effectively distinguish TF/TR/CR families with common conserved domains. Our systematic performance evaluations indicate that PlantTFcat annotates known TF/TR/CR families with high coverage and sensitivity.

Conclusions

PlantTFcat provides an analysis tool to identify and categorize plant TF/TR/CR genes on a genomic scale. PlantTFcat is freely available to the public at http://plantgrn.noble.org/PlantTFcat/.