Identification of lipoxygenase (LOX) genes from legumes and their responses in wild type and cultivated peanut upon Aspergillus flavus infection

Understanding the role of oxylipins in Cannabis to enhance cannabinoid production

Phytocannabinoids are medically important specialized defense compounds that are sparsely distributed among plants, yet Cannabis sativa can synthesize unprecedented amounts of these compounds within highly specialized surface cell factories known as glandular trichomes. The control mechanisms that allow for this high level of productivity are poorly understood at the molecular level, although increasing evidence supports the role of oxylipin metabolism in phytocannabinoid production. Oxylipins are a large class of lipid-based oxygenated biological signaling molecules. Although some oxylipins are known to participate in plant defense, roles for the majority of the ca. 600 plant oxylipins are largely unknown. In this review, we examine oxylipin gene expression within glandular trichomes and identify key oxylipin genes that determine the fate of common lipid precursors. Mechanisms by which oxylipins may be interacting with phytocannabinoid metabolism, as well as specialized plant metabolism more broadly, are discussed and a model summarizing these contributions proposed.

Comprehensive proteomic profiling of Cucumber mosaic virus infection: identifying key proteins and pathways involved in resistance and susceptibility in melon

BACKGROUND

Melon (Cucumis melo L.) is the model species of the Cucurbitaceae family and an important crop. However, its yield is primarily affected by viruses. Cucumber mosaic virus (CMV) is particularly significant due to its broad host range, capable of infecting over 100 plant families. Resistance to CMV in the melon accession Songwhan Charmi (SC) is controlled by the recessive gene cmv1, which encodes the Vacuolar Protein Sorting 41, involved in vesicle transport to the vacuole. cmv1 restricts the virus to the bundle sheath cells and impedes viral access to the phloem, preventing a systemic infection. This phenotype depends on the viral movement protein (MP). However, little is known about the broader cellular changes that CMV triggers in melon or the specific biological responses that facilitate or restrict the virus entry into the phloem in susceptible and resistant varieties.

RESULT

We profiled the proteomes of CMV-resistant or susceptible melon genotypes inoculated with CMV-LS or FNY strains. Analysis of co-abundance networks revealed the rewiring of central biological pathways during different stages of CMV infection. Upon inoculation, resistant varieties do not trigger any signalling event to the new leaves. Local infection triggers a general depletion in proteins related to translation, photosynthesis and intracellular transport, whereas only in resistant varieties CMV triggers an increase in lipid modification and phloem proteins. During the systemic infection of susceptible melon plants, there is a strong increase in proteins associated with stress responses, such as those involved in the ER-associated degradation (ERAD) and phenylpropanoid pathways, along with a decrease in translation and photosynthesis. Key hub proteins have been identified in these processes.

CONCLUSIONS

This study is the first comprehensive high-throughput proteomic analysis of CMV-infected melon plants, providing a novel and detailed understanding of the proteomic changes associated with CMV infection, highlighting the differential responses between resistant and susceptible genotypes and identifying key proteins that could be potential targets for future research and CMV management strategies.

OsLOX1 positively regulates seed vigor and drought tolerance in rice

The lipoxygenase (LOX) gene family is widely distributed in plants, and its activity is closely associated with seed viability and stress tolerance. In this study, we cloned the rice(Oryza sativa)lipoxygenase gene OsLOX1, a key participant in the 13-lipoxygenase metabolic pathway. Our primary focus was to investigate its role in mediating responses to drought stress and seed germination in rice. Histochemical staining and qPCR analysis indicated that the expression level of OsLOX1 was relatively high in leaves and early germinating seeds. Our findings revealed that mutant lines with CRISPR/Cas9-induced knockout of OsLOX1 exhibited reduced tolerance to drought stress compared with the wild-type. This was accompanied by elevated levels of H2O2 and malondialdehyde, and a decrease in the expression levels of genes associated with antioxidant enzymes. Furthermore, knockout of OsLOX1 reduced the longevity of rice seeds increased H2O2 and MDA levels, and decreased the activities of the antioxidant enzymes superoxide dismutase and catalase, compared with the wild-type. These findings demonstrated that OsLOX1 positively regulated rice seed vigor and drought stress.

Competing endogenous RNA networks and ferroptosis in cancer: novel therapeutic targets

As a newly identified regulated cell death, ferroptosis is a metabolically driven process that relies on iron and is associated with polyunsaturated fatty acyl peroxidation, elevated levels of reactive oxygen species (ROS), and mitochondrial damage. This distinct regulated cell death is dysregulated in various cancers; activating ferroptosis in malignant cells increases cancer immunotherapy and chemoradiotherapy responses across different malignancies. Over the last decade, accumulating research has provided evidence of cross-talk between non-coding RNAs (ncRNAs) and competing endogenous RNA (ceRNA) networks and highlighted their significance in developing and progressing malignancies. Aside from pharmaceutical agents to regulate ferroptosis, recent studies have shed light on the potential of restoring dysregulated ferroptosis-related ceRNA networks in cancer treatment. The present study provides a comprehensive and up-to-date review of the ferroptosis significance, ferroptosis pathways, the role of ferroptosis in cancer immunotherapy and chemoradiotherapy, ceRNA biogenesis, and ferroptosis-regulating ceRNA networks in different cancers. The provided insights can offer the authorship with state-of-the-art findings and future perspectives regarding the ferroptosis and ferroptosis-related ceRNA networks and their implication in the treatment and determining the prognosis of affected patients.

Genome-Wide Characterization of the Phenylalanine Ammonia-Lyase Gene Family and Their Potential Roles in Response to Aspergillus flavus L. Infection in Cultivated Peanut (Arachis hypogaea L.)

Phenylalanine ammonia-lyase (PAL) is an essential enzyme in the phenylpropanoid pathway, in which numerous aromatic intermediate metabolites play significant roles in plant growth, adaptation, and disease resistance. Cultivated peanuts are highly susceptible to Aspergillus flavus L. infection. Although PAL genes have been characterized in various major crops, no systematic studies have been conducted in cultivated peanuts, especially in response to A. flavus infection. In the present study, a systematic genome-wide analysis was conducted to identify PAL genes in the Arachis hypogaea L. genome. Ten AhPAL genes were distributed unevenly on nine A. hypogaea chromosomes. Based on phylogenetic analysis, the AhPAL proteins were classified into three groups. Structural and conserved motif analysis of PAL genes in A. hypogaea revealed that all peanut PAL genes contained one intron and ten motifs in the conserved domains. Furthermore, synteny analysis indicated that the ten AhPAL genes could be categorized into five pairs and that each AhPAL gene had a homologous gene in the wild-type peanut. Cis-element analysis revealed that the promoter region of the AhPAL gene family was rich in stress- and hormone-related elements. Expression analysis indicated that genes from Group I (AhPAL1 and AhPAL2), which had large number of ABRE, WUN, and ARE elements in the promoter, played a strong role in response to A. flavus stress.

