A Predictive Model for Time-to-Flowering in the Common Bean Based on QTL and Environmental Variables

Identification of superior haplotypes for flowering time in pigeonpea through candidate gene-based association study of a diverse minicore collection

KEY MESSAGE

In current study candidate gene (261 genes) based association mapping on 144 pigeonpea accessions for flowering time and related traits and 29 MTAs producing eight superior haplotypes were identified. In the current study, we have conducted an association analysis for flowering-associated traits in a diverse pigeonpea mini-core collection comprising 144 accessions using the SNP data of 261 flowering-related genes. In total, 13,449 SNPs were detected in the current study, which ranged from 743 (ICP10228) to 1469 (ICP6668) among the individuals. The nucleotide diversity (0.28) and Watterson estimates (0.34) reflected substantial diversity, while Tajima's D (-0.70) indicated the abundance of rare alleles in the collection. A total of 29 marker trait associations (MTAs) were identified, among which 19 were unique to days to first flowering (DOF) and/or days to fifty percent flowering (DFF), 9 to plant height (PH), and 1 to determinate (Det) growth habit using 3 years of phenotypic data. Among these MTAs, six were common to DOF and/or DFF, and four were common to DOF/DFF along with the PH, reflecting their pleiotropic action. These 29 MTAs spanned 25 genes, among which 10 genes clustered in the protein-protein network analysis, indicating their concerted involvement in floral induction. Furthermore, we identified eight haplotypes, four of which regulate late flowering, while the remaining four regulate early flowering using the MTAs. Interestingly, haplotypes conferring late flowering (H001, H002, and H008) were found to be taller, while those involved in early flowering (H003) were shorter in height. The expression pattern of these genes, as inferred from the transcriptome data, also underpinned their involvement in floral induction. The haplotypes identified will be highly useful to the pigeonpea breeding community for haplotype-based breeding.

Plant mitochondrial introns as genetic markers - conservation and variation

Plant genomes are comprised of nuclear, plastid and mitochondrial components characterized by different patterns of inheritance and evolution. Genetic markers from the three genomes provide complementary tools for investigations of inheritance, genetic relationships and phenotypic contributions. Plant mitochondrial genomes are challenging for universal marker development because they are highly variable in terms of size, gene order and intergenic sequences and highly conserved with respect to protein-coding sequences. PCR amplification of introns with primers that anneal to conserved, flanking exons is effective for the development of polymorphic nuclear genome markers. The potential for plant mitochondrial intron polymorphisms to distinguish between congeneric species or intraspecific varieties has not been systematically investigated and is possibly constrained by requirements for intron secondary structure and interactions with co-evolved organelle intron splicing factors. To explore the potential for broadly applicable plant mitochondrial intron markers, PCR primer sets based upon conserved sequences flanking 11 introns common to seven angiosperm species were tested across a range of plant orders. PCR-amplified introns were screened for indel polymorphisms among a group of cross-compatible Citrus species and relatives; two Raphanus sativus mitotypes; representatives of the two Phaseolus vulgaris gene pools; and congeneric pairs of Cynodon, Cenchrus, Solanum, and Vaccinium species. All introns were successfully amplified from each plant entry. Length polymorphisms distinguishable by gel electrophoresis were common among genera but infrequent within genera. Sequencing of three introns amplified from 16 entries identified additional short indel polymorphisms and nucleotide substitutions that separated Citrus, Cynodon, Cenchrus and Vaccinium congeners, but failed to distinguish Solanum congeners or representatives of the Phaseolus vulgaris major gene pools. The ability of primer sets to amplify a wider range of plant species' introns and the presence of intron polymorphisms that distinguish congeners was confirmed by in silico analysis. While mitochondrial intron variation is limited in comparison to nuclear introns, these exon-based primer sets provide robust tools for the amplification of mitochondrial introns across a wide range of plant species wherein useful polymorphisms can be identified.

The Domestication of the Seasons: The Exploitation of Variations in Crop Seasonality Responses by Later Prehistoric Farmers

Discussions of domestication and its associated genetic traits has focused upon dispersal mechanisms, for example in the rachis structure of cereals. Here we consider another trait of importance to domesticated crops – their seasonality response. We particularly consider flowering response to changing daylength, and how that relates to domesticated species that are carried by farmers considerable distances from their regions of origin. We consider the implications of cultivating crops in regions whose seasonality patterns contrast with those of the crop’s region of origin, and consider the consequent genetic changes. Reference is made to the particular case of barley, and a discussion of archaeogenetic research into its photoperiod response mechanism is provided.

