Evolution, diversity, and function of the disease susceptibility gene Snn1 in wheat.

Septoria nodorum blotch (SNB), caused by Parastagonospora nodorum, is a disease of durum and common wheat initiated by the recognition of pathogen-produced necrotrophic effectors (NEs) by specific wheat genes. The wheat gene Snn1 was previously cloned, and it encodes a wall-associated kinase that directly interacts with the NE SnTox1 leading to programmed cell death and ultimately the development of SNB. Here, sequence analysis of Snn1 from 114 accessions including diploid, tetraploid, and hexaploid wheat species revealed that some wheat lines possess two copies of Snn1 (designated Snn1-B1 and Snn1-B2) approximately 120 kb apart. Snn1-B2 evolved relatively recently as a paralog of Snn1-B1, and both genes have undergone diversifying selection. Three point mutations associated with the formation of the first SnTox1-sensitive Snn1-B1 allele from a primitive wild wheat were identified. Four subsequent and independent SNPs, three in Snn1-B1 and one in Snn1-B2, converted the sensitive alleles to insensitive forms. Protein modeling indicated these four mutations could abolish Snn1-SnTox1 compatibility either through destabilization of the Snn1 protein or direct disruption of the protein-protein interaction. A high-throughput marker was developed for the absent allele of Snn1, and it was 100% accurate at predicting SnTox1-insensitive lines in both durum and spring wheat. Results of this study increase our understanding of the evolution, diversity, and function of Snn1-B1 and Snn1-B2 genes and will be useful for marker-assisted elimination of these genes for better host resistance.

Genetic mapping and genomic prediction of sclerotinia stem rot resistance to rapeseed/canola (Brassica napus L.) at seedling stage.

KEY MESSAGE: GWAS detected ninety-eight significant SNPs associated with Sclerotinia sclerotiorum resistance. Six statistical models resulted in medium to high predictive ability, depending on trait, indicating potential of genomic prediction for disease resistance breeding. The lack of complete host resistance and a complex resistance inheritance nature between rapeseed/canola and Sclerotinia sclerotiorum often limits the development of functional molecular markers that enable breeding for sclerotinia stem rot (SSR) resistance. However, genomics-assisted selection has the potential to accelerate the breeding for SSR resistance. Therefore, genome-wide association (GWA) mapping and genomic prediction (GP) were performed using a diverse panel of 337 rapeseed/canola genotypes. Three-week-old seedlings were screened using the petiole inoculation technique (PIT). Days to wilt (DW) up to 2 weeks and lesion phenotypes (LP) at 3, 4, and 7 days post-inoculation (dpi) were recorded. A strong correlation (r = - 0.90) between DW and LP_4dpi implied that a single time point scoring at four days could be used as a proxy trait. GWA analyses using single-locus (SL) and multi-locus (ML) models identified a total of 41, and 208 significantly associated SNPs, respectively. Out of these, ninety-eight SNPs were identified by a combination of the SL model and any of the ML models, at least two ML models, or two traits. These SNPs explained 1.25-12.22% of the phenotypic variance and considered as significant, could be associated with SSR resistance. Eighty-three candidate genes with a function in disease resistance were associated with the significant SNPs. Six GP models resulted in moderate to high (0.42-0.67) predictive ability depending on SSR resistance traits. The resistant genotypes and significant SNPs will serve as valuable resources for future SSR resistance breeding. Our results also highlight the potential of genomic selection to improve rapeseed/canola breeding for SSR resistance.

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 genetics and physiology of seed dormancy, a crucial trait in common bean domestication.

