A comprehensive draft genome sequence for lupin (Lupinus angustifolius), an emerging health food: insights into plant-microbe interactions and legume evolution.

Lupins are important grain legume crops that form a critical part of sustainable farming systems, reducing fertilizer use and providing disease breaks. It has a basal phylogenetic position relative to other crop and model legumes and a high speciation rate. Narrow-leafed lupin (NLL; Lupinus angustifolius L.) is gaining popularity as a health food, which is high in protein and dietary fibre but low in starch and gluten-free. We report the draft genome assembly (609 Mb) of NLL cultivar Tanjil, which has captured >98% of the gene content, sequences of additional lines and a dense genetic map. Lupins are unique among legumes and differ from most other land plants in that they do not form mycorrhizal associations. Remarkably, we find that NLL has lost all mycorrhiza-specific genes, but has retained genes commonly required for mycorrhization and nodulation. In addition, the genome also provided candidate genes for key disease resistance and domestication traits. We also find evidence of a whole-genome triplication at around 25 million years ago in the genistoid lineage leading to Lupinus. Our results will support detailed studies of legume evolution and accelerate lupin breeding programmes.

Transcriptome sequencing of different narrow-leafed lupin tissue types provides a comprehensive uni-gene assembly and extensive gene-based molecular markers.

Narrow-leafed lupin (NLL; Lupinus angustifolius L.) is an important grain legume crop that is valuable for sustainable farming and is becoming recognized as a human health food. NLL breeding is directed at improving grain production, disease resistance, drought tolerance and health benefits. However, genetic and genomic studies have been hindered by a lack of extensive genomic resources for the species. Here, the generation, de novo assembly and annotation of transcriptome datasets derived from five different NLL tissue types of the reference accession cv. Tanjil are described. The Tanjil transcriptome was compared to transcriptomes of an early domesticated cv. Unicrop, a wild accession P27255, as well as accession 83A:476, together being the founding parents of two recombinant inbred line (RIL) populations. In silico predictions for transcriptome-derived gene-based length and SNP polymorphic markers were conducted and corroborated using a survey assembly sequence for NLL cv. Tanjil. This yielded extensive indel and SNP polymorphic markers for the two RIL populations. A total of 335 transcriptome-derived markers and 66 BAC-end sequence-derived markers were evaluated, and 275 polymorphic markers were selected to genotype the reference NLL 83A:476 × P27255 RIL population. This significantly improved the completeness, marker density and quality of the reference NLL genetic map.

Comparative proteomic analysis of genotypic variation in germination and early seedling growth of chickpea under suboptimal soil-water conditions.

Protein expression patterns in imbibed seeds of three cultivars of chickpea (Cicer arietinum L.) with different rates of germination under limiting water supply in soil (>10% water holding capacity) were compared. A large number of soluble proteins expressed earlier and at higher levels in cv Rupali seeds compared to two other genotypes that germinated less rapidly (KH850) or not at all (KJ850). Among the proteins identified were those with chaperone-like functions, including LEA and HSP proteins and proteins involved in metabolism of reactive oxygen species (ROS). Only NAD-malate dehydrogenase was identified as an early, differentially abundant enzyme of the TCA cycle, but in cv Rupali, expression of phospho-enol-pyruvate carboxykinase rose very rapidly to a high level, indicating that an anaplerotic C input to the TCA cycle may have been important. Proteinase inhibitors were more highly expressed in the genotype that did not germinate compared to cv Rupali. Clustering analysis of proteomic data indicated a link between groups of proteins, implying a common regulatory mechanism possibly at the transcriptional level. The chaperone-like proteins and enzymes of ROS homeostasis provide a useful starting point for molecular genetic analysis that may well identify other important genes for the early germination trait.

Development of an assay to evaluate differences in germination rate among chickpea genotypes under limited water content.