Dynamic QTL-based ecophysiological models to predict phenotype from genotype and environment data

BACKGROUND

Predicting the phenotype from the genotype is one of the major contemporary challenges in biology. This challenge is greater in plants because their development occurs mostly post-embryonically under diurnal and seasonal environmental fluctuations. Most current crop simulation models are physiology-based models capable of capturing environmental fluctuations but cannot adequately capture genotypic effects because they were not constructed within a genetics framework.

RESULTS

We describe the construction of a mixed-effects dynamic model to predict time-to-flowering in the common bean (Phaseolus vulgaris L.). This prediction model applies the developmental approach used by traditional crop simulation models, uses direct observational data, and captures the Genotype, Environment, and Genotype-by-Environment effects to predict progress towards time-to-flowering in real time. Comparisons to a traditional crop simulation model and to a previously developed static model shows the advantages of the new dynamic model.

CONCLUSIONS

The dynamic model can be applied to other species and to different plant processes. These types of models can, in modular form, gradually replace plant processes in existing crop models as has been implemented in BeanGro, a crop simulation model within the DSSAT Cropping Systems Model. Gene-based dynamic models can accelerate precision breeding of diverse crop species, particularly with the prospects of climate change. Finally, a gene-based simulation model can assist policy decision makers in matters pertaining to prediction of food supplies.

Improving Association Studies and Genomic Predictions for Climbing Beans With Data From Bush Bean Populations

Common bean (Phaseolus vulgaris L.) has two major origins of domestication, Andean and Mesoamerican, which contribute to the high diversity of growth type, pod and seed characteristics. The climbing growth habit is associated with increased days to flowering (DF), seed iron concentration (SdFe), nitrogen fixation, and yield. However, breeding efforts in climbing beans have been limited and independent from bush type beans. To advance climbing bean breeding, we carried out genome-wide association studies and genomic predictions using 1,869 common bean lines belonging to five breeding panels representing both gene pools and all growth types. The phenotypic data were collected from 17 field trials and were complemented with 16 previously published trials. Overall, 38 significant marker-trait associations were identified for growth habit, 14 for DF, 13 for 100 seed weight, three for SdFe, and one for yield. Except for DF, the results suggest a common genetic basis for traits across all panels and growth types. Seven QTL associated with growth habits were confirmed from earlier studies and four plausible candidate genes for SdFe and 100 seed weight were newly identified. Furthermore, the genomic prediction accuracy for SdFe and yield in climbing beans improved up to 8.8% when bush-type bean lines were included in the training population. In conclusion, a large population from different gene pools and growth types across multiple breeding panels increased the power of genomic analyses and provides a solid and diverse germplasm base for genetic improvement of common bean.

Genetic dissection of days to flowering via genome-wide association studies in Turkish common bean germplasm

Common bean is a nutrient-dense legume crop serving as a source of food for millions of people. Characterization of unexplored common bean germplasm to unlock the phenotypic and genetic variations is still needed to explore the breeding potential of this crop. The current study aimed to dissect the genetic basis having association for days to flowering (DF). A total of 188 common bean accessions collected from 19 provinces of Turkey were used as plant material under five environments and two locations. Analysis of variance (ANOVA) revealed that genotypes and genotype by environment interaction have significant effects on DF. A total of 10 most stable accessions were evaluated from stability analysis. Overall maximum (75) and minimum (54) DF were observed for Hakkari-51 and Mus-46 accessions, respectively. The implemented constellation plot divided studied germplasm according to their DF and growth habit. A total of 7900 DArTseq markers were used for association analysis. Mixed linear model using the Q + K Model resulted a total of 18 DArTseq markers from five environments. DArT-8668385 marker identified in Bolu during 2016 was also associated with DF in Sivas during 2017. Combined data of five years resulted a total of four markers (DArT-22346534, DArT-3369768, DArT-3374613, and DArT-3370801) having significant association ( p  <  0.01 ) for DF. DArT-22346534 present on Pv 08 accounted a maximum of 9.89% variation to the studied trait. A total of four putative candidate genes were predicted from sequences reflecting homology to identified four DArTseq markers. We envisage that exploitation of identified DArTseq markers will hopefully beneficial for the development of new common bean varieties having better adaptation ability to changing climatic conditions.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01029-8.