BACKGROUND: Physical seed dormancy is an important trait in legume domestication. Although seed dormancy is beneficial in wild ecosystems, it is generally considered to be an undesirable trait in crops due to reduction in yield and / or quality. The physiological mechanism and underlying genetic factor(s) of seed dormancy is largely unknown in several legume species. Here we employed an integrative approach to understand the mechanisms controlling physical seed dormancy in common bean (Phaseolus vulgaris L.). RESULTS: Using an innovative CT scan imaging system, we were able to track water movements inside the seed coat. We found that water uptake initiates from the bean seed lens. Using a scanning electron microscopy (SEM) we further identified several micro-cracks on the lens surface of non-dormant bean genotypes. Bulked segregant analysis (BSA) was conducted on a bi-parental RIL (recombinant inbred line) population, segregating for seed dormancy. This analysis revealed that the seed water uptake is associated with a single major QTL on Pv03. The QTL region was fine-mapped to a 118 Kb interval possessing 11 genes. Coding sequence analysis of candidate genes revealed a 5-bp insertion in an ortholog of pectin acetylesterase 8 that causes a frame shift, loss-of-function mutation in non-dormant genotype. Gene expression analysis of the candidate genes in the seed coat of contrasting genotypes indicated 21-fold lower expression of pectin acetylesterase 8 in non-dormant genotype. An analysis of mutational polymorphism was conducted among wild and domesticated beans. Although all the wild beans possessed the functional allele of pectin acetylesterase 8, the majority (77%) of domesticated beans had the non-functional allele suggesting that this variant was under strong selection pressure through domestication. CONCLUSIONS: In this study, we identified the physiological mechanism of physical seed dormancy and have identified a candidate allele causing variation in this trait. Our findings suggest that a 5-bp insertion in an ortholog of pectin acetylesterase 8 is likely a major causative mutation underlying the loss of seed dormancy during domestication. Although the results of current study provide strong evidences for the role of pectin acetylesterase 8 in seed dormancy, further confirmations seem necessary by employing transgenic approaches.

Computational identification of receptor-like kinases "RLK" and receptor-like proteins "RLP" in legumes.

BACKGROUND: In plants, the plasma membrane is enclosed by the cell wall and anchors RLK and RLP proteins, which play a fundamental role in perception of developmental and environmental cues and are crucial in plant development and immunity. These plasma membrane receptors belong to large gene/protein families that are not easily classified computationally. This detailed analysis of these plasma membrane proteins brings a new source of information to the legume genetic, physiology and breeding research communities. RESULTS: A computational approach to identify and classify RLK and RLP proteins is presented. The strategy was evaluated using experimentally-validated RLK and RLP proteins and was determined to have a sensitivity of over 0.85, a specificity of 1.00, and a Matthews correlation coefficient of 0.91. The computational approach can be used to develop a detailed catalog of plasma membrane receptors (by type and domains) in several legume/crop species. The exclusive domains identified in legumes for RLKs are WaaY, APH Pkinase_C, LRR_2, and EGF, and for RLP are L-lectin LPRY and PAN_4. The RLK-nonRD and RLCK subclasses are also discovered by the methodology. In both classes, less than 20% of the total RLK predicted for each species belong to this class. Among the 10-species evaluated ~ 40% of the proteins in the kinome are RLKs. The exclusive legume domain combinations identified are B-Lectin/PR5K domains in G. max, M. truncatula, V. angularis, and V. unguiculata and a three-domain combination B-lectin/S-locus/WAK in C. cajan, M. truncatula, P. vulgaris, V. angularis. and V. unguiculata. CONCLUSIONS: The analysis suggests that about 2% of the proteins of each genome belong to the RLK family and less than 1% belong to RLP family. Domain diversity combinations are greater for RLKs compared with the RLP proteins and LRR domains, and the dual domain combination LRR/Malectin were the most frequent domain for both groups of plasma membrane receptors among legume and non-legume species. Legumes exclusively show Pkinase extracellular domains, and atypical domain combinations in RLK and RLP compared with the non-legumes evaluated. The computational logic approach is statistically well supported and can be used with the proteomes of other plant species.

Genetic analysis of threshability and other spike traits in the evolution of cultivated emmer to fully domesticated durum wheat.

