Glycoside Hydrolase 3 (GH3), a member of the Auxin-responsive gene family, is involved in plant growth, the plant developmental process, and various stress responses. The GH3 gene family has been well-studied in Arabidopsis thaliana and Zea mays. However, the evolution of the GH3 gene family in Oryza species remains unknown and the function of the GH3 gene family in Oryza sativa is not well-documented. Here, a systematic analysis was performed in six Oryza species/subspecies, including four wild rice species and two cultivated rice subspecies. A total of 13, 13, 13, 13, 12, and 12 members were identified in O. sativa ssp. japonica, O. sativa ssp. indica, Oryza rufipogon, Oryza nivara, Oryza punctata, and Oryza glumaepatula, respectively. Gene duplication events, structural features, conserved motifs, a phylogenetic analysis, chromosome locations, and Ka/Ks ratios of this important family were found to be strictly conservative across these six Oryza species/subspecies, suggesting that the expansion of the GH3 gene family in Oryza species might be attributed to duplication events, and this expansion could occur in the common ancestor of Oryza species, even in common ancestor of rice tribe (Oryzeae) (23.07~31.01 Mya). The RNA-seq results of different tissues displayed that OsGH3 genes had significantly different expression profiles. Remarkably, the qRT-PCR result after NaCl treatment indicated that the majority of OsGH3 genes play important roles in salinity stress, especially OsGH3-2 and OsGH3-8. This study provides important insights into the evolution of the GH3 gene family in Oryza species and will assist with further investigation of OsGH3 genes’ functions under salinity stress.

Rice (Oryza sativa L.) is one of the most important staple food crops worldwide, while its growth and productivity are threatened by various abiotic stresses, especially salt stress. Unraveling how rice adapts to salt stress at the transcription level is vital. It can provide valuable information on enhancing the salt stress tolerance performance of rice via genetic engineering technologies. Here, we conducted a meta-analysis of different rice genotypes at the seedling stage based on 96 public microarray datasets, aiming to identify the key salt-responsive genes and understand the molecular response mechanism of rice under salt stress. In total, 5559 genes were identified to be differentially expressed genes (DEGs) under salt stress, and 3210 DEGs were identified during the recovery process. The Gene Ontology (GO) enrichment results revealed that the salt-response mechanisms of shoots and roots were different. A close-knit signaling network, consisting of the Ca²⁺ signal transduction pathway, the mitogen-activated protein kinase (MAPK) cascade, multiple hormone signals, transcription factors (TFs), transcriptional regulators (TRs), protein kinases (PKs), and other crucial functional proteins, plays an essential role in rice salt stress response. In this study, many unreported salt-responsive genes were found. Besides this, MapMan results suggested that TNG67 can shift to the fermentation pathway to produce energy under salt stress and may enhance the Calvin cycle to repair a damaged photosystem during the recovery stage. Taken together, these findings provide novel insights into the salt stress molecular response and introduce numerous candidate genes for rice salt stress tolerance breeding.

Plant cell walls play a fundamental role in several ways, providing structural support for cells, resistance against pathogens and facilitating the communication between cells. The glycosyltransferase family 8 (GT8) is involved in the formation of the plant cell wall. However, the evolutionary relationship and the functional differentiation of this important gene family remain obscure in Gramineae crop genomes. In the present investigation, we identified 269 GT8 genes in the seven Gramineae representative crop genomes, namely, 33 in Hordeum vulgare, 37 in Brachypodium distachyon, 40 in Oryza sativa ssp. japonica, 41 in Oryza rufipogon, 36 in Setaria italica, 37 in Sorghum bicolor, and 45 in Zea mays. Phylogenetic analysis suggested that all identified GT8 proteins belonged to seven subfamilies: galacturonosyltransferase (GAUT), galacturonosyltransferase-like (GATL), GATL-related (GATR), galactinol synthase (GolS), and plant glycogenin-like starch initiation proteins A (PGSIP-A), PGSIP-B, and PGSIP-C. We estimated that the GAUT subfamily might be further divided into four subgroups (I–IV) due to differentiation of gene structures and expression patterns. Our orthogroup analysis identified 22 orthogroups with different sizes. Of these orthogroups, several orthogroups were lost in some species, such as S. italica and Z. mays. Moreover, lots of duplicate pairs and collinear pairs were discovered among these species. These results indicated that multiple duplication modes led to the expansion of this important gene family and unequal loss of orthogroups and subfamilies might have happened during the evolutionary process. RNA-seq, microarray analysis, and qRT-PCR analyses indicated that GT8 genes are critical for plant growth and development, and for stresses responses. We found that OsGolS1 was significantly up-regulated under salt stress, while OsGAUT21, OsGATL2, and OsGATL5 had remarkable up-regulation under cold stress. The current study highlighted the expansion and evolutionary patterns of the GT8 gene family in these seven Gramineae crop genomes and provided potential candidate genes for future salt- and cold- resistant molecular breeding studies in O. sativa.