Ancient relaxation of an obligate short-day requirement in common bean through loss of CONSTANS-like gene function

Common bean (Phaseolus vulgaris L.) is a major global food staple and source of dietary protein that was domesticated independently in Mexico and Andean South America. Its subsequent development as a crop of importance worldwide has been enabled by genetic relaxation of the strict short-day requirement typical of wild forms, but the genetic basis for this change is not well understood. Recently, a loss of photoperiod sensitivity was shown to result from mutations in the phytochrome photoreceptor gene Ppd/PHYA3 that arose independently within the two major domesticated lineages. Here, we define a second major photoperiod sensitivity locus, at which recessive alleles associate with deleterious mutations affecting the CONSTANS-like gene COL2. A wider survey of sequence variation in over 800 diverse lines, including wild, landrace, and domesticated accessions, show that distinct col2 haplotypes are associated with early flowering in Andean and Mesoamerican germplasm. The relative frequencies and distributions of COL2 and PHYA3 haplotypes imply that photoperiod adaptation developed in two phases within each gene pool: an initial reduction in sensitivity through impairment of COL2 function and subsequent complete loss through PHYA3. Gene expression analyses indicate that COL2 functions downstream of PHYA3 to repress expression of FT genes and may function in parallel with PvE1, the bean ortholog of a key legume-specific flowering repressor. Collectively, these results define the molecular basis for a key phenological adaptation, reveal a striking convergence in the naturally replicated evolution of this major crop, and further emphasize the wider evolutionary lability of CONSTANS effects on flowering time control.

Genome-wide association studies of Ca and Mn in the seeds of the common bean (Phaseolus vulgaris L.)

SNP markers linked to genes controlling Ca and Mn uptake were identified in the common bean seeds using DArT-based association mapping (AM). The Ca concentration in the seeds varied between 475 and 3,100 mg kg^-1 with an average of 1,280.9 mg kg^-1 and the Mn concentration ranged from 4.87 to 27.54 mg kg^-1 with a mean of 11.76 mg kg^-1. A total of 19,204 SNP markers were distributed across 11 chromosomes that correspond to the haploid genome number of the common bean. The highest value of ΔK was determined as K = 2, and 173 common bean genotypes were split into two main subclusters as POP1 (Mesoamerican) and POP2 (Andean). The results of the UPGMA dendrogram and PCA confirmed those of STRUCTURE analysis. MLM based on the Q + K model identified a large number of markers-trait associations. Of the 19,204 SNPs, five (on Pv2, 3, 8, 10 and 11) and four (on Pv2, 3, 8 and 11) SNPs were detected to be significantly related to the Ca content of the beans grown in Bornova and Menemen, respectively in 2015. In 2016, six SNPs (on Pv1-4, 8 and 10) were identified to be significantly associated with the Ca content of the seeds obtained from Bornova and six SNPs (on Pv1-4, 8 and 10) from Menemen. Eight (on Pv3, 5 and 11) and four (on Pv2, 5 and 11) SNPs had a significant association with Mn content in Bornova in 2015 and 2016, respectively. In Menemen, eight (on Pv3, 5, 8 and 11) and 11 (on Pv1, 2, 5, 10 and 11) SNPs had a significant correlation with Mn content in 2015 and 2016, respectively.

Genomic Prediction of Agronomic Traits in Common Bean (Phaseolus vulgaris L.) Under Environmental Stress

In plant and animal breeding, genomic prediction models are established to select new lines based on genomic data, without the need for laborious phenotyping. Prediction models can be trained on recent or historic phenotypic data and increasingly available genotypic data. This enables the adoption of genomic selection also in under-used legume crops such as common bean. Beans are an important staple food in the tropics and mainly grown by smallholders under limiting environmental conditions such as drought or low soil fertility. Therefore, genotype-by-environment interactions (G × E) are an important consideration when developing new bean varieties. However, G × E are often not considered in genomic prediction models nor are these models implemented in current bean breeding programs. Here we show the prediction abilities of four agronomic traits in common bean under various environmental stresses based on twelve field trials. The dataset includes 481 elite breeding lines characterized by 5,820 SNP markers. Prediction abilities over all twelve trials ranged between 0.6 and 0.8 for yield and days to maturity, respectively, predicting new lines into new seasons. In all four evaluated traits, the prediction abilities reached about 50-80% of the maximum accuracies given by phenotypic correlations and heritability. Predictions under drought and low phosphorus stress were up to 10 and 20% improved when G × E were included in the model, respectively. Our results demonstrate the potential of genomic selection to increase the genetic gain in common bean breeding. Prediction abilities improved when more phenotypic data was available and G × E could be accounted for. Furthermore, the developed models allowed us to predict genotypic performance under different environmental stresses. This will be a key factor in the development of common bean varieties adapted to future challenging conditions.

