The Vigna unguiculata Gene Expression Atlas (VuGEA) from de novo assembly and quantification of RNA-seq data provides insights into seed maturation mechanisms

The pattern of genetic variability in a core collection of 2,021 cowpea accessions

Cowpea is a highly drought-adapted leguminous crop with great promise for improving agricultural sustainability and food security. Here, we report analyses derived from array-based genotyping of 2,021 accessions constituting a core subset of the world's largest cowpea collection, held at the International Institute of Tropical Agriculture (IITA) in Ibadan, Nigeria. We used this dataset to examine genetic variation and population structure in worldwide cowpea. We confirm that the primary pattern of population structure is two geographically defined subpopulations originating in West and East Africa, respectively, and that population structure is associated with shifts in phenotypic distribution. Furthermore, we establish the cowpea core collection as a resource for genome-wide association studies by mapping the genetic basis of several phenotypes, with a focus on seed coat pigmentation patterning and color. We anticipate that the genotyped IITA Cowpea Core Collection will serve as a powerful tool for mapping complex traits, facilitating the acceleration of breeding programs to enhance the resilience of this crop in the face of rapid global climate change.

The Content of Anthocyanins in Cowpea (Vigna unguiculata (L.) Walp.) Seeds and Contribution of the MYB Gene Cluster to Their Coloration Pattern

The intensively pigmented legumes belonging to Phaseolus and Vigna spp. are valued as an essential component of healthy nutrition due to their high content of flavonoids. In this context, we used the accessions of Vigna unguiculata with different colors of seed coats from the N.I. Vavilov All-Russian Institute of Plant Genetic Resources collection as the main object of this research. We applied confocal laser scanning microscopy, biochemical analysis, and wide in silico and molecular genetic analyses to study the main candidate genes for anthocyanin pigmentation within the MYB cluster on chromosome 5. We performed statistical data processing. The anthocyanin content ranged from 2.96 mg/100 g DW in reddish-brown-seeded cowpea accessions to 175.16 mg/100 g DW in black-seeded ones. Laser microscopy showed that the autofluorescence in cowpea seeds was mainly caused by phenolic compounds. The maximum fluorescence was observed in the seed coat, while its dark color, due to the highest level of red fluorescence, pointed to the presence of anthocyanins and anthocyanidins. Genes of the MYB cluster on chromosome 5 demonstrated a high homology and were segregated into a separate clade. However, amplification products were not obtained for all genes because of the truncation of some genes. Statistical analysis showed a clear correlation between the high content of anthocyanins in cowpea seeds and the presence of PCR products with primers Vigun05g0393-300-1.

Parallel selection of loss-of-function alleles of Pdh1 orthologous genes in warm-season legumes for pod indehiscence and plasticity is related to precipitation

Pod dehiscence facilitates seed dispersal in wild legumes but results in yield loss in cultivated legumes. The evolutionary genetics of the legume pod dehiscence trait remain largely elusive. We characterized the pod dehiscence of chromosome segment substitution lines of Glycine max crossed with Glycine soja and found that the gene underlying the predominant quantitative trait locus (QTL) of soybean pod-shattering trait was Pod dehiscence 1 (Pdh1). A few rare loss-of-function (LoF) Pdh1 alleles were identified in G. soja, while only an allele featuring a premature stop codon was selected for pod indehiscence in cultivated soybean and spread to low-precipitation regions with increased frequency. Moreover, correlated interactions among Pdh1's haplotype, gene expression, and environmental changes for the developmental plasticity of the pod dehiscence trait were revealed in G. max. We found that orthologous Pdh1 genes specifically originated in warm-season legumes and their LoF alleles were then parallel-selected during the domestication of legume crops. Our results provide insights into the convergent evolution of pod dehiscence in warm-season legumes, facilitate an understanding of the intricate interactions between genetic robustness and environmental adaptation for developmental plasticity, and guide the breeding of new legume varieties with pod indehiscence.

A comparative analysis of stably expressed genes across diverse angiosperms exposes flexibility in underlying promoter architecture

Promoters regulate both the amplitude and pattern of gene expression-key factors needed for optimization of many synthetic biology applications. Previous work in Arabidopsis found that promoters that contain a TATA-box element tend to be expressed only under specific conditions or in particular tissues, while promoters which lack any known promoter elements, thus designated as Coreless, tend to be expressed more ubiquitously. To test whether this trend represents a conserved promoter design rule, we identified stably expressed genes across multiple angiosperm species using publicly available RNA-seq data. Comparisons between core promoter architectures and gene expression stability revealed differences in core promoter usage in monocots and eudicots. Furthermore, when tracing the evolution of a given promoter across species, we found that core promoter type was not a strong predictor of expression stability. Our analysis suggests that core promoter types are correlative rather than causative in promoter expression patterns and highlights the challenges in finding or building constitutive promoters that will work across diverse plant species.