An assay system that provides rapid and reproducible germination under low soil water content (<10% water holding capacity (WHC)) was developed and used to compare how chickpea (Cicer arietinum L.) genotypes complete germination, without the technical difficulties of accurately controlling water levels. The system consisted of small plastic containers (50mm×50mm×60mm) filled with river sand and tightly closed (but not sealed) to minimise water loss and maintain constant soil water content during germination. Seed size influenced germination performance at low WHC. Small seeds within a single genotype germinated successfully and entered into the early stages of seedling growth, but germination of large seeds was inhibited, failing to germinate at 5% WHC. Small seeds were more efficient in remobilising seed reserves to seedling tissues than larger seeds. Under optimal WHC, the germination rate and subsequent radicle growth was similar among genotypes but at low WHC, there was variation despite seeds being of comparable size and imbibing equally. This suggests that the physiological threshold of threshold water potential for initiation of germination reflects genotypic differences. The assay system provides a suitable experimental tool to examine gene expression in contrasting genotypes during germination and early stages of seedling growth with a view to identifying the genes involved in superior performance under water limited field conditions.

Development of genomic resources for the narrow-leafed lupin (Lupinus angustifolius): construction of a bacterial artificial chromosome (BAC) library and BAC-end sequencing.

BACKGROUND: Lupinus angustifolius L, also known as narrow-leafed lupin (NLL), is becoming an important grain legume crop that is valuable for sustainable farming and is becoming recognised as a potential human health food. Recent interest is being directed at NLL to improve grain production, disease and pest management and health benefits of the grain. However, studies have been hindered by a lack of extensive genomic resources for the species. RESULTS: A NLL BAC library was constructed consisting of 111,360 clones with an average insert size of 99.7 Kbp from cv Tanjil. The library has approximately 12 × genome coverage. Both ends of 9600 randomly selected BAC clones were sequenced to generate 13985 BAC end-sequences (BESs), covering approximately 1% of the NLL genome. These BESs permitted a preliminary characterisation of the NLL genome such as organisation and composition, with the BESs having approximately 39% G:C content, 16.6% repetitive DNA and 5.4% putative gene-encoding regions. From the BESs 9966 simple sequence repeat (SSR) motifs were identified and some of these are shown to be potential markers. CONCLUSIONS: The NLL BAC library and BAC-end sequences are powerful resources for genetic and genomic research on lupin. These resources will provide a robust platform for future high-resolution mapping, map-based cloning, comparative genomics and assembly of whole-genome sequencing data for the species.

Identification and characterisation of seed storage protein transcripts from Lupinus angustifolius.

BACKGROUND: In legumes, seed storage proteins are important for the developing seedling and are an important source of protein for humans and animals. Lupinus angustifolius (L.), also known as narrow-leaf lupin (NLL) is a grain legume crop that is gaining recognition as a potential human health food as the grain is high in protein and dietary fibre, gluten-free and low in fat and starch. RESULTS: Genes encoding the seed storage proteins of NLL were characterised by sequencing cDNA clones derived from developing seeds. Four families of seed storage proteins were identified and comprised three unique α, seven ß, two γ and four δ conglutins. This study added eleven new expressed storage protein genes for the species. A comparison of the deduced amino acid sequences of NLL conglutins with those available for the storage proteins of Lupinus albus (L.), Pisum sativum (L.), Medicago truncatula (L.), Arachis hypogaea (L.) and Glycine max (L.) permitted the analysis of a phylogenetic relationships between proteins and demonstrated, in general, that the strongest conservation occurred within species. In the case of 7S globulin (ß conglutins) and 2S sulphur-rich albumin (δ conglutins), the analysis suggests that gene duplication occurred after legume speciation. This contrasted with 11S globulin (α conglutin) and basic 7S (γ conglutin) sequences where some of these sequences appear to have diverged prior to speciation. The most abundant NLL conglutin family was ß (56%), followed by α (24%), δ (15%) and γ (6%) and the transcript levels of these genes increased 103 to 106 fold during seed development. We used the 16 NLL conglutin sequences identified here to determine that for individuals specifically allergic to lupin, all seven members of the ß conglutin family were potential allergens. CONCLUSION: This study has characterised 16 seed storage protein genes in NLL including 11 newly-identified members. It has helped lay the foundation for efforts to use molecular breeding approaches to improve lupins, for example by reducing allergens or increasing the expression of specific seed storage protein(s) with desirable nutritional properties.

Macromolecular composition of phloem exudate from white lupin (Lupinus albus L.).

