Global transcriptome analysis reveals resistance genes in the early response of common bean (Phaseolus vulgaris L.) to Colletotrichum lindemuthianum.

BACKGROUND: Disease can drastically impair common bean (Phaseolus vulgaris L.) production. Anthracnose, caused by the fungal pathogen Colletotrichum lindemuthianum (Sacc. and Magnus) Briosi and Cavara, is one of the diseases that are widespread and cause serious economic loss in common bean. RESULTS: Transcriptome analysis of the early response of common bean to anthracnose was performed using two resistant genotypes, Hongyundou and Honghuayundou, and one susceptible genotype, Jingdou. A total of 9,825 differentially expressed genes (DEGs) responding to pathogen infection and anthracnose resistance were identified by differential expression analysis. By using weighted gene coexpression network analysis (WGCNA), 2,051 DEGs were found to be associated with two resistance-related modules. Among them, 463 DEGs related to anthracnose resistance were considered resistance-related candidate genes. Nineteen candidate genes were coexpressed with three resistance genes, Phvul.001G243600, Phvul.001G243700 and Phvul.001G243800. To further identify resistance genes, 46 candidate genes were selected for experimental validation using salicylic acid (SA) and methyl jasmonate (MeJA). The results indicated that 38 candidate genes that responded to SA/MeJA treatment may be involved in anthracnose resistance in common bean. CONCLUSIONS: This study identified 38 resistance-related candidate genes involved in the early response of common bean, and 19 resistance-related candidate genes were coexpressed with anthracnose resistance genes. This study identified putative resistance genes for further resistance genetic investigation and provides an important reference for anthracnose resistance breeding in common bean.

Genome-wise association study identified genomic regions associated with drought tolerance in mungbean (Vigna radiata (L.) R. Wilczek).

KEY MESSAGE: A total of 282 mungbean accessions were resequenced to identify genome-wide variants and construct a highly precise variant map, and drought tolerance-related loci and superior alleles were identified by GWAS. Mungbean (Vigna radiata (L.) R. Wilczek) is an important food legume crop that is highly adapted to drought environments, but severe drought significantly curtails mungbean production. Here, we resequenced 282 mungbean accessions to identify genome-wide variants and constructed a highly precise map of mungbean variants. A genome-wide association study was performed to identify genomic regions for 14 drought tolerance-related traits in plants grown under stress and well-watered conditions over three years. One hundred forty-six SNPs associated with drought tolerance were detected, and 26 candidate loci associated with more than two traits were subsequently selected. Two hundred fifteen candidate genes were identified at these loci, including eleven transcription factor genes, seven protein kinase genes and other protein coding genes that may respond to drought stress. Furthermore, we identified superior alleles that were associated with drought tolerance and positively selected during the breeding process. These results provide valuable genomic resources for molecular breeding and will accelerate future efforts aimed at mungbean improvement.

Fast integration and accumulation of beneficial breeding alleles through an AB-NAMIC strategy in wheat.

Wheat (Triticum aestivum) is among the most important staple crops for safeguarding the food security of the growing world population. To bridge the gap between genebank diversity and breeding programs, we developed an advanced backcross-nested association mapping plus inter-crossed population (AB-NAMIC) by crossing three popular wheat cultivars as recurrent founders to 20 germplasm lines from a mini core collection. Selective backcrossing combined with selection against undesirable traits and extensive crossing within and between sub-populations created new opportunities to detect unknown genes and increase the frequency of beneficial alleles in the AB-NAMIC population. We performed phenotyping of 590 AB-NAMIC lines and a natural panel of 476 cultivars for six consecutive growing seasons and genotyped these 1066 lines with a 660K SNP array. Genome-wide association studies of both panels for plant development and yield traits demonstrated improved power to detect rare alleles and loci with medium genetic effects in AB-NAMIC. Notably, genome-wide association studies in AB-NAMIC detected the candidate gene TaSWEET6-7B (TraesCS7B03G1216700), which has high homology to the rice SWEET6b gene and exerts strong effects on adaptation and yield traits. The commercial release of two derived AB-NAMIC lines attests to its direct applicability in wheat improvement. Valuable information on genome-wide association study mapping, candidate genes, and their haplotypes for breeding traits are available through WheatGAB. Our research provides an excellent framework for fast-tracking exploration and accumulation of beneficial alleles stored in genebanks.