A genome-wide study of the lipoxygenase gene families in Medicago truncatula and Medicago sativa reveals that MtLOX24 participates in the methyl jasmonate response

BACKGROUND

Lipoxygenase (LOX) is a multifunctional enzyme that is primarily related to plant organ growth and development, biotic and abiotic stress responses, and production of flavor-associated metabolites. In higher plants, the LOX family encompasses several isozymes with varying expression patterns between tissues and developmental stages. These affect processes including seed germination, seed storage, seedling growth, fruit ripening, and leaf senescence. LOX family genes have multiple functions in response to hormones such as methyl jasmonate (MeJA) and salicylic acid.

RESULTS

In this study, we identified 30 and 95 LOX homologs in Medicago truncatula and Medicago sativa, respectively. These genes were characterized with analyses of their basic physical and chemical properties, structures, chromosomal distributions, and phylogenetic relationships to understand structural variations and their physical locations. Phylogenetic analysis was conducted for members of the three LOX subfamilies (9-LOX, type I 13-LOX, and type II 13-LOX) in Arabidopsis thaliana, Glycine max, M. truncatula, and M. sativa. Analysis of predicted promoter elements revealed several relevant cis-acting elements in MtLOX and MsLOX genes, including abscisic acid (ABA) response elements (ABREs), MeJA response elements (CGTCA-motifs), and antioxidant response elements (AREs). Cis-element data combined with transcriptomic data demonstrated that LOX gene family members in these species were most likely related to abiotic stress responses, hormone responses, and plant development. Gene expression patterns were confirmed via quantitative reverse transcription PCR. Several MtLOX genes (namely MtLOX15, MtLOX16, MtLOX20, and MtLOX24) belonging to the type I 13-LOX subfamily and other LOX genes (MtLOX7, MtLOX11, MsLOX23, MsLOX87, MsLOX90, and MsLOX94) showed significantly different expression levels in the flower tissue, suggesting roles in reproductive growth. Type I 13-LOXs (MtLOX16, MtLOX20, MtLOX21, MtLOX24, MsLOX57, MsLOX84, MsLOX85, and MsLOX94) and type II 13-LOXs (MtLOX5, MtLOX6, MtLOX9, MtLOX10, MsLOX18, MsLOX23, and MsLOX30) were MeJA-inducible and were predicted to function in the jasmonic acid signaling pathway. Furthermore, exogenous MtLOX24 expression in Arabidopsis verified that MtLOX24 was involved in MeJA responses, which may be related to insect-induced abiotic stress.

CONCLUSIONS

We identified six and four LOX genes specifically expressed in the flowers of M. truncatula and M. sativa, respectively. Eight and seven LOX genes were induced by MeJA in M. truncatula and M. sativa, and the LOX genes identified were mainly distributed in the type I and type II 13-LOX subfamilies. MtLOX24 was up-regulated at 8 h after MeJA induction, and exogenous expression in Arabidopsis demonstrated that MtLOX24 promoted resistance to MeJA-induced stress. This study provides valuable new information regarding the evolutionary history and functions of LOX genes in the genus Medicago.

Identification and characterization of lipoxygenase (LOX) genes involved in abiotic stresses in yellow horn

Lipoxygenase (LOX) gene plays an essential role in plant growth, development, and stress response. 15 LOX genes were identified, which were unevenly distributed on chromosomes and divided into three subclasses in this study. In promoter region analysis, many cis-elements were identified in growth and development, abiotic stress response, hormonal response, and light response. qRT-PCR showed that the LOX gene showed tissue specificity in seven tissues, especially XsLOX1, 3, and 7 were relatively highly expressed in roots, stems, and axillary buds. The different expression patterns of LOX genes in response to abiotic stress and hormone treatment indicate that different XsLOX genes have different reactions to these stresses and play diversified roles. This study improves our understanding of the mechanism of LOX regulation in plant growth, development, and stress and lays a foundation for further analysis of biological functions.

Recombinant Soybean Lipoxygenase 2 (GmLOX2) Acts Primarily as a ω6(S)-Lipoxygenase

The lipoxygenase (LOX) cascade is a source of bioactive oxylipins that play a regulatory role in plants, animals, and fungi. Soybean (Glycine max (L.) Merr.) LOXs are the classical models for LOX research. Progress in genomics has uncovered a large diversity of GmLOX isoenzymes. Most of them await biochemical investigations. The catalytic properties of recombinant soybean LOX2 (GmLOX2) are described in the present work. The GmLOX2 gene has been cloned before, but only for nucleotide sequencing, while the recombinant protein was not prepared and studied. In the present work, the recombinant GmLOX2 behavior towards linoleic, α-linolenic, eicosatetraenoic (20:4), eicosapentaenoic (20:5), and hexadecatrienoic (16:3) acids was examined. Linoleic acid was a preferred substrate. Oxidation of linoleic acid afforded 94% optically pure (13S)-hydroperoxide and 6% racemic 9-hydroperoxide. GmLOX2 was less active on other substrates but possessed an even higher degree of regio- and stereospecificity. For example, it converted α-linolenic acid into (13S)-hydroperoxide at about 98% yield. GmLOX2 showed similar specificity towards other substrates, producing (15S)-hydroperoxides (with 20:4 and 20:5) or (11S)-hydroperoxide (with 16:3). Thus, the obtained data demonstrate that soybean GmLOX2 is a specific (13S)-LOX. Overall, the catalytic properties of GmLOX2 are quite similar to those of GmLOX1, but pH is optimum.

Exogenous Methyl Jasmonate Mediated MiRNA-mRNA Network Improves Heat Tolerance of Perennial Ryegrass

