Comprehensive genomic analyses of Vigna unguiculata provide insights into population differentiation and the genetic basis of key agricultural traits.

Vigna unguiculata is an important legume crop worldwide. The subsp. sesquipedalis and unguiculata are the two major types grown; the former is mainly grown in Asia to produce fresh pods, while the latter is mainly grown in Africa to produce seeds. Here, a chromosome-scale genome for subsp. sesquipedalis was generated by combining high-fidelity (HiFi) long-read sequencing with high-throughput chromosome conformation capture (Hi-C) technology. The genome size for all contigs and N50 were 594 and 18.5 Mb, respectively. The Hi-C interaction map helped cluster 91% of the contigs into 11 chromosomes. Genome comparisons between subsp. sesquipedalis and unguiculata revealed extensive genomic variations, and some variations resulted in gene loss. A germplasm panel with 315 accessions of V. unguiculata was resequenced, and a genomic variation map was constructed. Population structure and phylogenetic analyses suggested that subsp. sesquipedalis originated from subsp. unguiculata. Highly differentiated genomic regions were also identified, and a number of genes functionally enriched in adaptations were located in these regions. Two traits, pod length (PL) and pod width (PW), were observed for this germplasm, and genome-wide association analysis of these traits was performed. The quantitative trait loci (QTLs) for these two traits were identified, and their candidate genes were uncovered. Interestingly, genomic regions of PL QTLs also showed strong signals of artificial selection. Taken together, the results of this study provide novel insights into the population differentiation and genetic basis of key agricultural traits in V. unguiculata.

Insights into Temperature and Hypoxia Tolerance in Cowpea Weevil via HIF-1.

Cowpea weevil (Callosobruchus maculatus) is a major pest that leads to severe damage of the stored leguminous grains. Several management approaches, including physical barriers, biological or chemical methods, are used for controlling bruchid in cowpea. These methods usually target the metabolically active state of weevil. However, it becomes less effective at early stages as egg, larva, or pupa under low temperature and oxygen conditions. Since hypoxia-inducible factor-1 (HIF-1) is known to coordinate multiple gene responses to low oxygen or low temperature signals, we examined the HIF-1α gene expression under low temperature and hypoxic treatments. At -20 °C, it took 4 h to reduce the survival rate for eggs, larvae, and pupae down to 10%, while at 4 °C and 15 °C, the survival rate remained higher than 50% even after 128 h as HIF-1α gene expression peaked at 15 °C. Moreover, HIF-1 protein offers a valuable target for early stage pest control complementary to traditional methods. In particular, HIF-1 inhibitor camptothecin (CPT), one of the five HIF-1 inhibitors examined, achieved a very significant reduction of 96.2% and 95.5% relative to the control in weevil survival rate into adult at 4 °C and 30 °C, respectively. Our study can be used as one model system for drug development in virus infections and human cancer.

Transcriptomic profiling and analysis of differentially expressed genes in asparagus bean (Vigna unguiculata ssp. sesquipedalis) under salt stress.

Asparagus bean (Vigna unguiculata ssp. sesquipedalis) is a warm season legume which is widely distributed over subtropical regions and semiarid areas. It is mainly grown as a significant protein source in developing countries. Salinity, as one of the main abiotic stress factors, constrains the normal growth and yield of asparagus bean. This study used two cultivars (a salt-sensitive genotype and a salt-tolerant genotype) under salt stress vs. control to identify salt-stress-induced genes in asparagus bean using RNA sequencing. A total of 692,086,838 high-quality clean reads, assigned to 121,138 unigenes, were obtained from control and salt-treated libraries. Then, 216 root-derived DEGs (differentially expressed genes) and 127 leaf-derived DEGs were identified under salt stress between the two cultivars. Of these DEGs, thirteen were assigned to six transcription factors (TFs), including AP2/EREBP, CCHC(Zn), C2H2, WRKY, WD40-like and LIM. GO analysis indicated four DEGs might take effects on the "oxidation reduction", "transport" and "signal transduction" process. Moreover, expression of nine randomly-chosen DEGs was verified by quantitative real-time-PCR (qRT-PCR) analysis. Predicted function of the nine tested DEGs was mainly involved in the KEGG pathway of cation transport, response to osmotic stress, and phosphorelay signal transduction system. A salt-stress-related pathway of "SNARE interactions in vesicular transport" was concerned. As byproducts, 15, 321 microsatellite markers were found in all the unigenes, and 17 SNP linked to six salt-stress induced DEGs were revealed. These candidate genes provide novel insights for understanding the salt tolerance mechanism of asparagus bean in the future.

The genome assembly of asparagus bean, Vigna unguiculata ssp. sesquipedialis.

Asparagus bean (Vigna. unguiculata ssp. sesquipedialis), known for its very long and tender green pods, is an important vegetable crop broadly grown in the developing Asian countries. In this study, we reported a 632.8 Mb assembly (549.81 Mb non-N size) of asparagus bean based on the whole genome shotgun sequencing strategy. We also generated a linkage map for asparagus bean, which helped anchor 94.42% of the scaffolds into 11 pseudo-chromosomes. A total of 42,609 protein-coding genes and 3,579 non-protein-coding genes were predicted from the assembly. Taken together, these genomic resources of asparagus bean will help develop a pan-genome of V. unguiculata and facilitate the investigation of economically valuable traits in this species, so that the cultivation of this plant would help combat the protein and energy malnutrition in the developing world.