Development of a Model for Genomic Prediction of Multiple Traits in Common Bean Germplasm, Based on Population Structure.

Due to insufficient identification and in-depth investigation of existing common bean germplasm resources, it is difficult for breeders to utilize these valuable genetic resources. This situation limits the breeding and industrial development of the common bean (Phaseolus vulgaris L.) in China. Genomic prediction (GP) is a breeding method that uses whole-genome molecular markers to calculate the genomic estimated breeding value (GEBV) of candidate materials and select breeding materials. This study aimed to use genomic prediction to evaluate 15 traits in a collection of 628 common bean lines (including 484 landraces and 144 breeding lines) to determine a common bean GP model. The GP model constructed by landraces showed a moderate to high predictive ability (ranging from 0.59-0.88). Using all landraces as a training set, the predictive ability of the GP model for most traits was higher than that using the landraces from each of two subgene pools, respectively. Randomly selecting breeding lines as additional training sets together with landrace training sets to predict the remaining breeding lines resulted in a higher predictive ability based on principal components analysis. This study constructed a widely applicable GP model of the common bean based on the population structure, and encouraged the development of GP models to quickly aggregate excellent traits and accelerate utilization of germplasm resources.

QTL mapping and identification of genes associated with the resistance to Acanthoscelides obtectus in cultivated common bean using a high-density genetic linkage map.

BACKGROUND: Common bean (Phaseolus vulgaris L.) is an important agricultural product with large nutritional value, and the insect pest Acanthoscelides obtectus (Say) seriously affects its product quality and commodity quality during storage. Few researches on genes of bruchid resistance have investigated in common bean cultivars. RESULTS: In this study, a bruchid-resistant cultivar black kidney bean and a highly susceptible accession Longyundou3 from different gene banks were crossed to construct a recombinant inbred line population. The genetic analysis indicated a quantitative inheritance of the bruchid resistance trait controlled by polygenes. A high-density genetic map of a total map distance of 1283.68 cM with an average interval of 0.61 cM between each marker was constructed using an F6 population of 157 recombinant inbred lines. The map has 3106 bin markers, containing 2,234,769 SNPs. Using the high-density genetic map, a new quantitative trait locus for the resistance to Acanthoscelides obtectus was identified on chromosome 6. New molecular markers based on the candidate region were developed, and this locus was further delimited to an interval of 122.3 kb between SSR markers I6-4 and I6-16 using an F2 population. This region comprised five genes. Phvul.006G003700, which encodes a bifunctional inhibitor, may be a potential candidate gene for bruchid resistance. Sequencing analysis of candidate gene identified a 5 bp insertion-deletion in promoter of gene Phvul.006G003700 between two parents. Expression analysis of candidate gene revealed that the expression level of Phvul.006G003700 in bruchid-resistant parent was markedly higher than that in bruchid-susceptible parent both in dry seeds and leaves. CONCLUSIONS: A high-density genetic linkage map was constructed utilizing whole-genome resequencing and one new QTL for bruchid resistance was identified on chromosome 6 in common bean cultivar. Phvul.006G003700 (encoding a bifunctional inhibitor) may be a potential candidate gene. These results may form the basis for further research to reveal the bruchid resistance molecular mechanism of common bean.

The aquaporin gene PvXIP1;2 conferring drought resistance identified by GWAS at seedling stage in common bean.

