Synteny mapping between common bean and soybean reveals extensive blocks of shared loci

Artificial selection of mutations in two nearby genes gave rise to shattering resistance in soybean

Resistance to pod shattering is a key domestication-related trait selected for seed production in many crops. Here, we show that the transition from shattering in wild soybeans to shattering resistance in cultivated soybeans resulted from selection of mutations within the coding sequences of two nearby genes - Sh1 and Pdh1. Sh1 encodes a C2H2-like zinc finger transcription factor that promotes shattering by repressing SHAT1-5 expression, thereby reducing the secondary wall thickness of fiber cap cells in the abscission layers of pod sutures, while Pdh1 encodes a dirigent protein that orchestrates asymmetric lignin distribution in inner sclerenchyma, creating torsion in pod walls that facilitates shattering. Integration analyses of quantitative trait locus mapping, genome-wide association studies, and allele distribution in representative soybean germplasm suggest that these two genes are primary modulators underlying this domestication trait. Our study thus provides comprehensive understanding regarding the genetic, molecular, and cellular bases of shattering resistance in soybeans.

Fine-mapping and evolutionary history of R-BPMV, a dominant resistance gene to Bean pod mottle virus in Phaseolus vulgaris L

KEY MESSAGE

R-BPMV is located within a recently expanded TNL cluster in the Phaseolus genus with suppressed recombination and known for resistance to multiple pathogens including potyviruses controlled by the I gene. Bean pod mottle virus (BPMV) is a comovirus that infects common bean and legumes in general. BPMV is distributed throughout the world and is a major threat on soybean, a closely related species of common bean. In common bean, BAT93 was reported to carry the R-BPMV resistance gene conferring resistance to BPMV and linked with the I resistance gene. To fine map R-BPMV, 182 recombinant inbred lines (RILs) derived from the cross BAT93 × JaloEEP558 were genotyped with polymerase chain reaction (PCR)-based markers developed using genome assemblies from G19833 and BAT93, as well as BAT93 BAC clone sequences. Analysis of RILs carrying key recombination events positioned R-BPMV to a target region containing at least 16 TIR-NB-LRR (TNL) sequences in BAT93. Because the I cluster presents a suppression of recombination and a large number of repeated sequences, none of the 16 TNLs could be excluded as R-BPMV candidate gene. The evolutionary history of the TNLs for the I cluster were reconstructed using microsynteny and phylogenetic analyses within the legume family. A single I TNL was present in Medicago truncatula and lost in soybean, mirroring the absence of complete BPMV resistance in soybean. Amplification of TNLs in the I cluster predates the divergence of the Phaseolus species, in agreement with the emergence of R-BPMV before the separation of the common bean wild centers of diversity. This analysis provides PCR-based markers useful in marker-assisted selection (MAS) and laid the foundation for cloning of R-BPMV resistance gene in order to transfer the resistance into soybean.

The genome-wide characterization of WOX gene family in Phaseolus vulgaris L. during salt stress

The WUSCHEL-Related Homeobox (WOX) family is a type of homeobox transcription factor superfamily and its members perform many functions ranging from plant embryonic growth to organ formation in plants. Although the WOX proteins have been identified and characterized in many plant species, genome-wide identification and characterization of WOX proteins in the Phaseolus vulgaris genome has been performed for the first time in this study. Accordingly, 18 WOX proteins were identified using bioinformatics tools and biochemical/physicochemical properties of these proteins were investigated. Phvul-WOX genes were found to be categorized into three major phylogenetic groups according to the phylogenetic analysis and a total of five segmental duplication events were detected after duplication analysis. Moreover, the Phvul-WOX genes were found to be expressed in different plant tissues at different levels and some stress-related miRNAs have been found to target the Phvul-WOX genes based on miRNA analysis. Additionaly, MDA content, total protein level and catalase enzyme activity analyses were conducted in two P. vulgaris cultivars namely Yakutiye cv. and Zulbiye cv. subjected to 150 mM salt stress. Next, these cultivars were used for screening the expression levels of Phvul-WOX-1, Phvul-WOX-9, Phvul-WOX-11, Phvul-WOX-15 and Phvul-WOX-16 genes in response to salt stress. The insights gained from this study may be of assistance to the researchers who work in this area.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-022-01208-1.

Phylogenetic analysis of phytochrome A gene from Lablab purpureus (L.) Sweet

BACKGROUND

Phytochromes are the best characterized photoreceptors that perceive Red (R)/Far-Red (FR) signals and mediate key developmental responses in plants. It is well established that photoperiodic control of flowering is regulated by PHY A (phytochrome A) gene. So far, the members of PHY A gene family remains unexplored in Lablab purpureus, and therefore, their functions are still not deciphered. PHYA3 is the homologue of phytochrome A and known to be involved in dominant suppression of flowering under long day conditions by downregulating florigens in Glycine max. The present study is the first effort to identify and characterize any photoreceptor gene (PHYA3, in this study) in Lablab purpureus and decipher its phylogeny with related legumes.

RESULTS

PHYA3 was amplified in Lablab purpureus cv GNIB-21 (photo-insensitive and determinate) by utilizing primers designed from GmPHYA3 locus of Glycine max. This study was successful in partially characterizing PHYA3 in Lablab purpureus (LprPHYA3) which is 2 kb longer and belongs to exon 1 region of PHYA3 gene. Phylogenetic analysis of the nucleotide and protein sequences of PHYA genes through MEGA X delineated the conservation and evolution of Lablab purpureus PHYA3 (LprPHYA3) probably from PHYA genes of Vigna unguiculata, Glycine max and Vigna angularis. A conserved basic helix-loop-helix motif bHLH69 was predicted having DNA binding property. Domain analysis of GmPHYA protein and predicted partial protein sequence corresponding to exon-1 of LprPHYA3 revealed the presence of conserved domains (GAF and PAS domains) in Lablab purpureus similar to Glycine max.

CONCLUSION

Partial characterization of LprPHYA3 would facilitate the identification of complete gene in Lablab purpureus utilizing sequence information from phylogenetically related species of Fabaceae. This would allow screening of allelic variants for LprPHYA3 locus and their role in photoperiod responsive flowering. The present study could aid in modulating photoperiod responsive flowering in Lablab purpureus and other related legumes in near future through genome editing.

Genomics Associated Interventions for Heat Stress Tolerance in Cool Season Adapted Grain Legumes

Cool season grain legumes occupy an important place among the agricultural crops and essentially provide multiple benefits including food supply, nutrition security, soil fertility improvement and revenue for farmers all over the world. However, owing to climate change, the average temperature is steadily rising, which negatively affects crop performance and limits their yield. Terminal heat stress that mainly occurred during grain development phases severely harms grain quality and weight in legumes adapted to the cool season, such as lentils, faba beans, chickpeas, field peas, etc. Although, traditional breeding approaches with advanced screening procedures have been employed to identify heat tolerant legume cultivars. Unfortunately, traditional breeding pipelines alone are no longer enough to meet global demands. Genomics-assisted interventions including new-generation sequencing technologies and genotyping platforms have facilitated the development of high-resolution molecular maps, QTL/gene discovery and marker-assisted introgression, thereby improving the efficiency in legumes breeding to develop stress-resilient varieties. Based on the current scenario, we attempted to review the intervention of genomics to decipher different components of tolerance to heat stress and future possibilities of using newly developed genomics-based interventions in cool season adapted grain legumes.

Identification and Characterization of SRS Genes in Phaseolus vulgaris Genome and Their Responses Under Salt Stress

SHI-Related Sequence (SRS) transcription factors comprise a protein family with important roles in growth and development. However, the genome-wide study of the SRS protein family has not yet been carried out in the common bean. For this reason, the SRS family has been characterized in depth at both gene and protein levels and several bioinformatics methods have been used. As a result, 10 SRS genes have been identified and their proteins have been phylogenetically categorized into three major groups within the common bean. By investigating duplications that play a major role in the development of gene families, 19 duplication events have been identified in the SRS family (18 segmental and 1 tandem). In addition, using available RNAseq data, comparative expression analysis of Pvul-SRS genes was performed and expression changes in Pvul-SRS-1, 2, 4, 6, 7, and 10 genes were observed under both salt and drought stress. Five Pvul-SRS genes were selected based on RNAseq data (Pvul-SRS-1, 2, 4, 6, and 10) and screened with RT-qPCR in two common bean cultivars (Yakutiye 'salt-resistant' and Zulbiye 'salt-susceptible' cv.). These genes also showed different levels of expression between two common bean cultivars under salt stress conditions and this may explain the responses of Pvul-SRS genes against abiotic stress. In summary, this work is the first study in which in silico identification and characterization of Pvul-SRS genes have been examined at gene expression level. The results could therefore provide the basis for future studies of functional characterization of Pvul-SRS genes.

