A genetic linkage map of tetraploid orchardgrass (Dactylis glomerata L.) and quantitative trait loci for heading date

Metabolomic Study of Dactylis glomerata Growing on Aeolian Archipelago (Italy)

The Aeolian Islands (Italy) are a volcanic archipelago in the Tyrrhenian Sea comprising seven main islands, among which are two active volcanoes. The peculiar geological features and the wide variety of environments and soils have an important impact on native plants, and in particular, the Aeolian populations of Dactylis glomerata (a perennial cool-season bunchgrass) exhibit remarkable phenotypic variability. Considering that environmental drivers also strongly affect the production of plant metabolites, this work aimed at comparing the metabolomic profiles of D. glomerata (leaves) harvested at different altitudes on four islands of the Aeolian archipelago, namely: Lipari, Vulcano, Stromboli and Panarea. Samples were analyzed by ¹H NMR profiling, and data were treated by PCA. Samples collected on Stromboli were very different from each other and from the samples collected in the other islands. Through an Orthogonal Partial Least Squares (OPLS) model, using altitude as the y variable, it emerged that the concentration of proline, glycine betaine, sucrose, glucose and chlorogenic acid of D. glomerata growing on Stromboli decreased at increasing altitude. Conversely, increasing altitude was associated with an increment in valine, asparagine, fumaric acid and phenylalanine.

Identification of Flowering Regulatory Networks and Hub Genes Expressed in the Leaves of Elymus sibiricus L. Using Comparative Transcriptome Analysis

Flowering is a significant stage from vegetative growth to reproductive growth in higher plants, which impacts the biomass and seed yield. To reveal the flowering time variations and identify the flowering regulatory networks and hub genes in Elymus sibiricus, we measured the booting, heading, and flowering times of 66 E. sibiricus accessions. The booting, heading, and flowering times varied from 136 to 188, 142 to 194, and 148 to 201 days, respectively. The difference in flowering time between the earliest- and the last-flowering accessions was 53 days. Furthermore, transcriptome analyses were performed at the three developmental stages of six accessions with contrasting flowering times. A total of 3,526 differentially expressed genes (DEGs) were predicted and 72 candidate genes were identified, including transcription factors, known flowering genes, and plant hormone-related genes. Among them, four candidate genes (LATE, GA2OX6, FAR3, and MFT1) were significantly upregulated in late-flowering accessions. LIMYB, PEX19, GWD3, BOR7, PMEI28, LRR, and AIRP2 were identified as hub genes in the turquoise and blue modules which were related to the development time of flowering by weighted gene co-expression network analysis (WGCNA). A single-nucleotide polymorphism (SNP) of LIMYB found by multiple sequence alignment may cause late flowering. The expression pattern of flowering candidate genes was verified in eight flowering promoters (CRY, COL, FPF1, Hd3, GID1, FLK, VIN3, and FPA) and four flowering suppressors (CCA1, ELF3, Ghd7, and COL4) under drought and salt stress by qRT-PCR. The results suggested that drought and salt stress activated the flowering regulation pathways to some extent. The findings of the present study lay a foundation for the functional verification of flowering genes and breeding of new varieties of early- and late-flowering E. sibiricus.

Development of EST-SSR Markers Linked to Flowering Candidate Genes in Elymus sibiricus L. Based on RNA Sequencing