Genetic mutations in genes governing wheat threshability were critical for domestication. Knowing when these genes mutated during wheat evolution will provide more insight into the domestication process and lead to further exploitation of primitive alleles for wheat improvement. We evaluated a population of recombinant inbred lines derived from a cross between the durum variety Rusty and the cultivated emmer accession PI 193883 for threshability, rachis fragility, and other spike-related traits. Quantitative trait loci (QTL) associated with spike length, spikelets per spike, and spike compactness were primarily associated with known genes such as the pleiotropic domestication gene Q. Interestingly, rachis fragility was not associated with the Q locus, suggesting that this trait, usually a pleiotropic effect of the q allele, can be influenced by the genetic background. Threshability QTL were identified on chromosome arms 2AS, 2BS, and 5AL corresponding to the tenacious glume genes Tg2A and Tg2B as well as the Q gene, respectively, further demonstrating that cultivated emmer harbors the primitive non-free-threshing alleles at all three loci. Genetic analysis indicated that the effects of the three genes are mostly additive, with Q having the most profound effects on threshability, and that free-threshing alleles are necessary at all three loci to attain a completely free-threshing phenotype. These findings provide further insight into the timeline and possible pathways of wheat domestication and evolution that led to the formation of modern day domesticated wheats.

Mapping and characterization of two stem rust resistance genes derived from cultivated emmer wheat accession PI 193883.

KEY MESSAGE: Two stem rust resistance genes identified on chromosome arms 2BL and 6AL of the cultivated emmer wheat accession PI 193883 can be used for protecting modern varieties against Ug99 strains. The wheat research community consistently strives to identify new genes that confer resistance to stem rust caused by the fungal pathogen Puccinia graminis f. sp. tritici Eriks & E. Henn (Pgt). In the current study, our objective was to identify and genetically characterize the stem rust resistance derived from the cultivated emmer accession PI 193883. A recombinant inbred line population developed from a cross between the susceptible durum wheat line Rusty and PI 193883 was genotyped and evaluated for reaction to Pgt races TTKSK, TRTTF, and TMLKC. Two QTLs conferring resistance were identified on chromosome arms 2BL (QSr.fcu-2B) and 6AL (QSr.fcu-6A). The stem rust resistance gene (Sr883-2B) underlying QSr.fcu-2B was recessive, and based on its physical location it is located proximal to the Sr9 region. QSr.fcu-6A was located in the Sr13 region, but PI 193883 is known to carry the susceptible haplotype S4 for Sr13, indicating that the gene underlying QSr.fcu-6A (Sr883-6A) is likely a new allele of Sr13 or a gene residing close to Sr13. Three IWGSC scaffold-based simple sequence repeat (SSR) and two SNP-based semi-thermal asymmetric reverse PCR (STARP) markers were developed for the Sr883-2B region, and one STARP marker was developed for Sr883-6A. Sr883-2B was epistatic to Sr883-6A for reaction to TTKSK and TRTTF, and the two genes had additive effects for TMLKC. These two genes and the markers developed in this research provide additional resources and tools for the improvement in stem rust resistance in durum and common wheat breeding programs.

White seed color in common bean (Phaseolus vulgaris) results from convergent evolution in the P (pigment) gene.

The presence of seed color in common bean (Phaseolus vulgaris) requires the dominant-acting P (pigment) gene, and white seed is a recessive phenotype in all domesticated races of the species. P was classically associated with seed size, thus describing it as the first genetic marker for a quantitative trait. The molecular structure of P was characterized to understand the selection of white seeds during bean diversification and the relationship of P to seed weight. P was identified by homology searches, a genome-wide association study (GWAS) and gene remodeling, and confirmed by gene silencing. Allelic variation was assessed by a combination of resequencing and marker development, and the relationship between P and seed weight was assessed by a GWAS study. P is a member of clade B of subclass IIIf of plant basic helix-loop-helix (bHLH) proteins. Ten race-specific P alleles conditioned the white seed phenotype, and each causative mutation affected at least one bHLH domain required for color expression. GWAS analysis confirmed the classic association of P with seed weight. In common bean, white seeds are the result of convergent evolution and, among plant species, orthologous convergence on a single transcription factor gene was observed.

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 QTL alleles for quality-related traits in spring wheat revealed by RIL population derived from supernumerary × non-supernumerary spikelet genotypes.

