Marker-Assisted Molecular Profiling, Deletion Mutant Analysis, and RNA-Seq Reveal a Disease Resistance Cluster Associated with Uromyces appendiculatus Infection in Common Bean Phaseolus vulgaris L.

Common bean (Phaseolus vulgaris L.) is an important legume, useful for its high protein and dietary fiber. The fungal pathogen Uromyces appendiculatus (Pers.) Unger can cause major loss in susceptible varieties of the common bean. The Ur-3 locus provides race specific resistance to virulent strains or races of the bean rust pathogen along with Crg, (Complements resistance gene), which is required for Ur-3-mediated rust resistance. In this study, we inoculated two common bean genotypes (resistant "Sierra" and susceptible crg) with rust race 53 of U. appendiculatus, isolated leaf RNA at specific time points, and sequenced their transcriptomes. First, molecular markers were used to locate and identify a 250 kb deletion on chromosome 10 in mutant crg (which carries a deletion at the Crg locus). Next, we identified differential expression of several disease resistance genes between Mock Inoculated (MI) and Inoculated (I) samples of "Sierra" leaf RNA within the 250 kb delineated region. Both marker assisted molecular profiling and RNA-seq were used to identify possible transcriptomic locations of interest regarding the resistance in the common bean to race 53. Identification of differential expression among samples in disease resistance clusters in the bean genome may elucidate significant genes underlying rust resistance. Along with preserving favorable traits in the crop, the current research may also aid in global sustainability of food stocks necessary for many populations.

New genomic regions associated with white mold resistance in dry bean using a MAGIC population.

Dry bean (Phaseolus vulgaris L.) production in many regions is threatened by white mold (WM) [Sclerotinia sclerotiorum (Lib.) de Bary]. Seed yield losses can be up to 100% under conditions favorable for the pathogen. The low heritability, polygenic inheritance, and cumbersome screening protocols make it difficult to breed for improved genetic resistance. Some progress in understanding genetic resistance and germplasm improvement has been accomplished, but cultivars with high levels of resistance are yet to be released. A WM multiparent advanced generation inter-cross (MAGIC) population (n = 1060) was developed to facilitate mapping and breeding efforts. A seedling straw test screening method provided a quick assay to phenotype the population for response to WM isolate 1980. Nineteen MAGIC lines were identified with improved resistance. For genome-wide association studies (GWAS), the data was transformed into three phenotypic distributions-quantitative, polynomial, and binomial-and coupled with ∼52,000 single-nucleotide polymorphisms (SNPs). The three phenotypic distributions identified 30 significant genomic intervals [-log10 (P value) ≥ 3.0]. However, across distributions, four new genomic regions as well as two regions previously reported were found to be associated with resistance. Cumulative R2 values were 57% for binomial distribution using 13 genomic intervals, 41% for polynomial using eight intervals, and 40% for quantitative using 11 intervals. New resistant germplasm as well as new genomic regions associated with resistance are now available for further investigation.

High Level of Nonsynonymous Changes in Common Bean Suggests That Selection under Domestication Increased Functional Diversity at Target Traits.

Crop species have been deeply affected by the domestication process, and there have been many efforts to identify selection signatures at the genome level. This knowledge will help geneticists to better understand the evolution of organisms, and at the same time, help breeders to implement successful breeding strategies. Here, we focused on domestication in the Mesoamerican gene pool of Phaseolus vulgaris by sequencing 49 gene fragments from a sample of 45 P. vulgaris wild and domesticated accessions, and as controls, two accessions each of the closely related species Phaseolus coccineus and Phaseolus dumosus. An excess of nonsynonymous mutations within the domesticated germplasm was found. Our data suggest that the cost of domestication alone cannot explain fully this finding. Indeed, the significantly higher frequency of polymorphisms in the coding regions observed only in the domesticated plants (compared to noncoding regions), the fact that these mutations were mostly nonsynonymous and appear to be recently derived mutations, and the investigations into the functions of their relative genes (responses to biotic and abiotic stresses), support a scenario that involves new functional mutations selected for adaptation during domestication. Moreover, consistent with this hypothesis, selection analysis and the possibility to compare data obtained for the same genes in different studies of varying sizes, data types, and methodologies allowed us to identify four genes that were strongly selected during domestication. Each selection candidate is involved in plant resistance/tolerance to abiotic stresses, such as heat, drought, and salinity. Overall, our study suggests that domestication acted to increase functional diversity at target loci, which probably controlled traits related to expansion and adaptation to new agro-ecological growing conditions.

Genome-wide association studies identifies seven major regions responsible for iron deficiency chlorosis in soybean (Glycine max).

