Strategies to Increase Prediction Accuracy in Genomic Selection of Complex Traits in Alfalfa (Medicago sativa L.)

Anther-specific expression of MsMYB35 transcription factor in alfalfa (Medicago sativa L.) and its crucial role in pollen development

Alfalfa (Medicago sativa L.) is a high-quality forage crop and an essential resource for livestock. Understanding the molecular mechanisms underlying male sterility in alfalfa is pivotal for the development of superior forage varieties. Despite the critical role of anther development in plant reproduction, its molecular regulation-particularly the involvement of transcription factors in M. sativa-remains insufficiently explored. This study bridges this gap by isolating and characterizing an R2R3-MYB transcription factor, MsMYB35, and unveiling its regulatory role in anther development. Quantitative RT-PCR (qRT-PCR) revealed that MsMYB35 is predominantly expressed during early anther development and is homologous to AtMYB35. MsMYB35 was found to localize in both the cytoplasm and nucleus. DNA affinity purification sequencing (DAP-seq) identified 3647 target genes of MsMYB35, with enrichment analysis uncovering three recognition motifs. Integrated DAP-seq and RNA-seq analyses revealed that MsMYB35 directly regulates two key anther development-related genes. Functional analyses showed that overexpression of MsMYB35 promotes anther development, while silencing MsMYB35 leads to defective anther sacs and wrinkled pollen grains. Proper MsMYB35 expression ensures the formation of viable and fertile pollen grains, solidifying its role as a critical regulator of anther development. These findings provide a novel perspective on the molecular mechanisms regulating anther development in M. sativa and offer valuable insights for improving molecular breeding and hybrid seed production strategies. By advancing the fundamental understanding of transcriptional regulation in anther development, this study sets the stage for innovative approaches to alfalfa crop improvement.

Alfalfa Mosaic Virus and White Clover Mosaic Virus Combined Infection Leads to Chloroplast Destruction and Alterations in Photosynthetic Characteristics of Nicotiana benthamiana

Alfalfa mosaic virus (AMV) is one of the most widely distributed viruses; it often exhibits combined infection with white clover mosaic virus (WCMV). Even so, little is known about the effects of co-infection with AMV and WCMV on plants. To determine whether there is a synergistic effect of AMV and WCMV co-infection, virus co-infection was studied by electron microscopy, the double-antibody sandwich enzyme-linked immunosorbent assay (DAS-ELISA), and real-time fluorescence quantitative PCR (RT-qPCR) of AMV and WCMV co-infection in Nicotiana benthamiana. Meanwhile, measurements were carried out on the photosynthetic pigments, photosynthetic gas exchange parameters, and chlorophyll fluorescence parameters. The results showed that the most severe disease development was induced by AMV and WCMV co-infection, and the disease grade was scale 7. N. benthamiana leaves induced mottled yellow-green alternating patterns, leaf wrinkling, and chlorosis, and chloroplasts were observed to be on the verge of disintegration. The relative accumulation of AMV CP and WCMV CP was significantly increased by 15.44-fold and 10.04-fold upon co-infection compared to that with AMV and WCMV single infection at 21 dpi. In addition, chlorophyll a, chlorophyll b, total chlorophyll, the net photosynthetic rate, the water use efficiency, the apparent electron transport rate, the PSII maximum photochemical efficiency, the actual photochemical quantum yield, and photochemical quenching were significantly reduced in leaves co-infected with AMV and WCMV compared to AMV- or WCMV-infected leaves and CK. On the contrary, the carotenoid content, transpiration rate, stomatal conductance, intercellular CO2 concentration, minimal fluorescence value, and non-photochemical quenching were significantly increased. These findings suggest that there was a synergistic effect between AMV and WCMV, and AMV and WCMV co-infection severely impacted the normal function of photosynthesis in N. benthamiana.

Genome-wide family prediction unveils molecular mechanisms underlying the regulation of agronomic traits in Urochloa ruziziensis