KEY MESSAGE: A population developed from an exotic line with supernumerary spikelets was genetically dissected for eight quality traits, discovering new genes/alleles with potential use in wheat breeding programs. Identifying new QTLs and alleles in exotic germplasm is paramount for further improvement of quality traits in wheat. In the present study, an RIL population developed from a cross of an elite wheat line (WCB414) and an exotic genotype with supernumerary spikelets (SS) was used to identify QTLs and new alleles for eight quality traits. Composite interval mapping for 1,000 kernels weight (TKW), kernel volume weight (KVW), grain protein content (GPC), percent of flour extraction (FE) and four mixograph-related traits identified a total of 69 QTLs including 19 stable QTLs. These QTLs were located on 18 different chromosomes (except 4D, 5D, and 6D). Thirteen of these QTLs explained more than 15% of phenotypic variation (PV) and were considered as major QTLs. In this study, we identified 11 QTLs for TKW (R (2) = 7.2-17.1 %), 10 for KVW (R (2) = 6.7-22.5%), 11 for GPC (R (2) = 4.7-16.9%), 6 for FE (R (2) = 4.8-19%) and 31 for mixograph-related traits (R (2) = 3.2-41.2%). In this population, several previously identified QTLs for SS, nine spike-related and ten agronomic traits were co-located with the quality QTLs, suggesting pleiotropic effects or close linkage among loci. The traits GPC and mixogram-related traits were positively correlated with SS. Indeed, several loci for quality traits were co-located with QTL for SS. The exotic parent contributed positive alleles that increased PV of the traits at 56% of loci demonstrating the suitability of germplasm with SS to improve quality traits in wheat.

Application of in silico bulked segregant analysis for rapid development of markers linked to Bean common mosaic virus resistance in common bean.

BACKGROUND: Common bean was one of the first crops that benefited from the development and utilization of molecular marker-assisted selection (MAS) for major disease resistance genes. Efficiency of MAS for breeding common bean is still hampered, however, due to the dominance, linkage phase, and loose linkage of previously developed markers. Here we applied in silico bulked segregant analysis (BSA) to the BeanCAP diversity panel, composed of over 500 lines and genotyped with the BARCBEAN_3 6K SNP BeadChip, to develop codominant and tightly linked markers to the I gene controlling resistance to Bean common mosaic virus (BCMV). RESULTS: We physically mapped the genomic region underlying the I gene. This locus, in the distal arm of chromosome Pv02, contains seven putative NBS-LRR-type disease resistance genes. Two contrasting bulks, containing BCMV host differentials and ten BeanCAP lines with known disease reaction to BCMV, were subjected to in silico BSA for targeting the I gene and flanking sequences. Two distinct haplotypes, containing a cluster of six single nucleotide polymorphisms (SNP), were associated with resistance or susceptibility to BCMV. One-hundred and twenty-two lines, including 115 of the BeanCAP panel, were screened for BCMV resistance in the greenhouse, and all of the resistant or susceptible plants displayed distinct SNP haplotypes as those found in the two bulks. The resistant/susceptible haplotypes were validated in 98 recombinant inbred lines segregating for BCMV resistance. The closest SNP (~25-32 kb) to the distal NBS-LRR gene model for the I gene locus was targeted for conversion to codominant KASP (Kompetitive Allele Specific PCR) and CAPS (Cleaved Amplified Polymorphic Sequence) markers. Both marker systems accurately predicted the disease reaction to BCMV conferred by the I gene in all screened lines of this study. CONCLUSIONS: We demonstrated the utility of the in silico BSA approach using genetically diverse germplasm, genotyped with a high-density SNP chip array, to discover SNP variation at a specific targeted genomic region. In common bean, many disease resistance genes are mapped and their physical genomic position can now be determined, thus the application of this approach will facilitate further development of codominant and tightly linked markers for use in MAS.

