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

BACKGROUND: Alfalfa (Medicago sativa L.) is the primary forage legume crop species in the United States and plays essential economic and ecological roles in agricultural systems across the country. Modern alfalfa is the result of hybridization between tetraploid M. sativa ssp. sativa and M. sativa ssp. falcata. Due to its large and complex genome, there are few genomic resources available for alfalfa improvement. RESULTS: A de novo transcriptome assembly from two alfalfa subspecies, M. sativa ssp. sativa (B47) and M. sativa ssp. falcata (F56) was developed using Illumina RNA-seq technology. Transcripts from roots, nitrogen-fixing root nodules, leaves, flowers, elongating stem internodes, and post-elongation stem internodes were assembled into the Medicago sativa Gene Index 1.2 (MSGI 1.2) representing 112,626 unique transcript sequences. Nodule-specific and transcripts involved in cell wall biosynthesis were identified. Statistical analyses identified 20,447 transcripts differentially expressed between the two subspecies. Pair-wise comparisons of each tissue combination identified 58,932 sequences differentially expressed in B47 and 69,143 sequences differentially expressed in F56. Comparing transcript abundance in floral tissues of B47 and F56 identified expression differences in sequences involved in anthocyanin and carotenoid synthesis, which determine flower pigmentation. Single nucleotide polymorphisms (SNPs) unique to each M. sativa subspecies (110,241) were identified. CONCLUSIONS: The Medicago sativa Gene Index 1.2 increases the expressed sequence data available for alfalfa by ninefold and can be expanded as additional experiments are performed. The MSGI 1.2 transcriptome sequences, annotations, expression profiles, and SNPs were assembled into the Alfalfa Gene Index and Expression Database (AGED) at http://plantgrn.noble.org/AGED/ , a publicly available genomic resource for alfalfa improvement and legume research.

Transgene silencing of sucrose synthase in alfalfa (Medicago sativa L.) stem vascular tissue suggests a role for invertase in cell wall cellulose synthesis.

BACKGROUND: Alfalfa (Medicago sativa L.) is a widely adapted perennial forage crop that has high biomass production potential. Enhanced cellulose content in alfalfa stems would increase the value of the crop as a bioenergy feedstock. We examined if increased expression of sucrose synthase (SUS; EC 2.4.1.13) would increase cellulose in stem cell walls. RESULTS: Alfalfa plants were transformed with a truncated alfalfa phosphoenolpyruvate carboxylase gene promoter (PEPC7-P4) fused to an alfalfa nodule-enhanced SUS cDNA (MsSUS1) or the ß-glucuronidase (GUS) gene. Strong GUS expression was detected in xylem and phloem indicating that the PEPC7-P4 promoter was active in stem vascular tissue. In contrast to expectations, MsSUS1 transcript accumulation was reduced 75-90 % in alfalfa plants containing the PEPC7-P4::MsSUS1 transgene compared to controls. Enzyme assays indicated that SUS activity in stems of selected down-regulated transformants was reduced by greater than 95 % compared to the controls. Although SUS activity was detected in xylem and phloem of control plants by in situ enzyme assays, plants with the PEPC7-P4::MsSUS1 transgene lacked detectable SUS activity in post-elongation stem (PES) internodes and had very low SUS activity in elongating stem (ES) internodes. Loss of SUS protein in PES internodes of down-regulated lines was confirmed by immunoblots. Down-regulation of SUS expression and activity in stem tissue resulted in no obvious phenotype or significant change in cell wall sugar composition. However, alkaline/neutral (A/N) invertase activity increased in SUS down-regulated lines and high levels of acid invertase activity were observed. In situ enzyme assays of stem tissue showed localization of neutral invertase in vascular tissues of ES and PES internodes. CONCLUSIONS: These results suggest that invertases play a primary role in providing glucose for cellulose biosynthesis or compensate for the loss of SUS1 activity in stem vascular tissue.

An RNA-Seq based gene expression atlas of the common bean.

