Last-Generation Genome-Environment Associations Reveal the Genetic Basis of Heat Tolerance in Common Bean (Phaseolus vulgaris L.)

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.

Genomic evidence for climate-linked diversity loss and increased vulnerability of wild barley spanning 28 years of climate warming

The information on how plant populations respond genetically to climate warming is scarce. Here, landscape genomic and machine learning approaches were integrated to assess genetic response of 10 wild barley (Hordeum vulgare ssp. spontaneum; WB) populations in the past and future, using whole genomic sequencing (WGS) data. The WB populations were sampled in 1980 and again in 2008. Phylogeny of accessions was roughly in conformity with sampling sites, which accompanied by admixture/introgressions. The 28-y climate warming resulted in decreased genetic diversity, increased selection pressure, and an increase in deleterious single nucleotide polymorphism (dSNP) numbers, heterozygous deleterious and total deleterious burdens for WB. Genome-environment associations identified some candidate genes belonging to peroxidase family (HORVU2Hr1G057450, HORVU4Hr1G052060 and HORVU4Hr1G057210) and heat shock protein 70 family (HORVU2Hr1G112630). The gene HORVU2Hr1G120170 identified by selective sweep analysis was under strong selection during the climate warming of the 28-y, and its derived haplotypes were fixed by WB when faced with the 28-y increasingly severe environment. Temperature variables were found to be more important than precipitation variables in influencing genomic variation, with an eco-physiological index gdd5 (growing degree-days at the baseline threshold temperature of 5 °C) being the most important determinant. Gradient forest modelling revealed higher predicted genomic vulnerability in Sede Boqer under future climate scenarios at 2041-2070 and 2071-2100. Additionally, estimates of effective population size (Ne) tracing back to 250 years indicated a forward decline in all populations over time. Our assessment about past genetic response and future vulnerability of WB under climate warming is crucial for informing conservation efforts for wild cereals and rational use strategies.

Genome-Wide Association Study Reveals a QTL Region for Ashy Stem Blight Resistance in PRA154 Andean Common Bean

Ashy stem blight (ASB) caused by Macrophomina phaseolina (Tassi) Goidanich affects the common bean (Phaseolus vulgaris L.) at all growing stages. Higher levels of resistance were observed in Andean common beans, but specific resistant quantitative trait loci (QTLs) conferring resistance to this pathogen have not been reported in this gene pool. The objectives of this research were to: (i) conduct a genome-wide association study (GWAS) and QTL mapping for resistance in the Andean breeding line PRA154; and (ii) identify single nucleotide polymorphism (SNP) markers and candidate genes for ASB resistance. Phenotyping was conducted under greenhouse conditions by inoculating the 107 F6:7 recombinant inbred lines (RILs) derived from the cross between the susceptible cultivar 'Verano' and the partial-resistant breeding line PRA154 twice with the M. phaseolina isolate PRI21. Genotyping was performed with 109,040 SNPs distributed across all 11 P. vulgaris chromosomes. A novel major QTL was located between 28,761,668 and 31,263,845 bp, extending 2.5 Mbp on chromosome Pv07, and the highest significant SNP markers were Chr07_28761668_A_G, Chr07_29131720_G_A, and Chr07_31263845_C_T with the highest LOD (more than 10 in most of the cases) and R-squared values, explaining 40% of the phenotypic variance of the PRI21 isolate. The gene Phvul.007G173900 (methylcrotonyl-CoA carboxylase alpha chain and mitochondrial 3-methylcrotonyl-CoA carboxylase 1 [MCCA]) with a size of 10,891 bp, located between 29,131,591 and 29,142,481 bp on Pv07, was identified as one candidate for ASB resistance in PRA154, and it contained Chr07_29131720_G_A. The QTL and genetic marker information could be used to assist common bean breeders to develop germplasm and cultivars with ASB resistance through molecular breeding.

High Genetic Diversity in the Himalayan Common Bean (Phaseolus vulgaris) Germplasm with Divergence from Its Center of Origin in the Mesoamerica and Andes

The common bean is found in the Himalayan region of Pakistan with substantial morphological variability. Genetic diversity within any crop species is a precursor for genetic improvement; however, little is known about common bean genetic diversity in this region. We explored the genetic diversity in the common bean from the Himalayan region (Khyber Pakhtunkhwa, Gilgit-Baltistan, Kashmir) of Pakistan. Microsatellite genotyping was carried out for 147 samples with 40 simple sequence repeat (SSR) markers. The results revealed a clear divergence of the Pakistani population from the primary gene pool (with FST values of 0.2 with Andes and 0.27 with Mesoamerica). However, within the Himalayan germplasm, no clear evidence of spatial structure was observed (with the maximum FST values of only 0.025), probably due to the dispersal of seeds by human activity within the region. This was further elucidated by the discriminant analyses of principal components. Considering the diversity parameters, high genotypic diversity was observed for the indigenous lines (0.990), comparable to the primary gene pool (0.976 for Mesoamerica and 0.976 for Andes populations). A high genotypic diversity was observed within the Himalayan population (ranging from 0.500 for Upper Dir to 0.952 for Mansehra). Gene diversity across loci varied between 0.28 for Chitral to 0.38 for Kurram. Our results suggested a divergent and independent evolution of the Himalayan population, which might have led to the diversification of the common bean germplasm in the region postintroduction into the region. The diversity observed could also be exploited in future breeding programs for the development and introduction of climate-resilient varieties.

Comprehensive Time-Course Transcriptome Reveals the Crucial Biological Pathways Involved in the Seasonal Branch Growth in Siberian Elm (Ulmus pumila)

Timber, the most prevalent organic material on this planet, is the result of a secondary xylem emerging from vascular cambium. Yet, the intricate processes governing its seasonal generation are largely a mystery. To better understand the cyclic growth of vascular tissues in elm, we undertook an extensive study examining the anatomy, physiology, and genetic expressions in Ulmus pumila. We chose three robust 15-year-old elm trees for our study. The cultivars used in this study were collected from the Inner Mongolia Autonomous Region in China and nurtured in the tree farm of Shandong Normal University. Monthly samples of 2-year-old elm branches were taken from the tree from February to September. Marked seasonal shifts in elm branch vascular tissues were observed by phenotypic observation: In February, the cambium of the branch emerged from dormancy, spurring growth. By May, elms began generating secondary xylem, or latewood, recognized by its tiny pores and dense cell structure. From June to August, there was a marked increase in the thickness of the secondary xylem. Transcriptome sequencing provides a potential molecular mechanism for the thickening of elm branches and their response to stress. In February, the tree enhanced its genetic responses to cold and drought stress. The amplified expression of CDKB, CYCB, WOX4, and ARF5 in the months of February and March reinforced their essential role in the development of the vascular cambium in elm. Starting in May, the elm deployed carbohydrates as a carbon resource to synthesize the abundant cellulose and lignin necessary for the formation of the secondary wall. Major genes participating in cellulose (SUC and CESA homologs), xylan (UGD, UXS, IRX9, IRX10, and IRX14), and lignin (PAL, C4H, 4CL, HCT, C3H, COMT, and CAD) biosynthetic pathways for secondary wall formation were up-regulated by May or/and June. In conclusion, our findings provided a foundation for an in-depth exploration of the molecular processes dictating the seasonal growth of elm timber.

Identification of Two Flip-Over Genes in Grass Family as Potential Signature of C4 Photosynthesis Evolution

