Phosphorus (P) is an essential element for plant growth. Crop production of common bean (Phaseolus vulgaris), the most important legume for human consumption, is often limited by low P in the soil. Functional genomics were used to investigate global gene expression and metabolic responses of bean plants grown under P-deficient and P-sufficient conditions. P-deficient plants showed enhanced root to shoot ratio accompanied by reduced leaf area and net photosynthesis rates. Transcript profiling was performed through hybridization of nylon filter arrays spotted with cDNAs of 2,212 unigenes from a P deficiency root cDNA library. A total of 126 genes, representing different functional categories, showed significant differential expression in response to P: 62% of these were induced in P-deficient roots. A set of 372 bean transcription factor (TF) genes, coding for proteins with Inter-Pro domains characteristic or diagnostic for TF, were identified from The Institute of Genomic Research/Dana Farber Cancer Institute Common Bean Gene Index. Using real-time reverse transcription-polymerase chain reaction analysis, 17 TF genes were differentially expressed in P-deficient roots; four TF genes, including MYB TFs, were induced. Nonbiased metabolite profiling was used to assess the degree to which changes in gene expression in P-deficient roots affect overall metabolism. Stress-related metabolites such as polyols accumulated in P-deficient roots as well as sugars, which are known to be essential for P stress gene induction. Candidate genes have been identified that may contribute to root adaptation to P deficiency and be useful for improvement of common bean.

Cytosine DNA methylation is one avenue for passing information through cell divisions. Here, we present epigenomic analyses of soybean recombinant inbred lines (RILs) and their parents. Identification of differentially methylated regions (DMRs) revealed that DMRs mostly cosegregated with the genotype from which they were derived, but examples of the uncoupling of genotype and epigenotype were identified. Linkage mapping of methylation states assessed from whole-genome bisulfite sequencing of 83 RILs uncovered widespread evidence for local methylQTL. This epigenomics approach provides a comprehensive study of the patterns and heritability of methylation variants in a complex genetic population over multiple generations, paving the way for understanding how methylation variants contribute to phenotypic variation.

• MicroRNAs (miRNAs) play a pivotal role in post‐transcriptional regulation of gene expression in plants. Information on miRNAs in legumes is as yet scarce. This work investigates miRNAs in an agronomically important legume, common bean (Phaseolus vulgaris). • A hybridization approach employing miRNA macroarrays – printed with oligonucleotides complementary to 68 known miRNAs – was used to detect miRNAs in the leaves, roots and nodules of control and nutrient‐stressed (phosphorus, nitrogen, or iron deficiency; acidic pH; and manganese toxicity) common bean plants. • Thirty‐three miRNAs were expressed in control plants and another five were only expressed under stress conditions. The miRNA expression ratios (stress:control) were evaluated using principal component and hierarchical cluster analyses. A group of miRNAs responded to nearly all stresses in the three organs analyzed. Other miRNAs showed organ‐specific responses. Most of the nodule‐responsive miRNAs showed up‐regulation. miRNA blot expression analysis confirmed the macroarray results. Novel miRNA target genes were proposed for common bean and the expression of selected targets was evaluated by quantitative reverse transcriptase–polymerase chain reaction. • In addition to the detection of previously reported stress‐responsive miRNAs, we discovered novel common bean stress‐responsive miRNAs, for manganese toxicity. Our data provide a foundation for evaluating the individual roles of miRNAs in common bean.

Phosphorus (P), an essential element for plants, is one of the most limiting nutrients for plant growth. A few transcription factor (TF) genes involved in P‐starvation signalling have been characterized for Arabidopsis thaliana and rice. Crop production of common bean (Phaseolus vulgaris L.), the most important legume for human consumption, is often limited by low P in the soil. Despite its agronomic importance, nothing is known about transcriptional regulation in P‐deficient bean plants. We functionally characterized the P‐deficiency‐induced MYB TF TC3604 (Dana Farber Cancer Institute, Common Bean Gene Index v.2.0), ortholog to AtPHR1 (PvPHR1). For its study, we applied RNAi technology in bean composite plants. PvPHR1 is a positive regulator of genes implicated in P transport, remobilization and homeostasis. Although there are no reports on the regulatory roles of microRNAs (miRNA) in bean, we demonstrated that PvmiR399 is an essential component of the PvPHR1 signalling pathway. The analysis of DICER‐like1 (PvDCL1) silenced bean composite plants suppressed for accumulation of PvmiR399 and other miRNAs suggested that miR399 is a negative regulator of the ubiquitin E2 conjugase: PvPHO2 expression. Our results set the basis for understanding the signalling for P‐starvation responses in common bean and may contribute to crop improvement.

Phosphorus (P) deficiency is widespread in regions where the common bean (Phaseolus vulgaris), the most important legume for human consumption, is produced, and it is perhaps the factor that most limits nitrogen fixation. Global gene expression and metabolome approaches were used to investigate the responses of nodules from common bean plants inoculated with Rhizobium tropici CIAT899 grown under P-deficient and P-sufficient conditions. P-deficient inoculated plants showed drastic reduction in nodulation and nitrogenase activity as determined by acetylene reduction assay. Nodule transcript profiling was performed through hybridization of nylon filter arrays spotted with cDNAs, approximately 4,000 unigene set, from the nodule and P-deficient root library. A total of 459 genes, representing different biological processes according to updated annotation using the UniProt Knowledgebase database, showed significant differential expression in response to P: 59% of these were induced in P-deficient nodules. The expression platform for transcription factor genes based in quantitative reverse transcriptase-polymerase chain reaction revealed that 37 transcription factor genes were differentially expressed in P-deficient nodules and only one gene was repressed. Data from nontargeted metabolic profiles indicated that amino acids and other nitrogen metabolites were decreased, while organic and polyhydroxy acids were accumulated, in P-deficient nodules. Bioinformatics analyses using MapMan and PathExpress software tools, customized to common bean, were utilized for the analysis of global changes in gene expression that affected overall metabolism. Glycolysis and glycerolipid metabolism, and starch and Suc metabolism, were identified among the pathways significantly induced or repressed in P-deficient nodules, respectively.

Small RNAs are ubiquitous, versatile repressors and include (1) microRNAs (miRNAs), processed from mRNA forming stem-loops; and (2) small interfering RNAs (siRNAs), the latter derived in plants by a process typically requiring an RNA-dependent RNA polymerase. We constructed and analyzed an expression atlas of soybean (Glycine max) small RNAs, identifying over 500 loci generating 21-nucleotide phased siRNAs (phasiRNAs; from PHAS loci), of which 483 overlapped annotated protein-coding genes. Via the integration of miRNAs with parallel analysis of RNA end (PARE) data, 20 miRNA triggers of 127 PHAS loci were detected. The primary class of PHAS loci (208 or 41% of the total) corresponded to NB-LRR genes; some of these small RNAs preferentially accumulate in nodules. Among the PHAS loci, novel representatives of TAS3 and noncanonical phasing patterns were also observed. A noncoding PHAS locus, triggered by miR4392, accumulated preferentially in anthers; the phasiRNAs are predicted to target transposable elements, with their peak abundance during soybean reproductive development. Thus, phasiRNAs show tremendous diversity in dicots. We identified novel miRNAs and assessed the veracity of soybean miRNAs registered in miRBase, substantially improving the soybean miRNA annotation, facilitating an improvement of miRBase annotations and identifying at high stringency novel miRNAs and their targets.