BACKGROUND: Common bean (Phaseolus vulgaris) is grown throughout the world and comprises roughly 50% of the grain legumes consumed worldwide. Despite this, genetic resources for common beans have been lacking. Next generation sequencing, has facilitated our investigation of the gene expression profiles associated with biologically important traits in common bean. An increased understanding of gene expression in common bean will improve our understanding of gene expression patterns in other legume species. RESULTS: Combining recently developed genomic resources for Phaseolus vulgaris, including predicted gene calls, with RNA-Seq technology, we measured the gene expression patterns from 24 samples collected from seven tissues at developmentally important stages and from three nitrogen treatments. Gene expression patterns throughout the plant were analyzed to better understand changes due to nodulation, seed development, and nitrogen utilization. We have identified 11,010 genes differentially expressed with a fold change ≥ 2 and a P-value < 0.05 between different tissues at the same time point, 15,752 genes differentially expressed within a tissue due to changes in development, and 2,315 genes expressed only in a single tissue. These analyses identified 2,970 genes with expression patterns that appear to be directly dependent on the source of available nitrogen. Finally, we have assembled this data in a publicly available database, The Phaseolus vulgaris Gene Expression Atlas (Pv GEA), http://plantgrn.noble.org/PvGEA/ . Using the website, researchers can query gene expression profiles of their gene of interest, search for genes expressed in different tissues, or download the dataset in a tabular form. CONCLUSIONS: These data provide the basis for a gene expression atlas, which will facilitate functional genomic studies in common bean. Analysis of this dataset has identified genes important in regulating seed composition and has increased our understanding of nodulation and impact of the nitrogen source on assimilation and distribution throughout the plant.

An RNA-Seq transcriptome analysis of orthophosphate-deficient white lupin reveals novel insights into phosphorus acclimation in plants.

Phosphorus, in its orthophosphate form (P(i)), is one of the most limiting macronutrients in soils for plant growth and development. However, the whole-genome molecular mechanisms contributing to plant acclimation to P(i) deficiency remain largely unknown. White lupin (Lupinus albus) has evolved unique adaptations for growth in P(i)-deficient soils, including the development of cluster roots to increase root surface area. In this study, we utilized RNA-Seq technology to assess global gene expression in white lupin cluster roots, normal roots, and leaves in response to P(i) supply. We de novo assembled 277,224,180 Illumina reads from 12 complementary DNA libraries to build what is to our knowledge the first white lupin gene index (LAGI 1.0). This index contains 125,821 unique sequences with an average length of 1,155 bp. Of these sequences, 50,734 were transcriptionally active (reads per kilobase per million reads ≥ 3), representing approximately 7.8% of the white lupin genome, using the predicted genome size of Lupinus angustifolius as a reference. We identified a total of 2,128 sequences differentially expressed in response to P(i) deficiency with a 2-fold or greater change and P ≤ 0.05. Twelve sequences were consistently differentially expressed due to P(i) deficiency stress in three species, Arabidopsis (Arabidopsis thaliana), potato (Solanum tuberosum), and white lupin, making them ideal candidates to monitor the P(i) status of plants. Additionally, classic physiological experiments were coupled with RNA-Seq data to examine the role of cytokinin and gibberellic acid in P(i) deficiency-induced cluster root development. This global gene expression analysis provides new insights into the biochemical and molecular mechanisms involved in the acclimation to P(i) deficiency.

Legume genomics: understanding biology through DNA and RNA sequencing.

BACKGROUND: The legume family (Leguminosae) consists of approx. 17 000 species. A few of these species, including, but not limited to, Phaseolus vulgaris, Cicer arietinum and Cajanus cajan, are important dietary components, providing protein for approx. 300 million people worldwide. Additional species, including soybean (Glycine max) and alfalfa (Medicago sativa), are important crops utilized mainly in animal feed. In addition, legumes are important contributors to biological nitrogen, forming symbiotic relationships with rhizobia to fix atmospheric N2 and providing up to 30 % of available nitrogen for the next season of crops. The application of high-throughput genomic technologies including genome sequencing projects, genome re-sequencing (DNA-seq) and transcriptome sequencing (RNA-seq) by the legume research community has provided major insights into genome evolution, genomic architecture and domestication. SCOPE AND CONCLUSIONS: This review presents an overview of the current state of legume genomics and explores the role that next-generation sequencing technologies play in advancing legume genomics. The adoption of next-generation sequencing and implementation of associated bioinformatic tools has allowed researchers to turn each species of interest into their own model organism. To illustrate the power of next-generation sequencing, an in-depth overview of the transcriptomes of both soybean and white lupin (Lupinus albus) is provided. The soybean transcriptome focuses on analysing seed development in two near-isogenic lines, examining the role of transporters, oil biosynthesis and nitrogen utilization. The white lupin transcriptome analysis examines how phosphate deficiency alters gene expression patterns, inducing the formation of cluster roots. Such studies illustrate the power of next-generation sequencing and bioinformatic analyses in elucidating the gene networks underlying biological processes.