In flowering plants, C4 photosynthesis is superior to C3 type in carbon fixation efficiency and adaptation to extreme environmental conditions, but the mechanisms behind the assembly of C4 machinery remain elusive. This study attempts to dissect the evolutionary divergence from C3 to C4 photosynthesis in five photosynthetic model plants from the grass family, using a combined comparative transcriptomics and deep learning technology. By examining and comparing gene expression levels in bundle sheath and mesophyll cells of five model plants, we identified 16 differentially expressed signature genes showing cell-specific expression patterns in C3 and C4 plants. Among them, two showed distinctively opposite cell-specific expression patterns in C3 vs. C4 plants (named as FOGs). The in silico physicochemical analysis of the two FOGs illustrated that C3 homologous proteins of LHCA6 had low and stable pI values of ~6, while the pI values of LHCA6 homologs increased drastically in C4 plants Setaria viridis (7), Zea mays (8), and Sorghum bicolor (over 9), suggesting this protein may have different functions in C3 and C4 plants. Interestingly, based on pairwise protein sequence/structure similarities between each homologous FOG protein, one FOG PGRL1A showed local inconsistency between sequence similarity and structure similarity. To find more examples of the evolutionary characteristics of FOG proteins, we investigated the protein sequence/structure similarities of other FOGs (transcription factors) and found that FOG proteins have diversified incompatibility between sequence and structure similarities during grass family evolution. This raised an interesting question as to whether the sequence similarity is related to structure similarity during C4 photosynthesis evolution.

Identification of QTLs and Putative Candidate Genes for Plant Architecture of Lotus Revealed by Regional Association Mapping

The lotus (Nelumbo Adans.) is one of the most economically relevant ornamental aquatic plants. Plant architecture (PA) is an important trait for lotus classification, cultivation, breeding, and applications. However, the underlying genetic and molecular basis controlling PA remains poorly understood. In this study, an association study for PA-related traits was performed with 93 genome-wide microsatellite markers (simple sequence repeat, SSR) and 51 insertion-deletion (InDel) markers derived from the candidate regions using a panel of 293 lotus accessions. Phenotypic data analysis of the five PA-related traits revealed a wide normal distribution and high heritability from 2013 to 2016, which indicated that lotus PA-related traits are highly polygenic traits. The population structure (Q-matrix) and the relative kinships (K-matrix) of the association panels were analyzed using 93 SSR markers. The mixed linear model (MLM) taking Q-matrix and K-matrix into account was used to estimate the association between markers and the traits. A total of 26 markers and 65 marker-trait associations were identified by considering associations with p < 0.001 and Q < 0.05. Based on the significant markers, two QTLs on Chromosome 1 were identified, and two candidate genes were preliminarily determined. The results of our study provided useful information for the lotus breeding aiming at different PA phenotypes using a molecular-assisted selection (MAS) method and also laid the foundation for the illustration of the molecular mechanism underlying the major QTL and key markers associated with lotus PA.

The genome size, chromosome number and the seed adaption to long-distance dispersal of Ipomoea pes-caprae (L.)

Ipomoeapes-caprae (L.) (IPC) is a common species in tropical and subtropical coastal areas and one of the world's most widely distributed plants. It has attracted researchers for its outstanding biological, ecological and medicinal values. It has been reported that the genetic diversity of IPCs located on different continents is very low because of their frequent gene flow. During the long journey of evolution, every aspect of the plant morphologies has evolved to the best adaptivity to the environment, seeking their survival and progeny expansion. However, the fundamental genetic characteristics of IPC and how their seed adapted to the success of population expansion remain unknown. In this study, the fundamental genetic characteristics, including the genome size and the chromosome number of IPC, were investigated. The results showed that IPC's genome size is approximately 0.98-1.08 GB, and the chromosome number is 2n=30, providing the basic information for further genome analysis. In order to decipher the long-distance dispersal secret of this species, the fruit and seed developments, seed morphology, and seed germination were extensively investigated and described. The results showed an exquisite adaptive mechanism of IPC seeds to fulfil the population expansion via ocean currents. The large cavity inside the seeds and the dense tomenta on the surface provide the buoyancy force for the seeds to float on the seawater. The hard seed coats significantly obstructed the water absorption, thus preventing the seed from germination during the dispersal. Meanwhile, the fully developed embryos of IPC also have physiological dormancy. The physical and physiological characteristics of IPC seeds provide insight into the mechanism of their long-distance dispersal across the oceans. Moreover, based on morphological observation and semi-section microscopy, the development pattern of IPC glander trichomes was described, and their physiological functions were also discussed.

Delineating the role of plant stature towards heat stress tolerance in field pea (Pisum sativum L.)

Terminal heat stress severely affects field pea production in tropical climates. Identifying and characterizing marker-trait(s) remain vital for breeding heat-tolerant cultivars of field pea. Field pea genotypes are highly variable for plant stature; however, the significance of plant stature for yield stability under high-temperature conditions is not yet well understood. The study aimed to investigate the sensitivity and significance of plant stature toward yield sustainability of field pea under high-temperature environments. A panel of 150 diverse genotypes with variable plant statures [dwarf (<50 cm), semi-dwarf (50-80 cm), medium-tall (80-150 cm)] were grown under late sowing-induced high-temperature environments for two consecutive years (2017-2019). During the first year of the experiment, the late sown crops (15 and 30 days) were exposed to high-temperatures at flowering (+3.5 to +8.1 °C) and grain-filling (+3.3 to +6.1 °C) over timely sown crops. Likewise, elevated temperature during flowering (+3.7 to +5.2 °C) and grain filling (+5.4 to +9.9 °C) were recorded in late-sown environments (delayed by 27 and 54 days) in the next year. Medium-tall genotypes had longer grain-filling duration (7-10%), higher pod-bearing nodes (8-18%) and yield (22-55%), and lower yield losses (13-18%) over semi-dwarf and dwarf genotypes under high-temperature environments. Significant associations of plant height with yield, yield loss, and heat-susceptibility index in high-temperature environments suggested higher heat tolerance capacity of tall-type plants compared to dwarf and semi-dwarf types. GGEbiplot analysis revealed that the heat-tolerant genotypes were all medium tall-type (mean = 108 cm), while the heat-susceptible genotypes were mostly dwarf in stature. Hence, tall-type genotypes had better adaptability to high-temperature environments. Henceforth, the breeding approach for high-temperature tolerance in field pea may be designed by embracing tall-type backgrounds over dwarf plant to develop climate resilient cultivars.

Evolution and Expression Patterns of the Fructose 1,6-Bisphosptase Gene Family in a Miracle Tree (Neolamarckia cadamba)

Neolamarckia cadamba (N. cadamba) is a fast-growing tree species with tremendous economic and ecological value; the study of the key genes regulating photosynthesis and sugar accumulation is very important for the breeding of N. cadamba. Fructose 1,6-bisphosptase (FBP) gene has been found to play a key role in plant photosynthesis, sugar accumulation and other growth processes. However, no systemic analysis of FBPs has been reported in N. cadamba. A total of six FBP genes were identifed and cloned based on the N. cadamba genome, and these FBP genes were sorted into four groups. The characteristics of the NcFBP gene family were analyzed such as phylogenetic relationships, gene structures, conserved motifs, and expression patterns. A cis-acting element related to circadian control was first found in the promoter region of FBP gene. Phylogenetic and quantitative real-time PCR analyses showed that NcFBP5 and NcFBP6 may be chloroplast type 1 FBP and cytoplasmic FBP, respectively. FBP proteins from N. cadamba and 22 other plant species were used for phylogenetic analyses, indicating that FBP family may have expanded during the evolution of algae to mosses and differentiated cpFBPase1 proteins in mosses. This work analyzes the internal relationship between the evolution and expression of the six NcFBPs, providing a scientific basis for the evolutionary pattern of plant FBPs, and promoting the functional studies of FBP genes.

Improving Coconut Using Modern Breeding Technologies: Challenges and Opportunities