A High Density Genetic Map Derived from RAD Sequencing and Its Application in QTL Analysis of Yield-Related Traits in Vigna unguiculata.

Cowpea [Vigna unguiculata (L.) Walp.] is an annual legume of economic importance and widely grown in the semi-arid tropics. However, high-density genetic maps of cowpea are still lacking. Here, we identified 34,868 SNPs (single nucleotide polymorphisms) that were distributed in the cowpea genome based on the RAD sequencing (restriction-site associated DNA sequencing) technique using a population of 170 individuals (two cowpea parents and 168 F2:3 progenies). Of these, 17,996 reliable SNPs were allotted to 11 consensus linkage groups (LGs). The length of the genetic map was 1,194.25 cM in total with a mean distance of 0.066 cM/SNP marker locus. Using this map and the F2:3 population, combined with the CIM (composite interval mapping) method, eleven quantitative trait loci (QTL) of yield-related trait were detected on seven LGs (LG4, 5, 6, 7, 9, 10, and 11) in cowpea. These QTL explained 0.05-17.32% of the total phenotypic variation. Among these, four QTL were for pod length, four QTL for thousand-grain weight (TGW), two QTL for grain number per pod, and one QTL for carpopodium length. Our results will provide a foundation for understanding genes related to grain yield in the cowpea and genus Vigna.

Network Analysis Reveals the Recognition Mechanism for Dimer Formation of Bulb-type Lectins.

The bulb-type lectins are proteins consist of three sequential beta-sheet subdomains that bind to specific carbohydrates to perform certain biological functions. The active states of most bulb-type lectins are dimeric and it is thus important to elucidate the short- and long-range recognition mechanism for this dimer formation. To do so, we perform comparative sequence analysis for the single- and double-domain bulb-type lectins abundant in plant genomes. In contrast to the dimer complex of two single-domain lectins formed via protein-protein interactions, the double-domain lectin fuses two single-domain proteins into one protein with a short linker and requires only short-range interactions because its two single domains are always in close proximity. Sequence analysis demonstrates that the highly variable but coevolving polar residues at the interface of dimeric bulb-type lectins are largely absent in the double-domain bulb-type lectins. Moreover, network analysis on bulb-type lectin proteins show that these same polar residues have high closeness scores and thus serve as hubs with strong connections to all other residues. Taken together, we propose a potential mechanism for this lectin complex formation where coevolving polar residues of high closeness are responsible for long-range recognition.

Development of 12 chloroplast microsatellite markers in Vigna unguiculata (Fabaceae) and amplification in Phaseolus vulgaris.

PREMISE OF THE STUDY: Vigna unguiculata is an economically important legume, and the complexity of its variability and evolution needs to be further understood. Based on publicly available databases, we developed chloroplast microsatellite primers to investigate genetic diversity within V. unguiculata and its related species Phaseolus vulgaris. • METHODS AND RESULTS: Twelve polymorphic chloroplast microsatellite markers were developed and characterized in 62 V. unguiculata individuals. The number of alleles per locus varied between two and four, the unbiased haploid diversity per locus ranged from 0.123 to 0.497, and the polymorphism information content varied from 0.114 to 0.369. In cross-species amplifications, nine of these markers showed polymorphism in 29 P. vulgaris individuals. • CONCLUSIONS: The newly developed chloroplast microsatellite markers exhibit variation in V. unguiculata as well as their transferability in P. vulgaris. These markers can be used to investigate genetic diversity and evolution in V. unguiculata and P. vulgaris.

Development and characterization of polymorphic microsatellite markers in Momordica charantia (Cucurbitaceae).

PREMISE OF THE STUDY: Polymorphic microsatellite markers were developed for Momordica charantia L. to investigate the genetic diversity and population structure within and between M. charantia and its four related species (Cucurbita pepo L., Luffa cylindrical L., Lagenaria siceraria L., and Cucumis sativus L.). • METHODS AND RESULTS: Using the Fast Isolation by AFLP of Sequence COntaining Repeats (FIASCO) method, 16 polymorphic microsatellite loci were identified in 36 individuals of M. charantia. Across all the M. charantia samples, the number of alleles per locus ranged from three to eight. Seven primers successfully amplified in the four related species. • CONCLUSIONS: These markers will be useful to study population ecology and population differentiation among M. charantia species and its related species.

Genetic analysis of salt stress responses in asparagus bean (Vigna unguiculata (L.) ssp. sesquipedalis Verdc.).

Salt stress responses of 23 asparagus bean cultivars were evaluated using 14-day-old seedlings after 15-day exposure to 75 mM NaCl in a hydroponics culture system. Salt-induced changes in plant growth and morphology, photosynthetic capacity, cell membrane integrity, and cellular protection enzyme systems as well as other physiological and biochemical traits were investigated to identify genotypic variability in salt response. This study also analyzed heredity parameters and correlation of the salt response index (SRI). Salt stress suppressed seedling growth and simultaneously reduced leaf area, content of chlorophyll and soluble proteins, net assimilation rate, as well as dry matter accumulation. In contrast, leaf blade cell membrane relative permeability, content of malondialdehyde, and antioxidant enzyme activity were elevated after salt treatment. Analysis of the heredity parameters has identified that the 18 investigated traits have different genotypic variance of the SRI values. Based on Ward's distances estimated for the sum of square variance in the SRI values, all the cultivars were classified into 2 discrete salt-tolerant and salt-sensitive clusters. Findings from this study will provide theoretical bases for identification and breeding of salt-tolerant cultivars in asparagus bean.