KEY MESSAGE: A whole-genome resequencing-derived SNP dataset used for genome-wide association analysis revealed 12 loci significantly associated with drought stress based on survival rate after drought stress at seedling stage. We further confirmed the drought-related function of an aquaporin gene (PvXIP1;2) located at Locus_10. A variety of adverse conditions, including drought stress, severely affect common bean production. Molecular breeding for drought resistance has been proposed as an effective and practical way to improve the drought resistance of common bean. A genome-wide association analysis was conducted to identify drought-related loci based on survival rates at the seedling stage using a natural population consisting of 400 common bean accessions and 3,832,340 SNPs. The coefficient of variation ranged from 40.90 to 56.22% for survival rates in three independent experiments. A total of 12 associated loci containing 89 significant SNPs were identified for survival rates at the seedling stage. Four loci overlapped in the region of the QTLs reported to be associated with drought resistance. According to the expression profiles, gene annotations and references of the functions of homologous genes in Arabidopsis, 39 genes were considered potential candidate genes selected from 199 genes annotated within all associated loci. A stable locus (Locus_10) was identified on chromosome 11, which contained LEA, aquaporin, and proline-rich protein genes. We further confirmed the drought-related function of an aquaporin (PvXIP1;2) located at Locus_10 by expression pattern analysis, phenotypic analysis of PvXIP1;2-overexpressing Arabidopsis and Agrobacterium rhizogenes-mediated hairy root transformation systems, indicating that the association results can facilitate the efficient identification of genes related to drought resistance. These loci and their candidate genes provide a foundation for crop improvement via breeding for drought resistance in common bean.

Transcriptome Analysis of Resistance to Fusarium Wilt in Mung Bean (Vigna radiata L.).

Fusarium wilt is a destructive soil-borne disease that threatens the production of mung bean. Mung bean lines Zheng8-4 and Zheng8-20 show high resistance and high susceptibility to Fusarium wilt, respectively. Transcriptome analysis was carried out to identify candidate genes involved in Fusarium wilt resistance using Zheng8-4 and Zheng8-20 at 0, 0.5, 1, 2, and 4 days post inoculation (dpi). Differential expression analysis showed that 3,254 genes responded to pathogen infection and were differentially expressed in the resistant and susceptible lines. Weighted gene co-expression network analysis (WGCNA) was also performed to identify five modules highly correlated with Fusarium wilt resistance, in which 453 differentially expressed genes (DEGs) were considered likely to be involved in Fusarium wilt resistance. Among these DEGs, we found 24 genes encoding resistance (R) proteins, 22 encoding protein kinases, 20 belonging to transcription factor families, 34 encoding proteins with oxidoreductase activity, 17 involved in stimulation/stress responses, and 54 annotated to pathogen resistance-related pathways. Finally, 27 annotated genes were further selected as candidate genes of Fusarium wilt resistance in mung bean. This study identifies novel potential resistance-related genes against Fusarium wilt and provides a theoretical basis for further investigation of Fusarium wilt resistance in mung bean breeding.

Genetic dissection of drought resistance based on root traits at the bud stage in common bean.

KEY MESSAGE: A whole-genome resequencing-derived SNP dataset used for genome-wide association analysis revealed 196 loci significantly associated with drought stress based on root traits. Candidate genes identified in the regions of these loci include homologs of known drought resistance genes in A. thaliana. Drought is the main abiotic constraint of the production of common bean. Improved adaptation to drought environments has become a main goal of crop breeding due to the increasing scarcity of water that will occur in the future. The overall objective of our study was to identify genomic regions associated with drought resistance based on root traits using genome-wide association analysis. A natural population of 438 common bean accessions was evaluated for root traits: root surface area, root average diameter, root volume, total root length, taproot length, lateral root number, root dry weight, lateral root length, special root weight/length, using seed germination pouches under drought conditions and in well-watered environments. The coefficient of variation ranged from 11.24% (root average diameter) to 38.19% (root dry weight) in the well-watered environment and from 9.61% (root average diameter) to 39.05% (lateral root length) under drought stress. A whole-genome resequencing-derived SNP dataset revealed 196 loci containing 230 candidate SNPs associated with drought resistance. Seventeen candidate SNPs were simultaneously associated with more than two traits. Forty-one loci were simultaneously associated with more than two traits, and eleven loci were colocated with loci previously reported to be related to drought resistance. Candidate genes of the associated loci included the ABA-responsive element-binding protein family, MYB, NAC, the protein kinase superfamily, etc. These results revealed promising alleles linked to drought resistance or root traits, providing insights into the genetic basis of drought resistance and roots, which will be useful for common bean improvement.