Comparison of stem morphology and anatomy of two alfalfa clonal lines exhibiting divergent cell wall composition.

BACKGROUND: In previous research, two alfalfa clonal lines (252 and 1283) were identified that exhibited environmentally stable differences in stem cell walls. Compared with stems of 1283, stems of 252 have a higher cell wall concentration and greater amounts of lignin and cellulose but reduced levels of pectic sugar residues. These results suggest greater deposition of secondary xylem and a reduction in pith in stems of 252 compared with 1283. RESULTS: The stem morphology and anatomy of first-cut and second-cut harvests of field-grown 1283 and 252 were examined. For both harvests, stems of 1283 were thicker and had a higher leaf/stem ratio compared with stems of 252. Stem cross-sections of both genotypes were stained for lignin, and the proportions of stem area that were pith and secondary xylem were measured using ImageJ. Stems of 252 exhibited greater deposition of secondary xylem and a reduction in pith proportion compared with stems of 1283 for the first-cut harvest, but this difference was not statistically significant for the second-cut harvest. CONCLUSION: The results indicate that the proportions of secondary xylem and pith are not environmentally stable in these two genotypes and hence cannot be the sole basis for the differences in cell wall concentration/composition.

White lupin cluster root acclimation to phosphorus deficiency and root hair development involve unique glycerophosphodiester phosphodiesterases.

White lupin (Lupinus albus) is a legume that is very efficient in accessing unavailable phosphorus (Pi). It develops short, densely clustered tertiary lateral roots (cluster/proteoid roots) in response to Pi limitation. In this report, we characterize two glycerophosphodiester phosphodiesterase (GPX-PDE) genes (GPX-PDE1 and GPX-PDE2) from white lupin and propose a role for these two GPX-PDEs in root hair growth and development and in a Pi stress-induced phospholipid degradation pathway in cluster roots. Both GPX-PDE1 and GPX-PDE2 are highly expressed in Pi-deficient cluster roots, particularly in root hairs, epidermal cells, and vascular bundles. Expression of both genes is a function of both Pi availability and photosynthate. GPX-PDE1 Pi deficiency-induced expression is attenuated as photosynthate is deprived, while that of GPX-PDE2 is strikingly enhanced. Yeast complementation assays and in vitro enzyme assays revealed that GPX-PDE1 shows catalytic activity with glycerophosphocholine while GPX-PDE2 shows highest activity with glycerophosphoinositol. Cell-free protein extracts from Pi-deficient cluster roots display GPX-PDE enzyme activity for both glycerophosphocholine and glycerophosphoinositol. Knockdown of expression of GPX-PDE through RNA interference resulted in impaired root hair development and density. We propose that white lupin GPX-PDE1 and GPX-PDE2 are involved in the acclimation to Pi limitation by enhancing glycerophosphodiester degradation and mediating root hair development.

Objective Phenotyping of Root System Architecture Using Image Augmentation and Machine Learning in Alfalfa (Medicago sativa L.).

