Genetic diversity and local adaption of alfalfa populations (Medicago sativa L.) under long-term grazing.

Genomic information on alfalfa adaptation to long-term grazing is useful for alfalfa genetic improvement. In this study, 14 alfalfa populations were collected from long-term grazing sites (> 25 years) across four soil zones in western Canada. Alfalfa cultivars released between 1926 and 1980 were used to compare degree of genetic variation of the 14 populations. Six agro-morphological and three nutritive value traits were evaluated from 2018 to 2020. The genotyping-by-sequencing (GBS) data of the alfalfa populations and environmental data were used for genotype-environment association (GEA). Both STRUCTURE and UPGMA based on 19,853 SNPs showed that the 14 alfalfa populations from long-term grazing sites had varying levels of parentages from alfalfa sub-species Medicago sativa and M. falcata. The linear regression of STRUCTURE membership probability on phenotypic data indicated genetic variations of forage dry matter yield, spring vigor and plant height were low, but genetic variations of regrowth, fall plant height, days to flower and crude protein were still high for the 14 alfalfa populations from long-term grazing sites. The GEA identified 31 SNPs associated with 13 candidate genes that were mainly associated with six environmental factors of. Candidate genes underlying environmental factors were associated with a variety of proteins, which were involved in plant responses to abiotic stresses, i.e., drought, cold and salinity-alkali stresses.

Transcriptomic analysis of differentially expressed genes in leaves and roots of two alfalfa (Medicago sativa L.) cultivars with different salt tolerance.

BACKGROUND: Alfalfa (Medicago sativa L.) production decreases under salt stress. Identification of genes associated with salt tolerance in alfalfa is essential for the development of molecular markers used for breeding and genetic improvement. RESULT: An RNA-Seq technique was applied to identify the differentially expressed genes (DEGs) associated with salt stress in two alfalfa cultivars: salt tolerant 'Halo' and salt intolerant 'Vernal'. Leaf and root tissues were sampled for RNA extraction at 0 h, 3 h, and 27 h under 12 dS m- 1 salt stress maintained by NaCl. The sequencing generated a total of 381 million clean sequence reads and 84.8% were mapped on to the alfalfa reference genome. A total of 237 DEGs were identified in leaves and 295 DEGs in roots of the two alfalfa cultivars. In leaf tissue, the two cultivars had a similar number of DEGs at 3 h and 27 h of salt stress, with 31 and 49 DEGs for 'Halo', 34 and 50 for 'Vernal', respectively. In root tissue, 'Halo' maintained 55 and 56 DEGs at 3 h and 27 h, respectively, while the number of DEGs decreased from 42 to 10 for 'Vernal'. This differential expression pattern highlights different genetic responses of the two cultivars to salt stress at different time points. Interestingly, 28 (leaf) and 31 (root) salt responsive candidate genes were highly expressed in 'Halo' compared to 'Vernal' under salt stress, of which 13 candidate genes were common for leaf and root tissues. About 60% of DEGs were assigned to known gene ontology (GO) categories. The genes were involved in transmembrane protein function, photosynthesis, carbohydrate metabolism, defense against oxidative damage, cell wall modification and protection against lipid peroxidation. Ion binding was found to be a key molecular activity for salt tolerance in alfalfa under salt stress. CONCLUSION: The identified DEGs are significant for understanding the genetic basis of salt tolerance in alfalfa. The generated genomic information is useful for molecular marker development for alfalfa genetic improvement for salt tolerance.

Tissue specific changes in elements and organic compounds of alfalfa (Medicago sativa L.) cultivars differing in salt tolerance under salt stress.

