Expanding the eukaryotic genetic code with a biosynthesized 21st amino acid
Explore the source record for details and available documents.
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
Anaerobic fungi are emerging biotechnology platforms with genomes rich in biosynthetic potential. Yet, the heterologous expression of their biosynthetic pathways has had limited success in model hosts like E. coli. We find one reason for this is that the genome composition of anaerobic fungi like P. indianae are extremely AT-biased with a particular preference for rare and semi-rare AT-rich tRNAs in E coli, which are not explicitly predicted by standard codon adaptation indices (CAI). Native P. indianae genes with these extreme biases create drastic growth defects in E. coli (up to 69% reduction in growth), which is not seen in genes from other organisms with similar CAIs. However, codon optimization rescues growth, allowing for gene evaluation. In this manner, we demonstrate that anaerobic fungal homologs such as PI.atoB are more active than S. cerevisiae homologs in a hybrid pathway, increasing the production of mevalonate up to 2.5 g/L (more than two-fold) and reducing waste carbon to acetate by ~90% under the conditions tested. This work demonstrates the bioproduction potential of anaerobic fungal enzyme homologs and how the analysis of codon utilization enables the study of otherwise difficult to express genes that have applications in biocatalysis and natural product discovery.
Understanding the interactions between plants and microorganisms can inform microbiome management to enhance crop productivity and resilience to stress. Here, we apply a genome-centric approach to identify ecologically important leaf microbiome members on replicated plots of field-grown switchgrass and miscanthus, and to quantify their activities over two growing seasons for switchgrass. We use metagenome and metatranscriptome sequencing and curate 40 medium- and high-quality metagenome-assembled-genomes (MAGs). We find that classes represented by these MAGs (Actinomycetia, Alpha- and Gamma- Proteobacteria, and Bacteroidota) are active in the late season, and upregulate transcripts for short-chain dehydrogenase, molybdopterin oxidoreductase, and polyketide cyclase. Stress-associated pathways are expressed for most MAGs, suggesting engagement with the host environment. We also detect seasonally activated biosynthetic pathways for terpenes and various non-ribosomal peptide pathways that are poorly annotated. Our findings support that leaf-associated bacterial populations are seasonally dynamic and responsive to host cues.
The evolution of pathway enzymes enhances the biosynthesis of high-value chemicals, crucial for pharmaceutical, and agrochemical applications. However, unpredictable evolutionary landscapes of pathway genes often hinder successful evolution. Here, the presence of complex epistasis is identifued within the representative naringenin biosynthetic pathway enzymes, hampering straightforward directed evolution. Subsequently, a biofoundry-assisted strategy is developed for pathway bottlenecking and debottlenecking, enabling the parallel evolution of all pathway enzymes along a predictable evolutionary trajectory in six weeks. This study then utilizes a machine learning model, ProEnsemble, to further balance the pathway by optimizing the transcription of individual genes. The broad applicability of this strategy is demonstrated by constructing an Escherichia coli chassis with evolved and balanced pathway genes, resulting in 3.65 g L -1 naringenin. The optimized naringenin chassis also demonstrates enhanced production of other flavonoids. This approach can be readily adapted for any given number of enzymes in the specific metabolic pathway, paving the way for automated chassis construction in contemporary biofoundries.
A series of porphyrins and catalysts has been prepared as a model for the origin of photosynthesis on the primordial earth. These compounds have been used to test the hypotheses that (1) the biosynthetic pathway to chlorophyll recapitulates the evolutionary history of photosynthesis, and (2) the proto-photosythetic function of biogenetic porphyrins (biosynthetic chlorophyll precursors) was the oxidation of organic molecules by photoexcited porphyrins with the attendant emission of molecular hydrogen. This paper describes experiments in which photoexcited biogenetic porphyrins oxidize ethylenediamine tetraacetic acid (EDTA). The concomitant reduction of protons to hydrogen gas occurs in the presence of a colloidal platinum catalyst. The addition of methyl viologen, a one-electron shuttle, increases the amount of molecular hydrogen generated during long irradiations and the quantum yield of hydrogen production. When the porphyrin and catalyst are held in association by molecular complexes, the increased efficiency of electron transfer produces higher yields of hydrogen gas.
The design and optimization of biosynthetic pathways for industrially relevant, non-model organisms is challenging due to transformation idiosyncrasies, reduced numbers of validated genetic parts and a lack of high-throughput workflows. Here we describe a platform for in vitro prototyping and rapid optimization of biosynthetic enzymes (iPROBE) to accelerate this process. In iPROBE, cell lysates are enriched with biosynthetic enzymes by cell-free protein synthesis and then metabolic pathways are assembled in a mix-and-match fashion to assess pathway performance. We demonstrate iPROBE by screening 54 different cell-free pathways for 3-hydroxybutyrate production and optimizing a six-step butanol pathway across 205 permutations using data-driven design. Observing a strong correlation (r = 0.79) between cell-free and cellular performance, we then scaled up our highest-performing pathway, which improved in vivo 3-HB production in Clostridium by 20-fold to 14.63 ± 0.48 g l -1 . Finally, we expect iPROBE to accelerate design–build–test cycles for industrial biotechnology.
Abstract Background The oleaginous yeast Rhodotorula toruloides is a promising chassis organism for the biomanufacturing of value-added bioproducts. It can accumulate lipids at a high fraction of biomass. However, metabolic engineering efforts in this organism have progressed at a slower pace than those in more extensively studied yeasts. Few studies have investigated the lipid accumulation phenotype exhibited by R. toruloides under nitrogen limitation conditions. Consequently, there have been only a few studies exploiting the lipid metabolism for higher product titers. Results We performed a multi-omic investigation of the lipid accumulation phenotype under nitrogen limitation. Specifically, we performed comparative transcriptomic and lipidomic analysis of the oleaginous yeast under nitrogen-sufficient and nitrogen deficient conditions. Clustering analysis of transcriptomic data was used to identify the growth phase where nitrogen-deficient cultures diverged from the baseline conditions. Independently, lipidomic data was used to identify that lipid fractions shifted from mostly phospholipids to mostly storage lipids under the nitrogen-deficient phenotype. Through an integrative lens of transcriptomic and lipidomic analysis, we discovered that R. toruloides undergoes lipid remodeling during nitrogen limitation, wherein the pool of phospholipids gets remodeled to mostly storage lipids. We identify specific mRNAs and pathways that are strongly correlated with an increase in lipid levels, thus identifying putative targets for engineering greater lipid accumulation in R. toruloides . One surprising pathway identified was related to inositol phosphate metabolism, suggesting further inquiry into its role in lipid accumulation. Conclusions Integrative analysis identified the specific biosynthetic pathways that are differentially regulated during lipid remodeling. This insight into the mechanisms of lipid accumulation can lead to the success of future metabolic engineering strategies for overproduction of oleochemicals.
Abstract Modification of lignin in feedstocks via genetic engineering aims to reduce biomass recalcitrance to facilitate efficient conversion processes. These improvements can be achieved by expressing exogenous enzymes that interfere with native biosynthetic pathways responsible for the production of the lignin precursors. In planta expression of a bacterial 3-dehydroshikimate dehydratase in poplar trees reduced lignin content and altered the monomer composition, which enabled higher yields of sugars after cell wall polysaccharide hydrolysis. Understanding how plants respond to such genetic modifications at the transcriptional and metabolic levels is needed to facilitate further improvement and field deployment. In this work, we acquired fundamental knowledge on lignin-modified poplar expressing 3-dehydroshikimate dehydratase using RNA-seq and metabolomics. The data clearly demonstrate that changes in gene expression and metabolite abundance can occur in a strict spatiotemporal fashion, revealing tissue-specific responses in the xylem, phloem, or periderm. In the poplar line that exhibited the strongest reduction in lignin, we found that 3% of the transcripts had altered expression levels and ~19% of the detected metabolites had differential abundance in the xylem from older stems. The changes affected predominantly the shikimate and phenylpropanoid pathways as well as secondary cell wall metabolism, and resulted in significant accumulation of hydroxybenzoates derived from protocatechuate and salicylate.
Coagulase-positive staphylococci, which frequently colonize the mucosal surfaces of animals, also cause a spectrum of opportunistic infections including skin and soft tissue infections, urinary tract infections, pneumonia, and bacteremia. However, recent advances in bacterial identification have revealed that these common veterinary pathogens are in fact zoonoses that cause serious infections in human patients. The global spread of multidrug-resistant zoonotic staphylococci, in particular the emergence of methicillin-resistant organisms, is now a serious threat to both animal and human welfare. Accordingly, new therapeutic targets that can be exploited to combat staphylococcal infections are urgently needed. Enzymes of the methylerythritol phosphate pathway (MEP) of isoprenoid biosynthesis represent potential targets for treating zoonotic staphylococci. Here we demonstrate that fosmidomycin (FSM) inhibits the first step of the isoprenoid biosynthetic pathway catalyzed by deoxyxylulose phosphate reductoisomerase (DXR) in staphylococci. In addition, we have both enzymatically and structurally determined the mechanism by which FSM elicits its effect. Using a forward genetic screen, the glycerol-3-phosphate transporter GlpT that facilitates FSM uptake was identified in two zoonotic staphylococci, Staphylococcus schleiferi and Staphylococcus pseudintermedius . A series of lipophilic ester prodrugs (termed MEPicides) structurally related to FSM were synthesized, and data indicate that the presence of the prodrug moiety not only substantially increased potency of the inhibitors against staphylococci but also bypassed the need for GlpT-mediated cellular transport. Collectively, our data indicate that the prodrug MEPicides selectively and robustly inhibit DXR in zoonotic staphylococci, and further, that DXR represents a promising, druggable target for future development.
The raw data (metagenomes and metatranscriptomes) for this study are available in the Joint Genomes Institute Genome Portal (https://genome.jgi.doe.gov/portal/ Project ID 503249) with projects designated by year and product type. The MAG genomes analyzed in this paper are available on NCBI, as bioproject PRJNA800073. Plants and microorganisms form beneficial associations. Understanding plant-microbe interactions will inform microbiome management to enhance crop productivity and resilience to stress. Here, we apply a genome-centric approach to identify ecologically important leaf microbiome members on field-grown switchgrass and miscanthus and to quantify their activities for switchgrass over two growing seasons. We integrate metagenome and metatranscriptome sequencing from 192 leaf samples collected over representative time points in crop phenology. We curated 40 medium- and high-quality metagenome-assembled-genomes (MAGs) and focused analysis on seasonal transcript recruitment to them. Classes represented by these focal MAGs (Actinomycetia, Alpha- and Gamma- Proteobacteria, and Bacteroidota) were active and had increases in transcripts for short-chain dehydrogenase, molybdopterin oxidoreductase, and polyketide cyclase in the late season. The majority of MAGs had activated stress-associated pathways, including trehalose metabolism, indole acetic acid degradation, betaine biosynthesis, and reactive oxygen species degradation, suggesting direct engagement with the host environment. We also detected seasonally activated biosynthetic pathways for terpenes (carotenoids and isoprenoids) and for various non-ribosomal peptide pathways that were poorly annotated. Overall, this study overcame laboratory and bioinformatic challenges associated with field-based leaf metatranscriptome analysis to inform both general and likely specialized activities of these phyllosphere populations. These activities collectively support that leaf-associated bacterial populations are seasonally dynamic, responsive to host cues, and interactively engage in feedback with the plant. This analysis represented quality filtering of metagenomes and metatranscriptomes (data-preparation folder), metagenome assemblies (metagenome-assembly folder) and metagenome-assembled genome binning, curation, refinement, annotation (mag-evaluation folder). Abundances of sequencing libraries were calculated based on reads mapped (mapping folder). Additionallly, analysis of our annotated results are also included (analysis folder).
Abstract Plants make complex and potent therapeutic molecules 1,2 , but sourcing these molecules from natural producers or through chemical synthesis is difficult, which limits their use in the clinic. A prominent example is the anti-cancer therapeutic paclitaxel (sold under the brand name Taxol), which is derived from yew trees (Taxusspecies) 3 . Identifying the full paclitaxel biosynthetic pathway would enable heterologous production of the drug, but this has yet to be achieved despite half a century of research 4 . WithinTaxus’ large, enzyme-rich genome 5 , we suspected that the paclitaxel pathway would be difficult to resolve using conventional RNA-sequencing and co-expression analyses. Here, to improve the resolution of transcriptional analysis for pathway identification, we developed a strategy we term multiplexed perturbation × single nuclei (mpXsn) to transcriptionally profile cell states spanning tissues, cell types, developmental stages and elicitation conditions. Our data show that paclitaxel biosynthetic genes segregate into distinct expression modules that suggest consecutive subpathways. These modules resolved seven new genes, allowing a de novo 17-gene biosynthesis and isolation of baccatin III, the industrial precursor to Taxol, inNicotiana benthamianaleaves, at levels comparable with the natural abundance inTaxusneedles. Notably, we found that a nuclear transport factor 2 (NTF2)-like protein, FoTO1, is crucial for promoting the formation of the desired product during the first oxidation, resolving a long-standing bottleneck in paclitaxel pathway reconstitution. Together with a new β-phenylalanine-CoA ligase, the eight genes discovered here enable the de novo biosynthesis of 3’-N-debenzoyl-2’-deoxypaclitaxel. More broadly, we establish a generalizable approach to efficiently scale the power of co-expression analysis to match the complexity of large, uncharacterized genomes, facilitating the discovery of high-value gene sets.
Abstract Streptomycetes are highly metabolically gifted bacteria with the abilities to produce bioproducts that have profound economic and societal importance. These bioproducts are produced by metabolic pathways including those for the biosynthesis of secondary metabolites and catabolism of plant biomass constituents. Advancements in genome sequencing technologies have revealed a wealth of untapped metabolic potential from Streptomyces genomes. Here, we report the largest Streptomyces pangenome generated by using 205 complete genomes. Metabolic potentials of the pangenome and individual genomes were analyzed, revealing degrees of conservation of individual metabolic pathways and strains potentially suitable for metabolic engineering. Of them, Streptomyces bingchenggensis was identified as a potent degrader of plant biomass. Polyketide, non-ribosomal peptide, and gamma-butyrolactone biosynthetic enzymes are primarily strain specific while ectoine and some terpene biosynthetic pathways are highly conserved. A large number of transcription factors associated with secondary metabolism are strain-specific while those controlling basic biological processes are highly conserved. Although the majority of genes involved in morphological development are highly conserved, there are strain-specific varieties which may contribute to fine tuning the timing of cellular differentiation. Overall, these results provide insights into the metabolic potential, regulation and physiology of streptomycetes, which will facilitate further exploitation of these important bacteria.
Nitrogenase is an active target of heterologous expression because of its importance for areas related to agronomy, energy, and environment. One major hurdle for expressing an active Mo-nitrogenase in Escherichia coli is to generate the complex metalloclusters (P- and M-clusters) within this enzyme, which involves some highly unique bioinorganic chemistry/metalloenzyme biochemistry that is not generally dealt with in the heterologous expression of proteins via synthetic biology; in particular, the heterologous synthesis of the homometallic P-cluster ([Fe 8 S 7 ]) and M-cluster core (or L-cluster; [Fe 8 S 9 C]) on their respective protein scaffolds, which represents two crucial checkpoints along the biosynthetic pathway of a complete nitrogenase, has yet to be demonstrated by biochemical and spectroscopic analyses of purified metalloproteins. Here, we report the heterologous formation of a P-cluster-containing NifDK protein upon coexpression of Azotobacter vinelandii nifD, nifK, nifH, nifM, and nifZ genes, and that of an L-cluster-containing NifB protein upon coexpression of Methanosarcina acetivorans nifB, nifS, and nifU genes alongside the A. vinelandii fdxN gene, in E. coli. Our metal content, activity, EPR, and XAS/EXAFS data provide conclusive evidence for the successful synthesis of P- and L-clusters in a nondiazotrophic host, thereby highlighting the effectiveness of our metallocentric, divide-and-conquer approach that individually tackles the key events of nitrogenase biosynthesis prior to piecing them together into a complete pathway for the heterologous expression of nitrogenase. As such, this work paves the way for the transgenic expression of an active nitrogenase while providing an effective tool for further tackling the biosynthetic mechanism of this important metalloenzyme.
The siderophore rhizoferrin (N 1 ,N 4 -dicitrylputrescine) is produced in fungi and bacteria to scavenge iron. Putrescine-producing bacterium Ralstonia pickettii synthesizes rhizoferrin and encodes a single nonribosomal peptide synthetase-independent siderophore (NIS) synthetase. From biosynthetic logic, we hypothesized that this single enzyme is sufficient for rhizoferrin biosynthesis. We confirmed this by expression of R. pickettii NIS synthetase in Escherichia coli, resulting in rhizoferrin production. This was further confirmed in vitro using the recombinant NIS synthetase, synthesizing rhizoferrin from putrescine and citrate. Heterologous expression of homologous lbtA from Legionella pneumophila, required for rhizoferrin biosynthesis in that species, produced siderophore activity in E. coli. Rhizoferrin is also synthesized by Francisella tularensis and Francisella novicida, but unlike R. pickettii or L. pneumophila, Francisella species lack putrescine biosynthetic pathways because of genomic decay. Francisella encodes a NIS synthetase FslA/FigA and an ornithine decarboxylase homolog FslC/FigC, required for rhizoferrin biosynthesis. Ornithine decarboxylase produces putrescine from ornithine, but we show here in vitro that FigA synthesizes N-citrylornithine, and FigC is an N-citrylornithine decarboxylase that together synthesize rhizoferrin without using putrescine. We co-expressed F. novicida figA and figC in E. coli and produced rhizoferrin. A 2.1 Å X-ray crystal structure of the FigC N-citrylornithine decarboxylase reveals how the larger substrate is accommodated and how active site residues have changed to recognize N-citrylornithine. FigC belongs to a new subfamily of alanine racemase-fold PLP-dependent decarboxylases that are not involved in polyamine biosynthesis. These data reveal a natural product biosynthetic workaround that evolved to bypass a missing precursor and re-establish it in the final structure.
Abstract Despite the widespread emergence of multidrug‐resistant nosocomial Gram‐negative bacterial infections and the major public health threat it brings, no new class of antibiotics for Gram‐negative pathogens has been approved over the past five decades. Therefore, there is an urgent medical need for developing effective novel antibiotics against multidrug‐resistant Gram‐negative pathogens by targeting previously unexploited pathways in these bacteria. To fulfill this crucial need, we have been investigating a series of sulfonyl piperazine compounds targeting LpxH, a dimanganese‐containing UDP‐2,3‐diacylglucosamine hydrolase in the lipid A biosynthetic pathway, as novel antibiotics against clinically important Gram‐negative pathogens. Inspired by a detailed structural analysis of our previous LpxH inhibitors in complex with K. pneumoniae LpxH ( Kp LpxH), here we report the development and structural validation of the first‐in‐class sulfonyl piperazine LpxH inhibitors, JH‐LPH‐45 ( 8 ) and JH‐LPH‐50 ( 13 ), that achieve chelation of the active site dimanganese cluster of Kp LpxH. The chelation of the dimanganese cluster significantly improves the potency of JH‐LPH‐45 ( 8 ) and JH‐LPH‐50 ( 13 ). We expect that further optimization of these proof‐of‐concept dimanganese‐chelating LpxH inhibitors will ultimately lead to the development of more potent LpxH inhibitors for targeting multidrug‐resistant Gram‐negative pathogens.
The Krebs cycle (citric acid or tricarboxylic acid cycle), the final common pathway in aerobic metabolism for the oxidation of carbohydrates, fatty acids and amino acids, is known to be irreversible. It liberates CO2 and generates NADH whose aerobic oxidation yields ATP but it does not operate in reverse as a biosynthetic pathway for CO2 assimilation. In 1966, our laboratory described a cyclic pathway for CO2 assimilation (Evans, Buchanan and Arnon 1966) that was unusual in two respects: (i) it provided the first instance of an obligate photoautotroph that assimilated CO2 by a pathway different from Calvin's reductive pentose phosphate cycle (Calvin 1962) and (ii) in its overall effect the new cycle was a reversal of the Krebs cycle. Named the 'reductive carboxylic acid cycle' (sometimes also called the reductive tricarboxylic acid cycle) the new cycle appeared to be the sole CO2 assimilation pathway in Chlorobium thiosulfatophilum (Evans et al. 1966) (now known as Chlorobium limicola forma thiosulfatophilum). Chlorobium is a photosynthetic green sulfur bacterium that grows anaerobically in an inorganic medium with sulfide and thiosulfate as electron donors and CO2 as an obligatory carbon source. In the ensuing years, the new cycle was viewed with skepticism. Not only was it in conflict with the prevailing doctrine that the 'one important property ... shared by all (our emphasis) autotrophic species is the assimilation of CO2 via the Calvin cycle' (McFadden 1973) but also some of its experimental underpinnings were challenged. It is only now that in the words of one of its early skeptics (Tabita 1988) 'a long and tortuous controversy' has ended with general acceptance of the reductive carboxylic acid cycle as a photosynthetic CO2 assimilation pathway distinct from the pentose cycle. (Henceforth, to minimize repetitiveness, the reductive pentose phosphate cycle will often be referred to as the pentose cycle and the reductive carboxylic acid cycle as the carboxylic acid cycle.) Aside from photosynthetic pathways which are the focus of this article, CO2 assimilation is also known to sustain autotrophic growth via the acetyl-CoA pathway (Wood et al. 1986). Our aim here is to discuss (i) the findings that led our group to the discovery of the reductive carboxylic acid cycle, (ii) the nature and resolution of the controversy that followed, and (iii) the possible evolutionary implications of the cycle as an ancient mechanism for photosynthetic CO2 assimilation that preceded the pentose cycle and served as a precursor of the Krebs cycle in aerobic metabolism.
Gmelina arborea Roxb is a fast-growing tree species of commercial importance for tropical countries due to multiple industrial uses of its wood. Wood is primarily composed of thick secondary cell walls of xylem cells which imparts the strength to the wood. Identification of the genes involved in the secondary cell wall biosynthesis as well as their cognate regulators is crucial to understand how the production of wood occurs and serves as a starting point for developing breeding strategies to produce varieties with improved wood quality, better paper pulping or new potential uses such as biofuel production. In order to gain knowledge on the molecular mechanisms and gene regulation related with wood development in white teak, a de novo sequencing and transcriptome assembly approach was used employing secondary cell wall synthesizing cells from young white teak trees. For generation of transcriptome, RNA-seq reads were assembled into 110,992 transcripts and 49,364 genes were functionally annotated using plant databases; 5071 GO terms and 25,460 SSR markers were identified within xylem transcripts and 10,256 unigenes were assigned to KEGG database in 130 pathways. Among transcription factor families, C2H2, C3H, bLHLH and MYB were the most represented in xylem. Differential gene expression analysis using leaves as a reference was carried out and a total of 20,954 differentially expressed genes were identified including monolignol biosynthetic pathway genes. The differential expression of selected genes (4CL, COMT, CCoAOMT, CCR and NST1) was validated using qPCR. We report the very first de novo transcriptome of xylem-related genes in this tropical timber species of commercial importance and constitutes a valuable extension of the publicly available transcriptomic resource aimed at fostering both basic and breeding studies.
Common bean (Phaseolus vulgaris L.) market classes have distinct seed coat colors, which are directly related to the diverse flavonoids found in the mature seed coat. To understand and elucidate the molecular mechanisms underlying the regulation of seed coat color, RNA-Seq data was collected from the black bean 5-593 and used for a differential gene expression and enrichment analysis from four different seed coat color development stages. 5-593 carries dominant alleles for 10 of the 11 major genes that control seed coat color and expression and has historically been used to develop introgression lines used for seed coat genetic analysis. Pairwise comparison among the four stages identified 6,294 differentially expressed genes (DEGs) varying from 508 to 5,780 DEGs depending on the compared stages. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed that phenylpropanoid biosynthesis, flavonoid biosynthesis, and plant hormone signal transduction comprised the principal pathways expressed during bean seed coat pigment development. Transcriptome analysis suggested that most structural genes for flavonoid biosynthesis and some potential regulatory genes were significantly differentially expressed. Further studies detected 29 DEGs as important candidate genes governing the key enzymatic flavonoid biosynthetic pathways for common bean seed coat color development. Additionally, four gene models, Pv5-593.02G016100, 593.02G078700, Pv5-593.02G090900, and Pv5-593.06G121300, encode MYB-like transcription factor family protein were identified as strong candidate regulatory genes in anthocyanin biosynthesis which could regulate the expression levels of some important structural genes in flavonoid biosynthesis pathway. These findings provide a framework to draw new insights into the molecular networks underlying common bean seed coat pigment development.