Saline–alkaline stress is one of the major damages that severely affects rice (Oryza sativa L.) growth and grain yield; however, the mechanism of the tolerance remains largely unknown in rice. Herein, we comparatively investigated the transcriptome and metabolome of two contrasting rice subspecies genotypes, Luohui 9 (abbreviation for Chao2R under study, O. sativa ssp. indica, saline–alkaline-sensitive) and RPY geng (O. sativa ssp. japonica, saline–alkaline-tolerant), to identify the main pathways and important factors related to saline–alkaline tolerance. Transcriptome analysis showed that 68 genes involved in fatty acid, amino acid (such as phenylalanine and tryptophan), phenylpropanoid biosynthesis, energy metabolism (such as Glycolysis and TCA cycle), as well as signal transduction (such as hormone and MAPK signaling) were identified to be specifically upregulated in RPY geng under saline–alkaline conditions, implying that a series of cascade changes from these genes promotes saline–alkaline stress tolerance. The transcriptome changes observed in RPY geng were in high accordance with the specifically accumulation of metabolites, consisting mainly of 14 phenolic acids, 8 alkaloids, and 19 lipids based on the combination analysis of transcriptome and metabolome. Moreover, some genes involved in signal transduction as hub genes, such as PR5, FLS2, BRI1, and NAC, may participate in the saline–alkaline stress response of RPY geng by modulating key genes involved in fatty acid, phenylpropanoid biosynthesis, amino acid metabolism, and glycolysis metabolic pathways based on the gene co-expression network analysis. The present research results not only provide important insights for understanding the mechanism underlying of rice saline–alkaline tolerance at the transcriptome and metabolome levels but also provide key candidate target genes for further enhancing rice saline–alkaline stress tolerance.

The plant hormones gibberellins (GAs) regulate plant growth and development and are closely related to the yield of cash crops. The GA oxidases (GAoxs), including the GA2ox, GA3ox, and GA20ox subfamilies, play pivotal roles in GAs’ biosynthesis and metabolism, but their classification and evolutionary pattern in Gramineae crops remain unclear. We thus conducted a comparative genomic study of GAox genes in six Gramineae representative crops, namely, Setaria italica (Si), Zea mays (Zm), Sorghum bicolor (Sb), Hordeum vulgare (Hv), Brachypodium distachyon (Bd), and Oryza sativa (Os). A total of 105 GAox genes were identified in these six crop genomes, belonging to the C₁₉-GA2ox, C₂₀-GA2ox, GA3ox, and GA20ox subfamilies. Based on orthogroup (OG) analysis, GAox genes were divided into nine OGs and the number of GAox genes in each of the OGs was similar among all tested crops, which indicated that GAox genes may have completed their family differentiations before the species differentiations of the tested species. The motif composition of GAox proteins showed that motifs 1, 2, 4, and 5, forming the 2OG-FeII_Oxy domain, were conserved in all identified GAox protein sequences, while motifs 11, 14, and 15 existed specifically in the GA20ox, C₁₉-GA2ox, and C₂₀-GA2ox protein sequences. Subsequently, the results of gene duplication events suggested that GAox genes mainly expanded in the form of WGD/SD and underwent purification selection and that maize had more GAox genes than other species due to its recent duplication events. The cis-acting elements analysis indicated that GAox genes may respond to growth and development, stress, hormones, and light signals. Moreover, the expression profiles of rice and maize showed that GAox genes were predominantly expressed in the panicles of the above two plants and the expression of several GAox genes was significantly induced by salt or cold stresses. In conclusion, our results provided further insight into GAox genes’ evolutionary differences among six representative Gramineae and highlighted GAox genes that may play a role in abiotic stress.

Adipose tissue composition contributes greatly to the quality and nutritional value of meat. Transcriptomic and lipidomic techniques were used to investigate the molecular mechanisms of the differences in fat deposition in Ningxiang pigs, Berkshires and F₁ offspring. Transcriptomic analysis identified 680, 592, and 380 DEGs in comparisons of Ningxiang pigs vs. Berkshires, Berkshires vs. F₁ offspring, and Ningxiang pigs vs. F₁ offspring. The lipidomic analysis screened 423, 252, and 50 SCLs in comparisons of Ningxiang pigs vs. Berkshires, Berkshires vs. F₁ offspring, and Ningxiang pigs vs. F₁ offspring. Lycine, serine, and the threonine metabolism pathway, fatty acid biosynthesis and metabolism-related pathways were significantly enriched in comparisons of Berkshires vs. Ningxiang pigs and Berkshires vs. F₁ offspring. The DEGs (PHGDH, LOC110256000) and the SCLs (phosphatidylserines) may have a great impact on the glycine, serine, and the threonine metabolism pathway. Moreover, the DEGs (FASN, ACACA, CBR4, SCD, ELOV6, HACD2, CYP3A46, CYP2B22, GPX1, and GPX3) and the SCLs (palmitoleic acid, linoleic acid, arachidonic acid, and icosadienoic acid) play important roles in the fatty acid biosynthesis and metabolism of fatty acids. Thus, the difference in fat deposition among Ningxiang pig, Berkshires, and F₁ offspring may be caused by differences in the expression patterns of key genes in multiple enriched KEGG pathways. This research revealed multiple lipids that are potentially available biological indicators and screened key genes that are potential targets for molecular design breeding. The research also explored the molecular mechanisms of the difference in fat deposition among Ningxiang pig, Berkshires, and F₁ pigs, and provided an insight into selection for backfat thickness and the fat composition of adipose tissue for future breeding strategies.