Legume-wide comparative analysis of pod shatter locus PDH1 reveals phaseoloid specificity, high cowpea expression, and stress responsive genomic context

Pod dehiscence is a major source of yield loss in legumes, which is exacerbated by aridity. Disruptive mutations in "Pod indehiscent 1" (PDH1), a pod sclerenchyma-specific lignin biosynthesis gene, has been linked to significant reductions in dehiscence in several legume species. We compared syntenic PDH1 regions across 12 legumes and two outgroups to uncover key historical evolutionary trends at this important locus. Our results clarified the extent to which PDH1 orthologs are present in legumes, showing the typical genomic context surrounding PDH1 has only arisen relatively recently in certain phaseoloid species (Vigna, Phaseolus, Glycine). The notable absence of PDH1 in Cajanus cajan may be a major contributor to its indehiscent phenotype compared with other phaseoloids. In addition, we identified a novel PDH1 ortholog in Vigna angularis and detected remarkable increases in PDH1 transcript abundance during Vigna unguiculata pod development. Investigation of the shared genomic context of PDH1 revealed it lies in a hotspot of transcription factors and signaling gene families that respond to abscisic acid and drought stress, which we hypothesize may be an additional factor influencing expression of PDH1 under specific environmental conditions. Our findings provide key insights into the evolutionary history of PDH1 and lay the foundation for optimizing the pod dehiscence role of PDH1 in major and understudied legume species.

Rhizosphere Priestia species altered cowpea root transcriptome and enhanced growth under drought and nutrient deficiency

Priestia species isolated from the cowpea rhizosphere altered the transcriptome of cowpea roots by colonization and enhanced nutrient uptake, antioxidant mechanisms, and photosynthesis, protecting cowpea from drought and nutrient deficiency. Cowpea is a significant grain legume crop primarily grown in sub-Saharan Africa, Asia, and South America. Drought and nutrient deficiency affect the growth and yield of cowpea. To address this challenge, we studied the phyto-beneficial effects of stress-tolerant rhizobacteria on the biomass yield of cowpea under water- and nutrient-deficit conditions. Among the bacteria isolated, two rhizobacillus genotypes, C8 (Priestia filamentosa; basonym: Bacillus filamentosus) and C29 (Priestia aryabhattai; basonym: Bacillus aryabhattai) were evaluated for the improvement of seed germination and growth of cowpea under stress. Our study revealed that C8 protected cowpea from stress by facilitating phosphorus and potassium uptake, protecting it from oxidative damage, reducing transpiration, and enhancing CO₂ assimilation. A 17% increase in root biomass upon C8 inoculation was concomitant with the induction of stress tolerance genes in cowpea roots predominantly involved in growth and metabolic processes, cell wall organization, ion homeostasis, and cellular responses to phosphate starvation. Our results indicate a metabolic alteration in cowpea root triggered by P. filamentosa, leading to efficient nutrient reallocation in the host plant. We propose inoculation with P. filamentosa as an effective strategy for improving the yield of cowpea in low-input agriculture, where chemical fertilization and irrigation are less accessible to resource-poor farmers.

A Genome-Wide Scan Divulges Key Loci Involved in Resistance to Aphids (Aphis craccivora) in Cowpea (Vigna unguiculata)

Cowpea aphids (Aphis craccivora Koch) double as a direct damaging pest and a virus vector to cowpea, threatening the economic yield of the crop. Given the multiple ecotypes, different alleles have been implicated in aphid resistance, necessitating the identification of key genes involved. The present study implemented a genome-wide scan using 365 cowpea mini-core accessions to decipher loci involved in resistance to aphid ecotype from Kano, Nigeria. Accessions were artificially inoculated with A. craccivora in insect-proof cages and damage severity assessed at 21 days after infestation. Significant phenotypic differences based on aphid damage severity were registered among the accessions. Skewed phenotypic distributions were depicted in the population, suggesting the involvement of major genes in the control of resistance. A genome-wide scan identified three major regions on chromosomes Vu10, Vu08 and Vu02, and two minor ones on chromosomes Vu01 and Vu06, that were significantly associated with aphid resistance. These regions harbored several genes, out of which, five viz Vigun01g233100.1, Vigun02g088900.1, Vigun06g224900.1, Vigun08g030200.1 and Vigun10g031100.1 were the most proximal to the peak single nucleotide polymorphisms (SNPs) positions. These genes are expressed under stress signaling, mechanical wounding and insect feeding. The uncovered loci contribute towards establishing a marker-assisted breeding platform and building durable resistance against aphids in cowpea.