Soil salinity is a global concern and often the primary factor contributing to land degradation, limiting crop growth and production. Alfalfa (Medicago sativa L.) is a low input high value forage legume with a wide adaptation. Examining the tissue-specific responses to salt stress will be important to understanding physiological changes of alfalfa. The responses of two alfalfa cultivars (salt tolerant 'Halo', salt intolerant 'Vernal') were studied for 12 weeks in five gradients of salt stress in a sand based hydroponic system in the greenhouse. The accumulation and localization of elements and organic compounds in different tissues of alfalfa under salt stress were evaluated using synchrotron beamlines. The pattern of chlorine accumulation for 'Halo' was: root > stem ~ leaf at 8 dSm-1, and root ~ leaf > stem at 12 dSm-1, potentially preventing toxic ion accumulation in leaf tissues. In contrast, for 'Vernal', it was leaf > stem ~ root at 8 dSm-1 and leaf > root ~ stem at 12 dSm-1. The distribution of chlorine in 'Halo' was relatively uniform in the leaf surface and vascular bundles of the stem. Amide concentration in the leaf and stem tissues was greater for 'Halo' than 'Vernal' at all salt gradients. This study determined that low ion accumulation in the shoot was a common strategy in salt tolerant alfalfa up to 8 dSm-1 of salt stress, which was then replaced by shoot tissue tolerance at 12 dSm-1.

Advancing crested wheatgrass [Agropyron cristatum (L.) Gaertn.] breeding through genotyping-by-sequencing and genomic selection.

Crested wheatgrass [Agropyron cristatum (L.) Gaertn.] provides high quality, highly palatable forage for early season grazing. Genetic improvement of crested wheatgrass has been challenged by its complex genome, outcrossing nature, long breeding cycle, and lack of informative molecular markers. Genomic selection (GS) has potential for improving traits of perennial forage species, and genotyping-by-sequencing (GBS) has enabled the development of genome-wide markers in non-model polyploid plants. An attempt was made to explore the utility of GBS and GS in crested wheatgrass breeding. Sequencing and phenotyping 325 genotypes representing 10 diverse breeding lines were performed. Bioinformatics analysis identified 827, 3,616, 14,090 and 46,136 single nucleotide polymorphism markers at 20%, 30%, 40% and 50% missing marker levels, respectively. Four GS models (BayesA, BayesB, BayesCπ, and rrBLUP) were examined for the accuracy of predicting nine agro-morphological and three nutritive value traits. Moderate accuracy (0.20 to 0.32) was obtained for the prediction of heading days, leaf width, plant height, clump diameter, tillers per plant and early spring vigor for genotypes evaluated at Saskatoon, Canada. Similar accuracy (0.29 to 0.35) was obtained for predicting fall regrowth and plant height for genotypes evaluated at Swift Current, Canada. The Bayesian models displayed similar or higher accuracy than rrBLUP. These findings show the feasibility of GS application for a non-model species to advance plant breeding.

Advancing Bromegrass Breeding Through Imaging Phenotyping and Genomic Selection: A Review.

Breeding forage crops for high yields of digestible biomass along with improved resource-use efficiency and wide adaptation is essential to meet future challenges in forage production imposed by growing demand, declining resources, and changing climate. Bromegrasses (Bromus spp.) are economically important forage species in the temperate regions of world, but genetic gain in forage yield of bromegrass is relatively low. In particular, limited breeding efforts have been made in improving abiotic stress tolerance and resource-use efficiency. We conducted a literature review on bromegrass breeding achievements and challenges, global climate change impacts on bromegrass species, and explored the feasibility of applying high-throughput imaging phenotyping techniques and genomic selection for further advances in forage yield and quality selection. Overall genetic gain in forage yield of bromegrass has been low, but genetic improvement in forage yield of smooth bromegrass (Bromus inermis Leyss) is somewhat higher than that of meadow bromegrass (Bromus riparius Rehm). This low genetic gain in bromegrass yield is due to a few factors such as its genetic complexity, lack of long-term breeding effort, and inadequate plant adaptation to changing climate. Studies examining the impacts of global climate change on bromegrass species show that global warming, heat stress, and drought have negative effects on forage yield. A number of useful physiological traits have been identified for genetic improvement to minimize yield loss. Available reports suggest that high-throughput imaging phenotyping techniques, including visual and infrared thermal imaging, imaging hyperspectral spectroscopy, and imaging chlorophyll fluorescence, are capable of capturing images of morphological, physiological, and biochemical traits related to plant growth, yield, and adaptation to abiotic stresses at different scales of organization. The more complex traits such as photosynthetic radiation-use efficiency, water-use efficiency, and nitrogen-use efficiency can be effectively assessed by utilizing combinations of imaging hyperspectral spectroscopy, infrared thermal imaging, and imaging chlorophyll fluorescence techniques in a breeding program. Genomic selection has been applied in the breeding of forage species and the applications show its potential in high ploidy, outcrossing species like bromegrass to improve the accuracy of parental selection and improve genetic gain. Together, these new technologies hold promise for improved genetic gain and wide adaptation in future bromegrass breeding.