Active breeding programs specifically for root system architecture (RSA) phenotypes remain rare; however, breeding for branch and taproot types in the perennial crop alfalfa is ongoing. Phenotyping in this and other crops for active RSA breeding has mostly used visual scoring of specific traits or subjective classification into different root types. While image-based methods have been developed, translation to applied breeding is limited. This research is aimed at developing and comparing image-based RSA phenotyping methods using machine and deep learning algorithms for objective classification of 617 root images from mature alfalfa plants collected from the field to support the ongoing breeding efforts. Our results show that unsupervised machine learning tends to incorrectly classify roots into a normal distribution with most lines predicted as the intermediate root type. Encouragingly, random forest and TensorFlow-based neural networks can classify the root types into branch-type, taproot-type, and an intermediate taproot-branch type with 86% accuracy. With image augmentation, the prediction accuracy was improved to 97%. Coupling the predicted root type with its prediction probability will give breeders a confidence level for better decisions to advance the best and exclude the worst lines from their breeding program. This machine and deep learning approach enables accurate classification of the RSA phenotypes for genomic breeding of climate-resilient alfalfa.

An integrative approach to genomic introgression mapping.

Near-isogenic lines (NILs) are valuable genetic resources for many crop species, including soybean (Glycine max). The development of new molecular platforms promises to accelerate the mapping of genetic introgressions in these materials. Here, we compare some existing and emerging methodologies for genetic introgression mapping: single-feature polymorphism analysis, Illumina GoldenGate single nucleotide polymorphism (SNP) genotyping, and de novo SNP discovery via RNA-Seq analysis of next-generation sequence data. We used these methods to map the introgressed regions in an iron-inefficient soybean NIL and found that the three mapping approaches are complementary when utilized in combination. The comparative RNA-Seq approach offers several additional advantages, including the greatest mapping resolution, marker depth, and de novo marker utility for downstream fine-mapping analysis. We applied the comparative RNA-Seq method to map genetic introgressions in an additional pair of NILs exhibiting differential seed protein content. Furthermore, we attempted to optimize the comparative RNA-Seq approach by assessing the impact of sequence depth, SNP identification methodology, and post hoc analyses on SNP discovery rates. We conclude that the comparative RNA-Seq approach can be optimized with sufficient sampling and by utilizing a post hoc correction accounting for gene density variation that controls for false discoveries.

Phenotyping seedlings for selection of root system architecture in alfalfa (Medicago sativa L.).

BACKGROUND: The root system architecture (RSA) of alfalfa (Medicago sativa L.) affects biomass production by influencing water and nutrient uptake, including nitrogen fixation. Further, roots are important for storing carbohydrates that are needed for regrowth in spring and after each harvest. Previous selection for a greater number of branched and fibrous roots significantly increased alfalfa biomass yield. However, phenotyping root systems of mature alfalfa plant is labor-intensive, time-consuming, and subject to environmental variability and human error. High-throughput and detailed phenotyping methods are needed to accelerate the development of alfalfa germplasm with distinct RSAs adapted to specific environmental conditions and for enhancing productivity in elite germplasm. In this study methods were developed for phenotyping 14-day-old alfalfa seedlings to identify measurable root traits that are highly heritable and can differentiate plants with either a branched or a tap rooted phenotype. Plants were grown in a soil-free mixture under controlled conditions, then the root systems were imaged with a flatbed scanner and measured using WinRhizo software. RESULTS: The branched root plants had a significantly greater number of tertiary roots and significantly longer tertiary roots relative to the tap rooted plants. Additionally, the branch rooted population had significantly more secondary roots > 2.5 cm relative to the tap rooted population. These two parameters distinguishing phenotypes were confirmed using two machine learning algorithms, Random Forest and Gradient Boosting Machines. Plants selected as seedlings for the branch rooted or tap rooted phenotypes were used in crossing blocks that resulted in a genetic gain of 10%, consistent with the previous selection strategy that utilized manual root scoring to phenotype 22-week-old-plants. Heritability analysis of various root architecture parameters from selected seedlings showed tertiary root length and number are highly heritable with values of 0.74 and 0.79, respectively. CONCLUSIONS: The results show that seedling root phenotyping is a reliable tool that can be used for alfalfa germplasm selection and breeding. Phenotypic selection of RSA in seedlings reduced time for selection by 20 weeks, significantly accelerating the breeding cycle.

Transcript profiling of common bean (Phaseolus vulgaris L.) using the GeneChip Soybean Genome Array: optimizing analysis by masking biased probes.