An integrated transcriptome mapping the regulatory network of coding and long non-coding RNAs provides a genomics resource in chickpea

Large-scale transcriptome analysis can provide a systems-level understanding of biological processes. To accelerate functional genomic studies in chickpea, we perform a comprehensive transcriptome analysis to generate full-length transcriptome and expression atlas of protein-coding genes (PCGs) and long non-coding RNAs (lncRNAs) from 32 different tissues/organs via deep sequencing. The high-depth RNA-seq dataset reveal expression dynamics and tissue-specificity along with associated biological functions of PCGs and lncRNAs during development. The coexpression network analysis reveal modules associated with a particular tissue or a set of related tissues. The components of transcriptional regulatory networks (TRNs), including transcription factors, their cognate cis-regulatory motifs, and target PCGs/lncRNAs that determine developmental programs of different tissues/organs, are identified. Several candidate tissue-specific and abiotic stress-responsive transcripts associated with quantitative trait loci that determine important agronomic traits are also identified. These results provide an important resource to advance functional/translational genomic and genetic studies during chickpea development and environmental conditions.

A full-length transcriptome and expression atlas of protein-coding genes and long non-coding RNAs is generated in chickpea. Components of transcriptional regulatory networks and candidate tissue-specific transcripts associated with quantitative trait loci are identified.

Gene activation via Cre/lox-mediated excision in cowpea (Vigna unguiculata)

Expression of Cre recombinase by AtRps5a_pro or AtDD45_pro enabled Cre/lox-mediated recombination at an early embryonic developmental stage upon crossing, activating transgenes in the hybrid cowpea and tobacco. Genetic engineering ideally results in precise spatiotemporal control of transgene expression. To activate transgenes exclusively in a hybrid upon fertilization, we evaluated a Cre/lox-mediated gene activation system with the Cre recombinase expressed by either AtRps5a or AtDD45 promoters that showed activity in egg cells and young embryos. In crosses between Cre recombinase lines and transgenic lines harboring a lox-excision reporter cassette with ZsGreen driven by the AtUbq3 promoter after Cre/lox-mediated recombination, we observed complete excision of the lox-flanked intervening DNA sequence between the AtUbq3_pro and the ZsGreen coding sequence in F₁ progeny upon genotyping but no ZsGreen expression in F₁ seeds or seedlings. The incapability to observe ZsGreen fluorescence was attributed to the activity of the AtUbq3_pro. Strong ZsGreen expression in F₁ seeds was observed after recombination when ZsGreen was driven by the AtUbq10 promoter. Using the AtDD45_pro to express Cre resulted in more variation in recombination frequencies between transgenic lines and crosses. Regardless of the promoter used to regulate Cre, mosaic F₁ progeny were rare, suggesting gene activation at an early embryo-developmental stage. Observation of ZsGreen-expressing tobacco embryos at the globular stage from crosses with the AtRps5a_proCre lines pollinated by the AtUbq3_prolox line supported the early activation mode.

Transcriptome atlas of Phalaenopsis equestris

The vast diversity of Orchidaceae together with sophisticated adaptations to pollinators and other unique features make this family an attractive model for evolutionary and functional studies. The sequenced genome of Phalaenopsis equestris facilitates Orchidaceae research. Here, we present an RNA-seq-based transcriptome map of P. equestris that covers 19 organs of the plant, including leaves, roots, floral organs and the shoot apical meristem. We demonstrated the high quality of the data and showed the similarity of the P. equestris transcriptome map with the gene expression atlases of other plants. The transcriptome map can be easily accessed through our database Transcriptome Variation Analysis (TraVA) for visualizing gene expression profiles. As an example of the application, we analyzed the expression of Phalaenopsis “orphan” genes–those that do not have recognizable similarity with the genes of other plants. We found that approximately half of these genes were not expressed; the ones that were expressed were predominantly expressed in reproductive structures.

Constraints and Prospects of Improving Cowpea Productivity to Ensure Food, Nutritional Security and Environmental Sustainability