Coconut (Cocos nucifera L.) is a perennial palm with a wide range of distribution across tropical islands and coastlines. Multitude use of coconut by nature is important in the socio-economic fabric framework among rural smallholders in producing countries. It is a major source of income for 30 million farmers, while 60 million households rely on the coconut industry directly as farm workers and indirectly through the distribution, marketing, and processing of coconut and coconut-based products. Stagnant production, inadequate planting materials, the effects of climate change, as well as pests and diseases are among the key issues that need to be urgently addressed in the global coconut industry. Biotechnology has revolutionized conventional breeding approaches in creating genetic variation for trait improvement in a shorter period of time. In this review, we highlighted the challenges of current breeding strategies and the potential of biotechnological approaches, such as genomic-assisted breeding, next-generation sequencing (NGS)-based genotyping and genome editing tools in improving the coconut. Also, combining these technologies with high-throughput phenotyping approaches and speed breeding could speed up the rate of genetic gain in coconut breeding to solve problems that have been plaguing the industry for decades.

Biochemical and Transcriptional Responses in Cold-Acclimated and Non-Acclimated Contrasting Camelina Biotypes under Freezing Stress

Cold-acclimated and non-acclimated contrasting Camelina (Camelina sativa L.) biotypes were investigated for changes in stress-associated biomarkers, including antioxidant enzyme activity, lipid peroxidation, protein, and proline content. In addition, a well-known freezing tolerance pathway participant known as C-repeat/DRE-binding factors (CBFs), an inducer of CBF expression (ICE1), and a cold-regulated (COR6.6) genes of the ICE-CBF-COR pathway were studied at the transcriptional level on the doubled-haploid (DH) lines. Freezing stress had significant effects on all studied parameters. The cold-acclimated DH34 (a freezing-tolerant line) showed an overall better performance under freezing stress than non-acclimated plants. The non-cold-acclimated DH08 (a frost-sensitive line) showed the highest electrolyte leakage after freezing stress. The highest activity of antioxidant enzymes (glutathione peroxidase, superoxide dismutase, and catalase) was also detected in non-acclimated plants, whereas the cold-acclimated plants showed lower enzyme activities upon stress treatment. Cold acclimation had a significantly positive effect on the total protein and proline content of stressed plants. The qRT-PCR analysis revealed significant differences in the expression and cold-inducibility of CsCBF1-3, CsICE1, and CsCOR6.6 genes among the samples of different treatments. The highest expression of all CBF genes was recorded in the non-acclimated frost-tolerant biotype after freezing stress. Interestingly a significantly higher expression of COR6.6 was detected in cold-acclimated samples of both frost-sensitive and -tolerant biotypes after freezing stress. The presented results provide more insights into freezing tolerance mechanisms in the Camelina plant from both a biochemical point of view and the expression of the associated genes.

Identification of a QTL region for ashy stem blight resistance using genome-wide association and linage analysis in common bean recombinant inbred lines derived from BAT 477 and NY6020-4

Ashy stem blight (ASB), caused by the fungus Macrophomina phaseolina (Tassi) Goidanich is an important disease of the common bean (Phaseolus vulgaris L.). It is important to identify quantitative trait loci (QTL) for ASB resistance and introgress into susceptible cultivars of the common bean. The objective of this research was to identify QTL and single nucleotide polymorphism (SNP) markers associated with ASB resistance in recombinant inbred lines (RIL) derived from a cross between BAT 477 and NY6020-4 common bean. One hundred and twenty-six F_6:7 RIL were phenotyped for ASB in the greenhouse. Disease severity was scored on a scale of 1-9. Genotyping was performed using whole genome resequencing with 2x common bean genome size coverage, and over six million SNPs were obtained. After being filtered, 72,017 SNPs distributed on 11 chromosomes were used to conduct the genome-wide association study (GWAS) and QTL mapping. A novel QTL region of ~4.28 Mbp from 35,546,329 bp to 39,826,434 bp on chromosome Pv03 was identified for ASB resistance. The two SNPs, Chr03_39824257 and Chr03_39824268 located at 39,824,257 bp and 39,824,268 bp on Pv03, respectively, were identified as the strongest markers associated with ASB resistance. The gene Phvul.003G175900 (drought sensitive, WD repeat-containing protein 76) located at 39,822,021 - 39,824,655 bp on Pv03 was recognized as one candidate for ASB resistance in the RIL, and the gene contained the two SNP markers. QTL and SNP markers may be used to select plants and lines for ASB resistance through marker-assisted selection (MAS) in common bean breeding.

Heat Stress Tolerance: A Prerequisite for the Selection of Drought- and Low Phosphorus-Tolerant Common Beans for Equatorial Tropical Regions Such as Ghana

Forty common bean accessions of multiple genetic background trait attribution regarding drought tolerance were selected based on mean yield performance from an earlier field test evaluation conducted using augmented RCBD. The various bean genotypes were further evaluated with phosphorus and water treatment interactions at two different levels for each factor. The experiment was conducted in a 2 × 2 × 40 factorial using RCBD with three replications under screen-house conditions at the CSIR-Crops Research Institute, Kumasi-Ghana. The objective was to select drought- and low phosphorus-tolerant common bean genotypes; which are suitable for tropical climatic conditions. The results showed that common bean with drought and heat trait tolerance survived, developed flowers and podded with seeds to physiological maturity, whilst genotypes with no heat trait tolerance had impaired reproductive structural development and growth disruption; thus, flowers could not develop into pods with seeds. This reproductive developmental anomaly was due to prevailing average daytime and nighttime high temperatures of 35.45 °C and 29.95 °C, respectively, recorded during the growth period, which reduced pollen fertility. Among the 478 experimental bean plants (two plants were missing) analyzed, 141 (29.5%) did not flower, 168 (35.18%) had their pods dropped whilst 99 (20.7%) podded with seeds to achieve physiological maturity. The podded-seed bean genotypes were of the SEF-line pedigrees, which were shown to be heat and drought-tolerant. Meanwhile, bean accessions with SMC, SMN and SMR code prefixes did not pod into seed despite possessing drought-tolerant traits. The effects of interactions between phosphorus and water treatments on the root characteristics of drought-tolerant common bean were as follows: root length, root surface area, average root diameter and root volume growth extensions doubled dimensionally under optimum conditions (P₂W₂) compared to stressed conditions (P₁W₁). The results from the present study identified four SEF-bean genotypes, namely, SEF15, SEF 47, SEF 60 and SEF 62, as superior yield performers, even under low soil phosphorus and in extreme high temperature conditions. Therefore, breeding for the selection of drought- and low-P-tolerant common bean for tropical agro-ecological environments must also consider concomitant heat stress tolerance.

Worldwide Selection Footprints for Drought and Heat in Bread Wheat (Triticum aestivum L.)

Genome–environment Associations (GEA) or Environmental Genome-Wide Association scans (EnvGWAS) have been poorly applied for studying the genomics of adaptive traits in bread wheat landraces (Triticum aestivum L.). We analyzed 990 landraces and seven climatic variables (mean temperature, maximum temperature, precipitation, precipitation seasonality, heat index of mean temperature, heat index of maximum temperature, and drought index) in GEA using the FarmCPU approach with GAPIT. Historical temperature and precipitation values were obtained as monthly averages from 1970 to 2000. Based on 26,064 high-quality SNP loci, landraces were classified into ten subpopulations exhibiting high genetic differentiation. The GEA identified 59 SNPs and nearly 89 protein-encoding genes involved in the response processes to abiotic stress. Genes related to biosynthesis and signaling are mainly mediated by auxins, abscisic acid (ABA), ethylene (ET), salicylic acid (SA), and jasmonates (JA), which are known to operate together in modulation responses to heat stress and drought in plants. In addition, we identified some proteins associated with the response and tolerance to stress by high temperatures, water deficit, and cell wall functions. The results provide candidate regions for selection aimed to improve drought and heat tolerance in bread wheat and provide insights into the genetic mechanisms involved in adaptation to extreme environments.