Heat stress can hinder the growth of perennial ryegrass (Lolium perenne L.). Methyl jasmonate (MeJA) applied exogenously can increase heat stress tolerance in plants; however, the regulatory mechanisms involved in heat tolerance mediated by MeJA are poorly understood in perennial ryegrass. Here, the microRNA (miRNA) expression profiles of perennial ryegrass were assessed to elucidate the regulatory pathways associated with heat tolerance induced by MeJA. Plants were subjected to four treatments, namely, control (CK), MeJA pre-treatment (T), heat stress treatment (H), and MeJA pre-treatment + heat stress (TH). According to the results, 102 miRNAs were up-regulated in all treatments, with 20, 27, and 33 miRNAs being up-regulated in the T, H, and TH treatment groups, respectively. The co-expression network analysis between the deferentially expressed miRNAs and their corresponding target genes showed that 20 miRNAs modulated 51 potential target genes. Notably, the miRNAs that targeted genes related to with regards to heat tolerance were driven by MeJA, and they were involved in four pathways: novel-m0258-5p mediated signal transduction, novel-m0350-5p mediated protein homeostasis, miR397-z, miR5658-z, and novel-m0008-5p involved in cell wall component, and miR1144-z and miR5185-z dominated chlorophyll degradation. Overall, the findings of this research paved the way for more research into the heat tolerance mechanism in perennial ryegrass and provided a theoretical foundation for developing cultivars with enhanced heat tolerance.

Genome-wide survey of catalase genes in Brassica rapa, Brassica oleracea, and Brassica napus: identification, characterization, molecular evolution, and expression profiling of BnCATs in response to salt and cadmium stress

Catalase (CAT, EC 1.11.1.6), one of the most important antioxidant enzymes, can control excess levels of H₂O₂ produced under oxidative stress in plants. In this study, 16, 8, and 7 CAT genes in the genome of Brassica napus, B. rapa, and B. oleracea were identified, respectively. Phylogenetic studies showed that CATs could be divided into two main groups, each containing specific monocotyledon and dicotyledon subgroups. Motifs, gene structure, and intron phase of CATs in B. napus, Brassica rapa, and Brassica oleracea are highly conserved. Analysis of codon usage bias showed the mutation pressure and natural selection of the codon usage of CATs. Segmental duplication and polyploid were major factors in the expansion of this gene family in B. napus, and genes have experienced  negative selection during evolution. Existence of hormones and stress-responsive cis-elements and identifying miRNA molecules affecting CATs showed that these genes are complexly regulated at the transcriptional and posttranscriptional levels. Based on RNA-seq data, CATs are divided into two groups; the first group has moderate and specific expression in flowers, leaves, stems, and roots, while the second group shows expression in most tissues. qRT-PCR analysis showed that the expression of these genes is dynamic and has a specific expression consistent with other CAT genes in response to salinity and cadmium (Cd) stresses. These results provide information for further investigation of the function of CAT genes in response to stresses and the development of tolerant  plants.

Genome-Wide Analysis of Lipoxygenase (LOX) Genes in Angiosperms

Lipoxygenases (LOXs) are enzymes that catalyze the addition of an oxygen molecule to unsaturated fatty acids, thus forming hydroperoxides. In plants, these enzymes are encoded by a multigene family found in several organs with varying activity patterns, by which they are classified as LOX9 or LOX13. They are involved in several physiological functions, such as growth, fruit development, and plant defense. Despite several studies on genes of the LOX family in plants, most studies are restricted to a single species or a few closely related species. This study aimed to analyze the diversity, evolution, and expression of LOX genes in angiosperm species. We identified 247 LOX genes among 23 species of angiosperms and basal plants. Phylogenetic analyses identified clades supporting LOX13 and two main clades for LOX9: LOX9_A and LOX9_B. Eudicot species such as Tarenaya hassleriana, Capsella rubella, and Arabidopsis thaliana did not present LOX9_B genes; however, LOX9_B was present in all monocots used in this study. We identified that there were potential new subcellular localization patterns and conserved residues of oxidation for LOX9 and LOX13 yet unexplored. In summary, our study provides a basis for the further functional and evolutionary study of lipoxygenases in angiosperms.

Genome-wide identification and analysis of LOX genes in soybean cultivar “Zhonghuang 13”

Lipoxygenases (LOXs; EC1.13.11.12) are a family of iron- or manganese-containing dioxygenases that catalyze the oxygenation of polyunsaturated fatty acids (PUFAs) and play important roles in plant growth, development, and stress response. In this study, a total of 36 LOX gene family members were identified and annotated in Zhonghuang 13, a soybean cultivar bred by Chinese scientists in 2001. Sanger sequencing of the GmLOX1-coding sequence and colorimetric assays for the GmLOX1 protein showed that Zhonghuang 13 possessed the GmLOX1 gene. These LOX genes are divided into three subfamilies: 9-LOX, type Ⅰ 13-LOX and type II 13-LOX. In the 13-LOX group, the number of GmLOX members was the highest. These GmLOX genes are unevenly distributed on chromosomes 3, 7, 8, 10, 11, 12, 13, 15, 16, 19, and 20. Most of the 13-LOX genes exist in the form of gene clusters, indicating that these genes may originate from tandem duplications. The analysis of duplicated gene pairs showed that GmLOX genes underwent purifying selective pressure during evolution. The gene structures and conserved functional domains of these genes are quite similar. Compared to the orthologous gene pairs of LOX genes between wild soybean (Glycine soja W05) and Zhonghuang 13, the sequences of most gene pairs are relatively conserved. Many cis-elements are present in the promoter region and are involved in stress response, growth and development, hormone response and light response. The tissue-specific gene expression of GmLOX genes was evaluated. Represented by GmLOX1, GmLOX2, and GmLOX3, which were expressed at extremely high levels in seeds, they showed the characteristics of specific expression. This study provides detailed information on soybean lipoxygenase gene family members in Zhonghuang 13, which lays a foundation for further research.

De novo genes in Arachis hypogaea cv. Tifrunner: systematic identification, molecular evolution, and potential contributions to cultivated peanut

De novo genes are derived from non-coding sequences, and they can play essential roles in organisms. Cultivated peanut (Arachis hypogaea) is a major oil and protein crop derived from a cross between Arachis duranensis and Arachis ipaensis. However, few de novo genes have been documented in Arachis. Here, we identified 381 de novo genes in A. hypogaea cv. Tifrunner based on comparison with five closely related Arachis species. There are distinct differences in gene expression patterns and gene structures between conserved and de novo genes. The identified de novo genes originated from ancestral sequence regions associated with metabolic and biosynthetic processes, and they were subsequently integrated into existing regulatory networks. De novo paralogs and homoeologs were identified in A. hypogaea cv. Tifrunner. De novo paralogs and homoeologs with conserved expression have mismatching cis-acting elements under normal growth conditions. De novo genes potentially have pluripotent functions in responses to biotic stresses as well as in growth and development based on quantitative trait locus data. This work provides a foundation for future research examining gene birth processes and gene function in Arachis and related taxa.