The alleles bc-ud and bc-ur (previously bc-4 gene), representing coding mutations within Vps4 AAA+ ATPase ESCRT protein, interact with other genes to condition resistance to BCMV and BCMNV in common bean.

Bean common mosaic virus (BCMV) and bean common mosaic necrosis virus (BCMNV) have a damaging impact on global common bean (Phaseolus vulgaris L.) cultivation, causing potential yield losses of over 80%. The primary strategy for controlling these viruses is through host plant resistance. This research aimed to identify and validate structural variations for the bc-ud gene as revealed by long-read sequencing, develop an efficient DNA marker to assist selection of bc-ud in snap and dry beans, and examine the interactions between the bc-ud allele and other BCMV resistance genes. A gene (Phvul.005G125100) model on chromosome Pv05, encoding a vacuolar protein-sorting 4 (Vps4) AAA+ ATPase endosomal sorting complexes required for transport (ESCRT) protein, was identified as the best candidate gene for bc-ud. An 84-bp repetitive insertion variant within the gene, exhibited 100% co-segregation with the bc-ud resistance allele across 264 common bean accessions. The 84-bp repetitive insertion was labeled with an indel marker IND_05_36225873 which was useful for tracking the bc-ud allele across diverse germplasm. A different single nucleotide polymorphism variant within the same candidate gene was associated with the bc-4 gene. Segregation in F2 populations confirmed bc-ud and bc-4 were alleles, so bc-4 was renamed bc-ur to fit gene nomenclature guidelines. The interactions of bc-ud and bc-ur with other resistance genes, such as bc-1 (receptor-like kinase on Pv03) and bc-2 (Vps4 AAA+ ATPase ESCRT protein on Pv11), validated gene combinations in the differential "host groups" effective against specific BCMV/BCMNV "pathogroups." These findings increase our understanding of the Bc-u locus, and enhance our ability to develop more resilient bean varieties through marker-assisted selection, reducing the impact of BCMV and BCMNV.

T and Z, partial seed coat patterning genes in common bean, provide insight into the structure and protein interactions of a plant MBW complex.

Flavonoids are secondary metabolites associated with plant seed coat and flower color. These compounds provide health benefits to humans as anti-inflammatory and antioxidant compounds. The expression of the late biosynthetic genes in the flavonoid pathway is controlled by a ternary MBW protein complex consisting of interfacing MYB, beta-helix-loop-helix (bHLH), and WD40 Repeat (WDR) proteins. P, the master regulator gene of the flavonoid expression in common bean (Phaseolus vulgaris L.), was recently determined to encode a bHLH protein. The T and Z genes control the distribution of color in bean seeds and flowers and have historically been considered regulators of the flavonoid gene expression. T and Z candidates were identified using reverse genetics based on genetic mapping, phylogenetic analysis, and mutant analysis. Domain and AlphaFold2 structure analyses determined that T encodes a seven-bladed ß-propeller WDR protein, while Z encodes a R2R3 MYB protein. Deletions and SNPs in T and Z mutants, respectively, altered the 3D structure of these proteins. Modeling of the Z MYB/P bHLH/T WDR MBW complex identified interfacing sequence domains and motifs in all three genes that are conserved in dicots. One Z MYB motif is a possible beta-molecular recognition feature (ß-MoRF) that only appears in a structured state when Z MYB is modeled as a component of a MBW complex. Complexes containing mutant T and Z proteins changed the interaction of members of the complex in ways that would alter their role in regulating the expression of genes in the flavonoid pathway.