BACKGROUND: Members of the legume genus Lupinus exude phloem 'spontaneously' from incisions made to the vasculature. This feature was exploited to document macromolecules present in exudate of white lupin (Lupinus albus [L.] cv Kiev mutant), in particular to identify proteins and RNA molecules, including microRNA (miRNA). RESULTS: Proteomic analysis tentatively identified 86 proteins from 130 spots collected from 2D gels analysed by partial amino acid sequence determination using MS/MS. Analysis of a cDNA library constructed from exudate identified 609 unique transcripts. Both proteins and transcripts were classified into functional groups. The largest group of proteins comprised those involved in metabolism (24%), followed by protein modification/turnover (9%), redox regulation (8%), cell structural components (6%), stress and defence response (6%) with fewer in other groups. More prominent proteins were cyclophilin, ubiquitin, a glycine-rich RNA-binding protein, a group of proteins that comprise a glutathione/ascorbate-based mechanism to scavenge oxygen radicals, enzymes of glycolysis and other metabolism including methionine and ethylene synthesis. Potential signalling macromolecules such as transcripts encoding proteins mediating calcium level and the Flowering locus T (FT) protein were also identified. From around 330 small RNA clones (18-25 nt) 12 were identified as probable miRNAs by homology with those from other species. miRNA composition of exudate varied with site of collection (e.g. upward versus downward translocation streams) and nutrition (e.g. phosphorus level). CONCLUSIONS: This is the first inventory of macromolecule composition of phloem exudate from a species in the Fabaceae, providing a basis to identify systemic signalling macromolecules with potential roles in regulating development, growth and stress response of legumes.

Consequences of transforming narrow leafed lupin (Lupinus angustifolius [L.]) with an ipt gene under control of a flower-specific promoter.

Phenotypes of five transgenic lines of narrow-leafed lupin (Lupinus angustifolius [L] cv Merrit) stably transformed with the isopentenyl pyrophosphate transferase (ipt) gene from Agrobacterium tumefaciens coupled to a flower-specific promoter (TP12) from Nicotiana tabacum [L.] are described. Expression of the transgene was detected in floral tissues and in shoot apical meristems on all orders of inflorescence. In each transgenic line there was significant axillary bud outgrowth at all nodes on the main stem with pronounced branch development from the more basal nodes in three of the lines. The lowest basal branches developed in a manner similar to the upper stem axillary branches on cv Merrit and bore fruits, which, in two lines, contained a significant yield of filled seeds at maturity. Senescence of the cotyledons was delayed in all lines with green cotyledons persisting beyond anthesis in one case. IPT expression increased cytokinin (CK) levels in flowers, meristem tissues and phloem exudates in a form specific manner, which was suggestive of localized flower and meristem production with significant long-distance re-distribution in phloem. The total number of fruits formed (pod set) on some transgenic lines was increased compared to cv Merrit. Grain size compared to cv Merrit was not significantly altered in transgenic lines.

Macromolecules in phloem exudates--a review.

Proteomic and transcriptomic analyses using the growing resources of genomic information have been applied to identification of macromolecules in exudates collected from phloem. Most of the analyses rely on collection of exudate following incisions made to the vasculature, but some limited data are available for exudates collected from excised aphid stylets. Species examined, to date, include a number of cereals (rice, barley, and wheat), a number of cucurbits, castor bean, members of the genus Lupinus, brassicas, and Arabidopsis. As many as 1,100 proteins, some hundreds of transcripts, and a growing number of small ribonucleic acids (RNAs), including micro-RNAs, have been identified across the species with a high degree of commonality. Questions relating to the nature and extent of contamination of sieve element contents with those of surrounding companion cells and nonvascular cells are addressed together with likely functions of identified macromolecules. The review considers likely translocation and systemic signaling functions among the macromolecular inventory of phloem exudates.

Uric acid deposits in symbiotic marine algae.

The symbiosis between cnidarians and dinoflagellate algae is not understood at the cell or molecular level, yet this relationship is responsible for the formation of thousands of square kilometres of coral reefs. We have investigated the nature of crystalline material prominent within marine algal symbionts of Aiptasia sp. anemones. This material, which has historically been considered to be calcium oxalate, is shown to be uric acid. We demonstrate that these abundant uric acid stores can be mobilized rapidly, thereby allowing the algal symbionts to flourish in an otherwise N-poor environment. This is the first report of uric acid accumulation by symbiotic marine algae. These data provide new insight and considerations for understanding the physiological basis of algal symbioses, and represent a new and previously unconsidered aspect of N metabolism in cnidarian, and a variety of other, marine symbioses.