Tropical forage grasses, particularly those belonging to the Urochloa genus, play a crucial role in cattle production and serve as the main food source for animals in tropical and subtropical regions. The majority of these species are apomictic and tetraploid, highlighting the significance of U. ruziziensis, a sexual diploid species that can be tetraploidized for use in interspecific crosses with apomictic species. As a means to support breeding programs, our study investigates the feasibility of genome-wide family prediction in U. ruziziensis families to predict agronomic traits. Fifty half-sibling families were assessed for green matter yield, dry matter yield, regrowth capacity, leaf dry matter, and stem dry matter across different clippings established in contrasting seasons with varying available water capacity. Genotyping was performed using a genotyping-by-sequencing approach based on DNA samples from family pools. In addition to conventional genomic prediction methods, machine learning and feature selection algorithms were employed to reduce the necessary number of markers for prediction and enhance predictive accuracy across phenotypes. To explore the regulation of agronomic traits, our study evaluated the significance of selected markers for prediction using a tree-based approach, potentially linking these regions to quantitative trait loci (QTLs). In a multiomic approach, genes from the species transcriptome were mapped and correlated to those markers. A gene coexpression network was modeled with gene expression estimates from a diverse set of U. ruziziensis genotypes, enabling a comprehensive investigation of molecular mechanisms associated with these regions. The heritabilities of the evaluated traits ranged from 0.44 to 0.92. A total of 28,106 filtered SNPs were used to predict phenotypic measurements, achieving a mean predictive ability of 0.762. By employing feature selection techniques, we could reduce the dimensionality of SNP datasets, revealing potential genotype-phenotype associations. The functional annotation of genes near these markers revealed associations with auxin transport and biosynthesis of lignin, flavonol, and folic acid. Further exploration with the gene coexpression network uncovered associations with DNA metabolism, stress response, and circadian rhythm. These genes and regions represent important targets for expanding our understanding of the metabolic regulation of agronomic traits and offer valuable insights applicable to species breeding. Our work represents an innovative contribution to molecular breeding techniques for tropical forages, presenting a viable marker-assisted breeding approach and identifying target regions for future molecular studies on these agronomic traits.

Simultaneous improvement of grain yield and grain protein concentration in durum wheat by using association tests and weighted GBLUP

KEY MESSAGE

Simultaneous improvement for GY and GPC by using GWAS and GBLUP suggested a significant application in durum wheat breeding. Despite the importance of grain protein concentration (GPC) in determining wheat quality, its negative correlation with grain yield (GY) is still one of the major challenges for breeders. Here, a durum wheat panel of 200 genotypes was evaluated for GY, GPC, and their derived indices (GPD and GYD), under eight different agronomic conditions. The plant material was genotyped with the Illumina 25 k iSelect array, and a genome-wide association study was performed. Two statistical models revealed dozens of marker-trait associations (MTAs), each explaining up to 30%. phenotypic variance. Two markers on chromosomes 2A and 6B were consistently identified by both models and were found to be significantly associated with GY and GPC. MTAs identified for phenological traits co-mapped to well-known genes (i.e., Ppd-1, Vrn-1). The significance values (p-values) that measure the strength of the association of each single nucleotide polymorphism marker with the target traits were used to perform genomic prediction by using a weighted genomic best linear unbiased prediction model. The trained models were ultimately used to predict the agronomic performances of an independent durum wheat panel, confirming the utility of genomic prediction, although environmental conditions and genetic backgrounds may still be a challenge to overcome. The results generated through our study confirmed the utility of GPD and GYD to mitigate the inverse GY and GPC relationship in wheat, provided novel markers for marker-assisted selection and opened new ways to develop cultivars through genomic prediction approaches.

Genome-wide genotyping data renew knowledge on genetic diversity of a worldwide alfalfa collection and give insights on genetic control of phenology traits

China's and Europe's dependence on imported protein is a threat to the food self-sufficiency of these regions. It could be solved by growing more legumes, including alfalfa that is the highest protein producer under temperate climate. To create productive and high-value varieties, the use of large genetic diversity combined with genomic evaluation could improve current breeding programs. To study alfalfa diversity, we have used a set of 395 alfalfa accessions (i.e. populations), mainly from Europe, North and South America and China, with fall dormancy ranging from 3 to 7 on a scale of 11. Five breeders provided materials (617 accessions) that were compared to the 400 accessions. All accessions were genotyped using Genotyping-by-Sequencing (GBS) to obtain SNP allele frequency. These genomic data were used to describe genetic diversity and identify genetic groups. The accessions were phenotyped for phenology traits (fall dormancy and flowering date) at two locations (Lusignan in France, Novi Sad in Serbia) from 2018 to 2021. The QTL were detected by a Multi-Locus Mixed Model (mlmm). Subsequently, the quality of the genomic prediction for each trait was assessed. Cross-validation was used to assess the quality of prediction by testing GBLUP, Bayesian Ridge Regression (BRR), and Bayesian Lasso methods. A genetic structure with seven groups was found. Most of these groups were related to the geographical origin of the accessions and showed that European and American material is genetically distinct from Chinese material. Several QTL associated with fall dormancy were found and most of these were linked to genes. In our study, the infinitesimal methods showed a higher prediction quality than the Bayesian Lasso, and the genomic prediction achieved high (>0.75) predicting abilities in some cases. Our results are encouraging for alfalfa breeding by showing that it is possible to achieve high genomic prediction quality.

Combining GS-assisted GWAS and transcriptome analysis to mine candidate genes for nitrogen utilization efficiency in Populus cathayana

BACKGROUND