Specific Roles of Lipoxygenases in Development and Responses to Stress in Plants

Lipoxygenases (LOXs), naturally occurring enzymes, are widely distributed in plants and animals. LOXs can be non-sulfur iron, non-heme iron, or manganese-containing dioxygenase redox enzymes. LOXs catalyze the oxidation of polyunsaturated fatty acids into fatty acid hydroperoxides. Linolenic acid, a precursor in the jasmonic acid (JA) biosynthesis, is converted to 12-oxo-phytodienoic acid through oxygenation with LOX, allene oxide synthase, and allene oxide cyclase. Moreover, JA participates in seed germination, fruit ripening, senescence, and many other physio-biochemical processes. LOXs also play crucial roles in defense responses against biotic stress, i.e., insects, pests, pathogenic attacks, and abiotic stress, such as wounding, UV-rays, extreme temperature, oxidative stress, and drought.

Identification of the LOX Gene Family in Peanut and Functional Characterization of AhLOX29 in Drought Tolerance

Lipoxygenases (LOXs) are a gene family of nonheme iron-containing dioxygenases that play important roles in plant development and defense responses. To date, a comprehensive analysis of LOX genes and their biological functions in response to abiotic stresses in peanut has not been performed. In this study, a total of 72 putative LOX genes were identified in cultivated (Arachis hypogaea) and wild-type peanut (Arachis duranensis and Arachis ipaensis) and classified into three subfamilies: 9-LOX, type I 13-LOX and type II 13-LOX. The gene structures and protein motifs of these peanut LOX genes were highly conserved among most LOXs. We found that the chromosomal distribution of peanut LOXs was not random and that gene duplication played a crucial role in the expansion of the LOX gene family. Cis-acting elements related to development, hormones, and biotic and abiotic stresses were identified in the promoters of peanut LOX genes. The expression patterns of peanut LOX genes were tissue-specific and stress-inducible. Quantitative real-time PCR results further confirmed that peanut LOX gene expression could be induced by drought, salt, methyl jasmonate and abscisic acid treatments, and these genes exhibited diverse expression patterns. Furthermore, overexpression of AhLOX29 in Arabidopsis enhanced the resistance to drought stress. Compared with wide-type, AhLOX29-overexpressing plants showed significantly decreased malondialdehyde contents, as well as increased chlorophyll degradation, proline accumulation and superoxide dismutase activity, suggesting that the transgenic plants exhibit strengthened capacity to scavenge reactive oxygen species and prevent membrane damage. This systematic study provides valuable information about the functional characteristics of AhLOXs in the regulation of abiotic stress responses of peanut.

Unraveling Microbial Volatile Elicitors Using a Transparent Methodology for Induction of Systemic Resistance and Regulation of Antioxidant Genes at Expression Levels in Chili against Bacterial Wilt Disease

Microbial volatiles benefit the agricultural ecological system by promoting plant growth and systemic resistance against diseases without harming the environment. To explore the plant growth-promoting efficiency of VOCs produced by Pseudomonas fluorescens PDS1 and Bacillus subtilis KA9 in terms of chili plant growth and its biocontrol efficiency against Ralstonia solanacearum, experiments were conducted both in vitro and in vivo. A closure assembly was designed using a half-inverted plastic bottle to demonstrate plant–microbial interactions via volatile compounds. The most common volatile organic compounds were identified and reported; they promoted plant development and induced systemic resistance (ISR) against wilt pathogen R. solanacearum. The PDS1 and KA9 VOCs significantly increased defensive enzyme activity and overexpressed the antioxidant genes PAL, POD, SOD, WRKYa, PAL1, DEF-1, CAT-2, WRKY40, HSFC1, LOX2, and NPR1 related to plant defense. The overall gene expression was greater in root tissue as compared to leaf tissue in chili plant. Our findings shed light on the relationship among rhizobacteria, pathogen, and host plants, resulting in plant growth promotion, disease suppression, systemic resistance-inducing potential, and antioxidant response with related gene expression in the leaf and root tissue of chili.

Genome-wide identification and expression pattern analysis of lipoxygenase gene family in banana

The LOX genes have been identified and characterized in many plant species, but studies on the banana LOX genes are very limited. In this study, we respectively identified 18 MaLOX, 11 MbLOX, and 12 MiLOX genes from the Musa acuminata, M. balbisiana and M. itinerans genome data, investigated their gene structures and characterized the physicochemical properties of their encoded proteins. Banana LOXs showed a preference for using and ending with G/C and their encoded proteins can be classified into 9-LOX, Type I 13-LOX and Type II 13-LOX subfamilies. The expansion of the MaLOXs might result from the combined actions of genome-wide, tandem, and segmental duplications. However, tandem and segmental duplications contribute to the expansion of MbLOXs. Transcriptome data based gene expression analysis showed that MaLOX1, 4, and 7 were highly expressed in fruit and their expression levels were significantly regulated by ethylene. And 11, 12 and 7 MaLOXs were found to be low temperature-, high temperature-, and Fusarium oxysporum f. sp. Cubense tropical race 4 (FocTR4)-responsive, respectively. MaLOX8, 9 and 13 are responsive to all the three stresses, MaLOX4 and MaLOX12 are high temperature- and FocTR4-responsive; MaLOX6 and MaLOX17 are significantly induced by low temperature and FocTR4; and the expression of MaLOX7 and MaLOX16 are only affected by high temperature. Quantitative real-time PCR (qRT-PCR) analysis revealed that the expression levels of several MaLOXs are regulated by MeJA and FocTR4, indicating that they can increase the resistance of banana by regulating the JA pathway. Additionally, the weighted gene co-expression network analysis (WGCNA) of MaLOXs revealed 3 models respectively for 5 (MaLOX7-11), 3 (MaLOX6, 13, and 17), and 1 (MaLOX12) MaLOX genes. Our findings can provide valuable information for the characterization, evolution, diversity and functionality of MaLOX, MbLOX and MiLOX genes and are helpful for understanding the roles of LOXs in banana growth and development and adaptations to different stresses.

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

