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

Molecular insights into Solanum sisymbriifolium's resistance against Globodera pallida via RNA-seq

BACKGROUND

The presence of potato cyst nematodes (PCN) causes a significant risk to potato crops globally, leading to reduced yields and economic losses. While the plant Solanum sisymbriifolium is known for its resistance to PCN and can be used as a trap crop, the molecular mechanisms behind this resistance remain poorly understood. In this study, genes differentially expressed were identified in control and infected plants during the early stages of the S. sisymbriifolium - G. pallida interaction.

RESULTS

Gene expression profiles were characterized for two S. sisymbriifolium cultivars, Melody and Sis6001, uninfected and infected by G. pallida. The comparative transcriptome analysis revealed a total of 4,087 and 2,043 differentially expressed genes (DEGs) in response to nematode infection in the cultivars Melody and Sis6001, respectively. Gene ontology (GO) enrichment analysis provided insights into the response of the plant to nematode infection, indicating an activation of the plant metabolism, oxidative stress leading to defence mechanism activation, and modification of the plant cell wall. Genes associated with the jasmonic and salicylic acid pathways were also found to be differentially expressed, suggesting their involvement in the plant's defence response. In addition, the analysis of NBS-LRR domain-containing transcripts that play an important role in hypersensitive response and programmed cell death led to the identification of ten transcripts that had no annotations from the databases, with emphasis on TRINITY_DN52667_C1_G1, found to be upregulated in both cultivars.

CONCLUSIONS

These findings represent an important step towards understanding the molecular basis underlying plant resistance to nematodes and facilitating the development of more effective control strategies against PCN.

Whole-genome re-sequencing, diversity analysis, and stress-resistance analysis of 77 grape rootstock genotypes

INTRODUCTION

Grape rootstocks play critical role in the development of the grape industry over the globe for their higher adaptability to various environments, and the evaluation of their genetic diversity among grape genotypes is necessary to the conservation and utility of genotypes.

METHODS

To analyze the genetic diversity of grape rootstocks for a better understanding multiple resistance traits, whole-genome re-sequencing of 77 common grape rootstock germplasms was conducted in the present study.

RESULTS

About 645 billion genome sequencing data were generated from the 77 grape rootstocks at an average depth of ~15.5×, based on which the phylogenic clusters were generated and the domestication of grapevine rootstocks was explored. The results indicated that the 77 rootstocks originated from five ancestral components. Through phylogenetic, principal components, and identity-by-descent (IBD) analyses, these 77 grape rootstocks were assembled into ten groups. It is noticed that the wild resources of V. amurensis and V. davidii, originating from China and being generally considered to have stronger resistance against biotic and abiotic stresses, were sub-divided from the other populations. Further analysis indicated that a high level of linkage disequilibrium was found among the 77 rootstock genotypes, and a total of 2,805,889 single nucleotide polymorphisms (SNPs) were excavated, GWAS analysis among the grape rootstocks located 631, 13, 9, 2, 810, and 44 SNP loci that were responsible to resistances to phylloxera, root-knot nematodes, salt, drought, cold and waterlogging traits.

DISCUSSION

This study generated a significant amount of genomic data from grape rootstocks, thus providing a theoretical basis for further research on the resistance mechanism of grape rootstocks and the breeding of resistant varieties. These findings also reveal that China originated V. amurensis and V. davidii could broaden the genetic background of grapevine rootstocks and be important germplasm used in breeding high stress-resistant grapevine rootstocks.

Deciphering evolutionary dynamics of WRKY genes in Arachis species

BACKGROUND

Cultivated peanut (Arachis hypogaea), a progeny of the cross between A. duranensis and A. ipaensis, is an important oil and protein crop from South America. To date, at least six Arachis genomes have been sequenced. WRKY transcription factors (TFs) play crucial roles in plant growth, development, and response to abiotic and biotic stresses. WRKY TFs have been identified in A. duranensis, A. ipaensis, and A. hypogaea cv. Tifrunner; however, variations in their number and evolutionary patterns across various Arachis spp. remain unclear.

RESULTS