GWAS of resistance to three bacterial diseases in the Andean common bean diversity panel.

Bacterial brown spot (BBS) caused by Pseudomonas syringae pv. syringae (Pss), common bacterial blight (CBB) caused by Xanthomonas axonopodis pv. phaseoli (Xap) and Xanthomonas fuscans subsp. fuscans (Xff), and halo bacterial blight (HBB), caused by Pseudomonas syringae pv. phaseolicola (Psph), are major bacterial diseases that severely affect common bean yields and global food security. Andean-origin dry beans, representing large-seeded market classes, are particularly susceptible. Using 140,325 SNPs, a multi-locus GWAS was conducted on subsets of the Andean diversity panel (ADP) phenotyped for BBS in South Africa, CBB in Puerto Rico, South Africa, and Zambia, and HBB in South Africa, through natural infection, artificial inoculation, or both. Twenty-four QTL associated with resistance were identified: nine for BBS, eight for CBB, and seven for HBB. Four QTL intervals on Pv01, Pv03, Pv05, and Pv08 overlapped with BBS and HBB resistance. A genomic interval on Pv01, near the fin gene, which determines growth habit, was linked to resistance to all three pathogens. Different QTLs were detected for BBS and CBB resistance when phenotyped under natural infection versus artificial inoculation. These results underscore the importance of combining phenotyping methods in multi-GWAS to capture the full genetic spectrum. Previously recognized CBB resistance QTL SAP6 and SU91 and HBB resistance QTL HB4.2, and HB5.1, were observed. Other common (MAF >0.25) and rare (MAF <0.05) resistance QTL were also detected. Overall, these findings enhance the understanding and utilization of bacterial resistance present in ADP for the development of common beans with improved resistance.

Artificial selection for determinate growth habit in soybean.

Determinacy is an agronomically important trait associated with the domestication in soybean (Glycine max). Most soybean cultivars are classifiable into indeterminate and determinate growth habit, whereas Glycine soja, the wild progenitor of soybean, is indeterminate. Indeterminate (Dt1/Dt1) and determinate (dt1/dt1) genotypes, when mated, produce progeny that segregate in a monogenic pattern. Here, we show evidence that Dt1 is a homolog (designated as GmTfl1) of Arabidopsis terminal flower 1 (TFL1), a regulatory gene encoding a signaling protein of shoot meristems. The transition from indeterminate to determinate phenotypes in soybean is associated with independent human selections of four distinct single-nucleotide substitutions in the GmTfl1 gene, each of which led to a single amino acid change. Genetic diversity of a minicore collection of Chinese soybean landraces assessed by simple sequence repeat (SSR) markers and allelic variation at the GmTfl1 locus suggest that human selection for determinacy took place at early stages of landrace radiation. The GmTfl1 allele introduced into a determinate-type (tfl1/tfl1) Arabidopsis mutants fully restored the wild-type (TFL1/TFL1) phenotype, but the Gmtfl1 allele in tfl1/tfl1 mutants did not result in apparent phenotypic change. These observations indicate that GmTfl1 complements the functions of TFL1 in Arabidopsis. However, the GmTfl1 homeolog, despite its more recent divergence from GmTfl1 than from Arabidopsis TFL1, appears to be sub- or neo-functionalized, as revealed by the differential expression of the two genes at multiple plant developmental stages and by allelic analysis at both loci.

Linked candidate genes of different functions for white mold resistance in common bean (Phaseolus vulgaris L) are identified by multiple QTL mapping approaches.