Elymus sibiricus L. is an important cold-season grass with excellent cold and drought tolerance, good palatability, and nutrition. Flowering time is a key trait that affects forage and seed yield. Development of EST-SSR (expressed sequence tag simple sequence repeat) markers based on flowering genes contributes to the improvement of flowering traits. In the study, we detected 155 candidate genes related to flowering traits from 10,591 unigenes via transcriptome sequencing in early- and late-flowering genotypes. These candidate genes were mainly involved in the photoperiodic pathway, vernalization pathway, central integrator, and gibberellin pathway. A total of 125 candidate gene-based EST-SSRs were developed. Further, 15 polymorphic EST-SSRs closely associated to 13 candidate genes were used for genetic diversity and population structure analysis among 20 E. sibiricus accessions, including two contrasting panels (early-flowering and late-flowering). Among them, primer 28366, designed from heading date 3a (HD3a), effectively distinguished early- and late-flowering genotypes using a specifically amplified band of 175 bp. The polymorphic information content (PIC) value ranged from 0.12 to 0.48, with an average of 0.25. The unweighted pair group method analysis (UPGMA) cluster and structure analysis showed that the 20 E. sibiricus genotypes with similar flowering times tended to group together. These newly developed EST-SSR markers have the potential to be used for molecular markers assisted selection and germplasm evaluation of flowering traits in E. sibiricus.

Evaluation of genetic integrity of pearl millet seeds during aging by genomic-SSR markers

Seed is an important way to store germplasm resources but its genetic integrity will decrease during long-term preservation. So, it's essential to update seeds according to the aging level of different species. Pearl millet [Cenchrus americanus (L.) Morrone syn., Pennisetum glaucum (L.) R. Br.] is a crucial forage grass, biofuel plant and important crops in the world bringing huge economic and ecological benefits. However, there is no report about the impact of aging on genetic integrity of its seeds. In this study, four genetic diversity indexes (the percentage of polymorphic bands, PPB; the effective number of alleles, Ne; the Nei's gene diversity index, H; the Shannon's information index, I) and 20 pairs of genomic-SSR primers were used to certify the optimal sample volume of pearl millet for molecular study and found that the best sample volume was 60. After the artificial aging test, the germination rate and four genetic diversity parameters (the number of alleles, Na; Ne; H; I) were used to evaluate the change of genetic integrity at different aging levels. The results showed that the germination rate and these four genetic diversity parameters declined with the increase of aging levels. Furthermore, when the germination rate of pearl millet seeds went down to 68.23%, a significant difference in genetic integrity was observed with unaged seeds. In conclusion, the optimal sample size of pearl millet was 60 and the critical point of germination rate to renew germplasm resources was 68.23% and these finds might contribute to the scientific study and the safe conservation of pearl millet.

Identification and Distribution of NBS-Encoding Resistance Genes of Dactylis glomerata L. and Its Expression Under Abiotic and Biotic Stress

Orchardgrass (Dactylis glomerata L.) is drought resistant and tolerant to barren landscapes, making it one of the most important forages for animal husbandry, as well as ecological restoration of rocky landscapes that are undergoing desertification. However, orchardgrass is susceptible to rust, which can significantly reduce its yield and quality. Therefore, understanding the genes that underlie resistance against rust in orchardgrass is critical. The evolution, cloning of plant disease resistance genes, and the analysis of pathogenic bacteria induced expression patterns are important contents in the study of interaction between microorganisms and plants. Genes with nucleotide binding site (NBS) structure are disease-resistant genes ubiquitous in plants and play an important role in plant attacks against various pathogens. Using sequence analysis and re-annotation, we identified 413 NBS resistance genes in orchardgrass. Similar to previous studies, NBS resistance genes containing TIR (toll/interleukin-1 receptor) domain were not found in orchardgrass. The NBS resistance genes can be divided into four types: NBS (up to 264 homologous genes, accounting for 64% of the total number of NBS genes in orchardgrass), NBS-LRR, CC-NBS, and CC-NBS-LRR (minimum of 26 homologous genes, only 6% of the total number of NBS genes in orchardgrass). These 413 NBS resistance genes were unevenly distributed across seven chromosomes where chromosome 5 had up to 99 NBS resistance genes. There were 224 (54%) NBS resistance genes expressed in different tissues (roots, stems, leaves, flowers, and spikes), and we did not detect expression for 45 genes (11%). The remaining 145 (35%) were expressed in some tissues. And we found that 11 NBS resistance genes were differentially expressed under waterlogging stress, 5 NBS resistance genes were differentially expressed under waterlogging and drought stress, and 1 NBS resistance was is differentially expressed under waterlogging and heat stress. Most importantly, we found that 65 NBS resistance genes were significantly expressed in different control groups. On the 7th day of inoculation, 23 NBS resistance genes were differentially expressed in high resistance materials alone, of which 7 NBS resistance genes regulate the "plant-pathogen interaction" pathway by encoding RPM1. At the same time, 2 NBS resistance genes that were differentially expressed in the high resistance material after inoculation were also differentially expressed in abiotic stress. In summary, the NBS resistance gene plays a crucial role in the resistance of orchardgrass to rust.