WRKY TFs were identified and compared across different Arachis species, including A. duranensis, A. ipaensis, A. monticola, A. hypogaea cultivars (cv.) Fuhuasheng, A. hypogaea cv. Shitouqi, and A. hypogaea cv. Tifrunner. The results showed that the WRKY TFs underwent dynamic equilibrium between diploid and tetraploid peanut species, characterized by the loss of old WRKY TFs and retention of the new ones. Notably, cultivated peanuts inherited more conserved WRKY orthologs from wild tetraploid peanuts than their wild diploid donors. Analysis of the W-box elements and protein-protein interactions revealed that different domestication processes affected WRKY evolution across cultivated peanut varieties. WRKY TFs of A. hypogaea cv. Fuhuasheng and Shitouqi exhibited a similar domestication process, while those of cv. Tifrunner of the same species underwent a different domestication process based on protein-protein interaction analysis.

CONCLUSIONS

This study provides new insights into the evolution of WRKY TFs in Arachis spp.

Vibrio splendidus flagellin C binds tropomodulin to induce p38 MAPK-mediated p53-dependent coelomocyte apoptosis in Echinodermata

As a typical pathogen-associated molecular pattern, bacterial flagellin can bind Toll-like receptor 5 and the intracellular NAIP5 receptor component of the NLRC4 inflammasome to induce immune responses in mammals. However, these flagellin receptors are generally poorly understood in lower animal species. In this study, we found that the isolated flagellum of Vibrio splendidus AJ01 destroyed the integrity of the tissue structure of coelomocytes and promoted apoptosis in the sea cucumber Apostichopus japonicus. To further investigate the molecular mechanism, the novel intracellular LRR domain-containing protein tropomodulin (AjTmod) was identified as a protein that interacts with flagellin C (FliC) with a dissociation constant (K_d) of 0.0086 ± 0.33 μM by microscale thermophoresis assay. We show that knockdown of AjTmod also depressed FliC-induced apoptosis of coelomocytes. Further functional analysis with different inhibitor treatments revealed that the interaction between AjTmod and FliC could specifically activate p38 MAPK, but not JNK or ERK MAP kinases. We demonstrate that the transcription factor p38 is then translocated into the nucleus, where it mediates the expression of p53 to induce coelomocyte apoptosis. Our findings provide the first evidence that intracellular AjTmod serves as a novel receptor of FliC and mediates p53-dependent coelomocyte apoptosis by activating the p38 MAPK signaling pathway in Echinodermata.

Identification of miRNA and their target genes in Cestrum nocturnum L. and Cestrum diurnum L. in stress responses

MicroRNAs (miRNAs) are small, highly conserved non-coding RNA molecules and products of primary miRNAs that regulate the target gene expression. Homology-based approaches were employed to identify miRNAs and their targets in Cestrum nocturnum L. and Cestrum diurnum L. A total of 32 and 12 miRNA candidates were identified in C. nocturnum and C. diurnum. These miRNAs belong to 26 and 10 miRNA families and regulate 1024 and 1007 target genes in C. nocturnum, and C. diurnum, respectively. The functional roles of these miRNAs have not been earlier elucidated in Cestrum. MiR815a, miR849, miR1089 and miR172 have a strong propensity to target genes controlling phytochrome-interacting factor 1 (PIF1), ubiquitin-specific protease 12 (UBP12), leucine-rich repeat (LRR) protein kinase and GAI, RGA, SCR (GRAS) family transcription factor in C. nocturnum. While miR5205a, miR1436 and miR530 regulate PATATIN-like protein 6 (PLP6), PHD finger transcription factor and myb domain protein 48 (MYB48) in C. diurnum. Overall, these miRNAs have regulatory responses in biotic and abiotic stresses in both plant species. Eight putative miRNAs and their target genes were selected for qRT-PCR validation. The validated results suggested the importance of miR815a, miR849, miR5205a, miR1089, miR172, miR1436, and miR530 in exerting control over stress responses in C. nocturnum and C. diurnum.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-022-01127-1.

Transcriptome analysis reveals the mechanism of improving erect-plant-type peanut yield by single-seeding precision sowing

Background

