Genome-wide identification and classification of Lipoxygenase gene family and their roles in sorghum-aphid interaction

Integrated transcriptomic and pathway analyses of sorghum plants revealed the molecular mechanisms of host defense against aphids

Sugarcane aphid has emerged as a major pest of sorghum recently, and a few sorghum accessions were identified for resistance to this aphid so far. However, the molecular and genetic mechanisms underlying this resistance are still unclear. To understand these mechanisms, transcriptomics was conducted in resistant Tx2783 and susceptible BTx623 sorghum genotypes infested with sugarcane aphids. A principal component analysis revealed differences in the transcriptomic profiles of the two genotypes. The pathway analysis of the differentially expressed genes (DEGs) indicated the upregulation of a set of genes related to signal perception (nucleotide-binding, leucine-rich repeat proteins), signal transduction [mitogen-activated protein kinases signaling, salicylic acid (SA), and jasmonic acid (JA)], and plant defense (transcription factors, flavonoids, and terpenoids). The upregulation of the selected DEGs was verified by real-time quantitative PCR data analysis, performed on the resistant and susceptible genotypes. A phytohormone bioassay experiment showed a decrease in aphid population, plant mortality, and damage in the susceptible genotype when treated with JA and SA. Together, the results indicate that the set of genes, pathways, and defense compounds is involved in host plant resistance to aphids. These findings shed light on the specific role of each DEG, thus advancing our understanding of the genetic and molecular mechanisms of host plant resistance to aphids.

Analysis of plant expression profiles revealed that aphid attack triggered dynamic defense responses in sorghum plant