White mold (WM) is a major disease in common bean (Phaseolus vulgaris L.), and its complex quantitative genetic control limits the development of WM resistant cultivars. WM2.2, one of the nine meta-QTL with a major effect on WM tolerance, explains up to 35% of the phenotypic variation and was previously mapped to a large genomic interval on Pv02. Our objective was to narrow the interval of this QTL using combined approach of classic QTL mapping and QTL-based bulk segregant analysis (BSA), and confirming those results with Khufu de novo QTL-seq. The phenotypic and genotypic data from two RIL populations, 'Raven'/I9365-31 (R31) and 'AN-37'/PS02-029C-20 (Z0726-9), were used to select resistant and susceptible lines to generate subpopulations for bulk DNA sequencing. The QTL physical interval was determined by considering overlapping interval of the identified QTL or peak region in both populations by three independent QTL mapping analyses. Our findings revealed that meta-QTL WM2.2 consists of three regions, WM2.2a (4.27-5.76 Mb; euchromatic), WM 2.2b (12.19 to 17.61 Mb; heterochromatic), and WM2.2c (23.01-25.74 Mb; heterochromatic) found in both populations. Gene models encoding for gibberellin 2-oxidase 8, pentatricopeptide repeat, and heat-shock proteins are the likely candidate genes associated with WM2.2a resistance. A TIR-NBS-LRR class of disease resistance protein (Phvul.002G09200) and LRR domain containing family proteins are potential candidate genes associated with WM2.2b resistance. Nine gene models encoding disease resistance protein [pathogenesis-related thaumatin superfamily protein and disease resistance-responsive (dirigent-like protein) family protein etc] found within the WM2.2c QTL interval are putative candidate genes. WM2.2a region is most likely associated with avoidance mechanisms while WM2.2b and WM2.2c regions trigger physiological resistance based on putative candidate genes.

Integrating de novo QTL-seq and linkage mapping to identify quantitative trait loci conditioning physiological resistance and avoidance to white mold disease in dry bean.

White mold (WM), caused by the ubiquitous fungus Sclerotinia sclerotiorum, is a devastating disease that limits production and quality of dry bean globally. In the present study, classic linkage mapping combined with QTL-seq were employed in two recombinant inbred line (RIL) populations, "Montrose"/I9365-25 (M25) and "Raven"/I9365-31 (R31), with the initial goal of fine-mapping QTL WM5.4 and WM7.5 that condition WM resistance. The RILs were phenotyped for WM reactions under greenhouse (straw test) and field environments. The general region of WM5.4 and WM7.5 were reconfirmed with both mapping strategies within each population. Combining the results from both mapping strategies, WM5.4 was delimited to a 22.60-36.25 Mb interval in the heterochromatic regions on Pv05, while WM7.5 was narrowed to a 0.83 Mb (3.99-4.82 Mb) region on the Pv07 chromosome. Furthermore, additional QTL WM2.2a (3.81-7.24 Mb), WM2.2b (11.18-17.37 Mb, heterochromatic region), and WM2.2c (23.33-25.94 Mb) were mapped to a narrowed genomic interval on Pv02 and WM4.2 in a 0.89 Mb physical interval at the distal end of Pv04 chromosome. Gene models encoding gibberellin 2-oxidase proteins regulating plant architecture are likely candidate genes associated with WM2.2a resistance. Nine gene models encoding a disease resistance protein (quinone reductase family protein and ATWRKY69) found within the WM5.4 QTL interval are putative candidate genes. Clusters of 13 and 5 copies of gene models encoding cysteine-rich receptor-like kinase and receptor-like protein kinase-related family proteins, respectively, are potential candidate genes associated with WM7.5 resistance and most likely trigger physiological resistance to WM. Acquired knowledge of the narrowed major QTL intervals, flanking markers, and candidate genes provides promising opportunities to develop functional molecular markers to implement marker-assisted selection for WM resistant dry bean cultivars.

Whole-genome QTL analysis of Stagonospora nodorum blotch resistance and validation of the SnTox4-Snn4 interaction in hexaploid wheat.