In China, double-seed (DS) sowing (i.e., sowing two seeds per hole) has been conventionally performed towards the erect-plant-type peanuts to increase the low germination rate due to poor seed preservation conditions. However, the corresponding within-hole plant competition usually limits the subsequent plant growth and the final yield. We developed a high-yield cultivation system of single-seed (SS) precision sowing to solve this paradox, saving 20% of seeds and increasing yields by more than 10% relative to the conventional DS sowing.

Methods

To explore the mechanisms of these two different cropping patterns in peanut yields, we conducted transcriptomic and physiological comparisons in the seeding plant leaf and root tissues between SS precision sowing and standard DS sowing treatments.

Results

After assembly, each library contained an average of 43 million reads and generated a total of 523,800, 338 clean reads. After GO and Kyoto Encyclopedia of Genes and Genomes pathway analysis, we found the key genes for biotic and abiotic stress showed higher expression in roots of plants grown under the SS precision sowing treatment, including genes encoding disease resistance, oxidation-reduction, hormone related, and stress response transcription factors and signaling regulation proteins. In particular, the resveratrol synthesis genes related to stress and disease resistance appeared induced in roots under the SS sowing treatment.

Conclusion

These data indicated that Aspergillus flavus resistance and stress tolerance in roots under SS precision sowing were enhanced compared with roots under the DS sowing treatment. This work benefits the development of underground pods and thus increasing peanut yields.

Comprehensive analysis of coding sequence architecture features and gene expression in Arachis duranensis

Coding sequence (CDS) architecture affects gene expression levels in organisms. Codon optimization can increase the gene expression level. Therefore, understanding codon usage patterns has important implications for research on genetic engineering and exogenous gene expression. To date, the codon usage patterns of many model plants have been analyzed. However, the relationship between CDS architecture and gene expression in Arachis duranensis remains poorly understood. According to the results of genome sequencing, A. duranensis has many resistant genes that can be used to improve the cultivated peanut. In this study, bioinformatic approaches were used to estimate A. duranensis CDS architectures, including frequency of the optimal codon (Fop), polypeptide length and GC contents at the first (GC1), second (GC2) and third (GC3) codon positions. In addition, Arachis RNA-seq datasets were downloaded from PeanutBase. The relationships between gene expression and CDS architecture were assessed both under normal growth as well as nematode and drought stress conditions. A total of 26 codons with high frequency were identified, which preferentially ended with A or T in A. duranensis CDSs under the above-mentioned three conditions. A similar CDS architecture was found in differentially expressed genes (DEGs) under nematode and drought stresses. The GC1 content differed between DEGs and non-differentially expressed genes (NDEGs) under both drought and nematode stresses. The expression levels of DEGs were affected by different CDS architectures compared with NDEGs under drought stress. In addition, no correlation was found between differential gene expression and CDS architecture neither under nematode nor under drought stress. These results aid the understanding of gene expression in A. duranensis.

Genome-Wide Identification of Populus Malectin/Malectin-Like Domain-Containing Proteins and Expression Analyses Reveal Novel Candidates for Signaling and Regulation of Wood Development

Malectin domain (MD) is a ligand-binding protein motif of pro- and eukaryotes. It is particularly abundant in Viridiplantae, where it occurs as either a single (MD, PF11721) or tandemly duplicated domain (PF12819) called malectin-like domain (MLD). In herbaceous plants, MD- or MLD-containing proteins (MD proteins) are known to regulate development, reproduction, and resistance to various stresses. However, their functions in woody plants have not yet been studied. To unravel their potential role in wood development, we carried out genome-wide identification of MD proteins in the model tree species black cottonwood (Populus trichocarpa), and analyzed their expression and co-expression networks. P. trichocarpa had 146 MD genes assigned to 14 different clades, two of which were specific to the genus Populus. 87% of these genes were located on chromosomes, the rest being associated with scaffolds. Based on their protein domain organization, and in agreement with the exon-intron structures, the MD genes identified here could be classified into five superclades having the following domains: leucine-rich repeat (LRR)-MD-protein kinase (PK), MLD-LRR-PK, MLD-PK (CrRLK1L), MLD-LRR, and MD-Kinesin. Whereas the majority of MD genes were highly expressed in leaves, particularly under stress conditions, eighteen showed a peak of expression during secondary wall formation in the xylem and their co-expression networks suggested signaling functions in cell wall integrity, pathogen-associated molecular patterns, calcium, ROS, and hormone pathways. Thus, P. trichocarpa MD genes having different domain organizations comprise many genes with putative foliar defense functions, some of which could be specific to Populus and related species, as well as genes with potential involvement in signaling pathways in other tissues including developing wood.