Construction of the first high-density genetic linkage map and identification of seed yield-related QTLs and candidate genes in Elymus sibiricus, an important forage grass in Qinghai-Tibet Plateau

Background

Elymus sibiricus is an ecologically and economically important perennial, self-pollinated, and allotetraploid (StStHH) grass, widely used for forage production and animal husbandry in Western and Northern China. However, it has low seed yield mainly caused by seed shattering, which makes seed production difficult for this species. The goals of this study were to construct the high-density genetic linkage map, and to identify QTLs and candidate genes for seed-yield related traits.

Results

An F₂ mapping population of 200 individuals was developed from a cross between single genotype from “Y1005” and “ZhN06”. Specific-locus amplified fragment sequencing (SLAF-seq) was applied to construct the first genetic linkage map. The final genetic map included 1971 markers on the 14 linkage groups (LGs) and was 1866.35 cM in total. The length of each linkage group varied from 87.67 cM (LG7) to 183.45 cM (LG1), with an average distance of 1.66 cM between adjacent markers. The marker sequences of E. sibiricus were compared to two grass genomes and showed 1556 (79%) markers mapped to wheat, 1380 (70%) to barley. Phenotypic data of eight seed-related traits (2016–2018) were used for QTL identification. A total of 29 QTLs were detected for eight seed-related traits on 14 linkage groups, of which 16 QTLs could be consistently detected for two or three years. A total of 6 QTLs were associated with seed shattering. Based on annotation with wheat and barley genome and transcriptome data of abscission zone in E. sibiricus, we identified 30 candidate genes for seed shattering, of which 15, 7, 6 and 2 genes were involved in plant hormone signal transcription, transcription factor, hydrolase activity and lignin biosynthetic pathway, respectively.

Conclusion

This study constructed the first high-density genetic linkage map and identified QTLs and candidate genes for seed-related traits in E. sibiricus. Results of this study will not only serve as genome-wide resources for gene/QTL fine mapping, but also provide a genetic framework for anchoring sequence scaffolds on chromosomes in future genome sequence assembly of E. sibiricus.

Association of candidate genes with heading date in a diverse Dactylis glomerata population

Flowering occurs in response to cues from both temperature and photoperiod elicitors in cool-season, long-day forage grasses, and genes involved in sensing the elicitors and inducing downstream flowering responses have been associated with heading date and flowering time in perennial forage grasses as well as cereal grasses. In this study we test for association between orchardgrass (Dactylis glomerata L.) heading date and polymorphisms in the CONSTANS (DgCO1), FLOWERING TIME (DgFT1), a VRN1 like MADS-box (DgMADS), and PHOTOPERIOD (DgPPD1-like) containing genes. A diverse population of 150 genotypes was measured for heading date across three years, genotyped, and candidate genes sequenced. Although pairwise population kinship values were generally low, the genotypes fit into a two-group structure model. Linkage disequilibrium decayed rapidly, reaching r² levels below 0.2 within the 500bp of each gene. SNPs significantly associated with heading date were detected in equal-dose and tetraploid dosage models. The DgCO1 gene had the most significant polymorphisms and those with the largest effects, while DgMADS had several significant polymorphisms in its first intron with smaller effects. These polymorphisms can be used for further validation, selection, and development of breeding lines of orchardgrass.

Construction of high-density genetic linkage map and identification of flowering-time QTLs in orchardgrass using SSRs and SLAF-seq