Necrotrophic effectors (also known as host-selective toxins) are important determinants of disease in the wheat-Stagonospora nodorum pathosystem. To date, five necrotrophic effector-host gene interactions have been identified in this system. Most of these interactions have additive effects while some are epistatic. The Snn4-SnTox4 interaction was originally identified in a recombinant-inbred population derived from a cross between the Swiss winter wheat cultivars 'Arina' and 'Forno' using the S. nodorum isolate Sn99CH 1A7a. Here, we used a recombinant-inbred population consisting of 121 lines developed from a cross between the hexaploid land race Salamouni and the hexaploid wheat 'Katepwa' (SK population). The SK population was used for the construction of linkage maps and quantitative trait loci (QTL) detection using the Swiss S. nodorum isolate Sn99CH 1A7a. The linkage maps developed in the SK population spanned 3,228 centimorgans (cM) and consisted of 441 simple-sequence repeats, 9 restriction fragment length polymorphisms, 29 expressed sequence tag sequence-tagged site markers, and 5 phenotypic markers. The average marker density was 6.7 cM/marker. Two QTL, designated QSnb.fcu-1A and QSnb.fcu-7A on chromosome arms 1AS and 7AS, respectively, were associated with disease caused by the S. nodorum isolate Sn99CH 1A7a. The effects of QSnb.fcu-1A were determined by the Snn4-SnTox4 interaction and accounted for 23.5% of the phenotypic variation in this population, whereas QSnb.fcu-7A accounted for 16.4% of the phenotypic variation for disease but was not associated with any known effector sensitivity locus. The effects of both QTL were largely additive and collectively accounted for 35.7% of the total phenotypic variation. The results of this research validate the effects of a compatible Snn4-SnTox4 interaction in a different genetic background, and it provides knowledge regarding genomic regions and molecular markers that can be used to improve Stagonospora nodorum blotch resistance in wheat germplasm.

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.

Modified screening method of middle american dry bean genotypes reveals new genomic regions on Pv10 associated with anthracnose resistance.

Anthracnose, caused by the fungal pathogen Colletotrichum lindemuthianum (Sacc. & Magnus) Lams.-Scrib., is one of the most devastating diseases in dry bean (Phaseolus vulgaris L.) with seed yield losses up to 100%. Most anthracnose resistance genes thus far identified behave in a dominant manner and were identified by seedling screening. The Middle American Diversity Panel (MDP; n=266) was screened with a modified greenhouse screening method to evaluate the response to anthracnose race 73. Thirty MDP genotypes exhibited resistance to the race of which 16 genotypes were not known to contain anthracnose resistance genes to race 73. GWAS with ~93,000 SNP markers identified four genomic regions, two each on Pv01 and Pv10, associated race 73 resistance. A likelihood-ratio-based R2 analysis indicated the peak four SNP markers are responsible for 26% of the observed phenotypic variation, where one SNP, S10_072250, explains 23% of the total variation. SNP S10_072250 is associated with a new region of anthracnose resistance and is in an intron of a ZPR1-like gene. Further greenhouse testing of the 16 resistant lines without previously known resistance to race 73 revealed various levels of resistance under various levels of disease pressure. Disease resistance was further characterized in the field using four representative genotypes. GTS-900 and Remington exhibited field resistance while Merlot and Maverick were susceptible. Field testing with two different fungicide regimes revealed the resistant genotypes had no significant disease differences. The results suggest resistance to anthracnose may differ at various growth stages and that breeders have been selecting for major genes at early seedling stages while ignoring the effect of alternative genes that may be active at later stages. The newly identified resistant lines may be related to Age Related Resistance (ARR) and could be exploited as parental sources of anthracnose resistance in addition to already known major genes. The physical localization of the multiple regions of resistance confirms the presence of two clusters of disease resistance genes on Pv01 and identifies two new regions of anthracnose resistance on Pv10 possibly associated with ARR. Future research should look at the mode of inheritance of this resistance and its effect when combined with other anthracnose resistance loci.