This report shows detailed characterization of LOX gene family in sorghum and provides new insight of sorghum LOX genes in genetic structure and their roles in plant response to infestation by sugarcane aphids. Lipoxygenases (LOXs) are monomeric, nonheme iron-containing dioxygenases that initiate the fatty acid oxidation pathway creating oxylipins and plant hormone jasmonate both have a key role in plant development and defense. To date, a comprehensive and systematic analysis of sorghum LOXs is still deficient. Thus, we performed a genome-wide analysis of the sorghum LOXs genome and identified nine LOXs genes. Detailed examination of protein sequences and phylogenetic analysis categorized the sorghum LOXs into two subclasses, 9-LOXs (SbLOX1, SbLOX3, SbLOX4, SbLOXm, and SbLOXo), 13-LOXs (SbLOX9, SbLOX5, and SbLOX2), and the unclassified SbLOX8. This classification was further supported by sequence similarity/identity matrix and subcellular localization analysis. The lipoxygenase domains, motifs, and vital amino acids were highly conserved in all sorghum LOX genes. In silico analysis of the promoter region of SbLOXs identified different hormones responsive cis-elements. Furthermore, to explore the roles of sorghum LOXs during sugarcane aphid feeding and exogenous MeJA application, expression analysis was conducted for all the eight LOXs in resistant (Tx2783) and susceptible (Tx7000) sorghum lines, respectively. As detailed in this report, the data generated from both genome-wide identification and expression analysis of lipoxygenase genes suggest the putative functions of two 13-LOXs (SbLOX9 and SbLOX5) and three 9-LOXs (SbLOX1, SbLOX3, and SbLOXo) in biosynthesis of jasmonic acid, green leaf volatiles and death acids, and all of them are involved in defense-related functions in plants. Furthermore, this report represents the first genome-wide analysis of the LOX gene family in sorghum, which will facilitate future studies to characterize the roles of each individual LOXs gene in aphid resistance and defense responses to other stresses.

Magnaporthe oryzae systemic defense trigger 1 (MoSDT1)-mediated metabolites regulate defense response in Rice

BACKGROUND

Some of the pathogenic effector proteins play an active role in stimulating the plant defense system to strengthen plant resistance.

RESULTS

In this study, ultra-high performance liquid chromatography-quadrupole time-of-flight mass spectrometry (UHPLC/Q-TOF-MS) was implemented to identify altered metabolites in transgenic rice containing over-expressed M. oryzae Systemic Defense Trigger 1 (MoSDT1) that was infected at three-time points. The characterized dominating metabolites were organic acids and their derivatives, organic oxygen compounds, lipids, and lipid-like molecules. Among the identified metabolites, shikimate, galactinol, trehalose, D-mannose, linolenic acid, dopamine, tyramine, and L-glutamine are precursors for the synthesis of many secondary defense metabolites Carbohydrate, as well as amino acid metabolic, pathways were revealed to be involved in plant defense responses and resistance strengthening.

CONCLUSION

The increasing salicylic acid (SA) and jasmonic acid (JA) content enhanced interactions between JA synthesis/signaling gene, SA synthesis/receptor gene, raffinose/fructose/sucrose synthase gene, and cell wall-related genes all contribute to defense response in rice. The symptoms of rice after M. oryzae infection were significantly alleviated when treated with six identified metabolites, i.e., galactol, tyramine, L-glutamine, L-tryptophan, α-terpinene, and dopamine for 72 h exogenously. Therefore, these metabolites could be utilized as an optimal metabolic marker for M. oryzae defense.

Identification, Phylogeny, and Comparative Expression of the Lipoxygenase Gene Family of the Aquatic Duckweed, Spirodela polyrhiza, during Growth and in Response to Methyl Jasmonate and Salt

Lipoxygenases (LOXs) (EC 1.13.11.12) catalyze the oxygenation of fatty acids and produce oxylipins, including the plant hormone jasmonic acid (JA) and its methyl ester, methyl jasmonate (MeJA). Little information is available about the LOX gene family in aquatic plants. We identified a novel LOX gene family comprising nine LOX genes in the aquatic plant Spirodela polyrhiza (greater duckweed). The reduced anatomy of S. polyrhiza did not lead to a reduction in LOX family genes. The 13-LOX subfamily, with seven genes, predominates, while the 9-LOX subfamily is reduced to two genes, an opposite trend from known LOX families of other plant species. As the 13-LOX subfamily is associated with the synthesis of JA/MeJA, its predominance in the Spirodela genome raises the possibility of a higher requirement for the hormone in the aquatic plant. JA-/MeJA-based feedback regulation during culture aging as well as the induction of LOX gene family members within 6 h of salt exposure are demonstrated.

Comparative Genomic and Secretomic Analysis Provide Insights Into Unique Agar Degradation Function of Marine Bacterium Vibrio fluvialis A8 Through Horizontal Gene Transfer

Agarose-oligosaccharide production from agar degradation by agarase exhibits lots of advantages and good application prospects. In this study, a novel agar-degrading bacterium Vibrio sp. A8 was isolated from a red algae in the South China Sea. The whole genome sequencing with comparative genomic and secretomic analysis were used to better understand its genetic components about agar degradation. This strain exhibited good agarase production in artificial seawater after culture optimization. The complete genome (4.88 Mb) of this strain comprised two circular chromosomes (3.19 and 1.69 Mb) containing 4,572 protein-coding genes, 108 tRNA genes and 31 rRNA genes. This strain was identified as Vibrio fluvialis A8 by comparative genomic analysis based on genome phylogenetic tree and average nucleotide identity (ANI) similarity. Different from other 20 similar strains including three strains of the same species, V. fluvialis A8 possessed unique agar degradation ability with four β-agarases (GH50) and one α-1,3-L-NA2 hydrolase (GH117) due to the horizontal gene transfer. Secretomic analysis showed that only β-agarase (gene 3152) was abundantly expressed in the secretome of V. fluvialis A8. This agarase had a good substrate specificity and wide work conditions in complex environments, suggesting its potential application for agarose-oligosaccharide production.

Comprehensive identification and analysis of DELLA genes throughout the plant kingdom

BACKGROUND

DELLAs play key roles in plant gibberellin signaling pathways and are generally important in plant development and growth. However, DELLAs in many plant taxa have not yet been systematically analyzed.

RESULTS