Genotyping-by-Sequencing Enhances Genetic Diversity Analysis of Crested Wheatgrass [Agropyron cristatum (L.) Gaertn.].

Molecular characterization of unsequenced plant species with complex genomes is now possible by genotyping-by-sequencing (GBS) using recent next generation sequencing technologies. This study represents the first use of GBS application to sample genome-wide variants of crested wheatgrass [Agropyron cristatum (L.) Gaertn.] and assess the genetic diversity present in 192 genotypes from 12 tetraploid lines. Bioinformatic analysis identified 45,507 single nucleotide polymorphism (SNP) markers in this outcrossing grass species. The model-based Bayesian analysis revealed four major clusters of the samples assayed. The diversity analysis revealed 15.8% of SNP variation residing among the 12 lines, and 12.1% SNP variation present among four genetic clusters identified by the Bayesian analysis. The principal coordinates analysis and dendrogram were able to distinguish four lines of Asian origin from Canadian cultivars and breeding lines. These results serve as a valuable resource for understanding genetic variability, and will aid in the genetic improvement of this outcrossing polyploid grass species for forage production. These findings illustrate the potential of GBS application in the characterization of non-model polyploid plants with complex genomes.

Genotyping-by-sequencing data of 272 crested wheatgrass (Agropyron cristatum) genotypes.