Providing safe and secure food for an increasing number of people globally is challenging. Coping with such a human population by merely applying the conventional agricultural production system has not proved to be agro-ecologically friendly; nor is it sustainable. Cowpea (Vigna unguiculata (L) Walp) is a multi-purpose legume. It consists of high-quality protein for human consumption, and it is rich in protein for livestock fodder. It enriches the soil in that it recycles nutrients through the fixation of nitrogen in association with nodulating bacteria. However, the productivity of this multi-functional, indigenous legume that is of great value to African smallholder farmers and the rural populace, and also to urban consumers and entrepreneurs, is limited. Because cowpea is of strategic importance in Africa, there is a need to improve on its productivity. Such endeavors in Africa are wrought with challenges that include drought, salinity, the excessive demand among farmers for synthetic chemicals, the repercussions of climate change, declining soil nutrients, microbial infestations, pest issues, and so forth. Nevertheless, giant strides have already been made and there have already been improvements in adopting sustainable and smart biotechnological approaches that are favorably influencing the production costs of cowpea and its availability. As such, the prospects for a leap in cowpea productivity in Africa and in the enhancement of its genetic gain are good. Potential and viable means for overcoming some of the above-mentioned production constraints would be to focus on the key cowpea producer nations in Africa and to encourage them to embrace biotechnological techniques in an integrated approach to enhance for sustainable productivity. This review highlights the spectrum of constraints that limit the cowpea yield, but looks ahead of the constraints and seeks a way forward to improve cowpea productivity in Africa. More importantly, this review investigates applications and insights concerning mechanisms of action for implementing eco-friendly biotechnological techniques, such as the deployment of bio inoculants, applying climate-smart agricultural (CSA) practices, agricultural conservation techniques, and multi-omics smart technology in the spheres of genomics, transcriptomics, proteomics, and metabolomics, for improving cowpea yields and productivity to achieve sustainable agro-ecosystems, and ensuring their stability.

Genetic, anatomical, and environmental patterns related to pod shattering resistance in domesticated cowpea [Vigna unguiculata (L.) Walp]

Pod shattering, which causes the explosive release of seeds from the pod, is one of the main sources of yield losses in cowpea in arid and semi-arid areas. Reduction of shattering has therefore been a primary target for selection during domestication and improvement of cowpea, among other species. Using a mini-core diversity panel of 368 cowpea accessions, four regions with a statistically significant association with pod shattering were identified. Two genes (Vigun03g321100 and Vigun11g100600), involved in cell wall biosynthesis, were identified as strong candidates for pod shattering. Microscopic analysis was conducted on a subset of accessions representing the full spectrum of shattering phenotypes. This analysis indicated that the extent of wall fiber deposition was highly correlated with shattering. The results from this study also demonstrate that pod shattering in cowpea is exacerbated by arid environmental conditions. Finally, using a subset of West African landraces, patterns of historical selection for shattering resistance related to precipitation in the environment of origin were identified. Together, these results shed light on sources of resistance to pod shattering, which will, in turn, improve climate resilience of a major global nutritional staple.

Fab Advances in Fabaceae for Abiotic Stress Resilience: From 'Omics' to Artificial Intelligence

Legumes are a better source of proteins and are richer in diverse micronutrients over the nutritional profile of widely consumed cereals. However, when exposed to a diverse range of abiotic stresses, their overall productivity and quality are hugely impacted. Our limited understanding of genetic determinants and novel variants associated with the abiotic stress response in food legume crops restricts its amelioration. Therefore, it is imperative to understand different molecular approaches in food legume crops that can be utilized in crop improvement programs to minimize the economic loss. 'Omics'-based molecular breeding provides better opportunities over conventional breeding for diversifying the natural germplasm together with improving yield and quality parameters. Due to molecular advancements, the technique is now equipped with novel 'omics' approaches such as ionomics, epigenomics, fluxomics, RNomics, glycomics, glycoproteomics, phosphoproteomics, lipidomics, regulomics, and secretomics. Pan-omics-which utilizes the molecular bases of the stress response to identify genes (genomics), mRNAs (transcriptomics), proteins (proteomics), and biomolecules (metabolomics) associated with stress regulation-has been widely used for abiotic stress amelioration in food legume crops. Integration of pan-omics with novel omics approaches will fast-track legume breeding programs. Moreover, artificial intelligence (AI)-based algorithms can be utilized for simulating crop yield under changing environments, which can help in predicting the genetic gain beforehand. Application of machine learning (ML) in quantitative trait loci (QTL) mining will further help in determining the genetic determinants of abiotic stress tolerance in pulses.

Developing a rapid and highly efficient cowpea regeneration, transformation and genome editing system using embryonic axis explants