In our study, we searched for DELLA genes across 58 green plant genomes and found 181 DELLAs. Structure analysis showed some DELLA domains do not contain "D-E-L-L-A" sequences and instead contain similar domains, including DGLLA and DSLLH domains. "VHYNP" motifs in plant DELLAs comprise 23 types of sequences, while some DELLAs did not contain GRAS domains. In grape, we found that the DELLA protein GSVIVT01015465001 contains an F-box domain, while apple DELLA proteins MDP0000220512 and MDP0000403162 contain a WW domain and a BCIP domain, respectively. These DELLAs can be divided into 22 homologous groups and 17 orthologous groups, and 35 paralogous genes were identified. In total, 35 positively selected genes (PSGs) and 121 negatively selected genes (NSGs) were found among DELLAs based on selective pressure analysis, with an average K_s of NSGs that was significantly higher than that of PSGs (P < 0.05). Among the paralogous groups, CBI and Fop were significantly positively correlated with GC, GC1, GC2, GC12, and GC3, while CAI was significantly positively correlated with GC, GC1, GC12, and GC. The paralogous groups with ω values exceeding 1 had significantly higher K_a values. We also found some paralogous groups with ω values exceeding 1 that differed in their motifs.

CONCLUSIONS

This study provides helpful insights into the evolution of DELLA genes and offers exciting opportunities for the investigation of DELLA functions in different plants.

High-density SNP map facilitates fine mapping of QTLs and candidate genes discovery for Aspergillus flavus resistance in peanut (Arachis hypogaea)

Two novel resistant QTLs mapped and candidate genes identified for Aspergillus flavus resistance in cultivated peanut using SLAF-seq. Aflatoxin contamination in peanuts caused by Aspergillus flavus is a serious food safety issue for human health around the world. Host plant resistance to fungal infection and reduction in aflatoxin are crucial for mitigating this problem. Identification of the resistance-linked markers can be used in marker-assisted breeding for varietal development. Here we report construction of two high-density genetic linkage maps with 1975 SNP loci and 5022 SNP loci, respectively. Two consistent quantitative trait loci (QTL) were identified as qRAF-3-1 and qRAF-14-1, which located on chromosomes A03 and B04, respectively. QTL qRAF-3-1 was mapped within 1.67 cM and had more than 19% phenotypic variance explained (PVE), while qRAF-14-1 was located within 1.34 cM with 5.15% PVE. While comparing with the reference genome, the mapped QTLs, qRAF-3-1 and qRAF-14-1, were located within a physical distance of 1.44 Megabase pair (Mbp) and 2.22 Mbp, harboring 67 and 137 genes, respectively. Among the identified candidate genes, six genes with the same function were found within both QTLs regions. In addition, putative disease resistance RPP13-like protein 1 (RPP13), lipoxygenase (Lox), WRKY transcription factor (WRKY) and cytochrome P450 71B34 genes were also identified. Using microarray analysis, genes responded to A. flavus infection included coding for RPP13, pentatricopeptide repeat-containing-like protein, and Lox which may be possible candidate genes for resistance to A. flavus. The QTLs and candidate genes will further facilitate marker development and validation of genes for deployment in the molecular breeding programs against A. flavus in peanuts.

Transcript Abundance Patterns of 9- and 13-Lipoxygenase Subfamily Gene Members in Response to Abiotic Stresses (Heat, Cold, Drought or Salt) in Tomato (Solanum lycopersicum L.) Highlights Member-Specific Dynamics Relevant to Each Stress

Lipoxygenases (LOXs; EC 1.13.11.12) catalyze the oxygenation of fatty acids to produce oxylipins including the jasmonate family of plant hormones. The involvement of jasmonates in plant growth and development and during abiotic stress has been documented, however, the response and regulation of each member of the LOX gene family under various abiotic stresses is yet to be fully deciphered. Previously, we identified fourteen members of the tomato LOX gene family, which were divisible into nine genes representing the 9-LOX family members and five others representing the 13-LOX family members based on the carbon oxidation position specificity of polyunsaturated fatty acids. Here, we have determined the transcript abundance patterns of all the 14 LOX genes in response to four independent abiotic stresses, namely, heat, cold, drought and salt. Our results show that each of these stresses leads to a time-dependent, variable or indifferent response of specific and different set(s) of LOX gene members of both subfamilies, differentiating functional relevance of the 14 LOX genes analyzed. Out of the 14 gene members, three LOX genes were expressed constitutively or were non-responsive to either heat (SlLOX9), cold (SlLOX9) or salt (SlLOX4) stress. An in-silico LOX gene promoter search for stress-responsive elements revealed that only some but not all of the LOX genes indeed are decorated with specific and known stress responsive cis-acting elements. Thus, these data implicate some other, yet to be discovered, cis-acting elements present in the LOX gene family members, which seemingly regulate tomato responses to defined abiotic stresses presented here.

Genome Editing, Gene Drives, and Synthetic Biology: Will They Contribute to Disease-Resistant Crops, and Who Will Benefit?

Genetically engineered crops have been grown for more than 20 years, resulting in widespread albeit variable benefits for farmers and consumers. We review current, likely, and potential genetic engineering (GE) applications for the development of disease-resistant crop cultivars. Gene editing, gene drives, and synthetic biology offer novel opportunities to control viral, bacterial, and fungal pathogens, parasitic weeds, and insect vectors of plant pathogens. We conclude that there will be no shortage of GE applications totackle disease resistance and other farmer and consumer priorities for agricultural crops. Beyond reviewing scientific prospects for genetically engineered crops, we address the social institutional forces that are commonly overlooked by biological scientists. Intellectual property regimes, technology regulatory frameworks, the balance of funding between public- and private-sector research, and advocacy by concerned civil society groups interact to define who uses which GE technologies, on which crops, and for the benefit of whom. Ensuring equitable access to the benefits of genetically engineered crops requires affirmative policies, targeted investments, and excellent science.

A halotolerant growth promoting rhizobacteria triggers induced systemic resistance in plants and defends against fungal infection

A halotolerant rhizobacteria, Klebsiella species (referred to MBE02), was identified that had a growth stimulation effect on peanut. To gain mechanistic insights into how molecular components were reprogrammed during the interaction of MBE02 and peanut roots, we performed deep RNA-sequencing. In total, 1260 genes were differentially expressed: 979 genes were up-regulated, whereas 281 were down-regulated by MBE02 treatment as compared to uninoculated controls. A large component of the differentially regulated genes were related to phytohormone signalling. This included activation of a significant proportion of genes involved in jasmonic acid, ethylene and pathogen-defense signalling, which indicated a role of MBE02 in modulating plant immunity. In vivo and in vitro pathogenesis assays demonstrated that MBE02 treatment indeed provide fitness benefits to peanut against Aspergillus infection under controlled as well as field environment. Further, MBE02 directly reduced the growth of a wide range of fungal pathogens including Aspergillus. We also identified possible molecular components involved in rhizobacteria-mediated plant protection. Our results show the potential of MBE02 as a biocontrol agent in preventing infection against several fungal phytopathogens.