Crested wheatgrass [Agropyron cristatum L. (Gaertn.)] is an important cool-season forage grass widely used for early spring grazing. However, the genomic resources for this non-model plant are still lacking. Our goal was to generate the first set of next generation sequencing data using the genotyping-by-sequencing technique. A total of 272 crested wheatgrass plants representing seven breeding lines, five cultivars and five geographically diverse accessions were sequenced with an Illumina MiSeq instrument. These sequence datasets were processed using different bioinformatics tools to generate contigs for diploid and tetraploid plants and SNPs for diploid plants. Together, these genomic resources form a fundamental basis for genomic studies of crested wheatgrass and other wheatgrass species. The raw reads were deposited into Sequence Read Archive (SRA) database under NCBI accession SRP115373 (https://www.ncbi.nlm.nih.gov/sra?term=SRP115373) and the supplementary datasets are accessible in Figshare (10.6084/m9.figshare.5345092).

RNA-Seq Analysis of Plant Maturity in Crested Wheatgrass (Agropyron cristatum L.).

Crested wheatgrass (Agropyron cristatum L.) breeding programs aim to develop later maturing cultivars for extending early spring grazing in Western Canada. Plant maturity is a complex genetic trait, and little is known about genes associated with late maturity in this species. An attempt was made using RNA-Seq to profile the transcriptome of crested wheatgrass maturity and to analyze differentially expressed genes (DEGs) between early and late maturing lines. Three cDNA libraries for each line were generated by sampling leaves at the stem elongation stage, spikes at the boot and anthesis stages. A total of 75,218,230 and 74,015,092 clean sequence reads were obtained for early and late maturing lines, respectively. De novo assembly of all sequence reads generated 401,587 transcripts with a mean length of 546 bp and N50 length of 691 bp. Out of 13,133 DEGs detected, 22, 17, and eight flowering related DEGs were identified for the three stages, respectively. Twelve DEGs, including nine flowering related DEGs at the stem elongation stage were further confirmed by qRT-PCR. The analysis of homologous genes of the photoperiod pathway revealed their lower expression in the late maturing line at the stem elongation stage, suggesting that their differential expression contributed to late maturity in crested wheatgrass.

RNA-Seq analysis of gene expression for floral development in crested wheatgrass (Agropyron cristatum L.).

Crested wheatgrass [Agropyron cristatum L. (Gaertn.)] is widely used for early spring grazing in western Canada and the development of late maturing cultivars which maintain forage quality for a longer period is desired. However, it is difficult to manipulate the timing of floral transition, as little is known about molecular mechanism of plant maturity in this species. In this study, RNA-Seq and differential gene expression analysis were performed to investigate gene expression for floral initiation and development in crested wheatgrass. Three cDNA libraries were generated and sequenced to represent three successive growth stages by sampling leaves at the stem elongation stage, spikes at boot and anthesis stages. The sequencing generated 25,568,846; 25,144,688 and 25,714,194 qualified Illumina reads for the three successive stages, respectively. De novo assembly of all the reads generated 311,671 transcripts with a mean length of 487 bp, and 152,849 genes with an average sequence length of 669 bp. A total of 48,574 (31.8%) and 105,222 (68.8%) genes were annotated in the Swiss-Prot and NCBI non-redundant (nr) protein databases, respectively. Based on the Kyoto Encyclopedia of Genes and Genome (KEGG) pathway database, 9,723 annotated sequences were mapped onto 298 pathways, including plant circadian clock pathway. Specifically, 113 flowering time-associated genes, 123 MADS-box genes and 22 CONSTANS-LIKE (COL) genes were identified. A COL homolog DN52048-c0-g4 which was clustered with the flowering time genes AtCO and OsHd1 in Arabidopsis (Arabidopsis thaliana L.) and rice (Oryza sativa L.), respectively, showed specific expression in leaves and could be a CONSTANS (CO) candidate gene. Taken together, this study has generated a new set of genomic resources for identifying and characterizing genes and pathways involved in floral transition and development in crested wheatgrass. These findings are significant for further understanding of the molecular basis for late maturity in this grass species.

Fermentation, degradation and microbial nitrogen partitioning for three forage colour phenotypes within anthocyanidin-accumulating Lc-alfalfa progeny.

BACKGROUND: Alfalfa has the disadvantage of having a rapid initial rate of protein degradation, which results in pasture bloat, low efficiency of protein utilisation and excessive nitrogen (N) pollution into the environment for cattle. Introducing a gene that stimulates the accumulation of monomeric/polymeric anthocyanidins might reduce the ruminal protein degradation rate (by fixing protein and/or direct interaction with microbes) and additionally reduce methane emission. The objectives of this study were to evaluate in vitro fermentation, degradation and microbial N partitioning of three forage colour phenotypes (green, light purple-green (LPG) and purple-green (PG)) within newly developed Lc-progeny and to compare them with those of parental green non-transgenic (NT) alfalfa. RESULTS: PG-Lc accumulated more anthocyanidin compared with Green-Lc (P < 0.05), with LPG-Lc intermediate. Volatile fatty acids and potentially degradable dry matter (DM) and N were similar among the four phenotypes. Gas, methane and ammonia accumulation rates were slower for the two purple-Lc phenotypes compared with NT-alfalfa (P < 0.05), while Green-Lc was intermediate. Effective degradable DM and N were lower in the three Lc-phenotypes (P < 0.05) compared with NT-alfalfa. Anthocyanidin concentration was negatively correlated (P < 0.05) with gas and methane production rates and effective degradability of DM and N. CONCLUSION: The Lc-alfalfa phenotypes accumulated anthocyanidin. Fermentation and degradation parameters indicated a reduced rate of fermentation and effective degradability for both purple anthocyanidin-accumulating Lc-alfalfa phenotypes compared with NT-alfalfa.

Interference of condensed tannin in lignin analyses of dry bean and forage crops.

Legumes with high concentrations of condensed tannin (pinto bean [Phaseolus vulgaris L.], sainfoin [Onobrychis viciifolia Scop.], and big trefoil [Lotus uliginosus Hoff.]), were compared to a selection of forages, with low or zero condensed tannin (smooth bromegrass [ Bromus inermis Leyss], Lotus japonicus [Regel] K. Larsen, and alfalfa [Medicago sativa L.]), using four methods to estimate fiber or lignin. Protocols were validated by using semipurified condensed tannin polymers in adulteration assays that tested low-lignin tissue with polyphenolic-enriched samples. The effect on lignin assay methods by condensed tannin concentration was interpreted using a multivariate analysis. There was an overestimation of fiber or lignin in the presence of condensed tannin in the acid detergent fiber (ADF) and Klason lignin (KL) assays compared to that in the thioglycolic acid (TGA) and acid detergent lignin (ADL) methods. Sulfite reagents (present in TGA lignin method) or sequential acidic digests at high temperatures (ADF followed by ADL) were required to eliminate condensed tannin. The ADF (alone) and KL protocols are not recommended to screen nonwoody plants, such as forages, where condensed tannin has accumulated in the tissue.

Genetic variation and correlation of agronomic traits in meadow bromegrass (Bromus riparius Rehm) clones

Meadow bromegrass (Bromus riparius Rehm.) é uma gramínea forrageira para pastejo recentemente introduzida no oeste do Canadá. Sua produção de matéria verde aliada a uma rebrota rápida a transformou numa forrageira bastante utilizada no pastejo de gado de corte naquela região. Existe pouca informação sobre o comportamento da espécie quando submetida a trabalhos de melhoramento genético. O principal objetivo deste estudo foi, portanto, estimar a variabilidade genética total, a herdabilidade no sentido amplo, bem como as correlações fenotípicas e genotípicas. Caracteres agronômicos de quarenta e quatro clones de Meadow bromegrass foram avaliadas em dois locais. A variação genética, para produção matéria seca, produção de sementes, índice de fertilidade, índice de colheita, altura das plantas, diâmetro das plantas, proteina bruta, fibra em detergente ácido e fibra em detergente neutro, foi estatisticamente significativa. As estimativas de herdabilidade foram superiores a 50 por cento, para todos os caracteres avaliados. Estas estimativas superaram seu erro padrão em, no mínimo, 3,5 vezes. As correlações genotípicas e fenotípicas entres todos os caracteres também foram estimadas, sendo as mesmas similares tanto em sinal como em magnitudes. Estas correlações demostram que é possível melhorar simultaneamente a produção de forragem e a de sementes. Estes resultados indicam que é possível desenvolver cultivares com alta produtividade e altos teores de proteina bruta, associados a baixo teores de fibra.

Expression of anthocyanins and proanthocyanidins after transformation of alfalfa with maize Lc.

Three anthocyanin regulatory genes of maize (Zea mays; Lc, B-Peru, and C1) were introduced into alfalfa (Medicago sativa) in a strategy designed to stimulate the flavonoid pathway and alter the composition of flavonoids produced in forage. Lc constructs included a full-length gene and a gene with a shortened 5'-untranslated region. Lc RNA was strongly expressed in Lc transgenic alfalfa foliage, but accumulation of red-purple anthocyanin was observed only under conditions of high light intensity or low temperature. These stress conditions induced chalcone synthase and flavanone 3-hydroxylase expression in Lc transgenic alfalfa foliage compared with non-transformed plants. Genotypes containing the Lc transgene construct with a full-length 5'-untranslated region responded more quickly to stress conditions and with a more extreme phenotype. High-performance liquid chromatography analysis of field-grown tissue indicated that flavone content was reduced in forage of the Lc transgenic plants. Leucocyanidin reductase, the enzyme that controls entry of metabolites into the proanthocyanidin pathway, was activated both in foliage and in developing seeds of the Lc transgenic alfalfa genotypes. Proanthocyanidin polymer was accumulated in the forage, but (+)-catechin monomers were not detected. B-Peru transgenic and C1 transgenic populations displayed no visible phenotypic changes, although these transgenes were expressed at detectable levels. These results support the emerging picture of Lc transgene-specific patterns of expression in different recipient species. These results demonstrate that proanthocyanidin biosynthesis can be stimulated in alfalfa forage using an myc-like transgene, and they pave the way for the development of high quality, bloat-safe cultivars with ruminal protein bypass.