Cowpea (Vigna unguiculata (L.) Walp.) is one of the most important legume crops planted worldwide, but despite decades of effort, cowpea transformation is still challenging due to inefficient Agrobacterium-mediated transfer DNA delivery, transgenic selection and in vitro shoot regeneration. Here, we report a highly efficient transformation system using embryonic axis explants isolated from imbibed mature seeds. We found that removal of the shoot apical meristem from the explants stimulated direct multiple shoot organogenesis from the cotyledonary node tissue. The application of a previously reported ternary transformation vector system provided efficient Agrobacterium-mediated gene delivery, while the utilization of spcN as selectable marker enabled more robust transgenic selection, plant recovery and transgenic plant generation without escapes and chimera formation. Transgenic cowpea plantlets developed exclusively from the cotyledonary nodes at frequencies of 4% to 37% across a wide range of cowpea genotypes. CRISPR/Cas-mediated gene editing was successfully demonstrated. The transformation principles established here could also be applied to other legumes to increase transformation efficiencies.

Developing a rapid and highly efficient cowpea regeneration, transformation and genome editing system using embryonic axis explants

Cowpea (Vigna unguiculata (L.) Walp.) is one of the most important legume crops planted worldwide, but despite decades of effort, cowpea transformation is still challenging due to inefficient Agrobacterium-mediated transfer DNA delivery, transgenic selection and in vitro shoot regeneration. Here, we report a highly efficient transformation system using embryonic axis explants isolated from imbibed mature seeds. We found that removal of the shoot apical meristem from the explants stimulated direct multiple shoot organogenesis from the cotyledonary node tissue. The application of a previously reported ternary transformation vector system provided efficient Agrobacterium-mediated gene delivery, while the utilization of spcN as selectable marker enabled more robust transgenic selection, plant recovery and transgenic plant generation without escapes and chimera formation. Transgenic cowpea plantlets developed exclusively from the cotyledonary nodes at frequencies of 4% to 37% across a wide range of cowpea genotypes. CRISPR/Cas-mediated gene editing was successfully demonstrated. The transformation principles established here could also be applied to other legumes to increase transformation efficiencies.

Arachis hypogaea gene expression atlas for fastigiata subspecies of cultivated groundnut to accelerate functional and translational genomics applications

Spatio-temporal and developmental stage-specific transcriptome analysis plays a crucial role in systems biology-based improvement of any species. In this context, we report here the Arachis hypogaea gene expression atlas (AhGEA) for the world's widest cultivated subsp. fastigiata based on RNA-seq data using 20 diverse tissues across five key developmental stages. Approximately 480 million paired-end filtered reads were generated followed by identification of 81 901 transcripts from an early-maturing, high-yielding, drought-tolerant groundnut variety, ICGV 91114. Further, 57 344 genome-wide transcripts were identified with ≥1 FPKM across different tissues and stages. Our in-depth analysis of the global transcriptome sheds light into complex regulatory networks namely gravitropism and photomorphogenesis, seed development, allergens and oil biosynthesis in groundnut. Importantly, interesting insights into molecular basis of seed development and nodulation have immense potential for translational genomics research. We have also identified a set of stable expressing transcripts across the selected tissues, which could be utilized as internal controls in groundnut functional genomics studies. The AhGEA revealed potential transcripts associated with allergens, which upon appropriate validation could be deployed in the coming years to develop consumer-friendly groundnut varieties. Taken together, the AhGEA touches upon various important and key features of cultivated groundnut and provides a reference for further functional, comparative and translational genomics research for various economically important traits.

Identification of QTL for perenniality and floral scent in cowpea (Vigna unguiculata [L.] Walp.)

Perennial habit and floral scent are major traits that distinguish domesticated cowpeas from their wild relatives. However, the genetic basis of these two important traits remains largely unknown in cowpea. Plant longevity, a perenniality-related trait, and floral scent, an outcrossing trait, were investigated using a RIL population derived from a cross between a domesticated and a wild cowpea. QTL analysis revealed three significant loci, one on chromosome 8 associated with plant longevity and two, on chromosomes 1 and 11, for floral scent. Genes within the QTL regions were identified. Genes encoding an F-box protein (Vigun08g215300) and two kinases (Vigun08g217000, Vigun08g217800), and involved in physiological processes including regulation of flowering time and plant longevity, were identified within the perenniality QTL region. A cluster of O-methyltransferase genes (Vigun11g096800, Vigun11g096900, Vigun11g097000, Vigun11g097600, and Vigun11g097800) was identified within the floral scent QTL region. These O-methyltransferase cowpea genes are orthologs of the Arabidopsis N-acetylserotonin O-methyltransferase (ASMT) gene, which is involved in the biosynthesis of melatonin. Melatonin is an indole derivative, which is an essential molecule for plant interactions with pollinators. These findings lay the foundation for further exploration of the genetic mechanisms of perenniality and floral scent in cowpea. Knowledge from this study can help in the development of new extended-growth cycle lines with increased yield or lines with increased outcrossing for population breeding.