Gracilaria dura extract confers drought tolerance in wheat by modulating abscisic acid homeostasis

Water stress severely reduces the production of wheat. Application of seaweed extracts have started to show promise in protecting plants from environmental stresses as they contain several biostimulants. However, the modes of action of these biostimulants are not clear. Here, we investigated the role of Gracilaria dura (GD), a red alga, in conferring stress tolerance to wheat during drought under glasshouse and agro-ecological conditions by integrating molecular studies with physiological and field investigations. GD-sap application conferred drought tolerance (as the biomass increased by up to 57% and crop yield by 70%), via facilitating physiological changes associated to maintaining higher water content. GD-sap application significantly increased ABA accumulation (2.34 and 1.46 fold at 4 and 6 days of drought, respectively) due to enhanced expression of biosynthesis genes. This followed an activation of ABA response genes and physiological processes including reduced stomatal opening, thus reducing water loss. Moreover, GD-sap application enhanced the expression of stress-protective genes specifically under water stress. Treatment with fluridone, an ABA inhibitor, further support the role of ABA in GD-sap mediated drought tolerance in wheat. The findings of this study provide insights into the functional role of GD-sap in improving drought tolerance and show the potential to commercialize GD-sap as a potent biostimulant for sustainable agriculture in regions prone to drought.

The soybean Rhg1 amino acid transporter gene alters glutamate homeostasis and jasmonic acid-induced resistance to soybean cyst nematode

Rhg1 (resistance to Heterodera glycines 1) is an important locus that contributes to resistance against soybean cyst nematode (SCN; Heterodera glycines Ichinohe), which is the most economically damaging disease of soybean worldwide. Simultaneous overexpression of three genes encoding a predicted amino acid transporter, an α-soluble N-ethylmaleimide-sensitive factor attachment protein (α-SNAP) and a predicted wound-induced protein resulted in resistance to SCN provided by this locus. However, the roles of two of these genes (excluding α-SNAP) remain unknown. Here, we report the functional characterization of Glyma.18G022400, a gene at the Rhg1 locus that encodes the predicted amino acid transporter Rhg1-GmAAT. Although the direct role of Rhg1-GmAAT in glutamate transport was not demonstrated, multiple lines of evidence showed that Rhg1-GmAAT impacts glutamic acid tolerance and glutamate transportation in soybean. Transcriptomic and metabolite profiling indicated that overexpression of Rhg1-GmAAT activated the jasmonic acid (JA) pathway. Treatment with a JA biosynthesis inhibitor reduced the resistance provided by the Rhg1-containing PI88788 to SCN, which suggested that the JA pathway might play a role in Rhg1-mediated resistance to SCN. Our results could be helpful for the clarification of the mechanism of resistance to SCN provided by Rhg1 in soybean.

Complete Genome Sequence Analysis and Characterization of Selected Iron Regulation Genes of Pasteurella Multocida Serotype A Strain PMTB2.1

Although more than 100 genome sequences of Pasteurella multocida are available, comprehensive and complete genome sequence analysis is limited. This study describes the analysis of complete genome sequence and pathogenomics of P. multocida strain PMTB2.1. The genome of PMTB2.1 has 2176 genes with more than 40 coding sequences associated with iron regulation and 140 virulence genes including the complete tad locus. The tad locus includes several previously uncharacterized genes such as flp2, rcpC and tadV genes. A transposable phage resembling to Mu phages was identified in P. multocida that has not been identified in any other serotype yet. The multi-locus sequence typing analysis assigned the PMTB2.1 genome sequence as type ST101, while the comparative genome analysis showed that PMTB2.1 is closely related to other P. multocida strains with the genomic distance of less than 0.13. The expression profiling of iron regulating-genes of PMTB2.1 was characterized under iron-limited environment. Results showed significant changes in the expression profiles of iron-regulating genes (p < 0.05) whereas the highest expression of fecE gene (281 fold) at 30 min suggests utilization of the outer-membrane proteins system in iron acquisition at an early stage of growth. This study showed the phylogenomic relatedness of P. multocida and improved annotation of important genes and functional characterization of iron-regulating genes of importance to the bacterial growth.

Genome-wide identification of LRR-containing sequences and the response of these sequences to nematode infection in Arachis duranensis

BACKGROUND

Leucine-rich repeat (LRR)-containing genes are involved in responses to various diseases. Recently, RNA-seq data from A. duranensis after nematode (Meloidogyne arenaria) infection were released. However, the number of LRR-containing genes present in A. duranensis and the response of LRR-containing genes to nematode infection are poorly understood.

RESULTS

In this study, we found 509 amino acid sequences containing nine types of LRR domains in A. duranensis. The inferred phylogenetic relationships revealed that the nine types of LRR domains had two originations. The inferred selective pressure was mainly consistent with LRR domains undergoing purifying selection. Twenty-one LRR-containing genes were associated with possible resistance to nematode infection after 3, 6, and 9 days. Among them, Aradu.T5WNW, Aradu.JM17V, and Aradu.MKP1A were up-regulate at these three time points, while Aradu.QD5DS and Aradu.M0ENQ were up-regulated 6 and 9 days after nematode infection. The expression of the above mentioned five genes was significantly and negatively correlated with the number of LRR8 domain, indicating that fewer LRR8 domains are associated with the promotion of LRR-containing genes that resist nematode infection. Patterns of co-expression and cis-acting elements indicated that WRKY possibly regulate the responses of LRR-containing genes to nematode infection and that expansin genes may work together with LRR-containing genes in response to nematode infection.

CONCLUSIONS

We identified the number and type of LRR-containing genes in A. duranensis. The LRR-containing genes that were found appear to be involved in responses to nematode infection. The number of LRR8 domains was negatively correlated with expression after nematode infection. The WRKY transcription factor may regulate resistance to nematode infection based on LRR-containing genes. Our results could improve the understanding of resistance to nematodes and molecular breeding in peanuts.

Comparative analysis of selection mode reveals different evolutionary rate and expression pattern in Arachis duranensis and Arachis ipaënsis duplicated genes