Synthesis, transport and accumulation of quinolizidine alkaloids in Lupinus albus L. and L. angustifolius L.

Each of the principal quinolizidine alkaloids (QA) found in both xylem and phloem exudates together with extracts from all component organs collected from bitter (cv. Lupini) and sweet (cv. Ultra) cultivars of Lupinus albus L. were quantified by gas chromatographic analyses throughout reproductive development. In addition to establishing the major translocated QA species estimates for fluxes of QA to developing fruits based on their sap composition and water economy showed that around half of the QA that accumulated in fruit tissues was due to synthesis in situ and half to translocation principally by phloem. Detailed analyses of QA in transport fluids and component organs were extended to reciprocal homo- and hetero-grafts using bitter (cv. Fest) and sweet (cv. Danja) cultivars of L. angustifolius L. These data confirmed that the majority of QA were synthesized in shoot tissues. In both lupin species feeding and analysis of deuterated QA (lupanine and 13-hydroxylupanine) were used as tracers to demonstrate direct redistribution of alkaloids by translocation from mature leaves in phloem.

Alternative splicing of the Vupur3 transcript in cowpea produces multiple mRNA species with a single protein product that is present in both plastids and mitochondria.

A heterogeneous population of cDNAs (designated Vupur3) encoding phosphoribosylglycinamide formyltransferase (GART; EC 2.1.2.2) was isolated from a cowpea (Vigna unguiculata L. Walp.) nodule library. Three classes of cDNA with the same ORF, but differing in their 3'-UTRs, were identified. Southern analysis and sequencing of genomic DNA confirmed that these differences result from alternative splicing of the primary transcript of a single Vupur3 gene. Alternative splicing does not appear to play a role in the production of soybean (Glycine max Merrill.) pur3 transcripts. The presence of the protein product of the Vupur3 gene, GART, in plastids and mitochondria was confirmed by immunoblotting with antibodies raised against the recombinant protein. The antibodies recognised two proteins with apparent molecular masses of 27 and 27.5 kDa in both mitochondria and plastids. All Vupur3 transcripts have two in-frame start codons that are active in wheatgerm in vitro transcription / translation experiments suggesting a mechanism by which the gene product could be targeted to two organelles. Like other genes encoding enzymes for purine synthesis, Vupur3 is expressed in nodules before nitrogen fixation begins but in contrast to these genes its expression does not increase markedly after nitrogen fixation begins.

Adenylate gradients and Ar:O2 effects on legume nodules. II. Changes in the subcellular adenylate pools.

Central infected zone tissue of soybean (Glycine max L. Merr.) nodules was fractionated into separate subcellular compartments using density gradient centrifugation in nonaqueous solvents to better understand how exposure to Ar:O(2) (80:20%, v/v) atmosphere affects C and N metabolism, and to explore a potential role for adenylates in regulating O(2) diffusion. When nodules were switched from air to Ar:O(2), adenylate energy charge (AEC) in the plant cytosol rose from 0.63 +/- 0.02 to 0.73 +/- 0.02 within 7 min and to 0.80 +/- 0.01 by 60 min. In contrast, AEC of the mitochondrial compartment of this central zone tissue remained high (0.80 +/- 0.02 to 0.81 +/- 0.02) following Ar treatment while that of the bacteroid compartment was unchanged, at 0.73 +/- 0.02, after 7 min, but declined to 0.57 +/- 0.03 after 60 min. These results were consistent with a simulation model that predicted Ar:O(2) exposure would first reduce ATP demand for ammonia assimilation and rapidly increase cytosolic AEC, before the Ar:O(2)-induced decline mediated by a decrease in nodule O(2) permeability reduces bacteroid AEC. The possibility that adenylates play a key, integrating role in regulating nodule permeability to oxygen diffusion is discussed.

Adenylate gradients and Ar:O(2) effects on legume nodules: I. Mathematical models.