Identification and Expression Analysis of the Genes Involved in the Raffinose Family Oligosaccharides Pathway of Phaseolus vulgaris and Glycine max

Raffinose family oligosaccharides (RFO) play an important role in plants but are also considered to be antinutritional factors. A profound understanding of the galactinol and RFO biosynthetic gene families and the expression patterns of the individual genes is a prerequisite for the sustainable reduction of the RFO content in the seeds, without compromising normal plant development and functioning. In this paper, an overview of the annotation and genetic structure of all galactinol- and RFO biosynthesis genes is given for soybean and common bean. In common bean, three galactinol synthase genes, two raffinose synthase genes and one stachyose synthase gene were identified for the first time. To discover the expression patterns of these genes in different tissues, two expression atlases have been created through re-analysis of publicly available RNA-seq data. De novo expression analysis through an RNA-seq study during seed development of three varieties of common bean gave more insight into the expression patterns of these genes during the seed development. The results of the expression analysis suggest that different classes of galactinol- and RFO synthase genes have tissue-specific expression patterns in soybean and common bean. With the obtained knowledge, important galactinol- and RFO synthase genes that specifically play a key role in the accumulation of RFOs in the seeds are identified. These candidate genes may play a pivotal role in reducing the RFO content in the seeds of important legumes which could improve the nutritional quality of these beans and would solve the discomforts associated with their consumption.

Characterization of Vascular plant One-Zinc finger (VOZ) in soybean (Glycine max and Glycine soja) and their expression analyses under drought condition

Vascular plant one-zinc-finger (VOZ) transcription factors regulate plant growth and development under drought conditions. Six VOZ transcription factors encoding genes exist in soybean genome (both in Glycine max and Glycine soja). Herein, GmVOZs and GsVOZs were identified through in silico analysis and characterized with different bioinformatics tools and expression analysis. Phylogenetic analysis classified VOZ genes in four groups. Sequence logos analysis among G. max and G. soja amino acid residues revealed higher conservation. Presence of stress related cis-elements in the upstream regions of GmVOZs and GsVOZs highlights their role in tolerance against abiotic stresses. The collinearity analysis identified 14 paralogous/orthologous gene pairs within and between G. max and G. soja. The Ka/Ks values showed that soybean VOZ genes underwent selection pressure with limited functional deviation arising from whole genome and segmental duplication. The GmVOZs and GsVOZs were found to express in roots and leaves at seedling stage. The qRT-PCR revealed that GmVOZs and GsVOZs transcripts can be regulated by abiotic stresses such as polyethylene glycol (PEG). The findings of this study will provide a reference to decipher physiological and molecular functions of VOZ genes in soybean.

Shoot-derived miR2111 controls legume root and nodule development

Legumes control their nodule numbers through the autoregulation of nodulation (AON). Rhizobia infection stimulates the production of root-derived CLE peptide hormones that are translocated to the shoot where they regulate a new signal. We used soybean to demonstrate that this shoot-derived signal is miR2111, which is transported via phloem to the root where it targets transcripts of Too Much Love (TML), a negative regulator of nodulation. Shoot perception of rhizobia-induced CLE peptides suppresses miR2111 expression, resulting in TML accumulation in roots and subsequent inhibition of nodule organogenesis. Feeding synthetic mature miR2111 via the petiole increased nodule numbers per plant. Likewise, elevating miR2111 availability by over-expression promoted nodulation, while target mimicry of TML induced the opposite effect on nodule development in wild-type plants and alleviated the supernodulating and stunted root growth phenotypes of AON-defective mutants. Additionally, in non-nodulating wild-type plants, ectopic expression of miR2111 significantly enhanced lateral root emergence with a decrease in lateral root length and average root diameter. In contrast, hairy roots constitutively expressing the target mimic construct exhibited reduced lateral root density. Overall, these findings demonstrate that miR2111 is both the critical shoot-to-root factor that positively regulates root nodule development and also acts to shape root system architecture.

Allelic characterization and protein structure analysis reveals the involvement of splice site mutation for growth habit differences in Lablab purpureus (L.) Sweet

BACKGROUND

Interrelationship between growth habit and flowering played a key role in the domestication history of pulses; however, the actual genes responsible for these traits have not been identified in Indian bean. Determinate growth habit is desirable due to its early flowering, photo-insensitivity, synchronous pod maturity, ease in manual harvesting and short crop duration. The present study aimed to identify, characterize and validate the gene responsible for growth habit by using a candidate gene approach coupled with sequencing, multiple sequence alignment, protein structure prediction and binding pocket analysis.

RESULTS

Terminal flowering locus was amplified from GPKH 120 (indeterminate) and GNIB-21 (determinate) using the primers designed from PvTFL1y locus of common bean. Gene prediction revealed that the length of the third and fourth exons differed between the two alleles. Allelic sequence comparison indicated a transition from guanine to adenine at the end of the third exon in GNIB 21. This splice site single-nucleotide polymorphism (SNP) was validated in germplasm lines by sequencing. Protein structure analysis indicated involvement of two binding pockets for interaction of terminal flowering locus (TFL) protein with other proteins.

CONCLUSION

The splice site SNP present at the end of the third exon of TFL locus is responsible for the transformation of shoot apical meristem into a reproductive fate in the determinate genotype GNIB 21. The splice site SNP leads to absence of 14 amino acids in mutant TFL protein of GNIB 21, rendering the protein non-functional. This deletion disturbed previously reported anion-binding pocket and secondary binding pocket due to displacement of small β-sheet away from an external loop. This finding may enable the modulation of growth habit in Indian bean and other pulse crops through genome editing.

Orthology and synteny analysis of receptor-like kinases "RLK" and receptor-like proteins "RLP" in legumes

BACKGROUND

Legume species are an important plant model because of their protein-rich physiology. The adaptability and productivity of legumes are limited by major biotic and abiotic stresses. Responses to these stresses directly involve plasma membrane receptor proteins known as receptor-like kinases and receptor-like proteins. Evaluating the homology relations among RLK and RLP for seven legume species, and exploring their presence among synteny blocks allow an increased understanding of evolutionary relations, physical position, and chromosomal distribution in related species and their shared roles in stress responses.

RESULTS

Typically, a high proportion of RLK and RLP legume proteins belong to orthologous clusters, which is confirmed in this study, where between 66 to 90% of the RLKs and RLPs per legume species were classified in orthologous clusters. One-third of the evaluated syntenic blocks had shared RLK/RLP genes among both legumes and non-legumes. Among the legumes, between 75 and 98% of the RLK/RLP were present in syntenic blocks. The distribution of chromosomal segments between Phaseolus vulgaris and Vigna unguiculata, two species that diverged ~ 8 mya, were highly similar. Among the RLK/RLP synteny clusters, seven experimentally validated resistance RLK/RLP genes were identified in syntenic blocks. The RLK resistant genes FLS2, BIR2, ERECTA, IOS1, and AtSERK1 from Arabidopsis and SLSERK1 from Solanum lycopersicum were present in different pairwise syntenic blocks among the legume species. Meanwhile, only the LYM1- RLP resistant gene from Arabidopsis shared a syntenic blocks with Glycine max.

CONCLUSIONS

The orthology analysis of the RLK and RLP suggests a dynamic evolution in the legume family, with between 66 to 85% of RLK and 83 to 88% of RLP belonging to orthologous clusters among the species evaluated. In fact, for the 10-species comparison, a lower number of singleton proteins were reported among RLP compared to RLK, suggesting that RLP positions are more physically conserved compared to RLK. The identification of RLK and RLP genes among the synteny blocks in legumes revealed multiple highly conserved syntenic blocks on multiple chromosomes. Additionally, the analysis suggests that P. vulgaris is an appropriate anchor species for comparative genomics among legumes.