Our results reveal that K_s is a determining factor affecting selective pressure and different evolution and expression patterns are detected between PSGs and NSGs in wild Arachis duplicates. Selective pressure, including purifying (negative) and positive selection, can be detected in organisms. However, studies on comparative evolutionary rates, gene expression patterns and gene features between negatively selected genes (NSGs) and positively selected genes (PSGs) are lagging in paralogs of plants. Arachis duranensis and Arachis ipaënsis are ancestors of the cultivated peanut, an important oil and protein crop. Here, we carried out a series of systematic analyses, comparing NSG and PSG in paralogs, using genome sequences and transcriptome datasets in A. duranensis and A. ipaënsis. We found that synonymous substitution rate (K_s) is a determining factor affecting selective pressure in A. duranensis and A. ipaënsis duplicated genes. Lower expression level, lower gene expression breadth, higher codon bias and shorter polypeptide length were found in PSGs and not in NSGs. The correlation analyses showed that gene expression breadth was positively correlated with polypeptide length and GC content at the first codon site (GC1) in PSGs and NSGs, respectively. There was a negative correlation between expression level and polypeptide length in PSGs. In NSGs, the K_s was positively correlated with expression level, gene expression breadth, GC1, and GC content at the third codon site (GC3), but selective pressure was negatively correlated with expression level, gene expression breadth, polypeptide length, GC1, and GC3 content. The function of most duplicated gene pairs was divergent under drought and nematode stress. Taken together, our results show that different evolution and expression patterns occur between PSGs and NSGs in paralogs of two wild Arachis species.

WRKY transcription factors in legumes

BACKGROUND

WRKY transcription factors, so named because of the WRKYGQK heptapeptide at the N-terminal end, are widely distributed in plants and play an important role in physiological changes and response to biotic and abiotic stressors. Many previous studies have focused on the evolution of WRKY transcription factors in a given plant; however, little is known about WRKY evolution in legumes. The gene expression pattern of duplicated WRKY transcription factors remains unclear.

RESULTS

We first identified the WRKY proteins in 12 legumes. We found that the WRKYGQK heptapeptide tended to mutate into WRKYGKK. The Q site in WRKYGQK preferentially mutated, while W, K, and Y were conserved. The phylogenetic tree shows that the WRKY proteins in legumes have multiple origins, especially group IIc. For example, WRKY64 from Lupinus angustifolius (LaWRKY64) contains three WRKY domains, of which the first two clustered together in the N-terminal WRKY domain of the group I WRKY protein, and the third WRKY domain grouped in the C-terminal WRKY domain of the group I WRKY protein. Orthologous WRKY genes have a faster evolutionary rate and are subject to constrained selective pressure, unlike paralogous WRKY genes. Different gene features were observed between duplicated WRKY genes and singleton WRKY genes. Duplicated Glycine max WRKY genes with similar gene features have gene expression divergence.

CONCLUSIONS

We analyzed the WRKY number and type in 12 legumes, concluding that the WRKY proteins have multiple origins. A novel WRKY protein, LaWRKY64, was found in L. angustifolius. The first two WRKY domains of LaWRKY64 have the same origin. The orthologous and paralogous WRKY proteins have different evolutionary rates. Duplicated WRKY genes have gene expression divergence under normal growth conditions in G. max. These results provide insight into understanding WRKY evolution and expression.

Transcriptome analysis of the tea oil camellia (Camellia oleifera) reveals candidate drought stress genes

BACKGROUND

The tea-oil camellia (Camellia oleifera) is the most important oil plant in southern China, and has a strong resistance to drought and barren soil. Understanding the molecular mechanisms of drought tolerance would greatly promote its cultivation and molecular breeding.

RESULTS

In total, we obtained 76,585 unigenes with an average length of 810 bp and an N50 of 1,092 bp. We mapped all the unigenes to the NCBI 'nr' (non-redundant), SwissProt, KEGG, and clusters of orthologous groups (COG) databases, where 52,531 (68.6%) unigenes were functionally annotated. According to the annotation, 46,171 (60.8%) unigenes belong to 338 KEGG pathways. We identified a series of unigenes that are related to the synthesis and regulation of abscisic acid (ABA), the activity of protective enzymes, vitamin B6 metabolism, the metabolism of osmolytes, and pathways related to the biosynthesis of secondary metabolites. After exposed to drought for 12 hours, the number of differentially-expressed genes (DEGs) between treated plants and control plants increased in the G4 cultivar, while there was no significant increase in the drought-tolerant C3 cultivar. DEGs associated with drought stress responsive pathways were identified by KEGG pathway enrichment analysis. Moreover, we found 789 DEGs related to transcription factors. Finally, according to the results of qRT-PCR, the expression levels of the 20 unigenes tested were consistent with the results of next-generation sequencing.

CONCLUSIONS

In the present study, we identified a large set of cDNA unigenes from C. oleifera annotated using public databases. Further studies of DEGs involved in metabolic pathways related to drought stress and transcription will facilitate the discovery of novel genes involved in resistance to drought stress in this commercially important plant.

Comparative analysis of NBS-LRR genes and their response to Aspergillus flavus in Arachis

Studies have demonstrated that nucleotide-binding site-leucine-rich repeat (NBS-LRR) genes respond to pathogen attack in plants. Characterization of NBS-LRR genes in peanut is not well documented. The newly released whole genome sequences of Arachis duranensis and Arachis ipaënsis have allowed a global analysis of this important gene family in peanut to be conducted. In this study, we identified 393 (AdNBS) and 437 (AiNBS) NBS-LRR genes from A. duranensis and A. ipaënsis, respectively, using bioinformatics approaches. Full-length sequences of 278 AdNBS and 303 AiNBS were identified. Fifty-one orthologous, four AdNBS paralogous, and six AiNBS paralogous gene pairs were predicted. All paralogous gene pairs were located in the same chromosomes, indicating that tandem duplication was the most likely mechanism forming these paralogs. The paralogs mainly underwent purifying selection, but most LRR 8 domains underwent positive selection. More gene clusters were found in A. ipaënsis than in A. duranensis, possibly owing to tandem duplication events occurring more frequently in A. ipaënsis. The expression profile of NBS-LRR genes was different between A. duranensis and A. hypogaea after Aspergillus flavus infection. The up-regulated expression of NBS-LRR in A. duranensis was continuous, while these genes responded to the pathogen temporally in A. hypogaea.