Forest trees such as poplar, shrub willow, et al. are essential natural resources for sustainable and renewable energy production, and their wood can reduce dependence on fossil fuels and reduce environmental pollution. However, the productivity of forest trees is often limited by the availability of nitrogen (N), improving nitrogen use efficiency (NUE) is an important way to address it. Currently, NUE genetic resources are scarce in forest tree research, and more genetic resources are urgently needed.

RESULTS

Here, we performed genome-wide association studies (GWAS) using the mixed linear model (MLM) to identify genetic loci regulating growth traits in Populus cathayana at two N levels, and attempted to enhance the signal strength of single nucleotide polymorphism (SNP) detection by performing genome selection (GS) assistance GWAS. The results of the two GWAS analyses identified 55 and 40 SNPs that were respectively associated with plant height (PH) and ground diameter (GD), and 92 and 69 candidate genes, including 30 overlapping genes. The prediction accuracy of the GS model (rrBLUP) for phenotype exceeds 0.9. Transcriptome analysis of 13 genotypes under two N levels showed that genes related to carbon and N metabolism, amino acid metabolism, energy metabolism, and signal transduction were differentially expressed in the xylem of P. cathayana under N treatment. Furthermore, we observed strong regional patterns in gene expression levels of P. cathayana, with significant differences between different regions. Among them, P. cathayana in Longquan region exhibited the highest response to N. Finally, through weighted gene co-expression network analysis (WGCNA), we identified a module closely related to the N metabolic process and eight hub genes.

CONCLUSIONS

Integrating the GWAS, RNA-seq and WGCNA data, we ultimately identified four key regulatory genes (PtrNAC123, PtrNAC025, Potri.002G233100, and Potri.006G236200) involved in the wood formation process, and they may affect P. cathayana growth and wood formation by regulating nitrogen metabolism. This study will provide strong evidence for N regulation mechanisms, and reliable genetic resources for growth and NUE genetic improvement in poplar.

Alfalfa genomic selection for different stress-prone growing regions

Alfalfa (Medicago sativa L.) selection for stress-prone regions has high priority for sustainable crop-livestock systems. This study assessed the genomic selection (GS) ability to predict alfalfa breeding values for drought-prone agricultural sites of Algeria, Morocco, and Argentina; managed-stress (MS) environments of Italy featuring moderate or intense drought; and one Tunisian site irrigated with moderately saline water. Additional aims were to investigate genotype × environment interaction (GEI) patterns and the effect on GS predictions of three single-nucleotide polymorphism (SNP) calling procedures, 12 statistical models that exclude or incorporate GEI, and allele dosage information. Our study included 127 genotypes from a Mediterranean reference population originated from three geographically contrasting populations, genotyped via genotyping-by-sequencing and phenotyped based on multi-year biomass dry matter yield of their dense-planted half-sib progenies. The GEI was very large, as shown by 27-fold greater additive genetic variance × environment interaction relative to the additive genetic variance and low genetic correlation for progeny yield responses across environments. The predictive ability of GS (using at least 37,969 SNP markers) exceeded 0.20 for moderate MS (representing Italian stress-prone sites) and the sites of Algeria and Argentina while being quite low for the Tunisian site and intense MS. Predictions of GS were complicated by rapid linkage disequilibrium decay. The weighted GBLUP model, GEI incorporation into GS models, and SNP calling based on a mock reference genome exhibited a predictive ability advantage for some environments. Our results support the specific breeding for each target region and suggest a positive role for GS in most regions when considering the challenges associated with phenotypic selection.

Application of machine learning to explore the genomic prediction accuracy of fall dormancy in autotetraploid alfalfa

Fall dormancy (FD) is an essential trait to overcome winter damage and for alfalfa (Medicago sativa) cultivar selection. The plant regrowth height after autumn clipping is an indirect way to evaluate FD. Transcriptomics, proteomics, and quantitative trait locus mapping have revealed crucial genes correlated with FD; however, these genes cannot predict alfalfa FD very well. Here, we conducted genomic prediction of FD using whole-genome SNP markers based on machine learning-related methods, including support vector machine (SVM) regression, and regularization-related methods, such as Lasso and ridge regression. The results showed that using SVM regression with linear kernel and the top 3000 genome-wide association study (GWAS)-associated markers achieved the highest prediction accuracy for FD of 64.1%. For plant regrowth height, the prediction accuracy was 59.0% using the 3000 GWAS-associated markers and the SVM linear model. This was better than the results using whole-genome markers (25.0%). Therefore, the method we explored for alfalfa FD prediction outperformed the other models, such as Lasso and ElasticNet. The study suggests the feasibility of using machine learning to predict FD with GWAS-associated markers, and the GWAS-associated markers combined with machine learning would benefit FD-related traits as well. Application of the methodology may provide potential targets for FD selection, which would accelerate genetic research and molecular breeding of alfalfa with optimized FD.