Genetic diversity and population structure of wild and cultivated Crotalaria species based on genotyping-by-sequencing

Crotalaria is a plant genus that is found all over the world, with over 700 species of herbs and shrubs. The species are potential alternative food and industrial crops due to their adaptability to different environments. Currently, information on the genetic diversity and population structure of these species is scanty. Genotyping-by-sequencing (GBS) is a cost-effective high-throughput technique in diversity evaluation of plant species that have not been fully sequenced. In the current study, de novo GBS was used to characterize 80 Crotalaria accessions from five geographical regions in Kenya. A total of 9820 single nucleotide polymorphism (SNP) markers were obtained after thinning and filtering, which were then used for the analysis of genetic diversity and population structure in Crotalaria. The proportion of SNPs with a minor allele frequency (maf) > = 0.05 was 45.08%, while the Guanine-Cytosine (GC) content was 0.45, from an average sequence depth of 455,909 reads per base. The transition vs transversion ratio was 1.81 and Heterozygosity (He) ranged between 0.01–0.07 in all the sites and 0.04 to 0.52 in the segregating sites. The mean Tajima’s D value for the population was -0.094, suggesting an excess of rare alleles. The fixation index (Fst) between the different populations based on the Wright Fst (1943) ranged from 0.0119 to 0.066 for the Eastern-Western and Nairobi-Western populations. Model based techniques of population structure analysis including structure, k-means and cross-entropy depicted eight clusters in the study accessions. Non-model based techniques especially DAPC depicted poor population stratification. Correspondence Analysis (CA), Principal coordinate analyses (PCoA) and phylogenetic analysis identified a moderate level of population stratification. Results from this study will help conservationists and breeders understand the genetic diversity of Crotalaria. The study also provides valuable information for genetic improvement of domesticated species.

Genome–Environment Associations, an Innovative Tool for Studying Heritable Evolutionary Adaptation in Orphan Crops and Wild Relatives

Leveraging innovative tools to speed up prebreeding and discovery of genotypic sources of adaptation from landraces, crop wild relatives, and orphan crops is a key prerequisite to accelerate genetic gain of abiotic stress tolerance in annual crops such as legumes and cereals, many of which are still orphan species despite advances in major row crops. Here, we review a novel, interdisciplinary approach to combine ecological climate data with evolutionary genomics under the paradigm of a new field of study: genome–environment associations (GEAs). We first exemplify how GEA utilizes in situ georeferencing from genotypically characterized, gene bank accessions to pinpoint genomic signatures of natural selection. We later discuss the necessity to update the current GEA models to predict both regional- and local- or micro-habitat–based adaptation with mechanistic ecophysiological climate indices and cutting-edge GWAS-type genetic association models. Furthermore, to account for polygenic evolutionary adaptation, we encourage the community to start gathering genomic estimated adaptive values (GEAVs) for genomic prediction (GP) and multi-dimensional machine learning (ML) models. The latter two should ideally be weighted by de novo GWAS-based GEA estimates and optimized for a scalable marker subset. We end the review by envisioning avenues to make adaptation inferences more robust through the merging of high-resolution data sources, such as environmental remote sensing and summary statistics of the genomic site frequency spectrum, with the epigenetic molecular functionality responsible for plastic inheritance in the wild. Ultimately, we believe that coupling evolutionary adaptive predictions with innovations in ecological genomics such as GEA will help capture hidden genetic adaptations to abiotic stresses based on crop germplasm resources to assist responses to climate change.

“I shall endeavor to find out how nature’s forces act upon one another, and in what manner the geographic environment exerts its influence on animals and plants. In short, I must find out about the harmony in nature” Alexander von Humboldt—Letter to Karl Freiesleben, June 1799.

Single trait versus principal component based association analysis for flowering related traits in pigeonpea

Pigeonpea, a tropical photosensitive crop, harbors significant diversity for days to flowering, but little is known about the genes that govern these differences. Our goal in the current study was to use genome wide association strategy to discover the loci that regulate days to flowering in pigeonpea. A single trait as well as a principal component based association study was conducted on a diverse collection of 142 pigeonpea lines for days to first and fifty percent of flowering over 3 years, besides plant height and number of seeds per pod. The analysis used seven association mapping models (GLM, MLM, MLMM, CMLM, EMLM, FarmCPU and SUPER) and further comparison revealed that FarmCPU is more robust in controlling both false positives and negatives as it incorporates multiple markers as covariates to eliminate confounding between testing marker and kinship. Cumulatively, a set of 22 SNPs were found to be associated with either days to first flowering (DOF), days to fifty percent flowering (DFF) or both, of which 15 were unique to trait based, 4 to PC based GWAS while 3 were shared by both. Because PC1 represents DOF, DFF and plant height (PH), four SNPs found associated to PC1 can be inferred as pleiotropic. A window of ± 2 kb of associated SNPs was aligned with available transcriptome data generated for transition from vegetative to reproductive phase in pigeonpea. Annotation analysis of these regions revealed presence of genes which might be involved in floral induction like Cytochrome p450 like Tata box binding protein, Auxin response factors, Pin like genes, F box protein, U box domain protein, chromatin remodelling complex protein, RNA methyltransferase. In summary, it appears that auxin responsive genes could be involved in regulating DOF and DFF as majority of the associated loci contained genes which are component of auxin signaling pathways in their vicinity. Overall, our findings indicates that the use of principal component analysis in GWAS is statistically more robust in terms of identifying genes and FarmCPU is a better choice compared to the other aforementioned models in dealing with both false positive and negative associations and thus can be used for traits with complex inheritance.

Identification of Candidate Genes Regulating Drought Tolerance in Pearl Millet

Pearl millet is an important crop of the arid and semi-arid ecologies to sustain food and fodder production. The greater tolerance to drought stress attracts us to examine its cellular and molecular mechanisms via functional genomics approaches to augment the grain yield. Here, we studied the drought response of 48 inbreds representing four different maturity groups at the flowering stage. A set of 74 drought-responsive genes were separated into five major phylogenic groups belonging to eight functional groups, namely ABA signaling, hormone signaling, ion and osmotic homeostasis, TF-mediated regulation, molecular adaptation, signal transduction, physiological adaptation, detoxification, which were comprehensively studied. Among the conserved motifs of the drought-responsive genes, the protein kinases and MYB domain proteins were the most conserved ones. Comparative in-silico analysis of the drought genes across millet crops showed foxtail millet had most orthologs with pearl millet. Of 698 haplotypes identified across millet crops, MyC2 and Myb4 had maximum haplotypes. The protein–protein interaction network identified ABI2, P5CS, CDPK, DREB, MYB, and CYP707A3 as major hub genes. The expression assay showed the presence of common as well as unique drought-responsive genes across maturity groups. Drought tolerant genotypes in respective maturity groups were identified from the expression pattern of genes. Among several gene families, ABA signaling, TFs, and signaling proteins were the prospective contributors to drought tolerance across maturity groups. The functionally validated genes could be used as promising candidates in backcross breeding, genomic selection, and gene-editing schemes in pearl millet and other millet crops to increase the yield in drought-prone arid and semi-arid ecologies.

Genome-Wide Association Study and Genomic Prediction for Bacterial Wilt Resistance in Common Bean (Phaseolus vulgaris) Core Collection