The style and substance of plant flavonoid decoration; towards defining both structure and function

Over 8000 different flavonoids have been described and a considerable number of new flavonoid structures are being elucidated every year. The advent of metabolomics alongside the development of phytochemical genetics - wherein the genetic basis underlying the regulation of the levels of plant metabolites is determined - has provided a massive boost to such efforts. That said our understanding of the individual function(s) of the vast majority of the metabolites that constitute this important class of phytochemicals remains unknown. Here we review what is known concerning the major decorative modifications of flavonoids in plants, namely hydroxylation, glycosylation, methylation and acylation. Our major focus is with regard to the in planta function of these modified compounds, however, we also highlight the demonstrated bioactive roles which they possess. We additionally performed a comprehensive survey of the flavonoids listed in the KNApSAcK database in order to assess the frequency of occurrence of each type of flavonoid modification. We conclude that whilst considerable research has been carried out regarding the biological roles of flavonoids most studies to date have merely provided information on the compound class or sub-classes thereof as a whole with too little currently known on the specific role of individual metabolites. We, therefore, finally suggest a framework based on currently available tools by which the relative importance of the individual compounds can be assessed under various biological conditions in order to fill this knowledge-gap.

Pod indehiscence is a domestication and aridity resilience trait in common bean

Plant domestication has strongly modified crop morphology and development. Nevertheless, many crops continue to display atavistic characteristics that were advantageous to their wild ancestors but are deleterious under cultivation, such as pod dehiscence (PD). Here, we provide the first comprehensive assessment of the inheritance of PD in the common bean (Phaseolus vulgaris), a major domesticated grain legume. Using three methods to evaluate the PD phenotype, we identified multiple, unlinked genetic regions controlling PD in a biparental population and two diversity panels. Subsequently, we assessed patterns of orthology among these loci and those controlling the trait in other species. Our results show that different genes were selected in each domestication and ecogeographic race. A chromosome Pv03 dirigent-like gene, involved in lignin biosynthesis, showed a base-pair substitution that is associated with decreased PD. This haplotype may underlie the expansion of Mesoamerican domesticates into northern Mexico, where arid conditions promote PD. The rise in frequency of the decreased-PD haplotype may be a consequence of the markedly different fitness landscape imposed by domestication. Environmental dependency and genetic redundancy can explain the maintenance of atavistic traits under domestication.

High-resolution DNA methylome reveals that demethylation enhances adaptability to continuous cropping comprehensive stress in soybean

BACKGROUND

Continuous cropping stress involves such factors as biological barriers, allelopathic autotoxicity, deterioration of soil physicochemical properties, and soil fertility imbalance and is regarded as a kind of comprehensive stress limiting soybean yield and quality. Genomic DNA methylation is an important regulatory mechanism for plants to resist various environmental stresses. Therefore, it is especially worthwhile to reveal genomic methylation characteristics under stress and clarify the relationship between DNA methylation status and continuous cropping stress adaptability in soybean.

RESULTS

We generated a genome-wide map of cytosine methylation induced by this kind of comprehensive stress in a tolerant soybean variety (Kang Xian 2, KX2) and a sensitive variety (He Feng, HF55) using whole-genome bisulfite sequencing (WGBS) technology. The expression of DNA demethylase genes was detected using real-time quantitative PCR (qRT-PCR). The functions of differentially methylated genes (DMGs) involved in stress response in biochemical metabolism and genetic information transmission were further assessed based on Gene Ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. The results showed that genomic DNA demethylation was closely related to continuous cropping comprehensive stress adaptability in soybean, which was further verified by the increasing expression of DNA demethylases ROS1 and DML. The demethylation of mCpG and mCpHpG (mCpApG preferred) contexts was more critical, which mainly occurred in gene-regulatory regions at the whole-chromosome scale. Moreover, this kind of stress adaptability may be related to various stress responders generated through strengthened glucose catabolism and amino acid and fatty acid anabolism, as well as fidelity transmission of genetic information.

CONCLUSIONS

Genomic DNA demethylation was closely associated with continuous cropping comprehensive stress adaptability, highlighting the promising potential of screening continuous cropping-tolerant cultivars by DNA methylation index and further exploring the application of DNA demethylases in soybean breeding.

Genome wide association study discovers genomic regions involved in resistance to soybean cyst nematode (Heterodera glycines) in common bean

