Background: Heat and drought are serious threats for crop growth and development. As the sixth largest cereal crop in the world, pearl millet can not only be used for food and forage but also as a source of bioenergy. Pearl millet is highly tolerant to heat and drought. Given this, it is considered an ideal crop to study plant stress tolerance and can be used to identify heat-resistant genes. Results: In this study, we used Pacbio sequencing data as a reference sequence to analyze the Illumina data of pearl millet that had been subjected to heat and drought stress for 48 h. By summarizing previous studies, we found 26,299 new genes and 63,090 new transcripts, and the number of gene annotations increased by 20.18%. We identified 2792 transcription factors and 1223 transcriptional regulators. There were 318 TFs and 149 TRs differentially expressed under heat stress, and 315 TFs and 128 TRs were differentially expressed under drought stress. We used RNA sequencing to identify 6920 genes and 6484 genes differentially expressed under heat stress and drought stress, respectively. Conclusions: Through Pacbio sequencing, we have identified more new genes and new transcripts. On the other hand, comparing the differentially expressed genes under heat tolerance with the DEGs under drought stress, we found that even in the same pathway, pearl millet responds with a different protein.

High-temperature stress negatively affects the growth and development of plants, and therefore threatens global agricultural safety. Cultivating stress-tolerant plants is the current objective of plant breeding programs. Pearl millet is a multi-purpose plant, commonly used as a forage but also an important food staple. This crop is very heat-resistant and has a higher net assimilation rate than corn under high-temperature stress. However, the response of heat resistant pearl millet has so far not been studied at the transcriptional level. In this study, transcriptome sequencing of pearl millet leaves exposed to different lengths of heat treatment (1 h, 48 h and 96 h) was conducted in order to investigate the molecular mechanisms of the heat stress response and to identify key genes related to heat stress. The results showed that the amount of heat stress-induced DEGs in leaves differs with the length of exposure to high temperatures. The highest value of DEGs (8286) was observed for the group exposed to heat stress for 96 h, while the other two treatments showed lower DEGs values of 4659 DEGs after 1 h exposure and 3981 DEGs after 48 h exposure to heat stress. The DEGs were mainly synthesized in protein folding pathways under high-temperature stress after 1 h exposure. Moreover, a large number of genes encoding ROS scavenging enzymes were activated under heat stress for 1 h and 48 h treatments. The flavonoid synthesis pathway of pearl millet was enriched after heat stress for 96 h. This study analyzed the transcription dynamics under short to long-term heat stress to provide a theoretical basis for the heat resistance response of pearl millet.

Orchardgrass (Dactylis glomerata L.), is a well-known perennial forage species; however, rust diseases have caused a noticeable reduction in the quality and production of orchardgrass. In this study, genetic diversity was assessed and the marker-trait associations for rust were examined using 18 EST-SSR and 21 SCoT markers in 75 orchardgrass accessions. A high level of genetic diversity was detected in orchardgrass with an average genetic diversity index of 0.369. For the EST-SSR and SCoT markers, 164 and 289 total bands were obtained, of which 148 (90.24%) and 272 (94.12%) were polymorphic, respectively. Results from an AMOVA analysis showed that more genetic variance existed within populations (87.57%) than among populations (12.43%). Using a parameter marker index, the efficiencies of the EST-SSR and SCoT markers were compared to show that SCoTs have higher marker efficiency (8.07) than EST-SSRs (4.82). The results of a UPGMA cluster analysis and a STRUCTURE analysis were both correlated with the geographic distribution of the orchardgrass accessions. Linkage disequilibrium analysis revealed an average r² of 0.1627 across all band pairs, indicating a high extent of linkage disequilibrium in the material. An association analysis between the rust trait and 410 bands from the EST-SSR and SCoT markers using TASSEL software revealed 20 band panels were associated with the rust trait in both 2011 and 2012. The 20 bands obtained from association analysis could be used in breeding programs for lineage selection to prevent great losses of orchardgrass caused by rust, and provide valuable information for further association mapping using this collection of orchardgrass.

To select the most stable reference genes in annual ryegrass (Lolium multiflorum), we studied annual ryegrass leaf tissues exposed to various abiotic stresses by qRT-PCR and selected 11 candidate reference genes, i.e., 18S rRNA, E2, GAPDH, eIF4A, HIS3, SAMDC, TBP-1, Unigene71, Unigene77, Unigene755, and Unigene14912. We then used GeNorm, NormFinder, and BestKeeper to analyze the expression stability of these 11 genes, and used RefFinder to comprehensively rank genes according to stability. Under different stress conditions, the most suitable reference genes for studies of leaf tissues of annual ryegrass were different. The expression of the eIF4A gene was the most stable under drought stress. Under saline-alkali stress, Unigene14912 has the highest expression stability. Under acidic aluminum stress, SAMDC expression stability was highest. Under heavy metal stress, Unigene71 expression had the highest stability. According to the software analyses, Unigene14912, HIS3, and eIF4A were the most suitable for analyses of abiotic stress in tissues of annual ryegrass. GAPDH was the least suitable reference gene. In conclusion, selecting appropriate reference genes under abiotic stress not only improves the accuracy of annual ryegrass gene expression analyses, but also provides a theoretical reference for the development of reference genes in plants of the genus Lolium.

Cis-regulatory elements (CREs) are essential for regulating gene expression, yet their evolutionary dynamics in plants remain elusive. Here we constructed a single-cell chromatin accessibility atlas for Oryza sativa from 103,911 nuclei representing 126 cell states across nine organs. Comparative genomics between O. sativa and 57,552 nuclei from four additional grass species (Zea mays, Sorghum bicolor, Panicum miliaceum and Urochloa fusca) revealed that chromatin accessibility conservation varies with cell-type specificity. Epidermal accessible chromatin regions in the leaf were less conserved compared to other cell types, indicating accelerated regulatory evolution in the L1-derived epidermal layer of O. sativa relative to other species. Conserved accessible chromatin regions overlapping the repressive histone modification H3K27me3 were identified as potentially silencer-like CREs, as deleting these regions led to up-regulation of gene expression. This study provides a comprehensive epigenomic resource for the rice community, demonstrating the utility of a comparative genomics approach that highlights the dynamics of plant cell-type-specific CRE evolution.

Drought stress is a global problem, and the lack of water is a key factor that leads to agricultural shortages. MicroRNAs play a crucial role in the plant drought stress response; however, the microRNAs and their targets involved in drought response have not been well elucidated. In the present study, we used Illumina platform (https://www.illumina.com/) and combined data from miRNA, RNA, and degradome sequencing to explore the drought- and organ-specific miRNAs in orchardgrass (Dactylis glomerata L.) leaf and root. We aimed to find potential miRNA–mRNA regulation patterns responding to drought conditions. In total, 519 (486 conserved and 33 novel) miRNAs were identified, of which, 41 miRNAs had significant differential expression among the comparisons (p < 0.05). We also identified 55,366 unigenes by RNA-Seq, where 12,535 unigenes were differently expressed. Finally, our degradome analysis revealed that 5950 transcripts were targeted by 487 miRNAs. A correlation analysis identified that miRNA ata-miR164c-3p and its target heat shock protein family A (HSP70) member 5 gene comp59407_c0 (BIPE3) may be essential in organ-specific plant drought stress response and/or adaptation in orchardgrass. Additionally, Gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) analyses found that “antigen processing and presentation” was the most enriched downregulated pathway in adaptation to drought conditions. Taken together, we explored the genes and miRNAs that may be involved in drought adaptation of orchardgrass and identified how they may be regulated. These results serve as a valuable genetic resource for future studies focusing on how plants adapted to drought conditions.