Regulatory Mechanisms of the Resistance to Common Bacterial Blight Revealed by Transcriptomic Analysis in Common Bean (Phaseolus vulgaris L.).

Common bean blight (CBB), primarily caused by Xanthomonas axonopodis pv. phaseoli (Xap), is one of the most destructive diseases of common bean (Phaseolus vulgaris L.). The tepary bean genotype PI 319443 displays high resistance to Xap, and the common bean genotypes HR45 and Bilu display high resistance and susceptibility to Xap, respectively. To identify candidate genes related to Xap resistance, transcriptomic analysis was performed to compare gene expression levels with Xap inoculation at 0, 24, and 48 h post inoculation (hpi) among the three genotypes. A total of 1,146,009,876 high-quality clean reads were obtained. Differentially expressed gene (DEG) analysis showed that 1,688 DEGs responded to pathogen infection in the three genotypes. Weighted gene coexpression network analysis (WGCNA) was also performed to identify three modules highly correlated with Xap resistance, in which 334 DEGs were likely involved in Xap resistance. By combining differential expression analysis and WGCNA, 139 DEGs were identified as core resistance-responsive genes, including 18 genes encoding resistance (R) proteins, 19 genes belonging to transcription factor families, 63 genes encoding proteins with oxidoreductase activity, and 33 plant hormone signal transduction-related genes, which play important roles in the resistance to pathogen infection. The expression patterns of 20 DEGs were determined by quantitative real-time PCR (qRT-PCR) and confirmed the reliability of the RNA-seq results.

Marker-Trait Association Analysis of Seed Traits in Accessions of Common Bean (Phaseolus vulgaris L.) in China.

Seed weight and seed size are the key agronomic traits that determine yield in common bean. To investigate the genetic architecture of four seed traits (100-seed weight, seed length, seed width, and seed height) of common bean in China, marker-trait association analysis of these seed traits was performed in a nationwide population of 395 common bean accessions using 116 polymorphic SSR markers. The four seed traits were evaluated in six trials across three environments. Seed size varied among the environments. Population structure was evaluated based on SSR markers and phaseolin, which divided the accessions into two main subpopulations representing the two known gene pools. Seed weight and seed size had a strong relationship with population clustering. In addition, in a Genome-wide association studies (GWAS), 21 significantly associated markers were identified for the four seed traits with two models, namely, general linear model (GLM) and mixed linear model (MLM). Some markers had pleiotropic effects, i.e., controlled more than one trait. The significant quantitative trait loci identified in this study could be used in marker-assisted breeding to accelerate the genetic improvement of yield in common bean.

Resequencing of 683 common bean genotypes identifies yield component trait associations across a north-south cline.

We conducted a large-scale genome-wide association study evaluation of 683 common bean accessions, including landraces and breeding lines, grown over 3 years and in four environments across China, ranging in latitude from 18.23° to 45.75° N, with different planting dates and abiotic or biotic stresses. A total of 505 loci were associated with yield components, of which seed size, flowering time and harvest maturity traits were stable across years and environments. Some loci aligned with candidate genes controlling these traits. Yield components were observed to have strong associations with a gene-rich region on the long arm of chromosome 1. Manipulation of seed size, through selection of seed length versus seed width and height, was deemed possible, providing a genome-based means to select for important yield components. This study shows that evaluation of large germplasm collections across north-south geographic clines is useful in the detection of marker associations that determine grain yield in pulses.

Morphology independent triplet formation in pentalene films: Singlet fission as the triplet formation mechanism.