Common bean (Phaseolus vulgaris L.) is an important high protein crop grown worldwide. North Dakota and Minnesota are the largest producers of common beans in the USA, but crop production is threatened by soybean cyst nematode (SCN; Heterodera glycines Ichinohe) because most current cultivars are susceptible. Greenhouse screening data using SCN HG type 0 from 317 plant introductions (PI's) from the USDA core collection was used to conduct a genome wide association study (GWAS). These lines were divided into two subpopulations based on principal component analysis (Middle American vs. Andean). Phenotypic results based on the female index showed that accessions could be classified as highly resistant (21% and 27%), moderately resistant (51% and 48%), moderately susceptible (27% and 22%) and highly susceptible (1% and 3%) for Middle American and Andean gene pools, respectively. Mixed models with two principal components (PCs) and kinship matrix for Middle American genotypes and Andean genotypes were used in the GWAS analysis using 3,985 and 4,811 single nucleotide polymorphic (SNP) markers, respectively which were evenly distributed across all 11 chromosomes. Significant peaks on Pv07, and Pv11 in Middle American and on Pv07, Pv08, Pv09 and Pv11 in Andean group were found to be associated with SCN resistance. Homologs of soybean rhg1, a locus which confers resistance to SCN in soybean, were identified on chromosomes Pv01 and Pv08 in the Middle American and Andean gene pools, respectively. These genomic regions may be the key to develop SCN-resistant common bean cultivars.

Genome-wide Estimation of Evolutionary Distance and Phylogenetic Analysis of Homologous Genes

Homologous genes, including paralogs and orthologs, are genes that share sequence homologies within or between different species. Homologous genes originate from a common origin through speciation, genetic duplication or horizontal gene transfer. Estimation of the sequence divergence of homologous genes help us to understand divergence time, which makes it possible to understand the evolutionary patterns of speciation, gene duplication and gene transfer events. This protocol will provide a detailed bioinformatics pipeline on how to identify the homologous genes, compare their sequence divergence and phylogenetic relationships, focusing on homologous genes that show syntenic relationships using soybean (Glycine max) and common bean (Phaseolus vulgaris) as example species.

Evolutionary Divergence of TNL Disease-Resistant Proteins in Soybean (Glycine max) and Common Bean (Phaseolus vulgaris)

Disease-resistant genes (R genes) encode proteins that are involved in protecting plants from their pathogens and pests. Availability of complete genome sequences from soybean and common bean allowed us to perform a genome-wide identification and analysis of the Toll interleukin-1 receptor-like nucleotide-binding site leucine-rich repeat (TNL) proteins. Hidden Markov model (HMM) profiling of all protein sequences resulted in the identification of 117 and 77 regular TNL genes in soybean and common bean, respectively. We also identified TNL gene homologs with unique domains, and signal peptides as well as nuclear localization signals. The TNL genes in soybean formed 28 clusters located on 10 of the 20 chromosomes, with the majority found on chromosome 3, 6 and 16. Similarly, the TNL genes in common bean formed 14 clusters located on five of the 11 chromosomes, with the majority found on chromosome 10. Phylogenetic analyses of the TNL genes from Arabidopsis, soybean and common bean revealed less divergence within legumes relative to the divergence between legumes and Arabidopsis. Syntenic blocks were found between chromosomes Pv10 and Gm03, Pv07 and Gm10, as well as Pv01 and Gm14. The gene expression data revealed basal level expression and tissue specificity, while analysis of available microRNA data showed 37 predicted microRNA families involved in targeting the identified TNL genes in soybean and common bean.

Patterns and Consequences of Subgenome Differentiation Provide Insights into the Nature of Paleopolyploidy in Plants

Polyploidy is an important feature of plant genomes, but the nature of many polyploidization events remains to be elucidated. Here, we demonstrate that the evolutionary fates of the subgenomes in maize (Zea mays) and soybean (Glycine max) have followed different trajectories. One subgenome has been subject to relaxed selection, lower levels of gene expression, higher rates of transposable element accumulation, more small interfering RNAs and DNA methylation around genes, and higher rates of gene loss in maize, whereas none of these features were observed in soybean. Nevertheless, individual gene pairs exhibit differentiation with respect to these features in both species. In addition, we observed a higher number of chromosomal rearrangements and higher frequency of retention of duplicated genes in soybean than in maize. Furthermore, soybean "singletons" were found to be more frequently tandemly duplicated than "duplicates" in soybean, which may, to some extent, counteract the genome imbalance caused by gene loss. We propose that unlike in maize, in which two subgenomes were distinct prior to the allotetraploidization event and thus experienced global differences in selective constraints, in soybean, the two subgenomes were far less distinct prior to polyploidization, such that individual gene pairs, rather than subgenomes, experienced stochastic differences over longer periods of time, resulting in retention of the majority of duplicates.

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.

Current understanding of the pathways of flavonoid biosynthesis in model and crop plants

Flavonoids are a signature class of secondary metabolites formed from a relatively simple collection of scaffolds. They are extensively decorated by chemical reactions including glycosylation, methylation, and acylation. They are present in a wide variety of fruits and vegetables and as such in Western populations it is estimated that 20-50 mg of flavonoids are consumed daily per person. In planta they have demonstrated to contribute to both flower color and UV protection. Their consumption has been suggested to presenta wide range of health benefits. Recent technical advances allowing affordable whole genome sequencing, as well as a better inventory of species-by-species chemical diversity, have greatly advanced our understanding as to how flavonoid biosynthesis pathways vary across species. In parallel, reverse genetics combined with detailed molecular phenotyping is currently allowing us to elucidate the functional importance of individual genes and metabolites and by this means to provide further mechanistic insight into their biological roles. Here we provide an inventory of current knowledge of pathways of flavonoid biosynthesis in both the model plant Arabidopsis thaliana and a range of crop species, including tomato, maize, rice, and bean.

Genetic Architecture of Flooding Tolerance in the Dry Bean Middle-American Diversity Panel

Flooding is a devastating abiotic stress that endangers crop production in the twenty-first century. Because of the severe susceptibility of common bean (Phaseolus vulgaris L.) to flooding, an understanding of the genetic architecture and physiological responses of this crop will set the stage for further improvement. However, challenging phenotyping methods hinder a large-scale genetic study of flooding tolerance in common bean and other economically important crops. A greenhouse phenotyping protocol was developed to evaluate the flooding conditions at early stages. The Middle-American diversity panel (n = 272) of common bean was developed to capture most of the diversity exits in North American germplasm. This panel was evaluated for seven traits under both flooded and non-flooded conditions at two early developmental stages. A subset of contrasting genotypes was further evaluated in the field to assess the relationship between greenhouse and field data under flooding condition. A genome-wide association study using ~150 K SNPs was performed to discover genomic regions associated with multiple physiological responses. The results indicate a significant strong correlation (r > 0.77) between greenhouse and field data, highlighting the reliability of greenhouse phenotyping method. Black and small red beans were the least affected by excess water at germination stage. At the seedling stage, pinto and great northern genotypes were the most tolerant. Root weight reduction due to flooding was greatest in pink and small red cultivars. Flooding reduced the chlorophyll content to the greatest extent in the navy bean cultivars compared with other market classes. Races of Durango/Jalisco and Mesoamerica were separated by both genotypic and phenotypic data indicating the potential effect of eco-geographical variations. Furthermore, several loci were identified that potentially represent the antagonistic pleiotropy. The GWAS analysis revealed peaks at Pv08/1.6 Mb and Pv02/41 Mb that are associated with root weight and germination rate, respectively. These regions are syntenic with two QTL reported in soybean (Glycine max L.) that contribute to flooding tolerance, suggesting a conserved evolutionary pathway involved in flooding tolerance for these related legumes.

Molecular characterization of dihydroneopterin aldolase and aminodeoxychorismate synthase in common bean—genes coding for enzymes in the folate synthesis pathway

Common beans (Phaseolus vulgaris) are excellent sources of dietary folates, but different varieties contain different amounts of these compounds. Genes coding for dihydroneopterin aldolase (DHNA) and aminodeoxychorismate synthase (ADCS) of the folate synthesis pathway were characterized by PCR amplification, BAC clone sequencing, and whole genome sequencing. All DHNA and ADCS genes in the Mesoamerican cultivar OAC Rex were isolated and compared with those genes in the genome of Andean genotype G19833. Both genotypes have two functional DHNA genes and one pseudo gene. PvDHNA1 and PvDHNA2 proteins have similar secondary structures and conserved residues as DHNA homologs in Staphylococcus aureus and Arabidopsis. Sequence analysis and synteny mapping indicated that PvDHNA1 might be a duplicated and transposed copy of PvDHNA2. There is only one ADCS gene (PvADCS) identified in the bean genome and it is identical in OAC Rex and G19833. PvADCS has the conserved motifs required for catalytic activity similar to other plant ADCS homologs. DHNA and ADCS gene-specific markers were developed, mapped, and compared to their physical locations on chromosomes 1 and 7, respectively. The gene-specific markers developed in this study should be useful for detection and selection of varieties with enhanced folate contents in bean breeding programs.

Natural rhizobial diversity helps to reveal genes and QTLs associated with biological nitrogen fixation in common bean

Common bean is one of the most important crops for human feed, and the most important legume for direct consumption by millions of people, especially in developing countries. It is a promiscuous host legume in terms of nodulation, able to associate with a broad and diverse range of rhizobia, although the competitiveness for nodulation and the nitrogen fixation capacity of most of these strains is generally low. As a result, common bean is very inefficient for symbiotic nitrogen fixation, and nitrogen has to be supplied with chemical fertilizers. In the last years, symbiotic nitrogen fixation has received increasing attention as a sustainable alternative to nitrogen fertilizers, and also as a more economic and available one in poor countries. Therefore, optimization of nitrogen fixation of bean-rhizobia symbioses and selection of efficient rhizobial strains should be a priority, which begins with the study of the natural diversity of the symbioses and the rhizobial populations associated. Natural rhizobia biodiversity that nodulates common bean may be a source of adaptive alleles acting through phenotypic plasticity. Crosses between accessions differing for nitrogen fixation may combine alleles that never meet in nature. Another way to discover adaptive genes is to use association genetics to identify loci that common bean plants use for enhanced biological nitrogen fixation and, in consequence, for marker assisted selection for genetic improvement of symbiotic nitrogen fixation. In this review, rhizobial biodiversity resources will be discussed, together with what is known about the loci that underlie such genetic variation, and the potential candidate genes that may influence the symbiosis' fitness benefits, thus achieving an optimal nitrogen fixation capacity in order to help reduce reliance on nitrogen fertilizers in common bean.

Further evidence for the genetic basis of qualitative traits and their linkage relationships in Dolichos bean (Lablab purpureus L.)

An investigation on inheritance of qualitative traits in dolichos bean revealed biallelic monogenic control of photoperiodinduced sensitivity to flowering time and flower colour in F₂ and F₂ generations. While, growth habit and pod curvature are each controlled by two genes that exhibit classical complementary epistasis, raceme emergence was controlled by two genes that displayed classical inhibitory epistasis. The dominant alleles, at two different unlinked pairs of genes are necessary for plants to exhibit indeterminate growth habit and bear straight pods. Any other combination of alleles at the two pairs of genes result in plants displaying determinate growth habit and bearing curved pods.While, the genes controlling growth habit, PSFT and raceme emergence are linked. Those controlling flower colour and pod curvature are segregated independent of each other. These results are discussed in relation to strategies for breeding dolichos bean.

Comparative Transcriptome Analysis of Resistant and Susceptible Common Bean Genotypes in Response to Soybean Cyst Nematode Infection

Soybean cyst nematode (SCN; Heterodera glycines Ichinohe) reproduces on the roots of common bean (Phaseolus vulgaris L.) and can cause reductions in plant growth and seed yield. The molecular changes in common bean roots caused by SCN infection are unknown. Identification of genetic factors associated with SCN resistance could help in development of improved bean varieties with high SCN resistance. Gene expression profiling was conducted on common bean roots infected by SCN HG type 0 using next generation RNA sequencing technology. Two pinto bean genotypes, PI533561 and GTS-900, resistant and susceptible to SCN infection, respectively, were used as RNA sources eight days post inoculation. Total reads generated ranged between ~ 3.2 and 5.7 million per library and were mapped to the common bean reference genome. Approximately 70-90% of filtered RNA-seq reads uniquely mapped to the reference genome. In the inoculated roots of resistant genotype PI533561, a total of 353 genes were differentially expressed with 154 up-regulated genes and 199 down-regulated genes when compared to the transcriptome of non- inoculated roots. On the other hand, 990 genes were differentially expressed in SCN-inoculated roots of susceptible genotype GTS-900 with 406 up-regulated and 584 down-regulated genes when compared to non-inoculated roots. Genes encoding nucleotide-binding site leucine-rich repeat resistance (NLR) proteins, WRKY transcription factors, pathogenesis-related (PR) proteins and heat shock proteins involved in diverse biological processes were differentially expressed in both resistant and susceptible genotypes. Overall, suppression of the photosystem was observed in both the responses. Furthermore, RNA-seq results were validated through quantitative real time PCR. This is the first report describing genes/transcripts involved in SCN-common bean interaction and the results will have important implications for further characterization of SCN resistance genes in common bean.

Integration of Lupinus angustifolius L. (narrow-leafed lupin) genome maps and comparative mapping within legumes

Narrow-leafed lupin (Lupinus angustifolius L.) has recently been considered a reference genome for the Lupinus genus. In the present work, genetic and cytogenetic maps of L. angustifolius were supplemented with 30 new molecular markers representing lupin genome regions, harboring genes involved in nitrogen fixation during the symbiotic interaction of legumes and soil bacteria (Rhizobiaceae). Our studies resulted in the precise localization of bacterial artificial chromosomes (BACs) carrying sequence variants for early nodulin 40, nodulin 26, nodulin 45, aspartate aminotransferase P2, asparagine synthetase, cytosolic glutamine synthetase, and phosphoenolpyruvate carboxylase. Together with previously mapped chromosomes, the integrated L. angustifolius map encompasses 73 chromosome markers, including 5S ribosomal DNA (rDNA) and 45S rDNA, and anchors 20 L. angustifolius linkage groups to corresponding chromosomes. Chromosomal identification using BAC fluorescence in situ hybridization identified two BAC clones as narrow-leafed lupin centromere-specific markers, which served as templates for preliminary studies of centromere composition within the genus. Bioinformatic analysis of these two BACs revealed that centromeric/pericentromeric regions of narrow-leafed lupin chromosomes consisted of simple sequence repeats ordered into tandem repeats containing the trinucleotide and pentanucleotide simple sequence repeats AGG and GATAC, structured into long arrays. Moreover, cross-genus microsynteny analysis revealed syntenic patterns of 31 single-locus BAC clones among several legume species. The gene and chromosome level findings provide evidence of ancient duplication events that must have occurred very early in the divergence of papilionoid lineages. This work provides a strong foundation for future comparative mapping among legumes and may facilitate understanding of mechanisms involved in shaping legume chromosomes.

Comparative analysis of genome-wide Mlo gene family in Cajanus cajan and Phaseolus vulgaris

The Mlo gene was discovered in barley because the mutant 'mlo' allele conferred broad-spectrum, non-race-specific resistance to powdery mildew caused by Blumeria graminis f. sp. hordei. The Mlo genes also play important roles in growth and development of plants, and in responses to biotic and abiotic stresses. The Mlo gene family has been characterized in several crop species, but only a single legume species, soybean (Glycine max L.), has been investigated so far. The present report describes in silico identification of 18 CcMlo and 20 PvMlo genes in the important legume crops Cajanus cajan (L.) Millsp. and Phaseolus vulgaris L., respectively. In silico analysis of gene organization, protein properties and conserved domains revealed that the C. cajan and P. vulgaris Mlo gene paralogs are more divergent from each other than from their orthologous pairs. The comparative phylogenetic analysis classified CcMlo and PvMlo genes into three major clades. A comparative analysis of CcMlo and PvMlo proteins with the G. max Mlo proteins indicated close association of one CcMlo, one PvMlo with two GmMlo genes, indicating that there was no further expansion of the Mlo gene family after the separation of these species. Thus, most of the diploid species of eudicots might be expected to contain 15-20 Mlo genes. The genes CcMlo12 and 14, and PvMlo11 and 12 are predicted to participate in powdery mildew resistance. If this prediction were verified, these genes could be targeted by TILLING or CRISPR to isolate powdery mildew resistant mutants.

Construction of high resolution genetic linkage maps to improve the soybean genome sequence assembly Glyma1.01

BACKGROUND

A landmark in soybean research, Glyma1.01, the first whole genome sequence of variety Williams 82 (Glycine max L. Merr.) was completed in 2010 and is widely used. However, because the assembly was primarily built based on the linkage maps constructed with a limited number of markers and recombinant inbred lines (RILs), the assembled sequence, especially in some genomic regions with sparse numbers of anchoring markers, needs to be improved. Molecular markers are being used by researchers in the soybean community, however, with the updating of the Glyma1.01 build based on the high-resolution linkage maps resulting from this research, the genome positions of these markers need to be mapped.

RESULTS

Two high density genetic linkage maps were constructed based on 21,478 single nucleotide polymorphism loci mapped in the Williams 82 x G. soja (Sieb. & Zucc.) PI479752 population with 1083 RILs and 11,922 loci mapped in the Essex x Williams 82 population with 922 RILs. There were 37 regions or single markers where marker order in the two populations was in agreement but was not consistent with the physical position in the Glyma1.01 build. In addition, 28 previously unanchored scaffolds were positioned. Map data were used to identify false joins in the Glyma1.01 assembly and the corresponding scaffolds were broken and reassembled to the new assembly, Wm82.a2.v1. Based upon the plots of the genetic on physical distance of the loci, the euchromatic and heterochromatic regions along each chromosome in the new assembly were delimited. Genomic positions of the commonly used markers contained in BARCSOYSSR_1.0 database and the SoySNP50K BeadChip were updated based upon the Wm82.a2.v1 assembly.

CONCLUSIONS

The information will facilitate the study of recombination hot spots in the soybean genome, identification of genes or quantitative trait loci controlling yield, seed quality and resistance to biotic or abiotic stresses as well as other genetic or genomic research.

Divergence in Enzymatic Activities in the Soybean GST Supergene Family Provides New Insight into the Evolutionary Dynamics of Whole-Genome Duplicates

Whole-genome duplication (WGD), or polyploidy, is a major force in plant genome evolution. A duplicate of all genes is present in the genome immediately following a WGD event. However, the evolutionary mechanisms responsible for the loss of, or retention and subsequent functional divergence of polyploidy-derived duplicates remain largely unknown. In this study we reconstructed the evolutionary history of the glutathione S-transferase (GST) gene family from the soybean genome, and identified 72 GST duplicated gene pairs formed by a recent Glycine-specific WGD event occurring approximately 13 Ma. We found that 72% of duplicated GST gene pairs experienced gene losses or pseudogenization, whereas 28% of GST gene pairs have been retained in the soybean genome. The GST pseudogenes were under relaxed selective constraints, whereas functional GSTs were subject to strong purifying selection. Plant GST genes play important roles in stress tolerance and detoxification metabolism. By examining the gene expression responses to abiotic stresses and enzymatic properties of the ancestral and current proteins, we found that polyploidy-derived GST duplicates show the divergence in enzymatic activities. Through site-directed mutagenesis of ancestral proteins, this study revealed that nonsynonymous substitutions of key amino acid sites play an important role in the divergence of enzymatic functions of polyploidy-derived GST duplicates. These findings provide new insights into the evolutionary and functional dynamics of polyploidy-derived duplicate genes.

Proteny: discovering and visualizing statistically significant syntenic clusters at the proteome level

Background: With more and more genomes being sequenced, detecting synteny between genomes becomes more and more important. However, for microorganisms the genomic divergence quickly becomes large, resulting in different codon usage and shuffling of gene order and gene elements such as exons.

Results: We present Proteny, a methodology to detect synteny between diverged genomes. It operates on the amino acid sequence level to be insensitive to codon usage adaptations and clusters groups of exons disregarding order to handle diversity in genomic ordering between genomes. Furthermore, Proteny assigns significance levels to the syntenic clusters such that they can be selected on statistical grounds. Finally, Proteny provides novel ways to visualize results at different scales, facilitating the exploration and interpretation of syntenic regions. We test the performance of Proteny on a standard ground truth dataset, and we illustrate the use of Proteny on two closely related genomes (two different strains of Aspergillus niger) and on two distant genomes (two species of Basidiomycota). In comparison to other tools, we find that Proteny finds clusters with more true homologies in fewer clusters that contain more genes, i.e. Proteny is able to identify a more consistent synteny. Further, we show how genome rearrangements, assembly errors, gene duplications and the conservation of specific genes can be easily studied with Proteny.

Availability and implementation: Proteny is freely available at the Delft Bioinformatics Lab website http://bioinformatics.tudelft.nl/dbl/software.

Contact:t.gehrmann@tudelft.nl

Supplementary information: Supplementary data are available at Bioinformatics online.

Genome-wide identification of evolutionarily conserved alternative splicing events in flowering plants

Alternative splicing (AS) plays important roles in many plant functions, but its conservation across the plant kingdom is not known. We describe a methodology to identify AS events and identify conserved AS events across large phylogenetic distances using RNA-Seq datasets. We applied this methodology to transcriptome data from nine angiosperms including Amborella, the single sister species to all other extant flowering plants. AS events within 40–70% of the expressed multi-exonic genes per species were found, 27,120 of which are conserved among two or more of the taxa studied. While many events are species specific, many others are shared across long evolutionary distances suggesting they have functional significance. Conservation of AS event data provides an estimate of the number of ancestral AS events present at each node of the tree representing the nine species studied. Furthermore, the presence or absence of AS isoforms between species with different whole genome duplication (WGD) histories provides the opportunity to examine the impact of WDG on AS potential. Examining AS in gene families identifies those with high rates of AS, and conservation can distinguish ancient events vs. recent or species specific adaptations. The MADS-box and SR protein families are found to represent families with low and high occurrences of AS, respectively, yet their AS events were likely present in the MRCA of angiosperms.

Evolutionary patterns and coevolutionary consequences of MIRNA genes and microRNA targets triggered by multiple mechanisms of genomic duplications in soybean

The evolutionary dynamics of duplicated protein-encoding genes (PEGs) is well documented. However, the evolutionary patterns and consequences of duplicated MIRNAs and the potential influence on the evolution of their PEG targets are poorly understood. Here, we demonstrate the evolution of plant MIRNAs subsequent to a recent whole-genome duplication. Overall, the retention of MIRNA duplicates was correlated to the retention of adjacent PEG duplicates, and the retained MIRNA duplicates exhibited a higher level of interspecific preservation of orthologs than singletons, suggesting that the retention of MIRNA duplicates is related to their functional constraints and local genomic stability. Nevertheless, duplication status, rather than local genic collinearity, was the primary determinant of levels of nucleotide divergence of MIRNAs. In addition, the retention of duplicated MIRNAs appears to be associated with the retention of their corresponding duplicated PEG targets. Furthermore, we characterized the evolutionary novelty of a legume-specific microRNA (miRNA) family, which resulted from rounds of genomic duplication, and consequent dynamic evolution of its NB-LRR targets, an important gene family with primary roles in plant-pathogen interactions. Together, these observations depict evolutionary patterns and novelty of MIRNAs in the context of genomic duplication and evolutionary interplay between MIRNAs and their PEG targets mediated by miRNAs.

Intra- and interchromosomal rearrangements between cowpea [Vigna unguiculata (L.) Walp.] and common bean (Phaseolus vulgaris L.) revealed by BAC-FISH

Cowpea (Vigna unguiculata) is an annual legume grown in tropical and subtropical regions, which is economically relevant due to high protein content in dried beans, green pods, and leaves. In this work, a comparative cytogenetic study between V. unguiculata and Phaseolus vulgaris (common bean) was conducted using BAC-FISH. Sequences previously mapped in P. vulgaris chromosomes (Pv) were used as probes in V. unguiculata chromosomes (Vu), contributing to the analysis of macrosynteny between both legumes. Thirty-seven clones from P. vulgaris 'BAT93' BAC library, corresponding to its 11 linkage groups, were hybridized in situ. Several chromosomal rearrangements were identified, such as translocations (between BACs from Pv1 and Pv8; Pv2 and Pv3; as well as Pv2 and Pv11), duplications (BAC from Pv3), as well as paracentric and pericentric inversions (BACs from Pv3, and Pv4, respectively). Two BACs (from Pv2 and Pv7), which hybridized at terminal regions in almost all P. vulgaris chromosomes, showed single-copy signal in Vu. Additionally, 17 BACs showed no signal in V. unguiculata chromosomes. The present results demonstrate the feasibility of using BAC libraries in comparative chromosomal mapping and karyotype evolution studies between Phaseolus and Vigna species, and revealed several macrosynteny and collinearity breaks among both legumes.

Molecular marker applications in plants

Individuals within a population of a sexually reproducing species will have some degree of heritable genomic variation caused by mutations, insertion/deletions (INDELS), inversions, duplications, and translocations. Such variation can be detected and screened using molecular, or genetic, markers. By definition, molecular markers are genetic loci that can be easily tracked and quantified in a population and may be associated with a particular gene or trait of interest. This chapter will review the current major applications of molecular markers in plants.

Genetic control of flowering time in legumes

The timing of flowering, and in particular the degree to which it is responsive to the environment, is a key factor in the adaptation of a given species to various eco-geographic locations and agricultural practices. Flowering time variation has been documented in many crop legumes, and selection for specific variants has permitted significant expansion and improvement in cultivation, from prehistoric times to the present day. Recent advances in legume genomics have accelerated the process of gene identification and functional analysis, and opened up new prospects for a molecular understanding of flowering time adaptation in this important crop group. Within the legumes, two species have been prominent in flowering time studies; the vernalization-responsive long-day species pea (Pisum sativum) and the warm-season short-day plant soybean (Glycine max). Analysis of flowering in these species is now being complemented by reverse genetics capabilities in the model legumes Medicago truncatula and Lotus japonicus, and the emergence of genome-scale resources in a range of other legumes. This review will outline the insights gained from detailed forward genetic analysis of flowering time in pea and soybean, highlighting the importance of light perception, the circadian clock and the FT family of flowering integrators. It discusses the current state of knowledge on genetic mechanisms for photoperiod and vernalization response, and concludes with a broader discussion of flowering time adaptation across legumes generally.

Remnants of the Legume Ancestral Genome Preserved in Gene-Rich Regions: Insights from Lupinus angustifolius Physical, Genetic, and Comparative Mapping

The narrow-leafed lupin (Lupinus angustifolius) was recently considered as a legume reference species. Genetic resources have been developed, including a draft genome sequence, linkage maps, nuclear DNA libraries, and cytogenetic chromosome-specific landmarks. Here, we used a complex approach, involving DNA fingerprinting, sequencing, genetic mapping, and molecular cytogenetics, to localize and analyze L. angustifolius gene-rich regions (GRRs). A L. angustifolius genomic bacterial artificial chromosome (BAC) library was screened with short sequence repeat (SSR)-based probes. Selected BACs were fingerprinted and assembled into contigs. BAC-end sequence (BES) annotation allowed us to choose clones for sequencing, targeting GRRs. Additionally, BESs were aligned to the scaffolds of the genome sequence. The genetic map was supplemented with 35 BES-derived markers, distributed in 14 linkage groups and tagging 37 scaffolds. The identified GRRs had an average gene density of 19.6 genes/100 kb and physical-to-genetic distance ratios of 11 to 109 kb/cM. Physical and genetic mapping was supported by multi-BAC-fluorescence in situ hybridization (FISH), and five new linkage groups were assigned to the chromosomes. Syntenic links to the genome sequences of five legume species (Medicago truncatula, Glycine max, Lotus japonicus, Phaseolus vulgaris, and Cajanus cajan) were identified. The comparative mapping of the two largest lupin GRRs provides novel evidence for ancient duplications in all of the studied species. These regions are conserved among representatives of the main clades of Papilionoideae. Furthermore, despite the complex evolution of legumes, some segments of the nuclear genome were not substantially modified and retained their quasi-ancestral structures. Cytogenetic markers anchored in these regions constitute a platform for heterologous mapping of legume genomes.

Development and validation of DNA markers linked to Sdvy-1, a common bean gene conferring resistance to the yellowing strain of Soybean dwarf virus

The yellowing strain of Soybean dwarf virus (SbDV-YS) causes yellowing and yield loss in common bean (Phaseolus vulgaris). The most effective control is achieved through breeding for resistance. An indeterminate climbing cultivar with a white seed coat, 'Oofuku', is resistant to SbDV-YS in inoculation tests. We crossed 'Oofuku' with an elite cultivar, 'Taisho-Kintoki', which is SbDV-YS-susceptible, determinate dwarf with a red-purple seed coat, and performed amplified-fragment-length polymorphism analysis of F3 lines. From nucleotide sequences of the resistant-specific fragments and their flanking regions, we developed five DNA markers, of which DV86, DV386, and DV398 were closely linked to Sdvy-1, a resistance gene. Using the markers, we developed 'Toiku-B79' and 'Toiku-B80', the near-isogenic lines (NILs) incorporating Sdvy-1 in the background of 'Taisho-Kintoki'. The NILs had similar growth habit, maturity date and seed shape to those of 'Taisho-Kintoki'. The quality of boiled beans was also similar, except that the NILs had more seed coat cracking than 'Taisho-Kintoki'. The NILs showed no SbDV-YS infection in inoculation tests. We suggest that Sdvy-1 is a useful source of SbDV-YS resistance in common bean.

Analysis of BAC-end sequences in common bean (Phaseolus vulgaris L.) towards the development and characterization of long motifs SSRs

The increasing volume of genomic data on the Phaseolus vulgaris species have contributed to its importance as a model genetic species and positively affected the investigation of other legumes of scientific and economic value. To expand and gain a more in-depth knowledge of the common bean genome, the ends of a number of bacterial artificial chromosome (BAC) were sequenced, annotated and the presence of repetitive sequences was determined. In total, 52,270 BESs (BAC-end sequences), equivalent to 32 Mbp (~6 %) of the genome, were processed. In total, 3,789 BES-SSRs were identified, with a distribution of one SSR (simple sequence repeat) per 8.36 kbp and 2,000 were suitable for the development of SSRs, of which 194 were evaluated in low-resolution screening. From 40 BES-SSRs based on long motifs SSRs (≥ trinucleotides) analyzed in high-resolution genotyping, 34 showed an equally good amplification for the Andean and for the Mesoamerican genepools, exhibiting an average gene diversity (H E) of 0.490 and 5.59 alleles/locus, of which six classified as Class I showed a H E ≥ 0.7. The PCoA and structure analysis allowed to discriminate the gene pools (K = 2, FST = 0.733). From the 52,270 BESs, 2 % corresponded to transcription factors and 3 % to transposable elements. Putative functions for 24,321 BESs were identified and for 19,363 were assigned functional categories (gene ontology). This study identified highly polymorphic BES-SSRs containing tri- to hexanucleotides motifs and bringing together relevant genetic characteristics useful for breeding programs. Additionally, the BESs were incorporated into the international genome-sequencing project for the common bean.

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

BACKGROUND

Common bean (Phaseolus vulgaris) is grown throughout the world and comprises roughly 50% of the grain legumes consumed worldwide. Despite this, genetic resources for common beans have been lacking. Next generation sequencing, has facilitated our investigation of the gene expression profiles associated with biologically important traits in common bean. An increased understanding of gene expression in common bean will improve our understanding of gene expression patterns in other legume species.

RESULTS

Combining recently developed genomic resources for Phaseolus vulgaris, including predicted gene calls, with RNA-Seq technology, we measured the gene expression patterns from 24 samples collected from seven tissues at developmentally important stages and from three nitrogen treatments. Gene expression patterns throughout the plant were analyzed to better understand changes due to nodulation, seed development, and nitrogen utilization. We have identified 11,010 genes differentially expressed with a fold change ≥ 2 and a P-value < 0.05 between different tissues at the same time point, 15,752 genes differentially expressed within a tissue due to changes in development, and 2,315 genes expressed only in a single tissue. These analyses identified 2,970 genes with expression patterns that appear to be directly dependent on the source of available nitrogen. Finally, we have assembled this data in a publicly available database, The Phaseolus vulgaris Gene Expression Atlas (Pv GEA), http://plantgrn.noble.org/PvGEA/ . Using the website, researchers can query gene expression profiles of their gene of interest, search for genes expressed in different tissues, or download the dataset in a tabular form.

CONCLUSIONS

These data provide the basis for a gene expression atlas, which will facilitate functional genomic studies in common bean. Analysis of this dataset has identified genes important in regulating seed composition and has increased our understanding of nodulation and impact of the nitrogen source on assimilation and distribution throughout the plant.

In silico analysis of high affinity potassium transporter (HKT) isoforms in different plants

BACKGROUND

High affinity potassium transporters (HKTs) are located in the plasma membrane of the vessels and have significant influence on salt tolerance in some plants. They exclude Na(+) from the parenchyma cells to reduce Na(+) concentration. Despite many studies, the underlying regulatory mechanisms and the exact functions of HKTs within different genomic backgrounds are relatively unknown. In this study, various bioinformatics techniques, including promoter analysis, identification of HKT-surrounding genes, and construction of gene networks, were applied to investigate the HKT regulatory mechanism.

RESULTS

Promoter analysis showed that rice HKTs carry ABA response elements. Additionally, jasmonic acid response elements were detected on promoter region of TmHKT1;5. In silico synteny highlighted several unknown and new loci near rice, Arabidopsis thaliana and Physcomitrella patent HKTs, which may play a significant role in salt stress tolerance in concert with HKTs. Gene network prediction unravelled that crosstalk between jasmonate and ethylene reduces AtHKT1;1 expression. Furthermore, antiporter and transferase proteins were found in AtHKT1;1 gene network. Interestingly, regulatory elements on the promoter region of HKT in wild genotype (TmHKT1;5) were more frequent and variable than the ones in cultivated wheat (TaHKT1;5) which provides the possibility of rapid response and better understanding of environmental conditions for wild genotype.

CONCLUSION

Detecting ABA and jasmonic acid response elements on promoter regions of HKTs provide valuable clues on underlying regulatory mechanisms of HKTs. In silico synteny and pathway discovery indicated several candidates which act in concert with HKTs in stress condition. We highlighted different arrangement of regulatory elements on promoter region of wild wheat (TmHKT1;5) compared to bread wheat (TaHKT1;5) in this study.

Fine mapping of Co-x, an anthracnose resistance gene to a highly virulent strain of Colletotrichum lindemuthianum in common bean

The Co - x anthracnose R gene of common bean was fine-mapped into a 58 kb region at one end of chromosome 1, where no canonical NB-LRR-encoding genes are present in G19833 genome sequence. Anthracnose, caused by the phytopathogenic fungus Colletotrichum lindemuthianum, is one of the most damaging diseases of common bean, Phaseolus vulgaris. Various resistance (R) genes, named Co-, conferring race-specific resistance to different strains of C. lindemuthianum have been identified. The Andean cultivar JaloEEP558 was reported to carry Co-x on chromosome 1, conferring resistance to the highly virulent strain 100. To fine map Co-x, 181 recombinant inbred lines derived from the cross between JaloEEP558 and BAT93 were genotyped with polymerase chain reaction (PCR)-based markers developed using the genome sequence of the Andean genotype G19833. Analysis of RILs carrying key recombination events positioned Co-x at one end of chromosome 1 to a 58 kb region of the G19833 genome sequence. Annotation of this target region revealed eight genes: three phosphoinositide-specific phospholipases C (PI-PLC), one zinc finger protein and four kinases, suggesting that Co-x is not a classical nucleotide-binding leucine-rich encoding gene. In addition, we identified and characterized the seven members of common bean PI-PLC gene family distributed into two clusters located at the ends of chromosomes 1 and 8. Co-x is not a member of Co-1 allelic series since these two genes are separated by at least 190 kb. Comparative analysis between soybean and common bean revealed that the Co-x syntenic region, located at one end of Glycine max chromosome 18, carries Rhg1, a major QTL contributing to soybean cyst nematode resistance. The PCR-based markers generated in this study should be useful in marker-assisted selection for pyramiding Co-x with other R genes.

A reference genome for common bean and genome-wide analysis of dual domestications

Common bean (Phaseolus vulgaris L.) is the most important grain legume for human consumption and has a role in sustainable agriculture owing to its ability to fix atmospheric nitrogen. We assembled 473 Mb of the 587-Mb genome and genetically anchored 98% of this sequence in 11 chromosome-scale pseudomolecules. We compared the genome for the common bean against the soybean genome to find changes in soybean resulting from polyploidy. Using resequencing of 60 wild individuals and 100 landraces from the genetically differentiated Mesoamerican and Andean gene pools, we confirmed 2 independent domestications from genetic pools that diverged before human colonization. Less than 10% of the 74 Mb of sequence putatively involved in domestication was shared by the two domestication events. We identified a set of genes linked with increased leaf and seed size and combined these results with quantitative trait locus data from Mesoamerican cultivars. Genes affected by domestication may be useful for genomics-enabled crop improvement.

Comparative mapping of the wild perennial Glycine latifolia and soybean (G. max) reveals extensive chromosome rearrangements in the genus Glycine

Soybean (Glycine max L. Mer.), like many cultivated crops, has a relatively narrow genetic base and lacks diversity for some economically important traits. Glycine latifolia (Benth.) Newell & Hymowitz, one of the 26 perennial wild Glycine species related to soybean in the subgenus Glycine Willd., shows high levels of resistance to multiple soybean pathogens and pests including Alfalfa mosaic virus, Heterodera glycines Ichinohe and Sclerotinia sclerotiorum (Lib.) de Bary. However, limited information is available on the genomes of these perennial Glycine species. To generate molecular resources for gene mapping and identification, high-density linkage maps were constructed for G. latifolia using single nucleotide polymorphism (SNP) markers generated by genotyping by sequencing and evaluated in an F2 population and confirmed in an F5 population. In each population, greater than 2,300 SNP markers were selected for analysis and segregated to form 20 large linkage groups. Marker orders were similar in the F2 and F5 populations. The relationships between G. latifolia linkage groups and G. max and common bean (Phaseolus vulgaris L.) chromosomes were examined by aligning SNP-containing sequences from G. latifolia to the genome sequences of G. max and P. vulgaris. Twelve of the 20 G. latifolia linkage groups were nearly collinear with G. max chromosomes. The remaining eight G. latifolia linkage groups appeared to be products of multiple interchromosomal translocations relative to G. max. Large syntenic blocks also were observed between G. latifolia and P. vulgaris. These experiments are the first to compare genome organizations among annual and perennial Glycine species and common bean. The development of molecular resources for species closely related to G. max provides information into the evolution of genomes within the genus Glycine and tools to identify genes within perennial wild relatives of cultivated soybean that could be beneficial to soybean production.

High density SNP and SSR-based genetic maps of two independent oil palm hybrids

BACKGROUND

Oil palm is an important perennial oil crop with an extremely long selection cycle of 10 to 12 years. As such, any tool that speeds up its genetic improvement process, such as marker-assisted breeding is invaluable. Previously, genetic linkage maps based on AFLP, RFLP and SSR markers were developed and QTLs for fatty acid composition and yield components identified. High density genetic maps of crosses of different genetic backgrounds are indispensable tools for investigating oil palm genetics. They are also useful for comparative mapping analyses to identify markers closely linked to traits of interest.

RESULTS

A 4.5 K customized oil palm SNP array was developed using the Illumina Infinium platform. The SNPs and 252 SSRs were genotyped on two mapping populations, an intraspecific cross with 87 palms and an interspecific cross with 108 palms. Parental maps with 16 linkage groups (LGs), were constructed for the three fruit forms of E. guineensis (dura, pisifera and tenera). Map resolution was further increased by integrating the dura and pisifera maps into an intraspecific integrated map with 1,331 markers spanning 1,867 cM. We also report the first map of a Colombian E. oleifera, comprising 10 LGs with 65 markers spanning 471 cM. Although not very dense due to the high level of homozygosity in E. oleifera, the LGs were successfully integrated with the LGs of the tenera map. Direct comparison between the parental maps identified 603 transferable markers polymorphic in at least two of the parents. Further analysis revealed a high degree of marker transferability covering 1,075 cM, between the intra- and interspecific integrated maps. The interspecific cross displayed higher segregation distortion than the intraspecific cross. However, inclusion of distorted markers in the genetic maps did not disrupt the marker order and no map expansion was observed.

CONCLUSIONS

The high density SNP and SSR-based genetic maps reported in this paper have greatly improved marker density and genome coverage in comparison with the first reference map based on AFLP and SSR markers. Therefore, it is foreseen that they will be more useful for fine mapping of QTLs and whole genome association mapping studies in oil palm.

Systems and evolutionary characterization of microRNAs and their underlying regulatory networks in soybean cotyledons

MicroRNAs (miRNAs) are an emerging class of small RNAs regulating a wide range of biological processes. Soybean cotyledons evolved as sink tissues to synthesize and store seed reserves which directly affect soybean seed yield and quality. However, little is known about miRNAs and their regulatory networks in soybean cotyledons. We sequenced 292 million small RNA reads expressed in soybean cotyledons, and discovered 130 novel miRNA genes and 72 novel miRNA families. The cotyledon miRNAs arose at various stages of land plant evolution. Evolutionary analysis of the miRNA genes in duplicated genome segments from the recent Glycine whole genome duplication revealed that the majority of novel soybean cotyledon miRNAs were young, and likely arose after the duplication event 13 million years ago. We revealed the evolutionary pathway of a soybean cotyledon miRNA family (soy-miR15/49) that evolved from a neutral invertase gene through an inverted duplication and a series of DNA amplification and deletion events. A total of 304 miRNA genes were expressed in soybean cotyledons. The miRNAs were predicted to target 1910 genes, and form complex miRNA networks regulating a wide range of biological pathways in cotyledons. The comprehensive characterization of the miRNAs and their underlying regulatory networks at gene, pathway and system levels provides a foundation for further studies of miRNAs in cotyledons.

Annotation and sequence diversity of transposable elements in common bean (Phaseolus vulgaris)

Common bean (Phaseolus vulgaris) is an important legume crop grown and consumed worldwide. With the availability of the common bean genome sequence, the next challenge is to annotate the genome and characterize functional DNA elements. Transposable elements (TEs) are the most abundant component of plant genomes and can dramatically affect genome evolution and genetic variation. Thus, it is pivotal to identify TEs in the common bean genome. In this study, we performed a genome-wide transposon annotation in common bean using a combination of homology and sequence structure-based methods. We developed a 2.12-Mb transposon database which includes 791 representative transposon sequences and is available upon request or from www.phytozome.org. Of note, nearly all transposons in the database are previously unrecognized TEs. More than 5,000 transposon-related expressed sequence tags (ESTs) were detected which indicates that some transposons may be transcriptionally active. Two Ty1-copia retrotransposon families were found to encode the envelope-like protein which has rarely been identified in plant genomes. Also, we identified an extra open reading frame (ORF) termed ORF2 from 15 Ty3-gypsy families that was located between the ORF encoding the retrotransposase and the 3'LTR. The ORF2 was in opposite transcriptional orientation to retrotransposase. Sequence homology searches and phylogenetic analysis suggested that the ORF2 may have an ancient origin, but its function is not clear. These transposon data provide a useful resource for understanding the genome organization and evolution and may be used to identify active TEs for developing transposon-tagging system in common bean and other related genomes.