Sorghum [Sorghum bicolor (L.) Moench] is one of the most important cereal crops grown worldwide but is often attacked by greenbug (aphid). In response to aphid attack, host plant initiates a large transcriptional reorganization, leading to activation of the host defense genes in aphid-attacked plants. In this study, our objective was to analyze defensive responses of sorghum against aphid and identify aphid resistance genes in sorghum. For the experiments, seedlings developed from an aphid resistant germplasm line (PI 550607) were divided into two groups, then, one group was infested with greenbug ((Schizaphis graminum Rondani) and the other group was used as control (un-infested). In addition, seedlings of sorghum cultivar Tx 7000, a susceptible genotype, prepared under the same conditions, were used as a genetic control. Those plant samples were used to develop transcriptional profiles using the microarray method, from which 26.1% of the 1,761 cDNA sequences spotted on the microarray showed altered expression between two treatments at 4 days after infestation. Sequence annotation and molecular analysis revealed that many differentially expressed genes (DEGs) were related to direct host defense or signal transduction pathways, which regulate host defense. In addition to common responsive genes, unique transcripts were identified in response to greenbug infestation specifically. Later, a similar transcriptional profiling was conducted using the RNA-seq method, resulted in the identification of 2,856 DEGs in the resistant line with a comparison between infested and non-infested at 4 days and 4,354 DEGs in the resistant genotype compared to the susceptible genotype at 4 days. Based on the comparative analysis, the data of RNA-seq provided a support for the results from the microarray study as it was noticed that many of the DEGs are common in both platforms. Analysis of the two differential expression profiles indicate that aphid triggered dynamic defense responses in sorghum plants and sorghum plant defense against aphid is a complex process involving both general defense systems and specific resistance mechanisms. Finally, the results of the study provide new insights into the mechanisms underlying host plant defense against aphids and will help us design better strategies for effectively controlling aphid pest.

Genome-wide identification, evolution and expression pattern analysis of the GATA gene family in Sorghum bicolor

The GATA family of transcription factors is zinc finger DNA binding proteins involved in a variety of biological processes, including plant growth and development and response to biotic/abiotic stresses, and thus play an essential role in plant response to environmental changes. However, the GATA gene family of Sorghum (SbGATA) has not been systematically analyzed and reported yet. Herein, we used a variety of bioinformatics methods and quantitative Real-Time Polymerase Chain Reaction (qRT-PCR) to explore the evolution and function of the 33 SbGATA genes identified. These SbGATA genes, distributed on 10 chromosomes, are classified into four subfamilies (I-IV) containing one pair of tandem duplications and nine pairs of segment duplications, which are more closely related to the monocot Brachypodium distachyon and Oryza sativa GATA genes. The physicochemical properties of the SbGATAs are significantly different among the subfamilies, while the protein structure and conserved protein motifs are highly conserved in the subfamilies. In addition, the transcription of SbGATAs is tissue-specific during Sorghum growth and development, which allows for functional diversity in response to stress and hormones. Collectively, our study lays a theoretical foundation for an in-depth analysis of the functions, mechanisms and evolutionary relationships of SbGATA during plant growth and development.

Transcriptomic and metabolomic profiling reveal the mechanism underlying the inhibition of wound healing by ascorbic acid in fresh-cut potato

Ascorbic acid (AsA) inhibits wound healing in fresh-cut potatoes (FCP); however, the comprehensive regulatory mechanisms of the chemical during wound healing remain unclear. Here, physiobiochemical, transcriptomic, and metabolomic analyses were performed. In total, 685 differentially expressed genes (DEGs) and 1921 differentially accumulated metabolites (DAMs) were identified between control and AsA-treated samples. The level of the majority of DEGs expression and DAMs abundance in AsA-treated samples were similar to data of newly cut samples. The collective data indicated that the AsA treatment inhibited wound healing in FCPs by regulating glutathione metabolism, enhancing starch metabolism, and inhibiting phenylalanine metabolism, sucrose degradation, and fatty acid synthesis. Major genes and metabolites affected by AsA treatment included StGST, StPAL, StPHO1 and StLOX5, and starch, sucrose, and linoleic acid. AsA treatment increased starch content and amylase and lipoxygenase activity and decreased free fatty acid level. Our research provides fundamental insights into wound healing mechanisms in FCP.

Genome-wide identification and characterization of lipoxygenase genes related to the English grain aphid infestation response in wheat

44 wheat LOX genes were identified by silico genome-wide search method. TaLOX5, 7, 10, 24, 29, 33 were specifically expressed post aphid infestation, indicating their participation in wheat-aphid interaction. In plants, LOX genes play important roles in various biological progresses including seed germination, tuber development, plant vegetative growth and most crucially in plant signal transduction, stress response and plant defense against plant diseases and insects. Although LOX genes have been characterized in many species, the importance of the LOX family in wheat has still not been well understood, hampering further improvement of wheat under stress conditions. Here, we identified 44 LOX genes (TaLOXs) in the whole wheat genome and classified into three subfamilies (9-LOXs, Type I 13-LOXs and Type II 13-LOXs) according to phylogenetic relationships. The TaLOXs belonging to the same subgroup shared similar gene structures and motif organizations. Synteny analysis demonstrated that segmental duplication events mainly contributed to the expansion of the LOX gene family in wheat. The results of protein-protein interaction network (PPI) and miRNA-TaLOXs predictions revealed that three TaLOXs (TaLOX20, 22 and 37) interacted mostly with proteins related to methyl jasmonate (MeJA) signaling pathway. The expression patterns of TaLOXs in different tissues (root, stem, leaf, spike and grain) under diverse abiotic stresses (heat, cold, drought, drought and heat combined treatment, and salt) as well as under diverse biotic stresses (powdery mildew pathogen, Fusarium graminearum and stripe rust pathogen) were systematically analyzed using RNA-seq data. We obtained aphid-responsive candidate genes by RNA-seq data of wheat after the English grain aphid infestation. Aphid-responsive candidate genes, including TaLOX5, 7, 10, 24, 29 and 33, were up-regulated in the wheat aphid-resistant genotype (Lunxuan144), while they were little expressed in the susceptible genotype (Jimai22) during late response (48 h and 72 h) to the English grain aphid infestation. Meanwhile, qRT-PCR analysis was used to validate these aphid-responsive candidate genes. The genetic divergence and diversity of all the TaLOXs in bread wheat and its relative species were investigated by available resequencing data. Finally, the 3D structure of the TaLOX proteins was predicted based on the homology modeling method. This study not only systematically investigated the characteristics and evolutionary relationships of TaLOXs, but also provided potential candidate genes in response to the English grain aphid infestation and laid the foundation to further study the regulatory roles in the English grain aphid infestation of LOX family in wheat and beyond.

Genome-wide studies of PAL genes in sorghum and their responses to aphid infestation

Phenylalanine ammonia-lyase (PAL, EC 4.3.1.25) plays a crucial role in plant adaptation to biotic and abiotic stresses. However, the current knowledge about PAL proteins in sorghum is essentially lacking. Thus, in this study we aimed to analyze the PAL family genes in sorghum using a genome-wide approach and to explore the role of PAL genes in host plant resistance to aphids via SA-mediated defense signaling. Here, we report gene structural features of 8 PAL (SbPAL) genes in sorghum (Sorghum bicolor), their phylogeny, protein motifs and promoter analysis. Furthermore, we demonstrated that the SbPAL genes were induced by sugarcane aphid (SCA) infestation and SbPAL exhibited differential gene expression in susceptible and resistant genotypes. PAL activity assays further validated upregulated expression of the SbPAL genes in a resistant genotype. In addition, exogenous application of SA reduced plant damage and suppressed aphid population growth and fecundity in susceptible genotype, suggesting that those SbPAL genes act as positive regulator of the SA-mediated defense signaling pathway to combat aphid pests in sorghum. This study provides insights for further examination of the defense role of PAL in sorghum against other pests and pathogens.

Revealing Differential Expression of Phytohormones in Sorghum in Response to Aphid Attack Using the Metabolomics Approach

Sorghum (Sorghum bicolor) is an important multipurpose crop grown worldwide, but like many other crops, it is often threatened by insect pests. Sugarcane aphid (SCA, Melanaphis sacchari Zehntner), for example, is one of the most severe pests in sorghum, which causes plant damage and yield loss. The main objective of this study was to assess the effect of phytohormones on host plant resistance to aphid attack. Two sorghum genotypes, BTx623 (susceptible) and Tx2783 (resistant), were selected for a comparative analysis of differential expression of a group of phytohormones in response to aphid infestation. The quantification of phytohormones through LC-MS demonstrated higher levels of jasmonic acid (JA), salicylic acid (SA), abscisic acid (ABA), and auxins in the resistant genotype infested with SCA. The PCA plot supports the strong differential responses between resistant and susceptible genotypes, indicating a positive correlation between JA and ABA and a negative correlation between SA and auxins. Similarly, RT-PCR results of the phytohormones-related marker genes showed higher expression in the resistant genotype compared to the susceptible one. Furthermore, to corroborate the role of phytohormones in plant defense, the susceptible genotype was treated with SA, JA, and ABA. The exogenous application of SA and JA + ABA significantly reduced plant mortality, aphid number, and damage in the susceptible genotype, suggesting a strong correlation between phytohormones and plant survival. Our findings indicate that phytohormones play positive roles in plant defense against aphids and provide new insights into the molecular mechanisms operating in plants for self-protection. These findings could also stimulate further research into the mystery about the regulation of phytohormone production during plant interaction with aphids.

Genome-wide identification and expression pattern analysis of lipoxygenase gene family in turnip (Brassica rapa L. subsp. rapa)

Turnip (Brassica rapa L. subsp. rapa) is an important crop with edible and medicinal values, and various stresses, especially salt stress and drought stress, seriously threaten the yield of turnips. LOXs play important roles in regulating plant growth and development, signal transduction, and biotic and abiotic stress responses through secondary metabolites produced by the oxylipin metabolic pathway, and although the turnip genome has been published, however, the role of LOX family genes in various abiotic stress responses has not been systematically studied in turnips. In this study, a total of 15 LOX genes (BrrLOX) were identified in turnip, distributed on six chromosomes. Phylogenetic tree analysis classified these LOX genes into two classes: three 9-LOX proteins and 12 13-LOX type II proteins. Gene duplication analysis showed that tandem and segmental duplication were the main pathways for the expansion of the BrrLOX gene family. The Ka and Ks values of the duplicated genes indicate that the BrrLOX gene underwent strong purifying selection. Further analysis of the cis-acting elements of the promoters suggested that the expression of the BrrLOX gene may be influenced by stress and phytohormones. Transcriptome data analysis showed that 13 BrrLOX genes were expressed at one or more stages of turnip tuber development, suggesting that LOX genes may be involved in the formation of turnip fleshy roots. The qRT-PCR analysis showed that four stresses (salt stress, drought stress, cold stress, and heat stress) and three hormone treatments (methyl jasmonate, salicylic acid, and abscisic acid) affected the expression levels of BrrLOX genes and that different BrrLOX genes responded differently to these stresses. In addition, weighted gene co-expression network analysis (WGCNA) of BrrLOX revealed seven co-expression modules, and the genes in these co-expression modules are collectively involved in plant growth and development and stress response processes. Thus, our results provide valuable information for the functional identification and regulatory mechanisms of BrrLOX in turnip growth and development and stress response.

Genome-wide association mapping of resistance to the sorghum aphid in Sorghum bicolor

Since 2013, the sorghum aphid (SA), Melanaphis sorghi (Theobald), has been a serious pest that hampers all types of sorghum production in the U.S. Known sorghum aphid resistance in sorghum is limited to a few genetic regions on SBI-06. In this study, a subset of the Sorghum Association Panel (SAP) was used along with some additional lines to identify genomic regions that confer sorghum aphid resistance. SAP lines were grown in the field and visually evaluated for SA resistance during the growing seasons of 2019 and 2020 in Tifton, GA. In 2020, the SAP accessions were also evaluated for SA resistance in the field using drone-based high throughput phenotyping (HTP). Flowering time was recorded in the field to confirm that our methods were sufficient for identifying known quantitative trait loci (QTL). This study combined phenotypic data from field-based visual ratings and reflectance data to identify genome-wide associated (GWAS) marker-trait associations (MTA) using genotyping-by-sequencing (GBS) data. Several MTAs were identified for SA-related traits across the genome, with a few common markers that were consistently identified on SBI-08 and SBI-10 for aphid count and plant damage, as well as loci for reflectance-based traits on SBI-02, SBI-03, and SBI-05. Candidate genes encoding leucine-rich repeats (LRR), Avr proteins, lipoxygenases (LOXs), calmodulins (CAM) dependent protein kinase, WRKY transcription factors, flavonoid biosynthesis genes, and 12-oxo-phytodienoic acid reductase were identified near SNPs that had significant associations with different SA traits. In this study, flowering time-related genes were also identified as a positive control for the methods. The total phenotypic variation explained by significant SNPs across SA-scored traits, reflectance data, and flowering time ranged from 6 to 61%, while the heritability value ranged from 4 to 69%. This study identified three new sources of resistant lines to sorghum aphid. These results supported the existing literature, and also revealed several new loci. Markers identified in this study will support marker-assisted breeding for sorghum aphid resistance.

Genome-Wide Analysis of AP2/ERF Superfamily Genes in Contrasting Wheat Genotypes Reveals Heat Stress-Related Candidate Genes

The AP2/ERF superfamily is one of the largest groups of transcription factors (TFs) in plants, which plays important roles in regulating plant growth and development under heat stress. A complete genome-wide identification, characterization, and expression analysis of AP2/ERF superfamily genes focusing on heat stress response were conducted in bread wheat. This study identified 630 putative AP2/ERF superfamily TF genes in wheat, with 517 genes containing well-defined AP2-protein domains. They were classified into five sub-families, according to domain content, conserved motif, and gene structure. The unique genes identified in this study were 112 TaERF genes, 77 TaDREB genes, four TaAP2 genes, and one TaRAV gene. The chromosomal distribution analysis showed the unequal distribution of TaAP2/ERF genes in 21 wheat chromosomes, with 127 pairs of segmental duplications and one pair of tandem duplication, highly concentrated in TaERF and TaDREB sub-families. The qRT-PCR validation of differentially expressed genes (DEGs) in contrasting wheat genotypes under heat stress conditions revealed that significant DEGs in tolerant and susceptible genotypes could unequivocally differentiate tolerant and susceptible wheat genotypes. This study provides useful information on TaAP2/ERF superfamily genes and reveals candidate genes in response to heat stress, which forms a foundation for heat tolerance breeding in wheat.

Genome-Wide Identification, Characterization and Expression Analysis of Lipoxygenase Gene Family in Artemisia annua L

Lipoxygenase (LOX) is a ubiquitous oxygenase found in animals and plants and plays a pivotal role in diverse biological processes, including defense and development. Artemisinin, which can only be obtained from Artemisia annua L., is the most effective therapeutic drug for malaria without serious side effects. This study identified and analyzed LOX gene family members in the A. annua genome at the chromosomal level. Twenty LOX genes with various molecular weights, isoelectric points, and amino acid numbers were identified and named AaLOX, which were located in the cytoplasm or chloroplast. The average protein length of all AaLOX was 850 aa. Phylogenetic tree analysis revealed that the AaLOX was divided into two major groups, 9-LOX and 13-LOX. The exon numbers ranged from 1 to 12, indicating that different AaLOX genes have different functions. The secondary structure was mainly composed of alpha helix and random coil, and the tertiary structure was similar for most AaLOX. Upstream promoter region analysis revealed that a large number of cis-acting elements were closely related to plant growth and development, light response, hormone, and other stress responses. Transcriptome data analysis of different tissues suggested that the gene family was differently expressed in the roots, stems, leaves, and flowers of two A. annua strains HAN1 and LQ9. qRT-PCR confirmed that AaLOX5 and AaLOX17 had the highest expression in flowers and leaves. This study provides a theoretical basis for the further functional analysis of the AaLOX gene family.

Genome-wide characterization of the sorghum JAZ gene family and their responses to phytohormone treatments and aphid infestation

Jasmonate ZIM-domain (JAZ) proteins are the key repressors of the jasmonic acid (JA) signal transduction pathway and play a crucial role in stress-related defense, phytohormone crosstalk and modulation of the growth-defense tradeoff. In this study, the sorghum genome was analyzed through genome-wide comparison and domain scan analysis, which led to the identification of 18 sorghum JAZ (SbJAZ) genes. All SbJAZ proteins possess the conserved TIFY and Jas domains and they formed a phylogenetic tree with five clusters related to the orthologs of other plant species. Similarly, evolutionary analysis indicated the duplication events as a major force of expansion of the SbJAZ genes and there was strong neutral and purifying selection going on. In silico analysis of the promoter region of the SbJAZ genes indicates that SbJAZ5, SbJAZ6, SbJAZ13, SbJAZ16 and SbJAZ17 are rich in stress-related cis-elements. In addition, expression profiling of the SbJAZ genes in response to phytohormones treatment (JA, ET, ABA, GA) and sugarcane aphid (SCA) was performed in two recombinant inbred lines (RILs) of sorghum, resistant (RIL 521) and susceptible (RIL 609) to SCA. Taken together, data generated from phytohormone expression and in silico analysis suggests the putative role of SbJAZ9 in JA-ABA crosstalk and SbJAZ16 in JA-ABA and JA-GA crosstalk to regulate certain physiological processes. Notably, upregulation of SbJAZ1, SbJAZ5, SbJAZ13 and SbJAZ16 in resistant RIL during JA treatment and SCA infestation suggests putative functions in stress-related defense and to balance the plant defense to promote growth. Overall, this report provides valuable insight into the organization and functional characterization of the sorghum JAZ gene family.

The Responses of the Lipoxygenase Gene Family to Salt and Drought Stress in Foxtail Millet (Setaria italica)

Plant lipoxygenases (LOXs), a kind of non-heme iron-containing dioxygenases, participate plant physiological activities (especially in response to biotic and abiotic stresses) through oxidizing various lipids. However, there was few investigations on LOXs in foxtail millet (Setaria italica). In this study, we identified the LOX gene family in foxtail millet, and divided the total 12 members into three sub-families on the basis of their phylogenetic relationships. Under salt and drought stress, LOX genes showed different expression patterns. Among them, only SiLOX7 showed up-regulated expression in Yugu1 (YG1) and Qinhuang2 (QH2), two stress-tolerant varieties, indicating that SiLOX7 may play an important role in responses to abiotic stress. Our research provides a basis for further investigation of the role of LOX genes in the adaptation to abiotic stresses and other possible biological functions in foxtail millet.

VqBGH40a isolated from Chinese wild Vitis quinquangularis degrades trans-piceid and enhances trans-resveratrol

Resveratrol (3,5,4'-trihydroxy-stilbene) is a phytoalexin that can prevent plants from pathogen attacks. Piceid is the glycosylation product of resveratrol and the main storage form of stilbenes in grapevines. Here, we reported the function of a β-glycoside hydrolase gene, VqBGH40a, from the Chinese wild grapevine Vitis quinquangularis accession Danfeng-2 in the regulation of plant resistance to powdery mildew (Uncinula necator). VqBGH40a belonging to β-glycoside hydrolase family 1 encoded 506 amino acids and was located on the cytomembrane. Its optimal induction condition was 28 or 30℃, for 4 h, with 0.1 mM IPTG in a prokaryotic expression system. Enzyme activity detection showed that purified VqBGH40a could hydrolyze trans-piceid to form trans-resveratrol in vitro. VqBGH40a was transiently overexpressed in Danfeng-2 leaves and then artificially inoculated with powdery mildew showed that VqBGH40a protein could hydrolyze trans-piceid in vivo. Additionally, a comparative family analysis between VqBGH40a and 38 VviBGHs was performed. Overall, these results demonstrate that VqBGH40a can hydrolyze trans-piceid, enhance trans-resveratrol content, and participate in the defense mechanism of grapevine against powdery mildew.