Singlet fission (SF), a spin-allowed multiexciton generation process, experienced renewed interest in the last decade due to its potential to increase the efficiency of photovoltaic devices. The hurdles now lie in the limited range of SF-capable materials and demanding morphology requirement for an efficient fission process. Although primary fission to yield triplet pair (1TT) can occur independently of film morphology in intramolecular singlet fission (iSF) materials, the separation of the 1TT state has been shown to be highly dependent on the packing motif and morphologies. In this work, we have demonstrated that both iSF and triplet pair separation processes took place irrelevant of molecular order and/or film morphology in a series of pentalene compounds. With the >180% fission efficiency, the suitable triplet energy levels, and the long lifetime of the triplet excitons, these iSF systems can be integrated into practical photovoltaic application.

Salt and drought stress and ABA responses related to bZIP genes from V. radiata and V. angularis.

Mung bean and adzuki bean are warm-season legumes widely cultivated in China. However, bean production in major producing regions is limited by biotic and abiotic stress, such as drought and salt stress. Basic leucine zipper (bZIP) genes play key roles in responses to various biotic and abiotic stresses. However, only several bZIP genes involved in drought and salt stress in legumes, especially Vigna radiata and Vigna angularis, have been identified. In this study, we identified 54 and 50 bZIP proteins from whole-genome sequences of V. radiata and V. angularis, respectively. First, we comprehensively surveyed the characteristics of all bZIP genes, including their gene structure, chromosome distribution and motif composition. Phylogenetic trees showed that VrbZIP and VabZIP proteins were divided into ten clades comprising nine known and one unknown subgroup. The results of the nucleotide substitution rate of the orthologous gene pairs showed that bZIP proteins have undergone strong purifying selection: V. radiata and V. angularis diverged 1.25 million years ago (mya) to 9.20 mya (average of 4.95 mya). We also found that many cis-acting regulatory elements (CAREs) involved in abiotic stress and plant hormone responses were detected in the putative promoter regions of the bZIP genes. Finally, using the quantitative real-time PCR (qRT-PCR) method, we performed expression profiling of the bZIP genes in response to drought, salt and abscisic acid (ABA). We identified several bZIP genes that may be involved in drought and salt responses. Generally, our results provided useful and rich resources of VrbZIP and VabZIP genes for the functional characterization and understanding of bZIP transcription factors (TFs) in warm-season legumes. In addition, our results revealed important and interesting data - a subset of VrbZIP and VabZIP gene expression profiles in response to drought, salt and ABA stress. These results provide gene expression evidence for the selection of candidate genes under drought and salt stress for future study.

Absence of Intramolecular Singlet Fission in Pentacene-Perylenediimide Heterodimers: The Role of Charge Transfer State.

A new class of donor-acceptor heterodimers based on two singlet fission (SF)-active chromophores, i.e., pentacene (Pc) and perylenediimide (PDI), was developed to investigate the role of charge transfer (CT) state on the excitonic dynamics. The CT state is efficiently generated upon photoexcitation. However, the resulting CT state decays to different energy states depending on the energy levels of the CT state. It undergoes extremely rapid deactivation to the ground state in polar CH2Cl2, whereas it undergoes transformation to a Pc triplet in nonpolar toluene. The efficient triplet generation in toluene is not due to SF but CT-mediated intersystem crossing. In light of the energy landscape, it is suggested that the deep energy level of the CT state relative to that of the triplet pair state makes the CT state actually serve as a trap state that cannot undergoes an intramolecular singlet fission process. These results provide guidance for the design of SF materials and highlight the requisite for more widely applicable design principles.

Quantitative trait locus mapping under irrigated and drought treatments based on a novel genetic linkage map in mungbean (Vigna radiata L.).