Common bean (Phaseolus vulgaris) is one of the major legume crops cultivated worldwide. Bacterial wilt (BW) of common bean (Curtobacterium flaccumfaciens pv. flaccumfaciens), being a seed-borne disease, has been a challenge in common bean producing regions. A genome-wide association study (GWAS) was conducted to identify SNP markers associated with BW resistance in the USDA common bean core collection. A total of 168 accessions were evaluated for resistance against three different isolates of BW. Our study identified a total of 14 single nucleotide polymorphism (SNP) markers associated with the resistance to BW isolates 528, 557, and 597 using mixed linear models (MLMs) in BLINK, FarmCPU, GAPIT, and TASSEL 5. These SNPs were located on chromosomes Phaseolus vulgaris [Pv]02, Pv04, Pv08, and Pv09 for isolate 528; Pv07, Pv10, and Pv11 for isolate 557; and Pv04, Pv08, and Pv10 for isolate 597. The genomic prediction accuracy was assessed by utilizing seven GP models with 1) all the 4,568 SNPs and 2) the 14 SNP markers. The overall prediction accuracy (PA) ranged from 0.30 to 0.56 for resistance against the three BW isolates. A total of 14 candidate genes were discovered for BW resistance located on chromosomes Pv02, Pv04, Pv07, Pv08, and Pv09. This study revealed vital information for developing genetic resistance against the BW pathogen in common bean. Accordingly, the identified SNP markers and candidate genes can be utilized in common bean molecular breeding programs to develop novel resistant cultivars.

Genome-Wide Identification of the LHC Gene Family in Kiwifruit and Regulatory Role of AcLhcb3.1/3.2 for Chlorophyll a Content

Light-harvesting chlorophyll a/b-binding (LHC) protein is a superfamily that plays a vital role in photosynthesis. However, the reported knowledge of LHCs in kiwifruit is inadequate and poorly understood. In this study, we identified 42 and 45 LHC genes in Actinidia chinensis (Ac) and A. eriantha (Ae) genomes. Phylogenetic analysis showed that the kiwifruit LHCs of both species were grouped into four subfamilies (Lhc, Lil, PsbS, and FCII). Expression profiles and qRT-PCR results revealed expression levels of LHC genes closely related to the light, temperature fluctuations, color changes during fruit ripening, and kiwifruit responses to Pseudomonas syringae pv. actinidiae (Psa). Subcellular localization analysis showed that AcLhcb1.5/3.1/3.2 were localized in the chloroplast while transient overexpression of AcLhcb3.1/3.2 in tobacco leaves confirmed a significantly increased content of chlorophyll a. Our findings provide evidence of the characters and evolution patterns of kiwifruit LHCs genes in kiwifruit and verify the AcLhcb3.1/3.2 genes controlling the chlorophyll a content.

Whole Transcriptome Sequencing Unveils the Genomic Determinants of Putative Somaclonal Variation in Mint (Mentha L.)

Mint (Mentha L., Lamiaceae) is a strongly scented herb of the family Lamiaceae that is grown mostly by clonal propagation, making it a valuable species for the study of somaclonal variation and its phenotypic consequences. The recent introduction of a few species of mint in South America, followed by a presumably rampant propagation, make this region particularly ideal for studying the extent of somaclonal genetic diversity. Hence, the objective of this work was to offer a preliminary characterization of somaclonal genetically coding diversity of the mint in the northern Andes in order to address the question of whether somaclonal variants may have emerged despite relatively recent introductions in a region where mint is not native. A total of 29 clonally propagated specimens, collected in mint export farms in the province of Antioquia, a major region for mint production in the northwest Andes of Colombia, were genotyped using RNA sequencing (RNA-Seq). SNP calling was carried out from the leaves’ transcriptome profiles of each plant by combining the GATK4 and TRINITY protocols, obtaining a total of 2033 loci across 912 transcripts with a minimum read depth of 20X and 4% of missing data. Unsupervised machine learning algorithms considered the K-means, AGNES and UPGMA approaches, all of which suggested three genetic clusters for M. spicata and a unique cluster for M. × piperita. The results indicate that at least two different origins of M. spicata reached the eastern region of the Antioquia province, clonally propagated in the locality ever since for local consumption and export. One of these ancestries had more population structure, possibly due to environmental or anthropological pressures that intervened in the fragmentation of this genetic group or to a higher somaclonal mutation rate. This work offers a first step into the study of the accumulation and transmission of presumably quasi-neutral somatic mutations at coding regions in an herbaceous clonally propagated scented species such as mint, likely favored by an expected population expansion after its Andean introduction. These ad hoc hypotheses warrant further study as part of future research.

Global DNA Methylation and mRNA-miRNA Variations Activated by Heat Shock Boost Early Microspore Embryogenesis in Cabbage (Brassica oleracea)

Microspore culture, a type of haploid breeding, is extensively used in the cultivation of cruciferous crops such as cabbage. Heat shock (HS) treatment is essential to improve the embryo rate during the culture process; however, its molecular role in boosting early microspore embryogenesis (ME) remains unknown. Here we combined DNA methylation levels, miRNAs, and transcriptome profiles in isolated microspores of cabbage ‘01-88’ under HS (32 °C for 24 h) and normal temperature (25 °C for 24 h) to investigate the regulatory roles of DNA methylation and miRNA in early ME. Global methylation levels were significantly different in the two pre-treatments, and 508 differentially methylated regions (DMRs) were identified; 59.92% of DMRs were correlated with transcripts, and 39.43% of miRNA locus were associated with methylation levels. Significantly, the association analysis revealed that 31 differentially expressed genes (DEGs) were targeted by methylation and miRNA and were mainly involved in the reactive oxygen species (ROS) response and abscisic acid (ABA) signaling, indicating that HS induced DNA methylation, and miRNA might affect ME by influencing ROS and ABA. This study revealed that DNA methylation and miRNA interfered with ME by modulating key genes and pathways, which could broaden our understanding of the molecular regulation of ME induced by HS pre-treatment.

Genetic and Genomic Resources for Soybean Breeding Research

Soybean (Glycine max) is a legume species of significant economic and nutritional value. The yield of soybean continues to increase with the breeding of improved varieties, and this is likely to continue with the application of advanced genetic and genomic approaches for breeding. Genome technologies continue to advance rapidly, with an increasing number of high-quality genome assemblies becoming available. With accumulating data from marker arrays and whole-genome resequencing, studying variations between individuals and populations is becoming increasingly accessible. Furthermore, the recent development of soybean pangenomes has highlighted the significant structural variation between individuals, together with knowledge of what has been selected for or lost during domestication and breeding, information that can be applied for the breeding of improved cultivars. Because of this, resources such as genome assemblies, SNP datasets, pangenomes and associated databases are becoming increasingly important for research underlying soybean crop improvement.

Iconography of Beans and Related Legumes Following the Columbian Exchange

Common bean (Phaseolus vulgaris L.), maize, and squash were described by explorers as early as 1492. The illustration of common bean recognized as the first in Europe is in Fuchs' Di Historias Stirpium, published in 1542 and a half-century after beans were observed in the Caribbean. Besides herbals and herbarium specimens, the sources of information on the introduction of New World crops are paintings and illustrations. Two early sources of images of maize and squash are the Grandes Heures d'Anne de Bretagne and the Loggia di Amore e Psiche in the Villa Farnesina, Rome. The former was illustrated between 1507 and 1508 and has an image identified as the common bean. The Villa Farnesina Loggia was decorated in 1515-1518, with festoons containing three instances of bean pods. Our first objective was to evaluate these images to determine whether they represented depictions of common bean earlier than the illustration by Fuchs. Neither image appears to be a common bean based on a combination of botanical characters and size. Folio 194 of the Grandes Heures d'Anne de Bretagne is most likely a Vigna species in the Ceratotropis subgenus. In the Loggia, one set of pods appears to be a species in the Mimosoideae subfamily and the second and third sets of pods most closely resemble Canavalia gladiata. Neither image likely represents common beans and are probably Old-World species. Secondly, illustrations of common beans from ten early herbals were analyzed for traits that are characteristic of the centers of domestication and races of common beans. Our objective was to characterize the diversity observed among herbals and determine whether beans from both centers of domestication were present. We potentially identified both Middle American, race Mesoamerica and Andean, race Nueva Granada types. We posit that both Middle American and Andean types were in the Caribbean at the time of the Columbian exchange and that beans from both centers were informally introduced into Europe early on. This review of 16th-century manuscripts and illustrations has provided some answers to the questions of what and when common beans reached Europe and provide new hypotheses for researchers studying the origins, diversity, and distribution of this crop.