Seed Coat Pattern QTL and Development in Cowpea (Vigna unguiculata [L.] Walp.)

The appearance of the seed is an important aspect of consumer preference for cowpea (Vigna unguiculata [L.] Walp.). Seed coat pattern in cowpea has been a subject of study for over a century. This study makes use of newly available resources, including mapping populations, a reference genome and additional genome assemblies, and a high-density single nucleotide polymorphism genotyping platform, to map various seed coat pattern traits to three loci, concurrent with the Color Factor (C), Watson (W), and Holstein (H) factors identified previously. Several gene models encoding proteins involved in regulating the later stages of the flavonoid biosynthesis pathway have been identified as candidate genes, including a basic helix-loop-helix gene (Vigun07g110700) for the C locus, a WD-repeat gene (Vigun09g139900) for the W locus and an E3 ubiquitin ligase gene (Vigun10g163900) for the H locus. A model of seed coat development, consisting of six distinct stages, is described to explain some of the observed pattern phenotypes.

Gene Expression Maps in Plants: Current State and Prospects

For many years, progress in the identification of gene functions has been based on classical genetic approaches. However, considerable recent omics developments have brought to the fore indirect but high-resolution methods of gene function identification such as transcriptomics, proteomics, and metabolomics. A transcriptome map is a powerful source of functional information and the result of the genome-wide expression analysis of a broad sampling of tissues and/or organs from different developmental stages and/or environmental conditions. In plant science, the application of transcriptome maps extends from the inference of gene regulatory networks to evolutionary studies. However, only some of these data have been integrated into databases, thus enabling analyses to be conducted without raw data; without this integration, extensive data preprocessing is required, which limits data usability. In this review, we summarize the state of plant transcriptome maps, analyze the problems associated with the combined analysis of large-scale data from various studies, and outline possible solutions to these problems.

The genome of cowpea (Vigna unguiculata [L.] Walp.)

Cowpea (Vigna unguiculata [L.] Walp.) is a major crop for worldwide food and nutritional security, especially in sub-Saharan Africa, that is resilient to hot and drought-prone environments. An assembly of the single-haplotype inbred genome of cowpea IT97K-499-35 was developed by exploiting the synergies between single-molecule real-time sequencing, optical and genetic mapping, and an assembly reconciliation algorithm. A total of 519 Mb is included in the assembled sequences. Nearly half of the assembled sequence is composed of repetitive elements, which are enriched within recombination-poor pericentromeric regions. A comparative analysis of these elements suggests that genome size differences between Vigna species are mainly attributable to changes in the amount of Gypsy retrotransposons. Conversely, genes are more abundant in more distal, high-recombination regions of the chromosomes; there appears to be more duplication of genes within the NBS-LRR and the SAUR-like auxin superfamilies compared with other warm-season legumes that have been sequenced. A surprising outcome is the identification of an inversion of 4.2 Mb among landraces and cultivars, which includes a gene that has been associated in other plants with interactions with the parasitic weed Striga gesnerioides. The genome sequence facilitated the identification of a putative syntelog for multiple organ gigantism in legumes. A revised numbering system has been adopted for cowpea chromosomes based on synteny with common bean (Phaseolus vulgaris). An estimate of nuclear genome size of 640.6 Mbp based on cytometry is presented.

Legume Cytosolic and Plastid Acetyl-Coenzyme-A Carboxylase Genes Differ by Evolutionary Patterns and Selection Pressure Schemes Acting before and after Whole-Genome Duplications

Acetyl-coenzyme A carboxylase (ACCase, E.C.6.4.1.2) catalyzes acetyl-coenzyme A carboxylation to malonyl coenzyme A. Plants possess two distinct ACCases differing by cellular compartment and function. Plastid ACCase contributes to de novo fatty acid synthesis, whereas cytosolic enzyme to the synthesis of very long chain fatty acids, phytoalexins, flavonoids, and anthocyanins. The narrow leafed lupin (Lupinus angustifolius L.) represents legumes, a plant family which evolved by whole-genome duplications (WGDs). The study aimed on the contribution of these WGDs to the multiplication of ACCase genes and their further evolutionary patterns. The molecular approach involved bacterial artificial chromosome (BAC) library screening, fluorescent in situ hybridization, linkage mapping, and BAC sequencing. In silico analysis encompassed sequence annotation, comparative mapping, selection pressure calculation, phylogenetic inference, and gene expression profiling. Among sequenced legumes, the highest number of ACCase genes was identified in lupin and soybean. The most abundant plastid ACCase subunit genes were accB. ACCase genes in legumes evolved by WGDs, evidenced by shared synteny and Bayesian phylogenetic inference. Transcriptional activity of almost all copies was confirmed. Gene duplicates were conserved by strong purifying selection, however, positive selection occurred in Arachis (accB2) and Lupinus (accC) lineages, putatively predating the WGD event(s). Early duplicated accA and accB genes underwent transcriptional sub-functionalization.