KEY MESSAGE: A novel genetic linkage map was constructed using SSR markers and stable QTLs were identified for six drought tolerance related-traits using single-environment analysis under irrigation and drought treatments. Mungbean (Vigna radiata L.) is one of the most important leguminous food crops. However, mungbean production is seriously constrained by drought. Isolation of drought-responsive genetic elements and marker-assisted selection breeding will benefit from the detection of quantitative trait locus (QTLs) for traits related to drought tolerance. In this study, we developed a full-coverage genetic linkage map based on simple sequence repeat (SSR) markers using a recombinant inbred line (RIL) population derived from an intra-specific cross between two drought-resistant varieties. This novel map was anchored with 313 markers. The total map length was 1010.18 cM across 11 linkage groups, covering the entire genome of mungbean with a saturation of one marker every 3.23 cM. We subsequently detected 58 QTLs for plant height (PH), maximum leaf area (MLA), biomass (BM), relative water content, days to first flowering, and seed yield (Yield) and 5 for the drought tolerance index of 3 traits in irrigated and drought environments at 2 locations. Thirty-eight of these QTLs were consistently detected two or more times at similar linkage positions. Notably, qPH5A and qMLA2A were consistently identified in marker intervals from GMES5773 to MUS128 in LG05 and from Mchr11-34 to the HAAS_VR_1812 region in LG02 in four environments, contributing 6.40-20.06% and 6.97-7.94% of the observed phenotypic variation, respectively. None of these QTLs shared loci with previously identified drought-related loci from mungbean. The results of these analyses might facilitate the isolation of drought-related genes and help to clarify the mechanism of drought tolerance in mungbean.

Genome-Wide Association Study Identifies NBS-LRR-Encoding Genes Related with Anthracnose and Common Bacterial Blight in the Common Bean.

Nucleotide-binding site and leucine-rich repeat (NBS-LRR) genes represent the largest and most important disease resistance genes in plants. The genome sequence of the common bean (Phaseolus vulgaris L.) provides valuable data for determining the genomic organization of NBS-LRR genes. However, data on the NBS-LRR genes in the common bean are limited. In total, 178 NBS-LRR-type genes and 145 partial genes (with or without a NBS) located on 11 common bean chromosomes were identified from genome sequences database. Furthermore, 30 NBS-LRR genes were classified into Toll/interleukin-1 receptor (TIR)-NBS-LRR (TNL) types, and 148 NBS-LRR genes were classified into coiled-coil (CC)-NBS-LRR (CNL) types. Moreover, the phylogenetic tree supported the division of these PvNBS genes into two obvious groups, TNL types and CNL types. We also built expression profiles of NBS genes in response to anthracnose and common bacterial blight using qRT-PCR. Finally, we detected nine disease resistance loci for anthracnose (ANT) and seven for common bacterial blight (CBB) using the developed NBS-SSR markers. Among these loci, NSSR24, NSSR73, and NSSR265 may be located at new regions for ANT resistance, while NSSR65 and NSSR260 may be located at new regions for CBB resistance. Furthermore, we validated NSSR24, NSSR65, NSSR73, NSSR260, and NSSR265 using a new natural population. Our results provide useful information regarding the function of the NBS-LRR proteins and will accelerate the functional genomics and evolutionary studies of NBS-LRR genes in food legumes. NBS-SSR markers represent a wide-reaching resource for molecular breeding in the common bean and other food legumes. Collectively, our results should be of broad interest to bean scientists and breeders.

MicroRNAs associated with drought response in the pulse crop common bean (Phaseolus vulgaris L.).

Drought stress significantly reduces common bean yields. Recently, some drought-related miRNAs were found in various plants. However, reports of miRNAs involved in drought stress in common bean are limited. Here, we obtained four sRNA samples from drought-tolerant and -sensitive cultivars of common bean that experienced with or without drought treatment. A total of 49 novel miRNAs and 120 known miRNAs were detected. Under drought treatment, 9 and 7 known miRNAs were down and up-regulated, respectively, and 5 and 3 of the novel miRNAs were increased and decreased, respectively. Among these miRNAs, four miRNAs shared the same pattern of expression between Long 22-0579 and Naihua. Target genes of these miRNAs included transcription factors, protein kinases, and nuclear transcription factors. Finally, we verified all of the differentially expressed miRNAs by RT-qPCR, and we identified 16 miRNAs that are potentially associated with the drought stress response. These miRNAs and target genes will be useful in future basic studies and in applied studies investigating how miRNA regulation can be used to enhance drought resistance in plant species.