Insights into the plateau adaptation of Salvia castanea by comparative genomic and WGCNA analyses

INTRODUCTION

Salvia castanea, a wild plant species is adapted to extreme Qinghai-Tibetan plateau (QTP) environments. It is also used for medicinal purposes due to high ingredient of tanshinone IIA (T-IIA). Despite its importance to Chinese medicinal industry, the mechanisms associated with secondary metabolites accumulation (i.e. T-IIA and rosmarinic acid (RA)) in this species have not been characterized. Also, the role of special underground tissues in QTP adaptation of S. castanea is still unknown.

OBJECTIVES

We explored the phenomenon of periderm-like structure in underground stem center of S. castanea with an aim to unravel the molecular evolutionary mechanisms of QTP adaptation in this species.

METHODS

Morphologic observation and full-length transcriptome of S. castanea plants were conducted. Comparative genomic analyses of S. castanea with other 14 representative species were used to reveal its phylogenetic position and molecular evolutionary mechanisms. RNA-seq and WGCNA analyses were applied to understand the mechanisms of high accumulations of T-IIA and RA in S. castanea tissues.

RESULTS

Based on anatomical observations, we proposed a "trunk-branches" developmental model to explain periderm-like structure in the center of underground stem of S. castanea. Our study suggested that S. castanea branched off from cultivated Danshen around 16 million years ago. During the evolutionary process, significantly expanded orthologous gene groups, 24 species-specific and 64 positively selected genes contributed to morphogenesis and QTP adaptation in S. castanea. RNA-seq and WGCNA analyses unraveled underlying mechanisms of high accumulations of T-IIA and RA in S. castanea and identified NAC29 and TGA22 as key transcription factors.

CONCLUSION

We proposed a "trunk-branches" developmental model for the underground stem in S. castanea. Adaptations to extreme QTP environment in S. castanea are associated with accumulations of high secondary metabolites in this species.

Whole-Genome Sequencing and Annotation of the Yeast Clavispora santaluciae Reveals Important Insights about Its Adaptation to the Vineyard Environment

Clavispora santaluciae was recently described as a novel non-Saccharomyces yeast species, isolated from grapes of Azores vineyards, a Portuguese archipelago with particular environmental conditions, and from Italian grapes infected with Drosophila suzukii. In the present work, the genome of five Clavispora santaluciae strains was sequenced, assembled, and annotated for the first time, using robust pipelines, and a combination of both long- and short-read sequencing platforms. Genome comparisons revealed specific differences between strains of Clavispora santaluciae reflecting their isolation in two separate ecological niches—Azorean and Italian vineyards—as well as mechanisms of adaptation to the intricate and arduous environmental features of the geographical location from which they were isolated. In particular, relevant differences were detected in the number of coding genes (shared and unique) and transposable elements, the amount and diversity of non-coding RNAs, and the enzymatic potential of each strain through the analysis of their CAZyome. A comparative study was also conducted between the Clavispora santaluciae genome and those of the remaining species of the Metschnikowiaceae family. Our phylogenetic and genomic analysis, comprising 126 yeast strains (alignment of 2362 common proteins) allowed the establishment of a robust phylogram of Metschnikowiaceae and detailed incongruencies to be clarified in the future.

Inheritance of Yield Components and Morphological Traits in Avocado cv. Hass From "Criollo" "Elite Trees" via Half-Sib Seedling Rootstocks

Grafting induces precocity and maintains clonal integrity in fruit tree crops. However, the complex rootstock × scion interaction often precludes understanding how the tree phenotype is shaped, limiting the potential to select optimum rootstocks. Therefore, it is necessary to assess (1) how seedling progenies inherit trait variation from elite 'plus trees', and (2) whether such family superiority may be transferred after grafting to the clonal scion. To bridge this gap, we quantified additive genetic parameters (i.e., narrow sense heritability-h ², and genetic-estimated breeding values-GEBVs) across landraces, "criollo", "plus trees" of the super-food fruit tree crop avocado (Persea americana Mill.), and their open-pollinated (OP) half-sib seedling families. Specifically, we used a genomic best linear unbiased prediction (G-BLUP) model to merge phenotypic characterization of 17 morpho-agronomic traits with genetic screening of 13 highly polymorphic SSR markers in a diverse panel of 104 avocado "criollo" "plus trees." Estimated additive genetic parameters were validated at a 5-year-old common garden trial (i.e., provenance test), in which 22 OP half-sib seedlings from 82 elite "plus trees" served as rootstocks for the cv. Hass clone. Heritability (h ²) scores in the "criollo" "plus trees" ranged from 0.28 to 0.51. The highest h ² values were observed for ribbed petiole and adaxial veins with 0.47 (CI 95%0.2-0.8) and 0.51 (CI 0.2-0.8), respectively. The h ² scores for the agronomic traits ranged from 0.34 (CI 0.2-0.6) to 0.39 (CI 0.2-0.6) for seed weight, fruit weight, and total volume, respectively. When inspecting yield variation across 5-year-old grafted avocado cv. Hass trees with elite OP half-sib seedling rootstocks, the traits total number of fruits and fruits' weight, respectively, exhibited h ² scores of 0.36 (± 0.23) and 0.11 (± 0.09). Our results indicate that elite "criollo" "plus trees" may serve as promissory donors of seedling rootstocks for avocado cv. Hass orchards due to the inheritance of their outstanding trait values. This reinforces the feasibility to leverage natural variation from "plus trees" via OP half-sib seedling rootstock families. By jointly estimating half-sib family effects and rootstock-mediated heritability, this study promises boosting seedling rootstock breeding programs, while better discerning the consequences of grafting in fruit tree crops.

Rootstock-Mediated Genetic Variance in Cadmium Uptake by Juvenile Cacao (Theobroma cacao L.) Genotypes, and Its Effect on Growth and Physiology

Grafting typically offers a shortcut to breed tree orchards throughout a multidimensional space of traits. Despite an overwhelming spectrum of rootstock-mediated effects on scion traits observed across several species, the exact nature and mechanisms underlying the rootstock-mediated effects on scion traits in cacao (Theobroma cacao L.) plants often remain overlooked. Therefore, we aimed to explicitly quantify rootstock-mediated genetic contributions in recombinant juvenile cacao plants across target traits, specifically cadmium (Cd) uptake, and its correlation with growth and physiological traits. Content of chloroplast pigments, fluorescence of chlorophyll a, leaf gas exchange, nutrient uptake, and plant biomass were examined across ungrafted saplings and target rootstock × scion combinations in soils with contrasting levels of Cd. This panel considered a total of 320 progenies from open-pollinated half-sib families and reciprocal full-sib progenies (derived from controlled crosses between the reference genotypes IMC67 and PA121). Both family types were used as rootstocks in grafts with two commercial clones (ICS95 and CCN51) commonly grown in Colombia. A pedigree-based best linear unbiased prediction (A-BLUP) mixed model was implemented to quantify rootstock-mediated narrow-sense heritability (h 2) for target traits. A Cd effect measured on rootstocks before grafting was observed in plant biomass, nutrient uptake, and content of chloroplast pigments. After grafting, damage to the Photosystem II (PSII) was also evident in some rootstock × scion combinations. Differences in the specific combining ability for Cd uptake were mostly detected in ungrafted rootstocks, or 2 months after grafting with the clonal CCN51 scion. Moderate rootstock effects (h 2> 0.1) were detected before grafting for five growth traits, four nutrient uptake properties, and chlorophylls and carotenoids content (h 2 = 0.19, 95% CI 0.05-0.61, r = 0.7). Such rootstock effects faded (h 2< 0.1) when rootstock genotypes were examined in soils without Cd, or 4 months after grafting. These results suggest a pervasive genetic conflict between the rootstock and the scion genotypes, involving the triple rootstock × scion × soil interaction when it refers to Cd and nutrient uptake, early growth, and photosynthetic process in juvenile cacao plants. Overall, deepening on these findings will harness early breeding schemes of cacao rootstock genotypes compatible with commercial clonal scions and adapted to soils enriched with toxic levels of Cd.