Identification of Candidate Genes Controlling Black Seed Coat and Pod Tip Color in Cowpea (Vigna unguiculata [L.] Walp)

Seed coat color is an important part of consumer preferences for cowpea (Vigna unguiculata [L.] Walp). Color has been studied in numerous crop species and has often been linked to loci controlling the anthocyanin biosynthesis pathway. This study makes use of available resources, including mapping populations, a reference genome, and a high-density single nucleotide polymorphism genotyping platform, to map the black seed coat and purple pod tip color traits, with the gene symbol Bl, in cowpea. Several gene models encoding MYB domain protein 113 were identified as candidate genes. MYB domain proteins have been shown in other species to control expression of genes encoding enzymes for the final steps in the anthocyanin biosynthesis pathway. PCR analysis indicated that a presence/absence variation of one or more MYB113 genes may control the presence or absence of black pigment. A PCR marker has been developed for the MYB113 gene Vigun05g039500, a candidate gene for black seed coat color in cowpea.

The RNA-Seq-based high resolution gene expression atlas of chickpea (Cicer arietinum L.) reveals dynamic spatio-temporal changes associated with growth and development

Chickpea is one of the world's largest cultivated food legumes and is an excellent source of high-quality protein to the human diet. Plant growth and development are controlled by programmed expression of a suite of genes at the given time, stage, and tissue. Understanding how the underlying genome sequence translates into specific plant phenotypes at key developmental stages, information on gene expression patterns is crucial. Here, we present a comprehensive Cicer arietinum Gene Expression Atlas (CaGEA) across different plant developmental stages and organs covering the entire life cycle of chickpea. One of the widely used drought tolerant cultivars, ICC 4958 has been used to generate RNA-Seq data from 27 samples at 5 major developmental stages of the plant. A total of 816 million raw reads were generated and of these, 794 million filtered reads after quality control (QC) were subjected to downstream analysis. A total of 15,947 unique number of differentially expressed genes across different pairwise tissue combinations were identified. Significant differences in gene expression patterns contributing in the process of flowering, nodulation, and seed and root development were inferred in this study. Furthermore, differentially expressed candidate genes from "QTL-hotspot" region associated with drought stress response in chickpea were validated.

Assembled genomic and tissue-specific transcriptomic data resources for two genetically distinct lines of Cowpea ( Vigna unguiculata (L.) Walp)