BACKGROUND: Common bean (Phaseolus vulgaris L.) and soybean (Glycine max) both belong to the Phaseoleae tribe and share significant coding sequence homology. This suggests that the GeneChip(R) Soybean Genome Array (soybean GeneChip) may be used for gene expression studies using common bean. RESULTS: To evaluate the utility of the soybean GeneChip for transcript profiling of common bean, we hybridized cRNAs purified from nodule, leaf, and root of common bean and soybean in triplicate to the soybean GeneChip. Initial data analysis showed a decreased sensitivity and accuracy of measuring differential gene expression in common bean cross-species hybridization (CSH) GeneChip data compared to that of soybean. We employed a method that masked putative probes targeting inter-species variable (ISV) regions between common bean and soybean. A masking signal intensity threshold was selected that optimized both sensitivity and accuracy of measuring differential gene expression. After masking for ISV regions, the number of differentially-expressed genes identified in common bean was increased by 2.8-fold reflecting increased sensitivity. Quantitative RT-PCR (qRT-PCR) analysis of 20 randomly selected genes and purine-ureide pathway genes demonstrated an increased accuracy of measuring differential gene expression after masking for ISV regions. We also evaluated masked probe frequency per probe set to gain insight into the sequence divergence pattern between common bean and soybean. The sequence divergence pattern analysis suggested that the genes for basic cellular functions and metabolism were highly conserved between soybean and common bean. Additionally, our results show that some classes of genes, particularly those associated with environmental adaptation, are highly divergent. CONCLUSIONS: The soybean GeneChip is a suitable cross-species platform for transcript profiling in common bean when used in combination with the masking protocol described. In addition to transcript profiling, CSH of the GeneChip in combination with masking probes in the ISV regions can be used for comparative ecological and/or evolutionary genomics studies.

Knockdown of CELL DIVISION CYCLE16 reveals an inverse relationship between lateral root and nodule numbers and a link to auxin in Medicago truncatula.

The postembryonic development of lateral roots and nodules is a highly regulated process. Recent studies suggest the existence of cross talk and interdependency in the growth of these two organs. Although plant hormones, including auxin and cytokinin, appear to be key players in coordinating this cross talk, very few genes that cross-regulate root and nodule development have been uncovered so far. This study reports that a homolog of CELL DIVISION CYCLE16 (CDC16), a core component of the Anaphase Promoting Complex, is one of the key mediators in controlling the overall number of lateral roots and nodules. A partial suppression of this gene in Medicago truncatula leads to a decrease in number of lateral roots and a 4-fold increase in number of nodules. The roots showing lowered expression of MtCDC16 also show reduced sensitivity to phytohormone auxin, thus providing a potential function of CDC16 in auxin signaling.

Molecular analysis of SCARECROW genes expressed in white lupin cluster roots.

The Scarecrow (SCR) transcription factor plays a crucial role in root cell radial patterning and is required for maintenance of the quiescent centre and differentiation of the endodermis. In response to phosphorus (P) deficiency, white lupin (Lupinus albus L.) root surface area increases some 50-fold to 70-fold due to the development of cluster (proteoid) roots. Previously it was reported that SCR-like expressed sequence tags (ESTs) were expressed during early cluster root development. Here the cloning of two white lupin SCR genes, LaSCR1 and LaSCR2, is reported. The predicted amino acid sequences of both LaSCR gene products are highly similar to AtSCR and contain C-terminal conserved GRAS family domains. LaSCR1 and LaSCR2 transcript accumulation localized to the endodermis of both normal and cluster roots as shown by in situ hybridization and gene promoter::reporter staining. Transcript analysis as evaluated by quantitative real-time-PCR (qRT-PCR) and RNA gel hybridization indicated that the two LaSCR genes are expressed predominantly in roots. Expression of LaSCR genes was not directly responsive to the P status of the plant but was a function of cluster root development. Suppression of LaSCR1 in transformed roots of lupin and Medicago via RNAi (RNA interference) delivered through Agrobacterium rhizogenes resulted in decreased root numbers, reflecting the potential role of LaSCR1 in maintaining root growth in these species. The results suggest that the functional orthologues of AtSCR have been characterized.

Functional analysis and development of a CRISPR/Cas9 allelic series for a CPR5 ortholog necessary for proper growth of soybean trichomes.