Phylogenomic and Microsynteny Analysis Provides Evidence of Genome Arrangements of High-Affinity Nitrate Transporter Gene Families of Plants

Nitrate transporter 2 (NRT2) and NRT3 or nitrate-assimilation-related 2 (NAR2) proteins families form a two-component, high-affinity nitrate transport system, which is essential for the acquisition of nitrate from soils with low N availability. An extensive phylogenomic analysis across land plants for these families has not been performed. In this study, we performed a microsynteny and orthology analysis on the NRT2 and NRT3 genes families across 132 plants (Sensu lato) to decipher their evolutionary history. We identified significant differences in the number of sequences per taxonomic group and different genomic contexts within the NRT2 family that might have contributed to N acquisition by the plants. We hypothesized that the greater losses of NRT2 sequences correlate with specialized ecological adaptations, such as aquatic, epiphytic, and carnivory lifestyles. We also detected expansion on the NRT2 family in specific lineages that could be a source of key innovations for colonizing contrasting niches in N availability. Microsyntenic analysis on NRT3 family showed a deep conservation on land plants, suggesting a high evolutionary constraint to preserve their function. Our study provides novel information that could be used as guide for functional characterization of these gene families across plant lineages.

Genome-Environment Association Analysis for Bio-Climatic Variables in Common Bean (Phaseolus vulgaris L.) from Brazil

Abiotic stress is a limiting factor for common bean (Phaseolus vulgaris L.) production globally. The study of the genotypic, phenotypic, and bio-climatic variables in a broad set of accessions may assist the identification of genomic regions involved in the climatic adaptation of the common bean. We conducted a genotyping-by-sequencing analysis using 28,823 SNPs on 110 georeferenced common bean accessions from Brazil to discover associations between SNPs and bio-climatic indexes. The population structure analysis clustered the accessions into two groups corresponding to the Andean and Mesoamerican gene pools. Of the 19 bioclimatic variables, 17 exhibited a significant association with SNPs on chromosomes Pv01, Pv02, Pv03, Pv04, Pv06, Pv09, Pv10, and Pv11 of common bean. Ten candidate genes were associated with specific bio-climatic variables related to temperature and precipitation. The candidate genes associated with this significant Pv09 region encode a Platz transcription factor family protein previously reported to be an essential regulator of drought stress. The SNP markers and candidate genes associated with the bio-climatic variables should be validated in segregating populations for water stress, which could further be used for marker-assisted selection. As a result, bean breeding programs may be able to provide advances in obtaining drought-tolerant cultivars.

Genome-Wide Association Study and Genomic Prediction for Soybean Cyst Nematode Resistance in USDA Common Bean (Phaseolus vulgaris) Core Collection

Soybean cyst nematode (SCN, Heterodera glycines) has become the major yield-limiting biological factor in soybean production. Common bean is also a good host of SCN, and its production is challenged by this emerging pest in many regions such as the upper Midwest USA. The use of host genetic resistance has been the most effective and environmentally friendly method to manage SCN. The objectives of this study were to evaluate the SCN resistance in the USDA common bean core collection and conduct a genome-wide association study (GWAS) of single nucleotide polymorphism (SNP) markers with SCN resistance. A total of 315 accessions of the USDA common bean core collection were evaluated for resistance to SCN HG Type 0 (race 6). The common bean core set was genotyped with the BARCBean6K_3 Infinium BeadChips, consisting of 4,654 SNPs. Results showed that 15 accessions were resistant to SCN with a Female Index (FI) at 4.8 to 9.4, and 62 accessions were moderately resistant (10 < FI < 30) to HG Type 0. The association study showed that 11 SNP markers, located on chromosomes Pv04, 07, 09, and 11, were strongly associated with resistance to HG Type 0. GWAS was also conducted for resistance to HG Type 2.5.7 and HG Type 1.2.3.5.6.7 based on the public dataset (N = 276), consisting of a diverse set of common bean accessions genotyped with the BARCBean6K_3 chip. Six SNPs associated with HG Type 2.5.7 resistance on Pv 01, 02, 03, and 07, and 12 SNPs with HG Type 1.2.3.5.6.7 resistance on Pv 01, 03, 06, 07, 09, 10, and 11 were detected. The accuracy of genomic prediction (GP) was 0.36 to 0.49 for resistance to the three SCN HG types, indicating that genomic selection (GS) of SCN resistance is feasible. This study provides basic information for developing SCN-resistant common bean cultivars, using the USDA core germ plasm accessions. The SNP markers can be used in molecular breeding in common beans through marker-assisted selection (MAS) and GS.

Sequencing and Comparative Genomic Analysis of a Highly Metal-Tolerant Penicillium janthinellum P1 Provide Insights Into Its Metal Tolerance

Heavy metal pollution is a global knotty problem and fungi hold promising potential for the remediation of wastewater containing heavy metals. Here, a new highly chromium-tolerance species, Penicillium janthinellum P1, is investigated. The genome of P1 was sequenced and assembled into 30 Mb genome size containing 10,955 predicted protein-coding genes with a GC content of 46.16% through an integrated method of Illumina short-read sequencing and single-molecule real-time Pacific Biosciences sequencing platforms. Through a phylogenetic analysis with model species of fungi, the evolutionary divergence time of Penicillium janthinellum P1 and Penicillium oxalicum 114-2 was estimated to be 74 MYA. 33 secondary metabolism gene clusters were identified via antiSMASH software, mainly including non-ribosomal peptide synthase genes and T1 polyketide synthase genes. 525 genes were annotated to encode enzymes that act on carbohydrates, involving 101 glucose-degrading enzymes and 24 polysaccharide synthase. By whole-genome sequence analysis, large numbers of metal resistance genes were found in strain P1. Especially ABC transporter and Superoxide dismutase ensure that the P1 fungus can survive in a chromium-polluted environment. ChrA and ChrR were also identified as key genes for chromium resistance. Analysis of their genetic loci revealed that the specific coding-gene arrangement may account for the fungus’s chromium resistance. Genetic information and comparative analysis of Penicillium janthinellum are valuable for further understanding the mechanism of high resistance to heavy metal chromium, and gene loci analysis provides a new perspective for identifying chromium-resistant strains.

Gene Flow in Phaseolus Beans and Its Role as a Plausible Driver of Ecological Fitness and Expansion of Cultigens