Orchardgrass (Dactylis glomerata L.) is one of the most economically important perennial, cool-season forage species grown and pastured worldwide. High-density genetic linkage mapping is a valuable and effective method for exploring complex quantitative traits. In this study, we developed 447,177 markers based on SLAF-seq and used them to perform a comparative genomics analysis. Perennial ryegrass sequences were the most similar (5.02%) to orchardgrass sequences. A high-density linkage map of orchardgrass was constructed using 2,467 SLAF markers and 43 SSRs, which were distributed on seven linkage groups spanning 715.77 cM. The average distance between adjacent markers was 0.37 cM. Based on phenotyping in four environments, 11 potentially significant quantitative trait loci (QTLs) for two target traits–heading date (HD) and flowering time (FT)–were identified and positioned on linkage groups LG1, LG3, and LG5. Significant QTLs explained 8.20–27.00% of the total phenotypic variation, with the LOD ranging from 3.85–12.21. Marker167780 and Marker139469 were associated with FT and HD at the same location (Ya’an) over two different years. The utility of SLAF markers for rapid generation of genetic maps and QTL analysis has been demonstrated for heading date and flowering time in a global forage grass.

Assessing and Broadening Genetic Diversity of Elymus sibiricus Germplasm for the Improvement of Seed Shattering

Siberian wild rye (Elymus sibiricus L.) is an important native grass in the Qinghai-Tibet Plateau of China. It is difficult to grow for commercial seed production, since seed shattering causes yield losses during harvest. Assessing the genetic diversity and relationships among germplasm from its primary distribution area contributes to evaluating the potential for its utilization as a gene pool to improve the desired agronomic traits. In the study, 40 EST-SSR primers were used to assess the genetic diversity and population structure of 36 E. sibiricus accessions with variation of seed shattering. A total of 380 bands were generated, with an average of 9.5 bands per primer. The polymorphic information content (PIC) ranged from 0.23 to 0.50. The percentage of polymorphic bands (P) for the species was 87.11%, suggesting a high degree of genetic diversity. Based on population structure analysis, four groups were formed, similar to results of principal coordinate analysis (PCoA). The molecular variance analysis (AMOVA) revealed the majority of genetic variation occurred within geographical regions (83.40%). Two genotypes from Y1005 and ZhN06 were used to generate seven F₁ hybrids. The molecular and morphological diversity analysis of F₁ population revealed rich genetic variation and high level of seed shattering variation in F₁ population, resulting in significant improvement of the genetic base and desired agronomic traits.

Combining Drought Survival via Summer Dormancy and Annual Biomass Productivity in Dactylis glomerata L

Under Mediterranean climates, the best strategy to produce rain-fed fodder crops is to develop perennial drought resistant varieties. Summer dormancy present in native germplasm has been shown to confer a high level of survival under severe drought. Nevertheless it has also been shown to be negatively correlated with annual biomass productivity. The aim of this study was to analyze the correlations between summer dormancy and annual biomass productivity related traits and to identify quantitative trait loci (QTL) for these traits in a progeny of a summer dormant cocksfoot parent (Kasbah) and a summer active parent (Medly). A total of 283 offspring and the parents were phenotyped for summer dormancy, plant growth rate (PGR) and heading date in Morocco and for maximum leaf elongation rate (LERm) in France. The individuals were genotyped with a total of 325 markers including 59 AFLP, 64 SSR, and 202 DArT markers. The offspring exhibited a large quantitative variation for all measured traits. Summer dormancy showed a negative correlation with both PGR (-0.34 p < 0.005) and LERm (-0.27 p < 0.005). However, genotypes with both a high level of summer dormancy and a high level of PGR were detected in the progeny. One genetic map per parent was built with a total length of 377 and 423 cM for Kasbah and Medly, respectively. Both different and co-localized QTL for summer dormancy and PGR were identified. These results demonstrate that it should be possible to create summer dormant cocksfoot varieties with a high annual biomass productivity.