Genome-wide identification and characterization of nonspecific lipid transfer protein (nsLTP) genes in Arachis duranensis

Nonspecific lipid transfer proteins (nsLTPs) play vital roles in lipid metabolism, cell apoptosis and biotic and abiotic stresses in plants. However, the distribution of nsLTPs in Arachis duranensis has not been fully characterized. In this study, we identified 64 nsLTP genes in A. duranensis (designated AdLTPs), which were classified into six subfamilies and randomly distributed along nine chromosomes. Tandem and segmental duplication events were detected in the evolution of AdLTPs. The K_s and ω values differed significantly between Types 1 and D subfamilies, and eight AdLTPs were under positive selection. The expression levels of AdLTPs were changed after salinity, PEG, low-temperature and ABA treatments. Three AdLTPs were associated with resistance to nematode infection, and DOF and WRI1 transcription factors may regulate the AdLTP response to nematode infection. Our results may provide valuable genomic information for the breeding of peanut cultivars that are resistant to biotic and abiotic stresses.

Identification and expression analysis of WRKY gene family under drought stress in peanut (Arachis hypogaea L.)

WRKY transcription factors play crucial roles in regulation mechanism leading to the adaption of plants to the complex environment. In this study, AhWRKY family was comprehensively analyzed using bioinformatic approaches in combination with transcriptome sequencing data of the drought-tolerant peanut variety ‘L422’. A total of 158 AhWRKY genes were identified and named according to their distribution on the chromosomes. Based on the structural features and phylogenetic analysis of AhWRKY proteins, the AhWRKY family members were classified into three (3) groups, of which group II included five (5) subgroups. Results of structure and conserved motifs analysis for the AhWRKY genes confirmed the accuracy of the clustering analysis. In addition, 12 tandem and 136 segmental duplication genes were identified. The results indicated that segmental duplication events were the main driving force in the evolution of AhWRKY family. Collinearity analysis found that 32 gene pairs existed between Arachis hypogaea and two diploid wild ancestors (Arachis duranensis and Arachis ipaensis), which provided valuable clues for phylogenetic characteristics of AhWRKY family. Furthermore, 19 stress-related cis-acting elements were found in the promoter regions. During the study of gene expression level of AhWRKY family members in response to drought stress, 73 differentially expressed AhWRKY genes were obtained to have been influenced by drought stress. These results provide fundamental insights for further study of WRKY genes in peanut drought resistance.

Evolutionary balance between LRR domain loss and young NBS-LRR genes production governs disease resistance in Arachis hypogaea cv. Tifrunner

BACKGROUND

Cultivated peanut (Arachis hypogaea L.) is an important oil and protein crop, but it has low disease resistance; therefore, it is important to reveal the number, sequence features, function, and evolution of genes that confer resistance. Nucleotide-binding site-leucine-rich repeats (NBS-LRRs) are resistance genes that are involved in response to various pathogens.

RESULTS

We identified 713 full-length NBS-LRRs in A. hypogaea cv. Tifrunner. Genetic exchange events occurred on NBS-LRRs in A. hypogaea cv. Tifrunner, which were detected in the same subgenomes and also found in different subgenomes. Relaxed selection acted on NBS-LRR proteins and LRR domains in A. hypogaea cv. Tifrunner. Using quantitative trait loci (QTL), we found that NBS-LRRs were involved in response to late leaf spot, tomato spotted wilt virus, and bacterial wilt in A. duranensis (2 NBS-LRRs), A. ipaensis (39 NBS-LRRs), and A. hypogaea cv. Tifrunner (113 NBS-LRRs). In A. hypogaea cv. Tifrunner, 113 NBS-LRRs were classified as 75 young and 38 old NBS-LRRs, indicating that young NBS-LRRs were involved in response to disease after tetraploidization. However, compared to A. duranensis and A. ipaensis, fewer LRR domains were found in A. hypogaea cv. Tifrunner NBS-LRR proteins, partly explaining the lower disease resistance of the cultivated peanut.