Iron deficiency chlorosis (IDC) is a yield limiting problem in soybean (Glycine max (L.) Merr) production regions with calcareous soils. Genome-wide association study (GWAS) was performed using a high density SNP map to discover significant markers, QTL and candidate genes associated with IDC trait variation. A stepwise regression model included eight markers after considering LD between markers, and identified seven major effect QTL on seven chromosomes. Twelve candidate genes known to be associated with iron metabolism mapped near these QTL supporting the polygenic nature of IDC. A non-synonymous substitution with the highest significance in a major QTL region suggests soybean orthologs of FRE1 on Gm03 is a major gene responsible for trait variation. NAS3, a gene that encodes the enzyme nicotianamine synthase which synthesizes the iron chelator nicotianamine also maps to the same QTL region. Disease resistant genes also map to the major QTL, supporting the hypothesis that pathogens compete with the plant for Fe and increase iron deficiency. The markers and the allelic combinations identified here can be further used for marker assisted selection.

Genetic architecture of chalcone isomerase non-coding regions in common bean (Phaseolus vulgaris L.).

Sequence data for 2 non-coding regions of the chalcone isomerase gene were analyzed to study the genetic architecture of common bean (Phaseolus vulgaris L.). One region corresponded to the first 596 nucleotides (nt) of the 5'-untranslated region (UTR). The other region was the 710 nt intron 3. Data were collected from 67 genotypes representing both landraces and cultivars from the geographical range of the cultivated form of the species. Variability in the 5'-UTR region was represented by single nucleotide polymorphisms (SNPs), whereas intron 3 variation was due to a collection of SNPs and insertion-deletion events. Diversity was greater in the 5'-UTR (pi = 0.0175) than in intron 3 (pi = 0.0089). For each region, diversity was greater for genotypes of Middle American than Andean origin. A single recombination event was observed, and the hybridization pattern necessary to derive the recombinant genotypes supported the previous observation of an ancestral gene pool from which modern domesticated genotypes are derived. For both regions, a strongly supported Andean group was observed, whereas the presence of 2 Middle American subgroups was also supported. Although a significantly positive Tajima's D statistic was observed for the 5'-UTR for all genotypes, we conclude that this is more likely the result of a strong demographic effect and that balancing selection is occurring only among the Middle American genotypes.

Sequence diversity analysis of dihydroflavonol 4-reductase intron 1 in common bean.

Variation in common bean (Phaseolus vulgaris L.) was investigated by sequencing intron 1 of the dihydroflavonol 4-reductase (DFR) gene for 92 genotypes that represent both landraces and cultivars. We were also interested in determining if introns provide sufficient variation for genetic diversity studies and if the sequence data could be used to develop allele-specific primers that could differentiate genotypes using a standard PCR assay. Sixty-nine polymorphic sites were observed. Nucleotide variation (pi/bp) was 0.0481, a value higher than that reported for introns from other plant species. Tests for significant deviation from the mutation drift model were positive for the population as a whole, the cultivar and landrace subsets, and the Middle American landrace set. Significant linkage disequilibrium extended about 300 nucleotides. Twenty haplotypes were detected among the cultivated genotypes. Seven recombination events were detected for the whole population, and six events for the landraces. Recombination was not observed among the landraces within either the Middle American or Andean gene pools. Evidence for hybridization between the two gene pools was discovered. Five allele-specific primers were developed that could distinguish 56 additional genotypes. The allele-specific primers were used to map duplicate DFR genes on linkage group B8.

A molecular marker tightly linked to P, a gene required for flower and seedcoat color in common bean (Phaseolus vulgaris L.), contains the Ty3-gypsy retrotransposon Tpv3g.

In common bean (Phaseolus vulgaris L.), the expression of color in flower and seedcoat tissues requires the dominant allele of the P gene. The fully recessive p allele completely suppresses color expression in these tissues, whereas in specific genetic backgrounds, the p(gri) allele potentiates a grayish-white seedcoat and pale violet (nearly white) flowers with two violet dots on the banner petals. As a first step to gaining a better understanding of this important gene, we phenotypically scored an F2 population segregating for P and p(gri) and subsequently screened contrasting bulk DNA samples with oligonucleotide primers to uncover random amplified polymorphic DNA (RAPD) fragments. OU3(2300), an RAPD marker linked in coupling phase to the dominant allele, mapped 1.3 cM from P. The core 'BAT93' x 'Jalo EEP558' recombinant inbred population was scored, and the marker mapped to linkage group B7. The segregating fragment was cloned, sequenced, and shown to possess significant homology to the Ty3-gypsy class of retrotransposons. We have named the element Tpv3g. It is estimated that about 100 copies of the element are present in the common bean genome. Phylogenetic analysis placed Tpv3g in the class A group of plant retrotransposons.