Hairy root transgene expression analysis of a secretory peroxidase (PvPOX1) from common bean infected by Fusarium wilt.

Plant peroxidases (POXs) are one of the most important redox enzymes in the defense responses. However, the large number of different plant POX genes makes it necessary to carefully confirm the function of each paralogous POX gene in specific tissues and disease interactions. Fusarium wilt is a devastating disease of common bean caused by Fusarium oxysporum f. sp. phaseoli. In this study, we evaluated a peroxidase gene, PvPOX1, from a resistant common bean genotype, CAAS260205 and provided direct evidence for PvPOX1's role in resistance by transforming the resistant allele into a susceptible common bean genotype, BRB130, via hairy root transformation using Agrobacterium rhizogenes. Analysis of PvPOX1 gene over-expressing hairy roots showed it increased resistance to Fusarium wilt both in the roots and the rest of transgenic plants. Meanwhile, the PvPOX1 expressive level, the peroxidase activity and hydrogen peroxide (H2O2) accumulation were also enhanced in the interaction. The result showed that the PvPOX1 gene played an essential role in Fusarium wilt resistance through the occurrence of reactive oxygen species (ROS) induced hypersensitive response. Therefore, PvPOX1 expression was proven to be a valuable gene for further analysis which can strengthen host defense response against Fusarium wilt through a ROS activated resistance mechanism.

Genome-Wide Investigation of WRKY Transcription Factors Involved in Terminal Drought Stress Response in Common Bean.

WRKY transcription factor plays a key role in drought stress. However, the characteristics of the WRKY gene family in the common bean (Phaseolus vulgaris L.) are unknown. In this study, we identified 88 complete WRKY proteins from the draft genome sequence of the "G19833" common bean. The predicted genes were non-randomly distributed in all chromosomes. Basic information, amino acid motifs, phylogenetic tree and the expression patterns of PvWRKY genes were analyzed, and the proteins were classified into groups 1, 2, and 3. Group 2 was further divided into five subgroups: 2a, 2b, 2c, 2d, and 2e. Finally, we detected 19 WRKY genes that were responsive to drought stress using qRT-PCR; 11 were down-regulated, and 8 were up-regulated under drought stress. This study comprehensively examines WRKY proteins in the common bean, a model food legume, and it provides a foundation for the functional characterization of the WRKY family and opportunities for understanding the mechanisms of drought stress tolerance in this plant.

The iSelect 9 K SNP analysis revealed polyploidization induced revolutionary changes and intense human selection causing strong haplotype blocks in wheat.

A Chinese wheat mini core collection was genotyped using the wheat 9 K iSelect SNP array. Total 2420 and 2396 polymorphic SNPs were detected on the A and the B genome chromosomes, which formed 878 haplotype blocks. There were more blocks in the B genome, but the average block size was significantly (P < 0.05) smaller than those in the A genome. Intense selection (domestication and breeding) had a stronger effect on the A than on the B genome chromosomes. Based on the genetic pedigrees, many blocks can be traced back to a well-known Strampelli cross, which was made one century ago. Furthermore, polyploidization of wheat (both tetraploidization and hexaploidization) induced revolutionary changes in both the A and the B genomes, with a greater increase of gene diversity compared to their diploid ancestors. Modern breeding has dramatically increased diversity in the gene coding regions, though obvious blocks were formed on most of the chromosomes in both tetraploid and hexaploid wheats. Tag-SNP markers identified in this study can be used for marker assisted selection using haplotype blocks as a wheat breeding strategy. This strategy can also be employed to facilitate genome selection in other self-pollinating crop species.