Mathematical models were developed to test the likelihood that large cytosolic adenylate concentration gradients exist across the bacteria-infected cells of legume nodules. Previous studies hypothesized that this may be the case to account for the unusually low adenylate energy charge (AEC; 0.65) measured in the plant fraction of metabolically active nodules (M.M. Kuzma, H. Winter, P. Storer, I. Oresnik, C.A. Atkins, D.B. Layzell [1999] Plant Physiol 119: 399-407). Simulations coupled leghemoglobin-facilitated O(2) diffusion into the infected cell, through bacteroid nitrogenase activity, with the ATP demand for transport and ammonia assimilation in the plant fraction of ureide- and amide-producing nodules. Although large cytosolic adenylate gradients were predicted to exist in both nodule types, amide nodules were predicted to have steeper AEC gradients (0.82-0.52) than ureide nodules (0.82-0.61). The differences were attributed to an additional ATP demand for Asn synthesis in the amide nodule. Simulations for nodules transferred to an Ar:O(2) atmosphere predicted a major reduction in the magnitude of adenylate gradients and an increase in the AEC of the plant fraction. Results were consistent with a number of experimental studies and were used to propose an experimental test of the models.

Differential expression of a novel gene during seed triacylglycerol accumulation in lupin species ( Lupinus angustifolius L. and L. mutabilis L.).

Seed triacylglycerols (TAGs) are stored as energy reserves and extracted for various end-product uses. In lupins, seed oil content varies from 16% in Lupinus mutabilisto 8% in L. angustifolius. We have shown that TAGs rapidly accumulate during mid-stages of seed development in L. mutabilis compared to the lower seed oil species, L. angustifolius. In this study, we have targeted the key enzymes of the lipid biosynthetic pathway, acetyl-CoA carboxylase (ACCase) and diacylglycerol acyltransferase (DAGAT), to determine factors regulating TAG accumulation between two lupin species. A twofold increase in ACCase activity was observed in L. mutabilis relative to L. angustifolius and correlated with rapid TAG accumulation. No difference in DAGAT activity was detected. We have identified, cloned and partially characterised a novel gene differentially expressed during TAG accumulation between L. angustifolius and L. mutabilis. The gene has some identity to the glucose dehydrogenase family previously described in barley and bacteria and the significance of its expression levels during seed development in relation to TAG accumulation is discussed. DNA sequence analysis of the promoter in both L. angustifolius and L. mutabilis identified putative matrix attachment regions and recognition sequences for transcription binding sites similar to those found in the Adh1 gene from Arabidopsis. The identical promoter regions between species indicate that differential gene expression is controlled by alternative transcription factors, accessibility to binding sites or a combination of both.

Effect of short-term N(2) deficiency on expression of the ureide pathway in cowpea root nodules.

Root systems of 28-d-old cowpea (Vigna unguiculata L. Walp cv Vita 3: Bradyrhizobium sp. strain CB756) plants bearing nitrogen-fixing nodules in sand culture were exposed to an atmosphere of Ar:O(2) (80:20, v/v) for 48 h and then returned to air. Root systems of control plants were maintained in air throughout. Nodules were harvested at the same times in control and Ar:O(2)-treated root systems. Activities of two enzymes of de novo purine synthesis, glycinamide ribonucleotide transformylase (GART; EC 2.1.2.2), aminoimidazole ribonucleotide synthetase (AIRS; EC 6.3.3.1), uricase (EC 1.7.3.3), and phosphoenolpyruvate carboxylase (PEPC; EC 4.1.1.31) were measured together with the protein level of each using immune-specific polyclonal antibodies. AIRS activity and protein both declined to very low levels within 6 h in Ar:O(2) together with a decline in transcript level of pur5, the encoding gene. GART activity, protein, and transcript (pur3) levels were relatively stable. Uricase activity declined in Ar:O(2) as rapidly as AIRS activity but the protein was stable. PEPC activity showed evidence of increased sensitivity to inhibition by malate but the protein level was stable. The data indicate that the flux of fixed N from bacteroids (N(2)-fixing nodule bacteria) is in some way associated with transcriptional control over pur5 and possibly also catabolism of AIRS protein. In contrast, there is limited posttranslational control over GART and PEPC and close posttranslational control over uricase activity. The significance of these different levels of regulation is discussed in relation to the overall control of enhanced expression of plant enzymes in the cowpea symbiosis.