Cowpea ( Vigna unguiculata (L.) Walp) is an important legume crop for food security in areas of low-input and smallholder farming throughout Africa and Asia. Genetic improvements are required to increase yield and resilience to biotic and abiotic stress and to enhance cowpea crop performance. An integrated cowpea genomic and gene expression data resource has the potential to greatly accelerate breeding and the delivery of novel genetic traits for cowpea. Extensive genomic resources for cowpea have been absent from the public domain; however, a recent early release reference genome for IT97K-499-35 ( Vigna unguiculata v1.0, NSF, UCR, USAID, DOE-JGI, http://phytozome.jgi.doe.gov/) has now been established in a collaboration between the Joint Genome Institute (JGI) and University California (UC) Riverside. Here we release supporting genomic and transcriptomic data for two transformable cowpea varieties, IT97K-499-35 and IT86D-1010. The transcriptome resource includes six tissue-specific datasets for each variety, with particular emphasis on reproductive tissues that extend and support the V. unguiculata v1.0 reference. Annotations have been included in our resource to allow direct mapping to the v1.0 cowpea reference. The resource described here is supported by downloadable raw and assembled sequence data.

Assembled genomic and tissue-specific transcriptomic data resources for two genetically distinct lines of Cowpea ( Vigna unguiculata (L.) Walp)

Cowpea ( Vigna unguiculata (L.) Walp) is an important legume crop for food security in areas of low-input and smallholder farming throughout Africa and Asia. Genetic improvements are required to increase yield and resilience to biotic and abiotic stress and to enhance cowpea crop performance. An integrated cowpea genomic and gene expression data resource has the potential to greatly accelerate breeding and the delivery of novel genetic traits for cowpea. Extensive genomic resources for cowpea have been absent from the public domain; however, a recent early release reference genome for IT97K-499-35 ( Vigna unguiculata v1.0, NSF, UCR, USAID, DOE-JGI, http://phytozome.jgi.doe.gov/) has now been established in a collaboration between the Joint Genome Institute (JGI) and University California (UC) Riverside. Here we release supporting genomic and transcriptomic data for two transformable cowpea varieties, IT97K-499-35 and IT86D-1010. The transcriptome resource includes six tissue-specific datasets for each variety, with particular emphasis on reproductive tissues that extend and support the V. unguiculata v1.0 reference. Annotations have been included in our resource to allow direct mapping to the v1.0 cowpea reference. The resource described here is supported by downloadable raw and assembled sequence data.

QTL mapping and transcriptome analysis of cowpea reveals candidate genes for root-knot nematode resistance

Cowpea is one of the most important food and forage legumes in drier regions of the tropics and subtropics. However, cowpea yield worldwide is markedly below the known potential due to abiotic and biotic stresses, including parasitism by root-knot nematodes (Meloidogyne spp., RKN). Two resistance genes with dominant effect, Rk and Rk2, have been reported to provide resistance against RKN in cowpea. Despite their description and use in breeding for resistance to RKN and particularly genetic mapping of the Rk locus, the exact genes conferring resistance to RKN remain unknown. In the present work, QTL mapping using recombinant inbred line (RIL) population 524B x IT84S-2049 segregating for a newly mapped locus and analysis of the transcriptome changes in two cowpea near-isogenic lines (NIL) were used to identify candidate genes for Rk and the newly mapped locus. A major QTL, designated QRk-vu9.1, associated with resistance to Meloidogyne javanica reproduction, was detected and mapped on linkage group LG9 at position 13.37 cM using egg production data. Transcriptome analysis on resistant and susceptible NILs 3 and 9 days after inoculation revealed up-regulation of 109 and 98 genes and down-regulation of 110 and 89 genes, respectively, out of 19,922 unique genes mapped to the common bean reference genome. Among the differentially expressed genes, four and nine genes were found within the QRk-vu9.1 and QRk-vu11.1 QTL intervals, respectively. Six of these genes belong to the TIR-NBS-LRR family of resistance genes and three were upregulated at one or more time-points. Quantitative RT-PCR validated gene expression to be positively correlated with RNA-seq expression pattern for eight genes. Future functional analysis of these cowpea genes will enhance our understanding of Rk-mediated resistance and identify the specific gene responsible for the resistance.

A compendium of transcription factor and Transcriptionally active protein coding gene families in cowpea (Vigna unguiculata L.)

BACKGROUND

Cowpea (Vigna unguiculata (L.) Walp.) is the most important food and forage legume in the semi-arid tropics of sub-Saharan Africa where approximately 80% of worldwide production takes place primarily on low-input, subsistence farm sites. Among the major goals of cowpea breeding and improvement programs are the rapid manipulation of agronomic traits for seed size and quality and improved resistance to abiotic and biotic stresses to enhance productivity. Knowing the suite of transcription factors (TFs) and transcriptionally active proteins (TAPs) that control various critical plant cellular processes would contribute tremendously to these improvement aims.

RESULTS

We used a computational approach that employed three different predictive pipelines to data mine the cowpea genome and identified over 4400 genes representing 136 different TF and TAP families. We compare the information content of cowpea to two evolutionarily close species common bean (Phaseolus vulgaris), and soybean (Glycine max) to gauge the relative informational content. Our data indicate that correcting for genome size cowpea has fewer TF and TAP genes than common bean (4408 / 5291) and soybean (4408/ 11,065). Members of the GROWTH-REGULATING FACTOR (GRF) and Auxin/indole-3-acetic acid (Aux/IAA) gene families appear to be over-represented in the genome relative to common bean and soybean, whereas members of the MADS (Minichromosome maintenance deficient 1 (MCM1), AGAMOUS, DEFICIENS, and serum response factor (SRF)) and C2C2-YABBY appear to be under-represented. Analysis of the AP2-EREBP APETALA2-Ethylene Responsive Element Binding Protein (AP2-EREBP), NAC (NAM (no apical meristem), ATAF1, 2 (Arabidopsis transcription activation factor), CUC (cup-shaped cotyledon)), and WRKY families, known to be important in defense signaling, revealed changes and phylogenetic rearrangements relative to common bean and soybean that suggest these groups may have evolved different functions.

CONCLUSIONS

The availability of detailed information on the coding capacity of the cowpea genome and in particular the various TF and TAP gene families will facilitate future comparative analysis and development of strategies for controlling growth, differentiation, and abiotic and biotic stress resistances of cowpea.