The genetic control of leaf allometry in the common bean, Phaseolus vulgaris

BACKGROUND

To maximize photosynthetic efficiency, plants have evolved a capacity by which leaf area scales allometrically with leaf mass through interactions with the environment. However, our understanding of genetic control of this allometric relationship remains limited.

RESULTS

We integrated allometric scaling laws expressed at static and ontogenetic levels into genetic mapping to identify the quantitative trait loci (QTLs) that mediate how leaf area scales with leaf mass and how such leaf allometry, under the control of these QTLs, varies as a response to environment change. A major QTL detected by the static model constantly affects the allometric growth of leaf area vs. leaf mass for the common bean (Phaseolus vulgaris) in two different environments. The ontogenetic model identified this QTL plus a few other QTLs that determine developmental trajectories of leaf allometry, whose expression is contingent heavily upon the environment.

CONCLUSIONS

Our results gain new insight into the genetic mechanisms of how plants program their leaf morphogenesis to adapt to environmental perturbations.

Characterization of QTL and Environmental Interactions Controlling Flowering Time in Andean Common Bean (Phaseolus vulgaris L.)

Genetic variation for response of flowering time to photoperiod plays an important role in adaptation to environments with different photoperiods, and as consequence is an important contributor to plant productivity and yield. To elucidate the genetic control of flowering time [days to flowering (DTF); growing degree days (GDD)] in common bean, a facultative short-day plant, a quantitative trait loci (QTL) analysis was performed in a recombinant inbred mapping population derived from a cultivated accession and a photoperiod sensitive landrace, grown in different long-day (LD) and short-day (SD) environments by using a multiple-environment QTL model approach. A total of 37 QTL across 17 chromosome regions and 36 QTL-by-QTL interactions were identified for six traits associated with time to flowering and response to photoperiod. The DTF QTL accounted for 28 and 11% on average of the phenotypic variation in the population across LD and SD environments, respectively. Of these, a genomic region on chromosome 4 harboring the major DTF QTL was associated with both flowering time in LD and photoperiod response traits, controlling more than 60% of phenotypic variance, whereas a major QTL on chromosome 9 explained up to 32% of flowering time phenotypic variation in SD. Different epistatic interactions were found in LD and SD environments, and the presence of significant QTL × environment (QE) and epistasis × environment interactions implies that flowering time control may rely on different genes and genetic pathways under inductive and non-inductive conditions. Here, we report the identification of a novel major locus controlling photoperiod sensitivity on chromosome 4, which might interact with other loci for controlling common bean flowering time and photoperiod response. Our results have also demonstrated the importance of these interactions for flowering time control in common bean, and point to the likely complexity of flowering time pathways. This knowledge will help to identify and develop opportunities for adaptation and breeding of this legume crop.

Parallel origins of photoperiod adaptation following dual domestications of common bean

Common bean (Phaseolus vulgaris L.) is an important grain legume domesticated independently in Mexico and Andean South America approximately 8000 years ago. Wild forms are obligate short-day plants, and relaxation of photoperiod sensitivity was important for expansion to higher latitudes and subsequent global spread. To better understand the nature and origin of this key adaptation, we examined its genetic control in progeny of a wide cross between a wild accession and a photoperiod-insensitive cultivar. We found that photoperiod sensitivity is under oligogenic control, and confirm a major effect of the Ppd locus on chromosome 1. The red/far-red photoreceptor gene PHYTOCHROME A3 (PHYA3) was identified as a strong positional candidate for Ppd, and sequencing revealed distinct deleterious PHYA3 mutations in photoperiod-insensitive Andean and Mesoamerican accessions. These results reveal the independent origins of photoperiod insensitivity within the two major common bean gene pools and demonstrate the conserved importance of PHYA genes in photoperiod adaptation of short-day legume species.