Developments in genomic and genome editing technologies have facilitated the mapping, cloning, and validation of genetic variants underlying trait variation. This study combined bulked-segregant analysis, array comparative genomic hybridization, and CRISPR/Cas9 methodologies to identify a CPR5 ortholog essential for proper trichome growth in soybean (Glycine max). A fast neutron mutant line exhibited short trichomes with smaller trichome nuclei compared to its parent line. A fast neutron-induced deletion was identified within an interval on chromosome 6 that co-segregated with the trichome phenotype. The deletion encompassed six gene models including an ortholog of Arabidopsis thaliana CPR5. CRISPR/Cas9 was used to mutate the CPR5 ortholog, resulting in five plants harboring a total of four different putative knockout alleles and two in-frame alleles. Phenotypic analysis of the mutants validated the candidate gene, and included intermediate phenotypes that co-segregated with the in-frame alleles. These findings demonstrate that the CPR5 ortholog is essential for proper growth and development of soybean trichomes, similar to observations in A. thaliana. Furthermore, this work demonstrates the value of using CRISPR/Cas9 to generate an allelic series and intermediate phenotypes for functional analysis of candidate genes and/or the development of novel traits.

Molecular cloning, characterization and regulation of two different NADH-glutamate synthase cDNAs in bean nodules.

NADH-dependent glutamate synthase (NADH-GOGAT) is a key enzyme in primary ammonia assimilation in Phaseolus vulgaris nodules. Two different types of cDNA clones of PvNADH-GOGAT were isolated from the nodule cDNA libraries. The full-length cDNA clones of PvNADH-GOGAT-I (7.4 kb) and PvNADH-GOGAT-II (7.0 kb), which displayed an 83% homology between them, were isolated using cDNA library screening, 'cDNA library walking' and RT-PCR amplification. Southern analysis employing specific 5' cDNA probes derived from PvNADH-GOGAT-I and PvNADH-GOGAT-II indicated the existence of a single copy of each gene in the bean genome. Both these proteins contain approximately 100 amino acid sequences theoretically addressing each isoenzyme to different subcellular compartments. RT-PCR analysis indicated that PvNADH-GOGAT-II expression is higher than PvNADH-GOGAT-I during nodule development. Expression analysis by RT-PCR also revealed that both of these genes are differentially regulated by sucrose. On the other hand, the expression of PvNADH-GOGAT-I, but not PvNADH-GOGAT-II, was inhibited with nitrogen compounds. In situ hybridization and promoter expression analyses demonstrated that the NADH-GOGAT-I and -II genes are differentially expressed in bean root and nodule tissues. In silico analyses of the NADH-GOGAT promoters revealed the presence of potential cis elements in them that could mediate differential tissue-specific, and sugar and amino acid responsive expression of these genes.

Localization of superoxide dismutases and hydrogen peroxide in legume root nodules.

Superoxide dismutases (SODs) catalyze the dismutation of superoxide radicals to O2 and H2O2 and thus represent a primary line of antioxidant defense in all aerobic organisms. H2O2 is a signal molecule involved in the plant's response to pathogen attack and other stress conditions as well as in nodulation. In this work, we have tested the hypothesis that SODs are a source of H2O2 in indeterminate alfalfa (Medicago sativa) and pea (Pisum sativum) nodules. The transcripts and proteins of the major SODs of nodules were localized by in situ RNA hybridization and immunogold electron microscopy, respectively, whereas H2O2 was localized cytochemically by electron microscopy of cerium-perfused nodule tissue. The transcript and protein of cytosolic CuZnSOD are most abundant in the meristem (I) and invasion (II) zones, interzone II-III, and distal part of the N2-fixing zone (III), and those of MnSOD in zone III, especially in the infected cells. At the subcellular level, CuZnSOD was found in the infection threads, cytosol adjacent to cell walls, and apoplast, whereas MnSOD was in the bacteroids, bacteria within infection threads, and mitochondria. The distinct expression pattern of CuZnSOD and MnSOD suggests specific roles of the enzymes in nodules. Large amounts of H2O2 were found at the same three nodule sites as CuZnSOD but not in association with MnSOD. This colocalization led us to postulate that cytosolic CuZnSOD is a source of H2O2 in nodules. Furthermore, the absence or large reduction of H2O2 in nodule tissue preincubated with enzyme inhibitors (cyanide, azide, diphenyleneiodonium, diethyldithiocarbamate) provides strong support to the hypothesis that at least some of the H2O2 originates by the sequential operation of an NADPH oxidase-like enzyme and CuZnSOD. Results also show that there is abundant H2O2 associated with degrading bacteroids in the senescent zone (IV), which reflects the oxidative stress ensued during nodule senescence.