CONCLUSIONS

Although relaxed selection acted on NBS-LRR proteins and LRR domains, LRR domains were preferentially lost in A. hypogaea cv. Tifrunner compared to A. duranensis and A. ipaensis. The QTL results suggested that young NBS-LRRs were important for resistance against diseases in A. hypogaea cv. Tifrunner. Our results provid insight into the greater susceptibility of A. hypogaea cv. Tifrunner to disease compared to A. duranensis and A. ipaensis.

Molecular and transcriptional characterization of phosphatidyl ethanolamine-binding proteins in wild peanuts Arachis duranensis and Arachis ipaensis

BACKGROUND

Phosphatidyl ethanolamine-binding proteins (PEBPs) are involved in the regulation of plant architecture and flowering time. The functions of PEBP genes have been studied in many plant species. However, little is known about the characteristics and expression profiles of PEBP genes in wild peanut species, Arachis duranensis and Arachis ipaensis, the diploid ancestors of cultivated peanuts.

RESULTS

In this study, genome-wide identification methods were used to identify and characterize a total of 32 peanut PEBP genes, 16 from each of the two wild peanut species, A. duranensis and A. ipaensis. These PEBP genes were classified into 3 groups (TERMINAL FLOWER1-like, FLOWERING LOCUS T-like, and MOTHER OF FT AND TFL1-like) based on their phylogenetic relationships. The gene structures, motifs, and chromosomal locations for each of these PEBPs were analyzed. In addition, 4 interchromosomal duplications and 1 tandem duplication were identified in A. duranensis, and 2 interchromosomal paralogs and 1 tandem paralog were identified in A. ipaensis. Ninety-five different cis-acting elements were identified in the PEBP gene promoter regions and most genes had different numbers and types of cis-elements. As a result, the transcription patterns of these PEBP genes varied in different tissues and under long day and short day conditions during different growth phases, indicating the functional diversities of PEBPs in different tissues and their potential functions in plant photoperiod dependent developmental pathways. Moreover, our analysis revealed that AraduF950M/AraduWY2NX in A. duranensis, and Araip344D4/Araip4V81G in A. ipaensis are good candidates for regulating plant architecture, and that Aradu80YRY, AraduYY72S, and AraduEHZ9Y in A. duranensis and AraipVEP8T in A. ipaensis may be key factors regulating flowering time.

CONCLUSION

Sixteen PEBP genes were identified and characterized from each of the two diploid wild peanut genomes, A. duranensis and A. ipaensis. Genetic characterization and spatio-temporal expression analysis support their importance in plant growth and development. These findings further our understanding of PEBP gene functions in plant species.

Systematic Identification and Evolution Analysis of Sox Genes in Coturnix japonica Based on Comparative Genomics

Coturnix japonica (Japanese quail) has been extensively used as a model animal for biological studies. The Sox gene family, which was systematically characterized by a high-mobility group (HMG-box) in many animal species, encodes transcription factors that play central roles during multiple developmental processes. However, genome-wide investigations on the Sox gene family in birds are scarce. In the current study, we first performed a genome-wide study to explore the Sox gene family in galliform birds. Based on available genomic sequences retrieved from the NCBI database, we focused on the global identification of the Sox gene family in C. japonica and other species in Galliformes, and the evolutionary relationships of Sox genes. In our result, a total of 35 Sox genes in seven groups were identified in the C. japonica genome. Our results also revealed that dispersed gene duplications contributed the most to the expansion of the Sox gene family in Galliform birds. Evolutionary analyses indicated that Sox genes are an ancient gene family, and strong purifying selections played key roles in the evolution of CjSox genes of C. japonica. More interestingly, we observed that most Sox genes exhibited highly embryo-specific expression in both gonads. Our findings provided new insights into the molecular function and phylogeny of Sox gene family in birds.