The genus Phaseolus, native to the Americas, is composed of more than eighty wild species, five of which were domesticated in pre-Columbian times. Since the beginning of domestication events in this genus, ample opportunities for gene flow with wild relatives have existed. The present work reviews the extent of gene flow in the genus Phaseolus in primary and secondary areas of domestication with the aim of illustrating how this evolutionary force may have conditioned ecological fitness and the widespread adoption of cultigens. We focus on the biological bases of gene flow in the genus Phaseolus from a spatial and time perspective, the dynamics of wild-weedy-crop complexes in the common bean and the Lima bean, the two most important domesticated species of the genus, and the usefulness of genomic tools to detect inter and intraspecific introgression events. In this review we discuss the reproductive strategies of several Phaseolus species, the factors that may favor outcrossing rates and evidence suggesting that interspecific gene flow may increase ecological fitness of wild populations. We also show that wild-weedy-crop complexes generate genetic diversity over which farmers are able to select and expand their cultigens outside primary areas of domestication. Ultimately, we argue that more studies are needed on the reproductive biology of the genus Phaseolus since for most species breeding systems are largely unknown. We also argue that there is an urgent need to preserve wild-weedy-crop complexes and characterize the genetic diversity generated by them, in particular the genome-wide effects of introgressions and their value for breeding programs. Recent technological advances in genomics, coupled with agronomic characterizations, may make a large contribution.

Harnessing Crop Wild Diversity for Climate Change Adaptation

Warming and drought are reducing global crop production with a potential to substantially worsen global malnutrition. As with the green revolution in the last century, plant genetics may offer concrete opportunities to increase yield and crop adaptability. However, the rate at which the threat is happening requires powering new strategies in order to meet the global food demand. In this review, we highlight major recent 'big data' developments from both empirical and theoretical genomics that may speed up the identification, conservation, and breeding of exotic and elite crop varieties with the potential to feed humans. We first emphasize the major bottlenecks to capture and utilize novel sources of variation in abiotic stress (i.e., heat and drought) tolerance. We argue that adaptation of crop wild relatives to dry environments could be informative on how plant phenotypes may react to a drier climate because natural selection has already tested more options than humans ever will. Because isolated pockets of cryptic diversity may still persist in remote semi-arid regions, we encourage new habitat-based population-guided collections for genebanks. We continue discussing how to systematically study abiotic stress tolerance in these crop collections of wild and landraces using geo-referencing and extensive environmental data. By uncovering the genes that underlie the tolerance adaptive trait, natural variation has the potential to be introgressed into elite cultivars. However, unlocking adaptive genetic variation hidden in related wild species and early landraces remains a major challenge for complex traits that, as abiotic stress tolerance, are polygenic (i.e., regulated by many low-effect genes). Therefore, we finish prospecting modern analytical approaches that will serve to overcome this issue. Concretely, genomic prediction, machine learning, and multi-trait gene editing, all offer innovative alternatives to speed up more accurate pre- and breeding efforts toward the increase in crop adaptability and yield, while matching future global food demands in the face of increased heat and drought. In order for these 'big data' approaches to succeed, we advocate for a trans-disciplinary approach with open-source data and long-term funding. The recent developments and perspectives discussed throughout this review ultimately aim to contribute to increased crop adaptability and yield in the face of heat waves and drought events.

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.

Allelic Diversity at Abiotic Stress Responsive Genes in Relationship to Ecological Drought Indices for Cultivated Tepary Bean, Phaseolus acutifolius A. Gray, and Its Wild Relatives

Some of the major impacts of climate change are expected in regions where drought stress is already an issue. Grain legumes are generally drought susceptible. However, tepary bean and its wild relatives within Phaseolus acutifolius or P. parvifolius are from arid areas between Mexico and the United States. Therefore, we hypothesize that these bean accessions have diversity signals indicative of adaptation to drought at key candidate genes such as: Asr2, Dreb2B, and ERECTA. By sequencing alleles of these genes and comparing to estimates of drought tolerance indices from climate data for the collection site of geo-referenced, tepary bean accessions, we determined the genotype x environmental association (GEA) of each gene. Diversity analysis found that cultivated and wild P. acutifolius were intermingled with var. tenuifolius and P. parvifolius, signifying that allele diversity was ample in the wild and cultivated clade over a broad sense (sensu lato) evaluation. Genes Dreb2B and ERECTA harbored signatures of directional selection, represented by six SNPs correlated with the environmental drought indices. This suggests that wild tepary bean is a reservoir of novel alleles at genes for drought tolerance, as expected for a species that originated in arid environments. Our study corroborated that candidate gene approach was effective for marker validation across a broad genetic base of wild tepary accessions.

Downstream Signalling from Molecular Hydrogen

Molecular hydrogen (H2) is now considered part of the suite of small molecules that can control cellular activity. As such, H2 has been suggested to be used in the therapy of diseases in humans and in plant science to enhance the growth and productivity of plants. Treatments of plants may involve the creation of hydrogen-rich water (HRW), which can then be applied to the foliage or roots systems of the plants. However, the molecular action of H2 remains elusive. It has been suggested that the presence of H2 may act as an antioxidant or on the antioxidant capacity of cells, perhaps through the scavenging of hydroxyl radicals. H2 may act through influencing heme oxygenase activity or through the interaction with reactive nitrogen species. However, controversy exists around all the mechanisms suggested. Here, the downstream mechanisms in which H2 may be involved are critically reviewed, with a particular emphasis on the H2 mitigation of stress responses. Hopefully, this review will provide insight that may inform future research in this area.

Modern Strategies to Assess and Breed Forest Tree Adaptation to Changing Climate

Studying the genetics of adaptation to new environments in ecologically and industrially important tree species is currently a major research line in the fields of plant science and genetic improvement for tolerance to abiotic stress. Specifically, exploring the genomic basis of local adaptation is imperative for assessing the conditions under which trees will successfully adapt in situ to global climate change. However, this knowledge has scarcely been used in conservation and forest tree improvement because woody perennials face major research limitations such as their outcrossing reproductive systems, long juvenile phase, and huge genome sizes. Therefore, in this review we discuss predictive genomic approaches that promise increasing adaptive selection accuracy and shortening generation intervals. They may also assist the detection of novel allelic variants from tree germplasm, and disclose the genomic potential of adaptation to different environments. For instance, natural populations of tree species invite using tools from the population genomics field to study the signatures of local adaptation. Conventional genetic markers and whole genome sequencing both help identifying genes and markers that diverge between local populations more than expected under neutrality, and that exhibit unique signatures of diversity indicative of "selective sweeps." Ultimately, these efforts inform the conservation and breeding status capable of pivoting forest health, ecosystem services, and sustainable production. Key long-term perspectives include understanding how trees' phylogeographic history may affect the adaptive relevant genetic variation available for adaptation to environmental change. Encouraging "big data" approaches (machine learning-ML) capable of comprehensively merging heterogeneous genomic and ecological datasets is becoming imperative, too.

Predicting Thermal Adaptation by Looking Into Populations' Genomic Past

Molecular evolution offers an insightful theory to interpret the genomic consequences of thermal adaptation to previous events of climate change beyond range shifts. However, disentangling often mixed footprints of selective and demographic processes from those due to lineage sorting, recombination rate variation, and genomic constrains is not trivial. Therefore, here we condense current and historical population genomic tools to study thermal adaptation and outline key developments (genomic prediction, machine learning) that might assist their utilization for improving forecasts of populations' responses to thermal variation. We start by summarizing how recent thermal-driven selective and demographic responses can be inferred by coalescent methods and in turn how quantitative genetic theory offers suitable multi-trait predictions over a few generations via the breeder's equation. We later assume that enough generations have passed as to display genomic signatures of divergent selection to thermal variation and describe how these footprints can be reconstructed using genome-wide association and selection scans or, alternatively, may be used for forward prediction over multiple generations under an infinitesimal genomic prediction model. Finally, we move deeper in time to comprehend the genomic consequences of thermal shifts at an evolutionary time scale by relying on phylogeographic approaches that allow for reticulate evolution and ecological parapatric speciation, and end by envisioning the potential of modern machine learning techniques to better inform long-term predictions. We conclude that foreseeing future thermal adaptive responses requires bridging the multiple spatial scales of historical and predictive environmental change research under modern cohesive approaches such as genomic prediction and machine learning frameworks.