Regulatory patterns of a large family of defensin-like genes expressed in nodules of Medicago truncatula.

Root nodules are the symbiotic organ of legumes that house nitrogen-fixing bacteria. Many genes are specifically induced in nodules during the interactions between the host plant and symbiotic rhizobia. Information regarding the regulation of expression for most of these genes is lacking. One of the largest gene families expressed in the nodules of the model legume Medicago truncatula is the nodule cysteine-rich (NCR) group of defensin-like (DEFL) genes. We used a custom Affymetrix microarray to catalog the expression changes of 566 NCRs at different stages of nodule development. Additionally, bacterial mutants were used to understand the importance of the rhizobial partners in induction of NCRs. Expression of early NCRs was detected during the initial infection of rhizobia in nodules and expression continued as nodules became mature. Late NCRs were induced concomitantly with bacteroid development in the nodules. The induction of early and late NCRs was correlated with the number and morphology of rhizobia in the nodule. Conserved 41 to 50 bp motifs identified in the upstream 1,000 bp promoter regions of NCRs were required for promoter activity. These cis-element motifs were found to be unique to the NCR family among all annotated genes in the M. truncatula genome, although they contain sub-regions with clear similarity to known regulatory motifs involved in nodule-specific expression and temporal gene regulation.

A re-sequencing based assessment of genomic heterogeneity and fast neutron-induced deletions in a common bean cultivar.

A small fast neutron (FN) mutant population has been established from Phaseolus vulgaris cv. Red Hawk. We leveraged the available P. vulgaris genome sequence and high throughput next generation DNA sequencing to examine the genomic structure of five P. vulgaris cv. Red Hawk FN mutants with striking visual phenotypes. Analysis of these genomes identified three classes of structural variation (SV); between cultivar variation, natural variation within the FN mutant population, and FN induced mutagenesis. Our analyses focused on the latter two classes. We identified 23 large deletions (>40 bp) common to multiple individuals, illustrating residual heterogeneity and regions of SV within the common bean cv. Red Hawk. An additional 18 large deletions were identified in individual mutant plants. These deletions, ranging in size from 40 bp to 43,000 bp, are potentially the result of FN mutagenesis. Six of the 18 deletions lie near or within gene coding regions, identifying potential candidate genes causing the mutant phenotype.

A standardized method for analysis of Medicago truncatula phenotypic development.

Medicago truncatula has become a model system to study legume biology. It is imperative that detailed growth characteristics of the most commonly used cultivar, line A17 cv Jemalong, be documented. Such analysis creates a basis to analyze phenotypic alterations due to genetic lesions or environmental stress and is essential to characterize gene function and its relationship to morphological development. We have documented morphological development of M. truncatula to characterize its temporal developmental growth pattern; developed a numerical nomenclature coding system that identifies stages in morphological development; tested the coding system to identify phenotypic differences under phosphorus (P) and nitrogen (N) deprivation; and created visual models using the L-system formalism. The numerical nomenclature coding system, based on a series of defined growth units, represents incremental steps in morphological development. Included is a decimal component dividing growth units into nine substages. A measurement component helps distinguish alterations that may be missed by the coding system. Growth under N and P deprivation produced morphological alterations that were distinguishable using the coding system and its measurement component. N and P deprivation resulted in delayed leaf development and expansion, delayed axillary shoot emergence and elongation, decreased leaf and shoot size, and altered root growth. Timing and frequency of flower emergence in P-deprived plants was affected. This numerical coding system may be used as a standardized method to analyze phenotypic variation in M. truncatula due to nutrient stress, genetic lesions, or other factors and should allow valid growth comparisons across geographically distant laboratories.